Trayectoria de la bola en la bolera americana

Sea una bola de masa m, radio R y momento de inercia, Ic=2mR2/5. En un instante dado t, el punto de contacto P de la esfera con la pista se encuentra en la posición (x, y). La velocidad del centro de masas es V , y la velocidad angular de rotación ω alrededor de un eje que pasa por el centro de masas.

La velocidad del punto P de contacto de la bola con la pista es

v p = V + ω × R v p =( V x V y 0 )( i ^ j ^ k ^ ω x ω y ω z 0 0 R )=( V x R ω y ) i ^ +( V y +R ω x ) j ^ { v x = V x R ω y v y = V y +R ω x

Ecuaciones del movimiento

Las fuerzas sobre la bola son:

Sean Fx y Fy las componentes de esta fuerza

Las ecuaciones del movimiento son

Relacionamos las componentes de la velocidad del centro de masas con las componentes de la velocidad angular de rotación. Teniendo en cuenta que el momento de inercia Ic=2mR2/5.

R d ω x dt = 5 2 d V y dt R d ω y dt = 5 2 d V x dt

Hemos calculado las componentes de la velocidad del punto P de contacto con la pista, derivando con respecto del tiempo

{ d v x dt = d V x dt R d ω y dt = 7 2 d V x dt = 7 2 F x m d v y dt = d V y dt +R d ω x dt = 7 2 d V y dt = 7 2 F y m

Cuando el punto P desliza sobre la pista, la fuerza de rozamiento es de signo contrario a la velocidad del punto P y vale, F r =μN v ^ =μmg v ^ , siendo μ el coeficiente cinético y v ^ el vector unitario en la dirección de la velocidad del punto de contacto P. Las ecuaciones son

{ d v x dt = 7 2 μg v x v x 2 + v y 2 d v y dt = 7 2 μg v y v x 2 + v y 2

La ecuación de una línea recta de pendiente m que pasa por el origen es y=mx, dy=m·dx, o bien dy/y=dx/x.

Las componentes de la velocidad del punto de contacto P están relacionadas, vx=cvy, donde la constante c se determina a partir de las condiciones iniciales , en el intente t=0, la velocidad del centro de masas es (V0x, V0y) y la velocidad angular de rotación (ω0x, ω0y)

c= v 0x v 0y = V 0x R ω 0y V 0y +R ω 0x

Integramos la primera ecuación diferencial, sustituyendon vx=c·vy

d v x dt = 7 2 μg c v y c 2 v y 2 + v y 2 v x = v 0x 7 2 μg c c 2 +1 t

Integramos la segunda ecuación diferencial, de forma similar

d v y dt = 7 2 μg v y c 2 v y 2 + v y 2 d v y dt = 7 2 μg 1 c 2 +1 v y = v 0y 7 2 μg 1 c 2 +1 t

La velocidad del punto P de contacto se anula cuando vx y vy se anulen simultáneamente, en el instante t tal que

v 0x = 7 2 μg c c 2 +1 t v 0y = 7 2 μg 1 c 2 +1 t

Elevando al cuadrado y sumando

v 0x 2 + v 0y 2 = 49 4 μ 2 g 2 c 2 +1 c 2 +1 t 2 t c = 2 7 v 0 μg

A partir de este instante, la bola rueda sin deslizar

Movimiento del centro de masas

La recta Y' es la dirección constante de la velocidad v del punto P y por tanto, de la fuerza de rozamiento. forma un ángulo c=tanθ con el eje Y

Establecemos un sistema de coordenadas X',Y' giradas un ángulo θ respecto de X,Y , tal como se muestra en la figura. El movimiento del centro de masas de la bola es composición de dos movimientos: uniforme a lo largo del eje X' y uniformemente acelerado a lo largo del eje Y'

{ a x ' =0 a y ' =μg { V x ' = V 0x ' V y ' = V 0y ' μgt { x'= V 0x ' t y'= V 0y ' t 1 2 μg t 2

La posición y velocidad del centro de masas en el instante t=tc es

{ V cx ' = V 0x ' V cy ' = V 0y ' μg t c { x c ' = V 0x ' t c y c ' = V 0y ' t c 1 2 μg t c 2

A partir de este instante, la bola rueda sin deslizar y el centro de masas se mueve a lo largo de una línea recta con velocidad constante

{ x'= x c ' + V cx ' (t t c ) y'= y c ' + V cy ' (t t c )

Transformación de coordenadas

La relación entre la posición (x, y) del centro de masas C de la bola en el Sistema de Referencia X,Y y la misma posición (x', y') en el Sistema de Referencia X',Y' se puede deducir de la figura

{ x=x'cosθ+y'sinθ y=x'sinθ+y'cosθ

Teniendo en cuenta las relaciones trigonométricas

cosθ= 1 1+ tan 2 θ ,sinθ= tanθ 1+ tan 2 θ

Obtenemos las relaciones

{ x= 1 c 2 +1 ( x'+cy' ) y= 1 c 2 +1 ( cx'+y' )

Trayectoria

En el instante t=0, la velocidad del centro de masas es (V0x, V0y), y la velocidad angular inicial de rotación alrededor de un eje que pasa por el centro de masas es (ω0x, ω0y).

Para dibujar la trayectoria de la bola en la pista seguimos estos pasos

Ejemplo 1

function bolera
    V0x=0; %velocidades iniciales
    V0y=18*30.48/100;
    w0x=0;
    w0y=4;
    R=0.354*30.48/100; %radio
    mu=0.09; %coeficiente de rozamiento
    v0x=V0x-R*w0y;
    v0y=V0y+R*w0x;
    c=v0x/v0y; %tangente de th
    tc=2*sqrt(v0x^2+v0y^2)/(7*9.8*mu); %tiempo deslizando

%Desliza
    Vp0x=(V0x-c*V0y)/sqrt(c^2+1);
    Vp0y=(c*V0x+V0y)/sqrt(c^2+1);
    t=linspace(0,tc,50);
    xp=Vp0x*t;
    yp=Vp0y*t-mu*9.8*t.^2/2;
    [x,y]=gira(xp,yp);
    hold on
    plot(x,y)
%Rueda sin deslizar
    Vcx=Vp0x;
    Vcy=Vp0y-mu*9.8*tc;
    xc=Vp0x*tc;
    yc=Vp0y*tc-mu*9.8*tc^2/2;
    t=linspace(tc,2*tc, 50);
    xp=xc+Vcx*(t-tc);
    yp=yc+Vcy*(t-tc);
    [x,y]=gira(xp,yp);
    plot(x,y)
    [x,y]=gira(xc,yc);
    line([x,x],[0,y],'lineStyle','--')
    line([0,x],[y,y],'lineStyle','--')
    hold off
    grid on
    xlabel('x')
    ylabel('y')
    title('Movimiento de la bola')
    
    function [x,y]=gira(xp,yp)
       x=(xp+c*yp)/sqrt(c^2+1); 
       y=(-c*xp+yp)/sqrt(c^2+1); 
    end 
end

En la figura, se distingue las dos trayectorias que sigue la bola: cuando el punto P de contacto de la bola con la pista desliza y cuando la bola rueda sin deslizar. Se señala la posición de la bola (xc=0.1099 , yc=8.3836) m en el instante tc=1.7828 s

Ejemplo 2

Cambiamos las velocidades iniciales

En la figura, se distingue las dos trayectorias que sigue la bola: cuando el punto P de contacto de la bola con la pista desliza y cuando la bola rueda sin deslizar. Se señala la posición de la bola (xc=0.3690, yc=8.9070) m en el instante tc=1.9192 s

Para que se pueda apreciar la trayectoria, las escalas en el eje X e Y son distintas

Referencias

D. C. Hopkins, J. D. Patterson. Bowling frames: Paths of a bowling ball. Am. J. Phys. 45 (3) March 1977, pp. 263-266