Lie Symmetries , Optimal System and Invariant Reductions to a Nonlinear Timoshenko System

Lie symmetries and their Lie group transformations for a class of Timoshenko systems are presented. The class considered is the class of nonlinear Timoshenko systems of partial differential equations (PDEs). An optimal system of one-dimensional sub-algebras of the corresponding Lie algebra is derived. All possible invariant variables of the optimal system are obtained. The corresponding reduced systems of ordinary differential equations (ODEs) are also provided. All possible non-similar invariant conditions prescribed on invariant surfaces under symmetry transformations are given. As an application, explicit solutions of the system are given where the beam is hinged at one end and free at the other end.


Introduction
The problem of transverse vibration of beams is of importance in many engineering problems.Some of the early studies were based on the Euler-Bernoulli model, which takes into account the bending moment and the lateral displacement.Later models were based on adding shear or rotary inertia effect.In [1,2], Timoshenko proposed taking into consideration the shear, as well as the rotation effects, which proved to be suitable for non-slender beams and high frequency vibrations.Dimplekumar et al. [3] discuss the mathematical modeling for the mechanics of a solid using the distribution theory of Schwartz to the beam bending differential equations; the governing differential equations of a Timoshenko beam with jump discontinuities in slope, deflection, flexural stiffness and shear stiffness were obtained in this paper in the space of generalized functions.Mustafa and Messaoudi [4] and Fatiha [5] studied stability of the following Timoshenko system with the nonlinear frictional damping term in one equation: where the functions ϕ, ψ depend on (t, x) ∈ (0, ∞) × (0, L) and model the vertical displacement of a beam and the rotation angle of a filament, respectively.The shear angle is ψ − ϕ x , and L denotes the length of the beam.The physical parameters appearing in the system are ρ 1 , the mass density per unit length, ρ 2 , the polar moment of inertia of a cross-section, E, Young's modulus of elasticity, I, the moment of inertia of the cross-section, and k, the shear modulus.In this paper the beam is assumed to be uniform, that is the physical parameters ρ 1 , ρ 2 , EI and k are all positive constants, and χ is an arbitrary nonlinear function of the damping term ψ t .
We have solved this system by determining its symmetries and performing a Lie theoretic analysis as detailed below.
The classification of group invariant solutions of differential equations by means of the optimal systems is one of the main applications of Lie group analysis to differential equations.The main idea behind the method is discussed in detail in Ovsiannikov [6], Ibragimov [7], Olver [8] and Hydon [9].We can always construct a family of group invariant solutions corresponding to a subgroup of a symmetry group admitted by a given differential equation.Since there is an infinite number of such subgroups, it is not possible to list all of the group invariant solutions.An effective and systematic way of classifying these group invariant solutions is to obtain optimal systems of subalgebras of the Lie symmetry algebra.This leads to non-similar invariant solutions under symmetry transformations.
Recall that one calls a list { ṽα } α∈A a one-dimensional optimal system, if it satisfies the conditions: (1) completeness, i.e., any one-dimensional subalgebra is equivalent to some ṽα ; (2) inequivalence, i.e., ṽα and ṽβ are inequivalent for distinct α and β.In this paper, we have used all of the basic invariants to determine the conjugacy classes of one-dimensional subalgebras and at the same time shown that these representatives are comprehensive and mutually inequivalent.This is done by using the formula given in Section 3 of [10].This formula is more or less a direct consequence of definitions.For the convenience of the reader, a detailed explanation is given in Section 2. The idea of using invariant functions to determine conjugacy classes is also discussed in [11,12].
Here is a brief outline of the paper.In Section 2, a formula for computing invariants in the adjoint representation is given.In Section 3, Lie symmetries and their Lie group transformations for the Timoshenko system are presented.In Section 4, an optimal system of one-dimensional subalgebras of the corresponding Lie algebra is derived.In Section 5, all possible invariant variables of the optimal system are presented; moreover, the corresponding reduced systems of ordinary differential equations (ODEs) are also given.As an illustration, some invariant solutions are given explicitly by solving the reduced systems of ODEs.Furthermore, all possible non-similar invariant conditions prescribed on invariant surfaces under symmetry transformations are given.A hinged-free beam has been considered in Example 1, with a constant torque control at the hinged end and a linear force control at the free end.

Invariants in the Adjoint Representations
The following formula is stated without a proof in [10].Let X 1 , ..., X N be a basis of the Lie algebra L and ω 1 , ..., ω N be the corresponding dual basis of L * .Let X be an element of L and Y = ∑ N i=1 x i X i be a general element of L.Then, X L , the fundamental vector field corresponding to X in the adjoint representation, is given by: To see this, recall that a vector field V on an open set Ω of R N assigns to each point p of Ω a vector V(p) = (V 1 (p), ..., V N (p)).Therefore, if e 1 , ..., e N is a basis of R N and ω 1 , ..., ω N is the dual basis, then: The vector field V gives rise to the differential operator of taking the directional derivative of any function f in the direction of V(p).We identify V with the operator ∑ N j=1 ω j (V(p)) ∂ ∂x j | p .Specializing to the adjoint representation, the fundamental vector field X L is given by: Therefore, if Y = ∑ N i=1 x i X i then:

Application to Computations of Invariants
Let f be a function on L, which is invariant in the adjoint representation of L. This means that f (exp ad tX)(Y) = f (Y) for all X in the Lie algebra L and all reals t.Therefore, if we write Letting X run through the basis {X 1 , ..., X N }, we get the system: System (2) needs to be solved only for a set of generators of the Lie algebra L. The solution of this linear system of partial differential equations (PDEs) gives the invariants in the adjoint representations.
For a solvable Lie algebra L, it is best to work with the central series for the commutator L of L, which is a nilpotent subalgebra.Therefore, one can find a chain of subalgebras, each an ideal in its successor and each of codimension one.Using this chain, one is essentially dealing with invariants of a single vector field, namely if we have found the basic invariants for a subalgebras S, which is an ideal of codimension one in S, then S S operates on invariants of S, and this gives the invariants of S. Having obtained the invariants of L , one can again reach L by a series of ideals, each of codimension one, to determine the basic invariants of L. Such a series also gives a factorization of the group G = exp(Ad X) : X ∈ L in terms of one-dimensional subgroups.
The advantage of working with invariant functions is that if f is an invariant of L in the adjoint representation, then the adjoint group operates on the level sets, as well as sub-level sets of f , and one can systematically use invariants to work on level and sub-level sets of low dimensions to determine all of the conjugacy classes of one-dimensional subalgebras.This procedure will be used in Section 4 to determine one-dimensional optimal systems of the nonlinear damped Timoshenko system.This in turn will be used in Section 5 to determine optimal reduction.

Lie's Infinitesimal Method
Lie symmetry method [13,14], described extensively in the literature, is invoked in the sequel to apply Lie symmetry analysis for Timoshenko System (1).
Using the invariance condition of the system of PDEs (1): where X [2] is the second prolongation of the vector field differential operator: Comparing coefficients of the various derivatives of the dependent variables ϕ and ψ yields an over-determined linear PDE system.Carrying out the Janet basis of this over-determined system in the degree reverse lexicographical ordering as ψ > ϕ > x > t and η 2 > η 1 > ξ 2 > ξ 1 by using the command "JanetBasis" involved in the Maple package "Janet" [15] leads to the following determining equations: The command "Denominators" involved in the Maple package "Janet" returns the functions by which the Janet basis algorithm had to divide.These functions may give rise to new cases.The command "Denominators" tells us that Janet basis of the determining Equations ( 4) is produced when k = ρ 1 EI ρ 2 .However, Janet basis of the determining equations for this case is equivalent to (4).The solution to the determining Equation ( 4) is: where f (t, x) and g(t, x) satisfy the following system: Differentiating the second equation of System (6) with respect to ψ t gives: Since we are concerned with nonlinear damping term, χ (ψ t ) = 0, which implies g = g(x) and c 0 = 0. Solving System (6) yields: Thus, the Lie point symmetry generators admitted by System (1) are given by: In order to obtain the group transformations, which are generated by the resulting infinitesimal symmetry generators (9), we need to solve the following initial value problem: The one parameter group G i ( ) = e X i generated by X i for i = 1, ..., 8 is as follows: Theorem 1.If ϕ = f (t, x) and ψ = g(t, x) is a solution of Timoshenko System (1), then so is: where { i } 8 i=1 are arbitrary real numbers.
Proof.The eight-parameter group: generated by X i for i = 1, ..., 8, can be given by the composition of the transformations (11) as follows: (13)

Optimal System of One-Dimensional Sub-Algebras of the Nonlinear Damped Timoshenko System
In this section, we give the one-dimensional optimal system for the algebra with basis (9).In order to find the optimal system, one needs to classify the one-dimensional sub-algebras under the action of the adjoint representation.We follow the algorithm explained by Ibragimov [7] and Olver [8].
The non-zero commutators of the Lie algebra L 8 with basis (9) are given by: Recall that the adjoint representation is given by: The Lie algebra L 8 is solvable, and the adjoint table is given in Table 1 below: The adjoint group is generated by exp(adX) : X ∈ L .Using the solvability of L, this group consists of the elements: A = Ad(e 1 X 1 ).Ad(e 2 X 2 ).Ad(e 3 X 3 ).Ad(e 4 X 4 ).Ad(e 5 X 5 ).Ad(e 6 X 6 ).Ad(e 7 X 7 ).Ad(e 8 X 8 ).
Therefore, A is given by: Theorem 2. An optimal system of one-dimensional sub-algebras of L 8 with basis ( 9) is provided by the following generators: Proof.Let X and X be two elements in the Lie algebra L 8 with basis (9) given by X = ∑ 8 i=1 a i X i and X = ∑ 8 i=1 ãi X i .For simplicity, we will write X and X as row vectors of the coefficients on the form X = ( a 1 a 2 . . .a 8 ) and X = ( ã1 ã2 . . .ã8 ).Then, in order that X and X are in the same conjugacy class, we must have X = XA, where A is given by (15).Therefore, the theorem is proven by solving the system:  (17) for { i } 8 i=1 in terms of {a i } 8 i=1 in order to get the simplest values of { ãi } 8 i=1 .The results are presented for different cases in a tree diagram, where each of its vertices is an invariant function, and its leaves are given completely.Moreover, one can verify that all of the parameters α, β and γ appearing in each case are invariants.Therefore, the inequivalence and completeness conditions are satisfied.
All joint invariants are obtained by solving the following system of PDEs, which is given by using Formula (2): By solving this system, we obtain Φ(a 1 , a 2 , a 3 , a 4 , a 5 , a 6 , a 7 , a 8 ) = F(a 1 , a 2 , a 8 ), where Fis an arbitrary function of a 1 , a 2 , a 8 .Hence, the basic invariants of Timoshenko System (1) are a 1 , a 2 and a 8 ; this means that the first three vertices of the tree will be these invariants in any order.In our case, we will consider the order a 1 , a 2 , a 8 .The full details for each leaf are given as follows: a 2 , 7 = − a6 a 2 and 8 = − a 7 a 2 to have ã3 = ã4 = ã5 = ã6 = ã7 = 0: the conjugacy class is a 2 to have ã3 = ã5 = ã6 = ã7 = 0.The conjugacy class is EI a 4 to make ã3 = ã4 = ã7 = 0: the conjugacy class is When a 1 = 0, a 2 = 0, a 8 = 0, we need to solve system (18) again to see what are the invariants as well the next vertices of the tree.

Case 10
7 and 2 = a 6 a 7 to make ã3 = ã6 = 0: the conjugacy class is In case a 1 = a 2 = a 8 = a 7 = 0, solving system (18) yields Φ(a 1 , a 2 , a 3 , a 4 , a 5 , a 6 , a 7 , a 8 ) = F(a 4 , a 6 ) which implies, after considering a 4 and a 7 as vertices, that one can consider the invariant a 6 as a new vertex.

Optimal Reductions and Invariant Solutions
It is known that the invariant solutions for PDEs can be determined by two procedures, namely the invariant form method and the direct substitution method [14].The idea of looking for group invariant solutions generalizes quite naturally to PDEs with any number of independent and dependent variables.A one parameter group that acts nontrivially on one or more independent variables can be used to reduce the number of independent variables by one.
In this section, we focus on the invariant form method, which requires that at least one of the infinitesimals ξ 1 and ξ 2 is not equal to zero [9,14].Basically, the reason is that a vector field of the type X = η 1 ∂ ∂ϕ + η 2 ∂ ∂ψ cannot leave the submanifold of codimension two defined by the equations ϕ = f (t, x) and ψ = g(t, x) invariant, otherwise the dimension of the solution space would go down.This means the transversality condition in the sense of [16] will not be satisfied.For this reason, only the generators X 1 , ..., X 8 from system (16) are considered in this section.
We solve the invariance surface conditions explicitly by solving the corresponding characteristic equation given by: dt to get the corresponding invariants, which are used to reduce the number of independent variables by one.The procedure is explained in detail in the following example.Moreover, all possible invariant variables of the optimal system (16) and their corresponding reductions are given in Table 2.

Generators in (16) Invariant Variables
The Reduced System As another application for the optimal system based on the definition of the invariant boundary value problem given in [14,17], we classify the non-similar invariant conditions prescribed on invariant surfaces under symmetry transformations for Timoshenko system (1) as given in Table 3.This is achieved by finding invariant conditions as arbitrary functions of the invariants up to the first order and invariant surfaces as arbitrary functions of the invariants of order zero, which depend on the independent variables for each generator in the optimal system (16).Table 3. Invariant conditions prescribed on invariant surfaces and theirs reductions.
Example 1. Reduction and invariant solution using X 6 : Consider the generator X 6 = X 1 + α X 5 = ∂ ∂t − αx ∂ ∂ϕ + α ∂ ∂ψ where α ∈ R, from the optimal system (16).Solving the corresponding characteristic equations of the first prolongation: gives the following invariants up to the first order: Since the order of the Timoshenko system (1) is two, then the invariant conditions have to be functions of the first order invariants of the form: B Γ ( I 1 , I 2 , I 3 , I 4 , I 5 , I 6 , I 7 ) = 0, Γ = 1, ..., 4. (21) The invariant conditions (21) are prescribed on the invariant surfaces: which are invariants of order zero that depend on the independent variables.For α = 0, if we restrict the invariant conditions (21) to be linear, then it will take the form: where A Γ,i (x) are arbitrary functions for Γ = 1, ..., 4 and i = 1, ..., 7.
Moreover, the invariants of order zero I 1 , I 2 and I 3 give the invariant variables: The reduction of Timoshenko System (1) with the boundary conditions (23) prescribed on the surfaces x = 0 and x = L using the invariant variables (24) is the system of ODEs of the form: with general boundary conditions of the form: For instance, we consider a beam that models small motions of a hinged arm, which is hinged at the origin and free at its other end.This case was considered in [18] with two control functions.The control functions are the force f (t) applied at the free end and a torque τ(t) applied at the hinged end.The associated boundary conditions: ϕ(t, 0) = 0, ψ x (t, 0) = τ(t), ϕ x (t, L) − ψ(t, L) = f (t), ψ x (t, L) = 0, (27) can be obtained from Equation (23) for the non-zero values A 1,2 = 1 at x = 0; A 2,1 = −r, A 2,7 = 1 at x = 0; A 3,1 = −q, A 3,3 = −1, A 3,5 = 1 at x = L; A 4,7 = 1 at x = L with force f (t) = q − 2 αt and torque τ(t) = r.