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

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

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.     

Fast Multipole BEM for Simulation of 2-D Solids Containing Large Numbers of Cracks

The fast multipole method was used to solve the traction boundary integral equation for 2-D crack analysis, The use of both multipole and local expansions reduces both the computational complexity and the memory requirement to O(N). The multipole expansion uses a complex Taylor series expansion to reduce the number of multipole moments, The generalized minimum residual method solver (GMRES) was selected as the iterative solver, An improved preconditioner for GMRES was developed which uses less CPU time and less memory. A new initial candidate vector for the iterative solver was developed to further improve the efficiency, The numerical examples apply the method to the analysis of cracks in infinite 2-D space with the largest model having 900 000 degrees of freedom.     

一种新的用于二维弹性静力学的快速多极边界元法

快速多极边界元法(fastmultipole BEM)是近几年发展起来的边界元新型算法。本文提出了一种新型的适合二维弹性静力学问题的快速多极边界元格式，并用于含有多个夹杂的二维复合材料的应力分析。数值结果表明这种方法非常适合解决大规模问题。     

三维弹塑性摩擦接触多极边界元法

三维弹塑性摩擦接触问题是多重非线性问题，对接触区表面和塑变区的离散，须划分大量单元进行大规模运算才能获得接触位移、面力及应力场的准确信息。传统边界元法由于离散自由度所需内存大，CPU计算时间冗长，完整解题运算规模受到限制。本文在三维弹性多极留数边界元法的基础上，开发研制三维弹性数学规划型摩擦接触多极边界元法及源程序，建立三维弹塑性摩擦接触多极边界元法并研制其源程序，更新了课题组开发的原传统的三维弹塑性摩擦接触边界元法。数值试验表明，本法使计算机内存量减少近百倍，从而使细划分单元的大规模运算成为可能，并提高了计算精度。     

三维快速多极虚边界元配点法

以三维弹性力学问题为研究背景,提出了一种三维快速多极虚边界元配点法的求解思想,即将三维快速多极展开的基本思想和广义极小残值法运用于求解传统虚边界元配点法方程.文中将三维弹性问题的基本解推导为适合于虚边界元快速多极算法的展开格式,经数值计算格式的演变,使求解方程的计算量和储存量与所求问题的计算自由度数成线性比例,以达到数值模拟大规模自由度问题的目的.算例说明了该方法的可行性、计算效率和计算精度.此外,该方法的思想具有一般性,应用上具有扩展性.     

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 …     

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.     

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

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

Multi-Variable Non-Singular BEM for 2-D Potential Problems

A multi-variable non-singular boundary element method (MNBEM) is presented for 2-D potential problems. This method is based on the coincident collocation of non-singular boundary integral equations(BIEs) of the potential and its derivatives, where the nodal potential derivatives are considered independent of the nodal potential and flux. The system equation is solved to determine the unknown boundary potentials and fluxes, with high accuracy boundary nodal potential derivatives obtained from the solution at the same time. A modified Gaussian elimination algorithm was developed to improve the solution efficiency of the final system equation. Numerical examples verify the validity of the proposed algorithm.     

压力容器开孔结构分析的快速多极边界元研究

压力容器开孔结构的应力高度集中。为精确模拟该结构的应力状况,该文提出一种三维高阶快速多极边界元法。在三维弹性力学边界元法的基础上,推导出二阶单元的基本解快速多极展开格式。该算法通过多极展开概念,大大降低了对存储量的要求,并且不损失精度。使用高阶快速多极边界元法分析含多个开孔的压力容器整体结构,所得应力结果与大规模高阶有限元法的结果吻合得很好。研究结果表明,高阶快速多极边界元法易于分析此类大规模问题,并具有很高的数值计算精度,满足工程设计的要求。     

磁多极场的场参数理论

通过对磁多极场的分析,其对称性可分为两大类:一类是二、六、十、十四、…极场,另一类是四、八、十二、十六、…极场.根据磁多极场的对称性,导出了磁多极场磁感应三个分量的泰勒级数展开式,定义了磁多极场的场参数.结果表明,通过y轴上By分量的数个参数,就可以近似计算得到整个空间的磁感应值.     

求解三维准静态矢量磁位的快速多极方法

将快速多极算法（FMM）应用于三维准静态电磁场矢量磁位的求解，首先根据计算精度的要求把连续分布的场源进行离散化处理，然后通过静电类比分析，将求解三维准静态矢量磁位的问题转化为多体问题，进而利用快速多极方法来计算三维空间中载流导体产生的矢量磁位，可以将计算量由O（N^2）降低为O（N）次运算，大大提高了计算速度．算例的计算结果表明，当取剖分体积单元的边长等于0．25倍透入深度时，采用FMM方法计算的电流密度不均匀分布载流导体在其自身所在空间的磁矢位与精确解的相对误差小于0．005，而其在自身所在空间以外的磁矢位的FMM计算结果，具有更高的精度．经过积分方程离散和静电模拟分析，应用FMM算法可正确地计算三维空间载流导体的矢量磁位，计算误差可通过剖分密度进行控制．提出的方法扩展了FMM算法在准静态矢量磁位数值计算领域中的应用，为芯片上互连电感参数的计算奠定了基础．     

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

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

三维电大尺寸复杂群目标的单站RCS的快速多极子分析

用快速多极子算法(FMM)和共轭梯度法(CG)求解三维电大尺寸复杂群目标的电磁散射特性。对单站雷达散射截面(RCS)的预估,更采用了物理光学电流近似和相位修正的继承迭代法两项措施进一步加快了求解过程。该方法具有节省内存,计算量小,迭代速度快且精确度高的特点,特别适于准确分析多个电大尺寸目标间的相互影响。用三维计算实例验证了该方法在解决电大尺寸复杂群目标电磁散射分析方面的有效性和优越性。     

静态电场分析的多极理论—边界元耦合法

从Ｌａｐｌａｃｅ方程的积分解出发,推导出静态电场分析的多极理论－边界元耦合的公式,给出了多极理论－边界元耦合法的计算公式及其使用条件．通过对两个实例的计算表明,所建立的方法是可行的,且具有比其它数值方法更高的计算精度,可以很方便地应用于静态电场的分析与设计中．这种方法也可用于求解Ｐｏｉｓｏｎ方程的边值问题．     

Taylor级数多极边界元法远场影响的误差估计

快速多极方法能够有效地提高边界元法的计算效率．求解的计算量和内存量与问题的自由度数N成正比．求解的精度与传统边界元法相比有所下降．分析了Taylor级数多极边界元法的计算精度和远场影响系数的误差．研究了核函数r的Taylor级数展开性质,推导了三维弹性问题基本解的误差估计公式．说明了影响多极边界元法计算精度的因素．数值算例显示了误差估计公式的正确性和有效性．     

层状地基与弹性薄板相互作用的边界元解

将无限大薄板的基本解作为薄板边界积分方程的核函数,对薄板的内部和边界进行离散,并假定薄板内部和边界上的节点与地基反力的分布情况,得到薄板的边界元方程组;同时基于层状地基的解析层元解,通过Guass-Legendre积分得到地基柔度矩阵;结合地基与薄板接触面上的位移协调条件,得到层状地基与薄板共同作用问题总的边界元法方程组;求解该方程组,得到层状地基与薄板共同作用问题的解答.基于本文理论,编制了相应的FORTRAN程序,通过与已有文献结果对比验证本文理论及程序的正确性,数值分析结果表明:方形基础薄板情况下,离板中心越近,垂直于坐标轴y(x)方向、距离相等的2条线段的竖向位移差越小,且该位移差随着板-土刚度比减小而减小;随着板长宽比的增大,板中心点与长边中点位移差变化不明显,而短边中心与边界角点的位移差也有相类似的规律.