论文部分内容阅读
本文首先将计算地震波区域划分成许多四边形单元，然后基于二维弹性波动方程，从其弱形式出发，通过对方程进行空间离散，在每个单元上把近似解表示成截断的Lagrange正交多项式展开，时间域中引入时间辛算法；最后用Galerkin方法求解正交问题的差分格式，得到全局保结构解。由于在空间和时间域均采用了保辛算法，因此较非辛算法具有较好的长时程模拟能力。