用边界元方法和复合形法求解三维结构下限安定载荷.pdfVIP

第20卷第3期 I。 程 )3 哗 VOl20No．3 2003年6月 ENGINEERJNGMECl4ANICS Jun2003 文章编号：1000—4750(2003)03-0013—07 用边界元方法和复合形法求解 三维结构的下限安定载荷 张晓峰，刘应华，岑章志 (清华人学工程力学系，北京100084) 摘要：基于安定分析的静力定理，建立r用常规边界元方法进行二维理想弹颦性结构安定分析的整套求解算 法。F限安定分析所需的弹性应力场直接由边界元方法求出，所需的白平衡应力场由一组带有待定系数的自平 衡廊力场基矢量的线性组合进行模拟，这些自平衡应力场基矢量通过边界元弹塑性迭代计算获取。安定分析问 题最终被归结为系列未知变量较少的非线性数学规划子问题并通过复合形法直接求解。计算结果表明了算法 的和效性。 关键词：边界元法：安定分析：自平衡成力场：非线性规划；复合形法 125 中图分类号：TB 文献标识码：A 1 引言 难。由于结构安定分析最终归结为求解一个数学规 划的极值问题，而这个数学规划格式含有多个变量 结构安定分析是塑性力学最有实用意义的分 和约束，且一般来说为非线性，求解的规模往往很 支之一，其宗旨是确定各类工程结构的安定载荷， 大。同时，在将非线性屈服条件线性化的过程中， 为工程设计和安全评估提供必要的理论依据。在循 为保证一定的精度，需要引入一定数量的约束和变 环加载情况下，安定载荷是控制工程结构塑性损伤 量。而约束和变量的数目越多，数学规划的规模将 累积的重要指标。当循环载荷的幅值小于安定载荷 呈指数上升趋势，变成一大规模的数学规划问题， 时，结构虽然可能在加载初期出现局部变形，但经 这就是所谓的维数障碍。如何建立一种有效的求解 过一定次数的载荷循环后塑性损伤不再积累，结构 算法，在保证计算精度的同时，又能大大降低求解 将呈现弹性行为。反之，当载荷幅值大于安定载荷 的计算规模，是安定分析方法能够真正应用于工程 时，每个载荷循环都将产生塑性变形，结构将进入 实际问题的关键所在。 塑性损伤不断积累的逐步破坏过程，最终导致循环 迄今为止，安定分析的数值计算方法主要是基 塑性(低周疲劳)或渐增塑性(棘轮效应)破坏。安定载于应用范围已经相当广泛的有限元方法。采用边界 荷已经成为近代工程强度设计规范中塑性失效准 元方法进行极限和安定分析，早在1983年就由 则的重要判据。 Maiertll等提出了采用常规边界元方法进行极限和 在实际的工程问题中，往往无法事先知道循环 安定分析的方案， 载荷随时间变化的历史，而只能知道载荷的变化范 Galerkin边界元进行极限和安定分析的方案。由于 围，因此用弹塑性增量方法进行分析是非常困难 没有克服数值计算上的困难，上述方案并未给出具 的。而进行安定分析可以避开加载历史，直接求解 体的计算方法和有关的计算结果。边界元方法，由 安定载荷，这对于工程结构的设计和分析是十分有 利的。但在实际的数值计算时却遇到了很大的困