我正在尝试对地震波方程进行建模，因此一直在阅读有关离散化方案及其稳定性的信息。我最近看到一篇关于“波动方程的 Galerkin FEM 方法”的有见地的论文http://www.math.tamu.edu/~bangerth/publications/2010-waves.pdf。在本文档和许多其他文档中，对于 newmark beta 方案，波动方程的最终离散形式是加速度、速度和位移。这也是我在其他文件中看到的方式。

但是，我还看到了另一种形式的 Newmark 波动方程方案，即以下形式：（Pg 130）

http://www.lmn.pub.ro/~daniel/ElectromagneticModelingDoctoral/Books/Computational%20EM/Bondeson-%20Computational%20%20Electromagnetics.pdf

现在，我对波动方程如何具有两种不同形式的 Newmark 离散化感到非常困惑。如果您能对此有所了解或指出正确的方向，我将不胜感激。

最终目的是在 fenics 和 DEAL 2 中对其进行建模