Computer-aided diagnosis of occlusal disorders and design of dental restorations requires an automated evaluation of jaw occlusion and chewing ability. This requires simulation of the motion of the jaws and characterization of contacts between the surfaces of teeth. We propose approaches to evaluation of the load on teeth and of the grinding process. These characteristics are derived in interactive time, and are based on distance maps and topological structure of the contact zones. The proposed approaches are general and usable in applications where modeling of contact between objects with complex geometry is required.