固体力学中快速多极边界元法研究进展

和快速多极方法相结合，使边界元法处理大规模工程与科学问题变得十分有效,首先概括介绍了快速多极边界元法，接着介绍其精度和效率的验证以及与常规边界元法的比较，并给出在微机机群上的快速多极边界元并行算法,给出了快速多极边界元法的一些应用，其中包括：复合材料的二维、三维模拟，含大量裂纹的二维弹性固体及其疲劳裂纹扩展的模拟．此外还介绍了用于弹塑性问题的快速多极边界元新方法．     

基于积分核级数展开的多极边界元法及其截断误差分析

对于二维Helmholtz方程问题,本文提出一种基于积分核级数展开的多极边界元方法,推导证明了二维Helmholtz方程的多极展开定理,给出了多极边界元法计算公式和计算过程,分析了截断误差,说明截断误差可由截断项数控制,并给出一个可广泛应用于实际计算的截断项数的近似表达式。     

快速多极边界元方法在二维声散射问题中的应用

快速多极算法(FMM)是求解边界元方法(BEM)在大尺度情况下的一种非常有效的算法.研究了快速多极算法在二维声散射问题的边界积分方程求解中的应用.给出了积分核函数以及其共轭积分算子核函数的多极展开式,局部展开式以及相应展开系数之间的转化关系.分别应用两种不同的层级树结构的FMM来进行求解,并对两种树结构下的求解效率进行了对比.数值算例表明用快速多极算法求解该问题时在存储量和计算量上比直接求解方法效率更高.     

FMBEM求解直圆柱线性波绕射问题

运用快速多极子边界元法(Fast Multipole Boundary Element Method,FMBEM)求得单圆柱在线性波浪中的绕射问题的数值解.所谓的快速多极子边界元法就是采用快速多极子法(Fast Multipole Method,FMM)加速传统边界元法的求解速度.在文中通过求解二维的Helmholtz方程证明FMBEM法具有高精度和高效率,适用于求解大规模的数值问题.另外,给出了单圆柱线性平面波绕射问题中相关水动力学系数的数值计算结果.     

二维各向异性位势薄体问题的虚边界元法

传统边界元法分析各向异性薄体问题时涉及奇异边界积分和拟奇异边界积分的处理,估计这些积分具有相当的难度而且耗时.提出了求解二维各向异性位势薄体问题的虚边界元方法,给出了求解此类问题的新途径,同时拓展了虚边界元法的应用范围.数值算例表明,虚边界元法可有效求解二维各向异性位势薄体问题,且方法简单、精度高、易于程序设计.     

三维位势快速多极虚边界元最小二乘法

将快速多极展开法(FMM)和广义极小残值法(GMRES)结合于三维位势问题的虚边界元最小二乘法,使求解方程的计算量和储存量与所求问题的计算自由度数成线性比例;欲达到数值模拟大规模自由度问题的目的.基于位势问题虚边界元最小二乘法的数值求解格式,将对角化和指数展开系数的概念引入到常规的快速多极展开法中,将三维位势问题的基本解推导为更适合于快速多极算法的展开格式,并用广义极小残值法求解方程组,旨在达到进一步提高效率且仍保证较高计算精度的目的.数值算例说明了该方法的可行性,及计算效率和计算精度.     

二维位势问题快速多极虚边界元解

快速多极虚边界元法是近期发展起来的一种数值算法;其对大规模复杂问题的计算,能在保证求解精度的前提下,使计算量和存储量均比常规虚边界元法具有在数量级上的减少.本文给出了二维位势问题快速多极虚边界元法的求解思想,并进行了数值论证;由文中的数值结果可知,本文方法具有可行性,且有较好的计算精度.     

三维快速多极边界元高性能并行计算

该文实现了快速多极边界元法的一种高性能并行计算。其并行求解器基于自适应新版本快速多极边界元算法,采用三维二次等参元和等精度积分格式,并通过实测的任务量进行分布式并行环境下的合理负载划分。数值算例表明,该求解器在保持高次边界元高精度优点的基础上,对于几何形状不规则的结构仍能保持较好的并行效率,和传统边界元法相比使解题规模有了数量级的提高。这种并行计算为边界元法在大规模复杂工程问题中的应用提供了有效方案。     

快速多极边界元法中的优化数值技术应用

快速多极边界元法是近几年发展起来的边界元新型数值算法,利用多极边界元法解题的关键和难点是求解大规模稀疏矩阵方程组.引入最优化数值技术很好地解决了这一问题,并通过数值实验验证,该方法可节约求解时间,从而为求解大规模问题奠定了理论基础.     

用线性边界元分析二维静电场问题

从二维拉普拉斯方程的积分解和线性插值函数出发，推导出二维静电场边值问题的线性单元的边界元方程及电位、场强计算公式。编制用线性边界元法分析了二维静电场边值问题计算程序，并对几个算例进行了计算。所获得的计算结果与相关文献中所报道的数据吻合较好，表明该方法的有效性；同时还分析计算了十字形内导体圆外导体特种同轴传输线的电位、场强和特性阻抗。     

Multipole BEM for 3-D Elasto-Plastic Contact with Friction

The analysis of 3-D elasto-plastic contact with friction is a highly nonlinear problem. The elements in the contact and plastic zones should be refined to obtain accurate information about the real size,displacement, and traction in the contact zone. However, the increase in the number of degrees of freedom is limited when traditional boundary element method (BEM) is used with the larger memory size and long CPU time required for the solution procedure. This paper describes the additional mathematical friction model to the 3-D elastic multipole BEM to develop a 3-D elasto-plastic contact multipole BEM. Numerical tests show that with this new method, the needed computer memory size is only 2% of the traditional BEM model with friction, which erases large-scale computing with refined meshes and improves the computational accuracy.     

Taylor级数多极边界元法及其在轧制工程中的应用

通过基本解的多极展开与边界元线性方程组的隐式求解方法(GMRES)相结合,开发出了快速多极边界元法。Taylor级数多极边界元法更新了传统边界元法的求解模式,大大提高了计算效率,扩大了边界元法的求解规模。介绍了Taylor级数多极边界元法的发展历史和现状,给出了Taylor级数多极边界元法的基本思想、基本原理和分类,给出了基本解的Taylor展开方法和边界积分的基本实现步骤。将该方法应用于轧制工程中,通过轧辊弹性变形和HC轧机辊系接触和变形的数值解析,说明了Taylor级数多极边界元法适合于大规模轧制工程     

Large Scale Analysis of Mechanical Properties in 3-D Fiber-Reinforced Composites Using a New Fast Multipole Boundary Element Method

Fiber-reinforced composites are commonly used in various engineering applications. The mechanical properties of such composites depend strongly on micro-structural parameters. This paper presents a new boundary element method (BEM) for numerical analysis of the mechanical properties of 3-D fiber-reinforced composites. Acceleration of the BEM is achieved by means of a fast multipole method (FMM), in allowing large scale simulations of a finite elastic domain containing up to 100 elastic fibers to be performed on one personal computer. The maximum number of degrees of freedom can reach a value of over 250 000. The effects of several key micro-structural parameters on the local stress fields and on the effective elastic moduli of fiber-reinforced composites are evaluated. The numerical results are compared with analytical predictions and good agreement is observed. The results show that the fast multipole BEM could be a prom- ising tool for further understanding of the mechanical behavior of such composites.     

Fast Multipole BEM for 3-D Elastostatic Problems with Applications for Thin Structures

The fast multipole method (FMM) has been used to reduce the computing operations and memory requirements in large numerical analysis problems. In this paper, the FMM based on Taylor expansions is combined with the boundary element method (BEM) for three-dimensional elastostatic problems to solve thin plate and shell structures. The fast multipole boundary element method (FM-BEM)requires O(N) operations and memory for problems with N unknowns. The numerical results indicate that for the analysis of thin structures, the FM-BEM is much more efficient than the conventional BEM and the accuracy achieved is sufficient for engineering applications.     

Mathematical Programming Solution for the Frictional Contact Muitipole BEM

This paper presents a new mathematical model for the highly nonlinear problem of frictional contact. A programming model, multipole boundary element method (BEM), was developed for 3-D elastic contact with friction to replace the Monte Carlo method. A numerical example shows that the optimization programming model for the point-to-surface contact with friction and the fast optimization generalized minimal residual algorithm (GMRES(m)) significantly improve the analysis of such problems relative to the conventional BEM.     

Mathematical Programming Solution for the Frictional Contact Multipole BEM

IntroductionElastic friction contact problems require accuratetracking of the movement of objects before and aftercontact and the interaction during contacts and correctsimulation of the frictional behavior between the con-tact surfaces. The boundary element method (BEM)[1,2]is well suited to accurately describe the variation of thefrictional contact conditions since the highly nonlineareffects only occur on the contact surface. For nonlinear frictional contact, various approacheshave been …     

二维位势自适应边界元中的一种误差估计方法

在前期线弹性研究工作的基础上，针对二维位势问题，提出一种新的误差指示来探测位势边界元分析的误差分布。这种误差估计是超奇位势导数边界积分方程的一种残余的度量。数值算例表明本文提出的误差指示成功地跟踪描述了真实误差，进而有效地引导了一个h型网格优化过程。     

快速多极边界元法在薄板结构中的应用

基于Taylor级数多极展开研究了边界元快速多极算法(FM—BEM)，并将它应用于薄板结构。算例分析表明FM—BEM的计算时间和存储空间明显少于常规边界元迭代解法。随着问题规模的增大，这种优势将更加突出。