Objective To reconstruct three-dimensional finite element method models (FEM), including FEM for lower dentition and bilateral condyles for analyzing the characteristics of temporomandibular joint (TMJ) stress distribution result from alteration of dental occlusion . Methods A human cadaveric skull was used for CT scan at the plane parallel to the Frankfort Horizontal plane. Geometrical information of mandibular CT images was extracted for preparation of being imported into software. Occlusion surface of its plaster casts was then scanned in a 3D fashion (accuracy:0.5mm). The related parameters of both cast and CT images were converted into a coordinate system. A simulated disc was added on each of the condyles. Constraints were applied to the upper surface of the disc and the mandibular angle. Results After meshing the FEM containing bilateral condyles?simulated discs, and lower dentition was established. Conclusion By using the methods developed for this study, the FEM was established to study the result of TMJ stress distribution and its alterations accompanied with changes in occlusion.