Modelling of Auxetic Woven Structures for Composite Reinforcement

The current research is focused on the design and development of auxetic woven structures. Finite element analysis based on computational modeling and prediction of axial strain as well as Poisson’s ratio was carried out. Further, an analytical model was used to calculate the same parameters by a foldable zig-zag geometry. In the analytical model, Poisson’s ratio is based on the crimp percentage, bending modulus, yarn spacing, and coefficient of friction. In this yarn, properties and fabric parameters were also considered. Experimental samples were evaluated for the actual performance of the defined auxetic material. Auxetic fabric was developed with foldable strips created in a zig-zag way in the vertical (warp) direction. It is based on the principle that when the fabric is stretched, the unfolding of the folds takes place, leading to an increase in transverse dimensions. Both the analytical and computational models gave close predictions to the experimental results. The fabric with foldable strips created in a zig-zag way in the vertical (warp) direction produced negative Poisson’s ratio (NPR), up to 8.7% of axial strain, and a maximum Poisson’s ratio of −0.41 produced at an axial strain of around 1%. The error percentage in the analytical model was 37.14% for the experimental results. The computational results also predict the Poisson’s ratio with an error percentage of 22.26%. Such predictions are useful for estimating the performance of auxetic woven structures in composite reinforcement. The auxetic structure exhibits remarkable stress-strain behavior in the longitudinal as well as transverse directions. This performance is useful for energy absorption in composite reinforcement.


Introduction
Auxetic structures belong to a class of extraordinary materials that become thicker in the perpendicular direction when it is stretched longitudinally. In other words, these materials exhibit negative Poisson's ratio (NPR) and therefore show perceptible improvement over conventional materials in toughness, resilience, shear resistance, and acoustic properties mainly due to their special structure and associated deformation mechanics [1][2][3][4][5]. Research has shown that the auxeticity of a structure depends on its geometrical arrangement and the way it deforms under the action of external force. Zero Auxeticity describes the state when the material is stretched in the longitudinal direction and there is no change in the transverse direction, that is, neither contraction nor expansion. In various studies in the literature, it has been reported that auxeticity can be introduced on all scales making it scale-independent. Auxetics can be broadly categorized into two categories, natural and synthetic, based on their origin [6][7][8][9][10][11][12][13]. Certain natural auxetic materials have been reported like the cancellous bone in humans [14], cow teat skin, cat skin, and salamander

Analytical Model
To develop the analytical model, a foldable zig-zag geometry as shown in Figure 1 was considered for the formation of a woven construction. The basic unit cell design that is to be geometrically analyzed is shown in Figure 2.
Textiles 2021, 1, FOR PEER REVIEW 3 load-bearing capacity and energy absorption in a composite reinforced with auxetic woven structures.

Analytical Model
To develop the analytical model, a foldable zig-zag geometry as shown in Figure 1 was considered for the formation of a woven construction. The basic unit cell design that is to be geometrically analyzed is shown in Figure 2. The combination of plain weave and float is used to generate the desired auxetic design as shown in Figure 2. In this design, when the structure is pulled along the direction of the warp, it tries to open the fold and expands the fabric in the transversal direction. It is assumed that two different types of yarns are used in the weft with elasticity, El and Es, and elasticity of warp yarn be Ep.  ven structures.

Analytical Model
To develop the analytical model, a foldable zig-zag geometry as shown in Figure 1 was considered for the formation of a woven construction. The basic unit cell design that is to be geometrically analyzed is shown in Figure 2. The combination of plain weave and float is used to generate the desired auxetic design as shown in Figure 2. In this design, when the structure is pulled along the direction of the warp, it tries to open the fold and expands the fabric in the transversal direction. It is assumed that two different types of yarns are used in the weft with elasticity, El and Es, and elasticity of warp yarn be Ep.  The combination of plain weave and float is used to generate the desired auxetic design as shown in Figure 2. In this design, when the structure is pulled along the direction of the warp, it tries to open the fold and expands the fabric in the transversal direction. It is assumed that two different types of yarns are used in the weft with elasticity, E l and E s , and elasticity of warp yarn be E p .
In the unit cell of the design, there are 18 warp and 17 weft yarns. Considering each pick, there are two portions of different weave patterns, for example, the first weft is interlaced with nine warp yarns in a plain weave pattern, and the rest is with 1-up-4-down Textiles 2022, 2 4 twill (float). The weft has two yarns that are elastic and stiffer yarn. The elastic yarn will be responsible for the extension in the floating part initially. The pull-out force of yarn from a plain-woven fabric is much more compared to that of 1 up 4 down twill weaves, as the interlacements in the case of plain weave are more. The stiffer weft yarn would require force to overcome the friction resistance. The stress at both sides of the fabric is equal and varying as the rate of displacement of the ends of the fabric is constant. Solving for the first pick, there are two segments-plain weave and float-the plain weave is toward the left side of the line of center. The force on both sides (connecting point) of the portion is different due to the non-symmetry of the structure. The force acting on the right side of the plain weave part would be more than on the left side. This would be balanced by the frictional force at each contact, which would accumulate in one direction due to differences in the end forces. The frictional force at each contact point in plain weave and float is given by Pan [50] as shown in Equations (1)- (3).
where P is the compression force acting on the point of intersection of warp and weft yarns, θ is the weave angle, and σ is the applied stress, µ is coefficient of friction, where pw stands for plain weave, and f stands for float region. Figure 3 refers to a symbolic representation of interlaced (rigid region) and noninterlaced (loose region) areas. The float portion would be experiencing an extra force in the left that would not be balanced by its frictional force as shown in Figure 3. Figure 4 demonstrates a line diagram for the weave, where the zig-zag solid line denotes the point of interlacement of warp and weft, and the space between them constitutes the float area. The dotted line refers to the line connecting the midpoints of the twill line. This will move the center of the float towards the line of center. From Figure 4, the floating part gets another extra force towards the left due to the difference of accumulated friction force of the plain weave and float part; point A moves toward the line of center. After analyzing the single weft, while stretching, the center of the float will deflect toward the side (w.r.t line of center) where there is a significant portion of plain weave surrounding the float. Point B will move towards the line of the center from the left and point C will travel towards the line of the center from the right. All the weft in between point A, B, C will follow the course according to the argument as mentioned earlier and the center of float in each weft will try to come closer towards each other in the axial direction. The virtual line AB and BC would be formed where all the center of floats will lie on that line, and their tendency of motion will force the acute angle between the lines to overlap with the line of center. This is due to couple formation in the virtual lines AB and BC individually which restricts the motion of point B in the y-axis (weft direction). Point B transfers its component force to the lines AB and BC in upward and downward directions, respectively.
Textiles 2021, 1, FOR PEER REVIEW 4 In the unit cell of the design, there are 18 warp and 17 weft yarns. Considering each pick, there are two portions of different weave patterns, for example, the first weft is interlaced with nine warp yarns in a plain weave pattern, and the rest is with 1-up-4-down twill (float). The weft has two yarns that are elastic and stiffer yarn. The elastic yarn will be responsible for the extension in the floating part initially. The pull-out force of yarn from a plain-woven fabric is much more compared to that of 1 up 4 down twill weaves, as the interlacements in the case of plain weave are more. The stiffer weft yarn would require force to overcome the friction resistance. The stress at both sides of the fabric is equal and varying as the rate of displacement of the ends of the fabric is constant. Solving for the first pick, there are two segments-plain weave and float-the plain weave is toward the left side of the line of center. The force on both sides (connecting point) of the portion is different due to the non-symmetry of the structure. The force acting on the right side of the plain weave part would be more than on the left side. This would be balanced by the frictional force at each contact, which would accumulate in one direction due to differences in the end forces. The frictional force at each contact point in plain weave and float is given by Pan [50] as shown in Equations (1)-(3).

=
(1) = 2 sin (For uniform plain weave) (2) where P is the compression force acting on the point of intersection of warp and weft yarns, θ is the weave angle, and σ is the applied stress, is coefficient of friction, where pw stands for plain weave, and f stands for float region. Figure 3 refers to a symbolic representation of interlaced (rigid region) and non-interlaced (loose region) areas. The float portion would be experiencing an extra force in the left that would not be balanced by its frictional force as shown in Figure 3.  Figure 4, the floating part gets another extra force towards the left due to the difference of accumulated friction force of the plain weave and float part; point A moves toward the line of center. After analyzing the single weft, while stretching, the center of the float will deflect toward the side (w.r.t line of center) where there is a significant portion of plain weave surrounding the float. Point B will move towards the line of the center from the left and point C will travel towards the line of the center from the right. All the weft in between point A, B, C will follow the course according to the argument as mentioned earlier and the center of float in each weft will try to come closer towards each other in the axial direction. The virtual line AB and BC would be formed where all the center of floats will lie on that line, and their tendency of motion will force the acute angle between the lines to overlap with the line of center. This is due to couple formation in the virtual lines AB and BC individually which restricts the motion of point B in the y-axis (weft direction). Point B transfers its component force to the lines AB and BC in upward and downward directions, respectively.   The extension in the plain weave that would contribute to the overall extension of fabric in the axial direction would be due to the stiffer part (de-crimping) and the extension in the elastic yarn would be the same as of stiffer yarn due to interlacement present in the plain weave. The extension in that part would be according to the saw teeth model as shown in Figure 5, which is a representation of the unit cell of interlacement of warp and weft. where applied forces are not large, considering the bending energy. This analysis is useful for understanding the initial extension and de-crimping of the fabric under uniaxial load. All these notations refer to Figure 5, where the p is the spacing between the yarns, V is the inter-yarn pressure at crossover points (here V = P), h is the crimp height, it is otherwise called crimp amplitude where the amplitude is the maximum displacement of yarn from the central axis. It may be mentioned that crimp in the yarn of woven fabric makes a waveform, B is the bending rigidity of the yarn, l refers to the length of yarn between two consecutive crossover points, 1 and 2 represent warp and weft, respectively. Accordingly, to this model, the extension in the warp direction is given by The extension in the plain weave that would contribute to the overall extension of fabric in the axial direction would be due to the stiffer part (de-crimping) and the extension in the elastic yarn would be the same as of stiffer yarn due to interlacement present in the plain weave. The extension in that part would be according to the saw teeth model as shown in Figure 5, which is a representation of the unit cell of interlacement of warp and weft. where applied forces are not large, considering the bending energy. This analysis is useful for understanding the initial extension and de-crimping of the fabric under uniaxial load.  The extension in the plain weave that would contribute to the overall extension of fabric in the axial direction would be due to the stiffer part (de-crimping) and the extension in the elastic yarn would be the same as of stiffer yarn due to interlacement present in the plain weave. The extension in that part would be according to the saw teeth model as shown in Figure 5, which is a representation of the unit cell of interlacement of warp and weft. where applied forces are not large, considering the bending energy. This analysis is useful for understanding the initial extension and de-crimping of the fabric under uniaxial load. All these notations refer to Figure 5, where the p is the spacing between the yarns, V is the inter-yarn pressure at crossover points (here V = P), h is the crimp height, it is otherwise called crimp amplitude where the amplitude is the maximum displacement of yarn from the central axis. It may be mentioned that crimp in the yarn of woven fabric makes a waveform, B is the bending rigidity of the yarn, l refers to the length of yarn between two consecutive crossover points, 1 and 2 represent warp and weft, respectively. Accordingly, to this model, the extension in the warp direction is given by All these notations refer to Figure 5, where the p is the spacing between the yarns, V is the inter-yarn pressure at crossover points (here V = P), h is the crimp height, it is otherwise called crimp amplitude where the amplitude is the maximum displacement of yarn from the central axis. It may be mentioned that crimp in the yarn of woven fabric makes a waveform, B is the bending rigidity of the yarn, l refers to the length of yarn between two consecutive crossover points, 1 and 2 represent warp and weft, respectively. Accordingly, to this model, the extension in the warp direction is given by The final extension in the plain weave part is multiplied by the number of plain weave interlacements (n) The extension in the floating part would be because the elastic yarn is in taut condition and the stiffer yarn is in slack condition. The float is extended till the stiffer yarn is in a taut condition.
where ∆F f the difference in the force of friction of the two portions, L is the float length, E l is the elasticity of the weft yarn, δx 2 is the final extension in the float region. Now, the total extension would be the sum of extension in the plain weave and float part. The strain in the axial direction is total extension divided by the length of the fabric (in this case unit cell). Now, considering the forces that are acting on warp yarns, as there is no external stress applied in the transversal direction, the frictional force, and the transferred force due to the restricted motion of Point B would be in action. In the warp yarn also, there is a tightly woven and loosely woven part. Solving for the warp containing point A and C, the extension in plain weave is The float will also extend according to the same argument but with a different compression force. Similarly, we add the extension of both the sections and dividing by the total length of warp yarn. The strain in the transverse direction would be positive and finally, the Poisson's ratio would be negative until the tension (applied force) becomes too large to decrease the angle of the weave to zero. After this point, Poisson's ratio becomes positive because the yarn extension and strength come into the role. The above model is valid for the tension less than the critical point.

Material
Auxetic fabric is developed with foldable strips created in a zig-zag way in a vertical (warp) direction. It is based on the principle that when the fabric is stretched unfolding of the folds takes place, leading to an increase in transverse dimensions. Foldable geometry was realized into the fabric using a combination of loose and tight weaves to create differential shrinkage [14][15][16][17]. The unit cell of interlacement pattern of two single-layered fabrics based on parallel in-phase zigzag alternately folded stripes running along the warp direction as shown in Figure 6. The yarn used as weft was elastic (core spun spandex) with a fineness of 38 tex and non-elastic yarn (cotton) with a fineness of 20 tex. The yarn in the warp was non-elastic yarn (cotton) with a fineness of 2/30 tex. These yarns are used to exploit the differential shrinkage effect. The thread density in both directions is 16/cm. Mainly three types of yarns are used to exploit the auxetic nature of the developed fabric. The stretch properties The yarn used as weft was elastic (core spun spandex) with a fineness of 38 tex and non-elastic yarn (cotton) with a fineness of 20 tex. The yarn in the warp was non-elastic yarn (cotton) with a fineness of 2/30 tex. These yarns are used to exploit the differential shrinkage effect. The thread density in both directions is 16/cm. Mainly three types of yarns are used to exploit the auxetic nature of the developed fabric. The stretch properties of these yarns are presented in Table 1.

Test Method
Testing of all the specimen samples is performed on a video extensometer shown in Figure 7. The video extensometer is designed to accurately measure specimen strain during a materials test without the need to contact the specimen. The non-contacting measuring process leads to no mechanical influence on the specimen, ease of use, and reproducibility. The developed fabric has extensibility in one direction therefore test was carried out along this direction however the fabric can be tested along the warp direction but due to the limited extensibility in this direction, the fabric may break at smaller strain and cannot be tested over a wider range of strain. The tensile test was conducted on an Instron 5982 tensile machine. The capacity of the load cell used was 100 kN. The tensile speed and the gauge length were set as 50 mm/min and 100 mm respectively. Fabric strips of dimension (100 mm * 50 mm) were cut for each sample. Four points were marked by keeping the scale at the center of the fabric at a 25 mm distance to simplify the recording of the information of the fabric deformation during the tensile test. The distances of the two marks in the tensile and the traversal directions were scanned by a camera to record the extension for each sample until the sample broke and the software of the machining process the data to calculate the engineering strain of fabric in both tensile directions and the traversal directions. The fabric parameters are given in Table 2.  The tensile test was conducted on an Instron 5982 tensile machine. The capacity of the load cell used was 100 kN. The tensile speed and the gauge length were set as 50 mm/min and 100 mm respectively. Fabric strips of dimension (100 mm × 50 mm) were cut for each sample. Four points were marked by keeping the scale at the center of the fabric at a 25 mm distance to simplify the recording of the information of the fabric deformation during the tensile test. The distances of the two marks in the tensile and the traversal directions were scanned by a camera to record the extension for each sample until the sample broke and the software of the machining process the data to calculate the engineering strain of fabric in both tensile directions and the traversal directions. The fabric parameters are given in Table 2.

Computational Model
For the computational modeling of the foldable woven fabric, the unit cell of the model is firstly created on Texgen software, according to the weave design. Texgen can create a wide variety of 2-D and 3-D weaves [1]. Here, the yarns are considered as single unit, that is, solid volumes instead of considering the fiber in the yarns as individual units.
Geometry: The unit cell shown in Figure 8, is created by varying the float length along the warp and weft direction. Here, we use a combination of the plain weave along with the 4 up 1 down twill weave. Here, in each weft, 9 warp yarns are interlaced in plain weave and 10 warp yarns are interlaced in a 4/1 twill pattern, forming a foldable geometry. The twill weave forms the folded part with a long float length, which unfolds when extended in the axial direction. The unit cell created in Texgen, when joined, can produce full-scale fabric. This sort of portrayal is significantly more helpful for computational analysis of material properties.

FEM:
The prepared unit cell is imported to Ansys Workbench as a Stereo Lithography file to FEM Modeler to generate the initial geometry. Then these data are transferred to the transient structural feature. Then the type of material is assigned to the yarns, that is, the weft yarns were kept elastic and warp yarns are kept non-elastic.
Contact: Initially the contact is created automatically on Ansys mechanical, analyzing the face-to-face detection shown in Figure 9 but the contacts are modified according to weave design. The contact method is changed to frictional with the coefficient of friction equal to 0.1. The formulation is set to Augmented Lagrange as the penetration matter in this model. The initial contact condition is analyzed through the contact tool and accordingly pinball radius is adjusted manually.

FEM:
The prepared unit cell is imported to Ansys Workbench as a Stereo Lithography file to FEM Modeler to generate the initial geometry. Then these data are transferred to the transient structural feature. Then the type of material is assigned to the yarns, that is, the weft yarns were kept elastic and warp yarns are kept non-elastic.
Contact: Initially the contact is created automatically on Ansys mechanical, analyzing the face-to-face detection shown in Figure 9 but the contacts are modified according to weave design. The contact method is changed to frictional with the coefficient of friction equal to 0.1. The formulation is set to Augmented Lagrange as the penetration matter in this model. The initial contact condition is analyzed through the contact tool and accordingly pinball radius is adjusted manually. is, the weft yarns were kept elastic and warp yarns are kept non-elastic.
Contact: Initially the contact is created automatically on Ansys mechanical, analyzing the face-to-face detection shown in Figure 9 but the contacts are modified according to weave design. The contact method is changed to frictional with the coefficient of friction equal to 0.1. The formulation is set to Augmented Lagrange as the penetration matter in this model. The initial contact condition is analyzed through the contact tool and accordingly pinball radius is adjusted manually.

Mesh generation:
After the contact point generation, the mesh was generated as shown in Figure 10. The mesh size taken here was medium for easy processing and accuracy and element size was 0.125 mm by default.
Textiles 2021, 1, FOR PEER REVIEW 10 Mesh generation: After the contact point generation, the mesh was generated as shown in Figure 10. The mesh size taken here was medium for easy processing and accuracy and element size was 0.125 mm by default. Boundary conditions: The boundary conditions are imposed, as shown in Figure 11, on the unit cell, keeping in mind that it is a repeating unit and not a full-scale fabric. The following kinematic boundary condition will be applied as in Equations (9) and (10).
where v1, v2, v3, and v4 are displacement vectors of points on the edges of 1, 2, 3, and 4, respectively, and x1, x2, x3 and x4 are the absolute value of distances corresponding to displacement vectors v1, v2, v3, and v4. E stands for macroscopic strain tensor of the fabric [51][52][53][54]. Boundary conditions: The boundary conditions are imposed, as shown in Figure 11, on the unit cell, keeping in mind that it is a repeating unit and not a full-scale fabric. The following kinematic boundary condition will be applied as in Equations (9) and (10).
where v 1 , v 2 , v 3 , and v 4 are displacement vectors of points on the edges of 1, 2, 3, and 4, respectively, and x 1 , x 2 , x 3 and x 4 are the absolute value of distances corresponding to displacement vectors v 1 , v 2 , v 3 , and v 4 . E stands for macroscopic strain tensor of the fabric [51][52][53][54].

Simulation:
The left boundary was fixed along the x-direction and a directional displacement was applied along the right boundary, with the displacement rate of 2 mm/min. For the tensile test environment simulation, z and y directional displacements were constrained on both left and right boundaries. For the mechanical analysis of the model, elastic strain in the transverse direction in Figure 12.
Textiles 2021, 1, FOR PEER REVIEW 11 Simulation: The left boundary was fixed along the x-direction and a directional displacement was applied along the right boundary, with the displacement rate of 2 mm/min. For the tensile test environment simulation, z and y directional displacements were constrained on both left and right boundaries. For the mechanical analysis of the model, elastic strain in the transverse direction in Figure 12.

Results and Discussion
The prepared fabric samples were tested for the tensile deformation under a video extensometer as shown in Figure 7. The Poisson's ratios at different axial strain were plotted. This experimental deformation behavior was compared with the results obtained from analytical and computational model results. The graphical representation of the NPR effect is shown as Poisson's ratio vs. axial strain values.

Experimental Results
The poison's ratio was measured for the auxetic foldable design according to the procedure explained in Section 2.2. On analysis of Figure 13, it may be observed that fabric

Results and Discussion
The prepared fabric samples were tested for the tensile deformation under a video extensometer as shown in Figure 7. The Poisson's ratios at different axial strain were plotted. This experimental deformation behavior was compared with the results obtained from analytical and computational model results. The graphical representation of the NPR effect is shown as Poisson's ratio vs. axial strain values.

Experimental Results
The poison's ratio was measured for the auxetic foldable design according to the procedure explained in Section 2.2. On analysis of Figure 13, it may be observed that fabric with foldable strips created in a zig-zag way in the vertical (warp) direction produced NPR up to 8.7% of axial strain and maximum Poisson's ratio of (−0.41) produced at an axial strain of around 1%. Although the NPR effect was lost when the strain percentage exceeded 9%. The increase in Poisson's ratio with the increase of strain is due to the force of up-crimping in the traverse direction. This up-crimping is increasing as the applied load increases, resulting in the shrinking behavior of the sample in the transverse direction.

Analytical Model Results
The axial and transverse strain obtained by solving Equations [4][5][6][7][8] at tension 0.13296 N is given in Table 3. So, the maximum strain in the transverse direction is 0.000834 and the strain in the axial direction is 0.0014837. The Poisson's ratio calculated from the above figures is −0.5623. This value is calculated for a particular point, as the tension increases the sin(θ) decreases leading to a change in Poisson's ratio. Drawing a continuous graph from the analytical model is difficult due to the non-linear dependence of θ on applied force. The above value calculated is at a particular value of applied force in which the maximum Poisson's ratio was obtained experimentally.

Computational Model Results
The Poisson's ratio obtained through simulation data has a good relationship with the experimental data. Under unidirectional displacement, a series of quasi-uniform wave forms was observed in the transverse direction. Figure 14 shows that the value of Poisson's ratio is dependent on the extent of axial strain, that is, as the axial strain increases the NPR effect decreases. The trend observed here is close to the experiment values, the variation is due to the difference between contact simulations and actual contact conditions. From Figure 12, it was observed that the float region undergoes maximum elonga-

Analytical Model Results
The axial and transverse strain obtained by solving Equations [4][5][6][7][8] at tension 0.13296 N is given in Table 3. So, the maximum strain in the transverse direction is 0.000834 and the strain in the axial direction is 0.0014837. The Poisson's ratio calculated from the above figures is −0.5623. This value is calculated for a particular point, as the tension increases the sin(θ) decreases leading to a change in Poisson's ratio. Drawing a continuous graph from the analytical model is difficult due to the non-linear dependence of θ on applied force. The above value calculated is at a particular value of applied force in which the maximum Poisson's ratio was obtained experimentally.

Computational Model Results
The Poisson's ratio obtained through simulation data has a good relationship with the experimental data. Under unidirectional displacement, a series of quasi-uniform wave forms was observed in the transverse direction. Figure 14 shows that the value of Poisson's ratio is dependent on the extent of axial strain, that is, as the axial strain increases the NPR effect decreases. The trend observed here is close to the experiment values, the variation is due to the difference between contact simulations and actual contact conditions. From Figure 12, it was observed that the float region undergoes maximum elongation, due to the unrestricted motion of the elastic yarn. Moreover, the variation in the result is due to the weft yarn arrangement, that is, for actual fabric construction, alternate elastic and non-elastic yarns are used, whereas for the FEM simulation, consecutive elastic yarns are considered as the weft.

Comparison of the computational and experimental results
From Figure 14, it was observed that the analytical model results are in good agreement with the experimental results despite having the limitation to predict the trend. The error percentage of the results of this model is 37.14% for experimental results. The computational results also predict the Poisson's ratio with an error percentage of 22.26%. The reason for this error is the difference between the behavior yarn at interlacement points and this increases with the increase in the number of interlacements.

Conclusions
The analytical and computational model was built to predict the variation of Poisson's ratio of the auxetic fabric with the axial strain. The model was based on the unit cell of the foldable stripes created in a zig-zag way. In the analytical model, Poisson's ratio is based on the crimp percentage, bending modulus, yarn spacing, and coefficient of friction. In this yarn, properties and fabric, parameters were also considered. Also, successful computational modeling of the same was done. Both the analytical and computational models gave close predictions to the experimental results. It was observed that fabric with foldable strips created in a zig-zag way in vertical (warp) direction produced a negative Poisson's ratio (NPR) of 8.7%. Poisson's ratio of (−0.41) was produced at an axial strain of around 1%. The NPR effect was lost when the strain percentage exceeded 9%. The error percentage in the analytical model was 37.14% for experimental results. The computational results also predict the Poisson's ratio with an error percentage of 22.26%. Such predictions are useful for estimating the performance of auxetic woven structures in composite reinforcement. The auxetic structure exhibits remarkable stress-strain behavior in

Comparison of the Computational and Experimental Results
From Figure 14, it was observed that the analytical model results are in good agreement with the experimental results despite having the limitation to predict the trend. The error percentage of the results of this model is 37.14% for experimental results. The computational results also predict the Poisson's ratio with an error percentage of 22.26%. The reason for this error is the difference between the behavior yarn at interlacement points and this increases with the increase in the number of interlacements.

Conclusions
The analytical and computational model was built to predict the variation of Poisson's ratio of the auxetic fabric with the axial strain. The model was based on the unit cell of the foldable stripes created in a zig-zag way. In the analytical model, Poisson's ratio is based on the crimp percentage, bending modulus, yarn spacing, and coefficient of friction. In this yarn, properties and fabric, parameters were also considered. Also, successful computational modeling of the same was done. Both the analytical and computational models gave close predictions to the experimental results. It was observed that fabric with foldable strips created in a zig-zag way in vertical (warp) direction produced a negative Poisson's ratio (NPR) of 8.7%. Poisson's ratio of (−0.41) was produced at an axial strain of around 1%. The NPR effect was lost when the strain percentage exceeded 9%.
Textiles 2022, 2 13 The error percentage in the analytical model was 37.14% for experimental results. The computational results also predict the Poisson's ratio with an error percentage of 22.26%. Such predictions are useful for estimating the performance of auxetic woven structures in composite reinforcement. The auxetic structure exhibits remarkable stress-strain behavior in longitudinal as well as transverse directions. This performance is useful for energy absorption in a composite reinforcement.