不同介质组合结构虚边界元分析研究.pdfVIP

强度与环境2000增刊 1 Shucldu+eEnvhumientEneermg,Su碑dement2000川 不同介质的组合结构虚边界元分析 许 强 戴月辉 (上海铁道大学土木工程系，上海，加0331) 文摘 对不同介质的组合结构，关于其精合面附近的分析计算的研究，既具有理论意义，又具 有应用价值。本文依据虚边界元法[[,.21P出了不同介质组合结构弹性力学间题的求解方程。通 过本文的数值算例结果与解析解和有限元解的比较，证明了本文算法的有效诊.- 关.词 性力学虚边界元最小二乘阵 I 引言 在工程实际中，随着生产的需要和材料科学的发展，常使用两种或两种以上不同性质的 材料通过契合、胶合、焊接等手段组成的结构。对这些不同介质结合物体的结合面附近域进 行准确的分析计算是值得研究的间题，且具有应用价值。如不同材料的胶合结构，其结合面 上的接触应力不应超过材料间的允许粘接强度，否则将会引起裂缝，导致应力和位移的重分 布。对于这些问题的计算，倘若采用边界元直接法的韧合法求解，则由于解的奇异性使得结 合面附近区域内位移和应力计算精度很低;倘若采用文[3」的修正边界元法求解，则由于 Hentenyi，基本解表达式异常繁琐和该方法应用的局限性(对三维问题，仅适用于结合面为 平面;二维问题 ，仅适合于结合面为直线)，使其应用不易推广，另外，边“界层效应”也使该 法的有效性降低;再者，倘若采用常规有限元法求解，则由于位移是基本未知量，位移在单元 间是连续的，而应力在单元间不是连续的。为了获得较好的解，本文依据虚边界元法思 想1,27，从一般角度推导出不同介质组合结构弹性力学间题的求解方程。本文方法可避免 上述方法的不足。通过数值结果的比较，证明了本文算法的有效性。 2 不同介质弹性力学问题的虚边界元法 为了叙述方便，现以两个不同介质的组合物体为研究对象，当然，本文的思想同样适用 于。个不同介质的组合物体。如图1所示，设D1,l22分别为各自介质所占的弹性区域，且211 nn2二曰(空集);设弹性模量和泊松比为E:,v1;s1Us}为。:的边界，S2ns,为。2的边 界，s,为两种不同弹性介质的公共边界(即结合面);设sI.,s0;为对应介质物休的已知面力 许强，男,1951年生，安徽安庆，博士，副教授,(200331〕上海铁道大学土木工程系。 收摘Fl期初摘99-01-I8,修改稿99-03-26 110 和已知位移边界，则对n,有s;二S,Us0.，对rat有s:二s2Us}2，且均有s,:fls0.=必(空 集);设0,,n:的虚边界分别为S,和S2，其上 作用有未知待求的分布体积力(简称虚体力)， 分别记为。cu(s)(E任s,)和。(2)(F)(E任 s12). 根据文 【1〕的思想，可将习嵌‘入一个具有 空洞为日的‘无限域中其〔材料与r2的‘相同)，在 无限域中作一有限虚边界S;(s,DS;Us,),s; 包含的区域记为niA :)口)‘。每一。所‘对应 图1 两种不同介质的组合结构 的无限域在分析计算时只与口有‘关。 由弹性力学理论可知，任一给定的弹性力 学问题的解是唯一存在的，且具有连续性;而Kelvin解在无限域内满足泛定方程;由文 【1] 中的式(11)一(14)可知，峋(:)、。，(二).ui(x).p;(二)在伏内满足连续性条件。于是，可以通 过满足原问题的边界条件和结合面s。上的应力、位移连续性条件来求了0(E),了2)(F)。在 具体实施中，对已知边界SS2采用虚边界元数值积分法I[]，对结合面凡采用虚边界元配点 法1[21;再由最小二乘法，可建立关于虚体力函数g(,)(E),q(2)($)的方差泛函如下 I[qu，q2):二f，a[UP(x)一u)(x)]2dS·fS，R[pu{(x)一，}1)(x)]2dS ·fvia[u2;)(x)一I2“()(X)12dS·丁fs、,}，P，lp，l2〔()(x)一、2()(X)]2dS +艺a[u()(xj)一。,(2)(.j)]dS+艺PI[P尸,(])(}j)一P尸,(2)(.j)]2