EL MÉTODO DE LOS

ELEMENTOS FINITOS

Por: [Link]. Ing. Edwin Félix Achá Colque

 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Contenido
 Conceptos de elasticidad  Tipos de Estructuras
 Conceptos energéticos 1)Estructuras unidimensionales
 Conceptos matemáticos a) Barras y estructuras articuladas
 Descripción del MEF b) Vigas y estructuras reticuladas
1) División de elementos finitos 2) Estructuras bidimensionales
2) Vector de desplazamiento a) Elasticidad bidimensional
3) Matriz de rigidez b) Elemento triangular de tres nodos
4) Matriz completa de rigidez de la estructura c) Elemento rectangular de cuatro nodos
5) Respuesta de la estructura 3) Placas delgadas
 Análisis del MEF 4) Estructuras tridimensionales
1) Condiciones de las funciones de desplazamiento a) Elemento tetraédrico de cuatro nodos
2) Criterios de convergencia b) Elemento hexaédrico de ocho nodos
3) Equilibrio de la estructura
4) Estabilidad de los elementos

 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Contenido x1
P2 P3 F2 y
y F2
P1 4 5 6 2
1 F2 x
y1

F1 y
1
F1
1 x
1 2 3 F1x
y
y

q(x)
1

x
2
y2
x x M0 y1
x
1 2
q x1
x2
2
y
1
y1

P1 v2 y1
P2
2 x1
v1y1 u 2 x1
u1x1
4
2
1 x

2
1
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Contenido
w4
z 4
v4
u4

w3
w1 3 v3
1
v1
w2 u3
u1
2 v2
y
u2

z w8 w7
8 7
v8 w6 v7
w5
u8 u7
5 6
v5 w4 v6 w3
u5 4 u6
v3
v4 3
w1 w2
u4 u3
1
2 v2
v1
u1 u2
y

x
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Contenido
EL MÉTODO DE LOS
ELEMENTOS FINITOS

BASES TEORICAS
Por: [Link]. Ing. Edwin Félix Achá Colque
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Introducción
El método de los elementos finitos MEF es una herramienta poderosa para la solución numérica de un amplio
rango de problemas:
• Análisis Estructural (Edificios, aeronaves, puentes, etc.)
• Campos del flujo de calor
• Fluidos
• Magnetismo
• Filtraciones
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Desplazamientos 𝑢
Deformaciones 𝜀
Esfuerzos 𝜎

La teoría de la elasticidad determina el comportamiento de una estructura ante las cargas mediante la
resolución de sistemas de ecuaciones diferenciales llamadas ecuaciones de campo o de gobierno
complementados por las condiciones de contorno o de sustentación.
Los métodos tradicionales de cálculo emplean la resistencia de materiales que añade hipótesis simplificatórias
a la teoría de la elasticidad para resolver las ecuaciones de gobierno planteadas, con lo cual se obtiene soluciones
que se consideran exactas.

Ejemplo:

“Las caras permanecen planas antes y después de la deformación” (Navier-Bernouilli)

 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Estructuras Discretas y Continuas

Los métodos tradicionales de cálculo normalmente solo se pueden aplicar a estructuras articuladas y
reticuladas (Discretas). En estructuras más complejas como estructuras continuas es necesario aplicar métodos
aproximados como el Método de Elementos Finitos MEF (FEM).
Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

P P

a)
P P

b)
Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

a) b) c)
Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Método de los Elementos Finitos

Es dividir la estructura en partes, llamados ELEMENTOS FINITOS, que están interconectados entre si, el método
reduce el problema elástico a la resolución de un sistema de ecuaciones algebraicas,

Cuchillo Tipo carga

Elemento Finito Tipo Losa k

Carga Puntual Vertical dentro del
Elemento
l
j a11 x1  a12 x 2  a13 x3    a1n x n  b1
i a 21 x1  a 22 x 2  a 23 x3    a 2 n x n  b2
a31 x1  a32 x 2  a33 x3    a3n x n  b3
     
Carga Distribuida Verticalmente a n1 x1  a n 2 x 2  a n 3 x3    a nn x n  bn

Carga Distribuida dentro del Elemento

Elemento de Viga
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

El MEF se basa en una serie de conceptos que son:

- Conceptos de Elasticidad.
- Conceptos Energéticos.
- Conceptos Matemáticos.


Esfuerzo real de Fractura

U
Esfuerzo Ultimo B
f Esfuerzo de Fractura
u3
A
dFt 
h3 du3e3
Y
E
P

Ft 
P h2 du2e2

 u2
u1 h1du1e1

Región Elastica Fluencia Endurecimiento por deformación Estricción

 0 D C
dV

t
Comportamiento Elástico Comportamiento Plástico d t
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Elasticidad
Vector de Esfuerzos
Idealización de una fuerza
concentrada

Fuerza supeficial  F  0;  F
x y  0;  Fz  0
 M  0;  M
x y  0;  M z  0
Idealización de una carga
linealmente distribuida

n
Fn
F F Fn
  lim A0
A A
Ft t
Ft
  lim A0
A
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Elasticidad
Vector de Esfuerzos
 yz
 xz
z z
y
 xy
Fz F z  yx
x
Fy  zy
y y
Fx  zx
z
x
x
z
 zy  x 
 zx  yz  
 y
 xz y  z 
 yx y     
z  xy  xy 
x x  yz 
 
 zx 
y
x
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Elasticidad
Ecuaciones de Equilibrio
z
Fx 0
z  Fy  0
z
F
z 0
z
 zy
 xz   xz
 x   x  xy   xy
 zx qz
 y   y qy y  yz   yz  yz
 yx   yx qx  xz y
 yx z  xy
y z
y
z x
x x
 zx   zx x  zy   zy
y y  z   z
y x
x
x

 x  yx  zx  xy  y  zy  xz  yz  z
   qx  0    qy  0    qz  0
x y z x y z x y z

Estas ultimas 3 ecuaciones representan las ecuaciones de equilibrio del paralelepípedo elemental.
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Elasticidad
Vector de Deformación
Deformación unitaria normal

s's
 prom 
s
x 
  
 nt   '  y
2 B Lim
A a lo largo n
z 
    
C  A a lo largo t
Deformación unitaria cortante

 xy 
 yz 
 
 zx 
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Elasticidad
Ley de Hooke generalizada
x

z x 
1
E

 x  v y   z  
x

y
1
E

 y   y  v x   z  
y
1
E

 z   z  v x   y  
y
z

x En una dimensión

x 
1
 x 
E

z
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Elasticidad
Ley de Hooke generalizada

 yz  zx


 xy

L
1
 xy   xy
G
P
1
 yz   yz
G
1
 zx   zx
G
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Elasticidad
Relación Esfuerzo - Deformación (Ecuaciones de Lame )
 1 v v v 0 0 0 
x 
1
 
 x  v y   z   1

v

v
 x 
 
 v
 1 v v 0 0 0   x 
 
E 
 x   E E E   x   y  v v 1 v 0 0 0   y 
v  

 y   y  v x   z 
1
     v 1
   y 
 z 
 
E 
 0 0 0
1  2v
0
 
0   z 
 xy  1  v 1  2v  
E  y  E E E 2   xy 
 z   v 1   z  1  2v
v  yz   0 0   yz 

 z   z  v x   y 
1   0 0 0
 E E E     2  
E  zx   1  2v   zx 
 0 0 0 0 0
2 
1
 xy   xy 1
0 0
    D 
 xy   G
G

1     0 1   xy  Matriz Constitutiva
 yz   yz
 yz   0   yz 
G
 zx  
G    1 v v v 0 0 0 

1 1   zx   v
 1 v v 0 0 0 

0 0
 zx   zx  G   v

v 1 v 0 0 0 

G 1  2v
D  E
 0 0 0 0 0 
1  v 1  2v   2
1  2v 
 0 0 0 0 0 
 2 
 1  2v 
 0 0 0 0 0
2 
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Elasticidad
u v w
Relación Deformación - Desplazamiento x  ; y  ; z 
x y z
dy  dv
u v v w w u
B' ' '  xy   ;  yz   ;  zx  
z
v y x z y x z
1
 B' dv B' ' w  dw
v A'

w w
     
w  dw  x 0 0 0 0
B  x
    
A dy v  dv 
 x   0 0 0 0
y
   y   y 
 y   u   
z   0 0
z   v  0 0
z 
    
x
      
v  dv  xy   0   w  0
dv  yz   y x    y x 
C'' C'''        
 zx   0  0
  dw
 z y  z y 
z w  dw  
C'
dz  dw  0
  0

2  z
w  dw  dz C x   z x 
u  A'
    u
dz
u    v    w
v w
A v
 w y
  el operador matricial
x
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Elasticidad
Problema Elástico
El problema elástico que consiste en determinar en cada uno de los puntos de la estructura las 15 magnitudes siguientes:
- 6 componentes del estado de esfuerzos  x ,  y ,  z , xy , yz , zx
- 6 componentes del estado de deformaciones  x ,  y ,  z ,  xy ,  yz ,  zx
- 3 componentes del desplazamiento u , v, w

Estas 15 magnitudes están relacionadas entre si mediante 15 ecuaciones que son el fundamento de la teoría de elasticidad:

- 3 ecuaciones de equilibrio - 6 ecuaciones que relacionan - 6 ecuaciones de compatibilidad

esfuerzos y deformaciones de las deformaciones

 x  yx  zx  x    2u  x   y   z
u
   qx  0 x 
x
x y z
 y   x    2u  y   z y 
v
y
 xy  y  zy
   qy  0  y   x   y    2u  z z 
w
x y z z
 xy  G xy  xy 
u v

 xz  yz  z y x
x

y

z
 qz  0  yz  G yz v w
 yz  
 zx  G zx z y
w u
 zx  
x z
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Elasticidad
Problema Elástico
Al resolver este sistema de 15 ecuaciones en derivadas parciales hay que tener en cuenta las condiciones de contorno que
comprenden dos tipos de prescripciones:
a) Desplazamientos conocidos debido a la actuación de enlaces.
b) Esfuerzos conocidos originados en las proximidades de las fuerzas aplicadas.

A veces interesa expresar las ecuaciones de equilibrio en función de los desplazamientos y a las ecuaciones diferenciales
resultantes se las denomina ecuaciones diferenciales de gobierno.
qx  f x 
PA PB

A x 
d  du 
A B

 EA   q x 
x dx

N
qx 
N  dN
dx  dx 

dx

La solución del problema elástico es única y solamente puede hallarse de forma exacta en casos particulares mediante la Teoría
de la Elasticidad. Por el contrario, la Resistencia de Materiales resuelve el problema elástico de forma más simple utilizando,
usualmente, métodos aproximados a partir de hipótesis simplificadoras.
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Energéticos
Trabajo de un vector
z

P1 F
P2  


dr w F  dr
r (t ) P1
P2
 
r (t )  dr
 
y w   M  d

Trabajo Externo e Interno

La mayoría de los métodos de energía se basan en el principio de la conservación de la energía, que establece que el
trabajo externo hecho por todas las fuerzas externas que actúen dentro de la estructura, se transforma en trabajo interno
de energía de deformación, que se desarrolla cuando la estructura se deforma.
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Energéticos
Trabajo sobre cuerpo deformable

B
Ft  dFt
A
dFt
L

dL Ft
d t

0 D C
 t
Ft  dFt d t

F
U  w   Ft d t
0
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Energéticos
Trabajo sobre cuerpo rígido

P2
P2

P1
P1
Pi
Pi
n
w   wi  0  0
i

i 1

P3
P3

“En un sistema de fuerzas en equilibrio actuando sobre un cuerpo rígido, es nulo el trabajo de este sistema cuando el
cuerpo sufre un desplazamiento”.
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Energéticos
Principio del trabajo virtual

z Fi F F
F
F N
V N M
y M M
M
T
F i N
T T
N
V M
T
V
V M V
T
y V T
dy dy
dy N
dy N

dWe  dWr  dWd dWr  0 dWe  dWd

Esta expresión dice que el trabajo virtual total efectuado por todas las fuerzas que actúan sobre el elemento durante su
desplazamiento virtual (esto estando el trabajo virtual debido al desplazamiento de cuerpo rígido como el trabajo virtual
debido a la deformación del elemento) es igual al trabajo virtual realizado por esas mismas fuerzas solamente durante la
deformación virtual del elemento.
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Energéticos
Principio del trabajo virtual

 dW   dW
e d

Wext  Wint
Esta ecuación representa el principio del trabajo virtual y se puede expresar como:

”Si a una estructura deformable que esta en equilibrio bajo la acción de un sistema de cargas se le da una pequeña
deformación virtual, entonces el trabajo virtual efectuado por las fuerzas externas (o cargas) es igual al trabajo virtual
realizado por las fuerzas internas (o resultantes de esfuerzos).”
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Energéticos
Principio del trabajo virtual
Trabajo virtual de las cargas externas
  n
Wext   F  dr   F cos dr   Fi i
i

Trabajo virtual de las cargas internas

Trabajo virtual de la fuerza normal Trabajo virtual del momento flector Trabajo virtual de la fuerza cortante Trabajo virtual del momento torsor
d i d d

d i
T di
dd
V

M M  dM

d d
V  dV
dy
dy T  dT
d

dWint  Vd
dWint  Ndy dWint  Md dWint  Td
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Coordenadas naturales
Las coordenadas adimensionales que definen la posición de un punto se denominan coordenadas naturales. Según que estas
coordenadas se determinen mediante relaciones de longitudes (coordenadas de longitud), de áreas (coordenadas de área)
o de volúmenes (coordenadas de volumen), es decir son adimensionales.

a) Coordenadas de Longitud

L
L1 x2  x 1 1
1 P 2 1  1  1
1 1
L L 1 1 x 2  1 1 1 x  1 1 x  x1 x 2  x1 1 L 1
 L 2 L 2 L L
x1 L2 L1 L x  x1  2 L  x1 1 x x  x1 1 2 x1 x2  2 x1 x2  2
2  2 2 
L L x2  x1 x 2  x1 L L
x
x2
1 1 1 1

x x1 x2  2
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Coordenadas naturales
b) Coordenadas de Área
Sea un triangulo con vértices en los puntos 1, 2, 3, y un punto P cualquiera en el interior del triangulo definido por sus
coordenadas cartesianas x e y. El punto divide el triangulo en tres superficies de áreas A1,A2 ,A3 . La posición del punto P se
determina en un sistema de coordenadas naturales llamadas coordenadas de área, mediante las siguientes relaciones.
y 3
y 3 y1 x1
y x2
2A  2  y1 x 2  y 2 x3  y 3 x1  y 2 x1  y 3 x 2  y1 x3
y3 x3
A2 A2 y1 x1
A1 A1
y3 y x
P P y x2
2 2
2 A1  2  yx2  y 2 x3  y 3 x  y 2 x  y 3 x 2  yx3
A3 y A3 y3 x3
y2 y x
1
1 y1
y1 x1
x1 x y x
x 2 A2   y1 x  yx3  y 3 x1  yx1  y 3 x  y1 x3
x3 y3 x3
A1 A A y1 x1
1  ; 2  2 ; 3  3 x
A A A y1 x1
x2
y x2
1   2   3  1 2 A3  2  y1 x 2  y 2 x  yx1  y 2 x1  yx2  y1 x
y x
y1 x1
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Coordenadas naturales
b) Coordenadas de Área
A1 yx2  y 2 x3  y3 x  y 2 x  y3 x2  yx3  y  y2 x  x2  x3 y   y 2 x3  y3 x2 
1    3
A y1 x2  y 2 x3  y3 x1  y 2 x1  y3 x2  y1 x3 A
1 y3  y 2 x 2  x3 y 2 x3  y 3 x 2 x
1
2 
A2

y1 x  yx3  y3 x1  yx1  y3 x  y1 x3  y  y3 x  x3  x1 y   y3 x1  y1 x3 
 1  2  y1  y3 x3  x1 y3 x1  y1 x3 y
A y1 x2  y 2 x3  y3 x1  y 2 x1  y3 x2  y1 x3 A A
3 y 2  y1 x1  x2 y1 x2  y 2 x1 1
A3 y1 x2  y 2 x  yx1  y 2 x1  yx2  y1 x  y  y1 x  x1  x2 y   y1 x2  y2 x1 
3    2
A y1 x2  y 2 x3  y3 x1  y 2 x1  y3 x2  y1 x3 A

1
x y3  y 2 x 2  x3 y 2 x3  y 3 x 2 1 x1 x2 x3 1 x x1 x2 x3 1 x x1 x2 x3 1
1 1
y  A y1  y 3 x3  x1 y 3 x1  y1 x3 2  A y1 y2 y3  2 y  A y1 y2 y3  2 1
y  A y1 y2 y3  2
y1 x 2  y 2 x3  y 3 x1  y 2 x1  y 3 x 2  y1 x3 A
1 3
A
1 y 2  y1 x1  x 2 y1 x 2  y 2 x1 3 1 1 1 3 1 1 1 1 1 1 1 3

x x1 x2 x3 1
y  y1 y2 y3  2
1 1 1 1 3
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Coordenadas naturales
c) Coordenadas de Volumen
z 4
4 4
V1 V V V
1  ; 2  2 ; 3  3 ; 3  4
4 A A A A
P
P P P

1   2   3   4  1
3 3 3
1 1 1 P

2 2 2 1 2

y Volumen 1 Volumen 2 Volumen 3 Volumen 4

x
Sean los vectores , que mediante un sexto del triple producto escalar definen el volumen del los tetraedros dado tres vectores.
1      
V  12 13 14 
6   

12  x 2  x1 , y 2  y1 , z 2  z1  V1 
1      
12 1P 14 
 6   
13  x3  x1 , y 3  y1 , z 3  z1  1      
V2  12 13 1P 
 6   
14  x 4  x1 , y 4  y1 , z 4  z1  1      
V3  1P 13 14 
6 

1P  x  x1 , y  y1 , z  z1   
1              
V4  12 1P    13 1P   14 1P  
6       
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Coordenadas naturales
x 2  x1 y 2  y1 z 2  z1
c) Coordenadas de Volumen x  x1 y  y1 z  z1
x 2  x1 y 2  y1 z 2  z1
V x 4  x1 y 4  y1 z 4  z1 1
Tetraedro 1: 1  1   x  x1 y  y1 z  z1
V x 2  x1 y 2  y1 z 2  z1 V
x 4  x1 y 4  y1 z 4  z1
z 4
x3  x1 y 3  y1 z 3  z1
x 4  x1 y 4  y1 z 4  z1

V1  x  x1  y 4  y1 z 2  z1    y 2  y1 z 4  z1    y  y1 x 2  x1 z 4  x1   x 4  x1 z 2  z1 

P  z  z1 x 4  x1  y 2  y1   x 2  x1  y 4  y1 
3
1 V1  x y 4  y1 z 2  z1    y 2  y1 z 4  z1   yx 2  x1 z 4  x1   x 4  x1 z 2  z1  
 zx 4  x1  y 2  y1   x2  x1  y 4  y1   x1  y 4  y1 z 2  z1    y 2  y1 z 4  z1  
 y1 x2  x1 z 4  x1   x 4  x1 z 2  z1   z1 x 4  x1  y 2  y1   x2  x1  y 4  y1 
2
y

Ax  By  Cz  D
1 
V
x donde:
A   y 2  y1 z 3  z1    y 3  y1 z 2  z1 
B  x3  x1 z 2  z1   x 2  x1 z 3  z1 
C  x 2  x1  y 3  y1   x3  x1  y 2  y1 
D   x1  y 2  y1 z 3  z1    y 3  y1 z 2  z1   y1 x3  x1 z 2  z1   x 2  x1 z 3  z1  
 z1 x 2  x1  y 3  y1   x3  x1  y 2  y1 
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Coordenadas naturales
x 2  x1 y 2  y1 z 2  z1
c) Coordenadas de Volumen
x3  x1 y 3  y1 z 3  z1
x 2  x1 y 2  y1 z 2  z1
Tetraedro 2: V x  x1 y  y1 z  z1 1
2  2   x3  x1 y 3  y1 z 3  z1
V x 2  x1 y 2  y1 z 2  z1 V
z x  x1 y  y1 z  z1
x3  x1 y 3  y1 z 3  z1
4

x 4  x1 y 4  y1 z 4  z1

V 2  x  x1  y 2  y1 z 3  z1    y3  y1 z 2  z1    y  y1 x3  x1 z 2  z1   x 2  x1 z 3  z1 

 z  z1 x 2  x1  y3  y1   x3  x1  y 2  y1 
P

3
1

V 2  x y 2  y1 z 3  z1    y3  y1 z 2  z1   yx3  x1 z 2  z1   x2  x1 z 3  z1  

 zx 2  x1  y3  y1   x3  x1  y 2  y1   x1  y 2  y1 z 3  z1    y3  y1 z 2  z1  
2
y
 y1 x3  x1 z 2  z1   x 2  x1 z 3  z1   z1 x 2  x1  y3  y1   x3  x1  y 2  y1 

Ex  Fy  Gz  H
2 
V
x donde:

E   y 2  y1 z 3  z1    y 3  y1  z 2  z1 
F  x3  x1 z 2  z1   x 2  x1 z 3  z1 
G   x 2  x1  y 3  y1   x3  x1  y 2  y1 
H   x1  y 2  y1 z 3  z1    y 3  y1 z 2  z1   y1 x3  x1 z 2  z1   x 2  x1 z 3  z1  
 z1 x 2  x1  y 3  y1   x3  x1  y 2  y1 
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Coordenadas naturales
x  x1 y  y1 z  z1
c) Coordenadas de Volumen x3  x1 y 3  y1 z 3  z1
x  x1 y  y1 z  z1
V x 4  x1 y 4  y1 z 4  z1 1
Tetraedro 3: 3  3   x3  x1 y 3  y1 z 3  z1
V x 2  x1 y 2  y1 z 2  z1 V
x 4  x1 y 4  y1 z 4  z1
z 4
x3  x1 y 3  y1 z 3  z1
x 4  x1 y 4  y1 z 4  z1

V 3  x  x1  y3  y1 z 4  z1    y 4  y1 z 3  z1    y  y1 x 4  x1 z 3  z1   x3  x1 z 4  z1 

P
 z  z1 x3  x1  y 4  y1   x 4  x1  y3  y1 
3

V 3  x y3  y1 z 4  z1    y 4  y1 z3  z1   yx4  x1 z3  z1   x3  x1 z 4  z1  

1

2  zx3  x1  y4  y1   x4  x1  y3  y1   x1  y3  y1 z 4  z1    y4  y1 z3  z1  

y  y1 x4  x1 z3  z1   x3  x1 z 4  z1   z1 x3  x1  y4  y1   x4  x1  y3  y1 

Ix  Jy  Kz  L
3 
V
x donde:

I   y 3  y1 z 4  z1    y 4  y1 z 3  z1 
J  x 4  x1 z 3  z1   x3  x1 z 4  z1 
K  x3  x1  y 4  y1   x 4  x1  y 3  y1 
L   x1  y 3  y1 z 4  z1    y 4  y1 z 3  z1   y1 x 4  x1 z 3  z1   x3  x1 z 4  z1  
 z1 x3  x1  y 4  y1   x 4  x1  y 3  y1 
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Coordenadas naturales
x2  x y2  y z2  z
c) Coordenadas de Volumen x3  x y3  y z3  z
x2  x y2  y z2  z
V x4  x y4  y z4  z 1
Tetraedro 4: 4  4   x3  x y3  y z3  z
V x 2  x1 y 2  y1 z 2  z1 V
x4  x y4  y z4  z
z 4
x3  x1 y 3  y1 z 3  z1
x 4  x1 y 4  y1 z 4  z1

P V 4  x y 2 z 4  z 3   y3 z 2  z 4   y 4 z 3  z 2   yx 2 z 3  z 4   x3 z 4  z 2   x 4 z 2  z 3  
3  zx2  y 4  y3   x3  y 2  y 4   x4  y3  y 2   x2  y3 z 4  y 4 z 3   x3  y 4 z 2  y 2 z 4  
 x4  y 2 z 3  y3 z 2 
1

2
y
Mx  Ny  Oz  P
4 
V

x donde:
M  y 2 z 4  z 3   y3 z 2  z 4   y 4 z 3  z 2 
N  x 2  z 3  z 4   x3  z 4  z 2   x 4  z 2  z 3 
O  x 2  y 4  y 3   x3  y 2  y 4   x 4  y 3  y 2 
P  x 2  y 3 z 4  y 4 z 3   x3  y 4 z 2  y 2 z 4   x 4  y 2 z 3  y 3 z 2 
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Coordenadas naturales
c) Coordenadas de Volumen
Expresando de forma matricial:
1 A B C D x
2 1 E F G H y

z 4
3 V I J K L z
4 M N O P 1

1
x A B C D 1
y 1 E F G H 2

P z V I J K L 3
3 1 M N O P 4
1

2
x x1 x2 x3 x4 1
y
y y1 y2 y3 y4 2

z z1 z2 z3 z4 3
x 1 1 1 1 1 4
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Jacobiano de transformación
Considerando las coordenadas rectangulares (x,y,z) de un punto P en función de las nuevas variables (u1,u2,u3)que
llamaremos coordenadas curvilíneas de P , en la forma:

x  x(u1 , u2 , u3 ); y  y(u1 , u2 , u3 ); z  z(u1 , u2 , u3 )

u1  u1 ( x, y, z); u2  u2 ( x, y, z); u3  u3 ( x, y, z)

Coordenadas curvilíneas ortogonales 

r Denominando como factores de
 u1 escala y siendo iguales a:
e1  
z
Linea u 3 z r
u1
u1
P  x, y , z 

dr   
u1  c1

r r r
u 2  c2 r h1  ; h2  ; h3 
u1 u 2 u 3
 u2
 
r u1   d r

e2 
r u1 
P 
u 3  c3 Linea u 2 r
Linea u1 u2
y

  
1 r 1 r 1 r
y
r   
 u3 e1  ;e2  ;e3 
x e3   h1 u1 h2 u 2 h3 u 3
r
u3
x
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Jacobiano de transformación
Elemento de Línea Elemento de Área Elemento de Volumen
 u3
z
h3du3e3

h3du3e3
u1
dz
ds P  x, y , z  
h2 du2e2
P

P h2 du2e2
u1 u2

h1du1e1  u2
dA3 u1 h1du1e1
dy dV
y
dx dA1  h2 h3 du 2 du3
dA2  h1 h3 du1 du3 dV  h1h2 h3 du1du2 du3
x
dA3  h1 h2 du1 du 2

ds 2  h1 du1  h2 du2  h3 du3

2 2 2 2 2 2

Considerando las ecuaciones:

x  x(u1 , u2 , u3 ); y  y(u1 , u2 , u3 ); z  z(u1 , u2 , u3 )

x y z
u1 u1 u1
 x, y , z  x y z  x, y , z 
J    J    h1 h2 h3
 u1 , u 2 , u 3  u 2 u 2 u 2  u1 , u 2 , u 3 
x y z
u 3 u 3 u 3
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Jacobiano de transformación

Elemento de Volumen z
w4
z
4 4
u3 v4
u4 2  0

h3du3e3
1  0
w3 3  0
w1 3 3
 v3
P h2 du2e2 1 1
v1
w2 u3
 u1
u1 h1du1e1 u2 4  0
2 v2 2
dV y y
u2

x
dV  h1h2 h3 du1du2 du3 x

  1   

1 3 
Plano de Proyección  , 
2
  1 

 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Transformación ortogonal de coordenadas en R3
y1
z y1
z 3  
   
         
P x  OP i   x1 i1  y1 j1  z1 k1   i  x1  i1  i   y1  j1  i   z1  k1  i 
2 P        
z1 3 2 x  x1 cos1  y1 cos 2  z1 cos 3
k z1
j1


   
         
j y  OP j   x1 i1  y1 j1  z1 k1   j  x1  i1  j   y1  j1  j   z1  k1  j 
y 3        
k1

y
o o 1
y  x1 cos 1  y1 cos  2  z1 cos  3
i
1
i1

2 
   
         
x z  OP k   x1 i1  y1 j1  z1 k1   k  x1  i1  k   y1  j1  k   z1  k1  k 
x        
1
x1
x1 z  x1 cos 1  y1 cos 2  z1 cos 3

       x cos 1 cos 2 cos 3 x1

OP  x i  y j  z k  x1 i1  y1 j1  z1 k1
y  cos 1 cos  2 cos  3 y1
z cos 1 cos 2 cos 3 z1
Matriz de rotación de coordenadas entre dos sistemas cartesianos ortogonales.
cos 1 cos 2 cos 3
M R  cos 1 cos  2 cos  3
cos 1 cos 2 cos 3
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Transformación ortogonal de coordenadas en R3
y1
z
P
z1
k
j1

j
k1

o y
i
i1

x1

z y1
3
2 P
3 2
z1

3 y
o 1

1
x 2
1
x1
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Integral a través de una superficie
z

   x, y, z  dS
S
n

Vector unitario normal exterior
dS
 z z   z z 
  x ,  y ,1   x ,  y ,1
n   k n   0, 0,1   1
2 2 k n
 z   z   z   z 
2 2 2
 z   z 
2
k      1      1      1
j  x   y   x   y   x   y 
i y
2
 z   z 
2

dy
dx
R
d S  ndS  
k d S  k n dS  dxdy dS 
dxdy
k n
dS        1dxdy
 x   y 

Plano xy: Plano xz: Plano yz:

2
 x   x 
2
 y   y 
2 2

   x, y, z  dS      x, y, z     x, y, z  dS      x, y, z        1dydz
2
 z   z 
2
      1dxdz
   x, y, z  dS      x, y, z 
S R
      1dxdy
 x   y  S R  x   z  S R  y   z 
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Integral a través de una superficie
z

   x, y, z  dS
S
n


dS
w4
z 4
p  x, y , z 
k v4
j u4
i y
w3
dx
dy R 3 v3
1
w2
x u3

2 v2
y
u2

x
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Sistema de Ecuaciones Lineales
a11 x1  a12 x 2  a13 x3    a1n x n  b1
a 21 x1  a 22 x 2  a 23 x3    a 2 n x n  b2
a31 x1  a32 x 2  a33 x3    a3n x n  b3
     
X  A 1 B
a n1 x1  a n 2 x 2  a n 3 x3    a nn x n  bn

La solución de sistemas casi de cualquier número de ecuaciones, es una realidad hoy en día gracias a las
computadoras, lo cual proporciona un atractivo especial a las técnicas de solución mediante métodos directos y
métodos iterativos.

Métodos directos se encuentran:

Métodos iterativos se encuentran:
• Eliminación de Gauss.
• Método de Jacobi.
• Eliminación de Jordán.
• Método de Gauss-Seidel.
• Factorización de matrices LU.
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Interpolación Polinomial
Una función de interpolación es aquella que pasa a través de puntos dados como datos, la enorme ventaja de
aproximar información discreta o funciones complejas con funciones analíticas sencillas, radica en su mayor facilidad
de evaluación y manipulación, situación necesaria en el campo de la ingeniería.

Por ejemplo los monomios (elementos finitos unidimensionales) se constituyen en la familia mas empleada, sus
combinaciones generan aproximaciones del tipo polinomial.

a0  a1 x  a2 x 2    an x n
No X y
1 0 2
2 2 6
3 4 4
4 8 7
5 6 5
6 7 8
7 11 5

Pero para algunos casos también se emplean binomios y trinomios (para elementos finitos bidimensionales y
tridimensionales).
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Aproximación Polinomial
Cuando se necesita un polinomio de interpolación donde los datos son separados de manera no uniforme, además
que los datos proporcionados son exactos y que representan valores de alguna función desconocida se pueden
emplear una aproximación polinomial.

Suponiendo que se dan 𝑁 + 1 puntos como:

x0 , x1 ,  x N
f 0 , f1 , f N

El polinomio de orden N que pasa a través de los N+1 puntos se puede escribir en una serie de potencias como:

g x   a0  a1 x  a2 x 2    a N x N

El ajuste de la serie de potencias a los N+1 puntos dados, da un sistema de ecuaciones lineales:

f 0  a0  a1 x0  a 2 x02    a N x0N
f1  a0  a1 x1  a 2 x12    a N x1N

f 0  a0  a1 x N  a 2 x N2    a N x NN

Donde 𝑎𝑖 los coeficientes pueden determinarse resolviendo la ecuaciones simultaneas. Esto será accesible cuando el
número de puntos sea reducido, pero por fortuna, se pueden aplicar métodos para determinar una interpolación
polinomial sin resolver las ecuaciones lineales cuando el número de puntos sea considerable, estos se describe en los
siguientes puntos.
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Formula de interpolación de Lagrange
La formula de interpolación de Lagrange que para un orden N de puntos es:

x  x1 x  x2 x  x N  x  x0 x  x2 x  x N  x  x0 x  x1 x  x N 1 

g x   f0  f1    fN
x0  x1 x0  x2 x0  x N  x1  x0 x1  x2 x1  x N  x N  x0 x N  x1 x N  x N 1 

El polinomio de interpolación de orden N que se ajusta a N+1 puntos es único. Estos es importante, ya que indica que
todos los polinomios de orden N que se ajustan a un conjunto dado de N+1 puntos son matemáticamente idénticos,
aun cuando sus formas sean distintas. No X y
0 0 -3
1 1 0
2 3 5
3 6 7

g x  
x  1x  3x  6  3  x  0x  3x  6 0  x  0x  1x  6 5  x  0x  1x  3 7 
0  10  30  6 1  01  31  6 3  03  13  6 6  06  16  3

g x   
3 3 3 2 276
x  x  x 3
90 90 90
 Método de los Elementos Finitos por [Link]. Ing. Edwin Félix Achá Colque

Conceptos de Matemáticos
Polinomios de Interpolación de Hermite
Un polinomio se puede ajustar no solo a los valores de la función sino también a las derivadas de los puntos. Los polinomios
ajustados a los valores de la función y su derivada se llaman polinomios de interpolación de Hermite.

Si se conocen los puntos x0 , x1 ,, x N y los valores de la función y de todas sus derivadas hasta de orden p f i , f i ' , f i  p  , i  0,1,, N
El número total de datos es K   p  1N  1 Un polinomio de orden K  1 será:
K 1
g x    a j x j
j 0

se puede ajustar a los K datos, donde es un coeficiente a j . Al igualar la ecuación anterior con los datos, obtenemos un
conjunto de K   p  1N  1 ecuaciones:
g  xi   f i , i  0,1, , N
g '  xi   f i ' , i  0,1,, N

g  p   xi   f i  p  , i  0,1,, N

Los valores de los coeficientes se pueden determinar resolviendo esta ecuación en forma exacta si K es pequeño.

Ejemplo en Vigas:
d 2v M  M 

dv

M
dx   v   dx dx
dx 2
EI dx EI  EI 