Three-dimensional finite element analysis of occlusal stress distribution in the human skull with premolar extraction

To analyze the effect of orthodontic treatment with premolar extraction on the stress distribution of the occlusal force in the human skull.

A three-dimensional finite element (3D FE) model was constructed based on computed tomography scan data, and it served as the pretreatment model. For the extraction model simulating postorthodontic occlusion, the first premolar was removed in the pretreatment model, and the anterior and posterior segments were repositioned. Stress distribution was evaluated by 3D FE analysis in both models under the simulation of 1000 N for occlusal forces and 400 N for masseter muscle force.

The occlusal stresses were concentrated at the alveolar bone near the teeth, the infrazygomatic crest, the frontal process, the temporal process of the zygomatic bone, the infraorbital rim, the pyriform aperture region, and the pterygoid plate in both models. The von Mises stress at the pterygoid plate area was lower in the extraction model (3.53 MPa) than in the pretreatment model (5.57 MPa), while the stress at the frontal process of the maxilla was higher in the extraction model (2.32 MPa) than in the pretreatment model (2.16 MPa).

The results indicated that the occlusal forces were transferred through the maxillonasal, maxillozygomatic, and maxillopterygoid stress trajectories and that stress distribution moved more "forward" with the orthodontic treatment with premolar extraction.