% 
%   Formule per la lubrificazione del cusncietto completo lubrificato (full Sommerfeld)
%
close all
clear all
%
%
%
ep  = 0:0.01:0.999;
for j=1:1:length(ep)
A(j) = 0.12e2 * ep(j) * pi / (ep(j) ^ 2 + 2) * ((1 - ep(j) ^ 2) ^ (-0.1e1 / 0.2e1));
B(j) = 0.4e1 * pi * (2 * ep(j) ^ 2 + 1) / (ep(j) ^ 2 + 2) * ((1 - ep(j) ^ 2) ^ (-0.1e1 / 0.2e1));
k(j) = 0.2e1 / 0.3e1 * (ep(j)^ (-0.1e1 / 0.2e1)) * sqrt(pi) * (2 * ep(j) ^ 2 + 1) * sqrt(0.3e1) * ((ep(j) ^ 2 + 2) ^ (-0.1e1 / 0.2e1)) * ((1 - ep(j) ^ 2) ^ (-0.1e1 / 0.4e1));
end;
%
scala_A = 0.5;
scala_B=0.5;
scala_k = 2;
%
figure(1)
set(1,'Color','w');
hold on
plot(ep,A*scala_A,'-r','LineWidth',1);
hold on
plot(ep,B*scala_B,'-g','LineWidth',1);
hold on
plot(ep,k*scala_k,'-.b','LineWidth',1);
%
title(" Full Sommerfeld - Functions ");
ylabel("  A(epsilon), B(epsilon), and k(epsilon) ");
xlabel("  x axis  ");
legend ("A(epsilon)","B(epsilon)","k(epsilon)", "location","outside");
grid "minor" on;
hold off
%
%
% 
%  Dati numerici in input
%
mu  = 0.07;
R = 0.0375;
epsilon  = 0.7;
delta = 0.07;
V = 1.176;
%
%
ho= -2 / (-1 + epsilon ^ 2) * delta / (epsilon ^ 2 + 2) * (1 - 2 * epsilon ^ 2 + epsilon ^ 4);
%
theta = 0:0.1:2*pi;
zeri = theta * 0;
h = delta*(1+epsilon*cos(theta));
for j = 1:length(theta)
pp(j) =  0.6e1 * mu * V * R * epsilon * ((3 * epsilon) + (epsilon ^ 2)* cos(theta(j)) + 0.2e1* cos(theta(j))) / delta ^ 2 / (epsilon ^ 2 + 2) / (0.1e1 + epsilon* cos(theta(j)))^ 3;
end;
for j=1:length(theta)
p(j) = 0.6e1 * mu * V * R * epsilon / delta ^ 2 / (epsilon ^ 2 + 2) * sin(theta(j)) * (epsilon * cos(theta(j)) + 0.2e1) / ((epsilon ^ 2) * cos(theta(j)) ^ 2 + 0.2e1 * epsilon * cos(theta(j)) + 0.1e1);
end;
%
%
scala_h = 10;
scala_pp = 1;
scala_p = 2;
%
figure(2)
set(2,'Color','w');
plot(theta,zeri,'-y','LineWidth',1);
hold on
plot(theta,h*scala_h,'-r','LineWidth',1);
hold on
plot(theta,pp*scala_pp,'-g','LineWidth',1);
hold on
plot(theta,p*scala_p,'-.b','LineWidth',1);
%
title(" Full Sommerfeld ");
ylabel("  h(x), pp(x), and p(x) ");
xlabel("  x axis  ");
legend ("Riferimento","Altezza","Gradienti","Pressioni (per u.d.l.)", "location","outside");
grid "minor" on;
hold off
%
%
figure(3)
%
scala_cerchio = 10;
scala_p = scala_p * 0.4;
scala_pp = scala_pp * 0.3;
cerchio = scala_cerchio *ones(1,length(theta));
polar(theta,cerchio,'-y');
hold on
polar(theta,(cerchio + pp*scala_pp),'-g');
hold on
polar(-theta,(cerchio + p*scala_p),'-.b');
axis("equal");
title(" Full Sommerfeld ");
ylabel("  pp(theta), and p(theta) ");
xlabel("  x axis  ");
legend ("Riferimento","Gradienti","Pressioni (per u.d.l.)", "location","outside");
grid "minor" on;