Mechanics, Fatigue, and Fracture of Structural JointsView this Special Issue
Mechanical Properties Optimization of Friction Stir Welded Lap Joints in Aluminium Alloy
This work focuses on the influence of the rotational and travel speed on the strength of AA 2024 T3 friction stir welded lap joints. Tensile tests were carried out on minispecimens extracted from different welding zones. A central composite design was applied to identify the relative importance of the variable factors’ effects and their interaction on yield/ultimate strength and elongation for both the heat affected-zone (HAZ) and nugget zone. Surface methods and gradient algorithms were used to optimize the yield strength of the joints. Shear and microhardness tests were executed to achieve a more complete mechanical characterization.
The feasibility of replacing the riveting process with friction stir welding (FSW) technology in the assembly of fuselage skin-stiffener panels was the objective of several scientific papers in the last decades [1–6]. The potential of FSW to result in significantly lower assembly times and manufacturing costs, yet higher productivity, cannot be overemphasized. However, commercial applications of this process require addressing issues such as strength analysis and design as well as optimizing the manufacturing process parameters for more reliable welds with minimum distortion [7, 8]. One of the many experimental strategies used to determine the process parameters as well as to optimize the process is defined “by trial” and relies on the researcher’s technical and theoretical process knowledge . This strategy consists of a sequential procedure wherein the value of just one parameter among the set of parameters identified in the previous experiment is varied. This specific parameter is the one that mostly influences the process response while keeping fixed the level of the other parameters . Such an approach has at least two disadvantages: it lacks objectivity in determining mathematical relations between process inputs and response variables and it also lacks efficiency and accuracy in optimizing the surface response through trials . Another experimental strategy is the “one factor per time approach,” wherein the researcher subsequently varies the levels of each factor within its range while keeping fixed at the base level the other factors . Obviously, this approach lacks information about the factors’ interactions, but it enables estimating main factors’ effects. In any event, the one factor per time approach is less accurate and less efficient than the statistical approach . Many studies of experimental strategies prove that the correct experimental strategy for problems involving more than one factor and data points affected by experimental errors must be taken from design of experiments (DOE) techniques if the aim is to remove objective results from experimental data . These techniques allow maximum conclusive information to be drawn from minimum work, time, energy, money, or other limited resources. So, DOE aims to maximize information per run by choosing a reduced number of input sample points. Moreover, it enables developing a robust process, i.e., a process that minimizes the influence of noise variables . Data organized by DOE allow the most powerful use of analysis of variances, since main factor and interaction effects can be estimated to determine their significance. Another considerable advantage of using these techniques involves combining DOE with response surface methodologies (RSMs), more efficiently to evaluate the effectiveness of investigative trade-offs . One of the most important advantages of using DOE techniques is extrapolation . This property is significant because it allows optimum prediction even if “optimum” lies outside the initial design range, and the number of runs must be reduced to achieve optimum. Good predictive properties allow DOE association with other complex mathematical methods without compromising accuracy and assuring the method of steepest ascent. This method is a viable technique for sequentially moving toward the optimum response . A comprehensive recitation of additional advantages of these DOE techniques is beyond the scope of this paper. Several papers deal with the influence of process parameters on properties of FSW joints [19–22]; in particular, it was adopted a DOE to investigate such influence . From the literature, it is also evident that the main process parameters that need to be investigated are the rotational speed and travel speed [24–26]. In particular, it was proved that the inverse of travel speed measurements was better than travel speed to fit experimental data . Therefore, this study examined the strength of FSW aluminium alloy lap joints in two different welding zones: HAZ and nugget zone. Aiming to test separately the tensile strength of HAZ and nugget zone, minispecimens from both these zones were subjected to shear and microhardness tests. The process parameters optimization issue was approached using central factorial design, response surface methods, and gradient algorithms.
2. Materials and Methods
AA 2024 T3 sheets were used as the bottom and the top sheet of the lap joints. The nominal weight percent composition (major alloying additions) of AA 2024 is 4.4% Cu, 1.5% Mg, 0.6% Mn, and the rest is aluminium. Sheets’ thickness is 1.27 mm. The ultimate tensile strength (UTS) of the AA 2024 T3 base metal is 475 MPa. As for the welding process parameters, the inverse of travel speed (TS−1) ranged from 0.155 to 0.533 s/mm and the rotational speed (RS) ranged from 950 to 2222.5 rpm. Plunge depth was fixed to 2.11 mm, tilt angle to zero degree, and travel angle to one degree.
Nine different TS−1-RS combinations were used for obtaining 13 welds in this study; five combinations were derived by 22 factorial design with five centre points, and the remaining four combinations were obtained by using the steepest ascent algorithm. One tool configuration equipped with a pin consisting of a threaded frustum of the cone was used; the lower and the upper diameters are 2.77 mm and 5.06 mm, respectively, with a length of 2.03 mm. The shoulder diameter is 12.04 with a concavity of 7°. Two different machines were used to perform FSW lap and overlap joints with single pass. A controlled numerical machine (CNM) was employed to perform the set of runs concerned with the central factorial design featured by high heat generation (hot runs); a vertical milling machine was used to carry out the set of runs concerned with the steepest algorithm featured by low heat generation (cold runs).
Both the top and bottom sheets, 152.40 mm long and 2.54 mm thick, were positioned as shown in Figure 1. Sheets were degreased prior to welding using acetone as a cleaner. After the FSW process, an optical specimen was taken out from each run; microscope observations of the cross sections were carried out to identify the exact position of the welding zones. At the same time, the cross sections of the welds were observed and analysed by image processing software to identify the position of the hook defects and measure the grain size. According to ASTM E 112, the grain size number is determined by using the general intercept method. For this purposes, standard metallographic polishing procedures were used with modified Keller’s reagent. Thanks to the aforementioned analyses, it was possible to remove three minitensile specimens from both the nugget and heat-affected zones for each run by using a minimilling machine.
Minitensile testing was carried out by following design of experiments. Figure 2 shows the minitensile specimen drawing with its geometrical dimensions. It is worth to note that smaller is the specimen higher are the mechanical properties recorded; experimental results of minitensile samples refer to local properties of the material which exceed the global ones . Each sample was pulled out from the lap joints; then, it was refined and polished on both sides by using abrasive papers of 30 µm and 12 µm removing exceeding material and leading the thickness from 1 mm to 0.50 mm. Abrasive papers of 9 µm and 3 µm were used to eliminate the surface scratches; finally, a diamond suspension of 1 µm was used to create a like-mirror surface. The waiting time between the welding and the minitensile testing was typically 150 h. The minitensile testing was performed through a minitensile testing machine by imposing a strain rate equal to 10−3 s−1. The ultimate strength, the yield strength, and the elongation of the weld were measured. The last property was evaluated by using the Epsilon ONE optical high-precision extensometer, which is a noncontact device able to measure accurately the strain of this type of minitensile specimen at the narrow section.
The mechanical characterization of the lap joints obtained by FSW technology also included shear testing and microhardness measurements. The strength of the lap joints loaded nominally in overlap shear was examined, even though in this case, 150 h elapsed before shear testing was carried out. All the specimens tested were 25.4 mm wide, 127 mm long, and 2.54 mm thick (Figure 3), manufactured using a shear machine, and ground on both sides by milling. Three specimens for each run were tested, and the failure loads were averaged. All overlap shear tests were performed using the same strain rate used in the tensile testing. The maximum loads (failure loads) were recorded for each specimen. Each weld was sectioned and polished so the hardness could be measured using a digital Vickers hardness tester. The centre of the welds was tested throughout the cross section with 1 mm spacing between data points; all welding zones down to the base material were involved.
As for the experiments design, the central composite design (CCD), able to fit the response surface through a quadratic regression model, was used in this study. CCD contains an embedded factorial or fractional factorial design with centre points that are augmented with a group of “star points.” The importance of the centre points is that they allow for a curvature estimation . The CCD could be implemented more easily with a sequential procedure: firstly, a factorial design with centre points should be performed to estimate the significance of the main factor effects and interaction effects, and above all to prove or to reject the presence of a curvature. Then, if the presence of the curvature is proved, the “star points” procedure can be performed. According to the aforementioned procedure, a 22 factorial design with five centre points and three repeated measurements was carried out to perform tensile testing for the nugget zone and for the heat-affected zone. The experimental table is summarized (Table 1).
To avoid systematic errors, the whole set of experiments has been randomized. To perform an ANOVA in every possible case, even in the case of heteroskedasticity of data (i.e., not constant variance), a weighted ANOVA was carried out . Such analysis requires repeated measurements for each run. Experimental data of yield strength, ultimate strength, and elongation are analysed and plotted separately. MINITAB, one of the most common statistical software packages, was used to define the DOE table and to perform a weighted ANOVA as well as residuals analysis. Experiments included repeats of the response, whose measurements allowed variability analysis in response data through a weighted regression. Such a method is suitable for handling data with different variances.
3. Results and Discussion
Experimental data of yield strength are shown in Table 2. An important point is that the column called WGTS_ADJ contains weights, i.e., an estimate of the reciprocal variance of the mean. These weights are due to the analysis of variability that is shown in Table 3. This analysis calculates and stores the standard deviations of the repeated responses detecting differences or dispersion effects across factor settings. This analysis even allows weights to be stored to perform a weighted ANOVA.
This analysis shows that only two effects could be considered significant: rpm with a value of 0.094 and the inverse of the travel speed with a value of 0.055; the interaction effect seems not to be noticeable. Neither curvature could be considered significant. The weighted ANOVA related to the experimental data of yield strength is shown in Table 4. The regression model has been reduced by removing the interaction effect since a previous analysis proved its negligibility; in fact, the value of F related to the interaction effect is 0.09 corresponding to a value of 0.782.
The analysis shows a curvature value of 0.145 that cannot be able to prove its presence and a main effect value of 0.040 that ensures a significant influence; specifically, the value related to rpm is 0.0045, and the value related to TS−1 is 0.0055. Both the values can be considered significant. Regression model coefficients are estimated as shown in Table 4. Aiming to prove the model adequacy, residuals analysis was performed that is not here reported for the sake of brevity. Such an analysis verified the hypothesis made on the residuals; i.e., residuals can be assumed to be normally distributed with null average and constant standard deviation. Moreover, the regression model fits the data very well since the high value related to the lack of fit (0.782) is just a consequence of the interaction term omission. Experimental data of ultimate strength are shown in Table 5. The analysis of variability related to the experimental data of ultimate strength is shown in Table 6.
It can be seen that only two effects could be considered significant: rpm with a value of 0.059 and the inverse of the travel speed with a value of 0.006; the interaction effect seems not to be noticeable. Either curvature could be considered significant. The weighted ANOVA related to the experimental data of ultimate strength is shown in Table 7. The regression model presented has been reduced by removing the interaction effect since a previous analysis proved its negligibility. In fact, the value of F related to the interaction effect is 0.41 corresponding to a value of 0.558.
The analysis shows a significant curvature value, whose value is less than 0.001, a main effect value less than 0.001 ensuring a huge significant influence; specifically, the value related to rpm is less than 0.001 and the value related to TS−1 is 0.001. Both the effects can be considered very significant. Regression model coefficients are estimated as shown in Table 7. It was also performed a residuals analysis that seems to verify the hypothesis made on the residuals, so they can be assumed to be normally distributed with null average and constant standard deviation. By the way, the regression model cannot fit well the data, even if the value related to the lack of fit is very high, since the adopted design cannot allow second-order estimate effects being significant in this case.
As for the nugget zone, a 22 factorial design with two centre points and three repeated measurements was carried out for the elongation. Experimental data of yield strength are shown in Table 8.
In this case, a simple ANOVA was carried out since the analysis of variability did not show any significant effect. Reliance on this consideration would deem the weights negligible. The ANOVA related to the experimental data of elongation is shown in Table 9.
The analysis shows a curvature value of 0.070 that proves its presence; the main effect value is 0.051; specifically, the value related to rpm is 0.037 which is significant. The value related to TS−1 is 0.198, which is not significant. The residuals can be assumed to be normally distributed, with null average and constant standard deviation. An important point is that the high value related to the curvature proves even a lack of fit of the regression model. Looking data of ultimate strength, yield strength, and elongation relative of the analysis of variance, it can be seen that the influence effects are not significant.
Regarding the optimization of the response surface related to the yield strength, which is the main mechanical property of interest, an algorithm of steepest ascent to the yield strength response of the nugget zone was proposed. Assuming that , where i = 1, 2, …, k identifies the factors and [−1; 0; 1] are the coded variables identifying the lower, the centre, and the upper level of the factors, respectively, the algorithm procedure can be described according to the following procedure: defining the step length for one of the process variables and the step length of the other variables can be defined by using equation (1). Finally, turning the step length ∆xi from the coded variable into natural variables.
This algorithm can be applied only when the regression model is of the first order and does not include interactions. Runs obtained by implementing the aforementioned algorithm are shown in Table 10:
The results are similar to those expected, and this ascertainment guarantees the reliability of the previous analysis. The average of the yield strength has increased from a value of 363.33 MPa (point: 1377.5 rpm; 0.267 s/mm) to values of 399.25 MPa (point: 1250 rpm; 0.241 s/mm) and 395.8 MPa (point: 1150 rpm; 0.213 s/mm). The region where the two last points are set is the optimum region; a new investigation using a central composite design would have been necessary to determine the exact position of the yield strength optimum point. Following the next steps of the steepest path, the yield strength decreases, which suggests that the steepest path is leaving from the optimum region. A similar behaviour is acted also by the ultimate strength, but steepest yield strength of the nugget zone path is even one of the improvement directions of yield and ultimate strength of the heat-affected zone. Seemingly, the mechanical properties of the HAZ improve markedly when the process gets colder and colder.
As for the shear testing, the results showed averaged failure loads, since three specimens were tested for each run. In Table 11, a summary of averaged failure loads is presented.
According to the tensile testing results, the runs corresponding to the first two steepest ascent algorithms show the best performances for yield strength, especially the run (1150 rpm; 0.213 s/mm) which had an averaged failure load of 10125.16 N. Those failures occur always in the nugget zone explaining the correlation between tensile and shear testing.
Vickers microhardness measurements, which were carried out from the centre of the welding along the cross section passing through HAZ, TMAZ, and nugget zone, are shown in Figure 4.
The results reveal a switch of the best mechanical behaviour zones, from welds showing better hardness performances in the nugget zone to welds showing better hardness performances in the HAZ. Looking at the welds carried out by setting rpm equal to 1050 and TS−1 equal to 0.184 s/mm and the ones obtained with rpm equal to 1150 and TS−1 equal to 0.213 s/mm (Figure 4), the switch in hardness performances confirms the results already obtained with tensile testing. In fact, as well as tensile testing has held up when the runs get very cold (namely, the following runs: rpm = 950; TS−1 = 0.155 s/mm, rpm = 1050; TS−1 = 0.184 s/mm), improved mechanical behaviour in the HAZ exceeds nugget zone performance. Thus, microhardness testing confirms the same trend of tensile testing.
Optical investigations included the study of the oxide film; this kind of defect, which is a feature for friction stir welding technology, is commonly called hook defect because of its shape. Hook defects must be taken seriously since they represent cross-section thinning, and for this reason, they may explain shear testing results . Hook’s distances are shown in Table 12.
This study does not contain statistical evidence for hook’s distance and shear tests correlation; in fact, a linear regression analysis showed a very low coefficient of determination (R2 = 0.102 and R2(adj) = 0) and a very high ANOVA value (F = 0.68; value = 0.44).
As for the grain size number of the nugget zone, the results are shown in Table 13.
It can be seen a decrease in grain size for runs featured by high heat generation. A regression analysis between grain size estimates and heat index was carried out as proof. Heat index is considered as a measurement of heat generated during welding, as the higher the index, the higher is the heat generated. It was defined by using the following equation:
The coefficient of determination presented a high value (R2 = 0.68, R2(adj) = 0.63), which suggests a good correlation between the two parameters, as well the analysis of variance as proof of a significant influence of the heat index on grain size showing a value of 0.012. In other words, these analyses seem to ensure this assertion: “As big is the heat generation as big is the grain size in the Nugget Zone.” The mechanical behaviour of the welds is explained by regression analysis between hardness Vickers tests and grain size. Such an analysis revealed a lack of correlation between the two properties, thus disagreeing with the Hall–Petch rule. So, the causes of the mechanical behaviour should be found in the differences of the weld microstructure; i.e., different temporal thermal profiles produce different precipitate compositions, where some are more hardened than others [22, 32, 33].
In this study, a procedure to optimize the mechanical behaviour of friction stir welded joints was developed. The yield strength was chosen as a factor of interest, and its response was optimized by using a response surface method. This method, which consisted of a central composite design and a subsequent steepest ascent algorithm, provided for optimal yield strength surface. Tensile tests on base metal specimens were conducted with benchmarking measures obtained from processed samples proving the high performance of FSW technology. The best performances, which were evaluated choosing the minimum value between HAZ and nugget zone ratio for both yield and ultimate strength, were found under these welding conditions: rpm = 1250 and TS−1 = 0.241 s/mm. Shear testing confirmed the results provided by tensile tests. Hook’s defect was studied, and the grain size of the nugget zone was estimated.
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
V. M. Magalhães, C. Leitão, and D. M. Rodrigues, “Friction stir welding industrialisation and research status,” Science and Technology of Welding and Joining, vol. 23, no. 5, pp. 400–409, 2017.View at: Publisher Site | Google Scholar
X. Liu, S. Lan, and J. Ni, “Analysis of process parameters effects on friction stir welding of dissimilar aluminum alloy to advanced high strength steel,” Materials & Design, vol. 59, pp. 50–62, 2014.View at: Publisher Site | Google Scholar
R. Nandan, T. Debroy, and H. Bhadeshia, “Recent advances in friction-stir welding-process, weldment structure and properties,” Progress in Materials Science, vol. 53, no. 6, pp. 980–1023, 2008.View at: Publisher Site | Google Scholar
Y. S. Sato, M. Urata, and H. Kokawa, “Parameters controlling microstructure and hardness during friction-stir welding of precipitation-hardenable aluminum alloy 6063,” Metallurgical and Materials Transactions A, vol. 33, no. 3, pp. 625–635, 2002.View at: Publisher Site | Google Scholar
A. Ali, M. W. Brown, C. A. Rodopoulos, and S. Gardiner, “Characterization of 2024-T351 friction stir welding joints,” Journal of Failure Analysis and Prevention, vol. 6, no. 4, pp. 83–96, 2006.View at: Publisher Site | Google Scholar
M. W. Mahoney, C. G. Rhodes, J. G. Flintoff, W. H. Bingel, and R. A. Spurling, “Properties of friction-stir-welded 7075 T651 aluminum,” Metallurgical and Materials Transactions A, vol. 29, no. 7, pp. 1955–1964, 1998.View at: Publisher Site | Google Scholar
R. Nandan, G. G. Roy, T. J. Lienert, and T. Debroy, “Three-dimensional heat and material flow during friction stir welding of mild steel,” Acta Materialia, vol. 55, no. 3, pp. 883–895, 2007.View at: Publisher Site | Google Scholar
H. J. Liu, H. Fujii, M. Maeda, and K. Nogi, “Tensile properties and fracture locations of friction-stir-welded joints of 2017-T351 aluminum alloy,” Journal of Materials Processing Technology, vol. 142, no. 3, pp. 692–696, 2003.View at: Publisher Site | Google Scholar
R. F. Gunst, “Response surface methodology: process and product optimization using designed experiments,” Technometrics, vol. 38, no. 3, pp. 284–286, 1996.View at: Publisher Site | Google Scholar
J. Antony, “Simultaneous optimisation of multiple quality characteristics in manufacturing processes using taguchi's quality loss function,” International Journal of Advanced Manufacturing Technology, vol. 17, no. 2, pp. 134–138, 2001.View at: Publisher Site | Google Scholar
T. V. Sibalija and V. D. Majstorovic, “An integrated approach to optimise parameter design of multi-response processes based on Taguchi method and artificial intelligence,” Journal of Intelligent Manufacturing, vol. 23, no. 5, pp. 1511–1528, 2010.View at: Publisher Site | Google Scholar
M. J. Anderson, “A new method for non-parametric multivariate analysis of variance,” Austral Ecology, vol. 26, no. 1, pp. 32–46, 2001.View at: Publisher Site | Google Scholar
W. B. Schaufeli, M. Salanova, V. González-romá, and A. B. Bakker, “The measurement of engagement and burnout: a two sample confirmatory factor analytic approach,” Journal of Happiness Studies, vol. 3, no. 1, pp. 71–92, 2002.View at: Publisher Site | Google Scholar
N. Alagumurthi, K. Palaniradja, and V. Soundararajan, “Optimization of grinding process through design of experiment (DOE)-a comparative study,” Materials and Manufacturing Processes, vol. 21, no. 1, pp. 19–21, 2006.View at: Publisher Site | Google Scholar
W. Chen, J. K. Allen, K.-L. Tsui, and F. Mistree, “A procedure for robust design: minimizing variations caused by noise factors and control factors,” Journal of Mechanical Design, vol. 118, no. 4, p. 478, 1996.View at: Publisher Site | Google Scholar
G. Vicente, A. Coteron, M. Martinez, and J. Aracil, “Application of the factorial design of experiments and response surface methodology to optimize biodiesel production,” Industrial Crops and Products, vol. 8, no. 1, pp. 29–35, 1998.View at: Publisher Site | Google Scholar
B. K. Rout and R. K. Mittal, “Parametric design optimization of 2-DOF R-R planar manipulator-A design of experiment approach,” Robotics and Computer-Integrated Manufacturing, vol. 24, no. 2, pp. 239–248, 2008.View at: Publisher Site | Google Scholar
A. E. Bryson and W. F. Denham, “A steepest-ascent method for solving optimum programming problems,” Journal of Applied Mechanics, vol. 29, no. 2, p. 247, 1962.View at: Publisher Site | Google Scholar
M. Jones, P. Heurtier, C. Desrayaud, F. Montheillet, D. Allehaux, and J. Driver, “Correlation between microstructure and microhardness in a friction stir welded 2024 aluminium alloy,” Scripta Materialia, vol. 52, no. 8, pp. 693–697, 2005.View at: Publisher Site | Google Scholar
G. Liu, L. E. Murr, C.-S. Niou, J. C. McClure, and F. R. Vega, “Microstructural aspects of the friction-stir welding of 6061-T6 aluminum,” Scripta Materialia, vol. 37, no. 3, pp. 355–361, 1997.View at: Publisher Site | Google Scholar
G. D’Urso and C. Giardini, “The influence of process parameters and tool geometry on mechanical properties of friction stir welded aluminum lap joints,” International Journal of Material Forming, vol. 3, no. 1, pp. 1011–1014, 2010.View at: Publisher Site | Google Scholar
Y. S. Sato, M. Urata, H. Kokawa, and K. Ikeda, “Hall-Petch relationship in friction stir welds of equal channel angular-pressed aluminium alloys,” Materials Science and Engineering: A, vol. 354, no. 1-2, pp. 298–305, 2003.View at: Publisher Site | Google Scholar
P. M. G. P. Moreira, T. Santos, S. M. O. Tavares, V. Richter-Trummer, P. Vilaça, and P. M. S. T. de Castro, “Mechanical and metallurgical characterization of friction stir welding joints of AA6061-T6 with AA6082-T6,” Materials & Design, vol. 30, no. 1, pp. 180–187, 2009.View at: Publisher Site | Google Scholar
L. Commin, M. Dumont, J.-E. Masse, and L. Barrallier, “Friction stir welding of AZ31 magnesium alloy rolled sheets: influence of processing parameters,” Acta Materialia, vol. 57, no. 2, pp. 326–334, 2009.View at: Publisher Site | Google Scholar
W. Zhang, Y. Shen, Y. Yan, and R. Guo, “Dissimilar friction stir welding of 6061 Al to T2 pure Cu adopting tooth-shaped joint configuration: microstructure and mechanical properties,” Materials Science and Engineering: A, vol. 690, pp. 355–364, 2017.View at: Publisher Site | Google Scholar
R. S. Mishra and Z. Y. Ma, “Friction stir welding and processing,” Materials Science and Engineering: R: Reports, vol. 50, no. 1-2, pp. 1–78, 2005.View at: Publisher Site | Google Scholar
A. Murphy, M. Price, R. Curran, and P. Wang, “Integration of strength and process modeling of friction-stir-welded fuselage panels,” Journal of Aerospace Computing, Information, and Communication, vol. 3, no. 4, pp. 159–176, 2006.View at: Publisher Site | Google Scholar
A. V. Sergueeva, J. Zhou, B. E. Meacham, and D. J. Branagan, “Gage length and sample size effect on measured properties during tensile testing,” Materials Science and Engineering: A, vol. 526, no. 1-2, pp. 79–83, 2009.View at: Publisher Site | Google Scholar
M. Ahmadi, F. Vahabzadeh, B. Bonakdarpour, E. Mofarrah, and M. Mehranian, “Application of the central composite design and response surface methodology to the advanced treatment of olive oil processing wastewater using Fenton's peroxidation,” Journal of Hazardous Materials, vol. 123, no. 1–3, pp. 187–195, 2005.View at: Publisher Site | Google Scholar
W. S. Cleveland and S. J. Devlin, “Locally weighted regression: an approach to regression analysis by local fitting,” Journal of the American Statistical Association, vol. 83, no. 403, pp. 596–610, 1988.View at: Publisher Site | Google Scholar
W. Li, J. Li, Z. Zhang, D. Gao, W. Wang, and C. Dong, “Improving mechanical properties of pinless friction stir spot welded joints by eliminating hook defect,” Materials & Design (1980–2015), vol. 62, pp. 247–254, 2014.View at: Publisher Site | Google Scholar
N. Dialami, M. Chiumenti, M. Cervera, and C. Agelet de Saracibar, “Challenges in thermo-mechanical analysis of friction stir welding processes,” Archives of Computational Methods in Engineering, vol. 24, no. 1, pp. 189–225, 2016.View at: Publisher Site | Google Scholar
M. M. Attallah, C. L. Davis, and M. Strangwood, “Microstructure-microhardness relationships in friction stir welded AA5251,” Journal of Materials Science, vol. 42, no. 17, pp. 7299–7306, 2007.View at: Publisher Site | Google Scholar