ANÁLISIS AVANZADO DE ESTRUCTURAS

Prof. David Sebastián Cotes Prieto

Ingeniero civil-Universidad Industrial de Santander

Magister en ingeniería civil - Énfasis en estructuras y materiales

[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

Ejemplos Método numérico→ Herramienta matemática utilizada

para determinar soluciones aproximadas
𝑑 2 𝑈𝑥 𝑥
𝐴𝐸 +𝑛=0 ✓ Método de Runge-Kutta → ODE
𝑑𝑥 2

𝑑 4 𝑈𝑦 (𝑥) ✓ Método de Diferencias Finitas→ ODE and PDE

𝐸𝐼 +𝑤 =0
𝑑𝑥 4
✓ Método de Elementos Finitos→ ODE and PDE
𝑑 2 𝑈𝑥 𝑑𝑈𝑥
𝑚 + c + kUx = 0 ✓ Método de Elementos de Contorno→ ODE and PDE
𝑑𝑡 2 𝑑𝑡
✓ Método de Volúmenes Finitos→ ODE and PDE
𝜕 4 𝑈𝑧 𝜕 4 𝑈𝑧 𝜕 4 𝑈𝑧 12(1 + 𝑣 2 )
+2 2 2+ =− 𝑤
𝜕𝑥 4 𝜕𝑥 𝜕𝑦 𝜕𝑦 4 𝐸𝑡 3 ¡DO NOT USE A STRUCTURAL ANALYSIS PROGRAM UNLESS YOU FULLY UNDERSTAND
THE THEORY AND APPROXIMATIONS USED WITHIN THE PROGRAM! Edward Wilson

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)

✓ Discretización → División del dominio (región geométrica) en subdominios.

✓ Elemento finito→ Subdominio rodeado por nodos. Las variables en su interior se
describen por funciones de forma o de prueba.
✓ Nodo → Punto de conexión entre elementos finitos
✓ Solución de prueba→ Función asumida que describen los desplazamientos al
interior del elemento finito
✓ Función de forma (interpolación)→ Función parte de la solución de prueba, que
cumple las BC, donde las constantes están en términos de los valores
(desplazamientos) en los nodos
✓ Funciones de prueba→ → Función parte de la solución de prueba, que cumple
las BC, donde las constantes son arbitrarias
✓ Malla → Grupo de elementos finitos y nodos
✓ Condiciones de contorno→ Valores específicos que debe satisfacer la ecuación
diferencial en ciertos puntos (fuerzas de superficie y apoyos)
✓ Error relativo→ Diferencia porcentual obtenida entre las soluciones para
diferentes enmallados propuestos
✓ Convergencia → Cuando el error relativo alcanza un valor objetivo.

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

Curva de convergencia Curva de exactitud

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)

2. Seleccionar el tipo de elementos finitos por usar(1D, 2D y/o 3D)

3. Definir las propiedades geométricas y materiales de los

elementos Preproceso

4. Discretizar el dominio (Enmallar)

5. Ensamblar los elementos finitos (definir conexiones)

6. Imponer las condiciones de contorno

7. Resolver las ecuaciones ensambladas Proceso

8. Reportar los resultados

Postproceso
9. Evaluar la convergencia

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

Deducción de las ecuaciones de los FE→ Solución de prueba del BVP

Ecuación diferencial de gobierno Subdominio Condiciones de borde
𝑑𝑈𝑥
𝑑𝑢 𝑑 𝑛 𝑢 𝑑 2 𝑈𝑥 𝑥1 ≤ 𝑥 ≤ 𝑥2 𝐴𝐸 (𝑥 ) = b
𝐹 𝑥, 𝑢, … =0 𝐴𝐸 +𝑛 =0 𝑑𝑥 1
𝑑𝑥 𝑑𝑥 𝑛 𝑑𝑥 2
𝑈𝑥 (𝑥2 ) = c

Solución de prueba 𝑢𝑀𝐸𝐹 Trial solution → 𝑢𝐹𝐸𝑀 𝑥1 , 𝑥2 … 𝑥𝑚 Tipos de condiciones de contorno

𝑢 ≈ 𝑢𝑀𝐸𝐹 = ෍ 𝑐𝑖 Φ𝑖 Trial functions
✓ Esenciales (EBC) → Relacionadas con las variables
→ Φ𝑖 (𝑥1 , 𝑥2 … 𝑥𝑚 ) independientes o sus derivadas y NO aparecen en la
forma débil. Usualmente son fuerzas
𝑢 ≈ 𝑢𝑀𝐸𝐹 = ෍ 𝑢𝑖 N𝑖 Nodal interpolation (shape) functions
𝑈𝑥 (𝑥2 ) = c
→ N𝑖 (x1 , x2 … xm )
✓ Naturales (NBC)→ Relacionadas con las derivadas
de la variable dependiente y SI aparecen en la
𝑢 ≈ 𝑢𝑀𝐸𝐹 = ෍ 𝑐𝑖 Φ𝑖 + ෍ 𝑢𝑗 N𝑗 forma débil. Usualmente son desplazamientos

𝑑𝑈𝑥
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 𝑑𝑥 𝑑𝑥 𝑛

Hay regiones del Método de Galerkin

dominio donde el error
U Propuesto por el ingeniero Ruso Boris Galerkin en 1915. Asume las funciones
puede ser mayor
de peso igual a las funciones de prueba (o a las funciones de forma).
U analítica
𝑢 ≈ 𝑢𝐹𝐸𝑀 = ෍ 𝑐𝑖 Φ𝑖 𝑢 ≈ 𝑢𝐹𝐸𝑀 = ෍ 𝑢𝑖 N𝑖

𝜕𝑢𝐹𝐸𝑀 𝜕𝑢𝐹𝐸𝑀
𝑊𝑖 = 𝑊𝑖 =
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

1. Escribir la forma fuerte integral ponderada

Pistas
𝑥2
න 𝑊𝑖 (𝑥) 𝑅(𝑥, … ) 𝑑𝑥 = 0
𝑥1
2. Asumir una solución de prueba que satisfaga el orden de diferenciabilidad requerido → La cantidad de constantes independientes
3. Escribir la solución de prueba en forma de funciones de forma o funciones de forma en el paso 2 se reducirá en el paso 3
de manera que se cumplan todas las BC.
según la cantidad de BC
𝑢 ≈ 𝑢𝑀𝐸𝐹 = ෍ 𝑐𝑖 Φ𝑖 𝑢 ≈ 𝑢𝑀𝐸𝐹 = ෍ 𝑢𝑖 N𝑖
4. Asumir las funciones de peso igual a las funciones de prueba o de forma → El tamaño del sistema de ecuaciones que
se debe resolver será igual al número de
𝑊𝑖 = Φ𝑖 𝑊𝑖 = N𝑖 constantes que resultan del paso 3
5. Reemplazar las funciones de prueba (o de forma) y las funciones de peso en la forma
fuerte integral ponderada.
𝑥2 𝑥2
න Φ𝑖 𝑅 𝑑𝑥 = 0 න N𝑖 𝑅 𝑑𝑥 = 0
𝑥1 𝑥1

6. Resolver el sistema de ecuaciones para obtener el valor de los coeficientes (c𝑖 o u𝑖 )

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

Compare your FEM solution to the analytic solution

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

Enfoque físico Energía potencial total πP

Enfoque matemático
𝒙𝟐 En mecánica, la forma débil (funcional o
𝒅𝒖 𝒅𝒏 𝒖
𝑰 𝒖 = න 𝑭 𝒙, 𝒖, ,… 𝒏 variacional) corresponde a la energía
𝒙𝟏 𝒅𝒙 𝒅𝒙 potencial total del sistema

П𝑝 = 𝐸𝑝𝑠 − 𝑊𝑃
𝑰 𝒖 = 𝑩 𝒖, 𝒖 + 𝑳(𝒖) 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

Método de Rayleigh-Ritz→ Minimizar Procedimiento

el funcional respecto a los
coeficientes 1. Escribir la forma débil basada en el funcional de energía potencial
П𝑝 = 𝐸𝑝𝑠 − 𝑊𝑃

2. Asumir una solución de prueba que satisfaga el orden de diferenciabilidad requerido

3. Escribir la solución de prueba en forma de funciones de forma o funciones de forma de
manera que se cumplan todas las EBC.

𝑢 ≈ 𝑢𝑀𝐸𝐹 = ෍ 𝑐𝑖 Φ𝑖 𝑢 ≈ 𝑢𝑀𝐸𝐹 = ෍ 𝑢𝑖 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

John William Strutt, 3rd Walther Heinrich Wilhelm 𝜕П𝑝 𝜕П𝑝

Baron Rayleigh (1842-1919) Ritz (1878-1909) =0 =0
𝜕c𝑖 𝜕u𝑖
6. Resolver el sistema de ecuaciones para obtener el valor de los coeficientes (c𝑖 o u𝑖 )

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

Compare your FEM solution to the analytic solution

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

Funciones de forma (interpolación) y discretización

u(x) u(x)
𝑢(0) = 0 Lk, uk+1
𝑑 2 𝑈𝑥
𝐴𝐸 +𝑛=0 0≤𝑥 ≤𝐿 𝑑𝑢 𝑢(𝑥) = 𝑐1 𝑥 + 𝑐0
𝑑𝑥 2 𝐴𝐸 (𝐿) = 0
𝑑𝑥 0, Uk
𝑢𝑘+1 − 𝑢𝑘
𝑐1 =
𝐿𝑘
ux(x)
𝑐0 = 𝑢𝑘

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

Funciones de forma (interpolación) y discretización

u(x) u(x) u(x) Lk, uk+1

1 1
𝑥 𝑥
𝑁𝑘+1 = 0, Uk
)+ (
𝑁𝑘 = 1 −
𝑢𝑘
( 𝐿𝑘 𝑢𝑘+1 𝐿𝑘
) =

x x x
Lk Lk

𝑢𝑀𝐸𝐹 = 𝑢𝑘 𝑁𝑘 + 𝑢𝑘+1 𝑁𝑘+1 𝑥 𝑥

𝑢𝑀𝐸𝐹 = 𝑢𝑘 1− + 𝑢𝑘+1
𝐿𝑘 𝐿𝑘

¡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

Procedimiento de ensamble y solución de un sistema discretizado

Los elementos finitos conectados comparten grados de libertan, por lo tanto, sus resultados de desplazamiento en los puntos de
conexión deben ser iguales. Al procedimiento que permite considerar esto se le denomina ensamble
1. Definir los nodos y elementos, usualmente usando números
2. Escribir las matrices y vectores característicos de cada

(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)

LK 5. Eliminar las filas asociadas a los desplazamientos nulos y a las

fuerzas desconocidas. Esto dará un sistema de la siguiente
forma con incógnitas solo de desplazamientos. La matriz [K]
resultante se conoce como matriz reducida
N+1

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.

Compare your solutions to the analytic solution (https://es.symbolab.com/solver/ordinary-differential-equation-calculator/)

5 kN/m 10 kN/m

3m 1.5 [m] 1.5 [m]

Figure A. Bar subjected to uniform axial load

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 𝑀𝑗𝑍
𝐹𝐿 = 𝑇 𝐹𝐺

αij 𝑈𝑥𝑖𝐿 𝑈𝑥𝑖𝐺

𝑐 𝑠 0 0 0 0
XG 𝑈𝑦𝑖𝐿 −𝑠 𝑐 0 0 0 0 𝑈𝑦𝑖𝐺
𝜃𝑖𝐿 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

Matrices de rigidez de elementos finitos en coordenadas globales

Matriz de rigidez del elemento Matriz de rigidez del elemento barras en coordenadas globales
cercha en coordenadas globales

𝑐 𝑠 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

Formulaciones de elementos finitos de barras rotulados

A partir de los elementos anteriormente formulados, se pueden deducir elementos con rótulas, utilizando el procedimiento conocido
como condensación estática

Ejemplo para viga conectada con articulación en el extremo izquierdo

K0 K1
𝑉𝑖 12 6𝐿 −12 6𝐿 𝑈𝑦𝑖 𝑈𝑦𝑖
Reordenar 𝑉𝑖 12 −12 6𝐿 6𝐿 Agrupar
0 𝐸𝐼 6𝐿 4𝐿2 −6𝐿 2𝐿2 𝜃𝑖 𝑉𝑗 𝐸𝐼 −12 12 −6𝐿 −6𝐿 𝑈𝑦𝑗 𝐾 𝐾1 𝑈1
= 3 𝐹1
𝑉𝑗 𝐿 −12 −6𝐿 12 −6𝐿 𝑈𝑦𝑗 = 3 = 0
𝑀𝑖 𝐿 6𝐿 −6𝐿 4𝐿2 2𝐿2 𝜃𝑖 0 𝐾2 𝐾3 𝑈2
𝑀𝑗 6𝐿 2𝐿2 −6𝐿 4𝐿2 𝜃𝑗 0 6𝐿 −6𝐿 2𝐿2 4𝐿2 𝜃𝑗
K2 K3
Separar
Despejar 𝐾3 𝑈2 = − 𝐾2 𝑈1 Reemplazar 𝐹1 = 𝐾0 𝑈1 − 𝐾1 𝐾3 −1
𝐾2 𝑈1
operaciones 𝐹1 = 𝐾0 𝑈1 + 𝐾1 𝑈2
−1 −1
𝑈2 = − 𝐾3 𝐾2 𝑈1 𝐹1 = ( 𝐾0 − 𝐾1 𝐾3 𝐾2 ) 𝑈1
0 = 𝐾2 𝑈1 + 𝐾3 𝑈2
3𝐸𝐼 3𝐸𝐼 3𝐸𝐼
0 −
Resolver la matriz 𝑉𝑖 𝐿3 𝐿3 𝐿2 𝑈𝑦𝑖 ???
𝐸𝐼 12 −12 6𝐿 6𝐿 1 0 0 0 0
condensada −12 12 −6𝐿 − −6𝐿 [6𝐿 − 6𝐿 2𝐿2 ] 0
3𝐸𝐼 3𝐸𝐼 3𝐸𝐼
𝜃𝑖 ???
𝐿3 4𝐿 2
𝑉𝑗 = 𝑈𝑦𝑗 + ? ? ?
6𝐿 −6𝐿 4𝐿2 2𝐿2 − 3 0 − 2
𝑀𝑗 𝐿 𝐿3 𝐿
3𝐸𝐼 3𝐸𝐼 3𝐸𝐼 𝜃𝑗 ???
0 − 2
𝐿2 𝐿 𝐿

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

general fuerzas uniformes

YG
𝑁𝑖𝑗 = − න 𝑛𝑖𝑗 𝑥 𝑑𝑥 − 𝑁𝑖 𝑁𝑖𝑗 = −𝑛𝑖𝑗 𝑥 − 𝑁𝑖
αij
𝑉𝑖𝑗 = − න 𝑤𝑖𝑗 𝑥 𝑑𝑥 + 𝑉𝑖 𝑉𝑖𝑗 = −𝑤𝑖𝑗 𝑥 + 𝑉𝑖
XG
𝑤𝑖𝑗 𝑥 2
𝑀𝑖𝑗 = − + 𝑉𝑖 𝑥 − 𝑀𝑖
𝑀𝑖𝑗 = න 𝑉𝑖𝑗 𝑥 𝑑𝑥 − 𝑀𝑖 2

𝐹𝐿 = 𝐾𝐿 𝑈𝐿 + 𝐹𝐿𝑒 Computacionalmente, es usual evaluar las fuerzas internas en

𝑁𝑖
puntos “p” específicos separados una distancia uniforme “s”.
𝑉𝑖
Estos puntos se conocen como estaciones de salida y permiten
𝑀𝑖 aproximar gráficamente los diagramas de fuerzas internas.
𝑈𝐿 = 𝑇 𝑈𝐺 𝐹𝐿 = 𝑁
𝑗
𝑉𝑗 𝐿
𝑠=
𝑀𝑗 𝑝−1

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:

a) W=50 kN/m, P=100 kN

b) W=0 kN/m, P=100 kN
w
c) W=50 kN/m, P=0 kN
d) W=0 kN/m, P varying from 0 to 100 kN in steps of 20 kN
P

5m
3m

Figure A. Simple frame

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)

Elementos finitos típicos Elementos estructurales

Figure 19. Shell element típicos
✓ Triángulo
✓ Rectángulo ✓ Losa maciza
El modelo matemático por resolver
✓ Muros de contención
es el de la mecánica del medio
✓ Zapatas
continuo

21
 Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes

Comportamiento en el plano (membrana)

Esfuerzos Deformaciones
dA
σyy

σyx
Configuración
deformada σxy ϵyy

σxx σxx

σxy
fv ϵxx
Configuración σyx
Y F
de referencia
σyy

Sólido con espesor

fs X constante t
σxy= σyx ϒxy

22
 Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes

Comportamiento en el plano (membrana)

Variables

✓ Esfuerzo → [F/ L2] ✓ Fuerza volumétrica → [F/L3]

σ𝑥𝑥 𝑓𝑣𝑥
{σ(𝑥, 𝑦)} = σ𝑦𝑦 {𝑓𝑣 (𝑥, 𝑦)} =
𝑓𝑣𝑦
σ𝑥𝑦

✓ Deformación → [L/L] or [rad] ✓ Fuerza de superficie→ [F/L2]

ϵ𝑥𝑥 𝑓𝑠𝑥
fv {ϵ(𝑥, 𝑦)} = ϵ𝑦𝑦 {𝑓𝑠 (𝑥, 𝑦)} =
𝑓𝑠𝑦
F γ𝑥𝑦
Y
✓ Desplazamiento [L] ✓ Fuerza puntual→ [F]

fs 𝑈𝑥 𝐹𝑥
{𝑈(𝑥, 𝑦)} = 𝑈 {𝐹} = 𝐹
X 𝑦 𝑦

23
 Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes

In plane (membrane) behavior of area element

Ecuaciones de equilibrio Ecuaciones constitutivas Ecuaciones cinemáticas

𝜕𝑈𝑥 𝜕𝑈𝑥 𝜕𝑈𝑦

𝜕σ𝑥𝑥 𝜕σ𝑥𝑦 𝜎𝑥𝑥 − ν 𝜎𝑦𝑦 𝜎𝑦𝑦 − ν 𝜎𝑥𝑥 𝜖𝑥𝑥 = 𝛾𝑥𝑦 = +
+ + 𝑓𝑣𝑥 = 0 τ𝑥𝑦 = τ𝑦𝑥 𝜖𝑥𝑥 = 𝜖𝑦𝑦 = 𝜕𝑥 𝜕𝑦 𝜕𝑥
𝜕𝑥 𝜕𝑦 𝐸 𝐸

𝜕σ𝑦𝑥 𝜕σ𝑦𝑦 𝜎𝑥𝑦 2(1 + ν)𝜎𝑥𝑦 𝜕𝑈𝑦

+ + 𝑓𝑣𝑦 = 0 𝛾𝑥𝑦 = = 𝜖𝑦𝑦 =
𝜕𝑥 𝜕𝑦 𝐺 𝐸 𝜕𝑦
𝜕
0
𝜕 𝜕 1 ν 0 ϵ𝑥𝑥 𝜕𝑥
0 σ𝑥𝑥 σ𝑥𝑥 ϵ𝑥𝑥 𝜕 𝑈𝑥
𝜕𝑥 𝜕𝑦 σ 𝑓𝑣𝑥 0 𝐸 ν 1 0 ϵ𝑦𝑦 = 0
𝑦𝑦 + = σ𝑦𝑦 = ϵ
𝜕 𝜕 σ 𝑓𝑣𝑦 0 1 − ν2 1 − ν 𝑦𝑦 γ𝑥𝑦 𝜕𝑦 𝑈𝑦
0 𝑥𝑦 𝜎𝑥𝑦 0 0 γ𝑥𝑦
𝜕𝑦 𝜕𝑥 2 𝜕 𝜕
𝜕𝑦 𝜕𝑥

𝐿 𝑇 σ + 𝑓𝑣 = 0 ϵ = [𝐿] 𝑈
σ = [𝐷] ϵ

24
 Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes

Formulación variacional general de los elementos de área

→ Ahora, recordemos que los
П𝑝 = 𝐸𝑝𝑠 − 𝑊𝑣 − 𝑊𝑠 − 𝑊𝑝 desplazamientos se asumen como

ϵ = [𝐿] 𝑈 {𝑈} = [N𝑖 ]{𝑈𝑖 }

1 𝑉 𝑇
𝑉
𝑇
𝐴
𝑇
П𝑝 = න ϵ σ 𝑑𝑉 − න 𝑈 𝑓𝑣 𝑑𝑉 − න 𝑈 𝑓𝑠 𝑑𝐴 − ෍ 𝑈𝑖 𝑇 {𝐹}
2 0 0 0
ϵ = [𝐿][N𝑖 ]{𝑈𝑖 }
𝑉
𝐸ϵ 2 → Para entender la expresión de arriba se puede hacer
П𝑝 = න − 𝐹𝑈 analogía con el caso 1D estudiado al inicio del curso → Reemplazando el funcional
0 2
П𝑝 =
𝐴 𝐴 𝐿
𝑡 𝑡 𝐴
П𝑝 = න ϵ 𝑇 [𝐷] ϵ 𝑑𝐴 − 𝑡 න 𝑈 𝑇
𝑓𝑣 𝑑𝐴 − 𝑡 න 𝑈 𝑇 𝑇
𝑓𝑠 𝑑𝐿 − ෍ 𝑈𝑖 {𝐹} න 𝑈 𝑇
𝑁𝑖 𝑇
𝐿 𝑇
𝐷 [𝐿][𝑁𝑖 ] {𝑈𝑖 }𝑑𝐴
2 0 0 0 2 0 𝑖
𝐴
Energía de deformación Potencial de trabajo de − 𝑡 න 𝑈𝑖 𝑇
𝑁𝑖 𝑇
𝑓𝑣 𝑑𝐴
las fuerzas volumétricas, 0
𝐿
de superficie y puntuales 𝑇 𝑇
−𝑡 න 𝑈𝑖 𝑁𝑖 𝑓𝑠 𝑑𝐿
0
→ Nótese que en el caso de membranas esto corresponde a un sistema 2x2
𝑇
donde las incógnitas son las funciones de desplazamientos en X y Y − ෍ 𝑈𝑖 {𝐹}

25
 Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes

Formulación variacional general de los elementos de área

→ Aplicando Rayleigh-Ritz → Para el término de energía de → Para el término de las fuerzas puntuales
deformación
𝜕П𝑝 𝜕
=0= 𝑇
𝜕𝑈𝑖 ෍ 𝑈𝑖 𝐹
𝜕 𝑡 𝐴 𝑇 𝑇 𝑇
𝜕𝑈𝑖
𝜕 𝑡 𝐴 න 𝑈 𝑁𝑖 𝐿 𝐷 [𝐿][𝑁𝑖 ] {𝑈𝑖 }𝑑𝐴
න 𝑈𝑖 𝑇 𝑁𝑖 𝑇 𝐿 𝐷 [𝐿][𝑁𝑖 ] {𝑈𝑖 }𝑑𝐴 𝜕𝑈𝑖 2 0 𝑖 → Para entender esta
𝜕𝑈𝑖 2 0 𝜕
𝐴 𝐹𝑈 = 𝐹 derivada se puede
𝜕 𝑡 𝐴 𝜕 𝜕𝑈𝑖
− 𝑡 න 𝑈𝑖 𝑇 𝑁𝑖 𝑇 𝑓𝑣 𝑑𝐴 න 𝑈 𝑇
𝑁𝑖 𝑇
𝐿 𝑇
𝐷 [𝐿][𝑁𝑖 ] {𝑈𝑖 }𝑑𝐴
hacer analogía con el
𝜕𝑈𝑖 0 2 0 𝜕𝑈𝑖 𝑖 caso 1D del resorte
𝐿
𝜕
− 𝑡 න 𝑈𝑖 𝑇 𝑁𝑖 𝑇 𝑓𝑠 𝑑𝐿 → Para entender esta 𝜕
𝜕𝑈𝑖 0 𝜕 𝐾𝑈 2 𝑇
= 𝐾𝑈 derivada se puede ෍ 𝑈𝑖 𝐹 = {𝐹}
𝜕 𝜕𝑈𝑖 2 𝜕𝑈𝑖
− ෍ 𝑈𝑖 𝑇 {𝐹} hacer analogía con el
𝜕𝑈𝑖 caso 1D del resorte
𝐴
→ La solución de los términos de fuerzas de
𝜕 superficie y volumétricas corresponden a las
𝑡න 𝑁𝑖 𝑇 [𝐷][𝑁𝑖 ]{𝑈𝑖 }𝑑𝐴
𝐹 = 𝐾 𝑈 + {𝐹𝑒 } 0 𝜕𝑈𝑖
fuerzas de empotramiento (ver Libro de
Seshu)
𝐴 𝐴
La clave en los elementos de área es 𝑇
determinar la matriz cinemática 𝑡 න 𝑁𝑖 𝐿 [𝐷][𝐿][𝑁𝑖 ] 𝑈𝑖 𝑑𝐴 = 𝑡 න 𝐵 𝑇 [𝐷] [𝐵] 𝑑𝐴 𝑈𝑖 = 𝐾 𝑈
𝑇
0 0

26
 Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes

Elemento triangular de deformación constante (CST)

Geometría del CST Uy3
3 x3, y3
Ux3
Elemento finito

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

Fuerzas externas en el CST

Fuerzas volumétricas Fuerzas de superficie Fuerzas nodales
Fy3

Fx3

Fy2

Fy1
Fx2
fs
fv
Fx1
𝑓𝑣𝑥 𝑓𝑠𝑥 𝐹𝑥
{𝑓𝑣 (𝑥, 𝑦)} = {𝑓𝑠 (𝑥, 𝑦)} = {𝑓𝑝 (𝑥, 𝑦)} = 𝐹
𝑓𝑣𝑦 𝑓𝑠𝑦 𝑦

28
 Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes

Funciones de forma del CST

La deformación se asume constante, por ende, los desplazamientos se asumen lineales 𝑈𝑥1
𝑈𝑦1
𝑈𝑥 = 𝑈𝑥1 𝑁1 + 𝑈𝑥2 𝑁2 + 𝑈𝑥3 𝑁3 𝑈𝑦 = 𝑈𝑦1 𝑁1 + 𝑈𝑦2 𝑁2 + 𝑈𝑦3 𝑁3 𝑈𝑥 𝑁1 0 𝑁2 0 𝑁3 0 𝑈𝑥2 {𝑈} = [N𝑖 ]{𝑈𝑖 }
𝑈𝑦 = 𝑈𝑦2
0 𝑁1 0 𝑁2 0 𝑁3
𝑈𝑥3
Z 𝑈𝑦3
1
Y
𝑥2 𝑦3 − 𝑥3 𝑦2 + y2 − y3 x + x3 − x2 y
Z 𝑁1 =
1 det[𝐽]
Y
X 𝑥3 𝑦1 − 𝑥1 𝑦3 + y3 − y1 x + x1 − x3 y
𝑁2 =
det[𝐽]
Z X 𝑥1 𝑦2 − 𝑥2 𝑦1 + y1 − y2 x + x2 − x1 y
1 𝑁3 =
Y det[𝐽]
𝑥1 − 𝑥3 𝑦1 − 𝑦3 𝑥13 𝑦13 𝜃𝑚𝑎𝑥 − 60° 60° − 𝜃𝑚𝑖𝑛
𝐽 = 𝑥 −𝑥 = 𝑠𝑘𝑒𝑤𝑛𝑒𝑠𝑠 = max ,
2 3 𝑦2 − 𝑦3 𝑥23 𝑦23 120° 60°
X El Jacobiano de transformación condiciona la exactitud de la solución numérica. Ojo con los elementos
muy oblicuos (0< skewness <0.5)

29
 Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes

Deducción de la matriz cinemática del CST

ϵ = [𝐿][N𝑖 ]{𝑈𝑖 }

𝜕
𝑈𝑥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

Ejemplo para el primer término ϵ = [𝐵] 𝑈𝑖

considerando la notación
xi -xj=xij=-xji
yi -yj=yij=-yji Matriz cinemática para el CST

𝜕𝑁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

Deducción de la matriz de rigidez para el CST

𝐴
𝐾 = 𝑡 න 𝐵 𝑇 [𝐷] [𝐵] 𝑑𝐴
0

𝐾 = 𝐴𝑡 𝐵 𝑇 𝐷 𝐵

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

Formulación general del CST sin incluir fuerzas de empotramiento

Uy3

𝐹 = 𝐾 𝑈 + {𝐹𝑒 }

𝐹𝑥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

Desventajas del elemento CST

Esfuerzos parásitos (Shear locking) Necesitdad de alta discretización

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

Out of plane behavior of area element

z
Stresses

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

Out of plane behavior of area element

Normal stresses
Bending around Y Bending around X
dy
𝑚𝑥 = න 𝜎𝑥 𝑧𝑑𝑧 𝑚𝑦 = න 𝜎𝑦 𝑧𝑑𝑧
z dx

σy If linear function is assumed for the stresses

t/2 y my 12 12
𝜎𝑥 = 𝑚𝑥 𝑧 3 𝜎𝑦 = 𝑚𝑦 𝑧 3
t/2 𝑡 𝑡
x
σx mx

Bending moments mx and my are

moments per unit length

40
 Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes

Out of plane behavior of area element

Shear stresses due to torsion Torsion

dy 𝜏𝑥𝑦 = 𝜏𝑦𝑥

z dx 𝑚𝑥𝑦 = න 𝜏𝑥𝑦 𝑧𝑑𝑧

myx If linear function is assumed for the stresses

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

Out of plane behavior of area element

Shear stresses due to torsion Shear in face x Shear in face y

dy
𝑣𝑥 = න 𝜏𝑥𝑧 𝑑𝑧 𝑣𝑦 = න 𝜏𝑦𝑧 𝑑𝑧
z dx

If parabolic stress is assumed

τyz 12𝑣𝑥 ‫𝑧𝑑𝑧 ׬‬ 12𝑣𝑦 ‫𝑧𝑑𝑧 ׬‬

t/2 y 𝜏𝑥𝑧 =− 𝜏𝑦𝑧 =−
𝑡3 𝑡3
t/2 vy
x vx τxz Shear forces vx and vy is a force per
σx
unit length

42
 Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes

Out of plane behavior of area element

Equilibrium equations

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

Out of plane behavior of area element

Deriving the first equilibrium equation respect X and the second equilibrium respect to Y and replace the result in the third equilibrium equation

𝜕 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

Out of plane behavior of area element

Kinematic equations
𝜕 2 𝑈𝑧
𝜕𝑈𝑥 𝜖𝑥 = − 𝜕𝑧
𝜕𝑥 2
𝜖𝑥 = ¿Geometric interpretation?
𝜕𝑥
𝜕 2 𝑈𝑧
𝜕𝑈𝑦 𝜖𝑥 = − 𝜕𝑧
𝜖𝑦 = 𝜕𝑦 2
𝜕𝑦
𝜕𝑈𝑥 𝜕𝑈𝑦 𝜕 2 𝑈𝑧
𝛾𝑥𝑦 = + 𝛾𝑥𝑦 = −2 𝜕𝑧
𝜕𝑦 𝜕𝑥 𝜕𝑥𝜕𝑦

𝜕𝑈𝑦 𝜕𝑈𝑧 𝜕𝑈𝑧 𝜕𝑈𝑦 𝜕𝑈𝑧

𝛾𝑦𝑧 = + =− 𝜕𝑈𝑦 = − 𝜕𝑧
𝜕𝑧 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑦
𝜕𝑈𝑧 𝜕𝑈𝑥 𝜕𝑈𝑧 𝜕𝑈𝑥 𝜕𝑈𝑧
𝛾𝑧𝑥 = + =− 𝜕𝑈𝑥 = − 𝜕𝑧
𝜕𝑥 𝜕𝑧 𝜕𝑥 𝜕𝑧 𝜕𝑥

45
 Método de elementos finitos-Análisis avanzado de estructuras
Prof. David Cotes

Out of plane behavior of area element

Replacing kinematic equations in constitutive equation Deriving respect to x and y, and replacing in the equilibrium equation
𝜕 2 𝑈𝑧 𝜕 2 𝑚𝑥 𝜕 2 𝑚𝑥𝑦 𝜕 2 𝑚𝑦
− 𝜕𝑧 +2 + + 𝑞𝑧 = 0
1 𝜈 0 𝜕𝑥 2 𝜕𝑥 2 𝜕𝑥𝜕𝑦 𝜕𝑦 2
𝜎𝑥
𝐸 𝜈 1 0 𝜕 2 𝑈𝑧
𝜎𝑦 = − 𝜕𝑧
1 − 𝜈2 1−𝜈 𝜕𝑦 2
𝜏𝑥𝑦 0 0
2 𝜕 2 𝑈𝑧
−2 𝜕𝑧 Differential equation of the Kirchhoff-Love plate
𝜕𝑥𝜕𝑦
𝜕 4 𝑈𝑧 𝜕 4 𝑈𝑧 𝜕 4 𝑈𝑧 12 1 + 𝑣 2
Multiplying by z and integrating respect to dz +2 2 2+ − 𝑞𝑧 = 0
𝜕𝑥 4 𝜕𝑥 𝜕𝑦 𝜕𝑦 4 𝐸𝑡 3
𝜕 2 𝑈𝑧 ✓ Second order partial differential equation
−
1 𝜈 0 𝜕𝑥 2 ✓ Neglects the strain due to shear force
𝑚𝑥 ✓ Recommended for plates where L/t>10
𝐸𝑡 3 𝜈 1 0 𝜕 2 𝑈𝑧
𝑚𝑦 = −
12 1 − 𝜈 2 1−𝜈 𝜕𝑦 2
𝑚𝑥𝑦 0 0
2 𝜕 2 𝑈𝑧
−2
𝜕𝑥𝜕𝑦

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

z Rectangular plate finite element

Trial solution
𝑈𝑧1 , 𝜃𝑥1 , 𝜃𝑦1 𝑈𝑧4 , 𝜃𝑥4 , 𝜃𝑦4 𝑈𝑧 = 𝑐1 + 𝑐2 𝑥 + 𝑐3 𝑦 + 𝑐4 𝑥 2 + 𝑐5 𝑥𝑦 + 𝑐6 𝑦 2 + 𝑐7 𝑥 3 + 𝑐8 𝑥 2 𝑦 + 𝑐9 𝑥𝑦 2
y +𝑐10 𝑦 3 + 𝑐11 𝑥 3 𝑦 + 𝑐12 𝑥𝑦 3

x Stiffness matrix
2a

2b 𝑈𝑧3 , 𝜃𝑥3 , 𝜃𝑦3

𝑈𝑧2 , 𝜃𝑥2 , 𝜃𝑦2

Sources:
1. Structural Analysis with the Finite Element Method - Eugenio Oñate
47
 Método de elementos finitos
Prof. David Cotes

THANKS FOR YOUR

ATTENTION…
ANY QUESTIONS?