岩石拉伸剪切破裂试验研究

岩石拉伸剪切破裂是一类特殊应力状态条件下的破裂形式,属于同时受垂直于破裂面的法向拉应力和平行于破裂面的剪应力作用的复合破裂模式。在研制的DSC-800电液伺服测控岩石拉伸剪切试验仪的基础上,进行了大量花岗闪长岩和砂岩的拉伸剪切试验,开展了配套的破裂断口三维激光扫描、扫描电子显微镜（SEM）、岩石物理力学性质试验、颗粒流离散元（PFC）数值模拟等相关试验,利用分形理论研究了岩石拉剪破裂面特征,研究了岩石拉剪-压剪全区破裂准则、剪切速率对岩石拉剪破裂强度的影响,采用颗粒流离散元研究了岩石拉剪破裂过程。研究结论如下:（1）岩石拉剪破裂面的宏观与微观分形维数即粗糙度随着拉应力的增加而增大;（2）岩石的微观断裂形式是拉伸破坏和剪切破坏的结合。当拉应力较小时,岩石的微观断裂形式主要表现为剪切破坏,并且随着拉应力的增加,岩石的拉伸破坏形式表现得更加明显;（3）岩石在拉伸剪切区的破裂拉应力与剪应力成线性负相关关系,在拉伸剪切应力区的岩石破裂线斜率比压缩剪切区大,岩石在拉伸剪切应力条件下比压缩剪切应力条件下容易破裂;（4）在岩石拉伸剪切条件下,剪切速率与剪切强度成非线性反相关关系,随着剪切速率的增加,岩石拉剪破裂面粗糙度增加;（5）建立了岩石拉伸剪切PFC数值试验模型,模拟了岩石拉伸剪切破裂过程中的力链演化以及剪切速率对拉剪破裂面粗糙度的影响,获得了与实验室试验一致的结果。     

Modelling and simulation of load transfer in reinforced soils: Part 1

Results of both triaxial and direct shear tests on reinforced soil samples performed by different investigators have shown that soil dilatancy and extensibility of the reinforcements have a significant effect on the generated tension forces in the inclusions. An appropriate soil--reinforcement load transfer model, integrating the effect of soil dilatancy and reinforcement extensibility is therefore needed to adequately predict forces in the inclusions under expected working loads. This paper present a load transfer model assuming an elastoplastic strain hardening behaviour for the soil and an elastic--perfectly plastic behaviour for the reinforcement. This model is used to analyse the response of the reinforced soil material under triaxial compression loading. A companion paper present the application of this model for numerical simulations of direct shear tests on sand samples reinforced with different types of tension resisting reinforcements. The model allows an evaluation of the effect of various parameters such as mechanical characteristics and dilatancy properties of the soil, extensibility of the reinforcements, and their inclination with respect to the failure surface, on the development of resisting tensile stresses in the reinforcements. A parametric study is conducted to evaluate the effect of these parameters on the behaviour of the reinforced soil material. An attempt is also made to verify the proposed model by comparing numerical predictions with available experimental results of both triaxial and direct shear tests on reinforced soil samples. This model can be used for analysis and design of reinforced soil walls with different types of tension resisting inclusions to predict tension forces under expected working loads.     

Numerical analysis of the response of reinforced soils to direct shearing: Part 2

A soil-reinforcement load transfer model was developed by the authors¹ to simulate the response of the reinforced soil material to triaxial compression and direct shearing. This paper presents the application of the proposed model for the numerical analysis of direct shear tests on sand samples reinforced with different types of tension resisting reinforcements. A parametric study is conducted to evaluate the effect of the mechanical characteristics and dilatancy properties of the soil, extensibility (elastic modulus) of the reinforcements, and their inclination with respect to the failure surface on the response of the reinforced soil material to direct shearing. An attempt is made to verify the proposed model by comparing numerical test simulations with experimental results reported by Jewell,² and Gray and Ohashi.³ Comparisons of predicted and experimental results illustrate that the model can provide adequate simulations of the response of the reinforced soil material to shearing. In particular, it allows an evaluation of the effect of soil dilatancy (or contractancy), and extensibility of the reinforcement on tension forces generated in inclusions during shearing.     

Granite rock fragmentation at percussive drilling – experimental and numerical investigation

The aim of this study is to numerically model the fracture system at percussive drilling. Because of the complex behavior of rock materials, a continuum approach is employed relying upon a plasticity model with yield surface locus as a quadratic function of the mean pressure in the principal stress space coupled with an anisotropic damage model. In particular, Bohus granite rock is investigated, and the material parameters are defined based on previous experiments. This includes different tests such as direct tension and compression, three‐point bending, and quasi‐oedometric tests to investigate the material behavior at both tension and confined compression stress states. The equation of motion is discretized using a finite element approach, and the explicit time integration method is employed. Edge‐on impact tests are performed, and the results are used to validate the numerical model. The percussive drilling problem is then modeled in 3D, and the bit‐rock interaction is considered using contact mechanics. The fracture mechanism in the rock and the bit penetration‐ resisting force response are realistically captured by the numerical model. Copyright © 2013 John Wiley & Sons, Ltd.     

强降雨条件下具有张裂缝边坡临界滑动场

降雨持时较长且雨强较大时,雨水的入渗不仅会增大孔隙水压力,且易使张裂缝充水形成静水压力,对边坡稳定不利。结合降雨条件下饱和-非饱和渗流分析,考虑降雨过程中的瞬态孔隙水压力场与瞬态强度场,并同时考虑张裂缝充水时的静水压力,对边坡临界滑动场法进行改进,提出降雨条件下具有张裂缝边坡临界滑动场数值模拟方法,且对其进行了验证。将该方法用于一个典型均质黏土边坡算例,结果表明,（1）文中方法可考虑张裂缝具体位置、深度及其充水状态下对边坡稳定性的影响,并能搜索出任意形状危险滑面,计算结果合理可靠;（2）降雨条件下张裂缝中静水压力对边坡稳定性及滑面形状有较大影响;（3）降雨条件下张裂缝位置对边坡稳定有较大影响,距坡肩越近,张裂缝对边坡稳定性影响越大;（4）张裂缝处在最不利位置且充水时存在一个最不利深度。     

含孔岩石加载过程的热辐射温度场变化特征

选择含圆孔岩石作为试样,利用单轴加载试验系统和红外热像仪,并结合数值模拟分析手段,对模型受力及破裂过程的热辐射时空演化特征进行了试验研究。结果表明：含圆孔岩石在加载过程中压、拉应力呈对称分布,导致热像的升温和降温区对称分布,即压应力区升温,拉应力区降温,应力场与红外辐射温度场之间呈很好的对应关系。试件加载后期发生破裂,破裂的性质与红外辐射温度变化有密切关系,剪性破裂辐射温度升高,而张性破裂辐射温度无明显变化。试件最终呈剪性破裂,而该区域的带状升温是岩石破裂失稳的重要红外前兆。试验结果再次揭示,利用热成像技术可以监测岩石的应力与灾变现象。     

An improved implicit numerical integration of a non‐associated,three‐invariant cap plasticity model with mixed isotropic–kinematic hardening for geomaterials

The Sandia GeoModel is a continuum elastoplastic constitutive model that captures many features of the mechanical response for geological materials over a wide range of porosities and strain rates. Among the specific features incorporated into the formulation are a smooth compression cap, isotropic/kinematic hardening, nonlinear pressure dependence, strength differential effect, and rate sensitivity. This study attempts to provide enhancements regarding computational tractability, domain of applicability, and robustness of the model. A new functional form is presented for the yield and plastic potential functions. This reformulation renders a more accurate, robust, and efficient model as it eliminates spurious solutions attributed to the original form. In addition, we achieve a high‐performance implementation, because the local iterative method is allowed to recast residual vectors with a uniform dimensionality. The model is also furnished with a smooth, elliptical tension cap to account for the tensile yielding. Moreover, an efficient algorithm is introduced, which decreases the computational cost by differentiating the updated shear yield surface from the cap surface based on the trial relative stress state. Finally, various numerical examples including a large‐scale boundary value problem are presented to demonstrate the fidelity of the modified model and to analyze its numerical performance. Copyright © 2015 John Wiley & Sons, Ltd.     

Experimental and numerical investigations on the performance of sacrificial piles in reducing local scour around pile groups

This paper presents both experimental study and numerical simulation to investigate the effects of an innovative sacrificial piles concept to reduce local scour surrounding service piles. These sacrificial piles are placed upstream of the pile groups, which serves to change the flow field and reduce the erosive force before it reaches the service piles. Four types of pile group arrangements were considered in the experimental and simulation analyses, including single pile, double pile in tandem arrangements, double pile with side-by-side arrangements, and pile array of 2 × 2 arrangements. Scaled experiments were conducted in flume for each group of design, and local scour around piles with and without sacrificial piles was monitored. It was found that the sacrificial piles effectively reduced the development of local scour surrounding the service piles, especially at their heels. The results also showed that the use of sacrificial piles reduced the time required for the scour depth to reach the equilibrium conditions than those without the sacrificial piles. To corroborate the experimental observation, numerical simulation with three-dimensional computational fluid dynamics model has also been carried out with commercial code with verified erosion model. The development of dynamic scour process and maximum scour depth were obtained by the computational simulations. Comparison of experimental and computational results both indicates that the use of sacrificial pile is effective in reducing the scour depth. The extent of reduction in the maximum scour depth varied with the arrangement of sacrificial piles.     

大型岩质滑坡地震变形破坏过程物理试验与数值模拟研究

以汶川地震触发的大光包岩质滑坡为例,结合野外现场调查,以其地质结果为背景,建立起破坏前的物理模型和三维数值模型,利用振动试验台和数值计算方法对滑坡变形破坏过程进行了研究。物理试验与数值模拟方法相互验证,取得了较为一致的结果。研究结果表明: 该滑坡的破坏模式为坡体顶部与中部拉张贯穿破坏中部沿层面滑移前缘剪切破坏,中部拉裂缝与主滑面首先形成滑动边界,前缘首先滑出; 滑坡变形过程中的加速度与速度响应研究表明其放大效应明显。同时,通过对比基岩与滑带加速度与速度放大系数,显示了结构面对斜坡变形破坏过程的控制作用。     

Numerical Study on Coalescence of Pre-Existing Flaw Pairs in Rock-Like Material

The present numerical study, which is an extension of our previous numerical analysis on cracking processes of a single pre-existing flaw, focuses on the coalescence of two pre-existing parallel open flaws in rock subjected to a uniaxial compressive loading. To facilitate a systematic investigation, the arrangements of the flaw pair are classified into 11 categories. Simulations engaging AUTODYN are conducted on each category. The numerical results are compared with some published physical experimental test results. Eleven typical coalescence patterns are obtained, which are in good agreement with the experimental results, which include two coalescence patterns obtained in flaw pair arrangements (II) and (VIII″) not being reported in previous studies. The information gathered in the simulations helps identify the type (tensile/shear) of each crack segment involved in the coalescence. Most of the coalescence cracks initiate at or around the flaw tips, except those in flaw pair arrangements (II) and (IX′) with a very short ligament length, in which the coalescence cracks initiate on the flaw surfaces away from the flaw tip regions. Based on the numerical simulation results, the properties of the 11 coalescence patterns are obtained. Except those in flaw pair arrangements (II) and (IX′), the other coalescence patterns can be interpreted with respect to the basic crack types—tensile wing crack, horsetail crack and anti-wing crack. In addition, based on the type of crack segments involved in coalescence, namely tensile and shear, the coalescence can be classified into T mode (tensile mode), S mode (shear mode) and TS mode (mixed tensile–shear mode).     

Three‐dimensional simulations of tensile cracks in geomaterials by coupling meshless and finite element method

Failure in geotechnical engineering is often related to tension‐induced cracking in geomaterials. In this paper, a coupled meshless method and FEM is developed to analyze the problem of three‐dimensional cracking. The radial point interpolation method (RPIM) is used to model cracks in the smeared crack framework with an isotropic damage model. The identification of the meshless region is based on the stress state computed by FEM, and the adaptive coupling of RPIM and FEM is achieved by a direct algorithm. Mesh‐bias dependency, which poses difficulties in FEM‐based cracking simulations, is circumvented by a crack tracking algorithm. The performance of our scheme is demonstrated by two numerical examples, that is, the four‐point bending test on concrete beam and the surface cracks caused by tunnel excavation. Copyright © 2014 John Wiley & Sons, Ltd.     

A boundary element solution to land subsidence above 3-D gas/oil reservoirs

A linear boundary element (BE) model is proposed for the uncoupied simulation of land subsidence due to gas, oil and hot water production over three-dimensional (3-D) arbitrarily shaped reservoirs. The pore pressure decline is assumed to be specified in advance, e.g. via a numerical model of flow. Use is made of the fundamental solution derived in 1885 by Boussinesq for a vertical load acting upon the traction-free surface of a semi-infinite medium. A straightforward application of Betti's (1872) reciprocal theorem allows for the development of a boundary integral whose numerical execution yields directly the downward settlement over the point of interest. The new procedure is applied to assess land sinking caused by an uniform pore pressure decline occurring within fields of elliptical shape and to explore the influence of the assumption of small reservoir thickness which underlies the ‘tension center’ or ‘strain nucleus’ approach previously developed by Geertsma in 1966. The results emphasize the numerical efficiency of the solution and the promising features of the BE method for the evaluation of ground subsidence in 3-D problems. The present model is based on the theory of the linear poroelasticity and is implemented for a mechanically homogeneous and isotropic half-space. It allows for any arbitrary geometry of the reservoir and for a non-uniform distribution of the pore pressure decline. It may easily be extended to other physical settings for which a vertical surface point load solution is available.     

Predicting the surface tension of aqueous 1:1 electrolyte solutions at high salinity

The surface tension of the air/water interface is a phenomenon of particular interest in the water-unsaturated zone of porous media because it influences the contact angle and consequently the capillary water volume. A mechanistic model based on the modified Poisson-Boltzmann equation and the Pitzer theory is described and used to predict, under isothermal and isobaric conditions, the surface tension of 1:1 electrolytes at high salinity. These theories enable the determination of the electrical potential at the air/water interface and the activity coefficient of the ionic species in the bulk pore water, respectively. Hydration free energies of the structure-making and structure-breaking ions that influence the surface tension at high salinity are taken into account. Structure-making ions flee the air/water surface because they can better organize the water dipoles in bulk water than at the interface. Structure-breaking ions are positively adsorbed at the air/water interface because the bulk water can better organize their hydrogen-bonding network without these ions. The resulting surface tension increases and decreases, respectively, compared to the surface tension of pure water. The predictions are in good agreement with the surface tension data of 1:1 electrolytes (NaCl, KCl, HCl, NaNO₃, KNO₃, HNO₃ aqueous solutions) and the optimized parameters depend on the effective electrostatic diameters of cations and on the hydration free energies of the ions at the interface.     

刚/柔性组合墙面加筋土挡墙内部破坏机制

基于离心模型试验成果,建立软土地基刚/柔性组合墙面加筋土挡墙离散-连续耦合数值模型,采用离散单元颗粒流程序PFC和有限差分程序FLAC分别模拟加筋土挡墙和软土地基,分析挡墙的变形、筋材拉力分布及内部破坏演化过程,并与刚性地基上挡墙情况进行比较,以探讨刚/柔性组合墙面加筋土挡墙的内部破坏机制。研究结果表明,数值模型计算结果与离心模型试验结果吻合,其内部潜在破坏面经过连接件锚固端位置,各层筋材拉力均在连接件锚固端位置最大;软土地基上挡墙和刚性地基上挡墙的内部破坏面均一致,且筋材均由下而上依次断裂;软土地基上挡墙内部破坏面与地基圆弧滑移面贯通,为复合滑动模式,而刚性地基上挡墙整体破坏模式为折线型,连接件埋深范围内加筋体沿其底部连接件水平滑出。     

A discussion on analytical and numerical modelling of the land subsidence induced by coal seam gas extraction

Coal seam gas (CSG) is an increasingly important source of natural gas all over the world. Although the influence of conventional oil and gas extraction on surface subsidence has been widely recognized and studied, few studies are carried out on the surface subsidence in coal seam gas fields and its impact on surface infrastructure and the environment. This paper discusses modelling of the surface subsidence associated with coal seam gas production by applying both analytical and numerical methods. By comparison of results from the numerical model and two analytical models, i.e. the disc-shaped reservoir model and the uniaxial compaction model, the analytical solutions cannot describe the complex process of water and gas extraction and have the limitations to predict the surface subsidence, while the numerical model can be better used in prediction of subsidence. After applying the numerical model in numerical analysis, the deformation characteristics of coupled fluid flow, and the effects of permeability change of coal seam, associated overlying and underlying layers, and depressurization rates on surface subsidence are investigated. The results demonstrate that the proposed model can simulate the production of water and gas from coal seams and the associated surface subsidence.     

基于离心机和数值模拟的软硬互层反倾层状岩质边坡变形特征分析

我国西部山区建设揭露了众多大型弯曲倾倒变形体,它们多具有软硬互层结构。为进一步探明软硬互层反倾岩质边坡的变形破坏规律,本研究融合大型土工离心机试验与自开发的可以综合考虑节理面拉剪和压剪破坏的Hoek-Brown与Mohr-Coulomb联合强度准则,对此类边坡进行试验与数值模拟分析。首先,结合监测点位移和应力曲线对边坡的变形破坏过程进行详述并验证了所提出的强度准则及所建立的数值模型的正确性;然后,基于此数值模型研究不同几何因素对此类边坡倾倒破坏特征的影响。结果表明:（1）节理单元采用Hoek-Brown与Mohr-Coulomb联合强度准则可以较准确地模拟软硬互层反倾边坡的层间错动以及岩层弯折;（2）此类坡体倾倒变形破坏全过程为:层间先出现相互错动,然后边坡自坡脚部位开始出现弯曲变形,随后坡体后缘出现拉张裂缝,与此同时边坡整体向临空面弯曲倾倒,最终形成2个或3个破坏面;（3）随着岩层倾角的增大,边坡一级破坏面逐渐向坡体深处发展;（4）随着硬/软岩层厚比的减小,坡顶竖向位移变小,且坡体滑动的整体性逐渐增强;（5）随着软/硬岩层厚比的增加,坡体破坏面逐渐由粗糙的“锯齿状”向平滑的“圆弧状”过渡。     

基于格子Boltzmann方法的饱和冻土孔隙成冰介观尺度模拟

路基冻胀问题严重影响寒区高速铁路的安全服役,而成冰相变过程是解释冻胀机制的关键。基于介观尺度的格子Boltzmann方法,将修正的孔隙水冻结温度算法与焓法固液相变格子Boltzmann模型相结合,模拟了悬浮液滴冻结和冻土孔隙水成冰两个过程,分别揭示了液态水在自由状态和孔隙束缚状态下冰水相变的细观机制。计算结果表明：土体孔隙中冰晶由中心向外生长的过程与悬浮在空气中的液滴冻结过程截然不同,并且孔隙水越接近颗粒表面,其冻结温度越低。相同粒径颗粒按照不同排列方式得到的冻结特征曲线（soil freezing characteristic curves,简称SFCC）具有明显差异;不同粒径的SFCC随着颗粒增大残余水含量逐渐变少,形态更加陡峭。通过与文献试验结果对比,验证了格子Boltzmann方法的有效性,表明该方法能够为研究多孔介质水气迁移与相变过程提供介观尺度的新手段。     

Application of the DDM and fracture mechanics model on the simulation of rock breakage by mechanical tools

A DDM (displacement discontinuity method) program coupled with a modified energy criterion is used to simulate the development of cracks and chips by indentation tools. In our analysis a cavity model is applied to represent the expansion of crushed rock to the surrounding rock and the cracks are formedin two-dimensional and quasi-static conditions. The model parameters, rock properties and load magnitudes are varied in the numerical calculations. The results show that chips are formed by multiple mechanisms of either tension or shear, or their combinations. The cracks may either propagate to the free surface to form chips or rest in the rock subsurface. The crack development is dependent upon rock and fracture properties, loading force and tool characteristics. The DDM is a convenient tool in the study of rock fragmentation and cracks.     

陡倾断层上下盘开挖引起地表变形的数值模拟分析

考虑地下开挖区与断层相对位置关系,采用简化数值模型对陡倾断层上下盘开挖引起的地表变形特征进行了数值模拟。结果表明,当断层位于地下开挖引起地表变形的压缩区且开挖区位于下盘时,地表出现不连续变形的可能性小;当断层位于地下开挖引起地表变形的压缩区且开挖区位于上盘时,随下向开挖进行地表断层处裂缝有减小趋势;而当断层位于地下开挖引起的地表变形拉张区时,断层出露处水平拉张明显,表现为拉张裂缝,随下向开挖深度增加、规模增大,不论开挖区位于上盘还是下盘,极可能导致地表出现正断层式的错动．