计算固体计算力学 - 第二章 非线性方程组的解法综述.pptVIP

博士研究生课程 计算固体力学 课程编号：017090 王生楠，谢伟 西北工业大学 航空学院; 非线性问题可分为三类：材料非线性、几何非线性和边界非线性。我们只讨论前两类问题。; 上述方程组可简单的用一个矢量方程表示 Ψ(a)=0，a为待求的未知量。; 2.1 直接迭代法 2.2 Newton-Raphson法(简称N-R法) 2.3 改进的Newton-Raphson法(简称M-N-R法) 2.4 增量法;2.1 直接迭代法;当设范数为;，;，;收敛性;2.2 Newton法;有限元分析中Jacobi矩阵代表切线刚度矩阵; 注意到在迭代的每一步，得到的近似解一般不会严格满足方程组;;;;;其中μn的作用是改变切线矩阵KT的主对角元素，使奇异性或病态得到改善。;;; 简化的Newton法求解过程如图：;;加速技术：修正牛顿法的过量修正技术－搜索办法;加速技术：Aitken加速技术; 如果吸取Newton法收敛快与简化Newton法每步计算工作量小的优点，???把M次简化Newton步组成一次Newton步。也就是说，在使用简化Newton法时，每经过M次迭代后，重新计算一次Jacobi矩阵KT; 为了避免Newton法每步计算偏导数，可以用函数值代替导数来建立新的计算方法，拟Newton法就是一种有效的算法。; 一个变量的非线性方程的拟Newton迭代过程如下图所示，Kn是割线的斜率。因此，N维方程组的拟Newton法可以看做是割线法的一种推广。; 使用某种算法的计算效率，除了与收敛速度有 关外，还与每一步迭代所花费的计算量有关。关于 每步的计算量，牛顿法最大，而修正牛顿法最小。 因此在实际问题的计算中判断使用哪种方法效率较 高，往往需要进行数值实验。总的看来，不同的算 法可能适用于不同的问题。选用哪种算法，与所研 究问题的性质 (例如，对线性的偏离程度)、规模 (离散的自由度总数) 以及容许误差等因素有关。;; 设“荷载”(R)在任一增量步的值为λ(R)，λ 为荷载增量因子，(R)为标准荷载列阵,则非线性方程Ψ(a)= 0可写为; 设荷载增量因子λ分别取如下值; 在一个自变量的情况下，求解非线性方程组的过程如下图所示。如果Δλm足够小，则认为所得解即为方程组的合理的近似解。;自修正算法，在增量法每一增量步进行自修正的迭代计算。其m增量步n次迭代的计算公式为; 所谓混合法是指，在增量法每一增量步进行自修正的迭代计算。其m增量步n次迭代的计算公式为;;;;;;;;本章结束！