Abstract: In order to describe the elasto-plastic dynamic response of fluid-saturated porous media, incremental elasto-plastic wave equations of fluid-saturated porous media are developed by the fundmental theory of continuum mechanics and in accordunce with the characteristic of fluid-saturated porous media. The above equations are divided in the space domain to get their Galerkin formula, and these formulas are divided in the time domain with the integral method, which consists of the central difference method and the Newmark constant average acceleration method to get the explicit time integral formulas for solving the above wave equations. On the basis of the integral formulas mentioned above, the time-domain explicit finite element method is developed for calculation and analysis of the elasto-plastic dynamic response of fluid-saturated porous media. In this method, the decoupling technique is adopted and it does not need to solve a set of linear equations in each time step, so the compuer memory space can be saved considerably and the calculation speed can be increased sharply by using it.