Experimental Investigation of Thermal Fatigue Die Casting Dies by Using Response Surface Modelling

Mechanical and thermal sequences impact largely on thermo-mechanical fatigue of dies in a die casting operations. Innovative techniques to optimize the thermo-mechanical conditions of samples are major focus of researchers. This study investigates the typical thermal fatigue in die steel. Die surface initiation and crack propagation were stimulated by thermal and hardness gradients, acting on the contact surface layer. A design of experiments (DOE) was developed to analyze the effect of as-machined surface roughness and die casting parameters on thermal fatigue properties. The experimental data were assessed on a thermo-mechanical fatigue life assessment model, being assisted by response surface methodology (RSM). The eminent valuation was grounded on the crack length, hardness properties and surface roughness due to thermal fatigue. The results were analyzed using analysis of variance method. Parameter optimization was conducted using response surface methodology (RSM). Based on the model, the optimal results of 26.5 μm crack length, 3.114 μm surface roughness, and 306 HV0.5 hardness properties were produced.


Introduction
Reducing process lead time and design time are important aspect to reduce total cost of die casting process.Minimizing the trial-and-error stage of the production can further assist in reducing the cost.The shape of the cavity and die geometry is directly related to soldering tendency of a die casting.The more complicated geometry a die has, the more likely soldering can occur.The dies of complicated geometry usually procure sharp angles, core pins, nooks, and the part that acts as hot spots.These hot spots can induce soldering, due to their higher temperatures than other areas [1].Moreover, prolonging die service life and preventing catastrophic die failures is essential because the die cost contributes considerably to the overall process cost.Die-casting dies are subjected to high mechanical and thermal loads, which can cause appalling or delayed die failure due to mount up damages.Most failures developed gradually and can be predicted.However, ongoing failures after few cycles can root for a larger economical losses due to wrecked tools, expensive down times, and disordered delivery schedules [2].
Thermal fatigue is one of the most common complications encountered in a die casting process.It is a result of the cyclic, rapid, non-uniform heating and cooling of dies.Thermal fatigue gets amplified by mechanical loadings and other damage mechanisms, such as erosion and corrosion.Thermal fatigue is one of the major foundations of poor-quality castings and die failures.The initial phase of the thermal fatigue initiates with the formation of micro-crack networks denoted as heat checking.Heat checking results in the deterioration of die cavity surfaces.Further damage is due to heat checking leads to clasping, defective castings, and subsequently die damage and failure [3][4][5][6][7][8].When an unconstrained object is slowly and uniformly heated or cooled, it will expand or contract proportionally to its coefficient of thermal expansion, and a temperature change will occur in relation to the mentioned temperature.However, if the body is constrained during the heating and cooling process, it will develop stresses and strains because it cannot reach its unconstrained dimensions properly.Similarly, if a body does not load and release through the same stress-strain path, repeated or cyclic heating and cooling might outcome in the build-up of inelastic deformation.If the time-varying stress within the body is tensile, accumulation may lead to a thermal fatigue.Even an unconstrained body may develop large thermal stresses if exposed to non-uniform or rapid heating and cooling.
Objects with sudden geometric and compositional variations retain internal constraints which are caused by temperature distribution.In a case of rapid surface temperature changes in a solid body, the temperature of a nearby surface layer increases and decreases promptly, whereas the rest of the body cannot react fast enough to the alterations in temperature.Different expansions and contractions within the body provide internal constraints, and various layers expand and contract at different degrees.Each region gets constrained by its neighbouring region, which results in a stress-strain field within the body [7][8][9].Thermally induced stresses are not the only factors involved in the thermal fatigue phenomenon.Die material properties at elevated temperatures, mechanical loads from filling and locking, residual stresses, and cavity surface conditions also affect a material's response to heat checking and a thermal fatigue.Tempering, decarburization, and phase or structure changes due to an exposure to elevated temperatures accelerate thermal fatigue.
To predict and avoid die failures and their consequences, thermal management of dies, the interactions between thermal and stress/strain fields, material properties, and the influence of all these factors on die service life must be understood.By investigating the effects of these factors on die life during the design stage of dies and the casting process, prolonged die service life can be succeeded.A reliable die service life prediction allows for accurate estimates of actual die and production cost, reduces undesired machine down times, and helps to achieve good production management during production processes [2,[9][10][11].Furthermore, thermo-mechanical fatigue (TMF) experiments on die casting require costly equipment, time consuming, and often conducted through thermal testing under the same operating conditions as that for die-casting dies; however, in fact, the temperature and number of cycles are fixed and kept constant [9][10][11][12].The damage on the die obtained by TMF experiments is less than that obtained in isothermal fatigue (IF) experiments [13].The purpose of this study is to fit a model of the die life by using an experimental data.Experimental procedures were designed to mimic the commercial die casting process.The thermal fatigue process was conducted in a cyclic manner with constant sample temperature during the test [14][15][16].Analysis of dies was performed by associating several of the previously mentioned factors that influence the temperature profile and the thermal gradients with the structural state, within the tool.The step was implemented by employing the temperature profiles obtained through a previous 1D thermal analysis and imposing them on the die cavity/casting interface via 2D thermal finite element models (FEMs), respectively.Thermal finite element analysis was conducted using ABAQUS software (Version 6.13.1, ABAQUS, Johnston, RI, USA, 2013), and this analysis was followed by a sequentially coupled structural analysis.The cooling effects of the lubricants were studied through experimental and numerical investigations.The effects of initial die cavity surface temperature and spray fluid density (spray volume per-unit surface area per-unit time as defined by Lee), have been a focus of several previous researchers [3][4][5][6][17][18][19][20].
In this study the effects on service life were investigated through response surface methodology (RSM), employing a three-factor Box-Behnken design.The thermal fatigue cycle parameters were evaluated with three levels for each factor.An experimental investigation was conducted to quantify the relationship between the input factors (as-machined surface roughness, R a1 , wall thickness, and immersion time) and responses (crack length, surface roughness due to thermal fatigue, R a2 , and hardness properties).Moreover, by combining existing and new data on heat transfer in die-casting dies and die material properties with RSM modelling, a method for predicting the onset of heat checking in H13 steel dies has been proposed.

Materials
The samples in this study were made from H13 tool steel .In the fatigue testing, the sample was dipped into molten A356 aluminium alloy was used in the study.Chemical configuration for both H13 tool steel and A356 aluminium alloy are shown in Tables 1 and 2 respectively.The geometry of the samples is shown in Figure 1 with different sample wall thicknesses to analyze thermal gradient, determine efficient cooling in the central part and thermal fatigue cracks and varied surface roughness due to thermal fatigue.the relationship between the input factors (as-machined surface roughness, Ra1, wall thickness, and immersion time) and responses (crack length, surface roughness due to thermal fatigue, Ra2, and hardness properties).Moreover, by combining existing and new data on heat transfer in die-casting dies and die material properties with RSM modelling, a method for predicting the onset of heat checking in H13 steel dies has been proposed.

Materials
The samples in this study were made from H13 tool steel .In the fatigue testing, the sample was dipped into molten A356 aluminium alloy was used in the study.Chemical configuration for both H13 tool steel and A356 aluminium alloy are shown in Tables 1 and 2 respectively.The geometry of the samples is shown in Figure 1 with different sample wall thicknesses to analyze thermal gradient, determine efficient cooling in the central part and thermal fatigue cracks and varied surface roughness due to thermal fatigue.

Experimental Apparatus and Work Procedure
A simulative laboratory experiment was conducted to generate controlled thermo-mechanical fatigue failures.Thermal wear set-up simulates an aluminium alloy die-casting conditions with testing cycle time of 24 s, 28 s and 32 s.The samples were dipped into molten aluminium at 700 °C and quenched in water at 32 °C for 7 s, 9 s and 11 s to generate thermal gradients during cooling process.In between the heating and water quenching, the samples were air-cooled at 28 °C for 5 s before repeating the next cycle.Thermal fatigue loading was achieved by 1850 cyclic movements.The inner wall of the samples was cooled by spraying water at room temperature for 3 s.

Experimental Apparatus and Work Procedure
A simulative laboratory experiment was conducted to generate controlled thermo-mechanical fatigue failures.Thermal wear set-up simulates an aluminium alloy die-casting conditions with testing cycle time of 24 s, 28 s and 32 s.The samples were dipped into molten aluminium at 700 • C and quenched in water at 32 • C for 7 s, 9 s and 11 s to generate thermal gradients during cooling process.In between the heating and water quenching, the samples were air-cooled at 28 • C for 5 s before repeating the next cycle.Thermal fatigue loading was achieved by 1850 cyclic movements.The inner wall of the samples was cooled by spraying water at room temperature for 3 s.
Samples surfaces were examined through optical microscope (MT8000 Series Metallurgical Microscopes, Microscope.com,Roanoke, VA, USA) and scanning electron microscope (SEM) (HITACHI Tabletop Microscope TM3030Plus, Hitachi High Technologies America, Inc., Schaumburg, IL, USA).For each immersion experiment scheme cycle, a characteristic testing sample from the final treated sample (when the sample was clamped on the equipment), was obtained.The transverse parts of the cycled samples were investigated for the exclusive evaluation of the detrimental penetration depth, microhardness, and microstructural amendments by using (SEM) and Vickers hardness tester (Instran ITW Test & Measurement Co., Ltd., Shanghai, China), respectively.The sample surface profile was measured with a two-dimensional surface profilometer MarSurf PS1 (Mahr, Göttingen, Germany).Surface roughness were measured using MAHR parameter MarSurf PS1 (Mahr, Göttingen, Germany), before and after the process.Thermally worn samples were investigated for crack analysis, metallographic study, and hardness properties.Metallographic study and crack analysis were performed along a longitudinal cross-section through SEM and Energy-dispersive X-ray spectroscopy (EDXS) analysis.The crack length on the sample surface was measured with ImageJ software (Version 2, National Institutes of Health, Bethesda, MD, USA, 2008); an image processing software available online [21].The crack tips were then marked manually.The hardness properties were measured across the sample thickness with Vickers indenter.

Factors and Response Surface Modeling
Among the factors that may affect the response, several were selected to be kept constant during the experiment.The selected factors included are material, cooling fluid, and the sample heat treatment.In addition, several important factors related to thermal fatigue were considered.These factors included as-machined surface roughness, sample wall thickness, and immersion time cycle, which is a parameter that depends on the amount of cracks and hardness induced in the material (Table 3).The maximum crack length, gradient hardness, and roughness after the thermal fatigue test were observed as responses in the model.Thermal-mechanical and temperature cycles are phase variations, which means that the minimum load coincides with the minimum temperature, and the maximum tensile load is applied with the maximum temperature [22,23].The investigated material was the tool steel used in the thermal fatigue wear test (H13 tool steel).The experimental plan and results for the crack length (CL s ) of samples, surface roughness due to thermal fatigue (R a2 ) and hardness properties obtained from the experimentation are presented in Table 4.

Effect of Surface Roughness, Wall Thickness and Immersion Time on Crack Spherical Shape (SPH)
The second-order polynomial model was used to approximate the relationship between the three factors and response crack with Equation ( 1) suitably.The model showed that crack increased with an increase in the immersion time and decreased with an increase in the wall thickness.Another observation indicated that these factors (as-machined surface roughness, wall thickness, and immersion time) were "significant" to crack, respectively [11].Table 5 shows the analysis of variance indicating that the model is adequate because the p-value of the square is much more than linear.
CL s = +2.3622+ 0.41111 × A + 0.086558B + 0.17442C + 7.7737 × 10 −4 AB + 6.1881 × 10 −3 AC − 5.1788 × 10 −3 BC, (1) where CL s is crack length, A is as-machined surface roughness, B is wall thickness, and C is immersion time.Crack SPH was obtained experimentally, and other values were predicted, in which the estimated regression coefficient for the second order predicted the crack and analysis of variance, as shown in Table 6.The small p-values for the linear term showed that their contribution was significant in the model.Moreover, the main effects were individually significant at the 0.05 significance level.Crack SPH was significant to the response model at α = 0.05.From the value of R 1 (96.17%), the fit of data can be measured from the estimated model. Figure 2 shows that the residuals follow a straight line, and the errors appear to be normally distributed.In such case, many parameters were carefully estimated to justify the test.Figure 3 shows that the response surface model determined above is relevant to the crack surface.The factors that are not on the plot are at their average level.The plots exhibited a square according to the surface 2FI model, and immersion time and R a1 contribute equally to reducing crack.Meanwhile, reduced crack length extends the lifetime of dies, but this factor cannot compensate for the reverse effect due to an increase in the other factors; wall thickness and R a1 .

Effect of As-Machined Surface Roughness, Wall Thickness and Immersion Time on R a2
To describe the effective factors, a second-order equation was established under the conditions of surface roughness due to thermal fatigue, R a2 , is expressed as follows: where R a2 is the surface roughness due to thermal fatigue.Table 7 indicates that the model is adequate because the p-value of the square is much more than linear.To assess the validity of the model (Equation ( 2)), a probability plot (Figure 6) was used to compare the measured and predicted y 2 .The R a1 , wall thickness, and immersion time data are plotted in Figure 4. Figure 4 show that most value matched one another well, except that the difference between the measured and predicted values exceeded 0.6 µm.The smallest value of R a2 was measured for H13 tool steel after thermal fatigue cycles (1850 cycles) with a constant temperature of 700 • C. The data on additional factors were used to generate the probability plot, except for two data points where the model overpredicted the measured data by over 0.6 µm.The additional test data fit the model reasonably and supports the validity of the model in predicting the values of R a2 .

Effect of As-Machined Surface Roughness, Wall Thickness and Immersion Time on R a2
To describe the effective factors, a second-order equation was established under the conditions of surface roughness due to thermal fatigue, Ra2, is expressed as follows: R a2 = 1.191 + 0.249A + 0.027B -4.274 × 10 -3 AB + 4.564 × 10 -3 AC -2.291 × 10 -3 BC, (2) where R a2 is the surface roughness due to thermal fatigue.Table 7 indicates that the model is adequate because the p-value of the square is much more than linear.To assess the validity of the model (Equation ( 2)), a probability plot (Figure 6) was used to compare the measured and predicted y2.The Ra1, wall thickness, and immersion time data are plotted in Figure 4. Figure 4 show that most value matched one another well, except that the difference between the measured and predicted values exceeded 0.6 μm.The smallest value of Ra2 was measured for H13 tool steel after thermal fatigue cycles (1850 cycles) with a constant temperature of 700 °C.The data on additional factors were used to generate the probability plot, except for two data points where the model overpredicted the measured data by over 0.6 μm.The additional test data fit the model reasonably and supports the validity of the model in predicting the values of Ra2.  Figure 5 shows a contour plot of Ra2, responding to three parameters.The increase in Ra1, and immersion time affected Ra2 dramatically.In Figure 6, the Ra2 reached the highest value when the Ra1 increases.
The response surface plots (i.e., graphing of Equation ( 2)) for Ra2 against Ra1 and wall thickness ratio, and for Ra2 against Ra1 and immersion time are shown in Figure 5a-d respectively.In Figure 5b, Figure 5 shows a contour plot of R a2 , responding to three parameters.The increase in R a1 , and immersion time affected R a2 dramatically.In Figure 6, the R a2 reached the highest value when the R a1 increases.The response surface plots (i.e., graphing of Equation ( 2)) for R a2 against R a1 and wall thickness ratio, and for R a2 against R a1 and immersion time are shown in Figure 5a-d respectively.In Figure 5b, the measured R a2 increased with R a1 and immersion time [24,25].Generally, the predicted R a2 for the different samples was lower than that for steel B and C, which is consistent with the measured results (Figure 5a,b).The effect of the difference in time and wall thickness on y 2 varied for the different H13 samples (Figure 5b,d).With R a1 versus immersion time increasing from 0% to 24%, R a2 increased for parameters A and C. With a decrease in R a2 , the effect on tool steel also decreased.the measured Ra2 increased with Ra1 and immersion time [24,25].Generally, the predicted Ra2 for the different samples was lower than that for steel B and C, which is consistent with the measured results (Figure 5a,b).The effect of the difference in time and wall thickness on y2 varied for the different H13 samples (Figure 5b,d).With Ra1 versus immersion time increasing from 0% to 24%, R a2 increased for parameters A and C. With a decrease in R a2 , the effect on tool steel also decreased.

Effect of As-Machined Surface Roughness, Wall Thickness and Immersion Time on Hardness Properties
The second-order equation was established to describe the influencing factors and conditions investigated in this study on the surface hardness.The second-order model can be expressed as where y3 is hardness.The model (Equation ( 3)) shows that surface hardness increased with an increase in surface roughness and wall thickness [26,27].Furthermore, the square of surface roughness, Ra1, and immersion time provided a good indication that as-machined surface roughness was a major factor

Effect of As-Machined Surface Roughness, Wall Thickness and Immersion Time on Hardness Properties
The second-order equation was established to describe the influencing factors and conditions investigated in this study on the surface hardness.The second-order model can be expressed as where y 3 is hardness.
The model (Equation ( 3)) shows that surface hardness increased with an increase in surface roughness and wall thickness [26,27].Furthermore, the square of surface roughness, R a1 , and immersion time provided a good indication that as-machined surface roughness was a major factor Metals 2017, 7, 191 9 of 13 in R a2 changes.The analysis of variance as shown in Table 8 indicates that the model was adequate because the p-value of the response surface quadratic model is significant.The hardness values for all the samples at different through-the-wall thickness increased slightly with the increasing temperature.However, compared with hardness due to thermal fatigue cycles, the HV 0.5 values were much lower, when the R a1 and immersion time increased.These tests require the error term to be normally and independently distributed with mean zero and variances.Figure 6 shows the normal probability, fitted values, and histogram of residuals, respectively.
Metals 2017, 7, 191 9 of 13 in Ra2 changes.The analysis of variance as shown in Table 8 indicates that the model was adequate because the p-value of the response surface quadratic model is significant.The hardness values for all the samples at different through-the-wall thickness increased slightly with the increasing temperature.However, compared with hardness due to thermal fatigue cycles, the HV0.5 values were much lower, when the Ra1 and immersion time increased.These tests require the error term to be normally and independently distributed with mean zero and variances.Figure 6 shows the normal probability, fitted values, and histogram of residuals, respectively.

Influence of Cooling Rate on Hardness Properties and Crack Length
In the experiment, the maximum and the variety of temperatures reached at the surface and the difference inside the sample was investigated only as a function of the cooling system.A cooling system brought the maximum temperature at the surface.Therefore, the effect of cooling system on crack length as a function of maximum temperature at the surface of the sample.The value assessments of maximum crack length decreased with the increase in the maximum temperature.The identical trend was observed in the experiment conducted for the thermal fatigue cracking parameters.The relationship between the cracking and microhardness measured at a distance equal to the average maximum crack length, followed the temperature trend and confirmed the observation made for different immersion times, as shown in Figure 8.The longer crack is mainly because of the higher temperature at the surface, the lower is the hardness next to crack.

Influence of Cooling Rate on Hardness Properties and Crack Length
In the experiment, the maximum and the variety of temperatures reached at the surface and the difference inside the sample was investigated only as a function of the cooling system.A cooling system brought the maximum temperature at the surface.Therefore, the effect of cooling system on crack length as a function of maximum temperature at the surface of the sample.The value assessments of maximum crack length decreased with the increase in the maximum temperature.The identical trend was observed in the experiment conducted for the thermal fatigue cracking parameters.The relationship between the cracking and microhardness measured at a distance equal to the average maximum crack length, followed the temperature trend and confirmed the observation made for different immersion times, as shown in Figure 8.The longer crack is mainly because of the higher temperature at the surface, the lower is the hardness next to crack.

Conclusions
In the following research study, a mathematical model with second order has been employed to anticipate the crack parameters including the immersion time, surface roughness, Ra1, and wall thickness process by evaluating the surface roughness due to thermal fatigue, Ra2, and the temperature respectively, which is based on response surface methodology (RSM).The mathematical model outcomes were later equated with the experimental findings.It was observed that the immersion time, Ra1, and wall thickness affected Ra2 and temperature distribution when milling the H13 tool steel.Some important conclusions derived from the study are concluded as following: 1.
The response surface methodology (RSM) resulted in an advantageous procedure for the surface roughness and temperature analysis.In addition, designing experiments are essential to produce statistics, which in turn is beneficial in expanding the calculating equations for surface roughness, crack length, and hardness properties.The investigation of variance for the second order for both the studied model displays that the immersion time is the most affected parameter which afflicted the hardness and the crack lengths followed by wall thickness.

2.
Both second order models were observed to be expedient in forecasting the key effects and the square effects of diverse dominant arrangements of the machining constraints.The process was found cost-effective in shaping the effect of several parameters in a methodical way.In addition, the process for the thermal fatigue cycle of H13 tool steel, the legitimacy of the process is typically restricted to the collection of factors measured during the investigation.

3.
The RSM model could effectively relate the machining parameters with the responses, crack, surface roughness due to thermal fatigue, and hardness properties.The optimal parameter setting resulted crack length of 26.5 μm, surface roughness of 3.114 μm, and hardness properties of 306 HV0.5.

4.
The results generated by the predicted model are equated with the experimental results.The observed percentage error is very low, which is only 2% for both the predicted models.

Conclusions
In the following research study, a mathematical model with second order has been employed to anticipate the crack parameters including the immersion time, surface roughness, R a1 , and wall thickness process by evaluating the surface roughness due to thermal fatigue, R a2 , and the temperature respectively, which is based on response surface methodology (RSM).The mathematical model outcomes were later equated with the experimental findings.It was observed that the immersion time, R a1 , and wall thickness affected R a2 and temperature distribution when milling the H13 tool steel.Some important conclusions derived from the study are concluded as following: 1.
The response surface methodology (RSM) resulted in an advantageous procedure for the surface roughness and temperature analysis.In addition, designing experiments are essential to produce statistics, which in turn is beneficial in expanding the calculating equations for surface roughness, crack length, and hardness properties.The investigation of variance for the second order for both the studied model displays that the immersion time is the most affected parameter which afflicted the hardness and the crack lengths followed by wall thickness.

2.
Both second order models were observed to be expedient in forecasting the key effects and the square effects of diverse dominant arrangements of the machining constraints.The process was found cost-effective in shaping the effect of several parameters in a methodical way.In addition, the process for the thermal fatigue cycle of H13 tool steel, the legitimacy of the process is typically restricted to the collection of factors measured during the investigation.

3.
The RSM model could effectively relate the machining parameters with the responses, crack, surface roughness due to thermal fatigue, and hardness properties.The optimal parameter setting resulted crack length of 26.5 µm, surface roughness of 3.114 µm, and hardness properties of 306 HV 0.5 .

4.
The results generated by the predicted model are equated with the experimental results.The observed percentage error is very low, which is only 2% for both the predicted models.

Figure 3 .
Figure 3. Contour plot of crack length responding to (a) Ra1 and wall thickness, (b) Ra1 and immersion time.(c,d) a 3D view of crack length interaction with the respective parameters.

Figure 3 .
Figure 3. Contour plot of crack length responding to (a) Ra1 and wall thickness, (b) Ra1 and immersion time.(c,d) a 3D view of crack length interaction with the respective parameters.

Figure 3 .
Figure 3. Contour plot of crack length responding to (a) R a1 and wall thickness, (b) R a1 and immersion time.(c,d) a 3D view of crack length interaction with the respective parameters.

Figure 5 .
Figure 5. Contour plot of surface roughness due to thermal fatigue, Ra2, responding to (a) Ra1 and immersion time, (b) Ra1 and wall thickness.(c,d) a 3D view of Ra2 interactions with the respective parameters.

Figure 5 .
Figure 5. Contour plot of surface roughness due to thermal fatigue, R a2 , responding to (a) R a1 and immersion time, (b) R a1 and wall thickness.(c,d) a 3D view of R a2 interactions with the respective parameters.

Figure 7
Figure 7 shows the hardness contours at three different parameters.The increase in Ra1 and immersion time clearly affected Ra2 dramatically.In Figure 5, Ra2 reached the highest value, when the Ra1 increases.Hardness properties increased more rapidly with increasing temperature in comparison to increasing time.All the empirical equations predicting change in the properties are based on the three factors, namely, as-machined surface roughness, wall thickness, and immersion time.

Figure 7
Figure7shows the hardness contours at three different parameters.The increase in R a1 and immersion time clearly affected R a2 dramatically.In Figure5, R a2 reached the highest value, when the R a1 increases.Hardness properties increased more rapidly with increasing temperature in comparison to increasing time.All the empirical equations predicting change in the properties are based on the three factors, namely, as-machined surface roughness, wall thickness, and immersion time.

Figure. 7 .
Figure.7.Contour plots of hardness properties responding to (a) Ra1 and immersion time, (b) Ra1 and wall thickness.(c,d) a 3D view of hardness properties, interaction with the respective parameters.

Figure 7 .
Figure 7. Contour plots of hardness properties responding to (a) R a1 and immersion time, (b) R a1 and wall thickness.(c,d) a 3D view of hardness properties, interaction with the respective parameters.

Table 1 .
Chemical composition of American Iron and Steel Institute (AISI) H13.

Table 2 .
Chemical composition of Al 356.

Table 1 .
Chemical composition of American Iron and Steel Institute (AISI) H13.

Table 2 .
Chemical composition of Al 356.

Table 3 .
Surface design factor levels.

Table 4 .
Outcome of defect thermal fatigue cycles.

Table 5 .
Response surface 2FI model for crack length.

Table 6 .
Analysis of variance for crack.

Table 7 .
Estimated regression coefficients for surface roughness due to thermal fatigue, R a2 .

Table 7 .
Estimated regression coefficients for surface roughness due to thermal fatigue, Ra2.

Table 8 .
Response Surface Quadratic model by the analysis of variance (ANOVA).

Table 8 .
Response Surface Quadratic model by the analysis of variance (ANOVA).