5-5.txt

是用apdl语言编的

/prep7
/VIEW, 1 ,1,1,1 
et,1,shell63
et,2,link8
r,1,0.25        !顶层与底层的厚度
r,2,0.6         !腹板的厚度
r,3,1            !V形支撑厚度
r,4,0.02,-0.005   !预应力钢筋的面积初始与初始应变
mp,ex,1,3.25e10
mp,prxy,1,0.2
mp,dens,1,2700
mp,ex,2,207e9
mp,prxy,2,0.3
mp,dens,2,7800

!桥面关键点  
k,1
k,2,4 
k,3,10   
k,4,16  
k,5,20  
k,6,24  
k,7,30   
k,8,36
k,9,40  
SAVE

!箱形梁底面关键点
k,10,4,-1.8   
k,11,10,-1.8 
k,12,16,-1.8
k,13,24,-1.8
k,14,30,-1.8 
k,15,36,-1.8
SAVE
SAVE

kgen,2,1,15,1,,,-16
kgen,2,1,15,1,,,-85
kgen,2,1,15,1,,,-154 
kgen,2,1,15,1,,,-170
SAVE

k,101,,,-36   
k,102,4,,-36
k,103,10,,-36  
k,104,16,,-36 
k,105,20,,-36 
k,106,24,,-36 
k,107,30,,-36  
k,108,36,,-36   
k,109,40,,-36
kgen,2,102,104,1,,-3,,8 
kgen,2,106,108,1,,-3,,7 
SAVE

kgen,2,101,115,1,,,-12,15   
kgen,2,116,130,1,,,-12,15
kgen,2,131,145,1,,,-50,15
kgen,3,146,160,1,,,-12,15
SAVE

!顶层面建模并网格划分
*do,i,1,8,1
a,i,i+15,i+16,i+1
*enddo
*do,i,16,23,1
a,i,i+85,i+86,i+1
*enddo
*do,i,101,108,1
a,i,i+15,i+16,i+1
*enddo
*do,i,116,123,1
a,i,i+15,i+16,i+1
*enddo
*do,i,31,38,1
a,i,i+100,i+101,i+1
*enddo
*do,i,31,38,1
a,i,i+115,i+116,i+1
*enddo
*do,i,146,153,1
a,i,i+15,i+16,i+1
*enddo
*do,i,161,168,1
a,i,i+15,i+16,i+1
*enddo
*do,i,46,53,1
a,i,i+130,i+131,i+1
*enddo
*do,i,46,53,1
a,i,i+15,i+16,i+1
*enddo
save
lesize,all,4
aatt,1,1,1
amesh,all


!lsel,s,loc,x,-0.1,0.1
!lsel,a,loc,x,41.9,42.1
!lesize,all,,,10
!allsel
!lplot
!lsel,s,loc,z,-0.1,0.1
!lsel,a,loc,z,-165,-166
!lesize,all,,,2
!lplot

asel,u,loc,y,-0.1,0.1	!排除顶面面元素
aplot

!底层
*do,i,10,11,1
a,i,i+15,i+16,i+1
*enddo
*do,i,13,14,1
a,i,i+15,i+16,i+1
*enddo
*do,i,25,26,1
a,i,i+85,i+86,i+1
*enddo
*do,i,28,29,1
a,i,i+85,i+86,i+1
*enddo
*do,i,110,111,1
a,i,i+15,i+16,i+1
*enddo
*do,i,113,114,1
a,i,i+15,i+16,i+1
*enddo
*do,i,125,126,1
a,i,i+15,i+16,i+1
*enddo
*do,i,128,129,1
a,i,i+15,i+16,i+1
*enddo
save
*do,i,40,41,1
a,i,i+100,i+101,i+1
*enddo
*do,i,43,44,1
a,i,i+100,i+101,i+1
*enddo
*do,i,40,41,1
a,i,i+115,i+116,i+1
*enddo
*do,i,43,44,1
a,i,i+115,i+116,i+1
*enddo
*do,i,155,156,1
a,i,i+15,i+16,i+1
*enddo
*do,i,158,159,1
a,i,i+15,i+16,i+1
*enddo
*do,i,170,171,1
a,i,i+15,i+16,i+1
*enddo
*do,i,173,174,1
a,i,i+15,i+16,i+1
*enddo
*do,i,55,56,1
a,i,i+130,i+131,i+1
*enddo
*do,i,58,59,1
a,i,i+130,i+131,i+1
*enddo
*do,i,55,56,1
a,i,i+15,i+16,i+1
*enddo
*do,i,58,59,1
a,i,i+15,i+16,i+1
*enddo
save

lesize,all,4
aatt,1,1,1
amesh,all

asel,u,loc,z,-170.1,0.1
/pnum,kp,1
kplot
aplot
save


!腹板
*do,i,2,4,1
a,i,i+8,i+23,i+15
*enddo
*do,i,6,8,1
a,i,i+7,i+22,i+15
*enddo
*do,i,17,19,1
a,i,i+8,i+93,i+85
*enddo
*do,i,21,23,1
a,i,i+7,i+92,i+85
*enddo
*do,i,102,104,1
a,i,i+8,i+23,i+15
*enddo
*do,i,106,108,1
a,i,i+7,i+22,i+15
*enddo
*do,i,117,119,1
a,i,i+8,i+23,i+15
*enddo
*do,i,121,123,1
a,i,i+7,i+22,i+15
*enddo
*do,i,32,34,1
a,i,i+8,i+108,i+100
*enddo
*do,i,36,38,1
a,i,i+7,i+107,i+100
*enddo
*do,i,32,34,1
a,i,i+8,i+123,i+115
*enddo
*do,i,36,38,1
a,i,i+7,i+122,i+115
*enddo
*do,i,147,149,1
a,i,i+8,i+23,i+15
*enddo
*do,i,151,153,1
a,i,i+7,i+22,i+15
*enddo
*do,i,162,164,1
a,i,i+8,i+23,i+15
*enddo
*do,i,166,168,1
a,i,i+7,i+22,i+15
*enddo
*do,i,47,49,1
a,i,i+8,i+138,i+130
*enddo
*do,i,51,53,1
a,i,i+7,i+137,i+130
*enddo
*do,i,47,49,1
a,i,i+8,i+23,i+15
*enddo
*do,i,51,53,1
a,i,i+7,i+22,i+15
*enddo
save

aplot
!!先对y方向杆设置划分数目为1和2份，然后在对Z!方向的杆设置单元长度为4
lsel,s,loc,z,-0.1,0.1
lsel,a,loc,z,-17,-15
lsel,a,loc,z,-86,-84
lsel,a,loc,z,-155,-153
lsel,a,loc,z,-171,-169
lesize,all,,,2

lsel,all
lsel,s,loc,z,-37,-35
lsel,a,loc,z,-61,-59
lsel,a,loc,z,-111,-109
lsel,a,loc,z,-135,-133
lesize,all,,,2

!esize,all,4
aatt,1,2,1
amesh,all
save

!!!!!!!!!!!!!!!!!!!!!!!!!!!!预应力杆设置
lsel,s,loc,x,3.9,4.1
lsel,a,loc,x,9,11
lsel,a,loc,x,15,17
lsel,a,loc,x,23,25
lsel,a,loc,x,29,31
lsel,a,loc,x,35,37
lsel,u,loc,y,-1.7,-0.1
lplot

lesize,all,4
latt,2,4,2
lmesh,all
save



!V型支座
asel,u,loc,y,-3.1,0.1
aplot
k,1001,4,-15,-48
k,1002,16,-15,-48
k,1003,24,-15,-48
k,1004,36,-15,-48
k,2001,4,-15,-122
k,2002,16,-15,-122
k,2003,24,-15,-122
k,2004,36,-15,-122

a,110,112,1002,1001
a,112,113,1003,1002
a,113,115,1004,1003
a,140,142,1002,1001
a,142,143,1003,1002
a,143,145,1004,1003


a,155,157,2002,2001
a,157,158,2003,2002
a,158,160,2004,2003
a,185,187,2002,2001
a,187,188,2003,2002
a,188,190,2004,2003

save

lesize,all,4
aatt,1,3,1
amesh,all
save
allsel
nummrg,all
numcmp,all
/eshape,1
eplot

!施加约束
nsel,s,loc,z,-0.1,0.1
nsel,a,loc,z,-171,-169
nsel,r,loc,y,-1.9,-1.7
nplot
d,all,uy
allsel
nsel,s,loc,y,-16,-14
nplot
d,all,all
allsel
gplot


/solu
acel,,9.8
pstres,on
allsel
solve
finish

/post1
plnsol,u,y
/eshape,0
/DSCALE,1,user,10
ETABLE, ,LS,1   
PLETAB,LS1,NOAV 
fini

!模态求解
/solu
antype,2
upcoord,1,on
pstres,on
modopt,subsp,10
mxpand,10,,,1
allsel
solve
FINISH  

!后处理
/POST1  
SET,LIST
SET,FIRST   
PLDISP,0
!SET,NEXT
!PLDISP,0
!SET,NEXT
!PLDISP,0
!finish