Based on the finite element method (FEM), the parametric variational principle (PVP) is combined with a numerical time-domain integral method to simulate the dynamic behavior of the pantograph-catenary system. Based on PVP, formulations for the nonlinear droppers in the catenary and for the contact between the pantograph and the contact wire are proposed. The formulations can accurately determine the tension state or compression state of the nonlinear droppers and the contact state between the pantograph and the contact wire. Based on the periodicity of the catenary and the precise integration method (PIM), a numerical timeintegration method is developed for the dynamic responses of the catenary. For this method, the matrix exponential of only one unit cell of the catenary is computed, which greatly improves the computational efficiency. Moreover, the validation shows that the formulations can compute the contact force accurately and represent the nonlinearity of the droppers, which demonstrates the accuracy and reliability of the proposed method. Finally, the dynamic behaviors of the pantograph-catenary system with different types of catenaries are simulated.