弹性力学问题的局部 Petrov-Galerkin方法

提出了弹性力学平面问题的局部Petrov-Galerkin方法,这是一种真正的无网格方法。这种方法采和移动最小二乘近似函数作为试函数,并且采用移动最小二乘近似函数的权函数作为加权残值法加权函数;同时这种方法只包含中心在所考虑点处的规则局部区域上以及局部边界上的积分,所得系统矩阵是一个带状稀疏矩阵,该方法可以容易推广到求解非线性问题以及非均匀介质的力学问题。还计算了两个弹性力学平面问题的例子,给出了位移和能量的索波列夫模及其相对误差。所得计算结果证明：该方法是一种具有收敛快、精度高、简便有效的通用方法;在工程中具有广阔的应用前景。     

弹性力学问题中一个新的边界积分方程——自然边界积分方程

位移导数边界积分方程一直存在着超奇异积分计算的障碍.该文提出以符号算子δij和εij作用于位移导数边界积分方程,施用一系列变换将边界位移、面力和位移导数转成为新的边界张量,从而得到一个新的边界积分方程——自然边界积分方程.自然边界积分方程的奇异性为强奇性,文中给出了相应的Cauchy主值积分算式.自然边界积分方程与位移边界积分方程联合可直接获取边界应力.几个算例表明了自然边界积分方程的正确性.     

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

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

弹性力学边界积分方程的一个发展

本文讨论二维弹性力学平面问题，独立于Ｒｉｚｚｏ型边界分方程，一类新型的边界积分方程，其边界场变量包含应力分量σｉｊｔｉｔｊ（其中ｔｉ是边界切向余弦）。该应力分量可直接用数值方法解边界积分方程求出，它比常规的边界元解提高一阶精度。文末的算例表明确定论的实用性和有效性。     

平面问题等价边界积分方程的三次边界轮廓法

基于弹性力学平面问题等的边界积分方程,给出了三次单元的边界轮廓法。根据平面问题解的复变函数表示,构造了三次形函数。给出了对于混合边值问题求解系统方程确定的边界轮廓方程配置和三次单元界轮廓法的实施。     

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

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

弹性力学平面问题的等价边界积分方程的边界轮廓法

基于边界积分方程中被积函数散度为零的特性,提出了弹性力学平面问题的等价边界积分方程的边界轮廓法,该方法无需进行数值积分,只需要计算单元两结点势函数值之差。实例计算说明,基于传统的边界积分方程的边界轮廓法所得到的面力结果是错误,而本文建立的边界轮廓法则可给出精确的结果。     

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

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

平面弹性力学充要的随机边界元

将平面弹性力学确定性的充分必要的边界积分方程推广到含材料常数随机的不确定问题中去，给出了位移的均值以及偏差的充分必要的边界积分方程。数值计算结果表明，和确定性的积分方程一样，习用的随机边界积分方程在退化尺度附近，无论是均值还是偏差都存在巨大的误差，而充要的随机边界积分方程则始终保持良好的精度     

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

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

A meshless method of a thin plate

A meshless approach to analysis of arbitrary Kirchhoff plates by the local boundary integral equation(LBIE) method is presented. The method combines the advantageous features of, all the three methods: the Galerkin finite element method (GFEM), the boundary element method (BEM) and the element-free Galerkin method (EFGM). It is a truly meshless method, which means that the discretization is independent of geometric subdivision into elements or cells, but is only based on a set of nodes (ordered or scattered) over a domain in question. It involves only boundary integration, however, over a local boundary centered at the node in question; It poses no difficulties in satisfying the essential boundary conditions while leading to banded and sparse system matrices using the moving least square (MLS) approximations. It is shown that high accuracy can be achieved for arbitrary geometries for clamped and simply-supported edge conditions. The method is found to be simple, efficient, and attractive. Project supported by the National Science Foundation of China (No. 19972019).     

Improved non-singular local boundary integral equation method

When the source nodes are on the global boundary in the implementation of local boundary integral equation method (LBIEM),singularities in the local boundary integrals need to be treated specially. In the current paper,local integral equations are adopted for the nodes inside the domain trod moving least square approximation (MLSA) for the nodes on the global boundary,thus singularities will not occur in the new al- gorithm.At the same time,approximation errors of boundary integrals are reduced significantly.As applications and numerical tests,Laplace equation and Helmholtz equa- tion problems are considered and excellent numerical results are obtained.Furthermore, when solving the Hehnholtz problems,the modified basis functions with wave solutions are adapted to replace the usually-used monomial basis functions.Numerical results show that this treatment is simple and effective and its application is promising in solutions for the wave propagation problem with high wave number.     

A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless method is a new numerical technique presented in recent years .It uses the moving least square (MLS) approximation as a shape function . The smoothness of the MLS approximation is determined by that of the basic function and of the weight function, and is mainly determined by that of the weight function. Therefore, the weight function greatly affects the accuracy of results obtained. Different kinds of weight functions, such as the spline function, the Gauss function and so on, are proposed recently by many researchers. In the present work, the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method. The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed. Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and a in Gauss and exponential weight functions are in the range of reasonable values, respectively, and the higher the smoothness of the weight function, the better the features of the solutions.     

Research on the companion solution for a thin plate in the meshless local boundary integral equation method

The meshless local boundary integral equation method is a currently developed numerical method, which combines the advantageous features of Galerkin finite element method ( GFEM ), boundary element method (BEM) and element free Galerkin method (EFGM), and is a truly meshless method possessing wide prospects in engineeringapplications. The companion solution and all the other formulas required in the meshless local boundary integral equation for a thin plate were presented, in order to make this method apply to solve the thin plate problem.     

AN IMPROVED HYBRID BOUNDARY NODE METHOD IN TWO-DIMENSIONAL SOLIDS

The hybrid boundary node method （HBNM） is a promising method for solving boundary value problems with the hybrid displacement variational formulation and shape functions from the moving least squares（MLS） approximation. The main idea is to reduce the dimensionality of the former and keep the meshless advantage of the latter. Following its application in solving potential problems, it is further developed and numerically implemented for 2D solids in this paper. The rigid movement method is employed to solve the hyper-singular integrations. Numerical examples for some 2D solids have been given to show the characteristics. The computation results obtained by the present method are in excellent agreement with the analytical solution. The parameters that influence the performance of this method are studied through numerical examples.     

泊松方程中域积分化为边界积分的方法

本文讨论了二维和三维泊松方程中域积分化为边界积分的方法。对于形如x~ig_x(y,z)、y~ig_x(x,z)和z~ig_z(x,y)的荷载给出了域积分转化为边界积分的正确公式。而对于复杂荷载,利用泰勒展开将域积分近似地转化为边界积分并给出了误差估计。计算结果表明利用本文方法可大大节省计算时间。因此,本文方法是一种十分有效的方法。     

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

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