复合型加载条件下锈岩断裂韧度试验研究

岩石断裂韧度表征其抵抗裂纹起裂和扩展的能力,作为其力学性能的一个重要指标,在岩石力学理论研究与岩体工程应用中有着不可替代的作用。由于岩石结构的复杂性,其破坏形式大多呈现为复合型断裂破坏,研究复合型加载条件下岩石断裂韧度具有重要的意义。为了研究锈岩的断裂力学性能,以便其在建筑工业中得到更为广泛的应用,通过18个中心裂纹圆盘（CSTBD）试件径向受压试验,进行了复合型加载条件下锈岩断裂韧度的试验研究。测得了锈岩在纯I型,纯II型以及复合型加载时的断裂韧度,并将试验结果与基于广义最大周向应力（GMTS）准则的理论值进行了对比分析,结果表明：锈岩的纯I型断裂韧度为1.01 MPa?m0.5,而纯II型断裂韧度为1.51 MPa?m0.5,是纯I型断裂韧度的1.49倍,这与基于GMTS准则的理论值1.34非常接近,而比基于最大周向应力（MTS）准则的理论值0.87大很多。裂纹尖端附近的T应力及断裂过程区裂纹尖端的临界距离rc对岩石类材料的开裂路径以及复合型断裂韧度都有较大的影响。考虑了T应力的广义最大周向应力（GMTS）准则能很好的预测试验结果。     

T应力对岩石裂纹扩展路径及起裂强度的影响研究

经典断裂判据在分析岩石起裂角度及强度时,往往只选择裂尖应力场展开式中的r1/2奇异应力项,而将高阶的O(r1/2)项以及非奇异应力项（T应力）忽略,造成理论预测并不能完全吻合实际的试验结果。在充分考虑非奇异应力项对裂纹扩展影响作用的基础上,利用最大周向应力判据重新研究了岩石类材料脆性破坏的I、II及I-II复合型裂纹扩展。研究结果表明：（1）纯I型裂纹状态下,T应力为负时（压缩）,裂纹的扩展是稳定的;而当T 应力为正时（拉伸）,只有当 ,裂纹才沿着裂纹方向扩展（经典断裂判据所预测的方向）,而其他情况下裂纹的扩展将发生偏转;（2）纯II型裂纹时,不仅沿着裂纹方向的T应力影响起裂角度及强度,垂直裂纹面上的N应力也同样具有重要的影响作用;（3）对I-II型裂纹而言,修正后的最大周向应力准则与试验结果吻合良好,且正的T应力增大开裂角,而负的T应力降低开裂角。通过控制T应力的大小,可以对裂纹扩展方向加以控制,从而使得裂纹扩展偏离最危险的方向,最大地限度阻止或延缓结构整体断裂的发生。     

拉剪应力状态下岩体裂隙扩展的断裂力学机制及物理模型试验

开挖卸荷岩体裂隙面通常处于拉剪应力状态。在裂隙应力和变形状态分析的基础上,采用线弹性断裂力学理论和物理模型试验研究了拉剪应力状态下裂隙扩展的力学机制。岩体裂隙在拉剪应力状态下沿裂隙面间的滑动抗剪摩擦力消失,裂隙起裂沿Ⅰ型张拉裂隙断裂韧度KⅠ最小的方向起裂,并最终发展与卸荷拉应力方向垂直;拉剪应力状态下岩体的总位移方向平行于拉应力方向,通过合理的位移假设,基于能量及线弹性断裂力学理论,求解了拉剪应力状态下分支裂隙扩展过程中尖端的动态应力强度因子和扩展长度判据;通过拉剪应力状态下单裂隙扩展物理模型试验验证了理论推导的正确性。     

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

Determination of Fracture Toughness of Anisotropic Rocks by Boundary Element Method

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     

Numerical Observation of Three-Dimensional Wing Cracking of Cracked Chevron Notched Brazilian Disc Rock Specimen Subjected to Mixed Mode Loading

The cracked chevron notched Brazilian disc (CCNBD) specimen has been suggested by International Society for Rock Mechanics for measuring mode I fracture toughness of rocks. Subsequently, this specimen geometry has been widely extended to conduct mixed mode fracture tests on rocks as well. A straight through crack front during the fracturing process upon the root of the chevron notch is assumed in the testing principle, but has never been thoroughly evaluated before. In this study, for the first time, the progressive rock fracture mechanism of the CCNBD rock specimen under mixed mode loading is numerically simulated. Specimens under representative mixed mode loading angles are modelled; and the assumption of the straight through crack front growth is critically assessed. The results show that not only the notch tip but also the saw-cut chevron notch cracks during the experiments, yielding a prominent twisted front, far from being straight. The crack front never grows up to the root of the notch ligament and the straight through crack front assumption is never satisfied in the realistic rock fracture progress of this chevron notched specimen subjected to mixed mode loads. In contrast, the fracture progress features typical three-dimensional wing cracking towards the loading ends. The numerically observed progressive fracture mechanism reveals that the measuring principle of mixed mode fracture tests employing CCNBD specimens is significantly violated and the measures of both modes I and II fracture toughness are uncertain.     

Mode I crack propagation in concrete under fatigue: microscopic observations and modelling

In this study, three‐point bending tests were carried out on notched beams to investigate mode I crack propagation in plain concrete under fatigue. The first part of the study focused on microscopic observations of the crack growth features. Microscopic observations were made using the replica method associated with scanning electron microscopy (SEM). Observations of fatigue crack growth both on the surface and inside the specimens are presented as a comparison between the observed crack lengths and those estimated by the compliance calibration method. In the second part, a finite element model of mode I crack propagation under fatigue is presented. According to the cohesive crack concept, a cohesive force distribution on the crack at various loading stages is assumed, according to both the stress‐crack opening relation worked out by Hordijk (1991; Thesis, Technische Universiteit) and a new proposed relation with hysteresis loop. Finite element computation is used to evaluate the crack extension in the bending beams. Numerical predictions are discussed in comparison with experimental results. Copyright © 2002 John Wiley & Sons, Ltd.     

Modeling of dynamic cohesive fracture propagation in porous saturated media

In this paper, a mathematical model is presented for the analysis of dynamic fracture propagation in the saturated porous media. The solid behavior incorporates a discrete cohesive fracture model, coupled with the flow in porous media through the fracture network. The double‐nodded zero‐thickness cohesive interface element is employed for the mixed mode fracture behavior in tension and contact behavior in compression. The crack is automatically detected and propagated perpendicular to the maximum effective stress. The spatial discretization is continuously updated during the crack propagation. Numerical examples from the hydraulic fracturing test and the concrete gravity dam show the capability of the model to simulate dynamic fracture propagation. The comparison is performed between the quasi‐static and fully dynamic solutions, and the performance of two analyses is investigated on the values of crack length and crack mouth opening. Copyright © 2010 John Wiley & Sons, Ltd.     

Discrete element model for crack propagation in brittle materials

We propose a discrete element model for brittle rupture. The material consists of a bidimensional set of closed‐packed particles in contact. We explore the isotropic elastic behavior of this regular structure to derive a rupture criterion compatible to continuum mechanics. We introduce a classical criterion of mixed mode crack propagation based on the value of the stress intensity factors, obtained by the analysis of two adjacent contacts near a crack tip. Hence, the toughness becomes a direct parameter of the model, without any calibration procedure. We verify the consistency of the formulation as well as its convergence by comparison with theoretical solutions of tensile cracks, a pre‐cracked beam, and an inclined crack under biaxial stress. Copyright © 2015 John Wiley & Sons, Ltd.     

考虑裂隙水压力的岩体压剪裂纹扩展规律研究

在库水位大幅度升、降变化时,容易导致岩体内增量裂隙压力的集中,使断裂面上有效应力降低,裂纹面尖端的应力强度因子增加。当达到临界强度因子时,可能使岩体内裂纹、裂隙贯通、扩展,形成连续的复式破坏面,从而使边坡稳定性降低,造成边坡的失稳。基于此,从断裂力学角度分析了裂隙水压力对裂纹强度因子的影响,对考虑裂隙水压力作用的Ⅰ、Ⅱ型复合裂纹扩展规律进行了研究,结果表明：Ⅰ、Ⅱ型复合裂纹的裂纹扩展角的变化,不仅与裂纹的闭合程度、斜裂纹倾角、双向应力大小有关,还与裂隙水压力的大小、裂纹面的摩擦系数有关;并且在相同情况下,未闭合裂纹的扩展角要大于闭合裂纹的扩展角;对于闭合裂纹,裂纹面摩擦系数越小,扩展角越大;最后,推导了基于摩尔-库仑准则考虑裂隙水压力的岩体断裂韧度KIc、KIIc和压剪状态下Ⅰ、Ⅱ型复合断裂判据。研究成果为分析水岩作用下裂隙岩体的失稳破坏提供了重要的参考。     

砂岩的I/III复合型断裂试验研究

为了了解砂岩的I/III复合型断裂力学性能,采用含有裂纹面垂直于底面但倾斜于正背面的三点弯曲（TPB）试件实现面外断裂。使用Abaqus通过有限元法研究了不同裂纹面倾斜角下沿着试件厚度的I型、II型与III型无量纲应力强度因子。通过讨论可知,在靠近试件厚度中点处I型与III型无量纲应力强度因子取得较大值,而II型无量纲应力强度因子取得接近于0的可以忽略的极小值,且靠近自由面处的应力强度因子不作考虑。故该种试件可用于测量材料的I/III复合型断裂性能。数值计算结果与常规TPB试件的解析公式及倾斜裂纹面TPB试件的近似计算公式的结果进行对比,证明了数值模拟结果的准确性并讨论了近似公式存在的问题。使用7组共28个该种三点弯曲试件,研究砂岩的I/III复合型断裂特性,结果表明：砂岩的I型与III型临界应力强度因子随着裂纹面倾斜角的增大先增大后减小,但其变化趋势不同;砂岩的I/III复合型等效断裂韧度随着裂纹面倾斜角的增大先增大后减小,在倾斜角接近20°时取得最大值。裂纹扩展时裂纹面将发生扭转,倾斜角θ越大则扭转角度越大。     

单轴压缩条件下含三叉裂隙类岩石试样力学特性的细观研究

三叉裂隙是自然界普遍存在的一种岩体缺陷形式,其对岩体的力学特性有重要影响。对含预制三叉裂隙的水泥砂浆试样进行室内单轴压缩试验,配合使用摄像机拍摄裂纹的起裂、扩展、贯通过程,通过数字图像技术处理获取试样的应变场云图,并结合PFC2D程序研究不同?、? 条件下试样的强度特征、裂纹模式和裂纹演化扩展规律。研究表明：三叉裂隙对试样单轴抗压强度有明显的削弱作用。当? 恒定为120°时,试样在? = 30°时达到最大抗压强度;当? 恒定为90°时,随?增大,试样抗压强度呈先减小后增大的趋势,且当? = 45°时达到最大抗压强度。试样产生的裂纹可分为3类,分别是张拉型裂纹(Ⅰ型裂纹)、剪切型裂纹(Ⅱ型裂纹)、混合型裂纹(Ⅲ型裂纹)。这3类裂纹通常从裂隙尖端开始产生,并且Ⅰ型裂纹沿加载方向扩展,通常未扩展至试样边界;Ⅱ型和Ⅲ型裂纹通常与加载方向呈一定角度扩展至试样边界。通过对裂纹的几何形态和组成宏观裂纹的微裂纹成分的分析,得知导致含三叉裂隙试样在单轴压缩条件下失效的是张拉破坏。数字图像技术得到的应变云图表明,当载荷达到一定阶段,裂隙尖端出现应力集中,微破裂开始发育并聚集成微破裂区,微破裂区扩大产生宏观裂纹。通过对主应变和剪应变云图分析,发现导致试样失效的是张拉破坏,剪应变在裂纹扩展过程中的影响较小。     

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

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

Mixed mode crack propagation in low brittle rock-like materials

Mixed mode fracture is quite common in rock structures. Numerous investigators have used the Brazilian disk specimens with a central crack for investigating modes I, II, and mixed fracture toughness in brittle materials. In this study, analytical, experimental, and numerical investigations were planned and performed on Central Straight Through Crack Brazilian Disk (CSCBD) specimens. Ranking of geometrical parameters effective on the value of stress intensity factors (SIFs) of CSCBD specimens were obtained using stochastic analysis. Furthermore, experimental tests were undertaken in order to evaluate the crack propagation in rock-like material of low brittleness. Finally, numerical modeling was performed to assess the effect of crack length on the failure mode of CSCBD specimens. Analytical analyses revealed that the inclination angle of the crack with respect to the diametrical load has the most important impact on the SIFs among the geometrical parameters of CSCBD specimen. Performed experimental and numerical analyses also confirmed the effect of inclination angle and crack length and their impact on the mode of failure of the tested specimen.     

单轴压缩非贯通节理岩石尖端塑性区及弹塑性断裂模型研究

基于Drucker-Prager（下简称D-P）准则,建立压缩载荷作用下的非贯通节理岩石的弹塑性断裂模型。针对节理岩石小范围屈服翼裂纹尖端塑性区,推导了D-P屈服准则的纯I、纯II及I、II复合型3种翼裂纹无量纲塑性区径长函数,并与Mises准则的塑性区进行对比;结果表明,D-P准则的I型和复合型塑性区较Mises屈服准则的塑性区大,且其II型及I、II复合型塑性区在翼裂纹上下表面不连续。进一步,引入断裂软化因子以表征节理岩石裂隙断裂扩展后的断裂软化规律,考虑非贯通节理岩石复合型断裂软化,是由于节理尖端翼裂纹应变能密度超过最小应变能密度导致其成核扩展引起的,提出用应变能密度的指数函数形式表征断裂软化变量的演化;塑性屈服函数采用Borja等的应力张量3个不变量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;建立节理岩石的弹塑性断裂本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性断裂模型的程序。以单轴压缩下非贯通节理岩石为例,分析岩石断裂韧度、节理摩擦系数和节理倾角等参数的影响,结果表明,所提出的弹塑性断裂模型与数值和试验结果比较吻合。     

条形药包爆破预制贯通裂纹动态断裂过程研究

采用动态焦散线试验和ABAQUS数值分析方法,对条形药包爆破载荷作用下不同角度预制贯通裂纹的扩展行为进行了研究。研究结果表明：在爆炸应力波的作用下,柱部区域和端部区域0°预制贯通裂纹远端翼裂纹背离炮孔方向扩展,而近端翼裂纹朝向炮孔扩展。预制贯通裂纹远端的应力集中程度较近端高,并且远端翼裂纹的止裂韧度更低,翼裂纹更易扩展。在爆炸应力波的作用下,柱部区域90°预制贯通裂纹由张开逐渐转为闭合,爆炸应力波在已闭合的预制贯通裂纹面发生透射,并在预制贯通裂纹尖端产生压剪应力集中,形成以II型断裂为主的I-II复合型裂纹,并近似垂直于预制贯通裂纹面扩展,随后,在自由面反射应力波的作用下,反翼裂纹沿预制贯通裂纹面起裂扩展;炮孔端部区域90°预制贯通裂纹处翼裂纹的起裂是由于爆炸应力波在预制贯通裂纹处产生反射拉伸波的结果,促使预制贯通裂纹端部产生拉应力集中,形成近似I型裂纹,随后,翼裂纹逐渐转向爆炸应力波传播方向扩展。     

基于扩展有限元的压剪复合断裂数值分析

针对岩石类材料压剪破坏,国内外学者做了大量真实岩石和模型材料的试验,并进行相应的数值模拟工作。但由于试验结果往往具有很强的离散性,且对影响试验结果的各因素缺乏系统地理论分析,加之传统常规有限元在处理裂纹这类强不连续问题时存在很大的困难,如在裂尖附近的高应力区需要剖分令人难以接受的密度的网格,同时在模拟裂纹生长时还需要对网格进行重新剖分,操作繁琐复杂,效率极低,因此,针对压剪复合断裂翼裂纹扩展规律的研究始终没有取得令人满意的成果。对原生裂纹的几何形态,主裂纹面与其间填充物的相互作用对翼裂纹起裂角的影响做了系统的理论分析,并利用扩展有限元对翼型裂纹扩展全过程进行数值模拟,与理论分析结果对比,吻合良好。相关成果为以后针对压剪复合断裂的研究提供了理论分析框架,同时对今后相关试验的开展具有理论指导作用     

Prediction of fracture trajectory in anisotropic rocks using modified maximum tangential stress criterion

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.     

圆孔缺陷与I型运动裂纹相互作用的试验研究

为了研究运动裂纹与圆形孔缺陷的相互作用机制,采用数字激光动态焦散线方法进行了含圆形孔缺陷的冲击试验。结果表明：在冲击载荷下,I型运动裂纹在与圆形孔贯通前,断裂面光滑,扩展路径平直;运动裂纹从圆形孔上端起裂后,断裂面凹凸不平,扩展路径也更为弯曲。当运动裂纹朝向邻近的圆形孔扩展时,圆形孔对运动裂纹的扩展速度和动态应力强度因子有抑制作用,且随着圆形孔直径的增大,这种抑制作用不断增强。当运动裂纹与圆形孔缺陷汇聚后,圆形孔阻碍了裂纹的继续扩展,裂纹尖端被钝化,钝裂纹的起裂韧度较尖裂纹提高9.58%～13.87%,裂纹再次起裂时的扩展速度和动态应力强度因子存在明显的跳跃,表明钝裂纹更难起裂,需要消耗的能量更多。     

Impact of water and nitrogen fracturing fluids on fracturing initiation pressure and flow pattern in anisotropic shale reservoirs

A numerical model is proposed to investigate the impact of water and nitrogen fracturing fluids on the fracturing initiation pressure and the flow pattern in anisotropic shale reservoirs. This model considers the anisotropy of shale deformation and permeability, the compressibility of fracturing fluid, and the fluid moving front. A crack initiation criterion is established with the stress intensity factor of mode I crack. Both crack initiation pressure and seepage area are verified and analyzed. These results show that shale deformation and permeability anisotropy, fracturing fluid compressibility, viscosity, and pressurization rate have significant impacts on fracturing initiation pressure and seepage area.