New enriched FEM that accurately simulates wave propagation in fractured media. Combines the Phantom Node Method to model fractures and the Generalized FEM. Eliminates spurious oscillations that appear in FEM models of fractured media. Critical time step size is the same order of magnitude of that of the underlying FEM. A specific kinetic energy-consistent lumped mass matrix is proposed.