Analisis estructural

TRABAJO DOMICILIARIO DE ANALISIS ESTRUCTURAL

GRADOS DE LIBERTAD: GDL=4

MATRIZ DE RIGIDEZ DE CADA ELEMENTO:

ELEMENTO AB
E=2*107 Ton/m2 SECCION AB= 0.3*0.4
IAB= 0.3*(0.4)312 = 0.0016

1 6m 2

1
1
2
2
1
1
2
2
K(1) = 4EI/L2EI/L2EI/L4EI/L =4(32000)/62(32000)/62(32000)/64(32000)/6 = 21333.310666.6710666.6721333.3

ELEMENTO BC
E=2*107 Ton/m2 SECCION AB= 0.3*0.45
IAB= 0.3*(0.45)312 = 0.0022781

2 4m 3
2
2
3
3
2
2

3
3
K(2) = 4EI/L2EI/L2EI/L4EI/L = 4(45562.5)/42(45562.5)/42(45562.5)/44(45562.5)/4 = 45562.522781.2522781.2545562.5

ELEMENTO CD
E=2*107 Ton/m2 SECCION AB= 0.3*0.3
IAB= 0.3*(0.3)312 =0.000675

3 5m 4
3
3
4
4
3
3

4
4
K(3) = 4EI/L2EI/L2EI/L4EI/L = 4(13500)/52(13500)/52(13500)/54(13500)/5 = 108005400540010800

MATRIZ RIGIDEZ DEL SISTEMA
KSISTEMA = K11K12K13K14K21K22K23K24K31K32K33K34K41K42K43K44
K11=21333.3, K12=10666.67, K13=0, K14= 0,
K21=10666.67, K22=21333.3+45562.5=66895.8, K23=22781.25,K24=0
K31=0, K32=22781.25, K33= 45562.5+10800=56362.5, K34=0
K41=0, K42=0, K43=5400, K44=10800

KSISTEMA = 21333.310666.670010666.6766895.822781.250022781.2556362.5540000540010800

VECTOR DE FUEZAS DEL SISTEMA
Fe =0000

VECTOR DE CARGAS DE EMPOTRAMIENTO DE CADA ELEMENTO
WL212
WL212

M°A = 12M°B = -12
Pab2L2
Pab2L2
F1 = 12-12

M°B = 4 M°C = - 4
F2 = 4-4
WL212
WL212

M°C= 6.25 M°D = -6.25

F3 = 6.25-6.25

VECTOR FUERZAS DE EMPOTRAMIENTO DEL SISTEMA

F° = 12-12+4-4+6.25-6.25 = 12-82.25-6.25

VECTOR DE FUERZAS TOTALES DEL SISTEMA

FT= Fe - F° = 0000 - 12-82.25-6.25 =-128-2.256.25
CALCULO DE LOS DESPLAZAMIENTOS

F=K*U Ke U=K-1*F

U = 21333.310666.670010666.6766895.822781.250022781.2556362.5540000540010800-1 * -128-2.256.25

U = -0.00071930.0003137-0.00023330.0006953

VECTOR DE FUERZAS INTERNAS DE CADA ELEMENTO
FT= F° + Ke * Ue

ELEMENTO AB
FT=12-12+21333.310666.6710666.6721333.3*-0.00071930.0003137 = 0-12.98ELEMENTO BC
FT=4-4+45562.522781.2522781.2545562.5*0.0003137-0.0002333 = 12.98-7.845
ELEMNTO CD
FT=6.25-6.25+108005400540010800*-10/3537727/37498 = 7.08540
HALLANDO LAS REACCIONES DE LA ESTRUCTURA

PROGRAMACION EN EL PROGRAMA MATLAB
fprintf(' PRACTICA DOMICILIARIA DE ESTRUCTURAL II\n\n')
fprintf('Nombre: Heber Huaman LLamoca CICLO: VIII SECCION: B\n\n')
fprintf('DATOSGENERALES DEL PROBLEMA\n\n')
fprintf('Ingreso de datos\n\n')
fprintf('Ingreso de longitud L(m) \n')
L1=input('ingrese (L1):');
L2=input('ingrese (L2):');
L3=input('ingrese (L3):');
fprintf('Ingreso de bases (m) \n')
B1=input('seccion (B1):');
B2=input('seccion (B2):');
B3=input('seccion (B3):');
fprintf('Ingreso las alturas\n')
H1=input('seccion (H1):');
H2=input('seccion(H2):');
H3=input('seccion (H3):');
E=input('Elasticidad (E):');
fprintf('Ingreso de cargas distribuidas (Ton/m)\n')
W1=input('carga (W1):');
W3=input('carga (W3):');
fprintf('Ingreso de cargas puntuales (Ton)\n')
P2=input('carga puntual (P2):');
fprintf('Hallando el momento de Inercia (Ton/m)\n')
I1=B1*(H1)^3/12
I2=B2*(H2)^3/12
I3=B3*(H3)^3/12
fprintf('Ingreso de Matrizrigidez del elemento(1)(Ton/m) \n')
K1=zeros(2,2);
K1(1,1)=4*E*I1/L1;
K1(1,2)= 2*E*I1/L1;
K1(2,1)=K1(1,2);
K1(2,2)=K1(1,1)
fprintf('Ingreso de Matriz rigidez del elemento(2)(Ton/m)\n')
k2=zeros(2,2);
K2(1,1)=4*E*I2/L2;
K2(1,2)= 2*E*I2/L2;
K2(2,1)=K2(1,2)
K2(2,2)=K2(1,1)
fprintf('Ingreso de Matriz rigidez del elemento(3)(Ton/m)\n')
k3=zeros(2,2);
K3(1,1)=4*E*I3/L3;...