Herein are provided methods, systems, computer-readable media, techniques and processes for occlusion estimation in dental prosthesis design. These include determining a first contact point between antagonist teeth, then simulating motion in order to determine additional contact points. Various of the determined contact points may not be used in the simulated motion, based on the simulated physics. When one or more predetermined stopping criteria have been met, then the contact points may be used to define the relative occlusal positions of antagonists. In some embodiments, sets of teeth in a bridge or crown, e.g., may be simulated as moving separately and a new state for each set of teeth may be determined.