*Article* **Dynamic Equivalent Resistance Model of Knitted Strain Sensor under In-Plane and Three-Dimensional Surfaces Elongation**

**Yutian Li <sup>1</sup> , Pibo Ma <sup>2</sup> , Mingwei Tian 1,\* and Miao Yu 1,3,\***


**Abstract:** The dynamic equivalent resistance is a major index that determines the sensing performance of knitted strain sensors, and has the characteristics of in-plane and three-dimensional curved strain sensing. Therefore, in addition to establishing the in-plane equivalent resistance, it is necessary to establish a three-dimensional equivalent resistance model to fully explain the surface sensing performance. This project establishes two equivalent resistance models of knitted strain sensors under in-plane deformation and one equivalent resistance model of three-dimensional curved surface strain. Based on the length of resistance and the geometric topological structure, an in-plane strain macro– micro equivalent resistance model and a topological equivalent resistance model are established, respectively. In addition, a three-dimensional curved surface equivalent resistance model is created based on the volume resistance. By comparing the theoretical model with the experimental data, the results prove that the proposed in-plane and three-dimensional models can be utilized to calculate the resistance change of knitted strain sensors. Length resistance, coil transfer, and curved surface deformation depth are the main factors that affect the equivalent resistance of knitted strain sensors.

**Keywords:** knitted sensor; equivalent resistance; topology model; volume resistance; three-dimensional

### **1. Introduction**

The application of wearable strain sensors in detecting body motion posture, exercise rehabilitation, and healthcare makes our life much convenient [1–4]. Knitting fabrics, due to their multifunctional characteristics in terms of elasticity, flexibility, and deformability, are suitable to be employed in wearable garments as body monitoring sensors for movement [5–7]. Compared to woven and non-woven fabrics, knitting fabrics have unique properties [8–12]. Most textile-based sensors are based on the change in resistance [13–16], which can be used to predict fabric deformation and mechanical properties [17–19]. Therefore, there is a huge need to study the basic electro-mechanical properties of conductive knitted fabric, and develop the application of knitted strain sensors [20–24].

In the study of conductive fabric, some authors have studied the equivalent resistance model of textile sensors in their original state and under uniaxial stretching [22,25,26]; for example, Zhao et al. calculated and predicted the resistance of woven conductive fabric through the radius of the warp yarn and the resistance of one unit of conductive yarn [27]. Li et al. [28,29] explored the electro-mechanical properties of conductive yarn and conductive knitted fabric, with results given in terms of length-related resistance. It was shown that the resistance related to the length plays a dominant role in the total equivalent resistance [30]. However, the resistance change of the loop transfer was not mentioned [31]. Wang et al. examined the electro-mechanical property of conductive elastic knitted fabric based on a loops structure under biaxial extensions. Tokarska evaluated the planar anisotropy of conductive knitted fabrics and determined the resistance value

**Citation:** Li, Y.; Ma, P.; Tian, M.; Yu, M. Dynamic Equivalent Resistance Model of Knitted Strain Sensor under In-Plane and Three-Dimensional Surfaces Elongation. *Polymers* **2022**, *14*, 2839. https://doi.org/10.3390/ polym14142839

Academic Editors: Subhadip Mondal and Jeong In Han

Received: 2 May 2022 Accepted: 8 July 2022 Published: 12 July 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

on the surface of the knitted fabrics [32,33]. However, the fact is that in the process of stress deformation of conductive knitted fabric, which is affected by many factors such as environmental friction, both macroscopic fabric deformation and microscopic loop resistance changes should be considered.

The strain knitted sensor must consider the superimposed resistance problem of mutual contact between the coils in the initial state [34,35]. Therefore, the superimposed resistance of the loops squeezing against each other also plays an important role when the resistance changes [36–39]. Few studies can take into account the close contact resistance of the knitted loops in the initial state, and no appropriate theoretical model is given for exploring the equivalent resistance model of the knitted sensor under the initial state strain [40]. To explore the in-plane tensile properties, several papers have described corresponding equivalent resistance models under uniaxial and biaxial elongation and verified the validity of the corresponding model [41–45]. However, there are limited studies on the equivalent resistance model based on three-dimensional curved surface strain sensing.

This paper proposes two equivalent resistance models of knitted strain sensors under an in-plane deformation. They are macro–micro combined models based on length resistance and a topological equivalent resistance network based on segmented loop structure. In addition, a three-dimensional curved surface equivalent resistance model based on volume resistivity is offered. Comparing the theoretical model and experimental results, it is proven that in-plane and three-dimensional curved surface resistance models can be used to calculate the knitted strain sensor network. Key factors affecting the in-plane knitted strain sensing performance are the length resistance and the loop transfer, as shown in Figure 1. In three-dimensional surface sensing, the key factor influencing the sensing performance is the depth of the surface deformation.

**Figure 1.** The concepts of length resistance, loop transfer, and in-plane deformation.

### **2. Materials and Methods**

### *2.1. Materials*

The knitted strain sensor samples used in the experiments included a commercial silver-coated nylon filament yarn (40 dtex/12 f) purchased from the Tianyin Textile Technology Company (Qingdao, China). The silver-coated nylon yarn's normalized electrical resistance was 14 Ω/cm, with a resistivity of 0.1 ± 0.05 Ω·m and a tensile strength of 5.6 cN/dtex. Ordinary nylon filament yarn (75 dtex/36 f) and nylon-wrapped spandex yarn (50 dtex spandex filament covered by 20 dtex nylon staple fiber) were purchased from the Kejia Textile Fiber Products Company (Nantong, China). The ordinary nylon yarn was used as received for making a knitted base fabric, and was treated as an insulating material. A single circular knitting (SM8-TOP2 MP2, Pitch 0.907 mm, E28, Diameter 15 inch, Santoni Spa, Italy) machine was used for making conductive knitted fabric. The silver-coated yarn was knitted into the fabric as plating stitches. The specimen's off-loom course density was 130 wales/5 cm, and its wale density was 85 courses/5 cm. The size of the sensors was 5 cm × 1.5 cm (126 courses × 26 wales).

### *2.2. Sensing Performance Test*

The load and strain of the knitted strain sensor were measured by an Instron Model 5966, in which two pairs of clamps were used to fix the fabric in the wale direction. After selecting this machine, a sensor was fixed in the bottom carrier with a bursting tool, and

a ball was installed on the top to perform a curved strain sensing test on the sensor, as shown in Figure 2. The resistance variation of the knitted sensor was tested by a resistance acquisition program using a two-wire method to connect with the sample. This program was a real-time resistance acquisition program specially designed for knitted sensors.

**Figure 2.** Experimental equipment.

In accordance with the ASTM D5035 standard test method for the breaking force and elongation of textile fabrics, a specimen was stretched in the wale direction. In the experiments, the speed of the clamps was 100 mm/min, the gauge length was 100 mm, and the pre-load was 0.1 N. The relationship between strain and resistance was tested, and repeatability was tested with 2800 stretching cycles at 250 mm/min speed. The in-plane strain-resistance properties were tested by stretching the samples to produce a 100% strain 30 times under the same environment and calculating their resistance ratios, and the testing of the resistance change of the strain sensor was carried out in 2 courses × 27 wales with 50% stretch and calculating the ratios. Three-dimensional surface sensing performance was tested with a set of 400 mm/min for tensile rate and 60% strain.

### *2.3. Theoretical Model*

In-plane elongation is generally divided into two categories, course- and wale-stretch elongation [46]. The longitudinal elongation of the knitted fabric is greater than the transverse direction, and the longitudinal direction of the conductive fabric is used as the sensor's sensing direction.

Alternatively, the three-dimensional curved surface deformation can be approximated as a bursting experimental principle. Through the vertical drop of the sphere, the fabric is required to produce strains in all directions in the three-dimensional space. These two strain methods include the mechanical behavior of the sensor in various spatial directions, which is the basic premise for further research on the in-plane and three-dimensional curved surface strain sensing characteristics of knitted strain sensors.

The change of loop form is one of the main changing characteristics of a knitted strain sensor structure under both in-plane and three-dimensional stretching. The lengths of the loop segments change with loop transmission. When the arc of a parallel and sinker loop is transferred to a leg of loop under the longitudinal tensile, the loop section is transferred to the loop column part. Its tensile resistance deformation process is assumed to be the following steps:

(1) In the unstretched state, fabric loops are tightly abutted, and the resistance of the circuit is low. When subjected to longitudinal stretching, according to the law of resistance as given in Equation (1), the overall circuit resistance increases due to the length change.

$$\mathbf{R} = \frac{\rho \mathbf{L}}{\mathbf{S}} \tag{1}$$

where ρ is the proportionality constant and known as the resistivity or the specific resistance of the material of the conductor or substance; L is the length of the substance; S is the cross-sectional area of the substance.


$$\mathcal{R}\_{\mathbb{C}} = \frac{\mathfrak{p}}{2} \sqrt{\frac{\pi \mathcal{H}}{\mathfrak{np}}} \tag{2}$$

where ρ (Ωm−<sup>1</sup> ), H (Nm−<sup>2</sup> ), n, and p (N) are the electrical resistivity, material hardness, number of contact points, and contact pressure between the conductive yarn, respectively.

### *2.4. Model 1: Macro–Micro Equivalent Resistance Models Based on Length Resistance*

From a macro perspective, it was assumed that the conductive area of knitted fabric was a whole yarn, the overall length was L0, the cross-sectional area of the fabric was S0, and the resistance was R0. When the fabric was stretched and deformed, the overall length changed. Therefore, the length, cross-sectional area, and resistance were L1, S<sup>1</sup> and R1, respectively, as shown in Figure 3a. Then, from a microscopic point of view, there was resistance on each course of the loop. Taking 2 courses × 2 wales as an example, the arc of a parallel and leg of loop column constituted the first red line of course resistance (R1), the sinker loop and the second-course arc of a parallel/leg of loop column constituted a second yellow line of resistance (R2), and the second-course sinker loop constituted a third blue line of resistance (R3), as shown in Figure 3b. The combination of the red and blue line segments was the yellow loop, and the resistance was Equation (3). According to the 2 courses × 2 wales circuit networks and the series-parallel resistance calculation method, the equivalent resistance equation was Equation (4).

$$\mathbf{R\_2} = \mathbf{R\_1} + \mathbf{R\_3} \tag{3}$$

$$\frac{1}{\mathcal{R}\_{(2,2)}} = 2\left(\frac{1}{\mathcal{R}\_1} + \frac{1}{\mathcal{R}\_2} + \frac{1}{\mathcal{R}\_3}\right), \mathcal{R}\_{(2,2)} = 2\left(\frac{\mathcal{R}\_1 \,^2 + \mathcal{R}\_3 \,^2 + 3\mathcal{R}\_1 \mathcal{R}\_3}{\mathcal{R}\_1 \,^2 \mathcal{R}\_3 + \mathcal{R}\_3 \mathcal{R}\_1}\right) \tag{4}$$

**Figure 3.** (**a**) macro structure equivalent resistance model; (**b**) micro structure equivalent resistance model; (**c**) topological structure equivalent resistance model; (**d**) three-dimensional surface equivalent resistance model.

### *2.5. Model 2: The Equivalent Resistance of the Knitted Sensor Based on the Topology Model*

Figure 3c shows the loop structure of the knitted strain sensor, and the equivalent resistance is represented by the topology model. The knitting loop contained by the arc of parallel, sinker loop, and loop column. The arc of parallel had the same resistance as the sinking loop, the arc of parallel resistance was Ra, the sinking arc was Ra/2, and the loop column resistance was Rb. Taking the structure of 2 courses × 2 wales as an example, in the initial state, the conductive knitted fabric was under the influence of spandex, the loops were in close contact, and the loop column and the loop column and the settlement arc and the settlement arc were close to each other. Therefore, the mutual contact resistance of the loop columns was 2Rb, and the mutual contact resistance of the sinking loop was Ra.

### *2.6. Model 3: The Equivalent Resistance of the Knitted Sensor Based on the Topology Model*

Depending on the characteristics of the three-dimensional curved surface strain sensor, the sensors strain surface area change was approximated as a frustum of a cone. As shown in Figure 3d, in the burst test, R was the upper-bottom surface cone radius, r was the bottom surface cone radius, and l was the generatrix of a cone. Therefore, the three-dimensional surface resistance was calculated according to the volume resistivity formula, which refers to material resistance of current per unit volume, and is used to characterize material electrical properties.

$$\mathbf{R\_{V}} = \frac{\rho \mathbf{D}}{\mathbf{S\_{V}}} \tag{5}$$

where R<sup>v</sup> is the three-dimensional curved surface deformation resistance value, S<sup>v</sup> is the curved surface deformation area, D is the sensor's burst depth, and ρ is volume resistivity in Ω/cm.

### **3. Results**

### *3.1. Calculation of the Macro–Micro Equivalent Resistance Model Based on Length Resistance*

Calculated from the macrostructure at a relaxed state (according to the law of resistance), when the fabric is in its unstretched state, its length, cross-sectional area, and resistance are L0, S0, and R0, respectively. The strain rate is ε for the in-plane stretching. Its length, cross-sectional area, and resistance are L1, S1, and R<sup>1</sup> after the in-plane stretching, as shown in Equation (6). Therefore, the theoretical model of the equivalent resistance ratio of the in-plane strain can be deduced from Equation (7). The calculation resistance ratio R <sup>0</sup>/R is assumed to be y<sup>1</sup> , where the relationship between y<sup>1</sup> and strain ε is:

$$\begin{cases} \begin{aligned} \mathcal{R}\_{0} &= \frac{\rho \mathcal{L}\_{0}}{\mathcal{S}\_{0}} = \frac{\rho \mathcal{L}\_{0}}{\pi \mathbf{r}^{2}}\\ \mathcal{L}\_{1} &= \varepsilon \mathcal{L}\_{0} + \mathcal{L}\_{0} \\ \mathcal{S}\_{1} &= \pi \left(\frac{\mathbf{r}}{\varepsilon + 1}\right)^{2} \\ \mathcal{R}\_{1} &= \frac{\rho \left(\varepsilon + 1\right) \mathcal{L}\_{0}}{\pi \left(\frac{\mathbf{r}}{\varepsilon + 1}\right)^{2}} = \frac{\rho \mathcal{L}\_{0} \left(\varepsilon + 1\right)^{3}}{\pi \mathbf{r}^{2}} \\ \mathcal{R}\_{0} &= \varepsilon^{3} + 3\varepsilon^{2} + 3\varepsilon + 1 \end{aligned} \tag{6}$$
 
$$\mathcal{Z}\_{1} = \mathcal{Z}\_{2}$$

$$\mathbf{y}\_1 = \varepsilon^3 + \mathbf{3}\varepsilon^2 + \mathbf{3}\varepsilon + \mathbf{1} \tag{7}$$

Calculated from the microstructure, according to the 2 courses × 2 wales equivalent circuit and by using the same method, Equation (8) shows that in the relationship among the equivalent resistances of 2 courses × 1 wale, 2 courses × 3 wales, 3 courses × 2 wales, and 4 courses × 2 wales, the series-parallel circuit's connection of equivalent resistance of knitting sensor exists in course and wale direction. Therefore, the M course × N wale equivalent resistance of the knitted sensor can be deduced from Equation (9); it should

be noted that M and N refer to the number of loop courses and wales in the sensing area, respectively, as shown in Figure 4.

 1 R(2,1) = <sup>1</sup> R1 + <sup>1</sup> R2 + <sup>1</sup> R3 , R(2,1) = R1 <sup>2</sup>+R<sup>3</sup> <sup>2</sup>+3R1R<sup>3</sup> R1 <sup>2</sup>R3+R<sup>3</sup> <sup>2</sup>R<sup>1</sup> 1 R(2,2) = 2 1 R1 + <sup>1</sup> R2 + <sup>1</sup> R3 , R(2,2) = 2 R1 <sup>2</sup>+R<sup>3</sup> <sup>2</sup>+3R1R<sup>3</sup> R1 <sup>2</sup>R3+R<sup>3</sup> <sup>2</sup>R<sup>1</sup> 1 R(2,3) = 3 1 R1 + <sup>1</sup> R2 + <sup>1</sup> R3 , R(2,3) = 3 R1 <sup>2</sup>+R<sup>3</sup> <sup>2</sup>+3R1R<sup>3</sup> R1 <sup>2</sup>R3+R<sup>3</sup> <sup>2</sup>R<sup>1</sup> 1 R(3,2) = 2 1 R1 + <sup>1</sup> R2 + <sup>1</sup> R2 + <sup>1</sup> R3 , R(3,2) = 2 R1 <sup>2</sup>+R<sup>3</sup> <sup>2</sup>+4R1R<sup>3</sup> R1 <sup>2</sup>R3+R<sup>3</sup> <sup>2</sup>R<sup>1</sup> 1 R(4,2) = 4 1 R1 + <sup>1</sup> R2 + <sup>1</sup> R2 + <sup>1</sup> R2 + <sup>1</sup> R3 , R(4,2) = 4 R1 <sup>2</sup>+R<sup>3</sup> <sup>2</sup>+5R1R<sup>3</sup> R1 <sup>2</sup>R3+R<sup>3</sup> <sup>2</sup>R<sup>1</sup> (8)

$$\frac{1}{\mathcal{R}\_{\text{(M,N)}}} = \mathcal{N}\left(\frac{1}{\mathcal{R}\_1} + (\mathcal{M} - 1)\frac{1}{\mathcal{R}\_2} + \frac{1}{\mathcal{R}\_3}\right), \mathcal{R}\_{\text{(M,N)}} = \mathcal{N}\left(\frac{\mathcal{R}\_1^2 + \mathcal{R}\_3^2 + (\mathcal{M} + 1)\mathcal{R}\_1\mathcal{R}\_3}{\mathcal{R}\_1^2\mathcal{R}\_3 + \mathcal{R}\_3^2\mathcal{R}\_1}\right) (\mathcal{M} \mathbf{>1}) \tag{9}$$

**Figure 4.** (**a**) 2 courses × 1 wale, 2 courses × 3 wales, 3 courses × 2 wales, and 4 courses × 2 wales equivalent circuits; (**b**) The equivalent circuit of M course × N wale.

When knitted fabrics are in-plane stretched an extending of each loop and yarn segment occurs, transferring among the arc of parallel, leg and sinker loop. Then, the adjacent loops become tightly interlinked with each other. Due to the large strain caused by the range of extension, the following assumptions are considered to simplify the structure of knitted loop under a large strain. Assuming that the resistance changes of R1, R2, and R<sup>3</sup> are approximately equal (for the simple calculation, the three resistances are approximately equal to R2), then the total resistance (R) of the network and the resistance after stretching (R 0 ) can be calculated as in Equation (10). The equivalent resistance ratio R 0/R is assumed to be y<sup>2</sup> , where the relationship between y<sup>2</sup> , strain ε, course, and wale is presented as Equation (11).

$$\begin{cases} \frac{1}{\mathcal{R}} = \frac{1}{\mathcal{N}} \left( \frac{1}{\mathcal{R}\_2} + \frac{1}{\mathcal{R}\_2} + \dots + \frac{\mathcal{M}}{\mathcal{R}\_2} \right) \\ \mathcal{R} = \mathcal{N} \frac{\mathcal{R}\_2}{\mathcal{M} + 1} = \mathcal{N} \frac{\frac{\rho \mathcal{L}}{\mathcal{S}}}{\mathcal{M} + 1} \\ \mathcal{L}' = (\varepsilon + 1) L\_\prime \ S' = \pi \left( \frac{\mathbf{r}}{\varepsilon + 1} \right)^2 \\ \mathcal{R}' = \frac{\rho \mathcal{L}'}{\mathcal{S}'} = \frac{\rho 2 \left( \varepsilon + 1 \right) \mathcal{L}}{\pi \left( \frac{\mathbf{r}}{\varepsilon + 1} \right)^2} \\ \frac{\mathcal{R}'}{\mathcal{R}} = \mathcal{N} \frac{2 \left( \varepsilon + 1 \right)^3}{\mathcal{M} + 1} \end{cases} \tag{10}$$

$$\mathbf{y\_2 = N} \frac{2(\varepsilon + 1)^3}{\mathbf{M} + 1} \tag{11}$$

### *3.2. Calculation of the Equivalent Resistance Based on the Topology Model*

According to Kirchhoff's current and voltage law, the equivalent resistance of the topology model and the corresponding equations can be calculated [47]. The total current (I1) in the circuit can be calculated by Matlab, and the equivalent resistance (R) of the knitted strain sensor can be defined as Equation (12), where V is the total voltage that the circuit loads; (I1) is the total current in the circuit.

$$\mathbf{R} = \frac{\mathbf{V}}{\mathbf{I}\_1} \tag{12}$$

Based on the geometric structure of the knitting loop, the knitting strain sensor starts with 2 courses × 2 wales and extends to M course × N wale. Kirchhoff's current and voltage law is used to calculate the total circuit I<sup>1</sup> of the knitted fabric circuit. Figure 5a is a loop circuit with 2 courses × 2 wales. The knitting sensor is connected to the power supply at both ends, and its voltage is V. The circuit equations calculated by Kirchhoff's current and voltage law are shown in Equation (13), and this circuit contains 5 loop currents. Solving the above equation with Matlab comprehensively calculates the equivalent resistance of the 2 courses × 2 wales loop circuits as Equation (14).

**Figure 5.** (**a**) 2 courses × 2 wales loop circuit; (**b**) 3 courses × 2 wales and 4 courses × 2 wales loop circuits.

The same method can be used to solve for the equivalent resistances of the 3 courses × 2 wales loop circuits and 4 courses × 2 wales loop circuits, with the circuit equation set calculated by Kirchhoff's current and voltage laws, as shown in Figure 5b and the Equations (15) and (16). The 3 courses × 2 wales circuits contain seven loop circuits, and the 4 courses × 2 wales circuits contain nine loop currents.

$$\begin{cases} 2\mathbf{I}\_{1}\mathbf{R}\_{\mathbf{a}} - \mathbf{I}\_{2}\mathbf{R}\_{\mathbf{a}} - \mathbf{I}\_{3}\mathbf{R}\_{\mathbf{a}} = V\\\ \mathbf{I}\_{2}(2\mathbf{R}\_{\mathbf{a}} + 3\mathbf{R}\_{\mathbf{b}}) - \mathbf{I}\_{1}\mathbf{R}\_{\mathbf{a}} - \mathbf{I}\_{4}\mathbf{R}\_{\mathbf{a}} - 2\mathbf{I}\_{3}\mathbf{R}\_{\mathbf{b}} = 0\\\ \mathbf{I}\_{3}(2\mathbf{R}\_{\mathbf{a}} + 3\mathbf{R}\_{\mathbf{b}}) - \mathbf{I}\_{1}\mathbf{R}\_{\mathbf{a}} - \mathbf{I}\_{5}\mathbf{R}\_{\mathbf{a}} - 2\mathbf{I}\_{2}\mathbf{R}\_{\mathbf{b}} = 0\\\ \mathbf{I}\_{4}(2\mathbf{R}\_{\mathbf{a}} + 3\mathbf{R}\_{\mathbf{b}}) - 2\mathbf{I}\_{5}\mathbf{R}\_{\mathbf{b}} - \mathbf{I}\_{2}\mathbf{R}\_{\mathbf{a}} = 0\\\ \mathbf{I}\_{5}(2\mathbf{R}\_{\mathbf{a}} + 3\mathbf{R}\_{\mathbf{b}}) - \mathbf{I}\_{3}\mathbf{R}\_{\mathbf{a}} - 2\mathbf{I}\_{4}\mathbf{R}\_{\mathbf{b}} = 0\end{cases} \tag{13}$$

$$\mathbf{R}\_{(2,2)} = \frac{2\mathbf{R}\_{\mathbf{a}} \left(\mathbf{R}\_{\mathbf{a}}\,^2 + 3\mathbf{R}\_{\mathbf{a}}\mathbf{R}\_{\mathbf{b}} + \mathbf{R}\_{\mathbf{b}}\,^2\right)}{3\mathbf{R}\_{\mathbf{a}}\,^2 + \mathbf{R}\_{\mathbf{b}}\,^2 + 4\mathbf{R}\_{\mathbf{a}}\mathbf{R}\_{\mathbf{b}}}\tag{14}$$

 2I1R<sup>a</sup> − I2R<sup>a</sup> − I3R<sup>a</sup> = *V* I2(2R<sup>a</sup> + 3Rb) − I1R<sup>a</sup> − I4R<sup>a</sup> − 2I3R<sup>b</sup> = 0 I3(2R<sup>a</sup> + 3Rb) − I1R<sup>a</sup> − I5R<sup>a</sup> − 2I2R<sup>b</sup> = 0 I4(2R<sup>a</sup> + 3Rb) − I2R<sup>a</sup> − I6R<sup>a</sup> − 2I5R<sup>b</sup> = 0 I5(2R<sup>a</sup> + 3Rb) − I3R<sup>a</sup> − I7R<sup>a</sup> − 2I4R<sup>b</sup> = 0 I6(2R<sup>a</sup> + 3Rb) − I4R<sup>a</sup> − 2I7R<sup>b</sup> = 0 I7(2R<sup>a</sup> + 3Rb) − I5R<sup>a</sup> − 2I6R<sup>b</sup> = 0 (15) 2I1R<sup>a</sup> − I2R<sup>a</sup> − I3R<sup>a</sup> = *V* I2(2R<sup>a</sup> + 3Rb) − I1R<sup>a</sup> − I4R<sup>a</sup> − 2I3R<sup>b</sup> = 0 I3(2R<sup>a</sup> + 3Rb) − I1R<sup>a</sup> − I5R<sup>a</sup> − 2I2R<sup>b</sup> = 0 I4(2R<sup>a</sup> + 3Rb) − I2R<sup>a</sup> − I6R<sup>a</sup> − 2I5R<sup>b</sup> = 0 I5(2R<sup>a</sup> + 3Rb) − I3R<sup>a</sup> − I7R<sup>a</sup> − 2I4R<sup>b</sup> = 0 I6(2R<sup>a</sup> + 3Rb) − I4R<sup>a</sup> − I8R<sup>a</sup> − 2I7R<sup>b</sup> = 0 I7(2R<sup>a</sup> + 3Rb) − I5R<sup>a</sup> − I9R<sup>a</sup> − 2I6R<sup>b</sup> = 0 I8(2R<sup>a</sup> + 3Rb) − I6R<sup>a</sup> − 2I9R<sup>b</sup> = 0 I9(2R<sup>a</sup> + 3Rb) − I7R<sup>a</sup> − 2I8R<sup>b</sup> = 0 (16)

The equivalent resistance of the 3 courses × 2 wales and 4 courses × 2 wales loop circuits can be comprehensively calculated by solving the above equation with Matlab, as in the Equations (17) and (18). The mathematical transformation of the equivalent resistance Equations (14), (17) and (18) for 2 courses × 2 wales, 3 courses × 2 wales, and 4 courses × 2 wales results in the following formulas being obtained:

$$\mathbf{R}\_{(3,2)} = \frac{2\mathbf{R}\_{\mathbf{a}} \left(\mathbf{R}\_{\mathbf{a}}{}^{3} + 6\mathbf{R}\_{\mathbf{a}}{}^{2}\mathbf{R}\_{\mathbf{b}} + 5\mathbf{R}\_{\mathbf{a}}\mathbf{R}\_{\mathbf{b}}{}^{2} + \mathbf{R}\_{\mathbf{b}}{}^{3}\right)}{4\mathbf{R}\_{\mathbf{a}}{}^{3} + \mathbf{R}\_{\mathbf{b}}{}^{3} + 10\mathbf{R}\_{\mathbf{a}}{}^{2}\mathbf{R}\_{\mathbf{b}} + 6\mathbf{R}\_{\mathbf{a}}\mathbf{R}\_{\mathbf{b}}{}^{2}} \tag{17}$$

$$\mathbf{R}\_{(4,2)} = \frac{2\mathbf{R}\_{\mathbf{a}} \left( \mathbf{R}\_{\mathbf{a}}{}^{4} + 10\mathbf{R}\_{\mathbf{a}}{}^{3}\mathbf{R}\_{\mathbf{b}} + 15\mathbf{R}\_{\mathbf{a}}{}^{2}\mathbf{R}\_{\mathbf{b}}{}^{2} + 7\mathbf{R}\_{\mathbf{a}}\mathbf{R}\_{\mathbf{b}}{}^{3} + \mathbf{R}\_{\mathbf{b}}{}^{4} \right)}{5\mathbf{R}\_{\mathbf{a}}{}^{4} + \mathbf{R}\_{\mathbf{b}}{}^{4} + 20\mathbf{R}\_{\mathbf{a}}{}^{3}\mathbf{R}\_{\mathbf{b}} + 21\mathbf{R}\_{\mathbf{a}}{}^{2}\mathbf{R}\_{\mathbf{b}}{}^{2} + 8\mathbf{R}\_{\mathbf{a}}\mathbf{R}\_{\mathbf{b}}{}^{3}} \tag{18}$$

 R(2,2) = 2Ra[(Ra+Rb) <sup>2</sup>+RaRb] (Ra+Rb) <sup>2</sup>+2Ra(Ra+Rb) = (R<sup>a</sup> + Rb) − Ra <sup>3</sup>+R<sup>a</sup> <sup>2</sup>Rb+3RaR<sup>b</sup> <sup>2</sup>+R<sup>b</sup> 3 3Ra <sup>2</sup>+4RaRb+R<sup>b</sup> 2 R(3,2) = <sup>2</sup> 3 (R<sup>a</sup> + Rb) − 2(Ra <sup>4</sup>−4R<sup>a</sup> <sup>3</sup>Rb+R<sup>a</sup> <sup>2</sup>R<sup>b</sup> <sup>2</sup>+4RaR<sup>b</sup> <sup>3</sup>+R<sup>b</sup> 4 ) 3(4Ra <sup>3</sup>+10R<sup>a</sup> <sup>2</sup>Rb+6RaR<sup>b</sup> <sup>2</sup>+R<sup>b</sup> 3) R(4,2) = <sup>1</sup> 2 (R<sup>a</sup> + Rb) − Ra <sup>5</sup>−15R<sup>a</sup> <sup>4</sup>Rb−19R<sup>a</sup> <sup>3</sup>R<sup>b</sup> <sup>2</sup>+R<sup>a</sup> <sup>2</sup>R<sup>b</sup> <sup>3</sup>+5RaR<sup>b</sup> <sup>4</sup>+R<sup>b</sup> 5 2(5R<sup>a</sup> <sup>4</sup>+20R<sup>a</sup> <sup>3</sup>Rb+21R<sup>a</sup> <sup>2</sup>R<sup>b</sup> <sup>2</sup>+8RaR<sup>b</sup> <sup>3</sup>+R<sup>b</sup> 4) (19)

It is impossible to find the rule of M course × 2 wales from this expression. Therefore, if Equation (19) is factorized, the denominator of the last term in the equation is much larger than the numerator. When there are more courses, the last term in the three expressions is approximately zero and can be ignored. Thus, the following relationship is obtained:

$$\begin{cases} \begin{array}{l} \mathcal{R}\_{(2,2)} \approx (\mathcal{R}\_{\mathbf{a}} + \mathcal{R}\_{\mathbf{b}})\\\mathcal{R}\_{(3,2)} \approx \frac{2}{3}(\mathcal{R}\_{\mathbf{a}} + \mathcal{R}\_{\mathbf{b}})\\\mathcal{R}\_{(4,2)} \approx \frac{1}{2}(\mathcal{R}\_{\mathbf{a}} + \mathcal{R}\_{\mathbf{b}})\\\mathcal{R}\_{(\mathbf{M},2)} \approx \frac{2}{3}(\mathcal{R}\_{\mathbf{a}} + \mathcal{R}\_{\mathbf{b}}) \end{array} \tag{20}$$

It is inferred that along the wale direction of the knitted sensor, the equivalent resistance decreases with the increase of courses. The circuit is a parallel circuit, and the equivalent resistance expression of M course × 2 wales is Equation (21). Figure 6 shows the equivalent resistances of the circuits with 2 courses × 3 wales, 2 courses × 4 wales, and 2 courses × 5 wales. The equivalent resistances are calculated by the same method, and the equations obtained are shown in Equations (22)–(24).

$$\mathbf{R}\_{\text{(M,2)}} = \frac{2}{\mathbf{M}} (\mathbf{R}\_{\text{\(\mathbf{a}\)}} + \mathbf{R}\_{\text{b}}) \tag{21}$$

 3I1R<sup>a</sup> − I2R<sup>a</sup> − I3R<sup>a</sup> − I4R<sup>a</sup> = *V* I2(2R<sup>a</sup> + 3Rb) − I1R<sup>a</sup> − 2I3R<sup>b</sup> − I5R<sup>a</sup> = 0 I3(2R<sup>a</sup> + 4Rb) − I1R<sup>a</sup> − 2I2R<sup>b</sup> − I6R<sup>a</sup> − 2I4R<sup>b</sup> = 0 I4(2R<sup>a</sup> + 3Rb) − 2I3R<sup>b</sup> − I1R<sup>a</sup> − I7R<sup>a</sup> = 0 I5(2R<sup>a</sup> + 3Rb) − I2R<sup>a</sup> − 2I6R<sup>b</sup> = 0 I6(2R<sup>a</sup> + 4Rb) − I3R<sup>a</sup> − 2I5R<sup>b</sup> − 2I7R<sup>b</sup> = 0 I7(2R<sup>a</sup> + 3Rb) − I4R<sup>a</sup> − 2I6R<sup>b</sup> = 0 (22) 4I1R<sup>a</sup> − I2R<sup>a</sup> − I3R<sup>a</sup> − I4R<sup>a</sup> − I5R<sup>a</sup> = *V* I2(2R<sup>a</sup> + 3Rb) − I1R<sup>a</sup> − 2I3R<sup>b</sup> − I6R<sup>a</sup> = 0 I3(2R<sup>a</sup> + 4Rb) − I1R<sup>a</sup> − 2I2R<sup>b</sup> − I7R<sup>a</sup> − 2I4R<sup>b</sup> = 0 I4(2R<sup>a</sup> + 4Rb) − 2I3R<sup>b</sup> − 2I5R<sup>b</sup> − I1R<sup>a</sup> − I8R<sup>a</sup> = 0 I5(2R<sup>a</sup> + 3Rb) − I1R<sup>a</sup> − I9R<sup>a</sup> − 2I4R<sup>b</sup> = 0 I6(2R<sup>a</sup> + 3Rb) − I2R<sup>a</sup> − 2I7R<sup>b</sup> = 0 I7(2R<sup>a</sup> + 4Rb) − I3R<sup>a</sup> − 2I6R<sup>b</sup> − 2I8R<sup>b</sup> = 0 I8(2R<sup>a</sup> + 4Rb) − I4R<sup>a</sup> − 2I7R<sup>b</sup> − 2I9R<sup>b</sup> = 0 I9(2R<sup>a</sup> + 3Rb) − I5R<sup>a</sup> − 2I8R<sup>b</sup> = 0 (23) 5I1R<sup>a</sup> − I2R<sup>a</sup> − I3R<sup>a</sup> − I4R<sup>a</sup> − I5R<sup>a</sup> − I6R<sup>a</sup> = *V* I2(2R<sup>a</sup> + 3Rb) − I1R<sup>a</sup> − I32R<sup>b</sup> − I7R<sup>a</sup> = 0 I3(2R<sup>a</sup> + 4Rb) − I1R<sup>a</sup> − I8R<sup>a</sup> − I22R<sup>b</sup> − I42R<sup>b</sup> = 0 I4(2R<sup>a</sup> + 4Rb) − I1R<sup>a</sup> − I9R<sup>a</sup> − I32R<sup>b</sup> − I52R<sup>b</sup> = 0 I5(2R<sup>a</sup> + 4Rb) − I1R<sup>a</sup> − I10R<sup>a</sup> − I42R<sup>b</sup> − I62R<sup>b</sup> = 0 I6(2R<sup>a</sup> + 3Rb) − I1R<sup>a</sup> − I11R<sup>a</sup> − I52R<sup>b</sup> = 0 I7(2R<sup>a</sup> + 3Rb) − I2R<sup>a</sup> − I82R<sup>b</sup> = 0 I8(2R<sup>a</sup> + 4Rb) − I3R<sup>a</sup> − I72R<sup>b</sup> − I92R<sup>b</sup> = 0 I9(2R<sup>a</sup> + 4Rb) − I4R<sup>a</sup> − I82R<sup>b</sup> − I102R<sup>b</sup> = 0 I10(2R<sup>a</sup> + 4Rb) − I5R<sup>a</sup> − I92R<sup>b</sup> − I112R<sup>b</sup> = 0 I11(2R<sup>a</sup> + 3Rb) − I6R<sup>a</sup> − I102R<sup>b</sup> = 0 (24)

Through the resistance formulas of 2 courses × 3 wales, 2 courses × 4 wales, and 2 courses × 5 wales (Equation (25)), the resistance expression law of 2 courses × N wale cannot be directly obtained, so the equivalent resistance of the knitted sensor topological model can instead be obtained through the corresponding equations of M course × N wale, as in Equation (26), and as shown in Figure 7 (the equivalent circuit of the knitted strain sensor with M course × N wale).

 R(2,3) = Ra(9R<sup>a</sup> <sup>4</sup>+94R<sup>a</sup> <sup>3</sup>Rb+271R<sup>a</sup> <sup>2</sup>R<sup>b</sup> <sup>2</sup>+260RaR<sup>b</sup> <sup>3</sup>+48R<sup>b</sup> 4 ) 9Ra <sup>4</sup>+16R<sup>b</sup> <sup>4</sup>+84R<sup>a</sup> <sup>3</sup>Rb+187R<sup>a</sup> <sup>2</sup>R<sup>b</sup> <sup>2</sup>+112RaR<sup>b</sup> 3 R(2,4) = 2Ra(6R<sup>a</sup> <sup>4</sup>+45R<sup>a</sup> <sup>3</sup>Rb+92R<sup>a</sup> <sup>2</sup>R<sup>b</sup> <sup>2</sup>+62RaR<sup>b</sup> <sup>3</sup>+8R<sup>b</sup> 4 ) 9Ra <sup>4</sup>+4R<sup>b</sup> <sup>4</sup>+60R<sup>a</sup> <sup>3</sup>Rb+95R<sup>a</sup> <sup>2</sup>R<sup>b</sup> <sup>2</sup>+40RaR<sup>b</sup> 3 R(2,5) = Ra(45R<sup>a</sup> <sup>6</sup>+690R<sup>a</sup> <sup>5</sup>Rb+3627R<sup>a</sup> <sup>4</sup>R<sup>b</sup> <sup>2</sup>+8348R<sup>a</sup> <sup>3</sup>R<sup>b</sup> <sup>3</sup>+8592R<sup>a</sup> <sup>2</sup>R<sup>b</sup> <sup>4</sup>+3520RaR<sup>b</sup> <sup>5</sup>+320R<sup>b</sup> 6 ) 27Ra <sup>4</sup>+64R<sup>b</sup> <sup>4</sup>+396R<sup>a</sup> <sup>6</sup>Rb+1929R<sup>a</sup> <sup>4</sup>R<sup>b</sup> <sup>2</sup>+3920R<sup>a</sup> <sup>3</sup>R<sup>b</sup> <sup>3</sup>+3232R<sup>a</sup> <sup>2</sup>R<sup>b</sup> <sup>4</sup>+896RaR<sup>b</sup> 5 (25)

 (N + 1)I1R<sup>a</sup> − I2R<sup>a</sup> − I3R<sup>a</sup> − I4R<sup>a</sup> − · · · − IN+1R<sup>a</sup> = *V* I2(2R<sup>a</sup> + 3Rb) − I1R<sup>a</sup> − I32R<sup>b</sup> − IN+2R<sup>a</sup> = 0 I3(2R<sup>a</sup> + 4Rb) − I1R<sup>a</sup> − I22R<sup>b</sup> − I42R<sup>b</sup> − IN+3R<sup>a</sup> = 0 . . . IN+1(2R<sup>a</sup> + 3Rb) − I1R<sup>a</sup> − I2N+1R<sup>a</sup> − IN2R<sup>b</sup> = 0 IN+2(2R<sup>a</sup> + 3Rb) − I2R<sup>a</sup> − I2N+2R<sup>a</sup> − IN+32R<sup>b</sup> = 0 IN+3(2R<sup>a</sup> + 4Rb) − I3R<sup>a</sup> − I2N+3R<sup>a</sup> − IN+22R<sup>b</sup> − IN+42R<sup>b</sup> = 0 . . . I2N+1(2R<sup>a</sup> + 3Rb) − IN+1R<sup>a</sup> − I3N+1R<sup>a</sup> − I2N2R<sup>b</sup> = 0 . . . <sup>I</sup>(M−2)N+<sup>2</sup> (2R<sup>a</sup> + 3Rb) − <sup>I</sup>(M−1)N+2R<sup>a</sup> − <sup>I</sup>(M−3)N+2R<sup>a</sup> − <sup>I</sup>(M−2)N+32R<sup>b</sup> = <sup>0</sup> <sup>I</sup>(M−2)N+<sup>3</sup> (2R<sup>a</sup> + 4Rb) − <sup>I</sup>(M−1)N+3R<sup>a</sup> − <sup>I</sup>(M−3)N+3R<sup>a</sup> − <sup>I</sup>(M−2)N+22R<sup>b</sup> − <sup>I</sup>(M−2)N+42R<sup>b</sup> = <sup>0</sup> . . . <sup>I</sup>(M−1)N+<sup>1</sup> (2R<sup>a</sup> + 3Rb) − <sup>I</sup>MN+1R<sup>a</sup> − <sup>I</sup>(M−2)N+1R<sup>a</sup> − <sup>I</sup>(M−1)N2R<sup>b</sup> = <sup>0</sup> <sup>I</sup>(M−1)N+<sup>2</sup> (2R<sup>a</sup> + 3Rb) − <sup>I</sup>(M−2)N+2R<sup>a</sup> − <sup>I</sup>(M−1)N+32R<sup>b</sup> = <sup>0</sup> <sup>I</sup>(M−1)N+<sup>3</sup> (2R<sup>a</sup> + 4Rb) − <sup>I</sup>(M−2)N+3R<sup>a</sup> − <sup>I</sup>(M−1)N+22R<sup>b</sup> − <sup>I</sup>(M−1)N+42R<sup>b</sup> = <sup>0</sup> . . . (26)

<sup>I</sup>MN+1(2R<sup>a</sup> + 3Rb) − <sup>I</sup>(M−1)N+1R<sup>a</sup> − <sup>I</sup>MN2R<sup>b</sup> = <sup>0</sup>

**Figure 6.** The equivalent resistances of the circuits with 2 courses × 3 wales, 2 courses × 4 wales, and 2 courses × 5 wales.

**Figure 7.** The equivalent circuit of the knitted strain sensor with M course × N wale.

According to the M course × N wale topology model equations on loop transfer and length resistance (take 2 courses × 2 wales as an example to calculate the equivalent resistance ratio during longitudinal strain), R<sup>a</sup> 0 , R<sup>b</sup> 0 , and the equivalent resistance ratio R 0 (2,2)/R(2,2) of 2 courses × 2 wales can be calculated, as shown in Equation (27).

 R(2,2) = 2Ra(R<sup>a</sup> <sup>2</sup>+3RaRb+R<sup>b</sup> 2 ) 3Ra <sup>2</sup>+R<sup>b</sup> <sup>2</sup>+4RaR<sup>b</sup> R 0 <sup>a</sup> = R<sup>a</sup> <sup>1</sup> <sup>−</sup> <sup>ε</sup> 2 R 0 <sup>b</sup> = Rb(1 + ε) R 0 (2,2) = 2Ra(1<sup>−</sup> <sup>ε</sup> 2 ) [Ra(1<sup>−</sup> <sup>ε</sup> <sup>2</sup> )]<sup>2</sup> <sup>+</sup>3Ra(1<sup>−</sup> <sup>ε</sup> <sup>2</sup> )Rb(1+ε)+[Rb(1+ε)]<sup>2</sup> <sup>3</sup>[Ra(1<sup>−</sup> <sup>ε</sup> <sup>2</sup> )]<sup>2</sup> +[Rb(1+ε)]2+4Ra(1<sup>−</sup> <sup>ε</sup> <sup>2</sup> )Rb(1+ε) R 0 (2,2) R(2,2) = (3R<sup>a</sup> <sup>2</sup>+4RaRb+R<sup>b</sup> 2 ) h 2Ra 3 ( ε <sup>2</sup> <sup>−</sup>1) 3 −R<sup>b</sup> 2 (1+ε) <sup>2</sup>+3RaRb( ε <sup>2</sup> <sup>−</sup>1)(1+ε) i 2Ra h Rb <sup>2</sup>(1+ε) <sup>2</sup>+3R<sup>a</sup> 2( ε <sup>2</sup> <sup>−</sup>1) 2 <sup>−</sup>4RaRb( ε <sup>2</sup> <sup>−</sup>1)(1+ε) i (Ra <sup>2</sup>+3RaRb+R<sup>b</sup> 2) (27)

### *3.3. Calculation of the Equivalent Resistance Model Based on the Three-Dimensional Curved Surface Strain*

According to the change of the sensor surface area, the relationship between the burst depth (D) and curved surface strain is established. The volume resistivity formula is used to calculate the resistance (Rv) before breaking and the equivalent resistance (R 0 v ) after the curved surface strain. The ratio of surface strain resistance is derived as shown in Equation (28). To avoid the burst depth of the formula being zero in an initial state, the burst depth (D) in the initial state is set to 1 mm, and bottom surface cone radius (r) is set to one-half of a stainless steel radius. Therefore, the resistance change ratio of the curved surface strain is defined as y<sup>v</sup> , and the relationship between y<sup>v</sup> and the burst depth (D) is expressed as Equation (29).

$$\begin{cases} \begin{aligned} \mathbf{R}\_{\mathbf{v}} &= \frac{\rho \mathbf{D}}{\mathbf{S}\_{\mathbf{v}}}\\ \mathbf{R}\_{\mathbf{v}}' &= \frac{\rho \mathbf{D}'}{\mathbf{S}\_{\mathbf{v}}'} = \frac{\rho \mathbf{D}'}{\pi \mathbf{I}'(\mathbf{R}+\mathbf{r}) + \pi \mathbf{r}^2} \\\ \mathbf{I}' &= \sqrt{(\mathbf{R}-\mathbf{r})^2 + \mathbf{D}'} \end{aligned} \tag{28}$$
 
$$\begin{aligned} \frac{\mathbf{R}\_{\mathbf{v}}'}{\mathbf{R}\_{\mathbf{v}}} &= \frac{\mathbf{S}\_{\mathbf{v}} \mathbf{D}'}{\pi \mathbf{I}' \mathbf{D}(\mathbf{R}+\mathbf{r}) + \pi \mathbf{r}^2 \mathbf{D}} = \frac{7.5 \mathbf{D}'}{0.88705 \sqrt{\mathbf{D}'^2 - 2.5 \mathbf{D}' + 4.84 + 0.12265625}} \\\ \mathbf{y}\_{\mathbf{v}} &= \frac{7.5 \mathbf{D}}{0.88705 \sqrt{\mathbf{D}^2 - 2.5 \mathbf{D} + 4.84} + 0.12265625} \end{aligned} \tag{29}$$

### *3.4. Fabric Cycle Test Results under In-Plane Stretching*

The initial strain responsiveness was tested by stretching the samples to produce a 100% strain at a speed of 300 mm/min. Following 2800 occurrences of in-plane stretching, 50 groups of cycle data were extracted according to the cyclic stretching data for the analysis of variance, as shown in Figure 8. The results show that the *p*-value of the in-plane stretching is 0.955 > 0.05, as shown in Table 1, indicating that the extracted actual resistance data have no significant positive influence on all cycle data. It proved that an equivalent resistance model established for the same knitted sensor is feasible.



*3.5. Comparison of the Experimental Data and the Equivalent Resistance Model Calculation Results during In-Plane Stretching*

The equivalent resistance model during in-plane stretching is composed of macro– micro models; that is, the combination of the Equations (7) and (11). Therefore, the in-plane tensile equivalent resistance model Equation (30) can be derived. Figure 9 illustrates the simulation results of the equivalent resistance and the experimental data of the inplane stretching, where the two curves Y and y are close to each other (within acceptable accuracy), and further analyzes the fit of the two curves. It can be seen from Table 2 that the *R* <sup>2</sup> value of the model is 0.975, meaning that the in-plane strain fitting resistance can explain 97.5% of the change of real resistance. When the *F* test was performed on the model, it was found that the *F* test (*F* = 39060.82, *p* = 0.00 < 0.05) was applied to the model, which indicated that the in-plane fitting resistance curve will have a significant positive impact on a real resistance curve. The equivalent resistance theoretical model of in-plane tensile properties can be used to predict the same knitted fabric resistance according to the given strain, and the same method can be used in future studies to predict the deformation of the fabric based on the actual resistance obtained.

$$\mathbf{y} = \mathbf{y}\_1 + \mathbf{y}\_2 = \varepsilon^3 + 3\varepsilon^2 + 3\varepsilon + \frac{2\mathbf{n}(\varepsilon + 1)^3}{\mathbf{m} + 1} + 1 \tag{30}$$

**Figure 9.** Comparison of the actual resistance and the in-plane tensile equivalent resistance model calculation results.

**Table 2.** Parameter estimates for in-plane tensile equivalent resistance model.


### *3.6. Comparison of the In-Plane Tensile Test Results and the Equivalent Resistance Calculation Results of the Topology Model*

Comparing the simulation results of the knitting sensor topology resistance model with the results of longitudinal stretching shows that the two curves are close to each other within acceptable accuracy, further verifying the fit of the two curves, as shown in Figure 10. It can be seen from Table 3 that the *R* <sup>2</sup> value of the model is 0.99, meaning that the topology model fitting resistance can explain 99% of the change of real resistance. When the *F* test was performed on the model, it was found that the *F* test (*F* = 43,544.39, *p* = 0.00<0.05) was applied to the model, which indicated that the topology model fitting the resistance curve will have a significant positive impact on the real resistance curve. The equivalent resistance theoretical model of the topology model can be used to predict the same knitted fabric resistance according to a given strain, and the same method can be used in future studies to predict the deformation of the fabric based on the actual resistance obtained.

Comparing the two equivalent resistance models, the macro–micro model establishes the equivalent resistance of the knitted sensor from two perspectives. When the selected strain sensor is larger in size, the macro–micro model can be used to calculate the equivalent resistance more easily and quickly, and the calculation efficiency is higher. The topology equivalent resistance model comprehensively considers the loop transfer and length resistance. When selecting a smaller sensor application, it is more accurate to use the topology resistance model to calculate the equivalent resistance. Compared with the macro–micro models, the topological model has a higher resistance prediction accuracy, but the calculation is more complicated.

**Figure 10.** Comparison of actual resistance and topology equivalent resistance model calculation results.

**Table 3.** Parameter Estimates for topology equivalent resistance model.


### *3.7. Comparison of Three-Dimensional Strain Equivalent Resistance Model and Test Results*

Equation (29) demonstrates a proportional relationship between burst depth and surface strain. Therefore, the relationship between the equivalent resistance ratio and the curved surface strain can be established, as shown in Figure 9. The red curve in Figure 9 aligns with the smooth curve of the burst test data, and the black curve is the predicted result of the equivalent resistance model.

From the analysis in Figure 9, two curves are close to each other within acceptable accuracy, which further verifies the fit of two curves, as shown in Table 4. The fitted resistance is used as an independent variable and actual resistance is used as the dependent variable for linear regression analysis. The accuracy of the fit *R* 2 is 0.867, which means that the fitted resistance can explain 86.7% of the actual resistance. The model passes the *F* test (*F* = 366274.180, *p* = 0.000 < 0.05), which means that the fitted resistance will affect the actual resistance. The specific analysis shows that the regression coefficient value of the fitting resistance is 0.888 (*t* = 605.206, *p* = 0.000 < 0.01), which means that the fitting resistance will have a significant positive influence on the actual resistance. Therefore, a three-dimensional curved surface equivalent resistance model can be used to predict the curved surface resistance change.

**Table 4.** Linear regression analysis of the three-dimensional equivalent resistance model and the actual resistance.


From the analysis in Figure 11, there is a power function relationship between equivalent resistance and burst depth. As the burst depth increases, the resistance change rate tends to be flat after rapid increases in resistance. This change rule better explains three-dimensional curved surface strain-resistance changes, which are consistent with the results of the curved surface resistance test. Therefore, a theoretical formula of equivalent resistance can explain the three-dimensional curved surface strain resistance change rule.

**Figure 11.** The theoretical model of equivalent resistance and the results of the burst test when the three-dimensional surface was strained.

### **4. Conclusions**

In this study, two in-plane equivalent resistance models of knitted strain sensors and a three-dimensional equivalent resistance were presented. It is verified that these three resistance models can predict the resistance change of knitted strain sensors. Based on the in-plane strain, a macro–micro equivalent resistance model and topological equivalent resistance model are established. The macro–micro equivalent resistance model is a combination of two structural models, based on length resistance from a macro point of view and loop transfer from a micro point of view. The topological equivalent resistance model is based on a single loop and considers the superimposed resistance between the loops.

Comparing the in-plane tensile experimental data with the resistance models, the topological resistance model has a higher prediction and fitting accuracy, but the calculation is more complex, which is suitable for small-area knitted sensor applications. The macro–micro equivalent resistance model is simple and quick to calculate, suitable for large-area knitting sensor applications. Based on volume resistivity, a three-dimensional curved surface equivalent resistance model is established, which has a significant positive correlation effect on the change of the curved surface strain resistance. Furthermore, the results strongly suggest that the resistance change occurrence is directly related to the length resistance, loop transfer, and burst depth of the knitted sensor.

**Author Contributions:** Conceptualization, Y.L.; methodology, M.T.; software, M.Y.; validation, Y.L. and P.M.; formal analysis, M.Y.; investigation, Y.L.; resources, P.M.; data curation, Y.L.; writing original draft preparation, Y.L.; writing—review and editing, Y.L.; visualization, M.Y.; supervision, M.T.; project administration, M.T.; funding acquisition, M.Y. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Opening Fund of the China National Textile and Apparel Council Key Laboratory of Flexible Devices for Intelligent Textile and Apparel, Soochow University, grant number SDHY2106.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Informed consent was obtained from all subjects involved in the study.

**Data Availability Statement:** Not applicable.

**Acknowledgments:** We wish to thank the Engineering Research Center for Knitting Technology of Jiangnan University for their timely help in knitting fabric samples.

**Conflicts of Interest:** The authors declare no conflict of interest.

### **References**

