Airstandard Cycle on PV diagram and calculating Efficiency

close all
clear all
clc

p1=101325;
t1=500;
gamma=1.4;
t3=2300;

bore=0.1;
stroke=0.1;
con_rod=0.15;
cr=12;

v_swept=(pi/4)*bore^2*stroke;
v_clearance=v_swept/(cr-1);

v1=v_swept+v_clearance;
v2=v_clearance;

p2=p1*cr^gamma;
t2=p2*v2*t1/(p1*v1);

constant_c=p1*v1^gamma;
v_compression=engine_kinematics (bore,stroke,con_rod,cr,180,0);

p_compression=constant_c./v_compression.^gamma;

v3=v2;
p3=p2*t3/t2;

constant_c=p3*v3^gamma;
v_expansion=engine_kinematics (bore,stroke,con_rod,cr,0,180);

p_expansion=constant_c./v_expansion.^gamma;

v4=v1;
p4=p3*(v3/v4)^gamma;

effi1=cr^(gamma-1);
effi2=1/(effi1);
E=1-effi2

figure(1);
hold on;

plot(v1,p1,\'*\',\'color\',\'r\');
plot(v_compression,p_compression);

plot(v2,p2,\'*\',\'color\',\'r\');
plot([v2 v3],[p2 p3],\'linewidth\',1,\'color\',\'r\');

plot(v3,p3,\'*\',\'color\',\'r\');
plot(v_expansion,p_expansion);

plot(v4,p4,\'*\',\'color\',\'r\');
plot([v4 v1],[p4 p1],\'linewidth\',1,\'color\',\'r\');

axis([0 8.9e-004 0 5699000])

hold off;