第7章-非线性方程组的求解.pptVIP

第7章 7.1 概述 7.2 增量法 增量法的实质是用分段线性的折线去代替非线性曲线。 增量法的特点是：不存在迭代收敛问题、计算误差在计算过程中会累积、不能确定所得解与精确解的误差大小。 根据每级荷载增量求位移增量时所采用的刚度矩阵的不同，增量法又可以分为： 始点刚度法 平均刚度法 中点刚度法 二. 始点刚度法 始点刚度法是在计算第 i 级荷载增量引起的位移增量时，采用第 i-1 级荷载增量末的刚度值来求解，即： 三. 中点刚度法 四. 平均刚度法 7.3 迭代法 7.4 考虑结构负刚度的一些算法 逐步搜索法 对于只要求极值点时适用。 计算步骤： 施加一荷载增量△P，计算发散时，退回到上一级荷载并改用较小步长如： △P /2、 △P /4 、 △P /8 等再进行计算，直至收敛。 2. 增加虚拟弹簧法 在适当地方增加虚拟的、有较大刚度的线弹性弹簧，使施加弹簧后的结构反应不出现下降段，则结构的实际反应等于带弹簧虚拟结构的反应减去线弹性弹簧的相应反应。 3. 位移控制法 分析过程中控制位移增量而不是荷载增量，即由给定的位移求相应的荷载，从而避免控制荷载时的一值多解。 将刚度矩阵重新排列，使得要作为控制加载的位移排到最后一项，同时将原刚度矩阵分块，将有限元方程重新改写而获得新的迭代方程。 3. 位移控制法 分析过程中控制位移增量而不是荷载增量，即由给定的位移求相应的荷载，从而避免控制荷载时的一值多解。 将刚度矩阵重新排列，使得要作为控制加载的位移排到最后一项，同时将原刚度矩阵分块，将有限元方程重新改写而获得新的迭代方程。 4. 其他方法 强制迭代法－用三角分解法分解刚度矩阵并判断结构刚度矩阵对角元中是否已出现负元素，若是，则施加负的荷载增量继续迭代计算。 弧长法－在荷载增量法中，是控制荷载因子的步长 ；在位移增量法中，是控制位移增量的步长；而在弧长法中则同时控制荷载因子和位移增量的步长。 q q 实际梁反应 弹簧反应 带弹簧梁反应 使用时，直接删除本页！ 精品课件，你值得拥有！ 精品课件，你值得拥有！ 使用时，直接删除本页！ 精品课件，你值得拥有！ 精品课件，你值得拥有！ 使用时，直接删除本页！ 精品课件，你值得拥有！ 精品课件，你值得拥有！ * 一. 非线性有限元方程组的特点 用有限元法进行结构非线性分析，最终的控制方程是一组非线性代数方程： 在线弹性结构分析中，为常量。但在非线性结构分析中，它随着节点位移的变化而变化。 二. 非线性有限元方程组常用的求解方法 增量法 始点刚度法； 中点刚度法； 平均刚度法 迭代法 切线刚度法－Newton－Raphson法； 割线刚度法－直接迭代法； 等刚度法－ Modified Newton－Raphson法。 混合法 增量法＋迭代法 一. 增量法求解的基本思路 将荷载分成多个增量后逐级施加； 在每级荷载增量内，结构的刚度 取为常值； 求每级荷载增量引起的位移增量； 累加各级荷载的位移增量，就可以求出任一荷载作用下的总位移。 增量法示意 精确解 增量解 精确解 始点刚度法示意 中点刚度法是在计算第i级荷载增量引起的位移增量时，先求出引起 的位移增量 ，将其累加至位移总量后，计算此位移状态下的刚度矩阵 ，再计算此级荷载增量下的位移增量，即有： 精确解 中点刚度法示意 平均刚度法在计算第i级荷载增量引起的位移增量时，先利用 求得初步位移量 ，将其累加至位移总量后，计算此位移状态下的刚度矩阵 ，再利用平均刚度 求此级荷载增量下的位移增量，即有： 平均刚度法示意 精确解 一. 迭代法求解的基本思路 根据作用于结构上的荷载，作一系列迭代计算 ； 在每次迭代中，结构刚度取为某一数值 ； 迭代后计算总荷载的不平衡力，并把它作为下一次迭代的荷载，以计算出附加的位移增量 ； 如此重复进行，直到计算结果达到预定精度即终止迭代过程 。 二. 迭代法的种类 割线刚度法（直接迭代法） ； 切线刚度法－Newton－Raphson 法； 等刚度迭代法－Modified Newton－Raphson 法。 1. 割线刚度法 割线刚度法（直接迭代法）是在迭代过程中采用割线刚度求出总荷载引起的总位移，再根据总位移调整割线刚度的一种迭代法。 割线刚度迭代法示意 K0 K1 K2 2. 切线刚度法-Newton-Raphson 法 切线刚度法是在迭代过程中采用变化的切线刚度求出总荷载引起的总位移，再根据总位移调整切线刚度的一种迭