Geometric interferences must be detected and reflected in the equations of motion. Global and local searches are carried out to find the exact location of a contact. NURBS and triangular patch methods are used to represent the surfaces. Contact penalty forces are generated for a contact.