共查询到20条相似文献，搜索用时 31 毫秒
1.
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. 相似文献
2.
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. 相似文献
3.
4.
A new criterion to predict crack propagation trajectory in anisotropic rocks with incorporating the concept of T-stress in formulating stress field near the crack tip was developed. The developed criterion along with enrichment functions and interaction integral in the extended finite element method (XFEM) framework made a sophisticated tool in modeling fracturing process in anisotropic media. Numerical results indicated that stress intensity factors considerably depend on orientation of anisotropy axes and ratio of the elastic modulus. The proposed formulation for anisotropic media provides a more accurate prediction of crack propagation trajectory compared with conventional methods, especially in mixed mode conditions. 相似文献
5.
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. 相似文献
6.
Summary This paper presents a systematic procedure for determining fracture toughness of an anisotropic marble using the diametral
compression test (Brazilian test) with a central crack on the discs. Additionally, a novel formulation to increase the accuracy
in Stress Intensity Factor (SIF) calculations using Boundary Element Method (BEM) is applied to determine the stress intensity
factors and the fracture toughness of anisotropic rocks under mixed-mode loading. The numerical results show that the SIFs
for both isotropic and anisotropic problems are in good agreement with those reported by previous authors. The marble with
clear white-black foliation from Hualien (in eastern Taiwan), was selected for the Brazilian tests. Diametral loading was
conducted on the Cracked Straight Through Brazilian Disc (CSTBD) specimens to evaluate their fracture toughness. In addition,
a new fracture criterion was developed to predict pure mode I, pure mode II or mixed mode (I–II) fracture toughness of the
anisotropic marble. The new fracture criterion is based on the examination of mode I, mode II and mixed mode (I–II) fracture
toughness for different crack angles and anisotropic orientation.
Authors’ address: Associate Professor Chao-Shi Chen, Department of Resources Engineering, National Cheng Kung University,
No. 1 Dasyue Rd., East District, Tainan 701, Taiwan 相似文献
7.
By using the numerical method to model the ter-rain effect on the magnetotelluric field,few resultshave been obtained. The finite element method(FEM) was used by Chouteau and Bouchard (1988)and Wannamaker et al .(1986) ,andthe boundary el-ement method (BEM) was used by Xu and Zhou(1997) and Xu (1995) to model 2Dtopographyinflu-ences on magnetotelluric surveys . The BEM methodwas also used to model the 3Dtopographic effect onmagnetotelluric deep sounding (Xu et al .,1997 ;Xu,1995) .In t… 相似文献
8.
三维地形大地电磁场的边界元模拟方法 总被引：1，自引：0，他引：1
提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果. 相似文献
9.
Abstract The calculations of unsteady flow to a multiple well system with the application of boundary element method (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem of parabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time variable. The image function of water head H can be solved by BEM. We derived the boundary integral equation of the transformed variable H and the discretization form of it, so that there is no need to discretize the boundaries of well walls and it becomes easier to solve the groundwater head H by numerical inversion. 相似文献
10.
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. 相似文献
11.
结合边界元法和离散元法,提出一种可以进行计算颗粒内部应力和破碎路径的方法。该方法利用离散元法求解颗粒的相互作用和每个颗粒上的荷载。然后利用边界元法计算颗粒的应力分布,为了实现动态平衡,将颗粒的加速度视为恒定大小的体力。但体力导致边界积分方程中出现域积分,故采用直线积分法将域积分转化为边界积分,以保证边界元法降维的优势。为了提高边界元的计算效率,对于几何形状相似的颗粒,以其中一个颗粒作为模板颗粒,只需要计算模板颗粒在局部坐标系中的系数矩阵,其他相似颗粒可以通过局部和全局坐标系之间的映射获得。在得到应力后,基于Hoek-Brown准则来判断颗粒是否破碎。此外,将破坏路径简化为直线,并采用最小二乘法拟合得到破坏路径。 相似文献
12.
目前隧道及大型地下工程往往在裂隙岩体中开挖,而裂隙与地下空间的距离及裂隙的扩展条件,制约着隧道及地下工程的稳定性。应用能考虑孔洞和裂纹问题的新型边界积分方程与无网格加辽金法结合,建立一种新型的边界无单元法。在该方法中基本的未知量是由边界上的面力和边界上位移密度函数构成的复变量边界函数 。文中应用的边界积分公式和Muskhelishvili的积分公式直接相关。将无网格构造方法引入新型的边界积分方程,建立了新型的边界无单元法。应用该方法详细分析了含隧道和裂纹间相互关系等问题,其数值结果与解析结果吻合很好,说明该方法的正确性和可行性。 相似文献
13.
14.
The paper presents closed‐form solutions for stress and displacement influence functions for stress discontinuity (SD) and displacement discontinuity (DD) elements, for a two‐dimensional plane‐strain elastic, transversely anisotropic medium. The solutions for SD elements are based on Kelvin's problem and for DD elements on the concept of dipoles. Stress and displacement influence functions are derived for the following elements: constant SD, linear SD, constant DD, linear DD, square root DD, parabolic DD, constant DD surface, and linear DD surface elements. The formulations are incorporated into FROCK, a hybridized boundary element method code, and are validated by providing comparisons between the results from FROCK and the finite element code ABAQUS. A limited parametric analysis shows the effects of slight anisotropy on the stress field around the tip of a crack and of the orientation of the crack with respect to the axes of elastic symmetry. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
15.
The complex variables boundary element method (CVBEM) is used to study interaction between a circular opening and fractures originating from its boundary in a piecewise homogeneous plane. A new complex hypersingular equation for piecewise homogeneous media with a circular opening is obtained. The equation is solved using the CVBEM technique with circular and straight boundary elements and polynomial approximations (with square root asymptotics for crack tip elements) for the unknown functions. The algorithm is verified through comparison with known semi‐analytical and numerical solutions that involve interaction between a circular opening and specific systems of cracks or other openings. New numerical results concerning the interaction of the circular opening with the cracks and circular inclusions are presented. The method is applied to an important problem in the petroleum industry: modelling propagation of hydraulic fractures in the vicinity of a borehole. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
16.
The finite‐element (FE) method is used for modeling geotechnical and pavement structures exhibiting significant non‐homogeneity. Property gradients generated due to non‐homogeneous distributions of moisture is one such example for geotechnical materials. Aging and temperature‐induced property gradients are common sources of non‐homogeneity for asphalt pavements. Investigation of time‐dependent behavior combined with functionally graded property gradation can be accomplished by means of the non‐homogeneous viscoelastic analysis procedure. This paper describes the development of a generalized isoparametric FE formulation to capture property gradients within elements, and a recursive formulation for solution of hereditary integral equations. The formulation is verified by comparison with analytical and numerical solutions. Two application examples are presented: the first describes stationary crack‐tip fields for viscoelastic functionally graded materials, and the second example demonstrates the application of the proposed procedures for efficient and accurate simulations of interfaces between layers of flexible pavement. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
17.
C. Fidelibus 《国际地质力学数值与分析法杂志》2007,31(11):1329-1348
This paper describes the essential features of a numerical technique for the simulation of the coupled fluid flow and deformation in a 2D assembly of poroelastic blocks and transmissive fractures. The boundary element method (BEM) is applied to each block to reduce Navier and diffusion equations to a set of integral equations involving block boundary terms, whereas a Galerkin weighted‐residuals finite element method (FEM) is applied to the fracture diffusion equations. In addition, fracture local equilibrium is rendered through spring‐like equations relating the stresses to the relative displacements of the fracture walls. A time‐marching process is implemented leading to an algebraic system where the right‐hand side vector is built based on the collected solutions of the previous time steps. The technique requires the meshing of the fracture network only. The accuracy of the results is adequate even with relatively coarse meshes without the resort to small time steps at the beginning of the simulation. It furnishes outputs that focus only on the salient features of the response. The efficiency of the technique is demonstrated through the illustration of the results of three examples. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
18.
基于边界元法的边坡矢量和稳定分析 总被引：4，自引：0，他引：4
矢量和法物理力学意义明确,计算简单,且能根据边坡当前的应力分布状态合理地评价其整体稳定性状态。其中边坡的应力状态通常是采用有限元法来求解。由于边界元法具有研究问题降阶、离散化带来的误差值仅产生在边界以及计算量小等优点,在工程中得到了广泛应用;对于平面问题,以源点作为原点,以所积分单元的切向和法向为坐标轴建立局部坐标系,对于线性单元可以得到所有积分的解析解。因此,可以得到计算区域内部任意点的场变量的解析解,这就保证了位于边界附近区域场变量的精度。利用边界元法得到二维边坡体内连续的应力分布状态,使用矢量和法对该边坡进行稳定性分析,并且与基于有限元的矢量和法、极限平衡法进行对比分析。边坡圆弧滑面和折线滑面的计算结果表明,基于边界元法得到的矢量和安全系数和基于有限元的矢量和法、极限平衡法基本一致;边界元法对应的矢量和安全系数对边界单元尺寸不敏感。 相似文献
19.
A fully coupled element‐free Galerkin model for hydro‐mechanical analysis of advancement of fluid‐driven fractures in porous media 
下载免费PDF全文
下载免费PDF全文 Hydraulic fracturing (HF) of underground formations has widely been used in different fields of engineering. Despite the technological advances in techniques of in situ HF, the industry uses semi‐analytical tools to design HF treatment. This is due to the complex interaction among various mechanisms involved in this process, so that for thorough simulations of HF operations a fully coupled numerical model is required. In this study, using element‐free Galerkin (EFG) mesh‐less method, a new formulation for numerical modeling of hydraulic fracture propagation in porous media is developed. This numerical approach, which is based on the simultaneous solution of equilibrium and continuity equations, considers the hydro‐mechanical coupling between the crack and its surrounding porous medium. Therefore, the developed EFG model is capable of simulating fluid leak‐off and fluid lag phenomena. To create the discrete equation system, the Galerkin technique is applied, and the essential boundary conditions are imposed via penalty method. Then, the resultant constrained integral equations are discretized in space using EFG shape functions. For temporal discretization, a fully implicit scheme is employed. The final set of algebraic equations that forms a non‐linear equation system is solved using the direct iterative procedure. Modeling of cracks is performed on the basis of linear elastic fracture mechanics, and for this purpose, the so‐called diffraction method is employed. For verification of the model, a number of problems are solved. According to the obtained results, the developed EFG computer program can successfully be applied for simulating the complex process of hydraulic fracture propagation in porous media. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
20.
In this paper, an elastostatic half-plane boundary element method (BEM) formulation was applied to analyze the stress behavior of underground pressure pipes, embedded in two-layer soils. In the use of this method, only the boundary of pipe and interfaces were required to be discretized. In this regard, first, a computer code was prepared based on a multi-region substructuring process in the BEM scheme. Then, the efficiency and applicability of the method as well as the prepared algorithm were verified by solving some practical examples and comparing the results with those of the published works. Finally, a parametric study was done to evaluate the effect of pipe depth and determine the soil stress distribution. The studies showed that the half-plane BEM was in good agreement with the existing solutions and its capability was very favorable for elastostatic problems including semi-infinite domain. It is obvious that this method can be practically used to analyze the geotechnical underground buildings in substituting the full-plane BEM formulation. 相似文献