桥梁结构几何非线性.ppt

11.5非线性方程组的求解11.5.1求解方法概述用有限元法进行结构非线性分析，其控制方程最终是一组非线性代数方程。非线性代数方程组的求解方法很多，其选择往往与物理问题的性质、特点、非线性程度、对计算结果的要求以及计算机的容量、计算速度等因素有关。这就要求研究者对非线性问题的求解过程以及程序设计有较全面的了解。以下介绍几种常用的求解方法。第80页,共101页，星期六，2024年，5月1)直接求解法直接求解法是基于全量列式的求解过程，应用最多的是直接迭代法，由虚功原理建立的非线性有限元平衡方程为：(11－80)当设定位移向量{}的初值{}后，改进的近似解可由下式得到：(11－81)整个迭代过程可用下式表示：(11－82)当迭代结果满足预定的收敛准则时，就得到了所要求的节点位移向量。第81页,共101页，星期六，2024年，5月图11.12(a)为取{}={0}时单自由度问题的迭代过程取得收敛的示意图。直接迭代法应用简单，运算速度一般也较快，可应用于具有轻微非线性的问题。这一求解过程的成功与否很大程度上取决于对初值位移{}的正确估计。图11.12(b)表示的是直接迭代法迭代过程发散时的情形。为改善收敛性和收敛速度，可以采用将荷载分成若干级的做法。图11.12直接迭代法收敛和发散过程第82页,共101页，星期六，2024年，5月2)增量法增量形式的有限元列式方法具有一个共同的特点：将整个荷载变形过程划分为一连串增量段，每一增量段中结构的荷载反应被近似地线性化。简单增量法将每一级增量荷载下直接求得的状态变量视作结构平衡状态，计算相应的切线刚度阵，进而作下一级荷载计算，并不断累加其位移增量。图11.13描述了简单增量法的求解过程。图11.13简单增量法的求解过程第83页,共101页，星期六，2024年，5月几何非线性问题的有限元分析最初多采用简单增量法进行，虽然这种求解方法对每一级荷载作用时的计算速度较快，但由于每一级荷载作用前结构并未精确地到达平衡位置，所求得的解答会随着增量过程的继续而越来越偏离真实的荷载－变形过程。为了保证计算精度，常常将增量区间划分得相当小。此外，为了评价解的精度，一般要对同一问题在进一步细分增量区间后再次求解，通过两次解的比较判定是否收敛。这样就需要消耗大量的计算时间。作为对这一方法的改进，可将不平衡力作为一种修正荷载并入下一级荷载增量。这就是有一阶自校正的增量法。一阶自校正增量法求解过程的示意图如图11.14所示。一阶自校正增量法具有较高的求解速度，同时也比简单增量法的计算精度高。这一方法在求解非线性问题特别是求解塑性问题时得到广泛的应用。第84页,共101页，星期六，2024年，5月图11.14一阶自校正增量法的收敛过程第85页,共101页，星期六，2024年，5月11.5.2Newton-Raphson法对于式(11－83)给出的单自由度非线性平衡方程(11－83)采用牛顿——下山法，将在点展开成Taylor级数，取线性近似公式：(11－84)求非线性方程式(11－83)的根可按如下公式进行迭代计算：(11－85)第86页,共101页，星期六，2024年，5月单自由度非线性刚度方程一般为：(11－86)或(11－87)式(11－87)代入式(11－84)，有：(11－88)对式(11－87)求导，并注意到得：(11－89)式(11－89)即为体系在处的切线刚度表达式，求相应的迭代公式为：第87页,共101页，星期六，2024年，5月(11－90)为失衡力，式(11－90)即为Newton-Raphson法)求解结构非线性问题的最简单形式，其收敛过程如图11.15(a)所示。法(图11.15(a)第88页,共101页，星期六，2024年，5月对于多自由度体系，同样可以导出相应的迭代公式：(11－91)