In　order　to　describe　the　elastoplastic　dynamic　response　of　fluid-saturated　porous　media,　incremental　elastoplastic　wave　equations　of　fluid-saturated　porous　media　are　developed　by　the　fundamental　theory　of　continuum　mechanics　and　in　accordance　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　elastoplastic　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　computer　memory　space　can　be　saved　considerably　and　the　calculation　speed　can　be　increased　sharply　by　using　it.