On dynamics of a rigid body moving on a horizontal plane by means of motion of an internal particle

Boris Bardin1 , Aleksandr Panev2

1, 2Moscow Aviation Institute (National Research University), Moscow, Russia

2Corresponding author

Vibroengineering PROCEDIA, Vol. 8, 2016, p. 135-141.
Received 31 July 2016; accepted 3 September 2016; published 7 October 2016

Copyright © 2016 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Abstract.

We deal with mechanical system consisting of a body moving on a horizontal plane and a particle which moves within the body in a circle. We suppose that dry friction acts between the plane and the body and assume that the body moves without jumps on the plane. We study the dynamics of the above system in detail. In particular, we show that, depending on parameters values, motion of the body approaches to one of three periodic modes.

Keywords: periodic motion, asymptotic stability, rigid body with moving particle, dry friction, mobile robots.

1. Introduction

Investigation of dynamic of rigid body, which moves by means of relative motions of its internal components is a relevant engineering problem. The advantage of the so-called vibrating mobile robots is that they are simple in design, i.e. do not require external elements such as wheels or track chains and can be produced as isolated capsules. This circumstance makes them resistant to aggressive external environment. That is why they can be used both on hard surfaces and in liquids. In papers [1-9] dynamics and control of a rigid body moving on a rough horizontal plane by means of relative motions of the internal components has been studied.

In this paper we consider the mechanical system consisting of a rigid body of mass M and a particle of mass m moving inside the body in a circle of radius R. The center of the circle is located in the body mass center and the angular velocity ω of the relative circular motion of the particle is a constant.

Fig. 1. Coordinate systems

We consider a rectilinear translation of the body in a rough horizontal plane. In this case the position of the body can be defined by one variable X, which is the X-coordinate of its center of mass in fixed coordinate system OXY (see Fig.1).

Let us introduce dimensionless time t so that the period of particle relative motion is equal to 2π. Without loss of generality we suppose that at moment t= 0 the particle, moving counterclockwise, passes the lowest point of the circle. Then the time t can be regarded as angle defining the position of the particle in moving (with the body) coordinate system O1ξη (see Fig. 1).

Taking into account gravity and dry friction forces we have the following equation of motion [7]:

(1)
d u d t = sin t + f c ,

where:

(2)
f c = - k μ + cos t s i g n   u , u 0 , - sin t , u = 0   a n d   sin t k μ + cos t , - k μ + cos t s i g n sin t , u = 0   a n d   sin t > k μ + cos t ,

and u= dx/dt. The dimensionless coordinate x and the parameter μ, are introduced by formulas:

(3)
x = R m X M + m ,           μ = M + m g R m ω 2 ,

where g is the gravitational acceleration, k is the coefficient of dry friction.

We also assume that the body can start moving from state of rest, i.e. k21/(µ2-1) as well as the body moves without jumps on the horizontal plane, i.e. μ1. In [7] the partial solution of Eq. (1) with initial condition u0=0 has been studied in detail for all possible parameters values. In this paper we study solutions of Eq. (1) with arbitrary initial conditions. In Sections 3 and 4 we show that the above partial solution plays an important role in dynamics of the system considered. In particular, it has properties of a limit cycle.

2. Zones of friction deceleration and periodic motions

First we note that at ut0 the right hand side of Eq. (1) is equal to zero at the moments of time, when the horizontal component of inertia force applied to the particle and the dry friction force have opposite directions, but their absolute values are equal. In the interval [0; 2π] such moments of time can be calculated as follows [7]:

(4)
t 1 = 2 a r c t g   1 - 1 - k 2 µ 2 - 1 k µ - 1 ,     t 2 = 2 a r c t g   1 + 1 - k 2 µ 2 - 1 k µ - 1 ,
t 3 = 2 π - t 2 ,       t 4 = 2 π - t 1 .

The positions of the particle corresponding to the moments t1, t2, t3, t4   are shown in Fig. 2. In the intervals t2, t3 and  (t4, t1) direction of the body’s acceleration is opposite to the velocity of the body. The above intervals we call upper and lower zones of friction deceleration respectively. They play a significant role for the analysis of the body motion. In particular, if the body stops, when the particle is in one of them, then the body remains at rest until the particle leaves the zone of friction deceleration. This phenomenon we call the body sticking.

Fig. 2. Zones of friction deceleration

It was shown in [7] that for any possible values of parameters k and μ there exists motion of body with 2π-periodically changing velocity. Depending on the parameters such a periodic motion can have a quite different nature. The area of possible values of parameters k and μ is separated into three domains I, II and III corresponding to different types of the periodic motions (see Fig. 3).

Domain I is defined by the following inequalities:

(5)
k 2 1 µ 2 - 1 ,     μ k 2 + 1 - k 2 μ 2 - 1 μ k k 2 + 1 - a r c c t g   k < 0 .

In domain I the solution of Eq. (1) with initial condition ut1=0 is periodic. This solution describes the motion of the body with stops both in the upper and in the lower zones of friction deceleration. Moreover, the coordinate xt is also a 2π-periodic function in this case.

Domain II is defined by inequalities:

(6)
μ k 2 + 1 - k 2 μ 2 - 1 μ k k 2 + 1 - a r c c t g   k 0 ,
cos t 1 + cos t 1 + sin t 1 + 1 - cos t 1 - π k μ 2 μ π k μ .

In this domain the solution of Eq. (1) with initial condition ut1=0 is also 2π-periodic. However, the coordinate xt is not periodic function in this case. The body moves with a stop only in the lower zone of friction deceleration such that xt+2π-xt>0.

Domain III is defined by the inequalities:

(7)
cos t 1 + cos t 1 + sin t 1 + 1 - cos t 1 - π k μ 2 μ π k μ ,     μ > 1 .

In domain III the solution of Eq. (1) with initial condition ut1=0 is not periodic, but it approaches a periodic limit cycle, which describes motion of the body without stopping in any zones of friction deceleration.

Fig. 3. Domains of periodic motion; horizontal axis – k, vertical axis – μ

3. Motion of the body in domains I and II

Let the parameters belong to domain I. Consider the solution of Eq. (1) with initial condition ut1=u1, where u1 is a real positive number. For positive velocity u the Eq. (1) reads:

(8)
d u d t = sin t - k μ + cos t .

That is:

(9)
u t = u 1 + t 1 t [ sin τ - k μ + cos τ ] d τ ,

Eq. (9) describes the motion of the body until the velocity u is positive, i.e. the motion in the interval [t1,t'], where ut'=0. It can be shown that sint-kµ+cost>0 in the interval [t1,t2] and sint-kµ+cost<0 in the interval [t2;t1+2π]. It means that, function ut increases in the interval [t1,t2] and decreases in the interval [t2;t1+2π].

On the other hand, from Eq. (9) we have:

(10)
u t 1 + 2 π = u 1 - 2 π k μ .

Thus, if u1>2πkμ, then in the interval [t1;t1+2π] the body moves in positive direction without stopping. The velocity of the body decreases by 2πkμ for the period 2π.

Let us now suppose that u1<2πkμ. In this case the body stops in the interval [t2;t1+2π], i.e. t'<t1+2π. Note that t'>t2, because function ut increases in the interval [t1,t2].

If the body stops in the upper or in the lower zone of friction deceleration, then it remains at rest till the moment t3 or t1+2π respectively. In both cases ut1+2π=0, that is for t> t' the body moves with 2π-periodically changing velocity.

It remains to consider the case when the body stops in the interval [t3,t4]. Thus, we suppose that t3<t'<t4 and consider the motion of the body for t> t'. In this case Eq. (1) takes the following form:

(11)
d u d t = sin t + k μ + cos t .

That is:

(12)
u t = t ' t [ sin τ + k μ + cos τ ] d τ .

It can be shown that sint+kμ+cost<0 in the interval t3,t4 and sint+kμ+cost>0 in the interval [t4;t1+2π]. It means, that function ut defined by Eq. (11) decreases in the interval [t',t4] and increases in the interval [t4;t1+2π]. Thus, after the stop at the moment t=t' the body starts to move in opposite direction, i.e. with negative velocity until it stops again at some moment t=t''>t4. By using Eq. (12) and taking in to account the condition ut''=0 we have:

(13)
t ' t 4 [ sin τ + k μ + cos τ ] d τ = - t 4 t ' ' [ sin τ + k μ + cos τ ] d τ .

On the other hand, there exists the periodic motion in domain I such that at the moment t=t3 the body starts to move in negative direction of axis OX and stops in lower zone of friction deceleration at the moment  t=t**, where t4<t**<t1+2π, (see [7] for more details). For such a motion we have ut3=ut**=0. Taking this fact into account by analogy with Eq. (13) we have:

(14)
t 3 t 4 [ sin τ + k μ + cos τ ] d τ = - t 4 t * * [ sin τ + k μ + cos τ ] d τ .

Now we recall that sint+kμ+cost<0 in the interval [t3,t4], that is why:

(15)
t ' t 4 [ sin τ + k μ + cos τ ] d τ > t 3 t 4 [ sin τ + k μ + cos τ ] d τ .

By taking into account Eqs. (13), (14) we can rewrite the inequality (15) in the following form:

(16)
t 4 t ' ' [ sin τ + k μ + cos τ ] d τ < t 4 t * * [ sin τ + k μ + cos τ ] d τ .

Since sint+kμ+cost>0 in the interval [t4;t1+2π], it follows from Eq. (16) that t''< t**<t1+2π, i.e. the body stops in the lower zone of friction deceleration and remains at rest till the moment t1+2π. That is, for t>t1 the body moves with 2π-periodically changing velocity.

Thus, if the body starts to move with positive velocity, then it always passes in periodic mode of motion for a finite period of time. It is easy to see that the same result takes place in the case of negative initial velocity.

By analogy it can be also shown that in domain II the motion of the body has similar properties.

4. Motion of the body in domain III. Stability of periodic motion

Let us suppose that the parameters belong to domain III. In this case there exists a periodic mode of motion such that the body stops and reverses its direction of motion twice in the period 2π [7]. The stops of the body take place outside of both upper and lower zone of friction deceleration. That is the sticking phenomenon does not take place and there exists two moments of time t1*[t1;t2] and t2*[t3;t4], such that ut1*=ut2* =0. The values t1* and t2* satisfy the following system of equations [7]:

(17)
- cos t 2 * + cos t 1 * - k μ t 2 * - t 1 * - k sin t 2 * - sin t 1 * = 0 , cos t 2 * - cos t 1 * - k μ t 1 * - t 2 * + 2 π + k sin t 1 * - sin t 2 * = 0 .

It is clear that if values t1* and t2* are found, than the periodic mode of motion is fully determined. Numerical investigations have shown that in domain III the above system of Eq. (17) has a unique solution.

In order to prove the asymptotic stability of the above periodic mode of motion, we consider solution of Eq. (1) with the initial condition ut1*=ε, where ε is a small value. Without loss of generality we assume that ε>0. We will construct such a solution by integrating Eq. (1) in finite intervals of time, where the velocity of the body keeps a definite direction, i.e. the function ut does not change sign within these intervals.

At the first step we integrate Eq. (1) with the initial condition ut1*=ε in the interval t1*, t2*+τ, where value τ satisfies the condition ut2*+τ=0. As a result, we have:

(18)
- cos t 2 * + τ + cos t 1 * - k μ t 2 * + τ - t 1 * - k sin t 2 * + τ - sin t 1 * = - ε .

Let us expand the left hand side of Eq. (18) in a power series in τ:

(19)
sin t 2 * τ - k μ τ - k cos t 2 * τ + O τ 2 = - ε .

By solving Eq. (19) with respect to τ, we have:

(20)
τ - ε sin t 2 * - k μ - k cos t 2 * +   O ε 2 .

At the second step by integrating Eq. (1) in the interval (t2*+τ, t1*+2π) with the initial condition ut2*+τ=0, we get:

(21)
u t 1 * + 2 π = - cos t 1 * + cos t 2 * + τ + k μ t 1 * + 2 π - t 2 * - τ + k sin t 1 * - sin t 2 * + τ .

Now we expand left and right hand side of Eq. (21) in a power series in τ:

(22)
u t 1 * + 2 π = - sin t 2 * - k μ - k cos t 2 * τ + O τ 2 .

By substituting Eq. (20) in Eq. (22), we get:

(23)
u t 1 * + 2 π ε sin t 2 * + k μ + k cos t 2 * sin t 2 * - k μ - k cos t 2 * + O ε 2 .

Consider the expression:

(24)
Δ = sin t 2 * + k μ + k cos t 2 * sin t 2 * - k μ - k cos t 2 * .

First, we note that kμ+cost2*>0, because we have supposed μ>1. Then, we recall that t2*[t3;t4], hence sint2*<0. Now it is easy to see that the following inequality holds:

(25)
sin t 2 * - k μ + cos t 2 * > sin t 2 * + k μ + cos t 2 *

The last inequality can be rewritten in the form:

(26)
< 1 .

The process of integration of Eq. (1) can be continued. Finally, we get the following asymptotic estimate:

(27)
u t 1 + 2 π n = ε n + O ε 2 .

It follows from Eq. (27) that u(t1*+2πn)0 as n, i.e. the solution with the initial condition ut1*=ε approaches the periodic solution with the initial condition ut1*=0 at t. It means that the above periodic solution is asymptotically stable.

5. Conclusions

In this paper an analytical study of dynamics of a rigid body with a moving internal particle is performed. The body moves on a rough horizontal plane without jumps. It appears that for any possible parameter values there exists a periodic mode of motion, which is a limit cycle. The above periodic mode is asymptotically stable. The stability has global character, i.e. for all initial conditions the motion of the body approaches the periodic mode. Depending on parameters a finite or infinite period of time is necessary to pass into the periodic mode of motion.

Acknowledgements

The work is carried out at the cost of the Grant of the Russian Scientific Foundation (Project No. 14-21-00068) at the Moscow Aviation Institute (National Research University).

References

  1. Chernous’ko F. L. Analysis and optimization of the motion of a body controlled by means of a movable internal mass. Journal of Applied Mathematics and Mechanics, Vol. 70, Issue 6, 2006, p. 819-842. [CrossRef]
  2. Chernous’ko F. L. The optimal periodic motions of a two-mass system in a resistant medium. Journal of Applied Mathematics and Mechanics, Vol. 72, Issue 2, 2008, p. 116-125. [CrossRef]
  3. Bolotnik N. N., Zeidis I. M., Zimmermann K., Yatsun S. F. Dynamics of controlled motion of vibration-driven systems. Journal of Computer and Systems Sciences International, Vol. 45, Issue 5, 2006, p. 831-840. [CrossRef]
  4. Bolotnik N. N., Figurina T. Yu. Optimal control of the rectilinear motion of a rigid body on a rough plane by means of the motion of two internal masses. Journal of Applied Mathematics and Mechanics, Vol. 72, Issue 2, 2008, p. 126-135. [CrossRef]
  5. Bolotnik N. N., Figurina T. Yu., Chernous’ko F. L. Optimal control of the rectilinear motion of a two-body system in a resistive medium. Journal of Applied Mathematics and Mechanics, Vol. 76, Issue 1, 2012, p. 1-14. [CrossRef]
  6. Vartholomeos P., Vlachos K., Papadopoulos E. Analysis and motion control of a centrifugal-force microrobotic platform. IEEE Transactions on Automation Science and Engineering, Vol. 10, Issue 3, 2013, p. 545-553. [CrossRef]
  7. Bardin B. S., Panev A. S. Periodic motions of the body with a movable internal mass on a horizontal surface. Trudy MAI, Vol. 84, 2015. [CrossRef]
  8. Xiong Z., Jian X. Locomotion analysis of a vibration-driven system with three acceleration-controlled internal masses. Advances in Mechanical Engineering, Vol. 7, Issue 3, 2015. [CrossRef]
  9. Sakharov A. V. Rotation of a body with two movable internal masses on a rough plane. Journal of Applied Mathematics and Mechanics, Vol. 79, Issue 2, 2015, p. 132-141. [CrossRef]