Pistola de aire comprimido

En primer lugar, consideramos un modelo simplificado del movimiento de la bala en tubo del cañón propulsada por un gas a presión. Supondremos:

El gas a alta presión está inicialmente contenido en un volumen Ax0 y tiene una temperatura T0. Se libera el émbolo de masa m que representa el proyectil y el gas se expande. El gas (en color amarillo) a alta presión p y temperatura T, se comporta como un gas ideal. Su presión disminuye al aumentar el volumen V=Ax, siendo A el área de la sección transversal del tubo y x la posición del proyectil.

pV=nRTF= nRT x

Supondremos que sobre el émbolo actúa únicamente la fuerza F=pA. Aplicamos la segunda ley de Newton para describir el movimiento del proyectil

m d 2 x d t 2 = nRT x

El gas se expande adiabáticamente, de modo que su energía interna U=ncT disminuye a medida que realiza un trabajo, ncdT=-F·dx. c es el calor específico a volumen constante, 3R/2 para los gases monoatómicos y 5R/2 para los diatómicos.

nc·dT= nRT x dx ln( T T 0 )= R c ln( x x 0 )

Donde T0 es la temperatura inicial del gas y Ax0 el volumen inicial del gas

T= ( x 0 x ) R/c T 0

Multiplicando por nR y teniendo en cuenta que el volumen V=Ax, obtenemos la relación entre presiones

p= ( x 0 x ) R/c+1 p 0

Velocidad del proyectil

Escribimos la ecuación diferencial del movimiento en otra equivalente de primer orden en v y x.

v dv dx = nR T 0 x 0 R/c m x (R/c+1)

Integramos la ecuación diferencial con la condición inicial: en la posición x0 la velocidad inicial de la bala es v=0.

v 2 = 2nc T 0 m ( 1 ( x 0 x ) R/c )

Cuando x se hace grande, la velocidad tiende hacia un valor límite constante

v = 2nc T 0 m

Supongamos un gas diatómico c=5R/2, cuya presión inicial p0=10 atm, para que la presión se reduzca a p=1 atm.

1= ( 1 z ) 7/5 10z= x x 0

z=5.2. Si n=0.1 moles de gas están a una temperatura inicial T0=300 K, la masa del proyectil es m=0.1 kg, la velocidad v del proyectil es

v 2 = 2·0.1·5·8.314·300 2·0.1 ( 1 5.2 5/2 )

La velocidad del proyectil es v=110.7 m/s, cercana a la velocidad límite constante

hold on
fplot(@(z) sqrt(3*(1-z.^(-3/2))),[1,9])
fplot(@(z) sqrt(5*(1-z.^(-5/2))),[1,9])
hold off
grid on
xlabel('x/x_0')
ylabel('v/v_\infty')
legend('monoatómico','diatómico','location','south')
title('Velocidad del proyectil')

Posición del proyectil

Integramos de nuevo, para obtener la posición x de la bala en función del tiempo, t

dx dt = 2nc T 0 m ( 1 ( x 0 x ) R/c ) x 0 x dx ( 1 ( x 0 x ) R/c ) = 2nc T 0 m 0 t dt

Esta integral se resuelve haciendo el cambio de variable

x 0 x = sin 2c/R θ

Con este cambio la integral se reduce a

2 x 0 c R π/2 φ dθ sin 2c/R+1 θ = 2nc T 0 m tsinφ= ( x 0 x ) R/(2c)

Resolveremos la integral para dos casos concretos: gas monoatómico c=3R/2 y gas diatómico c=5R/2

Gas monoatómico

3 x 0 π/2 φ dθ sin 4 θ = 3nR T 0 m tsinφ= ( x 0 x ) 1/3

La integración se realiza con el cambio de variable y la identidad siguientes:

y=cotθ1+ cot 2 θ= 1 sin 2 θ dy= dθ sin 2 θ dθ sin 4 θ = dy sin 2 θ = ( 1+ y 2 )dy=( y+ y 3 3 ) =cosθ 2 sin 2 θ+1 3 sin 3 θ

Deshacemos los cambios

cosφ 2 sin 2 φ+1 sin 3 φ = 1 x 0 3nR T 0 m tsinφ= ( x 0 x ) 1/3 z 2 z 4/3 ( 2 z 1/3 +1 )= 1 x 0 3nR T 0 m tz= x 0 x z 2 +3 z 4/3 4 = 3 nR T 0 m x 0 2 t

Gas diatómico

5 x 0 π/2 φ dθ sin 6 θ = 5nR T 0 m tsinφ= ( x 0 x ) 1/5

La integración se realiza con el cambio de variable y la identidad siguientes:

y=cotθ1+ cot 2 θ= 1 sin 2 θ dy= dθ sin 2 θ dθ sin 6 θ = dy sin 4 θ = ( 1+ y 2 ) 2 dy=( y+ 2 3 y 3 + y 5 5 ) = cosθ 3+4 sin 2 θ+8 sin 4 θ 15 sin 5 θ

Desahecemos los cambios

cosφ 3+4 sin 2 φ+8 sin 4 φ 3 sin 5 φ = 1 x 0 5nR T 0 m tsinφ= ( x 0 x ) 1/5 1 3 9 z 2 +15 z 8/5 +40 z 6/5 64 = 5 nR T 0 m x 0 2 tz= x 0 x

f=@(z) sqrt((z.^2+3*z.^(4/3)-4)/3);
g=@(z) sqrt((9*z.^2+15*z.^(8/5)+40*z.^(6/5)-64)/45);
hold on
fplot(f,[1,9])
fplot(g,[1,9])
hold off
grid on
legend('monoatómico','diatómico','location','southeast')
xlabel('x/x_0')
ylabel('\tau')
view([90 -90]) %intercambian los ejes X e Y
title('posición en función del tiempo')

Velocidad final del proyectil

Consideremos un modelo de pistola de aire comprimido, consistente en un recipiente de volumen V0 que contiene gas comprimido a la presión p0. Un dispositivo libera el proyectil de masa m que se mueve a lo largo del tubo de longitud L y sección A.

Las fuerzas sobre el proyectil son:

La ecuación del movimiento del proyectil es:

m d 2 x d t 2 =ApA p atm f mv dv dx =ApA p atm f

A medida que se mueve el proyectil, el volumen del gas se incrementa V(x)=V0+Ax

Transformación isoterma

Supondremos que la expansión es isotérmica: p(V0+Ax)=p0·V0

Integramos para obtener la velocidad final v del proyectil cuando sale del tubo, x=L, sabiendo que en la posición inicial x=0, el proyectil parte del reposo, v=0

v(L)= 2 m ( p 0 V 0 ln( 1+ AL V 0 ) p atm ALfL )

La velocidad del proyectil alcanza un máximo para la longitud del tubo Lmax

L max = V 0 A ( p 0 A p a A+f 1 )

pa=1.013e5; %presión atmosférica
p0=200*pa; %presión inicial
V0=1.6e-5; %volumen inicial
A=1.6e-5; %sección del tubo
m=5.2e-4; %masa proyectil
f=0.5*p0*A; %fuerza de rozamiento

Lmax=V0*(p0*A/(pa*A+f)-1)/A
v=@(x) sqrt(2*(p0*V0*log(1+A*x/V0)-pa*A*x-f*x)/m);
fplot(v,[0,1.5*Lmax])
grid on
xlabel('L')
ylabel('v')
title('Proceso isotermo')
Lmax =    0.9802

Transformación adiabática

Si suponemos que la expansión es adiabática

p ( V 0 +Ax) γ = p 0 V 0 γ

Integramos para obtener la velocidad final del proyectil cuando sale del tubo, x=L

v(L)= 2 m ( p 0 V 0 γ1 ( 1 ( V 0 V 0 +AL ) γ1 ) p atm ALfL )

La velocidad v del proyectil alcanza un máximo para la longitud del tubo Lmax

L max = V 0 A ( ( p 0 A p a A+f ) 1/γ 1 )

pa=1.013e5; %presión atmosférica
p0=200*pa; %presión inicial
V0=1.6e-5; %volumen inicial
A=1.6e-5; %sección del tubo
m=5.2e-4; %masa proyectil
f=0.5*p0*A; %fuera de rozamiento
gamma=7/3; %gamma (gas diatómico)

Lmax=V0*((p0*A/(pa*A+f))^(1/gamma)-1)/A;
v=@(x) sqrt(2*((p0*V0/(gamma-1))*(1-(V0/(V0+A*x))^(gamma-1))-pa*A*x-f*x)/m);
fplot(v,[0,1.5*Lmax])
grid on
xlabel('L')
ylabel('v')
title('Proceso adiabático')
Lmax =    0.3402

Comparación entre ambas transformaciones

Representamos la velocidad v del proyectil en función de la presión inicial p0. Si no tenemos en cuenta la fuerza de rozamiento f. Veremos que no hay diferencia apreciable entre las velocidades finales del proyectil para un tubo de L=0.88 m de longitud

pAtm=101325; %Pa
m=19.4/1000; %kg
L=0.8825;
A=2.87233e-4; %m2
V0=4.196e-3; %volumen
gamma=1.4;
p0=110000:10000:800000;

vF=sqrt((p0*V0*(1-(V0/(A*L+V0))^(gamma-1))/(gamma-1)-A*L*pAtm)*2/m);
hold on
plot(p0,vF,'r')
vF=sqrt((p0*V0*log(1+A*L/V0)-A*L*pAtm)*2/m);
plot(p0,vF,'b')
hold off
xlabel('p(Pa)')
ylabel('v_f(m/s)')
title('Pistola de aire comprimido')

Movimiento del gas en el tubo

En la página, La fuerza que ejerce la presión atmosférica hemos tenido en cuenta el movimiento de la columna de aire a medida que la pelota se desplaza por el tubo. Se trata de un problema de masa variable

En este apartado no hemos tenido en cuenta que cuando el proyectil se ecuentra en la posición x a lo largo del tubo, una porción mg de la masa Mg de gas se desplaza detrás del proyectil

m g (x)= Ax V 0 +Ax M g M g = p 0 V 0 R T 0

La ecuación del movimiento se complicaría

d( (m+ m g )v ) dt =( p p atm )Af

Referencias

Carl E Mungan. Internal ballistics of a pneumatic potato cannon. Eur. J. Phys. 30 (2009). pp. 453-457

Disponible en https://www.usna.edu/Users/physics/mungan/Publications/Publications.php#fndtn-panel120162017

Mark Denny. The Internal Balistics of an Air GunThe Physics Teacher, Vol 49, February 2011, pp. 81-83

Z. J. Rohrbach, T. R. Buresh, M. J. Madsen. Modeling the exit velocity of a compressed air cannon Am. J. Phys. 80(1) January 2012, pp. 24-26