On the Motion of the Oloid Toy
Total Page:16
File Type:pdf, Size:1020Kb
On the motion of the Oloid toy On the motion of the Oloid toy Alexander S. Kuleshov Mont Hubbard Dale L. Peterson Gilbert Gede [email protected] Abstract Analysis and simulation are performed for the commercially available toy known as the Oloid. While rolling on the fixed horizontal plane the Oloid moves very swinging but smooth: it never falls over its edges. The trajectories of points of contact of the Oloid with the supporting plane are found analytically. 1 Introduction Let us consider the motion of the Oloid on a fixed horizontal plane. The Oloid is a developable surface comprises of two circles of radius R whose planes of symmetry make a right angle between each other with the distance between the centers of the circles equals to their radius R. The resulting convex hull is called Oloid. The Oloid have been constructed for the first time by Paul Schatz [2,3]. The geometric properties of the surface of the Oloid have been discussed in the paper [1]. The Oloid is also used for technical applications. Special mixing-machines are constructed using such bodies [4]. In our paper we make the complete kinematical analysis of motion of this object on the horizontal plane. Further we briefly describe basic facts from Kinematics and Differential Geometry which we will use in our investigation. The Frenet - Serret formulas. Consider a particle which moves along a continuous differentiable curve in three - dimensional Euclidean Space 3. We can introduce the R following coordinate system: the origin of this system is in the moving particle, τ is the unit vector tangent to the curve, pointing in the direction of motion, ν is the derivative of τ with respect to the arc-length parameter of the curve, divided by its length and β is the cross product of τ and ν: β = [τ ν]. × Then the Frenet - Serret formulas for the derivatives of τ , ν and β are valid: dτ /dt = ks˙ν, dν/dt = ks˙τ + æs ˙β, dβ/dt = æs ˙ν. − − The tangent τ , the normal ν and the binormal β unit vectors are known as the Frenet - Serret frame. The Poisson formulas. Let e1, e2, e3 is any moving coordinate system and let ω be the angular velocity of this system. Then for the derivatives of e1, e2 and e3 the following Poisson formulas are valid: de /dt = [ω e ] , i = 1, 2, 3. i × i According to the Poisson formulas we can rewrite the Frenet - Serret formulas as follows: dτ /dt = [Ω τ ] , dν/dt = [Ω ν] , dβ/dt = [Ω β] , × × × 275 Proceedings of XXXIX International Summer School–Conference APM 2011 Figure 1: Kinematics of the moving lamina. where Ω = ks˙β + æs ˙τ is the angular velocity of the Frenet - Serret frame known also as the Darboux vector. In particular, for the plane curve we have æ = 0 and therefore Ω = ks˙β. (1) The law of composition of angular velocities. Consider three coordinate systems: the fixed system Oxyz and two moving systems Aξηζ and Cx1x2x3. Then the angular velocity of the system Cx1x2x3 with respect to Oxyz - system is the sum of the angular velocity of the Cx1x2x3 - system with respect to the Aξηζ - system and the angular velocity of the Aξηζ - system with respect to the Oxyz - system: ωOxyz = ωOxyz + ωAξηζ . Cx1x2x3 Aξηζ Cx1x2x3 We will call the angular velocity ωAξηζ of the moving coordinate system Cx x x Cx1x2x3 1 2 3 with respect to another moving system Aξηζ as ”the relative angular velocity”. The Oxyz angular velocity ωAξηζ will called ”the transfer angular velocity” and the angular velocity ωOxyz will called ”the absolute angular velocity”. Cx1x2x3 Angular velocity of the moving lamina. Let us consider the following problem: the rigid lamina bounded by the convex curve of the given curvature k(s) is rolling without slipping on the fixed convex curve of the given curvature K(s) (we assume that the both curves have the same arc parametrization s = s(t)). Let us find the angular velocity Ω of the moving lamina. For this purpose we introduce three coordinate system: the fixed system Oxyz, con- nected with the fixed curve, the moving system Aξηζ rigidly connected with the moving lamina and the Frenet - Serret frame Cτνβ with the origin at the point of contact of two curves (Fig. 1). We will use the law of composition of angular velocities. The absolute angular velocity of the Frenet - Serret frame with respect to the fixed coordinate system Oxyz is Ks˙β (according to (1)). But this angular velocity is the vector sum of relative and transfer angular velocities. The relative angular velocity of the Frenet - Serret frame is its angular velocity with respect to the Aξηζ - system. This angular velocity is equal ks˙β (since in the moving frame Aξηζ the motion of the Frenet - Serret frame Cτνβ is the motion along the bound of the lamina). The transfer angular velocity of the Frenet - Serret frame is the angular velocity of the Aξηζ - system with respect to Oxyz - system, i.e. it is the angular velocity of the moving lamina Ω. Finally we have Ω = (K k)s ˙β. − In particular when the lamina moves along the straight line, its angular velocity equals Ω = ks˙β. − 276 On the motion of the Oloid toy Two contact points. Let us consider the rolling of the body consisting of two similar symmetric laminas in perpendicular planes connected to each other along the common axis of symmetry.This body rolls without slipping on the fixed horizontal plane. Let G be the center of mass of the moving body and K1, K2 be two contact points of the body with the plane. From the rolling conditions for the two contact points follows that the angular velocity is always parallel to the line between the contact points i.e. ω −−−→K K = 0. × 1 2 h i Now we introduce four coordinate systems. The first coordinate system is the fixed coordinate system Oxyz with the origin at any point O of the fixed plane. The Oz - axis is directed upward. The second coordinate system is Gx1x2x3 system with Gx2 axis directed along the common axis of symmetry of the two laminas. The Gx3 - axis is perpendicular to the plane of the first lamina and the Gx1 - axis is perpendicular to the plane of the second lamina. The third coordinate system is the Frenet - Serret frame K1τνβ for the bound of the first lamina: the τ - vector is the tangent vector to the bound of the lamina, the ν - vector is in the plane of the lamina and the β - vector coincides in the unit vector of Gx3 axis (i.e. β is perpendicular to the plane of lamina). The fourth coordinate system is also the Frenet - Serret frame. The origin of this system is also in K1, the first unit vector τ 1 of this system coincides with the vector τ of the system K1τνβ. The second vector ν1 of this system is in the plane of motion and the third vector β1 is perpendicular to the plane of motion (i.e. the vector β1 coincides with the vector ez of the system Oxyz). Therefore we can conclude that the coordinate system K1τ 1ν1beta1 is the Frenet - Serret frame for the curve which is obtained by the motion of the contact point K1 (i.e. it is a ”trace” of the point of contact on the plane). Let us try to find the angular velocity of the moving body. We will use the law of composition of angular velocities. Let us consider the systems Oxyz, K1τνβ and Gx1x2x3. The absolute angular velocity of Gx1x2x3 system with respect to Oxuz system is the sum of angular velocity of Gx1x2x3 with respect to K1τνβ and the angular velocity of K1τνβ with respect to Oxyz. But the angular velocity of Gx1x2x3 system with respect to K1τνβ system can be found very easy: it is equal ks˙β. So we need calculate now the angular − velocity of K1τνβ system with respect to Oxyz system. Let us denote by ϕ the angle between two unit vectors β = e3 and β1 = ez. Then we can represent the angular velocity of K1τνβ system with respect to Oxyz system is the sum of the angular velocity of K1τνβ system with respect to K1τ 1ν1β1 system and the angular velocity of K1τ 1ν1β1 system with respect to Oxyz system. But the angular velocity of K1τνβ with respect to K1τ 1ν1β1 system is equalϕ ˙τ and the angular velocity of K1τ 1ν1β1 system with respect to Oxyz system is equal Ks˙β1. Therefore finally we obtain: ω =ϕ ˙τ ks˙β + Ks˙β . − 1 Since β = ν sin ϕ + β cos ϕ, we have − 1 1 ω =ϕ ˙τ + ks˙ sin ϕν ks˙ cos ϕβ + Ks˙β . 1 1 − 1 1 But from the rolling conditions we already obtained that the vector ω is proportional to −−−→K1K2 vector i.e. ω is always in the plane of motion ω = ω1τ 1 + ω2ν1. 277 Proceedings of XXXIX International Summer School–Conference APM 2011 This means that K = k cos ϕ or ρ cos ϕ = r, (2) where ρ and r are radii of curvature of the curve on the fixed plane (trajectory of the point K1) and the bound of the first lamina. Natural equations of a curve. Let us consider the planar curve, defined by its parametric equations: r = r (s)= x(s)ex + y(s)ey, with the arc-length s as a parameter. We denote by α(s) the angle between the unit tangent vector τ to the given curve dr dx dy τ = = e + e ds ds x ds y and the unit vector ex of the Ox - axis.