稀疏存储的显式有限元三角网格地震波数值模拟及其犘犕.pdfVIP

第 卷 第 期 地 球 物 理 学 报 ， ５６ ９ Ｖｏｌ．５６ Ｎｏ．９ 　 　 　 　 　 年 月 ， ２０１３ ９ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ｓｅ．２０１３ 　 　 　 ｐ 刘有山，滕吉文，刘少林等 稀疏存储的显式有限元三角网格地震波数值模拟及其 吸收边界条件 地球物理学报， ， ． ＰＭＬ ． ２０１３ （）： ， ： ／ ５６９ ３０８５３０９９ｄｏｉ１０．６０３８ ｊｇ ， ， ， ＬｉｕＹＳ Ｔｅｎ ＪＷ ＬｉｕＳＬ ｅｔａｌ．Ｅｘｌｉｃｉｔｆｉｎｉｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄｗｉｔｈｔｒｉａｎｌｅｍｅｓｈｅｓｓｔｏｒｅｄｂ ｓａｒｓｅｆｏｒｍａｔａｎｄｉｔｓ ｇ ｐ ｇ ｙ ｐ （ ）， ， （）： ， ： ｅｒｆｅｃｔｌｍａｔｃｈｅｄｌａｅｒｓａｂｓｏｒｂｉｎ ｂｏｕｎｄａｒｃｏｎｄｉｔｉｏｎ．犆犺犻狀犲狊犲犑．犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ２０１３ ５６９ ３０８５３０９９ｄｏｉ１０． ｐ ｙ ｙ ｇ ｙ 狆狔 ／ ６０３８ ｊｇ 稀疏存储的显式有限元三角网格地震波数值模拟 及其犘犕犔吸收边界条件 ， １２ １ １ １ 刘有山 ，滕吉文 ，刘少林 ，徐 涛 　 １中国科学院地质与地球物理研究所，北京 １０００２９ 　 ２中国科学院大学，北京 １０００４９ 　 摘 要 有限元法是复杂介质地震模拟的有力工具，它能比较客观地反映地震波的传播，比较细致地再现地震图 　 　 像 但是，为了获得较精确的结果，有限元法模拟地震波的传播需要的网格点数多，具有计算量大和消耗内存多的 ． 缺点 针对上述缺点，本文对刚度矩阵采用压缩存储行（ ）格式，以减少计算量并节省内存；采用集中质量矩阵 ． ＣＳＲ 得到对角的质量矩阵以提高有限元法（显式有限元）的计算效率；时间离散采用保能量的Ｎｅｗｍａｒｋ算法以提高有 限元法的计算精度；采用变分形式（弱形式）的 吸收边界条件对人工截断边界进行处理 通过与高精度的数值 ＰＭＬ ． 方法———谱元法的数值试验的对比表明，上述方法的引入可使有限元法在计算精度