横观各向同性饱和两相介质波动问题的时域显式有限元方法 被引量：2

1

作者 李亮 翟威 杜修力 《岩土工程学报》 EI CAS CSCD 北大核心 2012年第3期464-470,共7页

针对横观各向同性饱和两相介质的弹性波动问题,分别应用空间解耦技术和时域显式逐步积分计算格式完成对该问题的有限元空间离散和时间离散,建立了横观各向同性饱和两相介质波动问题计算分析的时域显式有限元方法。应用该方法进行了输入... 针对横观各向同性饱和两相介质的弹性波动问题,分别应用空间解耦技术和时域显式逐步积分计算格式完成对该问题的有限元空间离散和时间离散,建立了横观各向同性饱和两相介质波动问题计算分析的时域显式有限元方法。应用该方法进行了输入地震波作用下横观各向同性饱和两相介质动力反应的计算分析,并将计算结果与完全各向同性饱和两相介质的计算结果进行对比研究,并对各向异性系数取值对于横观各向同性饱和两相介质动力反应计算结果的影响进行了研究和讨论。计算结果表明,横观各向同性饱和两相介质与完全各向同性饱和两相介质的动力反应特性具有较为显著的差异,各向异性系数取值对于横观各向同性饱和两相介质动力反应计算结果具有较为重要的影响。本文的数值计算工作同时表明,本文建立的时域显式有限元方法是进行横观各向同性饱和两相介质动力反应计算分析的一种有效方法。 展开更多

关键词 横观各向同性 饱和两相介质:弹性波动 时域显式有限元法

在线阅读 下载PDF 职称材料

基于精细时程积分的饱和两相介质波动问题时域解法 被引量：2

2

作者 段雪铭 李亮 +1 位作者 杜修力 宋佳 《岩土力学》 EI CAS CSCD 北大核心 2015年第9期2702-2707,2715,共7页

针对u-p形式的饱和两相介质的波动方程,采用精细时程积分方法计算固相位移u,采用向后差分算法求解流体压力p,建立了基于精细时程积分技术的饱和两相介质波动问题的时域求解方法。针对标准算例,将该方法的计算结果与可视为标准结果的Zien... 针对u-p形式的饱和两相介质的波动方程,采用精细时程积分方法计算固相位移u,采用向后差分算法求解流体压力p,建立了基于精细时程积分技术的饱和两相介质波动问题的时域求解方法。针对标准算例,将该方法的计算结果与可视为标准结果的Zienkiewicz隐式算法的计算结果进行比较分析,二者符合较好,表明了该方法具有良好的计算精度。同时,该方法的计算过程为交替迭代求解,避免了在每个时间分析步上求解耦联方程组,因而具有较高的计算效率。该方法具有时域显式计算方法的基本特点,是进行饱和两相介质动力问题计算与分析的一种有效方法。 展开更多

关键词 饱和两相介质 u-p形式波动方程 精细时程积分 时域解法 向后差分

在线阅读 下载PDF 职称材料

集中力作用下两相饱和介质二维位移场Green函数 被引量：12

3

作者 丁伯阳 丁翠红 +1 位作者 陈禹 陶海冰 《应用数学和力学》 EI CSCD 北大核心 2004年第8期869-874,共6页

　由于工程场地的对称性,集中力作用下的位移场Green函数在土力学、地震工程学和动力基础方面的应用需以二维模型出现·　在理论推导上Green函数的二维模型要比三维模型复杂·　根据丁伯阳等人已得到的三维位移场中集中力作用... 　由于工程场地的对称性,集中力作用下的位移场Green函数在土力学、地震工程学和动力基础方面的应用需以二维模型出现·　在理论推导上Green函数的二维模型要比三维模型复杂·　根据丁伯阳等人已得到的三维位移场中集中力作用下两相饱和介质位移场Green函数,采用DeHoop与Manolis给出的沿x3方向在无穷域积分方法,得到了集中力作用下两相饱和介质二维位移场Green函数·　相比已有的工作,所得结果不仅简单。 展开更多

关键词 集中力 两相饱和介质 二维Green函数

在线阅读 下载PDF 职称材料

集中力作用下的两相饱和介质位移场Green函数 被引量：13

4

作者 丁伯阳 丁翠红 孟凡丽 《力学学报》 EI CSCD 北大核心 2001年第2期234-241,共8页

以复模两相饱和介质Ｂｉｏｔ动力学方程为基础，根据该方程Ｄ’Ａｌｅｍｂｅｒｔ解的Ｆｏｕｒｉｅｒ变换所属的Ｈｏｍｈｏｌｔｚ方程特性，由Ｂｉｏｔ方程解的相容性条件及δ函数性质较好地解决了快、慢纵波位势的耦合问题，较为简便地... 以复模两相饱和介质Ｂｉｏｔ动力学方程为基础，根据该方程Ｄ’Ａｌｅｍｂｅｒｔ解的Ｆｏｕｒｉｅｒ变换所属的Ｈｏｍｈｏｌｔｚ方程特性，由Ｂｉｏｔ方程解的相容性条件及δ函数性质较好地解决了快、慢纵波位势的耦合问题，较为简便地得到了两相饱和介质在集中力作用下低频（ω＜ωｃ）时的频域和时域的Ｇｒｅｅｎ函数． 展开更多

关键词 集中力 两相饱和介质 GREEN函数 位移场 Homholtz方程 Δ函数

在线阅读 下载PDF 职称材料

流固耦合饱和两相介质动力模型在ABAQUS中的实现 被引量：8

5

作者 李亮 崔智谋 +1 位作者 康翠兰 王相宝 《岩土工程学报》 EI CAS CSCD 北大核心 2013年第S2期281-285,共5页

采用以固相位移和孔压为基本未知量的u–p形式的方程组对于研究饱和土体的地震动力响应问题相对简单。利用通用有限元软件ABAQUS所提供的用户自定义单元(UEL)接口,编制相应的接口程序,开发了基于流固耦合两相介质动力模型的孔压单元。... 采用以固相位移和孔压为基本未知量的u–p形式的方程组对于研究饱和土体的地震动力响应问题相对简单。利用通用有限元软件ABAQUS所提供的用户自定义单元(UEL)接口,编制相应的接口程序,开发了基于流固耦合两相介质动力模型的孔压单元。通过算例对所开发单元的正确性进行了验证,结果表明,所开发的孔压单元能够对饱和土体的动力响应进行较为准确的计算,而ABAQUS中提供的CEP4P孔压单元动力反应计算结果与实际情况有较为明显的差异。 展开更多

关键词 饱和两相介质 流固耦合动力模型 ABAQUS 用户自定义单元(UEL)

在线阅读 下载PDF 职称材料

流体饱和两相多孔介质动力反应计算分析 被引量：3

6

作者 李亮 杜修力 +1 位作者 赵成刚 李立云 《岩土力学》 EI CAS CSCD 北大核心 2008年第1期113-118,共6页

基于流体饱和两相多孔介质的弹性波动方程组,运用显式逐步积分格式与局部透射人工边界相结合的时域显式有限元方法对该波动方程组进行求解,对两相介质在输入地震波作用下的弹性动力反应进行计算和分析;对在是否考虑孔隙流体渗流的两种... 基于流体饱和两相多孔介质的弹性波动方程组,运用显式逐步积分格式与局部透射人工边界相结合的时域显式有限元方法对该波动方程组进行求解,对两相介质在输入地震波作用下的弹性动力反应进行计算和分析;对在是否考虑孔隙流体渗流的两种情况下计算得到的两相介质弹性动力反应结果的差异进行对比研究,从而揭示孔隙流体渗流对两相介质动力反应性质的影响。计算结果表明:两相介质弹性动力反应时程的波形与入射地震波的波形相同,且弹性动力反应的峰值出现的时刻对应于入射地震波的峰值出现的时刻;孔隙流体的渗流将对两相介质的弹性动力反应性质产生显著的影响。数值计算同时表明,时域显式有限元方法是进行流体饱和两相多孔介质弹性动力反应计算分析的一种有效的方法。 展开更多

关键词 流体饱和两相多孔介质 弹性 动力反应 时域显式有限元方法 渗流影响

在线阅读 下载PDF 职称材料

两相饱和介质层与单相介质层应力-位移函数的传递与退化 被引量：3

7

作者 丁伯阳 陈樟龙 徐庭 《应用数学和力学》 CSCD 北大核心 2014年第2期162-180,共19页

根据Biot孔隙弹性介质动力控制方程,利用快、慢纵波的解耦,求得满足两相饱和介质位移-应力传播的一阶微分方程组.该方程组及传递函数能退化到单相介质的位移-应力传播微分方程组.利用界面应力-位移连续条件,分析了位移-应力从两相饱和... 根据Biot孔隙弹性介质动力控制方程,利用快、慢纵波的解耦,求得满足两相饱和介质位移-应力传播的一阶微分方程组.该方程组及传递函数能退化到单相介质的位移-应力传播微分方程组.利用界面应力-位移连续条件,分析了位移-应力从两相饱和介质向单相介质传播,构建了界面过渡传递矩阵.使原有的6×6阶应力-位移传递矩阵过渡为4×6阶矩阵,能与单相介质的4×4阶应力-位移传递矩阵结合.最后,采用经典的波传播模型对比验算了结果,它们一致吻合. 展开更多

关键词 两相饱和介质 单相介质 多层介质 应力位移函数 过渡传递矩阵 退化

在线阅读 下载PDF 职称材料

两相介质饱和土三维非轴对称稳态响应分析 被引量：8

8

作者 张玉红 黄义 《应用力学学报》 CAS CSCD 北大核心 2002年第3期85-89,共5页

用积分变换方法分析饱和土三维非轴对称波动方程 ,用刚度矩阵方法分析层状饱和土三维非轴对称稳态响应。以数值算例对比分析了单相土介质与两相饱和土介质三维非轴对称稳态动力响应。分析中未引入势函数及传统的Helmholtz分解 ,推导过... 用积分变换方法分析饱和土三维非轴对称波动方程 ,用刚度矩阵方法分析层状饱和土三维非轴对称稳态响应。以数值算例对比分析了单相土介质与两相饱和土介质三维非轴对称稳态动力响应。分析中未引入势函数及传统的Helmholtz分解 ,推导过程及解答均非常简捷 。 展开更多

关键词 两相介质 饱和土 三维 非轴对称 稳态响应 积分变换

在线阅读 下载PDF 职称材料

u–p格式饱和两相介质动力问题的显–隐式交替算法 被引量：4

9

作者 杜修力 宋佳 李亮 《岩土工程学报》 EI CAS CSCD 北大核心 2016年第3期395-403,共9页

针对以土骨架位移u和孔隙流体压力p为基本未知量的饱和两相介质动力方程,提出基于中心差分法和精细时程积分法相结合的时域显–隐交替算法,分析饱和两相介质动力反应。该算法对孔压方程非齐次项的计算分别采用Gauss积分,梯形积分和Lagra... 针对以土骨架位移u和孔隙流体压力p为基本未知量的饱和两相介质动力方程,提出基于中心差分法和精细时程积分法相结合的时域显–隐交替算法,分析饱和两相介质动力反应。该算法对孔压方程非齐次项的计算分别采用Gauss积分,梯形积分和Lagrange二次插值函数近似位移的方法近似模拟;同时提出一种基于精细积分方法的三次样条插值函数法,进而回避了对孔压方程耦合速度项积分时涉及的矩阵求逆运算问题。计算结果与Zienkiewicz显–隐式结果对比,两者结果基本一致,证明了方法的有效性。 展开更多

关键词 饱和两相介质 u–p动力方程 显–隐式交替算法 精细时程积分法

在线阅读 下载PDF 职称材料

实验室尺度下气水两相裂缝型介质弹性波速度的数值模拟分析 被引量：2

10

作者 段茜 刘向君 《石油物探》 EI CSCD 北大核心 2017年第3期338-348,共11页

岩心的非均质结构使得采用实验手段研究弹性波速度与岩石力学参数之间的关系时面临较大困难,且超声波测试工作量大,耗资耗时。基于弹性波动理论和Wood孔隙流体模量模型,采用交错网格有限差分方法对气水两相裂缝型介质的弹性波场进行了... 岩心的非均质结构使得采用实验手段研究弹性波速度与岩石力学参数之间的关系时面临较大困难,且超声波测试工作量大,耗资耗时。基于弹性波动理论和Wood孔隙流体模量模型,采用交错网格有限差分方法对气水两相裂缝型介质的弹性波场进行了数值模拟。计算了岩样在不同裂缝分布位置、孔隙度、裂缝产状、裂缝密度和含水饱和度情况下的纵、横波速度及波速比,讨论了这些因素对波速的影响。模拟分析结果表明,该数值模拟计算方法简单、精度高,对利用岩石波速来反演岩石的孔隙结构有一定的指导作用。 展开更多

关键词 弹性波速度 波动理论 气水两相 介质 裂缝

在线阅读 下载PDF 职称材料

饱和两相介质动力固结问题时域求解的精细时程积分方法 被引量：1

11

作者 李亮 高超 吴利华 《计算力学学报》 CAS CSCD 北大核心 2017年第1期27-34,共8页

针对u-p形式的饱和两相介质波动方程,采用精细时程积分方法计算固相位移u,采用向后差分算法求解流体压力p,建立了饱和两相介质动力固结问题时域求解的精细时程积分方法。针对标准算例,对该方法的计算精度进行了校核。开展了该方法相关... 针对u-p形式的饱和两相介质波动方程,采用精细时程积分方法计算固相位移u,采用向后差分算法求解流体压力p,建立了饱和两相介质动力固结问题时域求解的精细时程积分方法。针对标准算例,对该方法的计算精度进行了校核。开展了该方法相关算法特性的研究,对采用不同数值积分方法计算非齐次波动方程特解项计算精度的差异进行了对比研究,并对采用不同积分点数目的高斯积分法计算特解项条件下计算精度的差异进行了对比研究。研究结果表明,(1)该方法具有良好的计算精度。(2)计算非齐次波动方程特解项的数值积分方法中,梯形积分法的计算精度最差,高斯积分法、辛普生积分法和科茨积分法都具有较好的计算精度。(3)增加高斯积分点数目对于提高计算精度的作用并不显著。 展开更多

关键词 饱和两相介质 动力固结 u-p波动方程 精细时程积分 向后差分

在线阅读 下载PDF 职称材料

基于传递矩阵的两相介质波动问题时域显式有限元方法稳定性研究

12

作者 李亮 杜修力 《岩土力学》 EI CAS CSCD 北大核心 2014年第6期1768-1774,共7页

饱和两相介质波动问题的时域显式有限元方法具有条件稳定性。这类方法的稳定性可通过方法中动力反应时域递推计算格式传递矩阵传递因子的数值进行评估,传递因子的数值越小,方法的稳定性越好。针对饱和两相介质波动问题的两种时域显式有... 饱和两相介质波动问题的时域显式有限元方法具有条件稳定性。这类方法的稳定性可通过方法中动力反应时域递推计算格式传递矩阵传递因子的数值进行评估,传递因子的数值越小,方法的稳定性越好。针对饱和两相介质波动问题的两种时域显式有限元方法,基于其动力反应时域递推计算格式的传递矩阵的性质,在两相介质体系物理参数取值不同的条件下,对两种方法稳定性性质的差异进行了对比研究。结果表明:渗透系数的取值对两种时域显式有限元方法稳定性性质的差异具有较为显著的影响。当渗透系数取值较大时,两种方法的稳定性性质没有明显的差异;当渗透系数的取值较小时,两种方法的稳定性性质将表现出明显的差异,方法1的稳定性要优于方法2。 展开更多

关键词 饱和两相介质 波动问题 时域显式有限元方法 稳定性 传递矩阵 渗透系数

在线阅读 下载PDF 职称材料

固体-液固两相饱和介质耦合动力有限元分析

13

作者 陈震 赵顺波 《人民黄河》 CAS 北大核心 2012年第4期138-139,142,共3页

基于固体介质和液固两相饱和介质动力学方程,采用有限元方法模拟脉冲荷载作用下介质耦合对结构动力响应的影响,并采用黏弹性人工边界模拟地震波的远域能量逸散效应。计算结果表明:介质耦合对地震波的传播有一定的影响;黏弹性人工边界可... 基于固体介质和液固两相饱和介质动力学方程,采用有限元方法模拟脉冲荷载作用下介质耦合对结构动力响应的影响,并采用黏弹性人工边界模拟地震波的远域能量逸散效应。计算结果表明:介质耦合对地震波的传播有一定的影响;黏弹性人工边界可以削弱地震波在固定边界处产生的反射现象,较固定边界计算结果更符合实际情况。 展开更多

关键词 固体介质 液固两相饱和介质 结构-地基动力相互作用 黏弹性人工边界 有限元分析

在线阅读 下载PDF 职称材料

两相介质中波在平面界面上的反射 被引量：2

14

作者 徐长节 吴世明 +1 位作者 蔡袁强 陈云敏 《爆炸与冲击》 EI CAS CSCD 北大核心 1998年第1期8-14,共7页

借助Ｂｉｏｔ多孔介质中的波动方程，根据各种界面条件导出了波在两相介质各种界面上反射的一般计算公式。作为例子，数值计算了Ｐ１波入射于几种典型界面时，Ｐ１、Ｐ２和Ｓ波的反射系数与频率、入射角之间的关系。结果表明：各种波的... 借助Ｂｉｏｔ多孔介质中的波动方程，根据各种界面条件导出了波在两相介质各种界面上反射的一般计算公式。作为例子，数值计算了Ｐ１波入射于几种典型界面时，Ｐ１、Ｐ２和Ｓ波的反射系数与频率、入射角之间的关系。结果表明：各种波的反射系数与入射角、频率有关系，若把多孔介质当作单相均一固体处理，将会得到显著不同的结果。 展开更多

关键词 两相介质 反射 弹性波 平面界面 P波 波动方程

在线阅读 下载PDF 职称材料

基于显式有限元方法的两相介质弹塑性动力反应计算分析 被引量：2

15

作者 李亮 杜修力 +1 位作者 赵成刚 翟威 《计算力学学报》 EI CAS CSCD 北大核心 2011年第3期423-429,共7页

针对增量形式的流体饱和两相多孔介质弹塑性波动方程组,运用基于显式逐步积分格式的时域显式有限元方法对该波动方程组进行求解,并应用基于SMP破坏准则的弹塑性动力本构模型描述两相介质的动力反应性质,对两相介质在输入地震波作用下的... 针对增量形式的流体饱和两相多孔介质弹塑性波动方程组,运用基于显式逐步积分格式的时域显式有限元方法对该波动方程组进行求解,并应用基于SMP破坏准则的弹塑性动力本构模型描述两相介质的动力反应性质,对两相介质在输入地震波作用下的弹塑性动力反应进行计算和分析,将计算结果与相应的弹性动力反应的计算结果进行对比;对本文应用的弹塑性动力本构模型进行参数研究,揭示模型参数的取值对两相介质弹塑性动力反应计算结果的影响。计算结果表明:两相介质弹塑性位移反应与相应的弹性位移反应具有较为显著的差别,表现为二者的峰值与反应时程的波形均有比较明显的差异;本文应用的弹塑性本构模型中的塑性势参数的取值对于两相介质弹塑性位移反应的计算结果具有显著的影响,而初始加载和卸载与重新加载阶段的硬化参数的比值仅影响动力反应结束后的塑性变形的幅值。本文工作表明,时域显式有限元方法是进行流体饱和两相多孔介质弹塑性动力反应计算分析的一种有效的方法。 展开更多

关键词 流体饱和两相多孔介质 弹塑性 动力反应 时域显式有限元方法 基于SMP破坏准则的动力本构模型

在线阅读 下载PDF 职称材料

双重孔隙介质中两相驱替理论及其计算 被引量：4

16

作者 郭大立 《西南石油大学学报（自然科学版）》 CAS CSCD 1993年第S1期90-92,共3页

1 基本方程一般多孔介质中两相驱替理论已由Buckley和Leverett建立,对于裂缝——孔隙介质中的同类问题,自60年代初提出了双重孔隙介质中两相渗流的数学模型以来。

关键词 孔隙介质 驱替 两相渗流 多孔介质 饱和度分布 裂缝系统 Buckley 同类问题 数学模型 相对渗透率曲线

在线阅读 下载PDF 职称材料

基于显式有限元方法的两相多孔介质地震反应研究

17

作者 李亮 杜修力 +1 位作者 赵成刚 李立云 《应用力学学报》 EI CAS CSCD 北大核心 2007年第4期550-553,共4页

基于流体饱和两相多孔介质的弹性波动方程组,运用显式逐步积分格式与局部透射人工边界相结合的时域显式有限元方法对该波动方程组进行求解,对两相多孔介质在输入地震波作用下的弹性动力反应进行了计算和分析,以揭示两相多孔介质弹性地... 基于流体饱和两相多孔介质的弹性波动方程组,运用显式逐步积分格式与局部透射人工边界相结合的时域显式有限元方法对该波动方程组进行求解,对两相多孔介质在输入地震波作用下的弹性动力反应进行了计算和分析,以揭示两相多孔介质弹性地震反应的规律和性质。计算结果表明:两相介质弹性地震反应时程的波形与入射地震波的波形相同,且弹性地震反应的峰值出现的时刻对应于入射地震波的峰值出现的时刻。本文的数值计算同时表明了时域显式有限元方法在进行两相多孔介质地震反应计算分析时的有效性。 展开更多

关键词 流体饱和两相多孔介质 弹性 地震反应 时域显式有限元方法

在线阅读 下载PDF 职称材料

横观各向同性含液饱和多孔介质黏弹性瞬态反应

18

作者 王俊林 祝彦知 张天航 《岩土力学》 EI CAS CSCD 北大核心 2007年第7期1315-1322,共8页

在考虑横观各向同性含液饱和多孔介质固体骨架和流体可压缩性以及固体骨架的黏弹性特征下,基于横观各向同性含液饱和多孔介质u-w形式的三维动力控制方程,以固相位移u、液相相对位移w为基本未知量,综合运用Laplace变换、双重Fourier变换... 在考虑横观各向同性含液饱和多孔介质固体骨架和流体可压缩性以及固体骨架的黏弹性特征下,基于横观各向同性含液饱和多孔介质u-w形式的三维动力控制方程,以固相位移u、液相相对位移w为基本未知量,综合运用Laplace变换、双重Fourier变换等方法,在直角坐标系下通过引入中间变量,将六元2阶动力控制方程组化为两组各含4个未知变量的常微分方程组,给出了直角坐标系下横观各向同性含液饱和多孔介质三维黏弹性动力反应的积分形式一般解;作为理论推导的验证,通过引入初始条件和边界条件,对横观各向同性含液饱和多孔介质半空间黏弹性瞬态反应问题进行了求解。解答的退化验证表明,所推导的理论解是正确的。 展开更多

关键词 BIOT波动方程 横观各向同性 含液饱和介质 黏弹性 动力反应

在线阅读 下载PDF 职称材料

孔隙-裂隙双重变形介质中的油水两相流动理论研究

19

作者 马利成 唐汝众 +1 位作者 仝兴华 薛世峰 《合肥工业大学学报（自然科学版）》 CAS CSCD 北大核心 2007年第1期95-99,共5页

在连续介质力学理论架构下,文章建立了饱和油水两相渗流与变形孔隙-裂隙双重介质耦合作用的理论来模拟裂隙砂岩储层中的油-水渗流,提出了基于双孔隙固相系统变形与油水两相渗流的全耦合力学模型;在对力学模型给予适当简化的基础上,应用... 在连续介质力学理论架构下,文章建立了饱和油水两相渗流与变形孔隙-裂隙双重介质耦合作用的理论来模拟裂隙砂岩储层中的油-水渗流,提出了基于双孔隙固相系统变形与油水两相渗流的全耦合力学模型;在对力学模型给予适当简化的基础上,应用解耦的有限元数值解方法对裂隙储层的两相流体压力、饱和度以及储层变形等参数进行了数值模拟;结果表明,裂隙储层流体运动主要由裂隙变形、油/水的PVT特性所决定。 展开更多

关键词 微裂缝性岩层 孔隙-裂隙双重介质 饱和油水两相渗流 耦合理论模型

在线阅读 下载PDF 职称材料

考虑惯性耦合的两相介质动力反应显式数值计算方法 被引量：1

20

作者 徐静茹 李亮 宋佳 《防灾减灾工程学报》 CSCD 北大核心 2023年第5期1009-1015,共7页

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

关键词 饱和两相介质 动力反应 时域全显式数值算法 惯性耦合 u‑p形式波动方程

在线阅读 下载PDF 职称材料