Establishment of a finite element model of a neonate's skull to evaluate the stress pattern distribution resulting during nasoalveolar molding therapy of cleft lip and palate patients

Nasoalveolar Molding (NAM) is associated with ambivalent acceptance regarding effectiveness and unknown long-term results. Our purpose was to analyze the stress distribution patterns within the viscero- and neurocranium of neonates during the first phase of NAM therapy. A finite element (FE) model of a healthy four-week-old neonate was generated, derived from a computed tomography scan allowing the implementation of a bone-density-dependent material model. The influence of dental germs with variable material properties, the cleft width and area of expected force application were analyzed in a worst-case scenario. The resulting stress distribution patterns for each situation were analyzed using the software Ansys APDL. The established FE model was verified with a convergence analysis. Overall, stress patterns at the age of four weeks showed von Mises stress values below 60.000 Pa in the viscero- and neurocranium. The influences of the allocation of material properties for the dental germs, the area of force application, and the cleft width were negligible. A workflow to simulate the stress distribution and deformation in neonates attributable to various areas of force application has been established. Further analyses of the skulls of younger and older neonates are needed to describe the stress distribution patterns during NAM therapy.