一种新型的边界元法——边界轮廓法

利用传统边界元积分方程的被积函数的散度等于零的特性，提出一种新型的边界元法——边界轮廓法，使求解问题的维数降低两维。对线弹性平面问题，选择二次位移形函数，求得相应的位移和应力势函数，使二维问题的求解转化为边界点的数值计算，给出了边界点的位移和面力及域内点的应力和位移的计算公式。实例计算表明，该方法具有较高的精度。     

弹性力学平面问题的无奇异边界积分方程

将弹性力学平面问题归化成无奇异边界积分方程,避免了传统的边界元法中的柯西主值（CPV）积分和Hadamard-Finite-Parts(HFP)积分的计算,建立完整的数值求解体系。     

弹性力学中一种新的边界轮廓法

利用基本解的特性,将面力积分方程化成仅含有Ｃａｕｃｈｙ主值积分的形式,基于这种边界积分方程,提出了一种新的边界轮廓法,对于三维问题,该方法只须计算沿边界单元界线的线积分,对二维问题,则只需计算边界单元两点的热函数之差,无须进行数值积分计算,实例计算说明该方法是有效的。     

三维弹塑性有限变形问题边界元法中奇异积分的直接分析去奇解法

作为本文作者研究工作的继续,本文提出了处理三维弹塑性有限变形问题边界元法中二次元区域弱奇及Cauchy 主值奇异积分的二次极坐标变换—分析去奇法.该方法先通过适当的二次极坐标变换降低奇异积分的奇异性,然后利用Causs 散度定理去除Cauchy主值积分的奇异性.通过三维弹塑性及三维有限变形问题数值算例说明该方法具有良好的精度及数值稳定性,并且实施较方便.本文方法可直接推广应用于二阶以上高阶元离散模型奇异积分处理.     

用面力边界积分方程求解断裂力学J积分

证明面力边界积分方程被积函数的散度等于零，应用Stokes公式，对平面线弹性问题，将面力边界积分的求解转化为边界点的位移势函数的点值计算。应用边界积分方程的求解结果，推导出J积分亦可表示为边界点的积分势函数的点值计算，无需进行数值积分，实例计算说明该方法的有效性。     

完备型二次边界轮廓法的研究及求解断裂力学Ja,L,M积分

对线弹性平面问题的边界轮廓法，选用完备的二次位移形函数，使求问题的维数降低两维，给出了求解边界位移和面力以及内点应力的求解方法。证明平面弹怀断鲜明力学Ｊａ积分、Ｍ积分、Ｌ积分方程的被积函数的散度均等于零，将它们分别转化为边界点的位移和面力的线性迭加，无需计算数值积分，算例表明，本文方法具有较高的精度。     

弹性理论几类导数边界积分方程之间的变换关系

导数场边界积分方程通常难以应用，因为存在着超奇异主值积分的计算障碍。弹性理论中有几类不同的位移导数边界积分方程，本文采用算子δij和∈ij(排列张量)作用于这些导数边界积分方程，做一系列变换，原有的超奇异积分被正则化为强奇异积分获解。从而建立了这些位移导数边界积分方程之间的转换关系，它们均可以归结为自然边界积分方程。自然边界积分方程仅存在容易计算的Cauchy主值积分。自然边界积分方程分析可直接获得边界应力和位移导数。     

边界元法计算切口多重应力奇性指数

提出采用边界元法直接计算V形切口的多重应力奇性指数。首先在切口尖端挖出一微小扇形域,在该域边界列常规边界积分方程,后将扇形域内的位移场和应力场表示成关于切口尖端距离ρ的渐近级数展开式,回代入切口边界积分方程,离散后得到关于切口奇性指数的代数特征方程,从而求解获得V形切口的应力奇性指数。该法避免了常规边界元法和有限元法在切口尖端附近布置细密单元的缺陷,并可同时求得多阶应力奇性指数。     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

完备型二次边界轮廓法的研究及J积分和应力强度因子的计算

选择二次完全多项多作为位移形函数,对边界轮廓法作了进一步的发展,证明二维弹性断裂问题的Ｊ积分方程的被积分函数的散度等于零,将Ｊ积分化为边界点的势函数数值的计算,无需计算数值积分,算例表明,该方法较传统边界元法求得的结果精度更好。     

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

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

Stress singularity at a crack tip for various intermediate zones in bimaterial structures (mode III)

The influence of the geometry of a thin intermediate zone on the stress distribution has been investigated in the vicinity of a crack tip in a bimaterial structure. Corresponding modelling boundary value problems are reduced to functional-difference equations by the Mellin transform technique, and later to singular integral equations with fixed point singularities. It has been observed that the order of the stress singularity is essentially dependent on the model parameters. Numerical results concerning the stress singularity exponents and generalized stress intensity factors are presented.     

Evaluation of strongly singular domain integrals for internal stresses in functionally graded materials analyses using RIBEM

An accurate evaluation of strongly singular domain integral appearing in the stress representation formula is a crucial problem in the stress analysis of functionally graded materials using boundary element method.To solve this problem,a singularity separation technique is presented in the paper to split the singular integral into regular and singular parts by subtracting and adding a singular term.The singular domain integral is transformed into a boundary integral using the radial integration method.Analytical expressions of the radial integrals are obtained for two commonly used shear moduli varying with spatial coordinates.The regular domain integral,after expressing the displacements in terms of the radial basis functions,is also transformed to the boundary using the radial integration method.Finally,a boundary element method without internal cells is established for computing the stresses at internal nodes of the functionally graded materials with varying shear modulus.     

Investigation of the scattering of antiplane shear waves by two griffith cracks using the non-local theory

In this paper, the scattering of harmonic antiplane shear waves by two finite cracks is studied using the non-local theory. The Fourier transform is applied and a mixed boundary value problem is formulated. Then a set of triple integral equations is solved using a new method, namely Schmidt's method. This method is more exact and more reasonable than Eringen's for solving this kind of problem. The result of the stress near the crack tip was obtained. Contrary to the classical elasticity solution, it is found that no stress singularity is present at the crack tip, which can explain the problem of macroscopic and microscopic mechanics.     

The scattering of harmonic elastic anti-plane shear waves by a finite crack in infinitely long strip using the non-local theory

In this paper, the scattering of harmonic anti-plane shear waves by a finite crack in infinitely long strip is studied using the non-local theory. The Fourier transform is applied and a mixed boundary value problem is formulated. Then a set of dual integral equations is solved using the Schmidt method instead of the first or the second integral equation method. A one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress occurring at the crack tips. Contraty to the classical elasticity solution, it is found that no stress singularity is present at the crack tip. The non-local dynamic elastic solutions yield a finite hoop stress at the crack tip, thus allowing for a fracture criterion based on the maximum dynamic stress hypothesis. The finite hoop stress at the crack tip depends on the crack length, the width of the strip and the lattice parameter. Supported by the Post Doctoral Science Foundation of Heilongjiang Province, the Natural Science Foundation of Heilongjiang Province and the National Foundation for Excellent Young Investigators.     

关于边界元方法中柯西主值积分的探讨

本文对边界元方法中的各类积分根据其奇异性作分类，并对主值积分的收敛条件、变量替换等进行了讨论，又给出了变替换附加项显式。文中提供的主值积分配项消奇术在边界元方法中是有普遍意义的。     

Investigation of the Scattering of Harmonic Elastic Waves by Two Collinear Symmetric Cracks Using the Non-Local Theory

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｌａｓｔｆｏｕｒｄｅｃａｄｅｓｈａｖｅｗｉｔｎｅｓｓｅｄｔｈｅｉｎａｕｇｕｒａｔｉｏｎｏｆａｎｏｖｅｌｔｈｅｏｒｙｏｆｍａｔｅｒｉａｌｂｏｄｉｅｓ,ｎａｍｅｄｔｈｅｎｏｎ＿ｌｏｃａｌｍｅｃｈａｎｉｃｓ.ＴｈｉｓｗａｓｄｏｎｅｐｒｉｍａｒｉｌｙｄｕｅｔｏｔｈｅｅｆｆｏｒｔｓｏｆＥｄｅｌｅｎ[1],Ｅｒｉｎｇｅｎ[2 ],ＧｒｅｅｎａｎｄＲｉｖｌｉｎ[3].Ａｃｃｏｒｄｉｎｇｔｏｔｈｅｎｏｎ＿ｌｏｃａｌｔｈｅｏｒｙ ,ｔｈｅｓｔｒｅｓｓａｔａｐｏｉｎｔＸｉｎａｂｏｄｙｄｅｐｅｎｄｓｎｏ…     

两相材料V形切口应力强度因子边界元分析

建立了边界元法计算两相材料粘结V形切口奇异应力场的新途径。在V形切口尖端挖出一小扇形,将该扇形弧线边界的位移和面力表示为有限项奇性指数和特征角函数的线性组合,其组合系数即为广义应力强度因子,将该组合回代到在被挖去小扇形后的剩余结构内建立的边界积分方程,离散后可求解出组合系数,获得两相材料粘结V形切口尖端的应力强度因子。算例证明了本文方法的有效性。     

Analysis of a three-dimensional dissimilar material joint with one real singularity using a conservative integral

In the present study, a conservative integral based on the Betti reciprocal principle is formulated to determine the intensity of singularity at a vertex of the interface in three-dimensional dissimilar material joints with one real singularity. Eigenanalysis formulated using a three-dimensional finite element method (FEM) is used to calculate the order of stress singularity, angular functions of displacements and stresses. Models with various element sizes and various integral areas are used to investigate the effect of the integration area on the accuracy of the results. The results are compared with those obtained from the boundary element method (BEM) using a curve-fitting technique to calculate the intensity of singularity. In addition, models of various lengths and various material combinations are used to investigate the stress singularity characteristics in three-dimensional dissimilar material joints. The results of the present study indicate that the conservative integral can be used to determine the intensity of singularity in three-dimensional bi-material joints. The accuracy of the results can be improved by mesh refinement. Finally, the relationships among the intensity of singularity, the order of stress singularity and the model geometry are discussed.