A simplified procedure is presented for the determintion of rates of change of frequency response function matrix of system for the non-proportional changeable damping ,by means of the second order perturbation technique in conjunction with complex modal analysis. Consequently, when excitation information in frequency domain is known,the sensitivity of the response power spectra and the response time history can be obtained.