梁弯曲变形仿真的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命令流
!(1)工作环境设置/FILENAME,RC-BEAM !指定工作文件名/TITLE,ALAL YSIS OF A RC-BEAM !指定图形标题!(2)进行前处理器/prep7!(3)定义单元类型ET,1,LINK8 !定义钢筋单元ET,2,SOLID65 !定义混凝土单元ET,3,MESH200 !用于拉伸成体单元KEYOPT,3,1,6!(4)定义钢筋截面积r,1,28.3r,2,50.3r,3,314.1!(5)为solid65单元定义一个实参数组r,4,!(6)定义混凝土材料MP,EX,2,2.55E10MP,PRXY,2,0.3TB,CONC,2,1,9, !定义混凝土的破坏参数TBDA TA,,0.3,0.55,1.55E6,-1,,TBDA TA,,,,0.6!(7)定义钢筋材料模型及参数mp,ex,2,2e5 !纵向受拉钢筋材料mp,prxy,2,0.3tb,bkin,2,1,2,1tbdata,,350mp,ex,3,2e5 !横向箍筋，架立钢筋材料mp,prxy,3,0.25tb,bkin,3,1,2,1tbdata,,200!(8) 创建以及复制节点/pnum,node,1/pnum,elem,1n,1n,9,200fill,1,9ngen,11,9,1,9,1,,30ngen,11,99,1,99,1,,,-150/view,1,1,1,1!(9)建立箍筋单元type,1real,1mat,3!水平箍筋*do,i,11,16,1e,i,i+1e,i+(83-11),i+(83-11)+1*enddo!坚直箍筋*do,i,11,74,9e,i,i+9e,i+6,i+6+9*enddo!产生整个模型的箍筋egen,11,99,all!(10)建立架立筋以及纵筋单元!创建上部的架立钢筋单元*do,i,83,node(25,270,-1500+150),99 e,i,i+99e,i+6,i+6+99*enddo!纵向受拉钢筋单元的属性type,1real,3mat,2!创建纵筋单元*do,i,11,node(25,30,-1500+150),99 e,i,i+99e,i+3,i+3+99e,i+6,i+6+99*enddo/view,1,1,1,1/pnum,elem,0/pnum,node,0/eshape,1eplot!(11)建立混凝土剖面并划分网格k,1K,2,200,k,3,200,300k,4,,300a,1,2,3,4lsel,s,loc,y,0lsel,a,loc,y,300lesize,all,,,8lsel,alllsel,s,loc,x,0lsel,a,loc,x,200lesize,all,,,10type,3amesh,all!(12)拉伸形成混凝土单元type,2real,3mat,1extopt,esize,20extopt,aclear,1vext,all,,,,,-1500/pnum,mat,1/pnum,node,0/pnum,elem,0eplotallsel!(13)合并压缩节点编号nummrg,allnumcmp,alleplot!(14)施加支座约束nsel,s,loc,y,0nsel,r,loc,z,-1500+75d,all,uyd,all,ux!(15)施加对称面约束asel,s,loc,z,0da,all,symm!(16)施加载荷nsel,allnsel,s,loc,y,300nsel,r,loc,z,-450d,all,uy,-30allsel!(17)退出前处理器fini!(18) 进入求解器/solu!(19)求解器选项设置nlgeom,onnsubst,200outres,all,allneqit,50pred,oncnvtol,f,,0.05,2,0.5allsel!(20)求解并退出求解器solvefini!(21)进入通用处理器并读入最后一个子步的结果/post1set,last!(22)后处理操作plnsol,u,y,0,1/device,vector,1plcrack,0,0!(23)退出通用后处理器finish。
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计算命令流
!****************************************************************************** ***********************! 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在材料力学的弯曲变形教学中的运用作者：葛仁余来源：《科技视界》 2014年第15期葛仁余(安徽工程大学建筑工程学院，安徽芜湖 241000)【摘要】关于材料力学的弯曲变形问题，传统的教学方法就是积分法和叠加法。
本文在弯曲变形传统教学方法的基础上，设计了一种新的教学方法，即将传统教学方法和ANSYS有限元软件运用有机结合，将原本抽象的力学理论和枯燥乏味的概念直观化、形象化，使学生对弯曲变形获得感性认识，提高了学习兴趣。
【关键词】材料力学；弯曲变形；教学改革0 引言在材料力学的教学中，弯曲变形既是教学重点又是难点，当工程构件的受载情况不同时，获得的弯曲变形也就不同。
关于弯曲变形的转角和挠度的计算，在材料力学的教学中非常重要，教师在讲授时都做重点讲解，而学生在学习过程中，常常感到难以掌握，其中一个主要原因就是理论枯燥、难懂，且学生对弯曲变形感性认识欠缺。
为了攻克这一教学难点，笔者在材料力学的教学实践基础上，对弯曲变形的教学方法进行了研究，提出了ANSYS有限元软件应用与传统教学方法结合的一种新的教学方法。
实践证明，采用这种方法进行教学，有助于提高学生的形象思维能力，帮助学生理解教材内容，能激发学生学习材料力学的热情。
1 弯曲变形的传统教学方法材料力学教师在讲授梁的弯曲变形时，通常是先讲纯弯曲情况下，弯矩与曲率间的关系，再与高等数学中关于曲线的曲率公式联立，从而获得挠曲线微分方程近似表达式。
基于挠曲线近似微分方程，可由积分法获得转角公式和挠度公式，它们中的积分常数由边界条件获得，至此可以获得梁的弯曲变形完整的转角和挠度解析表达式，这种求解梁的弯曲变形的方法就叫积分法。
积分法求解单一荷载作用下梁的弯曲变形优势十分明显，但对于求多荷载作用下梁的弯曲变形问题，其计算过程十分繁杂，计算量也是较大。
所谓叠加法，就是在弯曲变形很小且材料服从胡克定律的情况下，挠曲线的近似常微分方程是线性的。
又因在小变形的前提下，计算弯矩时用梁变形前的位置，于是弯矩与载荷的关系也是线性的。
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=,C=27,ψ=45的命令如下：MP,EX,1,1E8MP,NUXY,1,TB,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命令得到更详细的解释;后继载荷步在后继载荷步中,用户可以随意杀死或重新激活单元;象上面提到的,要正确的施加和删除约束和结点载荷;用下列命令杀死单元：Command:EKILLGUI: Main Menu>Solution>-Load Step Opts-Other>Kill Elements用下列命令重新激活单元：Command: EALIVEGUI: Main Menu>Solution>-Load Step Opts-Other>Activate Elem第二个或后继载荷步：TIME,...ESEL,...EKILL,... 杀死选择的单元ESEL,...EALIVE,... 重新激活选择的单元...FDELE,... 删除不活动自由度的结点载荷D,... 约束不活动自由度...F,... 在活动自由度上施加合适的结点载荷DDELE,... 删除重新激活的自由度上的约束SAVESOLVE四、u /grid, keykey: “0”或“off”无网络“1”或“on” xy网络“2”或“x”只有x线“3”或“y”只有y线u xvar, nn: “0”或“1”将x轴作为时间轴“n”将x轴表示变量“n”“-1”u /axlab, axis, lab 定义轴线的标志axis: “x”或“y”lab: 标志,可长达30个字符u plvar, nvar, nvar2, ……,nvar10 画出要显示的变量作为纵坐标五、Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点缺省R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0”使用正负号“1”仅用绝对值六、VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体nv1:初始体号nv2:最终的体号ninc:体号之间的间隔kswp=0:只删除体kswp=1:删除体及组成关键点,线面如果nv1=all,则nv2,ninc不起作用七、VSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWPType,是选择的方式,有选择s,补选a,不选,全选all、反选inv等,其余方式不常用Item, Comp 是选取的原则以及下面的子项如 volu 就是根据实体编号选择,loc 就是根据坐标选取,它的comp就可以是实体的某方向坐标其余还有材料类型、实常数等MIN, VMAX, VINC,这个就不必说了吧,例：vsel,s,volu,,14vsel,a,volu,,17,23,2上面的命令选中了实体编号为 14,17,19,21,23的五个实体u rforce, nvar, node, item, comp, name 指定待存储的节点力数据nvar: 变量号node: 节点号item compF x,M x, y,zname: 给此变量一个名称,8个字符u add, ir, ia,ib,ic,name,--,--,facta, factb, factc 将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic：变量号name: 变量的名称Fini退出四大模块,回到BEGIN层/cle 清空内存,开始新的计算1．定义参数、数组,并赋值.2． /prep7进入前处理定义几何图形：关键点、线、面、体定义几个所关心的节点,以备后处理时调用节点号;设材料线弹性、非线性特性设置单元类型及相应KEYOPT设置实常数设置网格划分,划分网格根据需要耦合某些节点自由度定义单元表存盘3．/solu加边界条件设置求解选项定义载荷步求解载荷步4./post1通用后处理5./post26 时间历程后处理菜单命令7.参数化设计语言8.理论手册Fini退出四大模块,回到BEGIN层/cle 清空内存,开始新的计算1 定义参数、数组,并赋值.u dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组par: 数组名type： array 数组,如同fortran,下标最小号为1,可以多达三维缺省char 字符串组每个元素最多8个字符tableimax,jmax, kmax 各维的最大下标号var1,var2,var3 各维变量名,缺省为row,column,plane当type为table时2 /prep7进入前处理定义几何图形：关键点、线、面、体u csys,kcnkcn , 0 迪卡尔zuobiaosi1 柱坐标2 球4 工作平面5 柱坐标系以Y轴为轴心n 已定义的局部坐标系u numstr, label, value 设置以下项目编号的开始nodeelemkplineareavolu注意：vclear, aclear, lclear, kclear 将自动设置节点、单元开始号为最高号,这时如需要自定义起始号,重发numstru K, npt, x,y,z, 定义关键点Npt：关键点号,如果赋0,则分配给最小号u Kgen,itime,Np1,Np2,Ninc,Dx,Dy,Dz,kinc,noelem,imoveItime：拷贝份数Np1,Np2,Ninc：所选关键点Dx,Dy,Dz：偏移坐标Kinc：每份之间节点号增量noelem: “0”如果附有节点及单元,则一起拷贝;“1”不拷贝节点和单元imove：“0”生成拷贝“1”移动原关键点至新位置,并保持号码,此时itime,kinc,noelem被忽略注意：MAT,REAL,TYPE 将一起拷贝,不是当前的MAT,REAL,TYPEu A, P1, P2, ……… P18 由关键点生成面u AL, L1,L2, ……,L10 由线生成面面的法向由L1按右手法则决定,如果L1为负号,则反向;线需在某一平面内坐标值固定的面内u vsba, nv, na, sep0,keep1,keep2 用面分体u vdele, nv1, nv2, ninc, kswp 删除体kswp: 0 只删除体1 删除体及面、关键点非公用u vgen, itime, nv1, nv2, ninc, dx, dy, dz, kinc, noelem, imove 移动或拷贝体itime: 份数nv1, nv2, ninc：拷贝对象编号dx, dy, dz ：位移增量kinc: 对应关键点号增量noelem,：0：同时拷贝节点及单元1：不拷贝节点及单元imove： 0：拷贝体1：移动体u cm, cname, entity 定义组元,将几何元素分组形成组元cname: 由字母数字组成的组元名entity: 组元的类型volu, area, line, kp, elem, nodeu cmgrp, aname, cname1, ……,cname8 将组元分组形成组元集合aname: 组元集名称cname1……cname8: 已定义的组元或组元集名称u cmlist,nameu cmdele,nameu cmplot, label1定义几个所关心的节点,以备后处理时调用节点号;u n,node,x,y,z,thxy, thyz, thzx 根据坐标定义节点号如果已有此节点,则原节点被重新定义,一般为最大节点号;设材料线弹性、非线性特性u mp,lab, mat, co, c1,…….c4 定义材料号及特性lab: 待定义的特性项目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次项的系数u Tb, lab, mat, ntemp,npts,tbopt,eosopt 定义非线性材料特性表Lab: 材料特性表之种类Bkin: 双线性随动强化Bis 双线性等向强化Mkin: 多线性随动强化最多5个点Mis 多线性等向强化最多100个点Dp: dp模型Mat: 材料号Ntemp: 数据的温度数对于bkin: ntemp缺省为6mis ntemp缺省为1,最多20bis ntemp缺省为6,最多为6dp: ntemp, npts, tbopt 全用不上Npts: 对某一给定温度数据的点数u TBTEMP,temp,kmod 为材料表定义温度值temp: 温度值kmod: 缺省为定义一个新温度值如果是某一整数,则重新定义材料表中的温度值注意：此命令一发生,则后面的TBDATA和TBPT均指此温度,应该按升序若Kmod为crit, 且temp为空,则其后的tbdata数据为solid46,shell99,solid191中所述破坏准则如果kmod为strain,且temp为空,则其后tbdata数据为mkin中特性;u TBDATA, stloc, c1,c2,c3,c4,c5,c6给当前数据表定义数据配合tbtemp,及tb使用stloc: 所要输入数据在数据表中的初始位置,缺省为上一次的位置加1 每重新发生一次tb或tbtemp命令上一次位置重设为1,发生tb后第一次用空闲此项,则c1赋给第一个常数u tbpt, oper, x,y 在应力-应变曲线上定义一个点oper: defi 定义一个点dele 删除一个点x,y：坐标设置单元类型及相应KEYOPTu ET, itype, ename, kop1……kop6, inopr 设定当前单元类型Itype：单元号Ename：单元名设置实常数u Keyopt, itype, knum, valueitype: 已定义的单元类型号knum: 单元的关键字号value: 数值注意：如果 ,则必须使用keyopt命令,否则也可在ET命令中输入设置网格划分,划分网格映射网格划分1.面映射网格划分条件：a. 3或4条边b.面的对边必须划分为相同的单元或其划分与一个过渡形网格的划分相匹配c. 该面如有3条边,则划分的单元不必须为偶数,并且各边单元数相等d. mahkeye. mshpattern如果多于四条边,可将线合并成Lcomb可用amap命令,先选面,再选4个关键点即可指定面的对边的分割数,以生成过渡映射四边形网格,只适用于有四条边的面2. 体映射网格划分1若将体划分为六面体单元,必须满足以下条件a. 该体的外形为块状六面体、楔形或棱形五面体、四面体b. 对边必须划分为相同的单元数,或分割符合过渡网格形式c. 如果体是棱形或四面体,三角形面上的单元分割数必须是偶数2 当需要减少围成体的面数以进行映射网格划分时,可以对面相加或连接;如果连接而有边界线,线也必须连接在一起;3体扫掠生成网格步骤：a. 确定体的拓扑是否能够进行扫掠;侧面不能有孔；体内不能有封闭腔;源面与目标面必须相对b. 定义合适的单元类型c. 确定扫掠操作中如何控制生成单元层的数目 lesized. 确定体的哪一个边界面作为源面、目标面e. 有选择地对源面、目标面和边界面划分网格3. 关于连接线和面的一些说明连接仅是映射网格划分的辅助工具4. 用desize定义单元尺寸时单元划分应遵守的级别高：lesizekesizeesizedesize用smartzing定义单元尺寸时单元划分应遵守的级别高：lesizekesizesmartsizeu LESIZE,NL1,Size, Angsiz,ndiv,space,kforc,layer1,layer2,kyndiv 为线指定网格尺寸NL1: 线号,如果为all,则指定所有选中线的网格;Size: 单元边长,程序据size计算分割份数,自动取整到下一个整数Angsiz: 弧线时每单元跨过的度数Ndiv: 分割份数Space: “+”: 最后尺寸比最先尺寸“-“: 中间尺寸比两端尺寸free: 由其他项控制尺寸kforc 0: 仅设置未定义的线,1：设置所有选定线,2：仅改设置份数少的,3：仅改设置份数多的kyndiv: 0,No,off 表示不可改变指定尺寸1,yes,on 表示可改变u ESIZE,size,ndiv 指定线的缺省划分份数已直接定义的线,关键点网格划分设置不受影响u desize, minl, minh,……控制缺省的单元尺寸minl: n 每根线上低阶单元数缺省为3defa 缺省值stat 列出当前设置off 关闭缺省单元尺寸minh: n 每根线上高阶单元数缺省为2u mshape, key, dimension 指定单元形状key: 0 四边形2D,六面体3D1 三角形 2D, 四面体3DDimension: 2D 二维3D 三维u smart,off 关闭智能网格u mshkey, key 指定自由或映射网格方式key: 0 自由网格划分1 映射网格划分2 如果可能的话使用映射,否则自由即使自由smartsizing也不管用了u Amesh, nA1,nA2,ninc 划分面单元网格nA1,nA2,ninc 待划分的面号,nA1如果是All,则对所有选中面划分u SECTYPE, ID, TYPE, SUBTYPE, NAME, REFINEKEY定义一个截面号,并初步定义截面类型ID: 截面号TYPE: BEAM:定义此截面用于梁SUBTYPE: RECT 矩形CSOLID:圆形实心截面CTUBE: 圆管I: 工字形HREC: 矩形空管ASEC: 任意截面MESH: 用户定义的划分网格NAME: 8字符的截面名称字母和数字组成REFINEKEY: 网格细化程度：0~5对于薄壁构件用此控制,对于实心截面用SECDATA控制u SECDATA, VAL1, VAL2, …….VAL10 描述梁截面说明：对于SUBTYPE=MESH, 所需数据由SECWRITE产生,SECREAD读入u SECNUM,SECID 设定随后梁单元划分将要使用的截面编号u LATT, MAT, REAL, TYPE, --, KB, KE, SECNUM为准备划分的线定义一系列特性MAT: 材料号REAL: 实常数号TYPE: 线单元类型号KB、KE: 待划分线的定向关键点起始、终止号SECNUM: 截面类型号u SECPLOT,SECID,MESHKEY 画梁截面的几何形状及网格划分SECID:由SECTYPE命令分配的截面编号MESHKEY：0：不显示网格划分1：显示网格划分u /ESHAPE, SCALE 按看似固体化分的形式显示线、面单元SCALE: 0:简单显示线、面单元1：使用实常数显示单元形状u esurf, xnode, tlab, shape 在已存在的选中单元的自由表面覆盖产生单元xnode: 仅为产生surf151 或surf152单元时使用tlab: 仅用来生成接触元或目标元top 产生单元且法线方向与所覆盖的单元相同,仅对梁或壳有效,对实体单元无效Bottom产生单元且法线方向与所覆盖的单元相反,仅对梁或壳有效,对实体单元无效Reverse 将已产生单元反向Shape: 空与所覆盖单元形状相同Tri 产生三角形表面的目标元注意：选中的单元是由所选节点决定的,而不是选单元,如同将压力加在节点上而不是单元上u Nummrg,label,toler, Gtoler,actiontch 合并相同位置的itemlabel: 要合并的项目node: 节点, Elem,单元,kp: 关键点也合并线,面及点mat: 材料,type: 单元类型,Real: 实常数cp：耦合项,CE：约束项,CE: 约束方程,All：所有项toler: 公差Gtoler：实体公差Action: sele 仅选择不合并空合并注意：可以先选择一部分项目,再执行合并;如果多次发生合并命令,一定要先合并节点,再合并关键点;合并节点后,实体荷载不能转化到单元,此时可合并关键点解决问题;u Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneuunselectinve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元u Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点缺省R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0”使用正负号“1”仅用绝对值u NSLL,type, nkey 选择与所选线相联系的节点u nsla, type, nkey: 选择与选中面相关的节点type：s 选一套新节点r 从已选节点中再选a 附加一部分节点到已选节点u 从已选节点中去除一部分nkey: 0 仅选面内的节点1 选所有和面相联系的节点如面内线,关键点处的节点u esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元Type: S: 选择一组单元缺省R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组Stat: 显示当前选择状态Item： Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号u ALLSEL, LABT, ENTITY 选中所有项目LABT: ALL: 选所有项目及其低级项目BELOW: 选指定项目的直接下属及更低级项目ENTITY: ALL: 所有项目缺省VOLU:体高级AREA:面LINE :线KP:关键点ELEM:单元NODE:节点低级u Tshap,shape 定义接触目标面为2D、3D的简单图形Shape: line:直线Arc:顺时针弧Tria:3点三角形Quad:4点四边形………….根据需要耦合某些节点自由度u cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotznode1-node17: 待耦合的节点号;如果某一节点号为负,则此节点从该耦合组中删去;如果node1=all,则所有选中节点加入该耦合组;注意：1,不同自由度类型将生成不同编号2,不可将同一自由度用于多套耦合组u CPINTF, LAB, TOLER 将相邻节点的指定自由度定义为耦合自由度LAB：UX,UY,UZ,ROTX,ROTY,ROTZ,ALLTOLER: 公差,缺省为说明：先选中欲耦合节点,再执行此命令定义单元表说明：1,单元表仅对选中单元起作用,使用单元表之前务必选择一种类型的单元2,单元表各行为选中各单元,各列为每单元的不同数据u ETABLE, LAB, ITEM, COMP 定义单元表,添加、删除单元表某列LAB:用户指定的列名REFL, STAT, ERAS 为预定名称ITEM: 数据标志查各单元可输出项目COMP: 数据分量标志存盘u save, fname, ext,dir, slab 存盘fname : 文件名最多32个字符缺省为工作名ext: 扩展名最多32个字符缺省为dbdir: 目录名最多64个字符缺省为当前slab: “all”存所有信息“model”存模型信息“solv”存模型信息和求解信息3 /soluu /solu 进入求解器加边界条件u D, node, lab, value, value2, nend, ninc, lab2, lab3, ……lab6 定义节点位移约束Node : 预加位移约束的节点号,如果为all,则所有选中节点全加约束,此时忽略nend和ninc.Lab: ux,uy,uz,rotx,roty,rotz,allValue,value2: 自由度的数值缺省为0Nend, ninc: 节点范围为：node-nend,编号间隔为nincLab2-lab6: 将lab2-lab6以同样数值施加给所选节点;注意：在节点坐标系中讨论设置求解选项u antype, status, ldstep, substep, actionantype: static or 1 静力分析buckle or 2 屈曲分析modal or 3 模态分析trans or 4 瞬态分析status: new 重新分析缺省,以后各项将忽略rest 再分析,仅对static,full transion 有效ldstep: 指定从哪个荷载步开始继续分析,缺省为最大的,runn数指分析点的最后一步substep: 指定从哪个子步开始继续分析;缺省为本目录中,runn文件中最高的子步数action, continue: 继续分析指定的ldstep,substep说明：继续以前的分析因某种原因中断有两种类型singleframe restart: 从停止点继续需要文件：必须在初始求解后马上存盘单元矩阵或 .osav : 如果.esav坏了,将.osav改为.esavresults file: 不必要,但如果有,后继分析的结果也将很好地附加到它后面注意：如果初始分析生成了.rdb, .ldhi, 或rnnn 文件;必须删除再做后继分析步骤： 1进入anasys 以同样工作名2进入求解器,并恢复数据库3antype, rest4指定附加的荷载5指定是否使用现有的矩阵缺省重新生成kuse: 1 用现有矩阵6求解multiframe restart:从以有结果的任一步继续用不着u pred,sskey, --,lskey….. 在非线性分析中是否打开预测器sskey: off 不作预测当有旋转自由度时或使用solid65时缺省为offon 第一个子步后作预测除非有旋转自由度时或使用solid65时缺省为on -- ：未使用变量区lskey: off 跨越荷载步时不作预测缺省on 跨越荷载步时作预测此时sskey必须同时on注意：此命令的缺省值假定solcontrol为onu autots, key 是否使用自动时间步长key:on: 当solcontrol为on时缺省为onoff: 当solcontrol为off时缺省为off1: 由程序选择当solcontrol为on且不发生autots命令时在 .log文件中纪录“1”注意：当使用自动时间步长时,也会使用步长预测器和二分步长u NROPT, option,--,adptky 指定牛顿拉夫逊法求解的选项OPTION: AUT程序选择FULL:完全牛顿拉夫逊法MODI:修正的牛顿拉夫逊法INIT：使用初始刚阵UNSYM：完全牛顿拉夫逊法,且允许非对称刚阵ADPTKY:ON: 使用自适应下降因子OFF：不使用自适应下降因子u NLGEOM,KEYKEY: OFF:不包括几何非线性缺省ON：包括几何非线性u ncnv, kstop, dlim, itlim, etlim, cplim 终止分析选项kstop: 0 如果求解不收敛,也不终止分析1 如果求解不收敛,终止分析和程序缺省2如果求解不收敛,终止分析,但不终止程序dlim：最大位移限制,缺省为itlim: 累积迭代次数限制,缺省为无穷多etlim：程序执行时间秒限制,缺省为无穷cplim：cpu时间秒限制,缺省为无穷u solcontrol ,key1, key2,key3,vtol 指定是否使用一些非线性求解缺省值key1: on 激活一些优化缺省值缺省CNVTOL Toler=%Minref=对力和弯矩NEQIT 最大迭代次数根据模型设定在15~26之间ARCLEN 如用弧长法则用较更先进的方法PRED 除非有rotx,y,z或solid65,否则打开LNSRCH 当有接触时自动打开CUTCONTROL Plslimit=15%, npoint=13SSTIF 当NLGEOM,on时则打开NROPT,adaptkey 关闭除非：摩擦接触存在；单元12,26,48,49,52存在；当塑性存在且有单元20,23,24,60存在AUTOS 由程序选择off 不使用这些缺省值key2: on 检查接触状态此时key1为on此时时间步会以单元的接触状态据keyopt7的假定为基础当keyopt2=on 时,保证时间步足够小key3: 应力荷载刚化控制,尽量使用缺省值空：缺省,对某些单元包括应力荷载刚化,对某些不包括查nopl:对任何单元不包括应力刚化incp:对某些单元包括应力荷载刚化查vtol：u outres, item, freq, cname 规定写入数据库的求解信息item: all 所有求解项basic 只写nsol, rsol, nload, strsnsol 节点自由度rsol 节点作用荷载nload 节点荷载和输入的应变荷载strs 节点应力freq: 如果为n,则每n步包括最后一步写入一次none: 则在此荷载步中不写次项all: 每一步都写last: 只写最后一步静力或瞬态时为缺省定义载荷步u nsubst, nsbstp, nsbmx, nsbmn, carry 指定此荷载步的子步数nsbstp: 此荷载步的子步数如果自动时间步长使用autots,则此数定义第一子步的长度；如果solcontrol打开,且3D面-面接触单元使用,则缺省为1-20步；如果solcontrol打开,并无3D接触单元,则缺省为1子步；如果solcontrol关闭,则缺省为以前指定值；如以前未指定,则缺省为1nsbmx, nsbmn：最多,最少子步数如果自动时间步长打开u time, time 指定荷载步结束时间注意：第一步结束时间不可为“0”u f, node, lab, value, value2, nend, ninc 在指定节点加集中荷载node:节点号lab: Fx,Fy,Fz,Mx,My,Mzvalue: 力大小value2: 力的第二个大小如果有复数荷载nend,ninc：在从node到nend的节点增量为ninc上施加同样的力注意：1节点力在节点坐标系中定义,其正负与节点坐标轴正向一致u sfa, area, lkey, lab, value, value2 在指定面上加荷载area: n 面号all 所有选中号lkey: 如果是体的面,忽略此项lab: presvalue: 压力值u SFBEAM, ELEM, LKEY, LAB, VALI, VALJ, VAL2I, VAL2J, IOFFST, JOFFST 对梁单元施加线荷载ELEM: 单元号,可以为ALL,即选中单元LKEY: 面载类型号,见单元介绍;对于BEAM188,1为竖向；2为横向；3为切向VALI,VALJ: I, J节点处压力值VAL2I,VAL2J: 暂时无用IOFFST, JOFFST: 线载距离I, J 节点距离u lswrite, lsnum 将荷载与荷载选项写入荷载文件中lsnum ：荷载步文件名的后缀,即荷载步数当 stat 列示当前步数init 重设为“1”缺省为当前步数加“1”注意1. 尽量加面载,不加集中力,以免奇异点2. 面的切向荷载必须借助面单元求解载荷步u lssolve, lsmin, lsmax, lsinc 读入并求解多个荷载步lsmin, lsmax, lsinc ：荷载步文件范围4 /post1通用后处理u set, lstep, sbstep, fact, king, time, angle, nset 设定从结果文件读入的数据lstep :荷载步数sbstep：子步数,缺省为最后一步time：时间点如果弧长法则不用nset： data set numberu dscale, wn, dmult 显示变形比例wn: 窗口号或all,缺省为1dmult, 0或auto : 自动将最大变形图画为构件长的5%u pldisp, kund 显示变形的结构kund： 0 仅显示变形后的结构1 显示变形前和变形后的结构2 显示变形结构和未变形结构的边缘u get, par, node, n, u, xy,z 获得节点n的xy,z位移给参数par等价于函数 ux,uy,uznodex,y,z: 获得x,y,z节点号arnodex,y,z：获得和节点n相连的面注意：此命令也可用于/solu模块u fsum, lab, item 对单元之节点力和力矩求和lab: 空在整体迪卡尔坐标系下求和rsys 在当前激活的rsys坐标系下求和item: 空对所有选中单元不包括接触元求和cont: 仅对接触节点求和u PRSSOL, ITEM, COMP 打印BEAM188、BEAM189截面结果说明：只有刚计算完还未退出ANSYS时可用,重新进入ANSYS时不可用item comp 截面数据及分量标志S COMP X,XZ,YZ应力分量PRIN S1,S2,S3主应力SINT应力强度,SEQV等效应力EPTO COMP 总应变PRIN 总主应变,应变强度,等效应变EPPL COMP 塑性应变分量PRIN 主塑性应变,塑性应变强度,等效塑性应变u plnsol, item, comp, kund, fact 画节点结果为连续的轮廓线item: 项目见下表comp: 分量kund: 0 不显示未变形的结构1 变形和未变形重叠2 变形轮廓和未变形边缘fact: 对于接触的2D显示的比例系数,缺省为1item comp discriptionu x,y,z,sum 位移rot x,y,z,sum 转角s x,y,z,xy,yz,xz 应力分量1,2,3 主应力Int,eqv 应力intensity,等效应力epeo x,y,z,xy,yz,xz 总位移分量1,2,3 主应变Int,eqv 应变intensity,等效应变epel x,y,z,xy,yz,xz 弹性应变分量1,2,3 弹性主应变Int,eqv 弹性intensity,弹性等效应变eppl x,y,z,xy,yz,xz 塑性应变分量u PRNSOL, item, comp 打印选中节点结果item: 项目见上表comp: 分量u PRETAB, LAB1, LAB2, ……LAB9 沿线单元长度方向绘单元表数据LABn : 空：所有ETABLE命令指定的列名列名：任何ETABLE命令指定的列名u PLLS, LABI, LABJ, FACT, KUND 沿线单元长度方向绘单元表数据LABI:节点I的单元表列名LABJ:节点J的单元表列名FACT: 显示比例,缺省为1kund: 0 不显示未变形的结构1 变形和未变形重叠2 变形轮廓和未变形边缘5 /post26 时间历程后处理u nsol, nvar, node, item, comp,name在时间历程后处理器中定义节点变量的序号nvar：变量号从2到nv根据numvar定义node: 节点号item compu x, y,zrot x, y,zu ESOL, NVAR, ELEM, NODE, ITEM, COMP, NAME 将结果存入变量NVAR: 变量号,2以上ELEM: 单元号NODE: 该单元的节点号,决定存储该单元的哪个量,如果空,则给出平均值ITEM:COMP:NAME: 8字符的变量名, 缺省为ITEM加COMPu rforce, nvar, node, item, comp, name 指定待存储的节点力数据nvar: 变量号node: 节点号item compF x,M x, y,zname: 给此变量一个名称,8个字符u add, ir, ia,ib,ic,name,--,--,facta, factb, factc将ia,ib,ic变量相加赋给ir变量ir, ia,ib,ic：变量号name: 变量的名称u /grid, keykey: “0”或“off”无网络“1”或“on” xy网络“2”或“x”只有x线“3”或“y”只有y线u xvar, nn: “0”或“1”将x轴作为时间轴“n”将x轴表示变量“n”“-1”u /axlab, axis, lab 定义轴线的标志axis: “x”或“y”lab: 标志,可长达30个字符u plvar, nvar, nvar2, ……,nvar10 画出要显示的变量作为纵坐标u prvar, nvar1, ……,nvar6 列出要显示的变量6 PLOTCONTROL菜单命令u pbc, ilem, ……,key, min, max, abs 在显示屏上显示符号及数值item: u 所加的位移约束rot 所加的转角约束key: 0 不显示符号1 显示符号2 显示符号及数值u /SHOW, FNAME, EXT, VECT, NCPL 确定图形显示的设备及其他参数FNAME: X11:屏幕文件名：各图形将生成一系列图形文件JPEG: 各图形将生成一系列JPEG图形文件说明：没必要用此命令,需要的图形文件可计算后再输出7 参数化设计语言u do, par, ival, fval, inc 定义一个do循环的开始par: 循环控制变量ival, fval, inc：起始值,终值,步长正,负u enddo 定义一个do循环的结束u if,val1, oper, val2, base: 条件语句val1, val2: 待比较的值也可是字符,用引号括起来oper: 逻辑操作当实数比较时,误差为1e-10eq, ne, lt, gt, le, ge, ablt, abgtbase: 当oper结果为逻辑真时的行为lable: 用户定义的行标志stop: 将跳出anasysexit: 跳出当前的do循环cycle: 跳至当前do循环的末尾then: 构成if-then-else结构。
钢筋混凝土梁ansys分析附命令流
钢筋混凝⼟梁ansys分析附命令流钢筋混凝⼟⾮线性分析2015⼤作业上海交通⼤学陈明1、参数选择梁的截⾯宽度为200mm，上部配置2Φ8受压筋，混凝⼟的净保护层厚度为25 mm（从纵向钢筋外边缘算起），箍筋两端区采⽤8@100的双肢箍，中间区取8@200 双肢箍1）梁的截⾯⾼度选300mm；2）两加载间的距离选1000mm；3）混凝⼟选C30；4）纵向受拉钢筋配筋选218；2、描述选⽤的有限元模型及单元的特点采⽤ansys软件进⾏模拟计算，钢筋混凝⼟模型采⽤分离式模型，不考虑钢筋与混凝⼟之间的相对滑移。
混凝⼟采⽤solid65单元模拟，solid65⽤于模拟三维有钢筋或⽆钢筋的混凝⼟模型。
该单元能够计算拉裂和压碎。
在混凝⼟应⽤中，该单元的实体功能可以⽤于建⽴混凝⼟模型，同时，还可⽤加筋功能建⽴钢筋混凝⼟模型。
另外，该单元还可以应⽤于加强复合物和地质材料。
该单元由⼋个节点定义，每个节点有三个⾃由度:节点坐标系的x，y，z⽅向的平动。
⾄多可以定义三种不同规格的钢筋。
钢筋单元采⽤link180单元模拟，link180是⼀个适⽤于各类⼯程应⽤的三维杆单元。
根据具体情况，该单元可以被看作桁架单元、索单元、链杆单元或弹簧单元等等。
本单元是⼀个轴向拉伸⼀压缩单元，每个节点有三个⾃由度:节点坐标系的x，y，z⽅向的平动。
本单元是⼀种顶端铰接结构，不考虑单元弯曲。
本单元具有塑性、蠕变、旋转、⼤变形和⼤应变功能。
缺省时，当考虑⼤变形时任何分析中LINK180单元都包括应⼒刚化选项。
3、描述选⽤的混凝⼟与钢筋粘结滑移本构关系的具体形式、参数等。
钢筋的应⼒应变关系曲线考虑到极限塑性应变最⼤值为0.01,钢筋本构模型采⽤多线性模型kinh，初始弹性模量为Es=200000Mpa，强化系数为0.001。
混凝⼟的应⼒应变关系曲线混凝⼟选⽤各向同性的miso模型，当计⼊下降端时，程序报错，所以只取了前⾯的上升段，⽤5段折线模拟混凝⼟应⼒应变曲线。
ANSYS绘制弯矩、剪力、轴力图命令流完全教程
ANSYS绘制弯矩、剪力、轴力图命令流完全教程1.绘制弯矩图建立弯矩单元表。
例如梁单元i节点单元表名称为imom，j节点单元表名称为jmom，ETABLE,NI,SMISC,1 !单元I点轴力ETABLE,NJ,SMISC,7 !单元J点轴力ETABLE,QI,SMISC,2 !单元I点剪力ETABLE,QJ,SMISC,8 !单元J点剪力ETABLE,MI,SMISC,6 !单元I点弯矩ETABLE,MJ,SMISC,12 !单元J点弯矩plls,MI,MJ2.标注弯矩图PLOTCTRLS>>NUMBERING>>SVAL ON即可在画出弯矩图的同时在图上标出弯矩值的大小3.调整弯矩图如果弯矩图方向错误，则绘制弯矩图命令为plls,imom,jmom，-1同一个节点处两边的单元内力有细微差别，导致内力数字标注出现重影。
观察上面整体轴力图也可以发现，一段一段的，好像马赛克，其实上面整体弯矩图也是，不过不是很明显罢了。
这是EULER-BEONOULI梁理论以及ANSYS输出定义造成的（详细原因就不展开了，看看梁理论的书和ANSYS的说明吧）。
为了修正重影和节点两边内力值不一样的问题，遍制了宏文件ITFAVG.MAC命令文件内容如下：!---------------------------------------------------------------------!宏：ITFAVG.MAC（INTERNAL FORCE AVERAGE MACRO）!获取线性单元内力，并对单元边界处的内力进行平衡!输入信息!内力类型：MFORX,MFORY,MFORZ,MMOMX,MMOMY,MMOMZ*ASK,ITFTYPE,'PLEASE INPUT THE TYPE OF INTERNAL FORCE','MMOMY'!需处理的单元包*ASK,EASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF ELEMENTS TO BE PROCESSED!','EOUTER'!需处理的节点包*ASK,NASSEMBLY,'PLEASE INPUT THE COMPONENT NAME OF NODE TO BE PROCESSED!','NOUTER'!无需处理的节点包*ASK,UNASSEMBLY,'PLEASE INPUT THE COMPONENTNAME OF THE UNCHANGED NODE!(NONE IF THERE'S NO SUCH COMPONENT)','NONE'/POST1!输入信息：内力类型，欲处理单元的集合，欲处理节点的集合!ITFTYPE='MMOMY'!EASSEMBLY='EOUTER'!NASSEMBLY='NOUTER'!按内力类型确定ANSYS输出信息SMISC的编号*IF,ITFTYPE,EQ,'MFORX',THENITFINUM=1ITFJNUM=7*ELSEIF,ITFTYPE,EQ,'MFORY',THENITFINUM=2ITFJNUM=8*ELSEIF,ITFTYPE,EQ,'MFORZ',THENITFINUM=3ITFJNUM=9*ELSEIF,ITFTYPE,EQ,'MMOMX',THENITFINUM=4ITFJNUM=10*ELSEIF,ITFTYPE,EQ,'MMOMY',THENITFINUM=5ITFJNUM=11*ELSEIF,ITFTYPE,EQ,'MMOMZ',THENITFINUM=6ITFJNUM=12*ELSE*ENDIF!对不需平均的节点进行处理*IF,UNASSEMBLY,NE,'NONE',THEN!选出不进行处理的节点包并获取不进行处理节点的数目CMSEL,S,UNASSEMBLY*GET,UNNODNUM,NODE,0,COUNT!定义长度为UNNODNUM的数组（UNNOD），以存放选中单元的单元编号*DIM,UNNOD,ARRAY,UNNODNUM!将选中单元的编号按顺序存入数组UNNOD*DO,I,0,UNNODNUM-1,1UNNOD(I+1)=NDNEXT(I)*ENDDO*ELSEUNNODNUM=0*ENDIF!选出所需的单元和节点包CMSEL,S,EASSEMBLYCMSEL,S,NASSEMBLY!获得当前选中单元总数（存入变量SELELENUM）*GET,SELELENUM,ELEM,0,COUNT!定义长度为SELELENUM的数组（ELENUM），以存放选中单元的单元编号*DIM,ELENUM,ARRAY,SELELENUM!将选中单元的编号按顺序存入数组ELENUM*DO,I,0,SELELENUM-1,1ELENUM(I+1)=ELNEXT(I)*ENDDO!获得当前选中节点总数（存入变量SELNODNUM）*GET,SELNODNUM,NODE,0,COUNT!定义长度为SELNODNUM的数组（NODNUM），以存放选中单元的单元编号*DIM,NODNUM,ARRAY,SELNODNUM!将选中单元的编号按顺序存入数组NODNUM*DO,I,0,SELNODNUM-1,1NODNUM(I+1)=NDNEXT(I)*ENDDO!定义所需的线性单元内力ETABLE，节点I的内力存入数组ITNFI,!节点J的内力存入数组ITNFJETABLE,ITNFI,SMISC,ITFINUMETABLE,ITNFJ,SMISC,ITFJNUM!定义所需的结果数组，并将其置零ETABLE,ITNFINEO,SMISC,5SADD,ITNFINEO,ITNFI,,1ETABLE,ITNFJNEO,SMISC,11SADD,ITNFJNEO,ITNFJ,,1*DO,K,1,SELNODNUM,1!处理不需平均的节点INDEX=0*IF,UNNODNUM,GE,1,THEN*DO,J,1,UNNODNUM*IF,NODNUM(K),EQ,UNNOD(J),THENINDEX=1*ELSE*ENDIF*ENDDO*ELSE*ENDIF*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中，I节点（对线性单元而言）为节点K的单元编号*IF,NELEM(ELENUM(J),1),EQ,NODNUM(K),THENELEI=ELENUM(J)*EXIT*ELSE*ENDIF*ENDDO*DO,J,1,SELELENUM,1!选出和节点K相连的线性单元中，J节点（对线性单元而言）为节点K的单元编号*IF,NELEM(ELENUM(J),2),EQ,NODNUM(K),THENELEJ=ELENUM(J)*EXIT*ELSE*ENDIF*ENDDO*IF,INDEX,EQ,0,THEN*IF,ELEJ,NE,0,THEN !有可能出现ELEJ为0的情况!取出I节点为节点K的单元的I节点端的内力放入参数ETELEI *GET,ETELEI,ELEM,ELEI,ETAB,ITNFI!取出J节点为节点K的单元的J节点端的内力放入参数ETELEJ *GET,ETELEJ,ELEM,ELEJ,ETAB,ITNFJ!平均节点K的单元的I节点端的内力和节点K的单元的J节点端的内力ETAVE=(ETELEI+ETELEJ)/2!将平均后的内力存入结果数组中DETAB,ELEI,ITNFINEO,ETAVEDETAB,ELEJ,ITNFJNEO,ETAVE*ELSE*ENDIF*ELSE*ENDIF*ENDDO/UDOC,1,LOGO,OFFPLLS,ITNFINEO,ITNFJNEO!END OF ITFAVG.MAC（2）对体和面来说，ANSYS默认的结果输出格式是云图格式，而这种彩色云图打印为黑白图像时对比很不明显，无法表达清楚，对于发表文章非常不便。
ansys命令流
工程结构有限元数值仿真学院：力学与建筑工程学院专业： 14-硕结构工程学号：TSP140601031姓名：王卓第一次作业命令流：/BATCH/input,menust,tmp,'',,,,,,,,,,,,,,,,1 WPSTYLE,,,,,,,,0/PREP7ET,1,LINK180ET,2,LINK180 MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,2e5 MPDATA,PRXY,1,,0.3 ETDEL,2R,1, , ,RDEL,1R,1,400, , MPTEMP,,,,,,,,MPTEMP,1,0MPDE,EX,1MPDE,PRXY,1 MPDATA,EX,1,,2E+005 MPDATA,PRXY,1,,0.3 TYPE, 1MAT, 1REAL, 1ESYS, 0 SECNUM,TSHAP,LINEET,2,LINK180ETDEL,2K,1,,,,K,2,1500,,,K,3,,1500,,K,4,1500,1500,,K,5,250,2500,,K,6,1250,2500,,K,7,250,3500,,K,8,1250,3500,,K,9,-750,4000,,K,10,2250,4000,,K,11,-1750,4000,,K,11,-1750,4500,,K,12,-750,4500,,K,13,250,4500,,K,14,1250,4500,,K,15,2250,4500,,K,16,3250,4500,,K,16,3250,4500,, LSTR, 1, 2 LSTR, 2, 4 LSTR, 1, 3 LSTR, 3, 4 LSTR, 1, 4 LSTR, 4, 6 LSTR, 3, 5 LSTR, 5, 7 LSTR, 6, 8 LSTR, 7, 8 LSTR, 5, 6 LSTR, 4, 5 LSTR, 5, 8 LSTR, 8, 10 LSTR, 10, 16 LSTR, 16, 15 LSTR, 15, 14 LSTR, 14, 13 LSTR, 13, 12 LSTR, 12, 11 LSTR, 11, 9 LSTR, 9, 7 LSTR, 7, 13 LSTR, 8, 14 LSTR, 9, 13 LSTR, 9, 12 LSTR, 10, 14 LSTR, 10, 15 LSTR, 8, 13 TYPE, 1MAT, 1REAL, 1 ESYS, 0 SECNUM,LESIZE,ALL, , ,1, ,1, , ,1, FLST,2,28,4,ORDE,2 FITEM,2,2FITEM,2,-29LMESH,P51XFLST,2,26,4,ORDE,5 FITEM,2,1FITEM,2,-14FITEM,2,17FITEM,2,-27FITEM,2,29LMESH,P51XFINISH/SOLFLST,2,1,1,ORDE,1 FITEM,2,3/GOD,P51X, , , , , ,UX,UY,UZ, , , FLST,2,1,1,ORDE,1 FITEM,2,1/GOD,P51X, , , , , ,UY,UZ, , , , FLST,2,1,1,ORDE,1 FITEM,2,15/GOF,P51X,FY,-1000FLST,2,1,1,ORDE,1 FITEM,2,10/GOF,P51X,FY,-1000SOLVEFINISH/POST1A VPRIN,0, ,ETABLE, ,LS, 1A VPRIN,0, ,ETABLE, ,LS, 2A VPRIN,0, ,ETABLE, ,U,X轴力图：变形图：第二次作业命令流：/BATCH/input,menust,tmp,'',,,,,,,,,,,,,,,,1 WPSTYLE,,,,,,,,0/input,menust,tmp,'',,,,,,,,,,,,,,,,1 WPSTYLE,,,,,,,,0/NOPRKEYW,PR_SET,1KEYW,PR_STRUC,1KEYW,PR_THERM,0KEYW,PR_FLUID,0KEYW,PR_ELMAG,0KEYW,MAGNOD,0KEYW,MAGEDG,0KEYW,MAGHFE,0KEYW,MAGELC,0KEYW,PR_MULTI,0KEYW,PR_CFD,0/GO/PREP7ET,1,BEAM188ET,2,SHELL181/MPLIB,STATMPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,3e4MPDATA,PRXY,1,,0.2 MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,2.5e-9 MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,2,,2.1e5MPDATA,PRXY,2,,0.3 MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,2,,7.8e-9SECTYPE, 1, BEAM, RECT, ZHU, 0 SECOFFSET, CENTSECDATA,500,500,0,0,0,0,0,0,0,0,0,0 SECTYPE, 2, BEAM, RECT, BIANLIANG, 0 SECOFFSET, CENTSECDATA,300,500,0,0,0,0,0,0,0,0,0,0 SECTYPE, 3, BEAM, I, NEILIANG, 0 SECOFFSET, CENTSECDATA,200,200,500,16,16,10,0,0,0,0,0,0 SECTYPE, 3, BEAM, I, NEILIANG, 0 SECOFFSET, CENTSECDATA,200,200,500,16,16,10,0,0,0,0,0,0 /UI,BEAM,OFFsect,4,shell,,BANsecdata, 100,1,0.0,3secoffset,MIDseccontrol,,,, , , ,K,1,,,,FLST,3,1,3,ORDE,1FITEM,3,1KGEN,3,P51X, , , ,5400, , ,0FLST,3,3,3,ORDE,2FITEM,3,1FITEM,3,-3KGEN,4,P51X, , ,5400, , , ,0FLST,3,12,3,ORDE,2FITEM,3,1FITEM,3,-12KGEN,4,P51X, , , , ,5400, ,0LSTR, 1, 13LSTR, 13, 25LSTR, 25, 37LSTR, 2, 14LSTR, 14, 26LSTR, 26, 38LSTR, 3, 15LSTR, 15, 27LSTR, 27, 39LSTR, 4, 16LSTR, 16, 28LSTR, 28, 40LSTR, 5, 17LSTR, 17, 29LSTR, 29, 41LSTR, 6, 18LSTR, 18, 30LSTR, 30, 42LSTR, 7, 19LSTR, 19, 31LSTR, 8, 20 LSTR, 20, 32 LSTR, 32, 44 LSTR, 9, 21 LSTR, 21, 33 LSTR, 33, 45 LSTR, 10, 22 LSTR, 22, 34 LSTR, 34, 46 LSTR, 11, 23 LSTR, 23, 35 LSTR, 35, 47 LSTR, 12, 24 LSTR, 24, 36 LSTR, 36, 48 LSTR, 13, 14 LSTR, 14, 15 LSTR, 16, 17 LSTR, 17, 18 LSTR, 19, 20 LSTR, 20, 21 LSTR, 22, 23 LSTR, 23, 24 LSTR, 3, 6 LSTR, 6, 9 LSTR, 9, 12 LSTR, 2, 5 LSTR, 5, 8 LSTR, 8, 11 LSTR, 1, 4 LSTR, 4, 7 LSTR, 7, 10 GPLOTLSTR, 25, 26 LSTR, 26, 27 LSTR, 37, 38 LSTR, 38, 39 LSTR, 40, 41 LSTR, 41, 42 LSTR, 42, 33 LDELE, 60 LSTR, 13, 16 LSTR, 16, 19LSTR, 14, 17 LSTR, 17, 20 LSTR, 20, 23 LSTR, 15, 18 LSTR, 18, 21 LSTR, 21, 24 LSTR, 39, 42 LSTR, 42, 45 LSTR, 45, 48 LSTR, 27, 30 LSTR, 30, 33 LSTR, 33, 36 LSTR, 38, 41 LSTR, 41, 44 LSTR, 44, 47 LSTR, 37, 40 LSTR, 40, 43 LSTR, 43, 46 LSTR, 48, 47 LSTR, 47, 46 LSTR, 45, 44 LSTR, 44, 43 LSTR, 26, 29 LSTR, 29, 32 LSTR, 32, 35 LSTR, 36, 35 LSTR, 35, 34 LSTR, 30, 29 LSTR, 29, 28 LSTR, 33, 32 LSTR, 32, 31 LSTR, 25, 28 LSTR, 28, 31 LSTR, 31, 34 FLST,2,9,4,ORDE,2 FITEM,2,45FITEM,2,-53LDELE,P51XFLST,2,4,4FITEM,2,38FITEM,2,66FITEM,2,40FITEM,2,63AL,P51X FLST,2,4,4 FITEM,2,37 FITEM,2,63 FITEM,2,39 FITEM,2,60 AL,P51X FLST,2,4,4 FITEM,2,40 FITEM,2,67 FITEM,2,42 FITEM,2,64 AL,P51X FLST,2,4,4 FITEM,2,39 FITEM,2,64 FITEM,2,41 FITEM,2,61 AL,P51X FLST,2,4,4 FITEM,2,42 FITEM,2,68 FITEM,2,44 FITEM,2,65 AL,P51X FLST,2,4,4 FITEM,2,41 FITEM,2,65 FITEM,2,43 FITEM,2,62 AL,P51X FLST,2,4,4 FITEM,2,55 FITEM,2,72 FITEM,2,90 FITEM,2,85 AL,P51X FLST,2,4,4 FITEM,2,54 FITEM,2,85 FITEM,2,91 FITEM,2,94 AL,P51X FLST,2,4,4FITEM,2,73 FITEM,2,92 FITEM,2,86 AL,P51X FLST,2,4,4 FITEM,2,91 FITEM,2,86 FITEM,2,93 FITEM,2,95 AL,P51X FLST,2,4,4 FITEM,2,92 FITEM,2,74 FITEM,2,88 FITEM,2,87 AL,P51X FLST,2,4,4 FITEM,2,93 FITEM,2,87 FITEM,2,89 FITEM,2,96 AL,P51X FLST,2,4,4 FITEM,2,57 FITEM,2,69 FITEM,2,59 FITEM,2,75 AL,P51X FLST,2,4,4 FITEM,2,56 FITEM,2,75 FITEM,2,58 FITEM,2,78 AL,P51X FLST,2,4,4 FITEM,2,59 FITEM,2,70 FITEM,2,83 FITEM,2,76 AL,P51X FLST,2,4,4 FITEM,2,58 FITEM,2,76FITEM,2,79AL,P51XFLST,2,4,4 FITEM,2,83 FITEM,2,71 FITEM,2,81 FITEM,2,77AL,P51XFLST,2,4,4 FITEM,2,84 FITEM,2,77 FITEM,2,82 FITEM,2,80AL,P51X/UI,MESH,OFF FLST,5,18,4,ORDE,2 FITEM,5,1 FITEM,5,-18CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , , ,1 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 FLST,5,18,4,ORDE,2 FITEM,5,19 FITEM,5,-36CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , , ,1 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 TYPE, 1MAT, 1 REAL,ESYS, 0SECNUM, 1/UI,MESH,OFF FLST,5,4,4,ORDE,4 FITEM,5,56 FITEM,5,-57 FITEM,5,69 FITEM,5,78CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , , ,2 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 TYPE, 1MAT, 1 REAL,ESYS, 0 SECNUM, 1 CM,_Y,LINE LSEL, , , , 94 CM,_Y1,LINE CMSEL,S,_Y FLST,5,4,4,ORDE,4 FITEM,5,54 FITEM,5,-55 FITEM,5,72 FITEM,5,94CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , , ,2 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 FLST,5,10,4,ORDE,6 FITEM,5,56 FITEM,5,-57 FITEM,5,69 FITEM,5,-71FITEM,5,78 FITEM,5,-82CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , , ,2 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 FLST,5,10,4,ORDE,8 FITEM,5,54 FITEM,5,-55 FITEM,5,72 FITEM,5,-74 FITEM,5,88 FITEM,5,-89 FITEM,5,94 FITEM,5,-96CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , , ,2 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 FLST,5,10,4,ORDE,8 FITEM,5,37 FITEM,5,-38 FITEM,5,43 FITEM,5,-44 FITEM,5,60 FITEM,5,-62 FITEM,5,66 FITEM,5,-68CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , , ,2CMDELE,_Y CMDELE,_Y1 FLST,5,7,4,ORDE,6 FITEM,5,58 FITEM,5,-59 FITEM,5,75 FITEM,5,-77 FITEM,5,83 FITEM,5,-84 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,2, ,1, , , ,3 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 FLST,5,7,4,ORDE,4 FITEM,5,85 FITEM,5,-87 FITEM,5,90 FITEM,5,-93 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,2, ,1, , , ,3 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 FLST,5,7,4,ORDE,4 FITEM,5,39 FITEM,5,-42 FITEM,5,63 FITEM,5,-65 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,2, ,1, , , ,3CMDELE,_YCMDELE,_Y1FLST,5,18,5,ORDE,2FITEM,5,1FITEM,5,-18CM,_Y,AREAASEL, , , ,P51XCM,_Y1,AREACMSEL,S,_YCMSEL,S,_Y1AATT, 1, , 2, 0, 4 CMSEL,S,_YCMDELE,_YCMDELE,_Y1LESIZE,ALL, , ,20, ,1, , ,1,FLST,2,87,4,ORDE,4FITEM,2,1FITEM,2,-44FITEM,2,54FITEM,2,-96LMESH,P51XMSHKEY,0FLST,5,18,5,ORDE,2FITEM,5,1FITEM,5,-18CM,_Y,AREAASEL, , , ,P51XCM,_Y1,AREACHKMSH,'AREA'CMSEL,S,_YAMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2FINISH/SOLFLST,2,12,1,ORDE,12FITEM,2,1FITEM,2,62FITEM,2,123FITEM,2,184FITEM,2,245FITEM,2,306FITEM,2,367FITEM,2,428FITEM,2,489FITEM,2,550FITEM,2,611FITEM,2,672/GOD,P51X, , , , , ,UX,UY,UZ, , , ACEL,0,0,9.8,SOLVEFINISH/POST1A VPRIN,0, ,ETABLE, ,SMISC, 1A VPRIN,0, ,ETABLE, ,SMISC, 14A VPRIN,0, ,ETABLE, ,SMISC, 2A VPRIN,0, ,ETABLE, ,SMISC, 15 ANMODE,10,0.5, ,0 ANMODE,10,0.5, ,0 ANMODE,10,0.5, ,0 ANMODE,10,0.5, ,0结构图：柱轴力图：梁弯矩图：。
混凝土梁的ANSYS命令流
！ 水＊＊＊＊＊ ＊ ＊＊＊＊＊定义支座垫板的几何特性及材料特性＊＊＊＊＊＊ ＊ ＊
Ｅ ， Ｏ Ｉ４ Ｔ ３ Ｓ ＬＤ ５ ， Ｍ， ３３Ｏ １ ＰＥ， ． ０ Ｘ ， Ｅ
ＳＶ ＡＥ
！＊＊＊＊＊ 建立混凝土梁单元＊＊＊＊ ＊＊＊＊＊＊ ＊＊＊ ＊
附录 ２
Ｂ Ｏ Ｋ ０ Ｌ ０ ２ ， ． ＬＣ ， ０ ． ０ ０ ２ ， ， ５ ， ０ ＢＯ Ｋ ０ Ｌ０ ２ ， ５０ ０ ２ ＬＣ ， ， ５ ０ ３ ， ． ， ． ． ， ０ ＢＯ ＫＯ Ｌ ２ ， ５－ ． ０ ＬＣ ， ０ ５０ ３ ，０ １， ， ． ． ５ ＢＯ Ｋ ０ Ｌ０ ２ ， ５ ０２ ， ５ ＬＣ ， ， ５０ ３ ， ０ ０ ３ ， ． ． ． ． Ｖ Ｅ ， ． Ｓ ＬＡ Ｉ Ｌ ＶＬ ＥＡ Ｉ Ｇ Ｕ， Ｌ
ＳＶ ＡＥ
左翼缘构造筋
Ｌ 部构造钢筋划分网格 ［Ｅ ，， Ｅ，６ １， ＳｌＳ ＬＮ， １， １ Ｉ ６ ＬＥ， ＩＥ， ７１， Ｓ ＬＡ ＬＮ ，１， １ ， ７ ＬＴ ， ， ＡＴ ２ ３ ２ ， ＬＳ Ｚ ， ， ５ Ｅ ＩＥ Ａ Ｌ ００ Ｌ ．
ＬＥ ＩＡ Ｌ Ｍ ＳＩ Ｌ ， ｌ ｌＡＬ 名Ｅ ， Ｌ
Ｋ ２ ， ， ５ ０ １， ， ３ ５ ０ ０ ， ５ ４ ． ． ． Ｌ Ｔ ， ＳＲ ２ ， １ ２ ２ Ｌ Ｔ ， ＳＲ ２ ， ３ ２ ４
！＊＊＊＊＊下面纵向钢筋划分网格＊ ＊＊ ＊ ＊＊＊＊＊＊ ＊ ＊＊
ＬＥ ， ， Ｅ ， ３ １， Ｓ Ｌ Ｓ ＬＮ ， １， ２ Ｉ ５
Ｒ１ ，，
！ 定义混凝土的常数
Ｍ， ， ｅ０弹性模量要根据其他公式求 ＰＥ， ３ ｌ！ Ｘ１２ ４ ．
Ｍ ， Ｙ １０ ２ Ｐ Ｐ Ｘ ， ． Ｒ ， ０ Ｍ ， Ｓ ｌ２Ｅ 密度 ＰＤＮ ， ５３！ Ｅ ，
ansys命令流中文说明
ansys命令流中文说明展开全文KB、KE: 待划分线的定向关键点起始、终止号SECNUM: 截面类型号u SECPLOT,SECID,MESHKEY 画梁截面的几何形状及网格划分SECID:由SECTYPE命令分配的截面编号MESHKEY：0：不显示网格划分1：显示网格划分u /ESHAPE, SCALE 按看似固体化分的形式显示线、面单元SCALE: 0:简单显示线、面单元1：使用实常数显示单元形状u esurf, xnode, tlab, shape 在已存在的选中单元的自由表面覆盖产生单元xnode: 仅为产生surf151 或surf152单元时使用tlab: 仅用来生成接触元或目标元top 产生单元且法线方向与所覆盖的单元相同，仅对梁或壳有效，对实体单元无效Bottom产生单元且法线方向与所覆盖的单元相反，仅对梁或壳有效，对实体单元无效Reverse 将已产生单元反向Shape: 空与所覆盖单元形状相同Tri 产生三角形表面的目标元注意：选中的单元是由所选节点决定的，而不是选单元，如同将压力加在节点上而不是单元上u Nummrg,label,toler, Gtoler,action,switch 合并相同位置的itemlabel: 要合并的项目node: 节点， Elem,单元，kp: 关键点（也合并线，面及点）mat: 材料，type: 单元类型，Real: 实常数cp：耦合项，CE：约束项，CE: 约束方程，All：所有项toler: 公差Gtoler：实体公差Action: sele 仅选择不合并空合并switch: 较低号还是较高号被保留（low, high）注意：可以先选择一部分项目，再执行合并。
如果多次发生合并命令，一定要先合并节点，再合并关键点。
合并节点后，实体荷载不能转化到单元，此时可合并关键点解决问题。
u Lsel, type, item, comp, vmin, vmax, vinc, kswp 选择线type: s 从全部线中选一组线r 从当前选中线中选一组线a 再选一部线附加给当前选中组aunoneu(unselect)inve: 反向选择item: line 线号loc 坐标length 线长comp: x,y,zkswp: 0 只选线1 选择线及相关关键点、节点和单元u Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备Type: S: 选择一组新节点（缺省）R: 在当前组中再选择A: 再选一组附加于当前组U: 在当前组中不选一部分All: 恢复为选中所有None: 全不选Inve: 反向选择Stat: 显示当前选择状态Item: loc: 坐标node: 节点号Comp: 分量Vmin,vmax,vinc: ITEM范围Kabs: “0” 使用正负号“1”仅用绝对值u NSLL,type, nkey 选择与所选线相联系的节点u nsla, type, nkey: 选择与选中面相关的节点type：s 选一套新节点r 从已选节点中再选a 附加一部分节点到已选节点u 从已选节点中去除一部分nkey: 0 仅选面内的节点1 选所有和面相联系的节点（如面内线，关键点处的节点）u esel, type, item, comp, vmin, vmax, vinc, kabs 选择一组单元Type: S: 选择一组单元（缺省）R: 在当前组中再选一部分作为一组A: 为当前组附加单元U: 在当前组中不选一部分单元All: 选所有单元None: 全不选Inve: 反向选择当前组（？）Stat: 显示当前选择状态Item： Elem: 单元号Type: 单元类型号Mat: 材料号Real: 实常数号Esys: 单元坐标系号u ALLSEL, LABT, ENTITY 选中所有项目LABT: ALL: 选所有项目及其低级项目BELOW: 选指定项目的直接下属及更低级项目ENTITY: ALL: 所有项目（缺省）VOLU:体高级AREA:面LINE :线KP:关键点ELEM:单元NODE:节点低级u Tshap,shape 定义接触目标面为2D、3D的简单图形Shape: line:直线Arc:顺时针弧Tria:3点三角形Quad:4点四边形………….2.6 根据需要耦合某些节点自由度u cp, nset, lab,,node1,node2,……node17nset: 耦合组编号lab: ux,uy,uz,rotx,roty,rotznode1-node17: 待耦合的节点号。
ANSYS分析实例与工程应用命令流学习笔记
ANSYS分析实例与工程应用命令流学习笔记1大纲静力分析：2杆、3梁、5薄膜和板壳、4实体单元梁单元：简化计算，结构总体受力情况实体单元：较复杂的结构，局部细节的受力情况稳定性分析：6振动、模态分析：7简单振动和梁的振动、8膜板和实体振动2杆系结构的静力分析2.1铰接杆在外力作用下的变形二维杆单元LINK1*AFUN,DEG：三角函数默认为弧度，改为角度后处理：结构变形图、显示节点位移和杆件应力2.2人字形屋架的静力分析后处理：杆单元的轴力、轴向应力、轴向应变2.3超静定拉压杆的反力计算后处理：节点反力2.4平行杆件与刚性梁连接的热应力问题定义3点的UY为耦合自由度，即三者的UY位移相等温度（增量）后处理：寻找特定位置的节点和单元，并从单元表中提取它们的内力2.5端部有间隙的杆的热膨胀二维带厚度的平面应力单元PLANE42、二维接触单元CONTACT26温度（始、末）后处理：定义水平应力和铅直应力单元表，并提取3号单元的应力结果*Status,ParmFINISH定义数组变量，将计算结果通过数组变量输出到文件3梁的弯曲静力分析3.1单跨等截面超静定梁的平面弯曲二维弹性梁单元BEAM3后处理：定义以两端弯矩和剪力的单元表，并列出单元表数据并用单元表数据绘制剪力图和弯矩图更细的节点划分方案，更精细3.2四跨连续梁的内力计算体素建模：keypoint, line, area, volume便于细分单元3.3七层框架结构计算3.4工字形截面外伸梁的平面弯曲3.5矩形截面梁的纵横弯曲分析考虑应力强化效应后处理：迭代过程3.6空间刚架静力分析三维梁单元BEAM43.7悬臂梁的双向弯曲三维8节点耦合场实体单元SOLID5三维20节点固体单元SOLID92三维10节点耦合场实体单元SOLID98三维结构实体自适应单元SOLID147定义宏程序，对应四种工况，各种结果差别不大3.8圆形截面悬臂杆的弯扭组合变形三维直管单元PIPE16（只定义外直径，不定义内直径）3.9悬臂等强度梁的弯曲四边形壳单元SHELL63（这里用退化的三角形单元，并使用节点耦合自由度保证模型的对称变形）三维非对称锥形梁单元BEAM44（定义横截面主轴，单元宽度线性变化）计算结果都很好，但壳体单元更能模拟出等强度梁的实际几何形状，更直观，截面定义更简单。
ansys命令流实例
ansys命令流实例1.结构静力分析－－梁分析/PREP7 !进入前处理器ET,1,BEAM3 !定义单元类型R,1,0.25,0.0052,0.5 !定义实常数MP,EX,1,210E6MP,PRXY,1,0.3 !定义材料属性N,1,0N,2,1N,3,5N,4,7N,5,7N,6,9N,7,11N,8,11N,9,13N,10,14 !定义节点E,1,2E,2,3E,3,4E,5,6E,6,7E,8,9E,9,10 !生成单元CP,1,UX,4,5CP,2,UY,4,5CP,3,UX,7,8CP,4,UY,7,8 !耦合节点FINISH/SOLU !进入求解器D,2,UXD,2,UYD,3,UYD,6,UYD,9,UY !施加位移约束F,10,FY,-4 !施加集中约束SFBEAM,1,1,PRES,4,4SFBEAM,2,1,PRES,4,4 !施加均布力SOLVE !求解/POST1 !进入后处理器PLDISP !绘制结构变形图PRDISP !列出各节点的位移ETABLE,IMOMENT,SMISC,6ETABLE,JMOMENT,SMISC,12ETABLE,ISHEAR,SMISC,2ETABLE,JSHEAR,SMISC,8 !将节点弯矩、剪力制表PRETAB !列表显示单元的弯矩、剪力/TITLE,SHEAR FORCE DISTRIBUTION !设置剪力分布图的标题PLLS,ISHEAR,JSHEAR !绘制剪力分布图/TITLE,BENDING MOMENT IDSTRIBUTION !设置弯矩分布图的标题PLLS,IMOMENT,JMOMENT !绘制弯矩分布图流固耦合实例RAD=0.8 !底面半径H=1G=9.8OMEGAR=2ROU=1000 !定义参数变量/PREP7 !进入前处理器ET,1,FLUID79 !选择单元类型KEYOPT,1,3,1 !设置单元关键字MP,EX,1,2E9 !设置杨氏模量MP,DENS,1,ROU !设置材料密度K,1K,2,RADK,3,RAD,HK,4,,H !生成关键点A,1,2,3,4 !连接关键点生成面积LESIZE,ALL,,,10 !设置网格划分精度AMESH,ALL !将面积划分网络/SOLU !进入求解器DL,2,,UXDL,1,,UYNSEL,S,LOC,XDSYM,SYMM,XD,ALL,UXNSEL,ALL !施加位移约束ACEL,,GOMEGA,,OMEGAR !施加惯性力SOLVE !求解/POST1SET,LAST !进入通用后处理器PLNSOL,U,X,0,1 !绘制应力云图UCENT=UY(22)UEDGE=UY(12)UELEV=UEDGE-UCENT !提取节点位移结构静力分析－－壳结构内力分析LENGTH=100YOUNG=200000THICKNESS=2FORCE=1000DENSITY=9E-6 !将材质、载荷、板的几何尺寸等参数化/PREP7 !进入前处理器MP,EX,1,YOUNGMP,NUXY,1,0.3MP,DENS,1,DENSITY !定义材质ET,1,SHELL63 !定义单元类型R,1,THICKNESS,THICKNESS,THICKNESS,THICKNESS !定义实常数!构建结构的几何模型K,1,0,0K,2,LENGTH,0K,3,LENGTH,LENGTHK,4,0,LENGTH !定义关键点A,1,2,3,4!通过关键点生成面LSEL,ALLLESIZE,ALL,,,16AMESH,ALL !设定网格划分参数，划分网格FINISH/SOLU !进入求解器NSEL,S,LOC,X,0,0D,ALL,ALL,0 !选择X＝0的节点将其固定NSEL,S,LOC,X,LENGTH,LENGTHD,ALL,ALL,0 !选择X＝LENGTH的节点将其固定NSEL,S,LOC,X,0.5*LENGTH,0.5*LENGTHNSEL,R,LOC,Y,0.5*LENGTH,0.5*LENGTHF,ALL,FZ,FORCEALLSEL !捕捉板的中心点并在中心点处施加集中力荷载SOLVE !求解FINISH/POST1 !进入后处理器/DSC,,10PLNSOL,U,Z,0,1 !绘图显示板的竖向变形NSEL,ALL !提取板的最大竖向变形NSORT,U,Z,1,1 !将节点的位移绝对值以升序排序*GET,MAXDEFLECTION,SORT,0,MAX !提取位移最大值并赋给变量NSEL,S,LOC,X,0NSEL,A,LOC,X,LENGTH,LENGTH !选择固定边节点NSORT,S,EQV,1,1 !将节点等效应力的绝对值以升序排序*GET,MAXSTRESS,SORT,0,MAX !提取等效应力最大值并赋给变量结构静力分析－－平面桁架分析/PREP7 !进入前处理器ET,1,LINK1 !选择单元R,1,0.1 !定义实常数MP,EX,1,30E6MP,PRXY,1,0.3 !定义材料属性N,1,0N,2,4N,3,8N,4,12N,5,0,3N,6,4,3N,7,8,3 !生成节点E,1,2E,2,3E,3,4E,4,7E,3,7E,2,7E,2,6E,2,5E,1,5E,5,6E,6,7 !生成单元FINISH/SOLU !进入求解器D,1,UXD,1,UYD,5,UX !实加位移约束F,2,FY,-15F,3,FY,-15F,4,FY,-15 !施加集中力SOLVE !求解/POST1 !进入能用后处理器PRESOL,FORC !列表显示反力。
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命令流使用方法(中文)修改
Finish(退出四大模块，回到BEGIN层)/clear (清空内存，开始新的计算)1．定义参数、数组，并赋值.2．/prep7(进入前处理)定义几何图形：关键点、线、面、体定义几个所关心的节点，以备后处理时调用节点号。
设材料线弹性、非线性特性设置单元类型及相应KEYOPT设置实常数设置网格划分，划分网格根据需要耦合某些节点自由度定义单元表3．/solu加边界条件设置求解选项定义载荷步求解载荷步4./post1（通用后处理）5./post26 （时间历程后处理）6.PLOTCONTROL菜单命令7.参数化设计语言8.理论手册Finish(退出四大模块，回到BEGIN层)/clear (清空内存，开始新的计算)1.定义参数、数组，并赋值.dim, par, type, imax, jmax, kmax, var1, vae2, var3 定义数组par: 数组名type：array 数组，如同fortran,下标最小号为1，可以多达三维（缺省）char 字符串组（每个元素最多8个字符）tableimax,jmax, kmax 各维的最大下标号var1,var2,var3 各维变量名，缺省为row,column,plane(当type为table时) 2./prep7(进入前处理)2.1 设置单元类型及相应KEYOPTET, itype, ename, kop1……kop6, inopr 设定当前单元类型Itype：单元号Ename：单元名设置实常数Keyopt, itype, knum, valueitype: 已定义的单元类型号knum: 单元的关键字号value: 数值注意：如果,则必须使用keyopt命令，否则也可在ET命令中输入2.2 定义几个所关心的节点，以备后处理时调用节点号。
n,node,x,y,z,thxy, thyz, thzx 根据坐标定义节点号如果已有此节点，则原节点被重新定义，一般为最大节点号。
ANSYS简支梁计算命令流
ANSYS简支梁计算命令流finish/clear/PREP7et,1,link8 !定义link8单元et,2,solid65 !定义solid65单元keyopt,2,7,1r,1,314 !定义link8单元的面积r,2 !定义solid65的实常数号mp,ex,1,2e5 !定义link8单元的弹性模量mp,prxy,1,0.3 !定义link8单元的泊松系数mp,ex,2,4e4 !定义solid65单元的弹性模量mp,prxy,2,0.3 !定义solid65单元的泊松系数tb,concr,2 !定义2号为混凝土tbdata,,0.9,1,1.8,50 !定义混凝土的c1,c2,Rl,Ra!----------建立几何模型-------------------------blc4, , ,50,200,1500 !定义梁体/view,1,1,1,1 !定义ISO查看/ang,1vplot !绘制梁体kwpave,5 !工作平面移动到关键点5wpoff,25 !工作平面移动25mmwprot,0,0,90 !工作平面旋转vsbw,1 !分割梁体wpoff,0,40 !工作平面移动40mmwprot,0,90 !工作平面旋转vsbw,all !分割梁体wpoff,300 !再将梁体分割为三个区域wprot,0,0,90 !(原为控制单元密度而设)vsbw,allwpoff,0,0,900vsbw,allwpstyl !关闭工作平面显示nummrg,all !整理编号numcmp,all!------------划分单元网格---------------------------lsel,s,,,41,57,16 !定义line41,line57为新的选择集lsel,a,,,9 !定义line9也在选择集中latt,1,1,1 !定义选择集的属性lesize,9,,,20,,,,,1 !将line9划分为20个单元lesize,41,,,20,,,,,1 !将line41划分为20个单元lesize,57,,,60,,,,,1 !将line57划分为60个单元lmesh,all !对线划分单元!-------------allsel,all !新的选择集为所有的实体vsel,s,,,all !定义所有体为选择集vatt,2,2,2 !定义选择集的属性lesize,22,,,2,,,,,0 !定义体线的单元数以控制体单元lesize,23,,,2,,,,,0lesize,33,,,2,,,,,0lesize,28,,,4,,,,,0mshape,0,3d !将体划分单元的形状定位HEX mshkey,1 !采用MAPPED划分器vmesh,all !对体进行划分单元finish!--------------定义约束及荷载-------------------------/soluasel,s,,,2,22,20asel,a,,,14,17,3da,all,uy !对支撑端面施加Uy方向约束asel,s,,,3,23,20asel,a,,,15,19,4da,all,uz !对正对称面施加Uz方向约束asel,s,,,32,36,4asel,a,,,50,52,2asel,a,,,41,43,2da,all,ux !对侧对称面施加Ux方向约束asel,s,,,21,27,6asel,a,,,42,46,4asel,a,,,35,51,16sfa,all,1,pres,0.04 !对顶面施加面荷载(N/mm*mm) allsel,all!--------------定义非线性分析选项-------------------- cnvtol,f,,0.05,2 !定义收敛条件,使用缺省的VALUE nsubst,25 !定义子步数outres,all,all !输出每一子步的结果autots,1 !打开自动时间步控制lnsrch,1 !打开线性搜索ncnv,2 !如果不收敛时结束而不退出neqit,50 !每一子步中方程的迭代次数限值pred,on !打开预测器time,50 !定义时间solve !求解!--------------进入后处理过程----------------------- finish/post1pldisp,0 !绘制变形图(缺省为最后一步的) etable,sig,ls,1 !定义应力表plls,sigi,sigi,1 !绘制上述应力plnsol,s,z,0,1 !绘制Z方向的应力。
【干货】基于ANSYS的悬臂梁模态分析
【干货】基于ANSYS的悬臂梁模态分析1、连续系统的振动实际的振动系统都是连续体，它们具有连续分布的质量与弹性，因而又称连续系统或分布参数系统。
由于确定连续体上无数质点的位置需要无限多个坐标，因此连续体是具有无限多自由度的系统。
连续体的振动要用时间和空间坐标的函数来描述，其运动方程不再像有限多自由度系统那样是二阶常微分方程组，它是偏微分方程。
在物理本质上，连续体系统和多自由度系统没有什么差别，连续体振动的基本概念与分析方法与有限多自由度系统是完全类似的。
2、说明(1) 本章讨论的连续体都假定为线性弹性体，即在弹性范围内服从虎克定律。
(2) 材料均匀连续；各向同性。
(3) 振动满足微振动的前提。
3、梁的弯曲振动动力学方程考虑细长梁的横向弯曲振动梁参数：ρ单位体积梁的质量E弹性模量I截面对中性轴的惯性距S 梁横截面积外部力：m(x,t): 单位长度梁上分布的外力矩f(x,t): 单位长度梁上分布的外力假设：(1) 梁各截面的中心惯性轴在同一平面xoy内(2) 外载荷作用在该平面内(3) 梁在该平面作横向振动（微振）(4) 这时梁的主要变形是弯曲变形(5) 在低频振动时可以忽略剪切变形以及截面绕中性轴转动惯量的影响伯努利－欧拉梁（Bernoulli-Euler Beam）令：y(x,t):距原点x处的截面在t时刻的横向位移微段受力分析力平衡方程：4、悬臂梁的固有频率和模态函数5、两端固定杆的纵向模态分析问题描述：一悬臂梁截面为矩形，如图1所示，几何尺寸及材料特性如下，分析其前三阶固有频率及振型。
GUI操作如下：一、菜单建模分析过程第一步，清除内存准备分析1) 清除内存：选择菜单Utility Menu>File>Clear& Start New，单击OK按钮。
2) 更换工作文件名：选择菜单Utility Menu>File>ChangeJobname，输入vibration of cantilever，单击OK按钮。
实验三简支梁的变形Ansys分析
实验三简支梁的变形分析实验目的：了解和掌握简支梁的变形分析的方法和步骤。
实验内容：完成工字梁端面受力分析。
图1 是一工字梁受力简图，工字梁型号为32a，其跨度为1.0m，求其在集中力P1和P2作用下O点的挠度。
该工字梁材料的弹性模量为220GPa，泊松比为0.3，集中力P1=8000N，P2=5000N。
要求采用自顶向下的建模方法即采用3D体素建立起几何模型。
图1 工字梁受力简图该问题属于梁的弯曲问题。
由于该工字梁属于标准型材，在分析过程中选择图2所示的几何模型，通过选择相应的梁单元进行求解，这样就简化了3D模型复杂的建模过程和求解步骤。
图2 几何模型实验步骤：⒈定义工作文件名和工作标题1）选择Utility Menu→File→Change Jobname 命令，出现Change Jobname对话框，在[/FILNAM] Enter new jobname 输入栏中输入工作文件名exercise3，并将New log and error files 设置为Yes,单击【OK】按钮关闭该对话框。
2）选择Utility Menu→File→Change Title命令，出现Change Title对话框，在[/TITLE]Enter new title 栏中输入I BEAM SUBJECTED TO CONCENTRATED FORCE，单击【OK】按钮关闭该对话框。
⒉定义单元类型1）选择Main Menu→Preprocessor→ElementType→Add/Edit/Delete命令,出现ElementTypes 对话框，单击【Add】按钮, 出现Library of ElementTypes对话框。
2）在Library of ElementTypes复选框中选择Structural Beam 2D elastic 3,在Element type reference number 输入栏中输入1，如图3所示，单击【OK】按钮关闭该对话框。