Modeling fracture in multilayered teeth using the finite volume-based phase field method.

The present work, utilizing the finite volume-based phase field method (FV-based PFM), aims to investigate the initiation and propagation of cracks in the second molar of the left mandible under occlusal loading. By reconstructing cone beam computed tomography scans of the patient, the true morphology and internal mesostructure of the entire tooth are implemented into numerical simulations, including both 2D slice models and a realistic 3D model. Weibull functions are introduced to represent the tooth's heterogeneity, enabling the stochastic distribution characteristics of mechanical parameters. The results indicate that stronger heterogeneity leads to greater crack tortuosity, uneven damage distribution, and lower fracture stress. Additionally, different cusp angles (50° and 70°) and pre-existing fissure morphologies (i.e., U-shape, V-shape, IK-shape, I-shape, and IY-shape) also significantly affect the mechanical performance of the tooth. The study reveals that different cusp angles affect the location of crack initiation. Overall, this work demonstrates the utility of the FV-based PFM framework in capturing the complex fracture behavior of teeth, which can contribute to improved clinical treatment and prevention of tooth fractures. The insights gained from this study can inform the design of dental crown restorations and the optimization of cusp inclination and contact during clinical occlusal adjustments.