UNIVERSIDAD PRIVADA ANTENOR ORREGO
FACULTAD DE INGENIERÍA
ESCUELA DE INGENIERÍA CIVIL
APLICACIÓN DE LOS MÉTODOS NUMÉRICOS EN
INGENIERÍA CIVIL-TERMODINÁMICA
ALUMNOS: -Asmat Castillo, Jean (217619)
-Cachay Mejia, Mayra (220763)
-Centeno Alquizar, Kevin (217784)
-Cordova Infante, Eduardo (217899)
-Cusma Flores, Sayuri (201265)
-Macso Monteza, Solange (217116)
-Otiniano Oliva, Jhuliana (217660)
-Zapata Niño, Pedro (216201)
DOCENTE: ANGULO VILCA, PERCY
ASIGNATURA: MATEMÁTICA IV
CICLO: IV
TRUJILLO – 2020
1
ÍNDICE
1. INTRODUCCIÓN………………………………………………………..
03
2. PROBLEMÁTICA……………………………………………………….
04
3. OBJETIVOS……………………………………………………………..
05
4. FUNDAMENTO TEÓRICO……………………………………………. 06
4.1. CONCEPTOS Y DEFINICIONES BÁSICAS……………....... 06
4.2. MARCO TEÓRICO……………………………………………...
07
5. SOLUCIÓN DEL PROBLEMA………………………………………...
13
5.1. TOMA DE DATOS…………………………………………….... 13
5.2. FÓRMULAS MATEMÁTICAS…………………………………..
13
5.3. DESCRIPCIÓN Y FORMALIZACIÓN DEL PROBLEMA…… 14
5.4. RESULTADOS………………………………………………….. 15
5.5. CÓDIGO EN MATLAB…………………………………………..
18
6. CONCLUSIONES……………………………………………………….
20
7. RECOMENDACIONES………………………………………………... 21
8. BIBLIOGRAFÍA…………………………………………………………. 22
2
1. INTRODUCCIÓN
Desde hace varios años se conoce que el rol del ingeniero civil está basado
principalmente en el conocimiento de los fenómenos y mecanismos
relacionados a la naturaleza y la realidad; dando un servicio a la comunidad, es
por ello que utiliza su intelecto y todo lo que ha aprendido en su etapa de
formación profesional. Hoy en día, debido al avance de la tecnología, se utilizan
diversos métodos para dar solución a los problemas relacionados a la sociedad
y que estos involucren el tema constructivo; ya que esto va a generar un fuerte
impacto en la vida del ser humano.
Uno de los temas más vistos en la carrera de ingeniería civil son los métodos
numéricos, debido a que estos están relacionados al ámbito constructivo; es
decir, están vinculados a los fenómenos físicos de la ingeniería civil. Como
estudiantes de esta carrera nos sentimos dichosos de poder abordar los
distintos modelos matemáticos que se verán reflejados a lo largo de nuestra
carrera profesional, con el único objetivo de dar una solución a los diversos
problemas o situaciones que puedan ocurrir en nuestra sociedad.
La fábrica ORLANDO S.A.C tiene como objetivo principal la producción de
material constructivo destinado a las obras públicas y privadas, manteniendo
un control en la temperatura de los materiales que fabrica, debido a la
3
presencia de altas temperaturas en la zona en la que ubica. Además de estar
comprometida con el bienestar de sus consumidores y sus trabajadores.
En el presente trabajo de investigación, analizaremos una situación relacionada
a la termodinámica en la que emplearemos un método numérico para hallar
una solución aproximada, de acuerdo a lo aprendido en clase.
2. PROBLEMÁTICA
Es de suma importancia a nivel mundial llevar el control y análisis respectivo
durante una instalación de tuberías, dentro de ello los materiales o
substancias involucradas en los mismos. Según estadísticas de la MHIDAS
(Major Hazard Incident Data Service), se han podido registrar los accidentes
según el lugar o actividad en dónde estos se han producido, hemos visto
que entre los distintos orígenes donde se trabaja con tuberías, la mayoría
de accidentes han ocurrido en plantas de proceso que se dan por
combustión de substancias presentes en tanques de proceso (36.7%).
Es de importancia comentar el gran despliegue de tuberías, codos y
conexiones que se da en la mayoría de las plantas de proceso,
representando un riesgo potencial de fuga de sustancias inflamables.
Debido a esto, la fábrica ORLANDO S.A.C ha decidido llevar un control de
calidad de sus productos, para así en caso de encontrar algún desperfecto
poder modificar y mejorarlo; garantizando seguridad tanto a sus
consumidores como a sus trabajadores, para poder reducir perdidas
económicas y mortales.
Para dicho control se estima que si T(t) representa la temperatura del
cuerpo en el tiempo t, 𝑇𝑚 representa la temperatura del medio ambiente.
4
Entonces, si el cuerpo es un tubo de acero calentado a 400°C, se requiere
conocer su temperatura al cabo de 5 minutos. El tubo se encuentra en una
cámara a temperatura de 25°C. Para este caso, la constante de
proporcionalidad es -0.213.
De acuerdo al enunciado y a los datos mencionados, el problema a analizar
es el siguiente:
¿Será posible conocer la temperatura aproximada a la que llegará el
tubo de acero al ser calentado a 400°C durante 5 minutos?
3. OBJETIVOS
3.1 Objetivo general
Conocer la temperatura aproximada a la que llegará el tubo de
acero al ser calentado a 400°C durante 5 minutos.
3.2 Objetivos específicos
Calcular el control de calidad de los productos.
Usar el método de Euler para determinar la temperatura
aproximada a la que podría llegar el tubo de acero en 5
minutos al ser calentado a 400°C.
5
4. FUNDAMENTO TEÓRICO
4.1 CONCEPTOS Y DEFINICIONES BÁSICAS
MATEMÁTICA: Es una ciencia en la cual se emplea símbolos para crear una
teoría exacta basada en definiciones.
Esta ciencia impulsa a una persona a pensar de forma lógica y, por ende, a
desarrollar capacidades, solucionar problemas y tomar decisiones.
INGENIERÍA: La ingeniería es la aplicación de la ciencia en la conversión
óptima de los recursos naturales para uso de la humanidad.
Asociado a la ingeniería existe un monumental cuerpo de conocimientos
especiales, la preparación para la práctica profesional implica una amplia
capacitación en la aplicación de ese conocimiento.
INGENIERÍA CIVIL: La ingeniería civil es la rama de la ingeniería donde se
aplica métodos científicos, mayormente se aplican los conocimientos de la
física y la matemática.
Su objetivo es ejecutar edificaciones que sean útiles para la sociedad.
6
TEMPERATURA: La temperatura es la cantidad de energía cinética que hay
en un cuerpo.
A mayor temperatura tendremos más sensación de calor y a menor
temperatura la sensación de frio aumenta.
TERMODINÁMICA: Es una disciplina que se encarga de estudiar las acciones
mecánicas del calor y otras formas semejantes de energía, con esta disciplina
se analizan los efectos que producen los cambios de magnitudes tales como: la
temperatura, densidad, la presión, masa, volumen, etc.
Los procesos termodinámicos están determinados por tres leyes básicas.
La primera ley permite comprender cómo se conserva la energía.
La segunda ley es usada para conocer las condiciones necesarias para que
ocurra la transferencia de energía.
La tercera ley sirve para conocer el comportamiento de los sistemas en
equilibrio.
4.2. MARCO TEÓRICO
RECTA
Líneas compuestas por puntos que suceden de manera indefinida, no sufren
curvaturas, indican una única dirección y dos sentidos contrarios.
RECTA TANGENTE
Una recta se le puede decir que es tangente a una recta cuando la toca en un
solo punto y no la corta.
Recta Normal Recta Tangente
7
f (x)
f (a)
a x
DERIVADAS
La derivada es uno de los conceptos más importante en matemáticas. La
derivada es el resultado de un límite y representa la pendiente de la recta
tangente a la gráfica de la función en un punto.
f ( a+ h )−f (a)
f ' ( ∂ ) =lim
h→0 h
ECUACIÓN DIFERENCIAL
Se puede definir como la relación entre una variable x y la función buscada y(x)
es decir, es una ecuación matemática que vincula una función con sus
derivadas Solo las ecuaciones diferenciales más simples se pueden solucionar
por medio de fórmulas explícitas, sin embargo, se pueden establecer algunas
propiedades de las soluciones de una cierta ecuación diferencial sin encontrar
su forma exacta.
dy
=f ( x , y)
dx
ECUACIONES DIFERENCIALES ORDINARIAS
8
Es una ecuación que involucra a una variable independiente x, una función y(x)
y una y varias derivadas, se les denominar ordinarias porque aparecen
derivadas ordinarias de una a más variables dependientes respecto a una
única variable independiente.
dy
{
[ L ] = dt
=f (t , y )
y ( t 0 )= y 0
ECUACIÓN EN DERIVADAS PARCIALES
Una ecuación en derivadas parciales es una ecuación diferencial que contiene
una función multi variable y sus derivadas parciales. Lo que nos permite hablar
esta ecuación es una cantidad física como, por ejemplo, temperatura,
concentración de un contaminante, intensidad de una señal acústica,
deformación de una estructura, etc.
∂y ∂y
−k =0
∂x ∂x
MÉTODO DE EULER
El método de Euler consiste en encontrar la solución de ecuaciones
diferenciales de primer orden por medio de la aproximación a partir de un valor
dado al inicio.
Teniendo como datos iniciales:
dy
=f ( x , y ) , para , x 0 ≤ x ≤ x f
dx
y ( x 0 )= y 0
9
El desarrollo del método consistirá en lo siguiente:
Se hallará el h (ancho) ; dividiendo los intervalos que van desde x 0 hasta
x f entre la cantidad de sub intervalos ( n ) :
x f −x 0
h=
n
Obteniéndose un conjunto de n + 1 puntos del intervalo [ x ¿ ¿ 0 , x f ]¿
La condición inicial y ( x 0 )= y 0 representa el punto P0=(x 0 , y 0 ) por donde
pasa la curva solución de la ecuación inicial del problema, el cual lo
podemos denotar como F ( x )= y .
Teniendo el punto P0 se procede a realizar la primera derivada de F ( x )
dy
'
en el punto dado; obteniendo que F = dx | va a ser igual a f (x , y ).
P0
0 0
Con esto, se procede a hacer el siguiente gráfico. Se ha trazado una
recta que pasa por el punto ( x 0 , y 0 ) con una pendiente f (x 0 , y 0 ), además
de la curva solución.
y
( x 1 , y ( x 1))
Error
y1
Recta Secante
F ( x1 ) ( x1 , y1)
f (x 0 , y 0 )
10
y0
( x0 , y0 )
x0 x1 x
Según la gráfica y por fórmula de la pendiente, se obtiene lo siguiente:
y 1− y 0
=f (x0 , y 0 )
x 1−x 0
Despejando y 1:
y 1= y 0 + ( x 1−x 0 ) f ( x 0 , y 0 )
y 1= y 0 +h f ( x 0 , y 0 )
y n+ 1= y n +h f ( xn , y n)
Despejando x 1:
x 1−x 0
h=
n
x 1=x 0+ h∗n
x n+1=x n +h∗n o x n+1= y n+ h
LEY DE ENFRIAMIENTO DE NEWTON
La ley de enfriamiento de Newton hace referencia a que la pérdida de calor de
un cuerpo es proporcional a la diferencia de la temperatura del cuerpo entre la
del medio ambiente en el que se encuentra. La fórmula se expresa de la
siguiente forma:
11
dQ
=α S(T −T m )
dt
Donde:
α : Coeficiente de intercambio de calor
S: Área superficial del cuerpo expuesto al medio ambiente
Gracias al modelo matemático de Newton es posible saber los cambios de
temperatura con mucha precisión cuando la diferencia de temperatura entre
dos cuerpos no es grande. Para ello, se analizará la siguiente ecuación
diferencial:
dT
=α(T −T m )
dt
dT
=k (T −T m)
dt
Donde:
T: Temperatura del objeto al cual se analizará
T m: Temperatura del medio ambiente en el que se encuentra el objeto
t: Tiempo en el que ocurre el cambio de temperatura.
k: Constante de proporcionalidad.
12
5. SOLUCIÓN DEL PROBLEMA
5.1. TOMA DE DATOS
Ti 400 ° C
Tm 25 ° C
Tf ?
t 5 min
k −0.213
5.2. FÓRMULAS MATEMÁTICAS
REALIZACIÓN DEL MODELO MATEMÁTICO
Este método se puede realizar de la siguiente manera, de acuerdo al
establecimiento de la Ley de Newton enunciada:
13
La diferencia entre la temperatura del cuerpo y el medio circundante
T −T m
Como la disminución de temperatura es proporcional a la diferencia, se
tiene:
k (T −T m )
Lo anterior, representa al sistema en movimiento, por lo que existe un
motivo de cambio de la temperatura del cuerpo con respecto al tiempo
de enfriamiento:
ⅆT
ⅆt
Y el modelo matemático es:
ⅆT
=k (T −T m)
ⅆt
En el método, el intervalo se divide en sub intervalos, con un número de
pasos (n) que lo definimos de acuerdo con el contexto del problema.
El tamaño del paso lo definimos como:
b−a
h=
n
Aplicar el método de Euler, el cual tiene propiedades similares a las
aproximaciones sucesivas.
y i+1 = y i+ hf (x , y )
5.3. DESCRIPCION Y FORMALIZACION DEL PROBLEMA
14
Si T representa la temperatura corporal en la era t, 𝑇 representa la temperatura
ambiental. Entonces, si el cuerpo humano es un tubo de acero calentado a
400°C, es preciso conocer su temperatura al cabo de 5 min. El tubo está en
una cámara a temperatura de 25°C. Para esta situación, podemos trabajar con
la constante de proporcionalidad -0.213.
Organización de la información:
Se le llama t al tiempo, en minutos (min)
T m, es la temperatura de la cámara en °C
T, es la temperatura del tubo en °C
Sustituyendo en la fórmula de Euler:
T i+1=T i+ h[ k ( T i−T m ) ]
Generamos una tabla
T 0=T ( 0 )=400
T m=25
0 ≤ t ≤5
k =−0.213
h=0.5
5.4 RESULTADOS
Aplicandola f ó rmula del paso 3 :
ITERACION 01:
T ( 0 )=400
ITERACION 02:
15
T 1=T 0 + h [ k ( T 0 −T m ) ]
T ( 1 )=T 0 +0.5 [ −0.213(400−25) ]
T ( 1 )=360.0625
ITERACION 03:
T 2=T 1 +h [ k ( T 1−T m ) ]
T ( 2 ) =360.0625+ 0.5 [ −0.213(360.0625−25) ]
T ( 2 ) =324.3783
ITERACION 04:
T 3=T 2+ h [ k ( T 2 −T m ) ]
T ( 3 )=324.3783+0.5 [ −0.213(324.3783−25) ]
T ( 3 )=292.4945
ITERACION 05:
T 4=T 3 +h [ k ( T 3−T m ) ]
T ( 4 ) =292.4946+0.5 [ −0.213(292.4946−25) ]
T ( 4 ) =264.0063
ITERACION 06:
T 5=T 4 +h [ k ( T 4 −T m ) ]
T ( 5 ) =264.0064+0.5 [ −0.213( 264.0064−25) ]
T ( 5 ) =238.5521
16
ITERACION 07:
T 6=T 5+ h [ k ( T 5−T m ) ]
T ( 6 )=238.5522+0.5 [ −0.213(238.5522−25) ]
T ( 6 )=215.8088
ITERACION 08:
T 7=T 6+ h [ k ( T 6−T m ) ]
T ( 7 )=215.8089+0.5 [ −0.213(215.8089−25) ]
T ( 7 )=195.4877
ITERACION 09:
T 8=T 7+ h [ k ( T 7−T m ) ]
T ( 8 )=195.4878+0.5 [ −0.213(195.4878−25)]
T ( 8 )=177.3308
ITERACION 10:
T 9=T 8+ h [ k ( T 8−T m ) ]
T ( 9 )=177.3308+0.5 [ −0.213(177.3308−25) ]
T ( 9 )=161.1076
ITERACION 11:
T 10=T 9 + h [ k ( T 9 −T m ) ]
T ( 10 ) =161.1076+0.5 [ −0.213(161.1076−25) ]
T ( 10 ) =146.6121
I h T(final)
01 0 400
02 0.5 360.0625
03 1 324.3783
04 1.5 292.4945
05 2 264.0063
17
06 2.5 238.5521
07 3 215.8088
08 3.5 195.4877
09 4 177.3308
10 4.5 161.1076
11 5 146.6121
Cuando t=5 min. , entonces , latemperatura del tuboaproximada es 146.6121 ° C .
5.5 CODIGO EN MATLAB
function [x,y]= Euler(f,a,b,h,y0)
%Metodo de Euler para problemas de valor incial
%n: numero de iteraciones
%h: paso
%a,b: intervalo de x
%y0: valor incial
n=(b-a)/h;
x=[a:h:b]';
y=zeros(n+1,1);
y(1)=y0;
for k=1:n
y(k+1)=y(k)+h*feval(f,x(k),y(k));
end
end
Aplicándolo a nuestro problema:
>> f=@(x,y) -0.213*(y-25)
f=
function_handle with value:
@(x,y)-0.213*(y-25)
>> [x,y]= Euler(f,0,5,0.5,400)
18
x= 0
0.5000
1.0000
1.5000
2.0000
2.5000
3.0000
3.5000
4.0000
4.5000
5.0000
y=
400.0000
360.0625
324.3783
292.4945
264.0064
238.5522
215.8089
195.4877
177.3308
161.1076
146.6121
>> plot(x,y)
xlabel('Eje de abscisas')
ylabel('Eje de las ordenadas')
title('GRAFICA DE EULER')
hold on
19
6. CONCLUSIONES
Considerar que para solucionar cada problema de los procedimientos
numéricos se necesita tener orden pues la proporción de datos son
demasiados, además es necesario tener los programas para solucionar
cada procedimiento.
Concluimos que la Termodinámica es el análisis de las características
de sistemas de monumental escala en equilibrio en las que la
temperatura es una variable fundamental.
Los métodos numéricos, permite solucionar diversos problemas físicos
en forma eficiente.
Estos se pueden aplicar en distintos campos, para encontrar resultados
aproximados a sistemas complejos utilizando sólo las operaciones
matemáticas más simples.
20
7. RECOMENDACIONES
Identificar los datos para así incluirlos adecuadamente al método
elegido.
Generar una tabla con las variables de iteración para así obtener una
mayor organización, evitando errores.
Comprobar si el sistema planteado tiene una solución coherente para la
simulación del mismo para un análisis preciso de resultados.
Escoger un “n” que ayude a obtener una solución coherente al problema.
21
8. BIBLIOGRAFIA
Agrazal, C., Arauz, A., Bernal, J., Bonilla, N., Marciaga, F., Miranda, E.,
Vivar, L. (06 de diciembre de 2011). Método de Euler. Panamá:
Issuu.com. Recuperado de
https://issuu.com/luiseduardovivar/docs/monografia_metodo_de_euler
Muoz, E. (14 de junio de 2017). Ecuación diferencial (Ley de
enfriamiento de Newton). Eduardomuozfisica.
https://eduardomuozfisica.wordpress.com/2017/06/14/ecuacion-
diferencial-ley-de-enfriamiento-de-newton/
Pérez, J., Merino, M. (2010). Definición de Termodinámica.
Definicion.DE. https://definicion.de/termodinamica/
Vilas, I. (7 de octubre de 2016). Definición de la temperatura de diseño.
Arveng: training & engineering. https://arvengtraining.com/definicion-la-
temperatura-diseno
López, J. (junio de 2017). Análisis cuantitativo de riesgos de tuberías de
transporte de substancias peligrosas.
https://upcommons.upc.edu/bitstream/handle/2117/110231/TGF_Javier_
Lopez_Lopez.pdf
22
Angulo, P. (Diciembre de 2020). Método de Euler. [diapositiva]. Trujillo,
Perú, 22 diapositivas.
23
