PC051217

kel 3=ke_beam_m3 (EU, uj, ' 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 c Slo=[ -Pl/2 -PI*L1/

S2o=[ (-P1/2)+(p*L2/2) (-Pl*L2/8) ł (p * 1.2 " 2 / 12)    (-P1/2) + (p*L2

p*L2^A2/12)+(Pl*L2

s wektor obciążeń \

R=[0 0]';

\ agregacja \

Ro=[Sio(4 ; •

P=R-Ro; q=inv(K)* P

%wektor przemieszczeń koncow element Dl=[0    0    0 q(1) ] ' ;

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