含液多孔介质中失稳现象理论研究及应变局部化的有限元-无网格耦合方法

摘要

研究含液多孔介质的失稳现象，诸如基坑工程或隧洞的开挖过程中可能出现的地基沉降甚至坑壁垮塌；以及边坡或堤坝由地震或暴雨所引发的滑动破坏(通常称为滑坡现象或边坡失稳)具有很重要的工程和理论意义。本论文从理论分析和数值模拟两个方面致力于研究在静、动力荷载作用下含液(特别是饱和)多孔介质中驻波间断、颤振失稳和应变局部化等破坏现象和过程。 在多孔介质受到冲击或爆炸等高频模态占主导地位的脉动荷载作用时，人们需要研究介质中应力波的传播过程。工程中许多多孔材料为塑性应变软化材料。在应力波的传播过程中，介质内某处的受力状态将首先达到材料的极限承载能力，并伴随以在介质中局部狭窄区域内急剧发生非弹性应变为特征的应变局部化现象和承载能力的急剧下降。 本论文基于饱和与非饱和含液多孔介质的非线性动力-渗流耦合模型(广义Biot模型)。计及介质中流固两相的惯性耦合，具体考虑模拟介质压力相关弹塑性本构行为的非关联Drucker-Prager准则，忽略液相和固相颗粒的压缩性。详细分析了二维情况下含液饱和多孔介质在动力荷载作用下波传播过程的间断和失稳，导出了在波传播过程中产生驻波间断和动力颤振失稳的临界条件。驻波间断是由于应变软化导致材料失稳的结果，它并不一定意味着完全丧失波通过间断面在介质中继续传播的能力。颤振失稳则是因为模拟含液多孔介质固体骨架的非关联塑性本构行为所致，它可以先于驻波间断、即在塑性硬化阶段发生；但它仅可能在含液多孔介质中发生，对于固体材料即使为非关联塑性连续体也不可能发生颤振失稳现象。 波的逸散性意味着波的相速度随频率而变化。这一性质与正确模拟波在因应变软化引起的应变局部化区域中的传播密切相关。本论文基于上述耦合模型，对单轴应变一维情况讨论了非线性饱和一非饱和多孔介质中波传播过程的失稳现象和逸散性。分析了流固粘性耦合、流固惯性耦合、流固混合体的压缩性、孔隙饱和度及固体骨架材料在高应变速率下粘弹塑性本构行为等因素对失稳与逸散性的影响。由此所获得的结果和结论将为克服含液多孔介质在强动荷载作用下波传播过程数值模拟的困难提供理论基础和线索。 实验观察表明，在粘性土等多孔介质材料中因应变软化引起的在局部区域发生并急剧发展塑性变形的剪切带具有一定的宽度。此外，剪切带的萌发、发展直至最终形成是一个渐进破坏过程。为数值模拟和再现这一渐进破坏过程，本论文工作中作为正则化机制引入梯度塑性模型。对饱和多孔介质(也能作为退化情况用于固体材料)提出了一个归结为线性互补问题(LinearComplementaryProblem)求解过程的梯度塑性连续体有限元—无网格耦合方法。 为模拟材料的弹塑性本构行为，对固体材料和饱和多孔介质分别采用von-Mises准则和非关联Drucker-Prager准则。利用在积分点上定义的离散塑性乘子值和采用基于移动最小二乘(MovingLeast-Square)的无网格法插值近似假定塑性乘子场。而位移和压力场则利用定义在节点上的离散值采用有限元法插值近似。因而可充分发挥无网格法与有限元法的各自优势，而避免它们的各自缺点。 通过建立平衡方程的弱形式实现空间离散化，结合在积分点上逐点满足而不是积分意义下满足非局部本构方程和屈服条件，导出相应的线性互补问题标准型。并通过Lexico-Lemke算法求解。发展了一个基于向后欧拉返回映射积分方案和利用Newton-Raphson方法的全局迭代过程的一致性算法，使得空间离散的平衡方程和在每个积分点上的非局部本构方程和屈服准则在每次全局迭代中同时满足。值得强调指出，所提出方法在保证二阶收敛率的同时无需形成非局部一致性切线刚度矩阵；另外，对于非关联塑性模型，所导出的为线性互补问题求解的全局广义刚度阵仍保持对称。数值结果表明，所发展的模型和一致性算法能正确模拟由应变软化引起的以应变局部化为特征的渐进破坏过程。

目录

文摘

英文文摘

声明

第一章绪论

1.1课题的工程背景和理论意义

1.2土力学的发展历史

1.3多孔介质应变局部化问题的研究现状

1.3.1实验观测

1.3.2理论分析

1.3.3数值模拟

1.4无网格方法简介

1.4.1无网格法的发展历程

1.4.2无网格法评价

1.5本文的主要工作

参考文献

第二章 多孔介质的基本概念与数学模型

2.1多孔介质概述

2.2多孔介质的结构及其构造

2.3多孔介质的基本物理指标

2.4多孔介质的应力

2.4.1饱和多孔介质的有效应力

2.4.2非饱和多孔介质的有效应力

2.4.3非饱和多孔介质的净应力

2.5多孔介质中液相和气相的流动

2.5.1 Darcy定律

2.5.2渗透率

2.6多孔介质的数学模型

2.6.1饱和多孔介质的数学模型

2.6.2饱和多孔介质固结问题的数学模型(u-p模型)

2.6.3非饱和多孔介质的数学模型

参考文献

第三章经典塑性与梯度塑性模型

3.1应力不变量和应力空间

3.2经典塑性模型

3.2.1 Tresca屈服条件

3.2.2 vorl-Mises屈服条件

3.2.3 Mohr-Coulomb屈服条件

3.2.4 Drucker-Prager屈服条件

3.3梯度塑性模型

3.3.1梯度塑性屈服函数

3.3.2内部长度参数

3.4本构关系

3.4.1流动法则

3.4.2等效塑性应变

3.4.3向后欧拉积分算法

3.4.4一致性本构方程

参考文献

第四章含液多孔介质的失稳与逸散性分析

4.1引言

4.2含液体多孔介质的数学模型与控制方程

4.3固体骨架的非线性本构方程

4.4.波传播过程中的失稳和逸散性

4.4.1率无关弹塑性介质中考虑渗透粘性的波传播过程分析

4.4.2率相关弹塑性多孔介质中的波传播过程

4.5小结

参考文献

第五章饱和多孔介质的驻波间断和颤振失稳

5.1引言

5.2饱和多孔介质的控制方程

5.3间断面上的运动学协调条件和动量平衡条件

5.4波速特征方程

5.5波传播过程的间断面上的驻波失稳和颤振失稳

5.5.1波传播过程的间断性和驻波失稳

5.5.2惯性耦合对波速的影响

5.5.3波传播过程的颤振失稳

5.6小结

参考文献

第六章梯度塑性连续体的有限元一无网格耦合方法

6.1引言

6.2无网格方法

6.2.1移动最小二乘法基本原理

6.2.2权函数和影响域

6.2.3形函数的导数及其改进

6.2.4.本质边界条件的施加

6.2.5区域积分的实现

6.3梯度塑性连续体的控制方程

6.4控制方程的空间离散

6.5线性互补问题标准型的推导及其一致性迭代算法

6.5.1迭代求解方案的基本方程

6.5.2线性互补问题的迭代求解形式

6.6数值算例

6.7小结

参考文献

第七章饱和多孔介质梯度塑性模型的有限元一无网格耦合方法

7.1引言

7.2饱和多孔介质控制方程和梯度塑性本构模型

7.3控制方程的离散化

7.3.1空间离散

7.3.2时间离散

7.4线性互补问题标准型的推导及其一致性迭代算法

7.4.1迭代求解方案的基本方程

7.4.2线性互补问题及其迭代求解形式

7.5数值算例

7.6小结

参考文献

第八章程序实现

8.1数据结构说明

8.2程序模块说明

8.3程序流程

第九章总结与展望

9.1总结

9.2展望

论文创新点摘要

攻读博士学位期间发表学术论文情况

致谢