# Nonlinear Vibrational Analysis of Nanobeams Embedded in an Elastic Medium including Surface Stress Effects

Authors: S. Azizi; B. Safaei; A. M. Fattahi; M. Tekere
Journal: Advances in Materials Science and Engineering
 (2015)

---

## Abstract

Due to size-dependent behavior of nanostructures, the classical continuum models are not applicable for the analyses at this submicron size. Surface stress effect is one of the most important matters which make the nanoscale structures have different properties compared to the conventional structures due to high surface to volume ratio. In the present study, nonlinear free vibrational characteristics of embedded nanobeams are investigated including surface stress effects. To this end, a thin surface layer is assumed on the upper and lower surfaces of the cross section to separate the surface and bulk of nanobeams with their own different material properties. Based on harmonic balance method, closed-form analytical solution is conducted for nonlinear vibrations to obtain natural frequencies of embedded nanobeams with and without considerations of surface elasticity and residual surface tension effects corresponding to the various values of nondimensional amplitude, elastic foundation modulus, and geometrical variables of the system. Selected numerical results are given to indicate the influence of each one in detail.


---

## Body

## 1. Introduction

Because of advancing science and technology, there has been quite request for new materials with extraordinary properties for implementation in various emerging fields of development. Due to a wide range of applications nanostructures and materials at nanoscale have been attracting consequential interest in different aspects of nanoscience and nanotechnology.For example, Xu and his assistants [1] explored the PH-sensitive property of the single-walled carbon nanotubes (SWCNTs) based on electroactive group on SWCNT by different pulse voltammetry technique to show the broad biosensor application of carbon nanotube in electrochemical area. Wu et al. [2] investigated the resonant frequency shift of the fixed-free SWCNTs caused by the addition of a nanoparticle to the beam tip using a beam-bending model to explore the suitability of the SWCNT as a mass detector device. A nanocomposite film of embedded SWCNTs was prepared by Penza et al. [3] to develop highly sensitive microacoustic vapor sensors based on surface acoustic waves. They indicated that SWCNT nanocomposite vapor sensitivity results in being significantly enhanced with respect to traditional organic molecular cavities materials.The surface stress effect can be easily observed at the atomic scale due to high surface to volume ratio, and this has been clearly indicated and explained [4, 5]. The main reason for the phenomenon is based on different environment conditions for the atoms which are near free surface compared to the ones at the bulk of the material. For conventional structures in which the number of atoms near free surface is so negligible relative to the total number of atoms of the system, stress effects on surface physics are unimportant and can be neglected. However, in the case of nanostructures, through increasing ratio of surface area and volume these effects can make considerable influence.There are so many computational and experimental researches in which the influence of surface stress effects on behavior of the structures at nanoscale has been studied. Gurtin et al. [6] analyzed the effect of residual surface tension on the stress concentration around a nanosized spheroidal cavity in an isotropic medium based on surface elasticity. Radi [7] described the material behavior with the indeterminate theory of couple surface elasticity to investigate the problem of a stationary semi-infinite crack in an elastic solid. The effect of surface/interface stress on the macroscopic plastic behaviors of nanocomposites has been investigated by Zhang et al. [8] through a new second-order moment nonlinear micromechanics theory.Gurtin et al. [6] presented a model for vibrational analysis of microbeams including surface elasticity influences represented by applying distributed traction and compressive axial force. Then, for an isotropic surface, Gurtin and Ian Murdoch [9] proposed a continuum model of surface elasticity in which the Laplace-Young equation extended to solid materials. The model has the capability to incorporate the surface stress effects into the mechanical response of nanostructures which has been applied in many studies conducted for various problems about mechanical behavior of the structures at nanoscale [10–18].A new frontier of research in the area of computational nanomechanics is to study the behavior of structures at very small length scales. Then Ansari and Sahmani [19] proposed a nonclassical solution to analyze bending and buckling responses of nanobeams including surface stress effects. Explicit formulas are proposed relevant to each type of beam theory to evaluate the surface stress effects on the displacement profile and critical buckling load of the nanobeams. Ansari and Sahmani [20] presented the free vibration characteristics of nanoplates including surface stress effects are investigated based on the continuum modeling approach. To this end, Gurtin-Ian Murdoch continuum elasticity approach is incorporated into the different types of plate theory, namely, classical plate theory (CLPT) and first-order shear deformation theory (FSDT) to develop nonclassical continuum plate models for free vibration analysis of the nanoplates including surface stress effects.Karličić et al. [21] proposed the stochastic parametric vibrations of a Voigt-Kelvin nanobeam based on Eringen’s nonlocal elasticity theory of the Helmholtz and bi-Helmholtz type of kernel. The axial force consists of a constant part and a time-dependent stochastic function. By using the direct Lyapunov method, bounds of almost sure asymptotic stability of a viscoelastic nanobeam are obtained as a function of retardation time, variance of the stochastic force, geometric ratio, scale coefficients, and intensity of the deterministic component of axial loading.In this work, the surface stress effects including both surface elasticity and residual surface tension are incorporated into the nonlinear free vibrations of nanobeams embedded in an elastic medium based on the model of Gurtin et al. [6]. Selected numerical results are presented to indicate the influence of nondimensional amplitude, elastic foundation modulus, and geometrical variables of the system on the value of surface stress effects.


## 2. Surface Elasticity Model for Nanobeams

Consider a straight uniform beam with the lengthL and rectangular cross section of thickness 2
h. A coordinate system (
x
,
y, and z
) is introduced on the central axis of the beam, whereas x-axis is taken along the length of the beam, y-axis is taken in the width direction, and z-axis is taken along the depth (height) direction. Also, the origin of the coordinate system is selected at the left end of the beam (Figure 1). It is assumed that the deformations of the beam take place in the x
-
z plane, so the displacement components (
u
1
,
u
2, and u
3
) along the axes (
x
,
y, and z
) are only dependent on the x and z coordinates and time t.Figure 1
Geometry and properties of the considered problem.To represent the surface stress effects in this work, thin surface layers with surface elastic modulusE
s are assumed on the upper and lower surfaces of the cross section to separate the surfaces and the bulk of nanobeams. According to the Laplace-Young equation, the upper (
τ
u
) and lower (
τ
b
) residual surface tension can be considered as distributed loading on surfaces [6] which can be expressed in the following forms:(1a)
p
u
x
=
τ
u
b
∂
2
w
∂
x
2
,
(1b)
p
b
x
=
τ
b
b
∂
2
w
∂
x
2
,where w is the transverse displacement of the beam, b denotes the width of nanobeam, and τ
u and τ
b stand for the residual surface stresses on the upper and lower surfaces of nanobeam, respectively.Furthermore, based on the Winkler elastic foundation model [22] the normal pressure of elastic medium is simulated as a series of closely spaced, vertical linear elastic springs. Thereupon, the loading corresponding to this type of foundation can be expressed as(2)
p
e
x
=
-
K
w
,where K is the Winkler modulus parameter corresponding to the normal pressure.On the basis of energy method, the total potential energy of the system contains the elastic strain energy of the bulk and surfaces. Through implementing generalized Hamilton’s principle, the governing differential equation for nonlinear free vibrations of Euler-Bernoulli nanobeam embedded in an elastic medium is obtained as(3)
E
I
+
2
E
s
b
h
2
∂
4
w
x
,
t
∂
x
4
-
τ
u
+
τ
b
b
∂
2
w
x
,
t
∂
x
2
+
K
w
x
,
t
+
ρ
A
∂
2
w
x
,
t
∂
t
2
=
E
A
2
L
∫
0
L
∂
w
x
,
t
∂
x
2
d
x
∂
2
w
x
,
t
∂
x
2
,where E and E
s are the bulk and surface Young’s modulus, respectively.The transverse displacementw can be considered in the following generalized form which satisfies the simply supported boundary conditions:(4)
w
x
,
t
=
W
m
t
sin
⁡
m
π
x
L
.Substituting (4) into the governing equation of (3) yields(5)
d
2
W
m
t
d
t
2
+
m
π
L
4
E
I
ρ
A
+
m
π
L
4
2
E
s
b
h
2
ρ
A
+
m
π
L
2
τ
u
+
τ
b
b
ρ
A
+
K
ρ
A
W
m
t
+
1
4
E
ρ
m
π
L
4
W
m
3
t
=
0in which m is the vibrational mode number, ρ denotes the mass density, and A and I are cross-sectional area and moment of inertia of the nanobeam, respectively.The transformations are considered as follows:(6)
r
=
L
A
,
X
m
=
W
m
r
,
τ
=
τ
u
+
τ
b
,
μ
=
ω
t
,
ω
L
=
m
π
L
2
E
I
ρ
A
,
ω
K
=
K
ρ
A
,
ω
E
=
m
π
L
2
2
E
s
b
h
2
ρ
A
,
ω
τ
=
m
π
L
τ
u
+
τ
b
b
ρ
A
.Now (5) can be expressed in the following dimensionless nonlinear equation: (7)
ω
ω
l
2
X
¨
m
+
1
+
ω
E
ω
l
2
+
ω
τ
ω
l
2
+
ω
k
ω
l
2
X
m
+
α
X
m
3
=
0
,in which α
=
1
/
4.The harmonic balance method is utilized to solve the nonlinear equations given by (7). Based on this method, it is assumed that the solution can be represented as a truncated Fourier series of the form(8)
X
m
μ
=
A
m
cos
⁡
μ
+
A
m
+
n
cos
⁡
3
μ
.Applying the assumed solution to (7) in conjunction with expanding each term as a Fourier series leads to throwing out higher harmonics other than the ones included in the original assumed solution. After equating the coefficients like cosine harmonics (such as the cos
⁡
τ and cos
⁡
3
τ terms), one can arrive at a set of algebraic equations of this form (9)
-
ω
ω
l
2
A
m
+
1
+
ω
E
ω
l
2
+
ω
τ
ω
l
2
+
ω
k
ω
l
2
A
m
+
3
4
α
A
m
2
A
m
+
n
+
3
2
α
A
m
+
n
2
A
m
+
3
4
α
A
m
3
=
0
,
-
9
ω
ω
l
2
A
m
+
n
+
1
+
ω
E
ω
l
2
+
ω
τ
ω
l
2
+
ω
k
ω
l
2
A
m
+
n
+
3
2
α
A
m
2
A
m
+
n
+
3
4
α
A
m
+
n
3
+
1
4
α
A
m
3
=
0
.The preceding nonlinear equations can be solved for ω, A
m
+
n using the Newton-Raphson numerical scheme by choosing A
m as a determined parameter.


## 3. Numerical Results and Discussion

In this section, selected numerical results are presented to show the influence of different parameters on the nonlinear vibrational response of nanobeams embedded in an elastic medium including surface stress effects. Anodic alumina with two surface crystal orientations of1
0
0 and 1
1
1 is considered in the analytical solution. The material properties and related surface parameters for each crystal orientation are given in Table 1.Table 1
Mechanical properties of anodic alumina.
Crystal orientation
Young’s modulus (GPa)
Poisson’s ratio
Mass density (Kg/m3)
Surface elastic modulus (N/m)
Residual surface tension (N/m)
[1 0 0]
70
0.3
2700
-
7.9253
0.5689
[1 1 1]
70
0.3
2700
5.1882
0.9108


### 3.1. Effect of Elastic Foundation Modulus on Nonlinear Frequency Ratio

To explore the effect of surrounding elastic medium on the nonlinear vibrations of nanobeams, variation ofω
nonlinear  with  surface  stress  effects
/
ω
linear  with  surface  stress  effects with elastic foundation modulus for 1
0
0 and 1
1
1 crystal orientations is illustrated in Figures 2(a) and 2(b), respectively. In both cases it is assumed that L
=
30 nm and h
=
b
=
0.1
L. It can be observed that, for both crystal orientations, the difference between linear and nonlinear responses of nanobeams increases with increase of nondimensional amplitude of vibration which is more prominent for lower values of elastic foundation modulus. In other words, the elastic medium makes the effect of nonlinearity on the vibrational behavior of nanobeams less important.Figure 2
Effect of elastic foundation modulus on nonlinear frequency response for (a)1
0
0 crystal orientation, (b) 1
1
1 crystal orientation.
(a)
(b)


### 3.2. Effect of Length of Nanobeam on Nonlinear Frequency Ratio

The effect of length of nanobeams on the nonlinear free vibrational analysis is shown in Figure3. In this figure, the value of ω
nonlinear  with  surface  stress  effects
/
ω
linear  with  surface  stress  effects is plotted against the length of the embedded nanobeam corresponding to different elastic foundation modulus and crystal orientations. It can be seen that, for both crystal orientations and all values of elastic foundation moduli, the increase of length of nanobeams leads to decrease of the effect of nonlinearity inclusion. Also, it is shown that, for the higher value of elastic foundation moduli, the behaviors of 1
0
0 and 1
1
1 nanobeams are approximately the same. However, for the lower value of K, the difference between the responses of two nanobeams with different crystal orientations is considerable especially for lower values of length.Figure 3
Effect of length of the nanobeams on nonlinear frequency response.Additionally, for the case of1
1
1 crystal orientation, it is seen that, for the higher elastic medium (K
=
10
8), by increasing the length of nanobeam, the effect of surface stress decreases with lower slope, especially for shorter nanobeams.


### 3.3. Effect of Height of Nanobeam on the Frequency Ratio

To see the effect of height of nanobeams on the nonlinear vibrational response, the variation ofω
nonlinear  with  surface  stress  effects
/
ω
nonlinear  without  surface  stress  effects with respect to the height of embedded nanobeams relevant to various elastic foundation moduli is depicted in Figure 4 for both 1
0
0 and 1
1
1 crystal orientations.Figure 4
Effect of height of nanobeams on nonlinear frequency response corresponding to different values of elastic foundation modulus: (a)1
0
0 crystal orientation, (b) 1
1
1 crystal orientation.
(a)
(b)It is observed that, in the case having no elastic medium(
K
=
0
), the value of frequency ratio increases for all values of height. However, with the inclusion of elastic medium, the behavior of nanobeam changes totally as the height of nanobeam can make increase or decrease of the frequency ratio at different values. Furthermore, in the case of 1
0
0 crystal orientation, a critical value for the height of nanobeam is found corresponding to which the surface stress effect will be the same for all values of elastic foundation modulus.Also, for both crystal orientations, it is observed that, forK
=
10
7, by increasing the height of nanobeam, the effect of surface stress increases first and then decreases. However, for a stiffer elastic medium (K
=
10
8), the surface stress effect decreases continuously by increasing the height of nanobeam. In other words, for a medium value of elastic medium, in order to decrease the effect of surface stress, it is necessary to increase the height of nanobeam above a special value. However, for high elastic medium, any increase in the value of beam height leads to lower surface stress effect.Figure5 indicates the influence of height of nanobeams on the value of ω
nonlinear  with  surface  stress  effects
/
ω
nonlinear  without  surface  stress  effects relevant to different amounts of nondimensional amplitude of vibrations for 1
0
0 and 1
1
1 crystal orientations. It can be found that, in both types of crystal orientation, the height of nanobeam decreases the frequency ratio until a certain value of height. Beyond this certain height the value of nonlinear frequency ratio increases. It is seen that this pattern of variation of the frequency ratio is more sensible corresponding to higher values of nondimensional amplitude.Figure 5
Effect of height of the nanobeams on nonlinear frequency response corresponding to different values of nondimensional amplitude: (a)1
0
0 crystal orientation, (b) 1
1
1 crystal orientation.
(a)
(b)Also, it can be observed that, for the case of1
0
0 crystal orientation, the values of frequency ratio corresponding to a certain value of height are equal for all nondimensional amplitudes. In other words, there is an amount of the height for the nanobeam with 1
0
0 crystal orientation relevant to which the surface stress effect is the same for all values of nondimensional amplitude.


## 3.1. Effect of Elastic Foundation Modulus on Nonlinear Frequency Ratio

To explore the effect of surrounding elastic medium on the nonlinear vibrations of nanobeams, variation ofω
nonlinear  with  surface  stress  effects
/
ω
linear  with  surface  stress  effects with elastic foundation modulus for 1
0
0 and 1
1
1 crystal orientations is illustrated in Figures 2(a) and 2(b), respectively. In both cases it is assumed that L
=
30 nm and h
=
b
=
0.1
L. It can be observed that, for both crystal orientations, the difference between linear and nonlinear responses of nanobeams increases with increase of nondimensional amplitude of vibration which is more prominent for lower values of elastic foundation modulus. In other words, the elastic medium makes the effect of nonlinearity on the vibrational behavior of nanobeams less important.Figure 2
Effect of elastic foundation modulus on nonlinear frequency response for (a)1
0
0 crystal orientation, (b) 1
1
1 crystal orientation.
(a)
(b)


## 3.2. Effect of Length of Nanobeam on Nonlinear Frequency Ratio

The effect of length of nanobeams on the nonlinear free vibrational analysis is shown in Figure3. In this figure, the value of ω
nonlinear  with  surface  stress  effects
/
ω
linear  with  surface  stress  effects is plotted against the length of the embedded nanobeam corresponding to different elastic foundation modulus and crystal orientations. It can be seen that, for both crystal orientations and all values of elastic foundation moduli, the increase of length of nanobeams leads to decrease of the effect of nonlinearity inclusion. Also, it is shown that, for the higher value of elastic foundation moduli, the behaviors of 1
0
0 and 1
1
1 nanobeams are approximately the same. However, for the lower value of K, the difference between the responses of two nanobeams with different crystal orientations is considerable especially for lower values of length.Figure 3
Effect of length of the nanobeams on nonlinear frequency response.Additionally, for the case of1
1
1 crystal orientation, it is seen that, for the higher elastic medium (K
=
10
8), by increasing the length of nanobeam, the effect of surface stress decreases with lower slope, especially for shorter nanobeams.


## 3.3. Effect of Height of Nanobeam on the Frequency Ratio

To see the effect of height of nanobeams on the nonlinear vibrational response, the variation ofω
nonlinear  with  surface  stress  effects
/
ω
nonlinear  without  surface  stress  effects with respect to the height of embedded nanobeams relevant to various elastic foundation moduli is depicted in Figure 4 for both 1
0
0 and 1
1
1 crystal orientations.Figure 4
Effect of height of nanobeams on nonlinear frequency response corresponding to different values of elastic foundation modulus: (a)1
0
0 crystal orientation, (b) 1
1
1 crystal orientation.
(a)
(b)It is observed that, in the case having no elastic medium(
K
=
0
), the value of frequency ratio increases for all values of height. However, with the inclusion of elastic medium, the behavior of nanobeam changes totally as the height of nanobeam can make increase or decrease of the frequency ratio at different values. Furthermore, in the case of 1
0
0 crystal orientation, a critical value for the height of nanobeam is found corresponding to which the surface stress effect will be the same for all values of elastic foundation modulus.Also, for both crystal orientations, it is observed that, forK
=
10
7, by increasing the height of nanobeam, the effect of surface stress increases first and then decreases. However, for a stiffer elastic medium (K
=
10
8), the surface stress effect decreases continuously by increasing the height of nanobeam. In other words, for a medium value of elastic medium, in order to decrease the effect of surface stress, it is necessary to increase the height of nanobeam above a special value. However, for high elastic medium, any increase in the value of beam height leads to lower surface stress effect.Figure5 indicates the influence of height of nanobeams on the value of ω
nonlinear  with  surface  stress  effects
/
ω
nonlinear  without  surface  stress  effects relevant to different amounts of nondimensional amplitude of vibrations for 1
0
0 and 1
1
1 crystal orientations. It can be found that, in both types of crystal orientation, the height of nanobeam decreases the frequency ratio until a certain value of height. Beyond this certain height the value of nonlinear frequency ratio increases. It is seen that this pattern of variation of the frequency ratio is more sensible corresponding to higher values of nondimensional amplitude.Figure 5
Effect of height of the nanobeams on nonlinear frequency response corresponding to different values of nondimensional amplitude: (a)1
0
0 crystal orientation, (b) 1
1
1 crystal orientation.
(a)
(b)Also, it can be observed that, for the case of1
0
0 crystal orientation, the values of frequency ratio corresponding to a certain value of height are equal for all nondimensional amplitudes. In other words, there is an amount of the height for the nanobeam with 1
0
0 crystal orientation relevant to which the surface stress effect is the same for all values of nondimensional amplitude.


## 4. Conclusion

In this work, nonlinear vibrational characteristics of embedded nanobeams in an elastic medium are investigated including surface stress effects. To this end, both surface elasticity and residual surface tension effects are incorporated into the nonlinear vibrational analysis of embedded Euler-Bernoulli nanobeam.Analytical solution is presented for simply supported nanobeams and selected numerical results are given to show the influences of nondimensional amplitude, elastic foundation modulus, and geometrical variables of the system in detail corresponding to two types of crystal orientations. It is found that surface stress effect can have significant influence on the nonlinear vibrational behavior of nanobeams and this effect varies relevant to different parameters of the system.



---

*Source: 318539-2015-12-28.xml*