PC051212

%Piotr Pawlak grll Zadanie kontrolne nr 5 elear, clc Ll=4; %[m]

1,2=1; % [m]

L3=5;% [m]

E=200*10^A6;%[kN/m^A2]

11=1050*10"-8;%[m^A4]

12=2050*10^A-8; % [m^A4]

I3=750*10^A-8;%[m^A4]

EI1=E*I1;

EI2=E*I2;

EI3=E*I3; p=4.5;%[kN/m]

Pl=30;%[kN]

P2=25;%[kN]

%macierz sztywności elem. 1 kel_l=ke_beam_m3(Eli,LI,* 11¹);

%macierz sztywności elem. 2 kel_2=ke_beam_m3(EI2,L2,’11’) ;

%macierz sztywności elem. 3 kel_3=ke_beam_m3(EI3,L3,'11');

% agregacja macierzy sztywności układu

K=[kel_l(4,4)+kel_2(2,2) kel_2(2,4);kel_2 (4,2) kel_2(4,4)+kel_3 (2,2)]

%wektor sil przywezlowych od obciążeń przeslowych Slo= [ -Pl/2 -Pl*Ll/8 -Pl/2 Pl*Ll/8 ]';

S2o=[ (-P1/2) + (p*L2/2)    (-Pl*L2/8) + (p*L2^A2/12)    (-P1/2) + (p*L2/2)    (-

p*L2^A2/12)+(Pl*L2/8)    ]’;

S3o= [ -P2/2 -P2*L3/8 -P2/2 P2*L3/8 ]';

% wektor obciążeń węzłowych R=[0 0]