考虑惯性耦合的两相介质动力反应显式数值计算方法∗

基于包含惯性耦合（耦合质量密度系数）的 u?p 形式饱和两相介质弹性波动方程,开展了对饱和两相介质近场波动问题时域全显式数值计算方法的研究。通过对波动方程中的质量集中矩阵和流体压缩矩阵对角化处理,消除了方程中的动力耦联的影响。分别应用中心差分法和 Newmark 常平均加速度法求解固相位移和固相速度,基于向后差分法求解孔隙流体压力,建立了考虑惯性耦合的饱和两相介质的时域显式数值计算方法。
将该方法的数值解与相应的解析解进行对比,二者符合良好,验证了该方法的正确性。将该方法应用于饱和两相介质自由场的地震反应问题,研究了惯性耦合对两相介质动力反应的影响。基于本文研究的两种情况表明,惯性耦合对饱和两相介质动力反应的影响程度有限。