## Abstract

In this paper, a new method to simulate the structure and loop deformation behavior of double-bar reflex-lapping warp-knitted fabrics based on the structural characteristics is proposed. A simplified mass-spring model was built in which loops knitted by filaments were considered as particles with the uniform mass distribution connected by structure springs for overlaps and shear springs for underlaps. Deformation forces and direction on particles were analyzed to describe the displacement and deformation behavior of particles. A loop model with eight control points was established, and the relationship between control points and particles was studied combining the quadratic Bezier curves. The deformation simulation was implemented by a simulator program with C# and JavaScript via web technology on Visual Studio 2015. The stereoscopic sense of filaments was realized by changing the direction and intensity of the light. The results show that the fabric deformation and the loop shape can be accurately achieve using the simplified mass-spring model compared with the real sample.

## 1 Introduction

The double-bar reflex-lapping warp-knitted fabric is one of the most basic kinds in warp-knitted fabrics, which needs repeated proofing during development by lack of strong visibility and intuition of actual results. In recent years, the efficiency and accuracy of warp-knitted fabric design have been improved successfully because of the application of computer-aided design software. However, three-dimension simulation including structure deformation and single loop shape change for warp-knitted fabrics remains an ongoing research problem. Structural deformation is a basic attribute of warp-knitted fabric, and it is also the core of the simulation. The key influencing factors on deformation, such as density, chain notation, and lapping regulars, will result in limited production and massive resource consumption in continuous sampling. The accurate judgment of the appearance leads to the real state of fabric structure, and the pattern of the fabrics can be well displayed so that the design draft can be revised repeatedly until the requirements are met. Thus, many methods have been proposed in this regard.

In the study of fabric deformation, the mass-spring model was one of the most widely used and typical models that was developed by Pierce [1] in 1947. In this model, the mass of the fabric was assumed to be particles of uniform mass, and the connection of the two loops was assumed to be an elastic spring. Provot et al. [2-3] put forward an improved mass-spring model in which the fabric was regarded as a collection of several particles, and the mutual relationship among particles was attributed to the spring action. It was widely used in the dynamic simulation of fabrics [4, 5, 6]. Another mass-spring model with modulus and tensile data of each region in a piece of fabric was proposed by Ji et al. [7] to simulate the dynamic draping behavior of woven and knitted fabrics.

Choi and Ko [8] proposed a stable and practical method to deal with the post-buckling instability without any actual damping. Based on this method and the classic mass-spring model, Hu et al. [6] simplified the physical model for cloth by simplifying the calculation method of spring bending resistance and the structural force. The shear force of the spring was also combined. Liu and Song [9] also reduced some shear springs in the same direction. But it did not apply to the warp-knitted fabric with underlaps in multiple directions. Instead of explicit calculation of internal force and particle velocity, a constrained mass-spring model with reverse dynamics constraining for the internal force simulation was established by Wang et al. [5] All of the above methods reduced the amount of calculation and increase the simulation speed.

In the research of Mozafary [10-11], a new mass-spring model defined as loop mesh was applied for every single loop with eight mass particles of weft-knitted fabrics. A single loop was divided into seven segments, by which displacement and shape of each part were controlled separately on account of various stress conditions so that selvage curling of weft-knitted fabric could be well simulated. In order to model the fabric thickness, a mass-spring mesh with a layered structure [12] was proposed in 2012. A mass-spring model [13] was used to describe the deformation of the warp-knitted loop in detail, in which each segment of the loop was modeled as a spring with a definite constant. This method was used for the prediction of the elastic modulus of composites instead of the fabric deformation. The fabric deformation and even the loop shape could be simulated accurately through the preceding ways, but there also existed a problem of big calculation.

Sha [14, 15, 16] not only established the spring-mass model of weft-knitted fabric but also studied the loop details. In this model, a rectangle formed by every four adjacent particles represented a loop, and the loop shape was controlled by these four particles. Similarly, a mass-spring model with mesh was proposed by Peng et al. [17], where two grids corresponded to one loop. A loop was divided into two parts, upper and lower, which was controlled by six particles, and the force analysis and deformation calculation of the whole loop was implemented. However, in their research, the length of shear springs was fixed, which was not suitable for warp-knitted structures. Moreover, Arumugam et al. [18] focused on the relationship between shear performance and deformation of weft-knitted spacer fabrics. Although fabric deformation was involved, the details of loop deformation are ignored.

The relationship between yarn tension and loop deviation based on the structure characteristic such as underlaps was focused [19-20]. Nevertheless, only the deformation influence for the current loop was considered in this method, ignoring the chain reaction for neighboring loops. Moreover, although the transverse displacement of loops was well studied, no consideration of displacement was put in other directions.

A deformation simulation model was based on the mass-spring model for warp-knitted fabrics, and attention was only paid to Jacquard structures [21]. With the further study of warp-knitted jacquard fabric deformation, springs for ground pillars, jacquard inlays, and pattern inlays were generated separately [22]. However, both models were developed for Jacquard fabrics only and thus not suited for broad application. In addition, the problem of super-elasticity in knitted fabric simulation was effectively solved [4,23-24].

Moreover, the explicit Euler integration method had been widely used in numerical integration in recent years [25]. One of the most significant advances was implicit integration. However, implicit integration was a linear solution system, which was not easy to calculate. Therefore, some effective approximate methods were proposed [26-27]. The velocity-Verlet method was one of them that took advantage of force ranges and integral methods. Stable movement with a large time step was obtained [22].

Most deformation simulators were developed for weft-knitted fabrics. Different from weft-knitted fabrics, yarns of each bar in warp-knitted fabrics were fed as layers, which resulted in one needle being affected by one or more loops. Using yarns in the proper specification, fabrics knitted by the high-speed tricot warp-knitting machine with two ground bars have the following advantages: this kind of machine has large working width, which leads to a wide fabric widely applied in industry and clothing. Moreover, knitting the same type of loops, which meant that fabrics knitted only by open laps or only by close laps, achieved optimal performance of machines with high efficiency.

Motivated by the above requirements, this paper constructs a mass-spring model for ordinary warp-knitted fabric, especially for double-bar reflex-lapping warp-knitted fabrics. The characteristics of ordinary warp-knitted fabric are considered, and the force is updated several times in the order of bars for a mass-spring model related to one fabric. Although force on the current course is the main cause of displacement, loops on the next course also influence the shape and position, especially for the loop arc because of the intermeshing relationship. The shape change of each loop caused by a nonuniform force is also studied. Additionally, the deformation range changes accordingly with the process parameters such as density. The pattern effect is simulated in real-time once the parameters change. Thus, the fabric can be modified before it is put into production to save cost and time. Furthermore, this simulator relies on web service, which supports the 3D rendering technology of WebGL to implement the display of fabric simulation. Another advantage is that fabrics can be designed and simulated as long as users have internet access at the place, which is very convenient with a break of space and geographical restrictions.

## 2 Mass-spring modeling

Generally speaking, the deformation of fabrics knitted by two bars results from a superposition of the force of multiple loops rather than one single loop. Ideally, the warp-knitted fabric is assumed as flat with uniform quality, which can be represented by a group of evenly distributed particles linked by massless springs in a relaxed state. Thus, a simplified mass-spring model for double-bar reflex-lapping warp-knitted fabrics is presented in this paper. For an ideal warp-knitted flat fabric, three-dimensional bending deformation controlled by bending spring with low elastic coefficient can be ignored, which is removed in this improved mass-spring model.

### 2.1 Simplified particle system

It is assumed that the warp-knitted fabric is composed of *m* wales and *n* courses. Each particle is positioned at the bottom of the loop, that is, the intersection position of the current and the previous current loops. Thus, the particle composition of the fabric is represented with the specification of *m* * (*n*+1). It is revealed from Figure 1 that the particles are arranged from the right-bottom to left-top because yarns on the machine are equipped from right to left, and courses in the fabric are knitted from bottom to top. The finished wale density in 1 cm is denoted as wpc, and the density of the finished course is cpc. *Wg* is defined as the length between adjacent wales determined by finished horizontal density (wpc), while *hg* is the course distance determined by finished longitudinal density (cpc), and they are respectively calculated as shown in equations (1) and (2).

### 2.2 Mass-spring model for overlaps

Double-bar reflex-lapping warp-knitted fabrics are formed by interconnected loops whose loop arc of the current course is connected with the loop root of the next course. Therefore, the position of a loop is determined by two particles *P(m, n)* controlling the loop root and *P(m, n+1)* controlling the loop arc. Springs linking these two particles are referred to as massless structural springs in the longitudinal direction as shown in Figure 2. However, there is no structural spring in the transverse direction because neighboring wales are connected by underlaps rather than overlaps.

### 2.3 Mass-spring model for underlaps

Underlaps of double-bar reflex-lapping warp-knitted fabrics are used for the connection of overlaps between adjacent courses whose length depends on both loop height and the number of wales they across. Different lengths of underlaps are subjected to different forces, and the corresponding shear springs are shown in Figure 3. It is also revealed in Figure 3 that one particle with a loop is acted upon by forces of two underlaps in one single bar.

### 2.4 Structure deformation calculation

The warp-knitted fabric is simulated when laid on a flat platform so there is no consideration of outer forces, such as support forces, gravity, and so on. Only the internal forces between loops covering forces of structural springs, shear springs, and damping force are analyzed when studying factors for deformation. The forces can be calculated by the following equations [15,17]:

Where
*P _{m,n}* and

*P*

_{m,n+1},

*H*

_{m,n−m,n+1}is the length of this structural spring under forces, and

*k*

_{st}is the structural spring coefficient.

*P*,

_{m,n}*P*

_{m+d,n+1}·

*P*

_{m+d,n+1}is the particle connected by the underlap starting from the loop on the particle

*P*, and

_{m,n}*dis*is the number of wales crossed by this particular underlap.

*l*is the shear spring length under force, and

_{shear}*k*is the coefficient of the shear spring.

_{sh}*P*·

_{m,n}*P*(

_{m,n}*t*) is the position at current time t, and

*P*(

_{m,n}*t*− Δ

*t*) is the position at the last moment (

*t*− Δ

*t*).

*t*is the time step, and

*k*is damping coefficient. Comprehensively,

_{d}*F*is the sum of all internal forces the particle

_{m, n}*P*received.

_{m, n}The Verlet integration method is used to determine the particle position without calculation of particle velocity of the particle^{28}. As shown in equation (7), the position of the particle at the future time (*t* − Δ*t*) can be derived from the position at the last and the current moment, which leads to a more stable system. In equation (7), *M _{m,n}* is the particle mass.

During the movement, the kinetic energy of particles will go down continually because of the damping force and displacement will also decrease with it. The displacement threshold of particles is set as *L _{th}* in order to avoid overstretching of the spring and ensure the rationality of the deformation. Because excessive stretching of springs leads to yarn breakage, and excessive shift causes interpenetration between particles. If the particle displacement

*L*≤

*L*, it means that particles are in a stable state and the deformation is completed.

_{th}## 3 Geometric loop model

### 3.1 Bezier curve model for loops

The Bezier curve has a powerful function for drawing smooth curves by any given point, and the loop path can be limited within the appropriate ranges to generate a more reasonable loop shape. In this paper, the quadratic Bezier curve with three control points shown in Figure 4 is selected to calculate the coordinate due to less amount of calculation. The equation of the quadratic Bezier curve is shown as equation (8). *Q(t)* represents the dynamic point at the moment *t* on the curve, and *Q _{1}-Q_{3}* are the control points of the quadratic Bezier curve.

Several control points are required to determine the basic shape of loops in the warp-knitted fabric in this paper as shown in Figure 5. The green dots (*A _{0}-A_{7}*) represent control points, while the blue ones (

*M*-

_{1}*M*

_{7}) are the midpoints of each segment connected by two adjacent control points.

*A*and

_{1}-A_{2}*A*control two-loop limbs, and

_{5}-A_{6}*A*are for the loop arc.

_{2}-A_{3}-A_{4}-A_{5}*A*-

_{0}*A*and

_{1}*A*-

_{6}*A*are transition articulation curves linking the loops in adjacent courses. This figure also reveals that

_{7}*h*is the height of the loop limb, which is the relaxed vertical distance of the particles in the warp direction.

*w*is the width of the loop,

*b*is the height of the loop arc, and

*d*is the thick of the loop model. The loop model viewed from the front and the side are shown in Figure 5 (b) and (c) separately. Polylines connecting green dots are auxiliary lines to illustrate the position and size relationships between each controlling point. The frame of the actual stitch curves can be derived from points on these polylines. The overlap in warp-knitted fabric includes the loop limb and loop arc, which can be denoted as a smooth path composed of multiple Bezier curves as shown in Figure 6. The Bezier curves are generated automatically by the built-in geometry engine in the Three.js once the midpoints (M1-M7) are calculated.

The 3D loop model proposed in this paper consists of eight segments. Figure 7 is a plane graph of loops in two courses. *A _{i}* and

*B*are the control points, and

_{i}*M*is the midpoint of each segment. The underlaps are controlled by the last two points of the current loop (

_{i}*A*,

_{6}*A*) and the start two points of the loop in the next course (

_{7}*B*).

_{0}, B_{1}*M*-

_{1}*A*-

_{1}*M*is the first segment, and

_{2}*M*-

_{1}*A*-

_{2}*M*is the second segment. And so on, for each of the seven segments, a complete loop model is composed.

_{3}Take the second segment *M _{2}*-

*A*-

_{2}*M*for example shown in Figure 6 (b). Assuming that the coordinate of

_{3}*M*is (

_{2}*x*,

_{1}*y*,

_{1}*z*),

_{1}*A*is (

_{2}*x*,

_{2}*y*,

_{2}*z*),

_{2}*M*is (

_{3}*x*,

_{3}*y*,

_{3}*z*), and

_{3}*N*on the curve is (

*x*,

_{n}*y*,

_{n}*z*). In particular,

_{n}*y*coordinate of point

*N*is the same as that of point

*A*. Thus,

_{2}*M*-

_{2}*A*-

_{2}*M*is a Bezier curve through

_{3}*N*so that the coordinate of

*N*can be used in the quadratic Bezier curve equation. Calculation equations are as follows:

Solving *t _{n}* from equation (8) gives that

Substituting the value of *t _{n}* into the quadratic Bezier equation, values of

*x*and

_{n}*z*can be obtained as:

_{n}

Other curves are calculated in a similar way. And space relationship of loops without shape deformation is shown in Figure 8. It is a common fabric knitted by two guide bars. The red yarns are threaded in the needle on the front bar (GB1), whose overlaps are on the technical face, and the underlaps are on the technical back. The green stitches are knitted by GB2, which are wrapped in the middle of the fabric.

### 3.2 Loop deformation with control points

The loop model established above is ideal in a relaxed state. When subjected to a force, the loop shape will deform to a certain extent. It is assumed that the force only acts in the *x* and *y* directions without changing the *z*-axis coordinate. As shown in Figure 9(a), the basic position of a loop is determined by particles *P _{0}* (

*m*,

*n*) and

*P*(

_{1}*m*,

*n*+1), which represent one particle on the wale

*m*and the course

*n*and the other particle on the wale

*m*and course

*n*+1. The warp-knitted fabric has a structure with loop interlocking, and the loop arc of the current course is connected to the root of the next course. Therefore, the position of the loop root is controlled by particle

*P*(

_{0}*m*,

*n*), while the loop arc depends on

*P*(

_{1}*m*,

*n*+1). Figure 9(b) shows that the control points surrounding particles will move with the change of particles’ coordinates under forces. It is necessary to establish the relationship between the particles and the control points in order to better describe the shape change of the loop.

In Figure 9(a) and (b), *A _{0}*-

*A*are control points in the relaxed state, while

_{7}*A*’-

_{0}*A*’ are in one of the states under force.

_{7}*w*-

_{1}*w*are relaxed distances between control points and the centerline.

_{4}*h*and

*b*are the relaxed height of the loop limb and the loop arc. When the loop is deformed, an inclination angle

*α*is generated, which can be calculated by the coordinates of

*P*and

_{0}*P*as in equations (15)-(16). The length of

_{1}*P*-

_{0}*P*and

_{1}*P*remain unchanged as

_{1}M’*h*and

*b*.

*Δy*and

_{1}*Δy*are the vertical height of the loop limb and loop arc after deformation.

_{2}*Δx*is the relative displacements to

_{1}*P*horizontally, while

_{1}*Δx*is that to

_{2}*P*. The values can be obtained by equations (17)–(20).

_{0}

Values of the *z* axis of control points remain stable since only the position change in the X-Y plane is considered in this paper. But the need to pay attention to it is that the start point *A _{0}* and the endpoint

*A*change dynamically. Referring to Figure 9 (c) for specific definitions, the green segments are the start of the loop, and the yellow ones are the segments linking to the next loop. The direction of them is affected by the loop types and position relationship of loops connected to it. The coordinates of control points are shown in Table 1.

_{7}Point |
x |
y |
z |
---|---|---|---|

M’ | P_{1}(x)+b *cosα |
P_{1}(y)+b *sinα |
A_{3}(z) |

*A_{0}’ |
P_{0}(x) - w_{1} |
P_{0}(y)-w_{0} |
A_{0}(z) |

A_{1}’ |
P_{0}(x) - w_{1} |
P_{0}(y) |
A_{1}(z) |

A_{2}’ |
P_{1}(x) - w_{2} |
P_{1}(y) |
A_{2}(z) |

A_{3}’ |
P_{1}(x)+b *cosα-w_{3} |
P_{1}(y)+b *sinα |
A_{3}(z) |

A_{4}’ |
P_{1}(x)+b *cosα+w_{3} |
P_{1}(y)+b *sinα |
A_{4}(z) |

A_{5}’ |
P_{1}(x)+w_{2} |
P_{1}(y) |
A_{5}(z) |

A_{6}’ |
P_{0}(x)+w_{1} |
P_{0}(y) |
A_{6}(z) |

*A_{7}’ |
P_{0}(x)+(w +_{0}w)_{1} |
P_{0}(y) |
A_{7}(z) |

Note that

*A*and_{0}*A*with ‘*’ mean that they change dynamically._{7}

## 4 Results and discussion

The simulator proposed in this paper is designed by Visual Studio 2015. Using C#, a kind of computer language, the mass-spring model and warp-knitted loop models are established, and key parameters, such as the force, displacement of particles, and coordinate of loop control points, are calculated. Due to the control points of loops computed by two vertical neighboring particles, each loop can be deformed along with the displacement of these related particles. Once getting the point coordinates, smooth tubes are generated along a specific path with three.js technology. Ambient light and directional light are used to cause shadings so that the stereo feeling of yarns can be reflected, realizing the simulation of double-bar reflex-lapping warp-knitted fabrics with chemical filaments. The flow chart of deformation simulation is as Figure 10 shows.

In order to compare the influence of different chain notations on fabric structure deformation, two pieces of double-bar reflex-lapping fabrics using the same yarns are selected. The technical parameters are shown in Table 2. The particle displacement results are shown in Figure 11. Due to the lateral displacement of the warp-knitting bars, some laps of the pattern move more than one wales. Particularly on the leftmost yarn moves beyond the particle model, loops out of range need to be constrained on particles at the edge to ensure the model in regular.

No. |
Chain notation |
Yarn specification |
Threading |
---|---|---|---|

1 | GB1:1-0/1-2//; GB2:1-2/1-0//; |
A: nylon in 22dtex/1F, pink B: nylon in 22dtex/1F, blue |
GB1:1A; GB2:1B. |

2 | GB1:1A, 1*; GB2:2B. |

When the particles are in equilibrium as shown in Figure 11 (a) and (b), particles remain at the original position without offset, and loops are in an ideal upright state. However, Figure 11 (c) and (d) illustrate if there are some needles without yarn threading, which means some particles are in the forces of two loops and some are in the forces of only one loop, there will be an imbalance in the force on particles, and displacement will cause the deformation of the fabric. In Figure 11 (c), the red mesh represents the initial position without deformation. Based on this, the irregular blue mesh with white particles shows the deviation results.

Another example is taken to illustrate the effect of the parameter change on fabric deformation and the real-time simulation. A pattern repeat contains six courses and two wales knitted on a high-speed tricot warp-knitting machine HKS2 (E28) produced by KARL MAYER (Germany). The basic parameters are listed in Table 3.

Pattern width |
2 wales |
Pattern height |
4 courses |
---|---|---|---|

Finished horizontal density |
9 wpc |
Finished longitudinal density |
11.2 cpc |

Threading |
GB1:1*,1A; GB2:1*,1A; | ||

Run-in values |
GB1:1800mm/rack; GB2:1800mm/rack; | ||

Chain notation |
GB1:2-3/2-1/2-3/1-0/1-2/1-0//; GB2:1-0/1-2/1-0/2-3/2-1/2-3//; |
||

Materials |
A:76dtex/24F, polyester; |

After the technical parameters are set, the simulation is conducted according to the flow chart. Four figures show the fabrics in the view dimensions of 234.96 mm and 83.75 mm. The real sample and the simulation result of the one-one ratio are shown in Figure 12 (a) and (b). The fabric contains about 13 diamonds in the course direction and about 6.5 diamonds in the wale direction. To test the real-time simulation, one of the technical parameters is changed to resimulation and compared with real samples. Figure 12 (c) shows the simulation with finished horizontal density in 6 wpc and finished longitudinal density in 11.2 cpc, which contains about 8.3 diamonds in the course direction and about 6.5 diamonds in the wale direction. While Figure 12 (d) shows the finished density in 9 wpc and 11.2 cpc with about 20 diamonds in the course direction and about 6.5 diamonds in the wale direction. The final simulated deformation effect is displayed once the parameter is changed.

Figure 13 are the simulation effects of the same part with actual knitting parameters. The data in Table 4 are the average value of 50 measurements for different parts of the real sample, and the measurement values of simulation are generated automatically. The specific measuring position is marked in Figure 13. P_{1}-P_{4} are four vertices of the diamond, while <α and <β are two adjacent inside vertex angles of the diamond. As revealed from Table 4, the deviations between the simulation and the real sample are small, which indicate the deformation simulator has high accuracy.

Tested pairs |
Simulation |
Sample |
Deviation |
---|---|---|---|

P_{1}-P_{3} |
12.5 mm | 12.9 mm | 3.10% |

P_{2}-P_{4} |
17.8 mm | 18.1 mm | 1.58% |

<α |
118.127° | 118.902° | 0.65% |

<β |
64.514° | 64.973° | 0.70% |

## 5 Conclusion

A deformation simulation method with a simplified mass-spring model for double-bar reflex-lapping warp-knitted fabrics is proposed in this paper. Structural springs and shear springs are, respectively, denoted as yarns of loop limbs and underlaps based on the warp-knitted technical characteristics. The forces on particles are analyzed, and the particles’ motion state and deformations are calculated based on Verlet integral. Besides, a simple model for hyperelastic detection of particles is established in order to improve the accuracy and rationality of particle movement. The loop path is generated by quadratic Bezier curves, and dynamic models for loops and inlays are built. The control points of each loop in a deformation state are calculated according to the relationship of the particles and control points. The validity and correctness of the spring particle model to simulate the deformation behavior are verified by calculating the deformation ratio analysis. C# and JavaScript via web technology were used to implement the deformation simulation and display of warp-knitted fabrics. Furthermore, the simulation results change with different parameters, which is helpful for the designer to modify the parameters before knitting. It also has obviously exhibited time and cost savings and greater practicability.

## ACKNOWLEDGMENTS

This work was supported by the National Natural Science Foundation of China (61772238 & 61902150) and Taishan Industry Leading Talents (tscy20180224).

## References

[1] Peirce, F. T. (1947). Geometrical principles applicable to the design of functional fabrics. Textile Research Journal, 17(3), 123–147.10.1177/004051754701700301Search in Google Scholar

[2] Provot, X. (1995). Deformation constraints in a Mass-Spring Model to describe rigid cloth behavior. Proceedings of the 1995 Graphic Interface Conference, 147–155.Search in Google Scholar

[3] Provot, X. (1997). Collision and self-collision handing in cloth model dedicated to design garment. Proceedings of the 1997 Graphics Interface Conference, 177–189.10.1007/978-3-7091-6874-5_13Search in Google Scholar

[4] Hui, W. H., Guan Z. X. (2007). Research on cloth modeling simulation based on mass spring model. Computer Engineering and Design, (3), 664–667.Search in Google Scholar

[5] Wang, J. C., Hu, X.R., Wang H. (2008). Simulation technology of cloth modeling based on mass spring model. Journal of Wuhan University of Science and Engineering, (9), 5–8.Search in Google Scholar

[6] Hu, J. F., Yan, P. M., Zhou, J. W., et al. (2008). Real-time cloth simulation based on improved mass-spring model. Computer Simulation, (7), 184–187.Search in Google Scholar

[7] Ji, F., Li, R. Q., Qiu, Y. P. (2007). Simulate the Dynamic Draping Behavior of Woven and Knitted Fabrics. Journal of Industrial Textiles, 35(3), 201–215.10.1177/1528083706055753Search in Google Scholar

[8] Choi, K. J., Ko, H. S. (2002). Stable but responsive cloth. ACM Transactions on Graphics, 21(3), 604–611.10.1145/566654.566624Search in Google Scholar

[9] Liu, L. X., Song, Q. (2011). Cloth-object deformation simulation based on a simplified mass-spring model. Computer Simulation, (5), 406–409.Search in Google Scholar

[10] Mozafary, V., Payvandy, P. (2017). A novel method based on loop shape for simulating knitted fabric using mass spring model. Fibers and Polymers, 18(3), 533–541.10.1007/s12221-017-1066-0Search in Google Scholar

[11] Mozafary, V., Payvandy, P., Rezaeian, M. (2018). A novel approach for simulation of curling behavior of knitted fabric based on mass spring model. Journal of the Textile Institute, 109(12), 1620–1641.10.1080/00405000.2018.1453635Search in Google Scholar

[12] Güdükbay, U., Bayraktar, S., Koca, Ç., et al. (2014). Particle-based simulation of the interaction between fluid and knitwear. Signal Image and Video Processing, 8(3), 415–422.10.1007/s11760-012-0308-2Search in Google Scholar

[13] Dabiryan, H., Ashouri, M. H. (2020). Mesoscale model for elastic modulus of warp-knitted fabric reinforced composites based on spring model. Journal of Industrial Textiles, 50(1), 46–69.10.1177/1528083718821873Search in Google Scholar

[14] Sha, S., Jiang, G. M., Ma, P. B., (2015). et al. 3-D dynamic behaviors simulation of weft knitted fabric based on particle system. Fibers and Polymers, 16(8), 1812–1817.10.1007/s12221-015-5254-5Search in Google Scholar

[15] Sha, S., Ma, P. B., Chapman, L. P., et al. (2017). Three-dimensional modeling and simulation of deformation behavior of fancy weft knitted stitch fabric. Textile Research Journal, 87(14), 1742–1751.10.1177/0040517516658518Search in Google Scholar

[16] Sha, S., Jiang, G. M., Chapman, L. P., et al. (2017). Fast penetration resolving for weft knitted fabric based on collision detection. Journal of Engineered Fibers and Fabrics, 12(1), 50–58.10.1177/155892501701200106Search in Google Scholar

[17] Peng, J. J., Jiang, G. M., Xia, F. L., et al. (2020). Deformation and geometric modeling in three-dimensional simulation of fancy weft-knitted fabric. Textile Research Journal, 90(13–14), 1527–1536.10.1177/0040517519894749Search in Google Scholar

[18] Arumugam, V., Mishra, R., Militky, J., et al. (2016). In-plane shear behavior of 3D spacer knitted fabrics. Journal of Industrial Textiles, 46(3), 868–886.10.1177/1528083715601509Search in Google Scholar

[19] Zhang, L. Z., Jiang, G. M. (2010). Loop offset of warp-knitted structures in 3-D simulation. Journal of Textile Research, 31(3), 64–67.Search in Google Scholar

[20] Guo, C. X. (2018). The study of structure and deformation of warp-knitted mesh fabric with less guide bars. Wuhan Textile University, Wuhan, China.Search in Google Scholar

[21] Deng, W., Wang, T. T., Deng, Z. M. (2011). Mass-spring models used in jacquard fabric. Journal of Wuhan University of Science and Engineering, 24(3), 74–78.Search in Google Scholar

[22] Li, X. X., Zhang, A. J., Ma, P. B., et al. (2017). Structural deformation behavior of Jacquardtronic lace based on the mass-spring model. Textile Research Journal, 87(10), 1242–1250.10.1177/0040517516651103Search in Google Scholar

[23] Chen, Y., Wu, M. Z., Hu, X. R. (2013). An improved mass-spring cloth model and simulation. Journal of South China Normal University(Natural Science Edition), 45(5), 43–46.Search in Google Scholar

[24] Mozatary, V., Payvandy, P. (2016). Study and comparison techniques in fabric simulation using mass spring model. International Journal of Clothing Science and Technology, 28(5), 634–689.10.1108/IJCST-03-2015-0033Search in Google Scholar

[25] Baraff, D., Witkin, A. (1998). Large Steps in Cloth Simulation. Proceedings of Siggraph, Orlando, 43–54.10.1145/280814.280821Search in Google Scholar

[26] Ng, H. N., Grimsdale, R. L. (1995). Image Analysis Applications and Computer Graphics. Lecture Notes in Computer Science, Springer-verlag, Berlin, 1024, 124–331.10.1007/3-540-60697-1_95Search in Google Scholar

[27] Kubiak, B., Pietroni, N., Ganovelli, F., et al. (2007). A Robust Method for Real-time Thread Simulation. Proceedings of the ACM VRST, California, 85–88.10.1145/1315184.1315198Search in Google Scholar

[28] Mongus, D., Repnik, B., Mernik, M., et al. (2012). A hybrid evolutionary algorithm for tuning a cloth-simulation model. Applied Soft Computing, 12(8): 266–273.10.1016/j.asoc.2011.08.047Search in Google Scholar

**Published Online:**2021-12-22

© 2021 Haisang Liu et al., published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.