Modelling of mandible bone properties in the numerical analysis of oral implant biomechanics

Anchorage loss of the posterior teeth under different extraction patterns in maxillary and mandibular arches using clear aligner: a finite element study

BACKGROUND

Extracting the premolars is an effective strategy for patients with bimaxillary dentoalveolar protrusion. Clear aligners (CAs) close the extraction spaces through shortening the length of aligners. The contraction force generated by the terminal of aligners makes the posterior teeth tip mesially, which is known as the roller coaster effect. This phenomenon is even worse in the 2nd premolar extraction cases. Posterior anchorage preparation is commonly used to protect the angulation of molars, taking the form of presetting distal tipping value. However, the distal tipping design aggravates the anchorage loss of anterior teeth simultaneously. This study aimed to explore the different anchorage loss of the posterior teeth when the 1st or 2nd premolars were extracted using CAs, respectively in maxillary and mandibular arches, further providing guidance for anchorage preparation design in clinical practice.

METHODS

Two bimaxillary finite element models with different extraction patterns were established to simulate the anterior en-masse retraction process of the CAs. In Model 1, the maxillary and mandibular 1st premolars were extracted, while in Model 2, the 2nd premolars were extracted. Finite element analysis methods were utilized to analyze the tipping angle of the anterior and posterior teeth.

RESULTS

Compared between two models, the anterior teeth exhibited a greater lingual inclination tendency and the posterior teeth exhibited a slighter mesial tipping tendency in Model 1 regarding individual tooth. The closer to the extraction spaces, the greater the tip, and the distal tipping tendency of the 1st premolars was more evident than the mesial tipping tendency of the 1st molars in Model 2. Compared between the maxillary and mandibular arches, the mesial tipping tendency of individual posterior tooth was more evident in the maxilla. In addition, the highest hydrostatic stress of the periodontal ligaments was concentrated on the cervical and apical parts directly adjacent to the extraction spaces, and it exhibited relatively uniform distribution in Model 1.

CONCLUSIONS

The individual posterior tooth showed the same mesial tipping direction but to different degree when the 1st or the 2nd premolars were extracted during clear aligner treatment. Presetting anchorage preparation design for the posterior teeth is essential to alleviate the roller coaster effect, especially in the 2nd premolar extraction cases. Furthermore, larger anchorage preparation value should be proposed for the maxillary posterior teeth.

Comparison of Bulk Polymeric Resin Composite and Hybrid Glass Ionomer Cement in Adhesive Class I Dental Restorations: A 3D Finite Element Analysis

This study aimed to investigate the mechanical behavior of resin composites and hybrid glass ionomer cement in class I adhesive dental restorations under loading and shrinkage conditions. Three CAD models of a mandibular first molar with class I cavities were created and restored with different techniques: a bi-layer of Equia Forte HT with Filtek One Bulk Fill Restorative composite (model A), a single layer of adhesive and Filtek One Bulk Fill Restorative (model B), and a single layer of Equia forte HT (model C). Each model was exported to computer-aided engineering software, and 3D finite element models were created. Models A and B exhibited a similar pattern of stress distribution along the enamel-restoration interface, with stress peaks of 12.5 MPa and 14 MPa observed in the enamel tissue. The sound tooth, B, and C models showed a similar trend along the interface between dentine and restoration. A stress peak of about 0.5 MPa was detected in the enamel of both the sound tooth and B models. Model C showed a reduced stress peak of about 1.2 MPa. A significant stress reduction in 4 mm deep class I cavities in lower molars was observed in models where non-shrinking dental filling materials, like the hybrid glass ionomer cement used in model C, were applied. Stress reduction was also achieved in model A, which employed a bi-layer technique with a shrinking polymeric filling material (bulk resin composite). Model C's performance closely resembled that of a sound tooth.

Optimizing restorative procedure and material selection for pulpotomized primary molars: Mechanical characterization by 3D finite element analysis

Purpose

This study aimed to assess the stress distribution in pulpotomized primary molars with different types of restorative materials using 3D-finite element analysis (FEA), and provide valuable insights into the selection and application of restorative materials, with the ultimate goal of reducing the risk of pulpotomy failure and protecting residual dental tissue.

Methods

Four 3D models of pulpotomized primary molars with different restorative materials according to the material and its elastic modulus were analysed: resin composite, stainless steel crowns (SSCs), prefabricated zirconia crowns and endocrowns. The food layer was also designed before vertical and bucco-lingual forces were applied to simulate physiological masticatory conditions. The results were obtained by colorimetric graphs of the von Mises stresses (VMS) in the restoration and tooth remnant. The maximum shear stress on the bonding interfaces and pressure stress on the Mineral trioxide aggregate (MTA)-pulp interfaces were recorded.

Results

The results of the 3D-FEA showed that all restorative materials generated stresses and strains on the tooth structure after pulpotomy. In the resin composite group, the marginal enamel exhibited the highest stress peaks. In the zirconia crown and SSC groups, there was a concentration of stress at the dentin-restoration margin. The shear stress concentrations were mainly at the adhesive margins, with lower levels around endocrowns compared to other groups. MTA in the resin composite group experienced more VMS than in the other group. The resin composite group also generated relatively higher pressure stress values at the MTA-pulp interface compared to the other groups.

Significance

In the model of primary teeth following pulpotomy, the three types of restorations covering the occlusal surface can effectively reduce the stress on pulp capping materials under occlusal loads, thereby potentially decreasing the risk of pulpotomy failure. In addition, the group of endocrowns demonstrated reduced stress at the bonding interface and in the stress concentration zone near the dentist-restoration edge, making them more effective at protecting residual dental tissue.

Rigid flexible coupling contact mechanism for oral and maxillofacial skin and soft tissues

BACKGROUND AND OBJECTIVES

The existing medical clinical treatment institutions mostly use rigid structures to come into contact with flexible skin. The rigid flexible coupled contact biomechanical model for the skin is the first step that urgently needs to be considered in the process of medical clinical operations. However, there has been currently no effective biomechanical contact model available.

METHODS

Based on the principle of elastic interface deformation, the basic biomechanical characteristics of oral and maxillofacial skin and soft tissues were analyzed to address the unknown mechanism of rigid body and maxillofacial contact in oral imaging operations. A nonlinear characterization method for the mechanical properties of oral and maxillofacial skin soft tissues was proposed by deriving a general contact force model that takes into account energy dissipation. However, the problem of the inability to obtain analytical solutions for the parameters of the dynamic model exists. It is necessary to perform particle swarm parameter identification on different nonlinear contact models and verify the accuracy of the algorithm through numerical simulation. A maxillofacial contact experiment was conducted to verify the operation process of an oral imaging robot.

RESULTS

After experimental analysis, it was found that the comprehensive average error between the model and the actual contact force was 0.13325 N. The absolute error of the maximum deformation displacement was below 0.18 N, which verified the effectiveness and safety of the contact model in the contact process of the oral imaging robot system.

CONCLUSIONS

The results indicate that the output force of the model has been in good agreement with the actual contact force.

Biomechanical Modelling for Tooth Survival Studies: Mechanical Properties, Loads and Boundary Conditions-A Narrative Review

Recent biomechanical studies have focused on studying the response of teeth before and after different treatments under functional and parafunctional loads. These studies often involve experimental and/or finite element analysis (FEA). Current loading and boundary conditions may not entirely represent the real condition of the tooth in clinical situations. The importance of homogenizing both sample characterization and boundary conditions definition for future dental biomechanical studies is highlighted. The mechanical properties of dental structural tissues are presented, along with the effect of functional and parafunctional loads and other environmental and biological parameters that may influence tooth survival. A range of values for Young's modulus, Poisson ratio, compressive strength, threshold stress intensity factor and fracture toughness are provided for enamel and dentin; as well as Young's modulus and Poisson ratio for the PDL, trabecular and cortical bone. Angles, loading magnitude and frequency are provided for functional and parafunctional loads. The environmental and physiological conditions (age, gender, tooth, humidity, etc.), that may influence tooth survival are also discussed. Oversimplifications of biomechanical models could end up in results that divert from the natural behavior of teeth. Experimental validation models with close-to-reality boundary conditions should be developed to compare the validity of simplified models.

Effects of assessing the bone remodeling process in biomechanical finite element stability evaluations of dental implants

BACKGROUND AND OBJECTIVE

While an accurate assessment of the biomechanical stability of implants is essential in dental prosthesis planning and associated treatment assurance, the bone remodeling process is often ignored in biomechanical studies using finite element (FE) analysis. In this study, we aimed to analyze the significance of assessing the bone remodeling process in FE analysis for evaluating the biomechanical stability of dental implants. We compared the FE results considering the bone remodeling process with FE results simulated using commonly used conditions, with no considerations of the bone remodeling process.

METHODS

The mathematical model proposed by Komarova et al. was used to calculate cell population dynamics and changes in bone density at a discrete site. The model was implemented in the FE software ABAQUS, using the UMAT subroutine. Three-dimensional FE models were constructed for two types of bone (III and IV) and three values of implant diameter (4.0, 4.5, and 5.0 mm). An average biting force of 50 N in the vertical direction was applied during the bone remodeling process for 150 days. Afterwards, the maximum biting force of 200 N in the 30° oblique direction was applied to evaluate the stability of the implant systems.

RESULTS

To understand the impact of bone remodeling on the resultant mechanical responses, we focused on peri-implant cancellous bone based on two parameters: apparent density change and microstrain distribution. The bone density decreased by an average of 5.3 % after implantation, and it was the lowest on the 6^th day. The average density increases of the peri-implant cancellous bone were 264.4 kgm³ (bone type III) and 220.0 kgm³ (bone type IV) over 150 days. For the bone stability analysis, the maximum principal strain in the peri-implant bone was used to evaluate the bone stability. If the bone remodeling process is ignored, then the bone volume within the fatigue failure range of the microstrain differs significantly from that if the bone remodeling process is considered, i.e., 60 % higher for bone type III and 33.4 % lower for bone type IV than when the bone remodeling process is considered.

CONCLUSIONS

The FE result without considering the bone remodeling process could be considered a conservative criterion for bone type III. However, in bone type IV, the FE result without considering the bone remodeling process tends to underestimate the risks. The bone remodeling process is more affected by the initial bone quality than the implant diameter.

Evaluation of Viability and Cell Proliferation in Bone and Gingival on Dental Implant Fixtures with Active Sandblasted and Sandblasted Surfaces by the Cytotoxicity Test Method

In recent years, the technology of dental implants has advanced a lot and this has become one of the main reasons for dentists to refer to implants immediately after tooth extraction. Evaluation of cell survival in implantology studies to determine cell sensitivity determines the outcome of treatment. This study aimed to investigate the bone integration properties as well as the cytotoxicity between the implant surface and the jaw bone. In this study, cytotoxicity test was used to evaluate the toxicity and viability of sandblasted large grit acid etched, and sandblasted large grit acid etched active surfaces in 3A brand dental implants with G292 osteoblastic cells and human gingival fibroblasts cells are discussed. This operation was performed using a laboratory incubator of the German company MEMMERT for 24 hours, by neubauer lam cells counting for one hundred thousand cultured cells in each test at a temperature of 37 °C, a pressure of 1 atmosphere and 90% humidity. Based on the scanning electron microscopy images and the cytotoxicity test results, it can be seen that the bone graft of the implant, with the sandblasted large grit acid etched active surface treatment, is much better and also one week faster than the implant with the sandblasted large grit acid etched one. However, the viability of the implant with the sandblasted large grit acid etched active surface treatment for both G292 osteoblastic cells and human gingival fibroblasts cells samples was equal to 98.4% and 97.3%, respectively, and is lower than the sandblasted large grit acid etched surface treatment. The results show that the viability of the sandblasted large grit acid etched implant is about 1.5 to 2% higher than the sandblasted large grit acid etched active one, but the surface integrity of sandblasted large grit acid etched active is better than sandblasted large grit acid etched in all samples, and the treatment process is reduced by one week.

Influence of restorative material and cement on the stress distribution of endocrowns: 3D finite element analysis

Purpose

This study aimed to evaluate the influence of different types of restorative materials and resin cements on the stress distribution in the regions of the restoration, cement layer and dental remnant in endodontically treated posterior endocrowns.

Methods

A 3D finite element analysis (FEA) model of the first mandibular molar that was restored with an endocrown designed by computer-aided design (CAD) software was generated. Three kinds of restorative materials (Vita Enamic (VE), IPS e.max CAD (EMX) and Grandio blocs (GR)) and two types of cementing materials (NX3 and Maxcem Elite Chroma (MX)) were analysed with such a model. The food layer was also designed before vertical (600 N) forces were applied to simulate physiological masticatory conditions. Thermal expansion was used to simulate the polymerization shrinkage effects of cement layers. The results were obtained by colorimetric graphs of the maximum principal stress in the restoration and tooth remnant. The failure risk of the cement layer was also calculated based on the normal stress.

Results

The elastic modulus was positively correlated with the tensile stress peak values in the restoration, mainly at the intaglio surface. However, in the cervical enamel and cement layer, restorative material with a higher elastic modulus generated lower peak stress values. The cement with a higher elastic modulus resulted in higher stress peak values inside the cement layer. The combination of EMX (restorative material) and NX3 (cement material) in the cement layer resulted in the lowest failure risk.

Significance

The ceramic material EMX with a higher elastic modulus appeared to be more effective at protecting the cement layer and residual enamel tissue. Based on the analysis of the failure risk of the cement layer, the combination of EMX and NX3 was recommended as an optional material for endocrowns for endodontically treated posterior teeth.

Design of dental implant using design of experiment and topology optimization: A finite element analysis study

Ever since the introduction of topology optimization into the industrial and manufacturing fields, it has been a top priority to maximize the performance of any system by optimizing its geometrical parameters to save material while keeping its functionality unaltered. The purpose of this study is to design a dental implant macro-geometry by removing expendable material using topology optimization and to evaluate its biomechanical function. Three-dimensional finite element models were created of an implant embedded in cortical and cancellous bone. Parameters like the length and diameter of the implant and the bone quality (±20% variation in Young’s modulus, Poisson’s ratio and density for both cortical and cancellous bone) were varied to evaluate their effect on the principal stresses induced on the peri-implant bone tissues and the micromotion of the implant at 150 N applied load. Design optimization is used to select one suitable implant for each material property combination with optimum parameters that experiences the least von Mises stress and axial deformation, out of twenty implants with different length and diameter for each material property combination. Topology optimization was then used on the selected implants to remove the redundant material. The biomechanical functions of the implants with optimized parameter and volume were then evaluated. The finite element analyses estimated that a reduction of 32% to 45% in the implant volume is possible with the implant still retaining all of its functionality.

The use of different adhesive filling material and mass combinations to restore class II cavities under loading and shrinkage effects: a 3D-FEA

3D tooth models were virtually restored: flowable composite resin + bulk-fill composite (A), glass ionomer cement + bulk-fill composite (B) or adhesive + bulk-fill composite (C). Polymerization shrinkage and masticatory loads were simulated. All models exhibited the highest stress concentration at the enamel-restoration interfaces. A and C showed similar pattern with lower magnitude in A in comparison to C. B showed lower stress in dentine and C the highest cusps displacement. The use of glass ionomer cement or flowable composite resin in combination with a bulk-fill composite improved the biomechanical behavior of deep class II MO cavities.

Adhesive class I restorations in sound molar teeth incorporating combined resin-composite and glass ionomer materials: CAD-FE modeling and analysis

OBJECTIVES

To investigate the influence of different resin composite and glass ionomer cement material combinations in a "bi-layer" versus a "single-layer" adhesive technique for class I cavity restorations in molars using numerical finite element analysis (FEA).

MATERIALS AND METHODS

Three virtual restored lower molar models with class I cavities 4mm deep were created from a sound molar CAD model. A combination of an adhesive and flowable composite with bulk fill composite (model A), of a glass ionomer cement with bulk fill composite (model B) and of an adhesive with bulk fill composite (model C), were considered. Starting from CAD models, 3D-finite element (FE) models were created and analyzed. Solid food was modeled on the occlusal surface and slide-type contact elements were used between tooth surface and food. Polymerization shrinkage was simulated for the composite materials. Physiological masticatory loads were applied to these systems combined with shrinkage. Static linear analyses were carried out. The maximum normal stress criterion was adopted as a measure of potential damage.

RESULTS

All models exhibited high stresses principally located along the tooth tissues-restoration interfaces. All models showed a similar stress trend along enamel-restoration interface, where stresses up to 22MPa and 19MPa was recorded in the enamel and restoration, respectively. A and C models showed a similar stress trend along the dentin-restoration interface with a lower stress level in model A, where stresses up to 11.5MPa and 7.5MPa were recorded in the dentin and restoration, respectively, whereas stresses of 17MPa and 9MPa were detected for model C. In contrast to A and C models, the model B showed a reduced stress level in dentin, in the lower restoration layer and no stress on the cavity floor.

SIGNIFICANCE

FE analysis supported the positive effect of a "bi-layer" restorative technique in a 4mm deep class I cavities in lower molars versus "single-layer" bulk fill composite technique.

Periprosthetic biomechanical response towards dental implants, with functional gradation, for single/multiple dental loss

The differences in shape and stiffness of the dental implants with respect to the natural teeth (especially, dental roots) cause a significant alteration of the periprosthetic biomechanical response, which typically leads to bone resorption and ultimately implant loosening. In order to avoid such clinical complications, the implant stiffness needs to be appropriately adapted. In this study, hollow channels were virtually introduced within the designed implant screws for reduction of the overall stiffness of the prototype. In particular, two opposing radial gradients of increasing hollow channel diameters, i.e., outside to inside (Channel 1) and inside to outside (Channel 2) were considered. Two clinical situations of edentulism were addressed in this finite element-based study, and these include a) loss of the first molar, and b) loss of all the three molars. Consequently, two implantation approaches were simulated for multiple teeth loss - individual implantation and implant supported dental bridge. The effects of implant length, approach and channel distribution on the biomechanical response were evaluated in terms of the von Mises stress within the interfacial periprosthetic bone, under normal masticatory loading. The results of our FE analysis clearly reveal significant variation in periprosthetic bone stress between the different implant designs and approaches. An implant screw length of 11 mm with the Channel 2 configuration was found to provide the best biomechanical response. This study also revealed that the implant supported dental bridge approach, which requires lower bone invasion, results in favorable biomechanical response in case of consecutive multiple dental loss.

Influences of geometrical and mechanical properties of bone tissues in mandible behaviour - experimental and numerical predictions

The properties and geometry of bone in the mandible play a key role in mandible behaviour during a person's lifetime, and attention needs to be paid to the influence of bone properties. We analysed the effect of bone geometry, size and bone properties in mandible behaviour, experimenting on cadaveric mandibles and FE models. The study was developed using the geometry of a cadaveric mandible without teeth. Three models of cadaveric condyles were experimentally tested with instrumented with four rosettes, and a condyle reaction of 300 N. Four finite element models were considered to validate the experiments and analyse mandible behaviour. One numeric model was simulated with 10 muscles in a quasi-static condition. The experimental results present different condyle stiffness's, of 448, 215 and 254 N/mm. The values presented in the rosettes are influenced by bone geometry and bone thickness; maximum value was -600 με in rosette #4, and the maximum strain difference between mandibles was 111%. The numerical results show that bone density decreases and strain distribution increases in the thinner mandible regions. Nevertheless, the global behaviour of the structure remains similar, but presents different strain magnitudes. The study shows the need to take into account bone characteristics and their evolutions in order to improve implant design and fixation throughout the patient life. The change in bone stiffness promotes a change in maximum strain distribution with same global behaviour.

Investigation of the biomechanical behaviour of hindfoot ligaments

The aim of this work is to provide a computational tool for the mechanical characterization of the hindfoot ligaments. The investigation is performed by a coupled numerical and experimental approach. For this purpose, a numerical model that represents the complex structural configuration of the hindfoot and the typical features of the mechanical behaviour of the ligament tissue is developed. The geometrical analysis of the anatomical site is performed starting from the processing of computed tomography and magnetic resonance images. Accounting for morphometric measurements, the virtual solid model provides an averaged configuration of the hindfoot structure. In order to specify the mechanical behaviour of the ligament tissue, a fibre-reinforced visco-hyperelastic model is adopted. The formulation accounts for the anisotropic configuration, geometric non-linearity, non-linear elasticity and time-dependent phenomena. Numerical analyses are performed to evaluate the biological tissues and structure mechanics with regard to physiological boundary conditions, accounting for dorsiflexion and plantarflexion movements. In order to evaluate the reliability of the numerical model developed, the experimental data are compared with the numerical results. The numerical results are in agreement with the range of values obtained by experimental test confirming the accuracy of the procedure adopted.

Investigations on the viscoelastic behaviour of a human healthy heel pad: In vivo compression tests and numerical analysis

The aim of this study was to investigate the viscoelastic behaviour of the human heel pad by comparing the stress–relaxation curves obtained from a compression device used on an in vivo heel pad with those obtained from a three-dimensional computer-based subject-specific heel pad model subjected to external compression. The three-dimensional model was based on the anatomy revealed by magnetic resonance imaging of a 31-year-old healthy female. The calcaneal fat pad tissue was described with a viscohyperelastic model, while a fibre-reinforced hyperelastic model was formulated for the skin. All numerical analyses were performed to interpret the mechanical response of heel tissues, with loading conditions and displacement rate in agreement with experimental tests. The heel tissues showed a non-linear, viscoelastic behaviour described by characteristic hysteretic curves, stress–relaxation and viscous recovery phenomena. The reliability of the investigations was validated by the interpretation of the mechanical response of heel tissues under the application of three pistons with diameter of 15, 20 and 40 mm, at the same displacement rate of about 1.7 mm/s. The maximum and minimum relative errors were found to be less than 0.95 and 0.064, respectively.