l=40; %Celková délka vagonu
a=10;
%a=(2-(2^0.5))*l; %optimální rozvor
%a=20.91096; %rozvor kol, hodnota, kdy je pruhyb uprostřed nulový a=20.91096;
b=(l-a)/2; %délka převislých konců
m=a+b; %pomocná vzdálenost
q=1; %spoj zat
E=210000000; %modul pružnosti [MPa]
sirka=0.3;
vyska=0.1;
J=(sirka*(vyska^3))/12; %moment
G=q/(2*E*J); %ohybová tuhost
%%
x1=(0:0.001:b);
x2=(b:0.001:m);
x3=(m:0.001:l);
x=[x1,x2,x3];
%%
Tx1=-q*x1;
Tx2=(q/2)*(l-2*x2);
Tx3=q*(l-x3);
Tx=[Tx1,Tx2,Tx3];
%%
Mx1=-q*((x1.^2)/2);
Mx2=(q/2)*(l*x2-l*b-x2.^2);     %((q*l)/2)*(x2-b)-q*((x2.*x2)/2);
Mx3=-(q/2)*(x3-l).^2;           %-q*(((l-x3).*(l-x3))/2);
Mx=[Mx1,Mx2,Mx3];
%%
nulax=(0:0.01:l);
nulay=0*nulax;
%%
hold off
subplot(4,1,1);
plot (x,Tx)
hold on
plot (nulax,nulay, 'Color', 'red')

hold off
subplot(4,1,2);
plot (x,Mx)
hold on
plot (nulax,nulay, 'Color', 'red')
%%
format long;
U=-(b^4)/12;
V=((b^3)/12)*(-4*l-b);
W=((m^2)/12)*(2*l*m-6*b*l-m^2);
X=((m^2)/12)*(4*l*m-(6*l^2)-m^2);
Y=((b^2)*l)/2;
Z=(m/2)*((2*l^2)-l*m-2*b*l);

vektor=[U;V;W;X;Y;Z];
matice=[b,1,0,0,0,0;0,0,b,1,0,0;0,0,m,1,0,0;0,0,0,0,m,1;1,0,-1,0,0,0;0,0,1,0,-1,0];
C=matice\vektor;
C1=C(1);
C2=C(2);
C3=C(3);
C4=C(4);
C5=C(5);
C6=C(6);


%%
fi1=G*(((x1.^3)/3)+C1);
fi2=G*(((x2.^3)/3)-((l*x2.^2)/2)+b*l*x2+C3);
fi3=G*((l*l*x3)-(l*x3.^2)+((x3.^3)/3)+C5);
fi=-[fi1,fi2,fi3];
%%
w1=G*(((x1.^4)/12)+C1*x1+C2);
w2=G*(((x2.^4)/12)-((l*x2.^3)/6)+(b*l*x2.^2)/2+C3*x2+C4);
w3=G*(((l*l*x3.^2)/2)-(l*x3.^3)/3+((x3.^4)/12)+C5*x3+C6);

w=-[w1,w2,w3];

%%
hold off
subplot(4,1,3);
plot (x,fi)
hold on
plot (nulax,nulay, 'Color', 'red')


hold off
subplot(4,1,4);
plot (x,w)
hold on
plot (nulax,nulay, 'Color', 'red')

wlpul=-G*((-(l^4)/64)+((l^3)/8).*b+(l/2).*C3+C4);

ax=(0:0.001:l);
bx=(l-ax)/2;
mx=bx+ax;



Mlpulx=(q/8)*(l*l-4*l*bx);

Ux=-(bx.^4)/12;
Vx=((bx.^3)/12).*(-4*l-bx);
Wx=((mx.^2)/12).*(2*l.*mx-6.*bx*l-mx.^2);
Xx=((mx.^2)/12).*(4.*l.*mx-(6*l^2)-mx.^2);
Yx=((bx.^2)*l)/2;
Zx=(mx./2).*((2*l^2)-l.*mx-2.*bx*l);


C3x=(Vx-Wx)./(bx-mx);
C4x=Vx-bx.*C3x;
C1x=Yx+C3x;
C5x=C3x-Z;
C2x=Ux-bx.*C1x;
C6x=Xx-mx.*C5x;

wlpulx=-G*((-(l^4)/64)+((l^3)/8).*bx+(l/2).*C3x+C4x);