(1)转贴一个例题,来自中国有限元联盟论坛-运行没有问题/title,fixed-fixedconcretebeamexample/prep7et,1,65mp,ex,1,3e7!steelrebar(unitsarepounds,inches)MP,PRXY,1,0.3mp,ex,2,1e6,!concreteMP,PRXY,2,0.3mp,dens,2,.00025tb,concr,2tbdata,1,.3,.5,200,4000!shearcoeffs,tensileandcompressstrengthr,1,1,.03,0,0!mat1(steel),3percentreinforcementinxdirr,2,1,.01,0,0!mat1(steel),1percentr,3,1,.04,0,0!mat1(steel),4percentblock,,100,,5,,5block,,100,5,10,,5block,,100,10,15,,5!vovlap,allNUMMRG,KP,,,,LOWnumcmp,voluesize,5mat,2!concretematerialreal,1!rebarvmesh,1real,2vmesh,2real,3vmesh,3nsel,s,loc,xd,all,allnsel,s,loc,x,100d,all,allnsel,allfini/solunsel,s,loc,y,15sf,all,pres,100nsel,allOUTRES,ALL,ALL,nsub,10solvefini/post1/DEVICE,VECTOR,1set,lastplcrackfini(2)!ANSYSSOLID65环向布置钢筋的例子――运行没有问题!一个管道,环向配筋率为1%,纵向配筋率为0.5%,径向配筋率为0.1%!作者:陆新征清华大学土木系!FINISH/CLEAR/PREP7!*!单元属性ET,1,SOLID65!*KEYOPT,1,1,0!表示考虑大变形KEYOPT,1,5,0!表示只打印质心的线性解KEYOPT,1,6,0!表示只打印质心的线性解KEYOPT,1,7,1!表示开裂后考虑应力松弛,有助于计算收敛!*!实参数1:不同方向配筋R,1,2,.001,,,2,.01,RMORE,90,,2,.005,90,90,!材料属性!混凝土基本材料属性MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,30E2MPDATA,PRXY,1,,.2!屈服准则TB,MISO,1,1,5,!多线性等向强化,材料号为1,数据的温度数为1,对某一给定温度数据的点数为5。TBTEMP,0!为材料表定义温度值。TBPT,,0.0005,15TBPT,,0.001,21TBPT,,0.0015,24TBPT,,0.002,27TBPT,,0.003,24!破坏准则TB,CONC,1,1,9,!材料号为1,数据的温度数为1,对某一给定温度数据的点数为9。TBTEMP,0!为材料表定义温度值。TBDATA,,.5,.9,3,30,,!给当前数据表定义数据TBDATA,,,,1,,,MPTEMP,,,,,,,,MPTEMP,1,0!钢材基本属性MPDATA,EX,2,,200E3MPDATA,PRXY,2,,.27!屈服准则TB,BISO,2,1,2,!双线性等向强化,材料号为2,数据的温度数为1,对某一给定温度数据的点数为2。TBTEMP,0!为材料表定义温度值为0度。TBDATA,,310,2E3,,,,!给当前数据表定义数据。!管道内径和外径CYL4,0,0,3000,,,,10000!生成内径为2000,外径为3000,长度为10000的管道CYL4,0,0,2000,,,,10000VSBV,1,2!体相减,在体1中减去体2!定义局部柱坐标CSWPLA,11,1,1,1,!通过当前定义的工作平面的原点为中心建立局部坐标系,第一个1表示为柱坐标。后!两个1表示比率,一般不用KWPAVE,11!移动工作平面到关键点11WPRO,,-90.000000,!旋转工作平面VSBW,3!由工作平面剖分体3WPCSYS,-1,0!通过当前坐标系的R-0来定义工作平面,-1表示图形窗口编号,0为坐标系的编号。KWPAVE,1!移动工作平面到关键点1WPRO,,,-90.000000!旋转工作平面VSBW,ALL!工作平面剖分体ESIZE,500,0,!设置单元尺寸大小!注意:设定单元局部坐标VATT,1,1,1,11!给体分配属性,材料号为1,实常数为1,TYPE号为1,坐标系号为11。!*VSWEEP,ALL!划分网格,体扫掠划分/DEVICE,VECTOR,1!绘制等值线/ESHAPE,1.0/REPLO!注意:红色代表最大配筋方向,绿色代表其次,蓝色表示最小配筋方向(3)/clear――beam189单元的定性命令,能运行的通/filname,beam189test,1/PREP7ET,1,BEAM189!定义单元类型及其KEYOPT选项KEYOPT,1,7,2KEYOPT,1,8,3KEYOPT,1,9,3SECTYPE,1,BEAM,RECT,,0!定义第一种截面SECOFFSET,CENTSECDATA,300,500,0,0,0,0,0,0,0,0SECTYPE,2,BEAM,RECT,,0!定义第二种截面SECOFFSET,CENTSECDATA,500,500,0,0,0,0,0,0,0,0MPTEMP,,,,,,,,!定义第一种材料弹性模量,泊松比,密度MPTEMP,1,0MPDATA,EX,1,,2.06e5MPDATA,PRXY,1,,0.3MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,7.85e-9K,1,0,0,0,!定义关键点KEYPOINTK,2,6000,0,0,K,3,6000,4000,0,K,4,0,4000,0,K,5,0,0,-5000,K,6,6000,0,-5000,K,7,6000,4000,-5000,K,8,0,4000,-5000,LSTR,1,2!通过KEYPOINT定义线LINE,梁LSTR,2,3LSTR,3,4LSTR,4,1LSTR,1,5!通过KEYPOINT定义线LINE,柱LSTR,2,6LSTR,3,7LSTR,4,8LSEL,S,,,1,2,1,!给LINE赋属性材料,截面,单元类型LATT,1,,1,,6,,1LSEL,S,,,3,4,1,LATT,1,,1,,8,,1LSEL,S,,,5,,,LATT,1,,1,,8,,2LSEL,S,,,6,,,LATT,1,,1,,7,,2LSEL,S,,,7,,,LATT,1,,1,,6,,2LSEL,S,,,8,,,LATT,1,,1,,5,,2ALLSEL,ALLLESIZE,ALL,1000,,,,,,,1!指定划分单元长度LMESH,ALL!划分单元/soluNSEL,S,LOC,Z,-5000!按位置选择需要施加约束的节点D,ALL,,,,,,ALL,,,,,!柱底施加固定约束ACEL,0,0,9810,!施加重力荷载ESEL,S,SEC,,1!选择全部梁SFBEAM,ALL,1,PRES,-500,-500,,,,,!对全部梁施加压力PRESSUREALLSEL,ALLSOLVE!求解/post1!绘制位移结果PLNSOL,U,Z,0,1/ESHAPE,1!绘制应力结果(4)finish――――带有三个圆孔的平面支座的分析,可以运行的通/clear,start/filname,ex2_2/title,2dplanestressbracket/prep7blc4,0,0,80,100cyl4,80,50,50cyl4,0,20,20cyl4,0,80,20blc4,-20,20,20,60aadd,allcyl4,80,50,30cyl4,0,20,10cyl4,0,80,10asba,6,1asba,4,3asba,1,2et,1,plane82keyopt,1,3,3r,1,20mp,ex,1,200000mp,prxy,1,0.3aesize,all,7amesh,allfinish/soluantype,0dl,7,,all,0dl,8,,all,0dl,9,,all,0dl,10,,all,0dl,11,,all,0dl,12,,all,0dl,13,,all,0dl,14,,all,0fk,9,fy,-1000solvefinish/post1plnsol,s,eqvancntr,10,0.15finish(5)/title,bucklingofabar―――杆件的屈曲分析,可以运行的通/prep7et,1,plane42,,,3r,1,0.01mp,ex,1,2.1e11mp,prxy,1,0.3n,1n,21,,1fillngen,2,21,1,21,,0.01e,22,23,2,1egen,20,1,1finish/solud,1,all,,,22,21f,21,fy,-0.5,,42,21pstrss,onsolvefinish/soluantype,bucklebucopt,subsp,1mxpand,1solvefinishFinish―――――在接头固定处定义预紧,能运行的通/clear/prep7/title,SampleapplicationofPSMESH!标题,PSMESH运用示例et,1,92mp,ex,1,1e7MP,PRXY,1,0.3mp,alpx,1,1.3e-5mp,ex,2,3e7MP,PRXY,2,0.3mp,alpx,2,8.4e-6tref,70/foc,,-.09,.34,.42/dist,,.99/ang,,-55.8/view,,.39,-.87,.31/pnum,volu,1/num,1cylind,0.5,,-0.25,0,0,180cylind,0.5,,1,1.25,0,180cylind,0.25,,0,1,0,180wpoff,.05wpstyle,,,,,,,,0cylind,0.35,1,0,0.75,0,180wpoff,-.1wpstyle,,,,,,,,0cylind,0.35,1,0.75,1,0,180vglue,allnumc,allvplotmat,1smrt,offvmesh,4,5mat,2vmesh,1,3/pnum,mat,1eplotpsmesh,,example,,volu,1,0,z,0.5,,,,nCM,lines,LINE/dist,,1.1cmplot/solueqslve,pcg,1e-8asel,s,loc,yda,all,symmasel,alldk,1,uxdk,12,uxdk,1,uzsload,1,9,tiny,force,100,1,2/title,SampleapplicationofPSMESH-preloadonlysolve!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!Finally,weconstructtheactualsolutionofinterest.Wewantto!knowwhathappenstothepreloadinthebolt,andthestressfieldaround!it,whentheassemblytemperaturerisesto150F.!Boththepreloadandthestressesincreasebecause,forauniform!temperaturerise,thereisgreaterthermalexpansioninthealuminumplates!thaninthesteelbolt.Anymethodforapplyingpreloadtha