ESTRUCTURAS III

Análisis de estructuras mediante el método de los elementos finitos

INTRODUCCIÓN
Concepto
ESTRUCTURA es el conjunto
j de elementos conectados de forma
f que
q cumplan
p una función
f resistente,,
soportando un estado de cargas en unas condiciones aceptables

Etapas de un proyecto estructural

1
1. Establecimiento del plan de necesidades
2. Definición y análisis de las propuestas más adecuadas
3. Construcción de la solución elegida
Idealización: barras, placas, láminas, sólidos
Elección del material
Régimen de cargas: Estático o DINÁMICO
DISEÑO Y ANÁLISIS DE ESTRUCTURAS
Relaciones constitutivas o de comportamiento
1. Tipología Relaciones de equilibrio: global, parcial, diferencial
Relaciones de compatibilidad
2
2. Definición de un modelo estructural
3. Resolución de un problema matemático Métodos gráficos
4. Validación de los resultados Métodos analíticos
Métodos numéricos: MEF (elementos finitos)

Materiales

Tipologias Métodos
estructurales de análisis

Diseño Modelado Análisis Construcción

ESTRUCTURAS III: Análisis de estructuras mediante el MEF

MEF (Método de los elementos finitos): Método numérico de obtención de soluciones

aproximadas de EDP sujetas a unas determinadas condiciones de contorno

Cuarto curso, Optativa

Objetivos:
Desarrollar los fundamentos teóricos y prácticos del análisis de estructuras mediante el
método de los elementos finitos

 Comprender las bases teóricas del método

método, los elementos y las tipologías de aplicación
 Aplicar el MEF utilizando programas de diseño y cálculo de estructuras por ordenador
 Despertar una visión crítica del alumno sobre los resultados obtenidos
 Iniciar al alumno en el análisis no lineal de estructuras

Organización
 Tipos
Ti de
d clases:
l
Teóricas: 20 horas
Prácticas con ejercicios teóricos : 10 horas
P á ti
Prácticas con programas: 30 h
horas
ESTRUCTURAS III: Análisis de estructuras mediante el MEF

Bibliografía Básica y apuntes

 Concepts and applications of finite element analysis; Cook R., Malkus D., Plesha. M., John Wiley,1989.
 Cálculo de estructuras por el método de elementos finitos; E. Oñate, CIMNE, 1992.
 I t d ti to
Introduction t finite
fi it elements
l t in
i engineering;
i i Ch d
Chandrupatla
tl T.R.,
T R Belegundu
B l d A.
A Prentice
P ti Hall,
H ll 1997.
1997
 The finite element method (fifth ed.). Vol 1: The Basis, Vol2: Solid mechanics; Zienkiewicz, O.C.,
Taylor, R.L, Thomas Telford, 2000.
 y
Linear static finite element analysis. Online training;
g; Anderson,, W.J.,, Automated Analysis
y Corporation
p
5CDs, 1994.
 The Finite Element Method. Linear Static and Dynamic Finite Element Analysis; T.J. Hughes, Prentice-
Hall, 1987.
 Finite Element Procedures; K.J.
K J Bathe,
Bathe Prentice-Hall,
Prentice Hall 1996.
1996
 Introduction to finite element computations; Hinton, E., Owen, D.R.J., Pineridge Press, 1980.
 Cálculo de estructuras mediante el MEF; Romera L.E., Hernández S, Ediciones Torculo, . 2003.
ESTRUCTURAS III: Análisis de estructuras mediante el MEF

Sistema de Evaluación
La nota final se obtiene como promedio de los trabajos prácticos (2) y teóricos (1) realizados, entregados y
presentados tras su corrección, durante el cuatrimestre.

H
Horas de
d Consulta
C lt
En horas de trabajo, y durante las sesiones de practicas en el laboratorio de cálculo de estructuras por ordenador.
Durante el curso se fijarán las horas de tutorías.
Información Adicional:
Se presupone que los estudiantes han utilizado programas de cálculo de estructuras por métodos matriciales.
ftp://ceres.udc.es/asignaturas/ITS_Caminos/Optativas/Estructuras_III
https://campusvirtual.udc.es/moodle: 632011604 - Estructuras III

En la página están disponibles los temas teóricos desarrollados, prácticas y ejercicios adicionales resueltos.

P t l de
Portales d Internet
I t t
www.cimne.upc.es Centro internacional de métodos numéricos en la ingeniería
http://homepage.usask.ca/~ijm451/finite/fe_resources/node15.html International Finite Element Resources
http://www.6dof.com web profesional sobre FEM
http://femci.gsfc.nasa.gov/links.html Conexiones con paginas web FEM interesantes
http://www.cimne.upc.es/iacm/main.htm International Association for Computational Mechanics
http://www.civilfem.com/ Programa Civilfem (Ansys)
http://www.ansys.com/ Programa Ansys
http://www.mscsoftware.com/ Programas Nastran, Patran y Marc
http://robot-structures.com/usa/index.php Programa Robot (Ansys)
http://www.simulia.com Programa Abaqus
http://www comp engineering com/ Computers and engineering (Sap2000 y otros)
http://www.comp-engineering.com/
http://www.iberisa.com/ Programas y servicios
ESTRUCTURAS III: Análisis de estructuras mediante el MEF
Programa de Estructuras III

1 Introducción al método de los elementos finitos

2 Elementos unidimensionales: barra a axil

3 Elementos finitos en elasticidad bidimensional (I)

4 Introducción al programa Abaqus

5 Elementos finitos en elasticidad bidimensional (II)

6 Aplicación del MEF en problemas térmicos

7 Sólidos
lid ded revolución
l i
8 Elementos finitos tridimensionales

9 Elementos unidimensionales: elemento viga

10 Elementos placa

11 Elementos lámina

12 Elementos lámina axisimétrica

13 Estimación del error, aspectos computacionales y otros temas de interés

14 Introducción al análisis no lineal de estructuras mediante el MEF

15 Introducción al análisis dinámico de estructuras mediante el MEF

T1. Introducción

1.1. Simulación numérica

1.2. Tipos básicos de elementos finitos (EF)
1.2.1. Elementos unidimensionales
1.2.2. Elementos bidimensionales
1.2.3. Elementos tridimensionales
1.2.4. Elementos lámina
1 3 Formulación
1.3. F l ió mediante
di t ell Principio
P i i i ded Trabajos
T b j Virtuales
Vi t l (PTV)
1.3.1. Barra articulada
1.3.2. Barra de tres nudos a axil
1.3.3. Viga plana
1.4. Ejemplo de modelización: mensula
1 5 Ejemplos en teoría lineal
1.5.
1.6. Ejemplos en teoría no lineal
T1. Introducción: Simulación numérica
Proceso de simulación numérica Discretización
PROBLEMA
FÍSICO
Idealización e
identificación
MODELO
40 m
MATEMÁTICO
Métodos
de resolución
y
SOLUCIÓN
x

20 m

uk
uj
ui
u (x,y)

Y j

X i
T1. Introducción: Simulación numérica
Preproceso, análisis numérico y postproceso: proceso iterativo
T1. Introducción: Simulación numérica
Preproceso, análisis numérico y postproceso: proceso iterativo
T1. Introducción: Simulación numérica
Preproceso, análisis numérico y postproceso: proceso iterativo
T1. Introducción: Simulación numérica
Problema Físico Cambio del problema
físico

Modelo matemático
Ecuaciones diferenciales, obtenidas
mediante hipótesis:
- Geométricas
- Cinemáticas Mejora del modelo
- Comportamiento del material
- Propiedades del material
matemático
- Cargas actuantes
- Condiciones de contorno
...

Solución por EF
Solución por elementos finitos del modelo matemático
Selección de:
- Tipos de elementos finitos
- Conexiones entre elementos
- Densidad de la malla
- Parámetros de solución
- Modelización de las cargas
- Modelización de las - Modificación de la malla
condiciones de contorno - Modificación de parámetros
... de solución
...

Estimación de la precisión de la solución

por EF del modelo matemático

Interpretación de resultados Mejora del análisis

Mejora del diseño

Análisis de sensibilidad y
Optimización Estructural
T1. Introducción: Tipos de elementos finitos
1.2.1. Elementos unidimensionales

Barras articuladas planas y espaciales

Vigas planas y espaciales

Vigas planas curvas Gap o tope

T1. Introducción: Tipos de elementos finitos
1.2.2. Elementos bidimensionales
2 grados de libertad por nudo: ux, uy
T1. Introducción: Tipos de elementos finitos
1.2.3. Elementos tridimensionales 3 grados de libertad por nudo: ux, uy, uz
T1. Introducción: Tipos de elementos finitos
1.2.4. Elementos lámina
6 grados de libertad por nudo
T1. Introducción: Formulación de los EF mediante el PTV

En su versión de los movimientos virtuales el PTV establece que en un cuerpo en equilibrio bajo la acción de un
sistema de cargas externas, el trabajo virtual interno es igual al trabajo virtual de las cargas externas para cualquier
campo de movimientos virtuales u compatible.

Denominando u al campo virtual y  al campo de deformaciones virtuales asociado, la expresión general

resultante es

 
T
 d   u T b d   u T t d   u iT p i
  t i

En la expresión anterior el vector tensión  es la tensión producida por las cargas externas, previa a la aplicación
del movimiento virtual, mientras que u es el campo de movimientos virtuales aplicado, y  las deformaciones
virtuales que producen los movimientos virtuales.

Ku = F

s
2 3 u 
U Winterno  
Fext   
s (u)

 
e
1 u  
u
e (u) de (du)
T1. Introducción: Formulación de los EF mediante el PTV
1.3.1. Ejemplo
j p de cálculo de matrices de rigidez
g elementales: barra articulada

Se considera una barra articulada de longitud L, propiedades mecánicas A y del material E, constantes.
L, E, A F2
2
x u2
F1 1
u1

En ejes locales la expresión matricial de equilibrio y la matriz de rigidez elemental, de una barra a axil es:

F1  EA  1  1 u1 
K eu e  F e      
F2  L  1 1 u 2 

Cálculo directo

Se resuelve mediante la aplicación de las ecuaciones de equilibrio, comportamiento y compatibilidad

du L u 2  u1
- Compatibilidad: suponiendo deformación constante   
dx L L
- Comportamiento: u 2  u1
  E  E
L
u u
N  A  EA 2 1
L

- Equilibrio:
ilib i F1  EA  1  1 u1 
u 2  u1
F2  F1  N  EA      
L F2  L  1 1 u 2 
T1. Introducción: Formulación de los EF mediante el PTV
1.3.1. Formulación de un elemento finito barra articulada: Principio de los trabajos virtuales

Considerando la barra articulada aislada y en su estado deformado, sometida a la acción de la fuerza F = (F1, F2) en
sus extremos y con movimientos u = (u1, u2), si le aplicamos unos movimientos virtuales u1 y u2, la expresión del
PTV es:
Wf int  Wf ext     d  F1 u1  F2 u 2
e

Siendo  las deformaciones virtuales producidas por le campo de movimientos virtuales y  la tensión de la barra.

Desarrollando el término de trabajo interno:    d  

e
Le
  EAdx

Para resolver la integral anterior es necesario expresar las deformaciones en función de los desplazamientos mediante
la ecuación de compatibilidad:  =du/dx.

Suponiendo una variación lineal de los desplazamientos en la barra (coherente con la hipótesis de que las
deformaciones son constantes), expresamos el movimiento u de cualquier punto de la barra en función de los
movimientos nodales mediante las funciones de interpolación N1 y N2 siguientes.

x N2(x)
N1  1  N1(x) 1 2
x x
L u ( x )   N i u i  N1 ( x )u1  N 2 ( x )u 2  (1  )u1  u 2
x 2 i 1 L L
N2  1 x
L 1
T1. Introducción: Formulación de los EF mediante el PTV
1.3.1. Formulación de un elemento finito barra articulada
2
x x
u ( x )   N i u i  N1 ( x )u1  N 2 ( x )u 2  (1  )u1  u 2
i 1 L L
Luego: du u 2  u1
 
dx L
d( u ) d( N1u1  N 2 u 2 ) dN1 dN 2 u  u1
    u1  u 2  2
dx dx dx dx L

Sustituyendo las expresiones anteriores en el PTV se obtiene:

u2  u1 u2  u1
 e(
L L
)(
L
) EAdx  F1 u1  F2 u2

Integrando:
EA EA
u1 ( u1  u 2 )  u 2 ( u 2  u1 )  F1 u1  F2 u 2 u1 , u 2
L L

Puesto que la ecuación se debe cumplir para cualquier movimiento virtual compatible:

 EA
F1  L ( u1  u 2 ) F1  EA  1  1 u1 
      
F  EA ( u  u ) F2  L  1 1 u 2 
 2 L
2 1
T1. Introducción: Formulación de los EF mediante el PTV
1.3.2. Barra a axil de 3 nudos
(  1)
N1(x) N2(x) N3(x) N1 (  ) 
1(x1) 2 (x2) 3(x3) 2
x 1 1 1
N 2 ()  1   2
a1 a2 a3 x (  1)
1(x=-1) 2(x=0) 3(x=1) N 3 () 
2
u ( x )  c1  c 2 x  c 3 x 2
3
u ( x )   N i a i  N1 ( x )a 1  N 2 ( x )a 2  N 3 ( x )a 3
i 1

a e  a1,,a2 ,,a3 
du du d 2 du 2 d ( Na e ) 2  dN1 dN 2 dN 3  e
ε     e , , a
dx d dx Le d Le d L  d d d 
2  1  2 1  2 
ε  Be ae luego Be  B1 , B2 , B3   e  , ,  2 
L  2 2 
  D ε  EB e a e

 7 8 1  e
EA
K e   B e D B e Adx  e  8 16  8
T

3L
Le  1  8 7 
T1. Introducción: Formulación de los EF mediante el PTV
1.3.3. Viga plana
w1 w2
N1 () 
1
4

2  3   3  N 2 () 
1
4

1     2  3 
q1 q2
Le N1(x) N2(x)
dN1(x)/dx = 0
x 1 (x1) 2 (x2) 1
dN2(x)/dx = 0

1
1(x = -1) x=0 2(x = 1) 1(x = -1) x=0 2(x = 1)

N3(x) N4(x)
dN3(x)/dx = 0
dN4(x)/dx = 0
1
1
1(x = -1) x=0 2(x = 1) 1(x = -1) x=0 2(x = 1)

N 3 () 
1
4

2  3   3  N 4 () 
1
4

 1     2  3 

 6 
 (Le ) 2 
 
 3ξ  1   12 6L  12 6L 
1  e   2 
EI  6L 4L  6L 2L 
2
L  6 3ξ  1  6ξ 3  1 Le
K  EI y  
e
  e 2 , e , e 2 , e  d  3
1 
 6 ξ   (L ) L (L ) L  2 L  12  6L 12  6L 
 ((L ) 
e 2  
 6L 2L  6L 4L2 
2
 
 3  1
 Le 
T1. Introducción: Ejemplo de modelización
— Mensula:
El soporte metálico de la Figura esta cargado en un extremo y unido mediante dos tornillos a un pilar de acero.
Se desea calcular la flecha y el momento en la sección AA, utilizando distintas modelizaciones.

P = 1000 N
L = 27.5 cm
A
r = 0.5 cm
n E = 2x107 N/cm2
h t  = 0.3
r
A Límite elástico: y = 2400 Kg/cm2
L p h = 6 cm
y,v
espesor t = 0.4 cm

x,u
Modelo 1: Viga en voladizo empotrada, de longitud L+ r y con carga en un extremo.

Hipótesis consideradas:
- Comportamiento como viga incluyendo la deformación por cortante
- Material elástico lineal (no se producen plastificaciones para la carga aplicada)
- Pequeños desplazamientos y deformaciones
- Problema estático (se supone una aplicación lenta de la carga)
- No se consideran movimientos por holgura en los tornillos (empotramiento)

M
M A  PL  27500 Ncm  x max   1146 Kgg / cm 2
W
P( L  r ) 3 P( L  r )
p    0.053 cm
3EI 5
AG
6
T1. Introducción: Ejemplo de modelización
Modelo 2: Elasticidad 2D, estructura en tensión plana.

Hipótesis consideradas:
- Las únicas tensiones no nulas son: x, y, xy
- Material elástico lineal (no se producen plastificaciones para la carga aplicada)
- Pequeños
P ñ desplazamientos
d l i y deformaciones
d f i
- Problema estático (se supone una aplicación lenta de la carga)
- No se consideran movimientos por holguras de los tornillos (empotramiento)
  x  xy
M d l matemático:
Modelo t áti  x  y  0  n  0
Equilibrio interno y en el contorno: 
 
  xy  y  t  0
 x    0
1  0   x   x y
  E    
C
Comportamiento:
t i t  y    1
1   2 
0   y 
  0 0 (1   ) / 2   xy 
 xy   
u v u v
C
Compatibilidad:
tibilid d x  y   xy  
x y x y

Modelo de EF:
T1. Introducción: Ejemplo de modelización
Resultados máximos: p  0.062 cm  x max  1774 Kg / cm 2 M A  26600 Ncm

Ell desplazamiento
d l i vertical
i l máximo
i crece un 17% en ell modelo
d l bidimensional
bidi i l respecto all modelo
d l de
d barra
b a
flexión.

La tensión horizontal crece en un 55% debido al efecto de concentración de tensiones en la zona de acuerdo en las
proximidades
i id d ded la
l unión.
ió La tensión
ió máxima
á i ded Von Mises i / 2) es inferior
(18 8 Kg/cm
(1858 i f i a lal de
d plastificación
l ifi ió (2600
2
Kg/cm ) luego es correcta la hipótesis de material lineal.

Tanto los movimientos máximos (uymax = 0.22L), como las deformaciones máximas (max = x = 0.00049) cumplen
l hipótesis
la hi ót i de
d ser pequeños.
ñ
T1. Introducción: Ejemplos
T1. Introducción: Ejemplos

Uy max
(mm)
07
0.7
Sy max
0.6 (Kg/cm2)
0.5 2500
0.4 Mallas uniformes de
elementos C4 2000
0.3
Valor de referencia
0.2 Uy = 0.6571 mm 1500 Mallas uniformes de
elementos C4
0.1
1000 Valor de referencia
0 Sy = 2304 kg/cm2
0 200 400 600 800 1000 Nd
500

0
Nd
0 200 400 600 800 1000
T1. Introducción: Ejemplos
T1. Introducción: Ejemplos
T1. Introducción: Ejemplos
T1. Introducción: Ejemplos
T1. Introducción: Ejemplos
T1. Introducción: Ejemplos
T1. Introducción: Ejemplos
 T1. Introducción: Ejemplos
— Estructura de edificación:

Planta tipo
Alzado transversa (izda), Alzado frontal (c entro) y planta estandar (dc ha)

Modelo de láminas Modelo de barras

SC:
SC Centro
C t de
d esfuerzos
f cortantes
t t
CG1: Centro de gravedad de los muros
portantes
CG2: Centro de gravedad del forjado
CG3: Centro de gravedad de las cargas
permanentes y sobrecargas

Propiedades de las secciones

Modelos de lám inas y barras
T1. Introducción: Ejemplos

Esfuerzos de m em brana verticales en los m uros (dc ha), y esfuerzos de m em brana y de flexion en la
prim era planta (izda).
T1. Introducción: Ejemplos

z
r
j
y
x
T1. Introducción: Ejemplos en teoría no lineal
Los problemas estructurales se suelen suponer lineales y estáticos, es decir
- Materiales elásticos lineales
- Cargas de aplicación lenta
- Condiciones de contorno constantes
- Pequeños movimientos y pequeñas deformaciones
La modificación de las condiciones anteriores provoca comportamientos no lineales, como los siguientes:

No linealidad del material: Relación  no lineal. El material presenta elasticidad no lineal, plasticidad, fluencia,
termoplasticidad viscoelasticidad,
termoplasticidad, viscoelasticidad fractura...
fractura

No linealidad geométrica: Se producen grandes desplazamientos en la estructura, con o sin grandes deformaciones,
lo que provoca una flexibilización o una rigidización de la misma respecto al comportamiento lineal.

p P
Lineal No lineal
F
No lineal
u u Lineal

u
T1. Introducción: Ejemplos en teoría no lineal
Placa biempotrada lineal y con no linealidad geométrica

E = 2.1e8 KN/m2
= 0.3 Resultados teoría lineal: uymax = -5.28 cm.
Resultados teoría no lineal: uymax = -2.68
2 68 cm.
cm

Curva flecha central/fracción de carga aplicada, en el análisis no lineal

T1. Introducción: Ejemplos en teoría no lineal
Placa biempotrada lineal y con no linealidad geométrica

Teoría lineal: sólo flexión

Teoría no lineal: flexión y membrana

T1. Introducción: Ejemplos en teoría no lineal
Estructura de cables

Edificación
Cargas desplazamientos
p

Respuesta
Parámetro de carga elástica

Carga máxima

Límite elástico

Parámetro de desplazamiento
T1. Introducción: Ejemplos en teoría no lineal

1 T/m

5m

30 m

Modo de pandeo 1: cr1 = 0.289 Modo de pandeo 2: cr1 = 0.66

T1. Introducción: Ejemplos en teoría no lineal
No linealidad de contorno: Se debe a la alteración de las condiciones de contorno de la estructura en función de sus
d l
desplazamientos.
i En ell caso estático
i se denomina
d i contacto, y en ell caso dinámico
di i como en ell choque
h d dos
de d edificios
difi i
por sismo, impacto.

Detalle de deformaciones en la base

M ll y carga hidrostática
Malla hid táti Tensión Sy en la base
T1. Introducción: Ejemplos en teoría no lineal
T1. Introducción: Ejemplos en teoría no lineal
Proceso constructivo

Resultados del análisis lineal:

LC 1 (Carga perm.): Uy1 = -1.07 cm, M1 = 48.4 KNm, M2 = 51.7 KNm, M3 = 43.2 KNm
LC 2 (Sobrecarga ): Uy1 = -1.4 cm, M1 = 47.7 KNm, M2 = 49 KNm, M3 = 76 KNm

Resultados del análisis no lineal:

Fase 1 (Planta 1, t =1): M1 = 32.1 KNm
Fase 2 (Planta 2, t =2): M1 = 68.9 KNm
Fase 3 (Planta 3, t =3): Uy1 = -1.8 cm, M1 = 91.5 KNm, M2 = 43.6 KNm, M3 = 8.56 KNm
Fase 4 (Sobrecarga, t = 4): Uy1 = -3.2 cm, M1 = 133 KNm