ansys钢管混凝土拱桥

ansys钢管混凝土拱桥/prep7/title,theanalyseofsteelarchbridgeet,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.2asel,allaptn,alladele,5,8,1allselaadd,allsmrtsize,5amesh,allsecwrite,gg,sect,,1sectype,1,beam,meshsecoffset,cent,,,secread,gg,sect,,meshasel,allaclear,alladele,all,,,1/replotcyl4,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,,,1k,1,-0.5,0k,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,,,1k,1,-0.5,0k,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,,,1k,1,-0.5,0k,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.5,0k,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,,,1ldele,allkdele,allk,1,-0.5,0k,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.5,0k,13,-0.2,0.38k,14,-0.3,0.38k,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,2smrtize,5amesh,allsecwrite,hl4,sect,,1sectype,7,beam,meshsecoffset,cent,,,secread,hl4,sect,,meshasel,allaclear,alladele,all,,,1/replotet,8,beam44et,8,beam44mp,ex,8,,3.5e10k,1,-0.5,0k,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.5,0a,1,2,3,4,5,6,7,8,9,10,11,12smrtize,5amesh,allsecwrite,hl5,sect,,1sectype,8,beam,meshsecoffset,cent,,,secread,hl5,sect,,meshasel,allaclear,alladele,all,,,1/replotetdele,1/replotet,1,beam44mp,ex,1,2.1e11mp,dens,1,7800mp,prxy,1,0.3n,90000,0,0,30mp,dens,2,2600mp,prxy,2,0.1667et,3,beam44mp,ex,3,3.0e10mp,dens,3,2600mp,prxy,3,0.1667et,4,beam44mp,ex,4,3.0e10mp,dens,4,2600mp,prxy,4,0.1667et,5,beam44mp,ex,5,2.1e11mp,dens,5,7800mp,prxy,5,0.3et,6,link10mp,ex,6,1.9e11mp,dens,6,7800mp,prxy,6,0.3keyopt,6,3,0et,7,beam44mp,ex,7,3.0e10mp,dens,7,2600mp,prxy,7,0.1667mp,dens,8,2600mp,prxy,8,0.1667et,9,beam44mp,ex,9,3.0e10mp,dens,9,2600mp,prxy,9,0.1667et,10,link10mp,ex,10,2.0e11mp,dens,10,7800mp,prxy,10,0.3keyopt,10,3,0et,11,beam44mp,ex,11,2.85e10mp,dens,11,2500mp,prxy,11,0.1667n,90000,0,0,30n,1,0,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.025n,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,7e,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,16egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2egen,17,100,11,20,1,,,,,,0,5.3+1.5egen,2,100,171,180,1,,,,,,0,4.7+1.5/2+1.25/2cm,hl,elemesel,nonetype,3r,2real,2mat,3sectype,10,beam,rect,zlsecdata,0.4,1.2secnum,10e,1,101e,2,102e,3,103e,5,105e,12,112e,13,113e,15,115*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*enddoe,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*dim,z,array,16y(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.5y(11)=21.7y(12)=28.5y(13)=42.1y(14)=48.9y(15)=59.1y(16)=2.025fl=25ml=1.1k=log(ml+sqrt(mldim,ygd,array,23ygd(1)=0.4004e-02+1.29e-02dim,ygd,array,23ygd(1)=0.4004e-02+1.29e-02*do,i,1,16,1cs(i)=1-y(i)*2/125z(i)=25-fl/(ml-1)*(cosh(k*cs(i))-1)*enddoygd(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)=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.69-02ygd(15)=0.48248e-01+3.04e-02ygd(16)=0.35507e-01+2.24e-02ygd(17)=0.35507e-01+1.23e-02ygd(18)=0.65770e-01+0.40234e-01ygd(19)=0.69769e-01+0.41215e-01ngen,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,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)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,0ngen,2,3000,204,214,10,0,0,0ngen,2,3000,304,314,10,0,0,0ngen,2,3000,404,414,10,0,0,0ngen,2,3000,504,514,10,0,0,0ngen,2,3000,604,614,10,0,0,0ngen,2,3000,704,714,10,0,0,0ngen,2,3000,804,814,10,0,0,0ngen,2,3000,904,914,10,0,0,0ngen,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,3000,1804,1814,10,0,0,0ngen,2,3000,1904,1914,10,0,0,0ngen,2,3000,2004,2014,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)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,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,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)dgarea=109*3.1415926*0.25*0.007*0.007dgy1=0.3*1670*1e6dgl=dgyl*dgareaedg=1.95e11*dim,l0,arry,17*dim,l,array,17*dim,deltl,array,17*dim,istrn,array,17*dim,fdg,array,17*do,i,1,9,1l0(i)=z(i)l0(18-i)=z(i)*enddofdg(1)=dglfdg(2)=dglfdg(3)=dglfdg(4)=dglfdg(5)=dglfdg(6)=dglfdg(7)=dglfdg(8)=dglfdg(9)=dglfdg(10)=dglfdg(11)=dglfdg(12)=dglfdg(13)=dglfdg(14)=dglfdg(15)=dglfdg(16)=dglfdg(17)=dgltype,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+1000e,i*100+4+3000,i*100+4+6000e,i*100+10+4+3000,i*100+10+4+6000*enddocm,dg,elemesel,noneallselreal,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,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,-12.5,42.1,z(13)n,20024,-12.5,48.9,z(14)n,20025,-12.5,59.1,z(15)+ygd(18)n,20026,-12.5,65.9,z(15)+ygd(19)n,20027,-12.5,76.1,z(14)n,20028,-12.5,82.9,z(13)n,20029,-12.5,96.5,z(12)n,20030,-12.5,103.3,z(11)n,20000,0,62.5,z(15)type,5mat,5sectype,500,beam,ctube,windcsecdata,0.486,0.50secnum,500e,20002,10404,90000e,20002,10414,90000e,20002,10304,90000e,20002,10314,90000e,20004,10704,90000e,20004,10714,90000e,20004,10604,90000e,20004,10614,90000e,20000,20015,90000e,20000,20016,90000e,20000,20025,90000e,20000,20026,90000e,20007,11104,90000e,20007,11114,90000e,20007,11204,90000e,20007,11214,90000e,20009,11404,90000e,20009,11414,90000e,20009,11304,90000e,20009,11314,90000cm,fc,elemesel,nonen,30000,12.5,0,0n,30001,12.5,125,0n,30002,12.5,2.025,z(16)n,30003,12.5,125-2.025,z(16)n,40000,-12.5,0,0n,40001,-12.5,125,0n,40002,-12.5,2.025,z(16)n,40003,-12.5,125-2.025,z(16)*dim,zgzx,array,23*dim,ygzx,array,23zgzx(1)=30000zgzx(2)=30002zgzx(3)=10104zgzx(4)=10204zgzx(5)=10304zgzx(6)=10404zgzx(7)=10504zgzx(8)=10604zgzx(9)=10704zgzx(10)=10804zgzx(11)=20015zgzx(12)=10904zgzx(13)=20016zgzx(14)=11004zgzx(15)=11104zgzx(16)=11204zgzx(17)=11304zgzx(18)=11404zgzx(19)=11504zgzx(20)=11604zgzx(21)=11704zgzx(22)=30003zgzx(23)=30001ygzx(1)=40000ygzx(2)=40000ygzx(3)=10114ygzx(4)=10214ygzx(5)=10314ygzx(6)=10414ygzx(7)=10514ygzx(8)=10614ygzx(9)=10714ygzx(10)=10814ygzx(11)=20025ygzx(12)=10914ygzx(13)=20026ygzx(14)=11014ygzx(15)=11114ygzx(16)=11214ygzx(17)=11314ygzx(18)=11414ygzx(19)=11514ygzx(20)=11614ygzx(21)=11714ygzx(22)=40003ygzx(23)=40001n,90001,12.5,0,30n,90002,-12.5,0,30type,1mat,1r,2000secnum,1*do,i,1,22,1e,zgzx(i),zgzx(i+1),90001e,ygzx(i),ygzx(i+1),90002*enddocm,gzgg,elemesel,nonetype,2mat,2r,2001secnum,2*do,i,1,22,1e,zgzx(i),zgzx(i+1),90001e,ygzx(i),ygzx(i+1),90002*enddocm,gzhnt,elemesel,nonetype,9r,6real,6mat,9sectype,70,beam,rect,dizuosecdata,1.4,4secnum,70e,30001,30003,90001e,40001,40003,90002e,30000,30002,90001e,40000,40002,90002cm,gzdz,elemesel,nonengen,2,2000,1,16,1,0,-2.205,0ngen,2,100,1801,1816,1,0,2.025,0type,7mat,7r,7real,7sectype,60,beam,rect,gailsecdata,4,1.6secnum,60n,90005,0,0,20n,90006,0,125,20e,1901,1902,90006e,1902,1903,90006e,1903,1904,90006e,1904,1905,90006e,1905,1906,90006e,1901,1912,90006e,1912,1913,90006e,1913,1914,90006e,1914,1915,90006e,1915,1916,90006e,2001,2002,90005e,2002,2003,90005e,2003,2004,90005e,2004,2005,90005e,2005,2006,90005e,2001,2012,90005e,2012,2013,90005e,2013,2014,90005e,2014,2015,90005e,2015,2016,90005cm,gailiang,elemesel,nonen,90003,14.5,0,0n,90004,-14.5,0,0type,8mat,8r,8real,8sectype,50,beam,hrec,dunsecdata,3,4,0.8,0.8,0.8,0.8ngen,2,10000,2004,2014,10,0,0,-10ngen,2,10000,1904,1914,10,0,0,-10secnum,50e,2004,12004,90003e,2014,12014,90004e,1904,11904,90003e,1914,11914,90004cm,dunzhu,elemesel,nonetype,11mat,11r,1sectype,40,beam,rect,chengtaisecdata,12.7,9r,11real,11secnum,40ngen,2,20000,2004,2014,10,0,0,-11.5ngen,2,20000,1904,1914,10,0,0,-11.5ngen,2,30000,2004,2014,10,0,0,-13ngen,2,30000,1904,1914,10,0,0,-13e,12004,22004,90003e,11904,21904,90003e,12014,22014,90004e,11914,21914,90004e,22004,32004,90003e,21904,31904,90003e,22014,32014,90004e,21914,31914,90004cm,chtai,elemesel,nonetype,9mat,9r,9real,9n,90007,12.5,0,10n,90008,-12.5,0,10sectype,90,beam,rect,hxlsecdata,3,3secnum,90e,22004,22014,90007e,21904,21914,90008cm,hxl,elemesel,nonen,50000,6.15+1.85,2.7,-13n,50001,6.15+1.85,-2.7,-13n,50002,6.15+1.85+4.5,2.7,-13n,50003,6.15+1.85+4.5,-2.7,-13n,50004,6.15+12.7-1.85,2.7,-13n,50005,6.15+12.7-18.5,-2.7,-13nsel,s,,,50000,50005,1nsym,x,10,50000,50005,1allselngen,2,100,50000,50015,1,0,125,0ngen,10,1000,50000,51000,1,0,0,-2.5esel,nonetype,11mat,11r,11real,11sectype,20,beam,csolid,piersecdata,0.9secnum,20*do,i,0,8,1*do,k,0,5,1e,50000+k+idim,m,array,5m(1)=8dim,m,array,5m(1)=8m(2)=15m(3)=20m(4)=25m(5)=4000*dim,h,array,5h(1)=5h(2)=5h(3)=5h(4)=5h(5)=10*dim,c,array,4c(1)=m(1)*h(1)c(2)=m(2)*(h(1)+h(2))e,50000+k+10+i*1000,50000+k+10+(i+1)*1000,90003e,50000+k+100+i*1000,50000+k+100+(i+1)*1000,90003e,50000+k+110+i*1000,50000+k+110+(i+1)*1000,90003*enddo*enddocm,zhuang,elemesel,noned=1.8kf=0.9k0=1+1/dbl=0.6ll=3.6hl=3*(d+1)b=kf*k0*k*dc(3)=m(3)enddocm,spring1,elemenddocm,spring1,elemesel,nonengen,2,10000,53000,53115,1,1,0,0ngen,2,20000,53000,53115,1,0,1,0ngen,2,30000,53000,53115,1,0,0,1type,20real,30*do,i,0,5,1e,53000+i,53000+i+10000e,53000+i+10,53000+i+10+10000e(1)=0.5*c(1)*h(1)*b*1e6*do,i,2,4,1e(i)=0.5*(c(i-1)+c(i))*h(i)*b*1e6*enddoet,20,combin14r,20,e(1)r,30,e(2)r,40,e(3)r,50,e(4)r,60,m(5)*3.1415926*b*b/4*1e6ngen,2,10000,51000,51115,1,1,0,0ngen,2,20000,51000,51115,1,0,1,0ngen,2,30000,51000,51115,1,0,0,1type,20real,20*do,i,0,5,1e,51000+i,51000+i+10000e,51000+i+10,51000+i+10+10000e,51100+i,51100+i+10000e,51100+i+10,51100+i+10+10000e,51000+i,51000+i+20000e,51000+i+10,51000+i+10+20000e,51100+i,51100+i+20000e,51100+i+10,51100+i+10+20000e,53100+i,53100+i+10000e,53100+i+10,53100+i+10+10000e,53000+i,53000+i+20000e,53000+i+10,53000+i+10+20000e,53100+i,53100+i+20000e,53100+i+10,53100+i+10+20000*enddocm,spring?,elemesel,nonengen,2,10000,55000,55115,14,0,0ngen,2,20000,55000,55115,1,0,1,0type,20real,40*do,i,0,5,1e,55000+i,55000+i+10000e,55000+i+10,55000+i+10+10000e,55100+i,55100+i+10000e,55100+i+10,55100+i+10+10000e,55000+i,55000+i+20000e,55000+i+10,55000+i+10+20000e,55100+i,55100+i+20000e,55100+i+10,55100+i+10+20000*enddocm,spring3,elemesel,nonengen,2,10000,57000,57115,14,0,0ngen,2,20000,57000,57115,l,0,l,0type,20real,50*do,i,0,5,1e,57000+i,57000+i+10000e,57000+i+10,57000+i+10+10000e,57100+i,57100+i+10000e,57100+i+10,57100+i+10+10000e,57000+i,57000+i+20000e,57000+i+10,57000+i+10+20000e,57100+i,57100+i+20000e,57100+i+10,57100+i+10+20000*enddocm,spring4,elemesel,nonengen,2,10000,58000,58115,1,1,0,0ngen,2,20000,58000,58115,1,0,1,0type,20real,60*do,i,0,5,1e,58000+i,58000+i+10000e,58000+i+10,58000+i+10+10000e,58100+i,58100+i+10000e,58100+i+10,58100+i+10+10000e,58000+i,58000+i+20000e,58000+i+10,58000+i+10+20000e,58100+i,58100+i+20000e,58100+i+10,58100+i+10+20000*enddocm,spring5,elemesel,nonengen,2,10000,59000,59115,1,1,0,0ngen,2,20000,59000,59115,1,0,1,0ngen,2,30000,59000,59115,1,0,0,-1type,20real,60*do,i,0,5,1e,59000+i,59000+i+10000e,59000+i+10,59000+i+10+10000e,59100+i,59100+i+10000e,59100+i+10,59100+i+10+10000e,59000+i,59000+i+20000e,59000+i+10,59000+i+10+20000e,59100+i,59100+i+20000e,59100+i+10,59100+i+10+20000e,59000+i,59000+i+30000e,59000+i+10,59000+i+10+30000e,59100+i,59100+i+30000e,59100+i+10,59100+i+10+30000*enddocm,spring6,elemesel,noneallselnsel,s,loc,z,-50,-0.01nsel,r,loc,y,-10,10cm,zdnode,nodeesln,s,,allcmzdelem,elemallselnsel,u,,,zdnodensym,y,100000,allallselesel,u,,,zdelemesym,,100000,allallselxgarea=8*27*7*3.1415926*0.25*0.005*0.005exg=1.95e11*dim,fxg,array,2*dim,llast,array,2*dim,listrn,array,2*dim,detl,array,2xgyl=0.3*1860*1e6xgl=xgyl*xgareafxg(1)=xglfxg(2)=xgllzero=250*do,i,1,2,1detl(i)=lzero*fxg(i)/(exg*xgarea)llast(i)=lzero+detl(i)listrn(i)=detl(i)/llast(i)r,i+100,xgarea,listrn*enddotype,10mat,10real,101e,30001,130001real,102e,40001,140001cm,xigan,elemesel,noneallselcp,1,ux,2004,30000cp,1,uy,2004,30000cp,next,uz,2004,30000cp,next,roty,2004,30000cp,next,rotz,2004,30000cp,next,ux,1904,30001cp,next,uy,1904,30001cp,next,uz,1904,30001cp,next,roty,1904,30001cp,next,rotz,1904,30001cp,next,ux,2014,40000cp,next,uy,2014,40000cp,next,uz,2014,40000cp,next,roty,2014,40000cp,next,rotz,2014,40000cp,next,ux,1914,40001cp,next,uy,1914,40001cp,next,uz,1914,40001cp,next,roty,1914,40001cp,next,rotz,1914,40001cp,next,ux,101914,140001cp,next,uy,101914,140001cp,next,uz,101914,140001cp,next,roty,101914,140001cp,next,rotz,101914,140001cp,next,ux,101904,130001cp,next,uy,101904,130001cp,next,uz,101904,130001cp,next,roty,101904,130001cp,next,rotz,101904,130001allselnsel,s,,,30000,130000,100000nsel,a,,,40000,140000,100000nummrg,nodeallsel*do,i,1,3,1cp,next,all,1800+i,1900+i*enddoallsel*do,i,5,6,1cp,next,all,1800+i,1900+i*enddoallsel*do,i,1,2,1cp,next,all,1800+i+11,1900+i+11cp,next,all,i+11,2000+i+11,i+100000+11*enddoallsel*do,i,4,5,1cp,next,all,1800+i+11,1900+i+11cp,next,all,i+11,2000+i+11,i+100000+11*enddoallsel*do,i,1,3,1cp,next,all,i,2000+i,100000+i*enddoallsel*do,i,5,6,1cp,next,all,i,2000+i,100000+i*enddoallsel*do,i,1,3,1cp,next,alli,i,101800+i,101900+i*enddoallsel*do,i,5,6,1cp,next,all,i,101800+i,101900+i*enddoallsel*do,i,1,2,1cp,next,all,i,101800+i+11,101900+i+11*enddoallsel*do,i,4,5,1cp,next,all,i,101800+i+11,101900+i+11*enddoallselnsel,s,loc,z,-13nsel,r,loc,y,-3,3nsel,r,loc,x,8,20cp,next,all,allallselnsel,s,loc,z,-13nsel,r,loc,y,-3,3nsel,r,loc,x,-8,-20cp,next,all,allallselnsel,s,loc,z,-13nsel,r,loc,y,122,128nsel,r,loc,x,-8,-20cp,next,all,allallselnsel,s,loc,z,-13nsel,r,loc,y,122,128nsel,r,loc,x,8,20cp,next,all,allallselnsel,s,loc,z,-13nsel,r,loc,y,-122,-128nsel,r,loc,x,8,20cp,next,all,allallselnsel,s,loc,z,-13nsel,r,loc,y,-122,-128nsel,r,loc,x,-8,-20cp,next,all,allallsel*do,i,1,17,1cp,next,ux,i*100+4,i*100+4+3000cp,next,ux,i*100+4+6000,i*100+4+10000cp,next,uy,i*100+4,i*100+4+3000cp,next,uy,i*100+4+6000,i*100+4+10000cp,next,uz,i*100+4,i*100+4+3000cp,next,uz,i*100+4+6000,i*100+4+10000cp,next,ux,i*100+4+10,i*100+4+3000+10cp,next,ux,i*100+4+6000+10,i*100+4+10000+10cp,next,uy,i*100+4+10,i*100+4+3000+10cp,next,uy,i*100+4+6000+10,i*100+4+10000+10cp,next,uz,i*100+4+10,i*100+4+3000+10

cp,next,uz,i*100+4+6000+10,i*100+4+10000+10cp,next,ux,i*100+4+100000,i*100+4+3000+100000cp,next,ux,i*100+4+6000+100000,i*100+4+10000+100000cp,next,uy,i*100+4+100000,i*100+4+3000+100000cp,next,uy,i*100+4+6000+100000,i*100+4+10000+100000cp,next,uz,i*100+4+100000,i*100+4+3000+100000cp,next,uz,i*100+4+6000+100000,i*100+4+10000+100000cp,next,ux,i*100+4+100000+10,i*100+4+3000+100000+10cp,next,ux