\documentclass{jande}
\usepackage{jande}
\usepackage{rotating} 
\Volume[2012]{8}{2}
\Page{1}         

%%%%% to generate high-quality PDF
\usepackage[T1]{fontenc}
\usepackage{ae,aecompl}
%%%%% to generate high-quality PDF


\theoremstyle{definition}  
\theoremstyle{remark}      
\theoremstyle{plain}      
\newtheorem{theorem}{Theorem}%[section]
\newtheorem{proposition}[theorem]{Proposition}
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{coro}[theorem]{Corollary}
\newtheorem{conj}[theorem]{Conjecture}


\def\proof{{\bf Proof. }}


\theoremstyle{remark}
\newtheorem{remark}{Remark}%[section]

\theoremstyle{definition}
\newtheorem{example}{Example}%[section]
\newtheorem{definition}{Definition}%[section]



\titlemark{ Author 1, Author 2/Journal of Applied Nonlinear Dynamics Vol-number (year) page1--page2 }
\authormark{Author 1, Author 2/Journal of Applied Nonlinear Dynamics Vol-number (year) page1--page2}


\begin{document}

\title{\first{Title of Article }}


\setcounter{footnote}{1}

\author{\noindent\large Author 1\footnote{Corresponding author.\\\hspace*{5.5mm}Email address:
*@* }~, Author 2}




\address{\normalsize Department of Mechanical and Industrial Engineering, Southern Illinois University Edwardsville,\\[-2mm]
Edwardsville, IL62026-1805, USA}


\abstract{
\begin{table}[h!]
\vspace*{-5mm}
\doublerulesep 0.05pt
\tabcolsep 7.8mm
\vspace*{2mm}
\setlength{\tabcolsep}{7.5pt}
\hspace*{-2.5mm}\begin{tabular*}{16.5cm}{r|||||l}
%\hline
\multicolumn{2}{l}{\rule[-6pt]{16.5cm}{.01pt}}\\
\parbox[t]{6cm}{\small
\vspace*{.5mm}
\hfill {\bf Submission Info}\par
\vspace*{2mm}
\hfill Communicated by Referees\par
\hfill Received DAY MON YEAR \par
\hfill Accepted DAY MON YEAR\par
\hfill Available online DAY MON YEAR\par
\noindent\rule[-2pt]{6.3cm}{.1pt}\par
\vspace*{2mm}
\hfill {\bf Keywords}\par
\vspace*{2mm}
\hfill Keyword 1\par
\hfill Keyword 2\par
\hfill Keyword 3\par
\hfill Keyword 4}
&
\parbox[t]{9.85cm}{
\vspace*{.5mm}
{\normalsize\bf Abstract}\par
\renewcommand{\baselinestretch}{.8}
\normalsize \vspace*{2mm} {\small Each article should be preceded by
an abstract (8-11 lines long) that summarizes the content. The
abstract will appear online at www.****.com and be available with
unrestricted access. This allows unregistered users to read the
abstract as a teaser for the complete article.Each article should be preceded by
an abstract (8-11 lines long) that summarizes the content. The
abstract will appear online at www.****.com and be available with
unrestricted access. This allows unregistered users to read the
abstract as a teaser for the complete article. }
\par
\par
\par
\hfill{\scriptsize\copyright 2012 L\&H Scientific Publishing, LLC. All rights reserved.}}\\[-4mm]
&\\
\multicolumn{2}{l}{\rule[15pt]{16.5cm}{.01pt}}\\\end{tabular*}
\vspace*{-7mm}
\end{table}
}


\maketitle
\thispagestyle{first}


\renewcommand{\baselinestretch}{1}
\normalsize

\section{Section Heading}

\noindent Please note that the first line of text that follows a heading is not indented,
whereas the first lines of all subsequent paragraphs are.

Use the standard \verb|equation| environment to typeset your equations, e.g.
\begin{equation}
\ddot{x}+\delta\dot{x}-\alpha x+\beta x^3=Q_0\cos \Omega t
\end{equation}
where $\dot{{x}}=dx/dt$ is velocity. $Q_{0} $ and $\Omega $ are
excitation amplitude and frequency, respectively. $\delta $is
damping coefficient. $\alpha $ and $\beta$ are linear and nonlinear
stiffness coefficients of the Duffing oscillator. however, for
multiline equations we recommend to use the \verb|eqnarray|
environment\footnote{In physics texts please activate the class
option \texttt{vecphys} to depict your vectors in \textbf{\itshape
boldface-italic} type - as is customary for a wide range of physical
subjects.}.
\begin{eqnarray}
a \times b = c \nonumber\\
\vec{a} \cdot \vec{b}=\vec{c}
\label{eq:01}
\end{eqnarray}


\begin{theorem}
theorem
\end{theorem}



\begin{proposition}
theorem
\end{proposition}

\begin{lemma}
theorem
\end{lemma}

\begin{coro}
theorem
\end{coro}

\begin{conj}
theorem
\end{conj}

In mechanical engineering, in 1918, Duffing \cite{1} presented the hardening
spring model to describe the vibration of electro-magnetized vibrating beam.
Since then, the Duffing oscillator has been extensively used to describe
nonlinear structural vibrations in structural dynamics. In 1964, Hayashi \cite{2}
discussed the approximate periodic solutions and the corresponding stability
by the averaging method and harmonic balance method. In 1973, Nayfeh \cite{3}
used the perturbation method to approximate periodic motion of the Duffing
oscillators (also see, Nayfeh and Mook \cite{4}). In 1979, Holmes \cite{5} showed the
strange attractors of chaotic motions in nonlinear oscillators via the
Duffing oscillator with a twin-well potential. In 1980, Ueda \cite{6} used
numerical simulations to show chaotic motion via period-doubling of periodic
motions of Duffing oscillators. In 1997, Luo and Han \cite{7} analytically
presented the stability and bifurcation conditions of periodic motions of
the Duffing oscillator. The constant term of the analytical solution for the
steady-state motion of the Duffing oscillator was not considered. In 1996,
Luo and Han \cite{8} presented an improved solution of the Duffing oscillator
with a twin-well potential. For analytical prediction of chaos, in 1999, Luo
and Han \cite{9} investigated chaotic motions in nonlinear rod through the
Duffing oscillator. For the periodically forced Duffing oscillator with
damping, the analytical prediction of periodic solutions is still very
difficult. In this paper, the analytical solutions of periodic motions will
be investigated and the analytical route of periodic motions to chaos will
be of great interest.

\begin{figure}[h!]
\centering
\includegraphics{fig1.eps}
%\vspace*{-30mm}
\caption{The analytical prediction of periodic solutions based on two
harmonic terms (HB3): (a) constant term $a_0$; (b)-(d)
harmonic amplitudes $A_k (k=1,2,3)$; and (e)-(f)
harmonic phases $\varphi_k$ ($k=1,2)$ for right potential well.
($\delta=0.5,\alpha=-10.0,\beta=10.0,Q_0=10.0)$.}
\label{fig1}
\end{figure}




To look for approximate analytical solution of nonlinear oscillator, such an
issue started from Lagrange \cite{10} to investigate the three-body problem as a
perturbation of the two-body problem by the method of averaging. In the end
of the 19th century, Poincare \cite{11} further developed the perturbation theory
to investigate the motions of celestial bodies. In 1920, van der Pol \cite{12}
used the method of averaging to determine the periodic solutions of
oscillation systems in circuits. Until 1928, the asymptotic validity of the
method of averaging was not proved. In 1928, Fatou \cite{13} gave the proof of
the asymptotic validity through the solution existence theorems of
differential equations. In 1935, Krylov, Bogoliubov and Mitropolsky \cite{14}
further developed the method of averaging, and the detailed presentation was
given. In 1964, Hayashi \cite{2} presented the perturbation methods including
averaging method and principle of harmonic balance. In 1969, Barkham and
Soudack \cite{15} extended the Krylov-Bogoliubov method for the approximate
solutions of nonlinear autonomous second order differential equations (also
see, Barkham and Soudack \cite{16}). In 1987, a generalized harmonic balance
approach was used by Garcia-Margallo and Bejarano \cite{17} to determine
approximate solutions of nonlinear oscillations with strong nonlinearity. In
the same year, Rand and Armbruster \cite{18} used the perturbation method and
bifurcation theory to determine the stability of periodic solutions. In
1989, Yuste and Bejarano \cite{19} used the elliptic functions rather than
trigonometric functions to improve the Krylov-Bogoliubov method. In 1990,
Coppola and Rand \cite{20} used the averaging method with elliptic functions to
determine approximation of limit cycle. In 1997, Luo and Han \cite{7}
analytically studied the stability and bifurcations of periodic solutions of
Duffing oscillators through the first order harmonic balance method, and
provided the analytical conditions for the Hopf and saddle-node
bifurcations. To obtain accurate results of periodic solutions in nonlinear
vibration, many harmonic terms are included in the harmonic balance method.
In 2008, Peng et al \cite{21} presented the approximate period-1 solution for the
Duffing oscillator by the HB3 method compared with the fourth-order
Runge-Kutta method. In 2011, Luo and Huang \cite{22} further discussed a
generalized harmonic balance method to obtain the analytical solution of
period-1 motion. Luo and Huang \cite{23} also presented a generalized harmonic
balance method to determine period-$m$ solutions in nonlinear oscillators.



In this paper, the generalized harmonic balance method will be used to
investigate analytical periodic motions in the periodically forced Duffing
oscillator with a twin-well potential. The bifurcation tree from period-1
motions to chaos will be presented with varying parameters. The
corresponding unstable periodic motions in the Duffing oscillator will be
presented for a better understanding of nonlinear dynamics in such a Duffing
oscillator. Numerical illustrations of stable and unstable periodic motions
will be carried out.

\section{Section Heading}

From Eq.(1), the standard form is
\setcounter{equation}{2}
\begin{equation}
\ddot{x}+f(x,\dot{x},t)=0
\end{equation}

\subsection{Subsection Heading}
\noindent
The Fourier series expression of any periodic motion in nonlinear systems
needs infinite terms to give the exact solution of such a periodic motion.
In practice, it is impossible to do so. Thus, the truncated Fourier series
solutions will be used to give an approximate solution that can be close to
the exact solution. From such approximate, analytical solutions, the
equilibrium solution of coefficient dynamical system for the Fourier series
of the periodic motion can be obtained from Eq.(3) using Newton-Raphson
method, and the stability and bifurcation analysis of the such equilibrium
points can be completed through the eigenvalue analysis. The system
parameters are
\begin{equation}
\delta=0.5,\alpha=-10.0,\beta=10,Q_0=10.0
\end{equation}
The backbone curves of harmonic amplitude varying with excitation frequency
$\Omega$ are illustrated. The harmonic amplitude and phase are
defined by
\begin{equation}
A_{k/m}\equiv\sqrt{b^2_{k/m}+c^2_{k/m}},\varphi_{km}=\arctan\dfrac{c_{k/m}}{b_{k/m}}
\end{equation}
and the corresponding solution in Eq.(43) is
\begin{equation}
x^*(t)=a^{(m)}_0+\sum^N_{k=1}A_{k/m}\cos\left(\dfrac{k}{m}\Omega t-\varphi_{k/m}\right).
\end{equation}

In Luo and Han \cite{8}, one term harmonic term was considered for period-1
motions for the large and small orbit. In this paper, many harmonic terms
will be considered to achieve a more accurate prediction of the periodic
motions. For period-1 motion, the first three harmonic terms of the Fourier
series expansion (HB3) will be used to obtain the approximate periodic
solutions. The constant term $a^{(1)}_0\equiv a_0$ and the first three
harmonic amplitudes $A_k$ and phases $\varphi_k(k=1,2,3)$ versus excitation frequency are plotted in
Fig.1(a)-(g), respectively. A parameter map is presented in Fig.2.

The initial conditions for stable period-1 motion ($\Omega=2.75)$,
unstable period-1 motion and stable period-2 motion ($\Omega=2.753)$,
unstable period-1 motion, unstable period-2 motion and stable period-4
motion ($\Omega=2.7537)$ are listed in Table 1.

\begin{table}[h!]
\doublerulesep 0.1pt
\tabcolsep 7.8mm
\centering
\caption{\rm Input data for numerical simulations of periodic motions $(\delta=0.5,\alpha=-10.0,\beta=10.0,Q_0=10.0)$}
\vspace*{2mm}
\renewcommand{\arraystretch}{1.3}%%�����и�
\setlength{\tabcolsep}{22pt}
\footnotesize{\begin{tabular*}{16.5cm}{cccccc}
\hline\hline\hline
&  \raisebox{-2ex}[0pt][0pt]{$\Omega$}  &  \multicolumn{2}{c}{Initial conditions $(t=0.0)$}  &  \raisebox{-2ex}[0pt][0pt]{Stability}  &  \raisebox{-2ex}[0pt][0pt]{Period-$m$}\\\cline{3-4}
& & $x_0$  &  $\dot{x}_0$\\\hline
Fig.10(a)  &  2.75  &  $-0.724512$  &  0.251206  &  Stable(HB)  &  Period-1(HB5)\\
Fig.10(b)  &  2.75  &  $-0.724512$  &  0.251206  &  Stable(HB)  &  Period-1(HB5)\\
Fig.10(c)  &  2.75  &  $-0.724512$  &  0.251206  &  Stable(HB)  &  Period-1(HB5)\\
Fig.10(d)  &  2.75  &  $-0.724512$  &  0.251206  &  Stable(HB)  &  Period-1(HB5)\\
Fig.10(e)  &  2.75  &  $-0.724512$  &  0.251206  &  Stable(HB)  &  Period-1(HB5)\\
Fig.10(f)  &  2.75  &  $-0.724512$  &  0.251206  &  Stable(HB)  &  Period-1(HB5)
\\\hline\hline\hline\end{tabular*}
}
\renewcommand{\arraystretch}{1}
\end{table}



\begin{figure}[h!]
\centering
\includegraphics{fig2.eps}
\caption{A parameter map from the analytical prediction of periodic
solutions based on three harmonic terms (HB3): (a) Global view and (b)
zoomed view. ($\delta=0.5,\alpha=-10.0,\beta=10,Q_0=10.0)$.}
\label{fig1}
\end{figure}



\section{Conclusions}

\noindent
In this paper, analytical routines of period-1 motions to chaos in the
Duffing oscillator with a twin-potential well were discussed comprehensively
through the generalized harmonic balance method. The analytical solutions of
period-$m$ motions were developed by the Fourier series and the corresponding
Hopf bifurcations of periodic motions cause new periodic motions with
period-doubling. Three analytical routes of asymmetric period-1 motions to
chaos were developed. The approximate, analytical periodic solutions were
verified via numerical simulations, and the analytical, unstable periodic
motions were given as well. With exact unstable periodic motion, the
numerical simulations should stay with the analytical solution if without
any computational errors. The analytical routes with unstable periodic
motions can lead us to find unstable chaos.



\begin{thebibliography}{999}
\bibitem{1}Duffing, G. (1918), Erzwunge Schweingungen bei veranderlicher eigenfrequenz. Braunschweig: F. Viewig u. Sohn.
\bibitem{2}Hayashi, G. (1964), \textit{Nonlinear oscillations in Physical Systems, }McGraw-Hill Book Company: New York
\bibitem{3}Nayfeh, A.H. (1973), \textit{Perturbation Methods}, John Wiley: New York.
\bibitem{4}Nayfeh, A.H. and Mook, D.T. (1979), \textit{Nonlinear Oscillation}, John Wiley: New York.
\bibitem{5}Holmes, P.J. (1979), A nonlinear oscillator with strange attractor, \textit{Philosophical Transactions of the Royal Society} A, \textbf{292}, 419-448.
\bibitem{6}Ueda, Y. (1980), Explosion of strange attractors exhibited by the Duffing equations, \textit{Annuals of the New York Academy of Science}, \textbf{357}, 422-434.

\bibitem{7}
Luo, A.C.J. and Han, R.P.S.\textsc{ (1997), }A quantitative stability and bifurcation analyses of a generalized Duffing oscillator with strong nonlinearity, \textit{Journal of Franklin Institute}, \textbf{334B}, 447-459.

\bibitem{8}
Han, R.P.S. and Luo, A.C.J. (1996), Comments on ``Subharmonic resonances and criteria for escape and chaos in a driven oscillator'', \textit{Journal of Sound and Vibration}, \textbf{196}(2), 237-242.

\bibitem{9}
Luo, A.C.J. and Han, R.P.S. (1999), Analytical predictions of chaos in a nonlinear rod, \textit{Journal of Sound and Vibration}, \textbf{227}(2), 523-544.

\bibitem{10}
Lagrange, J.L. (1788), \textit{Mecanique Analytique} (2 vol.)(edition Albert Balnchard: Paris, 1965).

\bibitem{11}
Poincare, H. (1899), \textit{Methodes Nouvelles de la Mecanique Celeste}, Vol.3, Gauthier-Villars: Paris.

\bibitem{12}
van der Pol, B. (1920), A theory of the amplitude of free and forced triode vibrations, \textit{Radio Review}, \textbf{1,} 701-710, 754-762.

\bibitem{13}
Fatou, P. (1928), Sur le mouvement d'un systeme soumis `a des forces a courte periode, \textit{Bull. Soc. Math. }\textbf{56}, 98--139

\bibitem{14}
Krylov, N.M. and Bogolyubov, N.N. (1935), \textit{Methodes approchees de la mecanique non-lineaire dans leurs application a l'Aeetude de la perturbation des mouvements periodiques de divers phenomenes de resonance s'y rapportant}, Academie des Sciences d'Ukraine: Kiev. (in French).

\bibitem{15}
Barkham, P.G.D. and Soudack, A.C. (1969), An extension to the method of Krylov and Bogoliubov, \textit{International Journal of Control}, \textbf{10}, 377-392.

\bibitem{16}
Barkham, P.G.D. and Soudack, A.C.(1970), Approximate solutions of nonlinear, non-autonomous second-order differential equations, \textit{International Journal of Control}, \textbf{11}, 763-767.

\bibitem{17}
Rand, R.H. and Armbruster, D. (1987), \textit{Perturbation Methods, Bifurcation Theory, and Computer Algebra}. Applied Mathematical Sciences, no. \textbf{65}, Springer-Verlag: New York.

\bibitem{18}
Garcia-Margallo J. and Bejarano, J.D. (1987), A generalization of the method of harmonic balance, \textit{Journal of Sound and Vibration}, \textbf{116}, 591-595.

\bibitem{19}
Yuste, S.B. and Bejarano, J.D. (1989), Extension and improvement to the Krylov-Bogoliubov method that use elliptic functions, \textit{International Journal of Control}, \textbf{49}, 1127-1141.

\bibitem{20}
Coppola, V.T. and Rand, R.H. (1990), Averaging using elliptic functions: Approximation of limit cycle, \textit{Acta Mechanica}, \textbf{81}, 125-142.

\bibitem{21}
Peng, Z.K., Lang, Z.Q., Billings, S.A. and Tomlinson, G.R. (2008), Comparisons between harmonic balance and nonlinear output frequency response function in nonlinear system analysis, \textit{Journal of Sound and Vibration}, \textbf{311}, 56-73.

\bibitem{22}
Luo, A.C.J. and Huang, J.Z. (2011), Approximate solutions of periodic motions in nonlinear systems via a generalized harmonic balance, \textit{Journal of Vibration and Control}, \underline {in press}.

\bibitem{23}
Luo, A.C.J. and Huang, J.Z. (2011), Analytical dynamics of period-$m$ flows and chaos in nonlinear systems, \textit{International Journal of Bifurcation and Chaos}, \underline {in press}.

\bibitem{24} Luo, A.C.J. (2011), \textit{Regularity and Complexity in Dynamical Systems}, Springer: New York.
\end{thebibliography}

\end{document}