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:

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.

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 :

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

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:

Equation (4), then, provides the relationship between the applied force at the tip and the tip-displacement as:

Our lumped parameter model is based on assuming that the density function is given by

For fixed time, integrating the equation over gives:

Equation (7) gives the following nonlinear spring-mass equation

The corresponding autonomeous equation of (8) in which is denoted by, is given by:

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

where is a positive parameter.

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:

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.

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,:

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

sought to have the following form

Substituting (13) in (10) and identifying the term of the same power of we obtain the following system of initial value problems:

We will show that the solutions of Equations (14) and (15) are given by:

where and. The solution of Equation (10) will then be given by:

• The solution of Equation (14) has the following form:

Substituting in (15) and writing as the Fourier series:

the following equation is obtained

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.

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:

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.

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.

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.