A numerical algorithm for the dynamic contact analysis in flexible multibody systems is developed. The governing differential equations with algebraic constraints are derived using finite element approximations. The checking of contact status is made depending on the appropriate conditions such as normal distance, stick or slip, and so on. The contacting point in the target is obtained by a linear interpolation scheme. The resulting algorithm is applied to four example problems of different nature and shown applicable to realistic problems.