Process variables for PET in a typical melt-spinning operation.
Abstract
This chapter discusses an approach for process simulation in the design of melt spinning process for finding optimal design parameters concerning spinneret, quench air unit and other technical parameters for maximum throughput and quality. The property of as-spun fiber is a function of structural parameters at a given condition and orientation of the structural parameter and it is highly governed by stress level at freeze line. Thus, to define structural property and associated relationship, it requires to identify the process to control the variables (or factors) that affect the structural parameter as well as final fiber property. In addition, this chapter also provides a System-of-Systems (SOS) perspective on melt spinning process and its computer modeling along with mathematical equations for estimating spinline stress with a change in process variables. The spinline stress will be used as an input for a computer simulation to have process optimization by changing the necessary variables until it optimized.
Keywords
- melt spinning
- system-of-systems
- poly ethylene terephthalate
- simulation
- process parameters
- process optimization
1. Introduction
Modeling and Simulation (M&S) have become important tools for evaluating and scheming a melt spinning in a comprehensive arrangement of disciplines fluctuating from fiber spinning and engineering to melt spinning [1]. For example, in engineering scheme, modeling the parameters and simulation used to evaluate the effectiveness of a melt spinning process concept, verify whether all the functional design specifications are meet, or suggest modifications for improving the manufacturability of a product [2]. Melt spinning processes are considered as a system-of-systems process based on the thermoplastic polymer processing method as shown in Figure 1 [3]. For this process a system is required to ensure the polymer is melted above its melting temperature, and the melted polymer is then transported to another system where it is metered for constant mass through a spinneret into a quench air stream blowing across the spinline [4]. The Spinline-and-Free-line” forming can be considered as another system that is used to form spinline, cool and finally solidify at a distance from the spinneret called ‘freeze- line’ [5]. The solidified polymer filament is winding at a speed significantly greater than the extrusion velocity; this cause the final cross-sectional area is considerably smaller (about 100–200 times) than the initial extruded area. This “Winding” process can be considered as another system in the melt spinning process [2, 6]. The way we decompose the melt spinning process into SOS view is based on the “System Control” view-point. For examples: (i) the melted polymer process is considered as a system because it requires a temperature control to ensure the polymer melts above its melting temperature, and (ii) the winding process is considered as another system because it requires to control the winding speed ensuring that the winding speed is significantly greater than the extrusion velocity [7]. Concerning “Spinline-and-Free-line” forming process, note that without considering crystallization in the spinline, solidified filament orientation can be represented by the orientation in the spinning threadline frozen at or near the glass transition temperature of the polymer due to result of tensile deformation under stress in the flowing melt. So, stress level at freeze line can easily govern the orientation, in which the process variables such as extrusion temperature, extrusion velocity, take-up speed, quench air velocity and temperature affect the stress level at freeze line. This shows that molecular orientation in the as-spun filaments is a function of extensional deformation and cooling. Cooling rate across may vary across the filament and cause a change in oriented filament cross section [8]. Thus, “Spinline-and-Freeze-line” system requires a complex system control to ensure accurate orientation, glass transition temperature, and tensile deformation under stress. In practice, the automation in melt spinning can be tested using simulation approach for estimating and gaining insight of the changes of process variables in order to reduce the risk on uniformity of the fiber quality [9].
2. Spinline orientation
Quantification of melt spinning requires the development of mathematical equations describing processing-structural-property correlations. The mathematical equations will be used in the computer simulation for predicting the fiber-line stress using readily measurable process variables such as melt temperature, take-up speed, and extrusion velocity and quantitative correlations between response variable like fiberline stress and structural parameter of the solidified fiber [5]. Using a simple quantitative expressions critical operating conditions and material properties for the melt-spinning for Polyethylene terephthalate (PET) identified. The method used depends on the experimental observation that for vitrified polymers like PET spun at speeds below 3000 m/min, the molecular orientation of the as-spun filament is uniquely determined by stress at freeze line when spinline stress is at the glass transition temperature [9, 10].
Clearly, the spinning stage represents a non-isothermal, uniaxial elongational flow situation with variable physical properties [5]. In order to develop the governing mathematical equations, the following assumptions will be used through this chapter.
The process is operating under steady state conditions.
The temperature and velocity field are independent of radial position.
The elongational viscosity is independent of extension rate and is equal to the Trouton viscosity.
All the molecular motions at temperatures less than the glass transition temperature, Tg.
The effects of inertia, gravity, surface tension and air drag are negligible.
Dieswell has insignificant effect.
Figure 1 shows the coordinate system used for the analysis. Based on the above assumptions, the governing equations can be written as:
or alternatively
and
The appropriate boundary conditions are
Where L is the distance of the freeze line from the spinneret. The temperature, T(x), the velocity V(x), and the spin line tension F, are obtained by solving Eqs. (2), (3), and (4) subject to the conditions represented by Eqs. (5) and (6). In order to accomplish that, however, correlations for the temperature dependence of physical properties (ρ, Cp ‘
Rheological measurements of the shear viscosity of PET indicate that it behaves like a Newtonian liquid for shear rates up to about 200 S−l. Thus, the assumption of Newtonian behavior may not be too inappropriate for PET [9, 10, 11, 12]. The variations of zero shear viscosity of PET,
For density and specific heat, we assume the following linear variations with temperature [13].
For the heat transfer coefficient, the best available correlation is [9].
A proper choice of F is crucial for smooth operation of the numerical procedure; a fairly good starting estimate can be obtained from the following equation.
At intermediate take-up speeds (3000 m/min), melt-spinning of PET results in almost amorphous as-spun fibers. That is why molecular orientation of the as-spun fiber is expected to represent the ‘frozen-in’ stress in the spinline [1, 2, 14]. Molecular orientation which is measured by birefringence quantitatively related to process variable form a link between the latter and the as-spun fiber properties.
Experimental data shown in Figure 2, for melt-spun PET fibers under different process conditions demonstrate that the birefringence is directly proportional to the spinline stress for PET,
C being the stress optical coefficient [9]. The intermediate take-up speeds in a melt-spun PET fibers show amorphous at a slow crystallization behavior of PET, and the birefringence depends on the (experimentally measured) spinline stress. To show the relationship between as-spun fiber birefringence and the calculated stress at the freeze line, the stress at freeze line is taken to be the spinline stress at which the velocity is 95% of the take-up velocity [5, 15]. It is observed that a linear variation is obtained with stress optical coefficient approximately about 6 × 10−9 to 9 × 10−9 Pa −1, which is then compared with the values obtained by different experiments [5]. For other experiments, the spinline stress controls the orientation of the as-spun fiber, which in turn determines the mechanical properties and associated fiberline processing ability. Thus, it is clear that stress at freeze-line controls the orientation of the as-spun fiber, and in turns it can determine the mechanical properties [4, 13]. The sensitivity of the as-spun fiber properties is directly related to the sensitivity of freeze-line stress to the changes in process variables. Stress at freeze-line, on the other hand, can readily be obtained from the process simulation calculations, provided the operating conditions and the melt properties are available [16].
3. Sensitivity analysis
Sensitivity analysis provides an approach to look into the effect of a factor affecting the response variables and making effective changes in the process parameters for controlling process stability and product quality improvement [1, 6]. As shown in Table 1, the critical process variables have a significant effect on stress level at the freeze-line are quantified for a typical PET spinning process [17]. The spinline stress is calculated under different operating conditions and melt properties, which is then used in the assessment of stress sensitivity, are shown in Figures 3 and 4 below [9].
Process variable | Value |
---|---|
Intrinsic velocity | 0.611 dl g−1 |
Extrusion temperature | 285°C |
Glass transition temperature | 67°C |
Ambient air temperature | 25°C |
Quench air velocity | 48 m min−1 |
Spinneret hole diameter | 0.00038 m |
Take-up velocity | 1500 m min−1 |
Mass flow rate | 0.001 kg min−1 |
Spinline tensiona | 0.000972 N |
Freeze line location from spinneret | 0.37 m |
The results indicate that the changes in process variables like extrusion temperature, melt, take-up velocity and the melt throughput rate that can affect the response variables (e.g., stress at the freeze-line) [5, 8]. The slope of the curve provides a measure of the sensitivity of stress at the freeze-line to process variable changes. In comparison, stress at the freeze-line is moderately sensitive to changes in quench air velocity and relatively insensitive to alterations in the quench air temperature [9]. But it is more sensitive to a changes in the throughput rate as well as the extrusion temperature and is relatively unperturbed by changes in the melt intrinsic viscosity. But other changes exert only moderate influence in determining the freeze-line location [6]. Since the as-spun fiber quality is uniquely determined by its birefringence, which in turn is directly related to stress at the freeze-line, so the degree of sensitivity of the as-spun fiber properties changes in these process variables are the same as that of stress at the freeze-line [15, 18, 19].
4. Process implications
The information obtained from the sensitivity analysis provides good insight for identifying the critical process with associated level of risk [14]. The process for controlling the as-spun filament orientation requires critical process parameters like extrusion temperature, intrinsic viscosity of the polymer melt take-up velocity and melt flow rate. These parameters need to be taken into consideration, unless the control process quality is having issues as discussed below [5].
If there is an increase in temperature and/or a decrease in melt intrinsic viscosity, the spinline stress is being easily reduced, which affects the stretch ability as well as final fiber quality [20]. These higher sensitivity of the spinline stress in response to extrusion temperature has important implications for maintaining uniform product quality in a multiposition commercial plant with each position carrying multiple hole spinnerets [21, 22]. As shown in Figure 5, Curve 4, a variation of about 3°C in quench air temperature profile of the spinneret represents a 1% change in the extrusion temperature; but this could result in about 10% variation in the spinline stress by increasing the coefficient of variation of elongation [9]. So, having static mixers upstream of the spin pack is the best way for delivering a thermally homogeneous melt to the spinnerets [21, 22, 23]. Both the melt flow rate and the take-up velocity can also affect as-spun filament quality. This is due to, increase in take-up velocity at a constant flow rate results in increase orientation due to enhanced stretching. But for a constant take-up velocity, the orientation is reduced as the flow rate is enhanced. Because the filament denier increases rapidly with increase in flow rate, and even the spinline tension actually decreases proportionally with a net effect in a reduction in spinline stress as throughput is increased [14, 20]. By controlling such process variables (i.e., the melt flow rate and the take-up speed), it is possible to produce products over a range of denier and tenacity if needed [9, 17]. In addition, if the process is a batch process, the moisture content needs to be controlled less than 0.005% and even a moisture level of 0.01% could lower the polymer intrinsic viscosity from 0.6 to 0.56. In case if the change in intrinsic viscosity is 7%, which leads to a 35% change in spin line stress, this thereby affecting spinnability and product quality that desired [5].
In plants with continuous polymerization, variation in polymer intrinsic viscosity may also result from reactor perturbation [22]. Therefore, incorporation of an on-line viscometer to monitor the melt quality is likely to be a necessary instrumentation. There are different ways of increasing productivity. One approach is to use a process optimization discussed above, i.e., by changing the process variables to find optimum solution. The other approach is to increase in the number of holes in the spinneret by 30% [24]. If possible, as per hole layout guidelines depend on the investment to manipulate such multi-hole spinneret [9, 23, 24]. The location of the freeze-line has critical practical implications for spinline stability and product uniformity in multifilament spinning air turbulence [1, 9, 22]. An incorrect location of the freeze-line can lead to filament breakage and/or fused filaments resulting in process interruptions and poor product quality. Reducing the quench air temperature from 25 to 15°C show 40% change in the filaments stability by moving the freeze-line closer to the spinneret without a significant effect on the structural parameter like orientation, this makes it be fewer disturbances to air [9]. The take-up velocity is not the most dominant controlling variable for orientation development since its process implication does not have much impact on spun-fiber.
5. Stress-orientation relationship for PET
The Dutta Nadkarni simulation package is recommended for simulation of the melt spinning process [9, 20, 23]. Before we use the computer simulation for spinning process optimization, it is important to develop a way to quantify the correlation between the spinline stress and the as-spun filament orientation, such that the operator can confirm its validity over a wide range of spinning process parameters in current industrial operations is the data are collected from three different plants for computing the spinline stress, stress at the freeze-line [5, 25]. Table 2 summarizes the industrial data collected. Using the data birefringence’s, of the as-spun filaments orientation are measured experimentally (Figure 6).
Process parameter | Range | No. of variables |
---|---|---|
Intrinsic velocity | 0.6–0.65 | 2 |
Extrusion temperature | 275–290 | 3 |
Quenching air temperature | 15–25 | 7 |
Quench air velocity | 30–130 | 3 |
Spinneret hole diameter | 0.2–0.3 | 11 |
Take-up velocity | 500–1250 | 9 |
Throughput rate | 0.5–1.0 | 2 |
Filament diameter | 4.0–12.0 | 11 |
The above equation is used as an input to the melt spinning process for optimizing the relationship between birefringence’s and spinline stress for the melt-spinning of PET at speeds up to 3000 m min−1 with a correlation coefficient of 0.96.
6. Case study for process optimization
The use of Eq. (1) along with the recommended computer simulation package for a case study to: (i) Increase the production rate, (ii) Quality of the fiber, and (iii) a new desired product developed through process optimization in spinning without affecting the fiber line processing conditions [9, 23]. Maintain the denier and orientation of the as-spun filaments at the reference levels because for a given denier product, the throughput rate and the take-up speed have to be increased proportionally which avoid downstream process optimization at the higher productivity [15]. For cases to increase in productivity changing the polymer, viscosity is not desirable rather other process parameters are desirable, such as extrusion temperature, quench air velocity and quench air temperature [7, 26]. According to the process sensitivity analysis, the orientation of the as-spun filaments, as measured by birefringence is very sensitive to changes in the spinning temperature and since the orientation related to freeze-line stress, the spinning temperature is to be changed so as to maintain freeze line stress and spin line tension at the original values [21, 25]. Quench air velocity and temperature are the other process variables that can be changed but they do not show significantly affect the spinline stress, rather on the freeze-line location in order to ensure stability of the spinline to avoid denier variation and final fiber quality [3, 20, 27]. The production rate can also be improved by a proportional increase in the number of spinneret holes with no changes in throughput per hole, take-up speed and spinning temperature and even it is best option than process optimization through changing variable [11].
The reference process parameters of a spin- line whose productivity is to be improved are summarized in Table 3 [9] along with the process conditions determined for production rate increases of 10, 20 and 30%. The process variables provided in this table are representative of values used in a typical spinning operation, as averaged out from within the ranges given in Table 2. The computer simulation of Dutta and Nadkarni is used for determining the changes in the process variables when the throughput rate and take-up speed changes by 10, 20 and 30% and the effect on the response parameters, including spinline stress, spinline tension and freeze-line location, for increased productivity is shown in Figure 7 [11].
Parameter | Original condition | Increase in melt flow rate, take up velocity(W,VL) | ||
---|---|---|---|---|
10% | 20% | 30% | ||
Intrinsic velocity | 0.6 | 0.6 | 0.6 | 0.6 |
Extrusion temperature | 280 | 280 | 280 | 280 |
Quenching air temperature | 20 | 20 | 20 | 20 |
Quench air velocity | 100 | 100 | 100 | 100 |
Spinneret hole diameter | 0.2 | 0.2 | 0.2 | 0.2 |
Take-up velocity | 1000 | 1100 | 1200 | 1300 |
Throughput rate | 0.5 | 0.55 | 0.6 | 0.65 |
Spinline stress | 2.18 | 2.32 | 2.46 | 2.60 |
Spin line tension | 82.8 | 87.5 | 93.8 | 98.0 |
Freeze line location | 26.4 | 28.0 | 30.8 | 33.0 |
All process variables show proportional increments with increased throughput and take-up velocity [15]. Thus, to maintain the molecular orientation and the downstream process variables constant, it is essential to have the values of spinline stress, spinline tension and freeze-line location with proper adjustment of other process variables [25]. The sensitivity analysis indicates that increasing extrusion temperature decreases the spinline stress and spinline tension but extends the freeze-line location as shown in Figure 8 below.
When extrusion temperature is assumed to be equal to the spinneret exit temperature, there is a 10% increase in productivity. In order to return the spinline stress and tension from a 6% increase to the original level, the increase in the extrusion temperature should be approximately 1.87% [5, 7]. However, a 1.87% increases in extrusion temperature would result in a further 1.75% increase in Freeze-line location. The freeze-line location is thus sensitive to the combined changes in throughput rate and extrusion temperature. [17, 24] Referring to Table 4, although the spinline stress and tension can be brought back to the original values by a change in the extrusion temperature, the freeze-line location remains different from the original value by about 9.5%. The location of the freeze-line can affect the spinline stability and product uniformity [5, 25]. Hence its retention at the original value is desirable. A process variable is to be sought, which would move the location of the freeze line closer to the spinneret without any significant changes in either spinline stress or Freeze-line location. The quench air temperature is known to influence freeze-line location to a greater extent relative to spinline tension and spinline stress and can be effectively used for stabilizing the spinline. From Figure 9 it is seen that to move the freeze- line location by about 8.0%, the quench air temperature has to be decreased by about 45%. Furthermore, the effect of this process change on spinline stress and freeze line location is insignificant. The entire exercise is summarized in Table 4, which indicates the changes in process variables required to increase productivity while maintaining the denier, orientation and stability of the spinline [5].
Parameter | Reference condition intrinsic viscosity (IV) = 0.6 | 10% decrease in freeze line location, take up velocity(L,W) | 1.9% increase in extrusion temperature (TE) | 45% decrease in quench air temperature (Ta) |
---|---|---|---|---|
Intrinsic velocity | 0.6 | 0.6 | 0.6 | 0.6 |
Extrusion temperature | 280a | 280 | 285 | 285a |
Quenching air temperature | 20a | 20 | 20 | 11a |
Quench air velocity | 100 | 100 | 100 | 100 |
Spinneret hole diameter | 0.2 | 0.2 | 0.2 | 0.2 |
Take-up velocity | 1000a | 1100 | 1100 | 1100a |
Throughput rate | 0.5a | 0.55 | 0.55 | 0.55a |
Spin line stress | 2.18b | 2.32 | 2.20 | 2.19b |
Spin line tension | 82.8b | 87.9 | 83.3 | 82.9b |
Freeze line location | 26.4b | 28.0 | 28.9 | 26.8b |
7. Case study for product development of low pill fibers
Lowering the intrinsic viscosity or molecular weight by having drawn fiber denier of 1.2 held constant to make the downstream process variables at constant values [12, 28] (Table 5). Because of this the changeable variables, it is required to lower the draw ratio to 3.5 from reference value to prevent fiber breakage in downstream process and this change as-spun denier from 4.5 to 4.2. Due to 8.3% decreases in intrinsic viscosity a new pill fiber is developed. However, this phenomena (decrease in spun fiber diner) causes the decrease in Spinline stress and Spinline tension, which makes melt unspinnable increases the take-up velocity and the throughput rate to new values (22.5%, 14.2% respectively), and lowering extrusion temperature to 275°C. This is a way to increase the Spinline tress with as-spun denier at a value of 4.2. Freeze-line location has effect on the spinline stability and product uniformity stress [29]. However having 50% decreases in quench air temperature from reference value, to maintain the spinline stress and freeze-line location at the original values, we need to move the location of the freeze-line closer to the spinneret [9, 11, 12, 30].
Parameter | Reference condition intrinsic viscosity (IV) = 0.6 | 8.3% decrease in melt intrinsic viscosity (IV) | 22.5% increase in take up velocity (VL) 14.2% increase in melt flow rate (W) | 1.8% decrease in extrusion temperature (TE) | 50% decrease in quench air temperature (Ta) |
---|---|---|---|---|---|
Intrinsic velocity | 0.6a | 0.55 | 0.55 | 0.55 | 0.55a |
Extrusion temperature | 280a | 280 | 280 | 275 | 275a |
Quenching air temperature | 20a | 20 | 20 | 20 | 10a |
Quench air velocity | 100 | 100 | 100 | 100 | 100 |
Spinneret hole diameter | 0.2 | 0.2 | 0.2 | 0.2 | 0.2 |
Take-up velocity | 1000a | 1000 | 1225 | 1225 | 1225a |
Throughput rate | 0.5a | 0.5 | 0.571 | 0.571 | 0.571a |
Spin line stress | 2.18b | 1.63 | 2.05 | 2.16 | 2.18b |
Spin line tension | 82.8a | 61.8 | 72.5 | 76.6 | 76.9a |
Freeze line location | 26.4a | 26.3 | 29.6 | 29.3 | 26.5a |
As-spun fiber denier | 4.5a | 4.5 | 4.2 | 4.2 | 4.2a |
Drawn fiber denier | 1.2b | 1.2 | 1.2 | 1.2 | 1.2b |
Draw ratio | 3.75a | 3.75 | 3.5 | 3.5 | 3.5a |
8. Conclusion
The mechanical property of fiber is highly governed by the orientation level in the thread line. However, this structural parameter is heavily affected by the melting process as shown in the sensitivity analysis above with different level of effects. According to the above analysis, it is shown that the molecular orientation developed in a melt flow in spinning process is governed by the spin-line stress (stress at freeze line). A mathematical equation, which is an input for the simulation is developed according to a given quantitative data, shows a linear relationship between the orientation and stress levels. The simulation results are correlated with the variables associated with the spin-line stress to determine if the stated assumptions are in agreement with the simulation results. The combination of the process simulation and stress-orientation relationship gives us a procedure for identifying and simulating the important spinning process variables that would affect the spin-line stress, and hence the mechanical properties and fiber line processability of the as-spun filaments. Once the effect is simulated the process optimization is done by changing the process variables according to the final property desired. Therefore, the process optimization can be achieved using SOS perspective and computer simulation. This is a good advancement in melt spinning to achieve the process optimization. The future work will tell on computer simulation and modeling of manmade fibers and their parameters based on systems of system perspective.
Notations
initial cross sectional area
spinneret exit temperature
cross sectional area at (x) distance
spinneret temperature at (x) distance
quench air velocity
final velocity
quenches air temperature
final cross sectional area
distance from spinneret
freeze line location
glass transition temperature
birefringence
stress at freeze line
References
- 1.
Baran I, Cinar K, Ersoy N, Akkerman R, Hattel JH. A review on the mechanical modeling of composite manufacturing processes. Archives of Computational Methods in Engineering. 2017; 24 (2):365-395 - 2.
Barshop BA, Wrenn RF, Frieden C. Analysis of numerical methods for computer simulation of kinetic processes: Development of KINSIM—A flexible, portable system. Analytical Biochemistry. 1983; 130 (1):134-145 - 3.
Bungartz HJ, Schafer M. Modelling. Simulation. Optimization. 2014; 53 - 4.
Doufas AK, McHugh AJ, Miller C. Simulation of melt spinning including flow-induced crystallization: Part I. Model development and predictions. Journal of Non-Newtonian Fluid Mechanics. 2000; 92 (1):27-66 - 5.
Gupta VB, Kothari VK e. Manufactured fiber technology. In: Springer Science & Business Media. 2012 - 6.
Mattrand C, Beakou A, Charlet K. Numerical modeling of the flax fiber morphology variability. Composites Part A: Applied Science and Manufacturing. 2014; 63 :10-20 - 7.
Spagele T, Kistner A, Gollhofer A. Modelling, simulation and optimisation of a human vertical jump. Journal of Biomechanics. 1999; 32 (5):521-530 - 8.
Lin S, Ryu S, Tokareva O, Gronau G, Jacobsen MM, Huang W, et al. Predictive modelling-based design and experiments for synthesis and spinning of bio inspired silk fibres. Nature Communications. 2015; 6 (1):1-12 - 9.
Nadkarni, V.M., 1997. Computer simulation of melt-spinning. In Manufactured Fibre Technology (pp. 98–123). Springer, Dordrecht - 10.
Singha K, Maity S, Singha M, PAL S. Computer simulations of textile non-woven structures. Frontiers in Science. 2012; 2 (2):11-17 - 11.
Schutze M, Butler D, Beck BM. Modelling, Simulation and Control of Urban Wastewater Systems. Springer Science & Business Media; 2011 - 12.
Shengwang Y. Influence of TCS process conditions on pet FDY formation [j]. China Synthetic Fiber Industry. 2001; 5 - 13.
Cao C, Chung TS, Chen SB, Dong Z. The study of elongation and shear rates in spinning process and its effect on gas separation performance of poly (ether sulfone) (PES) hollow fiber membranes. Chemical Engineering Science. 2004; 59 (5):1053-1062 - 14.
Bell WP, Edie DD. Measured orientation and internal stress distributions in melt-spun fibers. Journal of Applied Polymer Science. 1987; 33 (4):1089-1102 - 15.
Kikutani, T., 2009. Structure development in synthetic fiber production. In Handbook of Textile Fibre Structure (pp. 157–180). Woodhead Publishing - 16.
Muller-Plathe F. Molecular dynamics simulation of gas transport in amorphous polypropylene. The Journal of Chemical Physics. 1992; 96 (4):3200-3205 - 17.
Hehenberger P, Vogel-Heuser B, Bradley D, Eynard B, Tomiyama T, Achiche S. Design, modelling, simulation and integration of cyber physical systems: Methods and applications. Computers in Industry. 2016; 82 :273-289 - 18.
Clancy TC, Mattice WL. Computer simulation of polyolefin interfaces. Computational and Theoretical Polymer Science. 1999; 9 (3–4):261-270 - 19.
Marla VT, Shambaugh RL. Three-dimensional model of the melt-blowing process. Industrial & Engineering Chemistry Research. 2003; 42 (26):6993-7005 - 20.
Bhuvanesh, Y.C. and Gupta, V.B., 1990. Computer Simulation of Melt Spinning of Poly (Ethylene Terephthalate) Using a Steady-State Model - 21.
Cicek, H., 1990. Computer Simulation of Melt Spinning - 22.
Elser M, Tschegg EK, Finger N, Stanzl-Tschegg SE. Fracture behaviour of polypropylene-fibre-reinforced concrete: Modelling and computer simulation. Composites Science and Technology. 1996; 56 (8):947-956 - 23.
Toth R, Coslanich A, Ferrone M, Fermeglia M, Pricl S, Miertus S, et al. Computer simulation of polypropylene/organoclay nanocomposites: Characterization of atomic scale structure and prediction of binding energy. Polymer. 2004; 45 (23):8075-8083 - 24.
Huang N, Yu X, Tang Z, Hou P, Xu X. Computer simulation of polyester high-speed thermal channel spinning. Macromolecular Materials and Engineering. 2002; 287 (5):330-338 - 25.
Gagon DK, Denn MM. Computer simulation of steady polymer melts spinning. Polymer Engineering & Science. 1981; 21 (13):844-853 - 26.
Liang X, Ding Y, Wang Z, Hao K, Hone K, Wang H. Bidirectional optimization of the melting spinning process. IEEE transactions on cybernetics. 2013; 44 (2):240-251 - 27.
Bhuvanesh YC, Gupta VB. Computer simulation of melt spinning of polypropylene fibers using a steady-state model. Journal of Applied Polymer Science. 1995; 58 (3):663-674 - 28.
Cooke WD. 21—Fibre fatigue and the pilling cycle part iii: Pill wear-off and fabric attrition. Journal of the Textile Institute. 1984; 75 (3):201-211 - 29.
Ukponmwan JO, Mukhopadhyay A, Chatterjee KN. Pilling. Textile Progress. 1998; 28 (3):1-57 - 30.
Sriprateep K, Bohez ELJ. A new computer geometric modelling approach of yarn structures for the conventional ring spinning process. The journal of the textile Institute. 2009; 100 (3):223-236