反复加载-涉及弹塑性问题

xiaoxu

借用网上的一个例子，想做个用刚性滚子在圆心处加力，对板材进行反复加载，借用的例子中，能进行一次加载和卸载的分析，可是我想在多加上几个加载和卸载的情况，以便分析多次加载对板材变形的影响，发现后加的几个加载和卸载不收敛，恳请高手指点。（存在塑性变形，非线性）
原版命令流如下：
/PREP7
ET,1,182,1,,2
ET,2,169
ET,3,171
R,2,0,0,1,0.01
RMORE,0.001,0.0001

MP,EX,1,160E9
MP,NUXY,1,0.33
MP,MU,1,0.3
TB,MISO,1,0,9
TBPT,,0.000625,100E6
TBPT,,0.0025,150E6
TBPT,,0.005,210E6
TBPT,,0.01,290E6
TBPT,,0.015,326E6
TBPT,,0.02,347E6
TBPT,,0.04,362.5E6
TBPT,,0.1,390E6
TBPT,,0.2,402.5E6
K,1
K,2,0.03
K,3,0.1
KGEN,2,1,3,1,,0.025,,3
A,1,2,5,4
A,2,3,6,5
！选择单元类型
TYPE,1
！选择材料模型
MAT,1
LESIZE,4,,,10
LESIZE,2,,,10
LESIZE,6,,,10
LESIZE,1,,,12
LESIZE,3,,,12
LESIZE,5,,,12,5
LESIZE,7,,,12,1/5
AMESH,ALL
！对称生成面
ARSYM,X,ALL     !以垂直x轴的平面为对称面
！进行节点、线、关键点等的压缩
NUMMRG,ALL
！不选择任何面
ASEL,NONE
！不选择任何线
LSEL,NONE
！定义局部坐标系,移动工作面
WPAVE,,0.1
CSWPLA,11,1  !定义柱坐标,设置工作平面坐标系
！生成关键点
K,20,0.075,180
K,21,0.075,-90
K,22,0.075,0     !生成圆弧面上的关键点
K,23,0,0         !生成控制关键点   
！生成曲线         
L,20,21
L,21,22
！定义单元类型
TYPE,2
！定义实常数
REAL,2
！对线划分网格
LMESH,ALL
！对关键点23进行网格划分
KMESH,23
！选择关键点
KSEL,S,KP,,23
！选择关键点上的节点
NSLK,S
！将关键点中最大的命名为N_LOAD
*GET,N_LOAD,NODE,,NUM,MAX
！选择所有
ALLSEL
！转换回到全局笛卡儿坐标系
CSYS,0
NSEL,S,LOC,Y,0.025
！定义单元类型
TYPE,3
！生成接触单元
ESURF
FINISH
/SOLU
！选择Y＝0处的节点
NSEL,S,LOC,Y,0
！给选择的节点进行约束
D,ALL,UY
！选择X＝0处的节点
NSEL,S,LOC,X,0
！从上面X＝0的节点当中选择Y＝0的节点
NSEL,R,LOC,Y,0
！定义约束
D,ALL,UX
！选择所有节点
NSEL,ALL
！进入求解器
！选择静态求解
ANTYPE,STATIC
！选择大变形分析
NLGEOM,ON
！选择自动时间步
SOLC,ON
！定义输入时间步
NSUBST,20,500,10
！定义输出文件
OUTRES,ALL,ALL
！施加载荷
F,N_LOAD,FY,-2000000*9.8
TIME,1              !设定载荷步终止时间
SOLVE
！撤销压力
F,N_LOAD,FY,0
TIME,2
SOLVE
FINISH
........
打算加上的部分：
F,N_LOAD,FY,-2000000*9.8
TIME,3            !设定载荷步终止时间
SOLVE
！撤销压力
F,N_LOAD,FY,0
TIME,4
SOLVE
F,N_LOAD,FY,-2000000*9.8
TIME,5             !设定载荷步终止时间
SOLVE
！撤销压力
F,N_LOAD,FY,0
TIME,6
SOLVE
FINISH
FINISH
....
这样加上以后发现不行，不能收敛，恳请有经验的指点。

xiaoxu

自己顶，大家来帮忙啊

xiaoxu

自己顶，大家来帮忙啊

xiaoxu

dajihao

eoa-kk

共同讨论讨论