# Experimental Study of Aircraft Achieving Dutch Roll Mode Stability without Weathercock Stability

Authors: Jingcheng Fu; Jun Huang; Lei Song; Daqing Yang
Journal: International Journal of Aerospace Engineering
 (2020)

---

## Abstract

Weathercock stability is usually considered essential to achieve normal flight, while the Dutch roll mode stability can still be achieved without weathercock stability which has been algebraically proved. This paper proposed a flight experiment to investigate the characteristics of an airplane with Dutch roll mode stability but no weathercock stability. Firstly, the algebraic analysis based on a standard lateral-directional mode approximation was made to demonstrate the effect of yawing stability derivativeCnβ on the Dutch roll mode characteristics. The flight experiment was organized after that using a model glider which was modified to have zero Cnβ but with marginal change on Cyβ. The convergence of Dutch roll mode in flight meets the algebraic and numerical analysis as expected. However, the difference of handling characteristics between the original and modified configurations indicates some other roles the weathercock stability plays in flight as well as some limitations of utilizing mode criterion in flight quality analysis.


---

## Body

## 1. Introduction

Weathercock stability, also known as directional static stability or yaw stiffness, is the ability of the aircraft to turn in to wind in order to maintain directional equilibrium. The static stability principle suggests that weathercock stability is essential for a stable airplane [1]. The derivative Cnβ that corresponds to the yawing moment due to sideslip is designed to be a positive value for most conventional airplanes as a result of weathercock stability assumption. The main contribution to Cnβ comes from vertical stabilizer. The fuselage and wing are also major contributors to the value of this derivative [2].From the perspective of lateral-directional stability study in flight dynamics, the perturbed state problems were first introduced in 1911 to investigate whether the response to the external disturbance was acceptable [3, 4]. After the airplane stability quartic in terms of its polynomial coefficients was rooted by Bairstow, the Routh stability criterion was the first used to study the dynamic stability behavior of an airplane; then, the term Dutch roll was used by Hunsaker to identify the oscillatory roots of the lateral-directional characteristic equations [5–8]. For the lateral-directional dynamic stability, a hypothetic small airplane with various modifications of the fin area and dihedral setting was used to investigate the motions caused by the lateral disturbance before 1938 [9]. As demonstrated in Ref. [9], with different combinations of Cnβ and Clβ, types of spiral and oscillatory divergence limits on these two derivatives were found.To make specifications of the flight quality for new military airplanes, MIL-8785B (ASG), which includes specifications of lateral-directional mode characteristics, was published in 1969 based on the data from typical existing aircrafts [10]. It was continuously validated and revised by the US military, institutes, and manufacturers; then, the final revision of MIL-F-8785C was issued in 1980 [11]. Later in 1995, this military specification was expanded and revised to MIL-STD-1797A and designated as a handbook for guidance only, where the flight quality specifications for the lateral-directional mode were involved identically [12].For the lateral-directional stability and control study from the designer’s perspective, in the preliminary configuration design stage, the tail-volume method is often used in the first-cut estimation of vertical stabilizer sizing. According to the statistics of the tail-volume range for specific categories of existing aircrafts, the appropriate tail volumes are determined, and the area of the vertical stabilizer can be obtained [2, 13]. As suggested in Ref. [14], for most airplanes, the derivative Cnβ has a lower limit of 0.001-0.0025.The vertical stabilizer appears indispensable for an aircraft. However, in aviation history before the stability augment system was invented, several aircrafts with flying wing configurations such as Ho.229 and N-9M appeared which have no vertical stabilizer to provide enough directional static stability. On the other hand, it has been proven that the lateral-directional stability of this configuration can be achieved by adjusting the spanwise dihedral layout without installing a stability augment system [15, 16]. Instead of focusing on the required lateral-directional stability derivative Clβ or Cnβ, a criterion of dynamic stability was included in the aerodynamic design for flying wing, which considered the Dutch roll mode characteristics as the critical element for flight qualities [17]. Moreover, a unified algebraic approach was developed to approximate aircraft lateral-directional modes and predict departure susceptibility [18]. And in Ref. [19], the relationship between vertical stabilizer and lateral-directional stability was studied and illustrated by some representative experiments. Above all, the study on flying wing aircrafts and the Dutch roll mode stability characteristics illustrates the possibility of normal flight with no sufficient weathercock stability. But the role that weathercock stability played in normal flight is still worth investigation.This paper is aimed at demonstrating that the weathercock stability is not essential to achieve Dutch roll mode stability for a conventional layout airplane but has its influence on flight quality. The demonstration is made from three aspects: algebraic analysis, numerical study, and flight experiments. The flight test result may indicate some other significance of the derivativeCnβ in flight quality which are not demonstrated in Dutch roll mode flight quality specifications.


## 2. The Essential Factor for Dutch Roll Mode Stability

To demonstrate the weathercock stability unessential and find the crucial factors to achieve the Dutch roll mode stability, the classic lateral-directional mode literal approximation was considered for theoretical research. In fact, such literal expression involves aerodynamic and inertial data and it provides valuable insight into the connection between the relative values of certain key derivatives and behavior of aircraft motion [18].Commonly, the derivation of the approximation contains omissions of relatively unimportant terms, to make the relationship between the cause and effect apparent. But for the purpose of studying the crucial factor for Dutch roll stability, some certain omissions need examination again in the developing of approximation. The simplified model for the approximation is obtained by separating slow mode from fast modes. That is, for lateral-directional motions, the fourth-order equation approximation is made by separating the spiral mode (first-order system) from the roll and Dutch roll mode (third-order system).To develop the Dutch roll mode natural frequency approximation, we designate the stability axis reference frame and assume a condition of level flight; the dynamic matrix from the lateral-directional linearization small disturbance equations is given by [20]
(1)Alat=Yβ/V00−1g/V0Lβ′Lp′Lr′0Nβ′Np′Nr′00100,where Li′ and Ni′ are standard primed derivatives, V0 refers to trimmed level flight velocity, and g refers the gravitational acceleration. Consider the standard factorization of the characteristic polynomial for lateral-directional modes,
(2)Plats=s+λSs+λRs2+2ζDRωDRs+ωDR2.The subscriptsS,R,DR refer to the spiral, roll, and Dutch roll modes, respectively. Letting
(3)Plats=s4+Bs3+Cs2+Ds+E,one has
(4)B=λR+λS+2ζDRωDR,C=ωDR2+2ζDRωDRλR+λS+λRλS,D=λR+λSωDR2+2ζDRωDRλRλS,E=λRλSωDR2.The fourth-order polynomial coefficients can be determined by derivatives in matrixAlat. Consider the four coefficients in the characteristic polynomial for the lateral-directional system matrix, they are
(5)B=−YβV0−Nr′−Lp′,C=Nβ′+Lp′YβV0+Nr′+YβNr′V0−Np′Lr′,D=Np′YβLr′V0+Lβ′−Lp′Nr′YβV0+Nβ′−Lβ′gV0,E=Nr′Lβ′−Nβ′Lr′gV0.First to consider first-order spiral mode. The time constantλS is customary assumed to be λS≪minωDR,λR. From that,
(6)D≅ωDR2λR.Consider the classic roll approximation which is given by(7)λR=−Lp′that is found by considering a single DOF roll subsidence model [20]. Accordingly, as long as the relationship λS≪minωDR,λR can be satisfied, from equations (5), (6), and (7), the classic Dutch roll natural frequency approximation yields
(8)ωDR2=Nβ′+YβNr′V0−Np′Lp′Lβ′+YβLr′V0+Lβ′Lp′gV0.In order to demonstrate that weathercock stability is not essential to Dutch roll stability and investigate the crucial factor, we move the focus on the right terms of equation (8) for the result. As can be seen in the right polynomial, the first three terms are positive. For most of the aircraft configurations and flight conditions, weathercock stability aerodynamic coefficient Nβ′ is relatively big among the right four terms, where the third term on the right side is usually very small. Therefore, when simply setting Nβ′ to zero, the second term on the right always remains positive to guarantee the residual Dutch roll stability, and this can be achieved by a relatively big Yβ.The vertical tail has been widely recognized as essential for most fixed-wing aircrafts to achieve lateral-directional stability, particularly in the Dutch roll mode stability. From the perspective of reducing sideslip and keeping steady, the yawing moment and side force produced by vertical tail both make effort. The related two aerodynamic coefficients areNβ′ and Yβ. This amounts to saying that although the weathercock stability can affect the Dutch roll mode, it is not essential for the Dutch roll stability, and the side-force derivative Yβ should be another crucial factor. In fact, with respect to the flight quality specifications, zero weathercock stability but with sufficient Nβ′, the lateral-directional mode characteristics can satisfy the flight quality specification.


## 3. Flight Experiment Verification with Zero Weathercock Stability

To verify the algebraic approximation result, a model glider with the conventional configuration was selected as the flight test object with reduced weathercock stability but maintained side force stability. The glider has a wingspan of 1.68 m, a length of 0.87 m, and a mass of 0.714 kg, and it is powered by an electric motor-driven propeller. The inertia moment parameters relative to lateral-directional motions can be found in Table1. To have zero Nβ′ but a marginal effect on Yβ, the aircraft was modified by reducing the area of the vertical stabilizer to less than one-fourth and attaching a fin that was almost three times larger in front of the center of gravity.Table 1
Inertia moments involving the lateral-directional stability for the modified configuration.
Parameters
Valueskg⋅m2
Ix
0.046
Iz
0.074
Ixz
0.001The original and modified aircrafts are shown in Figures1 and 2. The typical level flight condition for this glider is 18 m/s at 100 m in altitude.Figure 1
3-dimensional view of the original configuration.Figure 2
Contrast between the modified and original configurations from the side view.The aerodynamic parameters were obtained from a vortex lattice method program [21]. The vortex layout of the modified configuration is shown in Figure 3. In Table 2, several representative aerodynamic derivatives involving lateral-directional stability are listed to show the difference between original and modified configurations. Then, three groups of eigenvalues of system matrix Alat were obtained by substituting the aerodynamic derivatives and inertial parameters. Here, the standard primed derivatives Lβ′ and Nβ′ as well as other relevant derivatives are given by [20]
(9)Li′=Li+Ixz/IzNiIx−Ixz2/Iz,Ni′=Ni+Ixz/IxLiIz−Ixz2/Ix,Yβ=Cyβq∗SmV0,Lβ=Clβq∗Sb,Lp=Clpq∗Sb,Lr=Clrq∗Sbb2V0,Nβ=Cnβq∗Sb,Np=Cnpq∗Sbb2V0,Nr=Cnrq∗Sbb2V0.Figure 3
The aerodynamic calculation vortex lattice layout.Table 2
Stability derivatives of lateral-directional mode characteristics for the original and modified glider configurations.
Lateral-directional stability derivative
Original configuration
Modified configuration
Cnβ
0.0593
0.0001
Cyβ
-0.1835
-0.1807
Cnr
-0.0477
-0.0188The calculated three groups of lateral-directional mode eigenvalues and the corresponding parameters are shown in Table3. According to the Level 1 criterion in MIL-F-8785C, the time constant of roll mode τR should be less than 1.4 s, and the time τdouble to double the amplitude of the spiral mode to 20 degrees should be longer than 20 s. For the Dutch roll mode, its damping ratio ζDR, natural frequency ωDR, and their product ζDRωDR should be larger than 0.08, 0.4, and 0.15, respectively. All three lateral-directional mode characteristics of the modified configuration including Dutch roll mode satisfy the Level 1 specification.Table 3
Lateral-directional mode characteristics of the modified and original configurations.
Dutch roll mode
Configuration
λDR
ζDR
ζDRωDR
ωDR
Modified
−0.3772±1.535i
0.239
0.377
1.581
Original
−0.773±3.748i
0.202
0.773
3.826
Spiral mode
Roll mode
Configuration
λS
τdouble
λR
τR
Modified
-0.0368
N/A
-1.959
0.0557 s
Original
0.181
3.823
-18.120
0.0552 sAs can be seen in Table3, the Dutch roll mode frequency of modified configuration is smaller compared with that of original configuration. This can be attributed to the significantly reduced Cnβ. From equation (8), the reduction of Cnβ is equivalent to the decrease of Nβ′ on the right side, thence caused the oscillation frequency ωDR decrease. The trend of frequency between these configurations is reasonable from this perspective.Before the flight test, a lateral-directional motion simulation of the modified configuration is made andβ and ϕ after a 10° sideslip disturbance are shown in Figure 4.Figure 4
β and ϕ of the modified glider after 10° sideslip disturbance in lateral-directional simulation.The glider with the modified configuration was flight-tested by remote control with no stability augmentation. Several rudder doublet excitations were conducted, and the convergence of the Dutch roll motion demonstrates the consistency between the flight test and the numerical result. A representative sequential shooting of the aircraft response to the rudder doublet is shown in Figure5. This convergence indicates that the aircraft can achieve Dutch roll mode stability without the weathercock stability.Figure 5
Configuration modified aircraft response to the rudder doublet Dutch roll mode excitation.Two noticeable outcomes besides appeared in the flight test according to the controller’s feedback. Firstly, the modified glider requires additional rudder control input in turning compared with the original configuration. For conventional aircrafts, sideslip arises but within a limit at a relatively small value with no rudder control conducted in turning maneuvers. But for the modified glider, obvious sideslip was observed when conducting a turning maneuver if rudder control was not involved, and the nose heading did not change much until rudder control input is towards the same direction. This can be explained as no yawing moment to reduce the sideslip because of the lack of weathercock stability. The simulation of the sideslip and rolling angle and the rolling and yawing angular velocity after the glider encounter a 10° right rolling angle disturbance is shown in Figure6 to make a comparison between original and modified gliders on this issue. As can be seen in Figure 6, for the original aircraft, the 10° rolling angle soon results in a relatively big yawing angle velocity up to 15°/s and goes on increasing to the same direction but the sideslip angle remains under 2°. On the contrary, the modified one is in big sideslip angle up to 8 at the beginning but yawing angular velocity increases slowly and remains under 10°/s.Figure 6
Lateral-directional motions of modified and original gliders after 10° rolling angle disturbance.Secondly, the modified glider was easier to enter spins than the original configuration. The flat spin occurred twice in climbing soon after the aircraft took off. During the climbing and turn at a slow airspeed, the sideslip of the glider was difficult to realize and eliminate in time. Thus, with the increase in sideslip angle and big attack angle in climbing at low airspeed, one of the wings stalled; then, the glider entered a spin and quickly developed to a flat spin because there was no weathercock stability to reduce sideslip, as shown in Figure7. From the sequential photos, obvious left sideslip during the right flat spin can be observed and this was difficult to recover under conventional spin recover operation, that is, same aileron and opposite rudder. The reduced rudder should cause the result.Figure 7
Flat spin in the takeoff and climb phase.


## 4. Discussion on the Phenomena and Relevant Inspirations

Both algebraic analysis and flight experiment illustrate that the weathercock stability is not essential to the Dutch roll mode stability. However, in flight appeared the phenomena that additional rudder control requirement in turning maneuvers and easy occurrence of spins reflect the important roles the weathercock stability plays in flight. The convergence of the Dutch roll mode demonstrates that, within the flight quality judgment criterion based on the lateral-directional mode characteristic specifications defined in MIL-F-8785C or MIL-STD-1797A, the Dutch roll dynamic stability can be achieved with zero weathercock stability. But those two phenomena can bring bad effect in flight and they cannot be avoided by just satisfying the mode criterions in MIL-8785C and MIL-STD-1797A.In turning maneuvers, the additional rudder control requirement of the modified glider reflects thatCnβ plays a more important role in limiting the sideslip angle during the maneuvers rather than in the Dutch roll mode stability. In addition, with less or even zero weathercock stability, the characteristics of the rudder control requirement in coordinated turns should also be investigated and set a reference standard.Meanwhile, the lack of weathercock stability makes it more difficult to recover from the beginning of spin. Sufficient weathercock stability can restrict the sideslip from going serious when spin happens; thus, the aerodynamic state on wings will be improved. As can be inferred this will help to prevent the spin from developing and recover to normal flight.While neither MIL-F-8785C nor MIL-STD-1797A involves such relevant specifications on weathercock stability, the phenomena found in flight experiment cannot be demonstrated basically by satisfying the mode criterions according to flight quality specification. From the view of the aircraft design, the sufficient weathercock stability has its meanings, regardless by the means of aerodynamics or flight control system. Besides the way like arranging appropriate vertical stabilizer, the split spoilers on the wing tips controlled by flight control system according to the simultaneously detected sideslip angle can also make the aircraft possess sufficient weathercock stability when needed. The method from the view of flight control system has its special signification for flying wing configuration and this is worth further detailed research.The impact of gust wind on flight has been investigated from different aspects like overload control. For the impact of lateral wind gust, it is closely related to the residual lateral-directional stability characteristics of the aircraft. To investigate the influence of zeroedCnβ on the response of aircraft to the lateral gust wind, the numerical method was used. The simulation is based on the lateral-directional linearization model of the original and modified gliders. The lateral gust wind was simplified to a disturbance in sideslip angle at 10°. By comparing the motions after the sideslip angle disturbance between different models, the influence of Cnβ on the response to lateral gust wind can be studied. The four typical lateral-directional motion parameters of the response motions are shown in Figure 8.Figure 8
Lateral-directional motions of modified and original gliders after 10° sideslip angle disturbance.From Figure8, it apparently shows that the modified glider with zeroed Cnβ has less directional and rolling angular rate r and p response after the 10° sideslip disturbance with zero control input compared with those of the original glider. Both configurations have sideslip angle convergence, though the modified glider expresses lower oscillation frequency which is consistent with the previous numerical and experimental study. From this perspective, the zeroed Cnβ can help to reduce the impact of lateral gust wind especially on directional motions. For flying wing and other aircraft configurations with less weathercock stability, they can be less affected by the impact of lateral gust wind compared with bigger weathercock stability configurations.


## 5. Conclusion

In this paper, the algebraic analysis and flight experiment are used to illustrate that the weathercock stability is unessential to the Dutch roll mode stability but plays some roles in normal flight. First, an algebraic approach was developed to approximate Dutch roll mode natural frequency and the literal expression demonstrates that the side force derivativeCyβ can also guarantee the residual Dutch roll stability at zeroed Cnβ. Then, a flight experiment was organized based on a model glider, which was modified to have zero Cnβ but with a marginal change on Cyβ. The flight test showed the convergence of the Dutch roll after procedural mode excitation and some other phenomena that are thought as a result of lack of weathercock stability. The phenomena include the extra requirement of rudder input in the turning maneuvers compared with original configuration and easy occurrence of spin. The relationship between weathercock stability and the impact of lateral gust wind is also investigated based on the simulation of the two glider configurations. The result shows that with less weathercock stability the impact of lateral gust wind can be weakened.The convergence of Dutch roll in the flight experiment showed the consistence with the algebraic and numerical analysis. With the lateral-directional mode criterions in specifications like MIL-F-8785C, they can be satisfied without the weathercock stability. From the view of flight experiment, the side force derivativeCyβ also plays an important role in Dutch roll mode stability.Even though a residual Dutch roll mode stability can be achieved with zeroCnβ, the weathercock stability still shows its roles played in flight. The extra rudder control input in turning maneuvers and easier occurrence of spin demonstrated its effects. These problems cannot be demonstrated or avoided by merely satisfaction of mode criterion in flight quality specifications. Hence, the appropriate weathercock stability has its meanings no matter by means of aerodynamics or directional stability augmentation. Moreover, the reference criterion of weathercock stability in flight quality is worth further investigation with the development of flying wing and other unconventional aircraft design.



---

*Source: 8971275-2020-02-11.xml*