原创 :最近刚完成自己的论文模型部分,变截面梁模拟实际桥梁.
以下为建模的一段命令流:
/prep7
/title,the analyse of the 3rd Bridge
et,1,plane82
*do,i,1,54,1
*if,i,le,25,then
xp=i*3/2
*elseif,i,le,53,then
xp=25*3/2+(i-25)*4/2
*else
xp=94.5
*endif
height=(8/94.5**1.5)*(94.5-xp)**1.5+3.57
tb=1.2-0.0093122*xp
*if,i,le,24,then
tw=0.7
*elseif,i,le,42,then
tw=0.6
*else
tw=0.5
*endif
x1=7.0/2-tw-1.5
x2=7.0/2-tw
x3=7.0/2-tw-1.0
k,1,-x1,x1*0.02,0
k,2,-x2,x2*0.02,0
k,3,-3.5,0.07,0
k,4,-5,0.1,0
k,5,-6.7,0.134,0
k,6,-6.7,-0.016,0
k,7,-5,-0.2,0
k,8,-3.5,-0.73,0
k,9,-x2,x1*0.02-0.3-0.47,0
k,10,-x1,x1*0.02-0.3,0
k,11,-3.5,-(height-tb-0.5),0
k,12,-x2,-(height-tb-0.5),0
k,13,-x3,-(height-tb),0
k,14,-x3,-height,0
k,15,-x2,-height,0
k,16,-3.5,-height,0
k,21,x1,-x1*0.02,0
k,22,x2,-x2*0.02,0
k,23,3.5,-0.07,0
k,24,5,-0.1,0
k,25,6.7,-0.134,0
k,26,6.7,-0.284,0
k,27,5,-0.4,0
k,28,3.5,-0.87,0
k,29,x2,-0.30-x1*0.02-0.53,0
k,30,x1,-0.30-x1*0.02,0
k,31,3.5,-(height-tb-0.5),0
k,32,x2,-(height-tb-0.5),0
k,33,x3,-(height-tb),0
k,34,x3,-height,0
k,35,x2,-height,0
k,36,3.5,-height,0
!
*do,ip,1,4,1
a,ip,11-ip,10-ip,ip+1
a,20+ip,20+11-ip,20+10-ip,20+ip+1
*enddo
a,1,21,30,10
a,8,9,12,11
a,11,12,15,16
a,12,13,14,15
a,28,31,32,29
a,31,36,35,32
a,32,35,34,33
a,13,33,34,14
esize,,1
amesh,all
secwrite,girder%i%,sect,,1
sectype,1000+i,beam,mesh
secoffset,cent !坐标放在节点上
secread,girder%i%,'sect',,mesh
aclear,all
adele,all,,,1
*enddo
!define the section of 0# segment
k,1,-0.8003,0.016,0
k,2,-2.3,0.046,0
k,3,-3.5,0.07,0
k,4,-5,0.1,0
k,5,-6.7,0.134,0
k,6,-6.7,-0.016,0
k,7,-5,-0.2,0
k,8,-3.5,-0.73,0
k,9,-2.3,-0.7542,0
k,10,-0.8003,-0.4841,0
k,11,-3.5,-9.57,0
k,12,-2.3,-9.57,0
k,13,-1.3,-10.07,0
k,14,-1.3,-11.57,0
k,15,-2.3,-11.57,0
k,16,-3.5,-11.57,0
k,21,0.8003,-0.016,0
k,22,2.3,-0.046,0
k,23,3.5,-0.07,0
k,24,5,-0.1,0
k,25,6.7,-0.134,0
k,26,6.7,-0.284,0
k,27,5,-0.4,0
k,28,3.5,-0.87,0
k,29,2.3,-0.8462,0
k,30,0.8003,-0.5161,0
k,31,3.5,-9.57,0
k,32,2.3,-9.57,0
k,33,1.3,-10.07,0
k,34,1.3,-11.57,0
k,35,2.3,-11.57,0
k,36,3.5,-11.57,0
*do,ip,1,4,1
a,ip,11-ip,10-ip,ip+1
a,20+ip,20+11-ip,20+10-ip,20+ip+1
*enddo
a,1,21,30,10
a,8,9,12,11
a,11,12,15,16
a,12,13,14,15
a,28,31,32,29
a,31,36,35,32
a,32,35,34,33
a,13,33,34,14
esize,,1
amesh,all
secwrite,girder0,sect,,1
sectype,1000,beam,mesh
secoffset,origin
secread,girder0,'sect',,mesh
aclear,all
adele,all,,,1
!define the section of the main pier
rectng,-3.5,3.5,-1.25,1.25
cyl4,2.9234375,0,1.3765625,65.23848614,,-65.23848614
cyl4,-2.9234375,0,1.3765625,114.76151386,,-114.76151386+360
aadd,all
esize,1
smartsize,5
amesh,all
secwrite,pier0,sect,,1
sectype,500,beam,mesh
secoffset,origin
secread,pier0,'sect',,mesh
aclear,all
adele,all,,,1
!
etdele,1
et,1,beam188
et,2,beam44
keyopt,1,1,1
keyopt,1,2,0
keyopt,1,3,1
keyopt,1,4,2
keyopt,1,6,2
keyopt,1,7,2
keyopt,1,8,2
keyopt,1,9,2
!50# concrete
mp,ex,1,3.5e10
mp,dens,1,2600
mp,prxy,1,0.1667
MP,alpx,1,1.0e-5
!40# concrete
mp,ex,2,3.3e10
mp,dens,2,2600
mp,prxy,2,0.1667
MP,alpx,2,1.0e-5
*dim,coorx,array,245
*dim,coory,array,245
*vread,coorx,coor1,dat
(f8.3)
*vread,coory,coor2,dat
(f8.3)
*do,i,1,245,1
n,i,coorx(i),0.0,0.0
*enddo
n,10000,0,50
!0# segment
r,1000
real,1000
secnum,1000
*do,i,61,68,1
e,i,i+1,10000
e,i+116,i+117,10000
*enddo
!
r,1
real,1
secnum,1054
*do,i,1,6,1
e,i,i+1,10000
*enddo
!
*do,i,7,60,1
r,i
real,i
sectype,i,taper
secdata,1001+(60-i),coorx(i),0,0
secdata,1000+(60-i),coorx(i+1),0,0
secnum,i
e,i,i+1,10000
*enddo
!
*do,i,69,122,1
r,i
real,i
sectype,i,taper
secdata,1000+(i-69),coorx(i),0,0
secdata,1001+(i-69),coorx(i+1),0,0
secnum,i
e,i,i+1,10000
*enddo
!
*do,i,123,176,1
r,i
real,i
sectype,i,taper
secdata,1001+(176-i),coorx(i),0,0
secdata,1000+(176-i),coorx(i+1),0,0
secnum,i
e,i,i+1,10000
*enddo
!
*do,i,185,238,1
r,i
real,i
sectype,i,taper
secdata,1000+(i-185),coorx(i),0,0
secdata,1001+(i-185),coorx(i+1),0,0
secnum,i
e,i,i+1,10000
*enddo
r,301
real,301
secnum,1054
*do,i,239,244,1
e,i,i+1,10000
*enddo
!2#pier & 3# pier
!type,2
mat,2
n,500,105.75,-11.5+105.75*0.005,0
n,511,105.75,-49.56+105.75*0.005,0
n,550,114.25,-11.5+114.25*0.005,0
n,561,114.25,-11.5-38.10+114.25*0.005,0
n,600,305.75,-11.5+305.75*0.005,0
n,611,305.75,-11.5-45.61+305.75*0.005,0
n,650,314.25,-11.5+314.25*0.005,0
n,661,314.25,-11.5-45.65+314.25*0.005,0
fill,500,511
fill,550,561
fill,600,611
fill,650,661
r,500
real,500
secnum,500
*do,i,500,510,1
e,i,i+1,1
e,i+50,i+51,1
e,i+100,i+101,1
e,i+150,i+151,1
*enddo
csys,0
/eshape,1
eplot
d,511,all
d,561,all
d,611,all
d,661,all
d,1,uy,,,,,uz,rotx
d,245,uy,,,,,uz,rotx
cerig,62,500,all
cerig,68,550,all
cerig,178,600,all
cerig,184,650,all
allsel,all
nummrg,all
numcmp,all
etcontrol,set,
esel,r,mat,,2
cm,piere,elem
nsle,s
cm,piern,node
allsel,all |