Open Access
Issue
Wuhan Univ. J. Nat. Sci.
Volume 28, Number 2, April 2023
Page(s) 106 - 116
DOI https://doi.org/10.1051/wujns/2023282106
Published online 23 May 2023

© Wuhan University 2023

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

0 Introduction

Time-scale analysis originated from the work of Hilger [1], which unifies and generalizes the differential equation and difference equation theory within the framework of time-scale calculus, thus avoiding the repeated solution of the two different kinds of equations and revealing the essential differences between continuous and discrete systems [2-4]. Time scale theory has been widely applied and achieved many achievements in various fields [4-10]. In last two decades, some new advances have emerged on the study of time-scale dynamics and its symmetries, such as kinetic equations [11-13], optimal control problems [14,15], fractional variational problems [16-18], Noether theorems [18-22], Lie symmetries [23-25], Mei symmetries [26,27], canonical transformation and Hamilton-Jacobi method [28,29], time-delay dynamics [30], Herglotz variational problems [31], higher-order delta derivatives [32], etc. However, Ref. [33] indicated that, possibly affected by Bohner's original paper [11], most studies on time-scale variational problems were only shifted variational problems. In fact, the numerical calculation scheme for variational problems based on non-shifted action functionals on time scales has good properties [33,34]. Recently, three symmetries and conservation laws for non-shifted time-scales dynamic systems were studied in Refs. [25, 27, 35-37]. The study of non-shifted variational problems on time scales is a new but important research direction of analytical mechanics.

Conservation law plays an irreplaceable role in the solution and reduction of differential equations and stability analysis, so it has always been a research hotspot in the field of analytical mechanics [38-44]. Using symmetry to find conservation laws is an effective method. However, the exploration of symmetry and conservation laws theory on time scales is still very preliminary. Moreover, the studies are mainly limited to Lagrangian systems or Birkhoffian systems so far, and non-conservative forces or non-holonomic constraints are not considered. Since the phase space of mechanical system has natural symplectic structure, it is easier to describe mathematically than Lagrangian mechanics. For some cases, the symmetry that is difficult to find in configuration space can be found in phase space, and the corresponding conserved quantity can be derived by using Noether's theorem in the canonical form [45]. Here we will explore the non-shifted Hamiltonian dynamics, Noether symmetry and non-shifted conservation laws on time scales under Hamiltonian framework.

1 Non-Shifted Time-Scale Hamiltonian Dynamics

For time scale calculus and its basic operations, please refer to Refs. [2, 3, 12].

The time-scale non-shifted Hamilton action reads

S ( γ ) = t 1 t 2 [ p s ( t ) q s Δ ( t ) - H ( t , q s ( t ) , p s ( t ) ) ] Δ t (1)

where γ is a certain curve, H:Tk×Rn×RnR is the Hamiltonian on time scales, qs and ps are generalized coordinate and generalized momentum, qsΔ is the delta derivative of qs with respect to t, where s=1,2,,n. All functions belong to Crd1,Δ(T).

1.1 Hamilton System

The isochronous variational principle

δ S = 0 (2)

with endpoint conditions

δ q s | t = t 1 = δ q s | t = t 2 = 0 ,     δ p s | t = t 1 = δ p s | t = t 2 = 0 (3)

and commutative relations

δ q s Δ = ( δ q s ) Δ ,     δ p s Δ = ( δ p s ) Δ (4)

is called the time-scale non-shifted Hamilton principle.

By carrying out the variational operation of Eq. (2) and using the relations (4), we can get

- ( δ q s t 1 t H q s Δ τ ) | t 1 t 2 - ( δ p s t 1 t ( H p s - q s Δ ) Δ τ ) | t 1 t 2 + t 1 t 2 { δ q s Δ ( p s + t 1 σ ( t ) H q s Δ τ ) + δ p s Δ t 1 σ ( t ) ( H p s - q s Δ ) Δ τ } Δ t = 0 (5)

Substituting the endpoint conditions (3) into Eq. (5), we get

t 1 t 2 { δ q s Δ ( p s + t 1 σ ( t ) H q s Δ τ ) + δ p s Δ t 1 σ ( t ) ( H p s - q s Δ ) Δ τ } Δ t = 0 (6)

From the independence of δqsΔ, δpsΔ(s=1,2,,n), according to Dubois-Reymond lemma [11], we get

p s + t 1 σ ( t ) H q s Δ τ = C s , t 1 σ ( t ) ( H p s - q s Δ ) Δ τ = D s ,   s = 1,2 , , n (7)

where Cs and Ds are some constants. Taking the nabla derivative of (7), we have

p s + σ ( t ) H q s = 0 ,    σ ( t ) ( H p s - q s Δ ) = 0 (8)

So there are

q s Δ = H p s ,    p s = - σ ( t ) H q s (9)

Eqs.(9) are the Hamilton canonical equations for the non-shifted Hamilton system on time scales.

If take T=R, then σ(t)=t, Eqs.(9) become

q ˙ s = H p s ,    p ˙ s = - H q s (10)

Eqs.(10) are the classical Hamilton canonical equations.

1.2 General Holonomic Mechanical System in Phase Space

For a general holonomic mechanical system, we extend principle (2) as follows

t 1 t 2 { δ [ p s ( t ) q s Δ ( t ) - H ( t , q s ( t ) , p s ( t ) ) ] + Q s δ q s } Δ t = 0 (11)

where Qs=Qs (t,qk(t),pk(t)) are non-potential generalized forces.

Similar to the derivation of Eq. (6), from principle (11), we get

t 1 t 2 { δ q s Δ ( p s + t 1 σ ( t ) ( H q s - Q s ) Δ τ ) + δ p s Δ t 1 σ ( t ) ( H p s - q s Δ ) Δ τ } Δ t = 0 (12)

From the independence of δqsΔ, δpsΔ(s=1,2,,n), according to Dubois-Reymond lemma [11], we get

p s + t 1 σ ( t ) ( H q s - Q s ) Δ τ = C s ' , t 1 σ ( t ) ( H p s - q s Δ ) Δ τ = D s ' (13)

where Cs' and Ds' are some constants. Hence, we have

q s Δ = H p s ,    p s = - σ ( t ) H q s + σ ( t ) Q s (14)

Eqs.(14) are the time-scale dynamic equations of the general holonomic system. When Qs=0, Eqs. (14) are reduced to Eqs. (9), which are the non-shifted Hamilton canonical equations.

1.3 Nonholonomic Mechanical System in Phase Space

Consider the system is subject to g bilateral ideal nonholonomic constraints

f β ( t , q s , q s Δ ) = 0 (15)

and the virtual displacements δqs need to meet the conditions

F β s ( t , q j , q j Δ ) δ q s = 0 (16)

where β=1,2,,g. In general, Fβs is independent of fβqsΔ, and the constraints are of non-Chetaev. If Fβs=fβqsΔ, then the constraints are of Chetaev.

If the non-shifted Lagrangian is L=L(t,qs,qsΔ), then

p s = L q s Δ (17)

According to Eq. (17), qsΔ=qsΔ(t,qj,pj) can be solved, and then substituted into Eqs. (15) and (16), thus the constraints (15) and restriction conditions (16) can be written as

f ˜ β ( t , q s , p s ) = 0 (18)

F ˜ β s ( t , q j , p j ) δ q s = 0 (19)

By introducing the constraint multiplier λβ multiplied by each of Eqs. (19) and summing over β, and integrating the equation on the interval [t1,t2], and by integration by parts, we get

( δ q s t 1 t λ β F ˜ β s Δ τ ) | t 1 t 2 - t 1 t 2 δ q s Δ ( t 1 σ ( t ) λ β F ˜ β s Δ τ ) Δ t = 0 (20)

Taking into account conditions (3), we get

- t 1 t 2 δ q s Δ ( t 1 σ ( t ) λ β F ˜ β s Δ τ ) Δ t = 0 (21)

Add Eq. (21) to Eq. (12), and we get

t 1 t 2 { δ q s Δ ( p s + t 1 σ ( t ) ( H q s - Q s - λ β F ˜ β s ) Δ τ ) + δ p s Δ t 1 σ ( t ) ( H p s - q s Δ ) Δ τ } Δ t = 0 (22)

According to the Lagrange multiplier method, without loss of generality, choose the multiplier λβ such that Cβ=0(β=1,2,,g), and using Dubois-Reymond lemma [11], from Eq. (22), we get

p s + t 1 σ ( t ) ( H q s - Q s - λ β F ˜ β s ) Δ τ = C s ,     t 1 σ ( t ) ( H p s - q s Δ ) Δ τ = D s (23)

where Cg+1,,Cn, Ds are some constants. So there are

q s Δ = H p s ,     p s = - σ ( t ) H q s + σ ( t ) Q s + σ ( t ) λ β F ˜ β s (24)

Assuming that the system is non-singular, by using Eqs. (24) and (18), λβ can be solved as the function of qs, ps and t. Therefore, Eqs.(24) can be expressed as

q s Δ = H p s ,     p s = - σ ( t ) H q s + σ ( t ) ( Q s + Λ s ) (25)

where Λs=λβF˜βs are the constraint forces corresponding to the nonholonomic constraints (18). Eqs. (25) can be regarded as a holonomic system corresponding to the nonholonomic system determined by Eqs.(18) and (24). If the initial values of qs and ps satisfy Eq. (18), namely

f ˜ β ( t 0 , q s 0 , p s 0 ) = 0 (26)

then the solution of (25) is the desired solution of time-scale nonholonomic systems (18) and (24).

2 Noether Symmetry for Time-Scale Hamiltonian Dynamics

2.1 Noether Symmetry for Hamilton System

The infinitesimal transformations are

t ¯ = t + ε ξ 0 ( t , q j , p j ) ,    q ¯ s ( t ¯ ) = q s ( t ) + ε ξ s ( t , q j , p j ) ,    p ¯ s ( t ¯ ) = p s ( t ) + ε η s ( t , q j , p j ) (27)

where ξ0, ξs and ηs are the generating functions, ε is the infinitesimal parameter, and s,j=1,2,,n. Let the map tα(t)=t+εξ0+o(ε) be a strictly increasing Crd1,Δ(T) function, whose image is denoted as T, delta derivative is Δ¯, forward jump operator σ¯, and σ¯α=ασ.

Under the transformation (27), the Hamilton action (1) reads

S ( γ ¯ ) = α ( t 1 ) α ( t 2 ) [ p ¯ s ( t ¯ ) q ¯ s Δ ¯ ( t ¯ ) - H ( t ¯ , q ¯ s ( t ¯ ) , p ¯ s ( t ¯ ) ) ] Δ ¯ t ¯    = t 1 t 2 [ ( p s + ε η s ) q s Δ + ε ξ s Δ 1 + ε ξ 0 Δ - H ( t + ε ξ 0 , q s + ε ξ s , p s + ε η s ) ] ( 1 + ε ξ 0 Δ ) Δ t = t 1 t 2 { p s q s Δ - H ( t , q s , p s ) + ε ( p s ξ s Δ + q s Δ η s ) - ε ( H t ξ 0 + H q s ξ s + H p s η s + H ξ 0 Δ ) + o ( ε ) } Δ t (28)

Therefore, the nonisochronous variation Δ*S, namely the main-line part of difference S(γ¯)-S(γ) relative to ε, is

Δ * S = t 1 t 2 ε [ p s ξ s Δ + q s Δ η s - ( H t ξ 0 + H q s ξ s + H p s η s + H ξ 0 Δ ) ] Δ t (29)

By straightforward calculation, formula (29) can also be expressed as

Δ * S = t 1 t 2 ε { t ( p s ξ s σ - H ξ 0 σ ) + ( H - σ ( t ) H t ) ξ 0 + σ ( t ) ( q s Δ - H p s ) η s - ( p s + σ ( t ) H q s ) ξ s } Δ t σ ( t ) (30)

Formulas (29) and (30) are two mutually equivalent nonisochronous variational formulas of non-shifted Hamilton action on time scales.

Definition 1   For the time-scale non-shifted Hamilton system (9), the transformation (27) is said to be Noether symmetric, if and only if Δ*S=0.

By using Eqs. (29) and (30), we obtain:

Criterion 1 For the time-scale non-shifted Hamilton system (9), if the Noether identity

p s ξ s Δ + q s Δ η s - H t ξ 0 - H q s ξ s - H p s η s - H ξ 0 Δ = 0 (31)

holds, then the transformation (27) is Noether symmetric.

Criterion 2 If the generating functions ξ0, ξs and ηs solve the equation

t ( p s ξ s σ - H ξ 0 σ ) + ( H - σ ( t ) H t ) ξ 0 + σ ( t ) ( q s Δ - H p s ) η s - ( p s + σ ( t ) H q s ) ξ s = 0 (32)

then the transformation (27) is Noether symmetric for the time-scale Hamilton system (9).

Assume that H1 is another Hamiltonian on time scales, if, considering the first-order approximation, the transformation (27) satisfies the following relation

t 1 t 2 [ p s ( t ) q s Δ ( t ) - H ( t , q s ( t ) , p s ( t ) ) ] Δ t = α ( t 1 ) α ( t 2 ) [ p ¯ s ( t ¯ ) q ¯ s Δ ¯ ( t ¯ ) - H 1 ( t ¯ , q ¯ s ( t ¯ ) , p ¯ s ( t ¯ ) ) ] Δ ¯ t ¯ (33)

then the action (1) is called quasi-invariant. Both H1 and H satisfy the same time-scale Hamilton equations. So there is

H 1 ( t , q s ( t ) , p s ( t ) ) = H ( t , q s ( t ) , p s ( t ) ) - ε G Δ (34)

where G=G(t,qs,ps) is the gauge function.

Definition 2   For the time-scale non-shifted Hamilton system (9), the transformation (27) is said to be Noether quasi-symmetric, if and only if

Δ * S = - t 1 t 2 ε G Δ Δ t (35)

By using Eqs.(29) and (30), we obtain:

Criterion 3 For the time-scale non-shifted Hamilton system (9), if the generalized Noether identity

p s ξ s Δ + q s Δ η s - H t ξ 0 - H q s ξ s - H p s η s - H ξ 0 Δ + G Δ = 0 (36)

holds, then the transformation (27) is Noether quasi-symmetric.

Criterion 4 If the generating functions ξ0, ξs and ηs solve the equation

t ( p s ξ s σ - H ξ 0 σ + G σ ) + ( H - σ ( t ) H t ) ξ 0 + σ ( t ) ( q s Δ - H p s ) η s - ( p s + σ ( t ) H q s ) ξ s = 0 (37)

then the transformation (27) is Noether quasi-symmetric for the time-scale Hamilton system (9).

2.2 Noether Symmetry for General Holonomic Mechanical System in Phase Space

For the general holonomic system, if the following relation

t 1 t 2 [ p s ( t ) q s Δ ( t ) - H ( t , q s ( t ) , p s ( t ) ) ] Δ t = α ( t 1 ) α ( t 2 ) [ p ¯ s ( t ¯ ) q ¯ s Δ ¯ ( t ¯ ) - H 1 ( t ¯ , q ¯ s ( t ¯ ) , p ¯ s ( t ¯ ) ) ] Δ ¯ t ¯ + t 1 t 2 Q s δ q s Δ t (38)

is satisfied, then the action (1) is called generalized quasi-invariant.

Definition 3   For the time-scale non-shifted general holonomic system (14), the transformation (27) is said to be generalized Noether quasi-symmetric, if and only if

Δ * S + t 1 t 2 Q s δ q s Δ t = - t 1 t 2 ε G Δ Δ t (39)

Criterion 5 For the time-scale non-shifted general holonomic system (14), if the generalized Noether identity

p s ξ s Δ + q s Δ η s - H t ξ 0 - H q s ξ s - H p s η s - H ξ 0 Δ + Q s   ( ξ s - q s Δ ξ 0 ) + G Δ = 0 (40)

holds, then the transformation (27) is generalized Noether quasi-symmetric.

Criterion 6 If the generating functions ξ0, ξs and ηs solve the equation

t ( p s ξ s σ - H ξ 0 σ + G σ ) + ( H - σ ( t ) H t - σ ( t ) Q s q s Δ ) ξ 0 + σ ( t ) ( q s Δ - H p s ) η s - ( p s + σ ( t ) H q s - σ ( t ) Q s ) ξ s = 0 (41)

then the transformation (27) is generalized Noether quasi-symmetric for the time-scale general holonomic system (14).

2.3 Noether Symmetry for Nonholonomic Mechanical System in Phase Space

For the nonholonomic system, we have

Definition 4   For the corresponding holonomic system (25), the transformation (27) is said to be generalized Noether quasi-symmetric, if and only if

Δ * S + t 1 t 2 ( Q s + Λ s ) δ q s Δ t = - t 1 t 2 ε G Δ Δ t (42)

Criterion 7 For the corresponding holonomic system (25), if the generalized Noether identity

p s ξ s Δ + q s Δ η s - H t ξ 0 - H q s ξ s - H p s η s - H ξ 0 Δ + ( Q s + Λ s ) ( ξ s - q s Δ ξ 0 ) + G Δ = 0 (43)

holds, then the transformation (27) is generalized Noether quasi-symmetric.

Criterion 8 If the generating functions ξ0, ξs and ηs solve the equation

t ( p s ξ s σ - H ξ 0 σ + G σ ) + ( H - σ ( t ) H t - σ ( t ) ( Q s + Λ s ) q s Δ ) ξ 0 + σ ( t ) ( q s Δ - H p s ) η s - ( p s + σ ( t ) H q s - σ ( t ) ( Q s + Λ s ) ) ξ s = 0 (44)

then the transformation (27) is generalized Noether quasi-symmetric for the corresponding holonomic system (25).

The restriction conditions of Eqs. (19) on the generating functions are

F ˜ β s ( ξ s - q s Δ ξ 0 ) = 0 (45)

Definition 5   For the time-scale non-shifted nonholonomic system determined by (18) and (24), if and only if the formula (42) and restriction conditions (45) hold, then the transformation (27) is said to be generalized Noether quasi-symmetric.

Criterion 9 For the time-scale non-shifted nonholonomic system determined by (18) and (24), if the generalized Noether identity (43) and the restriction conditions (45) hold, then the transformation (27) is generalized Noether quasi-symmetric.

Criterion 10 If the generating functions ξ0, ξs and ηs solve the equation (44) and the restriction conditions (45), then the transformation (27) is generalized Noether quasi-symmetric for the time-scale nonholonomic system determined by (18) and (24).

3 Noether Theorems for Time-Scales Hamiltonian Dynamics

Noether symmetry is closely related to conservation laws. Here we establish and prove Noether's theorems for time-scale non-shifted holonomic and nonholonomic Hamiltonian dynamics.

3.1 Noether Theorems for Hamilton System

Theorem 1   For the time-scale non-shifted Hamilton system (9), if the transformation (27) is Noether symmetric, then

I N = p s ξ s σ - H ξ 0 σ + t 1 t { [ H - σ ( t ) H t ] ξ 0 } t (46)

is a non-shifted Noether conserved quantity.

Proof   Due to

t I N = p s ξ s + σ ( t ) p s ξ s Δ - H ξ 0 - σ ( t ) H ξ 0 Δ + ( H - σ ( t ) H t ) ξ 0 = σ ( t ) ( p s ξ s Δ + q s Δ η s - H t ξ 0 - H q s ξ s - H p s η s - H ξ 0 Δ ) + σ ( t ) ( - q s Δ + H p s ) η s + ( p s + σ ( t ) H q s ) ξ s (47)

Substituting the non-shifted Hamilton equations (9) and the Noether identity (31) into (47), we get

t I N = 0 (48)

Therefore, formula (46) is a non-shifted Noether conserved quantity.

Theorem 2   For the time-scale non-shifted Hamilton system (9), if the transformation (27) is Noether quasi-symmetric, then

I N = p s ξ s σ - H ξ 0 σ + G σ + t 1 t { [ H - σ ( t ) H t ] ξ 0 } t (49)

is a non-shifted Noether conserved quantity.

Proof   Taking the nabla derivative of (49), and using Eqs.(9) and (36), we get the result immediately.

Theorem 1   and 2 are Noether's theorems for time-scale non-shifted Hamilton system.

3.2 Noether Theorems for General Holonomic Mechanical System in Phase Space

For the general holonomic mechanical system, we have

Theorem 3   For the time-scale non-shifted general holonomic system (14), if the transformation (27) is generalized Noether quasi-symmetric, then

I N = p s ξ s σ - H ξ 0 σ + G σ + t 1 t { [ H - σ ( t ) ( H t + Q s q s Δ ) ] ξ 0 } t (50)

is a non-shifted Noether conserved quantity.

Proof   Due to

t I N = p s ξ s + σ ( t ) p s ξ s Δ - H ξ 0 - σ ( t ) H ξ 0 Δ + ( G σ ) + ( H - σ ( t ) H t - σ ( t ) Q s q s Δ ) ξ 0 = σ ( t ) ( p s ξ s Δ + q s Δ η s - H t ξ 0 - H q s ξ s - H p s η s - H ξ 0 Δ + Q s ( ξ s - q s Δ ξ 0 ) + G Δ ) + σ ( t ) ( - q s Δ + H p s ) η s + ( p s + σ ( t ) H q s - σ ( t ) Q s ) ξ s = 0 (51)

This completes the proof.

Theorem 3   is Noether's theorem for time-scale non-shifted general holonomic system under Hamiltonian framework.

3.3 Noether Theorems for Nonholonomic Mechanical System in Phase Space

For the nonholonomic mechanical system, we have

Theorem 4   For the corresponding holonomic system (25), if the transformation (27) is generalized Noether quasi-symmetric, then

I N = p s ξ s σ - H ξ 0 σ + G σ + t 1 t { [ H - σ ( t ) ( H t + Q s q s Δ + Λ s q s Δ ) ] ξ 0 } t (52)

is a non-shifted Noether conserved quantity.

Proof   Taking the nabla derivative of (52), and using Eqs.(25) and (43), we get the results.

Theorem 5   For the time-scale non-shifted nonholonomic system determined by (18) and (24), if the transformation (27) is generalized Noether quasi-symmetric, then formula (52) is a non-shifted Noether conserved quantity.

Theorem 5   and Theorem 4 are Noether's theorems for time-scale non-shifted nonholonomic system and its corresponding holonomic system under Hamiltonian framework.

4 Examples

Example 1 Consider a non-shifted holonomic non-conservative system on time scale T={2m:mN0}, and let the Lagrangian function be

L = 1 2 [ ( q 1 Δ ) 2 + ( q 2 Δ ) 2 ] - 1 2 q 2 2 (53)

The generalized forces are

Q 1 =    - q 2 Δ ,     Q 2 =   0 (54)

The generalized momenta and the Hamiltonian are

p 1 = L q 1 Δ = q 1 Δ , p 2 = L q 2 Δ = q 2 Δ , H = 1 2 ( p 1 2 + p 2 2 ) + 1 2 q 2 2 (55)

According to equation (14), the Hamilton equations are

q 1 Δ = p 1 , q 2 Δ = p 2 , p 1 = - σ ( t ) p 2 , p 2 = - σ ( t ) q 2 (56)

If we take T=R, then Eqs.(56) are reduced to

q ˙ 1 = p 1 , q ˙ 2 = p 2 , p ˙ 1 + p 2 = 0 , p ˙ 2 + q 2 = 0 (57)

This is the classic Hojman-Urrutia problem [46].

The generalized Noether identity (40) reads

p 1 ξ 1 Δ + p 2 ξ 2 Δ + η 1 q 1 Δ + η 2 q 2 Δ - q 2 ξ 2 - p 1 η 1 - p 2 η 2 - 1 2 ( p 1 2 + p 2 2 ) ξ 0 Δ - 1 2 q 2 2 ξ 0 Δ - p 2 ( ξ 1 - q 1 Δ ξ 0 ) + G Δ = 0 (58)

Since σ(t)=2t, Eq.(58) has a solution

ξ 0 1 = 0 , ξ 1 1 = 1 , ξ 2 1 = 0 , η 1 1 = 1 , η 2 1 = 0 , G 1 = q 2 (59)

ξ 0 2 = 0 , ξ 1 2 = 2 t , ξ 2 2 = 1 , η 1 2 = 0 , η 2 2 = 0 , G 2 = - 2 q 1 + q 2 t (60)

By Theorem 3, we obtain

I N 1 = p 1 + q 2 σ = c o n s t . (61)

I N 2 = 4 p 1 t + p 2 + 2 q 2 σ t - 2 q 1 σ = c o n s t . (62)

The conserved quantity (61) and (62) correspond to the quasi-symmetry (59) and (60). Assume the initial conditions q1(0)=1, q2(0)=1, p1(0)=1, p2(0)=1. Let m[0,6]. The numerical simulation results of conserved quantities (61) and (62) on t[2,64] are shown in Fig. 1.

thumbnail Fig. 1

The values of conserved quantities (61) and (62) on t[2,64]

From Fig. 1, it can be seen intuitively that conserved quantities obtained from formulae (61) and (62) are both constants, which shows the correctness of Theorem 3.

Example 2 Let us study Appell-Hamel problem [47] on time scales. The non-shifted Lagrangian and nonholonomic constraint are respectively

L = 1 2 m [ ( q 1 Δ ) 2 + ( q 2 Δ ) 2 + ( q 3 Δ ) 2 ] - m g q 3 (63)

f = q 3 Δ - [ ( q 1 Δ ) 2 + ( q 2 Δ ) 2 ] 1 / 2 = 0 (64)

The constraint (64) is of Chetaev type. The generalized momenta and the Hamiltonian are

p 1 = L q 1 Δ = m q 1 Δ , p 2 = L q 2 Δ = m q 2 Δ , p 3 = L q 3 Δ = m q 3 Δ (65)

H = 1 2 m ( p 1 2 + p 2 2 + p 3 2 ) + m g q 3 (66)

In canonical coordinates, the constraint (64) can be shown as

f ˜ = p 3 - ( p 1 2 + p 2 2 ) 1 / 2 = 0 (67)

The time-scale dynamical equations can be expressed as

q 1 Δ = 1 m p 1 , q 2 Δ = 1 m p 2 , q 3 Δ = 1 m p 3 , p 1 = - σ ( t ) λ p 1 ( p 1 2 + p 2 2 ) - 1 / 2 ,   p 2 = - σ ( t ) λ p 2 ( p 1 2 + p 2 2 ) - 1 / 2 , p 3 = σ ( t ) ( - m g + λ ) (68)

Taking the nabla derivative of (67), we get

p 3 - 1 2 ( p 1 2 + p 2 2 ) - 1 / 2 ( p 1 + p 1 ρ ) p 1 - 1 2 ( ( p 1 ρ ) 2 + p 2 2 ) - 1 / 2 ( p 2 + p 2 ρ ) p 2 = 0 (69)

Substituting (68) into (69), we can get

λ = 2 m g ( p 1 2 + p 2 2 ) 1 / 2 2 ( p 1 2 + p 2 2 ) 1 / 2 + p 1 ( p 1 + p 1 ρ ) ( p 1 2 + p 2 2 ) - 1 / 2 + p 2 ( p 2 + p 2 ρ ) ( ( p 1 ρ ) 2 + p 2 2 ) - 1 / 2 (70)

Therefore, the nonholonomic constraint forces are

Λ 1 = - λ p 1 ( p 1 2 + p 2 2 ) - 1 / 2 , Λ 2 = - λ p 2 ( p 1 2 + p 2 2 ) - 1 / 2 , Λ 3 = λ (71)

From (43), the generalized Noether identity for the system is

p 1 ξ 1 Δ + p 2 ξ 2 Δ + p 3 ξ 3 Δ + η 1 q 1 Δ + η 2 q 2 Δ + η 3 q 3 Δ - m g ξ 3 - 1 m p 1 η 1 - 1 m p 2 η 2 - 1 m p 3 η 3

- 1 2 m ( p 1 2 + p 2 2 + p 3 2 ) ξ 0 Δ - m g q 3 ξ 0 Δ + Λ 1 ( ξ 1 - q 1 Δ ξ 0 ) + Λ 2 ( ξ 2 - q 2 Δ ξ 0 ) + Λ 3 ( ξ 3 - q 3 Δ ξ 0 ) + G Δ = 0 (72)

The restriction condition (45) reads

- p 1 ( p 1 2 + p 2 2 ) - 1 / 2 ( ξ 1 - q 1 Δ ξ 0 ) - p 2 ( p 1 2 + p 2 2 ) - 1 / 2 ( ξ 2 - q 2 Δ ξ 0 ) + ξ 3 - q 3 Δ ξ 0 = 0 (73)

The equations (72) and (73) have a solution

ξ 0 1 = - 1 , ξ 1 1 = 0 , ξ 2 1 = 0 , ξ 3 1 = 0 , η 1 1 = 0 , η 2 1 = 0 , η 3 1 = 0 , G 1 = 0 (74)

From Theorem 5, we get the Noether conserved quantity

I N = 1 2 m ( p 1 2 + p 2 2 + p 3 2 ) + m g q 3 + t 1 t { λ σ 2 m [ ( p 1 p 1 ρ + p 2 p 2 ρ ) ( p 1 2 + p 2 2 ) - 1 / 2 - p 3 ρ ] + 1 2 σ g ( p 3 + p 3 ρ ) - m g q 3 Δ } t = c o n s t . (75)

If T=R, then σ(t)=t, and formula (75) becomes

I N = 1 2 m ( p 1 2 + p 2 2 + p 3 2 ) + m g q 3 = c o n s t . (76)

This is the classical conservation law of energy, which has been given in Ref.[47].

If T=Z, then σ(t)=t+1, ρ(t)=t-1, μ(t)=1, and formula (75) becomes

I N = 1 2 m [ p 1 2 ( t ) + p 2 2 ( t ) + p 3 2 ( t ) ] + m g q 3 ( t ) + τ = t 1 + 1 t { λ ( τ ) 2 m [ ( p 1 ( τ ) p 1 ( τ - 1 ) + p 2 ( τ ) p 2 ( τ - 1 ) ) ( p 1 2 ( τ ) + p 2 2 ( τ ) ) - 1 / 2 - p 3 ( τ - 1 ) ] + 1 2 g [ p 3 ( τ ) + p 3 ( τ - 1 ) ] - m g [ q 3 ( τ + 1 ) - q 3 ( τ ) ] } (77)

This is the Noether conserved quantity of the discrete version with unit step size.

5 Conclusion

The time-scale calculus provides an excellent mathematical tool for exploring the dynamics of continuous and discrete systems or their mixtures, and has attracted extensive attentions. In this paper, we proposed the time-scale non-shifted Hamilton principle and extended it to non-conservative systems, and derived the dynamic equations for non-shifted Hamilton systems, non-shifted general holonomic systems and non-shifted nonholonomic systems. We defined Noether symmetries and gave their criteria. We proved Noether's theorems for non-shifted Hamilton systems, non-shifted general holonomic systems and non-shifted nonholonomic systems, and obtained the non-shifted Noether conserved quantities. The ideas presented here can be applied to solving the symmetries of complex dynamics under time-scale framework, such as nonlinear problems.

References

  1. Hilger S. Analysis on measure chains — A unified approach to continuous and discrete calculus [J]. Results in Mathematics, 1990, 18: 18-56. [CrossRef] [Google Scholar]
  2. Bohner M, Peterson A. Dynamic Equations on Time Scales: An Introduction with Applications [M]. Boston: Birkhäuser, 2001. [Google Scholar]
  3. Bohner M, Georgiev S G. Multivariable Dynamic Calculus on Time Scales [M]. Berlin: Springer Verlag, 2016. [CrossRef] [Google Scholar]
  4. Bohner M, Peterson A. Advances in Dynamic Equations on Time Scales [M]. Boston: Birkhäuser, 2003. [CrossRef] [Google Scholar]
  5. Agarwal R, Bohner M, O'Regan D, et al. Dynamic equations on time scales: A survey [J]. Journal of Computational and Applied Mathematics, 2002, 141(1-2): 1-26. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  6. Atici F M, Biles D C, Lebedinsky A. An application of time scales to economics [J]. Mathematical and Computer Modelling, 2006, 43(7-8): 718-726. [Google Scholar]
  7. Bohner M, Fan M, Zhang J. Periodicity of scalar dynamic equations and applications to population models [J]. Journal of Mathematical Analysis and Applications, 2007, 330(1): 1-9. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  8. Ferreira R A C, Torres D F M. Remarks on the calculus of variations on time scales [J]. International Journal of Ecological Economics and Statistics, 2007, 9(F07): 65-73. [Google Scholar]
  9. Benkhettou N, Brito da Cruz A M C, Torres D F M. A fractional calculus on arbitrary time scales: Fractional differentiation and fractional integration [J]. Signal Processing, 2015, 107: 230-237. [CrossRef] [Google Scholar]
  10. Dryl M, Torres D F M. Direct and inverse variational problems on time scales: A survey [C]// Modeling, Dynamics, Optimization and Bioeconomics II. Proceedings in Mathematics and Statistics. Berlin: Springer-Verlag, 2017, 195: 223-265. [MathSciNet] [Google Scholar]
  11. Bohner M. Calculus of variations on time scales [J]. Dynamic Systems and Application, 2004, 13: 339-349. [Google Scholar]
  12. Martins N, Torres D F M. Calculus of variations on time scales with nabla derivatives [J]. Nonlinear Analysis: Theory, Methods and Application, 2019, 71(12): e763-e773. [Google Scholar]
  13. Jin S X, Zhang Y. Generalized Chaplygin equations for nonholonomic systems on time scales [J]. Chinese Physics B, 2018, 27(2): 020502. [Google Scholar]
  14. Hilscher R, Zeidan V. First order conditions for generalized variational problems over time scales [J]. Computers and Mathematics with Applications, 2011, 62(9): 3490-3503. [CrossRef] [MathSciNet] [Google Scholar]
  15. Girejko E, Malinowska A B, Torres D F M. Delta-nabla optimal control problems [J]. Journal of Vibration and Control, 2011, 17(11): 1634-1643. [Google Scholar]
  16. Bastos N R O, Ferreira R A C, Torres D F M. Discrete-time fractional variational problems [J]. Signal Processing, 2011, 91(3): 513-524. [CrossRef] [Google Scholar]
  17. Mekhalfi K, Torres D F M. Generalized fractional operators on time scales with application to dynamic equations [J]. The European Physical Journal Special Topics, 2017, 226(16-18): 3489-3499. [NASA ADS] [CrossRef] [Google Scholar]
  18. Tian X, Zhang Y. Fractional time-scales Noether theorem with Caputo Formula derivatives for Hamiltonian systems [J]. Applied Mathematics and Computation, 2021, 393: 125753. [Google Scholar]
  19. Bartosiewicz Z, Torres D F M. Noether's theorem on time scales [J]. Journal of Mathematical Analysis and Applications, 2008, 342(2): 1220-1226. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  20. Cai P P, Fu J L, Guo Y X. Noether symmetries of the nonconservative and nonholonomic systems on time scales [J]. Science China Physics, Mechanics and Astronomy, 2013, 56(4): 1017-1028. [CrossRef] [Google Scholar]
  21. Song C J, Zhang Y. Noether theorem for Birkhoffian systems on time scales [J]. Journal of Mathematical Physics, 2015, 56(10): 102701. [CrossRef] [MathSciNet] [Google Scholar]
  22. Anerot B, Cresson J, Belgacem K H, et al. Noether's-type theorems on time scales [J]. Journal of Mathematical Physics, 2020, 61(11): 113502. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  23. Zhai X H, Zhang Y. Lie symmetry analysis on time scales and its application on mechanical systems [J]. Journal of Vibration and Control, 2019, 25(3): 581-592. [CrossRef] [MathSciNet] [Google Scholar]
  24. Zhang Y. Adiabatic invariants and Lie symmetries on time scales for nonholonomic systems of non-Chetaev type [J]. Acta Mechamica, 2020, 128(1): 293-303. [CrossRef] [Google Scholar]
  25. Zhang Y, Tian X, Zhai X H, et al. Hojman conserved quantity for time scales Lagrange systems [J]. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(10): 2814-2822(Ch). [Google Scholar]
  26. Zhai X H, Zhang Y. Mei symmetry for time-scales Euler-Lagrange equations and its relation to Noether symmetry [J]. Acta Physica Polonica A, 2019, 136(3): 439-443. [CrossRef] [MathSciNet] [Google Scholar]
  27. Zhang Y. Mei's symmetry theorems for time scales nonshifed mechanical systems [J]. Theoretical and Applied Mechanics Letters, 2021, 11(5): 100286. [CrossRef] [Google Scholar]
  28. Zhai X H, Zhang Y. Hamilton-Jacobi method for mechanical systems on time scales [J]. Complexity, 2018, 2018: 8070658. [Google Scholar]
  29. Zhang Y, Zhai X H. Generalized canonical transformation for second-order Birkhoffian systems on time scales [J]. Theoretical and Applied Mechanics Letters, 2019, 9(6): 353-357. [CrossRef] [Google Scholar]
  30. Zhai X H, Zhang Y. Noether theorem for non-conservative systems with time delay on time scales [J]. Communications in Nonlinear Science and Numerical Simulation, 2017, 52: 32-43. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  31. Tian X, Zhang Y. Noether symmetry and conserved quantity for Hamiltonian system of Herglotz type on time scales [J]. Acta Mechanica, 2018, 229(9): 3601-3611. [CrossRef] [MathSciNet] [Google Scholar]
  32. Ferreira R A C, Malinowska A B, Torres D F M. Optimality conditions for the calculus of variations with higher-order delta derivatives [J]. Applied Mathematics Letters, 2011, 24(1): 87-92. [Google Scholar]
  33. Bourdin L. Nonshifted calculus of variations on time scales with Formula Formula [J]. Journal of Mathematical Analysis and Applications, 2014, 411(2): 543-554. [CrossRef] [MathSciNet] [Google Scholar]
  34. Cresson J, Malinowska A B, Torres D F M. Time scale differential, integral, and variational embeddings of Lagrangian systems [J]. Computers & Mathematics with Applications, 2012, 64(7): 2294-2301. [Google Scholar]
  35. Song C J, Cheng Y. Noether's theorems for nonshifted dynamic systems on time scales [J]. Applied Mathematics and Computation, 2020, 374: 125086. [CrossRef] [MathSciNet] [Google Scholar]
  36. Zhang Y. Nonshifted dynamics of constrained systems on time scales under Lagrange framework and its Noether's theorem [J]. Communications in Nonlinear Science and Numerical Simulation, 2022, 108: 106214. [Google Scholar]
  37. Zhang Y. Mei's symmetry theorems for non-migrated Birkhoffian systems on a time scale [J]. Acta Physica Sinica, 2021, 70 (24): 244501(Ch). [CrossRef] [Google Scholar]
  38. Mei F X. Analytical Mechanics II [M]. Beijing: Beijing Institute of Technology Press, 2013. [Google Scholar]
  39. José J V, Saletan E J. Classical Dynamics: A Contemporary Approach [M]. Cambridge: Cambridge University Press, 1998. [Google Scholar]
  40. Mei F X, Zhang Y F, Shang M. Lie symmetries and conserved quantities of Birkhoffian systems [J]. Mechanics Research Communications, 1999, 26(1): 7-12. [CrossRef] [MathSciNet] [Google Scholar]
  41. Fu J L, Fu L P, Chen B Y, et al. Lie symmetries and their inverse problems of nonholonomic Hamilton systems with fractional derivatives [J]. Physics Letters A, 2016, 380(1-2): 15-21. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  42. Fu J L, Chen L Q. Form invariance, Noether symmetry and Lie symmetry of Hamiltonian systems in phase space [J]. Mechanics Research Communications, 2004, 31(1): 9-19. [CrossRef] [MathSciNet] [Google Scholar]
  43. Chen X W, Mei F X. Constrained mechanical systems and gradient systems with strong Lyapunov functions [J]. Mechanics Research Communications, 2016, 76: 91-95. [Google Scholar]
  44. Zhang Y. Noether's symmetry and conserved quantity for a time-delayed Hamiltonian system of Herglotz type [J]. Royal Society Open Science, 2018, 5(10): 180208. [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  45. Mei F X, Wu H B, Li Y M. A Brief History of Analytical Mechanics [M]. Beijing: Science Press, 2019(Ch). [Google Scholar]
  46. Santilli R M. Foundations of Theoretical Mechanics II [M]. New York: Springer-Verlag, 1983. [CrossRef] [MathSciNet] [Google Scholar]
  47. Mei F X. Applications of Lie Groups and Lie Algebras to Constrained Mechanical Systems [M]. Beijing: Science Press, 1999(Ch). [Google Scholar]

All Figures

thumbnail Fig. 1

The values of conserved quantities (61) and (62) on t[2,64]

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.