Next Article in Journal
Synthesis of Calcium Peroxide Nanoparticles with Starch as a Stabilizer for the Degradation of Organic Dye in an Aqueous Solution
Previous Article in Journal
Design of New Dispersants Using Machine Learning and Visual Analytics
Previous Article in Special Issue
Recent Trends in Protective Textiles against Biological Threats: A Focus on Biological Warfare Agents
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Geometric Analysis of Three-Dimensional Woven Fabric with in-Plane Auxetic Behavior

1
School of Fashion and Textiles, The Hong Kong Polytechnic University, Hong Kong, China
2
Research Institute for Intelligent Wearable Systems, The Hong Kong Polytechnic University, Hong Kong, China
3
Department of Mechanical Engineering, Hakim Sabzevari University, Sabzevar 9617976487, Iran
*
Author to whom correspondence should be addressed.
Polymers 2023, 15(5), 1326; https://doi.org/10.3390/polym15051326
Submission received: 18 January 2023 / Revised: 1 March 2023 / Accepted: 2 March 2023 / Published: 6 March 2023

Abstract

:
Auxetic textiles are emerging as an enticing option for many advanced applications due to their unique deformation behavior under tensile loading. This study reports the geometrical analysis of three-dimensional (3D) auxetic woven structures based on semi-empirical equations. The 3D woven fabric was developed with a special geometrical arrangement of warp (multi-filament polyester), binding (polyester-wrapped polyurethane), and weft yarns (polyester-wrapped polyurethane) to achieve an auxetic effect. The auxetic geometry, the unit cell resembling a re-entrant hexagon, was modeled at the micro-level in terms of the yarn’s parameters. The geometrical model was used to establish a relationship between the Poisson’s ratio (PR) and the tensile strain when it was stretched along the warp direction. For validation of the model, the experimental results of the developed woven fabrics were correlated with the calculated results from the geometrical analysis. It was found that the calculated results were in good agreement with the experimental results. After experimental validation, the model was used to calculate and discuss critical parameters that affect the auxetic behavior of the structure. Thus, geometrical analysis is believed to be helpful in predicting the auxetic behavior of 3D woven fabrics with different structural parameters.

1. Introduction

Auxetic textile is a sub-class of meta-materials that are differentiated from conventional materials due to their negative Poisson’s ratio (NPR) [1]. So far, scientists have successfully developed fibers [2,3], polymers [4,5], yarns [6,7,8], and fabrics [9,10] with the uncommon property of an NPR. Textile fabrics with an NPR mean that they expand when they are subjected to tensile loading and shrink when compressive loading is applied. Due to this unusual response to the applied loading, auxetic fabrics show outstanding extensibility in both warp and weft directions, better drapability, and become more compact upon compression [1]. Therefore, auxetic fabrics are considered to be the first choice for applications in protective textiles, medical textiles, smart textiles, sportswear, and auxetic composites [11]. Researchers have conducted an immense amount of work to develop auxetic fabrics in every possible way. Therefore, auxetic fabrics could be realized in woven form two-dimensional (2D) auxetic woven fabrics [12,13,14], three-dimensional (3D) auxetic textiles [15,16,17], knitted form (warp-knitted auxetic fabrics [18,19], weft-knitted auxetic fabrics [20,21]), non-woven form [22,23], and auxetic laminated fabrics [24,25]. Mainly, there are two techniques for introducing auxetic property into textile fabrics. The first and more simple technique is to use auxetic yarns and a conventional weave pattern. This technique was first introduced by Miller et al. [6], and they used their invented double helix yarn (DHY) in weft direction of a plain-woven fabric. The reason for using the DHY in weft direction was to ensure the out-of-register position of each neighboring DHY strand to maximize the auxetic effect. When the fabric was subjected to tensile tension, the DHY strands overlap each other in the thickness direction, which causes an out-of-plane auxetic effect by increasing the thickness of the sample. In this study, the maximum achieved out-of-plane Poisson’s ratio (PR) was −0.1. Similarly, the same research group also developed auxetic plain-woven fabrics using helical auxetic yarn (HAY) in warp direction and multi-filament polyester yarn in weft direction [26]. This study concluded that the auxetic effect is associated with the weft yarn’s material and weave geometry. Recently, auxetic woven fabrics were produced using auxetic plied yarns (4 ply and 6 ply) and DHY [27]. For fabric production, auxetic yarns were used in the warp direction only, while conventional elastic yarn was used in the weft direction. Three weave designs (plain, twill, and satin) were selected, along with other design parameters of plied yarn, to evaluate their effect on the NPR of the woven fabrics. Although producing auxetic woven fabrics with this method is the easiest way, there are still many drawbacks that limit the practical application of such fabrics. Those drawbacks are: the NPR of auxetic yarn transfers partially to the fabric due to weaving constraints, the unstable auxetic effect, and the out-of-register orientation of auxetic yarns. Furthermore, the orientation of auxetic yarn inside the fabric structure is crucial and difficult, i.e., the yarns should be as straight as possible to obtain a better auxetic effect. Due to this limitation, this technique is only suitable for woven fabrics and cannot be applied to knitted fabrics.
The second technique of developing auxetic fabrics includes the use of non-auxetic yarns and the realization of special auxetic geometry. It is well known that the auxetic effect is associated with the geometrical arrangement of a structure’s unit cells such as re-entrant hexagonal, rotating squares, etc. [28]. Therefore, the idea of developing auxetic fabrics progressed when scientists gave considerable attention to configuring conventional yarns in a special geometrical arrangement using knitting or weaving technologies. Up to now, various weft-knitted [20] and warp-knitted [18] auxetic fabrics have been developed using this technique. However, the low structural stability and low elastic recovery of the knitted auxetic fabrics make them impractical for many applications such as protective textiles, etc. Recently, a considerable amount of work has been reported on the development of 2D auxetic woven fabrics. For the first time, a uni-stretch auxetic woven fabric was developed using conventional weaving techniques and non-auxetic yarns [9]. Auxetic behavior was introduced to fabrics through foldable structures, re-entrant hexagonal geometry, and a rotating rectangle structure. The basic mechanism involved in the auxetic effect is the different shrinkage properties of tight weaves (plain) and loose weaves (twill or satin), along with the use of elastic and non-elastic yarns. These fabrics showed a zero PR or a very low NPR in one direction only. Therefore, another study was conducted by Zulifqar et al. [29] on developing bi-stretch auxetic woven fabrics. The methodology of the study was based on a similar principle of differential shrinkage, but elastic yarns were used in both the warp and weft directions along with non-elastic yarns. Among all the developed auxetic structures based on different geometries, the zig-zag folded stripes showed the highest NPR of −0.15. The developed 2D auxetic woven fabrics have the potential to be used for clothing, fashion garments, etc. However, these fabrics cannot be used for high-performance applications because they experience more longitudinal deformation under tension, a low NPR, and a sharp decrease in the NPR under increased strain.
Nowadays, 3D auxetic fabrics have been the central focus for many researchers because they can exhibit more extraordinary mechanical performances compared to those of 2D auxetic fabrics. Ge and Hu [17], for the first time, reported a novel 3D auxetic fabric developed by combining knitting and non-weaving techniques. After fabrication, the fabric structure was converted into a soft composite for stability. Uniquely, this fabric structure was designed to show the NPR effect under compression. Hence, it was recommended for indentation-resistant and impact-resistant applications. Similarly, 3D multi-layer auxetic woven structures having an out-of-plane NPR were developed using the differential modulus of the yarns [30]. The out-of-plane NPR effect is triggered by the Z-direction binding yarns when they are stretched longitudinally. Due to the limitation of the out-of-plane auxetic effect, the authors used their developed fabrics as a reinforcement in polymer composite. Furthermore, the authors claimed improved impact-resistant properties for the 3D auxetic woven fabric reinforced composites. Not long ago, a new milestone was achieved by developing 3D auxetic woven fabrics with in-plane NPR effect [31]. The development is considered to be very significant in view of broader application aspects because it addresses the drawbacks associated with the existing 2D and 3D auxetic textiles. Furthermore, the fabrics that exhibit in-plane auxetic behavior at a higher tensile tension are useful for releasing the contact pressure when they are used in applications such as automotive seat belts, garment belts, etc. Although significant interests have been shown on the design, fabrication, and characteristics of 3D auxetic woven structure, the important information on the fundamental geometrical analytics of structure is yet lacking.
In this study, a geometrical model of 3D woven fabric has been proposed at the micro-level (yarn level) to establish a relationship between the geometrical parameters (in terms of yarn parameters) and the tensile deformation. Since the geometrical analysis is carried out at the micro-level, the geometrical arrangement and the physical properties of yarns are very crucial to facilitate the geometrical model. Thus, a 3D auxetic woven fabric was fabricated to carefully observe the arrangements of the warp and weft yarns, forming auxetic geometry. After observing the geometry in a free state, the fabric was then subject to tensile deformation to examine the behavior of the yarns together with the geometrical deformation of the whole structure. Here, it was found that the structure extends longitudinally and axially under tensile deformation due to its auxetic nature, which causes tension and compression to the weft yarns. Thus, by recognizing the geometry and physical behavior of yarns, semi-empirical equations were drawn based on the geometrical arrangement of the warp and weft yarns, including the tension and compression effect of the yarns. The study shows that the extracted results from the established semi-empirical equations are in good agreement with the experimental results. Therefore, it could be said that the proposed geometrical model will be a useful tool to design, predict, and optimize the auxetic behavior of the 3D auxetic woven fabric.

2. Methodology

A comprehensive methodology was adopted to propose an effective geometrical model for the 3D auxetic woven structure. First, the 3D auxetic woven fabric structure was designed and developed with appropriate warp and weft yarns. Then, the fabric was subjected to a tensile loading to observe its deformation behavior at the micro-level and macro-level. Finally, a geometrical model was proposed by establishing semi-empirical equations based on the warp and weft yarn’s geometry at the initial and deformed states.

2.1. Structure Design and Fabric Formation

A novel multi-layer 3D auxetic woven structure was previously designed by modifying the conventional 3D orthogonal through-the-thickness structure [31]. The auxetic effect was introduced by featuring an unusual lateral crimp to the warp yarns, which form a special geometry resembling a re-entrant hexagon, as shown in Figure 1a. The unusual lateral crimp of the warp yarns was induced by using two different weft yarn systems, i.e., a coarser, non-elastic yarn and a fine elastic yarn. A solid, coarser yarn was inserted with a 2/2 twill pattern as a binder (through the thickness), creating alternate empty spaces within the warp yarns. Meanwhile, multiple insertions of fine elastic yarn were made in stretched form, which shrank the structure in a relaxed state by forcing the warp yarns to crimp laterally, filling those empty spaces. As a result, an unusual lateral crimp, similar to a sinusoidal wave, was induced in the warp yarns. When this structure is stretched longitudinally, the warp yarns try to become straight, pushing the coarse binding yarns in a lateral direction. Thus, under tensile extension, the width of the structure also increases, and an auxetic effect is achieved.
For sample fabrication, the three types of polymer-based yarn were procured from Wai Hung Weaving Factory Limited with the desired properties. A single type of yarn made of multi-filament polyester with a 0.7 mm diameter was selected as a warp. Two types of yarn, namely coarse binding yarn and elastic yarn, made of polyester-wrapped polyurethane (PU) having diameters of 2.5 mm and 0.62 mm, respectively, were selected as the weft. In addition, other features of the warp yarn include stability and strength because it will bear tensile loading. Similarly, the coarse binding yarn was solid, but easier to bend to facilitate orthogonal shaping during weaving. A semi-automatic weaving machine with a maximum of 16 heald frames attached to the dobby shedding mechanism was used to produce the fabric sample. The required number of warp yarns were passed through 8 heald frames with a straight drawing-in draft to weave a two-layer 3D woven fabric. A special denting plan, i.e., two yarns per dent followed by an empty dent, was used. This denting sequence kept a gap between the warp yarns which facilitate the insertion of coarse binding yarn. Furthermore, elastic yarn was inserted in a stretch form to shrink the structure along the width direction when the sample was cut off the loom. The developed fabric sample is shown in Figure 1b.
It is worth mentioning that, in a previous study [31], the 3D auxetic woven structure was investigated with four highly influencing parameters: diameter of binding yarn, bending stiffness of binding yarn, repeats of elastic weft yarn, and stretch percentage of elastic weft yarn. Thereupon, it was concluded that the auxetic behavior of the 3D woven fabrics is primarily associated with the yarn’s properties that limit the weaving process. For example, the higher diameter and stiffness of binding yarn were difficult to weave to achieve a favorable 3D auxetic structure. Therefore, it is preferred to select and develop a fabric sample with the optimized parameters having the maximum auxetic effect for validating the geometric model of the structure because it is considered that a sample with maximum auxetic behavior indicates the best formation of auxetic geometry. The physical properties of the three types of yarn used to develop the sample are given in Table 1. The bending stiffness of the binding yarn was calculated with the beam deflection method reported by Msalilwa et al. [32].

2.2. Tensile Test and Measurement of Negative Poisson’s Ratio

To avoid possible slippage during the tensile test, aluminum tabs were mounted on both ends of each sample, covering the gripping area of the clamps. Epoxy EL2 and AT30 hardener were used to mount the aluminum tabs onto the specimen. The narrow-woven fabric samples of width 25 mm and length 150 mm were prepared according to the guidelines provided by tensile testing standard ASTM D5035-11 [33]. Further, the gauge length between the two jaws of the machine was kept at 75 mm. An Instron 5982 universal testing machine (UTM) with a maximum load capacity of 100 kN was used in this experiment. The samples were subjected to a strain-controlled tensile test, where the lower jaw was fixed, and the upper jaw moved at the rate of 30 mm per minute. The test was performed three times for each set of samples, and average values were reported.
The UTM can only provide raw data of the mechanical properties, therefore, a separate digital camera (Canon EOS 800D) was set up in front of the machine to record the change in dimensions of the fabric samples during the tensile deformation, as shown in Figure 2. Thus, each sample was marked with four dots (two longitudinally and two vertically) at a distance of 25 mm before clamping them onto the UTM. The video recording and the tensile test started simultaneously. First, the recorded videos were processed, and pictures were extracted at a rate of 1.5 s, which was equivalent to 1% tensile strain. Then, the extracted pictures were analyzed with screen ruler software to calculate the lateral and longitudinal deformation of the fabric. Finally, the PR for each sample was calculated using Equation (1).
υ = ε y ε x = - ( Y Y o ) / Y o ( X X o ) / X o = ( Y Y o ) X o ( X X o ) Y o
where εx is strain in longitudinal direction and εy is strain in the lateral direction; Y o and X o are the initial width and length as shown in Figure 1b; Y and X are width and length in the stretched state, respectively.

2.3. Yarns Analysis

In a previous study [31], it was found that the auxetic effect of the 3D woven fabric is highly associated with the extension of elastic weft yarn and compression of coarse binding yarn during the tensile deformation. Therefore, the inclusion of these two characteristics of yarns are essential for developing a geometrical model.
The extension of elastic weft yarn was simply characterized by performing a tensile test on a single yarn using the Instron 5566 machine. The yarn was fixed between the machine’s jaws, where the lower jaw wax fixed, and the upper jaw moved at a constant velocity of 30 mm/min. Similarly, the compression of coarse binding yarn was determined by performing a compression test. However, conventional flat compression could be applied here because the compression of the binding yarn was applied by the warp yarns. Hence, a customized assembly, as shown in Figure 3, was prepared to perform the test. The assembly was designed in such a way that the two pins, whose diameter was equal to that of the warp yarn, applied the compressive force on the binding yarn. So, the binding yarn was placed between the upper and lower parts of the assembly, perpendicular to the pins. Then, the assembly was placed between the two platens of the Instron 5566 machine. To perform the compression test, the upper platen was moved downward at a constant velocity of 1 mm/min, and the lower platen was fixed. A very sensitive load cell of capacity 50 N was used to obtain precise results.

2.4. Geometrical Analysis

The parametric terms used in this geometrical analysis are given below:
  • a o —initial diameter of binding yarn.
  • b o —diameter of warp yarn at the initial state and during the first stage of deformation.
  • b i i —diameter of warp yarn during the second stage of deformation.
  • l o —length of warp yarn between the centers of the two binding yarns at the initial state and during the first stage of deformation.
  • l i i —length of warp yarn between the centers of the two binding yarns during the second stage of deformation.
  • θ o —initial inclination angle of warp yarn.
  • θ i —inclination angle of warp yarn in a deformed state.
  • D o —sum of the radii of binding yarn and warp yarn at the initial state.
  • D i —sum of the radii of binding yarn and warp yarn during the first stage of deformation.
  • D i i —sum of the radii of binding yarn and warp yarn during the second stage of deformation.
  • a i —height of binding yarn after compression.
  • a i ' —width of binding yarn after compression.
A geometric model is proposed to predict the NPR effect of the 3D in-plane auxetic woven structure. The following assumptions are made based on the deformation behavior of the fabric during the tensile test:
  • Figure 4 represent the schematic diagram of the 3D auxetic woven structure while the detailed geometry of the unit cell of the structure is illustrated in Figure 5. It is assumed that all the repeating unit cells are similar in size and shape in the initial state and deformed symmetrically during extension.
  • In the geometrical model, the presentation of an elastic weft yarn has been omitted because the warp and binding yarn analysis is suitable enough to predict the auxetic behavior. However, the effect of the elastic weft yarn, i.e., causing compression to the binding yarn during tensile loading, was included.
  • The initial cross-section of the binding yarn is circular, which, during compression, changes to a race track shape due to the compression of the warp yarns caused by the restoring force of the elastic weft yarns. However, the cross-sectional area of the binding yarn remains constant all the time.
  • The warp yarns are crimped by the binding yarns at the initial state. It is assumed that the part of the warp yarn that is in contact with the binding yarn is circular, while the non-contact part is straight and tangent to the circular part.
  • There is no slippage at the contact points of binding and warp yarns during tensile deformation.
  • During tensile stretching, the structure deforms in two stages: decrimping (Figure 5a) and elongation (Figure 5b) of the warp yarns under tensile loading. Furthermore, it is also assumed that at the decrimping stage, the cross-section of binding yarn changes from a circular one to a race track one, while the diameter of warp yarn remains unchanged. However, in the elongation stage, there is no further change to the dimensions of the race-track-shaped binding yarn, while the diameter of the warp yarn changes due to elongation.
    • First stage: In the first stage of deformation, the warp yarns start decrimping at a constant rate and become fully straight at a particular tensile strain. In this stage, the length and diameter of the warp yarn remain constant, as shown in Figure 5a. Based on the above assumptions, the geometrical model of the first stage can be re-illustrated in detail, as shown in Figure 6.
      From Figure 6, the initial distance ( X o ) between the centers of two binding yarns is given by:
      X o = 2 p q ¯ + q s ¯  
      The length of p q ¯ and q s ¯ in terms of yarn’s parameters can be calculated by solving right angle triangles Δ o p q and Δ q s t , respectively.
      p q ¯ = D o sin θ o
      q s ¯ = q t ¯ cos θ o
      Whereas:
      q t ¯ = l o 2 q r ^
      q r ^ = D o θ o
      By solving Equations (4)–(6), the following relation can be obtained
      q s ¯ = ( l o 2 D o θ o ) cos θ o
      Substituting Equations (3) and (7) into Equation (2) gives the following relation
      X o = 2 D o sin θ o + ( l o 2 D o θ o ) cos θ o
      Likewise, the initial height ( Y o ) between the centers of two binding yarns is given as
      Y o = 2 o p ¯ s t ¯
      To solve Δ o p q and Δ q s t , the following relation can be made for o p ¯ and s t ¯
      o p ¯ = D o cos θ o
      s t ¯ = q t ¯ sin θ o
      Substituting Equations (5) and (6) into Equation (11) gives the following equation
      s t ¯ = ( l o 2 D o θ o ) sin θ o
      Substituting Equations (10) and (12) into Equation (9) gives the following relation
      Y o = 2 D o cos θ o ( l o 2 D o θ o ) sin θ o
      This relation can be used for calculating the initial height ( Y o ) using the inclination angle ( θ o ), however, Y o can also be calculated simply by using the following relation.
      Y o = a o 2 + b o
      Similarly, the distance between the centers of the two binding yarns at a deformed state ( X i ), is given below.
      X i = 2 p q ¯ + q s ¯ + 2 n o ¯
      Whereas:
      n o ¯ = ( a i ' a i ) / 2
      By solving the right triangles Δ o p q   and Δ q s t , the lengths p q ¯ and q s ¯ can be expressed with Equations (17) and (18), respectively.
      p q ¯ = D i sin θ i
      q s ¯ = [ l o 2 D i θ i ( a i ' a i ) ] cos θ i
      Substituting Equations (16)–(18) into Equation (15) gives the following relation
      X i = 2 D i sin θ i + [ l o 2 D i θ i ( a i ' a i ) ] cos   θ i + ( a i ' a i )
      The height, Y i ,   between the coarse binding yarns can express mathematically as,
      Y i = 2 o p ¯ s t ¯
      From Δ o p q and Δ q s t , o p ¯ and s t ¯ can be expressed in the following Equations
      o p ¯ = D i cos θ i
      s t ¯ = ( l o 2 D i θ i ( a i ' a i ) ) sin θ i
      Substituting Equations (21) and (22) into Equation (20) gives the following relation
      Y i = 2 D i cos θ i [ l o 2 D i θ i ( a i ' a i ) ] sin θ i
      As both the lateral and longitudinal strains in the initial and deformed states are known, the PR ( υ i ) for the first stage of deformation can be calculated. By putting Equations (8), (13), (19), and (23) into Equation (1), the following relation can be obtained
      υ i = 2 D i cos θ i l o 2 D i θ i a i + a i sin θ i 2 D o cos θ o + l o 2 D o θ o sin θ o × 2 D o sin θ o + l o 2 D o θ o cos θ o 2 D i sin θ i + l o 2 D i θ i a i + a i cos θ i + a i a i 2 D o sin θ o l o 2 D o θ o cos θ o × 2 D o cos θ o l o 2 D o θ o sin θ o
    • Second stage: As the warp yarns become decrimped and become fully straight in the first stage, they then elongate at a constant rate due to the applied tensile loading. In addition, the warp yarn’s diameter will also decrease due to elongation, as shown in Figure 7.
It should be noted that the height and width of the race-track-shaped binding yarns will remain constant, while according to assumption 6, the length and diameter of the warp yarn will change during the second stage of deformation. Therefore, the sum of the radii of warp yarn and binding yarn D i will change to D i i , which is given in Equation (25), and the length of warp yarn l o will change to l i i , which is written in Equation (26).
D i i = a i 2 + b i i 2
l i i = l o + Δ l
where Δ l is the change in the length of warp yarn, it can be calculated from the magnitude of tensile extension.
As the warp yarns are decrimped and straightened, the inclination angle ( θ i ) becomes “0”. By putting θ i = 0 into Equations (19) and (23), the following relations can be obtained, respectively
X i i = l i i
Y i i = 2 D i i
As the lateral and longitudinal strains at the second stage of deformation are known, the PR υ i i for the second stage of deformation can be calculated by putting Equations (8), (13), (27), and (28) into Equation (1)
υ i i = 2 D i i 2 D o cos θ o + l o 2 D o θ o sin θ o × 2 D o sin θ o + l o 2 D o θ o cos θ o l i i 2 D o sin θ o l o 2 D o θ o cos θ o × 2 D o cos θ o l o 2 D o θ o sin θ o .

3. Results and Discussion

After successfully developing the geometrical model, the first step is to validate it through experimental results. Hence, the model can be used to predict, present, and analyze critical findings based on different structural parameters of the 3D auxetic woven fabric.

3.1. Experimental Observation

3.1.1. Determining the Poisson’s Ratio during the First Stage of Deformation

To validate the geometrical model, the calculated tensile strain and PR results should be based on the same geometrical parameters as that of the developed 3D woven fabric sample. Considering that, all the parameters in the initial state are reasonably straightforward to obtain. For example, the radii of the coarse binding yarn ( a o 2 ) and warp yarn ( b o 2 ) can be taken from Table 1, the initial distance ( X o ) between two coarse binding yarns can be determined by dividing the total length of a sample by the total repeats of coarse binding yarn, and the initial height ( Y o ) can be calculated using Equation (14). When these parameters are known, the initial inclination angle ( θ o ) and length ( l o ) of the warp yarn can be calculated using Equations (8) and (13).
In contrast, other parameters, particularly the compression effect of coarse binding yarn in a deformed state, are comparatively complicated to calculate. According to assumption (3), the cross-section of binding yarn changes to a race track shape during tensile deformation due to the compression of warp yarns caused by the restoring force of the elastic weft yarns. Therefore, the magnitude of force applied by elastic weft yarn should be known. During tensile deformation, the force applied by elastic weft yarn increases because of the lateral expansion of the structure. So, the first step is to calculate the amount of lateral deformation of the structure without assuming the coarse binding yarn’s compression. For this purpose, Equations (8), (13), (19), and (23) can be used, while considering a i ' = a i = a o . When the percentage of lateral deformation of the structure is known, the magnitude of force applied by elastic weft yarn can be calculated from the force–strain curve of elastic weft yarn as shown in Figure 8. Eventually, the compression caused to the coarse binding yarn at a given force can be determined from the force–compressive extension curve of the coarse binding yarn, as shown in Figure 9.
After calculating the compression of coarse binding yarn in a deformed state, other parameters such as the final height ( a i )   and width ( a i ' ) of race-track-shaped binding yarn can be calculated based on assumption 3. Now, using Equation (19), the inclination angle ( θ i ) can be calculated at a certain tensile deformation ( X i ). Thus, all the variables are known, therefore, the PR ( υ i ) during the first stage of deformation can be determined using Equation (24). The calculated values of PR for the first stage of deformation are given in Table 2.

3.1.2. Determining the Poisson’s Ratio during the Second Stage of Deformation

The inclination angle ( θ i ) reaches 0 at a tensile strain of 8.07%, which indicates that the warp yarn is decrimped and becomes fully straight. At this point, the second stage of deformation starts, where no changes occur to the other parameters except to the length and diameter of the warp yarn. This specific point is determined as “critical strain,” and the PR is called “critical Poisson’s ratio”. After the critical strain, the warp yarns elongate at a rate of tensile strain, which is a known parameter. However, the decrease in the diameter of the warp yarn ( b i i ) at a particular tensile strain needs to be calculated. So, a customized tensile test was performed on the warp yarn, in which the cross head of UTM machine needed to be stopped after every 1% strain for 10 s to observe the diameter of the warp yarn using a thickness gauge. The variation trend of warp yarn’s diameter ( b i i ) as a function of tensile strain ( ε ) is shown in Figure 10. It can be observed that the change in b i i has a polynomial trend, therefore, according to a second-degree polynomial trend, the following equation can be established
b i i = m 1 ε i 2 + m 2 ε i + m 3
where ε i is the tensile strain, and m 1 , m 2 , and m 3 are the constants that can be determined from the experimental resutls of the warp yarn. Substituting Equation (30) in Equation (25) gives the following equation
D i i = a i 2 + m 1 ε 2 + m 2 ε + m 3 2
From Equation (31), a i is a known parameter, whereas the three other unknown constants were calculated experimentally from the relationship of warp yarn’s diameter and the tensile strain, which are: m 1 = 0.003 , m 2 = 0.0102 , and m 3 = 0.6608 .
When the sum of radii of the binding yarn and warp yarn D i i are known at a particular tensile extension, X i i , then the PR at the second stage of deformation can be calculated using Equation (29). Finally, the yarn’s parameters, the structural variables of the geometric model, and the calculated PR are given in Table 2.
It should be noted that if the tensile strain is less than the critical tensile strain, then υ i (Equation (24)) will be used to calculate the PR, otherwise υ i i Equation (29) will be used.
Figure 11 presents the PR curves plotted against the tensile strain of the developed 3D auxetic woven fabric and the calculated results obtained from the geometrical model. It can be seen that both the experimental and calculated PRs follow a similar trend. The agreement of calculated and experimental results confirms the reliability of the geometrical model. The PR first decreases, reaches a minimum at 8.07% tensile strain, and then follows a steadily increasing trend. From the analysis of the structure’s geometry, it is clear that the lateral crimp of warp yarn plays an important role in the auxetic behavior under the tensile extension. The warp yarn starts decrimp under the tensile deformation, and the PR decreases. The effect of decrimping continues until the warp yarn becomes fully straight, which takes effect at 8.07% tensile strain. Hereafter, there is no decrimping of warp yarn that indicates there will be no lateral extension with further tensile extension. In addition, the diameter of warp yarns decreases due to longitudinal extension. As a result, the PR increased beyond 8.07% tensile strain.

3.2. Prediction of Auxetic Behavior

3.2.1. Effect of Binding Yarn Compression

A previous study [31] claimed that the coarse binding yarn becomes compressed during tensile extension. From the experimental assessment of the structure, it is quite clear that compression occurs due to the restoring force of elastic yarn. However, what could not be answered is how much the compression of binding yarn effects the auxetic behavior. In comparison, the geometrical model can calculate the effect of binding yarn compression on the auxetic behavior of 3D woven fabric. Figure 12 shows the PR versus tensile strain curves. The two curves, with compression (W-C) and without compression (W/O-C), are calculated using the geometrical model by keeping all the parameters constant, except for the compression effect of the binding yarn. For WO-C, it is supposed that the diameter of the binding yarn ( a o ) does not change during the tensile strain, while for W-C, the diameter of the binding yarn decreases as the structure elongates in both the longitudinal and lateral directions. It is worth mentioning that the compression that occurred to the binding yarn is calculated experimentally, as explained in Section 3.1. From Figure 12, it can be observed that the auxetic behavior of W-C is 54% less than that of WO-C. This means that the binding yarn’s compressive stiffness is critical and must be considered actively when one is developing 3D auxetic woven fabric.

3.2.2. Effect of Binding Yarn Diameter

It is well known that [17] the diameters of two yarns (warp and binding yarns) significantly affect the NPR of 3D auxetic structures. Later on, it was proven again through an experimental study [31]. Since then, the effect of diameter of yarns on auxetic behavior has been understood, however, a scientific reason for the difference is still lacking. To provide a concrete reason, PR of the 3D woven structure was calculated using the geometrical model. Three different binding yarns diameter, i.e., 1.9 mm, 2.28 mm, and 2.75 mm, were used, while other parameters were kept constant. The PR–tensile strain curves of the 3D woven fabrics based on the different diameters of binding yarns are shown in Figure 13.
It can be found that the auxetic effect increases by increasing the diameter of the binding yarn, which validates the literature. For a solid reason, the numerical values generated by the geometrical model were analyzed. It is concluded that the initial lateral crimp percentage (%) of warp yarn decreases by increasing the diameter of the binding yarn. The crimp% can be calculated using Equation (32). For binding yarn diameters 1.9 mm, 2.28 mm, and 2.75 mm, the crimp% are 8.086, 8.011, and 7.949, respectively. Less lateral crimp in the warp yarn means the lateral expansion of auxetic structure will be triggered at a lower longitudinal strain, and the structure will reach its maximum lateral deformation under less strain. Since the longitudinal strain is inversely proportional to the NPR (Equation (1)), the auxetic behavior of the structure with the warp yarns being less laterally crimped is less notable, and vice versa.
C r i m p   % = l o X o X o × 100

3.2.3. Effect of Binding Yarn Spacing

Similar to the effect of binding yarn compression and binding yarn diameter, the spacing between the two binding yarns is also important to assess the auxetic behavior of the 3D auxetic woven structure. For this assessment, a similar methodological approach was applied, i.e., all the other parameters were kept constant, and the initial distance between the two binding yarns ( X o ), also referred to as binding yarn spacing, was changed for three measurements. The three selected values for X o were 2.8 mm, 3 mm, and 3.2 mm. The PR values were then calculated using the developed geometrical model that are presented in Figure 14.
It can be explained that the binding yarn spacing has a direct relationship with the NPR effect, i.e., the NPR increases by increasing the yarn spacing ( X o ) , and vice versa. After evaluating the measured values of each parameter, it was found that the binding yarn spacing also effected the lateral crimping of the warp yarn. Higher yarn spacing caused less crimping of the warp yarn, producing a high NPR effect at a lower longitudinal strain value. The calculated lateral crimp values are 9.74 mm, 8.01 mm, and 6.97 mm when X o are 2.8 mm, 3 mm, and 3.2 mm , respectively.

4. Conclusions

The 3D woven fabric was designed and fabricated with in-plane auxetic behavior. The auxetic effect was visualized by the unusual lateral crimping of warp yarns that form a special geometry, resembling a re-entrant hexagon. Based on the geometrical arrangement of yarns in the structure, a micro-level geometrical model was proposed in terms of the yarn parameters to establish a relationship between PR and tensile strain through semi-empirical equations. The following conclusions can be drawn according to the study.
  • As the geometrical model is based on the micro-level analysis of the structure, therefore, it is useful to predict the auxetic behavior with given yarn parameters. Additionally, it can help to find specific reasons for the NPR effect in relation to a particular parameter.
  • Good agreements are observed between the experimental results and the results obtained from the established empirical equations based on the geometrical analysis.
  • It was found from the geometric analysis that the compression that is applied to binding yarns during lateral expansion significantly affects the NPR of the structure, i.e., the NPR decreases by 62.5% after the compression of the coarse binding yarn. Therefore, the compressive stiffness of binding yarn should be considered carefully.
  • The diameter of binding yarn and the spacing between the two binding yarns affect the lateral crimp percentage of the warp yarn. Both the diameter and spacing of binding yarns have an inverse relationship with the crimp percentage of the warp yarn. For example, the PR of the samples are −2.1 and −3.0 when the binding yarn spacing is 2.8 mm and 3.2 mm, respectively.
  • The lateral crimp of warp yarn has a strong relationship with the NPR of the structure. The NPR increases if the crimping degree is lower and decreases if the crimping degree is higher. For example, when the lateral crimping of warp yarn is 8.09%, the PR is −2.0, while the PR for a lateral crimping of 7.95% is −3.78.
The retainability of the auxetic effect is of great significance to its practical application. Therefore, the auxetic behavior of the fabrics under repetitive tensile forces can be evaluated experimentally. Furthermore, the current study uses coarse yarns to produce the fabric. The auxetic fabric can be used for broader applications if the fabric is developed with finer yarns.

Author Contributions

M.Z.: Methodology, Data curation, Writing—original draft, Writing—review and editing. H.H.: Visualization, Supervision, Methodology, Resources, Project administration, Writing—review and editing. E.E.: Methodology, Validation, Writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Research Institute for Intelligent Wearable Systems of The Hong Kong Polytechnic University in form of an internal project (No. P0039471).

Acknowledgments

The authors would like to acknowledge the support from the Technician, Ng Shui-wing, of the weaving workshop at The Hong Kong Polytechnic University.

Conflicts of Interest

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

References

  1. Liu, Y.; Hu, H. A Review on Auxetic Structures and Polymeric Materials. Sci. Res. Essays 2010, 5, 1052–1063. [Google Scholar]
  2. Ravirala, N.; Alderson, K.L.; Davies, P.J.; Simkins, V.R.; Alderson, A. Negative Poisson’s Ratio Polyester Fibers. Text. Res. J. 2006, 76, 540–546. [Google Scholar] [CrossRef]
  3. Alderson, K.L.; Alderson, A.; Smart, G.; Simkins, V.R.; Davies, P.J. Auxetic Polypropylene Fibres:Part 1—Manufacture and Characterisation. Plast. Rubber Compos. 2002, 31, 344–349. [Google Scholar] [CrossRef]
  4. Alderson, K.; Evans, K. The Fabrication of Microporous Polyethylene Having a Negative Poisson’s Ratio. Polymer (Guildf.) 1992, 33, 4435–4438. [Google Scholar] [CrossRef]
  5. Caddock, B.D.; Evans, K.E. Microporous Materials with Negative Poisson’s Ratios. I. Microstructure and Mechanical Properties. J. Phys. D Appl. Phys. 1989, 22, 1877–1882. [Google Scholar] [CrossRef]
  6. Miller, W.; Hook, P.B.; Smith, C.W.; Wang, X.; Evans, K.E. The Manufacture and Characterisation of a Novel, Low Modulus, Negative Poisson’s Ratio Composite. Compos. Sci. Technol. 2009, 69, 651–655. [Google Scholar] [CrossRef]
  7. Gao, Y.; Chen, X.; Studd, R. Experimental and Numerical Study of Helical Auxetic Yarns. Text. Res. J. 2021, 91, 1290–1301. [Google Scholar] [CrossRef]
  8. Sloan, M.R.; Wright, J.R.; Evans, K.E. The Helical Auxetic Yarn—A Novel Structure for Composites and Textiles; Geometry, Manufacture and Mechanical Properties. Mech. Mater. 2011, 43, 476–486. [Google Scholar] [CrossRef]
  9. Zulifqar, A.; Hua, T.; Hu, H. Development of Uni-Stretch Woven Fabrics with Zero and Negative Poisson’s Ratio. Text. Res. J. 2018, 88, 2076–2092. [Google Scholar] [CrossRef]
  10. Ali, M.; Zeeshan, M.; Qadir, M.B.; Ahmad, S.; Nawab, Y.; Anjum, A.S.; Riaz, R. Development and Comfort Characterization of 2D Woven Auxetic Fabric for Wearable and Medical Textile Applications. J. Cloth. Text. 2018, 36, 199–214. [Google Scholar] [CrossRef]
  11. Wang, Z.; Hu, H. Auxetic Materials and Their Potential Applications in Textiles. Text. Res. J. 2014, 84, 1600–1611. [Google Scholar] [CrossRef]
  12. Cao, H.; Zulifqar, A.; Hua, T.; Hu, H. Bi-Stretch Auxetic Woven Fabrics Based on Foldable Geometry. Text. Res. J. 2019, 89, 2694–2712. [Google Scholar] [CrossRef]
  13. Ali, M.; Zeeshan, M.; Qadir, M.B.; Riaz, R.; Ahmad, S.; Nawab, Y.; Anjum, A.S. Development and Mechanical Characterization of Weave Design Based 2D Woven Auxetic Fabrics for Protective Textiles. Fibers Polym. 2018, 19, 2431–2438. [Google Scholar] [CrossRef]
  14. Zeeshan, M.; Ali, M.; Riaz, R.; Anjum, A.S.; Nawab, Y.; Qadir, M.B.; Ahmad, S. Optimizing the Auxetic Geometry Parameters in Few Yarns Based Auxetic Woven Fabrics for Enhanced Mechanical Properties Using Grey Relational Analysis. J. Nat. Fibers 2022, 19, 4594–4605. [Google Scholar] [CrossRef]
  15. Ahmed, H.I.; Umair, M.; Nawab, Y.; Hamdani, S.T.A. Development of 3D Auxetic Structures Using Para-Aramid and Ultra-High Molecular Weight Polyethylene Yarns. J. Text. Inst. 2021, 112, 1417–1427. [Google Scholar] [CrossRef]
  16. Khan, M.I.; Akram, J.; Umair, M.; Hamdani, S.T.; Shaker, K.; Nawab, Y.; Zeeshan, M. Development of Composites, Reinforced by Novel 3D Woven Orthogonal Fabrics with Enhanced Auxeticity. J. Ind. Text. 2019, 49, 676–690. [Google Scholar] [CrossRef]
  17. Ge, Z.; Hu, H. Innovative Three-Dimensional Fabric Structure with Negative Poisson’s Ratio for Composite Reinforcement. Text. Res. J. 2013, 83, 543–550. [Google Scholar] [CrossRef]
  18. Zhao, S.; Hu, H.; Kamrul, H.; Chang, Y.; Zhang, M. Development of Auxetic Warp Knitted Fabrics Based on Reentrant Geometry. Text. Res. J. 2020, 90, 344–356. [Google Scholar] [CrossRef]
  19. Wang, Z.; Hu, H. 3D Auxetic Warp-Knitted Spacer Fabrics. Phys. Status Solidi. 2014, 251, 281–288. [Google Scholar] [CrossRef]
  20. Glazzard, M.; Breedon, P. Weft-Knitted Auxetic Textile Design. Phys. Status Solidi 2014, 251, 267–272. [Google Scholar] [CrossRef]
  21. Boakye, A.; Chang, Y.; Rafiu, K.R.; Ma, P. Design and Manufacture of Knitted Tubular Fabric with Auxetic Effect. J. Text. Inst. 2018, 109, 596–602. [Google Scholar] [CrossRef]
  22. Bhullar, S.K.; Ko, J.; Cho, Y.; Jun, M.B.G. Fabrication and Characterization of Nonwoven Auxetic Polymer Stent. Polym. Plast. Technol. Eng. 2015, 54, 1553–1559. [Google Scholar] [CrossRef]
  23. Verma, P.; Shofner, M.L.; Lin, A.; Wagner, K.B.; Griffin, A.C. Inducing Out-of-Plane Auxetic Behavior in Needle-Punched Nonwovens. Phys. Status Solidi. 2015, 252, 1455–1464. [Google Scholar] [CrossRef]
  24. Yuping, C.; Liu, Y.; Hong, H. Deformation Behavior of Auxetic Laminated Fabrics with Rotating Square Geometry. Text. Res. J. 2022, 92, 4652–4665. [Google Scholar] [CrossRef]
  25. Shukla, S.; Behera, B.K. Auxetic Fibrous Structures and Their Composites: A Review. Compos. Struct. 2022, 290, 115530. [Google Scholar] [CrossRef]
  26. Wright, J.R.; Burns, M.K.; James, E.; Sloan, M.R.; Evans, K.E. On the Design and Characterisation of Low-Stiffness Auxetic Yarns and Fabrics. Text. Res. J. 2012, 82, 645–654. [Google Scholar] [CrossRef]
  27. Ng, W.; Hu, H. Woven Fabrics Made of Auxetic Plied Yarns. Polymers 2018, 10, 226. [Google Scholar] [CrossRef] [Green Version]
  28. Etemadi, E.; Gholikord, M.; Zeeshan, M.; Hu, H. Improved Mechanical Characteristics of New Auxetic Structures Based on Stretch-Dominated-Mechanism Deformation under Compressive and Tensile Loadings. Thin-Walled Struct. 2023, 184, 110491. [Google Scholar] [CrossRef]
  29. Zulifqar, A.; Hu, H. Development of Bi-Stretch Auxetic Woven Fabrics Based on Re-Entrant Hexagonal Geometry. Phys. Status Solidi. Basic Res. 2019, 256, 1800172. [Google Scholar] [CrossRef] [Green Version]
  30. Liaqat, M.; Samad, H.A.; Hamdani, S.T.A.; Nawab, Y. The Development of Novel Auxetic Woven Structure for Impact Applications. J. Text. Inst. 2017, 108, 1264–1270. [Google Scholar] [CrossRef]
  31. Zeeshan, M.; Hu, H.; Zulifqar, A. Three-Dimensional Narrow Woven Fabric with in-Plane Auxetic Behavior. Text. Res. J. 2022, 92, 4695–4708. [Google Scholar] [CrossRef]
  32. Msalilwa, L.R.; Kyosev, Y.; Rawal, A.; Kumar, U. Investigation of the Bending Rigidity of Double Braided Ropes. In Recent Developments in Braiding and Narrow Weaving; Kyosev, Y., Ed.; Springer International Publishing: Cham, Germany, 2016; pp. 47–57. [Google Scholar]
  33. ASTM D5035-11; Standard Test Method for Breaking Force and Elongation of Textile Fabrics (Strip Method). ASTM International: West Conshohocken, PA, USA, 2019.
Figure 1. Three dimensional auxetic woven structure: (a) schematic representation; (b) real fabric.
Figure 1. Three dimensional auxetic woven structure: (a) schematic representation; (b) real fabric.
Polymers 15 01326 g001
Figure 2. Testing setup on a universal testing machine for Poisson’s ratio measurement.
Figure 2. Testing setup on a universal testing machine for Poisson’s ratio measurement.
Polymers 15 01326 g002
Figure 3. Binding yarn compression test assembly.
Figure 3. Binding yarn compression test assembly.
Polymers 15 01326 g003
Figure 4. Schematic diagram of 3D auxetic woven structure: (a) top view; (b) top cross-sectional view [31].
Figure 4. Schematic diagram of 3D auxetic woven structure: (a) top view; (b) top cross-sectional view [31].
Polymers 15 01326 g004
Figure 5. Geometrical model of 3D auxetic structure: (a) decrimping of warp yarns under tensile load; (b) elongation of warp yarns under tensile load.
Figure 5. Geometrical model of 3D auxetic structure: (a) decrimping of warp yarns under tensile load; (b) elongation of warp yarns under tensile load.
Polymers 15 01326 g005
Figure 6. First stage of deformation of a unit-cell: (a) initial state; (b) deformed state.
Figure 6. First stage of deformation of a unit-cell: (a) initial state; (b) deformed state.
Polymers 15 01326 g006
Figure 7. Second stage of deformation of a unit cell: (a) decrimped and (b) elongation of warp yarn.
Figure 7. Second stage of deformation of a unit cell: (a) decrimped and (b) elongation of warp yarn.
Polymers 15 01326 g007
Figure 8. Tensile force–strain curve of elastic weft yarn.
Figure 8. Tensile force–strain curve of elastic weft yarn.
Polymers 15 01326 g008
Figure 9. Compressive force–extension curve of coarse binding yarn.
Figure 9. Compressive force–extension curve of coarse binding yarn.
Polymers 15 01326 g009
Figure 10. Trend of a diameter of warp yarn against tensile strain.
Figure 10. Trend of a diameter of warp yarn against tensile strain.
Polymers 15 01326 g010
Figure 11. Comparison of calculated and experimental results of Poisson’s ratio against tensile strain.
Figure 11. Comparison of calculated and experimental results of Poisson’s ratio against tensile strain.
Polymers 15 01326 g011
Figure 12. Effect of binding yarn’s compression on auxetic behavior of 3D woven structure.
Figure 12. Effect of binding yarn’s compression on auxetic behavior of 3D woven structure.
Polymers 15 01326 g012
Figure 13. Effect of binding yarn diameter on auxetic behavior of 3D woven structure.
Figure 13. Effect of binding yarn diameter on auxetic behavior of 3D woven structure.
Polymers 15 01326 g013
Figure 14. Effect of binding yarn spacing on auxetic behavior of 3D woven structure.
Figure 14. Effect of binding yarn spacing on auxetic behavior of 3D woven structure.
Polymers 15 01326 g014
Table 1. Properties of the yarns used for sample preparation.
Table 1. Properties of the yarns used for sample preparation.
YarnMaterialTypeDiameter
(mm)
Tensile Modulus
(MPa)
Tensile Strength (N)Bending Stiffness (×10−6 Nm2)
Warp yarnPolyester multi-filamentBraided yarn0.70153106-
Binding yarnPolyester-wrapped PUBraided yarn2.2882030.66
Elastic weft yarnPolyester-wrapped PUBraided yarn0.62314-
Table 2. Structural parameters and the calculated Poisson’s ratio results of 3D auxetic woven fabric.
Table 2. Structural parameters and the calculated Poisson’s ratio results of 3D auxetic woven fabric.
StateTensile Strain (%) a o   or   a i
(mm)
b o   or   b i i
(mm)
D o   or   D i   or   D i i
(mm)
X o   or   X i   or   X i i
(mm)
Y o   or   Y
(mm)
θ o   or   θ i
(rad)
l o   or   l i i
(mm)
PR
Initial state02.280.701.4903.0001.8400.4603.2420
Under tensile deformationFirst stage12.190.701.4473.0301.8450.4363.242−0.29
22.120.701.4103.0601.8620.4083.242−0.61
32.050.701.3773.0901.8900.3763.242−0.90
42.000.701.3513.1201.9330.3373.242−1.26
51.930.701.3163.1501.9740.2953.242−1.46
61.870.701.2863.1802.0370.2433.242−1.79
71.780.701.2453.2102.1080.1773.242−2.08
81.620.701.1583.2402.2260.0473.242−2.61
8.071.620.701.1133.2422.22603.242−2.59
Second stage91.620.621.1103.2702.22103.27−2.33
101.620.611.1083.3002.21603.3−2.07
111.620.601.1053.3302.21103.33−1.86
And so on
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zeeshan, M.; Hu, H.; Etemadi, E. Geometric Analysis of Three-Dimensional Woven Fabric with in-Plane Auxetic Behavior. Polymers 2023, 15, 1326. https://doi.org/10.3390/polym15051326

AMA Style

Zeeshan M, Hu H, Etemadi E. Geometric Analysis of Three-Dimensional Woven Fabric with in-Plane Auxetic Behavior. Polymers. 2023; 15(5):1326. https://doi.org/10.3390/polym15051326

Chicago/Turabian Style

Zeeshan, Muhammad, Hong Hu, and Ehsan Etemadi. 2023. "Geometric Analysis of Three-Dimensional Woven Fabric with in-Plane Auxetic Behavior" Polymers 15, no. 5: 1326. https://doi.org/10.3390/polym15051326

APA Style

Zeeshan, M., Hu, H., & Etemadi, E. (2023). Geometric Analysis of Three-Dimensional Woven Fabric with in-Plane Auxetic Behavior. Polymers, 15(5), 1326. https://doi.org/10.3390/polym15051326

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop