热弹性平面问题的规则化边界积分方程

对于热弹性平面问题,过去广泛集中在直接变量边界元法研究,本文研究间接变量规则化边界元法,建立了间接变量规则化边界积分方程。和直接边界元法相比,间接法具有降低密度函数的连续性要求、位移梯度方程中的热载荷体积分具有较弱奇异性等优点。数值实施中,用精确单元描述边界几何,不连续插值函数逼近边界量。算例表明,本文方法效率高,所得数值结果与精确解相当吻合。     

各向异性位势平面问题的规则化边界元法

基于转化域方程为边界积分方程的极限定理及一个新颖的基本解分解技术, 建立间接变量规则化边界积分方程, 它有效地避免了奇异积分的直接计算. 与已有方法比,该方法不将问题变换为各向同性的问题去处理, 因而无需反演运算, 也有别于Galerkin方法, 无需计算重积分. 可计算任意边界位势梯度, 而不仅限于法向通量. 针对椭圆边界的边值问题, 提交一种精确单元来描述边界几何. 数值算例表明, 所提算法稳定且效率高, 所得数值结果与精确解吻合较好.      

无奇异性边界积分方程法

边界元法中如何有效地处理奇异积分,一直是人们极为关心的课题。本文提出了建立由“线源”产生的边界积分方程,其优点是可任意阶地降低奇异性阶次。计算实例表明,本文所提出的方法优于常规的边界元法。尤其是提高了位移的计算精度。对边界单元网格不等长划分,更显示出该法的优越性。     

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

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

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

从理论上提出一种新的方法，归化出间接变量无奇异边界积分方 程. 采用Lagrange二次单元，建立一个数值求解框架系统. 此外，基于问题的计算区域的 特殊性，给出一种边界近似方法. 数值算例表明该方法所取得的数值结果与精确解相当接 近，特别是边界量的数值结果. 此外，该方法容易被推广到三维问题. 和已有的直接变量的情形相比较，具有优点：1)无需处理HFP积分. 大大降低处理问 题的复杂性，并提高了计算效率和解的精度；2)摆脱了问题的具体形式，进入纯代数操作. 这样做的好处是从理论上建立一种普遍适用的方法，不仅适用于弹性力学问题，同样可应用 于其它问题，如位势问题, Stokes问题等. 3)提供了一种计算CPV积分的方法.     

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

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

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

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

导数场边界积分方程分析近边界应力分布

研究二维弹性力学问题边界积分方程,通过分部积分变换消除了常规导数边界积分方程中的超奇异积分,获得仅含强奇异积分的应力自然边界积分方程.对于近边界应力的计算,进一步运用正则化算法解析计算其中的几乎强奇异积分.较常规边界元法相比,应力自然边界积分方程可以求解离边界更加接近的内点应力值.算例证明了文中方法的可应用性和有效性.     

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

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

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

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

ELASTIC ANALYSIS OF ORTHOTROPIC PLANE PROBLEMS BY THE SPLINE FICTITIOUS BOUNDARY ELEMENT METHOD

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｅｃｏｎｖｅｎｔｉｏｎａｌｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ (ＢＥＭ) ,ｓｉｎｇｕｌａｒｐｏｉｎｔｓａｒｅｄｉｓｔｒｉｂｕｔｅｄａｌｏｎｇｔｈｅｂｏｕｎｄａｒｙｏｆｔｈｅｄｏｍａｉｎｕｎｄｅｒｓｔｕｄｙ ,ｗｈｉｃｈｌｅａｄｓｔｏｓｉｎｇｕｌａｒｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎｓｄｕｅｔｏｔｈｅｓｉｎｇｕｌａｒｉｔｙｏｆｆｕｎｄａｍｅｎｔａｌｓｏｌｕｔｉｏｎｓ .Ｔｈｅｒｅｆｏｒｅ ,ｓｉｎｇｕｌａｒｉｎｔｅｇｒａｔｉｏｎｍｕｓｔｂｅｈａｎｄｌｅｄｗｈｅｎｓｏｌ…     

Seismic oscillations of an elastic rectangular structure protected by a viscoelastic stratum: in-plane problem

We study the problem of seismic protection of a rectangular elastic construction by mean of a viscoelastic stratum made of material of Kelvin–Voigt type. In particular we consider the in-plane problem. By using the standard boundary integral equations method (BIEs), the dynamic problem is reduced to a system of BIEs over the set of boundary lines. In order to reduce the dimension of the problem in its discrete form, for wave processes in the stratum and in the half-space foundation we use a special Green’s function satisfying the stress-free boundary conditions over the boundary of the half-space. A numerical algorithm is used to solve the problem after discretization and finally we discuss the physical meaning of the results and the efficiency of the seismic protection by the viscoelastic stratum.     

Boundary element analysis for elastic and elastoplastic problems of 2D orthotropic media with stress concentration

Both the orthotropy and the stress concentration are common issues in modern structural engineering. This paper introduces the boundary element method (BEM) into the elastic and elastoplastic analyses for 2D orthotropic media with stress concentration. The discretized boundary element formulations are established, and the stress formulae as well as the fundamental solutions are derived in matrix notations. The numerical procedures are proposed to analyze both elastic and elastoplastic problems of 2D orthotropic media with stress concentration. To obtain more precise stress values with fewer elements, the quadratic isoparametric element formulation is adopted in the boundary discretization and numerical procedures. Numerical examples show that there are significant stress concentrations and different elastoplastic behaviors in some orthotropic media, and some of the computational results are compared with other solutions. Good agreements are also observed, which demonstrates the efficiency and reliability of the present BEM in the stress concentration analysis for orthotropic media. The project supported by the Basic Research Foundation of Tsinghua University, the National Foundation for Excellent Doctoral Thesis (200025) and the National Natural Science Foundation of China (19902007). The English text was polished by Keren Wang.     

正交各向异性厚板的边界元解法

本文利用 Hormander 算子法和平面波分解法导出了计入剪切变形的正交各向异性厚板的基本解。建立了计入剪切变形的正交各向异性厚板的边界积分方程。文中详细地讨论了基本解的数值计算,并用边界元法分析了一些算例。     

THE MESHLESS VIRTUAL BOUNDARY METHOD AND ITS APPLICATIONS TO 2D ELASTICITY PROBLEMS

A novel numerical method for eliminating the singular integral and boundary effect is processed. In the proposed method, the virtual boundaries corresponding to the numbers of the true boundary arguments are chosen to be as simple as possible. An indirect radial basis function network (IRBFN) constructed by functions resulting from the indeterminate integral is used to construct the approaching virtual source functions distributed along the virtual boundaries. By using the linear superposition method, the governing equations presented in the boundaries integral equations (BIE) can be established while the fundamental solutions to the problems are introduced. The singular value decomposition (SVD) method is used to solve the governing equations since an optimal solution in the least squares sense to the system equations is available. In addition, no elements are required, and the boundary conditions can be imposed easily because of the Kronecker delta function properties of the approaching functions. Three classical 2D elasticity problems have been examined to verify the performance of the method proposed. The results show that this method has faster convergence and higher accuracy than the conventional boundary type numerical methods.     

NOVEL REGULARIZED BOUNDARY INTEGRAL EQUATIONS FOR POTENTIAL PLANE PROBLEMS

The universal practices have been centralizing on the research of regulariza-tion to the direct boundary integal equations (DBIEs). The character is elimination of singularities by using the simple solutions. However, up to now the research of regular ization to the first kind integral equations for plane potential problems has never been found in previous literatures. The presentation is mainly devoted to the research on the regularization of the singular boundary integral equations with indirect unknowns. A novel view and idea is presented herein, in which the regularized boundary integral equations with indirect unknowns without including the Cauchy principal value (CPV) and Hadamard-finite-part (HFP) integrals are established for the plane potential problems. With some numerical results, it is shown that the better accuracy and higher efficiency, especially on the boundary, can be achieved by the present system.     

横观各向同性压电材料中裂纹问题的边界元法

采用Somigiliana公式给出了三维横观各向同性压电材料中的非渗漏裂纹问题的一般解和超奇异积分方程,其中未知函数为裂纹面上的位移间断和电势间断．在此基础上,使用有限部积分和边界元结合的方法,建立了超奇异积分方程的数值求解方法,并给出了一些典型数值算例的应力强度因子和电位移强度因子的数值结果,结果令人满意．     

A multi-variable non-singular BEM in 2D elasticity

A new BEM algorithm is proposed in the paper. This algorithm is based on the coincident collocation of the non-singular boundary integral equations (BIEs) of displacement and its derivatives, and the basic variables of the new BEM system equation are the boundary nodal displacements, tractions and displacement derivatives. Due to the non-singular character of the BIEs, all integrals in the algorithm can be evaluated by the Gaussian quadrature. Furthermore, the conforming element can be used for this algorithm as the displacement derivative can be considered independent of the displacement. Once the new BEM system equation is solved, the boundary stress can be evaluated immediately from the obtained displacement derivatives. Numerical examples show the validity of the algorithm.     

The vibration of a half-space due to a buried mode I crack opening

The solution of a dynamic problem for calculation of a displacement field on a half-space surface caused by an internal mode I crack opening is presented. The problem is reduced to the system of boundary integral equations (BIEs). The equations of motion are solved with the use of Helmholtz potentials and applying Fourier integral transform. The effects of the crack size, the crack depth and the distance from the crack epicenter to the observation point on the parameters of elastic waves are investigated. It is established that the increasing of the defect size leads to narrowing bandwidth of elastic waves and to lowering of center frequency. The analysis given here can be used for identification of the crack growth during technical diagnostic of an industry objects and structural elements by AE method.