The method of analyzing flexible mechanisms with clearances was studied considering flexibility of beams in the mechanism using finite elements. Both ends of a beam were modeled as free following Dubowskys impact pair model. Instead some force constraints were imposed at imposed at the connections between adjoining links. Coulomb model has been developed using dry frictions in place of tangential damping forces in the impact pair model and the contact compliance and damping coefficient approximated in a form of root function were used. As examples, impacts of a rigid ball in a cylinder, impact beam model and four-bar mechanisms made up of three flexible links with clearance connections were simulated numerically. The results from examples showed similar but a little bit smaller magnitude of impact forces compared with published studies.