# Patent

## Title
Method for continuously and adaptively generating a speed setpoint for an aircraft to observe an RTA

## Abstract

Method for generating a flight speed profile of an aircraft (100) according to an RTA constraint, characterized in that it permanently calculates, at any point of the flight plan, a speed setpoint VΩ(X) of the aircraft determined on the basis of reference speed profiles, comprising at least a minimum speed profile Vmin(X), and a maximum speed profile Vmax(X). The method can also take into account an optimum speed profile Vopt(X) determined, for example, on the basis of a cost or performance index.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
This application claims priority to foreign France patent application No. 0902613, filed on May 29, 2009, the disclosure of which is hereby incorporated by reference in its entirety.
FIELD OF THE INVENTION
The present invention relates to a method for continuously and adaptively generating a speed setpoint for an aircraft to observe an RTA. It applies notably to the field of avionics, and more particularly to a flight management device of an aircraft, designated hereinafter by the acronym FMS, standing for Flight Management System.
BACKGROUND OF THE INVENTION
In the general context of air traffic, it is a fact that the traffic density is increasing year-on-year. At the same time, it is necessary for the environmental impact of the air traffic to be reduced. Finally, it is necessary for the accident rate to be kept at the current level, or even reduced. Respecting these principles and requirements notably requires the air traffic control organizations to exercise increasingly tight control on the flight profiles of the aircraft occupying the air space. At the same time, it is necessary for the operators of these aircraft to implement devices whose function is to satisfy the requests originating from the air traffic control organizations, while minimizing their impact, notably on the piloting procedures, the initial flight plan, and observance of the provisions in terms of aircraft operating cost.
Typically, the air traffic control organizations can impose passages at given points of the air space on the aircraft, at given times, or RTA (Requested Time of Arrival) for a scheduled time constraint or CTA (Controlled Time of Arrival) for a time constraint activated by the air traffic controller. Hereinafter, these two types of constraint will be designated RTA without distinction. To observe an RTA, the pilot of an aircraft must control the speed profile along the flight plan. To this end, there are devices known from the prior art that assist or even supplant the pilot; such devices are notably implemented in the FMSs of aircraft that are equipped therewith.
Typically, such devices determine a performance index, for example a cost index usually designated by the acronym CI, which represents the ratio between the cost in terms of flight time and the cost in terms of fuel. Other indices can be employed; these indices are determined iteratively, and define a proportionality coefficient between different speed profiles. For a given index, the proportionality coefficient is applied uniformly over the entire profile, and makes it possible to obtain, at any point of the flight plan, the speed setpoint to be applied. In a simplified manner, one process then consists in arbitrarily setting a value of the index, calculating the corresponding speed profile, and estimating the time of passage, or ETA (Estimated Time of Arrival), at the point for which the RTA is applied. The index is then corrected, and the process reiterates the calculation of the ETA according to the new speed profile resulting therefrom; this continues until the ETA is sufficiently close to the RTA according to a predetermined criterion corresponding, for example, to a desired maximum time difference. Such a method has a certain number of drawbacks:

    
    
        determining the speed profile with which to observe the RTA requires a number of iterations that may be high, the number of iterations being relatively unpredictable and variable from one calculation to another;
        in some cases, the convergence of the calculation of the index is difficult, because of complex behaviours of the function linking the performance index to the time of passage at the point at which the RTA is applied;
        if there is a drift during a flight relative to the calculated profile, for example because of wind gusts or even unexpected aircraft performance levels, it is necessary for the calculation of the speed profile to be repeated in full. It is also necessary to ensure a relatively short periodicity in the performance of the calculation, in order to keep to precise predictions;
        on approaching the point at which the RTA is applied, it becomes difficult to compensate sufficiently responsively for any drifts.
    
    


A first method, described in the patent U.S. Pat. No. 6,507,782, allows for an adaptation of the speed profile according to a local sensitivity of the time of passage to the speed variation. However, this sensitivity is globally consolidated to compensate for the total error on the time of passage. This type of method does not call into question the principle whereby the speed profile must be determined by successive iterations, and each iteration can be likened to an open-loop calculation of the impact of the speed profile adopted on the time of passage resulting therefrom. Furthermore, drifts exhibited during the flight are compensated according to the same method, by iterations over the entire speed profile ahead of the aeroplane.
A second method, described in the patent U.S. Pat. No. 5,121,325, consists in varying the precision margins on the observance of the time of passage, as a function of time, so as not to unnecessarily constrain the speed profile when the aircraft is far from the constraint associated with the RTA, and on the other hand to increase the precision on approaching this point. This method makes it possible to optimize the profile, and improve confidence in the observance of the time constraint imposed by an RTA. On the other hand, this method also relies on an open-loop speed prediction, the slaving to the time of passage being based only on a calculation of predictions relating to the flight plan as a whole.
SUMMARY OF THE INVENTION
One aim of the present invention is to mitigate at least the abovementioned drawbacks, by proposing a method for generating a speed setpoint that is permanently slaved to the time of passage resulting therefrom at the point at which the RTA is applied, the method requiring calculations, the periodicity of which lightens the calculation workload of the device implementing it, for example an FMS.
Another advantage of the invention is that it allows for the application of various strategies for generating the speed setpoint, to ensure sufficient margins for satisfying the time constraint, in particular in the case of unstable or imprecise wind, and optimize the speed profile according to the different flight phases and/or according to the distance relative to the point at which the RTA constraint is applied.
To this end, the subject of the invention is a method for continuously and adaptively generating a speed setpoint for the observance of an RTA, characterized in that it comprises the following steps:

    
    
        for a determined current position of the aircraft along a flight plan, determination of at least two reference speed profiles according to the curvilinear abscissa X of the aircraft along the flight plan, comprising a minimum speed profile Vmin(X) and a maximum speed profile Vmax(X),
        storage of a plurality of pairs (TVmin(Xi),Xi) and (TVmax(Xi),Xi) of values of times of passage of the aircraft at curvilinear abscissa points Xi subsequent to the current position as far as the abscissa of the point at which the RTA is applied for the reference speed profiles,
        application at any point of curvilinear abscissa X of the flight plan of a speed setpoint VΩ(X) of the aircraft determined on the basis of the speed profiles Vmin(X) and Vmax(X) weighted by proportionality ratios between T(X), TVmin(X) and TVmax(X);

T(X) designating the time of passage of the aircraft at the point of the flight plan of curvilinear abscissa X, TVmax(X) designating the time of passage at the curvilinear abscissa point X, for a flight according to the minimum speed profile Vmin(X) reaching the point at which the RTA is applied by observing the RTA constraint, and TVmin(X) designating the time of passage at the curvilinear abscissa point X, for a flight according to the maximum speed profile Vmax(X) reaching the point at which the RTA is applied by observing the RTA constraint, the values of TVmin(X) and TVmax(X) being determined by a linear interpolation based on the corresponding stored pairs of values for the closest points either side of the curvilinear abscissa point X.

    
    


In one embodiment of the invention, the method for generating a flight speed profile of an aircraft can be characterized in that, for the current position of the aircraft along the flight plan, an additional optimum speed reference speed profile Vopt(X) is determined according to at least one defined flight criterion, the speed profile of the aircraft being constrained to tend towards the optimum speed profile Vopt(X) by the addition to the RTA of a tolerance ΔT the sign of which is opposite to the sign of the difference between the time of passage of the aircraft at the current position and the corresponding time of passage TVopt(X), until the absolute value of the difference between T(X) and TVopt(X) is less than a maximum time difference δt, the RTA then being restored to its initial value.
In one embodiment of the invention, said at least one defined flight criterion can comprise a cost index.
In one embodiment of the invention, said at least one defined flight criterion can comprise a performance index.
In one embodiment of the invention, the method for generating a flight speed profile of an aircraft can be characterized in that the speed setpoint VΩ(X) is determined by the relation:




  
    
      V
      Ω
    
    ⁡
    
      (
      X
      )
    
  
  =
  
    
      
        
          (
          
            
              
                T
                
                  V
                  ⁢
                  
                    
                  
                  ⁢
                  max
                
              
              ⁡
              
                (
                X
                )
              
            
            -
            
              
                T
                
                  V
                  ⁢
                  
                    
                  
                  ⁢
                  min
                
              
              ⁡
              
                (
                X
                )
              
            
          
          )
        
        ⁢
        
          
            V
            min
          
          ⁡
          
            (
            X
            )
          
        
        ⁢
        
          
            V
            max
          
          ⁡
          
            (
            X
            )
          
        
      
      
        
          
            (
            
              
                
                  T
                  
                    V
                    ⁢
                    
                      
                    
                    ⁢
                    max
                  
                
                ⁡
                
                  (
                  X
                  )
                
              
              -
              
                T
                ⁡
                
                  (
                  X
                  )
                
              
            
            )
          
          ⁢
          
            
              V
              max
            
            ⁡
            
              (
              X
              )
            
          
        
        -
        
          
            (
            
              
                
                  T
                  
                    V
                    ⁢
                    
                      
                    
                    ⁢
                    min
                  
                
                ⁡
                
                  (
                  X
                  )
                
              
              -
              
                T
                ⁡
                
                  (
                  X
                  )
                
              
            
            )
          
          ⁢
          
            
              V
              min
            
            ⁡
            
              (
              X
              )
            
          
        
      
    
    .
  




In one embodiment of the invention, the method for generating a flight speed profile of an aircraft can be characterized in that the speed setpoint VΩ(X) is determined by the relation:




  
    
      
        V
        Ω
      
      ⁡
      
        (
        X
        )
      
    
    =
    
      
        
          (
          
            
              
                T
                
                  V
                  ⁢
                  
                    
                  
                  ⁢
                  max
                
              
              ⁡
              
                (
                X
                )
              
            
            -
            
              
                T
                min
              
              ⁡
              
                (
                X
                )
              
            
          
          )
        
        ⁢
        
          
            V
            
              V
              ⁢
              min
            
          
          ⁡
          
            (
            X
            )
          
        
        ⁢
        
          
            V
            max
          
          ⁡
          
            (
            X
            )
          
        
      
      
        
          
            
              
                
                  (
                  
                    
                      
                        T
                        
                          V
                          ⁢
                          
                            
                          
                          ⁢
                          max
                        
                      
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                    -
                    
                      
                        T
                        corrected
                      
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                  
                  )
                
                ⁢
                
                  
                    V
                    max
                  
                  ⁡
                  
                    (
                    X
                    )
                  
                
              
              -
            
          
        
        
          
            
              
                (
                
                  
                    
                      T
                      
                        V
                        ⁢
                        
                          
                        
                        ⁢
                        min
                      
                    
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                  -
                  
                    
                      T
                      corrected
                    
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                
                )
              
              ⁢
              
                
                  V
                  min
                
                ⁡
                
                  (
                  X
                  )
                
              
            
          
        
      
    
  
  ,




Tcorrected(X) being equal to T(X)−δt if T(X) is less than TVopt(X)−δt, to T(X)+δt if T(X) is greater than TVopt(X)+δt, or to T(X) otherwise.

In one embodiment of the invention, the method for generating a flight speed profile of an aircraft can be characterized in that the speed setpoint VΩ(X) is corrected by a value δV(X) equal to




  
    
      V
      opt
    
    ⁡
    
      (
      X
      )
    
  
  ⁢
  
    
      
        
          T
          
            V
            ⁢
            
              
            
            ⁢
            opt
          
        
        ⁡
        
          (
          X
          )
        
      
      -
      
        (
        
          
            T
            ⁡
            
              (
              X
              )
            
          
          -
          
            Δ
            ⁢
            
              
            
            ⁢
            T
          
        
        )
      
    
    
      
        T
        RTA
      
      -
      
        (
        
          
            T
            ⁡
            
              (
              X
              )
            
          
          -
          
            Δ
            ⁢
            
              
            
            ⁢
            T
          
        
        )
      
    
  




    
    
        if T(X) is less than TVopt(X)−ΔT, or to
    
    






  
    
      V
      opt
    
    ⁡
    
      (
      X
      )
    
  
  ⁢
  
    
      
        
          T
          
            V
            ⁢
            
              
            
            ⁢
            opt
          
        
        ⁡
        
          (
          X
          )
        
      
      -
      
        (
        
          
            T
            ⁡
            
              (
              X
              )
            
          
          +
          
            Δ
            ⁢
            
              
            
            ⁢
            T
          
        
        )
      
    
    
      
        T
        RTA
      
      -
      
        (
        
          
            T
            ⁡
            
              (
              X
              )
            
          
          +
          
            Δ
            ⁢
            
              
            
            ⁢
            T
          
        
        )
      
    
  




    
    
        if T(X) is greater than TVopt(X)+ΔT.
    
    


In one embodiment of the invention, the method for generating a flight speed profile of an aircraft can be characterized in that it generates a plurality of speed profiles so as to observe a plurality of successive time of passage constraints along the flight plan
In one embodiment of the invention, the method for generating a flight speed profile of an aircraft can be characterized in that a plurality of different speed profiles are generated for different flight phases for which the optimum speed profiles are determined on the basis of different flight criteria.
Another subject of the present invention is a flight management system of an aircraft, characterized in that it implements a method according to any one of the preceding claims.
In one embodiment of the invention, the flight management system of an aircraft can be characterized in that it calculates, at any instant of the flight of the aircraft, a speed setpoint to be applied to it.
In one embodiment of the invention, the flight management system of an aircraft can be characterized in that the speed setpoint to be applied is displayed on a display interface.
In one embodiment of the invention, the flight management system of an aircraft can be characterized in that the speed setpoint to be applied is transmitted to an automatic pilot system.



BRIEF DESCRIPTION OF THE DRAWINGS
Other features and advantages of the invention will become apparent upon reading the description, given by way of example, and in light of the appended drawings which represent:
FIG. 1, examples of curves representing the times of passage of an aircraft at curvilinear abscissae along its flight plan, for flights following different speed profiles;
FIG. 2, examples of curves representing the times of passage of an aircraft at curvilinear abscissae along its flight plan, for flights following different speed profiles, transposed so as to culminate in a given RTA at the curvilinear abscissa corresponding to the point of constraint;
FIG. 3, examples of curves representing the times of passage of an aircraft, in a case where the aircraft is converging towards an optimum profile according to one embodiment of the invention;
FIG. 4, examples of time-of-passage curves of an aircraft illustrating a speed correction calculation principle according to the invention;
FIGS. 5a and 5b, respectively a flow diagram illustrating the various strategies that can be applied in the climb phase according to different configurations, and the corresponding time-of-passage curves, according to an exemplary embodiment of the invention;
FIGS. 6a and 6b, respectively a flow diagram illustrating the various strategies that can be applied in the cruise phase according to different configurations, and the corresponding time-of-passage curves, according to an exemplary embodiment of the invention; and
FIG. 7, a flow diagram illustrating the various strategies that can be applied in the descent phase according to different configurations, according to an exemplary embodiment of the invention.



DETAILED DESCRIPTION
FIG. 1 shows examples of curves representing the times of passage of an aircraft at curvilinear abscissae along its flight plan, for flights following different speed profiles.
Hereinafter, the letter X designates the curvilinear abscissa of an aircraft 100, along a flight plan, that is to say, a distance between a point of the flight plan and a reference point of the flight plan. A prediction of the time of passage T(X) of the aircraft 100 at the curvilinear abscissa X is made at an abscissa point Xpred.
A first curve 1 represents the times of passage ETA_Vmin(X) of the aircraft 100 according to the curvilinear abscissa X, assuming that the latter is moving according to a minimum speed profile Vmin(X) corresponding to the lower limit of the flight envelope along the profile, possibly augmented by an operational margin in order, for example, to limit acceleration inertia at low speeds or even the temporary risks of excessively low speed generated, for example, by wind shifts. This minimum speed profile is known for the aircraft 100, and depends notably on its performance characteristics, according to the different flight configurations and in the different flight phases according to the flight plan, and on possible additional limitations specified by the pilot or the airline responsible for operating the aircraft.
A second curve 2 represents the times of passage ETA_Vmax(X) of the aircraft 100 according to the curvilinear abscissa X, assuming that the latter is moving according to a maximum speed profile Vmax(X) corresponding to the upper limit of the flight envelope along the profile, possibly reduced by an operational margin in order, for example, to limit the temporary risks of overspeed generated, for example, by wind shifts. Similarly, this maximum speed profile is known for the aircraft 100, and depends notably on its performance characteristics, according to the different flight configurations and in the different flight phases according to the flight plan, and on possible additional limitations specified by the pilot or the airline responsible for operating the aircraft.
A third curve 3 represents the times of passage ETA_Vopt(X) of the aircraft 100 according to the curvilinear abscissa X, assuming that the latter is moving according to an optimum speed profile Vopt(X). This optimum speed profile corresponds to optimum flight conditions that do not take into account the RTA constraint, according to determined criteria, for example flight comfort, fuel consumption, or even according to a performance criterion or an index representative of the time/consumed fuel trade-off, etc.
In parallel with the three curves 1, 2 and 3 representing the times of passage of the aircraft 100 according to the curvilinear abscissa X along the flight plan, examples of curves 11, 12 and 13 representing the speed profiles are given in the figure. The three time-of-passage curves tend to be increasing, and the gradient of a time-of-passage curve becomes all the shallower as the speed of the aircraft 100 increases.
It should be noted that, hereinafter, it is considered as a first approximation that the lateral trajectories of the aircraft 100 are identical, regardless of the speed profile concerned. The reference trajectory of the aircraft 100 can, for example, be chosen as the trajectory corresponding to the optimum speed profile Vopt(X). This approximation does not of course correspond to reality, in which the flight mechanics impose different trajectories, in the horizontal plane, depending on the flight speeds of the aircraft 100, primarily during turn phases or transitions between different legs of the flight plan. It is also possible to consider the trajectory to correspond to a known speed profile that makes it possible to observe the RTA, if such a speed profile is available. It is even possible to calculate the trajectory transitions on the basis of the speed of the profile, augmented by a tolerance margin that makes it possible, from one refresh to another in the calculations, to compensate for any drifts. In practice, the duly defined trajectory cannot in reality be implemented at maximum speed for example, but the maximum speed profile is intended only to provide a time-of-passage reference on the curvilinear abscissa of a common trajectory.
It should also be observed that, regardless of the reference speed profile concerned, that is to say Vmin(X), Vmax(X) or Vopt(X), the ground speed of the aircraft 100 depends only on the curvilinear abscissa X, and not on time. This approximation is based on the assumption that the properties of a stationary system are applicable. These stationarity properties are linked to the fact that the atmospheric model used is stationary. However, it is possible to consider as an approximation that the system is virtually stationary, if variations in meteorological data are introduced, but remain slow, which is generally the practical case, given the time-of-passage differences between the profile Vmin(X) and the profile Vmax(X).
The times of passage T(X) predicted at the curvilinear abscissa Xpred are stored, to be used as and when predictions are subsequently refreshed.
According to the method of the present invention, the FMS notes the time of passage of the aircraft 100 at particular points of the flight plan, for example at waypoints WPT.
FIG. 2 shows examples of curves representing the times of passage of an aircraft at curvilinear abscissae along its flight plan, for flights following different speed profiles, transposed so as to culminate in a given RTA at the curvilinear abscissa XRTA corresponding to the point of constraint. FIG. 2 is based on the same speed profiles as FIG. 1, and illustrates the principles described hereinbelow.
The time of passage at a current abscissa point Xcurr is denoted T(Xcurr). At the same time, the FMS evaluates TVmax(Xcurr), which corresponds to the time of passage of the aircraft 100 at the curvilinear abscissa point Xcurr, required for the latter to fly according to the maximum speed profile Vmax and for it to reach the curvilinear abscissa point XRTA at which the RTA is applied, at the required RTA time, hereinafter denoted TRTA. Similarly, the FMS evaluates TVmin(Xcurr), which corresponds to the time of passage of the aircraft 100 at the curvilinear abscissa point Xcurr, required for the latter to fly according to the minimum speed profile Vmin and for it to reach the curvilinear abscissa point XRTA at the required time TRTA. Finally, the FMS evaluates TVopt(Xcurr), which corresponds to the time of passage of the aircraft 100 at the curvilinear abscissa point Xcurr, required for the latter to fly according to the optimum speed profile Vopt(x) and for it to reach the abscissa point XRTA at the required time TRTA. In FIG. 2, these principles are illustrated by a translation of the curves of FIG. 1 representing the times of passage according to the different speed profiles, so that they all meet at the abscissa point XRTA and ordinate point TRTA.
It should be noted that the required times of passage T(X) according to the different speed profiles Vmin(X), Vmax(X) and Vopt(X) are calculated by the FMS, for example with periodic refreshes, or done at given points of the flight plan. Such a point is represented in FIGS. 1 and 2, and its curvilinear abscissa is denoted Xpred. In an exemplary embodiment, the values TVmin(X), TVmax(X) and TVopt(X) can be recalculated based on the profiles ETAVmin(X), ETAVmax(X) and ETAVopt(X) stored in the memory of the FMS or of a device linked to the FMS, for a set of discrete points of the trajectory, corresponding, for example, to the waypoints WPT of the flight plan. The memory of the device then contains a plurality of pairs of values (ETAVmin(Xi),Xi), (ETAVmax(Xi),Xi) and (ETAVopt(Xi),Xi). The values ETAVmin(X), ETAVmax(X) and ETAVopt(X) can then be calculated by interpolation between the stored values at the points preceding and following the current position. The robustness of the method according to the present invention makes it possible to employ a linear interpolation technique with sufficient accuracy.
The values TVmin(X), TVmax(X) and TVopt(X) are then established according to the following formulae:

    
    
        Time of passage required to observe the RTA by flying according to the speed profile Vmin:

TVmin(x)=ETAVmin(x)+TRTA−ETAVmin(XRTA)

        Time of passage required to observe the RTA by flying according to the speed profile Vopt:

TVopt(x)=ETAVopt(x)+TRTA−ETAVopt(XRTA)

        Time of passage required to observe the RTA by flying according to the speed profile Vmax:

TVmax(x)=ETAVmax(x)+TRTA−ETAVmax(XRTA).

    
    


By definition, the ground speed at a curvilinear abscissa point X can be expressed as:




  
    
      
        
          V
          ⁡
          
            (
            X
            )
          
        
        =
        
          
            
              d
              X
            
            
              d
              t
            
          
          .
        
      
    
    
      
        (
        1
        )
      
    
  




It is thus also possible to write:




  dt
  =
  
    
      dX
      
        V
        ⁡
        
          (
          X
          )
        
      
    
    .
  




By incorporating dt along the trajectory, from the current curvilinear abscissa Xcurr to the curvilinear abscissa XRTA, the following relation is obtained:




  
    
      
        
          T
          ⁡
          
            (
            
              X
              RTA
            
            )
          
        
        =
        
          
            T
            ⁡
            
              (
              
                X
                curr
              
              )
            
          
          +
          
            
              ∫
              
                X
                curr
              
              
                X
                RTA
              
            
            ⁢
            
              
            
            ⁢
            
              
                d
                X
              
              /
              
                
                  V
                  ⁡
                  
                    (
                    X
                    )
                  
                
                .
              
            
          
        
      
    
    
      
        (
        2
        )
      
    
  




Hereinafter, in the interests of simplicity, we will use φN to denote the flight time between two points at a speed VN(X), VN(X) representing Vmin(X), Vopt(X) or Vmax(X); thus:




  
    
      
        
          
            
              φ
              VN
            
            ⁡
            
              (
              
                
                  X
                  curr
                
                ,
                
                  X
                  RTA
                
              
              )
            
          
          =
          
            
              ∫
              
                X
                curr
              
              
                X
                RTA
              
            
            ⁢
            
              
            
            ⁢
            
              
                d
                X
              
              /
              
                
                  V
                  N
                
                ⁡
                
                  (
                  X
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        3
        )
      
    
  
  
    
      and
    
    
      
        
      
    
  
  
    
      
        
          
            T
            VN
          
          ⁡
          
            (
            
              X
              RTA
            
            )
          
        
        =
        
          
            
              T
              VN
            
            ⁡
            
              (
              
                X
                curr
              
              )
            
          
          +
          
            
              
                φ
                VN
              
              ⁡
              
                (
                
                  
                    X
                    curr
                  
                  ,
                  
                    X
                    RTA
                  
                
                )
              
            
            .
          
        
      
    
    
      
        (
        4
        )
      
    
  




According to the method of the invention, the FMS seeks to determine a speed profile VΩ(x) in which the times of passage TΩ(x) along the trajectory satisfy:




  
    
      
        {
        
          
            
              
                
                  
                    T
                    Ω
                  
                  ⁡
                  
                    (
                    
                      X
                      curr
                    
                    )
                  
                
                =
                
                  T
                  curr
                
              
            
          
          
            
              
                
                  
                    
                      T
                      Ω
                    
                    ⁡
                    
                      (
                      
                        X
                        RTA
                      
                      )
                    
                  
                  =
                  
                    T
                    RTA
                  
                
                ;
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  




    
    
        or even:

TΩ(XRTA)=TΩ(Xcurr)+φΩ(Xcurr, XRTA)   (6).

    
    


Furthermore, according to the definition of the time profiles required at Vmin(X) and Vmax(X) given hereinabove, it is possible to write the following relation:




  
    
      
        {
        
          
            
              
                
                  
                    T
                    
                      V
                      ⁢
                      
                        
                      
                      ⁢
                      min
                    
                  
                  ⁡
                  
                    (
                    
                      X
                      RTA
                    
                    )
                  
                
                =
                
                  T
                  RTA
                
              
            
          
          
            
              
                
                  
                    T
                    
                      V
                      ⁢
                      
                        
                      
                      ⁢
                      max
                    
                  
                  ⁡
                  
                    (
                    
                      X
                      RTA
                    
                    )
                  
                
                =
                
                  
                    T
                    RTA
                  
                  .
                
              
            
          
        
      
    
    
      
        (
        7
        )
      
    
  




It is then possible to combine the preceding relations, to obtain the following relation:




  
    
      
        
          
            
              
                
                  φ
                  Ω
                
                ⁡
                
                  (
                  
                    
                      X
                      curr
                    
                    ,
                    
                      X
                      RTA
                    
                  
                  )
                
              
              -
              
                
                  φ
                  min
                
                ⁡
                
                  (
                  
                    
                      X
                      curr
                    
                    ,
                    
                      X
                      RTA
                    
                  
                  )
                
              
            
            
              
                
                  T
                  
                    V
                    ⁢
                    
                      
                    
                    ⁢
                    min
                  
                
                ⁡
                
                  (
                  
                    X
                    curr
                  
                  )
                
              
              -
              
                T
                curr
              
            
          
          =
          
            
              
                
                  φ
                  Ω
                
                ⁡
                
                  (
                  
                    
                      X
                      curr
                    
                    ,
                    
                      X
                      RTA
                    
                  
                  )
                
              
              -
              
                
                  φ
                  max
                
                ⁡
                
                  (
                  
                    
                      X
                      curr
                    
                    ,
                    
                      X
                      RTA
                    
                  
                  )
                
              
            
            
              
                
                  T
                  
                    V
                    ⁢
                    
                      
                    
                    ⁢
                    max
                  
                
                ⁡
                
                  (
                  
                    X
                    curr
                  
                  )
                
              
              -
              
                T
                curr
              
            
          
        
        ,
      
    
    
      
        (
        8
        )
      
    
  




    
    
        which can also be written as follows:
    
    






  
    
      
        
          
            φ
            Ω
          
          ⁡
          
            (
            
              
                X
                curr
              
              ,
              
                X
                RTA
              
            
            )
          
        
        =
        
          
            
              
                
                  
                    
                      
                        (
                        
                          
                            
                              T
                              
                                
                                  _
                                  ⁢
                                  V
                                
                                ⁢
                                
                                  
                                
                                ⁢
                                max
                              
                            
                            ⁡
                            
                              (
                              
                                X
                                curr
                              
                              )
                            
                          
                          -
                          
                            T
                            curr
                          
                        
                        )
                      
                      ⁢
                      
                        
                          φ
                          min
                        
                        ⁡
                        
                          (
                          
                            
                              X
                              curr
                            
                            ,
                            
                              X
                              RTA
                            
                          
                          )
                        
                      
                    
                    -
                  
                
              
              
                
                  
                    
                      (
                      
                        
                          
                            T
                            
                              V
                              ⁢
                              
                                
                              
                              ⁢
                              min
                            
                          
                          ⁡
                          
                            (
                            
                              X
                              curr
                            
                            )
                          
                        
                        -
                        
                          T
                          curr
                        
                      
                      )
                    
                    ⁢
                    
                      
                        φ
                        max
                      
                      ⁡
                      
                        (
                        
                          
                            X
                            curr
                          
                          ,
                          
                            X
                            RTA
                          
                        
                        )
                      
                    
                  
                
              
            
            
              
                
                  T
                  
                    V
                    ⁢
                    
                      
                    
                    ⁢
                    max
                  
                
                ⁡
                
                  (
                  
                    X
                    curr
                  
                  )
                
              
              -
              
                
                  T
                  
                    V
                    ⁢
                    
                      
                    
                    ⁢
                    min
                  
                
                ⁡
                
                  (
                  
                    X
                    curr
                  
                  )
                
              
            
          
          .
        
      
    
    
      
        (
        9
        )
      
    
  




The relation (9) can be satisfied by a large number of speed profiles VΩ. It is possible, in one embodiment of the invention, to choose a speed profile that satisfies the relation (9) at any point of the trajectory. This is equivalent to considering that the FMS determines a speed profile such that the proportion ratios between the flight times to RTA, or: TRTA−TVmin(X), TRTA−TVmax(X) and TRTA−TVopt(X), remain constant along the entire trajectory.
Thus, according to this embodiment, the speed profile can be determined by satisfying the following relation:




  
    
      
        
          1
          
            
              V
              Ω
            
            ⁡
            
              (
              X
              )
            
          
        
        =
        
          
            
              
                
                  (
                  
                    
                      
                        T
                        
                          V
                          ⁢
                          
                            
                          
                          ⁢
                          max
                        
                      
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                    -
                    
                      T
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                  
                  )
                
                ⁢
                
                  1
                  
                    
                      V
                      min
                    
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                
              
              -
              
                
                  (
                  
                    
                      
                        T
                        
                          V
                          ⁢
                          
                            
                          
                          ⁢
                          min
                        
                      
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                    -
                    
                      T
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                  
                  )
                
                ⁢
                
                  1
                  
                    
                      V
                      max
                    
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                
              
            
            
              
                
                  T
                  
                    V
                    ⁢
                    
                      
                    
                    ⁢
                    max
                  
                
                ⁡
                
                  (
                  X
                  )
                
              
              -
              
                
                  T
                  
                    V
                    ⁢
                    
                      
                    
                    ⁢
                    min
                  
                
                ⁡
                
                  (
                  X
                  )
                
              
            
          
          .
        
      
    
    
      
        (
        10
        )
      
    
  




The speed setpoint to be applied at any point of the trajectory is then written:




  
    
      
        
          
            V
            Ω
          
          ⁡
          
            (
            X
            )
          
        
        =
        
          
            
              
                (
                
                  
                    
                      T
                      
                        V
                        ⁢
                        
                          
                        
                        ⁢
                        max
                      
                    
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                  -
                  
                    
                      T
                      
                        V
                        ⁢
                        
                          
                        
                        ⁢
                        min
                      
                    
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                
                )
              
              ⁢
              
                
                  V
                  min
                
                ⁡
                
                  (
                  X
                  )
                
              
              ⁢
              
                
                  V
                  max
                
                ⁡
                
                  (
                  X
                  )
                
              
            
            
              
                
                  (
                  
                    
                      
                        T
                        
                          V
                          ⁢
                          
                            
                          
                          ⁢
                          max
                        
                      
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                    -
                    
                      T
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                  
                  )
                
                ⁢
                
                  
                    V
                    max
                  
                  ⁡
                  
                    (
                    X
                    )
                  
                
              
              -
              
                
                  (
                  
                    
                      
                        T
                        
                          V
                          ⁢
                          
                            
                          
                          ⁢
                          min
                        
                      
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                    -
                    
                      T
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                  
                  )
                
                ⁢
                
                  
                    V
                    min
                  
                  ⁡
                  
                    (
                    X
                    )
                  
                
              
            
          
          .
        
      
    
    
      
        (
        11
        )
      
    
  




It is therefore possible, knowing the required time-of-passage and ground speed at maximum and minimum speed profiles, to deduce therefrom a speed setpoint which takes account of the current ETA. This setpoint can be used both in the calculation of the predictions along the flight plan, and in the guidance setpoints to ensure that the RTA is observed, in particular on approaching the latter.
The trend of the predicted ETA along the trajectory of the aircraft 100 can then be illustrated by the curve 200 represented in FIG. 2, starting from the aircraft 100 at the abscissa Xcurr and ending at the abscissa point XRTA and ordinate point TRTA.
FIG. 3 shows examples of curves representing the times of passage of an aircraft, in a case where the aircraft is converging towards an optimum profile, according to one embodiment of the invention.
It may in fact be advantageous to force the convergence of the flight time profile of the aircraft 100 towards the time profile corresponding to the optimum speed profile Vopt(X). The way this convergence is determined is based on the choice of tolerance margins that are to be maintained to ensure that the RTA is observed.
Furthermore, this makes it possible to guarantee the observance of the RTA, while approaching as far as possible the optimum speed condition for which the initial trajectory was calculated.
To this end, the formula (11) is used, but it is biased to converge according to the case towards RTA, (RTA+ΔT), or (RTA−ΔT), in which ΔT designates a chosen tolerance, depending on whether the time of passage of the aircraft is close to, greater than or less than the time of passage that the aircraft would present by following the optimum speed profile Vopt(X). In practice, the bias is applied in an opposite manner to the value of T(X).
δt designates the desired maximum time deviation relative to the optimum speed profile.
The determination of the setpoint speed is made on the basis of a corrected time of passage Tcorrected(X), as follows:
If T(X)<TVopt(X)−δt, then:

Tcorrected(X)=T(X)−ΔT   (12):

    
    
        Else, if T(X)>TVopt(X)+δt, then:

Tcorrected(X)=T(X)+ΔT   (13);

        Else:

Tcorrected(X)=T(X)   (14).

    
    


The speed setpoint taking into account the corrected time of passage can then be written according to the following relation:




  
    
      
        
          
            V
            Ω
          
          ⁡
          
            (
            X
            )
          
        
        =
        
          
            
              
                (
                
                  
                    
                      T
                      
                        V
                        ⁢
                        
                          
                        
                        ⁢
                        max
                      
                    
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                  -
                  
                    
                      T
                      
                        V
                        ⁢
                        
                          
                        
                        ⁢
                        min
                      
                    
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                
                )
              
              ⁢
              
                
                  V
                  min
                
                ⁡
                
                  (
                  X
                  )
                
              
              ⁢
              
                
                  V
                  max
                
                ⁡
                
                  (
                  X
                  )
                
              
            
            
              
                
                  
                    
                      
                        (
                        
                          
                            
                              T
                              
                                V
                                ⁢
                                
                                  
                                
                                ⁢
                                max
                              
                            
                            ⁡
                            
                              (
                              X
                              )
                            
                          
                          -
                          
                            
                              T
                              corrected
                            
                            ⁡
                            
                              (
                              X
                              )
                            
                          
                        
                        )
                      
                      ⁢
                      
                        
                          V
                          max
                        
                        ⁡
                        
                          (
                          X
                          )
                        
                      
                    
                    -
                  
                
              
              
                
                  
                    
                      (
                      
                        
                          
                            T
                            
                              V
                              ⁢
                              
                                
                              
                              ⁢
                              min
                            
                          
                          ⁡
                          
                            (
                            X
                            )
                          
                        
                        -
                        
                          
                            T
                            corrected
                          
                          ⁡
                          
                            (
                            X
                            )
                          
                        
                      
                      )
                    
                    ⁢
                    
                      
                        V
                        min
                      
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                  
                
              
            
          
          .
        
      
    
    
      
        (
        15
        )
      
    
  




In this way, the time profile of the flight of the aircraft 100 is forced to converge towards the profile corresponding to the optimum speed profile Vopt(X), until in the vicinity δt of the latter, then slaving is applied to the observance of the RTA according to the optimum speed profile Vopt(X).
In the example illustrated by FIG. 3, the aircraft 100, in its current curvilinear abscissa position Xcurr, is delayed relative to the assumed time of passage that it would have by following an optimum speed profile Vopt(X) denoted by the reference symbol 3. Thus, the time at the RTA point is reduced by a tolerance ΔT, and the method of the invention is implemented by targeting an RTA that is early relative to the current RTA. This RTA point corrected by the tolerance ΔT is targeted until the difference between the time of passage of the aircraft 100 and the time of passage that it would have by following the optimum speed profile 3 is below a chosen maximum time deviation δt. Then, it is the current RTA that is targeted.
According to one embodiment of the invention, it may be useful, in order to apply a strategy for optimizing the speed profile, to seek to achieve a given time of passage at an intermediate position along the flight plan. For example, it may be desirable to get to the top of climb, usually designated T/C, at a time TVmax (time of passage required to achieve the RTA at Vmin) in such a way as to accelerate as little as possible when climbing.
To get to the target position whose curvilinear abscissa is denoted Xtgt at the time Ttgt, the preceding relations are applied in the same way:

TVmin(x)=ETA—Vmin(x)+Ttgt−ETA—Vmin(Xtgt)

TVmax(x)=ETA—Vmax(x)+Ttgt−ETA—Vmax(Xtgt)

If the wish is to get to this intermediate point with a tolerance ΔT, this margin possibly being positive or negative depending on the case, then the current time can be corrected according to the following relation:

Tcorrected(X)=T(X)+ΔT. 

The speed setpoint to be followed to reach this intermediate waypoint is then given by the relation (15).
This speed setpoint must then be bounded by the flight range, therefore by Vmin(X) and Vmax(X). In this way, an intermediate objective, even it is not achievable, may be used to guide a more global speed strategy.
This principle can also be used to establish a succession of time-of-passage setpoints, and serve as a basis for observing multiple RTAs.
It should be recalled that, for the assumptions defining the speed profile to remain applicable, it is essential for the different speed profiles to be established on a common trajectory of the aircraft 100. The latter may therefore not be achievable in practice, according to the maximum speed profile Vmax(X), or may prove highly sub-optimal according to the minimum speed profile Vmin(X). However, this is not of great importance, since it is sufficient for the trajectory to remain practicable according to the speed profile determined by the method of the invention.
Furthermore, it should be recalled that it has hitherto been considered that, for the first RTA observance prediction calculation, and in the absence of an initial speed profile with which to observe the RTA, the trajectory would have had to be calculated by using the speed profile Vopt(X). Now, if the speed required at the start of the trajectory is too distant from Vopt(X), the consequence may be an impracticable or sub-optimal trajectory.
It is possible to correct, at least partially, this consequence, by anticipating the speed correction that will be necessary, and by biasing accordingly the speed value that will be used to calculate the lateral transitions of the trajectory. This is applied for the calculation of the lateral trajectory, after a first calculation of predictions according to the speed profile Vopt(X), which makes it possible to know the time TVopt(X) required to observe the RTA.
FIG. 4 shows examples of time-of-passage curves of an aircraft illustrating a speed correction calculation principle according to the invention, and provides a better understanding of the following description.
To this end, it is assumed that, according to each speed profile, the ETA changes linearly up to the point at which the RTA constraint is applied. According to this assumption, the position Xcapt of capture of the profile Vopt(X) can be approximated when the tolerance margin ΔT is applied to the RTA setpoint:




  
    
      
        
          
            Δ
            ⁢
            
              
            
            ⁢
            T
          
          =
          
            
              
              
                
                  T
                  ⁡
                  
                    (
                    X
                    )
                  
                
                -
                
                  
                    T
                    Vopt
                  
                  ⁡
                  
                    (
                    X
                    )
                  
                
              
              
            
            ⁢
            
              
                
                  X
                  capt
                
                -
                
                  X
                  RTA
                
              
              
                X
                -
                
                  X
                  capt
                
              
            
          
        
        ,
      
    
    
      
        (
        16
        )
      
    
  
  
    
      
        
          X
          capt
        
        =
        
          
            
              
                Δ
                ⁢
                
                  
                
                ⁢
                
                  T
                  ·
                  X
                
              
              +
              
                
                  X
                  RTA
                
                ⁡
                
                  (
                  
                    
                    
                      
                        T
                        ⁡
                        
                          (
                          X
                          )
                        
                      
                      -
                      
                        
                          T
                          
                            V
                            ⁢
                            
                              
                            
                            ⁢
                            opt
                          
                        
                        ⁡
                        
                          (
                          X
                          )
                        
                      
                    
                    
                  
                  )
                
              
            
            
              
                
                
                  
                    T
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                  -
                  
                    
                      T
                      
                        V
                        ⁢
                        
                          
                        
                        ⁢
                        opt
                      
                    
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                
                
              
              +
              
                Δ
                ⁢
                
                  
                
                ⁢
                T
              
            
          
          .
        
      
    
    
      
        (
        17
        )
      
    
  




It is then possible to linearly bias the speed of the profile Vopt(X), by a bias δV(X), which linearly approximates the speed correction that will be applied by the abovementioned relation (15) when calculating the speed setpoint in order to observe the RTA.

    
    
        If T(X)<TVopt(X)−Δthen:
    
    






  
    
      
        
          
            δ
            ⁢
            
              
            
            ⁢
            
              V
              ⁡
              
                (
                X
                )
              
            
          
          =
          
            
              
                V
                opt
              
              ⁡
              
                (
                X
                )
              
            
            ⁢
            
              
                
                  
                    T
                    Vopt
                  
                  ⁡
                  
                    (
                    X
                    )
                  
                
                -
                
                  (
                  
                    
                      T
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                    -
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                  
                  )
                
              
              
                
                  T
                  RTA
                
                -
                
                  (
                  
                    
                      T
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                    -
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                  
                  )
                
              
            
          
        
        ;
      
    
    
      
        (
        18
        )
      
    
  




    
    
        else, if T(X)>TVopt(X)+ΔT then:
    
    






  
    
      
        
          δ
          ⁢
          
            
          
          ⁢
          
            V
            ⁡
            
              (
              X
              )
            
          
        
        =
        
          
            
              V
              opt
            
            ⁡
            
              (
              X
              )
            
          
          ⁢
          
            
              
                
                  
                    T
                    Vopt
                  
                  ⁡
                  
                    (
                    X
                    )
                  
                
                -
                
                  (
                  
                    
                      T
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                    +
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                  
                  )
                
              
              
                
                  T
                  RTA
                
                -
                
                  (
                  
                    
                      T
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                    +
                    
                      Δ
                      ⁢
                      
                        
                      
                      ⁢
                      T
                    
                  
                  )
                
              
            
            .
          
        
      
    
    
      
        (
        19
        )
      
    
  




The relations (18) and (19) are only approximations, but they make it possible to limit the deviations between the speed profile that determines the trajectory, and the one that is actually applied in flight.
The remaining deviations can then be compensated for by the tolerance of the speed setpoint defined by the abovementioned relation (15), which naturally corrects the deviations and the biases on the predicted ETA, to make the predictions converge towards observance of the RTA.
Another advantage of the method according to the present invention is that it makes it possible to apply different strategies for determining speed setpoints along the flight plan, for example according to the different flight phases. It is, in practice, possible to determine strategies optimizing the use of the acceptable speed ranges according to the flight phase, or else according to portions of the flight plan on which the speed variations have a specific impact on an optimization criterion. To this end, it is possible to apply the speed setpoint formulae defined by the preceding relations, with judiciously selected margins and intermediate waypoints.
The description that follows presents possible choices of applicable strategies for determining speed setpoints, in the respectively typical climb, cruise and descent phases, according to the position of the RTA in the flight plan, and the time-related situation of the aircraft relative to reference time-related situations.
This description is given as an example, and is not limiting on the method according to the invention.
FIG. 5a presents a flow diagram illustrating the different strategies that can be applied in the climb phase according to different configurations, according to an exemplary embodiment of the invention.
FIG. 5b gives the time-of-passage curves as a function of the curvilinear abscissa, illustrating the different strategies that can be applied in the climb phase.
When the aircraft is in the climb phase, or CLB phase, it is, for example, preferable to adopt a strategy that makes it possible to accelerate as little as possible, high speeds being unfavourable to an effective climb. The reference speed profile may then correspond to the pair of optimum speeds determining the cross-over altitude.
In this case:

    
    
        if the RTA is situated ahead of the T/C point of abscissa XT/C and regardless of the time-related position of the aircraft, then the speed compensation can be distributed uniformly along the climb phase. The time of passage targeted for the point situated at the abscissa XRTA is the time TRTA. The tolerance margin ΔT can be set to half the current maximum time deviation δt, according to the time-related position of the aircraft relative to TVopt(X).
        otherwise, that is to say, if the RTA is situated at a point belonging to the cruise phase or the descent phase, then:
        
            if the time of passage T(X) of the aircraft 100 at the abscissa point X is less than the corresponding time TVopt(X) at the same point, in a first region 51 with reference to FIG. 5b, then it is possible for the aircraft 100 to fly according to the minimum speed profile Vmin(X) as long as T(X) remains less than TVopt(X), then to fly according to the optimum speed profile Vopt(X);
            if the time of passage T(X) is less than TVmax(XT/C)+TVopt(X)−TVopt(XT/C), in a second region 52 with reference to FIG. 5b, then the time-related situation indicates that, by climbing according to the speed profile Vopt(X), the aircraft 100 can ensure that the RTA is observed during the climb. In practice, the aircraft does not reach the point T/C later than the time TVmax. In this case, according to a first option, it is possible to fly according to the optimum speed profile Vopt(X). According to a second option, it is possible to seek to minimize the acceleration during the climb phase. Referring to FIG. 5b, the targeted waypoint then has for abscissa that of the T/C point according to a flight at the maximum speed profile Vmax(X), and for ordinate the corresponding time of passage. The tolerance ΔT that is applied can then be set to a negative value, the absolute value of which equals half the current maximum time deviation δt.
            if the time of passage T(X) is greater than TVmax(XT/C)+TVopt(X)−TVopt(XT/C), in a third region 53 with reference to FIG. 5b, then it is possible to seek to minimize the acceleration during the climb, but the situation then dictates flying faster than according to the optimum speed profile Vopt(X). Referring to FIG. 5b, the targeted waypoint then has for abscissa that of the T/C point according to a flight at the maximum speed profile Vmax(X), and for ordinate the corresponding time of passage TVmax(X). The tolerance ΔT that is applied can then be set at a negative value, the absolute value of which equals half the current maximum time deviation δt.
        
        
    
    


FIG. 6a is a flow diagram illustrating the different strategies that can be applied in the cruise phase according to different configurations, according to an exemplary embodiment of the invention.
FIG. 6b gives the time-of-passage curves as a function of the curvilinear abscissa, illustrating the different strategies that can be applied in the cruise phase.
When the aircraft is now in the cruise phase, or CRZ phase:

    
    
        if the RTA is situated ahead of the top of descent point, usually designated T/D, of abscissa XT/D, at a distance less than a determined distance from the aircraft 100, for example 60 Nm, and in the leg of the flight plan or current leg, then, regardless of the time-related position of the aircraft, the speed compensation can be distributed uniformly along the cruise phase. Referring to FIG. 6b, the targeted waypoint is the point situated at the abscissa XRTA and at the ordinate TRTA. The tolerance can be set to half the current tolerance, according to the time-related position of the aircraft 100 relative to TVopt(X);
        if the RTA is situated ahead of the T/D point of abscissa XT/D, for example at a distance greater than a determined distance of 60 Nm, or beyond the current leg, then regardless of the time-related position of the aircraft 100, the speed compensation can be distributed uniformly along the cruise phase. However, once the speed setpoint has been calculated, the latter is maintained along the current leg, over a distance at most equal to 60 Nm, before once again determining a new speed setpoint. Referring to FIG. 6b, the targeted time of passage is the point situated at the abscissa XRTA and at the ordinate TRTA. The tolerance ΔT can be set at half the current tolerance δt, according to the time-related position of the aircraft 100 relative to TVopt(X);
        if the RTA is situated beyond the T/D point of abscissa XT/D, then:
        
            if the time of passage T(X) of the aircraft 100 is less than the time TVopt(X), in a first region 61 with reference to FIG. 6b, then it is possible to seek to retain a speed margin in descent, therefore arrive at least at the time TVopt(XT/D) at the point T/D. However, in order to avoid an unnecessary overcompensation, it is possible not to seek to immediately join the flight profile corresponding to the optimum speed profile Vopt(X). Referring to FIG. 6b, the targeted waypoint then has for abscissa that of the T/D point and an ordinate equal to TVopt (XT/D). The applied tolerance ΔT can then be set at a positive value equal to half the current maximum time deviation δt;
            if the time of passage T(X) is less than TVmax(XT/D)+TVopt(X)−TVopt(XT/D), in a second region 62 with reference to FIG. 6b, then the time-related situation indicates that, by flying according to the speed profile Vopt(X), the aircraft can ensure that the RTA is observed while cruising. In practice, the aircraft does not reach the T/D point later than the time TVmax. Furthermore, it is preferable to favour a flight according to the optimum speed profile Vopt(X), rather than to risk an unnecessary speed overcompensation. In this case, it is possible to fly according to the speed profile Vopt(X).
            if the time of passage T(x) is greater than TVmax(XT/D)+TVopt(X)−TVopt(XT/D), in a third region 63 with reference to FIG. 6b, then it is possible to seek to provide a speed margin for the descent, therefore to minimize the acceleration, even though the time-related situation dictates authorizing flying faster than according to the optimum speed profile Vopt(X). The targeted waypoint then has for abscissa that of the T/D point according to a flight at the maximum speed profile Vmax(X), and for ordinate the corresponding time of passage. The applied tolerance ΔT can then be set at a negative value, the absolute value of which equals half the current maximum time deviation δt.
        
        
    
    


FIG. 7 is a flow diagram illustrating the different strategies that can be applied in the descent phase according to different configurations, according to an exemplary embodiment of the invention.
When the aircraft is in the descent phase, or, to use the usual designation, DES phase, it is possible to seek to avoid switching to underspeed, because the practicable flight gradients are shallower, and keeping to the descent profile is then no longer possible. However, it is necessary beforehand to define the speeds with which the theoretical speed profile must be determined. In practice, the descent predictions are made in two stages: firstly, the calculation of a theoretical descent profile, determined by starting from the arrival airport, and which determines the descent profile required to ensure approach and landing; then, a prediction calculation deriving from the current position and used to converge with the theoretical descent profile. Before determining the predictions to be applied during the descent, it is therefore necessary to determine the theoretical descent profile. This profile is generally based on a pair of setpoint speeds (one Mach setpoint and one air speed setpoint, usually designated by the acronym CAS, standing for Computed Air Speed) applicable to the entire descent as long as no speed constraint is imposed by the arrival procedure. The descent is then made at constant Mach, according to the Mach setpoint, for the high altitudes, then at constant CAS, according to the CAS setpoint, for the low altitudes. The transition altitude between these two setpoints is called cross-over altitude, for which the CAS and Mach speeds correspond to the same physical speed (usually designated by the acronym TAS, standing for True Air Speed).
In order to avoid making an iterative calculation on the speeds of the theoretical profile during the descent phase, it is possible to determine a performance index, designated PI, or any other type of index, for example a cost-type index. This index PI is then used to determine the speeds to be applied for the theoretical descent profile; this theoretical profile is not recalculated during the descent, except in the event of a modification of the flight plan.
To determine the PI of the descent, the first step is, for example, to calculate the theoretical descent profiles at Vopt, Vmin and Vmax. The required times of passage TVmin(X), TVopt(X) and TVmax(X) are then deduced therefrom. The performance index to be applied is then determined by comparing the current ETA with these required times of passage. A PI with a value of between −100 and 0 corresponds to the relative positioning of the ETA between TVmin(X) and TVopt(X), and is reflected in a speed that is situated in the same proportion between Vmin and Vopt. Similarly, a PI with a value of between 0 and 100 corresponds to the relative positioning of the ETA between TVopt(X) and TVmax(X), and is reflected in a speed that is situated in the same proportion between Vopt and Vmax. A negative tolerance ΔT can be applied to the ETA, equal to the current maximum time deviation δt, in order to impose shallower gradients on the theoretical profile.
Thus, if the current ETA minus ΔT is less than TVopt(X):




  PI
  =
  
    100
    ⁢
    
      
        T
        -
        
          Δ
          ⁢
          
            
          
          ⁢
          T
        
        -
        
          T
          Vopt
        
      
      
        
          T
          Vopt
        
        -
        
          T
          
            V
            ⁢
            
              
            
            ⁢
            min
          
        
      
    
  




Also, the CAS and Mach speeds applied during the descent will be determined with the same ratio:




  V
  =
  
    
      PI
      100
    
    ⁢
    
      
        (
        
          
            V
            min
          
          -
          
            V
            opt
          
        
        )
      
      .
    
  




    
    
        Similarly, if the current ETA minus ΔT is greater than TVopt(X):
    
    






  PI
  =
  
    100
    ⁢
    
      
        T
        -
        
          Δ
          ⁢
          
            
          
          ⁢
          T
        
        -
        
          T
          Vopt
        
      
      
        
          T
          
            V
            ⁢
            
              
            
            ⁢
            max
          
        
        -
        
          T
          
            V
            ⁢
            opt
          
        
      
    
  




Also, the CAS and Mach speeds applied during the descent will be determined with the same ratio:




  V
  =
  
    
      PI
      100
    
    ⁢
    
      (
      
        
          V
          max
        
        -
        
          V
          opt
        
      
      )
    
  




Once the theoretical descent profile has been calculated with this pair of optimum speeds, the predictions are applied relative to this theoretical reference profile, by using the following speed strategy.

    
    
        Regardless of the position of the RTA relative to the current point, and regardless of the time-related situation, the strategy in the descent phase consists in compensating for the time deviations noted in order to converge towards the time constraint. The targeted time of passage is the point situated at the abscissa XRTA and at the ordinate TRTA. The tolerance ΔT can be set at half the current maximum time deviation δt, according to the time-related position of the aircraft relative to TVopt(X).

## Claims

1. A method for generating a flight speed profile of an aircraft according to an RTA constraint, using a flight management system, said method comprising the following steps:
for a determined current position of the aircraft along a flight plan, determination of at least two reference speed profiles according to the curvilinear abscissa X of the aircraft along the flight plan, comprising a minimum speed profile V min(X) and a maximum speed profile Vmax(X),
storage of a plurality of pairs (Tvmin(Xi),Xi) and (Tvmax(Xi),Xi) of values of times of passage of the aircraft at curvilinear abscissa points Xi subsequent to the current position as far as the abscissa of the point at which the RTA isapplied for the reference speed profiles,
application at any point of curvilinear abscissa X of the flight plan of a speed setpoint VΩ(X) of the aircraft determined on the basis of the speed profiles V min(X) and V max(X) weighted by proportionality ratios between T(X), Tvmin(X) and Tvmax(X) with flight management system;
T(X) designating the time of passage of the aircraft at the point of the flight plan of curvilinear abscissa X, Tvmax(X) designating the time of passage at the curvilinear abscissa point X, for a flight according to the minimum speed profile Vmin(X) reaching the point at which the RTA is applied by observing the RTA constraint, and Tvmin(X) designating the time of passage at the curvilinear abscissa point X, for a flight according to the maximum speed profile Vmax(X) reaching the point at which the RTA is applied by observing the RTA constraint, the values of Tvmin(X) and Tvmax(X) being determined by a linear interpolation based on the corresponding stored pairs of values for the closest points either side of the curvilinear abscissa point X.

2. The method for generating a flight speed profile of an aircraft according to claim 1, wherein, for the current position of the aircraft along the flight plan, an additional optimum speed reference speed profile Vopt(X) is determined according to at least one defined flight criterion, the speed profile of the aircraft being constrained to tend towards the optimum speed profile Vopt(X) by the addition to the RTA of a tolerance ΔT the sign of which is opposite to the sign of the difference between the time of passage of the aircraft at the current position and the corresponding time of passage Tvopt(X), until the absolute value of the difference between T(X) and Tvopt(X) is less than a maximum time difference δt, the RTA then being restored to its initial value.

3. The method for generating a flight speed profile of an aircraft according claim 2, wherein said at least one defined flight criterion comprises a cost index.

4. The method for generating a flight speed profile of an aircraft according to claim 2, wherein said at least one defined flight criterion comprises a performance index.

5. The method for generating a flight speed profile of an aircraft according to claim 1, wherein the speed setpoint VΩ(X) is determined by the relation:




  
    
      V
      Ω
    
    ⁡
    
      (
      X
      )
    
  
  =
  
    
      
        
          (
          
            
              
                T
                
                  V
                  ⁢
                  
                    
                  
                  ⁢
                  max
                
              
              ⁡
              
                (
                X
                )
              
            
            -
            
              
                T
                
                  V
                  ⁢
                  
                    
                  
                  ⁢
                  min
                
              
              ⁡
              
                (
                X
                )
              
            
          
          )
        
        ⁢
        
          
            V
            min
          
          ⁡
          
            (
            X
            )
          
        
        ⁢
        
          
            V
            max
          
          ⁡
          
            (
            X
            )
          
        
      
      
        
          
            (
            
              
                
                  T
                  
                    V
                    ⁢
                    
                      
                    
                    ⁢
                    max
                  
                
                ⁡
                
                  (
                  X
                  )
                
              
              -
              
                T
                ⁡
                
                  (
                  X
                  )
                
              
            
            )
          
          ⁢
          
            
              V
              max
            
            ⁡
            
              (
              X
              )
            
          
        
        -
        
          
            (
            
              
                
                  T
                  
                    V
                    ⁢
                    
                      
                    
                    ⁢
                    min
                  
                
                ⁡
                
                  (
                  X
                  )
                
              
              -
              
                T
                ⁡
                
                  (
                  X
                  )
                
              
            
            )
          
          ⁢
          
            
              V
              min
            
            ⁡
            
              (
              X
              )
            
          
        
      
    
    .

6. The method for generating a flight speed profile of an aircraft according claim 1, wherein the speed setpoint VΩ(X) is determined by the relation:




  
    
      
        V
        Ω
      
      ⁡
      
        (
        X
        )
      
    
    =
    
      
        
          (
          
            
              
                T
                
                  V
                  ⁢
                  
                    
                  
                  ⁢
                  max
                
              
              ⁡
              
                (
                X
                )
              
            
            -
            
              
                T
                min
              
              ⁡
              
                (
                X
                )
              
            
          
          )
        
        ⁢
        
          
            V
            
              V
              ⁢
              min
            
          
          ⁡
          
            (
            X
            )
          
        
        ⁢
        
          
            V
            max
          
          ⁡
          
            (
            X
            )
          
        
      
      
        
          
            
              
                
                  (
                  
                    
                      
                        T
                        
                          V
                          ⁢
                          
                            
                          
                          ⁢
                          max
                        
                      
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                    -
                    
                      
                        T
                        corrected
                      
                      ⁡
                      
                        (
                        X
                        )
                      
                    
                  
                  )
                
                ⁢
                
                  
                    V
                    max
                  
                  ⁡
                  
                    (
                    X
                    )
                  
                
              
              -
            
          
        
        
          
            
              
                (
                
                  
                    
                      T
                      
                        V
                        ⁢
                        
                          
                        
                        ⁢
                        min
                      
                    
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                  -
                  
                    
                      T
                      corrected
                    
                    ⁡
                    
                      (
                      X
                      )
                    
                  
                
                )
              
              ⁢
              
                
                  V
                  min
                
                ⁡
                
                  (
                  X
                  )
                
              
            
          
        
      
    
  
  ,




Tcorrected(X) being equal to T(X)−δt if T(X) is less than Tvopt(X)−δt, to T(X)+δt if T(X) is greater than Tvopt(X)+δt, or to T(X) otherwise.

7. The method for generating a flight speed profile of an aircraft according to claim 1, wherein the speed setpoint VΩ(X) is corrected by a value δV(X) equal to




  
    
      V
      opt
    
    ⁡
    
      (
      X
      )
    
  
  ⁢
  
    
      
        
          T
          Vopt
        
        ⁡
        
          (
          X
          )
        
      
      -
      
        (
        
          
            T
            ⁡
            
              (
              X
              )
            
          
          -
          
            Δ
            ⁢
            
              
            
            ⁢
            T
          
        
        )
      
    
    
      
        T
        RTA
      
      -
      
        (
        
          
            T
            ⁡
            
              (
              X
              )
            
          
          -
          
            Δ
            ⁢
            
              
            
            ⁢
            T
          
        
        )
      
    
  




if T(X) is less than Tvopt(X)−ΔT, or to




  
    
      V
      opt
    
    ⁡
    
      (
      X
      )
    
  
  ⁢
  
    
      
        
          T
          Vopt
        
        ⁡
        
          (
          X
          )
        
      
      -
      
        (
        
          
            T
            ⁡
            
              (
              X
              )
            
          
          +
          
            Δ
            ⁢
            
              
            
            ⁢
            T
          
        
        )
      
    
    
      
        T
        RTA
      
      -
      
        (
        
          
            T
            ⁡
            
              (
              X
              )
            
          
          +
          
            Δ
            ⁢
            
              
            
            ⁢
            T
          
        
        )
      
    
  




if T(X) is greater than Tvopt(X)+ΔT.

8. The method for generating a flight speed profile of an aircraft according to claim 1, wherein said method generates a plurality of speed profiles so as to observe a plurality of successive time of passage constraints along the flight plan.

9. The method for generating a flight speed profile of an aircraft according claim 1, wherein a plurality of different speed profiles are generated for different flight phases for which the optimum speed profiles are determined on the basis of different flight criteria.

10. A flight management system of an aircraft, wherein the system implements the meyhod according to claim 1.

11. The flight management system of an aircraft according to claim 10, wherein the flight management system calculates, at any instant of the flight of the aircraft, a speed setpoint to be applied to the flight management system.

12. The flight management system of an aircraft according to claim 10, wherein the speed setpoint to be applied is displayed on a display interface.

13. The flight management system of an aircraft according to claim 10, wherein the speed setpoint to be applied is transmitted to an automatic pilot system.

14. A method of generating a flight speed profile of an aircraft according to an RTA constraint by a flight management system, said method comprising:
determining for a current position of the aircraft along a flight plan, at least two reference speed profiles according to a curvilinear abscissa X of the aircraft along the flight plan by the flight management system, comprising a minimum speed profile Vmin(X) and a maximum speed profile Vmax(X);
storing a plurality of pairs (Tvmin(Xi),Xi) and (Tvmax(Xi),Xi) values of times of passage of the aircraft at curvilinear abscissa points Xi subsequent to the current position along the abscissa at the point at which the RTA is applied for the reference speed profiles; and
applying at any point of the curvilinear abscissa X of the flight plan a speed setpoint VΩ(X) of the aircraft determined based on the speed profiles Vmin(X) and Vmax(X) weighted by proportionality ratios between T(X), Tvmin(X) and Tvmax(X) with the flight management system,
wherein T(X) designates the time of passage of the aircraft at the point of the flight plan of the curvilinear abscissa X, Tvmax(X) designates the time of passage at the curvilinear abscissa point X, for a flight according to the minimum speed profile Vmin(X) reaching the point at which the RTA is applied by observing the RTA constraint, and Tvmin(X) designates the time of passage at the curvilinear abscissa point X, for a flight according to the maximum speed profile Vmax(X) reaching the point at which the RTA is applied by observing the RTA constraint, the values of Tvmin(X) and Tvmax(X) being determined by a linear interpolation based on the corresponding stored pairs of values for the closest points on either side of the curvilinear abscissa point X.