L = [3 4] % [m]

L = 1×2
3 4

I = [6.75e-4 1.529e-5] % [m^4]

I = 1×2
10-3 ×
0.6750 0.0153

fr=[100*cosd(45); -100*cosd(45) ;0;] % [ kN ]

fr = 3×1
70.7107
-70.7107
0

E= [1.8e7 2e8] % kN/m2

E = 1×2
18000000 200000000

A = [0.0900 0.002364] % [m^2]

A = 1×2
0.0900 0.0024

K1=ke_portico(E(1),A(1),I(1),L(1),1)

K1 = 6×6
5400 0 -8100 -5400 0 -8100
0 540000 0 0 -540000 0
-8100 0 16200 8100 0 8100
-5400 0 8100 5400 0 8100
0 -540000 0 0 540000 0
-8100 0 8100 8100 0 16200

K2=[E(2)*A(2)/L(2) 0 0 -E(2)*A(2)/L(2) 0 0;
0 0 0 0 0 0;
0 0 0 0 0 0;
-E(2)*A(2)/L(2) 0 0 E(2)*A(2)/L(2) 0 0;
0 0 0 0 0 0;
0 0 0 0 0 0;]

K2 = 6×6
105 ×
1.1820 0 0 -1.1820 0 0
0 0 0 0 0 0
0 0 0 0 0 0
-1.1820 0 0 1.1820 0 0
0 0 0 0 0 0
0 0 0 0 0 0

Kg=zeros(9,9);
disp(Kg)

0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0

1
 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

%
Kg(1:3,1:3)=K1(1:3,1:3);
Kg(1:3,4:6)=K1(1:3,4:6);
%
Kg(4:6,1:3)=K1(4:6,1:3);
Kg(4:6,4:6)=K1(4:6,4:6)+K2(1:3,1:3);
Kg(4:6,7:9)=K2(1:3,4:6)

Kg = 9×9
105 ×
0.0540 0 -0.0810 -0.0540 0 -0.0810 0 0
0 5.4000 0 0 -5.4000 0 0 0
-0.0810 0 0.1620 0.0810 0 0.0810 0 0
-0.0540 0 0.0810 1.2360 0 0.0810 -1.1820 0
0 -5.4000 0 0 5.4000 0 0 0
-0.0810 0 0.0810 0.0810 0 0.1620 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

%
Kg(7:9,4:6)=K2(4:6,1:3);
Kg(7:9,7:9)=K2(4:6,4:6);
disp(Kg)

1.0e+05 *

0.0540 0 -0.0810 -0.0540 0 -0.0810 0 0 0

0 5.4000 0 0 -5.4000 0 0 0 0
-0.0810 0 0.1620 0.0810 0 0.0810 0 0 0
-0.0540 0 0.0810 1.2360 0 0.0810 -1.1820 0 0
0 -5.4000 0 0 5.4000 0 0 0 0
-0.0810 0 0.0810 0.0810 0 0.1620 0 0 0
0 0 0 -1.1820 0 0 1.1820 0 0
0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0

kr=Kg

kr = 9×9
105 ×
0.0540 0 -0.0810 -0.0540 0 -0.0810 0 0
0 5.4000 0 0 -5.4000 0 0 0
-0.0810 0 0.1620 0.0810 0 0.0810 0 0
-0.0540 0 0.0810 1.2360 0 0.0810 -1.1820 0
0 -5.4000 0 0 5.4000 0 0 0
-0.0810 0 0.0810 0.0810 0 0.1620 0 0
0 0 0 -1.1820 0 0 1.1820 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0

% Definir las filas y columnas que se eliminarán

2
filas_eliminar = [1, 2, 3, 7, 8, 9]; % Filas 1, 2, 3, 7, 8
columnas_eliminar = [1, 2, 3, 7, 8, 9]; % Columnas 1, 2, 3, 7, 8

kr(filas_eliminar, :) = [];
kr(:, columnas_eliminar) = [];

disp(kr)

1.0e+05 *

1.2360 0 0.0810
0 5.4000 0
0.0810 0 0.1620

disp(kr)

1.0e+05 *

1.2360 0 0.0810
0 5.4000 0
0.0810 0 0.1620

ur=kr\fr

ur = 3×1
10-3 ×
0.5915
-0.1309
-0.2957

u=[0;0;0;ur(1,1);ur(2,1);ur(3,1);0;0;0]

FG=Kg*u

FG = 9×1
-0.7985
70.7107
2.3955
70.7107
-70.7107
-0.0000
-69.9122
0
0

fl1=K1*[u(1,1);u(2,1);u(3,1);u(4,1);u(5,1);u(6,1)]

fl1 = 6×1
-0.7985
70.7107
2.3955
0.7985
-70.7107
-0.0000

fl2=K2*[u(4,1);u(5,1);u(6,1);u(7,1);u(8,1);u(9,1)]

3
fl2 = 6×1
69.9122
0
0
-69.9122
0
0