Parametric Optimization of FDM Process for Improving Mechanical Strengths Using Taguchi Method and Response Surface Method: A Comparative Investigation

: In the present study, a comparison of two widely used optimization approaches for fused deposition modeling (FDM), that is, Taguchi method in contrast with response surface method (RSM), was investigated. Four operating parameters, namely extrusion temperature, layer thickness, raster width, print speed, and their interaction terms, were identiﬁed as control variables with three levels, while tensile strength and compressive strength were selected responses. L27 orthogonal array and face-centered central composite design (FCCCD) were used for the experimental approach for Taguchi and RSM, respectively. The signal-to-noise (S/N) ratio and analysis of variance (ANOVA) were employed to ﬁnd the optimal FDM parameter combination as well as the main factor that affect the performance of the PLA samples. Based on experimental results, it was observed that conclusions about signiﬁcant ranking of parameters on FDM process from these two methods were different. However, both the Taguchi method and RSM succeed in predicting better results compared with the original groups. In addition, the optimum combinations for tensile strength and compressive strength obtained from the RSM were 2.11% and 8.15% higher than Taguchi method, respectively.


Introduction
Additive manufacturing (AM), also known as 3D printing, can produce objects through depositing materials (e.g., powdered metal or plastic) layer by layer based on the digital model [1].In the recent past, AM has become a rapidly emerging innovative approach in many industrial fields.Fused Deposition Modeling (FDM) is the most extensively used AM approach, which extrudes the melted filaments (e.g., plastic materials such as ABS or PLA) from the nozzle to form the component on the platform according to the prescribed manner.Due to low cost, fast fabrication, minimum energy consumption, and low material wastage, FDM technique accounts for a large proportion in the AM market.With the development of new materials and improvement of process parameters, FDM shows potential application in high-end fields such as aviation, medicine, electronics, etc. [2].However, a well-identified limitation of the FDM products is weak mechanical properties.The input parameters play a crucial role in determining characteristics of the part quantitatively or qualitatively.Therefore, the suitable setting of parameters associated with different performances remains a challenge and must be determined.
A considerable amount of research work has been carried out to optimize and improve parameters of FDM process, and various techniques such as the Taguchi method [3,4], response surface method (RSM) [5,6], full factorial design [7][8][9], fractional factorial design [10,11], and other advanced techniques [12][13][14][15] have been used for analyzing and optimizing.Typically, the most widely accepted approaches are the Taguchi method and RSM, providing efficient techniques for process design and product optimization.
The Taguchi method provides an effective and reliable approach in simplification of feasibility of study and experimental plan of different process parameters.This is especially important to AM where the cost of products produced is still high.Liu et al. [16] analyzed five parameters (build orientation, layer thickness, raster width, air gap, and raster angle) on impact, flexural, and tensile strengths of ABS built parts using gray Taguchi analysis.The author observed that build orientation was the most significant parameter, followed by layer thickness and raster angle.Rinanto et al. [17] investigated the tensile strength of PLA parts for a set of extrusion temperature, raster angle, and infill density via Taguchi approach.The conclusion of their results was that 210 • C extrusion temperature, 45 • raster angle, and 40% infill density maximized the performance.Deng et al. [18] designed the experiment using L9 orthogonal array for tensile properties of PEEK, and the preferable parameters were low layer thickness, high extrusion temperature, and print speed.Sood et al. [19] considered the effect of five factors namely raster width, layer thickness, air gap, raster angle, and build orientation, and their interactions on ABS dimensional accuracy via grey Taguchi method.Experimental outcomes indicated that optimal parameter settings for different mechanical characteristics were not the same.Chacón et al. [20] used a low cost printer to produce PLA samples and analyzed the important of different process parameters (layer thickness, build orientation, and print speed) and their interactions on flexural and tensile properties.The study acknowledged that the effects of build orientation and layer thickness were coupled in impacting the tensile property.
The RSM is considered another promising optimization method, which gives very low standard errors to experimental verification.The RSMs used most frequently are the central composite design (CCD) and face-centered central composite design (FCCCD).Srivastava et al. [21] determined the role of six process parameters on build time and support material via CCD.Among the analyzed parameters, raster width, contour width, layer thickness, and air gap were significant parameters, while raster angle and build orientation were deemed insignificant in comparison.Other than tensile property, Percoco et al. [22] also experimentally studied the impact of raster angle, immersion time and raster width on compressive properties of chemical dipped specimens using 3-level CCD.Their results reported that raster angle had a quite weak effect on the compressive behavior, which however increased with the increase of raster width.The functional relationships between process parameters (build orientation, layer thickness, raster angle, air gap and raster width) and strengths (tensile, flexural and impact) were determined by Sood et al. [5] and Panda et al. [23], respectively.Similar outcomes from FCCCD indicated that zero air gap and thick raster width improved the mechanical properties, while raster angle with small value was not preferable.
In fact, there are many other studies based on Taguchi and RSM that are not listed due to space limitations.Most conclusions about effects of process parameters in the review literature of FDM are also drawn according to the results of these two approaches [24][25][26].From all the referred studies, the main concerns through different design of experiment (DOE) approaches usually to be answered are: (1) Among all parameters, which one has the maximum/minimum impact on FDM performance?(2) Which level of a parameter is the most beneficial to improve product's behavior?Statistical tools such as S/N analysis, ANOVA, interaction charts, and main effect plots are commonly used.However, an easily overlooked issue in FDM research is whether the conclusions obtained by different methods are the same or consistent, which have great effects on FDM parameters chosen for users and manufacturers.However, few studies focus on the comparison of different DOE approaches.Rashed et al. [27] compared Taguchi design and full factorial design to determine the influences of parameters on surface roughness and mechanical strength of Nylon 6/66 copolymer.Some errors were found in determining the influences of parameters on impact strength when using Taguchi DOE.Kechagias and Vidakis [28] assessed the efficiency of full factorial design in contrast with Box-Behnken design (BBD) of RSM in the tensile strength of PA 12 specimens, proving that, for the quadratic regression modeling, these two approaches have similar efficiency.Only in the research of Tontowi et al. [29] were Taguchi method and RSM compared on dimensional accuracy and tensile strength by investigating raster angle, build orientation, and layer thickness.The authors recommended that RSM was a better prediction approach.However, for these existing limited articles, the above two concerned questions have not been all answered directly.To our knowledge, a comprehensive comparison between the Taguchi method and RSM has not been wellestablished in the literature, which remains an urgent and vital problem to be solved.
This study attempted to present the comparison between Taguchi method and RSM by analyzing FDM parameters for optimum performance of PLA samples in terms of tension and compression.Four process parameters varied on three levels, that is, extrusion temperature, layer thickness, print speed, and raster width, along with their interaction terms, were analyzed in this research.In the following, the description of experimental design through these two approaches to investigate the FDM process parameters was given.The optimal combinations with regard to mechanical characteristics, which were performed using tension and compression tests, were considered.Results were then validated by Taguchi method to determine the predictive model when compared to RSM.Statical data were discussed, and finally the study ended with the findings of the investigation.

Materials and Equipment
The FDM materials range includes engineering polymers (e.g., nylon and PETG), thermoplastics (e.g., PLA and ABS) and composite materials (e.g., carbon fiber composite and glass fiber composite), among which PLA is the most commonly used filament because of its biodegradability, non-toxicity, low glass transition temperature, low price, and moderate flexibility [30,31].In this research, the material used for sample fabrication is natural PLA filament (provided by Esun Industrial Co., Ltd., Shenzhen, China) with a diameter of 1.75 mm.The FUNMAT HT 3D printer with IntamSuite slice software (provided by Intamsys Technologies Co., Ltd., Shanghai, China) is used in the study to fabricate experimental prints.Taguchi and RSM analysis are conducted using Minitab software (Version: 21.1).MIT-20 microcomputer controlled electronic universal testing machine (provided by Sanfeng Instrument Co., Ltd., Changzhou, China) is used to measure mechanical strengths, as shown in Figure 1.
rameters on impact strength when using Taguchi DOE.Kechagias and Vidakis [28] assessed the efficiency of full factorial design in contrast with Box-Behnken design (BBD) of RSM in the tensile strength of PA 12 specimens, proving that, for the quadratic regression modeling, these two approaches have similar efficiency.Only in the research of Tontowi et al. [29] were Taguchi method and RSM compared on dimensional accuracy and tensile strength by investigating raster angle, build orientation, and layer thickness.The authors recommended that RSM was a better prediction approach.However, for these existing limited articles, the above two concerned questions have not been all answered directly.To our knowledge, a comprehensive comparison between the Taguchi method and RSM has not been well-established in the literature, which remains an urgent and vital problem to be solved.
This study attempted to present the comparison between Taguchi method and RSM by analyzing FDM parameters for optimum performance of PLA samples in terms of tension and compression.Four process parameters varied on three levels, that is, extrusion temperature, layer thickness, print speed, and raster width, along with their interaction terms, were analyzed in this research.In the following, the description of experimental design through these two approaches to investigate the FDM process parameters was given.The optimal combinations with regard to mechanical characteristics, which were performed using tension and compression tests, were considered.Results were then validated by Taguchi method to determine the predictive model when compared to RSM.Statical data were discussed, and finally the study ended with the findings of the investigation.

Materials and Equipment
The FDM materials range includes engineering polymers (e.g., nylon and PETG), thermoplastics (e.g., PLA and ABS) and composite materials (e.g., carbon fiber composite and glass fiber composite), among which PLA is the most commonly used filament because of its biodegradability, non-toxicity, low glass transition temperature, low price, and moderate flexibility [30,31].In this research, the material used for sample fabrication is natural PLA filament (provided by Esun Industrial Co., Ltd., Shenzhen, China) with a diameter of 1.75 mm.The FUNMAT HT 3D printer with IntamSuite slice software (provided by Intamsys Technologies Co., Ltd., Shanghai, China) is used in the study to fabricate experimental prints.Taguchi and RSM analysis are conducted using Minitab software (Version: 21.1).MIT-20 microcomputer controlled electronic universal testing machine (provided by Sanfeng Instrument Co., Ltd., Changzhou, China) is used to measure mechanical strengths, as shown in Figure 1.

Parameters Selection
There are various parameters in FDM process found to impact mechanical characteristics of printed parts, such as infill density, air gap, number of contours, platform temperature, etc.Meanwhile, FDM process deposits material in a criss-cross way which results in direction dependence, leading to products with anisotropy.It has been proved that the anisotropic behavior of FDM parts has close relationship with parameters such as build orientation [32], raster angle [33], and infill pattern [34].Therefore, to exclude the interference of anisotropic property, four process parameters in Figure 2 ((A) extrusion temperature, (B) layer thickness, (C) print speed and (D) raster width) were selected as the controllable variables, which are considered to have little effect on FDM anisotropy [33].The values for these parameters are determined based on practical experience in processing, as presented in Table 1.According to the literature, some prominent characteristics of the quality for FDM parts are mechanical strength, dimensional accuracy, surface roughness, and manufacturing time.Obviously, the mechanical strength is the most concerning and important property to the consumers, which is also selected as the response in the study.

Parameters Selection
There are various parameters in FDM process found to impact mechanical characteristics of printed parts, such as infill density, air gap, number of contours, platform temperature, etc.Meanwhile, FDM process deposits material in a criss-cross way which results in direction dependence, leading to products with anisotropy.It has been proved that the anisotropic behavior of FDM parts has close relationship with parameters such as build orientation [32], raster angle [33], and infill pattern [34].Therefore, to exclude the interference of anisotropic property, four process parameters in Figure 2 ((A) extrusion temperature, (B) layer thickness, (C) print speed and (D) raster width) were selected as the controllable variables, which are considered to have little effect on FDM anisotropy [33].The values for these parameters are determined based on practical experience in processing, as presented in Table 1.According to the literature, some prominent characteristics of the quality for FDM parts are mechanical strength, dimensional accuracy, surface roughness, and manufacturing time.Obviously, the mechanical strength is the most concerning and important property to the consumers, which is also selected as the response in the study.

Sample Preparation
Samples for strength tests were fabricated following ASTM standard, namely ASTM D638 for tensile strength and ASTM D695 for compressive strength, respectively, the dimensions of which are depicted in Figure 3.To exclude the interference of mechanical anisotropy, all samples were built at 0° raster angle with flat build orientation.A 0.6 mm nozzle diameter, 100% infill density, and 70 °C platform temperature was maintained for all prints.For a given group of parameter settings, each sample consisted of three identical specimens, which were conditioned in a laboratory environment at room temperature (24 °C) for 24 h before testing to achieve equilibrium humidity content.The mean of each experiment trial was taken as the represented value of respective strength.

Sample Preparation
Samples for strength tests were fabricated following ASTM standard, namely ASTM D638 for tensile strength and ASTM D695 for compressive strength, respectively, the dimensions of which are depicted in Figure 3.To exclude the interference of mechanical anisotropy, all samples were built at 0 • raster angle with flat build orientation.A 0.6 mm nozzle diameter, 100% infill density, and 70 • C platform temperature was maintained for all prints.For a given group of parameter settings, each sample consisted of three identical specimens, which were conditioned in a laboratory environment at room temperature (24 • C) for 24 h before testing to achieve equilibrium humidity content.The mean of each experiment trial was taken as the represented value of respective strength.

Taguchi Method
The Taguchi method uses a fractional factorial design called orthogonal array (OA) to determine the influence of parameters and optimum combinations with fewer experiment runs.For the Taguchi method used in many FDM studies, interaction effects between control factors are often regarded as less important and can be ignored.However, some FDM parameters are dependent on the same factor and should be considered coupled.For example, print time and raster width are affected by layer thickness.Infill density significantly impact the print speed, which can be changed by adjusting raster width and air gap [24].Therefore, it is necessary to analyze the role of the intersection effect of related parameters and in what amount it is.In this study, the effects of three interaction terms, that is, layer thickness and print speed (B × C), layer thickness and raster width (B × D), and print speed and raster width (C × D) are taken into consideration.To select a proper orthogonal array for the experimental design, the degrees of freedom should be calculated, which can be determined as [35]: DOF = (level number − 1) × factor number + (level number − 1) 2 × interaction number For this experiment, DOF = 2 × 4 + 4 × 3 = 20.Therefore, the orthogonal array L27(3 13 ) was used, consisting of 13 columns for assigning four single factors and three two-factor interactions.Table 2 describes the factors and levels of the Taguchi design.The experimental results of 81 specimens for mechanical strengths are tabulated in Table 3.

Taguchi Method
The Taguchi method uses a fractional factorial design called orthogonal array (OA) to determine the influence of parameters and optimum combinations with fewer experiment runs.For the Taguchi method used in many FDM studies, interaction effects between control factors are often regarded as less important and can be ignored.However, some FDM parameters are dependent on the same factor and should be considered coupled.For example, print time and raster width are affected by layer thickness.Infill density significantly impact the print speed, which can be changed by adjusting raster width and air gap [24].Therefore, it is necessary to analyze the role of the intersection effect of related parameters and in what amount it is.In this study, the effects of three interaction terms, that is, layer thickness and print speed (B × C), layer thickness and raster width (B × D), and print speed and raster width (C × D) are taken into consideration.To select a proper orthogonal array for the experimental design, the degrees of freedom should be calculated, which can be determined as [35]: DOF = (level number − 1) × factor number + (level number − 1) 2 × interaction number For this experiment, DOF = 2 × 4 + 4 × 3 = 20.Therefore, the orthogonal array L 27 (3 13 ) was used, consisting of 13 columns for assigning four single factors and three two-factor interactions.Table 2 describes the factors and levels of the Taguchi design.The experimental results of 81 specimens for mechanical strengths are tabulated in Table 3.

S/N Ratio Analysis
The signal to noise ratio (S/N) is used to estimate the performance characteristics deviating from the desired values.The output of the orthogonal arrays is optimized with respect to the S/N ratio of the response.For response of mechanical strength, mean square deviation (MSD) for 'the-larger-the-better' quality characteristic is applied.
S/N = −10 log 10 (MSD) where y i is the mean mechanical strength for i th set experimental sample, n is the number of experiments in the orthogonal array.The main effect plots for mechanical strengths are shown in Figure 4. On the one hand, the S/N ratio for tensile strength (Figure 4a) decreases with increasing extrusion temperature, while it increases with increasing layer thickness and raster width.It can also be noted that, with an increase in print speed, the S/N ratio response first decreases and then increases slightly.In addition, the significances of parameters without considering interaction effects are B > A > D > C. The optimum combination is A -1 B 1 C -1 D 1 .When interaction terms are not ignored, the important orders are The best level of partial parameters for optimum performance is A -1 B 1 D 1 .However, considering the interaction term (C × D) has a higher priority than C, the level of C should be given by comparing different combinations of C and D. According to the interaction plots in Figure 5, when D is set at a high level of 1, tensile strength decreases as C is changed from −1 to 1, meaning C -1 D 1 is the best combination to maximize the tensile property.That is, the optimum levels considering the interaction effect remain A -1 B 1 C -1 D 1 (i.e., extrusion temperature 205 • C, layer thickness 0.35 mm, print speed 40 mm/s and raster width 0.6 mm).On the other hand, for compressive strength (Figure 4b), the S/N ratio first increases and then decreases with increase in extrusion temperature and layer thickness.The opposite behavior is observed from print speed and raster width, respectively.The important factors sorted in descending order without interaction terms are: B > A > D > C. The optimum setting levels are A 0 B 0 C -1 D 1 .When interaction effects are taken into consideration, the order becomes B > A > D > (C × D) > (B × D) > C > (B × C).In this situation, the optimum levels of partial parameters for compression are A 0 B 0 D 1 .Since the interaction term (C × D) has a higher priority than C, therefore the optimum combination by considering the interaction effect from Figure 5 becomes A 0 B 0 C -1 D 1 , which is also consistent with the previously unconsidered one (i.e., extrusion temperature 215 • C, layer thickness 0.25 mm, print speed 40 mm/s and raster width 0.6 mm).
Furthermore, the interactions plot in Figure 5 shows that no trends run parallel in a coupled group, which is a sign of strong interactions among the different experimental factors investigated.An interesting phenomenon is that optimum combinations are the same when considering and not considering interaction terms, both for tensile strength and compressive strength.One possible explanation is the importance of interaction term is less than individual process parameter, therefore the results remain.However, this phenomenon seems more likely to occur in Taguchi design.As we know, the orthogonal array of Taguchi's design matrix has a number limitation of factors [36] (e.g., there are at most three two-factor interactions can be investigated in L 27 (3 13 )), thus it is impossible to account for all interaction effects.Since the chosen interaction terms are based on past experience and practical processing, it is likely to concentrate on unimportant interactions but neglects the significant ones.This phenomenon is more evident when compared with RSM, which will be discussed later.On the other hand, for compressive strength (Figure 4b), the S/N ratio first increases and then decreases with increase in extrusion temperature and layer thickness.The opposite behavior is observed from print speed and raster width, respectively.The important factors sorted in descending order without interaction terms are: B > A > D > C. The optimum setting levels are A0B0C-1D1.When interaction effects are taken into consideration, the order becomes B > A > D > (C × D) > (B × D) > C > (B × C).In this situation, the optimum levels of partial parameters for compression are A0B0D1.Since the interaction term (C × D) has a higher priority than C, therefore the optimum combination by considering the interaction effect from Figure 5 becomes A0B0C-1D1, which is also consistent with the previously unconsidered one (i.e., extrusion temperature 215 °C, layer thickness 0.25 mm, print speed 40 mm/s and raster width 0.6 mm).
Furthermore, the interactions plot in Figure 5 shows that no trends run parallel in a coupled group, which is a sign of strong interactions among the different experimental factors investigated.An interesting phenomenon is that optimum combinations are the same when considering and not considering interaction terms, both for tensile strength and compressive strength.One possible explanation is the importance of interaction term is less than individual process parameter, therefore the results remain.However, this phenomenon seems more likely to occur in Taguchi design.As we know, the orthogonal array of Taguchi's design matrix has a number limitation of factors [36] (e.g., there are at most three two-factor interactions can be investigated in L27(3 13 )), thus it is impossible to account for all interaction effects.Since the chosen interaction terms are based on past experience

ANOVA Analysis
The purpose of the analysis of variance (ANOVA) is to investigate which parameters significantly affected the quality characteristic and the relative percent influence, which has been applied at 95% confidence interval level.For importance and significance check, F value and P value given in ANOVA table are used.If P value is less than 0.05, significance of related term is established.
Observations are made from Table 4 that the significant ranking of process parameters in terms of F value by ANOVA for tension and compression is , respectively, which are just in agreement with results of S/N analysis before.That is, S/N analysis and ANOVA mutually verify the correctness of results of each other.In addition, for tensile strength, layer thickness (B) significantly affects the tensile strength of the specimens as P value is less than 0.05.However, for compressive strength, none of the factors are significant as their P values are all above 0.05.

Comparative Confirmation Test
It is found that optimum parameters and their levels just correspond to existing experiment of the orthogonal array in Table 3.That is, the optimum tension A -1 B 1 C -1 D 1 is No.7, while the optimum compression is No.13, respectively.Furthermore, the strength values of these two settings are just the maximum of each array, proving the optimum combinations obtained from the S/N ratio analysis are correct.In terms of percentage, the difference between predicted values and actual values in Table 5 are 4.22% in tension and 1.89% in compression, respectively.Thus, these results verified the precision and efficiency of the Taguchi method.

Experimental Design
The empirical model for the response surface method (RSM) is usually built by central composite design (CCD), which contains five levels for each factor: zero level (center point), ±1 level (factorial points), and ±α level (axial points, or star points).Owing to machine constraints, face-centered central composite design (FCCCD) where α = 1 is usually considered, ensuring the position of axial point on the centers of faces of the cube [23], which is also adopted in this study.In this situation, the three levels of factors for RSM are the same as the Taguchi method.These two approaches can be regarded as part of full fractional design (3 4 = 81 experimental runs required), which is easier for comparison on the same scale.Therefore, total 30 runs with 6 center points, 8 axial points and 16 factorial points, are used in this study, as shown in Table 6.The number of experiments required is a little more than Taguchi's.
The quadratic regression model is developed, predicting the dependent responses in terms of independent variables, along with their interactions: where y is the response; x i is i th factor with coded levels (A, B, C and D); β 0 is constant, β i , β ii and β ij are the coefficient values for linear, quadratic, and interaction terms, respectively; k is the number of design factors, which is 4 in this research.

Results and Discussion
The Anderson-Darling (AD) normality tests of the residual of the respective strength are applied.Figure 6 shows that both mechanical strengths through RSM approach fit the normality distribution, whose P values are higher than 0.05.In addition, the ANOVA from Table 7 indicates that, for tensile strength, factors with significant effects are ordered as: In the current case, P value of 0.0003 denotes extrusion temperature (A) as significant.For compressive strength, the important orders are changed as: D The parameters with significant effects are D, (B × D), (A × B), and D 2 , with P values lower than 0.05.
The conclusions about important orders of parameters are different from the Taguchi method.This can be explained as follows: On the one hand, most parameter settings of tested samples for Taguchi (27 groups) and RSM (30 groups) are different, and only four groups are the same according to the statistics, as shown in Figure 7. Therefore, the regression model for each approach, which is greatly affected by the experimental data, is different.On the other hand, the empirical model for RSM is usually full quadratic with square terms and all interaction terms, which cannot be achieved by Taguchi method.Since there are more terms to be considered in RSM, the importance of the factor is changed accordingly (For example, factor A appears in linear, square, and interaction forms in RSM).As a result, the significant rankings of parameters obtained through these two approaches are usually different, which is applicable to the optimum combination.It should be pointed out that the effects of process parameters have not been discussed in this study, unlike many other studies done.This is because the conclusion of the behavior of a specific process parameter is affected by various factors, such as the DOE approach (just as this study presents), materials, environmental conditions, experimental standard, etc.Furthermore, the coupled effect of different parameter combinations will also lead to significant discrepancies.For example, Oubalouch et al. [37] investigated layer thickness, raster angle, and infill density on tensile property of PLA parts using full factorial experiment.In the case of 10% infill density, the tensile strength was maximum at medium layer thickness.In contrast, the lower layer thickness, the higher strength was found for both 50% and 100% infill densities.Based on the same approach and material, In the study of Rajpurohit and Dave [38], the authors also determined the effect of parameters (raster width, raster angle, and layer thickness) and their interaction on tensile strength.However, experimental data showed that tensile strength and stiffness decreased with the increase of layer thickness, in spite of raster angle.Therefore, in the authors' view, discussion of the influence of parameters on FDM process without specifying the same standard/condition has minor significance for users and manufacturers, which may lead to different or even opposite conclusions.

Comparative Confirmation Test
Based on experimental data, the following regression models are developed, showing the interactions between the proposed independent variables: The objective of the coded regression model is to properly describe the interaction of factors affecting the mechanical strength at the concentration ranges analyzed.The measured tensile strength and compressive strength vary between 30.97 MPa-38.07Mpa and 19.23 Mpa-25.47Mpa, respectively, which agree with the desired values as presented in Figure 8.
Max T S (or C S ) Machines 2022, 10, x FOR PEER REVIEW 14 of 17 ( ) . .Optimization of the experimental data is performed by making mechanical strengths be the maximum and other coded variables in range.The optimum combinations obtained from formula (6) are A 0.27 B 1 C 1 D -1 for tension and A 1 B -1 C 1 D -1 for compression, respectively.In fact, A 1 B -1 C 1 D -1 just corresponds to the existing experiment of No. 6 in Table 6, whose value 25.2 MPa is the second largest compressive strength, a little lower than the maximum value 25.47 MPa of No. 5. Considering experiment and measurement errors, the predicted compressive strength can be approximately deemed as the best.In comparison, for tensile strength, the optimum settings A 0.27 B 1 C 1 D -1 correspond to no existing experiment, as 0.27 is an intermediate value.Based on the optimized data presented in Table 8, the experiment was conducted to validate the prediction.Furthermore, 38.7 MPa of tensile strength is achieved, which is the maximum compared with current value in Table 6.This indicates a good agreement of the predicted and experimental values under optimum conditions.In addition, a comparison of the results with previously conducted research shows that the performance of the optimized samples is further improved by RSM than Taguchi method, in terms of both tensile strength (38.7 MPa vs. 37.9 MPa) and compressive strength (25.2 MPa vs. 23.3MPa), indicating that RSM is a superior tool for optimizing FDM process parameters.As stated before, the regression model of RSM contains more terms than Taguchi, which facilitates fitting the characteristic curve more accurately.Moreover, unlike Taguchi method, the optimal solution of RSM can achieve any intermediate value in constraint.As a result, the optimum combination obtained by RSM is more likely to perform better, as demonstrated by the experiment.

Conclusions and Summary
In this work, the functional relationships between process parameters and mechanical strengths (tensile and compressive) for FDM process have been determined using the Taguchi method and RSM, respectively.The process parameters investigated are extrusion temperature, print speed, layer thickness and raster width.From this research, the following conclusions can be drawn: (1) The Taguchi method offered 27 runs for design of experiment while RSM suggested 30 runs.Therefore, the Taguchi method can decrease the number of tests compared to RSM, even considering interaction terms.However, limited by the form of orthogonal array, Taguchi method can only consider a specific number of interaction items besides linear terms (e.g., In L 27 (3 13 ), only three two-factor interactions can be investigated).In contrast, RSM can consider all interaction terms as well as quadratic terms.Therefore, the mathematical model for RSM performs better to fit the observed data.
(2) The results of significant ranking of FDM process parameters drawn by Taguchi method and RSM were different.Since most experimental data of parameter settings are different for these two approaches, the regression models are generally different, and so are the important orders of factors.Therefore, it is more meaningful to find optimum level combination by different DOE methods than to evaluate the influence of process parameter on FDM products, whose result is affected by various factors.
(3) At optimized conditions, both the Taguchi method and RSM succeed in achieving better specimen performance in terms of tension and compression, with small prediction errors.The optimum parameters from each response are as follows: Tensile strength: Taguchi method with 37.9 MPa-205 • C extrusion temperature, 0.35 mm layer thickness, 40mm/s print speed and 0.6mm raster width; RSM with 38.7 MPa-218 • C extrusion temperature, 0.35 mm layer thickness, 60mm/s print speed and 0.4 mm raster width.
In addition, RSM performed superiorly due to obtaining higher mechanical strengths.In fact, the optimum result of Taguchi method is essentially the permutation of existing levels for various parameters, which may lead to local optimal solution.The optimum result of RSM is more likely to achieve global optimal solution because intermediate values can be obtained (e.g., A 0.27 B 1 C 1 D -1 for tension).
Thus, it can be concluded that both the Taguchi method and RSM are robust statistical tools for FDM experimental design and process optimization.On the one hand, the Taguchi method can qualitatively show the effect of factors and determine the best level combination of parameters via S/N analysis with fewer experimental runs.However, the optimized solution obtained may miss the best result.On the other hand, RSM describes problems more quantitatively by fitting the full quadratic model via ANOVA analysis.By considering all degrees and values of factors, the global optimal condition is more likely to be achieved.However, the experiment may be time-consuming in the case of large-scale parameters.

Figure 7 .
Figure 7. Relationship of experimental groups of Taguchi and RSM.

Table 1 .
FDM parameters and their levels.

Table 1 .
FDM parameters and their levels.

Table 4 .
ANOVA for mechanical strengths of Taguchi method.

Table 5 .
Confirmation experiment for Taguchi method.

Table 8 .
Confirmation experiment for RSM.