岩体结构面新型本构模型研究

对于涉及需要单独考虑岩体结构面的工程岩体结构分析,采用能反映岩体结构面主要力学特性的合理的本构模型是取得合理解答的关键问题之一。针对在经典连续介质力学理论框架内建立岩体结构面本构模型的缺点,基于岩体结构面的实际受力变形特性,采用直接法建立了一种新型岩体结构面本构模型。所建立的模型依据岩体结构面切向应力变形曲线及剪胀曲线的实际特征,将其分为峰前线性段、峰前非线性段以及峰后软化段,并分别给出了用于描述岩体结构面变形和强度等主要力学特性的数学模型,进而推导建立了结构面各变形阶段的增量型本构模型。最后,编写相关计算程序,采用所建立的新型本构模型以及被广泛采用的Plesha模型对经典的岩体结构面直剪试验成果进行拟合分析。结果表明,所建立的新型本构模型能更为合理的描述岩体结构面的主要力学特性,且模拟能力优于Plesha模型     

锚固节理直剪力学行为的颗粒流数值分析

为探究锚固节理在直剪加载作用下的细观力学响应,利用颗粒流数值计算方法建立不同锚固角度下的锚固节理数值模型,并进行不同法向荷载下的直剪试验。之后,通过分析和对比剪切-位移曲线和峰值剪切强度来对锚固节理宏观力学性质进行研究。与此同时,基于微裂纹分布规律,从细观角度揭示了剪切荷载下不同锚固角度的岩石节理面的破坏特征;研究结果显示。（1）不同于非锚固节理,锚固节理的剪切-位移曲线在剪切后期呈现出一定的增长趋势,且粗糙度不同曲线变化特征存在一定的差异。（2）不同锚固角度下,锚固节理的破坏模式存在着较为明显的差异。当锚固角为90°时,锚固节理模型的破坏主要集中于锚杆与上下节理面接触位置,且破坏主要为挤压式破坏。随着锚固角度的减小,接触位置的破碎区域不断减小,节理破坏主要表现为沿着锚杆的轴向拉伸变形破坏及内部的拉伸破坏。（3）锚固角度的变化对锚固节理的抗强度的影响程度与节理面的粗糙度存在一定关联。具体而言,较为平直节理抗剪强度随着锚固角度的增大,增长趋势明显,而粗糙节理面随着锚固角度的增大峰值抗剪强度的变化趋势相对较为平缓。     

节理化岩体复合式滑移破坏的模型试验研究

针对节理化岩体的自身特点,提出沿岩体层面产生的应力集中以及层面与节理之间的残余完整岩桥破坏,将可能导致岩体产生复合式滑移破坏。并在室内利用岩石三轴刚性伺服机,通过模型试验对岩体的层面和一组与层面正交的节理的情况进行模拟研究,分析了结构面性状和围压对复合式滑移破坏的影响,给出了在贯通性结构面(层面)倾角、非贯通性结构面以及围压等因素共同作用下,节理化岩体发生复合式滑移破坏的条件以及破坏面特征。在与前人的试验资料进行对比分析的基础上得到可靠的结论,期望为在将来的实践工作中正确评估节理化岩体的破坏模式提供指导性意见。      

节理化岩体复合式滑移破坏的模型试验研究

针对节理化岩体的自身特点,提出沿岩体层面产生的应力集中以及层面与节理之间的残余完整岩桥破坏,将可能导致岩体产生复合式滑移破坏。在室内利用岩石三轴刚性伺服机,通过模型试验对岩体的层面和一组与层面正交的节理的情况进行模拟研究,分析了结构面性状和围压对复合式滑移破坏的影响,给出了在贯通性结构面（层面）倾角、非贯通性结构面的以及围压等因素共同作用下,节理化岩体发生复合式滑移破坏的条件以及破坏面特征。最后在与前人的试验资料进行对比分析的基础上得到可靠的结论,期望为在将来的实践工作中正确评估节理化岩体的破坏模式提供指导性意见。     

Cosserat介质理论针对节理岩体的应用

由于考虑了内在尺度的影响,Cosserat介质理论较传统连续介质力学更精确和一般化,在考虑具有微结构介质的力学行为时拥有优势.节理岩体的力学性质明显受结构面影响控制,传统连续介质理论对于认识结构应力、应变的真实状况是有限度的.Cosserat理论应用于岩土介质模拟时,根据具体情况有三种方式一种是将岩土介质视为粒状材料,第二种则是将节理岩体看作是层状材料,此种模型应用最为常见;另一种则是将高度节理化岩体看作是块体结构.与试验或其它数值计算的对比表明,采用层状和块状模型对节理岩体进行数值模拟得到的分析成果是合理可信的.空间模型的应用也作了讨论.     

金属矿开采节理演化诱发覆岩变形失稳机制研究

金属矿开采诱发覆岩变形是一个复杂的力学问题,在工程活动的扰动下,采空区围岩的原始应力状态将被打破,引起覆岩的应力重分布,从而达到新的应力平衡状态。当覆岩中存在一些软弱结构面时,这些软弱结构面将不断地发展、演化,并可能出现贯通,其对覆岩的变形和破坏具有明显的影响。针对节理条件这一变量,采用Udec二维离散元模型,分析了双岩桥模型和四岩桥模型中各岩桥破裂特征及机制,并对比了两组模型中顶板应力释放规律及特点。采用3组对比模型实验,分析了0°/90°,45°/45°,20°/70°节理条件下岩移特点及变形机制。基于数值模型,对镜儿泉矿区实际节理条件下的地表变形特征和沉陷范围进行了研究,模型结果与地表变形GPS监测结果具有很好的一致性。     

非均质共面断续节理岩体拉伸剪切破裂机制研究

工程开挖面附近卸荷扰动区的岩体,受结构面和拉应力共同影响作用,其变形和破坏具有拉剪复合特征。为研究节理岩体的拉剪力学特性,基于颗粒离散元法针对共面断续节理岩体开展了系列数值模拟研究。通过假设粒间接触的力学参数服从Weibull分布表征岩体的非均质性,探讨了非均质性、均质度、法向拉应力和节理连通率对节理岩体拉剪强度和破坏模式的影响。研究表明:拉剪应力条件下非均质性节理岩体主要沿阶梯型破裂面破坏,剪应力-水平位移曲线可以分为线性变形阶段、非线性变形阶段、峰值及峰后阶段;随均质度提高,节理岩体的剪切强度逐渐增加且提升幅度逐渐减弱,趋于均质岩体,岩体中微裂纹由弥散型分布向破裂面集中;节理岩体峰值剪切强度和法向拉应力的大小呈非线性负相关关系;岩体剪切强度随节理连通率增加而显著降低。     

节理倾角对黄土力学特性影响试验研究

黄土节理种类多样、角度多变,是影响黄土力学性质的重要因素.为研究节理对黄土力学性质影响,开展含预制节理黄土试样单轴与三轴压缩强度试验,分析了非贯通节理对试样剪切破坏模式影响,探讨了节理倾角对黄土强度与变形特性影响规律.结果表明:节理试样剪切破坏均为压剪破坏,破坏模式可分为节理面与破裂面贯通型和节理面与破裂面斜交型,压力...     

考虑截面规则性的柱状节理岩体变形及强度特性研究

柱状节理岩体的复杂结构使其现场力学参数的确定存在困难,可靠地估计柱状节理岩体强度及变形特性对工程安全至关重要。结合Voronoi图形及3D打印技术制作不规则柱状节理岩体光敏树脂模具,并制备了具有不同倾角的规则及不规则柱状节理岩体试样。基于开展的单轴压缩试验结果,分析不同截面柱状节理岩体的变形及强度特性差异。根据试样最终破坏形态,揭示柱状节理岩体破坏模式及机制。结合已有研究结果,描述节理面材料对柱状节理岩体破坏形态及强度的影响。基于传统节理系数方法,提出了更能反映柱状节理岩体柱体结构特征的修正节理系数方法。采用修正节理系数法建立规则及不规则柱状节理岩体试验结果与现场变形及强度参数间的经验关系,将所提出计算公式应用于白鹤滩水电站工程,预测结果与现场实测值及RMR法（RMR为岩石质量评分指标）、Q法（Q为岩体质量指标）的预测值进行比较。结果表明,对于现场变形和强度参数,所提出的修正节理系数法的预测结果与现场实测结果较为吻合。     

基于单元劈裂法的岩体黏结型结构面数值模拟

节理岩体中富含结构面等不连续体,其中很多结构面之间具有一定黏结强度,在冲击荷载作用下会发生脱结（debond）行为,从而消耗岩体应变能。在该结构面之间的黏结性质在很大程度上决定着岩体的破坏行为,发生脱结后,结构面之间不再具有黏结强度,岩体的力学性能发生弱化。为了对该类结构面进行模拟,在改进的Xu-Needleman势函数基础上推导出结构面单元,并将其嵌入到单元劈裂法中,模拟结构面的开裂过程,当结构面完全脱结后结构面单元就转化为一般节理单元。相应地,在数值实现过程中只是将结构面单元替换为一般节理单元即可。该结构面单元与单元劈裂法相结合,能够有效地模拟节理岩体的破坏过程。     

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.     

岩石节理表面几何特性的三维统计分析

在利用三维非接触式高精度激光表面形状测定仪精确测量岩石节理的表面粗糙形状并进行数值化表达的基础上,应用地理信息系统(GIS)技术,实现了岩石节理粗糙面的三维可视化。再对岩石节理粗糙面的三维几何特性参数（诸如表面粗糙度、裂隙张开度和粗糙高度等）进行统计分析并计算频数分布的特征数字。最后,根据频率直方图的分布趋势和频数分布的特征数字研究各几何特性参数的分布规律。结果表明,花岗岩节理粗糙表面起伏比砂岩大,但二者几何特性参数的分布规律相似,节理粗糙高度和裂隙张开度分布近似于正态分布规律,而节理粗糙面倾斜角分布近似于Γ分布规律。     

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.     

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.     

Observation of the Degradation Characteristics and Scale of Unevenness on Three-dimensional Artificial Rock Joint Surfaces Subjected to Shear

The present study explores the degradation characteristics and scale of unevenness (small-scale roughness) on sheared rock joint surfaces at a low-stress regime. While the degradation characteristics of unevenness and the normal stress are mutually interrelated, an understanding of the degradation patterns of the three-dimensional roughness of rock joints is one of the important components needed to identify asperity failure characteristics and to quantify the role of damaged unevenness in establishing a shear strength model. A series of direct shear tests was performed on three-dimensional artificial rock joint surfaces at different normal stress levels. After shearing, the spatial distributions and statistical parameters of degraded roughness were analysed for the different normal stress levels. The length and area of the degraded zones showed bell-shaped distributions in a logarithmic scale, and the dominant scale (or the most frequently occurring scale) of the damaged asperities (i.e., unevenness) ranged from approximately, 0.5 to 5.0 mm in length and 0.1–10 mm² in area. This scale of the damaged unevenness was consistent regardless of the level of normal stress. It was also found that the relative area of damaged unevenness on a given joint area, and thus the contribution of the mechanical asperity failure component to shear strength increased as normal stress increased.     

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

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

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

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

A constitutive model for a laboratory rock joint with multi-scale asperity degradation

This paper presents a joint constitutive model that considers separately the mechanical contribution of waviness and unevenness of a joint to shear behaviour. The critical asperities for waviness and unevenness are determined from geometric properties in a lab-scale joint. The wear process is employed to model the degradation in dilation and strength during shear. From dimensional analysis, asperity degradation constants are developed using geometric parameters including asperity angle, wavelength, and amplitude as well as rock strength and stress. The applicability of the proposed model was assessed by performing direct shear tests on three joint roughness coefficient (JRC) profiles and providing its correlation with experimental results. Additionally, experimental data taken from literature were used to validate the model’s performance.     

Towards an understanding of joint roughness

Summary It is argued that the currently available joint models are incapable of accurate predictions of joint shear behaviour without resorting to substantial levels of empiricism. This is because these models fail to adequately quantify joint roughness, or appreciate the importance of scale. A novel approach, which uses fractal geometry to investigate joint roughness, is described. This approach goes beyond describing the symptoms of roughness and seeks to find the cause. In the application of this approach, the concepts of fractal geometry, fractal dimension and self similarity are described and used as a framework to formulate a statistically based and practical model for the characterisation of rock joint roughness. Important relationships between the fractal dimension and the more useful statistical parameters of standard deviation of both asperity angles and asperity heights are derived. These relationships not only provide a useful, working method for quantifying joint roughness, but are also shown to provide a basis for understanding the Barton empirical JRC-JCS model. In addition, the fractal model is able to provide conceptual models for the effects of normal stress on the shear behaviour of joints and the scale-dependence of joints.     

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.