Piezoresistive Theory and Numerical Calculation for Carbon Nanotube Polymer Composite

A three-dimensional theory has been established for the piezoresistivity of carbon nanotube (CNT) polymer composites. Based on the Mori–Tanaka method in meso-mechanics theory and considering quantum tunneling effect between CNTs, an approach to calculate equivalent electrical conductivity of composites was proposed. On this basis, a piezoresistive theory, which incorporates the effect of composites’ geometric nonlinearity, was developed for CNT polymer composites. The theory is dependent only on some basic physical parameters of the materials. A finite element formula of the theory for the numerical calculation of piezoresistivity was presented from the analysis of both elastic and electric fields. Numerical simulations demonstrated that the results predicted by the theory were in good agreement with those of the experimental tests. Parameter sensitivity analysis revealed that when both the potential barrier height of the matrix and the initial average separation distance between CNTs increased, the piezoresistivity obviously increased. However, with the increase in aspect ratio and CNT conductivity, the piezoresistivity decreased gradually. A practical engineering application of this theory is also provided.


Introduction
Carbon nanotubes (CNTs) have extraordinary mechanical and electrical properties and have been widely used as reinforcements for the fabrication of composite materials in recent years [1][2][3][4][5].Most CNT composites are polymer-based.It is well-known that polymers are poor in terms of electrical conductivity, only about 10 −16 S/m∼10 −12 S/m [6], while the conductivity of carbon nanotubes is 1∼10 6 S/m [7].Therefore, carbon nanotubes can not only improve the mechanical properties of the composite, but also increase its electrical conductivity.Experiments have shown that when the volume fraction of carbon nanotubes exceeds a certain value, the electrical conductivity of the composite increases sharply by more than 10 orders of magnitude, and then increases slowly, showing obvious percolation-like properties [7,8].In the vicinity of the percolation threshold, CNT composites were found to present piezoresistive behaviors, that is, the electrical conductivity of the composites varied regularly with the variation of the strain [9][10][11].It was further shown by experimental tests that this piezoresistive property depends on two conductivity mechanisms: the conductive network formed by CNTs and the quantum tunneling effect between adjacent CNTs [7,[12][13][14][15].Due to the piezoresistive properties, CNT composites have various potential applications, such as in smart structures and in the construction of advanced sensors [16].Therefore, it is of significance to study the piezoresistive property of CNT composites.
A few experiments have been conducted to test the piezoresistivity of CNT composites.Most of them often display two different types of piezoresistive behaviors: (1) piezoresistivity increases linearly when the strain is below a critical value, and (2) beyond the critical strain, piezoresistivity increases exponentially, and the quantum tunneling effect plays a major role [17][18][19].Since the quantum tunneling conductivity between CNTs varies exponentially with the distance between CNTs, a small variation in composite deformation may lead to a significant change in the conductivity of the composites.More recently, some research in the literature has been devoted to establishing both numerical and analytical models to predict the piezoresistivity of CNT composites, which may provide useful tools for the design of highly sensitive resistance-type strain sensors made with various CNT composites.Alamusi [20] proposed a numerical model to predict the piezoresistive behaviors of CNT composites, where both the conductive network formed by CNTs and the quantum tunneling effect between adjacent CNTs were considered.Wang [15] developed a simple analytical model to predict the piezoresistivity of CNT composites using micromechanics.However, both Alamusi and Wang's models were one-dimensional and linear, and therefore unable to predict the nonlinear variation trend of piezoresistivity.In order to simulate the strong nonlinearity and asymmetry of piezoresistivity, Gong [21] introduced CNTs' piezoresistivity (effect of CNT deformation) into his model, believing that CNT deformation is a dominant mechanism for the nonlinearity and asymmetry of piezoresistivity of CNT composites.Hu [22] proposed a comprehensive multi-scale resistor network for the prediction of piezoresistivity numerically, in which three mechanisms were involved: (a) a change of conductive network, (b) the quantum tunneling effect, and (c) CNTs' piezoresistivity.This study shows that the contribution from CNTs' piezoresistivity is quite small and can be neglected.Nevertheless, this is inconsistent with Gong's model.
In fact, experiments [17][18][19] showed that it is only when the critical strain of composites is exceeded that the piezoresistivity will increase non-linearly.Generally, the critical strain is so large that the small deformation hypothesis no longer applies.The significant deformation of composites changes the conductivity network, as well as the quantum tunneling distance between CNTs, leading to the nonlinear variation of piezoresistivity.Therefore, in order to reflect the nonlinear variation of piezoresistivity, it is necessary to consider the effect of geometrical nonlinearity (large deformation) on the piezoresistivity of composites.Unfortunately, there are few theoretical or numerical studies involving geometrical nonlinearity.On the other hand, the models mentioned above only predict one-dimensional piezoresistivity, that is, the piezoresistive effect due to uniaxial tension or compression.In fact, composites are usually subjected to a three-dimensional stress state.Thus, the effect of a complex stress state should be included.
In this study, the piezoresistivity of CNT polymer composites is analytically investigated.Based on the theory of meso-mechanics, a new method to calculate the equivalent conductivity of CNT polymer composites is proposed, which includes the mechanisms of both the quantum tunneling effect and CNTs' conductive network.On this basis, a three-dimensional theory to predict the piezoresistivity of CNT polymer composites under a complex stress state is established, and a numerical approach using the finite element method is developed.For the theory, the effects of geometrical nonlinearity on CNTs conductivity network and quantum tunneling electrical conductivity are considered.With this theory, both the electrical conductivity and piezoresistivity of CNT composites can be calculated by the use of some basic physical parameters of CNTs and polymers.
Notation: Compact tensor and matrix notation are used in this paper.As general rules, scalars are denoted by italic light-face Greek or Latin letters (e.g., π or f ); for vectors, second-and fourth-order tensors are signified by italic boldface characters (e.g., σ or ε); and the matrix is denoted by characters within [ ] (e.g., [σ]).

Piezoresistivity of CNT Composites Based on Deformation
It is assumed that CNT composites are elastic.The governing equations of displacement u, strain ε, and stress σ are where f is the body force and D is the elasticity tensors.Note that in Equation (1b), the geometrical nonlinearity term (∇ ⊗ u) • (∇ ⊗ u) T is involved.Boundary conditions: where T and u are boundary tractions and displacements, respectively; n is the outer unit normal vector of boundary.
The potential energy functional of an elastic body is: Equation ( 4) can be rewritten as follows if the elastic body is discretized into a mesh of finite elements: The superscript 'e' represents the element, and the superscript 'T' denotes transpose.
To facilitate the finite element analysis, matrix formulation was adopted as follows.In this paper, the iso-parametric element of a 20-node hexahedron was used for the discretization.For a particular element, displacements can be represented as: where [u e ] is element displacement vector, [N] is the shape function matrix, and [a e ] is the column vector of nodal displacements for a particular element.Substitution of Equation ( 6) into (1b) yields where where λ ix j T = ∂u i ∂x j ∂v i ∂x j ∂w i ∂x j (j = 1, 2, 3.For the hexahedral element of 20 nodes, i = 1, 2, 3, . . ., 20). where , and N i is the element shape function [23].Substituting Equations ( 6) and (7) into Equation ( 5) yields: where [G e ] is the transformation matrix; )dv is the element stiffness matrix; and T T dS represents the equivalent nodal load of the element. Because Equation ( 12) is a nonlinear equation, where [K m ] is a function of nodal displacement a.By introducing the displacement boundary conditions and solving the system of nonlinear equations by the iterative method, a nodal displacement column vector a can be obtained.Then, the volumetric strain at each point was obtained as follows: As the elastic modulus of CNTs is much larger than that of polymers, the deformation of the composite mainly came from the polymer matrix.Therefore, the volume change of CNTs can be neglected.The volume fraction of CNTs can be represented as: where V 0 is the initial volume of the composite and f 0 is the initial volume fraction of the CNTs.
The average separation distance between adjacent CNTs was taken to be d a .Several works have shown that d a is a function of the CNT volume fraction [24].
where d cp and f cp are the average separation distance and CNT volume fraction corresponding to the percolation threshold, respectively.It is shown that d cp = 1.8nm [13].f cp can be approximated by the aspect ratio of the CNTs [25]: where: is the effective length of the CNT, and D is its diameter.
Materials 2023, 16, 7090 5 of 14 Physical studies have proven that there is a penetration effect between CNTs.Simmons derived the conductivity resulting from the perforation effect as [26]: where γ is the potential barrier height, representing the energy required for the electron transition through the matrix; m (m = 9.10938291 × 10 −31 kg) is the electron mass; e (e = 1.602176565 × 10 −19 C) is the electric charge on an electron; and h (h = 6.626068 × 10 −34 m 2 kg/s) is the Planck constant.The polymer matrix is usually an insulator, so the quantum tunneling conductivity between CNTs can be approximated to that of the matrix.Using the Mori-Tanaka method [27], the effective electrical conductivity of CNT composites σ e can be obtained as follows: where with where S 33 = 1 − 2S 11 , n = σ cnt /σ m and σ cnt is the conductivity of the CNTs.
Since the volume resistivity of CNT composites is ρ = 1 /σ e , the piezoresistivity of the CNT composites can be obtained by where ρ 0 is the initial resistivity and σ e0 is the initial conductivity.
Note that σ e is a function of composite volumetric strain, and thus, the piezoresistivity will vary with the deformation of the composites.However, for a structure, volume strain is not uniformly distributed; for example, in a beam, the lower part is characterized by expansion and the upper part by compression, so the application of Equation ( 22) in a practical structure is not convenient.Therefore, an alternative method for the calculation of piezoresistivity was developed as follows.

Piezoresistivity of CNT Composites Based on an Electric Field
Let the CNT composite be subjected to an electric field, as shown in Figure 1.The electric potential ϕ, current density J tensor, and electric field strength E tensor of each point inside the composite satisfy the following relationships: Boundary conditions: where q and ϕ are the surface current density and electric potentials on the boundary surface, respectively.σ e is the conductivity tensor, and for isotropic case, σ e = σ e δ, where δ is the unit tensor.The equivalent integral form of Equation ( 25) and boundary conditions ( 27) is where -node hexahedron iso-parametric element was used for finite element an The computational domain was divided into a number of hexahedral elements, an the corner points and midpoints of each side were taken as nodes.For a particular el ( e ), the electric potential could be determined by: x y z , ,...., , ....., where i N is the shape function and e i  is the nodal electric potential.
  CNT composite under an electric field.
The equivalent integral form of Equation ( 25) and boundary conditions ( 27) is: where v and v are the weight function.The variation in electric potential was selected as the weight function, i.e., v(x, y, z) = δϕ(x, y, z), and v = −v was adopted.Considering ( 23) and ( 24), Equation (28) was rewritten as where

−ϕqdS (30)
A 20-node hexahedron iso-parametric element was used for finite element analysis.The computational domain was divided into a number of hexahedral elements, and both the corner points and midpoints of each side were taken as nodes.For a particular element (e), the electric potential could be determined by: where where N i is the shape function and ϕ e i is the nodal electric potential. where After discretization and substituting Equations ( 33) and (31) into Equation (30), we obtained: where The nodal electric potential vector of each element [ϕ e ] can be written as: where [G e ] is the transformation matrix, [ϕ] denotes the electric potential of all nodes, and ϕ k (k = 1, 2, 3 . . ., 20) denotes a particular nodal electric potential.Then, Equation (35) becomes: where Substitution of Equation (37) into Equation (29) yields Equation ( 38) is a set of linear equations.In order to solve Equation (38), boundary conditions ( 26) must be considered.
The electric potential of all nodes can be obtained by solving Equation (38).The voltage between point a and b is: The volume resistivity, voltage U, and current I satisfy the following relationships: where C is a constant and depends only on the geometry.The piezoresistivity of CNT composites can be determined by a coupling analysis of both the mechanical and electric fields.As shown in Figure 1, for each step of mechanical loading, the finite element analysis as stated in Section 2 was carried out to obtain the elastic deformation.Then, both the volume strain θ and electrical conductivity σ e of each point were calculated according to Equations ( 13)-(19).Thereafter, the voltage between two points before and after loading with the same current I was obtained with the electric field.Therefore, from Equation (40), the piezoresistivity of CNT composites can be represented as: where U 0 ab is the initial voltage before loading between a and b, and U ab is the voltage after loading.Using the above coupling analysis, piezoresistivity can be determined for any position in a non-uniformly deformed structure.
It should be noted that Equation ( 41) is very important for practical applications.With this equation, we can simply measure the piezoresistivity of CNT composites by measuring the voltage before and after loading.

Method Validation
The piezoresistivity of multi-walled CNT composites under tension was studied in reference [19].The matrix material was PEO (polyethylene oxide), and the initial volume fractions of multi-walled CNTs were 0.56% and 1.44%, respectively.The basic physical parameters of the composite were as follows: a Poisson's ratio of 0.46 [19] and an elastic modulus of 7.86 MPA [19].The CNTs' average diameter D was 15 nm [11]; the average length L c was 5 µm [11], and the conductivity of the carbon nanotubes was 100 S/m [11] d cp = 1.8 nm.The barrier height of the polymers was usually 0.5-2.5 eV [7].Therefore, 0.5, 1.5, and 2.5 eV were adopted.The experimental specimen (a) and the finite element mesh (b) for the piezoresistivity of CNT composites are shown in Figure 2. The dimension of the specimens was 160 mm × 40 mm × 40 mm.The CNTs were initially uniformly distributed in a resin matrix.The external electric field was applied at point A and B, and the ends of the specimen were subjected to tensile loads.The iso-parametric element of hexahedron 20 nodes was used in this simulation.The domain was divided into 16 × 4 × 4 = 256 elements, in which there were a total of 1505 nodes.The piezoresistivity obtained according to Equation (41) was compared with the experimental results shown in Figures 3 and 4, respectively.
As shown in Figure 3, the piezoresistivity of f 0 = 0.56% was compared with the experimental results.When the strain increased from 0 to 0.6%, both the test results and the simulation results increased linearly, and geometric nonlinearity was negligible.We noted that when the barrier height was between 1.5 and 2.5 eV, the simulation results agreed well with the test results.When the strain was between 0.6% and 0.8%, the simulation results presented a similar non-linear growth trend as the test results, which was caused by the geometric nonlinearity of the specimens.However, when the strain was greater than 0.008, the piezoresistivity increased sharply and tended towards infinity, indicating that damages or fractures in the matrix probably occurred between CNTs.Although the simulation results also show a very large non-linear trend, there was still a significant difference between the simulation and the tests, since the present approach did not incorporate the effect of damages or fractures in the composites.  .The barrier height of the polymers was usually 0.5-2.5 eV [7].Therefore, 0.5, 1.5, and 2.5 eV were adopted.The experimental specimen (a) and the finite element mesh (b) for the piezoresistivity of CNT composites are shown in Figure 2. The dimension of the specimens was m 160mm 40mm 40m .The CNTs were initially uniformly distrib- uted in a resin matrix.The external electric field was applied at point A and B, and the ends of the specimen were subjected to tensile loads.The iso-parametric element of hexahedron 20 nodes was used in this simulation.The domain was divided into16 4 4 256 elements, in which there were a total of 1505 nodes.The piezoresistivity obtained according to Equation (41) was compared with the experimental results shown in Figure 3 and Figure 4, respectively.As shown in Figure 3, the piezoresistivity of 0 f = 0.56% was compared with the experimental results.When the strain increased from 0 to 0.6%, both the test results and the simulation results increased linearly, and geometric nonlinearity was negligible.We noted that when the barrier height was between 1.5 and 2.5 eV, the simulation results agreed well with the test results.When the strain was between 0.6% and 0.8%, the simulation results presented a similar non-linear growth trend as the test results, which was caused by the geometric nonlinearity of the specimens.However, when the strain was greater than 0.008, the piezoresistivity increased sharply and tended towards infinity, indicating that damages or fractures in the matrix probably occurred between CNTs.Although the simulation results also show a very large non-linear trend, there was still a significant difference between the simulation and the tests, since the present approach did not incorporate the effect of damages or fractures in the composites.
Figure 4 shows the results of piezoresistivity for f0 = 1.44%.Under the same experimental conditions, highly discrete results were obtained from two specimens, which may have been due to the nonuniform distribution of the CNTs.When the strain increased from 0 to 1%, both the test results and the simulation results increased linearly, and the geometric nonlinearity had little influence.When the strain was greater than 1%, geometric nonlinearity played an important role, leading to the non-linear growth of piezoresistivity.It can be observed that the simulation results reflected this growth trend well.With the very large strain of 5%, the simulation results were in good agreement with the experimental results.
The above analyses show that, although further studies are still needed, the present approach can provide a reasonable prediction of the piezoresistivity of CNT composites.

Analysis of Parameters Affecting Piezoresistivity
Theoretically, the parameters which may affect the conductivity of CNT composites include both the aspect ratio and the conductivity of CNTs, the potential barrier height of the matrix, and the initial average separation distance between adjacent CNTs.In the following section, their effects on the prediction of piezoresistivity are analyzed according to Equation (41).

Influence of Aspect Ratio 
The effects of the aspect ratio ( 200 333.33 , and 500 ) on the prediction of pie- zoresistivity are illustrated in Figure 5.It is shown that the influence of  was weak when the strain was low, and the effect of  became more remarkable with the increase in strain.With the increase in the aspect ratio, the piezoresistivity of the composites Figure 4 shows the results of piezoresistivity for f 0 = 1.44%.Under the same experimental conditions, highly discrete results were obtained from two specimens, which may have been due to the nonuniform distribution of the CNTs.When the strain increased from 0 to 1%, both the test results and the simulation results increased linearly, and the geometric nonlinearity had little influence.When the strain was greater than 1%, geometric nonlinearity played an important role, leading to the non-linear growth of piezoresistivity.It can be observed that the simulation results reflected this growth trend well.With the very large strain of 5%, the simulation results were in good agreement with the experimental results.
The above analyses show that, although further studies are still needed, the present approach can provide a reasonable prediction of the piezoresistivity of CNT composites.

Analysis of Parameters Affecting Piezoresistivity
Theoretically, the parameters which may affect the conductivity of CNT composites include both the aspect ratio and the conductivity of CNTs, the potential barrier height of the matrix, and the initial average separation distance between adjacent CNTs.In the following section, their effects on the prediction of piezoresistivity are analyzed according to Equation (41).

Influence of Aspect Ratio α
The effects of the aspect ratio (α = 200, 333.33, and 500) on the prediction of piezoresistivity are illustrated in Figure 5.It is shown that the influence of α was weak when the strain was low, and the effect of α became more remarkable with the increase in strain.With the increase in the aspect ratio, the piezoresistivity of the composites decreased gradually.This is because the aspect ratio affected the formation of the conductive network and the quantum tunneling effect, and hence affected the conductivity of the composites.Therefore, when the strain was large, the influence of the aspect ratio on the piezoresistivity of the composites was remarkable.
decreased gradually.This is because the aspect ratio affected the formation of the conductive network and the quantum tunneling effect, and hence affected the conductivity of the composites.Therefore, when the strain was large, the influence of the aspect ratio on the piezoresistivity of the composites was remarkable.

Influence of Conductivity of CNTs
As shown in Figure 6, the influence of the conductivity of the CNTs ( , 1000S / m and 10 000 m S/ ， ) on the piezoresistivity was analyzed.The figure shows that the influence of the conductivity of the CNTs on piezoresistivity was not obvious when the strain is small; when the strain increased, the piezoresistivity decreased with the increase in the CNTs' conductivity.Particularly, when 10,000 cnt  = , the decreasing trend became more remarkable.This is because the significant CNTs conductivity led to insensitivity of the composite's conductivity to strain, which weakened the composite's piezoresistivity.

Influence of Conductivity of CNTs
As shown in Figure 6, the influence of the conductivity of the CNTs (σ cnT = 100 S/m, 1000 S/m and 10,000 S/m) on the piezoresistivity was analyzed.The figure shows that the influence of the conductivity of the CNTs on piezoresistivity was not obvious when the strain is small; when the strain increased, the piezoresistivity decreased with the increase in the CNTs' conductivity.Particularly, when σ cnt = 10,000, the decreasing trend became more remarkable.This is because the significant CNTs conductivity led to insensitivity of the composite's conductivity to strain, which weakened the composite's piezoresistivity.
decreased gradually.This is because the aspect ratio affected the formation of the conduc-tive network and the quantum tunneling effect, and hence affected the conductivity of the composites.Therefore, when the strain was large, the influence of the aspect ratio on the piezoresistivity of the composites was remarkable.

Influence of Conductivity of CNTs
As shown in Figure 6, the influence of the conductivity of the CNTs ( , 1000S / m and 10 000 m

S/ ，
) on the piezoresistivity was analyzed.The figure shows that the influence of the conductivity of the CNTs on piezoresistivity was not obvious when the strain is small; when the strain increased, the piezoresistivity decreased with the increase in the CNTs' conductivity.Particularly, when 10,000 cnt  = , the decreasing trend became more remarkable.This is because the significant CNTs conductivity led to insensitivity of the composite's conductivity to strain, which weakened the composite's piezoresistivity.

Influence of Potential Barrier Height
Shown in Figure 7 are the piezoresistivity curves for different barrier heights (γ = 0.5 eV, 1.5 eV, and 2.5 eV).When the potential barrier height increased from 0.5 eV to 1.5 eV, the piezoresistivity increased obviously.However, the increase amplitude decreased when the barrier height increased from 1.5 eV to 2.5 eV.The potential barrier height represents the energy required for the electron to transition through the matrix.High potential barrier height values made electron tunneling more difficult, which significantly affected the quantum tunneling conductivity.Quantum tunneling conductivity was more sensitive to composite deformation than CNT conductivity.Therefore, piezoresistivity was sensitive to changes in the potential barrier height.
to 1.5eV , the piezoresistivity increased obviously.However, the increase amplitude de- creased when the barrier height increased from 1.5eV to 2.5eV .The potential barrier height represents the energy required for the electron to transition through the matrix.High potential barrier height values made electron tunneling more difficult, which significantly affected the quantum tunneling conductivity.Quantum tunneling conductivity was more sensitive to composite deformation than CNT conductivity.Therefore, piezoresistivity was sensitive to changes in the potential barrier height.

Influence of Initial Average Separation Distance d cp
Shown in Figure 8 are the curves of piezoresistivity for different initial average separation distances (d cp = 1.4 nm, 1.6 nm, and 1.8 nm).The influence of the initial average separation distance on the piezoresistivity was obvious.When d cp increased from 1.4 nm to 1.8 nm, the piezoresistivity increased.The increasing trend became more obvious when the strain grew larger.The initial average separation distance was related to the quantum tunneling conductivity, which was sensitive to composite deformation.Large values of d cp made quantum tunneling more difficult, which increased the resistivity of the composites.Therefore, after the deformation of the composites, the piezoresistivity for large d cp values was higher than that for small d cp values.

Application
In order to demonstrate the application of the present approach to engineering structures, a finite element analysis was conducted on the cantilever beam shown in Figure 1.The dimensions of beam were 800 mm × 200 mm × 200 mm, and there was a concentrated force P acting on its cantilever end.The finite element mesh is shown in Figure 9.The beam was composed of CNT composites, with an elastic modulus of

Application
In order to demonstrate the application of the present approach to engineering structures, a finite element analysis was conducted on the cantilever beam shown in Figure 1.The dimensions of beam were 800 mm × 200 mm × 200 mm, and there was a concentrated force P acting on its cantilever end.The finite element mesh is shown in Figure 9.The beam was composed of CNT composites, with an elastic modulus of E = 3 GPa and a Poisson's ratio of v = 0.3.The initial CNT volume fraction was f 0 = 1.21%; the average diameter of the CNTs = 15 nm; and the average length L c = 5 µm.The CNT conductivity was 10,000 S/m, γ = 0.5 eV, and d cp = 1.8 nm.
was 10 000 S/ m ， , A steady-state electric field was applied before loading, with the current flowing in at point A at the bottom and flowing out at point B .Two sets of measuring points, 2 N , were used to test the voltage changes in the tension and compres- sion zones, respectively, in order to determine the corresponding piezoresistive effect.Figure 10 shows the variation curve of relative voltage in the tensile zone with the variation in the strain.Under concentrated force, the strain increased with the loads, and the piezoresistive behavior of CNT composites increased significantly.When the strain was low, the relative change in voltage increased linearly.With the further increase in strain, the piezoresistive behavior gradually became non-linear.The reason for this is that when the tensile strain was high, it led to an increase in the spacing between CNTs.The tunneling conductivity changed exponentially due to the influence of geometric nonlinearity, resulting in a non-linear increase in the relative voltage variation at the macro level.A steady-state electric field was applied before loading, with the current flowing in at point A at the bottom and flowing out at point B. Two sets of measuring points, M 1 N 1 and M 2 N 2 , were used to test the voltage changes in the tension and compression zones, respectively, in order to determine the corresponding piezoresistive effect.
Figure 10 shows the variation curve of relative voltage in the tensile zone with the variation in the strain.Under concentrated force, the strain increased with the loads, and the piezoresistive behavior of CNT composites increased significantly.When the strain was low, the relative change in voltage increased linearly.With the further increase in strain, the piezoresistive behavior gradually became non-linear.The reason for this is that when the tensile strain was high, it led to an increase in the spacing between CNTs.The tunneling conductivity changed exponentially due to the influence of geometric nonlinearity, resulting in a non-linear increase in the relative voltage variation at the macro level.When the strain was small, the nonlinear term of volume strain (Equation ( 13)) was so small that it could be ignored, and therefore, the piezoresistive behavior was linear.When the strain was small, the nonlinear term of volume strain (Equation ( 13)) was so small that it could be ignored, and therefore, the piezoresistive behavior was linear.Figure 11 shows that the piezoresistive behavior in the compressive zone decreased with increasing compressive strain.When the compressive strain was low, its linear trend was more obvious.As the compressive strain increased, the relative change in voltage gradually became non-linear.Compared with the tensile zone, this trend was smaller, and the variation in piezoresistive behavior was weaker.This is because the tensile strain was greater than the compressive strain under the same load.Therefore, the piezoresistive behavior due to geometric nonlinearity in the compressive zone was weaker.Figure 11 shows that the piezoresistive behavior in the compressive zone decreased with increasing compressive strain.When the compressive strain was low, its linear trend was more obvious.As the compressive strain increased, the relative change in voltage gradually became non-linear.Compared with the tensile zone, this trend was smaller, and the variation in piezoresistive behavior was weaker.This is because the tensile strain was greater than the compressive strain under the same load.Therefore, the piezoresistive behavior due to geometric nonlinearity in the compressive zone was weaker.
Figure 11 shows that the piezoresistive behavior in the compressive zone decreased with increasing compressive strain.When the compressive strain was low, its linear trend was more obvious.As the compressive strain increased, the relative change in voltage gradually became non-linear.Compared with the tensile zone, this trend was smaller, and the variation in piezoresistive behavior was weaker.This is because the tensile strain was greater than the compressive strain under the same load.Therefore, the piezoresistive behavior due to geometric nonlinearity in the compressive zone was weaker.

v
and v are the weight function.The variation in electric potential was se as the weight function, i.e.,

Figure 2 .
Figure 2. Testing specimen (a) and finite element mesh (b) for the piezoresistivity of CNT composites.

Figure 2 .Figure 3 .
Figure 2. Testing specimen (a) and finite element mesh (b) for the piezoresistivity of CNT composites.

, 1 .
cp d Shown in Figure 8 are the curves of piezoresistivity for different initial average sep-6nm , and 1.8nm ).The influence of the initial aver- age separation distance on the piezoresistivity was obvious.When cp d increased from 1.4 nm to 1.8 nm, the piezoresistivity increased.The increasing trend became more obvious when the strain grew larger.The initial average separation distance was related to the quantum tunneling conductivity, which was sensitive to composite deformation.Large values of cp d made quantum tunneling more difficult, which increased the resistivity of the composites.Therefore, after the deformation of the composites, the piezoresistivity for large cp d values was higher than that for small cp d values.

.
The initial CNT volume fraction was 0 1.21% f = ; the average diameter of the CNTs = 15 nm; and the average length state electric field was applied before loading, with the current flowing in at point A at the bottom and flowing out at point B .Two sets of measuring points,

Figure 9 .
Figure 9. Finite element mesh and measurement points.

Figure 9 .
Figure 9. Finite element mesh and measurement points.

Figure 10 .
Figure 10.The relative variation in voltage with tensile strain (simulation results).

Figure 10 .
Figure 10.The relative variation in voltage with tensile strain (simulation results).