% Pevnostní výpočet - Diplomová práce - Růžička Pavel
close all; clear all; clc

% Výpočet sil dFz a dFx pro stanovení sil v Abaqusu


%VSTUPNÍ HODNOTY
Dp=12;                                          %prumer rotoru v palcich
rpm=8197;                                       %otacky za min
rho=1.225;                                      %hustota vzduchu
D=Dp*0.0254;                                    %prumer rotoru [m]
A=pi*(D^2)/4;                                   %obsah disku [m^2]
Vtip=pi*D*rpm/60;                               %rychlost na konci listu [m/s]
dr=0.08; 
r=0.2:dr:1;
y=(D/2)*r;
dy=dr*D;

c = dlmread('c_pevnost.txt','');
Cl = dlmread('Cl_pevnost.txt','');
Cd = dlmread('Cd_pevnost.txt','');
fir= dlmread('fir_pevnost.txt','');
c=c';
Cl=Cl';
Cd=Cd';
fir=fir';

fi=fir*180/pi();
omega=2*pi*rpm/60;
for i=1:length(y)
    Ut(i)=omega*y(i);
end

U=Ut;
for i=1:length(y)
    dL(i)=0.5*rho*((U(i))^2)*c(i)*Cl(i)*dy;
    dD(i)=0.5*rho*((U(i))^2)*c(i)*Cd(i)*dy;
end


for i=1:length(y)
    dFz(i)=dL(i)*cos(fir(i))-dD(i)*sin(fir(i));
    dFx(i)=dL(i)*sin(fir(i))+dD(i)*cos(fir(i));
end

dFz_list=dFz/2;
dFx_list=dFx/2;

Fz_list=sum(dFz_list);
Fx_list=sum(dFx_list);

y_mm=y*1000;
y_mm=round(y_mm,2);
pevnost_abaqus=[y_mm(:), dFz_list(:), dFx_list(:)];


%ZÁPIS DO EXCELU
    Excel_pevnost_abaqus = array2table(pevnost_abaqus, 'VariableNames', {'y [mm]', 'Fz [N]', 'Fx [N]'});

  % Uložení tabulky do Excelu
    writetable(Excel_pevnost_abaqus, 'Pevnost_data_abaqus_navrh_ot.xlsx');