rwesfw

clearvars;
clc;
close all

%% fasta v?rden
r0=0.025;   %Rod radius
r1=0.375;   %Wheel radius
g = 9.81;
h = linspace(0,1,1000);
frw = 1.5;
ns=3;
Mb=1500; %boat weight
Mm=95;    %Counter weight mass
eta = 0.88;

Fh=50;   %Handforce on wheel


Fb=Mb*9.81;   %force from boat
Fm=Mm*9.81;     %force from counterweight

%% Krafter i konstruktionen

s = 0.317;          %from cad
h0 = 1.807;         %from cad
fi = atand(s./(h0-h));
Fboat = Fb/4;
Frope = Fboat./(cosd(fi));
Fx = Fboat.*tand(fi);
Fres = sqrt((Fx+Frope).^2+Fboat.^2);
Frear = frw*Fres;
Frearmax = Frear(end);
Frepb=Frope(end);      %Force in a rope connected to boat
Frepm=Fm;         %Force in a rope connected to CW

Tb=Frope(end)*r0*4;       %Total torque on rod from boat
Tm=8*Fm*r0;     %Total torque from CW on rod
Tut=(Tb-Tm);      %Torque from hand on rod
Tin=Fh*r1;      %Torque in wheel

ireq=Tut/(Tin*eta);       %needed transmission
itot=ireq;        %actual transmission

%% Skillnaden i vikt n?r b?ten ?r i vattnet mot i luften
% mboatvek1=[0:15:1500];
% mboatvek2(1:200)=1500;
% mboatvek=[mboatvek1,mboatvek2];
% hboatvek=[0:0.00333:1];

%plot(hboatvek,mboatvek);

%% Tid f?r upphissning

nin=60;  %rpm
Power = Tin*nin*pi/30;

varvlift=h(end)/(2*pi*r0);   %Varv st?ngen beh?ver snurra

varvratt=varvlift*itot; %Varv ratten beh?ver snurra

timelift=varvratt/nin; %Tid f?r upphissning

%% Ber?kning av krafter p? st?ngen


sigmamax = 350e6;
sigmatill = sigmamax/3;
ri = r0 - 0.5*(16*Tin/(pi*sigmamax))^(1/3);


%% kraften p? ett lager z-led

Flager=Fm*2/sqrt(2);

%% key lengths

Win=0.008;
Sy=350e6;
D=0.03;
H=0.007;

Lt=4*Tin*ns/(D*Win*Sy);
Lp=4*Tin*ns/(D*H*Sy);

Ltr=4*Tut*ns/(D*Win*Sy);
Lpr=4*Tut*ns/(D*H*Sy);


%% bearing calculations

p=3; %single row deep grove ball bearing
nut=nin/itot;
Prod=Flager;

L10h = 15*4/60*52;

Crod=Prod*(L10h*60*nut/10^6)^(1/p);

Pin=Fh;

Cin=Pin*(L10h*60*nin/10^6)^(1/p);

%% Beam parameters

Lh = 2.55;
Lv = 3.4;
E = 190e9;
n = 3;

%% Beam dimensions

Ih = 2*frw*Lh*Frepb/(E*pi^2*n);
Iv = Lv*frw*Frepb/(2.05*E*pi^2*n);
Ihc = 8433;
Ivc = 8433;
hsafe = Ihc/(Ih*1e9);
vsafe = Ivc/(Iv*1e9);

disp(['Safety factor horizontal beam ', num2str(hsafe)]);
disp(['Safety factor vertical beam ', num2str(vsafe)]);

%% Pulley, boat

s = 0.317;          %from cad
h0 = 1.807;         %from cad
fi = atand(s./(h0-h));
Fboat = Frepb;
Frope = Fboat./(cosd(fi));
Fx = Fboat.*tand(fi);
Fres = sqrt((Fx+Frope).^2+Fboat.^2);
Frear = frw*Fres;
Frearmax = Frear(end);

taumax = 350e6;
d = sqrt((16*Frearmax*n)/(3*pi*taumax));
disp(['Minimum pulley diameter ', num2str(d*1000),' mm']);
d = 0.015;
disp(['Chosen pulley diameter ', num2str(d*1000), ' mm']);
np = r1*nut/(d/2);
Cp = sqrt(2)*Frearmax*(L10h*60*np/10^6)^(1/p);

figure('Name','fi and Fres as functon of height');
subplot(2,1,1)
plot(h,fi);
xlabel('Height [m]');
ylabel('\phi [°]');
title('\phi (h) ')
grid on
subplot(2,1,2)
plot(h,Fres/1000);
xlabel('Height [m]');
ylabel('F_{res} [kN]');
title('F_{res} (h)')
grid on

%% Energy

alphamax = varvratt*2*pi;
M = linspace(Tin,Tin,1000);
alpha = linspace(0,alphamax,1000);
Ehand = trapz(alpha,M);
Ecw = 8*Mm*g*(h(end)-h(1));
Ein = Ehand + Ecw;
Losses = Ehand*(1-eta);

Eut = Mb*g*(h(end)-h(1)) + Losses;

diff = Ein - Eut;

disp(['Energy check: ', num2str(diff), ' Joule']);
disp(['Lifting time: ', num2str(timelift),' min']);