Research Papers

Distinct Element Method Modeling of Carbon Nanotube Bundles With Intertube Sliding and Dissipation

[+] Author and Article Information
Igor Ostanin

Department of Civil Engineering,
University of Minnesota,
500 Pillsbury Drive SE,
Minneapolis, MN 55455
e-mail: ostan002@umn.edu

Roberto Ballarini

Department of Civil Engineering,
University of Minnesota,
500 Pillsbury Drive SE,
Minneapolis, MN 55455
e-mail: broberto@umn.edu

Traian Dumitrică

Department of Mechanical Engineering,
University of Minnesota,
111 Church Street SE,
Minneapolis, MN 55455
e-mail: dtraian@me.umn.edu

It was assumed that an average CNT in a bundle has three slipped nearest neighbors. Critical strain εc is used to calculate stress and its terms in Eq. (13).

1Corresponding author.

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received November 27, 2013; final manuscript received January 9, 2014; accepted manuscript posted February 3, 2014; published online February 3, 2014. Editor: Yonggang Huang.

J. Appl. Mech 81(6), 061004 (Feb 03, 2014) (10 pages) Paper No: JAM-13-1487; doi: 10.1115/1.4026484 History: Received November 27, 2013; Revised January 09, 2014; Accepted February 03, 2014

The recently developed distinct element method for mesoscale modeling of carbon nanotubes is extended to account for energy dissipation and then applied to characterize the constitutive behavior of crystalline carbon nanotube bundles subjected to simple tension and to simple shear loadings. It is shown that if these structures are sufficiently long and thick, then they become representative volume elements. The predicted initial stiffness and strength of the representative volumes are in agreement with reported experimental data. The simulations demonstrate that energy dissipation plays a central role in the mechanical response and deformation kinematics of carbon nanotube bundles.

Copyright © 2014 by ASME
Your Session has timed out. Please sign back in to continue.


Thostenson, E., Li, C., and Chou, T.-W., 2005, “Nanocomposites in Context,” Compos. Sci. Technol., 65(3-4), pp. 491–516. [CrossRef]
Tahhan, M., Truong, V.-T., Spinks, G. M., and Wallace, G. G., 2003, “Carbon Nanotube and Polyaniline Composite Actuators,” Smart Mater. Struct., 12(4), pp. 626–632. [CrossRef]
Hobbie, E. K., Simien, D. O., Fagan, J. A., Huh, J. Y., Chung, J. Y., Hudson, S. D., and Obrzut, J., 2010, “Wrinkling and Strain Softening in Single-Wall Carbon Nanotube Membranes,” Phys. Rev. Lett., 104(12), p. 125505. [CrossRef] [PubMed]
Suhr, J., Kotakar, N., Keblinski, P., and Ajayan, P., 2005, “Viscoelasticity in Carbon Nanotube Composites,” Nature Mater., 4, pp. 134–137. [CrossRef]
Xu, M., Futaba, D. N., Yamada, T., Yumura, M., and Hata, K., 2010, “Carbon Nanotubes With Temperature-Invariant Viscoelasticity From −196° to 1000 °C,” Science, 330(6009), pp. 1364–1368. [CrossRef] [PubMed]
Thess, A., Lee, R., Nikolaev, P., Dai, H., Petit, P., Robert, J., Xu, C., Lee, Y. H., Kim, S. G., Rinzler, A. G., Colbert, D. T., Scuseria, G., Tomanek, D., Fischer, J. E., and Smalley, R. E., 1996, “Crystalline Ropes of Metallic Carbon Nanotubes,” Science, 273(5274), pp. 483–487. [CrossRef] [PubMed]
Yu, M.-F., Files, B. S., Arepalli, S., and Ruoff, R. S., 2000, “Tensile Loading of Ropes of Single Wall Carbon Nanotubes and Their Mechanical Properties,” Phys. Rev. Lett., 84(24), pp. 5552–5555. [CrossRef] [PubMed]
Vigolo, B., Penicaud, A., Coulon, C., Sauder, C., Pailler, R., Journet, C., Bernier, P., and Poulin, P., 2000, “Macroscopic Fibers and Ribbons of Oriented Carbon Nanotubes,” Science, 290(5495), pp. 1331–1334. [CrossRef] [PubMed]
Zhu, H. W., Xu, C. L., Wu, D. H., Wei, B. Q., Vajtai, R., and Ajayan, P. M., 2002, “Direct Synthesis of Long Single-Walled Carbon Nanotube Strands,” Science, 296(5569), pp. 884–886. [CrossRef] [PubMed]
Foroughi, J., Spinks, G. M., Wallace, G. G., Oh, J., Kozlov, M. E., Fang, S., Mirfakhrai, T., Madden, J. D. W., Shin, M. K., Kim, S. J., and Baughman, R. H., 2011, “Torsional Carbon Nanotube Artificial Muscles,” Science, 334(6055), pp. 494–497. [CrossRef] [PubMed]
Beese, A. M., Sarkar, S., Nair, A., Naraghi, M., An, Z., Moravsky, A., Loutfy, R. O., Buehler, M. J., Nguyen, S. T., and Espinosa, H. D., 2013, “Bio-Inspired Carbon Nanotube Polymer Composite Yarns With Hydrogen Bond-Mediated Lateral Interactions,” ACS Nano, 7(4), pp. 3434–3446. [CrossRef] [PubMed]
Dalton, A. B., Collins, S., Munoz, E., Razal, J. M., Ebron, V. H., Ferraris, J. P., Coleman, J. N., Kim, B. G., and Baughman, R. H., 2003, “Super-Tough Carbon-Nanotube Fibers,” Nature, 423(6941), pp. 703. [CrossRef] [PubMed]
Behabtu, N., Young, C. C., Tsentalovich, D. E., Kleinerman, O., Wang, X., Ma, A. W. K., Bengio, E. A., Waarbeek, R. F., Jong, J. J., Hoogerwerf, R. E., Fairchild, S. B., Ferguson, J. B., Maruyama, B., Kono, J., Talmon, Y., Cohen, Y., Otto, M. J., and Pasquali, M., 2013, “Strong, Light, Multifunctional Fibers of Carbon Nanotubes With Ultrahigh Conductivity,” Science339(6116), pp. 182–186. [CrossRef] [PubMed]
Dumitrică, T., Hua, M., and Yakobson, B., 2006, “Symmetry-, Time-, and Temperature-Dependent Strength of Carbon Nanotubes,” Proc. Natl. Acad. Sci. U.S.A., 103(16), pp. 6105–6109. [CrossRef] [PubMed]
Zhang, D.-B., and Dumitrică, T., 2008, “Elasticity of Ideal Single-Walled Carbon Nanotubes Via Symmetry-Adapted Tight-Binding Objective Modeling,” Appl. Phys. Lett., 93, p. 031919. [CrossRef]
Zhang, D.-B., James, R., and Dumitrică, T., 2009, “Electromechanical Characterization of Carbon Nanotubes in Torsion Via Symmetry Adapted Tight-Binding Objective Molecular Dynamics,” Phys. Rev. B, 80(11), p. 115418. [CrossRef]
Nikiforov, I., Zhang, D.-B., James, R., and Dumitrică, T., 2010, “Wavelike Rippling in Multiwalled Carbon Nanotubes Under Pure Bending,” Appl. Phys. Lett., 96(12), p. 123107. [CrossRef]
Berhan, L., Yi, Y. B., Sastrya, A. M., Munoz, E., Selvidge, M., and Baughman, R., 2004, “Mechanical Properties of Nanotube Sheets: Alterations in Joint Morphology and Achievable Moduli in Manufacturable Materials,” J. Appl. Phys., 95(8), pp. 4335–4345. [CrossRef]
Buehler, M. J., 2006, “Mesoscale Modeling of Mechanics of Carbon Nanotubes: Self-Assembly, Self-Folding, and Fracture,” J. Mater. Res., 21(11), pp. 2855–2869. [CrossRef]
Cranford, S. W., and Buehler, M. J., 2010, “In Silico Assembly and Nanomechanical Characterization of Carbon Nanotube Buckypaper,” Nanotechnology, 21(26), p. 265706. [CrossRef] [PubMed]
Hahm, M. G., Wang, H., Jung, H., Hong, S., Lee, S. G., Kim, S. R., Upmanyu, M., and JungY. J., 2012, “Bundling Dynamics Regulates the Active Mechanics and Transport in Carbon Nanotube Networks and Their Nanocomposites,” Nanoscale, 4(11), pp. 3584–3590. [CrossRef] [PubMed]
Xie, B., Liu, Y., Ding, Y., Zheng, Q., and Xu, Z., 2011, “Mechanics of Carbon Nanotube Networks: Microstructural Evolution and Optimal Design,” Soft Matter, 7(21), pp. 10039–10047. [CrossRef]
Li, Y., and Kroger, M., 2012, “A Theoretical Evaluation of the Effects of Carbon Nanotube Entanglement and Bundling on the Structural and Mechanical Properties of Buckypaper,” Carbon, 50(5), pp. 1793–1806. [CrossRef]
Li, Y., and Kroger, M., 2012, “Viscoelasticity of Carbon Nanotube Buckypaper: Zipping–Unzipping Mechanism and Entanglement Effects,” Soft Matter, 8(30), pp. 7822–7830. [CrossRef]
Zhigilei, L. V., Wei, C., and Srivastava, D., 2005, “Mesoscopic Model for Dynamic Simulations of Carbon Nanotubes,” Phys. Rev. B, 71(16), p. 165417. [CrossRef]
Volkov, A. N., and Zhigilei, L. V., 2010, “Mesoscopic Interaction Potential of Carbon Nanotubes of Arbitrary Length and Orientation,” J. Phys. Chem. C, 114(12), pp. 5513–5531. [CrossRef]
Volkov, A. N., and Zhigilei, L. V., 2010, “Structural Stability of Carbon Nanotube Films: The Role of Bending Buckling,” ACS Nano, 4(10), pp. 6187–6195. [CrossRef] [PubMed]
Jacobs, W. M., Nicholson, D. A., Zemer, H., Volkov, A. N., and Zhigilei, L. V., 2012, “Acoustic Energy Dissipation and Thermalization in Carbon Nanotubes: Atomistic Modeling and Mesoscopic Description,” Phys. Rev. B, 86(16), p. 165414. [CrossRef]
Anderson, T., Akatyeva, E., Nikiforov, I., Potyondy, D., Ballarini, R., and Dumitrică, T., 2010, “Toward Distinct Element Method Simulations of Carbon Nanotube Systems,” ASME J. Nanotechnol. Eng. Med., 1(4), p. 041009. [CrossRef]
Ostanin, I., Ballarini, R., Potyondy, D., and Dumitrică, T., 2013, “A Distinct Element Method for Large Scale Simulations of Carbon Nanotube Assemblies,” J. Mech. Phys. Solids, 61(3), pp. 762–782. [CrossRef]
Filleter, T., Yockel, S., Naraghi, M., Paci, J. T., Compton, O. C., Mayes, M. L., Nguyen, S-B. T.,Schatz, G. C., and Espinosa, H. D., 2012, “Experimental-Computational Study of Shear Interactions Within Double-Walled Carbon Nanotube Bundles,” Nano Lett., 12(2), pp. 732–742. [CrossRef] [PubMed]
Zhang, Q., Lu, Y. C., Du, F., Dai, L., Baur, J., and Foster, D. C., 2010, “Viscoelastic Creep of Vertically Aligned Carbon Nanotubes,” J. Phys. D: Appl. Phys., 43(31), p. 315401. [CrossRef]
Carlson, A., and Dumitrică, T., 2007, “Extended Tight-Binding Potential for Modeling Intertube Interactions in Carbon Nanotubes,” Nanotechnology, 18(6), p. 065706. [CrossRef]
Johnson, S. M., Williams, J. R., and Cook, B. K., 2007, “Quaternion-Based Rigid Body Rotation Integration Algorithms for Use in Particle Methods,” Int. J. Numer. Meth. Eng., 74(8), pp. 1303–1313. [CrossRef]
Itasca Consulting Group Inc., 2008, “PFC3D (Particle Flow Code in 3 Dimensions),” Version 4.0, Itasca Consulting Group Inc., Minneapolis, MN.
Tang, Y., Ballarini, R., Buehler, M. J., and Eppell, S. J., 2010, “Deformation Micromechanisms of Collagen Fibrils Under Uniaxial Tension,” J. R. Soc. Interface, 7(46), pp. 839–850. [CrossRef] [PubMed]
Li, C., Liu, Y., Yao, X., Ito, M., Noguchi, T., and Zheng, Q., 2010, “Interfacial Shear Strengths Between Carbon Nanotubes,” Nanotechnology, 21(11), p. 115704. [CrossRef] [PubMed]
Moore, D. F., Ding, Y., and Wang, Z. L., 2004, “Crystal Orientation-Ordered ZnS Nanowire Bundles,” J. Am. Chem. Soc., 126(44), pp. 14372–14373. [CrossRef] [PubMed]
Ma, J., Liu, X., Cao, X., Feng, S., and Fleet, E. M., 2006, “Bundle of Nanobelts Up to 4 cm in Length: One-Step Synthesis and Preparation of Titanium Trisulfide (TiS3) Nanomaterials,” Eur. J. Inorganic Chem., 2006(3), pp. 519–522. [CrossRef]


Grahic Jump Location
Fig. 1

(a) Coarse graining of a CNT into a chain of cylindrical segments, representing inertial properties of a CNT. Segments are linked with parallel bond interfaces, representing elastic properties of a CNT surface, arising from covalent bonds between C atoms within the CNT surface. Segments of neighboring CNTs interact via mesoscopic vdW contact model, acting in parallel with viscous forces that damp relative translational motion of CNT segments. (b) Illustration of definition of generalized coordinates r, θ, φ for two interacting cylindrical segments of neighboring CNTs. n1 and n2 are axial directions of first and second cylinder.

Grahic Jump Location
Fig. 2

(a) Schematics of a crystalline CNT bundle (N = 2). (b) Geometry of the prerelaxed bundle (N = 4, M = 4). Color legend gives x-component of displacements developing during the relaxation. (c) Changes in different terms of CNT bundle energy during prerelaxation. The change in strain energy ΔUstr, vdW adhesion energy ΔUvdw, total energy ΔUtot, and kinetic energy Ukin are presented. (d) Gage and grip regions of a CNT bundle specimen.

Grahic Jump Location
Fig. 3

(a), (b) SSCs of bundles of different length, containing 37 tubes on a cross section (N = 4), for LCNT = 0.068μm (specimens 1–9, (a)) and LCNT = 0.136μm (specimens 10–18, (b)). (c) Uniaxial tensile strength of a bundle σuts as a function of the length factor M for LCNT = 0.068μm and LCNT = 0.136μm. (d), (e) SSCs of bundles of different thickness, with length factor M = 4, and LCNT = 0.068μm (specimens 19–24, (d)) and LCNT = 0.136μm (specimens 25–30, (e)). (f) Uniaxial tensile strength of a bundle σuts as a function of the thickness factor N for LCNT = 0.068μm and LCNT = 0.136μm. The slope of reference lines given in figures (a), (b), (d), (e) corresponds to Young's modulus E0 estimated for hexagonal arrangement of stretched noninteracting CNTs.

Grahic Jump Location
Fig. 4

(a) Work of external force Aext, as compared to changes in elastic strain energy ΔUstr, vdW adhesion energy ΔUvdw, and dissipated energy ΔQ during the test. (b) Terms of decomposition (Eq. (13)) during the simulation. Total stress response σtot is calculated as a sum of derivatives of traced energy terms (dashed line) and directly from force balance (solid line). (c) Development of a localized deformation in a CNT bundle. Visualization of CNT bundle geometry and magnitude of a slip vector (on the surface of a bundle and on a horizontal axial cross section of a bundle).

Grahic Jump Location
Fig. 5

Low damping (α = 0.2,β = 0.0) mechanical tests on CNT bundles. (a) Stress-strain curves for specimen 12 and specimen 16 in Table 2 indicate the absence of an RVE. (b) Magnitude of slip vector, visualized on a horizontal axial cross-section, indicates immediate localization of the deformation and brittle fracture of a specimen. (c) Magnitude of a slip vector, averaged over thin slices of a bundle (specimen 16) along the length for few different values of viscous and local damping.

Grahic Jump Location
Fig. 6

Large deformation of a CNT bundle. (a) Work of external force Aext, as compared to changes in elastic strain energy ΔUstr, vdW adhesion energy ΔUvdw, and dissipated energy ΔQ during the test. (b) Terms of decomposition (Eq. (13)) during the simulation. (c) CNT bundle geometry and magnitude of a slip vector on the surface of a bundle and on a horizontal axial cross section of a bundle.

Grahic Jump Location
Fig. 7

Shear test on a close-packed CNT assembly. (a) Problem geometry and boundary conditions. (b) Definition of the cross-sectional area of a specimen for given values of B and H. (c) Displacement field (x component of displacement) observed in a simple shear test. (d) External work as compared to potential energy terms and dissipation during the test. (e) Decomposition of stress into vdW adhesive, elastic, and dissipative terms. (f) The influence of specimen length on its shear modulus. Red circles are simulation results, and blue crosses are theoretical predictions for the shear modulus due to surface tension. On the inset—the effect of cross section size (B and H) on the shear modulus.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In