Conservación del momento lineal y angular en las colisiones de dos discos I

En el capítulo de dinámica hemos estudiado las colisiones unidimensionales elásticas e inelásticas tanto desde el punto de vista de un observador situado en el Sistema de Referencia del Laboratorio como del Sistema de Referencia del Centro de Masa. A continuación, procedimos con el estudio de las colisiones bidimensionales. En ambos casos, hemos aplicado el principio de conservación del momento lineal a un sistema aislado de dos partículas interactuantes y a continuación, hemos efectuado el balance energético de la colisión.

En la figura, se muestra un esquema de la colisión de dos discos de masas m₁ y m₂, y radios r₁ y r₂ respectivamente. El segundo disco está en reposo u₂=0, mientras que el primero lleva una velocidad u₁ antes del choque. El choque está caracterizado por el denominado parámetro de impacto b, que es la distancia entre la dirección de la velocidad u₁ del primer disco y el centro del segundo disco en reposo.

Después del choque, el primer disco se mueve con velocidad v₁ haciendo un ángulo φ₁ con eje horizontal (dirección de la velocidad u₁ del disco incidente) y a su vez, gira alrededor de un eje perpendicular al disco y que pasa por su centro con velocidad angular w_1. El segundo disco, se mueve con velocidad v₂ haciendo un ángulo φ₂ con el eje horizontal y a su vez, gira alrededor de de un eje perpendicular al disco y que pasa por su centro con velocidad angular w₂.

Tenemos que despejar seis incógnitas del sistema de ecuaciones que describe el choque entre dos discos:

• Los módulos de las velocidades de los c.m. de los discos después del choque v₁ y v₂, y sus direcciones φ₁ y φ₂
• Las velocidades angulares de rotación de cada uno de los discos alrededor de un eje perpendicular al disco y que pasa por el c.m. w₁  y w₂.

Para obtener las seis ecuaciones, aplicamos:

• El principio de conservación del momento lineal al sistema aislado formado por los dos discos, nos proporcionan dos ecuaciones.
• El principio de conservación del momento angular  cada uno de los discos respecto del punto de contacto, nos proporcionan dos ecuaciones.
• El balance energético de la colisión a través del coeficiente de restitución, nos proporciona la quinta ecuación
• Cuando los discos están en contacto, puede deslizar una superficie sobre la otra o no deslizar. Tenemos una ecuación más, la sexta, distinta para cada caso.

Choques frontales

Los choques frontales son las más fáciles de describir ya que solamente, precisan la aplicación del principio de conservación del momento lineal y la definición de coeficiente de restitución.

m₁u₁ =m₁v₁+m₂v₂                (1)

 -e·u₁=v₁-v₂                                    (2)

Choques oblicuos

Parámetro de impacto

Se denomina parámetro de impacto b, a la distancia entre la dirección de la velocidad u₁ del primer disco y el centro del segundo disco en reposo. La relación entre le parámetro de impacto b y el ángulo θ que forma la dirección de la velocidad u₁ del primer disco y la recta que pasa por los centros de ambos discos, cuando entran en contacto en el momento del choque, se puede apreciar en la figura.

b=(r1+r2)·senθ

Conservación del momento lineal

• a lo largo del eje X

m₁u₁ =m₁v_1x+m₂v_2x                (1)

• a lo largo del eje Y

0=m₁v_1y+m₂v_2y                       (2)

Conservación del momento angular

Las fuerzas de interacción entre los discos se aplican en el punto de contacto P. El momento de dichas fuerzas respecto de P es cero.

• Se conserva el momento angular del disco 1 respecto del punto de contacto P.

• Se conserva el momento angular del disco 2 respecto del punto de contacto P.

-m₁·r₁·u₁senq= I₁ω₁- m₁·r₁·(v_1xsenθ+ v_1ycosθ)                   (3)

0=I₂ω₂ +m₂·r₂·(v_2xsenθ+ v_2ycosθ)                                    (4)

Balance energético. Coeficiente de restitución

La definición de coeficiente de restitución es

                 (5)

Cuando los discos están en contacto

• No hay deslizamiento

Consideremos el caso de que no hay deslizamiento de un disco respecto del otro en el punto de contacto P. Las velocidades de los dos discos en el punto de contacto P serán iguales.

r₁ω₁+v_1xsenθ+ v_1ycosθ= -r₂ω₂ + v_2xsenθ+ v_2ycosθ                     (6)

• Hay deslizamiento

Las fuerzas sobre el disco azul en el punto de contacto son:

• La reacción N,
• La fuerza de rozamiento F

 Si las superficie lateral del disco azul desliza sobre la del rojo, la relación entre ambas fuerzas es F=μ·N
La fuerza N actuando durante el pequeño intervalo de tiempo Dt en el que los discos están en contacto modifica la componente del momento lineal del disco en la dirección de dicha fuerza. Lo mismo cabe decir de la fuerza F.

μ(v_1ysenθ-v_1xcosθ+u₁cosθ)=- v_1xsenθ-v_1ycosθ+u₁senθ              (6)
Las fuerzas sobre el disco rojo en el punto de contacto son:

• La reacción N,
• La fuerza de rozamiento F

iguales y de sentido contrario a las que se ejercen sobre el disco azul.


μ(v_2xcosθ-v_2ysenθ)=v_2xsenθ+v_2ycosθ                                    (7)

Resolución del sistema de ecuaciones

Llamamos M=m₁/m₂

El momento de inercia de un sólido rígido en rotación alrededor de un eje fijo es

I=kmr² con k=1/2 para un disco

El sistema de seis ecuaciones con seis incógnitas se escribe.

Mv_1x+v_2x= Mu₁+u₂         (1)

Mv_1y+v_2y=0                      (2)

v_1xsenθ+ v_1ycosθ-kr₁ω₁= u₁senq            (3)
v_2xsenθ+ v_2ycosθ +kr₂ω₂ =0                     (4)

-v_1xcosθ + v_1ysenθ+v_2xcosθ-v_2ysenθ=eu₁cosθ           (5)

v_1xsenθ+ v_1ycosθ - v_2xsenθ-v_2ycosθ +r₁ω₁+r₂ω₂ =0      (6)   No desliza
(- μcosθ+senθ)v_1x+(μsenθ +cosθ) v_1y=u₁(senθ- μcosθ)   (6) Desliza

Para resolver el sistema de ecuaciones lineales, cambiamos el orden de las ecuaciones para evitar que los elementos de la diagonal de la matriz cuadrada (más abajo) sean nulos

• No desliza

v_1xsenθ+ v_1ycosθ - v_2xsenθ- v_2ycosθ +r₁ω₁+r₂ω₂ =0      (6)  
Mv_1y+v_2y=0                                                                       (2)
Mv_1x+v_2x= Mu₁+u₂                                                           (1)
-v_1xcosθ + v_1ysenθ+v_2xcosθ-v_2ysenθ=eu₁cosθ                  (5)
v_1xsenθ+ v_1ycosθ-kr₁ω₁= -u₁senq                                      (3)
v_2xsenθ+ v_2ycosθ +kr₂ω₂ =0                                              (4)
En forma matricial escribimos el sistema

Se resuelve el sistema de ecuaciones se evalúa el primer y segundo miembro de la ecuación (6) caso "desliza"

A=μ(v_1ysenθ-v_1xcosθ+u₁cosθ)
B=- v_1xsenθ-v_1ycosθ+u₁senθ

Si A>B se mantiene el resultado, en caso contrario se resuelve el sistema alternativo

• Desliza

(- μcosθ+senθ)v_1x+(μsenθ +cosθ) v_1y=u₁(senθ- μcosθ)    (6)
Mv_1y+v_2y=0                                                                       (2)
Mv_1x+v_2x= Mu₁+u₂                                                           (1)
-v_1xcosθ + v_1ysenθ+v_2xcosθ-v_2ysenθ=eu₁cosθ                  (5)
v_1xsenθ+ v_1ycosθ-kr₁ω₁= -u₁senq                                      (3)
v_2xsenθ+ v_2ycosθ +kr₂ω₂ =0                                              (4)

En forma matricial escribimos el sistema

Comprobamos que la solución es correcta en este caso, verificándose la ecuación

μ(v_2xcosθ-v_2ysenθ)=v_2xsenθ+v_2ycosθ                                    (7)

Procedimiento de cálculo

Llamamos M a la matriz, X al vector columna de las incógnitas y B al vector columna de los términos independientes

M·X=B

Despejamos el vector columna de las incógnitas

X=M^-1·B

Para hallar los elementos del vector columna X de las incógnitas, calculamos la matriz inversa M^-1 y a continuación, la multiplicamos por el vector columna B de los términos independientes.

Ejemplo:

Comprobaremos,  la conservación del momento lineal y del momento angular de los discos, a partir de los datos suministrados por el programa interactivo.

Datos relativos a los discos

• Las masas de los discos, m₁=1, m₂=1, por lo que M=m₁/m₂=1
• Los radios de los discos  r₁=1,  r₂=1
• Discos ambos de acero, e=0.94 y m =0.10

Antes del choque

• Velocidad inicial u₁=3.5 del primer disco
• Parámetro de impacto b=1.5, por lo que senθ=1.5/(1+1), θ=48.60º

Después del choque

• Velocidad del primer disco v₁=2.40, dirección φ₁=39.75º (por encima de la horizontal)
• Velocidad del segundo disco v₂=2.26, dirección φ₂=-42.88º (por debajo de la horizontal)
• Velocidad angular del primer disco r₁·ω₁=-0.45, o bien w₁=-0.45 (gira en el sentido de las agujas del reloj)
• Velocidad angular del segundo disco r₂·ω₂=-0.45, o bien w₂=-0.45

1.- Conservación del momento lineal

El momento lineal inicial del primer disco (el segundo está inicialmente en reposo) es igual a la suma vectorial de los momentos lineales de los discos después del choque.m1u1=m1v1cosφ1+m2v2cosφ2 0=m1v1senφ1+m2v2senφ2

• Antes del choque

Eje X: 1·3.5=3.5
Eje Y: 0.0

• Después del choque

Eje X: 1·2.40·cos39.75+1·2.26·cos(-42.88)=3.49
Eje Y: 1·2.40·sen39.75+2·2.26·sen(-42.88)=0.003 » 0.0

El momento lineal se conserva
2.-Conservación del momento angular
Se conserva el momento angular del disco 1 respecto del punto de contacto P.

-m₁·r₁·u₁senq= I₁ω₁- m₁·r₁·v₁sen(θ+φ₁)
-1·1·3.5·sen48.60=1·1²·(-0.45)/2-1·1·2.40·sen(48.60+39.75)
(Las velocidades angulares son positivas en el sentido contrario al de las agujas del reloj)
Se conserva el momento angular del disco 2 respecto del punto de contacto P.

0=I₂ω₂ +m₂·r₂·v₂sen(θ+φ₁)
0=1·1²·(-0.45)/2-1·1·2.26·sen(48.60+-42.88)

3. Balance energético

Energía inicial
E_i=6.125
Energía final
E_f=5.535

 Energía perdida en la colisión Q=E_f-E_i=-0.590

Ejemplo 2º Choques frontales

Datos relativos a los discos

• Masas de los discos m₂=1, m₁=0.5, por lo que M=m₁/m₂=0.5
• Discos ambos de acero, e=0.94

Antes del choque

• Velocidad inicial u₁=3.5 del primer disco
• Parámetro de impacto b=0.0, por lo que θ=0º

Después del choque

• Velocidad del primer disco v₁=-1.03
• Velocidad del segundo disco v₂=2.26,

1. Conservación del momento lineal

• Antes del choque

Eje X: 1.3.5=3.5

• Después del choque

Eje X: -1·1.03+2·2.26=3.49

Se conserva el momento lineal

2. Balance energético

Antes del choque

Después del choque

 Energía perdida en la colisión Q=E_f-E_i=-0.243