经历剪切变形历史的岩石节理表面粗糙特性分析

基于两类岩石节理剪切试验，采用非接触式激光测距仪，实测了剪切前后节理表面粗糙数据，统计分析了岩石节理表面粗糙特性的三维粗糙参数的物理含义，并将这些参数与相应的二维参数进行了充分比较。结果表明，三维参数与二维参数间存在的近似线性关系，而多个三维参数相对比单一的二维参数更全面有效地描述了节理表面粗糙特性。     

一种新的岩石节理面三维粗糙度分形描述方法

研究并提出一种新的岩石节理面三维粗糙度分形描述方法。首先,基于激光扫描数据将节理表面离散成三角网,并建立与剪切方向相关的三维均方根抵抗角的计算方法。其次,运用分形数学理论,提出一种新的基于三维均方根抵抗角的节理面粗糙度分形描述方法。最后,采用新方法对天然玄武岩节理和花岗岩张拉型节理的粗糙特性进行分析。研究结果表明,提出的新方法能够较全面地反映节理面的三维几何形貌信息,并能描述节理粗糙度的各向异性特性。研究成果为进一步建立岩石节理面的三维剪切强度公式和剪切本构理论奠定了基础。     

一种新型粗糙节理面随机强度模型及其应用

天然岩体在长期地质作用下会生成各种节理裂隙等不连续面,而地下工程结构的稳定性一般取决于这些不连续面的强度。在众多因素中,表面形态对岩石节理面剪切强度具有决定性影响。为了系统研究岩石节理面剪切强度的确定方法,把岩石节理面概化为一系列高度不同的微长方体凸起组成的粗糙表面结构,且微长方体凸起有剪胀破坏和非剪胀破坏两种模式。综合微长方体凸起破坏规律,应用概率密度函数描述节理面表面起伏分布的影响,建立了粗糙节理面随机强度模型,推导了节理面剪切强度理论公式,提出了节理面强度的随机评价方法。基于随机强度模型和评价方法编制Matlab计算程序计算自然粗糙节理面的剪切强度,并将计算结果与试验结果进行比较分析。研究表明:粗糙节理面随机强度模型综合了粗糙节理面表面形态和法向应力对节理剪切强度的影响机制,理论计算值与试验数据吻合良好,可以较好的评价粗糙节理的峰值剪切强度和残余剪切强度。该随机模型可作为进一步深入研究的重要基础,分析结构面的连续剪切过程,建立更完善的节理面强度模型。     

不同节理粗糙度类岩石材料三轴压缩力学特性试验研究

粗糙度是影响节理岩体强度与变形特性的重要因素之一。首先使用3D 打印机制作模具,并浇筑形成不同粗糙度（节理粗糙度系数JRC=2、7、12、17、22）的节理岩石试样。采用GCTS高温高压动静岩石三轴试验系统,对含有不同粗糙度节理岩石试样进行了三轴压缩试验,获得了不同粗糙度节理岩石试样的三轴应力–应变曲线,分析了JRC对岩石三轴强度和变形特性的影响规律,在三轴加载过程中采用声发射测试系统,分析了不同粗糙度节理岩石试样的声发射特性。运用数字三维视频显微系统观察节理面形态,讨论了不同围压下节理岩石试样峰值强度与JRC之间的关系。研究结果表明,节理面的存在直接导致节理岩石试样强度的大幅度降低,JRC对岩石破坏裂纹的形态、数量和空间分布特征亦有很大的影响,随着JRC值的增大,岩石节理面的抗剪强度增大,岩石试样的三轴抗压强度也会增大,岩石试样由脆性破坏转变为延性破坏。     

基于3D打印的含复杂节理岩石直剪特性及破坏机制研究

岩体中节理的几何形态及力学特性是影响其剪切力学特性及破坏模式的重要因素之一。基于3D打印技术,建立了不同节理粗糙系数(JRC)的节理模型、几何形态节理模型和复杂裂隙网络物理模型,通过开展室内直接剪切试验分析了各组试件的剪切强度及破坏模式。结果表明节理模型的抗剪强度随JRC波动性较大,波动幅值越高,峰值剪切位移越低;平面形节理模型的峰值抗剪强度最低,矩形节理模型的峰值抗剪强度最高,正弦形和三角形节理试件的抗剪能力相近;离散裂隙网络模型和粗糙裂隙网络模型的峰值抗剪强度显著低于实心试件,考虑了节理粗糙性的裂隙网络模型抗剪强度高于直线型节理模型;实心试件破坏模式为典型脆性剪切破坏,裂隙网络模型的破坏模式相对复杂,沿着剪切方向主剪切裂面波动萌生,破断面由多个节理面的交叉点破坏与沿节理面的滑移构成。研究成果可以为3D打印技术的推广和复杂节理岩体剪切力学特性的室内试验研究提供参考。     

基于孔内影像的岩石节理粗糙度特征研究

目前对岩石结构表面粗糙度的研究往往只局限于地表,难以反映深部岩石节理粗糙度的特征。钻孔孔壁上节理轮廓线包含有三维信息的特点,本文开展了基于数字钻孔摄像技术的岩石节理粗糙度分形特征的研究,利用数字钻孔摄像系统获取地下深度岩石节理全景图,采用边缘检测技术从全景图中提取出节理轮廓线,对其进行空间变换和视距离变换得到地下岩石节理粗糙表面轮廓线的真实状况。与Barton提出的10条标准剖面曲线进行对比得到每条轮廓线的JRC值,并计算其分形维数,根据最小二乘法原理拟合出JRC与分形维数之间的关系为:(D)=JRC=-541.9x2+1362x-818.53。本文研究内容为描述地下深部天然节理的结构及其特征提供了基础,对更深入地研究地下深部岩石节理的表面空间状况有重要的意义。     

于肃北山花岗岩节理表面形态特性研究

花岗岩是中国高放废物处置的候选围岩,大量存在的节理面是核素迁移的主要通道,同时节理面对岩体的力学、水力性质起着支配作用。而节理面的表面形态与这些性质密切相关,因此需要对候选花岗岩体节理面展开多样本、大尺寸的现场表面形态研究。以中国高放废物深地质处置甘肃北山预选区某采石场中的天然节理面为研究对象,运用基于几何理论的测量方法,获得了大量不同尺度、不同方位的现场节理面的表面形态特征曲线。采用Barton直边法和统计参数法处理这些特征曲线,系统分析了特征曲线的尺度效应、各向异性特征以及与不同节理类型的关系。为进一步分析裂隙岩体的力学、水力特性建立了基础。     

甘肃北山花岗岩节理表面形态特性研究

花岗岩是中国高放废物处置的候选围岩,大量存在的节理面是核素迁移的主要通道,同时节理面对岩体的力学、水力性质起着支配作用。而节理面的表面形态与这些性质密切相关,因此需要对候选花岗岩体节理面展开多样本、大尺寸的现场表面形态研究。以中国高放废物深地质处置甘肃北山预选区某采石场中的天然节理面为研究对象,运用基于几何理论的测量方法,获得了大量不同尺度、不同方位的现场节理面的表面形态特征曲线。采用Barton直边法和统计参数法处理这些特征曲线,系统分析了特征曲线的尺度效应、各向异性特征以及与不同节理类型的关系。为进一步分析裂隙岩体的力学、水力特性建立了基础。     

三维交叉裂隙渗流传质特性数值模拟

对交叉裂隙渗流传质特性的定量描述是研究整个裂隙网络渗透传质特性的基础。为真实模拟水流及溶质在三维交叉裂隙中的运移过程,首先通过三维轮廓仪获取天然岩石裂隙表面的形貌数据,再应用三维重构技术生成相应的三维交叉裂隙模型,随后求解Navier-Stokes方程,假定溶质运移满足Fick定律,模拟水流和溶质在三维交叉裂隙中的运移过程。通过对比粗糙裂隙模型与平行平板模型的模拟结果发现：粗糙度对流体的分布及流动状态存在显著的影响;不同进、出口工况下的流体流动及溶质运移状态亦表明：裂隙交叉的几何形貌会显著地影响溶质混合行为。这些结果表明,目前被广泛采用的平行平板模型在评估岩体内特别是交叉口的物质运移特性时将导致较大的偏差,在将来的研究中有必要针对裂隙交叉口的几何特征建立修正的模型以提高评估的准确性。     

岩石裂隙的非饱和渗透特性及其演化规律研究

岩石裂隙的非饱和渗透特性是岩土、能源和环境等领域科学研究中的热点问题。采用三维激光扫描获取花岗岩裂隙的表面形貌特征,分析裂隙微观形貌特征对非饱和渗透特性的影响。研究在张拉、压缩、剪切等复杂荷载作用下裂隙开度分布的演化规律,建立复杂荷载作用下岩石裂隙非饱和毛细压力曲线演化模型。基于裂隙的微观形貌特征推导了岩石裂隙非饱和相对渗透系数模型,通过与试验数据对比,验证了模型的准确性和有效性,并在此基础上建立了复杂荷载作用下岩石裂隙非饱和相对渗透系数演化模型。研究成果对非饱和条件下裂隙岩体的水-力耦合机制研究具有一定指导意义。     

Friction Factor of Water Flow Through Rough Rock Fractures

Fluid flow through rock joints occurs in many rock engineering applications. As the fluid flows through rough-walled rock fractures, pressure head loss occurs due to friction drag of the wall and local aperture changes. In this study, the friction factor was experimentally investigated by performing flow tests through sandstone fractures with joint roughness coefficient ranging from 5.5 to 15.4 under changing normal stress from 0.5 to 3.5 MPa. According to the experimental results, the friction factor was formulated as a function of two-independent variables—Reynolds number and relative roughness. Relative roughness is defined as the ratio of maximum asperity height to equivalent hydraulic aperture. The experimental results show that the proposed predictor of the friction factor fits the data with a coefficient of determination R ² > 0.93. Sensitivity analyses indicate that in general, the proposed friction factor increases with the relative roughness of confined fractures. The large difference of friction factor induced by relative roughness occurs when the Reynolds number is lower than unity, especially for Re < 0.2. For Reynolds numbers greater than unity, the difference of friction factor induced by relative roughness is smaller. Inclusion of joint roughness in calculating the friction resistance to fluid flow in rough rock joints and the influence of normal stress to the joints is a major step towards more accurate predictions for fluid flow in underground joint networks. This study provides a significant improvement in fundamental understanding of fluid flow in the jointed strata.     

基于Plesha本构的岩石节理多层结构模型研究

在Dong等工作的基础上[1],研究了基于Plesha本构的岩石节理多层结构模型。模型将节理面分解为不同层次细观结构面,每层结构面代表一级粗糙度,力学响应发生在最底层（基本面）,上层结构面受力性状由其下层结构面平均化得到,依次进行得到节理宏观力学性状。节理破坏由下而上分层次进行,基本面破坏后,其上层结构面转化为基本面,基本面采用Plesha节理模型。模型能模拟弹性变形、滑动变形、磨损、剪断、压碎、分离等作用机理。定义的接触因子与绕流因子能考虑节理接触面积及粗糙度变化对渗流的影响。通过将粗糙度定义为等效起伏角,可以模拟节理循环剪切性状。采用ABAQUS的用户子程序UEL进行了模型验证与参数分析。     

Some recent advances in the modelling of soft rock joints in direct shear

This paper presents a review of recent developments made by the authors into the modelling of rock joints in direct shear. Careful observation of laboratory direct shear testing on concrete/rock joints containing two-dimensional roughness has allowed theoretical models of behaviour to be developed. The processes modelled include asperity sliding, asperity shearing, post-peak behaviour, asperity deformation and distribution of stresses on the joint interface. Model predictions compare extremely well with laboratory test results. These models were then applied to direct shear tests on rock/rock joints, and although behaviour in general was well predicted, the strength of rock/rock joints was over-predicted. Direct shear tests have also been carried out on samples containing both two- and three-dimensional roughness to test the accuracy of the two-dimensional approximation to roughness adopted in the theoretical models.     

Numerical Investigations of the Dynamic Shear Behavior of Rough Rock Joints

The dynamic shear behavior of rock joints is significant to both rock engineering and earthquake dynamics. With the discrete element method (DEM), the dynamic direct-shear tests on the rough rock joints with 3D (sinusoidal or random) surface morphologies are simulated and discussed. Evolution of the friction coefficient with the slip displacement shows that the 3D DEM joint model can accurately reproduce the initial strengthening, slip-weakening, and steady-sliding responses of real rock joints. Energy analyses show that the strengthening and weakening behavior of the rock joint are mainly attributed to the rapid accumulation and release of the elastic energy in the joint. Then, effects of the surface roughness and the normal stress on the friction coefficient and the micro shear deformation mechanisms, mainly volume change and asperity damage, of the rock joint are investigated. The results show that the peak friction coefficient increases logarithmically with the increasing surface roughness, but decreases exponentially with the increasing normal stress. In addition, the rougher rock joint exhibits both higher joint dilation and asperity degradation. However, high normal stress constrains the joint dilation, but promotes the degree of asperity degradation significantly. Lastly, the effects of the 3D surface morphology on the shear behavior of the rock joint are investigated with a directional roughness parameter. It is observed that the anisotropy of the surface roughness consequently results in the variation of the peak friction coefficient of the joint corresponding to different shearing directions as well as the micro shear deformation mechanisms, e.g., the extent of joint dilation.     

Underestimation of roughness in rough rock joints

Numerous studies have been made to improve Barton's shear strength model for the quantification of rock joints. However, in these previous studies, the roughness and shear strength of the rock joint have been underestimated especially for relatively high undulated profiles (joint roughness coefficient (JRC) >14). The main factors of roughness underestimation in rough rock joints are investigated for the proper quantification of rock joint roughness. The aliasing effect and the roughness characteristics are analyzed by using artificial joint profiles and natural rock joint profiles. A 3D camera scanner is adopted to verify the main source of underestimation when using conventional measurement methods. Shear strength tests are carried out by using two types of shear apparatus to study the roughness mobilization characteristics, which may also affect the roughness underestimation. The results of joint roughness assessment, such as aliasing and undulation of waviness, show that the roughness can be underestimated in relatively rough joint profiles (JRC>14). At least two components of roughness parameters are needed to properly represent the joint roughness, for example, the amplitude and the inclination angle of joint asperity. Roughness mobilization is affected by both the normal stress and the asperity scale. Copyright © 2008 John Wiley & Sons, Ltd.     

循环剪切荷载作用下含二阶起伏体模拟岩石节理力学特性研究

利用人工材料浇注含二阶起伏体的模拟岩石节理试样,进行常法向荷载循环剪切试验,研究节理剪切力学特性在循环剪切过程中的劣化规律。试验结果表明：二阶起伏体对节理循环剪切力学特性有重要影响,剪切强度、剪切刚度、剪胀角随剪切循环次数增大而衰减,衰减趋势随着二阶起伏度的增大而加快;法向应力、二阶起伏度较大时,二阶起伏体对剪切力学特性的影响主要体现在第1轮剪切循环中,在随后的剪切循环中影响不明显;法向应力、二阶起伏度较小时,二阶起伏体的影响在前几轮循环剪切过程中均有较清晰的体现。基于Hertz接触力学理论,提出了节理面微凸体球面接触细观模型,揭示了节理循环剪切宏观试验现象的力学机制     

Effects of Cyclic Loading on the Shear Behaviour of Infilled Rock Joints Under Constant Normal Stiffness Conditions

The variation of the shear strength of infilled rock joints under cyclic loading and constant normal stiffness conditions is studied. To simulate the joints, triangular asperities inclined at angles of 9.5° and 18.5° to the shear movement were cast using high-strength gypsum plaster and infilled with clayey sand. These joints were sheared cyclically under constant normal stiffness conditions. It was found that, for a particular normal stiffness, the shear strength is a function of the initial normal stress, initial asperity angle, joint surface friction angle, infill thickness, infill friction angle, loading direction and number of loading cycles. Based on the experimental results, a mathematical model is proposed to evaluate the shear strength of infilled rock joints in cyclic loading conditions. The proposed model takes into consideration different initial asperity angles, initial normal stresses and ratios of infill thickness to asperity height.     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran     

New Methodology to Characterize Shear Behavior of Joints by Combination of Direct Shear Box Testing and Numerical Simulations

A new procedure is presented, which combines big shear box tests on rocks and corresponding numerical simulations with explicit consideration of joint roughness to get deeper insight into the shear behavior of rock joints. The procedure consists of three parts: (1) constant normal load- or CNS-shear box tests with registration of shear- and normal-components of stress and displacements and deduction of basis rock mechanical parameters; (2) high resolution 3D-scanning of joint surface to deduce joint topography; and (3) set-up, run and evaluation of 3-dimensional numerical model with explicit duplication of joint roughness as back-analysis of shear box tests. The numerical back-analysis provides deeper insight into the joint behavior at the micro-scale. Several parameters can be deduced, like micro-slope angle distribution, aperture size distribution, local normal stress distribution and detailed analysis of dilation in relation to shear direction. The potential of the new procedure is illustrated exemplary by shear box tests on slate.