Damage Mechanics of Carbon Nanotubes
George Z. Voyiadjis ^{1}, Peter I. Kattan ^{2,*}

Boyd Professor, Department of Civil and Environmental Engineering, Louisiana State University, Baton Rouge, LA 70803, USA

Visiting Professor, Department of Civil and Environmental Engineering, Louisiana State University, Baton Rouge, LA 70803, USA
* Correspondence: George Z. Voyiadjis
Academic Editor: Zhao Yang Dong
Special Issue: Nanostructured Materials for Clean & Sustainable Energy Storage & Efficiency for Wellness & Green Buildings
Received: October 27, 2020  Accepted: April 28, 2021  Published: May 11, 2021
Journal of Energy and Power Technology 2021, Volume 3, Issue 2, doi:10.21926/jept.2102021
Recommended citation: Voyiadjis GZ, Kattan PI. Damage Mechanics of Carbon Nanotubes. Journal of Energy and Power Technology 2021;3(2):17; doi:10.21926/jept.2102021.
© 2021 by the authors. This is an open access article distributed under the conditions of the Creative Commons by Attribution License, which permits unrestricted use, distribution, and reproduction in any medium or format, provided the original work is correctly cited.
Abstract
A robust mathematical method for the characterization of damage in carbon nanotubes is presented the presentation here is limited to elasticity. In this regard, the second and third order elastic stiffnesses are employed. All this is based on damage mechanics. The hypotheses of elastic strain equivalence and elastic energy equivalence are utilized. A new damage variable is proposed that is defined in terms of the surface area. This is in contrast to the classical damage variable which is defined in terms of the crosssectional area. In the presentation, both the onedimensional case (scalars) and the threedimensional case (tensors) are illustrated.
Keywords
Carbon nanotubes; damage mechanics; material modelling; elasticity; elastic energy equivalence; elastic strain equivalence; damage variable
1. Introduction
Carbon nanotubes have been used recently in worldwide applications. Their research is taking off at exponential speeds [1,2,3,4,5]. Recent work has been done also on grapheme [6,7,8,9,10]. The aim of this presentation is to provide a way to characterize damage inc carbon nanotubes. The authors in this respect utilize the principles of continuum damage mechanics. The topic of continuous damage mechanics was first proposed by Kachanov [11] and further developed by Lee et al. [12], Voyaadjis and Kattan [13,14,15,16], Sidoroff [17], Krajcinovic [18], and others; and Kattan and Voyadjis [19,20,21].
In 1958, Kachanov [11] first proposed the concept of effective stress and introduced the theme of continuous damage mechanics. Rabotnov[22] and later others [14,15,16,20,21,23,24] followed closely behind. In this framework, a damage variable is introduced that is defined in terms of the crosssectional area. This damage variable has a zero value for the virgin undamaged material and a value of 1 for the case of complete rupture (total damage).
Research in damage mechanics has progressed very quickly in the past 50 years [12,13,17,19,25,26,27,28,29,30,31,32,33,34,35]. In addition, there are also some noteworthy researches on the corresponding topics of healing mechanics [24,36,37,38,39]. Some recent research work is considered to be a combined damage/repair model for different types of materials. Basaran has introduced a new way to define damage based on entropy generation [40,41]. This work consists of two main parts, which deal with the complete theoretical characteristics of carbon nanotube damage within the framework of continuous damage mechanisms. More recent work has been done on damage mechanics especially for rocks [42,43,44,45,46,47,48].
Onedimenstional damage in carbon nanotubes is presented in Section 1 in terms of scalars. To illustrate this method, both the hypothesis of elastic strain equivalence and the hypothesis of elastic energy equivalence have been utilized. Especially for nanomaterials like carbon nanotubes, a new damage variable is introduced that is defined based on the surface area as the surface area plays a major role in these types of materials.
Threedimenstional damage in carbon nanotubes is illustrated in Section 2 in terms of tensor. For this purpose, the two hypotheses of elastic strain equivalence and elastic energy equivalence are utilized. Furthermore, the damage tensor is defined in terms of the surface area. This damage tensor is a generalization of the new damage variable that was introduced in Section1 previously.
2. OneDimensional Formulation (Scalars)
The elastic stressstrain formula for carbon nanotubes is given by the following formula for the onedimensional case using scalars (uniaxial tension) [4]:
\[\sigma = E\epsilon + D \epsilon^{2} \tag{1} \]
where σ and ε are the stress and strain, respectively, while E and D are the secondorder elastic stiffness and thirdorder elastic stiffness, respectively. The purpose of this work is to show the damage process and the transition of the damage process within the framework of continuous damage mechanics.
In Continuum Damage Mechanics, the effective stress σ is given by:
\[\bar{\sigma} = \frac{\sigma}{1\phi} \tag{2} \]
where φ is the classical damage variable. This damage variable is defined in terms of the loss of crosssectional area due to damage. The effective stress is defined as acting in the virtual undamaged structure, while the actual stress acts in the damaged structure of the material. The value of the damage variable ranges from 0 (undamaged material state) to 1 (completely broken).
The purpose of this work is to show the material constants E and D conversion methods in the damage process, that is, to show the two relations of the following form $\overline{E}=E{f}_{\left(\varphi \right)}$ and $\overline{D}=D{g}_{\left(\varphi \right)}$, where f and g are unspecified functions. The authors will show the explicit expressions of the scalar functions where $\overline{E}$ and $\overline{D}$ are the secondorder effective lossless elastic stiffness and the thirdorder effective lossless elastic stiffness, respectively. For this purpose, two classical assumptions of elastic strain equivalent and elastic energy equivalent will be used.
2.1 Scalar Formula Using Elastic Strain Equivalence Assumption
The first one starts with the assumption of elastic strain equivalence. In this assumption, it is assumed that the elastic strain is the same in the deformed/damaged configuration and the virtual/undamaged configuration. This assumption can be written as:
\[\bar{\epsilon} = \epsilon \tag{3} \]
Write the formula (1) in the effective/undamaged configuration for the barred quantities, as shown below:
\[\bar{\sigma} = \bar{E} \bar{\epsilon} + \bar{D} \bar{\epsilon}^{2} \tag{4} \]
Then substitute the strain in equation (3) into equation (4), and substitute the stress in equation (2) into equation (4) to get:
\[\frac{\sigma}{1\phi} = \bar{E}\bar{\epsilon} + \bar{D} \bar{\epsilon}^{2} \tag{5} \]
Comparing equations (1) and (5), the following conversion equations for the two elastic stiffnesses of carbon nanotubes can be immediately obtained:
\[E = \bar{E}\left(1\phi\right) \tag{6} \] \[D = \bar{D}\left(1\phi\right) \tag{7} \]
It should be noted from equations (6) and (7) that the two elastic stiffness conversions involve the same damage variable. For more elaborate models, two different and independent damage variables should be usedone for each elastic stiffness. This avenue of research will be discussed in Section 1.3.
2.2 Scalar Formula Using Elastic Energy Equivalence Assumption
The sequel discusses how the two elastic stiffnesses can be converted under the assumption of elastic energy equivalent. In this assumption, it is assumed that the elastic energy is the same between the deformed/damaged configuration and the virtual/undamaged configuration. The mathematical formula of this hypothesis is as follows:
\[\bar{U} = U \tag{8} \]
where U is the elastic energy in the deformed/damaged configuration and $\overline{U}$ is the elastic energy in the fictitious/undamaged configuration.
Using equation (1), in the deformed/damaged configuration, the strain energy of carbon nanotubes is obtained as follows:
\[U = \int\sigma d \epsilon = \int \left( E\epsilon + D \epsilon^{2} \right) d\epsilon = \frac{1}{2} E \epsilon^{2} \frac{1}{3} D \epsilon^{3} + c_{1} \tag{9} \]
where C_{1} is a constant of integration to be determined. Similarly one can show that the elastic energy for carbon nanotubes can be written as follows in the fictitious/undamaged configuration:
\[\bar{U} = \frac{1}{2} \bar{E} \bar{\epsilon}^{2} + \frac{1}{3} \bar{D} \bar{\epsilon}^{3} + c_{2} \tag{10} \]
where C_{2} is a constant of integration to be determined. Substituting equations (9) and (10) into equation (8), and utilizing the fact that c_{1} = 0 and c_{2} = 0 based on the initial conditions in each configuration, one obtains:
\[\frac{1}{2} \bar{E} \bar{\epsilon}^{2} + \frac{1}{3} \bar{D} \bar{\epsilon}^{3}=\frac{1}{2} E \epsilon^{2} + \frac{1}{3} D \epsilon^{3} \tag{11} \]
In the next formula, one rewrites equation (4) in the following form:
\[\bar{E} = \frac{\bar{\sigma}}{\bar{\epsilon}}  \bar{D} \bar{\epsilon} \tag{12} \]
Similarly equation (1) is rewritten as follows:
\[E = \frac{\sigma}{\epsilon}  D \epsilon \tag{13} \]
Substituting equations (12) and (13) into equation (11) and simplifying, we can obtain:
\[\frac{1}{2} \bar{\sigma} \bar{\epsilon}  \frac{1}{6} \bar{D} \bar{\epsilon}^{3} = \frac{1}{2} \sigma \epsilon  \frac{1}{6} D \epsilon^{3} \tag{14} \]
If the simplified assumptions expressed in the following two equations are not made based on the above equation (14), it cannot be further proceeded:
\[\frac{1}{2} \bar{\sigma} \bar{\epsilon} = \frac{1}{2} \sigma \epsilon \tag{15} \] \[\frac{1}{6} \bar{D} \bar{\epsilon}^{3} = \frac{1}{6} D \epsilon^{3} \tag{16} \]
People immediately recognized the expression in equation (15) as a hypothesis of the elastic energy equivalence of linear elastic materials. Next, it is performed by substituting equation (2) into equation (15). After simplifying the obtained equation, the following elastic strain conversion equation is obtained:
\[\bar{\epsilon} = \epsilon (1\phi) \tag{17} \]
The next equation substitutes equation (17) into equation (16). After simplifying the obtained equation, the following thirdorder elastic stiffness conversion equation can be obtained:
\[D = \bar{D} (1\phi)^{3} \tag{18} \]
Next, substitute equations (17) and (18) into equation (11). After simplifying the obtained equation, the conversion equation of elastic stiffness can be obtained, as shown below:
\[E = \bar{E} (1\phi)^{2} \tag{19} \]
Therefore, it can be seen from the two conversion equations (18) and (19) that the two elastic stiffnesses of carbon nanotubes are converted in two completely different ways. However, the same single damage variable is still used in both cases. The sequel explores the possibility of using two different and independent damage variables for the two elastic stiffnesses of carbon nanotubes.
2.3 Elastic Stiffness Degradation Damage Variables and Surface Area Damage Variables
In addition, two other scalar damage variables will be studied. The first is the damage variable defined by the degradation of elastic stiffness, and the second is the damage variable defined by the thirdorder degradation of elastic stiffness [32,49,50]. These two damage variables will be compared with the classical damage variables.
One attempts to define two new distinct and independent damage variables $\ell $ and m defined as follows in terms of elastic stiffness reduction:
\[\ell = \frac{\bar{E} E}{E} \tag{20}\]
From the above two equations, you can immediately get:
\[ \bar{E}=E\,\left( 1+\ell \right) \tag{22}\] \[\bar{D}=D\,\left( 1+m \right) \tag{23}\]
From equations (20) and (21) it appears that the maximum values of the two damage variables ℓ and m are infinity. This is because that the values of E and D tend to zero. However, in practical applications, the stiffness does not reduce that much. In addition, if one looks at equation (2) for the effective stress, one notices that when the value of the damage variable Ⴔ is 1, the value of the effective stress explodes to infinity. Thus the damage variable will approach 1 as its maximum value. Again, it is seen from equation (2) that any value of the damage variable that is beyond 1 is meaningless. The maximum value of Ⴔ is 1. By analogy, the authors limit the maximum values of ℓ and m to 1 also.
First, use the elastic strain equivalent assumption proposed in the transformation equation (6). One substitutes equation (22) into equation (6). After simplifying the obtained equation, we can obtain:
\[ \left( 1\phi \right) \left( 1+\ell \right)=1 \tag{24} \]
The next equation (23) is substituted into equation (7). After simplifying the obtained equation, we can obtain:
\[ \left( 1\phi \right)\left( 1+m \right)=1 \tag{25} \]
Comparing equations (24) and (25), one can immediately conclude that for the special case of the elastic strain equivalence assumption, the two independent damage variables are the same. The conclusion is as follows:
\[ \ell =m \tag{26} \]
The next step is to explore the nature of the relationship between the two independent damage variables when using the elastic energy equivalent assumption, as shown in the transformation equations (18) and (19). This is a repeat of the previous derivation but for the hypothesis of elastic energy equivalence. For this purpose, equation (22) is substituted into equation (19). After simplifying the obtained equation, we can obtain:
\[ {{\left( 1\phi \right)}^{2}}\,\left( 1+\ell \right)=1 \tag{27} \]
Similarly, substituting equation (23) into equation (18) and simplifying the resulting expression, we can obtain:
\[ {{\left( 1\phi \right)}^{3}}\,\left( 1+m \right)=1 \tag{28} \]
Comparing the two expressions in equations (27) and (28), some simple algebraic operations can be performed to obtain the following relationship between the two independent damage variables:
\[ {{\left( 1+\ell \right)}^{\frac{1}{2}}}={{\left( 1+m \right)}^{\frac{1}{3}}} \tag{29} \]
Solving equation (29) for in terms of , one obtains:
\[ \ell ={{\left( 1+m \right)}^{\frac{2}{3}}}1 \tag{30} \]
It is immediately seen from equation (30) that both damage variables $\ell $ and m vanish for the virgin/undamaged material. However, when the value of m reaches 1, then the value of $\ell $ reaches the value 0.587. This means that the maximum damage to the elastic stiffness E is only about half of the maximum damage to the thirdorder elastic stiffness D. The relationship between these two damage variables of elastic stiffness degradation is illustrated further by plotting the expression of equation (30) as shown in Figure 1.
Figure 1 Relationship between the two Damage Variables l and m.
One can get further insight into the relationship between the two damage variables $\ell $ and m by writing an approximation to equation (30) that is applicable to small values of damage. This can be accomplished easily by writing the Taylor series expansion of the power expression in parenthesis of equation (30). The following Taylor series expansion will be utilized:
\[ ~{{\left( 1+m \right)}^{\frac{2}{3}}}\approx 1+\frac{2}{3}m\frac{1}{9}{{m}^{2}}+\ldots \ldots . \tag{31} \]
Only substituting the first two terms on the right side of equation (31) into equation (30), we can get:
\[ l\approx \frac{2}{3}m \tag{32} \]
Thus it is seen from equation (32) that for small values of damage, the reduction in the elastic stiffness E is about twothirds of the reduction in the thirdorder elastic stiffness D. It is also noted from Figure 1 that for higher values of damage (approaching 1), the reduction in the elastic stiffness E reduces to about onehalf of the reduction in the thirdorder stiffness D. This fact can also be noted from Figure 1 in the relevant applicable range of values.
In nanomaterials, surface area is critical. Therefore, it is important to introduce a new damage variable, which is defined in terms of surface area reduction. To this end, the new damage variable (called the surface area damage variable) is defined as follows:
\[ {{\phi }_{s}}=\frac{S\bar{S}}{S} \tag{33} \]
where S is the surface area in the deformed/damage configuration while $\overline{S}$ is the corresponding surface area in the fictitious/undamaged configuration.
If one looks at equation (2) for the effective stress, one deduces immediately that the value of the crosssectional area cannot be 1 or more. The effective stress explodes to infinity when Ⴔ is equal to 1. Thus, its maximum value is limited to unity. The same argument holds for the surface area damage variable. Although the equation for the effective stress in terms of the surface area damage variable is not shown in the manuscript, the effective stress also explodes to infinity when the surface area damage variable approaches 1. Thus, its value is limited to unity also.
Next one writes the following geometric relationships that are applicable for Euclidean shapes/solids [51]:
\[ S=P\cdot t \tag{34} \] \[ A=\alpha \cdot {{P}^{2}} \tag{35} \]
where A is the crosssectional area, S is the surface area, P is the perimeter, t is the thickness, and α is a shape constant (e.g. $\alpha =\frac{1}{16}$ for a square and $\alpha =\frac{1}{4\pi}$ for a circle).
Starting with the definition of the classical crosssectional area damage variable as follows:
\[ \phi =\frac{A\bar{A}}{A} \tag{36} \]
one substitutes for A from equation (35) into equation (36), and for $\overline{A}$ a similar expression but using barred quantities to obtain:
\[ \phi =\frac{\alpha \,{{P}^{2}}\alpha \,{{{\bar{P}}}^{2}}}{\alpha \,{{P}^{2}}} \tag{37} \]
Note from equation (37) that it is assumed that the shape constant is not damaged. This is a reasonable assumption and can simplify the equation well. Formula (37) is simplified as follows:
\[ \phi =1{{\left( \frac{{\bar{P}}}{P} \right)}^{2}} \tag{38} \]
Equivalently the above equation is rewritten as follows:
\[ \frac{{\bar{P}}}{P}=\sqrt{1\phi } \tag{39} \]
Next one substitutes equation (34) for S into equation (33), and similar equation for $\overline{S}$ but with barred quantities to obtain:
\[ {{\phi }_{s}}=\frac{P\cdot t\bar{P}\cdot t}{P\cdot t} \tag{40} \]
It is assumed that no damage occurs along the thickness t. In order to draw conclusions and simplify the equation, this is a reasonable and necessary assumption. Therefore, equation (40) is simplified as follows:
\[ \frac{{\bar{P}}}{P}=1{{\phi }_{s}} \tag{41} \]
Finally, compare equations (39) and (41) to arrive at the desired relationship:
\[ 1{{\phi }_{s}}=\sqrt{1\phi } \tag{42} \]
The above expression is the relationship between the crosssectional area damage variable φ and the surface area damage variable φ_{s}. This relationship is rewritten in its final form as follows based on equation (42):
\[ {{\phi }_{s}}=1\sqrt{1\phi } \tag{43} \]
It can be seen from equation (43) that for the original/undamaged material, both the crosssectional damage variable and the surface area damage variable disappear. Similarly, both variables reach the maximum value of 1 at the same time. This can be seen by plotting the relationship of equation (43), as shown in Figure 2.
Figure 2 Relationship between the CrossSectional Area Damage Variable and the Surface Area Damage Variable.
Indeed, the presentation is mainly theoretical but rigorous. These are two of the main features of this work. The figures are indicative of this issue.
3. ThreeDimensional Formulation (Tensors)
Finally, the whole process is extended to (using tensors) the threedimensional deformation and damage state of carbon nanotubes. Classic indicial symbols will be used to represent tensors. For this, the following generalization of equation (1) will be used [4].
\[ {{\sigma }_{ij}}={{E}_{ijkl}}\,{{\varepsilon }_{kl}}+\frac{1}{2}{{D}_{ijklmn}}\,{{\varepsilon }_{kl}}\,{{\varepsilon }_{mn}} \tag{44} \]
where σ_{ij} and ε_{kl} are the components of the secondrank stress and secondrank strain tensors, respectively. In equation (44), E_{ijkl} and D_{ijklmn} are the fourthrank secondorder elasticity tensor and the sixthrank thirdorder elasticity tensor, respectively. Tensorial generalizations of the transformation equations (6), (7), (18), and (19) will be formulated. Note that by comparing the tensorial equation (44) with the scalar equation (1), one notes that $E\equiv {E}_{1111}$ and $D\equiv E{D}_{1111}/2$.
Based on equation (44) and using the following derivative derived from the foundation of solid mechanics:
\[ {{\sigma }_{ij}}=\frac{dU}{d{{\varepsilon }_{ij}}} \tag{45} \]
where U is the strain energy in the deformed/damaged configuration, one obtains the stresses. Substituting equation (45) into equation (44) and solving for U(by integrating the relevant expression), one obtains:
\[ U=\frac{1}{2}{{E}_{ijkl}}{{\varepsilon }_{ij}}{{\varepsilon }_{kl}}+\frac{1}{6}{{D}_{ijklmn}}{{\varepsilon }_{ij}}{{\varepsilon }_{kl}}{{\varepsilon }_{mn}} \tag{46} \]
Next one derives the transformation equations for the fourthrank secondorder elasticity tensor E_{ijkl} and the sixthrank thirdorder elasticity tensor D_{ijklmn} twice – one time using the hypothesis of elastic strain equivalence and another time using the hypothesis of elastic energy equivalence.
3.1 Tensor Formula Using Elastic Strain Equivalence Assumption
Start with the assumption of elastic strain equivalence. The strain tensor is represented in the following form:
\[ {{{\bar{\varepsilon }}}_{ij}}={{\varepsilon }_{ij}} \tag{47} \]
The generalization of the effective stress conversion equation (equation (2)) is usually written in the following form within the framework of continuous damage mechanics [14,15,17].
\[ {{{\bar{\sigma }}}_{ij}}={{M}_{ijkl}}{{\sigma }_{kl}} \tag{48} \]
where ${\overline{\sigma}}_{ij}$ is the effective stress tensor (defined in the fictitious/undamaged configuration) and M_{ijkl} is the fourthrank damage effect tensor.
Next, use the barred quantity to write the nonlinear elastic constitutive equation of carbon nanotubes in virtual/lossless configuration (44):
\[ {{{\bar{\sigma }}}_{ij}}={{{\bar{E}}}_{ijkl}}\,{{{\bar{\varepsilon }}}_{kl}}+\frac{1}{2}{{{\bar{D}}}_{ijklmn}}\,{{{\bar{\varepsilon }}}_{kl}}\,{{{\bar{\varepsilon }}}_{mn}} \tag{49} \]
Substituting for the effective stress tensor ${\overline{\sigma}}_{ij}$ from equation (48) into equation (49), one obtains:
\[ {{M}_{ijkl}}{{\sigma }_{kl}}={{{\bar{E}}}_{ijkl}}\,{{{\bar{\varepsilon }}}_{kl}}+\frac{1}{2}{{{\bar{D}}}_{ijklmn}}\,{{{\bar{\varepsilon }}}_{kl}}\,{{{\bar{\varepsilon }}}_{mn}} \tag{50} \]
Premultiplying equation (50) by ${M}_{\left(pqij\right)}^{(1)}$, carrying out the relevant multiplications, contractions, and simplifications, one arrives at the following equation:
\[ {{\sigma }_{pq}}=M_{pqij}^{1}{{{\bar{E}}}_{ijkl}}{{{\bar{\varepsilon }}}_{kl}}+\frac{1}{2}M_{pqij}^{1}{{{\bar{D}}}_{ijklmn}}{{{\bar{\varepsilon }}}_{kl}}{{{\bar{\varepsilon }}}_{mn}} \tag{51} \]
Substituting equation (47) into equation (51) (this means replacing the barred strain tensor components with unbarred strain tensor components), one obtains:
\[ {{\sigma }_{pq}}=M_{pqij}^{1}{{{\bar{E}}}_{ijkl}}{{\varepsilon }_{kl}}+\frac{1}{2}M_{pqij}^{1}{{{\bar{D}}}_{ijklmn}}{{\varepsilon }_{kl}}{{\varepsilon }_{mn}} \tag{52} \]
Finally, the stress tensor equations (44) and (52) are compared, and two conversion equations for the fourth and sixth order elastic tensors are obtained:
\[ {{E}_{pqkl}}=M_{pqij}^{1}{{{\bar{E}}}_{ijkl}} \tag{53} \] \[ {{D}_{pqklmn}}=M_{pqij}^{1}{{D}_{ijklmn}} \tag{54} \]
The above two expressions represent the conversion equations of the fourth and sixth order elastic tensors of carbon nanotubes under the assumption that the elastic strains are equivalent. It can be seen from the above two equations that both elastic tensors are transformed in exactly the same way using the same transformation tensor ${M}_{\left(pqij\right)}^{(1)}$ in the two expressions. This is not the case when using more complex assumptions of elastic energy equivalence, as shown below.
3.2 Tensor Formula Using Elastic Energy Equivalence Assumption
Next, using the assumption of elastic energy equivalence, the conversion equations of the fourth and sixth order elastic tensors of carbon nanotubes are derived. The mathematical formulation of this hypothesis is as follows:
\[ \bar{U}=\int {{{\bar{\sigma }}}_{ij\,}}\,d{{{\bar{\varepsilon }}}_{ij}}=\int {{\sigma }_{ij\,}}\,d{{\varepsilon }_{ij}}=U \tag{55} \]
where U is the elastic energy in the deformed/damaged configuration and $\overline{U}$ is the elastic energy in the fictitious/undamaged configuration.
Using equation (44), in the deformed/damaged configuration, the strain energy of carbon nanotubes is obtained as follows:
\[ U=\int {{\sigma }_{ij}}\,d{{\varepsilon }_{ij}}=\int \left( {{E}_{ijkl}}\,{{\varepsilon }_{kl}}+\frac{1}{2}{{D}_{ijklmn}}\,{{\varepsilon }_{kl}}\,{{\varepsilon }_{mn}}\, \right)d{{\varepsilon }_{ij}} \tag{56} \]
Carrying out the above tensorial integration, and using the relations $d\left({\u03f5}_{ij}{\u03f5}_{kl}\right)=2{\u03f5}_{kl}d{\u03f5}_{ij}$ and $d\left({\u03f5}_{ij}{\u03f5}_{kl}{\u03f5}_{kmn}\right)=3{\u03f5}_{kl}{\u03f5}_{kmn}d{\u03f5}_{ij}$, one obtains:
\[ U=\frac{1}{2}{{E}_{ijkl}}\,{{\varepsilon }_{ij}}\,{{\varepsilon }_{kl}}+\frac{1}{6}{{D}_{ijklmn}}\,{{\varepsilon }_{ij}}\,{{\varepsilon }_{kl}}\,{{\varepsilon }_{mn}} \tag{57} \]
Set the integral constant to zero according to the initial conditions. Similarly, we can prove that the elastic energy of carbon nanotubes can be written as follows in a virtual/undamaged configuration (using the same form of equation (57), but with a barred quantity):
\[ \bar{U}=\frac{1}{2}{{{\bar{E}}}_{ijkl}}\,{{{\bar{\varepsilon }}}_{ij}}\,{{{\bar{\varepsilon }}}_{kl}}+\frac{1}{6}{{{\bar{D}}}_{ijklmn}}\,{{{\bar{\varepsilon }}}_{ij}}\,{{{\bar{\varepsilon }}}_{kl}}\,{{{\bar{\varepsilon }}}_{mn}} \tag{58} \]
The next one starts with the expression of the effective stress tensor of equation (49) and rewrites it in the following form:
\[ {{{\bar{E}}}_{ijkl}}\,{{{\bar{\varepsilon }}}_{kl}}={{{\bar{\sigma }}}_{ij}}\frac{1}{2}{{{\bar{D}}}_{ijklmn}}\,{{{\bar{\varepsilon }}}_{kl}}\,{{{\bar{\varepsilon }}}_{mn}} \tag{59} \]
Postmultiplying equation (59) by ${\overline{\u03f5}}_{\left(lp\right)}^{(1)}$ and carrying out the needed algebraic manipulations, one arrives at:
\[ {{{\bar{E}}}_{ijpl}}=\frac{1}{2}{{{\bar{D}}}_{ijrsmn}}\,{{{\bar{\varepsilon }}}_{rs}}\,{{{\bar{\varepsilon }}}_{mn}}\,\bar{\varepsilon }_{lp}^{1}{{{\bar{\sigma }}}_{ij}}\,\bar{\varepsilon }_{lp}^{1} \tag{60} \]
Equation (60) can be written in the deformed/damaged configuration using unbarred quantities as follows:
\[ {{E}_{ijpl}}=\frac{1}{2}{{D}_{ijrsmn}}\,{{\varepsilon }_{rs}}\,{{\varepsilon }_{mn}}\,\varepsilon _{lp}^{1}{{\sigma }_{ij}}\,\varepsilon _{lp}^{1} \tag{61} \]
Substituting equations (60) and (61) together with equations (57) and (58) into the elastic energy equivalence hypothesis (550), and performing the required algebraic operations, we can get:
\[ \frac{1}{4}{{{\bar{D}}}_{ijrsmn}}\,{{{\bar{\varepsilon }}}_{rs}}\,{{{\bar{\varepsilon }}}_{mn}}\,{{{\bar{\varepsilon }}}_{ij}}\,{{I}_{ll}}\frac{1}{2}{{{\bar{\sigma }}}_{ij}}\,{{{\bar{\varepsilon }}}_{ij}}{{I}_{ll}}+\frac{1}{6}{{{\bar{D}}}_{ijklmn}}\,{{{\bar{\varepsilon }}}_{ij}}\,{{{\bar{\varepsilon }}}_{kl}}\,{{{\bar{\varepsilon }}}_{mn}} \] \[= \frac{1}{4}{{D}_{ijrsmn}}\,{{\varepsilon }_{rs}} {{\varepsilon }_{mn}}\,{{\varepsilon }_{ij}}\,{{I}_{ll}}\frac{1}{2}{{\sigma }_{ij}}\,{{\varepsilon }_{ij}}{{I}_{ll}}+\frac{1}{6}{{D}_{ijklmn}}\,{{\varepsilon }_{ij}}\,{{\varepsilon }_{kl}}\,{{\varepsilon }_{mn}} \tag{62} \]
where I_{ll} is related to the secondrank identity tensor ${I}_{il}\equiv {\delta}_{ij}$.
In order to proceed further and simplify the equation, it is necessary to assume that equation (62) can be decomposed into the following two equivalent equations:
\[ \frac{1}{2}{{{\bar{\sigma }}}_{ij}}\,{{{\bar{\varepsilon }}}_{ij}}\,{{I}_{ll}}=\frac{1}{2}{{\sigma }_{ij}}\,{{\varepsilon }_{ij}}\,{{I}_{ll}} \tag{63} \] \[ \frac{5}{12}{{{\bar{D}}}_{ijrsmn}}\,{{{\bar{\varepsilon }}}_{rs}}{{{\bar{\varepsilon }}}_{mn}}{{{\bar{\varepsilon }}}_{ij}}=\frac{5}{12}{{D}_{ijrsmn}}\,{{\varepsilon }_{rs}}{{\varepsilon }_{mn}}{{\varepsilon }_{ij}} \tag{64} \]
where I_{ll} = 3 is used to obtain equation (64) above.
The expression represented by equation (63) is immediately regarded as the form of the elastic energy equivalence hypothesis for linear elastic materials. Substituting formula (48) into formula (63), we can obtain:
\[ {{M}_{ijkl}}\,{{\sigma }_{kl}}\,{{{\bar{\varepsilon }}}_{ij}}={{\sigma }_{ij}}\,{{\varepsilon }_{ij}} \tag{65} \]
Premultiplying equation (65) by ${\sigma}_{\left(pq\right)}^{(1)}$ and simplifying, one obtains:
\[ {{M}_{ijpq}}\,{{{\bar{\varepsilon }}}_{ij}}={{\varepsilon }_{pq}} \tag{66} \]
Solving the effective strain tensor equation (66), the transformation equation of the strain tensor can be obtained, as shown below:
\[ {{{\bar{\varepsilon }}}_{ij}}=M_{ijpq}^{1}\,{{\varepsilon }_{pq}} \tag{67} \]
Next, continue to derive the transformation equations of the two elastic tensors. Substituting the effective elastic strain tensor of equation (67) into equation (64), we can obtain:
\[ {{{\bar{D}}}_{ijrsmn}}\,M_{rspq}^{1}\,{{\varepsilon }_{pq}}\,M_{mnkl}^{1}\,{{\varepsilon }_{kl}}\,M_{ijab}^{1}\,{{\varepsilon }_{ab}}={{D}_{ijrsmn}}\,{{\varepsilon }_{rs}}{{\varepsilon }_{mn}}\,{{\varepsilon }_{ij}} \tag{68} \]
By performing the required tensor manipulation and contraction to solve equation (68), the conversion equation of the sixthlevel thirdorder elastic tensor will be obtained:
\[ {{D}_{abpqkl}}={{{\bar{D}}}_{ijrsmn}}\,M_{rspq}^{1}\,M_{mnkl}^{1}\,M_{ijab}^{1} \tag{69} \]
All that remains now is to find the transformation equation for the fourthlevel elasticity tensor. To this end, formula (48) can be rewritten as:
\[ M_{ijmn}^{1}\,{{{\bar{\sigma }}}_{ij}}={{\sigma }_{mn}} \tag{70} \]
Next, substituting equations (69) and (70) into equation (44), we get:
\[ M_{cdji}^{1}\,{{{\bar{\sigma }}}_{cd}}={{E}_{ijkl}}\,{{\varepsilon }_{kl}}+\frac{1}{2}{{{\bar{D}}}_{abpqrs}}\,M_{abij}^{1}\,M_{pqkl}^{1}\,M_{rsmn}^{1}\,{{\varepsilon }_{kl}}\,{{\varepsilon }_{mn}} \tag{71} \]
Premultiplying equation (71) by M_{efil} and solving for the effective stress tensor, one obtains:
\[ {{{\bar{\sigma }}}_{ef}}={{M}_{efij}}\,{{E}_{ijkl}}\,{{\varepsilon }_{kl}}+\frac{1}{2}{{{\bar{D}}}_{abpqrs}}\,{{I}_{efab}}\,M_{pqkl}^{1}\,M_{rsmn}^{1}\,{{\varepsilon }_{kl}}\,{{\varepsilon }_{mn}} \tag{72} \]
where I_{efab} represents the fourthrank identity tensor. Next, rewrite equation (67) in the following form:
\[ {{\varepsilon }_{rs}}={{M}_{rsij}}\,{{{\bar{\varepsilon }}}_{ij}} \tag{73} \]
Substituting equation (73) into equation (72) and simplifying the resulting expression after some tensor manipulation, we get:
\[ {{{\bar{\sigma }}}_{ef}}={{M}_{efij}}\,{{E}_{ijkl}}\,{{M}_{klpq}}\,{{{\bar{\varepsilon }}}_{pq}}+\frac{1}{2}{{{\bar{D}}}_{efabcd}}\,{{{\bar{\varepsilon }}}_{ab}}\,{{{\bar{\varepsilon }}}_{cd}} \tag{74} \]
Comparing the equation (74) of the effective stress tensor with the original equation (49) of the effective stress tensor, the required conversion equation of the fourlevel elastic tensor can be immediately obtained, as shown below:
\[ {{{\bar{E}}}_{efab}}={{E}_{ijkl}}\,{{M}_{efij}}\,{{M}_{klab}} \tag{75} \]
3.3 Surface Area Damage Effect Tensor
Next one attempts to generalize equation (43) for the surface area damage variable using tensors. Let m_{ijkl} represent the new surface area damage effect tensor (which is a fourthrank tensor) that corresponds to the scalar surface area damage variable φ_{s}. Then the scalar relation in equation (43) can be written as follows using tensors with indicial notation:
\[ {{M}_{ijkl}}={{m}_{ijmn}}{{m}_{mnkl}} \tag{76} \]
Now the effective stress tensor equation (48) can be generalized as follows:
\[ {{{\bar{\sigma }}}_{ij}}={{m}_{ijmn}}{{m}_{nmkl}}{{\sigma }_{kl}} \tag{77} \]
4. Conclusion and Discussion
The nonlinear elastic constitutive equation of carbon nanotubes is used to establish a damage mechanics model to characterize the damage of carbon nanotubes. To this end, several issues were discussed. The formula is executed twiceusing scalars and tensors. The two parts containing the scalar formula and the tensor formula are kept separate. In these two parts, the transformation relationship between linear elastic variable/tensor and thirdorder elastic variable/tensor is derived.
New damage variables are also derived based on the surface area. This damage variable is very useful for studying the damage of nanomaterials such as carbon nanotubes, because surface area is critical in these types of materials. The surface area damage variable is compared with the classical damage variable based on the crosssectional area. In addition, a tensor summary of the fourthlevel surface area damage effect tensor is also proposed.
This a theoretical model based on the surface area damage effect tensor or scalar. This damage variable is very useful for studying the damage of nanomaterials such as carbon nanotubes, because surface area is critical in these types of materials. The surface area damage variable is compared with the classical damage variable based on the crosssectional area. The theoretical formulation is rigorous, however, experimental data on damage are not available at this time to calibrate the model for specific materials.
The model is valid for any configuration as this formulation is fundamental for the nonlinear elastic constitutive equation of carbon nanotubes and is used to establish a damage mechanics model to characterize the damage in carbon nanotubes. All configurations are valid for using this formulation for carbon nanotubes (CNTs) or multiwalled carbon nanotubes (MWCNTs).
Author Contributions
Both authors contributed to all sections of this work.
Funding
This research was not funded.
Competing Interests
The authors have declared that no competing interests exist.
References
 Zhang L. On the mechanics of singlewalled carbon nanotubes. J Mater Process Technol. 2009; 209: 42234228. [CrossRef]
 Pantano A, Parks DM, Boyce MC. Mechanics of deformation of single and multi wall carbon nanotubes. J Mech Phys Solids. 2004; 52: 789821. [CrossRef]
 Guo X, Wang J, Zhang H. Mechanical properties of singlewalled carbon nanotubes based on higher order Cauchy–Born rule. Int J Solids Struct. 2006; 43: 12761290. [CrossRef]
 Tanaka K, Iijima S. Carbon nanotubes and graphene. 2nd ed. Amsterdam: Elsevier; 2014.
 Li C, Chou TW. A structural mechanics approach for the analysis of carbon nanotubes. Int J Solids Struct. 2003; 40: 24872499. [CrossRef]
 Chu Y, Gautreau P, Basaran C. Parity conservation in electronphonon scattering in zigzag graphene nanoribbon. Appl Phys Lett. 2014; 105: 113112. [CrossRef]
 Chu Y, Gautreau P, Ragab T, Basaran C. An accelerated algorithm for full band electron–phonon scattering rate computation. Comput Phys Commun. 2014; 185: 33923397. [CrossRef]
 Chu Y, Ragab T, Basaran C. The size effect in mechanical properties of finitesized graphene nanoribbon. Comput Mater Sci. 2014; 81: 269274. [CrossRef]
 Chu Y, Ragab T, Basaran C. Temperature dependence of Joule heating in zigzag graphene nanoribbon. Carbon. 2015; 89: 169175. [CrossRef]
 Chu Y, Ragab T, Gautreau P, Basaran C. Mechanical properties of hydrogen edge–passivated chiral graphene nanoribbons. J Nanomech Micromech. 2015; 5: 04015001. [CrossRef]
 Kachanov LM. Time of the rupture process under creep conditions. Izy Akad Nank SSR Otd Tech Nauk. 1958; 8: 2631.
 Hao L, Ke P, June W. An anisotropic damage criterion for deformation instability and its application to forming limit analysis of metal plates. Eng Fract Mech. 1985; 21: 10311054. [CrossRef]
 Voyiadjis GZ, Kattan PI. A plasticitydamage theory for large deformation of solids—I. Theoretical formulation. Int J Eng Sci. 1992; 30: 10891108. [CrossRef]
 Voyiadjis GZ, Kattan PI. Damage mechanics. Boca Raton, Florida: CRC Press; 2005. [CrossRef]
 Voyiadjis G, Kattan P. Advances in damage mechanics: Metals and metal matrix composites with an Introduction to Fabric Tensors. 2nd ed: Elsevier; 2006.
 Voyiadjis GZ, Kattan PI. A comparative study of damage variables in continuum damage mechanics. Int J Damage Mech. 2009; 18: 315340. [CrossRef]
 Sidoroff F. Description of anisotropic damage application to elasticity. In: Physical NonLinearities in Structural Analysis. Germany: Springer; 1981. [CrossRef]
 Krajcinovic D. Damage mechanics. Amsterdam: Elsevier; 1996.
 Kattan PI, Voyiadjis GZ. A plasticitydamage theory for large deformation of solids—II. Applications to finite simple shear. Int J Eng Sci. 1993; 31: 183199. [CrossRef]
 Kattan PI, Voyiadjis GZ. Decomposition of damage tensor in continuum damage mechanics. J Eng Mech. 2001; 127: 940944. [CrossRef]
 Kattan PI, Voyiadjis GZ. Damage mechanics with finite elements: Practical applications with computer tools. Germany: Springer Science & Business Media; 2001. [CrossRef]
 Rabotonov YN. Creep Rupture. Proceedings of the XII International Congress Applied Mechanics; 1986 August 2631; Stanford, California, U.S. Germany: Springer.
 Ladeveze P, Lemaitre J. Damage effective stress in quasi unilateral conditions. Proceedings of the 16th International congress of theoretical and applied mechanics; 1984 August 1925; Lyngby, Denmark. Amsterdam: Elsevier.
 Voyiadjis GZ, Kattan PI. Mechanics of damage processes in series and in parallel: A conceptual framework. Acta Mech. 2012; 223: 18631878. [CrossRef]
 Celentano DJ, Chaboche JL. Experimental and numerical characterization of damage evolution in steels. Int J Plast. 2007; 23: 17391762. [CrossRef]
 Doghri I. Mechanics of deformable solids: Linear and nonlinear, analytical and computational aspects. Berlin: SpringerVerlag; 2000. [CrossRef]
 Hansen N, Schreyer H. A thermodynamically consistent framework for theories of elastoplasticity coupled with damage. Int J Solids Struct. 1994; 31: 359389. [CrossRef]
 Kattan PI, Voyiadjis GZ. A coupled theory of damage mechanics and finite strain elastoplasticity—I. Damage and elastic deformations. Int J Eng Sci. 1990; 28: 421435. [CrossRef]
 Ladeveze P, Poss M, Proslier L. Damage and fracture of tridirectional composites. In: Progress in science and engineering of composites. Tokyo: Japan Society for Composite Materials; 1982.
 Lubineau G. A pyramidal modeling scheme for laminatesidentification of transverse cracking. Int J Damage Mech. 2010; 19: 499518. [CrossRef]
 Lubineau G, Ladevèze P. Construction of a micromechanicsbased intralaminar mesomodel, and illustrations in ABAQUS/Standard. Comput Mater Sci. 2008; 43: 137145. [CrossRef]
 Luccioni B, Oller S. A directional damage model. Comput Methods Appl Mech Eng. 2003; 192: 11191145. [CrossRef]
 Rice JR. Inelastic constitutive relations for solids: An internalvariable theory and its application to metal plasticity. J Mech Phys Solids. 1971; 19: 433455. [CrossRef]
 Voyiadjis GZ. Degradation of elastic modulus in elastoplastic coupling with finite strains. Int J Plast. 1988; 4: 335353. [CrossRef]
 Voyiadjis GZ, Kattan PI. A coupled theory of damage mechanics and finite strain elastoplasticity—II. Damage and finite strain plasticity. Int J Eng Sci. 1990; 28: 505524. [CrossRef]
 Voyiadjis GZ, Shojaei A, Li G. A thermodynamic consistent damage and healing model for self healing materials. Int J Plast. 2011; 27: 10251044. [CrossRef]
 Voyiadjis GZ, Shojaei A, Li G, Kattan P. Continuum damagehealing mechanics with introduction to new healing variables. Int J Damage Mech. 2012; 21: 391414. [CrossRef]
 Voyiadjis GZ, Shojaei A, Li G, Kattan PI. A theory of anisotropic healing and damage mechanics of materials. Proc Math Phys Eng Sci. 2012; 468: 163183. [CrossRef]
 Darabi MK, AlRub RKA, Little DN. A continuum damage mechanics framework for modeling microdamage healing. Int J Solids Struct. 2012; 49: 492513. [CrossRef]
 Basaran C, Yan CY. A thermodynamic framework for damage mechanics of solder joints. J Electron Packag. 1998; 120: 379384. [CrossRef]
 Basaran C, Nie S. An irreversible thermodynamics theory for damage mechanics of solids. Int J Damage Mech. 2004; 13: 205223. [CrossRef]
 Onodera S, Okabe T. Analytical model for determining effective stiffness and mechanical behavior of polymer matrix composite laminates using continuum damage mechanics. Int J Damage Mech. 2020; 29: 15121542. [CrossRef]
 AbuFarsakh G, Asfa A. A unified damage model for fibrous composite laminae subject to inplane stressstate and having multi materialnonlinearity. Int J Damage Mech. 2020; 29: 13291344. [CrossRef]
 Wang S, Xu W. A coupled elastoplastic anisotropic damage model for rock materials. Int J Damage Mech. 2020; 29: 12221245. [CrossRef]
 Feng W, Qiao C, Niu S, Yang Z, Wang T. An improved nonlinear damage model of rocks considering initial damage and damage evolution. Int J Damage Mech. 2020; 29: 11171137. [CrossRef]
 Brünig M, Michalski A. Numerical analysis of damage and failure behavior of concrete. Int J Damage Mech. 2020; 29: 570590. [CrossRef]
 Baranger E. Accounting for frictional contact in an anisotropic damage model based on compliance tensorial variables, illustration on ceramic matrix composites. Int J Damage Mech. 2019; 28: 11501169. [CrossRef]
 Seupel A, Kuna M. A gradientenhanced damage model motivated by engineering approaches to ductile failure of steels. Int J Damage Mech. 2019; 28: 12611296. [CrossRef]
 Nichols JM, Abell AB. Implementing the degrading effective stiffness of masonry in a finite element model. Proceedings from the 9th North American Masonry Conference; 2003 June 23; Clemson, South Carolina.
 Nichols JM, Totoev YZ. Experimental investigation of the damage mechanics of masonry under dynamic inplane loads. Proceedings from the 8th North American Masonry Conference; 1999 June 69; Austin, Texas.
 Zhang H, Wei D. Fracture and damage behaviors of concrete in the fractal space. J Mod Phys. 2010; 1: 48. [CrossRef]