Glass fiber-filled thermoplastics. II. Cavity filling and fiber orientation in injection molding

Mold filling of a rectangular cavity of three different thick nesses fed from a reservoir is studied for unfilled and glass fiber-filled polypropylene and polystyrene. The shapes of flow fronts studied by short-shots are affected predominantly by the thickness of the cavity with other parameters playing a less important role. Pressure drop versus volumetric flow rate inside the thinnest cavity is studied experimentally and predictions are made from a computer simulation of mold filling. The orientation of fibers in the cavity is examined using a reflect-type microscope and the orientation is found to depend on cavity thickness, melt temperature, fiber content, and to a lesser extent, on volumetric flow rate. In the thinnest cavity, where the flow is quasi-unidirectional, the fibers remain in the plane of flow oriented either along the flow direction or perpendicular to it, except in the region near the flow front, where they follow a “fountain” flow behavior.     

Numerical simulation of three‐dimensional fiber orientation in injection molding including fountain flow effect

It is essential to predict the nature of flow field inside mold and flow‐induced variation of fiber orientation for effective design of short fiber reinforced plastic parts. In this investigation, numerical simulations of flow field and three‐dimensional fiber orientation were carried out in special consideration of fountain flow effect. Fiber orientation distribution was described using the second‐order orientation tensor. Fiber interaction was modeled using the interaction coefficient C_I. Three closure approximations, hybrid, modified hybrid, and closure equation for C_I=0, were selected for determination of the fiber orientation. The fiber orientation routine was incorporated into a previously developed program of injection mold filling (CAMPmold), which was based on the fixed‐grid finite element/finite difference method assuming the Hele‐Shaw flow. For consideration of the fountain flow effect, simplified deformation behavior of fountain flow was employed to obtain the initial condition for fiber orientation in the flow front region. Comparisons with experimental results available in the literature were made for film‐gated strip and centergated disk cavities. It was found that the orientation components near the wall were were accurately predicted by considering the fountain flow effect. Test simulations were also carried out for the filling analysis of a practical part, and it was shown that the currently developed numerical algorithm can be effectively used for the prediction of fiber orientation distribution in complex parts.     

Fiber orientation and mechanical properties of short-fiber-reinforced injection-molded composites: Simulated and experimental results

A numerical simulation is presented that combines the flow simulation during injection molding with an efficient algorithm for predicting the orientation of short fibers in thin composite parts. Fiber-orientation state is represented in terms of a second-order orientation tensor. Fiber-fiber interactions are modeled by means of an isotropic rotary diffusion. The simulation predicts flow-aligned fiber orientation (shell region)near the surface with transversely aligned (core region) fibers in the vicinity of the mid-plane. The effects of part thickness and injection speed on fiber orientation are analyzed. Experimental measurements of fiber orientation in plaque-shaped parts for three different combinations of cavity thickness and injection speed are reported. It is found that gapwise-converging flow due to the growing layer of solidified polymer near the walls tends to flow-align the fibers near the entrance, whereas near the melt front, gapwise-diverging flow due to the diminishing solid layer tends to lign the fibers transverse to the flow. The effect of this gapwise-converging-diverging flow is found to be especially significant for thin parts molded at slower injection speeds, which have a proportionately thicker layer of solidified polymer during the filling process. If the fiber orientation is known, predictions of the anisotropic tensile moduli and thermal-expansion coefficients of the composite are obtained by using the equations for unidirectional composites and taking an orientation average. These predictions are found to agree reasonably well with corresponding experimental measurements.     

Prediction of fiber orientation in a rotating compressing and expanding mold

In the rotating/compressing/expanding mold (RCEM), one mold wall can expand, compress, and rotate during injection molding, thus offering opportunities to control the thermomechanical history of a polymer and its microstructure. A computer simulation of flow and fiber orientation in RCEM was developed. The predictive model extends the generalized Hele‐Shaw formulation to account for compression/expansion and rotation of the mold wall, and uses the Folgar–Tucker model for fiber orientation predictions. A 20% GF polypropylene was molded under various molding conditions. The predicted fiber orientation distributions were compared with experiments. The model compares favorably with experiments, provided that the fiber orientation equation is modified by a strain‐reduction factor that slows the transient development of fiber alignment. The effect of fountain flow on orientation must also be included to correctly predict fiber orientation near the mold walls, mainly for the case of stationary and linear motions of the mold surface. Compression or expansion of the mold has only a small effect on fiber orientation, but rotation of the mold dramatically changes the orientation, causing fibers to align in the tangential direction across the entire thickness of the molding. This rotation action perturbs the fountain flow and becomes the dominant factor affecting fiber alignment across the entire cavity thickness. POLYM. ENG. SCI., 2008. © 2008 Society of Plastics Engineers.     

短纤维-弹性体复合材料的注射充模I．数学物理模型的建立

研究了短纤维－弹性体复合材料注射充模过程中,流动前锋和上的喷泉效应,固化层及纤维用量与纤维之间相互作用等的影响,建立了薄壁拉伸试样型腔中短纤维－弹性体复合材料注射充模过程及其纤维取向分布的数学物理模型。     

Molecular orientation distributions during injection molding of liquid crystalline polymers: Ex situ investigation of partially filled moldings

The development of molecular orientation in thermotropic liquid crystalline polymers (TLCPs) during injection molding has been investigated using two‐dimensional wide‐angle X‐ray scattering coordinated with numerical computations employing the Larson–Doi polydomain model. Orientation distributions were measured in “short shot” moldings to characterize structural evolution prior to completion of mold filling, in both thin and thick rectangular plaques. Distinct orientation patterns are observed near the filling front. In particular, strong extension at the melt front results in nearly transverse molecular alignment. Far away from the flow front shear competes with extension to produce complex spatial distributions of orientation. The relative influence of shear is stronger in the thin plaque, producing orientation along the filling direction. Exploiting an analogy between the Larson–Doi model and a fiber orientation model, we test the ability of process simulation tools to predict TLCP orientation distributions during molding. Substantial discrepancies between model predictions and experimental measurements are found near the flow front in partially filled short shots, attributed to the limits of the Hele–Shaw approximation used in the computations. Much of the flow front effect is however “washed out” by subsequent shear flow as mold filling progresses, leading to improved agreement between experiment and corresponding numerical predictions. POLYM. ENG. SCI.,, 2011. © 2011 Society of Plastics Engineers     

Formation mechanism of high-pressure water penetration induced fiber orientation in overflow water-assisted injection molded short glass fiber-reinforced polypropylene

Recently, there has been growing interest in water-assisted injection molding (WAIM) not only for its advantages over gas-assisted molding (GAIM) and conventional injection molding (CIM), but also for its great potential advantages in industrial applications. To understand the formation mechanism of water penetration induced fiber orientation in overflow water-assisted injection molding (OWAIM) parts of short glass fiber-reinforced polypropylene (SGF/PP), in this work, the external fields and water penetration process within the mold cavity were investigated by experiments and numerical simulations. The results showed that the difference of fiber orientation distribution in thickness direction between WAIM moldings and CIM moldings was mainly ascribed to the great external fields generated by water penetration. Besides, fiber orientation depended on the position both across the part thickness and along the flow direction. Especially in the radial direction, fiber orientation varied considerably. The results also showed that the melt temperature is the principal parameter affecting the fiber orientation along the flow direction, and a higher melt temperature significantly facilitated more fibers to be oriented along the flow direction, which is quite different from the results as previously reported in short-shot water-assisted injection molding (SSWAIM). A higher water pressure, shorter water injection delay time, and higher melt temperature significantly induced more fibers to be orderly oriented in OWAIM moldings, which may improve their mechanical performances and broaden their application scope.     

Effects of molding conditions on the electromagnetic interference performance of conductive ABS parts

Polymers filled with conducting fibers to prevent electromagnetic interference (EMI) performance have recently received great attention due to the requirements of 3C (computer, communication, and consumer electronics) products. In the present article, the effect of fiber content and processing parameters, including melt temperature, mold temperature, and injection velocity, on the electromagnetic interference shielding effectiveness (SE) in injection molded ABS polymer composites filled with conductive stainless steel fiber (SSF) was investigated. The influence of fiber orientation and distribution resulting from fiber content and molding conditions on EMI performance was also examined. It was found from measured results that fiber content plays a significant role in influencing part EMI SE performance. SE value can reach the highest values of approximately 40 dB and 60 dB at 1000 MHz frequency for fiber content 7 wt % and 14 wt %, respectively, under the best choice of molding conditions. Higher melt and mold temperature would increase shielding effectiveness due to a more uniform and random fiber orientation. However, higher injection velocity leading to highly‐orientated and less uniform distribution of fiber reduces shielding effectiveness. Among all molding parameters, melt temperature affects SE performance most significantly. Its influence slightly decreases as fiber content increases. Injection speed plays a secondary importance in affecting SE values, and its influence increases as fiber content increases. Upon examination of fiber distribution via optical microscope and subsequent image analysis, it was found that the fiber becomes more densely and random distributed toward the last melt‐filled region, whereas fiber exhibits less concentration around the middle way of the flow path. This can be attributed to the combined effects of fountain flow, frozen layer thickness, and gapwise melt front velocity. The results indicate that molding conditions, instead of fiber content alone, are very important on the SE performance for injection molded SSF filled ABS composites. © 2005 Wiley Periodicals, Inc. J Appl Polym Sci 98: 1072–1080, 2005     

注塑成型熔接线截面熔接过程数值模拟

在多浇口和带嵌件注塑制品的成型过程中必然存在熔体的熔接过程,从而形成熔接线。熔接线沿厚度方向的熔接过程是影响该区域的力学强度以及纤维取向等制品性能的重要因素。本文采用有限元法针对注塑制品的典型截面建立数学模型,采用T6P3单元(速度二次插值,压力线性插值),数值模拟了注塑制品熔接线的截面熔接过程。通过等厚度截面和非等厚度截面两个算例,给出了两股熔体熔接过程中的截面速度场和压力场分布。讨论了熔接线区域的壁厚均匀程度对熔接过程的影响。该计算结果可以为制品力学性能以及纤维取向等数值模拟提供数据支持。     

Analytical solutions for fiber orientation in 2-D flows of dilute suspensions

Analytical solutions for fiber orientation in dilute solutions are presented for both 2-D planar and axisymmetric flow. The viscosity is assumed to be independent of position along the flow direction, but varies through the thickness of the mold cavity. For the particular case of a laminate of two fluids, an analytical solution is derived. Based upon the assumption that the fluid adjacent to the mold wall exhibits reduced viscosity due to non-isothermal considerations, the fluid kinematics simplify, and an analytical solution to Jeffery's orientation equation is derived. Since the fiber orientation analysis is based on Jeffery's equation, it is valid only for dilute suspensions, i.e. for fiber volume fraction < 1/(fiber aspect ratio)². Qualitatively, the analytical results obtained provide superior correlation with experimental results, indicating that non-isothermal fluid flow may play an important role in the development of the fiber orientation distribution in the molded component.     

Modeling fountain flow and filling front shape in reaction injection molding

A numerical model of the reaction injection molding process was developed to test front shape and flow approximations employed in previous models. The model was two-dimensional and simulated the flow, reaction, and heat transfer in the typically long axial dimension and the typically small thickness dimension of a mold. The filling front shape and the velocity profiles in the filling fluid were determined by numerical solution of the momentum equation with the appropriate stress boundary conditions using the method of Patankar (1980). The predicted temperature and conversion results agreed with calculations assuming that the front was flat perpendicular to the flow and that a parabolic velocity profile existed behind the fountain flow region at the front. Thus, simple assumptions about front shape and velocity in the thin dimension of a reaction injection mold can be employed without significant loss of accuracy in modeling reaction injection molding.     

Mechanical properties of glass fiber reinforced polypropylene injection molded with a rotation mold

A special mold (Rotation, Compression, and Expansion Mold) was used to impose a controlled shear action during injection molding of short glass fiber reinforced polypropylene discs. This was achieved by superimposing an external rotation to the pressure‐driven advancing flow front during the mold filling stage. Central gated discs were molded with different cavity rotation velocities, inducing distinct levels of fiber orientation through the thickness. The mechanical behavior of the moldings was assessed, in tensile and flexural modes on specimens cut at different locations along the flow path. Complete discs were also tested in four‐point flexural and in impact tests. The respective results are analyzed and discussed in terms of relationships between the developed fiber orientation level and the mechanical properties. The experimental results confirm that mechanical properties of the moldings depend strongly on fiber orientation and can thus be tailored by the imposed rotation during molding. POLYM. ENG. SCI. 46:1598–1607, 2006. © 2006 Society of Plastics Engineers.     

凸截面注射成型全三维充填和纤维取向模拟研究

选用15 %玻璃纤维增强聚酰胺6（PA6）复合材料,根据Folgar-Tucker理论模型,研究了凸形截面零件全三维注射成型充填流动行为和纤维取向分布规律。结果表明,全三维模拟技术可以有效反映注射成型充填的流动形态和喷泉效应;凸形截面表层纤维沿剪切方向取向分布,且制品大端部分的纤维取向程度比小端部分差;3个坐标方向上的纤维取向分布非常复杂,沿x轴负方向（即充填方向）,随充填长度的增加,纤维取向程度先增大后减小;沿y轴方向,随着距中心距离的增大,凸形截面小端部分的纤维取向程度由内向外逐渐增大,而大端部分的纤维取向程度先增大后减少;沿z轴方向,凸形截面小端部分的纤维取向由内向外呈增大趋势,而大端部分的纤维取向程度却逐渐减小。     

Finite element analysis for the flow characteristics along the thickness direction in injection molding

A numerical simulation of the filling stage along the thickness direction is proposed by combining the free-surface boundary condition with the relevant governing equations. The mathematical model is based on the equations of continuity, momentum, and energy along with the inelastic power-law model and relevant boundary conditions. Because of the significant implications for microstructure development in the product, the fountain effect at the advancing free surface is explicitly taken into consideration in the simulation. The model yields data on free-surface shape and velocity, pressure, temperature, and shear stress distributions within the mold cavity. The rearrangement of the velocity and temperature profiles in the vicinity of the melt front is considered in detail.     

Boundary element study on particle orientation caused by the fountain flow in injection molding

A boundary element model is used for simulating the fountain flow in injection molding. The particle orientation caused by the fountain effect is illustrated with a tracer-in-domain technique which is similar to the Marker-And-Cell technique. The efficiency and accuracy of the boundary element model are shown through this simulation, in determination of the free surface shape, and the particle orientation patterns. Clear and correct information is obtained to interpret the complicated effects upon formation of the orientation, especially regarding the “V” shape pattern in the skin areas. The simulation is in agreement with both Schmidt's experiment and some finite element analysis.     

Numerical analysis of injection molding of glass fiber reinforced thermoplastics. Part 2: Fiber orientation

Numerical predictions of fiber orientation during injection molding of fiber-filled thermoplastics are compared to measurements. The numerical work successfully describes the flow of fiber-filled plastic during injection molding, using finite-difference solutions for the transport equations and marker particles to track the flow front. The flow is modeled as a 2-D, non-isothermal, free-surface flow with a new viscosity model dependent upon temperature, pressure, and fiber concentration. The fiber orientation is based upon solution to the Fokker-Planke equation. The comparison demonstrates fair agreement between predicted fiber orientation and experimental results for slow and fast injection speeds. For the slow speed case at 10 and 20 wt% fibers, the numerical and experimental works show that the fibers are more random at the flow front than at the centerline, and that the fibers become more aligned as they flow from the gate to the midstream region. At fast injection speeds, the agreement between the numerical and experimental works is not as good as at slow injection speed. Possible explanations for the discrepancies are that the flow is assumed to be simple shear when injection molding is known to be a pressure-driven flow, the fibers have an initial orientation for the runner rather than the assumed random orientation, the fibers that were displayed from the camera were more oriented just behind the flow front (owing to the fast injection speed), and the orientation requires more than a 2-D video image to represent a 3-D fiber orientation.     

Computer simulation and analysis of fountain flow in filling process of injection molding

The viscous flow in the filling stage of injection molding can be described in terms of an one-dimensional fully developed main flow and a complex two-dimensional flow near the advancing front, which is often termed the fountain flow. The transport characteristics in the front region of the mold flow are of increasing importance in injection process of composite materials such as resin injection molding (RIM). By using of finite element method, the simulation of non-isothermal viscous flow between two isothermal parallel plates with the generalized newtonian fluid is presented in detail. The un-folding of the fluid particles towards the mold wall directly affects transport characteristics such as the distribution of temperature, the orientation and the concentration of molecule near the front in filling stage.     

Dynamics of semi-flexible and breakable fibers under Poiseuille flow

In the processing of fiber-reinforced polymer composites, especially in injection molding, the fiber's orientation, length, and distribution vary depending on the location of the channel flow and its properties, which affects the performance of final products. To investigate the intricate behavior of fiber suspensions under Poiseuille flow, we used the hybrid simulation approach, multiparticle collision dynamics–molecular dynamics (MPC-MD), which takes hydrodynamic interactions and fiber properties (strength, flexibility) into account. For non-breakable and rodlike fibers, fibers align well along the flow direction while showing more alignment near the wall. As fiber becomes breakable and/or flexible, the length and orientation of fibers strongly depend on their properties. The interesting phenomenon is specifically seen for breakable and semiflexible fibers, where the orientation of the fiber exhibits non-monotonic behavior depending on the flow rate. This complex behavior highlights the importance of comprehending the dynamics of many types of fibers and necessitates research into the best conditions for injection molding.     

Measurement and numerical simulation of three‐dimensional fiber orientation states in injection‐molded short‐fiber‐reinforced plastics

To determine three‐dimensional fiber orientation states in injection‐molded short‐fiber composites, a confocal laser scanning microscope (CLSM) is used. Since the CLSM optically sections the specimen, more than two images of the cross sections on and below the surface of the composite can be obtained. Three‐dimensional fiber orientation states can be determined by using geometric parameters of fiber images obtained from two parallel cross sections. For experiments, carbon‐fiber‐reinforced polystyrene is examined by the CLSM and geometric parameters of fibers on each cross‐sectional plane are measured by an image analysis. In order to describe fiber orientation states compactly, orientation tensors are determined at different positions of the prepared specimen. Three‐dimensional orientation states are obtained without any difficulty by determining the out‐of‐plane angles utilizing fiber images on two parallel planes acquired by the CLSM. Orientation states are different at different positions and show the shell–core structure along the thickness of the specimen. Fiber orientation tensors are predicted by a numerical analysis and the numerically predicted orientation states show good agreement with measured ones. However, some differences are found at the end of cavity. They may result from the fountain flow effects, which are not considered in the numerical analysis. © 2003 Wiley Periodicals, Inc. J Appl Polym Sci 88: 500–509, 2003