高拱坝枢纽工程的长期安全稳定运行不仅是实现地区经济稳定发展的保障，还关系到水电站下游千百万人民群众的生命财产安全。受原始高地应力与工程扰动，拱坝坝肩岩石内部节理裂隙发育。电站蓄水后，两岸坝肩不仅承受坝体的拱推力，也承受来自库区数百米高的水压力。在拱推力和高水压力的长期共同作用下，坝肩软弱层间带裂隙岩石强度参数降低，易导致高坝枢纽结构失稳、渗透破坏和岸坡滑移等安全稳定问题。在20世纪，国外曾相继发生过水库大坝失事事件，如1926年美国圣佛西斯坝由坝基岩石遇水软化崩解造成大坝溃决；1959年法国马尔帕塞拱坝由坝基与坝基岩石弱化导致大坝失稳和破坏，造成421人死亡，经济损失达300×10⁸法郎^[1]。因此，研究水岩耦合作用下的坝肩软弱层间带裂隙岩石强度参数弱化效应对高坝枢纽的长期安全稳定运行至关重要。

针对水对岩石强度的弱化效应，国内外学者开展了大量研究。Hadizadeh等^[2]探讨了水对岩石的软化作用。Heggheim等^[3]通过研究海水和不同浓度的盐水对灰岩力学性质的影响发现，岩石中的矿物成分通过与水中的离子发生反应可改变岩石的成分和结构，进而改变其力学性质。Risnes等^[4]也研究了灰岩与水的耦合作用，并获得了相应的力学性质变化情况。王思敬等^[5]综述了水库地区水岩作用类型和机制，提出岩体结构面上渗压主要通过降低有效法向应力降低结构面抗剪强度；结构面抗剪强度与结构面的摩擦系数和内聚力有关，并取决于结构面上的法向应力和渗压。文[6-10]提出水-岩作用是微细观结构变化导致宏观力学特性劣化的过程。更多学者主要从试验角度研究水对岩石强度的弱化效应。汤连生等^[11]在试验中获得岩石结构、矿物成分、水溶液的化学性质与流动速率等影响水岩相互作用的主要因素，并从岩石的损伤机制和定量分析2个角度研究了水岩耦合的水化学损伤。刘新荣等^[12]的试验研究表明：在水作用下，岩石的Poisson比增大，其黏聚力和抗剪强度均变小。文[13-14]进行了水-岩耦合作用下砂板岩的抗剪强度弱化效应试验，研究发现：水压对黏聚力的影响最大，而对内摩擦角的影响极小，对变形模量的影响在两者之间。杨宝全等^[15]研究了水对坝肩软弱岩石的软化作用，并对坝肩稳定性进行了分析。文[16-18]也通过试验研究了水-岩耦合作用下岩石的力学特性。

学者们关于不同岩石的水压弱化效应研究对理论探索与工程建设具有十分重要的意义。然而，在岩石强度参数弱化试验研究中，学者们普遍关注完整岩石强度的变化规律，而对裂隙岩石与水压的讨论相对较少。本文认为裂隙岩石结构面强度存在一个临界值，即水压存在一个临界水压力；当水压大于该临界水压力时，抗剪强度不会随着水压增大继续降低，而是相对恒定，该临界值仅与内摩擦角有关。为验证这一设想，本文首先开展了裂隙岩石结构面强度的水压弱化试验与结果分析，获得临界水压力的理论依据与数据支撑，其次通过相关文献数据验证了临界水压力的存在，最后基于水压的弱化效应，研究了考虑水压的岩石统计损伤本构模型。

1 裂隙岩石水压弱化试验设备与方法 1.1 岩样制备

试验采用MTS815 Flex Test GT电液伺服岩石力学测试系统。为尽量避免大理岩岩样之间的差异，试验采用的不同岩样均取自结构构造比较均匀的同一块大理岩，岩样的制备按照实验室规范完成。具体过程如下：先将现场采取的大理岩制备成直径50 mm、高度100 mm的岩样，再沿与中轴线同一平面的岩样两端各钻取直径5 mm、深度80 mm的2个渗流孔，岩样渗水孔位置如图 1所示，2个渗流孔呈中心对称分布，距岩样边界各10 mm。设置目的是使水通过进口端小孔进入岩样，并通过岩样中的裂隙网络流向出口端小孔，从而保证水在岩样中形成渗流场。

1.2 试验方案及过程

根据相关工程实际，岩石试验渗压确定为0.5、1.0、1.5、2.0 MPa共4级，文中后续所述的水压均指岩石所受的渗透水压。试验围压确定为3.0、6.0、9.0、12.0 MPa共4级，试验方案与文[13-15]中试验过程相似，主要分为4个步骤。

步骤1   模拟大理岩受力变形破裂过程，采用60.0~100.0 MPa围压对岩样进行全过程三轴试验以模拟形成裂隙岩石，此阶段持续至轴压不随应变增加而变化时，即进入残余强度阶段；在残余强度阶段，监测岩样的变形模量，使其变形模量和强度参数接近工程实际参数，此时岩样的力学性能主要由结构面控制，可近似为现场天然岩体^[15]。

步骤2   模拟河流下蚀、卸荷的过程，将围压由60.0 MPa缓慢降低至弱化试验所需围压；降低围压的过程近似于河流下蚀、卸荷的过程；从试验技术方面考虑，此时轴压略高于围压，大理岩试验围压分别为3.0、6.0、9.0、12.0 MPa。

步骤3   天然状态下岩体试验，将围压降低至上述对应值后，进行三轴压缩试验，按国标及规程要求的速率加载至即将屈服时转为侧向应变控制，至应力-应变曲线转平后结束，获得天然状态下裂隙岩石的强度参数与应力-应变曲线^[15]。

步骤4   不同水压下的水岩耦合试验，在前述围压和轴压不变的基础上向试样内施加渗透水压，保证试样两端水压一致，达到岩样饱和的目的；岩样饱和后，在该级围压下进行三轴试验，按0.5、1.0、1.5、2.0 MPa逐级升高水压，试验过程按照规范要求进行，当轴压不随应变增加而变化时结束；再将轴压降低至初始轴压，水压升高至下一级后重复试验，直至水压2.0 MPa试验完成后结束全部试验，获得各级水压下裂隙大理岩的强度和变形参数^[14]。试验结果主要适用于裂隙岩石结构面强度参数的在渗透水压作用下的弱化效应。

1.3 试验应力-应变曲线

按照上述试验步骤对大理岩进行水岩耦合压缩试验，分析不同压缩阶段的试验数据成果。以围压为6.0 MPa的结果为例，不同水压条件下大理岩的应力-应变曲线如图 2a所示。围压6.0 Mpa、水压1.5 MPa的大理岩弱化试验典型应力-应变曲线如图 2b所示。图 2中σ₁、σ₃分别为轴压、围压，ε₁、ε_v、ε₃分别为轴向应变、体积应变、横向应变。

根据岩石试验数据和应力-应变曲线，获得的裂隙大理岩强度参数如表 1所示。由表 1可知，岩石的峰值强度呈现明显的水压弱化效应，即每一级围压下的峰值强度均是天然状态(无水压)最大，水压越大，岩石峰值强度越低，呈现随水压升高而减小的规律，水压达2.0 MPa时岩石峰值强度最低；岩石峰值强度呈现明显的围压强化效应，峰值强度随围压升高而增大，围压达12.0 MPa时岩石峰值强度最大。

2 裂隙岩石弱化试验结果与分析 2.1 水压弱化效应

为分析裂隙大理岩抗剪强度参数随水压变化的规律，本文以孔隙水压力P_w为横坐标，分别以摩擦系数f、黏聚力c为纵坐标，绘制了图 3中的大理岩抗剪强度参数与孔隙水压力关系曲线。

由图 3可知，摩擦系数随孔隙水压力增大有一定的升降，但整体上数值变化较小，说明孔隙水压力对摩擦系数的影响较小；而黏聚力随着孔隙水压力增大呈线性减小的变化规律，孔隙水压力对黏聚力影响较明显。因此，孔隙水压力对黏聚力具有明显的弱化效应。

为综合分析孔隙水压力对摩擦系数和黏聚力的影响，采用抗剪强度τ进行评定。假定正应力分别为5.0、10.0、15.0、20.0 MPa，根据Mole-Coulomb准则计算出对应的剪应力，获得的不同正应力下大理岩抗剪强度与孔隙水压力关系曲线如图 4所示。

对所获得的试验结果进行弱化率分析。弱化率表示如下^[15]：

$ W=\frac{X_{\mathrm{n}}-X_{\mathrm{w}}}{X_{\mathrm{n}}} \times 100 \%. $

(1)

其中：W为强度参数弱化率，X_n、X_w分别为天然状态和不同孔隙水压力下岩石的强度参数。

由式(1)可得，抗剪强度参数(f、c)的弱化率分别为W_f、W_c，正应力为10.0 MPa时抗剪强度的弱化率为W_τ，抗剪强度参数(f、c)与正应力为10.0 MPa时抗剪强度的弱化率随孔隙水压力的变化如图 5所示。

由图 5可知，随着孔隙水压力增大，摩擦系数的弱化率虽有差异，但整体变化幅度处于20%以内，变化趋势较平稳，这也进一步说明孔隙水压力对摩擦系数的影响较小；随着孔隙水压力增大，黏聚力的弱化率增加幅度较大，最大弱化率接近100%，这说明孔隙水压力对黏聚力具有明显的弱化效应，当孔隙水压力增大至一定值时，岩石可能完全没有黏聚力；随着孔隙水压力增大，岩石抗剪强度的弱化率明显增大。

2.2 弱化机理

根据2.1节试验结果与分析可知，孔隙水压力对岩石抗剪强度和黏聚力具有明显的弱化效应，这是从宏观角度得出的试验结果。本文从微观角度对大理岩岩样进行了电镜扫描，大理岩岩样不同倍数下电镜扫描图如图 6所示。

由图 6可知，大理岩内部存在大小不同的孔隙，随着放大倍数增加，孔隙逐渐清晰。而渗透水压在岩石内部表现为孔隙水压力的“楔入”作用^[19]，水压力促进了裂纹的扩展过程，使岩石产生渐进性破坏。

岩石内存在裂隙发育，而在试验过程中，为模拟实际工程中岩石所处围压环境和卸荷条件，先将岩石在60.0 MPa围压作用下加载荷载至峰值再卸载，此过程中岩石内部已产生大量裂隙，更有利于孔隙水压力的“楔入”。根据Mole-Coulomb强度理论与有效应力原理，当孔隙水压力作用于岩石时，有效正应力和岩石的抗剪强度分别表示如下：

$ \sigma^{\prime}=\sigma-\alpha^{\prime} P_{\mathrm{w}}, $

(2)

$ \begin{gathered} \tau^{\prime}=\sigma^{\prime} \tan \varphi+c= \\ \left(\sigma-\alpha^{\prime} P_{\mathrm{w}}\right) \tan \varphi+c=\sigma \tan \varphi+c_{\mathrm{w}} . \end{gathered} $

(3)

其中：σ′为有效应力，σ为轴向应力，α′为有效应力系数，τ′为岩石的抗剪强度，φ为内摩擦角，c_w为水影响后岩石的黏聚力。根据试验数据规律可知，c_w与孔隙水压力呈线性减小的趋势，可表示如下：

$ c_{\mathrm{w}}=c-\alpha^{\prime} P_{\mathrm{w}} \tan \varphi . $

(4)

由2.1节可知，孔隙水压力对摩擦系数的影响甚微。由式(4)可知，孔隙水压力的存在会使黏聚力减少α′P_wtanφ，岩石的抗剪强度也随之降低。

岩石的孔隙水压力可分为贯通孔隙水压力和封闭孔隙水压力^[20]。而岩石内平均孔隙水压力表示如下：

$ \bar{P}_{\mathrm{w}}=\alpha_1 P_1+\alpha_2 P_2. $

(5)

其中：P₁为封闭孔隙水压力，P₂为贯通孔隙水压力，α₁、α₂为比例系数。

岩石总应力大小与孔隙水压力有直接关系。当岩石在水压为0时，贯通孔隙水压力为0，平均孔隙水压力仅由封闭孔隙水压力组成；岩石强度仅与自身含水率有关。而岩石的弱化效应是岩石处于外界水压环境中表现出的相关特征差异。

岩石在外界水压环境下的封闭孔隙水压力变化较小，与零水压环境相近，而贯通孔隙水压力由外界渗透水压产生。贯通孔隙水压力取决于岩石内部的裂纹发育，随着加载应力增大，岩石微裂纹逐渐扩展，水通过扩展的裂纹渗入岩石内部，形成贯通孔隙水压力，该部分水压力是影响平均孔隙水压力的关键因素^[21]。当外界水压越大时，水越容易渗入岩石内部，从而增大平均孔隙水压力，降低岩石抗剪强度，岩石抗压强度也随之降低。同时，水的“楔入”作用加速裂纹扩展，进而加速岩石破坏。

当岩石在围压环境中时，围压可抑制岩石内部裂纹的扩展，延缓岩石的变形。围压越大，在同等应力作用下扩展的裂纹越少，水更难渗入岩石内部，贯通孔隙水压力也越小，岩石峰值强度也越大。因此，即使在水压作用下，围压对岩石强度依然有强化效应。水压与围压对岩石强度存在互斥作用。

综上所述，本文可从岩石内部孔隙水压力角度解释水压的弱化效应和围压的强化效应。水压对岩石的作用表现在对岩石破裂面剪切强度的影响，其实质是通过影响岩石破裂面的黏聚力和内摩擦角，从而影响岩石抗剪强度。水压对岩石的影响还可能与水和岩石内部矿物的化学作用，以及岩石内部裂纹贯通性、裂纹走向等裂缝发育特征多种因素有关，尚待进一步研究。

2.3 弱化临界水压力

水压的存在使黏聚力减少α′P_wtanφ，而试验中也发现，水压增大会明显弱化黏聚力，当水压达到某一定值后，岩石可能完全失去黏聚力。式(4)可改写为：

$ c_{\mathrm{w}}=c-k P_{\mathrm{w}}, $

(6)

$ k=\alpha^{\prime} \tan \varphi . $

(7)

其中k为水压与黏聚力关系曲线的斜率。水压对内摩擦角的影响较小，甚至可以忽略^[14]，式(7)中本文认为内摩擦角为一定值。故弱化后的黏聚力与水压呈线性关系。

由图 3b可知，随着水压增大，水压对黏聚力的弱化呈线性减小，故当水压达到一定值后，黏聚力减小至0，岩石完全丧失黏聚力^[16]。此时，c_w为0，根据式(6)可得：

$ P_{\mathrm{w}}=\frac{c}{k}=\frac{c}{\alpha^{\prime} \tan \varphi}. $

(8)

定义c_w为0时的水压力为弱化临界水压力P_r，表示如下：

$ P_{\mathrm{r}}=\frac{c}{\alpha^{\prime} \tan \varphi}. $

(9)

根据式(6)和试验数据可知，当c_w为0时，岩石完全丧失黏聚力，水压增大对岩石黏聚力失去作用，而水压变化对摩擦系数的影响甚微。因此，当水压大于临界水压力时，随着水压继续增大，岩石强度的弱化效应不明显。临界水压力的物理意义在于判定岩石黏聚力是否完全丧失，即水压对岩石强度弱化的临界点。当水压小于临界水压力时，岩石强度的水压弱化效应明显，黏聚力随着水压增大呈线性减小；当水压等于临界水压力时，岩石完全丧失黏聚力；当水压大于临界水压力时，随着水压增大，岩石强度不会明显降低，此时岩石强度的水压弱化效应不明显。

引入文[22]关于砂板岩、绿片岩相同试验方案的数据，得到的不同岩石弱化后的黏聚力与水压的线性关系如图 7所示。根据式(8)计算所得的3种岩石的弱化临界水压力如表 2所示。由表 2可知，绿片岩的弱化临界水压力达4.11 MPa，而其余2种岩石的弱化临界水压力均小于3.00 MPa。试验过程中未设置超过2.0 MPa的水压，所以无法验证水压大于临界水压力的弱化规律。

文[14]中砂板岩抗剪强度参数随水压变化如表 3所示。由表 3可知，当水压为4.0 MPa时，黏聚力弱化为0，但内摩擦角受水压影响较小。通过曲线拟合，可以得到该砂板岩的弱化临界水压力为3.8 MPa，这说明确实存在弱化的临界水压力。当水压大于临界水压力时，黏聚力仍为0，内摩擦角保持恒定，此时砂板岩抗剪强度不会随着水压增大而明显降低。

文[16]中细砂岩抗剪强度参数随水压变化如表 4所示。由表 4可知，该细砂岩的弱化临界水压力为3.48 MPa，当水压大于3.48 MPa时(如试验水压为6.0、9.0 MPa)，黏聚力均为0，内摩擦角变化较小。假定岩石在正应力10.0 MPa作用下，细砂岩的抗剪强度随水压的变化规律如图 8所示。由图 8可知，当水压小于临界水压力时，随着水压增大，细砂岩抗剪强度降低；当水压大于临界水压力时，随着水压增大，细砂岩抗剪强度基本没有变化。

通过以上分析可以表明：岩石强度随着水压升高并非无止境降低，存在一临界水压力，当水压大于该临界值时，岩石强度不会随着水压增大而继续降低，即为一个定值，该定值仅与内摩擦角有关。

由于岩石强度的水压弱化效应还受岩石内部裂纹发育、矿物组成等其他因素的综合影响，在同一围压下，水压增大可能会进一步降低岩石的峰值强度，该部分研究有待更进一步的试验与分析。

3 考虑水压弱化的岩石统计损伤本构模型 3.1 统计损伤本构模型

岩石材料中含大量裂隙、微裂纹等缺陷，其分布、尺寸、形状和贯通长度均表现出较明显的随机性，并且这些随机分布的缺陷影响岩石的强度，使强度也具有随机性。因此可将岩石微元体强度看作随机变量，对其概率分布进行研究。

岩石属于脆性材料，其强度具有较大的分散性，可采用瑞典学者Waloddi Weibull提出的Weibull分布函数对岩石强度进行统计处理。该分布的密度函数表示如下^[23]：

$ \phi\left(F^*\right)=\frac{m}{F_0}\left(\frac{F^*}{F_0}\right)^{m-1} \cdot \exp \left[-\left(\frac{F^*}{F_0}\right)^m\right] . $

(10)

其中：m、F₀为分布参数，可反映材料的力学特性；F^*为微元体强度；ϕ(F^*)为强度为F^*时微元体破坏的概率，可反映岩石的损伤程度。

损伤变量定义是研究岩石等材料在荷载作用下损伤演化规律的核心问题。损伤变量D^*定义为损伤面积S^*与材料全面积S之比^[24]，即

$ D^*=\frac{S^*}{S}=\int_0^{F^*} \phi\left(F^*\right) \mathrm{d} F^*=1-\exp \left[-\left(\frac{F^*}{F_0}\right)^m\right], $

(11)

$ F^*=\alpha I_1^*+\left(J_2^*\right)^{1 / 2}-K, $

(12)

$ \alpha=\frac{2 \sin \varphi}{\sqrt{3}(3-\sin \varphi)}, $

(13)

$ K=\frac{6 c \cos \varphi}{\sqrt{3}(3-\sin \varphi)}. $

(14)

其中：α、K为Drucker-Prager破坏准则参数，与黏聚力和内摩擦角有关；I₁^*为应力张量第一不变量；J₂^*为应力偏量第二不变量。I₁^*、J₂^*分别表示如下：

$ I_1^*=\sigma_1^*+\sigma_2^*+\sigma_3^*, $

(15)

$ J_2^*=\frac{\left(\sigma_1^*-\sigma_2^*\right)^2+\left(\sigma_2^*-\sigma_3^*\right)^2+\left(\sigma_1^*-\sigma_3^*\right)^2}{6} . $

(16)

其中σ_i^*为名义应力σ_i对应的有效应力，i为1、2、3。σ_i^*表示如下：

$ \sigma_i^*=\frac{\sigma_i}{1-D^*} . $

(17)

三轴压缩实验中，假定岩石符合广义Hooke定律，表示如下：

$ \left\{\begin{array}{l} \sigma_1=E \varepsilon_1+\mu\left(\sigma_2+\sigma_3\right), \\ \sigma_2=E \varepsilon_2+\mu\left(\sigma_1+\sigma_3\right), \\ \sigma_3=E \varepsilon_3+\mu\left(\sigma_1+\sigma_2\right) . \end{array}\right. $

(18)

其中：ε_i为主应变，i=1，2，3；E为弹性模量；μ为Poisson比。

围压作用会使岩石内部的微裂纹状态改变，故引入损伤变量进行分析。三维各向同性线弹性体损伤本构模型表示如下^[25]：

$ \{\boldsymbol{\sigma}\}=[\boldsymbol{C}]\{\boldsymbol{\varepsilon}\}\left(1-D^*\right) . $

(19)

其中：{σ}为名义应力矢量，[C]为材料弹性矩阵；{ε}为应变矢量。

结合文[26]，岩石微元体强度服从式(10)的Weibull分布建立的岩石损伤软化统计本构模型，岩石统计损伤本构模型如下：

$ \left\{\begin{array}{l} \sigma_1=E \varepsilon_1 \exp \left[-\left(\frac{F^*}{F_0}\right)^m\right]+\mu\left(\sigma_2+\sigma_3\right), \\ \sigma_2=E \varepsilon_2 \exp \left[-\left(\frac{F^*}{F_0}\right)^m\right]+\mu\left(\sigma_1+\sigma_3\right), \\ \sigma_3=E \varepsilon_3 \exp \left[-\left(\frac{F^*}{F_0}\right)^m\right]+\mu\left(\sigma_1+\sigma_2\right) . \end{array}\right. $

(20)

对于式(20)中的2个分布参数m和F₀，本文将采用文[26]的方法并根据应力峰值σ_c、应变峰值ε_c确定参数值，相关计算表示如下：

$ F_0=F^* \cdot m^{\frac{1}{m}}, $

(21)

$ m=\frac{1}{\ln \left[E \varepsilon_{\mathrm{c}} /\left(\sigma_{\mathrm{c}}-2 \mu \sigma_2\right)\right]}. $

(22)

3.2 考虑水压的统计损伤本构模型

根据岩石在不同水压下的弱化试验，水压对岩石抗剪强度参数的影响主要表现为水压对黏聚力具有明显的弱化效应，而其对内摩擦角的弱化效应不明显，一般可忽略。因此，考虑水压情况下，式(14)的K值可表示如下：

$ \begin{gathered} K_{\mathrm{w}}=\frac{6 c_{\mathrm{w}} \cos \varphi}{\sqrt{3}(3-\sin \varphi)}=\frac{6\left(c-\alpha^{\prime} P_{\mathrm{w}} \tan \varphi\right) \cos \varphi}{\sqrt{3}(3-\sin \varphi)}= \\ \frac{6 c \cos \varphi}{\sqrt{3}(3-\sin \varphi)}-\frac{6 \alpha^{\prime} P_{\mathrm{w}} \sin \varphi}{\sqrt{3}(3-\sin \varphi)} . \end{gathered} $

(23)

由式(23)可知，水压会使K_w一定程度上降低。但根据2.3节提出的弱化临界水压力原理，当水压等于临界水压力时，黏聚力弱化为0，随着水压增大，黏聚力不再变化。因此，根据Drucker-Prager破坏准则，随着水压增大，岩石强度在一定程度上降低，但当水压大于临界水压力时，随着水压增大，岩石强度不会发生明显的弱化现象。

在常规的岩石三轴压缩试验中，围压σ₂=σ₃，应变ε₂=ε₃，当考虑水压影响后，F^*发生变化，因此，根据式(20)可得考虑水压的统计损伤本构模型

$ \sigma_1=E \varepsilon_1 \exp \left[-\left(\frac{F_{\mathrm{w}}^*}{F_0}\right)^m\right]+2 \mu \sigma_2 . $

(24)

其中式(25)的关键是如何通过试验数据求解考虑水压的微元体强度F_w^*。

根据Drucker-Prager破坏准则，材料的宏观强度F可表示如下：

$ F=\alpha I_1+J_2^{1 / 2}-K . $

(25)

其中：I₁为宏观方面的应力张量第一不变量，J₂为宏观方面的应力偏量第二不变量。

$ I_1=\sigma_1+\sigma_2+\sigma_3, $

(26)

$ J_2=\frac{\left(\sigma_1-\sigma_2\right)^2+\left(\sigma_2-\sigma_3\right)^2+\left(\sigma_1-\sigma_3\right)^2}{6} \text {. } $

(27)

根据损伤定义及式(17), 考虑水压的岩石微元体强度可表示为

$ F_{\mathrm{w}}^*=\frac{F}{1-D^*}=\frac{\left(\alpha I_1+J_2^{1 / 2}-K_{\mathrm{w}}\right) E \varepsilon_1}{\sigma_1-\mu\left(\sigma_2+\sigma_3\right)} . $

(28)

至此，考虑水压的统计损伤本构模型中相关参数均能通过试验数据求解，其中黏聚力受水压弱化的关系可由式(4)得出，式(4)中-α′tanφ可根据黏聚力与水压的关系曲线求得。

3.3 本构模型试验数据验证

根据岩石弱化试验结果，可获得大理岩在1.0 MPa水压下不同围压的力学参数，如表 5所示。将图 3中c、φ值和相关试验数据代入式(11)、(12)、(13)、(14)、(21)、(22)，可得大理岩在不同围压作用下的分布参数F₀、m，以及大理岩峰值应力时的损伤变量，大理岩在不同围压下的分布参数与损伤变量(水压1.0 MPa)如表 6所示。

为验证模型的合理性，将分布参数m和F₀代入式(24)，可获得不同围压时的统计损伤本构模型。以试验应变值为模型输入值，采用本构模型求得应力值，并与试验应力值比较。图 9为各级围压下模型曲线与试验曲线的对比。由图 9可知，本文模型曲线能较好地反映岩石破裂全过程。峰后试验曲线与模型曲线的差异主要由试验过程中加载方式和加载速度改变造成。

4 结论

本文通过采用岩石力学测试系统对裂隙岩石进行水岩耦合作用下的三轴压缩试验，分析了水压对裂隙岩石的强度弱化效应，并提出了弱化临界水压力，构建了考虑水压弱化的岩石统计损伤本构模型，主要结论如下：

1) 水压对岩石抗剪强度参数的影响主要表现为其对黏聚力的弱化效应，随着水压增大，黏聚力的弱化率增加幅度较大，最大弱化率接近100%，弱化后的黏聚力随着水压增大呈线性减小；水压对岩石的摩擦系数影响较小，最大弱化率不超过20%，且无明显的升降规律。

2) 黏聚力弱化为0时的水压力为弱化临界水压力：当水压小于临界水压力时，岩石强度的水压弱化效应明显，黏聚力随着水压增大呈线性减小，直至0；当水压大于临界水压力时，随着水压增大，岩石强度不会明显降低，此时岩石强度的水压弱化效应不明显。临界水压力作为水压对岩石强度及黏聚力弱化的界限值，可用于界定水压对岩石强度的影响规律。

3) 通过构建考虑水压弱化效应的岩石统计损伤本构模型，并采用试验数据验证了理论模型，该模型能实现水岩耦合作用下岩石破裂全过程的演变规律计算，证明了模型的有效性与合理性。