基于有限元与宽频快速多极边界元的二维流固耦合声场分析

采用有限元/快速多极边界元法进行水下弹性结构的辐射和散射声场分析。Burton-Miller法用于解决传统单Helmholtz边界积分方程在求解外边界值问题时出现的非唯一解的问题。该文采用GMRES和快速多极算法加速求解系统方程。针对传统快速算法在高频处效率低和对角式快速算法在低频处不稳定这一问题,该文通过结合这两种快速算法形成宽频快速算法来克服。同时该文通过观察不同参数条件设置下,宽频快速多极法得到的数值结果在计算精度和计算时间上的变化,得到最优的参数组合值。最后通过数值算例验证该文算法的正确性和有效性。     

大型浮体水弹性作用的频域分析

对大型浮体水弹性响应的频域计算方法做了综述和介绍。分别介绍了干模态法、湿模态法和直接计算方法,对于干模态法的五种结构弹性模态函数做了介绍,研究了五种模态函数下计算结果的收敛速度。对于水动力分析方法的计算量和存储量做了分析,介绍了降低计算量和存储量的一些新计算方法,实现了一种柱坐标下的多极子展开高阶边界元方法,积分方程中的固角系数和柯西主值积分均采用直接方法计算,应用该方法计算了波浪与大型弹性浮体的相互作用问题,计算结果与实验值得到很好的吻合。     

基于快速多极子边界元法的齿轮箱声场分析

齿轮箱是广泛应用的工程机械零部件,准确地模拟其辐射声场对后续的降噪优化设计有着重要作用。边界元方法非常适合分析此类无限域下的声辐射问题。但传统边界元方法有着计算效率低、内存占用高的缺点。该研究发展了宽频的快速多极子边界元方法,并运用该方法计算了齿轮箱在特定频率下的场点声压以及辐射声场。通过对比商用软件的分析结果,验证了所提快速边界元方法的准确性。此外,运用多核并行计算方法,对计算量较大的扫频分析进行加速计算,最终快速、准确地获取了齿轮箱辐射声场的扫频结果。     

激光干涉仪反射镜三维温度场的快速多极边界元分析

温差造成光刻机激光干涉仪反射镜热变形,从而影响光刻机的精度。该文将基于二次单元的快速边界元法用于激光干涉仪反射镜的大规模温度场模拟。不连续单元的引入可以有效处理角点问题;新型快速多极算法用于边界元法的加速求解。建立统一的二次单元多极展开格式以处理混合边界。数值算例分析了快速多极边界元法的计算精度和效率,并和常规算法比较;使用该算法对激光干涉仪反射镜进行了大规模温度场计算,并和有限元法比较。结果表明：基于二次单元的快速多极边界元法可以高精度求解大规模三维传热问题。     

求解具有弱界面的复合材料有效性能的近似模型与边界元法

基于中间层模型, 建立了用于描述具有弱界面的多夹杂复合材料的三层嵌入式模型。首先, 给出了计算多夹杂复合材料有效性能的公式, 随后将三层嵌入式模型分成两种体系得到了相应的细观力学近似方法。此外, 针对中间层模型进一步给出了快速多极边界元的基本列式。最后, 通过算例验证所提出的两种方法的正确性及有效性, 并进一步考察了界面性能对复合材料整体宏观性能的影响。比较结果发现, 细观力学近似模型形式简单, 易于编程计算, 适用于快速预测具有弱界面的复合材料有效弹性性能。而数值算法能够快速有效直观的考察复合材料内部受弱界面影响后的应力分布。算例分析结果表明, 当中间层厚度较大而弹性模量较低时, 会使得夹杂不再具有增强效果。并且, 随着中间层厚度的增加, 基体内承担的最大应力也随之增加, 从而容易导致复合材料在界面处产生破坏。     

An Adaptive Fast Multipole Boundary Element Method for Three-dimensional Potential Problems

An adaptive fast multipole boundary element method (FMBEM) for general three-dimensional (3-D) potential problems is presented in this paper. This adaptive FMBEM uses an adaptive tree structure that can balance the multipole to local translations (M2L) and the direct evaluations of the near-field integrals, and thus can reduce the number of the more costly direct evaluations. Furthermore, the coefficients used in the preconditioner for the iterative solver (GMRES) are stored and used repeatedly in the direct evaluations of the near-field contributions. In this way, the computational efficiency of the adaptive FMBEM is improved significantly. The adaptive FMBEM can be applied to both the original FMBEM formulation and the new FMBEM with diagonal translations. Several numerical examples are presented to demonstrate the efficiency and accuracy of the adaptive FMBEM for studying large-scale 3-D potential problems. The adaptive FMBEM is found to be about 50% faster than the non-adaptive version of the new FMBEM in solving the model (with 558,000 elements) for porous materials studied in this paper. The computational efficiencies and accuracies of the FMBEM as compared with the finite element method (FEM) are also studied using a heat-sink model. It is found that the adaptive FMBEM is especially advantageous in modeling problems with complicated domains for which free meshes with much more finite elements would be needed with the FEM.     

Fast single domain–subdomain BEM algorithm for 3D incompressible fluid flow and heat transfer

In this paper acceleration and computer memory reduction of an algorithm for the simulation of laminar viscous flows and heat transfer is presented. The algorithm solves the velocity–vorticity formulation of the incompressible Navier–Stokes equations in 3D. It is based on a combination of a subdomain boundary element method (BEM) and single domain BEM. The CPU time and storage requirements of the single domain BEM are reduced by implementing a fast multipole expansion method. The Laplace fundamental solution, which is used as a special weighting function in BEM, is expanded in terms of spherical harmonics. The computational domain and its boundary are recursively cut up forming a tree of clusters of boundary elements and domain cells. Data sparse representation is used in parts of the matrix, which correspond to boundary‐domain clusters pairs that are admissible for expansion. Significant reduction of the complexity is achieved. The paper presents results of testing of the multipole expansion algorithm by exploring its effect on the accuracy of the solution and its influence on the non‐linear convergence properties of the solver. Two 3D benchmark numerical examples are used: the lid‐driven cavity and the onset of natural convection in a differentially heated enclosure. Copyright © 2008 John Wiley & Sons, Ltd.     

A fast multipole boundary element method for 2D multi-domain elastostatic problems based on a dual BIE formulation

A new fast multipole formulation for the hypersingular BIE (HBIE) for 2D elasticity is presented in this paper based on a complex-variable representation of the kernels, similar to the formulation developed earlier for the conventional BIE (CBIE). A dual BIE formulation using a linear combination of the developed CBIE and HBIE is applied to analyze multi-domain problems with thin inclusions or open cracks. Two pre-conditioners for the fast multipole boundary element method (BEM) are devised and their effectiveness and efficiencies in solving large-scale problems are discussed. Several numerical examples are presented to study the accuracy and efficiency of the developed fast multipole BEM using the dual BIE formulation. The numerical results clearly demonstrate the potentials of the fast multipole BEM for solving large-scale 2D multi-domain elasticity problems. The method can be applied to study composite materials, functionally-graded materials, and micro-electro-mechanical-systems with coupled fields, all of which often involve thin shapes or thin inclusions.     

快速多极边界元法在大规模传热分析中的应用

该文将快速多极边界元法用于三维稳态传热问题的大规模数值计算。多极展开的引入使得该算法能够在单台个人电脑上完成30万自由度以上的传热边界元分析。统一展开的基本解能够处理混合边界。广义极小残差法作为快速多极边界元法的迭代求解器,数值算例分析了快速多极边界元法的计算效率。结果表明:快速多极边界元法的求解效率与常规算法相比有数量级的提高;在模拟复杂结构大规模传热问题上将具有良好的应用前景。