Finite element contact analysis as a critical technique in dental biomechanics: A review

Biomechanical impact of different isthmus positions in mandibular first molar root canals: a finite element analysis

OBJECTIVE

This study used image-based finite element analysis (FEA) to assess the biomechanical changes in mandibular first molars resulting from alterations in the position of the root canal isthmus.

METHODS

A healthy mandibular first molar, characterized by two intact root canals and a cavity-free surface, was selected as the subject. A three-dimensional model for the molar was established using scanned images of the patient's mandibular teeth. Subsequently, four distinct finite element models were created, each representing varied root canal morphologies: non-isthmus (Group A), isthmus located at the upper 1/3 of the root (Group B), middle 1/3 of the root (Group C), and lower 1/3 of the root (Group D). A static load of 200 N was applied along the tooth's longitudinal axis on the occlusal surface to simulate regular chewing forces. The biomechanical assessment was conducted regarding the mechanical stress profile within the root dentin. The equivalent stress (Von Mises stress) was used to assess the biomechanical features of mandibular teeth under mechanical loading.

RESULTS

In Group A (without an isthmus), the maximum stress was 22.2 MPa, while experimental groups with an isthmus exhibited higher stresses, reaching up to 29.4 MPa. All maximum stresses were concentrated near the apical foramen. The presence of the isthmus modified the stress distribution in the dentin wall of the tooth canal. Notably, dentin stresses at specific locations demonstrated differences: at 8 mm from the root tip, Group B: 13.6 MPa vs. Group A: 11.4 MPa; at 3 mm from the root tip, Group C: 14.2 MPa vs. Group A: 4.5 MPa; at 1 mm from the root tip, Group D: 25.1 MPa vs. Group A: 10.3 MPa. The maximum stress in the root canal dentin within the isthmus region was located either at the top or bottom of the isthmus.

CONCLUSION

A root canal isthmus modifies the stress profile within the dentin. The maximum stress occurs near the apical foramen and significantly increases when the isthmus is located closer to the apical foramina.

3D finite element analysis of stress distribution on the shape of resected root-end or with/without bone graft of a maxillary premolar during endodontic microsurgery

Background/purpose

Apical root resection pattern affects the stress distribution behavior in the apical region of the resected tooth. The purpose of the study was to compare the biomechanical responses of resected teeth between endodontic microsurgery (horizontal resection) and targeted endodontic microsurgery (round resection).

Materials and methods

Five different models were developed. The basic model without resection (NR) was regarded as the control model, and the others involved: horizontal resection without bone grafting (HN), horizontal resection with bone grafting (HG), round resection without bone grafting (RN), and round resection with bone grafting (RG) models. A static load of 100 N was applied to the buccal and palatal cusps of all the teeth in a 30° oblique direction. The maximum von-Mises stress and tooth displacement values were analyzed and compared.

Results

Both the HN and RN models exhibited lower stress distribution values on bone compared with the NR (control) model. Regarding maximum stress distribution at the root apex, the stress value of the RN model was slightly higher compared to the HN model, whereas the RG model displayed a slightly lower stress value in comparison with the HG model. For maximum tooth displacement value, there were no significant differences between the HN and RN models, as well as the HG and RG models.

Conclusion

The round resection pattern had comparable stress distribution behaviors at the root apex and tooth displacement values with the horizontal resection pattern. Targeted endodontic microsurgery might provide better biomechanical response of the resected tooth after root-end resection.

A Safe Technique for Excising the Perpendicular Plate of the Ethmoid Bone in Patients with Crooked Nose: A Finite Element Analysis

BACKGROUND

Correction of the crooked nose, especially the perpendicular plate of the ethmoid bone, has the potential to cause skull base injury. At present, the safe and effective method for perpendicular plate resection has not been clearly defined through biomechanics.

METHOD

CT scan data of 48 patients with crooked nose and deviated nasal septum were divided into C-type, angular deformity-type, and S-type based on the morphology of the 3D model. Different types of finite element models of the nasal bony septum and skull base were established. The osteotomy depth, angle, and force mode of the PPE resection were simulated by assembling different working conditions for the models. The von Mises stress of the anterior cranial fossa was observed.

RESULTS

When the osteotomy line length was 0.5 cm, the angle was at 30° to the Frankfurt plane, and 50 N·mm torque was applied, the von Mises stress of the skull base was minimal in the four models, showing 0.049 MPa (C-type), 0.082 MPa (S-type), 0.128 MPa (angular deformity-type), and 0.021 MPa (control model). The maximum von Mises stress values were found at the skull base when the osteotomy line was 1.5 cm, the angle was 50°, and the force was 10 N along the X-axis, showing 0.349 MPa (C-type), 0.698 MPa (S-type), 0.451 MPa (angular deformity-type), and 0.149 MPa (control model).

CONCLUSION

The use of smaller resection angle with the Frankfurt plane, conservative resection depth, and torsion force can better reduce the stress value at the skull base and reduce the risk of basicranial fracture. It is a safe and effective technique for perpendicular plate resection of the ethmoid bone in the correction of crooked nose.

LEVEL OF EVIDENCE IV

This journal requires that authors assign a level of evidence to each article. For a full description of these Evidence-Based Medicine ratings, please refer to the Table of Contents or the online Instructions to Authors www.springer.com/00266 .

A mechanical and three-dimensional finite element study of the optimum tooth sectioning depth during the extraction of low-level horizontally impacted mandibular third molar

The present study aims to determine the optimum sectioning depth for the extraction of low-level horizontally impacted mandibular third molar (LHIM3M) using mechanical and finite element analysis. One hundred and fifty extracted mandibular third molars were randomly divided into three groups: 1, 2 or 3 mm of tooth tissue was retained at the bottom of the crown. The breaking force of teeth was tested in a universal strength testing machine. The fracture surface was observed and the type of tooth breakage was recorded. According to the three groups, corresponding 3D finite element models were created. The breaking force obtained in the mechanical study was, respectively, applied and the stress and strain of the teeth and surrounding tissues were analysed. Breaking force decreased as sectioning depth increased. The 2 mm group produced the lowest rate of incomplete breakage (10%). In the 2 mm model, the stresses were evenly distributed in the tooth tissue at the bottom of the fissure, and the maximal stress was located in the tissue close to the root segment. The maximum values of stresses in the bone and of strains in the periodontal ligament of the second molar and bone were lower in the 1 mm model than in other models. Their distribution was similar in the three models. A sectioning depth of 1 mm group saves labour during the extraction of LHIM3M, compared to 2 and 3 mm; 2 mm might be the appropriate sectioning depth in terms of breakage shapes.

Clinical Efficacy Analysis of the Personalization of Prosthetic Abutments in Implant Supported Restorations in Comparison to Available Standard Titanium Abutments

Personalized medicine has become an important direction to offer better solutions for health problems. In implantology, this trend was materialized through customizing dental abutments to each clinical situation. The demands for better esthetics and function of implant-supported restorations have imposed a more personalized variety of prosthetic abutments. This retrospective study compared clinical efficiency of personalized implant abutments with standard implant abutments in multiple implant restorations. Clinical data of patients who were admitted in a private clinic between 2011 and 2022 and received dental implant treatments were collected. All complications and undesired events from the patients' medical record charts were statistically analyzed. The implants were loaded using either standard or customized abutments. For complete arch rehabilitations with the SKY Fast & Fixed protocol, standard titanium prosthetic abutments were used. Our results suggest that the abutments choice for patients has moved throughout the years more towards the use of customized abutments. The number of customized abutments (414) was higher compared with the number of standard abutments (293). In our database, the most used abutments for the anterior area implants were made of titanium and zirconia, whereas for the posterior area, the preferred abutments were mostly titanium. The standard abutments were used almost entirely for immediate loading and implantation in both anterior and posterior areas (Fast & Fixed protocol). Complications were encountered mainly in restorations with standard abutments (9.22%) compared to customized abutments (2.7%), with titanium abutments being the most reliable, having only 1.79% complications.

Fracture Resistance of Class II MOD Cavities Restored by Direct and Indirect Techniques and Different Materials Combination

This study aimed to evaluate the fracture resistance of class II MOD cavities restored using different techniques and materials. Sixty extracted maxillary molars were selected and standardized class II MOD cavities were prepared using a custom-made paralleling device. The specimens were divided into four groups based on the restoration technique used: Group 1 (direct resin composite), Group 2 (short-fiber-reinforced composite resin), Group 3 (composite polyethylene fiber reinforcement), and Group 4 (CAD/CAM resin inlays). Fracture resistance was assessed for each group after thermocycling aging for 10,000 cycles. The mode of fracture was assigned to five types using Burke's classification. To compare the fracture force among the tested materials, a paired sample t-test was performed. The significance level for each test was set at p < 0.05. Significant differences in fracture resistance were observed among the different restoration techniques. CAD/CAM inlays (2166 ± 615 N), short-fiber-reinforced composite resin (2471 ± 761 N), and composite polyethylene fiber reinforcement (1923 ± 492 N) showed superior fracture resistance compared to the group restored with direct resin composite (1242 ± 436 N). The conventional resin composite group exhibited the lowest mean fracture resistance. The choice of restoration material plays a critical role in the clinical survival of large MOD cavities. CAD/CAM inlays and fiber-reinforced composites offer improved fracture resistance, which is essential for long-term success in extensive restorations.

Evaluation of biomechanics using different traction devices in distalization of maxillary molar with clear aligners: a finite element study

The study aimed to mechanically evaluate the tooth displacement of molar distalization by clear aligners combined with micro-implant through different traction devices using finite element analysis. A three-dimensional finite element model of complete maxillary dentition was constructed. Simultaneously move the maxillary first and second molars 0.2 mm distally at the height of 4 mm and 6 mm of micro-implant, and 150 g force was applied to button, precision cut and angelbutton respectively. Initial tooth movement in six different conditions of anterior tooth and molars was analyzed and calculated with ANSYS software. All the upper anterior tooth exhibited uncontrolled labial tipping and intrusion upon the six conditions, and the central incisor showed the largest tendency of crown labial inclination. Among the absolute values of crown-root displacement difference of the anterior tooth in sagittal direction, the angelbutton was the smallest, which means the torque control ability was superior to others. However, button played a more accurate role in the sagittal and vertical control of canine. With the increase of micro-implant height, the torque control ability of anterior tooth was decreased, but the intrusion trend increased. The controlled distal inclination with extrusion of the first molar and uncontrolled distal inclination with intrusion of the second molar were observed, and the angelbutton had more effective horizontal and vertical control on molars, which was close to bodily movement than others. As a new type of traction device, angelbutton has excellent anchorage control effect in clear aligners therapy of molar distalization, which further realizes the accurate expression of orthodontic force.

Bonding Performance of Surface-Treated Zirconia Cantilevered Resin-Bonded Fixed Dental Prostheses: In Vitro Evaluation and Finite Element Analysis

Debonding of zirconia cantilevered resin-bonded fixed dental prostheses (RBFDPs) remains the main treatment complication, therefore, the present in vitro study aimed to evaluate the effect of different surface pretreatments on the bonding of zirconia RBFDPs. Eighty milled zirconia maxillary central incisors, with complementary zirconia cantilevered RBFDPs, were randomly subjected to four different surface pretreatments (n = 20): as-machined (AM); airborne-particle abraded (APA); coated with nanostructured alumina coating (NAC); incisor air-abraded and RBFDP coated (NAC_APA). After bonding, half of each group (n = 10) was stored in deionized water (150 days/37 °C), thermocycled (37,500 cycles, 5-55 °C), and cyclically loaded (50 N/1.2 × 10⁶). Load-bearing capacity (LBC) was determined using a quasi-static test. Additionally, finite element analysis (FEA) and fractography were performed. t-test and one-way ANOVA were used for statistical-analysis. Before aging, the NAC group provided superior LBC to other groups (p < 0.05). After aging, the AM specimens debonded spontaneously, while other groups exhibited comparable LBC (p ˃ 0.05). The FEA results correlated with the in vitro experiment and fractography, showing highly stressed areas in the bonding interface, cement layer, and in RBFDP's retainer wing and connector. The NAC RBFDPs exhibited comparable long-term bonding performance to APA and should be regarded as a zirconia pretreatment alternative to APA.

Analysis of initial stress distribution in palatal bone around the implant in lingual orthodontics for single and double palatal implant systems: a FEM study

Objective:

To analyze and compare the Von Mises stress and principal stress distribution in palatal bone around the palatal implant in lingual orthodontics (LiO) for single and double palatal implant systems with varying lengths of lever arm.

Methods:

Two groups were assessed: single (Group 1) and double (Group 2) palatal implant systems, which were further divided into two subgroups, based on lever arm length, for analyzing stress in the palatal bone around the implant. Hence, two 3D finite element models of bilateral maxillary first premolar extraction cases were constructed in each system. Lingual brackets (0.018-in slot) were positioned at the center of the clinical crown. In both systems, 150g of retraction force was applied, and ANSYS v. 12.1 software was used to analyze and compare stress in the palatal bone around the palatal implant.

Results:

In this study, higher stress was observed at the inner threaded interface of cortical bone. Magnitude of Von Mises stress was higher in Group 2 (0.63 MPa and 0.65 MPa) in comparison to Group 1 (0.29 MPa and 0.29 MPa). Similarly, magnitude of principal stress was higher in Group 2, in comparison to Group 1. Higher stress was observed in the apical region of the implant-bone interface of cancellous bone.

Conclusion:

This study concluded that the Von Misses stress as well as principal stress in the palatal bone were within the optimal limit in both groups. Finally, it can be concluded that both systems (single and double palatal implant) were safe for the patients in clinical use of 150g of retraction force.

Precision medicine using patient-specific modelling: state of the art and perspectives in dental practice

The dental practice has largely evolved in the last 50 years following a better understanding of the biomechanical behaviour of teeth and its supporting structures, as well as developments in the fields of imaging and biomaterials. However, many patients still encounter treatment failures; this is related to the complex nature of evaluating the biomechanical aspects of each clinical situation due to the numerous patient-specific parameters, such as occlusion and root anatomy. In parallel, the advent of cone beam computed tomography enabled researchers in the field of odontology as well as clinicians to gather and model patient data with sufficient accuracy using image processing and finite element technologies. These developments gave rise to a new precision medicine concept that proposes to individually assess anatomical and biomechanical characteristics and adapt treatment options accordingly. While this approach is already applied in maxillofacial surgery, its implementation in dentistry is still restricted. However, recent advancements in artificial intelligence make it possible to automate several parts of the laborious modelling task, bringing such user-assisted decision-support tools closer to both clinicians and researchers. Therefore, the present narrative review aimed to present and discuss the current literature investigating patient-specific modelling in dentistry, its state-of-the-art applications, and research perspectives.

Finite element analysis and nanomechanical properties of composite and ceramic dental onlays

The aim is to evaluate the mechanical properties of the composite and ceramic based indirect restorative materials used in dental treatments with scanning nanoindentation test (NT). Finite element analysis (FEA) was applied to investigate the stress distribution. Four hybrid composite materials; Indirect resin composite (IRC), Resin nanoceramic (RNC), Polymer infiltrated ceramic (PIC) and Zirconia-reinforced lithium-di-silicate (ZRC) were divided into two subgroups for NT (n = 20) and fracture test (n = 40). Statistical analyses were performed with independent t-test, ANOVA and post-hoc Tukey tests (p ≤ 0.05). The highest hardness, elasticity and fracture toughness were observed in ZRC (p = 0.001). Frequency of vertical root fractures in RNC and IRC were statistically lower than ZRC (p = 0.032). Reinforced CAD-CAM ceramics revealed higher mechanical properties compared with IRC materials. The FEA model for fracture mechanism of RNC demonstrated lowest stress values and uniform stress distribution amongst all groups, while ZRC and PIC presented the highest fracture toughness.

Experimental validation of a micro-CT finite element model of a human cadaveric mandible rehabilitated with short-implant-supported partial dentures

PURPOSE

This study aimed to address the predictive value of a micro-computed tomography (μCT)-based finite element (μFE) model of a human cadaveric edentulous posterior mandible, rehabilitated by short dental implants. Hereby, three different prosthetic/implant configurations of fixed partial dentures ("Sp"-3 splinted crowns on 3 implants, "Br" - Bridge: 3 splinted crowns on 2 implants, and "Si"- 3 single crowns) were analysed by comparing the computational predictions of the global stiffness with experimental data.

METHODS

Experimental displacement of the bone/implant/prosthesis system was measured under axial and oblique loads of 100 N using an optical deformation system (GOM Aramis) and the overall movement of the testing machine (Zwick Z030). Together with the measured machine force, an "Aramis" (optical markers) and "Zwick" (test machine) stiffness were calculated. FE models were created based on μCT-scans of the cadaveric mandible sample (n = 1) before and after implantation and using stl-files of the crowns. The same load tests and boundary conditions were simulated on the models and the μFE-results were compared to experimental data using linear regression analysis.

RESULTS

The regression line through a plot of pooled stiffness values (N/mm) for the optical displacement recording (true local displacement) and the test machine (machine compliance included) had a slope of 0.57 and a correlation coefficient R² of 0.82. The average pooled correlation of global stiffness between the experiment and FE-analysis (FEA) showed a R² of 0.80, but the FEA-stiffness was 7.2 times higher. The factor was highly dependent on the test configuration. Sp-configuration showed the largest stiffness followed by Br-configuration (17% difference in experiment and 21% in FEA).

CONCLUSIONS

The current study showed good qualitative agreement between the experimental and predicted global stiffness of different short implant configurations. It could be deduced that 1:1 splinting of the short implants by the crowns is most favorable for the stiffness of the implant/prosthesis system. However, in the clinical context, the absolute in silico readings must be interpreted cautiously, as the FEA showed a considerable overestimation of the values.

Polymer-Based Bioactive Luting Agents for Cementation of All-Ceramic Crowns: An SEM, EDX, Microleakage, Fracture Strength, and Color Stability Study

The aim of the study was to compare microleakage and fracture loads of all ceramic crowns luted with conventional polymer resins and polymeric bioactive cements and to assess the color stability of polymeric bioactive cements. Seventy-five extracted premolar teeth were tested for fracture loads and microleakage in all-ceramic crowns cemented with two types of polymeric bioactive cements and resin cements. In addition, the degree of color change for each cement with coffee was assessed. Thirty maxillary premolar teeth for fracture loads and thirty mandibular premolar teeth for microleakage were prepared; standardized teeth preparations were performed by a single experienced operator. All prepared specimens were randomly distributed to three groups (n = 20) based on the type of cement, Group 1: resin cement (Multilink N); Group 2: polymeric bioactive cement (ACTIVA); Group 3: polymeric bioactive cement (Ceramir). The cementation procedures for all cements (Multilink, ACTIVA, and Ceramir) were performed according to the manufacturers' instructions. All specimens were aged using thermocycling for 30,000 cycles (5-55 °C, dwell time 30 s). These specimens were tested using the universal testing machine for fracture strength and with a micro-CT for microleakage. For the color stability evaluation, the cement specimens were immersed in coffee and evaluated with a spectrometer. Results: The highest and lowest means for fracture loads were observed in resin cements (49.5 ± 8.85) and Ceramir (39.8 ± 9.16), respectively. Ceramir (2.563 ± 0.71) showed the highest microleakage compared to resin (0.70 ± 0.75) and ACTIVA (0.61 ± 0.56). ACTIVA cements showed comparable fracture loads, microleakage, and stain resistance compared to resin cements.

Mechanical failure of posterior teeth due to caries and occlusal wear- A modelling study

OBJECTIVES

The purpose of the present work is to explore the effect of occlusal wear and different types and degrees of caries on the mechanical performance and structural integrity of posterior teeth.

METHODS

Three-dimensional (3D) computational models with different combinations of caries parameters (caries location, caries size and caries induced pulp shrinkage) and occlusal wear factors (enamel thickness, marginal ridge height and cuspal slope) were developed and analyzed using the extended finite element method (XFEM) to identify the stress distribution, crack initiation load and ultimate fracture load values. The effect of a non-drilling conservative treatment using resin infiltration on the recovery of mechanical properties of carious molar teeth was also investigated.

RESULTS

Presence of fissural caries, worn proximal marginal ridge and decreased enamel thickness due to occlusal wear, imparted a significant negative effect on the crack initiation load value of the lower molar models. Accordingly, models with intact and strong proximal marginal ridge, generally exhibited higher crack initiation loading, regardless of caries size and location. Presence of fissure caries drastically decreased (55%-70%) the crack initiation load compared to sound teeth. The depth of the fissural lesion and the presence of proximal caries did not have a major effect on crack initiation load values. However, increasing the caries size resulted in lower final fracture load values in most of the cases. Accordingly, the groups with combined and connected large fissural and proximal lesions experienced the largest drop in the fracture load values compared to sound tooth models. The worst condition consisted of two connected large proximal and fissural caries with no proximal marginal ridge, in which the fracture load dramatically decreased to only 25% of that for sound teeth with intact marginal ridge. On the other hand, decreased cuspal slope due to occlusal wear and shrinkage of the pulp due to caries appeared to have a protective role and a direct relation with the fracture resistance of the tooth. Following the application of resin infiltration on the carious models, the crack initiation load and the fracture load could recover up to 75% and 90% of the values for the corresponding sound tooth models, respectively.

SIGNIFICANCE

Presence of fissural caries, if not treated (either with remineralization, resin infiltration or restoration), can be a major risk factor in the initiation of tooth fracture. When combined with decreased enamel thickness and loss of proximal marginal ridge due to mechanical or chemical wear, the weakening effect of the caries will be amplified specially in teeth with steep cuspal slopes. The application of a conservative treatment with resin infiltration can be an effective approach in prevention of further mechanical failure of demineralized enamel. The findings of this study emphasize the importance of early interventions in the management of caries for the prevention of future cuspal or tooth fracture especially in subjects with higher risk factors for tooth fracture such as caries, wear and bruxism.

Novel and simplified optimisation pathway using response surface and design of experiments methodologies for dental implants based on the stress of the cortical bone

Dental implants are widely used as a long-term treatment solution for missing teeth. A titanium implant is inserted into the jawbone, acting as a replacement for the lost tooth root and can then support a denture, crown or bridge. This allows discreet and high-quality aesthetic and functional improvement, boosting patient confidence. The use of implants also restores normal functions such as speech and mastication. Once an implant is placed, the surrounding bone will fuse to the titanium in a process known as osseointegration. The success of osseointegration is dependent on stress distribution within the surrounding bone and thus implant geometry plays an important role in it. Optimisation analyses are used to identify the geometry which results in the most favourable stress distribution, but the traditional methodology is inefficient, requiring analysis of numerous models and parameter combinations to identify the optimal solution. A proposed improvement to the traditional methodology includes the use of Design of Experiments (DOE) together with Response Surface Methodology (RSM). This would allow for a well-reasoned combination of parameters to be proposed. This study aims to use DOE, RSM and finite element models to develop a simplified optimisation analysis method for dental implant design. Drawing on data and results from previous studies, two-dimensional finite element models of a single Branemark implant, a multi-unit abutment, two prosthetic screws, a prosthetic crown and a region of mandibular bone were built. A small number of combinations of implant diameter and length were set based on the DOE method to analyse the influence of geometry on stress distribution at the bone-implant interface. The results agreed with previous studies and indicated that implant length is the critical parameter in reducing stress on cortical bone. The proposed method represents a more efficient analysis of multiple geometrical combinations with reduced time and computational cost, using fewer than a third of the models required by the traditional methods. Further work should include the application of this methodology to optimisation analyses using three-dimensional finite element models.

[Study of the effect of the occlusal surface of a natural tooth and fixed partial dentures supported by dental implants on the stress distribution by finite element analysis]

OBJECTIVE

To study the effect of the occlusal surface of a natural tooth, cement-retained and screw-retained fixed partial dentures supported by dental implants and the coefficient of friction on the stress distribution in the peri-implant bone at maximum and minimum principal stresses.

MATERIAL AND METHODS

Study of maximum and minimum principal stresses in models with natural teeth and artificial crowns supported by dental implants, taking into account the coefficient of friction using the finite element analysis.

RESULTS

In models represented by fixed partial dentures supported by dental implants and a natural tooth, the maximum tensile stresses arise in the cortical bone in the cervical region of the artificial crown, and the maximum compressive stresses occur both in the cortical layer in the cervical region of the artificial crown and in the cervical region of the tooth. In models with two fixed partial dentures supported by dental implants or two natural teeth, the stress distributions in the cortical layers in the upper and lower jaw are almost identical.

CONCLUSION

Modeling the antagonist and adding to the FEA model is important in order to determine the precise and realistic direction of the resulting force vector. Amplification of the number of contact areas should be considered when modeling the occlusal surface of artificial crowns supported by dental implants.

OpenMandible: An open-source framework for highly realistic numerical modelling of lower mandible physiology

OBJECTIVE

Computer modeling of lower mandible physiology remains challenging because prescribing realistic material characteristics and boundary conditions from medical scans requires advanced equipment and skill sets. The objective of this study is to provide a framework that could reduce simplifications made and inconsistency (in terms of geometry, materials, and boundary conditions) among further studies on the topic.

METHODS

The OpenMandible framework offers: 1) the first publicly available multiscale model of the mandible developed by combining cone beam computerized tomography (CBCT) and μCT imaging modalities, and 2) a C++ software tool for the generation of simulation-ready models (tet4 and hex8 elements). In addition to the application of conventional (Neumann and Dirichlet) boundary conditions, OpenMandible introduces a novel geodesic wave propagation - based approach for incorporating orthotropic micromechanical characteristics of cortical bone, and a unique algorithm for modeling muscles as uniformly directed vectors. The base intact model includes the mandible (spongy and compact bone), 14 teeth (comprising dentin, enamel, periodontal ligament, and pulp), simplified temporomandibular joints, and masticatory muscles (masseter, temporalis, medial, and lateral pterygoid).

RESULTS

The complete source code, executables, showcases, and sample data are freely available on the public repository: https://github.com/ArsoVukicevic/OpenMandible. It has been demonstrated that by slightly editing the baseline model, one can study different "virtual" treatments or diseases, including tooth restoration, placement of implants, mandible bone degradation, and others.

SIGNIFICANCE

OpenMandible eases the community to undertake a broad range of studies on the topic, while increasing their consistency and reproducibility. At the same time, the needs for dedicated equipment and skills for developing realistic simulation models are significantly reduced.

The effect of the design of a mandibular implant-supported zirconia prosthesis on stress distribution

STATEMENT OF PROBLEM

Prosthetic complications have been frequently reported in implant-supported complete-arch prosthesis. Prosthetic restorations designed with an all-on-four treatment concept and fabricated from zirconia ceramic may be used to overcome these problems.

PURPOSE

The purpose of this biomechanical study was to evaluate the effects of cantilever length and inclination of implant on the stress distribution in bone tissue, implant, and a monolithic zirconia ceramic-lithium disilicate glass-ceramic superstructure for all-on-four prosthesis.

MATERIAL AND METHODS

All-on-four mandibular prosthesis fabricated from a zirconia and lithium disilicate glass-ceramic (LDGC) superstructure was designed with cantilever lengths of either 5 mm or 9 mm and posterior implants with a distal tilt of either 15 or 30 degrees. Stresses were evaluated with a simulated application of a static load of 600 N.

RESULTS

Increasing implant inclination from 15 to 30 degrees led to a decrease in maximum principal stress (MaxPS) values of approximately 4 to 7 MPa in cortical bone around all implants except the right anterior implant in the designs with short cantilevers and an increase in MaxPS values (approximately 3 to 19 MPa) in the same places in the designs with the long cantilevers. Increasing cantilever length from 5 to 9 mm resulted in an increase in minimum principal stress (MinPS) values of approximately 3 to 13 MPa in the cortical bone surrounding all posterior implants. In the designs with the long cantilever, MaxPS values increased approximately 3 to 4 MPa in spongy bone adjacent to the right posterior implant. An increase in cantilever length also led to higher vMS values at the first and second implant grooves in the right posterior implant in the design with the 15-degree implant tilt. An increase in implant inclination in the design with the short cantilever resulted in lower vMS values at the apex and all grooves of the left posterior implant, whereas in the design with the long cantilever, an increase in implant inclination resulted in lower stress values in the first and second grooves of the same implant. An increase in implant inclination led to in an increase in vMS values in the core structure.

CONCLUSIONS

In zirconia ceramic restorations by using an all-on-four design with an LDGC superstructure, short cantilevers may be preferable because they result in a more favorable distribution of stress than long cantilevers. An increase in implant angulation from 15 to 30 degrees decreased MaxPS values in cortical bone.

Validated Finite Element Models of Premolars: A Scoping Review

Finite element (FE) models are widely used to investigate the biomechanics of reconstructed premolars. However, parameter identification is a complex step because experimental validation cannot always be conducted. The aim of this study was to collect the experimentally validated FE models of premolars, extract their parameters, and discuss trends. A systematic review was performed following Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. Records were identified in three electronic databases (MEDLINE [PubMed], Scopus, The Cochrane Library) by two independent reviewers. Twenty-seven parameters dealing with failure criteria, model construction, material laws, boundary conditions, and model validation were extracted from the included articles. From 1306 records, 214 were selected for eligibility and entirely read. Among them, 19 studies were included. A heterogeneity was observed for several parameters associated with failure criteria and model construction. Elasticity, linearity, and isotropy were more often chosen for dental and periodontal tissues with a Young’s modulus mostly set at 18–18.6 GPa for dentine. Loading was mainly simulated by an axial force, and FE models were mostly validated by in vitro tests evaluating tooth strains, but different conditions about experiment type, sample size, and tooth status (intact or restored) were reported. In conclusion, material laws identified herein could be applied to future premolar FE models. However, further investigations such as sensitivity analysis are required for several parameters to clarify their indication.

Biomechanical Evaluation of Initial Stability of a Root Analogue Implant Design with Drilling Protocol: A 3D Finite Element Analysis

Background: The aim of this study was to biomechanically evaluate the initial stability of a patient-specific root analogue implant (RAI) design with drilling protocol by comparing it to designs without drilling protocol through a 3D finite element analysis (FEA). Methods: A 3D surface model of an RAI for the upper right incisor was constructed. To evaluate the effect of root apex drilling, four modified RAI shapes were designed with the press-fit implantation method: Non-modified, wedge added at root surface, lattice added at root surface, and apex-anchor added at root apex (AA). Each model was subjected to an oblique load of 100 N. To simulate the initial stability of implantation, contact conditions at the implant–bone interface were set to allow for the sliding phenomenon with low friction (frictional coefficient 0.1–0.5). Analysis was performed to evaluate micro-displacements of the implants and peak stress on the surrounding bones. Results: Under all low frictional coefficient conditions, the lowest von Mises stress level on the cortical bone and fewest micro-displacements of the implant were observed in the AA design. Conclusion: In view of these results, the AA design proved superior in reducing the stress concentration on the supporting cortical bone and the micro-displacement of RAI.

Root canal reconstruction using biological dentin posts: A 3D finite element analysis

Background. Several types of post have been developed for clinical use. A biological dentin post obtained from an extracted tooth eliminates the problems arising from material differences and reduces the fracture rate in teeth undergoing root canal treatment. This study used finite element analysis to compare a biological dentin post with posts made of two different materials. Methods. Three 3D models of the upper central incisor were created, and stainless-steel, glass fiber and biological dentin posts were applied to these models. The restoration of the models was completed by applying a composite as the core structure and a ceramic crown as the superstructure. Using finite element stress analysis in the restoration models, a 100-N force was applied in the vertical and horizontal directions and at a 45º angle, and the suitability of the biological dentin post was evaluated by comparing the data. Results. Under the applied forces, the greatest stress accumulation was seen in the models with the stainless steel post. Because the stainless steel post was more rigid, stress forces accumulated on the surface instead of being transmitted to the tooth tissue. In the models with the glass fiber and biological dentin posts, the post material responded to the stratification in tandem with the dental tissue and did not cause excessive stress accumulation on the tooth or post surfaces. Conclusion. The results showed that biological dentin posts prevent the accumulation of stresses that might cause fractures in teeth undergoing root canal treatment. In addition, the physical compatibility and biocompatibility of a biological dentin post with the tooth imply that it is a good alternative to the types of post currently used.

Influence of block-out on retentive force of thermoplastic resin clasps: an in vitro experimental and finite element analysis

PURPOSE

The purpose of this study was to determine the effect of block-out preparation, used to eliminate the undercut area, on the retentive force and stress distribution of resin clasps.

METHODS

A total of 72 polyester and polyamide resin clasps were fabricated on a premolar abutment crown following six block-out preparations. A combination of two types of vertical block-outs and three types of horizontal block-outs (on the missing side) was used on the abutment tooth. Each clasp was subjected to an in vitro removal test using a universal testing machine. The retentive force and traces of the clasp on the abutment tooth were recorded and analyzed with one-way analysis of variance and post hoc comparisons (α=0.05). Non-linear finite element analysis was performed to assess the stress distributions of the resin clasps.

RESULTS

Resin clasps with a vertical block-out of 0.75mm undercut showed significantly higher retentive force than those with the 0.5mm undercut. Resin clasps with horizontal block-out showed significantly lower retentive force than those without horizontal block-out. There was no significant difference between the two thermoplastic resins. The maximum first principal stress of the resin clasp was concentrated under the shoulder of the clasp and strongly influenced by the width of horizontal block-out in the finite element analysis.

CONCLUSIONS

This in vitro experiment suggested that a horizontal block-out is necessary even for a 0.5-mm undercut. The design of the resin clasp should be considered from two aspects: retentive force and deformation risk.

Biomechanical characteristics of immediately loaded and osseointegration dental implants inserted into Sika deer antler

This study aimed to compare biomechanical characteristics of immediately loaded (IL) and osseointegrated (OS) dental implants inserted into Sika deer antler and lay a foundation for developing an alternative animal model for dental implants studies. Two implants per antler were inserted. One implant was loaded immediately via a self-developed loading device; the other was submerged and unloaded as control. IL implants were harvested after different loading periods. The unloaded implants were collected after OS and the shedding of antler. Specimens were scanned by µCT scanner and finite element models were generated. A vertical force of 10 N was applied on the implant. The mean values of maximum displacements, stresses and strains were compared. The results showed that the density of antler tissue around the implants dramatically increased as the loading time increased. After shedding the antler, 3 pairs of antlers were collected and the density of antler tissue remained in a similar value in all specimens. The maximum values of displacement and stresses in implant and stresses and strains in antler tissue were significantly different among OS models. In one antler, all the biomechanical parameters of IL model were significantly higher than those of OS model of the same animal (P < 0.05) and wider distributions were obtained from IL model. It can be concluded that implants inserted into Sika deer antler might not disturb the growth and calcification process of antler and the use of Sika deer antler model is a promising alternative for implant studies that does not require animal sacrifice.

Numerical fatigue analysis of premolars restored by CAD/CAM ceramic crowns

OBJECTIVES

The purpose of this study was to estimate the fatigue life of premolars restored with two dental ceramics, lithium disilicate (LD) and polymer infiltrated ceramic (PIC) using the numerical method and compare it with the published in vitro data.

METHODS

A premolar restored with full-coverage crown was digitized. The volumetric shape of tooth tissues and crowns were created in Mimics^®. They were transferred to IA-FEMesh for mesh generation and the model was analyzed with Abaqus. By combining the stress distribution results with fatigue stress-life (S-N) approach, the lifetime of restored premolars was predicted.

RESULTS

The predicted lifetime was 1,231,318 cycles for LD with fatigue load of 1400N, while the one for PIC was 475,063 cycles with the load of 870N. The peak value of maximum principal stress occurred at the contact area (LD: 172MPa and PIC: 96MPa) and central fossa (LD: 100MPa and PIC: 64MPa) for both ceramics which were the most seen failure areas in the experiment. In the adhesive layer, the maximum shear stress was observed at the shoulder area (LD: 53.6MPa and PIC: 29MPa).

SIGNIFICANCE

The fatigue life and failure modes of all-ceramic crown determined by the numerical method seem to correlate well with the previous experimental study.

A REVIEW OF FINITE ELEMENT APPLICATIONS IN ORAL AND MAXILLOFACIAL BIOMECHANICS

Finite element method (FEM) is preferred to carry out mechanical analyses for many complex biomechanical structures. For most of the biomechanical models such as oral and maxillofacial structures or patient-specific dental instruments, including nonlinearities, complicated geometries, complex material properties, or loading/boundary conditions, it is not possible to accomplish an analytical solution. The FEM is the most widely used numerical approach for such cases and found a wide range of application fields for investigating the biomechanical characteristics of oral and maxillofacial structures that are exposed to external forces or torques. The numerical results such as stress or strain distributions obtained from finite element analysis (FEA) enable dental researchers to evaluate the bone tissues subjected to the implant or prosthesis fixation from the viewpoint of (i) mechanical strength, (ii) material properties, (iii) geometry and dimensions, (iv) structural properties, (v) loading or boundary conditions, and (vi) quantity of implants or prostheses. This review paper evaluates the process of the FEA of the oral and maxillofacial structures step by step as followings: (i) a general perspective on the techniques for creating oral and maxillofacial models, (ii) definitions of material properties assigned to oral and maxillofacial tissues and related dental materials, (iii) definitions of contact types between tissue and dental instruments, (iv) details on loading and boundary conditions, and (v) meshing process.

A three-dimensional topology optimization model for tooth-root morphology

To obtain the root of a lower incisor through structural optimization, we used two methods: optimization with Solid Isotropic Material with Penalization (SIMP) and Soft-Kill Option (SKO). The optimization was carried out in combination with a finite element analysis in Abaqus/Standard. The model geometry was based on cone-beam tomography scans of 10 adult males with healthy bone-tooth interface. Our results demonstrate that the optimization method using SIMP for minimum compliance could not adequately predict the actual root shape. The SKO method, however, provided optimization results that were comparable to the natural root form and is therefore suitable to set up the basic topology of a dental root.

New Design for an Adjustable Cise Space Maintainer

Objective

The aim of this study is to present a new adjustable Cise space maintainer for preventive orthodontic applications.

Methods

Stainless steel based new design consists of six main components. In order to understand the major displacement and stress fields, structural analysis for the design is considered by using finite element method.

Results

Similar to major displacement at y-axis, critical stresses σ_x and τ_xy possess a linear distribution with constant increasing. Additionally, strain energy density (SED) plays an important role to determine critical biting load capacity.

Conclusion

Structural analysis shows that the space maintainer is stable and is used for maintaining and/or regaining the space which arouses early loss of molar tooth.

The ball-on-disk cyclic wear of CAD/CAM machinable dental composite and ceramic materials

The purpose of this study was to investigate the wear volume and the principal strain of machinable dental composite and ceramics in simulated mastication. A ball-on-disk wear test was performed for 3,000 cycles in water, using nine ball/disk combinations of three commercial CAD/CAM materials: feldspathic, lithium disilicate glass ceramics, and a highly loaded composite material (n = 7 for each combination). The wear volume was optically measured using a digital scanner and analyzed for statistical differences based on the materials (α = 0.05). We used non-linear finite element analysis to calculate the principal strain. The wear volume of the ball was significantly larger than that of the disk when hardness and fracture toughness of the former was lower than that of the latter and vice versa (P < 0.05). The lithium disilicate glass ceramic constantly showed lower wear volume than the opposing antagonist. Except for the same material pairs of feldspathic and composite, the ball or disk specimen that showed a larger wear in the occluding pair coincided with the one with higher maximum strain. It was not possible to predict the magnitude of wear, whereas the result suggested a strong association between the maximum strain and wear volume of the ceramic surface.

Occlusal loading during biting from an experimental and simulation point of view

OBJECTIVES

Occlusal loading during clenching and biting is achieved by the action of the masticatory system, and forms the basis for the evaluation of the functional performance of prosthodontic and maxillofacial components. This review provides an overview of (i) current bite force measurement techniques and their limitations and (ii) the use of computational modelling to predict bite force. A brief simulation study highlighting the challenges of current computational dental models is also presented.

METHODS

Appropriate studies were used to highlight the development and current bite force measurement methodologies and state-of-the-art simulation for computing bite forces using biomechanical models.

RESULTS

While a number of strategies have been developed to measure occlusal forces in three-dimensions, the use of strain-gauges, piezo-electric sensors and pressure sheets remain the most widespread. In addition to experimental-based measurement techniques, bite force may be also estimated using computational models of the masticatory system. Simulations of different bite force models clearly show that the use of three-dimensional force measurements enriches the evaluation of masticatory functional performance.

SIGNIFICANCE

Hence, combining computational modelling with three-dimensional force measurement techniques can significantly improve the evaluation of masticatory system and the functional performance of prosthodontic components.

Effects of implant tilting and the loading direction on the displacement and micromotion of immediately loaded implants: an in vitro experiment and finite element analysis

Purpose

The purpose of this study was to investigate the effects of implant tilting and the loading direction on the displacement and micromotion (relative displacement between the implant and bone) of immediately loaded implants by in vitro experiments and finite element analysis (FEA).

Methods

Six artificial bone blocks were prepared. Six screw-type implants with a length of 10 mm and diameter of 4.3 mm were placed, with 3 positioned axially and 3 tilted. The tilted implants were 30° distally inclined to the axial implants. Vertical and mesiodistal oblique (45° angle) loads of 200 N were applied to the top of the abutment, and the abutment displacement was recorded. Nonlinear finite element models simulating the in vitro experiment were constructed, and the abutment displacement and micromotion were calculated. The data on the abutment displacement from in vitro experiments and FEA were compared, and the validity of the finite element model was evaluated.

Results

The abutment displacement was greater under oblique loading than under axial loading and greater for the tilted implants than for the axial implants. The in vitro and FEA results showed satisfactory consistency. The maximum micromotion was 2.8- to 4.1-fold higher under oblique loading than under vertical loading. The maximum micromotion values in the axial and tilted implants were very close under vertical loading. However, in the tilted implant model, the maximum micromotion was 38.7% less than in the axial implant model under oblique loading. The relationship between abutment displacement and micromotion varied according to the loading direction (vertical or oblique) as well as the implant insertion angle (axial or tilted).

Conclusions

Tilted implants may have a lower maximum extent of micromotion than axial implants under mesiodistal oblique loading. The maximum micromotion values were strongly influenced by the loading direction. The maximum micromotion values did not reflect the abutment displacement values.

The effect of direct and indirect force transmission on peri-implant bone stress - a contact finite element analysis

In almost all finite element (FE) studies in dentistry, virtual forces are applied directly to dentures. The purpose of this study was to develop a FE model with non-linear contact simulation using an antagonist as force transmitter and to compare this with a similar model that uses direct force transmission. Furthermore, five contact situations were created in order to examine their influence on the peri-implant bone stresses, which are relevant to the survival rate of implants. It was found that the peri-implant bone stresses were strongly influenced by the kind of force transmission and contact number.

Three-dimensional finite element analysis of implant-assisted removable partial dentures

STATEMENT OF PROBLEM

Whether the implant abutment in implant-assisted removable partial dentures (IARPDs) functions as a natural removable partial denture (RPD) tooth abutment is unknown.

PURPOSE

The purpose of this 3-dimensional finite element study was to analyze the biomechanical behavior of implant crown, bone, RPD, and IARPD.

MATERIAL AND METHODS

Finite element models of the partial maxilla, teeth, and prostheses were generated on the basis of a patient's computed tomographic data. The teeth, surveyed crowns, and RPDs were created in the model. With the generated components, four 3-dimensional finite element models of the partial maxilla were constructed: tooth-supported RPD (TB), implant-supported RPD (IB), tooth-tissue-supported RPD (TT), and implant-tissue-supported RPD (IT) models. Oblique loading of 300 N was applied on the crowns and denture teeth. The von Mises stress and displacement of the denture abutment tooth and implant system were identified.

RESULTS

The highest von Mises stress values of both IARPDs occurred on the implants, while those of both natural tooth RPDs occurred on the frameworks of the RPDs. The highest von Mises stress of model IT was about twice that of model IB, while the value of model TT was similar to that of model TB. The maximum displacement was greater in models TB and TT than in models IB and IT. Among the 4 models, the highest maximum displacement value was observed in the model TT and the lowest value was in the model IB.

CONCLUSIONS

Finite element analysis revealed that the stress distribution pattern of the IARPDs was different from that of the natural tooth RPDs and the stress distribution of implant-supported RPD was different from that of implant-tissue-supported RPD. When implants are used for RPD abutments, more consideration concerning the RPD design and the number or location of the implant is necessary.

Reliability of FEA on the Results of Mechanical Properties of Materials

The present study evaluated the reliability of FEA on the results of different mechanical properties (E and v) of materials. Two 3D models of a maxillary canine with endodontic treatment, intracanal post, composite resin core and restored with porcelain-fused-to-metal crown were generated according to micro-CT images. Two groups with different E and ν values for porcelain, metal coping alloy, resin cement and composite resin were established. The materials' properties for group GL were based on literature data, while for group GIE the impulse excitation technique was used. A load of 180 N was applied at 45° on the incisal third of the lingual surface of the canine tooth. All models were supported by the periodontal ligament (x=y=z=0). The von Mises stress (VMS) was calculated. The stress values revealed differences between the groups for both VMS distribution and value. The porcelain (GL: 5.966 MPa; GIE: 7.478 MPa), metal coping (GL: 3.811 MPa; GIE: 0.973 MPa) and core (GL: 4.771 MPa; GIE: 0.026 MPa) were significantly affected. In conclusion, this study showed that the determination of mechanical properties (E and ν) of materials is essential for the reliability on the results of FEA.