单轴压缩下横观各向同性岩石破裂过程的数值模拟

采用基于细观损伤力学基础上开发的RFPA2D数值模拟软件,用2种不同的岩石材料来组成7个不同岩层倾角的横观各向同性的岩石试件,通过单轴加载数值模拟试验,模拟横观各向同性岩石渐进破裂的整个过程,分析了岩层与最大主应力之间的倾角和强度之间的关系,讨论了不同岩层倾角的横观各向同性岩体的不同破裂模式及其破坏准则。     

A modi?ed failure criterion for transversely isotropic rocks

A modified failure criterion is proposed to determine the strength of transversely isotropic rocks. Me-chanical properties of some metamorphic and sedimentary rocks including gneiss, slate, marble, schist, shale, sandstone and limestone, which show transversely isotropic behavior, were taken into consider-ation. Afterward, introduced triaxial rock strength criterion was modified for transversely isotropic rocks. Through modification process an index was obtained that can be considered as a strength reduction parameter due to rock strength anisotropy. Comparison of the parameter with previous anisotropy in-dexes in literature showed reasonable results for the studied rock samples. The modified criterion was compared to modified Hoek-Brown and Ramamurthy criteria for different transversely isotropic rocks. It can be concluded that the modified failure criterion proposed in this study can be used for predicting the strength of transversely isotropic rocks.     

层理面特性对砂岩断裂力学行为的影响研究

层理弱面对层状岩石的力学特性影响较显著，为了研究层理面特性对岩石断裂力学特性的影响，开展了具有不同层理方向的砂岩试样三点弯试验，探讨了砂岩断裂韧度及断裂模式的各向异性。之后基于有限元中的黏聚单元建立了数值模型，采用数值模拟方法研究了层理面强度对各层理角度试样断裂力学行为的影响规律。结果表明：层理方向影响下砂岩的断裂韧度及模式存在各向异性；同一层理方向试样的断裂韧度随层理面强度的增大而增大，且试样的层理面与加载方向夹角越小，断裂韧度受层理面强度变化影响越明显；试样的断裂模式不仅与层理面强度有关，还受层理倾角的控制，层理面与加载方向夹角θ = 0o试样断裂模式基本不受层理面强度影响，θ = 30o试样主要沿层理面张拉或剪切破坏，且沿层理面的破裂长度随层理面强度的降低逐渐增大；层理面强度较大时，θ = 45o试样主要沿层理面张拉破坏，θ = 60o～90o试样主要以贯穿层理的张拉破坏为主；层理面强度较小时，θ = 45o～90o试样均以沿层理面的剪切破坏为主，其中θ = 45o试样沿层理剪切长度最大。另外，通过数值模拟结果分析了层理面强度及方向对试样的起裂角及裂纹扩展路径产生的影响。该研究成果可作为层状岩石断裂力学理论的有益补充。     

Slope failure analysis considering anisotropic characteristics of foliated rock masses

Weak planes affect the strength and deformational behaviors of rock slopes, and the anisotropic characteristics of rock mass should be considered in slope stability analysis. Effects of joint plane orientations on failure mechanism and strength response of inherently anisotropic rock samples were firstly investigated. The specimens with various orientations of joints were evaluated under uniaxial compression, Brazilian tensile, and direct shear tests. By treating the foliated rock as transversely isotropic materials, the relevant elastic constants and strength parameters were obtained from experimental results. The slope damage zone was then investigated using Comsol Multiphysics code based on Hoffman criterion. It is indicated that the failure mechanism and strength response depend highly on the inclination of specimens with respect to the loading direction. For disks with the same inclination angle, the value of tensile strength has an increasing trend with the total fracture length. Numerical results show that partial slope mass failed in single slope and no large-scale landslide occurred. The failure pattern in numerical results agrees well with the field observations. The cooperation between the experimental results and the numerical results allows an in-depth analysis of the experimental results and thus better understanding the dominant effect of joints on the deformation and failure of rock mass.     

层状岩体地下洞室破坏模式数值模型及验证

基于层面法向局部坐标系中的横观各向同性弹性本构模型,根据层面产状方位角,推导出整体坐标系中弹性应力的转换方法。针对层状岩体的各向异性特征,将岩体的破坏模式细分为5种,分别对每一种破坏模式建立各自的判别函数和屈服准则,根据塑性流动正交法则推导出相应的应力修正迭代计算方法。通过单轴、三轴压缩数值试验,对层状岩体的强度和变形性质进行研究,分析层面倾角和围压对层状岩体力学行为的影响。将所建模型与FLAC3D自带的横观各向同性弹性模型以及遍布节理模型进行对比,验证模型的正确性及精度。同时,将模型应用于三向等压理想圆形隧洞的开挖计算中,分析层状岩体破坏方式受层面方位角影响的变化规律,分析结果表明层面对围岩破坏模式和塑性区扩展方向起到了控制作用。     

Three-dimensional discrete element simulation of indirect tensile behaviour of a transversely isotropic rock

This paper presents the development of a three-dimensional discrete element model using flat-joint and smooth-joint contact models to investigate the effect of anisotropy on the tensile behaviour of slate, a transversely isotropic rock, under Brazilian testing from both macro and microscales. The effect of anisotropy is further realised by exploring the influence of foliation orientations (β and ψ) on the tensile strength, fracture pattern, microcracking and stress distribution of the transversely isotropic rock. The variation of tensile strength with foliation orientation is presented. The cross-weak-plane fracture growth observed in laboratory is reproduced, and the criterion for which to form is also given from the aspect of foliation orientation. Furthermore, the proportional variations of microcracks well account for the effects of foliation orientation on the tensile strength and failure pattern. Finally, it is found that the existence of weak planes increases both the heterogeneity and the anisotropy of stress distributions within the transversely isotropic rock, with the degree of influence varying with the foliation orientation.     

On the strength of transversely isotropic rocks

Accurate prediction of strength in rocks with distinct bedding planes requires knowledge of the bedding plane orientation relative to the load direction. Thermal softening adds complexity to the problem since it is known to have significant influence on the strength and strain localization properties of rocks. In this paper, we use a recently proposed thermoplastic constitutive model appropriate for rocks exhibiting transverse isotropy in both the elastic and plastic responses to predict their strength and strain localization properties. Recognizing that laboratory‐derived strengths can be influenced by material and geometric inhomogeneities of the rock samples, we consider both stress‐point and boundary‐value problem simulations of rock strength behavior. Both plane strain and 3D loading conditions are considered. Results of the simulations of the strength of a natural Tournemire shale and a synthetic transversely isotropic rock suggest that the mechanical model can reproduce the general U‐shaped variation of rock strength with bedding plane orientation quite well. We show that this variation could depend on many factors, including the stress loading condition (plane strain versus 3D), degree of anisotropy, temperature, shear‐induced dilation versus shear‐induced compaction, specimen imperfections, and boundary restraints.     

改进的遍布节理模型及其在层状岩体地下工程中的应用

基于显式有限差分程序FLAC3D,针对层状岩体建立了可以考虑其横观各向同性变形特性的遍布节理模型。通过FLAC3D程序的预留接口导入程序,将该改进的遍布节理模型中植入FLAC3D动态链接库。在此基础上进行了层状岩体变形与强度各向异性特性的研究,最后将该模型应用于龙滩水电站巨型地下硐室群的层状岩体围岩变形及破裂特征分析。研究表明,围岩变形主要表现为岩层同性面内的变形,其左右边墙变形不对称性主要受断层切割控制;围岩破坏型式以剪切破坏为主,其中岩体整体破坏受断层控制,表现为中低应力条件下的拉剪或压剪破坏;而开挖引起的岩层破坏受制于陡倾角层状岩体结构,表现为层间错动引起的剪切破坏。     

Fracture evolution of transversely isotropic rocks with a pre-existing flaw under compression tests based on moment tensor analysis

Acta Geotechnica - A numerical approach based on the particle discrete element theory is adopted to study the fracture evolution of transversely isotropic rocks with a pre-existing flaw under...     

A three-dimensional numerical investigation of the fracture of rock specimens containing a pre-existing surface flaw

Three-dimensional surface crack initiation and propagation in two kinds of heterogeneous rocks were numerically investigated via parallel finite element analysis using a supercomputer. Numerically simulated rock specimens containing a pre-existing flaw were subjected to uniaxial compression until failure. The initiation and propagation of wing cracks, anti-wing cracks, and shell-like cracks were reproduced by numerical simulations. The numerically simulated results demonstrate that the further propagation of wing cracks and shell-like cracks stop due to their wrapping (curving) behavior in three-dimensional spaces, even if the applied loads continue to increase. Furthermore, rock heterogeneity could significantly influence crack propagation patterns and the peak uniaxial compressive strengths of rock specimens. Moreover, anti-wing cracks only appeared in relatively heterogeneous rocks, and the peak uniaxial compressive strengths of the specimens were observed to depend on the inclination of the pre-existing flaw. Finally, the mechanism of surface crack propagation is discussed in the context of numerically simulated anti-plane loading tests, wherein it was identified that Mode III loading (anti-plane loading) does not lead to Mode III fracture in rocks due to their high ratio of uniaxial compressive strength to tensile strength. This finding could explain the lateral growth of an existing flaw in its own plane, which is a phenomenon that has not been observed in laboratory experiments.     

Numerical Modelling of the Anisotropic Mechanical Behaviour of Opalinus Clay at the Laboratory-Scale Using FEM/DEM

The Opalinus Clay (OPA) is an argillaceous rock formation selected to host a deep geologic repository for high-level nuclear waste in Switzerland. It has been shown that the excavation damaged zone (EDZ) in this formation is heavily affected by the anisotropic mechanical response of the material related to the presence of bedding planes. In this context, the purpose of this study is twofold: (i) to illustrate the new developments that have been introduced into the combined finite-discrete element method (FEM/DEM) to model layered materials and (ii) to demonstrate the effectiveness of this new modelling approach in simulating the short-term mechanical response of OPA at the laboratory-scale. A transversely isotropic elastic constitutive law is implemented to account for the anisotropic elastic modulus, while a procedure to incorporate a distribution of preferentially oriented defects is devised to capture the anisotropic strength. Laboratory results of indirect tensile tests and uniaxial compression tests are used to calibrate the numerical model. Emergent strength and deformation properties, together with the simulated damage mechanisms, are shown to be in strong agreement with experimental observations. Subsequently, the calibrated model is validated by investigating the effect of confinement and the influence of the loading angle with respect to the specimen anisotropy. Simulated fracture patterns are discussed in the context of the theory of brittle rock failure and analyzed with reference to the EDZ formation mechanisms observed at the Mont Terri Underground Research Laboratory.     

Assessments of Strength Anisotropy and Deformation Behavior of Banded Amphibolite Rocks

Assessment of strength anisotropy in transversely isotropic rocks has been one of the most challenging subjects in rock engineering. However, far too little attention has been paid to banded amphibolite rocks. This study aim to evaluate strength and deformation anisotropy behavior of banded amphibolite rocks. The dynamic mechanical tests including ultrasonic pulse test, uniaxial compressive strength, Brazilian test and deformability test were performed on drilled rock samples as a function of foliation plane angle (β = 0°, 30°, 60° and 90°). The results obtained have shown that the dynamic mechanical properties of amphibolite rocks have different values concerning banding plane. Compression and shear waves taken parallel to the foliation plane show highest values than those obtained in the other directions. Under uniaxial test, the banded amphibolite has a U-shaped anisotropy with maximum strength at β = 90° and minimum strength is obtained when β = 30°. Strength anisotropic index ranges between 0.96 and 1.47. It seems that the high range value of anisotropic index is mainly due to slight undulation of foliation planes, that being not perfectly straight. The results of elastic deformation test show that there is no clear dependence on microstructures characteristics of subtype-amphibolite rocks that controlling modulus “shape-anisotropy”. However, in this study, Young modulus values of amphibolite rocks with β follow both types of shape-anisotropy, “U-shape” and “decreased order-shaped”. Thus, this study recommended that further research be undertaken regarding the role of modulus “shape-anisotropy” within the same lithotype.     

A new Cosserat-like constitutive model for bedded salt rocks

Salt rocks are commonly used as geologic host rocks for storage of gas and crude oil, and are being considered for the disposal of radioactive waste. Different from the salt rock domes in many countries, the salt rock formations in China are usually laminar with many alternating layers, i.e. rock salt, anhydrite, and/or mudstone. Considering the unique stratigraphic characteristics of these salt rocks, a new Cosserat-like medium constitutive model is proposed in order to facilitate efficient modeling of the mechanical behavior of these formations. In this model, a new representative volume element, containing two different layers, is employed to simulate the compatibility of the meso-displacement between two different layers and also the bending effect. A new method for the deformation and failure analysis of bedded salt rocks is derived therefrom. Having the macro-average stresses, the conventional stresses in the different layers can be obtained in sequence. The conventional stresses can then be utilized in a routine way for the strength and failure analysis. For the initial numerical modeling, the new Cosserat-like medium is reduced to a transversely isotropic one. The simplified constitutive model for layered media is then implemented into FLAC3D codes. A test sample validates that the results by using the numerical model are in good agreement with that by using the built-in model, and the mesh size for the new model is reduced greatly. Finally, an application for the stability of oil storage caverns in deep thinly bedded salt rocks is carried out. The effects on convergence of storage caverns and on the failure of surrounding rock due to the presence of the mudstone interlayers (hard phase) are discussed in detail. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical studies on the failure process and associated microseismicity in rock under triaxial compression

In this paper, firstly the mesoscopic elemental mechanical model for elastic damage is developed and implemented into the rock and tool interaction code (R-T^2D). Then the failure processes of a heterogeneous rock specimen subjected to a wide variety of confining pressures (0–80 MPa) are numerically investigated using the R-T^2D code. According to the simulated results, on the one hand, the numerical simulation reproduced some of the well-known phenomena observed by previous researchers in triaxial tests. Under uniaxial compression, rock failure is caused by a combination of axial splitting and shearing. Dilatancy and a post-failure stage with a descending load bearing capacity are the prominent characteristics of the failure. As the confining pressure increases, the extension of the failed sites is suppressed, but the individual failure sites become dense and link with each other to form a shear fracture plane. Correspondingly, the peak strength, the residual strength and the shear fracture plane angle increase, but the brittleness decreases. When the confining pressure is high enough, the specimen behaves in a plastic manner and a narrow shear fracture plane leads to its failure. The prominent characteristics are volume condensation, ductile cataclastic failure, and a constant load bearing capacity with increasing strain. On the other hand, the numerical simulation revealed some new phenomena. The highest microseismicity events occur in the post-failure stage instead of the maximal stress, and most of the microseismicity energies are released in the failure localization process. As the confining pressure increases, the microseismicity events in the non-linear deformation stage increase dramatically and the ratio between the energies dissipated at the non-linear deformation stage and those dissipated in the whole loading process increases correspondingly. Therefore, it is concluded that the developed mesoscopic elemental mechanical model for elastic damage is able to reproduce accurately the failure characteristics in loading rock specimens under triaxial conditions, and the numerical modelling can furthermore obtain some new clarifications of the rock fracture process.     

深埋顺层偏压隧道围岩破坏机理及规律研究——以郑万线某隧道为例

处于薄—中层倾斜层状岩体中的深埋隧道常会产生地质顺层偏压的问题,导致隧道局部塌方、偏压变形及支护结构破坏。本文以郑万线某隧道为例,采用理论分析、数值模拟方法对深埋顺层隧道的破坏机理及不同结构面参数下的破坏规律展开了研究。研究结果表明:（1）深埋顺层偏压隧道洞周围岩将根据其切向应力与结构面夹角的不同发生岩层拉裂破坏、结构面剪切破坏及岩体自身破坏,其中切向应力与结构面平行位置,即反倾侧拱腰及顺倾侧拱脚位置主要发生拉裂破坏,此处围岩塑性区范围最广,围岩位移最大,围岩处于极不稳定状态;（2）顺层偏压隧道的破坏规律与结构面强度参数有直接关系,围岩塑性区范围及围岩位移均随着结构面摩擦角的增大而降低,且降低趋势逐渐放缓,当结构面摩擦角达到岩体摩擦角后,结构面摩擦角继续增加对围岩稳定性影响较小;（3）围岩塑性区及围岩位移场偏压分布特征随结构面倾角的变化而整体旋转,且对于隧道底部而言,结构面最不利倾角为0°,此时隧底最大上鼓量大于其他倾角下的最大上鼓量;对于隧道拱部而言,最不利倾角为40°,此时洞周最大收敛值大于其他倾角下的最大收敛值,最不利位置位于反倾侧拱腰。     

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

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

Measurement of Indirect Tensile Strength of Anisotropic Rocks by the Ring Test

Summary This paper presents a new approach, combined with the Boundary Element Method (BEM) analysis and the diametrical compression on a thin disc with a small central hole, referred to as the ring test, for determining the indirect tensile strength of anisotropic rocks. The stress distribution around the hole can be successfully obtained by the proposed single-domain BEM. The complex variable function method was used for conveniently computing the tractions and displacements of a two-dimensional anisotropic body. If we assume that the tensile strength is given by the maximum absolute value of stress in the direction perpendicular to the loaded diameter at the intersection of loaded diameter and the hole, then from the failure load recorded by laboratory testing of ring (disc), the indirect tensile strength of rocks could be obtained. A marble from Hualien (Taiwan) with clearly black-white foliation, which was assumed to be transversely isotropic, was selected to conduct both ring tests and Brazilian tests for evaluating the tensile strength. The variation of the marble tensile strength with the inclination angle of foliation and with the hole size was also investigated. In general, the tensile strength of anisotropic rocks determined by ring test is not a constant, but depends on the elastic properties of rocks, the angle between the planes of rock anisotropy and the loading direction, the diameter of the central hole, and the contact condition of loading.     

基于连续-非连续方法的裂隙破坏与相互作用研究

岩体内裂隙等非连续结构面对岩体的强度及变形等力学特性有着显著的影响,研究岩体裂隙起裂、扩展、相互作用和贯通机制,对工程岩体力学行为的表征和工程性能的评价十分重要。本文基于连续介质力学模型的离散元方法,通过考虑裂隙分布、模型加载条件及其与裂隙产状的关系,建立了一系列裂隙力学计算模型,研究了不同模型裂隙扩展演化特征和岩体破裂机制,分析了岩体裂隙扩展规律及其对岩体破坏路径和强度的影响,研究结果表明:（1）裂隙岩体模型加载条件下的破坏起裂点、最终贯通破坏特征及损伤分布受控于裂隙的产状及其与最大主压应力取向角度大小及围压大小。（2）裂隙弱面走向与最大主压应力取向斜交时,裂隙弱面在加载条件下其端部裂隙扩展、贯通破坏表现比较明显,反之,当裂隙弱面走向与最大主压应力取向一致时,裂隙弱面被动影响裂隙模型内新生裂隙的萌生、扩展和贯通模式,自身未出现新的扩展破坏。（3）裂隙数目的增多和围压的增大会显著增加模型内部剪切裂缝的数量和模型破坏后的破碎程度,模型内部的损伤区域主要围绕破裂面呈滑移线型交叉分布,非破裂面区域损伤呈条带状X型分布。（4）裂隙弱面走向与最大主压应力取向斜交时,裂隙对岩体模型强度的弱化程度高于裂隙弱面走向与最大主压应力取向一致的情况,而裂隙模型破坏后的残余强度则正好相反。     

Failure pattern of brittle rock governed by initial microcrack characteristics

Fracturing processes of rock are simulated using a micro- and fracture mechanical-based numerical modeling approach. The numerical model considers material heterogeneity and initial microcrack distributions. The gradual formation of macroscopic fractures by coalesced microcracks is successfully reproduced. Distinct failure modes are observed in the model under different loadings. Agreement is shown between the numerical results and laboratory observations. The influence of microcrack orientations on the fracture patterns is quantified by the numerical models. Possible mechanisms describing the splitting failure of hard brittle rocks under uniaxial compression are proposed and discussed.

     

Determining elastic constants of transversely isotropic rocks using Brazilian test and iterative procedure

The elastic constants of rocks are the basic parameters for rock mechanics, and play a very important role in engineering design. There are many laboratory methods to determine the elastic constants of transversely isotropic rocks, and the Brazilian test is a popular method. This paper presented a method combination of the Brazilian test, back calculation, and iterative procedure to evaluate the five independent elastic constants of transversely isotropic rocks in laboratory. The strain data at the centre of discs were obtained using Brazilian test. The stresses at the centre of discs were computed using numerical programs. By using back calculation, the temporary elastic constants were computed after the stresses and stains were substituted into elastic mechanics equations. After iterative procedure, the convergent values of the elastic constants can be obtained. One numerical example and three experimental cases were proposed to show the applicability of this method. The convergent values of the five independent elastic constants can be obtained in no more than 10 iterative cycles. The results coming from numerical analysis method exhibited satisfactory outcome in accordance with those of generalized reduced gradient method. The merits of this method include convenient specimen preparation of the Brazilian test, simple iterative procedure, and readily available commercially numerical programs, so that this method can be easily popularized in research and engineering analysis. Copyright © 2007 John Wiley & Sons, Ltd.