应力、应变均已非常接近正确值，却不收敛？

Phil_Leng

　　我做本构的二次开发，写了一段程序。调试通过后，进行计算。算例如下：仅一个单元，一面固支，其对面受均布面力，即单轴均匀受压问题。在静力状态下，荷载分20步施加，其中第一步加至材料屈服。

　　材料屈服后，第三荷载步（总第四步）不收敛。将其中计算结果输出来查看后发现，应力和应变都曾经非常接近正确解：应力正确解为2060000，计算值为2059999.997，应变的误差也在这个量级。但ADINA并未认为已经收敛，而是继续往下迭代计算，导致应力和应变值超过正确解，并最终发散。

　　查看out文件，发现在该荷载步，不平衡能量、不平衡力的范数和不平衡位移的范数从第一次迭代开始就都是越来越大，和前三个荷载步完全相反。

　　我作了以下尝试：我将材料屈服后的19步，变成1步完全加上，计算是收敛的，这排除了荷载步过大的可能性；我又将屈服后的荷载步变为原来的十分之一，仍然是在总第四步不收敛；我又尝试了在ADINA中采取不同的方程解法，采用不同的收敛控制和不同的收敛精度，情况没有得到改变。

　　看起来，似乎是在计算到第四步时有什么地方出了错。但我写的程序，只提供本构关系，并不涉及荷载步和迭代方法，而且既然前三步已经正确了，为什么第四步就不能计算正确呢？

　　各位高手，如果知道这其中的原因，或者可能的原因，请不吝赐教！

　　谢谢！

feizhang

太高深，不懂！