Numerical simulation of blow molding—prediction of parison diameter and thickness distributions in the parison formation process

In our previous study, we calculated the time course of parison length in the parison formation stage, but it could predict only the parison area swell. The next target in our study is to calculate the parison diameter and thickness swell. Annular extrudate swell simulation is necessary for the understanding of various kinds of swelling ratios in blow molding. We have examined three kinds of swells (outer diameter, thickness, and area swells) obtained from simulation results of annular extrudate swell, using the Giesekus model, and have developed a method of predicting parison outer diameter and thickness swell values. The predicted values of parison outer diameters are discussed in comparison with experimental data, and reasonable results are obtained by the proposed method. This prediction method could also be applied to the parison formation process using a parison controller. As a result, it is possible to predict approximately the whole process of parison formation by numerical simulation.     

Modeling of complex parison formation in extrusion blow molding: Effect of medium to large die heads and fuel tank geometry

This work presents the effect of die geometry and die gap opening on the extrudate swell phenomenon, in complex parison formation using the vertical wall distribution system (VWDS) and partial wall distribution system (PWDS). The BlowParison^© software from IMI is used to predict the parison formation for a combined VWDS/PWDS system, accounting for swell, sag, and nonisothermal effects. This software couples a fluid mechanics approach to represent the die flow, with a solid mechanics approach to represent the parison behavior outside the die, and a mathematical swell model to account for the pronounced elongational and shear stresses at high Weissenberg numbers. The emphasis is placed on experimental validation of the predicted parison dimensions using four diverging die geometries and different sets of VWDS/PWDS profiles. The experimental and predicted weight profiles for a dissected fuel tank are also presented. Both experimental and simulation results suggest a strong dependence of extrudate swell to the die geometry in the die land zone. The results also demonstrate the validity of the numerical predictions for part design purposes given the multitude of experimental validations presented in this work. POLYM. ENG. SCI., 2009. Published by the Society of Plastics Engineers     

熔体挤出速度对共挤吹塑型坯离模膨胀影响的数值模拟

基于三维非等温黏弹性熔体多相分层流动有限元数值模拟技术,模拟研究了熔体挤出速度对多层共挤吹塑成型环坯离模膨胀和初始温度场的影响规律,揭示了型坯离模膨胀的产生机理。结果表明,多层共挤吹塑成型环坯离模膨胀是由熔体的二次流动诱发而产生,与熔体流出机头进入自由膨胀段的二次流动强度成正比,而其二次流动强度随着熔体挤出速度的增大而增强,因而导致环坯离模膨胀随着熔体挤出速度的增加而增大;多层共挤吹塑成型熔体的二次流动强度与其第二法向应力差成正比关联关系,这与Debbaut的试验研究结论完全吻合,表明二次流动是由第二法向应力差驱动而产生。     

Effect of die gap width on annular extrudates by the annular extrudate swell simulation in steady-states

We calculated the steady-state annular extrudate swell of polymer melts through flow geometries encountered in processes used to control parison thickness. A streamline-upwinding finite element method with an under-relaxation for the rate of deformation tensor was used. The Giesekus model was employed as the constitutive equation. An operation that widens the die gap is appropriate for the control of parison thickness corresponding to the change of die gap width. However, a control process that decreases the die gap width is not useful, because the parison thickness does not correspond to the die gap width. Furthermore, thickness swells change strikingly with the Weissenberg number. It is difficult to control the parison outer diameter in the case of a converging die, because the change of the outer diameter swell becomes large with increasing Weissenberg number. In the case of a diverging die, the changing value of the outer diameter swell is smaller than that in the case of a converging die.     

塑料挤出吹塑的机理问题

采用不同的方法对挤出吹塑过程的型坯成型、型坯吹胀与制品冷却三阶段的机理问题进行了研究.采用人工神经网络方法预测了受模口温度和挤出流率影响的型坯成型阶段的膨胀.利用建立起来的神经网络模型预示的膨胀与实验结果很吻合,且可在一定范围内,预示不同工艺条件下型坯的直径膨胀和壁厚膨胀,为型坯的直径和壁厚的在线控制提供了理论依据.基于薄膜近似和neo-Hookean本构关系,建立了描述型坯自由吹胀的数学模型,并通过实验方法获得了型坯吹胀的瞬态图象.     

A comprehensive experimental study and numerical modeling of parison formation in extrusion blow molding

Parison dimensions in extrusion blow molding are affected by two phenomena, swell due to stress relaxation and sag drawdown due to gravity. It is well established that the parison swell and sag are strongly dependent on the die geometry and the operating conditions. The availability of a modeling technique ensures a more accurate prediction of the entire blow molding process, as the proper prediction of the parison formation is the input for the remaining process phases. This study considers both the simulated and the experimental effects of the die geometry, the operating conditions, and the resin properties on the parison dimensions using high density polyethylene. Parison programming with a moving mandrel and the flow rate evolution in intermittent extrusion are also considered. The parison dimensions are measured experimentally by using the pinch-off mold technique on two industrial scale machines. The finite element software BlowParison® developed at IMI is used to predict the parison formation, taking into account the swell, sag, and nonisothermal effects. The comparison between the predicted parison/part dimensions and the corresponding experimental data demonstrates the efficiency of numerical tools in the prediction of the final part thickness and weight distributions. POLYM. ENG. SCI., 47:1–13, 2007. © 2006 Society of Plastics Engineers     

A modeling approach to the effect of resin characteristics on parison formation in extrusion blow molding

The most critical stage in the extrusion blow‐molding process is the parison formation, as the dimensions of the blow‐molded part are directly related to the parison dimensions. The swelling due to stress relaxation and sagging due to gravity are strongly influenced by the resin characteristics, die geometry, and operating conditions. These factors significantly affect the parison dimensions. This could lead to a considerable amount of time and cost through trial and error experiments to get the desired parison dimensions based upon variations in the resin characteristics, die geometry, and operating conditions. The availability of a modeling technique ensures a more accurate prediction of the entire blow‐molding process, as the proper prediction of the parison formation is the input for the remaining process phases. This study considers both the simulated and the experimental effects of various high‐density polyethylene resin grades on parison dimensions. The resins were tested using three different sets of die geometries and operating conditions. The target parison length was achieved by adjusting the extrusion time for a preset die gap opening. The finite element software BlowParison® was used to predict the parison formation, taking into account the swell and sag. Good agreements were found between the predicted parison dimensions and the experimental data. POLYM. ENG. SCI., 2009. Published by Society of Plastics Engineers     

Overall numerical simulation of extrusion blow molding process

This paper focuses on the overall numerical simulation of the parison formation and inflation process of extrusion blow molding. The competing effects due to swell and drawdown in the parison formation process were analyzed by a Lagrangian Eulerian (LE) finite element method (FEM) using an automatic remeshing technique. The parison extruded through an annular die was modeled as an axisymmetric unsteady nonisothermal flow with free surfaces and its viscoelastic properties were described by a K‐BKZ integral constitutive equation. An unsteady die‐swell simulation was performed to predict the time course of the extrudate parison shape under the influence of gravity and the parison controller. In addition, an unsteady large deformation analysis of the parison inflation process was also carried out using a three‐dimensional membrane FEM for viscoelastic material. The inflation sequence for the parison molded into a complex‐shaped mold cavity was analyzed. The numerical results were verified using experimental data from each of the sub‐processes. The greatest advantage of the overall simulation is that the variation in the parison dimension caused by the swell and drawdown effect can be incorporated into the inflation analysis, and consequently, the accuracy of the numerical prediction can be enhanced. The overall simulation technique provides a rational means to assist the mold design and the determination of the optimal process conditions.     

Rheological behavior and processability of neoprene and acrylic rubbers

Rheological behavior of neoprene and acrylic rubbers and their carbon blackfilled compounds has been studied using a Monsanto Processability Tester. These systems show pseudoplastic flow behavior. Die swell increases with increase in shear rate up to a limit beyond which it decreases. It also decreases with increase in filler loading and temperature. The maximum recoverable deformation has been calculated by assuming that the viscous response of the rubbers obeys the power law model and elastic behavior is described by Hooke's law and has been correlated with die swell. A linear generalized relationship has been obtained between the die swell and the recoverable deformation. The principal normal stress difference has been found to increase nonlinearly with shear stress. Activation energy of melt flow process increases with increase in shear rate up to a limit, after which it decreases for acrylic rubber and 10 phr filled neoprene rubber.     

Measurement and calculation of parison dimensions and bottle thickness distribution during blow molding

Experimental data are reported regarding the dynamics of the blow molding process, including parison formation, growth, and inflation. These data have been obtained with the aid of high speed cinematography and pinch mold experiments, in conjunction with two commercial blow molding polyethylene resins. It is shown that pinch mold experiments alone do not yield accurate data regarding thickness and diameter swell. Furthermore, the inflation process involves decreasing rates of inflation with time, as a result of the rise in viscosity due to the cooling of the parison during inflation. Mathematical procedures are proposed for a first-order estimation of parison length and swell as a function of time and the inflation behavior after clamping. In the absence of more dependable basic procedures, the proposed treatment is employed to estimate the effective transient swell functions for the parison using experimental data obtained under the specified conditions. The mathematical treatment is extended to determine the thickness distribution of the bottle. Good agreement is obtained between experimental and calculated results.     

挤出吹塑中型坯成型的有限元模拟——型坯尺寸预测

建立了塑料挤出吹塑中平直型机头内非牛顿粘弹性熔体的流动模型，并利用POLYFLOW有限元软件进行求解，预测了不同型坯长度和不同流量时的型坯尺寸分布。     

Predictions of flow behaviors and entrance pressure drop characteristics of a rubber compound in a capillary die using various rheological models

Rubber compounds have highly viscoelastic properties. The viscoelastic behaviors that have been exhibited during die extrusion include die swell and vortices in regions of sudden contraction. In this study, the application of rheological models to the capillary die extrusion process is investigated. Experiments and simulations were conducted using a fluidity tester and finite element analysis, respectively. The velocity distributions, velocity profiles, pressure drops, and vortices at the capillary die entrance were analyzed through computer simulations for various viscoelastic models [i.e., Phan‐Thien and Tanner (PTT), Giesekus, POMPOM, simplified viscoelastic, and generalized Newtonian models]. Different models exhibited different pressure drops and different velocity profiles in the capillary die. Only the full viscoelastic models (PTT, Giesekus, and POMPOM) predicted the vortex at the corner of the reservoir that is the capillary die entrance. However, the simplified viscoelastic and generalized Newtonian models did not predict the vortex. All the viscoelastic models studied in this article predicted the die swells in various ways, and these were compared with the experimental results. The PTT and simplified viscoelastic models exhibited good agreement with the experimental results of the die swells. POLYM. ENG. SCI., 54:2441–2448, 2014. © 2013 Society of Plastics Engineers     

发射药药料离模膨胀和流动均匀性模拟及在模具设计中的应用

为了分析挤出成型过程中模具结构参数对七孔硝基胍发射药离模膨胀率及流动均匀性的影响规律,采用计算流体力学方法,对挤出成型过程进行模拟计算,讨论了模具各结构参数重要性的主次关系;对七孔发射药制备模具进行了结构优化,并进行了实验验证。结果表明,模具收缩角对膨胀率和药料出口速度均匀性的影响最大,压缩段高度次之,成型段长度的影响最小。模具优化后流道出口端速度分布均匀性提高36.53%,表明该模拟计算的可靠性与实用性。     

Modeling parison formation in extrusion blow molding by neural networks

A series of experiments were carried out on the parison formation stage in extrusion blow molding of high‐density polyethylene (HDPE) under different die temperature, extrusion flow rate, and parison length. The drop time of parison when it reached a given length and its swells, including the diameter, thickness, and area swells, were determined by analyzing its video images. Two back‐propagation (BP) artificial neural network models, one for predicting the length evolution of parison with its drop time, the other predicting the swells along the parison, were constructed based on the experimental data. Some modifications to the original BP algorithm were carried out to speed it up. The comparison of the predicted parison swells using the trained BP network models with the experimentally determined ones showed quite a good agreement between the two. The sum of squared error for the predictions is within 0.001. The prediction of the parison diameter and thickness distributions can be made online at any parison length or any parison drop time within a given range using the trained models. The predicted parison swells were analyzed. © 2005 Wiley Periodicals, Inc. J Appl Polym Sci 96: 2230–2239, 2005     

Die swell from capillary die and slit die: A theoretical study

Assuming polymer fluids obey the CEF equation, equations concerning die swell from capillary die and slit die were derived. The die swell of polymer increased with increasing shear rate and recoverable shear. Our theoretical predictions conformed well with the experimental data of die swell for PP and PS.     

Integrated numerical modeling of the blow molding process

The numerical modeling of the extrusion blow molding of a fuel tank is considered in this work. The integrated process phases are consecutively simulated, namely, parison formation, clamping, and inflation, as well as part solidification, part deformation (warpage), and the buildup of residual stresses. The parison formation is modeled with an integral type viscoelastic constitutive equation for the sag behavior and a semi-empirical equation for the swell behavior. A nonisothermal viscoelastic formulation is employed for the clamping and inflation simulation, since parison cooling during extrusion strongly affects the inflation behavior. Once the parison is inflated, it solidifies while in the mold and after part ejection. Warpage and residual stress development of the part are modeled with a linear viscoelastic solid model. Numerical predictions are compared with experimental results obtained on an industrial scale blow molding machine. Good agreement is observed. A process optimization based on a desired objective function, such as uniform part thickness distribution and/or minimal part weight, is performed. The integrated clamping, inflation, and cooling stages of the process are considered. The optimization is done by the systematic manipulation of the parison thickness distribution. Iterations are performed employing a gradient based updating scheme for the parison thickness programming, until the desired objective of uniform part thickness is obtained.     

Noncontact measurement of parison thickness profiles affected by swell and sag in continuous extrusion blow molding

Optimization of final part thickness distributions is crucial in the extrusion blow molding process in order to minimize resin usage. Prediction of part thickness distributions from basic process and material parameters would be ideal. However, attempts to do so have been unsuccessful, largely because of the inability to predict parison thickness profiles. One must therefore resort to measurement of the parison thickness profile and estimation of the final part thickness distribution by computational methods. This paper describes a new technique for the noncontact estimation of parison thickness profiles in continuous extrusion blow molding. The method accounts for sag and requires no previous knowledge of rheological data. It can be employed on-line for the purposes of process monitoring and control. The approach is based on the measurement of the parison length evolution with time during extrusion, the parison diameter profile, the flow rate, and the melt temperature gradient along the length of the parison. These parameters are utilized in conjunction with a theoretical approach that describes the extrusion of a parison under the effects of swell, sag, and extrusion into ambient conditions. Results are presented for three resins of various molecular weight distributions. The degree of sag is minimal at the top and bottom of the parison, and reaches a maximum near the center of the parison. Results are also presented to demonstrate the versatility of the method under other process conditions, such as varying flow rate, die temperature, and die gap.     

Swell and distortions of high-density polyethylene extruded through capillary dies

The effect of varying the die entrance angle and the die length on extrudate swell and on the onset of extrudate distortion in capillary extrusion has been studied. Using theory from the literature, we have analyzed the contribution to the total pressure drop from the elongational and shear deformation in the entrance region, and from the capillary pressure drop in the land region of the die. From the contribution of the elongational deformation, we obtained an estimate for the elongational viscosity of the polymer. The same analysis was used to study the influence of the die geometry on the stick-slip instability. It is found that the elongational component at the inlet region mainly influences the extrudate distortions. The onset of the stick-slip instability occurs within 10% at a wall stress τ_w of 0.3MPa, where τ_w is calculated from expressions assuming fully developed flow. The variation around this average value is systematic with changes in die geometry, and the observed variations are probably due to the non-homogeneous pressure field in the die. We also propose a model for predicting extrudate swell. Input to the model are material parameters obtainable from oscillatoric measurements of the loss and storage modulus and residence times calculated from the geometry of the die. The swell model includes a fitting parameter that sets the overall scale of the swell.     

An experimental investigation of capillary extrudate swell in relation to parison swell behavior in blow molding

An experimental study was carried out to study and characterize the capillary extrudate swell and parison swell behavior in extrusion blow molding of two commercial blow molding grade high density polyethylene resins. The capillary extrudate swell behavior of these resins were determined employing a capillary rheometer and a special thermostatting chamber. Parison swell behavior was determined using an Impco A13-R12 reciprocating screw blow molding machine in conjunction with cinematography and pinch-off. The experimental conditions under which capillary extrudate and parison swell data can be related are elucidated. Excellent agreement is found between the area swell values determined on the basis of capillary and parison swell experiments.     

Model prediction of non-symmetric normal stresses under oscillatory squeeze flow

The non-symmetric responses of normal stresses in oscillatory squeeze flow have been investigated with model calculations. The simplest and most widely used constitutive equations were employed to predict the non-symmetric normal stresses, which is a distinctive feature of oscillatory squeeze flow. The model prediction was compared with experimental data of polymer solution in terms of stress shape, Lissajous plot, stress decomposition, and Fourier transformation. The upper-convected Maxwell, Giesekus, and exponential Phan-Thien Tanner models predicted the nonsymmetric characteristics of normal stresses under oscillatory squeeze flow. The predictions showed fairly good agreement with experimental data. However, the upper-convected Maxwell model showed unrealistic result in the Lissajous plot of [stress vs. strain] and [stress vs. strain rate]. From stress decomposition, it could be confirmed that the non-symmetric nature arises from the elastic contribution of the normal stress, which was verified in both experiment and model calculation. This study is expected to provide useful insights for further understanding of the nonlinear and non-symmetric characteristics of oscillatory squeeze flow.