This is a presentation of robust and accurate explicit time-stepping strategy for finite element modeling of elastic discontinuous wave propagation in strongly heterogeneous, multi-material and graded one-dimensional media. One of the major issues in FEM modeling is the existence of spurious numerical stress oscillations close to theoretical wave fronts due to temporal-spatial dispersion behavior of FE discretization. The numerical strategy presented for modeling of 1D discontinuous elastic waves is based on (a) pushforward-pullback local stepping — ensuring the elimination of dispersion due to different critical time step sizes of finite elements, (b) domain decomposition via localized Lagrange multipliers — to satisfy coupling kinematics and dynamic equations , (c) asynchronous time scheme — ensuring the correct information transfer of quantities for the case of integer ratios of time step size for all domain pairs. Dispersion behaviors, existence of spurious stress oscillations, and sensitivity of the dispersion to time step size are then suppressed. The proposed method is numerically tested with regard to the rectangular step pulse elastic propagation problem considering in-space varying Young’s modulus. To prove robustness and accuracy, a comparison with results from commercial software, an analytical solution, and experimental data from partial assembly of a split Hopkinson pressure bar (SHPB) setup is provided.