A new variational integrator is proposed to solve constrained mechanical systems. The main distinguishing feature of the present integrator comes from the distinct discretization of Lagrangians based on the Hamilton's principle in its most general form. Specifically, Hermite interpolation is used for discrete positions, which provides at least continuity for generalized coordinates. The velocities and momentums are interpolated using quadratic polynomials for the consistency, such that the kinematic relation between velocities and positions can be exactly satisfied. Meanwhile, the Gauss-Legendre quadrature rule is employed to guarantee the accuracy of discrete Lagrange equations. To tackle constrained mechanical systems, a coordinate partition approach is used to eliminate the constraint equations. The local incremental rotation vector is exploited to get rid of rotation singularities in spatial problems. Moreover, an adaptive stepsize strategy is implemented to improve the efficiency. The strengths of the new integrator lie in the accessible large step sizes in the simulation and its global second-order accuracy for positions as well as velocities. Several examples are performed and analyzed to validate its accuracy and capabilities.