1. Introduction
Plastic injection molding technology has been regarded as one of the most important technologies for rapid, efficient manufacturing. During the injection molding process, a plastic melt is injected into the void cavity by the screw of the injection molding machine, and subsequently solidified in the cavity after a certain cooling time. It is a complicated multi-physical process involving a thermal and mechanical evolution history. The level of molding defects, like warpage or shrinkage magnitude, length of weld line, etc., are associated with this complicated thermal and mechanical history. Therefore, it has attracted researchers’ attention to investigate the thermal and mechanical behavior in the injection molding process, and relevant process optimization methods to improve the molding quality for a long time.
In fact, the thermal and mechanical behavior of the plastic material depends on the process parameters, molding temperature, pressure, etc. As a special injection molding technology, dynamic injection molding technology (DIMT) integrates melt vibration technology [
1] and injection molding technology. It attempts to introduce the vibrational force field into plasticizing, filling, and packing stages to disturb the flow behavior [
2,
3], crystallization behavior, etc. The final internal stress field may be more proper to eliminate the defects of the product. The mechanism is that the thermal and mechanical evolution history is changed due to this external vibration.
Many studies have introduced the research work of the DIMT effect on the mechanical or electronic properties of plastics from the experimental perspective. For instance, Kikuchi et al. [
4] improved the mechanical strength of plastics by using DIMT, and observed the alternation of residual stress distribution of the injected samples from the birefringence pattern. Yan et al. [
5] observed the shear stress and apparent viscosity reduction with DIMT. Yang et al. [
6] undertook detailed research on the electrical conductivity of the carbon nanotube/polypropylene composite with the help of DIMT, and observed the perfect conductivity paths under DIMT which can improve the electrical conductivity. Additionally, Yin et al. [
7] investigated the pressure response and crystallization behavior, and found that the crystallization is increased under DIMT. Wang et al. [
8] improved the properties of CaCO
3/PP composite material by using DIMT in plasticizing. However, there is a lack of literature on DIMT from a simulation point of view to investigate the thermal and mechanical evolution history, revealing a mechanism for improving the molding quality of plastic products.
Perfect molding quality cannot be obtained without the proper molding process parameters. For thin-walled products, warpage is one main molding defects. The main work of many studies is to find the optimum manufacturing process conditions to minimize the warpage. The method is to find a bridge to combine the computer aided engineering (CAE) technology with the mathematical optimization theory. For instance, Ozcelik and Erzurumlu [
9] used artificial neural network (ANN)- and genetic algorithm (GA)-based optimization methods to minimize the warpage of the product. Kitayama et al. [
10] reduced the warpage of the product by optimizing the packing pressure curve. Zhou and Turng [
11] proposed a process parameters optimization method on warpage and volumetric shrinkage with the help of a Gaussian process surrogate model. Wu, Ku, and Pai [
12] reduced the warpage of the injection molding product with weld line constraints. Gao and Wang [
13,
14] employed an efficient adaptive optimization method to reduce the warpage defect, which is based on a kriging surrogate model and expected improvement (EI) function. Shi et al. [
15,
16] used an ANN model to optimize the warpage of the product. Farshi et al. [
17] minimized both the warpage and the volumetric shrinkage by using the sequential variable-size simplex algorithm. Moreover, considering multi-class design variables, including process parameters, geometry information of the gate location, and thickness of the part, Deng et al. [
18] employed a particle swarm optimization (PSO) algorithm for solving the multi-objective problem including warpage, weld line, and air traps. Zhao et al. [
19] solved the multi-objective optimization problem of an injection molding product, including warpage, shrinkage and sink marks. However, this literature is all on the basis of the conventional injection molding technology (CIMT).
In our previous work [
20], we proposed a numerical method to simulate the process of DIMT and performed the process parameter optimization to minimize the warpage of a typical product. This paper is aimed at trying to extend our previous work to make an in-depth and detailed study on the mechanical evolution behavior and reveal the advantage and disadvantage for improving the molding quality. The implementation of numerical DIMT by using the finite element method is introduced with the given assumptions and solution conditions in the subsequent section. Then, pressure response analysis is performed in order to discuss the influence of dynamic parameters on the pressure field. This can help us to realize the mechanical evolution rules and clarify the mechanism for improving the molding quality. The most significant thing is to further argue the actual effect on warpage reduction by using a kriging model-based sequential optimization method. It is meaningful for engineers to recognize the features of DIMT when they produce products with DIMT. These are the differences and innovative aspects of this paper compared with our previous study of the literature [
20].
2. Flow Analysis Computation Method of Dynamic Injection Molding Technology (DIMT)
2.1. Fundamental Physical Equations
During the whole injection molding processing, the flow process is the early stage, which can determine the molding quality in the late stage. It can be divided into the filling phase and packing phase. In the filling phase, the plastic melt is injected into the void mold cavity. If the mold cavity is almost filled with the plastic material, the packing process proceeds. During the packing phase, there is still a certain amount of plastic material being injected into the cavity, which ensures filling the gap due to material shrinkage from the melt state to solid state. Although the molding process of the product is a complicated physical process, it should also satisfy the fundamental physical governing equations.
Equations (1)–(3) are the complete governing equations to describe the physical change of the plastic melt during injection molding, including the change of mass, momentum and energy.
where
ρ,
cp,
β, and
k represent four physical properties, density, specific heat, coefficient of thermal expansion, and thermal conductivity, respectively.
v and
f represent the velocity field and body force per unit mass, respectively.
σ is the stress tensor that can be decomposed into an isotropic part and a non-isotropic part, as in Equation (4).
is the heat generation rate per unit mass.
κ is the expansion viscosity that can determine the heat generation level during expansion.
T and
pe represent the temperature and equilibrium pressure, respectively.
Obviously, it is difficult to solve the complete equations above, so we have to make some physical assumptions for simulating the plastic injection molding process. The main assumptions are:
- (a)
For thin-walled injection molded products, the flow behavior of the melt can be regarded as the general Hele-Shaw flow, which neglects the “fountain flow” phenomenon at the front of the melt. The finite mesh type is the mid-plane mesh with triangular elements.
- (b)
Considering the high viscosity effect of the plastic melt, inertia force and gravity are neglected. Therefore, the flow process of the melt is treated as a laminar flow, and regarded as the general Newtonian fluid. The melt is dominated by shear stress.
- (c)
No-slip assumption. This means that there is no relevant movement between the melt and the mold cavity surface, or between the melt and the solid plastic material surface.
- (d)
For the energy transport issue, there is no heat source in the melt in practice. We only consider heat conduction, but neglect heat convection along the thickness direction. In contrast, we consider heat convection and neglect heat conduction along the flow direction.
According to Assumption (b), the constitutive equation of the melt should coincide with the form of the Newtonian fluid, which reads as:
where
is the rate of strain tensor.
η(
T,
p,
) is the viscosity of the melt, and is a highly non-linear function on temperature
T, pressure
p, and shear rate
, which is defined as:
Many researchers made efforts in order to obtain a perfect model of
η(
T,
p,
) [
21,
22,
23], listed in
Table 1. However, the Cross-WLF model is suitable to be used as a computation model for most plastic materials. In this paper, Cross-WLF model is employed.
By using the assumptions above, the momentum equation and energy equation are simplified as:
In order to solve these governing equations, it should combine with the solution conditions include the conditions of pressure, velocity, continuity and temperature:
By means of integrating the continuous Equation (1) along with the thickness and Galerkin-weighted residual method, we can obtain Equation (13) completely for solving the flow process:
where:
where
ξ is the position between the melt and solid material along the thickness.
In addition, many properties of the plastic material are variable with temperature and pressure, like thermal conductivity and density. For describing the thermal conductivity evolution, we employ:
where
λ1 to
λ5 are constant parameters.
The pressure-volume-temperature relationship, which is usually called the pvT state equation, is the critical rule to characterize the change of density (specific volume) with respect to temperature and pressure. For the pvT state equation, there are several state models to describe the property change, like the Spencer-Gilmore (S-G) model [
24] and the two-domain Tait model [
25]. In this study, we employ the two-domain Tait model, as Equation (15) expresses:
where
C = 0.0894, and:
where
b1 to
b9 are material parameters,
Tt represents the transition temperature, which is a crystalline temperature for crystalline plastic and the glass transition temperature for amorphous plastic. As an explanation of the pvT relationship,
Figure 1 presents the specific volume change of PC Lexan 105 with respect to temperature and pressure.
2.2. The Finite Element Equation for Computing the Pressure Field
For simulating the flow process during injecting, we need to deduce the finite element discrete form of Equation (13). Consider that the pressure in any element is interpolated by:
where index
m indicates the node number of the current element.
Lm is the interpolation function, which satisfies:
where
A is the area of the triangle element. Then, we can obtain the finial Equation (22):
where
K,
P, and
Q represent the stiffness matrix, pressure vector, and the net flow rate vector, respectively. Subscript
G,
M, and
F represent the injecting gate, melt region and the flow front of the melt, respectively. Therefore,
QM = 0 in the melt, and
PF = 0 at the front of the melt.
For the conventional injection molding, flow rate
QG is usually a given constant during filling analysis, and
PG is a given constant or linear form during the packing phase. We can make a slight change of Equation (22):
for computing the pressure and velocity field in filling analysis with
QG given, and:
for computing the pressure and velocity field in packing analysis with
PG given.
For the dynamic filling phase, the flow rate can be deduced by the advancing displacement of the screw which is:
where
vave is the average advancing velocity of the screw,
Samp,
ffill, and
φfill are amplitude, frequency, and phase of a vibration filling. Then, the flow rate is:
where
V is the volume of the product,
Asc here is the area of the screw,
tfill is the approximate filling time.
ς equals (
Samp∙2π∙
ffill∙
tfill∙
Asc)/
V herein.
Although Equation (26) is available to be the dynamic filling process, there is a limitation that, in order to guarantee no melt reflux problems, (Samp∙2π∙ffill∙tfill∙Asc)/V ≤ 1 has to be satisfied. This indicates that high frequency or high amplitude vibration cannot be realized, since if Samp or ffill is higher, ffill or Samp will be so low that the vibration will be “no sense” in practice. Therefore, we prefer to choose another case of dynamic filling rather than Equation (26) in the following study.
Therefore, if we express the advancing velocity of the screw in a simple form as:
where
ε is called the amplitude factor, we can easily obtain the dynamic filling flow rate expression as:
Equation (28) avoids the constraint problem mentioned above, and can be realized easily as well. For the packing phase, making the similar manipulation the filling phase, the packing pressure curve is controlled by:
where
Pfinal represents the finial pressure during the filling phase,
Pstart is the packing pressure fraction at the beginning,
δ is a factor to control the trend of the packing curve,
θ,
fpack, and
φpack represent the amplitude, frequency, and phase of a vibration term, respectively.
Overall, with Equations (28) and (29), Equations (23) and (24) can be solved to investigate the pressure response analysis and the molding quality improvement analysis of plastic product of numerical DIMT, which will be argued in the following sections.
3. Pressure Response Analysis in the Cavity
Defects, like warpage, shrinkage, etc., depend on the mechanical and thermal evolution history. By integrating Equation (7) associated with the solution conditions, we can deduce the velocity formula as:
where
h is the thickness of the product.
Obviously, solving the velocity field requires us to firstly solve the pressure field. Moreover, computing the heat contribution from shear motion relies on the pressure field as well. Therefore, pressure in the melt is a very important mechanical quantity to lead the material motion and the thermal evolution. This indicates that pressure response analysis can help us to understand the mechanical and thermal evolving history.
Therefore, a finite element model of a low-density polyethylene (LDPE) circle disk product (depicted as
Figure 2) is modeled as the investigating objective, which is the same as the product of Yin’s experiment [
7]. The basic material properties of LDPE are listed as
Table 2. The injecting gate is located at the center of the disk, which is colored with pink in the finite element model.
In order to compare the simulation result of DIMT with the experimental observation [
7], three monitor points, except the gate, A, B, and C, are selected in
Figure 2b. The distance of OA, OB, and OC are 10.08, 29.40, and 46.34 mm, respectively. The melt temperature and mold temperature are 190 °C and 50 °C, respectively. Moreover, injecting time and packing time are both 2 s. Other dynamic parameters are listed in
Table 3.
By finite computation, we can obtain the pressure response of these four monitors both for the CIMT case and DIMT case. From
Figure 3 and
Figure 4, it is clarified that, compared with the experimental observation, our pressure results are efficient, although there is a certain difference between the computation and experiment. The value of the simulation results is slightly lower than the experimental value. This difference is mainly because the material property and processing parameters of computation cannot be exactly the same as the actual experiment. However, this indicates that our computation method can be used as a simulation tool for investigating the mechanical evolution and the product quality under DIMT further.
There are several interesting phenomena in
Figure 4. The first one is that the pressure response in the cavity has a fluctuation compared with the CIMT case. The second one is that this fluctuation amplitude gets larger and larger as time goes by for each monitor. The third one is that, although the observation of the pressure response rule for different monitors is the same, the further the monitor location is away from the gate, the smaller the fluctuation amplitude is. The fourth one is that, if we observe the fluctuation frequency, we find that this frequency seems like the one of the dynamic flow rate. If we make a fast Fourier transform (FFT) operation to this fluctuation, we can obtain the frequency domain as
Figure 5, and the result shows that this fluctuation frequency is definitely the same as that of the flow rate. These observations exist in the filling phase. In the packing phase, the mean difference between these monitors is the fluctuation amplitude.
Recalling the finite element equation for computing the mechanical field under DIMT Equation (23), we can find the explanation of all of these phenomena. Actually, the entire solution of the equations not only includes the solution of the CIMT part, but also superposes the solution of the “vibration term”. For the solution of the CIMT part, it gives people the main trend of the pressure response curve, whereas the solution of the “vibration term” makes a certain pressure revision on the basis of the CIMT solution. For example, due to the periodic behavior of the flow rate in the dynamic filling phase, the higher flow rate of the DIMT case than the constant one of the CIMT case induces a pressure ascent, while the lower flow rate induces a pressure descent. This is the main reason that the pressure of the DIMT case fluctuates around the pressure of the CIMT case. For the difference of the fluctuation amplitude between these monitors, it is because the pressure value depends on the flow length between the node and the flow front. The larger the flow length is, the larger the fluctuation amplitude presents.
5. Product Warpage Optimization Design of Different Structure Types
5.1. Optimization Model
For thin-walled products, warpage is one of the most important molding defects. In the following sections, the purpose is to minimize the warpage of thin-walled products and argue the actual effect of DIMT compared with CIMT. The warpage optimization model is proposed as:
where
(
x),
(
x), and
(
x) represent the deformation in the
x,
y, and
z direction of node
i.
Nnode and
Nvar represent the node number and the variable number, respectively.
In order to compare the improvement level of warpage, the process parameter optimizations for CIMT and DIMT are both executed. For CIMT, melt temperature x1 (Tmelt), mold temperature x2 (Tmold), filling time x3 (tfill), packing time x4 (tpack), initial packing fraction x5 (Pstart), and trend coefficient x6 (δ) are selected as the optimization parameters. For DIMT, except the previous six parameters, filling amplitude x7 (ε), filling frequency x8 (ffill), filling phase x9 (φfill), packing amplitude x10 (θ), packing frequency x11 (fpack), and packing phase x12 (φpack) are chosen, additionally.
The optimization objective in Equation (40) is a “black” function with respect to
x. In this paper, a Kriging surrogate model [
26,
27] is taken into consideration to evaluate the objective approximately, since it can give us more prediction information, and a maximum experiment improvement (MEI) function-based sequential optimization method [
28] is employed to solve Equation (40). The detailed contents of the optimization method and procedure can be found in
Appendix A. Then, the equal optimization issue changes to be:
with the convergence criterion:
where
and
are the maximum and minimum response values in the sample set of the
Ith searching for Equation (41), respectively. Herein,
oc is 0.001 in this paper.
5.2. Structure and Parameters
In this section, three types of warpage optimization issues are studied to argue the improvement capability of DIMT compared with CIMT.
(a) Small Size Product with No Holes
The first product model is a polycarbonate (PC, the type is X-1) light cover with no holes in it. The geometry and the finite element model of the product are sketched as shown in
Figure 12. The basic properties of the molding material are listed in in detail
Table 4. The value ranges of 12 process parameters are tabulated as
Table 5. For the optimization of CIMT, the design variables involved are the first six variables in
Table 5.
(b) Small Size Product with Multiple Holes
The second product model is a phone cover used in our previous work [
20], which is made by the blend of PC and acrylonitrile butadiene styrene (ABS) copolymers. Herein, it is only combined as a comparison with other issues. In this case, there are multiple holes in the product compared with issue (a).
Figure 13 is the geometry and the finite element model of the product. The basic properties of the molding material are listed in
Table 6. The value ranges of 12 process parameters are tabulated as
Table 7, in which the design variables involved in CIMT optimization issue are the first six variables.
(c) Large, Narrow Product
The third product is a large, narrow PC light cover. The geometrical structure information and finite element model of this PC cover is sketched in
Figure 14. The basic properties information of the molding material is listed in
Table 8. The value ranges of 12 process parameters are tabulated as
Table 9. For the CIMT case, the design variables take the first six variables in
Table 9.
5.3. Results and Discussions
By using MEI-based optimization method in
Appendix A, all of the issues of warpage optimization are solved. The optimization results of three CIMT cases are listed in
Table 10,
Table 12, and
Table 14, respectively. The optimization results of three DIMT cases are listed in
Table 11,
Table 13, and
Table 15, respectively. For issue (a), the warpage of the PC light cover under DIMT (0.168 mm) is smaller than that of the CIMT case (0.197 mm), which is reduced by 14.72%. For issue (b), the warpage of the PC/ABS phone cover under DIMT (0.094 mm) is smaller than that of the CIMT case (0.221 mm), which is reduced by 57.47%. However, for issue (c), the warpage of the PC car light cover under DIMT (0.375 mm) is almost the same as that of the CIMT case (0.388 mm). This shows that DIMT do not improve the molding quality better than CIMT for this large narrow product. The warpage contour figures of all of the results are organized in
Figure 15,
Figure 16 and
Figure 17, respectively.
For these three issues, the results show that DIMT can greatly reduce the warpage of the products, even though the structure and molding material are different. For issues (a) and (c), the dynamic parameters are all “efficient”. However, the
fpack and
φpack are both “inefficient” in issue (b) since the packing amplitude is 0, and the molding process yields into the CIMT case during packing phase. This is to say that, in issue (b), only dynamic filling parameters worked. The discussion of this optimization result is not the main work here, and has been stated in our previous work [
20].
From the warpage contours, it shows that, in issues (a) and (b), the distribution area of the smaller warpage region under DIMT are larger than that of the CIMT cases, and the distribution area of the larger warpage region under DIMT is smaller than that of the CIMT cases. Therefore, DIMT can not only reduce warpage much more than CIMT, but also improve the entire molding quality for the small size products. Unfortunately, it does not happen for issue (c), which makes the molding quality much worse than CIMT. We should recall the pressure response observation of the cavity in
Section 3. From the analysis in
Section 3, we should realize that the pressure volatility relies on the flow length. If the geometry size is large, at the location far away from the gate, the pressure will not be affected by the vibration. This causes a great difference between the region around the gate and the region away from the gate. It may not have a good mechanical history for improving the molding quality. Therefore, it fails.
Actually, the in-cavity residual stress is the final result of the thermal and mechanical evolution before ejection. The uneven distribution of the in-cavity residual stress is the only reason for the warpage after ejection. DIMT takes an external force into the flow history. For small size product, it can disturb the flow behavior and may eliminate the magnitude and distribution difference of the in-cavity residual stress, which results in a lower warpage level compared with CIMT. Whereas, for large-sized, narrow products, this contribution is limited around the gate, and may not be good for improving the molding quality, even making it worse. Therefore, using DIMT should be seriously considered to produce the polymer product.
6. Conclusions
In this paper, the numerical DIMT is proposed by using the finite element method. More detailed work has been performed for observing and analyzing the mechanical evolution behavior of DIMT. The effect of dynamic parameters on the transient pressure of the cavity and the mean pressure of the gate are analyzed for realizing the mechanical behavior of DIMT. It is also fundamental for explaining the mechanism in the following warpage optimization work.
Three kinds of optimization issues with different molding materials are constructed in this paper. By using the EGO sequential optimization method, warpage values of all of the products are efficiently reduced. The further analysis and discussions reveal that DIMT is not available for all types of structures, especially for the large-scale, narrow structures. The main reason for this is that the thermal and mechanical evolution history can determine the in-cavity residual stress in the product before being ejected out of the mold. For small-sized products, dynamic process parameters make the key contribution to disturb the in-cavity residual stress distribution, which can finally increase the region of lower warpage and reduce the warpage difference in the product. Whereas, for large and narrow products, the dynamic contribution concentrates around the gate. It makes a great mechanical difference for this kind of structure, and causes a significant difference in warpage in the product. Therefore, as a recommendation, for this type of product, DIMT should be seriously considered. Overall, the content of this study focuses on revealing the mechanical mechanism of DIMT, and may be a guide for the actual manufacturing with DIMT.