This Page Blank
This Page Blank
This Page Blank
This Page Blank
This Page Blank
This Page Blank
This Page Blank
This Page Blank
This Page Blank
This Page Blank
This Page Blank
This Page Blank
This Page Blank
This Page Blank
/* -------------------------------------------------------------------PROGRAM: index.html Ver: 1.0 Rev: 03/01/2010 DESCRIPTION: www.neatinfo.com main menu BY: Jan Zumwalt - www.zoomaviation.com -------------------------------------------------------------------COMMENTS: Practical calculation of aircraft performance Compiled and ran on the free Pellec C compiler http://www.smorgasbordet.com/pellesc/ -------------------------------------------------------------------Ver info: V1.0 users will note slight variations in output compared to the basic version of this program due to different round off error in math packages. */ #include #include /* -------------------------------------------------------------------This section is user variables that can be customized to a particular aircraft. See The book for descriptions. -------------------------------------------------------------------*/ const float altitude_ft const float air_density_slug const float pi const float vel_delta const float vel_stall_clean_mph const float cl_max_clean const float cl_max_flap const float gross_lb const float useful_load_lb const float plane_efficiancy const float bhp const float vel_max_mph const float prop_dia_in const float prop_dia_ft const float wing_span_ft const float prop_max_rpm // end of user editable custom
= 0.00; = 0.00237; = 3.14159; = 1.00; = 67.00; = 1.53; = 2.10; = 1500.00; = 600.00; = 0.744; = 150.00; = 180.00; = 72.00; = 72 / 12; = 20.83; = 2700.00; variables
void main() { float wing_load_lb_ft = float vel_stall_flaps_mph = float wing_area_ft = float wing_aspect = float wing_chord_ft = float wing_span_effective = float wing_chord_effective = float wing_load_effective = float drag_area_ft = // float cd_drag = float vel_sink_min_ft = sqrt(sqrt(drag_area_ft)); // float pwr_min_req_hp = sqrt(wing_load_effective); // float rate_sink_min_ft = float ld_max =
// // // // // // // // // // // // // // // //
Defines the value of Pi as fixed (sealevel) Defines the value of Pi as fixed airspeed increment for each iteration VS1
Defines the value of Pi as fixed
cl_max_clean * pow(vel_stall_clean _mph,2) / 391; sqrt(wing_load_lb_f t * 391 / cl_max_flap); gross_lb / wing_load_lb_ft; pow(wing_span_ft,2) / wing_area_ft; wing_span_ft / wing_aspect; wing_span_ft * sqrt(plane_effician cy); wing_area_ft / wing_span_effective ; gross_lb / wing_span_effective; .8 * bhp * 146625 / pow(vel_max_mph,3);
// // VS0 // // // // // // //
drag_area_ft / wing_area_ft; 11.29 * sqrt(wing_load_effect ive) /
//
.03922 * sqrt( sqrt(drag_area_ft)) * wing_load_effective * 33000 * pwr_min_req_hp / gross_lb; .8862 * wing_span_effective / sqrt(drag_area_ft);
// //
float float float float float float
drag_min cl_min_sink rate_climb_ideal prop_tip_mach prop_vel_ref static_thrust_ideal
printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t printf("\n\t
= = = = = =
gross_lb / ld_max; 3.07 * sqrt(drag_area_ft) / wing_chord_effective; 33000 * bhp / gross_lb; prop_max_rpm * prop_dia_ft * .05236 / 1100; 41.9 * pow(bhp / pow(prop_dia_ft,2),.3 3333); 10.41 * pow(bhp * prop_dia_ft,.66666);
wing_load_lb_ft vel_stall_flaps_mph wing_area_ft wing_aspect wing_chord_ft wing_span_effective wing_chord_effective wing_load_effective drag_area_ft cd_drag vel_sink_min_ft pwr_min_req_hp rate_sink_min_ft ld_max drag_min cl_min_sink rate_climb_ideal prop_tip_mach prop_vel_ref static_thrust_ideal
= = = = = = = = = = = = = = = = = = = =
%.02f", %.02f", %.02f", %.02f", %.02f", %.02f", %.02f", %.02f", %.02f", %.04f", %.02f", %.02f", %.02f", %.02f", %.02f", %.02f", %.02f", %.02f", %.02f", %.02f",
// // // // // //
wing_load_lb_ft); vel_stall_flaps_mph); wing_area_ft); wing_aspect); wing_chord_ft); wing_span_effective); wing_chord_effective); wing_load_effective); drag_area_ft); cd_drag); vel_sink_min_ft); pwr_min_req_hp); rate_sink_min_ft); ld_max); drag_min); cl_min_sink); rate_climb_ideal); prop_tip_mach); prop_vel_ref); static_thrust_ideal);
printf("\n\n"); printf("\n\t ----------------------------------------------------------------------------"); printf("\n\t airspeed \t climb rate \t prop eff \t sink rate \t rennolds num"); printf("\n\t v(mph) \t rc(fpm) \t eta \t rs(fpm) \t re=rho*v*c/mu"); printf("\n\t ----------------------------------------------------------------------------");
float eta float fp float rc float rc1 float rc2 float rcmax float rec float rsh float rmu float rs float sig // float t float t1 float t2 float v float vh float vmax float vt float wv2 while (rc { vh rsh rs vt t2 eta rc
= = = = = = = = = = = = = = = = = =
1; 0; 1; 0; 0; 0; 0; 0; 1; 0; pow(1 - altitude_ft / 145800,4.265); = 518.7 - 0.00356 * altitude_ft; .3333; 0; vel_stall_clean_mph ; 0; 0; 0; 0;
> 0) = = = = = = =
v / vel_sink_min_ft; .25 * (pow(vh,4) + 3) / vh; rsh * rate_sink_min_ft; v / prop_vel_ref; sqrt(1 + .23271 * pow(vt,3)); .92264 * vt * (pow( 1 + t2,t1) - pow(t2 - 1,t1)) * .85; rate_climb_ideal * eta - rs;
rc2 = rc; rec = sig * v * wing_chord_ft * 9324 / rmu; if (rc < 0) break; rcmax = fmax(rc,rcmax); vmax = fmax(v,vmax); printf("\n\t %.01f \t %.01f \t %.02f rec ); v = v + vel_delta * rc2 / (rc2 - rc1); } fp
\t
%.01f
\t
%.0f",v, rc, eta, rs,
= rcmax * useful_load_lb / 33000 / bhp * (1 - (vel_stall_flaps_mph / vmax));
wv2
= gross_lb * pow(v,2);
printf("\n\n\t performance parameter......... fp printf("\n\t kinetic energy parameter...... wv2 printf("\n\t maximum rate of climb.. ...... rcmax printf("\n\t maximum speed............... .. vmax printf("\n\t useful load lb....... useful_load_lb
= = = =
= %.04f",fp); %.02f",wv2); %.02f",rcmax); %.02f",vmax); %.02f",useful_load_ lb);
printf("\n\n\t +------------------------------------------------------+"); printf("\n\t | Thank you for using |"); printf("\n\t | Air-Performance 1.0 |"); printf("\n\t +------------------------------------------------------+"); printf("\n\n\t Press key to exit... "); while ((getchar()) != '\n'); printf("\n"); }
OUTPUT... ----------------------------------------------------------------------------airspeed climb rate prop eff sink rate rennolds r ennolds num v(mph) rc(fpm) eta rs(fpm) re=rho*v*c/mu ----------------------------------------------------------------------------67.0 1175.5 0.63 896.1 2561031 68.0 1196.3 0.63 891.1 2599255 69.0 1216.2 0.64 886.6 2637480 70.0 1235.3 0.64 882.6 2675704 71.0 1253.5 0.65 879.1 2713928 72.0 1271.0 0.65 876.1 2752153 170.0 183.4 0.82 2521.4 6498138 171.0 147.3 0.82 2559.0 6536362 172.0 110.7 0.82 2597.2 6574587 173.0 73.5 0.82 2635.9 6612811 174.0 35.8 0.82 2675.0 6651035 performance parameter......... fp kinetic energy parameter...... wv2 maximum rate of climb.. ...... rcmax maximum speed................. vmax useful load lb....... useful_load_lb
= = = = =
0.1206 45937500.00 1482.42 174.00 600.00
+------------------------------------------------------+ | Thank you for using | | Air-Performance 1.0 | +------------------------------------------------------+ Press key to exit...
This Page Blank
Appendix A (information not in the original book)
Appendix B (information not in the original book)
Aviation Math Symbols
Α, α Alpha
angle
Β, β Γ , γ ∆, δ Ε , ε Ζ , ζ Η , η Θ, θ Ι, ι Κ, κ Λ, λ Μ, µ Ν, ν Ξ, ξ Ο, ο Π, π Ρ , ρ Σ, σ Τ, τ Υ, υ Φ, φ Χ, χ Ψ, ψ Ω, ω ∠ ° ≅ ≤ ≥ ± ∞ ∑ o
C F o R o K o
Beta Gamma Delta
change or press/ratio
Epsilon Zeta Eta Theta
temp/ratio
Iota Kappa Lambda Mu Nu Xi Omicron Pi
3.141
Rho
density
Sigma
density/ratio
Tau Upsilon Phi
trig angle, i.e. sin, cos
Chi Psi Omega
angle degree approximate less than or equal more than or equal plus or minus infinity sum = (oF-32) * 5/9 = (oC * 5/9) + 32 = oF + 460 = oC + 273
ft/min = mph * 88 ft/sec = kt * 1.69 ft/sec = mph * 1.47 kt = mph * 0.87 fps
kt
=
mph
= kt * 1.15
1.69
Note: unit of measure may be in subscript. For example, a distance (X) may be given as Xft or Xin. a = AC = AR = b = C = CF CG = CL CLmax
acceleration aerodynamic center wing aspect ratio (no dim) wing span chord = coefficient of force (no dim) center of gravity = coefficient of lift (no dim) = max CL (no dim)
d = D = Dmin = Di = Dp = Dt = EAS = F = Fb = FR = ft = ft/m ft/s G = H = h = hp = IAS = k = KE = kt = L = L/D = L/Dmax m = MAC = mph = N = = ŋ = ŋP P = PA = PE = PR = q = RN = ROC = S = sl = SL = = Ŧ t = = ŦA TAS = TE = = ŦR u = V = VS = VX = VY = W = X = Y =
δ =
θ = σ =
distance drag minimum drag induced drag parasite drag total drag equivalent air speed force braking force friction feet = feet per minute = feet per second gravity – 32.2(ft/s2) head (total pressure) height horse power indicated air speed constant kinetic energy knot lift lift to drag ratio (no dim) = max lift/drag ratio (no dim) mass mean aerodynamic chord mile per hour weight on wheels efficiency propulsive efficiency air density(slugs) power available potential energy power required dynamic pressure renold number rate of climb surface area slug sea level thrust (lb) time thrust available true air speed total energy thrust required friction coefficient (no dim) velocity velocity at stall velocity best angle velocity best ROC weight distance or unknown height or unknown
P( hg ) PSL ( hg )
T ° SL T ° δ δ SL
* θ =
δ δ SL
*
T ° SL T °
DT(lb) = ½*p(alt)*V2(fps)*S(ft2)*Cd 2
σ ∗ V ( kts )
q =
295
F (lb)
a(fps2) =
=
Di min = drag induced at (L/D)max(lb)
=
D min
Di =
D min
2
W (lb) / 32.2( g )
F (lb) m( sg ) 2
=
DP =
2
V start ( fps ) + V end ( fps )
2
AR = AR =
span
span 2
=
chord ( mean )
area
( L / D ) max
F b = u * N
hp =
T (lb) * V ( kt ) 325
L (lb) 2
2
1 / 2 * p ( sg ) * V ( fps ) * S ( ft )
W (lb)
(sg) (sg) = m CL =
32( fps )
295 ∗ w(lb) 2
2
EAS = CAS ± ∆V”chart
CAS = IAS ± ∆V”chart”
CL =
V * 1 V 2
W (lb)
D min =
* [T (lb) − D (lb) − F R (lb) ]
W (lb)
2
Dt = Dp + Di
2 * d ( ft ) 32.2( g )
V * 2 V 1
2
σ ∗ V ( kt ) ∗ S ( ft )
2 ∗ W (lb)
V s(fps) =
2 C L max * p (alt ) * S ( ft )
d(ft) = Vav(fps) * t(s)
=
V start ( fps ) + V end ( fps ) 2
[V
2
σ d 1 _ takeoff ( ft ) * 1 σ 2
dland(ft) =
d2_land(ft) =
V tip tip(fps) =
W (lb) d 1 _ takeoff ( ft ) * 2 W 1 (lb)
[V
r ( ft ) * rpm
2
9.55 2
2 * a ( fps ) * d ( ft )
V DIST =
2
( fps ) ± wind ( fps ) ]
2
Takeoff
295 ∗ W (lb)
V Dmin Dmin =
2
C L * σ * S ( ft )
2
2 * a ( fps )
W (lb) d 1 _ Land ( ft ) * 2 W 1 (lb)
wind ( fps )
V 1 _ land ( fps )
=
d 1 _ Land ( ft ) * 1 +
=
d 1 _ land ( ft ) *
V IAS
( fps ) ± wind ( fps )] 2 * a( fps )
d2_takeoff(ft) =
=
TAS = EAS * density = 2
Takeoff
dtakeoff(ft) =
* t(s)
σ 1 σ 2
2
S(ft2) =
t(s) =
W 2 (lb) * L2 ( ft ) W 1 (lb) + W 2 (lb)
V ( fps ) 2 a ( fps )
1 σ
A1*V1
=
A2*V2
q1
=
2 ½ * P * V 1
H
=
=
A3*V3 2
2
alt + bas
hyp =
=
base cos ∠a
sin ∠a
=
2
V2(kt) =
A2 ( ft 2 ) 2
q2 =
=
σ ∗ V 2 ( kts)
295
p2 = H - q2
bas tan ∠a
A3(ft ) =
= hyp * Sin ∠α
bas = hyp * Cos ∠α
=
alt
A2 ( ft 2 ) ∗ V 2 (kt ) V 3 (kt ) 2
tan ∠b
q3 =
σ ∗ V 3 (kts)
= alt * tan ∠b = hyp * cos ∠a
295
p3 = H – q3
= hyp * sine ∠b
Tan ∠α ∠α =
295
A1 ( ft 2 ) ∗ V 1 ( kt )
2
tan ∠b
σ ∗ V 1 (kts)
P1 + q1
cos ∠b
sin ∠b
hyp + bas
bas
=
alt
base
2
alt =
=
=
2
alt
=
alt bas
∠a = 90 – b ***************************************
=
CG(ft)
=
A1 = V1 =
P1 =
A2 = V2 =
P2 =
A3 =
V3 =
σ = σ = P3 =
q1 =
∑ [W (lb) * Arm( ft )] ∑ [W (lb)]
************************************
q2 =
σ =
∑ Mom( ft / lb) ∑W (lb)
q3 =
Appendix C (information not in the original book)
Design Notes
Appendix D (information not in the original book)
Graphs & Charts
Log Paper