DEBER
ARMADURA PLANA
Objetivo
El presente trabajo tiene como objetivo mostrar la resolución de una armadura plana, mediante el
método de la rigidez y el uso del matlab
También se realizará una comparación entre los resultados obtenidos con los calculados con el
programa SAP 2000
CASO DE ESTUDIO.
Se considera una armadura plana como la las mostrada en la figura I, se considera que todos los
elementos tienen la misma sección transversal y el mismo módulo de elasticidad
FIGURA I
Caso de estudio
MODELACIÓN EN MATLAB
1) Obtención de K para cada miembro
Lo primero que se debe realizar es identificar el grado de indeterminación cinemática de la
estructura, las armaduras planas tienen 2 grados de libertad por nudo libre además el apoyo tipo
rodillo permite un desplazamiento adicional, se tiene 3 reacciones 2 para el apoyo fijo y una para el
tipo rodillo, los gdl se muestran en la figura anterior.
Una vez que tenemos identificado los gdl procedemos a numerar los elementos, hay que recordar
que vamos a pasar del sistema de coordenadas locales al sistema de coordenadas de nudos por lo
que debemos usar una matriz de trasformación. Para cada elemento se aplica la fórmula:
c1 cos( 1 )
s sen( )
T 1 1
c2 cos( 2 )
s 2 sen( 2 )
alfa = Ángulo del elemento en grados
beta1 = Ángulo en coordenadas locales del nudo 1 en grados
beta2 = Ángulo en coordenadas locales del nudo 2 en grados
La rigidez de un elemento armadura es EA/L, se considero EA igual a 1 por lo tanto k de cada
miembro es únicamente función de L para hallar la matriz de rigidez de los elementos en
coordenadas de nudos se aplica el triple producto mostrado a continuación:
K T k T T
El primer paso es obtener en el Matllab la matriz de rigidez de los miembros en coordenadas de
nudo.
En matlab (Obtención de K para cada miembro):
%ELEMENTO 1
Alfa=0;Beta1=0;Beta2=0;k=1;
T1=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]'
K1=T1*k*T1'
%ELEMENTO 2
Alfa=90;Beta1=0;Beta2=0;k2=1;
T2=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]'
K2=T2*k2*T2'
%ELEMENTO 3
Alfa=135;Beta1=30;Beta2=0;k3=1/sqrt(2);
T3=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]';
K3=T3*k3*T3'
%ELEMENTO 4
Alfa=45;Beta1=0;Beta2=0;k4=1/sqrt(2);
T4=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]';
K4=T4*k4*T4';
% ELEMENTO 5
Alfa=90;Beta1=30;Beta2=0;k5=1;
T5=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]';
K5=T5*k5*T5';
%ELEMENTO 6
Alfa=0;Beta1=0;Beta2=30;k6=1;
T6=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]';
K6=T6*k6*T6';
Matriz K6 obtenida
2) Obtención de Snt para la estructura
La matrices K de todos los elementos son cuadradas de 4x4, el siguiente paso es relacionar los
grados de libertad de la estructura con los desplazamientos D1, D2, D3 y D4 de los miembros para
ello usamos la matriz Gm
DGL + R
NUDO
ELEMENTO 6 ELEMENTO 5 ELEMENTO 4 ELEMENTO 3 ELEMENTO 2 ELEMENTO 1 GDL
1 2 3 4 5 6 7 8
S
D1 1 0 0 0 0 0 0 0
D2 0 1 0 0 0 0 0 0
D3 0 0 1 0 0 0 0 0
D4 0 0 0 1 0 0 0 0
D1 0 0 0 0 0 1 0 0
D2 0 0 0 0 0 0 1 0
D3 1 0 0 0 0 0 0 0
D4 0 1 0 0 0 0 0 0
D1 0 0 0 0 1 0 0 0
D2 0 0 0 0 0 0 0 1
D3 1 0 0 0 0 0 0 0
D4 0 1 0 0 0 0 0 0
D1 0 0 0 0 0 1 0 0
D2 0 0 0 0 0 0 1 0
D3 0 0 1 0 0 0 0 0
D4 0 0 0 1 0 0 0 0
D1 0 0 0 0 1 0 0 0
D2 0 0 0 0 0 0 0 1
D3 0 0 1 0 0 0 0 0
D4 0 0 0 1 0 0 0 0
D1 0 0 0 0 0 1 0 0
D2 0 0 0 0 0 0 1 0
D3 0 0 0 0 1 0 0 0
D4 0 0 0 0 0 0 0 1
Por ejemplo en el miembro 6 su D1 está relacionado con R6 su D2 con R7, su D3 con GDL5 y su D4
con R8, una vez que se tiene Gm´ debemos ubicar las rigideces obtenidas anteriormente en una
matriz KF que tiene una dimensión igual al número de miembros por 4, en este caso se tiene una
matriz de 24 x 24
Se puede obtener entonces la matriz de rigidez de la estructura
Snt=Gm*KF*Gm';
Hay que recordar que en esta matriz [Snt] de 8x8 están incluidos los 5 desplazamientos y las 3
restricciones en el siguiente orden:
1 2 3 4 5 6 7 8
1
2
Su SLR 3
[S] = 4
5
6
SRL SRR 7
8
snt
1 2 3 4 5 6 7 8
1 1,35 -0,35 -1,00 0,00 -0,13 0,00 0,00 0,48
2 -0,35 1,35 0,00 0,00 0,13 0,00 -1,00 -0,48
3 -1,00 0,00 1,35 0,35 0,00 -0,35 -0,35 0,00
4 0,00 0,00 0,35 1,35 -0,50 -0,35 -0,35 -0,87
5 -0,13 0,13 0,00 -0,50 1,05 -0,87 0,00 -0,18
6 0,00 0,00 -0,35 -0,35 -0,87 1,35 0,35 0,50
7 0,00 -1,00 -0,35 -0,35 0,00 0,35 1,35 0,00
8 0,48 -0,48 0,00 -0,87 -0,18 0,50 0,00 1,66
En matlab (Obtención de Snt):
Gm=0;% PEGO GM DEL EXCEL
Gm=Gm';
X=zeros(4);
KF=[K1 zeros(4) zeros(4) zeros(4) zeros(4) zeros(4)
zeros(4) K2 zeros(4) zeros(4) zeros(4) zeros(4)
zeros(4) zeros(4) K3 zeros(4) zeros(4) zeros(4)
zeros(4) zeros(4) zeros(4) K4 zeros(4) zeros(4)
zeros(4) zeros(4) zeros(4) zeros(4) K5 zeros(4)
zeros(4) zeros(4) zeros(4) zeros(4) zeros(4) K6]
Snt=Gm*KF*Gm';
3) Obtención de desplazamientos y reacciones
La fórmula que se aplica es la siguiente:
Delta=S\(P-SLR*DELTA_R)
El vector de cargas para este ejemplo es:
P= 0
El vector anterior indica que existe una carga puntual en el sentido del gdl 1 se representan
gráficamente en la siguiente figura:
P
30°
Se coloca el vector DeltaR
DELTA_R=[0 0 0]
El anterior vector indica que en R6,R7 y R8 se tienen desplazamientos 0, se obtiene también a
partir de Snt, las submatrices SLR SRR
Obtenemos finalmente estos desplazamientos
2,1357
0,5598
1,6955
-0,45
-0,02
En el método de la rigidez es fácil hallar las reacciones a partir de los
desplazamientos
R=SRL*Delta+SRR*DELTA_R
Tenemos las siguientes reacciones
R5=-0,42
R6=-1,00
R7= 1,15
En matlab (Desplazamientos y reacciones):
%desplazamientos en los grados de libertad de la estructura
P=[1 0 0 0 0]'; DELTA_R=[0 0 0]';
%S DE LA ESTR
S=Snt(1:5, 1:5);SLR=Snt(1:5, 6:8);SRL=Snt(6:8, 1:5);SRR=Snt(6:8, 6:8);
Delta=S\(P-SLR*DELTA_R);
R=SRL*Delta+SRR*DELTA_R
Dm
Finalmente comparamos lar reacciones obtenidas con los resultados calculados con SAP 2000 y
podemos verificar que los resultados son iguales
Reacciones obtenidas con SAP 2000
RUTINA
%ELEMENTO 1
Alfa=0;Beta1=0;Beta2=0;k=1;
T1=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]'
K1=T1*k*T1'
%ELEMENTO 2
Alfa=90;Beta1=0;Beta2=0;k2=1;
T2=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]'
K2=T2*k2*T2'
%ELEMENTO 3
Alfa=135;Beta1=30;Beta2=0;k3=1/sqrt(2);
T3=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]';
K3=T3*k3*T3'
%ELEMENTO 4
Alfa=45;Beta1=0;Beta2=0;k4=1/sqrt(2);
T4=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]';
K4=T4*k4*T4';
% ELEMENTO 5
Alfa=90;Beta1=30;Beta2=0;k5=1;
T5=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]';
K5=T5*k5*T5';
%ELEMENTO 6
Alfa=0;Beta1=0;Beta2=30;k6=1;
T6=[-cosd(Alfa-Beta1) -sind(Alfa-Beta1) cosd(Alfa-Beta2) sind(Alfa-Beta2)]';
K6=T6*k6*T6';
Gm=0;% PEGO GM DEL EXCEL
Gm=Gm';
X=zeros(4);
KF=[K1 zeros(4) zeros(4) zeros(4) zeros(4) zeros(4)
zeros(4) K2 zeros(4) zeros(4) zeros(4) zeros(4)
zeros(4) zeros(4) K3 zeros(4) zeros(4) zeros(4)
zeros(4) zeros(4) zeros(4) K4 zeros(4) zeros(4)
zeros(4) zeros(4) zeros(4) zeros(4) K5 zeros(4)
zeros(4) zeros(4) zeros(4) zeros(4) zeros(4) K6]
Snt=Gm*KF*Gm';
%desplazamientos en los grados de libertad de la estructura
P=[1 0 0 0 0]'; DELTA_R=[0 0 0]';
%S DE LA ESTR
S=Snt(1:5, 1:5);SLR=Snt(1:5, 6:8);SRL=Snt(6:8, 1:5);SRR=Snt(6:8, 6:8);
Delta=S\(P-SLR*DELTA_R);
R=SRL*Delta+SRR*DELTA_R