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

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

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

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

基于边界元法求解三维弹性摩擦接触问题

边界元法求解三维摩擦接触问题,其中一个关键点在于如何确定滑移方向。即当出现相对滑移时,滑移方向如何确定。当前常采用的方法是,粘结点利用切向面力得到滑移方向,滑移点利用切向相对位移得到滑移方向。不过该方法难以保证收敛性。针对这一问题,本文采用滑移方向预测技术得到滑移方向。即以后出现相对滑移时,滑移方向采用预测技术中得到的滑移方向。由于摩擦接触问题和历史加载相关,本文采用增量法求解。不同摩擦系数下的数值结果都证明了本文算法的有效性和收敛性及滑移方向预测技术的有效性。     

DYNAMIC INTERACTION OF PLANE WAVES WITH A UNILATERALLY FRICTIONALLY CONSTRAINED INCLUSION-TIME DOMAIN BOUNDARY ELEMENT ANALYSIS

A 2D time domain boundary element method (BEM) is developed to solve the transient scattering of plane waves by a unilaterally frictionally constrained inclusion. Coulomb friction is assumed along the contact interface. The incident wave is assumed strong enough so that localized slip and separation take place along the interface. The present problem is in effect a nonlinear boundary value problem since the mixed boundary conditions involve unknown intervals (slip, separation and stick regions). In order to determine the unknown intervals, an iterative technique is developed. As an example, we consider the scattering of a circular cylinder embeddedin an infinite solid.     

受有集中力的样条边界元法

本文应用Ｂｅｔｔｉ定理，建立弹性体边界和体内受有有限个集中力时的样条边界积分方程，将集中力的影响表征为边界积分方程中的自由项，从客观实际出发，对具有两套奇性交会的积分方程给出一种方便有效的处理方法，使得集中力下的边界元法得以实施，在样插值基础上，即使稀疏剖分也能给出很高精度的位移场。应力场和未知集中反力。     

Numerical comparison of two boundary element methods for plane harmonic functions

A direct boundary element method (BEM) has been studied in the paper based on a set of sufficient and necessary boundary integral equations (BIE) for the plane harmonic functions. The new sufficient and necessary BEM leads to accurate results while the conventional insufficient BEM will lead to inaccurate results when the conventional BIE has multiple solutions. Theoretical and numerical analyses show that it is beneficial to use the sufficient and necessary BEM, to avoid hidden dangers due to non-unique solution of the conventional BIE.     

各向异性体对瞬态SH波散射问题的边界元方法

本文利用Fourier变换和加权残数法建立了正交各向异性体反平面瞬态波散射问题的边界积分方程,构造了在边界上能满足波动方程的边界元函数,它比常用的二次元和高次元具有较好的适应性。数值结果表明利用该形函数计算高频情形下的散射问题具有很高的精度。     

单侧摩擦接触约束夹杂物与平面波的动力相互作用--时域边界元分析

应用时域边界元法研究了瞬态平面波对单侧摩擦约束夹杂物的散射问题，假设界面摩擦遵守库仑定理，当入射波足够强时界面会出现局部分离和滑移，由于边界上的区域(分离区、滑移区和粘着区)是未知的，所以该问题实际上是个复杂的边界非线性问题，为了确定未知区域，该文发展了一种有效的迭代技术，作为算例，计算了一个无限域中圆柱埋置夹杂物对瞬态平面波的散射问题。     

A boundary element method for calculating the SH waves scattering from arbitrarily shaped holes in an anisotropic medium

The scattering problem of elastic wave by arbitrarily shaped cavities in an infinite anisotropic medium is investigated by the boundary integral equation (BIE) method. The formulations of BIE are derived with the help of generalized Green's formula. The discretization of BIE is based upon constant elements. After confirmation of the accuracy of the present method, some numerical examples are given for various cavities in a full space, in which an isotropic body with a circular cylinder hole is used for comparison and good agreement is observed. It has been proved that the method developed in this paper is effective.     

各向异性体内含任意孔洞对反平面波散射的边界元方法

本文借助于广义格林公式导出了用位移表示的各向异性介质中SH波入射时的边界积分方程.根据本文作者在文献[8]给出的基本解,求解了各向异性介质中孔洞对SH波的散射问题.边界积分方程的离散基于常数元模式.文中给出了一个圆柱、一个椭圆柱和两个椭圆柱形式的孔洞周围的位移场和应力场的数值结果.最后,对入射波频率较高时的情形作了说明.     

CRACK FACE CONTACT PROBLEM ANALYSIS USING THE SCALED BOUNDARY FINITE ELEMENT METHOD

比例边界有限元侧面上有任意荷载时,将侧面载荷分解成关于径向方向局部坐标的多项式函数的和,推导给出了考虑侧面载荷存在的新型形函数,并基于该形函数推导了刚度矩阵和等效节点载荷列阵.首次对比例边界有限元法求解裂纹面接触问题进行了研究,运用Lagrange乘子引入接触界面约束条件,推导给出了比例边界有限元求解裂纹面接触问题的控制方程.将裂纹面单元分为非裂尖单元和含有侧面的裂尖单元.在非裂尖单元中的裂纹面,裂纹面作为多边形单元的边界,边界上的接触力可等效到节点上,通过在节点上构造Lagrange乘子,采用点对点接触约束进行处理.对于含有侧面的裂尖单元,在整个侧面上构造Lagrange乘子的插值场,采用边对边接触约束进行处理.对三个不同的接触约束状态下的算例进行了数值计算,通过与解析解及有限元软件ABAQUS计算结果的对比,验证了本文提出的比例边界有限元点对点和边对边接触求解裂纹面接触问题的精确性与有效性.     

A NOVEL BOUNDARY INTEGRAL EQUATION METHOD FOR LINEAR ELASTICITY--NATURAL BOUNDARY INTEGRAL EQUATION

The boundary integral equation (BIE) of displacement derivatives is put at a disadvantage for the difficulty involved in the evaluation of the hypersingular integrals. In this paper, the operators δij and εij are used to act on the derivative BIE. The boundary displacements, tractions and displacement derivatives are transformed into a set of new boundary tensors as boundary variables. A new BIE formulation termed natural boundary integral equation (NBIE) is obtained. The NBIE is applied to solving two-dimensional elasticity problems. In the NBIE only the strongly singular integrals are contained. The Cauchy principal value integrals occurring in the NBIE are evaluated. A combination of the NBIE and displacement BIE can be used to directly calculate the boundary stresses. The numerical results of several examples demonstrate the accuracy of the NBIE.     

流体流动的边界滑移问题研究进展

最近十几年来,随着现代微/纳米测试以及分子动力学模拟技术的出现和发展,流体流动的边界滑移问题研究获得了突 飞猛进的发展.边界滑移相关研究大体可分为3个方面: 实验、分子动力学模拟和理论数值分析,前两者主要以发现边界滑移现象、探索边界滑移的产生机理以及各因素对边 界滑移的影响规律为主要研究目的,而后者主要研究边界滑移的物理模型、相关问题的计算方法以及边界滑移对流体系 统流体动力学行为的影响.本文首先简要回顾了液体流动的边界滑移及其相关问题的早期研究历史,随后对边界滑移问题 的研究现状进行了综述,最后展望了该领域今后的研究重点及其应用前景.     

二维弹性接触问题中随机边界元法与可靠性分析

发展了二维弹性接触问题中的随机边界元法，推导并建立了相应的随机边界元基本方程，并将所发展的方法用于静强度的可靠性分析，讨论了其数值解技术。通过算例分析表明，本文发展的方法是可行的。     

Boundary element analysis of contact problems using artificial boundary node approach

An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.     

部分滑移边界条件下自由液面流动问题研究

数值方法进行相关问题的研究就以圆柱形储液罐为例, 考虑部分滑移边界条件, 对所得的高 度非线性微分方程进行了数值求解, 讨论了贮液腔体内液面接触线和液体高度线, 及部分滑 移条件对流体层微观半径的影响; 所得结论对自由液面晃动问题数值研究中的边界条件处理 有参考意义.     

Two-dimensional polynomial eigenstrain formulation of boundary integral equation with numerical verification

The low-order polynomial-distributed eigenstrain formulation of the boundary integral equation (BIE) and the corresponding definition of the Eshelby tensors are proposed for the elliptical inhomogeneities in two-dimensional elastic media. Taking the results of the traditional subdomain boundary element method (BEM) as the control, the effectiveness of the present algorithm is verified for the elastic media with a single elliptical inhomogeneity. With the present computational model and algorithm, significant improvements are achieved in terms of the efficiency as compared with the traditional BEM and the accuracy as compared with the constant eigenstrain formulation of the BIE.     

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations （BIE） and solved with the newly developed boundary point method （BPM）. The model is closely derived from the concept of the equivalent inclusion Of Eshelby tensors. Eigenstrains are iteratively determined for each short-fiber embedded in the matrix with various properties via the Eshelby tensors, which can be readily obtained beforehand either through analytical or numerical means. As unknown variables appear only on the boundary of the solution domain, the solution scale of the inhomogeneity problem with the model is greatly reduced. This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM. The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element （RVE）, showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations(BIE)and solved with the newly developed boundary point method(BPM).The model is closely derived from the concept of the equivalent inclusion of Eshelby tensors.Eigenstrains are iteratively determined for each short.fiber embedded in the matrix with various properties via the Eshelby tensors,which can be readily obtained beforehand either through analytical or numerical means.As unknown variables appear only on the boundary of the solution domain,the solution scale of the inhomogeneity problem with the model is greatly reduced.This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM.The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element(RVE),showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

应用边界元方法的正齿轮接触应力分析

本文应用边界元方法,采用在齿廓上按赫兹分布的载荷,直接计算了在接触区附近齿表面和齿体内的真实接触应力。并通过计算定量地分析了各齿轮参数及齿面摩擦对接触应力的影响。为高精度的齿轮设计、寿命预测和渗碳层厚度选择提供了可靠的依据。