Issue 
Int. J. Simul. Multidisci. Des. Optim.
Volume 12, 2021
Simulation and Optimization for Industry 4.0



Article Number  30  
Number of page(s)  9  
DOI  https://doi.org/10.1051/smdo/2021031  
Published online  09 November 2021 
Research Article
Experimental analysis and optimization of mechanical properties of FDMprocessed polylactic acid using Taguchi design of experiment
^{1}
Team of Applied Physics and New Technologies, Department of Physics − Polydisciplinary Faculty, Sultan Moulay Slimane University, Beni Mellal, Morocco
^{2}
Systems Engineering Laboratory, University Sultan Moulay Slimane, BeniMellal, Morocco
^{3}
Polytech Annecy Chambery − Univ. Savoie, SYMME, 74000 Annecy, France
^{*} email: mohamed.abouelmajd@usms.ac.ma
Received:
1
March
2021
Accepted:
26
October
2021
Fused deposition modeling (FDM) is one of the most used additive manufacturing processes in the current time. Predicting the impact of different 3D printing parameters on the quality of printed parts is one of the critical challenges facing researchers. The present paper aims to examine the effect of three FDM process parameters, namely deposition velocity, extrusion temperature, and raster orientation on the bending strength, stiffness, and deflection at break of polylactic acid (PLA) parts using Taguchi design of experiment technique. The results indicate that the temperature has the highest impact on the mechanical properties of PLA specimens followed by the velocity and the orientation. The optimum composition offering the best mechanical behavior was determined. The optimal predicted response was 159.78 N, 39.92 N/mm, and 12.55 mm for the bending strength, bending stiffness, and deflection at break, respectively. The R^{2} obtained from analysis of variance (ANOVA) showed good agreement between the experimental results and those predicted using a regression model.
Key words: Additive manufacturing / Taguchi method / desirability function / FDM process optimization
© M. Abouelmajd et al., Published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Additive manufacturing (AM), also called 3D printing, refers to a new construction technology of threedimensional objects by depositing layers one over the others with high precision. AM enables the production of complex tools, prototypes, and models offering the benefit of reasonable cost and short manufacturing time compared to the traditional methods such as machining and molding [1]. The ability to fabricate parts faster, cheaper, and without tooling increases the use of AM in medical sectors, functional prototype development, automobile fields, toolmaking, and many more [2–5]. Currently, AM technology has been developed to manufacture solid parts from a wide variety of materials including ceramics, thermoplastic polymers, and metals [6,7].
Among many AM processes, fused deposition modeling (FDM) is one of the technologies most commonly used by consumers in the current time. The FDM process is relatively easier to configure and use compared to the other AM technologies available on the market today. Also, the desktop 3D printers based on FDM process offer the advantage of building parts with short manufacturing duration, low cost, and important flexibility [8]. The small cost and the simplicity of operation make the home users able to prepare and fabricate parts of any idea easily.
The FDM process uses many types of thermoplastic polymers such as polylactic acid (PLA), Nylon, acrylonitrile butadiene styrene (ABS), polycarbonate, and others [9]. The majority of desktop 3D printers based on the FDM technique widely use PLA in the fabrication process. The PLA is a biodegradable polymer made from lactic acid derived from many recycled resources. It has a relatively low melting temperature that is between 190 °C and 220 °C [10,11].
Despite the wide number of FDM process benefits, it has a major disadvantage that is the process parameters affect the mechanical performance of printed pieces [12]. Also, the main part of home users has no way to examine the quality of their products. The influence of FDM process settings on the mechanical properties requires specific tests and analysis. Therefore, to make the 3D printed parts stronger and useful in functional fields, it is necessary to determine the optimal FDM process settings.
Many study were conducted in the last years to optimize the compressive performance [13,14], the tensile properties [15,16] and the flexural behavior of PLA built parts [17,18]. Moreover, this study investigates the characteristics of the FDM process parameters and their effects on the mechanical performance of printed parts. Three important mechanical properties namely limit load, bending stiffness, and flexural deflection at break characterize the mechanical response. The study was performed using Taguchi design to examine the effect of FDM process parameters on the response. The analysis of variance (ANOVA) technique was used to determine the percentage of contribution of each parameter in the resulting overall influence. The optimization of process parameters to improve the mechanical properties of printed parts was conducted using the composite desirability function.
2 Experimental procedure
Based on previously conducted research studies, three FDM process parameters were taken into consideration. These parameters are defined as follows:

Deposition velocity (V): it is the displacement speed of the extruder during the manufacturing period. The speed depends on the type of 3D printer machine.

Extrusion temperature (T): refers to the heating nozzle temperature, which ensures the transformation of thermoplastic polymer into semiliquid state. It is selected based on the material used and the printed part shape.

Raster orientation (O): also called printing angle, it represents the angle formed by the deposited filament and Xaxis of the build platform.
Three levels of each parameter were chosen according to the equipment used in this study and the past research works. In order to investigate the effect of FDM process factors on the response of printed parts and extract the response from the fewest number possible of tests, the experiments were conducted according to Taguchi design of experiments method. The FDM process parameters represent the input factors of the design, while the output is characterized by the mechanical response. Table 1 shows the factors considered and their respective levels. The other FDM process parameters such as infill density, infill pattern, build orientation, layer thickness, air gap and bed temperature were kept at their fixed level as mentioned in Table 2.
As shown in Figure 1, threedimensional models of parallelepiped parts of dimensions 80 mm long, 10 mm wide and 4 mm thick are modeled using a conception aided design (CAD) software. The parts dimensions are chosen according to the ISO 178:2019 standard [19], which describes the experimental procedure followed to determine the bending properties of thermoplastics. The test samples were printed according to the 3D models from the PLA material using WANHAO Duplicator 4S printer machine (see Fig. 2). Figure 3 shows part of the samples used in the mechanical test.
Flexural properties of structures are determined from the results of threepoint bending test. This mechanical test is commonly used because of the simplicity and the inexpensiveness of its procedure. As shown in the schema given in Figure 4, a polymer structure is deposited on two supports and a load is applied at the center of the structure. Every PLA specimen was tested under bending load until the fracture to determine its mechanical response. In this study, the threepoint bend test was achieved using Instron 5569 test machine with a crosshead velocity of 2 mm/min, and a support span (D) of 60 mm. In addition, the test machine has a maximum load cell of 10 KN.
The experiments were performed conforming to the orthogonal array of Taguchi Design L9 as shown in Table 3. In order to obtain a reasonable estimate of experimental error, three PLA specimens were tested for each test run and the mean values were taken into account.
Factors and Their Selected Levels.
Constant FDM Process Parameters.
Fig. 1
Dimensions of PLA flexural specimen. 
Fig. 2
The 3D printer machine WANHAO DUPLICATOR 4S. 
Fig. 3
Part of the samples used in mechanical test. 
Fig. 4
Schema of three point bending test. 
Orthogonal Array of Taguchi Design L9.
3 Results and discussion
3.1 Bending test results
The testing machine was equipped with an electronic extensometer to determine the applied load in the middle of PLA specimens during the mechanical test. The deflection resulting was measured by a displacement sensor. The data recorded from the threepoint bending test machine was collected using a computer system to plot the distribution of the applied load as a function of resulting deflection.
The conditions presented in Table 3 were taken into account to fabricate three samples for each test run. Thus, each test run was composed of three samples with the same printing conditions. The samples were tested under bending load, and the mean values between three specimens of the same printing settings were calculated. The LoadDeflection curves of the nine experiments are shown in Figure 5.
Three important mechanical properties have been determined from the LoadDeflection graphs, namely the limit load, bending stiffness, and maximum bending deflection. The mechanical properties are defined as follows:

Limit Load: characterize the PLA specimen strength; it is the maximal load supported before the fracture.

Bending Stiffness: it is the slope of the LoadDeflection curve in the domain of elastic deformations.

Deflection at break: refers to the maximum deflection subjected by the specimens.
Taguchi analysis was carried out to investigate the main effect of FDM process factors on the average mechanical properties of PLA parts. For the signaltonoise S/N ratio, the option larger is better was selected for the response that was to be maximized. In addition, the analysis of variance (ANOVA) was used to identify the statistical significance and the importance of the process parameters' effect on the mechanical behavior of specimens. For a confidence level of 95%, if the pvalue is less than 0.05, the effect of the corresponding factor is statistically significant.
The fracture surfaces of the specimens used in the bending test are illustrated in Figure 6. The surface of each sample is composed of two zones, a white zone formed by a section solicited in tension and a brown part that represents the compressed face. The observation of the fracture facets shows the ductile behavior of the different test pieces. The rapid rupture is produced after progressive damage of the specimens presented by the propagation of cracks, as indicated by the combination of rough and smooth fracture surfaces. The upper part of the specimen during the test was rigorous in comparison with the lower part which is relatively smooth.
Fig. 5
LoadDeflection curves of PLA structures under bending test. 
Fig. 6
Fracture surface of PLA samples. 
3.2 The limit load
A graphical comparison of limit load at the different experiments is depicted in Figure 7. In this figure, the average limit load of PLA specimens fabricated with the same settings was determined for each test run. The main effect plot for the means limit load based on Taguchi analysis is shown in Figure 8. The ANOVA results for limit bending load are listed in Table 4.
The comparison of the minimum and maximum average values of limit load showed deviation between PLA parts manufactured with different printing settings in the range of [160.68; 149.76] N. The limit load values obtained were in accordance with those found by Kumar et al. and Aloyaydi et al. [20,21]. The maximal value was obtained during the first test, while specimens used in Test 9 have the lowest average limit load. From the main effects of process parameters shown in Figure 8, it can be seen that the limit load decreases with decreasing deposition velocity, extrusion temperature, and raster orientation. Based on Taguchi analysis, the optimized process parameters for achieving maximal limit load are 30 mm/s deposition velocity, 190 °C extrusion temperature, and 0 degree raster orientation.
The ANOVA results show that the three factors considered have a statistically significant effect on the limit load (p < 0.05). In addition, it shows that the extrusion temperature has the highest impact on the limit load with a Percentage of Contribution (PoC) of 67%, followed by deposition velocity influence of 21%, while the lowest effect is characterized by the raster orientation with 11%. The Rsquared and adjusted Rsquared values of ANOVA for the limit load, 99.82%, and 99.27% respectively, show a good fitting of the model to the experimental results.
Fig. 7
Limit bending load of PLA specimens. 
Fig. 8
Main effect for means of limit load plot. 
ANOVA Analysis For Limit Bending Load.
3.3 The bending stiffness
The results of bending stiffness of the PLA products are presented in Figure 9. The average bending stiffness obtained ranges from 37.85N/mm to 40.19 N/mm. These results were in line with those reported in the literature for the PLA material manufactured using fused deposition modeling technology and under different printing conditions [22,23]. The maximum bending stiffness was reached at Test 1 and the minimum stiffness presents in specimens used in Test 8. It observed from the main effect plot for means bending stiffness (Fig. 10) that the lowest deposition velocity and extrusion temperature, 30 mm/s and 190 °C respectively, and the high level of raster orientation (60 degree) were the combination of the optimal levels to improve the bending stiffness of PLA parts.
The ANOVA for the bending stiffness was calculated and summarized in Table 5. The pvalues obtained from ANOVA analysis indicate that the impact of the three process parameters considered was statistically significant with a pvalue less than 0.05. Moreover, the extrusion temperature shows the highest influence on bending stiffness with a contribution ratio of 72%, followed by raster orientation and deposition velocity at reduced PoC of 14% and 13%, respectively. The high values of Rsquared (greater than 99%) show that the data fits adequately.
Fig. 9
Bending stiffness of PLA specimens. 
Fig. 10
Main effect for means of bending stiffness plot. 
ANOVA Analysis For Bending Stiffness.
3.4 The maximum deflection
The comparison of the average deflection at break of three repeated experiments during each test is depicted in Figure 11. It is clear from this figure that the mean deflection at break was raging between 8.39 mm obtained at Test number 3 and 13.85 mm achieved during the experiments of Test 4. These results were in agreement with the deflection at break reported in previous studies [24,25]. From the main effect plot for deflection at break shown in Figure 12, it is observed that the bending deflection at break will be maximal for deposition velocity of 70 mm/s, extrusion temperature of 190 °C, and raster orientation of 0 degree.
The ANOVA results presented in Table 6, show that the factors considered have a pvalue less than 0.05, which implies the significance of their effect on the maximum bending deflection. Also, it shows that the raster orientation and the extrusion temperature have an important impact on the maximum bending deflection with a contribution ratio of 38% and 37% respectively, followed by the influence of deposition velocity on bending deflection with 24% of contribution in the total deviations produced. The Rsquared values obtained from the ANOVA table, 99.66% of Rsquared, 98.65% of adjusted Rsquared indicate that the regression ANOVA model excellently fits the experimental data.
Fig. 11
Deflection at break of PLA specimens. 
Fig. 12
Main effect for means of maximum deflection plot. 
ANOVA Analysis For Deflection at Break.
3.5 Composite optimization
Taguchi analysis used to optimize the FDM process parameters remains a limited method. Each optimal combination obtained through this method will only improve just the corresponding mechanical property. For this reason, the composite desirability function was applied to perform multiobjective optimization. The composite desirability function is a suitable method to determine the optimal combination that maximizes all mechanical responses simultaneously together.
The approach based on the desirability function involves converting estimated response models (mechanical property) into individual desirability functions (d), which are then collected into a composite desirability function (D). The variation range of individual desirability (d) was from 0 to 1. The optimization results represent a completely undesirable value for a zero desirability (d = 0), while d = 1 represents a completely desirable response value. Based on the response specifications, three options of individual desirability functions are available:

Smaller is better: for a purpose to minimize the response.

Larger is better: for an objective to maximize the response.

Target: a specific value of response was the best.
In this paper, we intend to maximize the three mechanical responses (limit load, bending stiffness, and bending deflection at break). Thus, the second option of the individual desirability function (Larger is better) was chosen. To maximize the mechanical properties, the individual desirability (d) was calculated using equation (1).(1)
where y is the found value of the mechanical response, y_{low} and y_{high} are the minimum and maximum response values respectively, and w is the weight variable. Since the limit load, bending stiffness, and bending deflection at break are equally important, the value of weight variable w was taken equal to 1. The composite desirability function (D), as the geometric mean of the individual desirability functions (d_{i}), was calculated using equation (2).(2)here, n represents the number of responses.
In this study, the desirability function combined three mechanical responses into one composite desirability function (D), which is given by equation (3).(3)
The composite optimization for FDM parameters and the predicted mechanical properties are given in Figure 13. The dashed lines, blue, and horizontal represent the optimum levels that will maximize each mechanical response individually. An accordance was observed between the optimization obtained through the individual desirability function and the Taguchi analysis. The same optimal combinations were determined by the two approaches. The solid, red, and vertical lines indicate the optimum levels for the limit load, bending stiffness, and bending deflection at break simultaneously.
As shown in Figure 13, the optimum combination that will improve the three mechanical properties considered simultaneously was [30 mm/s, 190 °C, 30 °] for deposition velocity, extrusion temperature, and raster orientation respectively. The optimal predicted response of the limit load, bending stiffness, and deflection at break was 159.78 N, 39.92 N/mm, and 12.55 mm, respectively. The optimization was performed with excellent composite desirability that is relatively close to 1 (D = 0.85).
Fig. 13
Composite optimization plot. 
4 Conclusion
In this experimental investigation, efforts are made to investigate the influence of three FDM process parameters such as deposition velocity, extrusion temperature, and raster orientation on three important mechanical properties, namely limit load, bending stiffness, and bending deflection at break using Taguchi design of experiment method and analysis of variance technique. The deposition velocity is taken equal to 30 mm/s, 50 mm/s, and 70 mm/s. Three values of extrusion temperature were studied such as 190 °C, 200 °C, and 210 °C. Three levels of raster orientation were considered which are 0°, 30°, and 60°. The mechanical proprieties of PLA samples were determined from the experimental results of threepoint bending test. The Analysis of Variance was performed to assess the percentage of contribution of printing parameters in the effect on the mechanical behavior of printed parts. The ANOVA results show that the data fits adequately.
In addition, the ANOVA results show that the contribution ratio was related to the considered mechanical property. Generally, the extrusion temperature has the most important contribution ratio in the influence on the mechanical performance of printed structures. The optimization results indicate that a velocity of 30 mm/s, temperature of 190 °C, and orientation of 0° gives the maximum bending strength. On other hand, the highest bending stiffness was found at 30 mm/s of velocity, 190 °C of temperature, and 60° of orientation. The combination formed by a velocity of 70 mm/s, temperature of 190 °C, and orientation of 30° shows the possibility of producing parts with a high deflection at break. It is observed that the optimal printing parameters are determined according to the mechanical property. The optimal process parameters combination for the limit load was not necessarily the combination of printing settings that allows us to improve another mechanical property. Therefore, the multiobjective optimization was performed using the desirability function. The optimum printing parameters combination for improving the mechanical properties of PLA parts simultaneously was found to be [30 mm/s, 190 °C, 30 °].
References
 M. Chaudhari, B.F. Jogi, R. Pawade, Comparative study of part characteristics built using additive manufacturing (FDM), Proc. Manufactur. 20, 73–78 (2018) [Google Scholar]
 M.K. Thompson, G. Moroni, T. Vaneker, G. Fadel, R.I. Campbell, I. Gibson et al., Design for additive manufacturing: trends, opportunities, considerations, and constraints, CIRP Ann. 65, 737–60 (2016) [Google Scholar]
 J.S. Mohammed, Applications of 3D printing technologies in oceanography, Methods Oceanogr. 17, 97–117 (2016) [Google Scholar]
 Z. Liu, Y. Wang, B. Wu, C. Cui, Y. Guo, C. Yan, A critical review of fused deposition modeling 3D printing technology in manufacturing polylactic acid parts, Int. J. Adv. Manufactur. Technol. 102, 2877–89 (2019) [Google Scholar]
 M. Abouelmajd, A. Bahlaoui, I. Arroub, M. Lagache, S. Belhouideg, Mechanical characterization of PLA used in manufacturing of 3D printed medical equipment for COVID19 pandemic, 2020 IEEE 2nd International Conference on Electronics, Control, Optimization and Computer Science (ICECOCS), IEEE (2020), pp. 1–5 [Google Scholar]
 I. Gibson, D. Rosen, B. Stucker, M. Khorasani, Materials for Additive Manufacturing. Additive Manufacturing Technologies (Springer; 2021), pp. 379–428 [Google Scholar]
 J. Gardan, Smart materials in additive manufacturing: state of the art and trends, Virtual Phys. Prototyp. 14, 1–18 (2019) [Google Scholar]
 V.E. Kuznetsov, A.N. Solonin, O.D. Urzhumtsev, R. Schilling, A.G. Tavitov, Strength of PLA components fabricated with fused deposition technology using a desktop 3D printer as a function of geometrical parameters of the process, Polymers 10, 313 (2018) [Google Scholar]
 D. Bourell, J.P. Kruth, M. Leu, G. Levy, D. Rosen, A.M. Beese et al., Materials for additive manufacturing, CIRP Ann. 66, 659–81 (2017) [Google Scholar]
 M, Van den Eynde, P. Van Puyvelde, 3D Printing of Poly (lactic acid). Industrial Applications of Poly (lactic acid) (2017), pp. 139–58 [Google Scholar]
 X. Zhou, S.J. Hsieh, Y. Sun, Experimental and numerical investigation of the thermal behaviour of polylactic acid during the fused deposition process, Virtual Phys. Prototyp. 12, 221–33 (2017) [Google Scholar]
 J.R.C. Dizon, A.H. Espera Jr, Q. Chen, R.C. Advincula, Mechanical characterization of 3Dprinted polymers, Addit. Manufactur. 20, 44–67 (2018) [Google Scholar]
 S. Brischetto, R. Torre, Tensile and compressive behavior in the experimental tests for PLA specimens produced via fused deposition modelling technique, J. Compos. Sci. 4, 140 (2020) [Google Scholar]
 H.K. Dave, S.R. Rajpurohit, N.H. Patadiya, S.J. Dave, K.S. Sharma, S.S. Thambad et al., Compressive strength of PLA based scaffolds: effect of layer height, infill density and print speed, Int. J. Mod. Manuf. Technol. 11, 21–27 (2019) [Google Scholar]
 Y. Zhao, Y. Chen, Y. Zhou, Novel mechanical models of tensile strength and elastic property of FDM AM PLA materials: experimental and theoretical analyses, Mater. Des. 181, 108089 (2019) [Google Scholar]
 T. Yao, Z. Deng, K. Zhang, S. Li, A method to predict the ultimate tensile strength of 3D printing polylactic acid (PLA) materials with different printing orientations, Compos. B 163, 393–402 (2019) [Google Scholar]
 C. Abeykoon, P. SriAmphorn, A. Fernando, Optimization of fused deposition modeling parameters for improved PLA and ABS 3D printed structures, Int. J. Lightweight Mater. Manufact. 3, 284–97 (2020) [Google Scholar]
 S.R. Rajpurohit, H.K. Dave, Flexural strength of fused filament fabricated (FFF) PLA parts on an opensource 3D printer, Adv. Manufactur. 6, 430–41 (2018) [Google Scholar]
 I. Standard, ISO B. Plastics—Determination of flexural properties, ISO Geneva, Switzerland (2019) [Google Scholar]
 S. Kumar, R. Singh, M. Singh, T. Singh, A. Batish, Multi material 3D printing of PLAPA6/TiO_{2} polymeric matrix: Flexural, wear and morphological properties, J. Thermoplast. Compos. Mater. 0892705720953193 (2020) [Google Scholar]
 B.A. Aloyaydi, S. Sivasankaran, H.R. Ammar, Influence of infill density on microstructure and flexural behavior of 3D printed PLA thermoplastic parts processed by fusion deposition modeling, AIMS Mater. Sci. 6, 1033–48 (2019) [Google Scholar]
 J. Chacón, M.A. Caminero, E. GarcíaPlaza, P.J. Núnez, Additive manufacturing of PLA structures using fused deposition modelling: effect of process parameters on mechanical properties and their optimal selection, Mater. Des. 124, 143–57 (2017) [Google Scholar]
 J. Porter, T. Cain, S. Fox, P. Harvey, Influence of infill properties on flexural rigidity of 3Dprinted structural members, Virtual Phys. Prototyp. 14, 148–59 (2019) [Google Scholar]
 Z. Liu, Q. Lei, S. Xing, Mechanical characteristics of wood, ceramic, metal and carbon fiberbased PLA composites fabricated by FDM, J. Mater. Res. Technol. 8, 3741–51 (2019) [Google Scholar]
 S.Z. Gebrehiwot, L.E. Leal, J. Eickhoff, L. Rechenberg, The influence of stiffener geometry on flexural properties of 3D printed polylactic acid (PLA) beams, Progr. Additive Manufactur. 6, 71–81 (2021) [Google Scholar]
Cite this article as: Mohamed Abouelmajd, Ahmed Bahlaoui, Ismail Arroub, Maria Zemzami, Nabil Hmina, Manuel Lagache, Soufiane Belhouideg, Experimental analysis and optimization of mechanical properties of FDMprocessed polylactic acid using Taguchi design of experiment, Int. J. Simul. Multidisci. Des. Optim. 12, 30 (2021)
All Tables
All Figures
Fig. 1
Dimensions of PLA flexural specimen. 

In the text 
Fig. 2
The 3D printer machine WANHAO DUPLICATOR 4S. 

In the text 
Fig. 3
Part of the samples used in mechanical test. 

In the text 
Fig. 4
Schema of three point bending test. 

In the text 
Fig. 5
LoadDeflection curves of PLA structures under bending test. 

In the text 
Fig. 6
Fracture surface of PLA samples. 

In the text 
Fig. 7
Limit bending load of PLA specimens. 

In the text 
Fig. 8
Main effect for means of limit load plot. 

In the text 
Fig. 9
Bending stiffness of PLA specimens. 

In the text 
Fig. 10
Main effect for means of bending stiffness plot. 

In the text 
Fig. 11
Deflection at break of PLA specimens. 

In the text 
Fig. 12
Main effect for means of maximum deflection plot. 

In the text 
Fig. 13
Composite optimization plot. 

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