地震波动数值模拟及稳定性研究

唐 晖

（中国地震局工程力学研究所，哈尔滨150080）

地震波动数值模拟及稳定性研究

唐 晖

（中国地震局工程力学研究所，哈尔滨150080）

本文针对显式有限元－有限差分法进行了系统的研究，包括理论分析、算法改进以及实际应用。通过理论推导定量地研究了采用这种数值模拟方法对波动传播规律的影响，以及其在结合局部透射人工边界进行计算时的稳定条件。就采用该方法结合局部透射人工边界求解波源问题和散射问题的情况分别提出了改进计算的措施，并通过数值试验验证其有效性。提出了一种有限元和显式有限元－有限差分法耦合方法，并采用其对可控震源激震过程进行了数值模拟。利用在汶川地震中获得的数据，进行场地效应分析。主要工作如下：

（1）根据由李小军等人提出的一种求解动力方程的显式积分格式应用于波动有限元模拟时的递推形式，给出了波动传递函数的形式；以一维均匀离散体系为模型，推导了该情况下数值波动的传递函数，并进行数值试验，验证了传递函数求解的正确性；并通过分析求解所得的传递函数，探讨了这种显式积分格式对离散网格中波动传播规律包括截止频率、频散现象和能量耗散等方面的影响。

（2）针对显式有限元－有限差分法结合局部透射边界进行波动数值模拟的情况，推导并计算了循环系数，给出了计算过程的稳定条件，并通过数值试验对该判定准则进行验证；在给出的循环系数的基础上，定量地分析了积分格式的能耗特性对由边界条件引入的高频失稳的抑制和消除作用；并探讨了如何更有效地利用能耗特性保证计算稳定性，即实际计算模型中时间离散步距以及空间模型大小等参数如何取值既能确保计算稳定又能最大程度减小计算量。

（3）本文针对采用显式有限元－有限差分法结合局部透射人工边界求解波源问题和散射问题的情况分别提出了改进计算的措施。对于波源问题，建议针对每个边界节点可以采用不同的人工波速，该人工波速可以通过试算一定时间长度，然后通过相关处理获得。通过数值实验证明了，在波源问题的求解中采用这种方法选择人工波速可以在不增加计算失稳风险的同时提高数值计算精度。对于散射问题，本文给出了一种关于波动输入方法的改进措施，即在分区分离波场中入射波不再采用解析解而是采用离散数值解，并通过数值实验证明了采用该方法可以有效地抑制漂移失稳。

（4）本文给出了一种基于有限元和显式有限元－有限差分法相耦合的区域分裂法。将计算模型分为两个区域，对于广义结构及其附近区域采用有限元法模拟，而土体其他部分采用显式有限元－有限差分法，两个区域的模拟通过搭接区域实现数据的结合。通过分别代表波源问题和散射问题的数值试验都验证了该方法的有效性。采用提出的耦合方法针对可控震源激振过程建立的数值计算模型，并对不同的工况进行计算，研究分析了各因素包括接触土体密度、接触土体刚度、接触土层厚度、软、硬夹层以及底平板半径对可控震源能量输入大地体系的影响，为如何能高效地将激震能量输入土体提供了参考。

（5）利用自贡地形影响台阵在2008年汶川地震中获得的主震强震记录，采用水平／垂直谱比法（HVSR）对该山脊地形场地对地震动的放大效应进行了分析，并探讨了所获记录结果的可靠性；采用显式有限元－有限差分法结合局部透射人工边界的二维场地影响分析模型与方法，模拟了该场地的地震动响应。对模拟结果与HVSR法分析结果进行对比分析，表明数值模拟结果与观测记录分析结果具有较好的一致性，利用数值计算方法建立的简单二维分析模型基本上可以反映山脊地形对地震动的影响。

显式有限元－有限差分法；局部透射人工边界；传递函数；耦合法；场地效应

P315；

A；

10.3969／j.issn.0235－4975.2012.05.011

（作者电子信箱，唐晖：white＿hit＠163.com）

2012－05－02。