Abstract
Time finite element method was used for analysis of multibody dynamics. Hamilton's law played an important role in the present formulation. Constraint equations were introduced in the functional through Lagrange's multiplier. A simple example showed that the existing method did not yield the stable numerical solutions once the constraint equations existed. The way for satisfying the constraint equations were modified in this paper so that the stable numerical solutions were obtained. Numerical results presented herein showed the effectiveness of this modification.
