An Equivalent Constitutive Model of Cancellous Bone With Fracture Prediction

Inverse finite element analysis for an axisymmetric model of vertical tooth extraction

BACKGROUND AND OBJECTIVE

Tooth extraction is a common clinical procedure with biomechanical factors that can directly influence patient outcomes. Recent development in atraumatic extraction techniques have endeavoured to improve treatment outcomes, but the characterization of extraction biomechanics is sparse. An axisymmetric inverse finite element (FE) approach is presented to represent the biomechanics of vertical atraumatic tooth extraction in an ex-vivo swine model.

METHODS

Geometry and boundary conditions from the model are determined to match the extraction of swine incisors in a self-aligning ex vivo extraction experiment. Material parameters for the periodontal ligament (PDL) model are determined by solving an inverse FE problem using clusters of data obtained from 10 highly-controlled mechanical experiments. A seven-parameter visco-hyperelastic damage model, based on an Arruda-Boyce framework, is used for curve fitting. Three loading schemes were fit to obtain a common set of material parameters.

RESULTS

The inverse FE results demonstrate good predictions for overall force-time curve shape, peak force, and time to peak force. The fit model parameters are sufficiently consistent across all three cases that a coefficient-averaged model was taken that compares well to all three cases. Notably, the initial modulus ,u, converged across trials to an average value of 0.472 MPa with an average viscoelastic constant g of 0.561.

CONCLUSIONS

The presented model is found to have consistent parameters across loading cases. The capability of this model to represent the fundamental mechanical characteristics of the dental complex during vertical extraction loading is a significant advancement in the modelling of extraction procedures. Future work will focus on verifying the model as a predictive design tool for assessing new loading schemes in addition to investigating its applications to subject-specific problems.

Three-dimensional finite element analysis of the effect of alveolar cleft bone graft on the maxillofacial biomechanical stabilities of unilateral complete cleft lip and palate

BACKGROUND

The objective is to clarify the effect of alveolar cleft bone graft on maxillofacial biomechanical stabilities, the key areas when bone grafting and in which should be supplemented with bone graft once bone resorption occurred in UCCLP (unilateral complete cleft lip and palate).

METHODS

Maxillofacial CAD (computer aided design) models of non-bone graft and full maxilla cleft, full alveolar cleft bone graft, bone graft in other sites of the alveolar cleft were acquired by processing the UCCLP maxillofacial CT data in three-dimensional modeling software. The maxillofacial bone EQV (equivalent) stresses and bone suture EQV strains under occlusal states were obtained in the finite element analysis software.

RESULTS

Under corresponding occlusal states, the EQV stresses of maxilla, pterygoid process of sphenoid bone on the corresponding side and anterior alveolar arch on the non-cleft side were higher than other maxillofacial bones, the EQV strains of nasomaxillary, zygomaticomaxillary and pterygomaxillary suture on the corresponding side were higher than other maxillofacial bone sutures. The mean EQV strains of nasal raphe, the maximum EQV stresses of posterior alveolar arch on the non-cleft side, the mean and maximum EQV strains of nasomaxillary suture on the non-cleft side in full alveolar cleft bone graft model were all significantly lower than those in non-bone graft model. The mean EQV stresses of bilateral anterior alveolar arches, the maximum EQV stresses of maxilla and its alveolar arch on the cleft side in the model with bone graft in lower 1/3 of the alveolar cleft were significantly higher than those in full alveolar cleft bone graft model.

CONCLUSIONS

For UCCLP, bilateral maxillae, pterygoid processes of sphenoid bones and bilateral nasomaxillary, zygomaticomaxillary, pterygomaxillary sutures, anterior alveolar arch on the non-cleft side are the main occlusal load-bearing structures before and after alveolar cleft bone graft. Alveolar cleft bone graft mainly affects biomechanical stabilities of nasal raphe and posterior alveolar arch, nasomaxillary suture on the non-cleft side. The areas near nasal floor and in the middle of the alveolar cleft are the key sites when bone grafting, and should be supplemented with bone graft when the bone resorbed in these areas.

Prediction of fracture initiation and propagation in pelvic bones

The objective is developing an XFEM model that is capable of predicting different types of fracture in the pelvic bone under various loading conditions. Previously published mechanical and failure characteristics of cortical and cancellous tissues were implemented and assigned to an intact pelvic bone with specified cortical and cancellous tissues. Various loading conditions, including combined load directions, were applied to the acetabulum to model different types of fracture (e.g., anterior/posterior wall fracture and transverse fracture) in the pelvic bone. The predicated types of fracture and the maximum force at fracture were compared to those acquired from previously published experimental tests. Anterior/posterior wall fracture and transverse fracture were the most common types of fractures determined in the simulations. The XFEM simulations were able to predict similar fractures to those reported in the experimental tests. The maximum fracture force in the XFEM model was found to be 18.6 kN compared to 8.85 kN reported in the previous experimental tests. The results revealed that different types of fracture in the pelvic bones can be caused by the various loading conditions in unstable high-rate impact loads. Using proper mechanical and failure behaviors of cortical and cancellous tissues, XFEM modeling of pelvic bone is capable of predicting bone fracture. In future work, the XFEM models of cancellous and cortical tissues can be assigned to other bones in human body skeleton so that the failure mechanism in such bones can be investigated.

Advanced Modeling Methods-Applications to Bone Fracture Mechanics

PURPOSE OF REVIEW

The goal of this review is to summarize recent advances in modeling of bone fracture using fracture mechanics-based approaches at multiple length scales spanning nano- to macroscale.

RECENT FINDINGS

Despite the additional information that fracture mechanics-based models provide over strength-based ones, the application of this approach to assessing bone fracture is still somewhat limited. Macroscale fracture models of bone have demonstrated the potential of this approach in uncovering the contributions of geometry, material property variation, as well as loading mode and rate on whole bone fracture response. Cortical and cancellous microscale models of bone have advanced the understanding of individual contributions of microstructure, microarchitecture, local material properties, and material distribution on microscale fracture resistance of bone. Nano/submicroscale models have provided additional insight into the effect of specific changes in mineral, collagen, and non-collagenous proteins as well as their interaction on energy dissipation and fracture resistance at small length scales. Advanced modeling approaches based on fracture mechanics provide unique information about the underlying multiscale fracture mechanisms in bone and how these mechanisms are influenced by the structural and material constituents of bone at different length scales. Fracture mechanics-based modeling provides a powerful approach that complements experimental evaluations and advances the understanding of critical determinants of fracture risk.