Improving the validation of finite element models with quantitative full-field strain comparisons

Data-driven cranial suture growth model enables predicting phenotypes of craniosynostosis

We present the first data-driven pediatric model that explains cranial sutural growth in the pediatric population. We segmented the cranial bones in the neurocranium from the cross-sectional CT images of 2068 normative subjects (age 0-10 years), and we used a 2D manifold-based cranial representation to establish local anatomical correspondences between subjects guided by the location of the cranial sutures. We designed a diffeomorphic spatiotemporal model of cranial bone development as a function of local sutural growth rates, and we inferred its parameters statistically from our cross-sectional dataset. We used the constructed model to predict growth for 51 independent normative patients who had longitudinal images. Moreover, we used our model to simulate the phenotypes of single suture craniosynostosis, which we compared to the observations from 212 patients. We also evaluated the accuracy predicting personalized cranial growth for 10 patients with craniosynostosis who had pre-surgical longitudinal images. Unlike existing statistical and simulation methods, our model was inferred from real image observations, explains cranial bone expansion and displacement as a consequence of sutural growth and it can simulate craniosynostosis. This pediatric cranial suture growth model constitutes a necessary tool to study abnormal development in the presence of cranial suture pathology.

Effect of crown-to-implant ratio and crown height space on marginal bone stress: a finite element analysis

BACKGROUND

Crown-to-implant ratio and crown height space, associated with the use of short implants, have been related with marginal bone loss. However, it is unclear which of the two entities would play the most important role on the bone remodelling process. Using a finite element analysis, the present work aims to help clarifying how those two factors contribute for the stress generation at the marginal bone level. A numerical model (reference model), with a crown-to-implant ratio of 4, was double validated and submitted to a numerical calculation. Then, it was modified in two different ways: (a) by decreasing the prosthetic height obtaining crown-to-implant ratios of 3, 2.5 and 2 and (b) by increasing the implants length obtaining a crown-to-implant ratio of 2.08. The new models were also submitted to numerical calculations.

RESULTS

The reference model showed a marginal bone stress of 96.9 MPa. The increase in the implants' length did not show statistically significant differences in the marginal bone stress (p-value = 0.2364). The decrease in the prosthetic height was accompanied with a statistically significant decrease in the marginal bone stresses (p-value = 2.2e- 16).

CONCLUSIONS

The results represent a paradigm change as the crown height space appears to be more responsible for marginal bone stress than the high crown-to-implant ratios or the implants' length. New prosthetic designs should be attempted to decrease the stress generated at the marginal bone level.

An open-source tool for the validation of finite element models using three-dimensional full-field measurements

Three-dimensional (3D) full-field measurements provide a comprehensive and accurate validation of finite element (FE) models. For the validation, the result of the model and measurements are compared based on two respective point-sets and this requires the point-sets to be registered in one coordinate system. Point-set registration is a non-convex optimization problem that has widely been solved by the ordinary iterative closest point algorithm. However, this approach necessitates a good initialization without which it easily returns a local optimum, i.e. an erroneous registration. The globally optimal iterative closest point (Go-ICP) algorithm has overcome this drawback and forms the basis for the presented open-source tool that can be used for the validation of FE models using 3D full-field measurements. The capability of the tool is demonstrated using an application example from the field of biomechanics. Methodological problems that arise in real-world data and the respective implemented solution approaches are discussed.

Comparison of ultimate force revealed by compression tests on extracted first premolars and FEA with a true scale 3D multi-component tooth model based on a CBCT dataset

OBJECTIVE

The aim of this study was to develop a new method for creating a multi-component and true scale 3-dimensional (3D) model of a human tooth based on cone-beam computed tomography (CBCT) images.

MATERIALS AND METHODS

First maxillary premolar tooth model was reconstructed from a patient's CBCT images. The 2D serial sections were used to create the 3D model. This model was used for finite element analysis (FEA). Model validation was performed by comparing the ultimate compressive force (UF) obtained experimentally using a universal testing machine and from simulation. The simulations of three component-omitting models (silicone, cementum, and omitting both) were performed to analyze the maximum (max.) principal stress and stress distribution.

RESULTS

The simulation-based UF indicating tooth fracture was 637 N, while the average UF in the in vitro loading was 651 N. The discrepancy between the simulation-based UF and the experimental UF was 2.2%. From the simulation, the silicone-omitting models showed a significant change in max. principal stress, resulting in a UF error of 26%, whereas there was no notable change in the cementum-omitting model.

CONCLUSION

This study, for the first time, developed a true scale multi-component 3D model from CBCT for predicting stress distribution in a human tooth.

CLINICAL RELEVANCE

This study proposed a method to create 3D modeling from CBCT in a true scale and multi-component manner. The PDL-like component-omitting simulation led to a higher error value of UF, indicating the importance of multi-component tooth modeling in FEA. Tooth 3D modeling could help determine mechanical failure in dental treatments in a more precise manner.

An assessment of the role of the falx cerebri and tentorium cerebelli in the cranium of the cat (Felis silvestris catus)

The falx cerebri and the tentorium cerebelli are two projections of the dura mater in the cranial cavity which ossify to varying degrees in some mammalian species. The idea that the ossification of these structures may be necessary to support the loads arising during feeding has been proposed and dismissed in the past, but never tested quantitatively. To address this, a biomechanical model of a domestic cat (Felis silvestris catus) skull was created and the material properties of the falx and tentorium were varied for a series of loading regimes incorporating the main masticatory and neck muscles during biting. Under these loading conditions, ossification of the falx cerebri does not have a significant impact on the stress in the cranial bones. In the case of the tentorium, however, a localized increase in stress was observed in the parietal and temporal bones, including the tympanic bulla, when a non-ossified tentorium was modelled. These effects were consistent across the different analyses, irrespective of loading regime. The results suggest that ossification of the tentorium cerebelli may play a minor role during feeding activities by decreasing the stress in the back of the skull.

Modelling human skull growth: a validated computational model

During the first year of life, the brain grows rapidly and the neurocranium increases to about 65% of its adult size. Our understanding of the relationship between the biomechanical forces, especially from the growing brain, the craniofacial soft tissue structures and the individual bone plates of the skull vault is still limited. This basic knowledge could help in the future planning of craniofacial surgical operations. The aim of this study was to develop a validated computational model of skull growth, based on the finite-element (FE) method, to help understand the biomechanics of skull growth. To do this, a two-step validation study was carried out. First, an in vitro physical three-dimensional printed model and an in silico FE model were created from the same micro-CT scan of an infant skull and loaded with forces from the growing brain from zero to two months of age. The results from the in vitro model validated the FE model before it was further developed to expand from 0 to 12 months of age. This second FE model was compared directly with in vivo clinical CT scans of infants without craniofacial conditions (n = 56). The various models were compared in terms of predicted skull width, length and circumference, while the overall shape was quantified using three-dimensional distance plots. Statistical analysis yielded no significant differences between the male skull models. All size measurements from the FE model versus the in vitro physical model were within 5%, with one exception showing a 7.6% difference. The FE model and in vivo data also correlated well, with the largest percentage difference in size being 8.3%. Overall, the FE model results matched well with both the in vitro and in vivo data. With further development and model refinement, this modelling method could be used to assist in preoperative planning of craniofacial surgery procedures and could help to reduce reoperation rates.

Form-function relationships in dragonfly mandibles under an evolutionary perspective

Functional requirements may constrain phenotypic diversification or foster it. For insect mouthparts, the quantification of the relationship between shape and function in an evolutionary framework remained largely unexplored. Here, the question of a functional influence on phenotypic diversification for dragonfly mandibles is assessed with a large-scale biomechanical analysis covering nearly all anisopteran families, using finite element analysis in combination with geometric morphometrics. A constraining effect of phylogeny could be found for shape, the mandibular mechanical advantage (MA), and certain mechanical joint parameters, while stresses and strains, the majority of joint parameters and size are influenced by shared ancestry. Furthermore, joint mechanics are correlated with neither strain nor mandibular MA and size effects have virtually play no role for shape or mechanical variation. The presence of mandibular strengthening ridges shows no phylogenetic signal except for one ridge peculiar to Libelluloidea, and ridge presence is also not correlated with each other. The results suggest that functional traits are more variable at this taxonomic level and that they are not influenced by shared ancestry. At the same time, the results contradict the widespread idea that mandibular morphology mainly reflects functional demands at least at this taxonomic level. The varying functional factors rather lead to the same mandibular performance as expressed by the MA, which suggests a many-to-one mapping of the investigated parameters onto the same narrow mandibular performance space.

The biomechanical function of periodontal ligament fibres in orthodontic tooth movement

Orthodontic tooth movement occurs as a result of resorption and formation of the alveolar bone due to an applied load, but the stimulus responsible for triggering orthodontic tooth movement remains the subject of debate. It has been suggested that the periodontal ligament (PDL) plays a key role. However, the mechanical function of the PDL in orthodontic tooth movement is not well understood as most mechanical models of the PDL to date have ignored the fibrous structure of the PDL. In this study we use finite element (FE) analysis to investigate the strains in the alveolar bone due to occlusal and orthodontic loads when PDL is modelled as a fibrous structure as compared to modelling PDL as a layer of solid material. The results show that the tension-only nature of the fibres essentially suspends the tooth in the tooth socket and their inclusion in FE models makes a significant difference to both the magnitude and distribution of strains produced in the surrounding bone. The results indicate that the PDL fibres have a very important role in load transfer between the teeth and alveolar bone and should be considered in FE studies investigating the biomechanics of orthodontic tooth movement.