The GHM(Golla-Hughes -McTavish) model of ZN-1 Viscoelastic Material is combined with finite element method (FEM), which is the most frequently used method in Engineering. It can approach time-consuming iteration in solving model parameters and responses by introducing dissipation coordinates. The parameters of GHM model are determined by nonlinear curve fitting in complex frequency domain. This is converted into nonlinear optimization problem with constrained condition. The results show that the method proposed in the present paper to determine the parameters of GHM model is correct, simple and effective to employ GHM model to perform dynamic analysis.