长江科学院院报 ›› 2024, Vol. 41 ›› Issue (7): 132-138.DOI: 10.11988/ckyyb.20230149

• 岩土工程 • 上一篇    下一篇

考虑非饱和土基质吸力-回弹指数相关性的弹塑性本构模型及有限元实现

程昊1, 熊大生1, 唐辉明2, 张抒2   

  1. 1.中铁第四勘察设计院集团有限公司,武汉 430010;
    2.中国地质大学,武汉 430010
  • 收稿日期:2023-02-17 修回日期:2023-04-13 出版日期:2024-07-01 发布日期:2024-07-08
  • 通讯作者: 熊大生(1961-),男,湖南常德人,正高级工程师,主要从事土力学方面的研究。E-mail:1079778756@qq.com
  • 作者简介:程 昊(1988-),男,湖北武汉人,工程师,博士,主要从事非饱和土力学、本构模型、数值计算等方面的研究。E-mail:chenghao216@qq.com
  • 基金资助:
    中铁第四勘察设计院集团有限公司科研课题(2020k158);湖北巴东地质灾害国家野外科学观测研究站开放基金课题(BNORSG202213)

Swelling Index and Matric Suction Relations for Unsaturated Soil Elasto-plastic Constitutive Model in FEM

CHENG Hao1, XIONG Da-sheng1, TANG Hui-ming2, ZHANG Shu2   

  1. 1. China Railway Siyuan Survey and Design Group Co.,Ltd., Wuhan 430010, China;
    2. China University of Geosciences,Wuhan 430010, China
  • Received:2023-02-17 Revised:2023-04-13 Published:2024-07-01 Online:2024-07-08

摘要: 非饱和土压缩指数受基质吸力影响明显,回弹指数是土的临界状态本构模型中的重要变形参数。但在本构模型理论研究中,回弹指数随基质吸力的变化影响往往被忽略。为解决这一矛盾,以杭州地区重塑粉土为试样,使用GDS非饱和直剪系统的竖向加载功能,对试样进行50、100、200 kPa 3种定吸力条件下一维固结回弹试验,获取不同基质吸力非饱和回弹指数κ(s),并拟合得到κ(s)函数关系,对比已有试验研究,验证了拟合结果的合理性。进而将κ(s)表达式引入非饱和土扩展剑桥弹塑性本构模型,形成一种考虑回弹指数随基质吸力变化的非饱和土弹塑性本构模型。应用子增量步显示积分算法将该模型编辑成ABAQUS有限元UMAT子程序,最后基于ABAQUS有限元平台结合自编UMAT子程序进行固结回弹数值试验,模拟实际试验过程,对比两者试验结果,验证了该UMAT子程序的正确性。实现了一种考虑回弹指数随基质吸力变化的非饱和土弹塑性本构模型的数值化,为其进一步应用建立基础。

关键词: 非饱和土, 基质吸力, 回弹指数, 本构模型, 有限元模型

Abstract: The swelling index, a crucial parameter for soil deformation, holds significant importance in the critical state elasto-plastic constitutive model. Unlike the compressibility index, which has been extensively reported in literature, the variation of the swelling index with matric suction has often been overlooked in theoretical studies of unsaturated soil constitutive models. To address this gap, oedometer tests were conducted on remolded silt specimens sourced from Hangzhou, China. Under suction-controlled conditions, utilizing an oedometer equipped with a GDS direct shear chamber, the experimental program comprised three tests conducted at matric suction levels of 50, 100, and 200 kPa. The relationship between matric suction and the swelling index, denoted as κ(s), was established and fitted. Validation was performed by comparison with previous experimental studies. Furthermore, this κ(s) relationship was incorporated into an unsaturated soil-extended Cam-clay elasto-plastic model. An explicit integration with a substepping error control algorithm was implemented within the ABAQUS finite element method (FEM) UMAT subroutine. The numerical modeling results of the silt oedometer tests, coupled with the UMAT subroutine, exhibited high agreement with laboratory tests.

Key words: unsaturated soil, matric suction, swelling index, constitutive model, finite element method

中图分类号: