UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ

CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Aplicación de métodos numéricos a la resolución de problemas

elásticos planos.

DEFINICIÓN DEL PROBLEMA ELÁSTICO PLANO a) desplazamientos u = f (x, y ) v = g ( x, y )

b) deformaciones ε x (x, y ) ε y ( x, y ) γ xy (x, y )

c) tensiones σ x ( x, y ) σ y ( x, y ) τ xy (x, y )
u 
a) campo de desplazamientos: dado por el vector de desplazamientos: u= 
v 
 ∂u 
 
 ε x   ∂x 
  ∂v 
b) campo de deformaciones: dado por el vector de deformaciones: ε = ε y  = 
 
γ xy   ∂y 
  ∂u ∂v
 + 
 ∂y ∂x 

σ x 
 
c) campo de tensiones: dado por el vector de tensiones: σ = σ y 
τ xy 
 
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

 d 11 d 12 d 13 
d) relación tensión - deformación: según Hooke (o Lamé): σ = D⋅ε
D = d 21 d 22 d 23 
matriz D - matriz de características materiales o matriz de tensión deformación
 d 31 d 32 d 33 
- es simétrica, como consecuencia del Teorema de Maxwell – Betty
- asume distinta formulación para estado plano de tensión o estado plano de deformación

E
d.1) estado plano de tensión: σ z = 0 d 11 = d 22 =
1 −ν 2
d 12 = d 21 = ν ⋅ d 11
E
d 33 = G =
2 ⋅ (1 + ν )

E ⋅ (1 −ν )
d.2) estado plano de deformación: εz =0 d 11 = d 22 =
(1 +ν )⋅ (1 − 2 ⋅ν )
ν
d 12 = d 21 = ⋅ d 11
(1 −ν )
E
d 33 = G =
2 ⋅ (1 + ν )
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

e) planteo del equilibrio: realizado por medio de la aplicación del PRINCIPIO DE LOS TRABAJOS VIRTUALES.

ENUNCIADO: una estructura (o parte de ella) se encuentra en equilibrio bajo la acción de un sistema de fuerzas externas si,
al imponer unos desplazamientos virtuales, compatibles con las condiciones de vínculo, el trabajo realizado por las fuerzas
exteriores sobre los desplazamientos virtuales es igual al trabajo que realizan las tensiones sobre las deformaciones
producidas por los desplazamientos virtuales.

∫∫ (δε x ⋅ σ x + δε y ⋅ σ y + δγ xy ⋅τ xy )⋅ t ⋅ dA = ∫∫ (b x ⋅ δu + b y ⋅ δv )⋅ t ⋅ dA + ∫ (t x ⋅ δu + t y ⋅ δv )⋅ t ⋅ dL + ∑ (Pxi ⋅ δu i + Pyi ⋅ δv i )

A A L i

T.V. INTERNO T.V. EXTERNO

∫∫ (δε x ⋅ σ x + δε y ⋅ σ y + δγ xy ⋅τ xy )⋅ t ⋅ dA trabajo virtual interno realizado por las tensiones sobre las deformaciones
A
virtuales, que son producidas por los desplazamientos virtuales: el término
entre paréntesis es el trabajo virtual por unidad de volumen, y el trabajo total
es la integral del mismo sobre el volumen total del cuerpo.

∫∫ (b x ⋅ δu + b y ⋅ δv )⋅ t ⋅ dA trabajo virtual externo, realizado por las fuerzas por unidad de volumen, sobre los despla-
A
zamientos virtuales: término entre paréntesis, trabajo por unidad de volumen, trabajo total,
integral sobre el volumen.
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

∫∫ (δε x ⋅ σ x + δε y ⋅ σ y + δγ xy ⋅τ xy )⋅ t ⋅ dA = ∫∫ (b x ⋅ δu + b y ⋅ δv )⋅ t ⋅ dA + ∫ (t x ⋅ δu + t y ⋅ δv )⋅ t ⋅ dL + ∑ (Pxi ⋅ δu i + Pyi ⋅ δv i )

A A L i

T.V. INTERNO T.V. EXTERNO

∫ (t x ⋅ δu + t y ⋅ δv )⋅ t ⋅ dL trabajo virtual externo, realizado por las fuerzas actuando en el contorno, sobre los desplaza-
L
mientos virtuales. La integral es curvilínea porque las fuerzas actúan a lo largo del contorno,
y son fuerzas por unidad de área.

∑ (Pxi ⋅ δu i + Pyi ⋅ δv i ) trabajo virtual externo, realizado por las fuerzas puntuales sobre los desplazamientos virtua-
i
les. Se suma la contribución de todas las fuerzas puntuales actuando en el cuerpo.

e.1) formulación matricial del principio de los trabajos virtuales: denominando:

δu  vector de desplazamientos virtuales b x 
δu =   b=  vector de fuerzas másicas
δv  b y 
 δε x 
  t x 
δε =  δε y  vector de deformaciones virtuales t=  vector de fuerzas distribuidas sobre
t y  el contorno
δγ xy 
 
 Pxi 
vector de fuerzas puntuales
δu 
qi =   δu i =  
 Pyi  δv  i
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

∫∫ δε ⋅ σ ⋅ t ⋅ dA = ∫∫ δu ⋅ b ⋅ t ⋅ dA + ∫ δu ⋅ t ⋅ t ⋅ dL + ∑ δu i ⋅ q
T T T T

A A L i

T.V. INTERNO T.V. EXTERNO

Consideraciones:
i) sobre el espesor t : para el caso de estado plano de tensión, el espesor es el real; para el caso de estado plano de
deformación, se considera un espesor unitario.
ii) en la expresión matricial del PTV, las transpuestas de los vectores aparecen debido que el trabajo es una cantidad
escalar.
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Método de Elementos Finitos. Preliminares.

Problemas DISCRETOS vs. Problemas CONTINUOS.

sistemas DISCRETOS

sistemas CONTINUOS
T

concepto de DISCRETIZACIÓN
“hacer alguna aproximación de tal naturaleza que quepa esperar que la misma se
acerque tan estrechamente como se quiera a la solución continua verdadera, conforme u=0
crece el número de variables discretas (que son continuas en el problema original)”
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Formas de encarar la solución de los problemas DISCRETOS y CONTINUOS.

Problemas DISCRETOS
a) establecer relaciones fuerza – desplazamiento para los elementos
b) ensamblar el “conjunto” estableciendo el equilibrio local en los
puntos de unión de los elementos (nudos)
c) resolver para obtener los desplazamientos en los nudos
d) utilizar las relaciones a) para obtener los esfuerzos en elementos

Problemas CONTINUOS
a) dividir el continuo mediante “líneas” o superficies imaginarias en un número finito de elementos
b) Suponer que los elementos finitos están conectados entre sí mediante un número discreto de
puntos (nudos) situados en los contornos de los mismos
T
c) los desplazamientos de los nudos son las incógnitas fundamentales del problema (ídem anterior),
d) se establecen funciones que definan de manera única el campo de desplazamientos dentro de
cada elemento en función de los desplazamientos de los nudos
e) dichas funciones definen unívocamente entonces el campo de deformación, y consecuentemente,
junto con las tensiones iniciales (si existen) y las relaciones constitutivas, definen el campo de
tensión en el interior del elemento y sus contornos u=0

f) se determina un sistema de fuerzas concentradas en los nudos, tal que equilibre las tensiones en
los contornos y las cargas de volumen y superficie
g) ensamblar el “conjunto” estableciendo el equilibrio local en los puntos de unión de los elementos
(nudos)
h) resolver para obtener los desplazamientos en los nudos
i) utilizar las funciones (puntos d) y e)) para determinar tensiones y deformaciones, reacciones de
vínculo, etc.
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

DISCRETIZACIÓN de los sistemas CONTINUOS.

T i

e k

y y j

u=0

X X
v3
3 (x 3 , y 3 ) u3

Elementos Finitos. Definición.

y
v
Un elemento está formulado a partir de los nodos (numeración y coordenadas) y
las conectividades. Las conectividades están referidas a la numeración global v1
(numeración de todos los nodos de la estructura). u
numeración local (se numeran solo los nodos del elemento): 1, 2 y 3. v2
(x, y )
numeración global: i, j y k. u1

1 (x1 , y1 ) u2
e 2 (x 2 , y 2 )
x
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Elementos Finitos. Interpolación del campo de desplazamientos. Determinación de los

campos de deformaciones y tensiones.
Tal como habíamos visto anteriormente, el elemento finito se define por medio de sus coordenadas nodales y su conectivi-
dad. Veamos ahora como podemos introducir un nuevo sistema de coordenadas, que resulte útil para la definición de las
matrices de rigidez elementales, y como se realiza la conexión entre ese sistema, y el general, propio del la formulación del
problema.

3 (x 3 , y 3 )
y k
v
2 ( j = 0, k = 1)
u v

(x, y )
e u

1 (x1 , y1 ) ( j, k )
2 (x 2 , y 2 )
e 1 ( j = 1, k = 0 )
x j
COORDENADAS GENERALES 3 ( j = 0, k = 0 ) COORDENADAS NATURALES
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

UTILIDAD DE LAS FUNCIONES DE FORMA

El método de elementos finitos utiliza esta herramienta (la de las funciones de forma) para representar los desplazamientos
dentro del elemento en términos de los desplazamientos nodales del mismo.
u = F (u1,u 2 ,u 3 ) u = N1 ⋅ u1 + N 2 ⋅ u2 + N 3 ⋅ u3
v = G (v1 ,v 2 ,v 3 ) v = N1 ⋅ v1 + N 2 ⋅ v2 + N 3 ⋅ v3

En el caso del elemento finito de tipo triangular, los desplazamientos nodales u1 , u2 , u3 , en los nodos del elemento, definen
una superficie plana.
u3
El desplazamiento u , de un punto interior al elemento, de
coordenadas naturales ( j, k ) , está representado por la
3 ordenada u , medida desde la superficie del elemento hasta
u1 u
el plano definido por los valores nodales u1 , u2 , u3 .

1 ( j, k ) Evidentemente, el desplazamiento en coincidencia con un

u2 nodo, debe igualarse al valor nodal en el mismo.
Ni A través de ello, de las ecuaciones de interpolación, pueden
k
obtenerse los valores de las funciones de forma para las
2
coordenadas nodales (en el sistema natural).
j
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Evidentemente, las funciones de forma asumirán el valor unitario en correspondencia con el nodo, y cero en los demás.
Por ello, asumen también, geométricamente, superficies planas.

3
N3 = 1

3
N1 = 1

1
2 3
1
2

Ni
k N2 = 1
1
j 2 Cualquier combinación lineal de ellas, será
también una superficie plana. Por ej.:
N1 + N 2 + N 3 = 1
Representa un plano // al elemento, a una
unidad de altura.
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

DEDUCCIÓN DE LAS FUNCIONES DE FORMA EN COORDENADAS NATURALES

Habiendo establecido el elemento en el sistema natural de coordenadas, y con la definición anterior de las funciones de
forma, las mismas pueden obtenerse fácilmente:
N1 = j
N2 = k
N3 = 1 − j − k

En el problema bidimensional, las coordenadas ( x, y ) , se mapean en las coordenadas ( j, k ) , y las funciones de forma
se definen como funciones de las coordenadas del sistema natural.
k
INTERPOLACIÓN DE DESPLAZAMIENTOS
2 ( j = 0, k = 1)
Los desplazamientos dentro del elemento se interpolan mediante las
funciones de forma, en términos de los desplazamientos nodales:
u = N1 ⋅ u1 + N 2 ⋅ u2 + N 3 ⋅ u3
v = N1 ⋅ v1 + N 2 ⋅ v2 + N 3 ⋅ v3
e En términos de las funciones de forma:
u = (u1 − u3 ) ⋅ j + (u2 − u3 ) ⋅ k + u3
1 ( j = 1, k = 0 ) j v = (v1 − v3 ) ⋅ j + (v2 − v3 ) ⋅ k + v3
3 ( j = 0, k = 0 )
En forma matricial, las relaciones de interpolación se escriben como: u = N ⋅ a , donde:
N 0 N2 0 N3 0 u 
N= 1 , u =   , a = [u1 v1 u2 v3 ]
T
 v2 u3
0 N1 0 N2 0 N3 v 
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

CONCEPTO DE REPRESENTACIÓN ISOPARAMÉTRICA

De la misma forma que los desplazamientos en el interior del elemento se interpola en términos de los desplazamientos
nodales por medio de las funciones de forma, se pueden interpolar las coordenadas ( x, y ) de un punto interior, en términos
de las coordenadas de los nodos del elemento, haciendo uso también de las funciones de forma. Esta forma de proceder se
denomina REPRESENTACIÓN ISOPARAMÉTRICA, cuya principal utilidad es simplificar el desarrollo y retener la uniformidad
en la formulación de elementos mas complejos.
Entonces:
x = N1 ⋅ x1 + N 2 ⋅ x2 + N 3 ⋅ x3 x = ( x1 − x3 ) ⋅ j + ( x2 − x3 ) ⋅ k + x3
y = N1 ⋅ y1 + N 2 ⋅ y2 + N 3 ⋅ y3 y = ( y1 − y3 ) ⋅ j + ( y2 − y3 ) ⋅ k + y3

Se tienen entonces las expresiones en términos de las coordenadas naturales, tanto para las coordenadas generales como
para los desplazamientos.

DETERMINACIÓN DEL CAMPO DE DEFORMACIONES

Para obtener las deformaciones, es necesario utilizar las relaciones cinemáticas:
∂u
εx =
∂x De las relaciones de interpolación, se observa que tanto ( x, y ) , como (u, v ) , son funciones de las
∂v coordenadas naturales ( j, k ) . Por tanto:
εy =
∂y u = u ( x( j , k ), y ( j , k ))
∂u ∂v v = v( x( j , k ), y ( j , k ))
γ xy = +
∂y ∂x Utilizando la regla de la cadena para expresar las derivadas, se tiene:
∂u ∂u ∂x ∂u ∂y  ∂u   ∂x ∂y   ∂u 
= ⋅ + ⋅
∂j ∂x ∂j ∂y ∂j  ∂j   ∂j ∂j  ⋅  ∂x 
 =   
∂u ∂u ∂x ∂u ∂y  ∂u   ∂x ∂y   ∂u 
= ⋅ + ⋅
∂k ∂x ∂k ∂y ∂k  ∂k   ∂k ∂k   ∂y 
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

DETERMINACIÓN DEL CAMPO DE DEFORMACIONES (cont.)

La matriz cuadrada de 2x2, se denomina JACOBIANO DE TRANSFORMACIÓN, y se utiliza para vincular ambos sistemas de
coordenadas.
 ∂u   ∂x ∂y   ∂u   ∂x ∂y 
 ∂j   ∂j ∂j   ∂x   
 =  ⋅  ∂u  J =  ∂j ∂j 
 ∂u   ∂x ∂y     ∂x ∂y 
 ∂k   ∂k ∂k   ∂y   ∂k ∂k 
De la ecuación anterior, las derivadas necesarias para el cálculo de las deformaciones se obtienen como:
 ∂u   ∂u 
 ∂x  −1
 ∂j 
 ∂u  = J ⋅  
   ∂u 
 ∂y   ∂k 
De la misma forma, se pueden expresar las derivadas de los desplazamientos v .
 ∂v   ∂v 
 ∂x  −1
 ∂j 
 ∂v  = J ⋅  
   ∂v 
 ∂y   ∂k 
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

DETERMINACIÓN DEL CAMPO DE DEFORMACIONES (cont.)

De las ecuaciones anteriores, podemos determinar: x = (x1 − x3 ) ⋅ j + (x2 − x3 ) ⋅ k + x3
 ∂x ∂y   ∂x ∂y  y = ( y1 − y3 ) ⋅ j + ( y2 − y3 ) ⋅ k + y3
= (x1 − x3 ) = x13 = ( y1 − y3 ) = y13 

J =  ∂j ∂j  =  ∂j ∂j  x13 y13 
  =
 ∂x ∂y   ∂x ∂y y23 
= (x2 − x3 ) = x23 = ( y2 − y3 ) = y23   23
x
 ∂k ∂k   ∂k ∂k 
1  y23 − y13 
⋅ adj (J ) = ⋅ [cof (J )] =
1 1
J −1 = ⋅
T

det J det J det J − x23 x13 

det J = x13 ⋅ y23 − x23 ⋅ y13 = 2 ⋅ A

Por lo tanto, las derivadas necesarias para el cálculo de las deformaciones unitarias, resultan ser:
 ∂u   ∂u   ∂u ∂u 
y23 ⋅ − y13 ⋅
 ∂x  1  y23 − y13   ∂j  1  ∂j ∂k 
 ∂u  = ⋅ ⋅  = ⋅ 
  det J − x23 x13   ∂u  det J − x ⋅ ∂u + x ⋅ ∂u 
 ∂y   ∂k  
23
∂j
13
∂k 

 ∂v   ∂v   ∂v ∂v 
y ⋅ − y ⋅
 ∂x  1  y23 − y13   ∂j  1  ∂j ∂k 
23 13

 ∂v  = ⋅ ⋅   = ⋅  
  det J − x23 x13   ∂v  det J − x ⋅ ∂v + x ⋅ ∂v 
 ∂y   ∂k   ∂j ∂k 
23 13

Las deformaciones unitarias vendrán dadas entonces por:

UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

DETERMINACIÓN DEL CAMPO DE DEFORMACIONES (cont.)  ∂u   ∂u  

y ⋅
∂u
− y13 ⋅
∂u 
 ∂x  1  y23 − y13   ∂j  1  23 ∂j ∂k 
 ∂u  = ⋅ ⋅  = ⋅ 
∂u 1  ∂u ∂u    det J − x23 x13   ∂u  det J − x ⋅ ∂u x13 ⋅
∂u 
εx = = ⋅  y23 ⋅ − y13 ⋅   ∂y   ∂k  
23
∂j ∂k 
∂x det J  ∂j ∂k 
 ∂v   ∂v   ∂v ∂v 
y ⋅ − y13 ⋅
∂v 1  ∂v ∂v   ∂x  1  y23 − y13   ∂j  1  23 ∂j ∂k 
 ∂v  = ⋅ ⋅  = ⋅ 
εy = = ⋅  x13 ⋅ − x23 ⋅    det J − x23 x13   ∂v  det J − x ⋅ ∂v x13 ⋅
∂v 
∂y det J  ∂k ∂j   ∂y   ∂k  
23
∂j ∂k 

∂u ∂v 1  ∂u ∂u  1  ∂v ∂v 
γ xy = + = ⋅  x13 ⋅ − x23 ⋅  + ⋅  y23 ⋅ − y13 ⋅ 
∂y ∂x det J  ∂k ∂j  det J  ∂j ∂k 

Recordando las derivadas de las funciones de interpolación, respecto de las coordenadas naturales, se tiene:
∂u
= (u1 − u3 )
u = (u1 − u3 ) ⋅ j + (u2 − u3 ) ⋅ k + u3 → ∂ j
∂u
= (u2 − u3 )
∂k
∂v
= (v1 − v3 )
v = (v1 − v3 ) ⋅ j + (v2 − v3 ) ⋅ k + v3 → ∂ j
∂v
= (v2 − v3 )
∂k
Reemplazando: ∂u
⋅ ( y23 ⋅ (u1 − u3 ) − y13 ⋅ (u 2 − u3 ))
1
εx = =
∂x det J
∂v
⋅ ( x13 ⋅ (v2 − v3 ) − x23 ⋅ (v1 − v3 ))
1
εy = =
∂y det J
∂u ∂v
⋅ ( x13 ⋅ (u2 − u3 ) − x23 ⋅ (u1 − u3 ) + y23 ⋅ (v1 − v3 ) − y13 ⋅ (v2 − v3 ))
1
γ xy = + =
∂y ∂x det J
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

DETERMINACIÓN DEL CAMPO DE DEFORMACIONES (cont.)

Reordenando las ecuaciones, y escribiendo con la notación siguiente, se tiene:
∂u
⋅ ( y23 ⋅ u1 + y31 ⋅ u2 + u3 ⋅ y12 )
1
xij = − x ji εx = =
∂x det J
yij = − y ji ∂v
⋅ ( x32 ⋅ v1 + x13 ⋅ v2 + x21 ⋅ v3 )
1
εy = =
yij = yik − y jk ∂y det J
∂u ∂v
⋅ (x32 ⋅ u1 + y23 ⋅ v1 + x13 ⋅ u 2 + y31 ⋅ v2 + x21 ⋅ u3 + y12 ⋅ v3 )
1
γ xy = + =
∂y ∂x det J
Escribiendo en forma matricial, se tiene:

 y23 0 y31 0 y12 0

x21  ⋅ [u1 v1 u2
1 
ε = B ⋅a = ⋅ 0 v3 ]
T
x32 0 x13 0 v2 u3
det J 
 x32 y23 x13 y31 x21 y12 

La matriz B, que se denomina matriz gradiente, posee coeficientes constantes, e independientes de las variables de
posición. Por ello, el vector deformación, dentro del elemento, es CONSTANTE.

Según lo que se ha definido anteriormente, el campo de tensiones dentro del elemento, es: σ e = D ⋅ ε e = D ⋅ B e ⋅ a e
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Campos de desplazamientos y de tensiones.

 UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Elementos Finitos. Ecuación de Equilibrio Elemental.

Consideremos que en el elemento anterior actúan fuerzas repartidas por unidad de volumen b , y fuerzas repartidas en el
contorno (por unidad de superficie) t .
Asumamos ahora que el equilibrio solo se satisface en
los nodos del elemento, por lo que existirán fuerzas pun-
v 3 ,V3
tuales actuando en los mismos, que equilibrarán las fuer-
u3 , H 3
y zas debidas a la deformación y al resto de las fuerzas ex-
ty
ternas actuando en el mismo.
Para el cálculo de las fuerzas nodales, aplicamos el PTV
v1 ,V1 tx al elemento:
by v 2 ,V 2
∫∫ δε ⋅ σ ⋅ t ⋅ dAe = ∫∫ δu ⋅ b ⋅ t ⋅ dA + ∫ δu ⋅ t ⋅ t ⋅ dL + ...
T T T

Ae Ae Le

... + ∑3i =1 (δu i ⋅ H i + δv i ⋅ Vi )

bx

u2 , H 2
x

Las fuerzas nodales (o fuerzas nodales de equilibrio) se denotan por H i , Vi ,y δu i , δv i representan los desplaza-
mientos virtuales de los nodos del elemento.
Las fuerzas nodales de equilbrio se obtienen entonces despejando de la ecuación de los trabajos virtuales:
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

∫∫ δε ⋅ σ ⋅ t ⋅ dAe − ∫∫ δu ⋅ b ⋅ t ⋅ dA − ∫ δu ⋅ t ⋅ t ⋅ dL = [δa e ] ⋅ q e
T T T T
En la expresión anterior:
Ae Ae Le
v 3 ,V3 δu1   H1 
 δv  V 
u3 , H 3  1  1
y ty δu  H 
δa e =  2  qe =  2 
δv 2   V2 
δu 3  H 3 
   
v1 ,V1 δ
 3 
v  V3 
by tx
v 2 ,V 2
Teniendo en cuenta que: ue = N ⋅ae
εe = B ⋅ae
bx
u1 , H 1 [δu e ]T = [N ⋅ δa e ]T = [δa e ]T ⋅ N T
u2 , H 2
x
[δε e ]T = [B ⋅ δa e ]T = [δa e ]T ⋅ B T

Reemplazando en la ecuación de trabajos virtuales por medio de las ecuaciones anteriores:

∫∫ δa e ⋅ B ⋅ D ⋅ B ⋅ a ⋅ t ⋅ dAe − ∫∫ δa e ⋅ N ⋅ b ⋅ t ⋅ dA − ∫ δa e ⋅ N ⋅ t ⋅ t ⋅ dL = [δa e ] ⋅ q e
T T T T T T T

Ae Ae Le

δa e T ⋅ ∫∫ B T ⋅ D ⋅ B ⋅ a ⋅ t ⋅ dAe − δa e T ⋅ ∫∫ N T ⋅ b ⋅ t ⋅ dA − δa e T ⋅ ∫ N T ⋅ t ⋅ t ⋅ dL = [δa e ]T ⋅ q e
Ae Ae Le

Teniendo en cuenta que los desplazamientos virtuales son arbitrarios, se deduce que:
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

∫∫ B ⋅ D ⋅ B ⋅ a e ⋅ t ⋅ dAe − ∫∫ N ⋅ b ⋅ t ⋅ dA − ∫ N ⋅ t ⋅ t ⋅ dL = q e
T T T
Esta es la ecuación de equilibrio entre las fuerzas nodales
Ae Ae Le
de equilibrio y las fuerzas debidas a:

deformación del elemento

fuerzas másicas
fuerzas de superficie
La expresión anterior puede re expresarse de la siguiente forma, puesto que los desplazamientos nodales son valores cons-
tantes:
K e ⋅ae − fe = qe
 
 ∫∫ B T
⋅ D ⋅ B ⋅ t ⋅ dAe  ⋅ a e − ∫∫ N ⋅ b ⋅ t ⋅ dA − ∫ N ⋅ t ⋅ t ⋅ dL = q e
T T
fe = fb + ft
 Ae  Ae Le
f b = ∫∫ N T ⋅ b ⋅ t ⋅ dA , f t = ∫ N T ⋅ t ⋅ t ⋅ dL
Ae Le

Elementos Finitos. Ecuación de Equilibrio Global.

La suma de las fuerzas nodales de equilibrio extendida a todos y cada uno de los nodos de la discretización, debe ser
igual a la fuerza nodal exterior: ∑ qi = P j La sumatoria se extiende a todos los vectores de fuerzas nodales de equi-
i
librio de los distintos elementos que comparten el j-ésimo nodo. P j representa el vector de fuerzas puntuales exteriores
actuando en el j-ésimo nodo.
Como resultado del planteo del equilibrio global, se obtiene un sistema de ecuaciones de la forma: K ⋅ a − f = r
Tanto la matriz de rigidez global como los vectores de fuerzas se obtienen a partir del procedimiento de ensamblaje desa-
rrollado para el análisis matricial de estructuras de barras.
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

EJEMPLO DE CÁLCULO MEDIANTE ELEMENTOS FINITOS

CASO A RESOLVER: determinación de los desplazamientos, tensiones y deformaciones en una ménsula sometida a una
carga concentrada aplicada en el extremo.
y P = 5tn

20cm
100cm

PASO 1: discretización del dominio. Se consideran dos elementos finitos triangulares, identificando los nodos y numerando
los elementos.

y
3 4

2
1
x
1 2
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

PASO 2: identificación de coordenadas nodales, conectividades por elemento y generación de grados de libertad.
NODO X Y ELEM NODO i NODO j NODO k
1 0 0 1 1 2 3
2 100 0 2 2 4 3
3 0 20
notas:
4 100 20
•la numeración nodal es en sentido antihorario.
y v3 v 4 •los grados de libertad se numeran correlativamente.

u3 u4

3 v1 2 4 v2

u1 1 u2
x

1 2

y 6 8

5 7

2 2 4

1 1 3
x
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Generación de los grados de libertad por elemento

ELEM GDL 1 GDL 2 GDL 3 GDL 4 GDL 5 GDL 6
1 1 2 3 4 5 6
2 3 4 7 8 5 6

y 6 8

5 7

2 2 4

1 1 3
x

PASO 3: formación de la matriz de características materiales

kg
Para: E = 2.1e6
cm 2
ν = 0.25 E
d 11 = d 22 =
 d 11 d 12 d 13  1 −ν 2 2.24e6 5.60e5 0 
Estado plano de tensión: D = d 21 d 23  D =  0 
kg
d 22 d 12 = d 21 = ν ⋅ d 11 2.24e6
 cm 2
 d 31 d 32 d 33  E  sim 8.40e5
d 33 = G =
2 ⋅ (1 + ν )
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

PASO 4: formación de las matrices gradiente (matrices de deformación – desplazamiento)

y 6 8

5 7
La matriz gradiente se forma a partir de las coor-
2 2 4
denadas de los nodos que limitan el elemento.
1 1 3
x

Procedemos por submatrices, teniendo en cuenta que la matriz gradiente está compuesta de tres submatrices, una por
nodo:
b1 0 b2 0 b3 0 b1 = y 2 − y 3 , b2 = y 3 − y1 , b3 = y1 − y 2
1 
B = [B 1 B2 B3 ] = ⋅ 0 c1 0 c2 0 c 3  c1 = x 3 − x 2 , c 2 = x1 − x 3 , c 3 = x 2 − x1
2 ⋅ Ae 
c1 b1 c2 b2 c3 b3 

 − .01 0 .01 0 0 0
B 3 ] =  0 0 .05
b1 = y 2 − y 3 , b2 = y 3 − y1 , b3 = y1 − y 2
B = [B 1
1
1 B2 − .05 0 0
cm c1 = x 3 − x 2 , c 2 = x1 − x 3 , c 3 = x 2 − x1
− .05 − .01 0 .01 .05 0 

 0 0 .01 0 − .01 0 
B = [B 2 B 3 ] =  0 0 
1
2 B4 − .05 0 .05 0
cm
− .05 0 .05 .01 0 − .01
b2 = y 4 − y 3 , b4 = y 3 − y 2 , b3 = y 2 − y 4
c 2 = x3 − x 4 , c 4 = x 2 − x3 , c3 = x 4 − x 2
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

PASO 5: planteo de la ecuación de equilibrio elemental

La forma siguiente no es la que en general se si-
y 6 8
gue en la práctica. No obstante, se demuestra que
5 7 se llega al mismo resultado por ensamble directo.

2 2 4

1 1 3
x

Se plantea la ecuación: K e ⋅ a e − f e = q e para todos los elementos de la malla (en este caso dos).
 
1 K 1 =  ∫∫ B T ⋅ D ⋅ B ⋅ t ⋅ dA1  = B T ⋅ D ⋅ B ⋅ t ⋅ ∫∫ dA1 = B T ⋅ D ⋅ B ⋅ t ⋅ A1
 A1  A1

Utilizando la matriz gradiente calculada anteriormente, y procediendo por submatrices:

BT 1   B T 1 ⋅ D ⋅ B1 BT 1 ⋅ D ⋅ B 2 BT 1 ⋅ D ⋅ B 3 
   
K 1 = B T ⋅ D ⋅ B ⋅ t ⋅ A1 = B T 2  ⋅ D ⋅ [B 1 B2 B 3 ]⋅ t ⋅ A1 = B T 2 ⋅ D ⋅ B 1 BT 2 ⋅ D ⋅ B 2 B T 2 ⋅ D ⋅ B 3  ⋅ t ⋅ A1
B T 3  B T 3 ⋅ D ⋅ B BT 3 ⋅ D ⋅ B 2 BT 3 ⋅ D ⋅ B3 
   1 
De este modo, cada elemento ij-ésimo de la matriz de rigidez, tiene la forma siguiente: k ij = B i ⋅ D ⋅ B j ⋅ t ⋅ Ae
T

El vector de fuerzas nodales equivalentes tiene la forma siguiente:

 F1 
H 
q 1 = F2  Fi =  i 
F3   Vi 
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Los vectores de fuerzas (tanto de volumen como de superficie) son nulas, ya que no hay actuación de ninguna de ellas:

f b = ∫∫ N T ⋅ b ⋅ t ⋅ dA = 0 (**)
Ae

f t = ∫ N T ⋅ t ⋅ t ⋅ dL = 0 (*)
Le

Supongamos, como ejemplo, que actúa una fuerza de superficie, uniformemente distribuida a lo largo del lado que une los
nodos 1 y 2. El vector de carga tendrá la forma siguiente:
t x   0  Debido que el equilibrio se establece solo en los nodos, con conocer cuales son las fuerzas que actúan
t= = 
t y  t y  en los nodos debido a la carga distribuida, ahorramos el problema de resolver la integral (*).
Usando las ecuaciones de la estática:

 f t1 
 f t1 x   0   f t2x   0   f t 3 x  0 
f t = f t 2  , f t1 =  = , f2 =  = , f3 =  = 
y
f t 3   f t 1 y (

  t y ⋅ L)
12 ⋅ t 2 
  f t 2 y
 (
  t y ⋅ L )
12 ⋅ t 2 
  f t 3 y  0 
6

2 4

1 1 3
x

ty
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Supongamos ahora, que actúa una fuerza distribuida por unidad de volumen en el elemento. Una fuerza típica de esta clase
es el peso propio. En lugar de resolver la integral (**), un razonamiento de sentido común nos dice que si la fuerza está uni-
formemente distribuida, la resultante debe distribuirse por partes iguales en todos los nodos del elemento (recordar que el
equilibrio se plantea en los nodos). Para el caso del peso propio, el vector de carga toma la forma siguiente:

b x   0  2 4
b= = 
b y  − γ  1
1 3
x

 f b1 
 f b1x     f b2 x     f b3 x   
f b = f b 2  , f b1 = 
0 0 0
 = =   = =   =
3  3  3
, f , f
 f b1 y  (− γ ⋅ Ve )  f b 2 y  (− γ ⋅ Ve )  f b3 y  (− γ ⋅ Ve )
b2 3
f b3 
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Volviendo al planteo de la ecuación de equilibrio elemental, el vector de los desplazamientos nodales está dado por:
y v3 v4
 u1 
u3 u4 v 
 a1   
1
u 
v1 2 v2 a 1 = a 2  =  2 
a 3   2 
v
u1 1 u2 u 3 
x  
 v 3 

 B T 1 ⋅ D ⋅ B1 BT 1 ⋅ D ⋅ B 2 BT 1 ⋅ D ⋅ B 3   a 1   F1 
 T 
1 K 1 ⋅ a1 = q1 B 2 ⋅ D ⋅ B 1 BT 2 ⋅ D ⋅ B 2 B T 2 ⋅ D ⋅ B 3  ⋅ (t ⋅ A)1 ⋅ a 2  = F2 
B T 3 ⋅ D ⋅ B BT 3 ⋅ D ⋅ B 2 BT 3 ⋅ D ⋅ B 3  a 3  F3 
 1 

De la misma forma se plantea la ecuación de equilibrio para el elemento 2:

B T 2 ⋅ D ⋅ B 2 BT 2 ⋅ D ⋅ B 4 BT 2 ⋅ D ⋅ B 3  a 2  F2 
 T 
2 K 2 ⋅a2 = q2 B 4 ⋅ D ⋅ B 2 BT 4 ⋅ D ⋅ B 4 B T 4 ⋅ D ⋅ B 3  ⋅ (t ⋅ A)2 ⋅ a 4  = F4 
B T 3 ⋅ D ⋅ B BT 3 ⋅ D ⋅ B 4 BT 3 ⋅ D ⋅ B 3  a 3  F3 
 2 
 UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

PASO 6: planteo de la ecuación de equilibrio global

Se plantea la siguiente condición, para todos los nodos del elemento: ∑ qi = P j
i
Para ello, tienen que ponerse de manifiesto las fuerzas externas que actúan en los nodos de la malla de elementos finitos.
y
P3 x
2 P4 x
P3 y P4 y
1
x
P1x
P2 x

P1 y P2 y
Las fuerzas externas están dadas por las reacciones de vínculo que ofrece el empotramiento y por la carga externa que actúa
en el extremo del voladizo. Las restantes son nulas: P2 x = P2 y = P4 x = 0

∑i =1 q i = P1 → k 111 ⋅ a1 + k 112 ⋅ a 2 + k 113 ⋅ a 3 = P1

1
1

∑i =1 q i = P2 → k 1 21 ⋅ a1 + k 1 22 ⋅ a 2 + k 1 23 ⋅ a 3 + k 2 22 ⋅ a 2 + k 2 24 ⋅ a 4 + k 2 23 ⋅ a 3 = P2
2
2

∑i =1 q i = P3 → k 131 ⋅ a1 + k 132 ⋅ a 2 + k 133 ⋅ a 3 + k 2 32 ⋅ a 2 + k 2 34 ⋅ a 4 + k 2 33 ⋅ a 3 = P3

2
3

∑i =1 q i = P4 → k 2 42 ⋅ a 2 + k 2 44 ⋅ a 4 + k 2 43 ⋅ a 3 = P4
1
4
 UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

Expresando en forma matricial, y desarrollando: Puede verse que, procediendo por medio del método de en-
 k 111 k 112 k 113 0   a1   P1  samblaje directo, se llega al mismo sistema de ecuaciones
 1     
k 21 k 1 22 + k 2 22 k 1 23 + k 2 23 k 2 24  a 2  P2  que según el planteo que hemos utilizado aquí.
⋅ =
k 1 31 k 1 32 + k 2 32 k 1 33 + k 2 33 2  a   P3 

k 34

3 El método de ensamblaje directo se plantea de la siguiente
   
 0 k 2 42 k 2 43 k 2 44  a 4  P4  manera: la posición ij-ésima en la matriz de rigidez local

y (en numeración global) ocupa la posición ij-ésima en la ma-

triz de rigidez global.
P3 x
2 P4 x
P3 y P4 y
1
x
P1x
P2 x

P1 y P2 y
PASO 7: solución del sistema de ecuaciones
Planteo de las condiciones de vínculo: en correspondencia con el empotramiento, los desplazamientos son nulos, por tanto:
a1 = 0 Esta condición equivale a “eliminar” del sistema de ecuaciones, las filas y columnas 1 y 3, quedando entonces:
a3 = 0 k 1 22 + k 2 22 k 2 24  a 2  P2 
 2
⋅
2  a 
=   → K ⋅a = R
 k 42 k 44   4  P4 
La solución se obtiene como: (K −1
)
⋅ K ⋅ a = K −1 ⋅ R → I ⋅ a = K −1 ⋅ R → a = K −1 ⋅ R
 UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

El sistema a resolver, expresado en forma explícita, es el siguiente:

4.648e6 0 − 4.20e6 − 8.40e5  u 2   0  u 2  − 0.0024122
   v  
 1.1368e7 − 5.60e5 − 1.12e7   v 2   0  
 2  =  − 0.028475  cm
⋅ = →
 4.648e6 1.40e6  u 4   0  u 4   0.0031425 
         
 sim 1 . 1368e 7   4 
v − 5000   4 
v − 0 . 02906 
y
P3 x
2 P4 x
P3 y P4 y
1
x
P1x
P2 x

P1 y P2 y
PASO 8: determinación de las reacciones de vínculo
Para el cálculo, se utilizan las expresiones de equilibrio global, o lo que es lo mismo, con los valores calculados para los
desplazamientos, se resuelve el sistema de ecuaciones completo.
Operando con el sistema global, se obtiene:  P1x   25000 
P1 =   =   kg
P
 1y   6134 . 70 
 P3 x   − 25000 
P3 =   =   kg
P
  
3 y − 1134 . 70 
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

PASO 9: determinación de tensiones y deformaciones

Se utilizan las relaciones entre deformaciones y desplazamientos εe = B ⋅ae
Y las relaciones entre tensiones y deformaciones σe = D⋅B ⋅ae
Se calcula para cada uno de los elementos. La ubicación de los puntos de tensión (o deformación) está en el centro del
elemento.
 a1   − 0.0241
ε 1 = B ⋅ a1 = [B 1 B2 B 3 ]⋅ a 2  = 1e − 3 ⋅  0 
a 3  − 0.2848
1
 a1   − 54.032 
σ 1 = D ⋅ B ⋅ a1 = D ⋅ [B 1 B 3 ]⋅ a 2  =  − 13.508 
kg
B2 2
a 3  − 239.19 cm

a 2   0.3142 
ε 2 = B ⋅ a 2 = [B 2 B4 B 3 ]⋅ a 4  = 1e − 4 ⋅  − 0.2921
a 3  − 0.1286
2
a 2   54.032 
σ 2 = D ⋅ B ⋅ a 2 = D ⋅ [B 2 B 3 ]⋅ a 4  = − 47.839
kg
B4 2
a 3   − 10.806  cm
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

PASO 10: post procesamiento de los resultados: salidas gráficas (malla deformada, contornos de tensiones, etc.)

MODELO DE DEFORMACIÓN
UNIVERSIDAD TECNOLÓGICA NACIONAL FACULTAD REGIONAL PARANÁ
CÁTEDRA DE ELASTICIDAD Y PLASTICIDAD Preparó: Ing. Juan M. Franco

DISTRIBUCIÓN DE σ x