固体力学中快速边界元法计算精度的若干研究.pdfVIP

第十三届现代数学和力学学术会议暨钱伟长诞辰100 周年纪念大会文集 固体力学中快速边界元法计算精度的若干研究 姚振汉 （清华大学航天航空学院工程力学系，北京 100084 ） 摘要：作为有限元法的重要补充，边界元法已经发展成为工程与科学分析中的一种精确高 效的数值方法，特别是各种快速算法的引入，快速边界元法克服了常规边界元法难以计算大 规模问题的弱点，为边界元法在复杂工程与科学问题中的应用开辟了广阔的前景。 在固体力学快速边界元法的发展中，常常将复合材料模拟作为一类重要的应用实例，许多 作者都用快速边界元法来计算复合材料的胞元，从而得到了完整的弹性力学解，并由此计算 材料的等效弹性常数，将其与细观力学的理论结果进行比较。作者和合作者在这方面也做过 很多研究，图 1 所示为两个较大规模计算模型的计算结果示意图。 图 1 含200 根随机分布骨形纤维和弯曲纤维胞元的快速多极边界元计算结果界面应力云图 研究表明：由于边界元法高精度的特点，此类计算的精度明显高于有限元计算；与细观力 学结果相比，边界元模拟得到的不仅是等效材料常数而且是各个样本的完整弹性力学解，而 且精确考虑了各个夹杂的实际形状和夹杂之间的相互作用，因此精确计算的结果可以具有比 细观力学解更高的精度；通过和复合材料实测结果的比较，发现夹杂与基体之间的界面条件 不一定是理想的界面条件，因此在计算模型中引入实际的界面条件也是保证模拟结果符合实 际的关键之一。 本文在扼要介绍快速多极边界元法在复合材料模拟中若干应用的基础上，以弹性力学平面 问题采用自适应交叉逼近（ACA ）算法的快速边界元法为例，通过如下系统的数值实验对于 快速边界元法的计算精度做进一步的考核： （1）方板单向拉伸：主要研究对这样的简单问题，当离散成大规模问题时用这种算法能 否得到精确的解。 （2 ）含 5×5 圆孔的方板以位移加载方式双向均匀伸长：这是一个周期重复结构的问题， 单个子结构的边界元解可以作为考核计算精度的依据。 （3 ）含 100 个随机分布圆孔的方板以面力加载方式单向均匀拉伸：这和一个长纤维增强 复合材料截面内弹性特性的胞元计算问题相近。 （4 ）方板的一边固定、在其对边以面力加载方式均匀拉伸：在固定边的两端具有Williams 特征解的应力奇异性，随着单元的细分，将会显现这种应力分布特征。 数值实验表明，ACA 算法通过由边界元法所得方程矩阵子块的自适应低秩近似能显着压 第十三届现代数学和力学学术会议暨钱伟长诞辰100 周年纪念大会文集 缩存储量，并显着提高计算效率。不过低秩近似引起矩阵的失真，会引起附加的计算误差。 对于边界全部给定面力，刚体位移由给定节点位移来约束的问题，这种附加的误差有时会特 别严重。图 2 为第 1 组和第 3 数值实验中两个出现此类问题的算例的变形图，即边界位移图， 可以发现在约束刚体位移的节点处产生了较大的不应有的支承反力，从而会给应力计算结果 带来很大的误差。 图 2 方板和含 100 个随机分布圆孔的方板单向拉伸的有显着误差的变形图 最后对于快速边界元法在固体力学中的应用做简要的展望。边界元法和快速边界元法的研 究是为了给有限元等其他方法提供补充，主要是提供精度高于其他方法的分析方法，因此在 保证计算精度方面一定要认真对待。 关键词：边界元法；快速边界元法；快速多极算法；自适应交叉逼近算法