基于混合有限元格式的完美匹配层与多次透射公式人工边界-地震学报.pdfVIP

第 卷 第 期 地 震 学 报 ， ２９ ６ Ｖｏｌ．２９ Ｎｏ．６ 　 　 　 　 　 　 年 月 （ ） ＡＣＴＡＳＥＩＳＭＯＬＯＧＩＣＡＳＩＮＩＣＡ ， ２００７ １１ ６４３ ６５３ Ｎｏｖ． ２００７ 　 　 ～ 　 文章编号： （ ） ０２５３３７８２２００７０６０６４３１１ 基于混合有限元格式的完美匹配层与 多次透射公式人工边界比较研究 ）， ，） ） 李 宁１ 谢礼立１ ２ 翟长海１ 　 　 　 １）中国哈尔滨１５００９０哈尔滨工业大学土木工程学院 ２）中国哈尔滨１５００８０中国地震局工程力学研究所 　 摘要 介绍了完美匹配层（ ）人工边界可以吸收不同频率和任意角度入射波的原理以及 ＰＭＬ 　 人工边界的构造方法 在此基础上，将 人工边界应用于地震波动数值模拟的速度 ＰＭＬ ． ＰＭＬ  应力混合有限元格式中，探讨了 ＰＭＬ应用的可行性，并通过数值试验研究了 ＰＭＬ人工边界 的反射率，比较了 人工边界与多次透射公式（ ）人工边界应用于体波和面波模拟中 ＰＭＬ ＭＴＦ 数值反射的差异，对两种边界的透射效果进行了分析 结果表明，尽管数值离散后 人工 ． ＰＭＬ 边界不再保持完美匹配特性，但ＰＭＬ人工边界在近场波动数值模拟中可获得比 ＭＴＦ人工边 界更为理想的吸收效果，在角点透射、大角度掠射情形下尤为明显；ＰＭＬ人工边界在混合有 限元格式的数值算法中，未见失稳等不良反应，比ＭＴＦ人工边界有更好的稳定性；在合理选 择参数的情况下，ＰＭＬ人工边界的运算量可接受． 关键词 完美匹配层 多次透射公式 弹性波 人工边界 混合有限元 　　 　 　 　 　 中图分类号：ＴＵ３５２．１ 文献标识码：Ａ 　　　 引言 人工边界的构造是近场地震波数值模拟过程中最重要的问题之一 在一定计算区域 ． 内，针对如何让外行波的能量在有限域边界处衰减、吸收这一问题，国内外学者提出了基 于各类型波动方程、适用于不同情况的全局或局部人工边界 （ ）在电磁波数