(15) 比较零频动态解和静态解结果(文献结果复现)
下载示例: ex15.tar.gz
由于在求解动态解过程中使用了复数频率,因此零频(\(\omega=0\))的动态解可以 近似 静态位移。 以下进行相关测试。
备注
虚频率的大小直接影响近似效果,鼓励感兴趣的读者自行进行更多测试。
半空间模型
Zhu and Rivera (2002) 对该问题以半空间模型为例进行了测试。 以下对文中 Figure 2 结果进行复现。具体参数详见原文和示例脚本。
以上仅比较了爆炸源的结果,进一步地,我们可以将六种震源的结果都进行比较。 从下图可看出,这种近似方法对部分震源而言误差还是比较明显的, 或者说,要求更小的虚频率。但虚频率过小会引入数值不稳定的问题。
层状模型
仅更换模型,其它参数设置与以上一致。