Ejemplo de Un Informe para Métodos Númericos

ESCUELA COLOMBIANA DE INGENIERIA JULIO GARAVITO
ANÁLISIS NUMÉRICO
GRUPO 2
PROFESOR:
ERNESTO RIVEROS OSPINA
ESTUDIANTES:
SERGIO DAVID PARDO GUERRERO
CESAR IVÁN RAMÍREZ MUNEVAR
PROYECTO FINAL
MOVIMIENTO DE UNA PELOTA DE GOLF
BOGOTÁ DC
30 DE NOVIEMBRE DE 2018
1
Tabla de contenido
1. Introducción……………………………………………………………………………3
2. Objetivos……………………………………………………………………………….3
2.1. Objetivo general…………………………………………………………………..3
2.2. Objetivos específicos……………………………………………………………..3
3. Cuerpo del trabajo…………………………………………………………………….3
3.1. Marco Teórico……………………………………………………………………..3
3.2. Procedimiento……………………………………………………………………..6
3.3. Análisis y resultados………………………………………………………………8
4. Conclusiones…………………………………………………………………………11
5. Bibliografía……………………………………………………………………………12
2
1. Introducción
En el presente trabajo se analizó el movimiento semi-parabólico de una pelota de

golf en el espacio donde se tuvo en cuenta la influencia de una corriente de viento
y las fuerzas que esta ejerce; para ello se emplearon los métodos numéricos de
Runge-Kutta y Heun los cuales nos permitieron modelar la trayectoria de la pelota
sin necesidad de usar técnicas de calculo exactas.
Con el fin de optimizar los cálculos se realizó un programa aplicado a Excel el cual
efectúa las operaciones necesarias programando una única vez las ecuaciones
planteadas para los métodos mencionados anteriormente.
Con el fin de corroborar los resultados obtenidos se hace uso de un graficador 3D

el cual nos muestra visualmente la trayectoria de la pelota.
2. Objetivos
2.1. Objetivo general

• Interpretar el movimiento de una pelota de golf aplicando métodos
numéricos para dar solución a su trayectoria.
2.2. Objetivos específicos

• Aplicar técnicas de analisis numérico para dar solución a problemas
específicos.
• Utilizar la programación como herramienta para dar solución a
problemas del ámbito ingenieril fuera del área de sistemas.
• Aprovechar los ambientes de programación como un método de
optimización del tiempo buscando la efectividad.
• Comprender la influencia del viento en el movimiento de objetos dentro
del espacio material compuesto por aire.
3. Cuerpo del trabajo
3.1. Marco Teórico
Segunda Ley de Newton
La Segunda Ley de Newton también conocida como Ley Fundamental de la

Dinámica, es la que determina una relación proporcional entre fuerza y variación
de la cantidad de movimiento o momento lineal de un cuerpo. Dicho de otra forma,
la fuerza es directamente proporcional a la masa y a la aceleración de un cuerpo.
3
La Primera Ley de Newton nos dice que para que un cuerpo altere su movimiento
es necesario que exista algo que provoque dicho cambio. Ese algo es lo que
conocemos como fuerzas. Estas son el resultado de la acción de unos cuerpos
sobre otros.
La Segunda Ley de Newton se encarga de cuantificar el concepto de fuerza. Nos

dice que la fuerza neta aplicada sobre un cuerpo es proporcional a la aceleración
que adquiere dicho cuerpo. La constante de proporcionalidad es la masa del
cuerpo, de manera que podemos expresar la relación de la siguiente manera:
ΣF = m*a
Esta ley explica qué ocurre si sobre un cuerpo en movimiento (cuya masa no tiene
por qué ser constante) actúa una fuerza neta: la fuerza modificará el estado de
movimiento, cambiando la velocidad en módulo o dirección. En concreto, los
cambios experimentados en la cantidad de movimiento de un cuerpo son
proporcionales a la fuerza motriz y se desarrollan en la dirección de esta; esto es,
las fuerzas son causas que producen aceleraciones en los cuerpos.
Fuerza de resistencia del aire
El aire es una masa que al ponerse en movimiento genera el fenómeno que

conocemos como viento; el cual genera una fuerza de fricción que se opone al
movimiento de los objetos que se desplazan dentro de este, según el principio de
Bernoulli esta fuerza esta dad por:
𝑘 𝑣 ∗ |𝑣 |
𝐹𝑟 = ( )
𝑚 2𝑔
Donde k es el factor de forma del objeto que se desplaza, m es la masa de dicho

objeto, g la gravedad y v es la velocidad de la masa.
Método de Runge-Kutta de cuarto orden
Un miembro de la familia de los métodos Runge-Kutta usado ampliamente es el de

cuarto orden. Es usado tanto que a menudo es referenciado como «RK4» o como
«el método Runge-Kutta».
Definiendo un problema de valor inicial como:
4
Entonces el método RK4 para este problema está dado por la siguiente ecuación:
Donde
Así, el siguiente valor (yn+1) es determinado por el presente valor (yn) más el
producto del tamaño del intervalo (h) por una pendiente estimada. La pendiente es
un promedio ponderado de pendientes, donde K1 es la pendiente al principio del
intervalo, K2 es la pendiente en el punto medio del intervalo, usando K1 para
determinar el valor de y en el punto Xn+h/2 usando el método de Euler. K3 es otra
vez la pendiente del punto medio, pero ahora usando K2 para determinar el valor
de y; K4 es la pendiente al final del intervalo, con el valor de y determinado por K3.
Promediando las cuatro pendientes, se les asigna mayor peso a las pendientes en
el punto medio:
Método de Heun
Un método para mejorar la estimación de la pendiente involucra la determinación y

promediado de dos derivadas para el intervalo (una en el punto inicial y otra en el
punto final).
En el método de Euler, la pendiente al inicio del intervalo se usa para extrapolar

linealmente a yi+1.
En el método de Heun la pendiente calculada en la estimación previa no es para la

respuesta final, sino para una predicción intermedia. Esta ecuación es llamada
predictor. Mejora una estimación de yi+1 que permite el cálculo de una estimación
de la pendiente al final del intervalo.
Aquí, y0i+1 es el predictor, y es la misma ecuación de Euler para encontrar

yi+1. Ésta nos sirve para calcular la pendiente y'i+1.
5
Las dos pendientes se promedian en el intervalo:
Esta pendiente promedio se utiliza para extrapolar linealmente desde yi hasta yi+1
usando el método de Euler.
Esta ecuación es conocida como ecuación correctora. El método de Heun es un

procedimiento predictor – corrector.
Se puede conseguir una mejor precisión en el resultado si hacemos varios

procesos correctores, esto lo logramos tomando yi+1 y reemplazándolo por y0i+1
en la ecuación y así encontrar un nuevo yi+1, y se repite el proceso hasta donde se
desee.
3.2. Procedimiento
Para el desarrollo del problema se parte de la sumatoria de fuerzas sobre la pelota

para hallar una funcion de aceleracion en términos de la velocidad y el tiempo en
cada una de las direcciones x, y y z. Dicha aceleracion es la derivada de la
velocidad con respecto al tiempo y a su vez esta velocidad es la derivada de la
posición respecto al tiempo, con el fin de hallar la posición de la pelota de golf para
una velocidad y tiempo determinado se usan los métodos de Runge Kutta y Heun
consecutivamente aplicando el método de Runge Kutta de 4 orden como primera
integral de la aceleracion hallando así la velocidad, y después se aplica el método
de Heun para la segunda integral y así obtener la posición.
Con el fin de aplicar los conocimientos de programación aprendidos a lo largo del

curso se propone solucionar este problema por medio de programación en Visual
Basic aplicada a Excel, en esta herramienta tecnológica se programan las
funciones que nos darán las velocidades y posiciones, donde el usuario solo
aportara datos tales como el factor de forma (k), el lapso de tiempo (h) y
condiciones iniciales como la velocidad inicial de la pelota junto con su ángulo de
elevación y polar, la velocidad del viento y su ángulo polar, haciendo así de un
problema complejo uno sumamente intuitivo para el usuario en el cual solo tendrá
que interpretar los resultados.
6
Figura 1: Movimiento descrito por la pelota en el plano xy.
Sumatorias de fuerzas en las direcciones x, y y z:
𝑑𝑣𝑧 𝑘 𝑣𝑧 ∗ |𝑣𝑧 |
= −𝑔 − ( )
𝑑𝑡 𝑚 2𝑔
𝑑𝑣𝑥 𝑘 (𝑣𝑥 − 𝑢𝑥 ) ∗ |(𝑣𝑥 − 𝑢𝑥 )|

=− ( )
𝑑𝑡 𝑚 2𝑔
𝑑𝑣𝑦 𝑘 (𝑣𝑦 − 𝑢𝑦 ) ∗ |(𝑣𝑦 − 𝑢𝑦 )|

=− ( )
𝑑𝑡 𝑚 2𝑔
Donde, u es la velocidad del viento en la dirección x o y.
Las funciones presentadas anteriormente son funciones de la forma f (t, v) las

cuales se usarán para el método de Runge Kutta
Para el caso en la dirección z, definimos
K1 = f(vzi)
K2 = f(vzi+h*K1/2)
K3 = f(vzi+h*K2/2)
7
K4 = f(vzi+h*K3)
Vsig = Vi+(h/6)*(K1+2*K2+2*K3+K4)
Dado que la velocidad es la derivada de la posición para Heun se utiliza la

velocidad anterior y la velocidad siguiente como K1’ y K2’ en la formula siguiente:
Possig = Posi+( K1’ + K2’)*h/2
Para el caso de x y y las ecuaciones son similares salvo que la velocidad es

relativa y la gravedad no se tiene en cuenta.
En el archivo de Excel se incluye una lista de puntos la cual permite al usuario

graficar la trayectoria de la pelota de golf en cualquier programa de gráficas 3D
como GeoGebra.
3.3. Analisis y Resultados
Anexo a este documento se presenta el código del programa para su posterior

revisión. Para probar el programa se usaron los siguientes datos:
Figura 2: Datos iniciales utilizados para probar el programa.
8
Se obtuvieron las siguientes graficas de posición en los distintos planos xy, xz, y
yz:
Figura 3: Movimiento reflejado sobre el plano xz.
Considerando la velocidad inicial de la pelota y la del viento se observa que,

gracias a su dirección y magnitud, en el movimiento de la pelota se ve la oposición
que presenta el viento ante la trayectoria de la pelota que, si no estuviera presente
la influencia del viento, el movimiento debería ser netamente parabólico y no semi-
parabólico como observamos en la gráfica.
Figura 4: Movimiento reflejado sobre el plano yz.
Esta grafica nos presenta una vista frontal del golfista, donde se evidencia que la
desviación de la trayectoria de la pelota es correcta según la dirección del viento,
9
ya que para un viento a 240° y un tiro a 20° respecto del eje x positivo se
esperaría que la pelota desvíe su trayectoria en dirección del eje y negativo.
Figura 5: Movimiento reflejado sobre el plano xy.
Aquí observamos una vista de techo del tiro del golfista, con esta grafica
respaldamos lo anteriormente planteado: la pelota efectivamente se desvía hacia
el eje y negativo.
A continuación, se muestra una grafica 3D generada en GeoGebra donde se

evidencia la trayectoria de la pelota de forma más clara.
Figura 6: Movimiento de la pelota en tres dimensiones.
10
4. Conclusiones
• El uso de las nuevas tecnologías permite al ingeniero optimizar y facilitar
sus tareas, además de proporcionar herramientas para quienes en futuras
tareas se enfrenten a problemas similares.
• La materia de analisis numérico amplia nuestra lista de herramientas y nos
permite ver que no estamos limitados únicamente a los cálculos exactos,
que en la vida cotidiana por lo general no se presentan.
• A diferencia del problema básico que se presenta sobre movimiento
parabólico ideal, vemos que un caso más acercado a la realidad se da
cuando consideramos fuerzas externas y no únicamente la gravedad, en
este caso un corriente de viento desvía drásticamente lo que seria la
trayectoria ideal de un movimiento parabólico.
11
5. Bibliografía y fuentes
• https://es.wikipedia.org/wiki/Método_de_Runge-Kutta
• http://aprendeenlinea.udea.edu.co/lms/moodle/mod/page/view.php?id=2453
7
• https://www.ecured.cu/Segunda_Ley_de_Newton
• HUERTA, Sarrate-Ramos, Rodríguez-Ferrán. Métodos Numéricos,
Introducción Aplicaciones y Propagación. Ediciones UPC. Primera Edición.
1998. Capítulo 9. Páginas 201-202.
• Chapra, S.C. Canale, R.P. Métodos numéricos para ingenieros. McGraw-
Hill: México D.F.
12

Ejemplo de Un Informe para Métodos Númericos

Cargado por

Copyright:

Formatos disponibles

Ejemplo de Un Informe para Métodos Númericos

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Ejemplo de Un Informe para Métodos Númericos

Cargado por

Copyright:

Formatos disponibles

ESCUELA COLOMBIANA DE INGENIERIA JULIO GARAVITO

ERNESTO RIVEROS OSPINA

SERGIO DAVID PARDO GUERRERO

CESAR IVÁN RAMÍREZ MUNEVAR

MOVIMIENTO DE UNA PELOTA DE GOLF

En el presente trabajo se analizó el movimiento semi-parabólico de una pelota de

Con el fin de corroborar los resultados obtenidos se hace uso de un graficador 3D

2.1. Objetivo general

2.2. Objetivos específicos

3. Cuerpo del trabajo

3.1. Marco Teórico

Segunda Ley de Newton

La Segunda Ley de Newton también conocida como Ley Fundamental de la

La Segunda Ley de Newton se encarga de cuantificar el concepto de fuerza. Nos

Fuerza de resistencia del aire

El aire es una masa que al ponerse en movimiento genera el fenómeno que

Donde k es el factor de forma del objeto que se desplaza, m es la masa de dicho

Método de Runge-Kutta de cuarto orden

Un miembro de la familia de los métodos Runge-Kutta usado ampliamente es el de

Definiendo un problema de valor inicial como:

Un método para mejorar la estimación de la pendiente involucra la determinación y

En el método de Euler, la pendiente al inicio del intervalo se usa para extrapolar

En el método de Heun la pendiente calculada en la estimación previa no es para la

Aquí, y0i+1 es el predictor, y es la misma ecuación de Euler para encontrar

Esta ecuación es conocida como ecuación correctora. El método de Heun es un

Se puede conseguir una mejor precisión en el resultado si hacemos varios

Para el desarrollo del problema se parte de la sumatoria de fuerzas sobre la pelota

Con el fin de aplicar los conocimientos de programación aprendidos a lo largo del

Sumatorias de fuerzas en las direcciones x, y y z:

𝑑𝑣𝑥 𝑘 (𝑣𝑥 − 𝑢𝑥 ) ∗ |(𝑣𝑥 − 𝑢𝑥 )|

𝑑𝑣𝑦 𝑘 (𝑣𝑦 − 𝑢𝑦 ) ∗ |(𝑣𝑦 − 𝑢𝑦 )|

Donde, u es la velocidad del viento en la dirección x o y.

Las funciones presentadas anteriormente son funciones de la forma f (t, v) las

Para el caso en la dirección z, definimos

Dado que la velocidad es la derivada de la posición para Heun se utiliza la

Possig = Posi+( K1’ + K2’)*h/2

Para el caso de x y y las ecuaciones son similares salvo que la velocidad es

En el archivo de Excel se incluye una lista de puntos la cual permite al usuario

3.3. Analisis y Resultados

Anexo a este documento se presenta el código del programa para su posterior

Figura 2: Datos iniciales utilizados para probar el programa.

Figura 3: Movimiento reflejado sobre el plano xz.

Considerando la velocidad inicial de la pelota y la del viento se observa que,

Figura 4: Movimiento reflejado sobre el plano yz.

Figura 5: Movimiento reflejado sobre el plano xy.

A continuación, se muestra una grafica 3D generada en GeoGebra donde se

Figura 6: Movimiento de la pelota en tres dimensiones.

También podría gustarte