基于数值流形法的三维裂纹扩展研究

基于数值流形法（NMM）进行三维裂纹扩展分析研究,编写了相应的C++程序。充分发挥数值流形法在非连续变形分析领域的优势,不需要扩展有限元中的水平集和阶跃函数概念,应用数值流形法计算结果分析裂纹尖端线的破坏状况,对已有的非局部求迹法和三角形推进法进行简化和扩展,提出一种简化算法确定最终的新裂纹扩展面。应用简化算法对水平钱币型裂纹扩展问题进行数值模拟和对比分析研究。计算结果表明,基于数值流形法的三维裂纹扩展算法是可行的,采用简化处理之后,克服原有算法的不足,大大减少了新裂纹尖端线和新裂纹面的数目,降低了计算网格对新裂纹面的影响,提高了计算效率和扩大了应用范围。     

Influence of inclination angles and confining pressures on mechanical behavior of rock materials containing a preexisting crack

To deeply understand the cracking mechanical behavior of brittle rock materials, numerical simulations of a rock specimen containing a single preexisting crack were carried out by the expanded distinct element method (EDEM). Based on the analysis of crack tips and a comparison between stress- and strain-based methods, the strain strength criterion was adopted in the numerical models to simulate the crack initiation and propagation processes under uniaxial and biaxial compression. The simulation results indicated that the crack inclination angle and confining pressure had a great influence on the tensile and shear properties, peak strength, and failure behaviors, which also showed a good agreement with the experimental results. If the specimen was under uniaxial compression, it was found that the initiation stress and peak strength first decreased and then increased with an increasing inclination angle α. Regardless of the size of α, tensile cracks initiated prior to shear cracks. If α was small (such as α ≤ 30°), the tensile cracks dominated the specimen failure, the wing cracks propagated towards the direction of uniaxial compression, and the propagation of shear cracks was inhibited by the high concentration of tensile stress. In contrast, if α was large (such as α ≥ 45°), mixed cracks dominated the specimen failure, and the external loading favored the further propagation of shear cracks. Analyzing the numerical results of the specimen with a 45° inclination angle under biaxial compression, it was revealed that lateral confinement had a significant influence on the initiation sequence and the mechanical properties of new cracks.     

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

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

含折线型裂隙砂岩试件翼型裂纹起裂与扩展机制研究

岩体内部赋存的裂隙很多表现为折线型,为探究这类岩体的断裂机制,制备含折线型裂隙砂岩试件并对其进行单轴压缩试验。采用数字图像相关(DIC)方法计算加载过程中的变形场演化,根据新生裂纹两侧的位移差异识别裂纹类型;运用扩展有限元法(XFEM)模拟断裂过程,根据应力分布特征解释翼型裂纹起裂与扩展机制。DIC计算结果表明,新生裂纹处出现应变局部化带,裂纹两侧发生相对分离;含直线型和折线型裂隙砂岩试件的翼型裂纹分别萌生于预制裂隙端部以及折角处,这是因为裂隙几何形态会改变拉应力集中位置;含折线型裂隙砂岩试件的起裂应力小于含直线型裂隙砂岩试件,这是因为相同加载条件下前者的最大拉应力值更大;这2类试件的裂纹扩展均是由于裂纹尖端集中的拉应力引起的,裂纹依然呈张开状态;裂隙几何形态未改变试件的最终破坏模式,均表现为对角剪切破坏。     

基于数值流形法的重力坝多裂纹扩展研究

传统数值流形法（NMM）在处理非连续变形问题时,仅限于几何构型不发生破坏的情况。针对这一不足,通过在裂纹尖端附近的物理片上增加用于模拟应力奇异性的增强位移函数,进一步发展了可用于几何构型破坏的扩展的高阶NMM。然后,将其应用到重力坝由连续到非连续的破坏过程分析中。首先,针对一含单裂纹的重力坝模型进行了敏感性分析,结果表明,在不同的扩展长度或网格密度下,其扩展路径基本相同且与文献结果保持一致。进而在此模型基础上又开展了多裂纹扩展分析,结果仅一条主导裂纹发生扩展,与文献结果基本一致。最后,针对印度的Koyna重力坝,通过设置不同的漫顶高度研究了其裂纹扩展路径的变化。结果表明,随着漫顶高度的增大,裂纹扩展路径逐渐趋向于水平方向扩展,而且坝体抵抗破坏的能力逐渐减弱。总体表明,NMM在求解实际工程问题时具有很好的数值稳定性和鲁棒性。     

Crack propagation trajectories for rocks under mixed mode I-II fracture

Propagation of a crack in engineering materials including rocks can cause failure. Knowledge of the stress state under which a crack can propagate, and the trajectory it may follow during its growth are thus very important for the stability of rock masses/materials and for the safe design of structures in/on rocks. In this paper, the crack initiation angle and subsequent crack propagation path are experimentally investigated for limestone rock specimens. This investigation was conducted under various mixed mode I-II loading conditions, including pure mode-I and pure mode-II.This study includes conducting diametrical compression tests on notched Brazilian disk specimens. Moreover, the effect of confining pressure and temperature on crack initiation and propagation were also studied. The experimental results were compared with theoretical predictions of crack initiation angle. The results showed that limestone behaves in brittle fashion, and the effects of confining pressure and temperature on failure trajectories were not significant. Generally, the crack initiation angle can be predicted by the maximum tangential stress criterion. However, for notched Brazilian disk with high value of crack orientation with respect to loading direction, crack does not propagate from the tip of the crack. This important observation indicated that the tensile-strength failure can become more critical than the fracture-toughness failure.     

Experimental Analysis on Deformation Evolution and Crack Propagation of Rock Under Cyclic Indentation

Indentation-induced cracking may have many implications in rock drilling and cutting processes. In the present work, cyclic indentation tests have been conducted on Yunnan sandstones to investigate the deformation and failure process. Deformation fields on the sample surface have been obtained via the digital image correlation (DIC) technique. It has been found that the macroscopic force–indentation curve exhibits hysteresis and the residual displacements increase during the loading–unloading cycles. Furthermore, a median crack has been observed to nucleate, running almost parallel to the loading axis when the indentation force amounts to a critical value. Such a crack always results in the final failure of the sandstone samples. Experimental analysis reveals the deformation and failure mechanism of sandstone during cyclic indentation. By the DIC technique, the displacement distribution along the crack surface can be easily calculated and, hence, the fracture parameters, such as crack opening displacement and fracture toughness.     

牙轮钻头动静耦合碎岩机理及旋挖成桩应用

总结了在不同高温状态下冷却、不同加载速率及随机裂隙发育状态下花岗岩动态抗拉力学特性的变化规律;结合牙轮钻头孔底碎岩过程,分析了动静载荷耦合作用下岩石破碎的载荷一侵深特性曲线,认为动、静载荷耦合作用的加载点（即动载的施加点）应是在静载处于卸载阶段;并根据加载能量大小讨论了不同动静耦合工况下产生的岩石破裂深度及破碎体积,表明通过一定范围内增大静载荷及冲击力、预加静压对岩石进行预应力损伤、加载一卸载一加载的破碎循环模式,有利于高效碎岩及裂纹的发育。嵌岩桩基础工程实践表明,通过改造牙轮钻头等钻具结构形式及布齿方式,利用动静耦合加载方式及对钻头冷却处理,可实现牙轮钻头在微风化花岗岩高效钻进的目的,为牙轮钻头的旋挖钻进成桩应用提供了重要的技术支撑。     

冲击作用下的压头破岩机制研究

与静态岩石破碎过程相比,冲击作用下岩石的应力改变具有时间效应,应力波传播过程中表现出压、拉变化。基于损伤演化原理和有限元数值模拟方法,针对冲击荷载作用下的压头破岩机制进行了模拟分析。为排除边界上反射波的影响,黏弹性边界被纳入计算中。首先论证了黏弹性边界在均质和非均质介质中的计算精度,然后分析了冲击作用下不同均质度的岩石以及砂砾岩的响应规律,结果显示：在弹性情况下,压头与岩石接触边缘以及自由面附近是拉应力分布区,接触边缘拉应力最大。剪应力最高值并不位于接触面附近,而是离接触面有一定距离。较均质岩石主要呈现拉伸破坏模式,先出现赫兹裂缝,然后是径向裂缝和侧向裂缝,拉应力的产生成为诱发裂缝萌生和扩展的主因。当岩石均质度较低时,岩石的破坏形式呈现多元化,剪切破坏比重加大,表现为复杂的拉剪破坏模式。对于砾石粒径较大、含量较多的砂砾岩,砾石和基质的非均匀性不可忽略,冲击下破坏模式以绕粒环行和穿粒破坏为主。总体说来,对于岩石类准脆性材料,应力波传播过程中产生的拉应力是失稳的诱发和扩展的关键。     

Boundary settings for the seismic dynamic response analysis of rock masses using the numerical manifold method

Aiming to accurately simulate seismic dynamic response of rock masses using the numerical manifold method (NMM), boundary settings must be treated carefully. In this paper, 4 issues in boundary settings are investigated to improve the performance of NMM: (1) Nonreflecting boundaries including the viscous boundary and viscoelastic boundary are considered; (2) A free‐field boundary is incorporated into NMM to accurately simulate external source wave motion; (3) A seismic input boundary is considered, and the force input method is introduced; and (4) A static‐dynamic unified boundary is incorporated for the convenience of transforming displacement boundary into other types of boundaries, such as nonreflecting boundaries and seismic input boundary. Several benchmark problems are solved to validate the improved NMM. Simulation results agree well with analytical ones, indicating that the improved NMM is able to simulate seismic dynamic response of rock masses reliably and correctly.     

深埋隧洞开挖卸荷引起的围岩开裂特征及影响因素

通过理论分析和数值计算研究深埋地下隧洞开挖卸荷引起的围岩开裂,并分别对准静态和瞬态开挖缷荷引起的围岩开裂机制与开裂特征进行分析。采用双向受压条件下的压剪裂纹扩展模型和应力强度因子计算公式,分析了开挖面上岩体应力瞬态释放和围岩应力瞬态调整对围岩开裂过程的影响,并分别对准静态和瞬态卸荷引起的围岩开裂机制及影响因素进行探讨。研究结果表明,围岩开挖缷荷是深埋隧洞围岩发生开裂的重要原因之一,高地应力条件下围岩以剪切型断裂破坏为主;瞬态卸荷存在动态拉应力效应,开挖卸荷时间越短,引起的拉应力区及围岩开裂范围越大;围岩开裂深度及范围随着侧压力系数增加而增大,且开裂区近似成V型。深埋隧洞围岩开裂特征及影响因素研究,对进一步认识围岩开挖破坏的力学机制具有重要的意义。     

岩石材料裂纹扩展的改进近场动力学方法模拟

在传统键基近场动力学模型的基础上,提出可以反映岩石类材料应力-应变变化特点的改进型模型,弥补了传统近场动力学模型在模拟岩石类材料时无法反映岩石内应力随应变先增加再减小最后破坏这一特性的不足。运用改进型模型,对含有单裂纹试件在不同预制裂纹角度及拉伸荷载条件下的裂纹扩展过程进行了数值模拟,分析了开裂角度在不同预制裂纹角度及拉伸荷载条件下的变化。结果表明：预制角度及拉伸荷载对开裂角度都有影响,在试件周围拉荷载接近的情况下会获得较小的开裂角,并且在拉荷载一定的条件下开裂角会有先增加后减少的规律。同时模拟了含裂纹岩石在单轴压缩条件下的裂纹扩展情况,将数值模拟的结果与室内试验的结果对比,验证了所提出新模型的有效性。所提出的新模型可以较好地模拟岩石类材料在拉荷载条件下的力学特性,在岩石数值模拟方面有着较为广泛的应用前景。     

Frictional crack initiation and propagation analysis using the numerical manifold method

By employing both a physical mesh and a mathematical mesh to formulate a physical problem, the numerical manifold method (NMM) can lead to a very simple meshing task, which allows directly capturing the discontinuities across the crack surfaces without further incorporating unknowns to the related nodes through enrichment functions. These features enable the NMM to handle complex crack problems. In this study, based on the contact technique of the NMM and the incorporation of the Mohr–Coulomb crack initiation criterion, the effects of the friction and cohesion on the crack growth from a closed flaw (crack) under compression were investigated. A limited number of comparisons between the numerical results and the physical experiments show that with the Mohr–Coulomb crack initiation criterion, the NMM can not only accurately predict the pure tensile or pure shear crack growth, but the NMM can also satisfactorily predict the development of mixed shear–tensile crack types. Using a parametric analysis, the effects of the confining stress, the flaw inclination angle, the flaw friction angle and the material strengths on the crack development (crack initiation stress, crack initiation angle, crack type developed) have been investigated.     

一种测量岩石压缩裂纹扩展过程声波波速的连续测量方法

中低围压或单轴压缩荷载条件下,可以将脆性岩石材料的裂纹扩展过程划分为原生裂纹压密、新生裂纹起裂并稳定扩展、不稳定裂纹扩展和交互贯通以及峰后5个主要发展阶段。含预制裂纹试样的直接观测法、应变监测、声发射（AE）监测、声波波速测试、CT（computerized tomography）扫描、微电镜观察等间接监测方法被用来研究裂纹起裂和扩展过程。本文基于在中、低围压以及单轴压缩荷载条件下,脆性岩石材料首先产生张拉裂纹,并优先沿加载方向发育的规律,提出采用声波波速连续测量方法进行岩石加载过程波速的连续观测,有助于分析裂纹扩展过程中声波波速的变化,为研究裂纹扩展过程提供一种间接测量手段。     

弹塑性数值流形方法在边坡稳定分析中的应用

针对现有数值流形方法只能进行弹性计算的不足,建立了一个能够反映完整岩块弹塑性变形特征的本构模型,并借助VC++开发了内置该本构模型的弹塑性数值流形程序。利用该程序模拟了含单节理岩样的室内压缩试验,分析得到了其强度和变形特性,计算结果符合实际的物理现象,表明程序是正确有效的。考虑到数值流形方法本身能够有效模拟材料的不连续变形,新增的弹塑性分析功能又可以反映岩石的强度特性,将弹塑性数值流形程序应用于某含有不连续面的岩石边坡的稳定性分析。并结合锚杆单元的使用,对比分析了不同锚固方案的加固效果。程序提供的变形、应力等计算结果表明:预应力锚杆不仅可以防止不连续面发生剪切破坏,增强坡体的稳定性,限制塑性变形的发展;而且可以使不连续软弱层面对岩体变形的消极影响得以减弱,起到提高岩体整体性的作用。     

A relationship between tensile strength and loading stress governing the onset of mode I crack propagation obtained via numerical investigations using a bonded particle model

Numerical models based on the discrete element method are used to study the fracturing process in brittle rock‐like materials under direct and indirect tension. The results demonstrate the capacity of the model to capture the essential characteristics of fracture including the onset of crack propagation, stable and unstable crack growth, arrest and reinitiation of fracturing, and crack branching. Simulations of Brazilian indirect tension tests serve to calibrate the numerical model, relating macroscopic tensile strength of specimens to their micromechanical breakage parameters. A second suite of simulations reveals a linear relationship between the tensile strength of specimens and the loading stress for which mode I tensile crack propagation ensues. Based on these results, a crack initiation criterion for brittle materials is proposed, prescribing the stressing conditions required to induce tensile failure. Such a criterion, if broadly applicable, provides a practical means to rapidly assess the failure potential of brittle materials under tensile loads.     

柱状节理玄武岩隧洞破坏模式及其力学机制模拟

柱状节理岩体由于其内部赋存大量的隐节理面,开挖卸荷后极易出现隐节理面开裂松弛等特征,导致其破坏模式异于一般岩体。其破坏模式主要受到异常发育的节理面和较高地应力的共同作用影响。由于柱状节理岩体节理面发育,岩体结构控制型破坏占主要部分,包括单临空面节理面滑移（塌方）、多临空面节理面滑移（塌方）、与错动带、断层等弱面相组合的坍塌等破坏模式;应力控制型破坏主要为河谷侧顶拱喷层开裂;应力-结构面型破坏主要为岩性交界处的节理岩体塌落等。柱状节理岩体表层主要发生柱内竖直隐节理面和柱间节理面的拉破坏,而围岩内部的柱状节理主要发生柱间节理面的剪切破坏。因此,现场柱状节理的支护也应主要包括两个方面：以喷射混凝土钢纤维来阻止柱状节理岩体表层的张性开裂,以系统锚杆来控制柱状节理岩体内部的剪切破坏。     

Investigation on strength and stability of jointed rock mass using three‐dimensional numerical manifold method

This paper proposes a numerical model for jointed rock masses within the 3‐D numerical manifold method (NMM) framework equipped with a customized contact algorithm. The strength of rock sample containing a few sets of discontinuities is first investigated. The results of models with simple geometries are compared with the available analytical solutions to verify the developed computer code, whereas models with complex geometries are simulated to better understand the fundamental behavior and failure mechanism of jointed rock mass. Furthermore, the stability of jointed rock mass in an underground excavation is studied, where rock failure process is determined by the 3‐D NMM simulation. The simulation results provide valuable guidance on excavation process design and stabilization design in rock engineering practice. Copyright © 2012 John Wiley & Sons, Ltd.     

基于3D-DIC的砂岩裂纹扩展及损伤监测试验研究

裂纹监测对岩石损伤演化的认识至关重要。为研究岩石裂纹扩展及损伤变形特性,开展了含不同倾角（0°～90°）预制裂隙的标准细黄砂岩样的单轴压缩试验。利用三维数字图像相关技术（3D-DIC）获取岩样三维空间坐标下的应变分布,并结合声发射从光学与声学的角度监测了裂纹的扩展演化。由此提出了一种裂纹主应变的计算方法,定量表征岩石劣化的损伤变量D值。最后,探讨了由声发射与损伤变量D值确定岩样特征强度的影响因素。结论如下：（1）裂纹主应变反映了岩样受荷过程中同源裂纹在时间上的变化速率与空间上的扩展趋势,能较好地表征岩石的开裂行为;（2）声发射适用于确定岩样的起裂应力,不适用于损伤应力的确定,损伤变量D值所确定的起裂应力滞后于声发射,但适用于损伤应力特征值的确定;（3）结合声发射与DIC技术确定的归一化起裂应力范围为0.63～0.94、归一化损伤应力的范围为0.83～0.99;（4）预有裂隙会影响岩石的材料力学性能。随着倾角的增加,岩石的起裂应力、损伤应力及峰值应力呈增长的趋势,由于难以形成局部应变场聚集,裂纹的萌生与起裂更加困难。结果表明,3D-DIC技术的利用可以提高对岩石开裂行为的理解,对岩石的损伤监测与判识更有重要的意义。     

高渗压条件下压剪岩石裂纹断裂损伤演化机制研究

探讨了高渗透压作用下压剪岩石裂纹的起裂规律及分支裂纹尖端应力强度因子的演变规律,建立了高渗压下裂隙岩体发生拉剪破坏的临界水压力值和初裂强度判据,分析了不同渗压作用下裂纹的扩展情况表明,渗透压的存在加剧了分支裂纹的扩展,高渗透压作用下分支裂纹扩展由稳定扩展变成不稳定扩展,并导致分支裂纹尖端岩桥剪切破坏,同时考虑分支裂纹的相互作用,建立了高渗透压作用下压剪岩石裂纹体岩桥剪切贯通的断裂破坏力学模型,最后依据裂隙岩体的损伤力学效应研究了岩体的初始损伤及损伤演化柔度张量,提出了高渗压下压剪岩石裂纹渐进破坏的损伤演化方程。该理论为定量研究高渗压下裂隙岩体的失稳破坏提供了理论依据。