Open Access
Issue
Wuhan Univ. J. Nat. Sci.
Volume 29, Number 3, June 2024
Page(s) 273 - 283
DOI https://doi.org/10.1051/wujns/2024293273
Published online 03 July 2024

© Wuhan University 2024

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

The Gauss principle is a differential variational principle proposed by Gauss in 1829, a general analytical mechanics principle[1]. Chen pointed out[2] that taking the Gauss principle as the fundamental principle in terms of mechanics concepts seems most appropriate. According to Mei[3], the Gauss principle can be used as a basis for analytical dynamics. Udwadia and Kalaba[4] took the Gauss principle as a starting point to derive the basic equations of analytical mechanics by using matrix algebra operations and recommended its application to holonomic and non-holonomic mechanics, which reveals the broad applicability of Gauss principle in describing the motion of constrained mechanical systems. Of all differential variational principles, only the Gauss principle is a stationary principle, which shows that the variation of compulsion function in the sense of Gauss is equal to zero[1]. For a system of particles, applying Gauss minimum compulsion principle, its motion equation can be obtained directly by calculating the extremum of the compulsion function[5,6]. Because of this, the Gauss principle of least compulsion is widely used in dynamics modeling and in finding approximate solutions. For example, robot dynamics[5], multi-body system dynamics[6-14], elastic rod dynamics[15-18], and hybrid dynamics[19], etc. So far, there have been many achievements in the Gauss and least compulsion principles for constrained mechanical systems and their applications[20-28].

Using the analytical mechanics method to study the relative motion dynamics of complex systems can unify the expression forms and show the superiority of analytical mechanics in solving the dynamics problems of complex systems. These complex systems comprise the carrier body and carried bodies moving relative to the former[29]. There are many such systems in practical engineering, such as the relative motion and control of spacecraft, the relative motion of satellites[30-33] and so on. Whittaker analyzed a holonomic system subject to uniform rotation constraints and derived its Lagrange equations[34]. Lurie studied holonomic mechanics with relative motion[35]. Mei et al extended it to nonholonomic mechanics[36]. Since then, progress has been made in the variational principle, equations of motion, integral theory, and symmetry of relative motion dynamics[37-44]. The Gauss principle of relative motion dynamics is further studied in this paper. Section 1 introduces the establishment of the Gauss principle for relative motion dynamics by analyzing the virtual displacement of acceleration space. In Section 2, the compulsion function of relative motion dynamics is constructed, and it is proved that real motion causes the compulsion function to reach an extreme value under Gaussian variation. Section 3 gives the formulae of acceleration energy and corresponding compulsion function when the carried body is a rigid body whose relative motion is planar motion. In Section 4, we study Gauss principle of relative motion dynamics and give its Appell, Lagrange, and Nielsen forms in generalized coordinates. In Section 5, from Gauss principle we obtained, we deduce dynamical equations with relative motion. In Section 6, two examples are given. Section 7 is the conclusion of the article.

1 Gauss Principle of Relative Motion Dynamics

Study a system of particles that comprises a rigid body (carrier) and N particles (carried bodies). The carried bodies are moving relative to the carrier. The moving frame of reference Oxyz is attached to the carrier. We use n generalized coordinates qs to describe the configuration of relative motion s=1,2,,n. The acceleration aO of the point O in a fixed frame O1x1y1z1, and the angular velocity ω of moving frame are the given functions of time t. For the i-th particle, let mi be its mass and ri'=ri'(qs,t) its position vector relative to Oxyz. The dynamic equation of relative motion is

- m i r ¨ ˜ i ' + F i   + N i + F e i I + F c i I = 0 ,   i = 1,2 , , N (1)

where Fi, Ni, FeiI=-miaO-miω˙×ri'-miω×(ω×ri'), FciI=-2miω×r˙˜i' are the active force, the constraint force, the convective inertial force, and the Coriolis inertial force, respectively. r˙˜i' is the relative velocity, r¨˜i' is the relative acceleration.

By dotting the equation (1) with δGr¨˜i' and summing over i, we get

i = 1 N ( - m i r ¨ ˜ i ' + F i   + N i + F e i I + F c i I ) δ G r ¨ ˜ i ' = 0 (2)

where δG(·) stands for the Gaussian variation[3]. Within acceleration space, the condition of ideal constraints yields

i = 1 N N i δ G r ¨ ˜ i ' = 0 (3)

Thus, formula (2) becomes

i = 1 N ( - m i r ¨ ˜ i ' + F i   + F e i I + F c i I ) δ G r ¨ ˜ i ' = 0 (4)

Formula (4) is the Gauss principle of relative motion dynamics.

2 Gauss Principle of Least Compulsion for Relative Motion Dynamics

The compulsion function of relative motion is explained as

Z r = i = 1 N 1 2 m i ( r ¨ ˜ i '   - F i   + F e i I + F c i I m i ) 2 (5)

then we have

δ G Z r = i = 1 N m i ( r ¨ ˜ i ' - F i   + F e i I + F c i I m i ) δ G ( r ¨ ˜ i ' - F i   + F e i I + F c i I m i ) = i = 1 N m i ( r ¨ ˜ i ' - F i   + F e i I + F c i I m i ) δ G r ¨ ˜ i '   (6)

Therefore, the principle (4) becomes

δ G Z r = 0 (7)

If r¨˜i' is the relative acceleration in real motion and r¨˜i'+δGr¨˜i' is of possible motion of which the constraints admit, subsequently, the distinction between the compulsion functions is

Δ Z r = i = 1 N 1 2 m i { ( r ¨ ˜ i ' - F i + F e i I + F c i I m i ) 2 - ( r ¨ ˜ i ' + δ G r ¨ ˜ i ' - F i + F e i I + F c i I m i ) 2 } = - i = 1 N 1 2 m i ( δ G r ¨ ˜ i ' ) 2 - i = 1 N m i ( r ¨ ˜ i ' - F i + F e i I + F c i I m i ) δ G r ¨ ˜ i ' = - i = 1 N 1 2 m i ( δ G r ¨ ˜ i ' ) 2 < 0 (8)

Thus, equation (7) shows that, at any instant, the real motion of a relative motion dynamic system minimizes the compulsion function Zr under Gaussian variation when compared with possible motions with the same configuration and the same velocity but with different accelerations. Equation (7) can be called the Gauss principle of least compulsion for relative motion dynamics. When, aO=0,ω=0, principles (4) and (7) degenerate to the classical Gauss principle and the least compulsion principle on the absolute motion[3].

3 Calculation of Acceleration Energy and Compulsion Function

Expanding formula (5), we have

Z r = i = 1 N 1 2 m i r ¨ ˜ i ' r ¨ ˜ i ' - i = 1 N ( F i + F e i I + F c i I ) r ¨ ˜ i ' + (9)

where the ellipsis "" symbolizes the terms that are independent of relative acceleration.

Let Sr denote the acceleration energy of relative motion, i.e.,

S r = i = 1 N 1 2 m i r ¨ ˜ i ' r ¨ ˜ i ' (10)

The compulsion function gives

Z r = S r - i = 1 N ( F i   + F e i I + F c i I ) r ¨ ˜ i ' + (11)

Next, we study the calculation of acceleration energy of relative motion if the carried body is rigid. First, if the relative motion is translation, denote the center of mass of the carried body as C and its relative acceleration aCr, then

S r = i = 1 N 1 2 m i r ¨ ˜ i ' r ¨ ˜ i ' = i = 1 N 1 2 m i a i r 2 = 1 2 m a C r 2 (12)

where air is the relative acceleration of the i-th particle, and m=i=1Nmi is the total mass. Second, in the case of fixed-axis rotation for relative motion, denote the relative angular velocity of the carried body ωr, the relative angular acceleration εr, and the moment of inertia about the rotation axis Aξ as Jξ, then

S r = i = 1 N 1 2 m i r ¨ ˜ i ' r ¨ ˜ i ' = i = 1 N 1 2 m i a i r a i r = i = 1 N 1 2 m i ( ρ i ε r τ i + ρ i ω r 2 n i ) ( ρ i ε r τ i + ρ i ω r 2 n i ) (13)

where ρi is the distance between the i-th particle and Aξ axis, and unit vectors τi and ni are along tangential and principal normal directions, respectively. Expanding equation (13), we get

S r = i = 1 N 1 2 m i ρ i 2 ε r 2 + i = 1 N 1 2 m i ρ i 2 ω r 4 = 1 2 J ξ ε r 2 + (14)

Third, in the event that the relative motion is planar motion, denote the relative angular velocity of the carried body with planar motion as ωr, the relative angular acceleration as εr, the relative acceleration as aCr, then

S r = i = 1 N 1 2 m i r ¨ ˜ i ' r ¨ ˜ i ' = i = 1 N 1 2 m i a i r a i r = i = 1 N 1 2 m i [ a C r + ρ i ( ε r τ i + ω r 2 n i ) ] [ a C r + ρ i ( ε r τ i + ω r 2 n i ) ] = i = 1 N 1 2 m i [ a C r 2 + 2 a C r ρ i ( ε r τ i + ω r 2 n i ) + ρ i 2 ( ε r 2 + ω r 4 ) ] (15)

where ρi is the distance between the i-th particle and C, and τi and ni are tangential and normal unit vectors relative to C, respectively. Obviously, from the centroid coordinate formula, we have

i = 1 N m i ρ i a C r τ i = 0 , i = 1 N m i ρ i a C r n i = 0 (16)

Hence, we obtain

S r = 1 2 m a C r 2 + 1 2 J C ε r 2 (17)

where JC=i=1Nmiρi2 represents the moment of inertia. Equation (17) shows that the acceleration energy of the relative motion of the carried body with planar motion equals the sum of the acceleration energy of relative translation with and relative rotation around the center of mass. Let

F = i = 1 N F i , F e I = i = 1 N F e i I , F c I = i = 1 N F c i I (18)

represent the principal vectors of the active forces, the convective inertial forces, and the Coriolis inertial forces, respectively, and

M C = i = 1 N M C ( F i ) , M C e I = i = 1 N M C ( F e i I ) , M C c I = i = 1 N M C ( F c i I ) (19)

represent the principal moment about point C, then

i = 1 N ( F i   + F e i I + F c i I ) r ¨ ˜ i ' = i = 1 N ( F i   + F e i I + F c i I ) [ a C r + ρ i ( ε r τ i + ω r 2 n i ) ]

= a C r i = 1 N ( F i   + F e i I + F c i I ) + ε r i = 1 N ρ i τ i ( F i   + F e i I + F c i I ) + ω r 2 i = 1 N ρ i n i ( F i   + F e i I + F c i I ) = a C r ( F   + F e I + F c I ) + ε r ( M C   + M C e I + M C c I ) + (20)

Substituting formulas (17) and (20) into (11), we get

Z r = 1 2 m a C r 2 + 1 2 J C ε r 2 - a C r ( F   + F e I + F c I ) - ε r ( M C   + M C e I + M C c I ) + (21)

This formula calculates the compulsion function of relative motion for the carried body in planar motion.

If the relative motion is translation, then the compulsion function (21) gives

Z r = 1 2 m a C r 2 - a C r ( F   + F e I + F c I ) + (22)

If the relative motion is fixed axis rotation, then the compulsion function (21) provides

Z r = 1 2 J ξ ε r 2 - ε r ( M ξ   + M ξ e I + M ξ c I ) + (23)

4 Gauss Principle of Relative Motion Dynamics in Generalized Coordinates

Taking the relative derivative of ri'=ri'(qs,t), we get

r ˙ ˜ i ' = s = 1 n r i ' q s q ˙ s + r i ' t (24)

r ¨ ˜ i ' = s = 1 n r i ' q s q ¨ s + s = 1 n k = 1 n 2 r i ' q s q k q ˙ s q ˙ k + 2 s = 1 n 2 r i ' t q s q ˙ s + 2 r i ' t 2 (25)

Hence, we have

δ G r ¨ ˜ i ' = s = 1 n r i ' q s δ G q ¨ s (26)

Calculating the Gaussian variation of equation (11), we get

δ G Z r = δ G S r - i = 1 N ( F i + F e i I + F c i I ) δ G r ¨ ˜ i ' (27)

Notice that

δ G S r = s = 1 n S r q ¨ s δ G q ¨ s (28)

i = 1 N F i δ G r ¨ ˜ i ' = s = 1 n i = 1 N F i r i ' q s δ G q ¨ s = s = 1 n Q s δ G q ¨ s (29)

i = 1 N ( F e i I + F c i I ) δ G r ¨ ˜ i ' = i = 1 N [ - m i a o - m i ω ˙ × r i ' - m i ω × ( ω × r i ' ) - 2 m i ω × r ˙ ˜ i ' ] s = 1 n r i ' q s δ G q ¨ s = s = 1 n i = 1 N { - m i a o r i ' q s - m i ( ω ˙ × r i ' ) r i ' q s - m i [ ω × ( ω × r i ' ) ] r i ' q s - 2 m i ( ω × r ˙ ˜ i ' ) r i ' q s } δ G q ¨ s = s = 1 n [ - q s ( Π o + Π ω ) + Q s ω ˙ + Γ s ] δ G q ¨ s (30)

where

Π o = i = 1 N m i a o r i ' = m a o r c ' (31)

is the potential energy of a uniform force field[29], and

Π ω = 1 2 i = 1 N m i [ ω × ( ω × r i ' ) ] r i ' = - 1 2 i = 1 N m i ( ω × r i ' ) ( ω × r i ' ) = - 1 2 ω θ o ω (32)

is the potential energy of centrifugal forces, θo=i=1Nmi[(ri')2E-ri'ri'] is the inertia tensor. And

Q s ω ˙ = - i = 1 N m i ( ω ˙ × r i ' ) r i ' q s (33)

represents the generalized rotational inertia force and

Γ s = - i = 1 N 2 m i ( ω × r ˙ ˜ i ' ) r i ' q s (34)

is the generalized gyroscopic force. By substituting formulas (28), (29) and (30) into formula (27), we get

δ G Z r = s = 1 n [ S r q ¨ s - Q s + q s ( Π o + Π ω ) - Q s ω ˙ - Γ s ] δ G q ¨ s (35)

Therefore, principle (7) gives

s = 1 n [ S r q ¨ s - Q s + q s ( Π o + Π ω ) - Q s ω ˙ - Γ s ] δ G q ¨ s = 0 (36)

Equation (36) is the Appell form of the Gauss principle of relative motion dynamics in generalized coordinates.

Now, two alternative forms of the principle are derived: the Lagrange form and the Nielsen form. First of all, we give some formulas for the subsequent derivation. From equations (24) and (25), we can easily obtain

r ¨ ˜ i ' q ¨ s = r ˙ ˜ i ' q ˙ s = r i ' q s (37)

r ¨ ˜ i ' q ˙ s = 2 k = 1 n 2 r i ' q s q k q ˙ k + 2 2 r i ' t q s = r ˙ ˜ i ' q s (38)

d ˜ d t r i ' q s = r ˙ ˜ i ' q s (39)

By the relation between the absolute derivative and the relative derivative of a vector, for any vector A, there is

d d t A = d ˜ d t A + ω × A (40)

Thus, we have

d d t r ˙ ˜ i ' = d ˜ d t r ˙ ˜ i ' + ω × r ˙ ˜ i ' = r ˙ ˜ i ' + ω × r ˙ ˜ i ' (41)

d d t r i ' q s = d ˜ d t r i ' q s + ω × r i ' q s = r ˙ ˜ i ' q s + ω × r i ' q s (42)

Secondly, denote Tr as the kinetic energy of relative motion, i.e.,

T r = 1 2 i = 1 N m i r ˙ ˜ i ' r ˙ ˜ i ' (43)

then we have

T r q ˙ s = i = 1 N m i r ˙ ˜ i ' r ˙ ˜ i ' q ˙ s = i = 1 N m i r ˙ ˜ i ' r i ' q s (44)

By using equations (41) and (42), we can obtain

d d t T r q ˙ s = i = 1 N m i d d t r ˙ ˜ i ' r i ' q s + i = 1 N m i r ˙ ˜ i ' d d t r i ' q s = i = 1 N m i ( r ¨ ˜ i ' + ω × r ˙ ˜ i ' ) r i ' q s + i = 1 N m i r ˙ ˜ i ' ( r ˙ ˜ i ' q s + ω × r i ' q s ) = i = 1 N m i r ¨ ˜ i ' r i ' q s + i = 1 N m i r ˙ ˜ i ' r ˙ ˜ i ' q s (45)

and

T r q s = i = 1 N m i r ˙ ˜ i ' r ˙ ˜ i ' q s (46)

Therefore, we have

d d t T r q ˙ s - T r q s = i = 1 N m i r ¨ ˜ i ' r i ' q s (47)

From equations (10) and (37), we get

S r q ¨ s = i = 1 N m i r ¨ ˜ i ' r ¨ ˜ i ' q ¨ s = i = 1 N m i r ¨ ˜ i ' r i ' q s (48)

By comparing formula (47) and formula (48), principle (36) can be expressed as

s = 1 n [ d d t T r q ˙ s - T r q s - Q s + q s ( Π o + Π ω ) - Q s ω ˙ - Γ s ] δ G q ¨ s = 0 (49)

Equation (49) is the Lagrange form of the Gauss principle of relative motion dynamics in generalized coordinates.

Calculating the time derivative of Tr, we get

d d t T r = i = 1 N m i r ˙ ˜ i ' d d t r ˙ ˜ i ' = i = 1 N m i r ˙ ˜ i ' ( r ¨ ˜ i ' + ω × r ˙ ˜ i ' ) = i = 1 N m i r ˙ ˜ i ' r ¨ ˜ i ' (50)

Here, the following relationship is applied, i.e.,

r ˙ ˜ i ' ( ω × r ˙ ˜ i ' ) = ω ( r ˙ ˜ i ' × r ˙ ˜ i ' ) = 0 (51)

Taking the partial derivative of ddtTr with respect to q˙s, we get

q ˙ s d d t T r = i = 1 N m i r ¨ ˜ i ' r ˙ ˜ i ' q ˙ s + i = 1 N m i r ¨ ˜ i ' q ˙ s r ˙ ˜ i ' = i = 1 N m i r ¨ ˜ i ' r i ' q s + 2 i = 1 N m i r ˙ ˜ i ' r ˙ ˜ i ' q s (52)

From formula (52) and formula (46), we obtain

q ˙ s d d t T r - 2 T r q s = i = 1 N m i r ¨ ˜ i ' r i ' q s (53)

By comparing formula (53) and formula (48), principle (36) can also be expressed as

s = 1 n [ q ˙ s d d t T r - 2 T r q s - Q s + q s ( Π o + Π ω ) - Q s ω ˙ - Γ s ] δ G q ¨ s = 0 (54)

Equation (54) is the Nielsen form of the Gauss principle of relative motion dynamics in generalized coordinates.

5 Dynamical Equations of Relative Motion

For a holonomic system, δGq¨s is independent and arbitrary, so from principle (36), we get

S r q ¨ s = Q s - q s ( Π o + Π ω ) + Q s ω ˙ + Γ s (55)

This is the Appell equation of relative motion dynamics, and s=1,2,,n. From principle (49), we get

d d t T r q ˙ s - T r q s = Q s - q s ( Π o + Π ω ) + Q s ω ˙ + Γ s (56)

This is the Lagrange equation of relative motion dynamics. From principle (54), we get

q ˙ s d d t T r - 2 T r q s = Q s - q s ( Π o + Π ω ) + Q s ω ˙ + Γ s (57)

This is the Nielsen equation of relative motion dynamics. For a nonholonomic system, let g ideal two-sided nonholonomic constraints be

ϕ α ( t , q s , q ˙ s ) = 0 ,   α = 1,2 , , g (58)

By differentiating equation (58), we obtain

ϕ α q ˙ s q ¨ s + ϕ α q s q ˙ s + ϕ α t = 0 (59)

Then we have

ϕ α q ˙ s δ G q ¨ s = 0 (60)

From the Gauss principle (36) and formula (60) in Appell form, using the Lagrange multiplier method, we get

S r q ¨ s = Q s - q s ( Π o + Π ω ) + Q s ω ˙ + Γ s + μ α ϕ α q ˙ s (61)

where μα is the Lagrange multiplier, s=1,2,,n. Formula (61) is the Appell equation with multipliers in generalized coordinates for nonholonomic systems in relative motion. From the Lagrange form of Gauss principle (49) and formula (60), we get

d d t T r q ˙ s - T r q s = Q s - q s ( Π o + Π ω ) + Q s ω ˙ + Γ s + μ α ϕ α q ˙ s (62)

This is the Lagrange equation with multipliers in generalized coordinates for nonholonomic systems in relative motion, also known as Routh equation. From the Nielsen form of the Gauss principle (54) and formula (60), we get

q ˙ s d d t T r - 2 T r q s = Q s - q s ( Π o + Π ω ) + Q s ω ˙ + Γ s + μ α ϕ α q ˙ s (63)

This is the Nielsen equation with multipliers in generalized coordinates for nonholonomic systems in relative motion.

6 Examples

Example 1 A physical pendulum with mass m is suspended at point O on a given block AB, as shown in Fig. 1. Let block AB do circumferential translation with radius l0, whose motion is determined by the angle θ and known as θ=θ(t). The angle describes the position of the pendulum relative to ABφ, and the distance from its center of mass C to O is R, and its radius of gyration to C is ρ. Try to establish the dynamic equation of relative motion using the Gauss principle.

thumbnail Fig. 1 A physical pendulum in relative motion

In this example, the carrier is the block AB, and the carried body is the physical pendulum. The carrier's motion is circumferential translation, and the relative motion of the carried body is fixed axis rotation around the axis Oξ. The acceleration energy of the relative motion of the pendulum is

S r = 1 2 J ξ ε r 2 + = 1 2 ( J C + m R 2 ) ε r 2 + = 1 2 m ( ρ 2 + R 2 ) φ ¨ 2 + (64)

The active force is only gravity mg, and the moment to the axis Oξ is

M ξ = - m g R s i n φ (65)

Since the convected motion is translation, there is no Coriolis inertia force, and the convected inertia force is

F e I τ = m l 0 θ ¨ , F e I n = m l 0 θ ˙ 2 (66)

The moment of convected inertia force about Oξ is

M ξ e = - F e I n R s i n ( φ - θ ) - F e I r R c o s ( φ - θ ) (67)

Therefore, from formula (22), the compulsion function of relative motion is

Z r = 1 2 J ξ ε r 2 - ε r ( M ξ   + M ξ e I + M ξ c I ) + = 1 2 m ( ρ 2 + R 2 ) φ ¨ 2 + [ m g R s i n φ + m l 0 R θ ˙ 2 s i n ( φ - θ ) + m l 0 R θ ¨ c o s ( φ - θ ) ] φ ¨ + (68)

To calculate the Gaussian variation δGZr and set it to zero, we get

δ G Z r = m ( ρ 2 + R 2 ) φ ¨ δ G φ ¨ + m l 0 R θ ˙ 2 s i n ( φ - θ ) δ G φ ¨ + m l 0 R θ ¨ c o s ( φ - θ ) δ G φ ¨ + m g R s i n φ δ G φ ¨ = 0 (69)

Due to the arbitrariness of δGφ¨, we get

m ( ρ 2 + R 2 ) φ ¨ + m l 0 R θ ˙ 2 s i n ( φ - θ ) + m l 0 R θ ¨ c o s ( φ - θ ) + m g R s i n φ = 0 (70)

This represents the differential equation governing the relative motion of a physical pendulum. Equation (70) is consistent with the results obtained using the Lagrange equation in Ref. [29].

Example 2 As shown in Fig. 2, a uniform rod AB with mass m and length l has one end, A, sliding along the vertical fixed axis Oz and the other end, B, sliding along the horizontal axis Ox. In contrast, Ox rotates around Oz at a uniform angular velocity ω. The point B is connected to the spring BD, and D is fixed on the Ox axis. Let θ indicate the angle between AB and the plumb line, when θ=0, the spring has its original length. Suppose that the spring stiffness is k, friction is ignored, and 0θπ2, find the dynamic equation of relative motion.

thumbnail Fig. 2 A uniform rod AB in relative motion

In this scenario, the carrier rotates at a uniform angular speed around Oz. The relative motion of the carried body AB is planar. With θ the generalized coordinate, the acceleration energy of the relative motion of rod AB is

S r = 1 2 m a C r 2 + 1 2 J C ε r 2 (71)

where JC=112ml2. Since

x C = 1 2 l s i n θ ,   z C = 1 2 l c o s θ (72)

Taking the time derivative of (72) twice, we have

x ¨ C = 1 2 l θ ¨ c o s θ - 1 2 l θ ˙ 2 s i n θ ,   z ¨ C = - 1 2 l θ ¨ s i n θ - 1 2 l θ ˙ 2 c o s θ (73)

Hence, we have

a C r 2 = x ¨ C 2 + z ¨ C 2 = 1 4 l 2 θ ¨ 2 + (74)

Substituting equation (74) into equation (71) and noting that εr=θ¨, we get

S r = 1 2 m 1 4 l 2 θ ¨ 2 + 1 2 1 12 m l 2 θ ¨ 2 + = 1 6 m l 2 θ ¨ 2 + (75)

Suppose we take a small segment dli on AB at a distance li from end A; then its mass is dmi=mldli. The coordinates in the moving coordinate system Oxyz are

x i = l i s i n θ ,   y i = 0 ,   z i = ( l - l i ) c o s θ (76)

Then we have

x ˙ i = l i θ ˙ c o s θ ,   y ˙ i = 0 ,   z ˙ i = - ( l - l i ) θ ˙ s i n θ (77)

x ¨ i = l i ( θ ¨ c o s θ - θ ˙ 2 s i n θ ) ,   y ¨ i = 0 , z ¨ i = - ( l - l i ) ( θ ¨ s i n θ + θ ˙ 2 c o s θ ) (78)

Now let us calculate the relevant terms in the compulsion function formula (11), and we get

ω = ω k (79)

r ˙ ˜ i ' = l i θ ˙ c o s θ i - ( l - l i ) θ ˙ s i n θ k (80)

r ¨ ˜ i ' = l i ( θ ¨ c o s θ - θ ˙ 2 s i n θ ) i - ( l - l i ) ( θ ¨ s i n θ + θ ˙ 2 c o s θ ) k (81)

F i = - d m i g k (82)

F e i I = - d m i a e i = d m i x i ω 2 i = d m i l i ω 2 s i n θ i (83)

F c i I = - 2 d m i ω × r ˙ ˜ i ' = - 2 d m i l i ω θ ˙ c o s θ j (84)

In addition, the elastic force FBD of the spring and the relative acceleration of its action point B are

F B D = - x B i = - k l s i n θ i (85)

r ¨ ˜ B ' = l ( θ ¨ c o s θ - θ ˙ 2 s i n θ ) i (86)

Thus, we have

i = 1 N ( F i   + F e i I + F c i I ) r ¨ ˜ i ' = i = 1 N { ( - d m i g k i   + d m i l i ω 2 s i n θ i - 2 d m i l i ω θ ˙ c o s θ j ) [ l i ( θ ¨ c o s θ - θ ˙ 2 s i n θ ) i - ( l - l i ) ( θ ¨ s i n θ + θ ˙ 2 c o s θ ) k ] }      + ( - k l s i n θ i ) ( l θ ¨ c o s θ - l θ ˙ 2 s i n θ ) i = 0 l [ m l ω 2 l i 2 s i n θ ( θ ¨ c o s θ - θ ˙ 2 s i n θ ) + m l g ( l - l i ) ( θ ¨ s i n θ + θ ˙ 2 c o s θ ) ] d l i - k l 2 s i n θ ( θ ¨ c o s θ - θ ˙ 2 s i n θ ) = ( 1 3 m l 2 ω 2 - k l 2 ) ( θ ¨ c o s θ - θ ˙ 2 s i n θ ) s i n θ + 1 2 m g l ( θ ¨ s i n θ + θ ˙ 2 c o s θ ) (87)

By substituting equations (75) and (87) into equation (11), we get

Z r = 1 6 m l 2 θ ¨ 2 - ( 1 3 m l 2 ω 2 - k l 2 ) ( θ ¨ c o s θ - θ ˙ 2 s i n θ ) s i n θ - 1 2 m g l ( θ ¨ s i n θ + θ ˙ 2 c o s θ ) + (88)

To calculate the Gaussian variation δGZr and set it to zero, we get

δ G Z r = [ 1 3 m l 2 θ ¨ - ( 1 3 m l 2 ω 2 - k l 2 ) s i n θ c o s θ - 1 2 m g l s i n θ ] δ G θ ¨ = 0 (89)

Due to the arbitrariness of δGθ¨, we get

1 3 m l 2 θ ¨ - ( 1 3 m l 2 ω 2 - k l 2 ) s i n θ c o s θ - 1 2 m g l s i n θ = 0 (90)

i.e.,

θ ¨ - ( ω 2 - 3 k m ) s i n θ c o s θ - 3 g 2 l s i n θ = 0 (91)

This is the dynamic equation of the relative motion of rod AB. It is consistent with the results obtained using the Lagrange equation in Ref. [29].

7 Conclusion

Complex mechanical systems, including the carrier and the carried bodies, are ubiquitous, so their study is significant. Using the theory of analytical mechanics to study the relative motion dynamics of complex systems not only has the unity of expression form but also shows the superiority of analytical mechanics in solving the problems of complex system dynamics. Unlike other differential variational principles, such as d'Alembert-Lagrange's or Jourdain's principle, Gauss principle is an extreme value principle from which the motion of a system can be directly obtained. The work conducted in this article includes the following aspects:

① The Gauss principle of relative motion dynamics and its least compulsion principle were established. Based on the dynamic equation of relative motion and the concept of virtual displacement in acceleration space, the Gauss principle for relative motion dynamics was presented. The compulsion function of relative motion was constructed, and it was proved that real motion makes the compulsion function yield its extreme value under Gaussian variation.

② The formulation of acceleration energy and compulsion function of relative motion was presented. The acceleration energy and compulsion function of relative motion were obtained when the carried rigid body was in planar motion.

③ The Appell, Lagrange, and Nielsen forms in generalized coordinates of the Gauss principle of relative motion were derived. According to the above documents of the Gauss principle, using the Lagrange multiplier method, we established the Appell equation, the Lagrange equation, and the Nielsen equation with multipliers for relative motion dynamics.

References

  1. Mei F X, Wu H B, Li Y M. A Brief History of Analytical Mechanics[M]. Beijing: Science Press, 2019(Ch). [Google Scholar]
  2. Chen B. Analytical Dynamics[M]. 2nd Ed. Beijing: Peking University Press, 2012(Ch). [Google Scholar]
  3. Mei F X. Analytical Mechanics Ⅱ[M]. Beijing: Beijing Institute of Technology Press, 2013(Ch). [Google Scholar]
  4. Udwadia F E, Kalaba R E. Analytical Dynamics — A New Approach[M]. New York: Cambridge University Press, 2008. [Google Scholar]
  5. Попов Е П. Operating Robot Dynamics and Algorithm[M]. Yu L J, Chen X J trans. Beijing: Mechanical Industry Press, 1983(Ch). [Google Scholar]
  6. Liu Y Z, Pan Z K, Ge X S. Dynamics of Multibody Systems[M]. 2nd Ed. Beijing: Higher Education Press, 2014(Ch). [Google Scholar]
  7. Liu Y Z. Dynamic modeling of multi-body system based on Gauss's principle[J]. Chinese Journal of Theoretical and Applied Mechanics, 2014, 46(6): 940-945(Ch). [Google Scholar]
  8. Liu Y Z. Dynamical modeling of a net system of rods based on Gauss's principle[J]. Journal of Dynamics and Control, 2018, 16(4): 289-294(Ch). [Google Scholar]
  9. Yao W L, Liu Y P, Yang L S. Dynamic modeling of nonideal system based on Gauss's principle[J]. Chinese Journal of Theoretical and Applied Mechanics, 2020, 52(4): 945-953(Ch). [Google Scholar]
  10. Yao W L, Yang L S, Song K W, et al. Optimization method for dynamics of non-holonomic system based on Gauss's principle[J]. Acta Mechanica Sinica, 2020, 36(5): 1133-1141. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  11. Yao W L, Yang L S, Guo M M. Gauss optimization method for the dynamics of unilateral contact of rigid multibody systems[J]. Acta Mechanica Sinica, 2021, 37(3): 494-506. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  12. Orsino R M M. Extended constraint enforcement formulations for finite-DOF systems based on Gauss's principle of least constraint[J]. Nonlinear Dynamics, 2020, 101(4): 2577-2597. [CrossRef] [Google Scholar]
  13. Yao W L, Dai B Q. Gauss principle of least constraint in generalized coordinates and its generalization[J]. Mechanics in Engineering, 2014, 36(6): 779-782, 785(Ch). [Google Scholar]
  14. Yang L S, Yao W L, Xue S F. Application of particle swarm optimization on the multi-body system dynamics with singular positions[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2021, 57(5): 795-803(Ch). [Google Scholar]
  15. Xue Y, Liu Y Z, Chen L Q. On analytical mechanics for a super-thin elastic rod[J]. Chinese Journal of Theoretical and Applied Mechanics, 2005, 37(4): 485-493(Ch). [Google Scholar]
  16. Xue Y, Weng D W. Gauss principle for a super-thin elastic rod dynamics[J]. Acta Physica Sinica, 2009, 58(1): 34-39(Ch). [Google Scholar]
  17. Liu Y Z, Xue Y. Dynamical model of Cosserat elastic rod based on Gauss principle[J]. Acta Physica Sinica, 2015, 64(4): 044601(Ch). [CrossRef] [Google Scholar]
  18. Xue Y, Qu J L, Chen L Q. Gauss principle of least constraint for Cosserat growing elastic rod dynamics[J]. Applied Mathematics and Mechanics, 2015, 36(7): 700-709(Ch). [Google Scholar]
  19. Зегжда СА, Солтаханов ШХ, Юшков МП. Equations of Motion for Nonholonomic Systems and Variational Principles of Mechanics: A New Class of Control Problems[M]. Mei F X trans. Beijing: Beijing Institute of Technology Press, 2007(Ch). [Google Scholar]
  20. Li B L, Song F P. On Gauss's principle of least constraint for impulsive motion[J]. Journal of Hunan University, 1995, 22(4): 23-28(Ch). [Google Scholar]
  21. Ivanov A P. On the variational formulation of the dynamics of systems with friction[J]. Regular and Chaotic Dynamics, 2014, 19(1): 100-115. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  22. Wang L S, Pao Y H. Jourdain's variational equation and Appell's equation of motion for nonholonomic dynamical systems[J]. American Journal of Physics, 2003, 71(1): 72-82. [NASA ADS] [CrossRef] [Google Scholar]
  23. Yan C C. Hamilton's principle and Schrodinger's equation derived from Gauss's principle of least squares[J]. Foundations of Physics Letters, 2000, 13(1): 79-87. [CrossRef] [MathSciNet] [Google Scholar]
  24. Zegzhda S A, Soltakhanov S K. Application of the generalized Gaussian principle to the problem of damping vibrations of mechanical systems[J]. Journal of Computer and Systems Sciences International, 2010, 49(2): 186-191. [CrossRef] [MathSciNet] [Google Scholar]
  25. Lewis A D. The geometry of the Gibbs-Appell equations and Gauss's principle of least constraint[J]. Reports on Mathematical Physics, 1996, 38(1): 11-28. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  26. Kalaba R E, Natsuyama H H, Udwadia F E. An extension Gauss's principle of least constraint[J]. International Journal of General Systems, 2004, 33(1): 63-69. [CrossRef] [MathSciNet] [Google Scholar]
  27. Zhang Y, Chen X Y. The generalized Gauss principle for mechanical system with variable mass and its generalization to higher order nonholonomic systems[J]. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(10): 2883-2891(Ch). [Google Scholar]
  28. Zhang Y, Song C J, Zhai X H. Generalized Gauss principle of least compulsion for variable acceleration dynamical systems[J]. Chinese Journal of Theoretical and Applied Mechanics, 2023, 55(5): 1174-1180(Ch). [Google Scholar]
  29. Mei F X. Analytical Mechanical Special Problem[M]. Beijing: Beijing Institute of Technology Press, 1988(Ch). [Google Scholar]
  30. Markeyev A P. Approximate equations of rotational motion of a rigid body carrying a movable point mass[J]. Journal of Applied Mathematics and Mechanics, 2013, 77(2): 137-144. [CrossRef] [MathSciNet] [Google Scholar]
  31. Prioroc C L, Mikkola S. Simple algorithms for relative motion of satellites[J]. New Astronomy, 2015, 34: 41-46. [NASA ADS] [CrossRef] [Google Scholar]
  32. Vepa R. Feedback tracking control of optimal reference trajectories for spacecraft relative motion[J]. Advances in Space Research, 2022, 69(9): 3478-3489. [NASA ADS] [CrossRef] [Google Scholar]
  33. Burnett E R, Schaub H. Geometric perspectives on fundamental solutions in the linearized satellite relative motion problem[J]. Acta Astronautica, 2022, 190: 48-61. [NASA ADS] [CrossRef] [Google Scholar]
  34. Whittaker E T. A Treatise on the Analytical Dynamics of Particles and Rigid Bodies[M]. Fourth Ed. Cambridge: Cambridge University Press, 1952. [Google Scholar]
  35. Lure A I. Analytical Mechanics[M]. Moscow: GIFML, 1961 (in Russian). [Google Scholar]
  36. Mei F X, Wu H B. Dynamics of Constrained Mechanical Systems[M]. Beijing: Beijing Institute of Technology Press, 2009. [Google Scholar]
  37. Liu G L, Qiao Y F, Zhang J F, et al. Relative motion dynamics of variable mass non-holonomic mechanical systems[J]. Acta Mechanica Sinica, 1989, 21(6): 742-748(Ch). [Google Scholar]
  38. Chen X W, Luo S K. Integration method for the dynamics equation of relative motion of variable mass nonlinear nonholonomic system[J]. Applied Mathematics and Mechanics, 1998, 19(5): 479-488. [Google Scholar]
  39. Luo S K. Integral theory for the dynamics of nonlinear nonholonomic system in noninertial reference frames[J]. Applied Mathematics and Mechanics, 1993, 14(10): 907-918. [Google Scholar]
  40. Zhang X W. The universal forms of the dynamic equations of holonomic mechanical system in relative motion[J]. Acta Physica Sinica, 2006, 55(6): 2669-2675(Ch). [Google Scholar]
  41. Xie Y L, Jia L Q, Luo S K. Special Lie symmetry and Hojman conserved quantity of Appell equations in a dynamical system of relative motion[J]. Chinese Physics B, 2011, 20(1): 010203. [CrossRef] [MathSciNet] [Google Scholar]
  42. Xie Y L, Jia L Q, Yang X F. Lie symmetry and Hojman conserved quantity of Nielsen equation in a dynamical system of the relative motion[J]. Acta Physica Sinica, 2011, 60(3): 030201(Ch). [Google Scholar]
  43. Mei F X, Wu H B. Lagrange symmetry for a dynamical system of relative motion[J]. Acta Physica Sinica, 2009, 58(9): 5919-5922(Ch). [CrossRef] [MathSciNet] [Google Scholar]
  44. Zhang Y, Mei F X. Algebraic structure of the dynamical equations of holonomic mechanical system in relative motion[J]. Journal of Beijing Intitute of Technology (English Edition), 1998, 7(1): 12-18. [Google Scholar]

All Figures

thumbnail Fig. 1 A physical pendulum in relative motion
In the text
thumbnail Fig. 2 A uniform rod AB in relative motion
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.