非线性有限元 第2章非线性代数方程组解法.pdfVIP

第二章 非线性代数方程组的解法 在非线性力学中，有多种类型的非线性问题，如材料非线性、几何非线性、接触非线性 等。无论是哪一类非线性问题，经过有限元离散后,它们都归结为求解一个非线性代数方程 组：    ψ (  ) 0 1 1 2 n ψ (    ) 0 2 1 2 n  ψ (    ) 0 n 1 2 n 其中 , , , 是未知量，ψ , ψ , ,ψ 是 , , , 的非线性函数，现引用矢量记号 1 2 n 1 2 n 1 2 n δ [    ]T 1 2 n ψ [ψ ψ  ψ ]T 1 2 n 上述方程组可表示为 ψ(δ) 0 还可以将它改写为 ψ(δ) F (δ) R K (δ)δR 0 K (δ) 是一个n n 的矩阵，其元素k 是矢量δ的函数，R 为已知矢量。在位移有限元中, δ ij 代表未知的结点位移，F (δ) 是等效结点力，R 为等效结点荷载，方程ψ(δ) 0 表示结点 的平衡方程。 在线弹性有限元中，线性代数方程组 KδR 0 可以毫无困难地求解，但对非线性方程组ψ(δ) 0 则不行。一般来说，难以求得其精确解， 通常采用数值解法，把非线性问题转化为一系列线性问题。为了使这一系列线性解收敛于非 线性解，曾经有过许多方法，但这些解法都有一定的局限性。某一解法对某一类非线性问题 有效，但对另一类问题可能不合适。因而，根据问题性质正确选用求解方法成为非线性有限 元的一个极重要的问题。本章将介绍有限元分析中常见的各种求解非线性方程组的数值方 法。 2.1 迭代法 前面已经提到，目前求解非线性方程组的方法一般为线性化方法。若对总荷载进行线性 化处理，则称为迭代法。 2.1.1 直接迭代法 对非线性方程组 K (δ)δR 0 (2-1) 设其初始的近似解为 0 δ δ ，由此确定近似的K 矩阵 0 0 K K (δ ) 根据式〈2-1 〉可得出改进的近似解 1 0 1 δ (K ) R 重复这一过程，以第 i 次近似解求出第 i ＋1 次近似解的迭代公式为 7