Scribd
Upload
16 views

% Datos Delproblema: Ejercicio 5

This document contains the calculations to determine stresses in a structural problem. Stiffness matrices K1, K2, and K3 are formed using modulus of elasticity E, cross-sectional area A, and length values. These are summed to form the total stiffness matrix K. Load vector f and displacement vector u are determined by inverting the mass matrix M. Stresses EZ12, EZ23, and EZ24 are then calculated using u, E, and trigonometric functions of the beam angles.

Uploaded by

Rony Ynoquio Rosales
Copyright
© © All Rights Reserved
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
16 views

% Datos Delproblema: Ejercicio 5

This document contains the calculations to determine stresses in a structural problem. Stiffness matrices K1, K2, and K3 are formed using modulus of elasticity E, cross-sectional area A, and length values. These are summed to form the total stiffness matrix K. Load vector f and displacement vector u are determined by inverting the mass matrix M. Stresses EZ12, EZ23, and EZ24 are then calculated using u, E, and trigonometric functions of the beam angles.

Uploaded by

Rony Ynoquio Rosales
Copyright
© © All Rights Reserved
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
You are on page 1/ 4

EJERCICIO 5:

>> % DATOS DELPROBLEMA


>> E=210*10^9;
>> A=5*10^(-4);
>> P=2700;
>> Q=2750;
>> L1=10.4;
>> L2=10;
>> L3=14.14;
>> S1=sind(16);
>> C1=cosd(16);
>> D1=S1^2;
>> F1=C1^2;
>> G1=S1*C1;
>> S2=sind(180);
>> C2=cosd(180);
>> D2=S2^2;
>> F2=C2^2;
>> G2=S2*C2;
>> S3=sind(135);
>> C3=cosd(135);
>> D3=S3^2;
>> F3=C3^2;
>> G3=S3*C3;
>> a=[E*A]/L1;
>> b=[E*A]/L2;

>> c=[E*A]/L3;
>> %MATRIZ DE RIGIDEZ 1
>> K1= a*[F1 G1 -F1 -G1 0 0 0 0;G1 D1 -G1 -D1 0 0 0 0; -F1 -G1 F1 G1 0 0 0 0;
-G1 -D1 G1 D1 0 0 0 0; 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0; 0 0 0 0
0 0 0 0]

>> %MATRIZ DE RIGIDEZ 2


>> K2= b*[0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0; 0 0 F2 G2 -F2 -G2 0 0; 0 0 G2
D2 -G2 -D2 0 0; 0 0 -F2 -G2 F2 G2 0 0; 0 0 -G2 -D2 G2 D2 0 0; 0 0 0 0 0
0 0 0; 0 0 0 0 0 0 0 0]

>> K3= c*[0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0; 0 0 F3 G3 0 0 -F3 -G3; 0 0 G3


D3 0 0 -G3 -D3; 0 0 0 0 0 0 0 0; 0 0 0 0 0 0 0 0; 0 0 -F3 -G3 0 0 F3 G3; 0
0 -G3 -D3 0 0 G3 D3]

>> K=K1+K2+K3

>> M=1.0e+07 *[2.3542


M=

23542000

-1038000

-1038000

4480000

>> f=[2750; -2700]

f=
2750
-2700

>> u=inv(M)*f
u=

1.0e-03 *
0.0912
-0.5816

-0.1038; -0.1038

0.4480]

CALCULO DE ESFUERZOS
>> EZ12=(E/L1)*[-1 1]*[C1 S1 0 0; 0 0 C1 S1]*[0; 0; 1.0e-03 *0.0912;
1.0e-03 *-0.5816]
EZ =
-1.4668e+06
>> EZ23=(E/L2)*[-1 1]*[C2 S2 0 0; 0 0 C2 S2]*[ 1.0e-03 *0.0912; 1.0e03 *-0.5816; 0; 0]
EZ23 =
1.9152e+06
>> EZ24=(E/L3)*[-1 1]*[C3 S3 0 0; 0 0 C3 S3]*[ 1.0e-03 *0.0912; 1.0e03 *-0.5816; 0; 0]
EZ24 =
7.0655e+06

You might also like