SHPB冲击作用下层状千枚岩多尺度破坏机理研究

通过分离式霍普金森杆对层状千枚岩施加动态载荷,得到不同层理倾角下层状千枚岩的动态抗压强度与宏观破坏模式。采用三维激光仪获得断裂面细观形貌,引入分形几何定量计算断口面粗糙度;结合SEM观察到的微观尺度下不同层理倾角断口破坏机理,分析了不同层理倾角下层状岩石的动态破坏机制。研究结果表明:动态压缩下层理弱面对岩石的抗压强度影响较大;不同层理倾角千枚岩的断口形貌分形维数随层理倾角增大呈U型变化;从强度与裂纹扩展两方面考虑层理弱面对层状岩石破坏特征的影响,对于层理倾角为0°的试样,强度由岩石基质控制,但层理弱面仍对岩石破坏的裂纹分布与走向产生较大影响;对于层理倾角为22.5°的试样,强度与裂纹走向受岩石基质与层理弱面共同控制;对于层理倾角为45°~67.5°的试样,强度与裂纹走向受层理弱面控制;而对于层理倾角为90°的试样,动态抗压强度受岩石基质的影响较大,在层理弱面较早形成纵向宏观裂纹,导致该层理弱面角度下裂纹受层理弱面的影响较大。     

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.     

Predicting variability in the dynamic failure strength of brittle materials considering pre-existing flaws

We perform two-dimensional dynamic fracture simulations of a specimen in biaxial tension, incorporating various distributions of pre-existing microcracks. The simulations consider the spatial distribution of flaws while modeling the discrete failure processes of crack interactions and coalescence, and predict the macroscopic variability in failure strength. The model quantitatively predicts the effect (on the dynamic failure strength) of different shapes of the flaw size distribution function, the random spatial distribution of flaws, and the random local resistance to crack growth (i.e. strength) associated with each flaw. The effect of changing material volumes on the variability in failure strengths is also examined in relation to the flaw size distribution. The effect of loading rate on the variability in failure strengths is presented in a form that will enable improved constitutive modeling using non-local formulations at the continuum scale.     

岩石紧凑拉伸断裂过程及其尺寸效应的三维数值模拟

采用最近开发的三维岩石破裂过程分析软件RFPA3D模拟单边裂纹紧凑拉伸断裂过程。试验中五个不同尺寸的岩样具有相同的力学性质参数分布,模拟结果得到了裂纹扩展中的应力场、位移场和声发射的空间分布以及单边裂纹扩展贯通的过程。单边裂纹拉伸断裂的路径是一个复杂的空间三维曲面,三维裂纹比二维裂纹更为复杂。分析了岩石试样的峰值强度和试样尺寸之间的关系。随着岩样尺寸的增加,峰值强度逐渐减小,并且延性破坏特征更加明显,模拟结果满足岩石的尺寸效应规律。最后模拟了三组不同均匀性的试样拉伸破坏过程,结果表明细观上的非均匀性对岩石尺寸效应有很大影响,随着非均匀性的增加,岩石宏观强度随之提高,即使在均匀材料中一样存在尺寸效应。     

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).     

Loading rate and confining pressure effect on dilatancy,acoustic emission,and failure characteristics of fissured rock with two pre-existing flaws

Investigating the dilatancy, acoustic emission and failure characteristics of fissured rock are significant to ensure their geotechnical stability. In this paper, the uniaxial and triaxial compression experiments with AE monitoring under different loading rates were carried out on fissured rock specimens with the same geometrical distribution of two pre-existing flaws. The dilatancy and AE activity of these specimens were discussed, and the effects of the confining pressure and loading rate on the mechanical parameters and failure characteristics were analyzed. The results show that the exponential strength criterion is more suitable than the Mohr–Coulomb strength criterion to characterize the strength characteristics of fissured rock. The crack evolution and failure characteristics of fissured rock specimens are more complicated than those of intact rock specimens. The failure characteristics of the fissured rock follow the tensile shear coalescence model, crack branching occurs with increasing the loading rate, and the multi-section coalescence model is verified with increasing the confining pressure. The phenomena of stress drop and yield platform usually occur after the dilatancy onset, the specimen does not fail instantaneously, and the propagation and coalescence of cracks cause a sharp increase in the AE signals, circumferential strain, and volumetric strain.     

Experimental study on cracking behaviour of moulded gypsum containing two non-parallel overlapping flaws under uniaxial compression

Failure of rock mass that is subjected to compres-sive loads occurs from initiation, propagation, and linkage of new cracks from preexisting fissures. Our research inves-tigates the cracking behaviour and coalescence process in a brittle material with two non-parallel overlapping flaws using a high-speed camera. The coalescence tensile crack and tensile wing cracks were the first cracks to occur from the pre-existing flaws. The initiation stresses of the primary cracks at the two tips of each flaw were simultaneous and decreased with reduced flaw inclination angle. The following types of coalescence cracks were identified between the flaws: pri-mary tensile coalescence crack, tensile crack linkage, shear crack linkage, mixed tensile-shear crack, and indirect crack coalescence. Coalescence through tensile linkage occurred mostly at pre-peak stress. In contrast, coalescence through shear or mixed tensile-shear cracks occurred at higher stress. Overall, this study indicates that the geometry of preexisting flaws affect crack initiation and coalescence behaviour.     

基于TCK模型的非贯通节理岩体动态损伤本构模型

提出在岩体动态损伤本构模型中应同时考虑宏、细观缺陷;基于能量原理和断裂力学理论推导得出了同时考虑节理几何及力学特征的宏观损伤变量(张量)的计算公式;基于综合考虑宏、细观缺陷的复合损伤变量(张量)及完整岩石动态损伤Taylor-Chen-Kuszmaul(TCK)模型,建立了相应的单轴压缩下节理岩体动态损伤本构模型;利用该模型讨论了节理内摩擦角及节理长度对岩体动态力学特性的影响规律。研究表明,试件动态峰值强度随着节理内摩擦角的增大而增大,随着节理长度的增加而减小。     

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.     

Experimental investigation on strength and failure behavior of pre-cracked marble under conventional triaxial compression

Fractures in natural rocks have an important effect on the strength and failure behavior of rock mass, which are often evaluated in rock engineering practice. The theoretical evaluation of mechanical behavior of fractured rock mass has no satisfactory answer due to the role of confining pressure and crack geometry. Therefore, in this paper, conventional triaxial compression experiments were carried out to study the strength and failure behavior of marble samples with two pre-existing closed cracks in non-overlapping geometry. Based on the experimental results of a number of triaxial compression tests, the effect of crack coalescence on the axial supporting capacity and deformation property were investigated with different confining pressures. The results show that intact samples and flawed samples (marble with pre-existing cracks) have different deformation properties after peak stress, which change from brittleness to plasticity and ductility with the increase of confining pressure. The peak strength and failure mode are found depending not only on the geometry of flaw, but also on the confining pressure. The strength of flawed samples shows distinct non-linear behavior, which is in a better agreement with non-linear Hoek–Brown criterion than linear Mohr–Coulomb criterion. For a kind of rock that has been evaluated as a Hoek–Brown material, a new evaluation criterion is put forward by adopting optimal approximation polynomial theory, which can be used to confirm more precisely the strength parameters (cohesion and internal friction angle) of flawed samples. For intact samples, the marble leads to typical shear failure mode with a single fracture surface under different confining pressures, while for flawed samples, under uniaxial compression and a lower confining pressure (σ₃ = 10 MPa), tests for coarse and medium marble (the coarse and medium refer to the grain size) exhibit three basic failure modes, i.e., tensile mode, shear mode, and mixed mode (tensile and shear). Shear mode is associated with lower strength behavior. However, under higher confining pressures (σ₃ = 30 MPa), for coarse marble, the axial supporting capacity is not related to the geometry of flaw. The friction among crystal grains determines the strength behavior of coarse marble. For medium marble, the failure mode and deformation behavior are dependent on the crack coalescence in the sample. The present research provides increased understanding of the fundamental nature of rock failure under conventional triaxial compression.     

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.     

Modeling dynamic brittle behavior of materials with circular flaws or pores

Compressive failure of brittle materials is driven primarily by crack growth from pre-existing flaws in the material. These flaws, such as grain boundaries, pores, preexisting cracks, inclusions and missing grains, are randomly spaced and have a range of possible shapes and sizes. The current work proposes a micromechanics-based model for compressive dynamic failure of brittle materials with circular pore flaws, which incorporates both the number density and the size distribution of flaws. Results show that the distribution of flaw sizes is very important, particularly at moderate strain rate, since analyses based solely on the mean flaw size overpredict strength. Therefore, in order to increase dynamic strength at low to moderate strain rates, it is most effective to control the presence of large flaws. At very high strain rates, however, crack growth is activated even in small flaws and therefore controlling the total number density rather than the size of the flaws is effective for increasing dynamic strength. Finally, the model shows that neglecting very small flaws in the pore population may not have significant effects on the results in many cases, suggesting that the model is a useful tool for identifying a minimum resolution required for experimental characterization of microstructure.     

Particle flow study on strength and meso-mechanism of Brazilian splitting test for jointed rock mass

A discrete element method （DEM） called particle flow code （PFC2D） was used to construct a model for Brazilian disc splitting test in the present study. Based on the experimental results of intact Brazilian disc of rock-like material, a set of micro-parameters in PFC2D that reflected the macro-mechanical behavior of rock-like materials were obtained. And then PFC2D was used to simulate Brazilian splitting test for jointed rock mass specimens and specimen containing a central straight notch. The effect of joint angle and notch angle on the tensile strength and failure mode of jointed rock specimens was detailed analyzed. In order to reveal the meso-mechanical mechanism of crack coalescence, displacement trend lines were applied to analyze the displacement evolution during the crack initiation and propagation. The investigated conclusions can be described as follows. （1） The tensile strength of jointed rock mass disc specimen is dependent to the joint angle. As the joint angle increases, the tensile strength of jointed rock specimen takes on a nonlinear variance. （2） The tensile strength of jointed rock mass disc specimen containing a central straight notch distributes as a function of both joint angle and notch angle. （3） Three major failure modes, i.e., pure tensile failure, shear failure and mixed tension and shear failure mode are observed in jointed rock mass disc specimens under Brazilian test. （4） The notch angle roles on crack initiation and and joint angle play important propagation characteristics of jointed rock mass disc specimen containing a central straight notch under Brazilian test.     

基于CT扫描的花岗岩三维数值试件重构模型及应用

本文以花岗岩试件为研究对象,借助于先进的CT技术获取试件内部结构切片图像,利用数字图像处理技术实现了基于CT图像的花岗岩试件结构的识别和表征,基于矢量化方法和岩石破裂过程分析系统RFPA3D建立了花岗岩三维数值试件,对花岗岩试件单轴压缩破裂过程进行了数值模拟。通过数值试验结果与物理实验结果对比发现,基于CT扫描的花岗岩数值试件的破裂模式与室内实验得到的破裂模式相似,为深入研究岩石、复合材料等力学特征提供了一种可行方法。     

Foam mechanics: nonlinear response of an elastic 3D-periodic microstructure

The compressive response of a 3D open-cell foam with periodic tetrakaidecahedral cells is studied through combined theoretical and numerical efforts. Under compressive loading the response is characterized by an extended load plateau following the relatively sharp rise to a maximum load. Several processes of loading have been simulated numerically using appropriately nonlinear kinematics. The onset of failure under macroscopic loading conditions is shown to be the reason of the load plateau. A failure surface is defined in macroscopic stress space by the onset of the first buckling-type instability encountered along proportional load paths. The analysis is carried out through two methods. The first one consists in increasing specimen size with periodic boundary conditions leading to the termed microfailure surface. The second one consists in considering both periodic and nonperiodic displacements variations on a minimum unit cell. The resulting failure surfaces are shown to coincide. Moreover, the postbuckling analysis has been carried out for two particular loadings: the uniaxial compression and the uniaxial deformation.     

Tensile failure of micro-concrete: from mechanical tests to FE meso-model with the help of X-ray tomography

In this work, concrete is studied at meso-scale (aggregates, macro-pores and mortar matrix), where the local failure mechanisms are known to drive the macroscopic behaviour of the material. In order to highlight the impact of the mechanical and morphological properties of each phase (along with their interfaces), micro-concrete specimens are prepared with rather small dimensions compared to the size of the heterogeneities. X-ray tomography is used to reliably obtain the morphology of the heterogeneous meso-structure, which is then given as an input to a 3D FE meso-model with enhanced discontinuities. A uniaxial tensile numerical simulation is performed as a first application. To validate the numerical model, a uniaxial tensile test of the same micro-concrete specimen is performed inside the X-ray scanner and the in-situ evolution of the micro-structure is followed. Thus, both a direct validation of the model and a valuable insight of the 3D fracture mechanisms while the load progresses are obtained. After identification of the numerical parameters, comparison of experimental and numerical results reveals the capability of the meso-model to reproduce the actual material response (in terms of macroscopic strength, Young’s modulus and fracture patterns), with the explicit representation of the meso-scale heterogeneities being its key feature. To further challenge the meso-model, a new morphology coming from an X-ray scan of another characteristic micro-concrete specimen is introduced and its macroscopic behaviour is computed without a priori numerical identification. Starting from an X-ray scan in meso-scale, it is shown that the 3D meso-model is capable to predict the macroscopic behaviour and the failure patterns of the material.

     

各向异性编织CMC弯曲断裂失效模型

通过对缺口弯曲断裂试验的研究发现，编织CMC存在两种断裂失效模式：界面主导的和纤维束主导的失效模式．针对这两种模式建立了界面失效模型和纤维束断裂失效模型．界面失效模型用半经验方法处理弱界面和大编织角度的断裂强度问题；纤维束断裂失效模型采用解析方法处理较强界面和小角度的断裂强度问题，为含有缺口试件的弯曲强度预报打下了初步的理论基础．     

基于数字图像的不同倾角节理灰岩破裂模式及细观尺度破裂过程研究

为研究节理倾角对灰岩破裂模式及破裂过程的影响,使用数字图像对灰岩的细观非均匀性进行表征,采用岩石破裂过程分析系统(RFPA2D-DIP)对不同倾角节理灰岩的细观破裂过程及宏观破坏模式进行了研究.结果表明,细观结构对各倾角节理试样的力学特性和最终破坏模式有重要影响,含节理灰岩的弹性模量及抗压强度具有明显的各向异性,随着节...