基于双应力原理的饱和岩体本构模型构建与有限元程序开发研究.docxVIP

基于双应力原理的饱和岩体本构模型构建与有限元程序开发研究

一、绪论

1.1研究背景与意义

在各类大型工程建设中，如水利水电工程中的大坝建设、地下厂房开挖，交通工程中的隧道修建，以及采矿工程中的矿井开采等，饱和岩体作为重要的工程介质，其力学行为对工程的稳定性和安全性起着决定性作用。例如，在大坝基础中，饱和岩体的变形和强度特性直接影响大坝的整体稳定性；在隧道工程中，饱和岩体的力学响应关系到隧道的支护设计和施工安全。因此，深入研究饱和岩体的力学特性具有至关重要的工程意义。

双应力原理作为描述饱和岩体力学行为的重要理论，考虑了有效应力和孔隙压力对岩体力学特性的综合影响，能够更准确地揭示饱和岩体在复杂应力状态下的力学响应机制。基于双应力原理构建饱和岩体本构模型，能够为工程设计和分析提供更符合实际的理论依据，提高工程的可靠性和安全性。

有限元程序作为一种强大的数值模拟工具，在岩体工程分析中具有广泛的应用前景。通过开发基于双应力原理的饱和岩体本构模型的有限元程序，可以对复杂的岩体工程问题进行数值模拟分析，预测饱和岩体在不同工况下的力学行为，为工程决策提供科学依据，有效减少工程风险和成本。因此，开展基于双应力原理的饱和岩体本构模型及有限元程序开发的研究，具有重要的理论意义和工程实用价值。

1.2研究现状

1.2.1饱和岩体本构模型研究进展

饱和岩体本构模型的发展经历了多个阶段，从早期的线弹性模型到后来的弹塑性模型、粘弹性模型、损伤模型等，不断向更能准确描述饱和岩体复杂力学行为的方向发展。线弹性模型如胡克定律，简单地假设岩体的应力与应变成正比，虽然计算简便，但无法考虑岩体的非线性特性和塑性变形，在实际工程应用中具有很大的局限性。随着对岩体力学特性认识的深入，弹塑性模型逐渐得到发展，如Mohr-Coulomb模型和Drucker-Prager模型。Mohr-Coulomb模型基于Mohr强度理论，考虑了岩体的抗剪强度和摩擦特性，在岩土工程中得到了广泛应用；Drucker-Prager模型则是在Mohr-Coulomb模型的基础上进行了改进，考虑了中间主应力的影响，更适合描述岩土材料的屈服和破坏行为。然而，这些早期的弹塑性模型仍存在一些不足之处，如无法准确描述岩体的剪胀性和应变软化特性。

为了克服上述问题，学者们进一步提出了考虑剪胀性和应变软化的弹塑性模型。这些模型通过引入相关的参数和理论，能够更准确地描述岩体在加载和卸载过程中的力学行为。例如，一些模型通过引入剪胀角来考虑岩体的剪胀特性，通过定义软化函数来描述岩体的应变软化现象。此外，粘弹性模型和损伤模型也在饱和岩体本构模型研究中得到了广泛关注。粘弹性模型考虑了岩体的时间效应，能够描述岩体在长期荷载作用下的蠕变和松弛现象；损伤模型则考虑了岩体内部微裂纹和孔隙的发展对其力学性能的影响，通过引入损伤变量来描述岩体的损伤演化过程。

尽管目前已经提出了多种饱和岩体本构模型，但这些模型在描述饱和岩体力学行为时仍存在一些局限性。一方面，部分模型过于简化，无法全面考虑饱和岩体的复杂特性，如多相介质的相互作用、各向异性等；另一方面，一些模型虽然能够考虑更多的因素，但模型参数众多，确定困难，且计算过程复杂，不利于工程实际应用。因此，开发一种既能准确描述饱和岩体力学行为，又具有简单实用特点的本构模型，仍然是当前饱和岩体力学研究的重要任务之一。

1.2.2有限元程序在岩体工程中的应用现状

有限元程序在岩体工程领域得到了广泛应用，涵盖了岩土边坡稳定性分析、地下洞室开挖模拟、地基基础设计等多个方面。在岩土边坡稳定性分析中，有限元程序通过模拟边坡在自重、地下水、地震等荷载作用下的应力应变分布，预测边坡的潜在滑动面和安全系数，为边坡的加固设计提供依据。在地下洞室开挖模拟方面，有限元程序可以考虑洞室开挖过程中的卸荷效应、支护结构与岩体的相互作用等因素，分析洞室围岩的变形和稳定性，优化支护方案。在地基基础设计中，有限元程序能够模拟地基在建筑物荷载作用下的沉降和承载力，评估地基的稳定性，指导基础设计。

众多学者和工程师利用有限元程序对各类岩体工程问题进行了深入研究，并取得了丰硕的成果。例如，在某大型水利枢纽工程的地下厂房开挖分析中，通过有限元模拟准确预测了围岩的变形和应力分布，为支护设计提供了科学依据，确保了工程的安全施工。在某高速公路边坡稳定性分析中，利用有限元程序对比了不同加固方案下边坡的稳定性，选择了最优的加固措施，保障了公路的运营安全。

然而，有限元程序在岩体工程应用中也存在一些问题。首先，岩体材料的本构模型复杂多样，不同的本构模型对计算结果的影响较大，如何选择合适的本构模型仍是一个难题。其次，岩体工程中的边界条件和初始条件难以准确确定，这也会给计算结果带来一定的误差。此外，有限元计算需要大量的计