二维弹塑性问题的复变函数边界元法.pdfVIP

第 19卷第2期 同 济 大 学 学 报 Vo1．19№ 2 1口口1年 6月 JOURNAL 0F TONGJIUNIVERS1TY Jlln．1991 二维弹塑性问题的复变函数边界元法 唐寿高 (工程力学系) 摘 要 ’ 本文将复变函数连界元发展用于解二维弹塑性力学 丁口题． 给出基于初应力荷载增量法的增量形 式的弹塑性问题复变教边界积分方程殛所有基本关系式．编制了相应的计算机程序并培出算例。与 通常边界元解此类问题相比，本身法具有 公式统一化． 程序简单且多功能．计算机计算时阉短且精 发高等优点．计算表明，复变函数边界元特别适合用于需迭R计算的弹塑性等非线性问题。 关键词 ：弹塑性，复变函数，边界积分方程。 一 、 引 言 自从 Swedlow和Cruse[”等人于 1971年将边界元发展用于解弹塑性力学问题 以来 · 许多学者对这一领域作了富有成效的研究，如文献r2l所述。就作者看来，其中具有代表佳之 一 的是 Baaer~ee等人的工作。他于 1980年首次提出了基于初应力法的直接边界元法用于二 维弹塑性分析 ”【，所用初应力迭代步骤与 Zienkiewicz~‘在有限元中所用方法相类似。 后来 Baaerjee和Raveendrat 引进加速初应力迭代技术和二次单元模式改进 了原有的成果。 最 近，Henry轺Banerjee~ 通过数学处理得到的特别积分， 避免了由于考虑非线性影响而要 对区域内部计算的面积分或体积分。然而，上述方法仍有待于完善。如由于奇异积分，使解 在边界及附近塑性区内难以收敛，考虑非线性因素后程序编制繁复，且无通用性I必须对所 有边界划分单元，这样若采用迭代法求解，需耗费较多的计算机时间等等。 本文的尝试为解决上述问题提供了一个可供选择的途径。文中仍采用惯甩的初应力法 · 但由于核函数用复位势基本解来表示，所得增量形式的边界积分方程具有一般意义，所编程 序可处理任意多个问题的基本解。针对具体问题采用相应的复位势基本解，这样可省去许多 本来孺要划分的边界单元，使方程组阶数大为降低，有种子反复迭代计算。再则解在不需分 割单元的边界及近旁弹塑性区域 内可给出令人满意的结果。还有本方法公式简单，程序编制 容易且具统一性多功能，均比通常的边界元法优越。 需要指出的是早~_1982年MukherjeeE。】就已将复变数用于边界元间接法分析 蠕变和断裂 问题。他也利用MttskhelishviliE”弹性力学复变数理论 ，但由于是直接从非齐次双调和方程 变换得到边界积分方程，因而方程的意义及核函数等与文献[8]中的不同。此外，Mukherjee 的方法在实际运算时各主要量都是事先转变成实数米进行的 然而利用复变数的主要 目的之 率定收到 日期：1990年 3月8日 同 济 大 学 学 报 筇 10卷 一 却是相同的，即在高度应力集中的边界上免去单元划分 以便作精确的数值分析。Hromadka 辱 ’将复变数用于边界元 目的是对柯西积分进行离散，只解决了直接 由柯西积分给出解的较 为简单的标量场问题，而弹性力学等较为复杂的向量场问题的解是不能直接由柯西积分公式 给出的。在边界元中使用复变函数基本解的概念首先是由黄弘量川提出，并对二维弹性问慝 作了成功的尝试。本文表明，复变函数边界元不仅适用于线弹性分析，而且特别适台于需选 代计算的弹塑性辱非线性分析。 二、复变函数边界积分方程及其基本公式 基于加权残值理论的弹塑性力学平面问题增量形式的复变函数边界积分方程为 c。I(’)；()=frfUii(，口)．(口)一Pj(，口)(q)3dF(q) + ( ) · )d口(口)