Viscoelastic damping materials are widely used to reduce noise and vibration because of its low cost and easy implementation, for examples, on the body structure of passenger cars, air planes, electric appliances and ships. To design the damped structures, the material property such as elastic modulus and loss factor is essential information. The four-parameter fractional derivative model well describes the dynamic characteristics of the viscoelastic damping materials with respect to both frequency and temperature. However, the identification procedure of the four-parameter is very time-consuming one. In this study a new identification procedure of the four-parameters is proposed by using an FE model and a gradient-based numerical search algorithm. The identification procedure goes two sequential steps to make measured frequency response functions(FRF) coincident with simulated FRFs: the first one is a peak alignment step and the second one is an amplitude adjustment step. A numerical example shows that the proposed method is useful in identifying the viscoelastic material parameters of fractional derivative model.