Abstract

In this paper a system derived by an optimal control problem for the ball-plate dynamics is considered. Symplectic and Lagrangian realizations are given and some symmetries are studied. The image of the energy-Casimir mapping is described and some connections with the dynamics of the considered system are presented.

1. Introduction

In 1983 the problem of optimal rolling of a sphere over a horizontal plane without slipping was posed by Hammersley [1]. Using quaternion calculus of variations and optimal control theory, Arthurs and Walsh [2] obtained a system of differential equations that can be integrated analytically in terms of elliptic integrals.

A similar problem, named the ball-plate problem, was considered by Brockett and Dai [3]. They have considered the problem of optimal rolling of a ball without slipping between two horizontal plates, the lower plate fixed, and the upper plate mobile. More precisely, by moving of the upper plate, the ball is rolling from the prescribed initial state into some terminal state along a path which minimizes among all possible path which satisfy these boundary conditions, where is the velocity of the moving plate, is the time of transfer, and is a positive constant. The control is the velocity of the center of the ball. The state of such a system is described by the position of the center of the ball and by the orientation of the ball in the three-dimensional space. The nonslipping assumptions mean that the instantaneous velocity of the point of contact of the ball with the lower plate is equal to zero, and the velocity of the point of the contact with the moving plate is equal to the velocity of the moving plate. Such kind of motion is of great interest in mechanics, robotics, and control theory [47], making it a widely investigated problem.

The ball-plate problem has been formulated by Jurdjevic [8] as an optimal control problem on the Lie group . Among the results obtained by Jurdjevic in his study we state here the equations of motion, their integrals of motion, and the integrability properties. Hamilton-Poisson formulations, a study of stability of the equilibrium points and numerical integration of the system obtained by Jurdjevic were presented in [9, 10]. We mention other topics regarding the problem of optimal rolling sphere such as exponential map and the corresponding group of symmetries, Maxwell points [11], parametrization of extremal trajectories by using elliptic functions, and asymptotics of extremal trajectories [12, 13].

According to the Jurdjevic’s paper [8], in the following, the derivation of the equations of motion in the ball-plate problem is recalled.

Let us consider two right-handed orthonormal frames in the space . One is a fixed frame centered at a point in the lower plane, with perpendicular to this plane and pointing upward, and the other one is a moving frame fixed to the ball at its center. Let be the coordinates of a point on the surface of the ball relative to the fixed frame and let be the coordinates of the same point relative to the moving frame. Thus , where denotes the rotation matrix which transforms the coordinates relative to the moving frame onto the coordinates relative to the stationary frame. Also, if the center of the ball has the coordinates (the radius of the ball is ) and is the coordinates of , both relative to the fixed frame, then . It gives that the motion of the point is described by a curve relative to the fixed frame, or by a path in , where describes the motion of the center of the ball relative to and , and for some path in .

Identifying the tangent space of at any point as the set of all matrices , where is an antisymmetric matrix, it followsfor some functions The vector is the angular velocity generated by . Thus Denoting the velocity of the moving plate and using the nonslipping assumptions, the equations which describe the motion of the ball areConsidering the Lie group and the corresponding Lie algebra generated by the base , where(3) may be written as a controllable systemwhere and are left-invariant vector fields whose values at the identity are, respectively, equal to , .

The control problem is to transfer the ball from a given initial position and orientation to a prescribed final position and orientation along a path which minimizes the functional Using the Maximum Principle, Jurdjevic obtained that the controls which minimize are given bywhere and are the solutions of the systemand .

If , then and it is easy to see that the trajectory of the center of the ball is a circle if , respectively, and a line if . In the sequel we assume that

The paper is organized as follows. At the beginning of the second section we recall Hamilton-Poisson structures of system (8). We give two symplectic realizations of system (8) and their corresponding Lagrangian realizations. In the end of this section we determine some symmetries of the considered system. More precisely, using the Euler-Lagrange’s equations that give a Lagrangian realization, we obtain a one-parameter Lie group of transformations that leaves invariant these equations. The third section is dedicated to the energy-Casimir mapping (). We characterize the equilibrium points as critical points of . Moreover, we give a description of the image of the energy-Casimir mapping and we remark that the boundary of this set is the union of the image through of the stable equilibrium points. Also, we decompose this set into a disjoint union of the connected semialgebraic sets which are related to the image through of the stable and unstable equilibrium points, obtaining semialgebraic stratifications. From the topological point of view, we prove that the fibers corresponding to the strata of these semialgebraic stratifications are, respectively, stable equilibrium points, homoclinic orbits, and periodic orbits. Moreover, we give parametric representations of these fibers. Also, employing energetic methods based on Lyapunov functions, we prove some known stability results and new ones. In the last section of this paper, using the parametric representations of the aforementioned fibers, we mention parametric equations of the trajectory of the center of the ball.

2. Hamiltonian Structures, Lagrangian Realizations, and Symmetries

In this section, we give symplectic and Lagrangian realizations of system (8). Also, some symmetries of the considered system are pointed out.

We recall that the following constants of motion of system (8)were given in [8].

In [9], two Poisson structureswere considered. Moreover, system (8) has the following Hamilton-Poisson realizations and , with the corresponding Casimirs and , respectively. Also, system (8) admits a family of Hamilton-Poisson realizations , where and , , . In this configuration the function , , is a Casimir.

Since , we can conclude that system (8) is a bi-Hamiltonian system.

In the sequel we present two symplectic realizations of system (8).

Theorem 1. The Hamilton-Poisson mechanical system has a full symplectic realization , where and

Proof. Denotingthe mapping , , is a surjective submersion. We notice that . For the Hamiltonian the corresponding Hamilton’s equations areTaking into account relations (13), we havewhence (13) are mapped onto (8). Moreover, we havethat is, the canonical structure induced by is mapped onto the Poisson structure .
Therefore is a full symplectic realization of the Hamilton-Poisson mechanical system .

We remark that is a first integral for system (13).

Theorem 2. The Hamilton-Poisson mechanical system has a full symplectic realization , where and

Proof. The corresponding Hamilton’s equations areAs in the proof of Theorem 1, one obtains that the application ,is a surjective submersion, (17) are mapped onto (8), and the structure induced by is mapped onto the Poisson structure . Also, , which finishes the proof.

Remark 3. For the Poisson structure , the center of the Poisson algebra is generated by the Casimir , . We notice that is a first integral for system (17).

Remark 4. Considering the first integrals and , respectively, and , it is easy to see that system (13) and system (17) respectively are completely integrable in the sense of Arnold-Liouville.

We mention that, in the particular case , , a symplectic realization and its complete integrability were given in [14].

In the next two results, we give Lagrangian realizations of the aforementioned Hamiltonian mechanical systems.

Proposition 5. The Hamiltonian mechanical system considered in Theorem 1 has a Lagrangian realization on the tangent bundle given by the Euler-Lagrange’s equationsgenerated by the Lagrangian

Proof. For the Lagrangian (20) the Euler-Lagrange equations , become (19). Also, differentiating equations (13) one obtains (19). It is easy to see that , where , , which finishes the proof.

Analogously one obtains the following.

Proposition 6. The Hamiltonian mechanical system considered in Theorem 2 has a Lagrangian realization on the tangent bundle given by the Euler-Lagrange’s equationsgenerated by the Lagrangian

In the last part of this section we study some symmetries of system (8). For this purpose, we will use the Lagrangian realization given in Proposition 6.

A symmetry group of a system of differential equations is a group of transformations which maps any solution to another solution of the system. Following [15], let us search for a one-parameter Lie group of transformationswhich leaves (21) invariant. The corresponding infinitesimal generator is given by a vector field :with the property that the action of its second prolongation on (21) vanishes. The second prolongation of the vector field is given byA method to find is given, for example, in [15, 16]. More precisely, taking into account the chain rule for the computation of and replacing by using (21), the relations obtained by applying the second prolongation of on (21) become two equations in , which are all independent. These equations must be satisfied identically in , which leads to the finding of . Detailed utilization of the aforementioned method can be found, for example, in [17, 18].

Proposition 7. The infinitesimal generator of (23) is given bywhere

Proof. It is easy to see that the action of on (21) vanishes, where is given by (26).

Now, some variational symmetries of system (21) can be presented.

Remark 8. (i) For and , we have that represents a translation in the cyclic direction which is related to the conservation of given in Remark 3.
(ii) For and , we have that represents the time translation symmetry which generates the conservation of energy given by (16).
Also, since , , whereand is given by (22), and are variational symmetries.

Remark 9. Taking into account (26), the generators are , , and , which means that a three-parameter group of symmetries was obtained. Denoting , , and , the 3-dimensional Lie algebra of symmetries of (21) endowed with the standard Lie bracket vector fields is generated by the base . Since , , , this Lie algebra is of type V in Bianchi’s classification [19].

3. Properties of the Energy-Casimir Mapping

It is a well-known fact that the dynamics of a Hamilton-Poisson system is foliated by the symplectic leaves associated to the Poisson structure. In [8], Jurdjevic reduced the dynamics (8) to the equation for a pendulum on the level set ,where (see also [11, 12]).

In [12], Mashtakov and Sachkov gave a partition of the cylinder into subsets corresponding to motions of the pendulum (28) of the same type, and using this partition, they introduced elliptic coordinates rectifying the phase flow of the pendulum on a subset of full measure of the cylinder . Furthermore, parametrization of the solutions of (28) is obtained.

The qualitative information about the Hamiltonian mechanical system , , can be given by the so-called energy momentum mapping (see, e.g., [15]). A corresponding of the energy momentum mapping in the case of the Hamilton-Poisson system is the energy-Casimir mapping (see, e.g., [20]). In [20], Tudoran et al. formulated the following open questions. “Is there any connection between the dynamical properties of a given dynamical system and the geometry of the image of the energy-Casimir mapping, and if yes, how can one detect as many dynamical elements as possible (e.g., equilibria, periodic orbits, and homoclinic and heteroclinic connections) and dynamical behavior (e.g., stability, bifurcation phenomena for equilibria, periodic orbits, and homoclinic and heteroclinic connections) by just looking at the image of this mapping?” Particularly, for some dynamical systems, the answers to these questions were given (see, e.g., [2123]).

The goal of this section is to give an answer to the above questions in the case of system (8). Consequently, in this section we study some connections between the dynamics of the considered system and the corresponding energy-Casimir mapping. Our approach regards stability problem, homoclinic and periodic orbits, and their connections with the image of the energy-Casimir mapping and the fibers of the energy-Casimir mapping, respectively.

In the following, we consider system (8) with the Hamilton-Poisson realization , where is given by (10) and the Hamiltonian and a Casimir , given by (9). Consequently, the energy-Casimir mapping corresponding to the considered system is defined by ,Using the critical points of mapping (29), we present a characterization of the equilibrium points of system (8).

Proposition 10. Let . The following assertions are equivalent:(i) is an equilibrium point of system (8);(ii) is a critical point of energy-Casimir mapping (29);(iii).

Proof. Let us observe that system (8) can be written in the equivalent form , . Therefore the equilibria of system (8) are given by the condition .
On the other hand, a point is a critical point of if and only if the derivative of at is not surjective; that is, has rank less than two. Taking into account thatthe conclusion follows.

In order to present the connections between the equilibrium points and the image of the energy-Casimir mapping, we give a description of the image of the energy-Casimir mapping. Also, we present a study of the stability of equilibrium points.

By the image of the energy-Casimir mapping is understood by the setNow, we describe the image of the energy-Casimir mapping.

Proposition 11. The image of the energy-Casimir mapping (29) is given by

Proof. We are looking for such that the systemhas solutions. It is clear and . Denoting , , we obtainIf we choose such that and , it giveswhich implies the conclusion.

The graphical description of is given in Figure 1.

Remark 12. The curve of the equation , , is an arch of the parabola of equation . One gets the image of the energy-Casimir mapping (32) is the disjoint union of the following connected semialgebraic sets (Figure 2):Following, for example, [24], we obtain the following semialgebraic stratifications: and . Here, and are principal strata and, as we will see in the following, the superscripts and stand, respectively, for stable and unstable.

By Proposition 7, the union gives the equilibrium points of system (8), whereIn [10] the stability problem was treated by using the energy-Casimir method. However, the stability of the equilibrium points and was not established. In the following we complete the study of the stability. Also, using Lyapunov functions, we give another proof of the known stability results.

Proposition 13 (see, also, [10]). All the equilibrium points from the family are nonlinearly stable.

Proof. Considering the function ,we havefor some neighborhood of . Moreover, using (8), it followsfor all
By [25, 26], we deduce that all the equilibrium states from the family are nonlinearly stable.

Proposition 14 (see, also, [10]). The equilibrium points from are nonlinearly stable for and unstable otherwise.

Proof. The characteristic polynomial associated with the linear part of system (8) at the equilibrium is given byIf we notice that a root of is strictly positive, whence is an unstable equilibrium state. Therefore, all the equilibrium states from the family are unstable in the case .
If , then , whence it is nonlinearly stable.
Now, if , let us consider the function ,It follows and, using (8), for all . Therefore, in order to prove the nonlinear stability of the equilibrium point , it only remains to prove that the condition implies ; that is, the systemhas a unique solution, .
If , it is obvious that the above system has a unique solution, namely, . Thus this equilibrium state is nonlinearly stable.
Let now . From geometric point of view, the solution of system (43) is the curve of intersection between a circular cylinder and a sphere. This curve is reducing at a point if and only if the sphere and the cylinder are tangent. In this geometrical configuration, it is enough to prove that the circles obtained by the intersection of the aforementioned surfaces with the plane are tangent that is the following system has a unique solution:By straightforward computation, one obtains , , which finishes the proof.

The diagram of equilibrium points is given in Figure 3.

In the next result we give some connections between the image through the energy-Casimir mapping of the families of equilibrium states and the set .

Proposition 15. Let and and be the families of nonlinearly stable equilibrium points and unstable equilibrium points, respectively, and let be the energy-Casimir mapping (29) of system (8). Then,(i), where denotes the boundary of the set ;(ii), , , and ;(iii), where denotes the convex hull.

Proof. The conclusions result by simple computations.

Remark 16. We notice that if , then and belongs to a line. Also, if , then and belongs to a curve, namely, a parabola.

In the sequel we study the topology of the fibers of the energy-Casimir mapping (29). A fiber of the energy-Casimir mapping is the preimage of an element through , that is, By fixing an element from , it belongs to one stratum (36). Therefore the classification of the fibers will be done taking into account the stratification of . Since the dynamics has place to the intersection between the surfaces and , one obtains that these fibers are, respectively, a stable equilibrium point, a periodic orbit, a pair of homoclinic orbits, a pair of periodic orbits, and two stable equilibrium points (Figures 4, 5, 6, and 7).

Now we prove the observed results. The first proposition refers to the fibers related to the stable equilibrium states.

Proposition 17. Let . Then the corresponding fiber contains only stable equilibrium points; namely,(i)if , then ;(ii)if , then ;(iii)if , then .

Proof. (i) Using (36) and (45), the conclusion immediately follows.
(ii) In this case the fiber is the intersection between the cylinder and the sphere . Observing that the sphere and the cylinder are tangent (interior or exterior), it follows that contains only one point, . Now, by straightforward computation, one obtains the conclusion.

Just looking to the image of the energy-Casimir mapping (Figure 2) and taking into account the papers [20, 22, 23], we expect that there exist homoclinic orbits in dynamics (8). Moreover, Figure 6 suggests that the fibers corresponding to the stratum are homoclinic orbits. More precisely, we have the following.

Proposition 18. If , then consists of a pair of homoclinic orbits, namely, , wherewith , ,

Proof. Following [8], we denoteIt follows . Since , it gives . Now, taking into account that , one getsIf we consider such thatand we denote , then (48) becomesand using , we obtainwhere . It follows , , with , whence . Computing and , we obtain (46).
Now, let us observe that by fixing we have , where , . It is easy to see that , whence the curves given parametrically by (46) are homoclinic orbits.

In the papers [2023, 27] it was reported that the fibers corresponding to the principal strata of the image of the energy-Casimir mapping associated to some particular systems of differential equations are periodic orbits. In our case, Figures 5 and 7 suggest that the fibers corresponding to the principal strata are also periodic orbits. In the following two results we prove the aforementioned connections.

Proposition 19. The fiber corresponding to a point belonging to the principal stratum , wherewith , , , and the modulus of the Jacobi’s elliptic functions , is a periodic orbit.

Proof. Since , it followsThus there exists such thatFollowing the steps (47)–(50) from the proof of Proposition 18, (50) becomesTherefore the equation of a mathematical pendulum is obtained. Using Jacobi’s elliptic functions (see, e.g., [28]) it follows , where and the modulus . Then . Using the properties of the Jacobi’s elliptic functions one gets . Also, by (47), one obtains and , which finishes the proof.

Proposition 20. Let . Then represents a pair of periodic orbits, wherewith , , , and the modulus of the Jacobi’s elliptic functions .

Proof. Performing the steps (47)–(50), (50) becomeswhere . Because it gives . We again obtain the equation of a mathematical pendulum and following [28], it follows , where . By straightforward computation the parametrization (56) is obtained.

4. The Trajectory of the Center of the Ball

Using a connection with the elastica problem, Jurdjevic [8] gave a completely geometric analysis of the motion of the center of the ball; he described various qualitative types of rolling for a sphere along elastics of various forms: inflectional, noninflectional, a circle, and a straight line.

In [12], using the coordinates obtained by integration of the pendulum motion on the partition of a cylinder, a parametrization of extremal trajectories is obtained.

In this section, we classify the types of the trajectory of the center of the ball taking into account the fibers of the energy-Casimir mapping. Using the parametric representations of the studied fibers, we obtain the parametric representations of the trajectory of the center of the ball. These representations are similarly with the ones discussed in [8], respectively [12].

As we have seen in the Introduction, namely, (3) and (7), the controlled motion of the center of the ball is described by the following system:where and are given by system (8) and .

In the following, considering that belongs to a specific fiber , we will integrate system (58):

Proposition 21. Let and . Then the ball does not move.

Proof. By Proposition 17, it is obvious.

Proposition 22. Let and . Then the ball rolls along the line

Proof. By Proposition 17 (ii), it immediately follows.

Proposition 23. If , where , then the trajectory of the center of the ball is given by or , wherewith

Proof. Using (46) with , it is easy to integrate system (58).

Figure 8 presents the situation given in Proposition 23. The blue curve represents the trajectory , , and the red one the trajectory , where .

In the following results we are using the Jacobi’s epsilon function [28]:where is the modulus.

Proposition 24. Let and . Then the orbit of the center of the ball is , wherewith and the modulus of the Jacobi’s elliptic functions .

Proof. Replacing (52) with in (59) and using (62) and the property , the conclusion follows.

In Figure 9 we present some trajectories of the center of the ball when is fixed and increases such that , first near the stratum and last near the stratum .

Proposition 25. If , where , then the orbit of the center of the ball is given by (Figure 10, blue) or (Figure 10, red), wherewith and the modulus of the Jacobi’s elliptic functions .

Proof. Replacing (56) with in (59) and taking into account , making use of (62) and the property , we obtain (64).

In Figure 10 some trajectories of the center of the ball are presented when is fixed and increases such that , first near the stratum .

Remark 26. We notice that the relations (63) and (64) can be written by using elliptic integrals. Indeed, we have [28], where is the Jacobi amplitude and is the elliptic integral of the second kind with the modulus .
Using Wolfram Mathematica software, one can draw a variety of possible trajectories of the center of the ball by customizing and and taking belonging to different strata.

5. Conclusions

In our work some properties of the dynamic of the ball-plate problem are studied. We particularly refer here to the connections between the dynamical properties of the considered system of differential equations and the geometrical properties of the corresponding energy-Casimir mapping . Thus, the equilibrium points are the critical points of , and the boundary of the image of the energy-Casimir mapping is the image through of the stable equilibrium states. Also, using the image through of the equilibrium states, semialgebraic partitions of the image of the energy-Casimir mapping are obtained. In addition, the fibers corresponding to the strata of the presented semialgebraic stratifications are, respectively, stable equilibrium points, homoclinic orbits, and periodic orbits. More precisely, the fibers corresponding to the unstable stratum are a pair of homoclinic orbits and the fibers corresponding to the principal strata are periodic orbits. It is an open problem to establish a class of three-dimensional system of differential equations with the same properties. In the end of the paper, using the parametric representations of the aforementioned fibers, parametric equations of the trajectory of the center of the ball are mentioned.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.