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

英文篇名：Study on Stability of Explicit Finite Element Method for Wave Motion of Fluid-saturated Porous Media
中文刊名：西北地震学报
英文刊名：Northwestern Seismological Journal
作者：李亮 ; 杜修力 ; 李立云 ; 翟威
英文作者：LI Liang ; DU Xiu-li ; LI Li-yun ; ZHAI Wei(The Key Laboratory of Urban Security and Disaster Engineerin of Ministry of Education ; Beijing University of Technology ; Beijing 100124 ; China)
中文关键词：饱和两相介质 ; 波动问题 ; 时域显式有限元方法 ; 稳定性 ; 传递矩阵 ; 判别准则
英文关键词：Fluid-saturated porous media ; Wave motion ; Time-domain explicit finite element method ; Stability ; Transfer matrix ; Discriminant criterion
出版日期：2011-09-30
机构：北京工业大学城市与工程安全减灾教育部重点实验室;
年：2011
期：03
出版单位：西北地震学报

摘要

基于动力反应递推计算格式传递矩阵的性质,进行了饱和两相介质波动问题时域显式有限元方法稳定性问题的研究。定义了综合考虑各种影响因素的稳定性判别指标——传递因子;研究了稳定性影响因素,包括时间步长、空间步距和渗透系数取值的作用规律;给出了方法稳定性的实用判别准则。研究结果表明:时间步长、空间步距和渗透系数取值都对饱和两相介质波动问题的时域显式有限元方法的稳定性有较为显著的影响。随着时间步长的增大,方法的稳定性降低;随着空间步距的增大,方法的稳定性增加;当渗透系数取值增大时,方法的稳定性增加。
The stability of time-domain explicit finite element method for wave motion of fluid-saturated porous media is studied based on the transfer matrix of calculating formula of dynamic response.A new discriminant index-transfer factor,which considers different influence factors synthetically,is defined to describe the stability of the method.The effects of time step,space interval and the value of permeability coefficient on stability are analyzed and the practical stability criterion is put forward.Calculating results show that with increase of time step,the stability of the method is reduced,with the increase of space interval,the stability is enhanced,and with the increase of the value of permeability coefficient,the stability of the method will be enhanced also.

引文

[1]Biot M A.Theory of propagation of elastic waves in a fluid-saturated porous solid[J].The Journal of the Acoustical Socie-ty of America,1956,28(2):168-191.
    [2]Zienkiewicz O C,Shiomi T.Dynamic behaviour of saturatedporous media;the generalized Biot formulation and its numeri-cal solution[J].International Journal for Numerical and Ana-lytical Methods in Geomechanics,1984,8(1):71-96.
    [3]Prevost J H.Wave propagation in fluid-saturated porousmedia:an efficient finite element procedure[J].Soil Dynamicsand Earthquake Engineering,1985,4(4):183-202.
    [4]赵成刚,王进廷,史培新等.流体饱和两相多孔介质动力反应分析的显式有限元法[J].岩土工程学报,2001,23(2):178-182.
    [5]王进廷,杜修力,赵成刚.液固两相饱和介质动力分析的一种显式有限元法[J].岩石力学与工程学报,2002,21(8):1199-1204.
    [6]廖振鹏,刘晶波.离散网格中的弹性波动[J].地震工程与工程振动,1986,6(2):1-16.
    [7]王进廷,杜修力,张楚汉.瑞利阻尼介质有限元离散模型动力分析的数值稳定性[J].地震工程与工程振动,2002,22(6):18-24.
    [8]景立平.波动有限元方程显式逐步积分格式稳定性分析[J].地震工程与工程振动,2004,24(5):20-26.
    [9]刘祥庆,刘晶波,丁桦.时域逐步积分算法稳定性与精度的对比分析[J].岩石力学与工程学报,2007,26(Supp.1):3000-3008.
    [10]邱启荣.矩阵理论及其应用[M].北京:中国电力出版社,2008.