首页 | 本学科首页   官方微博 | 高级检索  
     检索      

各向异性介质中波动方程有限元法模拟及其稳定性
引用本文:周辉,何樵登.各向异性介质中波动方程有限元法模拟及其稳定性[J].地球物理学报,1997,40(6):833-841.
作者姓名:周辉  何樵登
作者单位:1. 青岛海洋大学地质地球物理研究所, 青岛 266003; 2. 长春地质学院地球物理系, 长春 130026
摘    要:利用六面体单元和三线性插值函数,推导出在任意弹性各向异性介质中三维三分量波动方程所满足的有限元方程.给出用时域有限差分求解该有限元方程时,时间、空间采样率和弹性参数应满足的稳定性条件.提出避免每一时间步都求解一大型线性代数方程组的方法.在开始求波场前,先求质量矩阵的逆的因子形式,以后每时间步只需作刚度矩阵的逆的因子形式与向量的乘积运算,可大大提高有限元素法在微机上进行波动方程数值模拟的运算速度.结合两个数值模拟算例,说明该方法的可行性,并与伪谱法作了计算速度的比较,有限元法的计算速度高于伪谱法的计算速度.

关 键 词:各向异性介质  波动方程  有限元法  稳定性条件  数值模拟  
收稿时间:1996-06-10

MODELING OF WAVE EQUATIONS IN ANISOTROPIC MEDIA USING FINITE ELEMENT METHOD AND ITS STABILITY C0NDITION
Zhou Hui, Xu Shi-zhe, Liu Bin.MODELING OF WAVE EQUATIONS IN ANISOTROPIC MEDIA USING FINITE ELEMENT METHOD AND ITS STABILITY C0NDITION[J].Chinese Journal of Geophysics,1997,40(6):833-841.
Authors:Zhou Hui  Xu Shi-zhe  Liu Bin
Institution:1. Institute of Geology &; Geophysics., Ocean University of Qingdao, Qingdao 266003, China; 2. Changchun University of Earth Sciences, Changchun 130026, China
Abstract:The finite element equati0n for arbitrary anisotropic elastic medium is derived when the domain is divided with hexahedron elements and linear interpolation funchons are used in three direchons. The stability condih0ns of the numerical method, which relthes the sampling rates of hme and space and elashc constalltS, are given in thes paper. A method for solving the finite element equahons is also discussed. The method avoids solving a linear system at each hme step, and only needs the product of the faCtor form of the inverse of the mass matrix and a displacement vector. So the computational rate of the method is very high. Two examples of mode1ing show that the method proposed in this paper is feasible.
Keywords:Anisotropic medium  Finite element method  Wave equation  Stability condition  Numeric modeling
本文献已被 CNKI 维普 等数据库收录!
点击此处可从《地球物理学报》浏览原始摘要信息
点击此处可从《地球物理学报》下载免费的PDF全文
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号