U.A.T.F.          Ing. Juan Carlos Barrios C.

  ©2011  Página: 1 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
Tema N° 1 
Introducción al Método de Elementos  Finitos 
1.1 Métodos matriciales 

En la actualidad el estudio de los problemas de la Ingeniería Civil es abordado por diferentes métodos 
y modelos matemáticos considerando el efecto de continuidad de los elementos estudiados tomando 
en cuenta las relaciones denominadas ecuaciones de gobierno. 

Por lo general, las estrategias de solución a estos problemas se muestran en el siguiente diagrama: 

 
Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil
 

En cada paso se producen simplificaciones por lo que puede decirse que la solución al problema físico 
en forma “exacta” no existe. 

En  los  problemas  enfocados  como  continuos  no  se  pierden  las  propiedades  físicas  en  dicho  medio 
continuo pudiendo valorar en cada punto del elemento de estudio, por lo es importante conocer los 
valores para resolver las ecuaciones integro – diferenciales. 

Existen varios métodos donde el procedimiento es sistemático, permitiéndonos discretizar y volver a 
un sistema de ecuaciones lineales algebraicas, cuyos resultados serán soluciones al problema. 

 
 U.A.T.F.          Ing. Juan Carlos Barrios C.  ©2011  Página: 2 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
1.2 Métodos matriciales 
e
La notación aplicada es:  Ij
    
e
R Vj
    j
  De la física tenemos:  e
  Ii
e
  Entonces:  ento
elem
Vi i
1
 

     
Siempre esta expresión puede escribirse matricialmente de la siguiente forma: 

 
Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil
1 1 1
 
1 1
O también: 
 
 
Ejemplo 2.‐  e
Fj
 
De la física tenemos:   σ E*ε 
e e j
A, E
Donde:      ó  
e
Fi L
Y:              ó     i

Por tanto:      

Despejando:       

Y para el caso general:     

Y                                           

Matricialmente se escribe: 

1 1
 
1 1

Por  tanto  se  puede  inferir  que  la  aplicación  de  estas  ecuaciones  pueden  ser  escritas    como  un 
sistema lineal de ecuaciones en forma matricial. 

 
 U.A.T.F.          Ing. Juan Carlos Barrios C.  ©2011  Página: 3 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
V2
Ejemplo 3.‐ 
2
Matricialmente el problema es dado por:  1
R
1
V1 I
 

R
I1

2
I
2
1 R3
Escribiendo para cada elemento tendremos: 
I3 I3
Elemento 1.‐  V
3 3
1 1
0

 
1 1

Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil

 
0
0
0 0 0
Elemento 2.‐ 
0 0 0
0 1 1
0
 
1 1
0

Elemento 3.‐ 
1 1
0
0 0 0 0  
1 1
0

Entonces  todos  los  Kij  se  obtienen  de  la  suma  de  las  matrices  individuales  y  a  este  proceso  se 
denomina ensamblaje de la matriz. 

Por tanto: 

   

Luego el sistema se coloca seleccionando los valores desconocidos y conocidos de la siguiente forma: 

Que al resolver se tiene:                              

 
 U.A.T.F.          Ing. Juan Carlos Barrios C.  ©2011  Página: 4 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
1.3 Análisis del problema como medio continuo.‐ 
En  el  método  de  elementos  finitos,  además  del  proceso  de  discretización  se  incorporan  2 
conceptos: El continuo ( Ω ) y su contorno ( Γ  ), puesto que las ecuaciones de gobierno deben 
ser considerados en todo su medio continuo y su comportamiento en su contorno, así en el 
proceso de discretización puede escribirse: 
 
En el elemento finito  e podemos escribir: 
                      
 
Donde      se  obtiene  al  imponer  al  elemento  las 
ecuaciones  integro  diferenciales  que  gobierna  el 

 
Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil
problema en su contorno y su continuo. 
continuo  
e
Es  importante  definir  el  dominio,  normalmente  se 
tiene un contorno fijo, aunque puede existir un análisis 
contorno con  mayores  criterios  donde  necesiten  varias 
q ecuaciones. 
 
 
1.4 Consideraciones de Derivación.‐ 

Si tenemos definida una función f(x, y), siempre podemos calcular:   ;   

Por tanto se define: 

Por tanto: 

 es como multiplicar un vector por un escalar 

Si: 
 
Entonces el producto escalar será: 
°  
Y el producto vectorial: 
           
Ahora bien, 

 
 U.A.T.F.          Ing. Juan Carlos Barrios C.  ©2011  Página: 5 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
Si se desea determinar el área de continuo (Ω) a partir de una función definida tal como se muestra 
en la figura siguiente, se verifica que: 
z

f(x,y) °° Ω ° Γ 
Ω
y
Desarrollando tenemos: 

° Γ 

 
 

Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil

Para determinar el área, se tiene lo siguiente: 

x
y  
2
= y2 (x)

, ,  
1 = y1 (x)
x
a b  

Por  tanto,  la  evaluación  del  continuo  depende  de  la  definición  de  las  funciones  de  contorno  que 
encierra dicha superficie. 

Ahora bien, si se tiene un dominio tal como se muestra en la figura siguiente, se establece: 

Por Pitágoras, se aproxima: 
y

Γ  
d Si queremos conocer: 
dy
dx

x , Γ 

Donde:    x = x(ξ)   y  y=y(ξ), en la que se verifica: 

                  

 
 U.A.T.F.          Ing. Juan Carlos Barrios C.  ©2011  Página: 6 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
Entonces: 

Γ  

Reemplazando en la integral: 

, Γ  

En el caso general, se tiene:                , por tanto 

 
        

Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil

 
Por tanto: 

Γ 1  

Entonces: 

, Γ 1  

Haciendo un zoom a la fracción de curva del contorno, observamos: 

Siendo: 
t    
d    
dy Entonces: 

                                      Por tanto:            
dx
Entonces: 

° Γ Ω 

Por tanto generalizando se dice: 

°° Ω ° Γ 
Ω

 
 
 
 
 U.A.T.F.          Ing. Juan Carlos Barrios C.  ©2011  Página: 7 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
1.5 Consideraciones de Integración Numérica 

En el método de elementos finitos normalmente se integra entre el rango ‐1, +1, es decir:  

f( )
 

Por ejemplo para una función de primer grado, 

Tenemos: 

 
-1 0 +1
2

Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil

 

Si f(ξ) es una función polinómica de grado n 

1 1 1
 
2 3 1 |

1 1
1 1 1 1  
2 1

Por tanto: 

Donde wi se denomina peso de la función. 

Por tanto para una función lineal:  0    2 

Asimilando para un polinomio de grado n, tenemos un sistema de ecuaciones como el siguiente: 

     
 

     

En el sistema no se conocen los pesos ni los puntos de observación, entonces: 

2 

 
 U.A.T.F.          Ing. Juan Carlos Barrios C.  ©2011  Página: 8 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
0 

2/3 

… 
1 1
 
1
Por tanto las incógnitas del sistema son: 
           …  
 
           …   
 
Para que el sistema tenga solución, se debe cumplir la siguiente relación: 

 
1 2  

Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil

Donde:   p= grado de la función polinómica 
    n = N° de puntos de observación (cuadratura de Gauss) 
Por tanto, para integrar n polinomio de grado n, debe evaluarse dicha función teniendo en cuenta su 
grado, entonces determinar los puntos de observación. Para ello se dispone de 2 métodos numéricos 
que son: 
Gauss Legendre: 
Donde necesitamos los n puntos expresados a partir de la relación: 
1
 
2
Gauss Lobato: 
    Que necesita n puntos a partir de la expresión: 
3
 
2
Debe  anotarse  que  para  una  función  cualquiera  entre  cualquier  parábola  cúbica  que  pasa  por  2 
puntos, se obtiene la misma sección. 
Asimismo,  según  Gauss  Lobato  (que  se  aplica  en  problemas  de  fluidos)  se  definen  2  puntos  y  los 
puntos medios son recomendados en función de ver su contorno. 
 
1.6 Tensor de esfuerzos  
 
El  esfuerzo  depende  de  la  condición  de  vector  y  la 
superficie donde actúa, por lo que el esfuerzo es un 
tensor de rango, cuya definición es: 
 
   P 
Esfuerzo  lim A0  
   A 
 
Este  concepto  se  aplica  tomando  en  cuenta  que  el 
material es continuo. 
 

 
 U.A.T.F.          Ing. Juan Carlos Barrios C.  ©2011  Página: 9 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
La nomenclatura esta dado por: 
 
  xy  xz    x  xy  xz   X x X xy X xz 
  x     
  yx  y  yz    yx  y  yz   Yyx Y y Yyz 
  zx  zy  z   zx  zy  z   Z zx Z zy Z z 
 
Donde: 
 
 i (i  x, y , z ) Esfuerzos normales
   ij (i, j  x, y, z) i  j Esfuerzos tan genciales
 
 
 

 
Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil
Todo cuerpo alcanza su estado de equilibrio, cumpliendo con las siguientes ecuaciones: 
 

0 

0 

0 

Donde:     Se denomina Vector de fuerza de masa (o gravitacional) 
 
El  esfuerzo  normal  en  un  plano,  dado  su 
estado  tensional  ortogonal  esta  dado  por  la 
siguiente relación: 
 
 
 
 
Donde: 
L, m, n se denominan cosenos directores 
 
1.7 Tensor de deformaciones.‐ 
El estado de deformación de un cuerpo infinitesimal está definido por el tensor: 
 
 x  xy  xz 
 
 y  yz    x ,  y ,  z ,  xy ,  xz ,  yz T
 
    yx
  zy  z 
   zx
 

 
 U.A.T.F.          Ing. Juan Carlos Barrios C.  ©2011  Página: 10 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
Se  establece  la  relación  entre  deformación  y  desplazamiento  a  partir  de  las  siguientes 
relaciones: 
u u v
  x   xy  
  x y x
v u w
  y   xz  
  y z x
  w v w
z   yz  
  x z y
 
Siendo u, v, w el campo de los desplazamientos en las direcciones x, y, z respectivamente. 
Para determinar la deformación en un punto cualquiera del plano, se establece a partir de las 

 
Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil
relaciones: 
 
2  
 
Tomando  en  cuenta  las  relaciones  de  deformación  de  un  nodo,  se  establece  que  para  los 
elementos  concurrentes  en  dicho  nodo  se  produce  una  condición  de  compatibilidad,  por 
tanto dichas ecuaciones se escriben como:  
 

2  

2  

2  

 
 
 

 
 U.A.T.F.          Ing. Juan Carlos Barrios C.  ©2011  Página: 11 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
1.8 La Ley de Hoocke como ecuación que gobierna los problemas elásticos estructurales‐ 
Se  conoce  como  Ley  Generalizada  de  Hoocke  a  la  expresión  que  relaciona  el  tensor  de  los 
esfuerzos  con  el  tensor  de  las  deformaciones.  Dependiendo  de  la  relación,  ésta  se  puede 
escribir de 2 maneras: 
 
a) Condición  de  flexibilidad.‐  que  relaciona  el  tensor  de  las  deformaciones  en  función  del 
tensor de esfuerzos, es decir:    C
 
    x  c1,1 c1, 2 c1,3 c1, 4 c1,5 c1,6   x 
     c2, 2 c2,3 c2, 4 c2,5 c2,6   y 
 y   
    z   c3,3 c3, 4 c3,5 c3,6   z 

 
Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil
     
  xy   c4, 4 c4,5 c4,6   xy 
 
 yz   c5,5 c5, 6   yz 
      
   zx   c1,6   zx 
Donde ci,j, se denomina: coeficiente de flexibilidad 
 
b) Condición de rigidez.‐ que relaciona el tensor de los esfuerzos en función del tensor de las 
deformaciones, es decir: 
  D
 
 x  d1,1 d1, 2 d1,3 d1, 4 d1,5 d1,6    x 
     d 2, 2 d 2,3 d 2, 4 d 2,5 d 2,6    y 
   y   
   
 z   d d d d 
3, 6    z 
 
3, 3 3, 4 3, 5
 
   xy   d 4, 4 d 4,5 d 4,6   xy 
   yz   d 5,5 d 5, 6   yz 
    
   zx   d1,6   zx 
 
Donde di,j, se denomina coeficiente de rigidez 
 
Para  el  caso  de  un  material  isotrópico  y  simétrico,  esta  relación  se  reduce  a  los  siguientes 
términos: 
    x  c1,1 c1, 2 c1,3   x 
    
   y 
c2 , 2 c2 , 3   y 
    z   c3,3   z 
   
 xy     xy 
  c4 , 4
   yz   c5,5   yz 
      
 zx   c1, 6   zx 
 
Estos  coeficientes  pueden,  a  su  vez,  ser  representado  en  función  de  2  constantes, 
denominadas:  E  =  Módulo  de  Young  y  v=  Constante  de  Poissón,  por  tanto  las  relaciones 
establecen: 

 
 U.A.T.F.          Ing. Juan Carlos Barrios C.  ©2011  Página: 12 
Guía Didáctica  CIV 313 Elementos Finitos 
 
 
 
x y z  xy
  x     xy 
  E E E G
x y z  yz
   x      yz 
  E E E G
  x y z  zx
 x      zx 
  E E E G
  E
G
  2(1  v)
 
Matricialmente, se escribe: 

 
Introducción del Método de Elementos Finitos aplicados a la Ingeniería Civil
 
x  1     x 
       
   y  1   y 
   
 z 1  1   z 
    
   xy  E  21      xy 
   yz   21      yz 
    
   zx   21     zx 
 
 
1.9 Otras relaciones constitutivas 
El  método  de  Elementos  Finitos  no  solamente  se  aplica  a  problemas  estructurales,  sino  a 
diversos casos de estudios de la ingeniería civil. Por ello, se presenta una tabla en la que se 
pueden identificar algunos criterios referidos a los diferentes problemas físicos que aborda el 
ingeniero civil y que conducen a los criterios constitutivos que se aplican: 
 
Problema Físico  Principio de  Variable de  Flujo  Constante del  Ecuación 
Conservación  Estado  Material  Constitutiva 
Deformación de  Equilibrio de  Desplazamiento  Esfuerzo o  Módulo de  Ley de Hoocke 
un cuerpo  Fuerzas  o Fuerzas  deformación  Young y Ctte. de 
elástico  Poisson 
Red eléctrica  Equilibrio de  Voltaje o  Flujo eléctrico  Conductividad  Ley de Kirchoff 
corriente  Amperaje  eléctrica 
Torsión  Energía  Función de  Razón de giro  Esfuerzo de corte  Ley de Hoocke 
Potencial  esfuerzos 
Transferencia de  Energía  Temperatura  Flujo de calor  Conductividad  Ley de Fourier 
calor  térmica 
Flujo de fluídos  Momento  Gasto hidráulico  Esfuerzo  Viscosidad  Ley de Stokes 
cortante 
Flujo a través de  Masa  Velocidad,  Razón de flujo  Permeabilidad  Ley de Darcy 
un medio poroso   gasto hidráulico 
Electrostática  Flujo eléctrico  Potencial  Flujo eléctrico  Permitividad  Ley de 
eléctrico  Coulomb 
Magnetostática  Potencial  Potencial  Flujo magnético  Permeabilidad  Ley de Maxwell 
magnético  magnético  magnética