第三章 回弹 /units,si /prep7 et,1,plane42 keyopt,1,3,2 mp,ex,1,220e9 mp,nuxy,1,0.3 TB,MISO,1,1,25, TBTEMP,0 TBPT,,0.000814,175.88e6, TBPT,,0.001531,324.98e6, TBPT,,0.002819,591.26e6, TBPT,,0.003846,801.457e6, TBPT,,0.004239,873.22e6, TBPT,,0.004594,924.399e6, TBPT,,0.004955,956.88e6, TBPT,,0.00531,976.99e6, TBPT,,0.005699,990.2138e6, TBPT,,0.006261,1000e6, TBPT,,0.006957,1008.34e6, TBPT,,0.007804,1015.5e6, TBPT,,0.008779,1021.185e6, TBPT,,0.009827,1026.47e6, TBPT,,0.010929,1031.289e6, TBPT,,0.012066,1035.444e6, TBPT,,0.0132,1039.693e6, TBPT,,0.014354,1044.886e6, TBPT,,0.015494,1049.796e6, TBPT,,0.016631,1053.951e6, TBPT,,0.017785,1058.106e6, TBPT,,0.019039,1062.827e6, TBPT,,0.020276,1067.643e6, csys,1 k,1,0.00512 k,2,0.0064 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 nu mmrg,all nu mcmp,all sav e cp,1,u y ,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,u x nplot allsel ou tres,all,all time,1 NSUBST,100 nlgeom,on d,166,u y ,-0.00128 solv e *get,def1,node,166,U,Y nsel,s,loc,y ,90 ddele,all,u y allsel f,166,fy ,0 solv e *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,0.1 ! 给定厚度 LOCAL,11,1,,,,,,,0.5 ! 定义区域坐标系统为圆柱坐标系，Y 轴是X 的0.5 倍。 K,1,2,90 K,2,2,0 L,1,2 LOCAL,12,1,,,,,,,0.8461585 ! 定义区域坐标系统为圆柱坐标系，Y 轴是X 的0.846 倍。 K,3,3.25,90 K,4,3.25,0.0 L,3,4 CSYS,0 !回到卡式坐标系下。...
