边界元法在颗粒增强复合材料模拟中的应用.pdf

.386 · 边界元法在颗粒增强复合材料模拟中的应用 王海涛，姚振汉 (清华大学z程力学系，北京 100084) 摘要:在本文中，边界元法用于二维和三维颖粒增强复合材料的数值模拟.快速多极算法作为边界元法的求解算法，从而使 边界元法能够对含有大f随机分布颗粒的复合材料进行大规模模拟.广义极小残差法((GMRES)用作迭代运算的求解器。数 值算例表明，边界元法能够有效进行颖粒增强复合材料的等效模量和应力分析. 关锐询:边界元;快速多极:颗粒;复合材料 APPLICATIONOFBOUNDARYELEMENTMETHODFORSIMULATION OFPARTICLE-REINFORCEDCOMPOSITEMATERIALS WANGHaitao,YAOZhenhan (DepartmentofEngineeringMechanics,TsinghuaUniversity,Beijing100084,China) Abstract:Thispaperpresentstheapplicationofboundaryelementmethod仍EM)forsimulationof2Dand3D particle-reinforcedcompositematerials.FastmultipolemethodisusedasafastsolverforBEM,makingBEM applicableforlargescalesimulationofcompositeswithalargenumberofrandomlydistributedparticles.General MinimumResidualMethod(GMRES)isusedforiteration.NumericalresultsshowthatBEM isefficientfor effectivemodulusandstressanalysesofcompositescontainingparticles. Keywords:boundaryelement;fastmultipole;particle;compositematerial 1引言 边界元法是处理连续介质问题的有效的数值模拟方法之一。与其他数值方法相比，边界元法只需要在 边界离散单元，降低了问题的维数和规模，同时精度很高，因此边界元法特别适合处理含有很多颗粒的复 合材料结构III.但是常规的边界元法不适合处理大规模问题，因为边界元法形成的系数矩阵是非对称的满 阵。常规求解方法往往需要ow3)的计算量和ow2)的存储量，其中N是问题的自由度。而快速多极算法 的出现，使二者的规模都下降到了线性的O(N)量级。 快速多极算法由Rokhlin2(1在1985年提出，用于快速求解势问题。由于快速多极算法在求解大规模问 题上的潜力，将快速多极算法用于边界元法的求解得到了不同领域学者的关注和研究。文献[3,41是)1,种 O(NlogN)和O(N)量级的快速多极边界元法在弹性力学问题中的应用。最近，新型快速多极算法1[51的出现， 又大幅度提高了边界元法的求解效率。 本文中，采用新型快速多极加速的边界元法与相似子域法[[1]结合，用于含有大量随机分布的圆形颗粒 的二维复合材料的等效模量的模拟计算;采用GMRES迭代求解的常规边界元法与相似子域法结合，用1一几 基金顶目:国家自然科学鉴金资助项目.教，部博士点专顶基金 作者简介:王海脚，男 山东人，协士生，主拼从事计抹力学及计体工扭研究; 姚报汉.男 江苏人，教长，博士生辱师，主要从事边异元法、高性能有皿元挂和无网格法的研究 ·387 · 含有球形颗粒的三维复合材料的应力分析。 2用干二维弹性力学的快速多极边界元法简介 用于描述弹性力学问题的微分方程可以通过Kelvin解和Betti功互等定理等效成积分方程。 厂面以二 维问题为例，忽略体积力的二维弹性力学的边界积分方程是: cO,(x)up(x)+工T.,,(x,Y)u,,Y()drY()二工U}(x,Y)t,,Y()drY() (1) 其中x和