Introducción al Análisis Numérico
Modelos Numéricos
TRABAJO PRÁCTICO Nº1
Lic. y Farm. COGGIOLA, Vivian (3D1 y 3D3)
Ing. PEREZ, Rita (3D2) 1
MODELO MATEMATICO PARA PROBLEMAS DE
INGENIERIA
Un modelo matemático es una formulación que expresa las
características de un sistema físico en término matemáticos.
Solución
EXACTA
M. Analítico
MODELO
MATEMATICO
(Modelar…) Solución
APROXIMADA
M. Numérico
Fórmula numérica –Método numérico-
Diferencia de
incrementos
DE VELOCIDAD
RESPECTO AL
TIEMPO
PASO )
𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑖𝑛𝑖𝑐𝑖𝑎𝑙
E.D.O)* PASO
𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑓𝑖𝑛𝑎𝑙
EJERCICIO Nº 1
Antes de empezar se pide analizar el caso y representar el diagrama de cuerpo libre
para modelar el problema.
a) Determine la EDO que representa el problema.
b) Determine la velocidad límite alcanzada por el cuerpo al equilibrarse las fuerzas (P)
y resistencia del aire ( R )
c) Calcular la velocidad del cuerpo a los 10s.
d) Analice resultados con diferentes pasos-
Resolución del punto a)E.D.O
“Un diagrama de cuerpo libre es un diagrama F= m*a (1) Por ley de Newton
que muestra el cuerpo elegido solo, “libre” de F=P+R (2) Por representación con el Diag de fuerzas
su entorno, con vectores que muestren las F= m*g +(-c.v) (3) Por sustitución P=m*g R= -c*v
magnitudes y direcciones de todas las fuerzas
aplicadas sobre el cuerpo por todos los F= mg- c*v(3)= m*a (1)
cuerpos que interactúan con él.” ¿ Cuánto vale a?
Fr= c*v =a
¿Positiva o negativa ? =a
Resistencia (R )
=a
a= aceleración =
Peso (P) 𝑑𝑣 𝑐
=𝑔 − ∗𝑣
𝑑𝑡 𝑚
Peso= P=m*g
¿Positiva o negativa ?
E.D.O
Resolución del punto b)velocidad
máxima
E.D.O
𝑑𝑣 𝑐
=𝑔 − ∗𝑣
𝑑𝑡 𝑚
𝑑𝑣
=0 𝑺𝒆 𝒂𝒍𝒄𝒂𝒏𝒛 ó 𝒍𝒂𝒗𝒆𝒍𝒐𝒄𝒊𝒅𝒂𝒅 𝒎 á 𝒙𝒊𝒎𝒂‼ ! m=68,1 kg
𝑑𝑡 g=9,8 m/ s²
c=12,5 m/s
V==53,3904 m/s
Resolución del punto c) velocidad del
cuerpo a los 10s
Resuelva numéricamente y calcule la velocidad del cuerpo a los 10 s utilizando un paso (h) = 2s.
E.D.O PASO
Tiempo Velocidad
inicial
0 0,000 =9,8 2-0 Primera
iteración
2 19,600 2-0
Segunda
iteración
E.D.O PASO
Tiempo Velocidad
inicial
Primera
0 0,000 =9,8 2-0 iteración
2 19,600 2-0
Segunda
4 32,005 iteración
6
¿Qué tan buena es la aproximación con
un paso = 2s?
Fórmula para hallar v(t) con la solución Fórmula para hallar v(t) con la solución
analítica numérica
¿Cómo podemos
Velocidad final a mejorar la
los 10 s=47,969 aproximación?
m/s
t V(t)
Resolución del punto d) 0
0,5
????
Calculamos con un paso= 1s, luego un paso = 0,5 s 1,0
1,5
t V(t) 2,0
0 ????? 2,5
1 3,0
2 3,5
3 4,0
4 4,5
5 5,0
6 5,5
7 6,0
8 6,5
9 7,0
¿Qué pasaría si utilizamos un paso = 0,1 s? 7,5
10 ??
8,0
8,5
9,0
9,5
10,0
Resolución del punto d)
Calculamos con un paso= 1s, luego un paso = 0,5 s y finalmente utilizando un paso = 0,1 s.
PASO 2,000 1,000 0,500 0,100
Tiempo (seg.) v analítica(m/s) v aprox. (m/s) v aprox. (m/s) v aprox. (m/s) v aprox. (m/s)
0 0,000 0,000 0,000 0,000 0,000
2 16,400 19,600 17,801 17,063 16,531
4 27,770 32,005 29,667 28,673 27,943
6 35,640 39,856 37,577 36,572 35,822
8 41,100 44,824 42,849 41,947 41,262
10 44,870 47,969 46,364 45,604 45,017
Excel
EJERCICIO Nº 2
Angulo α = 27 º
Usamos calculadora
en GRADOS
Antes de empezar se pide analizar el caso y representar el
diagrama de cuerpo libre para modelar el problema.
a) Determine la EDO que representa el problema.
b) Determine la velocidad límite alcanzada por el cuerpo al equilibrarse las fuerzas.
c) Calcular la velocidad del cuerpo a los t=10s.
d) Analice resultados con diferentes pasos-
Resolución del punto a)
𝑃𝑥
𝑃𝑦
Resolución del punto a)
= =
=
𝑃𝑥
𝑃𝑦
E.D.O
Resolución del punto b)
m=2060,35 kg
g=9,8 m/ s²
=23 º
μ=0,03
Ln(7331,8311)
Resolución del punto c)
Velocidad Velocidad
final inicial
PASO
E.D.O
Solución con un paso h= 1
t V(t)
¿CONDICION INICIAL?
0 ????
1
2
3
4
5
6
7
8
9
10
PASO 1,0
TIEMPO V APROX (m/s)
0 0,00
1 3,56
2 7,10
3 10,07
4 2,19
5 5,74
6 9,15
7 8,15
8 10,02
9 2,64
10 6,19
En Excel
PASO 1,0
TIEMPO V APROX (m/s)
0 0,00
1 3,56
2 7,10
3 10,07
4 2,19
5 5,74
6 9,15
7 8,15
8 10,02
9 2,64
10 6,19
PASO 0,5
TIEMPO V APROX
0 0,00
0,5 1,78
1 3,56
1,5 5,33
2 7,05
2,5 8,55
3 9,07
3,5 8,73
4 9,01
4,5 8,81
5 8,96
5,5 8,85
6 8,94
6,5 8,87
7 8,92
7,5 8,88
8 8,91
8,5 8,89
9 8,91
9,5 8,89
10 8,91
En Excel
Resolución del punto d)
PASO 1,0 0,5 0,1
TIEMPO V APROX (m/s) V APROX V APROX
0 0,00 0,00 0,00
1 3,56 3,56 3,55
2 7,10 7,05 6,98
3 10,07 9,07 8,79 Analice resultados con diferentes pasos-
4 2,19 9,01 8,90
5 5,74 8,96 8,90
6 9,15 8,94 8,90
7 8,15 8,92 8,90
8 10,02 8,91 8,90
9 2,64 8,91 8,90
10 6,19 8,91 8,90
EJERCICIO Nº 3
Antes de empezar se pide analizar el caso y representar el
diagrama de cuerpo libre para modelar el problema.
a) Determine la EDO que representa el problema.
b) ¿Cuál es el intervalo de tiempo que transcurre hasta que se detiene por completo?
c) Analice resultados con diferentes pasos-
Resolución del punto a)
𝑑𝑣 −𝑐𝑣 − 𝜇 𝑚𝑔 −𝑐𝑣 𝜇 𝑚𝑔 −𝑐𝑣
= = − = −𝜇 .𝑔
𝑑𝑡 𝑚 𝑚 𝑚 𝑚
Para deducir la fórmula del modelo numérico :
No olvides identificar en la fórmula 4 componentes:
𝑣 𝑡 =𝑣 𝑡 +
𝑖+ 1 𝑖 [ −𝒄𝒗 ( 𝒕 )
𝒎 ]
− μ 𝒈 .(𝒕 ¿ ¿ 𝒊+ 𝟏− 𝒕 𝒊 )¿ 1. Velocidad inicial .Es decir
2. Velocidad final .Es decir
3. E.D.O
4. Paso “h” Es decir
=0,27
¿Velocidad g=9,8
=104326
inicial v(0)?
Resolución del punto b)
¿Cuál es la velocidad inicial?
Tiempo (s) Velocidad (m/s)
0 ??? ¿Cuál es el paso en la tabla 1?
5
¿Cuántas iteraciones debemos hacer?
10
15
20 𝑣 ( 𝑡 𝑖 +1 ) =𝑣 (𝑡 ¿¿ 𝑖)+¿ ¿
25
=
30
35
40
45
Tabla 1
Resolución del punto b)
¿Cuál es la velocidad inicial?
Tiempo (s) Velocidad (m/s)
CUIDADO CON LAS
0 95,28
343 𝑘𝑚 343 𝟎𝟎𝟎 𝑚 UNIDADES!!!
5 82,05 𝑣𝑖= = =𝟗𝟓 , 𝟐𝟕𝟕𝟕𝟕=𝟗𝟓 , 𝟐𝟖 𝒎/𝒔 La velocidad se expresa en
1h 3600 𝑠
10
15
20 𝑣 ( 𝑡 𝑖 +1 ) =𝑣 (𝑡 ¿¿ 𝑖)+¿ ¿
25
30 = 82 ,049954335….
35
40
Usaremos un paso de t=5s
Resolución del punto b)
Tiempo (s) Velocidad (m/s)
¿Qué análisis podemos hacer con los valores
0 95,28 de la tabla 1?
5 82,05
10 68,82
15 55,59 𝑣 ( 𝑡 𝑖 +1 ) =𝑣 (𝑡 ¿¿ 𝑖)+¿ ¿
20 42,36
25 29,13
30 15,90
35 2,67
40 -10,56
45
Tabla 1
Tiempo (s) Velocidad (m/s)
0 95,28
5 82,05
10 68,82
15 55,59
20 42,36
25 29,13
¿Por qué usar un paso de 5s?
A fines de explorar el comportamiento del modelo.
30 15,90
¿Cuándo frena?
35 2,67 Con este paso de h=5 ,podemos afirmar que frena “entre” los 35 y 40
40 -10,56 segundos
¿Qué podemos hacer si queremos mejorar la aproximación?
Rpta:…………………………………………
Resolución del punto c)
Para investigar el comportamiento del modelo
numérico para hacer PREDICCIONES , ajustando
los valores de h , podemos recurrir a una planilla
de Excel y escribir las fórmulas TOMANDO UN
CRITERIO para la cantidad de cifras decimales
utilizadas. Por ejemplo utilizando la función de
redondeo REDONDEAR ()
PASO PASO 5 PASO 1
TIEMPO V APROX (m/s) V APROX (m/s) PASO PASO 5 PASO 1
0 95,2800 95,2800 TIEMPO V APROX (m/s) V APROX (m/s)
1 92,6340 21 39,7139
2 89,9880 22 37,0679
3 87,3420 23 34,4219 Con el paso h= 5
4 84,6960 24 31,7759 podemos decir que
5 82,0500 82,0500 25 29,1298 29,1298 frena entre los 35 y 40 s.
6 79,4039 26 26,4838
7 76,7579 27 23,8378
8 74,1119 28 21,1918
9 71,4659 29 18,5458
10 68,8199 68,8199 30 15,8998 15,8998
11 66,1739 31 13,2538
12 63,5279 32 10,6078
13 60,8819 33 7,9618
14 58,2359 34 5,3158
15 55,5899 55,5899 35 2,6698 2,6698 Con el paso h= 1
16 52,9439 36 0,0238 podemos decir que
17 50,2979 37 -2,6222 frena entre los 36 y 37s.
18 47,6519 38 -5,2682
19 45,0059 39 -7,9142
20 42,3599 42,3599 40 -10,5602 -10,5602
