K=zeros(2*14); %globalna macierz pozorna
d=zeros(13,1); % wektor średnic wału
Jd=zeros(13,1); % wektor momnetów bezwładności przekrojów
Ld=zeros(13,1); % długości elementów
P=zeros(13*2,1); %wektor sił uogólnionych (siły skupione i momenty) z uwzględnieniem war. brzegowych
k=zeros(13,4,4); % macierz sztywności elementów
% dane do programu wal
E=2.1e5;
%średnice elementów (w mm)W tym miejscu wpisujesz średnice swoich elementów
d(1)=25; d(2)=40; d(3)=52; d(4)=52; d(5)=78; d(6)=96; d(7)=96;
d(8)= 112; d(9)=96; d(10)=75; d(11)=75; d(12)=60; d(13)=52;
% długości elementów (w mm) W tym miejscu wpisujesz długość swoich elementów
Ld(1)=10; Ld(2)=115; Ld(3)=35; Ld(4)=35; Ld(5)=80; Ld(6)=35; Ld(7)=35;
Ld(8)=8; Ld(9)=242; Ld(10)=15; Ld(11)=15; Ld(12)=150; Ld(13)=35;
%
for i=1:l3, Jd(i)=pi*d(i)*d(i)*d(i)*d(i)/64; end
% 0bciazenie zewnetrzne w N lub Nmm Tu wpisujesz siły [N] i momenty[Mmm].
Siły są w numerach nieparzystych momenty w parzystych.
Dla przykładu testowego siła jest w węźle 4 to siłę wpisujemy w P(7).
Jak moment byłby w węźle 5 to byśmy wpisali go w P(10).
Tylko jest jeden haczyk 😊. W węźle jest łożysko i w tym miejscu nie ma ugięcia, dlatego od tej liczby jak wynika numeru węzła trzeba odjąć tyle ile łożysk jest po lewej stronie. Dla przykładu testowego siłę wpisujemy w
P(7-1). Jakby np. moment byłby w węźle 7, czyli numer przypada 14 to odejmujemy dwa i wartość wpisujemy w P(12)=…
P(6)=-603; P(7)=129; P(12)=1891; P(13)=122; P(26)=-3900; P(27)=-164;
%
for i=1:13
J=Jd(i); L=Ld(i);
k(i,:,:)=(2*E*J/L^3)*[6, 3*L, -6, 3*L; 3*L, 2*L^2, -3*L, L^2; -6, -3*L, 6, -3*L; 3*L, L^2, -3*L, 2*L^2];
end
%
for ip=0:2:(13-1)*2
for i = 1:4
for j = 1:4
ks=0.5*ip+1;
K(i+ip,j+ip) = K(i+ip,j+ip) + k(ks,i,j);
end
end
end
%
% rzeczywista macierz globalna - macierz pozorna po uwzględnieniu warunków brzegowych
% (poprzez- wykreslenie wierszy i kolumn)% x,y położenie łożysk
Dla przykładu pierwsze łożysko jest w węźle drugim to x=2*2-1=3, drugie łożysko jest w węźle 6 to y=2*6-1=11
K([1 21],:)=[];
K(:,[1 21])=[];
%rozwiazanie rownania macierzowego Kx=P
%x - ugiecia i katy ugiec w wezlach
x = K \ P; %lub alternatywnie x=inv(K)*P;
%