第二章非线性代数方程组的解法 - Read.DOC

第二章 非线性代数方程组的解法 在非线性力学中，有多种类型的非线性问题，如材料非线性、几何非线性、接触非线性等。无论是哪一类非线性问题，经过有限元离散后,它们都归结为求解一个非线性代数方程组： 其中是未知量，是的非线性函数，现引用矢量记号 上述方程组可表示为 还可以将它改写为 是一个的矩阵，其元素是矢量的函数，为已知矢量。在位移有限元中, 代表未知的结点位移，是等效结点力，为等效结点荷载，方程表示结点的平衡方程。 在线弹性有限元中，线性代数方程组 可以毫无困难地求解，但对非线性方程组则不行。一般来说，难以求得其精确解，通常采用数值解法，把非线性问题转化为一系列线性问题。为了使这一系列线性解收敛于非线性解，曾经有过许多方法，但这些解法都有一定的局限性。某一解法对某一类非线性问题有效，但对另一类问题可能不合适。因而，根据问题性质正确选用求解方法成为非线性有限元的一个极重要的问题。本章将介绍有限元分析中常见的各种求解非线性方程组的数值方法。 2.1 迭代法 前面已经提到，目前求解非线性方程组的方法一般为线性化方法。若对总荷载进行线性化处理，则称为迭代法。 2.1.1直接迭代法 对非线性方程组 (2-1) 设其初始的近似解为，由此确定近似的矩阵 根据式〈2-1〉可得出改进的近似解 重复这一过程，以第i次近似解求出第i＋1次近似解的迭代公式为 (2-2) 直到 (2-3) 变得充分小，即近似解收敛时，终止迭代。 在迭代过程中,得到的近似解一般不会满足〈2-1〉式,即 作为对平衡偏离的一种度量，称为失衡力。 图2-1 为凸曲线 图2-2 为凹曲线 对于一个单变量问题的非线性方程，直接迭代法的计算过程如图2-1和图2-2所示，它们分别给出为凸和凹曲线时的迭代过程。可以看出就是过曲线上点() 与原点的割线斜率。对于单变量问题，这一迭代过程是收敛的，但对多自由度情况，由于未知量通过矩阵耦合，迭代过程可能不收敛。 2.1.2 Newton—Raphson方法 Newton—Raphson方法是求解非线性方程组 (2-4) 的一个著名方法，简称Newton法。以下将介绍这种方法。 设为具有一阶导数的连续函数，是方程（2-4）的第i次近似解。若 希望能找到一个更好的、方程（2-4）的近似解为 (2-5) 将（2-5）代入（2-4），并在附近按一阶Taylor级数展开，则在处的线性近似公式为 其中 引入记号 假定为真实解，则由 解出修正量为 (2-6) 由于这样确定的仅考虑了Taylor级数的线性项，因而按式（2-6）和（2-5）求出的新解仍然是近似解。这样，Newton法的迭代公式可归纳为 (2-7) 对于单变量的非线性问题，其迭代过程见图2-3和2-4，可以看出是曲线上通过点的切线斜率 Newton法的收敛性是好的，但对某些非线性问题，如理想塑性和塑性软化问题，在迭代过程中可能是奇异或病态的，于是的求逆就会出现困难。为此，可引入一个阻尼因子，使矩阵或者成为非奇异的，或者使它的病态减弱。这儿为阶的单位矩阵。的作用是改变矩阵主对角线元素不占优的情况。当变大时，收敛速度变 图2-3 图2-4 慢，当→0时，收敛速度最快。引入后，将用下式代替（2-6） (2-8) 2.1.3修正的Newton-Raphson法 采用直接迭代法和Newton法求解非线性方程组时，在迭代过程的每一步都需要重新计算。如将Newton法迭代公式中的改用初始矩阵，就成了修正的Newton-Raphson法（简称修正的Newton法）。此时，仅第一步迭代需要完全求解一个线性方程组，并将三角分解后的存贮起来，以后的每一步迭代都采用公式 (2-9) 图2－5 这样，只需按式(2-9)右端的进行回代即可。 修正Newton法的每一步迭代所用的计算