PC051218

Imacierz sztywności elem. 3 Jcel 3=ke beam_m3 (El3, L3, ' 11')

% agregacja macierzy sztywno; K= [kel_l (4,4) +kel_:

* wektor sil przywezlowych ocl p*L2*2/12)+(P1*L2

R= [0 0] ';

h agregacja wektora obc Ro=(Sio(4)+S2o(2) S2c % rozwiązanie układu rownan i

P=R-Ro;

q=inv(K)*P

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

D2=[0 q(1) 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

2(2,4);kel_2(4,2) kel_2(4,4j*k«i_3'2,2) ];

elementu