print " D : ", AsUnit(Prop.D, 'in') print " Pitch : ", AsUnit(Prop.Pitch, 'in') Vmax = 50 h = 0 * FT N = npy.linspace(1000, 6800, 6) * 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 = 2414 * RPM print print "Static Thrust : ", AsUnit(Prop.T(N, 0 * FT / SEC, h), 'lbf') print "Measured Thrust : ", AsUnit(max(npy.array(Prop.ThrustData)[:, 1]), 'lbf') 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()
print " Pitch : ", AsUnit(Prop.Pitch, 'in') Vmax = 50 h = 0 * FT N = npy.linspace(1000, 6800, 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 = 6024 * RPM print 'Static Thrust : ', AsUnit(Prop.T(N, 0 * FT / SEC, h), 'lbf') print 'Measured Thrust : ', AsUnit(max(npy.array(Prop.ThrustData)[:, 1]), 'lbf') N = 6410 * 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()