# Minimum Porosity Formation in Pressure Die Casting by Taguchi Method

Authors: Quang-Cherng Hsu; Anh Tuan Do
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

Die casting process is significantly used in the industry for its high productivity and less postmachining requirement. Due to light weight and good formability, aluminum die casting plays an important role in the production of transportation and vehicle components. In the current study of die casting for Automobile starter motor casing, the following issues are focused: shot piston simulation, defect analysis, and, finally, the use of the Taguchi multiquality analytical method to find the optimal parameters and factors to increase the aluminium ADC10 die casting quality and efficiency. Experiments were conducted by varying molten alloy temperature, die temperature, plunger velocities in the first and second stage, and multiplied pressure in the third stage using L27 orthogonal array of Taguchi method. After conducting a series of initial experiments in a controlled environment, significant factors for pressure die casting processes are selected to construct an appropriate multivariable linear regression analysis model for developing a robust performance for pressure die casting processes. The appropriate multivariable linear model is a useful and efficient method to find the optimal process conditions in pressure die casting associated with the minimum shrinkage porosity percent.


---

## Body

## 1. Introduction

High pressure die casting for nonferrous casting applications is increasingly used in the foundries today as an economically viable casting process. High pressure die casting (HPDC) process has been widely used to manufacture a large variety of products with high dimensional accuracy and productivities. It has a much faster production rate in comparison to other methods and it is an economical and efficient method for producing components with low surface roughness and high dimensional accuracy. All major aluminium automotive components can be processed with this technology.High Pressure Die Casting process is rapid and depends on many factors. So, to capture the problem it requires a lot of time and experience including testing and simulation. The conventional trial and error based die design and process development is expensive and time consuming. Such a procedure also might lead to higher casting rejections. The HPDC castings production process has many defects, such as shrinkage porosity, misrun, cold-shut, blister, scab, hot-tear. Several previous studies of defects in aluminum alloy by the method of HPDC and disability solutions (Shen et al. 2007 [1], Dargusch et al. 2006 [2], Verran et al. 2006 [3], Mousavi Anijdan et al. 2006 [4], Tsoukalas et al. 2004, 2008 [5, 6]). However, the study to optimize aluminum alloy casting process in the condition of production casting factory is essential. This study focused on analysis of shrinkage porosity defect with mold design and put into production casting by foundry factory conditions.Shrinkage porosity is one of the most common defects leading to rejection of aluminium die casting, often only showing up after much value has been added to the casting via operations such as machining, polishing, and coating. The added value of the casting at the point of rejection can be very high. If you find out the causes and how to reduce the defects of castings will be of great significance in reducing the production cost of die casting. However, optimizing the conditions to render aluminium die castings of minimum porosity percent is costly and time consuming, because many experiments are necessary to find the optimal parameters.Taguchi method is one of the efficient problems solving tools to upgrade the performance of products and processes with a significant reduction in cost and time involved. Taguchi’s parameter design offers a systematic approach for optimization of various parameters with regard to performance, quality, and cost (Syrcos 2003 [7], Taguchi 1986 [8]).


## 2. Materials and Methods

The die casting part product of this study is provided through aluminium die casting factory, so the casting body no changes. A major factor in the successful development of castings is the design of the die and design of gates, biscuit, and runner system. A well-designed gating and runner system should avoid turbulence in metal flow and to reduce incidence of inclusions and air entrapment in the casting. The die design is required to avoid solidification related defects like shrinkage, micro-porosities, hot-tear and so forth. Die design process is very much dependent on the experience and skill of the design engineer. The die for this study is the result of collaboration between the foundry factory and Department of Mechanical Engineering-National Kaohsiung University of Applied Sciences. The casting with full of the gating, runner system and biscuit, is shown in Figure1. The die casting is designed in CATIA V5R19 software, shown in Figure 2. Moreover, the die casting material selection is very important. The nature of the material will directly affect the quality of the casting and die casting parameters configuration, this study selects casting material as the aluminium alloy ADC10. The chemical composition of the aluminum alloy used in the experimental procedure is given in Table 1.Table 1
Chemical composition of the alloy ADC10 used in the experiment.
Element
Si
Fe
Cu
Mg
Mn
Ni
Zn
Sm
wt%
7.5~9.5
1.3
3.0~4.0
0.1
0.5
0.5
3
0.35Casting image.
(a)
(b)
(c)
(d)Part product is designed by CATIA software.
(a)
(b)
(c)
(d)Shrinkage porosity formation in pressure die casting is the result of a so much number of parameters. Figure3 shows a cause and effect diagram that was constructed to identify the casting process parameters that may affect die casting porosity (Tsoukalas et al. 2004, 2008 [5, 6]). In this case, holding furnace temperature, die temperature, plunger velocity in the first stage, plunger velocity in the second stage, and multiplied pressure in the third stage were selected as the most critical in the experimental design. The other parameters were kept constant in the entire experimentation. The range of holding furnace temperature was selected as 640~700°C, the range of die temperature as 180~260°C, the range of plunger velocity in the first stage as 0.05~0.35 m/s and in the second stage as 1.5~3.5 m/s, and the range of multiplied pressure in the third stage was chosen as 200~280 bars. The selected casting process parameters, along with its ranges, are given in Table 2.Table 2
The parameter and its value​​at three levels.
Process parameters
Parameters range
Level 1
Level 2
Level 3
Holding furnace temperature (°C)
640~700
640
670
700
Die temperature (°C)
180~260
180
220
260
Plunger velocity, 1st stage (m/s)
0.05~0.35
0.05
0.2
0.35
Plunger velocity, 2nd stage (m/s)
1.5~3.5
1.5
2.5
3.5
Multiplied pressure (bars)
200~280
200
240
280Figure 3
Cause and effect diagram.Taguchi method based design of experiment has been used to study the effects of five casting process parameters (holding furnace temperature: A, die temperature; B, plunger velocity in the first stage; C, plunger velocity in the second stage; D, multiplied pressure in the third stage; E, on an important output parameter (Shrinkage porosity). For selecting appropriate orthogonal array, degree of freedom (number of fair and independent comparisons needed for optimization of process parameters is one less than the number of levels of parameter) of the array is calculated.In the experimental layout plan with five factors and three levels using L27 orthogonal array, 27 experiments were carried out to study the effect of casting input parameters, shown in Table 3. The input parameters are installed in the ProCAST software to conduct 27 simulation experiments.Table 3
Experimental layout using an L27 orthogonal array.
Trials
Holding furnace temperature A
Die  temperatureB
Plunger velocity 1st stage C
Plunger velocity 2nd stage D
Multiplied pressure E
1
1
1
1
1
1
2
1
1
2
2
2
3
1
1
3
3
3
4
1
2
1
2
2
5
1
2
2
3
3
6
1
2
3
1
1
7
1
3
1
3
3
8
1
3
2
1
1
9
1
3
3
2
2
10
2
1
1
2
3
11
2
1
2
3
1
12
2
1
3
1
2
13
2
2
1
3
1
14
2
2
2
1
2
15
2
2
3
2
3
16
2
3
1
1
2
17
2
3
2
2
3
18
2
3
3
3
1
19
3
1
1
3
2
20
3
1
2
1
3
21
3
1
3
2
1
22
3
2
1
1
3
23
3
2
2
2
1
24
3
2
3
3
2
25
3
3
1
2
1
26
3
3
2
3
2
27
3
3
3
1
3Computer simulation procedure-based process development and die design can be used for rapid process development and die design in a shorter time. Such a computer simulation based procedure, often using FINITE ELEMENT ANALYSIS based software systems, can improve the quality and enhance productivity of the enterprise by way of faster development of new product. Analysis software is used as a ProCAST commercial with finite element method analysis for a casting process. In this study, all parameters can be able to affect the analysis process, choice of material is aluminum alloy die casting ADC10, and cold chamber die casting method with molding material is H13. FEM based simulation software systems help the designer to visualize the metal flow in the die cavity, the temperature variation, the solidification progress, and the evolution of defects such as shrinkage porosity, cold-shut, hot-tear.ProCAST a FEM simulation-based virtual casting environment for analysis of the casting process is used as a tool for die design and process optimization. ProCAST with Visual-Viewer module can provide temperature field, thermal cracking, flow field, solidification time, and shrinkage analysis. This paper focused on the analysis of shrinkage porosity by ProCAST software base on parameters input from Table3.The analysis of defects simulated by ProCAST software with Visual-Viewer module can detect many types of disabilities casting. The defective products do not necessarily reflect the loss of the original function, for example, the internal pore trims acceptable. However, with large structural castings, defect analysis of this study focuses on maximum porosity in the selection casting, and the important parts of the casting shrinkage analysis (an important component), casting defect analysis are described as follows.The Solid Fraction. Solid fraction may be available shrinkage prediction casting position, the present study is in accordance with the theory prediction of defect, and ProCAST manual referred to in the final period of solidification. Shrinkage solid fraction prone is greater than 0.7, here as the reference value of 0.7 solid fractions. When the solid fraction area is below this value and the area around the solid phase rate is rather than this value, we can predict this area shrinkage porosity occurred.Maximum Porosity. The maximum porosity analysis using the Shrinkage Porosity function of the Visual-Viewer comes defined in the manual. According to the ProCAST user manual shrinkage definition and with the solid fraction, it can be used to analyse the basis of the maximum porosity.Shrinkage Analysis. For the amount of inspection shrinkage casting part used for the Visual-Viewer module function for quantitative analysis. In each experiment we took five elements with the coordinates determined at the important positions in the working conditions of automobile starter motor casing. Each experiment was repeated five times in order to reduce experimental errors, as shown in Figure 4. Data from 27 experiments with five sampling times in each simulation are summarized as in Table 4. From this table, we conducted quality characteristics analysis.Table 4
Shrinkage porosity results of the L27 array design. (Full table).
Trials
Holding Furnace Temperature A
Die temperature B
Plunger velocity1st stage C
PlungerVelocity2nd stage D
MultipliedPressure E
Shrinkage porosity (%)
Average
MSD
S
/
N
Repetition
Repetition
Repetition
Repetition
Repetition
1
2
3
4
5
1
640
180
0.05
1.5
200
1.973
1.988
1.984
1.978
1.954
1.9754
3.902346
−5.91326
2
640
180
0.2
2.5
240
1.912
1.893
1.948
1.908
1.912
1.9146
3.666021
−5.64195
3
640
180
0.35
3.5
280
1.846
1.853
1.814
1.813
1.828
1.8308
3.352095
−5.25316
4
640
220
0.05
2.5
240
1.862
1.867
1.861
1.878
1.871
1.8678
3.488716
−5.42666
5
640
220
0.2
3.5
280
1.793
1.804
1.761
1.764
1.775
1.7794
3.166541
−5.00585
6
640
220
0.35
1.5
200
1.916
1.908
1.912
1.887
1.903
1.9052
3.629888
−5.59893
7
640
260
0.05
3.5
280
1.786
1.797
1.754
1.757
1.768
1.7724
3.141679
−4.97162
8
640
260
0.2
1.5
200
1.909
1.903
1.889
1.898
1.893
1.8984
3.603973
−5.56782
9
640
260
0.35
2.5
240
1.852
1.859
1.82
1.819
1.836
1.8372
3.375568
−5.28347
10
670
180
0.05
2.5
280
1.799
1.81
1.767
1.77
1.781
1.7854
3.18793
−5.03509
11
670
180
0.2
3.5
200
1.957
1.968
1.971
1.951
1.965
1.9624
3.851068
−5.85581
12
670
180
0.35
1.5
240
1.857
1.864
1.825
1.824
1.841
1.8422
3.393965
−5.30707
13
670
220
0.05
3.5
200
1.924
1.935
1.923
1.928
1.909
1.9238
3.701079
−5.68328
14
670
220
0.2
1.5
240
1.803
1.814
1.771
1.774
1.785
1.7894
3.202229
−5.05452
15
670
220
0.35
2.5
280
1.751
1.744
1.713
1.731
1.725
1.7328
3.002778
−4.77523
16
670
260
0.05
1.5
240
1.796
1.807
1.764
1.767
1.778
1.7824
3.177227
−5.02048
17
670
260
0.2
2.5
280
1.738
1.726
1.724
1.745
1.715
1.7296
2.991629
−4.75908
18
670
260
0.35
3.5
200
1.881
1.878
1.898
1.893
1.903
1.8906
3.574461
−5.53211
19
700
180
0.05
3.5
240
1.876
1.871
1.895
1.887
1.893
1.8844
3.551052
−5.50357
20
700
180
0.2
1.5
280
1.755
1.743
1.741
1.762
1.732
1.7466
3.050725
−4.84403
21
700
180
0.35
2.5
200
1.921
1.901
1.956
1.916
1.928
1.9244
3.703644
−5.68629
22
700
220
0.05
1.5
280
1.708
1.696
1.694
1.715
1.685
1.6996
2.888753
−4.6071
23
700
220
0.2
2.5
200
1.864
1.898
1.881
1.876
1.885
1.8808
3.537532
−5.487
24
700
220
0.35
3.5
240
1.825
1.832
1.793
1.792
1.807
1.8098
3.275642
−5.15296
25
700
260
0.05
2.5
200
1.856
1.861
1.855
1.872
1.865
1.8618
3.466338
−5.39871
26
700
260
0.2
3.5
240
1.814
1.821
1.782
1.781
1.796
1.7988
3.235948
−5.10001
27
700
260
0.35
1.5
280
1.692
1.694
1.706
1.683
1.712
1.6974
2.881274
−4.59585Casting measurement area.
(a)
(b)
(c)
(d)Quality Characteristics. The parameter design study involves control and noise factors. The measure of interactions between these factors with regard to robustness is signal-to-noise (S/N) ratio. S/N characteristics formulated for three different categories are as follows: the bigger the better, and the smaller the better, the nominal the best. This paper focused on studying the effects of five input parameters (A,B,C,D,E) to defect shrinkage porosity in the process of casting, so the criteria “the smaller the better” is selected.The smaller the better (for making the system response as small as possible) is as follows:(1)SNS=-10log(1n∑i=1nyi2),
where n is the number of sampling (Each experiment was repeated five times sampling, so that n=5), yi: value of Shrinkage porosity at each time sampling.The responding graph shown in Figure5 learned that the best combination for this study with shrinkage porosity defect value minimum is A3B3C3D1E3.Figure 5
S
/
N Response graphs.Process Parameter Optimization Using MVLR. The objective of the process optimization is to select the optimal control variables in aluminium die casting process in order to obtain the minimum porosity. In this work, the fitness function used in the optimization procedure was based on the MVLR model.Multivariable Linear Regression Analysis. In most cases, the form of the relationship between the response and the independent variables is usually unknown. Multiple linear regression (MLR) is a method used to model the linear relationship between a dependent variable and one or more independent variables. MLR is based on least squares: the model is fitted such that the sum-of-squares of differences of observed and predicted values is minimized.Letx1;x2;…;xr be a set of r predictors believed to be related to a response variable Y. The linear regression model for the jth sample unit has the form
(2)Yj=β0+β1xj1+β2xj2+⋯+βrxjr+εj,
where ε is a random error and the βi,i=0,1,…,r are unknown regression coefficients.In this paper, there are five independent variables and one dependent variable. The relationships between these variables are of the following form:(3)F(x)=β0+β1A+β2B+β3C+β4D+β5E.In whichF(x): Dependence variable,A (°C): Holding furnace temperature,B (°C): Die temperature,C (m/s): Plunger velocity 1st stage,D (m/s): Plunger velocity 2nd stage,E (bars): Multiplied pressure during the third phase.The results after analysing by Intercooled Stada 8.2 Softwere as shown in Table5.Table 5
The results after analysing by Intercooled Stada 8.2 software.
(a)
. reg FxA B C D E
Source
SS
df
MS
Number of obs = 27F(5,21)=146.74
Model
0.155044498
5
0.0310089
Prob >F = 0.0000
Residual
0.004437799
21
0.000211324
R-squared = 0.9722 Adj R-squared = 0.9655
Total
0.159482297
26
0.006133935
Root MSE = 0.01454
(b)
Fx
Coef.
Std. Err.
t
P
>
|
t
|
[95% Conf. Interval]
A
−0.0008844
0.0001142
−7.74
0.000
−0.001122
−0.0006469
B
−0.00083
0.0000857
−9.69
0.000
−0.0010081
−0.0006519
C
−0.0305925
0.7284277
−0.042
0.000
−0.0780965
0.0169114
D
0.0175444
0.0034264
5.12
0.000
0.0104189
0.02467
E
−0.0020122
0.0000857
−23.49
0.000
−0.0021904
−0.0018341
_cons
3.054569
0.0820708
37.22
0.000
2.883893
3.225244The final MVLR model equation for porosity after substituting regression coefficients is as follows:(4)F(x)=3.054569-0.8844*10-3A-0.83*10-3B-0.03059C+0.01754D-0.00201E.


## 3. Results and Discussion

We applied multivariable linear regression analysis (MVLR) to seek the optimal parameter in the casting process of independent parameter variables in this study. A stationary point for the optimal performance was obtained by using the multivariable linear regression method in this linear regression equation, and the result is presented in Figure6. A very good fit was observed and substantiated by the coefficient of determination, R2=0.9722. That is, the R2 value indicates that the polynomial model explains almost 97.22% of variability in the casting process.Figure 6
Experimental and predicted values of shrinkage porosity percent.Figure6 shows the efficacy of the optimization scheme by comparing the MVLR results with the experimental values. There is a convincing agreement between experimental values and predicted values for shrinkage porosity percent.Matlab code for finding optimization shrinkage porosity value.Program in Matlab (see Algorithm1).Algorithm 1
c
lc;
clear all;
close all;
f =
@
(x) 3.054569-0.8844e-3*x(1)-0.83e-3*x(2)-...
0.03059*x(3)   +   0.01754*x(4)-0.00201*x(5);
options = optimset(‘GradObj', ‘on');
[x,fval,exitflag,output] =
…
fmincon(f,[670;220;0.2;2.5;240],
[
]
,
[
]
,
[
]
,
[
]
,[600;180;0.05;1.5;
200],[700;260;0.35;3.5;280],
[
]
,optimset(‘Display', ‘iter’));
x
fvalResults after running in Matlab are as follows:(5)x=700.0000sothat⟶A=700°C260.0000B=260°C0.3500C=0.35m/s1.5000D=1.5m/s280.0000E=280barfval=1.6725Shrinkageporosity:1.6725%.By the Program in Matlab, we are known as the best combination in the 27 experimental configurations.This result is similar to quality characteristics and is the best combination for this studyA3B3C3D1E3.


## 4. Conclusion

In this paper, the optimum process parameters values predicted for casting of minimum shrinkage porosity (1.6725%) and the best combination parameters given as follows:holding furnace temperature 700°C,die temperature 260°C,plunger velocity, 1st stage 0.35 m/s,plunger velocity, 2nd stage 1.5 m/s,multiplied pressure 280 bar.The model proposed in this paper gives satisfactory results in the optimization of pressure die casting process. The predicted values of the process parameters and the calculated are in convincing agreement with the experimental values.The experiments which are conducted to determine the best levels are based on “Orthogonal Arrays,” and are balanced with respect to all control factors and yet are minimum in number. This in turn implies that the resources (materials, saving time, and money) required for the experiments are also minimized.



---

*Source: 920865-2013-12-04.xml*