www.budowlancezasrance.fora.pl

Seba

clear;clc

L1=6;
L2=5;
L3=5;
EI=1500;
p=45;
P1=80;
P2=60;
R=[0]';
kel_1=ke_beam_m3(EI, L1, '01');
kel_2=ke_beam_m3(EI, L2, '10');
kel_3=ke_beam_m3(EI, L3, '11');
K=zeros(1,1);
K(1,1)=kel_1(3,3)+kel_2(2,2)+kel_3(4,4);
S1o=[-3*p*L1/8 -5*p*L1/8 p*L1^2/8]';
S2o=[-11*P2/16 -3*P2*L2/16 -5*P2/16]';
S3o=[P1/2 P1*L3/8 P1/2 (-P1*L3)/8]';

Ro=[S1o(3)+S2o(2)+S3o(4)]';

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

D1=[0 0 q(1)]';
D2=[0 q(1) 0]';
D3=[0 0 0 q(1)]';

S1=kel_1*D1+S1o
S2=kel_2*D2+S2o
S3=kel_3*D3+S3o