The step-by-step integration is the most frequently adopted way to obtain the dynamic responses of complex nonlinear structure system. The major topic of this study is to develop a new integration method which can have computational efficiency, unconditional stability, and favorable numerical dissipation, which can accurately integrate the low frequency modes while it can effectively filter out the spurious participation of high frequency modes. Some numerical examples and actual pseudodynamic tests were conducted to confirm the numerical properties of the proposed integration method, especially the characteristics of favorable numerical dissipation.