基于传递矩阵的两相介质波动问题时域显式有限元方法稳定-岩土力学.PDF

第 35 卷第 6 期 岩 土 力 学 Vol.35 No.6 2014 年 6 月 Rock and Soil Mechanics Jun. 2014 文章编号：1000－7598 －(2014) 06 －1768－07 基于传递矩阵的两相介质波动问题时域显式 有限元方法稳定性研究 李 亮，杜修力 （北京工业大学 城市与工程安全减灾教育部重点实验室，北京 100124） 摘 要：饱和两相介质波动问题的时域显式有限元方法具有条件稳定性。这类方法的稳定性可通过方法中动力反应时域递推 计算格式传递矩阵传递因子的数值进行评估，传递因子的数值越小，方法的稳定性越好。针对饱和两相介质波动问题的两种 时域显式有限元方法，基于其动力反应时域递推计算格式的传递矩阵的性质，在两相介质体系物理参数取值不同的条件下， 对两种方法稳定性性质的差异进行了对比研究。结果表明：渗透系数的取值对两种时域显式有限元方法稳定性性质的差异具 有较为显著的影响。当渗透系数取值较大时，两种方法的稳定性性质没有明显的差异；当渗透系数的取值较小时，两种方法 的稳定性性质将表现出明显的差异，方法 1 的稳定性要优于方法2 。 关 键 词：饱和两相介质；波动问题；时域显式有限元方法；稳定性；传递矩阵；渗透系数 中图分类号：TU 435 文献标识码：A Study of stability of explicit finite element method for wave motion of two-phase medium based on transfer matrix LI Liang ，DU Xiu-li (The Key Laboratory of Urban Security and Disaster Engineering of Ministry of Education, Beijing University of Technology, Beijing 100124, China) Abstract: The time-domain explicit finite element method for the wave motion of saturated two-phase medium is conditionally stable, and the stability of this kind of method can be evaluated by the transfer factor of transfer matrix of calculating formula. The smaller the value of transfer factor is, the better the stability of