]> Movimiento de un pistón

Movimiento de un pistón

El sistema biela-manivela de una máquina motriz (máquina de vapor, motor térmico) se compone de una biela AB cuyo extremo A llamado pie de biela, se desplaza a lo largo de una recta, mientras que el otro extremo B, llamado cabeza de biela, articulado en B con una manivela OB describe una circunferencia de radio OB. El pie de biela está articulado en una pieza denominada patín solidaria con el pistón que se desplaza entre dos guías. El pistón describe un movimiento oscilatorio que como vamos a ver no es armónico simple, aunque se puede aproximar bastante a éste.

Descripción del movimiento

Supongamos que la manivela tiene radio r, y la biela tiene una longitud l (l>2r). La manivela gira con velocidad angular constante ω, y el pistón oscila. La posición del pistón respecto del centro de la rueda es

x e =r·cosθ+ l 2 r 2 sin 2 θ +d

Si situamos el origen en la posición en la posición del pistón para θ=90º.

x O = l 2 r 2 +d

Posición del pistón

x= x e x O =r·cosθ+ l 2 r 2 sin 2 θ l 2 r 2

Si la manivela se mueve con velocidad angular ω constante, la posición del pistón en función del tiempo es

x=r·cos(ωt)+ l 2 r 2 sin 2 (ωt) l 2 r 2

x=r+l l 2 r 2

x=r+l l 2 r 2

En la figura, se representa la posición x del pistón en función del tiempo (color azul) y el MAS (color rojo)

x=r·sin(ω·t+π/2)=r·cos(ω·t)

Velocidad

Derivando la posición x con respecto al tiempo obtenemos la velocidad

v= dx dt =rωsin(ωt)( 1+ rcos(ωt) l 2 r 2 sin 2 (ωt) )

En la figura, se representa la velocidad v del pistón en función del tiempo (color azul) y el MAS (color rojo)

v=-r·ω·sin(ω·t)

Aceleración

Derivando la velocidad v con respecto al tiempo obtenemos la aceleración

a= dv dt =r ω 2 cos(ωt)( 1+ rcos(ωt) l 2 r 2 sin 2 (ωt) ) rωsin(ωt)( rωsin(ωt) l 2 r 2 sin 2 (ωt) +rcos(ωt) r 2 ωcos(ωt)sin(ωt) l 2 r 2 sin 2 (ωt) l 2 r 2 sin 2 (ωt) )= r ω 2 cos(ωt)+ r 2 ω 2 cos 2 (ωt)( l 2 r 2 sin 2 (ωt))+ r 2 ω 2 sin 2 (ωt)( l 2 r 2 sin 2 (ωt)) r 4 ω 2 cos 2 (ωt) sin 2 (ωt) ( l 2 r 2 sin 2 (ωt)) 3/2

Simplificando se llega al resultado

a= dv dt =r ω 2 ( cos(ωt)+ r( l 2 cos(2ωt)+ r 2 sin 4 (ωt)) ( l 2 r 2 sin 2 (ωt)) 3/2 )

En la figura, se representa la aceleración a del pistón en función del tiempo (color azul) y el MAS (color rojo)

a=-r·ω2·cos(ω·t)

Aceleración nula, máxima velocidad

La velocidad es máxima cuando la aceleración es cero. Para calcular los ángulos ω·t para los cuales la aceleración es cero, hay que resolver la ecuación

cos(ωt) ( l 2 r 2 sin 2 (ωt)) 3/2 =r( l 2 cos(2ωt)+ r 2 sin 4 (ωt))

Las operaciones que hay que realizar son las siguientes

Se sustituye sin2(ωt)=1-cos2(ωt) y cos(2ωt)= cos2(ωt) - sin2(ωt) =2 cos2(ωt)-1

cos(ωt) ( l 2 r 2 + r 2 cos 2 (ωt) ) 3/2 =r( 2( l 2 r 2 ) cos 2 (ωt)( l 2 r 2 )+ r 2 cos 4 (ωt) )

Se eleva al cuadrado ambos miembros

cos 2 (ωt) ( l 2 r 2 + r 2 cos 2 (ωt) ) 3 = r 2 ( 2( l 2 r 2 ) cos 2 (ωt)( l 2 r 2 )+ r 2 cos 4 (ωt) ) 2

y después de realizar algunas operaciones algebraicas se llega a la ecuación

r 4 cos 6 (ωt)+ r 2 ( l 2 3 r 2 ) cos 4 (ωt)( l 2 r 2 )( l 2 +3 r 2 ) cos 2 (ωt)+ r 2 ( l 2 r 2 )=0

Haciendo el cambio de variable z=cos2(ωt) se obtiene la ecuación cúbica

z 3 + ( l 2 3 r 2 ) r 2 z 2 ( l 2 r 2 )( l 2 +3 r 2 ) r 4 z+ ( l 2 r 2 ) r 2 =0

Para calcular las raíces de la ecuación cúbica

z3+az2+bz+c=0

Calculamos los valores de Q y R dados por

Q= a 2 3b 9 = 4 9 l 4 R= 2 a 3 9ab+27c 54 = 11 l 6 27 l 4 r 2 54 r 6

Si R2>Q3 tenemos una raíz real y dos complejas, en caso contrario, hay tres raíces reales. Supongamos que se cumple la primera condición. La condición R2>Q3  equivale a

27r4-33r2l2>5l4

Como r>l esta condición no se cumple

Como R2<Q3 entonces la ecuación tiene tres raíces reales

θ=arccos( R Q 3 )=arccos( 11 l 2 27 r 2 16 l 2 ) z 1 =2 Q cos( θ 3 ) a 3 = 4 l 2 3 r 2 cos( θ 3 ) l 2 3 r 2 3 r 2 z 2 =2 Q cos( θ+2π 3 ) a 3 = 4 l 2 3 r 2 cos( θ+2π 3 ) l 2 3 r 2 3 r 2 z 3 =2 Q cos( θ2π 3 ) a 3 = 4 l 2 3 r 2 cos( θ2π 3 ) l 2 3 r 2 3 r 2

Conocida las raíces de la ecuación cúbica  z se calcula el ángulo ωt.

ωt=arccos z

De las tres raíces, una es negativa (no se puede hallar la raíz cuadrada), otra es mayor que la unidad (la raíz es también mayor que la unidad y no se puede calcular en arco coseno) y la tercera, la única válida, está comprendida entre 0 y 1.

La aceleración es nula en dos instantes

ω t 1 =arccos z 3 ω t 2 =2πarccos z 3  

Ejemplo:

R=5.037         Q=7.11

Comprobamos que R2<Q3 la ecuación cúbica tiene tres raíces reales

z1=-5.172
z2
=4.028
z3
=0.1434

La velocidad es máxima (en módulo) o la aceleración es cero para

ω t 1 =arccos 0.1434 =1.81rad=67.7º ω t 2 =2πarccos 0.1434 =5.10rad=292.3º

Actividades

Se introduce

Se pulsa el botón titulado Empieza

Observamos el movimiento del pistón y la representación gráfica de su posición en función del tiempo.

Se activa la casilla titulada Gráfica

Se selecciona un botón de radio

Se pulsa el botón titulado Empieza

En estas dos últimas representaciones, se señala en el eje horizontal las posiciones angulares para las cuales la velocidad es máxima o la aceleración es cero.

stokesApplet aparecerá en un explorador compatible con JDK 1.1.

Referencias

Bacon R. H., The motion of a piston, Am. J. Phys.  (10) 1942, pp. 145-147