廖振鹏,
谢志南
中国地震局工程力学研究所, 中国地震局地震工程与工程振动重点实验室, 哈尔滨 150080
基金项目: 中央级公益性科研院所基本科研业务费专项（2017B08），国家自然科学基金（51808516，U2039209），黑龙江省"头雁"创新团队计划资助
详细信息
作者简介: 章旭斌, 男, 1986年生, 博士, 助理研究员, 主要研究方向为地震波动数值模拟.E-mail: zhangxvbin.hebi@163.com
中图分类号: P315 收稿日期:2020-11-03
修回日期:2021-06-28
上线日期:2021-10-10
Mechanism of high frequency instability and stable implementation for transmitting boundary—P-SV wave motion
ZHANG XuBin,LIAO ZhenPeng,
XIE ZhiNan
Institute of Engineering Mechanics, China Earthquake Administration, Key Laboratory of Earthquake Engineering and Engineering Vibration of China Earthquake Administration, Harbin 150080, China
MSC: P315
--> Received Date: 03 November 2020
Revised Date: 28 June 2021
Available Online: 10 October 2021
摘要
摘要:波动数值模拟的稳定性是获得可靠结果的前提.透射边界是一类具有高阶、高效等特点的人工边界，其引发的高频失稳是由内域格式和透射边界的不当耦合所致.本文针对P-SV波动有限元模拟中透射边界引发的失稳问题，基于GKS定理的群速度解释，通过对有限元和透射边界的频散分析揭示了数值失稳机理为透射边界和相邻内域格式支持了群速度指向内域的高频P波或SV波，波动能量将从边界进入内域引发数值失稳.同时，对比连续模型频散指出引发失稳的谐波是由有限元离散引入.本文采用修改的数值积分方法调整有限元刚度，以消除有限元中引发边界失稳的高频波动成分，从而稳定实现透射边界.理论分析和数值实验均表明本文稳定措施的有效性.
关键词: 波动模拟/
高频失稳/
GKS定理/
频散分析/
数值积分
Abstract:The stability of wave numerical simulation is the premise of obtaining reliable results. Transmitting boundary is a kind of artificial boundary with high order and high efficiency, but there exists high frequency instability induced by the improper coupling between the interior scheme and artificial boundary scheme. Aiming at the instability problem caused by transmitting boundary in the finite element simulation of P-SV wave motion, based on the group velocity interpretation of GKS theorem, the mechanism of instability is further clarified by the dispersion analysis of finite element and transmitting boundary scheme. The mechanism is that the interior scheme and artificial boundary scheme support high frequency P wave or SV wave whose group velocity points to interior domain, thus wave energy is allowed to enter interior domain spontaneously from artificial boundary to cause instability. At the same time, compared with the dispersion of the continuous model, it shows that the high frequency wave components leading to instability are introduced by the finite element discretization. In order to eliminate the above wave components, modified integration rules are used to change the stiffness of the finite element. Theory analysis and numerical experiments show that the proposed stabilization measures can ensure transmitting boundary to be implemented stably.
Key words:Wave simulation/
High frequency instability/
GKS theorem/
Dispersion analysis/
Numerical integration
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=452b1cab-d3de-407f-91a1-828a2596c478
