Modeling of the viscoelastic material (VEM) properties of composites consisting of a polymer matrix filled with nanoparticles is essential for predicting the behavior of these composites for different volumetric fractions and physical properties of the constituents. Examples of the nanomaterials include nanotubes and nanoparticles, such as those shown in Figure 11.1. Among the important nanomaterials that can significantly impact the damping characteristics of polymer composites are carbon nanotubes, carbon black (CB) particles, and piezoelectric particles (Aldraihem 2011; Aldraihem et al. 2007).
The modeling techniques of the viscoelastic properties of the polymeric composite materials aim primarily at reducing the heterogeneous composite medium to an equivalent homogenous, anisotropic continuum as illustrated in Figure 11.2. The development approaches of the equivalent properties of the homogenous medium from the geometrical and physical properties of the constituents of the microstructure belong to the well‐established fields of “micromechanics” and “homogenization.”
The most common approaches of micromechanics and homogenization can generally be classified as analytical or numerical methods as is briefly outlined in Figure 11.3. Among the adopted analytical methods are the Halpin–Tsai method, the self‐consistent method, the Mori–Tanaka model, and the double‐inclusion method. Also, the widely used numerical methods include the finite element‐based methods, the method of macroscopic degrees of freedom (DOF), generalized method of cells (GMC), high‐fidelity generalized method of cells (HFGMC), and the variational asymptotic method for unit cell homogenization (VAMUCH) (Sejnoha and Zeman 2013; Nemat‐Nasser and Hori 1999).
Excellent reviews of these analytical and numerical methods are presented by Tucker and Liang (1999), Nemat‐Nasser and Hori (1999), and Torquato (2001).
Besides these homogenization methods, bounding techniques for the stiffness of the composite have been developed by Voigt (1887), Reuss (1929), and Hashin and Shtrikman (1962).
In this chapter, particular emphasis is placed on the Mori–Tanaka method (MTM) (1973), which is based on Eshelby's equivalent inclusion technique (1957). In the MTM, it is assumed that the average strain in the individual inclusion is proportional to the average strain in the matrix by the concentration Eshelby's strain tensor that relates also the strain in the inclusion to the applied strain.
It is important to note here that as the particle‐filled composites are multi‐scale in nature, as the scale of the individual constituents is of a much lower order of magnitude than that of the entire composite structure.
Hence, all these previously mentioned “homogenization” methods, by developing the effective homogeneous medium, simplify computational effort as the focus shifts to analysis of the macroscopic scale of the entire composite rather than on the microscopic scale of the individual constituents, which is computationally exhaustive.
The particle‐filled polymer composite considered in this section, shown in Figure 11.4, consists of a polymer matrix with embedded N nanoparticle phases. The constituents of the composite are assumed to be perfectly bonded. Furthermore, the nanoparticle inclusions can be spheroidal in shape and assume any arbitrary orientation within the matrix.
In this section, all the nanoparticles are considered to be completely passive. But, in Section 11.3, some of these particles are imparted active capabilities such as in piezoelectric nanoparticles and others are utilized to make the polymer conductive. With such a combination, the matrix provides electric current paths and resistive loading in order to transmit and dissipate the electrical energy generated by the piezoelectric inclusions. In other words, the nanoparticle‐filled composite acts as flexible shunted piezoelectric networks similar to those discussed in Chapter 9.
In this section, the theory governing the prediction of the constitutive characteristics of nanoparticle‐filled composites is outlined. The composite is assumed to obey the conventional assumptions which are often used in the micromechanics analysis. Also, the composite is assumed to be homogeneous on the macro‐scale.
The effective viscoelastic properties are derived using the approach introduced by Dvorak and Benveniste (1992) along with the correspondence principle. For a macro‐scale homogeneous composite, a representative volume element (RVE) that contains N perfectly bonded phases can be selected. The volume averaged fields of the hybrid composite are given by:
and
where and denote the overall average stress and strain fields, respectively. Also, fr denotes the volume fraction of phase r. Furthermore, and denote the average stress and strain fields in phase r, respectively.
The correspondence constitutive equation of the overall composite is given by
where [c*] is the overall complex stiffness matrix in the principle material (1‐2‐3) coordinate system.
Similarly, the correspondence constitutive equation in each phase can be written as
where denotes the viscoelastic complex stiffness of phase r in the principle material coordinate system.
It has been shown that the overall average stress of the hybrid composite is given by (Dvorak and Benveniste 1992). Similarly, the overall average strain of the hybrid composite is given by . Thus, under uniform far field [ε0] or [σ0] the volume‐averages fields in phase r can be expressed as
and
where and are the correspondence concentration factors of phase r.
Substituting Eq. (11.1) into Eq. (11.5) and Eq. (11.2) into Eq. (11.6), yields:
and
Hence, Eqs. (11.7) and (11.8) yield the following expression:
where [I] denotes the identity matrix.
Combining Eqs. (11.4) and (11.5), gives:
Substituting Eq. (11.10) into Eq. (11.2), yields:
Comparing Eqs. (11.3) and (11.11), yields:
or
Using Eq. (11.9) yields:
or
i.e.
Substituting Eq. (11.13) into Eq. (11.13), gives the following expression for the overall viscoelastic stiffness and compliance in the principle material coordinate system:
In a similar manner, the corresponding expression for the overall viscoelastic compliance [s*] in the principle material coordinate system is given by:
where and denote the complex stiffness and compliance of the matrix, respectively. Note that Eqs. (11.14) and (11.15) satisfy the internal consistency relationships; that is, the effective stiffness [c*] and compliance [s*] are the inverse of each other.
The correspondence concentration factors, and , can be estimated by the various micromechanics tools such as the Eshelby's method, the self‐consistent method, and the Mori–Tanaka approach, and so on. Among all those, the Mori–Tanaka approach is known to be one of the most elegant and accurate methods that can be used to estimate the overall elastic properties of aligned as well as randomly oriented short fiber composite (Tucker III and Liang 1999; Tandon and Weng 1986). Moreover, the Mori–Tanaka approach yields identical results when either a uniform elastic strain [ε0] or a uniform stress [σ0] field is prescribed at the boundary of the composite.
Inside the matrix, the stress [σ1] is given by the following expression:
However, within the inclusion the stress [σr] is given by the following expression:
According to Eshelby's approach, the constrained strain [εC] is uniquely related to the transformation strain [εT] such that:
where [S*] is the “Eshelby tensor,” which is a function of the Poisson's ratio of the matrix and the shape and geometry of the inclusions.
Furthermore, Eshelby developed an equivalence theorem between the problems of homogeneous inclusion and an inhomogeneous inclusion of the same shape. Such an equivalence can best be understood by considering the bodies shown in Figure 11.6a,b for the homogeneous and inhomogeneous inclusions, respectively. These inclusions have stiffnesses and , respectively.
In Figure 11.6a, the homogeneous inclusion is subjected to a transformation strain [εT] whereas the inhomogeneous inclusion of Figure 11.6b has no transformation strain. The two bodies of Figure 11.6 are subject to a uniform applied strain [εA]. Eshelby developed an expression for the transformation strain [εT] that results in the same stress and strain distributions in both problems. Hence, the constitutive Eq. (11.21) for the homogeneous inclusion can be modified to account for the applied strain [εA] to yield:
Accordingly, the corresponding constitutive equation for the inhomogeneous inclusion can be written as follows:
The equivalence of the stresses between the two problems yields:
Substituting Eq. (11.22) into Eq. (11.25) gives the following expression:
Equation (11.26) gives the expression of the transformation strain [εT] which is necessary for the equivalence in terms of the applied strain transformation strain [εA].
Note that, in the far field, the average strains in the matrix and the rth inclusion can be written as:
and
Then, substituting Eq. (11.28) into Eq. (11.25), it reduces to:
or
that is,
Hence,
Combining Eqs. (11.22) and (11.29) yields:
Substituting Eqs. (11.27) and (11.28) into Eq. (11.30), it reduces to:
or
In the Mori–Tanaka assumption, it is assumed that when many identical particles (inclusions) are embedded into the composite, the average particle strain is given by:
that is, each particle experiences a far‐field strain that is equal to the average strain of the matrix. In other words, a direct comparison between Eqs. (11.31) and (11.32) gives the following expression of the strain concentration factor:
This completes the proof of Eq. (11.18).
Now for proving Eq. (11.16), two approaches can be adopted:
First Approach:
This approach is simple and aiming at showing that the two sides of the equation are equal. It begins by multiplying both sides by fs and summing up over the number of phases s = 1,…N inside the composite to yield the following:
that is, the right hand side of Eq. (11.34) is equal to an identity [I].
But, from Eq. (11.9), it is also seen that the left‐hand side of Eq. (11.34) is also equal to an identity I. Hence, Eq. (11.16) is valid.
Second Approach:
From Eq. (11.1), as:
Then, substituting Eq. (11.32) in this equation, yields:
Pre‐multiplying both sides by and using Eq. (11.32) gives the following expression:
From Eqs. (11.32) and (11.5), Eq. (11.36) reduces to:
Note that the Eshelby tensors for spheroidal inclusions, in an isotopic matrix, can be expressed as follows (Nemat‐Nasser and Hori 1999)
The elements of the Eshelby tensor are a function of the inclusion aspect ratio α = a3/a1 and the Poisson's ratio νm of the matrix as tabulated in Table 11.1. where g1 and g2 appearing in the table are defined as follows:
and
Table 11.1 Elements of the Eshelby tensor for common inclusion shapes (Aldraihem 2011). Reproduced with permission of Elsevier.
Shape | Thin disc (α = 0) | Sphere (α = 1) | Fiber (α → ∞) | Oblate spheroid (0 < α < 1, g = g1) and | Prolate spheroid (1 < α < ∞ , g = g2) |
S1111 | 0 | ||||
S1122 | 0 | ||||
S1133 | 0 | ||||
S3311 | 0 | ||||
S3333 | 1 | 0 | |||
S1313 | |||||
S1212 | 0 |
Two processes are needed to find the effective properties of a hybrid composite with any desired inclusion orientation distribution. The first process is the evaluation of the effective properties for the hybrid composite with unidirectional aligned inclusions. This step is carried out in the previous section. The second process involves averaging the properties obtained from the first process in terms of the orientation distribution. Advani and Tucker III (1987) suggested a straightforward scheme for orientation averaging in short fiber composites. Furthermore, their averaging scheme can deliver reliable and accurate predictions for thermoelastic properties (Gusev et al. 2002). According to the approach of Advani and Tucker III (1987), the orientation average of Γ(θ, ϕ) is denoted by 〈Γ〉 and is defined as follows:
where ϕ and θ denote the two Euler angles that are usually used to describe the spatial orientation of an inclusion as shown in Figure 11.7. The angles (θi, ϕi) and (θf, ϕf) denote the lower and upper limits, respectively.
The inclusion distribution region is defined by (−θ0 = θi) ≤ θ ≤ (θf = θ0) and (π/2 − ϕ0 = ϕi) ≤ ϕ ≤ (ϕf = π/2 + ϕ0), where θ0 and ϕ0 are prescribed angles. The symbol, Ψ(θ, ϕ), denotes the probability distribution function. This function must satisfy specific conditions leading to:
Using Eq. (11.41) into Eq. (11.40) yields the final expression of the orientation average:
In this equation, the function Γ(θ, ϕ) can be used to represent the transformed viscoelastic stiffness, compliance, or any other property tensor in the global (1′‐2′‐3′) coordinate system. An expression for the transformed viscoelastic stiffness, , is determined by employing the transformation (Wetherhold 1988):
Similarly, the transformed viscoelastic compliance obtained by
where the transformation matrices [R1] and [R2] are defined in Appendix 11.A, and superscript T denotes the transpose.
Substituting Eqs. (11.43) and (11.44) into Eq. (11.42), one obtains the orientation average of the viscoelastic stiffness:
and the viscoelastic compliance:
Expressions (11.45) and (11.46) do not guarantee satisfying the internal consistency relationships. In other words, the orientation averaged stiffness (11.45) and compliance (11.46) are not necessarily each other’s inverse for all orientation distribution.
For predictions of the elastic properties, Hine et al. (2002, 2004) have shown that the orientation averaged stiffness and compliance equations, which are similar to (11.45) and (11.46), provide significantly different results for the thermoelastic properties of short fiber composites. Furthermore, Hine and coworkers have demonstrated that the stiffness equation provides accurate predictions of the orientation averaged properties when compared with the results from their numerical simulation.
According to this discussion, Eq. (11.45) will be used to explicitly determine the orientation averaged viscoelastic properties of the hybrid composites. Furthermore, the viscoelastic properties will be determined by evaluating the limits of Eq. (11.45) as the angles θ0 and ϕ0 approach certain values. Three general orientations will be considered. In aligned orientations, all the inclusions are uni‐directionally aligned along the 1′‐axis. To achieve this case, the values of θ0 and ϕ0 are set to approach zero. In two‐dimensional random orientations, the inclusions are randomly oriented in the 1′‐2′ plane. Here, the values of θ0 and ϕ0 are set to approach π and zero, respectively. In three‐dimensional random orientations, the inclusions are randomly oriented in all directions. In this orientation, the values of θ0 and ϕ0 are set to approach π and π/2, respectively. For these three orientations, the viscoelastic stiffness tensor can be written as:
where the components of are listed in Table 11.2.
Table 11.2 Components of complex moduli for different orientations (Aldraihem 2011). Reproduced with permission of Elsevier.
Aligned orientation (θ0 → 0, ϕ0 → 0) | 2D random orientation (θ0 → π, ϕ0 → 0) | 3D random orientation (θ0 → π, ϕ0 → π/2) | |
Several properties can be extracted from the three‐dimensional moduli of the hybrid composite (11.47). For example, the complex modulus in the global coordinate system can be expressed as
and the loss factor is accordingly defined as
with modulus related to the compliance via the following expression:
and
where denotes the storage modulus, denotes the loss modulus, and subscript q denotes the contracted index.
Based on the theory presented in Section 11.2, the computation of the viscoelastic properties of particle‐filled polymer composites is carried out according to the flow chart displayed in Figure 11.8.
Table 11.3 lists the material properties of typical nanoparticles which are commonly used in manufacturing nanoparticle polymer composites.
Table 11.3 Material properties of the ingredients of the hybrid composite (Aldraihem 2011). Reproduced with permission of Elsevier.
Property | Particle | |||||
PZT‐5Ha | PCM51b | MWCNTc | CBd | DGEBA‐DDMe | Hercules 3501‐6f | |
127.2 | 129.22 | 999.66 | 9.2885 | 2.4409 | 7.4879 | |
80.2 | 86.4 | 369.74 | 3.9808 | 1.6962 | 4.5894 | |
84.67 | 83.06 | 369.74 | 3.9808 | 1.6962 | 4.5894 | |
117.44 | 116.9 | 999.66 | 9.2885 | 2.4409 | 7.4879 | |
23 | 28.83 | 314.96 | 2.6538 | 0.3723 | 1.4493 | |
23.5 | 21.41 | 314.96 | 2.6538 | 0.3723 | 1.4493 | |
k31 | 0.388 | 0.37 | ||||
k33 | 0.752 | 0.72 | ||||
k15 | 0.675 | 0.72 | ||||
ηm | – | – | 0.047 | 0.03 |
a Electro Ceramic Division, Morgan Matroc Inc., Bedford, OH.
b Noliac, http://www.noliac.com (PCM51 is also called NCE51).
The isotropic stiffness matrix for the HDPE polymer is given by:
where , , , , , , , and .
Note that where E*= complex modulus of HDPE, storage modulus, and η1= loss factor. Values of and η1 are displayed in Figure 11.5 as a function of frequency ω.
Similarly, for the CB nanoparticles, the stiffness matrix of the second phase of the composite is given by:
Note that is a matrix with real elements.
Figure 11.10 displays the values of the complex modulus of the CB‐filled HDPE polymer composite as predicted by the micromechanics approach outlined in Section 11.2 for CB volume fractions varying between 0 and 20%. These predictions are in excellent agreement with the corresponding experimental results shown in Figure 11.11.
It is important to note that plotting the storage modulus against the loss modulus results in a unified master line on which the characteristics of the CB/HDPE composite collapse for all the considered CB volume fractions. The equation of this line is:
The equation is valid for both the experimental results and micromechanics predictions. This leads to a loss factor η ≅1.
Furthermore, the displayed results indicate that increasing the CB concentration makes the composite stiffer, but at the expense of reducing its damping behavior as quantified by its loss factor.
Consider the Voigt and Reuss composites shown in Figure 11.12. Let and denote the stiffness tensors of the HDPE and CB, respectively. These tensors are as described in Example 11.1.
Then, the effective stiffness tensor of the Voigt model, which assumes the uniformity of the strain is given by:
where f1 and f2 are the volume fractions of HDPE and CB such that f1 + f2 = 1.
Similarly, the effective compliance tensor of the Reuss model, which assumes the uniformity of the stress is given by:
where and are the compliances of HDPE and CB.
Hence, the effective stiffness tensor of the Reuss model is given by:
Using Eqs. (11.52) through (11.54) yields the characteristics displayed in Figure 11.13 for frequencies of 0.01, 1, and 100 Hz. It can be seen that the viscoelastic characteristics of the CB/HDPE composite lie inside the Voigt and Reuss bounds.
In this section, the predictions of the complex modulus of CB/polymer composites are determined using the classical filler reinforcement methods, which are either analytical in nature or based on finite element modeling of the RVE. These predictions are compared with those generated by the Eshelby's strain tensor approach outlined in Section 11.2.
The analytical approach presented here is based on the filler reinforcement theory developed by Smallwood (1944) and Guth (1945). Such a theory is limited to nanoparticle fillers of diameters less than 100 nm, such as CB particles that usually have diameters of about 40 nm, as outlined in Figure 11.14.
The interaction between the CB particles and the host polymer is analytically quantified by the Smallwood–Guth model, which for solid spherical particles is given as follows:
where E0 is the modulus of the polymer matrix, Ev is the modulus of the CB/polymer composite, and v is the volume fraction of the filler in the composite. Equation (11.55) applies to both the storage and loss moduli of the polymer.
Details of the derivation of Eq. (11.55) are given in Appendix 11.B.
Figures 11.16 and 11.17 show the E′ − v characteristics and the Cole–Cole plot of E′ versus E″ following the use of the Smallwood–Guth Eq. (11.55). It can be seen that the equation results in collapsing the E′ − v and E′ − E″ characteristics to form nearly single lines especially for CB volume fractions that are less than or equal to 15%.
These figures also display comparisons with the predictions obtained by using Eshelby's strain tensor approach. Close agreement is obvious between the predictions of the two approaches particularly for volume fractions less than 15%.
Figure 11.19a shows a finite element model of the RVE that is developed in ANSYS using SOLID186 elements. Figure 11.19b shows a unit cell with the associated boundary conditions that is used in analysis in order to simplify the computational effort (Brinson and Lin 1998).
Figure 11.20 displays a comparison between the deflections w in the z direction as predicted by Smallwood (1944) reinforcement theory and the ANSYS finite element of the RVE for a CB/polymer composite with a CB volume fraction of 20%. The Smallwood prediction is based on Eq. (11.B.9) given in the Appendix 11.B. The figures indicate that the contours of the deflection v as predicted by Smallwood model are in close agreement with those predicted by the ANSYS model.
Figure 11.21 displays the contours of the strain and stress, along the y direction, as predicted by the ANSYS finite element of an RVE of the CB/polymer composite with a CB volume fraction of 20%.
Figure 11.22 displays the contours of the deflection, along the z direction, as predicted by the ANSYS finite element of a RVE of the CB/polymer composite with a CB volume fraction ranging from 0 to 16.6%.
The effect of the volume fraction of the filler on the mechanical properties of a CB/polymer composite is determined using the approach outlined by Barbero (2013). In that approach, the components of the stiffness matrix cαβ are determined by using the average stress and average strain such that:
where α and β = 1,…,6.
In Eq. (11.56), the average stresses are evaluated by computing the stress field over the entire RVE when a unit strain εβ is applied along the β direction such that:
where V is the volume of the VRE.
The modulus of elasticity ET and Poisson's ratio vT in the transverse directions (y and z) are extracted from (Barbero 2013) as follows:
and
Figure 11.23 displays comparisons between the effect of the filler volume fraction on the reinforcement, in the z direction, as predicted by the ANSYS Finite Element, Smallwood–Guth model, and Eshelby's method. The reinforcement effect is quantified by the dimensionless ratio of the elastic modulus Ez of the CB/polymer composite, in the z direction, relative to the storage modulus of the pristine polymer E0; that is, Ez/E0.
The figure suggests that the Smallwood–Guth model adequately predicts the reinforcement effect over a practical range of filler volume fraction. Such a prediction is also computationally simple and efficient when compared to the finite element method (FEM) and Eshelby's method.
CB/polymer composites have attracted the attention for several decades primarily because of for their important use in manufacturing automotive tires. But, in this section, the emphasis is placed on their novel utilization in vibration damping applications. These applications stem mainly from of the fact that embedding CB inside the polymer matrix renders it to be electrically conducting. The extent of the electrical conduction depends on the volume fraction of the CB filler into the composite.
An important metric that quantifies the conductivity of the polymer is the “Percolation Threshold.” This threshold defines the concentration of the fillers that makes the polymer conductive. Generally speaking, as the concentration of the CB increases, the conductivity of the composite increases or the resistivity decreases as shown in Figure 11.24.
The figure indicates that the rate of increase of the conductivity, or decrease of resistivity, is slow at low concentrations as shown in region A of Figure 11.24. The resistivity drops rapidly with further increase in the CB as it goes into region B where the rate of change increases by more than 10 orders of magnitude. Further increase of the CB content fails to improve the resistivity as seen during region C.
It is important to physically understand the underlying phenomena that are behind such changes in the conductivity of the polymer matrix. At low CB content, the gap between the CB particles, where the electrons are transmitted, is very large and the resistivity of the composite is approximately that of the polymer matrix. As the concentration increases, the percolation threshold is reached where the resistivity starts to decrease abruptly as a function of the CB loading. In this region, region B, the gap between the CB particles is close but not touching. As a result, the electron must overcome the potential barrier and cross the gap between the CB particles.
In region C, where the CB loading is higher than the percolation threshold, the CB particles form chains that act as pure resistive conduits for conducting electricity as indicated in Figure 11.25. Figure 11.25a displays a schematic drawing of the microstructure of a CB/polymer composite whereas Figure 11.25b shows a scanning electron microscope photograph of the microstructure emphasizing the existence of CB chains.
Ding et al. (2013) and Wang et al. (2005) suggested that the equivalent electrical circuits of the CB/polymer composite depend, in general, on the CB volume fraction and, in particular, on the percolation characteristics as shown in Figure 11.26. In the figure, Ra, Rc, Cc, and L denote the CB aggregate resistance, the contact resistance, the contact capacitance, and inductance, respectively.
The full practical potential of the CB/polymer composites is usually utilized when operating in the first part of zone C where the transmission of electricity is governed primarily by the ohmic mechanism (Ra). Further increase of the CB content makes the mixing and manufacturing of a uniform CB composite rather difficult. More importantly, the complexity of the resulting equivalent electrical circuit of the composites complicates their operation.
For practical considerations, the contact capacitance and inductance are usually very small as reported by Wang et al. (2005). Therefore, the emphasis is placed here only on modeling the CB/polymer composites as resistive elements.
The electrical resistance R, under no loading conditions, of CB/Polymer composites can be determined from:
where ρs, A, and lc denote the resistivity, cross sectional area, and thickness of the composite, respectively.
It is important to note that the electrical resistance of the composite is usually called “piezo‐resistance.” Such a resistance defines changes in the resistance of the CB/polymers due to the application of load.
Generally, the piezo‐resistance depends on the properties of the polymer matrix, filler properties and filler concentration, and applied load. An excellent account that describes the interactions between all these parameters and their effect on the piezo‐resistivity of the conducting polymer composite is given by Zhang et al. (2000, 2001). In their work, Zhang et al. developed a physics‐based mathematical model to predict the piezo‐resistivity of polymers impregnated with 11 different fillers.
The total resistance R of one conducting path is given by:
where Rc = resistance between two adjacent filler particles, Ra = resistance across filler particle, L = number of particles forming one path, and S = number of conducting paths.
If the inter‐particle separation is very large, no current flows. However, as the separation becomes adequately small, a tunneling current I will flow due to the application of a voltage V such that (Simmons 1963):
where m = electron mass, e = electron charge, hP = Planck constant, s = separation between two adjacent particles, and ϕ = height of potential barrier between adjacent particles.
Assuming a2 = the cross sectional area of the conducting particle, then the resistance Rc can be determined from:
where
Because the conductivity of the particles is very large compared with that between two adjacent particles, then Ra ≅ 0, reducing Eq. (11.61) to:
Equation (11.65) can be used to predict the resistance of the conducting polymer composite and it is clear that it varies exponentially with the separation distance s between the particles, which is a function of the applied load or strain experienced by the composite. Now, let us assume that the inter‐particle separation changes from s0 to s due to the application of stress, then the fractional resistance change (−ΔR/R0) can be predicted from:
where R0 is the original resistance. Note that s and s0 can be related to the strain ε and the stress σ by the following relationships:
where E = modulus of elasticity of the polymer matrix.
The initial separation distance s0 is estimated from Han and Choi (1998) as follows:
where D = particle diameter and vf = volume fraction of filler.
Accordingly, Eq. (11.66) reduces to:
Equation (11.69) predicts the piezo‐resistance changes of conducting polymer composite as function of the applied stress σ, modulus of elasticity of the polymer matrix E, filler particle diameter D, filler volume fraction v, and the parameter γ. Equation (11.69) can be rewritten as:or
where
For σ/E < 0.01, Eq. (11.71) can be simplified to:
Expanding the exponential in Eq. (11.72), yields:
or
where Rε is the resistance of the CB/polymer composites following the application of a strain ε. Equation (11.73) indicates that the relationship between Rε/R0 and the σ/E (or ε) is a straight line with a slope of as shown in Figure 11.27.
The value of can be determined, for CB particle‐filled polymer, using the following parameters: mass of the electron m = 9.180938E‐31 kg, hP = Planck constant = 6.626E‐34 m2 kg s−1, and ϕ = 0.05 EV.
These parameters yield a slope of the piezo‐resistive characteristics of 54.97.
The electrical resistivity ρs of the CB/polymer composites can be determined by rewriting Eq. (11.65) as follows:
In Eq. (11.74), it is assumed that Sa2 and LD denote the area and the length of the conducting path, respectively.
Figure 11.28 displays the resistivity of the CB/polymer composite as predicted by Eq. (11.74), with CB particles of an average size of 30 nm (Kaiser 1993), in comparison with the prediction of Schwartz et al. (2000). In their prediction approach, Schwartz et al. determined the electrical resistivity of a composite of CB in an insulating matrix in terms of the volume fraction of the CB and aggregate size and distribution. It is assumed that the resistance of a random lattice composite between sites different aggregates varies exponentially with the gap.
In this section, the CB/polymer is utilized as a means for vibration damping. Consider a rod consisting of an assembly of CB/polymer layers sandwiching, in a periodic manner, piezoelectric layers as shown in Figure 11.29.
The interaction between the dynamics of the piezoelectric layer that is shunted by the resistive CB/polymer layer is described using the constitutive equation of the piezo‐layer as given in Chapter 9 as follows:
where T3 and S3p are the stress and strain along the composite rod. Also, E3 and D3 denote the electric field and the electrical displacement, respectively. Furthermore, , , , and .
This form of constitutive relationship is used to generate the potential and kinetic energies of an element of the CB/polymer composite and shunted piezo‐network in terms of the DOF {Δe} = {ui uj uk}T and the electrical DOF Qe (the charge of the eth cell). Note that the longitudinal displacements of the composite uc (x) and piezo‐layer up (x), at any position x, are given by the following shape function
where {Nc} = {(1−x/lc) x/lc 0} and {Np} = {0 (1−x/lp) x/lp} where lc and lp are equal the length of the composite layer and the piezo‐layer, respectively.
The kinetic energy KE of the eth unit cell can now be determined from:
where
Also, ρl and Al denote the density and cross sectional area of the lth layer (l = p for the piezo‐layer and l = c for the CB/polymer layer, respectively).
The potential energy PE of the eth unit cell is given by:
where S3p = strain in the piezo‐layer = u,xp, = stress on piezo‐layer, S3c = strain in the CB/polymer layer = u,xc, and = stress on CB/polymer layer with denoting the complex modulus of the composite layer.
Substituting from Eq. (11.75) into Eq. (11.78) gives:
Note that: S3p = u,xp = {N,xp}{Δe}, S3c = u,xc = {N,xc}{Δe}, and D3 = Qe/Ap = constant along the piezo‐layer, then Eq. (11.79) reduces to
where = stiffness matrix of piezo‐layer and CB/polymer assembly.
The virtual work δW associated with the inductive‐resistive‐capacitive shunting network of the CB/polymer composite and the external load {Fi} is given by:
Now, the equations governing the dynamics of the/piezo‐patch system can be extracted using the Lagrangian dynamics to give
and
where Le, Re, and Ce are the components of the inductive‐resistive‐capacitive shunting network of the CB/Polymer composite. Note that Re = Rε where it is assumed here that the effect of the strain on the resistance of the CB/Polymer composite is negligible, that is, Re = Rε ≈ R0. This assumption is considered only to simplify the analysis and makes the finite element model linear to enable direct solution. But, when the resistance Rε ≠ R0, the analysis becomes iterative in nature and requires checking of the convergence at every time or frequency step.
These equations can be written in the following matrix form:
or
where , , ,
The equation of motion of the entire CB composite/shunted piezo‐network assembly can then be determined by assembling the mass and stiffness matrices of the individual cells to yield:
where [Mo], [Co], and [Ko] are the overall mass, damping, and stiffness matrices of the entire assembly. Also, {X} and {Fo} denote the structural and electrical DOF as well as load vectors of the entire assembly.
Note that bonding the shunted piezo‐layer to the CB composite has resulted in adding damping to the assembly as quantified by the damping matrix [Co]. It has also modified the stiffness and mass matrices.
As the CB/polymer composite can be treated electrically as a resistive element, then Le and Ce are set equal to zero in Eq. (11.81). Further, using the static condensation method to condense the DOF of the charge Qe leads to:
i.e. then:
Let e = {0 1 − 1}, then is given by:
Substituting Eq. (11.88) into Eq. (11.83) and pre‐multiplying by gives:
or
where and .
For the unit cell, shown in Figure 11.29, the nodal deflection vector {Δe} is defined as:
where uk, uj, and ui denote the bottom, internal, and top deflection vectors as shown in Figure 11.29b.
This vector is condensed to support the Bloch wave propagation (Hussein 2009). Hence, the displacements at the boundaries are related as follows:
where k and L denote the wave number and the length of the unit cell, respectively.
Hence, define an independent nodal deflection vector such that:
The deflection vectors {Δe} and are related as follows:
where is a transformation matrix such as:
Substituting Eqs. (11.94) and (11.95) into Eq. (11.90), it reduces to:
where , , , and .
Equation (11.95) is now cast in the following state‐space form (Meirovitch 2010):
where .
Let the state‐space solution be assumed as
where {Y}c = { {x} + i{z}}c and λ = iω.
Equation (11.98) reduces to the following eigenvalue problem:
In a compact form, Eq. (11.99) reduces to
where
Equating the real and imaginary coefficients on both sides of Eq. (11.100) yields:
Equations (11.102) and (11.103) can be rewritten in compact and standard eigenvalue problem form such that:
where [A]c = [M*]c − 1 [ [D*]Tc[M*]c − 1[D*]c ].
Note that all the entries of the matrix [A]c are functions of the dimensionless wave number kL. Therefore, the eigenvalues of the matrix [A]c can be determined for different values of the wave number kL.
Then, the eigenvalues λs are given by
where s = 1 to n.
The dispersion characteristics of the unit cell are constructed by plotting the resonant frequency ωs against the wave number kL. The resulted dispersion curves further define the zones of stop and pass bands of the cell as described in Chapter 10.
Figure 11.30 parts a–c displays the dispersion characteristics of the non‐conducting VEM, conducting VEM with Re = 10 kohm, and VEM with Re = 1 Mohm, respectively.
Comparing Figure 11.30a,b indicates that making the VEM layer conduct with the CB filler, with an effective resistance of Re = 10 kohm, has resulted in enhancing the stop band characteristics, particularly at the high frequency end. Further increase of the effective resistance by adding an additional external resistance to augment the piezo‐resistance of the VEM conducting layer to reach 1 Mohm has significantly improved the stop band characteristics and extended it to 1 MHz as indicated in Figure 11.30c.
The effect of such improvements in the stop band characteristics is demonstrated clearly in the frequency response behavior of the composite as shown in Figure 11.31a,b.
Figure 11.31a indicates that the embedding the CB filler in the pristine VEM layer has introduced two important effects when the effective piezo‐resistance is Re = 10 Kohm. The first effect appears in the form of the reinforcement of the polymer, as predicted by the reinforcement theory of Smallwood (1944), which has resulted in stiffening the VEM. Such a stiffening effect is manifested by the shift of the resonant frequencies to higher frequency ranges that, in turn, produced a significant attenuation of the transmitted vibration.
The second effect manifests itself by the vibration attenuation at high frequencies, particularly at 12.93 KHz, which is attributed to the energy dissipation due to shunting of the piezo‐layer by the piezo‐resistance of the conducting VEM/CB layer.
Increasing further the effective shunting resistance to Re = 1 Mohm, as indicated in Figure 11.31b resulted in complete blockage of the high frequency vibration around 12 KHz. Such a blockage is due to the enhanced stop band characteristics shown in Figure 11.30c.
Table 11.4 The physical and geometrical properties of the composite.
PROPERTY | Geometry | Piezo‐layer | Viscoelastic layer | ||||||||
lc = lp (m) | Ac = Ap (m2) | sE(m2 N−1) | d33(m V−1) | k33 | ρp (kg m−3) | Ev (N m−2) | α | ωn (rad s−1) | ζ | ρv (kg m−3) | |
VALUE | 0.1 | 0.01 | 16 × 10−12 | 296 × 10−12 | 0.37 | 7600 | 1E4 | 5 | 25,000 | 10 | 1100 |
The hybrid composite considered in this section, shown in Figure 11.32, consists of a polymer matrix, electrically conducting particles, and piezoelectric inclusions. All the ingredients of the composite are assumed to be perfectly bonded. Further, the inclusions are assumed to be spheroidal in shapes that are arbitrarily oriented within the matrix. The considered piezoelectric inclusions are assumed to be transversely isotropic with poling directions and electrodes in the three direction. The state of stress and electrode distributions of all the piezoelectric inclusions are assumed identical and uniform. The electrically conducting particles and the matrix are assumed to form resistive electric current paths R, which are adequate to load the piezoelectric inclusions (Aldraihem 2011; Aldraihem et al. 2007) as displayed in Figure 11.32b. Finally, the composite is assumed to be homogeneous on the macro‐scale and obeys the basic assumptions used in micromechanics analysis.
Consider a hybrid composite that is subjected to a far‐field stress in the composite global (1′‐2′‐3′) coordinate system. Then, each of the embedded piezoelectric inclusions generates an electric energy that is dissipated into the resistance of the surrounding conducting polymer matrix. As outlined in Chapter 9, a resistively shunted piezoelectric element behaves like a VEM with the nonzero compliances, in the principle material (1‐2‐3) coordinate system, which are defined as follows:
In Eqs. (11.106a)–(11.106f) and (11.107), denotes the elements of elastic compliance matrix at constant electric field. Also, k31, k33, and k15 denote the piezoelectric coupling factors, is the dielectric permittivity at constant stress, R is the shunt resistance, and ω is the frequency. Furthermore, a3 defines the inclusion half‐length and A is the electrodes area with normal along the three‐axis as shown in Figure 11.7.
For the considered shunting, Eqs. (11.106a)–(11.106f) indicates that the out‐of‐plane shear compliances exhibit open‐circuit conditions, the in‐plane shear compliance exhibits short‐circuit conditions, whereas the remaining compliances exhibit shunted‐circuit conditions. Hence, the shear compliances primarily behave in an elastic manner and have limited damping capacity.
Substituting Eq. (11.107) into Eqs. (11.106a)–(11.106f), yields the complex compliance matrix that can be separated into real and imaginary parts as follows:
where [s′] and [s″] denote the real and imaginary part of the compliance, respectively.
Equations (11.106a)–(11.106f) or (11.107) represent the compliances in the principle material (1‐2‐3) coordinate system. Moreover, the compliance can be used to extract the stiffness as the stiffness is inverse of the compliance; that is, [c*] = [s*]−1 and [cE] = [sE]−1.
In this section, the micromechanics model presented in Section 11.6.2 is utilized to predict the mechanical properties of a hybrid composite consisting of an epoxy matrix (Hercules 3501‐6) (Lesieutre et al. 1993), CB particles, and piezoelectric particles (Noliac® – PCM51). The physical properties of the three ingredients of the composite are listed in Table 11.3. The epoxy has a Poisson's ratio, v, of 0.38, and the CB particles are embedded with a volume fraction 0.20. It is assumed that the piezo‐particles have aligned orientations.
It should be mentioned that both the inclusions shape and shunting characteristics are unknown. Hence, in the model predictions the inclusions aspect ratio, α, and the shunt frequency, Ω, are investigated.
Figure 11.33a displays the damping ratios of the composite in direction 1 and 2, as function of the shunt frequency, Ω = R CTω for different values of the volume fraction vp of the piezo‐particles and for aspect ratio α = 100. The figures indicate that there are optimal values of ohm, that is, optimal shunting resistance, for each volume fractions vp at which the loss factor attain maximum values. It can be seen that these values of ohm are 1.68, 1.54, and 1.43 when the volume fractions vp are equal to 0.15, 0.35, and 0.55. The corresponding maximum loss factors are 9.67, 11.45, and 11.99, respectively. The values of the loss factors are normalized with respect to that of the epoxy matrix. This indicates that using the shunting of the piezo‐particles enhances the damping characteristics by almost an order of magnitude in both the one and two directions as the hybrid composite is transversally isotropic.
In direction 3, that is, during compression loading, the corresponding maximum loss factors are 0.81, 0.68, and 0.67 when the volume fractions vp are equal to 0.15, 0.35, and 0.55, respectively, as shown in Figure 11.33b. The corresponding optimal shunt frequencies ohm are 1.28, 1.17, and 1.035.
Figures 11.34 and 11.35 display the effects of the volume fraction of the piezo‐particles vp and aspect ratio α, respectively. The figures indicate that increasing the volume fraction and the aspect ratio result in a significant increase in the transverse loss factor. The opposite is true for the compressive loss factors.
The displayed results in Figures 11.33" through 11.35 indicate that the loss factors of the hybrid composite depend on the volume fraction, the orientation distribution, the aspect ratio (shape), and the shunt frequency‐resistance parameter. Table 11.5 summarizes a listing of the shunt frequency parameter ohm at which the loss factor ratio, η/ηm, attains maximum values for different orientation distributions and various aspect ratios. The listed values are calculated for piezo‐particles PCM51 with volume fraction of 35% and CB volume fraction of 25%.
Table 11.5 Values of Ωmax and ηmax for various orientations and aspect ratios when vc = 0.20 and vp = 0.35.
Inclusion aspect ratio (α) | Aligned orientation (θ0 → 0, ϕ0 → 0) |
2D random orientation (θ0 → π, ϕ0 → 0) |
3D random orientation (θ0 → π, ϕ0 → π/2) |
|||||||||||
0.0 | 1.16 | 0.62 | 1.13 | 2.36 | 1.75 | 3.71 | 1.16 | 2.42 | 2.03 | 4.41 | 1.78 | 3.71 | 2.02 | 4.42 |
0.5 | 1.07 | 1.18 | 1.06 | 0.92 | 1.03 | 1.23 | 1.08 | 0.92 | 1.16 | 0.71 | 1.01 | 1.10 | 1.10 | 0.71 |
1.0 | 1.07 | 1.53 | 1.07 | 0.84 | 1.11 | 1.37 | 1.07 | 0.84 | 1.34 | 0.75 | 1.06 | 1.23 | 1.28 | 0.75 |
100 | 1.58 | 11.45 | 1.13 | 0.68 | 2.12 | 8.45 | 1.16 | 0.68 | 2.43 | 7.52 | 2.05 | 8.47 | 2.44 | 7.52 |
1000 | 1.59 | 11.61 | 1.19 | 0.68 | 2.15 | 8.60 | 1.16 | 0.68 | 2.49 | 7.68 | 2.02 | 8.47 | 2.51 | 7.68 |
Figures 11.36–11.38 display the effect of the shunt frequency, piezo‐particle volume fraction, and piezo‐particles aspect ratio on the compliance coefficients of the hybrid composite, respectively.
This chapter has presented the prediction of the mechanical properties of composites consisting of a polymer matrix filled with nanoparticles as a function of different volumetric fractions and physical properties of the constituents. The presented prediction approaches are based on the MTM, which is rooted on Eshelby's equivalent inclusion technique, the classical filler reinforcement methods that are either analytical in nature, and finite element modeling using the RVE approach.
In this chapter, particular emphasis is placed on electrically conducting composites, such as CB/polymer composites, because of the possibility of integrating these composites with piezoelectric particles or elements in order to enhance their damping properties through appropriate electric shunting.
The transformation matrix shown in Eqs. (11.43) and (11.44) is defined by Wetherhold (1988):
with
where θ and φ are rotation Euler angles defined in Figure 11.3.
The transformation matrix [R2] shown in Eqs. (11.43) and (11.44) can be obtained by the relationship:
In this appendix, the reinforcement of a polymer by virtue of embedding solid particles in it is analyzed using the approach developed by Smallwood (1944) in an attempt to determine the effect of particle size and volume fraction on the modulus of the resulting composite.
Consider a rigid, spherical filler particle of radius R embedded in an isotropic polymer of known storage modulus as shown in Figure 11.B.1. The deflections and the resulting stresses near the particle are determined when the polymer is subjected to unidirectional tensile force T along the y direction.
The governing equilibrium equations as obtained from the theory of elasticity are given, for example by Barber (2010), as follows:
where u, v, w are the deflections along the x, y, z directions. Also, Δ denotes the cubical dilation (Δ = εxx + εyy + εzz) and ∇ is the Laplacian. Note that the equilibrium Eq. (11.B.1) are written in terms of the Lamé constants λ and μ that are related to Young's modulus E and Poisson's ratio, υ, such that (Reddy 2013):
Note that Lamé constants λ and μ are also given by:
where G is the shear modulus.
Equation (11.B.1) are subjected to the following boundary conditions:
and
The boundary conditions given by Eqs. (11.B.4) through (11.B.7), at r = ∞, indicate that far away from the particle, the deflections (u, v, w) are caused only by the tension T.
Solutions of Eq. (11.B.1) subject to the boundary conditions are as follows:
and
where , r2 = x2 + y2 + z2, and
Assuming that Poisson's ratio μ of the composite is close to that of the polymer, which is nearly 0.5. Eq. (11.B.2) can be rewritten as:
For such a limit (μ/λ ≅ 0), the parameters A, B, and C reduce to:
In order to quantify the effect of the filler particle, the strain energy function WT is determined as follows:
where N is the number of particles and V is the total volume of the polymer and the particles.
Then, the strain energy per unit volume is given by:
where v is the volume fraction of the particles.
Substituting Eqs. (11.B.3) and (11.B.12) into Eq. (11.B.14) yields:
The deflection w, at any point, due to the effect of all the filler particles is given by modifying Eq. (11.B.7) to take the following form:
where accounts for the effect of the particles and takes the following form (Smallwood 1944):
Then, the associated strain energy per unit volume is given by:
The equivalence of Eqs. (11.B.15) and (11.B.18) yields:
Expanding Eq. (11.B.19) into a first order Taylor series yields:
Song and Zheng (2016) summarized a list of other mathematical models that intend to improve the approximation of Smallwood Eq. (11.B.20) by adding higher order terms. In their summary, the coefficient of second order Taylor series expansion v2 is reported by to vary between 2.5 and 15.6 depending on the adopted mathematical model. In here, the approach developed by Smallwood (1944) and then modified by Guth (1945) yields the following commonly used reinforcement relationship:
Show that the electrical impedance Z of the circuit is given by:
where Z1 and Z2 denote the real and imaginary components of the impedance given by:
By eliminating the frequency ω from Z1 and Z2, show that the equation of the Cole–Cole plot, displayed in Figure P11.2, is given by:
where and .
The stiffness matrix CA of the polymer matrix A, in a 2D state of plane‐strain, is given by:
where and υA are the complex modulus and Poisson's ratio of the polymer, respectively. Assume that where EA is arbitrary and υA = 0.4995.
Show that the complex modulus stiffness matrix (CB) of component B corresponding to the arrangement shown in Figure P11.3 is given by:
The component f11 denotes the reaction force along the six and eight directions. It is determined by restraining the DOF 1, 2, 3, 4, 5, and 7 while the DOF 6 and 8 are subjected to unit displacements. This yields f11 as follows:
Also, f12 denotes the reaction force along the DOF 1 and 5, which is obtained by restraining the DOF 2, 3, 4, 6, 7, and 8 while subjecting the DOF 1 and 5 to unit displacements. This yields f12 as follows:
Finally, f33 denotes the reaction force along the DOF 5 and 7, which is obtained by restraining the DOF 1, 2, 3, 4, 6, and 8 while subjecting the DOF 5 and 7 to unit displacements. This yields f33 as follows:
where and
Assume that the stiffness matrix of the component B is set to represent an isotropic material, then, show that this implies that:
Let k22 = k21, k32 = k31, and c21 = c22 to ensure that s1 = s2 and let c21/k21 = 0.2E‐3, then show that if k21 = 8.45EA and k31 = − 0.304k21, the complex modulus of component B is given by:
Determine the equivalent dimensionless storage modulus E*/EA and loss factor η/ηA of the composite as function of the frequency and for volume fraction vB = 0.15 of the ingredient B. Compare the obtained results with those corresponding to positive k31 = 0.304 k21. Show that the storage moduli and loss factors in the 11, 12, and 33 directions are as displayed in Figure P11.4a–c, respectively. Observe that the storage moduli of the composite with negative k31 are lower than those of the composite with positive k31 whereas the loss factors of the composite with negative k31 are higher.
Note that for such a 2D composite, the associated Eshelby strain tensor is given (Marenić et al. 2015) as follows:
where
s1111 = A(3 + γ) = s2222, s1122 = A(1 − γ) = s2211, and s1212 = A(1 + γ) with A = 1/[8(1 − υ)] and γ = 2(1 − 2υ).
where α = H/L and β = h/L.
Assuming Young's modulus EA and Poisson's ratio υA of the polymer are such that EA = 10 MPa and υA = 0.49, whereas Young's modulus EB and Poisson's ratio υB of the inclusion are such that EB = 210 GPa and υB = 0.30. If α = 1, β = 1/15, and θ = 30°, then determine the equivalent dimensionless moduli cij/EA of the composite using ANSYS analysis of the RVE as shown in Figure P11.7. Show that c11/EA = 604.69, c21/EA = 96.35, c31/EA = 81.09, c12/EA = 5.36, c22/EA = 5.66, c32/EA = 5.39, c13/EA = 5.03, c23/EA = 4.98, and c33/EA = 5.40.
Assuming Young's modulus EA and Poisson's ratio υA of the polymer are such that EA = 10 MPa and υA = 0.49, whereas Young's modulus EB and Poisson's ratio υB of the inclusion are such that EB = 210 GPa and υB = 0.30. Note that the mass‐in‐inclusion is made of steel with diameter 20 nm. If α = 1, β = 1/15, and θ = 30°, then determine the equivalent dimensionless moduli cij/EA of the composite using ANSYS analysis of the RVE as shown in Figure P11.9. Show that c11/EA = 1294.24, c21/EA = 335.84, c31/EA = 190.57, c12/EA = 5.36, c22/EA = 5.66, c32/EA = 5.39, c13/EA = 5.03, c23/EA = 4.98, and c33/EA = 5.40.
Assuming Young's modulus EA and Poisson's ratio υA of the polymer are such that EA = 10 MPa and υA = 0.49, whereas Young's modulus EB and Poisson's ratio υB of the inclusion are such that EB = 210 GPa and υB = 0.30. If α = 1, β = 1/15, and θ = 30°, then determine the equivalent dimensionless moduli cij/EA of the composite using ANSYS analysis of the RVE as shown in Figure P11.11. Show that c11/EA = 555.94, c21/EA = 84.01, c31/EA = 93.15, c12/EA = 5.35, c22/EA = 5.66, c32/EA = 5.39, c13/EA = 5.05, c23/EA = 4.97, and c33/EA = 5.44.
Assuming Young's modulus EA and Poisson's ratio υA of the polymer are such that EA = 10 MPa and υA = 0.49, whereas Young's modulus EB and Poisson's ratio υB of the inclusion are such that EB = 210 GPa and υB = 0.30. Note that the mass‐in‐inclusion is made of steel with diameter 20 nm. If α = 1, β = 1/15, and θ = 30°, then determine the equivalent dimensionless moduli cij/EA of the composite using ANSYS analysis of the RVE as shown in Figure P11.13. Show that c11/EA = 630.04, c21/EA = 102.78, c31/EA = 85.15, c12/EA = 5.33, c22/EA = 5.64, c32/EA = 5.25, c13/EA = 5.01, c23/EA = 4.93, and c33/EA = 5.38.