Scribd
Upload
17 views6 pages

% Datos de Longitud de Barras % Barra 1-2 L1 ( (36+100+4) ) (0.5)

The document contains calculations to determine the stiffness matrix (k1) for bar 1-2 of a structural system. It first calculates the length (L) of each bar using the Pythagorean theorem. It then determines direction cosines (Cx, Cy, Cz) for each bar. Finally, it calculates the stiffness coefficient (a) for bar 1-2 and populates the 12x12 k1 matrix using the calculated values of a, L, and the direction cosines.

Uploaded by

Rony Ynoquio Rosales
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
17 views6 pages

% Datos de Longitud de Barras % Barra 1-2 L1 ( (36+100+4) ) (0.5)

The document contains calculations to determine the stiffness matrix (k1) for bar 1-2 of a structural system. It first calculates the length (L) of each bar using the Pythagorean theorem. It then determines direction cosines (Cx, Cy, Cz) for each bar. Finally, it calculates the stiffness coefficient (a) for bar 1-2 and populates the 12x12 k1 matrix using the calculated values of a, L, and the direction cosines.

Uploaded by

Rony Ynoquio Rosales
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 6

>> % Datos de Longitud de barras

% Barra 1-2
L1= [(36+100+4)]^(0.5)

% Barra 1-3
L2= [(40)]^(0.5)
% Barra 1-4
L3= [(4+16+16)]^(0.5)

L1 =

11.8322

L2 =

6.3246

L3 =

>> % Cx; Cy; Cz para Barra 1-2


>> Cx1= (0-6)/L1

Cx1 =

-0.5071

>> % Cx; Cy; Cz para Barra 1-3


>> Cy1= (14-4)/L1

Cy1 =

0.8452

>> Cz1= (10-8)/L1

Cz1 =

0.1690

>> % Cx; Cy; Cz para Barra 1-3


>> Cx2= (0-6)/L2

Cx2 =

-0.9487

>> Cy2= (2-4)/L2

Cy2 =

-0.3162

>> Cz2= (0)/L2

Cz2 =

>> % Cx; Cy; Cz para Barra 1-4


>> Cx3= (8-6)/L3

Cx3 =

0.3333

>> Cy3= (0-4)/L3

Cy3 =

-0.6667

>> Cz3= (12-8)/L3

Cz3 =

0.6667

>> % Datos del ejercicio


>> E= 250 *10^9

E=

2.5000e+11

>> A= 4*10^(-4);
>> P= 1200;
>> Q=1800;
>> % Desarrollo de la matriz k1 de la barra 1-2
>> a= (E*A)/L1

a=

8.4515e+06

>> b= (E*A)/L2

b=

1.5811e+07

>> c= (E*A)/L3

c=

1.6667e+07

>> % Matriz k1 de 1-2


>> k1= a*[(Cx1)^2 (Cx1)*(Cy1) Cx1*Cz1 -(Cx1)^2 -(Cx1)*(Cy1) -(Cx1*Cz1) 0 0 0 0 0 0;
(Cx1)*(Cy1) (Cy1)^2 (Cy1*Cz1) -(Cx1)*(Cy1) -(Cy1)^2 -(Cy1*Cz1) 0 0 0 0 0 0; (Cx1*Cz1)
(Cy1*Cz1) (Cz)^2 -(Cx1*Cz1) -(Cy1*Cz1) -(Cz)^2 0 0 0 0 0 0; -(Cx1)^2 -(Cx1)*(Cy1) -Cx1*Cz1
(Cx1)^2 (Cx1)*(Cy1) (Cx1*Cz1) 0 0 0 0 0 0; -(Cx1)*(Cy1) -(Cy1)^2 -(Cy1*Cz1) (Cx1)*(Cy1)
(Cy1)^2 (Cy1*Cz1) 0 0 0 0 0 0; -(Cx1*Cz1) -(Cy1*Cz1) -(Cz)^2 (Cx1*Cz1) (Cy1*Cz1) (Cz)^2 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 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 0 0 0 0 0 0]
Undefined function or variable 'Cz'.

>> k1= a*[(Cx1)^2 (Cx1)*(Cy1) Cx1*Cz1 -(Cx1)^2 -(Cx1)*(Cy1) -(Cx1*Cz1) 0 0 0 0 0 0;


(Cx1)*(Cy1) (Cy1)^2 (Cy1*Cz1) -(Cx1)*(Cy1) -(Cy1)^2 -(Cy1*Cz1) 0 0 0 0 0 0; (Cx1*Cz1)
(Cy1*Cz1) (Cz1)^2 -(Cx1*Cz1) -(Cy1*Cz1) -(Cz1)^2 0 0 0 0 0 0; -(Cx1)^2 -(Cx1)*(Cy1) -Cx1*Cz1
(Cx1)^2 (Cx1)*(Cy1) (Cx1*Cz1) 0 0 0 0 0 0; -(Cx1)*(Cy1) -(Cy1)^2 -(Cy1*Cz1) (Cx1)*(Cy1)
(Cy1)^2 (Cy1*Cz1) 0 0 0 0 0 0; -(Cx1*Cz1) -(Cy1*Cz1) -(Cz1)^2 (Cx1*Cz1) (Cy1*Cz1) (Cz1)^2 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 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 0 0 0 0 0 0]

k1 =

1.0e+06 *

Columns 1 through 9

2.1733 -3.6221 -0.7244 -2.1733 3.6221 0.7244

-3.6221 6.0368 1.2074 3.6221 -6.0368 -1.2074

-0.7244 1.2074 0.2415 0.7244 -1.2074 -0.2415

-2.1733 3.6221 0.7244 2.1733 -3.6221 -0.7244

3.6221 -6.0368 -1.2074 -3.6221 6.0368 1.2074

0.7244 -1.2074 -0.2415 -0.7244 1.2074 0.2415


0

Columns 10 through 12

You might also like