Programación y aproximación a una ecuación de orden inferior por medio del

software Matlab
Julio Fernán Castañeda Ruiz
Escuela Colombiana de Ingeniería Julio Garavito
Elementos Finitos
21 de mayo del 2020
Abstract
En el siguiente informe se presenta la programación de una ecuación de orden
inferior, la cual se busca por medio del software Matlab la aproximación más cercana
a dicha ecuación, corroborando por medio de graficas que tan cercana y correcta fue
la programación para la ecuación de orden inferior planteada.

Introducción
A medida del tiempo todo problema físico se ha tendido a modelar mediante
planteamientos matemáticos siendo estos llevados al planteamiento de ecuaciones, para sí
tener una aproximación más cercana a cómo se comporta la problemática o problema a
analizar. En este documento plantea la solución de una ecuación de orden inferior
desarrollando principalmente el planteamiento de esta y luego llevada a la programación en
el software Matlab donde se una aproximación cercana a la ecuación original.
Marco teórico
A continuación, se hace el planteamiento como el desarrollo matemático de la ecuación
propuesta como estudio para el aprendizaje del software Matlab. La ecuación de orden
inferior caracterizada de la siguiente manera:
Ecuación 1
Las condiciones de problema serán:
Ecuación 2
Dicha ecuación será solucionada aplicando métodos como: residuos ponderados, Galerkin y
Polinomios de Lagrange. Principalmente se contempla el hecho de los coeficientes de la
ecuación como valores de A, B Y C.
Partir del método de los residuos ponderados tenemos:
Ecuación 3
Seguido de esto:
Ecuación 4
Para tener una mejor simplificación de la ecuación por lo cual el coeficiente de A será
quien divida a toda la ecuación centrándose principalmente en el lado izquierdo de la
ecuación. Entonces:
Ecuación 5
De esto resulta
Ecuación 6
Luego de esto tomando la ecuación 5 se procede a separar términos, de la siguiente manera
Ecuación 7
Se resolverá la primera integral usando el método de integración por partes, donde:
Ecuación 8
De resolver la integral por partes resulta
Ecuación 9
Donde C = 0, por condiciones del problema. Luego de esto se sigue trabajando en la zona
izquierda de la ecuación, teniendo en cuenta los valores obtenidos de la integral.
Ecuación 10
El termino:
Ecuación 11
Ya que por las características del problema este no presentara variación por sus
condiciones.
En esta parte se usa el método de Galerkin para lo cual se propone:
Ecuación 12
Sustiyendo
Ecuación 13
Otro planteamiento propuesto por el método es:
Ecuación 14
Entonces en cada valor de la ecuación se evaluará el propuesto anteriormente
Ecuación 15
Ahora se separan términos multiplicándolos
Ecuación 16
Seguidamente de agruparan términos semejantes en la ecuación
Ecuación 17
Ecuación 18
NL = N1
Ecuación 19
Ecuación 20
NL = N2
Ecuación 21
Ecuación 22
Según el planteamiento general de elementos finitos se debe llevar a la ecuación
fundamental de resortes compuesta por matrices, a continuación, se muestra la matriz
elemental:
Matriz ecuación 23
Donde de las sustituciones anteriores tendremos

Para k(1,1) = k(2,2)

Ecuación 24
Para k(1,2) = k(2,1)
Ecuación 25
Y para finalizar la matriz de fuerzas dada por
Ecuación 26
Ecuación 27
Hasta este momento se tiene el planteamiento matemático resulto según los métodos
propuestos para su desarrollo. Los polinomios de Lagrange son expuestos dentro el
programa de simulación en Matlab desarrollando los factores de N1 y N2.
Procedimiento
Posteriormente se mostrarán pantallazos del programa Matlab en línea donde se desarrolló
la programación de la aproximación de la ecuación para finalmente llegar a la comparación
de las dos gráficas, una mediante la solución analítica dada por Matlab y otra con ayuda de
la programación.
Inicialmente se cuenta con un programa main el cual está conformado por dos
funciones(subprogramas) las cuales se llamarán dentro del main en el momento de ser
necesario, esto con el fin de hacer más organizada y cómoda la programación además de
disminuir recursos utilizados por la máquina.
A continuación, se hará una explicación general de la composición y funcionalidad del
programa.

Imagen 1. Parte 1 Interfaz Usuario-Maquina.

Esta sección del programa presenta una interfaz a través de la cuál el usuario ingresará los
valores correspondientes a la asignación de los coeficientes para cada término de la
ecuación lineal.
Como la ecuación característica presenta tres términos se usa un ciclo FOR hasta tres para
ir pidiendo uno a uno los coeficientes de dichos términos. Cada término es evaluado en
caso de que por error el usuario ingrese un valor que no sea numérico, de ser así, el
programa automáticamente arroja un mensaje pidiendo al usuario ingresar nuevamente el
valor.
 Imagen 2. Parte 2 Interfaz Usuario-Máquina.

En caso de que alguno de los valores sea negativo, el programa convierte dicho valor a
positivo con el fin de seguir los lineamientos declarados matemáticamente para cada
integral.

Imagen 3. Parte 3 Interfaz Usuario-Máquina.

En la imagen anterior se presenta la misma lógica anterior.

 Imagen 4. Parte 4 Interfaz Usuario-Máquina.

En esa parte del código se solicita al usuario la cantidad de elementos que desea para la
solución aproximada de la ecuación. Dicho valor también es analizado para descartar que
no sea un valor numérico ni menor o igual a cero, ya que se necesita al menos un elemento
para poder aplicar el método.

Imagen 5. Declaración de variables iniciales.

En esta imagen se declaran las variables iniciales y se construye las coordenadas nodales
teniendo en cuenta punto inicial y punto final. Adicionalmente se declara la matriz de
elementos y se realiza el llenado según el número de elementos.
 Imagen 6. Inicialización del sistema en ceros y definición de variables.

Como se observa en la imagen 6, se construyen inicialmente con ceros tres matrices para el
sistema a solucionar. Seguidamente se definen los términos como N1 y N2, hallados
matemáticamente para nuestras integrales usadas según el método de aproximación.
 Imagen 7. Hallazgo de la matriz elemental y la matriz T.

En la imagen 7 se realiza el llenado de la matriz K elemental que servirá para hallazgo de la

matriz K Global y la matriz T.

Imagen 8. Llamado a la función LlenadoKglob.

Imagen 9. Función Llenado Kglob.

Esta función realiza el llenado de la matriz Kglobal a través de un ciclo FOR, sin pasar por
alto que se debe inicializar con ceros antes.

Imagen 10. Penalización del primer y último elemento de las matrices.

En la parte inicial de las variables se tomó en cuenta el valor de melf, este se usa para hacer
penalización a las matrices T y la matriz de Kglobal. Luego se grafica la sección del
llenado de kglob como también la solución y grafica de esta.
 Imagen 11. Llamado a función y graficación de las soluciones.

Imagen 12. Función para solucionar la ecuación a través de Matlab.

Esta función realiza la solución de la ecuación planteada con ayuda del software
tomándolos valores del usuario, se evalúa y se vectoriza con el comando para su solución.
Y volviendo al final de código main imagen 11, se grafica para compararse con la
aproximación hallada dada por programación.
Conclusiones
Viendo la comparación de las dos graficas resulta de una buena aproximación con valores
de números de elementos como 1,2,3,7,11 para así, poder rectificar que la programación
desarrollada tuvo una buena aproximación mediante el software Matlab.
Cabe resaltar que el lenguaje de Matlab cuenta con gran parte de comandos y de esto
resulta el hecho de conocer estos comandos y entenderlos hará que la programación sea
más eficiente rápida o en caso contrario se cuente con dificultades a la hora de llevar a cabo
el modelo de una problemática propuesta en una aplicación ingenieril. Sin embargo, aplicar
programas al modelo de aplicación físicas junto con programación dará un resultado más
eficaz a la hora de saber si un elemento, pieza, equipo o sistema podrá aceptarse o porque
no mejorarse con ayuda del método de elementos finitos.