ABSTRACT

A lot of methods, such as Jacobian elliptic function analysis, are used to look for the explicit exact solution of Duffing differential equation. The key of the analysis is to construct quotient trigonometric function, and then nonlinear algebraic equation set theory and method are used for the solution of some kinds of nonlinear Duffing differential equation. In this paper, the exact solution of Duffing equation is obtained by using constant variation method, making use of the formula to solve cubic equations and general solution of the homogeneous equation of Duffing equation with appropriate Constant m and function $f(t)$ .

Keywords:

Duffing Equation, Exact Solution, Constant Variation Method

1. Introduction

The Duffing equation is a non-linear second-order differential equation. The equation describes the motion of a damped oscillator with a more complicated potential than in simple harmonic motion; in physical terms, it models, for example, a spring pendulum whose spring’s stiffness does not exactly obey Hooke’s law. It is also an example of a dynamical system that exhibits chaotic behavior.

Many of the dynamic behavior of mechanical power model can be described by a single degree of freedom oscillator [1] [2] [3] [4] [5]. The systems can be described by the following equation

$\ddot{x}=f(x)+\epsilon g(x,\stackrel{\dot{}}{x},t)x\in R$

Alex (2013) discussed the exact solution of a cubic-quintic Duffing oscillator of the form [6]

$\ddot{x}+Ax+B{x}^3+G{x}^5=0$

Alex’s idea is both simple and enlightening, he turn looking for exact solution into solving algebraic equations set. Without lose of generally, in this paper, we consider the following Duffing equation without external force and damping

$\ddot{x}-x+x^3=0$ (1.1)

Talay’s solve approach is much complex. He assumes that the Equation (1.1) has boundary conditions in the form

$\stackrel{\dot{}}{x}(0)=0$ and $x(\pm \infty )=0$

Given $x(0)=a$ , $\stackrel{\dot{}}{x}(t)=ay(t)$ , where a is an undetermined parameter, the Equation (1.1) becomes

$\ddot{y}-y+a^2{y}^3=0$ (1.2)

Equation (1.2) satisfies the following boundary conditions

$y(0)=1$ , $\stackrel{\dot{}}{y}(0)=0$ , $y(\pm \infty )=0$ (1.3)

Additionally, Talay assumed that the solution of Equation (1.2) can be expressed by index Series in the form

$y(t)={\displaystyle \sum c_m}{e}^{-mt}=y_0(t)+{\displaystyle \underset{k=1}{\overset{\infty }{\sum }}{y}_k}(t)$

According to the boundary conditions, he deduces

$y_0(t)=2e^{-t}-e^{-2t}$

$\begin{array}{l}{y}_1(t)=\hslash [e^{-2t}+{\gamma }_0(e^{-3t}-\frac{4}{5}{e}^{-4t}+\frac{1}{4}{e}^{-5t}-\frac{1}{35}{e}^{-6t})]+c_1{e}^{-t}+c_2{e}^t\\ =\hslash (\frac{5}{64}{e}^{-t}-e^{-2t}+\frac{35}{36}{e}^{-3t}-\frac{7}{4}{e}^{-4t}+\frac{35}{64}{e}^{-5t}-\frac{1}{16}{e}^{-6t})\end{array}$

In the last few years, many scientists work hard to look for exact solutions and to develop some efficient procedures. However, what the key of problem is, the different type of nonlinear differential equation can only be efficiently solved by different corresponding methods. For instant, the Jacobi elliptic function expansion method [7], the exp-function method [8], the F-expansion method [9], the tanh method [10], the auxiliary equation method [11], the simplest equation method [12], and linearized harmonic balance method [13] are reliable methods for obtaining exact solutions of specific types of nonlinear differential equations.

Keqiang Li (2013) studied the periodic orbits of the modified Duffing differential equation [14]. He paid his attention to the existence of multiple periodic solutions of asymptotically linear Duffing equations with resonance on the left side of the first eigenvalue using index theory and the three-critical-points theorem and obtained a new result. As the function in the right side of his equation is zero, how to look for the exact solution of the equation is a significative problem. Alex (2013) used Jacobi elliptic functions to derive the exact solution of the cubic-quintic Duffing oscillator [6]. He also proved that the exact angular frequency of this cubic-quintic Duffing equation is given in terms of the complete elliptic integral of the first kind. Vasile (2011) has presented a method for obtaining exact solutions for Duffing Equations [15]. Marinca (2011) has presented a method for obtaining exact solutions for Duffing and double-well Duffing equations [16]. What he used is traditional quotient trigonometric function expansion method, and his solution procedure is similar to the exp-function method. However, his result in third section has a distance to reach perfect.

Exact analytical solutions for a family of autonomous ordinary differential equations arising in connection with many important applied problems such as a cubic-quintic Duffing oscillator, Helmholtz?Duffing oscillator, and nonlinear Schrödinger equation are considered. Recently, some theory and applications are found in many contexts, such as, Jiang et al. (2012) presented a new approach via the well-known Poincaré-Birkhoff theorem to obtain the existence of periodic solutions to impulsive problems, they considered an impulsive Duffing equation, and then found the possibility of applying a generalized form of the Poincaré-Birkhoff theorem due to Ding to construct infinitely many periodic solutions of the impulsive Duffing equation even in a resonance case [17]. Jaume and Ana (2015) studied the periodic solutions of a kind of non-autonomous Duffing differential equation with forced pendulum. And, they obtained the condition for that the non-autonomous Duffing differential system has a periodic solution such that it tends to the periodic solution when forced pendulum satisfied with some conditions [18]. Gholam and Emmanuel (2015) have provided a set of inequalities that helps in finding transformation to obtain solutions to other Duffing equations from a known set of solutions to a related equation in a systematic manner [19]. Alexander (2018) has obtained the sufficient and necessary conditions for the existence of a positive periodic solution to a kind of non-autonomous Duffing type equation [20]. David and Douglas (2018) perform a complete analysis of the effect on the asymptotic dynamics of the insertion of delays into a sort of unforced Duffing type equation [21].

One of the most exciting recent advances of nonlinear science has been the development of some methods to look for exact solutions of nonlinear differential equations. It is important that many mathematical models are described by nonlinear differential equations and exact solutions are preferable instead of approximate ones. The aim of this paper is to put forward constant variation method to derive the exact solution of a Duffing equation of the form (1.1). The exact solution of Duffing equation is obtained by solving cubic equations and looking for appropriate constant m and function $f(t)$ .

2. Constant Variation Method

In general, the Duffing equation does not admit an exact symbolic solution. However, many approximate methods work well:

• Expansion in a Fourier series will provide an equation of motion to arbitrary precision.

• The cubic term also called the Duffing term can be approximated as small and the system treated as a perturbed simple harmonic oscillator.

• The Jacobi elliptic function expansion method yields a complicated but workable solution.

• Any of the various numeric methods such as Euler’s method and Runge-Kutta can be used.

In the special case of the undamped and unforced, an exact solution of Duffing equation has been obtained using our constant variation method

Now, let’s begin at the solution of the homogeneous linear equations. The homogeneous linear equations corresponding to Equation (1.1) is

$\ddot{x}-x=0$ (2.1)

Its general solution is

$x=c_1{e}^t+c_2{e}^{-t}$ (2.2)

Assume we have functions $c_1(t)$ and $c_2(t)$ , such that the solution of Equation (1.1) has the form as

$x=c_1(t)e^t+c_2(t)e^{-t}$ (2.3)

Taking Equation (2.3) into Equation (1.1), we obtain following differential equation

${\ddot{c}}_1(t)e^t+{\ddot{c}}_2(t)e^{-t}+2({\stackrel{\dot{}}{c}}_1(t)e^t-{\stackrel{\dot{}}{c}}_2(t)e^{-t})+{(c_1(t)e^t+c_2(t)e^{-t})}^3=0$ (2.4)

Because functions $c_1(t)$ and $c_2(t)$ are discretionary, hence, we let $c_1(t)=c_2(t)=c(t)$ .

Then, Equation (2.4) becomes

$\ddot{c}(t)(e^t+e^{-t})+2\stackrel{\dot{}}{c}(t)(e^t-e^{-t})+c^3(t){(e^t+e^{-t})}^3=0$ (2.5)

Therefore, we obtain a more compact differential equation form for unknown function $c(t)$ .

If we can find an undetermined constant m and an undetermined function $f(t)$ , such that they satisfy with following equation

$f(t)=mc(t)-c^3(t){(e^t+e^{-t})}^2$ (2.6)

Then, it is possible for us to reduce solving differential Equation (2.5) to looking for a couple of undetermined constant m and undetermined function $f(t)$ . Our main conclusions as follows.

Theorem 1. The Equation (2.5) has solutions if only and if the following equations set have solutions for given constant m and undetermined function $f(t)$ .

$\{\begin{array}{l}\{\ddot{c}(t)+2\stackrel{\dot{}}{c}(t)\frac{e^t-e^{-t}}{e^t+e^{-t}}+mc(t)=f(t)\\ c^3(t){(}^{e^t}-mc(t)+f(t)=0\end{array}$ (2.7)

Proof: Let $z=c(t)(e^t+e^{-t})$ , then Equation (2.6) leads to a cubic equation as

$z^3-mz+f(t)(e^t+e^{-t})=0$ (2.8)

Its solutions can be expressed as

$\begin{array}{l}z=\sqrt[3]{\frac{-f(t)(e^t+e^{-t})+\sqrt{f^2(t)(e^t+e^{-t})-\frac{4m^3}{27}}}{2}}\\ +\sqrt[3]{\frac{-f(t)(e^t+e^{-t})-\sqrt{f^2(t){(e^t+e^{-t})}^2-\frac{4m^3}{27}}}{2}}\end{array}$ (2.9)

Specially, if $m=0$ , we have the simplest representation as follows

$z=-\sqrt[3]{f(t)(e^t+e^{-t})}$ .

Notice that Equation (2.7) has denominators in the form as $(e^t+e^{-t})$ ,so let

$f(t)=\frac{-k^3}{{(e^t+e^{-t})}^4}$ (2.10)

Applying function (2.10) to equations set (2.7), we get $k^3=8k$ , it means $k=2\sqrt{2}$ or $k=0$ .

On the other hand, we can prove that $k^3$ is forever equal to 8k when $m=0$ .

Actually, if $m=0$ , then $f(t)=\frac{-k^3}{{(e^t+e^{-t})}^4}$ , meanwhile, we have

$z=\frac{k}{(e^t+e^{-t})}$

Let’s make a transformation: $z=c(t)(e^t+e^{-t})$ , then

$c(t)=\frac{k}{{(e^t+e^{-t})}^2}$

Applying $c(t)$ and its first derivative and second derivative to equations set (2.7), we have

$2k\frac{{(e^t+e^{-t})}^2-3{(e^t-e^{-t})}^2}{{(e^t+e^{-t})}^4}+4k\frac{e^t-e^{-t}}{{(e^t+e^{-t})}^3}\cdot \frac{e^t-e^{-t}}{e^t+e^{-t}}=\frac{k^3}{{(e^t+e^{-t})}^4}$ (2.11)

Equation (2.11) follows $k^3=8k$ immediately. This finishes the proof of theorem1.

We have showed a few representative integer parameters m and functions in Table 1.

In fact, m is not only integer, but it maybe any fractional as well. Some values of m and respective function are showed in Table 2.

3. Applications

In this section, we shall give some examples to explain our constant variation method.

Example 1. Consider following Duffing equation

$\ddot{x}-x+a{x}^3=0$ (3.1)

1) $a>0$ . Take time transformation $\tau =\sqrt{a}t$ for Equation (3.1), we have

$\frac{d^{2}x}{d{\tau }^2}-\frac{1}{a}x+x^3=0$

According to theorem 1, we can easily find a couple of m and $f(t)$ (see Table 2), for instant, take m = 1/2 and $f(\tau )=-\sqrt{2}\cdot \frac{12-{(e^{\frac{1}{\sqrt{a}}\tau }-e^{-\frac{1}{\sqrt{a}}\tau })}^2}{{(e^{\frac{1}{\sqrt{a}}\tau }+e^{-\frac{1}{\sqrt{a}}\tau })}^4}$ . In the light of theorem 1, one can get the solution of Equation (3.1) at once

$x=\frac{2\sqrt{2}}{\sqrt{a}\cosh \frac{\tau }{\sqrt{a}}}=\frac{2\sqrt{2}}{\sqrt{a}\cosh t},a>0$ (3.2)

2) $a<0$ . Take time transformation $\tau =\sqrt{-a}t$ for Equation (3.1), in same way we also get

$x=\frac{2\sqrt{2}}{\sqrt{-a}\sinh t},a<0$ (3.3)

Example 2. Vasile (2011) discussed the unforced and undamped double-well Duffing equation [16]

$\ddot{x}-{\omega }_0^{2}x+\beta x^3=0$ (3.4)

where $\beta$ is a positive constant which may not be a small value.

We can easily obtain exact solution of Equation (3.4) in the same way as example 1.

$x=\frac{2\sqrt{2}}{\sqrt{\beta }\cosh ({\omega }_{0}t)}$ (3.5)

Example 3. Jaume and Ana (2013; 2015) discussed following modified Duffing differential equation [15] [19]

$\ddot{y}+ay-\epsilon y^3=\epsilon h(y,\stackrel{\dot{}}{y})$ (3.6)

In the special case $h=0,a<0$ , we obtain one of the solutions of Equation (3.6) at once in the same way as example1.

$y=\frac{2\sqrt{2}}{\sqrt{-\epsilon }\cosh (\sqrt{-a}t)},\epsilon <0$ (3.7)

$y=\frac{2\sqrt{2}}{\sqrt{\epsilon }\sinh (\sqrt{-a}t)},\epsilon >0$ (3.8)

Example 4. By using Jacobian elliptic function, Alex (2013) derive the exact solution of a cubic-quintic Duffing oscillator of the form [14]

$\ddot{x}+Ax+B{x}^3+G{x}^5=0$ (3.9)

Let’s discuss the special case $G=0,A<0,B>0$ . Take time transformation $\tau =\sqrt{B}t(B>0)$ for Equation (3.9), then (3.9) can be turned into

$\frac{d^{2}x}{d{\tau }^2}+\frac{A}{B}x+x^3=0$

According the result in example 3, one can get the solution at once

$x=\frac{2\sqrt{2}}{\sqrt{B}\cosh (\sqrt{-A}t)},B>0$ and $x=\frac{2\sqrt{2}}{\sqrt{-B}\sinh (\sqrt{-A}t)},B<0$ (3.10)

Acknowledgements

This research was supported by National Natural Science Foundation of China (71171128) and Shanghai University’s horizontal project (Project name: company convertible bond pricing method and its application).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Nelson, O.N.F., Yu, Z., Dorian, B.P. and Wang, Y.Y. (2018) A New Method for the Exact Solution of Duffing Equation. Journal of Applied Mathematics and Physics, 6, 2718-2726. https://doi.org/10.4236/jamp.2018.612225

References