单轴压缩下断续节理岩体动态损伤本构模型

断续节理将对工程岩体的强度及变形等力学特性产生显著影响,损伤力学中视节理为岩体的一种宏观损伤,因而采用损伤张量来刻画其对岩体的影响。目前学术界提出了用节理的几何、强度及变形等3类参数来描述节理的物理力学性质,而目前的岩体损伤张量计算方法都只涉及前2类参数,均没有涉及其变形参数即法向及切向刚度。为此,在前人研究的基础上,基于断裂及损伤理论提出了考虑节理法向及切向刚度的单轴压缩下单条断续节理引起的损伤张量计算公式,进而通过考虑节理间相互作用给出了单组单排或多排节理岩体损伤张量计算公式。其次,以岩石细观动态损伤模型为基础,结合宏细观损伤耦合观点提出了一个能够同时考虑节理几何、强度及变形参数的断续节理岩体动态损伤本构模型。最后,利用该模型讨论了节理参数及载荷应变率等对岩体动态力学特性的影响,认为节理长度减小及摩擦角增大将导致岩体动态峰值强度及弹性模量增大;岩体动态峰值强度及弹性模量则随着节理法向及切向刚度的增大分别减小或增大;而当节理法向及切向刚度按照同一比例增大时,岩体动态峰值强度及弹性模量则是增大的。岩体动态峰值强度与载荷应变率呈正相关。     

综合考虑宏细观缺陷的岩体动态损伤本构模型

针对节理岩体同时含有节理、裂隙等宏观缺陷及微裂隙、微孔洞等细观缺陷的客观事实, 提出了在节理岩体动态损伤本构模型中应同时考虑宏细观缺陷的观点。为此, 首先对基于细观动态断裂机理的经典岩石动态损伤本构模型—TCK(Taylor-Chen-Kuszmaul)模型进行了阐述, 其次基于Lemaitre等效应变假设推导了综合考虑宏细观缺陷的复合损伤变量(张量), 进而在此基础上建立了相应的节理岩体动态损伤本构模型, 并利用该模型讨论了载荷应变率及节理条数对岩体动态力学特性的影响规律。结果表明, 在不同载荷应变率下试件在变形初始阶段是重合的, 而后随着应变的增加, 试件峰值强度、峰值应变及总应变均随载荷应变率的增加而增加; 随着节理条数的增加, 试件峰值强度逐渐降低, 但降低趋势逐渐变缓并趋于某一定值。上述研究结论与目前的理论及实验研究结果的基本规律是一致的, 说明了本模型的合理性。     

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

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

An empirical statistical constitutive relationship for rock joint shearing considering scale effect

The scale effect of rock joint shearing is of great significance in rock engineering. Most existing shear constitutive models could describe the pre- and post-peak deformation of rock joints, but only in one particular scale, that is, those existing models will fail to depict the rock joint shearing in different length scales. Therefore, this study aims to establish a constitutive relationship for rock joints with considering the scale effect. Based on the assumption of a random statistical distribution of rock material strength and statistical mesoscopic damage theory, damage variables are defined as the ratio of the number of damaged elements to the total number in the shear process. Together with the nonlinear relationship between the microelement failure and the joint scale, an empirical statistical constitutive relationship for joint is established. And then, the determination method of constitutive relationship parameters and the variation laws with the scale are discussed. Results show that the predicted results of the proposed empirical relationship not only agree well with the experimental results but also fully describe nonlinear deformation, pre-peak softening, post-peak softening, residual stage, and other mechanical properties of the shear deformation of joint with different dimensions, thereby demonstrating the rationality of the constitutive relationship. The physical meaning of the constitutive relationship parameters is clear, and the expressions of the constitutive relationship parameters can be deduced from the experimental results. In addition, the influence of scale effect on the shear deformation of rock joints can be quantified using parameters, which help accurately describe the action form of scale effect.     

甘肃北山地质处置库围岩节理抗剪强度经验估算

节理抗剪强度参数是地质处置库预选地段工程地质对比和围岩稳定性分析的重要指标。本文在旧井地段英云闪长岩节理分组的基础上,运用定向统计测量方法估测节理粗糙度系数,通过评价节理粗糙度系数的尺寸效应,确定节理抗剪强度经验估算有效长度,由JRC-JCS模型求得各组节理4个方向的抗剪强度参数,并评价了地质处置库围岩节理抗剪强度的各向异性特征。     

节理抗剪强度综合评价的试验研究

选取不同吻合度自然节理试样,进行干燥状态和饱和状态的JRC-JCS模型和JRC-JMS模型试验研究。结果表明,在节理粗糙度系数的定向统计测量和尺寸效应分析,以及剪切过程衰减折减的基础上,运用JRC-JCS模型综合评价节理抗剪强度可靠性较好,考虑吻合度系数的JRC-JMC模型比JRC-JCS模型具更好的预测能力。     

含节理岩体爆生裂纹扩展的动焦散模型实验研究

应用动态焦散线测试系统,模拟含节理岩体断裂爆破过程,进行了PMMA模型透射式动态焦散线实验,着重研究了爆炸初始裂纹与节理面不同夹角的情况下,裂纹尖端动态应力强度因子的变化规律,裂纹穿过节理面的扩展规律,以及裂纹扩展速度的变化规律。实验结果分析表明,爆生裂纹穿过节理面时,裂纹尖端的动态应力强度因子和裂纹扩展速度显著下降,穿过节理面后,应力强度因子又有所增强;裂纹穿过节理面时,裂纹会沿节理面偏离一段距离后沿初始裂纹方向继续扩展。研究结果可以为节理岩体的断裂爆破提供理论依据。     

JRC-JCS模型抗剪强度估算的平均斜率法

结构面对工程岩体的稳定性有着决定性的影响,结构面的抗剪强度参数是分析岩体稳定的关键的参数。本文回顾了结构面抗剪强度估算的JRC-JCS模型及其估算方法,提出了可快速估算结构面抗剪强度参数的平均斜率法,并通过对比分析阐明平均斜率法的优越性。     

Modeling 3D revolute joint with clearance and contact stiffness

The clearances in the kinematic joints are due to deformations, wear, and manufacturing errors; the accurate modeling of these effects in multibody analysis is a complex issue but in many practical applications, it is mandatory to take into them into account in order to understand the actual behavior of mechanical systems. In this paper, the authors present a general computer-aided model of a 3D revolute joint with clearance suitable for implementation in multibody dynamic solvers. While a perfect revolute joint imposes kinematic constraints, the proposed revolute joint with clearance leads to a force constraint. The revolute joint has been modeled by introducing a nonlinear equivalent force system, which takes into account the contact elastic deformations. The model depends on the structural and geometrical properties of materials in contact that have been investigated using finite element models. The purpose is to give a general approach to study the influence of actual joints on kinematic, dynamic, and structural behavior of mechanisms. The proposed model has been applied in dynamic simulations of a spatial slider-crank mechanism.     

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.     

DEM analysis of the effect of joint geometry on the shear behavior of rocks

In order to comprehensively investigate the effect of different joint geometries on the shear behavior of rocks, the Distinct Element Method (DEM) was utilized with a new bond contact model. A series of direct shear tests on coplanar and non-coplanar jointed rocks was simulated using the PFC2D software, which incorporates our bond contact model. Both coplanar jointed rocks with different joint persistence and non-coplanar ones with different joint inclinations were simulated and investigated numerically. The numerical results were compared and discussed with relevant laboratory tests as well as some reported numerical works. The results show that for coplanar jointed rocks, the peak shear stress decreases nonlinearly with the joint persistence, and the failure process can be divided into four stages: elastic shearing phase, crack propagation, failure of rock bridges, and residual phase. For non-coplanar jointed rocks, as the absolute value of the inclination angle of the rock joints increases, its shear strength increases, changing the failure patterns and the length of new fractures between existing cracks. When the absolute value increases from 15° to 30°, the average shear capacity increases the most as 39%, while the shear capacity increases the least as 2.9% when the absolute value changes from 45° to 60°. There is a good consistency of the failure patterns obtained from experiments and numerical tests. All these demonstrate that the DEM can be further applied to rock mechanics and practical rock engineering with confidence in the future.     

Effect of non-persistent joints distribution on shear behavior

The particle flow code 2D (PFC2D) is used to establish a coplanar, non-persistent joint model. Three joint distribution types, namely, both-side (type a), scattered (type b), and central (type c), are set according to their position. Numerical simulations of the direct shear test are conducted to investigate the effect of non-persistent joint distribution and connectivity on shear mechanical behavior. Simulation results are in good agreement with the analytical solutions to Jennings' criterion, and show: (1) type-c and type-b joints have high strength, whereas type-a joints have low strength. Shear strength and modulus increase with a decrease in joint persistency, and the shear displacement that correspond to shear strength increases with a decrease in persistency. (2) The brittle failure characteristics of the sample are evident when the intact rock bridge area is large. Reinforcement at both ends of the joint limits shear deformation, and shear strength can be effectively improved when joint persistency is large. The small-area dispersed reinforcement joint method cannot effectively improve shear strength. (3) The comprehensive shear strength parameters and the shear strength of the non-persistent joints can be predicted well using Jennings' criterion. Cohesion is the dominant factor that controls shear strength.     

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.     

On flaw size and distribution in lap joints

By using adhesive as the bonding substance between metals or polymeric materials, simple structural joints can be made to bear relatively high loads. Applications have increasingly been made in substituting adhesive joints for conventional mechanical fastenings, especially in the aircraft and aerospace industries where weight is a predominant factor. In order to design a most effective adhesive-bonded joint, an understanding of the stress distribution along the joint is as important as the physical properties of the bonding agent. One of the most common and widely used adhesive joints is the single lap joint.Recent investigations using various analytical models have revealed that the cause of failure in an idealized ‘defect free’ lap joint is primarily due to the localized effect of high stress concentration at the lap ends. With the presence of flaw like defects in the adhesive layer, the load transfer from adherend to adhesive is expected to be different from the idealized joint. In addition, localized stress concentrations induced by irregular adhesive defects that may be found in practical engineering applications can further reduce fracture strength of such an imperfect joint.This paper is intended to describe an investigation into the effect of internal adhesive flaw size and distribution on the fracture behaviour of adhesive-bonded lap joints. The finite element method is used to gain a quantitative understanding of the localized shear stress distributions due to the presence of the internal flaws along the bonding layer. It is observed that the reduction in the fracture strength is relatively small when a flaw is located in the central portion of the bonding length. However, a flaw located near the lap ends of the adhesive joint can cause marked reduction in the fracture strength, due to its interaction with the high stress concentration at the lap ends.     

A method to simplify the strength design of bolted joints—Case of connecting-rod bolts

This study shows the development of the strength-design system for general bolted joints. Specifically, a simplified method of the strength design for the connectingrod bolts is proposed. The connecting-rod joint is transformed into a single-bolted-joint model by using the ring model. The axial stress on the bolt and separation behavior at the contact plane are discussed. The strength design of the connecting-rod joint can be performed by applying the unified procedure to the single-bolted-joint model. The local deformation of the connecting-rod joint and the stress on the bolt were measured and the analytical results are compared with the experimental results. Both results show fairly good agreement.     

STRENGTH PROPERTIES OF JOINTED ROCK MASSES BASED ON THE HOMOGENIZATION METHOD

This paper aims to determine the strength properties of jointed rock masses by means of the homogenization method.To reflect the microstructure of jointed rock masses,a representative element volume (R...     

不同节理角度和地应力条件下岩石双孔爆破的数值模拟

在采矿工程、地下交通工程和水利水电工程等岩体开挖工程中,爆破法仍是一种主要的破岩方法. 实际天然岩体中存在的裂隙、节理等不连续面和所处的地应力场环境,都会对爆破荷载传播过程和最终的破碎效果产生重要影响. 本文把岩石爆破视为爆炸应力波动态作用和爆生气体压力准静态作用两个独立的先后作用过程,同时考虑初始应力场的静态作用,建立岩石爆破的力学模型. 基于渗流方程描述爆生气体在爆生裂纹中的传播过程,进而基于流固耦合理论实现爆生气体准静态压力对裂纹尖端的力学作用. 爆炸应力波主要在孔周引起压碎区和径向微裂纹区,随后爆生气体压力楔入径向裂纹尖端,促使裂纹进一步扩展,形成径向主裂纹. 数值模型可以再现孔周压碎区、径向微裂纹区、径向主裂纹区萌生、扩展的完整演化过程. 针对不同节理角度和地应力条件下岩石双孔爆破过程的数值模拟结果表明,初始地应力场的压应力作用不利于爆生裂纹的萌生与扩展,但节理的存在对裂纹的扩展具有明显的导向和促进作用,有利于爆生裂纹沿节理面方向的扩展.      

基于时域自适应算法的单向粘弹性节理岩体的等效分析

在粘弹性节理岩体的数值仿真中,必需合理地计及节理的影响.若将节理与岩体作为独立的材料组分考虑,在节理密集的情况下,时空两方面的计算开销可能是很难承受的.本文通过时域展开技术与一个简单的等效假定,得到了递推格式的粘弹性单向节理岩体的等效本构关系,及等效场的递推求解格式,并由此数值模拟了一个粘弹性节理岩体中的洞室,计算结果与考虑独立组分的ANSYS粘弹性分析进行了比较,二者的计算平均误差约为11.03%,但前者的计算效率为后者16.81倍.本文工作有可能为粘弹性节理岩体的数值模拟提供一条新的途径.     

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

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

Shear relaxation characteristics of rock joints under stepwise loadings

The stress relaxation characteristic of rock mass is an important aspect of rheology and has important practical significance for rock engineering. In order to investigate the relaxation characteristic of rock joints with different slope ratios and normal stresses, a series of shear stress relaxation tests were conducted on artifical rock joints poured by cement mortar. Test results show that the relaxation curves can be divided into three stages, i.e. instantaneous relaxation stage, attenuation relaxation stage, and stable relaxation stage. Furthermore, the nonlinear Maxwell relaxation equation was obtained by using the relation between the viscosity coefficient and time, and the theoretical curves based on the empirical equation agreed well with the test results. Moreover, the change law of the initial viscosity coefficient was investigated. Accordingly, a stress relaxation method, termed as relaxation stress peak method, was proposed to determine the long-term strength of rock joints.