JMMCE  Vol.8 No.2 , March 2020
Experimental Investigation of Laser Surface Hardening of AISI 4340 Steel Using Different Laser Scanning Patterns
Laser surface transformation hardening becomes one of the most modern processes used to improve fatigue and wear properties of steel surfaces. In this process, the material properties and the heating parameters are the factors that present the most significant effects on the hardened surface attributes. The control of these factors using predictive modeling approaches to achieve desired surface properties leads to conclusive results. However, when the dimensions of the surface to be treated are larger than the cross-section of the laser beam, various laser-scanning patterns are involved. This paper presents an experimental investigation of laser surface hardening of AISI 4340 steel using different laser scanning patterns. This investigation is based on a structured experimental design using the Taguchi method and improved statistical analysis tools. Experiments are carried out using a 3 kW Nd: YAG laser source in order to evaluate the effects of the heating parameters and patterns design parameters on the physical and geometrical characteristics of the hardened surface. Laser power, scanning speed and scanning patterns (linear, sinusoidal, triangular and trochoid) are the factors used to evaluate the hardened depth and the hardened width variations and to identify the possible relationship between these factors and the hardened zone attributes. Various statistical tools such as ANOVA, correlations analysis and response surfaces are applied in order to examine the effects of the experimental factors on the hardened surface characteristics. The results reveal that the scanning patterns do not modify the nature of the laser parameters’ effects on the hardened depth and the hardened width. But they can accentuate or reduce these effects depending on the type of the considered pattern. The results show also that the sinusoidal and the triangular patterns are relevant when a maximum hardened width with an acceptable hardened depth is desired.

1. Introduction

Laser surface transformation hardening process (LSTHP) is increasingly used to heat treat mechanical parts surfaces in order to enhance their mechanical properties. By applying concentrated, short and localized thermal inputs, LSTHP can improve the wear and fatigue resistance of the treated zones [1] [2]. Compared to other surface heat treatment processes, LSTHP presents several advantages such as selective and precise hardening zones, minimization of undesirable deformation and distortion, application of very fast thermal cycles to obtain autogenous hardening, treatment for parts with complex geometry through the use of robots, and possibility of full process automation and control [3] [4].

In LSTHP, the hardened surface characteristics depend essentially on the laser heating process parameters and the thermo-physical properties of the material to be treated. For an adequate use of this process, it is necessary to better understand the operating mode and to develop robust predictive models that can be used to control and adjust the process parameters in order to estimate precise hardened surface characteristics without using the traditional long and costly experimental procedures [3] [4] [5]. The produced models can be used to predict the hardness profile attributes for various heating parameters and conditions. Two kinds of models can be developed in this perspective: numerical models and empirical models. The numerical models combine laser spot form and motion, non-linear properties of the material and the heat flow equation using the finite element method (FEM). This approach uses a heat flow equation in order to determine the temperature distribution for each heat treatment phase. The hardness is then approximated by the Ashby and Easterling equations [3] [6]. The empirical models are based on experimental data that can be collected using structured design of experiments, and analyzed using statistical tools in order to extract the arithmetical relationship between the heating parameters and the hardness profile attributes.

In LSTHP of large surfaces, the simplest technique frequently used is called multi-track treatment. This technique consists to execute several passes one next to the other until the surface is completely covered. In some cases, the number of laser beam passes can be prohibitive. The reduction of the number of passes for a given surface can be, theoretically, achieved by increasing the diameter of the beam spot. This is an unrealistic and very expensive solution. In addition, the multi-track technique generates overheated zones between every two successive passes where the hardness decreases dramatically. Therefore, it is necessary to improve the technique or to look for other alternatives [7] [8]. Researchers such as Lakhkar et al. have improved this method while trying to minimize overheated areas between hardening tracks where hardness decreases, and have succeeded to obtain a hardness that varies between 48 and 58 HRC using 850 W laser power. They have found also that hardness variation can be controlled by the spacing between the tracks [7]. But this method is still unsatisfactory since the number of passes increases rapidly by increasing the surface to be treated, which prolongs the treatment time and increases the energy supply causing the deformation of parts. Xue et al. [9] [10] have worked on another treatment technique called “Dot matrix treatment”, which generally consists of treating just a certain percentage of the surface with a uniform distribution using a pulsed laser beam. This method has reduced the energy applied to the part to be treated, thus avoiding distortion of the part, but the risk of wear and tear in untreated areas always precedes the treated areas [9] [10]. Using this method, Straffelini et al. proved that tensile strength and wear resistance increased as the density of the matrix increased [11].

There is also another interesting treatment technique that can be used and which is based on the use of complex laser beam moving patterns during the laser hardening process. The use of certain types of trajectories allows mainly to treat the maximum of the surface with a single pass, this allows to optimize the laser treatment process of the relatively large geometries. To date, few pieces of research have addressed this method although the energy supply is enormously optimized and the treatment time is reduced. Sistaninia et al. are almost among the only researchers who have studied this method of treatment [12] - [17]. However, their efforts were focused on the development of analytical models that analyze the temperature distribution obtained by different beam displacement patterns, which explains the importance of an experimental analysis of the effect of this technique on the obtained hardness profile. Sistaninia and Farrahi have determined that for given laser power, a sinusoidal moving trajectory has a more uniform temperature distribution than a trochoid trajectory [12] [16],and Zhou et al. have shown that an increase in frequency leads to a drop in the maximum temperature [13] [17]. Therefore, an experimental investigation of interactions between the different phenomena during the hardening process is necessary to develop predictive models able to control the parameters of this heat treatment technique in order to produce the appropriate hardened surface characteristics.

The objective of this paper is to evaluate the effect of using different displacement trajectories during the laser heat treatment process of AISI 4340 steel plates. This is done by conducting an experimental investigation of the moving pattern effect and the process parameters’ effect on the hardness profile. This study allows to determine the trajectory that will allow to have the best hardness profile (best depth and hardness distribution) to optimize the process. The investigation is based on a structured design of experiments and various statistical analysis tools. Experiments are carried out by using a 3 kW Nd: YAG laser source where four moving patterns are considered: linear, sinusoidal, triangular and trochoid pattern. The obtained results permit to understand the effects of patterns as well as laser parameters on the hardened zone attributes.

2. Experimentation

2.1. Methodology

The experimental investigation is conducted on AISI 4340 steel 60 mm long, 6 mm wide, 5 mm thick plates hardened and tempered to 25 - 30 HRC (Figure 1). The parts were oil quenched after complete austenitization in the furnace, then quenched to blacken their color and to improve the absorptivity of the laser irradiation. The chemical composition of the material is shown in Table 1. To carry out the various experiments, an Nd: YAG laser source was used, powered by a YLS-3000-ST2 IPG fiber laser and mounted on a six-axis FANUC M-710IC robot. The ILVDC-Scanner system mounted on the laser head is used to perform the desired scanning patterns. The used laser source has a maximum power of 3000 W, a wavelength of 1070 nm and a focal length of 310 mm. The diameter of the focal spot used in this study is 0.52 mm, and it can be modified by means of a collimator fixed on the laser head. The main elements of the laser cell are presented in Figure 2. Figure 3 illustrates the four moving patterns considered in this study.

Figure 1. A schematic illustration of laser transformation hardening.

Figure 2. Experimental setup of steel plate’s laser surface hardening using laser cell.

Figure 3. Schematic illustration of the considered laser scanning patterns.

Table 1. Chemical composition of AISI 4340 steel.

The main parameters selected in the present study are laser power (P), Scanning speed (SS) and type of moving pattern. As shown in Figure 4, the measured quantities are the hardened depth (HD) and the hardened width (HW), which represent two properties to be optimized by the laser heat treatment process. After the laser processing, the samples are analyzed using a standard metallographic approach that consists of carefully cutting the plates using an Abrasi-Matic300 cutting machine, then coating and polishing the samples, then measuring the hardness of the parts and finally performing the chemical etching to measure the hardened depth and the hardened width. The chemical attack was performed using Nital’s chemical solution (95% ethanol and 5% nitric acid). The hardened depth and hardened width measurements are taken by a CLEMEX precision microscope, and the hardness profile is determined based on the micro-hardness measurements using the Clemex machine.

2.2. Preliminary Tests

Due to the lack of a database in the literature that provides the ranges of variation of the laser hardening process parameters, preliminary tests were carried out to define its intervals in order to ensure a hardened surface with a minimal melting zone. By varying one parameter at a time, the combination of maximum power and minimum scanning speed allows to have maximum hardened depth and maximum hardened width but with some surface melting. On the other hand, the combination of minimum power and maximum scanning speed allows to have a minimum hardening depth and width but without melting the surface. Consequently, the minimum and maximum limits for each processing parameter are defined based on these preliminary tests, where the variation range of the power is selected between 500 and 800 W and the speed varies between 20 and 40 mm/s.

2.3. Experimental Design

To establish experimental tests, it is essential to use an experimental design that

Figure 4. Cross-section of the steel plate with the variables to be measured.

can provide relevant results with the minimum number of tests to be performed, while respecting the ranges of variation already defined by the preliminary tests. Depending on the number of factors and their assigned levels, a systematic design based on factor design is an excellent tool to use and leads to an understanding of the material's behavior towards this processing technique. Due to their efficiency and simplicity, factorial designs are the most advised recommended designs to use, although the number of tests increases exponentially when an additional parameter or level is adopted at design, making the study more expensive and time-consuming [18]. However, fractional designs can collect statistically significant data on the studied phenomenon, on the effects of the process parameters and even on the interactions between the factors with the minimum number of tests [18]. Consequently, a Taguchi plan was established as an experimental design for this study. For each moving trajectory, 2 factors, 3-levels L9 matrix is used in order to determine the effect of power and scanning velocity on the hardened depth and hardened width of the steel plates. Indeed, the choice to introduce three levels for each parameter is intended to study the nonlinearity of the parameter effects and the interactions between them. The selected factors and their levels are identified in Table 2. The experimental design that best fit this problem is illustrated in Table 3.

2.4. Preliminary Tests and Repeatability

To improve the accuracy and validity of the results, four tests were performed under the same working conditions (power of 650 W and a scanning speed of 30 mm/s for a linear track), which also allows to estimate the overall and relative errors of measurements. These errors can be caused by uncontrollable factors such as a sudden change in ambient temperature during processing, non-uniformity of the chemical composition of the base material, incorrect calibration or use of measuring instruments, etc. The results presented in Table 4 show good precision of the measurement since the variations are less than 6%. These results then guarantee the validity of the measurement method and make it possible to carry out the final tests with confidence.

3. Results Analysis and Discussion

The use of different laser scanning patterns leads to different hardened surfaces

Table 2. Selected parameters and their levels for experimentation.

Table 3. Experimental design used for each laser scanning patterns.

Table 4. Results of the repetition tests.

topology in terms of hardened depth, hardened width and hardness homogeneity. A slight surface fusion is also present in some cases but in a negligible and acceptable gradation, since the selection of parameter variation ranges during the preliminary tests contributed significantly to minimize this phenomenon. Figure 5 shows typical hardness profiles achieved using various laser parameters applied to the traditional linear scanning pattern. Images 5a, 5b and 5c present the hardness profiles achieved using scanning velocity of 20 mm/s for power of 500, 650 and 800 W respectively. Images 5d, 5e and 5f illustrate the hardness profiles at a fixed power of 800 W and at scanning speed of 20, 30 and 40 mm/s, respectively. Figure 5(a) illustrates also the difference between the base metal and the hardened zone at the microscopic level where it is clear that the metallographic structure of the treated area is denser than that of the base metal. This density is the result of the immigration of carbon grains into the austenite during heating by the laser beam and its transformation into martensite during cooling. From these figures, it is clear as expected that the laser power has a positive effect on the hardness zone dimensions, while the increase of the scanning velocity has a negative effect. These observations confirm that both factors have direct effects on the hardness zone in terms of both hardened depth and hardened width. However, these effects that are the consequences of spatial

Figure 5. Typical hardened zones achieved using linear scanning pattern at various laser parameters.

temperature distribution can be completely transformed when using others scanning patterns. This transformation modifies considerably the thermal topology of the hardened surface.

Figure 6 shows typical hardness profiles achieved using various laser scanning patterns. The variation of HD and HW from a scanning pattern to another is directly related to pattern specifications such as amplitude, frequency and focal spot. The figure describes the variation of hardness profile for test 7 that uses a scanning velocity of 20 mm/s and laser power of 800 W. The scanning pattern geometry affects the exposure time of different portions of the part surface to the laser radiation. Depending on the scanning pattern, some zones of the surface can be exposed to the heat source repeatedly and for a longer time. An increase of the exposure interval gives more time for the carbon to migrate homogeneously into the austenite producing a larger treated surface. However, this time must not exceed certain limits to avoid melting zones and unwanted deformations of the treated surface.

Moreover, Figure 6 shows that the sinusoidal pattern presents the best compromise between hardened depth, hardened width and hardness distribution uniformity. Trochoid pattern produces high HD and HW but with a very irregular hardness distribution. The triangular pattern produces HD and HW with a uniform distribution but remains less than that achieved by sinusoidal pattern.

Figure 6. Typical hardened zones using a power of 800 W and a scanning speed of 20 mm/s for linear, sinusoidal, triangular and trochoid scanning patterns respectively.

Linear pattern presents, as anticipated, the maximum hardened depth and a good hardness distribution but with the lowest HW. In the linear scanning pattern, the laser power remains concentrated below the focal spot. HD is maximum but HW slightly exceeds the beam diameter. For the other patterns, the laser power is distributed over the entire amplitude of the trajectory increasing the width of the hardened surface. The increase of the HW results inexorably in a decrease of HD.

Table 5 presents the hardened depths and hardened widths collected from the experiments. These results are analyzed using several statistical indicators extracted from the analysis of variance (ANOVA) such as percent contribution of each parameter in HD and HW variations, the average effects graphs, the correlation between the measured characteristics and the laser parameters, and the graphs of the response surfaces. Generally, the parameter with the highest percentage contribution represents the most significant parameter. The total sum of the contributions of the considered parameters should reach 100%. The difference is usually attributed to uncontrolled factors or experimental errors.

The graphs of average effects consist of visualizing the effect of each parameter on the hardened zone characteristics. These graphs allow easy comparisons between the parameter’s effects and a better understanding of the connections between these parameters and the characteristics to evaluate. The response surface method (RSM) is adopted to explore the nature of relationships between the experimental factors and the hardness profiles.

The results of ANOVA performed on the four L9 are presented in Table 6. The analysis of the parameter contribution reveals different results between linear pattern and other scanning patterns. Generally, the laser power affects HD more than the scanning speed regardless of the scanning pattern. However, the scanning speed affects HW more than the power. An increase in power

Table 5. Experimental results achieved using the four performed L9.

Table 6. ANOVA analysis for the hardened depth and width for the four moving patterns.

mainly leads to an increase in HD and HW, while an increase in the scanning speed leads to a decrease in HW and a decrease in HD.

For linear pattern, power contribution represents more than 75% in the HD variation while it varies between 45% and 58% for the rest of the patterns. The contribution of the scanning speed in the HW variation is between 73% and 86% for sinusoidal, triangular and trochoid patterns. However, it does not exceed 50% for linear pattern. The contribution of the error is minimal in most cases and does not exceed 10%. This points out the accuracy of the achieved experimental results and presumes good experimental planning and measurement.

These observations are confirmed by the graphs of the average effects presented in Figures 7-10. According to these figures, it is clear that the increase of laser power affects positively the hardened depth and width, while the increase of scanning speed affects negatively these two characteristics for the four scanning patterns. The magnitude of these effects varies considerably from one pattern to another. These results are confirmed by the correlation analysis between laser parameters and hardened zone dimensions for the four scanning patterns presented in Table 7. Strong correlation can be observed between power and HD and between scanning speed and HW.

Figure 7. Main effects of laser parameters on HD and HW for the linear moving pattern.

Figure 8. Main effects of laser parameters on HD and HW for sinusoidal moving pattern.

Figure 9. Main effects of laser parameters on HD and HW for triangular moving pattern.

Figure 10. Mains effects of laser parameters on HD and HW for trochoid moving pattern.

Table 7. Correlation analysis between laser parameters and hardened zone dimensions for the considered scanning patterns.

ANOVA is also used to develop linear multiple regression models in order to establish a simplified empirical relationship to estimate or predict the hardened depth and the hardened width as a function of the laser parameters for the four moving patterns. Table 8 presents the achieved regression models for both HD and HW that include the process parameters and their interactions. Curves shown in Figures 11-14 present comparisons between measured and predicted hardened depth and hardened width. The predicted HD and HW values follow the measured values with an overall error of less than 10%.

Applying the same experimental design for the four scanning patterns provides the possibility to evaluate their direct effects on the hardened depth and the hardened width. The experimental results reveal that linear scanning pattern achieves the maximal HD and the minimal HW compared to other scanning patterns. In the case of test 7 with a laser power of 800 W and a scanning speed of 20 mm/s, the linear pattern produces a hardened zone characterized by HD/HW of 1011/1154 µm, while the sinusoidal pattern produces a HD/HW of 456/2808 µm. Under the same conditions and for the same test, the triangular

Figure 11. Scatter plot—Measured and predicted HD and HW for linear scaning pattern.

Figure 12. Scatter plot—Measured and predicted HD and HW for sinusoidal scaning pattern.

Figure 13. Scatter plot—Measured and predicted HD and HW for triangular scaning pattern.

Table 8. Regression equations derived from ANOVA for the four moving patterns.

Figure 14. Scatter plot—Measured and predicted HD and HW for trochoid scaning pattern.

trajectory reaches an HD/HW of 574/2516 µm, while the HD/HW of the trochoid trajectory reaches 802/2578 µm. Figure 15 presents an overview of the effect of each moving pattern on HD and HW. Figure 16 presents the micro-hardness measurement results in the form of hardness curves produced for test 7 achieved using maximum power and minimum scanning speed for each scanning pattern. The hardness curves confirm the same tendencies observed in Figure 15. When using a sinusoidal scanning pattern, the laser beam spends much more time in the peripheral areas of the surface to be treated compared to the triangular or trochoid pattern, which concentrates more time and power in the central zone that increases the hardened depth. Although the trochoid pattern presents good result, it generates a non-uniform HD. This lack of uniformity can cause several problems in the case of thin steel parts surface heat treatment.

When the dimensions of the surface to be treated are larger than the cross section of the laser beam, some specific scanning patterns can be used to overcome the problem caused by the multi-track approach. The objective of choosing these scanning patterns is to maximize the surface of the hardened zone. This consists of maintaining a maximum uniform hardened width while respecting a minimum limit of hardened depth. In this perspective, if a minimum threshold is set for HD, the HW/HD ratio can be used as an indicator to evaluate the hardened zone surface according to the laser parameters for the different scanning patterns. HW/HD ratio also provides useful information on how heat is absorbed by the surface. A priori, if HW/HD increases, the hardened surface increases on condition that HD is kept close to HD limit. Table 9 presents the HW/HD ratio for various combinations of laser parameters and scanning patterns. The global average HW/HD ratio is, of course, higher than the average calculated excluding tests where HD is lower than 200 μm. In this table, it can be observed that sinusoidal and triangular scanning patterns present the highest HW/HD ratio. Therefore, these two patterns allow to achieve the largest hardened surfaces. With high HD and small HW, the linear pattern produces the low and relatively constant HW/HD ratio.

Although producing a higher ratio and therefore a higher HW than linear pattern, the trochoid pattern presents practically constant HW/HD ratio. These observations are confirmed by the contour plot presented in Figure 17. These

Figure 15. Comparison of HD and HW achieved using various scanning patterns.

Figure 16. Micro-hardness measurements achieved using the considered scanning patterns (test 7).

Table 9. HW/HD ratio for various combinations of laser parameters and scanning patterns.

Figure 17. Contour plot of HW/HD ratio vs P and SS for (a) linear; (b) sinusoidal; (c) triangular and (d) trochoid scanning pattern.

graphs reveal that HW/HD ratio increases with SS and decreases with P. This is typical because HW/HD tends to growth when HD decreases. If the laser parameters can be set to control a constant minimal HD close to 200 μm, by adopting the maximum ratio, HW can reach a minimum of 385, 3670, 3145 and 990 μm for linear, sinusoidal, triangular and trochoid scanning pattern, respectively. If the total generated heating energy is adjusted to maintain a constant minimal HD, HW can reach higher values when using sinusoidal or triangular patterns.

4. Conclusion

This paper presents an experimental investigation of laser surface hardening of AISI 4340 steel using different laser scanning patterns. The experimentations are carried out using a 3 kW Nd: YAG laser system and considering four laser source moving patterns: linear, sinusoidal, triangular and trochoid pattern. The effects of laser process parameters, such as laser power and scanning speed and their interactions on the hardened depth and hardened width, for various laser scanning patterns, are evaluated and analyzed using structured experimental design and improved statistical analysis tools. Graphs of main effects, percent contributions, F-test, regression analysis, correlation analysis and response surface methodology are programmed, implemented and discussed for each scanning pattern. The results show that the laser power positively affects the hardened depth more than the hardened width, while the scanning speed negatively affects the hardened depth but more significantly the hardened width. However, the parameter’s interactions have little effect on the hardness profiles. The results show also that the sinusoidal pattern gives a maximum width with an acceptable depth and a uniform hardness distribution, whereas the linear pattern gives a maximum depth with a width. The trochoid pattern presents a relatively good compromise between hardened depth and hardened width, but with an irregular hardness distribution and a non-uniform hardened depth. Hence, if a maximum of hardened width is required, it is recommended to use the sinusoidal pattern. For an extensive evaluation of the laser scanning pattern effects, it will be interesting to develop a numerical predictive model to conduct numerical investigations using a 3D thermal analysis.

Cite this paper
Tarchoun, B. , Ouafi, A. and Chebak, A. (2020) Experimental Investigation of Laser Surface Hardening of AISI 4340 Steel Using Different Laser Scanning Patterns. Journal of Minerals and Materials Characterization and Engineering, 8, 9-26. doi: 10.4236/jmmce.2020.82002.

[1]   Shercliff, H. and Ashby, M. (1991) The Prediction of Case Depth in Laser Transformation Hardening. Metallurgical Transactions A, 22, 2459-2466.

[2]   Selvan, J.S., Subramanian, K. and Nath, A.K. (1999) Effect of Laser Surface Hardening on En18 (AISI 5135) Steel. Journal of Materials Processing Technology, 91, 29-36.

[3]   Ashby, M. and Easterling, K.E. (1984) The Transformation Hardening of Steel Surfaces by Laser Beams—I. Hypo-Eutectoid Steels. Acta Metallurgica, 32, 1935-1948.

[4]   Leung, M.K.H., Man, H.C. and Yu, J.K. (2007) Theoretical and Experimental Studies on Laser Transformation Hardening of Steel by Customized Beam. International Journal of Heat and Mass Transfer, 50, 4600-4606.

[5]   Steen, W., Watkins, K.G. and Mazumber, J. (2010) Laser Material Processing. 4th Edition, Springer, London, xvii, 558 p.

[6]   Miokovic, T., Schulze, V., Vöhringer, O. and Löhe, D. (2006) Prediction of Phase Transformations During Laser Surface Hardening of AISI 4140 Including the Effects of Inhomogeneous Austenite Formation. Materials Science and Engineering, 435-436, 547-555.

[7]   Lakhkar, R.S., Shin, Y.C. and Krane, M.J.M. (2008) Predictive Modeling of Multi-track Laser Hardening of AISI 4140 Steel. Materials Science and Engineering, 480, 209-217.

[8]   Martínez, S., Lesyk, D., Lamikiz, A., Ukar, E. and Dzhemelinsky, V. (2016) Hardness Simulation of Over-Tempered Area during Laser Hardening Treatment. Physics Procedia, 83, 1357-1366.

[9]   Woodard, P.R. and Dryden, J. (1998) Thermal Analysis of a Laser Pulse for Discrete Spot Surface Transformation Hardening. Journal of Applied Physics, 85, 2488-2496.

[10]   Xue, L., Islam, M.U. and McGregor, G. (1999) Dot Matrix Hardening of Steels Using a Fiber Optic Coupled Pulsed Nd:YAG Laser. Materials and Manufacturing Processes, 14, 53-65.

[11]   Straffelini, G., Menapace, C. and Molinari, A. (2002) Interpretation of Effect of Matrix Hardening on Tensile and Impact Strength of Sintered Steels. Powder Metallurgy, 45, 167-172.

[12]   Farrahi, G.H. and Sistaninia, M. (2009) Thermal Analysis of Laser Hardening for Different Moving Patterns. International Journal of Engineering Transactions A: Basics, 22, 169-180.

[13]   Zhou, H. (2011) Temperature Rise induced by a Rotating or Dithering Laser Beam. Advanced Studies in Theoretical Physics, 5, 443-468.

[14]   Komanduri, R. and Hou, Z.B. (2001) Thermal Analysis of the Laser Surface Transformation Hardening Process. International Journal of Heat and Mass Transfer, 44, 2845-2862.

[15]   Komanduri, R. and Hou, Z.B. (2004) Thermal Analysis of Laser Surface Transformation Hardening—Optimization of Process Parameters. International Journal of Machine Tools and Manufacture, 44, 991-1008.

[16]   Sistaninia, M. (2009) Laser Heat Treatment of Metals Using Rotating and Dithering Beams. In: Recent Advances in Technologies, IntechOpen.

[17]   Tan, T. and Zhou, H. (2013) Numerical Solution for a Transient Temperature Distribution on a Finite Domain Due to a Dithering or Rotating Laser Beam. International Journal of Operations Research and Information Systems, 4, 38.

[18]   Philip, J.R. (1995) Taguchi Techniques for Quality Engineering. McGraw-Hill, New York, 249.