Método Elemento Finito Cap 1

EL MÉTODO DEL ELEMENTO
FINITO
Ing. Rodolfo Alejandro García Pérez

CAPÍTULO 1. CONCEPTOS BÁSICOS
El método de elemento finito es un método numérico general para la aproximación de
soluciones de ecuaciones diferenciales (ordinarias y parciales). Muy utilizado en diversos
problemas principalmente de ingeniería y física – matemática aplicada. Dada la
imposibilidad práctica de encontrar una solución analítica de estos problemas, con
frecuencia, en la práctica ingenieril los métodos numéricos y, en particular, los elementos
finitos, se convierten en la única alternativa práctica de cálculo.
Una propiedad importante del método es la convergencia; si se consideran particiones de
elementos finitos sucesivos más finas, la solución numérica calculada converge rápidamente
hacía la solución exacta del sistema de ecuaciones.
1.1. INTRODUCCIÓN
El MEF o AEF (análisis de elemento finito), es una técnica computacional usada para obtener
soluciones aproximadas de problemas con valores límite en diversos campos de ingeniería.
Declarado de manera simple, los problemas con valores límites son problemas matemáticos
en los cuales una o más variables dependientes deben satisfacer una ecuación diferencial en
todas partes dentro de un domino conocido de variables independientes y satisfacer
condiciones específicas en el límite del dominio del problema. Los problemas de valores
límites también se llaman a veces problemas de campo. El campo es el domino de interés y
a menudo representa una estructura física. Las variables de campo son las variables
dependientes de interés gobernadas por la ecuación diferencial. Las condiciones límite son
los valores especificados de las variables de campo (o variables relacionadas con sus
derivadas) sobre los límites de campo. Dependiendo del tipo del problema físico que se
analiza, las variables de campo pueden incluir desplazamiento físico, temperatura, flujo de
calor, y velocidad de flujo, para mencionar solamente algunas variables que se presentan.
1.2. CONCEPTOS DEL MÉTODO DE ELEMENTO FINITO
Se introducirán las técnicas y terminología generales del análisis del elemento finito con
referencia a la figura 1.1, desarrollando el concepto básico del método. La figura mencionada
muestra el volumen de algún material o materiales que poseen propiedades físicas conocidas.
El volumen representa el dominio del problema de valor límite al ser resuelto Por
simplicidad, en este punto asumimos un caso bi – dimensional con una sola variable de
campo  ( x, y ) para ser determinada en un punto P( x, y ) tal que una ecuación gobernante
conocida (o ecuaciones) se satisfagan exactamente en cada uno de tales puntos pertenecientes
al dominio espacial de la variable  en cuestión. Note que esto implica que se obtiene una
solución matemática exacta; es decir, la solución es una expresión algebraica de forma –
cerrada de las variables independientes y tal vez algunos parámetros del dominio.
En los problemas prácticos, el dominio puede ser geométricamente complejo, como sucede
a menudo; como también resulta engorrosa la ecuación gobernante del problema y es muy
baja la probabilidad de obtener una solución exacta en forma – cerrada para la situación así
planteada. Por consiguiente, soluciones aproximadas basadas en técnicas numéricas y el
cómputo digital son más a menudo obtenidas en análisis de ingeniería para los problemas
complejos. El análisis de elemento finito es una técnica basada en un proceso de
discretización del dominio de definición del problema, que es muy poderosa para obtener
tales soluciones aproximadas con una exactitud muy buena cuando la discretización del
dominio es refinada.
Un pequeño elemento triangular que encierra un sub – dominio de tamaño finito del área del
interés se muestra en la figura 1.1.(b). Que este elemento no un elemento diferencial de
tamaño dx  dy hace del mismo un elemento finito; es decir, un elemento de dimensiones
finitas (no infinitesimales). Cuando nosotros tratamos este ejemplo como un problema bi –
dimensional, se supone que el espesor en la dirección z es constante y la dependencia según
esta dirección no se indica en la ecuación diferencial. Los vértices del elemento triangular se
numeran para indicar que estos puntos son los nodos del elemento finito. Un nodo es un punto
específico perteneciente al elemento finito en el que el valor de la variable de campo va a ser
calculado explícitamente. Los nodos exteriores son localizados en los límites del elemento
finito y podrían usarse para conectar un elemento dado a elementos finitos adyacentes. Los
nodos que no se ubican en los limites del elemento son nodos interiores y no pueden ser
conectados a los nodos de cualquier otro elemento. El elemento triangular de la figura 1.1 (b)
tiene sólo nodos en sus vértices, y por ello los mismos son nodos exteriores.
Si sólo se computan los valores de la variable de campo en los nodos del elemento ¿cómo se
obtienen los valores en otros puntos dentro de un elemento finito? La respuesta a esta crucial
pregunta contiene la esencia del método de elemento finto. Los valores de la variable de
campo calculados para los nodos se usan para aproximar los valores en los puntos no –
nodales (es decir, en todos los puntos al interior del elemento) mediante un procedimiento de
interpolación usando los valores nodales ya calculados en forma previa. Para el ejemplo del
elemento triangular con nodos en sus vértices, en el que los nodos son todos exteriores, la
variable de cualquier otro punto dentro del elemento y en sus aristas limitantes, se describe
por la relación aproximada.
 ( x, y) = N1 ( x, y)1 + N2 ( x, y)2 + N3 ( x, y)3
1 , 2 , 3 → Valores de variables de campo en los puntos nodales.

N1, N 2 , N3 → Funciones de interpolación, funciones de forma o funciones de mezcla.
En el ámbito del método del elemento finito, los valores nodales de la variable de campo se
tratan se tratan como contantes desconocidas que serán determinadas. Las funciones de
interpolación son más a menudo formas polinómicas de las variables independientes,
deducidas para satisfacer ciertas condiciones requeridas en los nodos. Estas condiciones son
discutidas en detalle en los capítulos subsecuentes. El punto relevante aquí es que las
funciones de interpolación son predeterminadas, con formas de funciones conocidas de las
variables independientes; y estas funciones describen la variación de la variable de campo
dentro del elemento finito.
Se dice que el elemento triangular descrito por la Ecuación (1.1) tiene 3 grados de libertad,
en virtud de que se requieren tres valores nodales de la variable de campo para describir a la
variable dependiente en todos los puntos que conforma el elemento. Éste sería el caso si la
variable de campo representa una cantidad escala, como la temperatura en un problema de
transmisión de calor. Si el dominio de Figura 1.1 representa un cuerpo sólido delgado, sujeto
a un estado de tensión plana; la variable de campo se adopta como el vector desplazamiento
interno y en este caso deben computarse los valores de dos componentes en cada nodo. En
este último caso, el elemento triangular con nodos en sus tres vértices tiene 6 grados de
libertad. En general, el número de grados de libertad asociados con un elemento finito es
igual al producto del número de nodos y el número de valores del campo variable (y
posiblemente sus derivadas) que deben calcularse en cada nodo.
En esta instancia surge otra pregunta ¿cómo este procedimiento basado en el elemento, puede
ser ahora aplicado sobre el dominio entero de interés? Cada elemento se conecta mediante
sus nodos exteriores a los otros elementos circundantes a él. Las ecuaciones de elemento
finito se formulan de modo que, en las conexiones nodales, el valor de la variable en cualquier
punto de conexión común es el mismo para cada elemento que concurre a dicha ubicación y
está conectado al nodo. Así, se asegura la continuidad de la variable de campo (variable
dependiente en el problema) en los puntos nodales. De hecho, las formulaciones de elemento
finito son tales que la continuidad del campo variable lo largo de los límites de elementos
aledaños que comparten dicho límite también se asegura. Este rasgo de imposición sobre el
método en cuanto a la continuidad de la variable que está siendo modelada y aproximada
mediante el método de elemento finito, evita la posibilidad físicamente inaceptable de la
existencia de huecos o discontinuidades en el dominio considerado.
Aunque la continuidad del campo variable de elemento a elemento adyacente es inherente a
la formulación del elemento finito, la continuidad de los gradientes (derivadas) del campo de
la variable dependiente generalmente no existe. Ésta es una observación crítica hacía el
método del elemento finito. En la mayoría de los casos, dichas derivadas son de más interés
que los valores del campo variable en sí – mismo. Por ejemplo, en los problemas estructurales
la variable de campo es el desplazamiento, pero el verdadero interés está más a menudo en
la deformación y la tensión. Como la deformación unitaria (al igual que la tensión interna)
está definida en términos de las primeras derivadas espaciales de las componentes de
desplazamiento, por lo que se refiere a esta variable mecánica no es continúa sobre la frontera
límite del elemento (exceptuando sus puntos nodales exteriores). Sin embargo, las
magnitudes de las discontinuidades de las derivadas de la función que se esté manipulando
pueden ser usadas para evaluar la exactitud de la solución y su convergencia hacia la solución
real (o más exacta) a medida que el número de elementos se incrementa.
1.2.1 Convergencia del método
El proceso de representar un dominio físico con elementos finitos es referido como un
proceso de discretización, y la serie resultante de elementos es conocida como la malla de
elemento finito.
Si las funciones de interpolación satisfacen ciertos requisitos matemáticos, una solución de

elemento finito para un problema particular converge hacia la solución exacta del problema
a medida que el modelo utilizado sea asociado a una malla de elementos de menor tamaño.
Los cambios incrementales decrecen con el proceso de refinamiento de malla y la
aproximación hacía la solución exacta se verifica que tiene forma asintótica.
Como se verá más adelante, el desplazamiento es la variable de campo primaria en los
problemas estructurales. En la mayoría de los problemas, sin embargo, estamos interesados
principalmente en las tensiones inducidas por las cargas especificadas. Las tensiones deben
ser calculadas a través de las relaciones tensión – deformaciones apropiadas (ley de Hooke
generalizada), y las componentes de deformación unitaria se obtienen de la solución del
campo de desplazamiento. De aquí, tensiones y deformaciones se refieren, por guardar esta
dependencia como variables derivadas.
La formulación del método del elemento finito para un problema dado es tal que la variable
de campo primaria es continua de elemento al elemento adyacente, pero las variables
derivadas (como la tensión interna, por ejemplo) no necesariamente poseen continuidad. En
el proceso límite de refinamiento de la malla, las variables derivadas se ponen cada vez más
cercanas al estado de continuidad si el número de elementos de la malla se incrementa.
¿Cómo aseguramos nosotros la exactitud de una solución obtenida mediante el método de

elemento finito para un problema con una solución teórica desconocida? La respuesta a esta
pregunta no es simple.
La convergencia de un método numérico (como el método del elemento finito) no está dada
por ninguna convicción de que los medios utilizados indiquen que la convergencia es
precisamente hacia la solución correcta (porque podría suceder que la aproximación sea hacia
un valor erróneo). Una persona que usa la técnica de análisis de elemento finito, debe
examinar la solución muy metódicamente por lo que se refiere a: (1) la convergencia
numérica, (2) la racionalidad (¿el resultado obtenido tiene sentido?, (3) si las leyes físicas del
problema están satisfechas (¿la estructura está en equilibrio? ¿existe balance entre el calor de
entrada y el calor de salida?, y (4) si los valores de las discontinuidades de las variables
derivadas en los límites inter – elementos son razonables. Deben proponerse muchas de tales
preguntas y deben examinarse prioritariamente las respuestas antes de aceptar los resultados
de un análisis de elemento finito como representativos de una solución correcta y útil para
los propósitos del diseño.
1.2.2. El método de diferencias finitas
Presentamos los conceptos básicos del método de diferencias finitas por propósitos de
comparación con el método de elemento finito.
El método de diferencias finitas está basado en la definición de la derivada de una función
cualquiera f (x), esto es
df ( x) f ( x + x) − f ( x)
= lim (1.2)
dx x →0 x
La solución de diferencias finitas sólo es mostrada en los puntos discretos de evaluación de
la función. La manera de variación de la función entre los puntos calculados no es conocida
en el método de diferencias finitas; pero, mostramos una curva de ajuste en trazo segmentado
que pasa a través de los puntos calculados. Uno puede, sin embargo, interpolar linealmente
los valores hallados para producir una aproximación a la curva de la solución exacta; pero,
la manera de efectuar este proceso de interpolación no es algo relevante, ni de determinación
“a priori” en el método de diferencias finitas.
Para contrastar el método de diferencias finitas con el método de elemento finito, notamos
que, en el método de elemento finito, la variación de la variable de campo en el dominio
físico es una parte integral del procedimiento.
En un problema estructural, por ejemplo, ambos métodos proporcionan las soluciones del
desplazamiento, pero la solución de elemento finito puede usarse para calcular las
componentes de tensión directamente (en términos de la primera derivada). Para obtener
datos de tensión en el método de diferencias finitas se requieren consideraciones adicionales
no inherentes al modelo matemático mismo.
Hay también ciertas similitudes entre los dos métodos. Los puntos de integración en el
método de diferencias finitas son análogos a los nodos en un modelo de elemento finito.
1.3. Procedimiento General
Ciertos pasos en la formulación de un análisis de elemento finito de un problema físico son
comunes a todos los análisis de característica discreta, sea este estructural, transmisión de
calor, flujo fluido, o algún otro problema
1.3.1. Pre – procesamiento
El paso de pre – procesamiento es, muy generalmente, descrito como el proceso de definir
el modelo e incluye las siguientes acciones a ejecutar.
• Definir el dominio geométrico del problema.

• Definir el tipo de elemento a ser utilizado.
• Definir las propiedades materiales de los elementos.
• Definir las propiedades geométricas de los elementos.
• Definir las conectividades de los elementos.
• Definir las restricciones físicas.
• Definir la solicitación o perturbación externa.
El paso de pre–procesamiento (la definición del modelo) es realmente crítico. Una solución
de elemento finito perfectamente calculada numéricamente no posee absolutamente ningún
valor si corresponde a un pésimo proceso de modelado, o a un problema erróneamente
planteado.
1.3.2. Procesamiento o solución
Durante la fase de solución, el software de elemento finito ensambla las ecuaciones
algebraicas gobernantes en forma matricial y calcula los valores desconocidos de las
variables primarias del campo involucrado en el problema. Los valores computados se usan
luego por substitución regresiva para evaluar las variables adicionales, derivadas o
secundarias, como las fuerzas de reacción, tensiones del elemento, flujo de calor, velocidades
y aceleraciones de movimiento, etc.
1.3.3. Post – procesamiento
El análisis y la evaluación de los resultados de la solución es denominado como post
procesamiento.
• Clasificar las tensiones internas del elemento en orden de magnitud.

• Verificar el equilibrio estático estructural.
• Calcular factores de seguridad de diseño.
• Bosquejar la deformación estructural producida.
• Producir esquemas dinámicos de la respuesta del modelo.
• Producir gráficas de color – codificado para el campo de temperaturas.
El objetivo más importante es aplicar de modo legítimo el juicio de ingeniería con el objetivo
de determinar las condiciones bajo las cuales los resultados de la solución obtenida son
físicamente razonables.
CAPÍTULO 2
ANÁLISIS UNI – DIMENSIONAL

Método Elemento Finito Cap 1

Cargado por

Copyright:

Formatos disponibles

Método Elemento Finito Cap 1

Cargado por

Información del documento

Descripción original:

Título 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

Método Elemento Finito Cap 1

Cargado por

Copyright:

Formatos disponibles

EL MÉTODO DEL ELEMENTO

Ing. Rodolfo Alejandro García Pérez

 ( x, y) = N1 ( x, y)1 + N2 ( x, y)2 + N3 ( x, y)3

1 , 2 , 3 → Valores de variables de campo en los puntos nodales.

Si las funciones de interpolación satisfacen ciertos requisitos matemáticos, una solución de

¿Cómo aseguramos nosotros la exactitud de una solución obtenida mediante el método de

• Definir el dominio geométrico del problema.

• Clasificar las tensiones internas del elemento en orden de magnitud.

También podría gustarte