An Analytical Model for the Tension-Shear Coupling of Woven Fabrics with Di ﬀ erent Weave Patterns under Large Shear Deformation

Featured Application: The work provides an analytical model which can be useful in the optimization of process parameters of dry fabrics in material-forming industries. Abstract: It is essential to accurately describe the large shear behavior of woven fabrics in the composite preforming process. An analytical model is proposed to describe the shear behavior of fabrics with di ﬀ erent weave patterns, in which tension-shear coupling is considered. The coupling is involved in two parts, the friction between overlapped yarns and the in-plane transverse compression between two parallel yarns. By introducing the concept of inﬂection points of a yarn, the model is applicable for fabrics with di ﬀ erent weave patterns. The analytical model is validated by biaxial tension-shear experiments. A parametric study is conducted to investigate the e ﬀ ects of external load, yarn geometry, and weave structure on the large shear behavior of fabrics. The developed model can reveal the physical mechanism of tension-shear coupling of woven fabrics. Moreover, the model has a high computational e ﬃ ciency due to its explicit expressions, thus beneﬁting the material design process.


Introduction
Woven fabrics and their composites have been widely used in engineering due to the outstanding mechanical performance [1,2]. Recently, wet compression molding of woven composites has attracted more attention because of its high efficiency in the manufacturing process [3]. The shear deformation, the dominant mode of a woven fabric in the forming process, is closely related to the yarn reorientation which is critical for the mechanical properties of a composite. It is noted that a blank holder is commonly used to provide a tensile force in a preforming process so that wrinkle onset can be postponed [4]. Moreover, the shear stiffness of a woven fabric can be significantly affected by applying a tensile force Figure 1 shows the configuration of a plain-woven representative volume element (RVE). The RVE geometric parameters, w 0 , g 0 , and N y , are, respectively, the yarn width, the yarn gap, and the yarn amount in the RVE. L R , the side length of the RVE, is defined as follows:

Geometrical Model
(1) Appl. Sci. 2020, 10, x FOR PEER REVIEW 3 of 15 Figure 1 shows the configuration of a plain-woven representative volume element (RVE). The RVE geometric parameters, 0 , 0 , and , are, respectively, the yarn width, the yarn gap, and the yarn amount in the RVE. , the side length of the RVE, is defined as follows: In order to focus on the tension-shear coupling, this work only considers the case that the warp and fill yarns have identical geometry. It should be noted that Equation (1) is applicable for fabrics with different weave patterns. The plain weave pattern shown in Figure 1 is only an example for illustration. The number of inflection points of a yarn, is introduced to characterize different weave patterns. In a woven fabric, a warp yarn passes over and under fill yarns in sequence, and vice versa. The ratio of the under and over areas of warp and fill yarns determines the weave pattern. In this work, when a yarn turns from under to over or from over to under one time, it is counted as one inflection point. The inflection point at the boundary of an RVE is counted as a half point. For example, as shown in Figure 2a, the inflection point number of a warp yarn is 2. The weave pattern of RVEs with the same yarn amounts can be distinguished by different yarn inflection point numbers. For example, as shown in Figure 2a,b, both the twill and plain RVEs have 4 warp and 4 fill yarns, but their inflection point numbers of the warp yarn are 2 and 4, respectively. The number of inflection points in an RVE determines the normal pressure between overlapped yarns induced by the tensile force. Furthermore, normal pressure is closely related to the tension-shear coupling which will be explained in detail hereinafter. Hence, it is believed that the inflection point number is critical to evaluate the effects of different weave patterns on the tension-shear coupling of fabrics.

Incremental Energy Equilibrium Theory
In forming manufacturing industries, the shape and stress state of a formed fabric depend on the forming process parameters, such as molding shape, forming rate, etc. Numerical simulations are generally used to optimize such forming process parameters. The key to such a simulation is a reliable In order to focus on the tension-shear coupling, this work only considers the case that the warp and fill yarns have identical geometry. It should be noted that Equation (1) is applicable for fabrics with different weave patterns. The plain weave pattern shown in Figure 1 is only an example for illustration.
The number of inflection points of a yarn, η y is introduced to characterize different weave patterns. In a woven fabric, a warp yarn passes over and under fill yarns in sequence, and vice versa. The ratio of the under and over areas of warp and fill yarns determines the weave pattern. In this work, when a yarn turns from under to over or from over to under one time, it is counted as one inflection point. The inflection point at the boundary of an RVE is counted as a half point. For example, as shown in Figure 2a, the inflection point number of a warp yarn is 2. The weave pattern of RVEs with the same yarn amounts can be distinguished by different yarn inflection point numbers. For example, as shown in Figure 2a,b, both the twill and plain RVEs have 4 warp and 4 fill yarns, but their inflection point numbers of the warp yarn are 2 and 4, respectively. The number of inflection points in an RVE determines the normal pressure between overlapped yarns induced by the tensile force. Furthermore, normal pressure is closely related to the tension-shear coupling which will be explained in detail hereinafter. Hence, it is believed that the inflection point number is critical to evaluate the effects of different weave patterns on the tension-shear coupling of fabrics.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 3 of 15 Figure 1 shows the configuration of a plain-woven representative volume element (RVE). The RVE geometric parameters, 0 , 0 , and , are, respectively, the yarn width, the yarn gap, and the yarn amount in the RVE. , the side length of the RVE, is defined as follows: In order to focus on the tension-shear coupling, this work only considers the case that the warp and fill yarns have identical geometry. It should be noted that Equation (1) is applicable for fabrics with different weave patterns. The plain weave pattern shown in Figure 1 is only an example for illustration. The number of inflection points of a yarn, is introduced to characterize different weave patterns. In a woven fabric, a warp yarn passes over and under fill yarns in sequence, and vice versa. The ratio of the under and over areas of warp and fill yarns determines the weave pattern. In this work, when a yarn turns from under to over or from over to under one time, it is counted as one inflection point. The inflection point at the boundary of an RVE is counted as a half point. For example, as shown in Figure 2a, the inflection point number of a warp yarn is 2. The weave pattern of RVEs with the same yarn amounts can be distinguished by different yarn inflection point numbers. For example, as shown in Figure 2a,b, both the twill and plain RVEs have 4 warp and 4 fill yarns, but their inflection point numbers of the warp yarn are 2 and 4, respectively. The number of inflection points in an RVE determines the normal pressure between overlapped yarns induced by the tensile force. Furthermore, normal pressure is closely related to the tension-shear coupling which will be explained in detail hereinafter. Hence, it is believed that the inflection point number is critical to evaluate the effects of different weave patterns on the tension-shear coupling of fabrics.

Incremental Energy Equilibrium Theory
In forming manufacturing industries, the shape and stress state of a formed fabric depend on the forming process parameters, such as molding shape, forming rate, etc. Numerical simulations are generally used to optimize such forming process parameters. The key to such a simulation is a reliable

Incremental Energy Equilibrium Theory
In forming manufacturing industries, the shape and stress state of a formed fabric depend on the forming process parameters, such as molding shape, forming rate, etc. Numerical simulations are generally used to optimize such forming process parameters. The key to such a simulation is a reliable constitutive relation that describes the mechanical behavior of an RVE of a fabric. In practice, blank-holders are often used to introduce tension forces along yarns so that the wrinkling phenomenon may be postponed. Besides, a fabric is majorly shear deformed in the forming process. Thus, it is critical to establish a constitutive relation for large shear deformation of a fabric RVE with the consideration of tension-shear coupling. Once the constitutive relationship is obtained, the numerical simulation can be realized subsequently. Thus, this work mainly investigates a fabric RVE subjected to a combined tension and shear load, as shown in Figure 3.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 4 of 15 constitutive relation that describes the mechanical behavior of an RVE of a fabric. In practice, blankholders are often used to introduce tension forces along yarns so that the wrinkling phenomenon may be postponed. Besides, a fabric is majorly shear deformed in the forming process. Thus, it is critical to establish a constitutive relation for large shear deformation of a fabric RVE with the consideration of tension-shear coupling. Once the constitutive relationship is obtained, the numerical simulation can be realized subsequently. Thus, this work mainly investigates a fabric RVE subjected to a combined tension and shear load, as shown in Figure 3. The RVE is applicable for different weave patterns. Thus, no specific weave pattern is illustrated in Figure 3. The RVE is firstly subjected to tension loads along yarns, . The tension forces remain constant throughout the load case. Then, the RVE undergoes shear deformation under the force at the node P of the RVE. The node P moves to P' due to the shear deformation. The shear force ℎ is the projection of to the boundary of the RVE. The relationship between ℎ and is given as follows: As mentioned in the above paragraph, the tension forces along yarns are applied before the shear deformation, meaning that a yarn has been stretched before the RVE is sheared. Thus, no tension input energy is involved in the shear deformation. The input energy generated by the node force will be dissipated by friction, transverse compression, thickness variation, yarn elongation, yarn sliding, and waviness angle variation. Among these dissipation energies, the former two are dominant. The yarn thickness may increase due to in-plane transverse compression, but the variation is small since the increase of yarn thickness will be restricted by overlapped yarns. The yarn elongation is negligible, since the longitudinal modulus of a yarn is much larger than the shear and transverse compressive properties of a fabric. The yarn sliding is also negligible, since a couple of balanced tension forces are applied to the end of a yarn. The waviness angle variation mainly results from yarn sliding and thickness variation, thereby also being insignificant. Thus, only the friction and transverse compression dissipation energies are considered. According to the energy equilibrium theory, the input energy should equal the sum of the dissipation energies. Hence, the following equation is obtained: where d is the input energy increment generated by the node force, . d and d are incremental energy dissipated by the respective friction and transverse compression of yarns.
Our goal is to express all items in Equation (3) as functions of shear force, shear angle, geometric parameters, and mechanical properties of yarns so that a shear stress-angle curve can be obtained eventually. The RVE is applicable for different weave patterns. Thus, no specific weave pattern is illustrated in Figure 3. The RVE is firstly subjected to tension loads along yarns, F t . The tension forces remain constant throughout the load case. Then, the RVE undergoes shear deformation under the force F node at the node P of the RVE. The node P moves to P' due to the shear deformation. The shear force F sh is the projection of F node to the boundary of the RVE. The relationship between F sh and F node is given as follows: As mentioned in the above paragraph, the tension forces along yarns are applied before the shear deformation, meaning that a yarn has been stretched before the RVE is sheared. Thus, no tension input energy is involved in the shear deformation. The input energy generated by the node force F node will be dissipated by friction, transverse compression, thickness variation, yarn elongation, yarn sliding, and waviness angle variation. Among these dissipation energies, the former two are dominant. The yarn thickness may increase due to in-plane transverse compression, but the variation is small since the increase of yarn thickness will be restricted by overlapped yarns. The yarn elongation is negligible, since the longitudinal modulus of a yarn is much larger than the shear and transverse compressive properties of a fabric. The yarn sliding is also negligible, since a couple of balanced tension forces are applied to the end of a yarn. The waviness angle variation mainly results from yarn sliding and thickness variation, thereby also being insignificant. Thus, only the friction and transverse compression dissipation energies are considered. According to the energy equilibrium theory, the input energy should equal the sum of the dissipation energies. Hence, the following equation is obtained: Appl. Sci. 2020, 10, 1551 where dW node is the input energy increment generated by the node force, F node . dW f ri and dW com are incremental energy dissipated by the respective friction and transverse compression of yarns. Our goal is to express all items in Equation (3) as functions of shear force, shear angle, geometric parameters, and mechanical properties of yarns so that a shear stress-angle curve can be obtained eventually.

Input Energy of Shear Force
As shown in Figure 3, when a fabric RVE undergoes shear deformation, the node P moves to P with displacement of S node . The yarn angle α decreases with the increase of shear angle ϕ, where a relation exists: α = π/2 − ϕ. Then, the incremental input energy is given by In order to characterize the shear behavior of the RVE, it is necessary to express F node and dS node as functions of the shear stress and shear angle. Firstly, with Equation (2), the following equation is given as: Furthermore, for the purpose of normalization, the shear stress τ sh instead of the shear force F sh is employed. Then, we have where L R is the side length of the RVE, and t is the yarn thickness. Secondly, dS node is the differentiation of the node displacement S node , whose relationship is given by From Figure 3, S node is given by Then, dS node is given by Finally, substituting Equations (5), (6), and (9) into Equation (4), the input energy is expressed in terms of the shear stress and the geometric parameters as follows:

Dissipation Energy of Friction
Figure 4a is a crossover point of two yarns. As shown in Figure 4b, a friction force is generated when a yarn rotates. An equivalent circle is used to approximate the contact area of overlapped yarns. It is reasonable to assume that the friction force is evenly distributed in the contact area. Imagine an infinitesimal sector whose side lengths are dr and rdφ, as shown in Figure 4c. The sector moves from A to A' with the yarn rotation. Then, the moment of the friction force to the point O is given by Appl. Sci. 2020, 10, 1551 6 of 16 0 in the following equation: Ideally, 0 should be calculated from the yarn bending stiffness and the bending angle. However, the measured bending stiffness has large scatter, and its value is sensitive to the load status of a yarn. Thus, it is better to curve-fit 0 by comparing prediction results to experimental data.   Figure 5 illustrates the contact force induced by the tension force along a yarn for a plain-woven fabric. is the waviness angle of a yarn.
is the length of the inflected part of a yarn. Then, is: For fabrics with different weave patterns, the illustration of the contact force is shown in Figure  6. Then, the contact force at one cross point is: Furthermore, the friction dissipation energy in an RVE is obtained as where R is the equivalent radius. The contact area varies with the shear angle until parallel yarns contact each other in in-plane transverse direction. Thus, R is expressed as where ϕ L is the shear-locking angle, meaning that parallel yarns start to contact each other in in-plane transverse direction. The shear-locking angle can be expressed as It should be noted that the shear angle keeps increasing after ϕ L due to in-plane transverse compression of yarns. In other words, the shear angle is actually not locked. However, following the nomenclature of textile forming, the shear-locking angle is still employed.
For simplicity, Coulomb's friction theory is employed. Before a yarn rotates, the static friction force can be approximately seen as linearly dependent on the shear angle. When a yarn rotates, the friction force F f ri is calculated from the friction factor µ and the contact force between overlapped yarns F c . The friction force is expressed as: where ϕ 0 is the critical shear angle when a yarn starts to rotate. Obviously, the friction force depends on the contact force between overlapped yarns, F c . F c consists of two parts: the one induced by the tension force along a yarn, F c t , and the initial contact force induced by the yarn bending, F c 0 , as shown in the following equation: Ideally, F c 0 should be calculated from the yarn bending stiffness and the bending angle. However, the measured bending stiffness has large scatter, and its value is sensitive to the load status of a yarn. Thus, it is better to curve-fit F c 0 by comparing prediction results to experimental data.
Appl. Sci. 2020, 10, 1551 7 of 16 Figure 5 illustrates the contact force induced by the tension force along a yarn for a plain-woven fabric. θ is the waviness angle of a yarn. L in is the length of the inflected part of a yarn. Then, F c t is: (a) Cross point in an RVE (b) Illustration of friction force (c) Calculation of friction energy  Figure 5 illustrates the contact force induced by the tension force along a yarn for a plain-woven fabric. is the waviness angle of a yarn.
is the length of the inflected part of a yarn. Then, is: For fabrics with different weave patterns, the illustration of the contact force is shown in Figure  6. Then, the contact force at one cross point is: For fabrics with different weave patterns, the illustration of the contact force is shown in Figure 6.   Figure 5 illustrates the contact force induced by the tension force along a yarn for a plain-woven fabric. is the waviness angle of a yarn.
is the length of the inflected part of a yarn. Then, is: For fabrics with different weave patterns, the illustration of the contact force is shown in Figure  6. Then, the contact force at one cross point is: Then, the contact force at one cross point is: The friction force can be calculated by substituting Equation (18) into Equation (16) and then into Equation (15).
Finally, the dissipation energy, dW f ri , is obtained by substituting Equations (15) and (13) into Equation (12). The expression of dW f ri is obtained:

Dissipation Energy of in-plane Transverse Compression
As shown in Figure 7a, parallel yarns start to contact each other at the shear-locking angle. As the shear angle increases, the yarns are transversely compressed under the in-plane transverse compression force F com . The points C and D move to C' and D', and the total width reduction of a yarn is ∆w, as illustrated in Figure 7b. Since no yarn sliding is considered in this work, the side length of the parallelogram ABCD remains constant, which is (w 0 + g 0 ).

Dissipation Energy of in-plane Transverse Compression
As shown in Figure 7a, parallel yarns start to contact each other at the shear-locking angle. As the shear angle increases, the yarns are transversely compressed under the in-plane transverse compression force . The points C and D move to C' and D', and the total width reduction of a yarn is ∆ , as illustrated in Figure 7b. Since no yarn sliding is considered in this work, the side length of the parallelogram ABCD remains constant, which is ( 0 + 0 ). In an infinitesimal load increment, can be seen as constant. Since the fill and warp yarns are consistent in geometry and mechanical properties, the dissipation energy of the in-plane transverse compression of a fabric RVE is where is the differentiation of ∆ . ∆ and are given by where is the in-plane transverse compressive strain of a yarn. From the view point of simplicity, it is reasonable to assume that the yarn compression is linear. Then, the compression force is calculated with the following equation: where is the transverse compressive modulus of a yarn. Please remember that our goal is to obtain the constitutive relation between the shear stress and shear angle of a fabric RVE. Thus, must be written as a function of the shear angle. is the shear-locking angle given by Equation (14). Obviously, when ≤ , the width is not compressed, and both and are zero. When > , is the following equation is given as: Then, is obtained as Furthermore, is reached: In an infinitesimal load increment, F com can be seen as constant. Since the fill and warp yarns are consistent in geometry and mechanical properties, the dissipation energy of the in-plane transverse compression of a fabric RVE is where dw is the differentiation of ∆w. ∆w and dw are given by where ε com is the in-plane transverse compressive strain of a yarn. From the view point of simplicity, it is reasonable to assume that the yarn compression is linear. Then, the compression force F com is calculated with the following equation: where E com is the transverse compressive modulus of a yarn. Please remember that our goal is to obtain the constitutive relation between the shear stress and shear angle of a fabric RVE. Thus, ε com must be written as a function of the shear angle. ϕ L is the shear-locking angle given by Equation (14). Obviously, when ϕ ≤ ϕ L , the width is not compressed, and both ε com and dε com are zero. When ϕ > ϕ L , is the following equation is given as: Then, ε com is obtained as Furthermore, dε com is reached: Obviously, when ϕ ≤ ϕ L , the dissipation energy is zero. When ϕ > ϕ L , dW com can be achieved by substituting Equations (25) and (26) into Equations (22) and (23) and then into Equation (20). Finally, the expression of dW com is achieved: Appl. Sci. 2020, 10, 1551 9 of 16

Model Summary
Both the input energy and dissipation energies of friction and transverse compression are expressed as functions of the shear stress and shear angle. Then, the relation between shear stress and shear angle can be obtained by solving the energy equilibrium, Equation (3). It is noted that the friction dissipation energy is sensitive to the number of inflection points in an RVE, indicating that the present model can consider the effect of different weave patterns on the shear deformation of fabrics. In addition, the tension-shear coupling mechanism appears from two aspects in the present model. Firstly, the effect of tension force on the shear property of a fabric is directly reflected by the shear friction dissipation energy. An increase in the tension force leads to the increase of contact force between overlapped yarns and then results in the increase of friction dissipation energy. Obviously, a larger dissipation energy means larger resistance in deformation. Secondly, the tension force will make a yarn shrink. Generally, the transverse compressive modulus will rise due to the shrinking effect. Furthermore, the increase of transverse compressive modulus also generates more dissipation energy. With the help of the analytical model developed in this work, the effect of tension on the shear behavior of fabrics can be analyzed explicitly.

Model Validation
In this part, the present model is validated by experimental results of a dry woven fabric under tension and shear loads. The experimental data for a carbon fiber plain-woven fabric subjected to biaxial tension and shear force were presented in Nosrat-Nezami et al. [22]. As shown in Figure 8, the experiment is constructed by adding a tensile device to a picture frame configuration, and a needle bar is set up between clamps to guarantee the free rotation of yarns. The geometric information of the test specimen is shown in Table 1. The friction factor between warp and fill yarns of a carbon fabric is 0.6 according to Najjar et al. [23] and Allaoui et al. [24]. Please note that the width and gap of warp and fill yarns are slightly different. The shear-locking angle is determined by geometry of the yarn family which is first in-plane transversely contacted. Thus, the width and gap of the fill yarn are used in the present model.
Obviously, when ≤ , the dissipation energy is zero. When > , d can be achieved by substituting Equations (25) and (26) into Equations (22) and (23) and then into Equation (20). Finally, the expression of d is achieved:

Model Summary
Both the input energy and dissipation energies of friction and transverse compression are expressed as functions of the shear stress and shear angle. Then, the relation between shear stress and shear angle can be obtained by solving the energy equilibrium, Equation (3). It is noted that the friction dissipation energy is sensitive to the number of inflection points in an RVE, indicating that the present model can consider the effect of different weave patterns on the shear deformation of fabrics. In addition, the tension-shear coupling mechanism appears from two aspects in the present model. Firstly, the effect of tension force on the shear property of a fabric is directly reflected by the shear friction dissipation energy. An increase in the tension force leads to the increase of contact force between overlapped yarns and then results in the increase of friction dissipation energy. Obviously, a larger dissipation energy means larger resistance in deformation. Secondly, the tension force will make a yarn shrink. Generally, the transverse compressive modulus will rise due to the shrinking effect. Furthermore, the increase of transverse compressive modulus also generates more dissipation energy. With the help of the analytical model developed in this work, the effect of tension on the shear behavior of fabrics can be analyzed explicitly.

Model Validation
In this part, the present model is validated by experimental results of a dry woven fabric under tension and shear loads. The experimental data for a carbon fiber plain-woven fabric subjected to biaxial tension and shear force were presented in Nosrat-Nezami et al. [22]. As shown in Figure 8, the experiment is constructed by adding a tensile device to a picture frame configuration, and a needle bar is set up between clamps to guarantee the free rotation of yarns. The geometric information of the test specimen is shown in Table 1. The friction factor between warp and fill yarns of a carbon fabric is 0.6 according to Najjar et al. [23] and Allaoui et al. [24]. Please note that the width and gap of warp and fill yarns are slightly different. The shear-locking angle is determined by geometry of the yarn family which is first in-plane transversely contacted. Thus, the width and gap of the fill yarn are used in the present model.  [22]. Table 1. Geometric information of the T700-12K carbon fiber plain-woven fabric specimen [22].  When a tensile force is applied to a yarn, due to the yarn shrinking and the out-of-plane compression between overlapped yarns, the fabric thickness decreases and the transverse compressive modulus of the yarn increases. However, the variations of thickness and transverse compressive modulus are difficult to be measured experimentally. According to the thickness-pressure relation shown in Figure 9, Ivanov and Lomov [25] suggested that the thickness variation and the transverse compressive of a yarn can be approximated as bilinear functions of the compressive force. Please note that, although a fabric is subjected to a tension force along the yarn direction, the yarn is actually compressed in the transverse direction due to the yarn shrinking and the out-of-plane compression between overlapped yarns. Thus, the transverse stress state of a yarn in this work is similar to the transvers compressive test shown in the work of Ivanov and Lomov [25]. Besides, the materials involved in Figures 8 and 9 are both carbon fabrics with a 12 K yarn bundle. Thus, the quasi-bilinear trend shown in Figure 9 is reasonably applicable for the validation in this work.

Specimen Size
100×100 mm 0.45 mm 4.11mm 0.78 mm 4.94 mm 0 mm When a tensile force is applied to a yarn, due to the yarn shrinking and the out-of-plane compression between overlapped yarns, the fabric thickness decreases and the transverse compressive modulus of the yarn increases. However, the variations of thickness and transverse compressive modulus are difficult to be measured experimentally. According to the thicknesspressure relation shown in Figure 9, Ivanov and Lomov [25] suggested that the thickness variation and the transverse compressive of a yarn can be approximated as bilinear functions of the compressive force. Please note that, although a fabric is subjected to a tension force along the yarn direction, the yarn is actually compressed in the transverse direction due to the yarn shrinking and the out-of-plane compression between overlapped yarns. Thus, the transverse stress state of a yarn in this work is similar to the transvers compressive test shown in the work of Ivanov and Lomov [25]. Besides, the materials involved in Figures 8 and 9 are both carbon fabrics with a 12 K yarn bundle. Thus, the quasi-bilinear trend shown in Figure 9 is reasonably applicable for the validation in this work. Figure 9. Thickness-compressive pressure of a dry carbon yarn [23].
For the cases in Nosrat-Nezami et al. [22], the maximum tensile load is 2.4 N/mm. If the contact force induced by yarn bending is ignored, the tensile load of 2.4 N/mm corresponds to a compressive pressure of about 0.05 MPa, according to Equations (16) and (18), which is located around the first linear segment in Figure 9. Thus, it is rational to assume that the thickness and transverse compressive moduli linearly vary with the tensile force. In addition, the bending stiffness also varies due to the yarn shrinking. Furthermore, the bending stiffness directly affects the initial contact force 0 induced by yarn waviness. Thus, it is also reasonable to set 0 as a linear function of the tension force. Then, the thickness , the compressive modulus , and the initial contact force 0 can be obtained by a linear curve-fitting process from the cases of no tension and 2.4N/mm tension. Besides, the shear angle when a yarn starts to rotate, 0 , can be approximately seen as a constant throughout the load process, since it is not sensitive to the contact force between overlapped yarns. Thus, one curve is enough to fit the 0 . Overall, the cases of no tension and 2.4 N/mm tension are employed to curve-fit model parameters, and the other cases are used for validation. All the curve-fitted parameters are shown in Table 2.  Figure 9. Thickness-compressive pressure of a dry carbon yarn [23].
For the cases in Nosrat-Nezami et al. [22], the maximum tensile load is 2.4 N/mm. If the contact force induced by yarn bending is ignored, the tensile load of 2.4 N/mm corresponds to a compressive pressure of about 0.05 MPa, according to Equations (16) and (18), which is located around the first linear segment in Figure 9. Thus, it is rational to assume that the thickness and transverse compressive moduli linearly vary with the tensile force. In addition, the bending stiffness also varies due to the yarn shrinking. Furthermore, the bending stiffness directly affects the initial contact force F c 0 induced by yarn waviness. Thus, it is also reasonable to set F c 0 as a linear function of the tension force. Then, the thickness t, the compressive modulus E com , and the initial contact force F c 0 can be obtained by a linear curve-fitting process from the cases of no tension and 2.4N/mm tension. Besides, the shear angle when a yarn starts to rotate, ϕ 0 , can be approximately seen as a constant throughout the load process, since it is not sensitive to the contact force between overlapped yarns. Thus, one curve is enough to fit the ϕ 0 . Overall, the cases of no tension and 2.4 N/mm tension are employed to curve-fit model parameters, and the other cases are used for validation. All the curve-fitted parameters are shown in Table 2. Table 2. Curve-fitting parameters for a carbon fiber plain-woven fabric.

Loading Cases No Tension 0.8 N/mm 1.6 N/mm 2.4 N/mm
0.25 0.83 1.41 2 Figure 10 compares the results of the present model and experimental data. The cases of no tension and 2.4 N/mm tension are used to curve-fit the parameters of the present model, t, ϕ 0 , F c 0 , and E com . The cases of 0.8 N/mm and 1.6 N/mm are used for validation. Equation (28) provides a method to quantitatively characterize the average error between experimental data and the predicted results.
where ϕ i is the evenly distributed shear angle point with an angle-interval of around 5 • . τ pr sh (ϕ i ) and τ ex sh (ϕ i ) are the corresponding predicted and measured shear stresses for each ϕ i . M is the number of collected data points of a curve. The average errors for the four load cases are shown in Figure 10. The average errors of the four cases are around 10%, which is good enough considering the scatter of measured data.

Parametric Study
It is essential to understand the effect of geometry and mechanical properties of yarns on the large shear behavior. In this work, a kind of carbon fiber woven fabric is taken as an example. It should be noted in Figure 10 that the unit of tension force along a yarn is N/mm, meaning that the force is homogenized by the side length of an RVE. Thus, the force along a yarn is dependent on the thickness of a fabric. Conversely, the parametric study requires that the tensile load is independent of thickness variation. Thus, a tensile stress with unit of MPa is used. As a reference, some basic parameters are reasonably set in Table 3. The sensitivity study is performed by varying one of these parameters while keeping the others constant.  Figure 11a shows the shear deformation curves of fabrics with different weave patterns. The weave pattern is characterized by ⁄ , the ratio between the inflection points of a yarn and the number of yarns in an RVE. The ratios are 1, 1/2, and 1/4 for plain woven, 2*2 twill, and 8 harness satin, respectively. As shown in Figure 11a, a higher ratio corresponds to a higher shear stiffness. The result is reasonable because a higher ratio means more interaction between fill and warp yarns. Figure 11b is the parametric study of the fabric geometry. It is found that the shear stiffness increases with the increase of yarn thickness but decreases with the increase of yarn gap. This is because the waviness angle becomes larger when a yarn becomes thicker. Thus, the contact force between overlapped yarns increases, though the tension force along the yarn remains unchanged. However, the increase of yarn gap reduces the waviness angle, thereby weakening the effect of the tension force on shear deformation. Besides, the shear-locking phenomenon is postponed with the increase of gap size, 25° for 0 = 0.5 , 34° for 0 = 1 , and 44° for 0 = 2 . The transverse compression of parallel yarns is one main source of the shear resistance. Thus, the rising shear-locking angle significantly reduces the shear rigidity of a fabric. For the no tension and 2.4 N/mm tension cases, the theoretical curves agree well with experimental results, proving that the present model can effectively describe the shear mechanism and the tension-shear coupling effect of fabrics under biaxial tension. For the cases of 0.8 N/mm and 1.6 N/mm, the deviations between predicted and experimental results are very small, verifying the prediction capability of the present model.
Comparing the curves of the four cases, it is found that the shear rigidity of a fabric significantly increases with the applied tension force. Thus, it is evident that the tension-shear coupling must be taken into account in the material forming process of fabrics. As shown in Figure 10, a stress-strain curve can be divided into two parts: the initial linear part and the nonlinear part. The initial part is induced by the static friction between overlapped yarns before a yarn rotates. The maximum value of the static friction force determines the initial shear stress value of the nonlinear part. A higher tension force leads to a larger maximum static friction force, thereby resulting in stiffer shear property of a fabric. Thus, the friction between overlapped yarns is one source of the tension-shear coupling. For the nonlinear part, the stiffness rapidly increases when the tension force varies from zero to 2.4 N/mm. This happens because the fiber density of a yarn becomes larger due to yarn shrinking induced by a tension force. A higher fiber density means a higher transverse compressive Young's modulus. Furthermore, higher transverse stiffness of yarns brings a larger shear resistance for the shear deformation of fabrics. Obviously, the transverse compression of yarns is one main source of the tension-shear coupling.
For general fabrics, there should be a flat part between the initial linear part and the nonlinear rising part. The flat part is induced by friction during shear deformation before transverse compression of parallel yarns occurs. However, the fill gap is zero for the case in Figure 10, meaning that parallel yarns contact each other when a yarn starts to rotate. Thus, no obvious flat part is observed in Figure 10.

Parametric Study
It is essential to understand the effect of geometry and mechanical properties of yarns on the large shear behavior. In this work, a kind of carbon fiber woven fabric is taken as an example. It should be noted in Figure 10 that the unit of tension force along a yarn is N/mm, meaning that the force is homogenized by the side length of an RVE. Thus, the force along a yarn is dependent on the thickness of a fabric. Conversely, the parametric study requires that the tensile load is independent of thickness variation. Thus, a tensile stress σ y with unit of MPa is used. As a reference, some basic parameters are reasonably set in Table 3. The sensitivity study is performed by varying one of these parameters while keeping the others constant. Table 3. Parameter values of a reference-woven fabric.

Parameters Value Parameters Value
Weave pattern plain Initial shear angle ϕ 0 2 • Yarn width w f & w w 5 mm Tension stress along yarns σ y 5 MPa Yarn gap g f & g w 0.5 mm Friction factor µ 0.6 Yarn thickness t f & t w 0.5 mm Compressive modulus of yarn E com 1 MPa Figure 11a shows the shear deformation curves of fabrics with different weave patterns. The weave pattern is characterized by η y /N y , the ratio between the inflection points of a yarn and the number of yarns in an RVE. The ratios are 1, 1/2, and 1/4 for plain woven, 2*2 twill, and 8 harness satin, respectively. As shown in Figure 11a, a higher ratio corresponds to a higher shear stiffness. The result is reasonable because a higher ratio means more interaction between fill and warp yarns. Figure 11b is the parametric study of the fabric geometry. It is found that the shear stiffness increases with the increase of yarn thickness but decreases with the increase of yarn gap. This is because the waviness angle becomes larger when a yarn becomes thicker. Thus, the contact force between overlapped yarns increases, though the tension force along the yarn remains unchanged. However, the increase of yarn gap reduces the waviness angle, thereby weakening the effect of the tension force on shear deformation. Besides, the shear-locking phenomenon is postponed with the increase of gap size, 25 • for g 0 = 0.5, 34 • for g 0 = 1, and 44 • for g 0 = 2. The transverse compression of parallel yarns is one main source of the shear resistance. Thus, the rising shear-locking angle significantly reduces the shear rigidity of a fabric. Figure 11c illustrates the effect of the yarn property and the tensile load on the shear behavior of a fabric. It should be noted that in real cases, the transverse compressive modulus E com is sensitive to the tensile load F t . However, in the sensitivity study, the two parameters are assumed to be independent so as to investigate the effect of each parameter on the shear behavior. It can be seen that the increase of tensile force along yarns results in a larger shear rigidity of a fabric. Moreover, each curve in Figure 11 can be divided into three parts, the initial linear part, the flat part, and the nonlinear rising part. Curves in Figure 11a-c indicate that the weave pattern, the yarn geometry, and the tensile load along yarns can only influence the maximum stress of the first part and the shear stress of the flat part but have no effect on the third part. However, as shown in Figure 11c, the growth of E com significantly increases the slope of the third part of a curve, thereby enhancing the shear rigidity of a fabric. This happens because the effects of the weave pattern, the geometry, and the tensile load on shear behavior are characterized by the friction force which controls the linear and the flat part. However, the third nonlinear rising part is mainly controlled by the transverse compression of yarns. load along yarns can only influence the maximum stress of the first part and the shear stress of the flat part but have no effect on the third part. However, as shown in Figure 11c, the growth of significantly increases the slope of the third part of a curve, thereby enhancing the shear rigidity of a fabric. This happens because the effects of the weave pattern, the geometry, and the tensile load on shear behavior are characterized by the friction force which controls the linear and the flat part. However, the third nonlinear rising part is mainly controlled by the transverse compression of yarns.

Conclusions
An analytical model is proposed to describe the tension-shear coupling effect of fabrics with different weave patterns under large shear deformation. Different weave patterns are characterized by the number of inflection points of yarns. In this work, the tension-shear coupling mechanism can be revealed from two aspects, the friction dissipation energy and the in-plane transverse compressive dissipation energy. On the one hand, the subjected tensile force along yarns increases the contact force between overlapped yarns, leading to the growth of the friction dissipation energy. Then, larger friction dissipation energy brings more resistance with respect to shear deformation, thereby increasing the shear rigidity of a fabric. On the other hand, the tensile load along a yarn makes the yarn shrink. The shrinking improves the transverse compressive modulus of the yarn. Besides, the transverse compressive modulus of a yarn is also improved due to the compressive force between overlapped yarns induced by the tension force. Then, the dissipation energy of the transverse compression of yarns becomes larger due to the increase of the mechanical properties of yarns, thus enhancing the shear rigidity of a fabric. In addition to the description of the tension-shear coupling mechanism, the present model has a high computational efficiency due to its explicit and analytical expressions.
The model is validated by experimental results of the shear behavior of a fabric under different biaxial tensile loads. Furthermore, parametric studies are carried out using the present model. In general, fewer inflection points, a larger gap-to-width ratio, and softer in-plane compressive moduli decrease the shear stiffness of a fabric, which can improve the formability of a fabric. When tensile force along yarns is utilized to reduce wrinkling, one should consider that the tensile force will also affect the shear property of a fabric and will further affect the yarn re-orientation. With the help of the present analytical model, engineers can easily optimize the design of fabrics with different weave patterns. Because of the high computational efficiency of the model, the optimization process can be significantly improved.

Conclusions
An analytical model is proposed to describe the tension-shear coupling effect of fabrics with different weave patterns under large shear deformation. Different weave patterns are characterized by the number of inflection points of yarns. In this work, the tension-shear coupling mechanism can be revealed from two aspects, the friction dissipation energy and the in-plane transverse compressive dissipation energy. On the one hand, the subjected tensile force along yarns increases the contact force between overlapped yarns, leading to the growth of the friction dissipation energy. Then, larger friction dissipation energy brings more resistance with respect to shear deformation, thereby increasing the shear rigidity of a fabric. On the other hand, the tensile load along a yarn makes the yarn shrink. The shrinking improves the transverse compressive modulus of the yarn. Besides, the transverse compressive modulus of a yarn is also improved due to the compressive force between overlapped yarns induced by the tension force. Then, the dissipation energy of the transverse compression of yarns becomes larger due to the increase of the mechanical properties of yarns, thus enhancing the shear rigidity of a fabric. In addition to the description of the tension-shear coupling mechanism, the present model has a high computational efficiency due to its explicit and analytical expressions.
The model is validated by experimental results of the shear behavior of a fabric under different biaxial tensile loads. Furthermore, parametric studies are carried out using the present model. In general, fewer inflection points, a larger gap-to-width ratio, and softer in-plane compressive moduli decrease the shear stiffness of a fabric, which can improve the formability of a fabric. When tensile force along yarns is utilized to reduce wrinkling, one should consider that the tensile force will also affect the shear property of a fabric and will further affect the yarn re-orientation. With the help of the present analytical model, engineers can easily optimize the design of fabrics with different weave patterns. Because of the high computational efficiency of the model, the optimization process can be significantly improved.