Metodo de Los Elementos Finitos

2019
METODO DE LOS ELEMENTOS FINITOS
El Método de los Elementos Finitos es un método numérico de resolución de problemas de

Mecánica de Sólidos que resulta de gran importancia por su utilidad práctica. Es una
herramienta de cálculo muy potente que permite al ingeniero estructuralista resolver infinidad
de problemas. Sin embargo, es un método que no proporciona la solución “exacta” a un
problema dado, sino que, en realidad, posibilita obtener una solución aproximada que, con el
juicio ingenieril que se le supone al calculista, puede ser más que suficiente para la resolución
de un problema práctico.
Su idea básica no puede ser más sencilla: dado un sólido, sometido a un sistema de cargas y
coaccionado por unas ligaduras, el método consiste en subdividir el sólido en pequeñas partes
(elementos) interconectadas entre sí a través de los nudos de los elementos, de manera que
suponemos que, el campo de desplazamientos en el interior de cada elemento, puede
expresarse en función de los desplazamientos que sufren los nudos del elemento
(desplazamientos nodales); posteriormente, se podrá determinar la matriz de rigidez de cada
elemento, las cuales una vez ensambladas (siguiendo los pasos del análisis matricial de
estructuras), permitirán la obtención de los desplazamientos en los nudos de cada elemento.
De esa manera, una vez conocidos dichos desplazamientos, podríamos determinar, de una
forma aproximada como ya se dijo antes, las tensiones y las deformaciones en el interior del
elemento.
El elemento finito lineal (barra biarticulada): planteamiento en coordenadas cartesianas y

obtención de su matriz de rigidez
En este apartado vamos a obtener la matriz de rigidez, en ejes locales, de una barra
biarticulada pero, ahora, utilizando las ideas básicas del Método de los Elementos Finitos. Se
trata, más bien, de un ejemplo de cómo se procede en este método porque, en realidad, la
matriz de rigidez que vamos buscando ya la hemos obtenido utilizando conceptos de
Resistencia de Materiales.
Denominemos u1 y u2 a los desplazamientos, según el eje x, que experimentan sus nudos y

admitamos que el campo de desplazamientos, en el interior del elemento, lo aproximamos
como:
u(x) = α0 + α1
donde las dos constantes α0 y α1 son desconocidas aunque las podríamos obtener imponiendo
que cuando x=x1, u(x1)=u1 y que cuando x=x2, u(x2)=u2. Es decir:
La ecuación anterior es un sistema de dos ecuaciones con dos incógnitas, en este caso. Su
solución proporcionaría:
Conviene hacer una precisión importante para, luego, poder comprobar la importancia del uso
de los elementos isoparamétricos que, posteriormente, veremos: para determinar el campo de
desplazamientos dentro del elemento (para lo que es necesario determinar los coeficientes α0
y α1) hemos tenido que resolver un sistema de ecuaciones. Esto puede no parecer una cosa
importante, sobre todo cuando se está analizando este elemento, pero cuando se desarrollen
elementos finitos con muchos más grados de libertad, y teniendo en cuenta que el proceso de
resolución para la obtención de los coeficientes de las funciones de forma habría que hacerlo
para todos y cada uno de los elementos de la malla, se puede convertir en un esfuerzo de
cálculo importante.
El campo de deformaciones dentro del elemento lo podremos obtener a partir del campo de
desplazamientos como:
que, como se observa, es constante a lo largo de todo el elemento. La tensión dentro del
mismo se podría obtener a partir de la deformación utilizando la ley de Hooke:
Elemento finito triangular de deformación constante: Matriz de rigidez
Si fuéramos capaces de dividir un sólido, considerado éste como un medio continuo, en

elementos de manera que lográsemos, para cada uno de ellos, podríamos pensar en seguir las
mismas pautas que hemos seguido en el problema precedente para resolver problemas de
Elasticidad mucho más complicados. Para ello, supongamos el elemento triangular de la Figura
11, que supondremos que trabaja en tensión plana, y que puede ayudarnos para resolver
problemas de sólidos bidimensionales como el que hemos considerado en las Figuras (1) y (2).
Vamos a emplear, como sistema de referencia global del elemento el formado por los ejes x-y
de la figura.
Planteamiento y formulación genéricos del Método de los elementos finitos
plantearemos las ecuaciones del Método de los Elemento Finitos de una forma general,
particularizando para el elemento triangular que hemos estudiado en el apartado anterior. Así,
consideremos el elemento triangular de la Figura 16, en el que sus nudos vienen definidos por
los números i,j,k.
SISTEMAS DISCRETOS Y SISTEMAS CONTINUOS
Al efectuar una clasificación de las estructuras, suelen dividirse en discretas o reticulares y

continuas. Las primeras son aquéllas que están formadas por un ensamblaje de elementos
claramente diferenciados unos de otros y unidos en una serie de puntos concretos, de tal
manera que el sistema total tiene forma de malla o retícula. La característica fundamental de
las estructuras discretas es que su deformación puede definirse de manera exacta mediante
un número finito de parámetros, como por ejemplo las deformaciones de los puntos de unión
de unos elementos y otros. De esta manera el equilibrio de toda la estructura puede
representarse mediante las ecuaciones de equilibrio en las direcciones de dichas
deformaciones.
Como contrapartida, en los sistemas continuos no es posible separar, a priori, el sistema en un
número finito de elementos discretos. Si se toma una parte cualquiera del sistema, el número
de puntos de unión entre dicha parte y el resto de la estructura es infinito, y es por lo tanto
imposible utilizar el mismo método que en los sistemas discretos, pues los puntos de unión
entre los distintos elementos, que allí aparecían de manera natural, no existen aquí.
Las estructuras continuas son muy frecuentes en ingeniería, como por ejemplo: bastidores de
máquinas, carrocerías de vehículos, losas de cimentación de edificios, vasijas de reactores,
elementos de máquinas (bielas, poleas, carcasas...), y para su análisis es necesario disponer de
un método que tenga en cuenta su naturaleza continua. Hasta la llegada del Método de los
Elementos Finitos (MEF), los sistemas continuos se abordaban analíticamente, pero por esa vía
sólo es posible obtener solución para sistemas con geometría muy sencilla, y/o con
condiciones de contorno simples. También se han utilizado técnicas de diferencias finitas, pero
éstas plantean problemas cuando los contornos son complicados.
CRITERIOS DE CONVERGENCIA
Antes de estudiar los criterios para garantizar la convergencia en el MEF es necesario definir
dicho concepto, en el ámbito del MEF. Se dice que un análisis por el MEF es convergente si al
disminuir el tamaño de los elementos, y por lo tanto aumentar el número de nudos y de
elementos, la solución obtenida tiende hacia la solución exacta. Hay que indicar que en el
análisis por el MEF, se introducen, además de la hipótesis de discretización, otras
aproximaciones, que son fuentes de error en la solución: integración numérica, errores de
redondeo por aritmética finita... El concepto de convergencia aquí analizado se refiere
solamente a la hipótesis de discretización, prescindiendo de los otros errores, que deben ser
estudiados aparte, y cuyo valor debe en todo caso acotarse. Las funciones de interpolación
elegidas para representar el estado de deformación de un medio continuo deben satisfacer
una serie de condiciones, a fin de que la solución obtenida por el MEF, converja hacia la
solución real.
Criterio 1 : Las funciones de interpolación deben ser tales que cuando los desplazamientos de
los nudos del elemento correspondan a un movimiento de sólido rígido, no aparezcan
tensiones en el elemento. Este criterio se puede enunciar también de forma más sencilla: las
funciones de interpolación deben ser capaces de representar los desplazamientos como sólido
rígido, sin producir tensiones en el elemento. Esta condición es evidente, pues todo sólido que
se desplaza como un sólido rígido, no sufre ninguna deformación ni por lo tanto tensión. Sin
embargo adoptando unas funciones de interpolación incorrectas, pueden originarse tensiones
al moverse como sólido rígido.
Criterio 2 : Las funciones de interpolación deben ser tales que cuando los desplazamientos de
los nudos correspondan a un estado de tensión constante, este estado tensional se alcance en
realidad en el elemento. Claramente, a medida que los elementos se hacen más pequeños, el
estado de tensiones que hay en ellos se acerca al estado uniforme de tensiones. Este criterio lo
10 Introducción que exige es que los elementos sean capaces de representar dicho estado de
tensión constante. Se observa que este criterio de hecho es un caso particular del criterio 1, ya
que un movimiento como sólido rígido (con tensión nula) es un caso particular de un estado de
tensión constante. En muchas ocasiones ambos criterios se formulan como un sólo criterio. A
los elementos que satisfacen los criterios 1 y 2 se les llama elementos completos.
Criterio 3 : Las funciones de interpolación deben ser tales que las deformaciones unitarias que
se produzcan en las uniones entre elementos deben ser finitas. Esto es lo mismo que decir que
debe existir continuidad de desplazamientos en la unión entre elementos aunque puede haber
discontinuidad en las deformaciones unitarias (y por lo tanto en las tensiones, que son
proporcionales a ellas).
ECUACIONES GENERALES
CAMPO DE DEFORMACIONES
El campo de deformaciones en un punto cualquiera del dominio está definido por un vector u
que tiene tantas componentes como deformaciones existen en el dominio.
La matriz de funciones de interpolación N tiene tres filas y tantas columnas como grados de
libertad haya entre todos los nudos del elemento. La estructura de esta matriz siempre es del
tipo:
DEFORMACIONES UNITARIAS
Las deformaciones unitarias en un punto cualquiera del elemento, con la suposición de

pequeñas deformaciones, son:
Se pueden poner en la forma matricial siguiente:
En esta expresión se identifica el operador matricial ∂ que permite pasar de las deformaciones
de un punto u a las deformaciones unitarias ε. Este operador tiene tantas filas como
deformaciones unitarias haya en el problema y tantas columnas como componentes tenga el
campo de desplazamientos u.
BIBLIOGRAFIA
http://bdigital.unal.edu.co/10002/6/958932276X.2002.pdf
http://ocw.uc3m.es/mecanica-de-medios-continuos-y-teoria-de-estructuras/ingenieria-
estructural/material-de-clase-1/apuntes/Capitulo_9.-
Introduccion_al_metodo_de_los_elementos_finitos.pdf
https://www.academia.edu/7224640/Metodo_de_los_elementos_finitos_para_Analisis_Estru
ctural

Metodo de Los Elementos Finitos

Cargado por

Copyright:

Formatos disponibles

Metodo de Los Elementos Finitos

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Metodo de Los Elementos Finitos

Cargado por

Copyright:

Formatos disponibles

2019

METODO DE LOS ELEMENTOS FINITOS

El Método de los Elementos Finitos es un método numérico de resolución de problemas de

El elemento finito lineal (barra biarticulada): planteamiento en coordenadas cartesianas y

Denominemos u1 y u2 a los desplazamientos, según el eje x, que experimentan sus nudos y

Elemento finito triangular de deformación constante: Matriz de rigidez

Si fuéramos capaces de dividir un sólido, considerado éste como un medio continuo, en

SISTEMAS DISCRETOS Y SISTEMAS CONTINUOS

Al efectuar una clasificación de las estructuras, suelen dividirse en discretas o reticulares y

Las deformaciones unitarias en un punto cualquiera del elemento, con la suposición de

También podría gustarte