弱形式时域完美匹配层——滞弹性近场波动数值模拟

摘要

本文旨在构建适用于滞弹性近场时域波动有限元模拟的高精度人工边界条件:完美匹配层(Perfectly Matched Layer:PML),其中阻尼介质时域本构基于广义标准线性体建立.与以往研究不同,本文采用复坐标延拓技术变换弱形式波动方程构建了可直接用有限元离散的弱形式时域PML,规避以往独立对无限域内波动方程及界面条件进行延拓可导致的PML场方程和界面条件匹配不合理引发数值失稳、计算精度低下等问题.其次,针对PML中多极点有理分式与频域函数乘积的傅里叶反变换难以计算的问题,利用PML精度对复坐标延拓函数中延拓参数微调不敏感这一特点,明确给出了参数微调准则以规避多重极点,进而利用有理分式分解给出了一种普适、简便的计算方法,极大地简化了PML计算.基于该方法可实现任意高阶PML.最后,将本文构建滞弹性PML与高阶勒让德谱元(高精度集中质量有限元)结合得到滞弹性近场波动谱元离散方案.基于算例验证了滞弹性PML的计算效率、精度及新离散方案的长持时稳定特性.新离散方案可应用于计入实际介质阻尼的地震波动正、反问题数值模拟,提高波形模拟的精度以及地下波速结构反演的精度和可靠性.