In this study, a numerical approach based on mid-point integrated finite elements and a viscous boundary is proposed for time-domain wave-propagation analyses in infinite poroelastic media. The proposed approach is accurate, efficient, and easy to implement in time-domain analyses. In the approach, an infinite domain is truncated at some distance. The truncated domain is represented by mid-point integrated finite elements with real element-lengths and a viscous boundary is attached to the end of the domain. Given that the dynamic behaviors of the proposed model can be expressed in terms of mass, damping, and stiffness matrices only, it can be implemented easily in the displacement-based finite-element formulation. No convolutional operations are required for time-domain calculations because the coefficient matrices are constant. The proposed numerical approach is applied to typical wave-propagation and soil-structure interaction problems. The model is verified to produce accurate and stable results. It is demonstrated that the numerical approach can be applied successfully to nonlinear soil-structure interaction problems.