Abstract

Configuration of carbon nanotube (CNT) has been the subject of research to perform theoretical development for analyzing nanocomposites. A new theoretical solution is developed to study curved nanotube structures subjected to mechanical loadings. A curved nanotube structure is considered. A nonlocal displacement-based solution is proposed by using a displacement approach of Toroidal Elasticity based on Eringen’s theory of nonlocal continuum mechanics. The governing equations of curved nanotube structures are developed in toroidal coordinate system. The method of successive approximation is used to discretize the displacement-based governing equations and find the general solution subjected to bending moment. The numerical results show that all displacement components increase with increasing the nonlocal parameter. The present theoretical study highlights the significance of the geometry and nonlocal parameter effects on mechanical behavior of nanotube structures.

1. Introduction

Due to superior properties of CNT which offer unprecedented potential to produce multifunctional nanocomposites, many investigations have been performed on different applications of these nanostructures such as nanomechanical [1, 2] and nanobiological [3, 4] applications and optic [5] and molecular engineering [68]. CNTs appear in different shapes. Each shape has its own special properties and potential applications. CNT with curved configuration is one of their shapes. Actually, a single nanotube naturally curves during growth if no external forces exist [9]. Continuum mechanics approach has been employed by many researchers to study mechanical behavior of curved nanostructures.

The continuum mechanics approach has been used to investigate several aspects of mechanical analysis of CNT such as the static [1013], the buckling [1416], free vibration [17, 18], wave propagation [1923], and thermomechanical analysis of CNT [24].

The effect of initial stress on the vibration of an individual multiwall carbon nanotube (MWCNT) with simply supported ends was investigated based on a laminated elastic beam model [12]. Wang et al. [14] proposed nonlocal elastic beam and shell models. They applied the proposed method to study the small-scale effect on buckling analysis of CNT under compression. Based on the continuum mechanics and the molecular dynamics simulation, the investigation on the wave propagation in a single-walled carbon nanotube (SWCNT) was performed [18]. Hu et al. [19] used nonlocal elastic beam theories to study the scale effect of the CNT on the wave dispersion. They compared the results obtained by nonlocal elastic Timoshenko beam theory with the results obtained by molecular dynamics simulations. Based on thermal elasticity mechanics and nonlocal elasticity theory, a single-elastic beam model was presented to investigate the thermal vibration of SWCNT [20]. Murmu and Pradhan [20] studied the buckling of SWCNT. They considered the effect of temperature change and surrounding elastic medium.

Compared with the classical continuum theory, the nonlocal elasticity theory uses the constitutive assumptions that the stress at a point is a function of strains at all points in the continuum. Nonlocal elasticity theory was used by many researchers to analyze nanostructures advanced by Eringen [25, 26]. Peddieson et al. [27] used a version of nonlocal elasticity theory to propose a nonlocal Bernoulli-Euler beam model. Zhang et al. [28] presented a nonlocal double-elastic beam model for the free vibrations of double-walled carbon nanotubes (DWCNTs) based on theory of nonlocal elasticity and by considering the effect of small length scale. The elastic buckling analysis of nanotubes was analyzed based on Eringen’s nonlocal elasticity Timoshenko beam theories [29, 30]. The small-scale effect was considered through Eringen’s nonlocal elasticity theory while the effect of transverse shear deformation was considered through Timoshenko beam theory. Frequency equations and modal shape functions of beam structures were derived based on a nonlocal Bernoulli-Euler beam [31].

The above review shows that little work has been done to study behavior of curved nanotube structures under mechanical loadings. A new displacement-based solution is proposed to study theoretically mechanical behavior of curved nanotube structures based on Toroidal Elasticity that includes the full three-dimensional constitutive relations. Governing equations of curved nanotube are developed by using the constitutive relations of Eringen. Then, the method of successive approximation is employed to find the general solution for mechanical behavior of curved nanotube structures subjected to pure bending moment. The effects of nonlocal parameter, nanotube length, and thickness on the mechanical behavior of nanotubes are studied.

2. Nonlocal Theory (Constitutive Relations)

According to Eringen [25, 26], the stress field at a point in an elastic continuum depends not only on the strain field at the point but also on strains at all other points of the body. Eringen attributed this fact to the atomic theory of lattice dynamics and experimental observations on phonon dispersion. Therefore, the following equivalent differential form is represented based on Eringen [25, 26]:where and are the Laplacian operator and nonlocal parameter, respectively. In addition, and are a material constant and internal characteristic length, respectively.

3. Governing Equations for Curved Nanotube Structures

The nonclassical toroidal coordinate system is shown in Figure 1. The curved nanotube has a bend radius and an annular cross section bounded by radii and (see Figure 1). A general point in a constant thickness curved nanotube is represented easily by the nonclassical toroidal coordinate system , , and , where and are polar coordinates in the plane of the nanotube cross section and defines the position of the nanotube cross section. is a nondimensional radial coordinate in the toroidal coordinate, and is a reference length to be specified later. For convenience, nondimensional displacements , , and are defined asNondimensional stress components are defined aswhere and are modulus of elasticity and Poisson ratio for nanotubes. The kinematics relations for the case of small displacements are [34]where , , and represent the displacement components in the , , and directions, respectively. From both considerations of (1) and Hooke’s law, it is written thatBy considering (5), (3), and strain displacement relations (4), the nondimensional stress components may be expressed in terms of the three nondimensional displacement components and nonlocal parameter aswhereThree equilibrium equations in stress form in the toroidal coordinate are represented as [34]where is the radial cylindrical coordinate of the general point. Substituting the stress-displacement relations (6) into the equilibrium equations (8), the government Navier equations in toroidal coordinate for the curved nanotube are obtained aswhereThe three Navier equations serve as the fundamental equations for the curved nanotube. These equations are composed of three parts. The first part is independent of . The second and the third parts are the linear and nonlinear parts of . As it is impossible to find an exact solution for the Navier equations, the method of successive approximation is used to obtain an approximate solution.

4. Method of Successive Approximation

The method of successive approximation is employed to solve the Navier equations in the toroidal coordinate system. The solution for each component of the displacement is assumed to be a series in terms of a small parameter . The parameter is chosen asTherefore the nondimensional displacement components take the form aswhere the indexes of the right hand side terms indicate the order of the displacement components in the series. By substituting the series form displacement functions (12) into (3), the corresponding nondimensional stress components are derived asThe functions of Navier equations are expanded asBy substituting the expanded Navier functions (14) into Navier equations (9), matching terms in , and setting each order of to zero, the governing equations for the different orders are obtained. The details to obtain the solution are explained in Appendix A.

5. General Solution for Curved Nanotubes Subjected to Pure Bending

A direct solution is formulated to avoid the complications that arise from the use of stress function. The complementary solution for the displacement is given in the following form:The quantities , , , , , and are constants determined from zeroth-order Navier equations (13) and stress free boundary conditions on the free surfaces. The results arewhere the quantities , , , , , , , , , , , , , , , and are constants to be found from the boundary conditions. The detailed solution is presented for the curved nanotube subjected to bending moments, , as shown in Figure 2.

Solutions are proposed successively for the various orders, starting from the zeroth order. The zeroth-order solution is taken from the complementary solution. The particular part is arranged using trial functions. The displacement functions for that order are obtained by summing the complementary and particular parts. The constants in the complementary solution of each order are determined separately from the boundary conditions. The total solution is the sum of the particular and complementary solution of all the orders times the appropriate power of . Any constants remaining in the total solution are determined from end or edge boundary conditions.

5.1. The Zeroth-Order Solution

Considering the boundary conditions on the inner and outer free surfaces, the displacement components are obtained asConstant of this order is to be determined later from the boundary conditions on the ends of the curved nanotube structure.

5.2. The First-Order Solution

The first-order displacement is obtained as

5.3. The Second-Order Solution

The second-order displacement is obtained aswhere

5.4. The Third-Order Solution

The third-order solution then takes the formwhere , , , , , and are free constants. These are determined from the boundary conditions and governing equations. These constants are functions of the material property, , the ratio , and the nonlocal parameter . The coefficients appearing in (21) are defined in Appendix B.

5.5. The Fourth-Order Solution

The fourth-order solution takes the formIt should be noted that , , , , , and are free constants obtained in Appendix C by applying the boundary conditions and governing equations. Finally, substituting the solution for the various orders into (12) yields the total solution for the problem.

6. Results and Discussion

The theoretical investigations are performed for the curved single-walled carbon nanotube (SWCNT) subjected to bending moment (see Figure 2). The mechanical properties of SWCNT are given in Table 1 [32]. Furthermore, the displacement components are normalized as , where and where and are outer and inner radius of the curved SWCNT. The innermost radius and thickness of SWCNT are assumed to be 8.5 nm and 0.34 nm, respectively [33]. In addition, is the length of the curved nanotube. is defined as an aspect ratio.

Tube length, thickness, bend radius, and nonlocal parameter have effects on the nondimensional displacement components. The numerical results for curved SWCNT under bending moment () at the point in which are presented in Tables 2, 3, 4, and 5.

The effects of thickness and nonlocal parameter () on the nondimensional displacement components are investigated while assuming and as shown in Table 2. It is found from Table 2 that the displacement components increase with increasing the aspect ratio () while curved SWCNT length is kept constant. The displacement components increase with the increase of the nonlocal parameter () while the length and thickness of curved SWCNT are kept constant as shown in Table 2.

Table 3 shows the effects of thickness and nonlocal parameter (μ) on the nondimensional displacement components for and . It is observed from Table 3 that the displacement components increase with increasing the aspect ratio () while curved SWCNT length is kept constant. Observing Table 3, the displacement components increase as the nonlocal parameter () increases where length and thickness of curved SWCNT are fixed.

Comparing Table 2 with Table 3, it is clear that the nondimensional displacement components increase with increasing the curved SWCNT length.

The effects of thickness and bend radius ratio () on the nondimensional displacement components are studied for and as shown in Table 4. Observing Table 4, the displacement components decrease with the increase of bend radius ratio ().

Table 5 shows the effects of thickness and bend radius ratio () on the nondimensional displacement components for and . It is observed from Table 5 that the displacement components increase with increasing the aspect ratio () while curved SWCNT length is kept constant. Table 5 shows that the displacement components decrease as the bend radius ratio () increases where length and thickness of curved SWCNT are fixed.

7. Conclusions

We have presented a new displacement-based solution by using the nonlocal Eringen theory and Toroidal Elasticity to study curved nanotube structures under mechanical loadings. Equations of motion were developed based on Eringen’s differential constitutive equations of nonlocal elasticity for curved nanotube structures. The successive approximation method was employed to simplify the governing equations for curved nanotube structures under pure bending moment. The numerical results demonstrate the effectiveness of length, thickness, bend radius, and nonlocal parameter on all displacement components of curved nanotubes. It was observed that for curved SWCNT the values of displacement component in the longitudinal direction () were greater than those for the displacement components in the circumferential () and radial () directions.

Appendices

A. Governing Equations for the Different Orders

The equations are expressed for the zeroth order asFor the first orderFor the second orderFor the third order And for the fourth order The Navier functions of order in (A.1)–(A.5) are defined asThe displacement functions on the right hand members with negative indexes are assumed to be zero. It is noted from (A.1)–(A.5) that the zeroth-order governing equations are homogenous equations. The right hand side terms for orders greater than zero contain portions from the lower orders, and the complication of the right hand side terms increases with order. There is only a complementary solution for the zeroth-order governing equations. For orders greater than the zeroth, the solutions are composed of a complementary and a particular part. The complementary part is determined from the homogeneous part of that order’s governing equations. The particular part is derived to satisfy the right hand side terms of the governing equations. The particular part is developed by using trial displacement functions with free constants, to match the right hand side terms of the governing equations of that order.

B. The Coefficients for the First up to the Third-Order Solution

The coefficients appearing in (21) are obtained as follows:where

C. The Coefficients for the Fourth-Order Solution

The coefficients appearing in (22) are obtained as follows:where In addition, The coefficients appearing in (C.3) are given as

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.