ANSYS空间梁格法分析连续斜梁桥命令流
连续梁桥ansys命令流(仅供参考)
!!连续梁桥/prep7et,1,4 !!!!定义梁单元et,2,21 !!!!定义mass21单元!!定义粱材料!!泊松比!!密度mp,ex,2,3.45e10 !!直线段梁材料和1M段梁材料mp,nuxy,2,0.2mp,dens,2,3302.153125mp,ex,3,3.45e10mp,nuxy,3,0.2mp,dens,3,3301.658695mp,ex,4,3.45e10mp,nuxy,4,0.2mp,dens,4,3299.906778mp,ex,5,3.45e10mp,nuxy,5,0.2mp,dens,5,3298.327219mp,ex,6,3.45e10mp,nuxy,6,0.2mp,dens,6,3292.351605mp,ex,7,3.45e10mp,nuxy,7,0.2mp,dens,7,3284.137255mp,ex,8,3.45e10mp,nuxy,8,0.2mp,dens,8,3271.802136mp,ex,9,3.45e10mp,nuxy,9,0.2mp,dens,9,3260.41903mp,ex,10,3.45e10mp,nuxy,10,0.2mp,dens,10,3248.193657mp,ex,11,3.45e10mp,nuxy,11,0.2mp,dens,11,3235.117644mp,ex,12,3.45e10mp,nuxy,12,0.2mp,dens,12,3221.585664mp,ex,13,3.45e10mp,nuxy,13,0.2mp,dens,13,3208.826871mp,ex,14,3.45e10mp,nuxy,14,0.2mp,dens,14,3194.279207mp,ex,15,3.45e10mp,nuxy,15,0.2mp,dens,15,3179.924673mp,ex,16,3.45e10mp,nuxy,16,0.2mp,dens,16,3166.445716mp,ex,17,3.45e10mp,nuxy,17,0.2mp,dens,17,3152.555731mp,ex,18,3.45e10mp,nuxy,18,0.2mp,dens,18,3138.312105mp,ex,19,3.45e10mp,nuxy,19,0.2mp,dens,19,3124.795334mp,ex,20,3.45e10mp,nuxy,20,0.2mp,dens,20,3110.7135mp,ex,21,3.45e10mp,nuxy,21,0.2mp,dens,21,3097.080875mp,ex,22,3.45e10mp,nuxy,22,0.2mp,dens,22,3083.186268mp,ex,23,3.45e10mp,nuxy,23,0.2mp,dens,23,3068.968824mp,ex,24,3.45e10mp,nuxy,24,0.2mp,dens,24,3055.612436mp,ex,25,3.45e10mp,nuxy,25,0.2mp,dens,25,3045.857147mp,ex,26,3.45e10mp,nuxy,26,0.2mp,dens,26,3035.174287mp,ex,27,3.45e10mp,nuxy,27,0.2mp,dens,27,3026.696551mp,ex,28,3.45e10mp,nuxy,28,0.2mp,dens,28,3015.795365mp,ex,29,3.45e10mp,nuxy,29,0.2mp,dens,29,3007.710181mp,ex,30,3.45e10mp,nuxy,30,0.2mp,dens,30,3000.513837mp,ex,31,3.45e10mp,nuxy,31,0.2mp,dens,31,2978.611375 mp,ex,32,3.45e10mp,nuxy,32,0.2mp,dens,32,2958.618861mp,ex,33,3.45e10mp,nuxy,33,0.2mp,dens,33,2937.888072mp,ex,34,3.45e10mp,nuxy,34,0.2mp,dens,34,2919.475751mp,ex,35,3.45e10mp,nuxy,35,0.2mp,dens,35,2903.359983 !!6700处mp,ex,36,3.45e10mp,nuxy,36,0.2mp,dens,36,3302.153125 !!可以不用mp,ex,37,3.45e10mp,nuxy,37,0.2mp,dens,37,3302.153125mp,ex,38,3.45e10 !!-700处梁mp,nuxy,38,0.2mp,dens,38,3180.578901mp,ex,39,3.45e10 !!边支点横隔板mp,nuxy,39,0.2mp,dens,39,2868.674818mp,ex,40,3.45e10 !!合拢段横隔板mp,nuxy,40,0.2mp,dens,40,2868.674818mp,ex,41,3.45e10 !!中支点横隔板mp,nuxy,41,0.2mp,dens,41,2757.470588mp,ex,1,3.25e10 !!!!定义墩材料属性mp,nuxy,1,0.2mp,dens,1,2650!!定义实常数编号，面积，IYY，IZZ，宽，高，，RMORE，，抗扭惯距!!主梁截面r,1,11.851,28.52,215.151,16.95,4,, !!直线段rmore,,56.8972r,2,11.862,28.57,215.458,16.95,4.001,, !!截面100处rmore,,56.905r,3,11.871,28.91,215.622,16.95,4.009,,rmore,,57.345r,4,11.903,29.03,215.78,16.95,4.025,,rmore,,57.938r,5,11.932,29.626,215.88,16.95,4.049,,rmore,,58.502r,6,12.043,30.402,217.093,16.95,4.081,,rmore,,59.768r,7,12.199,31.407,218.899,16.95,4.121,,rmore,,61.445r,8,12.441,32.716,221.806,16.95,4.169,,rmore,,63.652r,9,12.673,34.21,224.418,16.95,4.226,,rmore,,66.058r,10,12.932,35.94,227.279,16.95,4.290,,rmore,,68.801r,11,13.221,37.941,230.512,16.95,4.362,,rmore,,71.948r,12,13.534,40.233,233.872,16.95,4.442,,rmore,,75.432r,13,13.843,42.747,237.133,16.95,4.530,,rmore,,79.173r,14,14.213,45.688,241.054,16.95,4.627,,rmore,,83.484r,15,14.598,48.964,245.004,16.95,4.731,,rmore,,88.155r,16,14.979,52.547,248.817,16.95,4.843,,rmore,,93.111r,17,15.393,56.562,252.927,16.95,4.963,,rmore,,98.528r,18,15.842,61.069,257.347,16.95,5.092,,rmore,,104.47r,19,16.293,65.992,261.673,16.95,5.228,,rmore,,110.742r,20,16.791,71.484,266.426,16.95,5.372,,rmore,,117.648r,21,17.303,77.568,271.225,16.95,5.525,, rmore,,124.861r,22,17.858,84.282,276.363,16.95,5.685,, rmore,,132.631r,23,18.464,91.768,281.923,16.95,5.854,, rmore,,141.068r,24,19.072,99.936,287.364,16.95,6.03,, rmore,,149.875r,25,19.524,108.482,291.159,16.95,6.214,, rmore,,158.46r,26,20.084,118.046,295.597,16.95,6.407,, rmore,,167.857r,27,20.536,128.202,298.9,16.95,6.607,, rmore,,177.137r,28,21.148,139.669,303.86,16.95,6.816,, rmore,,187.495r,29,21.626,151.737,307.186,16.95,7.032,, rmore,,197.507r,30,22.07,164.734,310.057,16.95,7.257,, rmore,,207.681r,31,23.541,182.225,323.433,16.95,7.490,, rmore,,223.527r,32,25.066,201.932,336.467,16.95,7.730,, rmore,,241.046r,33,26.871,223.700,351.448,16.95,7.979,, rmore,,257.765r,34,28.707,247.928,365.663,16.95,8.235,, rmore,,275.922r,35,30.533,274.36,378.934,16.95,8.500,, rmore,,294.461r,46,14.58,35.599,238.702,16.95,4,, rmore,,238.702!!横隔板截面!!边支点r,36,35.376,54.393,339.84,16.95,4.00,, rmore,,128.747!!合拢段r,37,35.376,53.405,339.84,16.95,4.00,, rmore,,127.015!!中支点r,38,71.981,454.871,543.871,16.95,8.5,, rmore,,684.565!!主墩截面!!截面1（实心）r,39,39.932,82.876,212.138,8,5,, rmore,,203.709!!截面（1/2空心）r,40,23.863,72.571,168.024,8,5,, rmore,,168.491!!截面3(空心）r,41,14.236,52.366,114.529,8,5,, rmore,,116.391!!边墩截面!!截面1（实心1）r,42,36.113,66.217,177.703,7.7,4.7,, rmore,,165.119!!截面2（空心1/2）r,43,20.753,35.743,99.711,7,4,, rmore,,88.486!!截面3（空心）r,44,9.103,16.082,45.911,6.5,3.5,, rmore,,39.39!!截面4（实心2）r,45,22.673,23.005,79.361,6.5,3.5,, rmore,,61.527!!截面5（1/2）空心2r,49,16,21.958,68.708,6.5,3.5,, rmore,,55.682!!建立第一跨梁节点xl=atan(1.029/80) n,1,0,n,2,2.25,2.25*xln,3,6,6*xln,4,9,9*xln,5,11,11*xln,6,13,13*xl-0.0045 n,7,15,15*xl-0.0135 n,8,17,17*xl-0.023 n,9,19,19*xl-0.038 n,10,21,21*xl-0.062 n,11,23,23*xl-0.092 n,12,25,25*xl-0.1275 n,13,27,27*xl-0.167 n,14,29,29*xl-0.2115 n,15,31,31*xl-0.261 n,16,33,33*xl-0.314 n,17,35,35*xl-0.3715 n,18,37,37*xl-0.435 n,19,39,39*xl-0.502 n,20,41,41*xl-0.5725 n,21,43,43*xl-0.648 n,22,45,45*xl-0.728 n,23,47,47*xl-0.8115 n,24,49,49*xl-0.8995 n,25,51,51*xl-0.9925 n,26,53,53*xl-1.0895 n,27,55,55*xl-1.1905 n,28,57,57*xl-1.325 n,29,59,59*xl-1.433 n,30,61,61*xl-1.5155 n,31,63,63*xl-1.6325 n,32,65,65*xl-1.7535 n,33,67,67*xl-1.879 n,34,69,69*xl-2.007 n,35,71,71*xl-2.1325 n,36,73,73*xl-2.2635 n,37,75,75*xl-2.399 n,38,77,77*xl-2.533 n,39,79,79*xl-2.631xl2=atan(0.7/140)local,11,0,80,80*xl-2.631,0,xl2nsym,x,35,5,39,1 !!复制粱结点从5到39结点编号增加35 csys,0n,75,149.6,0.15852local,12,0,150,0.15852,0,xl2nsym,x,36,40,75,1local,13,0,220,-1.6027,0nsym,x,111,1,111,1csys,0n,500,80,80*xl-2.631n,501,220,-1.6027n,502,360,-1.6148csys,0 !!返回普通坐标!!建立墩结点!!PM112边墩截面n,300,0,-2.165 !!(实心）2米长n,301,0,-4.365n,302,0,-9.365 !!（1/2空心）5M长n,303,0,-13.365n,304,0,-17.365n,305,0,-21.365n,306,0,-24.93n,307,0,-28.93n,308,0,-30.63!!PM113主敦截面n,309,80,0.5-6.665-0.8n,310,80,0.5-3.5-6.665 !!(实心)3.5mn,311,80,0.5-8.5-6.665 !!(1/2空心）5mn,312,80,0.5-12.5-6.665 !!(空心）4Mn,313,80,0.5-16.5-6.665 !!n,314,80,0.5-18.594-6.665n,315,80,0.5-22.594-6.665n,316,80,0.5-24.594-6.665+0.4!(实心）2m!!pm114主墩截面n,350,220,0.5-6.665-0.8 !n,351,220,0.5-3.5-6.665 !!(1/2空心)5m n,352,220,0.5-8.5-6.665 !!(空心）5mn,353,220,0.5-12.5-6.665 !!(空心）6.064M n,354,220,0.5-16.5-6.665 !!(1/2空心）4M n,355,220,0.5-19.294-6.665 !!(实心）2M n,356,220,0.5-23.294-6.665n,357,220,0.5-25.294-6.665+0.4local,13,0,220,-1.6027,0nsym,x,17,300,316,1csys,0!!生成边跨1mat,39real,36e,1,2mat,38real,46e,2,3mat,2real,1e,3,4e,4,5mat,2real,2e,5,6*do,i,3,35,1mat,ireal,ie,3+i,4+i*enddo*do,i,3,35,1mat,ireal,ie,38+i,39+i*enddo*do,i,3,35,1 mat,ireal,ie,74+i,75+i *enddo*do,i,3,35,1 mat,ireal,ie,114+i,115+i *enddo*do,i,3,35,1 mat,ireal,ie,149+i,150+i *enddo*do,i,3,35,1 mat,ireal,ie,185+i,186+i *enddomat,41 real,38e,39,500mat,41 real,38e,500,74mat,2real,2e,41,40mat,2real,1e,40,75mat,40 real,37 e,75,111mat,2 real,1 e,111,76mat,2 real,2 e,76,77mat,41 real,38 e,110,501mat,41 real,38 e,501,221mat,2 real,2 e,188,187mat,2 real,1 e,187,222mat,41 real,37 e,222,186mat,2 real,1 e,186,151mat,2 real,2 e,151,152mat,41 real,38 e,185,502mat,41 real,38 e,502,150mat,39 real,36 e,112,113mat,38 real,46 e,113,114mat,2 real,1e,114,115 e,115,116mat,2 real,2e,116,117!!建立桥墩!!边墩1 mat,1 real,42 e,300,301mat,1 real,43 e,301,302mat,1 real,44 e,302,303mat,1 real,44e,303,304mat,1 real,44 e,304,305mat,1 real,44 e,305,306mat,1 real,49 e,306,307mat,1 real,45 e,307,308!!主墩1 mat,1 real,39 e,309,310mat,1 real,40 e,310,311mat,1 real,41 e,311,312mat,1 real,41 e,312,313mat,1 real,41 e,313,314mat,1 real,40 e,314,315mat,1 real,39 e,315,316!!主墩2 mat,1 real,39 e,350,351mat,1 real,40 e,351,352mat,1 real,41 e,352,353mat,1 real,41 e,353,354mat,1 real,41 e,354,355mat,1 real,40 e,355,356mat,1 real,39 e,356,357!!主墩3 mat,1 real,39 e,326,327mat,1 real,40 e,327,328mat,1 real,41 e,328,329mat,1 real,41 e,329,330mat,1 real,41 e,330,331mat,1 real,40 e,331,332mat,1 real,39 e,332,333!!边墩2 mat,1 real,42 e,317,318mat,1 real,43 e,318,319mat,1 real,44 e,319,320mat,1 real,44 e,320,321mat,1 real,44 e,321,322mat,1real,44e,322,323mat,1real,49e,323,324mat,1real,45e,324,325!!添加质量惯距!!两米段质量惯距r,100,,,,1293348.4 !!100米处和直线段的质量惯距r,101,,,,1296019.6r,102,,,,1297493r,103,,,,1301181.8r,104,,,,1311723r,105,,,,1326622r,106,,,,1348967r,107,,,,1370728r,108,,,,1395061r,109,,,,1422801r,110,,,,1452757r,111,,,,1483258r,112,,,,1519733r,113,,,,1558030r,114,,,,1597229r,115,,,,1640292r,116,,,,1687605r,117,,,,1736625r,118,,,,1790923r,119,,,,1848602r,120,,,,1911418r,121,,,,1980562r,122,,,,2052690r,123,,,,2118097r,124,,,,2192307r,125,,,,2263641r,126,,,,2350704r,127,,,,2432292r,128,,,,2516392r,129,,,,2679987r,130,,,,2853514r,131,,,,3048284r,132,,,,3252032r,133,,,,3462506!!边跨直线3米段r,148,,,,1293348.4*3/2!!-700处3.75米段的质量惯距r,134,,,,2725866.2!!边支点横隔r,135,,,,2246142.5!!合龙段横隔r,136,,,,833679/2!!中支点横隔r,137,,,,5292760!!添加边墩质量惯距从上到下!!变截面实心2米段r,138,,,,1163498!!1/2空心5米段r,139,,,,1794766!!空心段4米r,140,,,,657211!!1/2空心4米段r,141,,,,1085079!!实心2米段r,142,,,,1064567!!添加主墩质量惯性距从上到下!!实心3米段r,143,,,,2345361!!1/2空心5米段r,144,,,,2869095!!空心4米段r,145,,,,1769087!!1/2空心4米段r,146,,,,2869095*4/5!!实心2米段r,147,,,,2345361*2/3!!在结点处加入质量惯性距!!边跨合龙段!左边边跨type,2real,135e,2type,2real,134e,3type,2real,148e,4type,2real,100e,5!!右边边跨type,2real,135e,113type,2real,134e,114type,2real,148e,115type,2real,100e,116!!对第一个主跨的循环*do,i,6,39,1type,2real,94+ie,i*enddo!!对第二个主跨的循环*do,i,41,74,1type,2real,59+ie,i*enddo!!对第3个循环*do,i,77,110,1type,2real,33+ie,i*enddo!!对第4个循环*do,i,188,221,1 type,2real,(-88)+ie,i*enddo!!对第5个循环*do,i,152,185,1 type,2real,i-52e,i*enddo!!对第6个循环*do,i,117,150,1 type,2real,(-17)+ie,i*enddo!!添加质量惯距到桥墩!!左边墩type,2real,138e,301type,2real,139e,302type,2real,140e,303type,2real,140e,304type,2real,140e,305type,2real,140e,306type,2 real,141 e,307type,2 real,142 e,308!!右边墩type,2 real,138 e,318type,2 real,139 e,319type,2 real,140 e,320type,2 real,140 e,321type,2 real,140 e,322type,2 real,140 e,323type,2 real,141 e,324type,2real,142e,325!!左边主墩type,2real,143e,310type,2real,144e,311type,2real,145e,312type,2real,145e,313type,2r,400,,,,1769087*3/4 e,314type,2real,146e,315type,2real,147e,316!!中间主墩type,2real,143e,351type,2real,144e,352type,2real,145e,353type,2real,145e,354type,2r,400,,,,1769087*3/4 e,355type,2real,146e,356type,2real,147e,357!!右边主墩type,2real,143e,327type,2real,144e,328type,2real,145e,329type,2real,145e,330type,2r,400,,,,1769087*3/4e,331type,2real,146e,332type,2real,147e,333!!添加横隔板处质量惯性距!!中支点type,2real,137e,500type,2real,137e,501type,2real,137e,502!!合龙段1type,2real,136e,75type,2real,136e,111type,2r,411,,,,1293348.4*3/10 e,40type,2r,412,,,,1293348.4*3/10 e,76!!合龙段2type,2real,136e,222type,2real,136e,186type,2r,411,,,,1293348.4*3/10 e,187type,2r,412,,,,1293348.4*3/10 e,151!!偶合结点cp,1,uy,1,300cp,2,uz,1,300cp,3,rotx,1,300cp,4,uy,309,500cp,5,uz,309,500cp,6,rotx,309,500cp,7,ux,350,501cp,8,uy,350,501cp,9,uz,350,501cp,10,rotx,350,501 cp,11,roty,350,501cp,12,uy,326,502cp,13,uz,326,502cp,14,rotx,326,502cp,15,uy,317,112cp,16,uz,317,112cp,17,rotx,317,112!!!!加约束d,308,all,all,0d,316,all,all,0d,357,all,all,0d,333,all,all,0d,325,all,all,0allselfinish!!!!求解/soluantype,2 !模态分析acel,,9.8 !加载重力modopt,subsp,20 !提取前30阶模态mxpand,20,,,0LUMPM,1solvefinish/post1set,list。
ANSYS空间梁格法分析连续斜梁桥命令流
/prep7!DEFINE THE ELEMENTARY PARAMETERS*DIM,L,ARRAY,10*DIM,DISTC,ARRAY,10N=4NBOX=2L(1)=20L(2)=30L(3)=30L(4)=20NUML=NBOX+3NUMC=(L(1)+L(2)+L(3)+L(4))+1CITA=25/180*3.1415926DISTC(1)=2.9584DISTC(2)=4.0802DISTC(3)=0.00DISTC(4)=4.0802DISTC(5)=2.9584!DEFINE THE NODES OF BRIDGEN, 1, 0, 0, 0.00N, NUMC, NUMC-1, 0, 0.00FILL,1,NUMCN, 200+1,(DISTC(1)+DISTC(2))*SIN(CITA)/COS(CITA), 0.0,-(DISTC(1)+DISTC(2))N,200+NUMC,NUMC-1+(DISTC(1)+DISTC(2))* SIN(CITA)/COS(CITA),0.0,-(DISTC(1)+DISTC(2)) FILL,201,200+NUMCN, 400+1, DISTC(2)*SIN(CITA)/COS(CITA), 0.0, -DISTC(2)N,400+NUMC,NUMC-1+DISTC(2)*SIN(CITA)/ COS(CITA), 0.0,-DISTC(2)FILL,401,400+NUMCN, 600+1, -DISTC(2)*SIN(CITA)/COS(CITA), 0.0, DISTC(2)N,600+NUMC,NUMC-1-DISTC(4)*SIN(CITA)/C OS(CITA),0.0, DISTC(4)FILL,601,600+NUMCN, 800+1,-(DISTC(4)+DISTC(5))*SIN(CITA)/COS(CITA ),0.0, DISTC(4)+DISTC(5) N,800+NUMC,NUMC-1-(DISTC(4)+DISTC(5))* SIN(CITA)/COS(CITA),0.0, DISTC(4)+DISTC(5) FILL,801,800+NUMC! 定义纵梁单元材料、几何参数ET,1,BEAM4MP,EX , 1, 3.45E+10MP,NUXY, 1, 0.2000MP,DENS, 1, 2600MP,ALPX, 1, 1.00E-05R,1,0.3133,0.0871,0.0959,1.0,0.25,,RMORE ,,0.00265,,,,,ET,2,BEAM4MP,EX , 2, 3.45E+10MP,NUXY, 2, 0.2000MP,DENS, 2, 2600MP,ALPX, 2, 1.00E-05R,2,1.6048,0.5486,1.5400,2.0,1.0,,RMORE ,,0.036,,,,,ET,3,BEAM4MP,EX , 3, 3.45E+10MP,NUXY, 3, 0.2000MP,DENS, 3, 2600MP,ALPX, 3, 1.00E-05R,3,2.3625,0.9885,2.2600,2.0,1.0,,RMORE ,,0.0956,,,,! 定义横梁单元材料、几何参数ET,4,BEAM4MP,EX , 4, 3.45E+10MP,NUXY, 4, 0.2000MP,DENS, 4, 0MP,ALPX, 4, 1.00E-05R,4,0.500,0.4582,0.0416,1.0,1.0,,RMORE ,,0.00879,,,,,ET,5,BEAM4MP,EX , 5, 3.45E+10MP,NUXY, 5, 0.2000MP,DENS, 5, 0MP,ALPX, 5, 1.00E-05R,5,1.0250,0.2770,0.0495,1.0,1.0,,RMORE ,,0.0609,,,,,! 生成纵向梁格单元TYPE,1MAT,1REAL,1*DO,I,1,NUMC-1,1 E,200+I,200+I+1 *ENDDO*DO,I,1,NUMC-1,1 E,800+I,800+I+1 *ENDDOTYPE,2MAT,2REAL,2*DO,I,1,NUMC-1,1 E,400+I,400+I+1 *ENDDO*DO,I,1,NUMC-1,1 E,600+I,600+I+1 *ENDDOTYPE,3MAT,3REAL,3*DO,I,1,NUMC-1,1 E,I,I+1*ENDDO! 生成横向梁格单元TYPE,5MAT,5REAL,5E,200+1,400+1 E,400+1,1E,1,600+1E,600+1,800+1 TYPE,4MAT,4REAL,4*DO,I,2,L(1),1E,200+I,400+IE,400+I,IE,I,600+IE,600+I,800+I *ENDDOTYPE,5MAT,5REAL,5E,200+L(1)+1,400+L(1)+1E,400+L(1)+1,L(1)+1E,L(1)+1,600+L(1)+1E,600+L(1)+1,800+L(1)+1TYPE,4MAT,4REAL,4*DO,I,L(1)+2,L(1)+L(2),1E,200+I,400+IE,400+I,IE,I,600+IE,600+I,800+I*ENDDOTYPE,5MAT,5REAL,5E,200+L(1)+L(2)+1,400+L(1)+L(2)+1E,400+L(1)+L(2)+1,L(1)+L(2)+1E,L(1)+L(2)+1,600+L(1)+L(2)+1E,600+L(1)+L(2)+1,800+L(1)+L(2)+1TYPE,4MAT,4REAL,4*DO,I,L(1)+L(2)+2,L(1)+L(2)+L(3),1E,200+I,400+IE,400+I,IE,I,600+IE,600+I,800+I*ENDDOTYPE,5MAT,5REAL,5E,200+L(1)+L(2)+L(3)+1,400+L(1)+L(2)+L(3)+1 E,400+L(1)+L(2)+L(3)+1,L(1)+L(2)+L(3)+1E,L(1)+L(2)+L(3)+1,600+L(1)+L(2)+L(3)+1E,600+L(1)+L(2)+L(3)+1,800+L(1)+L(2)+L(3)+1 TYPE,4MAT,4REAL,4*DO,I,L(1)+L(2)+L(3)+2,L(1)+L(2)+L(3)+L(4),1 E,200+I,400+IE,400+I,IE,I,600+IE,600+I,800+I*ENDDOTYPE,5MAT,5REAL,5E,200+L(1)+L(2)+L(3)+L(4)+1,400+L(1)+L(2)+L(3) +L(4)+1E,400+L(1)+L(2)+L(3)+L(4)+1,L(1)+L(2)+L(3)+L(4) +1E,L(1)+L(2)+L(3)+L(4)+1,600+L(1)+L(2)+L(3)+L(4) +1E,600+L(1)+L(2)+L(3)+L(4)+1,800+L(1)+L(2)+L(3) +L(4)+1/SOLACEL,,9.81,0!******************* ********************* ! 定义边界条件!**************************************** D,400+1,UX,0,,,,UY,UZD,600+1,UX,0,,,,UY,UZD,21,UY,0D,421,UY,0D,621,UY,0D,51,UY,0D,651,UY,0D,451,UY,0D,81,UY,0D,681,UY,0D,481,UY,0D,400+101,UX,0,,,,UY,UZD,600+101,UX,0,,,,UY,UZANTYPE,0/STATUS,SOLUSOLVE/POST1/REP,FASTA VPRIN,0, ,ETABLE, ,SMISC, 2A VPRIN,0, ,ETABLE, ,SMISC,8 A VPRIN,0, , ETABLE, ,SMISC, 4 A VPRIN,0, , ETABLE, ,SMISC, 5 A VPRIN,0, , ETABLE, ,SMISC, 6 A VPRIN,0, , ETABLE, ,SMISC, 10 A VPRIN,0, , ETABLE, ,SMISC, 11 A VPRIN,0, , ETABLE, ,SMISC, 12 FINISH。
ansys分析连续梁桥
fini/clear/title,60+100+60 three span continous bridge/prep7!上顶板的关键点位置k,1,0,0k,2,-1.5k,3,-2.5,-0.15k,4,-2.95,-0.15k,5,-3.4,-0.15k,6,-3.8,-0.075k,7,-4.2,0k,8,-4.6,0.075k,9,-5,0.15k,10,-5.4,0.225k,11,-6.1,0.225!下底板的关键点位置k,12,-2.95,-4.5k,13,,-4.5!上顶板关键点15m的位置kgen,4,1,13,1,0,0,5,100!上顶板关键点58m的位置kgen,10,301,311,1,0,0,43/9,100!上顶板关键点60m的位置kgen,2,1,11,1,0,0,60,1300!上顶板关键点62m的位置kgen,2,1,11,1,0,0,62,1400!上顶板关键点62+43m的位置kgen,10,1401,1411,1,0,0,43/9,100!上顶板关键点110m的位置kgen,2,2301,2311,1,0,0,5,100c1=0.0016225 !边、中跨的二次抛物线系数*dim,x1,array,9 !定义边跨的九分点位置*do,i,1,9,1x1(i)=i*43/9*enddo*dim,x2,array,9 !定义中跨的九分点位置*do,i,1,9,1x2(i)=i*43/9*enddo!下底板边跨（58m）九分点关键点的y坐标位置*dim,yb,array,9*do,i,1,9,1yb(i)=c1*x1(i)**2*enddo!下底板中跨（105m）九分点关键点的y坐标位置*dim,ym,array,9*do,i,1,9,1ym(i)=-c1*(x2(i)-43)**2+3*enddo!生成下底边跨（60m）九分点的关键点*do,i,1,9,1kgen,2,312,313,1,0,-yb(i),x1(i),i*100*enddo!下底板60m，62m处关键点kgen,3,1212,1213,1,0,0,2,100!生成下底中跨（62+43m）九分点的关键点*do,i,1,9,1kgen,2,1412,1413,1,0,ym(i),x2(i),i*100*enddo!生成下底中跨110m处关键点kgen,2,2312,2313,1,0,0,5,100!边墩处的横隔板位置kgen,2,1,13,1,0,0,-0.75,50!连成板面!上顶板,板厚0.5met,1,shell63mp,dens,1,2500mp,ex,1,3.5e10mp,prxy,1,0.1667r,1,0.5*do,i,0,23,1a,1+i*100,1+(i+1)*100,2+(i+1)*100,2+i*100*enddoa,1,51,52,2aatt,1,1,1!上顶板的加腋，板厚0.8mr,2,0.8*do,i,0,23,1a,2+i*100,2+(i+1)*100,3+(i+1)*100,3+i*100*enddoa,2,52,53,3aatt,1,2,1asel,none!肋板与翼缘交界处,板厚0.95mr,3,0.95*do,i,0,23,1a,3+i*100,3+(i+1)*100,4+(i+1)*100,4+i*100a,4+i*100,4+(i+1)*100,5+(i+1)*100,5+i*100*enddoa,3,53,54,4a,4,54,55,5aatt,1,3,1asel,none!翼缘厚度渐变，下面定义翼缘板的渐变厚度*dim,hd,array,7hd(7)=0.2hd(6)=0.2hd(5)=1/5*0.75+0.2hd(4)=2/5*0.75+0.2hd(3)=3/5*0.75+0.2hd(2)=4/5*0.75+0.2hd(1)=0.95*do,i,1,6,1r,30+i,hd(i),hd(i),hd(i+1),hd(i+1) !定义渐变的翼缘厚度*enddo*do,i,0,23,1*do,k,5,10,1a,k+i*100,k+(i+1)*100,k+1+(i+1)*100,k+1+i*100a,k,k+50,k+1+50,k+1aatt,1,k+26,1asel,none*enddo*enddo!肋板厚度0.9mr,4,0.9*do,i,0,23,1a,4+i*100,4+(i+1)*100,12+(i+1)*100,12+i*100*enddoa,4,54,62,12aatt,1,4,1asel,none!下底板厚度0.4～1.0m按二次抛物线过度c2=0.00035693 !边跨底板厚度变化的二次抛物线系数*dim,h1,array,16 !边跨底板厚度h1(1)=0.4h1(2)=0.4h1(3)=0.4*do,i,4,13,1h1(i)=0.4+c2*((i-4)*43/9)**2*enddoh1(14)=0.4+c2*43**2h1(15)=0.4+c2*43**2!生成边跨的底板单元*do,i,1,14,1r,40+i,h1(i),h1(i),h1(i+1),h1(i+1)a,12+(i-1)*100,12+i*100,13+i*100,13+(i-1)*100aatt,1,40+i,1asel,none*enddoa,12,62,63,13 !横隔板底面aatt,1,41,1*dim,h2,array,11 !中跨底板厚度h2(11)=0.4h2(10)=0.4h2(9)=0.4+c2*(1*43/9)**2h2(8)=0.4+c2*(2*43/9)**2h2(7)=0.4+c2*(3*43/9)**2h2(6)=0.4+c2*(4*43/9)**2h2(5)=0.4+c2*(5*43/9)**2h2(4)=0.4+c2*(6*43/9)**2h2(3)=0.4+c2*(7*43/9)**2h2(2)=0.4+c2*(8*43/9)**2h2(1)=0.4+c2*(9*43/9)**2!生成中跨的底板单元*do,i,1,10,1r,60+i,h2(i),h2(i),h2(i+1),h2(i+1)a,12+(i+13)*100,12+(i+14)*100,13+(i+14)*100,13+(i+13)*100 aatt,1,60+i,1asel,none*enddoallselesize,0.8mshape,0mshkey,1amesh,allnsel,allnsym,x,10000,allesym,x,10000,alllocal,11,0,,,110csys,11allselnsym,z,100000,allesym,z,100000,allallselnummrg,allnumcmp,allcsys,0。
ANSYS工程结构数值分析---命令流实例
!------------------------------------------------------!EX8.26 钢筋混凝土简支梁数值分析!分离式模型,关闭压碎,keyopt(1)=0,keyopt(7)=1 !力加载,位移收敛准则，误差1.5%,1/4模型分析!--------------------------------------------finish/clear/config,nres,2000/prep7!1.定义单元与材料性质--------------------et,1,solid65,,,,,,,1 !K1=0,k7=1et,2,link8mp,ex,1,13585mp,prxy,1,0.2fc=14.3ft=1.43tb,concr,1tbdata,,0.5,0.95,ft,-1tb,miso,1,,11tbpt,,0.0002,fc*0.19tbpt,,0.0004,fc*0.36tbpt,,0.0006,fc*0.51tbpt,,0.0008,fc*0.64tbpt,,0.001,fc*0.75tbpt,,0.0012,fc*0.84tbpt,,0.0014,fc*0.91tbpt,,0.0016,fc*0.96tbpt,,0.0018,fc*0.99tbpt,,0.002,fctbpt,,0.0033,fc*0.85mp,ex,2,2.0e5mp,prxy,2,0.3tb,biso,2tbdata,,300,0pi=acos(-1)r,1,0.25*pi*22*22r,2,0.25*pi*22*22/2r,3,0.25*pi*10*10r,4,0.25*pi*10*10/2!2.创建几何模型blc4,,,150/2,300,2000/2*do,i,1,9wpoff,,,100vsbw,all*enddo wpcsys,-1 wpoff,,,50 vsbw,all wpcsys,-1 wprota,,-90 wpoff,,,30 vsbw,all wpoff,,,240 vsbw,all wpcsys,-1 wpoff,30 wprota,,,90 vsbw,all wpcsys,-1!3.划分钢筋网格elemsiz=50 lsel,s,loc,x,30 lsel,r,loc,y,30 cm,zj,linelatt,2,1,2 lesize,all,elemsiz lsel,s,loc,x,75 lsel,r,loc,y,30 cm,zjb,line latt,2,2,2 lesize,all,elemsiz lsel,s,loc,x,30 lsel,r,loc,y,270 cm,jlj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,tan1,z lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 lsel,u,loc,z,50 cm,gj,linelatt,2,3,2 lesize,all,elemsiz lsel,s,loc,z,0 lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 cm,gjb,line latt,2,4,2lesize,all,elemsizlsel,allcmsel,s,zjcmsel,a,zjbcmsel,a,jljcmsel,a,gjcmsel,a,gjbcm,gj,linelmesh,alllsel,all!4.划分混凝土网格vatt,1,,1mshkey,1esize,elemsizvmesh,allallsel,all!5.施加荷载和约束lsel,s,loc,y,0lsel,r,loc,z,900dl,all,,uyasel,s,loc,z,0da,all,symmasel,s,loc,x,75da,all,symm!-----------------p0=180000q0=p0/150/100asel,s,loc,z,0,50asel,r,loc,y,300sfa,all,1,pres,q0allsel,all!6.求解控制设置/soluantype,0nsubst,80outres,all,allautos,onneqit,50cnvtol,u,,0.015solve!7.进入POST1查看结果/post1set,lastesel,s,type,,2etable,saxl,ls,1plls,saxl,saxlesel,s,type,,1/device,vector,onplcrackeall!8.进入时程后处理/post26nsol,2,205,u,yprod,3,1,,,,,,p0/1000prod,4,2,,,,,,-1xvar,4plvar,3!EX2.6D 60kg/m钢轨实体建模FINISH/CLEAR/PREP7!1.定义参数H1=30.5H2=48.5H3=48.5B1=45.75B2=29.25B3=16.5B4=36.5B5=10XD1=1/3XD2=1/9XD3=1/3XD4=1/20R0=400R1=2R2=4R3=40R4=20R5=25R7=13R8=80R9=300R10=2FAI=43H=H1+2*H2+H3!2.创建轨底主要组成部分的线K,1K,2,B1+B2K,3,B1+B2,H1-B1*XD1-B2*XD2K,4,B1,H1-B1*XD1K,5,0,H1*DO,I,1,4L,I,I+1*ENDDO!3.计算确定轨头下三个圆弧的位置CTA=A TAN(XD3)BTA=(90-FAI)*ACOS(-1)/180-CTAA=1/TAN(CTA)C1=B3/2+R0-R6*SIN(CTA)-(R5-R6)*COS(BTA)C2=R5*SIN(BTA)+R6*COS(CTA)-R6*SIN(BTA)-H2C=(C1-C2/TAN(CTA))/(R0-R5)C3=A*C-SQRT(A*A-C*C+1)C3=C3/(A*A+1)REFA=ASIN(C3)C1=R6*SIN(CTA)+(R5-R6)*COS(BTA)+(R0-R5)*COS(REFA) !4.创建轨腰和轨头结构线ROUX=B3/2+R0-C1ROUY=ROUX*TAN(CTA)ANGI=ASIN((H1+H2)/R0)K,6,B3/2+R0-R0*COS(ANGI)K,7,B3/2+R0-R0*COS(REFA),H1+H2+R0*SIN(REFA) LARC,6,7,2,R0DX1=R6*(COS(BTA)-SIN(CTA))DY1=R6*(COS(CTA)-SIN(BTA))K,8,ROUX-DX1,H1+2*H2+ROUY-DY1K,9,ROUX,H1+2*H2+ROUYLARC,7,8,2,R5LARC,8,9,2,R6K,10,B4,H1+2*H2+B4*XD3K,11,B4-(H3-B4*XD3)*XD4,HL,9,10L,10,11ANGI=ASIN(B5/R9)K,12,0,HK,13,B5,H-(R9-R9*COS(ANGI)) LARC,12,13,1,R9DX1=B4-B5+R8*SIN(ANGI)ANGI1=ASIN(DX1/R8)Y1=R8*COS(ANGI)-R8*COS(ANGI1) K,14,B4,KY(13)-Y1LARC,13,14,1,R8!5.倒角或弧线连接LFILLT,1,2,R1LFILLT,2,3,R2LFILLT,3,4,R3LPTN,4,5LFILLT,17,18,R4LFILLT,8,9,R10LPTN,9,11LFILLT,20,22,R7LDELE,15,16,1,1LDELE,19,21,2,1LSYMM,X,ALLNUMMRG,ALLNUMCMP,ALLAL,ALLASUMVOFFST,1,300!EX2.9B 斜向圆台建模finish/clear/prep7r0=3r1=2h=4cyl4,,,r0wpoff,r1,,hcyl4,,,r1*do,i,1,4askin,i,i+4*enddova,all!斜向延伸创建圆台vgen,1,1,,,,3*r0,,,,1wpcsyscyl4,,,r0vext,7,,,r1,,h,r1/r0,r1/r0!EX2.9B 已知函数方程时的建模!玫瑰线1==================================== !极坐标方程为ρ=asin(4Φ/3)finish/clear/prep7a=20csys,1*do,i,0,360*3fei=i*acos(-1)/180r0=a*sin(fei*4/3)k,,r0,i*enddo*do,i,1,360*3l,i,i+1*enddo!玫瑰线2=================================== !极坐标方程为ρ=asin(4Φ)finish/clear/prep7a=20csys,1*do,i,0,360fei=i*acos(-1)/180r0=a*sin(fei*4)k,,r0,i*enddo*do,i,1,360l,i,i+1*enddo!长辐圆内旋轮线============================ !参数方程x=(a-b)cost+λcos(a-b)/bt!参数方程y=(a-b)sint-λsin(a-b)/btfinish/clear/prep7a=50b=5lmda=15*afun,deg*do,i,0,360x1=(a-b)*cos(i)+lmda*cos((a-b)/b*i)y1=(a-b)*sin(i)-lmda*sin((a-b)/b*i)k,,x1,y1*enddo*do,i,1,360l,i,i+1*enddo!椭圆抛物面================================ !方程为z=x*x/a/a+y*y/b/b!使用蒙皮方法finish/clear/prep7a=20b=10n=20*do,i,1,niz=i/10lsel,nonewpoff,,,izcyl4,,,1adele,alllsscale,all,,,a*sqrt(iz),b*sqrt(iz),,,,1*enddoallsel*do,i,1,4*(n-1),4askin,i,i+4askin,i+1,i+5askin,i+2,i+6askin,i+3,i+7*enddo!EX2.9C 旋转图素建模!类花瓶建模--------------------finish/clear/prep7K,1K,2,100K,3,300,400K,4,150,550k,5,130,880K,6,300,1150k,7,400,1100k,8,500,1150bsplin,allk,1001,,1000arotat,all,,,,,,1,1001,360,5!类轮建模----------------------finish/clear/prep7!先创建类轮结构的一个断面，然后旋转之blc4,,,4,6cyl4,2,6,1.5asba,1,2k,100,-3k,101,-3,10vrotat,all,,,,,,100,101!EX3.1 自定义箱形截面finish/clear/prep7k,1k,2,2k,3,2,2.2k,4,3,2.3k,5,3,2.5k,6,0,2.5a,1,2,3,4,5,6blc4,,0.2,1.7,2asba,1,2wpoff,1.7wprota,,,90asbw,allwpoff,,,0.3asbw,allwpoff,,0.2wprota,,90asbw,allwpoff,,,-2asbw,allarsym,x,allwpcsysaglue,allet,1,plane82lesize,all,,,1amesh,allsecwrite,myboxfinish/clear/prep7et,1,beam189mp,ex,1,3.0e10mp,prxy,1,0.167sectype,1,beam,mesh secread,mybox,,,meshk,1k,2,10k,100,5,5l,1,2lesize,all,,,20latt,1,,1,,100,,1lmesh,all/eshape,1Eplot!EX3.2 自定义多种材料截面finish/clear/prep7Ro=1.5Ri=1.0csys,1cyl4,,,ricyl4,,,roaptn,allwprota,,90asbwa,allwprota,,,90asbw,allwpcsyset,1,plane82mymat1=4mymat2=7mp,ex,mymat1,1.0mp,ex,mymat2,2.0asel,s,loc,x,0,riaatt,mymat1,,1asel,s,loc,x,ri,roaatt,mymat2,,1allselesize,0.25mshape,0,2dmshkey,1amesh,allsecwrite,mycsolid,sectfinish/clear/prep7et,1,beam189mym1=4mym2=7mp,ex,mym1,3.0e10mp,prxy,mym1,0.167 mp,ex,mym2,2.1e11mp,prxy,mym2,0.3 sectype,1,beam,mesh secread,mycsolid,sect,,mesh k,1k,2,,,10l,1,2lesize,all,,,20latt,,,1,,,,1lmesh,all/eshape,1/pnum,mat,1eplot!EX3.3A 单个方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8K,1K,2,,,1000L,1,2K,100,,500,500LATT,1,,1,,100,,1LGEN,4,1,,,500LESIZE,ALL,,,10LMESH,ALL/ESHAPE,1EPLOT!EX3.3B 同时使用KB和KE方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0L,1,2K,100,,DLK,200,DXC,-DLK,300,2*DXC,DLK,301,2*DXC+DLK,400K,500,8*DXCLGEN,5,1,,,DXCLATT,1,,1,,100,,1LSEL,S,,,2LATT,1,,1,,200,,1LSEL,S,,,3LATT,1,,1,,300,301,1LSEL,S,,,4LATT,1,,1,,400,,1LSEL,S,,,5LATT,1,,1,,500,,1LSEL,ALLLESIZE,ALL,,,50LMESH,ALL/ESHAPE,1EPLOT!EX3.3C 双方位关键点翻转示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,I SECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,,,L0K,3,DXCK,4,DXC,,L0/2K,5,DXC,,L0L,1,2L,3,4L,4,5K,100,,DLk,101,,-DLK,200,DXC,DLk,201,2*DXCK,202,DXC,-DLLATT,1,,1,,100,101,1lesize,all,,,100LSEL,S,,,2LATT,1,,1,,200,201,1lesize,all,,,50LSEL,S,,,3LATT,1,,1,,201,202,1lesize,all,,,50LSEL,ALLLMESH,ALL/ESHAPE,1EPLOT!EX3.4A 两种单元形状和两种网格划分比较finish/clear/prep7et,1,plane82k,1k,2,8k,3,7,6k,4,1,6a,1,2,3,4esize,1mshape,0mshkey,0!mshape,0!mshkey,1!mshape,1!mshkey,0!mshape,1!mshkey,1amesh,all!EX3.4B 中间节点位置控制网格划分比较finish/clear/prep7cyl4,,,4,,8,60lesize,all,,,2mshape,0mshkey,1mshmid,0!mshmid,1!mshmid,2amesh,all!EX3.5 线上单元尺寸设置示例!下边密上边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10lsel,s,loc,x,0lesize,all,,,9,1/8lsel,s,loc,x,10lesize,all,,,9,8lsel,allmshape,0mshkey,1amesh,all!中间密外边稀finish/clear/prep7et,1,plane82blc4,,,10,10lsel,s,tan1,ylesize,all,,,10,-1/5lsel,s,tan1,xlesize,all,,,9,-1/8lsel,allmshape,0mshkey,1amesh,all!EX3.6 合并线和连接线以进行映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,3mshape,0mshkey,1lccat,1,2lcomb,4,5amesh,all!EX3.7 合并线或连接线的网格划分设置finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,10lesize,6,,,8lesize,1,,,4lesize,2,,,3lesize,4,,,5lesize,5,,,2mshape,0mshkey,1lcomb,1,2!lccat,1,2!如采用lcatt,1,2则采用7个划分数lcomb,4,5amesh,all!EX3.7B 简化映射网格划分finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1al,allesize,,6mshape,0amap,1,2,5,3,4!EX3.8A 过渡四边形映射网格finish/clear/prep7k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,8lesize,3,,,3lesize,4,,,7lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.8B 过渡四边形映射网格finish/clear/prep7et,1,plane42k,1k,2,10,-1k,3,8,6k,4,1,3a,1,2,3,4lesize,1,,,11lesize,3,,,3lesize,4,,,2lesize,2,,,2mshape,0,2dmshkey,1amesh,all!EX3.9A 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,3lesize,5,,,15allselmshape,0,3dvmesh,all!EX3.9B 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,all,,,4lesize,7,,,12mshape,0,3dmshkey,1vmesh,all!EX3.9C 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lsel,s,,,1,8,7lesize,all,,,6lsel,s,,,2,7,5lesize,all,,,3lsel,s,,,4,5lesize,all,,,7lsel,s,,,3,6,3lesize,all,,,2lsel,s,,,9,12,1lesize,all,,,5allselmshape,0,3dmshkey,1vmesh,all!EX3.9D 过渡六面体映射网格划分finish/clear/prep7et,1,95blc4,,,8,8,8lesize,2,,,2lesize,4,,,2lesize,5,,,2lesize,1,,,3lesize,3,,,3lesize,7,,,6lesize,6,,,7lesize,11,,,7lesize,9,,,3lesize,10,,,3lesize,12,,,3allselmshape,0,3dmshkey,1vmesh,all!EX3.10 带两孔的长方体的扫掠网格划分finish/clear/prep7a=10r=2et,1,mesh200,6et,2,solid45blc4,,,2*a,a,acyl4,a/2,a/2,r,,,,awprota,,90cyl4,1.5*a,a/2,r,,,,-avsbv,1,2vsbv,4,3wprota,,,90wpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpoff,,,a/2vsbw,allwpcsyswpoff,,a/2,a/2vsbw,allwprota,,90vsbw,allwpcsysesize,1amap,105,15,16,26,63amap,107,16,13,60,26amap,108,13,14,28,60amap,103,15,14,28,63lesize,94,,,4lesize,79,,,5vsel,s,loc,x,0,avsweep,allasel,s,loc,y,aasel,r,loc,x,a,2*alccat,2,45lccat,57,71lccat,68,78lccat,65,72mshape,0,2dmshkey,1amesh,allvsel,s,loc,x,a,2*avsweep,allallsel/view,1,1,2,3Eplot!EX3.11 拉伸类命令生成体单元网格finish/clear/prep7et,1,82et,2,95blc4,,,4,4blc4,6,,4,4blc4,12,,4,4blc4,18,,4,4esize,1amesh,allesize,,8vrotat,1,,,,,,1,4,90vext,2,,,,,10,0.5,0.5voffst,3,10vdrag,4,,,,,,35!EX3.12 圆的网格划分finish/prep7et,1,plane82r0=10cyl4,,,r0cyl4,3*r0,,,,r0 wprota,,90asbw,allwprota,,,90asbw,allwpoff,,,3*r0asbw,allwpcsys,-1asel,s,loc,x,-r0,r0lsla,slesize,all,,,8mshape,0,2dmshkey,1amesh,allallselasel,s,loc,x,2*r0,4*r0 lesize,all,,,8lsel,r,length,,r0lesize,all,,,8,0.1,1 amesh,allallsel!EX3.13 圆环的网格划分finish/clear/prep7et,1,plane82r0=10cyl4,,,r0/3,,r0,90cyl4,2*r0,,r0/10,,r0,90 asel,s,loc,x,-r0,r0lsla,slesize,all,,,8lsel,r,length,,r0*2/3 lesize,all,,,3,,1 mshape,0,2dmshkey,1amesh,allasel,s,loc,x,2*r0,4*r0lesize,5,,,12lesize,7,,,6lsel,s,length,,r0*9/10lesize,all,,,7amesh,all!EX3.14 圆柱面和圆柱体的网格划分!圆柱面finish/clear/prep7r0=10h0=50et,1,shell63cyl4,,,r0adele,1cm,l1cm,linek,50k,51,,,h0l,50,51adrag,l1cm,,,,,,5lsel,s,loc,z,0lesize,all,,,6lsel,s,length,,h0lesize,all,,,8mshape,0,2dmshkey,1amesh,all!圆柱体finish/clear/prep7r0=10h0=50et,1,solid95cyl4,,,r0,,,,h0wprota,,90vsbw,allwprota,,,90vsbw,allmshape,0,3dmshkey,1lsel,s,loc,z,0lsel,s,length,,h0lesize,all,,,8vmesh,all!EX3.15 圆锥的网格划分!扫掠网格划分finish/clear/prep7cone,10,,,15,,90et,1,200,7et,2,95lsel,s,,,5,6lesize,all,,,12,0.5lsel,alllesize,3,,,6mshape,0,2dmshkey,1amesh,3vsweep,1,3,4vsymm,x,allvsymm,y,allvglue,all!六面体映射网格划分finish/clear/prep7cone,10,,,15,,90et,2,95lesize,all,,,8mshape,0,3dmshkey,1vmesh,allvsymm,x,allvsymm,y,allvglue,all!EX3.16 正多边形面的通用网格划分finish/clear/prep7ns=11rpr4,ns,,,10kp0=100+nsk,kp0*if,mod(ns,2),eq,0,then*do,i,1,ns/2l,kp0,2*i-1*enddo*else*do,i,1,nsl,kp0,i*enddo*endiflsel,s,,,ns+1,2*nscm,l1cm,lineallselasbl,1,l1cm*if,mod(ns,2),eq,0,then cmsel,s,l1cmlesize,all,,,8lsel,invelesize,all,,,4*elselesize,all,,,8*endifallselmshape,0mshkey,1amesh,all!EX3.17 球体及球面网格划分finish/clear/prep7r0=10sphere,,r0,,90vsbw,allvdele,2,,,1esize,2et,1,solid95mshape,0,3dmshkey,1vsymm,x,all vsymm,y,all vsymm,z,allvglue,allfinish/clear/prep7r0=10sphere,,r0,,90 vsbw,allvdele,2,,,1vdele,allasel,s,loc,x,0asel,a,loc,y,0asel,a,loc,z,0adele,all,,,1allselesize,2et,1,shell63 mshape,0,2d mshkey,1amesh,allarsym,x,allarsym,y,allarsym,z,all nummrg,all!EX3.18 孔板网格划分finish/clear/prep7a0=300b0=800r0=15blc4,,,a0,b0cyl4,a0/4,b0/8,r0 agen,2,2,,,a0/2 agen,2,2,3,1,,b0/8 agen,2,2,5,1,,b0*5/8 asel,s,,,2,9,1cm,a2cm,areaallselwprota,,-90*do,i,1,5wpoff,,,b0/16asbw,all*enddowpoff,,,b0*5/16*do,i,1,5wpoff,,,b0/16asbw,all*enddowprota,,,90*do,i,1,3wpoff,,,a0/4asbw,all*enddowpcsys,-1numcmp,alllsel,s,radius,,r0lesize,all,,,8lsel,invelesize,all,,,4lsel,allet,1,82mshape,0,2dmshkey,1asel,u,loc,y,b0/16,b0*5/16 asel,u,loc,y,b0*11/16,b0*15/16 lsla,slsel,r,tan1,xlesize,all,50,,,,1amesh,allallselamap,21,30,31,54,62 amap,22,31,32,54,64 amap,32,29,30,62,76 amap,33,29,32,64,76amap,42,34,35,76,82 amap,43,35,36,76,84 amap,41,33,34,53,82 amap,44,33,36,53,84amap,19,22,23,49,64 amap,20,23,24,49,60 amap,9,21,22,64,74 amap,31,21,24,60,74amap,38,26,27,74,84 amap,39,27,28,74,80 amap,37,25,26,50,84 amap,40,25,28,50,80amap,17,14,15,43,59 amap,18,15,16,43,63 amap,7,13,14,59,73 amap,8,13,16,63,73amap,34,18,19,73,79 amap,30,17,18,79,44 amap,35,19,20,73,83 amap,36,17,20,44,83amap,15,6,7,39,63amap,16,7,8,39,58amap,4,6,5,63,71amap,5,5,8,58,71amap,27,10,11,71,83 amap,28,11,12,71,78 amap,26,9,10,83,40 amap,29,9,12,40,78!EX3.19 角支架的网格划分finish/clear/prep7blc4,,,150,50blc4,100,,50,-50cyl4,,25,25cyl4,125,-50,25aadd,allnumcmp,allcyl4,,25,10cyl4,125,-50,10 asel,s,,,2,3cm,a1cm,area asel,allasba,1,a1cm lcomb,1,6lfillt,1,2,20 asbl,4,6adele,1,,,1lfillt,3,4,20al,18,19,20 aadd,all numcmp,all!---------------------- wprota,,90 wpoff,,,-25 asbw,all wpoff,,,75 asbw,all wprota,,,90 asbw,all wpoff,,,125 asbw,all wpcsys,-1 wpoff,25 wprota,,,90 asbw,all kwpave,18 asbw,all kwpave,3 wprota,,90 asbw,all kwpave,21 wprota,,-45 asbw,8 wpcsys,-1!-------------------- et,1,plane82 mshape,0,2d mshkey,1 esize,6lesize,33,,,6 lesize,37,,,6 lesize,42,,,6amesh,3,5,2amap,6,9,10,4,23amap,7,9,12,1,23amesh,11,13,2amesh,1,9,8amesh,2,4,2amap,16,3,18,26,28amap,12,14,15,5,28amap,15,13,14,6,28lcomb,19,22lcomb,27,46amesh,10,14,4!========================= !EX3.20 大板小孔的网格划分A finish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1wprota,,90,45asbw,allet,1,plane82mshape,0,2dmshkey,1lsel,s,loc,x,0,a0/20lesize,all,,,8lsel,alllesize,7,,,12,0.1lesize,8,,,12,0.1lesize,9,,,12,10amesh,allcsys,0arsym,x,allarsym,y,allnummrg,all!============================= !EX3.20 大板小孔的网格划分Bfinish/clear/prep7a0=100blc4,,,a0,a0cyl4,,,a0/100asba,1,2csys,1k,50,a0/20k,51,a0/20,90l,50,51asbl,all,1et,1,plane82mshape,0,2dmshkey,1lesize,5,,,8lesize,1,,,10lesize,4,,,8,5lesize,6,,,8,5amesh,1lesize,7,,,20,0.1lesize,8,,,20,0.1amap,2,50,51,2,4csys,0arsym,x,allarsym,y,allnummrg,all!EX3.21 Z形面网格划分及局部细分finish/clear/prep7blc4,,,15,10blc4,10,6,14,12aadd,allwprota,,-90wpoff,,,6asbw,allwpoff,,,4wprota,,,90wpoff,,,10asbw,allwpoff,,,5asbw,allwpcsys,-1et,1,82esize,2mshape,0,2dmshkey,1amesh,allkrefine,9,10,1,1,,off !不光滑处理（节点位置不变）aclear,allamesh,allkrefine,9,10,1,1,,smooth !光滑处理，节点位置可改变，但既有单元不删除aclear,allamesh,allkrefine,9,10,1,1 !光滑处理，可删除既有单元重细分，节点位置改变!EX3.22 分布小孔环形面的网格划分finish/clear/prep7r1=15r2=25r3=3n=8cyl4,,,r1,,r2,180/ncyl4,0.5*(r1+r2),,r3asba,1,2kl,2larc,6,8,4,0.5*(r1+r2)asbl,all,4lesize,4,,,4lesize,10,,,4lesize,6,,,8lesize,3,,,6lesize,8,,,4lesize,5,,,8lesize,9,,,4lesize,7,,,4et,1,plane82mshape,0,2damap,2,6,7,4,8amap,1,5,6,8,1!lrefine,5,6,1,1,smootharsym,y,allcsys,1agen,n,all,,,,360/nnummrg,all!EX3.23 同心圆轴体的网格划分finish/clear/prep7cylind,50,0,100,150cylind,40,0,200,250cylind,30,0,100,250cylind,20,0,50,250cylind,10,0,0,250vptn,allwprota,,,90vsbw,allwprota,,90vsbw,allet,1,95csys,1lsel,s,length,,50lesize,all,,,5lsel,invelsel,u,radius,,10,50lesize,all,,,4n1=6n2=6n3=8n4=10n5=12lsel,r,loc,x,0,10lesize,all,,,n1,,1lsel,s,radius,,10lesize,all,,,n1lsel,s,radius,,20lsel,s,radius,,30lesize,all,,,n3lsel,s,radius,,40lesize,all,,,n4lsel,s,radius,,50lesize,all,,,n5allselmshape,0,3dmshkey,1vmesh,all!EX3.24A 圆柱与长方体组合模型的网格划分!采用六面体映射网格划分finish/clear/prep7a0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,allaccat,4,6et,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.24B 圆柱与长方体组合模型的网格划分!采用过渡六面体映射网格划分finish/cleara0=30h1=15h2=25r0=7blc4,,,a0/2,a0/2,h1cyl4,,,r0,,,90,h1+h2vptn,alllsel,s,radius,,r0lsel,a,length,,r0lesize,all,,,6lsel,s,loc,z,0lsel,a,loc,z,h1lesize,all,,,7asel,s,loc,x,a0/2asel,a,loc,y,a0/2accat,allallselet,1,95esize,3mshape,0,3dmshkey,1vmesh,allasel,s,accaadele,alllsel,s,lccaldele,allallselvsymm,x,allvsymm,y,allnummrg,all!EX3.25-任意圈数圆柱形螺旋弹簧的扫略网格划分FINISH/CLEAR/PREP7!1.定义弹簧参数（同上）D=4C=8DZ=C*DT=DZ/6N=4.7*IF,T,LT,D,THEN*ENDIFTKPD=90!2.创建全部螺旋线CSYS,1TDEG=N*360TDEG1=MOD(TDEG,TKPD) N0=(TDEG-TDEG1)/TKPD+1 *AFUN,DEG*DO,I,1,N0CTA=(I-1)*TKPDZ=T/360*CTAK,I,DZ/2,CTA,Z*ENDDO*IF,TDEG1,LT,1.0E-2,THEN *ELSEN0=N0+1CTA=CTA+TDEG1Z=T/360*CTAK,N0,DZ/2,CTA,Z*ENDIF*DO,I,1,N0-1L,I,I+1*ENDDOCM,L1,LINE!3.在螺旋线端部创建簧丝截面CSYS,0WPOFF,DZ/2 WPROTA,,90CYL4,,,D/2wprota,,90asbw,allwprota,,,90asbw,allcm,a1cm,areaVDRAG,a1cm,,,,,,L1et,1,mesh200,7et,2,solid95cmsel,s,a1cmlsla,slesize,all,,,6amesh,allallselesize,d/2vsweep,all!EX3.26 具坑缺陷圆柱的网格划分finish/clear/prep7l0=50r0=10r1=5r2=1h0=5et,1,solid45cyl4,,,r0,,,,l0wpave,0,r0,l0/2wprota,,90con4,,,r1,r2,h0vsbv,1,2wprota,,90vsbw,allwpoff,,,r1+3vsbw,allwpoff,,,-2*(r1+3)vsbw,allwprota,,,90vsbw,allwpoff,,,r1+1vsbw,allwpoff,,,-2*(r1+1)vsbw,allwpcsys,-1wprota,,90vsbw,allwpoff,,,-r1vsbw,allwpoff,,,2*r1vsbw,allwpcsys,-1numcmp,allksel,s,loc,z,l0/2ksel,r,loc,x,0ksel,u,loc,y,0*get,kp1,kp,0,num,min kp2=kpnext(kp1) ksel,s,loc,z,l0/2+r2*get,kp3,kp,0,num,min allselk,1000,,-r0,l0/2+5l,kp3,1000asel,none*get,l1,line,0,num,max arotat,l1,,,,,,kp1,kp2 cm,a1cm,areaallselvsba,all,a1cm!vsel,s,loc,z,17,33!vsel,r,loc,x,-6,6!aslv,s!aplotaccat,128,61accat,160,138 accat,224,202 accat,183,35accat,31,134accat,151,164 accat,215,231 accat,82,195accat,130,5accat,162,145 accat,229,209 accat,192,89accat,12,123accat,117,155 accat,109,219 accat,14,178mshape,0,3d mshkey,1esize,1vmesh,all!EX3.27长方体开圆柱槽的网格划分finish/clear/prep7blc5,,,20,10,30wpoff,,5,5cyl4,,,8,,,,10wpoff,,,10cyl4,,,4,,,,12vsel,s,,,2,3cm,v1cm,voluallselwpcsys,-1vsbv,1,v1cmvoffst,15,5voffst,16,3voffst,17,15vptn,allwprota,,,90vsbw,allwpcsys,-1wpoff,,,5vsbw,allwpoff,,,10vsbw,allwpoff,,,12vsbw,allet,1,mesh200,7esize,2mshape,0,2dmshkey,1asel,s,loc,z,30lccat,6,74lccat,8,68amesh,allasel,s,loc,z,0lccat,3,73lccat,1,67amesh,allallselet,2,solid95mshape,0,3dmshkey,1vsweep,all!ex4.1累加节点自由度约束值finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,1lmesh,alld,1,ux,1e-2,,,,uyd,1,rotzd,2,alldlistnsel,s,d,u,0,0.1dofsel,s,uydscale,2.5dlistdofsel,alldcum,addd,1,ux,2e-2dlistdcum,ignod,2,uy,1.0dlist!ex4.2 对线施加约束并转换finish/clear/prep7et,1,95blc4,,,10,10,10dl,7,,ux,0.1dl,5,,alldl,11,6,symmdl,10,6,asymdl,6,,symmdllistesize,2vmesh,alldtrandlist!ex4.3 施加集中荷载与节点坐标系finish/clear/prep7et,1,beam4k,1k,2,5k,3,10l,1,2l,2,3local,12,0,,,,90nrotat,alldk,1,allfk,2,fy,-1000esize,1lmesh,allnrotat,alllplotfk,3,fy,1000f,6,fx,-1000sbctraneplot!ex4.4a 3d单元sf加载示例finish/clear/prep7et,1,95blc4,,,10,10,20esize,,4vmesh,allasel,s,loc,y,10sf,all,pres,1000asel,s,loc,z,20nsla,s,1sf,all,pres,1000!ex4.4b 2d单元sf加载示例finish/clear/prep7et,1,82blc4,,,100,200blc4,30,60,40,80asba,1,2wprota,,-90wpoff,,,60asbw,allwpoff,,,80asbw,allwprota,,,90wpoff,,,30asbw,allwpoff,,,40asbw,allwpcsys,-1esize,5amesh,all/psf,pres,norm,2sf,all,pres,100sfdele,all,presnsel,s,loc,x,0sf,all,pres,100nsel,s,loc,x,15,20esln,s,1nsel,r,loc,x,15sf,all,pres,110nsel,s,loc,x,40,60nsel,r,loc,y,10,30esln,s,1sf,all,pres,100lsel,s,loc,x,100nsll,s,1esln,snsel,s,loc,x,95sf,all,pres,-100eplot!ex4.5 节点号及其荷载函数finish/clear/prep7et,1,45blc4,,,10,10,20esize,5vmesh,all*dim,mypres,,100*do,i,1,100mypres(i)=i*10.0*enddosffun,pres,mypres(1)nsel,s,loc,y,10sf,all,pres,10sflist*do,i,1,100mypres(i)=i*50.0*enddonsel,s,loc,z,20sf,all,pres,0allselsflist!ex4.6a 2d平面单元plane82 finish/clear/prep7et,1,82blc4,,,10,50esize,2amesh,allsfe,1,4,pres,,100,50sfe,6,4,pres,,200,40nsel,s,loc,x,10esln,ssfe,all,2,pres,,-100/psf,pres,norm,2,0,1eplot!ex4.6b 3d体单元solid95 finish/clear/prep7et,1,95blc4,,,10,10,30esize,5vmesh,all/psf,pres,norm,2,0,1sfe,22,3,pres,,100sfe,22,6,pres,,100,110,120,130 sfelist!ex4.6c 3d壳单元shell63 finish/clear/prep7et,1,63wprota,,90blc4,,,10,10esize,5amesh,all/psf,pres,norm,2,0,1sfe,1,1,pres,,100sfe,1,3,pres,,100sfe,1,6,pres,,100sfe,4,1,pres,,100,110,120,130!ex4.7 在梁单元上施加荷载finish/clear/prep7et,1,beam3k,1k,2,10l,1,2esize,,10lmesh,all/pnum,elem,1sfbeam,3,1,pres,50,100 sfbeam,5,1,pres,100 sfbeam,7,1,pres,50,100,,,0.2,0.1sfbeam,9,1,pres,100,,,,0.4,-1sfbeam,3,2,pres,50,100,,,0.2,0.1!ex4.8a 在线上施加面荷载finish/clear/prep7et,1,82blc4,,,10,30esize,5amesh,all/psf,pres,norm,2sfl,4,pres,10,60sfl,2,pres,60sftraneplot!ex4.8b 在线上施加面荷载finish/clear/prep7et,1,shell63wprota,,90blc4,,,10,30esize,5amesh,allsfl,3,pres,100sftraneplot!ex4.9 利用荷载梯度在直角坐标系下的施加方法finish/clear/prep7et,1,82blc4,,,10,60esize,2amesh,all/psf,pres,norm,2sfgrad,pres,,y,0,-5nsel,s,loc,x,0nsel,r,loc,y,0,40sf,all,pres,600sfgrad,pres,,y,30,-20nsel,s,loc,x,10sf,all,pres,0allseleplot!ex4.10利用荷载梯度在柱坐标系下的施加方法finish/clear/prep7csys,1k,1,10,-90k,2,10,90k,3,10,90,30l,1,2l,2,3adrag,1,,,,,,2ldele,2,,,1numcmp,allet,1,63esize,2amesh,allsfgrad,pres,1,y,-90,1sfa,all,2,pres,400sftran!可考察下述命令及其施加方法!①!sfgrad,pres,1,y,270,1sfa,all,2,pres,400sftran!②!local,12,1cscir,12,1sfgrad,pres,12,y,270,1sfa,all,2,pres,400sftran!ex4.11 ldread命令的示例finish/clear/filname,ldtest1/prep7et,1,plane82mp,ex,1,2.1e5mp,prxy,1,0.3blc4,,,10,40esize,2amesh,alld,all,allsfl,3,pres,100/solusolve/post1prrsolfinish/filename,ldtest2/solulsclear,allnsel,s,loc,y,0d,all,allallselldread,reac,,,,,ldtest1,rst solve/post1plnsol,s,y!ex4.12 初应力荷载finish/clear/filname,colu1/prep7et,1,plane82mp,ex,1,2e5mp,nuxy,1,0.3blc4,,,1,10esize,2amesh,allnsel,s,loc,y,0d,all,uyd,1,uxnsel,s,loc,y,10sf,all,pres,-10。
ansys结构分析入门-命令流
1 ANSYS概述1.1 ANSYS简介ANSYS是一种广泛的商业套装工程分析软件。
所谓工程分析软件，主要是在机械结构系统受到外力负载所出现的反应，例如应力、位移、温度等，根据该反应可知道机械结构系统受到外力负载后的状态，进而判断是否符合设计要求。
一般机械结构系统的几何结构相当复杂，受的负载也相当多，理论分析往往无法进行。
想要解答，必须先简化结构，采用数值模拟方法分析。
由于计算机行业的发展，相应的软件也应运而生，ANSYS软件在工程上应用相当广泛，在机械、电机、土木、电子及航空等领域的使用，都能达到某种程度的可信度，颇获各界好评。
使用该软件，能够降低设计成本，缩短设计时间。
到80年代初期，国际上较大型的面向工程的有限元通用软件主要有：ANSYS, NASTRAN, ASKA, ADINA, SAP等。
以ANSYS为代表的工程数值模拟软件，是一个多用途的有限元法分析软件，它从1971年的2.0版本与今天的5.7版本已有很大的不同，起初它仅提供结构线性分析和热分析，现在可用来求结构、流体、电力、电磁场及碰撞等问题的解答。
它包含了前置处理、解题程序以及后置处理，将有限元分析、计算机图形学和优化技术相结合，已成为现代工程学问题必不可少的有力工具。
1.2ANSYS软件主要功能ANSYS软件是融结构、热、流体、电磁、声学于一体的大型通用有限元软件，可广泛的用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、生物医学、水利、日用家电等一般工业及科学研究。
该软件提供了不断改进的功能清单，具体包括：结构高度非线性分析、电磁分析、计算流体力学分析、设计优化、接触分析、自适应网格划分及利用ANSYS参数设计语言扩展宏命令功能。
1.3ANSYS软件主要特点主要技术特点：∙唯一能实现多场及多场耦合分析的软件∙唯一实现前后处理、求解及多场分析统一数据库的一体化大型FEA软件∙唯一具有多物理场优化功能的FEA软件∙唯一具有中文界面的大型通用有限元软件∙强大的非线性分析功能∙多种求解器分别适用于不同的问题及不同的硬件配置∙支持异种、异构平台的网络浮动，在异种、异构平台上用户界面统一、数据文件全部兼容∙强大的并行计算功能支持分布式并行及共享内存式并行∙多种自动网格划分技术∙良好的用户开发环境支持的图形传递标准：∙SAT∙Parasolid∙STEP与CAD软件的接口∙Unigraphics∙Pro/ENGINEER∙I-Deas∙Catia∙CADDS∙SolidEdge∙SolidWorks1.4 运行环境（ANSYS5.7）Computer: Pentium-class systemMemory (RAM): 64 MB以上Hard Disk: 500MB以上自由空间.Operating System: Microsoft Windows 2000, Windows NT 4.0 (SP 5 or higher) or Windows 98 Graphics: A Windows 2000/NT 4.0 or 98 supported Graphics Card, capable of 1024x768 in High Color (16-bit).A 17 inch monitor (or larger) compatible with the above mentioned card is recommended.2ANSYS 的基本使用2.1 ANSYS环境简介ANSYS有两种模式：一种是交互模式（Interactive Mode），另一个是非交互模式（Batch Mode）。
曲线梁桥ANSYS计算命令流
!****************************************************************************** ***********************! case2:无偏载(以跨径布置30m+40m+30m,桥宽8.5为例）! 上海城市设计研究院L1+L2+L3预应力混凝土曲线连续梁桥结构分析! 两端为抗扭支座，中间支座为点铰支座! 每次要记得修改横隔梁的参数，即Mass21单元的实常数!****************************************************************************** ***********************FINI/CLE/prep7!DEFINE THE ELEMENTARY PARAMETERS*DIM,L,ARRAY,10*DIM,H,ARRAY,10*DIM,CITA,ARRAY,10!*****以下参数均可修改***************N=3 !跨数L(1)=30 !第一跨L(2)=40 !第二跨L(3)=30 !第三跨e1=1.25 !1#墩处内支座到中心线的间距e2=1.25 !1#墩处外支座到中心线的间距e3=0 !2#墩处的支座偏心距（正的表示外偏）e4=0 !3#墩处的支座偏心距e5=1.25 !4#墩处内支座到中心线的间距e6=1.25 !4#墩处外支座到中心线的间距R=10000 !曲线桥半径H0=1.0 !梁底到截面形心处的高度M=16146 !mass21单元质量J=27246.38 !mass21单元转动惯量!*************************************LL=0.0*DO,I,1,NLL=LL+L(I)CITA(I)=L(I)/R/3.1415925*180*ENDDOCITA0=LL/R/3.1415925*180LOCAL,11,1,0,0,0,0,270,0CSYS,11!DEFINE THE NODES OF BRIDGEN, 1, R, -179.9999, H0N, L(1)+1, R, -179.9999+CITA(1), H0N, L(1)+L(2)+1, R, -179.9999+CITA(1)+CITA(2), H0 N, L(1)+L(2)+L(3)+1, R, -179.9999+CITA0, H0FILL,1,L(1)+1FILL,L(1)+1,L(1)+L(2)+1FILL,L(1)+L(2)+1,L(1)+L(2)+L(3)+1!DEFINE THE REFERENCE POINT OF MAIN DECKN,2000,0,0,0!***************************************************************** !DEFINE THE NODES OF RIGID BEAM OF MAIN GIRDERN, 2001, R-e1, -179.9999, H0N, 2002, R+e2, -179.9999, H0N, 2003, R+e3, -179.9999+CITA(1), H0N, 2004, R+e4, -179.9999+CITA(1)+CITA(2), H0N, 2005, R-e5, -179.9999+CITA0, H0N, 2006, R+e6, -179.9999+CITA0, H0NROTA T,ALL!*****************************************************! 定义主梁单元材料、几何参数! 混凝土标号为C40!**************梁段截面几何特性*******************ET,1,BEAM4MP,EX , 1, 3.25E+10MP,NUXY, 1, 0.2000MP,DENS, 1, 3950.50MP,ALPX, 1, 1.00E-05R,1,4.39129,2.4428,20.7089,8.5,2.0,, !跨中梁断面RMORE ,,5.0388,,,,,!**************刚臂截面几何特性*******************! 增设抗扭或偏心支座时采用ET,3,BEAM4MP,EX , 3, 3.25E+10MP,NUXY, 3, 0.2000MP,DENS, 3, 0.0MP,ALPX, 3, 1.0E-08R,3,100,1000,1000,0.5,0.5,,RMORE ,,1000,,,,,!****************生成主梁单元**********************TYPE,1MA T,1REAL,1*DO,I,1,LL,1E,I,I+1,2000*ENDDO!*********************************************************** ! 生成主梁刚臂单元!*********************************************************** TYPE,3MA T,3REAL,3E,2001,1E,1,2002*IF,e3,GT,0.0,THENE,2003,L(1)+1*ELSENSEL,U, , ,2003*ENDIF*IF,e4,GT,0.0,THENE,2004,L(1)+L(2)+1*ELSENSEL,U, , ,2004*ENDIFE,2005,LL+1E,LL+1,2006!**************************************************! 支墩处的横隔梁采用Mass21单元来模拟!**************************************************ET,4,MASS21R,100, M ,M ,M ,J ,,,TYPE,4REAL,100E,1E,L(1)+1E,L(1)+L(2)+1E,LL+1!****施加约束条件**********D,2001,UX,,,,,UY,UZD,2002,UZ*IF,e3,GT,0.0,THEND,2003,,,,,,,UZ*ELSED,L(1)+1,,,,,,,UZ*ENDIF*IF,e4,GT,0.0,THEND,2004,,,,,,,UZ*ELSED,L(1)+L(2)+1,,,,,,,UZ*ENDIFD,2005,UX,,,,,,UZD,2006,UZ!********************************************************** ! 对桥梁桥面各节点逐点加载!********************************************************** /SOLUANTYPE,0 ! 静力分析ACEL,0,0,0K=1*DO,I,1,LL+1,1TIME,KFDELE,ALL,ALLF,I,FZ,-1000SOLVEK=K+1*ENDDO!******定义各片纵梁的竖向挠度、弯矩、扭矩、剪力的包络图******** *DIM,AMAX,ARRAY,LL+1*DIM,AMIN,ARRAY,LL+1*DIM,JZMAX,ARRAY,LL+1*DIM,JZMIN,ARRAY,LL+1*DIM,DISP1,ARRAY,LL+1,LL+1*DIM,MZ1,ARRAY,LL,LL+1*DIM,MX1,ARRAY,LL,LL+1*DIM,QY1,ARRAY,LL,LL+1*DIM,DISP1MAX,ARRAY,LL+1*DIM,DISP1MIN,ARRAY,LL+1*DIM,MZ1MAX,ARRAY,LL*DIM,MZ1MIN,ARRAY,LL*DIM,MX1MAX,ARRAY,LL*DIM,MX1MIN,ARRAY,LL*DIM,QY1MAX,ARRAY,LL*DIM,QY1MIN,ARRAY,LL*DIM,KK,ARRAY,LL+1*DIM,REACTION,ARRAY,6,LL+1*DIM,REACTIONMAX,ARRAY,6*DIM,REACTIONMIN,ARRAY,6*DO,I,1,LL+1,1/POST1SET,I!*****************竖向位移***************DO,J,1,LL+1,1*GET,DISP1(J,I),NODE,J,U,Y*ENDDO!*****************竖向弯矩Mz************DO,J,1,LL,1*GET,MZ1(J,I),ELEM, J, SMISC,6*ENDDO!*****************扭矩Mx************DO,J,1,LL,1*GET,MX1(J,I),ELEM, J, SMISC,4*ENDDO!*****************剪力Qy************DO,J,1,LL,1*GET,QY1(J,I),ELEM, J, SMISC,2*ENDDO!*****************支座反力************GET,REACTION(1,I),NODE, 2001, RF,FZ*GET,REACTION(2,I),NODE, 2002, RF,FZ*IF,e3,GT,0.0,THEN*GET,REACTION(3,I),NODE, 2003, RF,FZ*ELSE*GET,REACTION(3,I),NODE, L(1)+1, RF,FZ*ENDIF*IF,e4,GT,0.0,THEN*GET,REACTION(4,I),NODE, 2004, RF,FZ*ELSE*GET,REACTION(4,I),NODE, L(1)+L(2)+1, RF,FZ *ENDIF*GET,REACTION(5,I),NODE, 2005, RF,FZ*GET,REACTION(6,I),NODE, 2006, RF,FZKK(I)=I*ENDDO!**************************公路I级汽车活载影响线加载效应计算*******************************! 车道宽度假定为3.75m! 车道布置：横向NN=2车道，不考虑横向折减系数;跨径小于150m，不需要考虑纵向折减。
史上最全的ANSYS命令流大全
《史上最全的ANSYS命令流查询与解释》【1】*************************************************************************************对ansys主要命令的解释1， /PREP7 ! 加载前处理模块2， /CLEAR,NOSTART ! 清除已有的数据, 不读入启动文件的设置(不加载初始化文件)初始化文件是用于记录用户和系统选项设置的文本文件/CLEAR, START !清除系统中的所有数据,读入启动文件的设置/, EX10.5 ! 定义工程文件名称/TITLE, EX10.5 SOLID MODEL OF AN AXIAL BEARING ! 指定标题4， F,2,FY,-1000 ! 在2号节点上施加沿着-Y方向大小为1000N的集中力6， FINISH ! 退出模块命令7， /POST1 ! 加载后处理模块8， PLDISP,2 ! 显示结构变形图,参数“2”表示用虚线绘制出原来结构的轮廓9， ETABLE,STRS,LS,1 ! 用轴向应力SAXL的编号”LS,1”定义单元表STRSETABLE, MFORX,SMISC,1 ! 以杆单元的轴力为内容, 建立单元表MFORXETABLE, SAXL, LS, 1 ! 以杆单元的轴向应力为内容, 建立单元表SAXLETABLE, EPELAXL, LEPEL, 1 ! 以杆单元的轴向应变为内容, 建立单元表EPELAXLETABLE,STRS_ST,LS,1 !以杆件的轴向应力“LS,1”为内容定义单元表STRS_STETABLE, STRS_CO, LS,1 !以杆件的轴向应力“LS,1”定义单元表STRS_COETABLE,STRSX,S,X ! 定义X方向的应力为单元表STRSXETABLE,STRSY,S,Y ! 定义Y方向的应力为单元表STRSY*GET,STRSS_ST,ELEM,STEEL_E, ETAB, STRS_ST !从单元表STRS_ST中提取STEEL_E单元的应力结果，存入变量STRSS_ST;*GET, STRSS_CO,ELEM,COPPER_E,ETAB,STRS_CO”从单元表STRS_CO中提取COPPER_E单元的应力结果，存入变量STRSS_CO10 FINISH !退出以前的模块11, /CLEAR, START ! 清除系统中的所有数据,读入启动文件的设置12 /UNITS, SI !申明采用国际单位制14 /NUMBER, 2 !只显示编号, 不使用彩色/NUMBER, 0 ! 显示编号, 并使用彩色15 /SOLU ! 进入求解模块:定义力和位移边界条件,并求解ANTYPE, STATIC ! 申明分析类型是静力分析(STATIC或者0)OUTPR, BASIC, ALL ! 在输出结果中, 列出所有荷载步的基本计算结果OUTPR,BASIC,ALL !指定输出所有节点的基本数据OUTPR,BASIC,LAST ! 选择基本输出选项,直到最后一个荷载步OUTPR,,1 ! 输出第1个荷载步的基本计算结果OUTPR,BASIC,1 ! 选择第1荷载步的基本输出项目OUTPR,NLOAD,1 ! 指定输出第1荷载步的内容OUTRES,ALL,0 !设置将所有数据不记录到数据库。
曲线梁桥ANSYS计算命令流
!****************************************************************************** ***********************! case2:无偏载(以跨径布置30m+40m+30m,桥宽8.5为例）! 上海城市设计研究院L1+L2+L3预应力混凝土曲线连续梁桥结构分析! 两端为抗扭支座，中间支座为点铰支座! 每次要记得修改横隔梁的参数，即Mass21单元的实常数!****************************************************************************** ***********************FINI/CLE/prep7!DEFINE THE ELEMENTARY PARAMETERS*DIM,L,ARRAY,10*DIM,H,ARRAY,10*DIM,CITA,ARRAY,10!*****以下参数均可修改***************N=3 !跨数L(1)=30 !第一跨L(2)=40 !第二跨L(3)=30 !第三跨e1=1.25 !1#墩处内支座到中心线的间距e2=1.25 !1#墩处外支座到中心线的间距e3=0 !2#墩处的支座偏心距（正的表示外偏）e4=0 !3#墩处的支座偏心距e5=1.25 !4#墩处内支座到中心线的间距e6=1.25 !4#墩处外支座到中心线的间距R=10000 !曲线桥半径H0=1.0 !梁底到截面形心处的高度M=16146 !mass21单元质量J=27246.38 !mass21单元转动惯量!*************************************LL=0.0*DO,I,1,NLL=LL+L(I)CITA(I)=L(I)/R/3.1415925*180*ENDDOCITA0=LL/R/3.1415925*180LOCAL,11,1,0,0,0,0,270,0CSYS,11!DEFINE THE NODES OF BRIDGEN, 1, R, -179.9999, H0N, L(1)+1, R, -179.9999+CITA(1), H0N, L(1)+L(2)+1, R, -179.9999+CITA(1)+CITA(2), H0 N, L(1)+L(2)+L(3)+1, R, -179.9999+CITA0, H0FILL,1,L(1)+1FILL,L(1)+1,L(1)+L(2)+1FILL,L(1)+L(2)+1,L(1)+L(2)+L(3)+1!DEFINE THE REFERENCE POINT OF MAIN DECKN,2000,0,0,0!***************************************************************** !DEFINE THE NODES OF RIGID BEAM OF MAIN GIRDERN, 2001, R-e1, -179.9999, H0N, 2002, R+e2, -179.9999, H0N, 2003, R+e3, -179.9999+CITA(1), H0N, 2004, R+e4, -179.9999+CITA(1)+CITA(2), H0N, 2005, R-e5, -179.9999+CITA0, H0N, 2006, R+e6, -179.9999+CITA0, H0NROTA T,ALL!*****************************************************! 定义主梁单元材料、几何参数! 混凝土标号为C40!**************梁段截面几何特性*******************ET,1,BEAM4MP,EX , 1, 3.25E+10MP,NUXY, 1, 0.2000MP,DENS, 1, 3950.50MP,ALPX, 1, 1.00E-05R,1,4.39129,2.4428,20.7089,8.5,2.0,, !跨中梁断面RMORE ,,5.0388,,,,,!**************刚臂截面几何特性*******************! 增设抗扭或偏心支座时采用ET,3,BEAM4MP,EX , 3, 3.25E+10MP,NUXY, 3, 0.2000MP,DENS, 3, 0.0MP,ALPX, 3, 1.0E-08R,3,100,1000,1000,0.5,0.5,,RMORE ,,1000,,,,,!****************生成主梁单元**********************TYPE,1MA T,1REAL,1*DO,I,1,LL,1E,I,I+1,2000*ENDDO!*********************************************************** ! 生成主梁刚臂单元!*********************************************************** TYPE,3MA T,3REAL,3E,2001,1E,1,2002*IF,e3,GT,0.0,THENE,2003,L(1)+1*ELSENSEL,U, , ,2003*ENDIF*IF,e4,GT,0.0,THENE,2004,L(1)+L(2)+1*ELSENSEL,U, , ,2004*ENDIFE,2005,LL+1E,LL+1,2006!**************************************************! 支墩处的横隔梁采用Mass21单元来模拟!**************************************************ET,4,MASS21R,100, M ,M ,M ,J ,,,TYPE,4REAL,100E,1E,L(1)+1E,L(1)+L(2)+1E,LL+1!****施加约束条件**********D,2001,UX,,,,,UY,UZD,2002,UZ*IF,e3,GT,0.0,THEND,2003,,,,,,,UZ*ELSED,L(1)+1,,,,,,,UZ*ENDIF*IF,e4,GT,0.0,THEND,2004,,,,,,,UZ*ELSED,L(1)+L(2)+1,,,,,,,UZ*ENDIFD,2005,UX,,,,,,UZD,2006,UZ!********************************************************** ! 对桥梁桥面各节点逐点加载!********************************************************** /SOLUANTYPE,0 ! 静力分析ACEL,0,0,0K=1*DO,I,1,LL+1,1TIME,KFDELE,ALL,ALLF,I,FZ,-1000SOLVEK=K+1*ENDDO!******定义各片纵梁的竖向挠度、弯矩、扭矩、剪力的包络图******** *DIM,AMAX,ARRAY,LL+1*DIM,AMIN,ARRAY,LL+1*DIM,JZMAX,ARRAY,LL+1*DIM,JZMIN,ARRAY,LL+1*DIM,DISP1,ARRAY,LL+1,LL+1*DIM,MZ1,ARRAY,LL,LL+1*DIM,MX1,ARRAY,LL,LL+1*DIM,QY1,ARRAY,LL,LL+1*DIM,DISP1MAX,ARRAY,LL+1*DIM,DISP1MIN,ARRAY,LL+1*DIM,MZ1MAX,ARRAY,LL*DIM,MZ1MIN,ARRAY,LL*DIM,MX1MAX,ARRAY,LL*DIM,MX1MIN,ARRAY,LL*DIM,QY1MAX,ARRAY,LL*DIM,QY1MIN,ARRAY,LL*DIM,KK,ARRAY,LL+1*DIM,REACTION,ARRAY,6,LL+1*DIM,REACTIONMAX,ARRAY,6*DIM,REACTIONMIN,ARRAY,6*DO,I,1,LL+1,1/POST1SET,I!*****************竖向位移***************DO,J,1,LL+1,1*GET,DISP1(J,I),NODE,J,U,Y*ENDDO!*****************竖向弯矩Mz************DO,J,1,LL,1*GET,MZ1(J,I),ELEM, J, SMISC,6*ENDDO!*****************扭矩Mx************DO,J,1,LL,1*GET,MX1(J,I),ELEM, J, SMISC,4*ENDDO!*****************剪力Qy************DO,J,1,LL,1*GET,QY1(J,I),ELEM, J, SMISC,2*ENDDO!*****************支座反力************GET,REACTION(1,I),NODE, 2001, RF,FZ*GET,REACTION(2,I),NODE, 2002, RF,FZ*IF,e3,GT,0.0,THEN*GET,REACTION(3,I),NODE, 2003, RF,FZ*ELSE*GET,REACTION(3,I),NODE, L(1)+1, RF,FZ*ENDIF*IF,e4,GT,0.0,THEN*GET,REACTION(4,I),NODE, 2004, RF,FZ*ELSE*GET,REACTION(4,I),NODE, L(1)+L(2)+1, RF,FZ *ENDIF*GET,REACTION(5,I),NODE, 2005, RF,FZ*GET,REACTION(6,I),NODE, 2006, RF,FZKK(I)=I*ENDDO!**************************公路I级汽车活载影响线加载效应计算*******************************! 车道宽度假定为3.75m! 车道布置：横向NN=2车道，不考虑横向折减系数;跨径小于150m，不需要考虑纵向折减。
连续梁桥ansys命令流
!平面杆件结构静力分析!==========前处理==========FINISH !退出PREP7处理器/CLEAR !清除当前数据/Prep7 !进入prep7处理器!创建线，赋予线属性*DO,I,1,5,1 !连续定义5个单元类型参考号，单元类型为BEAM44ET,I,BEAM44*ENDDOET,6,Link8 !定义第6个单元类型参考号，单元类型为Link8KEYOPT,3,7,1 !定义单元左端绕Z轴角位移自由度释放类型号为3KEYOPT,5,8,1 !定义单元右端绕Z轴角位移自由度释放类型号为5MP,EX,1,200000000 !定义弹性模量MP,prxy,1,0.16 !定义泊松比MP,DENS,1,78 !定义单元容重K,1,77.7875,202.6708 !定义关键点K,2,88.1393,202.6708 !定义关键点K,3,98.4911,202.6708 !定义关键点K,4,108.843,202.6708 !定义关键点K,5,119.1948,202.6708 !定义关键点K,6,129.5466,202.6708 !定义关键点K,7,139.8984,202.6708 !定义关键点K,8,150.2503,202.6708 !定义关键点K,9,160.6021,202.6708 !定义关键点R,1,0.00485883,0.00022855,7.4E-6,0.144,0.4 !定义实常数LSEL,none !默认选前一线段L,1,2 !连接两关键点为直线LATT,1,1,1 !设置上面连接线的单元属性：包括材料属性、实常数及BEAM44单元类型R,2,0.00485883,0.00022855,7.4E-6,0.144,0.4 !定义实常数LSEL,none !默认选前一线段L,2,3 !连接两关键点为直线LATT,1,2,1 !设置上面连接线的单元属性：包括材料属性、实常数及BEAM44单元类型R,3,0.00485883,0.00022855,7.4E-6,0.144,0.4 !定义实常数LSEL,none !默认选前一线段L,3,4 !连接两关键点为直线LATT,1,3,1 !设置上面连接线的单元属性：包括材料属性、实常数及BEAM44单元类型R,4,0.00485883,0.00022855,7.4E-6,0.144,0.4 !定义实常数LSEL,none !默认选前一线段L,4,5 !连接两关键点为直线LATT,1,4,1 !设置上面连接线的单元属性：包括材料属性、实常数及BEAM44单元类型R,5,0.00485883,0.00022855,7.4E-6,0.144,0.4 !定义实常数LSEL,none !默认选前一线段L,5,6 !连接两关键点为直线LATT,1,5,1 !设置上面连接线的单元属性：包括材料属性、实常数及BEAM44单元类型R,6,0.00485883,0.00022855,7.4E-6,0.144,0.4 !定义实常数LSEL,none !默认选前一线段L,6,7 !连接两关键点为直线LATT,1,6,1 !设置上面连接线的单元属性：包括材料属性、实常数及BEAM44单元类型R,7,0.00485883,0.00022855,7.4E-6,0.144,0.4 !定义实常数LSEL,none !默认选前一线段L,7,8 !连接两关键点为直线LATT,1,7,1 !设置上面连接线的单元属性：包括材料属性、实常数及BEAM44单元类型R,8,0.00485883,0.00022855,7.4E-6,0.144,0.4 !定义实常数LSEL,none !默认选前一线段L,8,9 !连接两关键点为直线LATT,1,8,1 !设置上面连接线的单元属性：包括材料属性、实常数及BEAM44单元类型!==========划分网格==========LSEL,ALL !选择所有线段LESIZE,ALL,,,1 !网格划分定义：一条线段一个单元LMESH,ALL !所有线段进行网格划分，形成单元!==========定义约束==========K,10,77.7875,202.6708 !定义关键点DK,10,UX,,,,UY,UZ !这是一个固定铰支座，节点的X、Y、Z三个方向的线位移被约束K,11,160.6021,202.6708 !定义关键点DK,11,UY,,,,UZ !这是一个竖向活动支座，Y向位移被约束，同时也约束了Z向的线位移!==========定义荷载==========SFBEAM,1,2,PRES,50,50 !施加单元Y向线荷载SFBEAM,2,2,PRES,50,50 !施加单元Y向线荷载SFBEAM,3,2,PRES,50,50 !施加单元Y向线荷载SFBEAM,4,2,PRES,50,50 !施加单元Y向线荷载SFBEAM,5,2,PRES,50,50 !施加单元Y向线荷载SFBEAM,6,2,PRES,50,50 !施加单元Y向线荷载SFBEAM,7,2,PRES,50,50 !施加单元Y向线荷载SFBEAM,8,2,PRES,50,50 !施加单元Y向线荷载NUMMRG,KP,0.005 !对所有关键点进行合并!以下为释放自由度ALLSEL,ALL !选择全部实体/SOLU !进入solu求解器SOLVE !结构力学求解计算/POST1 !进入后处理器PLDISP,1 !显示变形ETABLE,FXI,SMISC,1 !定义单元左端轴力表ETABLE,FXJ,SMISC,7 !定义单元右端轴力表ETABLE,FYI,SMISC,2 !定义单元左端Y向剪力表ETABLE,FYJ,SMISC,8 !定义单元右端Y向剪力表ETABLE,FZI,SMISC,3 !定义单元左端Z向剪力表ETABLE,FZJ,SMISC,9 !定义单元右端Z向剪力表ETABLE,MXI,SMISC,4 !定义单元左端绕X轴钮矩表ETABLE,MXJ,SMISC,10 !定义单元右端绕X轴钮矩表ETABLE,MYI,SMISC,5 !定义单元左端绕Y轴弯矩表ETABLE,MYJ,SMISC,11 !定义单元右端绕Y轴弯矩表ETABLE,MZI,SMISC,6 !定义单元左端绕Z轴弯矩表ETABLE,MZJ,SMISC,12 !定义单元右端绕Z轴弯矩表ETABLE,SMINI,NMISC,2ETABLE,SMINJ,NMISC,4PLLS,FXI,FXJ !显示轴力PLLS,FYI,FYJ !显示剪力PLLS,MXI,MXJ !显示扭矩PLLS,MYI,MYJ !显示绕Y轴的弯矩PLLS,SMINI,SMINJPLLS,MZI,MZJ !显示单元绕Z轴的弯矩SECTYPE, 1, BEAM, CSOLID SECOFFSET, CENTSECDATA,1。
梁ansys分析实例讲解课件
ansys软件操作流程
定义弹性模量、泊松比、密度等材料属性 建立几何体
在Model模块中选择“Model”选项卡
ansys软件操作流程
• 创建梁的几何体，输入梁的截面尺寸等信息
ansys软件操作流程
网格划分 在Model模块中选择“Mesh”选项卡 设置网格大小、网格类型等参数
定义约束和载荷
在简支梁的两个端点上定义约束和 载荷。
边界条件与载荷施加
固定约束
在简支梁的两个端点施加固定约 束，以模拟简支边界条件。
均布载荷
在简支梁的跨中施加均布载荷， 以模拟简支梁受到的集中力。
网格划分与求解
网格划分
对简支梁进行网格划分，可以选 择合适的网格密度以提高求解精
度。
求解设置
在“Solution”菜单下进行求解 设置，包括迭代次数、收敛准则
减少模型规模：减小 模型的规模，降低内 存需求；
如何设置ANSYS软 件中的单位制？
升级硬件：增加物理 内存或使用更好的计 算机配置。
a型时， 选择合适的单位制，如米制或英
制；
在模型树中选择模型名称，进入 Model模块，在Model模块中选 择“Model”选项卡，在弹出的 对话框中选择“Units”选项卡
求解设置
进行求解设置，包括迭代次数、收敛标准等。
求解过程
进行求解，得到桥梁结构的应力分布、位移分布等结果。
06
ansys软件操作流程及常见问题 解答
ansys软件操作流程
建立模型 启动ANSYS软件，选择Workbench模式
创建新的模型文件，命名并保存
ansys软件操作流程
Ansys命令流大全（整理）
Ansys命令流大全（整理）1、A,P1,P2,P3,P4,P5,P6,P7,P8,P9此命令用已知的一组关键点点（P1~P9）来定义面（Area），最少使用三个点才能围成面，同时产生转围绕些面的线。
点要依次序输入，输入的顺序会决定面的法线方向。
如果超过四个点，则这些点必须在同一个平面上。
Menu Paths:Main Menu>Preprocessor>Create>Arbitrary>Through KPs2、*ABBR，Abbr,String－－定义一个缩略语．Abbr:用来表示字符串＂String＂的缩略语，长度不超过８个字符．String：将由＂Abbr＂表示的字符串，长度不超过６０个字符．3、ABBRES，Lab，Fname，Ext－从一个编码文件中读出缩略语．Lab：指定读操作的标题，NEW：用这些读出的缩略语重新取代当前的缩略语（默认）CHANGE：将读出的缩略语添加到当前缩略语阵列，并替代现存同名的缩略语．Ext：如果＂Fname＂是空的，则缺省的扩展命是＂ABBR＂．4、ABBSA V，Lab，Fname，Ext－将当前的缩略语写入一个文本文件里Lab：指定写操作的标题，若为ALL，表示将所有的缩略语都写入文件（默认）5、add, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic：变量号name: 变量的名称6、Adele,na1,na2,ninc,kswp ！kswp=0时只删除掉面积本身，＝1时低单元点一并删除。
7、Adrag, nl1,nl2,nl3,nl4,nl5,nl6, nlp1,nlp2,nlp3,nlp4,nlp5,nlp6 ！面积的建立，沿某组线段路径，拉伸而成。
8、Afillt,na1,na2,rad ！建立圆角面积，在两相交平面间产生曲面，rad为半径。
(完整版)ANSYS命令流总结(全)
ANSYS结构解析单元功能与特征/POST1/可以构成一一些命令，一般是一种整体命令（ session),三十也有特别，比方是办理 ! 是说明说明符号，，与其余软件的说明是相同的， ansys 不作为命令读取，*此符号一般是 APDL 的表记符，也就是 ansys 的参数化语言，如 *do ,,,*enddo 等等NSEL 的意思是node select，即选择节点。
s 就是 select，选择。
DIM是定义数组的意思。
array 数组。
MP 命令用来定义资料参数。
K 是建立要点点命令。
K, 要点点编号 ,x 坐标 ,y 坐标， z 坐标。
K, NPT, X, Y , Z 是定义要点点， K 是命令， NPT 是要点点编号， XYZ 是坐标。
NUMMRG , keypoint 用这个命令，要保证要点点的地点完整相同，不过要点点号不一样样的才行。
这个命令关于重复的线面都可以用。
这个很简单，压缩要点。
Ngen 复制节点e,节点号码：这个命令式经过节点来形成单元NUMCMP,ALL ：压缩所有编号，这样你所有的线都会挨次次重新编号 ~你若是需要固定的线固定的标号NSUBST,100,500,50 ：经过指定子步数来设置载荷步的子步LNSRCH 线性搜寻是求解非线性代数方程组的一种技巧，此法会在一段区间内，以必定的步长逐渐搜寻根，对比常用的牛顿迭代法所要耗费的计算量大得多，但它可以防备在一些状况下牛顿迭代法出现的跳跃现象。
LNSRCH激活线性搜寻PRED 激活自由度求解展望NEQIT 指定一个荷载步中的最大子步数AUTOTS自动求解控制打开自动时间步长.KBC -指定阶段状也许用跳板装载里面一个负荷步骤。
SPLINE :P1， P2， P3，P4， P5， P6， XV1 ， YV1 ， ZV1 ， XV6 ，YV6 ， ZV6 （生成分段样条曲线）*DIM ， Par，Type ，IMAX ，JMAX ， KMAX ， Var1，Var2， Var3（定义载荷数组的名称）【注】 Par: 数组名Type： array 数组，仿佛fortran, 下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8 个字符）tableIMAX ， JMAX ， KMAX各维的最大下标号Var1， Var2，Var3 各维变量名，缺省为row,column,plane( 当 type 为 table 时 )/config 是设置 ansys 配置参数的命令格式为 /CONFIG, Lab, V ALUELab 为参数名称value 为参数值比方： /config ， MXEL ，10000 的意思是最大单元数为10000杆单元 : LINK1、 8、 10、 11、 180梁单元： BEAM3、 4、 23、 24，44， 54， 188， 189管单元 : PIPE16， 17， 18， 20， 59， 602D实体元 : PLANE2， 25， 42， 82， 83， 145，146， 182， 1833D实体元 : SOLID45， 46， 64，65， 72， 73，92， 95， 147，148， 185， 186，187， 191壳单元 : SHELL28， 41， 43， 51， 61， 63， 91， 93， 99， 143， 150， 181，208， 209弹簧单元 : COMBIN7， 14， 37，39， 40质量单元 : MASS21接触单元 : CONTAC12， 52， TARGE169， 170， CONTA171， 172， 173， 174， 175， 178矩阵单元 : MATRIX27， 50表面效应元 : SURF153， 154粘弹实体元 : VISCO88， 89， 106， 107， 108，超弹实体元 : HYPER56， 58， 74， 84， 86， 158耦合场单元 : SOLID5， PLANE13， FLUID29， 30，38， SOLID62， FLUID79， FLUID80，81，SOLID98， FLUID129， INFIN110 ， 111， FLUID116，130界面单元 : INTER192， 193， 194， 195显式动力解析单元 : LINK160， BEAM161， PLANE162， SHELL163， SOLID164， COMBI16杆单元单元名称简称节点数节点自由度特征备注LINK12D杆2Ux,Uy EPCSDGB常用杆元LINK83D杆Ux,Uy,Uz EPCSDGBLINK103D仅受拉EDGB模拟缆索的废弛及或仅受压杆缝隙LINK113D线性调理EGB模拟液压缸和大转器动LINK1803D有限应变杆EPCDFGB另可考虑粘弹塑性E- 弹性 (Elasticity)，P-塑性(Plasticity)，C-蠕变(Creep)，S-膨胀(Swelling)，D-大变形或大挠度deflection)， F- 大应变 (Large strain)或有限应变(Finite strain)，B-单元存亡(Birth and dead),G-化 (Stress stiffness)或几何刚度(Geometric stiffening)，A-自适应降落(Adaptive descent)等。
用ANSYS进行桥梁结构分析
用ANSYS进行桥梁结构分析谢宝来华龙海引言：我院现在进行桥梁结构分析主要用桥梁博士和BSACS，这两种软件均以平面杆系为计算内核，多用来解决平面问题。
近来偶然接触到ANSYS，发现其结构分析功能强大，现将一些研究心得写出来，并用一个很好的学习例子（空间钢管拱斜拉桥）作为引玉之砖，和同事们共同研究讨论，共同提高我院的桥梁结构分析水平而努力。
【摘要】本文从有限元的一些基本概念出发，重点介绍了有限元软件ANSYS平台的特点、使用方法和利用APDL语言快速进行桥梁的结构分析，最后通过工程实例来更近一步的介绍ANSYS进行结构分析的一般方法，同时进行归纳总结了各种单元类型的适用范围和桥梁结构分析最合适的单元类型。
【关键词】ANSYS有限元APDL结构桥梁工程单元类型一、基本概念有限元分析(FEA)是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。
还利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。
有限元模型是真实系统理想化的数学抽象。
真实系统有限元模型自由度(DOFs)用于描述一个物理场的响应特性。
节点和单元1、每个单元的特性是通过一些线性方程式来描述的。
2、作为一个整体，单元形成了整体结构的数学模型。
3、信息是通过单元之间的公共节点传递的。
4、节点自由度是随连接该节点单元类型变化的。
单元形函数1、FEA 仅仅求解节点处的DOF 值。
2、单元形函数是一种数学函数，规定了从节点DOF 值到单元内所有点处DOF 值的计算方法。
3、因此，单元形函数提供出一种描述单元内部结果的“形状”。
4、单元形函数描述的是给定单元的一种假定的特性。
5、单元形函数与真实工作特性吻合好坏程度直接影响求解精度。
6、DOF 值可以精确或不太精确地等于在节点处的真实解，但单元内的平均值与实际情况吻合得很好。
7、这些平均意义上的典型解是从单元DOFs 推导出来的（如，结构应力，热梯度）。
【ANSYS算例】3.3.7(3) 三梁平面框架结构的有限元分析(GUI)及命令流
【ANSYS 算例】3.3.7(3) 三梁平面框架结构的有限元分析针对【典型例题】3.3.7(1)的模型，即如图3-19所示的框架结构，其顶端受均布力作用，用有限元方法分析该结构的位移。
结构中各个截面的参数都为：113.010Pa E ，746.510m I ，426.810m A ，相应的有限元分析模型见图3-20。
在ANSYS 平台上，完成相应的力学分析。
图3-19 框架结构受一均布力作用（a ） 节点位移及单元编号 （b ） 等效在节点上的外力图3-20 单元划分、节点位移及节点上的外载解答 对该问题进行有限元分析的过程如下。
1．基于图形界面的交互式操作(step by step)(1) 进入ANSYS(设定工作目录和工作文件)程序 →ANSYS → ANSYS Interactive →Working directory (设置工作目录) →Initial jobname (设置工作文件名): beam3→Run → OK(2) 设置计算类型ANSYS Main Menu: Preferences… → Structural → OK(3) 选择单元类型ANSYS Main Menu: Preprocessor →Element Type →Add/Edit/Delete… →Add… →beam ：2D elastic 3 →OK (返回到Element Types 窗口) →Close(4) 定义材料参数ANSYS Main Menu:Preprocessor →Material Props →Material Models→Structural →Linear →Elastic→Isotropic: EX:3e11 (弹性模量) →OK →鼠标点击该窗口右上角的“ ”来关闭该窗口(5) 定义实常数以确定平面问题的厚度ANSYS Main Menu: Preprocessor →Real Constant s… →Add/Edit/Delete →Add →Type 1 Beam3→OK→Real Constant Set No: 1 (第1号实常数), Cross-sectional area:6.8e-4 (梁的横截面积) →OK →Close(6) 生成几何模型生成节点ANSYS Main Menu: Preprocessor →Modeling →Creat→Nodes→In Active CS→Node number 1 →X:0,Y:0.96,Z:0 →Apply→Node number 2 →X:1.44,Y:0.96,Z:0 →Apply→Node number 3 →X:0,Y:0,Z:0→Apply→Node number 4 →X:1.44,Y:0,Z:0→OK生成单元ANSYS Main Menu: Preprocessor →Modeling →Create →Element →Auto Numbered →Thru Nodes →选择节点1，2(生成单元1)→apply →选择节点1，3(生成单元2)→apply →选择节点2，4(生成单元3)→OK(7)模型施加约束和外载左边加X方向的受力ANSYS Main Menu:Solution →Define Loads →Apply →Structural →Force/Moment →On Nodes →选择节点1→apply →Direction of force: FX →V ALUE：3000 →OK→上方施加Y方向的均布载荷ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Beams →选取单元1(节点1和节点2之间)→apply →V ALI：4167→V ALJ：4167→OK左、右下角节点加约束ANSYS Main Menu:Solution →Define Loads →Apply →Structural →Displacement →On Nodes →选取节点3和节点4 →Apply →Lab:ALL DOF →OK(8) 分析计算ANSYS Main Menu:Solution →Solve →Current LS →OK →Should the Solve Command be Executed? Y→Close (Solution is done! ) →关闭文字窗口(9) 结果显示ANSYS Main Menu: General Postproc →Plot Results →Deformed Shape … →Def + Undeformed →OK (返回到Plot Results)(10) 退出系统ANSYS Utility Menu: File→Exit …→Save Everything→OK(11) 计算结果的验证与MA TLAB支反力计算结果一致。
ansys桥梁受力分析命令流
ANSYS桥梁受力分析命令流ANSYS是一款广泛应用于工程领域的有限元分析软件，它可以帮助工程师、设计师、科学家等进行各种力学或结构分析。
在工程领域中，桥梁是一种比较常见的结构，因此本文将介绍在ANSYS中进行桥梁受力分析的命令流。
1. 建立桥梁模型首先，我们需要建立桥梁模型。
在ANSYS中，可以使用多种模型建立工具，例如现在比较流行的三维CAD软件，如SolidWorks、CATIA等，均可以建立出桥梁模型。
建立好桥梁模型后，需要导入到ANSYS中进行后续分析。
2. 分配材料属性桥梁模型建立完后，需要根据实际情况为它分配材料属性。
以混凝土桥梁为例，我们可以使用“MP（材料属性）”命令为桥梁定义混凝土的相应属性。
比如，定义混凝土的弹性模量、泊松比等。
3. 设定支座约束桥梁的受力分析需要考虑桥梁的支座约束情况。
在ANSYS中，我们可以通过使用“SUPPORT”命令设定不同类型的支座约束。
4. 设定荷载桥梁承受的荷载对于力学分析至关重要，因此在进行桥梁受力分析时，必须设定荷载的种类和大小。
ANSYS提供了多种荷载设定方式，如集中荷载、分布荷载、自重荷载等。
5. 载荷应用载荷应用是桥梁受力分析中的一个关键步骤，它能够准确模拟桥梁所受荷载。
在ANSYS中，可以使用“D”命令定义荷载的应用方式，如定义一个位移荷载、力荷载等。
6. 求解过程ANSYS中提供了多种求解器，可以快速准确地解析桥梁模型受力情况。
在进行桥梁受力分析时，需要选择合适的求解器并进行计算。
ANSYS中提供了两种类型求解器，一种是基于传统矩阵进行计算的Direct Solver求解器，另一种是基于有限元模型优化的Iterative Solver求解器。
7. 结果分析在桥梁受力分析计算完成后，需要对应力和位移等结果进行分析并产生可视化的输出。
在ANSYS中，可以使用“POST1”命令来输出结果，并对结果进行后处理。
8.经过以上步骤进行受力分析后，可以得到桥梁受力的详细情况，也可以对桥梁的结构进行合理的优化。
ANSYS命令流详解（超全）
ANSYS命令流详解（超全）一、定义材料号及特性mp,lab, mat, co, c1,…….c4lab: 待定义的特性项目（ex,alpx,reft,prxy,nuxy,gxy,mu,dens）ex: 弹性模量nuxy: 小泊松比alpx: 热膨胀系数reft: 参考温度reft: 参考温度prxy: 主泊松比gxy: 剪切模量mu: 摩擦系数dens: 质量密度mat: 材料编号（缺省为当前材料号）c 材料特性值，或材料之特性，温度曲线中的常数项c1-c4: 材料的特性-温度曲线中1次项，2次项，3次项，4次项的系数二、定义DP材料：首先要定义EX和泊松比：MP，EX，MAT，……MP，NUXY，MAT，……定义DP材料单元表（这里不考虑温度）：TB，DP，MAT进入单元表并编辑添加单元表：TBDATA，1，CTBDATA，2，ψTBDATA，3，……如定义：EX=1E8，NUXY=0.3，C=27，ψ=45的命令如下：MP，EX，1，1E8MP，NUXY，1，0.3TB，DP，1TBDATA，1，27TBDATA，2，45这里要注意的是，在前处理的最初，要将角度单位转化到“度”，即命令：*afun,deg三、单元生死载荷步!第一个载荷步TIME,... !设定时间值（静力分析选项）NLGEOM,ON !打开大位移效果NROPT,FULL !设定牛顿－拉夫森选项ESTIF,... !设定非缺省缩减因子（可选）ESEL,... !选择在本载荷步中将不激活的单元EKILL,... !不激活选择的单元ESEL,S,LIVE !选择所有活动单元NSLE,S !选择所有活动结点NSEL,INVE !选择所有非活动结点（不与活动单元相连的结点）D,ALL,ALL,0 !约束所有不活动的结点自由度（可选）NSEL,ALL !选择所有结点ESEL,ALL !选择所有单元D,... !施加合适的约束F,... !施加合适的活动结点自由度载荷SF,... !施加合适的单元载荷BF,... !施加合适的体载荷SAVESOLVE请参阅TIME,NLGEOM,NROPT,ESTIF,ESEL,EKILL,NSLE,NSEL,D,F,SF和BF命令得到更详细的解释。
基于ANSYS 的连续刚构桥分析(操作篇 )
目录一、工程背景 (1)二、工程模型 (1)三、ANSYS分析 (2)(一)前处理 (2)(1)定义单元类型 (2)(2)定义材料属性 (3)(3)建立工程简化模型 (3)(4)有限元网格划分 (5)(二)模态分析 (5)(1)选择求解类型 (5)(2)建立边界条件 (6)(3)输出设置 (6)(4)求解 (6)(5)读取结果 (6)(6)结果分析 (8)(三)结构试验载荷分析 (8)(1)第二跨跨中模拟车载分析 (8)(2)边跨跨中模拟车载分析 (9)四、结果分析与强度校核 (10)(一)结果分析 (10)(二)简单强度校核 (10)参考文献 (11)连续刚构桥分析一、工程背景：随着我国经济的发展，对交通运输的要求也不断提高；高速路，高铁线等遍布全国，这就免不了要架桥修路。
截至2014年年底，我国公路桥梁总数已达75.71万座，4257.89万延米i。
进百万的桥梁屹立在我国交通线上，其安全便是头等大事。
随着交通运输线的再扩大，连续刚构桥跨越能力大，施工难度小，行车舒顺，养护简便，造价较低等优点将被广泛应用。
二、工程模型：现有某预应力混凝土连续刚构桥，桥梁全长为184m，宽13m，其中车行道宽11.5m，两侧防撞栏杆各0.75m主梁采用C50混凝土。
桥梁设计载荷为公路—— 级。
图2-1桥梁侧立面图上部结构为48m+88m+48m三跨预应力混凝土边界面连续箱梁。
箱梁为单箱双室箱形截面，箱梁根部高5m，中跨梁高2.2m，边跨梁端高2.2m。
箱梁顶板宽12.7m,底板宽8.7m,翼缘板悬臂长2.0m，箱梁高度从距墩中心3.0m处到跨中合龙段处按二次抛物线变化。
0号至3号块长3m（4x3m)，4、5号块长3.5m(2x3.5m)，6号块到合龙段长4m(6x4m)，合龙段长2m。
边跨端部设1.5m横隔板，墩顶0号块设两道厚1.2m横隔板。
0号块范围内箱梁底板厚度为0.90m，1号块范围内底板厚度由0.90m线性变化到0.557m，2号块到合龙段范围内底板厚度由0.557m 线性变化到0.3m。
个人总结ansys命令流
Q235 属性：弹性模量E=2.1e5 N/mm2密度=7.85e-6kg/mm3泊松比=0.3mp,ex,1,2.1e5mp,prxy,1,0.3mp,dens,1,7.85e-61,ksymm 镜像点2,arsym 镜像面3,kgen 复制点4.adele删除面6，kdist,k1,k2 测量两关键点的距离7,adele,a,,,1 删除area and below8,创建圆柱面：circle 创建圆然后创建直线然(轴线)利用拉伸命令创建圆柱面creat__areas__by Linesadrag线拉伸成面modeling>operate>extrude>lines>>along linesVDRAG 面拉伸成体modeling>operate>extrude>areas>>along lines！创建空心圆柱体CYLIND, RAD1, RAD2, Z1, Z2, THETA1, THETA2这个命令Main Menu>Preprocessor>Modeling>Create>Volumes>Cylinder>By DimensionsMain Menu>Preprocessor>Trefftz Domain>TZ Geometry>Create>Volume>Cylinder>By Dimensions9,aptn 分割面10,asbw 用工作平面切割面11.wpoffs12.wprotang 过圆外一点做圆的切线(0°或180°)14，nummrg将重复的点消除15，asba面减去面16,两个圆柱面的相贯线作法：做出两个相穿的圆柱面，利用APTN命令17，选择面，不选择一部分面asel,u,loc,z,kz(735)18.在工作平面上生成一个矩形面RECTING,X1,X2,Y1,Y2X1，X2——矩形在工作平面X方向坐标值的变化范围Y1,Y2——矩形在工作平面Y方向坐标值的变化范围18,圆阵列建立工作平面与圆柱的横截面平行，在工作平面情况下建立局部坐标系（柱坐标系），然后利用agen命令复制。
ANSYS算例 题斜支座的处理及命令流
ANSYS算例题斜支座的处理及命令流【ANSYS算例】5.3(8)平面问题斜支座的处理如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45º的斜支座，试用一个4节点矩形单元分析该结构的位移场。
(a)平面结构(b)有限元分析模型图5-7 带斜支座的平面结构基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。
解答给出的操作过程及命令流如下。
1 基于图形界面的交互式操作(step by step)(1) 进入ANSYS（设定工作目录和工作文件）程序→ANSYS→ANSYS Interactive→Working directory（设置工作目录）→Initial jobname:Plane support(设置工作文件名)→Run→OK(2) 设置计算类型ANSYS Main Menu：Preferences…→Structural→OK(3) 定义单元类型ANSYS Main Menu：Preprocessor →Element Type→Add/Edit/Delete... →Add…→Solid: Quad 4node 42→OK（返回到Element Types窗口）→Type 1 PLANE42 →Options…→K3:Plane strs w/thk(带厚度的平面应力问题)→OK→Close(4) 定义实常数ANSYS Main Menu：Preprocessor →Real Constants…→Add…→Type 1 PLANE42 →OK→输入Real Constants Set No: 1,THK:0.1 →Close (关闭Real Constants 窗口)(5) 定义材料参数ANSYS Main Menu：Preprocessor→Material Props →Material Models→Structural→Linear →Isotropic→EX: 2E5, PRXY:0.25(定义泊松比及弹性模量) →OK →Close（关闭材料定义窗口）(6) 构造平面模型生成节点ANSYS Main Menu：Preprocessor →Modeling→Create→Nodes→In Active CS→Node number:1，X，Y，Z Location in active CS:2，2，0 →Apply →同样依次输入其余3个节点坐标（最左端为起始点，坐标分别为（0，2，0）、（0，0，0）、（2，0，0））（若采用【典型例题】5.3(7)中直接法（2），可将4号节点THXY设置为斜方向45°，然后直接对该节点使用UY=0的约束即可。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
/prep7
!DEFINE THE ELEMENTARY PARAMETERS
*DIM,L,ARRAY,10
*DIM,DISTC,ARRAY,10
N=4
NBOX=2
L(1)=20
L(2)=30
L(3)=30
L(4)=20
NUML=NBOX+3
NUMC=(L(1)+L(2)+L(3)+L(4))+1
CITA=25/180*3.1415926
DISTC(1)=2.9584
DISTC(2)=4.0802
DISTC(3)=0.00
DISTC(4)=4.0802
DISTC(5)=2.9584
!DEFINE THE NODES OF BRIDGE
N, 1, 0, 0, 0.00
N, NUMC, NUMC-1, 0, 0.00
FILL,1,NUMC
N, 200+1,
(DISTC(1)+DISTC(2))*SIN(CITA)/COS(CITA), 0.0,-(DISTC(1)+DISTC(2))
N,
200+NUMC,NUMC-1+(DISTC(1)+DISTC(2))* SIN(CITA)/COS(CITA),0.0,-(DISTC(1)+DISTC(2)) FILL,201,200+NUMC
N, 400+1, DISTC(2)*SIN(CITA)/COS(CITA), 0.0, -DISTC(2)
N,
400+NUMC,NUMC-1+DISTC(2)*SIN(CITA)/ COS(CITA), 0.0,-DISTC(2)
FILL,401,400+NUMC
N, 600+1, -DISTC(2)*SIN(CITA)/COS(CITA), 0.0, DISTC(2)
N,
600+NUMC,NUMC-1-DISTC(4)*SIN(CITA)/C OS(CITA),0.0, DISTC(4)
FILL,601,600+NUMC
N, 800+1,
-(DISTC(4)+DISTC(5))*SIN(CITA)/COS(CITA ),0.0, DISTC(4)+DISTC(5) N,
800+NUMC,NUMC-1-(DISTC(4)+DISTC(5))* SIN(CITA)/COS(CITA),0.0, DISTC(4)+DISTC(5) FILL,801,800+NUMC
! 定义纵梁单元材料、几何参数
ET,1,BEAM4
MP,EX , 1, 3.45E+10
MP,NUXY, 1, 0.2000
MP,DENS, 1, 2600
MP,ALPX, 1, 1.00E-05
R,1,0.3133,0.0871,0.0959,1.0,0.25,,
RMORE ,,0.00265,,,,,
ET,2,BEAM4
MP,EX , 2, 3.45E+10
MP,NUXY, 2, 0.2000
MP,DENS, 2, 2600
MP,ALPX, 2, 1.00E-05
R,2,1.6048,0.5486,1.5400,2.0,1.0,,
RMORE ,,0.036,,,,,
ET,3,BEAM4
MP,EX , 3, 3.45E+10
MP,NUXY, 3, 0.2000
MP,DENS, 3, 2600
MP,ALPX, 3, 1.00E-05
R,3,2.3625,0.9885,2.2600,2.0,1.0,,
RMORE ,,0.0956,,,,
! 定义横梁单元材料、几何参数
ET,4,BEAM4
MP,EX , 4, 3.45E+10
MP,NUXY, 4, 0.2000
MP,DENS, 4, 0
MP,ALPX, 4, 1.00E-05
R,4,0.500,0.4582,0.0416,1.0,1.0,,
RMORE ,,0.00879,,,,,
ET,5,BEAM4
MP,EX , 5, 3.45E+10
MP,NUXY, 5, 0.2000
MP,DENS, 5, 0
MP,ALPX, 5, 1.00E-05
R,5,1.0250,0.2770,0.0495,1.0,1.0,,
RMORE ,,0.0609,,,,,
! 生成纵向梁格单元
TYPE,1
MAT,1
REAL,1
*DO,I,1,NUMC-1,1 E,200+I,200+I+1 *ENDDO
*DO,I,1,NUMC-1,1 E,800+I,800+I+1 *ENDDO
TYPE,2
MAT,2
REAL,2
*DO,I,1,NUMC-1,1 E,400+I,400+I+1 *ENDDO
*DO,I,1,NUMC-1,1 E,600+I,600+I+1 *ENDDO
TYPE,3
MAT,3
REAL,3
*DO,I,1,NUMC-1,1 E,I,I+1
*ENDDO
! 生成横向梁格单元TYPE,5
MAT,5
REAL,5
E,200+1,400+1 E,400+1,1
E,1,600+1
E,600+1,800+1 TYPE,4
MAT,4
REAL,4
*DO,I,2,L(1),1
E,200+I,400+I
E,400+I,I
E,I,600+I
E,600+I,800+I *ENDDO
TYPE,5
MAT,5
REAL,5
E,200+L(1)+1,400+L(1)+1
E,400+L(1)+1,L(1)+1
E,L(1)+1,600+L(1)+1
E,600+L(1)+1,800+L(1)+1
TYPE,4
MAT,4
REAL,4
*DO,I,L(1)+2,L(1)+L(2),1
E,200+I,400+I
E,400+I,I
E,I,600+I
E,600+I,800+I
*ENDDO
TYPE,5
MAT,5
REAL,5
E,200+L(1)+L(2)+1,400+L(1)+L(2)+1
E,400+L(1)+L(2)+1,L(1)+L(2)+1
E,L(1)+L(2)+1,600+L(1)+L(2)+1
E,600+L(1)+L(2)+1,800+L(1)+L(2)+1
TYPE,4
MAT,4
REAL,4
*DO,I,L(1)+L(2)+2,L(1)+L(2)+L(3),1
E,200+I,400+I
E,400+I,I
E,I,600+I
E,600+I,800+I
*ENDDO
TYPE,5
MAT,5
REAL,5
E,200+L(1)+L(2)+L(3)+1,400+L(1)+L(2)+L(3)+1 E,400+L(1)+L(2)+L(3)+1,L(1)+L(2)+L(3)+1
E,L(1)+L(2)+L(3)+1,600+L(1)+L(2)+L(3)+1
E,600+L(1)+L(2)+L(3)+1,800+L(1)+L(2)+L(3)+1 TYPE,4
MAT,4
REAL,4
*DO,I,L(1)+L(2)+L(3)+2,L(1)+L(2)+L(3)+L(4),1 E,200+I,400+I
E,400+I,I
E,I,600+I
E,600+I,800+I
*ENDDO
TYPE,5
MAT,5
REAL,5
E,200+L(1)+L(2)+L(3)+L(4)+1,400+L(1)+L(2)+L(3) +L(4)+1
E,400+L(1)+L(2)+L(3)+L(4)+1,L(1)+L(2)+L(3)+L(4) +1
E,L(1)+L(2)+L(3)+L(4)+1,600+L(1)+L(2)+L(3)+L(4) +1
E,600+L(1)+L(2)+L(3)+L(4)+1,800+L(1)+L(2)+L(3) +L(4)+1
/SOL
ACEL,,9.81,0
!******************* ********************* ! 定义边界条件
!**************************************** D,400+1,UX,0,,,,UY,UZ
D,600+1,UX,0,,,,UY,UZ
D,21,UY,0
D,421,UY,0
D,621,UY,0
D,51,UY,0
D,651,UY,0
D,451,UY,0
D,81,UY,0
D,681,UY,0
D,481,UY,0
D,400+101,UX,0,,,,UY,UZ
D,600+101,UX,0,,,,UY,UZ
ANTYPE,0
/STATUS,SOLU
SOLVE
/POST1
/REP,FAST
A VPRIN,0, ,
ETABLE, ,SMISC, 2
A VPRIN,0, ,
ETABLE, ,SMISC,8 A VPRIN,0, , ETABLE, ,SMISC, 4 A VPRIN,0, , ETABLE, ,SMISC, 5 A VPRIN,0, , ETABLE, ,SMISC, 6 A VPRIN,0, , ETABLE, ,SMISC, 10 A VPRIN,0, , ETABLE, ,SMISC, 11 A VPRIN,0, , ETABLE, ,SMISC, 12 FINISH。