[BACK]
 Computer Modeling in Engineering & Sciences

DOI: 10.32604/cmes.2021.015153

ARTICLE

Modeling Additional Twists of Yarn Spun by Lateral Compact Spinning with Pneumatic Groove

1Key Laboratory of Textile Science & Technology, Ministry of Education, College of Textiles, Donghua University, Shanghai, 201620, China
2Merchandising and Digital Retailing, University of North Texas, Texas, 76203, USA
*Corresponding Authors: Longdi Cheng. Email: ldch@dhu.edu.cn; Bugao Xu. Email: bugao.xu@unt.edu
Received: 25 November 2020; Accepted: 26 February 2021

Abstract: Compact spinning with pneumatic grooves is a spinning process to gather fibers by blended actions of airflow and mechanical forces. Modified from the ring spinning system, the lateral compact spinning with pneumatic grooves can improve yarn appearance and properties due to generated additional twists. In this study, we investigated additional twists of the lateral compact spinning with pneumatic grooves via a finite element (FE) method. An elastic thin rod was used to model a fiber to simulate its dynamic deformation in the three-dimensional space, and the space bar unit was used to simplify the fiber model for the dynamic analysis. The stiffness equation of the elastic rod element and the dynamic equation of the rigid body mass element were derived from the differential equation of the elastic thin rod. In the analysis of the nonlinear geometric displacement of the space elastic thin rod unit, the large deformation problem was solved with the stepwise loading successive approximation. The simulation results explained the mechanism of generating additional twists, and the experiment results proved the existence of additional twists. The study demonstrated that the FE model is effective for predicting additional twists of fiber bundles in the agglomeration zone, and for simulating the fiber motion in the compact spinning with pneumatic grooves.

Keywords: Compact spinning; additional twist; finite element; dynamic analysis

1  Introduction

Yarn twist, or the number of twists per unit length of yarn, directly affects yarn’s structure and mechanical properties, and fabric styles [1,2]. In studies of compact spinning, Zhou et al. [3,4] discovered twists in the condensing zone of pneumatic compact spinning with lattice aprons, and named these twists as “additional twists.” The reason why additional twists can be generated when a fiber bundle goes through the condensing zone is the force of lateral groove [5].

In order to study additional twists, Feldman [6] used material points as a fiber model to establish equations of the fiber motion in the condensing field and to obtain the trajectory of the fiber in a two-dimensional and non-uniform flow field. Bangert et al. [7] regarded fibers as rigid cylinders when describing motion and orientation of fibers under flow field forces. This was a relatively simple fiber model in the early time, which had drawbacks of no expression on the fiber rigidity. Jeffery [8] first proposed an orbital theory to solve the motion of a single rigid ellipsoid suspended in Newtonian fluid of a simple shear flow. Jeffery’s orbital theory has been widely used in subsequent fiber dynamics studies, and was enriched and perfected by Anezurowski et al. [9], Chiba et al. [10,11] and other researchers. In Jeffery’s orbital theory, the fiber was regarded as a rigid particle to represent the position and orientation of the fiber. However, this rigid particle was unable to reflect the flexibility and deformation of the fiber. Skjetne et al. [12] proposed a fiber model that consisted of a series of ellipsoids or spheres connected by “bead-bowl joints.” Nyland et al. [13] proposed a fiber model consisting of a plurality of rigid segments, which were connected by “bead-bowl joints.” Cheng [14,15] established a chain-ball model to analyze fiber motion with the movement patterns and deformations of fibers in incompressible viscous fluids. But Cheng’s model was too complex, and the viscoelastic property was not the main physical property of the fiber, although the research provided great inspiration and guidance for future studies on fiber models. Yamamoto et al. [16,17] proposed a bead-spring-chain model by arranging spheres connected by springs in series. Zeng et al. [18,19] treated a fiber as multiple elastic and flexible chains, and illustrated the elastic elongation of the fiber with changes in the mass point spacing of the chain. They also studied bending behaviors of the fiber with changes in the deflection of the chain. Wang [20] used a rigid micro-element segment as a fiber model to simulate the fiber motion in compact spinning with a lattice apron. The aforementioned fiber models regarded a fiber as multi-rigid chains to reflect the bending deformation characteristics, but these multi-rigid chains could not illustrate the elastic features of the fiber [21].

There is still a large gap between the performances of these existing fiber models and the actual mechanical behavior of a fiber. In this paper, we would like to introduce a continuous elastic thin rod to study the structure and dynamics of the fiber flexible body. The simultaneous partial differential equations of general nonlinear mechanics of the elastic thin rods are used to express the small deformation behavior of the rods. Based on the spatial elastic thin rod unit, we establish the fiber finite element (FE) model to explore the additional twists of lateral compact spinning with pneumatic grooves. The FE model of continuous elastic fine rods is used to simulate the fiber motion in compact spinning for both efficient calculations and effective visualization in the three-dimensional space. The FE method is applied to solving the dynamic problem of the large deformation of an elastic thin rod. The stiffness equation of the elastic rod element and the dynamic equation of the rigid body are derived from the dynamic analysis of a combination of mass units and elastic rod units. The nonlinear geometric large deformation problem of the spatial elastic thin rod unit is studied by solved the method of step-by-step loading and successive approximation. The goal of the research is to present an effective and feasible theoretical model and method for examining the twist formation process and mechanism of flexible fibers under the mechanical force of lateral grooves.

2  Methodology

2.1 Establishment of Fiber Finite Element Model

Fig. 1a shows a roller system that draws and gather fiber bundles to form a desired yarn. And Figs. 1b and 1c are the groove shapes of the intermediate compact spinning with pneumatic groove and the lateral compact spinning with pneumatic groove.

Figure 1: Diagram of compact spinning devices. (a) Gathering process of fiber bundle, (b) the intermediate compact spinning with pneumatic groove, (c) the lateral compact spinning with pneumatic groove

In order to solve the nonlinear large deformation problem and to realize the large draw ratio (fiber length/fiber diameter) of fibers, a finite element model is established with spatial elastic thin-rod units., as shown in Figs. 2 and 3. Fig. 2 presents the micro-segment rigid body of the FE model that shows the general motion of node i, while Fig. 3 displays elastic thin rods in the finite element model. In Fig. 2, O-xyz is an overall coordinate system, and O’-uvw is a connected moving coordinate system of a rigid body, or a local coordinate system. The origin of the coordinate is located at the centroid of the rigid body, and the direction of a coordinate axis is the direction of the principal axis of the mass center of the rigid body. The u axis is in the axial direction of the micro-segment, the v and w axes are located on the cross section of the micro-segment rigid body through the centroid.

Figure 2: Rigid body unit for general motion in space

Figure 3: Elastic thin rods in the finite element model [22]

To analyze the dynamics of the elastic thin rod, the elastic thin rod is divided into n independent micro-segments and n + 1 nodes. Each node relative to the overall coordinate system has six degrees of freedom. Six coordinate parameters are needed to determine the relative position of a node. Three of the six parameters describe the moving line displacement of the node, and the other three describe the cross-sectional rotational angular displacement of the node. When the mass of the elastic thin rod is considered, the mass of each micro-segment is regarded as a rigid body moving with its local coordinate system (see Fig. 2). When the elasticity of the elastic thin rod is considered, the connection between two adjacent nodes is regarded as one elastic rod unit that is massless (see Fig. 3). The external force on the elastic thin rod is simplified so that it is only exerted on the two end nodes.

2.2 Dynamics of Fiber Model

Any general motion of the rigid body in space can be decomposed into the translation with its centroid and the rotation relative to the centroid. According to the centroid motion theorem, the dynamic equation of the rigid body translation with centroid O’ is

ma=F(1)

where m is the mass of the rigid body, a is the acceleration relative to the whole coordinate system, and F is the external force on the rigid body.

According to the momentum quadrature theorem, the dynamic equation of the rigid body rotating around the centroid O’ can be derived as

Jα+ω×(Jω)=M(2)

where J is the inertance of the rigid body to the centroid, M is the main quadrature of the external force on the centroid of the rigid body, ω is the angular velocity of the rigid body, and α is the derivative of ω to time. By combining Eqs. (1) and (2), the dynamic equation of the rigid body can be written as:

[m00J]i{aα}i+{0ω̃Jω}i={FM}i(3)

In coordinate O-xyz, it can be simplified as:

[mb]6×6i{ab}6×1+{Rb}6×1i={Fb}6×1i(4)

where

{ab}={ax,ay,az,αxαyαz}T (5)

{Fb}={Fx,Fy,Fz,MxMyMz}T (6)

{ω}={ωx,ωy,ωz}T (7)

[m]=[m000m000m] (8)

Assume that the inertia of a micro rigid body relative to O’-uvw is

[J*]=[Ju000Jv000Jw].(9)

Then

[J]=[A]T[J*][A](10)

where the direction cosine matrix of O’uvw relative to O-xyz is

[A]=[cos(u,x)cos(u,y)cos(u,z)cos(v,x)cos(v,y)cos(v,z)cos(w,x)cos(w,y)cos(w,z)](11)

The influence of the axial force on bending deformation must be considered in a large deformation of the elastic thin bar. The bending deformation of the elastic thin rod unit is shown in Fig. 4. Considering the influence of the rod axial force on the bending moment, the bending moment M of the rod element is

M=Fviu+Fn(v-Δvi)-Mwi(12)

where v is the deflection of the rod, Fn is the rod axial force, and Fn = − Fui = Fuj.

Figure 4: Bending deformation of the elastic thin rod unit

Given that the displacement of the rod element is small and in the limits of elastic deformation, plastic (unrecoverable) deformation can be ignored. The deflection of the rod v is

v=Δvi+ul(Δvj-Δvi)(13)

where l is the length of the rod unit.

If the displacement of the rod element is relatively small in the local coordinate system O’-uvw and the first derivative of the deflection is small, the differential equation of the rod bending deflection is

d2vdu2=MEIw(14)

where E is the elastic modulus, and Iw is the section moment of inertia.

According to Eqs. (13) and (14) and the integral of u, the boundary condition of the rod element node and the balance equation of the rod element force, the matrix equation of the bending deformation of the elastic rod unit in the plane uO’v can be expressed as follows

{FviMwiFvjMwj}=[12EIwl3+Fnl6EIwl2-12EIwl3-Fnl6EIwl26EIwl24EIwl-6EIwl22EIwl-12EIwl3-Fnl-6EIwl212EIwl3+Fnl-6EIwl26EIwl22EIwl-6EIwl24EIwl]×{ΔviΔθwiΔvjΔθwj}(15)

Similarly, the bending deformation matrix equation of the elastic rod unit in the plane uO0w can be derived.

For the axial tension compression and torsion deformation of the rod unit, the following relationships can be derived

{FuiFuj}=EAl[1-1-11]{ΔuiΔuj} (16)

{MuiMuj}=GIpl[1-1-11]{ΔθuiΔθuj} (17)

where G is the shear elastic modulus, A is the cross-sectional area, and Ip is polar moment of inertia of cross section. According to Eqs. (16) and (17), the stiffness matrix [k*] of the spatial bar element relative to the local coordinate system can be obtained.

[k*]=[EAl00000-EAl00000012EIwl3+Fnl0006EIwl20-12EIwl3-Fnl0006EIwl20012EIvl3+Fnl0-6EIvl2000-12EIvl3-Fnl0-6EIvl20000GIpl00000-GIpl0000-6EIvl204EIvl0006EIvl202EIvl006EIwl20004EIwl0-6EIwl206EIl202EIwl-EAl00000EAl000000-12EIwl3-Fnl000-6EIwl2012EIwl3+Fnl000-6EIwl200-12EIvl3-Fnl06EIvl200012EIvl3+Fnl06EIvl20000-GIpl00000GIpl0000-6EIvl202EIvl0006EIvl204EIvl006EIwl20002EIwl0-6EIwl20004EIwl]

Fig. 5 shows the elastic bar element in the local coordinate system, where i and j are the node positions of the elastic thin rod time t. The coordinate system of node i is taken as the local coordinate system of the unit, and the coordinate plane O’vw is called the cross section of node i, the i~j0 position is the reference position of the rod unit without deformation.

Figure 5: Elastic bar element and local coordinate system

The relationship between the deformation displacement of the elastic rod element and the force of the rod end relative to the overall coordinate system O-xyz can be expressed as

{F}=[k]{Δ}(18)

where the end force and deformation displacement of the element relative to O-xyz are

{F}={Fi,Fj}T,{Fi}={Fxi,Fyi,Fzi,MxiMyiMzi}T,{Fj}={Fxj,Fyj,Fzj,MxjMyjMzj}T,{Δ}={Δi,Δj}T,{Δi}={Δxi,Δyi,Δzi,Δθxi,Δθyi,Δθzi,}T,{Δj}={Δxj,Δyj,Δzj,Δθxj,Δθyj,Δθzj,}T,[k]=[t]T[k*][t](19)

where [k*] is the stiffness matrix of the element relative to O’-uvw, [k] is the stiffness matrix of the element relative to O-xyz; [t] is a coordinate transformation matrix. [t] is the direction cosine matrix of O’uvw relative to O-xyz.

[t]=[[t][t][t][t]];[t]=[cos(u,x)cos(u,y)cos(u,z) cos (v,x) cos (v,y) cos (v,z)cos(w,x)cos(w,y)cos(w,z)]

Taking the coordinate system O’-uvw of node i as the local coordinate system of the unit (see Fig. 4), then we have

{Δi}={0},{Δj}={δe} (20)

{δe}={δx,δy,δz,δφx,δφy,δφz,}T. (21)

Let {δ}={0,δe}T, and then Eq. (10) can be written as

{F}12×1e=[k]12×12e[δ]12×1e(22)

For the dynamic equation of all the nodes, Eq. (4) can be written in a matrix form as

[mb]6(n+1)×6(n+1)i{a}6(n+1)×1+{Rb}6(n+1)×1i={Fb}6(n+1)×1i(23)

where the vacancy element is zero, and {a}6(n+1)×1 is the acceleration of all node degrees of freedom.

The dynamic equation of each elastic rod unit e can be written in a matrix form as

{F}6(n+1)×1e=[k]6(n+1)×6(n+1)e[δ]6(n+1)×1e(24)

{F}6(n+1)×1i in Eq. (13) is regarded as the force required for each micro-segment rigid body, and {F}6(n+1)×1e in Eq. (14) is regarded as the force required for each elastic rod unit. By combing all the node micro-segment rigid bodies and all the elastic rod elements for the FE analysis, the force required for each node can be found

{Q}=i=1n+1{Fb}i+e=1n{F}e(25)

where {Q} is the external force added to all nodes

Substituting Eqs. (23) and (24) from Eq. (25), the new equation can be written as

[mb]{a}+{Rb}+{P}={Q}(26)

where,

[mb]=i=1n+1[mb]6(n+1)×6(n+1)i,{Rb}=i=1n+1{Rb}6(n+1)×1i,and{P}=e=1n[k]6(n+1)×6(n+1)e[δ]6(n+1)×1e.

Eq. (26) is the differential equation describing motions of elastic thin rods for the FE dynamic analysis.

2.3 Simulation

2.3.1 Initial Parameters

The differential equation of the dynamic FE model is solved by a specially developed MATLAB program to simulate the additional twists of yarn in the condensing zone. In the simulation, the ramie fiber was taken as an example where the fiber diameter along the fiber length was assumed to be uniform, the elastic modulus E of the fiber was 2224.6 cN/tex, and the shear elastic modulus G of the fiber was 106.2 cN/tex, the fiber diameter was 0.02 mm, the fiber density was 1.54 g/cm3, and the fiber length in the agglomeration area was 40 mm. The other parameters used in the analysis included the size of the fiber AB division unit of 0.5 mm, the number of fiber AB division units of 200, and the time integral step of 10−6. The gathering roller was made of rubber and supported by a copper sleeve, and thus the friction coefficient μ was set at 0.64. The left lateral angle was 15, the initial position of the fiber in lateral entry compact spinning with pneumatic groove was (X, Y, Z)=(0,0,-1) mm, and the initial velocity of the fiber was 0.

2.3.2 Boundary Constraints

When an edge fiber was connected with central fibers and the lateral pneumatic grooves, the contact was treated as a flexible connection constraint with sufficient stiffness in the simulation. Each node was subjected to the normal elastic connection occurring on the central cylinder, and the wall surface was constrained by the reaction force Fn. Its size was related to the normal spacing ε between them. ε0 was a small parameter that tended to be zero (in this paper, we set ε0=7.5×10-5). When ε was greater than ε0, the connect reaction force was zero. When ε was less than ε0, the contact reaction force was calculated as follows:

Fn=cl(ε0-ε)(27)

where l was the unit length, and c was a sufficiently large stiffness factor (in this paper, c=2×104). To consider the sliding friction between the mutual connect, each node was subjected to frictional force:

Ff=fFn(28)

where f was the coefficient of sliding friction at the connect. The direction of friction was opposite to the direction of the node speed.

To implement the numerical calculation for the FE dynamic differential equations of the elastic thin rods, the implicit difference scheme Euler method was adopted. The calculation process was stable and the calculation results were reliable [23].

2.3.3 Solving Methods

2.3.4 Simulation

Fig. 6 shows the path of a fiber bundle from the entering point A, the gathering roller, and to the exiting point B in the condensing zone. According to the specific value, the fiber motion in the agglomeration area was calculated by using the MATLAB software.

Figure 6: Fiber simulation model

Fig. 7 shows the movement of the fiber bundle. Fig. 7a is the simulation result of intermediate compact spinning with pneumatic grooves, while Fig. 7b is the simulation result of lateral compact spinning with pneumatic grooves. Figs. 7c and 7d are the fiber motion in the YZ plane of the two spinning systems. Unlike fibers in the intermediate compact spinning (Fig. 7c), fibers of the lateral compact spinning with pneumatic grooves had additional twists because fiber 5 crossed fiber 4 (Fig. 7d). Because of the lateral direction of fiber bundle in the agglomeration area, fibers, especially boundary fibers, would move closer to the center. The fiber bundle became more inseparable, which would improve yarn’s properties.

Figure 7: Fiber simulation results. (a): Fiber motion of intermediate compact spinning with pneumatic; (b): fiber motion of lateral compact spinning with pneumatic grooves; (c): YZ plane of fiber motion of intermediate compact spinning with pneumatic; (d): YZ plane of fiber motion of lateral compact spinning with pneumatic [22]

3  Experiments

In order to determine the additional twist, yarns of the lateral compact spinning of pneumatic grooves were spun and compared with the intermediate compact spinning with pneumatic grooves yarns. The raw material was ramie roving. The roving linear density was 4.70 g/10 m, and the moisture regain of the roving was 12.0%. All the spinning was finished by the FZ501 spinning machine. The pressure of pneumatic grooves was −2600 (Pa), the twist was 680 (T/m), the spinning speed was 7000 (r/min), and the count of yarn was 36 Nm.

The tests were performed in the standard climate where the relative humidity was 65%±3% and the temperature was (20±2)C. Before the testing, the specimen conditioned for 48 h in a constant temperature and humidity laboratory. The yarn tension of different counts was adjusted by (0.5±0.1) cN/tex. The tested specimens were 2 kinds of bobbins that had 10 replications per bobbin. The yarn twist tester was the YG331A yarn twist test tester. According to GB/T2543.2-2001 [24], we used a method of the untwist-retwist method to test the yarn of a certain length under the specified tension, and measured the number of rounds when the yarn returned to its initial length after the untwisting and the reversed twisting. The experiment used the counterclockwise direction, the speed was 800 r/min, the length of the sample yarn was 500 mm, and the allowable elongation limit was 4.0 mm. Tab. 1 shows the results of yarn twists.

Table 1: Yarn twist

In the paired t-test on the twists of the two compact spinning systems, the P-value was 0.000 (< 0.05). Therefore, the yarn twists of the lateral compact spinning with pneumatic grooves were significantly different from those of the intermediate compact spinning with pneumatic grooves.

Fig. 8 shows yarn forms of the two different spinning methods. The quality of the yarn made from the lateral compact spinning with pneumatic grooves is better than that from the intermediate compact spinning. The yarn in Fig. 8a appears to be less hairy and tighter due to additional twists than the yarn in Fig. 8b. Therefore, the lateral compact spinning yarns possess better shape and mechanical performance.

Figure 8: Yarn of two spinning systems: (a): yarn of lateral entry compact spinning with pneumatic grooves; (b): yarn of intermediate compact spinning with pneumatic grooves

4  Conclusions

This paper provided a theoretical method for modeling the dynamic motion that simulates fibers driven by mechanical forces in the lateral compact spinning. In the study, we established a fiber FE model to simulate the large deformation process of the elastic thin rod element, and deduced the dynamic equilibrium equation of the spatial elastic thin rod unit in the model. We numerically obtained the movements of fibers under agglomeration and visually displayed additional twists generated in the lateral compact spinning with pneumatic grooves. The experiments demonstrated that the lateral compact spinning could produce smoother and tighter yarns that the intermediate compact spinning due to additional twists. Thus, the fiber FE model may be used to solve the large deformation problem of fiber motion in the three-dimensional space.

Funding Statement: This work is supported in part by China Scholarship Council and National Key R&D Program of China (Grant No. 2017YFB0309100).

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

References

1. Yu, C. W. (2009). Spinning. CN: China Textile Press.
2. Yu, C. W. (2010). Spinning. CN: China Textile Press.
3. Zhou, S. P., Wang, J., & Yang, J. P. (2005). Mechanical analysis of sections of straight sections in the compact field. Donghua University (Nature Science), 31(2), 20-23. [Google Scholar] [CrossRef]
4. Zhou, S. P., Wang, J., & Yang, J. P. (2005). Mechanical analysis of sections of straight sections in the compact field. Donghua University (Nature Science), 31(3), 10-14. [Google Scholar] [CrossRef]
5. Yang, X., Wang, J., & Yang, J. P. (2003). Development status and prospects of concentrated spinning. Donghua University (Nature Science), 329(2), 105-110. [Google Scholar] [CrossRef]
6. Feldman, L. (1996). Theoretical trajectory studies of light bodies in non-uniform two-dimensional flows. Textile Research Journal, 36(9), 809-813. [Google Scholar] [CrossRef]
7. Bangert, L. H., & Sagdeo, P. M. (1997). On fiber alignment using fluid-dynamic forces. Textile Research Journal, 47(12), 773-780. [Google Scholar] [CrossRef]
8. Jeffery, G. B. (1922). The motion of ellipsoid particles immersed in a viscous fluid. Proceedings of The Royal Society A-Mathematical Physical and Engineering Sciences, 102, 161-179. [Google Scholar] [CrossRef]
9. Anezurowski, E., & Mason, S. G. (1967). The kinetics of flowing dispersions: Equilibrium orientations for rods and dices (theoretical). Journal of Colloid and Interface Science, 23(4), 522-532. [Google Scholar] [CrossRef]
10. Chiba, K., & Nakamura, K. (1998). Numerical solution of fiber suspension flow through a complex channel. Journal of Nanonewton Fluid, 78(2–3), 167-185. [Google Scholar] [CrossRef]
11. Chiba, K., & Komatsu, T. (2007). Numerical simulation for orientation of thin disk particles in a Newtonian flow through an L-shape channel. Journal of Textile Engineering, 53(1), 31-35. [Google Scholar] [CrossRef]
12. Skjetne, P., Ross, R. F., & Klingenberg, D. J. (1997). Simulation of single fiber dynamics. Journal of Chemical Physics, 107(6), 2108-2121. [Google Scholar] [CrossRef]
13. Nyland, G. H., Skjetne, P., Mikkelsen, A., & Elgsaeter, A. (1996). Brownian dynamics simulation of needle chains. Journal of Chemical Physics, 105(3), 1198-1207. [Google Scholar] [CrossRef]
14. Cheng, D. W. (1987). Discrimination of the plane motion of curved fiber. Donghua University, 13(6), 73-77. [Google Scholar]
15. Cheng, D. W. (1998). Suspension dynamics of viscoelastic fibers in a uniform deformation flow field. Donghua University, 14(6), 89-96. [Google Scholar]
16. Yamamoto, S. M. (1993). A method for dynamic simulation of rigid and flexible fibers in a flow field. Journal of Chemical Physics, 98(1), 644-650. [Google Scholar] [CrossRef]
17. Yamamoto, S. M. (1994). The viscosity of dilute suspensions of rod like particles: A numerical simulation method. Journal of Chemical Physics, 100(4), 3317-3324. [Google Scholar] [CrossRef]
18. Zeng, Y. C., & Yu, C. W. (2003). A bead-elastic rod model for dynamic simulation of fibers in high-speed air flow. International Journal of Nonlinear Science and Numerical Simulation, 4(2), 201-202. [Google Scholar] [CrossRef]
19. Zeng, Y. C. (2003). Research and application of fiber movement in nozzle high-speed airflow, (Ph.D. Thesis). Donghua University, China.
20. Wang, Y. (2010). Study on aggregation mechanism and yarn performance of grid-type compact aggregate spinning, (Ph.D. Thesis). Donghua University, China.
21. Liu, Y. Z. (2006). Nonlinear mechanics of elastic thin rods. CN: Tsinghua University Press.
22. Lyu, J. D., Cheng, L. D., Xu, B. G., Hua, Z. H. (2021). Numerical simulation of fiber motion of lateral compact spinning with pneumatic groove. Autex Research Journal, (In press).
23. Jiang, J. F., Hu, L. J., Tang, J. (2004). Numerical analysis and matlab experiment. CN: Science Press.
24. Huo, S. H. (2001). GB/T 2543.2-2001. Textiles-determination of twist in yarn-part2: Untwist-retwist method. CN: China Standards Press.