ANSYS塑性变形模拟例子

/一个周边简支的圆盘，其中心受到一个冲杆的周期作用(假定冲杆是刚性的)，需要进行圆盘在冲杆的周期作用下的塑性分析。本实例的模型简图如图19.1所示，材料特性如下所示，塑性时的应力－应变关系如表19.1，载荷历史如表19.2所示。弹性模量：EX=70000，泊松比：NUXY=0.325/PREP7/TITLE,CircularPlateLoadedbyaCircularPunch-KinematicHardening!*下面定义建模分析时需要的参数EXX=70000RPL=65RPU=5H=6.5STS1=55STN1=STS1/EXXSTS2=112STN2=0.00575STS3=172STN3=0.02925STS4=241STN4=0.1NEX=15NET=2NEX1=nint(0.8*NET)NEX2=NEX-NEX1!*ET,1,42,,,1！定义单元PLANE42，设置为轴对称!*MP,EX,1,EXX！定义材料属性MP,NUXY,1,0.325!*TB,KINH,1,1,4,！定义多线性随动强化准则TBPT,,STN1,STS1TBPT,,STN2,STS2TBPT,,STN3,STS3TBPT,,STN4,STS4!*创建节点N,1,RPL,,,,,,N,2,0,,,,,,N,3,,H/2,,,,,!*创建关键点K,1,,-(H/2),,K,2,RPU,-(H/2),,K,3,RPL,-(H/2),,KGEN,2,ALL,,,,H,,3,0！复制并平移关键点!*创建直线段L,1,2L,4,5L,2,3L,5,6!*用关键点创建面A,1,2,5,4A,2,3,6,5!*对线进行网格控制LESIZE,1,,,NEX1LESIZE,2,,,NEX1LESIZE,3,,,NEX2,2.5LESIZE,4,,,NEX2,2.5LESIZE,5,,,NETLESIZE,6,,,NETLESIZE,7,,,NET!*AMESH,ALL！对面划分网格!*NUMMRG,NODE,,,,LOW！合并重复节点NUMCMP,NODE!*CP,1,UY,3,7,9！定义耦合自由度!*D,1,UY！定义位移约束条件NSEL,S,LOC,X,0D,ALL,UXNSEL,ALL!*下面开始求解/SOLU!*ANTYPE,0！指定分析类型OUTRES,ALL,ALL！设置输出控制PRED,ON！打开预测器!*施加并输出第一个载荷步F,3,FY,0TIME,1e-6AUTOTS,1NSUBST,1,,,1KBC,0LSWRITE,1,!*施加并输出第二个载荷步F,3,FY,-6000TIME,1NSUBST,10,,,1LSWRITE,2,!*施加并输出第三个载荷步F,3,FY,750TIME,2LSWRITE,3,!*施加并输出第四个载荷步F,3,FY,-6000TIME,3LSWRITE,4,!*求解所有4个载荷步LSSOLVE,1,4,1,FINISH!*进入通用后处理器/POST1SET,LAST选定最后一个载荷子步PLDISP,1绘制变形图!*PLNSOL,EPPL,EQV,1,1！绘制等效塑性应变PLNSOL,EPEL,EQV,1,1！绘制等效弹性应变!*进入时间－历程后处理/POST26ESOL,2,3,3,S,Y,SY_2！定义变量ESOL,3,3,3,EPEL,Y,EPELY_3ESOL,4,3,3,EPPL,Y,EPPLY_4ADD,5,3,4,,EPTLY！对变量进行运算!*/AXLAB,X,Y-Strain/AXLAB,Y,Y-Stress!*！绘制应力－应变曲线XVAR,5PLVAR,2,FINISH