有限元第七章材料非线性问题的有限元

1、第七章材料非线性问题的有限元法7.1引言前面各章所讲述的问题，都属于线性变形体系。所谓线性变形体系是指位移与载荷呈线性关系的体系，而且当载荷全部撤除后，体系将完全恢复原始状态。这种体系也称为线性弹性体系，它需满足下列条件：（1）材料的应力与应变关系满足虎克定律：（2）位移是微小的；（3）所有约束均为理想约束。在分析线性弹性体系时，可以按照体系变形前的几何位置和形状建立平衡方程，并且可以应用叠加原理。根据这种理论建立起来的方程是线性的，对于小应变和小位移的情形这种分析是适用的。实际结构的位移与载荷可以不呈线性关系，这样的体系称为非线性变形体系。如果体系的非线性是由于材料应力与应变关系的非线性引起

2、的，则称为材料非线性，如材料的弹塑性性质、松驰、徐变等。如果结构的变位使体系的受力发生了显著的变化，以至不能采用线性体系的分析方法时就称为几何非线性，如结构的大变形、大挠度的问题等。还有一类非线性问题是边界条件非线性，或状态非线性，如各种接触问题等。但本书只讨论前两类非线性问题的有限元解法，即材料非线性和几何非线性问题的有限元解法，对接触问题的有限元解法，读者可参考其它书籍。材料非线性问题的处理相对比较简单，通常不必修改整个问题的表达式，而只需将应力应变关系线性化，求解一系列的线性问题，并通过某种校正方法，最终将材料特性调整到满足给定的木构关系，从而获得了问题的解。对于几何非线性问题，那就需要

3、对公式进行根本的修改，这个问题将在后面详细讨论，不过应该指出，用于求解材料非线性问题的基本迭代方法也同样适用于几何非线性问题的求解。事实上，有些工程结构问题同时具有这两类非线性性质，它们可以统一地加以处理。本章将首先介绍用有限元方法处理非线问题的一般方法，然后讨论这些方法在非线性弹性、弹塑性和蠕变问题中的应用。在介绍弹塑性问题的处理方法前，为便于讨论，需扼要叙述一下Mises屈服准则和Prandtl-Reuss塑性流动理论，并据此写出弹塑性矩阵表达式。最后对平面刚架的极限分析做了简要介绍。7.2非线性问题的一般处理方法非线性问题用有限元法离散化应得到如下形式的一组代数方法：或写成 （7.1）其

4、中。虽然线性方程组直接求解并无困难，但对于方程组（7.1），单元刚变矩阵是单元节点位移向量的函数，直接求解就行不通。然而，下面介绍的非线性方程组的各种解法，仍以反复地求解线性方程组去获得满足一定精度要求的非线性方程组的解答。对于方程（7.1） （7.2）最简单的求解方法是直接迭代法。开始求解时先假定一组初始值代入上式的中，可求得改进了的一次近似值式中 重复上述过程，将迭代格式写成迭代一直进行到误差的某种范数小于预先规定的容许值er，即满足则停止。可以看出，该法的每一次迭代都需形成一次系数矩阵，并求解一次线性代数方程组。这里还隐含着一个假定，系数矩阵可以表示成的显函数，因此该法只适用于与变形历史

5、无关的非线性问题，例如非线性弹性问题及可利用形变理论分析的弹塑性问题。而对于依赖于变形历史的非线性问题，直接迭代法是不适用的，例如加载路径不断变化或涉及卸载及反复加载等必须利用增量理论分析的弹塑性问题。图7.1显示了单变量问题中这种迭代过程收敛和发散的可能性。通常，如曲线是凹的，则迭代发散。 图7.1 直接迭代法若已获得方程（7.1）的第n次近似解，为了求得改进的近似解，可利用仅保留线性项的Taylor级数展开式 （7.3）则有 （7.4）前式中为切线矩阵，即 （7.5）于是从式（7.3）可以得到式（7.4）中的 （7.6）重复上述迭代过程，直至达到所要求的精度。Newton-Raphson方

6、法的迭代过程示于图7.2中，通常迭代过程是收敛的。但当所选取的初始值偏离真实解较大时，正如图7.2（b）所表示的那样，发散也是可能的。 图7.2Newton-Raphson迭代法由式（7.3）看出，该法在每项迭代中必须重新形成一个系数阵并求解一次线性代数方程组。应该指出，如果原始的离散化方程组是通过变分原理导出的，则切线刚度矩阵总是对称的，而利用直接迭代法，系数矩阵的这种对称性不一定能保持。Newton-Raphson方法（简称修正的NR法）为克服Newton-Raphson方法中每次迭代都需形成一次系数矩阵和求解一次方程组的缺点，通常切线矩阵总是采用它的初始值，即令 （7.7）因此式（7.6

7、）现在变成 （7.8）这样，每次迭代求解的是系数矩阵相同的方程组。对于这种系数矩阵不变的线性方程组，如果我们在迭代的开始就将求逆或分解，则以后每次迭代只须进行一次回代过程。因而每次迭代的计算工作量大大地减小了，但收敛速度却变得较慢。不过从总的效果看，还是合算的。图7.3表示了这种方法的迭代过程。一种改进方案是在迭代若干次后就将切线矩阵修正一次，修正到当前的值，这一改进方案有时是非常有效的。 图7.3修正的Newton-Raphson迭代法为了便于理解，假定方程（7.1）式表达的是结构应力分析问题，其中代表结构的位移，代表结构的载荷。所谓增量法首先将载荷分为若干步：，相应的位移也分为若干步：。每

8、二步之间的增长量称为增量。增量解法的一般做法是假设第m步载荷和相应的位移为已知，而后让载荷增加为，再求解。如果每步载荷增量足够小，则解的收敛性是可以保证的。同时，可以得到加裁过程的中间结果。为了说明这一方法，可将式（7.1）改写成 （7.9）其中是描述载荷变化的参数。将上式对入求导，得到由此得出 （7.10）式中仍为切线矩阵。有多种方法可用来求典型常微分方程（7.10）的解。其中最简单的是Euler法，它可表达成 （7.11）其中下标指示增量加载的次数，即或在实际的数值计算中，由于并不用式（7.9），而是用其增量形式（7.10），所以解常会发生漂移现象。为了克服这一缺点，可将N-R法或修正的N

9、-R法用于每一增量步。如采用N-R法，则对于的m+1次增量步，和N-R法的第n+1次迭代表达式可写成（7.13）式中是时第n次改进的切线矩阵。由上式解出，则 （7.14）开始迭代时，取，然后连续地进行迭代，直至方程（7.9）在时能够在规定误差范围内被满足。由式（7.13）可以看出，当采用N-R法进行迭代时，每次迭代都需重新形成，并求解一次代数方程组，致使计算工作量很大。因此通常采用修正的N-R法，这时式（7.13）中的应代之为 （7.15）如果每一增量步采用N-R法进行一次迭代，由式（7.13）则有 （7.16）且假定在前一增量步结束时支配方程（7.9）是精确满足的，即则有 （7.17）这实际

10、上就就是Euler法，即式（7.11）。不过采用式（7.16）时，能将前一步支配方程的误差在本次增量步中加以校正，这是带自校正的增量法。采用此法解方程时，解的漂移不严重。本节所讨论的上述算法是目前用于求解离散的非线性方程组的常用算法。由于用有限元分析非线性问题计算工作量很大，且有时收敛很慢甚至会导致解的发散，因而引起许多计算工作者的关注。一些加速收敛的措施和方法，好的修正计算方案已相继提出。读者如有需要可查阅有关文献。在以上介绍的各种解法中，很难说哪一种算法最好。因为在某种情况下最经济有效的方法，在另一种情况下则不然，甚至解收敛很慢或不收敛。不过，要为一个通用程序编入一种解法时，增量法是合宜的

11、。因为只要选择足够小的增量步，解总是收敛的，如可能，还可在每一增量步中采用N-R法或修正的N-R法，使计算结果满足一定的要求。7.3非线性弹性力学问题如果我们只考虑小变形，则平衡方程在整个求解域内可写成 （7.18）上式中的积分运算实际上应逐个单元进行，并按单元集成法把它们对节点的平衡的贡献进行叠加。右端节点力向量也由单元集成法形成。几何关系可写成 （7.19）但此时物理关系是非线性的，一般可写成 （7.20）由式（7.18）至（7.20）可以导出与非线性方程组（7.1）相同的表达式因而上节中所讨论的各种解法，原则上都可在此应用。根据非线性弹性力学问题中非线性应力一应变关系的不同表达方式，可以

12、产生下列几种求解方案。直接迭代法如果材料的应力一应变关系能表达成于是由式（7.19）可将应力写成将上式代入式（7.18）后，得 （7.21）其中用迭代法求解方程（7.21）时，总是首先取，然后求出，由求出位移的第一次近似值。重复这一过程，迭代格式为 （7.22）切线刚度法（即N-R法）如果材料的应力一应变关系能表示成增量关系 （7.23）式中为切线弹性矩阵，可将式（7.18）改写为 （7.24）将对求导，注意到关系式（7.23），则有 （7.25）式中 （7.26）是切线刚度矩阵。在N-R法迭代中，可首先取，由式（7.19）求得，再由式（7.23）确定。将代入式（7.26）求出，并由式（7.2

13、4）求出，由。解出，于是得到位移的第一次近似值为。重复上述过程，可将迭代公式表示为 （7.27）式中 （7.28）我们知道，式（7.28）右端的第一项是与相等价的节点载荷，因而可理解为结构上的不平衡节点载荷向量，每迭代一次相当于对不平衡节点载荷作一次矫正。初应力法初应力法实质上是修正的N-R法的具体应用。如果材料的应力一应变关系可以表示成 （7.29）即应力分量能由给定的应变分量确定，我们就可以用具有初应力的线弹性应力一应变关系 （7.30）去代替上式。其中D采用非线性应力一应变关系式在时的切线弹性矩阵。这就是说，可用调整初应力值的方法，使在给定应变下，用式（7.30）算得的应力与用式（7.2

14、1）算出的结果相同。于是有 （7.31）其中。在一维情况下初应力的定义见图7.4 图7.4初应力法如令 （7.32）为结构的起始切线刚度矩阵，则将式（7.30）代入（7.18），得 （7.33）式中 （7.34）为与初应力等价的节点载荷。根据式（7.33），可以采用如下迭代步骤：首先由下式计算第一次近似位移由求出，用式（7.31）求出初应力由式（7.34）求出用对位移进行一次矫正于是位移的第二次近似值是继续进行这种矫正，直至收敛。由此可见，用初应力法迭代时，每一步都是根据真实应力与弹性应力之差决定初应力的，并从而决定不平衡的节点载荷，依次进行一次调整，使位移进一步逼近真实值。初应力法迭代过程在

15、单向应力状态下的应力变化示于图7.5中。 图7.5初应变法对于某些问题，例如蠕变问题，其应变由应力值决定 (7.35)这时，上式可用具有初应变的线弹性应力一应变关系 （7.36）来代替，即调整初应变值，使在给定应力下，用式（7.36）求得的应变与用式（7.35）求得的结果一样。于是有 （7.37）式中。在单向应力的情况下初应变的定义见图7.6。 图7.6将式（7.36）代入（7.18），则有 （7.38）其中 （7.39）为与初应变等价的节点载荷。见式（7.32）。根据式（7.37），可以采用与初应力法类似的迭代过程。先计算位移的第一次近似值由计算应力，由式（7.37）确定初应变，再由式（7.

16、39）计算R1，然后用求出，对位移进行一次矫正重复上述过程，直到收敛。可见，初应变法的每一步迭代是根据真实应变与弹性应变之差，去对位移进行调整的。在单向应力状态下初应变迭代过程示于图7.7中。 图7.77.4弹塑性应力一应变关系在弹塑性小变形情况下，弹性力学中的平衡方程和几何关系仍然成立，但描写应力一应变关系的物理方程却是非线性的。材料的塑性性质对于大多数金属材料，单向拉伸试验的应力应变曲线如图7.8所示。当应力小于屈服极限时，应力和应变呈线弹性关系（一般假定比例极限与屈服极限重合）；而应力达到屈服极限后，应力一应变显示出非线性关系。图7.8（a）所示为应变硬化材料，图7.8(b)所示为理想弹

17、塑性材料。由实验知，对于硬化材料，应力超过屈服极限时，要增加应变，应力也需增加。而对于理想弹塑性材料，应力达到时，材料则发生流动变形。然而无论哪种情况，当应力达到后，卸载都是弹性的，卸载路经沿直线，其斜率大致与初始加载路径相同。完全卸去载荷后，试件会留下残余变形。因此，弹塑性材料当应力达到或超出时，应力一应变之间并无一一对应关系，即应变不仅决定于当时的应力，而且还与整个加载的历史有关。这是与上节讨论的非线性弹性材料的区别所在。还应当指出，材料进入塑性后，经卸载再加载时，应力一应变关系沿CBD变化。对于硬化材料显然屈服应力提高了。材料进入塑性后，经卸载并反向加载，材料会再次屈服，对于各向同性材料

18、，再次屈服的应力的绝对值大体上与开始卸载时的应力相等。 (a) (b) 图7.8Mises屈服准则Mises屈服准则假定，材料在复杂应力状态下的形状改变能达到单向拉伸屈服时的形状改变能时，材料就开始屈服。于是，Mises屈服条件可写成式中为主应力，是单向拉伸屈服极限。如果定义等效应力则Mises屈服条件为 （7.40）在一般应力状态下，等效应力可表示成引进应力偏量 （7.41）式中是平均应力，则等效应力可用应力偏量表示为（7.42）若记则可将等效应力简洁地表示成 （7.43）对于大多数金属材料，实验已经证明，Mises准则能较好符合实际情况。因此，下面将结合Mises准则来讨论。首先，讨论复杂

19、应力状态下的强化规律。假定材料进入屈服后，载荷按小增量方式逐步加载，在一个载荷增量作用下应力和应变都会增加一微小的增量和。此时总应变增量可分成弹性的和塑性的两部分 （7.44）弹性应变增量在卸载后可完全恢复，而塑性应变增量在卸载后则不能恢复，如图7.9所示。 图7.9与等效应力对应，定义等效应变 （7.45）对于单向拉伸，于是等效应变恰等于拉伸应变。对应于塑性应变增量的等效应变称为塑性等效应变增量，记作。由于塑性变形的泊桑比，于是有(7.46)注意到塑性变形中的体积应变等于零，上式还可表达成 （7.47）若记 （7.48）则等效塑性应变增量可改写成 （7.49）对于一般应力状态，实验资料证明了

20、如下的应变强化规律：材料进入屈服以后进行卸载或部分卸载然后再加载，其新的屈服应力值仅与卸载前的等效塑性应变总量有关。这就是说，只有当等效应力适合 （7.50）时，重新屈服才会发生。这里的函数H反映了新的屈服应力对等效塑性应变总量的依赖关系。由于单向拉伸时，就是拉伸应力，就是拉伸塑性应变增量，所以式（7.50）中的函数关系可以通过单向拉伸的和之间的关系来确定，如图7.10所示。式（7.50）反映了屈服和强化之间的关系，称为等向强化材料的Mises准则。 图7.10Prandtl-Reuss塑性流动理论如果将等向强化Mises准则的式（7.50）写成 （7.51）则F可以看成n维应力空间的一个曲面

21、，称为屈服面。屈服面的位置决定于当时材料中的等效塑性应变总量。对于金属一类的材料，理论和实验的广泛研究表明，塑性应变增量和屈服面之间存在如下关系 （7.52）其中是一个待定的比例常数。上式可以解释为塑性应变增量“向量”垂直于n维应力空间的屈服面，如图7.11所示。因此，式（7.52）称为法向流动法则。 图7.11现在来决定式（7.52）中的常数。为此先计算。由式（7.42），并注意到，就有将上式代入（7.52），并注意到式（7.48），得上式等号两边的模应相等，于是有由式（7.49）和（7.43），从上式可得出由于在加载时应取正值，于是有这样，向流动法则最终化为 （7.53）应力一应变关系如前

22、所述，当应力产生一无限小增量时，总应变增量可分解成弹性的和塑性的两部分，即弹性应变增量与应力增量是线性关系，可写成 （7.54）式中De是弹性矩阵。用前乘上式的两边，得如果将强化材料的Mises准则（7.50）写成微分的形式 (7.55)并利用式（7.53），则上式可化为 由此可得到等效塑性应变增量 和总应变增量 的关系式将式（7.53）代入（7.54），并应用上式，得到记 （7.56） （7.57）于是得到增量形式的弹塑性应力一应变关系 （7.58）通常称为弹塑性矩阵。顺便指出，对于理想弹塑性材料，上述应力一应变关系（7.58）仍可应用，不过此时取。弹塑性矩阵表达式上节已导出了弹塑性矩阵的一

23、般表达式，为便于应用。下面对三维空间问题，轴对称问题和平面问题分别写出它的显式。1三维空间问题的弹塑性矩阵利用空间问题的弹性矩阵表达式（1.7），并注意到，容易看出 =其中G为剪切模量。于是注意到由式（7.56）容易得 (7.59)将式（1.7）和（7.59）代入式（7.57），得到空间问题的弹塑性矩阵为（7.60）式中 （7.61）2轴对称问题的弹塑性矩阵我们记在式（7.60）中划去最后二行二列，可以得到轴对称问题的弹塑性矩阵 （7.62）式中仍由式（7.61）决定，但等效应力为3平面问题的弹塑性矩阵对于平面应力问题，应力增量和应变增量分别为等效应力为由此从而应用平面应力问题的弹性矩阵（2.

24、6）式，有将上两式代入（7.56），得（7.63）式中 （7.64）由式（7.57），经整理，将平面应力问题的弹塑性矩写成（7.65）其中 （7.66）而Q的表达式（7.64）可改写为对于平面应变问题，材料进入塑性后，等效应力可写成这时的弹塑性矩阵可直接由式（7.65）得出，只要将近其中E换成换成。切线模量的计算上节导出了几种情况下的弹塑性矩阵，它们的表达式中都含有切线模量，它是表示材料硬化性能的参数。由式（7.55），微分形式的Mises屈服准则给出，由此可以看出是等效应力相对于塑性等效应变的变化率。它必须通过单向拉伸试验所给出的应力一应变由曲线来确定。下面推导的计算公式。设单向拉伸试验时，

25、材料进入塑性后的应力一应变关系以某种函数的形式给出，如图7.12所示。 （7.67） 图7.12由于单向拉伸试验给出的是全应变，为了建立与塑性应变之间的关系，必须将它分解为弹性部分与塑性部分之和，且弹性部分服从虎光定律，故有将其代入按试验曲线给定的函数关系式上式两边取微分，则有式中代表应力一应变曲线在点的斜率。上式经整理可得对于单向拉伸，显然即，故正是对应于复杂应力状态下的，代入上式则可得到 （7.68）这就是计算的公式。的计算步骤是，首先求出等效应变，然后在给定材料的应力一应变曲线上求得该点所对应的曲线斜率，最后代入式（7.68）就可求得。这里有几种情况应该注意，一种情况是，在弹性范围内不能

26、计算的值，因此时它没有塑性变形，故应取值，在编写程序时，可根据计算机的情况取一个适当大的数即可表示。另一种情况是，对理想弹塑性材料，当材料进入屈服时，由，此时可由式（7.68）知，应取O值。最后，如果材料硬化程度不剧烈，即材料的应力一应变曲线较平缓时，则可考虑用以下的近似公式计算： （7.69）求出后，就可以利用节的有关公式，计算不同情况下的弹塑性矩阵。7.5弹塑性问题的求解方法弹塑性问题的有限元求解方法完全是在线弹性问题有限元求解方法的基础上发展起来的，所以本章前面各章所介绍的线弹性问题有限元公式仍然适用。只是在塑性区内应力和应变不再成线性关系，必须采用应力微分和应变微分的关系式（7.58）

27、，由于式（7.58）右端的系数矩阵与当时的应力水平有关，所以这个关系是非线性的。用有限单元法求解弹塑性问题，也是用适当的方法将问题线性化，然后用一系列的线性解去逼近一个非线性问题的解。而每一个线性解的求解方法和过程和线弹性问题的处理非常类似。为达到线性化的目的，可以采用增量法（或称增量加载法）。一般的做法是按比例施加载荷，将结构的弹性极限载荷作为第一个增量，其余的载荷再分成若干等分。如果实际载荷不是按比例施加的，可根据实际情况确定载荷增量。当材料进入塑性后，只要载荷增量适当地小，由式（7.58），应力增量和应变增量的关系可近似地表示成 （7.70）且可认为仅与加载前的应力水平有关，而与应力和应

28、变的增量无关。这样式（7.70）就可视为线性的。增量切线刚度法首先，可以进行一次线弹性的分析，得到弹性极限载荷下结构的位移，应变和应力，分别记为，和，在此基础上将载荷分为n个增量，并相继作用于结构。作用载荷增量时，对于应力处于弹性状态的单元，单元刚度阵用 (7.71)计算。而对于进入塑性的单元，单元刚度阵用 (7.72)计算。中的应力取增量加载前的应力。结构总刚度矩阵的形成与前几章所介绍的线弹性问题 方法完全相同，也是利用单元定位向量，直接由单元刚度矩阵进行集成，由此求解可得到，从而求出和，这样第一次增量加载后的位移，应变和应力是在加载第k次载荷增量后，可通过求解 （7.73）得到，进而得重复

29、上述过程，直至第n个载荷增量。最后的位移、应变和应力就是弹塑性分析的结果。增量切线刚度法，由于每个载荷增量步都需要重新计算一次刚度矩阵，所以这个方法也称为变刚度法。增量切线刚度法在实施过程中，随着不断地增量加载，解会越来越偏离真实解，即出现所谓解的漂移现象。产生这种现象的原因除了应力一应变关系用增量形式代替微分形式外，还来自过渡单元的不准确的刚度矩阵的计算。必要时可以用下面介绍的方法加以校正。结构在逐步加载过程中，塑性是不断扩展的。有一些单元在某一增量载荷步前处于弹性区内，而在该增量载荷步后可能进入塑性区，这些单元称为过渡单元。对于过渡单元，其刚度矩阵如果简单地按式（7.71）或（7.72）计

30、算都会带来较大误差。通常在计算其刚度矩阵时，应该用下述带权平均弹塑性矩阵。如用表示过渡单元达到屈服时所需的等效应变增量，而用表示该单元下次增量加载所引起的等效应变增量，如图7.13所示，记 图7.13显然，对过渡单元有。定义加权平均弹塑性矩阵则对过渡单元应按下式计算刚度矩阵 （7.74）如果采用数值积分计算，则在过渡的积分点上应作上述考虑。通常对的估计，开始往往不够精确。一般第一次估计是根据上次增量载荷的结果推出，然后用解算的结果来修改，经过二三次的迭代就可得到比较精确的结果。应该指出的是，上述用迭代处理和校正过渡单元刚度矩阵的计算，虽然提高了求解的精度，但在迭代中却付出了成倍的计算时间。下面

31、介绍另一种较好的求解方案。上述增量切线刚度法，对每一增量加载步，被求解的线性代数方程组（7.73）可以作如下的变动： （7.75）式中是本次增量加载后结构中承受的总载荷向量。而右端最后一项积分是前一增量加载步末与结构中的应力等价的结点载荷。很显然，这是我们在7.2节中所提到的带自校正的增量法的方程。由于这一方程能对以前的载荷不平衡作一校正，使求解结果具有较好的精度。因每次增量载荷步，方程（7.75）都具有“自校正”的作用，有人建议可以不去考虑上述过渡单元的迭代计算。实践证明，带自校正的增量切线刚度法（或称一阶自校正法）是一种较好的求解方案。初应力法对于弹塑性问题，增量形式的应力一应变关系可以定

32、义为而其中由式（7.56）计算。对某单元开始进入塑性以后的每次加载，用增量代替微分 （7.76）位移增量应满足的平衡方程是 （7.77）由于式（7.76）第一式中的系数阵为弹性矩阵，式（7.77）中的刚度矩阵就是弹性计算中的刚度矩阵。而是与初应力等值的节点载荷，它是不平衡的矫正载荷。可以看出，式（7.77）中的矫正载荷决定于应变增量，而在求解前是未知的，因此对每个增量载荷步，一个迭代过程是必须的，以便同时求出位移增量和应变增量。第k次增量载荷的迭代公式是 （7.78）第一次迭代是在下作纯弹性计算。以后的迭代是根据前次迭代求出的和加载前的应力水平计算，然后按式计算。逐次进行迭代，直至收敛。值得注

33、意的是，对于过渡单元，初应力的计算不应计及总应变增量中在进入屈服之前的部分（参看图7.13），即矫正载荷应用下式计算 （7.79）初应变法对于弹塑性问题，增量形式的应力一应变关系可定义为 （7.80）而由式（7.53）和（7.55），有 （7.81）用增量代替微分，将式（7.80）和（7.81）线性化 （7.82）此时位移增量应满足的平衡方程为 (7.83)式中仍为弹性计算中的刚度矩阵，而 （7.84）是与初应变等价的节点载荷，或称为矫正载荷。由于矫正载荷与应力增量有关，而本身又是待确定的量，因此求解方程组（7.83）必须采用迭代方法。第k次增量载荷步的迭代公式为 （7.85）每一次加载时，首

34、先取作一次弹性计算，以后的迭代用上次迭代的结果和加载前的应力水平计算，逐次进行迭代，直至收敛。方法的比较1增量切线刚度法是在每次增量加载中用调整刚度矩阵的办法求非线性问题的近似解的。初应力法实质上是在每一增量步上确立结构中弹塑性的应力值和弹性解应力值之差，并不断按弹性方式重新分配这个差值，以恢复平衡。而初应变法则是不断调整初应变的过程。2增量切线刚度法由于在每次增量加载中都形成一次刚度矩阵，并求解一次，故计算时间一般比初应力法和初应变法长。因初应力法和初应变法在迭代中应用的是弹性刚度矩阵，如果求解一开始就将刚度矩阵进行三角分解，则在每次迭代中只需进行右端项的约化和回代，计算时间较节省。3可以证明，对一般强化材料，初应力法的迭代过程总是收敛的，而对初应变法，收敛的充分条件是。对理想塑性材料，初应力法和初应变法的迭代过程都是发散的。4如果要为一个通用程序配置一个求解方案的话，建议采用带校正的增量切线刚度法。通常该方法能得到较好的求解精度。如果将增量切线刚度法与初应力法相结合，则效果更好。算例带孔平板的拉伸问题承受单向拉伸的开孔板条的几何尺寸和载荷情况如图7.14所示。利用对称性，计算区域取板条的四分之一部分。采用Mises屈服准则，应变硬化材料采用的是常斜率的单轴硬化曲线。图7.14(b)和(c)分别表示理想塑性和应变硬化情况下不同载荷水平塑