TribMike

format long
clc
clear all

%%  Inlämningsuppgift Tribologi
% Givna konstanter, samband och värden
M0 = 6;
n0 = 1800;
k = 119;
R = 0.025;
alfaN = 20; % grader
beta = 18; % grader
sigmaTN = 2.9E8;

nmin = 0.25*n0;

n = nmin;
M = M0*(1-n/n0);
F = k*M;

Fh = M/R;
Fv = M*tand(alfaN)/(R*cosd(beta));
Fa = M*tand(beta)/R;

T = 80;
etac = 0.2;

%% Uppgift 3
% Gissade värden på bredd
bA = 0.01;
bB = 0.01;

% Krafter
FBa = F - Fa
FBv = (Fv*(0.03/2 + bA/2) - Fa*R)/(0.03 + bA/2 + bB/2);
FAv = Fv - FBv;
FAh = Fh*(0.03/2 + bA/2)/(0.03 + bA/2 + bB/2);
FBh = Fh - FAh;

FAr = sqrt(FAv^2 + FAh^2)
FBr = sqrt(FBv^2 + FBh^2)

% Gissade värden på spänningskoncentrationsfaktorer
alfaAv = 2;
alfaAb = 2;
alfaBb = 2;

% Lagerdiameter

dA = (1/(pi*sigmaTN)^2*((16*alfaAb*FAr*bA)^2 + 3*(alfaAv*16*M)^2))^(1/6)

dB = (16*alfaBb*FBr*bB/(pi*sigmaTN))^(1/3)


%% Uppgift 4
L10h = 20E3;

CA = ((L10h*60*nmin)/1E6)^(1/3) * FAr

CB = ((L10h*60*nmin)/1E6)^(1/3) * (0.56*FBr + 1*FBa)


%% Uppgift 5
f0 = 13;
C0 = 2.36E3;

f = f0*FBa/C0

% Linjär interpolation
f1 = 2.07;
f2 = 3.45;
Y1 = 1.31;
Y2 = 1.15;

Y = Y1 + (f-f1)*(Y2-Y1)/(f2-f1);

% Nya värden på spänningskoncentrationsfaktorer och bredd ger nya d och C
bAny = 0.0035;
bBny = 0.009;
alfaAvny = 1.8;
alfaAbny = 2.9;
alfaBbny = 1.9;

% Nya krafter pga nya alfa och b
FBvny = (Fv*(0.03/2 + bAny/2) - Fa*R)/(0.03 + bAny/2 + bBny/2);
FAvny = Fv - FBvny;
FAhny = Fh*(0.03/2 + bAny/2)/(0.03 + bAny/2 + bBny/2);
FBhny = Fh - FAhny;

FArny = sqrt(FAvny^2 + FAhny^2);
FBrny = sqrt(FBvny^2 + FBhny^2);


dAny = (1/(pi*sigmaTN)^2*((16*alfaAbny*FArny*bAnyny)^2 + 3*(alfaAvny*16*M)^2))^(1/6)
CAny = ((L10h*60*nmin)/1E6)^(1/3) * FArny

dBny = (16*alfaBbny*FBrny*bBny/(pi*sigmaTN))^(1/3)
CBny = ((L10h*60*nmin)/1E6)^(1/3) * (0.56*FBrny + Y*FBa)


% Testa igen med 619/7 och *6200

% Nya värden på spänningskoncentrationsfaktorer och bredd
bAnyny = 0.005;
bBny = 0.009;
alfaAvnyny = 1.7;
alfaAbnyny = 2.1;
alfaBbny = 1.9;

% Nya krafter pga nya alfa och b
FBvny = (Fv*(0.03/2 + bAnyny/2) - Fa*R)/(0.03 + bAnyny/2 + bBny/2);
FAvny = Fv - FBvny;
FAhny = Fh*(0.03/2 + bAnyny/2)/(0.03 + bAnyny/2 + bBny/2);
FBhny = Fh - FAhny;

FArny = sqrt(FAvny^2 + FAhny^2);
FBrny = sqrt(FBvny^2 + FBhny^2);

dAverk = 0.007;
dBverk = 0.010;

sigmaAjmf = sqrt((alfaAbnyny*32*FArny*bAnyny/2/(pi*dAverk^3))^2 + 3*(alfaAvnyny*16*M/(pi*dAverk^3))^2)
sigmaBjmf = (alfaBbny*32*FBrny*(bBny/2))/(pi*dBverk^3)