PC051219

^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(El3,L3,'11');

£ agregacja macierzy sztywno:

K=[kel_l(4,4)+kel_2(2,2) kel

^wektor sil przywezlowych od <

Slo=[ -Pl/2 -P1*L1/

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

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

% agregacja wektora obciążeń miedzywezl Ro=[Slo(4)+S2o(2) S2o(4)+S3o(2)] ’ ;

% rozwiązanie układu rownan metody prze

P=R-Ro;

q=inv(K)*P

%%%%%%%%%%%%%%%%%%%%%%

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

D2= [0 q (1) 0 q (2) ]

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

%wektor sil przywezlowych Sl*kel l*Dl+Slo

mieś zczen