blob: 9bdd8718c303fc4c6f84776478d07539c152856b (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
pathname=get_absolute_file_path('9_01.sce')
filename=pathname+filesep()+'9_01data.sci'
exec(filename)
x=poly(0,'x');
P=x-10*x+9.5;
t=roots(P);
V2=%pi*b^2*(Stroke+t)*(10^-6)/4;disp(V2,"V2=%pi*b^2*(Stroke+t)/4");
V3=V2/r;disp(V3,"V3=V2/r");
V5=V2;V4=V3;
Wcomp=(P2*V2-P3*V3)/(1-y);
disp(Wcomp,"Wcomp=","Wcomp=P2*V2-P3*V3/(1-y);","work done in compression cycle Wcomp:")
Wpower=(P5*V5-P4*V4)/(1-y);
disp(Wpower,"Wpower=","Wpower=P5*V5-P4*V4/(1-y);","work done in power stroke Wpower:")
Pa=6*n*nm*(rpm)*(Wpower-Wcomp)/120;
disp(Pa,"Pa=n*nm*(rpm)*(Wpower-Wcomp)/120","power available Pa:")
printf("\Answer:\n")
printf("\n\Power available from the engine propeller combination: %f J/s \n\n",Pa)
|
