Transient dynamic responses of a large-size finite element model can be nowadays solved in a reasonable computation time due to rapid improvement of both numerical schemes and computing resources.However, the increasing demands for accurate simulation and complicated modeling have led to larger and more complex finite element models, which consequently yield in much higher computation time.In addition, when structural optimizations include transient performance functions such as displacement, velocity, and acceleration, the optimizations often do not end within a reasonable process time because the large-size simulation should be repeated many times.