Application of an improved boundary element method on the problem of elastic torsion

An improved boundary clement method has been used in analyzing and calculating the problems of the torsion of a prismatic bar with elliptical cross-section. In this paper the calculated results correspond with the values of boundary element method. However, the quantity of data required by the improved boundary element method is much less than that required by boundary element method, and the calculating time will be greatly reduced. Therefore, the procedure of this paper is an economical and efficient numerical computational way for solving Poisson equation problem.     

确定裂纹体等效弹性模量的边界元方法

采用边界元方法计算含有序分布裂纹的裂纹体在压缩载荷作用下的等效弹性模量，利用一种能适当考虑裂纹有间相作用的自洽理论，建立了相应的迭代格式，通过算例研究了裂纹方向，裂纹面间摩擦系数对裂纹体等效弹性模量的影响。     

不同材料契合弹性力学问题的虚边界元最小二乘法

本文针对不同材料契合弹性力学问题,由虚边界元方法出发,建立了引入拉氏乘子的最小二科解法,对该类问题避免了采用Ｈｅｔｅｎｙｉ’ｓ基本解的麻烦和限制。数值结果表明本文算法的有效性和计算精度高。     

The boundary element method applied to elastic contact problems with friction

In the current paper the boundary integral equations (BIE) for elastic contact problems with friction are derived from the incremental virtual work principle. After introducing contact conditions of adhesion and slip into BIE all variants on boundary are made to discretize by quadratic isoparametric boundary element. In the current paper not only an auto-increment loading law is presented but also the iterative calculation laws for open, slip and adhesion condition are given. The results of numerical examples are satisfactory.     

The 3-D boundary element method of roller bearing by plate element analogue

The mill roller bearing is made up of an internal ring, middle rolls and an external ring, the analysis of which is a multi-bodies contact problem. In this paper, based on the three-dimensional elastic contact BEM without friction, and using the structural characteristics of roller bearings, middle rolls are described by elastic plate units of different shapes, which is placed on the internal ring. The discontinuous traction problem of the contact element can be dealt with by the traction sub-element method. Therefore the contact problem can be changed from one of multi-bodies to one of two-bodies. As an elastic deformation of the plate element, it can be calculated by an elastic contact formula and is regarded as an equivalent gap to be put into the total matrix equality. The face-to-face contact condition of displacement and traction on the element will serve as the convergence judging condition. So convergence solution can be obtained before any geometry inconsistent on the contact element may take place. Therefore, the three-dimensional distribution of load of the mill roller bearing can be achieved. It has been proved that the calculating model of this method is visualized and simple, has fast convergence solution and higher accuracy, making it an effective numerical method of designing and analyzing load characteristics of the mill roller bearing. Project supported by the National Natural Science Foundation of China (50075075).     

The elastic solution of concentrated force acting in orthogonal anisotropic half-plane and constant element fundamental formulae of boundary element method

In this paper, the elastic solutions of concentrated force acting in orthogonal anisotropic half-plane are derived by imaginal method and the formulae of coefficient matrix for constant element are put forward. To solve half-plane problems numerically by BEM, this paper provides the necessary formulae. Because the expressions of fundamental solutions are very simple, the object functions could be obtained for every integral of constant element and higher order element of indirect BEM. Thus, the procedure of integration could be avoided in calculation program.     

A boundary element method for calculating elastic waves scattered by axisymmetric inclusion

A Boundary Element Method (BEM) is described to compute the scattering of elastic waves by an axisymmetric inclusion in an infinite elastic medium. The boundary loads applied to the inclusion is expanded in terms of Fourier series in an infinite space. The boundary integral equation is solved in the general direction of the axisymmetric inclusion by BEM. The problem of the 3-D scattering of elastic waves is reduced to a 1-Done. According to the geometric features of the axisymmetric in clusion the ring shell elements are adopted in this method. A comparison is made with other BEM methods. The numerical results show this method can reduce the amount of calculation and enhance the speed of convergence. Supported by Foundation of Ph. D Program of State Education Commission of China     

风沙影响下铁路道碴变形模量的离散元数值分析

针对风沙区有碴铁路道床的结构特性,对其在细沙贯入下的有效变形模量研究有助于理解风沙影响下的道床动力特性。采用离散单元模型对道碴碎石和细沙颗粒进行数值建模,并对不同含沙率γ下的有效变形模量进行了数值分析。结果表明,在低含沙率下(γ<30%)沙石混合体的变形模量基本保持不变,在高含沙率下(γ>30%),变形模量随含沙率的增加呈线性降低,以上变形模量随含沙量的变化规律与试验结果相一致。基于离散单元模型的数值模拟,在细观尺度上对沙石混合体的力链强度、空间分布及配位数进行了分析,揭示了沙石混合体有效变形模量随含沙量变化的内在机理。本文工作对风沙区有碴铁路道床力学行为的研究具有一定的借鉴意义,有助于促进风沙影响下有碴铁路道床的结构设计和沙害治理。     

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

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

NEW BOUNDARY ELEMENT METHOD FOR TORSION PROBLEMS OF CYLINDER WITH CURVILINEAR CRACKS

The Saint-Venant torsion problems of a cylinder with curvilinear cracks were considered and reduced to solving the boundary integral equations only on cracks. Using the interpolation models for both singular crack tip elements and other crack linear elements, the boundary element formulas of the torsion rigidity and stress intensity factors were given. Some typical torsion problems of a cylinder involving a straight, kinked or curvilinear crack were calculated. The obtained results for the case of straight crack agree well with those given by using the Gauss-Chebyshev integration formulas, which demonstrates the validity and applicability of the present boundary element method.     

Shape sensitivity analysis of elastic structures by multi-domain boundary element method

Based on the direct differentiation of boundary element equations to design variables, an algorithm of the boundary element method which is suitable for the shape sensitivity analysis of serial multi-domain elastic structures is built. By the algorithm, it is not necessary to form and solve the equations for the structure as a whole. Each domain can be solved separately. Two computational examples are presented to demonstrate the accuracy and versatility of the algorithm.This work has been partially supported by National Natural Science Foundation of China     

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

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

损伤拉索的等效弹性模量及其参数分析

损伤拉索会出现线形松弛、应力水平降低的情况,必然会影响拉索的等效弹性模量。本文首先引入损伤程度、位置及范围3个参数,用以描述拉索损伤形态的特征,建立损伤拉索索力和线形计算公式,采用数值方法计算了损伤拉索弦向等效弹性模量精确数值,并和经典的等效弹性模量公式的计算结果进行了比较分析,分析了考虑损伤时两种不同计算方法结果的误差。计算表明,对于500m弦向长度以内的损伤拉索,拉索的弦向长度Lc越大,倾角越小,等效弹性模量的损失越大,并且应用割线模量公式计算的误差也越大,当Lc=500m时,损伤拉索相对误差值在2.5%~4.5%之间。弦向应变越小,等效弹性模量损失越大,弦向应变在[0.001,0.004]内,应用割线模量公式计算的相对误差小于3.5%。损伤程度及损伤范围对引用等效弹性模量公式的误差影响较大,倾角对等效弹模公式相对误差的影响也不容忽视。弦向长度、弦向应变、倾角和损伤程度参数都是通过改变拉索的松弛程度进而影响等效弹性模量的数值以及公式的误差。     

An effective boundary element method for analysis of crack problems in a plane elastic plate

A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.     

求解位势问题的虚边界元法

本文提出了求解位势问题的虚边界元法，建立了位势问题的虚边界元的离散方程式，推导了离散化求系数的积分解析式。该方法与传统边界元法相比具有不存在奇异积分和边界附近精度较高等优点，可用来计算真空静电场、稳定温度场、流体绕流、介质中的渗流等各类位势问题。大量算例均获得了满意的结果。     

Effective elastic properties of doubly periodic array of inclusions of various shapes by the boundary element method

A rectangular cell of known boundary conditions is cut out from a medium containing the doubly periodic array of inclusions. The stress and strain relationship of the rectangular cell is obtained by using the classical boundary element methods. By matching the boundary condition requirements, the effective elastic properties of the doubly periodic array of inclusions can then be calculated. Numerical examples from the sub-domain boundary element method and the single domain boundary element method are compared and discussed. However, the present method cannot be readily extended to domains having circular or curved boundary parts.     

岩石爆破开挖诱发振动的等效模拟方法

岩石爆破开挖过程中爆炸荷载的复杂性、岩石介质模型的多变性、本构方程的多样性、以及群孔爆破模拟计算工作量巨大等诸多因素使得准确模拟爆破诱发的围岩振动存在较大的困难,鉴于此提出一种等效数值模拟方法。根据爆破过程炮孔周围岩石破坏范围的空间分布特征,确定群孔起爆条件下爆炸荷载作用的等效弹性边界;通过分析炸药起爆后炮孔空腔动力膨胀、岩体裂纹扩展、炮孔堵塞物运动以及爆生气体的逸出,计算了等效边界上沿炮孔轴向变化的爆炸荷载。结合瀑布沟水电站1# 尾水隧洞爆破开挖监测,基于动力有限元法运用该等效技术模拟了围岩的质点振动速度。数值模拟与现场实测数据的对比表明该等效模拟技术用于岩石爆破开挖围岩响应的计算是合适的。计算结果同时还显示,等效边界邻近岩体力学参数的选取对计算结果有重要影响。     

Boundary element method for moving and rolling contact of 2D elastic bodies with defects

A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes. The project supported by the National Natural Science Foundation of China (19772025)     

Boundary element method for arbitrary elastic thin shell

In this paper, a boundary element scheme for arbitrary elastic thin shells is elaborated, Based on BEM of 3D linear elasticity and Kirchhoff's hypothesis, boundary integral equations for shells are deduced. As a result, only Kelvin's solution is used, the difficulty in finding a fundamental solution of arbitrary shells is successfully avoided. The project supported by National Natural Science Foundation of China     

A higher-order boundary element method for three-dimensional potential problems

A method of eliminating the singularities involved in boundary element methods for three-dimensional potential problems is presented and the non-singular expressions of integrals on an element on which the singular point is situated are given for linear and quadratic interpolation functions. Numerical examples are compared with analytical solutions to show that the higher-order interpolations have better precision.