Abstract
Objective The purpose of this study was to evaluate, by using finite element analysis, the influence of post location and occlusal loading location on the stress distribution pattern inside the root canals of the mandibular 1st molar.
Material and Methods Three different 3-D models of the mandibular 1st molar were established: no post (NP) – a model of endodontic and prosthodontic treatments; mesiobuccal post (MP) – a model of endodontic and prosthodontic treatments with a post in the mesiobuccal canal; and distal post (DP) – a model of endodontic and prosthodontic treatments with a post in the distal canal. A vertical force of 300 N, perpendicular to the occlusal plane, was applied to one of five 1 mm2 areas on the occlusal surface; mesial marginal ridge, distal marginal ridge, mesiobuccal cusp, distobuccal cusp, and central fossa. Finite element analysis was used to calculate the equivalent von Mises stresses on each root canal.
Results The DP model showed similar maximum stress values to the NP model, while the MP model showed markedly greater maximum stress values. The post procedure increased stress concentration inside the canals, although this was significantly affected by the site of the force.
Conclusions In the mandibular 1st molar, the distal canal is the better place to insert the post than the mesiobuccal canal. However, if insertion into the mesiobuccal canal is unavoidable, there should be consideration on the occlusal contact, making central fossa and distal marginal ridge the main functioning areas.
Finite element analysis; Molar; Post; Stress
Introduction
The tooth is a complex structure mainly composed of pulp, dentin, cementum, and enamel, which is surrounded by periodontal tissues including the periodontal ligament and alveolar bone19. It deals with the most important function of grinding and chewing of foods, which subsequently results in force generation and transmission from the crown to the alveolar bone. The natural biomechanical balance of the tooth is well suited to this purpose, with the material properties of its aforementioned components functioning together to achieve this end. Therefore, the tooth must withstand the stress exerted during the chewing process; if it does not, it will fracture, especially in the root, resulting in permanent loss of function12.
It is commonly accepted that severely damaged tooth, whether by dental caries or fracture, must undergo the root canal therapy. When an extensive amount of coronal structure is lost, the post would be the first choice of treatment to prevent the core material from being detached6,8,18,26. It is in agreement that the post, if needed in the posterior region, should be inserted in the largest and straightest canal; namely, the distal canal in the mandibular molars, and the palatal canal in the maxillary molars26. Root canal treatment, which is intended to save the tooth, may adversely cause iatrogenic damage, particularly during canal preparation, canal filling, and post preparation2,23,28. It is obvious that endodontic treatment, whether with the post or not, changes the balance of the tooth structure, so the different phenomenon would happen as a force transmission system.
The stress distribution pattern within the root is the most important determinant in the fracture mechanism12. Concentrated stress may initiate and promote crack propagation along weakened surfaces14. However, it is not easy to directly measure the stress exerted on the root. Finite element analysis is a numerical method to investigate the stress and strain exerted on relevant structures, and in turn predict where a fracture might occur30. Many authors used finite element analysis to estimate the effect of posts on the stress distribution patterns, mostly in incisor and premolar teeth1,3,4,15,16,21,24,27. However, few authors studied the molars in the same regard, and none of those focused on the root canals of molars6,17. Therefore, this study was designed to evaluate the influence of post location and loading location on the stress distribution pattern inside root canals of the mandibular 1st molar.
Material and methods
3-D model of the mandibular 1st molar
A 3-D model of the mandibular 1st molar having two roots was created by using Design Modeler™, a part of ANSYS™ v.15.0 software (ANSYS, Inc., Canonsburg, Pennsylvania, U.S.A.). External shapes of mesial and distal roots were generated based on the information obtained from the textbook, followed by the generation of five cusped coronal portion20. The canals, two on mesial and one oval shape on distal, were then fabricated inside the roots with average contour and curvature5,10,11. In the clinical view, the mesiobuccal canal and the mesiolingual canal had the same primary root canal curvature of 28 degrees, while the distal root canal had almost straight morphology showing slight curvature at its apical third in a distal direction5,20. In the proximal view, the primary root canal curvatures of mesiobuccal canal and mesiolingual canal were 21 degrees and 20 degrees each5. The canals were smoothly connected inside the coronal portion to form the pulp chamber in the quadrilateral cervical cross section, which is distally tapered from the wider buccolingual measurement of the mesial of the tooth. The roots were surrounded by 0.2 mm periodontal ligament, which was enclosed by supporting alveolar bone at an intact bone level 3 mm below the cementoenamel junction (CEJ). The tooth dimensions are shown in Figure 1 20.
The model was subsequently subjected to access opening, root canal shaping and filling, mimicking the clinical procedures. All three canals were shaped with .06 taper instruments. After these procedures, the two mesial canals measured 1.3 mm in diameter at canal orifice, while the oval shape of distal canal had 2.0 mm in major axis and 1.3 mm in minor axis at canal orifice. The uppermost portions, at 2.0 mm from the canal orifices, of canal filling materials were then removed26. To insert the post, either the distal or mesiobuccal canal was subjected to drilling of ideal angulation and length. The drilling depths from the CEJ were 7.0 mm for the mesiobuccal canal and 8.0 mm for the distal canal. A post of 1.1 mm in diameter and 11 mm in length was inserted into the prepared bed of the distal or mesiobuccal root canal, with a 0.1 mm layer of cement around the post, and the empty space inside the teeth was filled with core material. To mimic a tooth with severe loss of coronal dentin, most of the coronal component was removed so that only 2.0 mm of the coronal portion above the cementoenamel junction was left; the removed portions were replaced with the restorative material. Finally, a 1.0 mm axial reduction and a 1.5 mm occlusal reduction were conducted, and the crown was added to the prepared coronal portion (Figure 2).
Transparent images of three different models showing internal root canals (A: no post model; B: mesiobuccal post model; C: distal post model). (D-G) Series of procedures to generate the mandibular 1st molar; (G) five loading locations are indicated on the crown
Finite element analysis
The models were transferred to the static structural analysis system in ANSYS for finite element analysis. Zirconium oxide was chosen as the material for the prosthetic crown, and glass fiber for the post. The canal was filled using gutta-percha, and the post was cemented into the canal using Panavia. The core and coronal portion were constructed using resin-based composite material. All materials, other than the glass fiber post, were assumed to be homogeneous, isotropic, and linear elastic; the glass fiber post was considered as orthotropic, linear elastic material. The elastic moduli and Poisson’s ratio of the materials used in this study are shown in Figure 3 4,6,15,25.
The 3-D meshes were generated using 10-node tetrahedral elements, with three degrees of freedom per node. The accuracy of the element size was checked so that same element size of 0.78 mm in average was achieved for each model to eliminate the mesh-dependent variation on the result. Thus, the meshing procedure generated about 200,000 nodes and 130,000 elements. All assemblies were assumed to be fully bonded, and the bottom of the alveolar bone was fixed to prevent any rigid dynamic motions. To assess the effect of tooth angulation, the long axis of the tooth was tilted 9 degrees mesially and 8 degrees lingually from the axis of occlusal loading29. A vertical force of 300 N, perpendicular to the occlusal plane, was applied to one of five 1 mm2 areas on the occlusal surface of the crown: the mesial marginal ridge, distal marginal ridge, mesiobuccal cusp, distobuccal cusp, and central fossa7,9.
Results
The root canals were divided into three portions: the coronal third, middle third, and apical third. Figures 4-6 show the equivalent von Mises stress distributions on the internal surfaces of root canals under five different loading locations. Table 1 shows the estimated maximum equivalent von Mises stress on the internal surface of each root canal at each portion. The largest among the three stress values from each canal is highlighted and plotted in Figures 7-9 to represent the canals in each model.
Equivalent von Mises stress distributions on the internal surfaces of root canals of the no post model under five different loading locations (A: distobuccal cusp; B: mesiobuccal cusp; C: mesial marginal ridge; D: central fossa; E: distal marginal ridge)
Maximum values of Equivalent Von Mises stresses (MPa) on inner surfaces of root canals under five different loading locations (The largest among the three stress values from each canal is is highlighted)
Maximum values of Equivalent Von Mises stresses on the inner surfaces of root canals of the no post model under five different loading locations. MB: mesiobuccal canal, ML: mesiolingual canal, D: distal canal
Equivalent von Mises stress distributions on the internal surfaces of root canals of the distal post model under five different loading locations (A; distobuccal cusp; B: mesiobuccal cusp; C: mesial marginal ridge; D: central fossa; E: distal marginal ridge)
Equivalent von Mises stress distributions on the internal surfaces of root canals of the mesiobuccal post model under five different loading locations (A: distobuccal cusp; B: mesiobuccal cusp; C: mesial marginal ridge; D: central fossa; E: distal marginal ridge)
Maximum values of Equivalent Von Mises stresses on the inner surfaces of root canals of the distal post model under five different loading locations. MB: mesiobuccal canal, ML: mesiolingual canal, D: distal canal
Maximum values of Equivalent Von Mises stresses on the inner surfaces of root canals of the mesiobuccal post model under five different loading locations. MB: mesiobuccal canal, ML: mesiolingual canal, D: distal canal
NP Model
The locations in which the maximum stress values were determined in each canal varied as the area of force was changed. The mesial canals experienced more stress than the distal canal, except when the force was applied to distal marginal ridge. Most forces made the maximum stress occur within the apical third for the distal and mesiolingual canals. On the other hand, the location of maximum stress in the mesiobuccal canal depended on where the force was applied.
DP Model
The two mesial canals showed stress distribution patterns similar to those of the NP model. In the distal canal, force applied to the distal marginal ridge, distobuccal cusp, and central fossa caused maximum stress in the middle third, where the apex of the post was located. However, when force was applied to the two mesial portions, the location of maximum stress in the distal canal was similar to that of the NP model the NP model, even though the stress distribution pattern was different.
MP Model
The mesiolingual and distal canals showed stress distribution patterns similar to those in the NP model. Applying forces to the buccal or mesial portions of the crown significantly increased maximum stress values within the mesiobuccal canal, with the largest value occurring when loaded on the mesiobuccal cusp. Forces on the central fossa and distal marginal ridge changed the stress distribution patterns within the mesiobuccal canal, causing a slight increase in the maximum stress value. Regardless of the loading location, the maximum stress in the mesiobuccal canal occurred in the middle third.
Discussion
The mandibular 1st molar is the first permanent tooth to erupt, emerging at around six years old. For this reason, coronal loss because of dental caries or fracture is more likely in the mandibular 1st molar, and such coronal loss may in turn require root canal treatment followed by post and crown procedure. The morphology of the mandibular 1st molar does not vary as much as that of the mandibular 2nd molar20. Moreover, the external and internal root shapes of the mandibular 1st molar have been investigated and analyzed by many authors5,10,11,20. With such accumulated data, the “standard morphology” of the mandibular 1st molar for finite element analysis was created.
Most authors in this field have designed their models by scanning a real tooth and adding the inner shape of the root canals to the outer shape obtained4,6,15,16,24,27. In this study, however, both outer and inner shapes of the tooth were designed using information accumulated from many previous authors5,10,11,20. In this way, it is possible to use the model to evaluate the influence of different root shapes, i.e., different root lengths or curvatures, on the stress distribution pattern inside the root canals by simply altering the shape of root designed.
Procedures such as root canal preparation and filling, or post drilling, can cause damage, especially micro-cracks, on the inner surface of the root canals2,23,28. Furthermore, concentrated stress may initiate and promote crack propagation along these weakened surfaces14. However, no study has focused on the internal surfaces of root canals of molar teeth before our research. Dejak and Mlotkowski6 (2013) compared the equivalent stresses in molars restored with endocrowns as well as posts and cores. They found that the higher equivalent stresses occurred in molar restored with posts, but they mentioned no equivalent stresses inside the root canals. On the other hand, Lu, et al.17 (2013) analyzed the Von Mises stresses both at the root canal orifice and inside the root canal with different adhesives17. In their study, they reported stress concentration near the root canal and at the bottom of the root canal, irrespective of adhesive materials. However, they did not estimate the influence of the post on the stress distribution pattern inside root canals. Therefore, it was worth to numerically analyze the influence of the post insertion with variant loading locations on the stress distribution pattern and the stress concentration inside the root canals of the molar.
The diameter of the post drill was larger than that of the root canal, which resulted in the generation of internal and external line angles of dentin at the post apex (Figure 2). While the internal line angle could be rounded using a post drill with a round end, it was not possible to properly round off the external line angle; this caused the stress concentration on the middle third of post inserted canals. On the contrary, maximum stress was reduced in the coronal third under most of the five forces; it could therefore be said that post procedure moved stress from the coronal third to the middle third. This result is supported by the previous study of Cailleteau, et al.3 (1992), in which the maximum stresses were shifted to the apical end of the post by the post placement. Reinhardt, et al.22 (1983) also found the stress concentration in the small amount of dentin remaining near the post apex, which agrees with our result22. This result can also be supported by the study of Liu, et al.16 (2014), in which they reported the stress increase at the apex and the decrease in the coronal root.
The mesial and distal root canals display different stress transmission characteristics, mainly because of a difference in morphology and mesiolingual angulation of the tooth. In our study, the mesial canals transferred stress less proficiently and showed high stress concentration. It follows that there is a greater chance of crack generation and propagation in the mesial canals. In the distal canal of the NP model, maximum stress was mostly exerted on the apical third. This could be explained by the influence of the straightness of the distal canal and the decrease of the circumference of root canal as it goes to the apex. Also, the value of maximum stress was significantly affected by the loading locations, due to the forces in different locations generated at different moments on the distal canal. In the DP model, even though the post procedure increased the stress concentration in the middle third of the distal canal, the increments were not high enough to change the maximum stress pattern in the distal canal. This meant that the maximum stress pattern of the distal canal in the DP resembled that in the NP. On the other hand, the location of maximum stress varied in the mesiobuccal canal of the NP model, depending on where the force was applied. This would be explained by the fact that the curved mesiobuccal canal, would be in different moments and force environments with varying locations of the load. In the mesiobuccal canal of the MP model, the post procedure increased the stress concentration in the middle third to a much greater degree, especially when loaded on distobuccal cusp, mesiobuccal cusp, and mesial marginal ridge. This could be partially explained by the different distance from the axis of force to the apex of the post. Among the five loading locations, the force on mesiobuccal cusp, in which the distance from the axis of force to the apex of the post was the closest, made the largest maximum stress increment in the middle third of the mesiobuccal canal. This could be the explanation for the difference in maximum stress pattern on the mesiobuccal canal between the NP and the MP models.
Many authors have used finite element analysis to compare the effect of post material on stress distribution patterns in the roots; most have concluded that a glass fiber post with an elastic modulus similar to dentin results in a well-distributed stress pattern with minimal concentration of stress1,15,16,21,24,27. In our study, the glass fiber post was chosen because it is expected to induce a smaller stress concentration than the other post materials. The influence of post material, as well as post diameter and length, on the stress distribution pattern should be analyzed in future studies.
In this study, all materials, except the glass fiber post, were assumed to be homogenous, isotropic, and linearly elastic; however, this may not match clinical reality. The periodontal ligament is a vascular connective tissue composed principally of collagen fibers and tissue fluids19. This implies that the assumptions of homogeneity, isotropicity, and linear elasticity must have been erroneous. However, for simplicity, and utility on finite element analysis, many authors have capitulated to the assumption and focused on finding the most suitable elastic modulus values of periodontal ligament as a homogenous, isotropic, and linear elastic material3,6,16,24,25. However, recent authors have regarded the periodontal ligament as a non-linear, hyperelastic material, and have found variable elastic modulus values for periodontal ligament13. The influence of periodontal ligament as a non-linear hyperelastic material should be considered in further studies.
In this study, because of the curvature of the mesiobuccal canal, the post could not be positioned as deep as the half of the bone insertion. Instead, the mesiobuccal post bed was prepared to the half of the root, which was the 7 mm from the CEJ, so that as little dentin would be removed as possible. Similarly, the distal post was inserted into the 8 mm from the CEJ to get rid of the possible influence of post length on the results. Moreover, only the vertical force on five areas of the occlusal surface was considered. However, chewing food can generate force anywhere on the crown, and in any direction. Besides, the young are susceptible to dental caries, while periodontal disease and sequential bone loss are more common in the aged. Therefore, the influences of the quantitative surrounding bone level as well as the force in other directions, especially horizontal, should be further investigated.
Conclusion
Within the limitations of our study, we concluded the following for the mandibular 1st molar:
-
Stress distribution pattern was significantly influenced by force location.
-
The post procedure increased the stress concentration in the post-placed root canal, especially in the middle third.
-
The distal canal is a better place to insert the post than the mesiobuccal canal.
-
If insertion into the mesiobuccal canal is unavoidable, then the central fossa and distal marginal ridge should be the main functioning areas, avoiding buccal and mesial aspects on the crown from being loaded as possible.
Acknowledgements
Hong Gi Yoon and Hyun Keun Oh contributed equally to this research as first authors.
The authors deny any conflicts of interest related to this study.
References
- 1 - Al-Omiri MK, Rayyan MR, Abu-Hammad O. Stress analysis of endodontically treated teeth restored with post-retained crowns: a finite element analysis study. J Am Dent Assoc. 2011;142:289-300.
- 2 - Barreto MS, Moraes RA, Rosa RA, Moreira CH, Só MV, Bier CA. Vertical root fractures and dentin defects: effects of root canal preparation, filling, and mechanical cycling. J Endod. 2012;38:1135-9.
- 3 - Cailleteau JG, Rieger MR, Akin JE. A comparison of intracanal stresses in a post-restored tooth utilizing the finite element method. J Endod. 1992;18:540-4.
- 4 - Chen D, Wang N, Gao Y, Shao L, Deng B. A 3-dimensional finite element analysis of the restoration of the maxillary canine with a complex zirconia post system. J Prosthet Dent. 2014;112:1406-15.
- 5 - Cunningham CJ, Senia ES. A three-dimensional study of canal curvatures in the mesial roots of mandibular molars. J Endod. 1992;18:294-300.
- 6 - Dejak B, Mlotkowski A. 3D-Finite element analysis of molars restored with endocrowns and posts during masticatory simulation. Dent Mater. 2013;29:e309-17.
- 7 - Dinçer M, Isik Aslan B. Effects of thermoplastic retainers on occlusal contacts. Eur J Orthod. 2010;32:6-10.
- 8 - Fernandes AS, Shetty S, Coutinho I. Factors determining post selection: a literature review. J Prosthet Dent. 2003;90:556-62.
- 9 - Gibbs CH, Mahan PE, Lundeen HC, Brehnan K, Walsh EK, Sinkewiz SL, et al. Occlusal forces during chewing - influences of biting strength and food consistency. J Prosthet Dent. 1981;46:561-7.
- 10 - Gu Y, Lu Q, Wang H, Ding Y, Wang P, Ni L. Root canal morphology of permanent three-rooted mandibular first molars - part I: pulp floor and root canal system. J Endod. 2010;36:990-4.
- 11 - Gu Y, Lu Q, Wang P, Ni L. Root canal morphology of permanent three-rooted mandibular first molars: Part II - measurement of root canal curvatures. J Endod. 2010;36:1341-6.
- 12 - Hayashi M, Takahashi Y, Imazato S, Ebisu S. Fracture resistance of pulpless teeth restored with post-cores and crowns. Dent Mater. 2006;22:477-85.
- 13 Huang H, Tang W, Yan B, Wu B, Cao D. Mechanical responses of the periodontal ligament based on an exponential hyperelastic model: a combined experimental and finite element method. Comput Methods Biomech Biomed Engin. 2016;19:188-98.
- 14 - Kishen A, Kumar GV, Chen NN. Stress-strain response in human dentine: rethinking fracture predilection in postcore restored teeth. Dent Traumatol. 2004;20:90-100.
- 15 - 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-15.
- 16 - Liu S, Liu Y, Xu J, Rong Q, Pan S. Influence of occlusal contact and cusp inclination on the biomechanical character of a maxillary premolar: a finite element analysis. J Prosthet Dent. 2014;112:1238-45.
- 17 - Lu LW, Meng GW, Liu ZH. Finite element analysis of multi-piece post-crown restoration using different types of adhesives. Int J Oral Sci. 2013;5:162-6.
- 18 - Morgano SM. Restoration of pulpless teeth: application of traditional principles in present and future contexts. J Prosthet Dent. 1996;75:375-80.
- 19 - Natali AN, Pavan PG, Scarpa C. Numerical analysis of tooth mobility: formulation of a non-linear constitutive law for the periodontal ligament. Dent Mater. 2004;20:623-9.
- 20 - Nelson SJ, Ash MM Jr. Mandibular first molar. In: Nelson SJ. Wheeler’s dental anatomy, physiology, and occlusion. 9th ed. St. Louis: Saunders Elsevier; 2010. p. 189-99.
- 21 - Pegoretti A, Fambri L, Zappini G, Bianchetti M. Finite element analysis of a glass fibre reinforced composite endodontic post. Biomaterials 2002;23:2667-82.
- 22 - Reinhardt RA, Krejci RF, Pao YC, Stannard JG. Dentin stresses in post-reconstructed teeth with diminishing bone support. J Dent Res. 1983;62:1002-8.
- 23 - Rippe MP, Santini MF, Bier CA, Baldissara P, Valandro LF. Effect of root canal preparation, type of endodontic post and mechanical cycling on root fracture strength. J Appl Oral Sci. 2014;22:165-73.
- 24 - Roscoe MG, Noritomi PY, Novais VR, Soares CJ. Influence of alveolar bone loss, post type, and ferrule presence on the biomechanical behavior of endodontically treated maxillary canines: strain measurement and stress distribution. J Prosthet Dent. 2013;110:116-26.
- 25 - Ruse ND. Propagation of erroneous data for the modulus of elasticity of periodontal ligament and gutta percha in FEM/FEA papers: a story of broken links. Dent Mater. 2008;24:1717-9.
- 26 - Schwartz RS, Robbins JW. Post placement and restoration of endodontically treated teeth: a literature review. J Endod. 2004;30:289-301.
- 27 - Sorrentino R, Aversa R, Ferro V, Auriemma T, Zarone F, Ferrari M, et al. Three-dimensional finite element analysis of strain and stress distributions in endodontically treated maxillary central incisors restored with different post, core and crown materials. Dent Mater. 2007;23:983-93.
- 28 - Tamse A. Iatrogenic vertical root fractures in endodontically treated teeth. Endod Dent Traumatol. 1988;4:190-6.
- 29 - Tong H, Kwon D, Shi J, Sakai N, Enciso R, Sameshima GT. Mesiodistal angulation and faciolingual inclination of each whole tooth in 3-dimensional space in patients with near-normal occlusion. Am J Orthod Dentofacial Orthop. 2012;141:604-17.
- 30 - Wandscher VF, Bergoli CD, Oliveira AF, Kaizer OB, Souto Borges AL, Limberguer IF. Fatigue surviving, fracture resistance, shear stress and finite element analysis of glass fiber posts with different diameters. J Mech Behav Biomed Mater. 2015;43:69-77.
Publication Dates
-
Publication in this collection
2018
History
-
Received
2 Sept 2016 -
Reviewed
12 Nov 2016 -
Accepted
22 Dec 2016