1 files changed, 29 insertions, 0 deletions
diff --git a/1955/CH5/EX5.16/example16.sce b/1955/CH5/EX5.16/example16.sce
new file mode 100755
index 000000000..066571015
--- /dev/null
+++ b/1955/CH5/EX5.16/example16.sce
@@ -0,0 +1,29 @@
+clc
+clear
+//input data
+b1m=46//Rotor blade angle at entry at mean section in degree
+b2m=75//Rotor blade angle at exit at mean section in degree
+a1m=75//Nozzle angle at exit at mean section in degree
+DhDt=0.6//Hub to tip ratio
+N=7500//Mean rotor speed in rpm
+Dh=0.45//Hub diameter in m
+
+//calculations
+R=0.5//Degree of reaction as a1m=b2m
+a2m=b1m//Stator angle at exit at mean section in degree
+Dm=(Dh+(Dh/DhDt))/2//Mean diameter of turbine at mean section in m
+Um=(3.1415*DhDt*N)/60//Mean blade speed in m/s
+Ca=Um/(tand(a1m)-tand(b1m))//Axial velocity in m/s
+pi=Ca/Um//Flow coefficient
+psil=pi*(tand(b1m)+tand(b2m))//Blade loading coefficient
+a1h=atand(tand(a1m)*((Dm/2)/(Dh/2)))//Nozzle angle at inlet at root section in degree
+Uh=(3.14*Dh*N)/60//Blade speed at root section in m/s
+b1h=atand(tand(a1h)-(Uh/Ca))//Rotor blade angle at entry at root section in degree
+a2h=atand(tand(a2m)*((Dm/2)/(Dh/2)))//Stator angle at exit at root section in degree
+b2h=atand((Uh/Ca)+tand(a2h))//Rotor blade angle at exit at root section in degree
+pih=Ca/Uh//Flow coefficient at root section
+Rh=(pih/2)*(tand(b2h)-tand(b1h))//Degree of reaction at root section
+psilh=pih*(tand(b1h)+tand(b2h))//Blade loading coefficient at root section
+
+//output
+printf('Mean section\n    (a)Degree of reaction is %3.1f\n    (b)Blade loading coefficient is %3.2f\nRoot section    (a)Degree of reaction is %3.2f\n    (b)Blade loading coefficient is %3.2f',R,psil,Rh,psilh)