Noise-vibration-harshness (NVH) design optimization problems have become major concerns in engineering design phase. For complex vehicle body structures, the numerous design variables and constraints make the direct optimization design inefficient needing enormous iterations. In this work, to reduce computational burden, the screening analysis is firstly performed to detect the most important variables prior to the design optimization. Then the response surface approximate method coupled with design of experiment (DOE) technique is employed to replace time-consuming finite element vibration analysis. Subsequently, an optimization of NVH performance considering dynamic response is formulated and solved using the adaptive simulated annealing (ASA) algorithm. The methodology aims at improving the NVH behavior of body structure by simultaneously suppressing several resonance peaks in frequency domain. Finally, the proposed method and process are successfully illustrated through a real-world car body-in-white example. The results demonstrate the proposed method of incorporating response surface model with ASA algorithm is feasible and cost-efficient in solving vibration optimization problem.