A 3-D ellipsoidal flaw model for brittle fracture in compression

Engineering materials are rarely free of flaws. Mode I cracking from pre-existing flaws is the major cause of the brittle fracture in compression of materials such as concrete and rock. A 3-D ellipsoidal flaw model is used to show the significant influence of flaw geometry on crack initiation in uniform uniaxial, biaxial and triaxial compression. The model shows that the governing criterion for crack initiation may change from energy to stress with increasing crack size, and that for voids of similar size a spherical void is the most critical shape for crack initiation. The model thus provides a basis for a better understanding of both the phenomenon and the mechanism of brittle fracture in compression.     

An interacting micro-crack damage model for failure of brittle materials under compression

A model is developed for brittle failure under compressive loading with an explicit accounting of micro-crack interactions. The model incorporates a pre-existing flaw distribution in the material. The macroscopic inelastic deformation is assumed to be due to the nucleation and growth of tensile “wing” micro-cracks associated with frictional sliding on these flaws. Interactions among the cracks are modeled by means of a crack-matrix-effective-medium approach in which each crack experiences a stress field different from that acting on isolated cracks. This yields an effective stress intensity factor at the crack tips which is utilized in the formulation of the crack growth dynamics. Load-induced damage in the material is defined in terms of a scalar crack density parameter, the evolution of which is a function of the existing flaw distribution and the crack growth dynamics. This methodology is applied for the case of uniaxial compression under constant strain rate loading. The model provides a natural prediction of a peak stress (defined as the compressive strength of the material) and also of a transition strain rate, beyond which the compressive strength increases dramatically with the imposed strain rate. The influences of the crack growth dynamics, the initial flaw distribution, and the imposed strain rate on the constitutive response and the damage evolution are studied. It is shown that different characteristics of the flaw distribution are dominant at different imposed strain rates: at low rates the spread of the distribution is critical, while at high strain rates the total flaw density is critical.     

An experimental and numerical study of fracture coalescence in pre-cracked specimens under uniaxial compression

This study presents crack initiation, propagation and coalescence at or near pre-existing open cracks or flaws in a specimen under uniaxial compression. The flaw geometry in the specimen was a combination of a horizontal flaw and an inclined flaw underneath. This flaw geometry is different from those reported in the previous studies, where a pair of parallel flaws was used. Three materials were used, PMMA (Poly Methyl MethAcrylate), Diastone (types of molded gypsum), and Hwangdeung granite. Crack initiation and propagation showed similar and different patterns depending on the material. In PMMA, tensile cracks initiated at the flaw tips and propagated to the tip of the other flaw in the bridge area. The cracks then coalesced at a point of the inclined flaw, which is affected by the flaw inclination angle. For Diastone and Hwangdeung granite, tensile cracks were observed followed by the initiation of shear cracks. Coalescence occurred mainly through the tensile cracks or tensile and shear cracks. Crack coalescence was classified according to the crack coalescence types of parallel flaws for overlapping flaw geometry in the past works. In addition, crack initiation and coalescence stresses in the double-flawed specimens were analyzed and compared with those in the single-flawed specimen. Numerical simulations using PFC^2D (Particle Flow Code in two dimensions) based on the DEM (Discrete Element Method) were carried out and showed a good agreement with the experimental results in the coalescence characteristics in Hwangdeung granite. These experimental and numerical results are expected to improve the understanding of the characteristics of cracking and crack coalescence and can be used to analyze the stability of rock and rock structures, such as the excavated underground openings or slopes, tunneling construction, where pre-existing cracks or fractures play a crucial role in the overall integrity of such structures.     

考虑微裂纹相互作用的的岩石微观力学弹塑性损伤模型研究

本文建立基于微裂纹扩展的岩石弹塑性损伤微观力学模型。用自洽方法考虑裂隙间相互影响,压缩载荷下微裂纹尖端翼裂纹稳定扩展表征岩石的微观损伤,基于应变能密度准则用Newton迭代法求复合型断裂的翼裂纹扩展长度,并采用微裂隙统计的二参数Weibull函数模型反映绝对体积应变对微裂纹分布数目影响,进而用翼裂纹扩展所表征的应力释放体积和微裂纹数目来表示含有微裂隙的岩石损伤演化变量;宏观塑性屈服函数采用Voyiadjis等的等效塑性应变的硬化函数,反映塑性内变量对硬化函数的影响;建立岩石的弹塑性损伤本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性损伤模型的程序。从围压和微裂隙长度等因素分析弹塑性损伤模型的岩石的损伤和宏观塑性特性。     

含微裂纹材料的损伤理论

本文从含微裂纹材料的变形能出发引出了裂纹的方位张量。在考虑裂纹受压闭合与滑动摩擦的基础上,给出了损伤张量、损伤应变及有效弹性常数。文中给出了损伤机构离散化的方法,并对方位密度给出了演化方程。最后给出一个单向拉压的应力应变关系例子,并揭示了裂纹扩展时的应力突跌现象。     

Influence of flaw inclination angle and loading condition on crack initiation and propagation

With reference to the experimental observation of crack initiation and propagation from pre-existing flaws in rock specimens under compression, the influences of pre-existing flaw inclination angle on the cracking processes were analyzed by means of finite element method (FEM) and non-linear dynamics method. FEM analysis on the stress field distribution induced by the presence of a pre-existing flaw provided better understanding for the influence of flaw inclination angle on the initiation position and initiation angle of the potential cracks. Numerical analysis based on the non-linear dynamics method was performed to simulate the cracking processes. The resultant crack types, crack initiation sequences and the overall crack pattern were different under different loading conditions. Under a relatively low loading rate or a small magnitude of maximum loading pressure, tensile cracks would tend to initiate prior to shear cracks. In contrast, under a relatively high loading rate and a large magnitude of maximum loading pressure, shear cracks would tend to initiate prior to tensile cracks instead.     

压缩载荷下孔隙结构变化的CT实验研究

为了研究孔隙对岩石力学性能的影响,利用自制的孔隙物理模型,通过 单轴压缩和CT扫描实验研究了受载条件下孔隙率对岩石孔隙结构的演化及其对外部物理力 学性能的影响,得到了不同加载阶段和不同CT观察尺度下孔隙模型的裂纹扩展规律以及孔隙 和固体介质的损伤变化情况. 实验结果表明：孔隙模型在受载条件下裂纹主要发生在峰值荷 载之后,主裂纹大都集中在孔隙密集的地方且伴随许多细小裂纹的产生;峰值载荷前出现了 少数微裂纹,微裂纹的产生与演化主要发生在孔隙周边.     

双轴压缩载荷下裂纹动态传播的实验研究

本文研究了双轴压缩载荷下裂纹的起裂和非共面传播问题，用动态光弹性方法成功地记录下裂纹的动态扩展过程。获得了双轴压缩载荷下扩展裂纹的传播速度，临界开裂载荷以及开裂角第一系重要断裂参数。研究表明，在双轴压缩载荷下不闭合裂纹前缘应力场是一个混合型的奇异场，裂纹的扩展过程为稳态过程并且表现出显著的阶段性特点。     

EXPERIMENTAL STUDIES ON DYNAMIC PROPAGATION OF CRACKS UNDER BIAXIAL COMPRESSION

This paper deals with the initiation and non-coplaner propagation ofcracks under biaxial compression.The running process of the propagating cracks issuccessfully recorded with the method of dynamic photoclasticity.A series of importantfracture parameters such as cracks' propagating velocity,critical opening loading andinitial angle of cracks' growth under hiaxial compression is calculated.Results showthat the stress field in front of the non-closed crack tip under biaxial compression is amixed-mode singular field,the dynamic propagation of the new tensile cracks is juststable and shows distinct stages,but the unstable crack appears under high loads andresults in the rupture of the matrix.     

An application of a damage constitutive model to concrete at high temperature and prediction of spalling

A characteristic feature of concrete under uniaxial compression is the development of cracks parallel to the loading direction. A damage constitutive model proposed by Ortiz [Ortiz, M., 1985. A constitutive theory for the inelastic behaviour of concrete. Mech. Mater. 4, 67–93] can predict the transverse tensile stress responsible for these cracks by considering the interaction between the aggregate and the mortar and the development of damage in the latter. When concrete is exposed to high temperature, as is the case during fire, the failure mode is thermal spalling. In order to improve the prediction of the stresses involved in this failure Ortiz’s model is extended to account for the effects of high temperature. Published experimental results for uniaxial and biaxial compression at high temperatures are used to calibrate the temperature dependence of some of the material properties. The transient creep strain is accounted for by modifying the constrained thermal strain. The stress analysis is coupled with hygro-thermal analysis of heat, mass transfer and pore pressure build-up. The effect of pore pressure on the damage evolution is modeled by applying a body force in the stress analysis module proportional to the pressure gradient. A numerical example of concrete under fire is solved and the computed results are discussed. Spalling is predicted when the damage variable reaches its maximum value of unity. The predicted depth and time of spalling for a range of variation of permeability and initial liquid water content are presented. They are in good agreement with published experimental results.     

单轴荷载下饱水岩石静态和动态抗压强度的细观力学分析

由于单轴荷载下饱水岩石的动态力学特性与静态力学特性存在很大差异,从宏观上进行力学分析 存在局限性。根据岩石受压全应力应变曲线的细观机制,分析了静态及动态单轴荷载条件下孔隙水影响饱水 岩石裂纹扩展的情况。在静态单轴压缩条件下,初始裂隙受压使自由水产生孔隙水压力,自由水对翼裂纹有 向外挤压的应力,促进裂纹扩展。在动态单轴压缩条件下,自由水会产生粘结力,抑制裂纹扩展。根据翼裂纹 受压扩展原理,推导出饱水单轴条件下动态抗压强度、静态抗压强度的计算公式,在相同断裂韧度下,饱水岩 石静态抗压强度风干岩石静态抗压强度饱水岩石动态抗压强度。对自然风干和饱水砂岩进行单轴静态、 动态压缩实验,结果与理论模型的结果相符。     

Numerical study of failure behaviour of pre-cracked rock specimens under conventional triaxial compression

Macroscopic pre-existing flaws play an important role in evaluating the strength and the failure modes of a heterogeneous rock mass. Crack initiation, propagation and coalescence from macroscopic pre-existing flaws are considered in a 3-D numerical model (RFPA3D) to investigate their effects on the underlying failure modes of rock. A feature of the code RFPA3D is that it can numerically simulate the evolution of cracks in three-dimensional space, as well as the heterogeneity of the rock mass. Three types of flaw geometries were evaluated numerically against experimental results: Type A for intact specimen, and Types B and C for flawed cylindrical specimens with different macroscopic pre-existing flaws, respectively. The effect of confining pressure on the fracture evolution was also considered. Numerical results showed that both the ligament angle and the flaw angle of two pre-existing cracks can affect the uniaxial compressive strength of the specimen and the mechanism of fracture evolution. In addition, both the uniaxial compressive strength and the accumulated acoustic emission increase with increasing heterogeneity.     

单轴条件下砂岩三维破裂过程的CT观测

岩石破裂过程研究一直是岩石力学专家关注的重要问题。本文采用特制的三轴压力仪与医用西门子SOMATOM -plusCT扫描仪结合 ,对砂岩进行了室内单轴压缩试验。通过对砂岩的CT图像和密度损伤增量与应力关系曲线分析 ,结果显示 ,砂岩的破裂演化过程可分为初始损伤的压密、裂纹出现—扩展、裂纹归并—分岔、裂纹重分岔—扩展以及裂纹惯通—宏观破坏等五个阶段。在初始损伤的压密阶段 ,砂岩的密度损伤增量为正值 ,速率也为正 ;在裂纹出现—扩展阶段 ,砂岩出现局部密度损伤增量减小 ,并随应力增加而由正值转为负值 ,速率也由正变负 ;在裂纹归并—分岔阶段 ,砂岩的密度损伤增量全为负值 ,速率也变快 ;在裂纹重分岔—扩展阶段 ,砂岩的密度损伤增量继续变负 ,但速率变慢 ;在裂纹惯通—宏观破坏阶段 ,砂岩的密度损伤增量继续变负 ,速率变的更快。     

Numerical study on coalescence of two pre-existing coplanar flaws in rock

Crack propagation and coalescence processes are the fundamental mechanisms leading to progressive failure processes in rock masses, in which parallel non-persistent rock joints are commonly involved. The coalescence behavior of the latter, which are represented as pre-existing coplanar flaws (cracks), is numerically investigated in the present study. By using AUTODYN as the numerical tool, the present study systematically simulates the coalescence of two pre-existing coplanar flaws in rock under compression. The cumulative damage failure criterion is adopted in the numerical models to simulate the cumulative damage process in the crack initiation and propagation. The crack types (shear or tensile) are identified by analyzing the mechanics information associated with the crack initiation and propagation processes. The simulation results, which are generally in a good accordance with physical experimental results, indicate that the ligament length and the flaw inclination angle have a great influence on the coalescence pattern. The coalescence pattern is relatively simple for the flaw arrangements with a short ligament length, which becomes more complicated for those with a long ligament length. The coalescence trajectory is composed of shear cracks only when the flaw inclination angle is small (such as β ? 30°). When the pre-existing flaws are steep (such as β ? 75°), the coalescence trajectory is composed of tensile cracks as well as shear cracks. When the inclination angle is close to the failure angle of the corresponding intact rock material, and the ligament length is not long (such as L ? 2a), the direct shear coalescence is the more favorable coalescence pattern. In the special case that the two pre-existing flaws are vertical, the model will have a direct tensile coalescence pattern when the ligament length is short (L ? a), while the coalescence between the two inner flaw tips is not easy to achieve if the ligament length is long (L ? 2a).     

压剪条件下预埋椭圆裂纹三维扩展实验研究

试制出两种不同性质的脆性类岩石材料：透明的低温树脂材料和非透明的砂浆类材料。实验研究了这两种材料内置椭圆形裂隙的三维始裂、扩展过程及其破裂过程中声发射信号特征,并采用三维断裂理论确定了内部裂纹的三维始裂状态。实验结果表明,在单轴压缩荷载作用下,两种材料中内置裂隙初期破裂面形态基本一致,均为一对反对称的包裹状翼裂纹,但初始起裂位置有所不同。透明树脂试样中,初始的新生裂纹几乎在预制裂纹长轴端部附近同时起裂,但起裂点并不在裂隙长轴端部,而在砂浆类材料内的起裂位置则靠近短轴端部,亦不在短轴端部。相同条件下,两种材料的最终破裂状态明显不同,透明树脂试样是以宏观劈拉破裂为主,而砂浆类试样则以宏观剪切破裂为主,这可能与材料性质和三维裂隙的起裂位置不同有关。三维断裂分析结果与实验结果一致。二维简化模型仅为三维模型的主法平面上的特征解,无法反映出III型断裂模式在三维裂隙的起裂和扩展中的作用。三维断裂机制分析应加强对III型断裂模式及其相关的复合模式的深入分析。     

单轴压缩下页岩类脆性材料细观损伤演化的离散元分析

了解广泛存在的类似页岩的脆性材料各向异性对工程安全具有重要意义。本研究将页岩视为粘结颗粒材料,基于离散单元方法研究了横观各向同性脆性页岩的损伤演化。再现了不同层理角的页岩试样的破坏模式,并对比了实验和数值模拟的抗压强度和弹性模量。引入微裂纹的概念,通过定义裂纹密度函数,系统地研究了单轴压缩条件下,页岩层理角对细观结构的影响。此外基于平均配位数建立了配位数变化与细观损伤的联系,并根据配位数的变化与裂纹数量将加载过程分为三个阶段,分析了不同阶段配位数与裂纹数量的对应变化关系。研究表明,页岩的裂纹密度随着层理角的增加而增加,而试样的平均配位数在加载过程中先上升后剧烈下降,颗粒集合体在单轴压缩条件下的应力应变及裂纹数量曲线与平均配位数曲线有良好的一致性。该研究揭示了横观各向同性脆性岩石的破坏过程和内在机理,将为页岩类脆性材料的工程应用提供理论指导。     

近场动力学研究进展与岩石破裂过程模拟

近场动力学（Peridynamics,PD）作为一种新兴的非局部性理论,在非连续处不需要任何处理,能够很好表述模型从连续到非连续的过程。首先,在PD基本理论简介的基础上,系统回顾了PD的国内外研究现状。其次,采用键型PD理论对非均匀性的圆孔岩板单轴拉伸破裂过程进行了二维数值模拟,采用态型PD理论对单轴、常规三轴以及真三轴等不同压缩条件下的岩石破裂过程进行了三维数值模拟,并以加拿大Mine-by隧洞为例对现场岩体破裂过程进行了模拟,结果表明PD在岩石破裂过程模拟上具有较强适用性。最后,指出当前PD在岩石破裂过程模拟中存在的主要问题和未来值得开展的若干研究课题。     

基于3D⁃DIC技术的约束岩石裂缝扩展研究

采用 RMT-301B岩石与混凝土力学实验系统对施加预应力值为峰值强度 20 %、40 %、60 %的泥质白云 岩进行单轴压缩实验,并借助三维数字图像相关技术(3D-DIC)测试系统观测受约束岩石试样的变形和破坏形 态,研究泥质白云岩在约束条件下的位移场和 Tresca应变场的演化过程．结果表明：(1) 岩石的破坏过程是由 内向外扩展,当应力积累超出岩石承受范围后,岩石出现裂纹,随后裂纹尖端继续扩展,应力也发生了重分 布,岩石内部多组裂纹开始联合扩展、贯穿,直到试件发生完全破坏．(2) 含约束岩石在单轴压缩下一共有轴 向劈裂破坏、沿片理破坏、Y形破坏、共轭剪切破坏四种破坏形态,且以轴向劈裂破坏为主     

Discrete element modeling on the crack evolution behavior of brittle sandstone containing three fissures under uniaxial compression

Based on experimental results of brittle, intact sandstone under uniaxial compression, the micro-parameters were firstly confirmed by adopting particle flow code(PFC2D). Then, the validation of the simulated models were cross checked with the experimental results of brittle sandstone containing three parallel fissures under uniaxial compression. The simulated results agreed very well with the experimental results, including the peak strength, peak axial strain, and ultimate failure mode. Using the same microparameters, the numerical models containing a new geometry of three fissures are constructed to investigate the fissure angle on the fracture mechanical behavior of brittle sandstone under uniaxial compression. The strength and deformation parameters of brittle sandstone containing new three fissures are dependent to the fissure angle. With the increase of the fissure angle, the elastic modulus, the crack damage threshold,and the peak strength of brittle sandstone containing three fissures firstly increase and secondly decrease. But the peak axial strain is nonlinearly related to the fissure angle. In the entire process of deformation, the crack initiation and propagation behavior of brittle sandstone containing three fissures under uniaxial compression are investigated with respect tothe fissure angle. Six different crack coalescence modes are identified for brittle sandstone containing three fissures under uniaxial compression. The influence of the fissure angle on the length of crack propagation and crack coalescence stress is evaluated. These investigated conclusions are very important for ensuring the stability and safety of rock engineering with intermittent structures.