Método de los elementos finitos en la mecánica estructural

NOTA: Este artículo está siendo traducido del artículo en inglés, debe eliminarse la página que redirige a este artículo por tener un título repetido

Método de los Elementos Finitos (FEM) es una poderosa técnica originalmente desarrollada para soluciones numéricas de problemas complejos en Mecánica estructural, y sigue siendo el método de elección para sistemas complejos. En el FEM, el sistema estructural es modelado por un conjunto de elementos finitos apropiadamente interconectados en puntos llamados nodos. Los elementos deben de tener propiedades físicas tales como espesor, Coeficiente de dilatación, Densidad, Módulo de elasticidad, Modulo de cortante y Coeficiente de Poisson.

Contenido

Propiedades de los Elementos

  • Elementos rectos o curvados unidimencionales con propiedades físicas tales como rigidez axial, doblado, y torcional. Este tipo de elemento es adaptable para modelar cables, tirantes, trusses, vigas, refuerzos, redes y marcos. Los elementos rectos usualmente tienen dos nodos, uno en cada extremo, mientras los elementos curvados necesitaran al menos tres nodos incluyendo los nodos extremos. Los elementos están posicionados en el eje Centroidal de los miembros actuales.
  • Elementos bidimensionales para acción de Membrana ( esfuerzo en el plano, Deformación) en elplano y/o acción de doblado (placas y cascarones). ellos deben tener una variedad de formas tales como planas o curvadas Triángulos y Cuadriláteros. Los nodos son usualmente colocados en las esquinas de los elementos y, si es necesario para mayor precision, nodos adicionales pueden ser colocados a lo largo de los filos de los elementos o incluso dentro de los elementos. Los elementos son posicionados a la mitad de la superficie del espesor actual de la capa.
  • Elementos en forma de Toro para problemas ejesimetricos tales como delgados, placas gruesas, cascarones, y solidos. La sección transversal de los elementos es similar a los tipos previamente descritos: unidimencional para placas delgadas y cascarones, y bidimensional para solidos, y placas gruesas y cascarones.
  • Elementos tridimensionales para modelas solidos en 3D tales como componentes de Máquina , Presas, Terraplén o masas solidas. Formas de elementos comunes incluyen Tetraedros y hexahedros. Los nodos son colocados en los vertices y posiblemente en las caras de elementos o dentro del elemento.

Interconexion de elementos y desplazamientos

Los elementos están interconectados solo en los nodos exteriores, y ellos completamente deberían cubrir el dominio entero tan preciso como sea posible. Los nodos tendrán (vector) desplazamiento nodal o grados de libertad los cuales deben incluir traslaciones, rotaciones, y para aplicaciones especiales, Derivadas de alto orden de desplazamientos. Cuando los nodos se desplazan, ellos arrastraran los elementos a lo largo en una cierta manera dictada por la formulación del elemento. En otras palabras, desplazamientos de algún punto en el elemento serán interpolados desde los desplazamientos nodales, y esta es la principal razón para la naturaleza aproximada de la solución.

Consideraciones practicas

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

  • Las condiciones de simetría o asimetría sean explotadas en disposición de reducir el tamaño del dominio.
  • La compatibilidad de desplazamientos, incluyendo alguna discontinuidad requerida, es asegurada en los nodos, y preferiblemente, a lo largo del filo de los elementos también, particularmente cuando elementos adyacentes son de diferentes tipos, material o densidad. La compatibilidad de desplazamiento de muchos nodos puede ser usualmente impuesta mediante relaciones de restricción--Cuando tal característica no está disponible en el paquete de software, un modelo físico que impone las restricciones debe ser usado en su lugar.
  • El comportamiento de los elementos captura las acciones dominantes del actual sistema, ambos local y globalmente.
  • La malla de elementos es suficientemente fina en disposición de tener una precisión aceptable. Para evaluar la precisión, la malla es refinada hasta que los resultados importantes muestran pequeños cambios. Para mayor precisión, la Relación de aspecto de los elementos debería ser tan cercana a la unidad como sea posible, y pequeños elementos son usados sobre las partes de grandes Gradientes de esfuerzo.
  • Restricciones adecuadas de soporte son impuestas, con especial atención, a nodos en ejes de simetría.

Paquetes de software comercial de gran escala usualmente proveen facilidades para generar la malla, salidas gráficas de entradas y salidas, los cuales grandemente facilitan la verificación de ambos datos de entrada e interpretación de resultados.

Vision general teórica de Formulación de Desplazamientos-FEM: Desde elementos a sistema a solución

Mientras la teoría de FEM puede ser representada en diferentes perspectivas o énfasis, su desarrollo para Análisis estructural siguiendo los más tradicionales enfoques mediante el 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 como es aplicable a ambos comportamientos de material lineal y no lineal.

El principio de desplazamientos virtuales para el sistema estructural expresando la identidad matemática de trabajo virtual interno y externo:

\mbox{Trabajo virtual externo} = \int_{V}\delta\boldsymbol{\epsilon}^T \boldsymbol{\sigma} \, dV \qquad \mathrm{(1)}

El trabajo virtual interno en el lado derecho de la ecuación superior debe ser encontrado sumando el trabajo virtual en los elementos individuales--Este es el paso crucial donde necesitaremos funciones de desplazamiento escrita solo para los bastantes dominios pequeños que sobre el sistema completo. como se muestra en las subsecuentes secciones, Eq.(1) lleva a las siguientes ecuaciones de equilibrio governando para el sistema:

\mathbf{R} = \mathbf{Kr} + \mathbf{R}^o \qquad \qquad \qquad \mathrm{(2)}

donde

\mathbf{R} = vector de fuerzas nodales, representando fuerzas externas aplicadas a los nodos del sistema.
\mathbf{r} = vector de desplazamientos nodales del sistema, los cuales serán, por interpolación, desplazamientos de cedencia en algún punto de la malla de elementos finitos.
\mathbf{R}^o = vector de fuerzas nodales equivalentes, representando todos los efectos externos difernetes de las fuerzas nodales las cuales son todavía incluidas en el vector de fuerza nodal precedente R. Estos efectos externos deben incluir fuerzas superficiales distribuidas o concentradas, fuerzas de cuerpo, efectos térmicos, esfuerzos iniciales y restricciones.
\mathbf{K} = matriz de rigides del sistema, la cual debe ser establecida por ensanblando las matrices de riguides de los elementos :\mathbf{k}^e .

Una vez las restricciones en los soportes son representadas, los desplazamientos nodales son encontrados resolviendo el Sistema de ecuaciones lineales (2), simbólicamente:

\mathbf{r} = \mathbf{K}^{-1} (\mathbf{R}-\mathbf{R}^o ) \qquad \qquad \qquad \mathrm{(3)}

Subsecuentemente, las restricciones y esfuerzos en elementos individuales debe ser encontrada como sigue:

\mathbf{\epsilon} = \mathbf{Bq} \qquad \qquad \qquad \qquad \mathrm{(4)}
\mathbf{\sigma} = \mathbf{E}(\mathbf{\epsilon} - \mathbf{\epsilon}^o)+\mathbf{\sigma}^o = \mathbf{E}(\mathbf{Bq} - \mathbf{\epsilon}^o)+\mathbf{\sigma}^o\qquad \qquad \qquad \mathrm{(5)}

donde

\mathbf{q} = vector de desplazamientos nodales de elementos--un subconjunto de los vectores de desplazamiento del sistema r que pertenecen a los elementos bajo consideración.
\mathbf{B} = matriz de restricciones-desplazamientos que transforman los desplazamientos nodales q a restricciones en algún punto del elemento.
\mathbf{E} = matriz de elasticidad que trasnforma restricciones efectivas a esfuerzos a algún punto en el elemento.
\mathbf{\epsilon}^o = vector de restricciones iniciales en el elemento.
\mathbf{\sigma}^o = vector esfuerzos iniciales en el elemento.

Aplicando el Trabajo virtual ecuación (1) al sistema, podemos establecer las matrices de elementos \mathbf{B}, \mathbf{k}^e así como la técnica de ensamblar las matrices del sistema \mathbf{R}^o y \mathbf{K}. Otras matrices tales como \mathbf{\epsilon}^o , \mathbf{\sigma}^o , \mathbf{R} and \mathbf{E} pueden ser directamente creadas de los datos de entrada.

Interpolación o funciones de forma

Permite a \mathbf{q} ser el vector de desplazamientos nodales de un elemento típico. Los desplazamientos en algún punto del elemento deben ser encontrados por funciones de Interpolación , simbólicamente:

\mathbf{u} = \mathbf{N} \mathbf{q} \qquad \qquad \qquad \mathrm{(6)}

donde

\mathbf{u} = vector de desplazamientos en algún punto {x, y, z} del elemento.
\mathbf{N} = matrix de funciones forma actuando como funciones de Interpolación .

La ecuación (6) da subida a otras cantidades de grandes intereses:

  • Desplazamientos virtuales consistentes con desplazamientos virtuales nodales:  \delta \mathbf{u} = \mathbf{N} \delta \mathbf{q} \qquad \qquad \qquad \mathrm{(6b)}
  • restricciones en los elementos:\mathbf{\epsilon} = \mathbf{Du} = \mathbf{DNq} \qquad \qquad \qquad \qquad \mathrm{(7)}
donde \mathbf{D} = matriz de operadores diferenciales que convierten desplazamientos a restricciones usando teoría de elasticidad lineal . Eq.(7) muestra qeu la matriz B en (4) es
\mathbf{B} = \mathbf{DN} \qquad \qquad \qquad \qquad \mathrm{(8)}
  • Restricciones virtuales consistentes con el desplazamiento nodal virtual de los elementos:  \delta \boldsymbol{\epsilon} = \mathbf{B} \delta \mathbf{q} \qquad \qquad \qquad \qquad \mathrm{(9)}

Trabajo virtual interno en un elemento típico

Para un elemento típico de volumen Ve, el trabajo virtual interno debido a los desplazamientos virtuales es obtenido por sustitución de (5) y (9) dentro de (1):

\mbox{Trabajo virtual interno} = \int_{V^e}\delta\boldsymbol{\epsilon}^T \boldsymbol{\sigma} \, dV^e = \delta\ \mathbf{q}^T \int_{V^e} \mathbf{B}^T \big\{\mathbf{E}(\mathbf{Bq} - \mathbf{\epsilon}^o)+\mathbf{\sigma}^o\big\} \, dV^e \qquad \mathrm{(10)}

Matrices de Elementos

Primeeramente para la conveniencia de referencias, las siguientes matrices pertenecen a elementos típicos que deben ser definidos:

Matriz de rigidez de elementos  \mathbf{k}^e = \int_{V^e} \mathbf{B}^T \mathbf{E} \mathbf{B} \, dV^e \qquad \mathrm{(11)}
Vector de carga de elemento equivalente  \mathbf{Q}^{oe} = \int_{V^e} - \mathbf{B}^T \big( \mathbf{E}\mathbf{\epsilon}^o - \mathbf{\sigma}^o\big ) \, dV^e \qquad \mathrm{(12)}

Estas matrices son usualmente evaluadas numéricamente usando la Cuadratura gaussiana para Integración numérica. Su uso simplifica (10) a lo siguiente:

\mbox{Internal virtual work} = \delta\ \mathbf{q}^T \big( \mathbf{k}^e \mathbf{q} + \mathbf{Q}^{oe} \big)  \qquad \mathrm{(13)}

Trabajo virtual de elementos en términos de los desplazamientos nodales del sistema

desde que el vector de desplazamientos nodales q es un sujeto del sistema de desplazamientos nodales r (para la compatibilidad con los elementos adyacentes), podemos reemplazar q con r expandiendo el tamaño de las matrices de elementos con nuevas columnas y filas de ceros:

\mbox{Trabajo virtual interno} = \delta\ \mathbf{r}^T \big( \mathbf{k}^e \mathbf{r} + \mathbf{Q}^{oe} \big)  \qquad \mathrm{(14)}

Donde, por sinplicidad, usamos los mismos sinbolos para las matrices de elementos, las cuales ahora an expandido su tamaño así como adecuadamente reordenado sus filas y columnas.

Trabajo virtual del sistema

Sumando los trabajos virtuales internos (14) para todos los elementos dados en el lado derecho de (1):

\mbox{Trabajo virtual inerno del sistema} = \sum_{e} \delta\ \mathbf{r}^T \big( \mathbf{k}^e \mathbf{r} + \mathbf{Q}^{oe} \big)  = \delta\ \mathbf{r}^T \big( \sum_{e} \mathbf{k}^e \big)\mathbf{r} + \delta\ \mathbf{r}^T \sum_{e} \mathbf{Q}^{oe}    \qquad \mathrm{(15)}

Considerando ahora el lado izquierdo de (1), el trabajo virtual externo del sistema consiste de:

  • El trabajo echo por las fuerzas nodales R:  \delta\ \mathbf{r}^T \mathbf{R} \qquad \mathrm{(16)}
  • El trabajo echo por las fuerzas externas  \mathbf{T}^e en la parte  \mathbf{S}^e de los filos de los elementos o superficies, y po rlas fuerzas del cuerpo  \mathbf{f}^e
 \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:
 \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  -\delta\ \mathbf{q}^T \sum_{e} (\mathbf{Q}^{te} +  \mathbf{Q}^{fe}) \qquad \mathrm{(17a)}
Donde emos introducido matrices de elementos adicionales definidas abajo:
 \mathbf{Q}^{te} =  -\int_{S^e} \mathbf{N}^T \mathbf{T}^e \, dS^e  \qquad \mathrm{(18a)}
 \mathbf{Q}^{fe} =  -\int_{V^e} \mathbf{N}^T \mathbf{f}^e \, dV^e  \qquad \mathrm{(18b)}
Otra vez, la Integración numérica es conveniente para su evaluación. Un reemplazo similar de q en (17a) con r dados, reordenando y expandiendo después los vectores  \mathbf{Q}^{te}, \mathbf{Q}^{fe} :
 -\delta\ \mathbf{r}^T \sum_{e} (\mathbf{Q}^{te} +  \mathbf{Q}^{fe}) \qquad \mathrm{(17b)}

Ensamblaje de matrices de sistemas

Añadiendo (16), (17b) e igualando la suma a (15) da:  \delta\ \mathbf{r}^T \mathbf{R} -\delta\ \mathbf{r}^T \sum_{e} (\mathbf{Q}^{te} +  \mathbf{Q}^{fe}) =  \delta\ \mathbf{r}^T \big( \sum_{e} \mathbf{k}^e \big)\mathbf{r} + \delta\ \mathbf{r}^T \sum_{e} \mathbf{Q}^{oe}

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

 \mathbf{R} = \big( \sum_{e} \mathbf{k}^e \big)\mathbf{r} + \sum_{e} \big( \mathbf{Q}^{oe} + \mathbf{Q}^{te} +  \mathbf{Q}^{fe} \big)

Comparando con (2) muestra que:

  • La matriz de rigides del sistema es obtenida sumando las matrices de rigidez de los elementos:
 \mathbf{K} = \sum_{e} \mathbf{k}^e
  • El vector de fuerzas nodales equivalentes es obtenido sumando los vectores de carga de los elementos:
 \mathbf{R}^o = \sum_{e} \big( \mathbf{Q}^{oe} + \mathbf{Q}^{te} +  \mathbf{Q}^{fe} \big)

En la práctica, las matrices de elementos ninguna es expandida o reordenada. En lugar, la matriz de rigides del sistema  \mathbf{K} es ensanblada añadiendo coeficientes individuales  {k}_{ij}^e a Kkl donde los subindices ij, kl significan que los desplazamientos nodales de los elementos  {q}_{i}^e, {q}_{j}^e coinciden respectivamente con los desplazamientos nodales del sistema rk,rl. Similarmente,  \mathbf{R}^o es ensamblado añadiendo coeficientes individuales  {Q}_{i}^e to  {R}^o_{k} donde  {q}_{i}^e matches rk. Esta adición directa de  {k}_{ij}^e dentro de Kkl da al procedimiento el nombre de Método de rigidez directa (Método matricial de la rigidez).

Véase también


Wikimedia foundation. 2010.

Mira otros diccionarios:

  • Método de los elementos finitos — Solución de MEF en 2D para una configuración de un magnetostato, (las líneas muestran la dirección de la densidad de flujo calculada, y el color, su magnitud) …   Wikipedia Español

  • Método de elementos de frontera — El método de elementos de frontera, es también conocido como método de elementos de contorno o método de elementos de borde (es mejor conocido como BEM, del inglés boundary element method) es un método numérico para resolver ecuaciones en… …   Wikipedia Español

  • Análisis de elementos finitos — Saltar a navegación, búsqueda Visualización de como un vehículo se deforma en un impacto asimétrico, obtenido usando elementos finitos. El análisis por elementos finitos (FEA por sus siglas en inglés para: Finite Element Analysis) es una técnica… …   Wikipedia Español

  • Ingeniería mecánica — Este artículo o sección necesita referencias que aparezcan en una publicación acreditada, como revistas especializadas, monografías, prensa diaria o páginas de Internet fidedignas. Puedes añadirlas así o avisar …   Wikipedia Español

  • Elasticidad (mecánica de sólidos) — Una varilla elástica vibrando, es un ejemplo de sistema donde la energía potencial elástica se transforma en energía cinética y viceversa. En física e ingeniería, el término elasticidad designa la propiedad mecánica de ciertos materiales de… …   Wikipedia Español

  • Escuela Técnica Superior de Ingenieros Aeronáuticos — Saltar a navegación, búsqueda Escuela Técnica Superior de Ingenieros Aeronáuticos Acrónimo E.T.S.I.A. Fundac …   Wikipedia Español

  • Reacción mecánica — Saltar a navegación, búsqueda En ingeniería estructural e ingeniería mecánica, una reacción es una fuerza de sujeción de un elemento resistente al suelo u otro elemento de grandes dimensiones que sirve de soporte al elemento resistente. Métodos… …   Wikipedia Español

  • Resistencia de materiales — Saltar a navegación, búsqueda Para otros usos de este término, véase Resistencia. La resistencia de materiales clásica es una disciplina de la ingeniería mecánica y la ingeniería estructural que estudia los sólidos deformables mediante modelos… …   Wikipedia Español

  • Analisis modal utilizando FEM — Saltar a navegación, búsqueda El objetivo del Analisis modal en la mecánica estructural es determinar las formas modales naturales y frecuencias de un objeto o estructura durante vibración libre. Es común usar el Método de los elementos finitos… …   Wikipedia Español

  • ANSYS — ANSYS, Inc. Fue fundada en 1970 (Swanson Analysis Systems, Inc.) utilizando aproximadamente 1700 empleados. La mayoría con experiencia en elemento finito y dinámica de fluido computacional. ANSYS desarrolla, comercializa y presta soporte a la… …   Wikipedia Español

Compartir el artículo y extractos

Link directo
Do a right-click on the link above
and select “Copy Link”