二维弹性问题边界元法中边界层效应问题的变换法

基于间接规则化边界积分方程,有效估计奇异边界积分,准确求得边界量,为场变量的计算奠定了基础。在计算场变量时,针对二维弹性力学边界元法中出现的几乎奇异积分,本文采用一类非线性变量替换法,有效地改善了被积函数的震荡特性,从而消除了核积分的几乎奇异性;在不增加计算量的情况下,极大地改进了几乎奇异积分计算的精度,成功地求解了弹性体近边界点上的力学参量,避免了边界层效应。此外,本文引入一种精确几何单元逼近,对于圆弧边界,这样的插值逼近几乎是精确的,提高了计算精度。数值算例表明,本文算法稳定,效率高,并可达到很高的计算精度,即使场点非常靠近边界,如场点到积分单元的距离小到纳米级,仍可避免边界层效应现象。     

三维变系数热传导问题边界元分析中几乎奇异积分计算

在边界积分的数值计算过程中,当源点离积分单元很近时,边界积分就会具有几乎奇异性,此时不能直接用高斯数值积分公式计算几乎奇异积分。本文以三维非均质热传导问题为例,介绍了一种计算几乎奇异边界积分的新方法。首先,采用Newton-Raphson迭代算法确定积分单元上离源点最近的点;然后,将积分单元上任意一点的坐标在最近点处展开成泰勒级数,并计算源点到积分单元任意点的距离;最后,将距离函数代入几乎奇异边界积分中,并运用指数变换方法导出积分单元上几乎奇异积分的计算公式。文中给出了两个非均质热传导问题的算例来验证所述方法的正确性、有效性和稳定性。     

基于坐标变换精确计算近奇异积分的双向sinh变换法

提出一种与坐标变换相结合的双向sinh变换方法用于精确计算近奇异积分.首先利用坐标变换法对近奇异积分进行双向分离,再针对分离出的两个方向的积分形式,根据复变函数极点理论,构造双向sinh变换.与传统结合环向变换和极坐标变换的单向sinh变换以及迭代sinh变换方法相比,双向sinh变换法的计算精度更高,并且最近点的位置对近奇异积分计算精度的影响明显变小.     

高阶边界元奇异积分的一种通用高效计算方法

针对三维边界元法中曲面单元上的(弱、强、超)奇异积分提出了一种通用高效的计算方法。经极坐标变换,将奇异积分转化为常规积分;采用数值方法计算Cauchy主值积分和Hadamard有限项积分系数;引入保角变换和反曲变换消除因单元畸形或因积分点靠近单元边界而引起的周向积分奇异性。该方法可以统一处理(弱、强、超)奇异积分,并且只需要知道核函数的奇异阶数和少数几个点上的被积函数值,不依赖于积分和函数的具体选取;所需的积分点少,精度高,并且受单元畸形程度影响较小,稳定性好。采用该方法计算了声学和弹性力学中的典型奇异积分,并结合二阶Nystrm方法求解了弹性力学的边界积分方程,验证了方法的高精度和高效性。本文数值积分程序可向作者索取。     

基于球面细分法的高效高精度近奇异积分计算

精确高效地计算近奇异积分,对边界元法的成功实施至关重要,也是边界元法在实际工程计算中面临的主要障碍之一。论文提出了一种基于球面细分技术的近奇异积分计算方法,可以精确计算任意基本解类型、任意单元形状和任意源点位置的近奇异积分。该方法首先通过计算源点到单元的最近最远距离,来确定球面细分的初始半径和终止半径;然后通过一系列半径呈指数级增长的球面来分割积分单元,得到一系列三角形和四边形子单元;最后把细分后得到的子单元变成弧形状,即三角形和四边形子单元分别变成扇形和环形子单元。由于球面细分是直接在三维笛卡尔坐标系下进行的,所以它适用于任何类型的单元。此外,由于基本解主要是源点到场点距离的函数,因此在同等精度下,近奇异积分在子单元的环向上所需要的高斯积分点数将大大减少。在径向方向上,由于球半径系列呈指数级变化,各个子块可以做到等精度高斯积分。数值算例表明,与传统近奇异积分计算方法相比,论文提出的方法更加稳定,精度更高。     

A GENERAL ALGORITHM FOR CALCULATING NEARLY SINGULAR INTEGRALS OVER HIGH-ORDER ELEMENTS IN 3D BEM

三维边界元分析中,高阶几何单元上的几乎奇异积分计算是一个重要而且困难的问题,该文对此进行了研究。使用8节点四边形和6节点三角形曲面单元来描述几何边界;构造了新的距离函数;拓展原有的指数函数非线性变换到三维边界元法中,利用拓展的变换来消除被积函数的几乎奇异性。数值算例表明,该算法稳定,效率高,即使计算点到实际边界的距离很小,依然可获得令人满意的数值解。     

三维边界元法高阶元几乎奇异积分半解析法

分析了三维边界元法高阶曲面单元几何特征,定义接近度来表征源点与积分单元的接近程度.利用源点在积分单元上的垂足点建立局部极坐标系,构造与几乎奇异积分核函数具有相同奇异性的近似函数.从奇异积分核函数中扣除其近似函数,分离出积分核中主导的奇异函数部分,将奇异积分分解为规则核函数和奇异核函数两项积分.规则核函数积分应用常规Gauss数值积分计算,奇异核函数积分在局部极坐标系ρθ下分离积分变量ρ和θ,对ρ积分建立解析计算列式,对θ积分应用常规Gauss数值积分计算,从而对三维位势问题高阶边界单元几乎强奇异和几乎超奇异积分建立一种新的半解析算法.给出了若干温度场算例,采用边界元法高阶单元几乎奇异积分半解析法计算了近边界内点位势和位势梯度,并与线性单元正则化算法计算结果对比,结果证明提出的半解析法计算几乎奇异面积分和薄壁结构更加高效.      

薄体位势问题边界元法中的解析积分算法

薄体结构的数值分析是边界元法的难点问题之一。该文导出了一种完全解析积分算法，用这种算法计算了薄体平面位势问题边界元法中出现的几乎弱奇异、强奇异和超奇异积分。当边界离散为一系列线性单元，边界积分方程离散计算的积分可归纳为三种形式。对薄体问题，源点与积分单元距离通常相距很近，这些积分产生显著几乎奇异性，直接采用常规高斯积分不能有效计算。为此该文导出了这些几乎奇异积分的全解析计算公式。按源点与单元的距离是否为零，公式分两种情况。新算法采用全解析积分公式处理几乎奇异积分，首先精确计算出薄体问题边界未知位势和法向位势梯度，然后再进一步计算了域内点的物理参量。算例表明该文算法可处理狭长比为1．E-08的薄体问题，显示了边界元法分析薄体问题具有独特的优势。     

边界积分方程中近奇异积分计算的一种变量替换法

准确估计近奇异边界积分是边界元分析中一项很重要的课题,其重要性仅次于对奇异积分的处理. 近年来已发展了许多方法,都取得了一定程度的成功,但这个问题至今仍未得到彻 底的解决. 基于一种新的变量变换的思想和观点,提交了一种通用的积分变换法, 它非常有效地改善了被积函数的震荡特性,从而消除了积分的近奇异性,在不增加计算量的情况 下, 极大地改进了近奇异积分计算的精度. 数值算例表明,其算法稳定,效率高, 并可达到很高的计算精度,即使区域内点非常地靠近边界,仍可取得很理想的结果.     

三维Laplace方程边界元中线性单元的精确积分法

边界元中的边界积分计算影响计算精度和计算速度。非奇异积分一般采用数值积分，当配置点接近积分单元时，计算精度降低。未知函数线性插值得到的解是连续解，但计算难度增大。本文采用积分区域变换，将三维Laplace问题的二维积分化为一维积分，这样奇异积分和非奇异积分能采用精确积分的方法计算，使求解精度，计算速度都得到提高。     

结合材料界面端的三维应力奇异性

本文利用特殊有限元方法，开发了一个用来求解结合材料界面端三维应力奇异性问题的数值分析程序。该方法只需对界面端的角度方向进行离散即可求得应力奇异性。结合材料的应力奇异性取决于两种材料的材料常数和界面端形状。选用三个材料参数作为变量，用来研究结合材料三维应力奇异性随材料常数的变化规律。文中计算了几种重要而且常见的情况，并以此为基础建立了数据库。同时，还分析了应力奇异性随界面端形状的变化规律，并得到了应力函数的分布图。     

Thermoelectroelastic Green's function and its application for bimaterial of piezoelectric materials

Summary For a two-dimensional piezoelectric plate, the thermoelectroelastic Green's functions for bimaterials subjected to a temperature discontinuity are presented by way of Stroh formalism. The study shows that the thermoelectroelastic Green's functions for bimaterials are composed of a particular solution and a corrective solution. All the solutions have their singularities, located at the point applied by the dislocation, as well as some image singularities, located at both the lower and the upper half-plane. Using the proposed thermoelectroelastic Green's functions, the problem of a crack of arbitrary orientation near a bimaterial interface between dissimilar thermopiezoelectric material is analysed, and a system of singular integral equations for the unknown temperature discontinuity, defined on the crack faces, is obtained. The stress and electric displacement (SED) intensity factors and strain energy density factor can be, then, evaluated by a numerical solution at the singular integral equations. As a consequence, the direction of crack growth can be estimated by way of strain energy density theory. Numerical results for the fracture angle are obtained to illustrate the application of the proposed formulation. Received 10 November 1997; accepted for publication 3 February 1998     

Plane steady vibration of an elastic wedge

The problem of plane steady vibration of an elastic wedge of arbitrary angle (less than 180 degrees) subject to harmonic normal and shearing tractions on its faces is reduced to a system of singular integral equations by the superposition of two half-plane solutions. The integral equations have kernels with Cauchy singularities of a non-translation type, except for the 90 degree wedge. The locations of these singularity lines are shown graphically as a function of wedge angle.     

Isotropic/orthotropic two-layered strips containing cracks terminating at and going through an interface

Crack problems for isotropic/orthotropic two-layered strips have been investigated. A system of two singular integral equations can be derived by using Fourier integral transformation and boundary conditions of crack problems. After stress singularities at crack tips or other special points are determined for internal and edge cracks, and for cracks terminating at and going through the interface, the system of singular integral equations is solved numerically by Gauss-Jacobi or Gauss-Chebyshev integration formulas for stress intensity factors at the tips and other singular points of cracks. Finally, possible crack growth behavior for cracks approaching and going through the interface is discussed.     

THE INTERACTION BETWEEN MULTIPLE CURVED RIGID LINE AND ANTIPLANE CIRCULAR INCLUSION

The interaction between multiple curved rigid line and ciruclar inclusion in antiplane loading condition is considered in this paper. By utilizing the point force elementary solutions and taking density function of traction difference along curved rigid lines, a group of weakly singular integral equations with logarithmic kernels can be obtained. After the numerical solution of the integral equations, the discrete values of density functions of traction difference are obtainable. So the stress singularity coefficient at rigid line tips can be calculated, and two numerical examples are given.     

Meshless acoustic analysis using a weakly singular Burton-Miller boundary integral formulation

The Burton-Miller boundary integral formulation is solved by a complex variable boundary element-free method (CVBEFM) for the boundary-only meshless analysis of acoustic problems with arbitrary wavenumbers. To regularize both strongly singular and hypersingular integrals and to avoid the computation of the solid angle and its normal derivative, a weakly singular Burton-Miller formulation is derived by considering the normal derivative of the solid angle and adopting the singularity subtraction procedures. To facilitate the implementation of the CVBEFM and the approximation of gradients of the boundary variables, a stabilized complex variable moving least-square approximation is selected in the meshless discretization procedure. The results show the accuracy and efficiency of the present CVBEFM and reveal that the method can produce satisfactory results for all wavenumbers, even for extremely large wavenumbers such as k = 10 000.     

Method for determining power-type complex singularities in solutions of singular integral equations with generalized kernels and complex conjugate unknowns

We develop a method for determining power-type complex singularities of solutions for a class of one-dimensional singular integral equations with generalized kernels and complex conjugate unknown functions. By analyzing the characteristic part of a singular integral equation, we reduce the problem of determining the solution singularity exponents at the ends of the integration interval to two independent transcendental equations for these exponents. We show that the distribution of admissible singularity exponents is of continuous character. We present numerical results for a two-dimensional elasticity problem whose mathematical statement leads to a singular integral equation of the class under study. We also reveal the drawbacks of one classical approach to the determination of stress field singularities.     

Anti-plane shear crack normal to and terminating at the interface of two bonded piezoelectric ceramics

The electroelastic analysis of two bonded dissimilar piezoelectric ceramics with a crack perpendicular to and terminating at the interface is made. By using Fourier integral transform, the associated boundary value problem is reduced to a singular integral equation with generalized Cauchy kernel, the solution of which is given in closed form. Results are presented for a permeable crack under anti-plane shear loading and in-plane electric loading. Obtained results indicate that the electroelastic field near the crack tip in the homogeneous piezoelectric ceramic is dominated by a traditional inverse square-root singularity, while the electroelastic field near the crack tip at the interface exhibits the singularity of power law r^−α, r being distance from the interface crack tip and α depending on the material constants of a bi-piezoceramic. In particular, electric field has no singularity at the crack tip in a homogeneous solid, whereas it is singular around the interface crack tip. Numerical results are given graphically to show the effects of the material properties on the singularity order and field intensity factors.     

一种处理弹性动力边界元法中奇异积分的方法

利用边界元法求解瞬态弹性动力学问题时,时域基本解函数的分段连续性和奇异性为该问题的求解带来很大的困难。为了解决时域基本解中的奇异性问题,本文依据柯西主值的定义,对经过时间解析积分之后的时域基本解进行奇异值分解,将其分成奇异和正则积分两部分;其中正则部分可通过采用常规高斯积分方法来计算,而奇异部分具有简单的形式,可以利用解析积分计算。经过上述操作之后,就可以达到直接消除时域基本解中奇异积分的目的。和传统方法相比,本文方法并不依赖静力学基本解来消除奇异性,是一种直接求解方法。最后给定两个数值算例来验证本文提出方法的正确性和可行性,结果表明使用本文算法可以解决弹性动力学边界积分方程中的奇异性问题。     

Analytical removal of singularity and direct evaluation of singular integrals in 3-D elastoplastic finite deformation analysis by BEM

As a further development of the present authors' research work [1,2], in this paper a method of the so-called quadratic pentahedron polar co-ordinate transformation and analytical removal of singularity of Cauchy principal value singular integrals is proposed to evaluate the strongly singular integrals in the sense of Cauchy principal values and the weakly singular integrals over quadratic internal cells in 3-D elastoplastic finite deformation analysis by BEM. First, a quadratic pentahedron polar co-ordinate transformation technique is used to reduce the order of singularity of the singular integrals. Then, a form of Gauss' theorem is introduced to remove the singularity in the Cauchy principal value singular integrals analytically. Therefore, the evaluation of all those strongly and weakly singular integrals can be carried out by standard Gaussian quadrature accurately and efficiently. Numerical examples of the 3-D elastoplastic problem and 3-D finite deformation problem are given to demonstrate that the method possesses good accuracy and numerical stability, and is convenient to implement. The method in this paper can be applied extensively to evaluating the singular integrals over cubic and higher order elements.