Issue 
Int. J. Simul. Multidisci. Des. Optim.
Volume 12, 2021
Advances in Modeling and Optimization of Manufacturing Processes



Article Number  32  
Number of page(s)  9  
DOI  https://doi.org/10.1051/smdo/2021032  
Published online  17 November 2021 
Research Article
Numerical simulation and experimentation of endodontic file using Taguchi DoE
Vel Tech Rangarajan Dr. Sagunthala R&D Institute of Science and Technology, Chennai 600062, India
^{*} email: drsalunkhesachin@veltech.edu.in
Received:
7
January
2021
Accepted:
26
October
2021
The endodontic file is a tapered, needle shape body used for the preparation of curved human root canals. During the preparation, process files get failed due to the locking action offered by the canal wall. The present study aimed to find the fatigue life of endodontic files at 23°, 33° and 43° root canal curvature angles. Four brands of files were selected for the present study viz. Hyplex CM, ProTaper Next, Hero Shaper, ProFile Vortex. The strain life analysis done using ANSYS showed that the Hyplex CM file gives the highest fatigue life at 23°, 33° and 43° root canal curvature angles. Therefore, Hyplex CM file was selected for DoE Taguchi Optimization study. Each experimental reading was conducted on X Smart Plus experimental setup under purely rotary and combined reciprocatingrotary motion, at 500 rpm, 600 rpm and 700 rpm and 23°, 33° and 43° root canal curvature angles. Hyplex CM file given maximum fatigue life at 23° root canal curvature angle and 500 rpm speed of rotation. But if file motion is combined rotaryreciprocating, life reduces. The ANNOVA study showed that Pvalue and significance F are very small, which represented that the regression model is effective.
Key words: Endodontic file / fatigue life / root canal / DoE / ANOVA
© P.R. Lokhande et al., Published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The endodontic treatment aims to save the tooth which is infected, decayed or highly susceptible to future infection. The endodontic treatment is also known as root canal treatment [1]. The cleaning and shaping of the infected root canal are known as root canal preparation and filling and sealing of the prepared root canal using biocompatible material Guttapercha is known as root canal obturation, both root canal preparation and obturation comes under endodontic treatment and is known as root canal treatment [2,3]. The endodontic files are used to prepare the root canals these files are tapered drill bitshaped instruments having continuous helical grooves called flutes. Flutes help debris to come out from the infected root canal during the preparation process [4]. The endodontic file rotates in a curved root canal in presence of infected tissue, debris and oral fluids which offers locking pressure to the rotating file. The locking action of the canal on the endodontic file leads to fatigue failure [5]. This fatigue failure can be assessed using cyclic fatigue or torsion fatigue approach. These approaches help to find the life of the endodontic file. The life of the endodontic file helps the Dentist to make use of the file so that failure of the file in the patient jaw can be avoided during canal preparation.
The present study is divided into two parts viz. (i) numerical simulation of the endodontic file and (ii) optimization of the endodontic file using Taguchi Design of Experiments (DoE). The study used finite element analysis and cyclic fatigue test for assessment of fatigue life of endodontic files during the root canal preparation process when placed in 23°, 33° and 43° root canal curvature angles and under purely rotary, combined reciprocatingrotary motions. The Taguchi DoE optimization was used to find optimized root canal angle and speed of rotation giving maximum fatigue life of the file. Regression analysis and ANNOVA were used for statistical analysis of obtained results.
2 Literature review
Fatigue is progressive, confined structural destruction caused to Endodontic files subjected to periodic loading. If applied loads are higher than the design load, a microscopic crack begins forming stress concentrators such as slip bands which propagate suddenly to cause fatigue failure. The stress variations in the NiTi rotary instruments were conducted and it failed to evaluate the life of files. The study used the virtual model to estimate the life of the file and found the locations of failure [6]. The finite element method for analysis of cyclic performance NiTi rotary files made up of different materials was studied [7]. There is a need to measure the resistance offered by the file to avoid fractures [8]. The instrument fractures are dependent of root canal curvature angle, operating speed, torsion, compression and tension, Endodontic file properties, instrument manufacturing methods/defects, instruments geometry and material fatigue [9]. The stress state for the Endodontic file instrument is affected by many parameters viz. size, taper, crosssectional area, rigidity, flexibility, the direction of rotation of files and forces acting on the rotary instruments during working [10]. The study reported that different methodologies affected the fatigue life of the instruments and therefore the clinical outcome. It is suggested to carry the international standard for similarity in methodology and comparable results [11].
3 Materials and method
This section describes two methods viz. (i) Numerical simulation of the endodontic file and (ii) Optimization of the endodontic file using Taguchi Design of Experiments (DoE). The numerical simulation of the endodontic file helped to find out the fatigue life of endodontic files of four brands viz. Hyplex CM, ProTaper Next, Hero Shaper, ProFile Vortex. Out of these four brands of files, the file with the highest fatigue life (best file) was selected for DoE optimization. Taguchi DoE has given the best parameter combination that suits maximum fatigue life.
3.1 Numerical simulation of endodontic file
The cyclic fatigue life assessments of Endodontic files have been carried out using finite element analysis for rotary motion at 23°, 33° and 43° root canal curvature angles. The Endodontic file model of each brand viz. Hyplex CM, ProTaper Next, Hero Shaper, ProFile Vortex were imported to ANSYS to find out the cyclic fatigue by varying root canal curvature angles. ANSYS is the FEAbased analysis software that calculated the fatigue life of endodontic files using the strainlife fatigue analysis approach. In any patient, the root canal angle varies from 0° to 45°. The present study considered 23°, 33° and 43° root canal curvature angles.
3.1.1 Measurement of dimensions and modeling endodontic file
Profilometry extracts the topological data from the surface of the component. The main aim of profilometry is to obtain surface morphology, length, breadth, depth, crest, trough, pitch, orientation angle. The view of the endodontic file on the optical profilemeter is shown in Figure 1.
The CAD models of Hyplex CM endodontic file was needed to perform fatigue life assessment in ANSYS. The dimensions of the file were measured using optical profilemeter was used to model the endodontic file (shown in Fig. 1). These CAD models of file were used for assessment of fatigue life in ANSYS shown in Figure 2.
Fig. 1
Endodontic file profile view obtained on optical profile meter. 
Fig. 2
CAD model of hyplex CM endodontic file. 
3.1.2 Material property inputs and meshing of endodontic file
In order to give the input of material properties of nickeltitanium alloy to ANSYS software, various tests need to be performed to obtain the properties of the material. The past studies have been already investigated these material properties and used them for their study [5]. The material properties are shown in Tables 1 and 2.
Any continuous object has infinite degrees of freedom and it's not just possible to solve the problem in this format. Finite element analysis reduces the degree of freedom from infinite to finite with the help of discretization or meshing. The meshing divided the endodontic file CAD model into 47230 nodes and 13118 elements. Meshing used two types of elements viz. 10noded tetrahedral elements (Solid187) and 20noded hexahedral elements (Solid186). Looking at the shape of the tetrahedral element, it can be easily adjusted on the curved/complex profile of the endodontic file. Therefore, the tetrahedral element was selected for the meshing body of the endodontic file. On the other hand, the handle is simple in structure and therefore, was meshed using hexahedral elements. The mesh model of Hyplex CM file is shown in Figure 3.
Material property for nickeltitanium alloy.
Linear elastic properties for nickeltitanium alloy.
Fig. 3
Meshed model of hyplex CM endodontic file. 
3.1.3 Boundary conditions input for endodontic file
Defining boundary conditions and loading conditions means fixing the value of displacement and load on a specific area of the model. Thus, the model correctly constrained can give highly accurate results. The endodontic file was fixed at the handle end and a force of 0.24 N was applied to cause deflection of the tip by 1.88 mm (shown in Fig. 4). When the endodontic file was inserted into 43° root canal its tip gets deflected by 1.88 mm from the axis line. For simulation of this deflection in ANSYS 0.24 N force was applied to cause 1.88 mm deflection which constitutes 43° root canal angle. The experimental study used an acrylic block in which the root canal cavity was made to simulate the root canal but in ANSYS we avoiding this which saved the time required to build the root canal CAD model.
Fig. 4
Boundary conditions input in hyplex CM endodontic file. 
3.1.4 Strain life analysis using ANSYS
The ANSYS is used to perform the strainlife fatigue analysis. Strainlife analysis requires the input of nickeltitanium alloys material properties such as strength coefficient, strength exponent, ductility coefficient, ductility exponent, cyclic strength coefficient and cyclic strain hardening exponent. The strain analysis provides the fatigue life in terms of the number of reversals. The fatigue life of the endodontic file can be calculated using the following equation.(1)
The geometric nonlinearity considered in the static analysis: Under solver controls ‘Large Deflection’ option selected. Fully reversed type loading was used as shown in Figure 5. The various analysis parameters inputs for analysis used includes Fatigue Stress Concentration factor = 1; Type of approach used in fatigue analysis: Strainlife (The strain life curve for analysis is shown in Fig. 6); Stress component considered: Equivalent Stress; No mean stress theory was considered. The other strain life parameters are tabulated in Table 3.
Fig. 5
Reversible loading. 
Fig. 6
Strainlife curve. 
Strain life parameters considered for the nickeltitanium file material.
3.1.5 Fatigue life of hyplex CM endodontic file based on ANSYS report
The infinite fatigue life was defined at 10^{7} cycles or above. The strainlife approach measures life in terms of reversals. The number of reversals to failure is equal to onehalf of the fatigue life of the file. The minimum number of reversals observed for hyplex CM endodontic file is 5055 (shown in Fig. 7). So, the life of Hyplex CM endodontic file in the terms of the number of cycles is 2527 number of rotations.
Similarly, the analysis procedure described for Hyplex CM endodontic file was followed for fatigue analysis of remaining files viz. ProTaper Next, Hero Shaper and ProFile Vortex at 23°, 33° and 43° root canal curvature angle. The results were noted in Table 4 and it is found that Hyplex CM endodontic files show the highest fatigue life at 23°, 33° and 43° root canal curvature angle. Therefore, Hyplex CM file was selected for DoE Taguchi Optimization study.
Fig. 7
Fatigue life for hyplex CM endodontic file (isometric view). 
Effect of rotary motion on fatigue life of endodontic file.
Contributing factors and their levels.
3.2 DoE of the endodontic file using Taguchi method
The life of the endodontic file gets affected due to root canal preparation process parameters viz. root canal curvature angle, speed of endodontic file and type of endodontic file motion. Therefore, it had been decided to conduct Taguchi DoE optimization to find parameter combinations giving maximum fatigue life. The Hyplex CM nickeltitanium endodontic file of 25 mm length and 6% taper was selected. Numerical simulation study resulted that Hyplex CM endodontic files show highest fatigue life at 23°, 33° and 43° root canal curvature angle. Therefore, Hyplex CM file was selected for DoE Taguchi Optimization study.
This section describes fatigue life assessment using XSmart Plus and implementation of Taguchi DoE optimization.
3.2.1 Fatigue life assessment using XSmart Plus
The cyclic fatigue life assessments of Endodontic files had been carried out using X Smart Plus experimental setup (shown in Fig. 8) for rotary and combined rotaryreciprocating motions at 23°, 33° and 43° root canal curvature angles. During each experimental run, the endodontic file was inserted into the artificial root canal cavity in the acrylic block to measure the fatigue life. The nickeltitanium alloy files have made it easy to prepare the root canal by eliminating the difficulties available with stainlesssteel files such as canal transportation, ledges, zips and perforation [12] (Peters 2004). Therefore, nickeltitanium alloy endodontic files were used in this study. To conform with standard conditions for each experimental run new endodontic file (25 mm & 6% taper) was used.
Fig. 8
X Smart Plus experimental setup. 
3.2.2 Implementation of Taguchi DoE optimization
In this study, three root canal preparation process parameters including root canal curvature angle, speed of endodontic file and type of endodontic file motion were selected as contributing factors. These process parameters were selected after consultation with the dental practitioners. A deign of experiments was used to study the interaction between contributing factors (input variables) and fatigue life of file (output variable). The full factorial DoE was performed, which helps to adjust the input to get optimized output. The orthogonal array L18 was selected consisting of 18 experimental runs. Three levels of root canal curvature angle, speed of endodontic file and two levels for type of endodontic file motion were considered. During the preparation of root canal cavity endodontic file operates under either purely rotary or combined rotaryreciprocating motion, therefore two types of motions parameters were selected. The speeds of endodontic files were selected as per the manufacturer's catalog. The set of experimental runs shown in Table 6 corresponds to Taguchi Optimization Method.
The Taguchi method converts the experimental values to the signaltonoise (S/N) ratio [13], the noise representing undesired values (Standard deviation) and the signal representing the desired values (Mean) for the output characteristics.
The optimization of root canal preparation parameters was done using Taguchi analysis for better fatigue life of the endodontic file. In the present study, three root canal preparation process parameters including root canal curvature angle, speed of endodontic file and type of endodontic file motion were selected as contributing factors. As per L18 orthogonal array (shown in Tab. 7) experimentation was planned and executed. The experimental fatigue life values are then converted into an S/N ratio and the higher the better is used to determine the output characteristics.
The endodontic file operates in the curved root canal cavity during the preparation process. Therefore, root canal curvature angle plays a vital role as it causes the stress state in file material. In another way speed of file and type of motion also causes stress state. As mentioned earlier Taguchi method was used for better fatigue life of the endodontic file. To get an exact idea about the role of these three parameters in fatigue failure of the file, the % contribution needs to be calculated.
Experimental design with contributing factors.
Experimental design using L18 orthogonal array.
Regression analysis.
4 Result and discussion
The present study successfully conducted a numerical simulation of the endodontic file and optimization of the endodontic file using Taguchi Design of Experiments (DoE). The study has used finite element analysis and cyclic fatigue test for assessment of fatigue life of endodontic files during the root canal preparation process when placed in 23°, 33° and 43° root canal curvature angles and under purely rotary, combined reciprocatingrotary motions. Root canal curvature angles viz. 23°, 33° and 43° were selected since for any patient the angle varies from 13° to 45°. Therefore, to find cyclic fatigue resistance, three anles steps 23°, 33° and 43° were selected. The Taguchi DoE optimization was used to find optimized root canal angle and speed of rotation giving maximum fatigue life of the file. Regression analysis and ANNOVA were used for statistical analysis of obtained results.
The FEA of the endodontic files indicated that Hyplex CM file works longer during the root canal preparation process under curved canals of 23°, 33° and 43°. Therefore, Hyplex CM endodontic file was for later study viz. optimization using Taguchi Design of Experiments (DoE). The optimized root canal preparation process parameters were estimated using the Taguchi method. The Taguchi DoE helped to find the optimized parameters which give maximum fatigue life. Each experimental run was conducted using X Smart Plus experimental setup as shown in Figure 8. Regression analysis was performed to estimate the relationship between contributing and output variables to assess the closeness of experimentation with the developed model. The % contribution of parameters viz. root canal curvature angle, speed of endodontic file and type of endodontic file motion were decided using ANNOVA.
In the root canal preparation process, regression analysis was performed to estimate the relationship between contributing and output variables and to assess the closeness of experimentation with the developed model. In this study three root canal preparation process parameters viz. root canal curvature angle, speed of endodontic file and type of endodontic file motion were selected as contributing factors and fatigue life of endodontic file as output. The linear regression model is used to check if the data of fatigue life signify or not the output characteristic.
Rsquare in the regression table indicates the effectiveness of the developed model. The value of Rsquare if equals to 1 then it indicates the 100% effectiveness of model. In the present study Rsquare is more than 80%, which represents that model is effective to predict relationship between contributing factors (viz. root canal curvature angle, speed of endodontic file and type of endodontic file motion) and response (viz. life of endodontic file).
The ANNOVA was conducted for analyzing obtained data. Significance F and Pvalue in ANNOVA table, ideally, should be as small as possible. For example, in Table 9, the Pvalue is 0.36%, regression model and intercept in Table 10 are 99.64% accurate. Significance F is also very small, therefore, regression is good. The present ANNOVA analysis was conducted at 95% confidence interval. The obtained Pvalue is higher than the accepted confidence interval.
Fig. 9
Fractured samples during root canal preparation process. 
Fig. 10
Simple linear regression curve. 
Analysis of variance (regression and residual).
Analysis of variance (intercept and variable).
5 Conclusion
The present study concludes that the strain life analysis done using ANSYS showed that Hyplex CM file gives the highest fatigue life at 23°, 33° and 43° root canal curvature angles. Taguchi DoE study performed, as per L18 orthogonal array, showed that under purely rotary motion Hyplex CM file gives maximum fatigue life at 23° root canal curvature angle and 500 rpm speed of rotation. If file motion is combined rotaryreciprocating, life reduces drastically, since debries inside the root canal apply more resistance to the reciprocating motion of file than purely rotary. Concerning regression analysis done, Rsquare is 85.33%, which represents that model is effective to predict the relationship between contributing factors (viz. root canal curvature angle, speed of endodontic file and type of endodontic file motion) and response (viz. life of endodontic file). On the other hand, ANNOVA showed that Pvalue and significance F are very small, which represented that the regression model is effective.
6 Future scope
The fatigue failure of endodontic file may be affected by many parameters such as size, taper, cross sectional area, rigidity, flexibility, direction of rotation of files and forces acting on the rotary instruments during working. The same can be studied in future to cover more research insights.
References
 R.R. Ghorpade, K. Sundaram, V. Hegde, Deformation studies of endodontic obturator tip using FEM approach, Int. J. Eng. Invent. 1, 19–26 (2012) [Google Scholar]
 A. Ciftci, D.A. Vardarlı, I.S. Sonmez, Coronal microleakage of four endodontic temporary restorative materials: an in vitro study, Oral Surg. Oral Med. Oral Pathol. Oral Radiol. Endodontol. 108, 067–070 (2009) [Google Scholar]
 S. Moradi, M. Lomee, M. Gharechahim, Comparison of fluid filtration and bacterial leakage techniques for evaluation of microleakage in endodontics, Dental Res. J. 12, 109–112 (2015) [Google Scholar]
 S.E.D.M. Saber, M.M. Nagy, E. Schafer, Comparative evaluation of the shaping ability of Wave One, R eciproc and One Shape single‐file systems in severely curved root canals of extracted teeth, Int. Endodontic J. 48, 109–114 (2015) [CrossRef] [Google Scholar]
 Y. Shen, H.M. Zhou, Y.F. Zheng, B. Peng, M. Haapasalo, Current challenges and concepts of the thermomechanical treatment of nickeltitanium instruments, J. Endodontics 39, 163–172 (2013) [CrossRef] [Google Scholar]
 A. Scattina, M. Alovisi, D.S. Paolino, D. Pasqualini, N. Scotti, G. Chiandussi, E. Berutti, Prediction of cyclic fatigue life of nickeltitanium rotary files by virtual modelling and finite elements analysis, J. Endodont. 41, 867–870 (2015) [Google Scholar]
 D. Montalvao, Q. Shengwen, M. Freitas, A study on the influence of NiTi MWire in the flexural fatigue life of endodontic rotary files by using Finite Element Analysis, Mater. Sci. Eng. C 40, 172–179 (2014) [CrossRef] [Google Scholar]
 E. Pedulla, G.R.M. La Rosa, S. Boninelli, O.G. Rinaldi, E. Rapisarda, H.C. Kim, Influence of different angles of file access on cyclic fatigue resistance of Reciproc and Reciproc Blue instruments, J. Endodont. 44, 849–855 (2018) [CrossRef] [Google Scholar]
 A. Jamleh, A. Alfadley, K. Alfouzan, Vertical force induced with waveone and waveone gold systems during canal shaping, J. Endodont. 44, 412–415 (2018) [Google Scholar]
 Y. Shen, C. Tra, A. Hieawy, Z. Wang, M. Haapasalo, Effect of torsional and fatigue preloading on HyFlex EDM files, J. Endodont. 44, 643–647 (2018) [CrossRef] [Google Scholar]
 T. Ozyurek, G. Uslu, U. Inan, A comparison of the cyclic fatigue resistance of used and new glide path files, J. Endodont. 43, 477–480 (2017) [CrossRef] [Google Scholar]
 O.A. Peters, Current challenges and concepts in the preparation of root canal systems: a review, J. Endodont. 30, 559–567 (2004) [CrossRef] [Google Scholar]
 V.K. Vankanti, V. Ganta, Optimization of process parameters in drilling of GFRP composite using Taguchi method, J. Mater. Res. Technol. 3, 35–41 (2014) [CrossRef] [Google Scholar]
Cite this article as: Pravin R. Lokhande, Sachin S. Salunkhe, Sethuraman Balaguru, Numerical simulation and experimentation of endodontic file using Taguchi DoE, Int. J. Simul. Multidisci. Des. Optim. 12, 32 (2021)
All Tables
All Figures
Fig. 1
Endodontic file profile view obtained on optical profile meter. 

In the text 
Fig. 2
CAD model of hyplex CM endodontic file. 

In the text 
Fig. 3
Meshed model of hyplex CM endodontic file. 

In the text 
Fig. 4
Boundary conditions input in hyplex CM endodontic file. 

In the text 
Fig. 5
Reversible loading. 

In the text 
Fig. 6
Strainlife curve. 

In the text 
Fig. 7
Fatigue life for hyplex CM endodontic file (isometric view). 

In the text 
Fig. 8
X Smart Plus experimental setup. 

In the text 
Fig. 9
Fractured samples during root canal preparation process. 

In the text 
Fig. 10
Simple linear regression curve. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.