Coupled analysis of injection molding filling and fiber orientation, including in‐plane velocity gradient effect

On injection molding of short fiber reinforced plastics, fiber orientation during mold filling is determined by the flow field and the interactions between the fibers. The flow field is, in turn, affected by the orientation of fibers. The Dinh and Armstrong rheological equation of state for semiconcentrated fiber suspensions was incorporated into the coupled analysis of mold filling flow and fiber orientation. The viscous shear stress and extra shear stress due to fibers dominate the momentum balance in the coupled Hele-Shaw flow approximation, but the extra in-plane stretching stress terms could be of the same order as those shear stress terms, for large in-plane stretching of suspensions of large particle number. Therefore, a new pressure equation, governing the mold filling process, was derived, including the stresses due to the in-plane velocity gradients. The mold filling simulation was then performed by solving the new pressure equation and the energy equation via a finite element/finite difference method, as well as evolution equations for the second-order orientation tensor via the fourth-order Runge-Kutta method. The effects of stresses due to the in-plane velocity gradient on pressure, velocity, and fiber orientation fields were investigated in the center-gated radial diverging flow in the cases of both an isothermal Newtonian fluid matrix and a nonisothermal polymeric matrix. In particular, the in-plane velocity gradient effect on the fiber orientation was found to be significant near the gate, and more notably for the case of a nonisothermal polymer matrix.