El método de elementos finitos (MEF) es una poderosa técnica desarrollada originalmente para la solución numérica de problemas complejos en mecánica estructural , y sigue siendo el método de elección para sistemas complejos. En el MEF, el sistema estructural se modela mediante un conjunto de elementos finitos apropiados interconectados en puntos discretos llamados nodos. Los elementos pueden tener propiedades físicas como espesor, coeficiente de expansión térmica , densidad , módulo de Young , módulo de corte y relación de Poisson .
El origen del método finito se remonta al análisis matricial de estructuras [1] [2] , donde se introdujo el concepto de enfoque matricial de desplazamiento o rigidez. Los conceptos de elementos finitos se desarrollaron basándose en métodos de ingeniería en la década de 1950. El método de los elementos finitos obtuvo su verdadero impulso en las décadas de 1960 y 1970 por parte de John Argyris y sus colaboradores; en la Universidad de Stuttgart , por Ray W. Clough ; en la Universidad de California, Berkeley , por Olgierd Zienkiewicz y sus compañeros de trabajo Ernest Hinton , Bruce Irons ; [3] en la Universidad de Swansea , por Philippe G. Ciarlet ; en la Universidad de París ; en la Universidad de Cornell , por Richard Gallagher y colaboradores. Los trabajos originales como los de Argyris [4] y Clough [5] se convirtieron en la base de los métodos actuales de análisis estructural de elementos finitos.
Elementos unidimensionales rectos o curvos con propiedades físicas como rigidez axial, de flexión y de torsión. Este tipo de elemento es adecuado para modelar cables, tirantes, cerchas, vigas, rigidizadores, rejillas y pórticos. Los elementos rectos suelen tener dos nodos, uno en cada extremo, mientras que los elementos curvos necesitarán al menos tres nodos, incluidos los nodos finales. Los elementos se ubican en el eje centroidal de los miembros reales.
Los elementos están interconectados sólo en los nodos exteriores y en conjunto deberían cubrir todo el dominio con la mayor precisión posible. Los nodos tendrán desplazamientos nodales (vectoriales) o grados de libertad que pueden incluir traslaciones, rotaciones y, para aplicaciones especiales, derivadas de desplazamientos de orden superior. Cuando los nodos se desplazan, arrastrarán los elementos de cierta manera dictada por la formulación del elemento. En otras palabras, los desplazamientos de cualquier punto del elemento se interpolarán a partir de los desplazamientos nodales, y esta es la razón principal de la naturaleza aproximada de la solución.
Desde el punto de vista de la aplicación, es importante modelar el sistema de manera que:
Los paquetes de software comerciales a gran escala a menudo proporcionan funciones para generar la malla y la visualización gráfica de la entrada y la salida, lo que facilita enormemente la verificación tanto de los datos de entrada como de la interpretación de los resultados.
Si bien la teoría de FEM se puede presentar en diferentes perspectivas o énfasis, su desarrollo para el análisis estructural sigue el enfoque más tradicional a través del principio de trabajo virtual o el principio de energía potencial total mínima . El enfoque del principio de trabajo virtual es más general ya que es aplicable a comportamientos de materiales tanto lineales como no lineales. El método del trabajo virtual es una expresión de conservación de energía : para sistemas conservadores, el trabajo agregado al sistema por un conjunto de fuerzas aplicadas es igual a la energía almacenada en el sistema en forma de energía de deformación de los componentes de la estructura.
El principio de desplazamientos virtuales para el sistema estructural expresa la identidad matemática del trabajo virtual externo e interno:
En otras palabras, la suma del trabajo realizado sobre el sistema por el conjunto de fuerzas externas es igual al trabajo almacenado como energía de deformación en los elementos que componen el sistema.
El trabajo interno virtual en el lado derecho de la ecuación anterior se puede encontrar sumando el trabajo virtual realizado en los elementos individuales. Esto último requiere que se utilicen funciones de fuerza-desplazamiento que describan la respuesta de cada elemento individual. Por tanto, el desplazamiento de la estructura se describe por la respuesta de elementos individuales (discretos) colectivamente. Las ecuaciones se escriben sólo para el pequeño dominio de elementos individuales de la estructura en lugar de una única ecuación que describa la respuesta del sistema como un todo (un continuo). Esto último daría lugar a un problema intratable, de ahí la utilidad del método de los elementos finitos. Como se muestra en las secciones siguientes, la ecuación ( 1 ) conduce a la siguiente ecuación de equilibrio gobernante para el sistema:
dónde
Una vez que se tienen en cuenta las restricciones de los apoyos, los desplazamientos nodales se encuentran resolviendo el sistema de ecuaciones lineales ( 2 ), simbólicamente:
Posteriormente, las deformaciones y tensiones en elementos individuales se pueden encontrar de la siguiente manera:
dónde
Aplicando la ecuación de trabajo virtual ( 1 ) al sistema, podemos establecer las matrices de elementos , así como la técnica de ensamblaje de las matrices del sistema y . Otras matrices como , y son valores conocidos y se pueden configurar directamente a partir de la entrada de datos .
Sea el vector de desplazamientos nodales de un elemento típico. Los desplazamientos en cualquier otro punto del elemento se pueden encontrar mediante el uso de funciones de interpolación como, simbólicamente:
dónde
La ecuación ( 6 ) da lugar a otras cantidades de gran interés:
donde = matriz de operadores diferenciales que convierten desplazamientos en deformaciones utilizando la teoría de la elasticidad lineal . La ecuación ( 7 ) muestra que la matriz B en ( 4 ) es
Para un elemento típico de volumen , el trabajo virtual interno debido a los desplazamientos virtuales se obtiene sustituyendo ( 5 ) y ( 9 ) en ( 1 ):
Principalmente para facilitar la referencia, ahora se pueden definir las siguientes matrices pertenecientes a elementos típicos:
Estas matrices generalmente se evalúan numéricamente utilizando la cuadratura gaussiana para la integración numérica . Su uso se simplifica ( 10 ) a lo siguiente:
Dado que el vector de desplazamiento nodal q es un subconjunto de los desplazamientos nodales del sistema r (para compatibilidad con elementos adyacentes), podemos reemplazar q con r expandiendo el tamaño de las matrices de elementos con nuevas columnas y filas de ceros:
donde, por simplicidad, utilizamos los mismos símbolos para las matrices de elementos, que ahora tienen un tamaño ampliado y filas y columnas adecuadamente reorganizadas.
Sumar el trabajo virtual interno ( 14 ) para todos los elementos da el lado derecho de ( 1 ):
Considerando ahora el lado izquierdo de ( 1 ), el trabajo virtual externo del sistema consta de:
donde hemos introducido matrices de elementos adicionales definidas a continuación:
Nuevamente, la integración numérica es conveniente para su evaluación. Una sustitución similar de q en ( 17a ) por r da, después de reorganizar y expandir los vectores :
Sumando ( 16 ), ( 17b ) e igualando la suma a ( 15 ) se obtiene:
Dado que los desplazamientos virtuales son arbitrarios, la igualdad anterior se reduce a:
La comparación con ( 2 ) muestra que:
En la práctica, las matrices de elementos no se expanden ni se reordenan. En cambio, la matriz de rigidez del sistema se ensambla agregando coeficientes individuales donde los subíndices ij, kl significan que los desplazamientos nodales del elemento coinciden respectivamente con los desplazamientos nodales del sistema . De manera similar, se ensambla agregando coeficientes individuales a donde coincide . Esta adición directa de into le da al procedimiento el nombre de Método de rigidez directa .