Hybrid Integration of Taguchi Parametric Design , Grey Relational Analysis , and Principal Component Analysis Optimization for Plastic Gear Production

The identification of optimal processing parameters is an important practice in the plastic injection moulding industry because of the significant effect of such parameters on plastic part quality and cost. However, the optimization design of injection moulding process parameters can be difficult because more than one quality characteristic is used in the evaluation.This study systematically develops a hybrid optimization method for multiple quality characteristics by integrating the Taguchi parameter design, grey relational analysis, and principal component analysis. A plastic gear is used to demonstrate the efficiency and validity of the proposed hybrid optimization method in controlling all influential injection moulding processing parameters during plastic gear manufacturing. To minimize the shrinkage behaviour in tooth thickness, addendum circle, and dedendum circle of moulded gear, the optimal combination of different process parameters is determined.The case study demonstrates that the proposed optimization method can produce plastic-moulded gear with minimum shrinkage behaviour of 1.8%, 1.53%, and 2.42% in tooth thickness, addendum circle, and dedendum circle, respectively; these values are less than the values in the main experiment. Therefore, shrinkage-related defects that lead to severe failure in plastic gears can be effectively minimized while satisfying the demand of the global plastic gear industry.


Introduction
Injection moulding is a complex process because of its requirements for numerous delicate adjustments.The quality of an injection-moulded part significantly depends on the selection of appropriate materials, parts, mould designs, and processing parameters.A review of previous works shows that the setting of processing parameters significantly affects the quality of plastic parts [1,2].The selection of injection moulding process parameters previously involves a trial-and-error method [3].However, obtaining an optimal parameter setting for complex-manufacturing processes is difficult because the trial-and-error method is a "one change at a time" test [4].This tuning exercise is repeated until the quality of the moulded part is found satisfactory, thus incurring high production costs and long setup times [5].Moreover, the adjustments and modifications of processing parameters rely heavily on the experience and intuition of the moulding personnel [6].Nevertheless, the growing demand in the industry for expert moulding personnel exceeds the supply; moreover, amateur moulding personnel require more than 10 years of experience to be considered experts [7].
In a fiercely competitive market, the trial-and-error approach cannot meet the challenges of globalization, particularly at the point where the disadvantages outweigh the advantages.Thus, the interest in substituting the trial-anderror method with a fast and reliable optimization approach such as the Taguchi method continues to increase.The Taguchi method considerably reduces the time and effort needed to obtain optimal process conditions and to determine the important factors that affect plastic part quality.The application of the Taguchi method has attracted considerable attention in the literature over the past 20 years and has been widely applied to various fields, such as manufacturing 2 Chinese Journal of Engineering systems [8], mechanical component design [9], and parameter optimization [10,11].Li et al. [12] exploited the Taguchi method to eliminate the weld line of the right door of a copy machine by optimizing the melt temperature, injection speed, and injection pressure.Wu and Liang [13] further improved the weld line strength of an injection-moulded part by optimizing the mould temperature, packing pressure, injection acceleration, and packing time.By focusing on warpage in the thin-walled plastic part, many researchers used the experimental design of the Taguchi method to optimize the injection moulding parameters; packing pressure was found to be the most influential factor on the warpage problem [14][15][16].Oktem et al. [17] improved the shrinkage and warpage in a plastic part by approximately 2.17% and 0.7%, respectively, by determining the optimal packing time, packing pressure, injection time, and cooling time by using the Taguchi method.
A number of studies have shown that the Taguchi method is a robust experimental design for obtaining the best combination of factors or levels with the lowest cost solution to achieve the quality requirements of plastic parts.Nevertheless, the Taguchi method exhibits a number of limitations.The optimization design of injection moulding process parameters can be difficult because more than one quality characteristic is used in the evaluation.Problems arise when the optimal process parameters contradict one another because of different mechanisms [18].To address this problem, many studies devised new experiment design methodologies that optimize multiple quality characteristics simultaneously while providing accurate results by integrating the Taguchi method with other techniques.One such technique is grey relational analysis (GRA).This technique was first proposed by Deng [19] to optimize the multiresponse problem by using the grey relational coefficient and grey relational grade.The grey relational coefficient can express the relationship between the desired and actual experiment results, whereas the grey relational grade is simultaneously computed to correspond to each quality characteristic.
Few works have studied the optimization of injection moulding process parameters by using the integration of Taguchi method and GRA.Fung [20] adopted the integration method to improve the wear volume losses in both the parallel and perpendicular directions and inspected the extent in which process parameters influence each quality characteristic by using the comparability sequence with a large value of grey relational grade.Another two studies conducted by Fung et al. [21] and Yang [22] used the same procedure to optimize the processing parameters for multiple quality characteristics.The former focused on the concurrent improvement of yield stress and elongation of the PC/ABS blend, whereas the latter examined the effects of processing parameters on the mechanical and tribological properties of PC composites, including ultimate stress, surface roughness, and friction coefficient.
Another technique that can be applied with the Taguchi method in addressing multiple quality characteristic problems is principal component analysis (PCA).By using PCA, a set of original responses is transformed into a set of uncorrelated components to find the optimal factor or level combination.PCA does nothing when responses are uncorrelated; thus, the best results are obtained when responses or quality characteristics are highly correlated (either positively or negatively) [23].The application of PCA involves a series of steps that are capable of solving the weakness of the standalone Taguchi method, which requires engineering judgment to handle multiple quality characteristics because the judgment of an engineer increases uncertainty during the decision-making process [24].Fung and Kang [25] applied PCA to the Taguchi method to improve friction properties, including friction coefficients and surface roughness, in different sliding directions (P-type and AP-type).The authors extracted the principal components and the coefficient of determination to establish a comprehensive index and obtain the final optimal parameter setting.
Although both techniques are widely integrated with the Taguchi method to overcome multiple quality characteristic problems, both integration techniques have several limitations.In combining the Taguchi method and GRA, a problem arises when calculating the value of the grey relational grade in GRA.An engineering judgment or subjective estimation is required to determine the weighting values for each quality characteristic, thus resulting in an increase in uncertainty during the decision-making process.This approach cannot objectively reveal the relative importance of multiple quality characteristic performance.For the integration of the Taguchi method and PCA, the primary problem is the difficulty of interpretation because the original variables are substituted and the principal components are heavily affected by the scaling of the variables.In both methods, no scientific analysis has been performed in terms of performance under optimal process conditions to verify the accuracy of the optimization results.Therefore, the effectiveness and reliability of the integration of the Taguchi method with GRA or PCA alone remain uncertain.Further validation experiments on the optimization results are needed to prove the effectiveness of the integrated approach.
This study develops a hybrid optimization method by integrating the robust parameter design of Taguchi, GRA, and PCA to solve problems with multiple quality characteristics.An injection-moulded plastic gear is used to demonstrate the efficiency and validity of the proposed hybrid optimization method in controlling all influential injection moulding processing parameters during plastic gear manufacturing.By focusing on minimizing the shrinkage behaviour in the tooth thickness, addendum circle, and dedendum circle of moulded gear, the optimal combination of four process parameters, namely, melt temperature, packing pressure, packing time, and cooling time, at three different levels is determined by using the proposed hybrid optimization method.

Taguchi/GRA/PCA Hybrid Method for the Optimization of Injection Moulding Process Parameters
2.1.Taguchi Method.The Taguchi method is used to design experiments based on the orthogonal arrays (OAs).Taguchi OAs are highly fractional orthogonal designs that can be used to study the whole parameter space with a small number of experiments.In analyzing the results, the Taguchi method uses a statistical measure of performance known as signal-tonoise (S/N) ratio.The S/N ratio is a measure of performance to develop products or processes that are insensitive to noise factors in a controlled manner [26].Noise factors are uncontrollable factors that influence product or process uncertainty.These factors include humidity and weather.Depending on the objective, three different methods can be used to calculate the S/N ratio in the Taguchi method: (1) smaller-the-better quality characteristic: (2) bigger-the-better quality characteristic: (3) nominal-the-better quality characteristic: where   is the th test at the th trial,  is the total number of tests, and  is the standard deviation.

GRA.
In the analysis of the processing parameters of the plastic gear in the case study, an appropriate mathematic model is established to study the relationship between target values and the quality characteristics of the plastic gear obtained from the experiment.The primary concern is to analyze the differences among the quality characteristics of plastic gear because of various processing parameters and to understand the relationship between quality characteristics and target values.GRA is a method that measures the correlation degree among factors based on the similarity or difference among factors.GRA is characterized by small data requirements and multifactor analysis.The procedure of the GRA is presented below.Grey relational generation involves data preprocessing and calculation according to the quality characteristics.The computing method of the grey relational generation is as follows.
(1) The-larger-the-better (the higher the target value, the better): (2) The-smaller-the-better (the smaller the target value, the better): (3) The-nominal-the-better characteristic (if the target is a specific value, set the target value as OB): ()  () is the measurement of the quality characteristic, max  ()  () is the largest value of  ()  (), and min  ()  () is the smallest value of  ()   ().We determine the difference sequence Δ 0 () and the minimum value Δ min and maximum value Δ max in the difference sequence, including the parameter values required to set the reference sequence: We set the identification coefficient , which generally has a value of 0.5.After the data preprocessing, a grey relational coefficient is calculated to express the relationship between the ideal and actual normalized experimental results.The grey relational coefficient can be expressed as follows: The average of the grey relational coefficient is then calculated to obtain the grey relational grade.The grey relational grade is defined as follows: However, the effect of each factor on the system is not exactly the same in real applications.Thus, ( 9) can be modified as follows: where   represents the normalized weighting value of factor . Given the same weights, ( 9) and ( 10) are equal.In GRA, the grey relational grade is used to show the relationship among sequences.If two sequences are identical, the value of the grey relational grade is equal to one.The grey relational grade also indicates the degree of influence that the comparability sequence can exert over the reference sequence.Therefore, if a particular comparability sequence is more important than the other comparability sequences to the reference sequence, the grey relational grade for that comparability sequence and reference sequence will be higher than other grey relational grades [27].In this study, the corresponding weighting values   are obtained from PCA.

PCA.
PCA was developed by Pearson [28] and Hotelling [29].This approach explains the structure of variancecovariance by the linear combinations of each quality characteristic.The procedures are described as follows.
(4) Principal Components.The uncorrelated principal component is formulated as follows: where  1 is the first principal component,  2 is the second principal component, and so on.The entire technical line of the hybrid Taguchi/GRA/PCA process optimization method for plastic injection moulding is summarized and illustrated in Figure 1.

Implementation of the Proposed Hybrid
Taguchi/GRA/PCA Optimization Procedures for Plastic Gear validity of the proposed hybrid optimization method.With the capability to run without lubrication and corrosion, plastic gears have been widely used in the automotive industry, office machines, household utensils, food and textile machinery, and a host of other applications [30].Plastic-moulded gears offer some material properties that are not achievable with metal-based gears, including unique advantages in weight, noise, modulus, self-lubrication, chemical resistance, and low cost.However, the practicability of injection moulding in producing low-cost plastic-moulded gears remains restricted by the occurrence of shrinkage in the final parts.
As one of the most common and prominent defects of plastic, severe shrinkage leads to the deflection or warpage of the moulded part and negatively affects the dimensional stability and accuracy of the involute profile, concentricity, roundness, tooth spacing uniformity, and gear size.Sever shrinkage defect is often forgotten and regarded lightly in many works related to plastic-moulded gears.However, this defect critically affects the quality of the final moulded gear, noise, and vibration and shortens the service life of the gear because of different damage mechanisms, such as tooth fatigue, creep, excessive wear, and plastic deformation.Hence, shrinkage minimization is the optimization objective of this paper.This paper focuses on the optimization of processing parameters via the integration of the Taguchi/GRA/PCA to minimize the shrinkage behaviour of plastic injection-moulded gears.To inspect shrinkage behaviour on the dimensional stability of the moulded gear, three quality characteristics are selected: tooth thickness, addendum circle, and dedendum circle.

Selection of Taguchi OA.
In this study, the crystalline thermoplastic polypropylene (PP) is specified for the gear.
The PP used in this study is obtained from Propylene Malaysia Sdn.Bhd., Malaysia.The general properties of PP are shown in Table 1.The spur gear design used is compliant with the American Gears Manufacturers Association standards.The detailed geometry and specifications of the gear are shown in Figure 2.
The PP gears are injected by using a Battenfeld TM750/210 injection moulding machine.The experiment is conducted with four controllable three-level processing parameters: melt temperature, packing pressure, packing time, and cooling time.Other processing parameters, such as mould temperature (25 ∘ C), injection pressure (80 bar), and stroke distance (60 mm), are kept constant during the experiment.The selected processing parameters and their levels are shown in Table 2.The L 9 (3 4 ) OA (Table 3) is conducted to study the four processing parameters.Nine trials of PP gears with five repetitions are produced by using the OA.

Shrinkage Measurement. Rax Vision DC 3000 Mitutoyo
profile projector is used to inspect the accuracy of the specific profile of the involute gear teeth after moulding.In this study, the profile projector is used to measure 2D tooth thickness, as well as addendum and dedendum circles, by using the coordinates of selected points along the gear profile.With large magnifications and micrometer readouts,  the profile projector can ensure fairly accurate measurements compared with a Vernier caliper and micrometer.In the process of addendum and dedendum circle measurement, 20 edge points of each tooth gear are set in the profile projector.
To ensure the integrity of the comparison procedure, 5 injection-moulded gears from the same batch are measured to determine the repeatability of the part geometry.The relative shrinkage of the selected quality characteristics is calculated by using the following equation: where  is shrinkage,  is the reading of the diameter measurement by using the profile projector, and   is the mould dimension.

S/N Analysis.
For constant processing parameters, the average value of five repeated results for shrinkage behaviour in tooth thickness, addendum circle, and dedendum circle is calculated and is considered the final result.In processing optimization, all the results of shrinkage behaviour in tooth thickness, addendum circle, and dedendum circle are transformed into S/N ratio.The conversion of the results into S/N ratio involves a series of calculations of the mean squared deviation.For this case, the smaller-the-better category is used to characterize the shrinkage behaviour of tooth thickness, addendum circle, and dedendum circle (see (1)).The final measured results and S/N ratios for the three quality characteristics are shown in Table 4.

Grey Generation of Raw Data.
In the GRA, the experimental results for the S/N ratios of shrinkage behaviour in tooth thickness, addendum circle, and dedendum circle in Table 4 are first normalized according to the-smaller-thebetter characteristic of the sequence by using (5).The values of tooth thickness, addendum circle, and dedendum circle are set as the reference sequence  () 0 (),  = 1,2,3 and the comparability sequences  ()   (),  = 1, 2, 3, . . ., 9,  = 1, 2, 3. Table 5 lists all the sequences after data preprocessing.According to Deng (1989) [19], a larger value of the normalized results corresponds to better performance and the maximum normalized results that are equal to one indicate the best performance.
The same calculating method is performed for  = 1, . . ., 9, and the results of all Δ 0 for  = 1, . . ., 9 are listed in Table 6.By investigating the data presented in Table 6, Δ max () and Δ min () can be expressed as follows:

Computation of the Grey Relational Coefficient of Response
Variables.The grey relational coefficients for each quality characteristic have been calculated by substituting the distinguishing coefficient  = 0.5 by using (8).Examples on grey relational coefficient  1() are provided as follows:  Thus,  1() = (0.3333, 0.7991, 0.9536),  = 1, 2, 3.A similar procedure is applied for  = 1, . . ., 9. Table 7 lists the grey relational coefficient for each trial of the L 9 OA.

Computation of the Contribution of Respective Quality
Characteristics by Using PCA.In optimizing a problem concerning multiple quality characteristics or performances, an engineering judgment or subjective estimation is needed to determine the weighting values for each quality characteristic.The use of the conventional method to determine the values for each quality characteristics is heavily reliant on experience and trial-and-error, thus resulting in an increase in uncertainty during the decision-making process.To reveal the relative importance for each quality characteristic in GRA objectively, the PCA is introduced.The PCA is adopted to determine the corresponding weighting values for each quality characteristic.The elements of the array for the multiple quality characteristics listed in Table 7 represent the grey relational coefficient of each quality characteristic.These data are used to evaluate the correlation coefficient matrix and to determine the corresponding eigenvalues from (13) (Table 8).The eigenvector corresponding to each eigenvalue is listed Chinese Journal of Engineering in Table 9, and the square of the eigenvector can represent the contribution of the corresponding quality characteristic to the principal component.The contribution of shrinkage behaviour in the tooth thickness, addendum circle, and dedendum circle of the PP-moulded gear is shown in Table 10; these contributions are listed as 0.2878, 0.3224, and 0.3898, respectively.Moreover, the variance contribution for the first principal component characterizing the three quality characteristics is as high as 71.34%.Therefore, for this study, the squares of the respective eigenvectors are selected as the weighting values of the related quality characteristic.
By using the same procedure, the grey relational grade of the comparability sequence for  = 1-9 can also be obtained and is presented in Table 11.The processing parameters were optimized with respect to a single grey relational grade rather than complicated multiple quality characteristics.

Optimal Combination of Injection Moulding Processing
Parameters and Their Levels.To determine the optimal combination of injection moulding processing parameters for shrinkage behaviour in the tooth thickness, addendum circle, and dedendum circle of the studied moulded gear, the average grey relational grade for each injection moulding parameter level is calculated by employing the main effect analysis of the Taguchi method.This process is performed by sorting the grey relational grades corresponding to the levels of the injection moulding parameters in each column of the OA and then taking the average of parameters with the same levels.For instance, for factor A (Table 3), experiments 1, 2, and 3 are set to level 1.Therefore, by using the data listed in Table 11, the average grey relational grade for  1 can be calculated as follows: The average grey relational grade for  2 and  3 is calculated as follows: By using a similar method, calculations are performed for each injection moulding parameter level, and the main effect analysis is constructed (Table 12 and Figure 3).
Considering that the grey relational grade represents the level of correlation between the reference and comparability sequences, a larger grey relational grade indicates that the comparability sequence exhibits a stronger correlation with the reference sequence.A larger grey relational grade results in better multiple quality characteristics.Figure 3 clearly shows that the multiple quality characteristics of the PPmoulded gear are significantly affected by the adjustments of the processing parameters.The increment of melt temperature, packing pressure, and cooling time initially increases the grey relational grades, thus resulting in the reduction of shrinkage behaviour in the tooth thickness, addendum circle, and dedendum circle of the moulded gear.However, the relational grades for the three quality characteristics are lower when the melt temperature, packing pressure, and cooling time are 240 ∘ C, 100%, and 50 s, respectively.By contrast, the increment of packing time reduces the value of the grey relational grade, thus increasing the shrinkage behaviour in tooth thickness, addendum circle, and dedendum circle.In this case, the best combination of processing parameters and levels can easily be obtained from the main effect analysis by selecting the level of each parameter with the highest grey relational grade.A 2 , B 2 , C 1 , and D 2 show the largest value of grey relational grade for factors A, B, C, and D, respectively (Figure 3).Thus, the optimal parameter setting that statistically results in the minimum shrinkage behaviour in the tooth thickness, addendum circle, and dedendum circle of the PP-moulded gear is predicted to be A 2 , B 2 , C 1 , and D 2 .The melt temperature is 220 ∘ C, packing pressure is 80%, packing time is 5 s, and cooling time is 40 s.

Effects of Processing Parameters on Quality Characteristics.
To examine the extent in which injection moulding parameters significantly influence the performance of moulded gear, ANOVA is performed on the Taguchi method for the grey relational grade of nine comparability sequences (Table 11).The computed quantity of degrees of freedom (DOF), sum of square (), variance (), and percentage contribution () are presented in Table 13.
In this case, the percentage contribution of each processing parameter is directly calculated from S because the DOF for the error term is equal to zero [31].The significance of each processing parameter in the shrinkage behaviour in the tooth thickness, addendum circle, and dedendum circle of the moulded gear can be determined by the percentage contribution.Roy [32] suggested an alternative by using the 10% rule; that is, a parameter is considered insignificant when its influence is less than 10% of the highest parameter influence.From the results of ANOVA in Table 13, melt temperature appears to be the most decisive processing parameter in reducing the shrinkage behaviour in the tooth thickness, addendum circle, and dedendum circle of the moulded gear with the highest percentage contribution of 71.47%, thus outweighing the other process variables.The analysis also reveals that packing pressure and cooling time are significant because their percentages are more than 10% of the highest parameter influence (7.15%).Cooling time and packing pressure achieve 14.92% and 13.06%, respectively.Packing time results in only 0.55%, which is significantly less than 10% of the highest parameter influence.In this case, the packing time is considered insignificant in the shrinkage behaviour of the tooth thickness, addendum circle, and dedendum circle of the moulded gear.

Verification Test.
Once the optimal levels of the injection moulding process parameters are identified, the subsequent step is to verify the improvements in the quality characteristics by using this optimal combination.The verification test can be used to assess the accuracy of the proposed hybrid Taguchi/GRA/PCA optimization method.An experimental verification test is conducted by using the same procedures as previous runs under the optimal process conditions, namely, A 2 , B 2 , C 1 , and D 2 , to produce the plastic-moulded gear.Table 14 lists the results of five repetitions of the verification tests by using the optimal process conditions obtained by the proposed hybrid optimization method.After optimization, the minimum shrinkage behaviour in tooth thickness, addendum circle, and dedendum circle is averaged as 1.8%, 1.53%, and 2.42%, respectively.The shrinkage behaviour for the three quality characteristics decreases with respect to those attained in the main experiment presented in Table 4.

Conclusion
In this study, a hybrid Taguchi/GRA/PCA optimization method for injection moulding process parameters has been developed systematically to overcome the shortcomings of individual methods in multiple quality characteristics problems.The problem that arises when calculating the weighting value for each quality characteristic in GRA has been addressed by integrating the PCA to determine the grey relational grades.Thus, the conventional approach of using engineering judgment or subjective estimation to determine the weighting values for each quality characteristics can be avoided.A plastic gear is employed in this study as a case study to demonstrate the efficiency and validity of the proposed hybrid optimization method in controlling all influential injection moulding processing parameters during plastic gear manufacturing.The shrinkage behaviour in tooth thickness, addendum circle, and dedendum circle is chosen as the optimization objectives.Through a series of analyses and optimizations of selected multiple quality characteristics for the case of a plastic gear, the minimum shrinkage behaviour in tooth thickness, addendum circle, and dedendum circle after optimization by the proposed hybrid optimization method is averaged as 1.8%, 1.53%, and 2.42%, respectively, which are less than the values in the main experiment.

Figure 3 :
Figure 3: Effect of injection moulding parameter levels on multiple quality characteristics.

Table 1 :
General properties of PP.
3.1.Determination of Quality Characteristics.In this paper, a plastic gear is used to demonstrate the efficiency and

Table 2 :
Processing parameters and levels studied.

Table 4 :
Average shrinkage and S/N ratios for each trial.

Table 7 :
Calculated grey relational coefficient for nine comparability sequences.

Table 8 :
Eigenvalues and explained variation for principal components.

Table 9 :
Eigenvectors for principal components.

Table 10 :
Contribution of each individual quality characteristic for the principal component.

Table 11 :
Grey relational grade and its order.

Table 12 :
Main effects table for grey relational grades.

Table 13 :
ANOVA table for the grey relational grade of nine comparability sequences.

Table 14 :
Results of the verification test.