Ejemplo 4.
1 (Página 104)
Introducción al estudio del Elemento Finito en Ingeniería
Segunda Edición
TITUPATHI R. CHANDRUPATLA – ASHOK D. BELEGUNDU
PRENTICE HALL
Considere la armadura de cuatro barras mostrada en la figura E4.1a. Para todos los
elementos E=29.5E6 psi y Ae=1in2
- determine la matriz de rigidez elemental para cada elemento
- ensamble la matriz de rigidez estructural K para toda la armadura
- Encuentre los desplazamientos nodales
- Calcule las fuerzas de reacción
- Recupere los esfuerzos para cada elemento.
Solución:
1) Convenciones.
Inicialmente se opta un sistema de coordenada cartesiana ortogonal y el sistema de
orientación local para cada elemento, definido por su nudo inicial y final como se
muestra.
1 de 21
Luego se divide la estructura en una serie de elementos finitos considerando sus puntos
extremos como nudos, se enumera nudos y barras. Cada miembro de la estructura debe
estar completamente definido por su nudo inicial y final dentro del sistema.
Según la numeración de los nudos, los grados de libertad son.
Tomando estas convenciones, se ordena los argumentos.
2 de 21
2. Argumentos
MEF/ARMADURAS/DATOS
2.1 Nudos
Cada fila representa un punto y las columnas son:
Columna 1: coordenada "x" global del nudo
Columna 2: corrdenada "y" global del nudo
NODE :=
1 2
1 0 0
2 40 0
3 40 30
4 0 30
2.2 Propiedades de los elementos
las propiedqades de la sección trasversal de los elementos. cada fila representa una
propiedad distinta y las columnas son:
Columna 1: Área de la sección trasversal del elemento
Columna 2: Módulo de elasticidad del material
PROP :=
1 2
1 1 2.95·107
2.3 Elementos
Cada fila representa una barra, contiene la información de la conectividad del elemento en el sistem
cada columna representa:
Columna 1: nudo inicial del elemento
Columna 2: nudo final del elemento
Columna 3: número de propiedad del elemento
MEMB :=
1 2 3
1 1 2 1
2 3 2 1
3 1 3 1
4 4 3 1
2.4 Restricciones/Apoyos
Cara fila representa un apoyo de la estructura, las columnas informan el comportamiento para cada
grado de libertad, la convensión es:
• "1" para los grados de libertad de desplazamiento restringido.
• "0" para los grados de libertad donde exista desplazamiento libre.
Cada columna representa:
3 de 21
MEF/ANÁLISIS DE ARMADURAS
Columna 1: número del nudo donde existe el apoyo
Columna 2: "ux?" información del desplazamiento en la dirección "x" global
Columna 3: "uy?" información del desplazamiento en la dirección "y" global
SUPP :=
1 2 3
1 1 1 1
2 2 0 1
3 4 1 1
2.5 cargas
Se admite cargas en los nudos y cada columna representa:
Columna 1: número del nudo en que actúa la carga
Columna 2: carga en la dirección "x" global
Columna 3: carga en la dirección "y" global
se debe ingresar considerando la orientación global del sistema
NLF :=
1 2 3
1 2 2·104 0
2 3 0 -2.5·104
MEF/ARMADURAS/DATOS
4 de 21
MEF/ANÁLISIS DE ARMADURAS
3. Formando Vector de Cargas
MEF/ARMADURAS/DATOS
MEF/ARMADURAS/VECTOR DE CARGAS
Ordenado las cargas según los grados de libertad.
F := f (x , y) ← 0
F ← matrix( 2 ⋅ rows( NODE) , 1 , f )
for i ∈ 1 .. rows( NLF )
n ← NLF
i, 1
fx ← NLF
i, 2
fy ← NLF
i, 3
F ← fx
2⋅ n− 1
F ← fy
2n
F
⎛ 0 ⎞
⎜ 0
⎟
⎜ ⎟
⎜ 2 × 104 ⎟
⎜ ⎟
F= ⎜ 0 ⎟
⎜ 0 ⎟
⎜ 4
⎟
⎜ −2.5 × 10 ⎟
⎜ 0 ⎟
⎜ ⎟
⎝ 0 ⎠
MEF/ARMADURAS/VECTOR DE CARGAS
5 de 21
MEF/ANÁLISIS DE ARMADURAS
4. Matriz de Rigidez de Elementos
MEF/ARMADURAS/DATOS
MEF/ARMADURAS/VECTOR DE CARGAS
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
4.1 Matriz de rigidez para el elemento
Elm := 1
4.1.1 la longitud del elemento y los
cosenos directores son
⎛ Le ⎞
⎜ λxe ⎟
⎜ ⎟
⎜ λye ⎟ := ni ← MEMBElm , 1
⎜ Ae ⎟
⎜ ⎟ nf ← MEMB
Elm , 2
⎝ Ee ⎠
prop ← MEMB
Elm , 3
xi ← NODE
ni , 1
yi ← NODE
ni , 2
xf ← NODE Con el programa se obtiene para el elemento
nf , 1
Elm = 1
yf ← NODE
nf , 2 • Longitud del elemento
Le ←
2
( xf − xi) + ( yf − yi)
2 Le = 40
xf − xi • Cosenos directores
λxe ←
Le
λxe = 1
yf − yi
λye ←
Le λye = 0
Ae ← PROP
prop , 1 • Propiedades
Ee ← PROP Ae = 1
prop , 2
T 7
( Le λxe λye Ae Ee ) Ee = 2.95 × 10
4.2.2 La matriz de rigidez en coordenada local está dado por:
Ee⋅ Ae ⎛ 1 −1 ⎞
k := ⋅⎜ ⎟
Le ⎝ −1 1 ⎠
29500000⋅ 1 ⎛ 1 −1 ⎞
k = ⋅⎜ ⎟
40 ⎝ −1 1 ⎠
⎛ 7.375 × 105 −7.375 × 105 ⎞
k=⎜ ⎟
⎜ 5 5 ⎟
⎝ −7.375 × 10 7.375 × 10 ⎠
6 de 21
MEF/ANÁLISIS DE ARMADURAS
4.2.3 La matriz de transformación de desplazamientos está dado por:
⎛ λxe λye 0 0 ⎞
L := ⎜ ⎟
⎝ 0 0 λxe λye ⎠
⎛1 0 0 0⎞
L= ⎜ ⎟
⎝0 0 1 0⎠
4.2.4 Matriz de rigidez del elemento respecto al sistema global es.
T
K := L ⋅ k ⋅ L
T
⎛ 1 0 0 0 ⎞ ⎡ 29500000⋅ 1 ⎛ 1 −1 ⎞⎤ ⎛ 1 0 0 0 ⎞
K = ⎜ ⎟ ⋅⎢ ⋅⎜ ⎟⎥ ⋅ ⎜ ⎟
⎝0 0 1 0⎠ ⎣ 40 ⎝ −1 1 ⎠⎦ ⎝ 0 0 1 0 ⎠
⎛⎜ 1 0 −1 0 ⎞
⎟
0 ⎟ 29.5 × 106
K= ⎜
0 0 0
⋅
⎜ −1 0 1 0⎟ 40
⎜0 ⎟
⎝ 0 0 0⎠
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
7 de 21
MEF/ANÁLISIS DE ARMADURAS
4. Matriz de Rigidez de Elementos
MEF/ARMADURAS/DATOS
MEF/ARMADURAS/VECTOR DE CARGAS
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
4.1 Matriz de rigidez para el elemento
Elm := 2
4.1.1 la longitud del elemento y los
cosenos directores son
⎛ Le ⎞
⎜ λxe ⎟
⎜ ⎟
⎜ λye ⎟ := ni ← MEMBElm , 1
⎜ Ae ⎟
⎜ ⎟ nf ← MEMB
Elm , 2
⎝ Ee ⎠
prop ← MEMB
Elm , 3
xi ← NODE
ni , 1
yi ← NODE
ni , 2
xf ← NODE Con el programa se obtiene para el elemento
nf , 1
Elm = 2
yf ← NODE
nf , 2 • Longitud del elemento
Le ←
2
( xf − xi) + ( yf − yi)
2 Le = 30
xf − xi • Cosenos directores
λxe ←
Le
λxe = 0
yf − yi
λye ←
Le λye = −1
Ae ← PROP
prop , 1 • Propiedades
Ee ← PROP Ae = 1
prop , 2
T 7
( Le λxe λye Ae Ee ) Ee = 2.95 × 10
4.2.2 La matriz de rigidez en coordenada local está dado por:
Ee⋅ Ae ⎛ 1 −1 ⎞
k := ⋅⎜ ⎟
Le ⎝ −1 1 ⎠
29500000⋅ 1 ⎛ 1 −1 ⎞
k = ⋅⎜ ⎟
30 ⎝ −1 1 ⎠
⎛ 9.833 × 105 −9.833 × 105 ⎞
k=⎜ ⎟
⎜ 5 5 ⎟
⎝ −9.833 × 10 9.833 × 10 ⎠
8 de 21
MEF/ANÁLISIS DE ARMADURAS
4.2.3 La matriz de transformación de desplazamientos está dado por:
⎛ λxe λye 0 0 ⎞
L := ⎜ ⎟
⎝ 0 0 λxe λye ⎠
⎛ 0 −1 0 0 ⎞
L= ⎜ ⎟
⎝ 0 0 0 −1 ⎠
4.2.4 Matriz de rigidez del elemento respecto al sistema global es.
T
K := L ⋅ k ⋅ L
T
⎛ 0 −1 0 0 ⎞ ⎡ 29500000⋅ 1 ⎛ 1 −1 ⎞⎤ ⎛ 0 −1 0 0 ⎞
K = ⎜ ⎟ ⋅⎢ ⋅⎜ ⎟⎥ ⋅ ⎜ ⎟
⎝ 0 0 0 −1 ⎠ ⎣ 30 ⎝ −1 1 ⎠⎦ ⎝ 0 0 0 −1 ⎠
⎛⎜ 0 0 ⎞⎟
0 0
1 0 −1 ⎟ 29.5 × 106
K= ⎜
0
⋅
⎜0 0 0 0 ⎟ 30
⎜0 ⎟
⎝ −1 0 1 ⎠
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
9 de 21
MEF/ANÁLISIS DE ARMADURAS
4. Matriz de Rigidez de Elementos
MEF/ARMADURAS/DATOS
MEF/ARMADURAS/VECTOR DE CARGAS
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
4.1 Matriz de rigidez para el elemento
Elm := 3
4.1.1 la longitud del elemento y los
cosenos directores son
⎛ Le ⎞
⎜ λxe ⎟
⎜ ⎟
⎜ λye ⎟ := ni ← MEMBElm , 1
⎜ Ae ⎟
⎜ ⎟ nf ← MEMB
Elm , 2
⎝ Ee ⎠
prop ← MEMB
Elm , 3
xi ← NODE
ni , 1
yi ← NODE
ni , 2
xf ← NODE Con el programa se obtiene para el elemento
nf , 1
Elm = 3
yf ← NODE
nf , 2 • Longitud del elemento
Le ←
2
( xf − xi) + ( yf − yi)
2 Le = 50
xf − xi • Cosenos directores
λxe ←
Le
λxe = 0.8
yf − yi
λye ←
Le λye = 0.6
Ae ← PROP
prop , 1 • Propiedades
Ee ← PROP Ae = 1
prop , 2
T 7
( Le λxe λye Ae Ee ) Ee = 2.95 × 10
4.2.2 La matriz de rigidez en coordenada local está dado por:
Ee⋅ Ae ⎛ 1 −1 ⎞
k := ⋅⎜ ⎟
Le ⎝ −1 1 ⎠
29500000⋅ 1 ⎛ 1 −1 ⎞
k = ⋅⎜ ⎟
50 ⎝ −1 1 ⎠
⎛ 5.9 × 105 −5.9 × 105 ⎞
k=⎜ ⎟
⎜ 5 5 ⎟
⎝ −5.9 × 10 5.9 × 10 ⎠
10 de 21
MEF/ANÁLISIS DE ARMADURAS
4.2.3 La matriz de transformación de desplazamientos está dado por:
⎛ λxe λye 0 0 ⎞
L := ⎜ ⎟
⎝ 0 0 λxe λye ⎠
⎛ 0.8 0.6 0 0 ⎞
L= ⎜ ⎟
⎝ 0 0 0.8 0.6 ⎠
4.2.4 Matriz de rigidez del elemento respecto al sistema global es.
T
K := L ⋅ k ⋅ L
T
⎛4 3 0 0⎞ ⎛4 3 0 0⎞
⎜5 5 ⎟ ⎡ 29500000⋅ 1 ⎛ 1 −1 ⎞⎤ ⎜ 5 5 ⎟
K = ⎜ ⎟ ⋅⎢ ⋅⎜ ⎟⎥ ⋅ ⎜ ⎟
⎜0 0 4 3⎟ ⎣ 50 ⎝ −1 1 ⎠⎦ ⎜ 0 0 4 3 ⎟
⎝ 5 5⎠ ⎝ 5 5⎠
⎛⎜ 0.64 0.48 −0.64 −0.48 ⎞
⎟
−0.48 −0.36 ⎟ 29.5 × 106
K= ⎜
0.48 0.36
⋅
⎜ −0.64 −0.48 0.64 0.48 ⎟ 50
⎜ −0.48 ⎟
⎝ −0.36 0.48 0.36 ⎠
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
11 de 21
MEF/ANÁLISIS DE ARMADURAS
4. Matriz de Rigidez de Elementos
MEF/ARMADURAS/DATOS
MEF/ARMADURAS/VECTOR DE CARGAS
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
4.1 Matriz de rigidez para el elemento
Elm := 4
4.1.1 la longitud del elemento y los
cosenos directores son
⎛ Le ⎞
⎜ λxe ⎟
⎜ ⎟
⎜ λye ⎟ := ni ← MEMBElm , 1
⎜ Ae ⎟
⎜ ⎟ nf ← MEMB
Elm , 2
⎝ Ee ⎠
prop ← MEMB
Elm , 3
xi ← NODE
ni , 1
yi ← NODE
ni , 2
xf ← NODE Con el programa se obtiene para el elemento
nf , 1
Elm = 4
yf ← NODE
nf , 2 • Longitud del elemento
Le ←
2
( xf − xi) + ( yf − yi)
2 Le = 40
xf − xi • Cosenos directores
λxe ←
Le
λxe = 1
yf − yi
λye ←
Le λye = 0
Ae ← PROP
prop , 1 • Propiedades
Ee ← PROP Ae = 1
prop , 2
T 7
( Le λxe λye Ae Ee ) Ee = 2.95 × 10
4.2.2 La matriz de rigidez en coordenada local está dado por:
Ee⋅ Ae ⎛ 1 −1 ⎞
k := ⋅⎜ ⎟
Le ⎝ −1 1 ⎠
29500000⋅ 1 ⎛ 1 −1 ⎞
k = ⋅⎜ ⎟
40 ⎝ −1 1 ⎠
⎛ 7.375 × 105 −7.375 × 105 ⎞
k=⎜ ⎟
⎜ 5 5 ⎟
⎝ −7.375 × 10 7.375 × 10 ⎠
12 de 21
MEF/ANÁLISIS DE ARMADURAS
4.2.3 La matriz de transformación de desplazamientos está dado por:
⎛ λxe λye 0 0 ⎞
L := ⎜ ⎟
⎝ 0 0 λxe λye ⎠
⎛1 0 0 0⎞
L= ⎜ ⎟
⎝0 0 1 0⎠
4.2.4 Matriz de rigidez del elemento respecto al sistema global es.
T
K := L ⋅ k ⋅ L
T
⎛ 1 0 0 0 ⎞ ⎡ 29500000⋅ 1 ⎛ 1 −1 ⎞⎤ ⎛ 1 0 0 0 ⎞
K = ⎜ ⎟ ⋅⎢ ⋅⎜ ⎟⎥ ⋅ ⎜ ⎟
⎝0 0 1 0⎠ ⎣ 40 ⎝ −1 1 ⎠⎦ ⎝ 0 0 1 0 ⎠
⎛⎜ 1 0 −1 0 ⎞
⎟
0 ⎟ 29.5 × 106
K= ⎜
0 0 0
⋅
⎜ −1 0 1 0⎟ 40
⎜0 ⎟
⎝ 0 0 0⎠
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
13 de 21
MEF/ANÁLISIS DE ARMADURAS
4-1. Matriz de Rigidez Ensamblado
MEF/ARMADURAS/DATOS
MEF/ARMADURAS/VECTOR DE CARGAS
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
MEF/ARMADURAS/MATRIZ RIGIDEZ ENSAMBLADO
Programa que ensambla la matriz de rigidez de la estructura "K"
K := N ← 2 rows( NODE)
f (x , y) ← 0
K ← matrix( N , N , f )
for i ∈ 1 .. rows( MEMB)
ni ← MEMB
i, 1
nf ← MEMB
i, 2
prop ← MEMB
i, 3
xi ← NODE
ni , 1
yi ← NODE
ni , 2
xf ← NODE
nf , 1
yf ← NODE
nf , 2
2 2
Le ← ( xf − xi) + ( yf − yi)
xf − xi
λxe ←
Le
yf − yi
λye ←
Le
Ae ← PROP
prop , 1
Ee ← PROP
prop , 2
Ee⋅ Ae ⎛ 1 −1 ⎞
k← ⋅⎜ ⎟
Le ⎝ −1 1 ⎠
⎛ λxe λye 0 0 ⎞
L← ⎜ ⎟
⎝ 0 0 λxe λye ⎠
T
Ke ← L ⋅ k ⋅ L
"Columna 2*ni-1"
Ke
1, 1
K ←K +
2⋅ ni− 1 , 2⋅ ni− 1 2⋅ ni− 1 , 2⋅ ni− 1 2
K ←K + Ke
2ni , 2ni− 1 2ni , 2ni− 1 2, 1
K ←K + Ke
2⋅ nf − 1 , 2⋅ ni− 1 2⋅ nf − 1 , 2⋅ ni− 1 3, 1
K ←K + Ke
2nf , 2ni− 1 2nf , 2ni− 1 4, 1
"Columna 2*ni"
14 de 21
MEF/ANÁLISIS DE ARMADURAS
Ke
2, 2
K ←K +
2ni , 2⋅ ni 2ni , 2⋅ ni 2
K ←K + Ke
2nf − 1 , 2ni 2nf − 1 , 2ni 3, 2
K ←K + Ke
2nf , 2ni 2nf , 2ni 4, 2
"Columna 2*nf-1"
Ke
3, 3
K ←K +
2nf − 1 , 2nf − 1 2nf − 1 , 2nf − 1 2
K ←K + Ke
2nf , 2nf − 1 2nf , 2nf − 1 4, 3
"Columna 2*nf"
Ke
4, 4
K ←K +
2nf , 2nf 2nf , 2nf 2
T
K←K + K
K
⎛ 22.68 5.76 −15 0 −7.68 −5.76 0 0⎞
⎜ 5.76 4.32 0 0
⎟
−5.76 −4.32 0 0
⎜ ⎟
⎜ −15 0 15 0 0 0 0 0⎟
⎜ 0 0 0 20 0 −20 0 0 ⎟ 29.5⋅ 106
K= ⎜ ⎟⋅
−5.76 0 0 22.68 5.76 −15 0
⎜ −7.68 ⎟ 600
⎜ −5.76 −4.32 0 −20 5.76 24.32 0 0 ⎟
⎜ 0 0 0 0 −15 0 15 0 ⎟
⎜ ⎟
⎝ 0 0 0 0 0 0 0 0⎠
MEF/ARMADURAS/MATRIZ RIGIDEZ ENSAMBLADO
15 de 21
MEF/ANÁLISIS DE ARMADURAS
5. Estableciendo las Condiciones de Frontera
MEF/ARMADURAS/DATOS
MEF/ARMADURAS/VECTOR DE CARGAS
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
MEF/ARMADURAS/MATRIZ RIGIDEZ ENSAMBLADO
MEF/ARMADURAS/ESTABLECE CONDICIONES DE FRONTERA
Modificando la matriz de rigidez "K" según las restricciones en los apoyos.
Km := Km ← K
for i ∈ 1 .. rows( SUPP)
n ← SUPP
i, 1
ux ← SUPP
i, 2
uy ← SUPP
i, 3
Km ← ∞ if ux = 1
2n− 1 , 2n− 1
Km ← ∞ if uy = 1
2n , 2n
Km
⎛⎜ 2.034 × 10302 −15 −7.68 −5.76 ⎞
5.76 0 0 0 ⎟
⎜ 302 ⎟
⎜ 5.76 2.034 × 10 0 0 −5.76 −4.32 0 0 ⎟
⎜ −15 0 15 0 0 0 0 0 ⎟
⎜ 302
⎟ 6
Km = ⎜
0 0 0 2.034 × 10 0 −20 0 0 ⎟ ⋅ 29.5⋅ 10
⎜ −7.68 −5.76 0 0 22.68 5.76 −15 0 ⎟ 600
⎜ ⎟
⎜ −5.76 −4.32 0 −20 5.76 24.32 0 0 ⎟
⎜ 302 ⎟
−15 2.034 × 10
⎜ 0 0 0 0 0 0
⎟
⎜ 302 ⎟
⎝ 0 0 0 0 0 0 0 2.034 × 10 ⎠
MEF/ARMADURAS/ESTABLECE CONDICIONES DE FRONTERA
16 de 21
MEF/ANÁLISIS DE ARMADURAS
5-1 Obteniendo Desplazamiento en Nudos
MEF/ARMADURAS/DATOS
MEF/ARMADURAS/VECTOR DE CARGAS
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
MEF/ARMADURAS/MATRIZ RIGIDEZ ENSAMBLADO
MEF/ARMADURAS/ESTABLECE CONDICIONES DE FRONTERA
MEF/ARMADURAS/DESPLAZAMIENTO EN NUDOS
Formando la matriz aumentada.
⎛⎜ 1 × 10302 2.832 −7.375 0 −3.776 −2.832 0 0 0 ⎞⎟
⎜ 302 ⎟
⎜ 2.832 1 × 10 0 0 −2.832 −2.124 0 0 0 ⎟
⎜ −7.375 0 7.375 0 0 0 0 0 0.2 ⎟
⎜ 302
⎟
augment( Km , F) = ⎜
0 0 0 1 × 10 0 −9.833 0 0 0 ⎟ 5
⋅ 10
⎜ −3.776 −2.832 0 0 11.151 2.832 −7.375 0 0 ⎟
⎜ ⎟
⎜ −2.832 −2.124 0 −9.833 2.832 11.957 0 0 −0.25 ⎟
⎜ 302 ⎟
−7.375 1 × 10
⎜ 0 0 0 0 0 0 0
⎟
⎜ 302 ⎟
⎝ 0 0 0 0 0 0 0 1 × 10 0 ⎠
el sistema de ecuaciones en su forma escalonada reducida
⎛⎜ 1 0 0 0 0 0 0 0 0 ⎞
⎟
⎜0 1 0 0 0 0 0 0 0 ⎟
⎜0 0 1 0 0 0 0 0 0.027 ⎟
⎜0 0 0 1 0 0 0 0 0
⎟
rref ( augment( Km , F) ) = ⎜ ⎟
⎜0 0 0 0 1 0 0 0 5.65 × 10
− 3⎟
⎜ ⎟
⎜0 0 0 0 0 1 0 0 −0.022 ⎟
⎜0 0 0 0 0 0 1 0 0 ⎟
⎜0 ⎟
⎝ 0 0 0 0 0 0 1 0 ⎠
los desplazamiento en los nudos son(según los grados de libertad):
〈rows( Km) + 1〉
Q := rref ( augment( Km , F) )
T (
Q = 0 × 10
0
0 × 10
0 −3
27.12 × 10 0 × 10
0
5.65 × 10
−3
−22.25 × 10
−3
0 × 10
0
0 × 10
0 )
MEF/ARMADURAS/DESPLAZAMIENTO EN NUDOS
17 de 21
MEF/ANÁLISIS DE ARMADURAS
6 Reacciones en los Apoyos
MEF/ARMADURAS/DATOS
MEF/ARMADURAS/VECTOR DE CARGAS
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
MEF/ARMADURAS/MATRIZ RIGIDEZ ENSAMBLADO
MEF/ARMADURAS/ESTABLECE CONDICIONES DE FRONTERA
MEF/ARMADURAS/DESPLAZAMIENTO EN NUDOS
MEF/ARMADURAS/REACCIONES EN APOYOS
Sustituyendo datos en la ecuación fundamental del elemento finito
R := K⋅ Q − F
T (
R = −1.583 × 10
4
3.125 × 10
3
0 2.187 × 10
4
0 0 −4.167 × 10
3
0 )
Ordenando las reacciones:
Columna 1: número de nudo
Columna 2: Reacción en la dirección "x" global
Columna 3: Reacción en la dirección "y" global
Ro := f ( x , y) ← 0
Ro ← matrix( rows( SUPP) , 3 , f )
for i ∈ 1 .. rows( SUPP)
n ← SUPP
i, 1
ux ← SUPP
i, 2
uy ← SUPP Reacciones en los nudos son:
i, 3
Ro ←R if ux = 1
i, 2 2n− 1 ⎛ 1 −1.583 × 104 3.125 × 103 ⎞
⎜ ⎟
Ro ←R if uy = 1 ⎜ ⎟
i, 3 2n Ro = 2 4
⎜ 0 2.187 × 10 ⎟
←n
Ro
i, 1 ⎜ 3 ⎟
⎝ 4 −4.167 × 10 0 ⎠
Ro
MEF/ARMADURAS/REACCIONES EN APOYOS
18 de 21
7 Esfuerzo Axial en los Elementos
MEF/ARMADURAS/DATOS
MEF/ARMADURAS/VECTOR DE CARGAS
MEF/ARMADURAS/MATRIZ RIGIDEZ C/E
MEF/ARMADURAS/MATRIZ RIGIDEZ ENSAMBLADO
MEF/ARMADURAS/ESTABLECE CONDICIONES DE FRONTERA
MEF/ARMADURAS/DESPLAZAMIENTO EN NUDOS
MEF/ARMADURAS/REACCIONES EN APOYOS
MEF/ARMADURAS/ESFUERZO AXIAL
Programa que obtiene los desplazamientos en todas las barras
σ := for i ∈ 1 .. rows( MEMB)
ni ← MEMB
i, 1
nf ← MEMB
i, 2
prop ← MEMB
i, 3
xi ← NODE
ni , 1
yi ← NODE
ni , 2
xf ← NODE
nf , 1
yf ← NODE
nf , 2
2 2
Le ← ( xf − xi) + ( yf − yi)
xf − xi
λxe ←
Le
yf − yi
λye ←
Le
Ee ← PROP
prop , 2
⎛ Q2ni−1 ⎞
⎜ ⎟ ⎛⎜ 2 × 104 ⎞⎟
⎜ Q2ni ⎟
q←⎜ ⎟ ⎜ 4⎟
⎜ Q2nf −1 ⎟ ⎜ −2.187 × 10 ⎟
σ=
⎜ Q ⎟ ⎜ 3⎟
⎝ 2nf ⎠ ⎜ −5.208 × 10 ⎟
⎜ 4.167 × 103 ⎟
σ ←
Ee
⋅ ( −λxe −λye λxe λye ) ⋅ q ⎝ ⎠
i Le
σ
MEF/ARMADURAS/ESFUERZO AXIAL
19 de 21
MEF/ANÁLISIS DE ARMADURAS
8- Comparando los resultados con el programa Sap2000 11.0.4.
8.1 Desplazamientos en los nudos.
Resultado con sap2000 11.0.4
Resultado con mathCAD.
8.2 reacciones en los apoyos.
Resultado con sap2000 11.0.4
Resultado con mathCAD
20 de 21
8.3 Esfuerzo axial en los elementos.
Resultado con sap2000 11.0.4
Resultado con mathCAD.
Conclusión:
Se Observa que los cálculos obtenidos son idénticos a los de sap2000 11.0.4
21 de 21
