A Variational Formulation for Coupled Single-Walled Carbon Nanotubes Undergoing Vibrations in the Presence of an Axial Magnetic Field

A Variational Formulation for Coupled Single-Walled Carbon Nanotubes Undergoing Vibrations in the Presence of an Axial Magnetic Field

Sarp Adali

Discipline of Mechanical Engineering, University of KwaZulu-Natal, Durban 4041, South Africa

Corresponding Author Email: 
adali@ukzn.ac.za
Page: 
1180-1188
|
DOI: 
https://doi.org/10.18280/mmep.100410
Received: 
12 November 2022
|
Revised: 
18 January 2023
|
Accepted: 
25 January 2023
|
Available online: 
30 August 2023
| Citation

© 2023 IIETA. This article is published by IIETA and is licensed under the CC BY 4.0 license (http://creativecommons.org/licenses/by/4.0/).

OPEN ACCESS

Abstract: 

A variational formulation is presented for a system consisting of two coupled single-walled carbon nanotubes subjected to forced vibrations, longitudinal magnetic field, and axial compression, based on the nonlocal elasticity theory. The variational principle for the double nanotube system is derived, followed by the application of Hamilton's principle to express kinetic and potential energies. Subsequently, the time-independent scenario is investigated, and governing equations for the freely vibrating system are provided. The variational formulation for this case is established, and expressions for Rayleigh quotients concerning the vibration frequency and buckling load are derived. The Rayleigh quotient for the frequency demonstrates that the magnetic field increases the vibration frequency of the coupled nanotube system. Nonlocal effects appear in both the numerator and the denominator of the Rayleigh quotient, influencing the frequency increase or decrease depending on the relative values of various problem parameters. In contrast, the magnetic field reduces the buckling load, as evidenced by its negative contribution to the numerator of the Rayleigh quotient for buckling. The effect of the nonlocal parameter on buckling, however, cannot be inferred directly from the Rayleigh quotient. In this study, involving a system of two coupled partial differential equations, it is crucial to derive variationally consistent boundary conditions. Utilizing the formulated variational principle, variationally consistent natural boundary conditions are established in terms of moment and shear force expressions. It is revealed that the Pasternak interlayer between the nanotubes results in coupled boundary conditions when a shear force and/or a moment is specified at the boundaries of the nanotube system.

Keywords: 

carbon nanotubes, magnetic field, Winkler-Pasternak interlayer, variational formulation, variationally consistent boundary conditions, vibrations

1. Introduction

Magnetic fields play a significant role in various nanotechnology applications, particularly in nano and micro electromechanical systems and nanosensors. Carbon nanotubes (CNTs) exposed to magnetic fields exhibit unique features, which can be harnessed for numerous nanotechnology applications as reported by Kibalchenko et al. [1]. The influence of magnetic fields on the vibration frequencies of CNTs was discussed by Zhang et al. [2], suggesting that magnetic fields can be utilized to control the dynamic characteristics of CNTs. Moreover, magnetic fields affect wave propagation in nanotubes with implications for fluid conveyance in CNTs, as noted by Arani et al. [3]. Alazwari et al. [4] observed that exposure to a magnetic field induces a hardening effect, increasing system stiffness and resulting in higher frequencies. This effect can be employed for resonance control in nanosensors, as reported by Kiani [5].

In the realm of biomedicine, magnetic fields offer advantages such as the development of non-invasive and harmless medical instruments like MRIs, as mentioned by Samadishadlou et al. [6]. Another medical application involves the use of CNTs with ferromagnetic nanoparticles for hyperthermia studies when exposed to a magnetic field, as studied by Raniszewski et al. [7]. Several studies by Bellucci et al. [8], Klinovaja et al. [9], Kibalchenko et al. [10], Kiani [11], and Fedorov et al. [12] have investigated the effect of magnetic fields on CNTs and CNT devices, providing valuable insights for the development of various nanodevices. For instance, Popov et al. [13] devised a magnetic field sensor, while Mandal et al. [14] and Pal et al. [15] explored magnetic nanomotors.

Recent investigations have delved into the magnetic properties of nanotubes and the behavior of nanotubes in a magnetic field. Bellucci et al. [16] examined the impact of magnetic fields on the transport properties of nanotubes, while Kibalchenko et al. [17] investigated their electronic properties. Furthermore, Li et al. [18] and Wang et al. [19] studied the influence of magnetic fields on the dynamics of multiwalled nanotubes. The mechanics of CNTs in a magnetic field have been extensively researched, with studies by Wang et al. [20, 21], Kiani [22], and Narendar et al. [23] focusing on the effect of a longitudinal magnetic field on wave propagation in single-walled carbon nanotubes (SWCTs).

More recent studies on wave propagation in SWCTs exposed to magnetic fields include works by Arani et al. [24, 25] involving fluid-conveying single and double-walled nanotubes, and Li et al. [26], which considered wave propagation in viscoelastic single-walled nanotubes while incorporating the surface effect using strain gradient theory. Zhen and Zhou [27] investigated wave propagation in fluid-conveying nanotubes under the influence of magnetic fields, considering thermal and surface effects. Vibrations of carbon nanotubes exposed to magnetic fields were studied for double-walled nanotubes by Kiani [28, 29] and Murmu et al. [30], and for double single-walled carbon nanotube systems by Murmu et al. [31] and Arani et al. [32]. Alazwari et al. [4] examined vibrations of functionally graded nanobeams in a magnetic field, while Kiani [33] investigated vibration and buckling of CNTs in a three-dimensional magnetic field.

Jena et al. [34] and Esen et al. [35] studied the effect of magnetic fields on vibrations and buckling of nanobeams. Nonlinear vibrations of nanobeams and nanotubes in a magnetic field were investigated by Chang [36], Ebrahimi and Hosseini [37], and Anh and Hieu [38]. Jalaei et al. [39] examined the dynamic stability of a Timoshenko nanobeam in a magnetic field, while Bahaadini et al. [40] focused on viscoelastic nanotubes conveying magnetic nanoflow in a magnetic field. These studies primarily explored wave propagation and linear/nonlinear vibrations of single-walled nanotubes and nanobeams under specific boundary conditions.

Complex mechanical systems formed by multiple carbon nanotubes, particularly those involving double single-walled nanotubes (DSWNTs), have garnered significant attention in recent years due to their distinct mechanical properties from double-walled nanotubes. A plethora of research has focused on studying the vibrational behavior of DSWNTs and double nanobeams connected by an elastic layer in the presence of a magnetic field, as evidenced by the works of Kiani [28, 29], Murmu et al. [30, 31], Arani et al. [32], Nasirshoaibi et al. [41], and Stamenkovic et al. [42]. Ebrahimi and Dabbagh [43] extended these findings by investigating the effect of magnetic fields on the propagation of acoustical waves in rotary double-nanobeam systems.

The current study aims to develop a variational formulation for a double single-walled carbon nanotube (SWCNT) system subjected to forced vibrations, an axial magnetic field, and axial compression, based on nonlocal elasticity theory. To date, a rigorous derivation of variational principles and associated boundary conditions for such double CNT systems exposed to magnetic fields remains unavailable. By providing a variational framework, the present study paves the way for future investigations into the vibration and buckling behavior of CNTs in magnetic fields.

Variational formulations have been previously explored in numerous cases involving carbon nanotubes and nanobeams, including those by Yang and Lim [44], Baretta et al. [45], Barretta et al. [46], and Chalamel [47]. Ike [48] presented a formulation based on the system's energy for flexural torsional buckling of open cross-section columns, and proposed a variational formulation for a Timoshenko beam undergoing bending in a subsequent study [49].

Deriving variationally consistent boundary conditions is crucial for complex mechanical systems. Several studies have focused on this topic, such as those by Shi and Voyiadjis [50] for shear deformable beams, Robinson and Adali [51] for carbon nanotubes subjected to uniformly and triangularly distributed axial loads, and Barretta et al. [52] for Timoshenko nanobeams based on nonlocal strain gradient theory. Further examples include research by Yu et al. [53] on size-dependent beams, Xu and Zheng [54] on Timoshenko nanobeams, and Pinnola et al. [55] on nonlocal gradient elastic beams. Barretta et al. [56] obtained variationally consistent boundary conditions for nonlocal strain gradient Timoshenko nanobeams.

Adali [57, 58] derived variational principles and variationally consistent boundary conditions for multi-walled carbon nanotubes experiencing buckling and vibrations, respectively. Xu and Deng [59] extended these findings to natural and geometric boundary conditions for multi-walled nanotubes subject to buckling and vibrations based on strain gradient theory. Additional research on this topic includes work by Kucuk et al. [60] on vibrating multi-walled nanotubes modeled as Timoshenko beams and Adali [61] on buckling. Adali [62] also formulated variational principles and variationally consistent boundary conditions for a double Rayleigh beam system based on local elasticity theory.

In this study, we present a variational principle for a double nanotube system subjected to a magnetic field, based on nonlocal elastic theory. By identifying the kinetic and potential energies in the variational expression, we derive Hamilton's principle. We then focus on the freely vibrating system, formulating the corresponding variational principle and deriving Rayleigh quotients for frequency and buckling load. The presence of a magnetic field is shown to increase the frequency, while decreasing the buckling load. Finally, we derive natural and geometric boundary conditions and provide expressions for shear force and moment.

The remainder of this paper is structured as follows: In Section 2, we develop the variational principle for a double SWCNT system connected by a Winkler-Pasternak interlayer, subjected to forced vibrations, axial magnetic field, and compressive force. We then present the variational formulation of the problem, derive Hamilton's principle based on the variational formulation, and obtain expressions for Rayleigh quotients for vibration frequency and buckling load. In Section 3, we derive variationally consistent boundary conditions in terms of moment and shear force expressions, which are essential for variational and approximate solution methods, particularly those involving the Rayleigh-Ritz method. Examples of such methods applied to Kirchhoff plates can be found in Ike [63] and Nwoji et al. [64], while applications to nanotubes and nanobeams are presented in studies by Ansari et al. [65], Chakraverty and Behera [66], Behera and Chakraverty [67], and Fakher and Hosseini-Hashemi [68]. In Section 4, we provide a summary and conclude the paper.

2. Governing Equations

Double single-walled nanotube system under consideration is subject to a longitudinal magnetic field Hx acting in the axial direction and this results in the Lorentz force $\eta A_i H_x^2$ acting in the transverse direction as noted in Stamenkovic et al. [42]. In the expression $\eta A_i H_x^2$, $\eta$ is the magnetic field permeability and Ai is the cross sectional area of the ith nanotube. Nanotubes 1 and 2 are subject to axial compressive forces P1 and P2 and continuous transverse forces fi(x, t) (i=1, 2), respectively, with the external forces acting on the beams in the time domain $t_1 \leq t \leq t_2$.

The nanotubes are connected by an elastic layer between them which is modelled as a combination of a Winkler layer with an elastic modulus of k0 and a Pasternak (shear) layer with an elastic modulus of G0 as shown in Figure 1.

Figure 1. Double single-walled carbon nanotubes with Winkler and shear layer subject to axial magnetic field

The deflections of the nanotubes are given by w1(x, t) and w2(x, t) with the elastic stiffnesses denoted by EIi, cross-sectional areas by Ai and the densities by ρi. The equations governing the forced vibrations of the double nanotubes subject to a magnetic field are given in Stamenkovic et al. [42] based on the Euler–Bernoulli beam and nonlocal elasticity theories detailed in Eringen [69, 70]. These equations can be expressed as follows:

$\begin{gathered}L_1\left(w_1\right)-\mu N_1\left(w_1\right)+K\left(w_1, w_2\right)=f_1(x, t)- \mu \frac{\partial^2 f_1(x, t)}{\partial x^2}\end{gathered}$          (1)

$\begin{gathered}L_2\left(w_2\right)-\mu N_2\left(w_2\right)-K\left(w_1, w_2\right)=f_2(x, t)- \mu \frac{\partial^2 f_2(x, t)}{\partial x^2}\end{gathered}$           (2)

where, the differential operators $L_i\left(w_i\right)$ and the coupling operators $K_i\left(w_1, w_2\right)$ are given by:

$\begin{gathered}L_i\left(w_i\right)=E I_i \frac{\partial^4 w_i}{\partial x^4}+P_i \frac{\partial^2 w_i}{\partial x^2}+\rho_i A_i \frac{\partial^2 w_i}{\partial t^2}- \eta A_i H_x^2 \frac{\partial^2 w_i}{\partial x^2}\end{gathered}$              (3)

$N_i\left(w_i\right)=P_i \frac{\partial^4 w_i}{\partial x^4}+\rho_i A_i \frac{\partial^4 w_i}{\partial x^2 \partial t^2}-\eta A_i H_x^2 \frac{\partial^4 w_i}{\partial x^4}$                  (4)

$\begin{gathered}K\left(w_1, w_2\right)=k_0\left(w_1-w_2\right)-G_0\left(\frac{\partial^2 w_1}{\partial x^2}-\frac{\partial^2 w_2}{\partial x^2}\right)- \\ \mu k_0\left(\frac{\partial^2 w_1}{\partial x^2}-\frac{\partial^2 w_2}{\partial x^2}\right)+\mu G_0\left(\frac{\partial^4 w_1}{\partial x^4}-\frac{\partial^2 w_2}{\partial x^4}\right)\end{gathered}$                (5)

It is noted that in Eqs. (1) and (2), $\mu=\left(e_0 a\right)^2$ is the nonlocal parameter as defined in Eringen [69, 70]. The derivations of the Eqs. (1)-(5) are given in Stamenkovic et al. [42] with the formulations of the magnetic field acting on the nanotubes based on the works of Murmu et al. [30, 31].

3. Variational Formulation

In the present section, the variational formulation of the problem is developed. In order to formulate the variational principle applicable to the system of Eqs. (1)-(2), we first introduce the variational functionals V1(w1, w2) and V2(w1, w2) as follows

$V\left[w_1, w_2\right]=V_1\left[w_1, w_2\right]+V_2\left[w_1, w_2\right]$            (6)

where, V(w1, w2) is the variational functional to be determined. The variational functionals V1(w1, w2) and V2(w1, w2) are expressed as follows:

$\begin{gathered}V_i\left[w_1, w_2\right]=\Phi_i\left[w_i\right]-\mu \Psi_i\left[w_i\right]+\Lambda_W\left[w_1, w_2\right]+ \Lambda_P\left[w_1, w_2\right]-w_i f_i+\mu w_i \frac{\partial^2 f_i}{\partial x^2}\end{gathered}$          (7)

where the functionals $\Phi_i\left[w_i\right]$, $\Psi_i\left[w_i\right]$, $\Lambda_W\left[w_1, w_2\right]$ and $\Lambda_P\left[w_1, w_2\right]$ are defined as follows:

$\Phi_i\left[w_i\right]=\frac{1}{2} \int_{t_1}^{t_2} \int_0^L\left[\begin{array}{c}E_i I_i\left(\frac{\partial^2 w_i}{\partial x^2}\right)^2-P_i\left(\frac{\partial w_i}{\partial x}\right)^2 \\ -\rho_i A_i\left(\frac{\partial w}{\partial t}\right)^2+\eta A_i H_x^2\left(\frac{\partial w}{\partial x}\right)^2\end{array}\right] d x d t$            (8)

$\Psi_i\left[w_i\right]=\frac{1}{2} \int_{t_1}^{t_2} \int_0^L\left[\begin{array}{c}P_i\left(\frac{\partial^2 w_i}{\partial x^2}\right)^2+\rho_i A_i \frac{\partial^2 w_i}{\partial x^2} \frac{\partial^2 w_i}{\partial t^2} \\ -\eta A_i H_x^2\left(\frac{\partial^2 w_i}{\partial x^2}\right)^2\end{array}\right] d x d t$             (9)

$\Lambda_W\left[w_1, w_2\right]=\frac{k_0}{2} \int_{t_1}^{t_2} \int_0^L\left[\begin{array}{c}\left(w_1-w_2\right)^2 \\ +\mu\left(\frac{\partial w_1}{\partial x}-\frac{\partial w_2}{\partial x}\right)^2\end{array}\right] d x d t$           (10)

$\Lambda_P\left[w_1, w_2\right]=\frac{G_0}{2} \int_{t_1}^{t_2} \int_0^L\left[\begin{array}{c}\left(\frac{\partial w_1}{\partial x}-\frac{\partial w_2}{\partial x}\right)^2 \\ +\mu\left(\frac{\partial^2 w_1}{\partial x^2}-\frac{\partial^2 w_2}{\partial x^2}\right)^2\end{array}\right] d x d t$            (11)

Thus we have:

$\begin{gathered}V_1\left(w_1, w_2\right)=\Phi_1\left[w_1\right]-\mu \Psi_1\left[w_1\right]+\Lambda\left[w_1, w_2\right]- w_1\left(f_1-\mu \frac{\partial^2 f_1}{\partial x^2}\right)\end{gathered}$            (12)

$\begin{gathered}V_2\left(w_1, w_2\right)=\Phi_2\left[w_2\right]-\mu \Psi_2\left[w_2\right]+\Lambda\left[w_1, w_2\right]- w_2\left(f_2-\mu \frac{\partial^2 f_2}{\partial x^2}\right)\end{gathered}$            (13)

where, $\Lambda\left[w_1, w_2\right]$ is defined as follows:

$\Lambda\left[w_1, w_2\right]=\Lambda_W\left[w_1, w_2\right]+\Lambda_P\left[w_1, w_2\right]$           (14)

It is observed that the Euler-Lagrange equations of the variational functional V(w1, w2) given by Eq. (6) corresponds to the governing equations of the double nanotube system given by Eqs. (1)-(2). This can be verified by taking the first variation of the functional $\delta V\left[w_1, w_2\right]$ with respect to w1 and w2.

4. Hamilton’s Principle

The Hamilton’s principle can be expressed as:

$\begin{gathered}\int_{t_1}^{t_2}\left[\delta K E(t)-\left(\delta W_E(t)+\delta P E_1(t)+\right.\right. \left.\left.\delta P E_2(t)\right)\right] d t=0\end{gathered}$           (15)

In the present problem, the functionals KE(t), WE(t), PE1(t) and PE2(t) are given by:

$\begin{gathered}K E(t)=\frac{1}{2} \sum_{i=1}^2 \int_0^L\left[\rho_i A_i\left(\frac{\partial w_i}{\partial t}\right)^2+\right. \left.\mu \rho_i A_i \frac{\partial^2 w_i}{\partial x^2} \frac{\partial^2 w_i}{\partial t^2}\right] d x\end{gathered}$             (16)

$W_E(t)=-\sum_{i=1}^2 \int_0^L\left(f_i-\mu \frac{\partial^2 f_i}{\partial x^2}\right) w_i(x, t) d x$               (17)

$P E_1(t)=\frac{1}{2} \sum_{i=1}^2 \int_0^L\left[\begin{array}{c}E_i I_i\left(\frac{\partial^2 w_i}{\partial x^2}\right)^2-P_i\left(\frac{\partial w_i}{\partial x}\right)^2 \\ +\eta A_i H_x^2\left(\frac{\partial w_i}{\partial x}\right)^2\end{array}\right] d x$             (18)

$P E_2(t)=\frac{1}{2} \sum_{i=1}^2 \int_0^L\left[\begin{array}{c}k_0\left(w_1-w_2\right)^2+\mu k_0\left(\frac{\partial w_1}{\partial x} \frac{\partial w_2}{\partial x}\right)^2 \\ +G_0\left(\frac{\partial w_1}{\partial x}-\frac{\partial w_2}{\partial x}\right)^2 \\ +\mu G_0\left(\frac{\partial^2 w_1}{\partial x^2}-\frac{\partial^2 w_2}{\partial x^2}\right)^2\end{array}\right] d x$                 (19)

In Eqs. (16)-(19), KE(t) is the kinetic energy of the double nanotube system, WE(t) is the work done by external forces and PE1(t) is the potential energy of deformation which includes the contributions by the compressive forces Pi and the magnetic field $\eta A_i H_x^2$. PE2(t) is the potential energy due to the Winkler and Pasternak layers between the nanotubes.

5. Free Vibrations

Variational principle for a freely vibrating double beam system is formulated next. In this case the external forces fi(x,t)=0 with i=1, 2. For freely vibrating beams, the deflection can be expressed as:

$w_i(x, t)=W_i(x) e^{i \omega t}$               (20)

where, ω is the vibration frequency. The differential equations for the freely vibrating double nanotube system can be obtained from the system of Eqs. (1)-(5) and can be expressed as:

$L_{F V 1}\left(W_1\right)-\mu N_{F V 1}\left(W_1\right)+K\left(W_1, W_2\right)=0$               (21)

$L_{F V 2}\left(W_2\right)-\mu N_{F V 2}\left(W_2\right)-K\left(W_1, W_2\right)=0$                (22)

In Eqs. (21)-(22), the differential operators $L_{F V i}\left(W_i\right)$ and $N_{F V i}\left(W_i\right)$ are given by:

$\begin{gathered}L_{F V i}\left(W_i\right)=E_i I_i \frac{d^4 W_i}{d x^4}+P_i \frac{d^2 W_i}{d x^2}-\omega^2 \rho_i A_i W_i- \eta A_i H_x^2 \frac{d^2 W_i}{d x^2}\end{gathered}$             (23)

$N_{F V i}\left(W_i\right)=P_i \frac{d^4 W_i}{d x^4}-\omega^2 \rho_i A_i \frac{d^2 W_i}{d x^2}-\eta A_i H_x^2 \frac{d^4 W_i}{d x^4}$             (24)

The coupling operator $K\left(W_1, W_2\right)$ is defined as:

$\begin{gathered}K\left(W_1, W_2\right)=k_0\left(W_1-W_2\right)-G_0\left(\frac{d^2 W_1}{d x^2}-\frac{d W_2}{d x^2}\right)- \mu k_0\left(\frac{d^2 W_1}{d x^2}-\frac{d^2 W_2}{d x^2}\right)+\mu G_0\left(\frac{d^4 W_1}{d x^4}-\frac{d^2 W_2}{d x^4}\right)\end{gathered}$         (25)

For the freely vibrating case, the variational functional $V_{F V}\left[w_1, w_2\right]$ is of the form:

$V_{F V}\left[W_1, W_2\right]=V_{F V 1}\left[W_1, W_2\right]+V_{F V 2}\left[W_1, W_2\right]$           (26)

with the variational functionals $V_{F V 1}\left[W_1, W_2\right]$ and $V_{F V 2}\left[W_1, W_2\right]$ given by:

$\begin{aligned} V_{F V i}\left[W_1, W_2\right]= & \Phi_{F V i}\left[W_i\right]-\mu \Psi_{F V i}\left[W_i\right] +\Lambda_{F V}\left[W_1, W_2\right]\end{aligned}$          (27)

Functionals $\Phi_{F V i}\left(W_i\right)$, $\Psi_{F V i}\left[W_i\right]$ and $\Lambda_{F V}\left[W_1, W_2\right]$ in Eq. (27) are given by:

$\Phi_{F V i}\left[W_i\right]=\frac{1}{2} \int_0^L\left[\begin{array}{c}E_i I_i\left(\frac{d^2 W_i}{d x^2}\right)^2-P_i\left(\frac{d W_i}{d x}\right)^2 \\ -\omega^2 \rho_i A_i W_i^2+\eta A_i H_x^2\left(\frac{d W_i}{d x}\right)^2\end{array}\right] d x$           (28)

$\Psi_{F V i}\left[W_i\right]=\frac{1}{2} \int_0^L\left[\begin{array}{c}P_i\left(\frac{d^2 W_i}{d x^2}\right)^2+\omega^2 \rho_i A_i\left(\frac{d W_i}{d x}\right)^2 \\ -\eta A_i H_x^2\left(\frac{d^2 W_i}{d x^2}\right)^2\end{array}\right] d x$             (29)

$\Lambda_{F V}\left[W_1, W_2\right]=\Lambda_{F V 1}\left[W_1, W_2\right]+\Lambda_{F V 2}\left[W_1, W_2\right]$                (30)

where,

$\begin{gathered}\Lambda_{F V 1}\left[W_1, W_2\right]=\frac{k_0}{2} \int_0^L\left[\left(W_1-W_2\right)^2+\mu\left(\frac{d W_1}{d x}-\right.\right. \left.\left.\frac{d W_2}{d x}\right)^2\right] d x\end{gathered}$                  (31)

$\Lambda_{F V 2}\left[W_1, W_2\right]=\frac{G_0}{2} \int_0^L\left[\begin{array}{c}\left(\frac{d W_1}{d x}-\frac{d W_2}{d x}\right)^2 \\ +\mu\left(\frac{d^2 W_1}{d x^2}-\frac{d^2 W_2}{d x^2}\right)^2\end{array}\right] d x$                     (32)

5.1 Rayleigh quotient for frequency

Next Rayleigh quotient is given for the freely vibrating double beam system. For this purpose, we introduce the functionals $Y_{F V i}\left(W_i\right)$ and $M_{F V i}\left(W_i\right)$ given by:

$\begin{gathered}Y_{F V i}\left(W_i\right)=\left(P_i+\eta A_i H_x^2\right) \int_0^L\left[\left(\frac{d W_i}{d x}\right)^2\right. +\left.\mu\left(\frac{d^2 W_i}{d x^2}\right)^2\right] d x\end{gathered}$              (33)

$M_{F V i}\left(W_i\right)=\int_0^L E_i I_i\left(\frac{d^2 W_i}{d x^2}\right)^2 d x$                 (34)

Using the variational principle given in Eq. (26), Rayleigh quotient for the vibration frequency ω can now be formulated and is given by the expression:

$\omega^2=\frac{\sum_{i=1}^2\,\,\left[M_{F V i}\,\left(W_i\right)-Y_{F V i}\,\left(W_i\right)\right]+\Lambda_{F V}\,\left[W_1, W_2\right]}{\sum_{i=1}^2 \int_0^L \rho_i A_i\left[W_i{ }^2-\mu\left(\frac{d W_i}{d x}\right)^2\right] d x}$                (35)

Frequency ω is computed by minimizing Eq. (35) with respect to Wi. In Eq. (35), $M_{F V i}\left(W_i\right)$ is defined by Eq. (34), $Y_{F V i}\left(W_i\right)$ by Eq. (33) and $\Lambda_{F V}\left[W_1, W_2\right]$ by Eq. (30). It is observed that the effect of the compressive forces Pi and the magnetic field $\eta A_i H_x^2$ is to reduce the frequency. This is due to the fact that the functional $Y_{F V i}\left(W_i\right)$ given by Eq. (33) is a positive definite functional and appears with a negative sign in the numerator of the frequency expression (35). The nonlocal effect appears in the denominator of Eq. (35) with a minus sign and in the numerator with a negative and a positive sign (in $\Lambda_{F V}\left[W_1, W_2\right]$) and as such its effect on the frequency depends on the relative magnitudes of the problem parameters.

5.2 Rayleigh quotient for buckling load

Next, Rayleigh quotient for buckling load is formulated. For this purpose, we introduce the functionals:

$Z_{F V i}\left(W_i\right)=\int_0^L \rho_i A_i\left[W_i^2-\mu\left(\frac{d W_i}{d x}\right)^2\right] d x$             (36)

$B_{F V i}\left(W_i\right)=\int_0^L\left[\left(\frac{d W_i}{d x}\right)^2+\mu\left(\frac{d^2 W_i}{d x^2}\right)^2\right] d x$               (37)

Rayleigh quotient for the buckling load for the case $P_1=P_2=P$ can be expressed as:

$P=\frac{\sum_{i=1}^2\,\left[\begin{array}{c}M_{F V i}\,\left(W_i\right)-\omega^2 Z_{F V i}\,\left(W_i\right) -\eta H_x^2 A_i B_{F V i}\,\left(W_i\right) +\Lambda_{F V}\left(W_1, W_2\right)\end{array}\,\,\,\right]}{\sum_{i=1}^2 B_{F V i}\,\,\left(W_i\right)}$              (38)

where, $Z_{F V i}\left(W_i\right)$ is given by Eq. (36) and $B_{F V i}\left(W_i\right)$ by Eq. (37). The buckling load P is computed by minimizing Eq. (38) with respect to Wi. It is observed that the effect of the magnetic field Hx on the buckling load is to reduce the buckling load. This is due to the fact that the functional $B_{F V i}\left(W_i\right)$ given by Eq. (37) is a positive definite functional and appears with a negative sign in the numerator of the buckling expression (38). The effect of the nonlocal parameter η on the buckling load cannot be deduced from the buckling expression (38) and can be positive or negative depending on the relative values of various quantities.

6. Boundary Conditions

In this section, natural and geometric boundary conditions are derived using the variational formulation (Eq. (25)) of the freely vibrating double beam system. The first variations of $V_{F V}\left(W_1, W_2\right)$ with respect to Wi, denoted as $\delta_{W_i} V_{F V}$, can be obtained by integration by parts. We note that:

$\begin{gathered}\delta_{w_1} V_{F V}=\delta_{w_1} V_{F V 1}+\delta_{w_1} V_{F V 2}=  \int_0^L D_1\left(W_1, W_2\right) \delta W_1 d x+\partial \Omega_1(0, L)\end{gathered}$             (39a)

$\begin{gathered}\delta_{w_2} V_{F V}=\delta_{w_2} V_{F V 1}+\delta_{w_2} V_{F V 2}= \int_0^L D_2\left(W_1, W_2\right) \delta W_2 d x+\partial \Omega_2(0, L)\end{gathered}$           (39b)

where

$\begin{gathered}D_i\left(W_1, W_2\right)=L_i\left(W_i\right)-\mu N_i\left(W_i\right)+ (-1)^{i+1} K\left(W_1, W_2\right)\end{gathered}$              (40)

In Eqs. (38)-(39), $\partial \Omega_i(0, L)$ denotes the boundary terms with $\delta_{W_i} V_{F V i}\left(W_1, W_2\right)$ given by:

$\begin{aligned} \delta_{W_i} V_{F V i}\left(W_1, W_2\right)= & \int_0^L\left[\begin{array}{c}\frac{\partial V_{F V i}}{\partial W_i} \frac{d}{d x}\left(\frac{\partial V_{F V i}}{\partial W_{i x}}\right) \\ +\frac{d^2}{d x^2}\left(\frac{\partial V_{F V i}}{\partial W_{i x x}}\right)\end{array}\right] \delta W_i d x+ \partial \Omega_i(0, L)\end{aligned}$             (41)

and $\partial \Omega_i(0, L)$ is defined as:

$\partial \Omega_i(0, L)=\left.\left(Q_i \delta W_i\right)\right|_{x=0}^{x=L}+\left.\left(M_i \delta W_{i x}\right)\right|_{x=0}^{x=L}$              (42)

where,

$\begin{gathered}Q_i=-\left(P_i-\eta A_i H_x^2+\mu \omega^2 \rho_i A_i\right) \frac{d W_i}{d x}+\mu k_0\left(\frac{d W_1}{d x}-\right. \left.\frac{d W_2}{d x}\right)+\mu G_0\left(\frac{d W_1}{d x}-\frac{d W_2}{d x}\right)\end{gathered}$               (43)

$\begin{gathered}M_i=E_i I_i+\mu\left(P_i-\eta A_i H_x^2\right) \frac{d^2 W_i}{d x^2}+\mu G_0\left(\frac{d^2 W_1}{d x^2}-\right. \left.\frac{d^2 W_2}{d x^2}\right)\end{gathered}$               (44)

Here Qi is the shear force and Mi is the moment expression. Thus, the boundary conditions at x=0 and x=L are given by:

$\begin{gathered}W_i \text { or } Q_i(x)=-\left(P_i-\eta A_i H_x^2+\mu \omega^2 \rho_i A_i\right) \frac{d W_i}{d x}+ \mu k_0\left(\frac{d W_1}{d x}-\frac{d W_2}{d x}\right)-\mu G_0\left(\frac{d W_1}{d x}-\frac{d W_2}{d x}\right) \text { specified }\end{gathered}$              (45)

$\begin{gathered}W_{i x} \text { or } M_i(x)=E_i I_i+\mu\left(P_i-\eta A_i H_x^2\right) \frac{d^2 W_i}{d x^2}+ G_0\left(\frac{d^2 W_1}{d x^2}-\frac{d^2 W_2}{d x^2}\right) \text { specified }\end{gathered}$              (46)

with i=1, 2. It is observed that due to Pasternak interlayer between the nanotubes, the boundary conditions are coupled if the shear force and/or moment is specified at x=0 or x=L as given by Eqs. (45)-(46).

Derivations of the variational expressions (26)-(29) and the boundary conditions (45), (46) are given next. We note that:

$\begin{gathered}\int_0^L E_i I_i \frac{d^4 W_i}{d x^4} \delta W_i d x=B_{1 i F V}\left(W_i, \delta W_i\right)+ \delta\left[\frac{1}{2} \int_0^L E_i I_i\left(\frac{d^2 W_i}{d x^2}\right)^2 d x\right]\end{gathered}$            (47)

where,

$\begin{gathered}B_{1 F V i}\left(W_i, \delta W_i\right)=\left.\left(E_i I_i \frac{d^3 W_i}{d x^3} \delta W_i\right)\right|_{x=0} ^{x=L}- \left.\left(E_i I_i \frac{d^2 W_i}{d x^2} \delta\left(\frac{d W_i}{d x}\right)\right)\right|_{x=0} ^{x=L}\end{gathered}$            (48)

Similarly,

$\begin{gathered}\int_0^L\left(P_i-\eta A_i H_x^2\right) \frac{d^2 W_i}{d x^2} \delta W_i d x=B_{2 i F V}\left(W_i, \delta W_i\right)- \delta\left[\frac{1}{2} \int_0^L\left(P_i-\eta A_i H_x^2\right)\left(\frac{d W_i}{d x}\right)^2 d x\right]\end{gathered}$             (49)

where,

$B_{2 F V i}\left(W_i, \delta W_i\right)=\left.\left(\left(P_i-\eta A_i H_x^2\right) \frac{d W_i}{d x} \delta W_i\right)\right|_{x=0} ^{x=L}$                (50)

Eqs. (45)-(47) and (49) indicate the variational expressions in Eqs. (34)-(35). Similarly, Eqs. (46), (48), and (50) have the boundary terms which appear in the boundary conditions shown in Eqs. (43)-(44).

7. Conclusions

A system of two partial differential equations govern the vibrations of a system of double single-walled carbon nanotubes with the connection between the nanotubes modelled as an elastic interlayer. The double nanotube system is subject to a longitudinal magnetic field and compressive forces. In the present study, the constitutive equations are based on nonlocal Euler-Bernouilli beam theory and the interlayer between the nanotubes is defined in the form of a Winkler-Pasternak layer with the expressions involving the displacements (Winkler) and second derivatives of displacements (Pasternak). As the constitutive relations are based on nonlocal elasticity, the formulation leads to fourth order derivatives describing the interlayer. The variational formulation of the problem is obtained for the coupled differential equations with the system subject to forced vibrations. For use in the approximate and variational methods of solutions of the problem, vibration frequency and the buckling load are expressed in the form of Rayleigh quotients. A study of the Rayleigh quotients indicates that the effect of the magnetic field is to reduce the frequency and the buckling load. An important part of the present study is the derivation of variationally consistent boundary conditions. Natural boundary conditions are obtained involving the moment and shear force expressions. It is observed that the presence of a Pasternak interlayer between the nanotubes leads to coupled boundary conditions. The results presented in the study can be used in a number of approximate and numerical methods of solution, especially in the presence of different boundary conditions.

It is noted that variational formulation of a problem defined by differential equations provide the weak formulation of the problem and the approximations based on this formulation are the weak solutions. However, weak formulations can provide accurate solutions of the differential equations in the absence of an analytical solution as noted in Nicolescu and Bobe [71]. Future work on the present study could involve obtaining the approximate solutions of the nanotube system subject to a combination of clamped, free, simply supported and rotational/torsional boundary conditions using the variational formulation of the present study. This would require choosing suitable trial functions for the approximate solution of the specific problem based on the applicable boundary conditions. Another direction in further work can involve the study of the second variation of the variational functional and determining the conditions under which it is positive-definite indicating a minimum. Moreover, present formulation is based on Euler-Bernoulli beam theory and further work can extend the results formulation based on Timoshenko beam theory in order to take the shear effects into consideration.

Acknowledgment

The research reported in this paper was supported by research grants from the University of KwaZulu-Natal (UKZN) and from National Research Foundation (NRF) of South Africa. The author gratefully acknowledges the supports provided by UKZN and NRF.

  References

[1] Kibalchenko, M., Payne, M.C., Yates J.R. (2011). Magnetic response of single-walled carbon nanotubes induced by an external magnetic field. ACS Nano, 5(1): 537–545. https://doi.org/10.1021/nn102590b

[2] Zhang, D., Lei, Y., Shen, Z. (2018). Effect of longitudinal magnetic field on vibration response of double-walled carbon nanotubes embedded in viscoelastic medium. Acta Mechanica Solida Sinica, 31: 187–206. https://doi.org/10.1007/s10338-018-0006-x

[3] Arani, A.G., Roudbari, M.A., Amir, S. (2016). Longitudinal magnetic field effect on wave propagation of fluid-conveyed SWCNT using Knudsen number and surface considerations. Appl Math Model. 40(3): 2025–2038. https://doi.org/10.1016/j.apm.2015.09.055

[4] Alazwari, M.A., Esen, I., Abdelrahman, A.A., Abdraboh, A.M., Eltaher, M.A. (2022). Dynamic analysis of functionally graded (FG) nonlocal strain gradient nanobeams under thermo-magnetic fields and moving load. Advances in Nano Research, 12(3): 231-251. https://doi.org/10.12989/ANR.2022.12.3.231.

[5] Kiani, K. (2014). Magnetically affected single-walled carbon nanotubes as nanosensors. Mechanics Research Communications, 60: 33-39. https://doi.org/10.1016/j.mechrescom.2014.05.005

[6] Samadishadlou, M., Farshbaf, M., Annabi, N., Kavetskyy, T., Khalilov, R., Saghfi, S., Akbarzadeh, A., Mousavi, S. (2018) Magnetic carbon nanotubes: Preparation, physical properties, and applications in biomedicine. Artificial Cells, Nanomedicine, and Biotechnology, 46(7): 1314-1330. https://doi.org/10.1080/21691401.2017.1389746

[7] Raniszewski, G., Miaskowski, A., Wiak, S. (2015). The application of carbon nanotubes in magnetic fluid hyperthermia. Journal of Nanomaterials, 2015: Article ID 527652. https://doi.org/10.1155/2015/527652.

[8] Bellucci, S., González, J., Guinea, F., Onorato, P., Perfetto, E. (2007). Magnetic field effects in carbon nanotubes. Journal of Physics: Condensed Matter, 19(39): 1-10. https://doi.org/10.1088/0953-8984/19/39/395017

[9] Klinovaja, J., Schmidt M.J., Braunecker B., Loss D. (2011). Carbon nanotubes in electric and magnetic fields. PhysicsxReviewsxB, 84(8): 1-22. https://doi.org/10.1103/PhysRevB.84.085452

[10] Kibalchenko, M., Payne, M.C., Yates, J.R. (2011). Magnetic response of single-walled carbon nanotubes induced by an external magnetic field. ACS Nano, 5(1): 537-545. https://doi.org/10.1021/nn102590b

[11] Kiani, K. (2015). Vibrations and instability of pretensioned current-carrying nanowires acted upon by a suddenly applied three-dimensional magnetic field. Materials Chemistry Physics, 62:531-541. https://doi.org/10.1016/j.matchemphys.2015.05.080

[12] Fedorov, G., Tselev, A., Jimenez, D., Latil, S., Kalugin, N.G., Barbara, P., Smirnov, D., Roche, S. (2007). Magnetically induced field effect in carbon nanotube devices. Nano Letters, 7(4): 960-964. https://doi.org/10.1021/nl063029v

[13] Popov, A.M., Lebedeva, I.V., Knizhnik, A.A., Lozovik, Y.E., Poklonski, N.A., Siahlo, A.I., Vyrko, S.E., Ratkevich, S.V. (2014). Force and magnetic field sensor based on measurement of tunneling conductance between ends of coaxial carbon nanotubes. Computational Materials Science, 92: 84-91. https://doi.org/10.1016/j.commatsci.2014.05.020

[14] Mandal, P., Chopra, V., Ghosh, A. (2015). Independent positioning of magnetic nanomotors. ACS Nano, 9(5): 4717-4725. https://doi.org/10.1021/acsnano.5b01518

[15] Pal, M., Somalwar, N., Singh, A., Bhat, R., Eswarappa, S.M., Saini, D.K., Ghosh, A. (2018). Maneuverability of magnetic nanomotors inside living cells. Advanced Materials, 30(22): 1-7. https://doi.org/10.1002/adma.201800429

[16] Bellucci, S., Gonzalez, J., Guinea, F., Onorato, P., Perfetto, E. (2007). Magnetic field effects in carbon nanotubes, Journal of Physics: Condensed Matter, 19: 395017. https://doi.org/10.1088/0953-8984/19/39/395017

[17] Kibalchenko, M., Payne, M.C., Yates, J.R. (2011). Magnetic response of single-walled carbon nanotubes induced by an external magnetic field. ACS Nano, 5(1): 537-545. https://doi.org/10.1021/NN102590B

[18] Li, S., Xie, H.J., Wang, X. (2011). Dynamic characteristics of multi-walled carbon nanotubes under a transverse magnetic field. Bulletin of Materials Science, 34: 45-52. https://doi.org/10.1007/s12034-011-0033-4

[19] Wang, X., Shen, J.X., Liu, Y., Shen, G.G., Lu, G. (2012). Rigorous van der Waals effect on vibration characteristics of multi-walled carbon nanotubes under a transverse magnetic field. Applied Mathematical Modelling, 36: 648-656. https://doi.org/10.1016/j.apm.2011.07.017

[20] Wang, H., Dong, K., Men, F., Yan, Y.J., and Wang, X. (2010). Influences of longitudinal magnetic field on wave propagation in carbon nanotubes embedded in elastic matrix. Applied Mathematical Modelling, 34: 878–889. https://doi.org/10.1016/j.apm.2009.07.005 

[21] Wang, B., Deng, Z., Ouyang, H., Wang, Y. (2015). Terahertz wave propagation in a fluid-conveying single-walled carbon nanotube with initial stress subjected to temperature and magnetic fields. Acta Mechanica, 226(9): 3031-3043. https://doi.org/10.1007/s00707-015-1367-6

[22] Kiani, K. (2012). Transverse wave propagation in elastically confined single-walled carbon nanotubes subjected to longitudinal magnetic fields using nonlocal elasticity models. Physica E: Low-Dimensional SystemsxandxNanostructures, 45: 86-96. https://doi.org/10.1016/j.physe.2012.07.015

[23] Narendar, S., Gupta, S., Gopalakrishnan, S. (2012). Wave propagation in single-walled carbon nanotube under longitudinal magnetic field using nonlocal Euler–Bernoulli beam theory. Applied Mathematical Modelling, 36(9): 4529-4538. https://doi.org/10.1016/j.apm.2011.11.073

[24] Arani, A.G., Roudbari, M.A. Amir, S. (2016). Longitudinal magnetic field effect on wave propagation of fluid-conveyed SWCNT using Knudsen number and surface considerations. Applied Mathematical Modelling, 40(3): 2025–2038. https://doi.org/10.1016/j.apm.2015.09.055

[25] Arani, A.G., Jalilvand, A., Kolahchi, R. (2014). Wave propagation of magnetic nanofluid-conveying double-walled carbon nanotubes in the presence of longitudinal magnetic field. Proceedings of the Institute of Mechanical Engineers, Part N, Journal of Nanoengineering and Nanosystems, 228(2): 82-92. https://doi.org/10.1177/1740349913488575

[26] Li, L., Hu, Y., Ling, L. (2016). Wave propagation in viscoelastic single-walled carbon nanotubes with surface effect under magnetic field based on nonlocal strain gradient theory. Physica E: Low-dimensional Systems and Nanostructures, 75: 118-124. https://doi.org/10.1016/j.physe.2015.09.028

[27] Zhen, Y., Zhou, L. (2017). Wave propagation in fluid-conveying viscoelastic carbon nanotubes under longitudinal magnetic field with thermal and surface effect via nonlocal strain gradient theory. Modern Physics Letters B, 31(8): 1750069.

https://doi.org/10.1142/s0217984917500695

[28] Kiani, K. (2014). Longitudinally varying magnetic field influenced transverse vibration of embedded double-walled carbon nanotubes. International Journal of Material Science, 87: 179–199. https://doi.org/10.1016/j.ijmecsci.2014.04.018

[29] Kiani, K. (2013). Characterization of free vibration of elastically supported double-walled carbon nanotubes subjected to a longitudinally varying magnetic field. Acta Mechanica, 224(12): 3139-3151. https://doi.org/10.1007/s00707-013-0937-8]

[30] Murmu, T., McCarthy, M.A., Adhikari, S. (2012). Vibration response of double-walled carbon nanotubes subjected to an externally applied longitudinal magnetic field: a nonlocal elasticity approach. Journal of Sound and Vibration, 331(23): 5069-5086. https://doi.org/10.1016/j.jsv.2012.06.005

[31] Murmu, T., McCarthy, M.A., Adhikari, S. (2012). Nonlocal elasticity based magnetic field affected vibration response of double single-walled carbon nanotube systems. Journal of Applied Physics, 111: 113511. http://doi.org/10.1063/1.4720084

[32] Arani, G.A., Amir, S., Dashti, P., Yousefi, M. (2014). Flow induced vibration of double bonded visco-CNTs under magnetic fields considering surface effect. Computational Materials Science, 86: 144-154. https://doi.org/10.1016/j.commatsci.2014.01.047

[33] Kiani, K. (2014). Vibration and instability of a single-walled carbon nanotube in a three-dimensional magnetic field. Journal of Physics and Chemistry of Solids, 75(1): 15-22. https://doi.org/10.1016/j.jpcs.2013.07.022

[34] Jena, S.K., Chakraverty, S., Malikan, M. (2020). Vibration and buckling characteristics of nonlocal beam placed in a magnetic field embedded in Winkler–Pasternak elastic foundation using a new refined beam theory: an analytical approach. The European Physical Journal Plus, 135: 164. https://doi.org/10.1140/epjp/s13360-020-00176-3

[35] Esen, I., Abdelrhmaan, A.A., Eltaher, M.A., (2021). Free vibration and buckling stability of FG nanobeams exposed to magnetic and thermal fields. Engineering with Computers. https://doi.org/10.1007/s00366-021-01389-5

[36] Chang, T.P. (2016). Nonlinear vibration of single-walled carbon nanotubes under magnetic field by stochastic finite element method. International Journal of Structural Stability and Dynamics, 16(8): 1-18. https://doi.org/10.1142/S0219455415500467

[37] Ebrahimi, F., Hosseini, S.H.S. (2021). Nonlinear vibration and dynamic instability analysis nanobeams under thermo-magneto-mechanical loads: a parametric excitation study. Engineering with Computers, 7: 395-408. https://doi.org/10.1007/s00366-019-00830-0

[38] Anh, N.D., Hieu, D.V. (2021). Nonlinear vibration of nonlocal strain gradient nanotubes under longitudinal magnetic field. Vietnam Journal of Mechanics, 43(1): 55–77. https://doi:10.15625/0866-7136/15467

[39] Jalaei, M.H., Arani, A.G., Nguyen-Xuan, H. (2019). Investigation of thermal and magnetic field effects on the dynamic instability of FG Timoshenko nanobeam employing nonlocal strain gradient theory. International Journal of Mechanical Science, 161: 105043. https://doi.org/10.1016/j.ijmecsci.2019.105043

[40] Bahaadini, R., Hosseini, M., Amiri, M. (2021). Dynamic stability of viscoelastic nanotubes conveying pulsating magnetic nanoflow under magnetic field. Engineering with Computers, 37: 2877-2889. https://doi.org/10.1007/s00366-020-00980-6

[41] Nasirshoaibi, M., Mohammadi, N., Nasirshoaibi, M. (2015). Forced transverse vibration of a closed double single-walled carbon nanotube system containing a fluid with effect of compressive axial load. Shock and Vibration, 2015: Article ID 435284. https://doi.org/10.1155/2015/435284

[42] Stamenkovic, M., Karlicic, D., Goran, J., Kozic, P. (2016). Nonlocal forced vibration of a double single-walled carbon nanotube system under the influence of an axial magnetic field. Journal of Mechanics of Materials and Structures, 11(3): 279-307. https://doi.org/10.2140/jomms.2016.11.279

[43] Ebrahimi. F., Dabbagh, A. (2021). Magnetic field effects on thermally affected propagation of acoustical waves in rotary double-nanobeam systems. Waves in Random and Complex Media, 31(1): 25-45. https://doi.org/10.1080/17455030.2018.1558308

[44] Yang, Y., Lim, C.W. (2011). A variational principle approach for buckling of carbon nanotube on nonlocal Timoshenko beam models. Nano, 6(4): 363-377. https://doi.org/10.1142/S1793292011002676

[45] Barretta, R., Feo, L., Luciano, R., Marotti de Sciarra, F. (2015). Variational formulations for functionally graded nonlocal Bernoulli–Euler nanobeams. Composite Structures, 129: 80-89. https://doi.org/10.1016/j.compstruct.2015.03.033

[46] Barretta, R., Marotti de Sciarra, F. (2019). Variational nonlocal gradient elasticity for nano-beams. International Journal of Engineering Science, 143: 73-91. https://doi.org/10.1016/j.ijengsci.2019.06.016

[47] Challamel, N. (2013). Variational formulation of gradient or/and nonlocal higher-order shear elasticity beams. Composite Structures, 105: 351-368. https://dx.doi.org/10.1016/j.compstruct.2013.05.026

[48] Ike, C.C. (2018). Energy formulation for flexural-torsional buckling of thin-walled column with open cross- section. Mathematical Modeling of Engineering Problems, 5(2): 58-66. https://doi.org/10.18280/mmep.050202

[49] Ike, C.C. (2019). Timoshenko beam theory for the flexural analysis of moderately thick beams–Variational formulation, and closed form solution. Italian Journal of Engineering Science, 63(1): 34-45.

[50] Shi, G., Voyiadjis, G.Z. (2011). A sixth-order theory of shear deformable beams with variationally consistent boundary conditions. ASME Journal of Applied Mechanics, 78: 021019-1-11. https://doi.org/10.1115/1.4002594

[51] Robinson, M.T.A., Adali, S. (2016). Variational solution for buckling of nonlocal carbon nanotubes under uniformly and triangularly distributed axial loads. Composite Structures, 156: 101-107. https://doi.org/10.1016/j.compstruct.2016.01.026

[52] Barretta, R., Faghidian, S.A., Marotti de Sciarra, F., Pinnola, F.P. (2021). Timoshenko nonlocal strain gradient nanobeams: Variational consistency, exact solutions and carbon nanotube Young moduli. Mechanics of Advanced Materials and Structures, 28(15): 1523-1536. https://doi.org/10.1080/15376494.2019.1683660

[53] Yu, Y.J., Zhang, K., Deng, Z.C. (2019). Buckling analyses of three characteristic-lengths featured size-dependent gradient-beam with variational consistent higher order boundary conditions. Applied Mathematical Modelling, 74: 1-20. https://doi.org/10.1016/j.apm.2019.04.019

[54] Xu, X., Zheng, M. (2019). Analytical solutions for buckling of size-dependent Timoshenko beams. Applied Mathematics and Mechanics-Engl. Ed., 40: 953–976. https://doi.org/10.1007/s10483-019-2494-8

[55] Pinnola, F.P., Faghidian, S.A., Barretta, R., Marotti de Sciarra, F. (2020). Variationally consistent dynamics of nonlocal gradient elastic beams. International Journal of Engineering Science, 149: 103220. https://doi.org/10.1016/j.ijengsci.2020.103220

[56] Barretta, R., Faghidian, S.A., Marotti de Sciarra, F., Pinnola, F.P. (2021). Timoshenko nonlocal strain gradient nanobeams: Variational consistency, exact solutions and carbon nanotube Young moduli. Mechanics of Advanced Materials and Structures, 28(15): 1523-1536, https://doi.org/10.1080/15376494.2019.1683660

[57] Adali, S. (2008). Variational principles for multiwalled carbon nanotubes undergoing buckling based on nonlocal elasticity theory. Physics Letters A, 372: 5701-5705. https://doi.org/10.1016/j.physleta.2008.07.003

[58] Adali, S. (2009). Variational principles for transversely vibrating multiwalled carbon nanotubes based on nonlocal Euler-Bernoulli beam model. Nano Letters, 9: 1737–1741. https://doi.org/10.1021/nl8027087

[59] Xu, X.J, Deng, Z.C. (2014). Variational principles for buckling and vibration of MWCNTs modeled by strain gradient theory. Applied Mathematics and Mechanics, English Ed., 35(9): 1115-1128. https://doi.10.1007/s10483-014-1855-6

[60] Kucuk, I., Sadek, I.S., Adali, S. (2010). Variational principles for multiwalled carbon nanotubes undergoing vibrations based on nonlocal Timoshenko beam theory. Journal of Nanomaterials, 2010: 461252. https://doi.org/10.1155/2010/461252

[61] Adali, S. (2012). Variational formulation for buckling of multiwalled carbon nanotubes modelled as nonlocal Timoshenko beams. Journal of Theoretical and Applied Mechanics, 50: 321-333.

[62] Adali, S. (2023). Variational principles for a double Rayleigh beam system undergoing vibrations and connected by a nonlinear Winkler-Pasternak elastic layer. Nonlinear Engineering, 2(1). https://doi.org/10.1515/nleng-2022-0259

[63] Ike, C.C. (2018). Flexural analysis of rectangular Kirchhoff plate on Winkler foundation using Galerkin-Vlasov variational method. Mathematical Modeling of Engineering Problems, 5(2): 83-92. https://doi.org/10.18280/mmep.050205

[64] Nwoji, C.U., Onah, H.N., Mama, B.O., Ike, C.C. (2018). Ritz variational method for bending of rectangular Kirchhoff plate under transverse hydrostatic load distribution. Mathematical Modeling of Engineering Problems, 5(1): 1-10. https://doi.org/10.18280/mmep.050101

[65] Ansari, R., Sahmani, A., Rouhi, H. (2011). Axial buckling analysis of single-walled carbon nanotubes in thermal environments via the Rayleigh–Ritz technique. Computational Materials Science, 50(10): 3050-3055. https://doi.org/10.1016/j.commatsci.2011.05.027

[66] Chakraverty, S., Behera, L. (2015). Free vibration of non-uniform nanobeams using Rayleigh–Ritz method. Physica E: Low-dimensional Systems and Nanostructures, 67: 38-46, https://doi.org/10.1016/j.physe.2014.10.039.

[67] Behera, L., Chakraverty, S. (2017). Static analysis of nanobeams using Rayleigh–Ritz method. Journal of Mechanics of Materials and Structures, 12(5): 603-616. https://doi:10.2140/jomms.2017.12.603

[68] Fakher, M., Hosseini-Hashem, S. (2017). Bending and free vibration analysis of nanobeams by differential and integral forms of nonlocal strain gradient with Rayleigh-Ritz method. Materials Research Express, 4(12): 125025. https://doi.org/10.1088/2053-1591/aa9dd6

[69] Eringen, A.C. (1983). On differential equations of nonlocal elasticity and solutions of screw dislocation and surface waves. Journal of Applied Physics, 54(9): 4703-4710. https://doi.org/10.1063/1.332803

[70] Eringen, A.C. (2002). Nonlocal continuum field theories, Springer, New York, 2002.

[71] Nicolescu, A.E., Bobe, A. (2021). Weak solution of longitudinal waves in carbon nanotubes. Continuum Mechanics and Thermodynamics, 33: 2065-2073. https://doi.org/10.1007/s00161-021-01001-7