WO2009083507A1 - Procede pour l'aide a l'atterrissage d'aeronef utilisant un gps et un mls dans le cadre d'une approche axiale calculee - Google Patents
Procede pour l'aide a l'atterrissage d'aeronef utilisant un gps et un mls dans le cadre d'une approche axiale calculee Download PDFInfo
- Publication number
- WO2009083507A1 WO2009083507A1 PCT/EP2008/068084 EP2008068084W WO2009083507A1 WO 2009083507 A1 WO2009083507 A1 WO 2009083507A1 EP 2008068084 W EP2008068084 W EP 2008068084W WO 2009083507 A1 WO2009083507 A1 WO 2009083507A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- aircraft
- mls
- point
- coordinates
- azimuth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Ceased
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/005—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 with correlation of navigation data from several sources, e.g. map or contour matching
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
- G01S19/15—Aircraft landing systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/51—Relative positioning
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/04—Control of altitude or depth
- G05D1/06—Rate of change of altitude or depth
- G05D1/0607—Rate of change of altitude or depth specially adapted for aircraft
- G05D1/0653—Rate of change of altitude or depth specially adapted for aircraft during a phase of take-off or landing
- G05D1/0676—Rate of change of altitude or depth specially adapted for aircraft during a phase of take-off or landing specially adapted for landing
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft
- G08G5/20—Arrangements for acquiring, generating, sharing or displaying traffic information
- G08G5/22—Arrangements for acquiring, generating, sharing or displaying traffic information located on the ground
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft
- G08G5/50—Navigation or guidance aids
- G08G5/54—Navigation or guidance aids for approach or landing
Definitions
- Each station transmits a narrow beating beam sweeping back and forth at constant angular velocity the coverage space along the angular coordinate considered.
- An aircraft antenna and receiver receive the flying beam a first time in the forward scan and a second time in the reverse scan. It is thus possible to determine the azimuth angle and the elevation angle by the following linear relationship:
- ⁇ ( ⁇ - T o ). ⁇ (1) where ⁇ is the azimuth angle or the elevation angle, T is a time interval between the reception of the forward and return passages of the beating beam,
- T 0 the value of the time interval T for a zero angle ⁇ and v the scanning angular velocity.
- T 0 and v are constants defined by international standards on MLS.
- the calculated axial approach microwave landing system is an MLS in which the azimuth station is not placed at the end of the runway. but deported to one side of the airstrip.
- the displacement of the azimuth station can in particular be used in two cases. In the first case, the azimuth station is located near the site station for simplicity of MLS deployment. This configuration is encountered mainly in the case of tactical equipment deployed on summary and unprepared terrain. In the second case, the MLS-cc is used for the approach on a secondary track not equipped with an MLS but located in the coverage area of the runway equipped with MLS.
- the receiver of the aircraft measures an angle, called real azimuth angle, which does not correspond to the azimuth angle in the conventional sense of the term, called angle of virtual azimuth. It is therefore necessary to calculate the virtual azimuth angle to be able to provide the pilot information recentered with respect to the runway axis. To evaluate it, it is necessary to calculate the position of the aircraft in a Cartesian coordinate system centered on a reference point on the ground, called the ground point. This point on the ground is for example one of the two stations or the point of intersection between the runway axis and a straight line perpendicular to the runway axis and passing through the site station. This intersection point is called MLS datum point.
- the calculation of the position of the aircraft is performed by a system of 3 equations with 3 unknowns, parameterized by the real azimuth angle, the elevation angle and a distance between the aircraft and the ground point. Since these equations are non-linear, iterative algorithms are used to solve the system. Classically, the iterative algorithms are of the Gauss-Seidel or Newton-Raphson type. By using a satellite positioning system, it is possible to use the MLS datum point as a ground point for determining the distance with the aircraft. However, there are risks associated with iterative algorithm behaviors in this situation, as discussed in DO-226.
- the iterative algorithms are executed by the receiver of the aircraft, for example a multimode receiver, called MMR, acronym for the English expression "MultiMode Receiver”.
- MMR comprises a radio frequency chain, a digital signal processor, called DSP, acronym for the English expression “Digital Signal Processor”, a receiver of global positioning system, better known by the acronym “GPS” for the Anglo-Saxon expression “Global Positioning System” and a microprocessor.
- DSP digital signal processor
- GPS Global Positioning System
- the radio frequency chain receives signals from different systems, in particular the guidance signals from the azimuth station and the site station.
- the processor processes the signals from the radio frequency chain in order to extract the azimuth and site angles as well as the auxiliary words contained in the guidance signals.
- the GPS receiver can be internal or external. It transmits the positions of the aircraft.
- the microprocessor performs several functions including a calculation of the distance between the MLS datum point and the aircraft, a calculation of the position of the aircraft from the iterative algorithms, a calculation of deviations between the position of the aircraft and an axis. optimal descent, equipment management and communication interfacing with an aircraft bus link, for example an ARINC bus, abbreviation for the "Aeronautical Radio Incorporated" series of standards.
- the microprocessor includes an algorithm management device for detecting and processing the divergence, non-convergence or false convergence of an algorithm, managing the initialization, sequencing and combination of the algorithms.
- An object of the invention is in particular to overcome all or part of the aforementioned drawbacks.
- the subject of the invention is a method for assisting the approach and landing of aircraft using a satellite positioning system and an MLS system in the context of a calculated axial approach.
- MLS system comprising an airstrip, a site station and an azimuth station, the landing strip defining a direct Cartesian landmark ⁇ R c ) including as axes a longitudinal axis of the airstrip, called first axis (x) , a transverse axis of the landing runway passing through the site station, called second axis (y) and an axis perpendicular to the first and second axes, called third axis (z), a point of intersection between the first axis (x ) and the second axis (y) defining the center of the marker (R c ), called the runway point, the site station comprising a site antenna and the azimuth station comprising an azimuth antenna, each antenna comprising a phase center, the center of antenna phase if te and azimuth being located in a plane parallel to a plane ⁇ , y) at a distance dz from the track point along the third axis (z), characterized in that it comprises the following steps: create
- the invention has the particular advantage that it makes it possible to reduce the complexity of calculating the position of the aircraft, since the calculation no longer uses a Gauss-Seidel or Newton-Raphson type algorithm, but a simple search algorithm for the aircraft. roots of a polynomial of the second degree, as described in the DO-198 standard. Therefore, the algorithm is less complex to certify and the problems related to the use of iterative algorithms are removed. In particular, divergence and stationary problems, error propagation and multiple solutions are eliminated. In addition, the algorithm is less expensive in terms of computing load.
- FIG. 1 an exemplary configuration of FIG. MLS with calculated axial approach
- FIG. 4 an example of steps implemented by the method according to the invention.
- FIG. 5 an example of sub-steps for determining a distance between a reference point and the aircraft
- FIG. 6 an example of an airstrip, reference point, azimuth station and aircraft configuration for determining the distance between the reference point and the aircraft;
- FIG. 7 an example of sub-steps for determining a distance between an intermediate point and the aircraft;
- - Figure 8 a geographical reference and a geocartésien landmark;
- Figure 1 shows an example of MLS configuration as part of a calculated axial approach.
- the MLS includes a runway center runway 1 x oriented from an end of runway 2 threshold to the runway start threshold 3, a site station 4 and an azimuth station 5.
- the runway 1 defines a direct Cartesian coordinate system R c of axes the longitudinal axis of the track, called the first axis x, a transverse axis of the landing runway 1 passing through the site station 4, called second axis y and an axis perpendicular to the first and second axis, called third axis z.
- the second axis is oriented from the site station 4 to the landing strip 1.
- a point of intersection between the first axis x and the second axis y defines the center of the reference mark R c and is called MLS datum point 6. It generally corresponds to the touchdown point of the aircraft on the runway 1.
- the site station 4 includes a site antenna 7 and the azimuth station 5 comprises an antenna azimuth 8.
- Each antenna comprises a phase center from which the coordinates of the antenna can be determined. For the rest of the description, we will assimilate the position of each antenna to that of its phase center.
- the antennas site 7 and azimuth 8 emit guidance signals by swinging beams 9 and 10 for vertical and lateral guidance of an aircraft 1 1 during an approach phase and / or landing.
- FIG. 3 represents a real azimuth angle ⁇ R and a virtual azimuth angle ⁇ v .
- the real azimuth angle ⁇ R is located in a Cartesian coordinate system R "centered on the azimuth antenna 8 and with axes an axis x", the axis y and a z axis ", the axes x" and z “ being parallel to the axes x and z, respectively, and is defined by the angle between the plane ⁇ ", z”) and a line passing through the azimuth antenna 8 and the aircraft 11.
- the virtual azimuth angle ⁇ v is located in the Cartesian coordinate system R c and defined by the angle between the plane (x, z) and a line passing through the aircraft 1 1 and a virtual point 8v of the plane ( ⁇ , z). to an azimuth antenna in the context of an MLS system, usually located on the first x-axis at a distance of 300 m after the end-of-runway threshold 2.
- Cartesian coordinate system R c centered on the MLS datum point 6
- Cartesian coordinate system R c ' centered on the site antenna 7
- Cartesian coordinate system R "centered on the azimuth antenna 8 are considered.
- the geocartean landmark R geocart ( ⁇ , X, Y, z) is a landmark in which the center O is close to the Earth's center of mass, the OZ axis is close to the Earth's axis of rotation, and the OXZ plane is close to the original meridian plane.
- Earth is represented by an ellipsoid, each point M ( ⁇ , ⁇ , h) being identified with respect to this ellipsoid, ⁇ is the longitude, ie the angle between the plane OXZ and the meridian plane containing the point concerned, ⁇ denotes the latitude, the angle between the plane OXY and the normal to the ellipsoid passing through the point concerned and h denotes the ellipsoidal height.
- the ellipsoid representing the Earth is, for example, the ellipsoid WGS84 of half a major axis a and half minor axis b. From the semi-major axis a and the half minor axis b, the eccentricity e and the large normal v ( ⁇ ) of the ellipsoid are defined by the relations:
- FIG. 4 shows an example of steps implemented by the method according to the invention.
- the calculated axial approach MLS configuration described with reference to FIG. 1 and in which the site antenna 7 and the azimuth antenna 8 are located in the same horizontal plane parallel to the plane ⁇ , y) are considered. and located above it at a distance dz.
- a reference point P ref of the same coordinates as those of the azimuth antenna 8 or the site antenna 7 is created.
- a distance p between the reference point P D f and the aircraft 1 1 is determined by a satellite positioning system.
- a third step 43 the real azimuth angle ⁇ R between the plane ⁇ ", z") and the line passing through the azimuth antenna 8 and the aircraft 11 is determined.
- the elevation angle ⁇ between the plane ⁇ , y) and the line passing through the site antenna 7 and the aircraft 11 is determined.
- the position of the aircraft 11 is determined with respect to the reference point P re f from the distance p, the real azimuth angle ⁇ R and the elevation angle ⁇ p .
- Such a method simplifies the determination of the position of the aircraft 1 1 with respect to the reference point P r ef.
- the position of the aircraft 1 1 can be determined via a simplified algorithm.
- Such a simplified algorithm is for example described in standard DO-198 and partially set out in the appendix. It is essentially limited to the determination of roots of a polynomial of the second degree and therefore has several advantages.
- a first advantage is the reduction of the complexity of the calculation of the position of the aircraft 1 1. Therefore, the algorithm is less complex to certify, requires less computing resources than the algorithms of the state of the art and can be executed faster.
- a second benefit is the removal of the iterative character the calculation of the position of the aircraft 1 1. As a result, the risks of multiple solutions, error propagation, divergence and stationary are eliminated. From the position of the aircraft 1 1 with respect to the reference point P ref , the virtual azimuth angle ⁇ v can be determined by the following relation:
- the coordinates y Az and z AZ can be determined from the auxiliary word A1 contained in the guidance signals.
- the distance D may be determined by the auxiliary word A3.
- the step 42 of determining the distance p between the reference point P ref and the aircraft 11 can be carried out according to at least two embodiments.
- FIGS. 5 and 6 A first embodiment is illustrated in FIGS. 5 and 6.
- a distance MLS between the intermediate point A and the aircraft 11 is determined.
- the determination of the distance MLS between the intermediate point A and the aircraft 1 1 comprises the following substeps, illustrated by FIG. 7.
- coordinates ⁇ M , ⁇ M , h M ) of the aircraft 1 1 are determined, for example in the reference R geogra ( ⁇ , ⁇ , h).
- These coordinates ⁇ M , ⁇ M , h M ) are determined by a satellite positioning system, for example the GPS system, the aircraft 1 1 comprising a receiver processing the signals received from the satellites to determine its position.
- the coordinates ⁇ MLS , ⁇ MLS , h MLS ) of the MLS datum point 6 are determined, for example in the reference R geogra ⁇ , ⁇ , h).
- These coordinates ⁇ MLS , ⁇ MLS , h MLS ) can in particular be determined from the signals emitted by the swinging beams 9 and 10.
- these coordinates are contained in the auxiliary words B40. and B41 signals.
- a third step 4223 the coordinates ( ⁇ ⁇ , ⁇ ⁇ , h ⁇ ) of the intermediate point A are determined from the coordinates ⁇ MLS , ⁇ ML s > h ML s) of the ML S datum point 6 by the following relations :
- the distance MLS between the intermediate point A and the aircraft 1 1 is calculated from the coordinates ( ⁇ M , ⁇ M , h M ) of the aircraft 1 1 and coordinates (A ⁇ , ⁇ ⁇ , h ⁇ ) of the intermediate point A.
- the coordinates ( ⁇ M , ⁇ M , h M ) and (A ⁇ , ⁇ ⁇ , h ⁇ ) are converted into coordinates (x M , Y M , Z M ) and (x ⁇ , Y ⁇ , Z ⁇ ) in the reference R geocart by the following formulas:
- d ML s ⁇ IU M - X A Y + (Y M - Y A Y + (Z M - Zj (6)
- the determination of the coordinates of the aircraft 1 1, corresponding to the step 4221 can also be performed after or during the step 4222 and / or step 4223.
- the step 4224 of The calculation of the distance MLS must be performed as soon as possible after the step 4221 for determining the coordinates of the aircraft 11 so as to obtain a distance MLS practically in real time.
- a third step 423a an azimuth angle ⁇ ref between the aircraft 1 1 and a plane parallel to the plane ( ⁇ , z) passing through the reference P ref is determined.
- the reference point P ref has the same coordinates as those of the azimuth antenna 8.
- the azimuth angle ⁇ ref can then be directly determined, this angle being equal to the real azimuth angle ⁇ R.
- a distance dy between the MLS datum point 6 and the reference point P ref along the second axis y is determined.
- the distance dy is determined from the auxiliary words transmitted to the aircraft 1 1 by the signals emitted by the site 4 and azimuth stations 5.
- the distance dy is contained in the bits 21 to Of the auxiliary word A1.
- the distance p is calculated from the azimuth angle ⁇ ref and the distances d MLS and dy by the following relation:
- the coordinates ⁇ AZ , AZ AZ , AZ AZ ) of the azimuth antenna 8, indicated by the point S AZ are determined in FIG. the Cartesian coordinate system R c , for example by the auxiliary words transmitted to the aircraft 1 1 by the site 4 and azimuth stations 5.
- the coordinate x AZ can be determined by the auxiliary word A1, in particular bits 31 to 43
- the coordinate y ⁇ Z can be determined by the bits 21 to 30 of the auxiliary word A1
- the z coordinate AZ can be determined by the bits
- the point R has coordinates in the reference R seocart : (9) where the coordinates ⁇ R , ⁇ R , h R ) can be obtained by auxiliary words
- the vector RS ⁇ Z can be determined in the reference frame R t by the relations:
- angle ⁇ is the angle between the first axis x and the true north, this angle being obtainable by the auxiliary word A1 or B41,
- a third step 423b the coordinates ( ⁇ M , ⁇ M , h M ) of the aircraft 1 1, identified by the point M, are determined, for example in the reference R geogra ⁇ , ⁇ , h).
- These coordinates ( ⁇ M , ⁇ M , h M ) may in particular be determined by a satellite positioning system, for example the GPS, the aircraft 1 1 comprising a receiver processing the signals received from the satellites to determine its position.
- steps 421b and 422b can be performed before, after, or simultaneously with steps 423b and 424b.
- the method described above may in particular be applied by a device on board an aircraft 1 1.
- This device is for example a multimode receiver, commonly called “multimode receiver” and abbreviated "MMR".
- a multimode receiver on board has the function of accurately guiding the aircraft 1 1 during a final approach phase and landing. It may in particular indicate to a pilot of the aircraft 1 1 that it deviates from the optimal trajectory of descent.
- the MMR can normally operate in ILS, VOR, MLS and GPS differential modes, with "VOR” being the acronym for "Very high frequency Omnidirectional Range".
- ILS, VOR, and MLS modes the MMR receives a radio frequency signal from an ILS antenna, a VOR antenna, or MLS antennas, as appropriate.
- the radiofrequency signal is decoded and the deviations of the aircraft 1 1 from the optimum descent trajectory are calculated and transmitted to a flight management system of the aircraft 11.
- the MMR is used in MLS mode.
- the MMR comprises radio frequency signal acquisition means receiving the guidance signals of the swinging beams 9 and 10, radio frequency signal processing means and means for determining the position of the aircraft 1 1 with respect to the reference point P r according to the process described above.
- the MMR may also include a receiver satellite positioning system, for example a GPS receiver.
- the GPS receiver can also be external to the MMR, an interface between the MMR and the GPS receiver being for example made by means of a serial bus, for example an ARINC bus.
- the MMR comprises means for determining the virtual azimuth angle ⁇ v from the position of the aircraft
- the coordinate y M of the aircraft is thus easily calculable from the sole knowledge of p and ⁇ .
- the second solution, x M -x s sin 2 ⁇ + -, corresponds to a
- the first solution, z M -yJp 2 -y M - ⁇ M , corresponds to a position of the aircraft below the azimuth antenna. This solution is not possible.
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Computer Networks & Wireless Communication (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Traffic Control Systems (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Navigation (AREA)
Abstract
L'invention concerne un procédé pour l'aide à l'atterrissage d'aéronef utilisant un GPS et un MLS dans le cadre d'une approche axiale calculée. Selon l'invention, le procédé utilise des coordonnées de l'antenne azimut et/ou de l'antenne site comme point de référence pour le calcul d'une position de l'aéronef dans un repère centré sur la piste d'atterrissage. Cette position de l'aéronef est ensuite utilisée pour déterminer un angle d'azimut entre un axe longitudinal de la piste d'atterrissage et l'aéronef. Le procédé selon l'invention permet de simplifier le calcul de la position de l'aéronef et le calcul de l'angle d'azimut.
Description
Procédé pour l'aide à l'atterrissage d'aéronef utilisant un GPS et un MLS dans le cadre d'une approche axiale calculée
L'invention concerne un procédé pour l'aide à l'approche et à l'atterrissage d'aéronef utilisant un système d'atterrissage hyperfréquence et un système de positionnement par satellites. L'invention concerne plus particulièrement un procédé pour l'aide à l'approche et à l'atterrissage d'aéronef utilisant un système d'atterrissage hyperfréquence dans le cadre d'une approche axiale calculée.
Le système d'atterrissage hyperfréquence, appelé MLS, acronyme de l'expression anglo-saxonne "Microwave Landing System", est un système d'aide à l'approche et à l'atterrissage aux instruments destiné à fournir à un aéronef sa position en coordonnées sphériques dans un repère lié à la piste d'atterrissage, c'est-à-dire un angle d'azimut, un angle de site et une distance entre la piste d'atterrissage et l'aéronef. La distance entre la piste et l'aéronef est fournie par un équipement auxiliaire de mesure de distance, connu sous l'acronyme DME pour l'expression anglo-saxonne "Distance Measurement Equipment" et fonctionnant sur une fréquence différente ou par un système de positionnement par satellites. Le MLS a été développé pour pallier les inconvénients du système d'atterrissage aux instruments, appelé ILS, acronyme de l'expression anglo-saxonne "Instrument Landing System", et retenu par l'organisation de l'aviation civile internationale (OACI) pour succéder à l'ILS. Le MLS permet d'effectuer des approches courbes et segmentées, des atterrissages de catégories I, II et III et une augmentation de cadence d'atterrissage des aéronefs grâce à une diminution de l'espacement entre les aéronefs.
Le MLS, tel que normalisé par l'OACI, transmet des signaux de guidage latéral, c'est-à-dire un angle d'azimut, et de guidage vertical, c'est-à- dire un angle de site, en utilisant une technique de faisceaux battants à temps référencé et un signal multiplexe dans le temps. L'utilisation d'un signal multiplexe dans le temps permet la transmission des signaux de guidage latéral et vertical sur le même canal de radiofréquence sans créer d'interférences entre les signaux de guidage latéral et les signaux de guidage vertical. Les signaux de guidage sont émis sur une fréquence voisine de 5 Giga Hertz (GHz) par une station azimut et une station site. La station azimut est placée en fin de piste tandis que la station site est située sur le côté de la
piste, à environ 300 mètres (m) du seuil de début de piste. Chaque station transmet un faisceau battant étroit balayant en aller et retour à vitesse angulaire constante l'espace de couverture suivant la coordonnée angulaire considérée. Une antenne et un récepteur de bord de l'aéronef reçoivent le faisceau battant une première fois lors du balayage aller et une deuxième fois lors du balayage retour. Il est ainsi possible de déterminer l'angle d'azimut et l'angle de site par la relation linéaire suivante :
Θ = (τ - To ).^ (1 ) où θ est l'angle d'azimut ou l'angle de site, T un intervalle de temps entre la réception des passages aller et retour du faisceau battant,
T0 la valeur de l'intervalle de temps T pour un angle θ nul et v la vitesse angulaire de balayage. T0 et v sont des constantes définies par les normes internationales sur le MLS.
Le système d'atterrissage hyperfréquence à approche axiale calculée, appelé MLS-cc, acronyme de l'expression anglo-saxonne "Microwave Landing System - Computed Centerline", est un MLS dans lequel la station azimut n'est pas placée en bout de piste mais déportée sur un côté de la piste d'atterrissage. Le déplacement de la station azimut peut notamment être utilisé dans deux cas de figure. Dans le premier cas, la station azimut est située à proximité de la station site à des fins de simplicité de déploiement du MLS. Cette configuration est rencontrée principalement dans le cas d'équipements tactiques déployés sur des terrains sommaires et non préparés. Dans le deuxième cas, le MLS-cc est utilisé pour l'approche sur une piste secondaire non équipée d'un MLS mais située dans la zone de couverture de la piste équipée du MLS. Dans les deux cas, du fait du déport de la station azimut, le récepteur de l'aéronef mesure un angle, appelé angle d'azimut réel, qui ne correspond pas à l'angle d'azimut au sens conventionnel du terme, appelé angle d'azimut virtuel. Il est donc nécessaire de calculer l'angle d'azimut virtuel pour pouvoir fournir au pilote une information recentrée par rapport à l'axe de piste. Pour l'évaluer, il est nécessaire de calculer la position de l'aéronef dans un repère cartésien centré sur un point de référence au sol, appelé point au sol. Ce point au sol
est par exemple l'une des deux stations ou le point d'intersection entre l'axe de piste et une droite perpendiculaire à l'axe de piste et passant par la station site. Ce point d'intersection est appelé MLS datum point.
Le calcul de la position de l'aéronef est réalisé par un système de 3 équations à 3 inconnues, paramétré par l'angle d'azimut réel, l'angle de site et une distance entre l'aéronef et le point au sol. Ces équations étant non linéaires, des algorithmes itératifs sont utilisés pour résoudre le système. Classiquement, les algorithmes itératifs sont du type Gauss-Seidel ou Newton-Raphson. En utilisant un système de positionnement par satellites, il est possible d'utiliser le MLS datum point comme point au sol pour la détermination de la distance avec l'aéronef. Cependant, il existe des risques liés aux comportements des algorithmes itératifs dans cette situation, comme exposé dans la norme DO-226. Dans le but d'optimiser la rapidité de convergence vers la position de l'aéronef et la précision de cette position, deux algorithmes itératifs peuvent être utilisés successivement, le premier afin de converger rapidement autour de la position de l'aéronef, le second afin d'obtenir une meilleure précision de cette position. Une telle combinaison d'algorithmes présente cependant l'inconvénient de complexifier la détermination de la position de l'aéronef, rendant difficile l'implantation des algorithmes itératifs et leur validation. Par ailleurs, ces algorithmes sont lents à exécuter, coûteux en termes de charge de calcul et risqués en comportement, en particulier du fait des risques de solutions multiples, de divergence et de stationnante. Enfin, ces algorithmes dégradent la précision de calcul par leur nature itérative, en particulier la propagation d'erreurs.
Les algorithmes itératifs sont exécutés par le récepteur de l'aéronef, par exemple un récepteur multimode, appelé MMR, acronyme de l'expression anglo-saxonne "MultiMode Receiver". Le MMR comprend une chaîne radiofréquence, un processeur de signaux numériques, appelé DSP, acronyme de l'expression anglo-saxonne "Digital Signal Processor", un récepteur de système de positionnement global, mieux connu sous l'acronyme "GPS" pour l'expression anglo-saxonne "Global Positioning System" et un microprocesseur. La chaîne radiofréquence reçoit des signaux provenant de différents systèmes, en particulier les signaux de guidage
provenant de la station azimut et de la station site. Le processeur traite les signaux provenant de la chaîne radiofréquence afin d'extraire les angles d'azimut et de site mais aussi des mots auxiliaires contenus dans les signaux de guidage. Le récepteur GPS peut être interne ou externe. Il transmet les positions de l'avion. Le microprocesseur assure plusieurs fonctions dont un calcul de la distance entre le MLS datum point et l'aéronef, un calcul de la position de l'aéronef à partir des algorithmes itératifs, un calcul de déviations entre la position de l'avion et un axe de descente optimal, une gestion de l'équipement et un interfaçage de communication avec une liaison d'un bus avion, par exemple un bus ARINC, abréviation pour la série de normes "Aeronautical Radio Incorporated". Le microprocesseur comprend un dispositif de gestion des algorithmes permettant de détecter et traiter la divergence, la non-convergence ou la fausse convergence d'un algorithme, de gérer l'initialisation, le séquencement et la combinaison des algorithmes.
Un but de l'invention est notamment de pallier tout ou partie des inconvénients précités. A cet effet, l'invention a pour objet un procédé pour l'aide à l'approche et à l'atterrissage d'aéronef utilisant un système de positionnement par satellites et un système MLS dans le cadre d'une approche axiale calculée, le système MLS comportant une piste d'atterrissage, une station site et une station azimut, la piste d'atterrissage définissant un repère cartésien direct {Rc) comprenant comme axes un axe longitudinal de la piste d'atterrissage, appelé premier axe (x), un axe transversal de la piste d'atterrissage passant par la station site, appelé deuxième axe (y ) et un axe perpendiculaire aux premier et deuxième axes, appelé troisième axe ( z ), un point d'intersection entre le premier axe (x ) et le deuxième axe ( y ) définissant le centre du repère {Rc), appelé point de piste, la station site comprenant une antenne site et la station azimut comprenant une antenne azimut, chaque antenne comprenant un centre de phase, les centres de phase des antennes site et azimut étant situés dans un plan parallèle à un plan {χ,y ) à une distance dz du point de piste selon le troisième axe ( z ), caractérisé en ce qu'il comporte les étapes suivantes :
créer un point de référence de mêmes coordonnées que celles du centre de phase de l'antenne azimut et/ou de l'antenne site, déterminer une distance p entre le point de référence et l'aéronef par le système de positionnement par satellites, déterminer un angle d'azimut ( ΘR ) entre une droite passant par le centre de phase de l'antenne azimut et l'aéronef et un plan parallèle à un plan (χ,z ) passant par le centre de phase de l'antenne azimut, déterminer un angle de site { φ) entre une droite passant par le centre de phase de l'antenne azimut et l'aéronef et le plan parallèle au plan {x,y ) passant les centres de phase des antennes site et azimut, déterminer, à partir de la distance p entre le point de référence et l'aéronef, de l'angle d'azimut et de l'angle de site, la position de l'aéronef par rapport au point de référence.
L'invention a notamment pour avantage qu'elle permet de réduire la complexité du calcul de la position de l'aéronef, le calcul n'utilisant plus d'algorithme du type Gauss-Seidel ou Newton-Raphson mais un algorithme de simple recherche des racines d'un polynôme du second degré, tel que décrit dans la norme DO-198. Par conséquent, l'algorithme est moins complexe à certifier et les problèmes liés à l'utilisation d'algorithmes itératifs sont supprimés. En particulier, les problèmes de divergence et de stationnante, de propagation d'erreurs et de solutions multiples sont supprimés. En outre, l'algorithme est moins coûteux en termes de charge de calcul.
L'invention sera mieux comprise et d'autres avantages apparaîtront à la lecture de la description détaillée de modes de réalisation donnés à titre d'exemple, description faite en regard de dessins annexés qui représentent : - la figure 1 , un exemple de configuration de MLS à approche axiale calculée ;
- la figure 2, un angle de site entre une antenne site et un aéronef ;
- la figure 3, des angles d'azimut entre, d'une part, l'antenne azimut et l'aéronef et, d'autre part, un axe de piste d'atterrissage et l'aéronef ;
- la figure 4, un exemple d'étapes mises en œuvre par le procédé selon l'invention ;
- la figure 5, un exemple de sous-étapes pour la détermination d'une distance entre un point de référence et l'aéronef ; - la figure 6, un exemple de configuration de piste d'atterrissage, de point de référence, de station azimut et d'aéronef pour la détermination de la distance entre le point de référence et l'aéronef ;
- la figure 7, un exemple de sous-étapes pour la détermination d'une distance entre un point intermédiaire et l'aéronef ; - la figure 8, un repère géographique et un repère géocartésien ;
- la figure 9, un autre exemple de sous-étapes pour la détermination de la distance entre le point de référence et l'aéronef.
La figure 1 présente un exemple de configuration de MLS dans le cadre d'une approche axiale calculée. Le MLS comprend une piste d'atterrissage 1 d'axe de piste x orienté d'un seuil de fin de piste 2 vers le seuil de début de piste 3, une station site 4 et une station azimut 5. La piste d'atterrissage 1 définit un repère cartésien direct Rc d'axes l'axe de piste longitudinal, appelé premier axe x , un axe transversal de la piste d'atterrissage 1 passant par la station site 4, appelé deuxième axe y et un axe perpendiculaire aux premier et deuxième axes, appelé troisième axe z . Le deuxième axe y est orienté de la station site 4 vers la piste d'atterrissage 1. Un point d'intersection entre le premier axe x et le deuxième axe y définit le centre du repère Rc et est appelé MLS datum point 6. Il correspond généralement au point de toucher des roues de l'aéronef sur la piste d'atterrissage 1. La station site 4 comprend une antenne site 7 et la station azimut 5 comprend une antenne azimut 8. Chaque antenne comprend un centre de phase à partir duquel les coordonnées de l'antenne peuvent être déterminées. Pour la suite de la description, on assimilera donc la position de chaque antenne à celle de son centre de phase. Les antennes site 7 et azimut 8 émettent des signaux de guidage par des faisceaux battants 9 et 10 pour un guidage vertical et latéral d'un aéronef 1 1 lors d'une phase d'approche et/ou d'atterrissage. Les faisceaux battants 9 et 10 permettent également la transmission de mots auxiliaires contenant, par exemple, des coordonnées de l'antenne site 7 et de l'antenne azimut 8.
La figure 2 représente un angle de site φ entre l'antenne site 7 et l'aéronef 1 1. L'angle de site φ est repéré dans un repère cartésien Rc' centré sur l'antenne site 7 et d'axes un axe x , l'axe y et un axe z , les axes x et z étant parallèles respectivement aux axes x et z . Il est défini par l'angle entre le plan {χ ,y ) et une droite passant par l'antenne site 7 et l'aéronef 1 1.
La figure 3 représente un angle d'azimut réel ΘR et un angle d'azimut virtuel θv . L'angle d'azimut réel ΘR est repéré dans un repère cartésien R" centré sur l'antenne azimut 8 et d'axes un axe x" , l'axe y et un axe z" , les axes x" et z" étant parallèles respectivement aux axes x et z . Il est défini par l'angle entre le plan {χ",z") et une droite passant par l'antenne azimut 8 et l'aéronef 11. L'angle d'azimut virtuel θv est repéré dans le repère cartésien Rc et défini par l'angle entre le plan (x,z ) et une droite passant par l'aéronef 1 1 et un point virtuel 8v du plan (χ,z ). Ce point virtuel correspond à une antenne azimut dans le cadre d'un système MLS, habituellement située sur le premier axe x à une distance de 300 m après le seuil de fin de piste 2.
Pour la suite de la description, on considère en plus du repère cartésien Rc centré sur le MLS datum point 6, du repère cartésien Rc' centré sur l'antenne site 7 et du repère cartésien R" centré sur l'antenne azimut 8, un repère géocartésien Rgeocart{θ,X,Y,z) et un repère géographique R geogÀλ >Φ> h), illustrés à la figure 8.
Le repère géocartésien Rgeocart(θ,X,Y,z) est un repère dans lequel le centre O est proche du centre de masse de la Terre, l'axe OZ est proche de l'axe de rotation terrestre et le plan OXZ est proche du plan méridien origine.
Le repère géographique Rgeogra(λ,φ,h) est un repère dans lequel la
Terre est représentée par un ellipsoïde, chaque point M(λ,φ,h) étant repéré par rapport à cet ellipsoïde, λ désigne la longitude, soit l'angle entre le plan OXZ et le plan méridien contenant le point concerné, φ désigne la latitude, soit l'angle entre le plan OXY et la normale à l'ellipsoïde passant par le point concerné et h désigne la hauteur ellipsoïdale. Chaque point M(λ,φ,h) définit un vecteur Ux parallèle au plan P tangent à l'ellipsoïde passant par le projeté orthogonal de M(λ,φ,h) et orienté vers le nord vrai, un vecteur ûφ
parallèle au plan P et orienté vers l'Est et un vecteur ûh orthogonal au plan
P . L'ellipsoïde représentant la Terre est par exemple l'ellipsoïde WGS84 de demi grand axe a et de demi petit axe b . A partir du demi grand axe a et du demi petit axe b , on définit l'excentricité e et la grande normale v(φ) de l'ellipsoïde par les relations :
La figure 4 montre un exemple d'étapes mises en œuvre par le procédé selon l'invention. Pour ce procédé, on considère la configuration de MLS à approche axiale calculée décrite en référence à la figure 1 et dans laquelle l'antenne site 7 et l'antenne azimut 8 sont situées dans un même plan horizontal parallèle au plan {χ,y ) et situé au-dessus de celui-ci à une distance dz . Dans une première étape 41 , un point de référence Préf de mêmes coordonnées que celles de l'antenne azimut 8 ou de l'antenne site 7 est créé. Dans une deuxième étape 42, une distance p entre le point de référence Préf et l'aéronef 1 1 est déterminée par un système de positionnement par satellites. Dans une troisième étape 43, l'angle d'azimut réel ΘR entre le plan {χ",z") et la droite passant par l'antenne azimut 8 et l'aéronef 11 est déterminé. Dans une quatrième étape 44, l'angle de site <p entre le plan {χ ,y) et la droite passant par l'antenne site 7 et l'aéronef 1 1 est déterminée. Dans une cinquième étape 45, la position de l'aéronef 11 est déterminée par rapport au point de référence Préf à partir de la distance p , de l'angle d'azimut réel ΘR et de l'angle de site <p .
Un tel procédé permet de simplifier la détermination de la position de l'aéronef 1 1 par rapport au point de référence Préf. En effet, la position de l'aéronef 1 1 peut être déterminée par l'intermédiaire d'un algorithme simplifié. Un tel algorithme simplifié est par exemple décrit dans la norme DO-198 et exposé partiellement en annexe. Il se limite essentiellement à la détermination de racines d'un polynôme du second degré et présente par conséquent plusieurs avantages. Un premier avantage est la réduction de la complexité du calcul de la position de l'aéronef 1 1. Par conséquent, l'algorithme est moins complexe à certifier, demande moins de ressources de calcul que les algorithmes de l'état de l'art et peut être exécuté plus rapidement. Un deuxième avantage est la suppression du caractère itératif
du calcul de la position de l'aéronef 1 1 . Par conséquent, les risques de solutions multiples, de propagation d'erreurs, de divergence et de stationnante sont supprimés. A partir de la position de l'aéronef 1 1 par rapport au point de référence Préf, l'angle d'azimut virtuel θv peut être déterminé par la relation suivante :
Où {xM , yM , zM ) sont les coordonnées de l'aéronef 1 1 dans le repère Rc ,
(xAZ , y Az , zAZ ) sont les coordonnées de l'antenne azimut 8 dans le repère Rc , Dse est la distance entre le seuil de fin de piste 2 et le MLS datum point 6.
Les coordonnées y Az et zAZ peuvent être déterminées à partir du mot auxiliaire A1 contenu dans les signaux de guidage. La distance Dse peut être déterminée par le mot auxiliaire A3.
L'étape 42 de détermination de la distance p entre le point de référence Préf et l'aéronef 1 1 peut être réalisée selon au moins deux modes de réalisation.
Un premier mode de réalisation est illustré par les figures 5 et 6.
Dans une première étape 421 a, un point intermédiaire A est créé. Le point intermédiaire A a, dans le repère Rc , les mêmes coordonnées que le MLS datum point 6 suivant le premier axe x et le deuxième axe y et la même coordonnée que le point de référence Préf suivant le troisième axe z . Le point intermédiaire A est donc situé au-dessus du MLS datum point 6, la distance dz séparant ces deux points.
Dans une deuxième étape 422a, une distance dMLS entre le point intermédiaire A et l'aéronef 1 1 est déterminée. Selon un mode particulier de réalisation, la détermination de la distance dMLS entre le point intermédiaire A et l'aéronef 1 1 comprend les sous-étapes suivantes, illustrées par la figure 7. Dans une première étape 4221 , des coordonnées {λM ,φM ,hM ) de l'aéronef 1 1 sont déterminées, par exemple dans le repère Rgeogra (Λ,φ,h). Ces coordonnées {λM ,φM ,hM ) sont déterminées par un système de positionnement par satellites, par exemple le système GPS, l'aéronef 1 1
comprenant un récepteur traitant les signaux reçus des satellites pour déterminer sa position. Dans une deuxième étape 4222, les coordonnées {λMLS,φMLS,hMLS ) du MLS datum point 6 sont déterminées, par exemple dans le repère Rgeogra{λ,φ,h). Ces coordonnées {λMLS,φMLS,hMLS ) peuvent notamment être déterminées à partir des signaux émis par les faisceaux battants 9 et 10. En particulier, selon les normes définies par l'OACI, ces coordonnées sont contenus dans les mots auxilaires B40 et B41 des signaux. Dans une troisième étape 4223, les coordonnées (λΛ,φΛ,hΛ ) du point intermédiaire A sont déterminées à partir des coordonnées {ΛMLS,ΦMLs> h MLs ) du MLS datum point 6 par les relations suivantes :
Dans une quatrième étape 4224, la distance dMLS entre le point intermédiaire A et l'aéronef 1 1 est calculée à partir des coordonnées (λM ,φM ,hM ) de l'aéronef 1 1 et des coordonnées (AΛ,φΛ,hΛ ) du point intermédiaire A. Dans une forme particulière de réalisation, les coordonnées (λM,φM,hM ) et (AΛ,φΛ,hΛ ) sont converties en coordonnées (xM , YM , ZM ) et (xΛ,YΛ,ZΛ ) dans le repère Rgeocart par les formules suivantes :
XM = (V(ΦM ) + K )C0S ΦM C0S λM γM = (V(ΦM ) + K ) COS φM sin λM (5)
ZM = ïb 2/a 2 HΦM
h M\ sin ΦM Les mêmes formules permettent de déterminer les coordonnées (xΛ,YΛ,ZΛ ) en remplaçant λM , φM e\ hM par λA , φA et hA .
La distance dMLS peut ensuite être calculée par la relation suivante : dMLs = ΛIUM - XA Y + (YM - YA Y + (ZM - Zj (6)
A noter que la détermination des coordonnées de l'aéronef 1 1 , correspondant à l'étape 4221 , peut également être effectuée après ou pendant l'étape 4222 et/ou l'étape 4223. Dans tous les cas, l'étape 4224 de calcul de la distance dMLS doit être accomplie le plus rapidement possible après l'étape 4221 de détermination des coordonnées de l'aéronef 1 1 afin d'obtenir une distance dMLS pratiquement en temps réel.
Dans une troisième étape 423a, un angle d'azimut θref entre l'aéronef 1 1 et un plan parallèle au plan (χ,z ) passant par le point de
référence Préf est déterminé. Avantageusement, le point de référence Préf a les mêmes coordonnées que celles de l'antenne azimut 8. L'angle d'azimut θref peut alors être directement déterminé, cet angle étant égal à l'angle d'azimut réel ΘR . Dans une quatrième étape 424a, une distance dy entre le MLS datum point 6 et le point de référence Préf suivant le deuxième axe y est déterminée. Selon un mode particulier de réalisation, la distance dy est déterminée à partir des mots auxiliaires transmis à l'aéronef 1 1 par les signaux émis par les stations site 4 et azimut 5. En particulier, la distance dy est contenue dans les bits 21 à 30 du mot auxiliaire A1.
Dans une cinquième étape 425a, la distance p est calculée à partir de l'angle d'azimut θref et des distances dMLS et dy par la relation suivante :
P = ^MLS 2 - {dy- C0S θref Ϋ + <ty- sin θ ref - (7) A noter que l'ordre des étapes 422a, 423a et 424a n'a pas d'importance, les opérations liées à ces étapes pouvant être effectuées dans un ordre différent ou simultanément.
Un deuxième mode de réalisation est illustré par les figures 8 et 9. Dans une première étape 421 b, les coordonnées {χAZ ,yAZ ,zAZ ) de l'antenne azimut 8, repérée par le point SAZ , sont déterminées dans le repère cartésien Rc , par exemple par les mots auxiliaires transmis à l'aéronef 1 1 par les stations site 4 et azimut 5. Selon les normes établies par l'OACI, la coordonnée xAZ peut être déterminée par le mot auxiliaire A1 , en particulier les bits 31 à 43, la coordonnée yΛZ peut être déterminée par les bits 21 à 30 du mot auxiliaire A1 et la coordonnée zAZ peut être déterminée par les bits
57 à 63 du mot auxiliaire A1 .
Dans une deuxième étape 422b, les coordonnées {χAZ ,yAZ ,zAZ ) dans le repère Rc sont converties en coordonnées (xΛZ , YΛZ , ZΛZ ) dans le repère Rgeocart . La conversion peut notamment être effectuée en décomposant le vecteur OSΛZ en un vecteur OR et un vecteur RSΛZ où le point R est le MLS datum point 6 :
~δsZ = ÔR + RS^ (8)
Le point R a pour coordonnées dans le repère Rseocart :
(9)
où les coordonnées {λR,φR,hR ) peuvent être obtenus par des mots auxiliaires
B40 et B41 provenant des faisceaux battants 9 et 10.
Le vecteur RSΛZ peut être déterminé dans le repère R t par les relations :
Où l'angle ψ est l'angle entre le premier axe x et le nord vrai, cet angle pouvant être obtenu par le mot auxiliaire A1 ou B41 ,
Et où les vecteurs ûλ , ûφ et ûh peuvent être déterminés dans le repère
Kocan Par les relations :
'— sin λ,
Dans une troisième étape 423b, les coordonnées (λM ,φM ,hM ) de l'aéronef 1 1 , repéré par le point M , sont déterminées, par exemple dans le repère Rgeogra{λ,φ,h). Ces coordonnées (λM ,φM ,hM ) peuvent notamment être déterminées par un système de positionnement par satellites, par exemple le GPS, l'aéronef 1 1 comprenant un récepteur traitant les signaux reçus des satellites pour déterminer sa position.
Dans une quatrième étape 424b, les coordonnées (AM,φM,hM ) dans le repère Rgeogra{λ,φ,h) sont converties en coordonnées (xM, YM , ZM ) dans le repère Rgeocart - L'aéronef 1 1 a pour coordonnées dans le repère
χ M = {V{ΦM ) + K )cos φM cos λM γM = HΦM ) + hM ) cos φM ήnλM (12)
Dans une cinquième étape 425b, la distance p est calculée à partir des coordonnées {xAZ , YAZ ,ZAZ ) et {xM, YM , ZM ) par la relation suivante : p = Vfc - x ΛZ y + (γ M - Y ΛZ y + (z M - zΛZ y (13)
A noter que les étapes 421 b et 422b peuvent être effectuées avant, après ou simultanément aux étapes 423b et 424b.
Ce deuxième mode de réalisation présente l'avantage de ne pas utiliser de mesure d'angle d'azimut pour le calcul de la distance p entre le point de référence Prét et l'aéronef 11. Il en résulte une plus grande précision du calcul de la distance p et par conséquent du calcul de la position de l'aéronef 1 1.
Le procédé décrit ci-dessus peut notamment être appliqué par un dispositif embarqué à bord d'un aéronef 1 1. Ce dispositif est par exemple un récepteur multimode de bord, communément appelé "multimode receiver" et abrégé par "MMR". Un récepteur multimode de bord a pour fonction de guider avec précision l'aéronef 1 1 pendant une phase d'approche finale et d'atterrissage. Il peut notamment indiquer à un pilote de l'aéronef 1 1 qu'il s'écarte de la trajectoire optimale de descente. Le MMR peut normalement opérer en modes ILS, VOR, MLS et GPS différentiel, "VOR" étant l'acronyme de "Very high frequency Omnidirectional Range". En modes ILS, VOR et MLS, le MMR reçoit un signal radiofréquence d'une antenne ILS, d'une antenne VOR ou des antennes MLS, selon le cas. Le signal radiofréquence est décodé et les déviations de l'aéronef 1 1 par rapport à la trajectoire optimale de descente sont calculées et transmises à un système de gestion de vol de l'aéronef 11.
Selon l'invention, le MMR est utilisé en mode MLS. Le MMR comprend des moyens d'acquisition de signaux radiofréquences recevant les signaux de guidage des faisceaux battants 9 et 10, des moyens de traitement des signaux radiofréquences et des moyens pour déterminer la position de l'aéronef 1 1 par rapport au point de référence Prét selon le procédé décrit ci-dessus. Le MMR peut également comprendre un récepteur
de système de positionnement par satellites, par exemple un récepteur GPS. Selon une autre forme de réalisation, le récepteur GPS peut également être externe au MMR, une interface entre le MMR et le récepteur GPS étant par exemple réalisée au moyen d'un bus série, par exemple un bus ARINC. Dans une forme particulière de réalisation, le MMR comprend des moyens pour déterminer l'angle d'azimut virtuel θv à partir de la position de l'aéronef
11. Le bus ARINC peut être utilisé pour échanger des données entre le MMR et un système de gestion de vol de l'aéronef 1 1. En particulier, la position de l'aéronef 1 1 , l'angle d'azimut virtuel θv , l'angle de site φ et/ou la distance p peuvent être transmis à un pilote automatique de l'aéronef 1 1.
ANNEXE
Cette annexe présente la problématique de calcul de la position d'un aéronef et un algorithme simplifié permettant ce calcul inspiré de la norme DO-198.
Présentation du problème et formalisation géométrique
La fonction MLS-cc est utilisée pour des atterrissages sur des pistes équipées de stations MLS situées dans des positions non optimales (déport de la station azimut ou station azimut d'une autre piste d'atterrissage).
La connaissance de l'angle d'azimut et de l'angle de site obtenus à partir des informations transmises par les stations MLS mais aussi de la distance par rapport à un point au sol (utilisation d'un DME ou d'un GPS) permet de calculer la position de l'aéronef, la position des stations les unes par rapport aux autres étant connue.
On définit les notations suivantes dans un repère cartésien direct :
{ xM ,yM , zM ) : coordonnées de l'aéronef, - ( xΛZ , yAZ , zΛZ ) : coordonnées de l'antenne azimut,
- ( xs, ys , zs ) : coordonnées de l'antenne site,
- ( xD,yD ,zD ) : coordonnées du DME ou du point de référence
Préf, p : distance entre l'aéronef et le DME ou le point de référence, - θ : angle d'azimut, en radians, φ : angle de site, en radians. L'angle φ correspond à l'angle entre l'aéronef et l'axe de piste et non à l'angle entre l'aéronef et le plan de descente tel que mesuré en MLS. Une correction, consistant à ajouter la pente, est donc nécessaire pour obtenir φ à partir de la donnée fournie par le MLS.
Il est possible de calculer p , θ et φ à partir des relations suivantes : p = \(χ M - χ D y + (y - y D y + (z - ZD y (0.1 )
θ = (0.2)
Les équations ci-dessus sont équivalentes à :
UM = zs + tan XM - Xs )2 + (y M - y s )2 (1 .3)
- La résolution du système d'équations non linéaires (1.1 ) à (1.3) est complexe et nécessite l'utilisation d'un algorithme itératif, par exemple défini dans la norme DO-226. Toutefois, dans certains cas de figure, le problème peut être simplifié. Pour cela, il faut des caractéristiques précises de géométrie au sol. Selon ces caractéristiques, l'une des deux antennes des stations MLS et le point de référence sont situés au même endroit et les antennes site et azimut sont situées dans un même plan horizontal.
Simplification du problème et hypothèses nécessaires
L'origine du repère cartésien est fixée au point où se trouve l'antenne azimut. On a donc {xAZ,yAZ ,zAZ ) = {θ,O,θ) . La distance entre l'aéronef et la station azimut est déterminée par les coordonnées {xM,yM,zM ) de l'aéronef obtenues par un récepteur de système de positionnement par satellites à bord de l'aéronef et les coordonnées {χAZ,yAZ,zAZ ) de l'antenne azimut transmises par la station site ou la station azimut dans les mots auxiliaires B40 et B41.
Pour la suite de la description, on considérera que l'antenne azimut et le point de référence sont situés au même endroit. On obtiendrait cependant des résultats similaires en considérant que l'antenne site et le point de référence sont situés au même endroit. On a donc (xAZ,yAZ ,zAZ ) = (xD,yD,zD) . A partir des équations (1.1 ) à (1.3), on peut écrire :
2 2 2 2
XM ~ P ~ y M ~ ZM (2.1) y M =-tanθy]xM 2 + Z M (2.2) zM =zs+ tan φ^{xM - xs f + {yM - ys f (2.3)
La combinaison de (2.1 ) et (2.2) donne : yM =pύnθ (2.4)
La coordonnée yM de l'aéronef est donc facilement calculable à partir de la seule connaissance de p et θ.
Du fait que les antennes site et azimut sont situées dans le même plan horizontal, on a zs=0. La combinaison des équations (2.1) et (2.3) donne : xM 2 -2xMxssm2 φ + sm2 φ\d2 + yM 2 -2yMys)- p2 cos2 θ COS2^ = O
(2.5) où d désigne la distance entre l'antenne site et l'antenne azimut, soit : d =^xs 2+ ys 2
Trouver la coordonnée xM de l'aéronef se réduit donc à résoudre une équation du second degré (2.5) où a = \, b = -2xssm2 φ et c = sin2 φ(d2 +yM 2 -2yMys)- p2 cos2 θ cos2 φ. Le déterminant est :
Δ = 4x2sin4^-4|sin2(%/2 + yM 2 -2yMys)- p2 cos2 θ cos2 φ\
La résolution de cette équation aboutit à des solutions multiples puisque les solutions sont de la forme : xM = -xs sin2 φ± (2.6)
La première solution, xM =-xsήn2 φ , est toujours négative, la coordonnée xs étant toujours positive. Cette solution correspond à une situation où l'aéronef est situé derrière l'antenne azimut, c'est-à-dire en dehors du secteur couvert par le balayage des stations site et azimut. Cette solution n'est donc pas possible. La seconde solution, xM =-xssin2 φ+ — , correspond à une
situation où l'aéronef est situé devant l'antenne azimut. C'est la solution à considérer dans le cas d'un aéronef effectuant une approche axiale calculée.
La coordonnée zM de l'aéronef peut alors être déterminée par l'équation (2.1 ), soit :
La première solution, zM = -yJp2 -yM -χM , correspond à une position de l'aéronef en dessous de l'antenne azimut. Cette solution n'est donc pas possible.
La deuxième solution, zM = +-Jp2 - yM -χM , correspond à une position de l'aéronef au dessus de l'antenne azimut. C'est la solution à considérer. II est donc toujours possible de déterminer la position de l'aéronef par les relations suivantes :
Claims
1. Procédé pour l'aide à l'approche et à l'atterrissage d'aéronef (1 1 ) utilisant un système de positionnement par satellites et un système MLS dans le cadre d'une approche axiale calculée, le système MLS comportant une piste d'atterrissage (1 ), une station site (4) et une station azimut (5), la piste d'atterrissage (1 ) définissant un repère cartésien direct
(Rc ) comprenant comme axes un axe longitudinal de la piste d'atterrissage
(1 ), appelé premier axe (x ), un axe transversal de la piste d'atterrissage (1 ) passant par la station site (4), appelé deuxième axe (y ) et un axe perpendiculaire aux premier et deuxième axes, appelé troisième axe ( z ), un point d'intersection entre le premier axe (x ) et le deuxième axe ( y ) définissant le centre du repère {Rc ), appelé point de piste (MLS datum point
6), la station site (4) comprenant une antenne site (7) et la station azimut (5) comprenant une antenne azimut (8), chaque antenne comprenant un centre de phase, les centres de phase des antennes site (7) et azimut (8) étant situés dans un plan parallèle à un plan {χ,y ) à une distance dz du point de piste (MLS datum point 6) selon le troisième axe ( z ), caractérisé en ce qu'il comporte les étapes suivantes :
(41 ) créer un point de référence (Préf) de mêmes coordonnées que celles du centre de phase de l'antenne azimut (8) et/ou de l'antenne site (7),
(42) déterminer une distance p entre le point de référence (Préf) et l'aéronef (1 1 ) par le système de positionnement par satellites,
(43) déterminer un angle d'azimut { ΘR ) entre une droite passant par le centre de phase de l'antenne azimut (8) et l'aéronef (1 1 ) et un plan parallèle à un plan (x,z ) passant par le centre de phase de l'antenne azimut
(8),
(44) déterminer un angle de site ( φ) entre une droite passant par le centre de phase de l'antenne azimut (8) et l'aéronef (1 1 ) et le plan parallèle au plan {χ,y) passant les centres de phase des antennes site (7) et azimut (8),
(45) déterminer, à partir de la distance p entre le point de référence (Préf) et l'aéronef (1 1 ), de l'angle d'azimut et de l'angle de site, la position de l'aéronef (11 ) par rapport au point de référence (Préf).
2. Procédé selon la revendication 1 , caractérisé en ce que la détermination de la distance p entre le point de référence (Préf) et l'aéronef
(1 1 ) comprend les étapes suivantes : créer, dans le repère Rc , un point intermédiaire (A) de mêmes coordonnées {λΛ ,φΛ,hΛ ) que les coordonnées {λMLS , φMLS , hMLS ) du point de piste (MLS datum point 6) suivant le premier axe ( x ) et le deuxième axe ( y ) et de même cordonnée que le point de référence (Préf) suivant le troisième axe ( z ), déterminer une distance dMLS entre le point intermédiaire (A) et l'aéronef (1 1 ), déterminer un angle d'azimut θref entre le plan ( x, z ) et une droite passant par le point de référence (Préf) et l'aéronef, déterminer une distance dy entre le point de piste (MLS datum point 6) et le point de référence (Préf) suivant le deuxième axe { y ), calculer la distance p entre le point de référence (Prêt) et l'aéronef
(1 1 ) par la relation p = ^dMLS 2 - {dy. cos θref f + dy. sin θref
3. Procédé selon la revendication 2, caractérisé en ce que la détermination de la distance dMLS entre le point intermédiaire (A) et l'aéronef (1 1 ) comprend les étapes suivantes : déterminer des coordonnées (λM,φM,hM ) de l'aéronef (1 1 ), déterminer des coordonnées {λMLS , φMLS , hMLS ) du point de piste
(MLS datum point 6), déterminer des coordonnées {λΛ,φΛ,hΛ ) du point intermédiaire (A) à partir des coordonnées {λMLS,φMLS,hMLS ) du point de piste (MLS datum point
6), calculer la distance dMLS entre le point intermédiaire (A) et l'aéronef (1 1 ) à partir des coordonnées {λM ,φM ,hM ) de l'aéronef (1 1 ) et celles (ΛΛ,φΛ,hΛ ) du point intermédiaire (A).
4. Procédé selon la revendication 3, caractérisé en ce que les coordonnées {λMLS , φMLS , hMLS ) du point de piste (MLS datum point 6) sont déterminées à partir de signaux émis par la station site (4) ou la station azimut (5).
5. Procédé selon la revendication 1 , caractérisé en ce que la détermination de la distance p entre le point de référence (Préf) et l'aéronef
(1 1 ) comprend les étapes suivantes : déterminer, dans le premier repère Rc , les coordonnées {χAZ,yAZ,zAZ ) du point de référence (Préf), convertir les coordonnées {χΛZ,yΛZ,zΛZ ) du point de référence (Prêt) du premier repère Rc dans un deuxième repère, appelé repère géocartésien Rgeocart , déterminer, dans un troisième repère, appelé repère géographique R geogra > des coordonnées [λM,φM,hM ) de l'aéronef (1 1 ), convertir les coordonnées {λM ,φM ,hM ) de l'aéronef (1 1 ) du troisième repère Rgeogra dans le deuxième repère Rgeocart , calculer la distance p entre le point de référence (Prêt) et l'aéronef (1 1 ) par les coordonnées {{xAZ, YAZ,ZAZ ) et {xM ,YM,ZM )) du point de référence (Prêt) et de l'aéronef (1 1 ) dans le deuxième repère Rgeocart .
6. Procédé selon la revendication 5, caractérisé en ce que les coordonnées {χΛZ,yΛZ,zΛZ ) du point de référence (Préf) dans le premier repère Rc sont déterminées par des données provenant de signaux émis par le MLS.
7. Procédé selon l'une des revendications 3 à 6, caractérisé en ce que les coordonnées (ΛM ,ψM ,hM ) de l'aéronef (1 1 ) sont déterminées par un système de positionnement global (GPS).
8. Récepteur multimode pouvant équiper un aéronef (1 1 ) pour l'aide à l'approche et à l'atterrissage à approche axiale calculée comportant des moyens d'acquisition de signaux radiofréquences et des moyens de traitement des signaux radiofréquences, caractérisé en ce qu'il comprend, en outre, des moyens pour déterminer une position d'un aéronef (1 1 ) par le procédé selon l'une des revendications précédentes.
9. Récepteur multimode selon la revendication 8, caractérisé en ce que le récepteur multimode comprend un récepteur de système de positionnement par satellites.
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| AT08866617T ATE554432T1 (de) | 2007-12-21 | 2008-12-19 | Verfahren zur flugzeuglandeunterstützung mithilfe von gps und mls im falle einer berechneten axial- annäherung |
| ES08866617T ES2385421T3 (es) | 2007-12-21 | 2008-12-19 | Procedimiento de ayuda al aterrizaje de una aeronave utilizando GPS y MLS en el marco de una aproximación axial calculada |
| EP08866617A EP2225622B1 (fr) | 2007-12-21 | 2008-12-19 | Procede pour l'aide a l'atterrissage d'aeronef utilisant un gps et un mls dans le cadre d'une approche axiale calculee |
| US12/809,636 US8290644B2 (en) | 2007-12-21 | 2008-12-19 | Method for aircraft landing assistance using GPS and MLS in case of calculated axial approach |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR07/09035 | 2007-12-21 | ||
| FR0709035A FR2925712B1 (fr) | 2007-12-21 | 2007-12-21 | Procede pour l'aide a l'atterrissage d'aeronef utilisant un gps et un mls dans le cadre d'une approche axiale calculee. |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2009083507A1 true WO2009083507A1 (fr) | 2009-07-09 |
Family
ID=39639341
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/EP2008/068084 Ceased WO2009083507A1 (fr) | 2007-12-21 | 2008-12-19 | Procede pour l'aide a l'atterrissage d'aeronef utilisant un gps et un mls dans le cadre d'une approche axiale calculee |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US8290644B2 (fr) |
| EP (1) | EP2225622B1 (fr) |
| AT (1) | ATE554432T1 (fr) |
| ES (1) | ES2385421T3 (fr) |
| FR (1) | FR2925712B1 (fr) |
| WO (1) | WO2009083507A1 (fr) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2946447B1 (fr) | 2009-06-05 | 2011-07-22 | Thales Sa | Procede pour permettre un atterrissage sur piste decalee |
| US9934692B2 (en) | 2012-02-22 | 2018-04-03 | Honeywell International Inc. | Display system and method for generating a display |
| US8949027B2 (en) * | 2012-07-10 | 2015-02-03 | Honeywell International Inc. | Multiple truth reference system and method |
| CN103646384B (zh) * | 2013-12-20 | 2016-06-22 | 江苏大学 | 一种遥感扫描成像平台飞行速度的优化方法 |
| US20170206663A1 (en) * | 2014-07-17 | 2017-07-20 | Nec Corporation | Airspace information processing device, airspace information processing method, and non-transitory computer-readable medium storing airspace information processing program |
| CN105571593B (zh) * | 2015-11-27 | 2018-08-31 | 中国电子科技集团公司第二十研究所 | 一种基于mls的地理位置信息获取方法 |
| US10152195B2 (en) | 2015-12-14 | 2018-12-11 | Honeywell International Inc. | Aircraft display system pertaining to energy management |
| US10654589B2 (en) * | 2017-03-27 | 2020-05-19 | Honeywell International Inc. | Avionic display systems and methods for generating vertical situation displays including instability prediction and avoidance symbology |
| US10839701B2 (en) | 2018-06-05 | 2020-11-17 | Honeywell International Inc. | Methods and systems for stabilized approach energy management |
| US10854091B2 (en) | 2018-07-03 | 2020-12-01 | Honeywell International Inc. | Energy management visualization methods and systems |
| US11275171B2 (en) * | 2018-11-26 | 2022-03-15 | Sr Technologies, Inc. | Location using round trip time vectors using an airborne platform |
| FR3103178B1 (fr) * | 2019-11-20 | 2022-12-30 | Thales Sa | Procede d'aide a l'atterrissage d'un aeronef sur une piste d'atterrissage |
| US11842629B2 (en) * | 2020-12-10 | 2023-12-12 | Honeywell International Inc. | Dynamic radar vectoring guidance methods and systems |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH02287900A (ja) * | 1989-04-28 | 1990-11-27 | Japan Radio Co Ltd | Mls機上装置 |
| EP0793116A2 (fr) * | 1996-02-28 | 1997-09-03 | Nec Corporation | Appareil pour détecter le passage d'un objet volant par un espace surveillé |
| EP1589351A1 (fr) * | 2004-04-22 | 2005-10-26 | Airbus France | Procédé et dispositif d'aide à l'atterrissage d'un aéronef sur une piste d'atterrissage |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4454510A (en) * | 1978-12-18 | 1984-06-12 | Crow Robert P | Discrete address beacon, navigation and landing system (DABNLS) |
| US6239745B1 (en) * | 1999-07-30 | 2001-05-29 | Rockwell Collins, Inc. | Satellite landing system having instrument landing system look alike guidance |
-
2007
- 2007-12-21 FR FR0709035A patent/FR2925712B1/fr not_active Expired - Fee Related
-
2008
- 2008-12-19 WO PCT/EP2008/068084 patent/WO2009083507A1/fr not_active Ceased
- 2008-12-19 EP EP08866617A patent/EP2225622B1/fr active Active
- 2008-12-19 ES ES08866617T patent/ES2385421T3/es active Active
- 2008-12-19 US US12/809,636 patent/US8290644B2/en not_active Expired - Fee Related
- 2008-12-19 AT AT08866617T patent/ATE554432T1/de active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH02287900A (ja) * | 1989-04-28 | 1990-11-27 | Japan Radio Co Ltd | Mls機上装置 |
| EP0793116A2 (fr) * | 1996-02-28 | 1997-09-03 | Nec Corporation | Appareil pour détecter le passage d'un objet volant par un espace surveillé |
| EP1589351A1 (fr) * | 2004-04-22 | 2005-10-26 | Airbus France | Procédé et dispositif d'aide à l'atterrissage d'un aéronef sur une piste d'atterrissage |
Non-Patent Citations (3)
| Title |
|---|
| D. VAN WILLIGEN ET AL.: "Radio Navigation - Air traffic systems", 26 September 2006 (2006-09-26), Faculty Electrical Engineering, Mathematics and Computer Science, Delft University of Technology, pages 51 - 64, XP002490240, Retrieved from the Internet <URL:http://www.pn.ewi.tudelft.nl/education/et4-022/notes/> [retrieved on 20080729] * |
| HALL, JOHN W., HATZIS, PATRICIA M., POWELL, FREDERIC D.: "Examination of RTCA/DO-198 Position Reconstruction Algorithms for Area Navigation with the Microwave Landing System", DEFENSE TECHNICAL INFORMATION CENTER OAI-PMH REPOSITORY, June 1990 (1990-06-01), US, pages 1 - 47, XP009104064 * |
| POWELL F D: "A rotated-coordinate Gaussian position reconstruction algorithm for the microwave landing system", NAVIGATION. JOURNAL OF THE INSTITUTE OF NAVIGATION USA, vol. 37, no. 3, November 1990 (1990-11-01), USA, pages 273 - 283, XP009104065, ISSN: 0028-1522 * |
Also Published As
| Publication number | Publication date |
|---|---|
| EP2225622A1 (fr) | 2010-09-08 |
| US20100286852A1 (en) | 2010-11-11 |
| FR2925712B1 (fr) | 2010-01-01 |
| ATE554432T1 (de) | 2012-05-15 |
| FR2925712A1 (fr) | 2009-06-26 |
| EP2225622B1 (fr) | 2012-04-18 |
| US8290644B2 (en) | 2012-10-16 |
| ES2385421T3 (es) | 2012-07-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2225622B1 (fr) | Procede pour l'aide a l'atterrissage d'aeronef utilisant un gps et un mls dans le cadre d'une approche axiale calculee | |
| EP0925515B1 (fr) | Procede d'aide a la navigation d'un mobile vers une cible egalement mobile | |
| EP2648017B1 (fr) | Système embarqué d'aide au pilotage d'un aéronef, basé sur un système GNSS, à architecture redondante et dissimilaire pour niveau d'integrité élevé | |
| EP0327447B1 (fr) | Système de radionavigation | |
| EP0875002B1 (fr) | Systeme d'aide au pilotage d'aeronefs a l'aide d'un viseur tete haute | |
| EP2353024B1 (fr) | Procede de geo-localisation d'un objet par multitelemetrie | |
| EP2902797B1 (fr) | Procédé et système radiofréquence de détermination, par couple d'engins spatiaux, de la position angulaire relative entre plusieurs engins spatiaux distants | |
| EP2490042B1 (fr) | Procédé et système de détermination des paramètres de navigation d'un aéronef | |
| WO2007082929A1 (fr) | Dispositif de controle de position(s) relative(s) par analyse de signaux bi-frequences, pour un engin spatial d'un groupe d'engins spatiaux en formation | |
| EP3698173B1 (fr) | Système de positionnement avec moyens de génération de signaux gnss et câble rayonnant | |
| EP2374022B1 (fr) | Dispositif d'hybridation en boucle fermee integre par construction | |
| WO2008040896A2 (fr) | Procédé de localisation d'un véhicule | |
| EP1460504B1 (fr) | Système de pilotage d'un aéronef, au moins pour piloter l'aéronef lors d'une approche de non précision en vue d'un atterrissage | |
| EP1205732B1 (fr) | Centrale inertielle de navigation comportant un récepteur GPS intégré | |
| EP2765390A1 (fr) | Système et procédé d'aide à la navigation d'un aéronef | |
| WO2011039457A1 (fr) | Procédé et dispositif de localisation dynamique d'un mobile | |
| WO2006067058A1 (fr) | Dispositif de determination autonome des coordonnees geographiques absolues d'un mobile evoluant en immersion | |
| EP1459979B1 (fr) | Procédé et dispositif pour déterminer au moins une information de position verticale d'un aéronef | |
| EP4492097A1 (fr) | Procede de determination d'une trajectoire d'un porteur mobile en temps reel et en temps differe | |
| CA2257349C (fr) | Recepteur de signaux de satellites avec controle d'integrite et exclusion d'axes defectueux | |
| CA2257194A1 (fr) | Recepteur de signaux de satellites avec detecteur d'incoherence entre mesures phase de code et de porteuse | |
| EP3534172A1 (fr) | Systeme de geolocalisation, aeronef et procede de geolocalisation associe | |
| EP3871005B1 (fr) | Procédé de correction d'une position précédemment estimée d'un véhicule | |
| FR3079608A1 (fr) | Procede de geolocalisation d'un drone, systeme de geolocalisation et produit programme d'ordinateur associes | |
| FR3152598A1 (fr) | Procédé de localisation tridimensionnelle d’un aéronef au sein d’une scène d’atterrissage comprenant au moins une piste d’atterrissage |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 08866617 Country of ref document: EP Kind code of ref document: A1 |
|
| DPE1 | Request for preliminary examination filed after expiration of 19th month from priority date (pct application filed from 20040101) | ||
| WWE | Wipo information: entry into national phase |
Ref document number: 2008866617 Country of ref document: EP |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 4195/CHENP/2010 Country of ref document: IN |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 12809636 Country of ref document: US |









