The suppression of aerodynamic response of long-span suspension bridges during erection and after completion by using single TMD and multi TMD is presented in this paper. An advanced finite-element-based aerodynamic model that can be used to analyze both flutter instability and buffeting response in the time domain is also proposed. The frequency-dependent flutter derivatives are transferred into a time-dependent rational function, through which the coupling effects of three-dimensional aerodynamic motions under gusty winds can be accurately considered. The modal damping of a structure-TMD system is analyzed by the state-space approach. The numerical examples are performed on the Akashi Kaikyo Bridge with a main span of 1990 m. The bridge is idealized by a three-dimensional finite-element model consisting of 681 nodes. The results show that when the wind velocity is low, about 20 m/s, the multi TMD type 1 (the vertical and horizontal TMD with 1% mass ratio in each direction together with the torsional TMD with ratio of 1% mass moment of inertia) can significantly reduce the buffeting response in vertical, horizontal and torsional directions by 8.6-13%. When the wind velocity increases to 40 m/s, the control efficiency of a multi TMD in reducing the torsional buffeting response increases greatly to 28%. However, its control efficiency in the vertical and horizontal directions reduces. The results also indicate that the critical wind velocity for flutter instability during erection is significantly lower than that of the completed bridge. By pylon-to-midspan configuration, the minimum critical wind velocity of 57.70 m/s occurs at stage of 85% deck completion.