非饱和土弹塑性本构模型数值实现及其在斜坡变形分析中的应用

论文摘要

土由固体颗粒和孔隙组成,当土中孔隙被液体充满时为饱和土,而当土中孔隙未被液体充满即同时存在液体和气体时为非饱和土。自然界中,非饱和土广泛分布于地下潜水位以上的表层土体中。当发生降雨或地下水位升降等水文环境变化时,非饱和土含水量会随之发生变化,引起土饱和度的变化,进而导致土的力学特性发生改变,一方面土饱和度的增加会造成土有效应力减小导致土强度降低,另一方面土饱和度的增加会改变土的状态使其发生自身材料属性的改变例如屈服面的扩张或收缩等现象。在工程实践中斜坡变形破坏是无法回避的问题,而土质斜坡由于地形等条件的特殊性,造就了其表层大量分布非饱和土的现象。大量工程实践表明,降雨或地下水升降等水文环境变化是造成斜坡变形乃至失稳破坏的重要因素,因此研究斜坡在渗流应力耦合作用下的变形演化规律具有重要意义。目前,已有的研究渗流应力耦合作用下斜坡变形破坏规律多基于有效应力原理,而考虑饱和度变化导致的土材料属性的改变的文献较少,因此研究非饱和土的本构关系,特别是研究非饱和土材料力学属性随饱和度的变化规律即本构模型,并应用其研究斜坡变形演化规律具有重大理论意义和工程实践指导价值。而在非饱和土本构模型研究方面,研究非饱和土本构关系的文献较多,但研究非饱和土本构关系数值实现以及考虑非饱和土水力滞回效应的文献较少,因此研究考虑水力滞回效应的非饱和土本构模型数值实现方法具有更重要的基础研究意义。鉴于非饱和土应力应变关系在斜坡变形演化规律研究中的重要性,以及现有研究的不足,本文考虑饱和度对土材料属性的改变,结合非饱和土水力滞回效应,应用基于误差控制的子增量步显示积分算法对已有非饱和土扩展剑桥模型进行Abaqus有限元二次开发,编制模型Umat用户材料子程序。在此基础上基于非饱和土扩展剑桥模型对饱和度的响应原理,建立一种非饱和土扩展边界面模型,并同样应用基于误差控制的子增量步显示积分算法进行Abaqus有限元二次开发,编制了非饱和土扩展边界面模型的Umat用户材料子程序。随后应用非饱和土扩展剑桥模型和非饱和土扩展边界面模型对比模拟研究了理想斜坡模型在强降雨入渗条件下,在瞬态渗流应力耦合作用下斜坡变形演化规律,并与传统修正剑桥模型的模拟结果进行对比。获得了有限元模拟非饱和土质斜坡变形演化规律的新方法,并验证了模拟方法的可靠性。主要研究成果及结论如下:(1)非饱和土扩展剑桥模型的显式数值实现结合考虑水力滞回效应的非饱和土扩展剑桥弹塑性本构模型的屈服面方程、强化准则、流动和法则,利用塑性一致性条件,将饱和度视为一应变分量,推导得到扩展剑桥模型的一致性切线模量矩阵。在此基础上,应用基于误差控制的自动子增量步算法,每个子增量内部利用改进欧拉法一步显示求解本构常微分方程。以简洁明了的算法实现了非饱和土扩展剑桥模型Abaqus本构Umat子程序二次开发。应用非饱和土扩展剑桥模型进行多种条件下数值实验,验证了模型程序的正确性,并探究了模型适用性,最终得到该模型对非饱和土多项力学行为的模拟适用性良好的结论,而在干湿循环过程中由于模型考虑了水力滞回效应因此模拟干湿循环是非饱和土力学特性有较好的适用性。(2)非饱和土扩展剑桥模型的建立和显式数值实现以非饱和土扩展剑桥模型的强化准则为基础,建立一种适用于非饱和土的扩展边界面本构模型。为进行有限元数值实现,结合塑性一致性条件,将饱和度视为一应变分量,推导该模型一致切线模量。应用基于误差控制的自动子增量步显式积分算法,与非饱和扩展剑桥模型一样,无需进行迭代,在子增量步内部使用改进欧拉法对边界面本构常微分方程进行显示一步求解。从而实现非饱和土扩展边界面模型的Abaqus有限元二次开发。应用非饱和土扩展边界面模型Umat子程序进行多种条件下的数值试验,验证了模型程序的正确性,并探究了模型性质及适用性,非饱和土扩展边界面模型适用于模拟多种非饱和土力学特性,由于该模型考虑了边界面效应和水力滞回效应,因此对非饱和土在应力循环荷载和干湿循环荷载条件下力学响应的模拟有更强的适用性。(3)滑带土非饱和直剪试验参数拟合取三峡库区黄土坡滑坡滑带土进行非饱和定吸力直剪试验,应用非饱和土扩展剑桥模型及非饱和土扩展边界面模型对滑带土定吸力直剪试验进行拟合,并获取相应的模型参数。从拟合结果可以看出,扩展剑桥模型模拟结果有明显的弹性段,塑性段主要表现为应变硬化。而从扩展边界面模型模拟结果看,应力应变曲线没有纯弹性段,整个试验过程均为弹塑性变形,曲线形态初段为应变硬化,后段为应变软化,相比于扩展剑桥模型扩展边界面模型拟合结果更接近于实验曲线。(4)强降雨入渗条件下理想斜坡变形演化研究对理想斜坡模型,采用滑带土直剪试验拟合参数,应用非饱和土扩展剑桥模型和扩展边界面模型使用瞬态渗流应力耦合分析过程模拟强降雨入渗条件下的变形演化特征,并与传统饱和土修正剑桥模型应用渗流应力耦合分析的模拟结果进行对比。通过对比发现,由于受饱和度变化影响,在相同降雨时长非饱和土扩展边界面模型位移和扩展剑桥模型位移均明显大于传统修正剑桥模型,并且相对修正剑桥模型,应变分布明显受降雨入渗影响,在坡表更为集中。三种模型模拟斜坡变形结果均呈现出随时间呈阶段性演化的特点,即初始变形阶段、稳定变形阶段和加速变形阶段直至破坏。其中扩展剑桥模型模拟的位移结果演化阶段划分更为明显,而扩展边界面模型稳定变形和加速变形阶段分异不明显,主要是由于扩展边界面模型材料一直处于弹塑性变形,没有截然的弹塑性界面,因此该模型更适合模拟流动性强的粘土材料。

论文目录

作者简历

摘要

abstract

第一章绪论

1.1 选题的来源、目的和意义

1.2 选题的国内外研究现状

1.2.1 非饱和土弹塑性本构模型研究现状

1.2.2 非饱和土弹塑性模型数值模拟研究现状

1.2.3 土水特征曲线研究现状

1.2.4 渗流应力耦合分析斜坡稳定性研究现状

1.3 存在的问题

1.4 主要研究内容、技术路线和创新点

1.4.1 主要研究内容

1.4.2 技术路线

1.4.3 创新点

第二章基本理论

2.1 非饱和土弹塑性本构模型基本理论

2.1.1 土的压缩性

2.1.2 非饱和土巴塞罗那模型

2.1.3 非饱和土Wheeler模型

2.1.4 非饱和土土水特征曲线

2.1.5 非饱和土有效应力原理及达西渗流定律

2.2 渗流应力耦合有限元计算基本原理

2.2.1 渗流应力耦合有限元法计算基本方程

2.2.2 渗流应力耦合基本方程解法及本构模型计算流程

2.3 Abaqus及相关子程序简介

2.3.1 Abaqus有限元软件简介

2.3.2 Umat和 Usfld子程序简介

2.4 本章小结

第三章非饱和土扩展剑桥弹塑性本构模型数值实现

3.1 扩展剑桥非饱和土本构模型简介

3.1.1 屈服面方程

3.1.2 流动法则与强化准则

3.2 非饱和土扩展剑桥弹塑性模型有限元实现

3.2.1 非饱和土扩展剑桥本构模型一致性切线模量矩阵

3.2.2 基于误差控制的自动子增量步显示积分算法

3.2.3 饱和度的获取

3.3 模型验证

3.3.1 饱和土三轴排水压缩模拟试验

3.3.2 三轴固结及体应变规律

3.3.3 不排水条件下三轴压缩模拟试验

3.3.4 不同恒定饱和度条件下三轴压缩模拟试验

3.3.5 三轴湿陷模拟试验

3.3.6 三轴干湿循环模拟试验

3.3.7 非饱和土三轴定吸力剪切模拟试验

3.4 本章小结

第四章非饱和土扩展边界面弹塑性本构模型及其数值实现

4.1 非饱和土扩展边界面弹塑性模型

4.1.1 边界面方程及映射规则

4.1.2 流动法则

4.1.3 强化准则及加载条件

4.2 非饱和土扩展边界面弹塑性模型有限元实现

4.2.1 非饱和土扩展边界面弹塑性本构模型一致切线模量矩阵

4.2.2 非饱和土扩展边界面弹塑性本构模型编程要点

4.3 模型验证及模型性质探究

4.3.1 饱和土三轴排水压缩模拟试验

4.3.2 三轴固结及体应变规律

4.3.3 不排水动三轴模拟

4.3.4 不排水三轴压缩模拟

4.3.5 不同恒定饱和度条件下三轴压缩模拟

4.3.6 三轴湿陷模拟试验

4.3.7 三轴干湿循环模拟试验

4.3.8 非饱和土定吸力三轴剪切模拟试验

4.4 本章小结

第五章非饱和土弹塑性模型在斜坡变形研究中的应用

5.1 滑带土非饱和直剪试验对比研究

5.1.1 试样采集

5.1.2 滑带土基本物理参数及土水特征曲线

5.1.3 非饱和直剪试验方法及模拟

5.1.4 非饱和直剪试验结果拟合及对比

5.2 降雨入渗条件下斜坡变形研究

5.2.1 理想斜坡模型设置

5.2.2 降雨入渗条件下斜坡变形演化

5.3 本章小结

第六章结论与展望

6.1 结论

6.2 展望

致谢

参考文献

文章来源

类型: 博士论文

作者: 程昊

导师: 唐辉明

关键词: 非饱和土,本构模型,有限元二次开发,斜坡变形演化

来源: 中国地质大学

年度: 2019

分类: 基础科学,工程科技Ⅱ辑

专业: 地质学,建筑科学与工程

单位: 中国地质大学

基金: 国家重点研发计划项目(NO.2017YFC1501305)

分类号: TU43

总页数: 136

文件大小: 8529K

下载量: 508

非饱和土弹塑性本构模型数值实现及其在斜坡变形分析中的应用

论文摘要

论文目录

文章来源

相关论文文献

猜你喜欢