PC051216

/ r

Kel J = Ke DedlllJUJ % agregacja macierzy sztywności układu

K=[kel 1 (4,4)+kel_2(2,2) kel_2 (2,4) ;ke1_2 (4,2)    kel_2(4

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

S2o=[    (-P1/2)+(p*L2/2

p*L2^A2/12)+(P1*L2 S3o= ,^r -P2/2 -P2*L3

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

h agregacja wektora obciążeń mi Ro=[Slo(4)+S2o (2

% rozwiązanie układu rownan met

P=R-Ro;

q=inv(K)*P

%wektor przemieszczeń koncow elementu Dl= [0    0    0 q(l) ] ' ;

D2=[0 q(l) 0 q(2)]';

D3=[0 q(2) 0 0]*;

%wektor sil przywezlowych Sl=kel_l*Dl+Slo

S2=kel_2*D2+S2o S3=kel 3*D3+S3o