The correlation between the initial time value of transient iterative parameters and the blowing pressure in the numerical simulation process of blowing balloon is investigated by POLYFLOW. The results show that: 1) As the blow molding pressure increases, the boundary value of the iterative time step decreases rapidly at first and then slowly. At the end of the first step of iterative calculation for each boundary value, the balloon parison is in the mold core cavity. 2) If the initial time value of transient iterative parameters is smaller than the boundary value of the iterative time step, the balloon parison is still in the mold core cavity at the end of the first iteration. However, if the iterative calculation continues, the calculation process may be interrupted when the time step is smaller than the initial time value of the transient iterative parameters, which makes the blow molding simulation of balloon unable to continue. 3) It is suggested that the initial time value of transient iterative parameters is one order of magnitude smaller than the boundary value of the iterative time step to complete smoothly the simulation of blow molding balloon.
To form the balloon, which is an important part of the balloon dilating catheter, various numerical methods are employed. Fu et al. used the MSC. Marc to conduct isothermal numerical simulation of the primary and secondary stretching processes of single-cavity balloon made of Polyamide 12, and the influence of stretching velocity on balloon wall thickness is analyzed [
POLYFLOW is widely used in the design of extrusion [
In this work, POLYFLOW is employed to simulate the molding process of single balloon. The correlations between the initial value of transient time step and the pressure of balloon blow molding are studied, which is expected to provide useful suggestion for the choice of the initial time step.
In order to define boundary conditions for contact detection, a distinction must be established between regions where contact has occurred and regions of the free surface where contact has not yet been detected [
A normal force can be applied on an internal free surface to take the inflation into account. The displacement of the free surface can be described by the following dynamic and kinematic conditions prior to contact. The normal force f must be equal to zero, or more generally, be prescribed a given value g .
f = g (1)
where g is nonzero when an external force is applied on the free surface.
In the case of a steady-state problem, the normal velocity must be equal to zero because no mass crosses the free surface.
v ⋅ n = 0 (2)
where v is the velocity vector and n is the normal vector to the free surface.
Once contact occurs, it is impossible for the free surface to displace in the direction normal to the mold surface, and the velocity of the fluid becomes identical to that of the mold. At the same time, the force condition that was previously prescribed is replaced by an essential condition on the velocity in the tangential direction.
Once a solid domain is slightly penetrated by a node, a penalty term can be used to modify the discrete momentum equation of the node. In order to model the physics accurately, the penalty term can be different in the normal and tangential directions.
In the normal direction, the condition that the fluid velocity must be equal to the wall velocity is enforced with the penalty coefficient k:
f n = − k ( v − v w a l l ) n (3)
The default value of k in POLYDATA is 109, which is appropriate for most cases. When the contact detection is modeled, the time-marching parameters must be modified.
In the tangential direction, the velocity is prescribed through a linear law of the “slip” type, given by
f s = − F s l i p ( v s − v s , w a l l ) (4)
where F s l i p is the slip coefficient. The default value of F s l i p in POLYDATA is 109, which usually corresponds to a condition of zero relative velocity in the tangential direction.
The geometric model of balloon blow molding simulation is shown in
The geometric model in
Calculation parameters are set in POLYDATA when the mesh file is imported. The boundary conditions are as follows. The mold core is a solid and the balloon parison is fluid. The inner wall of the parison is free boundary, which bears blow molding pressure. The outer wall of parison and the inner wall of mold core are the contact boundaries. The normal velocity and the tangential force on the end faces of parison are 0.
Under a given blow molding pressure Pb, the boundary value of the initial time value of transient iterative parameters ti is studied by judging the location of parison at the end of the first iterative time step. Different boundary values of the initial time of transient iterative parameters can be obtained by changing Pb.
When ti is reasonably set, the blow molding process of a balloon is carried out by POLYFLOW, then the normal blow molding result of a balloon can be obtained, as shown in
If ti is set beyond a reasonable range, abnormal molding phenomenon occurs. E.g. when ti is equal to 2.63e-5s for Pb = 1.0 MPa, the balloon parison flies out of the mold core cavity at the end of the first iterative calculation, as shown in
When ti is equal to 2.62e-5s for Pb = 1.0 MPa, the balloon parison is still in the mold core cavity at the end of the first iterative calculation, as shown in
When other technological parameters are fixed, under a given blow molding pressure of balloon, the outward expanding velocity of balloon parison is a constant. According to the linear law, it can be known that the product of the expanding velocity of balloon parison and ti is the expanding distance of balloon parison at the end of the first iteration.
When ti is greater than the boundary value, since the outward expansion distance of balloon parison is greater than the radial size of the mold core cavity at the end of the first iteration, the balloon parison can fly out of the mold core cavity (see
When ti is equal to t0, since the outward expansion distance of balloon parison is smaller than the radial size of the mold core cavity at the end of the first iteration, the balloon parison is still in the mold core cavity (see
However, if the iterative calculation continues, the calculation process may be interrupted when the time step is smaller than ti, which makes the blow molding simulation of balloon unable to continue.
In fact, when ti is one order of magnitude smaller than t0 in the process of setting software parameters, the iterative calculation of blow molding balloon can proceed smoothly. And the phenomenon that the error is reported when the time step of transient iterative calculation is smaller than ti does not occur.
Under any given blow molding pressure, a t0 can be obtained. In order to obtain the relative position relationship between the parison and the mold cavity at the end of the first iteration under different blow molding pressure conditions, the statistics of different given initial values of iterative time step are shown in
| Blow molding pressure (MPa) | Given initial values of iterative time step when the parison is inside the cavity after the first iteration (×10−5 s) | Given initial values of iterative time step when the parison is outside the cavity after the first iteration (×10−5 s) |
|---|---|---|
| 0.1 | 8.29 | 8.30 |
| 0.2 | 5.73 | 5.74 |
| 0.3 | 4.78 | 4.79 |
| 0.4 | 4.14 | 4.15 |
| 0.5 | 3.70 | 3.71 |
| 0.6 | 3.38 | 3.39 |
| 0.7 | 3.13 | 3.14 |
| 0.8 | 2.93 | 2.94 |
| 0.9 | 2.76 | 2.77 |
| 1.0 | 2.62 | 2.63 |
| 1.1 | 2.50 | 2.51 |
| 1.2 | 2.39 | 2.40 |
| 1.3 | 2.29 | 2.30 |
| 1.4 | 2.21 | 2.22 |
|---|---|---|
| 1.5 | 2.14 | 2.15 |
| 1.6 | 2.07 | 2.08 |
| 1.7 | 2.01 | 2.02 |
| 1.8 | 1.95 | 1.96 |
| 1.9 | 1.90 | 1.91 |
| 2.0 | 1.85 | 1.86 |
| 2.1 | 1.80 | 1.81 |
| 2.2 | 1.76 | 1.77 |
| 2.3 | 1.72 | 1.73 |
| 2.4 | 1.69 | 1.70 |
| 2.5 | 1.65 | 1.66 |
| 2.6 | 1.62 | 1.63 |
| 2.7 | 1.59 | 1.60 |
| 2.8 | 1.56 | 1.57 |
| 2.9 | 1.53 | 1.54 |
| 3.0 | 1.51 | 1.52 |
| 3.1 | 1.48 | 1.49 |
| 3.2 | 1.46 | 1.47 |
| 3.3 | 1.44 | 1.45 |
| 3.4 | 1.42 | 1.43 |
| 3.5 | 1.40 | 1.41 |
| 3.6 | 1.38 | 1.39 |
| 3.7 | 1.36 | 1.37 |
| 3.8 | 1.34 | 1.35 |
| 3.9 | 1.32 | 1.33 |
| 4.0 | 1.31 | 1.32 |
| 4.1 | 1.29 | 1.30 |
| 4.2 | 1.27 | 1.28 |
| 4.3 | 1.26 | 1.27 |
| 4.4 | 1.25 | 1.26 |
| 4.5 | 1.23 | 1.24 |
| 4.6 | 1.22 | 1.23 |
| 4.7 | 1.20 | 1.21 |
| 4.8 | 1.19 | 1.20 |
| 4.9 | 1.18 | 1.19 |
| 5.0 | 1.17 | 1.18 |
In this paper, the numerical simulation of a single balloon molding process is studied by POLYFLOW in a mold core cavity with L = 15.0 mm and D = 7.0 mm. And the following conclusions are obtained:
1) As Pb increases, t0 decreases rapidly at first and then slowly. At the end of the first step of iterative calculation for each t0, the balloon parison is in the mold core cavity.
2) If ti is smaller than t0, the balloon parison is still in the mold core cavity at the end of the first iteration. However, if the iterative calculation continues, the calculation process may be interrupted when the time step is smaller than t0, which makes the blow molding simulation of balloon unable to continue.
3) It is suggested that ti is one order of magnitude smaller than t0 to complete smoothly the simulation of blow molding balloon.
Project Supported by Scientific Research Fund of Chaohu University (KYQD- 201908).
The authors declare no conflicts of interest regarding the publication of this paper.
Bian, J.Y., Jin, G.B., Qin, P. and Dai, G.H. (2020) Influence of Initial Value of Transient Time Step on Numerical Simulation of Blow Molding Balloon. American Journal of Computational Mathematics, 10, 266-274. https://doi.org/10.4236/ajcm.2020.102014