搜索
fini/clear/title,concrete-filledtubearchbridge/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','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,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,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.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!!!!!!建立横梁标准截面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-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,-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,allaclear,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,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.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,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+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!!其余端横梁2～17号egen,17,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!中间纵梁1～16号*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*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＝（1－2x/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.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,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-02ygd(7)=0.56390e-01+4.95e-02ygd(8)=0.63070e-01+5.15e-02ygd