Fast vertical trajectory prediction method for air traffic management, and relevant ATM system
IPC분류정보
국가/구분
United States(US) Patent
등록
국제특허분류(IPC7판)
G08G-005/00
G05D-001/06
출원번호
US-0324047
(2011-12-13)
등록번호
US-9031720
(2015-05-12)
우선권정보
IT-RM2010A0672 (2010-12-20)
발명자
/ 주소
Accardo, Domenico
Moccia, Antonio
Grassi, Michele
Tancredi, Urbano
Caminiti, Lucio
Fiorillo, Luigi
Leardi, Alberto
Maresca, Giuseppe
출원인 / 주소
Selex Sistemi Integrati S.p.A.
대리인 / 주소
Ladas & Parry LLP
인용정보
피인용 횟수 :
0인용 특허 :
10
초록▼
A method and system for the prediction of aircrafts vertical trajectory, in particular for Air Traffic Management, using the following flight calculation modules: Take-off; Climb; Cruise; Descent; and Landing, corresponding to the relevant flight phases, in which: the calculation of the predicted ai
A method and system for the prediction of aircrafts vertical trajectory, in particular for Air Traffic Management, using the following flight calculation modules: Take-off; Climb; Cruise; Descent; and Landing, corresponding to the relevant flight phases, in which: the calculation of the predicted aircraft trajectory is effected by using a set of TEM equations using, as output variables, the vertical rate of climb or descent, the true air speed, the energy share factor, the thrust and the drag, the mass of the aircraft modeled as point-mass, and using, as input variables, the Mach number depending on true air speed and temperature and altitude, the gravity acceleration, and the fuel flow, and the flight path angle;the calculation of the predicted aircraft trajectory for Cruise phase, wherein only the mass is varies, is performed by using the following analytical solution to said set of TEM equations.
대표청구항▼
1. A computer-implemented method for controlling an aircraft's flight plan and guiding an aircraft based on prediction of the aircraft's vertical trajectory, comprising the following flight calculation modules: Takeoff; Climb; Cruise; Descent; and Landing, corresponding to relevant flight phases, wh
1. A computer-implemented method for controlling an aircraft's flight plan and guiding an aircraft based on prediction of the aircraft's vertical trajectory, comprising the following flight calculation modules: Takeoff; Climb; Cruise; Descent; and Landing, corresponding to relevant flight phases, wherein the following steps are executed: calculating by a computer processor a predicted aircraft trajectory using the following total energy model (TEM) equations: VRCD=(T-D)mgTASESF{M}mⅆTASⅆt=(T-D)-mgVRCDTAS=(T-D)(1-ESF{M})m.=-fⅆhⅆt=TASsinγ solving said equations for VRCD, TAS, and m where VRCD is vertical rate of climb or descent; TAS is true air speed, ESF is the energy share factor, T is thrust and D drag, m the mass of the aircraft modeled as point-mass, {M} is the Mach number depending on TAS and temperature and altitude, g is gravity acceleration, and f is fuel flow, and γ is the flight path angle; calculating the predicted aircraft trajectory for Cruise phase, wherein only the mass is varying, by using the following analytical solution to said TEM equations: t-t0=distTAS=1k9·k10·[tan-1(k10k9·mfin)-tan-1(k10k9·min)] solved for mass mfin at the end of the cruise phase as a function of initial mass min, and wherein t is elapsed flight time, k9 and k10 are constant terms predefined according to the aircraft sending the predicted aircraft trajectory to an Air Traffic Management (ATM) controller, which decides the flight plan accordingly;communicating, by the ATM controller the flight plan to the aircraft;guiding the aircraft to follow the flight plan; the TEM equation being solved:each time a flight plan is needed or changed;each time, within a fixed flight plan, the difference between the actual position of the aircraft, given by a radar detection, and the predicted position is greater than a pre-defined value. 2. Method according to claim 1, wherein for the Take-off phase, which is divided into ground roll, transition and initial climb phases, the calculation of the predicted aircraft trajectory is performed by using the following analytical solutions to said TEM equations: tTR=tLO+RTRVLOγTR;x(tTR)=xTR;h(tTR)=hTR;V(tTR)=VLO=1.2·(Vstall)TOⅆxⅆt=V1-(T-D)2W2ESF2·Cpow,red2 that are solved by the ground travelled distance x, wherein tTR is the transition phase time, tLO is the exact time of lift-off, xTR, the travelled distance at tTR, γTR the travelled angle during transition, VLO the lift-off TAS, h is the altitude, hTR the altitude at tTR, V is the current TAS, (Vstall)TO is the stall speed of the aircraft with gear down, W is the aircraft weight, Cpow,red is the pre-defined coefficient of power reduction. 3. Method according to claim 1, wherein for the Landing phase, the calculation of the predicted aircraft trajectory is effected by using the following analytical solutions: VC=(h-h0)VCA,screen+(hscreen-h)VC0hscreen-ho for the glide approach, and ⅆhⅆt=Vflaresinγ for the flare, wherein h is the altitude, h0 is the altitude for the beginning of glide approach, hscreen is the decision altitude, VCA, screen is the CAS speed to be reached at hscreen, VC0 is the CAS speed at h0, VC the CAS speed at h; Vflare the CAS speed during flare and gamma is the slope of the vertical trajectory during flare. 4. Method according to claim 1, wherein the integration of the TEM equations for the calculation of predicted trajectory is made by using a pair of maximum integration pitches for speed and height, in order to address the minimum computational load at an acceptable accuracy level, the maximum integration pitches pair being determined by performing the following steps: Performing simulations, according to said TEM equations, of climb, descent, and cruise phases for uniformly distributed set of pairs of speed and height pitches ranging from a minimum values pair to a maximum values pair:Assuming the minimum values pair as the most accurate values pair;For each simulation, comprising climb, descent, and cruise phases, carrying out contour plots reporting the percent RMS error of each pair of speed and height pitches with respect to said minimum values pair;Choosing the optimal pitches pair as the pair representing the point that has an error of less than a pre-defined threshold value and it is also the most distant from said minimum values pair. 5. Method according to claim 1, wherein for all the flight phases except Cruise, the TEM equations are integrated and, for any i-th, i being a positive integer number, integration step: one checks that the calculated performance status is within a predefined target PS, comprised of a target CAS and target altitude h, calculated on the basis of a pre-defined flight envelope;if the calculated performance status is outside the flight envelope, performing the following steps: substituting said calculated performance status with a corrected performance status that is nearest to the boundaries of the flight envelope and to which a safe margin distance from these boundaries is added, in order to avoid that in the subsequent calculation it goes outside the flight envelope;proceeding to the i+1-th step of integration starting with the corrected CAS and altitude h. 6. Method according to claim 1, wherein the effect of the wind is taken into account by adding the following equations to said TEM equations: GSLong=WSLong+TAS2-WSLat2-VRCD22ψ=β-arctg(-WSLatTASLong)=β-arctg(-WSLatGSLong-WSLong)==β-arctg(-WSLatTAS2-WSLat2-VRCD22) And solving for GSLong, that is the horizontal component of the aircraft ground speed, and for □ that is the heading angle, wherein WSLong is the horizontal wind speed, WSLat is the lateral component of wind speed, beta is the course angle. 7. Method according to claim 2, wherein for the Landing phase, the calculation of the predicted aircraft trajectory is effected by using the following analytical solutions: VC=(h-h0)VCA,screen+(hscreen-h)VC0hscreen-h0 for the glide approach, and ⅆhⅆt=Vflaresinγ for the flare, wherein h is the altitude, h0 is the altitude for the beginning of glide approach, hscreen is the decision altitude, VCA, screen is the CAS speed to be reached at hscreen, VC0 is the CAS speed at h0, VC the CAS speed at h; Vflare the CAS speed during flare and □ is the slope of the vertical trajectory during flare. 8. Method according to claim 1, characterized in that the TEM equations are solved: before the take-off phase of the aircraft. 9. Method according to claim 2, wherein the integration of the TEM equations for the calculation of predicted trajectory is made by using a pair of maximum integration pitches for speed and height, in order to address the minimum computational load at an acceptable accuracy level, the maximum integration pitches pair being determined by performing the following steps: Performing simulations, according to said TEM equations, of climb, descent, and cruise phases for uniformly distributed set of pairs of speed and height pitches ranging from a minimum values pair to a maximum values pair:Assuming the minimum values pair as the most accurate values pair;For each simulation, comprising climb, descent, and cruise phases, carrying out contour plots reporting the percent RMS error of each pair of speed and height pitches with respect to said minimum values pair;Choosing the optimal pitches pair as the pair representing the point that has an error of less than a pre-defined threshold value and it is also the most distant from said minimum values pair. 10. Method according to claim 2, wherein for all the flight phases except Cruise, the TEM equations are integrated and, for any i-th, i being a positive integer number, integration step: one checks that the calculated performance status is within a predefined target PS, comprised of a target CAS and target altitude h, calculated on the basis of a pre-defined flight envelope; if the calculated performance status is outside the target PS, performing the following steps: substituting said calculated performance status with a corrected performance status determined by assigning a safe margin distance to the nearest point within the flight envelope; proceeding to the i+1-th step of integration starting with the corrected performance status, i.e. corrected CAS and altitude h. 11. Method according to claim 2, wherein the effect of the wind is taken into account by adding the following equations to said TEM equations: GSLong=WSLong+TAS2-WSLat2-VRCD22ψ=β-arctg(-WSLatTASLong)=β-arctg(-WSLatGSLong-WSLong)==β-arctg(-WSLatTAS2-WSLat2-VRCD22) And solving for GSLong, that is the horizontal component of the aircraft ground speed, and for □ that is the heading angle, wherein WSLong is the horizontal wind speed, WSLat is the lateral component of wind speed, beta is the course angle. 12. Method according to claim 3, wherein the integration of the TEM equations for the calculation of predicted trajectory is made by using a pair of maximum integration pitches for speed and height, in order to address the minimum computational load at an acceptable accuracy level, the maximum integration pitches pair being determined by performing the following steps: Performing simulations, according to said TEM equations, of climb, descent, and cruise phases for uniformly distributed set of pairs of speed and height pitches ranging from a minimum values pair to a maximum values pair:Assuming the minimum values pair as the most accurate values pair;For each simulation, comprising climb, descent, and cruise phases, carrying out contour plots reporting the percent RMS error of each pair of speed and height pitches with respect to said minimum values pair;Choosing the optimal pitches pair as the pair representing the point that has an error of less than a pre-defined threshold value and it is also the most distant from said minimum values pair. 13. Method according to claim 3, wherein for all the flight phases except Cruise, the TEM equations are integrated and, for any i-th, i being a positive integer number, integration step: one checks that the calculated performance status is within a predefined target PS, comprised of a target CAS and target altitude h, calculated on the basis of a pre-defined flight envelope;if the calculated performance status is outside the target PS, performing the following steps: substituting said calculated performance status with a corrected performance status determined by assigning a safe margin distance to the nearest point within the flight envelope; proceeding to the i+1-th step of integration starting with the corrected performance status, i.e. corrected CAS and altitude h. 14. Method according to claim 3, wherein the effect of the wind is taken into account by adding the following equations to said TEM equations: GSLong=WSLong+TAS2-WSLat2-VRCD22ψ=β-arctg(-WSLatTASLong)=β-arctg(-WSLatGSLong-WSLong)==β-arctg(-WSLatTAS2-WSLat2-VRCD22) and solving for GSLong, that is the horizontal component of the aircraft ground speed, and for ψ that is the heading angle, wherein WSLong is the horizontal wind speed, WSLat is the lateral component of wind speed, beta is the course angle. 15. Method according to claim 4, wherein for all the flight phases except Cruise, the TEM equations are integrated and, for any i-th, i being a positive integer number, integration step: one checks that the calculated performance status is within a predefined target PS, comprised of a target CAS and target altitude h, calculated on the basis of a pre-defined flight envelope;if the calculated performance status is outside the target PS, performing the following steps:substituting said calculated performance status with a corrected performance status determined by assigning a safe margin distance to the nearest point within the flight envelope;proceeding to the i+1-th step of integration starting with the corrected performance status, i.e. corrected CAS and altitude h. 16. Method according to claim 4, wherein the effect of the wind is taken into account by adding the following equations to said TEM equations: GSLong=WSLong+TAS2-WSLat2-VRCD22ψ=β-arctg(-WSLatTASLong)=β-arctg(-WSLatGSLong-WSLong)==β-arctg(-WSLatTAS2-WSLat2-VRCD22) and solving for GSLong, that is the horizontal component of the aircraft ground speed, and for Ψ that is the heading angle, wherein WSLong is the horizontal wind speed, WSLat is the lateral component of wind speed, beta is the course angle. 17. Method according to claim 5, wherein the effect of the wind is taken into account by adding the following equations to said TEM equations: GSLong=WSLong+TAS2-WSLat2-VRCD22ψ=β-arctg(-WSLatTASLong)=β-arctg(-WSLatGSLong-WSLong)==β-arctg(-WSLatTAS2-WSLat2-VRCD22) and solving for GSLong, that is the horizontal component of the aircraft ground speed, and for □that is the heading angle, wherein WSLong is the horizontal wind speed, WSLat, is the lateral component of wind speed, beta is the course angle. 18. A method according to claim 1, wherein the results of the solutions of said equations are displayed graphically. 19. Method according to claim 1, wherein the results obtained are employed in air traffic control. 20. System for Air Traffic Management, comprising a computer processor for the calculation of the predicted trajectory and a controller of the flight providing to said computer processor natural language commands for the trajectory the aircraft has to follow, wherein said computer processor executes a translation of the natural language commands into input numerical values, and in that said computer processor executes the calculation of the method according to claim 1. 21. Computer processor, which comprises code means that execute, when run, the method according to claim 1. 22. A nontransitory tangible memory medium, readable by a computer, storing a program, wherein the program is the computer program according to claim 20.
연구과제 타임라인
LOADING...
LOADING...
LOADING...
LOADING...
LOADING...
이 특허에 인용된 특허 (10)
Roberts, Alison Laura Udal; Pember, Stephen James, Air traffic control.
Karthikeyan, Tharangini; Kumaar, Kavitha; Palaniswamy, Kalpana, System and method for informing an aircraft operator about a temporary flight restriction in perspective view.
Lewis, Michael S.; Holland, Reza T., Systems and methods for real-time conflict-checked, operationally preferred flight trajectory revision recommendations.
※ AI-Helper는 부적절한 답변을 할 수 있습니다.