非线性代数方程组解法.pdfVIP

第二章 非线性代数方程组的解法 从数学的角度，非线性问题的求解可以归纳为偏微分方程的边值问题和初值问题。大部 分静力问题属于边值问题，动力响应分析、热传导和流体力学问题则是初值问题。以下 2.1 至 2.3 节为非线性方程的边值解法，节 2.4 简单介绍非线性方程的初值解法。 在非线性力学中，有多种类型的非线性问题，如材料非线性、几何非线性、接触非线性 等。无论是哪一类非线性问题，经有限元离散后，它们都归结为求解一个非线性代数方程组： ψ (δ δ L)δ 0 1 1 2 n δ δ ψ (δ L) 0 2 1 2 n LL δ δ ψ (δ L) 0 1n 2 n , , δ ,δ Lδ ψ ψ , ,ψ , L , , δ ,δ Lδ 其中 是未知量， 是 的非线性函数，现引用矢量记号 1 2 n 1 2 n 1 2 n δ δ δ [δ L ]T 1 2 n ψ ψ ψ ψ [ L ]T 1 2 n 上述方程组可表示为 ψ(δ) 0 还可以将它改写为 ( ) ( ) ψ(δ) ≡F δR −K ≡ δ δR −0 (K )δ n n × R 是一个 的矩阵，其元素 是矢量 的函数， 为已知矢量。在位移有限元中, k δ δ ij (F )δ R ψ(δ) 0 代表未知的结点位移， 是等效结点力， 为等效结点荷载，方程 表示结点 的平衡方程。 在线弹性有限元中，线性代数方程组 K R 0 δ− ψ(δ) 0 则不行。一般来说，难以求得其精确解， 可以毫无困难地求解，但对非线性方程组 通常采用数值解法，把非线性问题转化为一系列线性问题。为了使这一系列线性解收敛于非 线性解，曾经有过许多方法，但它们都有一定的局限性。某一解法对某一类非线性问题有效， 但对另一类问题可能不合适。因而，根据问题性质正确选用求解方法成为非线性有限元的一 个极重要的问题。本章将介绍有限元分析中常见的各种求解非线性方程组的数值方法。 2.1 迭代法 前面已经提到，目前求解非线性方程组的方法一般为线性化方法。若对总荷载进行线性 化处理，则称为迭代法。 2.1.1 直接迭代法 对非线性方程组 (Kδ)δR −0 (2.1) δ δ 0 K 矩阵 设其初始的近似解为 ，由此确定近似的