Estimation of Wing Loading and Thrust Loading-2 Topics
Estimation of Wing Loading and Thrust Loading-2 Topics
Tulapurkara
Chapter-4
Chapter 4
Estimation of wing loading and thrust loading- 2
Lecture 10
Topics
4.4 Selection of wing loading based on prescribed flight speed (Vp)
4.4.1 Alternate break-up of drag polar
Example 4.2
Example 4.3
Example 4.4
4.5 Selection of wing loading based on absolute ceiling (Hmax)
Example 4.5
Hence,
1
ρ Vp2 SCD
tv = 2
W
1
Denoting ρ Vp2 as qp and (W/S) as p yields:
2
CD qp
tv = (4.13)
p
The drag coefficient (CD) is generally written as:
CD = CD0 (M) + K(M) CL2 (4.14)
where, CD0 = Parasite drag coefficient. CD0 is nearly constant for M < Mcrit and
then it becomes a function of M.
K CL2 = CDi = Lift-dependent drag coefficient. K is nearly constant for M < Mcrit, and
then it becomes a function of M.
4.4.1 Alternate break-up of drag polar
In the remark at the end of section 4.2, it is mentioned that in the approach of
Ref.1.6 the changes in drag polar, as a result of the changes in wing loading are
taken care of by an alternate representation of the drag polar. This
representation is explained in this section.
Do
where, CDO = parasite drag coefficient , Do = parasite drag
1 2
ρV S
2
1
K= , e = Oswald efficiency factor; A = Aspect ratio of wing
πAe
The expressions for CD0 and K for subsonic airplanes with jet engine, turbo-prop
engine and piston engine are given by Eqs.(3.57),(3.58) and (3.59). These
expressions include influence of parameters like wing aspect ratio and wing
sweep. These expressions are expected to give reasonable drag polars for the
present stage of the preliminary design process. The steps to obtain the alternate
break up of drag polar are as follows.
(1) Noting that the parasite drag (D0) is the sum of parasite drags of the major
components plus the effect of interference, it can be expressed as:
D0 = (D0)wing + (D0)fuselage +( D0)nacelle + (D0)ht +( D0)vt + (D0)etc + ( D0)int (4.16)
where, (D0)etc is the parasite drag of components like landing gear, external fuel
tanks and armament mounted on the airplane. (D0)int is the incremental parasite
drag due to the interference between flows past wing and fuselage, tail and
fuselage etc.
CD0 = D0 / (0.5 V2 S)
Hence,
1 2 1 1
ρV S ρ V 2 Sfuse ρ V 2 Snac
CD0 = 2 (CD0 ) w + 2 (CD0 )fuse + 2 (CDo )nac
1 2 ρ V 2
S 1 2
ρV S ρV S
2 2
1 1 1
ρ V 2 Sht ρ V 2 Svt ρ V 2 Setc
+ 2 (CDo )ht + 2 (CDo ) vt + 2 (CDo )etc + (CDo )int
1 2 1 2 1 2
ρV S ρV S ρV S
2 2 2
Note that (CDo)etc and Setc may include more than one items.
Sfuse S S S
Hence, (CDo )fuse + CDo nac nac + CDo etc etc +CDint is written as CDo other other i.e.
S S S S
Sfuse S S S
(CDo )fuse + CDo nac nac + CDo etc etc +CDint = CDo other other (4.20)
S S S S
The parasite drag coefficient can also be expressed as:
S wet
CDo = Cfe (4.21)
S
where, Cfe = equivalent skin friction drag coefficient and
Swet = wetted area of airplane
Example 4.2
Consider a high subsonic jet airplane with an initial estimate of gross weight (W0)
as 60,000 kgf and wing loading of 5500 N/m2. Estimate Cfe.
Solution :
In this case the wing area would be :
60000 × 9.81
S = W / (W/S) = = 107.02m2
5500
(4) Using the definition of Cfe in Eq.(4.21), the right hand side of Eq.(4.18) is
expressed as :
Swet w
CDo W 1+
Sht Svt
+ = K t CDo W = K t Cfe (4.22)
S S S
(Swet)w = wetted area of wing.
(Swet)w 2 Sexposedwing {1 + 1.2 (t/c)w }. The quantity Sexposed wing can be obtained as
the wing parameters and the fuselage width are roughly known from the
preliminary three view stage. This is illustrated in example 4.3.
(S wet )w
The quantity Kt Cfe is denoted by f1 i.e.
S
(S wet )w
F1 = K t Cfe (4.23)
S
Sother
(5) The quantity (CDo)other of Eq.(4.20) is expressed as :
S
Sother
CDo other = (CDo -F1 ) (4.24)
S
It may be pointed out that Sother remains same but ‘S’ changes as the wing
loading changes. Let the current value S be denoted as Scurrent and ‘S’ based on
preliminary three-view stage as Sold .
Rewriting Eq.(4.24) yields:
b b
Further, S = c r +c t = c r 1+ λ
2 2
2×107.02
Hence, c r = = 5.47m .
31.55×1.24
ct = 1.31 m.
Hence, equation for local chord (c) is
cr - c t
c = cr - y = 5.47 - 0.264 y
b/2
Half of the width of fuselage = 3.79/2 = 1.895 m.
31.55
Hence, semispan of exposed wing = -1.895 = 13.89 m
2
Root chord of exposed wing = 5.47 – 0.264 x 1.895 = 4.97 m
Tip chord of exposed wing = 1.31 m
Consequently, Sexposed wing = 13.89 (4.97+1.31) = 87.23
(Swet)w = 2 x Sexposedwing x {1 + 0.2(t/c)}
= 2 x 87.23 {1 + 0.2 x (0.14)} = 203.76 m2
Kt = 1 + 0.31 + 0.21 = 1.52
(S wet )w 203.76
Hence, F1 = K t Cfe = 1.52×0.003055× = 0.00884
Sold 107.02
CDo -F1 0.0168 - 0.00864
F2 = = = 1.447×10-6 m2 /N
W/S old 5500
1 0.2078
= 1.0447 + = 0.0444
3.14×9.3 cos2 25o
0.0444
Hence, F3 =
q2
Finally, the drag polar in alternate form is:
0.0444 2
CD = 0.00884 + 1.447×10-6 p + p
q2
Answer: The drag polar in alternate form is:
0.0444 2
CD = 0.00884 + 1.447×10-6 p + p
q2
Having estimated the drag polar the wing loading can be optimised for the
prescribed flight speed Vp at altitude Hp.
Equation(4.13) can be rewritten as :
Tvp CDqp qp
tv = = = {F1 +F2 p+F3 p2 }
W p p
F
t v = qp 1 + F2 + F3 p (4.31)
p
The wing loading (pv) which will minimize tp , for chosen Vp , is obtained by
Or pv = F1 / F3 = qp F1 π Ae (4.32a)
.
A typical plot of tv vs p is shown in Fig.4.5. The value of optimum wing loading
from consideration of Vp is marked in the figure. The minimum thrust loading is
indicated by (tv)min. Allowing 5% increase (Δt) above this minimum, and solving,
Remark:
To select an engine, its rating needs to be prescribed. Further, the rating
of an engine is the sea level static take-off thrust produced by it. Hence, to get
the required rating of the engine, which can give tv at the chosen altitude, the
characteristics. If the ratio of the sea level static thrust (Ts.l)static to the thrust
available at Vp and Hp be Rtslbhvp , then the required thrust loading is given by:
Ts.l.
= R tslbhvp t v (4.33)
W
For example, let tv be 0.06 when Vp corresponds to a Mach number of 0.8 and
the flight altitude be 35000 ft (10670 m). Now, from Fig.4.22, which shows typical
characteristics of a turbofan engine, the sea level static thrust for cruise rating is
Example 4.4
Consider the airplane in example 4.3 and obtain (a) the optimum wing loading for
cruising at M = 0.8 at h = 11 km. (b) the range of wing loadings if an increase by
5 % is permitted above the minimum value of thrust required (c) the thrust
loading from cruising consideration, if the ratio of sea level static take-off thrust to
that at cruising condition, Rtslbhvp , be 5.
Solution :
The data are :
Mcr = 0.8, hcr = 11 Km
0.0444 2
Drag polar is CD = 0.00884 + 1.447 x 10-6 p + p
q2
At 11 km altitude, speed of sound (a) is 295.1 m/s
Vc = Vp = 0.8 x 295.1 = 236.10 m/s
ρ at 11 km is 0.364 kg/m3
1
qp = ×0.364×236.12 = 10145.3 N/m2
2
Hence, F3 = 0.0444/10145.32 = 4.314 x 10-10 m2/N2
Hence, optimum wing loading from consideration of Vp is
0.00884
pv = F1 /F3 = = 4526.9 N/m2
4.314×10-10
F
t
v min = qp 1 +F2 +F3p
p
F
= qp 2 1 +F2
p
0.00884
= 10145.3 2× +1.447×10-6 = 0.0543
4526.9
(b) If a 5% increase is permitted over ( tv )min then tv can be :
0.0543 x 1.05 = 0.05701.
From Eq.(4.31):
F
0.05701 = qp 1 +F2 +F3p
p
0.00884
Or 0.05701 = 10145.3 +1.447×10 -6 + 4.314×10 -10 p
p
Simplifying,
p2 – 9671.2 p + 20.49 x 106 = 0
Or p = 3135 or 6536.4 N/m2
Thus, if the wing loading is between 3135 to 6536 N/m2 then tvp will be with 5%
of tv .
min
(c) The ratio Rtslphp is given as 5.0. Hence, thrust loading for this airplane would
T
be = 5.0×0.0543 = 0.2715 .
W
Answers: (a) Optimum wing loading from consideration of Vp is 4527 N/m2.
(b) Range of wing loading when an increase by 5 % is allowed
T
over : 3135 to 6536 N/m2. (c) = 0.2715 .
W
Remarks:
(i) During the discussion in section 4.4.2 and in example 4.4, the prescribed flight
speed (Vp) is considered as Vcr. However, if the prescribed speed is above Vcr
then, in the case of high subsonic jet airplanes, the value of CDO may be higher
due to compressibility effects. This increase in CDO should be taken into account
and the value of F1 increased suitably. Section 3.3.5 of Ref.3.3 gives guidelines
for estimating changes in CDO and K in the transonic Mach number range.
airplane must have the ability to climb to service ceiling. The rate of climb at hcr is
generally 500 ft/min (or 152.4 m/min or 2.54 m/s) for jet airplanes. Hence, the
thrust to weight ratio required at hcr is
Vc
t = + tlevel
Vcr
In example 4.4, this would give :
2.54
t = +0.0543 = 0.0651
236.1
Taking Rtslbvp = 5.0
T
= 5×0.0651 = 0.326
W
This value of T/W is typical of airplanes with engines with bypass ratio of 6.
Dmin 1
tHmax = =
W (L / D)max
For a parabolic polar i.e. CD = CD0 + K CL2 , the following relations are already
known.
(CD) (L/D)max = 2C D0
Hence,
1 2CD0
tHmax = = = 4 CD0 K (4.34)
(L/D)max CD0 / K
on alternate break up of drag polar in Section 4.4.1). The variation of t Hmax with p
F1
t H = 2 qHmax ( +F2 ) (4.36)
p
and B in Fig.4.6. Allowing tHmax as deviation from the optimum value, the
permissible upper and lower limits on p are obtained as shown in Fig.4.6.
Example 4.5
For the airplane in example 4.4 obtain the optimum wing loading from the
consideration of absolute ceiling. It is prescribed that the absolute ceiling be 12
km and velocity be corresponding to (L/D)max. Further, obtain the range of wing
loadings if it is permitted that t Hmax can vary between 5% from that
corresponding to the aforesaid optimum.
Solution:
From Eqs.(4.35) and (4.36) :
tHmax = 4K(F1 +F2p) (E4.5.1)
F
tH = 2 qHmax 1 +F2 (E4.5.2)
p
From example 4.4, F1= 0.00884 and F2 = 1.447 x 10-6. CDo = 0.0168, K = 0.0444 .
From example 4.2, (W/S)old = 5500 N/m2. It is prescribed that qHmax corresponds
to V = V(L/D)max
(CL )(L/D)max = CDO / K = 0.0168 / 0.0444 = 0.615
W/S 5500
qHmax = = = 8943.1 N/m2
(CL )(L/D)max 0.615
Solving Eqs.(E 4.5.1) and (E 4.5.2) gives pHmax = 5500 N/m2 as it should be.
t Hmax corresponding to pHmax is :
tHmax1 = 0.05187
tHmax2 = 0.05733
Solving Eq.(E.4.5.1) for the above values of tHmax1 and tHmax2 gives:
Solving Eq.(E.4.5.2) for the above values of tHmax1 and tHmax2 gives:
Hence, to satisfy both the above criteria the range of permissible wing loading is
given by :
(a)higher of the lower limits from the two considerations and (b) lower of the
upper limits from the two considerations.
From Eqs.(E4.5.3) and (E4.5.4) the range of values is which satisfies these
limitations is :
5028 p 6084 N/m2
Answer: The range of wing loading to satisfy Hmax requirement is :
5028 to 6084 N/m2.
Remark:
This example again bring out the fact that the alternate representation of drag
polar permits examination of the effect of wing loading on Hmax.