湿陷性黄土渠道浸水变形的数值模拟.pdfVIP

水 利 学 报 2005 年 3 月 SHUILI XUEBAO 第 36 卷 第 3 期 文章编号: 0559-9350(2005)03-0309-05 湿陷性黄土渠道浸水变形的数值模拟 1,2 1,2 1 陈铁林 ，沈珠江 ，杨代泉 (1.南京水利科学研究院 土工研究所，江苏 南京 210024；2.清华大学 水利水电工程系，北京 100084) 摘要：天然土体一般具有原始结构，在其结构被破坏前后会表现出不同的力学特性，这些特性很难用传统的弹塑 性理论来进行模拟。本文提出能够反映胶结和吸力两个因素影响的结构性黄土本构模型，采用损伤函数和屈服函 数反映天然黄土的湿陷变形。这一模型把变形过程中的天然黄土看作不同大小土块的集合体，总的变形由土块的 弹性变形、土块之间滑动引起的塑性变形和土块破碎引起的损伤变形三部分组成。塑性变形常用屈服函数描述， 损伤变形则引入一种类似的损伤函数加以描述。本文在严格的非饱和土力学理论的基础上，开发相应的有限元程 序，模拟一个黄土渠道典型断面的湿陷变形，分析了黄土渠道浸水变形的基本特性。 关键词：湿陷性；本构模型；非饱和土 中图分类号：TU444 文献标识码：A 黄土渠道的浸水变形对西北地区的灌溉事业造成很大危害，从20世纪50年代开始就受到广泛的注意。 经过地质学家、水利工程师和岩土工程师的长期努力，人们对黄土湿陷变形的机理已经有了比较深入的了 解[1]。但是，到目前为止，黄土渠道的岩土工程设计主要还是依靠经验，除了偶见一般性的边坡稳定分析 [2，3] [4] 和开挖过程的有限元模拟外 ，没有看到专门研究浸水引起的边坡失稳和变形分析的成果。沈珠江 曾用 有限元分析过羊毛湾大坝坝基黄土的浸水变形问题，但是当时未严格按非饱和土固结理论进行计算，工作 比较粗糙。有的研究者试图在弹塑性理论框架内建立适用于黄土湿陷变形的本构模型，但由于忽略了黄土 具有结构性的本质，因此是不成功的。例如著名的Alonso非饱和土本构模型，虽然通过建立 LC(Loading-Collapse)屈服面可以模拟湿陷变形，但算出的变形量随围压增大不断增大，不能得出有峰值 的湿陷曲线。天然黄土能保持大孔隙，既与颗粒之间存在一定胶结有关，也与黄土在比较干燥条件下形成 从而颗粒之间存在相当大的吸力有关。胶结矿物有水溶性，半水溶性和水稳性之分，浸水时吸力的丧失与 水溶性和半水溶性胶结的破坏是黄土湿陷性的根本原因。为此，近年来胶结和吸力两个因素同时存在对黄 [5] 土湿陷性影响的研究已经受到了重视 。 本文在严格的非饱和土力学理论的基础上，采用能够反映胶结和吸力两个因素影响的结构性黄土本构 模型，开发相应的有限元程序，模拟一个黄土渠道典型断面的湿陷变形，以期客观地反映黄土渠道浸水变 形的基本特性。 1 结构性黄土本构模型 该模型把黄土看作由胶结不均匀块体组成的结构体，在荷载和浸水的作用下，大土块逐渐碎裂成小块， 并最终粉碎成与重塑土类似的完全损伤土。土体的总应变由三部分组成，即土块的弹性应变、土块相互滑 移产生的塑性应变和土块碎裂引起的损伤应变，为 收稿日期：2004-01-18 基金项目：水利部重点基金项目(SZ9807)；中国博士后科学基金资助项目 作者简介：陈铁林(1970-)，男，安徽天长人，博士，主要从事岩土材料本构模型和土工数值计算研究。 1 水 利 学 报 2005 年 3 月 SHUILI XUEBAO 第 36 卷 第 3 期