非线性方程组解法局部弧长法.pdfVIP

第　29　卷　第　1　期 力　　学　　学　　报 V o l. 29,N o. 1 年　 　月　　　 　 　 　 　 1997　 1 A CTA M ECHAN ICA S IN ICA 　Jan. , 1997 非线性方程组的解法: 局部弧长法 段云岭 (清华大学土木系, 北京 100084) 摘要 描述了一个新的非线性方程组的求解方法——局部弧长法. 该方法是在弧长法的基础上发展 起来的适合于材料非线性有限元分析的数值解法. 其约束方程充分利用了结构中破坏区域内的非 线性变形信息, 有效地解决了材料非线性分析中的稳定性与收敛性问题. 数值计算表明, 该方法不 仅适合于求解结构的极限承载能力, 也适合于求解结构达到极限承载参力以后的荷载- 变形的全 过程. 关键词　非线性有限元, 应变软化材料, 裂缝开展, 局部弧长法, 极限承载力, 控制单元 引　言 在对诸如岩石、混凝土和砖石等应变软化材料组成的结构应用有限元方法进行非线性分析 时, 一旦结构接近于极限承载能力或者以后, 其计算分析的稳定性与收敛性相当差以至于无法进 行下去. 虽然已经有许多“自动”或者“自适应”的非线性方程组的数值解法被提出来, 但是到目前 为止尚未有一种方法在材料非线性分析时是有效且可靠的. ( [ 1 ] [ 2 ] [ 3 ] ) 在几何非线性分析中, 弧长法 R ik s 于1979年提出, 后经R amm 和 C risfield 发展 是处 理结构的软化现象或者 现象的有效算法. 该方法同时控制结构的变形与荷载, 使分析 Snap Back 沿着结构的平衡路径进行. 但是, 无论是R amm 的正交平面约束的弧长法还是C risfield 的圆弧 [ 4 ] 约束的弧长法都无法求解材料非线性问题 . 例如当结构接近于极限承载力时, R amm 方法中的 正交平面往往与结构的平衡路径不相交而导致迭代失败; 一旦结构内材料达到极限强度, C risfield方法中的圆弧约束的二次方程常常因没有实根而无法继续求解. 后来, 不少学者[ 5 ] 以及 C risfield 本人[ 6, 7 ] 曾经对弧长法在材料非线性中的应用做了大量的尝试, 如加入线性搜索因子或 者各式各样的加速因子, 但收效甚微. 利用结构内部局部变形信息的尝试是 [ 8 ] 和 [ 9 ] 首先进行的. 应用裂缝 de Bo rst Ro ts de Bo rst 表面的相对滑移信息、 利用裂缝的张开位移分别对具有一条预制裂缝的构件进行了成功的 Ro ts 分析. 可惜这种方法只能用于含有一条裂缝的结构并且结构的承载力是由这条裂缝控制的. 本文所提出的方法是在 的思路上发展起来的. 在荷载- 变形空间, 利用一阶线性 de Bo rst 约束方程强迫迭代过程沿着一个圆弧路径趋向于结构的平衡状态, 且在约束方程中仅使用破坏 区的相对位移以消除来自非破坏区域的变形信息, 从而保证了数值解法的稳定性和收敛性. 1　方法原理 在利用有限单元法对结构进行非线性分析时, 其求解过程可以概述为寻找一组位移向量 U 和荷载向量参考系数 以使得结构内的不平衡节点力向量为零.