请教塑性分析中：卸载后的残余变形；载荷步的结束时间；子载荷步；收敛

xarrnous

看了之前讨论的讨论，包括到处都能查到的那个圆盘的例子，想自己试试，但是发现有问题，希望大家能指点一下。

具体的，我做了一个长方体，xyz为10*10*50，Solid45，划分单元时单元尺寸设为1（global），弹性的材料属性 E=1000，柏松比=0.3，塑形模型选择的nonlinear-rate independent-kinematic hardening plasticity-Mises plasticity-bilinear,屈服应力10，屈服之后的模量E'=100

第一步，过原点的10*10面上，在节点上施加固定约束；过（0,0,50)的10*10面上，在节点上施加z=1的位移约束。（我认为z=0.5的时候，应变是0.01，应力是10，恰屈服，所以z=1的时候就应该是有z方向上0.5的塑形变形了），time at end of load=0.01, DOF solution predictior设为 on for all substep，Automatic time stepping设为prog chosen，Integration Pt设为No，然后Write LS File，1；

第二步，过原点的10*10面上，保留固定约束；过（0,0,50)的10*10面上的位移约束删除，time at end of load=0.1，然后Write LS File，2。其他不变。

结果就是ERROR: Solution not converged at time 1.12E-022).   

有几个问题：
1，当然是这个收敛问题怎么解决，按理说模型很简单，我觉得就是哪里没设置对，但是自己想不出来；
2，既然是rate independent, 那么time at end of load是不是就没什么意义了，我设1或是设100都可以，只要后面比前面大就行；
3，我这么算之后，第二步中能不能看到塑形变形，因为没结果，所以也不清楚到底这样做对不对了；
4，最后一个，如果说再加上第3个载荷步，重新在z正方向上加载，屈服极限是不是就已经变了，比如第二步之后过（0,0,50)的正方面已经有一个0.5的残余变形了，那么重新加载后，是不是在达到1.0（相对与最初，第一步而言）之前，弹性模量就按照1000计算了。

问题比较多，先谢谢大家了。

xarrnous

自己顶一个。之前的帖子我确实认真看过了，但是本人愚钝，遇到问题还是解决不了，还麻烦各位高手指点一下。

xarrnous

弄出来了，给大家汇报一下，可能是载荷变化太大的原因。最后我加了4个LS，分别是z方向位移1，0.95，0.8，以及只有固定端固定，另外一端自由，然后就收敛了。算出来的塑性von Mises应变恒定为0.01484（固定端附近由于柏松比引起的收缩也受到了限制，所以实际上附近的应力会稍大一点），弹性应变逐渐变小。值得说一下的是，当自由约束之后，在固定端附近仍有很小的弹性应变，0.002335，精度对我来说够了。

zhangxw023

好贴，顶一个。本人目前也在做一个卸载分析，不过我想问下楼主，你是四个载荷步一起算的，还是一个一个计算的？

zhangxw023

还有你这个塑性变形0.01484是怎么查看的？谢谢