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) .
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 [
x ¨ = f ( x ) + ε g ( x , x ˙ , t ) x ∈ R
Alex (2013) discussed the exact solution of a cubic-quintic Duffing oscillator of the form [
x ¨ + A x + 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
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
x ˙ ( 0 ) = 0 and x ( ± ∞ ) = 0
Given x ( 0 ) = a , x ˙ ( t ) = a y ( t ) , where a is an undetermined parameter, the Equation (1.1) becomes
y ¨ − y + a 2 y 3 = 0 (1.2)
Equation (1.2) satisfies the following boundary conditions
y ( 0 ) = 1 , y ˙ ( 0 ) = 0 , y ( ± ∞ ) = 0 (1.3)
Additionally, Talay assumed that the solution of Equation (1.2) can be expressed by index Series in the form
y ( t ) = ∑ c m e − m t = y 0 ( t ) + ∑ k = 1 ∞ y k (t)
According to the boundary conditions, he deduces
y 0 ( t ) = 2 e − t − e − 2 t
y 1 ( t ) = ℏ [ e − 2 t + γ 0 ( e − 3 t − 4 5 e − 4 t + 1 4 e − 5 t − 1 35 e − 6 t ) ] + c 1 e − t + c 2 e t = ℏ ( 5 64 e − t − e − 2 t + 35 36 e − 3 t − 7 4 e − 4 t + 35 64 e − 5 t − 1 16 e − 6 t )
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 [
Keqiang Li (2013) studied the periodic orbits of the modified Duffing differential equation [
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 [
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 ) .
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
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
c ¨ 1 ( t ) e t + c ¨ 2 ( t ) e − t + 2 ( c ˙ 1 ( t ) e t − 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
c ¨ ( t ) ( e t + e − t ) + 2 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 ) = m c ( 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 ) .
{ { c ¨ ( t ) + 2 c ˙ ( t ) e t − e − t e t + e − t + m c ( t ) = f ( t ) c 3 ( t ) ( e t + e − t ) 2 − m c ( t ) + f ( t ) = 0 (2.7)
Proof: Let z = c ( t ) ( e t + e − t ) , then Equation (2.6) leads to a cubic equation as
z 3 − m z + f ( t ) ( e t + e − t ) = 0 (2.8)
Its solutions can be expressed as
z = − f ( t ) ( e t + e − t ) + f 2 ( t ) ( e t + e − t ) − 4 m 3 27 2 3 + − f ( t ) ( e t + e − t ) − f 2 ( t ) ( e t + e − t ) 2 − 4 m 3 27 2 3 (2.9)
Specially, if m = 0 , we have the simplest representation as follows
z = − f ( t ) ( e t + e − t ) 3 .
Notice that Equation (2.7) has denominators in the form as ( e t + e − t ) ,so let
f ( t ) = − k 3 ( e t + e − t ) 4 (2.10)
Applying function (2.10) to equations set (2.7), we get k 3 = 8 k , it means k = 2 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 ) = − k 3 ( e t + e − t ) 4 , meanwhile, we have
z = k ( e t + e − t )
Let’s make a transformation: z = c ( t ) ( e t + e − t ) , then
c ( t ) = k ( e t + e − t ) 2
Applying c ( t ) and its first derivative and second derivative to equations set (2.7), we have
2 k ( e t + e − t ) 2 − 3 ( e t − e − t ) 2 ( e t + e − t ) 4 + 4 k e t − e − t ( e t + e − t ) 3 ⋅ e t − e − t e t + e − t = k 3 ( e t + e − t ) 4 (2.11)
Equation (2.11) follows k 3 = 8 k immediately. This finishes the proof of theorem1.
We have showed a few representative integer parameters m and functions in
In fact, m is not only integer, but it maybe any fractional as well. Some values of m and respective function are showed in
In this section, we shall give some examples to explain our constant variation method.
Example 1. Consider following Duffing equation
x ¨ − x + a x 3 = 0 (3.1)
1) a > 0 . Take time transformation τ = a t for Equation (3.1), we have
d 2 x d τ 2 − 1 a x + x 3 = 0
According to theorem 1, we can easily find a couple of m and f ( t ) (see
m | f ( t ) |
---|---|
0 | − 16 2 ( e t + e − t ) 4 |
1 | − 2 2 ⋅ 4 − ( e t − e − t ) 2 ( e t + e − t ) 4 |
2 | 4 2 ( e t − e − t ) 2 ( e t + e − t ) 4 |
3 | 2 2 ⋅ 4 + 3 ( e t − e − t ) 2 ( e t + e − t ) 4 |
4 | 8 2 ⋅ 2 + ( e t − e − t ) 2 ( e t + e − t ) 4 |
5 | 2 2 ⋅ 12 + 5 ( e t − e − t ) 2 ( e t + e − t ) 4 |
6 | 4 2 ⋅ 8 + 3 ( e t − e − t ) 2 ( e t + e − t ) 4 |
7 | 2 2 ⋅ 20 + 7 ( e t − e − t ) 2 ( e t + e − t ) 4 |
8 | 16 2 ⋅ 3 + ( e t − e − t ) 2 ( e t + e − t ) 4 |
9 | 2 2 ⋅ 28 + 9 ( e t − e − t ) 2 ( e t + e − t ) 4 |
10 | 4 2 ⋅ 16 + 5 ( e t − e − t ) 2 ( e t + e − t ) 4 |
11 | 2 2 ⋅ 36 + 11 ( e t − e − t ) 2 ( e t + e − t ) 4 |
12 | 8 2 ⋅ 10 + 3 ( e t − e − t ) 2 ( e t + e − t ) 4 |
m | f ( t ) |
---|---|
1 2 | − 2 ⋅ 12 − ( e t − e − t ) 2 ( e t + e − t ) 4 |
1 3 | − 2 2 3 ⋅ 20 − ( e t − e − t ) 2 ( e t + e − t ) 4 |
1 4 | − 2 2 2 ⋅ 14 − ( e t − e − t ) 2 ( e t + e − t ) 4 |
1 5 | − 2 2 5 ⋅ 36 − ( e t − e − t ) 2 ( e t + e − t ) 4 |
1 6 | − 2 3 ⋅ 44 − ( e t − e − t ) 2 ( e t + e − t ) 4 |
1 7 | − 2 2 7 ⋅ 52 − ( e t − e − t ) 2 ( e t + e − t ) 4 |
1 8 | − 2 4 ⋅ 60 − ( e t − e − t ) 2 ( e t + e − t ) 4 |
1 9 | − 2 2 9 ⋅ 68 − ( e t − e − t ) 2 ( e t + e − t ) 4 |
x = 2 2 a cosh τ a = 2 2 a cosh t , a > 0 (3.2)
2) a < 0 . Take time transformation τ = − a t for Equation (3.1), in same way we also get
x = 2 2 − a sinh t , a < 0 (3.3)
Example 2. Vasile (2011) discussed the unforced and undamped double-well Duffing equation [
x ¨ − ω 0 2 x + β x 3 = 0 (3.4)
where β 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 = 2 2 β cosh ( ω 0 t ) (3.5)
Example 3. Jaume and Ana (2013; 2015) discussed following modified Duffing differential equation [
y ¨ + a y − ε y 3 = ε h ( y , 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 = 2 2 − ε cosh ( − a t ) , ε < 0 (3.7)
y = 2 2 ε sinh ( − a t ) , ε > 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 [
x ¨ + A x + B x 3 + G x 5 = 0 (3.9)
Let’s discuss the special case G = 0 , A < 0 , B > 0 . Take time transformation τ = B t ( B > 0 ) for Equation (3.9), then (3.9) can be turned into
d 2 x d τ 2 + A B x + x 3 = 0
According the result in example 3, one can get the solution at once
x = 2 2 B cosh ( − A t ) , B > 0 and x = 2 2 − B sinh ( − A t ) , B < 0 (3.10)
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).
The authors declare no conflicts of interest regarding the publication of 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