弹性力学问题的虚边界元－最小二乘法及误差评估

本文从［１］提出的虚边界原理出发,采用最小二乘法建立满足弹性力学问题边界条件的边界积分方程,再用线性虚边界元将其离散化。然后详细地研究了这些离散化的边界积分方程的解折特性。文中引用了误差分析的拉依达（ｐａИＴａ）准则,用来衡量解的误差水平,取得了理想的效果。编制了微机程序,程序中采用高斯积分格式,并考虑了虚,实边界对称条件的具体处理。本文方法不仅可以成功地处理边界条件连续的情况,而且对边界条件不连续的情况也能得出满意的结果。数值算例表明,程序可靠,虚边界变动时算法稳定,具有较高的处理精度。     

应力强度因子计算的边界元法

本文利用非连续元离散边界的积分方程，推导了奇异积分的具体表达式，将非连续边界元和多域缩聚法用于二维弹性断裂应力强度因子计算，得到了合理的计算结果。     

车内声场的数学模型建立

本文首先利用Ｈｅｌｍｈｏｌｔｚ方程和Ｇｒｅｅｎ定理推导出适合多种边界条件的车内声场边界积分方程，然后利用边界元数值分析技术离散方程，得到已知某一封闭空间边界的振动特性求解其内部声压的边界元数学模型。作为验证，本文还对两个实例进行了试验，结果表明边界元计算值与理论值和试验实测值吻合良好。     

板弯曲问题三维边界元分析

本文抛弃以往解板弯曲问题的假设,直接从三维弹性力学微分方程出发,依据三给弹性力学问题的Ｋｅｌｖｉｎ解,应用最小二乘法建立了三维虚边界元法解板弯曲问题的一般方法,文中给出了具有各种约束的矩形板的数值算例,以作为本方法的作用,本文方法与边界元直接法相比,优点在于需处理奇积分,且系数阵是对称的;再者,本文方法是思想简单,且程序实现容易,易于被工程界接受。     

板弯曲问题三维虚边界元分析

本文抛弃以往解板弯曲问题的假设,直接从三维弹性力学微分方程出发,依据三维弹性力学问题的Kelvin解,应用最小二乘法建立了三维虚边界元法解板弯曲问题的一般方法。文中给出了具有各种约束的矩形板的数值算例,以作为本方法的应用。本文方法与边界元直接法相比,优点在于无需处理奇异积分,且系数阵是对称的：再者,本文方法思想简单,且程序实现容易,易于被工程界接受。     

瞬态弹性动力学问题的半解析边界元法

本文用半解析有限元法对边界积分方程作离散化处理,通过引入基本解函数和半解析半离散试函数的二次半解析过程,使三维弹性动力学问题简化为一维数值计算。文中又采用移动边界元法来模拟波在半无限介质中传播的表面积分问题,分析计算了各种瞬态波在介质内传播,绕射及地面运动问题。计算结果表明,半解析边界元法不仅计算精度高,而且工作量大大降低,具有较高的经济效益与应用价值。     

薄壁杆件翘曲剪应力的边界元精确积分解法

用非连续边界元对薄壁杆件的约束扭转进行了分析,推导出了求解边界点二次翘曲函数值的边界积分方程,给出了边界积分方程数值求解时积分计算的精确表达式。数值算例表明:利用边界积分方程方法分析薄壁杆件的约束扭转问题时效率和精度高,同时采用精确积分可以有效的处理"边界层效应"问题。     

地面谐运动激励下任意外形水坝稳态响应的双边界元解

本文基于有限水深带形域势流问题的基本解和二维线弹性力学问题的Kelvin解,建立了坝库系统在谐激励下稳态响应的双边界积分方程.推导过程中,利用了Nardini和Brebbia方法将分布惯性力项的体积分化为相应的边界积分.然后通过边界元离散技术,针对两个不同型式的坝体计算了作用在界面上的水动压力分布,其中一个算例的结果和已有的有限元解作了比较.     

正交各向异性厚板振动分析的边界积分方程法

本文采用正交各向异性厚板静力问题的基本解作为边界积分方程的核函数,利用加权残数法建立了正交各向异性厚板振动分析的边界积分方程。文中详细地讨论了边界积分方程的数值处理过程并给出了若干数值算例以论证本文方法的正确性。      

半平面固支边复变基本解及其应用

本文运用Ｒｉｅｍａｎｎ－Ｓｃｈｗａｒｚ对称原理导出直边界为固支的半无限平面问题的复变基本解。应用所求得的基本解建立了相 的复变边界积分方程法并给出数据算例。     

一种建立分区变分原理的新方法

提出了一种建立弹性理论分区变分原理的新方法。放松了分区交界面上位移、应力连续的条件，证明了弹性理论分区求解体系的微分形式与积分形式的等价关系。本文以微分形式为前提，利用这种等价关系,在统一的构架下,导出了分区广义虚功方程和弹性理论分区变分原理。变分原理是积分形式的一种表现形式。讨论了积分形式的物理含义，提出了广义虚函数的概念。广义虚函数具有任意性、虚拟性。     

平面电磁弹性固体的虚边界元——最小二乘配点法

从电磁弹性固体平面问题的基本方程出发,依据弹性力学虚边界元法的基本思想,利用电磁弹性固体平面问题的基本解,提出了电磁弹性固体平面问题的虚边界元——最小二乘配点法。电磁弹性固体的虚边界元法在继承传统边界元法优点的同时,有效地避免了传统边界元法的边界积分奇异性的问题。由于仅在虚实边界选取配点,此方法不需要网格剖分,并且不用进行积分计算。最后给出了一些具体算例,并和已有的解析解进行了对比,结果表明提出的虚边界元方法有很高的精度。     

Virtual boundary element-integral collocation method for the plane magnetoelectroelastic solids

This paper presents a virtual boundary element—integral collocation method (VBEM) for the plane magnetoelectroelastic solids, which is based on the basic idea of the virtual boundary element method for elasticity and the fundamental solutions of the plane magnetoelectroelastic solids. Besides sharing all the advantages of the conventional boundary element method (BEM) over domain discretization methods, it avoids the computation of singular integral on the boundary by introducing the virtual boundary. In the end, several numerical examples are performed to demonstrate the performance of this method, and the results show that they agree well with the exact solutions. The method is one of the efficient numerical methods used to analyze megnatoelectroelastic solids.     

A new boundary integral equation method for analysis of cracked linear elastic bodies

Abstract

A novel integral equation method is developed in this paper for the analysis of two‐dimensional general anisotropic elastic bodies with cracks. In contrast to the conventional boundary integral methods based on reciprocal work theorem, the present method is derived from Stroh's formalism for anisotropic elasticity in conjunction with Cauchy's integral formula. The proposed boundary integral equations contain boundary displacement gradients and tractions on the non‐crack boundary and the dislocations on the crack lines. In cases where only the crack faces are subjected to tractions, the integrals on the non‐crack boundary are non‐singular. The boundary integral equations can be solved using Gaussian‐type integration formulas directly without dividing the boundary into discrete elements. Numerical examples of stress intensity factors are given to illustrate the effectiveness and accuracy of the present method.     

Regularized boundary integral equations for two-dimensional crack problems in multi-field media

A systematic procedure is followed to develop a set of regularized boundary integral equations for modeling cracks in 2D linear multi-field media. The class of media treated is quite general and includes, as special cases, anisotropic elasticity, piezoelectricity and magnetoelectroelasticity. Of particular interest is the development of a pair of weakly-singular, weak-form integral equations for ‘generalized displacement’ and ‘generalized stress’; these serve as the basis for a weakly-singular symmetric Galerkin boundary element method.     

Error estimation using hypersingular integrals in boundary element methods for linear elasticity

A natural measure of the error in the boundary element method rests on the use of both the standard boundary integral equation (BIE) and the hypersingular BIE (HBIE). An approximate (numerical) solution can be obtained using either one of the BIEs. One expects that the residual, obtained when such an approximate solution is substituted to the other BIE is related to the error in the solution. The present work is developed for vector field problems of linear elasticity. In this context, suitable ‘hypersingular residuals’ are shown, under certain special circumstances, to be globally related to the error. Further, heuristic arguments are given for general mixed boundary value problems. The calculated residuals are used to compute element error indicators, and these error indicators are shown to compare well with actual errors in several numerical examples, for which exact errors are known. Conclusions are drawn and potential extensions of the present error estimation method are discussed.     

Integrationsverfahren zur Galerkin-Randeiementmethode für dreidimensionale Elastizitätsprobleme

The mixed boundary value problem in three-dimensional linear elasticity is solved via a system of singular boundary integral equations. This procedure is an alternative to the finite element method and has the main advantage that expensive volume mesh generation is omitted and only a surface mesh is sufficient. The integral equations are discretized by the Galerkin-type boundary element method, which has essential advantages compared to the widely used collocation method. At present the Galerkin method is almost never used in engineering, because this method leads to an unacceptably high effort for the computation of singular double integrals if traditional integration methods are used. The main result of this paper is a new method for the computation of such singular double integrals. The integration procedure leads to simple regular integrand functions also in the case of curved boundary elements. This result simplifies the implementation of the Galerkin-type boundary element method and makes this method applicable in mechanical engineering. Furthermore, the integration of regular double integrals is explained. Numerical tests for model problems in linear elasticity are discussed. Quadrature and discretization errors are analyzed.     

Numerical stability in regular BEM schemes for elastostatic problems

The regular boundary element method is employed for the static analysis of boundary value problems of elasticity. This method allows one to reduce a given boundary value problem to a system of regular integral equations of the first kind with respect to source functions not located on the boundary. This paper is concerned with the numerical stability analysis of regular boundary element methods. In particular, the existence and stability of approximate solutions for integral equations of the first kind with continuous kernels are discussed. The special regularization technique for treating such a class of integral equations is developed. Numerical examples illustrate proposed algorithms and demonstrate their advantages.     

A new boundary integral equation for notch problem of antiplane elasticity

In this paper the notch problem of antiplane elasticity is discussed and a new boundary integral equation is formulated. In the problem, the distributed dislocation density is taken to be the unknown function. Unlike the usual choice, the resultant force function is taken as the right hand term of the integral equation; therefore, a new boundary integral equation for the notch problem of antiplane elasticity with a weaker singular kernel (logarithmic) is obtained. After introducing a particular fundamental solution of antiplane elasticity, the notch problem for the half-plane is discussed and the relevant boundary integral equation is formulated. The integral equations derived are compact in form and convenient for computation. Numerical examples demonstrated that high accuracy can be achieved by using the new boundary equation.