Método de elementos finitos en mecánica estructural.

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 .

Historia

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.

• Elementos bidimensionales que resisten sólo fuerzas en el plano por acción de membrana ( tensión plana , deformación plana ), y placas que resisten cargas transversales por acción de corte transversal y flexión (placas y carcasas ). Pueden tener una variedad de formas, como triángulos y cuadriláteros planos o curvos . Los nodos generalmente se colocan en las esquinas del elemento y, si es necesario para una mayor precisión, se pueden colocar nodos adicionales a lo largo de los bordes del elemento o incluso dentro del elemento. Los elementos se sitúan en la mitad de la superficie del espesor de capa real.
• Elementos en forma de toro para problemas de simetría axial como membranas, placas gruesas, cáscaras y sólidos. La sección transversal de los elementos es similar a los tipos descritos anteriormente: unidimensional para placas delgadas y conchas, y bidimensional para sólidos, placas gruesas y conchas.
• Elementos tridimensionales para modelar sólidos tridimensionales como componentes de máquinas , presas , terraplenes o masas de suelo. Las formas de elementos comunes incluyen tetraédricos y hexaédricos . Los nodos se colocan en los vértices y posiblemente en las caras del elemento o dentro del elemento.

Interconexión y desplazamiento de elementos.

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 una 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.

Consideraciones prácticas

Desde el punto de vista de la aplicación, es importante modelar el sistema de manera que:

• Se aprovechan las condiciones de simetría o antisimetría para reducir el tamaño del modelo.
• La compatibilidad de desplazamiento, incluida cualquier discontinuidad requerida, se garantiza en los nodos, y preferiblemente también a lo largo de los bordes del elemento, particularmente cuando los elementos adyacentes son de diferentes tipos, materiales o espesores. La compatibilidad de los desplazamientos de muchos nodos generalmente se puede imponer mediante relaciones de restricción.
• Los comportamientos de los elementos deben capturar las acciones dominantes del sistema real, tanto a nivel local como global.
• La malla del elemento debe ser lo suficientemente fina para producir una precisión aceptable. Para evaluar la precisión, la malla se refina hasta que los resultados importantes muestren pocos cambios. Para una mayor precisión, la relación de aspecto de los elementos debe ser lo más cercana posible a la unidad y se utilizan elementos más pequeños sobre las partes con mayor gradiente de tensión .
• Se imponen restricciones de soporte adecuadas prestando especial atención a los nodos en los ejes de simetría.

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.

Descripción teórica de la formulación por desplazamiento FEM: de los elementos al sistema y a la solución

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 difícil de resolver, 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

${\displaystyle \mathbf {R} }$= vector de fuerzas nodales, que representa fuerzas externas aplicadas a los nodos del sistema.

${\displaystyle \mathbf {K} }$= matriz de rigidez del sistema, que es el efecto colectivo de las matrices de rigidez de los elementos individuales  : . ${\displaystyle \mathbf {k} ^{e}}$

${\displaystyle \mathbf {r} }$= vector de los desplazamientos nodales del sistema.

${\displaystyle \mathbf {R} ^{o}}$= vector de fuerzas nodales equivalentes, que representa todos los efectos externos distintos de las fuerzas nodales que ya están incluidos en el vector de fuerzas nodales anterior R. Estos efectos externos pueden incluir fuerzas superficiales distribuidas o concentradas, fuerzas corporales, efectos térmicos, tensiones y deformaciones iniciales.

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

${\displaystyle \mathbf {q} }$= vector de desplazamientos nodales: un subconjunto del vector de desplazamiento del sistema r que pertenece a los elementos considerados.

${\displaystyle \mathbf {B} }$= matriz deformación-desplazamiento que transforma los desplazamientos nodales q en deformaciones en cualquier punto del elemento.

${\displaystyle \mathbf {E} }$= matriz de elasticidad que transforma deformaciones efectivas en tensiones en cualquier punto del elemento.

${\displaystyle \mathbf {\epsilon } ^{o}}$= vector de deformaciones iniciales en los elementos.

${\displaystyle \mathbf {\sigma } ^{o}}$= vector de tensiones iniciales en los elementos.

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 . ${\displaystyle \mathbf {B} }$ ${\displaystyle \mathbf {k} ^{e}}$ ${\displaystyle \mathbf {R} ^{o}}$ ${\displaystyle \mathbf {K} }$ ${\displaystyle \mathbf {\epsilon } ^{o}}$ ${\displaystyle \mathbf {\sigma } ^{o}}$ ${\displaystyle \mathbf {R} }$ ${\displaystyle \mathbf {E} }$

Funciones de interpolación o forma.

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: ${\displaystyle \mathbf {q} }$

dónde

${\displaystyle \mathbf {u} }$= vector de desplazamientos en cualquier punto {x,y,z} del elemento.

${\displaystyle \mathbf {N} }$= matriz de funciones de forma que sirven como funciones de interpolación .

La ecuación ( 6 ) da lugar a otras cantidades de gran interés:

• Desplazamientos virtuales que son función de desplazamientos nodales virtuales:

• Deformaciones en los elementos que resultan de desplazamientos de los nodos del elemento:

  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 ${\displaystyle \mathbf {D} }$

• Deformaciones virtuales consistentes con los desplazamientos nodales virtuales del elemento:

Trabajo virtual interno en un elemento típico.

Para un elemento típico de volumen , el trabajo virtual interno debido a los desplazamientos virtuales se obtiene sustituyendo ( 5 ) y ( 9 ) en ( 1 ): ${\displaystyle V^{e}}$

Matrices de elementos

Principalmente para facilitar la referencia, ahora se pueden definir las siguientes matrices pertenecientes a elementos típicos:

Matriz de rigidez del elemento

Vector de carga de elemento equivalente

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:

Elemento de trabajo virtual en términos de desplazamientos nodales del sistema.

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.

trabajo virtual del sistema

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:

• El trabajo realizado por las fuerzas nodales R :

• El trabajo realizado por fuerzas externas por parte de los bordes o superficies de los elementos, y por las fuerzas del cuerpo. ${\displaystyle \mathbf {T} ^{e}}$ ${\displaystyle \mathbf {S} ^{e}}$ ${\displaystyle \mathbf {f} ^{e}}$

  ${\displaystyle \sum _{e}\int _{S^{e}}\delta \ \mathbf {u} ^{T}\mathbf {T} ^{e}\,dS^{e}+\sum _{e}\int _{V^{e}}\delta \ \mathbf {u} ^{T}\mathbf {f} ^{e}\,dV^{e}}$

  La sustitución de ( 6b ) da:

  ${\displaystyle \delta \ \mathbf {q} ^{T}\sum _{e}\int _{S^{e}}\mathbf {N} ^{T}\mathbf {T} ^{e}\,dS^{e}+\delta \ \mathbf {q} ^{T}\sum _{e}\int _{V^{e}}\mathbf {N} ^{T}\mathbf {f} ^{e}\,dV^{e}}$

  o

  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 : ${\displaystyle \mathbf {Q} ^{te},\mathbf {Q} ^{fe}}$

Montaje de matrices del sistema.

Sumando ( 16 ), ( 17b ) e igualando la suma a ( 15 ) se obtiene: ${\displaystyle \delta \ \mathbf {r} ^{T}\mathbf {R} -\delta \ \mathbf {r} ^{T}\sum _{e}\left(\mathbf {Q} ^{te}+\mathbf {Q} ^{fe}\right)=\delta \ \mathbf {r} ^{T}\left(\sum _{e}\mathbf {k} ^{e}\right)\mathbf {r} +\delta \ \mathbf {r} ^{T}\sum _{e}\mathbf {Q} ^{oe}}$

Dado que los desplazamientos virtuales son arbitrarios, la igualdad anterior se reduce a: ${\displaystyle \delta \ \mathbf {r} }$

${\displaystyle \mathbf {R} =\left(\sum _{e}\mathbf {k} ^{e}\right)\mathbf {r} +\sum _{e}\left(\mathbf {Q} ^{oe}+\mathbf {Q} ^{te}+\mathbf {Q} ^{fe}\right)}$

La comparación con ( 2 ) muestra que:

• La matriz de rigidez del sistema se obtiene sumando las matrices de rigidez de los elementos:

  ${\displaystyle \mathbf {K} =\sum _{e}\mathbf {k} ^{e}}$

• El vector de fuerzas nodales equivalentes se obtiene sumando los vectores de carga de los elementos:

  ${\displaystyle \mathbf {R} ^{o}=\sum _{e}\left(\mathbf {Q} ^{oe}+\mathbf {Q} ^{te}+\mathbf {Q} ^{fe}\right)}$

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 . ${\displaystyle \mathbf {K} }$ ${\displaystyle {k}_{ij}^{e}}$ ${\displaystyle {K}_{kl}}$ ${\displaystyle {q}_{i}^{e},{q}_{j}^{e}}$ ${\displaystyle {r}_{k},{r}_{l}}$ ${\displaystyle \mathbf {R} ^{o}}$ ${\displaystyle {Q}_{i}^{e}}$ ${\displaystyle {R}_{k}^{o}}$ ${\displaystyle {q}_{i}^{e}}$ ${\displaystyle {r}_{k}}$ ${\displaystyle {k}_{ij}^{e}}$ ${\displaystyle {K}_{kl}}$

Ver también

• Método de elementos finitos
• Método de flexibilidad
• Método de rigidez de la matriz
• Análisis modal usando FEM
• Lista de paquetes de software de elementos finitos
• Análisis estructural
• Trabajo virtual
• elemento finito de intervalo

Referencias

1. Análisis matricial de estructuras enmarcadas, tercera edición de Jr. William Weaver, James M. Gere, Springer-Verlag New York, LLC, ISBN  978-0-412-07861-3 , 1966
2. Teoría del análisis estructural matricial, JS Przemieniecki, McGraw-Hill Book Company, Nueva York, 1968
3. Hinton, Ernesto; Irons, Bruce (julio de 1968). "Suavizado de mínimos cuadrados de datos experimentales utilizando elementos finitos". Cepa . 4 (3): 24-27. doi :10.1111/j.1475-1305.1968.tb01368.x.
4. Argyris, JH y Kelsey, S. Teoremas de energía y análisis estructural Publicaciones científicas de Butterworth, Londres, 1954
5. Clough, RW, "El elemento finito en el análisis de tensiones planas". Actas, Segunda Conferencia de la ASCE sobre Computación Electrónica, Pittsburgh, septiembre de 1960