伽辽金边界元法在平面问题极限分析中的应用.pdf

文章编号:1006-7329(2000)06-0112一05 伽辽金边界元法在平面问题极限分析中的应用’ 赵亚楠，张晓峰，刘应华，岑章志 (清华大学 工程力学系，北京 100084) 摘妥:基于板限分析的下限定理，利用伽辽金边界元方法对二维结构进行离嵌，建立了结 构桩限(下限)分析的计算格式。利用弹皇性咐f葬法中同-4-k步上不同迭代少的应力 且作为基夫全构过了自平衡应力场，然后通过复合形法立接术解该林钱性况划问趁，得到 了二维结构权限或待下限来于。计葬结果表明了本文来用的方法的有效性。 关.词:知辽全边界元法;权限分衍，基失t，复合形法 中.分维号.095i.a 文欲坛识码:A 在结构分析和安全评估中，极限分析方法得到了广泛的应用.作为求解弹塑性间题的简化方 法，极限分析一般不盆要事先知道裁荷变化的历史，应用更简便.理论上这种方法可以遵免费时的 弹塑性增t计算，适用于对工程结构的极限承载能力进行求解，但在实际的数值计算时却遇到了很 大的困难。由于结构极限分析的教学规划格式含有多个变t和约束，且一般来说为非线性，求解的 规棋往往很大。在将非线性屈服条件线性化的过程中，为保证一定的精度，需要引入一定数量的约 束和变f。而约束和变tI的数目越多，数学规划的规模将呈指数上升趋势，变成一大规棋的数学规 划问压，这峨是所谓的维致障碍。目前很多研究都国绕粉发展极限分析方法而展开，其中大量的工 作二盆是为了克服维教障碍问题。 应用边界元方法进行极限和安定分析，早在1983年就由Maierc]等提出了采用常规边界元进行 极限和安定分析的方案，1992年Panaecao」给出应用伽辽金边界元进行极限和安定分析的方案。由 于没有克服效值计算上的困难，上述方案并未给出具体的计算方法和有关的计算结果。伽辽金边界 元方法，由于采用解析积分方案和对称的系数矩阵使该方法有很好的计算精度和计算效率。特别是 由内点应力积分公式直接得到的应力比有限元方法由位移计算的应力有较高的精度，便于构造虚 拟弹性应力场和自平衡应力场，采用伽辽金边界元方法进行极限下限分析是很有利的。本文基于极 限分析的Melan定理，采用伽辽金边界元方法对结构进行离欲，并利用弹塑性应力模拟法构造了自 平衡应力场，将间厄最终归结为一个数学规划何题并直接求解此非线性规划何题。计算结果表明了 该方法的有效性。 1 基本理论 1.1扭爪分析的下限定理 极限分析的下限定理1[31可以叙述为:“一个满足平衡条件井且到处都不破坏材料屈服条件的应 力场，是一个静力容许应力场，与静力容许应力场对应的外载荷是极限截荷的下限;最高的下限就 是极限载荷.即其狱学表达格式为: max,B s.t.KWI(二)+A!W) 。 Vx 〔n ((l1ba)) .收摘日翻,2000-0s-01 若金项目.份华大举班翻研究签金，.家自艘科学班金 作者简介:赵蓝摘。9”一〕.女，天津人，砚士生.主要从事计算力掌研究。 一-一~一~--~~~~，~--气口巴盯，，尸下，，-一 ii一居丁 三E 呈6A -一-kA*T-:lt-iTli-t*)c.k丝fInp7AAIX#*toAm一_ 一-」丝- P;;.;=o Vx任一 1c) Prini二0 VxE一 id) 其中目标函数9为待求极限载荷乘子，峙”为给定基准载荷作用下的虚拟弹性应力场，P,;(x)为自 平衡应力场，抓·〕为屈服函数。约束条件((16)式为屈服条件，(lc)式与((1d)式为自平衡应力场在域 内和边界上所必须满足的条件，而自平衡应力场是主要的自变量。 1.2弹盆性伽辽全对称型边界元法 翻辽金对称型弹塑性边界元法4〔〕的建立是从贝蒂功互等定理出发，以实际结构的受力状态为 实状态，结构在 “广义载荷..(任惫