Next Article in Journal
Exploring Alternatives to Polyacrylamide: A Comparative Study of Novel Polymers in the Flocculation and Dewatering of Iron Ore Tailings
Next Article in Special Issue
MC-Injection Molding with Liquid Silicone Rubber (LSR) and Acrylonitrile Butadiene Styrene (ABS) for Medical Technology
Previous Article in Journal
Tetramine-Based Hyperbranched Polyimide Membranes with Rigid Crosslinker for Improved Gas Permeability and Stability
Previous Article in Special Issue
Development of an Injection Mold with High Energy Efficiency of Vulcanization for Liquid Silicone Rubber Injection Molding of the Fisheye Optical Lens
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integration of Multivariate Statistical Control Chart and Machine Learning to Identify the Abnormal Process Parameters for Polylactide with Glass Fiber Composites in Injection Molding; Part I: The Processing Parameter Optimization for Multiple Qualities of Polylactide/Glass Fiber Composites in Injection Molding

Department of Materials Science and Engineering, National Taiwan University of Science and Technology, Taipei 10607, Taiwan
*
Author to whom correspondence should be addressed.
Polymers 2023, 15(14), 3018; https://doi.org/10.3390/polym15143018
Submission received: 1 May 2023 / Revised: 8 June 2023 / Accepted: 4 July 2023 / Published: 12 July 2023
(This article belongs to the Special Issue Injection Molding Process in Polymer Processing)

Abstract

:
This paper discusses the mixing of polylactide (PLA) and glass fiber which use injection molding to produce a functional composite material with glass fiber properties. The injection molding process explores the influence of glass fiber ratio, melt temperature, injection speed, packing pressure, packing time and cooling time on the mechanical properties of composite. Using the orthogonal table planning experiment of the Taguchi method, the optimal parameter level combination of a single quality process is obtained through main effect analysis (MEA) and Analysis of variance (ANOVA). Then, the optimal parameter level combination of multiple qualities is obtained through principal component analysis (PCA) and data envelopment analysis (DEA), respectively. It is observed that if all the quality characteristics of tensile strength, hardness, impact strength and bending strength are considered at the same time, the optimal process conditions are glass fiber addition 20 wt %, melt temperature 185 °C, injection speed 80 mm/s, holding pressure 60 MPa, holding time 1 s and cooling time 15 s, and the corresponding mechanical properties are tensile strength 95.04 MPa, hardness 86.52 Shore D, impact strength 4.4408 J/cm2, bending strength 119.89 MPa. This study effectively enhances multiple qualities of PLA/GF composite.

1. Introduction

Nowadays, the importance of improving component mechanical and thermal qualities, incorporating them into production systems, and employing eco-friendly materials is rising [1]. On the contrary, every industry must manage its environmental impact from conception to final disposal. A material that is thermo-mechanically durable, lightweight, and biodegradable is required. The demand for poly(lactic acid) (PLA) composites reinforced with natural and synthetic fibers has surged as a result of this sectoral development [2].
PLA can be biodegradable, low-pollution and non-toxic, high mechanical strength, biocompatibility, application in the fields of medicine, packaging, and daily necessities, which can greatly reduce the impact of resources on environmental damage, and become a new type of polymer material with great development potential [3]. PLA polymers can be produced through direct lactic acid poly-condensation and also via ring-opening polymerization of lactide, a cyclic dimer of lactic acid [4]. Because the PLA products produced by traditional processing technology have a slow crystallization rate and low crystallinity in the molding process. PLA also has poor heat resistance, with a heat deflection temperature (HDT) of 55–65 °C, which seriously limits the application range of PLA at higher temperatures [5], such as disposable heat-resistant products, utensils and other food packaging containers [6]. PLA stands out with its eco-friendly properties so it has been recently more preferable to synthetic polymers in some sectors. There are some issues such as poor heat resistance of PLA, poor crystallization rate, deformation and material brittleness of disposable tableware during transportation, so, there is an urgent need to modify PLA to improve its heat resistance and all other properties as well. The use of biodegradable polymer PLA is limited, so it is necessary to enhance its properties through composite reinforcements.
Glass fiber (GF) is effectively to modify the substrate, and enhance the mechanical properties as well. After recycling, it can be reused and meets environmental protection requirements. It is a well-known reinforcement material in the composites sector as well [7]. Since it comprises inorganic components, it has excellent dimensional stability, transparency, and mechanical appropriateness [8]. SiO2, Al2O3, Fe2O3, CaO, MgO, Na2O, K2O, and BaO could be identified as the primary chemical constituents of glass fiber [9]. Because GF has a low thermal conductivity, it can dissipate heat similarly to asbestos and other organic fibers [10]. This makes it a popular insulation material in many industries. GF reinforced plastics continue to be the most prominent and dominant material in the industry, despite advancements and innovations in the plastic-composite field.
Wang et al. [11] showed not only GF but also enhanced crystallization led to the outstanding mechanical performance of PLA/GF composites. While GF shows heat treatment can remarkably improve thermal stability, in particular for PLA/GF composites. Sun et al. [12] presented that the defect in fiber GF reinforced PLA is significantly improved by modification and the mechanical properties are improved by about 40%. The main reason is the unification of the surface polarity of fibers and PLA, as well as the connection established by the functional groups. Meanwhile, the surface modification of GF/PLA composites can also improve their thermal and degradation properties. Wang et al. [13] demonstrated that the outstanding mechanical properties arises from the strengthening effect of the GF network skeleton that shows good bonding with PLA matrix. GF led to simultaneously enhanced strength, rigidness and toughness of PLA. Thermal analysis showed that GF led to increased heat deflection temperature of PLA. Leu et al. [14] used injection molding for PLA/maleated styrene-ethylene/butylene-styrene/organo-montmorillonite to improve the mechanical properties of composite materials. Ma et al. [15] discussed preparation and foaming extrusion behavior of PLA/polybutylene succinate (PBS)/montmorillonoid nanocomposite. The compatibility between PLA and PBS, and the elongation at break and impact strength of composite materials can be improved. Shen et al. [16] prepared biocomposite of PLA/reinforced hydroxyapatite (HA)/carbon fiber from hot pressing a prepreg which was manufactured by solvent impregnation process in order to improve composite mechanical properties. Bledzki et al. [17] used injection molding to make tensile and impact test pieces to verify the improved mechanical properties of composites of PLA man-made cellulose and abaca fibers. An et al. [18] investigated the cutting characteristics of GF reinforced plastics with respect to tool materials and geometries to improve tensile strength, impact and flexural strength, corrosion resistance and non-conductive properties. Bigg et al. [19] presented that GF can be effective increase the mechanical properties of materials such as strength and stiffness in thermoplastic composite materials. Jaszkiewicz et al. [20] applied GF/abaca fiber/man-made cellulose to prepare composite materials with PLA and polypropylene (PP) respectively, showing that GF can enhance the impact strength of PLA and PP.
PLA has excellent processing performance and can be processed by extrusion, film blowing and injection molding [21,22,23]. The extrusion process has some disadvantages because of nozzle radius limits, reduces the final quality, limits the accuracy and speed when compared to other processes, consistent pressure of material is required in order to increase quality of finish. In film blowing method, it is difficult to control accurately the thickness of a blown film and method is quite complicated as well, and there are a number of factors that can go wrong, the manufacturing cost for blown films is high and not environment-friendly.
Oppositely, the injection molding is advantageous as compared to other techniques because its low-waste process, it minimizes molding costs, highly repeatable way of production with high precision. Injection molding can produce a huge amount of parts per hour from a wide range of other materials, injection molding technology can limit the waste by recycling wherever possible, planning production runs to maximize efficiency, and conserving energy. It has become the most important production technology in polymer plastics and composite plastic materials [24,25]. The materials and process parameters are the important factors affecting product quality.
Therefore, this paper will discuss the mixing of PLA and glass fiber which use injection molding to produce a functional composite material with glass fiber properties. The influence of processing parameters will be discussed on various qualities through the injection molding process. The orthogonal table in Taguchi method will be used to plan the experiment. Through the MEA and ANOVA to obtain the process optimization parameters of a single quality. In response to the multiple quality characteristics of this study, PCA is applied to reduce the dimensionality of the relevant quality characteristics into independent linear combinations, and DEA method calculates the objective optimal weight of the original data to obtain relative efficiency. Afterwards, the optimal combination of processing parameter factor levels will be found and the confirmation experiments will be conducted to verify the optimized results. PLA based glass fibers that are added in higher content to produces a desirable characteristic so that the treated fibers produced desirable reinforcement effects.

2. Experimental Methodology

2.1. Injection Molding Process

Plastic injection molding is one of the most widely used plastic fabrication processes for plastic mass production with numerous shapes and complicated geometries. It has preliminarily been estimated that over 30% of the polymers that are processed as well as consumed are produced by the injection-molded process [26,27].
In the injection molding manufacturing, while there are a number of parameters that must be determined, some have been recognized as the important process parameters in relation to product quality. As the most popular plastic molding processing method at present, Bozzelli [28] proposed that melt temperature, injection pressure, injection speed and cooling time are important factors for plastic thin shell injection molding. Jansen et al. [29] pointed out that the impact on the shrinkage, the biggest ones are melt temperature, holding pressure and injection speed. Shokri et al. [30] showed that the properties of fiber-reinforced thermoplastic injection molding products depend on the influence of packing pressure on fiber orientation. Kamaruddin et al. [31] presented that melt temperature, high injection pressure, low packing pressure, long holding time and cooling time can effectively reduce the shrinkage behavior of injection molded products.
Kuo et al. [32] indicated that cooling time, mold temperature, melt temperature, injection speed, injection pressure, packing pressure, and packing time are the key factors for plastic LCD light-guide plates.
The related research concerning about the process factors include the mold temperature, the melt temperature, the packing pressure, the packing time and cooling time [26]. The current manufacturing application determining the injection molding process parameters involves a combination of the use of the machine operation handbook and accompany with the adaptations through trial and error from experienced plastic engineers [26]. In order to guarantee that the optimal process parameters have been selected, the demand to establish these optimal parameters has given rise to this research.

2.2. Process Optimization

In traditional experiments, when the process parameters increase, the number of experiments will increase. In order to solve this problem, Karna et al. [33] used the Taguchi method, the robust design of the orthogonal table, the S/N ratio and ANOVA to study the impact of process parameters on the product. Liu et al. [34] used Taguchi method to analyze the parameter optimization of thin shell parts in the injection molding process, showing that the melt temperature and injection pressure are the most important processing parameters. Ghani et al. [35] using the Taguchi method in the high-speed milling process, through the S/N ratio and ANOVA, the optimal process parameters are optimized. However, the Taguchi method is only suitable for the improvement of a single quality. In actual industry, it needs to be combined with other analysis methods to achieve the goal of multi-quality process optimization. For example, Su et al. [36] applied principal component analysis method to reduce the dimension and complexity and solved multi-quality problems. Antony [37] used the PCA method, combined with the quality loss function to effectively improve and take into account the effect of multi-quality. Shih et al. [38] presented the inert gas shielded welding process to weld the foamed aluminum plate. Taguchi method combined with the PCA method showing that the current, welding speed and the gap between the workpieces are important control factors in the process. The optimal parameters of the process could improve the multi-quality characteristics of the aluminum foam board. Jeyapaul et al. [39] aimed at the operation of the gear processing machine with six control factors. It showed that compared with Taguchi method, the genetic algorithm and DEA method are used for the optimal factor level combination S/N ratio of the qualities, and the total expected improvement is 4.1498 db and 11.2506 db, respectively. Al-Refaie et al. [40] studied the improvement of the quality of the hard disk drive with controllable factors. Compared with Taguchi method, when PCA method and DEA method used to optimize the quality process parameters, the total expected improvement of the optimal factor level combination S/N ratio are 4.1498 db and 11.2506 db, respectively
Therefore, this paper will use the Taguchi method, and combine with PCA and DEA to achieve the goal of the optimizing multiple qualities.

2.3. Materials

Manufacturer: Nytex Composites Co., Ltd. New Taipei City, Taiwan. Product number: GG-0010N (TY11512706, 10% Glass fiber), GG-0015N (TY11512707, 15% Glass fiber), GG-0020N (TY11512708, 20% Glass fiber).
The material properties are shown in Table 1.

2.4. Experimental Methodolofy

This section will introduce the injection molding, material analysis, possible reaction of the composite and experimental scheme.

2.4.1. Injection Molding

Injection molding machines perform a wide range of mechanical movements with differing characteristics. Mold opening is a low-force high-speed movement, and mold closing a high-force low-speed movement. Plasticizing involves high torque and low rotational speed, while injection requires high force and medium speed. Injection molding machine consist of three major components i.e., (1) Screw motor drive (2) Reciprocating screw and barrel, (3) Heaters, thermocouple, and ring plunger.
The operation principle of the injection molding is very simple, where plastic material is heated above its melting point, resulting in the conversion of the solid polymer to a molten fluid with a reasonably low viscosity. It is then forced into a closed mold that defines the shape of the article to be produced. The operation elements are shown in Figure 1. The injection samples are shown in Figure 2.
The plastic material from the feeding hopper enters the barrel, mixed by the screw, sent to the front end of the heating tube along the spiral groove, and is heated by the peripheral heater. The screw rotates to fully mix the plastic material so that the plastic is in a molten state. When the screw rotates, the screw retreats due to the reaction force (back pressure) of the plastic material. At this time, use the limit switch to constrain the amount of retreat, stop the screw rotation at a certain position, then close the mold into the injection stage. Meanwhile, the hydraulic cylinder of the injection device exerts injection force on the screw, and the screw becomes an injection plunger. Under high pressure, the completely melted plastic material at the front end of the barrel is injected into the mold from the nozzle. After the material in the cavity cools down, the mold is opened and eject the finished product. The injection molding machine can form plastic products with complex shapes, precise dimensions or dense texture with metal inserts at one time.

2.4.2. Materials Analysis

The instrument used is differential scanning calorimeter (DSC) for thermal properties. The instruments to measure mechanical properties such as tensile strength, Shore hardness, impact strength, and bending strength. The model used is MTS 810, the maximum displacement range: ±75 mm, the maximum test load: ±100 kN. Comply with ASTM D790 standard, observe the strength change of tension and bending. According to the ASTM D2204-00 standard, the composite material studied is more than 90 Shore A, using D type Shore hardness tester. The impact test is to determine the toughness of the material. The model of Izod impact testing machine used in this research is Yasuda Seiki N0158, which measures the impact energy of materials according to ASTM D256 standard.
The possible reaction between PLA and GF to synthesize PLA/GF composite is given in Figure 3.
A coupling agent versatile molecule, was employed to modify the fiber surface which generate a chemical bond between the siloxy group and the alkyl group. Silane coupling agents transformed fibers by a multi-step process that included bonding, condensation, and hydrolysis. The hydrolysis of siloxy groups resulted in the formation of silanol. The hydrophobicity of the molecule was increased by its ability to interact with the hydroxyl group of cellulose during the condensation process, and the opposite side of the molecule reacted with the PLA matrix to establish a bond (Figure 4). The enhancement of interfacial characteristics was credited for a boost in tensile strength and flexibility. Another purpose of silane is to serve as a surface protective layer by penetrating the pores of the fiber surface.

2.4.3. Scheme of Experiment and Processing

In this section, the material properties of the composite material are analyzed to set up the range of processing parameters. The L18 orthogonal table is used to plan the experiment, combined with DEA and PCA respectively to achieve the optimization of multiple qualities. Then, the optimized parameter combination is implemented in the confirmation experiment to verify the feasibility and reproducibility of the optimized parameters. The planning process of this experiment is shown in Figure 4.

3. Taguchi and Other Statistical Techniques

This study uses PCA and DEA to optimize the process parameters of PLA and GF composites used in injection molding machines.

3.1. Taguchi Method [41]

3.1.1. Orthogonal Table

The orthogonal table is expressed as La (bc) represents the orthogonal table, a is the column number (experiment times), b is the level number, and c is the row number.

3.1.2. Signal-to-Noise (S/N) Ratio

The quality discussed in this study is the larger the mechanical properties of tensile strength, Shore hardness, impact strength, and bending strength, the larger the better (LTB). The S/N ratio of the maximum characteristic is defined as:
S / N L T B = 10 l o g 10 ( 1 n i = 1 n 1 y i 2 )
where MSD is the mean square deviation from the target value, n is the total number of measurements, and y i is the quality measurement value.

3.1.3. Main Effects Analysis (MEA)

Find the average response value of each factor level and the main effect value ΔFi from the experimental data, and then make a response table for the MEA of each factor. When the main effect value of a factor is larger, it means that the factor has a greater impact on the system. On the contrary, the smaller it is, such as Equations (2) and (3).
F ¯ i = 1 m j = 1 m η j
Δ F = max { F ¯ 1 , F ¯ 2 , F ¯ 3 , , F ¯ n } min { F ¯ 1 , F ¯ 2 , F ¯ 3 , , F ¯ n }
where m is the number of level i in the factor row of the orthogonal table, ηj is the S/N ratio produced by each j level column, n is the level of the factor.

3.1.4. ANOVA

ANOVA analyzes the contribution of each factor to determine the importance of each factor:
I.
Degree of freedom (DOF)
(1)
degrees of freedom for each factor
DOF factor = ( level   number ) 1
(2)
total number of degrees of freedom
DOF total = n × r 1 = L 1
where n is the number of experimental groups, r is the number of repeated experiments, and L is the total number of experiments.
(3)
error degrees of freedom
DOF error = DOF total DOF factor
II.
Total sum of squares (SST), the total variation
SS T = i = 1 n ( η i η ¯ i ) 2 = i = 1 n η i 2 CF
where n is the total number of experimental observations, ηi is the S/N ratio of each group of experiments, and is the average of overall S/N ratio.
CF is the correction factor, defined as:
CF = 1 n i = 1 n ( η i ) 2
III.
The sum of squares of each factor (SS), the variation of each factor. if a factor has p levels, and each level has m observations, then the sum of squares is:
SS A = 1 m A 1 2 + A 2 2 + A 3 2 + + A p 2 CF
IV.
Error sum of squares (SSerror):
SS error = SS T SS factor
V.
Mean square, MS, the variance:
MS = SS DOF factor
VI.
Error mean square (MSE)
MSE = SS error DOF error
VII.
F-ratio indicates the relationship between the factor effect and the error variation. When the F value is larger, it means that the factor has a more important influence on the system, and it is used to arrange the important order of the factors.
F = MS MSE
VIII.
Pure sum of square (SS′)
SS factor = SS factor DOF factor × MSE
IX.
Percent contribution (ρ), the relative ability to reduce variation for factors.
ρ = SS factor SS T × 100 %

3.1.5. Confidence Interval (CI)

To evaluate each observation value effectively, it is necessary to calculate its CI.
C I S / N = F α ; 1 , V 2 × MSE × 1 n e f f + 1 r
where F α ; 1 , V 2 is the F value with a significant error α, v2 is the degree of freedom of the combined error variance, MSE is the combined mean square error, r is the number of confirmation experiments, and neff is the effective observation value.
n eff = total   number   of   experiments 1 + ( sum   of   degrees   of   freedom   for   factors   to   evaluate   the   mean )
Calculate the 95% confidence interval to verify the validity of the confirmed experimental mean under the predicted optimal parameter conditions, as sown in Equation (18).
S N C I S / N μ confirmation S N + C I S / N
where μ confirmation is the mean value of the confirmation experiment.
And
S N = T ¯ + i = 1 n ( F i T ¯ )
where T ¯ is the total average of S/N ratio, Fi is the S/N ratio of significant factor level.

3.2. PCA [36]

The steps to use PCA are described as follows:
Step 1. List the quality data of each group of experiments, and obtain the S/N ratio of its quality characteristics for PCA.
Step 2. Use Equation (20) to normalize the data of each quality characteristic, so that the data is between 0 and 1
X normalization = x i ( j ) min [ x i ( j ) ] max [ x i ( j ) ] min [ x i ( j ) ]
Step 3. The normalized data is obtained to obtain the correlation coefficient matrix
ρ xy = ( x i x ¯ ) ( y i y ¯ ) ( x i x ¯ ) 2 ( y i y ¯ ) 2
where ρ x y is the correlation coefficient of x to y, and x ¯ is the average value of item x.
Step 4. Use the correlation coefficient matrix to obtain its eigenvalues, which are the principal components, and the corresponding eigenvectors. The variation of the i-th principal component is shown in Equation (22).
ρ i = λ i i = 1 λ i
where ρ i is the variance of the i-th principal component in the total variation, and λi is the eigenvalue of the i-th principal component.
Step 5: Using the eigenvectors corresponding to the eigenvalues of the principal components and the normalized matrix X, the score of the principal components can be obtained from Equation (22).
Y i = XV i

3.3. Data Envelopment Analysis (DEA) [40]

DEA is a fractional mathematical programming technique for evaluating the relative efficiency of decision making unit (DMU) with multiple inputs and multiple outputs. It combines various inputs and various outputs for a DMU into one performance measure (called relative efficiency).

3.3.1. Charnes, Cooper and Rhodes (CCR) Input-Oriented Model

Based on the current output level, discuss how much “input” should be used to be an efficient DMU, and establish an evaluation model for DMUk:
Maximize h k = r = 1 s μ r y rk i = 1 m V i X ik
Subject   to r = 1 s μ r y rj i = 1 m V i X ij 1 ,   j = 1 , , n μ r 0 ;   r = 1 , , s v i 0 ;   i = 1 , , m
where μ r , v i are the weights of the r-th output item and the i-th input item, respectively.
Equation (25) indicates that the “output combination” of any DMU cannot be greater than its “input combination”.
Set i = 1 m V i X ik = 1 , Equations (24) and (25) can be changed to Equations (26) and (27).
Maximize   h k = r = 1 s μ r y rk
Subject   to i = 1 m V i X ik = 1 r = 1 s μ r y rj i = 1 m V i X ij 0 ,   j = 1 , , n μ r ,   v i 0 ,   r , i

3.3.2. Cross-Efficiency Analysis Model

The cross-evaluation measure was introduced by Sexton, et al. [42]. Let Eoj denotes the cross-efficiency of DMUj calculated according to the optimal weights of DMUo. For each Eoj, it is the (weighted output)/(weighted input) obtained by substituting the u ro * and v i o * corresponding to the o-th evaluated unit into the observed value of the j-th evaluated unit, as shown in Equation (28).
E oj = r = 1 s u ro * y rj i = 1 m v io * x ij E oj 1 ,   o j
This uses DEA in a peer-evaluation instead of a self-evaluation calculated by CCR model.
Let the mean of cross-efficiencies for DMUj expressed as:
e j = o j E oj ( n 1 ) j = 1 , , n
The ordinal value is to rank the ej values such that the smallest ej value obtains one whereas the largest ej value gets n.
Let A O V f g is the average of the ordinal values for level g of factor f . From calculating A O V f g value for each factor level. The optimal factor level, g * , is chosen as the level that maximizes the value of A O V f g , denoted by
g * = g | max g { AOV f g f
Cross-efficiency maximizes self-evaluation efficiency and minimizes peer-evaluation efficiency.

3.4. Materials Analysis

The DSC is used to measure the melting point of the composites’ material. The sample 2.0 mg is placed in the sample pan. The operating condition rises from 20 °C to 270 °C at a heating rate of 20 °C/min as shown in Figure 5, Figure 6 and Figure 7. The melting point of the material is about 152 °C, which is close to the melting point (155 °C) provided by the manufacturer, so that processing temperature of the composite material should not be lower than this melting point. Verify that the recommended injection temperature provided by the manufacturer is 170 °C~195 °C, which can be used as the melt temperature factor level in the orthogonal table. Kumar and Prakash [43] explained the DSC analysis of pure PLA and composites of PLA. They discussed the thermal characterizations of the composites. There were two peaks at 60.06 °C for glass transition temperature (Tg) and 147.71 °C melt temperature (Tm), with Delta values 0.6354 J/g and 28.2 J/g was observed for pure PLA as explained in literature. When these peak values observed in 20% PLA composites with glass fibers, it was increased to 68.69 °C and 152.35 °C with Delta values 11.387 J/g and 20.371 J/g. Overall, these results explained that PLA composites marks an enhanced thermal behavior and these results are consistent with the literature [4]. The use of other material to synthesize PLA composite raises the polymer breakdown temperature. The differential scanning calorimetry (DSC) curves showed the same behavioral properties as explained in present articles.

3.5. Injection Molding Process Parameter Selection

This project is to use the water circulation to cool the injection molding test sample mold. This cooling method is especially suitable for molds with simple shapes and can achieve a uniform cooling effect. By ensuring that the mold is cooled evenly, we can ensure that the quality and dimensions of the product meet the requirements. The parameters that affect the finished workpiece in the injection molding process are speed, temperature, pressure and time [26,28,29,30,31]. Because the speed affects the amount of cavity filler, the temperature affects the shear viscosity of the material, the pressure affects the volumetric shrinkage, and the time depends on the size of the injection molding equipment and the residence time of the material. RTP Company has confirmed that the glass fiber content reinforced polylactic acid compound improves the mechanical properties of polylactic acid. Refer to the machine operation handbook as well, so the glass fiber, so the glass fiber content, melt temperature, injection speed, holding pressure, holding time, and cooling time are set as the control parameters of the injection molding machine. Then the experiments were actually tried out, and the other levels that could result in deviations in the quality of composite material were tried to find, thereby identifying a suitable working range. Finally, for the composite material injection molding processing parameters, the factors that were actually controllable by the injection-molding machine were chosen.
When the temperature is lower than 175 °C, due to high viscosity by the incomplete melting of the material, the nozzle will be stuck. When the temperature is higher than 195 °C, the injection molded test piece will be coked and carbonized, so the processing temperature range is set at 175 °C~195 °C. The control factors and their levels of this experiment are as shown in Table 2.
In this study, the level value of the control factor was applied to the L18 (36) orthogonal table for experimental planning. Each group had five test pieces, a total of 90 experimental data. The MEA and ANOVA were used to obtain the optimal process parameters for each quality.

4. Experiment results

4.1. Experimental Data and Corresponding S/N Ratios

The results for the three iterations of the 18 experiments over 5 iterations with averages, and S/N ratios of five quality characteristics are shown in Table 3.

4.2. Single Quality Optimization Analysis

4.2.1. Tensile Strength Test Data Analysis

(1)
MEA
From the S/N ratio obtained from the experiment as shown in Table 3, the main effect of each control factor is calculated, and the response graph is drawn, as shown in Figure 8. It shows that the best factor level selection is A3 (glass fiber 20%), B2 (melt temperature: 185 °C), C1 (injection speed: 40 mm/s), D2 (packing pressure: 60 MPa), E2 (packing pressure Time: 1 s), F3 (cooling time: 20 s). According to the amount of change in the graph, it is judged that the control factor A has the greatest influence on this quality characteristic, followed by D, C, B, E, and F.
(2)
ANOVA
From ANOVA, the larger the F value is, the greater the contribution is, and it is expressed as a significant factor. Generally, the F value less than 5 is regarded as a factor with a relatively low contribution and its error is incorporated into the combined error. The ANOVA of tensile strength as shown in Table 4. The most significant factor is A (glass fiber), followed by D (packing pressure), C (injection speed), B (melt temperature).
In order to effectively evaluate each observation value and calculate its confidence interval, the expected mean value of the calculation confirmation experiment is:
CI S / N = F α ; 1 , V 2 × MSE × 1 n eff + 1 r = 5.12 × 0.071785 × 1 18 1 + 8 + 1 5 = 0.5072
Its S N = 39.93545   db , 95% confidence interval is 39.4282 ≤ μconfirmation ≤ 40.4427.

4.2.2. Hardness Test Data Analysis

(1)
MEA
From the S/N ratio obtained from the experiment as shown in Table 4, the main effect of each control factor was calculated, and the response graph was drawn as shown in Figure 9. It shows that the optimal factor levels are A3 (glass fiber: 20%), B2 (melt temperature: 185 °C), C2 (injection speed: 60 mm/s), D2 (packing pressure: 60 MPa), E2 (packing time: 1 s), F2 (cool time: 15 s). According to the amount of change in the graph, it can be judged that the control factor A has the greatest influence on this quality characteristic, followed by F, C, E, B, and D.
(2)
ANOVA
It can be seen from Table 5 that the most significant factor is A (glass fiber), followed by F (cooling time), and C (injection speed). Since the F values of E, B, and D are less than 5, the contribution is considered relatively low factor into the combined error.
In order to effectively evaluate each observation value and calculate its confidence interval, the expected mean value of the calculation confirmation experiment is:
CI S / N = F α ; 1 , V 2 × MSE × 1 n eff + 1 r = 4.84 × 0.000837 × 1 18 1 + 6 + 1 5 = 0.0488
Its S N = 39.93545   db , 95% confidence interval is 38.7035 ≤ μconfirmation ≤ 38.8011.

4.2.3. Impact Strength Test Data Analysis

(1)
MEA
From the S/N ratio obtained from the experiment as shown in Table 3, the main effect of each control factor is calculated, and the response graph is drawn, as shown in Figure 10. It shows that the best factor levels are A3 (glass fiber 20%), B2 (melt temperature 185 °C), C3 (injection speed 80 mm/s), D2 (packing pressure 60 MPa), E2 (packing time 1 s), F3 (cooling time 20 s). According to the variation of the graph, it can be observed that the control factor E has the greatest influence on this quality characteristic, followed by A, D, C, B, and F.
(2)
ANOVA
From ANOVA Table 6, it shows that the most significant factor is E (packing time), followed by A (glass fiber), D (packing pressure), C (injection speed), and B (melt temperature).
In order to effectively evaluate each observation value and calculate its confidence interval, the expected mean value of the calculation confirmation experiment is:
CI S / N = F α ; 1 , V 2 × MSE × 1 n eff + 1 r = 5.59 × 0.269678 × 1 18 1 + 10 + 1 5 = 1.1058
Its S N = 12.13171   db , 95% confidence interval is 11.0259 ≤ μconfirmation ≤ 13.2375.

4.2.4. Bending Strength Experiment Data Analysis

(1)
MEA
From the S/N ratio obtained from the experiment as shown in Table 3, the main effect of each control factor was calculated, and the response graph was drawn, as shown in Figure 11. It shows that the optimal factor level selection is A3 (glass fiber is 20%), B3 (melt temperature 195 °C), C2 (injection speed 60 mm/s), D3 (holding pressure 70 MPa), E3 (holding time 1.5 s), F2 (cooling time 15 s). According to the variation of the graph, it can be observed that the control factor A has the greatest influence on this quality characteristic, followed by C, D, F, B, and E.
(2)
ANOVA
From ANOVA Table 7, it shows that the most significant factor is A (glass fiber), followed by C (injection speed), D (packing pressure), and F (cooling time).
In order to effectively evaluate each observation value and calculate its confidence interval, the expected mean value of the calculation confirmation experiment is:
CI S / N = F α ; 1 , V 2 × MSE × 1 n eff + 1 r = 5.12 × 0.344485 × 1 18 1 + 8 + 1 5 = 1.1111
Its S N = 41.82273   db , 95% confidence interval is 40.7116 ≤ μconfirmation ≤ 42.9339.

4.3. Multiple-Quality Optimization Analysis

In this section, the Taguchi method is used in conjunction with PCA and DEA to obtain multiple quality optimization process parameters.

4.3.1. PCA

Step 1. From Table 4, normalize the S/N ratio data of each quality according to Equation (20), as shown in Table 8.
Step 2. Calculate the correlation coefficient matrix of the normalized data according to Equation (21), as shown in Table 9.
Step 3. Use the correlation coefficient matrix to calculate the eigenvalues and the eigenvectors, such as in Table 10 and Table 11; According to Equation (22), the variation of each principal component in the total variation is obtained.
Step 4. Combine the normalized data in Table 8 and the eigenvectors in Table 11, and calculate the total scores of the principal components according to Equation (23), as shown in Table 12.
Step 5. Multi-quality optimal parameter combination. The principal component total scores corresponding to the various control factors are shown in Table 13.
The best combination of parameters is A2 (glass fiber: 15%), B2 (melt temperature: 185 °C), C1 (injection speed: 40 mm/s), D2 (packing pressure: 60 MPa), E1 (packing time: 0.5 s), F3 (cooling time: 20 s).

4.3.2. DEA

Step 1. According to Equations (26) and (27), the relative efficiency is calculated from Table 4, as shown in Table 14 and the optimal weight of output and input is shown in Table 15.
Step 2. According to Equation (28), Table 14 and Table 15 are sorted by cross efficiency, and calculate the level value of the corresponding control factor in the orthogonal table, as shown in Table 16.
Table 16 shows that the best parameter combinations are A3 (glass fiber: 20%), B2 (melt temperature: 185 °C), C3 (injection speed: 80 mm/s), D2 (packing pressure: 60 MPa), E2 (packing time: 1 s), F2 (cooling time: 15 s).

5. Discussions

5.1. S/N Ratio Additive Model

Use the S/N ratio addition model to predict the S/N ratio of the best combination to verify the rationality of the confirmation experimental data.
(1)
PCA
The best combined S/N ratio addition model of PCA is shown in Table 17. For example, the S/N ratio addition model of the strength quality of the optimal factor level combination is calculated as follows:
ρ A 2 B 2 C 1 D 2 E 1 F 3 = 38.44403 + ( 38.4606 38.44403 ) + ( 38.68471 38.4403 ) + ( 38.71073 38.44403 ) + ( 38.69628 38.44403 ) + ( 38.25426 38.44403 ) + ( 38.51775 38.44403 ) = 39.10419
Similarity, the best combined S/N ratio addition model of DEA is shown in Table 18.
(2)
DEA
Use the S/N ratio addition model to predict the S/N ratio of the best combination to verify the rationality of the confirmation experimental data.
ρ A 3 B 2 C 3 D 2 E 2 F 2 = 38.44403 + ( 39.17581 38.44403 ) + ( 38.68471 38.44403 ) + ( 38.4098 38.44403 ) + ( 38.69628 38.44403 ) + ( 38.55264 38.44403 ) + ( 38.47925 38.44403 ) = 39.77835
Similarity, the prediction of all qualities is shown in Table 18.

5.2. S/N Ratio Additive Model Comparison

From Table 19, it can be seen that the S/N ratio of the optimal factor level combination of DEA in total qualities can be improved by 5.537101 db compared with the PCA expectation, so it can be predicted that the optimal factor level combination of multiple qualities is A3B2C3D2E2F2.

5.3. Confirmation Experiment and Comparison

The best processing parameters are actually processed the test pieces on the injection molding machine, and carry out the confirmation experiment. Each group of experiments is performed 5 times, as shown in Table 20 and Table 21, and the comparison is as follows.
(1)
The S/N ratio of the confirmation experiment of the two methods are similar to those predicted by the S/N ratio additive model.
(2)
The average confirmation experiment data of DEA: tensile strength 95.03775 MPa, hardness 86.52 Shore D, impact strength 4.4408 J/cm2, bending strength 119.889 MPa.
(3)
The average confirmation experiment data of PCA: tensile strength 94.03601 MPa, hardness 86.28 Shore D, impact strength 3.285046 J/cm2, bending strength 98.21989 MPa.
(4)
The Taguchi method combined with DEA, the obtained optimal combination of process parameters has the characteristics of better and multi-quality considerations.
The comparison of the multiple quality confirmation experiment group with single quality best experiment group from Taguchi experiment is shown in Table 22. It is observed that the optimal combination of process parameters obtained from DEA can meet the goal of the best multi-quality optimization.

6. Conclusions

In this paper, polylactide with glass fiber composites were synthesized via injection molding process and optimized with process parameters. First, the Taguchi orthogonal table is used to conduct experiments, and the optimal parameters of the single-quality process are obtained through MEA and ANOVA. Then, the PCA and DEA was combined to get the optimal process parameters for multiple qualities, and five confirmation experiments are carried out respectively to verify the ability of multi-quality consideration. The optimal process conditions are found to be glass fiber addition of 20%, melt temperature of 185 °C, injection speed of 80 mm/s, holding pressure of 60 MPa, retaining time of 1 s, and cooling time of 15 s. The associated mechanical properties are tensile strength of 95.04 MPa, hardness of 86.52 Shore D, impact strength of 4.4408 J/cm2, and bending strength of 4.4408 J/cm2. This research successfully boosts several properties of the PLA/GF composite. The composite material used in this study, the degradability of polylactic acid and the recyclability of glass fiber can reduce environmental pollution, and the mechanical properties can also be enhanced at the same time, that non-decomposable plastic materials cannot achieve.

Author Contributions

C.-H.H.: Conceptualization, Methodology, Resources, Writing. C.-C.H.: Supervision, Writing–Editing, Visualization. C.-F.J.K.: Conceptualization, Methodology, Writing—Original Draft. N.A.: Methodology, Editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology of Taiwan under Grant No. 110-2622-E-011-012. And the APC was funded by bcc830e18ae179cd.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request.

Acknowledgments

This research was supported by the Ministry of Science and Technology of the Republic of China under Grant No. 110-2622-E-011-012.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Isaac, C.W. Crashworthiness performance of green composite energy absorbing structure with embedded sensing device providing cleaner environment for sustainable maintenance. Sustain. Mater. Technol. 2020, 25, 00196. [Google Scholar] [CrossRef]
  2. Çavdar, A.D.; Boran, S. Dogal liflerin otomotiv sanayinde kullanımı. Kast. Univ. J. For. Fac. 2016, 16, 253–263. [Google Scholar]
  3. Zhao, X.; Liu, J.; Li, J.; Liang, X.; Zhou, W.; Peng, S. Strategies and techniques for improving heat resistance and mechanical performances of poly(lactic acid) (PLA) biodegradable materials. Int. J. Biol. Macromol. 2022, 218, 115–134. [Google Scholar] [CrossRef] [PubMed]
  4. Ashothamana, A.; Sudhab, J.; Senthilkumar, N. A comprehensive review on biodegradable polylactic acid polymer matrix composite material reinforced with synthetic and natural fiber. Mater. Today Proc. 2023, 80, 2829–2839. [Google Scholar] [CrossRef]
  5. Nagarajan, V.; Mohanty, A.K.; Misra, M. Perspective on polylactic acid (PLA) based sustainable materials for durable applications: Focus on toughness and heat resistance. ACS Sustain. Chem. Eng. 2016, 4, 2899–2916. [Google Scholar] [CrossRef]
  6. Jandas, P.J.; Mohanty, S.; Nayak, S.K. Surface treated banana fiber reinforced poly (lactic acid) nanocomposites for disposable applications. J. Clean. Prod. 2013, 52, 392–401. [Google Scholar] [CrossRef]
  7. Joshi, S.V.; Drzal, L.T.; Mohanty, A.K.; Arora, S. Are natural fiber composites environmentally superior to glass fiber reinforced composites? Compos. Part A Appl. Sci. Manuf. 2004, 35, 371–376. [Google Scholar] [CrossRef]
  8. Jang, J.; Im, H.G.; Lim, D.; Bae, B.S. Preparation of high-performance transparent glass-fiber reinforced composites based on refractive index-tunable epoxy-functionalized siloxane hybrid matrix. Compos. Sci. Technol. 2021, 201, 108527. [Google Scholar] [CrossRef]
  9. Sathishkumar, T.P.; Satheeshkumar, S.; Naveen, J. Glass fiber-reinforced polymer composites—A review. J. Reinf. Plast. Compos. 2014, 33, 1258–1275. [Google Scholar] [CrossRef]
  10. Al-Homoud, M.S. Performance characteristics and practical applications of common building thermal insulation materials. Build. Environ. 2005, 40, 353–366. [Google Scholar] [CrossRef]
  11. Wang, G.; Zhang, D.; Li, B.; Wan, G.; Zhao, G.; Zhang, A. Strong and thermal-resistance glass fiber-reinforced polylactic acid (PLA) composites enabled by heat treatment. Int. J. Biol. Macromol. 2019, 129, 448–459. [Google Scholar] [CrossRef] [PubMed]
  12. Sun, Y.; Zheng, Z.; Wang, Y.; Yang, B.; Wang, J.; Mu, W. PLA composites reinforced with rice residues or glass fiber—A review of mechanical properties, thermal properties, and biodegradation properties. J. Polym. Res. 2022, 29, 422. [Google Scholar] [CrossRef]
  13. Wang, G.; Zhang, D.; Wan, G.; Li, B.; Zhao, G. Glass fiber reinforced PLA composite with enhanced mechanical properties, thermal behavior, and foaming ability. Polymer 2019, 181, 121803. [Google Scholar] [CrossRef]
  14. Leu, Y.Y.; Mohd, Z.A.; Chow, W.S. Mechanical thermal and morphological properties of injection molded polylactic acid/sebs-g-mah/organo-montmorillonite nanocomposites. J. Appl. Polym. Sci. 2012, 124, 1200–1207. [Google Scholar] [CrossRef]
  15. Ma, P.; Wang, X.; Liu, B.; Li, Y.; Chen, S.; Zhang, Y.; Xu, G. Preparation and foaming extrusion behavior of polylactide/polybutylene ccinate/montmorillonoid nanocomposite. J. Cell. Plast. 2012, 48, 191–205. [Google Scholar] [CrossRef]
  16. Shen, L.H.; Yang, H.J.; Ying, J.; Qiao, F.; Peng, M. Preparation and mechanical properties of carbon fiber reinforced hydroxyapatite/polylactide biocomposites. J. Mater. Sci. Mater. Med. 2009, 20, 2259–2263. [Google Scholar] [CrossRef]
  17. Bledzki, A.K.; Jaszkiewicz, A.; Scherzer, D. Mechanical properties of pla composites with man-made cellulose and abaca fibres. Compos. Part A Appl. 2009, 40, 404–412. [Google Scholar] [CrossRef]
  18. An, S.O.; Lee, E.S.; Noh, S.L. A study on the cutting characteristics of glass fiber reinforced plastics with respect to tool materials and geometries. J. Mater. Process. Technol. 1997, 68, 60–67. [Google Scholar] [CrossRef]
  19. Bigg, D.M. The impact behavior of thermoplastic sheet composites. J. Reinf. Plast. Compos. 1994, 13, 339–354. [Google Scholar] [CrossRef]
  20. Jaszkiewicz, A.A.; Bledzki, A.K.; Franciszczak, P. Improving the mechanical performance of pla composites with natural, man-made cellulose and glass fibers—A comparison to pp counterparts. Polimery 2013, 58, 435–442. [Google Scholar] [CrossRef]
  21. Durante, M.; Formisano, A.; Boccarusso, L.; Langella, A.; Carrino, L. Creep behaviour of polylactic acid reinforced by woven hemp fabric. Compos. Part B 2017, 124, 16–22. [Google Scholar] [CrossRef]
  22. Li, M.X.; Kim, S.H.; Choi, S.W.; Goda, K.; Lee, W.I. Effect of reinforcing particles on hydrolytic degradation behavior of poly (lactic acid) composites. Compos. Part B 2016, 96, 248–254. [Google Scholar] [CrossRef]
  23. Luo, J.; Meng, X.; Gong, W.; Jiang, Z.; Xin, Z. Improving the stability and ductility of polylactic acid via phosphite functional polysilsesquioxane. RSC Adv. 2019, 9, 25151–25157. [Google Scholar] [CrossRef] [PubMed]
  24. Meinig, L.; Boldt, R.; Spoerer, Y.; Kuehnert, I.; Stommel, M. Correlation between processing parameters, morphology, and properties of injection-molded polylactid acid (PLA) specimens at different length scales. Polymers 2023, 15, 721. [Google Scholar] [CrossRef]
  25. Setyarini, P.H.; Anggarwati, M.P.; Wahyudi, S.; Sulistyarini, D.H. Injection Molding Characterization of PLA and Chitosan Mixtures for Biomaterial Applications. Adv. Sci. Technol. 2023, 122, 27–32. [Google Scholar]
  26. Mukras, S.M.; Omar, H.M.; Al-Mufadi, F.A. Experimental-based multi-objective optimization of injection molding process parameters. Arab. J. Sci. Eng. 2019, 44, 7653–7665. [Google Scholar] [CrossRef]
  27. Johnnaber, F. Injection Molding Machine: A Use’s Guide; Hanser Publication: Munich, Germany, 2008. [Google Scholar]
  28. Bozzelli, J.W.; Cardinal, J.; Fierens, B. Pressure Loss in Thin Wall Moldings; William Andrew Publishing: Norwich, NY, USA, 2001; pp. 121–126. [Google Scholar]
  29. Jansen, K.M.B.; Van Dijk, D.J.; Husselman, M.H. Effect of processing conditions on shrinkage in injection molding. Polym. Eng. Sci. 1998, 38, 838–846. [Google Scholar] [CrossRef] [Green Version]
  30. Shokri, P.; Bhatnagar, N. Effect of the post-filling stage on fiber orientation at the mid-plane in injection molding of reinforced thermoplastics. Phys. Procedia 2012, 25, 79–85. [Google Scholar] [CrossRef] [Green Version]
  31. Kamaruddin, S.; Khan, Z.A.; Foong, S.H. Application of Taguchi method in the optimization of injection Moulding parameters for manufacturing products from plastic blend. Int. J. Eng. Technol. 2010, 2, 574–580. [Google Scholar] [CrossRef] [Green Version]
  32. Kuo, C.F.J.; Su, T.L. Optimization of injection molding processing parameters for LCD light-guide plates. J. Mater. Eng. Perform. 2007, 16, 539–548. [Google Scholar] [CrossRef]
  33. Karna, S.K.; Singh, R.V.; Sahai, R. Application of Taguchi Method in Indian Industry. Int. J. Emerg. Technol. Adv. Eng. 2012, 2, 387–391. [Google Scholar]
  34. Liu, S.J.; Hsu, C.H.; Chang, C.Y. Parametric characterization of the thin-wall injection molding of thermoplastic composites. J. Reinf. Plast. Compos. 2002, 21, 1027–1041. [Google Scholar] [CrossRef]
  35. Ghani, J.A.; Choudhury, I.A.; Hassan, H.H. Application of Taguchi method in the optimization of end milling parameters. J. Mater. Process. Technol. 2004, 145, 84–92. [Google Scholar] [CrossRef]
  36. Su, C.T.; Tong, L.I. Multi-response robust design by principal component analysis. Total Qual. Manag. 1997, 8, 409–416. [Google Scholar] [CrossRef]
  37. Antony, J. Multi-response optimization in industrial experiments using Taguchi’s quality loss function and principal component analysis. Qual. Reliab. Eng. Int. 2000, 16, 3–8. [Google Scholar] [CrossRef]
  38. Shih, J.-S.; Tzeng, Y.-F.; Yang, J.-B. Principal component analysis for multiple quality characteristics optimization of metal inert gas welding aluminum foam plate. Mater. Des. 2011, 32, 1253–1261. [Google Scholar] [CrossRef]
  39. Jeyapaul, R.; Shahabudeen, P.; Krishnaiah, K. Simultaneous optimization of multi-response problems in the Taguchi method using genetic algorithm. Int. J. Adv. Manuf. Technol. 2006, 30, 9–10. [Google Scholar] [CrossRef]
  40. Al-Refaie, A.; Al-Tahat, M.D. Solving the multi-response problem in Taguchi method by benevolent formulation in DEA. J. Intell. Manuf. 2011, 22, 505–521. [Google Scholar] [CrossRef]
  41. Taguchi, G. Introduction to Quality Engineering: Designing Quality into Products and Processes; Asian Productivity Organization: Tokyo, Japan, 1986. [Google Scholar]
  42. Silkman, R.H. Measuring Efficiency: An Assessment of Data Envelopment Analysis; Jossey-Bass: Hoboken, NJ, USA, 1986. [Google Scholar]
  43. Kumar, A.J.; Prakash, M. Thermal properties of basalt/Cissus-quadrangularis hybrid fiber reinforced polylactic acid biomedical composites. J. Therm. Anal. Calorim. 2020, 141, 717–725. [Google Scholar] [CrossRef]
Figure 1. The operation principle of the injection molding machine.
Figure 1. The operation principle of the injection molding machine.
Polymers 15 03018 g001
Figure 2. The injection samples.
Figure 2. The injection samples.
Polymers 15 03018 g002
Figure 3. Possible reaction between polylactic acid (PLA) and glass fiber.
Figure 3. Possible reaction between polylactic acid (PLA) and glass fiber.
Polymers 15 03018 g003
Figure 4. Scheme of experiment and processing.
Figure 4. Scheme of experiment and processing.
Polymers 15 03018 g004
Figure 5. PLA/10%GF (type: GG-0010N) DSC diagram.
Figure 5. PLA/10%GF (type: GG-0010N) DSC diagram.
Polymers 15 03018 g005
Figure 6. PLA/15%GF (type: GG-0015N) DSC diagram.
Figure 6. PLA/15%GF (type: GG-0015N) DSC diagram.
Polymers 15 03018 g006
Figure 7. PLA/20%GF type: GG-0020N) DSC diagram.
Figure 7. PLA/20%GF type: GG-0020N) DSC diagram.
Polymers 15 03018 g007
Figure 8. Response graph of Tensile strength.
Figure 8. Response graph of Tensile strength.
Polymers 15 03018 g008
Figure 9. Hardness factor response graph.
Figure 9. Hardness factor response graph.
Polymers 15 03018 g009
Figure 10. Impact strength response graph.
Figure 10. Impact strength response graph.
Polymers 15 03018 g010
Figure 11. Bending factor response graph.
Figure 11. Bending factor response graph.
Polymers 15 03018 g011
Table 1. PLA/GF material properties.
Table 1. PLA/GF material properties.
TypeGG-0010NGG-0015NGG-0020N
Category TY11512706TY11512707TY11512708
Raw material propertiesRatioRatioRatio
Filling contents (%)101520
Mold shrinkage (%)0.080.070.055
Melting point (°C)155155155
Specific weight1.3021.3521.373
Table 2. Injection molding machine control factors and their level values.
Table 2. Injection molding machine control factors and their level values.
FactorABCDEF
Level GF (%)Melt Temperature
(°C)
Injection Speed (mm/s)Packing Pressure (MPa)Packing Time (s)Cololing Time (s)
11017540500.510
2151856060115
32019580701.520
Table 3. L18 orthogonal array of experimental data.
Table 3. L18 orthogonal array of experimental data.
Exp.
No.
Tensile Strength Shore HardnessImpact StrengthBending Strength
Mean
(Mpa)
S/N Ratio
(db)
Mean
(Shore D)
S/N Ratio
(db)
Mean
(Mpa)
S/N Ratio
(db)
Mean (Mpa)S/N Ratio
(db)
173.5237.3283.5238.432.849.0566.8136.49
279.3837.9985.3638.623.8211.6473.5737.33
374.3537.4284.1638.503.3010.3669.9436.89
489.8139.0685.3238.623.3110.4181.8538.26
579.9738.0585.0838.592.728.65104.7840.40
681.2738.1985.838.662.678.4794.1039.46
790.0139.0886.6438.753.009.51114.2941.15
894.5839.5185.5238.643.4310.6489.0238.98
990.9339.1785.2838.613.2210.09113.0441.06
1074.1837.4084.9638.583.4110.5676.2337.63
1184.2038.5083.638.443.6811.2965.0236.25
1275.2237.5284.238.502.658.4566.9536.51
1376.8437.7185.1638.602.427.5993.9439.45
1486.9738.7885.4438.633.139.9085.1038.59
1588.5338.9485.1238.603.089.77102.6840.22
1690.5939.1385.5638.643.5911.0694.1039.46
1791.7739.2586.2438.713.139.90105.7440.47
1888.0438.8985.8838.673.4310.67117.2741.38
Table 4. ANOVA of tensile strength.
Table 4. ANOVA of tensile strength.
Source of VarianceDOFSSMSF-RatioSS′Contribution
(%)
A26.5749423.28747178.86536.49157369.43184
B20.5360760.2680386.4301370.4527064.842008
C20.7580980.3790499.0932560.6747297.216688
D20.8343860.41719310.008310.7510168.032637
E20.3263810.1631913.914890.2430122.599183
F20.1112560.0556281.3344990.0278870.29827
Error50.2084230.041685-0.7086397.579378
Combined
error
90.6460610.071785-0.97953810.47683
Total179.349562--9.349562100%
Table 5. Hardness ANOVA table.
Table 5. Hardness ANOVA table.
Source of Variance DOFSSMSF-RatioSSContribution (%)
A20.0783810.0391950.177950.07681961.88925
B20.0006450.0003220.41269−0.00092−0.73912
C20.0098790.0049396.3243270.0083176.700536
D20.00005830.00002920.037343−0.0015−1.21148
E20.0046010.00232.9454480.0030392.448299
F20.0266540.01332717.063460.02509220.21547
Error50.0039050.000781-0.01327710.69704
Combined error110.0092090.000837-0.01389511.19475
Total170.124123--0.124123100%
Table 6. ANOVA of impact strength.
Table 6. ANOVA of impact strength.
Source of VariationDOFSSMSF-RatioSSContribution (%)
A25.1948632.59743120.265454.93852224.29646
B21.8068170.9034087.0484921.5504767.628006
C21.9880540.9940277.7555111.7317138.519656
D22.0133871.0066937.8543341.7570468.644286
E27.4352313.71761629.005257.1788935.31859
F21.2468920.6234464.8641950.9905514.873297
error50.6408520.12817-2.17889710.7197
combined error71.8877440.269678-3.16944815.593
Total1720.3261--20.3261100%
Table 7. Bending ANOVA table.
Table 7. Bending ANOVA table.
Source of VariationDOFSSMSF-RatioSS′Contribution
(%)
A240.4993920.2497109.810240.1305879.54693
B21.2178930.6089463.3021980.849081.683048
C22.3331721.1665866.3261721.964363.893759
D22.3072241.1536126.2558151.9384113.842324
E20.9604430.4802222.6041490.5916311.172732
F22.2087771.1043895.9888871.8399653.647182
error50.9220320.184406-3.1349086.214022
combined error93.1003680.344485-4.5756199.069802
Toatl1750.44894--50.44894100%
Table 8. Normalization of quality data.
Table 8. Normalization of quality data.
ItemNormalization
Exp. No.Tensile Strength
(db)
Hardness
(db)
Impact Strength
(db)
Bending Strength
(db)
1000.3602880.046529
20.3055140.59510110.209861
30.0438140.2116060.683120.124504
40.7966520.5833440.6951650.39112
50.3356280.5063420.2608250.80949
60.3989770.7349690.2177050.626838
70.80587710.4746390.955724
810.6464520.7518090.532859
90.8454590.5704360.6162770.937851
100.0361970.468590.7339540.269505
110.539950.0298640.9125180
120.0928010.2233510.2126930.050202
130.1772340.53142500.624045
140.669610.6210670.56990.456059
150.7400980.5196410.5386670.77504
160.8301430.6568110.8554350.625622
170.882720.8743930.5700410.823749
180.7174390.7603870.7595011
Table 9. Correlation coefficient matrix.
Table 9. Correlation coefficient matrix.
Correlation CoefficientTensile StrengthHardnessImpact StrengthBending Strength
tensile strength10.6377120.3641630.631945
hardness0.63771210.0250610.804754
impact strength0.3641630.0250611−0.15557
bending strength0.6319450.804754−0.155571
Table 10. Eigenvalues and variances.
Table 10. Eigenvalues and variances.
Principal ComponentEigenvaluesVariance (%)Variance Accumulation (%)
12.397259.931559.9315
21.172929.3232389.25473
30.27646.91017396.1649
40.15343.835096100
Table 11. Eigenvectors.
Table 11. Eigenvectors.
Principal Component EigenvalueEigenvector
PC1PC2PC3PC4
tensile strength0.33270.6983−0.30610.555
hardness0.4779−0.63810.13170.5891
impact strength−0.2964−0.3228−0.89430.0906
bending strength−0.75710.03040.29870.5803
Table 12. The principal component scores.
Table 12. The principal component scores.
PC No.PC1PC2PC3PC4MPCI
1−0.14202−0.11489−0.308310.059643−0.13782
2−0.06924−0.48281−0.846760.732517−0.21349
3−0.18104−0.32116−0.559270.283114−0.23046
40.04166−0.02844−0.671891.0757380.011456
5−0.33653−0.14831−0.027510.977937−0.20957
6−0.05513−0.2416−0.032791.037881−0.06634
7−0.11825−0.19952−0.253971.633971−0.08426
80.015380.059314−0.734141.3131570.026238
9−0.338820.055965−0.454671.405343−0.16417
10−0.18560−0.50246−0.525240.519025−0.27496
11−0.076560.06343−0.977410.399939−0.07948
120.03657−0.14485−0.174210.231483−0.02372
13−0.15953−0.196370.202140.77356−0.10956
140.00539−0.09881−0.496611.053788−0.01965
15−0.251880.034907−0.408331.215434−0.12232
16−0.13713−0.09654−0.745751.288207−0.11262
17−0.08107−0.10051−0.418781.534682−0.04814
18−0.38013−0.19898−0.499991.495233−0.26337
Notes: PC: principal component. MPCI: multiple performance characteristic index.
Table 13. Total scores of the principal component.
Table 13. Total scores of the principal component.
FactorABCDEF
Level
1−0.15999−0.11796−0.09008−0.10849−0.05989−0.11126
2−0.08600−0.09068−0.15066−0.07241−0.15403−0.13492
3−0.10772−0.14506−0.11297−0.17281−0.13979−0.10753
Optimal combinationA2B2C1D2E1F3
Table 14. The relative efficiency of each DMUj.
Table 14. The relative efficiency of each DMUj.
DMUjInputOutputCCR
Relative Efficiency
x1jTensile Strength
(y1j)
Hardness (y2j)Impact Strength
(y3j)
Bending Strength
(y4j)
EO
DMU1173.5247783.522.847666.815450.963989
DMU2179.3810885.363.82867773.575611
DMU3174.3587784.163.301269.940870.977845
DMU4189.8110885.323.31756981.859090.991989
DMU5179.9776985.082.720092104.7840.981994
DMU6181.2717585.82.67270894.103770.990305
DMU7190.0175386.643.008923114.29641
DMU8194.5827885.523.43344689.029341
DMU9190.9370885.283.220646113.04141
DMU10174.1863584.963.41707776.236440.988585
DMU11184.2073383.63.68006265.027890.991002
DMU12175.2298184.22.65612366.955240.971837
DMU13176.8438285.162.42113893.948010.982918
DMU14186.9790385.443.13836985.100430.989086
DMU15188.5363785.123.089477102.68130.985056
DMU16190.5979585.563.59756994.101141
DMU17191.7759986.243.138585105.74531
DMU18188.0410585.883.4328117.27631
Table 15. The optimal weight of each DMUj.
Table 15. The optimal weight of each DMUj.
DMUjInputOutput
v*1ju*1ju*2ju*3ju*4j
DMU1100.01154200
DMU210.00004620.0107760.0192730.0000372
DMU3100.0109480.0170940
DMU410.0008280.0100680.0176540
DMU5100.01154200
DMU6100.01154200
DMU710.0033150.0063940.0218960.000715
DMU810.0062340.0017920.040550.001324
DMU910.0057720.0020660.0467610.001311
DMU10100.0109480.0170940
DMU1110.00377100.183010
DMU12100.01154200
DMU13100.01154200
DMU1410.0002080.0106540.0193390
DMU1510.0002080.0106540.0193390
DMU1610.0011580.0094060.0212170.000149
DMU1710.0036160.0059630.023340.000762
DMU1810.0008410.0098090.0207680.000105
Table 16. DEA Cross-efficiency sorting corresponds to the control factor level.
Table 16. DEA Cross-efficiency sorting corresponds to the control factor level.
FactorABCDEF
Level
15.39.09.89.27.58.0
27.710.88.511.711.210.5
315.58.710.27.79.810.0
Optimal combinationA3B2C3D2E2F2
Table 17. The best combined S/N ratio addition model of PCA.
Table 17. The best combined S/N ratio addition model of PCA.
Best CombinationTensile StrengthHardnessImpact StrengthBending Strength
A238.460638.620839.13508539.40227
B238.6847138.608910.3398938.67565
C138.7107338.5718110.0884138.79719
D238.6962838.6063110.3322438.46438
E138.2542638.594899.02038738.56648
F338.5177538.5816910.1621738.53941
ρ A 2 B 2 C 1 D 2 E 1 F 3 39.1041938.565349.61189938.98086
Table 18. The best combined S/N ratio addition model of DEA.
Table 18. The best combined S/N ratio addition model of DEA.
Optimal
Combination
Tensile StrengthHardnessImpact StrengthBending Strength
A339.1758138.6747810.3146240.42158
B238.6847138.608910.3398938.67565
C338.409838.6124410.168838.50768
D238.6962838.6063110.3322438.46438
E238.5526438.6262710.5491839.04502
F238.4792538.657949.98166739.38358
ρ A 3 B 2 C 3 D 2 E 2 F 2 39.7783538.7675512.2201241.03337
Table 19. DEA addition model improvement.
Table 19. DEA addition model improvement.
MethodDEA
S/N (db)
PCA
S/N (db)
Improvement
S/N (db)
Quality
tensile strength39.7783539.104190.67416
hardness38.7675538.565340.20221
impact strength12.220129.6118992.608221
bending strength41.0333738.980862.05251
Table 20. PCA’s confirmation experiment.
Table 20. PCA’s confirmation experiment.
Group12345AverageLTB
S/N (db)
Quaty
tensile strength94.1407493.8675393.9130794.1862794.0724394.0360139.46586
hardness86.486.286.286.286.486.2838.71819
impact strength3.2846153.2846153.4493853.2846153.1223.28504610.31786
bending strength97.2073198.2198998.2198998.2198999.2324798.2198939.84344
Table 21. PCA’s confirmation experiment.
Table 21. PCA’s confirmation experiment.
Group12345AverageLTB
S/N (db)
Quality
tensile strength94.5277795.0514195.3701495.5522894.6871495.0377539.5577
hardness86.286.886.6878686.5238.74209
impact strength4.3583084.5230774.4412314.5230774.3583084.440812.94564
bending strength119.484120.4966120.4966119.484119.484119.88941.57537
Table 22. The comparison of multiple quality confirmation experiment group with single quality best experiment group.
Table 22. The comparison of multiple quality confirmation experiment group with single quality best experiment group.
Quality Tensile Strength
(MPa)
Hardness
(Shore D)
Impact Strength (J/cm2)Bending Strength
(MPa)
Group
PCA confirmation experimental group94.0360186.283.28504698.21989
DEA confirmation experimental group95.0377586.524.4408119.889
Taguchi group 894.5827885.523.43344689.02934
Taguchi group 790.0175386.643.008923114.2964
Taguchi group 279.3810885.363.82867773.57561
Taguchi Group 1888.0410585.883.4328117.2763
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hsiao, C.-H.; Huang, C.-C.; Kuo, C.-F.J.; Ahmad, N. Integration of Multivariate Statistical Control Chart and Machine Learning to Identify the Abnormal Process Parameters for Polylactide with Glass Fiber Composites in Injection Molding; Part I: The Processing Parameter Optimization for Multiple Qualities of Polylactide/Glass Fiber Composites in Injection Molding. Polymers 2023, 15, 3018. https://doi.org/10.3390/polym15143018

AMA Style

Hsiao C-H, Huang C-C, Kuo C-FJ, Ahmad N. Integration of Multivariate Statistical Control Chart and Machine Learning to Identify the Abnormal Process Parameters for Polylactide with Glass Fiber Composites in Injection Molding; Part I: The Processing Parameter Optimization for Multiple Qualities of Polylactide/Glass Fiber Composites in Injection Molding. Polymers. 2023; 15(14):3018. https://doi.org/10.3390/polym15143018

Chicago/Turabian Style

Hsiao, Chi-Hao, Chang-Chiun Huang, Chung-Feng Jeffrey Kuo, and Naveed Ahmad. 2023. "Integration of Multivariate Statistical Control Chart and Machine Learning to Identify the Abnormal Process Parameters for Polylactide with Glass Fiber Composites in Injection Molding; Part I: The Processing Parameter Optimization for Multiple Qualities of Polylactide/Glass Fiber Composites in Injection Molding" Polymers 15, no. 14: 3018. https://doi.org/10.3390/polym15143018

APA Style

Hsiao, C. -H., Huang, C. -C., Kuo, C. -F. J., & Ahmad, N. (2023). Integration of Multivariate Statistical Control Chart and Machine Learning to Identify the Abnormal Process Parameters for Polylactide with Glass Fiber Composites in Injection Molding; Part I: The Processing Parameter Optimization for Multiple Qualities of Polylactide/Glass Fiber Composites in Injection Molding. Polymers, 15(14), 3018. https://doi.org/10.3390/polym15143018

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop