GUÍA DIDÁCTICA DE CIV 313

Introducción a la aplicación del Método de Elementos Finitos JCB ©2020

U.A.T.F. Ingeniería Civil Ing. Juan Carlos Barrios C.

Ejemplo 1. – Resolver la siguiente estructura sometida a fuerzas axiales. Asumir EA = Constante.

La discretización del modelo será:

JCB ©2020
Cálculo de las rigideces, fuerzas equivalentes y fuerzas nodales para cada elemento.

Elemento 1

𝐸𝐴 (1) 1 0.7 ∗ 2 1 (1)

𝑓1
−1 〈𝐹𝑒〉(1) = { }
𝑘 (1) = ( ) [ ] 2 1 {𝑓}(1) = { }
2 −1 1 (1)

Introducción a la aplicación del Método de Elementos Finitos

𝑓2
Elemento 2

𝐸𝐴 (2) 1 0.7 ∗ 2 1 (2)

𝑓1
−1 〈𝐹𝑒〉(2) = { }
𝑘 (2) = ( ) [ ] 2 1 {𝑓}(2) = { }
2 −1 1 (2)
𝑓2
Elemento 3

𝐸𝐴 (3) 1 0.7 ∗ 2.7 1 (3)

𝑓1
−1 〈𝐹𝑒〉(3) = { }
𝑘 (3) = ( ) [ ] 2 1 {𝑓}(3) = { }
2.7 −1 1 (3)
𝑓2
Elemento 4

𝐸𝐴 (4) 1 0.7 ∗ 1.6 1 (4)

𝑓1
−1 〈𝐹𝑒〉(4) = { }
𝑘 (4) = ( ) [ ] 2 1 {𝑓}(4) = { }
1.6 −1 1 (4)
𝑓2

Tomando en cuenta la condición de equilibrio, tenemos:

(1)
Nudo 1: 𝑓1 =𝑅
(1) (2)
Nudo 2: 𝑓2 + 𝑓1 = 1.20
(2) (3)
Nudo 3: 𝑓2 + 𝑓1 = 1.50
(3) (4)
Nudo 4: 𝑓2 + 𝑓1 = 1.00
(4)
Nudo 5: 𝑓2 = 2.40

1
 GUÍA DIDÁCTICA DE CIV 313
Introducción a la aplicación del Método de Elementos Finitos JCB ©2020
U.A.T.F. Ingeniería Civil Ing. Juan Carlos Barrios C.

Por lo que la ecuación será: [𝐾]{𝑎} − {𝐹𝑒} = {𝑓}

0.5000 −0.500 0.000 0.000 0.000 𝑢1 0.700 𝑅

−0.500 1.000 −0.500 0.000 0.000 𝑢2 1.400 1.2
𝐸𝐴 0.000 −0.500 0.8704 0.3704 0.000 𝑢3 − 1.605 = 1.5
0.000 0.000 −0.3704 0.9954 −0.625 𝑢4 1.505 1.0
[ 0.000 0.000 0 −0.625 0.625 ] {𝑢5} {0.560} {2.4}
Por las condiciones de borde, tenemos restringido el nudo 1, es decir: u1 =0

Esto hace que se elimine la fila y columna 1, obteniéndose un sistema reducido, el cual al
resolverse obtenemos:
𝑢2 22.4200
𝑢3 1 39.6400
{ } { }
𝑢4 𝐸𝐴 54.3955
𝑢5 59.1315
Estos resultados podrán permitir dibujar el estado deformado de la estructura, bajo el sistema
de cargas.

Las deformaciones unitarias, se calcular para cada elemento finito, con la siguiente relación:

JCB ©2020
[𝜀] = [𝐵]{𝑎} Donde: {𝐵} = [
−1 1
]
𝐿 𝐿
Entonces:

Elemento 1

Introducción a la aplicación del Método de Elementos Finitos

1 1 1 0 11.21
[𝜀] = [− ] { }=
2 2 𝐸𝐴 22.42 𝐸𝐴

Elemento 2
1 1 1 22.42 8.61
[𝜀] = [− ] { }=
2 2 𝐸𝐴 39.64 𝐸𝐴
Elemento 3
1 1 1 39.64 5.465
[𝜀] = [− ] { }=
2.7 2.7 𝐸𝐴 54.3955 𝐸𝐴
Elemento 4
1 1 1 54.3955 2.96
[𝜀] = [− ] { }=
1.6 1.6 𝐸𝐴 59.1315 𝐸𝐴

El tensor de los esfuerzos, se calcula según la expresión siguiente:

[𝜎] = [𝐸][𝜀]

Por tanto, los resultados serán:

Elemento 1
11.21
𝜎=
𝐴

Elemento 2
8.61
𝜎=
𝐴

2
 GUÍA DIDÁCTICA DE CIV 313
Introducción a la aplicación del Método de Elementos Finitos JCB ©2020
U.A.T.F. Ingeniería Civil Ing. Juan Carlos Barrios C.

Elemento 3
5.465
𝜎=
𝐴

Elemento 4
2.96
𝜎=
𝐴

Ejemplo 2. – Considere la barra de la figura siguiente. Determinar los desplazamientos nodales,

los esfuerzos en los elementos y las reacciones en los soportes.

JCB ©2020
Cálculo de las matrices de rigidez de cada elemento

Elemento 1

Introducción a la aplicación del Método de Elementos Finitos

Elemento 2

3
 GUÍA DIDÁCTICA DE CIV 313
Introducción a la aplicación del Método de Elementos Finitos JCB ©2020
U.A.T.F. Ingeniería Civil Ing. Juan Carlos Barrios C.

Elemento 3

Ensamblando la matriz de rigidez general, obtenemos:

JCB ©2020
Introducción a la aplicación del Método de Elementos Finitos
Los vectores de fuerza y desplazamiento nodal, están dados por:
𝑅1 0
300000 𝑢2
{𝑎} = { }
{𝐹} = { } 𝑢3
0
𝑅4 0

Las condiciones de borde, establecen que los desplazamientos en los nodos 1 y 4 son
nulos, por tanto, permite eliminar las filas y columnas 1 y 4, por lo que el sistema reducido
a resolver, proporciona los siguientes resultados:
81
𝑢2
130
{𝑎} = { } =
9
𝑢3
{ 26 }
Las reacciones podrán calcularse a partir de la siguiente relación:

−1000000 81 9 −2700000
𝑅1 = ∗ +0∗ = 𝑁
3 130 26 13

81 800000 9 −1200000
𝑅4 = 0 ∗ − ∗ = 𝑁
130 3 26 13

4
 GUÍA DIDÁCTICA DE CIV 313
Introducción a la aplicación del Método de Elementos Finitos JCB ©2020
U.A.T.F. Ingeniería Civil Ing. Juan Carlos Barrios C.

Los esfuerzos nodales están dados por:

[𝜎] = [𝐸][𝐵]{𝑎}

Elemento 1
0
200000 10800 𝑁
[𝜎] = [−1 1] { 81 } =
150 13 𝑚𝑚2
130
Elemento 2
81
200000 130 −4800 𝑁
[𝜎] = [−1 1] =
150 9 13 𝑚𝑚2
{ 26 }
Elemento 3

JCB ©2020
9
200000 −3000 𝑁
[𝜎] = [−1 1] {26} =
300 13 𝑚𝑚2
0

Introducción a la aplicación del Método de Elementos Finitos

Ejemplo 3. – Determinar la matriz de la siguiente estructura:

• Determinación de la geometría de la sección del elemento:

La variación de la altura de la sección es: ℎ(𝑥) = 𝑎𝑜 + 𝑎1 𝑥 + 𝑎2 𝑥 2

Por tanto, para: x=0 2ℎ = 𝑎0

x=L ℎ = 𝑎0 + 𝐿𝑎1 + 𝐿2 𝑎2

x = 2L 2ℎ = 𝑎0 + 2𝐿𝑎1 + 4𝐿2 𝑎2

resolviendo el sistema, obtenemos:

5
 GUÍA DIDÁCTICA DE CIV 313
Introducción a la aplicación del Método de Elementos Finitos JCB ©2020
U.A.T.F. Ingeniería Civil Ing. Juan Carlos Barrios C.

𝑎0 = 2ℎ −2ℎ ℎ
𝑎1 = 𝑎2 =
𝐿 𝐿
Por tanto:
2ℎ ℎ
ℎ(𝑥) = 2ℎ − 𝑥 + 2 𝑥2
𝐿 𝐿
• Determinación de la rigidez. – Que esta dado por la siguiente relación.

𝑥2
Donde: 1
(𝑒) 𝑇 𝐵 = [−1 1]
𝑘 = ∫ 𝐵 𝐷𝐵𝐴𝑑𝑥 𝐿
𝑥1
Por tanto:
𝑥2
1 −1 1
𝑘 (𝑒) = ∫ [ ] 𝐸 ∗ [−1 1] ∗ 2 ∗ 𝑏 ∗ ℎ(𝑥)𝑑𝑥 =
𝐿 𝐿
𝑥1 1

𝑥2
2𝐸𝑏ℎ 1 −1 2𝑥 𝑥 2

JCB ©2020
(𝑒)
𝑘 = 2 [ ] ∫ (2 − + 2 ) 𝑑𝑥 =
𝐿 −1 1 𝑥1 𝐿 𝐿

2𝐸𝑏ℎ 1 −1 4
𝑘 (𝑒) = 𝐿2
[ ]3𝐿 =

Introducción a la aplicación del Método de Elementos Finitos

−1 1
Finalmente:

8𝐸𝑏ℎ 1 −1
𝑘 (𝑒) = 3𝐿
[ ]
−1 1

TALLER DE PRÁCTICA NO 1

Comprende en la resolución de un conjunto de ejercicios que deberán enviar al

siguiente correo, durante UNA semana a partir de la presente fecha 17 Marzo de 2020
en forma digital (Formato PDF) al siguiente correo

[email protected]
El estudiante podrá enviar su práctica, cualquier momento hasta fecha 23 Marzo de 020
hasta Hrs 14:00

Incluir como identificador el carimbo que se encuentra en el link del presente grupo.

6
 GUÍA DIDÁCTICA DE CIV 313
Introducción a la aplicación del Método de Elementos Finitos JCB ©2020
U.A.T.F. Ingeniería Civil Ing. Juan Carlos Barrios C.

Para las siguientes barrar que se muestran en las figuras P3-4 a P3-8 que se muestra,
determinar, el estado de deformación, las reacciones y tensores de deformación
unitaria y esfuerzos.

JCB ©2020
Introducción a la aplicación del Método de Elementos Finitos
Resuelva los ejercicios 3.5 y 3.11 del libro de Tirupathi R. Chandrupatla & Belegundu,
Introducción al Método de Elementos Finitos en la Ingeniería, 2da Edición.

Potosí, 17 de Marzo de 2020