Zum Hauptinhalt springen

Finite element analysis of a one-piece zirconia implant in anterior single tooth implant applications.

Talmazov, Georgi ; Veilleux, Nathan ; et al.
In: PLoS ONE, Jg. 15 (2020), Heft 2, S. e0229360
Online academicJournal

Finite element analysis of a one-piece zirconia implant in anterior single tooth implant applications  Introduction

This study evaluated the von Mises stress (MPa) and equivalent strain occurring around monolithic yttria-zirconia (Zir) implant using three clinically simulated finite element analysis (FEA) models for a missing maxillary central incisor. Two unidentified patients' cone-beam computed tomography (CBCT) datasets with and without right maxillary central incisor were used to create the FEA models. Three different FEA models were made with bone structures that represent a healed socket (HS), reduced bone width edentulous site (RB), and immediate extraction socket with graft (EG). A one-piece abutment-implant fixture mimicking Straumann Standard Plus tissue level RN 4.1 X 11.8mm, for titanium alloy (Ti) and Zir were modeled. 178 N oblique load and 200 N vertical load were used to simulate occlusal loading. Von Mises stress and equivalent strain values for around each implant model were measured. Within the HS and RB models the labial-cervical region in the cortical bone exhibited highest stress, with Zir having statistically significant lower stress-strain means than Ti in both labial and palatal aspects. For the EG model the labial-cervical area had no statistically significant difference between Ti and Zir; however, Zir performed better than Ti against the graft. FEA models suggest that Ti, a more elastic material than Zir, contributes to the transduction of more overall forces to the socket compared to Zir. Thus, compared to Ti implants, Zir implants may be less prone to peri-implant bone overloading and subsequent bone loss in high stress areas especially in the labial-cervical region of the cortical bone. Zir implants respond to occlusal loading differently than Ti implants. Zir implants may be more favorable in non-grafted edentulous or immediate extraction with grafting.

Recent demands in esthetics have driven research and clinical applications of monolithic zirconia (Zir) and yttria-stabilized tetragonal zirconia polycrystal ceramics or Y-TZP not only as restorative materials, but to also fabricate dental implants. Zir has proven to have similar biological properties and similar short-term survival rate compared to the conventional titanium alloy, Ti-6Al-4V, implants (Ti).[[1]–[3]] Beside esthetics, Zir as a restorative material has proven to be superior to other ceramic materials[[3]–[8]] in terms of mechanical and biological properties including high fracture toughness, high elastic modulus, low thermal conductivity, as well as low dental plaque affinity. One-piece Zir dental implants have been shown to have clinically acceptable mechanical properties even for a reduced diameter implant design.[[7]] Placement of Zir dental implants has been advocated for anterior esthetic zone in a conventional healed edentulous site,[[1]] an immediate extraction site,[[9]] and immediate provisionalization.[[9]] Little information currently exists in the literature on the biomechanical relationship between a one-piece Zir dental implant and its peri-implant osseous structure especially in esthetic zone situations in contrast to the availability of clinical data for Ti. There seemed to be a gap in the literature on Zir dental implants in various clinical situations.

Finite element analysis (FEA) has been applied to understand the bone-implant interface, allowing engineers and clinicians to evaluate implant materials, designs, and their effects toward the surrounding osseous structures.[[10]–[14]] Implant designs have been developed based on FEA modeling that aids our understanding of the bone-implant interface mechanics especially crestal bone stress distribution where bone loading forces are the highest.[[15]] FEA has long been used to understand the functional loading effects on a single tooth implant[[17]] in terms of different implant designs, implant-abutment connections, and restorative designs.[[18]] Simulation of various bone types as well as clinical scenarios such as immediate extraction sockets or grafted sites can also be generated and analyzed using FEA models.[[20]–[21]] Older FEA studies draw their conclusion from models that have been created using simple geometry (ie. cylinders, rectangular blocks) to simplify computations of the mechanical properties of peri-implant bone and haven't really studied how implants behave in periodontal conditions. Recent studies utilize the advantages of cone beam computed tomography (CBCT) to create high resolution STL meshes used to model the heterogenous peri-implant osseous microstructure and environment improving the accuracy of the model.[[21]] An FEA study comparing a Ti-Ti two-piece abutment/implant, a Zir abutment/Ti implant, and a one-piece Zir implant, demonstrated that the one-piece Zir implant generated lower stresses in the peri-implant bone region. Marcián et al. concluded that one-piece Zir implants appeared to have favorably less cervical bone stress, implying peri-implant bone preservation in type III bone when compared to Ti control.[[22]] With esthetic consideration, Zir implant has been indicated in different types of clinical scenarios including healed sockets, immediate implant placement, as well as immediate loading protocol.[[1], [3], [6], [9]] Note also that anterior edentulous site often has reduced bone width resulting in more challenge in implant placement and in long term maintaining peri-implant bone.[[23]]

This study applied FEA to examine the effects of different clinical scenarios that had not been studied biomechanically. This includes a healed socket (HS), reduced bone width edentulous site (RB), and extraction socket with bone grafting (EG) on the peri-implant osseous structure around a one-piece Zir implant. This study utilizes high resolution segmentation from CBCT to create its simulated osseous structures for the various scenarios studied. The rationale of the study was that the different mechanical properties of Ti and Zir may influence the distribution of von Mises stress and equivalent strain on heterogenous peri-implant osseous structures. The null hypothesis then is that Zir implants loaded with different forces in certain situations would transmit those forces into the surrounding peri-implant bone with similar magnitude and distribution as the Ti implants. By accepting or rejecting the null hypothesis that there is no difference between Ti and Zir on peri-implant bone and graft, the goal of the study is to aid the clinician in establishing a rationale for use of Zir implants in various clinical scenarios a patient may present with in the anterior maxilla.

Materials and methods

The first step of the workflow (Fig 1) to create the in silico three dimensional FEA models was selecting CBCT scans from Virginia Commonwealth University (VCU) School of Dentistry internal radiology database with patient identifiers removed prior to the utilization of the scans. Two CBCT datasets were selected based on the appropriate mesio-distal and facio-lingual dimensions available for an esthetic zone single tooth implant, an edentulous maxillary right central incisor and a maxillary right central incisor that was deemed to be extracted and replaced with a single tooth implant.[[25]] The CBCT scan protocol used was iCAT FLX V10 (Imaging Sciences International LLC, Hatfield, PA) with standard implant scan parameters (16 deep 10 high cm volume, 0.3 mm voxel size, 8.9-second scan time, 3.7-second exposure time, 120 kVP, 5 mA, and 501.3 mGy/cm2).[[26]] The CBCT scans were imported into 3D Slicer (https://www.slicer.org/) to generate the initial 3D model[[27]] surfaces of cortical and trabecular bone through the use of threshold operations and the Grow-Cut algorithm. The maxillary anterior region of interest, area approximate to the maxillary right central incisor was isolated.[[25], [28]] The final mesh model was exported as a Standard Tessellation Language (STL) file for final optimization before assembled in the FEA.

Graph: Fig 1 FEA modeling workflow.

The osseous model was subsequently re-meshed and optimized for FEA using MeshLab(http://www.meshlab.net/). While Meshmixer 3.0 software (Autodesk Meshmixer) was used to remove invalid geometry using the Analysis Tool.[[29]] Three FEA osseous models were generated this way representing the healed socket with ideal osseous dimension for a single tooth implant (HS), a reduced bone width model (RB), and an extraction site with grafting (EG). These FEA models were built mimicking the clinical scenarios commonly found for a single tooth replacement for a maxillary first central incisor.[[30]] The HS and RB models were created from the partially edentulous CBCT dataset, while the EG model was created from the dentate CBCT dataset.

The mesh was optimized using MeshLab in order to improve accuracy and speed of the computational analysis.[[32]] Within MeshLab, each model mesh was simplified by using quadric edge collapse decimation function to reduce the total number of initial mesh nodes to 543,804 for the HS, 704,960 for the RB, and 766,888 for the EG models. Cortical bone mesh was modified to match the appropriate bone width for the chosen implant using Meshmixer. Building upon the HS model, the Siebert Class 1 defect mesh in the RB model was created by modifying the cortical bone layer in Autodesk Meshmixer. The residual ridge defect commonly seen clinically on the labial aspect of the implant was created (Fig 2). The EG model or immediate implant placement and grafting socket model was developed from the dentate CBCT dataset.

Graph: Fig 2 FEA model, loading protocols, measurement sites, and implant design.

A tissue-level root-form endosseous dental implant (Straumann Standard Plus Tissue Level RN) was used as an implant model similar to previous studies.[[33]] The model was generated using FreeCAD. (https://www.freecadweb.org/) [[35]] The implant had a 0.2 mm thread depth and 1.2mm pitch, carrying a 3.9 mm inner diameter, and osseous height of 10 mm and 1.8 mm platform height, from bone crest, for the tissues. It was modeled with cervical diameter of 4.8 mm (regular neck, RN, Straumann), fixture diameter of 4.1 mm, and a 5.5 mm in height anti-rotation abutment component. The analysis was focused on the one-piece Zir implant and therefore the abutment and implant were unified as one solid unit.[[36]] For the HS and RB models, the implant fixture was optimally positioned by taking into consideration buccal-palatal width, angulation, and expected soft tissue thickness & emergence profile (Fig 2). For the EG model, the implant fixture was placed within the tooth socket to match the angulation and emergence profile of other models while keeping the implant engagement to the palatal bone socket wall.[[25], [28], [37]] The remaining void in the socket between the implant and the labial bone was filled with mesh pieces to simulate the grafted material (Fig 2).

Each bone model underwent a final mesh refinement in 3-Matic (Materialise NV) for further triangle reduction and removal of unconnected, floating bone segments. The wrap operation in 3-Matic was used to remove sharp edges within each surface, a requirement to allow the mesh to be used in a FEA. Each model's bone parts and implant were then imported and combined in Solidworks (DS Solidworks Corp). The physical properties were applied and the simulations were carried out in Solidworks.[[38]] Cortical, trabecular, and morselized cancellous bone material properties were applied to the appropriate components (Table 1). The implant materials used were Yttria-stabilized Zirconia (Y-TZP) and Titanium alloy (Ti-6Al-4V), created in two separate study simulations for each model. Fixations were applied at the palatal direction of the abutment portion of the model. Then for each study comparing the two materials, two loading conditions were applied: 200 N force directed down the long axis and 178 N force directed along the implant labial axial wall, implant-crown margins, and incisal surface (Fig 3A and 3B). The force values were simulated normal biting force.[[40]] The oblique loading is a component of axial loading and perpendicular to long axis force application resulting in torque. The three regions (Fig 3A and 3B) were identified to simulate the sum of contact points for a single unit crown in an Angle's Class I occlusion. All the models have an assumption of bonded contact at all bone interfaces as well as at the implant-bone interfaces presuming complete osseointegration. The von Mises stress and equivalent strain distribution along the implant-bone interface were quantitatively reviewed. Von Mises stress and equivalent strain means were obtained from measurements of 10 (n = 10) fixated sensors across all models within the designated regions (Figs 2 and 3). From the sampled data in cortical and trabecular bone at implant interface across the 3 models, 2 loading scenarios, and 2 implant materials means and respective standard deviations were calculated in Excel (Microsoft Excel). To compare the biomechanical significance between Ti and Zir implants one-tailed p-values with α = 0.05 were computed from the 10 sample points obtained from regions 1–4 (Fig 2) for von Mises stress and equivalent strain.

Graph: Fig 3 FEA model and loading.(A) HS Model Stress and Strain Distributions, (B) RB Model Stress and Strain Distributions, and (C) EG Model Stress and Strain Distributions.

Graph

Table 1 List of materials used in this finite element analysis, their properties and numeric values used for the simulation*.

Y-TZPTi-6Al-4VCortical BoneTrabecular BoneMorselized Cancellous Bone53
Elastic Modulus (MPa)22000052104800.3171700042000910053
Poisson's Ratio0.31520.3170.4260.4100.253
Mass Density (kg/m^3)6000524428.7871900827023401
Tensile Strength (MPa)7455210507524
Compressive Strength (MPa)220052848–1080513041043
Yield Strength (MPa)30052827.3708874964.83501

1 *The empty cells indicate missing values that were not found in the literature and were not required for the FEA to compute. When available, the values were used to improve the accuracy of the model.

Results

Data measurements along the implants surfaces at the designated regions (Fig 2) from the FEA analysis showed that overall a Zir implant has higher internal von Mises stress (MPa) and lower equivalent strain when compared to its Ti counterpart (Table 2 and Table 3; Fig 4). Concurrently, for the HS model there was p-value significance of less than 0.01 for the 178 N oblique loading. For 200 N axial loading, Zir had significantly higher stress labially with p-value of 0.03 and less than 0.01 palatally. Strain was significantly lower for Zir on labial and palatal aspects with p-values less than 0.01. The RB model's Ti implant surface had significantly higher equivalent strain for 178 N oblique and 200 N axial loading, with labial and palatal p-value comparisons less than 0.01. For this model the labial internal von Mises stress difference was not significant for either loading conditions, p-values greater than 0.05. While the palatal for both loading conditions was greater for Zir, p-value less than 0.01 for 178 N oblique and 0.01 for 200 N axial. For the EG model, 178 N oblique loading showed higher von Mises stress for Zir where labial p-value was less than 0.01 and palatal was 0.04. Long axial loading had p-value less than 0.01 labially but Ti had higher stress in that region. Palatal stress difference was not significant (p-value 0.07). With respect to internal strain, Ti exhibited higher strain for both loading condition and it was significantly different compared to Ti with p-values less than 0.01 for all fixture regions.

Graph: Fig 4 Model loading at 178 N oblique for HS (4A), RB (4B), and EG (4C); at 200 N vertical for HS (4D), RB (4E), and EG (4F); and the summary of the two loading protocols (4G) and (4H).

Graph

Table 2 Stresses and strains computed from three FEA models at 178 N oblique loading.

HS MODEL
Stress (MPa)Strain
Ti-6Al-4VY-TZPP-valueTi-6Al-4VY-TZPP-value
LabialImplant (S4)39.18±9.8443.18±8.130.0004LabialImplant (S4)3.51E-04±2.59E-051.80E-04±6.55E-064E-10
Cortical (S4)13.3±5.9110.4±3.820.002Cortical (S4)7.87E-04±1.53E-046.09E-04±8.40E-052E-05
Trabecular (S2)3.90±1.024.05±1.160.01Trabecular (S2)1.39E-03±5.01E-041.43E-03±5.40E-040.01
PalatalImplant (S3)22.1±7.3625.58±7.332E-05PalatalImplant (S3)1.94E-04±1.73E-051.06E-04±4.63E-063E-09
Cortical (S3)6.15±4.514.28±2.470.009Cortical (S3)3.31E-04±1.05E-042.36E-04±5.49E-050.0001
Trabecular (S1)0.22±0.06380.19±0.08090.2Trabecular (S1)1.18E-04±3.38E-051.12E-04±9.81E-060.3
RB MODEL
Stress (MPa)Strain
Ti-6Al-4VY-TZPP-valueTi-6Al-4VY-TZPP-value
LabialImplant (S4)45.93±9.2147.2±10.30.08LabialImplant (S4)3.89E-04±2.46E-051.89E-04±1.22E-051E-12
Cortical (S4)48.6±18.845.9±20.70.04Cortical (S4)1.71E-03±6.51E-041.46E-03±9.04E-040.008
Trabecular (S2)0.311±0.2250.227±0.130.01Trabecular (S2)1.77E-04±9.77E-051.30E-04±5.61E-050.003
PalatalImplant (S3)20.51±10.622.689±10.40.0003PalatalImplant (S3)1.92E-04±3.86E-059.97E-05±1.63E-052E-07
Cortical (S3)10.24±7.219.11±5.220.07Cortical (S3)6.86E-04±2.57E-045.98E-04±1.75E-040.004
Trabecular (S1)0.38±0.1290.36±0.1340.01Trabecular (S1)1.50E-04±4.92E-051.65E-04±4.72E-050.0002
EG MODEL
Stress (MPa)Strain
Ti-6Al-4VY-TZPP-valueTi-6Al-4VY-TZPP-value
LabialImplant (S2)55.41±4.5657.34±5.208E-06LabialImplant (S2)4.75E-04±2.51E-052.36E-04±1.27E-052E-13
Cortical (S3)3.44±0.5443.17±0.9970.1Cortical (S3)1.44E-04±2.06E-051.31E-04±4.36E-050.09
Graft (S2)0.566±0.1240.429±0.08441E-06Graft (S2)4.44E-03±7.55E-043.36E-03±5.10E-041E-07
PalatalImplant (S1)40.46±1.9542.18±4.130.04PalatalImplant (S1)3.24E-04±8.39E-061.68E-04±8.17E-067E-19
Cortical (S1)11.5±5.656.83±2.870.0002Cortical (S1)5.24E-04±3.11E-053.18E-04±2.80E-058E-17

Graph

Table 3 Stresses and strains computed from three FEA models at 200 N loading along implant axis.

HS MODEL
Stress (MPa)Strain
Ti-6Al-4VY-TZPP-valueTi-6Al-4VY-TZPP-value
LabialImplant (S4)17.88±2.5218.81±2.960.03LabialImplant (S4)1.41E-04±9.37E-066.89E-05±5.71E-064E-12
Cortical (S4)9.48±2.538.26±2.150.0003Cortical (S4)5.12E-04±5.31E-054.50E-04±4.32E-055E-06
Trabecular (S2)4.32±0.9664.12±1.050.002Trabecular (S2)1.54E-03±5.52E-041.47E-03±5.52E-044E-05
PalatalImplant (S3)12.0±3.512.7±3.50.006PalatalImplant (S3)1.03E-04±6.54E-065.11E-05±2.71E-062E-11
Cortical (S3)8.44±3.777.76±2.940.04Cortical (S3)4.48E-04±6.88E-054.15E-04±5.20E-050.0001
Trabecular (S1)0.85±0.5020.78±0.5720.02Trabecular (S1)5.12E-04±6.46E-054.88E-04±9.00E-050.01
RB MODEL
Stress (MPa)Strain
Ti-6Al-4VY-TZPP-valueTi-6Al-4VY-TZPP-value
LabialImplant (S4)23.15±9.9023.11±10.620.5LabialImplant (S4)2.24E-04±4.25E-051.05E-04±2.11E-051E-08
Cortical (S4)39.9±16.537.7±17.60.04Cortical (S4)1.30E-03±5.28E-041.13E-03±7.10E-040.01
Trabecular (S2)0.394±0.2080.273±0.110.004Trabecular (S2)2.17E-04±9.33E-051.42E-04±3.78E-050.001
PalatalImplant (S3)12.02±3.2312.82±3.680.01PalatalImplant (S3)1.02E-04±6.02E-065.15E-05±4.00E-061E-13
Cortical (S3)8.22±2.797.50±2.270.004Cortical (S3)4.57E-04±9.74E-054.18E-04±7.65E-050.0001
Trabecular (S1)0.78±0.2710.73±0.2830.002Trabecular (S1)4.13E-04±8.60E-053.86E-04±8.65E-050.00004
EG MODEL
Stress (MPa)Strain
Ti-6Al-4VY-TZPP-valueTi-6Al-4VY-TZPP-value
LabialImplant (S2)18.31±1.0918.13±1.028E-05LabialImplant (S2)1.54E-04±6.95E-067.27E-05±3.19E-067E-14
Cortical (S3)0.639±0.1960.609±0.2640.2Cortical (S3)2.69E-05±8.52E-062.56E-05±1.17E-050.2
Graft (S2)0.148±0.03180.105±0.02093E-07Graft (S2)1.16E-03±1.96E-048.26E-04±1.29E-044E-08
PalatalImplant (S1)15.80±1.8416.24±1.160.07PalatalImplant (S1)1.36E-04±1.14E-056.51E-05±3.18E-063E-10
Cortical (S1)3.82±3.222.13±1.920.001Cortical (S1)1.50E-04±7.34E-068.17E-05±5.53E-069E-15

The FEA models' simulated bone samples of von Mises stress and equivalent strain were produced by loading and implant fixture of Ti and Zir vertically at 200 N and obliquely at 178 N. The distributions of von Mises stress and equivalent strain mean, standard deviation, and p-values were computed (Table 2 and Table 3; and Fig 4). From the HS model under 178 N oblique loading the labial aspect of the cortical bone showed a mean von Mises stress of 13.3±5.91 MPa for Ti and 10.4±3.82 MPa for Zir while on the palatal aspect 6.15±4.51 MPa and 4.28±2.47 MPa were observed respectively for Ti and Zir, both regions with a p-value of less than 0.01. For HS model with 178 N oblique loading labial area equivalent strain mean was found to be 7.87x10-4±1.53x10-4 for Ti while 6.09x10-4±8.40x10-5 for Zir. Whilst palatal region for Ti exhibited mean equivalent strain of 3.31x10-4±1.05x10-4 and 2.36x10-4±5.49x10-5 for Zir. Labial and palatal equivalent strain p-values were less than 0.01. Trabecular bone exhibited overall lower mean von Mises stress and equivalent strain compared to its cortical counterpart (Tables 2 and 3). Trabecular bone von Mises stress p-value for the labial aspect was 0.01 and palatal was 0.25 –implying significant difference between Ti and Zir, favoring Ti, in compression but no significance in tension. Trabecular bone equivalent strain labial region p-value was 0.01 implying statistical significance favoring Ti, while the palatal region p-value was 0.28 showing no statistical significance. With respect to long axial loading at 200 N, the mean von Mises stress at the labial cortical interface was 9.48±2.53 MPa for Ti and 8.26±2.15 MPa for Zir while on the palatal cortical interface was 8.44±3.77 MPa for Ti and 7.76±2.94 MPa for Zir. Labial p-value between the two materials was found to be less than 0.01and the palatal p-value was 0.04 showing statistical significance in von Mises stress between Ti and Zir. Analysis of the cortical labial equivalent strain showed a mean of 5.12x10-4±5.31x10-5 for Ti and 4.50x10-4±4.32x10-5 for Zir while the palatal aspect had a mean of 4.48x10-4±6.88x10-5 and 4.15x10-4±5.20x10-5, respectively. The equivalent strain p-values were less than 0.01 for labial and palatal aspects–implying statistical significance in both regions between Ti and Zir. With respect to the trabecular bone von Mises stress labial p-value was calculated to be less than 0.01 and palatal 0.02, while the equivalent strain p-values were both less than 0.01. There is statistical significance between Ti and Zir within trabecular bone at the labial and palatal aspects of the HS model that favor Zir.

The RB model showed cortical mean von Mises stress on the labial aspect under 178 N oblique loading of 48.6±18.8 MPa for Ti and 45.9±20.7 MPa for Zir, while on the palatal aspect those values were 10.24±7.21 MPa and 9.11±5.22 MPa respectively. The equivalent strain was found to be 1.71x10-3±6.51x10-4 and 1.46x10-3±9.04x10-4 respectively for the palatal region, and 6.86x10-4±2.57x10-4 and 5.98x10-4±1.75x10-4 respectively. The calculated von Mises stress p-value for the labial region was 0.04 and the palatal 0.07, while the equivalent strain comparison has a p-value less than 0.01 for the labial and palatal. For the trabecular bone, the von Mises stress p-value between Ti and Zir was calculated to be 0.01 for the labial and palatal regions–Ti had significantly higher stress than Zir. P-value calculated for trabecular bone equivalent strain was less than 0.01 for labial and palatal–Ti had significantly higher strain than Zir on the labial aspect but the opposite occurred for the palatal. With respect to the 200 N long axial loading condition, the von Mises stress in labial cortical bone was 39.9±16.5 MPa for Ti and 37.7±17.6 MPa for Zir, while the palatal region showed 8.22±2.79 MPa and 7.50±2.27 MPa, respectively. Associated p-values for stress within the labial region was found to be 0.04 and less than 0.01 for the palatal–Ti showed significantly higher stress in both regions. Equivalent strain wise a similar pattern appeared–labial cortical mean was 1.30x10-3±5.28x10-4 for Ti and 1.13x10-3±7.10x10-4 for Zir, while palatal was 4.57x10-4±9.74x10-5 for Ti and 4.18x10-4±7.65x10-5 for Zir. P-values were found to be 0.01 and less than 0.01 respectively. In trabecular bone, the von Mises stress as well as equivalent strain were significantly higher for Ti in both labial and palatal aspects–p-values both less than 0.01 for stress and strain.

The EG model overall favored Zir where Ti exhibited higher mean values in von Mises stress as well as equivalent strain. Under oblique loading at 178 N the labial cortical bone adjacent to the graft interface had a mean von Mises stress of 3.44±0.544 MPa for Ti and 3.17±0.997 MPa for Zir; however, with a p-value of 0.09 this was not statistically significant. The palatal region interfacing the implant surface had a mean von Mises stress of 11.50±5.65 MPa for Ti and 6.83±2.87 MPa for Zir–p-value of less than 0.01 was calculated signifying statistical significance. The stress analysis at the implant-graft interface showed a mean of 0.566±0.124 MPa for Ti and 0.429±0.0844 MPa, calculated p-value was less than 0.01. Equivalent strain wise the labial aspect had a mean of 1.44x10-4±0.206x10-4 for Ti and 1.31x10-4±0.436x10-4, and palatally 5.24x10-4±0.311x10-4 for Ti and 3.18x10-4±0.280x10-4 for Zir. Calculated p-values from the sample points were 0.09 and less than 0.01, respectively. At the graft-implant interface the equivalent strain for Ti was 4.44x10-3±0.755x10-3 and 3.36x10-3±0.51x10-3 for Zir with p-value of less than 0.01. For 200 N long axial loading the same pattern of stress and strain was observed such that Zir was favored over Ti.

Discussion

Zir dental implants have recently been introduced and advocated for the esthetic zone.[[1], [3], [6]] This study was one of the first to apply FEA models to represent the use of a Zir single implant in three common clinical scenarios in the esthetic zone; the healed socket or edentulous area (HS model), in the reduced bone width area (RB model), and in the extraction socket with grafting (EG model). Since the modulus of elasticity of Zir is approximately twice that of Ti, it was hypothesized that the stress/strain distribution might have been different. The results suggest that the stress/strain distributions were very similar between the two materials but were significantly different. Zir demonstrated lower von Mises stress and equivalent strain mean values in all measurements across all models' cortical bone. These results are similar to a previous FEA study comparing one-piece Zir implants with two-piece Zir/Ti abutment/implants, and two-piece Ti/Ti abutment/implants.[[22]] The higher elastic modulus of Zir provides for a "stiffer" material that dissipates more stress within its internal structure and equivalently lower displacement (ie. strain) of surrounding biologic structures–reflected in the implant's high internal von Mises stress and low equivalent strain (Fig 4). This shows that use of Zir in the three scenarios studied can potentially improve clinical outcomes and reduce peri-implant bone loss especially in the cervical area where bone resorption occurs more readily.[[21]] A recent FEA implant model study pointed out a phenomenon that zirconia received higher stress value compared to titanium.[[42]] The authors further suggested that there may be a direct relationship between Young's modulus of the material and the stress transferred to the implant, that a more rigid implant absorbs more stress.[[42]] This may partly explain why in our study, there was less stress in the peri-implant bone in the Zir implant compared to the Ti implant.

In this study, the CBCT scans were used to provide microstructures of the trabecular bone mimicking the natural bone similar to others.[[21], [43]–[46]] These CBCT based FEA models allow von Mises stress and equivalent stain distributions to be simulated more realistically than other geometrically simplified FEA models.[[10], [21], [46]] More importantly, the three FEA models, the healed edentulous site (HS model), reduced bone width edentulous site (RB model), and immediate implant placement with grafting (EB model), allow for a direct comparison between Zir and Ti implants as well as the examination of Zir in common scenarios. Zir implants provide an overall more favorable peri-implant bone stress/strain distribution. This suggests that in non-grafted edentulous single tooth implant site, Seibert defects that experience compression, and the grafted extraction socket site, Zir implants are mechanically superior to Ti implants in terms of occlusal load distribution in the oblique and axial direction. The EG model results suggest Zir implants may allow less graft resorption/labial bone resorption for immediate placement as well as their use in sites with periodontal defects. However, Zir implants and Ti implants in the immediate placement cases may also be selected based on the other clinical scenarios/esthetic considerations. A recent systematic review suggest that Zir implants have similar short time survival rate as Ti implants, ~92%.[[1]] Another systematic review showed Zir may be an alternative implant material for Ti.[[3]] On the contrary, one systematic review pointed out the low evidence level of Zir implant studies and advised caution in clinical use of Zir implants.[[6]] In the immediate implant scenario higher insertion torque, representing good primary stability, may be easier to achieve in the Ti implant due to their lower elastic modulus. Therefore, fracture of Zir implants can be an issue especially in a two-piece or small diameter Zir implants when higher torque is applied.[[7]]

This study shares similar limitations with other FEA studies. These limitations include the limited study designs, material properties, occlusal loads, individual variations.[[10]] This study only targets three scenarios designed to compare two materials with one-piece implant designs. The results may not be applicable to two-piece implants.[[22], [46]] The material properties used here are the ones commonly prescribed in FEA studies.[[22]] The particular properties may be an average and can vary in individual patients. The occlusal loads applied here in the FEA models is to provide the general idea of how stress/strain distributions would be.[[47]] The models unfortunately may not represent an individual's physiologic bite and consequently occlusal loading.[[48]] Assessing various patient-specific bone properties and occlusal loading is unrealistic. Finally, while the CBCT scans from a human volunteer were chosen here, the scans may not be a representative of a wide-range of clinical cases.[[43]–[46]] It is possible to define the properties of bone directly based off of the CT intensity and use more subjects; however, the technique requires more specific computational modeling and the additional value of the information for stress evaluation is not clear.[[49]]

Zir implants demonstrate a promising potential when compared to Ti implants, and based on the conducted FEA analysis the null hypothesis has been rejected. However, clinicians should consider other clinical parameters, such as esthetic demands, insertion torque, clinical occlusion, parafunctional habits, position of the site, type of edentulous areas, etc.[[50]–[54]] A prospective clinical study with a large population will be needed to prove that Zir implants are a good clinical alternative to Ti implants. Furthermore, a supplemental in vitro analysis to support the computational models is subsequently employed to verify the in silico results–FEA models are meant to show how materials behave and not be used for quantitative experimental data.

Conclusions

In general Zir implants behave more favorable then Ti implants in terms of peri-implant stress distributions. Three different FEA models, healed edentulous site (HS), vertical periodontal defect under compression (RB), and immediate extraction with bone grafting site (EG), mimicking the common clinical scenarios suggested the following conclusion. Due to the stiffness of the material and its intrinsically high elastic modulus, Zir implants transmit less von Mises stress and induce lower equivalent strain to the peri-implant bone compared to Ti implants. This was statistically significant when cervical cortical bone (HS and RB models) and graft (EG model) were evaluated for one-piece Zir implant compared to its Ti counterpart. Therefore, the peri-implant bone surrounding Zir implants may be less prone to mechanically induced biologic peri-implant bone resorption. Zir implants may be considered not only due to its esthetic properties, but also due to the stress modulation properties of the material.

Supporting information

S1 Data.

(XLSX)

The authors thank Dr. Jennifer Wayne and her laboratory for the assistance in the FEA modeling software. Thanks also to the faculty and staff of the VCU Center of Digital Dentistry.

List of abbreviations

• Zir

  • zirconia

• FEA

  • finite element analysis

• CBCT

  • cone-beam computed tomography

• Ti

  • titanium alloy

• STL

  • standard tessellation language

• HS

  • healed socket with ideal osseous dimension for a single tooth implant

• RB

  • reduced bone width model

• EG

  • extraction site with grafting

• Y-TZP

  • yttria-stabilized Zirconia

• Ti6Al4V

  • Titanium alloy Ti6Al4V

• N

• Newton

Footnotes 1 The authors have declared that no competing interests exist. References Hashim D, Cionca N, Courvoisier DS, Mombelli A. A systematic review of the clinical survival of zirconia implants. Clin Oral Investig. 2016; 20: 1403–1417. doi: 10.1007/s00784-016-1853-9, 27217032 2 Denry I, Kelly JR. State of the art of zirconia for dental applications. Dent Mater. 2008; 24: 299–307. doi: 10.1016/j.dental.2007.05.007, 17659331 3 Andreiotelli M, Wenz HJ, Kohal R-J. Are ceramic implants a viable alternative to titanium implants? A systematic literature review. Clin Oral Implants Res2009; 20Suppl 4:32–47. 4 Piconi C, Maccauro G. Zirconia as a ceramic biomaterial. Biomaterials1999; 20:1–25. doi: 10.1016/s0142-9612(98)00010-6, 9916767 5 Scarano A, Piattelli M, Caputi S, Favero GA, Piattelli A. Bacterial adhesion on commercially pure titanium and zirconium oxide disks: an in vivo human study. J Periodontol2004; 75:292–6. doi: 10.1902/jop.2004.75.2.292, 15068118 6 Depprich R, Naujoks C, Ommerborn M, Schwarz F, Kübler NR, Handschel J. Current findings regarding zirconia implants. Clin Implant Dent Relat Res2014; 16:124–37. doi: 10.1111/j.1708-8208.2012.00454.x, 24533568 7 Karl M, Scherg S, Grobecker-Karl T. Fracture of Reduced-Diameter Zirconia Dental Implants Following Repeated Insertion. Int J Oral Maxillofac Implants.2017; 32:971–5. doi: 10.11607/jomi.5592, 28906500 8 Sanon C, Chevalier J, Douillard T, Cattani-Lorente M, Scherrer SS, Gremillard L. A new testing protocol for zirconia dental implants. Dent Mater2015; 31:15–25. doi: 10.1016/j.dental.2014.09.002, 25262212 9 Grassi FR, Capogreco M, Consonni D, Bilardi G, Buti J, Kalemaj Z. Immediate occlusal loading of one-piece zirconia implants: five-year radiographic and clinical evaluation. Int J Oral Maxillofac Implants2015; 30:671–80. doi: 10.11607/jomi.3831, 26009919 Geng JP, Tan KB, Liu GR. Application of finite element analysis in implant dentistry: a review of the literature. J Prosthet Dent2001; 85:585–98. doi: 10.1067/mpr.2001.115251, 11404759 Kwon B-G, Kim S-G. Finite Element Analysis of Different Bone Substitutes in the Bone Defects Around Dental Implants. Implant Dent2006; 15:254–64. doi: 10.1097/01.id.0000219864.33618.8b, 16966899 Schrotenboer J, Tsao Y-P, Kinariwala V, Wang H-L. Effect of Platform Switching on Implant Crest Bone Stress: A Finite Element Analysis. Implant Dent2009; 18:260–9. doi: 10.1097/ID.0b013e31819e8c1d, 19509536 Hsu M-L, Chang C-L. Application of Finite Element Analysis in Dentistry. Finite Element Analysis2010 p 44–60. Balkaya MC. Investigation of Influence of Different Implant Size and Placement on Stress Distribution With 3-Dimensional Finite Element Analysis. Implant Dent2014; 1. Holmgren EP, Seckinger RJ, Kilgren LM, Mante F. Evaluating Parameters of Osseointegrated Dental Implants Using Finite Element Analysis: A Two-Dimensional Comparative Study Examining the Effects of Implant Diameter, Implant Shape, and Load Direction. J Oral Implantol1998; 24:80–8. doi: 10.1563/1548-1336(1998)0240080:EPOODI2.3.CO;2, 9835834 Chun H-J, Chun H-J., Cheong S-Y., Han J-H., Heo S-J., Chung J-P., et al. Evaluation of design parameters of osseointegrated dental implants using finite element analysis. J Oral Rehabil2002; 29:565–74. doi: 10.1046/j.1365-2842.2002.00891.x, 12071926 Papavasiliou G, Kamposiora P, Bayne SC, Felton DA. Three-dimensional finite element analysis of stress-distribution around single tooth implants as a function of bony support, prosthesis type, and loading during function. J Prosthet Dent1996; 76:633–40. doi: 10.1016/s0022-3913(96)90442-4, 8957790 Petrie CS, Williams JL. Comparative evaluation of implant designs: influence of diameter, length, and taper on strains in the alveolar crest. A three-dimensional finite-element analysis. Clin Oral Implants Res2005; 16:486–94. doi: 10.1111/j.1600-0501.2005.01132.x, 16117775 Himmlová L, Dostálová T 'jana, Kácovský A, S. Influence of implant length and diameter on stress distribution: A finite element analysis. J Prosthet Dent2004; 91:20–5. doi: 10.1016/j.prosdent.2003.08.008, 14739889 Pellizzer EP, Lemos CAA, Almeida DAF, de Souza Batista VE, Santiago Júnior JF, Verri FR. Biomechanical analysis of different implant-abutments interfaces in different bone types: An in silico analysis. Materials Science and Engineering C2018; 90:645–50. doi: 10.1016/j.msec.2018.05.012, 29853134 Marcián P, Wolff J, Horáčková L, Kaiser J, Zikmund T, Borák L. Micro finite element analysis of dental implants under different loading conditions. Comput Biol Med.2018; 96:157–65. doi: 10.1016/j.compbiomed.2018.03.012, 29587150 Çaglar A, Bal BT, Karakoca S, Aydın C, Yılmaz H, Sarısoy S. Three-dimensional finite element analysis of titanium and yttrium-stabilized zirconium dioxide abutments and implants. Int J Oral Maxillofac Implants2011; 26:961–9., 22010077 Huynh-Ba G, Pjetursson BE, Sanz M, Cecchinato D, Ferrus J, Lindhe J, et al. Analysis of the socket bone wall dimensions in the upper maxilla in relation to immediate implant placement. Clin Oral Implants Res2010; 21:37–42. doi: 10.1111/j.1600-0501.2009.01870.x, 20070745 Gluckman H, Pontes CC, Du Toit J. Radial plane tooth position and bone wall dimensions in the anterior maxilla: A CBCT classification for immediate implant placement. J Prosthet Dent2018; 120:50–6. doi: 10.1016/j.prosdent.2017.09.005, 29195817 Bencharit S, Byrd WC, Hosseini B. Immediate placement of a porous-tantalum, trabecular metal-enhanced titanium dental implant with demineralized bone matrix into a socket with deficient buccal bone: a clinical report. J Prosthet Dent2015; 113:262–9. doi: 10.1016/j.prosdent.2014.09.022, 25702965 Bencharit S, Staffen A, Yeung M, Whitley D 3rd, Laskin DM, Deeb GR. In Vivo Tooth-Supported Implant Surgical Guides Fabricated With Desktop Stereolithographic Printers: Fully Guided Surgery Is More Accurate Than Partially Guided Surgery. J Oral Maxillofac Surg2018; 76:1431–9. doi: 10.1016/j.joms.2018.02.010, 29550378 Pianykh OS. Digital Imaging and Communications in Medicine (DICOM): A Practical Introduction and Survival Guide. Springer Science & Business Media; 2009.p 3–14,363–390. Edelmann AR, Hosseini B, Byrd WC, Preisser JS, Tyndall DA, Nguyen T, et al. Exploring Effectiveness of Computer-Aided Planning in Implant Positioning for a Single Immediate Implant Placement. J Oral Implantol2016; 42:233–9. doi: 10.1563/aaid-joi-D-15-00056, 26652644 Wang W, Gao C-H, Lin S-F, Yin G-F. Modeling Technologies of FEA Grid Services and Grid Information Management Structure for FEA Grid Portal. 2009 WRI World Congress on Software Engineering. 2009. doi: 10.1109/wcse.2009.110 Schropp L, Wenzel A, Kostopoulos L, Karring T. Bone healing and soft tissue contour changes following single-tooth extraction: a clinical and radiographic 12-month prospective study. Int J Periodontics Restorative Dent2003; 23:313–23., 12956475 Koopaie M, Kolahdouz S. Three-dimensional simulation of human teeth and its application in dental education and research. Med J Islam Repub Iran2016; 30:461., 28491836 Spradley JP, Pampush JD, Morse PE, Kay RF. Smooth operator: The effects of different 3D mesh retriangulation protocols on the computation of Dirichlet normal energy. Am J Phys Anthropol2017; 163:94–109. doi: 10.1002/ajpa.23188, 28218399 Chun H-J, Park D-N, Han C-H, Heo S-J, Heo M-S, Koak J-Y. Stress distributions in maxillary bone surrounding overdenture implants with different overdenture attachments. J Oral Rehabil2005; 32:193–205. doi: 10.1111/j.1365-2842.2004.01407.x, 15707430 Degerliyurt K, Simsek B, Erkmen E, Eser A. Effects of different fixture geometries on the stress distribution in mandibular peri-implant structures: a 3-dimensional finite element analysis. Oral Surg Oral Med Oral Pathol Oral Radiol Endod2010; 110:e1–11. Gayer D, O'Sullivan C, Scully S, Burke D, Brossard J, Chapron C. FreeCAD visualization of realistic 3D physical optics beams within a CAD system-model. Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VIII2016. doi: 10.1117/12.2233268 Sevimay M, Turhan F, Kiliçarslan MA, Eskitascioglu G. Three-dimensional finite element analysis of the effect of different bone quality on stress distribution in an implant-supported crown. J Prosthet Dent2005; 93:227–234. doi: 10.1016/j.prosdent.2004.12.019, 15775923 Hosseini B, Byrd WC, Preisser JS, Khan A, Duggan D, Bencharit S. Effects of Antibiotics on Bone and Soft-Tissue Healing Following Immediate Single-Tooth Implant Placement Into Sites With Apical Pathology. J Oral Implantol2015; 41:e202–11. doi: 10.1563/aaid-joi-D-14-00075, 25076118 Akin JE. Finite Element Analysis Concepts: Via SolidWorks. World Scientific; 2010. Shih R. Introduction to Finite Element Analysis Using SOLIDWORKS Simulation 2016. SDC Publications; 2016. Ferrario VF, Sforza C, Serrao G, Dellavia C, Tartaglia GM. Single tooth bite forces in healthy young adults. J Oral Rehabil2004; 31:18–22. doi: 10.1046/j.0305-182x.2003.01179.x, 15125591 Ferrario VF, Sforza C, Zanotti G, Tartaglia GM. Maximal bite forces in healthy young adults as predicted by surface electromyographyJ Dent. 2004; 32:451–7. doi: 10.1016/j.jdent.2004.02.009, 15240063 Pérez-Pevida E, Brizuela-Velasco A, Chávarri-Prado D, Jiménez-Garrudo A, Sánchez-Lasheras F, Solaberrieta-Méndez E, et al. Biomechanical Consequences of the Elastic Properties of Dental Implant Alloys on the Supporting Bone: Finite Element Analysis. Biomed Res Int. 2016; 2016:1850401. doi: 10.1155/2016/1850401, 27995137 HR Ali S, Ali SHR. 3D-Design of Artificial Mandible Using CBCT, CAD and FEA. Current Trends in Biomedical Engineering & Biosciences2017; 8. doi: 10.19080/ctbeb.2017.08.555732 Zhang Y, Tehrani JN, Wang J. A Biomechanical Modeling Guided CBCT Estimation Technique. IEEE Trans Med Imaging2017; 36:641–52. doi: 10.1109/TMI.2016.2623745, 27831866 Bujtar P, Simonovics J, George S, Varadi K, Christopher A, Szucs A, et al. Emerging computer bioengineering technologies 1: Scanning (CBCT versus MSCT), modelling (bio-CAD) and functional simulation (FEA): Steps towards validation. Br J Oral Maxillofac Surg2014; 52: e71. Lee K-S, Shin S-W, Lee S-P, Kim J-E, Kim J-H, Lee J-Y. Comparative Evaluation of a Four-Implant-Supported Polyetherketoneketone Framework Prosthesis: A Three-Dimensional Finite Element Analysis Based on Cone Beam Computed Tomography and Computer-Aided Design. Int J Prosthodont2017; 30:581–5. doi: 10.11607/ijp.5369, 29095963 Borges Radaelli MT, Idogava HT, Spazzin AO, Noritomi PY, Boscato N. Parafunctional loading and occlusal device on stress distribution around implants: A 3D finite element analysis. J Prosthet Dent. 2018; doi: 10.1016/j.prosdent.2017.12.023, 29724560 Winter W, Klein D, Karl M. Effect of model parameters on finite element analysis of micromotions in implant dentistry. J Oral Implantol2013; 39:23–9. doi: 10.1563/AAID-JOI-D-11-00221, 23402356 Dragomir-Daescu D, Op Den Buijs J, McEligot S, Dai Y, Entwistle RC, Salas C, et al. Robust QCT/FEA models of proximal femur stiffness and fracture load during a sideways fall on the hip. Ann Biomed Eng2011; 39:742–55. doi: 10.1007/s10439-010-0196-y, 21052839 Liddelow G, Klineberg I. Patient-related risk factors for implant therapy. A critique of pertinent literature. Aust Dent J2011; 56:417–26. doi: 10.1111/j.1834-7819.2011.01367.x, 22126353 Ubhi K. Recognizing biological risk factors for implant failure. Dental Nursing2007; 3:336–42. Nassani M, El-Anwar M, El-Zawahry M, Ibraheem E, ElGabry H. New dental implant selection criterion based on implant design. Eur J Dent2017; 11:186. doi: 10.4103/1305-7456.208432, 28729790 Bryington M, De Kok IJ, Thalji G, Cooper LF. Patient Selection and Treatment Planning for Implant Restorations. Dent Clin North Am2014; 58:193–206. doi: 10.1016/j.cden.2013.09.009, 24286653 Brodt MD, Swan CC, Brown TD. Mechanical behavior of human morselized cancellous bone in triaxial compression testing. J Orthop Res1998; 16:43–9. doi: 10.1002/jor.1100160108, 9565072

By Georgi Talmazov; Nathan Veilleux; Aous Abdulmajeed and Sompop Bencharit

Reported by Author; Author; Author; Author

Titel:
Finite element analysis of a one-piece zirconia implant in anterior single tooth implant applications.
Autor/in / Beteiligte Person: Talmazov, Georgi ; Veilleux, Nathan ; Abdulmajeed, Aous ; Bencharit, Sompop
Link:
Zeitschrift: PLoS ONE, Jg. 15 (2020), Heft 2, S. e0229360
Veröffentlichung: Public Library of Science (PLoS), 2020
Medientyp: academicJournal
ISSN: 1932-6203 (print)
DOI: 10.1371/journal.pone.0229360
Schlagwort:
  • Medicine
  • Science
Sonstiges:
  • Nachgewiesen in: Directory of Open Access Journals
  • Sprachen: English
  • Collection: LCC:Medicine ; LCC:Science
  • Document Type: article
  • File Description: electronic resource
  • Language: English

Klicken Sie ein Format an und speichern Sie dann die Daten oder geben Sie eine Empfänger-Adresse ein und lassen Sie sich per Email zusenden.

oder
oder

Wählen Sie das für Sie passende Zitationsformat und kopieren Sie es dann in die Zwischenablage, lassen es sich per Mail zusenden oder speichern es als PDF-Datei.

oder
oder

Bitte prüfen Sie, ob die Zitation formal korrekt ist, bevor Sie sie in einer Arbeit verwenden. Benutzen Sie gegebenenfalls den "Exportieren"-Dialog, wenn Sie ein Literaturverwaltungsprogramm verwenden und die Zitat-Angaben selbst formatieren wollen.

xs 0 - 576
sm 576 - 768
md 768 - 992
lg 992 - 1200
xl 1200 - 1366
xxl 1366 -