fini
/clear
/title,three span continous bridge
/prep7
rad=200
h=1.6
a=18.5
b=12.8
b0=0.50!中隔板厚度
a0=3.60!箱室宽度
csys,1
k,,rad-b/2
k,,rad-a0/2-b0/2
k,,rad+a0/2+b0/2
k,,rad+b/2
k,,rad-b/2-0.2,,h
k,,rad-a0/2-b0/2,,h
k,,rad+a0/2+b0/2,,h
k,,rad+b/2+0.2,,h
k,,rad-a/2,,h
k,,rad+a/2,,h
kgen,2,1,10,1,,11.55/rad*180/3.14,,10
kgen,9,11,20,1,,1.25/rad*180/3.14,,10
kgen,3,91,100,1,,11.55/rad*180/3.14,,10
kgen,5,111,120,1,,1.25/rad*180/3.14,,10
et,1,shell63
mp,dens,1,2500
mp,ex,1,3.5e10
mp,prxy,1,0.1667
r,1,0.97
a,1,2,3,4,8,7,6,5
aatt,1,1,1
esize,1
amesh,all
asel,none
r,2,2
a,51,52,53,54,58,57,56,55
aatt,1,2,1
esize,1
amesh,all
asel,none
r,3,0.3
a,101,102,103,104,108,107,106,105
aatt,1,3,1
esize,1
amesh,all
asel,none
r,4,0.97
a,151,152,153,154,158,157,156,155
aatt,1,4,1
esize,1
amesh,all
asel,none
r,5,0.7
*do,i,1,15
a,5+(i-1)*10,5+10*i,1+10*i,1+(i-1)*10
a,8+(i-1)*10,8+10*i,4+10*i,4+(i-1)*10
aatt,1,5,1
esize,1
amesh,all
asel,none
*enddo
r,6,0.22
*do,k,1,3,1
a,k,k+1,k+1+10,k+10
aatt,1,6,1
esize,1.0
amesh,all
asel,none
*enddo
!底板厚度渐变
*dim,hd,array,5
hd(1)=0.42-(0.42-0.22)/4*4
hd(2)=0.42-(0.42-0.22)/4*3
hd(3)=0.42-(0.42-0.22)/4*2
hd(4)=0.42-(0.42-0.22)/4*1
hd(5)=0.42
*do,i,1,4,1
r,10+i,hd(i),hd(i),hd(i+1),hd(i+1) !定义渐变的底板厚度
*enddo
*do,i,1,4,1
*do,k,1,3,1
a,k+i*10,k+1+i*10,k+1+(i+1)*10,k+(i+1)*10
aatt,1,10+i,1
esize,1.0
amesh,all
asel,none
*enddo
*enddo
*do,i,5,8,1
*do,k,1,3,1
a,k+1+(i+1)*10,k+(i+1)*10,k+i*10,k+1+i*10
aatt,1,19-i,1
esize,1.0
amesh,all
asel,none
*enddo
*enddo
*do,i,10,11,1
*do,k,1,3,1
a,k+(i-1)*10,k+1+(i-1)*10,k+1+i*10,k+i*10
aatt,1,6,1
esize,1.0
amesh,all
asel,none
*enddo
*enddo
*do,i,11,14,1
*do,k,1,3,1
a,k+i*10,k+1+i*10,k+1+(i+1)*10,k+(i+1)*10
aatt,1,i,1
esize,1.0
amesh,all
asel,none
*enddo
*enddo
r,7,0.5
*do,i,1,15
a,6+(i-1)*10,2+(i-1)*10,2+10*i,6+10*i
a,7+(i-1)*10,3+(i-1)*10,3+10*i,7+10*i
aatt,1,7,1
esize,1.0
amesh,all
asel,none
*enddo
r,8,0.23
*do,i,1,15
a,9+(i-1)*10,5+(i-1)*10,5+10*i,9+10*i
a,5+(i-1)*10,6+(i-1)*10,6+10*i,5+10*i
a,6+(i-1)*10,7+(i-1)*10,7+10*i,6+10*i
a,7+(i-1)*10,8+(i-1)*10,8+10*i,7+10*i
a,8+(i-1)*10,10+(i-1)*10,10+10*i,8+10*i
aatt,1,8,1
esize,1.0
amesh,all
asel,none
*enddo
allsel,all
nsym,y,100000,all
esym,y,100000,all
!支座
nsel,s,,,54
nsel,a,,,61
nsel,a,,,150
nsel,a,,,157
nsel,a,,,100054
nsel,a,,,100061
nsel,a,,,100150
nsel,a,,,100157
d,all,uy,,,,,,uz,rotx,rotz
allsel,all
numcmp,all
nummrg,all
!工况1 一期恒载
/solu
acel,0,0,9.8
solve |