Modified bone density-dependent orthotropic material model of human mandibular bone

Heterogeneous material models for finite element analysis of the human mandible bone - A systematic review

Subject-specific finite element (FE) modeling of the mandible bone has recently gained attention for its higher accuracy. A critical modeling factor is including personalized material properties from medical images especially when bone quality has to be respected. However, there is no consensus on the material model for the mandible that realistically estimates the Young's modulus of the bone. Therefore, the present study aims to review FE studies employing heterogeneous material modeling of the human mandible bone, synthesizing these models, investigating their origins, and assessing their risk of bias. A systematic review using PRISMA guidelines was conducted on publications before 1st July 2024, involving PubMed, Scopus, and Web of Science. The search string considered (a) anatomical site (b) modeling strategy, and (c) metrics of interest. Two inclusion and five exclusion criteria were defined. A review of 77 FE studies identified 12 distinct heterogeneous material models, built based on different in vitro or computational methodologies leading to varied performance and highly deviated range of estimated Young's modulus. They are proposed for bones from five different anatomical sites than mandible and for both trabecular and cortical bone domains. The original studies were characterized with a low to medium risk of bias. This review assessed the current state of material modeling for subject-specific FE models in the craniomaxillofacial field. Recommendations are provided to support researchers in selecting density-modulus relationships. Future research should focus on standardizing experimental protocols, validating models through combined simulation and experimental approaches, and investigating the anisotropic behaviour of the mandible bone.

Homogenized finite element simulations can predict the primary stability of dental implants in human jawbone

Adequate primary stability is a pre-requisite for the osseointegration and long-term success of dental implants. Primary stability depends essentially on the bone mechanical integrity at the implantation site. Clinically, a qualitative evaluation can be made on medical images, but finite element (FE) simulations can assess the primary stability of a bone-implant construct quantitatively based on high-resolution CT images. However, FE models lack experimental validation on clinically relevant bone anatomy. The aim of this study is to validate such an FE model on human jawbones. Forty-seven bone biopsies were extracted from human cadaveric jawbones. Dental implants of two sizes (Ø3.5 mm and Ø4.0 mm) were inserted and the constructs were subjected to a quasi-static bending-compression loading protocol. Those mechanical tests were replicated with sample-specific non-linear homogenized FE models. Bone was modeled with an elastoplastic constitutive law that included damage. Density-based material properties were mapped based on μCT images of the bone samples. The experimental ultimate load was better predicted by FE (R2 = 0.83) than by peri-implant bone density (R2 = 0.54). Unlike bone density, the simulations were also able to capture the effect of implant diameter. The primary stability of a dental implant in human jawbones can be predicted quantitatively with FE simulations. This method may be used for improving the design and insertion protocols of dental implants.

Finite element analysis in implant dentistry: State of the art and future directions

OBJECTIVE

To discuss the state of the art of Finite Element (FE) modeling in implant dentistry, to highlight the principal features and the current limitations, and giving recommendations to pave the way for future studies.

METHODS

The articles' search was performed through PubMed, Web of Science, Scopus, Science Direct, and Google Scholar using specific keywords. The articles were selected based on the inclusion and exclusion criteria, after title, abstract and full-text evaluation. A total of 147 studies were included in this review.

RESULTS

To date, the FE analysis of the bone-dental implant system has been investigated by analyzing several types of implants; modeling only a portion of bone considered as isotropic material, despite its anisotropic behavior; assuming in most cases complete osseointegration; considering compressive or oblique forces acting on the implant; neglecting muscle forces and the bone remodeling process. Finally, there is no standardized approach for FE modeling in the dentistry field.

SIGNIFICANCE

FE modeling is an effective computational tool to investigate the long-term stability of implants. The ultimate aim is to transfer such technology into clinical practice to help dentists in the diagnostic and therapeutic phases. To do this, future research should deeply investigate the loading influence on the bone-implant complex at a microscale level. This is a key factor still not adequately studied. Thus, a multiscale model could be useful, allowing to account for this information through multiple length scales. It could help to obtain information about the relationship among implant design, distribution of bone stress, and bone growth. Finally, the adoption of a standardized approach will be necessary, in order to make FE modeling highly predictive of the implant's long-term stability.

Influence of surface coatings on the stress distribution by varying friction contact at implant-bone interface using finite element analysis

The present work aims to evaluate the effect of various surface coatings of titanium dental implants by varying the friction coefficient (µ) at the interface between the dental implant and jawbone using finite element analysis (FEA) methods and to provide a comparative analysis between the various surface coatings and implant designs. An accurate model of the dental implant prosthetics consisting of the hard (cortical) and the soft (cancellous) bone, with the various titanium dental implant designs was modelled using a 3D CAD software, and the FE mesh model was generated using HyperMesh 13.0. Three coatings having different coefficient of friction values were selected: Titanium Nitride (TiN) with a friction coefficient of 0.19, Titanium Oxide (TiO₂) with a friction coefficient of 0.30 and Zirconium Nitride (ZrN) with a coefficient of friction of 0.49. The non-linear static stress analysis was conducted under three different loading conditions (vertical, lateral and oblique loading) using a CAE solver. The present study showed that surface coatings with high friction coefficients generated lower stresses in the cancellous bone while generating higher stresses in the cortical bone. However, for dental implants having microthreads in their neck region, surface coatings with a high coefficient of friction generated lower stresses at the interface between the cortical bone and the implant. The FEA results indicate that selecting suitable surface coatings would significantly decrease the stresses developed at the bone-implant interface, and future studies should conduct in vivo trials to validate the FEA results obtained.

Variation in stress distribution modified by mandibular material property: a 3D finite element analysis

BACKGROUND

Temporomandibular joint disorder (TMD) is a common oral and maxillary facial disease. Finite element method (FEM) has been widely used in TMD studies. Material assignment significantly affects FEM results. The differences in the methods of material assignment used in previous studies have not been comprehensively assessed for further calculations.

METHODS

The mandible material modelling approaches were of four types, namely: uniform modelling with (A) cortical bone; and (B) cancellous bone; (C) semi-uniform modelling with division of cortical and cancellous bone; and (D) non-uniform modelling with Computed tomography (CT) gray value related modulus. Meanwhile, the Young's modulus of values ranging from 20 to 300 GPa were considered for the teeth. Ten modellings were used to analyze and discuss the differences in contact pressure and contact force.

RESULTS

(1) The increase in teeth elastic modulus increased the maximum contact pressure on the alveolar bone and contact force on teeth, but induced insignificant stress variation on the temporomandibular joint; (2) The location of the maximum contact pressure was steady for all four modelling approaches of the mandibular material. However, the maximum contact pressure and contact force exhibited an insignificant difference.

CONCLUSIONS

Teeth with a higher elastic modulus significantly enhanced the stress concentration in the alveolar bone; in contrast, it induced minor variations in the temporomandibular joint stress states. The extreme stress regions predicted by the four mandibular models were consistent with the actual damaged regions. However, non-uniform modellings based on CT values could better describe the mechanical properties of the human bone, which should be primarily considered.

An approach for simultaneous reduction and fixation of mandibular fractures

This article presents a new approach for the design of a flexible V-shaped miniplate for mandibular fractures, which combines simultaneous fracture reduction and fixation. A Computerized Tomography (CT) based finite element model was developed to assess the reliability of this design. Muscle and mastication forces were included to replicate post-surgery loading. The V-plate is compared with a standard, linear miniplate, typically used for mandibular fixation. The results indicate that the proposed design can support the fracture while inducing limited fracture displacement, in addition to reducing the duration of the surgery due to fracture reduction by tightening the wire.

Mapping cortical bone stiffness and mineralization from endosteal to periosteal surfaces of bovine mid-diaphyseal femur

INTRODUCTION

While bone literature abounds with correlations of mechanical stiffness to mineralization, such correlations are reported without relating the findings to specific intracortical locations. This study reports on mapping of stiffness and mineralization distributions in ring-shaped cortical bone samples sliced from mid-diaphyseal bovine femur. Stiffness and mineralization measurements were conducted at points across the intracortical thickness along radial lines emanating from the inner (endosteal) surface to the outer (periosteal) surface. Measurements were taken along approximately 4 mm distance of cortical bone thickness.

MATERIALS AND METHODS

Three experimental techniques were employed: Vickers microhardness (HV), energy-dispersive X-ray (EDX) spectroscopy, and computed tomography (CT). Stiffness values were extracted from the Vickers microhardness tests. Elemental mineralization values (calcium %wt. and phosphorus %wt.) were determined from EDX data. All measurements were repeated on three different femur bones taken from different bovines (collected fresh from butcher).

RESULTS

The study plots stiffness values and elemental mineralization (calcium %wt. and phosphorus %wt.) versus cortical thickness. Both stiffness and Ca %wt. and P %wt. are found to track and to linearly increase when plotted along the radial distance. The stiffness and mineralization trends collected from Vickers and EDX measurements were verified by employing the CT number (Hounsfield units, HU) via CT scans of the same bone samples. Data fitting via statistical methods revealed that all correlations were statistically significant.

CONCLUSION

Starting from endosteal to periosteal surfaces of mid-diaphyseal bovine femur, it was found that stiffness, mineralization, and HU values all exhibit increasing and correlating trends.

Design and Development of a Novel Oral Care Simulator for the Training of Nurses

OBJECTIVE

A Novel Oral Care Simulator was designed and developed to measure and visualise the facial and lingual forces exerted on teeth by the action of tooth brushing, considering the irregular geometry and structural composition of human dentition and the emulation of the realistic biomechanical deflection of the teeth.

METHOD

FEA simulations were carried out on a central incisor under facial loading and an appropriate force sensing mechanism was designed. An anatomically accurate mandibular jaw and 16 teeth were 3D printed, on which 16 force sensing structures were embedded. The signals from the sensors were amplified using a multichannel signal amplifier built using instrumentation amplifiers which were then visualised through a GUI.

RESULTS

The developed simulator is capable of indicating the magnitude of a force upto 15 N exerted on to the facial and lingual surfaces of teeth at a frequency of 60 Hz and above and it is capable of alerting the user if the force exceeds a pre-specified threshold.

CONCLUSION

The designed force sensing mechanism considers the irregular geometry and structural composition of human dentition in measuring the facial and lingual forces. It provides a reliable feedback by indicating the force and emulating the realistic biomechanical deflection of teeth.

SIGNIFICANCE

Nurses who care for the disabled, elderly and sick have explicitly stated the requirement for a simulator to train themselves on brushing the teeth of their subjects as their incorrect technique can cause longterm dental damage, for which a device has not been developed to date.

Implant-supported overdentures with different clinical configurations: Mechanical resistance using a numerical approach

STATEMENT OF PROBLEM

Implant-supported overdentures (IODs) are a treatment option for patients with complete edentulism. However, this treatment increases the possibilities of peri-implant complications, characterized by inflammation or partial loss of surrounding hard and soft tissues.

PURPOSE

The purpose of this finite element analysis study was to evaluate the mechanical performance of different bar-IOD designs under different clinical configurations by comparing the stress and strain distribution on the bone during secondary stabilization.

MATERIAL AND METHODS

A finite element model of the mandible representing a patient with complete edentulism was developed. Different designs of bar-IODs were modeled and compared. The parameters studied were the material properties (cobalt-chromium, zirconium dioxide, titanium grade 5, and titanium grade 4), diameter and bar-IOD cross-sectional shape, tilt of the posterior implants (30 degrees), presence of a distal extension cantilever in the bar-IODs (12 mm), and number of implants (4 or 6). Two different mastication loading conditions were analyzed. One- and 2-way ANOVAs and the Tukey honestly significant differences post hoc test (α=.05) were used to determine the significant von Mises stress and strain values in the bone.

RESULTS

The 4 materials tested in the bar-IOD did not have a significant mechanical effect on the bone (P<.05). A smaller diameter and structure of the bar-IOD led to significantly higher bone stress (P<.001). A distal extension cantilever led to an increased stress concentration (model M1 versus model M3: P<.001), which reached 50% in the event of tilting of the posterior implants (model M2 versus model M4: P<.001). Tilting of the posterior implants alone, without extension, had a nonsignificant effect (model M3 versus model M4: P=.999). Model M5 supported with 6 implants reduces the stress transferred to the bone compared with model M3 supported with 4 implants (P<.05).

CONCLUSIONS

Distal extensions in bar-IODs, the tilt of the posterior implants, and the low amount of material in the cross-sectional area in the bar-IOD were the most influential parameters on the mechanical resistance of dental implants in the mandibular bone.

Stress Distribution Patterns within Viscero- and Neurocranium during Nasoalveolar Molding: a Finite Element Analysis

Background

The purpose of this study was to evaluate the stress distribution patterns within the viscero- and neurocranium of neonates during nasoalveolar molding.

Methods

Finite element models of 3 different healthy neonates at different times of life (date of birth, 4 weeks, and 3.5 months) were generated on the basis of computed tomography scans. A validated workflow, including segmentation, meshing, setting of boundary conditions, and implementation of a bone density-dependent material model, was carried out for each model. A small and a large unilateral alveolar and hard palatal cleft were virtually cut in each model. The stress distribution pattern in each model was then analyzed by using Ansys APDL.

Results

Convergence analysis validated the results. The virtual experiments at the date of birth showed a stress pattern above a previously defined threshold value of 30,000 Pa in the ipsilateral naso-orbital-complex, frontal sinus, and the anterior fossa of the base of the skull, with von Mises values > 35,000 Pa. Stress patterns at the age of 4 weeks and 3.5 months showed reduced von Mises values at < 15,000 Pa.

Conclusions

Nasoalveolar molding therapy is a safe presurgical treatment modality without significant influence on the viscero- and neurocranium of neonates. Treatment, considering the stress distribution at the naso-orbital-complex and anterior fossa of the base of the skull, should begin in the second week of life, and treatment initiation of preterm infants should be adapted respectively.

Finite element analysis of dental implants with validation: to what extent can we expect the model to predict biological phenomena? A literature review and proposal for classification of a validation process

A literature review of finite element analysis (FEA) studies of dental implants with their model validation process was performed to establish the criteria for evaluating validation methods with respect to their similarity to biological behavior. An electronic literature search of PubMed was conducted up to January 2017 using the Medical Subject Headings “dental implants” and “finite element analysis.” After accessing the full texts, the context of each article was searched using the words “valid” and “validation” and articles in which these words appeared were read to determine whether they met the inclusion criteria for the review. Of 601 articles published from 1997 to 2016, 48 that met the eligibility criteria were selected. The articles were categorized according to their validation method as follows: in vivo experiments in humans (n = 1) and other animals (n = 3), model experiments (n = 32), others’ clinical data and past literature (n = 9), and other software (n = 2). Validation techniques with a high level of sufficiency and efficiency are still rare in FEA studies of dental implants. High-level validation, especially using in vivo experiments tied to an accurate finite element method, needs to become an established part of FEA studies. The recognition of a validation process should be considered when judging the practicality of an FEA study.

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.