基于薄层-有限元耦合的饱和两相介质波动问题深度解析与应用.docxVIP

基于薄层-有限元耦合的饱和两相介质波动问题深度解析与应用

一、引言

1.1研究背景与意义

饱和两相介质广泛存在于自然界和工程实际中，如海洋底部的饱和沉积土层、地下含水层以及一些建筑工程中的地基等。在这些领域中，饱和两相介质的波动问题对于工程的安全与稳定具有至关重要的影响。例如，在地震发生时，饱和土体中的地震波传播特性会直接影响到地面建筑物和地下结构的动力响应，进而关系到这些结构的安全性。如果不能准确地分析饱和两相介质中的波动问题，就难以对工程结构在地震等动力荷载作用下的性能进行有效评估，可能导致工程结构在灾害发生时遭受严重破坏，威胁人民生命财产安全。

传统的单相介质理论在描述饱和两相介质的力学行为时存在局限性，因为它无法考虑孔隙流体与固相骨架之间的相互作用。而实际的饱和两相介质中，孔隙流体的存在会显著影响介质的力学响应。例如，在饱和土体中，孔隙水的流动会产生渗流力，这种力会改变土体的有效应力分布，进而影响土体的变形和强度特性。因此，需要一种更准确的理论和方法来研究饱和两相介质的波动问题。

薄层-有限元耦合分析方法结合了薄层法和有限元法的优点，为解决饱和两相介质波动问题提供了新的途径。薄层法在处理层状介质问题时具有高效性和准确性，能够快速得到层状介质的动力响应特性。有限元法则具有强大的适应性，能够处理复杂的几何形状和边界条件。将这两种方法耦合起来，可以充分发挥它们的优势，更精确地模拟饱和两相介质的波动特性。例如，在分析饱和土-桩体系的动力相互作用时，利用薄层法可以快速计算出饱和土层的动力阻抗，再结合有限元法对桩体进行精细建模，能够准确地得到桩体在饱和土体中的动力响应。这种耦合分析方法在岩土工程、地震工程等领域具有广阔的应用前景，有望为相关工程的设计、施工和安全评估提供更可靠的理论支持和技术手段。

1.2国内外研究现状

国外学者在饱和两相介质波动理论和数值方法方面开展了大量的研究工作。1956年，Biot建立了流体饱和两相多孔介质波动方程，为该领域的研究奠定了理论基础。Biot假定两相介质具有孔隙相互连通且孔隙中的流体具有可压缩性的特性，两相介质具有统计各向同性，固相骨架具有线弹性的本构关系，由此建立了通过惯性项和粘性项耦合起来的弹性波动方程。此后，众多学者基于Biot理论进行了深入研究。Ishihara将Biot理论中两相介质的压缩性和方程中的系数关联起来，推导出饱和介质波动方程的近似形式。1980年，Zienkiewicz等首先提出u-p形式的波动方程，即以固相位移u和孔隙水压力p为基本变量，它适用于如地震荷载一类的情况，为数值求解提供了更便利的形式。1991年，苏联学者Nikolaevsky利用热力学第一定律建立起了两相多孔弹性介质的动力方程，该方法采用的β1具备可测定的量，因此相对于Biot的方程更加方便，从新的理论角度丰富了饱和两相介质波动理论。

在数值模拟方面，有限差分方法是最早用于处理饱和两相介质波动问题的数值方法之一。1974年，Garg等利用有限差分方法对一维弹性或非弹性饱和多孔介质中波的传播问题进行了研究，开启了数值求解该问题的先河。Biot动力方程的第一个时域有限元数值解法是1972年由Ghaboussi和Wilson提出的，该方法在Biot动力方程的基础上建立变分公式从而得到有限元方程，并通过有限元方程分析大坝的地震瞬态反应及饱和两相介质的动荷载作用下的瞬态反应，推动了有限元法在该领域的应用。Prevost提出了一种对时间步长的积分采用隐—显式算法的有限元分析方法，用于分析饱和两相介质的动力响应问题，进一步改进了有限元求解算法。1992年，Hirai建立了一种分析饱和多孔弹性层状半空间在Rayleigh波作用下的土层反应的有限元数值方法，拓展了有限元法在特定波动问题中的应用。1997年和1999年，Khalili等人先后建立了针对频域内分析两相介质无限域中波传播问题的一维和二维无限元方法，为解决无限域波动问题提供了新的思路。

国内学者在该领域也取得了丰硕的研究成果。苗天德等在1995年根据混合物理论，讨论了饱和多孔介质波动的一般方程，从混合物理论的角度完善了饱和两相介质波动理论体系。章根德在1993年介绍了流体、固体混合物连续介质理论及其应用，推动了相关理论在国内的传播与应用。崔杰在1987年及1995年利用有限元方法分别建立了分析地震波作用下一维及二维饱和土动力响应的模型，为国内有限元法求解饱和土波动问题提供了重要的研究基础。赵成刚等于2003年在国际上率先采用饱和多孔介质Biot动力理论模拟局部不均匀场地的地震波散射问题，并得到了多个解析解，还最早建立对两相饱和多孔介质波动