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.

finite element analysis; mechanical properties; research

Resumo

O presente estudo avaliou a confiabilidade do MEF nos resultados de diferentes propriedades mecânicas (E e v) dos materiais. Dois modelos 3D de um canino superior com tratamento endodôntico, pino intrarradicular, núcleo em resina composta e restaurado com coroa metalocerâmica foram gerados de acordo com imagens de micro-CT. Dois grupos com diferentes valores de E e v para a cerâmica, coping metálico, cimento resinoso e resina composta foram criados. As propriedades dos materiais para o grupo GL foram baseadas em dados obtidos da literatura, enquanto para o grupo GIE a técnica de excitação por impulso foi utilizada. Carga de 180 N foi aplicada a 45° na superfície lingual do terço incisal do canino. Todos os modelos foram suportados pelo ligamento periodontal (x=y=z=0). Tensões equivalentes de von Mises (VMS) foram calculadas. Os valores de tensões revelaram diferença entre os grupos, ambos para a distribuição e valor das VMS. A cerâmica (GL: 5,966 MPa; GIE: 7,478 MPa), coping metálico (GL: 3,811 MPa; GIE: 0,973 MPa) e núcleo (GL: 4,771 MPa; GIE: 0,026 MPa) foram significativamente afetados. Em conclusão, este estudo mostrou que a determinação das propriedades mecânicas (E and v) dos materiais é essencial para a confiabilidade dos resultados do MEF.

Introduction

Finite element analysis (FEA) has been suggested as an efficient method to simulate a real clinical scenario and solve complex problems by evaluations of inner structural performance ^{1 2}. Thus, it has also been considered as an important method to foresee long-term failures ³.

FEA is a computational analysis that transforms a physical problem into a virtual model that presents mathematical equations to the finite element software. The virtual model containing several structures with specific mechanical properties is divided into small fragments to generate the finite element mesh ⁴. As result, the performance of materials, techniques, displacement and stress distribution under load are numerically simulated in the model ⁵.

The method is considered as a fast and low-cost technique that does not require a large sample size and avoids human-related failures ^{6 7}. Conversely, the reliability of results obtained by numeric methods should be carefully interpreted because it depends on the accurate simulation of the structure's geometry, mechanical properties of the materials, interface, support and load ⁸.

Young modulus (E) and Poisson ratio (v) of each simulated structure are the main mechanical properties to be used in static FEA. In general, the properties suggested by the literature are usually applied in the studies since complex and time-consuming laboratorial techniques are required to establish those values ⁹. Thus, most FEA studies are based on mechanical properties of materials with specific composition and unknown brand, which may influence the results of the research.

Dynamic impulse excitation technique ^{10 11} has been indicated to determine the mechanical properties of the materials since it determines the Young modulus and Poisson ratio of different materials by using nondestructive approaches. Although the pulse-echo ultrasound method is also suitable to determine E and v^{10 11}, it is a high cost not commonly used alternative.

The aim of this study was to evaluate the reliability of FEA on the results of different mechanical properties (E and v) of materials, using literature data and the impulse excitation technique. The work hypothesis assumed no difference in stress value and distribution between the models.

Material and Methods

Finite Element Model

Micro computed tomography (micro-CT) slices (SkyScan 1174v2; Bruker-microCT, Kontich, Belgium) were used to determine the 3D geometry of a maxillary human canine. The 3D model was generated using the Simpleware 4.1 software (Simpleware Ltd, Exeter, UK) and micro-CT images. All simulated structures (enamel, dentin, pulp and periodontal ligament) were included in the solid model and the 0.2-mm-thick periodontal ligament was designed by Boolean operations ¹².

According to a previous study ¹³, the following scenarios were simulated with CAD structures (SolidWorks Corp., Concord, MA, USA): (i) endodontic treatment with 14 mm working length; (ii) intracanal restoration with a prefabricated parallel sided fiber post (14 mm long and 1.5 mm diameter), preserving 4 mm of gutta-percha at the apical region; (iii) 50-μm-thick cementation line ; (iv) composite resin core to restore crown structure and simulated tooth preparation with 1 mm reduction in buccal and lingual surfaces and 2 mm in the incisal edge; (v) prosthetic restoration with a porcelain-fused-to-metal crown using a metal coping (0.5 mm thick and 0.5 mm porcelain) in the buccal and lingual surfaces and 1.5 mm on the incisal edge. The CAD images were transferred to the Simpleware 4.1 to create the 3D model.

Properties of Materials

Two groups were created with different mechanical property values of the materials (Young modulus {E} and Poisson ratio, {v}): GL - mechanical properties of all materials were based on literature data (Table 1); and GIE - mechanical properties of resin cement, composite resin core, metal coping and porcelain were determined by the impulse excitation technique while the properties of the remaining materials were similar to GL (Table 2).

Table 1
Mechanical properties of the group GL materials

• Refs: ¹⁴,¹⁵,¹⁶,¹⁷,¹⁸

Table 2
Mechanical properties of the group GIE materials

• Refs:¹⁷,¹⁸

Sonelastic (ATCP Physical Engineering, Ribeirão Preto, SP, Brazil) equipment was used to determine the mechanical properties of the materials using the impulse excitation technique. Three samples (4 mm wide, 2 mm high and 25 mm long) were fabricated for each material: 1 - porcelain (IPS InLine PoM, Ivoclar Vivadent AG, Schaan, Liechtenstein), 2 - resin cement (Rely X, 3M/ESPE, St Paul, MN, USA), 3 - metal alloy (Ni-Cr) (Vera Bond II, AalbaDent Inc., Fairfield, CA, USA), and 4 - composite resin (Z250, 3M/ESPE). After positioning in a specific holder, each sample was submitted to a short impact by a pulsator using mechanical vibration. A transducer captured the acoustic response and translated it into an electric sign in order to read the resonance frequencies. The Young modulus was determined by bending excitation of the sample while the Poisson ratio was obtained by torsion excitation. The steps were repeated for each sample until the last two measurement values became homogeneous, with a difference equal to or less than 2% from each other.

All materials were assumed as homogeneous, isotropic and linearly elastic; except for the glass fiber post (GFP) (Table 3), which was assumed as orthotropic, homogeneous and linearly elastic ¹⁹. GFP was assumed as orthotropic since its mechanical properties are different along the fiber (x) and normal directions (y and z). Thus, the mechanical characteristics of these materials were represented by the Young modulus in the 3 axes (Ex, Ey, Ez), while Poisson ratio (vxz, vxz, vyz) and shearing modulus (Gxz, Gxz, Gyz) in the orthogonal planes (xy, xz and yz) (Table 3) ¹⁹. All structures of the models were assumed as perfectly joined, which means no failure at bonding and interposition.

Table 3
Orthotropic mechanical properties of the glass fiber post

Finite Element Mesh

The finite element mesh was obtained with the Simpleware 4.1 software using linear tetrahedral elements type C3D4. All models presented 234,824 nodes and 1,243,290 elements. The mesh refinement was based on convergence analysis of 6% ²⁰.

Boundary Conditions and Load

The finite element mesh was transferred to the finite element software (Abaqus 6.10-EF1, Dassault Systèmes Simulus Corp., Providence, RI, USA) to simulate a static occlusal load of 180 N applied at 45° on the incisal third of the lingual surface of the canine tooth in all models ¹³. As boundary condition, the nodes of the periodontal ligament were fixed on the three axes of Cartesian plane (x, y and z), assuming values of x = y = z = 0 ²¹.

Results Calculation

Results were obtained with the finite element software using stress values (hot colors for higher stress and cold colors for lower stress) and numerical analysis. After loading, the maximum von Mises stress (VMS) was calculated. This criterion allows an efficient evaluation of failure pattern and damage occurrence under complex stress ^{2 22}.

Results

The qualitative analysis of stress values revealed difference in von Mises stress distribution between groups GL and GIE (Fig. 1). Higher von Mises stress concentration at the loading area was found for both groups. GIE concentrated more stress in the porcelain and transferred lower stress to the composite resin core than GL. Better von Mises stress distribution was observed in root dentin for group GL, although the numerical analysis showed similar maximum stress for both groups (GL: 1.614 MPa and GIE: 1.975 MPa). The remaining structures (resin cement, dentin, periodontal ligament, gutta percha and glass fiber post) exhibited similar VMS values and distribution (Table 4).

Figure 1
von Mises stress (VMS) distribution in the canine restored with a glass fiber post.

Table 4
von Mises stress values (VMS) of the groups

In quantitative analysis, the porcelain (GL: 5.966 MPa and GIE: 7.478 MPa), metal coping (GL: 3.811 MPa and GIE: 0.973 MPa) and composite resin core (GL: 4.771 MPa and GIE: 0.026 MPa) showed the highest differences in von Mises stress values (Fig. 1) between the groups of study.

Discussion

The present study evaluated the reliability of FEA on the results of different mechanical properties (E and v) of materials, using literature data and the impulse excitation technique. The work hypothesis was rejected since the stress values showed that different Young modulus and Poisson ratio changed von Mises stress distribution and value (Table 4).

In the present work the groups GL and GIE used the same dental materials, but they showed different mechanical properties. This difference in mechanical properties of the materials (E and v) showed that materials with higher Young modulus absorb higher stress; which agrees with literature ^{2 23}. It was found that the von Mises stress value of GIE porcelain was 1.253 times higher than GL and VMS of the GL metal coping was 3.917 times higher than GIE. Thus, comparing the results between the groups, it may be assumed that the GIE porcelain is more protective than GL since it avoided stress concentration on fragile structures like the composite resin core (GL: 4.771 MPa and GIE: 0.026 MPa). Considering the individual results of the groups, it may be suggested that the risk of core fracture in GL was 183.5 times higher than in GIE (Table 5).

Table 5
Ratio between the von Mises stress values (VMS) of the groups

The FEA reliability depends on using an appropriate software to solve the mathematical problem, the structures geometry, material properties, quality of finite element mesh, boundary and loading conditions ^{4 24}. Over the years, there was a wide range of commercial brands and mechanical properties for the same dental material in order to provide longevity, stability and better clinical performance. Thus, it is expected that the materials' properties used at the beginning of this century be different from currently used dental materials. Therefore, the authors suggest that it is important to determine the mechanical properties of materials used in the study to be able to extrapolate the results to a clinical situation.

It is important to highlight that the present research does not intend to criticize the previous FEA studies but enlighten the scientific community about the relevance of determining the Young modulus and Poisson ratio for each new dental material in order to avoid inaccurate analyses and false results. Finally, the findings of this study indicate that the definition of mechanical properties (E and v) of materials is essential for the FEA results reliability.

Acknowledgements

The authors wish to thank the support by FAPESP (Grant # 2013/23905-9), Laboratory of Biomechanical Studies in Prosthodontics and Implants at the Department of Dental Materials and Prosthodontics, and the Endodontics Research Laboratory at the Department of Restorative Dentistry, School of Dentistry of Ribeirão Preto, University of São Paulo, Brazil.

References

• ¹
  Santos-Filho PC, Veríssimo C, Soares PV, Saltarelo RC, Soares CJ, Marcondes Martins LR. Influence of ferrule, post system, and length on biomechanical behavior of endodontically treated anterior teeth. J Endod 2014;40:119-123.
• ²
  Veríssimo C, Simamoto Júnior PC, Soares CJ, Noritomi PY,. Santos-Filho PC Effect of the crown, post, and remaining coronal dentin on the biomechanical behavior of endodontically treated maxillary central incisors. J Prosthet Dent 2014;111:234-246.
• ³
  Vafaee F, Khoshhal M, Rezaei A, Sooltani F, Jalalzadeh M, Yalpaniyan A, et al.. Multiple in vitro analyses of fracture resistance in maxillary central incisors restored with fiber posts. Iran Endod J 2010;5:125-130.
• ⁴
  Tajima K, Chen KK, Takahashi N, Noda N, Nagamatsu Y, Kakigawa H. Three-dimensional finite element modeling from CT images of tooth and its validation. Dent Mater J 2009;28:219-226.
• ⁵
  Ausiello P, Apicella A, Davidson CL. Effect of adhesive layer properties on stress distribution in composite restorations - a 3D finite element analysis. Dent Mater 2002;18:295-303.
• ⁶
  Aggarwal S, Garg V. Finite element analysis of stress concentration in three popular brands of fiber posts systems used for maxillary central incisor teeth. J Conserv Dent 2011;14:293-296.
• ⁷
  Oliveira JS, Chaves CAL, Silva-Sousa YTC, Gomes EA. The finite element methods on oral rehabilitation: new trend for endodontically treated teeth. Braz Dent Sci 2014;17:12-19.
• ⁸
  Murakami N, Wakabayashi N. Finite element contact analysis as a critical technique in dental biomechanics: a review. J Prosthodont Res 2014;58:92-101.
• ⁹
  Della Bona A, Borba M, Benetti P, Duan Y, Griggs JA.Three-dimensional finite element modelling of all-ceramic restorations based on micro-CT. J Dent 2013;41:412-419.
• ¹⁰
  Borba M, de Araújo MD, de Lima E, Yoshimura HN, Cesar PF, Griggs JA, et al.. Flexural strength and failure modes of layered ceramic structures. Dent Mater 2011;27:1259-1266.
• ¹¹
  Passos SP, Kimpara ET, Bottino MA, Rizkalla AS, Santos GC Jr. Effect of ceramic thickness and shade on mechanical properties of a resin luting agent. J Prosthodont 2014;23:462-466.
• ¹²
  Eraslan O, Aykent F, Yucel MT, Akman S. The finite element analysis of the effect of ferrule height on stress distribution at post-and-core-restored all-ceramic anterior crowns. Clin Oral Investig 2009;13:223-227.
• ¹³
  Gomes ÉA, Gueleri DB, da Silva SR, Ribeiro RF, Silva-Sousa YT. Three-dimensional finite element analysis of endodontically treated teeth with weakened radicular walls restored with different protocols. J Prosthet Dent 2015;114:383-389.
• ¹⁴
  Akkayan B. An in vitro study evaluating the effect of ferrule length on fracture resistance of endodontically-treated teeth restored with fiber-reinforced and zirconia dowel systems. J Prosthet Dent 2004;92:155-162.
• ¹⁵
  Silva NR, Bonfante E, Rafferty BT, Zavanelli RA, Martins LL, Rekow ED, et al.. Conventional and modified veneered zirconia vs. metalloceramic: fatigue and finite element analysis. J Prosthodont 2012;21:433-439.
• ¹⁶
  Asmussen E, Peutzfeld A, Sahafi A. Finite element analysis of stresses in endodontically treated, dowel-restored teeth. J Prosthet Dent 2005;94:321-329.
• ¹⁷
  Reinhardt RA, Krejci RF, Pao YC, Stannard JG. Dentin stresses in post-reconstructed teeth with diminishing bone support. J Dent Res 1983;62:1002-1008.
• ¹⁸
  Friedman CM, Sandrik JL, Heuer MA, Rapp GW. Composition and mechanical properties of gutta-percha endodontic points. J Dent Res. 1975;54:921-925.
• ¹⁹
  Lanza A, Aversa R, Rengo S, Apicella D, Apicella A 3D FEA of cemented steel, glass and carbon posts in a maxillary incisor. Dent Mater 2005;21:709-715.
• ²⁰
  Pessoa RS, Vaz LG, Marcantonio E Jr, Vander Sloten J, Duyck J, Jaecques SV. Biomechanical evaluation of platform switching in different implant protocols: computed tomography-based three-dimensional finite element analysis. Int J Oral Maxillofac Implants 2010;25:911-919.
• ²¹
  Rocha EP, Anchieta RB, Freitas A C Jr, de Almeida EO, Cattaneo PM, Chang Ko C. Mechanical behavior of ceramic veneer in zirconia-based restorations: a 3-dimensional finite element analysis using microcomputed tomography data. J Prosthet Dent 2011;105:14-20.
• ²²
  Dejak B, Młotkowski A. Finite element analysis of strength and adhesion of cast posts compared to glass fiber-reinforced composite resin posts in anterior teeth. J Prosthet Dent 2011;105:115-126.
• ²³
  Okamoto K, Ino T, Iwase N, Shimizu E, Suzuki M, Satoh G, et al.. Three-dimensional finite element analysis of stress distribution in composite resin cores with fiber posts of varying diameters. Dent Mater J 2008;27:49-55.
• ²⁴
  Marghalani TY, Hamed MT, Awad MA, Naguib GH, Elragi AF. Three-dimension finite element analysis of custom made ceramic dowel made using CAD/CAM technology. J Prosthodont 2012;21:440-450.

Publication Dates

• Publication in this collection
  Nov-Dec 2015

History

• Received
  04 Aug 2015
• Accepted
  01 Dec 2015