APLICACIÓN DE MÉTODOS NUMÉRICOS EN ANÁLISIS DE ELEMENTOS
FINITOS (AEF) APLICADA EN ESTRUCTURAS
Esaú Álvarez Y. (1), Edu Ruiz P. (1)
• Universidad Tecnológica del Perú – Sede: Chimbote.
Código: U22245644, U22226630
• RESUMEN
En el presente trabajo se dará a conocer como aplicamos los métodos numéricos a través de
análisis de elementos finitos, implementadas en estructuras. A continuación, se dará un
breve resumen.
Este proyecto consiste en dividir una estructura compleja en partes más pequeñas
llamadas elementos finitos, lo que permite analizar su comportamiento de forma
aproximada pero eficiente. Dividiremos en 3 elementos la estructura para llegar a
tener como resultados sus desplazamientos, reacciones y esfuerzos internos. Se
utiliza este enfoque principalmente para estudiar cómo se deforman o fallan las
estructuras cuando se les aplican cargas, temperaturas u otras condiciones. Al
emplear métodos numéricos como “Runge-Kutta” o diferencias finitas, es posible
resolver problemas que serían muy difíciles de abordar con métodos exactos.
Además, el AEF resulta útil en el diseño y evaluación de estructuras en ingeniería
civil, mecánica o aeroespacial, ya que ayuda a optimizar materiales, reducir costos
y garantizar la seguridad de los sistemas analizados.
• INTRODUCCIÓN
El Análisis de Elementos Finitos (AEF) surge como respuesta a la necesidad de
resolver problemas estructurales complejos que no podían abordarse fácilmente con
métodos analíticos tradicionales. A medida que las estructuras en ingeniería se
volvieron más sofisticadas, fue necesario utilizar métodos numéricos para predecir
con mayor precisión su comportamiento ante cargas, deformaciones y otros
factores. El AEF permite dividir una estructura en elementos más simples,
facilitando el análisis computacional.
a. Descripción del proyecto.
Plan de trabajo propuesto
Para llevar a cabo el proyecto sobre la aplicación de métodos numéricos en el
Análisis de Elementos Finitos (AEF) aplicado en estructuras, se propuso el
siguiente plan de trabajo dividido en cinco etapas principales:
• Investigación teórica (1 semana): Revisión de fundamentos sobre métodos
numéricos y conceptos básicos del AEF, así como su aplicación en ingeniería
estructural.
• Selección del caso de estudio (1 semana): Elección de una estructura simple
(viga, pórtico o placa) para aplicar el análisis numérico.
• Modelado y discretización (1 semana): Representación del modelo estructural
en elementos finitos, asignación de propiedades y condiciones de frontera.
• Aplicación de métodos numéricos (2 semanas): Implementación de técnicas
como diferencias finitas o métodos iterativos para resolver el sistema de
ecuaciones generado.
• Análisis de resultados y conclusiones (1 semana): Evaluación del
comportamiento estructural, comparación con resultados teóricos o software, y
redacción del informe final.
Duración total aproximada: 6 semanas
b. Objetivos.
• Objetivos principales:
• Modelar y resolver el comportamiento estructural de sistemas
complejos.
• Optimizar el diseño estructural mediante simulaciones confiables.
• Objetivos secundarios:
• Reducir el tiempo y costo en pruebas físicas.
• Validar y comparar soluciones analíticas y experimentales.
c. Alcances y limitaciones.
Alcances:
El proyecto permite comprender y aplicar los fundamentos del Análisis de
Elementos Finitos (AEF) utilizando métodos numéricos en estructuras simples. Se
logrará modelar un caso estructural básico, discretizarlo y resolverlo mediante
técnicas numéricas, lo cual brinda una base sólida para estudios más avanzados.
Además, se podrá comparar el comportamiento teórico con resultados aproximados
obtenidos numéricamente, fortaleciendo la comprensión del análisis estructural.
Limitaciones:
Debido al tiempo limitado del proyecto y a la disponibilidad de recursos, el análisis
se enfocará únicamente en estructuras de baja complejidad. No se utilizarán
herramientas de simulación avanzada ni software especializado de alto costo, lo que
limita la precisión de los resultados. La experiencia previa en programación y
métodos numéricos también puede restringir el nivel de detalle alcanzado en la
modelación.
• FUNDAMENTO TEÓRICO
• Métodos numéricos en la ingeniería estructural
Los métodos numéricos permiten resolver problemas complejos de la
ingeniería estructural que no pueden abordarse con soluciones analíticas
exactas. Estos métodos se basan en aproximaciones sucesivas para encontrar
soluciones aceptables de ecuaciones diferenciales que modelan el
comportamiento físico de estructuras sometidas a cargas. En el contexto del
AEF, los métodos numéricos facilitan el cálculo de desplazamientos,
esfuerzos y deformaciones en estructuras con geometrías irregulares,
materiales no homogéneos o condiciones de carga complejas, que serían
imposibles de abordar con técnicas clásicas.
• Fundamento del análisis por elementos finitos (AEF)
El AEF es una técnica numérica que permite analizar el comportamiento
estructural dividiendo el dominio continuo de una estructura en un número
finito de subdominios llamados elementos. Cada elemento está conectado a
otros mediante nodos, y dentro de ellos se asume una forma simple para el
campo de desplazamientos, que se representa con funciones de forma. Este
método transforma el problema continuo en un sistema de ecuaciones
algebraicas que puede resolverse numéricamente, proporcionando
información detallada sobre la respuesta estructural bajo diferentes
condiciones.
• Formulación de la matriz de rigidez
La matriz de rigidez es el núcleo del análisis en AEF, ya que relaciona
directamente las fuerzas aplicadas con los desplazamientos en la estructura.
Se construye a partir de la formulación variacional o del principio de los
trabajos virtuales, usando funciones de forma e integración numérica sobre
cada elemento. Luego, estas matrices elementales se ensamblan en una
matriz global del sistema. Este procedimiento requiere el uso de métodos
numéricos, especialmente para la integración, ya que los elementos pueden
tener geometría compleja o propiedades variables.
• Resolución del sistema de ecuaciones
Una vez ensamblada la matriz global del sistema y aplicadas las condiciones
de frontera, se obtiene un sistema de ecuaciones lineales:
K = matriz de rigidez global
u = vector de desplazamientos nodales desconocidos
F = vector de fuerzas nodales conocidas
Este sistema puede ser muy grande, por lo que se requieren métodos
numéricos eficientes para resolverlo, como el método de eliminación de
Gauss, factorización LU o métodos iterativos como el de Jacobi o Gauss-
Seidel, dependiendo del tamaño y naturaleza del problema.
• Método de eliminación gussiana
El método de eliminación gaussiana es una técnica numérica utilizada para
resolver sistemas de ecuaciones lineales. Consiste en transformar el sistema
original en una forma escalonada (triangular superior) mediante operaciones
elementales sobre las filas de la matriz aumentada. Una vez en esa forma, se
aplica sustitución hacia atrás para encontrar las incógnitas. Es un método
directo y fundamental en álgebra lineal. Se utiliza ampliamente en
ingeniería para resolver matrices provenientes del análisis estructural.
• MÉTODO Y MATERIALES
Para dar solución a los errores de una simulación de comportamiento de estructuras
bajo diferentes condiciones, así como también modelar y analizar estructuras
complejas sometidas a cargas estáticas y dinámicas. Se implementarán métodos
numéricos como el Método de eliminación gaussiana qué sirve para resolver
sistemas de ecuaciones lineales, así como también, la integración de gauss para
integrar numéricamente funciones en un dominio determinado.
Además, se aplicará el método de Newton para optimizar funciones utilizando la
matriz de Hessian. Estas técnicas permitirán tener el mínimo porcentaje de error al
momento de calcular la simulación de las estructuras de acuerdo a las condiciones
que están sometidas, siendo cargas estáticas o dinámicas.
PROCEDIMIENTO EXPERIMENTAL:
Se analizó una barra sometida a compresión axial con el objetivo de calcular su
matriz de rigidez. Se definieron las propiedades físicas del elemento: área de
sección transversal, módulo de elasticidad y longitud. Luego, se aplicó la
formulación de rigidez axial (k = AE/L) para establecer la matriz de rigidez 2×2 del
sistema. Estos cálculos se realizaron en Excel, permitiendo visualizar los
coeficientes que relacionan las fuerzas internas con los desplazamientos en los
nodos. Se completó la simulación ensamblando la matriz global del sistema y
verificando su simetría. Este procedimiento permite comprender cómo los
parámetros físicos de una barra afectan su rigidez estructural, base esencial del
análisis por elementos finitos.
Figura 1. Barra de acero para someterlo a diferentes propiedades
físicas [1].
• RESULTADOS
En primer lugar, hacemos el recopilado de datos que nos proporciona la imagen:
Elementos Longitud (m) Área (m2) E(kg/m2) A*E/L
1 3 0.005 21000000000 35000000
2 4 0.002 21000000000 10500000
3 5 0.003 21000000000 12600000
Luego generamos más matrices de rigidez locales con los resultados de la tabla anterior:
Matriz de rigidez
local 1 2 3 4
35000000 -35000000 1
-35000000 35000000 2
K1
3
4
Matriz de rigidez
local 1 2 3 4
1
10500000 -10500000 2
K2
-10500000 10500000 3
4
Matriz de rigidez local 1 2 3 4
1
2
K3
12600000 -12600000 3
-12600000 12600000 4
Sumamos las matrices locales para generar una global:
Matriz de rigidez
global 1 2 3 4
35000000 -35000000 0 0 1
-35000000 45500000 -10500000 0 2
K
0 -10500000 23100000 -12600000 3
0 0 -12600000 12600000 4
Por otro lado, recopilamos los siguientes datos de la imagen:
Vectores desplazamiento Vector Fuerzas
0 1 R1 1
d2 2 800 2
d= F=
d3 3 -900 3
0 4 R4 4
Para hallar los vectores de desplazamiento haremos uso de la sub matriz de rigidez
global y los vectores de fuerza:
Matriz de rigidez
global 1 2 3 4
35000000 -35000000 0 0 1
-35000000 45500000 -10500000 0 2
K
0 -10500000 23100000 -12600000 3
0 0 -12600000 12600000 4
45500000 -10500000 d2 800
x =
-10500000 23100000 d3 -900
A = Sub matriz de rigidez B = Vectores de
global desplazamiento C = Vectores de fuerza
2x2 2x1 2x1
Mediante la siguiente formula se halla los vectores de desplazamiento:
A^-1*A*B=A^-1*C
B=A^-1*C
Vectores de desplazamiento:
d2 9.59821E-06
d3 -3.45982E-05
Con los datos obtenidos completamos los datos que nos faltó (R1 y R4):
0 1 -335.9375 1
9.59821E-06 2 800 2
d F
-3.45982E-05 3 -900 3
0 4 435.9375 4
Resumen de resultados recopilados
Desplazamientos
d2 = 0.009598214 (mm)
-
d3 = 0.034598214 (mm)
Reacciones internas
R1 = -335.9375 (kg)
R4 = 435.9375 (kg)
Esfuerzos Internos
Elemento 1 335.9375 (kg) Tracción
Elemento 2 -464.0625 (kg) Compresión
Elemento 3 435.9375 (kg) Tracción
Representación grafica de los esfuerzos internos:
• CONCLUSIONES
• Al desarrollar nuestro proyecto, comprobamos que el método de elementos
finitos permite analizar estructuras con formas irregulares, distintos
materiales y cargas complejas, algo que sería muy difícil o imposible con
métodos tradicionales.
• Aprendimos que una malla más fina suele ofrecer mayor precisión, aunque
requiere más recursos. Además, al elegir adecuadamente el tipo de elemento
(barra, placa o sólido), logramos que el modelo represente de forma más
realista el comportamiento de la estructura.
• Mientras que el análisis estático considera todo en reposo, el análisis
dinámico toma en cuenta la respuesta de la estructura a lo largo del tiempo.
Para ello, se utilizan métodos especiales como el de Newmark, que integran
las ecuaciones en el tiempo. Aunque no se aplicó en el proyecto por su
complejidad, es un enfoque más realista y muy útil en situaciones como
terremotos o vibraciones.
• Una de las principales utilidades de nuestro proyecto fue la posibilidad de
aplicar distintos tipos de carga (constantes, variables y dinámicas) bajo una
cinemática básica, y observar cómo respondía la estructura. Esto es muy
relevante en la práctica, ya que en la realidad las cargas suelen cambiar con
el tiempo o presentarse de forma repentina, y no son siempre estáticas.
• Otro aspecto destacado es que el método no solo sirve para analizar
estructuras, sino también para optimizarlas. Con estas herramientas, es
posible realizar ajustes en el diseño para hacerlo más liviano o económico
sin comprometer su resistencia, lo cual resulta muy útil en campos como la
ingeniería civil, mecánica y aeroespacial.
• REFERENCIA BIBLIOGRÁFICA
• Logan, Daryl L.
A First Course in the Finite Element Method (6th Edition).
Cengage Learning, 2016.
• Zienkiewicz, Olek C., y Taylor, Robert L.
The Finite Element Method: Its Basis and Fundamentals (7th Edition).
Elsevier, 2013.
• Cook, Robert D., Malkus, David S., Plesha, Michael E., Witt, Robert J.
Concepts and Applications of Finite Element Analysis (4th Edition).
John Wiley & Sons, 2002.
• Chandrupatla, T. R., y Belegundu, A. D.
Introduction to Finite Elements in Engineering (4th Edition).
Pearson, 2011.
• Segerlind, Larry J.
Applied Finite Element Analysis (2nd Edition).
John Wiley & Sons, 1984.
• "Finite Element Analysis of Solids and Structures" – Sudip S. Bhattacharjee
(2021).
• "Método de los Elementos Finitos" – Pilar Ariza Moreno y Andrés Sáez Pérez
(1999).
• "Cálculo de Estructuras por el Método de Elementos Finitos. Análisis Estático
Lineal. Volumen 1. Sólidos" – Eugenio Oñate (2016).
• "Método de los Elementos Finitos: Introducción a ANSYS" – Pilar Ariza
Moreno y Andrés Sáez Pérez (1999).
• "Introduction to Finite Element Methods" – Christian Clason (2017).
• M. Cerrolaza (2007): El Método de los Elementos Finitos en Ingeniería y
Ciencias Aplicadas: Teoría y Programas, CDCH-UCV, Caracas, 2007.
• J. J. Benito Muñoz, R. Álvarez Cabal, F. Ureña Prieto, E. Salete Casino, E.
Aranda Ortega (2014): Introducción al Método de los Elementos Finitos,
UNED, Madrid, 2014.
• O. C. Zienkiewicz, R. L. Taylor (1994): El Método de los Elementos Finitos,
volúmenes 1 y 2. CIMNE-Mc Graw Hill, 1994.
• Courant, R. (1943). Variational methods for the solution of problems of equilibrium
and vibrations. Bulletin of the American Mathematical Society, 49(1), 1–23.
• Turner, M. J., Clough, R. W., Martin, H. C., & Topp, L. J. (1956). Stiffness and
deflection analysis of complex structures. Journal of the Aeronautical Sciences,
23(9), 805–823.
• González Carbonell, R. (2015). Una introducción al análisis por elementos finitos:
Aplicaciones y ejemplos.
• Bhattacharjee, S. S. (2021). Finite Element Analysis of Solids and Structures.
• Beas Bernuy, H. (2016). Análisis no lineal de estructuras de concreto armado
mediante el método de los elementos finitos. Tesis de pregrado, UNI.
• Chinchay Poma, L. (2017). Modelación del concreto sometido a compresión
uniaxial y triaxial utilizando el método de elementos finitos. Tesis de pregrado,
PUCP.
• Li, Q., et al. (2022). Compressive Behavior of Concrete-Filled Steel Tubular
Columns Using Finite Element Analysis. Engineering Structures.
• Xiong, Y. (s.f.). Finite Element Modelling of Sandwich Structures under Axial
Loading. Abaqus Documentation.
