Fracture parameters evaluation for the cracked nonhomogeneous enamel based on the finite element method and virtual crack closure technique

To accurately solve the fracture parameters of enamel, we have established computational nonhomogeneous enamel models and constructed the fracture element of enamel dumb nodes, based on the enamel mineral concentration, nonhomogeneous mechanical properties, and virtual crack closure technique. Through the commercial finite element software ABAQUS and the fracture element of the enamel dumb nodes, we have established the user subroutines UMAT and UEL, which enabled solving of the energy release rates of the nonhomogeneous enamel structure with cracks. The stress intensity factors of central cracks, three-point bend and compact stretched enamels, and double-edge notched stretched enamels are determined. By comparing them with analytical solutions, we have proved that the fracture element of the enamel dumb nodes is highly accurate, simple, and convenient. In addition, the cracks can be other elements rather than singular or special elements; they show versatility and other advantages. The stress intensity factor of the dental enamel can be solved more realistically. Thus, a new numerical method for prevention and treatment of dental diseases is provided.