如何模拟膨胀土区域的非饱和瞬变渗流与变形的耦合相互作用？

如何模拟膨胀土区域的非饱和瞬变渗流与变形的耦合相互作用？

How to simulate the interaction between variably saturated transient flow and the deformation of expansive soils？

问题描述：

膨胀土是国际上普遍认为难对付的一种含水量变化可引起胀缩变形的特殊非饱和黏土。因膨胀土胀缩变形带来铁路、公路、基坑、建筑物一系列的工程问题或事故。膨胀土因为其微观结构随着含水量而变化，力学行为比一般非饱和黏性土更复杂。膨胀土因失水收缩其表面一般具有裂隙性特征。降雨入渗一般是造成膨胀土工程问题的诱因，入渗过程中膨胀土中的水分迁移的模拟是一个难题。水分迁移影响膨胀土体的强度和变形特性参数。膨胀土区域的工程力学分析是复杂非饱和瞬变渗流与变形强耦合作用问题，其合理模拟与评估是工程合理设计、防灾减灾的关键。工程问题的模拟需要功能强大的软件来完成。在已有自主开发的具备渗流与变形耦合模拟、土与结构相互作用模拟能力的大型岩土工程有限元分析的基础上，发展膨胀土工程的模拟方法，开发满足膨胀土区域边坡、路基、基坑和建筑地基等工程分析需求的软件，一方面可直接满足重大工程的研究需求，另一方面通过对典型的膨胀土工程的研究，归纳总结出简单适用的膨胀土地区水分迁移、作用在结构上的膨胀力的计算评估方法，成为膨胀土工程区域合理设计、科学防灾减灾中亟待解决的问题。

问题背景：

膨胀土是一种含水量变化可引起胀缩变形的国际上公认难以对付的特殊黏土。美国因膨胀土工程问题造成的损失年均在20亿美元以上，已超过洪水、飓风、地震和龙卷风所造成的经济损失的总和。我国膨胀土分布面积超过10万平方公里，因膨胀土问题导致的铁路与公路路基与边坡工程、基坑工程、基础工程事故极为常见，损失情况近年来比美国更为严重。随着我国交通工程和城市开发的高速发展，膨胀土边坡、基坑工程因降雨诱发的破坏损失越来越大。

然而，学术界大多成果集中在机理和理论模型研究方面、缺乏对复杂工程的模拟方法和膨胀土宏观力学效果的研究，尚未有膨胀土胀缩变形对支护结构的膨胀力计算的实用方法；工程界通常仅采取要求勘察提供膨胀力或对强度参数进行折减，其他方面依然按照普通黏土工程来分析。随着我国公路、铁路和城市建设的高速发展，发生事故的概率越来越高。

最新进展：

膨胀土的基本特性及其描述方法的研究已经从定性向定量发展。膨胀土的胀缩性、裂隙性、超固结性是土体含水量变化及其循环变化的外在表现。胀缩的机理、主要影响因素、理论模型的研究相对较充分。膨胀土中裂隙发展的规律及其描述方法、含裂隙膨胀土的渗流模拟方法、膨胀土中的侧压力变化规律还处于定性研究阶段，其定量模拟方法有待于进一步发展。膨胀土的应力-应变模型，大多为基于特定的试验约束条件的经验模型，反映含水量变化导致变形的本构模型有待于进一步发展。

膨胀土工程的力学分析问题，是含裂隙多孔介质非饱和瞬变渗流与应力变形耦合问题，常常还包含土与结构的相互作用。将水分迁移导致裂隙的产生和胀缩变形对力学行为的影响，加入到有限元数值模拟方法和软件中，便能实现膨胀土工程的定量力学模拟。裂缝产生形成的孔隙-裂隙介质，是耦合模拟中对渗流模拟方法的新挑战，而胀缩变形的模拟，则可采用初应变法来实现。提高大规模耦合计算的收敛性与迭代效率，是膨胀土工程数值模拟问题的关键点。基于自主开发的软件系统发展和实现膨胀土工程模拟方法，是科学解决我国量大面广的膨胀土工程防灾减灾的关键。

重要意义：

突破膨胀土区域非饱和瞬变渗流与变形耦合相互作用模拟问题，实现复杂条件下膨胀土工程中水分迁移和应力变形的定量模拟，一方面可以为膨胀土区域的重大工程的设计方案研究提供直接的计算工具，另一方面可逐步取得对典型工程中水分迁移和膨胀土对结构作用力的规律性认识，并进一步提炼出简单实用的结构上膨胀力的计算公式，便于膨胀土工程设计中的推广应用。将对膨胀土工程的合理设计、膨胀土灾害的防治提供更加科学的方法。模拟方法和软件的发展，也可为我国水利水电、采矿、交通、建筑等设计科研单位的工程设计和研究提供可持续改进的自主可控分析软件。