非线性有限元课程.doc

1 非线性问题的求解方法 无论是哪一类非线性问题，经过有限元离散后,它们都归结为求解一个非线性代数方程组： …… 其中：,……是未知量，,……是,……的非线性函数，现引用矢量记号 ? 上述方程组可表示为： 还可以将它改写为： 是一个的矩阵，其元素是矢量的函数，为已知矢量。在位移有限元中，代表未知的结点位移，是等效结点力，为等效结点荷载，方程表示结点的平衡方程。 在线弹性有限元中，线性代数方程组可以毫无困难地求解，但对非线性方程组则不行。一般来说，难以求得其精确解，通常采用数值解法，把非线性问题转化为一系列线性问题。为了使这一系列线性解收敛于非线性解，曾经有过许多方法，但这些解法都有一定的局限性。某一解法对某一类非线性问题有效，但对另一类问题可能不合适。因而，根据问题性质正确选用求解方法成为非线性有限元的一个极重要的问题。 1.1 直接迭代法 目前求解非线性方程组的方法一般为线性化方法。若对总荷载进行线性化处理，则称为迭代法。 对非线性方程组： 1-1 设其初始的近似解为，由此确定近似的矩阵，根据式1-1可得出改进的近似解。重复这一过程，以第次近似解求出第次近似解的迭代公式为： 1-2 直到Δ变得充