Ten years ago, an original semi-recursive formulation for the dynamic simulation of large-scale multibody systems was presented by García de Jalón et al. (Advances in Computational Multibody Systems, pp. 1-23, 2005). By taking advantage of the cut-joint and rod-removal techniques through a double-step velocity transformation, this formulation proved to be remarkably efficient. The rod-removal technique was employed, primarily, to reduce the number of differential and constraint equations. As a result, inertia and external forces were applied to neighboring bodies. Those inertia forces depended on unknown accelerations, a fact that contributed to the complexity of the system inertia matrix. In search of performance improvement, this paper presents an approximation of rod-related inertia forces by using accelerations from previous time-steps. Additionally, a mass matrix partition is carried out to preserve the accuracy of the original formulation. Three extrapolation methods, namely, point, linear Lagrange and quadratic Lagrange extrapolation methods, are introduced to evaluate the unknown rod-related inertia forces. In order to assess the computational efficiency and solution accuracy of the presented approach, a general-purpose MATLAB/C/C++ simulation code is implemented. A 15-DOF, 12-rod sedan vehicle model with MacPherson strut and multi-link suspension systems is modeled, simulated and analyzed.