This study presents computational fluid dynamics (CFD) simulation of forced vibration of a bridge deck for extraction of amplitude-dependent aerodynamic damping in the vicinity of vortex lock-in region. The coordinate system fixed on the oscillating body is used to bypass the requirement of mesh updating. The three-dimensional flow around the vibrating deck is assumed to be homogeneous in spanwise direction and simulated using the efficient high-order spectral element/Fourier method. The aerodynamic damping is extracted from the simulated lift and then used to estimate the steady-state vortex-induced vibration (VIV) amplitude for given wind velocity and structural parameters. A good match of VIV amplitudes between simulations and wind tunnel tests validates the accuracy of simulation. The VIV mechanisms are discussed from the vortex shedding modes, distributions of dynamic pressure and energy transfer from flow to bridge deck. The leading-edge vortex shedding rather than Karman vortex shedding is identified to be responsible for the VIV. The hysteresis of VIV is observed which is related to the mode change of impinging leading-edge vortex and leading-edge vortex shedding. Finally, a novel forced vibration technique with a varying amplitude is proposed that permits extraction of aerodynamic damping by a single CFD run with improved efficiency.