考虑结构性的土体小应变本构模型研究

摘要

在岩土工程中绝大部分土体的变形都在小应变范围内；而在小应变范围内土体具有较高的刚度、逐渐屈服和结构性等变形特性。以剑桥模型为代表的传统土体本构模型难以描述土体的小应变特性。在这种情况下,采用传统土体本构模型分析基坑开挖、地铁施工等岩土工程问题时,预测计算结果将与实际情况相差较大。同时土体是天然地质体,土体结构对包括小应变在内的整个变形过程有强烈的影响。因此,建立考虑结构性的土体小应变本构模型具有重要的理论意义和现实意义。论文以国家自然科学基金项目“土体的小应变本构模型研究”(50309017)和重庆市自然科学基金项目“基于结构扰动的土体本构模型研究”(CSTC2005BB6114)为依托,从土体小应变变形特性的本质因为-土体的结构性入手,对土体的小应变本构模型进行研究。
　　 论文的主要研究工作与结论如下:
　　 ①对土体小应变情况下的变形特性及其本构模型、土体结构性及其本构模型和修正剑桥模型等研究进展进行分析,提出了论文研究思路--基于土体的结构性,对修正剑桥模型进行扩展,从而建立土体的小应变本构模型。
　　 ②基于非连续介质细观力学,采用理论分析和数值试验两种手段,分别对土体的小应变变形的结构性机理进行分析。
　　 建立了土体的细观模型,利用弹性接触力学理论讨论了胶结结构对土体小应变特性的影响规律。研究表明:在土体的整个变形过程中存在一个由胶结结构控制的变形区域；在该变形区域内,微量的胶结结构(胶结厚度为土体颗粒直径的万分之一)也会使土体宏观的模量增加10余倍。
　　 在颗粒流PFC2D平台上,建立了土体及其三轴试验的颗粒流模型。分别讨论了接触刚度、粘结强度、粘结刚度与半径等细观结构参数对小应变条件下变形特性影响。研究表明:在小应变范围内细观颗粒滑动趋势较小,摩擦力还未得到充分发挥,土体细观颗粒间的摩擦对小应变特性影响较小。细观颗粒间的胶结厚度、胶结半径和胶结模量等对小应变特性影响较大；随着细观颗粒间的胶结厚度、胶结半径和胶结模量的增加,土体宏观变形模量将增大。
　　 在小应变条件,微量胶结也可以使土体颗粒间点接触转换为面接触,从而使土体表现出较高的宏观模量；而当应变进一步增加时,胶结破损,土体细观颗粒间的面接触转换为点接触,从而使土体宏观模量迅速降低。
　　 土体宏观强度参数的粘聚力和摩擦系数两个方面,在小应变范围内主要是粘聚力发挥作用,而摩擦力相对影响较小。
　　 ③土体结构的模型及其本构关系研究
　　 对土体的物理结构、宏观变形特性等进行调查,提出了土体结构的材料结构模型和分布结构模型。
　　 将土体变形表现出的超固结划分为材料结构超固结与应力超固结这两种不同性质的超固结；从而使材料结构与应力超固结统一起来。然后对修正剑桥模型进行扩展,模型屈服面的尺寸由材料结构超固结参数和应力超固结参数两者之和共同决定,而屈服面尺寸大小的变化由塑性体积硬化与材料结构损伤共同描述；从而建立了一个可以统一考虑超固结和结构性的土体本构模型。该模型对于描述土体结构性,特别是在大应变范围内的结构性具有较高的适用性。
　　 基于土体的分布结构模型概念,假设土体的一维塑性硬化元件的摩擦滑动元件启动应力参数、塑性模量相互独立,并分别服从某种形式的统计分布；利用统计学和塑性力学理论,推导了基于分布结构模型概念的土体一维弹塑性本构模型。该一维弹塑性模型能够反映土体小应变时的逐渐屈服特征。
　　 ④考虑土体的材料结构,利用复合体损伤理论,提出了弹性变形部分的计算方法。在塑性变形部分模型基于土体的分布结构模型,认为塑性变形的内变量不是单一、固定的值,而是在一定区间上遵循统计分布密度函数分布的；扩展修正剑桥模型单一内变量的单个屈服面至尺寸连续的屈服面场,并采用相关联的流动法则,从而提出了可以反映土体小应变特性的土体本构模型。
　　 ⑤在有限差分FLAC2D平台上,对论文小应变本构模型进行了fish语言程序化,然后对天府基坑工程案例进行了数值分析。分析表明,基坑开挖工程中绝大部分土体都处于小应变变形情况,土体小应变本构模型有利于准确预测岩土工程的变形；相比修正剑桥模型,本文小应变本构模型更能反映工程土体的变形情况特性。

目录

文摘

英文文摘

1 绪论

1.1 问题提出

1.2 小应变范围的说明

1.3 国内外研究现状

1.3.1 土体小应变试验研究及其本构模型研究进展

1.3.2 土体结构性及其本构模型研究进展

1.3.3 剑桥模型及其描述小应变变形的不足

1.3.4 研究进展小结

1.4 主要研究内容与技术路线

1.4.1 主要研究内容

1.4.2 研究技术路线

2 土体小应变变形的结构性机理研究

2.1 引言

2.2 土体的细观结构调查

2.3 胶结结构对小应变特性的影响机理的理论分析

2.3.1 细观分析模型

2.3.2 胶结对小应变变形特性影响的理论推导

2.3.3 胶结对小应变变形特性影响讨论

2.4 土体及其三轴试验的颗粒流模型的建立

2.4.1 土体颗粒及粒团的颗粒流模拟

2.4.2 土体细观颗粒间相互作用的模拟

2.4.3 土体的三轴试验颗粒流模型的建立

2.5 基于颗粒流模型的土体小应变变形机理分析

2.5.1 颗粒间摩擦对小应变特性的影响分析

2.5.2 颗粒间接触刚度对小应变特性的影响分析

2.5.3 颗粒间接触粘结强度对小应变特性的影响分析

2.5.4 颗粒间平行粘结的对小应变特性的影响分析

2.5.5 围压对小应变特性的影响分析

2.6 本章结论

3 土体结构模型及其本构关系研究

3.1 土体结构的模型描述

3.1.1 材料结构模型

3.1.2 分布结构模型

3.1.3 土体分布结构随机变量的分布函数

3.1.4 分布函数与概率密度函数

3.2 基于材料结构模型的土体本构关系

3.2.1 材料结构超固结与应力超固结概念的提出

3.2.2 一般弹塑性公式

3.2.3 结构性土体的塑性变形

3.2.4 模型参数与算例

3.2.5 考虑材料结构本构模型的讨论

3.3 基于分布结构模型的一维弹塑性本构关系

3.3.1 塑性硬化元件的逐渐屈服分析

3.3.2 一维逐渐屈服弹塑性本构模型

3.3.3 算例

3.4 本章结论

4 考虑结构性的土体小应变本构模型

4.1 多个屈服面同时激活条件下的弹塑性表达

4.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.4 模型参数及其讨论

4.5 算例

4.6 本章结论

5 土体小应变模型在工程中的应用研究

5.1 论文模型程序化

5.1.1 有限差分的计算方法

5.1.2 自定义本构模型的执行方法

5.1.3 基于fish语言的自定义小应变本构模型算法

5.2 工程实例

5.2.1 工程概况

5.2.2 数值模型

5.2.3 结果分析

5.2.4 讨论

5.3 本章结论

6 结论及建议

6.1 结论

6.2 建议

致谢

参考文献

附录

A.作者在攻读学位期间发表的论文目录

B.作者在攻读学位期间参加科研项目