Partial Integrability Conditions and an Integrating Algorithm for Fully Rheonomous Affine Constraints

Abstract

In this paper, integrability conditions and an integrating algorithm of fully rheonomous affine constraints (FRACs) for the partially integrable case are studied. First, some preliminaries on the FRACs are illustrated. Next, necessary and sufficient conditions on the partially integrable case for the FRACs are derived. Then, an integrating algorithm to calculate independent first integrals of the FRACs for the partially integrable case is derived. Moreover, the existence of an inverse function utilized in the algorithm is proven. After that, an example is presented for evaluation of the effectiveness of the proposed method. As a result, it turns out that the proposed integrating algorithm can easily calculate independent first integrals for given partially integrable FRACs, and thus this new algorithm is expected to be applied to various research fields.

Share and Cite:

Kai, T. (2014) Partial Integrability Conditions and an Integrating Algorithm for Fully Rheonomous Affine Constraints. Circuits and Systems, 5, 133-141. doi: 10.4236/cs.2014.56015.

1. Introduction

A lot of researches on nonholonomic systems have been actively done in the research fields of nonlinear control theory and robotics so far [1] -[5] . Moreover, sub-Riemannian geometry has been also studied in the research fields of differential geometry and control theory [6] [7] . In these two research fields, the common key is the existence of constraints, which is represented as differential equations. The concepts “integrability” and “nonintegrability” of such constraints are quite important since they determine various characteristics of systems. Hence, investigation of integrability and nonintegrability is an essential task in the first stage.

There are some classes in constraints, and the simplest class of constraints is called linear constraints:, which has been mainly studied so far. The class of the linear constraints covers wide-ranging mechanical systems such as mobile caars, acrobatic robots, and biomimetic robots. As a wider class of constraints than the linear constraints, the author has focused and researched scleronomous affine constraints:, and A-rheonomous affine constraints: from the viewpoints of mathematics and control theory [8] -[12] . It must be noted that in the field of analytical mechanics, the terminology “rheonomous” means “time-varying,” and the opposite word of it is “scleronomous.” The affine constraints can be found in physical and mechanical systems: space robots with initial angular momenta, a ball on a rotating table, a ship on a running river, and so on. These work have made it possible to treat such constraints systematically, however, we are still interested in fully rheonomous affine constraints: as a much wider class of constraints than the scleronomous affine constraints. In [13] , the author has derived a complete integrability condition for the fully rheonomous affine constraints. Moreover, in [14] , an integrating algorithm to calculate first integrals of the completely integrable fully rheonomous affine constraints is developed. The word “complete integrability” means all the constraints are integrable. However, there is another case called “partial integrability”, which means the situation where integrable constraints and nonintegrable ones are mixed. Researches on partially integrable fully rheonomous affine constraints have not been done.

Hence, the aim of this paper is to develop partial integrability conditions and an integrating algorithm for the fully rheonomous affine constraints. This paper is organized as follows. First, in Section 2, some preliminaries on the fully rheonomous affine constraints are shown. Next, in Section 3, necessary and sufficient conditions on partial integrability for the fully rheonomous affine constraints are derived. Then, an integrating algorithm to calculate independent first integrals for partially integrable rheonomous affine constraints is developed in Section 4. After that, Section 5 illustrates in order to verify the effectiveness and the availability of the new results.

2. Preliminaries

2.1. Fully Rheonomous Affine Constraints (FRACs)

This section presents some preliminaries on fully rheonomous affine constraints are presented. See [13] for more details. First, this subsection gives the definition of fully rheonomous affine constraints and explains their geometric representation. Denote the time variable by and a time interval by. Let be an n-dimensional configuration manifold and be a local coordinate of. Associated with, we refer as a tangent vector field. A set of of differential equations in the form:

(1)

is called fully rheonomous affine constraints (FRACs). Note that all the coefficients explicitly depend on the time variable. We now rewrite (1) as

(2)

where a rheonomous affine term is a vector-valued function whose i-th entry is, and a rheonomous velocity coefficient matrix is a matrix-valued function whose ij-th entry is. In this paper, we assume the following sufficient condition on independency of the FRACs (2):

(3)

Next, a geometric representation method of the FRACs (2) is explained. From (3), we see that the row vectors of in (2) are independent of each other. Hence, we consider vector fields which are independent of each other and annihilators of the row vectors of, and denote them by as time-varying vector fields on. Furthermore, we also denote a space spanned by, that is, a time-varying distribution on by

(4)

Since the basial vectors of: are independent of each other, is a nonsingular distribution, that is,

(5)

holds. A curve on: is said to satisfy the FRACs (2) if for a time-varying vector field on: and the generalized velocity of:,

(6)

holds. We call a rheonomous affine vector field and it satisfies the equation:

(7)

This definition is a natural extension of the one for the scleronomous affine constraints that do not contain the time variable explicitly [8] . Geometric representation of the FRACs is defined as follows and can allow us to analyze them geometrically and derive geometric properties.

 Definition 1: The FRACs (2) are geometrically represented by a pair, where is an m-dimensional time-varying distribution defined by (4) and is called a rheonomous affine vector and satisfies (7).                                                                                  □

2.2. Rheonomous Bracket

Next, this subsection explains a new operator for the fully rheonomous affine constraints (2), called the rheonomous bracket. The rheonomous bracket is originally introduced in order to analyze the A-rheonomous affine constraints in [10] -[12] , and plays important roles in derivation of a complete integrability condition and an integrating algorithm. The rheonomous bracket is fundamentally defined based on the normal Lie bracket which is an operator for two vector fields:

(8)

The definition of the rheonomous bracket for the FRACs (2) is as follows [10] -[12] .

 Definition 2 [10] -[12] : For the vector fields defined on on the geometric representation of the FRACs (2):, the rheonomous bracket is an operator: that satisfies the following three properties:(a) For a rheonomous affine vector field, (9)holds.(b) is defined as a set of vector fields that consists of and iterated rheonomous brackets of and does not contain. For a rheonomous affine vector field and a vector field, (10)holds.(c) For two vector fields, (11)holds.                                                                                   □

For the rheonomous bracket, it is noted that the rheonomous affine vector field is perceived as special, and this yields an additional term of a time differential of a vector field as the property (b). It must be also noted that from Definition 2 the rheonomous bracket is equivalent to the normal Lie bracket for scleronomous affine constraints, that is, constraints that do not contain the time variable explicitly. The following proposition shows that the rheonomous bracket has some important characteristics in common with the normal Lie bracket [10] -[12] .

 Proposition 1: For the vector fields on the geometric representation of the FRACs (2): and the set of iterated vector fields of them:, the following properties (a), (b), and (c) hold.(a) Bilinearlity: (12)(b) Skew-symmetry: (13)(c) Jacobi’s identity: (14)□

From Proposition 1, it can be confirmed that we only have to consider the iterated rheonomous brackets in the form:

(15)

in checking a complete integrability condition for the FRACs, which will be shown in the next subsection. Furthermore, the Philip Hall basis [2] , which is a systematic method to generate iterated Lie brackets with an order efficiently, can be also constructed for the rheonomous bracket as follows [10] -[12] .

 Algorithm 1: For iterated rheonomous brackets (15) of the geometric representation of the FRACs (2):, we define the length of (15) as, that is, the number of vector fields in the iterated rheonomous bracket. In addition, the symbol means the magnitude relation for two iterated rheonomous brackets. Then, the Philip Hall basis for the rheonomous bracket can be constructed by the next rules.(a) are the first elements of and.(b) If, then.(c) if and only if and, either or holds or with and.                                             □

3. Partial Integtability Conditions for FRACs

In this section, a partial integrability condition for the FRACs (2) is derived. In the simplest terms, the partially integrable case means that some constraints are integrable and the others are nonintegrable. If the FRACs in (2) are integrable, that is, there exist independent first integrals of (2), then they are said to be k-th order partially integrable. Now, we define a smallest and involutive time-varying distribution which contains and iterated rheonomous brackets of them, and satisfies, that is, is spanned by all the rheonomous brackets of with the exception of. Then, necessary and sufficient conditions on partilly integrability for the FRACs (2) are given as the next theorem.

 Theorem 1: For the fully rheonomous affine constraints defined on an n-dimensional manifold (2) and a time interval, the following statements (a) and (b) are equivalent to each other. If they hold, the fully rheonomous affine constraints (2) are said to be k-th order partially integrable.(a) There exist independent first integrals of the fully rheonomous affine constraints (2).(b) For a smallest and involutive time-varying distribution, (16)holds.

(Proof) Let us consider the -dimensional product space, where is the space of the time variable, and the local coordinate of is given by. On, the FRACs (2) can be represented by Pfaffian equations of differential forms:

(17)

Since the rheonomous affine vector field of geometric representation satisfies (7), vector fields on which annihilate (17) are given by

(18)

We now define an involutive distribution defined on which contains and iterated Lie brackets for. Therefore, by Frobenius’ theorem [1] [2] [15] [16] and the results [13] [14] , It can be confirmed that a necessary and sufficient condition on k-th order partial integrability for (17) is given by

(19)

Calculating the iterated Lie brackets for, then we obtain

(20)

We can easily see that is independent of and the iterated Lie brackets (20). Hence, we can change the necessary and sufficient condition (19) into the condition such that and the iterated Lie brackets which consist of span an k-dimensional space. From the calculations (18) and (20), we only have to consider on instead of on, and iterated rheonomous brackets for on instead of Lie brackets for on. Consequently, a necessary and sufficient condition on k-th order partially integrability for the FRACs (2) is that and the iterated rheonomous brackets for span an k-dimensional space, that is, (16) holds.

From Theorem 1, it turns out that the partial integrability condition for the FRACs (2) is quite simple and is a natural extension of the case of the complete integrability case [13] [14] . Moreover, we can see that the rheonomous bracket plays a significant role in the condition (16).

4. Integrating Algorithm for Partially Integrable FRACs

4.1. Integrating Algorithm

The result in Section 3 shows that if the FRACs (2) are k-th order partially integrable, there exist independent first integrals of them. For reduction of the dimension of a given configuration manifold subject to partially integrable FRACs, we need the explicit forms of independent first integrals of them. The author has developed algorithms to calculate independent first integrals of them in the case of complete integrability [14] . However, the algorithm cannot be utilized for the case of partial integrability. Hence, in this subsection, an integrating algorithm for partially integrable FRACs (2) is constructed.

Assume that the FRACs (2) satisfy Theorem 1, that is, they are k-th order partially integrable. Then, we can find vector fields that satisfy

(21)

In addition, set linearly independent vector fields for, that is

(22)

holds. Then, let us denote flows (1-parameter local transformation groups) of and by and with parameters and, respectively. We set an initial point at the initial time as. We also consider vector fields defined on the expanded configuration manifold (18), and then their flows on are obtained as

(23)

Calculating the composite mapping of flows (23) yields

(24)

where

(25)

is the composite mapping of flows. From (24), we can see that the projection of onto is equivalent to. Hence, by applying the idea of the integrating algorithm for scleronomous linear constraints defined on [1] [2] to and considering projection of it onto, Consequently, we have the following algorithm to calculate independent first integrals of partially integrable FRACs.

 Algorithm 2: For the case where the A-rheonomous affine constraints (2) are -th order partially integrable, we can obtain independent first integrals of them by the following procedure.Step 1: Set vector fields of geometric representation for the FRACs (2).Step 2: By calculating iterated rheonomous brackets of, derive linearly independent vector fields that satisfy (26)Step 3: For, derive linearly independent vector fields that satisfy (27)Step 4: Calculate flows (1-parameter local transformation groups) of with parameters, respectively:.Step 5: Combine flows calculated in Step 4 as (28)Step 6: Set and calculate the inverse function. Then, the last components of are first integrals of k-th order partially integrable FRACs (2).                                  □

4.2. Theoretical Analysis of Integrating Algorithm

Next, this subsection gives theoretical analysis on the integrating algorithm for partially integrable FRACs (Algorithm 2). In Algorithm 2, we have to derive the inverse function of the combined mapping: in order to calculate independent first integrals. However, the existence of the inverse function is not guaranteed. In general, it is quite difficult to calculate an inverse function of a given function. For Algorithm 2, the next proposition guarantees the existence of the inverse mapping.

 Proposition 2 Assume that the FRACs (2) are k-th order partially integrable. Then, there exists an time interval and is a diffeomorphism at any time. That is to say, there exists its inverse mapping.

(Proof) Set. Calculating the partial differential of (24) with the chain rule of differential calculation, we have

(29)

Substituting into (29), we obtain

(30)

that is to say,

(31)

Since are linearly independent of each other,

(32)

holds. Therefore, it turns out that is a diffeomorphism by the implicit function theorem [15] [16] . Since the projection of onto is equivalent to, is also a diffeomorphism. Consequently, the proposition is proven.

From Proposition 2, we can see that similar to the case of complete integrability [14] , the existence of the inverse mapping is also guaranteed. For an example of calculation of independent first integrals for the partially integrable FRACs, see the next section.

5. Example

Finally, this section deals with an example to evaluate the new results. Let us consider a 3-dimensional configuration manifold:

(33)

with, and a fully rheonomous affine constraints on:

(34)

with. We here consider a time interval. Then, it turns out that Assumption 1 holds for (34). One geometric representation for (34) can be obtained as follows:

(35)

Calculating an iterated rheonomous bracket for and above, we obtain

(36)

Since holds, we have

(37)

and then it can be confirmed that

(38)

holds. From Theorem 1, it turns out that the FRACs (34) are 2nd order partially integrable, that is, there exists one independent first integrals of (34).

Next, we shall calculate the first integrals of (34) according to Algorithm 1. Set and a new vector field that satisfy (26) as

(39)

For the vector fields, we calculate their flows as

(40)

where is the initial point at the initial time. Combining the flows (40) like (28), we have

(41)

By solving the equation, we calculate the inverse mapping of (41) as

(42)

Consequently, we can obtain one independent first integrals of (34) as the last component of (43):

(43)

6. Conclusions

This paper has studied necessary and sufficient conditions for the fully rheonomous affine constraints (FRACs) in the case of partial integrability, and develops an integrating algorithm in order to calculate independent first integrals. It turns out that the new proposed integrating algorithm is expected to be a powerful mathematical tool and can be utilized in various research fields.

Future work includes applications of the new results to various systems, extensions of sub-Riemannian geometry, and analysis of more general classes of constraints.

Conflicts of Interest

The authors declare no conflicts of interest.

 [1] Isidori, A. (1995) Nonlinear Control Systems. 3rd Edition, Springer-Verlag, London. http://dx.doi.org/10.1007/978-1-84628-615-5 [2] Sastry, S.S. (1999) Nonlinear Systems. Springer-Verlag, New York. http://dx.doi.org/10.1007/978-1-4757-3108-8 [3] Cortés, J. (2002) Geometric, Control and Numerical Aspects of Nonholonomic Systems. Springer-Verlag, Berlin Heidelberg. http://dx.doi.org/10.1007/b84020 [4] Bloch, A.M. (2003) Nonholonomic Mechanics and Control. Springer-Verlag, New York. http://dx.doi.org/10.1007/b97376 [5] Bullo, F. and Rewis, A.D. (2004) Geometric Control of Mechanical Systems. Springer Science+business Media, Inc. [6] Montgomery, R. (2002) A Tour of Subriemannian Geometries, Their Geodesics and Applications. American Mathematical Society. [7] Calin, O. and Change, D.C. (2009) Sub-Riemannian Geometry: General Theory and Examples. Cambridge University Press, Cambridge. http://dx.doi.org/10.1017/CBO9781139195966 [8] Kai, T. and Kimura, H. (2006) Theoretical Analysis of Affine Constraints on a Configuration Manifold—Part I: Integrability and Nonintegrability Conditions for Affine Constraints and Foliation Structures of a Configuration Manifold. Transactions of the Society of Instrument and Control Engineers, 42, 212-221. [9] Kai, T. (2011) Integrating Algorithms for Integrable Affine Constraints. IEICE Transactions on Fundamentals of Electronics, Communications and Computer Sciences, E94-A, 464-467. [10] Kai, T. (2012) Mathe-matical Modelling and Theoretical Analysis of Nonholonomic Kinematic Systems with a Class of Rheonomous Affine Constraints. Applied Mathematical Modelling, 36, 3189-3200. http://dx.doi.org/10.1016/j.apm.2011.10.015 [11] Kai, T. (2012) Theoretical Analysis for a Class of Rheonomous Affine Constraints on Configuration Manifolds—Part I: Fundamental Properties and Integrability/Nonintegrability Conditions. Mathemati-cal Problems in Engineering, 2012, Article ID: 543098. [12] Kai, T. (2012) Theoretical Analysis for a Class of Rheonomous Affine Constraints on Configuration Manifolds—Part II: Foliation Structures and Integrating Algorithms. Mathematical Problems in Engineering, 2012, Article ID: 345942. [13] Kai, T. (2013) On Integrability of Fully Rheonomous Affine Constraints. International Journal of Modern Nonlinear Theory and Application, 2, 130-134. http://dx.doi.org/10.4236/ijmnta.2013.22016 [14] Kai, T. (2013) An Integrating Algorithm and Theoretical Analysis for Fully Rheonomous Affine Constraints: Completely Integrable Case. Applied Mathematics, 4, 1720-1725. http://dx.doi.org/10.4236/am.2013.412235 [15] Nomizu, S. and Kobayashi, K. (1996) Foundations of Differential Geometry Volume I. Wiley-Inter-science. [16] Nomizu, S. and Kobayashi, K. (1996) Foundations of Differential Geometry Volume II. Wiley-Inter-science.