考虑损伤的节理本构模型

本文在弹塑性损伤的理论框架内,讨论了节理等地质间断面的本构模型。这个模型能够反映节理面的损伤弱化,扩容和弹性刚度劣化等复杂特性。这个模型的另一优点是,塑性变形增量与屈服面是非正交的,但本构矩阵具有对称性。这种对称性在岩石力学的理论研究和数值分析中是至关紧要的。     

基于复合损伤的节理岩体动态本构模型研究

基于运用霍普金森压杆（SHPB）装置对节理岩体动载试验得出的数据,从节理面倾角、贯通度、厚度、组数、填充物及应变率等不同方面分析各因素对节理岩体力学特性的影响,通过对节理岩体在高应变率下的损伤机制和破坏形式进行分析的基础上,基于复合损伤理论,对广义Bingham模型进行改进,构造了节理岩体材料在不同应变率下动态响应的本构模型,对模型计算和试验结果的比较表明,该模型能够很好地描述节理岩体动荷载下初始弹性变形阶段、稳态塑性变形阶段和加速变形破坏阶段的应力-应变关系,并且理论与试验结果吻合较好,从而证明了该模型的正确性和合理性。     

节理岩体本构模型与强度理论的统计断裂力学分析

本文提出了节理岩体的三维本构模型、强度理论及破坏概率的表达式。这里把节理岩体视作等效的各向异性连续介质,其具有裂隙张量S_(ifst)和柔度张量C_(ifst)。为了获得本构模型、强度判据及破坏概率,这里引用了能量平衡方程ω=ω_0+ω_c、断裂力学能量判据、岩体节理统计理论及可靠性理论,并引入了剩余剪应力τ_n,有效应力σ_(if)及最可能出现的最大节理半径a_m等概念。这套理论的初步应用是令人满意的。     

基于摩尔-库仑准则的断续节理岩体复合损伤本构模型

针对工程实际中断续节理裂隙岩体的损伤本构模型,假设岩石微元强度服从Weibull随机分布,以摩尔-库仑破坏准则作为描述微元强度的表示方法,推导出细观损伤变量。利用能量和断裂力学理论,综合考虑节理几何特征及力学特性,推导宏观损伤变量计算公式。基于Lemaitre应变等效假设,考虑宏细观缺陷耦合作用,推导出复合损伤变量,建立了基于摩尔-库仑准则的宏细观缺陷耦合作用的断续裂隙岩体损伤本构模型。研究结果表明：（1）采用摩尔-库仑准则作为描述微元强度的统计分布变量建立的损伤模型能够较好地反映岩石内部缺陷分布和变形特征,该模型真实地反映岩石微元强度受应力状态的影响。（2）该模型建立的理论曲线与断续节理岩体的试验曲线吻合较好。（3）节理裂隙岩体宏观损伤变量及峰值强度随节理倾角的变化规律与综合考虑宏细观耦合作用下的损伤变量及裂隙岩体峰值强度随节理倾角的变化规律基本一致。（4）宏细观耦合作用下的等效弹性模量与节理贯通率呈非线性负相关;在节理倾角一定的情况下,损伤变量与节理长度呈非线性正相关;在贯通率较小时,岩体的宏观损伤变量与内摩擦角的关系呈线性负相关变化,贯通率达到一定程度,线性关系变成非线性关系。     

双轴应力下非贯通节理岩体压缩损伤本构模型

针对目前节理岩体损伤变量定义中大多仅考虑节理长度、倾角等几何性质,而未考虑节理抗剪强度等力学性质的不足,基于断裂力学中的由于单个节理存在引起的附加应变能增量与损伤力学中的损伤应变能释放量相关联的观点,推导出了在双轴应力下含单条非贯通闭合节理岩体的损伤变量计算公式,并根据断裂力学理论对双轴压缩荷载下的单个节理尖端应力强度因子计算方法进行了研究,得出了应力强度因子KⅠ、KⅡ的计算公式;同时考虑多节理间的相互作用,给出了单组单排及单组多排非贯通节理尖端应力强度因子计算公式,由此建立了相应的节理岩体双轴压缩损伤本构模型,并利用该模型进行了相应的算例分析。结果表明：对含单条非贯通闭合节理的岩体而言,当节理倾角小于其内摩擦角时,岩体强度与完整岩石相同,岩体损伤为0,而后随着节理倾角增加,岩体强度、损伤随节理倾角的变化分别呈开口向上及向下的抛物线,当节理倾角约为60°时,岩体损伤最大,强度最低。随着节理长度增加,岩体损伤增加,而随着节理内摩擦角的增加,岩体损伤则减小;对含单组单排非贯通闭合节理的岩体而言,当节理总长度一定时,随着单条节理长度的减小及节理条数的增加,岩体损伤则逐渐减小,但其减小幅度与节理条数并不呈线性关系。     

单轴压缩荷载下非贯通闭合节理岩体损伤本构模型

目前损伤力学已被认为是研究节理岩体力学行为的有效工具,但是在目前的节理岩体损伤变量定义中大多仅考虑节理几何特征而未考虑节理内摩擦角等力学参数,这显然不能很好地反映节理岩体的力学特征。为此,拟推导出一个能够综合考虑节理几何及力学参数的损伤变量（张量）,并由此建立单轴压缩荷载下岩体损伤本构模型。首先,基于断裂力学的由于单个节理存在引起的附加应变能增量与损伤力学的损伤应变能释放量相关联的观点,推导出了含非贯通节理岩体的损伤变量计算公式;其次,根据断裂力学理论对单轴压缩荷载下的单个节理尖端应力强度因子计算方法进行了研究,得出了应力强度因子K_Ⅰ、K_Ⅱ的计算公式;同时考虑多节理间的相互作用给出了单组单排及多排非贯通节理应力强度因子计算公式。最后,利用该模型对含单条非贯通节理的岩体在单轴压缩荷载作用下的峰值强度及损伤变量进行了分析计算。结果表明,当节理倾角小于其内摩擦角时,岩体强度与完整岩石相同,岩体损伤为零,而后随着节理倾角增加,岩体强度、损伤随节理倾角的变化分别呈开口向上及向下的抛物线,当节理倾角约为60°时,岩体损伤最大,强度最低。随着节理长度增加,岩体损伤增加,而随着节理内摩擦角的增加,岩体损伤则减小。     

节理卸荷岩体非线性本构关系再研究

对地质材料模型、石膏模型、砂浆模型、进行了模型试验,通过对模型试验得出的数据进行分析、拟合,得出初步的节理卸荷岩体的增量本构。针对其不足,从参数选择、物理意义等方面进一步进行修正得出一改进的增量本构关系,并在有限元分析程序ADINA中接入该本构关系,对清江隔河岩电站厂房高边坡三维模型进行有限元分析。     

考虑损伤门槛的统计损伤本构模型研究

在统计损伤理论的基础上,从统计损伤本构模型的一般性公式出发,考虑到岩石的受力状态,提出统计损伤本构模型应考虑损伤门槛的影响,建立了考虑损伤门槛的统计损伤本构方程。以石膏角砾岩的常规三轴试验结果作为实例,进行了验证。理论计算结果与试验结果对比表明,考虑损伤门槛的统计损伤本构模型是更加合理的。     

裂隙岩体中弹性波传播特性试验及宏细观损伤本构模型研究

以不同加载方式下裂隙岩体中弹性波传播特性试验为基础,首先通过对比试验中弹性波波速、波幅与岩体宏观损伤的关系,给出了采用弹性波波幅来表征的裂隙岩体宏观损伤变量;其次利用统计强度理论定义了岩石的细观损伤变量,并基于连续损伤理论建立了宏、细观裂隙耦合的裂隙岩体损伤本构模型;最后结合试验数据对该模型进行了验证与分析。研究表明:完整岩体中细观裂隙的起裂和扩展对岩体中传播弹性波的波幅和波速的影响较小,而裂隙岩体中宏观裂隙的张开和闭合对弹性波波幅和波速的影响较大。裂隙岩体在外力作用下时,弹性波波速和波幅的变化规律类似,相较于弹性波波速,弹性波波幅对裂隙岩体宏观损伤的变化表现的更敏感,可用来定义裂隙岩体的宏观损伤变量。不同加载方式下裂隙岩体的损伤特性与宏、细观裂隙以及其所处的应力状态相关,基于宏、细观裂隙耦合的裂隙岩体损伤本构模型可以较好地反映不同加载方式下裂隙岩体的损伤力学特性。在三轴压缩试验中,宏观裂隙对岩体轴向压应力方向上损伤的影响主要体现在试验的中后期阶段,在三轴拉压试验中,宏观裂隙对岩体轴向拉应力方向上损伤的影响贯穿于试验全过程中。     

节理岩体爆破数值模型及模拟研究

运用多刚体系统动力学及接触冲击理论建立了节理岩体弹塑性动态本构模型。基于DYNA－2D程序原理及其框架,利用所建模型,模拟了节理岩体的爆破过程,量化了节理岩体在爆炸载荷作用下不同时刻的应力场。     

基于宏细观损伤耦合的非贯通裂隙岩体本构模型

针对非贯通裂隙岩体工程结构中的受荷岩体,提出受荷细观损伤与裂隙宏观损伤的概念。以完整岩石的初始损伤状态作为基准损伤状态,综合考虑裂隙宏观缺陷的存在、微裂纹细观缺陷在受荷下的损伤扩展以及宏细观缺陷在受荷过程中的耦合,基于Lemaitre应变等效假设,推导了考虑宏细观缺陷耦合的复合损伤变量,并给出同时考虑试件尺寸、裂隙几何与力学特性的宏观损伤变量的计算公式,从而建立了基于宏细观缺陷耦合的非贯通裂隙岩体在荷载作用下的损伤本构模型。用宏细观损伤耦合的本构模型来描述非贯通裂隙岩体在受荷过程中的细观损伤演化与宏观损伤行为,与非贯通裂隙岩体实际受荷情况符合较好。研究结果表明：（1）完整岩样和裂隙岩样的应力-应变行为在峰值强度之前差异较大,峰值强度以后差异逐渐减小,最后趋于一致,二者具有相近的残余强度;（2）裂隙岩体强度随裂隙贯通率的增加而增大,随裂隙倾角的变化具有明显的各向异性,同时还与裂隙面的内摩擦角有关;（3）裂隙倾角为90°时,裂隙岩样的峰值强度最高;张开型裂隙岩样的裂隙倾角为45°时,峰值强度最低;（4）非贯通裂隙岩体工程结构中的受荷岩体,其力学性能由受荷细观损伤与裂隙宏观损伤及其耦合效应所决定,基于宏细观损伤耦合的复合损伤变量可以较好地反映非贯通裂隙岩样的力学特性。     

节理岩体等效损伤流变模型初步研究

对节理岩体等效损伤流变模型进行了初步研究,假设岩石为各向同性体各向同性损伤,节理面的法向和切向损伤不同,分别建立了岩石和节理面流变损伤演变函数。假设材料无损,计算应力采用有效应力,推导了节理岩体等效损伤流变模型有限元计算公式,编制了相应的有限元程序。算例分析表明,采用的节理岩体等效损伤流变模型是可行的,可较好地反映节理岩体的损伤演变过程。     

A 3D model of a jointed rock mass and its deformation properties

Conclusion A geometric model of a jointed rock mass and some formulae of strain and deformation parametersE, G and were given. They may be useful for evaluating the deformation parameters of rock masses and classifying the rock mass quatitatively.     

节理岩体各向异性等效弹性参数研究

采用复合材料力学分析方法,针对规则节理岩体考虑坐标转换及刚度系数,建立各向异性等效弹性参数的解析公式。通过算例证明了解析公式的可行性,同时对规则节理岩体等效参数规律性进行了研究,分析在不同的节理刚度比值下等效弹性参数的变化规律。结合白鹤滩坝址区柱状节理岩体特征,探讨了其等效弹性参数取值和各向异性特征。研究表明：所采用方法在刚度取值相等时岩体可退化为各向同性介质,且该方法计算结果与实测结果较为吻合;同时发现柱体偏转对等效弹性参数的影响较大,偏转后水平向和铅直向弹性模量分别降低约54%和17%,根据曲线规律性,偏转角度越大,降低程度越高,等效弹性参数的曲线变化反映了其具有明显的各向异性特征。     

节理岩体强度特征研究综述

随着岩体工程建设规模的加大,对节理岩体的研究提出了更高的要求,需要对节理岩体强度及其相关问题有全新的认识。在认真总结国内外节理岩体理论与试验研究的基础上,以节理岩体强度为主线,详细论述了理论强度准则、经验强度准则及其与强度有关的节理岩体试验、本构、数值模拟等内容,指出了节理岩体需进一步研究的几个关键问题。     

Microplane damage model for jointed rock masses

The paper presents a new microplane constitutive model for the inelastic behavior of jointed rock masses that takes into account the mechanical behavior and geometric characteristics of cracks and joints. The basic idea is that the microplane modeling of rock masses under general triaxial loading, including compression, requires the isotropic rock matrix and the joints to be considered as two distinct phases coupled in parallel. A joint continuity factor is defined as a microplane damage variable to represent the stress‐carrying area fraction of the joint phase. Based on the assumption of parallel coupling between the rock joint and the rock matrix, the overall mechanical behavior of the rock is characterized by microplane constitutive laws for the rock matrix and for the rock joints, along with an evolution law for the microplane joint continuity factor. The inelastic response of the rock matrix and the rock joints is controlled on the microplane level by the stress–strain boundaries. Based on the arguments enunciated in developing the new microplane model M7 for concrete, the previously used volumetric–deviatoric splits of the elastic strains and of the tensile boundary are avoided. The boundaries are tensile normal, compressive normal, and shear. The numerical simulations demonstrate satisfactory fits of published triaxial test data on sandstone and on jointed plaster mortar, including quintessential features such as the strain softening and dilatancy under low confining pressure, as well as the brittle–ductile transition under higher confining pressure, and the decrease of jointed rock strength and Young's modulus with an increasing dip angle of the joint. Copyright © 2014 John Wiley & Sons, Ltd.     

循环冻融条件下节理岩体损伤破坏试验研究

通过循环冻融和相似材料试验,研究了节理岩体在冻融条件下的损伤破坏机制及其相应的力学特性。通过对经历冻融循环后的试件损伤破坏模式的观察和单轴压缩试验,重点研究了节理倾角、节理贯通度、节理条数、节理充填物厚度、节理充填物类型、试件饱和度、冻融循环次数等对试件冻融损伤破坏模式、单轴抗压强度和弹性模量的影响。研究发现：节理存在及其物理力学性质对岩体的冻融损伤破坏模式及强度均有很大影响。节理倾角通过影响冻融裂纹的起裂位置进而影响其破坏模式和强度;随着贯通度的增加,试件表面裂纹由稀变密;随着节理条数增加,试件受冻融影响明显加剧;随着节理充填物厚度增加,试件冻融损伤程度先增加后减小;节理充填物类型对试件冻融损伤程度也有一定影响;随着试件饱和度的增加,试件冻融损伤程度先减小后增加;随着冻融循环次数的增加,节理试件表面因冻胀引起的裂纹逐渐增多、变宽,且其抗冻融特性较完整试件差。上述研究成果对寒区岩体工程建设及安全运营具有重要的参考价值。     

充水节理岩体中大断面隧洞开挖的数值模型研究

采用非饱和的渗流应力耦合模型分析了糯扎渡水电站2#导流洞的开挖过程,导流洞部分洞身穿过节理带,洞顶水头较高。详细地讨论了高水头和节理带对施工过程的影响。节理带采用描述高密度平行节理组的各向异性节理本构来逼近;考虑开挖引起的介质变形对渗透系数的影响;考虑排水引起的饱和度变化对渗透系数的影响。所用的非饱和瞬态耦合模型可以模拟出开挖引起的EDZ区域孔隙水压力急剧升高、有效应力减小、渗透系数动态的变化以及排水对洞室稳定性的提高。数值模拟的计算结果与国外类似试验的一般性观测结论相吻合,因此,可以用来评价水位以下隧洞施工方法和施工速度的合理性和经济性。