otto cycle pv diagram

Objective: % function function compressionvolume = otto_cycle_function(stroke,bore,cond_rod,cr,v_clearance,start_angle,end_angle) a=stroke/2; R=cond_rod/a; teta = linspace(start_angle,end_angle,180); term1=0.5*(cr-1);term2=R+1-cosd(teta);term3=(R^2-sind(teta).^2).^0.5;compressionvolume = v_clearance.*(1+(term2-term3).*term1);end…