非线性方程组解法.pptVIP

**非线性方程组解法材料非线性几何非线性边界非线性以上三类问题的有限元问题都将得到待解的非线性方程组。非线性方程组一般都采用线性化方法。一般线弹性有限元的总体刚度矩阵为[K]是一个常数矩阵，上述方程是一个线性方程组。一般非线弹性或弹塑性有限元的总体刚度矩阵为[K(U)]不是一个常数矩阵，上述方程组成为非线性方程组。可记为简例材料非线性图示简单拉伸问题中如果，本构方程是非线性弹性的设杆端位移为u1，则这里我们得到的是非线性方程。εσ非线性方程解法迭代法（总载荷法）Newton－Raphson法（切线法）修正Newton－Raphson法（初始切线法）拟Newton－Raphson法(割线法）增量法（逐步法）Euler－Cauchy法混合法(逐步迭代法）Euler－Newton法迭代法（总载荷法）Newton－Raphson法（切线法）在{U}n处展开，并用un表示{U}n设则得到逐步迭代就可得到解。求解的几何意义在一维的情况下是明确的，[K(u)]T表示该处的切线，{F}-{R(un)}称为不平衡力，不平衡力逐步减小表示解是收敛的，反之，解不收敛。Puiui+1ui+2R(ui)K(ui)本方法每步均需求切线的刚度矩阵，计算量大，收敛速度快。如果每次都使用初始刚度矩阵，则称为修正Newton－Raphson法（初始切线法）。计算量小，但速度慢。Puiui+1ui+2R(ui)K(u0)增量法（逐步法）用载荷增量进行线性化处理。Euler-Cauchy法如图所示，采用下式步骤进行求解:：Puiui+1ui+2PiK(ui)Pi+1在每一增量步中都按线性弹性问题求解，虽然也将物体变得刚硬一些，会出现“漂移”，但当载荷增量取得甚小时，分段的线性解就接近非线性解。Puiui+1ui+2PiK(ui)Pi+1混合法（逐步迭代法）用混合使用增量法和迭代法求解。一般采用在一个增量步中采用迭代方法。Euler-Newton法如用{P}i，{P}i+1表示第i步和第i+1步{P}的初值和终值，{U}i，{U}i+1表示第i步和第i+1步{U}的初值和终值，在该增量步中采用下式迭代Pi+1PiKi0+1迭代收敛判据在使用迭代法和混合法时，必须给出收敛判据，否则无法中止迭代。当判据给得不恰当时，就不太精确或计算时间过长，甚至计算失败。