Guía01 2 2022 1

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 6

Universidad de Sucre & Departamento de Matemáticas

Curso: Métodos Numéricos – Semestre 02 - 2023

Guía de clase 2. Ingeniería, modelación matemática y métodos numéricos

Modelación matemática: es un proceso dinámico en el que se identifican los factores más relevantes de un fenómeno natural o artificial, los cuales
se expresan en una relación funcional que debe ser contrastada con la realidad, para realizar ajustes que permitan una mejor predicción del
fenómeno.

El proceso de modelación matemática no termina con la formulación funcional de la relación entre factores, 𝑦 = 𝑓(𝑥1 , 𝑥2 , … , 𝑥𝑛 ), es cíclico por
cuanto siempre surgen nuevos datos y/o contextos que permiten generalizar el modelo matemático.
El modelo matemático de un fenómeno es esquema simplificado e idealizado de aquel, constituido por símbolos y operaciones matemáticas. Es
una traducción de la realidad física en términos matemáticos, es decir, una forma de representar cada uno de los factores que intervienen en un
cierto proceso físico mediante objetos matemáticos y sus interrelaciones.

Ejemplo 1 Modelo matemático para el periodo de oscilación de un péndulo.

𝐿
𝑇 = 2𝜋√
𝑔

𝐿 𝜃0 𝑑𝜃 𝐿 1 2 𝜃0 1∙3 2 𝜃0 1∙3∙5 2 𝜃0
𝑇 = 2𝜋√ ∫ = 2𝜋√ (1 + ( ) 𝑠𝑖𝑛2 + ( ) 𝑠𝑖𝑛4 + ( ) 𝑠𝑖𝑛6 +∙ ∙ ∙)
𝑔 0 𝑔 2 2 2∙4 2 2∙4∙6 2
√𝑠𝑖𝑛2 (𝜃0 ) − 𝑠𝑖𝑛2 (𝜃)
2 2

𝐼𝑜 1 2 𝜃
2 0
1∙3 2 𝜃
4 0
1∙3∙5 2 𝜃0
𝑇 = 2𝜋√ (1 + ( ) 𝑠𝑖𝑛 +( ) 𝑠𝑖𝑛 +( ) 𝑠𝑖𝑛6 +∙ ∙ ∙)
𝑀𝑔𝐿 2 2 2∙4 2 2∙4∙6 2

El modelo matemático aún es incompleto, falta incorporar en la expresión funcional, 𝑻 = 𝒇(𝑰𝟎 , 𝑴, 𝑳, 𝜽𝟎 ), el efecto de la fricción y la temperatura.
Ejemplo 2 Modelo matemático para la distribución de la temperatura sobre una placa conductora de calor, en régimen estacionario.
Para modelar el proceso de distribución de temperatura se tiene que resolver una ecuación diferencial parcial; esto puede hacerse de forma
analítica ó numérica, mediante un sistema de ecuaciones, 𝐴𝑈 = 𝑏, que puede tener muchas incógnitas. En la secundaria a veces resolvíamos
sistemas de ecuaciones de tres incógnitas y hacer las operaciones para resolverlas era tedioso, como aquí son muchísimas incógnitas, desarrollamos
un método numérico, que implementado en una computadora, permite obtener una buena aproximación a la matriz 𝑈.

Para que la computadora haga este trabajo, se desarrollan algoritmos y software. Por eso se habla de modelación matemática en ambiente
computacional ó computación científica.

La computación científica: la respuesta al anhelo ancestral de predecir ó controlar a la naturaleza.


Solución del problema 1.3: en vez de la relación lineal de la ecuación (1.7), elija modelar la fuerza hacia arriba sobre el paracaidista
como una relación de segundo orden, 𝐹𝑈 = −𝑐 ′ 𝑣 2 , donde 𝑐 ′ =coeficiente de arrastre de segundo orden (kg/m), a. Con el empleo del
cálculo, obtenga la solución de forma cerrada para el caso en que al inicio el saltador se encuentra en reposo, (𝑣 = 0 en 𝑡 = 0), b.
Repita el cálculo numérico del ejemplo 1.2 con los mismos valores de condición inicial y de parámetros. Utilice un valor de 0.225 kg/m
para 𝑐 ′ .
a. Solución analítica:
𝑑𝑣 𝑑𝑣 𝑐′ 𝑑𝑣
𝑚 ∙ 𝑎(𝑡) = 𝑚 ∙ 𝑑𝑡 = 𝑚𝑔 − 𝑐 ′ 𝑣 2 , 𝑑𝑡 = 𝑔 − 𝑚 𝑣 2 ↔ 𝑐′
= 𝑑𝑡
𝑔− 𝑣 2
𝑚

𝑚
𝑐 ′ 𝑑𝑣 = 𝑚 ∫
𝑣 𝑡 𝑣 𝑣
𝑑𝑣 𝑑𝑣 𝑚 𝑣 𝑑𝑣 𝑚 𝑣
∫ ′
= ∫ 𝑑𝑡 → 𝑡 = ∫ 𝑚𝑔 ′ 𝑚𝑔 = ′
∫ = 𝑡𝑎𝑛ℎ −1 ( )
0 𝑐 2 0 0 − 𝑣2 𝑐 0 ′ − 𝑣2 𝑐 0 𝑚𝑔 2 𝑚𝑔 𝑚𝑔
𝑔− 𝑣 𝑐 ′ 𝑐 (√ ′ ) − 𝑣 2 ′
𝑐√ ′ √ ′
𝑚 𝑐 𝑐 𝑐

𝑚𝑔 𝑚𝑔 𝑚𝑔
𝑐′√ 𝑐′ 𝑣 𝑐′√ 𝑐′ 𝑣 𝑚𝑔 𝑐′√ 𝑐′
↔ 𝑡 = tanh−1 ↔ tanh 𝑡 = → 𝑣(𝑡) = √ ′ tanh 𝑡
𝑚 𝑚𝑔 𝑚 𝑚𝑔 𝑐 𝑚
√ √ ′
( 𝑐′ ) ( ) 𝑐 ( )

Para los valores, 𝑚 = 68.1, 𝑐 ′ = 0.225, 𝑔 = 9.81, 𝒗(𝒕) = 𝟓𝟒. 𝟒𝟗𝟎𝐭𝐚𝐧𝐡(𝟎. 𝟏𝟖𝟎 ∗ 𝒕), modelo matemático para el paracaidista.
b. Solución numérica:

𝑑𝑣 ∆𝑣 𝑣(𝑡 ) − 𝑣(𝑡𝑖 )
𝑚𝑎(𝑡) = 𝑚 ≈𝑚 = 𝑚 𝑖+1 = 𝑚𝑔 − 𝑐′ 𝑣2(𝑡𝑖)
𝑑𝑡 ∆𝑡 𝑡𝑖+1 − 𝑡𝑖
𝑣(𝑡𝑖+1) − 𝑣 (𝑡𝑖 ) 𝑐′
= 𝑔 − (𝑣(𝑡𝑖 ))2
𝑡𝑖+1 − 𝑡𝑖 𝑚
De la igualdad anterior se obtiene una expresión que permite estimar la velocidad:

𝑐′ 2
𝑣 (𝑡𝑖+1 ) = 𝑣 (𝑡𝑖 ) + [𝑔 − (𝑣(𝑡𝑖 )) ] (𝑡𝑖+1 − 𝑡𝑖 ), 𝑖 = 0,1,2, . . , 𝑁
𝑚
Ahora, al empezar los cálculos con 𝑡0 = 0, 𝑣(𝑡0 ) = 0 y (𝑡𝑖+1 − 𝑡𝑖 ) = 2, constante, se obtienen los siguientes valores de velocidad
0.225
𝑡1 = 2, 𝑣 (2) = 𝑣(0) + [9.81 − 𝑣(0)2 ] (2) = 0.000 + 19.620 = 19.620
68.1
0.225
𝑡2 = 4, 𝑣(4) = 𝑣(2) + [9.81 − 𝑣(2)2 ] (2) = 19.620 + 17.076 = 36.696
68.1
0.225
𝑡3 = 6, 𝑣(6) = 𝑣(4) + [9.81 − 𝑣(4)2 ] (2) = 36.696 + 10.722 = 47.418
68.1
Las velocidades siguientes se estiman de modo similar.

Estos resultados pueden obtenerse de modo más fácil y en menor tiempo, mediante un programa informático, codificado en un software para
computo científico. A continuación, se muestra el programa VP, codificado en el software Octave, que permite aproximar las velocidades para una
sucesión de tiempos contenidos en [0, T].
Actividades para el estudiante: resuelva de modo analítico y numérico, los ejercicios 1.13 y 1.14, propuestos en la página 17 del libro guía.

Bibliografía

1. Chapra Steven, R. Canales. Métodos numéricos para ingenieros, edición 7 ma. Editorial Mc. Graw-Hill.
2. María Lucía Brito-Vallina, Isidro Alemán-Romero, Elena Fraga-Guerra, José Luís Para-García, Ruth Irene Arias-de Tapia. Papel de la modelación
matemática en la formación de los ingenieros. Ingeniería Mecánica vol.14 no.2 La Habana Mayo-ago. 2011
3. C. Henry Edward, David E. Penney. Ecuaciones diferenciales y problemas con valores en la frontera, cómputo y modelado. Edición 4 ta.,
editorial Prentice-Hall.

También podría gustarte

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy