A　precise　time-integration　method　for　the　time-domain　solution　of　the　dynamic　consolidation　problem　of　fluid-saturated　porous　media　is　proposed　on　the　basis　of　u-p　dynamic　formulation.　The　solid-phase　displacement　u　is　computed　by　the　precise　time-integration　algorithm,　and　the　pore　fluid　pressure　p　is　computed　by　the　backward　difference　algorithm.　The　proposed　method　is　verified　by　a　standard　numerical　example.　The　characteristic　of　the　proposed　method　is　also　investigated.　The　computational　accuracy　is　compared　when　several　numerical　integral　methods　are　employed　to　calculate　the　special　solution　term　of　u-p　dynamic　formulation.　The　computational　accuracy　of　Gauss　integral　method　is　also　compared　when　different　numbers　of　integral　points　are　employed.　It　is　shown　that　(1)　the　proposed　precise　time-integration　method　has　fine　computational　accuracy.　(2)　Gauss　integral　method,　Simpson　integral　method　and　Coates　integral　method　all　have　fine　computational　accuracy,　while　trapezoidal　integral　method　has　poor　computational　accuracy　when　used　to　compute　the　special　solution　term.　(3)　the　effect　of　the　number　of　Gauss　integral　points　on　computational　accuracy　is　insignificant.　©　2017,　Editorial　Office　of　Chinese　Journal　of　Computational　Mechanics.　All　right　reserved.