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

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

论文摘要

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

论文目录

  • 0 引言
  • 1 无限域滞弹性波动方程
  • 2 弱形式时域滞弹性PML
  •   2.1 利用虚位移原理建立弱形式无限域波动方程
  •   2.2 将弱形式波动方程变换至频域,并对其做复坐标延拓得到弱形式频域PML
  •   2.3 对弱形式频域PML进行坐标变换将其转换为实坐标系下弱形式频域PML
  • 3 任意多极点有理分式与频域函数乘积傅里叶反变换的简化计算及滞弹性PML时空离散
  • 4 数值算例
  •   4.1 全空间均匀介质波动模拟
  •   4.2 半空间均匀介质波动模拟
  •   4.3 层状半空间波动模拟
  • 5 结论
  • 附录A二维弱形式频域PML的详尽表达式
  • 文章来源

    类型: 期刊论文

    作者: 谢志南,郑永路,章旭斌,唐丽华

    关键词: 滞弹性,广义标准线性体,完美匹配层,勒让德谱元,波动模拟

    来源: 地球物理学报 2019年08期

    年度: 2019

    分类: 基础科学,工程科技Ⅰ辑

    专业: 地质学,地球物理学,矿业工程

    单位: 中国地震局工程力学研究所中国地震局地震工程与工程振动重点实验室,新疆维吾尔自治区地震局

    基金: 国家重点研发计划(2018YFC1504004),国家自然科学基金(51678539),地震星火计划(XH19051),黑龙江省级资助(GX16C006)联合资助

    分类号: P315;P631.4

    页码: 3140-3154

    总页数: 15

    文件大小: 660K

    下载量: 98

    相关论文文献

    • [1].DP钢板非线性滞弹性力学行为及对回弹预测的应用研究[J]. 塑性工程学报 2018(02)
    • [2].Q&P钢板滞弹性变形行为研究[J]. 塑性工程学报 2016(06)
    • [3].饱和岩石滞弹性弛豫机理的实验研究[J]. 实验力学 2011(03)
    • [4].风的垂直切变对滞弹性大气中静力适应过程的影响[J]. 地球物理学报 2018(01)
    • [5].阻尼减振材料滞弹性位移场模型参数寻优及其应用[J]. 振动与冲击 2011(08)
    • [6].下地幔中R_(s/p)区域模型的形成机理、地震波速异常的起因和不连续界面的性质[J]. 西南大学学报(自然科学版) 2010(09)
    • [7].建立η坐标系下的假不压模式及其数值试验[J]. 大气科学 2018(01)
    • [8].锆基非晶合金准静态循环加载过程滞弹性行为的反速率依赖性[J]. 稀有金属材料与工程 2017(09)

    标签:;  ;  ;  ;  ;  

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

    猜你喜欢