% Loops and Basic Programming
% Dr Owen Williams
% University of Wales TSD

clear;clc

% Calculate the Real Cylinder Volume

Bore=98;                %Cylinder Bore (mm)
Stroke=66;              %Piston Stroke (mm)
Rod=120;                %Con Rod Length (mm)
Clearance=58.8;         %Clearance Volume (cc)

Vc=Clearance*1000                                       %Clearance Volume (mm^3)
Rad=Stroke/2                                            %Crank Radius (mm)
Comp=((Stroke*((pi)*(Bore.^2))/4)+Vc)/Vc                %Compression Ratio
Crank=-180:1:180;                                       %Crank Angle Range
s=Rad*cosd(Crank)+((Stroke.^2)-(Rad^2)*sind(Crank).*sind(Crank)).^(1/2);  %s  
Vol=(Vc+((((pi)*(Bore.^2))/4))*(Stroke+Rad-s))/1000;    %Cylinder Volume (cc)

figure(1)
plot(Crank,Vol)

% Calculation Loop for Non-Combustion Cylinder Pressure
Pcomp=1             % Ambient Air Pressure (BAR)

for i=1:1:360;
    if (i<180)
    n=1.33;
    else 
    n=1.27;
    end
    Pcomp(i+1)=Pcomp(i)*(Vol(i)/Vol(i+1))^n;
end

figure(2)
plot(Crank,Pcomp)