Next Article in Journal
An Agent-Specific Stochastic Model of Generalized Reaching Task Difficulty
Next Article in Special Issue
A Microcomputed Tomographic Evaluation of Root Canal Morphology of Maxillary Second Premolars in a Pakistani Cohort
Previous Article in Journal
Community Composition and Ex Situ Cultivation of Fungi Associated with UNESCO Heritage Monuments in the Bay of Naples
Previous Article in Special Issue
Functional Evaluation of a Novel Multi-Axial Alveolar Distractor—Preliminary In Vivo Animal Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computerized Generation and Finite Element Stress Analysis of Endodontic Rotary Files

by
Victor Roda-Casanova
1,*,
Álvaro Zubizarreta-Macho
2,3,
Francisco Sanchez-Marin
1,
Óscar Alonso Ezpeleta
4,
Alberto Albaladejo Martínez
3 and
Agustín Galparsoro Catalán
2
1
Mechanical Engineering and Construction Department, Universitat Jaume I, 12071 Valencia, Spain
2
Department of Dentistry, Alfonso X el Sabio University, 28691 Madrid, Spain
3
Department of Orthodontics, University of Salamanca, 37008 Salamanca, Spain
4
Department of Endodontics, University of Zaragoza, 50009 Zaragoza, Spain
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(10), 4329; https://doi.org/10.3390/app11104329
Submission received: 20 April 2021 / Accepted: 8 May 2021 / Published: 11 May 2021
(This article belongs to the Special Issue Applied Simulation and Experiment Research in Dentistry)

Abstract

:
Introduction: The finite element method has been extensively used to analyze the mechanical behavior of endodontic rotary files under bending and torsional conditions. This methodology requires elevated computer-aided design skills to reproduce the geometry of the endodontic file, and also mathematical knowledge to perform the finite element analysis. In this study, an automated procedure is proposed for the computerized generation and finite element analysis of endodontic rotary files under bending and torsional conditions. Methods: An endodontic rotary file with a 25 mm total length, 0.25 mm at the tip, 1.20 mm at 16 mm from the tip, 2 mm pitch and squared cross section was generated using the proposed procedure and submitted for analysis under bending and torsional conditions by clamping the last 3 mm of the endodontic rotary file and applying a transverse load of 0.1 N and a torsional moment of 0.3 N · cm . Results: The results of the finite element analyses showed a maximum von Mises stress of 398 MPa resulting from the bending analysis and a maximum von Mises stress of 843 MPa resulting from the torsional analysis, both of which are next to the encastre point. Conclusions: The automated procedure allows an accurate description of the geometry of the endodontic file to be obtained based on its design parameters as well as a finite element model of the endodontic file from the previously generated geometry.

1. Introduction

Endodontic rotary files have experienced continuous development, both in the chemical composition of alloys and geometrical design, which has improved their mechanical behavior during operation [1]. However, unexpected fractures of endodontic rotary files within the root canal system are still a concern and remain a challenge for clinicians despite continuous enhancements in their geometrical design and manufacture processes to reduce the incidence of fractures [2], which ranges from 0.09 % to 5 % [3,4].
Previous studies have suggested that flexural failure (bending) and torsional failure are the main causes for endodontic rotary file fractures [5,6,7]. Torsional failure is caused by the blockage of the endodontic files during rotational movement [5]; however, endodontic rotary instrument failures are mainly caused by cyclic bending fatigue, which occurs when a endodontic rotary instrument rotates in a curved root canal [8]. In addition, endodontic rotary files can be simultaneously subjected to cumulative torsional and bending stress during root canal shaping, which can lead to the fracture of endodontic rotary files and therefore affect the prognosis of root canal treatment as the fractured fragment blocks the access to the apex, preventing root canal system disinfection [9,10].
In addition, root canal dentine removal during the retrieval procedures of separated files can weaken the tooth structure up to 4.25 ± 0.68 mm 3 [11]. Moreover, fractures of endodontic instruments on teeth with periapical pathology significantly decrease the success rate of root canal treatment [12]. Therefore, improvements of the chemical composition of the alloy and geometrical design may increase torsional and bending resistance, preventing the fracture of the endodontic rotary instruments [13]. Furthermore, some geometrical factors have been highlighted as influencing the instrument’s performance, such as taper [14], apical diameter [7,14], cross section [15], flute length, helix angle and pitch [16]; however, most studies were conducted experimentally using custom-made cyclic fatigue devices, hardly reproducing the clinical setting and the absence of regulation prevents data comparison [17].
However, finite element (FE) analysis is an alternative that allows the mechanical behavior and stress distribution of endodontic rotary instruments subjected to different conditions to be simulated through mathematical analysis [18]. Moreover, FE analyses have been previously used in endodontics to assess both the mechanical behavior of endodontic instruments [19] and stress distribution during root canal treatments [16].
Two different approaches can be found in recent works when developing FE models for the stress analysis of endodontic files, which differ in the way in which the endodontic file geometry is obtained. In both cases, the obtained geometry must be converted to a finite element model by means of a FE analysis pre-processor, in order to proceed to its analysis. In addition, elevated solid mechanics knowledge is required in both cases, in order to define a finite element model that can accurately reproduce the physics of the problem.
In the most conventional approach [1,20,21], the generation of the geometry of the file is conducted through computer-aided design (CAD) software. This approach requires elevated CAD knowledge to accurately reproduce the geometric characterization of the endodontic rotary file parameters. In fact, there are several cases in which the generated geometries cannot accurately reproduce the geometry of the actual endodontic file, especially at the transition between the active part and the shaft.
In another approach [22], the geometry of the endodontic file is obtained from a tomography scanner. This approach requires expensive equipment in order to obtain the computer representation from the tomography scanner and elevated knowledge to manipulate the obtained geometry files.
Both processes are costly in terms of time and have to be accomplished for each assigned case of design of various geometries and cases of investigation. Therefore, generating an automated procedure that simplifies the procedures described above and enables the performance of FE analyses without the need for this knowledge will help us to better understand the mechanical behavior of endodontic rotary files under various conditions. Thus, the aims of this work are as follows:
  • To develop an automated procedure to obtain an accurate description of the geometry of the endodontic file from its design parameters;
  • To develop an automated procedure to obtain an FE model of the endodontic file from the generated geometry;
  • To demonstrate the capabilities of the proposed procedures.

2. Materials and Methods

This section describes a new method developed for the computerized generation of endodontic files and the obtention of FE models for their stress analysis. The geometry considered for the endodontic file is shown in Figure 1a. It is defined with respect to a coordinate system S f , whose z f axis is alined with the theoretical axis of rotation of the endodontic file and its origin O f is at the tip of the endodontic file. The geometry of the endodontic file is fully parametrized by the following set of parameters:
  • The diameter of the shaft, d s h ;
  • The diameter of the tip of the active part, d a p ;
  • The length of the active part, L a p ;
  • The total length of the file, L t o t a l ;
  • The pitch of the active part, p z , which represents the distance between adjacent points in the cutting edge.
Three different parts can be identified for this endodontic file: the active part, the shaft and the transition part. The active part contains the cutting edges and has an arbitrary cross section (squared, triangular, etc.). The shaft has a circular cross section and connects the active part with the support of the endodontic file (which is not included in this model). Finally, the transition part allows for a smooth transition from the cross section of the active part to the cross section of the shaft. Figure 2 shows an example of the cross section of the endodontic file at each one of these parts.
The computerized generation of this endodontic file is achieved as the intersection of two different geometries: the raw material (Figure 1b) and the machined material (Figure 1c). The raw material corresponds to a geometry of the endodontic file before the cutting edges are machined over it, and its generation is described in Section 2.1. The geometry of the machined material is an auxiliary geometry that is used to generate the cutting edges of the endodontic file, and its generation is described in Section 2.2. The intersection of both geometries and the obtention of the final geometry of the endodontic file is described in Section 2.3.
Having a computational description of the geometry of the endodontic file allows for an automated generation of FE models for its stress analysis. Section 2.4 provides further details on how to obtain these FE models.

2.1. Definition of the Geometry of the Raw Material for the Endodontic File

The raw material (Figure 1b) has an axisymmetric geometry that consists of two different parts:
  • A cylindrical part, which corresponds to the shaft, characterized by its length ( L t o t a l L a p ) and its diameter ( d s h );
  • A truncated cone, which corresponds to the active part. It is characterized by its length ( L a ) and its diameters at the front ( d a p ) and back ends ( d s h ).
The position of a point P on the surface of the raw material is given by the parametric equation r f ( θ , v ) , which provides the coordinates of point P in coordinate system S f from its parametric coordinates θ 0 , 360 and v 0 , L t o t a l :
r f ( θ , v ) = λ ( v ) 2 · cos θ + φ ( v ) λ ( v ) 2 · sin θ + φ ( v ) v T
where angle φ ( v ) is defined as a function of v; assuming a constant pitch p z along the endondontic file, it can be defined as
φ ( v ) = v p z · 2 π rad
On the other hand, λ ( v ) represents the diameter of the raw material at coordinate v, and it is defined as
λ ( v ) = d a p + v · d s h d a p L a p , if 0 v < L a p d s h , if L a p v < L t o t a l

2.2. Definition of the Geometry of the Machined Material for the Endodontic File

The machined material (Figure 1c) is an auxiliary geometry that is used for the generation of the geometry of the endodontic file, and it is obtained from the transformation (rotation, translation and scalation) of an arbitrary cross section.
For the definition of this geometry, let us consider an auxiliary coordinate system S g , whose z g axis is parallel to the z f axis of coordinate system S f , as indicated in Figure 3 left. The origin of coordinate system S g is located on axis z f at a distance v from O f , and S g is rotated with respect to S f an angle φ ( v ) .
The cross section of the machined material is defined with respect to coordinate system S g , and it is circumscribed into a normalized circumference (d = 1 mm), as illustrated in Figure 3 right for the particular case of a squared cross section. The position of any point P over the contour of the file section is defined by a parametric function r g ( θ ) that provides, for a given parametric coordinate θ 0 , 360 , the homogeneous coordinates [23] of point P in coordinate system S g :
r g ( θ ) = x g ( θ ) y g ( θ ) z g ( θ ) 1 T
where z g ( θ ) = 0 , since the cross section is contained in the x g y g plane. In the particular case of squared transversal cross section, x g ( θ ) and y g ( θ ) can be defined as
x g ( θ ) = d · 2 4 · tan ( 90 θ ) , if 45 θ < 135 1 , if 135 θ < 225 tan ( 270 θ ) , if 225 θ < 315 1 , otherwise
y g ( θ ) = d · 2 4 · 1 , if 45 θ < 135 tan ( 180 θ ) , if 135 θ < 225 1 , if 225 θ < 315 tan θ , otherwise
where d = 1 mm is the diameter of the normalized circumference. The definitions of x g ( θ ) and y g ( θ ) for other transversal cross sections can be trivially derived.
The position of a point P on the surface of the machined material is given by the parametric equation r f ( θ , v ) , which provides the homogeneous coordinates of point P in coordinate system S f from its parametric coordinates θ 0 , 360 and v 0 , L t o t a l :
r f ( θ , v ) = M f g ( v ) · r g ( θ )
where r g ( θ ) is given by Equation (4), v is a distance coordinate (see Figure 3) and M f g ( v ) is a transformation matrix defined as
M f g ( v ) = λ ( v ) · cos φ ( v ) λ ( v ) · sin φ ( v ) 0 0 λ ( v ) · sin φ ( v ) λ ( v ) · cos φ ( v ) 0 0 0 0 1 v 0 0 0 1
The advantage of using homogeneous coordinates is that the transformation matrix M f g ( v ) allows three simultaneous transformations of the transversal cross section of the machined material to be performed:
  • The translation of the cross section a distance v along the z f axis;
  • The rotation of the cross section at an angle φ ( v ) around the z f axis;
  • The scaling of the cross section in the x f y f plane by a factor λ ( v ) .
The angle φ ( v ) is defined as a function of v and, assuming a constant pitch p z along the endondontic file, this angle can be defined as
φ ( v ) = v p z · 2 π rad
Note that the equality φ ( v ) = φ ( v ) must be satisfied to guarantee the robustness of the proposed formulation. On the other hand, the scale factor λ ( v ) is defined as
λ ( v ) = d a p + v · d s h d a p L a p , if 0 v < L a p d a p + v · d s h d a p L a p + a s h · v L a p 2 , if L a p v < L t o t a l
In the interval corresponding to the active part, it is fulfilled that λ ( v ) = λ ( v ) . In the transition part and the shaft, λ ( v ) > λ ( v ) , because a quadratic term is included in the definition of λ ( v ) . This quadratic term allows for a smooth transition from the machined material to the raw material, and it is characterized by the parabola coefficient a s h .

2.3. Definition of the Geometry of the Endodontic File

The final geometry of the endodontic file is obtained from the intersection of the geometries of the raw material and machined material, as described above. In order to compute this intersection, it is necessary to compute the inner point of these geometries for the same θ and v coordinates. For such a purpose, the distance from any point in the cross section to the z f axis needs to be calculated. Distances d ( v ) and d ( θ , v ) are defined for the raw material and the machined material (respectively) as
d ( v ) = λ ( v ) 2
d ( θ , v ) = r f ( θ , v ) r f ( θ , v ) · k f
where k f is the unit vector in the direction given by the z f axis of the coordinate system S f .
Note that d ( v ) does not depend on θ because the raw material has a circular cross section and, in consequence, all the points within the section are at the same distance of z f . According to the definition of the scaling factors for the geometries of the raw material and the machine material, the following statements are fulfilled (see Figure 2):
  • In the active part of the endodontic file, d ( θ , v ) d ( v ) ;
  • In the shaft of the endodontic file, d ( θ , v ) > d ( v ) ;
  • In the transition part of the endodontic file, the relation between d ( θ , v ) and d ( v ) cannot be known in advance.
Taking this into account, the geometry of the endodontic file is represented in coordinate system S f by function q f ( θ , v ) as
q f ( θ , v ) = r f ( θ , v ) , if d ( θ , v ) < d ( v ) r f ( θ , v ) , otherwise

2.4. Development of an Fe Model for Stress Analysis of Endodontic Files

This section describes a new procedure for the automated development of FE models for the stress analysis of endodontic files. This procedure covers the typical steps in the generation of an FE model for stress analysis, which include the discretization of the volume of the endodontic file into finite elements, the definition of the material properties for the finite elements of the model and the definition of the loading and boundary conditions of the model.
In the first step, the volume of the endodontic file is discretized into tetrahedral elements by means of a constrained Delaunay tetrahedralization [24]. The constrained Delaunay tetrahedralization requires the boundary surfaces of the endodontic file to be represented by a set of triangles, which will constitute the external faces of the tetrahedrons once the mesh is generated. This boundary representation is achieved following these steps:
  • The external surfaces that define the volume of the endodontic file are represented analytically considering the parametric equation q f ( θ , v ) (Equation (11)), which provides the Cartesian coordinates of a point in the surface of the endodontic file from its parametric coordinates θ 0 , 360 and v 0 , L t o t a l ;
  • The analytical determination of the position of the nodes in the surface of the endodontic file is performed, taking into account the number of desired elements in the longitudinal ( n v ) and angular directions ( n θ ). Let us reiterate that all nodes of the FE mesh are determined analytically using Equation (11);
  • Discretization of the surface of the endodontic file into triangular faces is accomplished using the nodes determined previously.
Using this boundary representation of the endodontic file, the tetrahedral mesh is constructed using the constrained Delaunay tetrahedralization implemented in the TetGen library [25]. This tetrahedralization does not change the position of the nodes defined for the endodontic file surface, and the resulting tetrahedrons keep the faces provided as input information.
When the FE mesh is generated, mechanical properties are assigned to its elements. The mechanical response of a material in its elastic range can be defined by the elastic modulus E and the Poisson’s coefficient ν . Assuming the endondontic file is manufactured using a NiTi alloy, these material parameters are established as E = 36 GPa and ν = 0.3 [1]. Of course, the user can specify other values.
In the last step, the loading and boundary conditions for the FE model are established. Endodontic files are usually tested for bending and torsional loads, and the typical procedure for these tests is described in [26] and summarized in Figure 4. In these analyses, the endodontic file is approached to a cantilever beam. A fixed support is defined at a distance d e n c from the tip of the endodontic file, in such a way that all the nodes below that dimension have all their degrees of freedom constrained. Typical magnitudes for d e n c are between 3 mm [2,20,22,27,28] and 4 mm [29], although other values could be considered [1,21].
Loads are specified at the other part of the endodontic file. For such a purpose, a rigid surface is defined comprising all the nodes within the end face. A reference node is defined at the center of this face, whose movements are coupled to the movements of the rigid surface by means of a multi-point constraint. Loads are applied directly over this reference node. Two different load cases are derived to study the bending (Figure 4a) and torsional (Figure 4b) response of the endodontic file. In the first of them (Figure 4a), a tangential load F t is introduced over the reference node; in the second (Figure 4b), a torsional moment M z is considered over this reference node.
This procedure provides an input file for the FE solver that contains the entire definition of the FE model. It must be emphasized that this implementation allows for the development of fully parametrized FE models, in which the geometry, the mesh density, the material parameters and the loading conditions can be specified by the user without using any other third-party software.

3. Results

The developed procedure has been implemented in Matlab R2018 [30] (MathWorks Inc., Natick, MA, USA). A graphical interface, which is shown in Figure 5, has been programmed to ease the specification of the endodontic file geometry and the configuration of the FE model. After completing all the parameters that define the model, this software provides an input file for the FE solver, which contains the entire definition of the FE model.
The performance of the proposed procedure and its computational implementation is illustrated with the generation and stress analysis of an endodontic file, which is defined by the set of parameters shown in Table 1.
Two different load cases are considered in this study, according to Bonessio et al. [27] and De Arruda Santos et al. [28]. In the first load case, the bending behavior of the endodontic file is studied by considering a transversal force F t = 0.1 N . In the second load case, the torsional behavior of the endodontic file is studied by considering a torsional moment M z = 0.3 N · cm . In both cases, the length of the encastre is d e n c = 3 mm .

3.1. Mesh Convergence Study

A mesh convergence study has been conducted to establish the optimum size and type for the elements of the mesh. In this study, several finite element models are obtained for the endodontic file, in which the density of the mesh is progressively increased by varying the magnitudes of n v and n θ . Four different average element sizes have been considered in the analysis, which vary from 0.05 mm to 0.20 mm , and both linear and quadratic tetrahedral elements have been taken into account. Figure 6 shows one of the meshes used in this study, in which the endodontic file is meshed using quadratic tetrahedrons with an average element size of 0.10 mm ( n θ = 24 and n v = 250 ). The resulting FE models were analyzed using ABAQUS [31] (Dassault Systèmes Simulia Corp, Providence, RI, USA).
As suggested by Żmudzki [32,33], the convergence of the mesh can be assessed by means of an error energy norm E. In this error estimator, the stress vector σ ^ and the averaged stress vector σ * are used to determine a stress error vector Δ σ :
Δ σ = σ * σ ^
Then, for a given element i of the mesh and following the ideas presented by Zienkiewicz [34], an element energy error e i is determined as:
e i = 1 2 · V i Δ σ T · D 1 · Δ σ · d V i
where D 1 is the material elasticity matrix and V i denotes the volume of element i. The energy error e is estimated for the whole finite element model as the sum of the error contributions of all the elements in the mesh:
e = i = 1 N e e i
where N e denotes the number of elements in the mesh. Finally, and considering the total strain energy of the finite element model U, the error energy norm E is calculated as:
E = e U + e
Figure 7a shows the evolution of the error energy norm E as a function of the average element size, determined for each load case (torsion and bending) and element type (linear and quadratic). In all the cases, it can be observed that reducing the average element size also decreases the error in energy norm. For any load case and average element size, the error energy norm is significantly smaller when the problem is solved using quadratic elements than when it is solved using linear elements.
Figure 7b shows the maximum element energy error ( max ( e i ) ) determined for each load case (torsion and bending) and element type (linear and quadratic), as a function of the average element size. The trends observed in this figure are similar to the ones observed in Figure 7a. In general, the maximum element energy error is reduced as the average element size is decreased, and quadratic elements demonstrate a better performance than linear elements. In all the cases the element with the maximum element energy error is found in the vicinity of the encastre, indicating that this is the region of the finite element model where mesh refinement is more needed.
Although there is no universal agreement over which is the error in energy norm that can be used as a threshold to validate the convergence of a given finite element mesh, an error in energy norm of 15 % has been used in some works [35] for such a purpose. Taking this into consideration, and examining the results shown in Figure 7a, it could be said that in this case the mesh convergence is achieved when the finite element model is built using quadratic elements and the average element size is equal or less than 0.10 mm .

3.2. Interpretation of the Results of the Stress Analysis

Figure 8 shows a von Mises stress plot along the endodontic file, which is obtained from a torsional analysis in which the finite element model is built with quadratic elements with an average element size of 0.05 mm . The maximum von Mises stress reached in this analysis is 843 MPa , and it is produced in the vicinity of the encastre (but not in the encastre itself, as it can be observed in the detail of Figure 8). The figure also shows the stress distribution along the cross section where maximum von Mises stress takes place. In both cases, the stress results follow the expected distribution for this type of analyses. The rotation at the free end of the file 0.39 rad .
Figure 9 shows a von Mises stress plot along the endodontic file, which is obtained from a bending analysis in which the finite element model is built with quadratic elements with an average element size of 0.05 mm . The maximum von Mises stress reached in this analysis is 613 MPa , and it is produced in the vicinity of the encastre. The figure also shows the stress distribution along the cross section where the maximum stress is found. In both cases, the stress results follow the expected distribution for this type of analyses. The deflection at the free end of the file is 3.21 mm .
Unlike in the torsional analysis, in the bending analysis the maximum von Mises stress occurs in a node where a boundary condition is applied and, as a consequence, a numerical singularity is produced that leads to unrealistic stress values. This issue is illustrated in Figure 10a, where a detail of the von Mises stress plot over the area where the bending stresses reach their maximum value is shown. This figure also displays the von Misses stress values at the nodes of the edge where the maximum stress takes place. It can be observed that the maximum value of 613 MPa , produced at the node in the corner of the encastre section, is unreasonably high considering the stress magnitude at the surrounding nodes.
According to many researchers [36,37], these stress results at singularity points cannot be considered to evaluate the strength of the endodontic file. To get around this issue, Żmudzki [32] proposes to exclude the stress results at these points and then extrapolate the extreme value from the stress values in the remaining nodes. Following this idea, Figure 10b shows a plot of the evolution of the von Mises stress along the observed edge. A polynomial approximation to the stress values has been obtained, discarding the maximum value of 613 MPa . Following this approximation, the predicted von Misses stress value at the encastre is 398 MPa , leaving the maximum von Mises stress of the model at 474 MPa .

3.3. Parametric Study

The developed procedure allows for a simple and fast generation of finite element models of endodontic files, and, in consequence, it is especially suitable for performing parametric and optimization studies. In this section, the performance of the proposed procedure is illustrated with three parametric studies, in which the pitch of the active part, the length of the active part and the diameter of the tip were varied. In each parametric study, only one parameter was varied and the rest were kept constant and equal to those shown in Table 1.
Figure 11 shows the results obtained from the parametric study in which the pitch of the active part was varied. Four different magnitudes were considered for the pitch of the active part, which were between 1 mm and 4 mm . In each case, the endodontic file has been conveniently rotated so the orientation of the cross section at the encastre coincides with that shown in Figure 9. Figure 11a shows the results for the bending analysis. It can be observed that both the maximum von Mises stress and tip deflection reduced as the pitch of the active part increased. The maximum von Mises stress decreased from 517 MPa (when p z = 1 mm ) to 451 MPa (when p z = 4 mm ). The tip deflection decreased from 3.38 mm (when p z = 1 mm ) to 3.17 mm (when p z = 4 mm ).
On the other hand, Figure 11b shows the results for the torsional analysis. In this case, minor variations are observed for the maximum von Mises stress and the tip rotation, indicating that these parameters may be independent of the pitch of the active part (at least in the observed range).
Figure 12 shows the results obtained from the parametric study in which the length of the active part was varied. Five different magnitudes were considered for the length of the active part, which were between 14 mm and 18 mm . Figure 12a shows the results for the bending analysis. It can be observed that both the maximum von Mises stress and tip deflection increased as the length of the active part increased. The maximum von Mises stress increased from 399 MPa (when L a p = 14 mm ) to 554 MPa (when L a p = 18 mm ). The tip deflection increased from 2.44 mm (when L a p = 14 mm ) to 4.05 mm (when L a p = 18 mm ).
On the other hand, Figure 12b shows the results for the torsional analysis. As happened in the bending analysis, both the maximum von Mises stress and tip rotation increased as the length of the active part increased. The maximum von Mises stress increased from 690 MPa (when L a p = 14 mm ) to 960 MPa (when L a p = 18 mm ). The tip rotation increased from 0.30 rad (when L a p = 14 mm ) to 0.52 rad (when L a p = 18 mm ).
Finally, Figure 13 shows the results obtained from the parametric study in which the diameter of the tip of the active part was varied. Four different magnitudes were considered for the diameter of the tip of the active part, which were between 0.20 mm and 0.35 mm . Figure 13a shows the results for the bending analysis. It can be observed that both the maximum von Mises stress and tip deflection decreased as the diameter of the tip of the active part increased. The maximum von Mises stress decreased from 630 MPa (when d a p = 0.20 mm ) to 295 MPa (when d a p = 0.35 mm ). The tip deflection decreased from 4.20 mm (when d a p = 0.20 mm ) to 2.01 mm (when d a p = 0.35 mm ).
On the other hand, Figure 13b shows the results for the torsional analysis. As happened in the bending analysis, both the maximum von Mises stress and tip rotation decreased as the diameter of the tip of the active part increased. The maximum von Mises stress decreased from 1130 MPa (when d a p = 0.20 mm ) to 500 MPa (when d a p = 0.35 mm ). The tip rotation decreased from 0.51 mm (when d a p = 0.20 mm ) to 0.26 mm (when d a p = 0.35 mm ).

4. Discussion

FE models have been proposed to analyze the mechanical performance of endodontic rotary files, isolating the variables independently under controlled settings [38] to facilitate the comparison between the mechanical behavior and stress distribution of existing endodontic rotary files and even to allow the design of new endodontic rotary instruments by reducing time and costs [27]. In addition, He and Ni used FE methods to improve the geometrical design of existing endodontic rotary instruments and therefore the mechanical behavior under bending and torsional loads to reduce the prevalence of failure [2].
He and Ni also analyzed the influence of the helix angle, taper and flute length on the bending flexibility and torsional stiffness of endodontic files (V-Taper) through numerical simulations and concluded that the geometric features influence on the mechanical performance of endodontic files under bending and torsional conditions [2]. Arbab-Chirani et al. compared numerically the mechanical behavior of Hero, HeroShaper, ProFile, Wtwo and ProTaper F1 endodontic rotary NiTi instruments under bending and torsional conditions and showed that the different designs for tapers, pitch and cutting blades influence on the bending and torsional mechanical behavior [39]. Baek et al. evaluated the effect from pitch and cross-section design on torsional stiffness of NiTi endodontic rotary instruments and evidenced that torsional deformation and fracture of NiTi rotary files might be reduced by reducing the pitch and increasing the cross-sectional areas rather than the center core area [40].
Although metallurgical properties and metallurgical treatments have been shown to influence the properties of NiTi endodontic rotary instruments, geometrical design has been highlighted to be directly related to the stiffness [41], clinical efficiency and cutting performance of endodontic rotary files [29,42]. The cross-section design of endodontic rotary files has experienced a continuous development; Versluis et al. reported that increasing the cross-sectional surface area lead to a decrease in the flexibility of endodontic rotary files by a stress accumulation at the upper and lower edges of the cross section during bending tests, compared with smaller cross-sectional surface areas [17].
In addition, cross-sectional morphology has been also widely analyzed, and Galal and Hamdy reported that a convex triangular cross-section design showed maximum von Misses stresses upon bending tests ( 1962 MPa ), followed by a triangular cross-section design ( 1535 MPa ), parallelogram cross-section design ( 870 MPa ) and rectangular cross-section design ( 814 MPa ); however, a convex triangular cross-section design showed maximum von Misses stresses upon torsional tests ( 4321 MPa ), followed by a triangular cross-section design ( 4285 MPa ), rectangular cross-section design ( 2659 MPa ) and parallelogram cross-section design ( 2307 MPa ).
Additionally, the stress distribution pattern was localized at the cutting edges under bending tests and more concentrated in the base of the flutes in the first 2 mm of the triangular cross-section design under torsional tests [43]. These results agree with the findings reported by Tsao et al., who showed a maximum von Misses stress of 1170 MPa at the cutting edges of a triangular cross-section design endodontic rotary file after applying a transverse load of 1 N at the tip of the endodontic rotary file during bending tests [44]; additionally, Prados-Privado et al. showed a maximum von Misses stress of 1146.6 MPa at the convex triangular cross-section design of a Wave One endodontic reciprocating file, followed by the S-shaped cross-section design of Reciproc ( 942.48 MPa ), the parallelogram cross-section design of WaveOne Gold ( 847.35 MPa ) and the S-shaped cross-section design of Reciproc Blue ( 614.21 MPa ) under bending conditions.
However, torsional resistance tests showed a maximum von Misses stress of 990.14 MPa at the WaveOne endodontic reciprocating files, followed by Reciproc ( 678.69 MPa ), Reciproc Blue ( 669.75 MPa ) and WaveOneGold ( 481 MPa ) [20]; however, these results are not aligned with those obtained in the present study, possibly because the authors used a different cross-section design (squared) and also a lower transverse load ( 0.1 N ) .
However, FE methods require elevated CAD knowledge to accurately reproduce the geometric characterization of the endodontic rotary file parameters as well as numerical analysis knowledge to simulate the mechanical behavior and stress distribution during bending and torsional resistance tests. In addition, some drawbacks related to FE models have been highlighted as they used to be difficult and expensive to build, complex and time-consuming; furthermore, the formulations can be complicated and they have difficulty in dealing with complex geometry problems [29]. Therefore, the analytical model proposed in this study represents a solution to generate FE models repeatably.

5. Conclusions

In this work, a new procedure for the computerized generation of the geometry of endodontic rotary files and the automated creation of finite element models for the stress analysis of endodontic files is proposed. This new procedure allows us to create fully parametrized finite element models of the endodontic files directly from their design and material parameters.
The proposed procedure has been used to implement computer software, and its performance has been demonstrated by means of numerical examples. In general, it can be said that the proposed procedure can help to increase the accuracy with which the geometry of the finite element models of the rotary files can be generated and, at the same time, simplifies the analysis process and reduces time consumption.

Author Contributions

Conceptualization, V.R.-C. and Á.Z.-M.; methodology, V.R.-C. and F.S.-M.; software, V.R.-C. and F.S.-M.; investigation, V.R.-C. and Á.Z.-M.; resources, Ó.A.E., A.A.M. and A.G.C.; writing—original draft preparation, V.R.-C. and Á.Z.-M.; writing—review and editing, V.R.-C. and Á.Z.-M.; visualization, V.R.-C.; project administration, Ó.A.E., A.A.M. and A.G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available on request due to restrictions (privacy and ethical).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Basheer Ahamed, S.; Vanajassun, P.; Rajkumar, K.; Mahalaxmi, S. Comparative Evaluation of Stress Distribution in Experimentally Designed Nickel-titanium Rotary Files with Varying Cross Sections: A Finite Element Analysis. J. Endod. 2018, 44, 654–658. [Google Scholar] [CrossRef]
  2. He, R.; Ni, J. Design improvement and failure reduction of endodontic files through finite element analysis: Application to V-Taper file designs. J. Endod. 2010, 36, 1552–1557. [Google Scholar] [CrossRef]
  3. Parashos, P.; Gordon, I.; Messer, H. Factors influencing defects of rotary nickel-titanium endodontic instruments after clinical use. J. Endod. 2004, 30, 722–725. [Google Scholar] [CrossRef] [PubMed]
  4. Spili, P.; Parashos, P.; Messer, H. The impact of instrument fracture on outcome of endodontic treatment. J. Endod. 2005, 31, 845–850. [Google Scholar] [CrossRef] [PubMed]
  5. Sattapan, B.; Nervo, G.; Palamara, J.; Messer, H. Defects in rotary nickel-titanium files after clinical use. J. Endod. 2000, 26, 161–165. [Google Scholar] [CrossRef] [Green Version]
  6. Barbosa, F.; Ponciano Gomes, J.; de Araújo, M. Fractographic Analysis of K3 Nickel-Titanium Rotary Instruments Submitted to Different Modes of Mechanical Loading. J. Endod. 2008, 34, 994–998. [Google Scholar] [CrossRef]
  7. Plotino, G.; Grande, N.; Cordaro, M.; Testarelli, L.; Gambarini, G. A Review of Cyclic Fatigue Testing of Nickel-Titanium Rotary Instruments. J. Endod. 2009, 35, 1469–1476. [Google Scholar] [CrossRef]
  8. Inan, U.; Gonulol, N. Deformation and Fracture of Mtwo Rotary Nickel-Titanium Instruments After Clinical Use. J. Endod. 2009, 35, 1396–1399. [Google Scholar] [CrossRef] [PubMed]
  9. McGuigan, M.; Louca, C.; Duncan, H. Clinical decision-making after endodontic instrument fracture. Br. Dent. J. 2013, 214, 395–400. [Google Scholar] [CrossRef] [Green Version]
  10. Mcguigan, M.; Louca, C.; Duncan, H. The impact of fractured endodontic instruments on treatment outcome. Br. Dent. J. 2013, 214, 285–289. [Google Scholar] [CrossRef] [Green Version]
  11. Yang, Q.; Shen, Y.; Huang, D.; Zhou, X.; Gao, Y.; Haapasalo, M. Evaluation of Two Trephine Techniques for Removal of Fractured Rotary Nickel-titanium Instruments from Root Canals. J. Endod. 2017, 43, 116–120. [Google Scholar] [CrossRef]
  12. Strindberg, L. The Dependence of the Results of Pulp Therapy on Certain Factors: An Analytic Study Based on Radiographic and Clinical Follow-up Examinations. Acta Dontol. Scand. 1956, 14 (Suppl. 21), 1–175, N.G. Mauritzons Boktryckerl. [Google Scholar]
  13. Kim, H.; Kim, H.; Lee, C.; Kim, B.; Park, J.; Versluis, A. Mechanical response of nickel-titanium instruments with different cross-sectional designs during shaping of simulated curved canals. Int. Endod. J. 2009, 42, 593–602. [Google Scholar] [CrossRef] [PubMed]
  14. Faus-Llácer, V.; Kharrat, N.H.; Ruiz-Sánchez, C.; Faus-Matoses, I.; Zubizarreta-Macho, Á.; Faus-Matoses, V. The Effect of Taper and Apical Diameter on the Cyclic Fatigue Resistance of Rotary Endodontic Files Using an Experimental Electronic Device. Appl. Sci. 2021, 11, 863. [Google Scholar] [CrossRef]
  15. Sekar, V.; Kumar, R.; Nandini, S.; Ballal, S.; Velmurugan, N. Assessment of the role of cross section on fatigue resistance of rotary files when used in reciprocation. Eur. J. Dent. 2016, 10, 541–545. [Google Scholar] [CrossRef] [Green Version]
  16. Kwak, S.; Ha, J.H.; Lee, C.J.; El Abed, R.; Abu-Tahun, I.; Kim, H.C. Effects of Pitch Length and Heat Treatment on the Mechanical Properties of the Glide Path Preparation Instruments. J. Endod. 2016, 42, 788–792. [Google Scholar] [CrossRef] [PubMed]
  17. Hülsmann, M. Research that matters: Studies on fatigue of rotary and reciprocating NiTi root canal instruments. Int. Endod. J. 2019, 52, 1401–1402. [Google Scholar] [CrossRef]
  18. Zienkiewicz, O.; Taylor, R.; Zhu, J. The Finite Element Method: Its Basis and Fundamentals; Elsevier Science: Amsterdam, The Netherlands, 2005. [Google Scholar]
  19. Berutti, E.; Chiandussi, G.; Gaviglio, I.; Ibba, A. Comparative analysis of torsional and bending stresses in two mathematical models of nickel-titanium rotary instruments: ProTaper versus ProFile. J. Endod. 2003, 29, 15–19. [Google Scholar] [CrossRef]
  20. Prados-Privado, M.; Rojo, R.; Ivorra, C.; Prados-Frutos, J. Finite element analysis comparing WaveOne, WaveOne Gold, Reciproc and Reciproc Blue responses with bending and torsion tests. J. Mech. Behav. Biomed. Mater. 2019, 90, 165–172. [Google Scholar] [CrossRef]
  21. Scattina, A.; Alovisi, M.; Paolino, D.; Pasqualini, D.; Scotti, N.; Chiandussi, G.; Berutti, E. Prediction of cyclic fatigue life of nickel-titanium rotary files by virtual modeling and finite elements analysis. J. Endod. 2015, 41, 1867–1870. [Google Scholar] [CrossRef] [Green Version]
  22. De Arruda Santos, L.; De Azevedo Bahia, M.; De Las Casas, E.; Buono, V. Comparison of the mechanical behavior between controlled memory and superelastic nickel-titanium files via finite element analysis. J. Endod. 2013, 39, 1444–1447. [Google Scholar] [CrossRef] [PubMed]
  23. Lou, Y.; Li, Z.; Hervé, J. Projective, Affine and Euclidean Geometric Transformations and Mobility in Mechanisms; LAP Lambert Academic Publishing: Chisinau, Moldova, 2011; pp. 1–18. [Google Scholar]
  24. Shewchuk, J. General-dimensional constrained delaunay and constrained regular triangulations, I: Combinatorial properties. Discret. Comput. Geom. 2008, 39, 580–637. [Google Scholar] [CrossRef]
  25. Si, H. TetGen, a Delaunay-Based Quality Tetrahedral Mesh Generator. ACM Trans. Math. Softw. 2015, 41, 1–36. [Google Scholar] [CrossRef]
  26. ISO 3630-1:2008-Dentistry-Root-Ccanal Instruments-Part 1: General Requirements and Test Methods; Standard, International Organization for Standardization: Geneva, Switzerland, 2008.
  27. Bonessio, N.; Pereira, E.; Lomiento, G.; Arias, A.; Bahia, M.; Buono, V.; Peters, O. Validated finite element analyses of WaveOne Endodontic Instruments: A comparison between M-Wire and NiTi alloys. Int. Endod. J. 2015, 48, 441–450. [Google Scholar] [CrossRef]
  28. De Arruda Santos, L.; López, J.; De Las Casas, E.; De Azevedo Bahia, M.; Buono, V. Mechanical behavior of three nickel-titanium rotary files: A comparison of numerical simulation with bending and torsion tests. Mater. Sci. Eng. C 2014, 37, 258–263. [Google Scholar] [CrossRef]
  29. Montalvão, D.; Alçada, F.; Braz Fernandes, F.; De Vilaverde-Correia, S. Structural characterisation and mechanical fe analysis of conventional and m-wire Ni-Ti alloys used in endodontic rotary instruments. Sci. World J. 2014, 2014. [Google Scholar] [CrossRef] [Green Version]
  30. MATLAB. 9.4.0.813654 (R2018a); The MathWorks Inc.: Natick, MA, USA, 2018. [Google Scholar]
  31. Abaqus Theory Guide; Dassault Systèmes Simulia Corp: Johnston, RI, USA, 2015.
  32. Zmudzki, J.; Walke, W.; Chladek, W. Influence of model discretization density in FEM numerical analysis on the determined stress level in bone surrounding dental implants. Adv. Soft Comput. 2008, 47, 559–567. [Google Scholar] [CrossRef]
  33. Zmudzki, J.; Walke, W.; Chladek, W. Stresses present in bone surrounding dental implants in FEM model experiments. J. Achiev. Mater. Manuf. Eng. 2008, 27-1, 71–74. [Google Scholar]
  34. Zienkiewicz, O.C.; Zhu, J.Z. A simple error estimator and adaptive procedure for practical engineerng analysis. J. Achiev. Mater. Manuf. Eng. 1987, 24, 337–357. [Google Scholar] [CrossRef]
  35. Patil, H.; Jeyakarthikeyan, P. Mesh convergence study and estimation of discretization error of hub in clutch disc with integration of ANSYS. IOP Conf. Ser. Mater. Sci. Eng. 2018, 402, 012065. [Google Scholar] [CrossRef]
  36. Stolk, J.; Verdonschot, N.; Huiskes, R. Management of stress fields around singular points in a finite element analysis. Comput. Methods Biomech. Biomed. Eng. 2001, 3, 57–62. [Google Scholar]
  37. Chen, G.; Pettet, G.; Pearcy, M.; McElwain, D. Comparison of two numerical approaches for bone remodelling. Med. Eng. Phys. 2007, 29, 134–139. [Google Scholar] [CrossRef] [PubMed]
  38. Versluis, A.; Kim, H.C.; Lee, W.; Kim, B.M.; Lee, C.J. Flexural Stiffness and Stresses in Nickel-Titanium Rotary Files for Various Pitch and Cross-sectional Geometries. J. Endod. 2012, 38, 1399–1403. [Google Scholar] [CrossRef] [PubMed]
  39. Arbab-Chirani, R.; Chevalier, V.; Arbab-Chirani, S.; Calloch, S. Comparative analysis of torsional and bending behavior through finite-element models of 5 Ni-Ti endodontic instruments. Oral Surg. Oral Med. Oral Pathol. Oral Radiol. Endodontol. 2011, 111, 115–121. [Google Scholar] [CrossRef] [PubMed]
  40. Baek, S.H.; Lee, C.J.; Versluis, A.; Kim, B.M.; Lee, W.; Kim, H.C. Comparison of torsional stiffness of nickel-titanium rotary files with different geometric characteristics. J. Endod. 2011, 37, 1283–1286. [Google Scholar] [CrossRef] [PubMed]
  41. Hamdy, T.M.; Galal, M.; Ismail, A.G.; Abdelraouf, R.M. Evaluation of Flexibility, Microstructure and Elemental Analysis of Some Contemporary Nickel-Titanium Rotary Instruments. Open Access Maced. J. Med Sci. 2019, 7, 3647–3654. [Google Scholar] [CrossRef] [Green Version]
  42. Omar, N.; Ismail, A.; Galal, M.; Zaazou, M.H.; Mohamed, M.D.A. A comparative finite analysis of the mechanical behavior of ProTaper NEXT and WaveOne rotary files. Bull. Natl. Res. Cent. 2019, 43, 1–6. [Google Scholar] [CrossRef] [Green Version]
  43. Galal, M.; Hamdy, T.M. Evaluation of stress distribution in nickel-titanium rotary instruments with different geometrical designs subjected to bending and torsional load: A finite element study. Bull. Natl. Res. Cent. 2020, 44, 1–11. [Google Scholar] [CrossRef]
  44. Tsao, C.; Liou, J.; Wen, P.; Peng, C.; Liu, T. Study on bending behaviour of nickel-titanium rotary endodontic instruments by analytical and numerical analyses. Int. Endod. J. 2013, 46, 379–388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Definition of the geometry of the endodontic file (a), the raw material (b) and the machined material (c).
Figure 1. Definition of the geometry of the endodontic file (a), the raw material (b) and the machined material (c).
Applsci 11 04329 g001
Figure 2. Cross-sections of the endodontic file at different positions along the z g axis: active part (a), transition area (b) and shaft (c).
Figure 2. Cross-sections of the endodontic file at different positions along the z g axis: active part (a), transition area (b) and shaft (c).
Applsci 11 04329 g002
Figure 3. (Left) Definition of coordinate systems S f and S g and (Right) definition of a sample squared cross section for the endodontic file.
Figure 3. (Left) Definition of coordinate systems S f and S g and (Right) definition of a sample squared cross section for the endodontic file.
Applsci 11 04329 g003
Figure 4. Definition of loading and boundary conditions.
Figure 4. Definition of loading and boundary conditions.
Applsci 11 04329 g004
Figure 5. Graphical interface for the developed computer program.
Figure 5. Graphical interface for the developed computer program.
Applsci 11 04329 g005
Figure 6. Finite element mesh for the endodontic file; the average element size is 0.10 mm ( n θ = 24 and n v = 250 ).
Figure 6. Finite element mesh for the endodontic file; the average element size is 0.10 mm ( n θ = 24 and n v = 250 ).
Applsci 11 04329 g006
Figure 7. Results for the convergence study: (a) error energy norm of the whole model and (b) maximum element energy error.
Figure 7. Results for the convergence study: (a) error energy norm of the whole model and (b) maximum element energy error.
Applsci 11 04329 g007
Figure 8. Stress results for torsional analysis (undeformed plot).
Figure 8. Stress results for torsional analysis (undeformed plot).
Applsci 11 04329 g008
Figure 9. Stress results for bending analysis (undeformed plot).
Figure 9. Stress results for bending analysis (undeformed plot).
Applsci 11 04329 g009
Figure 10. Stress results for bending analysis: (a) detail of the region where maximum von Mises stress takes place and (b) evolution of the von Mises stresses along the edge.
Figure 10. Stress results for bending analysis: (a) detail of the region where maximum von Mises stress takes place and (b) evolution of the von Mises stresses along the edge.
Applsci 11 04329 g010
Figure 11. Results of the parametric study of the pitch of the active part: (a) bending analysis and (b) torsional analysis.
Figure 11. Results of the parametric study of the pitch of the active part: (a) bending analysis and (b) torsional analysis.
Applsci 11 04329 g011
Figure 12. Results of the parametric study of the length of the active part: (a) bending analysis and (b) torsional analysis.
Figure 12. Results of the parametric study of the length of the active part: (a) bending analysis and (b) torsional analysis.
Applsci 11 04329 g012
Figure 13. Results of the parametric study of the diameter of the tip of the active part: (a) bending analysis and (b) torsional analysis.
Figure 13. Results of the parametric study of the diameter of the tip of the active part: (a) bending analysis and (b) torsional analysis.
Applsci 11 04329 g013
Table 1. Definition of cases of study.
Table 1. Definition of cases of study.
NameMagnitude
Cross sectionSquared
Diameter of the shaft, d s h 1.20 mm
Diameter of the tip, d a p 0.25 mm
Length of the active part, L a p 16.00 mm
Total length of the file, L t o t a l 25.00 mm
Taper, c 6 %
Pitch of the active part, p z 2.00 mm
Parabola coefficient transition area, a s h 0.10 mm 1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Roda-Casanova, V.; Zubizarreta-Macho, Á.; Sanchez-Marin, F.; Alonso Ezpeleta, Ó.; Albaladejo Martínez, A.; Galparsoro Catalán, A. Computerized Generation and Finite Element Stress Analysis of Endodontic Rotary Files. Appl. Sci. 2021, 11, 4329. https://doi.org/10.3390/app11104329

AMA Style

Roda-Casanova V, Zubizarreta-Macho Á, Sanchez-Marin F, Alonso Ezpeleta Ó, Albaladejo Martínez A, Galparsoro Catalán A. Computerized Generation and Finite Element Stress Analysis of Endodontic Rotary Files. Applied Sciences. 2021; 11(10):4329. https://doi.org/10.3390/app11104329

Chicago/Turabian Style

Roda-Casanova, Victor, Álvaro Zubizarreta-Macho, Francisco Sanchez-Marin, Óscar Alonso Ezpeleta, Alberto Albaladejo Martínez, and Agustín Galparsoro Catalán. 2021. "Computerized Generation and Finite Element Stress Analysis of Endodontic Rotary Files" Applied Sciences 11, no. 10: 4329. https://doi.org/10.3390/app11104329

APA Style

Roda-Casanova, V., Zubizarreta-Macho, Á., Sanchez-Marin, F., Alonso Ezpeleta, Ó., Albaladejo Martínez, A., & Galparsoro Catalán, A. (2021). Computerized Generation and Finite Element Stress Analysis of Endodontic Rotary Files. Applied Sciences, 11(10), 4329. https://doi.org/10.3390/app11104329

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop