文档介绍：第三章回弹/units,si /prep7 et,1,plane42 keyopt,1,3,2 mp,ex,1,2 .160687961 e 11 mp,nuxy,1, TB,MISO,1,1,25, TBTEMP,0 TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,1000e6, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, TBPT,,,, csys,1 k,1, k,2, kgen,36,1,2,1,,10 *do,j,1,69,2 a,j,j+1,j+3,j+2 *enddo a,71,72,2,1 save lesize,all,,,4 amesh,all nummrg,all numcmp,all save cp,1,uy,166,170,171,172,173 finish /solu nsel,s,loc,y,270 d,all,all nplot allsel nplot nsel,s,loc,y,90 d,all,ux nplot allsel outres,all,all time,1 NSUBST,100 nlgeom,on d,166,uy,- solve *get,def1,node,166,U,Y nsel,s,loc,y,90 ddele,all,uy allsel f,166,fy,0 solve *get,def2,node,166,U,Y finish /post1 ! 进入后处理器 set,list ! 显示所有计算的结果 set,1,last,1 ! 读取第一个加载步最后子步的结果 pldisp,1 ! 同时显示变形前与变形后结构的形状。 plnsol,s,eqv,1,1 ! 显示节点的等效应力 set,2,last,1 ! 读取第二个加载步最后子步的结果 pldisp,1 ! 同时显示变形前与变形后结构的形状。 plnsol,s,eqv,1,1 ! 显示节点的等效应力自适应划分/PREP7 SMRT,OFF ANTYPE,STATIC ET,1,PLANE42,,,3 ! 定义平面应力单元有单元厚度 MP,EX,1,210E9 ! 定义材料属性 MP,NUXY,1,.3 R,1, ! 给定厚度 LOCAL,11,1,,,,,,, ! 定义区域坐标系统为圆柱坐标系, Y 轴是 X的 倍。 K,1,2,90 K,2,2,0 L,1,2 LOCAL,12,1,,,,,,, ! 定义区域坐标系统为圆柱坐标系, Y 轴是 X的 倍。 K,3,,90 K,4,, L,3,4 CSYS,0 ! 回到卡式坐标系下。 L,2,4 L,1,3 AL,2,4,1,3 ! 创建面积 DL,4,1,SYMM ! 创建对称面约束 DL,3,1,SYMM SFL,2,PRES,-10E6 ! 定义线上的分布载荷 MSHK,2 ! 混合网格划分 MSHA,0,2D ! 四边形网格 SAVE ! 备份数据库 FINISH ADAPT,7,7,,,1 ! 求解次数为 nsoln=4 次, 能量误差百分比 stargt=7, 单元尺寸改变最大比例因子 facmx=1 /POST1 EPLOT ! 绘制单元 PRERR ! 打印能量误差百分比 NSEL,S,LOC,Y, NSEL,R,LOC,X, *GET,MNODE,NODE,,NUM,MAX *G