Aae-2. Mef
Aae-2. Mef
Aae-2. Mef
[email protected]
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Contenido específico
2.7 Transformación de coordenadas y ensamblaje de matrices
2.1 Introducción métodos numéricos para la solución de características
ecuaciones diferenciales
2.8 Aplicación computacional del método de elementos finitos en
2.2 Introducción al método de elementos finitos cerchas y pórticos planos
2.3 Formas fuerte y débil de la ecuación diferencial de 2.9 Método de Rayleigh-Ritz y su aplicación en el elemento
gobierno membrana triangular de deformación constante (CST)
2.4 Descripción de los métodos de residuos ponderados y 2.10 Aplicación computacional del método de elementos finitos en
variacionales para solución por elementos finitos membranas
2.5 Método de peso residual de Galerkin y su aplicación en 2.11 Método de Rayleigh-Ritz y su aplicación en el elemento placa
barras en carga axial de Kirchhoff-Love
2.6 Método de Rayleigh-Ritz y su aplicación en barras en 2.12 Aplicación computacional del método de elementos finitos en
carga axial y barras en flexión placas
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Introducción al Método de Elementos Finitos
Ecuación diferencial de gobierno→ Ecuación diferencial de un modelo matemático que permite relacionar las variables dependientes
con las variables independientes (usualmente posición y tiempo). También se denominan ecuaciones de campo
1
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Introducción al Método de Elementos Finitos
Método de Elementos Finitos (FEM or FEA) → Método numérico para resolver ecuaciones Ray Clough (1920-
2016) y Olgierd Cecil
diferenciales con condiciones de contorno (BVP) a través de la discretización del dominio. Zienkiewicz (1921-
Usualmente, la solución es un conjunto de polinomios. 2009)
Sources:
1. http://eccomas.org/vnews/1069/in-memoriam-ray-clough
2. https://www.nap.edu/read/13338/chapter/69 2
3. https://www.youtube.com/watch?v=bYgktqSQACQ
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Introducción al Método de Elementos Finitos
Ecuación diferencial
𝑑𝑃
=𝑅 0 ≤ 𝜃 ≤ 2𝜋
𝑑𝜃 Dominio Nodo
Solución analítica
𝑃 = 𝜃𝑅 0 ≤ 𝜃 ≤ 2𝜋 Elemento
finito
𝑃 = 2𝜋𝑅
Solución numérica
Función de
𝜋 forma
𝑃 = 𝑛2𝑅𝑠𝑖𝑛
𝑛
Discretización
3
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Introducción al Método de Elementos Finitos
1.05 0.2
0.18
1 0.16
Solución numérica
0.14
Error relativo
0.95 0.12
0.1
0.9 0.08
0.06
0.85 0.04
0.02
0.8 0
0 20 40 60 80 0 20 40 60 80
n n
Sources:
1. Author
4
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Elementos Finitos
Real
Geometrías típicas de elementos finitos
Barras→1D Sólidos→3D
FEM
Área→2D
Sources:
1. The Finite Element Method in Engineering-Rao
5
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Elementos finitos
Combinaciones de elementos finitos
Sources:
1. https://perguntas.habitissimo.com.br/pergunta/galpao-pre-moldado
2. BEHAVIOR OF REINFORCED CONCRETE FRAME WITH SHORT COLUMN EFFECT UNDER LATERAL CYCLIC LOADING 6
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Procedimiento general del MEF
1. Plantear el problema físico (Ecuación de gobierno y condiciones
Software comercial de MEF
de contorno)
ENTRE LOS PASOS 1 Y 2 YACE EL MARAVILLO MUNDO DE LA DEDUCCIÓN DE LAS ECUACIONES DE UN ELEMENTO FINITO.
VAMOS A APRENDER A HACER ESTO PARA ALGUNOS ELEMENTOS Y ASÍ RECONOCER LAS HIPÓTESIS INHERENTES A ELLOS
7
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
𝑑𝑈𝑥
EL OBJETIVO PRIMORDIAL DEL MÉTODO ES CALCULAR LOS COEFICIENTES 𝐴𝐸 (𝑥 ) = b
𝑑𝑥 1
8
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Formulaciones fuerte y débil del BVP
𝑑𝑢 𝑑 𝑛 𝑢 𝑑𝑢𝑀𝐸𝐹 𝑑 𝑛 𝑢𝑀𝐸𝐹 Forma residual de la ecuación diferencial
𝐹 𝑥, 𝑢, … =0 𝐹 𝑥, 𝑢𝑀𝐸𝐹 , … ≈0
𝑑𝑥 𝑑𝑥 𝑛 𝑑𝑥 𝑑𝑥 𝑛
𝑑𝑢𝑀𝐸𝐹 𝑑 𝑛 𝑢𝑀𝐸𝐹
→ → 𝑅 𝑥, 𝑢𝑀𝐸𝐹 , … ≈0
2
𝑑 𝑢 𝑑 2 𝑢𝑀𝐸𝐹 𝑑𝑥 𝑑𝑥 𝑛
𝑎2 + 𝑎0 = 0 𝑎2 + 𝑎0 ≈ 0
𝑑𝑥 2 𝑑𝑥 2
Forma fuerte→ Forma natural de la ecuación diferencial Forma débil→ Forma integral reducida de la ecuación diferencial,
escrita a través del residual. también conocida como forma variacional
✓ La solución de prueba debe satisfacer el grado de diferencial ✓ El orden de diferenciabilidad requerido de la ED se reduce
requerido por la ecuación de gobierno ✓ La solución de prueba debe satisfacer solo las EBC
✓ La solución de prueba debe satisfacer todas las condiciones ✓ Las NBC se satisfacen solo en la forma de la ecuación, pero no en la
de contorno (NBC y EBC) solución de prueba.
𝑑𝑢𝐹𝐸𝑀 𝑑𝑛 𝑢𝐹𝐸𝑀 ✓ Se calcula usando integración por partes o con el enfoque del
𝑅 𝑥, 𝑢𝐹𝐸𝑀 , … 𝑛
≈0 funcional de energía
𝑑𝑥 𝑑𝑥
Métodos disponibles
Métodos disponibles
→ Collocation → Galerkin
→ Collocation → Least-squares method → Subdomains → Rayleigh-Ritz
→ Subdomains → Galerkin → Least-squares method
9
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Métodos de Galerkin basado en la forma fuerte
En general, estos métodos buscan que el residual integrado (suma de todos los errores) sea cero. El error se pondera a través de unas funciones
adicionales llamadas funciones de peso W(x). A esta forma se le conoce como forma fuerte integral ponderada.
𝑥2
𝑑𝑢𝑀𝐸𝐹 𝑑𝑛 𝑢𝑀𝐸𝐹
න 𝑾 𝒙 𝑅 𝑥, 𝑢𝑀𝐸𝐹 , … 𝑑𝑥 = 0
𝑥1 𝑑𝑥 𝑑𝑥 𝑛
𝜕𝑢𝐹𝐸𝑀 𝜕𝑢𝐹𝐸𝑀
𝑊𝑖 = 𝑊𝑖 =
U MEF 𝜕c𝑖 𝜕𝑢𝑖
𝑥2 𝑥2
න Φ𝑖 𝑅 𝑑𝑥 = 0 න N𝑖 𝑅 𝑑𝑥 = 0
𝑥1 𝑥1
x
Sources:
1. https://en.wikipedia.org/wiki/Boris_Galerkin
10
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Procedimiento con el método de Galerkin
11
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Exercises
Objective: To apply Finite Element Analysis based on Galerkin method
The bars subjected to the loads shown in the figures have AE=1 and EI=1. Calculate the displacements using the Galerkin approach
5 kN/m 15 kN w w
3m
3m L
L
E
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Métodos de Rayleigh-Ritz basado en la forma débil (variacional)
Funcional (cuadrático)→ Función de funciones
П𝑝 = 𝐸𝑝𝑠 − 𝑊𝑃
𝑰 𝒖 = 𝑩 𝒖, 𝒖 + 𝑳(𝒖) Desplazamiento
U
𝑰 𝒖, 𝒘 = 𝑩 𝒖, 𝒘 + 𝑳(𝒘) Energía de Potencial de trabajo de
deformación las fuerzas externas
𝟎 = 𝑩 𝒖, 𝒘 + 𝑳(𝒘)
1 𝑁2
𝐸𝑝𝑠𝑁 = න 𝑑𝑥 𝑊𝑃 = 𝑃𝑖 𝑈𝑖
2 𝐴𝐸
𝟏
𝑰 𝒖, 𝒖 = П𝒑 = 𝑩 𝒖, 𝒖 + 𝑳(𝒖) RECORDAR QUE ESTA TÉCNICA
𝟐
1 𝑀𝑧2 REDUCIRÁ EL GRADO DE
𝐸𝑝𝑠𝑀 = න 𝑑𝑥 𝑊𝑃 = න 𝑤𝑖 𝑥 𝑈𝑖 (𝑥)𝑑𝑥 DIFERENCIABILIDAD REQUERIDO
2 𝐼𝑧 𝐸
Sources:
1. https://phys.libretexts.org/Under_Construction/Purgatory/2%3A_Applying_Models_to_Mechanical_Phenomena/2.5%3A_Force_and_Potential_Energyht
12
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Método de Rayleigh-Ritz
𝑢 ≈ 𝑢𝑀𝐸𝐹 = 𝑐𝑖 Φ𝑖 𝑢 ≈ 𝑢𝑀𝐸𝐹 = 𝑢𝑖 N𝑖
4. Reemplazar las funciones de prueba (o de forma) y en el funcional de energía potencial
5. Minimizar (derivar) el funcional respecto a los coeficientes
Sources:
1. https://en.wikipedia.org/
13
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Exercises
Objective: To apply Finite Element Analysis based on Rayleigh-Ritz method
The bars subjected to the loads shown in the figures have AE=1 and EI=1. Calculate the displacements using the Rayleigh-Ritz approach
5 kN/m 15 kN w w
3m
3m L
L
E
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
x x
𝑥 𝑥 NOTA
LK 𝑢𝑀𝐸𝐹 = 𝑢𝑘 1− + 𝑢𝑘+1 Para un elemento finito
𝐿𝑘 𝐿𝑘
en formato general los
Nodos 𝑢𝑀𝐸𝐹 = 𝑢𝑘 𝑁𝑘 + 𝑢𝑘+1 𝑁𝑘+1 vectores y matrices
1 k k+1 N+1 obtenidos se conocen
x
𝑢𝑀𝐸𝐹 = 𝑢𝑖 N𝑖 como vectores y
Elemento 1 Elemento k Elemento N Elementos finitos matrices características
14
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
1 1
𝑥 𝑥
𝑁𝑘+1 = 0, Uk
)+ (
𝑁𝑘 = 1 −
𝑢𝑘
( 𝐿𝑘 𝑢𝑘+1 𝐿𝑘
) =
x x x
Lk Lk
¡IMPORTANTE!
→ La solución de prueba es una superposición de las funciones de forma multiplicadas por las soluciones de la ecuación diferencial
(desplazamientos) en los nodos
→ Nótese que para obtener cada función de interpolación se activa de a uno cada grado de libertad
15
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Exercises
Objective: To apply Finite Element Analysis based on Rayleigh-Ritz method
For a general linear-elastic bar subjected to axial load with sectional stiffness AE, calculate the characteristic matrices assuming linear shape
functions and implementing the functional form using the Rayleigh-Ritz method:
Ni n(x) Nj
𝑑 2 𝑈𝑥 𝑥 1 𝑁2 𝑥 𝑥
𝐴𝐸 +𝑛=0 𝐸𝑝𝑠𝑁 = න 𝑑𝑥 𝑈𝑥 = 𝑈𝑖 1 − + 𝑈𝑗
𝑑𝑥 2 2 𝐴𝐸 𝐿 𝐿
0≤𝑥≤𝐿 𝑈𝑥 = 𝑈𝑖 𝑁1 + 𝑈𝑗 𝑁2
E
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Exercises
Objective: To apply Finite Element Analysis based on Rayleigh-Ritz method
For a general prismatic Euler Bernoulli with sectional stiffness IE, calculate the characteristic matrices assuming cubic shape functions and using
the functional form using and the Rayleigh-Ritz Method
w(x)
Mi Mj
Vi Vj
𝑑 4 𝑈𝑦 1 𝑀𝑧2
𝐸𝐼 + 𝑤(𝑥) = 0 𝐸𝑝𝑠𝑀 = න 𝑑𝑥 𝑈𝑦 = 𝑈𝑦𝑖 𝑁1 + 𝜃𝑖 𝑁2 + 𝑈𝑦𝑖 𝑁3 + 𝜃𝑗 𝑁4
𝑑𝑥 4 2 𝐼𝑧 𝐸
0≤𝑥≤𝐿 3𝑥 2 2𝑥 3 2𝑥 2 𝑥 3 3𝑥 2 2𝑥 3 𝑥2 𝑥3
𝑁1 = 1 − 2 + 3 𝑁2 = 𝑥 − + 2 𝑁3 = 2 − 3 𝑁4 = − + 2
𝐿 𝐿 𝐿 𝐿 𝐿 𝐿 𝐿 𝐿
E
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Formulaciones de elementos finitos de barras
𝐿
𝐴𝐸 𝐴𝐸 − න 𝑁1𝑁 𝑛 𝑥 𝑑𝑥
0 0 − 0 0 0
𝐿
𝐴𝐸 𝐴𝐸 𝑛𝐿 𝐿 𝐿 𝐿
− − න 𝑁1𝑁 𝑛 𝑥 𝑑𝑥 = − 12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼 න 𝑁1𝐵 𝑤 𝑥 𝑑𝑥
𝑁𝑖 𝑈𝑥𝑖 2
𝑁𝑗 = 𝐿 𝐿
𝑈𝑥𝑗 + 0
𝑁𝑖 0
𝐿3 𝐿2
0 − 3
𝐿 𝐿2 𝑈𝑥𝑖 0
𝐴𝐸 𝐴𝐸 𝐿
𝑛𝐿 𝐿
− − න 𝑁2𝑁 𝑛 𝑥 𝑑𝑥 = − 𝑉𝑖 6𝐸𝐼 4𝐸𝐼 6𝐸𝐼 2𝐸𝐼 𝑈𝑦𝑖
𝐿 𝐿 2 0 0 − 2 න 𝑁2𝐵 𝑤 𝑥 𝑑𝑥
0 𝑀𝑖 𝐿2 𝐿 𝐿 𝐿 𝜃𝑖 0
𝑁𝑗 = 𝑈𝑥𝑗 +
Solo para n uniforme 𝐴𝐸 𝐴𝐸 𝐿
Elemento cercha 𝑉𝑗 − 0 0 0 0 − න 𝑁2𝑁 𝑛 𝑥 𝑑𝑥
𝐿 𝐿 𝑈𝑦𝑗
0
𝑀𝑗 12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼 𝜃𝑗 𝐿
0 − − 0 −
𝐿3 𝐿2 𝐿3 𝐿2 න 𝑁3𝐵 𝑤 𝑥 𝑑𝑥
6𝐸𝐼 2𝐸𝐼 6𝐸𝐼 4𝐸𝐼 0
𝐿
𝑤𝐿 0 0 − 2 𝐿
𝐿2 𝐿 𝐿 𝐿 න 𝑁4𝐵 𝑤 𝑥 𝑑𝑥
12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼 න 𝑁1𝐵 𝑤 𝑥 𝑑𝑥 =
− 3 0 2 0
𝐿3 𝐿2 𝐿 𝐿2 𝐿
𝑉𝑖 6𝐸𝐼 4𝐸𝐼 6𝐸𝐼 2𝐸𝐼 𝑈𝑦𝑖 𝐵
𝑤𝐿2 Frame element
− 2 න 𝑁2 𝑤 𝑥 𝑑𝑥 =
𝑀𝑖 𝐿2 𝐿 𝐿 𝐿 𝜃𝑖 0 12
𝑉𝑗 = 𝑈𝑦𝑗 +
12𝐸𝐼 6𝐸𝐼 12𝐸𝐼 6𝐸𝐼 𝐿
𝑤𝐿
𝑀𝑗 − 3 − 2 − 2 න 𝑁3𝐵 𝑤 𝑥 𝑑𝑥 =
𝐿 𝐿 𝐿3 𝐿 𝜃𝑗 2 LAS BARRAS AQUÍ DEDUCIDAS ESTÁN
0
6𝐸𝐼 2𝐸𝐼 6𝐸𝐼 4𝐸𝐼 𝐿
− 2 𝐵
𝑤𝐿2 ESCRITAS EN UN SISTEMA DE
𝐿2 𝐿 𝐿 𝐿 න 𝑁4 𝑤 𝑥 𝑑𝑥 = −
0 12 REFERENCIA DONDE “X” VA A LO
Solo para w LARGO DE LA BARRA (SISTEMA LOCAL)
Elemento viga de Euler-Bernoulli con conexiones rígidas
uniforme
16
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
(planteamiento)
ux(x) elemento
Ensamblaje
3. Sumar los términos de las matrices de rigidez y vectores de
Se debe fuerzas de empotramiento que relacionan los mismos GDL La
ensamblar matriz [K] resultante se conoce como matriz completa.
4. Imponer las condiciones de contorno naturales (fuerzas) y
esenciales (desplazamientos)
Solución
1 k k+1
x 𝐹 = 𝐾 𝑈 + 𝐹𝑒
Elemento 1 Elemento k Elemento N 6. Resolver los desplazamientos usando técnicas numéricas o
algebraicas. Una usada comúnmente es la técnica de la matriz
Nota: Todo esto aplica para elementos formulados bajo inversa. Usando los desplazamientos, se pueden obtener las
el mismo sistema de referencia demás fuerzas 𝑈 = 𝐾 −1 𝐹 − 𝐹 𝑒
17
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Exercises
Objective: Apply Finite Element Analysis based on Galerkin method
1. For the BVPs shown below, find a solution using the FEM implementing the weak form and finite elements with linear shape function. Use
AE=1 kN and two finite elements.
5 kN/m 10 kN/m
Figure C. Structure
E
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Transformación de coordenadas en MEF
Procedimiento que permite reescribir los vectores y matrices característicos de un EF de un sistema coordenado (local) a otro
(global)
𝑁𝑖 𝑐 𝑠 0 0 0 0 𝐹𝑋𝑖
𝑉𝑖 −𝑠 𝑐 0 0 0 0 𝐹𝑌𝑖
YG 𝑀𝑖 0 0 1 0 0 0 𝑀𝑖
𝑁𝑗 = 𝐹𝑋𝑗
0 0 0 𝑐 𝑠 0
𝑉𝑗 0 0 0 −𝑠 𝑐 0 𝐹𝑌𝑗
𝑀𝑗 0 0 0 0 0 1 𝑀𝑗𝑍
𝐹𝐿 = 𝑇 𝐹𝐺
18
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Transformación de coordenadas en MEF
𝑁𝑖 𝐹𝑋𝑖 𝑈𝑥𝑖𝐿 𝑈𝑥𝑖𝐺
TRANSFORMACIÓN DE
𝑐 𝑠 0 0 0 0 𝑐 𝑠 0 0 0 0 FUERZAS
𝑉𝑖 −𝑠 𝑐 0 0 0 0 𝐹𝑌𝑖 𝑈𝑦𝑖𝐿 −𝑠 𝑐 0 0 0 0 𝑈𝑦𝑖𝐺
𝑀𝑖 𝑀𝑖 𝐹𝐿 = 𝑇 𝐹𝐺
0 0 1 0 0 0 𝜃𝑖𝐿 0 0 1 0 0 0 𝜃𝑖𝐺
𝑁𝑗 = 𝐹𝑋𝑗 𝑈𝑥𝑗𝐿 = 𝑈𝑥𝑗𝐺
0 0 0 𝑐 𝑠 0 0 0 0 𝑐 𝑠 0
𝑉𝑗 0 0 0 −𝑠 𝑐 0 𝐹𝑌𝑗 𝑈𝑦𝑗𝐿 0 0 0 −𝑠 𝑐 0 𝑈𝑦𝑗𝐺 TRANSFORMACIÓN DE
𝑀𝑗 0 0 0 0 0 1 𝑀𝑗𝑍 𝜃𝑗𝐿 0 0 0 0 0 1 𝜃𝑗𝐺 DESPLAZAMIENTOS
𝑈𝐿 = 𝑇 𝑈𝐺 𝑈𝐿 = 𝑇 𝑈𝐺
𝐹𝐿 = 𝑇 𝐹𝐺
𝑇 [𝐹𝐺𝑒 ] = 𝑇 𝑇 𝐹𝐿𝑒
𝐾𝐺 = 𝑇 𝐾𝐿 𝑇
𝐹𝐿 = 𝐾𝐿 𝑈𝐿 + 𝐹𝐿𝑒
Matriz de rigidez Matriz de
𝑇 𝐹𝐺 = 𝐾𝐿 𝑇 𝑈𝐺 + 𝐹𝐿𝑒 Matriz de rigidez Fuerzas E en Fuerzas E en
en coordenadas transformación
en coordenadas globales locales
globales locales
𝑇 −1 𝑇 𝐹𝐺 = 𝑇 −1 𝐾𝐿 𝑇 𝑈𝐺 + 𝑇 −1 𝐹𝐿𝑒
𝐹𝐺 = 𝑇 𝑇
𝐾𝐿 𝑇 𝑈𝐺 + 𝑇 𝑇
𝐹𝐿𝑒 𝐹𝐺 = [𝐾𝐺 ] 𝑈𝐺 + 𝐹𝐺𝑒
19
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
𝑐 𝑠 0 0 𝐾𝐺 = 𝑇 𝑇
𝐾𝐿 𝑇
𝑇 =
0 0 𝑐 𝑠
𝐴𝐸 𝐴𝐸
− 𝐴𝐿2 2 𝐴𝐿2 𝐴𝐿2 2 𝐴𝐿2
𝐾𝐿 = 𝐿 𝐿 𝑐 + 12𝑠 2 − 12 𝑠𝑐 −6𝐿𝑠 − 𝑐 + 12𝑠 2 − − 12 𝑠𝑐 −6𝐿𝑠
𝐴𝐸 𝐴𝐸 𝐼 𝐼 𝐼 𝐼
−
𝐿 𝐿 𝐴𝐿2 𝐴𝐿2 2 𝐴𝐿2 𝐴𝐿2 2
− 12 𝑠𝑐 𝑠 + 12𝑐 2 6𝐿𝑐 − − 12 𝑠𝑐 − 𝑠 + 12𝑐 2 6𝐿𝑐
𝐼 𝐼 𝐼 𝐼
𝑇 𝐸𝐼
𝐾𝐺 = 𝑇 𝐾𝐿 𝑇 −6𝐿𝑠 6𝐿𝑐 4𝐿2 6𝐿𝑠 −6𝐿𝑐 2𝐿2
2
𝐿3 𝐴𝐿 2 𝐴𝐿2 2
𝐴𝐿 2 𝐴𝐿2
− 𝑐 + 12𝑠 2 − − 12 𝑠𝑐 6𝐿𝑠 𝑐 + 12𝑠 2 − 12 𝑠𝑐 6𝐿𝑠
𝐼 𝐼 𝐼 𝐼
𝑐2 𝑠𝑐 −𝑐 2 −𝑠𝑐
𝐴𝐸 𝑠𝑐 𝑠2 −𝑠𝑐 −𝑠 2 𝐴𝐿2 𝐴𝐿2 2 𝐴𝐿2 𝐴𝐿2 2
− − 12 𝑠𝑐 − 𝑠 + 12𝑐 2 −6𝐿𝑐 − 12 𝑠𝑐 𝑠 + 12𝑐 2 −6𝐿𝑐
𝐿 −𝑐 2 −𝑠𝑐 𝑐2 𝑠𝑐 𝐼 𝐼 𝐼 𝐼
−𝑠𝑐 −𝑠 2 𝑠𝑐 𝑠2 −6𝐿𝑠 6𝐿𝑐 2𝐿2 6𝐿𝑠 −6𝐿𝑐 4𝐿2
20
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
21
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Fuerzas internas en barras
Método de integración Método de integración para
Coordenadas locales
20
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Exercise
Objective: To apply Finite Element Analysis on plane structures
1. For the truss structure shown in Figure, calculate 2. For the beam structure shown in Figure, calculate displacements at all joints
displacements at all joints using FEA. The bars in the using FEA. The beam has rectangular cross section of 30x40 cm and Young modulus
truss have cross section with diameter 50 mm and of 20 GPa. Compare your solution to a finite element analysis software solution.
Young modulus of 200 GPa. Compare your solution to 10 kN/m 10 kN/m
a finite element analysis software solution.
2
2m
3m 3m 3m 3m
10 kN/m 10 kN/m
1m
50 kN
3
1
3m 3m 3m 3m
E
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Exercise
Objective: To apply Finite Element Analysis on plane structures
1. For the truss structure shown in Figure, calculate 2. For the beam structure shown in Figure, calculate displacements in all joints
displacements in all joints using FEA. The bars in the using FEA. The bars in the beam have rectangular cross section of 30x45 cm and
truss have cross section with diameter 50 [mm] and Young modulus of 20 [GPa]. Compare your solution to a finite element analysis
Young modulus of 70 [GPa]. Compare your solution to software solution.
a finite element analysis software solution.
1 30 [kN]
2 [m] 9 kN/m
25 [kN]
19 kN*m 10 kN*m
3
1 [m]
2 [m]
7 [m] 5[m]
2
Figure A. Truss Figure B. Beam
E
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Exercise
Objective: Apply Finite Element Analysis on plane structures
1. For the plane frame structure shown in Figure, calculate displacements at all joints and diagrams of internal forces using FEA. The beams and
columns are made of concrete with E=24 GPa. Columns are 400x400 mm. Beams are 300x400 mm. Consider self-weight of beams and columns
assuming a unit weight of 24 kN/m3. Consider the following scenarios for w and P:
5m
3m
E
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Elementos de área
z
Comportamiento en el plano (membrana)
y
Elementos finitos típicos Elementos estructuras
típicos
fz ✓ CST
✓ Triángulo de 6 nodos ✓ Muros de corte
✓ Cuadrilátero de 4 nodos ✓ Ménsulas
✓ Cuadrilátero de 8 nodos ✓ Vigas de acople
fx ✓ Cuadrilátero de 9 nodos
fy
t ✓ Cuadrilátero con rotación en los
nodos (drilling)
x Comportamiento fuera del plano (placa)
21
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
σyx
Configuración
deformada σxy ϵyy
σxx σxx
σxy
fv ϵxx
Configuración σyx
Y F
de referencia
σyy
22
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
fs 𝑈𝑥 𝐹𝑥
{𝑈(𝑥, 𝑦)} = 𝑈 {𝐹} = 𝐹
X 𝑦 𝑦
23
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
𝐿 𝑇 σ + 𝑓𝑣 = 0 ϵ = [𝐿] 𝑈
σ = [𝐷] ϵ
24
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
25
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
26
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Nodo Uy2
Uy1 x2, y2
2 Ux2
Y
Y
1 x1, y1 Ux1
✓ Elemento triangular
✓ 3 nodos
X
✓ 6 GDL
X ✓ Espesor t
27
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Fx3
Fy2
Fy1
Fx2
fs
fv
Fx1
𝑓𝑣𝑥 𝑓𝑠𝑥 𝐹𝑥
{𝑓𝑣 (𝑥, 𝑦)} = {𝑓𝑠 (𝑥, 𝑦)} = {𝑓𝑝 (𝑥, 𝑦)} = 𝐹
𝑓𝑣𝑦 𝑓𝑠𝑦 𝑦
28
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
29
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
𝜕
𝑈𝑥1
𝑈𝑥1
0
𝑈𝑦1 𝑈𝑦1
ϵ𝑥𝑥 𝜕𝑥 ϵ𝑥 𝑦23 0 𝑦31 0 𝑦12 0
𝜕 𝑁1 0 𝑁2 0 𝑁3 0 𝑈𝑥2 1 𝑈𝑥2
ϵ𝑦𝑦 = 0 ϵ𝑦 = 0 𝑥32 0 𝑥13 0 𝑥21
𝜕𝑦 0 𝑁1 0 𝑁2 0 𝑁3 𝑈𝑦2 γ𝑥𝑦 det[𝐽] 𝑈𝑦2
γ𝑥𝑦 𝑥32 𝑦23 𝑥13 𝑦31 𝑥21 𝑦12
𝜕 𝜕 𝑈𝑥3 𝑈𝑥3
𝜕𝑦 𝜕𝑥 𝑈𝑦3 𝑈𝑦3
𝜕𝑁1 𝜕 𝑥2 𝑦3 − 𝑥3 𝑦2 + y2 − y3 x + x3 − x2 y
+0= = y2 − y3 = 𝑦23
𝜕𝑥 𝜕𝑥 det[𝐽]
30
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
𝐾 = 𝐴𝑡 𝐵 𝑇 𝐷 𝐵
1 ν 0
𝐸 1 𝑦23 0 𝑦31 0 𝑦12 0
ν 1 0
𝐷 = 1−ν 𝐵 = 0 𝑥32 0 𝑥13 0 𝑥21
1 − ν2 det[𝐽] 𝑥 𝑦23 𝑥13 𝑦31 𝑥21 𝑦12
0 0 32
2
𝑇 1 ν 0
1 𝑦23 0 𝑦31 0 𝑦12 0 𝐸 1 𝑦23 0 𝑦31 0 𝑦12 0
ν 1 0
𝐾 = 𝐴𝑡 0 𝑥32 0 𝑥13 0 𝑥21 1 − ν det[𝐽] 0 𝑥32 0 𝑥13 0 𝑥21
det[𝐽] 𝑥 𝑦23 𝑥13 𝑦31 𝑥21 𝑦12 1 − ν2 𝑥32 𝑦23 𝑥13 𝑦31 𝑥21 𝑦12
32 0 0
2
31
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
𝐹 = 𝐾 𝑈 + {𝐹𝑒 }
𝐹𝑥1 𝑈𝑥1
𝐹𝑦1 𝑇 1 ν 0 𝑈𝑦1
1 𝑦23 0 𝑦31 0 𝑦12 0 𝐸 1 𝑦23 0 𝑦31 0 𝑦12 0
𝐹𝑥2 ν 1 0 𝑈𝑥2
= 𝐴𝑡 0 𝑥32 0 𝑥13 0 𝑥21 1 − ν det[𝐽] 0 𝑥32 0 𝑥13 0 𝑥21
𝐹𝑦2 det[𝐽] 𝑥 1 − ν2 𝑈𝑦2
32 𝑦23 𝑥13 𝑦31 𝑥21 𝑦12 0 0 𝑥32 𝑦23 𝑥13 𝑦31 𝑥21 𝑦12
𝐹𝑥3 2 𝑈𝑥3
𝐹𝑦3 𝑈𝑦3
32
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Sources:
1. Three dimensional static and dynamic structural analysis – Edward Wilson
33
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Exercise
Objective: To apply Finite Element Analysis on membranes using CST
1. For the wall structure shown in Figure, calculate displacements at all joints using FEA. The wall is made of concrete with E=20 [GPa], v=0.2
and has thickness t=150 [mm]. Compare your solution to a FEA software solution. Use the following approaches.
a. Euler-Bernoulli beam
b. CST membrane elements. Use N=2
c. CST membrane elements. Use N=8
100 kN 100 kN
100 kN 100 kN
Figures A, B and C.
Y Wall real, CST model
Y and EB model
Y
2m 9m
X X X
150 mm
Sources:
1. https://www.structuralguide.com/shear-walls/
E
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Exercise
Objective: To apply Finite Element Analysis on membranes using CST
1. Concrete corbel shown in Figure is made of concrete with E=24 [GPa], v=0.2 and has thickness t=30 [cm]. Calculate displacements and
reactions. Use one, two and three CST element according to Figures
50 kN 50 kN 50 kN
50 cm 50 cm 50 cm
Y Y Y
40 cm 40 cm 40 cm
X
X X
E
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
dy
y
fsz
x F z dx
fvz
𝑞𝑧 = න 𝑓𝑣𝑧 𝑑𝑧
t σy
τyz
t/2 y
τyx
t/2
x τxz τxy
tdA σx
Solid with thickness t
39
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
40
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
dy 𝜏𝑥𝑦 = 𝜏𝑦𝑥
12
t/2 𝜎𝑥𝑦 = 𝑚𝑥𝑦 𝑧
𝑡3
y τyx
t/2
x mxy τxy Torsional moment mxy is a moment per
unit length
41
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
dy
𝑣𝑥 = න 𝜏𝑥𝑧 𝑑𝑧 𝑣𝑦 = න 𝜏𝑦𝑧 𝑑𝑧
z dx
42
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
12 12 12𝑣𝑥 𝑧𝑑𝑧
𝜕σ𝑥 𝜕τ𝑥𝑦 𝜕τ𝑥𝑧 𝜕𝑚𝑥 𝑧 𝜕𝑚 𝑧
𝑥𝑦 𝑡 3 𝜕 𝜕𝑚𝑥 𝜕𝑚𝑥𝑦
න 𝑡3 + − 𝑡 3
𝑑𝑧 = 0 + − 𝑣𝑥 = 0
+ + +𝑓𝑣𝑥 = 0 𝜕𝑥 𝜕𝑦
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧
12 12 12𝑣𝑦 𝑧𝑑𝑧
𝜕τ𝑦𝑥 𝜕σ𝑦 𝜕τ𝑦𝑧 𝜕𝑚𝑥𝑦 𝑧 𝜕𝑚 𝑧
𝑦 𝑡3 𝜕 𝜕𝑚𝑥𝑦 𝜕𝑚𝑦
+ + + 𝑓𝑣𝑦 = 0 න 𝑡3 + − 𝑡3 𝑑𝑧 = 0 + − 𝑣𝑦 = 0
𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦
𝜕𝑣𝑥 𝜕𝑣𝑦
𝜕 𝜕 𝜕𝑣𝑥 𝜕𝑣𝑦
𝜕τ𝑧𝑥 𝜕τ𝑧𝑦 𝜕σ𝑧 න 𝜕𝑧 + 𝜕𝑧 +𝑓 𝑑𝑧 = 0 + + 𝑞𝑧 = 0
+ + +𝑓𝑣𝑧 = 0 𝜕𝑥 𝜕𝑦 𝑣𝑧 𝜕𝑥 𝜕𝑦
𝜕𝑥 𝜕𝑦 𝜕𝑧
43
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
𝜕 2 𝑚𝑥 𝜕 2 𝑚𝑥𝑦 𝜕 2 𝑚𝑦
+2 + + 𝑞𝑧 = 0
𝜕𝑥 2 𝜕𝑥𝜕𝑦 𝜕𝑦 2
Constitutive equations
𝜎𝑥 − ν𝜎𝑦
𝜖𝑥 = 1 𝜈 0
𝐸 𝜎𝑥 𝜖𝑥
𝐸 𝜈 1 0
𝜎𝑦 − ν𝜎𝑥 𝜎𝑦 = 𝜖𝑦
𝜖𝑦 = 1 − 𝜈2 1−𝜈
𝐸 𝜏𝑥𝑦 0 0 𝛾𝑥𝑦
2
𝐸
𝜏𝑥𝑦 = 𝛾
2(1 + 𝑣) 𝑥𝑦
𝜏𝑧𝑥 = 𝐺𝛾𝑧𝑥
Kirchhoff-Love plates vs Reissner Mindlin plates
𝜏𝑦𝑧 = 𝐺𝛾𝑦𝑧
44
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
45
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
46
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
Exercise
Objective: Apply Finite Element Analysis on plates using Galerkin
1. A solid slab is simply supported in all its edges. The slab is made of concrete with E=20 GPa and v=0.2. Using the FEM with strong
formulation and Galerkin method, calculate de displacement at mid-span. Consider the trial solution presented in the equation.
10 kN/m
z
𝜋𝑥 𝜋𝑦
𝑈𝑍 = 𝑐11 sin sin
𝐿𝑥 𝐿𝑦
4m
y
20 cm
x 5m
E
Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes
x Stiffness matrix
2a
Sources:
1. Structural Analysis with the Finite Element Method - Eugenio Oñate
47
Método de elementos finitos
Prof. David Cotes