Research Papers

Three-Dimensional Viscoelastic Simulation for Injection/Compression Molding Based on Arbitrary Lagrangian Eulerian Description

[+] Author and Article Information
Wei Cao, Shaozhen Hua, Shixun Zhang, Yaming Wang, Haimei Li, Chuntai Liu, Changyu Shen

National Engineering Research Center of
Mold and Die,
Zhengzhou University,
Zhengzhou, Henan 450002, China

Tao Wang

Beijing Institute of Aeronautical Materials,
Beijing 100095, China

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received May 23, 2015; final manuscript received December 27, 2015; published online February 3, 2016. Assoc. Editor: Zdravko Terze.

J. Comput. Nonlinear Dynam 11(5), 051004 (Feb 03, 2016) (8 pages) Paper No: CND-15-1139; doi: 10.1115/1.4032384 History: Received May 23, 2015; Revised December 27, 2015

Different from conventional injection molding (CIM), injection/compression molding (ICM) evolves boundary variation in gapwise direction. In order to describe melt flow characteristics in ICM correctly, a new material derivative based on arbitrary Lagrangian Eulerian (ALE) description was introduced to modify the material derivatives in the governing and constitutive equations. To avoid large amount of calculation and weak stability of integral numerical method, an iterative approach employing twofold iterations was proposed to decouple the interdependence between velocity, stress, and temperature. The initial values of material parameters in constitutive equations were obtained or fitted by rheological experiments. The ICM experiments for an iso-thick and a var-thick rectangular panel were carried out to validate the proposed method and find the special characteristics of ICM. In addition, the photoelastic tests on a quarter of spherical part processed by ICM were conducted to identify the relationship between residual flow-induced stress distributions and flow fields. Both simulations and experiments show that the pressure profile displays a plateau during compression, temperature decreases with time according to exponential law, large flow-induced stress originates in thick transitional region, flow start, and flow end areas, and gravity has significant effect on meltfront for thick part ICM. The good agreement between experiments and simulations indicates that the current method can properly describe the flow characteristics of ICM.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.


Wu, C.-H. , and Chen, W.-S. , 2006, “ Injection Molding and Injection Compression Molding of Three-Beam Grating of DVD Pickup Lens,” Sens. Actuators A, 125(2), pp. 367–375. [CrossRef]
Sortino, M. , Totis, G. , and Kuljanic, E. , 2014, “ Comparison of Injection Molding Technologies for the Production of Micro-Optical Devices,” Procedia Eng., 69(1), pp. 1296–1305. [CrossRef]
Lawal, A. , and Kalyon, D. M. , 2000, “ Compressive Squeeze Flow of Generalized Newtonian Fluids With Apparent Wall Slip,” Int. Polym. Process., 15(1), pp. 63–71. [CrossRef]
Lee, S. J. , Denn, M. M. , Crochet, M. J. , Metzner, A. B. , and Riggins, G. J. , 1984, “ Compressive Flow Between Parallel Disks II. Oscillatory Behavior of Viscoelastic Materials Under a Constant Load,” J. Non-Newtonian Fluid Mech., 14, pp. 301–325. [CrossRef]
Isayev, A. I. , and Azari, A. D. , 1986, “ Viscoelastic Effect in Compression Molding of Elastomers: Shear-Free Squeezing Flow,” Rubber Chem. Technol., 59(5), pp. 868–882.
Isayev, A. I. , Zhang, Y. , and Zook, C. , 1999, Advances in the Flow and Rheology of Non-Newtonian Fluids, Elsevier Science, Amsterdam, The Netherlands.
Kim, I . H. , Park, S. J. , Chung, S. T. , and Kwon, T. H. , 1999, “ Numerical Modeling of Injection/Compression Molding for Center-Gated Disk (Part I: Injection Molding With Viscoelastic Compressible Fluid Model),” Polym. Eng. Sci., 39(10), pp. 1930–1942. [CrossRef]
Kim, N. H. , 2009, “ Injection-Compression and Co-Injection Molding of Amorphous Viscoelastic Simulation and Experiment,” Ph.D. dissertation, University of Akron, Akron, OH.
Wang, J. T. , 1997, CAE and Intelligent Processing of Polymeric Materials, ASME, New York.
Kim, I . H. , Park, S. J. , Chung, S. T. , and Kwon, T. H. , 1999, “ Numerical Modeling of Injection/Compression Molding for Center-Gated Disk (Part II: Effect of Compression Stage),” Polym. Eng. Sci., 39(10), pp. 1943–1951. [CrossRef]
Lee, Y. B. , Kwon, T. H. , and Yoon, K. , 2002, “ Numerical Prediction of Residual Stresses and Birefringence in Injection/Compression Molded Center-Gated Disk (Part II: Effects of Processing Conditions),” Polym. Eng. Sci., 42(11), pp. 2273–2292. [CrossRef]
Chen, S.-C. , Chen, Y.-C. , Peng, H.-S. , and Huang, L.-T. , 2002, “ Simulation of Injection Compression Molding Process (Part 3: Effect of Process Conditions on Part Birefringence),” Adv. Polym. Technol., 21(3), pp. 177–187. [CrossRef]
Chen, S. C. , Chen, Y. C. , and Peng, H. S. , 2000, “ Simulation of Injection-Compression Molding Process (II: Influence of Process Characteristics on Part Shrinkage),” J. Appl. Polym. Sci., 75(13), pp. 1640–1654. [CrossRef]
Hu, S. T. , Chiu, H. S. , Chien, C. C. , Yu, C. K. , and Chang, R. Y. , 2010, “ True 3D Numerical Simulation in Injection Compression Molding (ICM),” SPE ANTEC, Orlando, FL, May 16–20, 2010, SPE ANTEC Technical Paper No. 1, pp. 716–720.
Li, Y. , Zhang, Y. , and Li, D. , 2010, “ Shrinkage Analysis of Injection-Compression Molding for Transparent Plastic Panel by 3D Simulation,” Appl. Mech. Mater., 44–47, pp. 1029–1033. [CrossRef]
Kim, N. H. , and Isayev, A. I. , 2013, “ Birefringence in Injection-Compression Molding of Amorphous Polymers: Simulation and Experiment,” Polym. Eng. Sci., 53(8), pp. 1786–1808. [CrossRef]
Braess, H. , and Wriggers, P. , 2000, “ Arbitrary Lagrangian Eulerian Finite Element Analysis of Free Surface Flow,” Comput. Methods Appl. Mech. Eng., 190(1–2), pp. 95–109. [CrossRef]
Ho, J.-Y. , Park, J. M. , Kang, T. G. , and Park, S. J. , 2012, “ Three-Dimensional Numerical Analysis of Injection-Compression Molding Process,” Polym. Eng. Sci., 52(4), pp. 901–911. [CrossRef]
Dwiyantoro, B. A. , 2013, “ A Numerical Study of an Injection-Compression Molding Process by Using a Moving Grid,” Appl. Mech. Mater., 249–250, pp. 472–476.
Mahmadi, K. , and Aquelet, N. , 2014, “ Delayed Mesh Relaxation for Multi-Material ALE Formulation,” Int. J. Heat Fluid Flow, 46(4), pp. 102–111. [CrossRef]
Zatloukal, M. , 2003, “ Differential Viscoelastic Constitutive Equations for Polymer Melts in Steady Shear and Elongational Flows,” J. Non-Newtonian Fluid Mech., 113(2–3), pp. 209–227. [CrossRef]
Verbeeten, W. M. H. , Peters, G. W. M. , and Baaijens, F. P. T. , 2004, “ Numerical Simulations of the Planar Contraction Flow for a Polyethylene Melt Using the XPP Model,” J. Non-Newtonian Fluid Mech., 117(117), pp. 73–84. [CrossRef]
Chang, H. H. , Hieber, C. A. , and Wang, K. K. , 1991, “ A Unified Simulation of the Filling and Postfilling Stages in Injection Molding (Part I: Formulation),” Polym. Eng. Sci., 31(2), pp. 116–124. [CrossRef]
Fan, Y. , Tanner, R. I. , and Phan-Thien, N. , 1999, “ Galerkin/Least-Square Finite-Element Methods for Steady Viscoelastic Flows,” J. Non-Newtonian Fluid Mech., 84(84), pp. 233–256. [CrossRef]
Guénette, R. , and Fortin, M. , 1995, “ A New Mixed Finite Element Method for Computing Viscoelastic Flows,” J Non-Newtonian Fluid Mech., 60(1), pp. 27–52. [CrossRef]
Silva, L. , Valette, R. , Laure, P. , and Coupez, T. , 2012, “ A New Three-Dimensional Mixed Finite Element for Direct Numerical Simulation of Compressible Viscoelastic Flows With Moving Free Surfaces,” Int. J. Mater. Form., 5(1), pp. 55–72. [CrossRef]
Hirt, C. W. , Amsden, A. A. , and Cook, J. L. , 1997, “ An Arbitrary Lagrangian–Eulerian Computing Method for All Flow Speeds,” J. Comput. Phys., 135(2), pp. 203–216. [CrossRef]
Lai, H. E. , and Wang, P. J. , 2008, “ Study of Process Parameters on Optical Qualities for Injection-Molded Plastic Lenses,” Appl. Opt., 47(12), pp. 2017–2027. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

Diagram of a point transformation in ALE scheme

Grahic Jump Location
Fig. 2

Diagram of polymer melt flow during compression

Grahic Jump Location
Fig. 3

Mixed finite elements used in this study with (a) P2(10) element: four vertices + six edge midpoints, (b) P1(4) element: four vertices in which •—value of u,T and ×—value of p, τ

Grahic Jump Location
Fig. 4

Diagram of the dimension of the studied parts: (a) iso-thick part and (b) var-thick part. O—sensor location and ◀—polymer entrance.

Grahic Jump Location
Fig. 5

Simulated meltfronts of iso-thick panel at times: (a) 3 s, (b) 5 s, (c) 7 s, and (d) 9 s

Grahic Jump Location
Fig. 6

Measured pressures of sensors 2 and 4 change with time during iso-thick part ICM

Grahic Jump Location
Fig. 7

Simulated meltfronts of var-thick panel at times: (a) 3 s, (b) 5 s, (c) 8 s, and (d) 10 s

Grahic Jump Location
Fig. 8

Measured sensors' pressures of iso-thick panel injection molding with injection velocity (for screw): 60 mm/s, filling/packing switch: 30 mm, packing pressure: 80 MPa, and melt temperature: 300 °C

Grahic Jump Location
Fig. 9

Comparison between the measured and simulated sensors' pressures of (a) iso-thick panel and (b) var-thick panel during filling and compressing stages. − (thick and discrete line) represents experimental data and − (thin and smooth line) denotes the simulated data.

Grahic Jump Location
Fig. 10

Comparison between the measured and simulated sensors' temperatures of (a) iso-thick panel and (b) var-thick panel during filling and compressing stages. − (thick and discrete line) represents experimental data and − (thin and smooth line) denotes the simulated data.

Grahic Jump Location
Fig. 11

Dimension of one-quarter of spherical part

Grahic Jump Location
Fig. 12

Comparison between the simulated (a) first normal stress difference and (b) photoelasticpictures



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In