Next Article in Journal
The Fuzzy u-Chart for Sustainable Manufacturing in the Vietnam Textile Dyeing Industry
Previous Article in Journal
A Case Study on Iteratively Assessing and Enhancing Wearable User Interface Prototypes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lie Symmetry Classification of the Generalized Nonlinear Beam Equation

1
Department of Mathematics, East China University of Science and Technology, Shanghai 200237, China
2
School of Data Science and Engineering, East China Normal University, Shanghai 200062, China
*
Author to whom correspondence should be addressed.
Symmetry 2017, 9(7), 115; https://doi.org/10.3390/sym9070115
Submission received: 22 March 2017 / Revised: 25 June 2017 / Accepted: 6 July 2017 / Published: 11 July 2017

Abstract

:
In this paper we make a Lie symmetry analysis of a generalized nonlinear beam equation with both second-order and fourth-order wave terms, which is extended from the classical beam equation arising in the historical events of travelling wave behavior in the Golden Gate Bridge in San Francisco. We perform a complete Lie symmetry group classification by using the equivalence transformation group theory for the equation under consideration. Lie symmetry reductions of a nonlinear beam-like equation which are singled out from the classification results are investigated. Some classes of exact solutions, including solitary wave solutions, triangular periodic wave solutions and rational solutions of the nonlinear beam-like equations are constructed by means of the reductions and symbolic computation.

1. Introduction

In this paper, we study the group properties of higher-order nonlinear wave-type equations. As a basic model, we consider the fourth-order generalized nonlinear beam equation (GNBE) or nonlinear wave equation of the form
u t t = - [ K ( u ) u x x x ] x + [ D ( u ) u x ] x + F ( u ) ,
where K = K ( u ) , D = D ( u ) and F = F ( u ) are arbitrary smooth functions ( hereafter the subscripts t and x denote differentiation with respect to these variables). Equation (1) generalizes a wide range of the known nonlinear wave equations arising in applications. The case with K identically zero and D ( u ) 0 for almost all u, is the case of second-order nonlinear wave equation which has already been studied in many practical contexts including shallow water waves theory, dynamics of a finite nonlinear string and elastic-plastic materials, etc. (see Refs. [1], pp. 50–52 and [2]). Hereinafter we assume that K ( u ) is not identically zero, so that the governing equation is of the fourth order, including a fourth-order wave term when K is non-negative.
The simplest equation of the form Equation (1), with K ( u ) = 1 , D ( u ) = 0 and F ( u ) = 0 , is the one dimensional linear beam equation or the fourth-order linear wave equation [3,4]
u t t = - u x x x x .
when K ( u ) = 1 , D ( u ) = 0 and F ( u ) = - k u + + 1 , u + = max { u , 0 } , Equation (1) is reduced to the following classical nonlinear beam equation or the fourth-order nonlinear wave equation of the form [5]
u t t = - u x x x x - k u + + 1 ,
which is an normalization of the original beam equations
u t t = - u x x x x - k u + + W ( x ) + ε f ( t , x )
arising in the historical events of travelling wave behavior in the Golden Gate Bridge in San Francisco [6] and has been proposed as a model for a suspension bridge. When k ( u ) = 1 , D ( u ) = 0 , Equation (1) become
u t t = - u x x x x + F ( u ) ,
which is a slight generalization of the classical nonlinear beam Equation (3) and is also presented by McKenna and Walter [7] in studying travelling wave solutions.
Up to now, the mathematic structure and properties of the fourth-order beam Equations (2)–(4) have been widely investigated. For instance, the mathematical state of the art concerning the operator and semigroup theory for Equation (2) on a bounded x-interval has been studied in [4]. In [5,8], the existence of travelling wave solutions and standing wave oscillations of the nonlinear beam Equation (3) have been considered. Chen and McKenna gave a variational proof of the existence of travelling wave solutions for Equation (4) via the Mountain Pass Lemma and showed the existence of at least one nontrivial solution for the equation under consideration in [9,10]. In particular, for F ( u ) = e u - 1 - 1 , they did also obtain numerical solutions by applying the Mountain Pass algorithm to a finite subinterval of R. They claimed that solutions seem to exist in the range 0 < c < 2 . As the wave speed approaches 2 , the solutions became highly oscillatory in nature, whereas when c approaches 0, they appear to go to infinity in amplitude. In [11,12,13], Equation (4) has been studied numerically by using either continuation methods or variational numerical methods in order to gain more information on the structure of the equation solutions set. In [14] Humphreys and McKenna considered the existence of multiple periodic solutions for the nonlinear beam Equation (4), while in [15,16] Doole and Hogan transformed the beam Equation (4) into an ordinary differential equation and treated it by using dynamical systems method. Recently, there are also some researches devoted to symmetry group structure and exact solutions of the fourth-order nonlinear beam Equation (4) [17,18].
However, the generalized nonlinear beam Equation (1) has been much less extensively investigated. It was only a few researches [19] that were devoted to qualitative mathematic properties such as quasi-periodic solutions of Equation (1) with K ( u ) = D ( u ) = 1 , while the symmetry group properties and corresponding algebraic structure as well as explicit exact solutions of Equation (1) and its invariant models still remain open. Therefore, the aim of the present work is to find all possible Lie symmetries, which Equation (1) can admit depending on the function triplets ( K , D , F ) , i.e., to solve the so-called group classification problem, which was formulated and solved for a class of nonlinear heat equations in the pioneering work by Ovsiannikov in 1959 [20] and now is the core stone of modern group analysis [21,22]. This problem for the second-order wave equation was probably first solved by Barone et al. in [23] and subsequently was extended to other general forms by many authors in the last two decades [2,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43], but were all limited to second-order cases.
Ovsiannikov’s method (also referred as the Lie-Ovsiannikov method) of Lie symmetry classification of differential equations [21] is based on the classical Lie scheme and a set of equivalence transformations of a given equation [44]. The formal application of this method to equations containing several arbitrary functions (Equation (1) contains three arbitrary functions) usually leads to a large number of equations admitting nontrivial Lie algebras of invariance [44].
Recent years, this method has been extended by many authors, in which they proposed a numbers of novel techniques, such as algebraic methods based on subgroup analysis of the equivalence group [45,46,47,48] and their generalizations [49,50,51,52,53,54], local transformations and form-preserving transformations [44,55,56,57,58], to solve group classification problem for numerous nonlinear partial differential equations. In this paper we extend the classical Lie-Ovsiannikov method based on equivalence transformations to the generalized nonlinear beam equation. We first carry out group classification of Equation (1) under the usual equivalence group. Then similar reductions of the classification models are performed and invariant solutions are also constructed. It is found that some similarity solutions are solutions with physical interest, including solitary wave solutions, triangular periodic wave solutions and rational solution.
The rest of paper is organized as follows. In Section 2, we derive the equivalence group and perform the group classification related to Equation (1), i.e., all possible Lie symmetries, which this equation can admit depending on the form the functions K , D and F, are found. In Section 3, the symmetry reductions and some exact solutions are constructed for particular case of Equation (1) that are likely to be useful in applications. Section 4 some concluding remarks are reported.

2. Lie Symmetry Classification

Background and procedures of the modern Lie group theory are well described in literature [21,22,35,37,50,59]. Without going into the details of the theory, we present only the results below.
Let
v = τ ( t , x , u ) t + ξ ( t , x , u ) x + ϕ ( t , x , u ) u .
be a vector field or infinitesimal operator on the space of independent and dependent variables t , x , u . A local group of transformations G is a symmetry group of Equation (1) if and only if
pr ( 4 ) v | ( ) = 0 ,
whenever = u t t + [ K ( u ) u x x x ] x - [ D ( u ) u x ] x - F ( u ) = 0 for every generator of G, where pr ( 4 ) v is the fourth-order prolongation of v (Here we do not write the explicit prolongation formulas to avoid tediousness, one can see the book [59] or the paper [54] for details). Expanding Equation (6) we get
ϕ t t + [ K ( u ) u x u x x x + K ( u ) u x x x x - D ( u ) u x 2 - D ( u ) u x x - F ( u ) ] ϕ + [ K ( u ) u x x x - 2 D ( u ) u x ] ϕ x - D ( u ) ϕ x x + K ( u ) u x ϕ x x x + K ( u ) ϕ x x x x = 0 ,
which must be satisfied whenever Equation (1) is satisfied. Substituting the formulae of ϕ t t , ϕ x , ϕ x x , ϕ x x x and ϕ x x x x into Equation (1) we get an equation of t , x , u and the derivatives of τ , ξ , ϕ , u . Replacing u t t by the right side of Equation (1) whenever it occurs, and equating the coefficients of the various independent monomials in the partial derivatives of u to zero, we obtain the determining equations
τ x = τ u = ξ t = ξ u = ϕ u u = 0 , 2 ( 2 ϕ x u - 3 ξ x x ) K + ϕ x K u = 0 , 2 ( τ t - 2 ξ x ) K + ϕ K u = 0 , ( ϕ x u - ξ x x ) K u = 0 , ( 4 ϕ x x x u - ξ x x x x ) K + ϕ x x x K u + ( ξ x x - 2 ϕ x u ) D - 2 ϕ x D u = 0 , 2 ϕ t u - τ t t = 0 , 2 ( 3 ϕ x x u - 2 ξ x x x ) K + 2 ( ξ x - τ t ) D - ϕ D u = 0 , ϕ t t + ϕ x x x x K - ϕ x x D + ( ϕ u - 2 τ t ) F - ϕ F u = 0 , K u ( ϕ u + 2 τ t - 4 ξ x ) + K u u ϕ = 0 , D u ( 2 ξ x - 2 τ t - ϕ u ) - D u u ϕ + K u ( 3 ϕ x x u - ξ x x x ) = 0 .
Investigating the compatibility of system Equation (8) we find that the last two equations of system Equation (8) are two identities (substituting the third, the fourth and the seventh equations of system Equation (8) into the last two one can yield this conclusion). Furthermore, with the aid of K ( u ) 0 , the second equation and the derivative of the third equation with respect to x imply that ξ x x = 2 ϕ x u . Thus, the determining Equation (8) is reduced to
τ x = τ u = ξ t = ξ u = ϕ u u = 0 , τ t t = 2 ϕ t u , ϕ x u K u = 0 , ξ x x = 2 ϕ x u , ϕ K u + 2 ( τ t - 2 ξ x ) K = 0 , ϕ x x x K u - 2 ϕ x D u + 2 ϕ x x x u K = 0 , ϕ D u + 2 ( τ t - ξ x ) D + 2 ϕ x x u K = 0 , ϕ t t + ϕ x x x x K - ϕ x x D + ( ϕ u - 2 τ t ) F - ϕ F u = 0 .
The first two equations of system Equation (9) do not contain arbitrary elements. Integration of them yields
τ = τ ( t ) , ξ = ξ ( x ) , ϕ = ϕ 1 ( t , x ) u + ϕ 0 ( t , x ) , ϕ 1 ( t , x ) = 1 2 τ t + α ( x ) ,
Thus, group classification of Equation (1) reduces to solving the rest equations of of system Equation (9).
Splitting the rest of the system Equation (9) with respect to the arbitrary elements and their non-vanishing derivatives gives the equations τ t = 0 , ξ x = 0 , ϕ = 0 for the coefficients of the operators from A ker of Equation (1). As a result, we obtain the following assertion.
Proposition 1.
The Lie algebra of the kernel of principal groups of Equation (1) is a two-dimensional algebra A ker = t , x .
In order to make the classification as simple as possible, we next look for equivalence transformations of class Equation (1), and then solve system Equation (9) under these transformations (Here we do not give detailed statements about equivalence transformation to avoid tediousness, one can see the book [21] or the paper [54] for details). An equivalence transformation is a nondegenerate change of the variables t, x and u taking any equation of the form Equation (1) into an equation of the same form, generally speaking, with different K ( u ) , D ( u ) and F ( u ) . The set of all equivalence transformations forms the equivalence group G . To find the connected component of the unity of G , we have to investigate Lie symmetries of the system that consists of Equation (1) and some additional conditions, i.e.,
u t t = - K u u x u x x x - K u x x x x + D u u x 2 + D u x x + F , K t = 0 , K x = 0 , D t = 0 , D x = 0 , F t = 0 , F x = 0 .
That is to say we must seek for an operator of the Lie algebra A of G in the form
X = τ ( t , x , u ) t + ξ ( t , x , u ) x + ϕ ( t , x , u ) u + π ( t , x , u , K , D , F ) K + ρ ( t , x , u , K , D , F ) D + θ ( t , x , u , K , D , F ) F .
Here u and K , D , F are considered as different variables: u is on the space ( t , x ) and K , D , F is on the extended space ( t , x , u ) . The coordinates τ , ξ , ϕ of the operator Equation (12) are sought as functions of t, x, u while the coordinates π , ρ , θ are sought as functions of t, x, u and K , D , F . Applying
pr ( 4 ) X = v + ϕ t t u t t + ϕ x u x + ϕ x x u x x + ϕ x x x u x x x + ϕ x x x x u x x x x + π t K t + π x K x + π u K u + ρ t D t + ρ x D x + ρ u D u + θ t F t + θ x F x
to Equation (11) we get the infinitesimal criterion
ϕ t t = - u x u x x x π u - K u ϕ x u x x x - K u u x ϕ x x x - π u x x x x - K ϕ x x x x - u x 2 ρ u + 2 D u u x ϕ x + ρ u x x + D ϕ x x + θ π x = 0 , π t = 0 , ρ x = 0 , ρ t = 0 , θ x = 0 , θ t = 0 ,
which must be satisfied whenever Equation (11) is satisfied. Substituting the formulae of ϕ t t , ϕ x , ϕ x x , ϕ x x x , ϕ x x x x , π x , π t , π u , ρ x , ρ t , ρ u , θ x , and θ t into Equation (13) we get equations of t, x, u, f, and the partial derivatives of ξ , τ , ϕ , π , ρ , θ . Replacing u t t , K x , K t , D x , D t , F x , and F t by the right hand side of Equation (11) whenever they occur, and equating the coefficients of various independent monomials to zero, we obtain
τ x = 0 , τ u = 0 , ξ t = 0 , ξ u = 0 , ϕ u u = 0 , τ t t - 2 ϕ t u = 0 , K u ( ϕ x u - ξ x x ) = 0 , K ( 4 ξ x - 2 τ t ) - π = 0 , K u ( 4 ξ x - 2 τ t ) - π u - K u π K - D u π D - F u π F = 0 , 2 D u ( ϕ u + τ t - ξ x ) - K u ( 3 ϕ x x u - ξ x x x ) - ρ u - K u ρ K - D u ρ D - F u ρ F = 0 , 2 D ( τ t - ξ x ) - K ( 6 ϕ x x u - 4 ξ x x x ) + ρ = 0 , D ( 2 ϕ x u - ξ x x ) - K u ϕ x x x - K ( 4 ϕ u x x x - ξ x x x x ) + 2 D u ϕ x = 0 , ϕ x x D - ϕ x x x x K + θ - ϕ t t - F ( ϕ u - 2 τ t ) = 0 , π t - K u ϕ t = 0 , π x - K u ϕ x = 0 , ρ t - D u ϕ t = 0 , ρ x - D u ϕ x = 0 , θ t - F u ϕ t = 0 , θ x - F u ϕ x = 0 ,
which can be reduced to
τ = c 1 t + c 4 , ξ = c 2 x + c 5 , ϕ = c 3 u + c 6 , π = 2 ( 2 c 2 - c 1 ) K ( u ) , ρ = 2 ( c 2 - c 1 ) D ( u ) , θ = ( c 3 - 2 c 1 ) F ( u ) ,
where c 1 , c 2 , c 3 , c 4 , c 5 and c 6 are arbitrary constants.
Thus the Lie algebra of G for class Equation (11) is
A = t , x , u , t t - 2 K K - 2 D D - 2 F F , x x + 4 K K + 2 D D , u u + F F .
Continuous equivalence transformations of class Equation (11) are generated by the operators from A . In fact, G contains the following continuous transformations:
t ˜ = t ε 1 + ε 4 , x ˜ = x ε 2 + ε 5 , u ˜ = u ε 3 + ε 6 , K ˜ = K ε 1 - 2 ε 2 4 , D ˜ = D ε 1 - 2 ε 2 2 , F ˜ = F ε 1 - 2 ε 3 ,
where ε 1 , ε 2 , , ε 6 are arbitrary constants.
Solve the system of determining Equation (9) under the above equivalence group G , we can obtain sixteen inequivalent equations of class Equation (1) with respect to the transformations from G .
Theorem 1.
A complete set of G -inequivalent extensions of A max A ker for Equation (1) is exhausted by ones given in Table 1.
Proof. 
To obtain the classification result we need to solve the system Equation (9) using the compatibility method [49,50]. The basic idea of this method is based on the fact that the substitution of the coefficients of any operator from A max \ A ker into the classifying equations results in nonidentity equations for arbitrary elements (see [49,50] for more details and exhaustive examples of applications). In our case the procedure of looking for the possible cases mostly depends on the fifth equation of system Equation (9). For any symmetry operator, equations ϕ K u + 2 ( τ t - 2 ξ x ) K = 0 give some equations on K of the general form
( a u + b ) K u + c K = 0 ,
where a, b and c are constants. For all operators from A max the number k of such independent equations is not greater than 2; otherwise they form an incompatible system on K. k is an invariant value for the transformations from G . Therefore, there exist three inequivalent cases for the value of k: (i) k = 0 : K ( u ) is arbitrary; (ii) k = 1 : K ( u ) = e μ u or K ( u ) = u μ ( μ 0 ) mod G , and (iii) k = 2 : K ( u ) = 1 mod G . Therefore, to complete the classification we have to consider all possible cases of the values of K ( u ) . We attempted to present our calculations in reasonable detail so that verification would be feasible.
Case 1: k = 0 . In this case, the fifth equation of system Equation (9) imply that ϕ = 0 and τ t - 2 ξ x = 0 . With this condition, the system Equation (9) can be reduced to
τ t t = 0 , ξ x x = 0 , 2 ( τ t - 2 ξ x ) = 0 , ( τ t - ξ x ) D = 0 , τ t F = 0 .
If D 0 or F 0 , which is corresponding to the Case 1. If D = F = 0 , we immediately arrive at Case 2 of Table 1 from system Equation (15).
Case 2: k = 1 . Here K { e μ u , u μ , μ 0 } mod G ^ and there exists v A max with ϕ 0 , otherwise there is no additional extension of the maximal Lie invariance algebra in comparison with the case k = 0 .
Case 2.1: Let us investigate the first possibility K = e μ u . In this case, the third equation and the fourth equation imply ϕ x u = 0 , ξ x x = 0 . Thus, from the fifth equation of system Equation (9), we can get ϕ x K u = 0 , which implies ϕ x = 0 , and thus leads to α ( x ) = c 3 , ϕ 0 ( t , x ) = ϕ 0 ( t ) , where c 3 is a constant. Furthermore, substituting K = e μ u into the fifth equation of system Equation (9) and using Equation (10), we can obtain ϕ 1 = 0 , μ ϕ 0 + 2 ( τ t - 2 ξ x ) = 0 , and thus leads to 1 2 τ t + c 3 = 0 , ϕ t 0 = 0 . Therefore, the system Equation (9) can be reduced to
τ t = - 2 c 3 , ξ x x = 0 , ϕ = c 4 μ ϕ + 2 ( τ t - 2 ξ x ) = 0 , ϕ D u + 2 ( τ t - ξ x ) D = 0 , - 2 τ t F - ϕ F u = 0 .
where c 4 is a constant. From the last two equations of system Equation (16), we can obtain D = d e ν u and F = f e γ u mod G , where d , f are two arbitrary constants. Substituting them into the system Equation (16) we can obtain
τ t = - 2 c 3 , ξ x x = 0 , ϕ = c 4 , μ ϕ + 2 ( τ t - 2 ξ x ) = 0 , d [ ν ϕ + 2 ( τ t - ξ x ) ] = 0 , f [ ϕ γ + 2 τ t ] = 0 .
The last two equations of the above system are two classifying conditions, which can be decomposed into four cases:
( i ) d = 0 , f = 0 ; ( i i ) d = 0 , f 0 ; ( i i i ) d 0 , f 0 ; ( i v ) d 0 , f = 0 .
(i) For this case the system Equation (17) can be reduced to
τ t = - 2 c 3 , ξ x x = 0 , ϕ = c 4 , μ ϕ + 2 ( τ t - 2 ξ x ) = 0 .
Solving this system, we can obtain τ = - 2 c 3 t + c 1 , ξ = ( μ 4 c 4 - c 3 ) x + c 2 , ϕ = c 4 , which is corresponding to Case 3.
(ii) For d = 0 , f 0 , the system Equation (17) can be reduced to
τ t = - 2 c 3 , ξ x x = 0 , ϕ = c 4 , μ ϕ + 2 ( τ t - 2 ξ x ) = 0 , γ ϕ + 2 τ t = 0 ,
from which we can obtain c 3 = γ 4 c 4 , τ = - γ 2 c 4 t + c 1 , ξ = μ - γ 4 c 4 x + c 2 , ϕ = c 4 , which is corresponding to Case 4.
(iii) When d 0 , f 0 , the system Equation (17) can be reduced to
τ t = - 2 c 3 , ξ x x = 0 , ϕ = c 4 , μ ϕ + 2 ( τ t - 2 ξ x ) = 0 , ν ϕ + 2 ( τ t - ξ x ) = 0 , ϕ γ + 2 τ t = 0 .
Solving the above system, we can obtain γ = 2 ν - μ , c 3 = γ 4 c 4 , τ = - γ 2 c 4 t + c 1 , ξ = μ - γ 4 c 4 x + c 2 , ϕ = c 4 , which is corresponding to Case 5.
(iv) For d 0 , f = 0 , the system Equation (17) can be reduced to
τ t = - 2 c 3 , ξ x x = 0 , ϕ = c 4 , μ ϕ + 2 ( τ t - 2 ξ x ) = 0 , ν ϕ + 2 ( τ t - ξ x ) = 0 ,
from which we can obtian c 3 = 2 ν - μ 4 c 4 , τ = μ - 2 ν 2 c 4 t + c 1 , ξ = μ - ν 2 c 4 x + c 2 , ϕ = c 4 , which is corresponding to Case 5.
Case 2.2: Consider the case K = u μ ( μ 0 ) . The third equation of system Equation (9) implies μ ϕ x u u μ - 1 = 0 , thus we have ϕ x u = 0 . Then, from the fourth equation and the fifth equation we can obtain ξ x x = 0 , ϕ x = 0 , which implies α ( x ) = c 3 , ϕ 0 ( t , x ) = ϕ 0 ( t ) , ξ ( x ) = c 1 x + c 2 , and thus the system Equation (9) can be reduced to
ξ ( x ) = c 1 x + c 2 , ϕ = ( 1 2 τ t + c 3 ) u + ϕ 0 ( t ) , μ ϕ + 2 ( τ t - 2 ξ x ) u = 0 , ϕ D u + 2 ( τ t - ξ x ) D = 0 , 1 2 τ t t t u + ϕ t t 0 + ( ϕ u - 2 τ t ) F - ϕ F u = 0 .
The second and the third equations of system Equation (18) implies
[ ( μ 2 + 2 ) τ t + μ c 3 - 4 ξ x ] u + μ ϕ 0 ( t ) = 0 , ( μ 2 + 2 ) τ t t = 0 .
Thus we have ϕ 0 ( t ) = 0 , and there exist two cases: (I) τ t t 0 or (II) τ t t = 0 .
(I) For τ t t 0 , from system Equation (19) we have μ = - 4 , c 1 = - c 3 . Thus the system Equation (18) can be reduced to
ϕ = ( 1 2 τ t + c 3 ) u , ϕ D u + 2 ( τ t + c 3 ) D = 0 , 1 2 τ t t t u + ( ϕ u - 2 τ t ) F - ϕ F u = 0 .
Substituting the first equation of system Equation (20) into the second equation of system Equation (20) and then take the derivative for it with respect to the variable t, we have
( u D u + 4 D ) τ t t = 0 , c 3 ( u D u + 2 D ) = 0 .
Solving the first equation of the above system we can obtain D ( u ) = d u - 4 , where d is an arbitrary constant. Substituting it into the second equation of the system Equation (21), we can obtain a classifying condition d c 3 = 0 , which imply there exist two cases: (i) c 3 = 0 or (ii) c 3 0 .
(i) For c 3 = 0 , we have c 1 = 0 , and thus the last equations of system Equation (20) can be reduced to
u F u + 3 F - τ t t t τ t u = 0 ,
which implies τ t t t τ t = ω , where ω is a constant. When ω = 0 , from Equation (22) we have F ( u ) = f u - 3 , where f is a constant. Furthermore, τ t t t = 0 and c 1 = 0 imply τ = c 6 t 2 + c 5 t + c 4 , ξ = c 2 , and thus ϕ = ( c 6 t + c 5 2 ) u , which is corresponding to Case 6. When ω > 0 , we can get τ = c 4 + c 5 e ω t + c 6 e - ω t , ξ = c 2 , ϕ = 1 2 ( c 5 ω e ω t - c 6 ω e - ω t ) u , and F ( u ) = 1 4 ω u + f u - 3 . which is corresponding to Case 7. When ω < 0 , we can get τ = c 4 + c 5 sin ( - ω t ) + c 6 cos ( - ω t ) , ξ = c 2 , ϕ = 1 2 ( c 5 - ω cos ( - ω t ) - c 6 - ω sin ( - ω t ) ) u , and F ( u ) = 1 4 ω u + f u - 3 . Which is corresponding to Case 8.
(ii) For c 3 0 , we have d = 0 and D ( u ) = 0 , and thus the last equations of system Equation (20) can be reduced to
1 2 τ t ( u F u + 3 F ) - 1 2 τ t t t u + c 3 ( u F u - F ) = 0 ,
Taking the derivative of Equation (23) with respect to the variable t, we have
u F u + 3 F - τ t t t t τ t t u = 0 .
which implies τ t t t t τ t t = λ .
When λ = 0 , we have τ t t t t = 0 and F ( u ) = f u - 3 , where f is an arbitrary constant. Substituting F ( u ) = f u - 3 into Equation (23), it can be reduced to
4 c 3 f u - 3 + 1 2 τ t t t u = 0 .
which implies τ t t t = 0 and f = 0 . Then we can get τ = c 6 t 2 + c 5 t + c 4 , ξ = - c 3 x + c 2 , and thus ϕ = ( c 6 t + c 5 2 + c 3 ) u , which is corresponding to Case 9.
When λ 0 , we can get τ t t t t 0 , and the Equation (23) can be separated into
1 2 τ t ( u F u + 3 F ) - 1 2 τ t t t u = 0 , c 3 ( u F u - F ) = 0 ,
From the second equation of Equation (24), we can get F ( u ) = f u , where f is an arbitrary constant. Substituting F ( u ) = f u into the first equation of Equation (24), we can obtain
4 f u - τ t t t τ t u = 0 .
which implies τ t t t τ t = 4 f and f 0 . When f > 0 , we can get τ = c 4 + c 5 e 2 f t + c 6 e - 2 f t , ξ = - c 3 x + c 2 , ϕ = ( c 5 f e 2 f t - c 6 f e - 2 f t + c 3 ) u , and F ( u ) = f u ( f > 0 ) , which is corresponding to Case 10. When f < 0 , we can get τ = c 4 + c 5 sin ( 2 - f t ) + c 6 cos ( 2 - f t ) , ξ = - c 3 x + c 2 , ϕ = ( c 5 - f cos ( 2 - f t ) - c 6 - f sin ( 2 - f t ) + c 3 ) u , and F ( u ) = f u ( f < 0 ) , which is corresponding to Case 11.
(II) τ t t = 0 implies τ t is a constant. Let 1 2 τ t + c 3 c , the system Equation (18) can be reduced to
τ t t = 0 , ξ ( x ) = c 1 x + c 2 , ϕ = c u , μ ϕ + 2 ( τ t - 2 ξ x ) u = 0 , ϕ D u + 2 ( τ t - ξ x ) D = 0 , ϕ F u - ( ϕ u - 2 τ t ) F = 0 .
As one can see, for any symmetry operator, equations ϕ D u + 2 ( τ t - ξ x ) D = 0 and ϕ F u + ( ϕ u - 2 τ t ) F = 0 give some equations (not greater than 2 each) on D and F of the general form
c u D u + r D = 0 , c u F u + s F = 0 ,
where c, r and s are constants. Solving this system, we can obtain D = d u ν and F = f u γ mod G , where d , f are two arbitrary constants. Substituting them into the system Equation (25) we can get
τ t t = 0 , ξ ( x ) = c 1 x + c 2 , ϕ = c u , c μ + 2 ( τ t - 2 ξ x ) = 0 , d [ c ν + 2 ( τ t - ξ x ) ] = 0 , f [ c γ - ( c - 2 τ t ) ] = 0 .
The last two equations of the above system are two classifying conditions, which can be decomposed into four cases:
( i ) d = 0 , f = 0 ; ( i i ) d = 0 , f 0 ; ( i i i ) d 0 , f 0 ; ( i v ) d 0 , f = 0 .
(i) For this case, the system Equation (26) can be reduced to
τ t t = 0 , ξ ( x ) = c 1 x + c 2 , ϕ = c u , c μ + 2 ( τ t - 2 ξ x ) = 0 .
Solving the above system, we can obtain τ = ( 2 c 1 - μ 2 c ) t + c 0 , ξ = c 1 x + c 2 , ϕ = c u , which is corresponding to Case 12.
(ii) When d = 0 , f 0 , the system Equation (26) can be reduced to
τ t t = 0 , ξ ( x ) = c 1 x + c 2 , ϕ = c u , c μ + 2 ( τ t - 2 ξ x ) = 0 , c γ - ( c - 2 τ t ) = 0 ,
from which we can obtain τ = 1 - γ 2 c t + c 0 , ξ = 1 - γ + μ 4 c x + c 2 , ϕ = c u , which is corresponding to Case 13.
(iii) For d 0 , f 0 , the system Equation (26) can be reduced to
τ t t = 0 , ξ ( x ) = c 1 x + c 2 , ϕ = c u , c μ + 2 ( τ t - 2 ξ x ) = 0 , c ν + 2 ( τ t - ξ x ) = 0 , c γ - ( c - 2 τ t ) = 0 .
Solving the above system, we can obtain γ = 1 - μ + 2 ν , τ = μ - 2 ν 2 c t + c 0 , ξ = μ - ν 2 c x + c 2 , ϕ = c u , which is corresponding to Case 14.
(iv) When d 0 , f = 0 , the system Equation (26) can be reduced to
τ t t = 0 , ξ ( x ) = c 1 x + c 2 , ϕ = c u , c μ + 2 ( τ t - 2 ξ x ) = 0 , c ν + 2 ( τ t - ξ x ) = 0 ,
from which we can obtain τ = μ - 2 ν 2 c t + c 0 , ξ = μ - ν 2 c x + c 2 , ϕ = c u , which is corresponding to Case 14.
Case 3: k = 2 . The assumption of two independent equations of form of the fifth equation of system Equation (9) for K yields K = const, i.e. K = 1 mod G . From the fifth equation of system Equation (9) we have τ t = 2 ξ x . Thus, we can get τ t t = 0 , ξ x x = 0 , ϕ t u = ϕ x u = 0 , which implies system Equation (9) can be reduced to
τ t t = ξ x x = ϕ t u = ϕ x u = 0 , τ t - 2 ξ x = 0 , 2 ϕ x D u = 0 , ϕ D u + 2 ( τ t - ξ x ) D = 0 , ϕ F u - ( ϕ u - 2 τ t ) F + ϕ x x D - ϕ x x x x - ϕ t t = 0 .
In a way similar to the proof in case 2.2, from the fourth equation of system Equation (27) we obtain the following different values of D ( u ) : D ( u ) = , or d e ν u , or d u ν ( d 0 , ν 0 ) , or d. Furthermore, taking derivative of the last equation of system Equation (27) with respect to u, we can get ϕ F u u + 2 τ t F u = 0 , which give some equations (not greater than 3) on F of the general form
( m u + n ) F u u + w F u = 0 ,
for any symmetry operator, where m, n and w are constants. Solving this equation up to G , we obtain the following different values of F ( u ) : F ( u ) = f e γ u , or f u γ ( f 0 , γ 0 ) , or f ln u mod G . Therefore, to complete the classification we have to consider all possible combinations of the values of D ( u ) and F ( u ) . Substituting all the different combinations of D ( u ) and K ( u ) to the system Equation (27), we can find three nontrivial inequivalent combinations: (I) D ( u ) = d e ν u ( d 0 , ν 0 ) , F ( u ) = f e γ u ( f 0 , γ 0 ) ; (II) D = d u ν ( d 0 , ν 0 ) , F = f u γ ( f 0 , γ 0 ) and (III) D = d , F = f u .
(I) For the first combination, substituting D ( u ) = d e ν u ( d 0 , ν 0 ) into the fourth equation of the system Equation (27), we can obtain
ϕ x D u = 0 , ϕ t D u = 0 .
Thus we have ϕ = c u + c 3 , and the system Equation (27) can be reduced to
τ t t = 0 , ξ x x = 0 , ϕ = c u + c 3 , τ t - 2 ξ x = 0 , ( c u + c 3 ) ν + 2 ( τ t - ξ x ) = 0 , ( c u + c 3 ) γ - ( c - 2 τ t ) = 0 .
The above system implies τ = - c 3 ν t + c 1 , ξ = - 1 2 c 3 ν x + c 2 , ϕ = c 3 and γ = 2 ν , where c i ( i = 1 , , 3 ) are arbitrary constants, which is corresponding to Case 5 with μ = 0 .
(II) For the second combination, substituting D = d u ν ( d 0 , ν 0 ) into the fourth equation of the system Equation (27), we can obtain
ϕ x D u = 0 , ϕ t D u = 0 .
Thus we have ϕ = c u + c 3 , and the system Equation (27) can be reduced to
τ t t = 0 , ξ x x = 0 , ϕ x = 0 , ϕ = c u + c 3 , τ t - 2 ξ x = 0 , c 3 ν + [ c ν + 2 ( τ t - ξ x ) ] u = 0 , ( c u + c 3 ) γ - ( c - 2 τ t ) u = 0 .
The above system implies τ = - c ν t + c 1 , ξ = - 1 2 c ν x + c 2 , ϕ = c u and γ = 1 + 2 ν , where c i ( i = 1 , , 3 ) are arbitrary constants, which is corresponding to Case 14 with μ = 0 .
(III) For the last combination, substituting them into the original system Equation (27), we can obtain
τ t t = 0 , ξ x x = 0 , ϕ = c u + b ( t , x ) , τ t - 2 ξ x = 0 , d ( ξ x - τ t ) = 0 , b t t + b x x x x - d b x x + f [ ( c - 2 τ t ) u - ( c u + b ( t , x ) ) ] = 0 .
The fifth equation of the above system is a classifying condition, which leads to two cases: (i) d = 0 ; (ii) d 0 .
(i) For d = 0 , the last equation of the system Equation (28) can be reduced to
b t t + b x x x x - f b ( t , x ) - 2 τ t f u = 0 ,
from which for f = 0 we can obtian τ = 2 c 3 t + c 1 , ξ = c 3 x + c 2 , ϕ = c u + b ( t , x ) , and b ( t , x ) is satisfied with the equation:
b t t + b x x x x = 0
which is corresponding to Case 15. For f 0 , we can obtain τ = c 1 , ξ = c 2 , ϕ = c u + b ( t , x ) , and b ( t , x ) is satisfied with the equation:
b t t + b x x x x - f b ( t , x ) = 0
which is corresponding to Case 16 with d = 0 .
(ii) For d 0 , the system Equation (28) implies τ = c 1 , ξ = c 2 , ϕ = c u + b ( t , x ) , where c i ( i = 1 , , 3 ) are arbitrary constants, and b ( t , x ) is satisfied with the equation:
b t t - d b x x + b x x x x - f b ( t , x ) = 0 .
which is corresponding to Case 16.   ☐

3. Symmetry Reduction and Exact Solutions

The Lie symmetry operators in the Table 1 found as a result of the group classification problem can be applied to construction exact solutions of the corresponding equations by means of symmetry reduction. The method of symmetry reduction with respect to subalgebras of Lie invariance algebras is well-known and quite algorithmic to use in most cases, we refer to the standard textbooks on the subject [21,59]. Below, we perform symmetry reduction for one classification model.
We consider Equation (1) with the fixed values of the parameter-functions K = u μ , D = u ν and F = u 1 - μ + 2 ν , i.e.,
u t t = - [ u μ u x x x ] x + [ u ν u x ] x + u 1 - μ + 2 ν ,
where μ , ν are two arbitrary constants. As shown in case 14 with d = f = 1 of the Table 1, the Equation (29) admits the three-dimensional Lie invariance algebra g generated by the operators
Q 1 = t , Q 2 = x , Q 3 = 1 2 ( μ - 2 ν ) t t + 1 2 ( μ - ν ) x x + u u ,
which is equivalent to
Q 1 = t , Q 2 = x , Q 3 = t t + μ - ν μ - 2 ν x x + 2 μ - 2 ν u u .
Because these operators satisfy the commutation relations
[ Q 1 , Q 3 ] = Q 1 , [ Q 2 , Q 3 ] = μ - ν μ - 2 ν Q 2 ,
and thus the corresponding symmetry algebra g is a realization of the algebra A 3 , 5 α [60], where 0 < | μ - ν μ - 2 ν | < 1 . A complete list of inequivalent non-zero one-dimensional subalgebras of g is exhausted by the algebras [60]
Q 1 , Q 2 , Q 3 , Q 1 + α Q 2 ,
where α = ± 1 . Lie reduction of Equation (29) to ordinary differential equations (ODEs) can be made with the one-dimensional subalgebra Q 1 , Q 2 , Q 3 and Q 1 + α Q 2 . The associated ansatz and the reduced ODEs are listed in Table 2.
Solving the reduced ODEs in Table 2, we can obtain some exact solutions of the corresponding Equation (29) by using ansatz. For example, solving the ODE corresponding to the case 4 of Table 2 with μ = - 1 , ν = - 1 by using the extended tanh function method [61,62], we can obtain the following exact travelling wave solutions for the nonlinear beam-like Equation (29) with μ = - 1 , ν = - 1 :
Case 1. Solitary wave solutions:
u = 12 R tanh 2 [ - R ( x - α t ) ] - 8 R + 1 , u = 12 R coth 2 [ - R ( x - α t ) ] - 8 R + 1 ;
where R < 0 .
Case 2. Rational solutions:
u = - 12 ( x - α t ) - 2 + 1 .
Case 3. Traingular periodic wave solutions:
u = - 12 R tan 2 [ R ( x - α t ) ] - 8 R + 1 , u = - 12 R cot 2 [ R ( x - α t ) ] - 8 R + 1 ,
where R > 0 .

4. Conclusions and Discussion

In this paper we present a complete group classification of the generalized nonlinear beam Equation (1) by using the classical Lie-Ovsiannikov method based on equivalence transformations. The main results on classification are collected in Table 1 where we list sixteen inequivalent cases of extensions with the corresponding Lie invariance algebras. It should be noted that the right-hand side of Equation (1) is exactly the same as the generalized thin film Equation (1.1) in [58]. Compare the classification results of Equation (1.1) in [58], we know that the generalized nonlinear beam Equation (1) have more inequivalent classification cases than the generalized thin film equation (eight cases). Furthermore, the maximal extension of Lie invariance algebra for the generalized nonlinear beam Equation (1) is five dimension, while for the generalized thin film equation it is four dimension. Therefore, the classification results of the two equations are different.
For one classification models from the list we perform corresponding Lie symmetry reduction which are presented in Table 2. This enabled to obtain some exact solutions for the equation under consideration, including solitary wave solutions, triangular periodic wave solutions and rational solutions. These results may lead to further applications in physics and engineering such as tests in numerical solutions of Equation (1) and as trial functions for application of variational approach in the analysis of different perturbed versions of Equation (1). Other topics including nonclassical symmetry, non-Lie exact solutions and physical applications of class Equation (1) will be studied in subsequent publication.

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China (Grant No. 11501204), the Natural Science Foundation of Shanghai (Grant No. 15ZR1408300), and Shanghai Key Laboratory of Intelligent Information Processing (Grant No. IIPL-2014-001).

Author Contributions

All the authors D.J.H., X.X.L. and S.C.Y.(check all the computations) make contributions to Section 2 and Section 3. D.J.H make contribution to Section 1 and Section 4.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ames, W.F. Nonlinear Partial Differential Equations in Engineering; Academic: New York, NY, USA, 1972; Volume II, pp. 50–52. [Google Scholar]
  2. Ames, W.F.; Adams, E.; Lohner, R.G. Group properties of utt = [f(u)ux]x. Int. J. Non-Linear Mech. 1981, 16, 439–447. [Google Scholar] [CrossRef]
  3. Galaktionov, V.A. The formation of shocks and fundamental solution of a fourth-order quasilinear Boussinesq-type equation. Nonlinearity 2009, 22, 239–257. [Google Scholar] [CrossRef]
  4. Favini, A.; Goldstein, G.R.; Goldstein, J.A.; Romanelli, S. Classification of general Wentzell boundary conditions for fourth order operators in one space dimension. J. Math. Anal. Appl. 2007, 335, 219–235. [Google Scholar] [CrossRef]
  5. McKenna, P.J.; Walter, W. Travelling waves in a suspension bridges. SIAM J. Appl. Math. 1990, 50, 703–715. [Google Scholar] [CrossRef]
  6. Ammann, O.H.; Karman, T.V.; Woodruff, G.B. The Failure of the Tacoma Narrow Bridge; Federal Works Agency: Washington, DC, USA, 1941. [Google Scholar]
  7. Champreys, A.R.; McKenna, P.J.; Zegeling, P.A. Solitary waves in nonlinear beam equations: Stability, fission and fusion. Nonlinear Dynam. 2000, 21, 31–53. [Google Scholar] [CrossRef]
  8. McKenna, P.J.; Walter, W. Nonlinear oscillations in a suspension bridge. Arch. Ration. Mech. Anal. 1987, 98, 167–177. [Google Scholar] [CrossRef]
  9. Chen, Y.; McKenna, P.J. Traveling waves in a nonlinearly suspended beam: Theoretical results and numerical observations. J. Differ. Equ. 1997, 136, 325–355. [Google Scholar] [CrossRef]
  10. Chen, Y.; McKenna, P.J. Traveling waves in a nonlinearly suspended beam: Some computational results and four open questions. Phil. Trans. R. Soc. Lond. A 1997, 355, 2175–2184. [Google Scholar] [CrossRef]
  11. Choy, Y.S.; Jen, K.S.; McKenna, P.J. The structure of the solution set for periodic oscillations in a suspension bridge model. IMA J. Appl. Math. 1991, 47, 283–306. [Google Scholar] [CrossRef]
  12. Humphreys, L.D. Numerical mountain pass solutions of a suspension bridge equation. Nonlinear Anal. TMA 1997, 35, 1811–1826. [Google Scholar] [CrossRef]
  13. Lazer, A.C.; McKenna, P.J. Large Amplitude periodic oscillation in suspension bridges: Some new connections with nonlinear analysis. SIAM Rev. 1990, 32, 537–578. [Google Scholar] [CrossRef]
  14. Humphreys, L.D.; McKenna, P.J. Multiple periodic solutions for a nonlinear suspension bridge equation. IMA J. Appl. Math. 1999, 63, 37–49. [Google Scholar] [CrossRef]
  15. Doole, S.H.; Hogan, S.J. The nonlinear dynamics of suspension bridges under harmonic forcing. Appl. Nonlinear Math. Rep. 1996, 76, 127–128. [Google Scholar]
  16. Doole, S.H.; Hogan, S.J. A piecewise linear suspension bridge model nonlinear dynamics and orbit continuation. Dynam. Stabil. Syst. 1996, 11, 19–47. [Google Scholar] [CrossRef]
  17. Bruzón, M.S.; Camacho, J.C.; Gandarias, M.L. Similarity reductions of a nonlinear model for vibrations of beams. PAMM Proc. Appl. Math. Mech. 2007, 7, 2040063–2040064. [Google Scholar] [CrossRef]
  18. Camacho, J.C.; Bruzón, M.S.; Ramírez, J.; Gandarias, M.L. Exact travelling wave solutions of a beam equation. J. Nonlinear Math. Phys. 2011, 18, 33–49. [Google Scholar] [CrossRef]
  19. Gao, Y.X. Quasi-periodic Solutions of the General Nonlinear Beam Equations. Commun. Math. Res. 2012, 28, 51–64. [Google Scholar]
  20. Ovsiannikov, L.V. Group properties of the nonlinear heat-conduction equation. Dokl. Akad. Nauk SSSR 1959, V.125, 492–495. (In Russian) [Google Scholar] [PubMed]
  21. Ovsiannikov, L.V. Group Analysis of Differential Equations; Academic Press: New York, NY, USA, 1982. [Google Scholar]
  22. Bluman, G.; Cheviakov, A.; Anco, S. Applications of Symmetry Methods to Partial Differential Equations; Springer: Berlin, Germany, 2010. [Google Scholar]
  23. Barone, A.; Esposito, F.; Magee, C.G.; Scott, A.C. Theory and applications of the sine-Gordon equation. Riv. Nuovo Cimento 1971, 1, 227–267. [Google Scholar] [CrossRef]
  24. Arrigo, D.J. Group properties of uxx u y m uyy = f(u). Int. J. Non-Linear Mech. 1991, 26, 619–629. [Google Scholar] [CrossRef]
  25. Pucci, E.; Salvatori, M.C. Group properties of a class of semilinear hyperbolic equations. Int. J. Non-Linear Mech. 1986, 21, 147–155. [Google Scholar] [CrossRef]
  26. Torrisi, M.; Valenti, A. Group properties and invariant solutions for infinitesimal transformations of a nonlinear wave equation. Int. J. Non-Linear Mech. 1985, 20, 135–144. [Google Scholar] [CrossRef]
  27. Donato, A. Similarity analysis and nonlinear wave propagation. Int. J. Non-Linear Mech. 1987, 22, 307–314. [Google Scholar] [CrossRef]
  28. Ibragimov, N.H.; Torrisi, M.; Valenti, A. Preliminary group classification of equations vtt = f(x, vx)vxx + g(x, vx). J. Math. Phys. 1991, 32, 2988–2995. [Google Scholar] [CrossRef]
  29. Ibragimov, N.H. (Ed.) Lie Group Analysis of Differential Equations—Symmetries, Exact Solutions and Conservation Laws, V.1.; CRC Press: Boca Raton, FL, USA, 1994. [Google Scholar]
  30. Oron, A.; Rosenau, P. Some symmetries of the nonlinear heat and wave equations. Phys. Lett. A 1986, 118, 172–176. [Google Scholar] [CrossRef]
  31. Chikwendu, S.C. Non-linear wave propagation solutions by Fourier transform perturbation. Int. J. Non-Linear Mech. 1981, 16, 117–128. [Google Scholar] [CrossRef]
  32. Gandarias, M.L.; Torrisi, M.; Valenti, A. Symmetry classification and optimal systems of a non-linear wave equation. Int. J. Non-Linear Mech. 2004, 39, 389–398. [Google Scholar] [CrossRef]
  33. Pucci, E. Group analysis of the equation utt + λuxx = g(u, ux). Riv. Mat. Univ. Parma 1987, 12, 71–87. [Google Scholar]
  34. Bluman, G.W.; Cheviakov, A.F. Nonlocally related systems, linearization and nonlocal symmetries for the nonlinear wave equation. J. Math. Anal. Appl. 2007, 333, 93–111. [Google Scholar] [CrossRef]
  35. Bluman, G. W.; Kumei, S. Symmetries and Differential Equations; Springer: Berlin, Germany, 1989. [Google Scholar]
  36. Bluman, G.W.; Temuerchaolu; Sahadevan, R. Local and nonlocal symmetries for nonlinear telegraph equation. J. Math. Phys. 2005, 46, 023505. [Google Scholar] [CrossRef]
  37. Huang, D.J.; Ivanova, N.M. Group analysis and exact solutions of a class of variable coefficient nonlinear telegraph equations. J. Math. Phys. 2007, 48, 073507. [Google Scholar] [CrossRef]
  38. Huang, D.J.; Zhou, S.G. Group properties of generalized quasi-linear wave equations. J. Math. Anal. Appl. 2010, 366, 460–472. [Google Scholar] [CrossRef]
  39. Huang, D.J.; Zhou, S.G. Group-theoretical analysis of variable coefficient nonlinear telegraph equations. Acta Appl. Math. 2012, 117, 135–183. [Google Scholar] [CrossRef]
  40. Lahno, V.; Zhdanov, R.; Magda, O. Group classification and exact solutions of nonlinear wave equations. Acta Appl. Math. 2006, 91, 253–313. [Google Scholar] [CrossRef]
  41. Sophocleous, C.; Kingston, J.G. Cyclic symmetries of one-dimensional non-linear wave equations. Int. J. Non-Linear Mech. 1999, 34, 531–543. [Google Scholar] [CrossRef]
  42. Suhubi, E.S.; Bakkaloglu, A. Group properties and similarity solutions for a quasi-linear wave equation in the plane. Int. J. Non-Linear Mech. 1991, 26, 567–584. [Google Scholar] [CrossRef]
  43. Vasilenko, O.F.; Yehorchenko, I.A. Group classification of multidimensional nonlinear wave equations. Proc. Inst. Math. NAS Ukr. 2001, 36, 63–66. [Google Scholar]
  44. Cherniha, R.; Serov, M.; Rassokha, I. Lie symmetries and form-preserving transformations of reaction-diffusion-convection equations. J. Math. Anal. Appl. 2008, 342, 1363–1379. [Google Scholar] [CrossRef]
  45. Basarab-Horwath, P.; Lahno, V.I.; Zhdanov, R.Z. The structure of Lie algebras and the classification problem for partial differential equations. Acta Appl. Math. 2001, 69, 43–94. [Google Scholar] [CrossRef]
  46. Zhdanov, R.Z.; Lahno, V.I. Group classification of heat conductivity equations with a nonlinear source. J. Phys. A 1999, 32, 7405–7418. [Google Scholar] [CrossRef]
  47. Gazeau, J.P.; Winternitza, P. Symmetries of variable coefficient Korteweg-de Vries equations. J. Math. Phys. 1992, 33, 4087–4102. [Google Scholar] [CrossRef]
  48. Popovych, R.O.; Kunzinger, M.; Eshraghi, H. Admissible point transformations and normalized classes of nonlinear Schrödinger equations. Acta Appl. Math. 2010, 109, 315–359. [Google Scholar] [CrossRef]
  49. Nikitin, A.G.; Popovych, R.O. Group classification of nonlinear Schrödinger equations. Ukr. Math. J. 2001, 53, 1053–1060. [Google Scholar] [CrossRef]
  50. Popovych, R.O.; Ivanova, N.M. New results on group classification of nonlinear diffusion-convection equations. J. Phys. A 2004, 37, 7547–7565. [Google Scholar] [CrossRef]
  51. Ivanova, N.M.; Popovych, R.O.; Sophocleous, C. Group analysis of variable coefficient diffusion-convection equations. I. Enhanced group classification. Lobachevskii J. Math. 2010, 31, 100–122. [Google Scholar] [CrossRef]
  52. Huang, D.J.; Yang, Q.M.; Zhou, S.G. Lie symmetry classification and equivalence transformation of variable coefficient nonlinear wave equations with power nonlinearities. Chin. J. Contemp. Math. 2012, 33, 205–214. [Google Scholar]
  53. Huang, D.J.; Yang, Q.M.; Zhou, S.G. Conservation law classification of variable coefficient nonlinear wave equation with power Nonlinearity. Chin. Phys. B 2011, 20, 070202. [Google Scholar] [CrossRef]
  54. Huang, D.J.; Ivanova, N.M. Algorithmic framework for group analysis of differential equations and its application to generalized Zakharov-Kuznetsov equations. J. Differ. Equ. 2016, 260, 2354–2382. [Google Scholar] [CrossRef]
  55. Cherniha, R.; King, J.R. Lie symmetries of nonlinear multidimensional reaction-diffusion systems: I. J. Phys. A 2000, 33, 267–282. [Google Scholar] [CrossRef]
  56. Cherniha, R.; King, J.R. Lie symmetries of nonlinear multidimensional reaction-diffusion systems: II. J. Phys. A 2003, 36, 405–425. [Google Scholar] [CrossRef]
  57. Cherniha, R.; King, J.R. Lie symmetries and conservation laws of nonlinear multidimensional reaction-diffusion systems with variable diffusivities. IMA J. Appl. Math. 2006, 71, 391–408. [Google Scholar] [CrossRef]
  58. Cherniha, R.; Myroniuk, L. Lie Symmetries and Exact Solutions of the Generalized Thin Film Equation. J. Phys. Math. 2010, 2, P100508:1–P100508:19. [Google Scholar] [CrossRef]
  59. Olver, P.J. Application of Lie Groups to Differential Equations; Springer: New York, NY, USA, 1986. [Google Scholar]
  60. Patera, J.; Winternitz, P. Subalgebras of real three- and four-dimensional Lie algebras. J. Math. Phys. 1977, 18, 1449–1455. [Google Scholar] [CrossRef]
  61. Malfliet, W.; Hereman, W. The tanh method. I. Exact solutions of nonlinear evolution and wave equations. Phys. Scr. 1996, 54, 563–568. [Google Scholar] [CrossRef]
  62. Fan, E.G. Extended tanh-function method and its applications to nonlinear equations. Phys. Lett. A 2000, 277, 212–218. [Google Scholar] [CrossRef]
Table 1. Group classification of class Equation (1).
Table 1. Group classification of class Equation (1).
N K ( u ) D ( u ) F ( u ) Basis of A max
1 t , x
200 t , x , x x + 2 t t
3 e μ u 00 t , x , - 2 t t - x x , μ 4 x x + u
4 e μ u 0 f e γ u ( f 0 ) t , x , - γ 2 t t + μ - γ 4 x x + u
5 e μ u d e ν u ( d 0 ) f e ( 2 ν - μ ) u t , x , μ - 2 ν 2 t t + μ - ν 2 x x + u
6 u - 4 d u - 4 f u - 3 t , x , t t + 1 2 u u , t 2 t + t u u
7 u - 4 d u - 4 1 4 ω u + f u - 3 ( ω > 0 ) t , x , e ω t t + ω 2 e ω t u u , e - ω t t - ω 2 e - ω t u u
8 u - 4 d u - 4 1 4 ω u + f u - 3 ( ω < 0 ) t , x , sin ( - ω t ) t + - ω 2 cos ( - ω t ) u u ,
cos ( - ω t ) t - - ω 2 sin ( - ω t ) u u
9 u - 4 00 t , x , t t + 1 2 u u , t 2 t + t u u , - x x + u u
10 u - 4 0 f u ( f > 0 ) t , x , e 2 f t t + f e 2 f t u u , e - 2 f t t - f e - 2 f t u u ,
- x x + u u
11 u - 4 0 f u ( f < 0 ) t , x , sin ( 2 - f t ) t + - f cos ( 2 - f t ) u u ,
cos ( 2 - f t ) t - - f sin ( 2 - f t ) u u , - x x + u u
12 u μ 00 t , x , 2 t t + x x , - μ 2 t t + u u
13 u μ 0 f u γ ( f 0 ) t , x , 1 - γ 2 t t + 1 - γ + μ 4 x x + u u
14 u μ d u ν ( d 0 ) f u 1 - μ + 2 ν t , x , 1 2 ( μ - 2 ν ) t t + 1 2 ( μ - ν ) x x + u u
15100 t , x , 2 t t + x x , u u , v = b ( x , t ) u
where b t t + b x x x x = 0
161d f u t , x , u u , v = b ( x , t ) u
where b t t - d b x x + b x x x x - f b ( t , x ) = 0
Here, for Case 16, d, f can not be zero simultaneously.
Table 2. Reduced ODEs for Equation (29) ( α = ± 1 ) .
Table 2. Reduced ODEs for Equation (29) ( α = ± 1 ) .
N Snbalgebra Ansatz for u y Reduced ODE
1 < Q 1 > h ( y ) x - μ h μ - 1 h h - h μ h + ν h ν - 1 h 2 + h ν h + h 1 - μ + 2 ν = 0
2 < Q 2 > h ( y ) t h = h 1 - μ + 2 ν
3 < Q 3 > h ( y ) t 2 μ - 2 ν x t - μ - ν μ - 2 ν ( μ - ν ) 2 y 2 h + ( μ - ν ) ( 2 μ - 3 ν - 4 ) y h + 2 ( 2 - μ + 2 ν ) h =
( μ - 2 ν ) 2 ( - μ h μ - 1 h h - h μ h + ν h ν - 1 h 2 + h ν h + h 1 - μ + 2 ν )
4 < Q 1 + α Q 2 > h ( y ) x - α t h = - μ h μ - 1 h h - h μ h + ν h ν - 1 h 2 + h ν h + h 1 - μ + 2 ν

Share and Cite

MDPI and ACS Style

Huang, D.; Li, X.; Yu, S. Lie Symmetry Classification of the Generalized Nonlinear Beam Equation. Symmetry 2017, 9, 115. https://doi.org/10.3390/sym9070115

AMA Style

Huang D, Li X, Yu S. Lie Symmetry Classification of the Generalized Nonlinear Beam Equation. Symmetry. 2017; 9(7):115. https://doi.org/10.3390/sym9070115

Chicago/Turabian Style

Huang, Dingjiang, Xiangxiang Li, and Shunchang Yu. 2017. "Lie Symmetry Classification of the Generalized Nonlinear Beam Equation" Symmetry 9, no. 7: 115. https://doi.org/10.3390/sym9070115

APA Style

Huang, D., Li, X., & Yu, S. (2017). Lie Symmetry Classification of the Generalized Nonlinear Beam Equation. Symmetry, 9(7), 115. https://doi.org/10.3390/sym9070115

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop