function [Tmax,Txy,Txz,a_rez,b_rez,deleni] = napeti_krut_obdelnik(a_rez,b_rez,deleni,E,poisson,M,n)
% VYPOČTY K BAKALÁŘSKÉ PRÁCI VALENTA TOMÁŠ
%
% napeti v krutu obdelnikoveho profilu

x = -a_rez/2:deleni:a_rez/2; % [m]
y = -b_rez/2:deleni:b_rez/2; % [m]

sizeM = size(M);
% % urcuje zda chceme vykreslit rez nebo ne a v jakem miste
% % se pro vypocet cele kliky
% kresleni = menu('Budete chtít vykreslit smyková napětí v řezu?','Ano - zadejte pozici řezu','Ne');
% if kresleni == 1
%     fprintf('\nZadejte číslo řezu pro vykreslení z intervalu  < 1 , %d >\n',sizeM(1));
%     pozice_rezu_kresleni = input('   : ');
% else
%     pozice_rezu_kresleni = 0;
% end
[X,Y] = meshgrid(x,y);

G = E/(2*(1+poisson));
sum_J=0;
for k=1:2:n
    for l=1:2:n
        sum_J = sum_J + 1/(k^2*l^2*(k^2*b_rez^2+l^2*a_rez^2));
    end
end
clear k l
J = ((2^8*a_rez^3*b_rez^3)/pi^6)*sum_J; % moment prurezu 

sum_Txy=0;
sum_Txz=0;
for k=1:2:n % vypocet sumy ve vzorci na str 58 tau pro obdelnik
    for l=1:2:n
        sum_Txy = sum_Txy +(((-1)^((k+l)/2))/(k*(k^2*b_rez^2+l^2*a_rez^2))).*cos((k.*pi.*X)./a_rez).*sin((l.*pi.*Y)./b_rez);
        sum_Txz = sum_Txz +(((-1)^((k+l)/2))/(l*(k^2*b_rez^2+l^2*a_rez^2))).*sin((k.*pi.*X)./a_rez).*cos((l.*pi.*Y)./b_rez);
    end
end
clear k l

for i=1:sizeM(1)
    alpha(i) = M(i,1)/(G*J);
    
    Txy = (((alpha(i)*G*2^5*a_rez^2*b_rez)/pi^3)*sum_Txy)/1e6; % smykove napeti [MPa]
    Txz = (((-alpha(i)*G*2^5*a_rez*b_rez^2)/pi^3)*sum_Txz)/1e6; % smykove napeti [MPa]
    
    T = sqrt(Txy.^2+Txz.^2); % celkove smykove napeti [MPa]
    Tmax(i) = max(max(T)); 

    
% -------- VYKRESLENI --------
%     if pozice_rezu_kresleni == i
%         figure('Position',[300 200 700 650])
%         surf(X,Y,T,'EdgeColor','interp')
%         fontsize = 15;
%         title('Smykove napeti ','FontSize',fontsize)
%         xlabel('a [m]','FontSize',fontsize)
%         ylabel('b [m]','FontSize',fontsize)
%         zlabel('T [MPa]','FontSize',fontsize)
%         set(gca, 'DataAspectRatio', [diff(get(gca, 'XLim')) diff(get(gca, 'XLim')) diff(get(gca, 'ZLim'))],'FontSize',fontsize)
%         print '-depsc' 'c:/Users/Tomáš Valenta/Škola/Vejška/_bachelory thesis/zprava/obrazky/12_smykove_napeti_T.eps' 
%         
%         figure('Position',[300 200 700 650])
%         surf(X,Y,Txy,'EdgeColor','interp')
%         title('Txy','FontSize',fontsize)
%         xlabel('a [m]','FontSize',fontsize)
%         ylabel('b [m]','FontSize',fontsize)
%         zlabel('Txy [MPa]','FontSize',fontsize)
%         set(gca, 'DataAspectRatio', [diff(get(gca, 'XLim')) diff(get(gca, 'XLim')) diff(get(gca, 'ZLim'))],'FontSize',fontsize)
%         print '-depsc' 'c:/Users/Tomáš Valenta/Škola/Vejška/_bachelory thesis/zprava/obrazky/12_smykove_napeti_Txy.eps'
%         
%         figure('Position',[300 200 700 650])
%         surf(X,Y,Txz,'EdgeColor','interp')
%         title('Txz','FontSize',fontsize)
%         xlabel('a [m]','FontSize',fontsize)
%         ylabel('b [m]','FontSize',fontsize)
%         zlabel('Txz [MPa]','FontSize',fontsize)
%         set(gca, 'DataAspectRatio', [diff(get(gca, 'XLim')) diff(get(gca, 'XLim')) diff(get(gca, 'ZLim'))],'FontSize',fontsize)
%         print '-depsc' 'c:/Users/Tomáš Valenta/Škola/Vejška/_bachelory thesis/zprava/obrazky/12_smykove_napeti_Txz.eps'
%     end
 end

% save Tmax.txt Tmax -ascii % uklada v [Mpa]
% save Txy.txt Txy -ascii % uklada v [Mpa]
% save Txz.txt Txz -ascii % uklada v [Mpa]

end

