水是影响岩体工程安全性的重要因素之一,水与岩体之间不仅存在力学作用,还存在复杂的水-岩物理化学作用,许多岩体工程的破坏和失稳均与水-岩作用密切相关。随着我国大规模的岩体工程建设,岩体工程在荷载和水-岩作用下的应力-渗流-化学(M-H-C)耦合问题日趋重要。因此,开展岩体剪切特性的水-岩作用效应、岩石M-H-C耦合本构模型以及M-H-C耦合数值模拟的研究,对于受水-岩作用影响的岩体工程如边坡、坝基、地下洞室及核废料地下储存、溶浸采矿、油气田开采等问题,均具有重要的理论意义与应用价值。本文围绕存在水-岩作用的岩体工程稳定性问题,采用试验研究、理论分析与数值模拟相结合的方法,就水-岩作用对岩体剪切特性的影响效应开展了系统的试验研究,分析了水-岩作用对岩体剪切特性的影响规律,探讨了其作用机理；建立了岩石水化学损伤动态演化数学模型,开展了岩石水化学损伤动态演化的数值模拟研究；探讨了岩石M-H-C耦合的基本理论与相关问题；提出了岩石M-H-C耦合弹性本构模型和M-H-C耦合弹塑性本构模型,分别编制了其数值模拟程序,建立了相应的数值模拟方法,同时进行了工程算例分析。论文的主要研究内容如下：1.开展了不同水化学溶液侵蚀条件下砂岩的抗剪强度试验,分析了水化学作用对砂岩微细观结构的影响效应,探讨了砂岩的水化学损伤机理。在此基础上,分析了水化学作用对砂岩抗剪强度的影响规律,同时通过引入岩石水化学损伤度,定量表达了砂岩抗剪强度参数随水化学损伤的变化过程。2.开展了具不同含水率的砂岩软弱结构面剪切蠕变试验,分析并探讨了含水率对砂岩软弱结构面剪切蠕变特性及长期抗剪强度的影响规律与机理。研究表明,在相同应力水平下,随着含水率增大,其蠕变变形和蠕变速率显著增大,达到稳态阶段所需时间明显延长,初始蠕变强度和长期抗剪强度逐渐降低,且长期抗剪强度随含水率变化符合指数函数关系。3.根据试验结果,对黏弹塑性元件组合蠕变模型进行了改进,建立了能考虑含水率影响效应的砂岩软弱结构面剪切蠕变模型,同时采用进化优化算法对该蠕变模型中的多参数进行了反演。4.采用化学动力学方法建立了岩石水化学损伤动态演化数学模型。在此基础上,对工程时间尺度内重庆砂岩在水化学环境下的损伤演化过程进行了数值模拟,分析了水溶液酸碱度、离子种类与离子浓度对重庆砂岩损伤演化的影响效应与规律。5.提出了由水化学损伤动态演化模型、岩石固相基质本构方程、渗流方程、孔隙压力变化方程等构成的岩石M-H-C耦合弹性本构模型。6.基于所提出的M-H-C耦合弹性本构模型和模块化思想,利用多物理场耦合分析软件COMSOL的流-固耦合模块和水文地球化学分析软件Phreeqc,分别对其进行二次开发,编制相关的接口程序,并通过岩石水化学损伤度将Phreeqc的水化学损伤计算过程与COMSOL的应力-水力耦合计算过程有机融合,构建起岩石M-H-C耦合的弹性数值模拟程序与方法。7.通过引入化学软化概念,在D-P屈服准则中加入化学软化力,提出了考虑水化学作用效应的岩石屈服准则及加卸载准则,采用关联流动法则,建立了相应的弹塑性本构模型,进而构建起包括水化学损伤动态演化模型、渗流方程、弹性常数和渗透系数的水化学损伤演化方程等的岩石M-H-C耦合弹塑性本构模型。8.基于所提出的岩石M-H-C耦合弹性数值模拟方法,通过嵌入上述M-H-C耦合弹塑性本构模型,构建了岩石M-H-C耦合弹塑性数值模拟程序与方法。9.采用所提出的M-H-C耦合弹性和弹塑性数值模拟方法,对某一存在水-岩化学作用的隧道稳定性分别开展了弹性和弹塑性计算分析,模拟了隧道在不同酸碱度、不同离子种类与离子浓度水环境下的损伤演化过程,分析了水环境酸碱度、离子种类与离子浓度对其稳定性的影响效应,并就弹性和弹塑性结果进行了比较分析。