The determination of the interference geometry between two arbitrary objects is an essential problem encountered in the simulation of continuous contact dynamics and haptic interactions. In these applications, with known material properties, the interbody contact force is only a function of the interference geometry between two objects. Here a theoretical basis and algorithms for the calculation of the interference geometry, such as overlap region, contact area and normal, and interference volume, are presented. Two methods to obtain the contact area and normal are analyzed: an area-weighted method and a best-fitting method. The geometric properties of the area-weighted method are presented and the degenerate cases related to both methods are discussed. Methods to calculate the application point of an interbody contact force are discussed. Some numerical simulation results are presented based on the implementation of the geometric algorithms, which are verified by comparison with hand calculation. The continuity of contact normal and its application point are demonstrated for a case in which the contacting objects smoothly move with respect to each other in the simulation.