弱形式时域完美匹配层.PDF

第 卷 第 期 地 球 物 理 学 报 ， ６０ １０ Ｖｏｌ．６０ Ｎｏ．１０ 　 　 　 　 　 年 月 ， ２０１７ １０ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ｏｃｔ．２０１７ 　 　 　 谢志南，章旭斌 弱形式时域完美匹配层 地球物理学报， （ ）： ， ： ／ ．２０１７． ． ６０１０ ３８２３３８３１ｄｏｉ１０．６０３８ ｊｇ ， （ ）， （ ）： ＸｉｅＺＮ Ｚｈａｎ ＸＢ．２０１７．Ｗｅａｋｆｏｒｍｔｉｍｅｄｏｍａｉｎ ｅｒｆｅｃｔｌｍａｔｃｈｅｄｌａｅｒ．犆犺犻狀犲狊犲犑．犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ６０１０ ３８２３ ｇ ｐ ｙ ｙ 狆狔 ， ： ／ ３８３１ｄｏｉ１０．６０３８ ｊｇ 弱形式时域完美匹配层 谢志南，章旭斌 中国地震局工程力学研究所，中国地震局地震工程与工程振动重点实验室，哈尔滨　１５００８０ 摘要 应用高精度人工边界条件可有效提升近场波动数值模拟计算效率 完美匹配层是吸收层形式高精度人工边 ． 　 界条件，匹配层内场方程和界面条件通常分别采用复坐标延伸技术变换强形式无限域内波动方程和界面条件得 到，亦曾将无限域界面条件当作匹配层界面条件 场方程和界面条件构建过程相互独立，可能出现匹配不合理而引 ． 发数值失稳、计算精度低下等问题 本文提出采用复坐标延伸技术变换弱形式无限域波动方程以构建完美匹配层 ． 的方法 弱形式波动方程耦合了波动方程及界面条件，进而规避了变换后所得场方程与界面条件之间的匹配不合 ． 理问题 新方法可直接建立弱形式匹配层，在此基础上亦可给出强形式匹配层 弱形式便于有限元离散，强形式便 ． ． 于有限差分离散 基于弱形式完美匹配层，结合勒让德谱元建立了弹性介质近场波动谱元模拟方案 利用算例验证 ． ． 了新方案的精度及数值稳定性 本文工作可直接推广至多相耦合介质近场波动数值模拟 ． ． 关键词 完美匹配层；人工边界条件；有限元；谱单元；地震波动数值模拟 　 ： ／ 中图分类号 ， 收稿日期 ， 收修定稿 ｄｏｉ１０．６０３８ Ｐ３１５Ｐ６３１ ２０１７０１０９２０１７０５０５ ｊｇ