100

CAPITULO IV

MODELOS MATEMATICOS Y SU APLICACIÓN

PARA OBTENER DATOS BALISTICOS

4.1. MODELO ADIMENSIONAL

El presente modelo matemático pretende determinar la distribución de la

presión alrededor del proyectil en función del tiempo, así como también la
determinación de la velocidad, aceleración y distancia atravesada.

La determinación del movimiento del proyectil requiere del cálculo de la presión

en la base del proyectil, presión efectiva (PB) . En este modelo, PB se determina a
partir de la presión media.

Las ecuaciones analíticas para determinar las presiones base y “breech” desde
la presión media son llamadas ecuaciones de gradiente.

4.1.1. Gradiente de Lagrange

Este considera la ecuación de continuidad a través de la cual se puede

determinar la distribución de velocidad y la energía almacenada en un
propelente gaseoso y sólido. A través de la distribución de velocidad y ecuación
de momentum se determina la distribución de presión dada en términos de la
presión media.

La ecuación de continuidad se escribe de la siguiente manera:

∂ρ ∂ ( ρ v )
+ = 0 (4.1)
∂t ∂x
 101

donde ρ es la sensidad, v la velocidad, x la distancia y t el tiempo.

Desarrollando la segunda derivada parcial:

∂ρ ∂v ∂ρ
+ρ +v = 0 (4.2)
∂t ∂x ∂x

Si se asume δρ/δx=0, es decir, fluido incompresible:

∂v −1 ∂ ρ
= ⋅ = ft (4.3)
∂x ρ ∂t

integrando:

∫ dv = ∫ f dx
t

v = f t ⋅ x + f t1 (4.4)

A x = 0, v = 0, entonces ft1 = 0

En un tiempo t, x=xp y v=vp, donde xp es la posición de la base del

proyectil y vp es la velocidad del proyectil.

entonces:
vp
vp = f p ⋅ xp o ft =
xp

por lo tanto,
 102

vp
v = ⋅x (4.5)
xp

Es necesario conocer la energía almacenada en el movimiento de los

gases del propelente y del propente sólido.
La energía cinética se define como:

v 2 ⋅ dm
d Ep = y dm = ρ ⋅ A ⋅ dx
2

donde “A” es el área del orificio y m es la masa (se aume que el diámetro
de la cámara es igual al área de la perforación)

entonces:
x
ρ ⋅ A ⋅ v2
Ep = ∫
0
p
2
dx (4.6)

sustituyendo la ecuación 4.5 en la 4.6

xp
ρ ⋅ A ⋅ v 2p ⋅ x 2
Ep = ∫ 0 2 ⋅ x 2p
dx (4.7)

ρ ⋅ A ⋅ v 2p ⋅ x p
Ep = (4.8)
6

Si, Ct = ρ ⋅ A ⋅ xp (4.9)

donde Ct es la masa total de propelente

 103

entonces:
C t ⋅ v 2p
Ep = (4.10)
6

Consideremos ahora la ecuación de momentum:

∂v ∂v − 1 ∂P
+ v⋅ = ⋅ (4.11)
∂t ∂x ρ ∂x

x
usando la ecuación (4.5), v = vp ⋅
xp

∂v x ⋅ ap x ⋅ v 2p
= − (4.12)
∂t xp x 2p

∂v vp
= (4.13)
∂x xp

donde ap es la aceleración del proyectil.

sustituyendo 4.12 y 4.13 en 4.11 se obtiene:

x −1 ∂ P
⋅ ap = ⋅ (4.14)
xp ρ ∂x

integrando desde x = 0 a x

− [ρ ⋅ a p ⋅ x 2 ]
P( x) − P( 0 ) = (4.15)
2 ⋅ xp

Si, C t = ρ ⋅ A ⋅ xp y P ( 0 ) = PBR
 104

donde PBR es la presión “breech” (en el cierre)

Ct ⋅ a p ⋅ x 2
P( x) = PBR − (4.16)
2 ⋅ A ⋅ x 2p

entonces:

ap =
[P B
− br ] ⋅ A
(4.17)
mp
donde PB es la presión en la base del proyectil, br es la presión producto
del engarzamiento del proyectil, y mp es la masa del proyectil.

entonces:
C t ⋅ [PB − br ] ⋅ x 2
P( x) = PBR − (4.18)
2 ⋅ m p ⋅ x 2p

a x = xp

C t ⋅ [PB − br ]
PB = PBR − (4.19)
2 ⋅ mp

La presión media se define como:

xp

∫ P dx
PM = 0
(4.20)
xp

sustituyendo la ecuación 4.18 en la 4.20

 105

xp
⎡ C t ⋅ [PB − br ] ⋅ x 2 ⎤
∫0 ⎢ PBR − 2 ⋅ m ⋅ x 2 ⎥
⎣⎢ ⎥⎦
=
p p
PM (4.21)
xp

e integrando:

C t ⋅ [PB − br ]
PM = PBR − (4.22)
6 ⋅ mp

resolviendo simultáneamente las ecuaciones 4.22 y 4.19

C t ⋅ [PB − br ] C ⋅ [PB − br ] C t ⋅ [PB − br ]

PBR = PM + ⇒ PB = PM + t −
6 ⋅ mp 6 ⋅ mp 2 ⋅ mp
C t ⋅ [PB − br ] ⎡ Ct ⎤ C t ⋅ br
⇒ PB = PM − (*) ⇒ PB ⋅ ⎢1 + ⎥ = PM +
3 ⋅ mp ⎢⎣ 3 ⋅ m p ⎥⎦ 3 ⋅ mp

se obtiene:
C t ⋅ br
PM +
3 ⋅ mp
PB = (4.23)
Ct
1+
3 ⋅ mp

y de la ecuación 4.19:
C t ⋅ [PB − br ]
PBR = PB + (4.24)
2 ⋅ mp
 106

Ahora, buscando PM y PB como función de PBR, sin barrera de resistencia

De la ecuación 4.19, considerando br = 0:

PB ⋅ 2 ⋅ m p = 2 ⋅ m p ⋅ PBR − C t ⋅ PB ⇒ PB ⋅ [2 ⋅ m p + C t ] = 2 ⋅ m p ⋅ PBR

se obtiene:
PBR
PB = (4.25)
Ct
1+
2 ⋅ mp

se tiene de (*):

C t ⋅ PB ⎡ Ct ⎤
PB = PM − ⇒ PB ⋅ ⎢1 + ⎥ = PM
3 ⋅ mp ⎢⎣ 3 ⋅ m p ⎥⎦

entonces:
Ct
1+
3 ⋅ mp
PM = ⋅ PBR (4.26)
Ct
1+
2 ⋅ mp
 107

4.2. CÁLCULOS TERMODINÁMICOS

Determinación de la temperatura del gas a través de la ecuación de energía y

pérdida de energía y presión media a partir de la temperatura del gas y ecuación de
estado con propiedades termodinámicas, geometría y velocidad de quemado del
propelente.

4.2.1. Ecuación de energía

Q = U + E pt + pérdidas (4.27)

donde, Q es la energía entregada por el quemado de propelente, U es la

energía interna del gas propelente y Ept es el trabajo interno entregado por el
proyectil.

Se consideran las siguientes pérdidas:

Ep = Energía usada para mover el propelente sólido y el gas propelente.

Epr = Energía usada para rotar el proyectil.
Ebr = Energía usada para desengarzar el proyectil y superar la fricción.
Er = Energía usada a causa del retroceso.
Ed = Energía usada para mover el aire frente al proyectil.
Eh = Energía usada para calentar las paredes de la cámara.
entonces:

Q = U + E pt + E p + E pr + Ebr + Er + Ed + Eh (4.28)

Si el propelente se quema en un recipiente cerrado, entonces los gases

producidos tendrán una única temperatura llamada temperatura de llama
adiabática, TF.
 108

La energía entregada por un propelente es:

TF

Q = m g ⋅ ∫ C v dT (4.29)
0

donde, mg es la masa de gas quemado y Cv es el calor específico a volumen

constante. Para más de un propelente se asume que la energía total es igual
a la suma de las energías proporcionadas por los propelentes individuales, es
decir:
T Fi

Q = ∑ m gi ⋅ ∫ C vi dT
i
(4.30)
0

La energía interna de los gases del propelente será en general menor

que la energía total causada por el trabajo entregado por el proyectil y sus
pérdidas.
T

U = ∑ m gi ⋅ ∫ C vi dT
i
(4.31)
0

por lo tanto:

T Fi T

∑ m gi ⋅ ∫ C vi dT
i
= ∑ m gi ⋅ ∫ C vi dT + E pt + pérdidas
i
(4.32)
0 0

ó
T Fi

∑ m gi ⋅ ∫ C vi dT
i
= E pt + pérdidas (4.33)
0

Si se asume que Cvi es constante en un relativamente pequeño rango de

temperatura de T a TFi:
 109

TFi

∑ m gi ⋅ C vi ⋅ ∫ dT
i
= E pt + pérdidas (4.34)
T

integrando:

∑m gi
⋅ C vi ⋅ [TFi − T ] = E pt + pérdidas (4.35)
i

despejando la temperatura T:

∑m gi
⋅ C vi ⋅ TFi − E pt − pérdidas
T = i
(4.36)
∑ m gi ⋅ C vi
i

Para poder encontrar el valor de T a partir de la ecuación 4.36 es

necesario conocer Cvi, TFi, mgi, Ept y las pérdidas.

4.2.2. Evaluación de pérdidas

Se puede utilizar algunas correlaciones de mecánica racional para

evaluar las pérdidas.

a) La energía perdida a causa de la traslación del proyectil:

m p ⋅ v p2
E pt =
2

donde, mp y vp son la masa y la velocidad del proyectil respectivamente.

 110

b) La energía perdida debido a la rotación del proyectil:

I p ⋅ ω p2
E pr =
2

m p ⋅ d b2
donde el momento de inercia es, I p = y la velocidad angular es:
8
2 ⋅ π ⋅ Tw ⋅ v p
ωp =
Db

entonces:
π 2 ⋅ m p ⋅ v p2 ⋅ Tw2
E pr =
4

c) La energía perdida debido al gas propelente y el propelente no quemado

en movimiento es:

c t ⋅ v p2
Ep = donde ct es la masa total de propelente.
6

d) La energía perdida por el trabajo para vencer la resistencia debido al

engarzamiento y fricción es:

t
E br = A ⋅ ∫ br dx = A ⋅ ∫ br ⋅ v p dt
0

e) La energía debida al retroceso es:

m rp ⋅ v rp2
Er = donde, mrp y vrp son la masa y la velocidad de
2
las partes retrocedidas respectivamente.
 111

f) La energía perdida debido a la resistencia del aire:

t
E d = A ⋅ ∫ Pg dx = A ⋅ ∫ Pg v p dt donde Pg es la presión promedio
0

del gas en el frente del proyectil.

g) La energía perdida debido a la transferencia de calor en las paredes de la

cámara está dada por:

t
E h = ∫ Q dt
0

donde Q es el flujo de calor en las paredes de la cámara, con

Q = Aw ⋅ h ⋅ [T − Tc ] , siendo Aw el área de las paredes de la cámara más el

área del tubo del arma explotada por el gas propelente, h es el coeficiente
de transferencia de calor y Tc la temperatura de la pared de la cámara.

v 0 ⋅ π ⋅ Db
Con Aw = + 2 ⋅ A + π ⋅ Db ⋅ s , donde Db es el diámetro interior
A
del tubo.

De Nordheim, Soodak&Nordheim, h = λ ⋅ C p ⋅ ρ ⋅ v + h0 , donde

λ = [13.2 + 4 ⋅ log10 ⋅ [100 − Db ]] , asumiendo Db en metros, Cp es la

−2

capacidad calórica promedio a presión constante de los gases del

vp
propelente, ρ es la densidad media del gas, v = , h0 es el coeficiente de
2
transferencia de calor convectivo libre.
 112

La temperatura en la pared de la cámara es:

E h + f ⋅ E br
Tc = + T0
C pw ⋅ ρ w ⋅ Aw ⋅ Dw

Donde, Cpw es la capacidad calórica en la pared del tubo, ρw es la

densidad de la pared del tubo, Dw es el espesor de la pared del tubo a
través del cual fluye el calor, T0 es la temperatura inicial de la pared del
tubo y f es la fracción de energía por la barrera de fricción que va en las
paredes del tubo.

Cada propelente tiene una única característica llamada energía

específica o ímpetu definido como:
TFi
Fi = R⋅ (4.37)
Mwi

Donde, Fi es el ímpetu, TFi la temperatura de llama adiabática, Mwi

es el peso molecular de los gases del propelente y R es la constante
universal de los gases.

Por termodinámica:

R ⎡ C pi ⎤
= C pi + C vi = C vi ⋅ ⎢ − 1⎥ (4.38)
Mwi ⎣ C vi ⎦

y definiendo Gama como:

C pi
Γ = (4.39)
C vi
 113

entonces:

= C vi ⋅ [Γ − 1]
R
(4.40)
Mwi

sustituyendo 4.40 en 4.37:

Fi
C vi = (4.41)
TFi ⋅ [Γi − 1]

y finalmente sustituyendo 4.41 en 4.36, tenemos:

m gi ⋅ Fi
∑i Γi − 1
− E pt − pérdidas
T = (4.42)
m gi ⋅ Fi
∑ T ⋅ [Γ
i − 1]
Fi i

4.2.3. Energía del propelente

La energía del propelente por unidad de masa es:

TF
Q
Ct
= ∫C
0
v
dT (4.43)

Asumiendo que la capacidad calórica a volumen constante es

aproximadamente la misma que a baja temperatura, se tiene:
 114

Q R ⋅ TF F
= Cv ⋅ TF = = (4.44)
Ct Mw⋅ [Γ − 1] Γ −1

4.2.4. Ecuación de estado de covolumen

La ecuación de estado del covolumen es:

mg ⋅ R ⋅ T
P ⋅ [V F − m g ⋅ b] = (4.45)
Mw

Donde, P es la presión del gas del propelente, VF es el volumen

disponible para el gas propelente, b es el covolumen (una medida de volumen
ocupada por moléculas de gas propelente) y VF el volumen de la cámara menos
el propelente sólido remanente más el volumen abierto por el movimiento del
proyectil, es decir:

⎡C mg ⎤
= V0 − ⎢ − + A⋅s
ρ ⎥⎦
VF (4.46)
⎣ρ

Donde, C es la masa inicial de propelente, ρ la densidad del propelente

sólido, A es el área de la perforación y s la distancia que el proyectil ha
atravesado.

En términos del Ímpetu (o energía específica), la ecuación de estado es:

mg ⋅ F ⋅ T
P =
[V − m g ⋅ b] ⋅ TF
(4.47)
F
 115

Si tenemos más de un propelente se asume la ley de Dalton de

presiones parciales:
Fi ⋅ m gi ⋅ T
PM = ∑P
i
i = ∑i V c ⋅ T Fi
(4.48)

con:

⎡C mgi ⎤
Vc = V0 − ∑ ⎢ i − ⎥ − ∑ mgi ⋅ bi + A ⋅ s (4.49)
i
⎣ ρi ρi ⎦ i

4.2.5. Obtención de la presión generada en función del tiempo

Es necesario obtener Fi, bi, TFi, Γi y mgi para realizar los cálculos que
permitan obtener las variables que persigue obtener la balística interior. Existen
códigos computacionales termodinámicos que permiten calcular las
propiedades del propelente. Hoy en día, la mayor parte de los datos de balística
interior son adquiridos por estos métodos. Dichos códigos requieren de datos
tales como la composición química y los calores de formación de cada
constituyente.

Históricamente las propiedades

termodinámicas se han obtenido
de medidas experimentales con
la bomba manométrica. Una
bomba manométrica es un
recipiente cerrado muy resisten-
te a las altas presiones genera-
das cuando un propelente es
quemado en su interior.
Figura N°4.1 “Esquema bomba manométrica”

Este equipo permite medir la presión generada en función del tiempo.

 116

mg ⋅ R ⋅ T
Si, P ⋅ [VF − mg ⋅ b ] = , entonces:
Mw

VF ⎡ 1 ⎤ ⎡ T ⎤
= ⋅ R +b (4.50)
m g ⎢⎣ P ⎥⎦ ⎢⎣ Mw ⎥⎦

Si se realiza medidas con diferentes cargas de propelente se puede

graficar VF/mgi en función de 1/Pmax. En esta gráfica el intercepto en la abscisa
es el covolumen y la pendiente es el ímpetu.
 117

4.3. VOLUMEN DE QUEMADO

Para determinar mgi, la masa del propelente quemada en función del tiempo, se
asumirán diferentes geometrías de propelente a través de las cuales se evalúa la
distancia de propelente que se ha quemado. Asumiendo que el propelente se quema
en forma normal a todas las superficies (capas paralelas) se puede obtener la masa
de propelente quemada.

Ejemplo 4.1 Desarrollar una expresión que permita obtener la fracción de la masa
quemada y distancia de quemado de un propelente cuya geometría de grano
corresponde a la de tipo “cordita”

2
⎡D ⎤
Un grano antes de ser quemado, tiene un volumen inicial, V0 = π ⋅ ⎢ ⎥ ⋅ L
⎣2⎦

Figura N°4.2 “Dimensiones características del grano de propelente cordita”

después de quemada una distancia x, su volumen es:

⎡D − 2 ⋅ x⎤
2

Vx = π ⋅ ⎢
⎣ 2 ⎥⎦ ⋅ (L − 2 ⋅ x ) (4.51)
 118

L – 2x x

D – 2x

Figura N°4.3 “Dimensiones características después de quemado una distancia x”

entonces, la fracción de masa quemada zi es:

Vx V0 − Vx
zi = 1 − =
V0 V0

Si la masa original del propelente es Ci, entonces la cantidad de propelente

quemado en el tiempo que el propelente tiene una distancia de quemado x es:

mgi = C i ⋅ z i

Ejemplo 4.2 Desarrollar una expresión que permita obtener la fracción de la masa
quemada y distancia de quemado de un propelente cuya geometría de grano
corresponde a la de tipo “uniperforada” , explicando en concepto de desmenbramiento
para una geometría de grano heptatubular.

Desarrollo:
⎡ ⎡D ⎤
2
⎡d ⎤ ⎤
2

a) Su volumen inicial es: V0 = ⎢π ⋅ ⎢ ⎥ −π ⋅ ⎢ ⎥ ⎥ ⋅L

⎣⎢ ⎣2⎦ ⎣ 2 ⎦ ⎦⎥
 119

L – 2x

D D – 2x

Figura N°4.4 “Dimensiones características del grano de propelente uniperforada”

su volumen después del quemado de una distancia x es:

⎡ ⎡ D − 2 ⋅ x⎤ 2 ⎡d + 2 ⋅ x ⎤ ⎤
2

Vx = ⎢π ⋅ ⎢
⎣ 2 ⎥⎦ − π ⋅ ⎢⎣ 2 ⎥⎦ ⎥ ⋅ (L − 2 ⋅ x) (4.52)
⎣ ⎦

b) El desmembrado ocurre cuando las perforaciones de un grano multiperforado

se queman hasta el punto en que ellos se tocan. La trama se define de
manera que las distancias quemadas sean iguales en el desmembramiento.

Figura N°4.5 “Vista frontal del grano de propelente heptatubular”

 120

Más allá del punto de desmembramiento, el área superficial y por tanto el

volumen, para el mismo paso de incremento el quemado continua creciendo.
Después del desmembramiento el área superficial y el volumen decrecerán
rápidamente. Esto ocurre después que el 80% del propelente ha sido
quemado. Una granulación de siete perforaciones tendrá una velocidad más
alta que una granulación de una perforación.

4.3.1. Ecuaciones que rigen la velocidad de quemado

Por definición, el volumen del propelente está dado por:

V = ∫ S dx
0
(4.53)

donde, V es el volumen del propelente, s el área superficial del propelente

y x la distancia dentro del grano.

S
dx

Figura N°4.6 Avance quemado grano propelente después de un tiempo t

m = ρ p ⋅V (4.54)

donde, m es la masa de propelente en el volumen V y ρp la densidad del

propelente sólido.
 121

entonces:
x

m = ρ p ⋅ ∫ S dx (4.55)
0

diferenciando:

dm dx
= ρp ⋅S ⋅ (4.56)
dt dt

dx
definiendo la velocidad de quemado como r =
dt
dm
r = dt (4.57)
ρp ⋅S

de la ecuación de energía, asumiendo mezcla completa:

i
[ i
] p
[ pi
]
mi Cv Tf − T − mp Cv Tf − T − ma Cv [Ta − T ] = L
a
( 4.58)

donde:
mi = Masa del gas iniciador.
mp = Masa del gas propelente.
ma = Masa del aire.
Tfi = Temperatura de llama del gas iniciador.
Tfp = Temperatura de llama del gas propelente.
Ta = Temperatura inicial del aire.
Cvi = Calor específico del gas iniciador.
Cvp = Calor específico del gas propelente.
Cva = Calor específico del aire.
T = Temperatura de la mezcla de gases con la energía perdida.
 122

L = Energía perdida.

despejando T y sustituyendo en la ecuación de estado de covolumen:

⎡ mp mi ⎡ 1 ⎤ ⎡1 ⎤ ⎤
P ⎢V0 − 0
− + mp ⎢
o
− b p ⎥ + mi ⎢ − bi ⎥ − m a ba ⎥ =
⎢⎣ ρp ρi ⎣⎢ ρ p ⎦⎥ ⎣ ρi ⎦ ⎥⎦
(4.59)
m p ⋅ C v ⋅ T f + mi ⋅ C v ⋅ T f + m p ⋅ C v ⋅ T f
m ⋅ Rh ⋅
p p i i p p

m p ⋅ C v + mi ⋅ C v + m a ⋅ C v
p i a

donde:
P = Presión en la bomba
V0 = Volumen inicial de la bomba vacía.
m = mi+mp+ma
mpo = masa inicial de propelente.
mio = masa inicial del iniciador.
bp = covolumen del propelente.
bi = covolumen del iniciador.
ba = covolumen del aire.

m p ⋅ R hp + m i ⋅ R hi + ma ⋅ R ha
Rh = (4.60)
m

R R R
donde: R hp = ; R hi = y R ha =
Mw p Mw i Mw a

Mwp = Peso molecular del gas propelente.

Mwi = Peso molecular del gas iniciador.
Mwa = Peso molecular del aire.
 123

resolviendo 4.60 para mp:

− B − B2 − 4 ⋅ A ⋅ C
mp = (4.61)
2⋅ A

donde:
⎡ 1 ⎤
A = P⋅⎢ − b p ⎥ ⋅ Cv p − R hp ⋅ Cv p ⋅ Tfp
⎢⎣ ρ p ⎥⎦

⎡ 1 ⎤ ⎡ ⎡1 ⎤ ⎤
[
B = mi C v + m a C v P ⎢ ]
− b p ⎥ + C v P ⎢V0 −
m p 0 mi 0
− + m i ⎢ − bi ⎥ − m a ba ⎥
⎢⎣ ρ p ⎥⎦ ⎢⎣ ρp ρi ⎣ ρi ⎦ ⎥⎦
i a p

p p
[
− C v T f [m i ⋅ R hi + m a R ha ] − R hp m i C v T f + m a C v Ta − L
i a
]

⎡ ⎡ 1 ⎤ ⎤
C = ⎢m i ⋅ ⎢
mp mi
[
− bi ⎥ + V0 − 0 − o − ma ⋅ ba ⎥ ⋅ P ⋅ m i ⋅ Cv i + ma ⋅ Cv a ]
⎢⎣ ⎢⎣ ρ p ⎥⎦ ρp ρi ⎥⎦
[
− [mi ⋅ R hi + ma ⋅ R ha ] ⋅ mi ⋅ Cv i ⋅ Tff + ma ⋅ Cv < ⋅ Ta − L ]

si mp= f (P)

una alternativa de función de energía perdida, L, es:

L = a + k ⋅ [P − Po ]

siendo:
Po = Presión a la cual comienza el cálculo, usualmente la
presión calculada para el quemado de todo el iniciador.
a = Energía perdida sobre Po.
P = Presión medida.
 124

∂H = 1 + k ⋅ [Pmax − Po ]

donde:
δH = Energía total perdida a Pmax.
Pmax = Presión máxima observada.

La energía perdida a Pmax es la energía que se pierde para ir de la

presión teórica (Pteo) a la presión máxima.

La energía sin pérdida de calor está dada por:

E Pteo = m p ⋅ Cv p ⋅ Tfp + m i ⋅ Cv i ⋅ Tfii + ma ⋅ Cv a ⋅ Ta

de la ecuación 4.58 con L = 0

[
E pteo = T0 ⋅ m p ⋅ Cv p + mi ⋅ Cv i + ma ⋅ Cv a ]

donde T0 es la temperatura de la mezcla gaseosa sin pérdidas de

calor.

[
E pmax = Energía con pérdida de calor = T p ⋅ m p ⋅ Cv p + m i ⋅ Cv i + ma ⋅ Cv a ]
siendo, Tp es la temperatura del gas a Pmax.

entonces:
[ ]
E pteo − E pmax = [T0 − T p ]⋅ m p ⋅ Cv p + m i ⋅ Cv i + ma ⋅ Cv a = ∂H
 125

Usando la ecuación de estado

[T − Tp ] =
[Pteo − Pmax ] ⋅V
m ⋅ Rh
0

entonces:

∂H =
[Pteo − Pmax ] ⋅V ⋅ [m ⋅ C v p + m i ⋅ Cv i + m a ⋅ C v a ]
m ⋅ Rh
p

donde, V = V0 − mi ⋅ bi − m p ⋅ bp − ma ⋅ ba

implica que:
∂H − a
k=
Pmax − P0

La ecuación 4.61 puede ser diferenciada para obtener dm/dt, donde

dm dP
= f (P , ) y con la forma de la función resuelta para la superficie S, se
dt t
obtiene r.

Es importante hacer notar que la curva presión-tiempo puede

diferenciarse numéricamente para obtener dP/dt. De manera equivalente, se
puede obtener mt a Pt y si en la función se incluye la distancia quemada, xt, y
equivalentemente se puede obtener mt+δt a P t+δt para la distancia quemada r t+δt.
entonces:
x t + δt − x t
r =
δt
 126

4.3.2. Efectos de la temperatura

Si las velocidades de quemado de un propelente son medidas a altas

temperaturas, las velocidades de quemado serán altas. La diferencia entre un
quemado a temperatura ambiente y 60ºC, genera diferencias del orden del 10-
20% en la velocidad de quemado. Esta sería la causa de la presión máxima
trasera para un arma

4.3.2.1. Retroceso

En general, el retroceso es ignorado debido a que es un efecto muy

pequeño en el ciclo de balística interior. El efecto es pequeño debido a
que la masa que retrocede es muchos órdenes de magnitud más
pesada que el proyectil, por tanto, la masa que retrocede obtiene una
muy pequeña velocidad, y consecuentemente una muy pequeña
energía cinética.

Para las situaciones en que esto no es correcto, un modelo simple

propuesto para el retroceso es:

⎡ A ⎤ ⎡ RP ⎤
a rp = ⎢ ⎥ ⋅ ⎢ PBR − − br ⎥
⎢⎣ m rp ⎥⎦ ⎣ A ⎦

donde, arp y mrp son la aceleración y la masa de las partes

retrocedidas y RP es la fuerza resistiva al movimiento de retroceso.

RP0
Se debe tener en cuenta que arp=0 mientras PBR 〉 + br , donde
A
RP0 es la fuerza requerida para resistir la barra precargada.
 127

RP

RP0

tR0 Tiempo

Figura N°4.7 “Gráfico energía retroceso vs tiempo”

RP está determinada por la linealidad interpolada entre RP0 y la

línea de equilibrio.