CC BY
6Russian Journal of Nonlinear Dynamics, 2019, vol. 15, no. 2, pp. 171-178. Full-texts are available at http://nd.ics.org.ru DOI: 10.20537/nd190206
NONLINEAR PHYSICS AND MECHANICS
MSC 2010: 70E18, 70F40
The Rolling of a Homogeneous Ball with Slipping on a Horizontal Rotating Plane
T. B. Ivanova
This paper is concerned with the rolling of a homogeneous ball with slipping on a uniformly rotating horizontal plane. We take into account viscous friction forces arising when there is slipping at the contact point. It is shown that, as the coefficient of viscosity tends to infinity, the solution of the generalized problem on each fixed time interval tends to a solution of the corresponding nonholonomic problem.
Keywords: rotating surface, turntable, nonholonomic constraint, rolling ball, sliding, viscous friction
Introduction
This paper investigates the motion of a dynamically symmetric (in particular, homogeneous) ball rolling with .slipping on a uniformly rotating horizontal plane.
In the case where there is no slipping, the motion of the ball is subject to an inhomogeneous nonholonomic constraint, which corresponds to the condition that the velocities of the contacting points on the surface of the ball and the rotating plane be the same. In this case, the system is linear and admits explicit integration (see, e.g., [1-3] and references therein). The trajectories of the center of mass of the ball relative to a fixed reference system are circles, and the ball moves with constant angular velocity, which depends only on the angular velocity of the surface on which the ball rolls, and on the radius of inertia of the ball.
Received April 12, 2019 Accepted May 28, 2019
This work is supported by the Russian Science Foundation under grant 19-71-30012 and performed in Steklov Mathematical Institute of Russian Academy of Sciences.
Tatiana B. Ivanova tbesp@rcd.ru
Steklov Mathematical Institute, Russian Academy of Sciences, ul. Gubkina 8, Moscow, 119991 Russia
On the other hand, when the ball moves, for example, on a fast rotating plane, slipping may be observed at the contact point. In this case, at the point of contact with the plane, from the plane the ball is acted upon not only by the normal reaction, but also by the friction force.
In this paper, we assume that the friction force is viscous, that is, it depends linearly on the velocity of the contact point relative to the rotating plane and has an opposite direction.
Taking into account viscous friction forces leaves the problem linear and, as in the nonholo-nomic case, the system can be explicitly integrated. In addition, as the coefficient of viscosity tends to infinity, the solution on each fixed time interval tends to a solution of the corresponding nonholonomic problem.
We note that a similar result was proved in [4] for a Chaplygin sleigh executing inertial motion on a horizontal plane, and in [5] for a ball moving on the surface of a rotating cylinder. The correctness of such a passage to the limit from dissipative holonomic systems to conservative nonholonomic systems in general form is discussed in [6, 7].
The theoretical results obtained supplement the models of motion of a ball without slipping on a rotating plane [2]. The problem under consideration is of great importance for applications. On the one hand, the dynamics of the ball with slipping and viscous friction is described by a linear system of differential equations which is investigated by standard methods of linear algebra and of the theory of ordinary differential equations. On the other hand, the motion of the ball with slipping on a plane can be observed in an experiment. Using high-speed cameras and modern methods of processing experimental results, one can verify the adequacy of the friction models used.
1. Equations of motion of a ball moving with slipping on a rotating plane
Consider the motion of a homogeneous ball of mass m and radius R with slipping on a horizontal rough plane. The plane rotates with constant angular velocity Q about the axis coinciding with the external normal to the plane (Fig. 1).
Fig. 1. A schematic model of a ball on a rotating plane. In the left figure, the axis Ox is perpendicular to the figure plane and is directed toward the reader.
To describe the dynamics of the ball, we define a fixed (inertial) coordinate system Oxyz in such a way that the plane Oxy is parallel to the plane of motion of the center of mass of the ball and the axis Oz is directed along the external normal 7 to the plane of motion, g = (0,0,1).
In the chosen coordinate system, the position of the center of mass of the ball is defined by the radius vector r = (x,y,z), the free-fall acceleration vector has the form g = (0,0, -g), and the angular velocity of the rotating plane on which the ball rolls is O = (0,0, Q) (Fig. 1).
In this paper, we assume that the ball moves without loss of contact with the plane (z = R), with slipping at the point of contact. At the contact point, the ball is acted upon by the viscous friction force:
F f = —¡ou, (1.1)
where ¡0 is the coefficient of viscous friction, u is the velocity of the contact point on the ball relative to the rotating plane.
Let v denote the velocity of the center of mass of the ball relative to the fixed coordinate system:
v = r. (1.2)
Changes in the linear and angular momenta of the ball relative to its center are described by the Newton-Euler equations:
mV = 0u + mg + N, lu = ¡0Ry x u, (1.3)
where I is the central tensor of inertia of the ball, u is the angular velocity of the ball, and N = (0,0, N) is the reaction force acting from the plane and directed along the normal to it. The velocity u of the contact point located on the ball is given by
u = v + Rg x u — O x r. (1.4)
In what follows, we will write all variables and equations in dimensionless form. To this end, we take the radius of the ball R as the unit of distance, the quantity 1/S as the unit of time, and the mass of the ball m as the unit of mass, that is, in Eqs. (1.1)-(1.4) we make the following change of variables:
r V u w
-^r, tn^t, U, TZ^UJ,
R Rli Rli Si /-,
(1.5)
g N ¡¡0
Substituting the expression for u (1.4) into the system (1.2)-(1.3), we obtain a closed system of nine differential first-order equations (in dimensionless variables):
r = v,
V = —¡o (v + Y x u — y x r)+ g + N, (1.6)
^ x v — up + 7(7 • a;) + r — 7(7 • r)j,
W Jfe
where k = I/(mR2) and k = 2/5 for a homogeneous ball.
According to the condition of motion without loss of contact, Z = 0. Then, projecting the first equation of (1.6) onto the axis Oz, we find vz = 0, vz = 0. Thus, from the projection of the second equation of (1.6) onto the axis Oz we find the reaction force
N = g.
The projection of the third equation of (1.6) onto the axis Oz has the form
Wz = 0,
that is, the vertical projection of the angular velocity of the ball remains unchanged.
We finally find that the system of equations (1.6) reduces to a linear system of six differential first-order equations relative to the vector of the variables
Z = (x, y, Vx, Vy, Ux, Uy),
which can be represented in matrix form
(
0 I 0
Z = Lz, L= JJ -J I -J J
\k k J
(1.7)
where L is a 6x6 block matrix, I is a 2x2 identity matrix, and Jj = —£ij3-
Remark 1. Since the vertical projection of the angular velocity uz = (w • 7) remains unchanged and does not appear in the system (1.7), the behavior of the other quantities describing the motion of the ball does not depend on its value. This implies, in particular, that the motion of a rubber ball with the additional constraint uz = (w • 7) = 0 [8] is a particular case of the system of interest. The situation where the phase flow of one (nonholonomic) system is a particular case of a more general (Hamiltonian) system on a specific integral surface is also considered in [9].
2. Explicit integration
The eigenvalues Xj of the matrix L satisfy the equation
X2P4 (X) = 0, P4 (X) =
X
+ 4
X
+ 4
X
+
J J \J J \JJ
where we have introduced the following notation to abbreviate the formulae:
Jo(k +1) Jo 4k2
J=
2k
>0, e = ^ =
j2 jo(k + 1)2
> 0.
(2.1)
(2.2)
Equation (2.1) has two zero roots which correspond to two linearly independent eigenvectors forming the kernel of the matrix L:
ei = (1, 0, 0, 0,1, 0), e2 = (0,1, 0, 0, 0,1).
In this case there exists a family of fixed points of the system (1.7)
zc = xcei + yce2 = (xc,yc, 0, 0,xc,yc), (2.3)
where xc and yc are independent parameters defining this family.
For fixed (xc,yc) a fixed point of this family corresponds to a state of the system in which the center of mass of the ball does not move relative to the fixed reference system. In this case, the plane rotates under it so that the ball rolls without slipping. That is to say, the ball has a constant angular velocity that (in dimensional form) is related to the angular velocity Q of rotation of the plane by
xc Q ycQ
IT' W2 = 1T
Ui
4
3
2
The nonzero roots of the characteristic polynomial (2.1), i.e., the roots of the polynomial P4(A), are pairwise complex conjugate and have the form
Ai,2 = H ( — 1 + Vli 'it) = OL\ ± i/3,
A3,4 = ¡j, 1 — \/l ± is) = «2 ± if],
(2.4)
where f is a complex unit. The real parts a1,a2 and the imaginary part / of the roots of the polynomial P4 (A) have the form
(2.5)
+ £2 - 1.
We see that, according to (2.2) and (2.5), the following inequalities are satisfied (see Fig. 2):
a2 < 0 < ai. (2.6)
Thus, the fixed point zc (2.3) is an unstable focus.
Fig. 2. Dependence of the real parts (ai, a2) and the imaginary part (fi) of the eigenvalues (2.4) on the coefficient vo at k = 2/5.
According to (2.4), the solution z(t) of the system (1.7) is the sum:
z(t) = zo + z+ + z _,
zo = ^oei + yoe2, z+ = e"1t (clei^t + C2e_^ , z_ = e"2t (c3eifit + C4,
(2.7)
where x0,y0 are the constants defined from the initial conditions, Cj,j = 1... 4 are the complex constant vectors defined in terms of the eigenvectors of the matrix L and from the initial conditions z(0). Within the framework of this paper we do not present them explicitly, since they are cumbersome in form. We present only the expressions for x0 ,y0, which are expressed in terms of the initial values of linear vx(0),vy(0) and angular (0),wy(0) velocities:
kux(0) - Vy(0) kujyj0) - vx(0)
0 =-k-' Vo =-k-' ( ^
According to (2.6) and (2.7), the solution of the system (1.7) is the sum of the constant term (z0) and the oscillating terms (z+ and z_). The frequency of oscillations of these terms
is the same and equals /, the amplitudes are monotonie increasing (eai t ) and decreasing (ea2t) functions.
A characteristic time dependence of the increasing x+ and decreasing x- terms in the solution for the function x(t) (the first coordinate of the vector z) at different coefficients j0 is presented in Fig. 3. It can be seen that under these initial conditions the decreasing term (x-) gives a contribution to the general solution only at the initial stage, for the time that does not exceed the period of one complete oscillation of the coordinate x(t). We show that the solution (2.7) is analogous in the general case as well.
To this end we calculate how many times N the amplitude of oscillations ea2t of the term z-in one period T = 2n// will decrease:
N = exp ^ - j . According to (2.5), the ratio in the exponent is
|a2| Vl + VI + e2 + 1
/
VVl +1~2 -T
> 1.
(2.9)
(2.10)
Therefore, for one period of oscillation the amplitude of the term z_ decreases N > exp(2^) & 535 times. Thus, as shown in the example in Fig. 3, the term z_ decreases to insignificant values in a time interval less than one period of oscillations.
Thus, under given initial conditions the solution of the system (1.7) is defined, except for a short initial period, by the increasing terms z+. This corresponds to the motion of the center of mass of the ball in a monotonically untwisting spiral around the projection of point zo onto the plane Oxy. Examples of trajectories of the center of mass of the ball for different coefficients of friction are given in Fig. 4.
The motion described above is similar to that of the ball on the plane if one incorporates the moment of resistance to rolling which is proportional to the angular velocity of the ball [2]. But in this case the untwisting of the spiral occurs much faster (see also [11]).
-20
x-
0.2
0.1
\
\ mo = 0.1
/¿o = l
n
\
4 8 t
Fig. 3. A characteristic dependence of the increasing (x+) and decreasing (x_) terms of the solution for the function x(t), which are constructed for = 0.1 and = 1, k = 2/5 and the initial values of the coordinates and velocities x(0) = 1, y(0) = 1, ^(0) = 0, vy(0) = 0, ^(0) = 0, (0) = 0.
Fig. 4. Trajectories of the center of mass of a ball moving with slipping on a rotating plane, under the initial conditions x(0) = 1, y(0) = 1, vx (0) = 0, vy(0) = 0, (0) = 0, wy(0) = 1, for different coefficients of friction. The heavy dot denotes the position of the point (x0,y0) (2.8), and z(0) is the projection of the vector z onto the plane Oxy at the initial time instant.
On the other hand, as the coefficient of friction j0 increases, the velocity of untwisting of the spiral decreases appreciably (see Fig. 4). That is, one can assume that in the limiting case j0 ^^ the trajectory will correspond to the nonholonomic rolling model in which the ball moves in closed circles with a constant frequency determined only by the radius of inertia of the ball [2]. We consider this limiting case in more detail in the next section.
3. The limiting case of the friction coefficient
In this section we analyze the solution of the system (1.7) for j0 — m. We note that, according to (2.2), e — 0 as j0 — <x>.
Let us expand the expression in brackets (2.4) in a power series of the small parameter e and represent the roots A of the characteristic equation in the form
Substituting (2.2) into (3.1), we obtain
ik k3 ^ / 1
Al'2 ±k + l + №(fc + l)3
Jo(fc +1) , Jk_ fc3 ,n(l_
k fc + 1 jo(,fc + l)3 + Uo
(3.1)
(3.2)
Thus, for large j0 the general solution (2.7) can now be represented as
+ exp (m^w) cos (ibTT') +1)3sin (kTl*.
(3.3)
where D1, D2, D3 are the real constant vectors defined in terms of the eigenvectors of the matrix L and from the initial conditions z(0).
Thus, on each fixed time interval, as fi0 the solution (3.3) tends to the function
z(t) = z0 + D2 cos + Di sin
which describes exactly the rolling of a homogeneous ball in the case of a nonholonomic dissipation-free model of motion [1, 2, 10].
Thus, the nonholonomic constraint that arises in the case of a ball rolling on a rotating plane can be realized by means of viscous friction forces and corresponds to the case where the coefficient of viscous friction is equal to infinity.
The author extends her gratitude to A. V. Borisov, I. S.Mamaev and A. A. Kilin for the formulation of the problem and fruitful discussions of the results obtained.
References
[1] Earnshaw, S., Dynamics, or An Elementary Treatise on Motion, 3rd ed., Cambridge: Deighton, 1844.
[2] Borisov, A. V., Ivanova, T.B., Karavaev, Yu. L., and Mamaev, I. S., Theoretical and Experimental Investigations of the Rolling of a Ball on a Rotating Plane (Turntable), Eur. J. Phys, 2018, vol. 39, no. 6, 065001, 13 pp.
[3] Fufaev, N. A., A Sphere Rolling on a Horizontal Rotating Plane, J. Appl. Math. Mech., 1983, vol.47, no. 1, pp. 27-29; see also: Prikl. Mat. Mekh, 1983, vol.47, no. 1, pp.43-47.
[4] Fufaev, N.A., On the Possibility of Realizing a Nonholonomic Constraint by Means of Viscous Friction Forces, J. Appl. Math. Mech., 1964, vol. 28, no. 3, pp. 630-632; see also: Prtkl. Mat. Mekh, 1964, vol.28, no. 3, pp. 513-515.
[5] Kolesnikov, S. L., The Problem of a Homogeneous Heavy Ball Rolling with Slipping in a Vertical Cylinder, in Collection of scientific and Methodical Papers on Theoretical Mechanics: Vol. 17, Moscow: Vysshaya Shkola, 1986, pp. 118-121 (Russian).
[6] Karapetyan, A. V. and Rumyantsev, V. V., Stability of Conservative and Dissipative Systems, Itogi Nauki Tekh. Ser. Obshch. Mekh., vol.6, Moscow: VINITI, 1983 (Russian).
[7] Eldering, J., Realizing Nonholonomic Dynamics as Limit of Friction Forces, Regul. Chaotic Dyn., 2016, vol.21, no. 4, pp. 390-409.
[8] Borisov, A. V. and Mamaev, I. S., Conservation Laws, Hierarchy of Dynamics and Explicit Integration of Nonholonomic Systems, Regul. Chaotic Dyn., 2008, vol. 13, no. 5, pp. 443-490.
[9] Borisov, A.V., Kilin, A.A., and Mamaev, I.S., The Hamilton Principle and the Rolling Motion of a Symmetric Ball, Dokl. Phys., 2017, vol. 62, no. 6, pp. 314-317; see also: Dokl. Akad. Nauk, 2017, vol.474, no. 5, pp. 558-562.
[10] Gersten, J., Soodak, H., and Tiersten, M.S., Ball Moving on Stationary or Rotating Horizontal Surface, Am. J. Phys., 1992, vol.60, no. 1, pp. 43-47.
[11] Kyeong Min Kim, Donggeon Oh, Junghwan Lee, Young-Gui Yoon, Chan-Oung Park. Dynamics of Cylindrical and Spherical Objects on a Turntable HAL Id: hal — 01761333, version 1, https://hal.archives-ouvertes.fr/hal-01761333 (2018).
