1、斜拉桥的算例/prep7!*定义节点*!*塔*!*塔上部*do,i,54.5,33.5,-0.5k,110-i*2,-(82.1495+i)/7.403322,i,0*enddo!*塔下部*!*远朝!*上*do,i,33,0,-0.5k,110-i*2,-15.6213*i/33.5,i,7.4839*i/33.5-7.4839*enddo!*下*do,i,-0.5,-2.5,-0.5k,110-i*2,-2.7978*i/6,i,-i/2-7.4839*enddo!*近朝!*上*do,i,33,0,-0.5k,182-i*2,-15.6213*i/33.5,i,-7.4839*i/33.5+
2、7.4839*enddo!*下*do,i,-0.5,-2.5,-0.5k,182-i*2,-2.7978*i/6,i,i/2+7.4839*enddo!*桥面(配跨为曲线)*k,192,-2,0,-6k,193,0,0,-6k,194,2,0,-6k,195,9,0,-6*do,i,196,205k,i,(i-195)*9+9,0,-6*enddok,206,105,0,-6k,207,-2,0,6k,208,0,0,6k,209,2,0,6k,210,9,0,6*do,i,211,220k,i,(i-210)*9+9,0,6*enddok,221,105,0,6!*背索锚碇*k,222,-3
3、1,0,14k,223,-31,0,-14!*定义单元*!上部塔的截面几何特性*dim,BBS,43 *dim,HHS,43*dim,SSS,43*dim,IIYYS,43*dim,IIZZS,43*do,i,1,43BBS(i)=1.5*(i-1)/42+3HHS(i)=(i-1)/21+3SSS(i)=BBS(i)*HHS(i) IIzzS(i)=BBS(i)*HHS(i)*3/12IIyyS(i)=HHS(i)*BBS(i)*3/12*enddo!*平均*dim,BS,42*dim,HS,42*dim,SS,42*dim,IYYS,42*dim,IZZS,42*do,i,1,42hS(i
4、)=(hhS(i)+hhS(i+1)/2bS(i)=(bbS(i)+bbS(i+1)/2SS(i)=(SSS(i)+SSS(i+1)/2IYYS(i)=(IIYYS(i)+IIYYS(i+1)/2IZZS(i)=(IIZZS(i)+IIZZS(i+1)/2*enddo!下部塔的截面几何特性*dim,BBX,80*dim,HHX,80*dim,SSX,80*dim,IIYYX,80*dim,IIZZX,80*do,i,1,80BBX(i)=2.5HHX(i)=(i-1)*4.7164/79+5SSX(i)=BBX(i)*HHX(i) IIzzX(i)=BBX(i)*HHX(i)*3/12IIyy
5、X(i)=HHX(i)*BBX(i)*3/12*enddo!*平均*dim,BX,79*dim,HX,79*dim,SX,79*dim,IYYX,79*dim,IZZX,79*do,i,1,79BX(i)=(BBX(i)+BBX(i+1)/2HX(i)=(HHX(i)+HHX(i+1)/2SX(i)=(SSX(i)+SSX(i+1)/2IYYX(i)=(IIYYX(i)+IIYYX(i+1)/2IZZX(i)=(IIZZX(i)+IIZZX(i+1)/2*enddo!*定义单元类型-1-42-塔上部*do,i,1,42et,i,beam4mp,ex,i,3.4e10 MP,DENS,i,2.5
6、e3 mp,prxy,i,0.3r,i,SS(i),IZZS(i),IYYS(i),bS(i),hS(i)*enddo!*定义单元类型-43-115-塔下部*do,i,43,114et,i,beam4mp,ex,i,3.4e10 MP,DENS,i,2.5e3 mp,prxy,i,0.3r,i,SX(i-42),IZZX(i-42),IYYX(i-42),bX(i-42),hX(i-42)*enddo!*定义单元类型-115-桥面-1*et,115,shell63MP,EX,115,3.4e10 MP,DENS,115,2.5e3mp,prxy,115,0.3r,115,1,!*定义单元类型-
7、116-桥面-2*et,116,shell63MP,EX,116,3.4e10 MP,DENS,116,2.5e3mp,prxy,116,0.3r,116,4,!*定义单元类型-117-桥面-3*et,117,shell63MP,EX,117,3.4e10 MP,DENS,117,2.5e3mp,prxy,117,0.3r,117,1,!*定义单元类型-118-119-背索*dim,cable_strain_bs,2cable_strain_bs(1)=3.978873577e-3 !2500cable_strain_bs(2)=3.978873577e-3 !2500et,118,link1
8、0MP,EX,118,2e11 MP,DENS,118,7.85e3 mp,prxy,118,0.3R,118,0.031416,cable_strain_bs(1)et,119,link10MP,EX,119,2e11 MP,DENS,119,7.85e3 mp,prxy,119,0.3R,119,0.031416,cable_strain_bs(2)!*定义单元类型-120-129-主索*dim,cable_strain_zs,11cable_strain_zs(1)=2.419155135e-3 !380cable_strain_zs(2)=2.419155135e-3 !380cabl
9、e_strain_zs(3)=2.2918311805e-3 !360cable_strain_zs(4)=2.2918311805e-3 !360cable_strain_zs(5)=2.228169203e-3 !350cable_strain_zs(6)=2.228169203e-3 !350cable_strain_zs(7)=1.909859317e-3 !300cable_strain_zs(8)=1.65521141e-3 !260cable_strain_zs(9)=1.65521141e-3 !260cable_strain_zs(10)=1.40056340e-3 !220
10、cable_strain_zs(11)=1.273239545e-3 !200*do,i,120,130 et,i,link10MP,EX,i,2e11 MP,DENS,i,7.85e3 mp,prxy,i,0.3R,i,0.007853982,cable_strain_zs(i-119) !定义索面积 直径=0.1 米，初始应变*enddo!*定义桥面与塔的连接*et,131,beam4mp,ex,131,3.4e10 MP,DENS,131,2.5e3 mp,prxy,131,0.3r,131,16,21.3333,21.3333,4,4!*定义塔上部-beam4-1*do,i,1,42l
11、,i,i+1lsel,u,type,1,131latt,i,i,i,lesize,all,1,1lmesh,all,allalls*enddo!*塔下部-beam4-2*!*近朝*do,i,43,114l,i,i+1lsel,u,type,1,131latt,i,i,i,lesize,all,1,1lmesh,all,allalls*enddo!*远朝*l,43,116lsel,u,type,1,131latt,43,43,43,lesize,all,1,1lmesh,all,allalls*do,i,116,186l,i,i+1lsel,u,type,1,131latt,i-72,i-72,
12、i-72,lesize,all,1,1lmesh,all,allalls*enddo!*定义索-link10*!*背索*l,9,222lsel,u,type,1,131latt,118,118,118,lesize,all,1,1lmesh,all,allallsl,9,223lsel,u,type,1,131latt,119,119,119,lesize,all,1,1lmesh,all,allalls!*主索*do,i,1,11l,i*4+1,206-il,i*4+1,221-ilsel,u,type,1,131latt,i+119,i+119,i+119,lesize,all,1,1lm
13、esh,all,allalls*enddo!*定义桥面与塔的连接*l,110,193l,182,208lsel,u,type,1,131latt,131,131,131,lesize,all,1,1lmesh,all,allalls!*定义桥面-area-4*wproat,0,-90,0cyl4,-5.6695,-70,76,90,64,118.09989655,0wproat,0,90,0a,188,192,207,191lsel,u,type,1,131aatt,115,115,115,!lesize,all,2,2lesize,213,2,2lesize,214,2,2lesize,21
14、5,2,2lesize,216,2,2lesize,217,4,2lesize,218,2,2lesize,219,4,2allsamesh,alla,192,193,208,207a,193,194,209,208lsel,u,type,1,131aatt,116,116,116,lesize,all,2,1allsamesh,all*do,i,1,12a,i+193,i+194,i+209,i+208*enddolsel,u,type,1,131aatt,117,117,117,lesize,all,2,1allsamesh,all!*活载*!*不计冲击力-验算最大竖向挠度*!sfa,al
15、l,1,PRES,1834!f,196,fy,-114286!f,197,fy,-114286!*f,196,fy,-160000f,197,fy,-160000finish!*加约束*/SOLUd,115,alld,187,alld,188,alld,189,alld,235,uyd,239,uyd,280,uyd,292,uyd,735,uyd,739,uy!*重力加速度*acel,12sfa,all,1,PRES,2567!acel,10solvefinish/POST1 !进入后处理模块 POST1 esel,s,type,118,130etable,link-fx,smisc,1esel,s,type,1,114etable,MXI,smisc,4esel,s,type,1,114etable,MXJ,smisc,10esel,s,type,1,114etable,MYI,smisc,5esel,s,type,1,114etable,MYJ,smisc,11esel,s,type,1,114etable,MZI,smisc,6esel,s,type,1,114etable,MZJ,smisc,12esel,s,type,all
