Scribd
Upload
14 views5 pages

Tarea II

Uploaded by

BRYAM FELIPE SARI DOMINGUEZ
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
14 views5 pages

Tarea II

Uploaded by

BRYAM FELIPE SARI DOMINGUEZ
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

clc

clear

E = [ 2e5 2e5 2e5 2e5 8e5 2e5 2e5 2e5 2e5 ] ; % ( N / mm2 )
A = [ 150 200 250 200 250 200 100 100 120 ] ; % ( mm2 )
L = [3000 3000 5000 4000 5000 4000 4000 5000 3000 ] ; % (mm)
theta = [ 0 0 53.13 90 126.87 90 0 36.87 90 ] ; % (grados)

Operador para obtener EA/L para cada elemento

EAL = E.*A./L ;
disp(EAL)

1.0e+04 *

1.0000 1.3333 1.0000 1.0000 4.0000 1.0000 0.5000 0.4000 0.8000

Matrices locales de cada elemento

Ke1 = EAL ( 1 ) * [ 1 -1 ; -1 1 ] ;
Ke2 = EAL ( 2 ) * [ 1 -1 ; -1 1 ] ;
Ke3 = EAL ( 3 ) * [ 1 -1 ; -1 1 ] ;
Ke4 = EAL ( 4 ) * [ 1 -1 ; -1 1 ] ;
Ke5 = EAL ( 5 ) * [ 1 -1 ; -1 1 ] ;
Ke6 = EAL ( 6 ) * [ 1 -1 ; -1 1 ] ;
Ke7 = EAL ( 7 ) * [ 1 -1 ; -1 1 ] ;
Ke8 = EAL ( 8 ) * [ 1 -1 ; -1 1 ] ;
Ke9 = EAL ( 8 ) * [ 1 -1 ; -1 1 ] ;

Multiplicación de matrices

K1=[T_cel(theta(1))]'*Ke1*[T_cel(theta(1))]

K1 = 4×4
10000 0 -10000 0
0 0 0 0
-10000 0 10000 0
0 0 0 0

disp(K1)

10000 0 -10000 0
0 0 0 0
-10000 0 10000 0
0 0 0 0

K2=[T_cel(theta(2))]'*Ke2*[T_cel(theta(2))]

K2 = 4×4
104 ×
1.3333 0 -1.3333 0
0 0 0 0
-1.3333 0 1.3333 0

1
0 0 0 0

disp(K2)

1.0e+04 *

1.3333 0 -1.3333 0
0 0 0 0
-1.3333 0 1.3333 0
0 0 0 0

K3=[T_cel(theta(3))]'*Ke3*[T_cel(theta(3))]

K3 = 4×4
103 ×
3.6000 4.8000 -3.6000 -4.8000
4.8000 6.4000 -4.8000 -6.4000
-3.6000 -4.8000 3.6000 4.8000
-4.8000 -6.4000 4.8000 6.4000

disp(K3)

1.0e+03 *

3.6000 4.8000 -3.6000 -4.8000


4.8000 6.4000 -4.8000 -6.4000
-3.6000 -4.8000 3.6000 4.8000
-4.8000 -6.4000 4.8000 6.4000

K4=[T_cel(theta(4))]'*Ke4*[T_cel(theta(4))]

K4 = 4×4
0 0 0 0
0 10000 0 -10000
0 0 0 0
0 -10000 0 10000

disp(K4)

0 0 0 0
0 10000 0 -10000
0 0 0 0
0 -10000 0 10000

K5=[T_cel(theta(5))]'*Ke5*[T_cel(theta(5))]

K5 = 4×4
104 ×
1.4400 -1.9200 -1.4400 1.9200
-1.9200 2.5600 1.9200 -2.5600
-1.4400 1.9200 1.4400 -1.9200
1.9200 -2.5600 -1.9200 2.5600

disp(K5)

1.0e+04 *

2
1.4400 -1.9200 -1.4400 1.9200
-1.9200 2.5600 1.9200 -2.5600
-1.4400 1.9200 1.4400 -1.9200
1.9200 -2.5600 -1.9200 2.5600

K6=[T_cel(theta(6))]'*Ke6*[T_cel(theta(6))]

K6 = 4×4
0 0 0 0
0 10000 0 -10000
0 0 0 0
0 -10000 0 10000

disp(K6)

0 0 0 0
0 10000 0 -10000
0 0 0 0
0 -10000 0 10000

K7=[T_cel(theta(7))]'*Ke7*[T_cel(theta(7))]

K7 = 4×4
5000 0 -5000 0
0 0 0 0
-5000 0 5000 0
0 0 0 0

disp(K7)

5000 0 -5000 0
0 0 0 0
-5000 0 5000 0
0 0 0 0

K8=[T_cel(theta(8))]'*Ke8*[T_cel(theta(8))]

K8 = 4×4
103 ×
2.5600 1.9200 -2.5600 -1.9200
1.9200 1.4400 -1.9200 -1.4400
-2.5600 -1.9200 2.5600 1.9200
-1.9200 -1.4400 1.9200 1.4400

disp(K8)

1.0e+03 *

2.5600 1.9200 -2.5600 -1.9200


1.9200 1.4400 -1.9200 -1.4400
-2.5600 -1.9200 2.5600 1.9200
-1.9200 -1.4400 1.9200 1.4400

K9=[T_cel(theta(9))]'*Ke9*[T_cel(theta(9))]

K9 = 4×4

3
0 0 0 0
0 4000 0 -4000
0 0 0 0
0 -4000 0 4000

disp(K9)

0 0 0 0
0 4000 0 -4000
0 0 0 0
0 -4000 0 4000

Generación de la matriz de rigidez

K=zeros(12,12);

%Fila 1
K(1:2,1:2)=K1(1:2,1:2)+K3(1:2,1:2);
K(1:2,3:4)=K1(1:2,3:4);
K(1:2,5:6)=K(1:2,5:6);
K(1:2,7:8)=K3(1:2,3:4);
K(1:2,9:10)=K(1:2,9:10);
K(1:2,11:12)=K(1:2,11:12);

%FILA 2
K(3:4,1:2)=K1(3:4,1:2);
K(3:4,3:4)=K1(3:4,3:4)+K2(1:2,1:2)+K4(1:2,1:2);
K(3:4,5:6)=K2(1:2,3:4);
K(3:4,7:8)=K4(1:2,3:4);
K(3:4,9:10)=K(3:4,9:10);
K(3:4,11:12)=K(3:4,11:12);

%FIla 3
K(5:6,1:2)=K(5:6,1:2);
K(5:6,3:4)=K2(3:4,1:2);
K(5:6,5:6)=K2(3:4,3:4)+K5(3:4,3:4)+K6(1:2,1:2);
K(5:6,7:8)=K5(3:4,1:2);
k(5:6,9:10)=K6(1:2,3:4);
K(5:6,11:12)=K(5:6,11:12);

%FIla 4
K(7:8,1:2)=K3(3:4,1:2);
K(7:8,3:4)=K4(3:4,1:2);
K(7:8,5:6)=K5(1:2,3:4);
K(7:8,7:8)=K3(3:4,3:4)+K4(3:4,3:4)+K5(1:2,1:2)+K7(1:2,1:2)+K8(1:2,1:2);
k(7:8,9:10)=K7(1:2,3:4);
K(7:8,11:12)=K8(1:2,3:4);

%Fila 5
K(9:10,1:2)=K(9:10,1:2);
K(9:10,3:4)=K(9:10,3:4);

4
K(9:10,5:6)=K6(3:4,1:2);
K(9:10,7:8)=K7(3:4,1:2);
k(9:10,9:10)=K6(3:4,3:4)+K7(3:4,3:4)+K9(1:2,1:2);
K(9:10,11:12)=K9(1:2,3:4);

%Fila 6
K(11:12,1:2)=K(11:12,1:2);
K(11:12,3:4)=K(11:12,3:4);
K(11:12,5:6)=K(11:12,5:6);
K(11:12,7:8)=K8(3:4,1:2);
k(11:12,9:10)=K9(3:4,1:2);
K(11:12,11:12)=K8(3:4,3:4)+K9(3:4,3:4);

disp(K)

1.0e+04 *

1.3600 0.4800 -1.0000 0 0 0 -0.3600 -0.4800 0 0 0


0.4800 0.6400 0 0 0 0 -0.4800 -0.6400 0 0 0
-1.0000 0 2.3333 0 -1.3333 0 0 0 0 0 0
0 0 0 1.0000 0 0 0 -1.0000 0 0 0
0 0 -1.3333 0 2.7733 -1.9200 -1.4400 1.9200 0 0 0
0 0 0 0 -1.9200 3.5600 1.9200 -2.5600 0 0 0
-0.3600 -0.4800 0 0 -1.4400 1.9200 2.5560 -1.2480 0 0 -0.2560 -0.
-0.4800 -0.6400 0 -1.0000 1.9200 -2.5600 -1.2480 4.3440 0 0 -0.1920 -0.
0 0 0 0 0 0 -0.5000 0 0 0 0
0 0 0 0 0 -1.0000 0 0 0 0 0 -0.
0 0 0 0 0 0 -0.2560 -0.1920 0 0 0.2560 0.
0 0 0 0 0 0 -0.1920 -0.1440 0 0 0.1920 0.

You might also like