UNIVERSIDAD NACIONAL DE INGENIERA
FACULTAD DE INGENIERA MECNICA
Juan Alberto Segura Celis
ARMADURAS
Universidad Nacional de Ingeniera
NDICE
Enunciado del Problema..........................................................................2
Anlisis.....................................................................3
Tabla de Conectividad: Coordenadas y GDL...........................................3
Matriz de Rigidez Elemental de cada Elemento.......................................4
Ensamblado de la Matriz de Rigidez Estructural........................5
Desplazamientos Nodales.......................................................................5
Esfuerzos de cada Elemento....................................................................5
Fuerzas de Reaccin.........................................................................6
Diagrama de Flujo....................................................................................7
Caculo por Elementos Finitos
Universidad Nacional de Ingeniera
ARMADURA PLANA
PROBLEMA:
Considere la armadura de cuatro barras mostrada en la figura siguiente. Para
todos los elementos = 20.4 104 y = 645.162
a.
b.
c.
d.
e.
Determinar la matriz de rigidez elemental para cada elemento.
Ensamble la matriz de rigidez estructural K para toda la armadura.
Usando el mtodo de eliminacin, halle los desplazamientos nodales.
Recupere los esfuerzos para cada elemento.
Calcule las fuerzas de Reaccin.
GRFICO:
Caculo por Elementos Finitos
Universidad Nacional de Ingeniera
SOLUCIN
1. ANLISIS
2. TABLA DE CONECTIVIDAD: Coordenadas y GLD
NODO
1
2
3
4
X(mm)
0
1016
1016
0
Y(mm)
0
0
762
762
Elemento NODOS
GDL
Le
Ae
Ee
(1) (2) 1 2 3 4 (mm) (mm2) (N/mm2)
1
1
2 1 2 3 4 1016 645.16 20.4*104
3
5 6 3 4
1 2 5 6 1270 645.16 20.4*104
7 8 5 6 1016 645.16 20.4*104
Caculo por Elementos Finitos
762
645.16 20.4*104
Universidad Nacional de Ingeniera
3. MATRIZ DE RIGIDEZ ELEMENTAL DE CADA ELEMENTO ( )
(1)
= cos() =
(2)
= cos() =
=
[ 2
2 ] (3)
GDL
le
1 2 3 4
1016
5 6 3 4
762
-1
1 2 5 6
1270
0.8
0.6
7 8 5 6
1016
Elementos
1
1
4
20.4
10
645.16
0
1 =
[
1
1016
0
2
3 4
0 1 0 1
0 0 0 2
]
0 1 0 3
0 0 0 4
5
6
0 0
4
20.4
10
645.16
0 1
2 =
[
0 0
762
0 1
3
0
0
0
0
4
0
1
]
0
1
5
6
3
4
1
2
5
6
. 64
. 48 .64 .48
20.4 104 645.16 . 48
. 36 .48 .36
3
=
[
]
.64
.48
. 64
. 48
1270
.48 .36 . 48
. 36
7
1
4
20.4
10
645.16
0
4 =
[
1
1016
0
Caculo por Elementos Finitos
8
5 6
0 1 0 7
0 0 0 8
]
0 1 0 5
0 0 0 6
1
2
5
6
Universidad Nacional de Ingeniera
4. ENSAMBLADO DE LA MATRIZ DE RIGIDEZ ESTRUCTURAL (K)
La matriz de rigidez estructural K se ensambla a partir de ki, sumando las
contribuciones de cada elemento teniendo en cuenta su conectividad.
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
20.4 104 645.16
0
0
0
20
0
20
0
0
=
7.68 5.76
0
0
22.68 5.76 15 0
15240
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]
:
1
2
3
4
5
6
7
8
5. DESPLAZAMIENTOS NODALES
Las filas y columnas correspondientes a los grados de libertad que
corresponden a los soportes fijos se borran de la matriz, entonces la ecuacin del
elemento finito reducido ser:
=
(4)
: : : .
3
89181.82
0
0
20.4 104 645.16 15
[ 0 22.68 5.76 ] [5 ] = [
]
0
15240
111477.27
0
5.76 24.32
6
3
0.688
[5 ] = [ 0.143 ]
6
0.565
Luego, el vector desplazamiento nodal de la armadura ser:
[] = [0 0 0.688 0
0.143
0.565 0
0]
6. ESFUERZO DE CADA ELEMENTO
=
(5)
Elementos
GDL
1 2 3 4
-1
5 6 3 4
0.8
0.6
1 2 5 6
[0 0 0.143 0.565]
7 8 5 6
[0 0 0.143 0.565]
Caculo por Elementos Finitos
[0
0 0.688
[0.143 0.565
0]
0.688
0]
Universidad Nacional de Ingeniera
1 =
2 =
3
=
4
20.4 104
1016
20.4 104
762
20.4 104
1270
20.4 104
1016
[1 0 1
[0 1
0] [0 0
1 = 138.14
0]
0.688
0 1][0.143 0.565 0.688 0] ; 2 = 151.26
[0.8 0.6 0.8 0.6][0
[1 0 1
0][0
0.143
0.565] ; 3 = 36
0.143
4 = 28.71
0.565]
Luego, el vector esfuerzo de la armadura ser:
28.71]
[] = [138.14 151.26 36
7. FUERZAS DE REACCIN
Para hallar las fuerzas de reaccin en los apoyos, resolvemos la ecuacin
original del elemento finito mostrado a continuacin, tomando en cuenta
solamente las filas de K que correspondan a los grados de libertad 1, 2, 4, 7, 8
(grados de libertad que corresponden a los apoyos).
=
:
=0
: :
:
0
0
22.68 5.76 15 0 7.68 5.76 0 0
0.688
0
0 5.76 4.32 0 0
20.4 104 (645.16) 5.76 4.32
0
=
0
0
0
20
0
20
0 0
0.143
15240
0
0
0
0
15
0
15 0
0.565
[ 0
0
0
0
0
0
0 0]
0
[ 0 ]
1
70651.12
2
13965.44
= 4 = 97586.8
7
18524.22
]
[8 ] [
0
Caculo por Elementos Finitos
Universidad Nacional de Ingeniera
8. DIAGRAMA DE FLUJO DEL PROGRAMA
INICIO
Leer
datos de
entrada.
Para i=1
hasta N de
nodos
Ingresar coordenadas
de los nodos.
Calcular rea, N de filas de
cond_contorno(CC1)
Para i=1 hasta
2veces N de
nodos
Cont=0
Para j=1 hasta N de
filas de
cond_contorno(CC1)
Caculo por Elementos Finitos
Universidad Nacional de Ingeniera
SI
Si
i=CC(i,
Cont=1,
C2=CC1(i,2)
C1=CC1(i,1)
Si
SI
CC(i,1)=C1;
CC(i,2)=C2
N
CC(i,1)=0;
CC(i,2)=0
Para i=1
hasta N
elementos
Calcula Le, l, m, las
posiciones de la matriz de
rigidez global y su valor.
Caculo por Elementos Finitos
Universidad Nacional de Ingeniera
Para i=1 hasta
2veces N elementos.
SI
Si
i=CC(i,
1)
Q(i,1)=CC(i,2)
NO
Acumulamos
fuerzas
(FC=[FC;F(i)])
Para
j=1;2*Nnodo
s
SI
Si jCC(j,1)
acuh=[acuh,Kij(i,j)]
acumula filas
acuv=[acuv;acuh];
acumula columnas
Calcula los
desplazamientos generales
Q1=acuv\FC;
Caculo por Elementos Finitos
Universidad Nacional de Ingeniera
5
Para i=1;
2N nodos
Si
i==CC(i,1
)
Calcula las reacciones
r=Kij(i,1:2*nd)*Q-F(i,1);
R=[R;r i];
Para i=1 hasta N
de elementos
Calcula esfuerzos
Imprime
Desplazamientos,
reaciones y esfuerzos
Caculo por Elementos Finitos
10
