1. Introduction
The graphene-based resonator and its application to mass sensing based on nonlinear waves have been poorly studied numerically [1] . Some researchers use discrete atomic or Monte Carlo approach for numerical simu- lation and some use local or nonlocal continuum mechanics approaches, however, their models are based on linear material constitutive equation for graphene ([1] [2] ). It is, however, well-known that graphene behaves nonlinearly even for small strains and there is no obvious yield point or a linear portion on it’s stress-strain curve. In fact, it is proved experimentally and theoretically in [3] that the mechanical behaviour of a single layer of graphene sheet can be accuartely modeled by a continuum nonlinear constitutive equation ([4] -[6] ). This consti- tutive equation in it’s one dimensional form is:
(1)
where
is the axial strain,
the axial stress,
the Young’s modulus,
the third-order
elastic stiffness constant, and
the ultimate yield stress of the graphene. It appears that recent studies in literature have not incorporated the constant
into their models for the vibration analysis of graphene layers. The main objective of this work is to model and understand how graphene behaves in free and forced axial vibrations and to calculate the nonlinear resonance frequencies based on Equation (1). To initiate this study, a simplified nonlinear spring model is derived based on the lumped parameter method. We show that the third- order elastic stiffness
plays an important role in modeling the patterns of graphene in axial vibration. Within a range of the initial energy, we show that there exist periodic solutions similar to the ones obtained using the corresponding linear models and that the free oscillations are nearly sinusoidal. However, as the initial energy approaches a threshold level, the limiting free oscillations deviate drastically from the sinusoidal oscil- lations predicted by linear models. Our initial results provide some quantitative regimes in which a grap-hene resonator can operate near harmonic and non-harmonic motions. The initial results of this project provide some insight information and data on the patterns of axial vibration of a graphene monolayer which can be useful for design of graphene-based resonators. By extending this simple nonlinear spring-mass model to more realistic models, it is possible to provide new design guide to help make more efficient resonators and wave guides, shorten the design cycle and provide more accurate assessment of the mechanical behavior of these devices. In Section 2, we derive the nonlinear spring lumped parameter model from the nonlinear wave equation of a graphnene sheet under axial vibration; in Section 3, we study the existence of periodic solutions by using phase plane analysis and perturbation techniques; in Section 4, we compute the approximate analytical solutions of free vibrations using the two-scales splitting method and obtain the associated natural frequencies and ampli- tudes and compare to numerical results; in Section 5, we compute numerical solutions of forced vibrations and obtain frequency sweeps.
2. The Nonlinear Lumped Parameter Model
A graphene sheet with uniform cross-section in axial vibration with fixed-free ends can be modeled by sub- stituting (1) into the standard balance of momentum equation
to obtain the following nonlinear wave equation subject to initial and boundary conditions :
(2)
Here, we use
for second order time derivative of
and
for spatial derivative of
. The core- sponding steady state problem with a concentrated load of magnitude
at the tip
is
(3)
where
is the Dirac delta function. Assuming that
, the exact solution of (3) can be found by integrating (3) and applying the boundary conditions. First, we integrate from 0 to
and then from
to
and using
, for
we get:
(4)
Equation (4), then, provides the relationship between the applied force
at the tip and the tip-displacement
as:
(5)
Our lumped parameter model is based on assuming that the density function is given by
(6)
For fixed time
, integrating the equation
over
gives:
(7)
Equation (7) gives the following nonlinear spring-mass equation
(8)
The corresponding autonomeous equation of (8) in which
is denoted by
, is given by:
(9)
is the lumped-mass at the tip of the sheet,
is the first order stiffness and
is the thrid-order stiffness
constant in (1). Using the change of variable
we obtain the equivalent non-dimensional equation
(10)
where
is a positive parameter.
3. Existence of Periodic Solutions of Free Vibration
We will show that for given initial conditions
and
, Equation (10) has periodic solutions for certain range of
. To determine the ranges of
for which existence of periodic solutions occur, we examine the phase diagrams associated with the Equation (10) defined by:
(11)
where
.
We make the following observations:
1) The
-intercepts associated with (11) are
,
2)
represents the energy at initial position
,
3) when
, the phase diagrams are the circles
with center
and radius
.
We prove that for any
, there exists
such that for
Equation (10) has a periodic solution and for
there exists no periodic solutions of (10). Since periodic solutions of (10) correspond to closed curves of the phase diagram, we need to examine the x-intercepts of (11) and their dependence on the equation parameters. The x-intercepts of (11) are the zeros of:
![]()
which is an even function. Therefore it is enough to consider
, for
. Some properties of
are:
;
for
and
for
, since
; and
. Based on these properties we can distinguish the following three cases (corresponding to Figures 1(a)-(c)):
Case (a): ![]()
No periodic solution.
Case (b): ![]()
Only one
-intercept and there is a periodic solution.
Case (c):
Two
-intercepts and there is a periodic solution.
We conclude that the bifurcation point for a given
, is
, see Figure 2.
Furthermore, for each
when periodic solution exists, we determine the frequency and the period numerically (see Figure 3).
It is demonstrated in Figure 2 that at a lower energy level
the free vibration is approximately sinu- soidal, however at a higher level of energy
the free vibration deviate drastically from the sinusoidal pattern which has not been captured by previous models that do not include the third order elastic constant
. When periodic solutions exist, Figure 3 indicates that our model shows that at each fixed energy level, the frequency and period of a given graphene sheet depend nonlinearly on the parameter
which depend on the material elastic constants as well as the length of the sheet.
4. Double Scales Analytical Approximations of Free Vibration
Multiple scales method is often used to solve nonlinear equations with small parameters in nonlinear vibrations. Double scales are used herein to find an approximate solution of the first order for Equation (10). The solution is then compared to results obtained by numerical integration using Matlab. The new time scales are
,
where
represents the fast time and
represents the slow time. The derivative with respect to
will be written as function of the derivative with respect to
,
:
(12)
Instead of determining the solution as a function of
, we determine it as a function of
and
. To this end, we change the independent variable in Equation (10) from
to
and
. A solution of the equation is
(a)
(b)
Figure 3. (a) Frequency diagram; (b) Period diagram.
sought to have the following form
(13)
Substituting (13) in (10) and identifying the term of the same power of
we obtain the following system of initial value problems:
(14)
(15)
We will show that the solutions of Equations (14) and (15) are given by:
(16)
where
and
. The solution of Equation (10) will then be given by:
(17)
• The solution of Equation (14) has the following form:
(18)
Substituting
in (15) and writing
as the Fourier series:
![]()
the following equation is obtained
(19)
To avoid unbounded solutions, we set the secular terms of the
-Equation (19), containing
and
to zero. This gives the system:
![]()
whose solution gives
and
. The solution for
is then obtained by returning to
the original time using
.
To find
, we need to solve the linear differential equation:
![]()
and obtain (16), where
and
are determined easily from the initial conditions. Using
instead of
and
instead of
and replacing
and
by its values, the expression in Equation (17) is verified.
Remark
The solution for
shows an odd multiple of the
, this can be seen clearly in the expression for
. These frequencies are the harmonics of the main mode or frequency. It is a typicall feature of nonlinear differential equations that the harmonics are related directly to the nonlinear terms. Our expressions are verified numerically by calculating the solution in the frequency domain using the fourier transform and comparing with the analytical results. The results in time and frequency domains are shown below in Figure 4 and Figure 5.
5. Nonlinear Vibration under Harmonic Excitation
In this section we characterize the nonlinear spring Equation (10) by a harmonic excitation and studying the system’s nonliear responses. The equation of motion is given by:
(20)
where
,
is a small real number called detuning parameter. The frequency
of the excitation is near the resonnance of the coresponding linear frequency of the equation
. We present the numerical solutions in time and frequency domains and demonstrate the use of the frequency sweep method in detecting the nonlinear resonnance of the system.We solve Equation (20) using Matlab solver to obtain numerical results in the time domain. FFT algorithm is then applied to the time signal to find the frequencies of the solutions. The expected frequency corresponds to the excitation
, the nonlinear resonance and some harmonics. The double scales method can be used to find analytical approximate solution of Equation (20) similar to the autonomeous system case of Section 4. We present our numerical results in Figure 6.
We use the frequency sweep method to detect nonlinear resonnance of the nonlinear system by direct intergration. The method begins by defining a grid of frequencies around the linear resonnace and intergrate the
system at each point of frequency. The maximum displacement of the solution is then plotted against the fre- quency mesh. The curve shows a peak corresponding to the nonlinear resonnance of the system. The numerical results show the dependence of the nonlinear frequency on the magnitude of excitation and on the parameter
. Figure 7 and Figure 8 show the numerical results for some values of the system parameters.
6. Conclusion
A simplistic nonlinear spring model is derived from the axial wave equation of a graphene sheet based on the
quadratic constitutive stress-strain equation. Using phase plane analysis, existence of periodic wave solutions
and bifurcation points depending on the parameter
are verified for free vibrations. Perturbation
method of time scales depending on
is used to study axial vibrations subject to harmonic excitation. The results are compared with the corresponding linear spring model which does not include the third order elastic
constant
. It is demonstrated through our analysis and numerical solutions that the bifurcation
parameter
critically affects the solutions quantitatively and numerically, therefore we conclude that the third order elastic constant
in the continuum mechanics based modeling of graphene should be included in further study of the dynamic behavior if higher accuracy of solutions are desired. In future studies we plan to examine the axial vibrations corresponding to the full model (2) numerically using finite differences, finite element and numerical bifurcation techniques. In addition, we plan to examine the vertical vibrations using nonlinear beam and plate equations.
Acknowledgements
The paper’s first co-author acknowledges the funding provided by the NPRP grant 08-777-1-141 from the Qatar National Research Fund (a member of Qatar Foundation) to Prof. Prabir Daripa of Texas A& M University at College Station, TX 77842, USA while working on this project.