The oscillation susceptibility of the ADMIRE aircraft along the path of longitudinal flight equilibriums is analyzed numerically in the general and in a simplified flight model. More precisely, the longitudinal flight equilibriums, the stability of these equilibriums, and the existence of bifurcations along the path of these equilibriums are researched in both models. Maneuvers and appropriate piloting tasks for the touch-down moment are simulated in both models. The computed results obtained in the models are compared in order to see if the movement concerning the landing phase computed in the simplified model is similar to that computed in the general model. The similarity we find is not a proof of the structural stability of the simplified system, what as far we know never been made, but can increase the confidence that the simplified system correctly describes the real phenomenon.
1. Introduction
Frequently, we describe the evolution of real phenomena by systems of ordinary differential equations. These systems express physical laws and geometrical connections, and often they are obtained by neglecting some influences and quantities, which are assumed insignificant with respect to the others. If the obtained simplified system correctly describes the real phenomenon, then it has to be topologically equivalent to the system in which the small influences and quantities (which have been neglected) are also included. Furthermore, the simplified system has to be structurally stable. Therefore, when a simplified model of a real phenomenon is build up, it is desirable to verify the structural stability of the system.
Interest in oscillation susceptibility of an aircraft is generated by crashes of high-performance fighter airplanes, such as the YF-22A and B-2, due to the oscillations that were not predicted during the aircraft development [1]. Flying qualities and oscillation prediction are based on linear methods and their quasilinear extensions [2]. These analyses cannot, in general, predict the presence or the absence of oscillations, because of the large variety of nonlinear interactions that have been identified as factors contributing to oscillations. Some of these factors include pilot behavioral transitions, actuator rate limiting [3–5], and changes in aircraft dynamics caused by transitions in operating conditions [6], gain scheduling, and switching [7]. The oscillation susceptibility analysis in a nonlinear model involves the computation of nonlinear phenomena including bifurcations (Hopf or fold bifurcations) that leads sometimes to large changes in the stability of the aircraft.
Oscillation susceptibility analysis means the evaluation of the oscillation potential of a given aircraft: identify characteristics of the pilot-aircraft interaction that may result in oscillation, demonstrate the potential for oscillation by analysis and simulations using appropriate piloting tasks and test maneuvers, distinguish aircraft configurations that are less susceptible to oscillations from those that have high oscillations potential, and suggest to reduce and/or eliminate oscillation susceptibility [1].
As an example in [1] the X-15 aircraft oscillation caused by the rate limiting and an F/A-18 aircraft oscillation caused by nonlinear category III triggers are presented. The limit cycle amplitudes are computed for the longitudinal flight equations of motion and large jump in limit cycle amplitude indicating a significant change in the vehicle stability is revealed.
Our aim in this paper is to analyze numerically the oscillation susceptibility of the ADMIRE aircraft in a longitudinal flight in a quasilinear (simplified) and a nonlinear (general) flight model in landing phase, when the Automatic flight Control System (AFCS) is decoupled. The equations governing such a flight and the conditions which assure the existence of such a flight are presented. The equilibriums flights are analyzed numerically, from the point of view of bifurcations which can appear due to the changes of the elevator deflection. Appropriate piloting tasks and maneuvers for the touch-down moment are established. The behavior of the aircraft is simulated in both models. The computed results obtained in the models are compared in order to see if the simplified model correctly describes the real flight. This is not a proof of the structural stability of the simplified system but can increase the trust that the simplified system correctly describes the real phenomenon.
2. The General Nonlinear Model
The system of differential equations [8, 9], which describes the motion around the center of gravity of a rigid aircraft, with respect to an xyz body-axis system, where xz is the plane of symmetry, is
V∘V·cosα·cosβ-β∘·cosα·sinβ-α∘·sinα·cosβ=r·sinβ-q·sinα·cosβ-gV·sinθ+Xm·V,V∘V·sinβ+β∘·cosβ=p·sinα·cosβ-r·cosα·cosβ+gV·sinφ·cosθ+Ym·V,V∘V·sinα·cosβ-β∘·sinα·sinβ+α∘·cosα·cosβ=-p·sinβ+q·cosα·cosβ+gV·cosφ·cosθ+Zm·V,Ix·p∘-Ixz·r∘=(Iy-Iz)·q·r+Ixz·p·q+L,Iy·q∘=(Iz-Ix)·p·r-Ixz·(p2-r2)+M,Iz·r∘-Ixz·p∘=(Ix-Iy)·p·q-Ixz·q·r+N,φ∘=p+q·sinφ·tanθ+r·cosφ·tanθ,θ∘=q·cosφ-r·sinφ,ψ∘=q·sinφ+r·cosφcosθ.
The state parameters of this system are forward velocity V, angle of attack α, sideslip angle β, roll rate p, pitch rate q, yaw rate r, Euler roll angle φ, Euler pitch angle θ, and Euler yaw angle ψ. The constants Ix, Iy, and Iz are the moments of inertia about the x-, y-, and z-axis, respectively, Ixz is the product of inertia, g is the gravitational acceleration, and m is the mass of the vehicle. The aero dynamical forces X,Y,Z and moments L,M,N are functions of the state parameters and the control parameters: δa is the aileron deflection, δe is the elevator deflection, and δr is the rudder deflection (the body flap, speed break, δc, δca are available as additional controls but, for simplicity, they are set to 0 in the analysis to follow).
Definition 2.1.
A flight with constant forward velocity V is defined as a flight for which V=const (i.e., V∘=0).
Proposition 2.2.
In a flight with constant forward velocity Vthe following equalities hold:
-β∘·cosα·sinβ-α∘·sinα·cosβ=r·sinβ-q·sinα·cosβ-gV·sinθ+Xm·V,β∘·cosβ=p·sinα·cosβ-r·cosα·cosβ+gV·sinφ·cosθ+Ym·V,-β∘·sinα·sinβ+α∘·cosα·cosβ=-p·sinβ+q·cosα·cosβ+gV·cosφ·cosθ+Zm·V,
Proof.
Replacing V∘ by 0 in the system (2.1), then (2.2) is obtained.
Proposition 2.3.
If in a flight with constant forward velocity V one has β≡(2n+1)·π/2, then the following equalities hold:
(-1)n·r-gV·sinθ+Xm·V≡0,gV·sinφ·cosθ+Ym·V≡0,(-1)n+1·p+gV·cosφ·cosθ+Zm·V≡0.
Proof.
Replacing β∘=0 and β≡(2n+1)·π/2 in (2.2), then (2.3) is obtained.
Proposition 2.4.
If in a flight with constant forward velocity V one has β≠(2n+1)·π/2, then the following equality holds:
g[sinβ·sinφ·cosθ-cosα·cosβ·sinθ+sinα·cosβ·cosφ·cosθ]+Ym·sinβ+Xm·cosα·cosβ+Zm·sinα·cosβ≡0.
Proof.
Equation (2.4) is the solvability (compatibility) condition of system (2.2) when β≠(2n+1)·π/2.
Proposition 2.5.
If β≠(2n+1)·π/2 and equality (2.4) holds, then the system (2.2) can be solved with respect to α∘, β∘, obtaining the explicit system of differential equations, which describes the motion of the aircraft in a flight, with constant forward velocity V:
α∘=q-p·cosα·tanβ-r·sinα·tanβ+gV·cosβ·[cosφ·cosθ·cosα+sinθ·sinα]+1cosβ·[Zm·V·cosα-Xm·V·sinα],β∘=p·sinα-r·cosα+1cosβ·gV·sinφ·cosθ+1cosβ·Ym·V,Ix·p∘-Ixz·r∘=(Iy-Iz)·q·r+Ixz·p·q+L,Iy·q∘=(Iz-Ix)·p·r-Ixz·(p2-r2)+M,Iz·r∘-Ixz·p∘=(Ix-Iy)·p·q-Ixz·q·r+N,φ∘=p+q·sinφ·tanθ+r·cosφ·tanθ,θ∘=q·cosφ-r·sinφ,ψ∘=q·sinφ+r·cosφcosθ.
Proof.
System (2.5) is obtained solving system (2.2) with respect to α∘, β∘ and replacing in system (2.1), then (2.1)1, (2.1)2, and (2.1)3 with the above obtained α∘ and β∘.
Definition 2.6.
A longitudinal flight is defined as a flight for which
β≡p≡r≡φ≡ψ≡0,δa=δr=0.
Proposition 2.7.
A longitudinal flight is possible if and only if Y=L=N=0 for β=p=r=φ=ψ=0 and δa=δr=0.
Proof.
This result is obtained from (2.1) taking into account Definition 2.6.
Proposition 2.8.
The explicit system of differential equations which describes the motion of the aircraft in a longitudinal flight is
V∘=g·sin(α-θ)+Xm·cosα+Zm·sinα,α∘=q+gV·cos(θ-α)-Xm·V·sinα+Zm·V·cosα,q∘=MIy,θ∘=q.
Proof.
This result is obtained from (2.1) taking into account Definition 2.6.
Remark 2.9.
In system (2.7) X,Z,M depend only on α,q,θ, and δe. These dependences are obtained replacing in the general expression of the aerodynamic forces and moments: β=p=r=φ=ψ=0 and δa=δr=0.
Proposition 2.10.
The explicit system of differential equations which describes the motion of the aircraft in a longitudinal flight with constant forward velocity V is
α∘=q+gV·cos(θ-α)-Xm·V·sinα+Zm·V·cosα,q∘=MIy,θ∘=q.
Proof.
This system is obtained from (2.7) taking into account V∘=0.
Proposition 2.11.
A longitudinal flight with constant forward velocity is possible if the following equalities hold:
Y=L=N=0forβ=p=r=φ=ψ=0,δa=δr=0,g·sin(α-θ)+Xm·cosα+Zm·sinα=0.
Proof.
This result is obtained from Proposition 2.7 and system (2.7), taking into account the fact that V∘ is equal to zero.
Remark 2.12.
Notice that in (2.8) and (2.10) X,Z,M depend on α,q,θ,δe, and V. Taking into account (2.10), the system (2.8) can be written as
α∘=q+gV·cos(θ-α)-gV·sin(θ-α)·tanα+Zm·V·1cosα,q∘=MIy,θ∘=q.
Remark 2.13.
In system (2.11), the functions Z=Z(α,q,θ;δe,V) and M=M(α,q,θ;δe,V) are considered known; δe and V are parameters.
The system (2.11) describes the motion around the center of gravity of an aircraft in a longitudinal flight with constant forward velocity V and defines the general nonlinear model.
3. The Simplified Model of the ADMIRE Aircraft
The ADMIRE aircraft is an Aero Data Model In a Research Environment. To describe the flight of this vehicle with constant forward velocity V, the following explicit system of differential equations is employed:
System (3.1) can be obtained from (2.5) substituting the general aero dynamical forces and moments with those corresponding to the ADMIRE aircraft [10], assuming that α and β are small and making the following approximations:
The simplified system which governs the longitudinal flight with constant forward velocity V of the ADMIRE aircraft is
α∘=q+gV·cosθ+zα·α+zδe·δe,q∘=mα·α+mq·q+gV·(mα∘¯·cosθ-c2a·a2·sinθ)+mδe·δe,θ∘=q.
Proof.
System (3.3) is obtained from the system (3.1) for β=p=r=φ=0 and δa=δr=δc=δca=0 and defines the simplified nonlinear model of the motion around the center of gravity of the ADMIRE aircraft in a longitudinal flight with constant forward velocity V.
In system (3.3) g,V,zα,zδe,mα,mq,mα∘¯,c2,a2,a,mδe are considered constants (see Table 1).
The values of the parameters used in the simplified system (3.3).
Parameter
Value
Units
zα
-1.598075
rad/s
zδe
-0.52089
rad/s
mα
1.72514652738
rad/s2
mq
-22.61196
rad/s2
a
-0.485
s-1
mδe
-9.972922
—
a2
11.964
s-2
V
84.5
m/s
g
9.81
m/s2
mα̇¯
-5.26416
rad/s2
c2
-0.029
—
—
—
—
Remark 3.3.
According to the simplified nonlinear model the equilibriums in a longitudinal flight with constant forward velocity are the solutions of the nonlinear system of equations:
q+gV·cosθ+zα·α+zδe·δe=0,mα·α+mq·q+gV·(mα∘¯·cosθ-c2a·a2·sinθ)+mδe·δe=0,q=0.
System (3.4) defines the equilibriums manifold of the longitudinal flight with constant forward velocity V.
Proposition 3.4.
System (3.4) implies that α satisfies
A·α2+B·δe·α+C·δe2+D=0,
where A, B, C, D are given by
A=(mα-mα∘¯·zα)2+c22a2·a22·zα2,B=2·(mα-mα∘¯·zα)·(mδe-mα∘¯·zδe)+2·c22a2·a22·zα·zδe,C=(mδe-mα∘¯·zδe)2+c22a2·a22·zδe2,D=-g2V2·c22a2·a22.
Proof.
Equation (3.5) is obtained replacing q=0 in (3.4)1 and (3.4)2 and eliminating θ between the so-obtained equations.
For the numerical values given in Table 1, (3.5) has real solutions if and only if δe∈[δe,̲δe¯] where δe̲=-0.04678233231992[rad]=-2.681° and δe¯=-δe̲.
The computed α1(δe), θ1(δe), α2(δe), and θ2(δe) solutions are represented on Figures 1 and 2.
The α1(δe) and α2(δe) coordinates of the equilibriums on the manifold ℳV.
The θ1(δe)+2kπ and θ2(δe)+2kπ coordinates of the equilibriums on the manifold ℳV.
The equilibrium manifold ℳV is the union of the following two pieces:
The eigenvalues of the linearized system at an equilibrium are the solutions of the equation:
-λ3+(zα+mq)·λ2+(mα+a23-zα·mq)·λ+mα·gV·sinθ-zα·a23=0,
where a23=-(g/V)·[mα∘¯·sinθ-(c2/a)·a2·cosθ].
Proof.
The proof is given by computation.
Proposition 3.6.
For δe∈(δe̲,δe¯) the equilibriums of 𝒫1 are exponentially stable and those of 𝒫2 are unstable.
Proof.
These results were obtained computing the eigenvalues of the linearized system at the equilibriums of 𝒫1 and at the equilibriums of 𝒫2. More precisely, it was obtained that the eigenvalues are negative real numbers at the equilibriums of 𝒫1, and two of the eigenvalues are negative and the third is positive at the equilibriums of 𝒫2.
Proposition 3.7.
At the equilibriums, which correspond to δe̲ and δe¯, two of the eigenvalues are negative and one eigenvalue is equal to zero. Consequently, δe̲ and δe¯ are nonhyperbolic equilibrium points (turning points).
Proof.
The proof is given by computation.
Remark 3.8.
Transfers between two equilibriums which belong to a conex part of 𝒫1={(α1(δe),0,θ1(δe)):δe∈(δe̲,δe¯)} are possible by small changes of the elevator deflection δe (Figure 3). On the other hand, a small change of the elevator deflection δe transfers an equilibrium which belongs to 𝒫2={(α2(δe),0,θ2(δe)):δe∈(δe̲,δe¯)} into an equilibrium which belongs to 𝒫1 (Figure 4).
A transfer between two equilibriums, which belong to 𝒫1 (simulated in simplified model): δe′=-0.03866[rad]=2.216°→δe′′=0.03026[rad]=1.734°. α1′=0.07866974[rad]=4.509°; q1′=0[rad/s]; θ1′=-0.428832[rad]=24.582°→α1′′=0.065516[rad]=3.755°; q1′′=0[rad/s]; θ1′′=-0.698066[rad]=-40.016°.
A transfer of an equilibrium, which belongs to 𝒫2 into an equilibrium, which belongs to 𝒫1 (simulated in the simplified model): δe′=-0.03866[rad]=2.216°→δe′′=0.03026[rad]=1.734°. α2′=0.064883[rad]=3.719°; q2′=0[rad/s]; θ2′=0.7674624[rad]=43.994°→α1′′=0.065516[rad]=3.755°; q1′′=0[rad/s]; θ1′′=-0.6980667[rad]=-40.016° instead of α2′=0.064883[rad]=3.719°; q2′=0[rad/s]; θ2′=0.7674624[rad]=43.994°→α2′′=0.046845[rad]=2.685°; q2′′=0[rad/s]; θ2′′=1.03669718[rad]=59.428°.
Remark 3.9.
The behavior of the ADMIRE aircraft changes when the maneuver δe′→δe′′ is so that δe′∈(δe̲,δe¯) and δe′′∉(δe̲,δe¯). Computation shows that after such a maneuver α and q oscillate with the same period and θ tends to +∞ or -∞ (Figures 5 and 6).
Oscillations when δe′=-0.04678[rad]=-2.68°→δe′′=-0.05[rad]=-2.866° and the starting point is α1′=0.0869742[rad]=4.985°; q1′=0[rad/sec]; θ1′=0.159329[rad]=9.133°.
Oscillations when δe′=-0.04678[rad]=-2.681°→δe′′=0.048[rad]=2.751° and the starting point is α1′=0.086974288[rad]=4.985°; q1′=0[rad/sec]; θ1′=0.1593297[rad]=9.133°.
Since the nonhyperbolic equilibriums at δe̲ and δe¯ seem to be fundamental for the above behavior, we are going to prove that δe̲ is a saddle point bifurcation for the system (3.3). An analogous proof holds for δe¯.
Proposition 3.10.
The nonhyperbolic equilibrium at δe̲ is a saddle point bifurcation.
Proof.
Let us pose x=(α,q,θ) and write the system (3.3) as ẋ=f(x,δe). Moreover, set Pe̲=(α1(δe̲),0,θ1(δe̲))=(α2(δe̲),0,θ2(δe̲)) the corresponding equilibrium point. As it has been already checked that Pe̲ is a nonhyperbolic equilibrium point, it is sufficient to verify that w·Dδef≠0 and w·Dx2f(v,v)≠0, where v and w are right and left eigenvalues corresponding to the zero eigenvalues, respectively, and the derivatives are computed at Pe̲, δe̲ [11, page 148]. From (3.3) it follows that v=((1/zα)·(g/V)·sinθ,0,1)T and w=(mα/(zα·mq-mα),-zα/(zα·mq-mα),1) are right and left eigenvectors, respectively, corresponding to the zero eigenvalue and
w·Dδef=mα·zδe-zα·mδezα·mq-mα,w·Dx2f(v,v)=1zα·mq-mα·[-gV·cosθ·mα+gV·(mα̇¯·cosθ+c2a·a2·sinθ)·zα].
For the considered numerical data, given in Table 1, we have w·Dδef=-0.489≠0and w·Dx2f(v,v)=2.255≠0 . Therefore, Pe̲ is a saddle-node bifurcation.
Remark 3.11.
Along the path of longitudinal flight equilibriums only saddle-point bifurcation exists. There is no Hopf bifurcation.
Remark 3.12.
In both cases, that is, decrease of the parameter δe under δe̲ or increase of the parameter δe over δe¯, we have a loss of steady state and occurrence of an oscillatory flight. Computation shows that this loss is not catastrophic, because if δe is reset, then a stable equilibrium is recovered, as it is illustrated in Figure 7.
Resetting δe from δe′′=0.048[rad]=2.751°<δe̲ to δe′=-0.04678[rad]=-2.681° after 3000 [s] of oscillations a stable equilibrium is recovered.
Remark 3.13.
It is important to remark that the saddle-point bifurcation phenomenon occurring at δe̲ is of practical interest. Using this phenomenon it is possible to transfer the vehicle from a stable descending longitudinal flight equilibrium with constant forward velocity (θ<0 and θ-α<0) into a state which is appropriate for the touch-down moment.
Definition 3.14.
In a longitudinal flight with constant forward velocity an equilibrium is a descending flight if θ-α<0.
Such an equilibrium is, for instance, α1′=0.078669740237840[rad]=4.509°; q1′=0[rad/s]; θ1′=-0.428832005303479[rad]=-24.582° and corresponds to δe′=-0.03866[rad]=-2.216°. Moreover, it is stable descent flight equilibrium. By the maneuver δe′→δe′′ with δe′′=-0.1[rad]=5.732°<δe̲ the above state is transferred into a state which is appropriate for the touch-down moment θ′′-α′′=0 and θ′′>0 (θ′′ small) (Figure 8).
Figure 8 shows that the maneuver δe′=-0.03866[rad]=-2.216°→δe′′=-0.1[rad]=-5.13° made 20.16 [s] before the touch-down moment, transfers the aircraft in 20.16 seconds from the stable descending flight equilibrium α1′=4.509°; q1′=0; θ1′=-24.5° into the state α1′′=θ1′′=6.94°, which is appropriate for the touch-down moment. This flare maneuver has to be made when the aircraft is at H=840 [m] altitude.
Transfer from α1′=0.0786697[rad]=4.509°; q1′=0[rad/s]; θ1′=-0.428832[rad]=-24.58° into a state which is appropriate for the touch-down moment: θ′′-α′′=0 and θ′′>0 (θ′′ small).
Remark 3.15.
It is important to underline that the transfer period t depends strongly on the initial value δe′ and on the final value δe′′ of the elevator deflection.
If δe′ is fixed and δe′′ decreases, then t decreases and α′′=θ′′ increases.
For instance, when δe′=-0.03866[rad]=-2.216°, we have
δe′′=-0.049[rad]=-2.80°,t=239.72[s],α′′=θ′′=0.089[rad]=5.1°,δe′′=-0.1[rad]=-5.73°,t=20.15[s],α′′=θ′′=0.1211[rad]=6.94°,δe′′=-0.2[rad]=-11.46°,t=8.15[s],α′′=θ′′=0.184[rad]=10.5°.
Moreover, if the time interval between the moment of the elevator deflection change (flare) and the real touch-down moment is larger than the above presented periods of transitions, then the touch-down moment could be catastrophic, due to the oscillations.
Conclusion
In the simplified model it can be shown numerically that there exists a range [δe̲,δe¯] of the values of the elevator deflection δe such that to a value δe from [δe̲,δe¯] a set of equilibriums corresponds. Some of these equilibriums are stable and some of them are unstable. If δe′,δe′′∈(δe̲,δe¯) a change δe′→δe′′ of the elevator deflection δe transfers the vehicle into a stable equilibrium corresponding to δe′′. If δe′∈(δe̲,δe¯) and δe′′<δe̲-k or δe′′>δe¯+k with k>0, then the change δe′→δe′′ of the elevator deflection leads to an oscillatory movement of the vehicle , which is noncatastrophic from mathematical point of view. This is due to a saddle node bifurcation and can be useful for preparing the touch-down. Hopf bifurcations are not present.
4. General Constant forward Velocity Longitudinal Flight Model of the ADMIRE AircraftProposition 4.1.
The general constant forward velocity longitudinal flight model of the ADMIRE aircraft is defined by the system of differential equations:
α∘=q+gV·cos(θ-α)-gV·sin(θ-α)·tanα+[zα·α+zδe·δe]·1cosα,q∘=mα·α+mq·q+gV·(mα∘¯·cosθ-c2a·a2·sinθ)+mδe·δe,θ∘=q.
Proof.
System (4.1) is obtained from the system (2.11), substituting the general aero dynamical forces and moments with those corresponding to the ADMIRE aircraft.
Remark 4.2.
The simplified model (3.3) is obtained from (4.1), approximating cosα with 1 (small angle of attack α≈0). If the simplified system (3.3) describes correctly the real phenomenon, then it has to be topologically equivalent to the system (4.1) [12]. Furthermore, the simplified system (3.3) has to be structurally stable. As far we know, the structural stability of the system (3.3) never was proved. What we intend to prove numerically in this paragraph is that the steady states (stable and unstable) of system (4.1) are exactly the same as those of system (3.3). Moreover, we show that the behavior of the aircraft described by (4.1) is similar to that described by (3.3). This is not a proof of the structural stability of the simplified system, but it is an increase of the believe that the simplified system (3.3) describes correctly the real phenomenon and we are not in the case reported in [13]. In [13] it was shown that the simplified system of differential equations which governs the motion of the automatic—landing flight—experiment (ALFLEX) reentry vehicle, is neither structurally stable, nor topologically equivalent to the general system governing the same motion. In other words, the general and the simplified mathematical models of ALFLEX give different images of the same reality.
Remark 4.3.
According to the general nonlinear model, the equilibriums in a longitudinal flight with constant forward velocity Vare solutions of the nonlinear system of equations:
q+gV·cos(θ-α)-gV·sin(θ-α)·tanα+[zα·α+zδe·δe]·1cosα=0,mα·α+mq·q+gV·(mα∘¯·cosθ-c2a·a2·sinθ)+mδe·δe=0,q=0.
System (4.2) defines the equilibrium manifold of the longitudinal flight with constant forward velocity V in the general model.
System (4.2) implies that α, θ satisfy
gV·cos(θ-α)-gV·sin(θ-α)·tanα+[zα·α+zδe·δe]·1cosα=0,mα·α+gV·(mα∘¯·cosθ-c2a·a2·sinθ)+mδe·δe=0.
Since (4.3)1 can be written in the form
1cosα·[gV·cosθ+zα·α+zδe·δe]=0,
the system (4.3) has the same solutions as the system
gV·cosθ+zα·α+zδe·δe=0,mα·α+gV·(mα∘¯·cosθ-c2a·a2·sinθ)+mδe·δe=0.
It follows in this way.
Proposition 4.4.
In the case of the longitudinal flight with constant forward velocity of the ADMIRE aircraft the equilibrium manifold ℳV in the general model is the same as in the simplified model.
Proposition 4.5.
The eigenvalues of the linearized system at equilibrium are the solutions of the equation:
-λ3+(a11+mq)·λ2+(mα+a23-a11·mq)·λ+mα·a13-a11·a23=0,
where
a11=gV·sin(θ-α)+gV·cos(θ-α)·tanα-gV·sin(θ-α)·1cos2α+zαcosα+[zα·α+zδe·δe]·sinαcos2α,a13=-gV·sin(θ-α)-gV·cos(θ-α)·tanα,a23=gV·[-mα·¯·sinθ-c2a·a2·cosθ]
or
a11=zαcosα+gV·sinθcosα-gV·sin(θ-α)cos2α+[zα·α+zδe·δe]·sinαcos2α,a13=-gV·sinθcosα,a23=-gV·[mα∘¯·sinθ+c2a·a2·cosθ]
and α, q, and θ are the coordinates of the equilibrium.
Proof.
The proof is given by computation.
Remark 4.6.
Computing the roots of (4.6) at the equilibriums of 𝒫1 we find that these are negative real numbers for δe∈(δe̲,δe¯), and for δe=δe̲ or δe=δe¯ one of the roots is equal to zero.
Computing the roots of (4.6) at the equilibriums of 𝒫2 we find that two of them are negative real numbers and one is strictly positive for δe∈(δe̲,δe¯), and for δe=δe̲ or δe=δe¯ one of the roots is equal to zero.
We conclude the following.
Proposition 4.7.
An equilibrium is stable in the general model if and only if it is stable in the simplified model. Moreover, δe̲, δe¯ are saddle-points in both models.
Remark 4.8.
The same transfer maneuver between the same equilibriums which belong to a conex part of 𝒫1={(α1(δe),0,θ1(δe)):δe∈(δe̲,δe¯)}, simulated already in simplified model, is simulated now in the framework of the general model in Figure 9.
A transfer between two equilibriums, which belong to 𝒫1 (simulated in the general model): δe′=-0.03866[rad]=2.216°→δe′′=0.03026[rad]=1.734°. α1′=0.07866974[rad]=4.509°; q1′=0[rad/s]; θ1′=-0.428832[rad]=24.582°→α1′′=0.065516[rad]=3.755°; q1′′=0[rad/s]; θ1′′=-0.698066[rad]=-40.016°.
The same transfer maneuver from an unstable equilibrium which belongs to 𝒫2={(α2(δe),0,θ2(δe)):δe∈(δe̲,δe¯)} into a stable equilibrium which belongs to 𝒫1, simulated already in simplified model, is simulated now in the framework of the general model in Figure 10.
A transfer of an equilibrium, which belongs to 𝒫2, into an equilibrium which belongs to 𝒫1 (simulated in the general model): δe′=-0.03866[rad]=2.216°→δe′′=0.03026[rad]=1.734°. α2′=0.064883[rad]=3.719°; q2′=0[rad/s]; θ2′=0.7674624[rad]=43.994°→α1′′=0.065516[rad]=3.755°; q1′′=0[rad/s]; θ1′′=-0.6980667[rad]=-40.016° instead of α2′=0.064883[rad]=3.719°; q2′=0[rad/s]; θ2′=0.7674624[rad]=43.994°→α2′′=0.046845[rad]=2.685°; q2′′=0[rad/s]; θ2′′=1.03669718[rad]=59.428°.
Remark 4.9.
The above simulations show that from the point of view of this type of transfer, the results obtained in general model are similar to those obtained in the simplified model. More about transfer maneuvers can be found in [14].
Remark 4.10.
The behavior of the ADMIRE aircraft, when the maneuver δe′→δe′′ is so that δe′∈(δe̲,δe¯) and δe′′∉(δe̲,δe¯), simulated in general model, is represented in Figures 11 and 12 and shows that this is similar to that obtained in the simplified model.
Oscillations when δe′=-0.04678[rad]=-2.681°→δe′′=-0.05[rad]=-2.866° and the starting point is α1′=0.0869742[rad]=4.985°; q1′=0[rad/sec]; θ1′=0.159329[rad]=9.133°.
Oscillations when δe′=-0.04678[rad]=-2.681°→δe′′=-0.05[rad]=-2.866° and the starting point is α1′=0.086974288[rad]=4.985°; q1′=0[rad/sec]; θ1′=0.1593297[rad]=9.133°.
Remark 4.11.
Moreover, the simulation in general model shows that the loss of stability, due to maneuver, is noncatastrophic (from mathematical point of view), because if δe is reset, then a stable equilibrium is recovered; see Figure 13.
Resetting δe from δe′′=0.048[rad]=2.751°<δe̲ to δe′=-0.04678[rad]=-2.681° after 3000 [s] of oscillations a stable equilibrium is recovered.
Remark 4.12.
The bifurcation phenomenon at δe̲, present also in the general model, can be used for the transfer of the aircraft from a stable descending longitudinal flight equilibrium into a state which is appropriate for the touch-down moment. The simulation of such a transfer is presented in Figure 14.
Transfer from α1′=0.0786697[rad]=4.509°; q1′=0[rad/s]; θ1′=-0.428832[rad]=-24.58° into a state which is appropriate for the touch-down moment θ′′-α′′=0 and θ′′>0 (θ′′ small).
Remark 4.13.
The flare maneuver has to be made also when the aircraft is at 840 [m] altitude, but the transfer period, computed in the general model can be different, as it is shown in Table 2.
The transfer periods of the transitions corresponding to the maneuver δe'=-0.03866[rad]=-2.216°→δe''.
δe′′[rad]
-0.048
-0.05
-0.06
-0.08
-0.1
-0.5
Transfer period [s]
333
170
70
30
20
4
5. Conclusions
(i) Numerical computation shows that there exists a range [δe̲,δe¯] of the elevator deflection values δe, such that to a value δe from [δe̲,δe¯] the same set of equilibriums corresponds in both models.
(ii) An equilibrium is stable in the simplified model if and only if it is stable also in the general model.
(iii) If δe′,δe′′∈[δe̲,δe¯], a change δe′→δe′′ of the elevator deflection δe transfers the vehicle into a stable equilibrium corresponding to δe′′ in both models.
(iv) If δe′∈[δe̲,δe¯] and δe′′<δe-k or δe′′>δe¯+k with k>0, then the change δe′→δe′′ of the elevator deflection leads to a noncatastrophic oscillatory movement of the vehicle, in both models. This is due to the saddle node bifurcation in both models and can be used for preparing the touch-down in both models.
(v) In both models there is no Hopf bifurcation.
(vi) The similarity of the computed results increases the confidence that the simplified model correctly describes the real longitudinal flight, but it is not a proof of the structural stability of the simplified model. Such a proof would be welcome.
Acknowledgments
The authors thank the anonymous referees for their valuable comments, which led to the improvement of the manuscript. St. Balint and A. M. Balint thank the Romanian National Authority for Research for supporting the research under the Grant ID 354 no.7/2007.
MehraR. K.PrasanthR. K.Bifurcation and limit cycle analysis of nonlinear pilot induced oscillationsProceedings of the AIAA Atmospheric Flight Mechanics Conference and ExhibitAugust 1998Boston, Mass, USAAIAA-1998-424919961–4OhioWright-Patterson AFBWL-TR-96KishB. A.A limited flight test investigation of pilot-induced oscillation due to rate limitingProceedings of the AIAA Guidance, Navigation, and Control ConferenceAugust 1997New Orleans, La, USAAIAA Paper 97-3703KlydeD. H.McRuerD. T.MyersT. T.Pilot-induced oscillation analysis and prediction with actuator rate limiting19972018189EID2-s2.0-003078088510.2514/2.3998MehraR. K.Global stability and control analysis of aircraft at high angles of attack1977ONR-CR215-248IonitaA.HalanayA.Delay induced oscillationsProceedings of the AIAA Atmospheric Flight Mechanics ConferenceAugust 1997New Orleans, La, USAAIAA-97-3502ShammaJ. S.AthansM.Guaranteed properties of gain scheduled control for linear parameter-varying plants1991273559564MR1104514ZBL0754.93022EtkinB.ReidL.1996New York, NY, USAJohn Wiley & SonsCookM.1997New York, NY, USAJohn Wiley & SonsBalintSt.BalintA. M.IonitaA.2009224Journal of Aerospace Engineering. In pressGuckenheimerJ.HolmesP.198342New York, NY, USASpringerxvi+453Applied Mathematical SciencesMR709768KuznetsovY. A.19981122ndNew York, NY, USASpringerxx+591Applied Mathematical SciencesMR1711790KaslikE.ekaslik@yahoo.comBalintS.balint@balint.uvt.roStructural stability of simplified dynamical system governing motion of ALFLEX reentry vehicle2007204215219EID2-s2.0-2804444697010.1061/(ASCE)0893-1321(2007)20:4(215)BalintŞ.BrăescuL.KaslikE.20081Cambridge, UKCambridge Scientificii+206Mathematical Problems in Engineering and Aerospace SciencesMR2494388