Some features of the ASCE Shopping cart and login features of the website will be down for maintenance on Sunday, June 16th, 2024, beginning at 12:00 A.M. ET and ending at 6:00 A.M. ET. During this time if you need immediate assistance at 1-800-548-2723 or [email protected].

Free access
Technical Papers
Jun 22, 2017

Geometrically Nonlinear Static Analysis of an Embedded Multiwalled Carbon Nanotube and the van der Waals Interaction

Publication: Journal of Nanomechanics and Micromechanics
Volume 7, Issue 4

Abstract

On the basis of Reissner’s mixed variational theorem (RMVT), rather than the principle of virtual displacement (PVD), the authors presented a nonlocal Timoshenko beam theory (TBT) for the geometrically nonlinear static analysis of multiwalled carbon nanotubes (MWCNT) embedded in an elastic medium. The embedded MWCNT was subjected to mechanical loads on its outer-most surface, with combinations of free, simply supported, and clamped edge conditions. The van der Waals interaction between any pair of walls constituting the MWCNT was considered, and the interaction between the MWCNT and its surrounding medium was simulated using the Pasternak-type foundation model. In the formulation, the governing equations of a typical wall and the associated boundary conditions were derived, in which von Kármán geometrical nonlinearity was considered. Eringen’s nonlocal elasticity theory was used to account for the small-length scale effect. The deformations induced in the embedded MWCNT were obtained using the differential quadrature method and a direct iteration approach. In the numerical examples, solutions of the RMVT-based nonlocal TBT converged rapidly, and the convergent solutions of its linear counterpart closely agreed with the analytical and numerical solutions of the PVD-based nonlocal beam theories available in the literature.

Introduction

From the time carbon nanotubes (CNT) were first discovered by Iijima (1991), the development of nanoscience and nanotechnologies has proceeded at a fast rate, as has the application of CNT to various advanced industries because of their superior mechanical, electrical, thermal, and chemical properties. To assess the performance of CNT and extend their lifetimes, a number of studies carried out the structural analyses of nonembedded and embedded single- and multiwalled (SW and MW) CNT. Basically, there are four kinds of simulation methods used to conduct such analyses: (1) atomistic modeling, (2) hybrid atomistic-continuum mechanics, (3) experimental methods, and (4) continuum mechanics, with the first three approaches being especially time-consuming and expensive. However, classical continuum mechanics-based theories are not able to capture the small-length scale effect and, therefore, some modified theories have been developed, such as the strain gradient elasticity, coupled-stress elasticity, surface elasticity, and Eringen’s nonlocal elasticity theories (ENET) (Eringen 1972, 1983, 2002), with ENET being the most popular and widely used in various nonlocal continuum mechanics theories. Some comprehensive literature surveys with regard to the vibration characteristics of CNT and the various simulation methods for the analysis of these also have been carried out (Arash and Wang 2012; Brischetto 2014; Eltaher et al. 2016; Hu et al. 2012; Gibson et al. 2007; Wang et al. 2016).
Some theoretical and experimental investigations (Fu et al. 2006; Walgraef 2007) have examined the behaviors of CNT under large external loads, in which their deformations are nonlinear, rather than linear. However, there remain relatively few articles carrying out the large deflection and large amplitude vibration analyses of CNT, as compared with those examining their linear counterparts. Among the various ENET-based nonlocal beam theories used for the analysis of CNT, the nonlocal Timoshenko beam theory (TBT) commonly is used in the open literature and has been shown to produce satisfactory results as seen in applications, such as the nonlocal classical beam theory (CBT), TBT, Reddy’s beam theory (RBT), Levinson’s beam theory (LBT), generalized shear deformation beam theory (GSDBT), and sinusoidal-shear deformation beam theory (SSDBT) (Aydogdu 2008, 2009; Murmu and Pradhan 2009a, b; Reddy 2007; Reddy and Pang 2008; Thai 2012; Thai and Vo 2012; Wang et al. 2008). The literature survey in the current work thus will focus on articles related to the geometrically nonlinear analysis of SW/MW CNT that are nonembedded and embedded in an elastic medium, and different boundary conditions.
On the basis of the principle of virtual displacement (PVD), Yang et al. (2010) developed a nonlocal TBT combined with the von Kármán geometrical nonlinearity (VKGN) for the nonlinear free vibration analysis of SWCNT, with combinations of simply-supported and clamped edges, in which the effects of transverse shear deformation and rotary inertia were considered. Within the framework of the nonlocal CBT and VKGN, Ghayesh (2014) investigated the nonlinear dynamic responses of SWCNT with various boundary conditions and subjected to the primary and superharmonic excitations. A parametric study with regard to the effects of the length-to-diameter ratio, forcing amplitude and excitation frequency on the dynamic responses of SWCNT also was undertaken. Ke et al. (2009) extended the work of Yang et al. (2010) to the large amplitude vibration analysis of double-walled (DW) CNT embedded in an elastic medium. In their formulation, the nonlocal TBT-based motion equations were derived using Hamilton’s principle; the interactions between the DWCNT and its surrounding medium were considered, and the van der Waals (vdW) interaction between adjacent walls were simulated using the Winkler-type foundation model with different spring constants. The differential quadrature (DQ) method (Bert and Malik 1997; Du et al. 1994; Wu and Lee 2001) in conjunction with a direct iterative method then was used to solve the resulting nonlinear equations. This issue then was reexamined by Ansari et al. (2012), in which a two-parameter Pasternak foundation model was used to simulate the interaction between the DWCNT and its surrounding medium. Moreover, an alternative approach, the strain gradient theory, also has been used to study various structural analyses of nanoscale beams, such as linear bending analysis (Lim and Wang 2007) and wave propagation analysis (Lim et al. 2015), in which the local higher-order strain gradients were considered in a global sense.
Because almost all of the two-dimensional (2D) nonlocal beam theories used for the analysis of SW/MW CNT were derived on the basis of the PVD rather than Reissner’s mixed variational theorem (RMVT) (Reissner 1984, 1986), Wu and Lai (2015a, b) and Wu and Liou (2016) thus developed an RMVT-based nonlocal TBT for the linear bending, free vibration, and buckling analyses of SWCNT, embedded and non embedded in an elastic medium, and with different boundary conditions, in which the DQ and differential reproducing kernel methods (Wang et al. 2010) were used. In the PVD-based beam theories, the generalized displacement components (i.e., the normal displacement and rotation) were regarded as the primary variables; however, except for the generalized displacement components, the generalized resultant forces (i.e., shear force and bending moment) also were regarded as the primary variables in the RMVT-based beam theories. In Wu and Lai (2015a, b), the highest differential orders of the Euler-Lagrange equations in the RMVT-based strong and weak formulations were shown to be one-half lower than those in the PVD-based ones. The performance of the RMVT-based nonlocal TBT was shown to be superior to that of the PVD-based one, in that the convergence rate of the RMVT-based nonlocal TBT was faster than that of the PVD-based nonlocal one. Moreover, the predictions of the generalized force resultants obtained using the RMVT-based nonlocal TBT were more accurate than those obtained using the PVD-based nonlocal one. Finally, the RMVT-based nonlocal TBT was less time-consuming than the PVD-based one. Because of these benefits of the RMVT-based nonlocal TBT, and the fact that there are relatively few articles regarding the nonlinear analysis of MWCNT as compared with those related to the linear analysis of these, the authors of this paper extended the earlier work regarding the linear static analysis of embedded SWCNT to the geometrically nonlinear static analysis of embedded MWCNT with combinations of free, simply-supported, and clamped edge conditions, in which the effect of vdW interaction was considered. In the formulation, the Euler-Lagrange equations of a typical wall of the MWCNT and the associated boundary conditions were derived by using the stationary principle of Reissner’s energy functional, in which the ENET and VKGN were used. The vdW interaction forces between any pair of walls constituting the MWCNT and the Pasternak-type foundation model used to simulate the interaction between the MWCNT and its surrounding medium were taken into account. A parametric study related to the influence of some crucial effects on the nonlinear static behavior of the embedded MWCNT was conducted, such as the small-length scale effect, vdW interaction, aspect ratios, different boundary conditions, and the stiffness of the foundation.

Geometrically Nonlinear RMVT-Based Local TBT

Kinematic and Kinetic Assumptions

In this work, the authors considered a moderately thick MWCNT with a length L and the thickness of the kth-wall of the MWCNT hk (k=1-Nl), resting on a two-parameter Pasternak foundation model, and under the mechanical loads on the outer-most surface, in which Nl denotes the total number of walls constituting the CNT. The vdW interaction forces between any pair of walls were taken into account. A set of Cartesian coordinates (x, y, z) was located at the center of the MWCNT. The configuration, coordinates, and kinematics of an SWCNT are given in Fig. 1, which also can be used for a typical wall of the MWCNT by replacing the applied external load and recovery forces exerted from the foundation with the corresponding vdW interaction forces exerted from other walls. The external and internal loading system for a triple-walled (TW) CNT is shown in Fig. 2.
Fig. 1. Configuration, coordinates, and kinematics of an SWCNT [reprinted from Composite Structures, 122, Chih-Ping Wu and Wei-Wen Lai, “Reissner’s mixed variational theorem-based nonlocal Timoshenko beam theory for a single-walled carbon nanotube embedded in an elastic medium and with various boundary conditions,” pp. 390–404, Copyright (2015), with permission from Elsevier]
Fig. 2. External and internal force systems applied to a TWCNT
In the TBT, the shear deformation effect was considered to be a constant through the thickness coordinate of each wall, and the related displacement field of the kth-wall of the MWCNT is given as follows:
u1k(x,z)=uk(x)zϕk(x)
(1)
u2k(x,z)=0
(2)
u3k(x,z)=wk(x)
(3)
where uik(x,z) (i=13) = displacement components of the kth-wall of the MWCNT in the x, y, and z directions, respectively; uk(x) and wk(x) = midplane displacement components of the kth-wall in the x and z directions; and ϕk(x) = total rotation of the kth-wall in the x-z plane. The deformation of a section in the x-z plane is shown in Fig. 1.
The strain-displacement relations of the kth-wall of the MWCNT accounting for the VKGN are given by
ϵxk=uk,xzϕk,x+(1/2)(wk,x)2
(4)
γxzk=ϕk+wk,x
(5)
ϵyk=ϵzk=γyzk=γxyk=0
(6)
where ϵxk, ϵyk, ϵzk, γxzk, γyzk, and γxyk = strain components of the kth-wall of the MWCNT; and f,x=f/x, in which f=uk, wk, and ϕk.
The nonzero stress components of the kth-wall of the MWCNT in the local TBT are given by
σxk=Ek[uk,xzϕk,x+(1/2)(wk,x)2]
(7)
τxzk=kckGk(ϕk+wk,x)
(8)
where Ek and Gk = Young’s and the shear moduli of the kth-wall, respectively; kck = shear stress correction factor of the kth-wall, which accounts for the nonuniform shear stress distribution through the thickness; and the value of kck will depend on the cross section and Poisson’s ratio, and is reviewed by Cowper (1966).
The axial force (Nk), moment (Mk), and shear force (Qk) resultants of the kth-wall of the MWCNT are defined by
Nk=AkσxkdAk
(9)
Mk=AkσxkzdAk
(10)
Qk=AkτxzkdAk
(11)
where the sign conventions of Nk, Mk, and Qk can refer to Wu and Lai (2015a).
Using Eqs. (9)(11), the authors can express the stress components (σxk and τxzk) in terms of the axial force, bending moment, and shear force of the kth-wall of the MWCNT as follows:
σxk=NkAkzMkIk
(12)
τxzk=QkAk
(13)
where Ak and Ik = cross sectional area and the moment inertia of the kth-wall of the MWCNT.

vdW Interaction

In this work, a formula of the vdW interaction force derived by He et al. (2005a, b), was used. In He et al., for an MWCNT, the vdW interaction force (qij) of the ith-wall exerted from the jth-wall of a MWCNT is expressed as qij=C¯ijΔwij, in which Δwij is the relative displacement between the ith- and jth-walls, and C¯ij is the vdW interaction coefficient. According to He et al., C¯ij is given in the following form:
C¯ij=(1,001πϵσ123a4Eij131,120πϵσ69a4Eij7)Rj
(14)
where vdW interaction force is defined as the repulsive force when the value of Cij is positive, whereas it is the attractive force when the value of Cij is negative; Rj = midsurface radius of the jth-wall of the MWCNT; a = C–C bond length; and a=1.42, ϵ=2.968  meV=4.7488×1022  J, σ=0.3407  nm, and
Eijm=(Ri+Rj)m0π/2(1Kijcos2θ)m/2dθ
(15)
Kij=4RiRj(Ri+Rj)2
(16)
Because the vdW interaction coefficient C¯ij with a unit of N/(nm)3 is derived for a cylindrical shell model, it cannot be directly applied to the current nonlocal TBT. The equivalent vdW interaction coefficient Cij with a unit N/(nm)2 for a beam model, therefore, is derived as follows:
40π/2C¯ijΔwijRisin2θdθ=(πRiC¯ij)Δwij=CijΔwij
(17)
where Cij is obtained as Cij=(πRi)C¯ij; and the corresponding schematic diagram is given in Fig. 3.
Fig. 3. Distribution of the vdW interaction force of the ith-wall exerted from the jth-wall of an MWCNT and its equivalent vdW interaction force for a beam model

Stationary Principle of Reissner’s Energy Functional

The stationary principle of Reissner’s energy functional was used to derive the governing equations of the geometrically nonlinear RMVT-based local TBT, for which the generalized displacement variables (uk, wk, and ϕk) and the generalized force resultant ones (Nk, Mk, and Qk) were selected as the primary variables subject to variation, and the corresponding energy functional (ΠR) for the geometrically nonlinear RMVT-based local TBT is written in the form of
ΠR=k=1NlΠRk
(18)
where ΠRk = Reissner strain energy functional of the kth-layer of the MWCNT in the RMVT-based nonlocal TBT, and is given as follows:
ΠRk=0LAk[σxkϵxk+τxzkγxzk(σxk)2/(2Ek)(τxzk)2/(2kckGk)]dAkdx+δkNl0L(1/2)[kw(wk)2+kG(wk,x)2]dx+0L(qvdWk+δkNlq)wkdx+N¯k0uk0N¯kLukL+M¯k0ϕk0M¯kLϕkLV¯k0wk0+V¯kLwkL=0L{Nkuk,x+(Nk/2)(wk,x)2[(Nk)2/(2AkEk)]+Mkϕk,x[(Mk)2/(2EkIk)]+Qk(ϕkwk,x)[(Qk)2/(2kckGkAk)]}dx+δkNl0L(1/2)[kw(wk)2+kG(wk,x)2]dx+0L(qvdWk+δkNlq)wkdx+N¯k0uk0N¯kLukL+M¯k0ϕk0M¯kLϕkLV¯k0wk0+V¯kLwkL
(19)
where kw and kG= Winkler stiffness and shear modulus of the surrounding elastic medium, respectively; and N¯k0, N¯kL, M¯k0, M¯kL, V¯k0, and V¯kL = applied axial forces, moments, and shear forces of the kth-wall of the MWCNT at the edges, and are shown in Fig. 1; δkl = Kronecker delta symbol, in which δkl=1 when k=l, whereas δkl=0 when kl; and q and qvdWk = external load applied on the outer-most wall (k=Nl) and the vdW interaction force of the kth-wall exerted from other walls, respectively, the positive direction of which is defined to be downward. The term qvdWkδwk is given as follows:
qvdWkδwk=[l=1(lk)Nl(Cklwl)(l=1(lk)NlCkl)wk]δwk
(20)
where Ckl = vdW interaction coefficient.
Performing the first-order variation of the previously noted RMVT-based energy functional, and then using integration by parts, the authors obtain
δk=1NlΠRk=k=1Nl{0L{[uk,x+(1/2)(wk,x)2Nk/(AkEk)]δNk+(Nk,x)δuk+[φk,x(Mk/EkIk)]δMk+(Mk,x+Qk)δφk+[(φkwk,x)Qk/(kckGkAk)]δQk+[(Nkwk,x),x+Qk,x+qvdWk+δkNl(q+kwwkkGw,xxk)]δwk}dx+[(Qk0Nk0wk0,xδkNlkGwk0,xV¯k0)δwk0(QkLNkLwkL,xδkNlkGwkL,xV¯kL)δwkL(Mk0M¯k0)δφk0+(MkLM¯kL)δφkL(Nk0N¯k0)δuk0+(NkLN¯kL)δukL]}
(21)
Imposing the stationary principle of the previously noted energy functional (i.e., δk=1NlΠRk=0) leads to the governing equations of the geometrically nonlinear RMVT-based local TBT and the possible boundary conditions as follows:
Governing equations in the interior domain of the kth-wall of the MWCNT (0<x<L)
δuk:Nk,x=0
(22)
δwk:Qk,xNk,xwk,xNkwk,xx=qvdWkδkNlqδkNl(kwwkkGwk,xx)
(23)
δϕk:Mk,x+Qk=0
(24)
δNk:Nk+AkEkuk,x+(1/2)(wk,x)2=0
(25)
δMk:Mk+EkIkϕk,x=0
(26)
δQk:Qk+kckGkAk(ϕkwk,x)=0
(27)
where k=1,2,,Nl.
Possible boundary conditions at the edges of the MWCNT
Atx=0,eitheruk0=u¯k0orNk0=N¯k0
(28a)
eitherwk0=w¯k0orQk0Nk0(wk0,x)δkNlkG(wk0,x)=V¯k0
(28b)
eitherϕk0=ϕ¯k0orMk0=M¯k0
(28c)
Atx=L,eitherukL=u¯kLorNkL=N¯kL
(29a)
eitherwkL=w¯kLorQkLNkL(wkL,x)δkNlkG(wkL,x)=V¯kL
(29b)
eitherϕkL=ϕ¯kLorMkL=M¯kL
(29c)

Geometrically Nonlinear RMVT-Based Nonlocal TBT

Eringen’s Nonlocal Constitutive Relations

In the local elasticity theory, the stress components induced at a particular material point of the loaded elastic body depend only on the strain components induced at that point, whereas in the nonlocal elasticity theory, these will depend on the strain components induced at all the material points of the continuum because of the small-length scale effect. According to Eringen (1972, 1983), the nonlocal constitutive behavior of an elastic body can be written as follows:
(1μ2)σij=cijklϵkl
(30)
where μ = nonlocal parameter; μ=(e0a0)2, in which a0 is the internal characteristic length and e0 is a constant used to adjust the nonlocal continuum model by matching its results with some of the reliable results obtained by experiments or other models; e0a02.0  nm for an SWCNT (Wang and Wang 2007); cijkl = stiffness coefficient of the elastic beam-like solid; and σij and ϵkl = stress and strain components, respectively. The classical continuum mechanics should be recovered when the nonlocal parameter (μ) approaches zero.

Strong Formulation

Using Eq. (30), the constitutive equations of the kth-wall of an MWCNT for the nonlocal Timoshenko beam theory can be written as follows:
σxkμσx,xxk=Ekϵxk
(31)
τxzkμτxz,xxk=kckGkγxzk
(32)
Using Eqs. (31) and (32) in conjunction with Eqs. (9)(11), the axial force, bending moment, and shear force resultants of the kth-wall of the MWCNT in terms of the generalized displacement components can be expressed as
NkμNk,xx=AkEkuk,x+(1/2)(wk,x)2
(33)
MkμMk,xx=EkIkϕk,x
(34)
QkμQk,xx=kckGkAk(ϕkwk,x)
(35)
Eqs. (33)(35) combined with Eqs. (22)(24) constitute the set of governing equations of the geometrically nonlinear RMVT-based nonlocal TBT, and Eqs. (28a)(28c) and (29a)(29c) are the possible boundary conditions at the edges. These equations constitute the strong formulation of the geometrically nonlinear nonlocal TBT, and are a well-posed boundary value problem for the geometrically nonlinear bending analysis of an embedded MWCNT with assorted boundary conditions.

Weak Formulation

In contrast to the previously noted strong formulation, a weak formulation for the finite element (FE) analysis of the current issue can be obtained using the Galerkin approach, and is given as
k=1Nl{0L((Nk,x)δuk+[Qk,xNk,xwk,xNkwk,xx+qvdWk+δkNlq+δkNl(kwwkkGwk,xx)δwk]+(Mk,x+Qk)δϕk+{(Nk/AkEk)+μ(Nk,xx/AkEk)+[uk,x+(1/2)(wk,x)2]}δNk+[(Mk/EkIk)+μ(Mk,xx/EkIk)+ϕk,x]δMk+[(Qk/kckGkAk)+μ(Qk,xx/kckGkAk)+(ϕkwk,x)]δQk)dx+[(Qk0Nk0wk0,xδkNlkGwk0,xV¯k0)δwk0(QkLNkLwkL,xδkNlkGwkL,xV¯kL)δwkL(Mk0M¯k0)δϕk0+(MkLM¯kL)δϕkL(Nk0N¯k0)δuk0+(NkLN¯kL)δukL]}=0
(36)
Substituting the FE interpolation of each primary variable in Eq. (36), performing the integration by parts, and then integration through the element domain, the authors finally obtained the corresponding equilibrium equations for the FE method following a standard process. The relevant detailed derivation regarding the Galerkin approach can refer to Eltaher et al. (2013, 2014), and is omitted in this paper. In this work, the authors adopted a strong formulation-based numerical technique, namely the DQ method, to analyze the current issue. Because the DQ method commonly is used for the strong formulation-based numerical analysis, the detailed process can refer to Wu and Lai (2015a), and it is not repeated in this paper.

Applications

In this section, the DQ method was used to construct the weighting coefficients for the first- and higher-order derivatives of the primary variables in the spatial (x) domain, such that
Fr(x)xr|x=xi=j=1npAij(r)Fj
(37)
where r1; and Fj = (uk)j, (wk)j,   (ϕk)j, (Nk)j, (Mk)j, and (Qk)j which are the nodal values of in- and out-of-plane deflections, rotation, axial force, bending moment, and shear force components of the kth-wall of the MWCNT at x=xj.
Applying the DQ formulae [i.e., Eq. (37)] to the strong formulation of the RMVT-based nonlocal TBT for the nonlinear bending problem of an embedded MWCNT, which consists of the governing Eqs. (22)(24) and (33)(35) associated with the appropriate boundary conditions given in Eqs. (28a)(28c) and (29a)(29c), the authors obtained a set of simultaneous nonlinear algebraic equations. In these, four different boundary conditions were considered: (1) the simple-simple (SS), (2) clamped-simple (CS), (3) clamped-clamped (CC), and (4) clamped-free (CF) supports.
Applying the DQ formulae to the governing equations at the sampling nodes in the interior domain, the following is obtained:
δ(uk)i:j=1npAij(1)(Nk)j=0
(38)
δ(wk)i:j=1npAij(1)(Qk)j[j=1npAij(1)(Nk)j][j=1npAij(1)(wk)j](Nk)i[j=1npAij(2)(wk)j]=(qvdWk)iδkNlq(xi)δkNl[kw(wk)ikGj=1npAij(2)(wk)j]
(39)
δ(ϕk)i:j=1npAij(1)(Mk)j+(Qk)i=0
(40)
δ(Nk)i:(Nk)iμj=1npAij(2)(Nk)jAkEkj=1npAij(1)(uk)j(AkEk/2)[j=1npAij(1)(wk)j][j=1npAij(1)(wk)j]=0
(41)
δ(Mk)i:(Mk)iμj=1npAij(2)(Mk)jEkIkj=1npAij(1)(ϕk)j=0
(42)
δ(Qk)i:(Qk)iμj=1npAij(2)(Qk)jkckGkAk(ϕk)i+kckGkAkj=1npAij(1)(wk)j=0
(43)
where i=2,3,,(np1).
For a typical kth-wall, Eqs. (38)(43) represent 6(np2) equations in terms of 6np unknowns, whereas there are only three boundary conditions at each boundary point, and, thus, the authors needed to supplement three more governing equations at each of these. Four different cases of boundary conditions were considered and given as follows:
1.
SS supports
(uk)1=(wk)1=(Mk)1=0and(uk)np=(wk)np=(Mk)np=0
(44)
2.
CS supports
(uk)1=(wk)1=(ϕk)1=0and(uk)np=(wk)np=(Mk)np=0
(45)
3.
CC supports
(uk)1=(wk)1=(ϕk)1=0and(uk)np=(wk)np=(ϕk)np=0
(46)
4.
CF supports
(uk)1=(wk)1=(ϕk)1=0and(Nk)np=(Mk)np=[(Qk)npkGj=2npAnpj(1)(wk)j]=0
(47)
Finally, the sets of resulting equations considering the previously noted boundary conditions of Cases 1–4 are 6(np1)Nl simultaneous algebraic equations in terms of 6(np1)Nl unknowns, when k=1Nl, which represents a standard nonlinear boundary value problem. Because a direct iterative method was used to solve these nonlinear equations, they are rewritten in matrix form as follows:
{[K11(m)K12(m)K21(m)K22(m)]+[G11(m1)0G21(m1)0]}[X1(m)X2(m)]=(P0)
(48)
where superscript m=mth-iteration; coefficients G11(m1) and G21(m1) are in terms of the determined variables X1(m1) and X2(m1); G11(0)=G21(0)=0, which means the linear solutions will be used as the initial guess of the unknowns to obtain the nonlinear ones, i.e., the convergent solutions obtained from the direct iterative method with an allowable error a less than 105, in which the allowable error a is defined as a=[X(m)X(m1)][X(m)X(m1)]/[X(m)][X(m)], in which the symbol, , denotes the inner product operation for two vectors; and X(m)=[X1(m)X2(m)]T.
Once Eq. (48) is solved, the nonlinear solutions of the nodal generalized displacement and force resultant components for each layer can be obtained.

Illustrative Examples

Pressure-Loaded, Nonembedded Nanobeams

Reddy and Pang (2008) and Reddy (2007) investigated the linear static behavior of a simply-supported rectangular nanobeam subject to a uniformly-distributed load [q¯(x)=q0] by using the PVD-based nonlocal RBT and TBT, respectively, in which the area cross section of the nanobeam (A) and its associated area moment of inertia (I) were A=bh and I = bh3/12. This issue also was studied by Thai (2012) and Thai and Vo (2012), using Thai’s nonlocal RBT and the nonlocal SBT, respectively. These solutions in the linear bending case available in the literature thus are used as the benchmark solutions to validate the accuracy and convergence rate of the present RMVT-based nonlocal TBT, and a comparison between the linear and nonlinear solutions also was carried out.
For comparison purposes, the dimensionless deflection and loading intensity are defined as w¯=w/h and q¯=(q0L4)/(Eh5). The material properties of the nanobeam are the Young’s modulus (E)=1  TPa and Poisson’s ratio (ν)=0.3, and the shear correction factor (kc)=5/6.
Table 1 shows results of the convergence study for the DQ solutions of the RMVT-based nonlocal TBT for the dimensionless midspan deflection (w¯c) of a simply-supported, nanobeam subjected to a uniformly-distributed load [q(x)=q0], in which b=h=1  nm, L=10  nm, and the nonlocal parameter (μ)=0, 1, and 2 (nm)2. In the implementation of the DQ method, the total number of sampling nodes (np) was taken to be 9, 11, 13, and 21, and the node distribution was selected as the roots of an np-order Chebyshev polynomial (Wu and Lai 2015a).
Table 1. Results of the Convergence Study for the DQ Solutions of RMVT-Based Nonlocal TBT for the Dimensionless Midspan Deflection (|w¯c|) of a Simply-Supported Nanobeam without Foundation Models, and under a Uniformly Distributed Load
μ(nm)2Theoriesnpq¯=2 q¯=4 q¯=6 q¯=8
LinearNonlinearLinearNonlinearLinearNonlinearLinearNonlinear
0Present90.32030.26330.64060.41560.96090.51961.28120.6000
110.32030.26330.64060.41560.96090.51961.28120.6000
130.32030.26330.64060.41560.96090.51961.28120.6000
210.32030.26330.64060.41560.96090.51961.28120.6000
CBT0.3125N/A0.6250N/A0.9375N/A1.2500N/A
TBT0.3203N/A0.6406N/A0.9609N/A1.2812N/A
TBT0.3203N/A0.6406N/A0.9609N/A1.2812N/A
RBT0.3203N/A0.6406N/A0.9609N/A1.2812N/A
SSDBT0.3203N/A0.6406N/A0.9608N/A1.2811N/A
1Present90.35030.27680.70060.42791.05090.53021.40120.6091
110.35030.27690.70060.42801.05090.53031.40120.6093
130.35030.27690.70060.42801.05090.53031.40120.6093
210.35030.27690.70060.42801.05090.53031.40120.6093
CBT0.3425N/A0.6850N/A1.0275N/A1.3700N/A
TBT0.3503N/A0.7006N/A1.0509N/A1.4012N/A
TBT0.3509N/A0.7019N/A1.0528N/A1.4037N/A
RBT0.3509N/A0.7019N/A1.0528N/A1.4037N/A
SSDBT0.3509N/A0.7018N/A1.0527N/A1.4036N/A
2Present90.38030.28880.76060.43851.14090.53921.52120.6169
110.38030.28890.76060.43861.14090.53941.52120.6171
130.38030.28890.76060.43861.14090.53941.52120.6171
210.38030.28890.76060.43861.14090.53941.52120.6171
CBT0.3725N/A0.7450N/A1.1175N/A1.4900N/A
TBT0.3803N/A0.7606N/A1.1409N/A1.5212N/A
TBT0.3816N/A0.7631N/A1.1145N/A1.5262N/A
RBT0.3816N/A0.7631N/A1.1145N/A1.5262N/A
SSDBT0.3815N/A0.7631N/A1.1446N/A1.5261N/A
Table 1 shows that in the linear bending cases, the present DQ solutions converged rapidly, and the convergent solutions were in excellent agreement with the analytical solutions available in the literature using assorted PVD-based nonlocal theories, such as PVD-based nonlocal CBT (Reddy 2007), TBT (Reddy and Pang 2008), Reddy’s and Thai’s RBTs (Reddy 2007; Thai 2012), and SBT (Thai and Vo 2012), in which the convergent solutions were obtained when np=11. The solutions of the dimensionless deflection obtained using CBT were lower than those obtained using other nonlocal theories, which means CBT always underestimates the deflections because the effect of transverse shear deformations are discarded in the nonlocal CBT. The deviations of the solutions of the midspan deflection obtained using the nonlocal CBT and RBT were approximately 2.5%. In the linear bending cases, the variation of the dimensionless midspan deflection with the dimensionless loading intensity (q¯) was linear, although it appeared to be a higher-order polynomial function with a decreasing slope in the nonlinear bending cases. The deviation between the linear and nonlinear solutions became greater when the loading intensity was large. The ratios of nonlinear-to-linear solutions of dimensionless midspan deflection were 0.82, 0.65, 0.54, and 0.47 when the dimensionless loading intensity was q¯=2, 4, 6, and 8, respectively. The geometrically nonlinear effect stiffened the nanobeam and decreased the deflection of the nanobeam. In addition, the dimensionless midspan deflection increased when the value of the nonlocal parameter became greater, which means the small-length scale effect will soften the nanobeam.

Pressure-Loaded, Nonembedded SWCNT

In this section, the authors investigated the static behavior of zigzag (8, 0), (9, 0), and (11, 0) SWCNT with combinations of free, simply-supported, and clamped edges, and subjected to a uniformly distributed load, in which the geometric parameters and material properties of these SWCNT are given by Yang et al. (2010) as follows:
For thezigzag(8,0)SWCNT,L=4.8659  nm,2R=0.626  nm,h=0.34  nm,E=1.1556  TPa,andυ=0.19
(49a)
For thezigzag(9,0)SWCNT,L=4.8749  nm,2R=0.705  nm,h=0.34  nm,E=1.1572  TPa,andυ=0.19
(49b)
Forthezigzag(11,0)SWCNT,L=4.8857  nm,2R=0.861  nm,h=0.34  nm,E=1.1621  TPa,andυ=0.19
(49c)
The dimensionless deflection (w¯) is defined as w¯=w/D=w/(2R), in which R and D are the midsurface radius and diameter of the SWCNT considered. The area of the cross section (A) and area moment of inertia (I) are A=2πRh, and I=(π/64)[(2R+h)4(2Rh)4], respectively, and the shear correction factor kc (Cowper 1966) is taken as
kc=6(1+ν)(1+cr)2(7+6ν)(1+cr2)2+(20+12ν)cr2
(50)
in which cr=(2Rh)/(2R+h).
Table 2 shows the DQ solutions of the RMVT-based nonlocal TBT for the dimensionless midspan deflections of zigzag (8, 0), (9, 0), and (11, 0) SWCNT with SS, SC, CC, and CF edge conditions, in which the total number of sampling nodes is taken to be np=21 and the nonlocal parameters are μ=0, 1, and 2  nm2. Table 2 shows that on the basis of the same load intensity, the midspan deflections obtained using the nonlinear formulation were always less than those obtained using the linear formulation for the cases of SS, CS, and CC edges, whereas they were nearly identical to each other for the CF case. This means the effect of geometrical nonlinearity stiffened these SWCNT for the cases of SS, CS, and CC edges, whereas it was very minor for the enhancement of the gross stiffness of these SWCNT for the CF case. The deviations between the nonlinear and linear deflections tended to increase when the load intensity became greater. The results also showed that the midspan deflections of these SWCNT for different boundary conditions were CF>SS>CS>CC, in which > means larger. The midspan deflections of these SWCNT increased when the value of the nonlocal parameter became greater in the cases of SS and CS, whereas those for the SWCNT with CC edges did not have much variation in the range of μ=02  nm2, which means the small-length scale effect on the deflections of these SWCNT with different boundary conditions was SS>CS>CC. However, contrary to the former three cases, the midspan deflections of these SWCNT decreased when the value of the nonlocal parameter became greater in the cases of CF edges. Moreover, on the basis of the same values of nonlocal parameters and load intensity, and the same boundary conditions, the midspan deflections of these SWCNTs were sorted as zigzag (8, 0) > zigzag (9, 0) > zigzag (11, 0), which also meant that the midspan deflections decreased when the radius of the SWCNT became larger.
Table 2. DQ Solutions of RMVT-Based Nonlocal TBT for the Dimensionless Midspan Deflections (|w¯c|) of Zigzag (8, 0), (9, 0), and (11, 0) SWCNTs with Various Boundary Conditions, and under a Uniformly Distributed Load
Zigzag (m, n)BCsμ(nm)2q0=0.5 (N/m)q0=1.0 (N/m)q0=2.0 (N/m)q0=3.0 (N/m)
LinearNonlinearLinearNonlinearLinearNonlinearLinearNonlinear
(8, 0)SS00.12840.12510.25690.23480.51370.40250.77060.5242
10.17670.16560.35330.29230.70660.46311.06000.5800
20.22490.19970.44980.33290.89960.50051.34930.6135
CS00.05880.05860.11760.11620.23520.22490.35280.3225
10.07170.07110.14350.13880.28690.25740.43040.3540
20.08470.08320.16930.15890.33870.28090.50800.3732
CC00.03330.03330.06660.06640.13310.13200.19970.1962
10.03330.03320.06660.06620.13310.13070.19970.1921
20.03330.03320.06660.06610.13310.12950.19970.1884
CF00.43290.43290.86580.86581.73151.73152.59732.5973
10.38460.38460.76930.76931.53861.53862.30792.3079
20.33640.33640.67280.67281.34571.34572.01852.0185
(9, 0)SS00.08590.08480.17180.16380.34370.29620.51550.3992
10.11750.11370.23500.21060.47000.35400.70490.4563
20.14910.13990.29810.24730.59620.39240.89440.4919
CS00.04040.04040.08090.08040.16170.15790.24260.2306
10.04900.04880.09810.09640.19620.18460.29430.2614
20.05760.05710.11530.11150.23060.20640.34590.2835
CC00.02340.02340.04680.04670.09360.09320.14040.1390
10.02340.02340.04680.04670.09360.09270.14040.1374
20.02340.02340.04680.04660.09360.09220.14040.1359
CF00.28900.28900.57800.57801.15611.15611.73411.7341
10.25740.25740.51490.51491.02981.02981.54471.5447
20.22590.22590.45170.45170.90350.90351.35521.3552
(11, 0)SS00.04310.04300.08630.08500.17250.16360.25880.2330
10.05830.05770.11650.11240.23310.20740.34960.2837
20.07340.07200.14680.13730.29360.24150.44040.3195
CS00.02150.02150.04310.04300.08620.08550.12920.1269
10.02580.02580.05160.05140.10330.10120.15490.1483
20.03010.03000.06020.05960.12040.11580.18060.1668
CC00.01300.01300.02610.02600.05210.05200.07820.0779
10.01300.01300.02610.02600.05210.05190.07820.0776
20.01300.01300.02610.02600.05210.05180.07820.0773
CF00.14450.14450.28890.28890.57780.57780.86680.8668
10.12930.12930.25870.25870.51730.51730.77600.7760
20.11420.11420.22840.22840.45680.45680.68510.6851

Pressure-Loaded TWCNT Embedded in an Elastic Medium

In this section, the authors studied the nonlinear static behavior of a TWCNT with fully simply-supported and clamped edges, and under a uniform load, in which the vdW interaction forces of any pair of walls constituting the TWCNT were considered and the Pasternak-type foundation model was used. As noted previously, a formula derived by He et al. (2005a, b) for the vdW interaction forces between any pair of walls was modified and adopted in this work, which is given in Eqs. (14) and (17). The geometric parameters and material properties of the TWCNT were R1R2R3=(0.68  nm)(11.52), L/D3=5, E=1  TPa, and υ=0.19, and the effective thickness for each wall was h=0.34  nm. Ri and Di were the midsurface radius and diameter of the ith-wall; Di=2Ri and i=13 counted from the inner-most wall to the outer-most wall. The dimensionless load intensity (q¯), deflection (w¯), axial-spring constant (Kw), and shear modulus of the medium (KG) are defined as q¯=q0L4/E(R3)5, w¯=w/D3, Kw=kwL4/EI3, and KG=kGL2/EI3, in which I3 denoted the area moment of inertia of the third wall, which was the outer-most wall of the TWCNT.
Fig. 4 shows the DQ solutions for the deflection functions of the third wall of the TWCNT along the length direction with different boundary conditions and nonlocal parameters by using the linear and nonlinear formulations of the RMVT-based nonlocal TBT, in which Kw=KG=0, μ=0 and 2  nm2, and q¯=200. Fig. 4 shows that the maximum deflection occurred at the midspan in the cases of SS and CC. The effects of geometrical nonlinearity and small-length scale on the deflection functions of the TWCNT with different boundary conditions were SS>CC, in which > means more significant. In addition, the geometrical nonlinearity effect on the deflection functions of the TWCNT was more significant than the small-length scale effect on these.
Fig. 4. DQ solutions for the deflection functions of the third wall of an isolated TWCNT along the length direction for different boundary conditions: (a) SS; (b) CC supports
Fig. 5 shows the variations of the dimensionless midspan deflection of the third wall of the TWCNT with increasing load intensity and different nonlocal parameters for different boundary conditions, in which Kw=KG=0, μ=0 and 2  nm2, and q¯=0200. Again, it shows that in the linear cases, the variation of the dimensionless midspan deflection with the dimensionless loading intensity (q¯) was linear, whereas it appeared to be a higher-order polynomial function with a decreasing slope in the nonlinear cases. The small-length scale effect on the midspan deflections of the TWCNT in the linear cases was more significant than that in the nonlinear cases.
Fig. 5. Variations of the dimensionless midspan deflection of the third wall of an isolated TWCNT with increasing load intensity for different boundary conditions: (a) SS; (b) CC supports
Figs. 6 and 7 show the variations of the dimensionless midspan deflection of the third wall of the TWCNT with increasing load intensity and different values of Kw and KG for different boundary conditions, in which μ=1nm2, q¯=0200, KG=0, and Kw=25, 50, and 100 in Fig. 6; Kw=50 and KG=0, 5, and 10 in Fig. 7. Figs. 6 and 7 show that the midspan deflections of the embedded TWCNT tended to decrease when the values of Kw and KG became greater, which means the interaction effect between the TWCNT and its surrounding medium increased the gross stiffness of the TWCNT.
Fig. 6. Variations of the dimensionless midspan deflection of the third wall of a TWCNT resting on the Winkler-type foundation model with increasing load intensity for different boundary conditions: (a) SS; (b) CC supports
Fig. 7. Variations of the dimensionless midspan deflection of the third wall of a TWCNT resting on the Pasternak-type foundation model with increasing load intensity for different boundary conditions: (a) SS; (b) CC supports

Concluding Remarks

In this article, the authors developed the RMVT-based nonlocal TBT for the nonlinear static analysis of MWCNTs with combinations of free, simply-supported, and clamped edge conditions embedded in an elastic medium. The strong formulations of the RMVT-based nonlinear nonlocal TBT and its associated possible conditions were derived using the stationary principle of Reissner’s energy functional, in which the vdW interaction forces between any pair of walls constituting the MWCNT were considered and the interaction between the MWCNT and its surrounding medium was simulated using the Pasternak-type foundation model. The DQ method was used to discretize the generalized displacement and force resultant components, and a direct iteration method was used to solve the resulting system equations. The results showed that in the linear cases, the variation of the dimensionless midspan deflection with the dimensionless loading intensity was linear, whereas it appeared to be a higher-order polynomial function with a decreasing slope in the nonlinear cases. The geometrical nonlinear effect on the deflection functions of MWCNTs was more significant than the small-scale effect on these. The geometrical nonlinear effect on the static behavior of MWCNTs with different boundary conditions was significant for the SS, CS, and CC cases, rather than the CF case, and were SS>CS>CC, in which > means more significant.

Acknowledgments

This work was supported by the Ministry of Science and Technology of the Republic of China through Grant MOST 103-2221-E-006-064-MY3.

References

Ansari, R., Gholami, R., and Darabi, M. A. (2012). “Nonlinear free vibration of embedded double-walled carbon nanotubes with layerwise boundary conditions.” Acta Mech., 223(12), 2523–2536.
Arash, B., and Wang, Q. (2012). “A review on the application of nonlocal elastic models in modeling of carbon nanotubes and graphenes.” Comput. Mater. Sci., 51(1), 303–313.
Aydogdu, M. (2008). “Vibration of multi-walled carbon nanotubes by generalized shear deformation theory.” Int. J. Mech. Sci., 50(4) 837–844.
Aydogdu, M. (2009). “A general nonlocal beam theory: Its application to nanobeam bending, buckling and vibration.” Physica E, 41(9), 1651–1655.
Bert, C. W., and Malik, M. (1997). “Differential quadrature: A powerful new technique for analysis of composite structures.” Compos. Struct., 39(3–4), 179–189.
Brischetto, S. (2014). “A continuum elastic three-dimensional model for natural frequencies of single-walled carbon nanotubes.” Compos. Part B. Eng., 61, 222–228.
Cowper, G. (1966). “The shear coefficient in Timoshenko’s beam theory.” J. Appl. Mech., 33(2), 335–340.
Du, H., Lim, M., and Lin, R. (1994). “Application of generalized differential quadrature method to structural problems.” Int. J. Numer. Methods Eng., 37(11), 1881–1896.
Eltaher, M. A., Alshorbagy, A. E., and Mahmoud, F. F. (2013). “Vibration analysis of Euler-Bernoulli nanobeams by using finite element method.” Appl. Math. Modell., 37(7), 4787–4797.
Eltaher, M. A., Khairy, A., Sadoun, A. M., and Omar, F. A. (2014). “Static and buckling analysis of functionally graded Timoshenko nanobeams.” Appl. Math. Comput., 229, 283–295.
Eltaher, M. A., Khater, M. E., and Emam, S. A. (2016). “A review on nonlocal elastic models for bending, buckling, vibrations, and wave propagation of nanoscale beams.” Appl. Math. Modell., 40(5–6), 4109–4128.
Eringen, A. C. (1972). “Nonlocal polar elastic continua.” Int. J. Eng. Sci., 10(1), 1–16.
Eringen, A. C. (1983). “On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves.” J. Appl. Phys., 54(9), 4703–4710.
Eringen, A. C. (2002). Nonlocal continuum field theories, Springer, New York.
Fu, Y. M., Hong, J. W., and Wang, X. Q. (2006). “Analysis of nonlinear vibration for embedded carbon nanotubes.” J. Sound Vibr., 296(4–5), 746–756.
Ghayesh, M. H. (2014). “Nonlinear size-dependent behaviour of single-walled carbon nanotubes.” Appl. Phys. A, 117(3), 1393–1399.
Gibson, R. F., Ayorinde, E. O., and Wen, Y. F. (2007). “Vibrations of carbon nanotubes and their composites: A review.” Compos. Sci. Technol., 67(1), 1–28.
He, X. Q., Kitipornchai, S., and Liew, K. M. (2005a). “Buckling analysis of multi-walled carbon nanotubes: a continuum model accounting for van der Waals interaction.” J. Mech. Phys. Solids, 53(2), 303–326.
He, X. Q., Kitipornchai, S., Wang, C. M., and Liew, K. M. (2005b) “Modeling of van der Waals force for infinitesimal deformation of multi-walled carbon nanotubes treated as cylindrical shells.” Int. J. Solids Struct., 42(23), 6032–6047.
Hu, Y. G., Liew, K. M., and Wang, Q. (2012). “Modeling of vibration of carbon nanotuibes.” Procedia Eng., 31, 343–347.
Iijima, S. (1991). “Synthesis of carbon nanotubes.” Nature, 354(6348), 56–58.
Ke, L. L., Xiang, Y., Yang, J., and Kitipornchai, S. (2009). “Nonlinear free vibration of embedded double-walled carbon nanotubes based on nonlocal Timoshenko beam theory.” Comput. Mater. Sci., 47(2), 409–417.
Lim, C. W., and Wang, C. M. (2007). “Exact variational nonlocal stress modeling with asymptotic higher-order strain gradients for nanobeams.” J. Appl. Phys., 101(5), 054312.
Lim, C. W., Zhang, G., and Reddy, J. N. (2015). “A higher-order nonlocal elasticity and strain gradient theory and its applications in wave propagation.” J. Mech. Phys. Solids, 78, 298–313.
Murmu, T., and Pradhan, S. C. (2009a). “Small-scale effect on the vibration of nonuniform nanocantilever based on nonlocal elasticity theory.” Physica E, 41(8), 1451–1456.
Murmu, T., and Pradhan, S. C. (2009b). “Thermo-mechanical vibration of a single-walled carbon nanotube embedded in an elastic medium based on nonlocal elasticity theory.” Comput. Mater. Sci., 46(4), 854–859.
Reddy, J. N. (2007). “Nonlocal theories for bending, buckling and vibration of beams.” Int. J. Eng. Sci., 45(2–8), 288–307.
Reddy, J. N., and Pang, S. D. (2008). “Nonlocal continuum theories of beams for the analysis of carbon nanotubes.” J. Appl. Phys., 103(2), 023511.
Reissner, E. (1984). “On a certain mixed variational theorem and a proposed application.” Int. J. Numer. Methods Eng., 20(7), 1366–1368.
Reissner, E. (1986). “On a mixed variational theorem and on shear deformable plate theory.” Int. J. Numer. Methods Eng., 23(2), 193–198.
Thai, H. T. (2012). “A nonlocal beam theory for bending, buckling, and vibration of nanobeams.” Int. J. Eng. Sci., 52, 56–64.
Thai, H. T., and Vo, T. P. (2012). “A nonlocal sinusoidal shear deformation beam theory with application to bending, buckling, and vibration of nanobeams.” Int. J. Eng. Sci., 54, 58–66.
Walgraef, D. (2007). “On the mechanics of deformation instabilities in carbon nanotubes.” Eur. Phys. J., 146(1), 443–457.
Wang, C. M., Kitipornchai, S., Lim, C. W., and Eisenberger, M. (2008). “Beam bending solutions based on nonlocal Timoshenko beam theory.” J. Eng. Mech., 475–481.
Wang, K. F., Wang, B. L., and Kitamura, T. (2016). “A review on the application of modified continuum models in modeling and simulation of nanostructures.” Acta Mech. Sinica, 32(1), 83–100.
Wang, Q., and Wang, C. M. (2007). “The constitutive relation and small scale parameter of nonlocal continuum mechanics formodelling carbon nanotubes.” Nanotechnology, 18(7), 075702.
Wang, Y. M., Chen, S. M., and Wu, C. P. (2010). “A meshless collocation method based on the differential reproducing kernel interpolation.” Comput. Mech., 45(6), 585–606.
Wu, C. P., and Lai, W. W. (2015a). “Free vibration of an embedded single-walled carbon nanotube with various boundary conditions using the RMVT-based nonlocal Timoshenko beam theory and DQ method.” Physica E, 68, 8–21.
Wu, C. P., and Lai, W. W. (2015b). “Reissner’s mixed variational theorem-based nonlocal Timoshenko beam theory for a single-walled carbon nanotube embedded in an elastic medium and with various boundary conditions.” Compos. Struct., 122, 390–404.
Wu, C. P., and Lee, C. Y. (2001). “Differential quadrature solution for the free vibration analysis of laminated conical shells with variable stiffness.” Int. J. Mech. Sci., 43(8), 1853–1869.
Wu, C. P., and Liou, J. Y. (2016). “RMVT-based nonlocal Timoshenko beam theory for the stability analysis of embedded single-walled carbon nanotube with various boundary conditions.” Int. J. Struct. Stab. Dyn., 16(10), 1550068.
Yang, J., Ke, L. L., and Kitipornchai, S. (2010). “Nonlinear free vibration of single-walled carbon nanotubes using nonlocal Timoshenko beam theory.” Physica E, 42(5), 1727–1735.

Information & Authors

Information

Published In

Go to Journal of Nanomechanics and Micromechanics
Journal of Nanomechanics and Micromechanics
Volume 7Issue 4December 2017

History

Received: Oct 17, 2016
Accepted: Apr 21, 2017
Published online: Jun 22, 2017
Discussion open until: Nov 22, 2017
Published in print: Dec 1, 2017

Permissions

Request permissions for this article.

Authors

Affiliations

Chih-Ping Wu [email protected]
Professor, Dept. of Civil Engineering, National Cheng Kung Univ., Tainan 70101, Taiwan, ROC (corresponding author). E-mail: [email protected]
Zong-Li Hong [email protected]
Graduate Student, Dept. of Civil Engineering, National Cheng Kung Univ., Tainan 70101, Taiwan, ROC. E-mail: [email protected]
Yung-Ming Wang [email protected]
Associate Professor, Dept. of Civil Engineering, National Cheng Kung Univ., Tainan 70101, Taiwan, ROC. E-mail: [email protected]

Metrics & Citations

Metrics

Citations

Download citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by

View Options

Media

Figures

Other

Tables

Share

Share

Copy the content Link

Share with email

Email a colleague

Share