我想解决弹性动力学方程，以解决冲击快速缩放的问题（尽管没有冲击）。

$\rho\ddot{\mathbf{u}} - \nabla \cdot \sigma(\mathbf{u}) = \mathbf{f}$

在结构动力学中，该方程通常使用 Newmark 方案求解，该方案具有用于空间域离散化的拉格朗日元素。Newmark 方案是隐式的，但我相信考虑到快速的规模，使用显式方案解决会更便宜。为了避免从二阶导数项中反转质量矩阵，我已经看到一种常见的技术是将质量矩阵集中到对角矩阵中。这种技术的准确性如何？

其他方法涉及使用 DG 进行空间离散化，它给出了一个可以轻松反转的块对角矩阵，这是否比上述方法更便宜、更准确？解决这个问题的正确方法是什么？谢谢。