UNIVERSIDAD NACIONAL DEL

CALLAO
FACULTAD DE INGENIERÍA
MECÁNICA Y DE ENERGÍA

“RESOLUCIÓN DE LA 2da PRÁCTICA DE

LABORATORIO”
CÁLCULO NUMÉRICO Y COMPUTACIONAL
BAYONA HUERTA JOSE
LEVITA AMBROSIO CRISTOPHER
SOTO TAIPE CRISTIAN
TARAZONA SALAZAR JAVIER
TRUJILLO TABOADA BRUNO
Callao, junio, 2018
PERÚ
 RESOLUCIÓN DE LA SEGUNDA PRÁCTICA DE LABORATORIO DE CÁLCULO NUMÉRICO Y
COMPUTACIONAL

1. Halle diferentes polinomios de Lagrange para aproximar V (25).

X (metros) 10 15 22 30 35 41
V (velocidad) 5 7 10 12 15 17

SOLUCIÓN:

Lo primero será abrir en MATLAB el programa “AITKEN”.

Modificamos el algoritmo:

X
(m
V
et
(V
Im
ro
el
ag
s)oc
en
id
(x
ad
1)
)

Corremos( ) y nos dará el siguiente resultado:

Donde:

Polinomio de LaGrange
xk fk x k −25 Grado 1 Grado 2 Grado 3 Grado 4 Grado 5
x 0=10 5 -15 11.0000 11.357142 P03=¿ 11.035 P04=¿10.796 P05=¿ 10.590
x 1=15 7 -10 11.2857142 P13=¿10.92 P14=¿ 10.637 P15=¿10.369 0
25 8
x 2=22 10 -3 P23=¿10.75 P24=¿ 10.34 P25=¿9.9422 0 0
0 6
x 3=30 12 5 9.00000000 7.7878787 0 0 0
x 4=35 15 10 11.6666666 0 0 0 0
x 5=41 17 16 0 0 0 0 0

Por lo tanto, los valores del polinomio de LaGrange para V (25) son:

P23 (25)=10.750000000000000
P13 (25)=10.928571428571429
P24 (25)=10.346153846153847
P03 (25)=11.035714285714288
P14 (25)=10.637362637362639
P25 (25)=9.942215679057785
P04 (25)=10.796703296703299
P15 (25)=10.369998422630003
P05 (25)=10.590233196345253

2. Resolver:
y ' =xycos ( xy ) −sin ⁡(x + y )
a) Usar Taylor de orden 2.
b) Usar el método de Runge-Kutta de cuarto orden.

SOLUCIÓN:

a) Usando Taylor de orden 2:

Daremos valores iniciales de la forma:

De tal manera que:

a=0 ; b=1; h=0.5 ; y ( 0 )=1

Colocamos en Matlab la parte izquierda de la ecuación como se muestra a
continuación:
 Abrimos nuestro programa “Taylor2”:

Teniendo en cuenta la forma de nuestro programa como se muestra en la figura,

colocaremos nuestros datos:

Finalmente, al correr el comando, Matlab nos dará los siguientes resultados:

b) Usando el método de Runge-Kutta de cuarto orden:

Daremos valores iniciales de la forma:

De tal manera que:

a=0 ; b=2; h=0.5 ; y ( 0 )=1

Colocamos en Matlab la parte izquierda de la ecuación como se muestra a
continuación:

Abrimos nuestro programa “Rungekutta”:

Teniendo en cuenta la forma de nuestro programa como se muestra en la figura,
colocaremos nuestros datos:
 Dándonos como resultado:

3. Resolver:
x '' −2 x ' + y=cos ⁡(t)−sin ⁡(t )
y ' ' + 8 x ' −x=t−1
a) Utilice ode23.
b) Utilice el método de Runge-Kutta para sistemas de cuarto orden.

SOLUCIÓN:

a) Usando “ode23”:
Lo primero será convertir esas ecuaciones diferenciales en un sistema de
ecuaciones:

z 1=x → z '1=x ' =z 2

z 2=x ' → z '2=x '' =2 z 2−z 3 +cos ( t )−sen(t)
z 3= y → z '3= y ' =z 4
 z4 = y ' → z '4 = y '' =z 1−8 z 2 +t−1

Obtendremos el siguiente sistema:

z '1=z 2
z '2=2 z2 −z3 +cos ( t ) −sen (t)
z '3=z 4
z '4 =z1−8 z 2 +t−1

Luego, guardamos el sistema en la función “Ed2” tal como se indica:

Una vez hecho esto nos dirigimos a la ventana de comandos y colocamos:

[ x , y ] =ode 23(' Ed 2' , [ 1 3 ] ,[1.5 ; 2 ; 2.5 ; 3])

Dónde:

 [ 13 ]: Es el dominio en el que estamos trabajando.

 [1.5; 2 ; 2.5 ; 3]: Son los valores que les dimos a las imágenes de
x ( 1 ) , x' ( 1 ) , y (1) y y ' (1) respectivamente
Finalmente, Matlab nos presentará el siguiente resultado:
 b) Usando el método de Runge-Kutta para sistemas de cuarto orden:
4. Presentar una aplicación de elemento finito en una dimensión.
4.1 Definición:

Hasta aquí hemos empleado métodos por diferencias finitas para resolver ecuaciones
diferenciales parciales. En estos métodos, el dominio de la solución se divide en una
malla con puntos discretos o nodos (figura 31.1b). Entonces, se aplica la EDP en cada
nodo, donde las derivadas parciales se reemplazan por diferencias finitas divididas.
Aunque tal aproximación por puntos es conceptualmente fácil de entender, tiene
varias desventajas. En particular, es difícil de aplicar a sistemas con una geometría
irregular, con condiciones de frontera no usuales o de composición heterogénea.

Figura 1: a) Empaque irregular; b) Sistema modelado con la técnica de diferencias finitas; c) Sistema
modelado por el método de elementos finitos.
Fuente:[ CITATION Cha07 \l 3082 ]
El método del elemento finito ofrece una alternativa que es más adecuada para tales
sistemas. A diferencia de las técnicas por diferencias finitas, la técnica del elemento
finito divide el dominio de la solución en regiones con formas sencillas o “elementos”
(figura 31.1c). Se puede desarrollar una solución aproximada de la EDP para cada uno
de estos elementos. La solución total se genera uniendo, o “ensamblando”, las
soluciones individuales, teniendo cuidado de asegurar la continuidad de las fronteras
entre los elementos. De este modo, la EDP se satisface por secciones.
Como se observa en la figura 31.1c, el uso de elementos, en lugar de una malla
rectangular, proporciona una mejor aproximación para sistemas con forma irregular.
Además, se pueden generar continuamente valores de las incógnitas a través de todo
el dominio de la solución en lugar de puntos aislados.

Dicho todo esto, en resume, podemos decir que el Método del Elemento Finito (MEF)
está pensado para ser usado en computadoras y permite resolver ecuaciones
diferenciales asociadas a un problema físico sobre geometrías complicadas. El MEF se
usa en el diseño y mejora de productos y aplicaciones industriales, así como en la
simulación de sistemas físicos y biológicos complejos. La variedad de problemas a los
que puede aplicarse ha crecido enormemente, siendo el requisito básico que
las ecuaciones constitutivas y ecuaciones de evolución temporal del problema sean
conocidas de antemano.

4.2 Problema de aplicación:

Resolver la ecuación diferencial que representa la forma unidimensional de la ecuación

de Poisson:

d2T
=−f (x) (1)
dx2
Para una barra de 10 cm con condiciones en la frontera de T (0 ; t)=40 y
T (10 ; t)=200 y una fuente de calor uniforme con f (x)=10. Utilice cuatro
elementos del mismo tamaño con longitud igual 2.5 cm.

Figura 2: a) Barra larga y delgada sujeta a condiciones de frontera fijas y una fuente de calor; b)
Representación del elemento finito de cuatro elementos de la misma longitud y cinco nodos
Fuente:[ CITATION Cha07 \l 3082 ]

4.3 Procesos

Los pasos que se siguen para resolver el problema serán los siguientes:

a) Discretización:
 Este paso consiste en dividir el dominio de la solución en elementos finitos.

En nuestro problema, una configuración simple para modelar el sistema consiste

en una serie de elementos de igual longitud (figura 2b). Así, el sistema se trata
como cuatro elementos de igual longitud y cinco nodos.

b) Ecuaciones de los elementos:

Este paso consiste en desarrollar ecuaciones para aproximar la solución de cada

elemento y consta de dos pasos. Primero, se debe elegir una función apropiada
con coeficientes desconocidos que aproximará la solución. Segundo, se evalúan los
coeficientes de modo que la función aproxime la solución de manera óptima.

Existen dos métodos para conocer estas ecuaciones, el primero se conoce como
“el método directo” y el segundo como “método de los residuos ponderados”. El
ultimo es el más usado en el área de ingeniería y para aplicarlo recurrimos a la
siguiente ecuación:

(2)

En nuestro problema, el término de la fuente de calor en el primer renglón de la

x 2−x
ecuación (2) se evalúa sustituyendo la ecuación N 1= , e integrando para
x 2−x 1
obtener:
2.5

∫ 10 2.5−x
2.5
dx=12.5
0

x−x 1
Del mismo modo, la ecuación N 2= se sustituye en el término de la fuente
x 2−x 1
de calor del segundo renglón de la ecuación (2), el cual también se integra para
obtener:
2.5

∫ 10 x−0
2.5
dx=12.5
0

Estos resultados junto con los valores de los otros parámetros se emplean para
sustituirse en la ecuación (2) y así obtener:

−dT
0.4 T 1−0.4 T 2= ( x ) +12.5
dx 1
y

dT
−0.4 T 1 +0.4 T 2= ( x ) +12.5
dx 2
c) Ensamble:

Una vez obtenidas las ecuaciones de elementos individuales, éstas deben unirse o
ensamblarse para caracterizar el comportamiento de todo el sistema. El proceso
de ensamble está regido por el concepto de continuidad. Es decir, las soluciones
de elementos contiguos se acoplan, de manera que los valores de las incógnitas (y
algunas veces las derivadas) en sus nodos comunes sean equivalentes. Así, la
solución total será continua[ CITATION Cha07 \l 3082 ].
Antes de que se ensamblen las ecuaciones del elemento se debe establecer un
esquema de numeración global que especifique la topología o el arreglo espacial
del sistema como en la tabla (1), esto define la conectividad de los elementos de la
malla. Debido a que este caso es unidimensional, el esquema de numeración es
tan predecible que parece trivial. Sin embargo, en problemas para dos y tres
dimensiones, tal esquema es el único medio para especificar qué nodos
pertenecen a que elementos.
Tabla 1:Topología del sistema para el esquema de segmentación del elemento finito de la figura (2b)

NUMERO DE NODOS
ELEMENTO
LOCAL GLOBAL
1 1
1
2 2
1 2
2
2 3
1 3
3
2 4
1 4
4
2 5

Una vez que se han especificado la topología, la ecuación del elemento (2) se
puede escribir para otros elementos usando las coordenadas del sistema. Después,
esta se agrega una por una para ensamblar la matriz de todo del sistema, el
proceso se ilustra a continuación:
d) Condiciones de frontera:
Observe que conforme se ensamblan las ecuaciones se cancelan las condiciones de
frontera interna. Así, el resultado final de {F } en la ecuación (2) e tiene
condiciones de frontera solo para el primero y el ultimo nodo. Ya que T 1y T 5están
dados, dichas condiciones de fronteras naturales en los extremos de la barra,
dT ( x 1)/dx y dT ( x 5)/dx representan incógnitas. Por lo tanto, las ecuaciones re
expresan como sigue:

e) Solución:

De lo anterior se obtiene:

dT T 2=173.75 T 3=245
( x ) =66
dx 1
 T 4=253.75 dT
( x )=−34
dx 5

f) Procesamiento posterior:
Finalmente, los resultados se representan gráficamente en la figura (3) se
muestran los resultados del método del elemento finito, junto con la solución
exacta. Observe que el cálculo del elemento finito capta la tendencia general de la
solución exacta y, además, da una coincidencia exacta en los nodos. Sin embargo,
existe una discrepancia en el interior de cada elemento debido a la naturaleza
lineal de las funciones de forma.

Figura 3: Resultados al aplicar el método del elemento finito a una barra calentada. También se
muestra la solución exacta
Fuente: [ CITATION Cha07 \l 3082 ]