两相介质波动问题显式有限元方法稳定性的研究.pdf

第十一届土力学及岩土工程学术会议 中国·兰州 2011 年 8 月 16—19 日 两相介质波动问题显式有限元方法稳定性研究 李 亮，杜修力，李立云，翟 威 （北京工业大学 城市与工程安全减灾教育部重点实验室，北京市朝阳区平乐园 100 号 100124 ） 摘 要：针对饱和两相介质波动问题的时域显式有限元方法，基于动力反应递推计算格式的传递矩阵的性质，进行了方法稳 定性问题的研究。定义了综合考虑各种影响因素的稳定性判别指标—传递因子，用以表征方法的稳定性。研究了稳定性影响 因素，包括时间步长、空间步距和渗透系数取值的作用规律，并给出了方法稳定性的实用判别准则。研究结果表明：时间步 长、空间步距和渗透系数取值都对饱和两相介质波动问题的时域显式有限元方法的稳定性有较为显著的影响。随着时间步长 的增大，方法的稳定性降低；随着空间步距的增大，方法的稳定性增加；当渗透系数取值增大时，方法的稳定性增加。 关键词：饱和两相介质；波动问题；时域显式有限元方法；稳定性；传递矩阵 中图分类号：TU435 文献标识码：A 文章编号： Study on Stability of Explicit Finite Element Method for Wave Motion of Fluid-saturated Porous Media LI Liang, DU Xiu-li, LI Li-yun, ZHAI Wei (The Key Laboratory of Urban Security and Disaster Engineering ( Beijing University of Technology ), Ministry of Education, Beijing 100124, China) Abstract: In this paper, the stability of explicit finite element method for wave motion of fluid-saturated porous media is studied based on the transfer matrix of calculating formula of the method. Transfer factor is defined to describe the stability of the method. The effect of time step, space interval and the value of permeability coefficient on stability is analyzed and practical stability criterionis is put forward. Calculating results show that with the increase of time step, the stability of the method will be deteriorated; with the increase of space interval, the stability of the method will be enhanced; with the increase of the value of permeability coefficient, the stability of the method will be enhanced. Key words: Fluid-saturated porous media; Wave motion; Time-domain explicit finite element method; Stability; Transfer matrix 对于少数边界条件特殊，荷载形式简