Evaluating the stress intensity factors of anisotropic bimaterials using boundary element method

This paper presents a boundary element method (BEM) procedure for a linear elastic fracture mechanics analysis in two‐dimensional anisotropic bimaterials. In this formulation, a displacement integral equation is only collocated on the uncracked boundary, and a traction integral equation is only collocated on one side of the crack surface. A fundamental solution (Green's function) for anisotropic bimaterials is also derived and implemented into the boundary integral formulation so that except for the interfacial crack part, the discretization along the interface can be avoided. A special crack‐tip element is introduced to capture the exact crack‐tip behavior. A computer program using FORTRAN has been developed to effectively calculate the stress intensity factors of an anisotropic bimaterial. This BEM program has been verified to have a good accuracy with previous studies. In addition, a central cracked bimaterial Brazilian specimen constituting cement and gypsum is prepared to conduct the Brazilian test under diametral loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation. Copyright © 2007 John Wiley & Sons, Ltd.     

Boundary element analysis of stresses in an axisymmetric soil mass undergoing consolidation

A time domain boundary element method (BEM) for evaluating stresses in an axisymmetric soil mass undergoing consolidation has been developed. Previous BEM work on axisymmetric poroelasticity for boundary displacements and pore pressures is extended to permit the computation of stresses at both boundary and interior points. The stress formulation preserves the surface-only discretization. The boundary displacement integral equation is progressively differentiated to obtain the related stress and strain integral equations. Explicit expressions for the steady-state axisymmetric fundamental solutions are derived in this process. The transient components of the integrands are obtained directly from the transformation of the three-dimensional kernels into a cylindrical system. Numerical implementation of these integral equations is carried out within a general purpose BEM computer code and several illustrative examples are presented to validate the method.     

Modeling crack propagation path of anisotropic rocks using boundary element method

In a cracked material, the stress intensity factors (SIFs) at the crack tips, which govern the crack propagation and are associated with the strength of the material, are strongly affected by the crack inclination angle and the orientation with respect to the principal direction of anisotropy. In this paper, a formulation of the boundary element method (BEM), based on the relative displacements of the crack tip, is used to determine the mixed‐mode SIFs of isotropic and anisotropic rocks. Numerical examples of the application of the formulation for different crack inclination angles, crack lengths, and degree of material anisotropy are presented. Furthermore, the BEM formulation combined with the maximum circumferential stress criterion is adopted to predict the crack initiation angles and simulate the crack propagation paths. The propagation path in cracked straight through Brazilian disc specimen is numerically predicted and the results of numerical and experimental data compared with the actual laboratory observations. Good agreement is found between the two approaches. The proposed BEM formulation is therefore suitable to simulate the process of crack propagation. Copyright © 2008 John Wiley & Sons, Ltd.     

Fracture mechanics analysis of multiple cracks in anisotropic media

This paper presents a single‐domain boundary element method (BEM) for linear elastic fracture mechanics analysis in the two‐dimensional anisotropic material. In this formulation, the displacement integral equation is collocated on the un‐cracked boundary only, and the traction integral equation is collocated on one side of the crack surface only. A special crack‐tip element was introduced to capture exactly the crack‐tip behavior. A computer program with the FORTRAN language has been developed to effectively calculate the stress intensity factors of an anisotropic material. This BEM program has been verified having a good accuracy with the previous researches. Furthermore, by analyzing the different anisotropic degree cracks in a finite plate, we found that the stress intensity factors of crack tips had apparent influence by the geometry forms of cracks and media with different anisotropic degrees. Copyright © 2010 John Wiley & Sons, Ltd.     

Coupling of BEM with a large displacement and rotation algorithm

In stability analysis of rock blocks, the deformability of the blocks can conveniently be simulated using the boundary element method (BEM). However, all boundary conditions are given as stresses. Thus, the displacement solution is not unique. In this paper, an algorithm is proposed to remove rigid body motions in the solution of the boundary form of Somigliana identity discretized by the direct BEM formulation. The algorithm is applied to the calculation of the normal stiffness of rock blocks and coupled with BS3D, large displacement and rotation algorithm for the general stability of rock blocks. Copyright © 2010 John Wiley & Sons, Ltd.     

基于B-DEM的颗粒聚合体应力计算和破碎路径模拟

结合边界元法和离散元法,提出一种可以进行计算颗粒内部应力和破碎路径的方法。该方法利用离散元法求解颗粒的相互作用和每个颗粒上的荷载。然后利用边界元法计算颗粒的应力分布,为了实现动态平衡,将颗粒的加速度视为恒定大小的体力。但体力导致边界积分方程中出现域积分,故采用直线积分法将域积分转化为边界积分,以保证边界元法降维的优势。为了提高边界元的计算效率,对于几何形状相似的颗粒,以其中一个颗粒作为模板颗粒,只需要计算模板颗粒在局部坐标系中的系数矩阵,其他相似颗粒可以通过局部和全局坐标系之间的映射获得。在得到应力后,基于Hoek-Brown准则来判断颗粒是否破碎。此外,将破坏路径简化为直线,并采用最小二乘法拟合得到破坏路径。     

Some aspects of coupled FEBEM analysis of underground openings

The finite element method (FEM) and the boundary element method (BEM) are two well established numerical methods used for the analysis of underground openings. The advantages of both the methods are utilized by adopting FEBEM in which finite elements are coupled with boundary elements. A coupling procedure is presented in this paper. In using FEBEM, the effect of the location of interface boundary between finite element and boundary element regions, effect of Poisson's ratio and effect of stress ratio are discussed. It is shown that Poisson's ratio and stress ratio have significant effect on the accuracy of the results. Different discretization schemes are discussed to study their effect on accuracy and computation time. The use of different material properties in the FE region is presented. A comparative study is made with FEM for all the cases. It is shown that use of FEBEM is more advantageous than FEM.     

The integration of the Melan solution over a line segment—implementation for boundary element modelling of surface excavations

In this paper, we have deduced the integration solution of the Melan problem (concentrated force at a point in the infinite half-plane) over a line segment directly from the two-dimensional result by the complex variable function method. The BEM dependent on these integration solutions is shown to be valuable for modelling surface and near-surface excavations. Examples are given of the comparison with Mindlin and analytic solutions.     

A direct boundary element method for plane strain poroelasticity

This paper presents a direct boundary element method (BEM), formulated in the Laplace transform space, for plane strain poroelasticity. The paper expands on work by Cheng and Liggett by recasting the theoretical foundation of BEM within the framework of Rice and Cleary's formulation of the Biot theory of poroelasticity. Furthermore, the numerical algorithm is generalized to deal with both interior and exterior domain problems, and a method for indirectly calculating the Cauchy principal value of the singular integrals is presented. Formulae for the stress and flux inside the domain are also derived. Finally, the algorithm is validated by comparing the numerical results with the analytical solution of a borehole subject to a far-field deviatoric stress (exterior domain) and with the solution of Mandel's problem (interior domain). These two examples provide a critical test of the algorithm.     

基于边界元法的边坡矢量和稳定分析

矢量和法物理力学意义明确,计算简单,且能根据边坡当前的应力分布状态合理地评价其整体稳定性状态。其中边坡的应力状态通常是采用有限元法来求解。由于边界元法具有研究问题降阶、离散化带来的误差值仅产生在边界以及计算量小等优点,在工程中得到了广泛应用;对于平面问题,以源点作为原点,以所积分单元的切向和法向为坐标轴建立局部坐标系,对于线性单元可以得到所有积分的解析解。因此,可以得到计算区域内部任意点的场变量的解析解,这就保证了位于边界附近区域场变量的精度。利用边界元法得到二维边坡体内连续的应力分布状态,使用矢量和法对该边坡进行稳定性分析,并且与基于有限元的矢量和法、极限平衡法进行对比分析。边坡圆弧滑面和折线滑面的计算结果表明,基于边界元法得到的矢量和安全系数和基于有限元的矢量和法、极限平衡法基本一致;边界元法对应的矢量和安全系数对边界单元尺寸不敏感。     

A quadratic variation indirect boundary element method for traction boundary-value problems of two-dimensional elastostatics

Boundary integral equations for traction boundary-value problems of two-dimensional elastostatics are derived by the indirect boundary element method. Quadratic variation functions for the representation of geometry, fictitious forces and displacements over each boundary element are described. A system of equations approximating to the boundary integral equations is obtained by a Galerkin formulation in which the integral equation is written at Gauss integration points of elements. The method of computation of the Cauchy principal value is described. Examples of application to the analysis of stress and displacement around underground excavations demonstrate the accuracy and efficiency of the formulation.     

Simulation of sequential excavation with the Boundary Element Method

The paper deals with the application of the Boundary Element Method (BEM) to problems where the boundary geometry and conditions change with time. An example of this is the simulation of sequential excavation in underground construction. A novel approach is presented that is more efficient than the currently used methods. The method involves only one region, where the boundary conditions and the geometry are changed as time progresses. Test examples are presented and results are compared with the ones obtained from other methods. The examples show that the new method is more efficient than currently used methods.     

Elastoplastic deformation around underground openings under biaxial initial stress field

Estimation of elastoplastic deformation around an underground opening induced by the excavation of it, especially displacement and strain field in plastic region, is presented in this paper, as well as the formulation for calculating the displacement and strain in the plastic region around the underground opening by the coupled Boundary Element Method - Characteristics Method (BEM-CM). In this method, the non-associated flow rule is adopted to calculate the displacement and strain field in the plastic region, which is determined by the integration of the displacement along characteristics lines under the boundary condition of the elastic displacement on an elastoplastic interface analysed. It is shown that this method is one of the accurate and effective methods for estimating not only the shape and extent of the plastic region but also the state of the displacement and strain in the plastic region around the underground opening, comparing the theoretical solution with numerical results by this method for a circular opening under hydrostatic initial stress condition. Furthermore, this method is applied to rectangular and horse-shoe shaped openings and the characteristics of the strain field in the plastic region are discussed. © 1997 by John Wiley & Sons, Ltd.     

露天地下开采隔离层稳定性研究

露天地下开采隔离层稳定性分析是矿山中经常遇到的问题。以大宝山矿露天地下开采的工程实例,采用数值分析方法,对隔离层稳定性进行分析计算,发现空区顶板的拉应力是关系隔离层稳定性的关键因素。通过对隔离层的安全厚度数值计算,与5种理论计算方法结果进行了比较,综合求和归一法数据处理和多项式数值逼近得到了不同空区跨度与隔离层安全厚度关系。结果为露天地下开采隔离层稳定性分析的方法进行了验证和补充,也为矿山设计隔离层厚度提供了参考,对指导露天地下开采安全生产施工具有重要意义。     

Scattering of a plane wave by a lined cylindrical cavity in a poroelastic half-plane

The focus of this contribution is to develop a complex variable function method to solve the two-dimensional scattering of plane waves by a lined cylindrical cavity in the poroelastic half-plane. The poroelastic half-plane is based on Biot’s dynamic theory, and the governing equations are solved by reduction to three Helmholtz equations. The lining structure can be treated as an elastic material and decoupled into two Helmholtz equations. Here, the large circle assumption is applied to simulate the half-plane boundary. By using appropriate boundary conditions and continuity conditions, the unknown coefficients in the potentials can be determined. Selected numerical results are presented in this paper.     

基于间接边界元方法的SH波倾斜入射下边坡动力响应特征

地震作用下山体边坡的动力响应规律与地震波入射角度显著相关,本文基于间接边界元方法研究SH波倾斜入射下边坡的动力响应特征.基于间接边界元理论研发了用于边坡岩体动力响应分析的边界元程序,通过计算平直裂隙对弹性波的反射及散射波波场对程序进行验证,并详细探讨了不同类型的边坡在各种入射角度的应力波作用下的响应特性.SH波铅直向入射时,半圆形凹陷地形的谷肩位置地震动振幅最小,而半圆形凸起地形的最高点位置振幅最大,单面坡的坡肩位置地震动最大,楔形凸起地形的坡顶地震动最强烈;当SH波入射角增大,各类边坡的最大地震动力响应位置也有所偏移,地震放大系数亦随之变化;以樟木镇边坡为实例揭示了复杂形态边坡岩体中的局部凸起对地震动的放大作用.研究结果表明,边坡微地形特征对地震动力响应影响非常大,凸起山体坡顶的动力放大效应最为显著;同一边坡在不同入射角地震波作用下产生的动力放大系数以及发生的位置不同.本文的研究结果对评价坡体稳定性和边坡工程抗震防设等具有重要的指导意义.      

Elasto-plastic analysis of an underground opening by FEM and coupled FEBEM

The elasto-plastic analysis of a circular underground opening using the finite element method (FEM), the coupled finite element and the boundary element method (FEBEM) is presented. The coupling procedure and elasto-plastic formulation are discussed. The effect of in situ stress ratio on yielded zone, displaced shape and principal stresses is presented. The results of FEBEM analysis are compared with those obtained from FEM. The computation time and the number of iterations as required by FEBEM and FEM are compared.     

新型边界积分方程的无单元解及其应用

目前隧道及大型地下工程往往在裂隙岩体中开挖,而裂隙与地下空间的距离及裂隙的扩展条件,制约着隧道及地下工程的稳定性。应用能考虑孔洞和裂纹问题的新型边界积分方程与无网格加辽金法结合,建立一种新型的边界无单元法。在该方法中基本的未知量是由边界上的面力和边界上位移密度函数构成的复变量边界函数 。文中应用的边界积分公式和Muskhelishvili的积分公式直接相关。将无网格构造方法引入新型的边界积分方程,建立了新型的边界无单元法。应用该方法详细分析了含隧道和裂纹间相互关系等问题,其数值结果与解析结果吻合很好,说明该方法的正确性和可行性。     

基于边界元方法的次级断裂信息挖掘试验研究

为研究次级断裂发育状况,探讨了一种地震解释主断层与地质力学相结合的边界元模拟方法：首先在安德森断层分类体系下,依据断层与应力组合特征,完善次级断裂构造演化事件分类;然后应用应力张量降阶、线性叠加原理,简化应力与位移平衡方程,并在此基础上使用蒙特卡洛方法反演古应力场;最后求解应力、应变与不连续位移平衡方程,获取地层现今应力分布,分析次级断裂开度、相对密度等信息。运用该方法在杨柳矿区童亭背斜区域开展试验,结果表明区域内最小水平主应力以近北西向为主,主断层之间广泛发育交叉断裂带,次级断裂发育带连通性强,为杨柳矿区隐蔽灾害评估提供了参考资料。     

地面磁梯度测量在勘查地下污水管道中的应用

利用在辽宁某地区实测的垂直磁梯度数据,采用欧拉反褶积方法对地下异常情况进行处理和解释。正演研究表明,理论计算获得的磁梯度比实测磁总场的灵敏度高,结合磁梯度数据和欧拉反演方法可以进行较为复杂地形的解释。通过将计算获得的磁梯度异常与实测磁梯度异常对比,后者对地质体的分辨率更高。对地面垂直磁梯度数据进行反演,结果表明:该方法能够准确地确定地下污水管道的边界和埋深。