断裂问题的扩展有限元法研究

扩展有限元（extended finite element method,XFEM）是近年来发展起来的、在常规有限元框架内求解不连续问题的有效数值计算方法,其基于单位分解的思想,在常规有限元位移模式中加入能够反映裂纹面不连续性的跳跃函数及裂尖渐进位移场函数,避免了采用常规有限元计算断裂问题时需要对裂纹尖端重新加密网格造成的不便。在推导扩展有限元算法的基础上,分析了应力强度因子的J积分计算方法及积分区域的选取。采用XFEM对I型裂纹进行了计算,有限元网格独立于裂纹面,无需在裂纹尖端加密网格;分析了积分区域、网格密度对应力强度因子计算精度的影响,指出了计算应力强度因子的合适参数,验证了此方法的可靠性和准确性。     

岩体裂隙模拟的扩展有限元法应用研究

岩体中普遍存在着断层﹑节理和裂隙等结构面,这些结构面的存在和发展对岩体的整体强度﹑变形及稳定性有极大的影响。因此,研究岩体中原生结构面的萌生﹑发展以及贯通演化过程对评估岩体工程安全性和可靠性具有非常重要的理论与现实意义。扩展有限元法（XFEM）作为一种求解不连续问题的有效数值方法,模拟裂隙时独立于网格,因此,在模拟岩体裂隙扩展﹑水力劈裂等方面具有独特优势。针对扩展有限元法的基本理论及其在岩体裂隙扩展模拟中的应用展开了研究,建立了扩展有限元法求解岩体裂隙摩擦接触、岩体裂隙破坏等问题的数值模型,并将计算模型应用于岩质边坡稳定性分析和重力坝坝基断裂破坏等工程问题。     

Extended Finite Element Method for the Analysis of Discontinuities in Rock Masses

The strength and deformability of rock mass primarily depend on the condition of joints and their spacing and partially on the engineering properties of rock matrix. Till today, numerical analysis of discontinuities e.g. joint, fault, shear plane and others is conducted placing an interface element in between two adjacent rock matrix elements. However, the applicability of interface elements is limited in rock mechanics problems having multiple discontinuities due to its inherent numerical difficulties often leading to non-convergent solution. Recent developments in extended finite element method (XFEM) having strong discontinuity imbedded within a regular element provide an opportunity to analyze discrete discontinuities in rock masses without any numerical difficulties. This concept is based on partition of unity principle and can be used for cohesive rock joints. This paper summarizes the mathematical frameworks for the implementation of strong discontinuities in 3 and 6 nodded triangular elements and also provides numerical examples of the application of XFEM in one and two dimensional problems with single and multiple discontinuities.     

扩展有限元法的数值方面

扩展有限元法是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法,其原理是在裂尖附近用一些奇异函数和沿裂纹面用阶跃函数加强传统有限元的基,以考虑跨过裂纹的位移场的不连续,该加强策略允许计算网格独立于不连续体几何。讨论了扩展有限元法的一些数值方面,主要包括：水平集法确定界面和加强节点与加强方式、裂尖加强范围的选择、J积分区域的确定和积分方案等。     

模拟三维裂纹问题的扩展有限元法

扩展有限元法是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法,其计算网格与不连续面相互独立,因此模拟移动不连续面时无需对网格进行重新剖分。给出了模拟三维裂纹问题的扩展有限元法。在常规有限元位移模式中,基于单位分解的思想加进一个阶跃函数和二维渐近裂尖位移场,反映裂纹处位移的不连续性。用两个水平集函数表示裂纹。采用线性互补法求解裂纹面非线性接触条件,不需要迭代,提高了计算效率。采用两点位移外推法计算裂纹前缘应力强度因子。给出了3个三维弹性静力问题算例,其结果显示了所提方法能获得高精度的应力强度因子,并能有效地处理裂纹面间的接触问题,同时表明扩展有限元结合线性互补法求解不连续问题具有较好的前景。     

XFEM formulation with sub-interpolation,and equivalence to zero-thickness interface elements

This paper describes a particular formulation of the extended finite element method (XFEM) specifically conceived for application to existing discontinuities of fixed location, for instance, in geological media. The formulation is based on two nonstandard assumptions: (1) the use of sub-interpolation functions for each subdomain and (2) the use of fictitious displacement variables on the nodes across the discontinuity (instead of the more traditional jump variables). Thanks to the first of those assumptions, the proposed XFEM formulation may be shown to be equivalent to the standard finite element method with zero-thickness interface elements for the discontinuities (FEM+z). The said equivalence is theoretically proven for the case of quadrangular elements cut in two quadrangles by the discontinuity, and only approximate for other types of intersections of quadrangular or triangular elements, in which the XFEM formulation corresponds to a kinematically restricted version of the corresponding interface plus continuum scheme. The proposed XFEM formulation with sub-interpolation, also helps improving spurious oscillations of the results obtained with natural interpolation functions when the discontinuity runs skew to the mesh. A possible explanation for these oscillations is provided, which also explains the improvement observed with sub-interpolation. The paper also discusses the oscillations observed in the numerical results when some nodes are too close to the discontinuity and proposes the remedy of moving those nodes onto the discontinuity itself. All the aspects discussed are illustrated with some examples of application, the results of which are compared with closed-form analytical solutions or to existing XFEM results from the literature.     

Kinematic modelling of shear band localization using discrete finite elements

Modelling shear band is an important problem in analysing failure of earth structures in soil mechanics. Shear banding is the result of localization of deformation in soil masses. Most finite element schemes are unable to model discrete shear band formation and propagation due to the difficulties in modelling strain and displacement discontinuities. In this paper, a framework to generate shear band elements automatically and continuously is developed. The propagating shear band is modelled using discrete shear band elements by splitting the original finite element mesh. The location or orientation of the shear band is not predetermined in the original finite element mesh. Based on the elasto‐perfect plasticity with an associated flow rule, empirical bifurcation and location criteria are proposed which make band propagation as realistic as possible. Using the Mohr–Coulomb material model, various results from numerical simulations of biaxial tests and passive earth pressure problems have shown that the proposed framework is able to display actual patterns of shear banding in geomaterials. In the numerical examples, the occurrence of multiple shear bands in biaxial test and in the passive earth pressure problem is confirmed by field and laboratory observations. The effects of mesh density and mesh alignment on the shear band patterns and limit loads are also investigated. Copyright © 2003 John Wiley & Sons, Ltd.     

岩体结构面新型本构模型研究

对于涉及需要单独考虑岩体结构面的工程岩体结构分析,采用能反映岩体结构面主要力学特性的合理的本构模型是取得合理解答的关键问题之一。针对在经典连续介质力学理论框架内建立岩体结构面本构模型的缺点,基于岩体结构面的实际受力变形特性,采用直接法建立了一种新型岩体结构面本构模型。所建立的模型依据岩体结构面切向应力变形曲线及剪胀曲线的实际特征,将其分为峰前线性段、峰前非线性段以及峰后软化段,并分别给出了用于描述岩体结构面变形和强度等主要力学特性的数学模型,进而推导建立了结构面各变形阶段的增量型本构模型。最后,编写相关计算程序,采用所建立的新型本构模型以及被广泛采用的Plesha模型对经典的岩体结构面直剪试验成果进行拟合分析。结果表明,所建立的新型本构模型能更为合理的描述岩体结构面的主要力学特性,且模拟能力优于Plesha模型     

基于PFC2D数值试验的异性结构面剪切强度特性研究

天然岩体中广泛发育两侧岩性不同的异性结构面,开展异性结构面变形和强度特性研究旨在为岩体稳定性评价和利用提供依据。选取三峡库区侏罗系典型的砂岩-泥岩异性岩层,首先运用分形几何理论,定量计算了平直和4种不同不规则起伏形态结构面的粗糙度系数JRC值,然后基于PFC2D颗粒流程序,分别开展了以上5种形态异性结构面的数值剪切试验,获得了各形态结构面在不同正应力下的剪切应力-位移曲线。根据数值试验结果,采用巴顿的JRC-JCS模型分析了异性结构面强度特性,并与同性结构面强度性质进行对比研究。最后,在考虑异性结构面剪切破坏机制的基础上,引入强度因子的概念,提出了新的适用于异性结构面强度评价的两类改进巴顿准则。研究结果表明：异性岩体结构面抗剪强度介于相同粗糙度的两种同性结构面强度之间,在较低正应力下接近软岩同性结构面强度,符合Ⅰ型改进巴顿准则;在较高正应力下偏向硬岩同性结构面强度,符合Ⅱ型改进巴顿准则。实际工程中可利用改进准则并根据异性结构面应力状态对岩体稳定性进行评价,弥补了以往研究的不足。     

The role of non‐coaxiality in the simulation of strain localization based on classical and Cosserat continua

It is normally accepted that materials inside the shear band undergo severe rotation of the principal stress direction, which causes non‐coaxiality between the principal stress and principal plastic strain rate. However, classical plasticity flow theory implicitly assumes that the principal stress and the principal plastic strain rate are coaxial; thus, it may not correctly predict the onset of the shear band. In addition, classical continuum does not contain any internal length scales; as a result, it cannot provide a reasonable shear band thickness. In this study, the original vertex non‐coaxial plastic model based on the classical continuum is extended to the Cosserat continuum. The corresponding codes are implemented via the interface of the user defined element subroutine in ABAQUS. Through a simple shear test, the effectiveness of the user's codes is verified. Through a uniaxial compression test, the influence of non‐coaxiality on the onset, the orientation, and the thickness of the shear band is investigated. Results show that the onset of the shear localization is delayed, and the thickness of the shear band is widened when the non‐coaxial degree increases, while the orientation of the shear band is little affected by the non‐coaxial degree. In addition, it is found that the non‐coaxiality can weaken the micro‐polar effect to some extent; nonetheless, the Cosserat non‐coaxial model still has its advantage over the classical non‐coaxial model in capturing the pre‐bifurcation as well as the post‐bifurcation behaviors of strain localization. Copyright © 2016 John Wiley & Sons, Ltd.     

Frictional slip plane growth by localization detection and the extended finite element method (XFEM)

This paper is concerned with developing a numerical tool for detecting instabilities in elasto‐plastic solids (with an emphasis on soils) and inserting a discontinuity at these instabilities allowing the boundary value problem to proceed beyond these instabilities. This consists of implementing an algorithm for detection of strong discontinuities within a finite element (FE) framework. These discontinuities are then inserted into the FE problem through the use of a displacement field enrichment technique called the extended finite element method (XFEM). The newly formed discontinuities are governed by a Mohr–Coulomb frictional law that is enforced by a penalty method. This implementation within an FE framework is then tested on a compressive soil block and a soil slope where the discontinuity is inserted and grown according to the localization detection. Copyright © 2010 John Wiley & Sons, Ltd.     

含结构面岩体试样单轴强度与变形特征

结构面是岩体区别于岩石材料的一大特征,其产状、迹长、密度等参数对岩体的力学性质有着重要影响。本文利用FLAC3D对含结构面岩体试样的单轴压缩特性进行了较为系统的数值模拟研究。文中建立了含不同组贯通性结构面的岩体试样模型和含不同倾角及迹长的非贯通结构面岩体试样模型,对每个试样进行单轴压缩试验的数值模拟,结构体和结构面均采用Mohr-Coulomb剪切和拉伸破坏准则。模拟中用编制的伺服控制程序通过调节加载速度,控制试样内最大不平衡力,研究含结构面试样单轴压缩情况下的变形、强度及破坏方式等特征。模拟结果显示,含1-3组贯通性结构面试样呈现各向异性特征,而含4组贯通性结构面试件呈现各向同性特征。随着贯通性结构面数量的增多,同尺寸试件的变形强度参数劣化。含单组非贯通性结构面试件,其单轴压缩模拟试验的应力-应变曲线峰值后出现应力降。基于Mohr-Coulomb抗剪强度准则和损伤理论所得的解析解与数值模拟结果所得的非贯通性结构面试件的单轴压缩强度不符,说明用抗剪强度准则与损伤理论刻画非贯通结构面试样的强度并不合理。随着非贯通性结构面贯通率的增大,试件的变形、强度参数劣化。含单组结构面试件的破坏方式可分为结构面控制破坏,结构面部分控制破坏和结构面不控制破坏3种类型,而随着结构面组数的增多,结构面控制试样破坏的概率增加。     

基于岩体不连续面三维分形维岩体质量评价研究

岩体内部不连续面分布具有随机性、不规则性,但又具有统计自相似性。基于三维随机不连续面网络模拟技术,对岩体结构统计均质区内的不连续面进行计算机模拟,并应用分形理论计算岩体不连续面分布的分形维数。岩体不连续面分布分维数具有尺寸效应,随着岩体尺寸增大分维数减小,到一定程度趋于稳定。这一稳定的分维值称之为表征分维数,用此值描述不连续面分布特征。将岩体不连续面分布的表征分维数与按传统岩体分类标准所得岩体质量等级进行对比,提出以岩体不连续面分布表征分维数为指标的岩体质量分类方法。在岩体分类基础上,提出基于表征分维数的岩体等效抗剪强度指标的折减计算。     

扩展有限元法在岩石断裂力学中的应用

岩石裂纹的扩展是一个经典的不连续问题,常规有限元方法难以实现裂纹扩展过程的仿真模拟。扩展有限元法(XFEM)实现了计算网格与不连续面相互独立,因此模拟移动的不连续面时无需对网格进行重新剖分。本文介绍了XFEM基本原理和岩石断裂力学常用判据,尝试对岩石类材料单缝Ⅰ型三点弯曲、单缝剪切和双缝平板实验进行模拟。分析结果表明:扩展有限元模拟岩石类材料断裂问题不受网格划分限制,裂纹以实际应力场分布随机扩展;直观地给出岩样的微裂纹产生、演化,直至完全破坏的全过程,并与实验结果吻合。该方法能够应用到岩石断裂力学方面的研究,模拟岩石类材料的宏细观破坏过程,为解决复杂问题提供了方便的途径。     

The response of an elastic half‐space under a momentary shear line impulse

The response of an ideal elastic half‐space to a line‐concentrated impulsive vector shear force applied momentarily is obtained by an analytical–numerical computational method based on the theory of characteristics in conjunction with kinematical relations derived across surfaces of strong discontinuities. The shear force is concentrated along an infinite line, drawn on the surface of the half‐space, while being normal to that line as well as to the axis of symmetry of the half‐space. An exact loading model is introduced and built into the computational method for this shear force. With this model, a compatibility exists among the prescribed applied force, the geometric decay of the shear stress component at the precursor shear wave, and the boundary conditions of the half‐space; in this sense, the source configuration is exact. For the transient boundary‐value problem described above, a wave characteristics formulation is presented, where its differential equations are extended to allow for strong discontinuities which occur in the material motion of the half‐space. A numerical integration of these extended differential equations is then carried out in a three‐dimensional spatiotemporal wavegrid formed by the Cartesian bicharacteristic curves of the wave characteristics formulation. This work is devoted to the construction of the computational method and to the concepts involved therein, whereas the interpretation of the resultant transient deformation of the half‐space is presented in a subsequent paper. Copyright © 2003 John Wiley & Sons, Ltd.     

Study of evolution of shear band systems in sand retained by flexible wall

The present paper is concerned with the study of the evolution of regular patterns of shear bands observed experimentally (cf. Milligan, 1974), and provides an extension of the previous work published by the same authors. The purpose of this paper is to present an improved version of a simple theoretical model, derived basically from a classical equilibrium of the Coulomb wedge. This model constitutes a modified version of an extension of the classical Coulomb wedge analysis by assuming that soil parameters are varying during the deformation process and the initial configuration at which limit equilibrium occurs evolves toward a new equilibrium configuration. The application of the model in the analysis of shear band pattern observed in dredged model tests on cantilever walls provided realistic simulation of consecutive shear band formation. Copyright © 2001 John Wiley & Sons, Ltd.     

饱和砂土局部变形带模拟的有限元数值实现

基于有限变形理论,推导了Newton-Raphson 迭代算法在k+1步增量表达的矩阵形式,实现了饱和砂土变形局部化的有限元数值计算,得到了饱和砂土发生局部化变形的准则。基于Galerkin 方法,得到了位移场和应力场的空间离散化矩阵方程;由土体局部变形带的连续性条件,引入第1切线算子,推导出了砂土等颗粒状媒介发生局部化变形的必要条件。基于此核心算法,编制了有限元计算程序,模拟了饱和砂土在不排水条件下平面压缩过程中剪切带的形成与发展;通过比较分析,研究了有限元网格粗细对于土体局部变形带的影响,结果表明,网格粗细的病态依赖只是微小的,它只与变形条带的宽度有关,对于土体所表现出来的其他力学特性没有影响。     

Strain parameters of discontinuities in rock for finite element calculation

The bulk strain induced by excavation of rooms for storage of highly radioactive waste is usually unimportant, but the change in aperture and shear displacement of permeable fractures may be of practical importance. Hence, the normal and shear of such discontinuities have to be considered and are commonly predicted by using special deformation moduli termed “normal” and “shear” stiffnesses in numerical calculations. For the use of some finite element methods, these have to be converted to compression and shear moduli ,which, in turn, requires that the thickness of the discontinuity is known. The normal and shear strain can be expressed in terms of soil mechanical parameters, yielding the compression modulus M and the shear modulus G, with the required form for finite element calculations. By definition they are functions of the normal and shear stiffness K_n and K_s. It is concluded that calculation of the normal and shear strain of discontinuities with clastic fillings by the use of finite element methods and deformation moduli derived from stiffness numbers is very uncertain except when the geometry of the weaknesses can be accurately defined.     

Finite volume–based modeling of flow‐induced shear failure along fracture manifolds

In this paper, a numerical model to predict flow‐induced shear failure along pre‐existing fractures is presented. The framework is based on a discrete fracture representation embedded in a continuum describing the damaged matrix. A finite volume method is used to compute both flow and mechanical equilibrium, whereas specifically tailored basis functions are used to account for the physics at discontinuities. The failure criterion is based on a maximum shear strength limit, which changes with varying compressive stress on the fracture manifold. The displacements along fracture manifolds are obtained such that force balance is achieved under conditions, where shear stress of the failing fracture segment is constrained to the maximum shear strength at the segment. Simultaneously, the fluid pressure is computed independently of the shear slip. A relaxation model approach is used to obtain the maximum shear limit on the fracture manifold, which leads to grid convergence.     

A fracture mapping and extended finite element scheme for coupled deformation and fluid flow in fractured porous media

This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd.