Inspired by Curling

This problem is inspired by the motion of the curling stone on the ice. Since we make several simplifying assumptions our result does not apply directly to curling,

Consider a ring-shaped object sliding on a flat horizontal surface. The motion of the ring can be decomposed to a translational motion of the center with velocity u u and to a rotation around the center with angular velocity ω \omega . The ring is flat on the surface and the pressure between the ring and the surface is equally distributed along the ring. The local friction force can be described as d F s = μ m g 2 π r d dF_s=\frac{\mu mg}{2\pi r}d\ell , where μ \mu is the friction coefficient, m m is the mass of the ring, g g is the gravitational acceleration, r r is the radius of the ring and d d\ell is a short segment of the ring. The direction of the friction force is opposite to the velocity of that particular segment.

First we start the ring without rotation with velocity u u , and observe it to slide to a distance s s . Next we start it with the same initial velocity, but we also give it a spin, so that ω = u / r \omega =u/r . The disk slides to a distance s s' .

What is the ratio s / s s'/s ?

Note: One can prove that if ω = u / r \omega =u/r initially, than ω = u / r \omega =u/r for the entire course of motion.


The answer is 1.57.

This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try refreshing the page, (b) enabling javascript if it is disabled on your browser and, finally, (c) loading the non-javascript version of this page . We're sorry about the hassle.

1 solution

Laszlo Mihaly
Mar 8, 2018

When the ring is pushed without rotation the equation of motion is simply u ˙ = μ g \dot u =-\mu g . The motion is a simple deceleration from the initial velocity to zero velocity.

For the rotating ring let us consider a small segment of the ring, of length d = r d φ d\ell=r d\varphi at an angle φ \varphi from the x x axis, as indicated in the Figure.

First we determine the direction of the friction force acting on this segment. The x x and y y components of the velocity are

v x = u ω r sin φ v_x=u-\omega r \sin \varphi

v y = ω r cos φ v_y=\omega r \cos \varphi

and the components of the friction force are

d F x = μ m g 2 π v x v d φ = μ m g 2 π u ω r sin φ v d φ d F_x= - \frac{\mu mg}{2\pi} \frac{v_x}{|v|} d\varphi= - \frac{\mu mg }{2\pi} \frac{u-\omega r \sin \varphi}{|v|} d\varphi

d F y = μ m g 2 π v y v d φ = μ m g 2 π ω r cos φ v d φ d F_y= - \frac{\mu mg}{2\pi} \frac{v_y}{|v|} d\varphi= - \frac{\mu mg }{2\pi} \frac{\omega r \cos \varphi}{|v|} d\varphi

where the absolute value of the velocity is

v = v x 2 + v y 2 = ( u ω r sin φ ) 2 + ω r cos 2 φ = u 2 + ( ω r ) 2 2 u ω r sin φ |v|=\sqrt{v_x^2+v_y^2}=\sqrt{(u-\omega r \sin \varphi)^2+ \omega r \cos^2 \varphi}=\sqrt{u^2+(\omega r)^2-2u\omega r \sin \varphi}

The total force is obtained by integration

F x = μ m g 2 π 0 2 π u ω r sin φ u 2 + ( ω r ) 2 2 u ω r sin φ d φ F_x= - \frac{\mu mg}{2\pi} \int_{0}^{2\pi}{ \frac{u-\omega r \sin \varphi}{\sqrt{u^2+(\omega r)^2-2u\omega r \sin \varphi}} d\varphi}

F y = μ m g 2 π 0 2 π ω r cos φ u 2 + ( ω r ) 2 2 u ω r sin φ d φ F_y= - \frac{\mu mg}{2\pi} \int_{0}^{2\pi}{ \frac{\omega r \cos \varphi}{\sqrt{u^2+(\omega r)^2-2u\omega r \sin \varphi}} d\varphi}

One can divide the integral for F y F_y into four quadrants and show that each quadrant has a pair where the integrand is the same, except for an opposite in sign. Therefore F y = 0 F_y=0 and the center of mass will move along the x x axis. The equation of motion for the center of mass is m u ˙ = F x m\dot{u}=F_x or

u ˙ = μ g 2 π 0 2 π η sin φ 1 + η 2 2 η sin φ d φ \dot{u}= - \frac{\mu g}{2\pi} \int_{0}^{2\pi}{ \frac{\eta-\sin \varphi}{\sqrt{1+\eta^2-2\eta \sin \varphi}} d\varphi}

where we introduced the dimensionless quantity η = u ω r \eta=\frac{u}{\omega r} .

The torque around the center of mass is d τ = d F y r cos φ d F x r sin φ d \tau = dF_y r \cos \varphi -dF_xr\sin \varphi or

d τ = r μ m g 2 π ( cos φ v x v sin φ v y v ) d φ = r μ m g 2 π ω r c o s 2 φ u sin φ + ω r sin 2 φ v d φ = r μ m g 2 π ω r u sin φ v d φ d \tau= - \frac{ r \mu mg}{2\pi}\left( \cos \varphi \frac{v_x}{|v|}-\sin \varphi \frac{v_y}{|v|}\right) d\varphi = - \frac{r \mu mg }{2\pi} \frac{\omega r cos^2\varphi - u \sin \varphi +\omega r \sin^2\varphi }{|v|} d\varphi = - \frac{r \mu mg }{2\pi} \frac{\omega r - u \sin \varphi}{|v|} d\varphi

The total torque is

τ = r μ m g 2 π 0 2 π 1 η sin φ 1 + η 2 2 η sin φ d φ \tau= - \frac{r \mu mg }{2\pi} \int_0^{2 \pi} \frac{1-\eta \sin \varphi}{\sqrt{1+\eta^2-2\eta \sin \varphi}} d\varphi

The equation of motion for the rotation around the center is I ω ˙ = τ I\dot\omega= \tau , where I I is the moment of inertia, I = m r 2 I=mr^2 for a ring. We can express the equation of motion as

ω ˙ r = r 2 μ m g 2 π I 0 2 π 1 η sin φ 1 + η 2 2 η sin φ d φ = μ g 2 π 0 2 π 1 η sin φ 1 + η 2 2 η sin φ d φ \dot\omega r = - \frac{r^2 \mu mg }{2\pi I} \int_0^{2 \pi} \frac{1-\eta \sin \varphi}{\sqrt{1+\eta^2-2\eta \sin \varphi}} d\varphi = - \frac{\mu g }{2\pi } \int_0^{2 \pi} \frac{1-\eta \sin \varphi}{\sqrt{1+\eta^2-2\eta \sin \varphi}} d\varphi .

Notice that the integrand in the equation for the rotational motion can be easily transformed into the integrand in the equation for the translational motion using η = 1 / η \eta'=1/\eta , because

1 η sin φ 1 + η 2 2 η sin φ = 1 / η sin φ 1 + ( 1 / η ) 2 2 ( 1 / η ) sin φ \frac{1-\eta \sin \varphi}{\sqrt{1+\eta^2-2\eta \sin \varphi}}=\frac{1/\eta-\sin \varphi}{\sqrt{1+(1/\eta)^2-2(1/\eta) \sin \varphi}} .

Let us introduce the function

f ( η ) = 0 2 π 1 η sin φ 1 + η 2 2 η sin φ d φ f(\eta)= \int_0^{2 \pi} \frac{1-\eta \sin \varphi}{\sqrt{1+\eta^2-2\eta \sin \varphi}} d\varphi

This function can be expressed in terms of elliptic integrals of second kind. For η = 1 \eta=1 the integral can be done directly, and it yields f ( 1 ) = 2 / π = 0.636 f(1)=2/\pi=0.636 .

It is easy to see that the two equations of motions can be simplified to

u ˙ = μ g f ( η ) \dot u= -\mu g f(\eta)

ω ˙ r = μ g f ( 1 / η ) \dot \omega r = -\mu g f(1/\eta)

We started with the initial condition of η = 1 \eta=1 . It is clear from the equations of motion that this condition is stable. Although u u and ω \omega will change with time, but the rate of change is the same and the ratio η = u / ω r \eta=u/\omega r does not change. The motion is a simple deceleration with an "effective" gravitational acceleration of g = g f ( η = 1 ) = 0.636 g g'=g f(\eta=1)= 0.636 g . Therefore the distance traveled is 1 / 0.636 1/0.636 times the distance traveled in simple translational motion, s / s = 1 / 0.636 = 1.57 s'/s=1/0.636=1.57 .

Note #1: If ω > 0 \omega >0 and u > 0 u>0 but ω u / r \omega \ne u/r initially the motion is more complicated, but eventually it converges to ω = u / r \omega =u/r before the ring stops. The translational and the rotational motion always stops at the same time.

Note #2: For an arbitrary object with axial symmetry we have I = α m r 2 I=\alpha mr^2 . A curling stone, for example, touches the ice like a ring, but its moment of inertia is different from the moment of inertia of a ring. If we use this more general formula for the moment of inertia, the equations of motion are

u ˙ = μ g f ( η ) \dot u= -\mu g f(\eta)

ω ˙ r = μ g α f ( 1 / η ) \dot \omega r = -\frac{\mu g}{\alpha} f(1/\eta)

One may ask, is it possible to have a motion with constant η \eta ? To answer that question, let us look at the time derivative of η \eta and make it equal to zero:

η ˙ = d d t u ω r = u ˙ ω r u ω ˙ ω 2 r = 0 \dot \eta= \frac{d}{dt}\frac {u}{\omega r}= \frac{\dot u}{\omega r}-\frac{u\dot \omega}{\omega^2 r}=0

This yields the condition u ˙ ω ˙ r = u ω r \frac{\dot u}{\dot \omega r}=\frac{u}{\omega r} . From the equations of motions we get:

f ( η ) = η α f ( 1 / η ) f(\eta)=\frac{\eta}{\alpha}f(1/\eta)

First we note that for α = 1 \alpha=1 we indeed have a solution with η = 1 \eta=1 . This is what we had before. In the more general case we can solve the equation numerically, and we find that "stable" η \eta exists in the range of 0.5 < α < 2 0.5<\alpha<2 . At the lower end η \eta approaches infinity as α 1 / 2 \alpha \rightarrow 1/2 . For α 2 \alpha \rightarrow 2 , η 0 \eta \rightarrow 0

Note #3, about curling: The main difference between this problem and the curling stone is that the friction between the ice and the stone has a much more complicated behavior. Also, the curling stone has a height that is comparable to the radius of the circle where it touches the ice. As a result, when the stone decelerates there will be higher pressure in the front relative to the back side of the stone and the friction forces pointing sideways will not cancel each other. This will cause the stone to go sideways ( y y direction in our problem). There is a discussion of this in the New Yorker magazine and in Scientific American . Mark Shegelski, a curler and a physicist at the University of Northern British Columbia, published several papers on the subject matter.

We started with three points. Then moved to a continuous circular line. Would the next logical step be a circular disk, effectively moving from 1D to 2D?

Steven Chase - 3 years, 3 months ago

Log in to reply

Thanks for solving it and thanks for the comment. I think a circular disk is not much different from the ring. The next step could be to do this with a cylinder of finite height. In that case there will be a force in the y direction, creating a curling effect. BTW, with the usual friction law a clockwise rotation leads to a curling to the right, because the normal force and the friction is stronger at the front side of the motion. The curling stone on ice curls in the opposite direction, because the friction is actually less when the normal force is larger.

Laszlo Mihaly - 3 years, 3 months ago

Log in to reply

For the cylinder, would the contact area be a disk (the 2D projection of the cylinder)? If so, why does the height matter?

Steven Chase - 3 years, 3 months ago

Log in to reply

@Steven Chase I would take a thin-walled cylinder, with a contact area of a circle. As it is sliding and decelerates the friction force will create a torque around the center of mass. In the extreme case when the height is much larger than the diameter the cylinder would tip forward. If the height is less than that the torque created by friction can be balanced by a torque created by the normal force. Accordingly the normal force will be larger on the front side than on the back side.

Laszlo Mihaly - 3 years, 3 months ago

Log in to reply

@Laszlo Mihaly Interesting. So the load bearing profile would change over time. Sounds very non-trivial

Steven Chase - 3 years, 3 months ago

@Laszlo Mihaly I put one up which has a rod with the normal force distributed between the two end points.

Steven Chase - 3 years, 3 months ago

Log in to reply

@Steven Chase It looks very interesting! I think the sideways displacement (non-zero y coordinate) in your problem is coming from the fact that this is a rod with two support points, a very much asymmetric object. If the support is distributed along a ring or along the surface of circle the sideways motion is exactly 0.

Laszlo Mihaly - 3 years, 3 months ago

0 pending reports

×

Problem Loading...

Note Loading...

Set Loading...