弹性力学平面问题大单元边界元法.pdfVIP

弹性力学平面问题大单元边界元法 赵曼扬 (石油大学(北京)机电系，北京，102200) 摘要本文提出大单元(Macro--element)直接边 界元方法求解弹性力学平面问题。在应用边界元方法 2平面问题的边界元法 解材料不均匀或几何彤状豆栽荷情况复杂的弹性力 2．1弹性力学平面问题厦其边界条件 学平面问题时，先将求解域划分为几何形状同一的大 单元，而在太单元的边界划分边界单元，用边界元法 弹性力学平衡方程 求解。对同一材料性质的各丸单元只须计算一个大单 ％．，+bk一0 域n内 元的矩阵元素。建立的总刚度矩阵是分块带状稀疏 位移边界条件： 的。数值算侧表明这种方法的精度满意。 q=uj边界n 关键词大单元边界元法平面问题 应力边界条件： 1引言 p，=d，lHk=pJ 边界巳 用边界元法解弹性力学平面问题时，当材料不均 2．2弹性力学平面问题的边界元方法 匀时必须将求解域按相同材料性质翔j分为子域，然后 在各子域的边界上定义边界元。这是因为只有在材料 以基本解u《和Pt作为杈函数，用加权余量法可 性质相同的区域才船使用同一个基本解。另外当求解 以建立如下方程”。 域的形状比较复杂或者载荷比较复杂时也往往采用 这样的子域法。边界元方法的子域法的讨论见[1，z]。 J。(at¨+bk)“dn2J 4A一只)“‘dr— 一般情况下子域法对各子域的形状和大小以及边界 r】(m一“-)户AdF (4) 单元的划分并没有特别要求。本文提出的大单元是将 通过对(4)进行分步积分可得到积分方程： 材料性质相同的区域划分为几何形状完全相同的太 单元并在各大单元的边界上划分相同的边界单元。这 样在建立方程时，对同一材料性质的各大单元只须计 式中“=a．tn．，而ul表示位于i的点在1方向上的位 算一个大单元的矩阵元素。最后建立起一个分块稀疏 移。对二维问题各下标取1，2时方程(1)至(s)。各向 的方程矩阵而不是满阵。本文讨论了此方法的精度。 同性平面问题的基本解为： 当各子域具有相近变量数时其精度要好于各子域的 变量数有很大差别口j。而大单元法中各大单元的变量 数正是完全相同。数值算例表明大单元法可以取得较 (6) 好的结果。 (7) 靠—一石拇{毒l“一朋)△”+2毫毫j一(1—2小蠹m+毫帕) 式中△u是Kroneckerz▲。将基本解引人方程(5)，并通r，c‘=o．5}当i在外，c’一o。 过将源点移到边界上，得到边界积分方程 方程(8)或者(9)是边界元法解平面问题的基本 。·。：+f。。。(靠打= 方程·通过适当的方式将其离散并处理了奇异积分可 L丸。二：j+L仉。；。。 。。，以得到一个线性方程气。；， 。，。， 对不计体积力的情予…方掣为r ， 塞饕甚嚣篇篇豪嚣需篱字g宝霎耄，-…o…7●，●……7……l……J…1…Ⅱ