Journal of Energy and Power Technology (JEPT) is an international peer-reviewed Open Access journal published quarterly online by LIDSEN Publishing Inc. This periodical is dedicated to providing a unique, peer-reviewed, multi-disciplinary platform for researchers, scientists and engineers in academia, research institutions, government agencies and industry. The journal is also of interest to technology developers, planners, policy makers and technical, economic and policy advisers to present their research results and findings.

Journal of Energy and Power Technology focuses on all aspects of energy and power. It publishes not only original research and review articles, but also various other types of articles from experts in these fields, such as Communication, Opinion, Comment, Conference Report, Technical Note, Book Review, and more, to promote intuitive understanding of the state-of-the-art and technology trends.

Main research areas include (but are not limited to):
Renewable energies (e.g. geothermal, solar, wind, hydro, tidal, wave, biomass) and grid connection impact
Energy harvesting devices
Energy storage
Hybrid/combined/integrated energy systems for multi-generation
Hydrogen energy 
Fuel cells
Nuclear energy
Energy economics and finance
Energy policy
Energy and environment
Energy conversion, conservation and management
Smart energy system

Power Generation - Conventional and Renewable
Power System Management
Power Transmission and Distribution
Smart Grid Technologies
Micro- and nano-energy systems and technologies
Power electronic
Biofuels and alternatives
High voltage and pulse power
Organic and inorganic photovoltaics
Batteries and supercapacitors

Publication Speed (median values for papers published in 2023): Submission to First Decision: 5.1 weeks; Submission to Acceptance: 11.6 weeks; Acceptance to Publication: 7 days (1-2 days of FREE language polishing included)

Current Issue: 2024  Archive: 2023 2022 2021 2020 2019
Open Access Original Research

Damage Mechanics of Carbon Nanotubes

George Z. Voyiadjis 1,*, Peter I. Kattan 2

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

  2. 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): 021; 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 cross-sectional area. In the presentation, both the one-dimensional case (scalars) and the three-dimensional 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 cross-sectional 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].

One-dimenstional 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.

Three-dimenstional 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. One-Dimensional Formulation (Scalars)

The elastic stress-strain formula for carbon nanotubes is given by the following formula for the one-dimensional 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 second-order elastic stiffness and third-order 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 cross-sectional 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 E ¯ = E f ( ϕ ) and D ¯ = D g ( ϕ ) , where f and g are unspecified functions. The authors will show the explicit expressions of the scalar functions where E ¯ and D ¯ are the second-order effective lossless elastic stiffness and the third-order 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 used-one 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 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 C1 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 C2 is a constant of integration to be determined. Substituting equations (9) and (10) into equation (8), and utilizing the fact that c1 = 0 and c2 = 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 re-written 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 third-order 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 third-order 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 and m defined as follows in terms of elastic stiffness reduction:

\[\ell = \frac{\bar{E} -E}{E} \tag{20}\]

\[ m=\frac{\bar{D}-D}{D} \tag{21}\]

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 and m vanish for the virgin/undamaged material. However, when the value of m reaches 1, then the value of 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 third-order 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.

Click to view original image

Figure 1 Relationship between the two Damage Variables l and m.

One can get further insight into the relationship between the two damage variables 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 two-thirds of the reduction in the third-order 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 one-half of the reduction in the third-order 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 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 cross-sectional 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 cross-sectional area, S is the surface area, P is the perimeter, t is the thickness, and α is a shape constant (e.g. α = 1 16 for a square and α = 1 4 π for a circle).

Starting with the definition of the classical cross-sectional 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 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 re-written 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 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 cross-sectional area damage variable φ and the surface area damage variable φs. This relationship is re-written 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 cross-sectional 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.

Click to view original image

Figure 2 Relationship between the Cross-Sectional 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. Three-Dimensional Formulation (Tensors)

Finally, the whole process is extended to (using tensors) the three-dimensional 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 second-rank stress and second-rank strain tensors, respectively. In equation (44), Eijkl and Dijklmn are the fourth-rank second-order elasticity tensor and the sixth-rank third-order 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 E 1111 and D 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 fourth-rank second-order elasticity tensor Eijkl and the sixth-rank third-order elasticity tensor Dijklmn 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 σ ¯ i j is the effective stress tensor (defined in the fictitious/undamaged configuration) and Mijkl is the fourth-rank 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 σ ¯ i j 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} \]

Pre-multiplying equation (50) by M ( p q i j ) ( 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 ( p q i j ) ( 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 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 ( ϵ i j ϵ k l ) = 2 ϵ k l d ϵ i j and d ( ϵ i j ϵ k l ϵ k m n ) = 3 ϵ k l ϵ k m n d ϵ i j , 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} \]

Post-multiplying equation (59) by ϵ ¯ ( l p ) ( 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 Ill is related to the second-rank identity tensor I i l δ i j .

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 Ill = 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} \]

Pre-multiplying equation (65) by σ ( p q ) ( 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 sixth-level third-order 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 fourth-level 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} \]

Pre-multiplying equation (71) by Mefil 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 Iefab represents the fourth-rank 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 four-level 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 mijkl represent the new surface area damage effect tensor (which is a fourth-rank 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 twice-using 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 third-order 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 cross-sectional area. In addition, a tensor summary of the fourth-level 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 cross-sectional 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

  1. Zhang L. On the mechanics of single-walled carbon nanotubes. J Mater Process Technol. 2009; 209: 4223-4228. [CrossRef]
  2. Pantano A, Parks DM, Boyce MC. Mechanics of deformation of single- and multi- wall carbon nanotubes. J Mech Phys Solids. 2004; 52: 789-821. [CrossRef]
  3. Guo X, Wang J, Zhang H. Mechanical properties of single-walled carbon nanotubes based on higher order Cauchy–Born rule. Int J Solids Struct. 2006; 43: 1276-1290. [CrossRef]
  4. Tanaka K, Iijima S. Carbon nanotubes and graphene. 2nd ed. Amsterdam: Elsevier; 2014.
  5. Li C, Chou TW. A structural mechanics approach for the analysis of carbon nanotubes. Int J Solids Struct. 2003; 40: 2487-2499. [CrossRef]
  6. Chu Y, Gautreau P, Basaran C. Parity conservation in electron-phonon scattering in zigzag graphene nanoribbon. Appl Phys Lett. 2014; 105: 113112. [CrossRef]
  7. Chu Y, Gautreau P, Ragab T, Basaran C. An accelerated algorithm for full band electron–phonon scattering rate computation. Comput Phys Commun. 2014; 185: 3392-3397. [CrossRef]
  8. Chu Y, Ragab T, Basaran C. The size effect in mechanical properties of finite-sized graphene nanoribbon. Comput Mater Sci. 2014; 81: 269-274. [CrossRef]
  9. Chu Y, Ragab T, Basaran C. Temperature dependence of Joule heating in zigzag graphene nanoribbon. Carbon. 2015; 89: 169-175. [CrossRef]
  10. Chu Y, Ragab T, Gautreau P, Basaran C. Mechanical properties of hydrogen edge–passivated chiral graphene nanoribbons. J Nanomech Micromech. 2015; 5: 04015001. [CrossRef]
  11. Kachanov LM. Time of the rupture process under creep conditions. Izy Akad Nank SSR Otd Tech Nauk. 1958; 8: 26-31.
  12. 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: 1031-1054. [CrossRef]
  13. Voyiadjis GZ, Kattan PI. A plasticity-damage theory for large deformation of solids—I. Theoretical formulation. Int J Eng Sci. 1992; 30: 1089-1108. [CrossRef]
  14. Voyiadjis GZ, Kattan PI. Damage mechanics. Boca Raton, Florida: CRC Press; 2005. [CrossRef]
  15. Voyiadjis G, Kattan P. Advances in damage mechanics: Metals and metal matrix composites with an Introduction to Fabric Tensors. 2nd ed: Elsevier; 2006.
  16. Voyiadjis GZ, Kattan PI. A comparative study of damage variables in continuum damage mechanics. Int J Damage Mech. 2009; 18: 315-340. [CrossRef]
  17. Sidoroff F. Description of anisotropic damage application to elasticity. In: Physical Non-Linearities in Structural Analysis. Germany: Springer; 1981. [CrossRef]
  18. Krajcinovic D. Damage mechanics. Amsterdam: Elsevier; 1996.
  19. Kattan PI, Voyiadjis GZ. A plasticity-damage theory for large deformation of solids—II. Applications to finite simple shear. Int J Eng Sci. 1993; 31: 183-199. [CrossRef]
  20. Kattan PI, Voyiadjis GZ. Decomposition of damage tensor in continuum damage mechanics. J Eng Mech. 2001; 127: 940-944. [CrossRef]
  21. Kattan PI, Voyiadjis GZ. Damage mechanics with finite elements: Practical applications with computer tools. Germany: Springer Science & Business Media; 2001. [CrossRef]
  22. Rabotonov YN. Creep Rupture. Proceedings of the XII International Congress Applied Mechanics; 1986 August 26-31; Stanford, California, U.S. Germany: Springer.
  23. Ladeveze P, Lemaitre J. Damage effective stress in quasi unilateral conditions. Proceedings of the 16th International congress of theoretical and applied mechanics; 1984 August 19-25; Lyngby, Denmark. Amsterdam: Elsevier.
  24. Voyiadjis GZ, Kattan PI. Mechanics of damage processes in series and in parallel: A conceptual framework. Acta Mech. 2012; 223: 1863-1878. [CrossRef]
  25. Celentano DJ, Chaboche J-L. Experimental and numerical characterization of damage evolution in steels. Int J Plast. 2007; 23: 1739-1762. [CrossRef]
  26. Doghri I. Mechanics of deformable solids: Linear and nonlinear, analytical and computational aspects. Berlin: Springer-Verlag; 2000. [CrossRef]
  27. Hansen N, Schreyer H. A thermodynamically consistent framework for theories of elastoplasticity coupled with damage. Int J Solids Struct. 1994; 31: 359-389. [CrossRef]
  28. Kattan PI, Voyiadjis GZ. A coupled theory of damage mechanics and finite strain elasto-plasticity—I. Damage and elastic deformations. Int J Eng Sci. 1990; 28: 421-435. [CrossRef]
  29. 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.
  30. Lubineau G. A pyramidal modeling scheme for laminates-identification of transverse cracking. Int J Damage Mech. 2010; 19: 499-518. [CrossRef]
  31. Lubineau G, Ladevèze P. Construction of a micromechanics-based intralaminar mesomodel, and illustrations in ABAQUS/Standard. Comput Mater Sci. 2008; 43: 137-145. [CrossRef]
  32. Luccioni B, Oller S. A directional damage model. Comput Methods Appl Mech Eng. 2003; 192: 1119-1145. [CrossRef]
  33. Rice JR. Inelastic constitutive relations for solids: An internal-variable theory and its application to metal plasticity. J Mech Phys Solids. 1971; 19: 433-455. [CrossRef]
  34. Voyiadjis GZ. Degradation of elastic modulus in elastoplastic coupling with finite strains. Int J Plast. 1988; 4: 335-353. [CrossRef]
  35. Voyiadjis GZ, Kattan PI. A coupled theory of damage mechanics and finite strain elasto-plasticity—II. Damage and finite strain plasticity. Int J Eng Sci. 1990; 28: 505-524. [CrossRef]
  36. Voyiadjis GZ, Shojaei A, Li G. A thermodynamic consistent damage and healing model for self healing materials. Int J Plast. 2011; 27: 1025-1044. [CrossRef]
  37. Voyiadjis GZ, Shojaei A, Li G, Kattan P. Continuum damage-healing mechanics with introduction to new healing variables. Int J Damage Mech. 2012; 21: 391-414. [CrossRef]
  38. 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: 163-183. [CrossRef]
  39. Darabi MK, Al-Rub RKA, Little DN. A continuum damage mechanics framework for modeling micro-damage healing. Int J Solids Struct. 2012; 49: 492-513. [CrossRef]
  40. Basaran C, Yan CY. A thermodynamic framework for damage mechanics of solder joints. J Electron Packag. 1998; 120: 379-384. [CrossRef]
  41. Basaran C, Nie S. An irreversible thermodynamics theory for damage mechanics of solids. Int J Damage Mech. 2004; 13: 205-223. [CrossRef]
  42. 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: 1512-1542. [CrossRef]
  43. Abu-Farsakh G, Asfa A. A unified damage model for fibrous composite laminae subject to in-plane stress-state and having multi material-nonlinearity. Int J Damage Mech. 2020; 29: 1329-1344. [CrossRef]
  44. Wang S, Xu W. A coupled elastoplastic anisotropic damage model for rock materials. Int J Damage Mech. 2020; 29: 1222-1245. [CrossRef]
  45. 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: 1117-1137. [CrossRef]
  46. Brünig M, Michalski A. Numerical analysis of damage and failure behavior of concrete. Int J Damage Mech. 2020; 29: 570-590. [CrossRef]
  47. 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: 1150-1169. [CrossRef]
  48. Seupel A, Kuna M. A gradient-enhanced damage model motivated by engineering approaches to ductile failure of steels. Int J Damage Mech. 2019; 28: 1261-1296. [CrossRef]
  49. 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.
  50. Nichols JM, Totoev YZ. Experimental investigation of the damage mechanics of masonry under dynamic in-plane loads. Proceedings from the 8th North American Masonry Conference; 1999 June 6-9; Austin, Texas.
  51. Zhang H, Wei D. Fracture and damage behaviors of concrete in the fractal space. J Mod Phys. 2010; 1: 48. [CrossRef]
Newsletter
Download PDF Download Citation
0 0

TOP