应力边界元法解平面热弹性问题

本文提出了求解平面热弹性问题的应力边界元法。利用应力法由平面热弹性问题的基本方程出发,简要地叙述了边界积分方程的建立,给出了位移单值条件。这种方法适用于应力边界值问题。作为数值计算例,计算了圆形区域和具有偏心圆孔的圆形区域的热应力,得到了满意的结果。应力边界元法也可应用于平板弯曲问题。     

直角域中椭圆形夹杂对SH波的散射与地震动

利用复变函数法和保角映射法研究了SH波入射到含有椭圆形夹杂的直角域中时的散射问题。首先,利用"镜像"叠加原理构造出满足直角域两个直线边界应力自由条件的等效入射及反射波场,通过将椭圆形夹杂的外域映射为单位圆外域构造出满足直角域自由表面应力自由条件的椭圆形夹杂的散射波场。其次,利用椭圆形夹杂边界处的应力和位移连续条件建立求解未知系数的积分方程组,并通过截断有限项求解。最后,给出了在不同参数条件下直角域水平边界处的地表位移幅值。通过算例可以看出:入射波数、入射角度、夹杂位置、介质参数等参数均对地表位移幅值有影响。     

不连续位移超奇异积分方程法解三维多裂纹问题

本文采用Ｂｅｔｉ互等功定理，导出了三维不连续位移基本解的一般形式，然后以该基本解为核函数，建立了求解三维多裂纹问题的超奇异边界积分方程组，并采用三角形单元变换技术和有限部分积分的方法给出了超奇异积分的数值解法，引入非协调单元处理技术解决了法向不确定的角点问题。最后，由裂纹面间断位移可直接求得裂纹前沿任意点的应力强度因子     

双相材料空间中受剪切载荷作用的平片裂纹问题的超奇异积分方程解法

本文利用三维断裂力学的超奇异积分方程的求解方法，对双相材料空间中垂直于界面的平片裂纹在剪切载荷作用下的问题作了研究，首先使用边元界法，在有限部分积分的意义下的问题归结为一组以裂纹面位移间断（位错）为未知函数的超奇异积分方程，然后使用有限部分的理论，并给出边界元法为其建立了数值法，在此基础上，讨论了用裂纹面位移间断计算应力强度因子的方法，最后以两个典型的平片裂纹问题的应力强度因子进行了计算，其数值结     

粘弹性人工边界的虚位移原理

该文将结构及其近场地基作为动力平衡系统,将在人工边界上的波动分解为自由波和散射波,并将输入地震波动转化为作用于人工边界上的等效荷载以实现波动输入。基于以上假设通过分析结构及其近场地基系统的动力平衡关系和自由场的传播机制,给出了自由场的位移表达式、速度表达式,以及在人工边界上由自由场产生的等效荷载一般表达形式,最后建立了粘弹性人工边界统一的动力学积分弱解形式,同时基于有限元程序自动生成系统(FEPG)开发了粘弹性边界条件元件程序。经过计算验证:该文建立的具有粘弹性人工边界的动力学问题的积分弱解方程粘弹性边界条件元件程序可靠、正确。利用这些元件程序,在前处理中可像加位移或应力边界条件一样简便快捷地施加粘弹性边界条件。     

半空间双相介质垂直界面裂纹附近圆形衬砌和半圆凹陷对SH波的散射EI北大核心CSCD

利用Green函数法、镜像法与多级坐标法,对半空间中半圆凹陷和圆形衬砌对SH波的散射进行分析,得到其稳态响应。利用镜像法得到了满足水平边界应力自由、垂直边界位移与应力连续的波函数解析表达式。根据垂直边界连续性条件,利用"裂纹切割法"和"契合法"建立起求解问题的第一类Fredholm型积分方程,得到了圆形衬砌周边的动应力集中系数与裂纹尖端动应力强度因子的解析表达式。数值算例分析了入射波数、衬砌深度、半圆凹陷大小、裂纹长度等对动应力集中系数、裂纹尖端动应力强度因子与地表位移的影响,并与已有文献进行比较。     

用边界元反分析构造平面残余应力场

基于固有应变概念,采用边界元方法,提出一种反方法构造连续的满足域内自平衡条件的平面残余应力场。考虑到反分析的稳定性,固有应变场用一系列光滑基函数(如多项式和三角函数)近似;为了识别由剪切固有应变引起的残余应力,求出对应于固有应变的位移特解与面力特解,将域内积分用双重互易边界元法转换为边界积分,保持了边界元法的优势;同时导出了灵敏度矩阵的显式表达,以提高反分析的效率。最后给出了两个算例验证方法的可行性。     

求解受钉传载荷正交各向异性板双边缺口边缘裂纹应力强度因子的复变－广义变分解法

本文以单边边缘裂纹二维应力场与位移场展开式为基础，采用广义变分方法研究受钉传载菏各向异性板双边缺口边缘裂纹应力强度因子。首先建立精确满足正文各向异性板基本微分方程、裂纹表面边界条件、钉载孔处位移单值条件与合力平衡条件的应力场和位移场的级数表达式。然后应用广义变分方法满足边界条件从而确定应力强度因子，在变分方程中只存在沿板边界的线积分、计算程序简单，输入数据很少，结果收敛迅速，并与已知结果相当吻合而且所需机时较少。     

半空间双相介质垂直界面裂纹附近圆形衬砌和半圆凹陷对SH波的散射

利用Green函数法、镜像法与多级坐标法,对半空间中半圆凹陷和圆形衬砌对SH波的散射进行分析,得到其稳态响应。利用镜像法得到了满足水平边界应力自由、垂直边界位移与应力连续的波函数解析表达式。根据垂直边界连续性条件,利用"裂纹切割法"和"契合法"建立起求解问题的第一类Fredholm型积分方程,得到了圆形衬砌周边的动应力集中系数与裂纹尖端动应力强度因子的解析表达式。数值算例分析了入射波数、衬砌深度、半圆凹陷大小、裂纹长度等对动应力集中系数、裂纹尖端动应力强度因子与地表位移的影响,并与已有文献进行比较。     

有摩擦弹性接触问题边界元分析的一种新方法

在将切向接触力与切向相对位移的关系表为带有惩罚因子的线性互补形式后,直接利用法向接触力与法向相对位移的互补关系,结合边界元技术,本文给出了一种新的求解有摩擦弹性接触问题的数学规划法。     

A nonlinear complementarity approach for elastoplastic problems by BEM without internal cells

A nonlinear complementarity approach is presented to solve elastoplastic problems by the boundary element method, in which the equations are formulated by stress equations and complementarity function obtained from the plasticity constitutive law. The domain integrals involved are transformed into boundary integrals by radial integration method, using compactly supported radial basis functions. Two numerical examples demonstrate the algorithm’s applicability and effectiveness.     

Boundary element analysis in thermoelasticity with and without internal cells

In this paper, a set of internal stress integral equations is derived for solving thermoelastic problems. A jump term and a strongly singular domain integral associated with the temperature of the material are produced in these equations. The strongly singular domain integral is then regularized using a semi‐analytical technique. To avoid the requirement of discretizing the domain into internal cells, domain integrals included in both displacement and internal stress integral equations are transformed into equivalent boundary integrals using the radial integration method (RIM). Two numerical examples for 2D and 3D, respectively, are presented to verify the derived formulations. Copyright © 2003 John Wiley & Sons, Ltd.     

Boundary-only element solutions of 2D and 3D nonlinear and nonhomogeneous elastic problems

This paper presents a robust boundary element method (BEM) that can be used to solve elastic problems with nonlinearly varying material parameters, such as the functionally graded material (FGM) and damage mechanics problems. The main feature of this method is that no internal cells are required to evaluate domain integrals appearing in the conventional integral equations derived for these problems, and very few internal points are needed to improve the computational accuracy. In addition, one of the basic field quantities used in the boundary integral equations is normalized by the material parameter. As a result, no gradients of the field quantities are involved in the integral equations. Another advantage of using the normalized quantities is that no material parameters are included in the boundary integrals, so that a unified equation form can be established for multi-region problems which have different material parameters. This is very efficient for solving composite structural problems.     

Plane stress and plate bending coupling in BEM analysis of shallow shells

In this paper the shear deformable shallow shells are analysed by boundary element method. New boundary integral equations are derived utilizing the Betti's reciprocity principle and coupling boundary element formulation of shear deformable plate and two‐dimensional plane stress elasticity. Two techniques, direct integral method (DIM) and dual reciprocity method (DRM), are developed to transform domain integrals to boundary integrals. The force term is approximted by a set of radial basis functions. Several examples are presented to demonstrate the accuracy of the two methods. The accuracy of results obtained by using boundary element method are compared with exact solutions and the finite element method. Copyright © 2000 John Wiley & Sons, Ltd.     

A boundary element solution to elasto-plastic torsion of solids of revolution

A boundary element method for the solution of problems of elastic and elasto-plastic torsion of solids of revolution is proposed. The displacement and stress field produced by a circumferential ring load in an infinite elastic body are derived and used as the fundamental solution to establish the governing integral equations. For the elasto-plastic case, linear boundary elements and bilinear quadrilateral internal cells are used for discretization of these equations. In order to carry out Gaussian integrations along boundary elements and over internal cells, a method ensuring sufficient accuracy for the calculation of singular integrals is proposed. Numerical results for several problems are given.     

Boundary element analysis of reinforced shear deformable shells

In this paper, stiffened shear‐deformable shells are analysed using the boundary element method. Coupled boundary integral equations are presented for describing curved shells under general loading conditions. The equations are based on boundary integral equations for plane stress and plate bending, with coupling terms arising from the curvature of the shell. Domain integrals are transformed into boundary integrals using the dual reciprocity technique. Stiffeners are modelled as curved beams, continuously attached to the shell. Numerical solutions calculated using the present method are compared with finite element results in two examples. Copyright © 2002 John Wiley & Sons, Ltd.     

Buckling analysis of shear deformable plates by boundary element method

In this paper, the derivation and numerical implementation of boundary integral equations for the buckling analysis of shear deformable plates are presented. Plate buckling equations are derived as a standard eigenvalue problem. The formulation is formed by coupling boundary element formulations of shear deformable plate and two dimensional plane stress elasticity. The eigenvalue problem of plate buckling yields the critical load factor and buckling modes. The domain integrals which appear in this formulation are treated in two different ways: initially the integrals are evaluated using constant cells, and next, they are transformed into equivalent boundary integrals using the dual reciprocity method (DRM). Several examples with different geometry, loading and boundary conditions are presented to demonstrate the accuracy of the formulation. Copyright © 2004 John Wiley & Sons, Ltd.     

A two-dimensional time-domain boundary element method for dynamic crack problems in anisotropic solids

A time-domain boundary element method (BEM) for transient dynamic crack analysis in two-dimensional, homogeneous, anisotropic and linear elastic solids is presented in this paper. Strongly singular displacement boundary integral equations (DBIEs) are applied on the external boundary of the cracked body while hypersingular traction boundary integral equations (TBIEs) are used on the crack-faces. The present time-domain method uses the quadrature formula of Lubich for approximating the convolution integrals and a collocation method for the spatial discretization of the time-domain boundary integral equations. Strongly singular and hypersingular integrals are dealt with by a regularization technique based on a suitable variable change. Discontinuous quadratic quarter-point elements are implemented at the crack-tips to capture the local square-root-behavior of the crack-opening-displacements properly. Numerical examples for computing the dynamic stress intensity factors are presented and discussed to demonstrate the accuracy and the efficiency of the present method.     

A boundary element method without internal cells for solving viscous flow problems

In this paper, a new boundary element method without internal cells is presented for solving viscous flow problems, based on the radial integration method (RIM) which can transform any domain integrals into boundary integrals. Due to the presence of body forces, pressure term and the non-linearity of the convective terms in Navier–Stokes equations, some domain integrals appear in the derived velocity and pressure boundary-domain integral equations. The body forces induced domain integrals are directly transformed into equivalent boundary integrals using RIM. For other domain integrals including unknown quantities (velocity product and pressure), the transformation to the boundary is accomplished by approximating the unknown quantities with the compactly supported fourth-order spline radial basis functions combined with polynomials in global coordinates. Two numerical examples are given to demonstrate the validity and effectiveness of the proposed method.     

The use of Stokes' fundamental solution for the boundary only element formulation of the three-dimensional Navier–Stokes equations for moderate Reynolds numbers

This paper presents a boundary element formulation for the permanent Navier–Stokes equations in which the well-known closed-form fundamental solution for the steady Stokes equations is employed. In this way, from the integral representation formulae for the Stokes' equations, an integral equation is found in which the original non-linear convective terms of the Navier–Stokes equations appear as a domain integral. Additionally, the method of dual reciprocity is used to transform the domain integral to boundary integrals (this method is closely related to the method of particular integrals also used in the literature to transform domain integrals to boundary integrals). Numerical results are presented for the three-dimensional internal flow in a cylindrical container with a rotating cover, in which the accuracy of the method is shown.