为探究炭质板岩在不同含水率条件下的蠕变特性，该文基于室内分级加载蠕变试验数据，通过非线性黏性元件与表征屈服特性的塑性元件并联，构建了可描述蠕变全过程中非线性加速蠕变阶段的黏塑性体；进而将该非线性黏塑性体与经典西原模型串联，结合弹性模量和黏性系数的软化规律引入了4个损伤因子，建立了能描述不同含水率条件下炭质板岩蠕变全过程的损伤蠕变本构模型（即改进的西原模型）；通过ABAQUS有限元软件进行用户自定义材料子程序（UMAT）二次开发，实现了不同含水率条件下炭质板岩三轴蠕变试验的数值模拟，并通过对比岩样蠕变试验与数值模拟结果，验证了模型的适用性。研究结果表明：所构建的改进西原模型能显著提升加速蠕变阶段的模拟精度，不同含水率工况下的分级加载蠕变试验曲线与数值模拟曲线吻合良好，相关系数均大于0.9；数值模拟云图准确复现了不同含水率岩样的蠕变变形演化过程，验证了所提考虑含水损伤的炭质板岩蠕变本构模型及UMAT子程序开发的正确性与有效性。该研究成果可为富水深埋炭质板岩隧道工程的长期稳定性评估与灾害预警提供理论支撑。

To explore the creep characteristics of carbonaceous schist under different moisture content conditions, based on the indoor graded loading creep test data, this paper constructed a viscoplastic body capable of describing the nonlinear accelerated creep stage throughout the entire creep process. This is achieved by parallelly connecting a nonlinear viscous element with a plastic element that characterizes yield behavior, based on data from indoor graded loading creep tests. Subsequently, this nonlinear viscoplastic body was integrated in series with the classical Nishihara model. By incorporating the softening patterns of elastic modulus and viscosity coefficient, four damage factors were introduced to establish a damage-based creep constitutive model (i.e., an improved Nishihara model) that describes the entire creep process of carbonaceous slate under different moisture content conditions. Through secondary development of a user-defined material subroutine (UMAT) in ABAQUS finite element software, numerical simulations of triaxial creep tests on carbonaceous slate under varying moisture content conditions were implemented. The applicability of the model was validated by comparing experimental creep data from rock samples with numerical simulation results. The research results show that the constructed improved Nishihara model can significantly improve the simulation accuracy in the accelerated creep stage. The graded loading creep test curves under different moisture content conditions are in good agreement with the numerical simulation curves, and the correlation coefficients are all greater than 0.9. The numerical simulation cloud map accurately reproduced the creep deformation evolution process of rock samples with different water contents, verifying the correctness and effectiveness of the proposed creep constitutive model of carbonaceous shale considering water damage and the development of the UMAT subroutine. This research achievement can provide theoretical support for the long-term stability assessment and disaster early warning of deep-water buried carbonaceous slate tunnel projects.