钢管混凝土拱桥ansys命令流很实用很经典Word文档下载推荐.docx

1、k,9,0.75,0.25k,10,0.5,0a,1,2,3,4,5,6,7,8,9,10secwrite,hl1,sect,1sectype,4,beam,meshhl1建立横梁标准截面b-bk,2,-0.5,1k,3,-0.75,1.25k,4,-0.75,1.5k,5,0.75,1.5k,6,0.75,1.25k,7,0.5,1k,8,0.5a,1,2,3,4,5,6,7,8secwrite,hl2,sect,1sectype,5,beam,meshhl2建立横梁标准截面c-ck,2,-0.5,1.04k,3,-0.75,1.04+0.25k,4,-0.75,1.04+0.5k,5,-0

2、.55,1.04+0.5k,6,-0.55,1.04+0.5+0.35k,7,0.55,1.04+0.5+0.35k,8,0.55,1.04+0.5k,9,0.75,1.04+0.5k,10,0.75,1.04+0.25k,11,0.5,1.04k,12,0.5k,13,-0.2,0.58k,14,-0.2,1.04+0.5-0.3k,15,0.2,1.04+0.5-0.3k,16,0.2,0.58a,1,2,3,4,5,6,7,8,9,10,11,12a,13,14,15,16asba,1,2secwrite,hl3,sect,1sectype,6,beam,meshhl3建立横梁标准截面d

3、-dldele,allkdele,allk,2,-0.5,1.23k,3,-0.75,1.23+0.25k,4,-0.75,1.23+0.5k,5,-0.55,1.23+0.5k,6,-0.55,1.23+0.5+0.35k,7,0.55,1.23+0.5+0.35k,8,0.55,1.23+0.5k,9,0.75,1.23+0.5k,10,0.75,1.23+0.25k,11,0.5,1.23k,13,-0.2,0.38k,14,-0.3,0.38+0.2k,15,-0.3,1.23+0.5-0.1-0.2k,16,-0.2,1.23+0.5-0.1k,17,0.2,1.23+0.5-0.1

4、k,18,0.3,1.23+0.5-0.1-0.2k,19,0.3,0.38+0.2k,20,0.2,0.38a,13,14,15,16,17,18,19,20secwrite,hl4,sect,1sectype,7,beam,meshhl4建立横梁标准截面e-esecwrite,hl5,sect,1sectype,8,beam,meshhl5etdele,1 !将辅助单元删除*定义所有材料特性et,1,beam44 !钢管特性mp,ex,1,2.1e11mp,dens,1,7800mp,prxy,1,0.3n,90000,0,0,30 !参考点et,2,beam44 !钢管内50混凝土特性m

5、p,ex,2,3.5e10mp,dens,2,2600mp,prxy,2,0.1667et,3,beam44 !纵梁30混凝土钢管特性mp,ex,3,3.0e10mp,dens,3,2600mp,prxy,3,0.1667et,4,beam44 !横梁30混凝土钢管特性mp,ex,4,3.0e10mp,dens,4,2600mp,prxy,4,0.1667et,5,beam44 !风撑特性mp,ex,5,2.1e11mp,dens,5,7800mp,prxy,5,0.3et,6,link10 !吊杆特性（钢绞线）mp,ex,6,1.9e11mp,dens,6,7800mp,prxy,6,0.3

6、keyopt,6,3,0 !只拉吊杆et,7,beam44 !盖梁30混凝土特性mp,ex,7,3.0e10mp,dens,7,2600mp,prxy,7,0.1667et,8,beam44 !墩柱40混凝土特性mp,ex,8,3.3e10mp,dens,8,2600mp,prxy,8,0.1667et,9,beam44 !横梁30混凝土特性mp,ex,9,3.0e10mp,dens,9,2600mp,prxy,9,0.1667et,10,link10 !系杆特性（钢绞线）mp,ex,10,2.0e11mp,dens,10,7800mp,prxy,10,0.3keyopt,10,3,0 !只拉

7、系杆et,11,beam44 !承台桩基基础25混凝土特性mp,ex,11,2.85e10mp,dens,11,2500mp,prxy,11,0.1667#开始建立模型0号横梁（端横梁）n,1,2.025n,2,6.25,2.025n,3,6.25+5.1,2.025n,4,6.25+5.1+0.2+1.9/2,2.025 !12.5m位置n,5,6.25+5.1+2.3,2.025n,6,6.25+5.1+2.3+2.23,2.025nsym,x,10,2,6,1type,4r,1real,1mat,4secnum,7单元1到10号为0号横梁e,1,2secnum,6e,2,3secnum,

8、5e,3,4e,4,5secnum,4e,5,6e,1,12e,12,13e,13,14e,14,15e,15,161号端横梁egen,2,100,1,10,1,0,4.7+1.25/2+1.5/2其余端横梁217号egen,17,100,11,20,1,0,5.3+1.518号端横梁egen,2,100,171,180,1,0,4.7+1.25/2+1.5/2cm,hl,elemesel,none*纵梁单元type,3r,2real,2mat,3sectype,10,beam,rect,zlsecdata,0.4,1.2secnum,100号纵梁e,1,101e,2,102e,3,103e,

9、5,105e,12,112e,13,113e,15,115中间纵梁116号*do,i,0,15,1e,101+i*100,201+i*100e,102+i*100,202+i*100e,104+i*100,204+i*100e,112+i*100,212+i*100e,114+i*100,214+i*100*enddo17号纵梁（端横梁）e,1701,1801e,1702,1802e,1703,1803e,1705,1805e,1712,1812e,1713,1813e,1715,1815cm,zl,elem*定义拱轴线*dim,y,array,16 !将纵桥向距离定义为数组*dim,cs,a

10、rray,16 !定义分值cs（12x/l）*dim,z,array,16 !将拱高度定义为数组吊杆上的点首先定义吊杆的y向（顺桥向）坐标y（1）=8.1y（2）=14.9y（3）=21.7y（4）=28.5y（5）=35.3y（6）=42.1y（7）=48.9y（8）=55.7y（9）=62.5风撑上的点首先定义风撑的y向（顺桥向）坐标y（11）=21.7y（12）=28.5y（13）=42.1y（14）=48.9y（15）=59.10号横梁的y向坐标y（16）=2.025fl=25 !矢高ml=1.1k=log（ml+sqrt（ml*ml-1） !k=ch-1（m）*do,i,1,16,1

11、cs（i）=1-y（i）*2/125z（i）=25-fl/（ml-1）*（cosh（k*cs（i）-1） !拱轴高度上述定义的拱轴线可参照姚玲森桥梁工程p1299*dim,ygd,array,23 !定义预工拱度（是在后续工况中反复调算的结果） !初始运行中可将预拱度全部设为0ygd（1）=0.4000e-02+1.29e-02ygd（2）=0.10685e-01+2.22e-02ygd（3）=0.19483e-01+3.03e-02ygd（4）=0.29219e-01+3.69e-02ygd（5）=0.39022e-01+4.23e-02ygd（6）=0.48243e-01+4.65e-02

12、ygd（7）=0.56390e-01+4.95e-02ygd（8）=0.63070e-01+5.15e-02ygd（9）=0.68014e-01+5.21e-02ygd（10）=0.71003e-01+5.15e-02ygd（11）=0.71739e-01+4.95e-02ygd（12）=0.69962e-01+4.65e-02ygd（13）=0.65514e-01+4.23e-02ygd（14）=0.58289e-01+3.69e-02ygd（15）=0.48248e-01+3.04e-02ygd（16）=0.35507e-01+2.24e-02ygd（17）=0.20415e-01+1.2

13、3e-02ygd（18）=0.65770e-01+0.40234e-01ygd（19）=0.69769e-01+0.41215e-01生成各吊杆的竖向节点（z向）ngen,2,10000,104,114,10,0,0,z（1）+ygd（1）ngen,2,10000,204,214,10,0,0,z（2）+ygd（2）ngen,2,10000,304,314,10,0,0,z（3）+ygd（3）ngen,2,10000,404,414,10,0,0,z（4）+ygd（4）ngen,2,10000,504,514,10,0,0,z（5）+ygd（5）ngen,2,10000,604,614,10,

14、0,0,z（6）+ygd（6）ngen,2,10000,704,714,10,0,0,z（7）+ygd（7）ngen,2,10000,804,814,10,0,0,z（8）+ygd（8）ngen,2,10000,904,914,10,0,0,z（9）+ygd（9）ngen,2,10000,1004,1014,10,0,0,z（8）+ygd（10）ngen,2,10000,1104,1114,10,0,0,z（7）+ygd（11）ngen,2,10000,1204,1214,10,0,0,z（6）+ygd（12）ngen,2,10000,1304,1314,10,0,0,z（5）+ygd（13）

15、ngen,2,10000,1404,1414,10,0,0,z（4）+ygd（14）ngen,2,10000,1504,1514,10,0,0,z（3）+ygd（15）ngen,2,10000,1604,1614,10,0,0,z（2）+ygd（16）ngen,2,10000,1704,1714,10,0,0,z（1）+ygd（17） ngen,2,3000,104,114,10,0,0,0 ngen,2,3000,204,214,10,0,0,0 ngen,2,3000,304,314,10,0,0,0 ngen,2,3000,404,414,10,0,0,0 ngen,2,3000,504

16、,514,10,0,0,0 ngen,2,3000,604,614,10,0,0,0 ngen,2,3000,704,714,10,0,0,0 ngen,2,3000,804,814,10,0,0,0 ngen,2,3000,904,914,10,0,0,0 ngen,2,3000,1004,1014,10,0,0,0ngen,2,3000,1104,1114,10,0,0,0ngen,2,3000,1204,1214,10,0,0,0ngen,2,3000,1304,1314,10,0,0,0ngen,2,3000,1404,1414,10,0,0,0ngen,2,3000,1504,151

17、4,10,0,0,0ngen,2,3000,1604,1614,10,0,0,0ngen,2,3000,1704,1714,10,0,0,0ngen,2,6000,104,114,10,0,0,z（1）+ygd（1）ngen,2,6000,204,214,10,0,0,z（2）+ygd（2）ngen,2,6000,304,314,10,0,0,z（3）+ygd（3）ngen,2,6000,404,414,10,0,0,z（4）+ygd（4）ngen,2,6000,504,514,10,0,0,z（5）+ygd（5）ngen,2,6000,604,614,10,0,0,z（6）+ygd（6）ng

18、en,2,6000,704,714,10,0,0,z（7）+ygd（7）ngen,2,6000,804,814,10,0,0,z（8）+ygd（8）ngen,2,6000,904,914,10,0,0,z（9）+ygd（9）ngen,2,6000,1004,1014,10,0,0,z（8）+ygd（10）ngen,2,6000,1104,1114,10,0,0,z（7）+ygd（11）ngen,2,6000,1204,1214,10,0,0,z（6）+ygd（12）ngen,2,6000,1304,1314,10,0,0,z（5）+ygd（13）ngen,2,6000,1404,1414,10

19、,0,0,z（4）+ygd（14）ngen,2,6000,1504,1514,10,0,0,z（3）+ygd（15）ngen,2,6000,1604,1614,10,0,0,z（2）+ygd（16）ngen,2,6000,1704,1714,10,0,0,z（1）+ygd（17）*do,i,1,9,1ngen,2,10000,i*100+4,i*100+14,10,0,0,z（i）+ygd（i） !吊杆节点jj=1*do,i,10,17,1ngen,2,10000,i*100+4,i*100+14,10,0,0,z（i-2*jj）+ygd（i） !jj=jj+1为了对以后吊杆转动进行放松，重复

20、吊杆节点ngen,2,3000,i*100+4,i*100+14,10,0,0,z（i）+ygd（i） !吊杆下节点ngen,2,3000,i*100+4,i*100+14,10,0,0,z（i-2*jj）+ygd（i） !*do,i,1,10,1ngen,2,6000,i*100+4,i*100+14,10,0,0,z（i）+ygd（i） !ngen,2,6000,i*100+4,i*100+14,10,0,0,z（i-2*jj）+ygd（i） !定义吊杆的参数（每根吊杆均采用109根直径7mm高强钢丝）dgarea=109*3.141593*0.25*0.007*0.007dgyl=0.3

21、*1670*1e6 !吊杆除应力采用0.3倍抗拉强度（n/mm*mm）dgl=dgyl*dgareaedg=1.95e11 !吊杆弹模*dim,L0,array,17 !定义吊杆的初始0应力长度数组*dim,L,array,17 !定义吊杆的受力后长度数组*dim,deltL,array,17 !定义吊杆的受力后长度增量数组*dim,istrn,array,17 !定义吊杆的应力数组*dim,fdg,array,17 !定义在吊杆上施加的力数组L0（i）=z（i） !0应力时单元长度L0（18-i）=z（i） !吊杆沿中轴对称 *do,i,1,17,1fdg（i）=dgl !各吊杆施加的预应力

22、type,6mat,6deltL（i）=fdg（i）*L0（i）/（edg*dgarea） !吊杆受力后的增长量L（i）=L0（i）+deltL（i）istrn（i）=deltL（i）/L（i） !吊杆受力后的应变r,i+1000,dgarea,istrn（i） !定义吊杆的实常数特性real,i+1000 !赋予吊杆的实常数特性e,i*100+4+3000,i*100+4+6000 !生成吊杆单元e,i*100+10+4+3000,i*100+10+4+6000 cm,dg,elem*real,1018 !定义一个全为零的实常数便于利用截面特性时，不产生冲突生成风撑上的节点n,20001,0

23、,21.7,z（11）n,20002,0,28.5,z（12）n,20003,0,42.1,z（13）n,20004,0,48.9,z（14）n,20005,0,59.1,z（15）n,20006,0,65.9,z（15）n,20007,0,76.1,z（14）n,20008,0,82.9,z（13）n,20009,0,96.5,z（12）n,20010,0,103.3,z（11）n,20011,12.5,21.7,z（11）n,20012,12.5,28.5,z（12）n,20013,12.5,42.1,z（13）n,20014,12.5,48.9,z（14）n,20015,12.5,59.1,z（15）+ygd（18）n,20016,12.5,65.9,z（15）+ygd（19）n,20017,12.5,76.1,z（14）n,20018,12.5,82.9,z（13）n,20019,12.5,96.5,z（12）n,20020,12.5,103.3,z（11）n,20021,-12.5,21.7,z（11）n,20022,-12.5,28.5,z（12）n,20023,-