and Au F., “ Assessment and Improvement of Precise Time Step Integration Method,” Computers & Structures, Vol. 84, No. 12, 2006, pp. 779–786. and Dokainish M., “ A Survey of Direct Time-Integration Methods in Computational Structural Dynamics-2. and Subbaraj K., “ A Survey of Direct Time-Integration Methods in Computational Structural Dynamics-1. Bathe K.-J., Finite Element Procedures, Prentice–Hall, Upper Saddle River, NJ, 1996, pp. 769–785. It has been demonstrated that the improved PTI methods are accurate and insensitive to integration step length and, more importantly, are highly efficient for large-scale structures subjected to high-frequency excitations, which greatly extend the applicability of PTI methods. The effects of tunable parameters of both kinds of algorithms on the numerical efficiency are also discussed. The transient responses of a cantilever beam and an aeroengine bladed disk are used as numerical examples to validate the accuracy and efficiency of the proposed methods. The theories and implementations of the improved PTI methods are introduced. In both kinds of methods, only the multiplications of sparse matrix and vectors are involved. To remove the calculations of inverse and/or exponential of large-scale matrix in the precise time integration (PTI) method for transient response analysis of structures, new PTI methods equipped with the restarted Krylov subspace method and real Leja points interpolation method to calculate exp ( A ) b are presented in this paper.