三维非均匀介质中弹性波传播的数值模拟.pdfVIP

第 卷 第 期 固体 力 学 学报 VOI .22 NO .4 22 4 2001年 12月 ACTA MECHANICA SOLIDA SINICA December 2001 三维非均匀介质中弹性波传播的数值模拟! 张剑锋 刘铁林 （大连理工大学工业装备结构分析国家重点实验室，大连， ） （辽宁工学院土木工程系，锦州， ） 116023 121001 摘 要 提出了一种三维非均匀介质中弹性波传播数值模拟的方法，文中称为三维格子 法. 该算法是二维格子法（一种二维非均匀介质中P-SV 波传播的数值模拟算法）向三维非均匀 介质情况的推广. 在空间离散上该文方法与有限元方法类似，容许根据连续体的形状和介质分 界面任意剖分网格，且自然满足自由表面边界条件 不同于常规有限差分法在各个节点上满足 . 动力学微分方程，该算法通过满足各节点周围格子的整体平衡（积分平衡方程）来对问题进行求 解 三维格子法所需的计算机内存及计算耗时与同阶精度的规则网格有限差分法相当 算例表 . . 明，该文提出的三维格子法具有较高的精度且可很好地模拟三维复杂形状地表对弹性波的反射 和绕射. 关键词 三维弹性波模拟，三维格子法，有限差分法 ! 引言 有限差分方法是一种非均匀介质中弹性波数值模拟的有效方法，但由于其使用规则网 格，使得在处理复杂形状地表和不规则介质界面时出现困难 理论上，有限元法或有限体积 . ［］ 1 法 可处理这类复杂形状地表和不规则介质界面情况，但有限元法巨大的计算机内存需求 和计算耗时限制了它的应用，特别是应用于三维非均匀介质情况；而即使是二维情况，有限 体积法的自由表面边界条件也没能很好的解决. 某些学者针对三维介质模型中复杂形状地 ［ ］ 2 ~ 4 表的处理提出了一些改进的有限差分方法 ，以便更好地模拟复杂形状地表对地震波的 反射和绕射 然而，上述方法引入的自由表面边界条件或是非常繁琐或是带来某些近似，且 . ［］ 3 文献［］提出的弯曲网格方法并不适应于地表弯曲程度很大的情况 另外，用阶梯状网格 2 . 来近似处理不规则介质界面将导致人工绕射现象，特别是在使用高阶有限差分算法的情况 ［］ 时这种现象更加严重5 . 文献［］和［］综合了有限元法 、有限体积法和交错网格有限差分法的基本思想，给出 6 7 了一种二维非均匀介质中P-SV 波传播的数值模拟算法———格子法；该方法可灵活处理任意 形状地表和不规则介质界面，且自然满足自由表面边界条件. 本文进一步将这一方法推广到 三维非均匀介质情况，给出了一个可精细模拟复杂形状地表和不规则介质界面的三维弹性 波模拟算法———三维格子法 文中采用了一种四面体和平行六面体混合格子型的算法以便