钢管混凝土拱桥ansys命令流很实用很经典.docx

1、钢管混凝土拱桥ansys命令流很实用很经典fini/clear/title,concrete-filled tube arch bridge/prep7截面1et,1,82cyl4,0.6,0.6,0.6,0.586cyl4,0.6,2.4,0.6,0.586rectng,0.35,0.364,0.8,2.2rectng,0.85,0.836,0.8,2.2allselaptn,alladele,5,8,1allselaadd,allsmrtsize,5amesh,allsecwrite,gg,sect,1sectype,1,beam,meshsecoffset,cent,secread,gg

2、,sect,meshasel,allaclear,alladele,all,1截面2cyl4,0.6,0.6,0.586cyl4,0.6,2.4,0.586rectng,0.364,0.836,0.8,2.2allselaadd,allsmrtsize,5amesh,allsecwrite,hnt,sect,1sectype,2,beam,meshsecoffset,cent,secread,hnt,sect,meshasel,allaclear,alladele,all,1建立横梁标准截面a-ak,1,-0.5k,2,-0.75,0.25k,3,-0.75,0.5k,4,-0.55,0.5k

3、,5,-0.55,0.85k,6,0.55,0.85k,7,0.55,0.5k,8,0.75,0.5k,9,0.75,0.25k,10,0.5,0a,1,2,3,4,5,6,7,8,9,10smrtsize,5amesh,allsecwrite,hl1,sect,1sectype,4,beam,meshsecoffset,cent,secread,hl1,sect,meshasel,allaclear,alladele,all,1建立横梁标准截面b-bk,1,-0.5k,2,-0.5,1k,3,-0.75,1.25k,4,-0.75,1.5k,5,0.75,1.5k,6,0.75,1.25k,

4、7,0.5,1k,8,0.5a,1,2,3,4,5,6,7,8smrtsize,5amesh,allsecwrite,hl2,sect,1sectype,5,beam,meshsecoffset,cent,secread,hl2,sect,meshasel,allaclear,alladele,all,1建立横梁标准截面c-ck,1,-0.5k,2,-0.5,1.04k,3,-0.75,1.04+0.25k,4,-0.75,1.04+0.5k,5,-0.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

5、.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,2smrtsize,5amesh,allsecwrite,hl3,sect,1sectype,6,beam,meshsecoffset,cent,secread,hl3,sect,meshasel,allaclear,alladele,all,1

6、建立横梁标准截面d-dldele,allkdele,allk,1,-0.5k,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,12,0.5k,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

7、-0.1k,17,0.2,1.23+0.5-0.1k,18,0.3,1.23+0.5-0.1-0.2k,19,0.3,0.38+0.2k,20,0.2,0.38a,1,2,3,4,5,6,7,8,9,10,11,12a,13,14,15,16,17,18,19,20asba,1,2smrtsize,5amesh,allsecwrite,hl4,sect,1sectype,7,beam,meshsecoffset,cent,secread,hl4,sect,meshasel,allaclear,alladele,all,1建立横梁标准截面e-ek,1,-0.5k,2,-0.5,1.23k,3,-

8、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,12,0.5a,1,2,3,4,5,6,7,8,9,10,11,12smrtsize,5amesh,allsecwrite,hl5,sect,1sectype,8,beam,meshsecoffset,cent,secread,hl5,sect,meshasel,allac

9、lear,alladele,all,1etdele,1 !将辅助单元删除*!定义所有材料特性et,1,beam44 !钢管特性mp,ex,1,2.1e11mp,dens,1,7800mp,prxy,1,0.3n,90000,0,0,30 !参考点et,2,beam44 !钢管内50混凝土特性mp,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,d

10、ens,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.3keyopt,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,

11、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 !只拉系杆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+

12、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,5e,3,4e,4,5secnum,4e,5,6secnum,7e,1,12secnum,6e,12,13secnum,5e,13,14e,14,15secnum,4e,15,16!1号端横梁egen,2,100,1,10,1,0,4.7+1.25/2+1.5/2!其余端横梁217号egen,17,

13、100,11,20,1,0,5.3+1.5!18号端横梁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,10!0号纵梁e,1,101e,2,102e,3,103e,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*10

14、0e,112+i*100,212+i*100e,114+i*100,214+i*100*enddo!17号纵梁（端横梁）e,1701,1801e,1702,1802e,1703,1803e,1705,1805e,1712,1812e,1713,1813e,1715,1815cm,zl,elemesel,none*!定义拱轴线*dim,y,array,16 !将纵桥向距离定义为数组*dim,cs,array,16 !定义分值cs（12x/l）*dim,z,array,16 !将拱高度定义为数组!吊杆上的点!首先定义吊杆的y向（顺桥向）坐标y(1)=8.1y(2)=14.9y(3)=21.7y(4

15、)=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.1!0号横梁的y向坐标y(16)=2.025fl=25 !矢高ml=1.1k=log(ml+sqrt(ml*ml-1) !k=ch-1(m)*do,i,1,16,1cs(i)=1-y(i)*2/125z(i)=25-fl/(ml-1)*(cosh(k*cs(i)-1) !拱轴高度*enddo!上述定义的拱轴线可参照姚玲森桥梁工程p1299*dim,yg

16、d,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-02ygd(7)=0.56390e-01+4.95e-02ygd(8)=0.63070e-01+5.15e-02ygd(9)=0.68014e-01+5.21e-02ygd(10)

17、=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.23e-02ygd(18)=0.65770e-01+0.40234e-01ygd(19)=0.69769e-01+0.41215e-01!生成各吊杆的竖向节点（z向）ngen,2

18、,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,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)+

19、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)ngen,2,10000,1404,1414,10,0,0,z(4)+ygd(14)ngen,2,10000,1504,1514,10,0,0,z(3)+ygd(15)nge

20、n,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,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

21、,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,1514,10,0,0,0ngen,2,3000,1604,1614,10,0,0,0ngen,2,3000,1704,1714,10,0,0,0ngen,2,6000,104,1

22、14,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)ngen,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,9

23、04,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,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,

24、z(2)+ygd(16)ngen,2,6000,1704,1714,10,0,0,z(1)+ygd(17)!*do,i,1,9,1!ngen,2,10000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i) !吊杆节点!*enddo!jj=1!*do,i,10,17,1!ngen,2,10000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!为了对以后吊杆转动进行放松，重复吊杆节点!*do,i,1,9,1!ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i)+ygd(

25、i) !吊杆下节点!*enddo!jj=1!*do,i,10,17,1!ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!*do,i,1,10,1!ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i) !吊杆节点!*enddo!jj=1!*do,i,10,17,1!ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i) !吊杆节点!jj=jj+1!*enddo!定义吊杆的参数（每根吊杆均采用109根直径7m

26、m高强钢丝）dgarea=109*3.141593*0.25*0.007*0.007dgyl=0.3*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 !定义在吊杆上施加的力数组*do,i,1,9,1L0(i)=z(i) !0应力

27、时单元长度L0(18-i)=z(i) !吊杆沿中轴对称 *enddo*do,i,1,17,1fdg(i)=dgl !各吊杆施加的预应力*enddotype,6mat,6*do,i,1,17,1deltL(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) !定义吊杆的实常数特性*enddo*do,i,1,17,1real,i+1000 !赋予吊杆的实常数特性e,i*100+4+3000,i*100+4+6000 !生成

28、吊杆单元e,i*100+10+4+3000,i*100+10+4+6000 *enddocm,dg,elemesel,noneallsel!*real,1018 !定义一个全为零的实常数便于利用截面特性时，不产生冲突!*!生成风撑上的节点n,20001,0,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

29、,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,-