h = 0 * FT
N = npy.linspace(1000, 12300, 5) * RPM
Alpha = npy.linspace(-25, 25, 41) * ARCDEG
V = npy.linspace(0, Vmax, 30) * FT / SEC
Prop.CoefPlot(Alpha, fig=1)
Prop.PTPlot(N, V, h, 'V', fig=2)
#
# N = npy.linspace(0, 13000,31)*RPM
# V = npy.linspace(0,Vmax,5)*FT/SEC
#
# Prop.PTPlot(N,V,h,'N', fig = 3)
Prop.PlotTestData(fig=4)
N = 12240 * RPM
print "Max " + AsUnit(Prop.MaxRPM(), 'rpm', '%3.0f') + " at " + AsUnit(
Prop.MaxTipSpeed, 'ft/s') + " tip speed "
print
print "Static Thrust : ", AsUnit(Prop.T(N, 0 * FT / SEC, h), 'lbf')
print "Measured Thrust : ", AsUnit(max(npy.array(Prop.ThrustData)[:, 1]),
'lbf')
N = 12660 * RPM
print
print "Static Torque : ", AsUnit(
Prop.P(N, 0 * FT / SEC, h) / N, 'in*ozf')
print "Measured Torque : ", AsUnit(max(npy.array(Prop.TorqueData)[:, 1]),
'in*ozf')
pyl.show()
#
# These are corrected for standard day
#
# RPM, Thrust
#Prop.ThrustData = [(8100 *RPM, 4 *LBF + 8*OZF),
# (9200 *RPM, 5 *LBF + 13*OZF),
# (11200 *RPM, 9 *LBF + 3*OZF)]
#
# RPM, Torque
#Prop.TorqueData = [(11000 *RPM, 114.768*IN*OZF)]
################################################################################
if __name__ == '__main__':
print "Max " + AsUnit(Prop.MaxRPM(), 'rpm', '%3.0f') + " at " + AsUnit(Prop.MaxTipSpeed, 'ft/s') + " tip speed ",
Vmax = 100
h=0*FT
N=npy.linspace(1000, 13000,4)*RPM
Alpha = npy.linspace(-25,25,31)*ARCDEG
V = npy.linspace(0,Vmax,11)*FT/SEC
# Prop.CoefPlot(Alpha,fig = 1)
Prop.PTPlot(N,V,h,'V', fig = 2)
N = npy.linspace(0, 13000,31)*RPM
V = npy.linspace(0,Vmax,5)*FT/SEC
Prop.PTPlot(N,V,h,'N', fig = 3)