A new time integration scheme for shock propagation in saturated soil

This paper presents a time integration scheme capable to simulate blast loading causing high frequency wave propagation in porous media using coarse meshes. The scheme is based on the partition of unity method. The discontinuity is imposed on the velocity field, while the displacement field is kept continuous. The developed time integration scheme is unconditionally stable and has controllable numerical dissipation in the high frequency range. The time scheme has been implemented in combination with Biot’s theory of wave propagation in two-phase media. The paper presents the basic steps of the integration scheme and 1D-examples to illustrate that the proposed time scheme is, in addition to being accurate and stable, highly effective for coarse meshes and enabling the control of numerical damping in a consistent way. This makes the developed scheme suitable for large scale finite element analysis. A comparative 2D-example of an induced ground shock due to an explosion in a tunnel is presented.