无限域横观各向同性岩体介质中折线裂纹边界元方法

岩体材料的各向异性导致其中的光滑裂纹在扩展后会变为折线裂纹。对于折线裂纹,其折线处裂纹面外法线方向不惟一,不能用连续单元离散。为此,在光滑裂纹问题方法的基础上,引入4种新的不连续单元来离散折线处裂纹面,建立了适用于折线裂纹问题的对偶边界元方法,该方法基于横观各向同性基本解。算例验证表明,该方法具有较好的精度。最后用该方法分析了横观各向同性岩体中的折线裂纹,得到了该类裂纹的应力强度因子。当裂纹面上作用法向均布力,横观各向同性岩体介质中的矩形光滑裂纹发生弯折时,折线两侧的裂纹面在张开时存在抑制效应,从而导致折线裂纹裂尖应力强度因子小于原光滑裂纹。同时还发现,随裂纹面的逐渐弯折,其裂纹面对各向同性面的倾角发生变化,因此,其裂尖断裂特性还受到岩体各向异性的影响。     

含多裂纹沥青路面的动力响应分析

考虑沥青材料的黏弹性,结合动力学和断裂力学理论,通过有限元数值模拟研究了含多裂纹沥青路面在移动交通荷载作用下的动力特性,考察了裂纹间距变化对动态应力强度因子以及裂纹尖端应力的影响。研究结果表明,在移动交通荷载作用下基层底面裂纹处于Ⅰ、Ⅱ型动态应力强度因子交替控制下;裂纹间距增大,对基底裂纹应力强度因子的影响逐渐减小;表面裂纹尖端在移动荷载距其较近时受压应力,且随裂纹间距增大,最大压应力先微弱减小后逐渐增大;在移动荷载经过表面裂缝时,裂缝尖端经历了方向相反的剪切作用,随着裂纹间距的增大,最大剪应力呈现先减小后增大的趋势。同时对加铺土工布这一阻裂措施利用数值模拟进行了验证。     

三维空间中周期裂纹阵受剪切作用时的断裂力学分析

本文研究了岩石断裂力学中感兴趣的三维空间中周期裂纹阵受剪切作用的问题。求出了以裂纹张开位移为未知量的奇异积分方程。对矩形裂纹面进行了数值计算,得到了应力强度因子沿裂纹边界的分布规律。     

裂纹面接触对中心裂纹圆盘应力强度因子影响分析

从理论和数值两个方面进行分析,发现受径向集中力和围压作用的中心裂纹圆盘(CCBD)试件裂纹面接触会对II型应力强度因子产生较大的影响。通过理论研究,分析CCBD受集中力和围压作用裂纹面接触时圆盘内部的应力场,采用断裂力学权函数理论,推导得出在集中力和围压共同作用下,考虑裂纹面闭合时应力强度因子的解析解。然后,使用ANSYS软件建立了相应的数值模型计算应力强度因子,并与理论解和相关文献进行对比验证,证明了理论公式的正确性。无论裂纹张合与否,所提出的解析公式都能计算出不同裂纹长度、加载角、围压和摩擦系数的应力强度因子。最后,利用公式分析摩擦系数对应力强度因子的影响,结果表明:随着摩擦系数的增大,I型裂纹的应力强度因子不变,II型裂纹的应力强度因子随之显著减小;当加载角较大时,裂纹面产生更为复杂的二次裂纹,故压剪断裂测试的推荐加载角范围为30°~50°。     

基于改进数值流形法的接触裂纹问题研究

接触裂纹问题在工程结构中较为常见。结合数值流形法在裂纹处理上的优势,分析了压剪荷载作用下的接触裂纹问题,模拟了压剪裂纹渐进扩展过程。为了减少由于裂纹尖端位置不同而产生的误差,对裂纹尖端附近一定范围内的每一个物理覆盖附加奇异覆盖函数项,并根据裂纹尖端位置和单元含奇异物理覆盖的数目进行分区积分。选取一个压剪破坏算例,分析了法向接触力对应力强度因子计算结果的影响,并模拟了其渐进破坏过程。计算结果表明,所提方法在压剪裂纹问题方面的可行性,与未细化和覆盖细化方法得到的结果相比,更能准确地描述裂纹扩展路径。法向接触力对II型应力强度因子的贡献为0,对I型应力强度因子的影响较大,相对误差随网格密度变化明显,且法向接触力对I型应力强度因子的影响要比直接施加内压时的影响大。     

基于3D-ILC球体45º三维双内裂纹复合断裂研究

断裂力学是各领域基础科学,材料内裂纹相互作用是断裂力学的重要研究领域,但目前球体内裂纹相互作用研究成果较少。基于三维激光疲劳内裂纹（3D-ILC）技术,在对球体试样表面无任何损伤的情况下,制作三维平行纯内裂纹,开展45o倾角平行双内裂纹单轴压缩试验,对双内裂纹的裂纹扩展过程、应力双折射规律、断口特征、破坏形态、相互作用等多方面开展研究,进行含双内裂纹I-II-III型圆球断裂数值模拟,分析KⅠ、KⅡ、KⅢ分布规律,与试验相符。结果表明：（1）球体双内裂纹在单轴压缩下之间发生屏蔽作用,主裂纹面穿过两预制内裂纹;（2）翼型裂纹光滑区呈现I-II型裂纹断口特征,翼型裂纹侧面撕裂区呈现“矛状”III型裂纹特征;（3）基于M积分计算的KI、KII、KIII分析与基于断裂特征的I、II、III型裂纹类型分析一致。3D-ILC技术可用于球体内裂纹相互作用断裂试验研究,球体内裂纹试验与理论分析为脆性材料内裂纹扩展及相互作用研究提供试验与理论参考。     

基于ANSYS的PDC钻头钎焊裂纹应力场和应力强度因子分析

金刚石复合片(PDC)钻头钎焊裂纹尖端两侧临近界面处应力高度集中,是导致其断裂失效的主要原因。基于ANSYS对PDC钻头钎焊裂纹进行数值模拟,针对焊缝裂纹形状、厚度、位置的不同,能对裂纹尖端应力场等值线图及应力强度因子进行数值计算。借助这一计算PDC焊缝裂纹应力强度因子的新方法,阐述了PDC断裂失效机理,用以指导钻头焊接工艺,稳定和提高PDC钻头焊接质量。      

闭合裂纹断裂的有效剪应力准则

将闭合裂纹表面的有效剪应力引入裂尖应力强度因子的计算,获得了无限大板和有限宽板含中心闭合裂纹在不同裂纹长度、倾角以及摩擦系数下裂纹尖端的应力强度因子值。引入等径向剪应力线 这一概念,建立了闭合裂纹断裂的有效剪应力准则：（1）岩板内闭合裂纹将沿着等 线上双剪应力的和最小的方向扩展;（2）裂纹尖端的应力强度因子KⅡ达到材料的临界值KⅡC,裂纹将开始扩展。该准则成功预测了闭合裂纹的临界起裂角 ,与各种经典复合型断裂准则计算Ⅱ型裂纹起裂角比较,结果较为接近。将其应用于闭合裂纹的断裂判定是安全的。     

冲击载荷作用下准脆性材料Ⅱ型裂纹扩展研究

冲击载荷作用下准脆性材料的动态断裂一直是关注的热点问题,Ⅱ型裂纹试样受冲击剪切时其裂纹扩展方向同材料力学性质和冲击速度等密切相关。应用岩石动态破裂过程动态分析系统软件,对单边平行双裂缝试样开展了冲击载荷作用下的裂纹动态扩展数值模拟,分别研究了不同材料力学性质、材料均质度、入射应力脉冲幅值和历时对II型裂纹动态扩展的影响。数值模拟结果表明,纯II型裂纹在动荷载作用下的扩展,不仅受到剪切损伤,而且还存在拉伸损伤;准脆性材料的非均匀性导致了主裂缝周围产生大量微裂纹的破坏,影响裂缝的分岔和内部的应力值;应力幅值和应力脉冲历时分别超过某一定值时,主裂缝将出现分叉现象,试样的破坏程度加剧,其研究结果对于深入揭示准脆性材料在动荷载作用下II型裂纹扩展的规律及准脆性材料的损伤断裂机制具有重要的参考价值。     

预裂爆破作用下岩石裂纹动态应力强度因子分析

通过有限元分析软件ANSYS数值模拟手段,分析了爆破荷载作用下,裂纹长度与类型、不同的装药量对裂纹尖端动态应力强度因子的影响以及预裂与光面爆破动态应力强度因子比较分析,计算得出：预裂爆破预裂缝的产生主要是从炮孔处产生的开口裂纹在冲击波以及爆生气体的作用下扩展形成的;随着药径与孔径比的增大,动态应力强度因子也逐渐增大,动态应力强度因子曲线形态不变;由于自由面的存在,光面爆破裂纹应力强度动态因子后续峰值较预裂爆破的大。     

横观各向同性地基上刚性矩形基础分析

对横观各向同性地基上刚性矩形板进行了求解。首先,利用表面受矩形均布荷载作用下的层状横观各向同性地基的位移解答,获得地基的柔度矩阵;然后,通过刚性矩形基础与层状横观各向同性地基的协调条件,建立刚性矩形基础与横观各向同性地基共同作用的方程,进而求得基础的地基反力。通过编制相应的程序,确定了合理的网格划分值;最后,进行算例分析,分析了地基横观各向同性性质、矩形刚性基础的长宽比以及地基分层性对地基反力的影响。分析表明：以上3种因素对地基反力有重要影响。     

Crack Initiation, Propagation and Coalescence in Rock-Like Material Containing Two Flaws: a Numerical Study Based on Bonded-Particle Model Approach

Cracking and coalescence behavior in a rectangular rock-like specimen containing two parallel (stepped and coplanar) pre-existing open flaws under uniaxial compression load has been numerically studied by a parallel bonded-particle model, which is a type of bonded-particle model. Crack initiation and propagation from two flaws replicate most of the phenomena observed in prior physical experiments, such as the type (tensile/shear) and the initiation stress of the first crack, as well as the coalescence pattern. Eight crack coalescence categories representing different crack types and trajectories are identified. New coalescence categories namely “New 1” and “New 2”, which are first observed in the present simulation, are incorporated into categories 3 and 4, and category 5 previously proposed by the MIT Rock Mechanics Research Group, respectively. The flaw inclination angle (β), the ligament length (L) (spacing between two flaws) and the bridging angle (α) (inclination of a line linking up the inner flaw tips, between two flaws) have different effects on the coalescence patterns, coalescence stresses (before, at or post the peak stress) as well as peak strength of specimens. Some insights on the coalescence processes, such as the initiation of cracks in the intact part of specimens at a distance away from the flaw tips, and coalescence due to the development and linkage of a number of steeply inclined to vertical macro-tensile cracks are revealed by the present numerical study.     

Elastic solutions for a transversely isotropic half-space subjected to buried asymmetric-loads

Elastic closed-form solutions for the displacements and stresses in a transversely isotropic half-space subjected to various buried loading types are presented. The loading types include finite line loads and asymmetric loads (such as uniform and linearly varying rectangular loads, or trapezoidal loads). The planes of transverse isotropy are assumed to be parallel to its horizontal surface. These solutions are directly obtained from integrating the point load solutions in a transversely isotropic half-space, which were derived using the principle of superposition, Fourier and Hankel transformation techniques. The solutions for the displacements and stresses in transversely isotropic half-spaces subjected to linearly variable loads on a rectangular region are never mentioned in literature. These exact solutions indicate that the displacements and stresses are influenced by several factors, such as the buried depth, the loading types, and the degree and type of rock anisotropy. Two illustrative examples, a vertical uniform and a vertical linearly varying rectangular load acting on the surface of transversely isotropic rock masses, are presented to show the effect of various parameters on the vertical surface displacement and vertical stress. The results indicate that the displacement and stress distributions accounted for rock anisotropy are quite different for those calculated from isotropic solutions. Copyright © 1999 John Wiley & Sons, Ltd.     

危岩应力强度因子的计算研究

危岩崩塌是山区主要地质灾害之一,其主控结构面的断裂失稳是导致危岩崩塌的根本原因。运用线弹性断裂力学的理论,引入应力集中系数法,建立了自重作用下的危岩应力强度因子的计算公式,可用于主控结构面与裂缝共线的危岩稳定性计算。同时还推导了基于ABAQUS 的J积分、位移外推法计算危岩应力强度因子的过程。通过实例表明,2种数值计算方法与解析结果十分接近,验证了ABAQUS计算危岩应力强度因子的合理性,最后模拟了危岩裂缝角度以及临空面宽度变化对危岩应力强度因子的影响。利用ABAQUS对危岩进行分析,并绘制应力强度因子在各因素影响下的趋势图,可以推测危岩失稳的临界条件。     

Surface load dynamic solution of saturated transversely isotropic multilayer half-space

This paper treats the dynamic response of a multilayered transversely isotropic fluid saturated poroelastic half-space under surface time-harmonic traction. The governing system of partial differential equations is uncoupled with the use of a set of physically meaningful and complete potential functions that decompose different body waves in a saturated poroelastic transversely isotropic medium. After expressing the equations in the Hankel-Fourier domain, a proper algebraic factorization is applied to generate reflection and transmission matrices for decomposed waves. All responses including displacements, stresses, and pore fluid pressure for both general patch load and point load are presented in the form of semi-infinite line integrals. The verification of the method is confirmed with the degeneration of the solutions presented here to the existing solutions for dried both homogeneous and multilayered elastic half-spaces as well as poroelastic half-space. Selected numerical results are depicted to investigate the effects of layering and pore pressure on responses of a transversely isotropic poroelastic medium. The load distribution effects are studied by comparison of the patch and point load responses. Also, resonance notion and effective parameters on this phenomenon such as layering system and anisotropy contrast are discussed. Significant influence of materials and layering configuration on number and amplitude of resonances depicted through the numerical evaluation.     

Elastic solutions for a transversely isotropic half-space subjected to a point load

We rederive and present the complete closed-form solutions of the displacements and stresses subjected to a point load in a transversely isotropic elastic half-space. The half-space is bounded by a horizontal surface, and the plane of transverse isotropy of the medium is parallel to the horizontal surface. The solutions are obtained by superposing the solutions of two infinite spaces, one acting a point load in its interior and the other being free loading. The Fourier and Hankel transforms in a cylindrical co-ordinate system are employed for deriving the analytical solutions. These solutions are identical with the Mindlin and Boussinesq solutions if the half-space is homogeneous, linear elastic, and isotropic. Also, the Lekhnitskii solution for a transversely isotropic half-space subjected to a vertical point load on its horizontal surface is one of these solutions. Furthermore, an illustrative example is given to show the effect of degree of rock anisotropy on the vertical surface displacement and vertical stress that are induced by a single vertical concentrated force acting on the surface. The results indicate that the displacement and stress accounted for rock anisotropy are quite different for the displacement and stress calculated from isotropic solutions. © 1998 John Wiley & Sons, Ltd.     

横观各向同性含液饱和多孔介质黏弹性瞬态反应

在考虑横观各向同性含液饱和多孔介质固体骨架和流体可压缩性以及固体骨架的黏弹性特征下,基于横观各向同性含液饱和多孔介质u-w形式的三维动力控制方程,以固相位移u、液相相对位移w为基本未知量,综合运用Laplace变换、双重Fourier变换等方法,在直角坐标系下通过引入中间变量,将六元2阶动力控制方程组化为两组各含4个未知变量的常微分方程组,给出了直角坐标系下横观各向同性含液饱和多孔介质三维黏弹性动力反应的积分形式一般解;作为理论推导的验证,通过引入初始条件和边界条件,对横观各向同性含液饱和多孔介质半空间黏弹性瞬态反应问题进行了求解。解答的退化验证表明,所推导的理论解是正确的。     

Elastic fields in two joined transversely isotropic media of infinite extent as a result of rectangular loading

This paper presents the closed‐form solutions for the elastic fields in two bonded rocks induced by rectangular loadings. Each of the two bonded rocks behaves as a transversely isotropic linear elastic solid of semi‐infinite extent. They are completely bonded together at a horizontal surface. The rectangular loadings are body forces along either vertical or horizontal directions and are uniformly applied on a rectangular area. The rectangular area is embedded in the two bonded rocks and is parallel to the horizontal interface. The classical integral transforms are used in the solution formulation, and the elastic solutions are expressed in the forms of elementary harmonic functions for the rectangular loadings. The stresses and displacements in the rocks induced by both the horizontal and vertical body forces are also presented. The numerical results illustrate the important effect of the anisotropic bimaterial properties on the stress and displacement fields. The solutions can be easily implemented for numerical calculations and applied to problems encountered in rock mechanics and engineering. Copyright © 2011 John Wiley & Sons, Ltd.