The methods of minimization or maximization of lift, drag and pitching moment are extended to multiple coupled control surfaces with polynomial dependence of any degree on the control surface deflections. This serves as an introduction to the problem of minimization of cruise drag with pitch trim and unchanged lift, i.e. same airspeed and altitude, that is considered for a flying-wing configuration, using ten distinct strategies. Among the first four non-optimal strategies, only strategy II of using the centerbody elevator alone leads to drag reduction, albeit with a large deflection and increased angle-of-attack. There is drag increase for the strategies: (I) of equal deflection of all control surfaces; (III) preferential deflection of inner control surfaces with an aeroelastic limit of 7.5º; (VI) equal contribution to lift for all control surfaces. The two optimal strategies using Lagrange multipliers after a few interactions give: (strategy IV) a drag reduction for a good initial condition like strategy II; (strategy V) a drag increase for a poor initial condition like strategy I. A sub-optimal strategy VII of using multiples of optimal deflections, with a multiplication factor determined by lift equilibrium, does not change the angle-of-attack, and thereby increases drag; a drag decrease can be obtained (strategy VIII) freeing the constraint on angle-of-attack. The strategy IX of deflecting all control surfaces in two groups to minimize trim drag, leads to a larger angle-of-attack, and thereby increases drag. The strategy X of using the two most effective control surfaces with opposite drag slopes gives the best compromise: (i) Lift balance with unchanged angle-of-attack; (ii) Second best drag reduction with small control surfaces deflections. The small control surfaces deflections justify the use of linear aerodynamics and reduce the risk of adverse aeroelastic effects.
Multiple control surfaces, Flying-wing configuration, Optimization, Cruise drag, Pitch trim
c¯ : Mean Aerodynamic Chord
CDi : Dimensionless drag coefficient of the control surface i
C′Di : Derivative of the drag coefficient of the control surface i
C′Di : Derivative of the drag coefficient with regard to the deflection of the control surface i
C′Mi : Derivative of the drag coefficient
CMi : Pitch control coefficient of the control surface i
C′Mi : Derivative pitch control coefficient with regard to the deflection of the control surface
D : Trim drag
Li : Moment arm of the control surface i
li : Given c.g. position of the control surface i
M : Pitching moment
q : Dynamic pressure
V : Airspeed
Si : Area of the control surface i
S0 : Total control surface area
S: Wing area
α: Angle of attack
α0: Cruise angle of attack
δi : Deflection of the control surface i
δ¯i : Optimal deflection of the control surface i
λ : Lagrange multiplier
ρ : Mass density of atmospheric air
AOA: Angle of Attack
BWB: Blended Wing Body
c.g.: Center of Gravity
FW: Flying-Wing
Ma: Mach Number
The minimization of cruise drag is an instance of the general problem of maximization or minimization control forces and moments (6), that is simpler assuming: (i) Decoupled control surfaces such that the deflection of one does not affect the others; (ii) Forces and moments specified by quadratic functions of the deflections of control surfaces that are sufficient to specify base values for equilibrium, first-order derivatives for stability and second-order derivatives for optimization. The maximization or minimization of control forces and moments can be extended to (i) Interacting control surfaces, such that the deflection of one affects the others (6.1); (ii) Forces and moments specified by polynomials of any degree of the deflection angles of control surfaces, for higher accuracy. The (i) Coupled and (ii) High-order controls (6.3) Increase the data needed to fully address (6.4) The performance, stability and optimization problems.
The main aim of the paper in (7-10) to compare ten pitch trim strategies for cruise drag minimization that is most relevant to reduce fuel consumption, emissions, and operational costs [1-13]. The case of a blended-wing body (BWB) configuration [14-26] is used to compare eight pitch trim strategies for cruise drag minimization: (i) The four simplest non-optimal pitch trim strategies, which give a first feeling for the flight physics of the problem (7); (ii) The preceding can serve as initial condition for the two optimal pitch trim strategies, which are recursive and inevitably more complex (8); (iii) The preceding suggest four further pitch trim strategies which, although sub-optimal, may prove quite effective and relatively simple (9-10).
The comparison of the ten pitch trim strategies (11) was initially motivated by (i) Reduction of cruise drag, but raises other relevant and related issues, e.g.: (ii) If unchanged lift is required, to keep the same airspeed and altitude as for the untrimmed flight condition, can this be achieved with the same untrimmed angle-of-attack (AOA) or is a higher value needed?; (iii) The largest control surface deflections required are sufficiently small to avoid concerns about adverse aeroelastic [27-31] or aerodynamic [32-57] effects, e.g. loss of control effectiveness [58-65] or control inversion, and shock wave formation and boundary layer separation?; (iv) The AOA of the aircraft and deflections of all control surfaces are within the range where linear or weakly non-linear methods can be applied with confidence? It should be borne in mind that there are many approaches to drag reduction [66,67], and the optimization of control surface deflections is just one of them. The comparison of the ten pitch trim strategies according to these criteria shows that the most complex or demanding is not necessary the best; it is possible to find a relatively simple strategy X, which achieves pitch trim drag reduction relative to the untrimmed flight condition, while also keeping the low untrimmed AOA and using moderate deflections of only the two most effective control surfaces (10).
Considering an aircraft with multiple control surfaces, the choice of deflections which maximize and minimize any force or moment is considered without making the assumptions that: (a) The controls are independent or decoupled; (b) The control forces and moments are specified by quadratic functions of the deflections. Using these assumptions (a,b) leads [68] to the simplest application of the method to find the extrema (maxima and minima) of the control forces and moments, whether they lie inside the interval of possible deflections or at the limit deflections. The methods presented next extend earlier results [68] and apply (6) to the case of coupled or interdependent controls (6.1) and to the case of control forces and moments specified by a polynomial or arbitrary degree of the deflections (6.2), thus removing one at a time of the two restrictions (a) and (b) above; both restrictions can be removed simultaneously in the general case of coupled controls specified by high-order polynomials (6.3), which requires more aerodynamic data (6.4). The forces and moments along the three axis are denoted (1b) by Cα with the index α running (1a) from one to six:
α=1,2,...,6: Cα={X,Y,Z,L,M,N} (1a,b)
the index i numbers the N control surfaces:
i=1,...,N=5: i={body flap,inner flap,middle flap,outer flap,rudder} (2a,b)
In the case of interdependent or coupled controls, the force and moment on one surface are affected at most by the deflections of all others:
Cαi=Cαi0+∑i=1Naαijδj+∑j,k=1Nbαijkδjδk (3)
Assuming still a polynomial of second degree: (i) The static values are unchanged; (ii) The slopes are specified by a matrix (4a) indicating the effect of the deflection of control surface δj on the forces and moments on control surface Cαi:
aαij=∂Cαi∂δj; 2bαijk=∂2Cαi∂δj∂δk (4a,b)
(iii) The curvatures (4b) involve a multiplicity with three indices.
The extremum cannot be found separately for each surface, but only for all together:
0=dCα=∑i=1NdCαi=∑i,j=1Ndδj(aαij+2 ∑k=1Nbαijk δ¯k) (5)
where the multiplicity of slopes (4b) is symmetric in the last two indices (6a)
bαijk=bαikj ⇔ ∂2Cαi∂δj∂δk=∂2Cαi∂δk∂δj (6a,b)
on account of (6b). The deflections for the extrema are given by (7)
δ¯i,α=−12∑j,k=1Naαij(bαijk)−1 (7)
where (bαijk)−1 is the inverse matrix of bαijk with regard to jk. The extremum will be a maximum/minimum if the second-order differential is positive/negative:
d2Cα=∑Ni=1d2Cαi=∑Ni,j,k=1(∂2Cαi∂δi∂δj) dδidδj (8a,b)=2∑Ni,j,k=1bαijkdδjdδk{ > 0 < 0minimum,maximum;
this is equivalent to the matrix bαijk being positive/negative definite in (j,k) for every i, i.e. having all eigenvalues real and positive (negative):
bαijk={is positive-definite implies δ¯i is local minimum, is negative-definite implies δ¯i is local maximum. (9a,b)
If these conditions [(9a) or (9b)] are not all met for some (i,α) or the values (7) lie outside the range of possible deflections, the maximum/minimum occur at the extreme deflections. The case of coupled controls (3) reduces to independent controls when for each α the matrix aαij and multiplicity bαijk are replaced by two vectors (aαi,bαi).
Next the assumption of decoupled controls is retained, but the control forces and moments are allowed to be polynomials of any degree M:
Cαi=∑m=1MA(m)αi(δi)m (10)
where: (i) The static values (11a), slopes (11b) and curvatures (11c) are retained:
A(0)αi≡Cαi(0)≡Cα0, A(1)αi=dCαidδi≡aαi, (11a-d) 2 A(2)αi=d2Cαidδi2≡2 bαi, m!A(m)αi≡dmCαidδim,
and derivatives of order m up to M are added. The extremum can be obtained independently for each decoupled surface. The total extremum is the sum of partial extrema. Thus only partial extrema need to be considered. The extremum corresponds to the vanishing of the derivative of a polynomial:
0=dCαidδi=∑m=1MA(m)αim(δ¯i)m−1 (12)
In the case of second-degree there would be only one root, whereas in the case of degree M there are M-1 roots. Of these roots, those with positive/negative second-order derivative can be local minima/maxima:
d2Cαidδ¯i2=∑m=2MA(m)αim(m−1)(δ¯i)m−2 (13a-b){> 0 < 0implies δ¯i is local minimum,implies δ¯i is local maximum;
If there are several local minima/maxima then the lowest/highest is chosen as the overall minimum/maximum.
It may happen that the second order derivative vanishes and some higher-order derivatives vanish also. Not all derivatives up to order M can vanish otherwise the polynomial would reduce to a constant. Suppose that the first derivative which does not vanish at δ¯i is of order s; if s is odd the δ¯i is an inflexion point, not a maximum or a minimum (14a):
Cαi′(δ¯i)=...=C(s−1)αi(δ¯i)=0≠C(s)αi(δ¯i) (14a,b){s is odds is evenimplies δ¯i is an inflexion-point, implies δ¯i is an extremum;
If s is even then δ¯i is an extremum (14b) and the sign decides if it is a maximum or minimum.
Cαi′(δ¯i)=...=C(2s−1)αi(δ¯i)=0≠C(2s)αi(δ¯i) (15a,b){>0 < 0implies δ¯i is local minimum,implies δ¯i is local maximum.
All local maxima/minima, arising at second-order derivative (13a,b) or at any derivative s of higher even order (15a,b) are considered for the overall maximum and minimum:
Cαimax=sup{Cαi(δ¯i)}≥inf{Cαi(δ¯i)}=Cαimin (16)
If the deflection δ¯i lies outside the range of possible deflections then the corresponding extremum is at one of the deflection limits. The extrema for each decoupled control specify by addition the overall extrema.
The optimization method applies most simply to decoupled controls specified by second-order polynomials. These two restrictions were relaxed to allow either coupled controls (6.1) or high-order polynomials (6.2). Next both constraints are removed, so that the controls can be coupled and specified by polynomials of any degree. If the highest degree is M, all polynomials can be written, possibly with some zero coefficients, in the form:
Cαi=Aαi+∑j=1NAαijδj+∑j,k=1NAαijkδjδk+...=∑m=0M∑j1,...,jm=1NAαij1...jmδj1...δjm (17)
where: (i) The static values are unchanged; (ii) The slopes (18a) and curvatures (18b) are as before in the coupled case (14a,b).
aαij=∂Cαi∂δj=Aαij, 2bαijk=∂2Cαi∂δj∂δk=2Aαijk, ∂mCαi∂δj1...∂δjm=m!Aαij1...jm (18a-c)
(iii) The higher-order derivatives (18c) may all be designated by the same letter A and distinguished by the number of indices (j1,...,jm) for 0≤m≤M.
Due to the coupling the extremum must be sought jointly for all controls together:
0=dCα=∑i=1NdCαi=∑m=1Mm∑i=1N∑j1,...,jm=1NAαij1...jmδ¯j1...δ¯jmdδj1 (19)
the roots of (19) specify the extrema. If the second variation is positive/negative the extremum is a minimum/maximum
d2Cα=∑i=1Nd2Cαi=∑m=1Mm(m−1)∑i=1N∑j1,...,jm=1NAαij1...jmδ¯j1...δ¯jmdδj1dδj2
{ > 0 < 0implies δ¯i is local minimum,implies δ¯i is local maximum; (20a-c)
this is equivalent to stating that the matrix (21a) is positive/negative definite, i.e. as all eigenvalues real and positive (negative)
Xj1j2=∑m=2Mm(m−1)∑i=1N∑j1,...,jm=1NAαij1...jmδ¯j3...δ¯jm
{positive-definite implies local minimum,negative-definite implies local maximum. (21a-c)
The matrix (21a) may not be definite, e.g. it may be indefinite, positive or negative semi-definite, i.e. have eigenvalues with opposite signs, or some zero eigenvalues; in this case δ¯i is not an extremum. If all eigenvalues of (21a) are zero, then it is necessary to proceed to higher order to know if δ¯i is an extremum.
If the differential (22a) of lowest order which does not vanish is of odd order (22b) then δ¯i is an inflexion, not an extremum:
dCαi(δ¯î)=...=d(s−1)Cαi(δ¯î)=0≠d(s)Cαi(δ¯î)=X: (22a-d)⎧⎩⎨⎪⎪s is odd implies δ¯i is inflexion point,s is even { > 0 < 0implies δ¯i is local minimum,implies δ¯i is local maximum;⎫⎭⎬⎪⎪
If the first non-vanishing differential is of even order then δ¯i is an extremum (22c,d) viz. a local minimum/maximum if it is always positive/negative. If X does not have a fixed sign then δ¯i is not an extremum. If the value δ¯i lies outside the range of possible deflections then the extremum is at one of the extreme deflections. The set of values δ¯i is then used to calculate the overall maximum and minimum in (17).
The number of data points needed Q0 is lowest (6) for the three forces and moments (total 6), for N decoupled control surfaces (6×N), for second-degree polynomials with 3 coefficients each:
Q0(N)=6×N×3=18N (23)
If the controls remain decoupled but the polynomials are (10) of degree M, the 3 coefficients are replaced by M + 1, and the number of data points need is:
Q1(N,M)=6×N×(M+1)=6N(M+1) (24)
this reduces to (23) for M = 2. In the case (3) of a polynomial of degree M = 2 and N coupled controls, there are 6N static values, 6N2 first-order derivatives and 6N3 second-order derivatives; the latter are symmetric (6a), reducing their number to 6N2(N+1)/2 in:
M=2: Q2(N)=6[N+N2+N2(N+1)/2]=3(2N+3N2+N3) (25)
In the general case (17) of polynomials of degree M and N coupled controls for each force and moment there is one static value, N2 slopes, N2(N+1)/2 curvatures and
m=1,2,...: N(N+m−1m)=(N+m−1)!Nm!(N−1)! (26a)=N2(N+1)(N+2)...(N+m−1)m!
coefficients of order m that are symmetric in all indices; the total number of data points is:
Q3(N,M)=6N{1+N+N(N+1)/2+...+N...(N+M−1)/M!} (26b)=6N{1+∑m=1M(N+m−1m)}
For polynomials of second degree M=2 the number (27a) of data points (26b) agrees with (25), and it increases for polynomials of third degree to (27b):
Q3(N,2)=Q2(N); Q3(N,3)=6N[1+N+N(N+1)/2+N(N+1)(N+2)/6] (27a-b) =6N+11N2+6N3+N4
As can be seen in the Table 1 the number of data points needed increases rapidly with the degree of the polynomial, more so in the coupled case.
The interaction between control surfaces farther apart would be weaker, and could be neglected to reduce the number of data points needed. Suppose that there are N control surfaces each interacting only with its neighbors: (i) The first and last surfaces have each one neighbor, and thus one coefficient of degree zero, two of degree one (δ1,δ2), three of degree two [δ12,δ22,δ1δ2], four of degree three [(δ1)3,(δ1)3δ2,δ1δ22,(δ2)3], and (m+1) of degree m; (ii) The (N−2) control surfaces not at either end have two neighbors, and hence one coefficient of zero degree, three of first degree (δ1,δ2,δ3), six of second degree [δ12,δ22,δ32,δ1δ2,δ1δ3,δ2δ3] and (m+1)(m+2)/2 coefficients of degree m; (iii) These coefficients add for all degrees from 0 to M, for each of three components of forces and moments:
N≥2: Q4(N,M)=6{2∑m=0M(m+1)+(N−2)∑m=0M(m+1)(m+2)/2} (28)=6(M+1)(M+2)+3(N−2)∑m=0M(m+1)(m+2)
This is the number of data points in the case of N control surfaces coupled only to the nearest neighbor, and polynomial of degree M. For two control surfaces the number of parameters is: (i) The same for neighboring interaction (28) or all interactions (26b); (ii) Larger than for decoupled controls (24):
Q4(2,M)=6(M+1)(M+2)=12∑m=0M(m+1)=12∑m=0M(m+1m) (29) =Q3(2,M)>12(M+1)=Q1(2,M)
For more than two control surfaces (30a), the number of parameters for neighboring interaction (28) is: (i) Larger than (24) for decoupled controls (30b):
N≥3: Q4(N,M)>(M+1)[12+6(N−2)]=6N(M+1)=Q1(N,M) (30a,b)
(ii) Smaller than (26b) for all controls coupled (30c,d):
N≥3: Q4(N,M) < 6(N−2)∑m=0M(m+1)(m+2)=6(N−2)∑m=0M(m+22) < 6N∑m=0M(m+2m)≤6N∑m=0M(N+m−1m) < Q3(N,M). (30c-d)
In the case of polynomials of degree two (31a) and three (31b) the number of parameters for neighboring interaction:
M=2; N≥2: Q4(N,2)=72+60(N−2)=60N−48 (31a)
M=3; N≥2: Q4(N,3)=120+120(N−2)=120(N−1) (31b)
is indicated in Table 1 in comparison with decoupled (24) and fully coupled (25,27) controls. The case of each control surface interacting only with its neighbor(s) reduces significantly the number of data points needed, and should account for control coupling with good accuracy. Since the method applies both to coupled and uncoupled controls, and polynomials of any degree, it can be illustrated for the simplest case, for which data is more readily available.
Before proceeding to implement optimal pitch trim strategies, it is useful to start with some simple initial pitch trim conditions, for three reasons: (i) They provide a more direct understanding of cruise equilibrium conditions; (ii) They serve as a starting point for the iterative optimisation process; (iii) They may suggest sub-optimal and non-optimal strategies which are may be simpler but still quite effective. For a given cruise condition, the effect on drag due to pitch trim is compared for three simple cases, viz.: (i) Equal deflection of all control surfaces (strategy I in 7.1); (ii) Deflection of elevator alone (strategy II corresponding to strategy A in [68]); (iii) Preferential deflection of inner control surfaces with moderate upper limit (strategy III in 7.2). The optimal deflection using Lagrange multipliers use two distinct initial conditions: (a) Deflection of elevator alone (strategy IV corresponding to strategy B in [68]); (b) Equal deflection of all control surfaces (strategy V in 7.3).
Consider a BWB aircraft configuration at a cruise weight W = 715 t, intermediate between initial and final cruise weights. At cruise altitude z = 35 kft the sound speed is c=295 m.s−1, and the cruise Mach number Ma=0.85 corresponds to an indicated airspeed U=Ma×c=251 m.s−1. At the tropopause the mass density is ρ = 0.223 kg m-3, and thus the dynamic pressure is q=12 ρ U2=7.025×103 kg m−1s−2. For a reference wing area S = 961.35 m2, this corresponds to a cruise lift coefficient:
CL=2W/ρ S U2=W/q S=0.10588. (32)
The lift coefficient CL, and also the drag CD and pitching moment CM coefficients, are given in Table 2 for the chosen BWB aircraft configuration in a typical cruise condition. Linear interpolation of the values of the lift coefficient between zero α = 0° and three degree α = 3° AOA, shows from Table 2 that (25) corresponds to a cruise AOA
α0 = 0.71789°. (33)
Further linear interpolation in the Table 2 shows that this AOA corresponds to a cruise drag coefficient (34a):
CD = 0.00569, CM = - 0.02441, (34a,b)
and pitching moment coefficient (34b). The total dynamic force q¯=q S=6.7535×106N, then leads to a cruise drag D=CD q¯=3.8427×104N and pitching moment M=CM q¯ c¯=−5.2699×106 N m. The cruise lift-to-drag ratio is L/D = CL/CD = 18.6.
There are five sets i = 1,...,5 control surfaces available for pitch control whose areas are indicated in Table 3 together with the fraction of wing area Si/S, and fraction of total control surface are Si/S0, where:
S0≡∑i=15 Si. (35)
The lift CLi, drag CDi and pitching moment CMi coefficients for each surface i = 1,...,5, and those for all surfaces together i = 0, are given in Table 4 for upward or downward deflections of three degrees δi = ± 3°, at angles of attack α = 0° and α = 3°. From these are interpolated the values at cruise AOA α0 = 0.71789° indicated in the same Table 4. The first and second-order derivatives of the aerodynamic coefficients are calculated applying finite differences to the values in Table 4 leading to the results in Table 5.
The preceding tables contain all the data needed to implement various trimming strategies. The pitching moment coefficient at cruise (34b), must be balanced by the opposite trimming pitching moment, e.g. adopting the strategy I of equal deflection δI of all trailing-edge control surfaces (36a):
CMR≡0.02441=CM0(δI)⇒δI=4.4680º≡δIi, (36a,b)
then linear interpolation in Table 5 leads to the deflection (36b) indicated in the first column of Table 6. Since the latter deflection is small, the use of linear aerodynamics and linear interpolation is justified. The deflection (36b) leads from Table 4 to a drag coefficient (37a):
CDI≡CDO(δI)=0.00659, Δ CDI≡CDI−CD=0.00089, (37a,b)
which exceeds the cruise drag coefficient in Table 1 by 89 drag counts. Thus there is a 15.6% drag increase:
Δ CDI/CD=0.156, (38)
by adopting the strategy I of pitch trim by equal deflection of all control surfaces. This is a relatively large cruise drag penalty to achieve pitch trim, so another alternative is sought.
The strategy I of equal deflections δI of all control surfaces is unlikely to be good, because it does not exploit the differences between the control surfaces, in their ability to produce pitching moment and reduce drag. It is seen from Table 5, that the centerbody elevator has: (i) One of the largest positive pitching moment slopes, hence can trim the aircraft with smaller deflections; (ii) Has the largest (in modulus) negative drag slope, so that drag is reduced from small deflections. Thus suggests that the preceding strategy I be compared next with the strategies of maximal deflection of the inboard control surfaces. The strategy II of deflecting only the innermost trailing-edge control surface, namely the central body elevator was considered as strategy A before [68], leading to a large elevator deflection shown in the second column of Table 6. A large elevator deflection could have undesirable consequences, e.g. shock wave formation and boundary layer separation, or exciting aeroelastic modes and causing earlier onset of buffet. Thus it may be necessary to limit the maximum deflection of all control surfaces to a given value, say 7.5º. In this case the elevator deflection is limited, and control surfaces further outboard must be used, in strategy III until the required pitching moment is obtained.
In the strategy III, the body elevator is limited to say, a seven-and-a-half degree deflection (39a) to avoid aerodynamic, aeroelastic or other problems:
δIII1 = 7.5º, CMIII1 = CM1 (δIII1) = 0.01828, (39a,b)
the resulting pitching moment coefficient (39b) falls short of the required value (36a) by (40):
Δ CMIII1≡CMR−CMIII1=0.00560. (40)
Thus the strategy III of limiting the deflections to 7.5°, requires more than the elevator for pitch trim, as could be expected from the second column of Table 6 showing that the elevator alone would require a 9.0° deflection for pitch trim. The short-fall (40) ≡ (41a) in pitch trim would be covered by the outer elevator:
0.00560=Δ CMIII1=CM2(δIII2)=δ=III212.055º, (41a,b)
deflecting by (41b), which exceeds the 7.5° limit. Enforcing the latter limit for the outer elevator (43a):
δIII 2 = 7.5°, CMIII2 = CM2 (δIII2) = -0.00400, (42a,b)
leads to a negative contribution to the pitching moment (42b), so that the shortfall is now higher:
Δ CM III 2=ΔCM III 1−CM III 2=0.00960. (43)
This shortfall is not compensated by the inner flap at the maximum allowed deflection:
δIII 3 = 7.5°, CM III 3 = CM 3 (δIII 3) = 0.00536, (44a,b)
and the remaining shortfall (45a):
0.00424 = Δ M III 2 – CM III 3 = Δ CM III 3 = CM III 4 = C4 (δIII 4), (45a)
is covered by deflecting the outer flap by (45b):
δIII 4 = 3.5102°; δIII 5 = 0°, (45b,c)
thus the ailerons need not be deflected. Thus the strategy III of preferential deflection of inboard control surfaces, not exceeding the 7.5° limit, leads in the third column of Table 6 to the maximum deflections of the elevators (39a,42a) and inner flap (44a), smaller deflection of the middle flap (45b) leaving only the ailerons (45c) unused.
The Table 5 specifies drag coefficients corresponding to these deflections (46a-e):
CD III 1=CD 1(7.5º)=0.00404+0.000016×7.52=0.00494, (46a)
CD III 2=CD 2(7.5º)=0.00531+0.0000022×7.52=0.00543, (46b)
CD III 3=CD 3(7.5º)=0.00531+0.000013×7.52=0.00604, (46c)
CD III 4=CD 4(3.5º)=0.004615, (46d)
CD III 5=CD 5(0º)=0.00569, (46e)
where curvature correction was applied, for the first three deflections, which are larger. They are multiplied by the fraction of the total control surface area in Table 3, to give the trimmed drag:
CD III=∑i=15CD III i Si/S0=0.00570, (47)
which leads to a trimmed drag a little higher than the untrimmed drag (48a):
Δ CD III ≡ CD III – CD = 0.00001, Δ CD III / CD = 0.0018, (48a,b)
corresponding to 1 drag counts or a 0.2% penalty. The three non-optimal cruise pitch trim strategies are compared in Table 6: (i) The strategy I of equal collective deflection of all trailing-edge control surfaces leads to a 89 counts of drag as penalty; (ii) The drag count penalty is reduced to 1 by the strategy III of preferential deflection of inboard control surfaces with a 7.5° degree limit; (iii) The most elementary strategy II of elevator deflection alone gives a cruise drag reduction of 16 counts due to pitch trim. The result that the classical strategy II is the best so far from the point of view of drag, need not to be general, and could be a consequence of the particular aerodynamic data used in the present calculation.
The optimal strategy, after convergence, should lead to the same optimal deflections, for any initial condition. If, however, as is always the case in practice, a finite number of iterations is used, the iterated 'optimal' deflections will depend on the initial condition. The optimal strategy IV with four iterations starting with the non-optimal strategy II of elevator deflection alone corresponding to strategy B in [68], is compared next with the optimal strategy V with six iterations starting with the non-optimal strategy I of equal control surface deflections (7.1). The former strategy IV was shown in Table 6 in [68], and the latter strategy is shown in Table 7, where the calculations were performed with more digits in the first than in the second case. The optimal deflections for the latter strategy V:
δVi≡[δ(6)i+δ(5)i]/2={4.125º,0.02º,1.88º,0.01º,−3.685º}, (49)
are smaller than for strategy IV.
The distinct deflection for strategies IV and V, lead to different drag coefficients for the latter, from Table 2, viz.:
CDV1≡CD1(4.125º)=0.00618, (50a)
CDV2≡CD2(0.002º)=0.00776, (50b)
CDV3≡CD3(1.68º)=0.00768, (50c)
CDV4≡CD4(0.01º)=0.00846, (50d)
CDV5≡CD5(−3.685º)=0.00803. (50e)
The corresponding total drag (79a).
CDV=0.00770, ACDV≡CDV−CD=0.00201, (51a,b)
is 201 drag counts above (51b) the untrimmed value, for a 35.3% cruise drag penalty. From the latter it is clear that the 'optimal' deflections after 4 accurate iterations with a good initial value, give better results than six less accurate iterations starting with a worse value. The comparison of drags obtained by the methods I to V could be misleading because they correspond to different lifts, and hence to distinct flight conditions, such as different speed, altitude or both. A valid comparison of all strategies requires that pitch trim be achieved with the same lift (8).
The three pitch trim strategies I, II and VI are considered with constant lift, thus leading to modification of the strategies: (I) Of equal deflection of all trailing edge surfaces (8.1); (II) Of deflection of the centerbody elevator only (8.2); (VI) equipartition of the cruise lift contribution from each control surface (8.3).
The strategy I of equal deflection of all control surfaces is applied using the following procedure, which finds directly the equilibrium of lift, drag and pitching moment in cruise. The lift coefficient depends on AOA α, and for pitch trim strategy I on the collective angle of deflection δI for all control surfaces. For a given weight, altitude and airspeed, the lift coefficient cannot change, i.e. it takes the value (32); then must hold for small angles:
0.10588 = CL = 0.04411 + α CL α + δI CLδ0, (52)
where the derivatives are calculated at zero angles α = 0° = δI, viz. (47a) from Table 2:
CL α=(∂ CL/∂ α)δ0=0=[CL(α=3º)−CL(α=0º)]/3º (53a)=(0.30224−0.04411)/3º=0.086043,
CL δ0=(∂ CL0/∂ δ0)α0=0=[CL0(δ0=+3º)−CL0(δ0=−3º)]/6º (53b)=(0.15977−0.49675)/6º=−0.056163,
and (47b) from Table 4. It follows that for constant lift, the AOA α and collective control surface deflection must be related by:
δI = 1.5320 αI - 1.0998, αI = 0.65273 δI + 0.71787, (54a,b)
allowing an infinite choice of one parameter.
In the same small angle approximation, the pitching moment coefficient is given by:
CM(α,δ)=CM(0,0)+α CM α+δI CM δI. (55)
At constant lift (54a,b) it depends only on one parameter, e.g. AOA:
CM(α)=−0.00819−0.02259 αI+0.028935 (1.5320 αI−1.0998), (56)
where CMα was calculated from Table 2 and CMδ from Table 4. For pitch trim the pitching moment coefficient must be zero, leading to the angle of attack (57a):
CM(αI)=0: αI=1.8406º, δI=1.7200º, (57a,b)
and collective control surface deflection (57b). The drag coefficient cannot be accurately predicted by a small angle approximation such as:
CD(α,δ)=CD(0,0)+α CDα+δ CDδ, (58a)
and thus a different method is used as shown in Table 8: (i) The values of the drag coefficient in Table 4 are interpolated for the AOA (57a), for fixed deflections δ=±3º; (ii) Then the values in Table 8 are extrapolated for the deflection (57b):
C∗DI=CD0(δI)=CD0(1.7200º)=0.00766. (58b)
Thus the strategy I with cruise lift equilibration leads to a drag penalty (59a) or 197 counts, or 34.5% in (59b):
ΔCDI=C∗DI−CD=0.00197, ΔCDI/CD=0.345. (59a,b)
Since the deflection of control surfaces upwards (51b) decreases lift, the AOA is increased (57a), and hence also the drag, as shown on Table 9. In the case of the strategy I of equal deflection of all control surfaces, the condition of cruise lift equilibration causes trim drag penalty of 197 drag counts. The opposite will be shown next for the strategy II of elevator deflection only, when the cruise lift equilibration condition decreases the pitch trim drag by 18 counts, as shown next.
The direct, non-iterative approach to simultaneous determination of AOA and control surface deflection, used for strategy I in 7.1, applies also to other strategies, as shown next. For comparison, consider the decrease in cruise drag due to pitch trim is calculated next for strategy II of deflection of the centerbody elevator alone. Using (53a) and:
CLδ1≡(∂CL1/∂δ1)α=0=−0.011245, (60)
from Table 4, the condition (52) of no lift change is now replaced by:
δII=7.6517 αII−5.4931. (61)
The larger values in (61) compared to (54a), are due to deflecting only one surface instead of all by the same amount. Using the analogue of (55) for pitch trim:
0=CM(αII)=−0.00819−0.02259 αII+0.0053067 (62)×(7.6517αII−5.4931),
leads to the AOA (63a):
αII=2.0727°, δII=10.367°, (63a,b)
and elevator deflection (63b), from (61).
The drag is again interpolated from Table 2 to the AOA (63a), for fixed deflection δII=±3° as shown in Table 6, leading to:
C∗DII2=CD1(δII)=0,00319, (64a)
which is very low because ∂CDII/∂δII < 0, leading to a large reduction for large δII. It is not reasonable to assume constant slope, i.e. neglect curvature effects, for large δII. Thus a curvature correction made to (64a), viz.:
C∗DII2=0.00319+12C′′D2 (δII)2=0.00491, (64b)
where C′′D2 = 0.000032 from Table 5. The non-deflected surfaces i = 2, 3, 4, 5 influence drag and thus the total drag is now:
C∗DII≡C∗DII1 S1/S0+C∗DIIi(1−S1/S0)=0.00551. (65a)
The strategy II of elevator deflection only, with cruise lift equilibrium, leads to a drag reduction due to pitch trim (65b) of 18 counts or -3.2% in (65c):
ΔC∗DII≡C∗DII−CD=−0.00018, ΔC∗DII/CD=−0.032. (65b,c)
After showing the effect of the cruise lift equilibrium condition on the worst (I) and best (II) non-optimal pitch trim strategies, it is worthwhile to consider a non-optimal pitch trim strategy VI, for which cruise lift equilibration is taken as the starting point.
A more direct strategy VI, is to conserve the total lift coefficient:
0.10588=CL=0.04411+α CLα+∑i=15δiC''Li, (66)
taking equal contributions from each control surface:
0.012338=0.2αCLα+δViC''Li. (67)
It follows that the deflections are related to the AOA by:
δVI1=−1.0724+1.49575α, (68a)
δVI2=−3.5001+4.88189α, (68b)
δVI3=−1.2609+1.7581α, (68c)
δVI4=−6.86971+9.5817α, (68d)
δVI5=−2.473+3.45105α, (68e)
using (59a) and taking C''Li at δi=0 from Table 5; the pitching moment:
CM=−0.00819−0.02259αIV+∑i=15δIViC′Mi, (69)
vanishes at the trim conditions:
0=CM=−0.00819−0.02259 αIV+ 0.00546 (1.49575αIV -1.072)+
+0.0015185(4.8819αIV−3.5001)+0.004075(1.7587αIV−1.2609)+ (70)
+0.0087615(9.5817αIV−6.8697)+0.003297(3.45105αIV−2.4743),
where C′Mi were taken at δi=0 from Table 5. The trim condition corresponds to (71a) the AOA:
αIV=2.2794º, δIV i={6.6990º,7.62770º,2.7478º,14.971º,5.3920º}, (71a,b)
and the corresponding (68a-e) deflections are (71b).
The drag coefficients are given in Table 8 for the AOA (71a) and two deflections δi=±3º, based on Table 4; extrapolations for the deflections (71b), leads to the drag coefficients:
CD VI 1≡CD1(δIV 1)=0.00529+12 C′′D1(δIV1)2=0.00601, (72a)
CDVI2≡CD2(δIV2)=0.00727+C′′D2(δIV2)2=0.00730, (72b)
CDVI3≡CD3(δIV3)=0.00809+C′′D3(δIV3)2=0.00818, (72c)
CDVI4≡CD3(δIV4)=0.00640+C′′D4(δIV4)2=0.01660, (72d)
CDVI5≡CD5(δIV5)=0.00913+13C′′D5(δIV5)2=0.00936, (72e)
where the curvature corrections were applied using values from Table 5:
C′′Di={0.0000323, 0.00000133, 0.0000257, 0.0000910, 0.0000156}. (73)
The large deflections (71b), lead to relatively large drag coefficients (72a-e) which are weighted by the relative control surface areas in Table 3, to yield the total drag coefficient:
CDVI=∑i=15CDVI(Si/S0)=0.01089; (74)
this leads to a trim drag penalty of 520 counts (75a) or 91% in (75b):
ΔCDVI≡CDVI−CD=520, ΔCDVI/CD=0.914. (75a,b)
As shown in Table 9 there (75a,b) can be a large difference in trim drag according to pitch trim strategy, suggesting that optimal deflections to minimize cruise drag should be considered.
The double iteration loops needed to make optimal deflections (8.1-8.2) compatible with cruise lift equilibrium (8.3), can be avoided by sub-optimal strategies. A sub-optimal strategy uses optimal values in a non-optimal way, e.g. in the sub-optimal strategy VII, the optimal deflections (8.1-8.2) are multiplied by a constant, which is specified by the condition of cruise lift equilibrium (9.1). The ideas arising from optimal strategies may be used in simplified strategies, which are much more effective with less labour, e.g. the strategy of deflecting differently two distinct group of control surfaces, in two ways (strategies VIII and IX) as to reduce drag while producing the required pitching moment for each of them (9.2). This idea is further refined an strategy X, which deflects not all control surfaces, but only the two most effective in producing pitching moment with drag reduction (9.3).
The optimal deflections can be the used in a non-optimal may to specify a simpler sub-optimal strategy VII, which multiplies the deflections by a constant k, leading to the condition of constant lift (46) in the form:
0.10588=CL=0.04411+αCLα+k∑i=15δiCLi; (76a)
using (69a) and the second set of optimal deflections (63) leads to:
∑i=15δVIIiC'Li=−0.047728; (76b)
(72a,b) together with (33a) specify the ratio of deflections k in terms of the AOA:
k=1.8027αVII−1.2942. (76c)
The corresponding pitching moment coefficient (49) is given by:
CM=−0.00819−0.02259αVII+k∑i=15δVIIiC′Mi, (77a)
where from (83a) and (77):
∑i=15δVIIiC′Mi=0.018152; (77b)
using (83b) the condition of pitch trim is:
CM=0=−0.00819−0.02259αVII+0.018152k. (77c)
The two conditions of (i) constant lift and (ii) pitch trim are satisfied substituting (76c) in (77c):
0.00819+0.02259αVII=0.018152(1.8027αVII−1.294), (78)
which specifies the AOA (89a):
αVII=3.1267º, k=4.3424, (79a,b)
and coefficient (89b) for the final deflection angles:
δVIIi=kδVi={17.912º, 0.0868º, 8.1637º, 0.0434º, −16.002º}. (80)
The drag coefficients are calculated from Table 4 for the AOA (89a) and listed in Table 8 for the deflections δ = ±3°. From this Table 8 the drag coefficients can be interpolated for the deflections (90), viz.:
CDVII1=CD1(17.912º)=−0.00137+0.000016×17.9122=0.00376, (81a)
CDVII2=CD2(0.0868º)=0.01055, (81b)
CDVII3=CD3(8.1637º)=0.00672+0.000013×8.16372=0.00758, (81c)
CDVII4=CD4(0.0434º)=0.01155, (81d)
CDVII5=CD5(−16.002º)=0.01237+0.000008×16.0022=0.01442; (81e)
the drag coefficient weighted by the fraction of control surface area in Table 3:
CDVII=∑i=11CDVIIiSi/S0=0.00635, (82)
leads to a 66 counts cruise drag penalty of 66 counts (93a):
ΔCDVII≡CDVII−CD=0.00066, ΔCDVII/CD=0.116, (83a,b)
corresponding to 11.6% more drag (93b). This is shown in Table 10, which has a similar format to Table 9. It is seen from Table 10 that the sub-optimal strategy VII, of non-optimal use of optimal deflections from strategy V, improves over strategy V, but not over the better initial condition in strategy IV. This result is not surprising, since it has been shown that pitch trim is compatible with cruise drag reduction. The optimal strategy IV may not be the best in practice, if parameter limitations prevent from reaching the optimum; in the present case the deflections are limited to small values and the aerodynamic coefficients, and their first-and second-order derivatives, have a limited range of variation. The best results obtained so far were for a simpler non-optimal strategy of elevator deflection only. This result could be due to the linearized aerodynamic data used. In general, a simple non-optimal strategy can give good results if it incorporates the lessons learned from optimal strategies, as suggested by strategy IV, which is further improved next.
It is seen from Table 10 that for, the best optimal strategy IV and the sub-optimal strategy VII, the inboard and outboard control surfaces deflect in opposite directions. This is confirmed from Table 11 which shows that: (i) The aileron is the only surface with positive drag slope, meaning that a drag reduction corresponds to a small downward deflection; (ii) Conversely, for all the other surfaces, the drag slope is negative, and thus drag reduces for small upward deflections:
δ−≡δ5 < 0 < δ4=δ3=δ2=δ1≡δ+. (84)
The strategy VIII deflects downwards by the same amount δ− all surfaces with positive drag slope and deflects upward by the same amount δ+ all surfaces with negative drag slope, so as to reduce the drag for all control surfaces as shown by the example (94).
The condition of constant lift (60) for (94) is:
0.06177=αCLα+δ−C′L5+δ+∑i=14C′Li, (85a)
where may be substituted (47a) and (81a), viz.:
0.06177=0.086043αVIII−0.0049865δ−−0.042775δ+. (85b)
The condition (63) of pitch trim is:
0=CM=−0.00819−0.02259α+δ−C''M5+δ+∑i=14C′Mi, (86a)
where (83a) may by substituted:
0.00819+0.02259αVIII=0.0032δ−+0.019815δ+. (86b)
If the strategy VIII is adopted of keeping the AOA at the initial untrimmed value αVII=0.71789º=α0, in (26a) then deflections satisfy:
0.4934δVIII−+0.31223δVIII+=0, (87a)
0.003279δVIII−+0.019642δVIII+=0.024407; (87b)
it follows that the deflections:
δVIII−=+24.727º < 0 < −2.8825º=δVIII+, (88a,b)
are large for the ailerons (88a) and smaller for the other control surfaces (88b), and of opposite sign to what was required (84) for low drag. Thus trying to achieve trim at untrimmed AOA (26a), leads to unfavourable control surface deflections, in the sense that pitch trim is obtained by increasing drag. This is confirmed from the drag coefficients in Table 4 for the AOA (26a), viz.:
CDVIII 1≡CD1(δVIII+)=CD1(2.8825º)=0.00661, (89a)
CDVIII 2≡CD2(δVIII+)=CD2(2.8825º)=0.00587, (89b)
CDVIII 3≡CD3(δVIII+)=CD3(−2.8825º)=0.00587, (89c)
CDVIII 5≡CD5(δVIII−)=CD5(24.727º)+0.000008×24.7272=0.00856, (89e)
which lead by (43a) and Table 2 to a total drag coefficient (90a):
CDVIII =0.00551, ΔCDVIII=CDVIII−CD=0.00082, (90a-c)ΔCDVIII/CD=0.144,
corresponds to a large drag increase of 82 counts in (90b) or 14.4% in (90c), as shown in Table 12.
Allowing an AOA larger than the untrimmed value (26a), could lead to control surface deflections satisfying (84), and hence less control drag; this might overcome the extra drag due to increased angle of attack. Thus (85b) and (86b) are now solved for the control surface deflections δ±, without restriction on AOA, leading to:
δ−=37.287−17.495αVII, (91a)
δ+=5.7751+4.0510αVII. (91b)
For pitch trim without change of AOA αVIII = 0.71789° in (26a), from (91a,b) would follow (88a,b). No deflection of the inner control surfaces would be needed for the AOA:
δ+=0:α+=5.7751/4.0510=1.4256º, (92a)
δ−(α+)=037.287−17.495x1.4256º=12.346º; (92b)
the positive aileron deflection would lead to a drag increase though. Instead, the strategy IX is adopted, of choosing the AOA so that no aileron deflection is needed:
δIX−=0: αIX=37.287/17.495=2.1313º; (93a)
the deflection of the inner surfaces is now positive:
δIX+≡δ+(αIX)=−5.7751+4.0510x2.1313=2.8588º, (93b)
ensuring a drag reduction due to control surfaces. To confirm this, the drag coefficients are calculated from Table 4 for the AOA (93a), and listed in Table 8 for the deflections δ = ± 3°. Interpolating for the deflection δB+ in (93b), leads to the drag coefficients:
CDIX1=CD1(2.8588º)=0.00723, (94a)
CDIX2=CD2(2.8588º)=0.00816, (94b)
CDIX3=CD3(2.8588º)=0.00782, (94c)
CDIX4=CD4(2.8588º)=0.00794, (94d)
CDIX5=CD5(0.0º)=0.008685; (94e)
using (43a) and Table 3 leads to the total drag coefficients (95a):
CDIX=0.00787, ΔCDIX=CDIX−CD=0.00218, ΔCDIX/CD=0.383, (95a,b,c)
corresponding to a drag increase of 218 counts (95b) or 38.3% in (95c).Comparison of strategy IX, with strategy VIII in the Table 10 or (90a,b,c) shows that the increase in drag due to a larger AOA, outweighed the drag decrease due to favourable deflection of control surfaces for pitch trim.
As a final strategy X, may be considered a further elaboration of strategies VIII-IX, as follows: (i) Of all control surfaces with negative drag slope C'Di < 0, which should have positive deflection for drag reduction, choose that with larger |CDi| = -C'Di drag benefit, or stronger pitching moment slope C'Mi > 0, or better combination of both specified by larger product-C'Di C'Mi; (ii) Of all control surfaces with positive drag slope C'Di > 0, which should have negative deflection for reduced drag, choose that with largest drag slope coefficient C'Di, or largest pitching moment coefficient C'Mi or largest product of the two C'Di C'Mi > 0. From Table 11 it follows that only the aileron is in class (iii) Of positive slope of the drag coefficient; the remaining four surfaces have negative drag slope, with the centerbody elevator having the second largest |C'D0|, and largest C'Mi and largest product |C'Di| CMi. Thus:
δVII1≡δ+>0=δVIII2=δVIII3=δVIII4>0=δVIII5≡δ−; (96)
the centerbody elevator is given positive deflection, the aileron is given negative deflection, and all other surfaces are not deflected.
Using (81a, 47a) the condition of cruise lift equilibrium (46) for (96) is:
0.06177=0.086043 αVIII−0.011505 δ+−0.0049805 δ−; (97a)
the corresponding condition of pitch trim, viz. CM = 0 in (49) is:
0.00819 =0.02259 αVIII+0.00546 δ++0.003297 δ−. (97b)
Retaining the untrimmed value of AOA (21a) leads to:
0.011505δ++0.0049865δ−=0, (98a)
0.011546δ++0.003297δ−=−0.00803. (98b)
Solving this leads to the deflections
−8.6295º=δ− < 0 < δ+=3.7402º. (99a,b)
In the case of strategy X, keeping the untrimmed AOA, gives the favourable deflections to reduce drag due to control surface deflection. Together with the low drag for small AOA, this looks promising for low cruise drag. To confirm this the drag coefficients are calculated from Table 5 for the deflections (96) and (99a,b), viz.:
CDX1=CD1(3.7402º)=0.00483, (100a)
i=2,3,4: CDXi=CDi(0.0º)=0.00569, (100b)
CDX5=CD5(−8.6295º)=0.00478+0.00008×8.62952=0.00597; (100c)
using (43a) and Table 3 the corresponding total drag is (101a):
CDX=0.00554, ΔCDX=CDX−CD=−0.00015, ΔCDX/CD=−0.026, (101a-c)
which represents a drag reduction (101b) of 15 counts or 2.6% in (101c).
The discussion of pitch trim in cruise starts from the remark that it is compatible with (10.1) drag reduction relative to the untrimmed condition. The favourable combination of pitch trim and drag reduction has been found both for optimal and non-optimal pitch trim strategies, both in cases of constant and non-constant lift (10.2).
It should not be too surprising that pitch trim can, in the right circumstances, reduce cruise drag. This is illustrated in Figure 1a, which shows at the top a wing profile at AOA, with a stable pitch down moment. The projection perpendicular to the airflow L1 is roughly proportional to drag. In order to trim the 'flying wing' a pitch up moment is produced, by deflecting the elevator up in Figure 1b. This reduces the projection orthogonal to the free stream L2 < L1, leading to lower drag. This simple diagram suggests pitch trim can lead to a cruise drag reduction, for moderate upward deflection of trailing-edge control surfaces. The reality may be less simple than the illustration in Figure 1, and actually relates to the slopes of aerodynamic derivatives for control surfaces in Table 10. It is seen that the derivative of the drag coefficient with regard to control surface deflection is negative in all cases, except for the aileron when it is positive. Thus, the situation in Figure 1 of upward deflection leads to drag reduction for all control surfaces except the aileron. This difference may be due to wing twist or downwash variation along the wing span. In contrast, the derivative of the pitching moment is positive for all control surfaces. The preceding results show that moderate control surface deflections are enough to trim the aircraft if the untrimmed AOA is moderate. Thus, the trim strategy is to achieve lift balance at the smallest possible AOA, so that the nose down pitching moment is moderate, and can be countered by moderate upward deflection of trailing-edge control surfaces.
The possibility of cruise drag reduction by pitch trim was already evident in strategy II of elevator deflection upwards. Since this reduces lift, the trim AOA has to be increased relative to the untrimmed value, as shown is Table 9. The optimisation strategy V with strategy II as initial condition ultimately suggests strategy X which achieves the lowest cruise drag with pitch trim in Table 10. The intermediate step was as shown in the strategy IX of minimization drag by deflecting upward by the same amount (i) All control surfaces with negative drag derivative and deflecting downward by the same amount (ii) All surfaces with positive drag derivative. The two collective deflections can be chosen to satisfy lift equilibrium and pitch trim. The strategy X is derived from strategies VIII-IX by choosing in each group (i) and (ii) the most effective control surface in terms of (a) Drag reduction, (b) Pitch trim or (a) combination of both. The ten strategies I to X are summarized in List 1; the two strategies which combine pitch trim and lift balance with the largest drag reduction are singled-out in Table 12.
The Table 13 ranks the strategies according to the most desirable features, viz.: (i) Lowest cruise drag; (ii) Smallest AOA; (iii) Smallest control deflections in the sense of smallest average of modulus of largest and smallest deflection δ¯≡(|δmin|+|δmax|)/2; (iv) Linear aerodynamics works best for smallest AOA α and smallest mean deflection δ¯; (v) Simplicity of the computational algorithm, viz. least number of variables and least iterations. The strategy X has been shown to have five relative advantages: (i) It has the second lowest trimmed cruise drag of all eight strategies and twelve sub-strategies considered: (ii) The trimmed AOA is unchanged from the untrimmed value, and hence is the lowest and small; (iii) The pitching moment is thus small, and leads to moderate (5th smallest) control surface deflections; (iv) The small AOA and control surface deflections imply that linear aerodynamics can give reliable results (3rd best); (v) The strategy is easy to apply (3rd best). The other pitch trim strategies considered in Table 13 do not have one or more of these attributes.
The reasons why strategy X is quite effective may be explained as follows: (i) The ailerons are the only control surfaces for which a positive deflection reduces drag; (ii) For the remaining control surfaces, a negative deflection reduces drag; (iii) Of the latter, the centerbody elevator provides the best combination of large pitching moment for trim and relatively large drag reduction; (iv) Positive deflection of the elevator and negative deflection of the aileron can balance each other for the effect on lift; (v) Thus the untrimmed AOA is unchanged by pitch trim, and there is no drag penalty; (vi) The unchanged drag due to AOA, and reduced drag to control surface deflection, yield a cruise drag reduction due to trim; (vii) The small AOA and moderate control surface deflections give confidence on the predictions of linear aerodynamics; (viii) This pitch trim strategy is simple to apply, since it involves only two control surface deflections, and no change in AOA.
The strategy II of deflection of centerbody elevator alone gave a slightly better cruise drag reduction due to trim (minus 18 drag counts) than the strategy X of opposite deflection of the aileron also (minus 15 drag counts). However, strategy X allows an unchanged, hence small AOA of 0.7° whereas strategy II increases AOA by a factor of three to 2.1°. For this reason, strategy X is considered preferable. It is also possible to explain why other strategies are less successful or even unfavourable. For example, the strategy I of equal deflection of all control surfaces, does not exploit the drag reduction due to the deflection of the ailerons opposite to the other control surfaces; thus it leads to high cruise drag due to trim. The same applies to strategy III, of preferential deflection of inner control surfaces, which do not use aileron deflection to reduce cruise drag and AOA. The methods which give intermediate results, like the early iterations, not fully converged, of the optimal strategies V, exhibit deflection of outer control surfaces opposite to inner control surfaces, but do not exploit them fully. It is clear from the present account that the choice of the best pitch trim strategy is guided by the aerodynamic derivatives of the particular aircraft configuration and the control surfaces available for pitch trim.
The preceding are examples of the qualitative explanation of the results obtained in Table 13 for each of the eight pitch trim strategies in List 1, viz.: (i) The strategy I of equal deflection of all control surfaces, causes the ailerons to be deflected as the inner control surfaces, thus increasing drag and trimmed AOA; (ii) The strategy II of deflection of centerbody elevator alone, is effective at reducing cruise drag, but also reduces lift, and requires increased AOA for trim, with no AOA reduction provided by the undeflected aileron; (iii) The strategy III of preferential deflection of inner control surfaces, also does not deflect the ailerons, and thus leads to an increased AOA for trim, which dominates cruise drag; (iv) The optimal strategies IV and V evolve from any initial condition in the right direction, deflecting in opposition inner and outer control surfaces, but the convergence is slow, and a few iterations do not reach the best results; (v) The strategy VI of equal contribution to lift from all control surfaces, does not use the ailerons to reduce lift and trimmed AOA, and hence also increases cruise drag; (vi) The sub-optimal strategy VII of multiplying the optimal deflection by a constant, determined from cruise lift equilibrium, improves on strategy for the same number of iterations, because it accelerates convergence; (vii) The strategies VIII-IX of deflecting the ailerons opposite to all other inner control surfaces is correct, but is spoiled by poor lift balance of the single aileron versus all the other inner control surfaces leading to a large trimmed AOA and hence increase in cruise drag; (viii) The strategy X works out best because it deflects the ailerons in opposition to centerbody elevator alone, achieving lift equilibrium with unchanged small AOA, and requiring moderate control surface deflections each of which contributes to cruise drag reduction. Thus a discerning use of simple pitch trim strategies turns out to be more effective than brute force recourse to computationally complex and slow converging optimal procedures.
In a conventional aircraft configuration the elevator or all-flying tail is the only effective pitch control surface. The elevator must be large enough to provide adequate pitch control authority at low speed, when large deflections are possible. At high speed, when aerodynamic and/or aeroelastic effects may limit deflections, the higher dynamic pressure may compensate, and not be a “sizing” criterion requiring large area, with aerodynamic, aeroelastic, weight or actuator problems. In any case the compromise in elevator sizing is straightforward once the required pitching moment and possible lever arm are known. The BWB or FW configuration has the potential to provide pitch trim with several combinations of control surface deflections. There are simple choices with obvious consequences but it is not clear what is the best compromise, e.g.: (i) The strategy II of deflecting the centerbody elevator only can lead to a large deflection; (ii) Using all control surfaces with the same deflection (strategy I) or same lift (strategy VI) does not exploit the most effective controls; (iii) Deflecting preferentially the inner control surfaces with an aeroelastic limit (strategy III) may not be the best compromise among the preceding choices.
The optimization methods may or may not do better than the simpler non-optimal strategies depending on whether the drag has a single global minimum or several local minima. If there are several local minima convergence will occur to that closest to the initial condition. The influence of the initial condition on the result of the optimization may not be obvious. It may be difficult to know if a local minimum is lower that neighbouring ones or is the global minimum. Also the optimization algorithm may require data, like second-order control derivatives which is not accurate, leading to a possible dilemma: (i) Too few iterations do not converge closely enough to the optimum; (ii) Too many iterations accumulate too large an error around the optimum. It may also occur that the local minimum is outside the acceptable range of control deflections, and the optimum is not a local minimum but rather the best value on a boundary. Thus, optimal methods have strengths and weakness which may be different from simple methods and a combination of the two may be more effective than either of them.
A careful examination of control derivatives is essential to understand how an optimization procedure is going to work, or what are the consequences of a simple strategy. In the case of pitch trim for cruise the starting ideas seems to be to select the most effective control surfaces, in the sense of combining the largest possible pitching moment and drag reduction with small deflections; it is also important to achieve lift balance, so that it is not necessary to increase the AOA for trimming purposes. A larger AOA increases drag, may require a larger pitching moment, and thus works against having selected trim with the “most effective” control surfaces; thus the most effective control surfaces should be chosen to (i) Provide pitching moment, (ii) Reduce drag and (iii) Balance lift. In the present example it was possible to achieve pitch trim while reducing drag relative to an untrimmed flight condition, also not changing AOA; the latter contributes to smaller control deflections, avoiding undesirable aeroelastic effects, and staying within the limits of linear or weakly non-linear aerodynamics and control methods.
This work was supported by FCT (Foundation for Science and Technology), through IDMEC (Institute of Mechanical Engineering), under LAETA, project UID/EMS/50022/2019. This work was started under the NACRE project of the European Union and benefited from comments of other partners in this activity.