Método de los Elementos Finitos

2022/23
Presentación de la asignatura
TFG ‐ Diseño generativo 3D MEF ‐ relleno de Metamaterial 1D‐3D ‐ Resultados MEF NASTRAN

https://twitter.com/gamosinosETSIAE
https://es.linkedin.com/in/gamosinos‐etsiae‐upm‐3967731b3
http://gamosinosetsiae.blogspot.com
https://www.youtube.com/channel/UC8hGeJo7aE1Jfnh‐TH7vrKw Optimización topológica funcional
en hueso 3D MEF
MEF 2022/23 ‐ 1
 Profesorado de MEF

• Profesorado:
Francisco Javier Montáns Leal, CU
Miguel Ángel Sanz, TU (coordinador MEF)
José María Benítez Baena, TU
Luis Saucedo Mora, CD
Laura Moreno Corrales, Ayu
Ismael Ben‐Yelun Insenser, Ayu

• Contacto:
Moodle y Teams
Localización: Departamento de Aeronaves y Vehículos
Espaciales, edificio C (Hangar), despachos: C113 y C114

• Clases magistrales y prácticas:

La clases magistrales en el aula; los ejercicios prácticos se
realizarán con Matlab, ADINA o Nastran, en el aula de
informática (virtual o presencial). Se propondrán ejercicios de
autoevaluación (Matlab Grader – Mathworks, Kahoot).

• Guía Docente y exámenes: web de la ETSIAE

(https://www.etsiae.upm.es/index.php?id=270&L=808%27)

MEF 2022/23 ‐ 2
 R.W. Clough (padre del MEF)

“When I arrived in the summer of 1952 to Boeing, Jon Turner asked me to work on the
vibration analysis of a delta wing structure. Because of its triangular plane form, this
problem could not be solved by procedures based on standard beam theory; so I spent
the summer of 1952 trying to formulate a delta wing model built up of an assembly of
one‐dimensional beams and struts.
However, the results of deflection analysis based on this type of mathematical model
were in very poor agreement with data obtained from laboratory tests of a scale model
of a delta wing. My final conclusion was that my summer’s work was a total failure,
however, at least I learned what did not work.”

F106 Delta Dart Avro 698

MEF 2022/23 ‐ 3
 El método de los elementos finitos

OBJETIVO: modelar el comportamiento de diferentes tipos de sólidos y estructuras, medios continuos

en general a través del método de cálculo más usado en la industria: el método de los elementos
finitos (resuelve numéricamente y de forma aproximada olver ecuaciones diferenciales)

Elementos sólidos: 203300

Elementos planos: 425700 tetras: 5400
quads: 417500 pentas: 26200
trias: 8200 hexas: 176700 MEF 2022/23 ‐ 4
 MEF: el ingeniero decide…
Nivel subestructura 2

¿placa o viga?
Nivel subestructura 1
Nivel elemento
individual
Elementos
Viga
Opción A Todo elementos Placa
Opción B

Elemento
Placa

MEF 2022/23 ‐ 5
 Elementos del MEF

Elementos 3D

Elementos 1D

Elementos 3D

Elementos 1D

Elementos 1D versus elementos 3D sólido en el estudio de

metamateriales de fabricación aditiva (TFG –TFM)
MEF 2022/23 ‐ 6
 Elementos estándar del MEF

Rod‐truss Beam Rod‐truss Beam

2 nodos 3 nodos

Tria 3 nodos Tria 6 nodos

Quad 9 nodos

Quad 4 nodos

Tetra 4 nodos Penta 6 nodos

Hexa 8 nodos Hexa 20 nodos
MEF 2022/23 ‐ 7
 Temario de la asignatura

Tema 1. Introducción al MEF en el cálculo de estructuras aeronáuticas.

1.1. Introducción a la mecánica computacional en medios continuos.

1.2. El cálculo matricial y el M.E.F en elasticidad lineal.
1.3. Conceptos de repaso del cálculo matricial de estructuras.

Bike frame en material compuesto con

elementos 2D shell y 3D solid (TFG)
MEF 2022/23 ‐ 8
 Temario de la asignatura

Tema 2. Herramientas del MEF en el cálculo de estructuras.

2.1. Funciones de Forma (espacios de Sovolev).

2.2. Integración numérica.
2.3. Imposición de condiciones de contorno: condiciones SPC y MPC.
2.3.1. Condensación estática de matrices.
2.3.2. Procedimiento de penalización.
2.3.3. Procedimiento de los multiplicadores de Lagrange.
2.3.4. Procedimientos con elementos master y slave.
2.4. Utilización de mallas no congruentes/no conformes.

Interacción fluido‐estructura 2D Estudio MEF del Springback en el conformado plástico (TFM – TFG)
MEF 2022/23 ‐ 9
 Temario de la Asignatura

Tema 3. Elementos finitos unidimensionales: elementos barra y viga

3.1. Aplicación del MEF a ecuaciones diferenciales de la Elasticidad.

3.2. Elementos unidimensionales de tipo barra.
3.2. Elementos unidimensionales de torsión (árbol de torsión).
3.3. Elementos unidimensionales de tipo viga (vigas de Bernoulli y de Thimoshenko).

Metamateriales de comportamiento auxético (TFG) Modos de Fallo en Metamateriales (TFG)

MEF 2022/23 ‐ 10
 Temario de la Asignatura

Tema 4. Elementos finitos del Medio Continuo: elementos sólidos

4.1. Elementos sólidos en el plano y formulación isoparamétrica: elementos 2D sólido.

4.2. Elementos sólidos 3D: elementos axilsimétricos de revolución y elementos 3D sólido.
4.3. Elementos de orden superior: lagrangianos y serendípitos.

Optimización topológica con

elementos sólidos 3D MEF
(TFG)
Diseño generativo 3D MEF (TFG)
MEF 2022/23 ‐ 11
 Temario de la Asignatura

Tema 5. Elementos finitos estructurales planos: elementos tipo shell (placa y

membrana)

5.1. Formulación de elementos planos: elemento membrana.

5.2. Formulación de elementos planos en flexión: elemento placa.
5.2. Aplicación operativa del elemento shell.

MEF Cajón HTP ‐horizontal tail plane‐ (TFG)

Cono de cola no presurizado en

material compuesto con elementos
2D Shell: resultados MEF (TFG)

MEF 2022/23 ‐ 12
 Temario de la Asignatura

Tema 6. Método de los elementos finitos en estructuras aeronáuticas: prácticas

computacionales con herramientas comunes de cálculo (Patran‐Nastran, Adina,
Hyperworks‐Optistruct).

Validación MEF y análisis no lineal de estructuras en grandes

desplazamientos (TFM)

Ala con comportamiento bimórfico (morphing) en

grandes desplazamientos (TFM)
MEF 2022/23 ‐ 13
 Algo de historia del MEF

• Cuando los aviones eran de tela, y luego de metal: ¿Se calculaban las
estructuras? ¿Cómo?
• “Lösung von Variationsproblemen” W. Ritz 1915…. (y Rayleigh en el s.XIX)
• "Analysis of Continuous Frames by Distributing Fixed‐End Moments”, Método
de Cross 1930.

• Cálculo matricial de estructuras:

• A method for the solution of oscillations problems by matrices, W. J. Duncan
and A. R. Collar, 1934, National Physics Laboratory (NPL‐UK), para estudios de
monoplanos.
• The direct stiffness method of structural analysis, Structural and Materials
Panel Paper AGARD. (M.J. Turner, 1950).

• Método de los Elementos Finitos

• “Rods and Plates. Series occurring in various questions concerning the elastic
equilibrium of rods and plates” B. Galerkin 1915.
• “Variational methods for the solution of problems of equilibrium,” R. Courant
1943.
• Ala delta en Boeing 1952: “Stiffness and deflection analysis of complex
Structures”, M. J. Turner, R. W. Clough, H. C. Martin, and L. J. Topp 1956.
• Clough, R.W., "The finite element method in plane stress analysis“. 1960
• J. H. Argyris, O.C. Zinkiewicz, T. Hughes…

• Códigos de Elementos Finitos

• NASA STRuctrual ANalysis ‐ NASTRAN (1968), ANSYS (1970), SAP (1973), LS
DYNA (1976), ABAQUS (1978), ADINA (1986), Hyperworks‐Optistruct (1995)…
La plataforma Sleippner‐A se hundió en 1991
debido a un error de cálculo con Nastran
MEF 2022/23 ‐ 14
 Los comienzos de los Elementos Finitos
(Boeing/Berkeley 1952‐1960)

Boeing Univ. Berkeley

Página 6
MEF 2022/23 ‐ 15
 Los comienzos de los Elementos Finitos
(Boeing/Berkeley 1952‐1960)

Página 16 Página 18 MEF 2022/23 ‐ 16

 Ordenadores y calculadoras científicas

UNIVAC I (1951): válvulas de vacío IBM PC (1981) 512Kb RAM

IBM 704 (1954) y 709 (1958): FORTRAN IBM Summit (2020) 200 petaFLOPS

MEF 2022/23 ‐ 17
 Ordenadores y calculadoras científicas

Regla de cálculo (hasta la aparición de las calculadoras…)

Olivetti Programma 101 (1965) Hewlett Packard HP9100A (1965) Hewlett Packard HP35 (1973)
MEF 2022/23 ‐ 18
 Bibliografía General MEF

• “Finite Element Procedures”, K. J. Bathe. Prentice Hall 2006.

• “The Finite Element Method Linear Static an Dynamic Analysis”, T.J.R.
Hughes. Dover 2005.
• “The Finite Element Method” O.C. Zienkiewicz & R. Taylor, Ed. Varios
editores y volúmenes.
• “Cálculo Matricial de Estructuras”, E. Alarcón, R. Álvarez, M.S. Gómez.
Editorial Reverte 1990.
• “Cálculo de estructuras por el método de los elementos finitos”, E. Oñate.
CIMNE 1995.
• “Theory of Matrix Structural Analysis”, Przemieniecki. Dover 1985
Capilarización por difusión
(reimpresión de la edición de 1967) MEF en masas tumorales (TFG)

• “A First Course in Finite Elements”. Jacob Fish, Ted Belytschko. Wiley 2007.
• “Finite Element Modelling for Stress‐Analysis”, R.D. Cook, Wiley 1995.

MEF 2022/23 ‐ 19
 Tema 1. Introducción al MEF a través del
cálculo matricial de estructuras
Miguel Angel Sanz Gómez

Los elementos finitos proporcionan una solución aproximada de un modelo matemático

que creemos aproxima la realidad: una realidad que desconocemos…
MEF 2022/23 ‐ 20
Guía Docente de RME cursos pasados

12‐ 18 horas de clase

MEF 2022/23 ‐ 21
 Repaso del Cálculo Matricial

• Vectores de fuerza y desplazamiento. Matriz de rigidez.

• Elemento barra axil.
• Elemento árbol de torsión.
• Elemento viga de flexión.
• Cambio de sistemas de referencias.
• Ensamblaje global.

q0 q0
3 e3 4 e4 5

2
e5
x e2 e8
q x   q0   6
L
L

L
2 8
e6

L/2
e1 e7
1 7

L L

1.5L
2L
2L

MEF 2022/23 ‐ 22
 Movimientos de deformación y de sólido
rígido

Movimientos de sólido rígido: NO cambia alguna dimensión del sólido cuando se mueve, NO se deforma
y
x 

u u Giros rígidos
Traslaciones rígidas


u u

Movimientos de deformación: cambia alguna dimensión del sólido cuando se mueve, se deforma

1
û3  ˆ1 û6  ˆ2
û1 Le
û5  wˆ2
û2  wˆ1 2
û 4

MEF 2022/23 ‐ 23
 Relación entre fuerzas internas del
elemento y desplazamientos nodales

• Elemento barra axil en el eje: movimiento de deformación


y 
x nodos elemento nodos
Fuerzas internas
 fˆ1  EA  1 1 uˆ1 
ˆf  Ke u
ˆ  ˆ     
e 1 1
f
 2  L   uˆ2 
Rigidez: reproducen las fuerzas internas en nodos del elemento
como función de los desplazamientos de esos nodos

• Elemento barra axil en el plano: movimientos de deformación y de sólido rígido

u2 y u2 x '
 fˆ1x  1 0 1 0  uˆ1x 
ˆ     u2 x
f
 1y  EA  0 0 0 0  uˆ1y 
ˆf  Ke u 
ˆ  fˆ   e  1 0 1 0  uˆ2x 
 2x  L    
 fˆ2y   0 0 0 0  uˆ2y 
 
u1 y

u1 x ' u1x

Cambio de base X X'

MEF 2022/23 ‐ 24
 Rotación vectorial y matricial: cambio de
base

Rotación Universal: entre vectores tridimensionales

 cos( x ', x ) cos( x ', y ) cos( x ', z ) 

u RX  X '  cos( y ', x ) cos( y ', y ) cos( y ', z )   R RX ' X  R -1  R T
 
 xx '  cos( z ', x ) cos( z ', y ) cos( z ', z ) 

Vector Tensor‐Matriz
 xx '
u'  RX  X ' u = R u K'  R T K R

Rotación Universal entre vectores en el plano

 cos  sin  0
 RX  X '    sin  cos  0  R
 
 0 0 1 
 cos   sin  0
R X ' X   sin  cos  0  R T
 
 0 0 1 

MEF 2022/23 ‐ 25
 Cambio de base de dos nodos

y' e1
 u    u  
1 e1
 x '    x  
u
ˆ1   u  0    u  0  
1    uˆ1 
x'
1
u  
e1
  y '  1  
y
 
e1 u
ˆ2   u    u   ˆ4 
u
e1  x ' y '   x '     x   xy
   
 4  u  0    u  0   u
2 e2
2
u '   y'  2  x 'y '  y  4 xy
  
 u' u
P
e3
G
AS

y ue'  T  ue = Tue  ue = T Tue'

M


X X '

3 x
Expresión sin condensar
Grados condensados por nodos  e1 
 u    u   
 ux '  
 cos   sin   0 0  
  x    1    ˆe1          
 uy 1   u2    uˆ1  RT 022   u 1     sin  cos  
  0 0    uy ' 1 
      22   
      e1    cos   sin    e1 
u u ˆ
u
  x     7    4 xy  0 RT 22   u 0 0
 22    ˆ2         u  
 u    u    0 0  sin  cos     x '  
  y  4    8    

    u  
    y ' 2 
ue TT
ue ' 

TT
MEF 2022/23 ‐ 26
 Cambio de base de la matriz de rigidez

y' e e
1   fˆ    uˆ  
  1x   e   1x   e
  fˆ1y  
x'  f1   u1 
1 e  uˆ1y   e
e1   ˆ    f   f  uˆ    u   u
e1   f2x    2    2x    2 
  fˆ    uˆ2y  
2 e2  4    2y    
2

 e
 k k12  k13 k14  
P   11    e
e3  K11 K12 
 k21 k22 

k
 23
k24  
  e
 k   K
G

k32  k33 k34  K21 K22 

AS

y 
M

  31 


  
 k41 k42  k43 k44  
  
3 x

 1 0   1 0   k k12  k13 k14  

 RT       11   
e T e 0    0 0  0
 
0 
 R 0   k21 k22 

k
 23
k24  

K  T K 'T      
 0 T   1 0  1 0  0 R  k k32  k33 k34  

R
         31   
  0 0  0
 
0 
  k 41 k 42 

k
 43
k44  

  
  
T
TT
MEF 2022/23 ‐ 27
 Ensamblaje: equilibrio de fuerzas nodales en
coordenadas globales
  fˆ    uˆ  

Fuerzas en los e1   ˆx   e1  x  
1 y
1  f1    fy    u1   uˆy  
   1     1 

x
R1 nodos   ˆ  
e1 f1e1
f
 2    fx   u
 2  ˆ 
 ux  
  fˆ    uˆy  
e1 f 2e1   y 2    2 
2 e2  4
R2 2 f1e 2 f 2e 2  4
  k k12  k13 k14  
  e1   11   
e2  K11 K12   k k22  k k24  
P       21   23 
e3
 K21 K22  k k32  k33 k34  
P   31   
f 2e3  k41 k 42  k 43 k 44  
G

 
AS


M

f1e3


e3
3 e
R3 3 fint  Ke ue
xy

Equilibrio en los nodos:  ˆfext

n
=  ˆfint
n _e
=  Ke uˆe = Ku Ensamblaje nodal

nodo 1: R1  q1  f1e1  e1
K11u e1
ˆ1  K12 ˆ4
u 
  R1  _ _ _ _  u
ˆ 
e2 e2       1
nodo 2: R2  q 2  f2e2  K11 ˆ 2  K12
u ˆ4
u R 2  = _ _ _ _ uˆ2 
 ˆf
ext  K
 ˆ
u
nodo 3: R 3  q 3  f3e3  e3
K11 ˆ 3  K1e23u
u ˆ4  R  _ _ _ _  u
ˆ  
  3    3
  F  _ _ _ _ uˆ
nodo 4: F  q 4  f4e1  f4e2  f4e3 e1
 K21u e1
ˆ1  K22u e2
ˆ 4  K21u e2
ˆ2  K22u e3
ˆ 4  K21u e3
ˆ 3  K22uˆ4    4  fint

↑
Equilibrio de fuerzas externas e internas en los nodos
MEF 2022/23 ‐ 28
 Ensamblaje: equilibrio de fuerzas nodales en
coordenadas globales
Elemento 1. Fuerzas nodales equivalente: fuerzas externas aplicadas en el interior de los elementos, no en los nodos

R1x  q11x  K11 1 1

K12 1
0 0 0 0 K13 1 
K14 
      u1x  … condensación por nodos 
1   
R1y  q1x  K 21
1 1 1 1
K 22 0 0 0 0 K 23 K 24 u1y 
 0      qe1   Ke1 
   0   0 0 0 0 0 0 0 0  u2x  R1  0 0 Ke121   u1  
     1   11    
 0   0   0 0 0 0 0 0 0 0  u2x  0   0   0 0 0 0  u2 
 0      0        
   0   0 0 0 0 0 0 0 0  u3x     0  
0 0 0 0
 u 3  
 0   0   0 0 0 0 0 0 0 0  u3x  f   e1   e1 0 Ke221  u  
   1   1    4 81 q 2 81 K22 0
 88  4 81 
 f4x  q2x  K 31 1
K 32 0 0 0 0 1
K 33 1  u
K 34  4x  
   1   1  
1 1 1  u 4y  
 f4y  q2y  K 41 K 42 0 0 0 0 K 43 K 44

 


Elemento 2. 
 e e
1  0   0   0  0 0
   e2   e2   e2
0 0 
e2

 u1 
 



 
fn  fneq =   Ke u n 
R1 R2  q  f  0 K11 0 K12  u2  
...      1    2     
f1e1 0 

0 u 
e1    0   0  0 0 0  3 
e1 f 2  f  qe 2  f e 2  0 Ke 2 0 Ke222  u  
 4   2  4   21  88  4 18

f1e 2
R2 2 f 2e 2  4 Elemento 3. 
 ˆf  K u
ˆ
e2

 0   0   0  0 0 0 0   u1 


ext 
    0   0  0 0 0  
0   
P  0       u2  
...      e 3    e 3     
f1e3
f 2e3 R
 3  q1  f3  0 0 K11
e3

 f  qe 3  f e 3  0 0 Ke 3
e3 
K12 
3 
u 3 
u 4 

 fint
e3 K22 
 4   2  4   21  88  18 

R3 3

! Todo descrito en coordenadas globales ¡

MEF 2022/23 ‐ 29
 Fuerzas nodales equivalentes

Fuerzas externas aplicadas en el

interior de los elementos
q Fuerzas estáticas de equilibrio en el elemento

q
M1 M2

R1x R2x
R1y R2y

-R1y
feq1 feq2 -R2y
Fuerzas estáticamente equivalentes:
-M1 -M2
trasladadas a los nodos
e
ˆf eq  -R1x -R2x
1 
e e
 . 
 . 
 
ensamblaje
 
  n n
ˆf  ˆf eq =
  Ke uˆn 
ˆf eq 
2 

MEF 2022/23 ‐ 30
 Cálculo de incógnitas

Solución ˆf  K u
ext
ˆ  ˆ  K1 f̂ ext
u

ˆf K es singular, simétrica y
int definida positiva

1 T
 int  uˆ K uˆ  0  uˆ   M , uˆ  0
2

Aplicación de las condiciones de contorno por condensación estática de matrices

R = grado de libertad impuesto o restringido (condición de contorno esencial ‐ Dirichlet)

L = grado de libertad no impuesto o libre (incógnita)

KLL no es singular
 K LL K LR   u L   f L 
K        f L  K LL u L  K LR u R  u L  K LL1  f L  K LR u R 
 RL K RR   u R  f R  incógnita

incógnita f R  K RL u L  K RR u R  !pueden ser 0 o no serlo!

MEF 2022/23 ‐ 31
 Matriz de rigidez directa (Turner 1956)

Matriz de rigidez completa de la estructura

(coordenadas globales)
 K111  K11
5 1
K12 5
K12 0 
• Simétrica: el orden de las cargas no altera el resultado  
• Cuadrada: dimensiones el número de nodos por los  K 21
1
K122  K114
 K11
2 2
K12 4
K12 
grados de libertad nodales.  K 521 K 221 K 222  K 522  K11
3 3
K12 
 
• Definida Positiva: los términos de la diagonal principal  0 K 421 K 521 4
K 22  K 322 
nunca son negativos o nulos: la energía de deformación
directa es siempre positiva.
• Los elementos de la diagonal son los términos directos 
Skyline
de rigidez entre fuerzas y desplazamientos en un
mismo grado de libertad o línea de acción. 11 12 . 0 0 0
• Los elementos fuera de la diagonal son las rigideces  21 22 23 . 0 0 
entre desplazamientos y fuerzas de diferentes líneas 
 . 32 33 34 . 0
de acción.  
• Singular: si hay un movimiento de sólido rígido, el vector 0 . 43 .44 .
de fuerzas en la estructura no se debe ver afectado. 0 0 . . ij 
 
Déficit de rango de la matriz: movimientos rígidos con 0 0 0 . ji ii 
autovalores imaginarios o nulos.
semi-bandwidth

Matriz en Banda: Bandwidth y Skyline

• La ordenación de términos en la matriz de rigidez Una correcta ordenación de nodos hará que los elementos
depende la forma de enumerar los nodos. no nulos se encuentren cerca de la diagonal principal.

MEF 2022/23 ‐ 32
 Elemento 1D barra axil y de torsión

• Elemento árbol de torsión

y
 x1  tˆ1  GJ  1 1 ˆ1 
fˆtor = K uˆ tor
e
   e  1 1   ˆ 
tˆ2  L    2 
z x

K 
 GJ
 x2
tor

• Elemento barra/árbol: rod (Nastran y Optistruct): los desplazamientos axiles y de torsión están
desacoplados en torsión libre
 AE AE 
 e 0  0 
ux1 1  L Le 
 uˆ1   fˆ1  fx GJ GJ 
  ˆ  
 0 0 
 uˆ2  x1  f2  tx1 e 
ˆ
u ... ˆf   e
K  L Le 
uˆ  u 2  fˆ  ... 2
 3 x  3  fx   AE 0
AE
0 
uˆ
 4   1  fˆ4  t 1  Le Le 
x x  GJ GJ 
 0  0 
 L e e 
L  4 4
MEF 2022/23 ‐ 33
 Matriz de rigidez de la viga en el plano
(barra‐viga)

Fuerzas internas sobre los nodos de los elementos Desplazamientos locales

y uˆ3  ˆ1z uˆ6  ˆ2z

x uˆ1  uˆ1x
uˆ5  wˆ2y
uˆ2  wˆ1y
uˆ4  uˆ2x

criterio de signos positivos

• Rigidez: utilizar las fuerzas internas en nodos del elemento como función de los desplazamientos
de esos nodos
 fˆ  uˆ 
 N 1x    1  1  u1x 
nodos
   fˆ  elemento nodos uˆ   
 V1y    2  2  w1y 
M  ˆ  ˆf    K  K    u ˆ   
z
 f3   
 1    
 11   12     1   uˆ3 
 1  
           1z  grados libres de
 N 2x   fˆ4  ˆf  K  K  u ˆ2   uˆ4   u2x  deformación en
V   2 
       21   22  
      w 
ˆ    flexión
 2 y  f
5 uˆ
 5  2y 
M 2z  ˆ     2z 
f
 6 
Condensación estática de matrices uˆ
 6 
MEF 2022/23 ‐ 34
 Matriz de rigidez de la viga en el plano
(barra‐viga)

Fuerzas internas sobre los nodos de los elementos Desplazamientos locales

y uˆ3  ˆ1z uˆ6  ˆ2z

x uˆ1  uˆ1x
uˆ5  wˆ2y
uˆ2  wˆ1y
uˆ4  uˆ2x
Pequeños desplazamientos
  AE   AE 
  0 0   0 0 
  L   L 
  0 12 EI 6EI   0 12 EI 6EI   uˆ  
  1   u  
  3 2   3 2   x  
  L L   L L    uˆ  
6EI 4EI  6EI 2EI   2 
  uy  
 K  K     0  0       
      uˆ3  
 11   12   L2 2
Ke      L   L L          z 1 
  K      AE   uˆ    u  
  K21   22  66   AE 0 0   0 0    
4  x  
 L 
  L   uˆ    uy  
 12EI 6EI   0 12EI 6EI   5    
 0      
  2   u6  
ˆ   z 2 
 L3 L2   L3
L    
 6EI 2EI   0 6EI 4EI 
 0    
 L2 L   L2 L 
MEF 2022/23 ‐ 35
 Matriz de rigidez de la viga en el espacio

y v

x
1 v : vector director de la viga contenido en el plano xy
z

  AE   AE 
  0 0 0 0 0   0 0 0 0 0 
  L   L 
  0 12EI z 6EI z   0 12EI z 6EI z  
0 0 0  0 0 0
  
  L3 L2   L3 L2  
12EI y 6EI y   12EI y 6EI y
 
 0   0   0 
  u  
  0 0
3
0
2  
0
3
0
2   x  
L L L L
 
GJ
 
GJ
  uy  
  0 0 0 0 0   0 0 0  0 0   u  
  L   L 
  6EI y 4EI y   6EI y 2EI y   z  
  0 0  0 0   0 0 0 0   x  
  L2 L   L2 L     
  6EI z 4EI z   6EI z 2EI z    y  
 K  K     0 0 0 0   0  0 0 0   z  
  11 
K  
e  12     
 L2 L   L2 L      1 
  K    AE   AE  
 K21   22  1212  0 0 0 0 0   0 0 0 0 0    ux  
 L   L    u  


12EI z

6EI z   0 12EI z

6EI z    y  
 0 0 0 0 0 0 0   u z  
 L3 L2   L3 L2  
 12EI y 6EI y   12EI y 6EI y 
    
 0 0  0 0 

 0

0 0 0  
  x  
L3 L2  y  
3 2
   L L  
 
GJ
0   0 GJ
0      
 0 0 0 0 0 0 0  

L   L    z 2 
6EI y 2EI y   6EI y 4EI y  
 0 0  0 0   0 0 0 0  
 L2 L   L2 L  
 6EI z 2EI z   6EI z 4EI z  
 0 0 0 0   0  0 0 0  
  L2 L   L2 L  
MEF 2022/23 ‐ 36
 Aplicación: estructuras estéreas 3D

Estructuras articuladas Estructuras reticuladas

z
y
x

z x

MEF 2022/23 ‐ 37
 Objetivos del Cálculo Matricial

¿Qué diferencia el Cálculo Matricial de Estructuras del MEF?

Deformada elástica y Skyline de la matriz de rigidez

con 39 grados de libertad

Deformada elástica y Skyline

de la matriz de rigidez con 12
grados de libertad

Diferencias

MEF 2022/23 ‐ 38