Tiro parabólico con resistencia del aire proporcional a v 2

LUIS ANGEL TOVAR MARTINEZ

1. Introducción
El tiro parabólico es uno de los problemas fundamentales de la mecánica clásica. En su formulación ideal,
donde se desprecia la resistencia del aire, el movimiento del proyectil puede describirse mediante ecuaciones
analíticas simples que conducen a una trayectoria parabólica perfectamente simétrica.
Sin embargo, en situaciones reales, la interacción del proyectil con el aire introduce una fuerza de arrastre
que modifica de manera significativa el movimiento. Cuando esta fuerza es proporcional al cuadrado de la
rapidez, el sistema dinámico resultante se vuelve no lineal y acoplado, lo que impide, en general, obtener
soluciones analíticas cerradas. En este contexto, los métodos numéricos se convierten en una herramienta
indispensable.
El objetivo de este trabajo es formular el modelo físico del tiro parabólico con resistencia del aire propor-
cional a v 2 , reescribirlo como un sistema de ecuaciones diferenciales ordinarias de primer orden y resolverlo
numéricamente mediante el método de Runge–Kutta de cuarto orden.

2. Motivación: caso ideal y caso real

En ausencia de resistencia del aire (c = 0), las ecuaciones del movimiento de un proyectil lanzado con
velocidad inicial v0 y ángulo θ son

1
x(t) = v0 cos θ t, y(t) = v0 sin θ t − gt2 .
2

Estas expresiones describen una trayectoria parabólica simétrica y permiten calcular de manera exacta el
alcance, la altura máxima y el tiempo de vuelo.
Al introducir la resistencia del aire, dicha simetría se rompe y la trayectoria deja de ser parabólica.
Además, la presencia del arrastre provoca disipación de energía y reduce significativamente las magnitudes
características del movimiento, lo que justifica el uso de métodos numéricos para su análisis.

3. Modelado físico
Sea la velocidad del proyectil
q
v(t) = (vx (t), vy (t)), v = ∥v∥ = vx2 + vy2 .

Las fuerzas que actúan sobre el proyectil son:

Gravedad:
Fg = (0, −mg).

1
 Arrastre cuadrático:
Fd = −c v v,

fuerza opuesta al movimiento cuya magnitud crece como v 2 .

Aplicando la segunda ley de Newton,

dv
m = Fg + Fd ,
dt
se obtiene un sistema dinámico no lineal. Durante la subida, tanto la gravedad como el arrastre actúan hacia
abajo, mientras que durante la caída el arrastre se opone a la gravedad.

4. Ecuaciones de movimiento
Separando la ecuación vectorial en componentes cartesianas se obtiene:

dvx dvy
m = −c v vx , m = −mg − c v vy .
dt dt

Además, por definición cinemática:

dx dy
= vx , = vy .
dt dt
El término q
v= vx2 + vy2

introduce no linealidad y acopla las ecuaciones, lo que impide resolver el sistema de manera analítica en el
caso general.

5. Sistema de EDO de primer orden

Definiendo el vector de estado
u(t) = (x(t), y(t), vx (t), vy (t)),

el sistema puede escribirse como:

dx
= vx ,
dt
dy
= vy , q
dt v= vx2 + vy2 .
dvx c
=− v vx ,
dt m
dvy c
= −g − v vy ,
dt m

Este sistema tiene la forma compacta

u′ (t) = F(t, u),

adecuada para la aplicación de métodos numéricos.

2
6. Condiciones iniciales y parámetros
El lanzamiento se realiza desde el origen:

x(0) = 0, y(0) = 0.

La velocidad inicial v0 y el ángulo de lanzamiento θ determinan:

vx (0) = v0 cos θ, vy (0) = v0 sin θ.

En la simulación se utilizaron los parámetros:

g = 9,81 m/s2 , m = 1,0 kg, c = 0,04 kg/m,

con v0 = 40 m/s y θ = 45◦ . :contentReference[oaicite:3]index=3

7. Método numérico: RK4

Debido a la no linealidad del sistema, se utiliza el método de Runge–Kutta de cuarto orden (RK4), el
cual presenta un error global de orden O(h4 ) y ofrece buena precisión para sistemas acoplados. :contentRe-
ference[oaicite:4]index=4
El esquema general del método es:

h
un+1 = un + (k1 + 2k2 + 2k3 + k4 ) ,
6

donde  
h h
k1 = F(tn , un ), k2 = F tn + , un + k1 ,
2 2
 
h h
k3 = F tn + , un + k2 , k4 = F(tn + h, un + hk3 ).
2 2
:contentReference[oaicite:5]index=5

8. Código implementado en Python

A continuación se muestra el núcleo del método RK4 utilizado para actualizar el estado u = (x, y, vx , vy )
en cada paso temporal. :contentReference[oaicite:6]index=6

# Estado u = [x , y , vx , vy ]
# Definimos F (t , u ) como el lado derecho del sistema :
# dx / dt = vx
# dy / dt = vy
# dvx / dt = -( c / m ) * v * vx
# dvy / dt = -g -( c / m ) * v * vy
# donde v = sqrt ( vx ^2 + vy ^2)

k1 = F ( tk , uk ) # pendiente al inicio

3
k2 = F ( tk + 0.5* h , uk + 0.5* h * k1 ) # pendiente al medio ( con k1 )
k3 = F ( tk + 0.5* h , uk + 0.5* h * k2 ) # pendiente al medio ( con k2 )
k4 = F ( tk + h , uk + h * k3 ) # pendiente al final

uk = uk + ( h /6.0) *( k1 + 2* k2 + 2* k3 + k4 ) # c o m b i n a c i n ponderada
tk = tk + h

La interpretación del método es que aproxima la integral del campo vectorial en el intervalo [tn , tn+1 ] mediante
un promedio ponderado de pendientes, capturando mejor la curvatura de la solución que métodos más simples.
:contentReference[oaicite:7]index=7

9. Discusión cualitativa del efecto del arrastre

Comparado con el caso ideal sin resistencia del aire (c = 0), la inclusión de un arrastre proporcional al
cuadrado de la rapidez produce modificaciones importantes en el movimiento del proyectil. Entre los efectos
más relevantes se encuentran:

Disminución del alcance horizontal R y, en general, del tiempo de vuelo T , debido a la pérdida
continua de energía mecánica causada por la fricción con el aire.

La trayectoria deja de ser parabólica y se vuelve asimétrica: la fase de subida es más suave, mientras
que la caída es más pronunciada.

Existe una mayor pérdida de rapidez al inicio del movimiento cuando la velocidad inicial v0 es
grande, ya que la fuerza de arrastre escala como v 2 .

El ángulo de lanzamiento que maximiza el alcance ya no es 45◦ y suele ser menor que 45◦ en presencia
de arrastre cuadrático.

Estos efectos muestran que el arrastre del aire rompe las simetrías presentes en el caso ideal y tiene un
impacto significativo en la dinámica del proyectil.

10. Significado del coeficiente de arrastre

El parámetro c > 0 mide la intensidad de la fuerza de resistencia del aire y depende de las propiedades
físicas tanto del objeto como del medio en el que se mueve. Entre los factores más importantes que influyen
en el valor de c se encuentran:

la forma geométrica del objeto,

el área frontal expuesta al flujo de aire,

la rugosidad de la superficie,

y la densidad del aire.

Objetos más aerodinámicos presentan valores menores del coeficiente de arrastre, mientras que objetos
irregulares o con mayor área frontal experimentan fuerzas de resistencia más intensas.

4
11. Influencia de la forma en el coeficiente de arrastre
El coeficiente de arrastre c (o, de manera más habitual, el coeficiente aerodinámico Cd ) cuantifica la
interacción entre el objeto y el aire y depende fuertemente de la forma geométrica del objeto lanzado.
Objetos con formas aerodinámicas permiten que el flujo de aire se deslice con menor turbulencia alrededor
del cuerpo, reduciendo la fuerza de arrastre. Por el contrario, objetos romos o con aristas pronunciadas generan
separación del flujo, lo que provoca un aumento significativo del arrastre.
De manera orientativa, algunos valores típicos del coeficiente aerodinámico Cd son los siguientes:

Forma del objeto Cd

Esfera lisa ≈ 0,47

Cilindro ≈ 1,0
Cubo ≈ 1,05
Proyectil aerodinámico 0,04 – 0,1

Estos valores ilustran la fuerte dependencia del arrastre con la forma del objeto y explican por qué
proyectiles con geometrías distintas presentan trayectorias muy diferentes aun bajo condiciones iniciales
similares.

12. Resultados: gráficas

12.1. Gráfica 1: Trayectoria y vs. x

Figura 1: Trayectoria del proyectil con fricción cuadrática

5
 Figura 2: Velicidad vs tiempo

12.2. Gráfica 2: Rapidez v(t) vs. tiempo

13. Conclusiones e interpretación de las gráficas

El tiro parabólico con resistencia del aire proporcional a v 2 conduce a un sistema no lineal y aco-
plado de ecuaciones diferenciales, para el cual no existe una solución analítica cerrada en general.
:contentReference[oaicite:8]index=8

En la gráfica de trayectoria (Figura 1) se observa que el movimiento deja de ser parabólico: la su-
bida es más suave y la caída es más pronunciada, produciendo asimetría en la trayectoria. Además,
el alcance, la altura máxima y el tiempo de vuelo se ven reducidos por el arrastre. :contentReferen-
ce[oaicite:9]index=9

En la gráfica de rapidez (Figura ??) se aprecia que al inicio la rapidez disminuye rápidamente por la
acción combinada de la gravedad y el arrastre; luego alcanza un mínimo cerca del punto más alto
(cuando la componente vertical se anula) y durante la caída vuelve a aumentar, pero sin recuperar
el valor inicial debido a la disipación de energía. :contentReference[oaicite:10]index=10

Este comportamiento confirma que la fricción del aire extrae energía mecánica del sistema, lo cual se
manifiesta como una disminución global de la rapidez respecto al caso ideal sin fricción. :contentRefe-
rence[oaicite:11]index=11

El método de Runge–Kutta de cuarto orden (RK4) permite aproximar la solución con buena precisión
y estabilidad, siendo adecuado para este tipo de problemas no lineales.