For efficient calculation of natural modes of structures, a numerical scheme which accelerates convergence of the subspace iteration method by employing accelerated starting Lanczos vectors was proposed in 2005. This paper is an extension of the study. The previous study simply showed feasibility of the proposed method by analyzing structures with smaller degrees of freedom. While, the present study verifies efficiency of the proposed method more rigorously by comparing closeness of conventional and accelerated starting vectors to genuine eigenvectors. This study also analyzes an example structure with larger degrees of freedom and more complex constraints in order to investigate applicability of the proposed method.