在研究非线性迭代如何做到和Abaqus完全一致，就做了一个简单的非线性弹簧，有两个地方Abaqus不一样，一直没有研究出来为什么？请教一下。

非线性弹簧本构关系如下：

也就是说U=1时，拉力为595,U=2时拉力为1128.

在弹簧右端加一个力为5000，初始增量为0.1

第一个增量步：

查看msg，第一个迭代步得到的AVERAGE FORCE =188，第一个力增量=500，

第一个问题：那么这个188怎么从500得到呢？

第二个增量步：

第一个增量步一个迭代步就收敛了，第二个增量步开始

在第二个迭代步，力增量为500，位移增量为0.92，那么Abaqus计算的K应该就是500/0.92=543.4783。

但这个与本构关系不一致。

因为位移肯定在2以内，2以内的本构关系算出的K为

0-1s内，K=595,

1-2s内，K=1128-595=533。

第二个问题：Abaqus的K与595和533都不一致，它怎么算出来的？

附上6.12-1版本的cae文件和msg，注意，无法分析成功，但通过msg可以查看前两个分析步的Abaqus结果。

SpringNonLinearCAE.rar

两个问题能解决任何一个都非常非常感谢。