华南理工大学学报(自然科学版) ›› 2026, Vol. 54 ›› Issue (2): 133-144.doi: 10.12141/j.issn.1000-565X.250055

• 土木建筑工程 • 上一篇    下一篇

不同含水率炭质板岩的损伤蠕变本构模型

胡涛涛1,2(), 刘可萌1, 赵玉龙1, 李浩1, 高咸超1, 王磊2,3()   

  1. 1.长安大学 公路学院,陕西 西安 710064
    2.长安大学 陕西省公路桥梁与隧道重点实验室,陕西 西安 710064
    3.西安科技大学 建筑与土木工程学院,陕西 西安 710054
  • 收稿日期:2025-03-05 出版日期:2026-02-25 发布日期:2025-07-18
  • 通信作者: 王磊 E-mail:tthu@chd.edu.cn;wl2013@xust.edu.cn
  • 作者简介:胡涛涛(1985—),男,博士,副教授,主要从事隧道工程研究。E-mail: tthu@chd.edu.cn
  • 基金资助:
    国家自然科学基金项目(52378388);陕西省重点研发计划项目(2025SF-YBXM-554);陕西省公路桥梁与隧道重点实验室(长安大学)水平提升项目(QLYSD2024K07);长安大学中央高校基本科研业务费专项资金项目(300102213211);长安大学青年学者学科交叉团队建设项目(300104240912)

Damage Creep Constitutive Model of Carbonaceous Slate with Different Moisture Content

HU Taotao1,2(), LIU Kemeng1, ZHAO Yulong1, LI Hao1, GAO Xianchao1, WANG Lei2,3()   

  1. 1.School of Highway,Chang’an University,Xi’an 710064,Shaanxi,China
    2.Key Laboratory for Bridge and Tunnel of Shaanxi Province,Chang’an University,Xi’an 710064,Shaanxi,China
    3.School of Architecture and Civil Engineering,Xi’an University of Science and Technology,Xi’an 710054,Shaanxi,China
  • Received:2025-03-05 Online:2026-02-25 Published:2025-07-18
  • Contact: WANG Lei E-mail:tthu@chd.edu.cn;wl2013@xust.edu.cn
  • Supported by:
    the National Natural Science Foundation of China(52378388);the Key Research and Development Program of Shaanxi Province(2025SF-YBXM-554)

摘要:

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

关键词: 炭质板岩, 含水率, 损伤, 改进西原模型, UMAT子程序

Abstract:

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.

Key words: carbonaceous slate, moisture content, damage, improved Nishihara model, UMAT subroutine

中图分类号: