CAPTULO 8

CONCEPTOS BSICOS DEL MTODO POR

ELEMENTO FINITO

8.1 GENERALIDADES
El mtodo del elemento finito (MEF en espaol o FEM en ingls) es un mtodo numrico
para la resolucin de ecuaciones diferenciales, utilizado en diversos problemas de
ingeniera y fsica.
El mtodo se basa en dividir el cuerpo, estructura o dominio (medio continuo) sobre el que
estn definidas ciertas ecuaciones integrales que caracterizan el comportamiento fsico
del problema (figura 8.1), en una serie de subdominios no intersectantes entre s
denominados elementos finitos. El conjunto de elementos finitos forma una particin del
dominio tambin llamada discretizacin.

FIGURA 8.1 Ejemplos de discretizacin

Dentro de cada elemento se distinguen una serie de puntos representativos llamados

nodos. Dos nodos son adyacentes s pertenecen al mismo elemento finito, adems, un
nodo sobre la frontera de un elemento finito puede pertenecer a varios elementos (figura
8.2). El conjunto de nodos considerando sus relaciones de adyacencia se conoce como

MODELADO DE PROCESOS DE MANUFACTURA

malla. Los clculos se realizan sobre una malla o discretizacin creada a partir del
dominio con programas generadores de mallas, en una etapa previa a los clculos que se
denomina pre-proceso. De acuerdo con estas relaciones de adyacencia o conectividad se
relaciona el valor de un conjunto de variables incgnitas definidas en cada nodo y
denominadas grados de libertad. El conjunto de relaciones entre el valor de una
determinada variable entre los nodos se puede escribir en forma de sistema de
ecuaciones lineales (o linealizadas), la matriz de dicho sistema de ecuaciones se llama
matriz de rigidez del sistema. El nmero de ecuaciones de dicho sistema es proporcional
al nmero de nodos.
Tpicamente, el mtodo del elemento finito se programa computacionalmente para
calcular el campo de desplazamientos y, posteriormente, a travs de relaciones
cinemticas y constitutivas, las deformaciones y tensiones respectivamente, cuando se
trata de un problema de mecnica de slidos deformables o ms generalmente un
problema de mecnica del medio continuo. El mtodo de los elementos finitos es muy
usado debido a su generalidad y a la facilidad de introducir dominios de clculo complejos
(en dos o tres dimensiones).

FIGURA 8.2 Discretizacin para aplicaciones geolgicas

Adems, el mtodo es fcilmente adaptable a problemas de difusin del calor, de

mecnica de fluidos para calcular campos de velocidades y presiones, o de campo
electromagntico. Dada la imposibilidad prctica de encontrar la solucin analtica de
estos problemas, con frecuencia, en la prctica ingenieril, los mtodos numricos y, en
particular, los elementos finitos se convierten en la nica alternativa prctica de clculo.

234

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

Una importante propiedad del mtodo es la convergencia, si se consideran particiones de

elementos finitos sucesivamente ms finas, la solucin numrica calculada converge
rpidamente hacia la solucin exacta del sistema de ecuaciones

8.2 HISTORIA
Diversos autores han considerado que Arqumedes (figura 8.3) utiliz un mtodo
semejante al del elemento finito para determinar el volumen de algunos slidos. Aunque l
calcul reas, longitudes y volmenes de objetos geomtricos, dividindolos en otros ms
sencillos y luego sumando sus contribuciones, el concepto de aproximacin variacional no
se observa por ningn lado. La relacin con la definicin de MEF es muy pobre. Se puede
argumentar que la medida del volumen (rea, longitud) de un objeto es una funcin
escalar de su geometra. Cambiando medida por energa y objetos por elementos en
las lneas anteriores, la descripcin se aproxima a lo establecido por el MEF la energa
del sistema es igual a la suma de la energa de cada elemento. Sin embargo,
Arqumedes necesitaba las definiciones de derivada para realizar sus clculos de energa
y el Clculo no fue inventado sino hasta 20 siglos despus.

FIGURA 8.3 Arqumedes de Siracusa

En 1941, Hrenikoff present una solucin para problemas elsticos usando el mtodo de
trabajo del marco. En un artculo publicado en 1943, Courant us interpolacin polinomial
por partes sobre subregiones triangulares para modelar problemas de torsin. Las ideas
bsicas del mtodo del elemento finito se originaron en el anlisis estructural de las
aeronaves. En el periodo de 1950-1962, Turner trabajando para Boeing formula y

235

MODELADO DE PROCESOS DE MANUFACTURA

perfecciona el Mtodo por Rigidez Directo. Turner y otros investigadores obtuvieron

matrices de rigidez para armaduras, vigas y otros elementos y presentaron sus resultados
en 1956. Clough fue el primero en acuar y emplear el trmino elemento finito en 1960.

En los primeros aos de la dcada de 1960, los ingenieros usaron el mtodo para obtener
soluciones aproximadas en problemas de anlisis de esfuerzos, flujo de fluidos,
transferencia de calor y otras reas. Un libro de Argyris, publicado en 1955, sobre
teoremas de energa y mtodos matriciales, ciment mtodos adicionales en los estudios
de elemento finito. El primer libro sobre elementos finitos por Zienkiewicz y Cheng fue
publicado en 1967. A finales de la dcada de 1960 y principios de la siguiente, el anlisis
por elemento finito se aplic a problemas no lineales y de grandes deformaciones. El libro
de Oden sobre continuos no lineales apareci en 1972.

8.3 DELIMITACIN DEL MTODO

El campo de la mecnica puede ser subdividido en tres reas (figura 8.4). La mecnica
terica estudia las leyes y principios fundamentales de la mecnica, pero considerando
solo su valor cientfico. La mecnica aplicada transfiere el conocimiento terico hacia
aplicaciones cientficas y de ingeniera, atendiendo principalmente la construccin de
modelos matemticos que representen fenmenos fsicos. La mecnica computacional
resuelve problemas especficos aplicando mtodos numricos implementados en
computadoras digitales (simulacin).

FIGURA 8.4 Ramas de la Mecnica

Mecnica computacional
Varias ramas de la mecnica computacional pueden distinguirse de acuerdo a la escala
fsica en que se enfocan (figura 8.5)

236

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

FIGURA 8.5 Ramas de la mecnica computacional

La nanomecnica trata con fenmenos de la materia a nivel molecular y atmico. Por sta
razn se encuentra ligada a la fsica y qumica de las partculas. La micromecnica, por su
parte, trabaja a nivel cristalogrfico y granular de la materia. Su principal aplicacin
tecnolgica es el diseo y fabricacin de materiales y microdispositivos.

La mecnica del continuo estudia los cuerpos a nivel macroscpico, utilizando modelos
continuos en los cuales la microestructura es considerada como homognea gracias a
promedios cualitativos. Las reas tradicionales de aplicacin son la mecnica de slidos y
de fluidos. Se incluyen las estructuras debido a que, obviamente, son fabricadas con
slidos. La mecnica computacional de slidos usa aproximaciones de ciencias aplicadas,
mientras que la mecnica estructural computacional hace nfasis en aplicaciones
tecnolgicas para el anlisis y diseo de estructuras.

La mecnica computacional de fluidos estudia problemas que involucran el movimiento y

equilibrio de lquidos y gases. Las reas beneficiadas con sus desarrollos tecnolgicos
son la hidrodinmica, aerodinmica, acstica, fsica de la atmsfera, combustin y
propulsin.

Finalmente, los sistemas identifican objetos mecnicos, naturales o artificiales, que

realizan una funcin fcilmente reconocible. Ejemplos de sistemas artificiales son todos
aquellos dispositivos hechos por el hombre. Los dems sistemas en los que el hombre no
tiene participacin se consideran como naturales.

237

MODELADO DE PROCESOS DE MANUFACTURA

Esttica vs. Dinmica

Los problemas del continuo pueden ser subdivididos de acuerdo a si sus efectos
inerciales son considerados o no (figura 8.6). En dinmica, la dependencia del tiempo
actual debe ser considerada explcitamente, porque el clculo de fuerzas inerciales
requiere tomar derivadas con respecto al tiempo actual.

FIGURA 8.6 Ramas de la mecnica del continuo

Los problemas en esttica tambin pueden depender del tiempo actual pero ignoran o
desprecian las fuerzas inerciales. Por consiguiente, los problemas estticos pueden ser
clasificados en estrictamente estticos y cuasiestticos. En los primeros, el tiempo no
necesita ser considerado explcitamente; cualquier tiempo registrado como respuesta, se
usa como parmetro de orden, si es que alguno se necesita. En los problemas
cuasiestticos, tales como la edificacin de cimientos, fluencia en metal, plasticidad o
fatiga cclica, una medida real es requerida pero las fuerzas inerciales an son
despreciadas.

Lineal vs. no lineal

Una clasificacin de los problemas estticos, particularmente relevante para estas notas,
se muestra en la figura 8.7. El anlisis esttico lineal involucra problemas estticos en los
cuales la respuesta es lineal en el sentido causa y efecto. Por ejemplo, si las fuerzas
aplicadas se duplican, los desplazamientos y esfuerzos internos tambin se duplican. Los
problemas que caen fuera de este dominio son clasificados como no lineales.

FIGURA 8.7 Divisin de problemas estticos

238

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

Mtodos de discretizacin
Finalmente, se pueden clasificar los anlisis estticos basndose en el mtodo de
discretizacin (figura 8.8) por medio del cual el modelo matemtico del continuo es
discretizado en el espacio, es decir, convertirlo en un modelo discreto con un nmero
finito de grados de libertad:

FIGURA 8.8 Mtodos de discretizacin

En los problemas lineales, el mtodo del elemento finito domina actualmente la escena en
lo relativo a discretizacin espacial. Los mtodos por fronteras finitas constituyen una
segunda alternativa en reas de aplicacin especficas. Para problemas no lineales, el
dominio de los mtodos por elemento finito es enorme.
Los mtodos por diferencias finitas aplicados en slidos y mecnica de estructuras han
desaparecido virtualmente debido a que son poco prcticos. Sin embargo, para mecnica
de fluidos estos mtodos son an importantes. Los mtodos por volmenes finitos, que se
relacionan directamente con la discretizacin de las leyes de conservacin, son
importantes en problemas difciles de la mecnica de fluidos; por ejemplo, en dinmica de
gases con nmero de Reynolds elevado.
Los mtodos espectrales se basan en correspondencias que transforman dimensiones
espaciales y/o temporales a espacios (por ejemplo, el dominio de frecuencias), donde el
problema es mucho ms sencillo de resolver.
De reciente aparicin se tienen los mtodos por mallado automtico. Estos combinan
tcnicas y herramientas del elemento finito, por ejemplo la formulacin variacional e
interpolacin, con caractersticas de diferencias finitas como el soporte no localizado.

239

MODELADO DE PROCESOS DE MANUFACTURA

Variantes del mtodo por elemento finito

Actualmente el trmino MEF identifica a una amplia gama de tcnicas (figura 8.9) que
comparten caractersticas descritas en los puntos anteriores. Dos subclasificaciones que
se ajustan a la mecnica de estructuras son:

FIGURA 8.9 Variantes del Elemento Finito

Utilizando la clasificacin anterior, podemos establecer el tema principal de estas notas:

Anlisis computacional por elemento finito de problemas estticos lineales. De las
variantes mostradas, se pondr especial atencin a la formulacin por desplazamiento y a
la solucin por rigidez. Esta combinacin recibe el nombre de Mtodo directo por rigidez o
DSM.

8.4 APLICACIONES DEL ELEMENTO FINITO

El objeto de estas notas es estudiar el Mtodo del elemento finito, pero qu es un
elemento finito? El concepto bsico puede ser parcialmente ilustrado a travs de un
antiguo problema: encontrar el permetro L de un crculo cuyo dimetro es d . Como

L d , esto equivale a obtener un valor numrico para . Se dibuja un crculo de radio

r y dimetro d 2r como se muestra en la figura 8.10a. Se inscribe un polgono regular
de n lados, donde n 8 en la figura 8.10b. Se renombran los lados del polgono como
elementos y los vrtices como nodos. Las etiquetas de los nodos son enteros que van de
1 a 8. Considrese un elemento tpico, el que une los nodos 4-5, como se muestra en la
240

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

figura 8.10c. Este es un caso del elemento genrico i j mostrado en la figura 8.10d. La
longitud del elemento es

Lij 2r sen( / n) 2r sen(180 / n) . Como todos los

elementos tienen la misma longitud, el permetro del polgono es Ln nLij , de donde la

aproximacin para resulta n Ln / d n sen(180 / n) .

FIGURA 8.10 El problema de encontrar el nmero tratado con conceptos de MEF: (a) objeto
continuo, (b) aproximacin discreta por polgonos regulares inscritos, (c) elemento
desconectado, (d) elemento genrico

Para facilitar la explicacin se ha implementado un sencillo algoritmo en Mathematica,

mismo que se anexa en las siguientes pginas (Archimedes_nb).

Algunas ideas del MEF, pueden identificarse gracias al ejemplo anterior. El crculo, visto
como un objeto matemtico, es reemplazado por polgonos. Estos constituyen la
aproximacin discreta del crculo. Los lados, renombrados como elementos, estn
completamente identificados por los nodos en sus extremos. Los elementos pueden
separarse desconectando sus nodos, un proceso llamado desensamble en el MEF.
Gracias a este proceso, un elemento genrico puede ser definido, independientemente
del crculo original, por el segmento que conecta dos nodos i y j . La propiedad relevante
del elemento, en este ejemplo, es la longitud de su lado Lij , misma que puede ser
calculada en el elemento genrico independientemente de los otros, una propiedad
conocida como soporte local en el MEF. La propiedad objetivo: el permetro del polgono,
es obtenido al reconectar n elementos y sumando su longitud; los pasos correspondientes
en el MEF son el ensamble y la solucin, respectivamente. Por supuesto que no existe
nada particular en el problema del crculo, pues la misma tcnica se puede utilizar para
obtener la longitud de curvas suaves.

241

MODELADO DE PROCESOS DE MANUFACTURA

8.5 TCNICAS PARA FORMULAR ECUACIONES ELEMENTALES

En general, hay diversos modos de enfocar el problema mediante elemento finito: (1)
formulacin directa, (2) formulacin mediante la energa potencial y (3) formulacin por
pesos residuales. Es importante notar que los pasos bsicos que se usan en el anlisis
por elementos finitos son los mismos, sin tener en cuenta cmo se genera el modelo.

242

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

Formulacin directa
De acuerdo con este mtodo, la matriz de rigidez y las ecuaciones de los elementos que
relacionan a las fuerzas nodales con los desplazamientos nodales son obtenidas usando
las condiciones de equilibrio de fuerzas para un elemento bsico. Este mtodo es muy
utilizado en problemas de elementos unidimensionales, debido a su sencillez de
implementacin.

243

MODELADO DE PROCESOS DE MANUFACTURA

Formulacin mediante pesos residuales

Este mtodo es usado para desarrollar las ecuaciones de los elementos cuando se conoce
la ecuacin diferencial que describe al fenmeno fsico. El mtodo ms usado es el de
Galerkin; ste permite obtener las mismas ecuaciones que en el mtodo de la energa, pero
existen caso especiales en donde la funcin de la energa potencial no es viable.

Formulacin usando la energa o trabajo

En el desarrollo de la matriz de rigidez y las ecuaciones de los elementos bidimensionales
y tridimensionales, es mucho ms sencillo aplicar el mtodo de la energa. El principio de
trabajo virtual (usando desplazamientos virtuales), el principio de energa potencial
mnima y el teorema de Castigliano son usados frecuentemente para encontrar las
ecuaciones de los elementos.
El principio de trabajo virtual es aplicable para cualquier tipo de material, el principio de
energa potencial mnima y el teorema de Castigliano son teoremas aplicables para
materiales elsticos. Sin embargo, el principio de trabajo virtual puede ser usado cuando
no existe funcin potencial. Los tres principios permiten obtener las mismas ecuaciones
elementales para materiales elsticos y lineales.

8.6 ETAPAS PARA SOLUCIONAR PROBLEMAS MEDIANTE MEF

Las siguientes etapas comprenden los pasos bsicos para formular problemas mediante
el mtodo del elemento finito.

Fase de preproceso
Crear y discretizar la solucin dominio en elementos finitos, esto es, subdividir el
problema en nodos y elementos.
Asumir una funcin forma a representar el comportamiento fsico de un elemento,
que es, una funcin continua aproximada que se asume para la solucin del
elemento.
Desarrollar las ecuaciones para el elemento
Armar los elementos a representar en el problema completo, construir la matriz
global de rigidez.
Aplicar condiciones de frontera, condiciones iniciales y cargas.
244

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

Fase de solucin
Resolver un conjunto de ecuaciones algebraicas lineales o no lineales simultneas
para obtener resultados globales; tal como valores de desplazamientos y
temperaturas en diferentes nodos.

Fase de postproceso
Con base en los resultados globales obtenidos, se puede volver en el proceso de
solucin para conocer los valores de desplazamiento o temperatura locales en
diversos nodos.
El siguiente ejemplo ilustra los pasos y el procedimiento a realizar en una formulacin
directa.

EJEMPLO 8.1. Sea un cuerpo en forma de trapecio (figura 8.11), con una seccin

transversal variable que soporta una carga P , segn lo mostrado en la figura. El cuerpo
est fijado en la base mayor y en el otro extremo, la base menor, lleva la carga P . El
ancho del cuerpo en la parte superior es w1 , y en la parte inferior es w2 ; su espesor es t
y su longitud es L . El mdulo de la elasticidad de la barra es E . Se debe determinar
cunto se deformar el cuerpo en varios puntos a lo largo de su longitud cuando se aplica
la carga P . No se tendr en cuenta el peso del cuerpo en el anlisis, si se asume que la
carga aplicada es considerablemente ms grande que el peso de la barra.

FIGURA 8.11 Dimensiones del trapecio

245

MODELADO DE PROCESOS DE MANUFACTURA

SOLUCIN

Fase de preprocesado
Se subdivide el trapecio en nodos y elementos, se simplificar el problema para resaltar
los pasos bsicos del anlisis por elemento finito, y se representar al mismo por medio
de un modelo que tiene 5 nodos y 4 elementos, como se indica en la figura 8.12.

FIGURA 8.12 Divisin en nodos y elementos

Luego se reconstruye al trapecio con cuatro subcuerpos rectangulares, cuyas secciones

sern el promedio de las secciones de los lmites de cada subcuerpo original. Se divide al
cuerpo en 4 subcuerpos, cada uno de ellos tendr una altura l , igual a L , de acuerdo
4
con la siguiente expresin:

w2 w1
wi w1
yi
L

Se calculan las sucesivas bases menores wi dentro del trapecio en funcin de y , luego
con el valor wi se obtienen las secciones a cada incremento de l , esto resulta en 5
secciones que se muestran en la figura 8.13.

246

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

FIGURA 8.13 Dimensiones de las bases menores

Se hace el promedio de estas dos secciones, una superior y otra inferior, luego se igualan
a la seccin promedio de un cuerpo rectangular de espesor t y el largo se determina
como L1 , L 2 , L3 y L 4 . Con esta simplificacin se logra transformar al cuerpo en forma
de trapecio (figura 8.14), en otro de geometra parecida pero formada de cuatro
subsecciones rectangulares, con 5 nodos y 4 elementos.

FIGURA 8.14 Trapecio aproximado por cuatro secciones rectangulares

247

MODELADO DE PROCESOS DE MANUFACTURA

8.7 ENERGA POTENCIAL Y EQUILIBRIO. MTODO DE RAYLEIGH - RITZ

En la mecnica de slidos, el problema consiste en determinar el campo de
desplazamientos del cuerpo (figura 8.15) que satisfaga las ecuaciones de equilibrio. Como
los esfuerzos estn relacionados con deformaciones unitarias que, a su vez, estn
relacionadas con desplazamientos, esto conduce a resolver ecuaciones diferenciales
parciales de segundo orden. A la solucin de este conjunto de ecuaciones se le llama una
solucin exacta. Tales soluciones exactas existen para geometras y condiciones de carga
simples, que pueden encontrarse con ayuda de la teora de la elasticidad. Para problemas
de geometras complejas y condiciones de frontera y de carga general, la obtencin de
tales soluciones es una tarea casi imposible. Los mtodos de solucin aproximada
usualmente emplean mtodos de energa potencial o con variacin, que imponen
condiciones menos estrictas sobre las funciones.

FIGURA 8.15 Cuerpo tridimensional

248

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

Resolucin numrica de sistemas lineales

La dificultad de resolver un sistema lineal radica en que, esencialmente, se debe invertir
una matriz. Una posible estrategia para simplificar los clculos consiste en sustituir la
matriz a invertir A por otra de fcil inversin y reescribir un sistema equivalente que se
resuelve de forma aproximada basndose en una solucin aproximada del sistema. El
esquema es el siguiente:
Dado el sistema Ax b escribir la matriz A en la forma A S P , donde la matriz S es
de fcil inversin. Entonces se puede reescribir el sistema como Sx b Px . Si x 0 es
una solucin aproximada, se busca la solucin del sistema Sx b Px0 , de fcil clculo:

x S 1 b Px0

(8.1)

En realidad, el resultado que se obtiene es una aproximacin x1 de la solucin del

sistema Ax b . A partir de ella se puede repetir el procedimiento y construir de forma
iterativa la sucesin

xk

que se espera converja a la solucin del sistema. La regla

iterativa ser

x k 1 S 1 b Px k

(8.2)

Y para conseguir la convergencia; si x es la solucin buscada, el error en el paso k 1

puede expresarse como

xk 1 x S 1P xk x S 1P x k x
y por tanto

x k 1 x S 1P

k 1

x0 x

(8.3)

por lo que si S 1P 1 el mtodo converge independientemente de la aproximacin de

partida x 0 .
Esta es una condicin suficiente para la convergencia, pero no es necesaria. Otra
condicin suficiente es que el autovalor mayor (en mdulo) de la matriz S 1P sea
estrictamente menor que 1.
A continuacin se muestra uno de los mtodos iterativos ms comunes y sencillos para
resolver sistemas lineales.
249

MODELADO DE PROCESOS DE MANUFACTURA

Mtodo de GaussJacobi
Consiste en despejar de la ecuacin i -sima la incgnita xi x x1, x2 ,..., xn , lo cual es
posible siempre que ningn elemento de la diagonal de A sea nulo (k , akk 0) . En
caso de que alguno lo fuera, se deben permutar las ecuaciones hasta conseguir una
diagonal libre de ceros. Con el sistema escrito en esta forma se utilizan las coordenadas
de la aproximacin x k en el lado derecho del sistema para obtener las coordenadas de

x k 1 en el lado izquierdo. Desarrollando para cada valor:

b a
a

x1k 1
x2k 1

1
b1 a12 x2k a13 x3k
a11
k
21 x1

a23 x3k

a1n xnk

a2 n xnk

22

xnk 1

1
bn an1x1k an 2 x2k
ann

ann 1xnk1

(8.4)

La descomposicin de A equivalente a este mtodo consiste en tomar como S la

diagonal de A y entonces P S A est formada por los elementos que no estn sobre
la diagonal (con signo cambiado) en su posicin y ceros en la diagonal.

EJEMPLO 8.2. Resuelva el siguiente sistema por el mtodo de Gauss Jacobi

8 2 0 x 8

2 8 2 y 4
0 2 8 z 8

SOLUCIN

La descomposicin A S P permite expresar el sistema como:

8 0 0 x 8 0 2 0 x

0 8 0 y 4 2 0 2 y
0 0 8 z 8 0 2 0 z

250

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

Y al invertir la matriz S , se logra la forma para el mtodo de Gauss Jacobi:

x 1 0
1 1
y 2 4
z
1 0

0 x

1 y
4

0 z

1
4

0
1
4

Si se inicia el mtodo con el vector x0 0, 0, 0 , las sucesivas iteraciones tienen por

34

resultado lo siguiente 0
3
4

1
1
2
1

1
1
4
1

34

0
3
4

63
64

0
63
64

1

0
1

que convergen a

Energa potencial
La energa potencial total de un cuerpo elstico se define como la suma de la energa
de deformacin unitaria total (U ) y el potencial de trabajo:
Energa Potencial Total = Energa de Deformacin + Potencial de Trabajo

WP

Para materiales elsticos lineales, la energa de deformacin por unidad de volumen en el

cuerpo es 12 . Para el cuerpo elstico que se muestra (figura 8.15), la energa de
deformacin unitaria total U est dada por

1
dV
2 V

siendo

[ x , y , z , yz , xz , xy ]

esfuerzos

[ x , y , z , yz , xz , xy ]

deformaciones

251

(8.5)

MODELADO DE PROCESOS DE MANUFACTURA

El potencial de trabajo WP est dado por

WP ufdV uTdS ui Pi
V

(8.6)

siendo

ui u , v, w

vector de desplazamientos

f f x , f y , f z

fuerza distribuida por unidad de volumen

T Tx , Ty , Tz

fuerza de cuerpo

Pi Px , Py , Pz

carga actuando en un punto i

El potencial total para el cuerpo elstico general mostrado en la figura 8.15 es

1
dV ufdV uTdS ui Pi
2 V
i
V
S

(8.7)

Se consideran sistemas conservativos, donde el potencial de trabajo es independiente de

la trayectoria, es decir, si el sistema se desplaza desde una configuracin dada y se trae
de regreso al estado inicial, las cargas efectan un trabajo nulo, independientemente de la
trayectoria. Entonces, el principio de la energa potencial se enuncia como sigue:

Principio de la energa potencial mnima

Para sistemas conservativos, de todos los campos de desplazamiento
cinemticamente admisibles, aquellos que corresponden a condiciones de
equilibrio extremizan la energa potencial total. Si la condicin es un mnimo,
el estado de equilibrio es estable.
Los desplazamientos cinemticamente admisibles son aquellos que satisfacen la
naturaleza de valor nico de los desplazamientos (compatibilidad) y las condiciones de
frontera. En los problemas donde los desplazamientos son las incgnitas, la
compatibilidad se satisface automticamente.
252

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

EJEMPLO 8.3. Se tiene un resorte (figura 8.16) cuya constante de elasticidad es k .

FIGURA 8.16 Resorte simple

Dos fuerzas externas f1 y f 2 estn aplicadas a los nodos 1 y 2, respectivamente. Los

desplazamientos de los nodos desde el estado no deformado son u1 y u2 . Considrese el
sentido positivo del eje x hacia la derecha.

Analizar la figura y determinar lo siguiente:

a. La energa potencial total
b. Las ecuaciones de equilibrio en forma matricial.

SOLUCIN

a. La energa potencial total se obtiene a partir de la expresin 8.7:

1
dV ufdV uTdS ui Pi
2 V
i
V
S

(8.7)

Y del anlisis se llega a la siguiente expresin:

1
dV ui Pi
2 V
i

(8.8)

Para calcular la energa asociada al resorte deformado, se utiliza la siguiente

grfica:

253

MODELADO DE PROCESOS DE MANUFACTURA

Ley de Hooke
rea bajo la curva
(representa energa)
Sustituyendo f en A

f k
1
A f
2
1
A k 2
2

Considrese que la deformacin del resorte es igual a la diferencia de

desplazamientos, es decir u2 u1 ; entonces se pueden sustituir las
expresiones obtenidas en la energa potencial total:

1
2
k u2 u1 f1u1 f 2u2
2

ENERGA POTENCIAL TOTAL

b. Para el equilibrio de este sistema, se necesita minimizar con respecto a u1 y u2

Las dos ecuaciones estn dadas por

0
ui

i 1, 2

(8.9)

que son en forma explcita

k u2 u1 f1 ku1 ku2 f1 0
u1

k u2 u1 f 2 ku1 ku2 f 2 0
u2

(8.10)

Estas ecuaciones de equilibrio pueden escribirse en la forma KU F , como

sigue:

k
k

k u1 f1

k u2 f 2

(8.11)

en donde K se conoce como matriz de rigidez elemental, y es la parte ms

importante del mtodo del elemento finito, pues al invertirla, y resolver
correctamente el sistema, se pueden conocer los desplazamientos, las
deformaciones, los esfuerzos, etc.

254

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

EJEMPLO 8.4.

Se tiene el siguiente sistema de resortes (figura 8.17)

FIGURA 8.17 Sistema de resortes

Analizar la figura y determinar lo siguiente:

a. La energa potencial total del sistema
b. Las ecuaciones de equilibrio en forma matricial

SOLUCIN

a) Para calcular la energa potencial total, es necesario manipular un poco la

expresin de la ley de Hooke. En la siguiente grfica se observa el
comportamiento de un resorte sometido a una carga de traccin

Ley de Hooke
rea bajo la curva
(representa energa)
Sustituyendo f en A

255

f k
1
A f
2
1
A k 2
2

MODELADO DE PROCESOS DE MANUFACTURA

Se aplica la ltima expresin en cada resorte y el resultado debe sumarse a la

energa generada por las cargas externas aplicadas en los nodos 1 y 3. De esta
forma, la energa potencial total es:

1
1
1
1
k112 k2 2 2 k3 32 k4 4 2 f1u1 f3u3
2
2
2
2

(8.12)

donde 1, 2 , 3 y 4 son los alargamientos de los cuatro resortes. Como

1 u1 u2 , 2 u2 u3 , 3 u3 u2 , 4 u5 u3 , se tiene

1
2 1
2 1
2
2 1
k1 u1 u2 k2 u2 u3 k3 u3 u2 k4 u5 u3 f1u1 f3u3
2
2
2
2
(8.13)

donde u1, u2 , u3 son los desplazamientos de los nodos 1, 2 y 3, respectivamente.

Para el equilibrio de este sistema con tres grados de libertad, se necesita

minimizar la funcin con respecto a u1, u2 , u3 . Las tres ecuaciones estn dadas
por

0
ui

i 1, 2,3

(8.14)

que son en forma explcita

k1 u1 u2 f1 0
u1

k1 u1 u2 k2 u2 u3 k3 u3 u2 0
u2

(8.15)

k3 u3 u2 k4 u5 u3 f3 0
u3

Esas ecuaciones de equilibrio pueden escribirse en la forma KU F , como sigue:

k1
k1
k k k k
1 1 2 3
0
k3

u1 f1
k3 u2 0
k3 k4 u3 f3
0

256

(8.16)

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

Por otra parte, s se escriben las ecuaciones del sistema considerando el equilibrio
de cada nodo por separado (figura 8.18).

FIGURA 8.18 Solucin considerando equilibrio nodal

resulta lo siguiente

k11 f1
k2 2 k11 k33 0

(8.17)

k3 3 k4 4 f3
que es precisamente el conjunto de ecuaciones representado anteriormente.

8.8 METODOLOGA PARA EL MODELADO NUMRICO

DE PROCESOS DE CORTE

En el proceso cognitivo de problemas complejos la aceptacin de modelos fsicos que

representen el nivel presente de conocimiento es muy importante. La complicacin
excesiva del modelo aumenta ligeramente el trabajo al definir los requisitos de entrada y
de hardware para obtener los resultados esperados, por lo que hace su efectiva utilizacin
imposible. De ah que el compromiso entre la desinformacin del problema y la
257

MODELADO DE PROCESOS DE MANUFACTURA

complicacin del modelo convertido sea esencial. En el siguiente paso se crea el modelo
matemtico entendido como una descripcin matemtica de las dependencias fsicas del
fenmeno observado. Esta descripcin se puede proporcionar mediante

ecuaciones

algebraicas, ecuaciones diferenciales ordinarias o parciales, ecuaciones diferenciales e

integrales o ecuaciones diferenciales con retardo. Estas ltimas cobran mayor relevancia
para el modelado de procesos donde debe considerarse la influencia que la historia del
proceso genera en el estado actual del sistema.

La siguiente etapa de modelado dinmico es el anlisis de las ecuaciones de movimiento

del sistema en estudio. El mtodo de anlisis depende del tipo de las ecuaciones del
modelo. Para ecuaciones diferenciales lineales se obtienen soluciones analticas. Sin
embargo, para los sistemas no-lineales ms simples no hay soluciones precisas, o la
determinacin de estas soluciones es muy inconveniente. En estos casos, se puede
aplicar una aproximacin mediante mtodos analticos con base en funciones
elementales, pero sobre todo los mtodos numricos predominan en el anlisis de
sistemas no-lineales.

En la ltima etapa del modelado, los parmetros para los que las caractersticas
dinmicas requeridas pueden ser obtenidas se seleccionan, a esto se le llama sntesis y
optimizacin. Aparte del anlisis terico mencionado anteriormente, se debe contar con
investigacin experimental adecuada. sta permite la identificacin de los parmetros
fsicos del modelo y la verificacin del modelo matemtico y los resultados obtenidos.

8.9 MTODO DE ELEMENTO FINITO APLICADO A PROCESOS

DE DEFORMACIN

La tcnica del elemento finito cuyo boom inici en la dcada de los sesenta y se logr
aplicar al clculo de anlisis estructuras mediante computadoras.

El concepto bsico del mtodo es la discretizacion. Un nmero de puntos finitos son

identificados en el dominio de la funcin, y el valor de la funcin y sus derivadas son
especificados cuando es conveniente. Estos puntos se conocen como nodos, el dominio

258

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

de la funcin es representada aproximadamente por una coleccin de subdominios

llamados elementos finitos.

Una vez que el cuerpo es dividido en elementos finitos, conectados mediante nodos, (este
proceso se conoce como mallado) se procede al ensamble de las ecuaciones en las
fronteras de cada elemento para que haya continuidad. En sistemas lineales los nodos
sern igual al nmero de variables o incgnitas, por lo tanto, entre ms fino sea el mallado
mayor cantidad de incgnitas se tendrn, de modo que es esencial el uso de una
computadora para resolverlas, en este caso son variables de campo.

El mtodo de elemento finito (MEF) puede ser aplicado a problemas que impliquen
deformaciones y esfuerzos, ya que permite resolver ecuaciones diferenciales parciales
para figuras arbitrarias, esto es debido a que las ecuaciones son muy difciles de resolver
mediante mtodos analticos, como la separacin de variables que solo abarca figuras
muy sencillas.

Independientemente de la naturaleza del problema, el MEF sigue los siguientes pasos:

1. Definicin del problema y su dominio
2. Discretizacion del dominio
3. Identificacin de las variables de estado
4. Formulacin del problema
5.

Establecimiento de los sistemas de referencia

6. Construccin de las formulaciones de aproximacin de los elementos

7. Determinacin de las ecuaciones al nivel de cada elemento
8. Transformacin de coordenadas
9. Ensamblaje de las ecuaciones de los elementos
10. Introduccin de las condiciones de contorno
11. Solucin del conjunto de ecuaciones simultneas resultantes
12. Interpretacin de los resultados

259

MODELADO DE PROCESOS DE MANUFACTURA

MEF aplicado al rango elstico

El MEF puede ser aplicado a muchos problemas que impliquen la solucin a ecuaciones
diferenciales por medio de mtodos numricos, en el caso del modelado del
comportamiento mecnico se tienen que resolver las siguientes ecuaciones: la
constitutiva, relacin entre deformacin y esfuerzo, las de equilibrio y las de velocidad, por
lo tanto se tienen 15 incgnitas, 6 debido a las componentes de las fuerzas y momentos,
6 esfuerzos (normales y tangenciales) y 3 componentes de velocidad. Aunque hay casos
donde se reducen las incgnitas de los problemas debido a que se pueden simplificar. El
comportamiento elstico lineal es utilizado en la industria porque muchos productos
trabajan bajo la consideracin de que no deben presentar una deformacin permanente o
un cambio de forma bajo la aplicacin de cargas externas, tambin es til cuando se
busca la solucin a problemas de esfuerzos y deformaciones a nivel microscpico, como
en las dislocaciones con ciertas condiciones particulares, en este parte slo analizaremos
el comportamiento elstico lineal.
La teora elstica lineal se basa en los siguientes cuatro puntos:

Las deformaciones son pequeas

El comportamiento del material es lineal

Los efectos dinmicos son despreciados

No se generan huecos o traslapes durante la deformacin

Como regla general se tiene que la deformacin debe de ser del orden de 10-2 de las
dimensiones del cuerpo, y el error asociado a este tipo de comportamiento no excedan el
1%. La decisin entre escoger un anlisis lineal depender de si los esfuerzos
sobrepasan el esfuerzo de cedencia o no, aunque si las deformaciones son pequeas el
resultado mediante anlisis lineal puede ser vlido, como en el caso del concreto. Por ello,
se presentan las ecuaciones que se resuelven por MEF en el rango elstico.

La idealizacin de que el comportamiento es esttico se debe a que la aceleracin cuando

se aplican las cargas es muy baja, por lo tanto no hay efectos inerciales. Un criterio para
determinar qu tan baja puede ser la aceleracin est determinado mediante:

f dlembert Ma
260

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

Si las f

dlembert

son pequeas en comparacin con las fuerzas inerciales, entonces es

aplicable.

La cuarta suposicin implica que el material es un medio continuo, sin traslapes o

cavidades en su interior.

Estas suposiciones tambin implican que el material no falla, evidentemente esto no es

cierto, en tales casos se deben considerar comportamientos no lineales.

Cinemtica
Si se define un vector con dos componentes una en la direccin i y la otra en la direccin
j, y utilizando coordenadas cartesianas, matricialmente se queda:

u x
u
u y
Se puede demostrar que en un volumen de control, bajo la suposicin de que existen
pequeos gradientes de desplazamiento, las deformaciones quedan definidas como:

xx

u y
u x
1 u y u x
, yy
, xy

x
y
2 x
y

As que las deformaciones se encuentran en trminos del campo de velocidades y se

puede escribir tensorialmente de la forma:

xx
u x

yy su
s ,
u y
xy

donde s 0

261

MODELADO DE PROCESOS DE MANUFACTURA

Esfuerzos y tracciones
Los esfuerzos se definen vectorialmente como

x xx i xy j , y xy i yy j , estos

son esfuerzos actuando en los planos x, y, respectivamente, en ocasiones es til

representar el tensor de esfuerzos de forma simtrica:

xx

xy

xy
yy

La traccin, al igual que los esfuerzos, es una fuerza aplicada sobre una superficie en
particular, mientras que los esfuerzos proveen informacin de la traccin en cualquier
punto.

De acuerdo a la imagen vector normal es n = nx i n y j , si se considera un espesor,

unitario el vector t queda como

t n

Ecuaciones de equilibrio
Si se considera que las fuerzas de traccin actan sobre el plano xy, y que puede que
existan fuerzas de cuerpo (debidas a campos gravitacionales, magnticos, etc)
definidas como:

b bx i by j
En un anlisis esttico la suma de fuerzas es igual a cero por lo que queda:

x y

b 0
x
y
262

CAPTULO 8. CONCEPTOS BSICOS DEL MTODO POR ELEMENTO FINITO

Escrito en forma matricial queda:

Ts b 0

Ecuaciones constitutivas
Las ecuaciones que relacionan los esfuerzos y las deformaciones se conocen como
ecuaciones constitutivas y en el rango elstico estn regidas por la Ley de Hook
generalizada.

D
Donde D es una matriz de 3 3 que depende en el caso de problemas bidimensionales si
se trata de un esfuerzo plano o de una deformacin plana.

Para esfuerzo plano:

Para deformacin plana:

1
E
D
1
1 2
0 0

E
(1 ) 1 2

1 / 2
0
0

1

0
0

/
2

0
0

FEM aplicado para deformacin plstica

Para el caso donde la deformacin se encuentra en el rango plstico los procesos son no
lineales, en el anlisis de conformado de metales como: forja, trefilado, laminacin en
caliente, etc., es normalmente aceptable despreciar las deformaciones elsticas frente a
las deformaciones plsticas o visco plsticas, el proceso es altamente considerado como
no-lineal en su comportamiento de deformacin y de la friccin de contacto, para la
mayora de los casos prcticos se admite un anlisis cuasi esttico, ya que los efectos
dinmicos pueden despreciarse a excepcin de los procesos de alta velocidad. En estos
casos, se analiza la evolucin de la deformacin de tal manera que en cada instante el
estado de un cuerpo en el curso de la deformacin satisfaga las condiciones de equilibrio.
263

MODELADO DE PROCESOS DE MANUFACTURA

La deformacin plstica generalmente sobrepasa a la deformacin elstica y las

idealizaciones de comportamiento rgido-plstico o rgido-viscoelstico son aceptables.

El problema asociado con la deformacin plstica es un problema de condiciones de

frontera. Para el caso de comportamiento rgido-visco plstico las condiciones de frontera
se establecen como sigue: en algn punto del proceso de distorsin cuasi esttica, la
forma de la figura, la distribucin interna de temperaturas y los valores de los parmetros
del material estn dados. As, el vector velocidad se encuentra definido en una superficie
Su a la cual se le aplica una carga F en cierta superficie SF. La solucin a este problema
es la distribucin de esfuerzos y velocidades que cumplan con las ecuaciones
constitutivas y las condiciones de frontera

Software libre
En los siguientes links podemos encontrar cdigos de uso libre para aplicar el MEF.
http://www.salome-platform.org/
http://www.dynaexamples.com/metalforming/metal-forming-iv

264