EP3205771B1 - Procédé de détermination de ripages d'un rail d'une voie ferrée en domaine absolu - Google Patents

Procédé de détermination de ripages d'un rail d'une voie ferrée en domaine absolu Download PDF

Info

Publication number
EP3205771B1
EP3205771B1 EP17155780.4A EP17155780A EP3205771B1 EP 3205771 B1 EP3205771 B1 EP 3205771B1 EP 17155780 A EP17155780 A EP 17155780A EP 3205771 B1 EP3205771 B1 EP 3205771B1
Authority
EP
European Patent Office
Prior art keywords
rail
sub
frequency
absolute
versine
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.)
Active
Application number
EP17155780.4A
Other languages
German (de)
English (en)
Other versions
EP3205771A1 (fr
Inventor
Alban LEYMARIE
Julien Faure
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
LEYFA MEASUREMENT
Original Assignee
Leyfa Measurement
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Leyfa Measurement filed Critical Leyfa Measurement
Publication of EP3205771A1 publication Critical patent/EP3205771A1/fr
Application granted granted Critical
Publication of EP3205771B1 publication Critical patent/EP3205771B1/fr
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E01CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
    • E01BPERMANENT WAY; PERMANENT-WAY TOOLS; MACHINES FOR MAKING RAILWAYS OF ALL KINDS
    • E01B35/00Applications of measuring apparatus or devices for track-building purposes
    • E01B35/06Applications of measuring apparatus or devices for track-building purposes for measuring irregularities in longitudinal direction
    • EFIXED CONSTRUCTIONS
    • E01CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
    • E01BPERMANENT WAY; PERMANENT-WAY TOOLS; MACHINES FOR MAKING RAILWAYS OF ALL KINDS
    • E01B33/00Machines or devices for shifting tracks, with or without lifting, e.g. for aligning track, for shifting excavator track
    • E01B33/02Machines or devices for shifting tracks, with or without lifting, e.g. for aligning track, for shifting excavator track for slewing, i.e. transversely shifting, in steps

Definitions

  • the present invention belongs to the field of roadway maintenance for guided transport systems. It relates more particularly to a method for determining shifts of a rail of a railway in absolute domain. The present invention finds a particularly advantageous application, although in no way limiting, in the maintenance of the railway tracks.
  • railway rails which provide both transverse guidance and vertical support trains, are subject to various constraints that can affect their structure, and therefore a fortiori their trajectory.
  • such displacements are determined by means of a method comprising a preliminary step of determining geometrical parameters characteristic of the curvature of said rail, as well as a step of calculating displacements to be imposed at points of said rail, called points. of control, in order to calibrate these geometric parameters on corresponding values of a reference profile, generally archived during the laying of the railway.
  • said arrows are determined along the rail by means of a measuring device, such as a lorry.
  • a measuring device such as a lorry.
  • the latter is comparable to a mechanical system comprising an input signal and an output signal which are respectively a representation of the intrinsic geometry of the rail along which it moves and said arrows raised, these input and output signals being mapped by means of a transfer function.
  • the present invention aims to overcome all or part of the disadvantages of the prior art, including those described above, by proposing a method that makes it possible to determine the shifts of a rail of a railway in absolute domain and which presents steps adapted to the calculation of the shifts from a faithful representation of the rail, as well as to the taking into account in the calculation of said shifts of defects of the rail distributed over an extended spectrum (wavelengths of the order from the meter to wavelengths of the order of a hundred meters) through the use of several transfer functions.
  • the invention relates to a method for determining shifts of a track rail in absolute domain, said rail being comparable to a curve sampled at successive control points in which arrows were obtained by means of a number N of transfer functions that are distinct from each other, N being greater than or equal to 1, and each arrow being associated with only one of said transfer functions so as to form N sets of arrows.
  • the method for determining shifts of an absolute domain track rail comprises one or more of the following characteristics, taken separately or in any technically possible combination.
  • said decomposition of the absolute profile estimated according to the intended route of the rail and the training of the rail is obtained from a predetermined frequency filtering of said absolute profile estimated according to a first frequency window and a second window, respectively. frequency, each of said frequency windows having a lower terminal and an upper terminal.
  • said first and second frequency windows are respectively a first wavelength interval and a second wavelength interval, the upper bound and the lower bound respectively of said first interval and said second interval. being equal to a terminal slaved to the curvilinear abscissa along the rail as well as to a predetermined value of slippage corridor.
  • a threshold value is associated with each frequency sub-window during said sub-step of frequency filtering, each of said frequency sub-windows further comprising a lower bound and an upper bound determining a frequency range. during which the variation of the module of said inverse transfer function associated with said frequency sub-window is bounded around 1 by said threshold value.
  • the respective lower and upper bounds of said frequency sub-windows are determined iteratively so that the meeting of said frequency sub-windows describe a continuous frequency spectrum similar to the spectrum of the rail.
  • said filtered intermediate absolute profiles are summed frequency component by frequency component during said substep of recomposition of said estimated absolute profile.
  • the figure 1 represents a flowchart of an example of implementation of a method for determining shifts of a rail of a railway in absolute domain.
  • the invention specifically targets a railway rail, but remains applicable to all types of rails, including those of guided transport systems intended to circulate on a monorail or multirail network.
  • a railway is a double rail network. Therefore, said invention applies to any one of the two rails of said railway, without loss of generality, and it being understood that once shifts determined for one of said two rails (preferably the rail of larger radius for the skilled person insofar as said rail of larger radius is in practice the curved guide rail), these shifts are also connectable to the other of said two rails thanks to the knowledge of the rail between said two rails .
  • each point of said rail curve is associated, in a manner known per se, with a Frenet marker comprising, on the one hand, a unit vector tangent to said curve at said point, and secondly a unit vector normal to said tangent vector.
  • the left and right sides of said rail are defined as being the sides situated on the left and the right of any guided transport system traveling along the rail in the direction of the mileage, from said first end to said second end.
  • the distance between two points of said curve is in the sense of the length of the segment connecting said two points.
  • said distance is differentiated from the curvilinear distance separating said two points, the latter being counted along the curve of the rail by means of the curvilinear abscissa s originating from one of said two points.
  • said signal R (X) does not correspond to an explicit representation of the rail curve, which can still be formulated from a point of mathematical view by observing that said function R (X) is multivalued, or else unequivocal (still differently, there is at least one line parallel to the Y-axis and intersecting the rail curve in at least 2 points) . It is known to those skilled in the art that such a representation is disadvantageous when it is important to study the geometry of the rail from a numerical point of view.
  • a parameterization of the rail curve is defined in the form of a signal (in other words a digital function) denoted Y t depending on the curvilinear abscissa along the rail curve.
  • Y t ⁇ ⁇ C s d s 2
  • Such a parameterization of the rail curve is well known to those skilled in the art and has the advantage of being explicit, that is to say that the function Y t ( s ) is monovalued, or univocal.
  • the function Y t ( s ) is defined only from Euclidean invariants which are respectively the curvilinear abscissa along the curve of the rail and the curvature of said curve.
  • This parameterization of the rail is therefore qualified as intrinsic absolute because it constitutes a representation of the geometry of the rail which is independent of its orientation in said longitudinal plane.
  • said intrinsic absolute parameterization is called the absolute profile of the rail.
  • said absolute profile is conventionally represented in an orthogonal coordinate system having for abscissa s and for ordinate Y t ( s ), and said absolute orthogonal reference.
  • an arrow is also defined at a given point of the rail, referred to as the point of the arrow, as being the distance between said point of the arrow and an auxiliary point of a chord, of predefined length, underlying a subset of said rail containing said point of the arrow.
  • said auxiliary point of the chord is the projection, on said chord, of the point of the arrow in the direction of the normal vector of the Frenet mark associated with said point of the arrow (which is still the direction along the ordinate of the mark of Frenet).
  • Said arrow is an algebraic measurement, that is to say a length with a positive or negative sign depending on whether the point of deflection is, conventionally, located on an element of the curve of the rail whose center of curvature is positioned respectively to the right or left of the rail.
  • said deflection also contributes to the characterization of the geometry of the rail insofar as its direction of variation, with respect to the position of the point of the arrow along the rail, is identical to that of the curvature, as described herein. -before.
  • the rope used to determine an arrow is of suitable length so that the two points of the rail positioned at the intersection of the rope and said rail, said points. adjacent, are equidistant from said point of the arrow of 10m (accordingly the rope has a length less than 20m).
  • Such a configuration is conventionally used in railway engineering, that is to say known to those skilled in the art, and is advantageous when, for example, a cord is used to measure the arrows in the field, the latter not being subjected to when no significant deformations under the effect of its weight.
  • chord length for example with a chord length of 10m and adjacent points positioned, in no way limiting, equidistant from the point of the arrow.
  • said point of the arrow is not equidistant from said adjacent points.
  • An arrow is measured (or, equivalently, read) by measuring means known to those skilled in the art, whose operation may be automatic, such as for example linear displacement sensors, or requiring human intervention.
  • the measurement of the deflection is therefore dependent on the geometry of the rail, but also above all on said measuring means used and their displacements along the rail so that it provides a parametrization of the rail described as relative, as opposed to the absolute profile described above.
  • the absolute profile and the arrows measured along the rail are matched in this order by means of a characteristic transfer function of said measurement means implemented.
  • the absolute profile of the rail constitutes an input signal of the measuring means, said transfer function performing mathematical operations from characteristic geometrical quantities of said absolute profile so as to provide an output signal of said measuring means, namely, said arrows.
  • said arrow signal is a function of the curvilinear abscissa along the rail curve.
  • Said arrows signal is conventionally represented in an orthogonal coordinate system having the abscissa s and ordinate the value of the arrow, and said orthogonal reference arrows.
  • a signal is a digital function that satisfies the necessary and sufficient conditions (regularity, periodicity) to harmonic development.
  • sinusoidal said harmonics ranging from low to high frequencies (respectively equivalently, large to small wavelengths) in a predetermined pitch, each of said sinusoidal harmonics being further characterized by its amplitude and phase.
  • the expression "frequential component" is also used instead of that of sinusoidal harmonic without this leading to confusion.
  • the fact of being in harmonic regime with respect to the description of signals is equivalent in a manner known to those skilled in the art to describe them in frequency regime.
  • the absolute profile of the rail (respectively the arrow signal) is a continuous function extending between said first and second ends of the rail, this function can also be defined as the restriction of a periodic function of period equal to the length of the rail.
  • the absolute profile of the rail (respectively the signal of arrows) satisfies the conditions necessary and sufficient for development in sinusoidal harmonics.
  • said measuring means matching the absolute profile of the rail with the arrows signal, are linear systems, continuous and invariant so that an output signal is a linear function of the signal of associated with it, and that the superposition principle applies.
  • Said auxiliary rope being defined so that its projection on the abscissa axis s in the absolute orthogonal coordinate system is of length L, and that its two ends project on the abscissa axis s at two respective abscissa points equal to curvilinear abscissa of said two adjacent points.
  • the absolute profile Y t ( s ) is an advantageous representation of the rail and its geometry. More particularly, the link between said absolute profile and the measurements of arrows along the rail, as explained above, means that "determination of shifts in absolute domain" the determination of shifts of said absolute profile Y t ( s ), these shifts being intended to vary the amplitude of said absolute profile Y t ( s ) so as to correct the defects of the rail.
  • shifting corridor corresponds to the maximum permissible deviation between shifts made respectively on the left and on the right of the rail.
  • Said corridor of shifts is typically a technical constraint set by a rail maintenance operator in order to respect, for example, template constraints and obstacles imposed on said rail.
  • the curve of the rail is sampled at successive control points in which arrows have previously been obtained by means of a number N of transfer functions that are distinct from one another, N being greater than or equal to 1, and each arrow being associated to only one of said transfer functions so as to form N sets of arrows.
  • these two transfer functions are respectively associated with two different measurement means, such as, for example, and by no means limiting, two lorries.
  • any means for measuring arrows is characterized by the length of rope used as well as by the relative position of an arrow point with respect to the points adjacent to said point of arrow.
  • the present invention is described with only theoretical limitation that an arrow point can not be confused with an adjacent point associated with it. This amounts to saying that, for each measuring means, said length of rope as well as said relative position of the arrow point, excluding adjacent points, are a priori left free.
  • there are configurations of the measuring means having known characteristics and, in some cases, sought by those skilled in the art.
  • the transfer function of said measuring means has a zero phase shift.
  • the transfer function of said measuring means has a non-zero phase shift.
  • the use of transfer functions having a zero phase shift or constant, of module with isolated zeros as well as decreasing slowly at long wavelengths is advantageous.
  • said control points are physically embodied on the railway by means of visually identifiable marks, such as, for example, terminals installed during the laying of said railway, and are intended to locate the locations. where the rail will be rectified by means of an arrow modification, thus of a shifting, said locations being therefore coincident with the arrow points as defined above.
  • the location of the control points is achieved by means of the measuring device described in the patent application. FR 14 50897 of the plaintiff.
  • the method of determining shifts in absolute domain comprises firstly a step 50, following the obtaining of the arrows at said control points of the rail, modeling the rail by a signal said absolute profile estimated by means of said arrows obtained.
  • a method of estimating the absolute profile of the rail adapted to the determination, for a transfer function associated with measuring means, of an inverse transfer function so that the application of said function reverse transfer to said arrows obtained provides a fair and accurate estimate of the absolute profile of the rail.
  • step 50 of modeling the rail with an estimated absolute profile is assumed available at any time a database comprising the control points along the rail curve, arrows at said points control and distributed according to N sets of arrows, N transfer functions as well as the estimated absolute profile of the rail.
  • N transfer functions are all distinct since they are associated with different measuring means, so that among said N transfer functions, there exists one and only one, so-called pass-through function, having a cut-off frequency higher than the respective cut-off frequencies of the other functions of transfer.
  • the only transfer function used in step 50 is also said pass-through function.
  • the improved operating conditions of the rail relate to the safety of movement of guided transports along the rail at prescribed speed as well as the comfort of passengers of said transport systems guided according to a set of standards, including in particular the standard NF EN 13848 -1.
  • the rail dressing meanwhile, is representative of rail defects that oppose the respect of safety conditions, and to a lesser extent the respect of user comfort.
  • said decomposition of the estimated absolute profile according to the intended route of the rail and the rail dressing is obtained from a predetermined frequency filtering of said absolute profile estimated according to a first frequency window respectively. and a second frequency window, each of said first and second frequency windows having a lower bound and an upper bound.
  • Said frequency filtering according to said first and second frequency windows thus corresponds respectively to a first and second band-pass filtering of the estimated absolute profile, so as to isolate the contributions of wavelengths, respectively contained in said frequency windows, to the amplitude. and at the phase of said estimated absolute profile.
  • said first and second frequency windows are respectively a first wavelength interval [3m, 70m] and a second wavelength interval [70m, 150m], so that the terminal of said first wavelength interval is equal to the lower limit of said second wavelength interval wavelengths.
  • the intended alignment of the rail is decomposed according to a sum of sinusoidal harmonics, the latter being of respective wavelengths between 3m and 70m (respectively between 70m and 150m), as well as of amplitudes and phases identical respectively to the amplitudes and phases of the sinusoidal harmonics of the estimated absolute profile having the same wavelengths. It is therefore understandable that the algebraic sum of the sinusoidal harmonics of the target layout and of the training corresponds to the development in sinusoidal harmonics of the estimated absolute profile.
  • the first wavelength interval can be separated into two sub-intervals which are respectively [3m, 25m] and [25m, 70m], so as to isolate with precision, if the maintenance of the rail requires it, sinusoidal harmonics of wavelengths between 3m and 25m representative of defects particularly problematic from the point of view of safety standards.
  • said first and second frequency windows are respectively a first wavelength interval and a second wavelength interval of which the upper and lower limits respectively are equal.
  • the frequency filtering of the estimated absolute profile is done by successive elements of the rail.
  • said elements are a priori of different geometrical nature, and moreover they each have their own defects, said slave terminal is a fortiori different depending on whether it is associated with such or such element of the rail.
  • the train of the rail is a signal function of the abscissa curvilinear and defined in pieces, each piece being associated with an element, so that said signal comprises discontinuities between two consecutive pieces. Therefore, it is known to those skilled in the art that said training can be made continuous by means of linear interpolation two consecutive pieces having a discontinuity, which is actually achieved.
  • step 100 is particularly advantageous insofar as it constitutes an adaptive band-pass filtering, that is to say that it makes it possible to select for each element of the rail frequency window terminals adapted to better isolate the rail defects.
  • the figure 2 represents another particular embodiment of the method of the figure 1 during which, during step 100, the target route is obtained by deconvolution of a rail outline, and the training is obtained by frequency filtering of intermediate training obtained as for him by difference between said estimated absolute profile and said intended route.
  • a rail outline is mapped to the absolute profile of said rail by means of a transfer function, said transfer function being associated with measuring means. used when laying the rail to raise said outline.
  • the measuring means used to provide a rail outline are pure 3-point type, with a length of 20m cord and the point of the arrow positioned equidistant from said adjacent points, so that the transfer function and the function of inverse transfer associated with said measuring means are known analytically.
  • the transfer function and the inverse transfer function associated with the means for measuring a rail outline are assumed to be known.
  • step 100 is broken down into several successive sub-steps.
  • step 100 comprises a sub-step 101 for determining the target plot by applying to a rail outline a reverse transfer function of the transfer function associated with said rail outline.
  • said target route is therefore different from the estimated absolute profile, and is representative of the intrinsic geometry of the rail during the laying of the railway.
  • Step 100 then comprises a sub-step 102 for determining an intermediate dressing by difference between said estimated absolute profile and said target plot.
  • Said intermediate dressing is a signal representative of the variations experienced by the rail over time due to its operation. In other words, he represents the variation of intrinsic geometry of the rail between the laying date of the rail and a date subsequent to said rail laying date. It should be noted, however, that the frequency spectrum of intermediate training is as extensive as that of the target route, or even the estimated absolute profile. This means in particular that said intermediate dressing comprises sinusoidal harmonics whose frequency is not representative of characteristic frequencies of rail defects as sought by those skilled in the art.
  • step 100 comprises a sub-step 103 for determining said rail dressing by frequency filtering said intermediate dressing.
  • said frequency filtering is effected by means of a frequency window so as to isolate rail defects considered to be of interest in order to correct the rail at the same time. medium of shifts.
  • said frequency filtering corresponds to a bandpass filtering.
  • said frequency window corresponds to a range of wavelengths [3m, 70m] or even [3m, 25m] as described above.
  • the upper limit of said frequency window is determined according to the same technical characteristics as those described above in said preferred embodiment of step 100.
  • step 100 is thus available the intended route of the rail representative rail track as desired under improved operating conditions. Also available is the rail dressing which reflects the defects of the rail, that is to say the deformations undergone by the latter over time and which, once added to said plot in the algebraic sense of the term, provide the profile absolute rail estimate as calculated in said step 50.
  • the method for determining shifts in absolute domain finally comprises a step 200 of determining at said control points of shifts so that at each checkpoint the sum of said shifting and the rail dressing is zero.
  • the determination of said shifts initially consists in creating a function called shifts function as being the opposite, in the algebraic sense of the term, of the shifts. rail dressing.
  • Said rail dressing being a function of the curvilinear abscissa
  • the function of shifts is obtained by multiplying said dressing by -1 so that at each point of the curvilinear abscissa along the rail, the sum of the dressing and the shifts function is zero.
  • the determination of the shifts at said control points is given by the value of said shifted function evaluated at the curvilinear abscissa of said control points.
  • Such a way of determining the shifts to said control points is advantageous because it is performed from a true geometric description of the rail, namely the estimated absolute profile of the latter. It is therefore an advantage over traditional methods of determining shifts, such as the Hallade method, in which shifts are calculated from arrow values which are a relative and not absolute parameterization of the rail. , said arrows being further obtained by applying at least one measuring means transfer function to the absolute profile of the rail.
  • the expression of said at least one transfer function used during said traditional methods of determining shifts is based on the measurement principle as theoretically implemented by said measuring means, moving away and the empirical reality resulting from the mechanical stresses experienced by said measuring means.
  • said at least one transfer function is that known to those skilled in the art, a device type 3 pure points.
  • step 200 the shifts obtained at the end of step 200 have been advantageously calculated from the estimated absolute profile describing accurately the rail as laid on the track, and not from its representation by means of arrows.
  • Such a way of proceeding makes it possible to gain in precision (of the order of the millimeter) in the determination of the values of the shifts.
  • the mathematical operation for determining the shifts from the estimated absolute profile is simple to perform insofar as it consists in taking the opposite of the rail dressing. This is again an advantage in terms of time-consuming iterative calculations made within traditional methods (a slip at a fixed control point of the rail being dependent on the values of the arrows at all the other control points located on either side of said fixed control point).
  • the determination of shifts as described in step 200 is advantageous because it makes it possible to correct the defects of the rail so that the position of the latter before and after correction of defects is not changed.
  • This allows in particular to respect the constraints of templates as well as obstacles imposed on the railway.
  • the periodicity of said sinusoidal harmonics implies that for the rail dressing to be canceled, the shifts are uniformly distributed to the left and right of the rail .
  • the shifts are algebraic distances, we deduce that the sum of the shifts to the right of the rail is equal to the inverse of the sum of the shifts to the left of the rail.
  • the process is adapted to have as many shifts to the right of the rail as to the left of the rail.
  • step 200 in another embodiment of step 200, and when it is not necessary to obtain a shifts function to determine shifts at points other than said control points by continuity of said function of shifts, the shifts are determined directly from the values of the rail dressing at said control points multiplied by -1.
  • This way of proceeding proves advantageous in the optics of a saving of time of computation since the function of shifts whole is not calculated.
  • an estimate of the shifting function can be obtained later, without making use of rail dressing, by interpolating said shifts determined at the control points of the rail (linear interpolation, splines, etc.). The quality of this estimate is conventionally dependent on the sampling rate of the rail at said control points.
  • the figure 3 represents a preferred embodiment of the method of figure 1 during which, in step 50, the estimated absolute profile of the rail is modeled by arrows taken along the rail by means of a number N of transfer functions, N being strictly greater than 1.
  • step 50 is broken down into several successive sub-steps.
  • step 50 comprises a substep 51a for estimating N intermediate absolute profiles by applying to each set of arrows a reverse transfer function of the transfer function associated with said set of arrows.
  • Each inverse transfer function is for example determined by means of the absolute profile estimation method described in the patent application. EP 2,806,065 .
  • Each of said N intermediate absolute profiles provides an approximation of the real absolute profile of the rail, the differences between these N intermediate absolute profiles resulting from the differences between the N transfer functions, and therefore N inverse transfer functions, to which they are respectively associated. That said N transfer functions, and thus also said N inverse transfer functions, are distinct from each other is for example due to the fact that they are associated with measuring means having distinct configurations, as mentioned above. before in the examples where the length of rope used and / or the position of the point of the arrow relative to said adjacent points vary.
  • Step 50 then comprises a sub-step 52a of frequency filtering of said N intermediate absolute profiles according to respectively N frequency sub-windows so as to obtain N filtered intermediate absolute profiles.
  • said frequency filtering according to said N frequency sub-windows corresponds to N bandpass filterings of the N intermediate absolute profiles. This is, for each intermediate absolute profile, to isolate the wavelength contributions contained in the frequency sub-window associated with said intermediate absolute profile.
  • the absolute profile of the rail can not be determined perfectly from a set of arrows associated with a transfer function, and therefore also with a reverse transfer function, insofar as said reverse transfer function has a bandwidth so that it behaves like a filter attenuating, or even cutting, certain wavelengths.
  • the inverse transfer function is thus unsuited to the restitution of geometric parameters and rail defects associated with said wavelengths too attenuated and cut, and therefore the end of the absolute profile of the rail in its entirety.
  • said inverse transfer function is adapted to correctly restore the sinusoidal harmonics (amplitudes and phases) of the absolute profile in a wavelength interval during which its module is variable. bounded around 1 in order to minimize any frequency distortion phenomenon. It is further known to those skilled in the art that it is advantageous to impose a monotonic stress on the variation of the modulus of the inverse transfer function on said wavelength range. This is not contemplated in the following description.
  • a threshold value V s i for i belonging to the discrete set [1, N], is associated with each frequency sub-window during said sub-step of frequency filtering so that the width of each frequency sub-window determines a frequency range (and therefore of equivalent manner a wavelength interval) during which the variation of the module of the inverse transfer function associated with said frequency sub-window is bounded around 1 by said threshold value so that the peak-to-peak amplitude of said module remains less than or equal to twice that threshold value.
  • the module of the inverse transfer function is in the interval 1 - V s i , 1 + V s i .
  • Such an implementation therefore amounts to applying, for each threshold value, a band-pass filtering according to said wavelength interval, the limits of said interval being determined by said threshold value.
  • the module of the inverse transfer function is in a range [0.8, 1.3].
  • each frequency sub-window advantageously makes it possible during said substep 52a to determine N filtered intermediate absolute profiles, each filtered intermediate absolute profile decomposing according to sinusoidal harmonics of frequencies (or wavelength equivalent) included in the frequency sub-window associated with said filtered intermediate absolute profile, these sinusoidal harmonics faithfully reproducing the amplitudes and phases of the harmonics of the same frequencies belonging to the absolute profile of the rail.
  • the respective lower and upper limits of said frequency sub-windows are determined iteratively by means of threshold values respectively associated with said sub-steps. frequency windows and so that the meeting of said frequency sub-windows describe a continuous frequency spectrum similar to the spectrum of the rail.
  • said N transfer functions are sorted in order of interest, said interest being a function of the measurement means most appropriate to the type of rail which is desired to obtain an estimated absolute profile.
  • This choice is left to the appreciation of the skilled person, and may for example consist of sorting said inverse transfer functions by decreasing chord length of the measuring means associated with them, or vice versa.
  • sorting said N transfer functions is equivalent to sorting in the same order said inverse transfer functions and said N intermediate absolute profiles obtained at the end of step 51a.
  • said N intermediate absolute profiles are sorted in order of decreasing interest, and are named in this order: first intermediate absolute profile, second intermediate absolute profile, etc.
  • a first filtered intermediate absolute profile is determined by filtering said first intermediate absolute profile by means of a first frequency sub-window.
  • the lower and upper limits of said first frequency sub-window are in turn determined by means of a first threshold value according to identical technical characteristics to those described above for the module of the inverse transfer function associated with said first profile.
  • absolute intermediary is bounded around 1.
  • a second filtered intermediate absolute profile is determined by filtering said second intermediate absolute profile by means of a second frequency sub-window.
  • the lower and upper limits of said second frequency sub-window are in turn determined by means of a second threshold value according to the same technical characteristics as those described above for the inverse transfer function module associated with said second profile.
  • intermediate absolute is bounded around 1, and moreover under the constraint that the lower bound of said second frequency sub-window is equal to the upper bound of said first frequency sub-window.
  • the lower and upper bounds of the other frequency sub-windows are then determined iteratively, according to the order of interest chosen, according to the same principle, so that the union of said frequency sub-windows describes a continuous frequency spectrum similar to the rail spectrum. .
  • the filtered intermediate absolute profiles respectively associated with said at least two frequency sub-windows comprise sinusoidal harmonics of identical respective frequencies since included in said common frequency range. Therefore, the amplitudes of said sinusoidal harmonics of respective common frequencies are added during the recomposition of the estimated profile, as described below, which leads to an energy contribution, at said common frequencies, overvalued for said estimated absolute profile.
  • said frequency sub-windows it is understood that it is also advantageous for said frequency sub-windows to be contiguous so that no sinusoidal harmonic, potentially carrying energy at a frequency belonging to the spectrum of the rail, is omitted during the recomposition of the profile. estimated absolute, as described below.
  • Step 50 then comprises a substep 53a of recomposition of said absolute profile estimated from said N filtered intermediate absolute profiles.
  • said filtered intermediate absolute profiles are summed frequency component by frequency component during said substep of recomposition of said estimated absolute profile.
  • the estimated absolute profile of the rail is obtained by summing all the sinusoidal harmonics of said N filtered intermediate absolute profiles. In this way, the estimated absolute profile obtained by summation is advantageously adapted to describe the real geometry of the rail.
  • the figure 4 represents another particular embodiment of the method of the figure 1 during which, in step 50, the estimated absolute profile of the rail is modeled by N sets of filtered arrows, N being strictly greater than 1.
  • step 50 is broken down into several successive sub-steps.
  • step 50 comprises a sub-step 51b of frequency filtering of said N sets of arrows respectively according to N frequency sub-windows so as to obtain N sets of filtered arrows.
  • said frequency filtering according to said N frequency sub-windows corresponds to N bandpass filterings N sets of arrows. This is for each set of arrows, to isolate the contributions of wavelengths contained in the frequency pane associated with said set of arrows.
  • Such frequency filtering is performed according to technical characteristics identical to those described above in sub-step 52a, namely that said sets of arrows of said substep 51b replace said intermediate absolute profiles of said substep 52a.
  • Step 50 then comprises a substep 52b for estimating N intermediate filtered absolute profiles by application to each set of filtered arrows of an inverse transfer function of the transfer function associated with the set of arrows whose filtering at step 51b provides said set of filtered arrows.
  • N filtered intermediate absolute profiles is carried out according to identical technical characteristics to those described above in sub-step 51a, namely that said filtered intermediate absolute profiles and filtered arrows sets of said substep 51b respectively replace said intermediate absolute profiles and sets of arrows of said substep 51b.
  • step 50 comprises a sub-step 53b of said recomposition absolute profile estimated from said N intermediate profiles absolute filtered.
  • Such recomposition of said estimated absolute profile is performed according to technical characteristics identical to those described above in the substep 53a.
  • the figure 5 represents another particular embodiment of the method of the figure 1 during which, during step 50, the estimated absolute profile of the rail is modeled by deconvolution of a filtered arrows signal obtained by recomposition of N sets of filtered arrows, N being strictly greater than 1.
  • step 50 is broken down into several successive sub-steps.
  • step 50 comprises a sub-step 51c of frequency filtering of said N sets of arrows respectively according to N frequency sub-windows so as to obtain N sets of filtered arrows.
  • Said sub-step 51c is strictly identical to said substep 51b described above.
  • each set of filtered arrows, and therefore also each transfer function associated with said set of filtered arrows is associated with a frequency range comprising a lower bound and an upper bound.
  • Step 50 then comprises a substep 52c for determining a filtered arrows signal by recomposing said N sets of filtered arrows.
  • a determination of said filtered arrows signal is carried out according to technical characteristics identical to those described above in sub-step 53a, namely that said filtered arrows signal and said N sets of filtered arrows of said substep 52c replace respectively said estimated absolute profile and said N filtered intermediate absolute profiles of said substep 53a.
  • the frequency spectrum of said filtered arrows signal covers the entire spectrum of the rail.
  • it is an artificial arrows signal insofar as it consists of an assembly of arrows signals respectively obtained from separate measuring means.
  • said filtered arrows signal in the frequency domain, it is theoretically possible to consider said filtered arrows signal as a signal obtained from theoretical measurement means whose input is the absolute profile of the rail.
  • Such theoretical measuring means establish the correspondence between said absolute profile of the rail and said filtered arrows signal by means of a theoretical transfer function.
  • step 50 comprises a sub-step 53c for determining said estimated absolute profile by applying to said filtered arrows signal a reverse transfer function of a theoretical transfer function admitting at input and output, respectively, said absolute profile of the rail and said filtered arrows signal.
  • the modules of the transfer functions associated with said sets of filtered arrows are restricted to the frequency ranges determined during said substep 51c so as to form respectively restricted modules.
  • Said restricted modules are then concatenated so as to form the module of said theoretical transfer function. It is thus understood that in this way the module of the theoretical transfer function is defined piece by piece, each piece being associated with only one of said ranges. frequency, and that it covers the entire spectrum of rail. Therefore, and in a second step, the module of the inverse transfer function is determined by inverting said restricted modules on their respective frequency ranges.
  • said estimated absolute profile is determined according to identical technical characteristics to those described in the patent application. EP 2,806,065 .
  • the invention has been described by considering a method of determining shifts of a rail of a railway.

Landscapes

  • Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Machines For Laying And Maintaining Railways (AREA)

Description

    DOMAINE TECHNIQUE
  • La présente invention appartient au domaine de la maintenance de voies de circulation pour des systèmes de transport guidés. Elle concerne plus particulièrement un procédé de détermination de ripages d'un rail d'une voie ferrée en domaine absolu. La présente invention trouve une application particulièrement avantageuse, bien que nullement limitative, dans la maintenance du tracé de voies ferrées.
  • ÉTAT DE LA TECHNIQUE
  • Au cours de leur utilisation, les rails de chemins de fer, qui assurent aussi bien le guidage transversal et le soutien vertical de trains, sont soumis à diverses contraintes pouvant affecter leur structure, et donc a fortiori leur trajectoire. En particulier, l'accroissement au cours du temps du trafic ferroviaire, mais aussi de la charge des convois, impactent de manière défavorable les profils notamment transversaux et longitudinaux des rails. C'est aussi le cas de facteurs externes à l'exploitation des rails, comme par exemple les mouvements de terrain ou encore des gradients de température élevés.
  • En résulte alors l'apparition de défauts, bien que faibles (de l'ordre de quelques millimètres) mais modifiant notamment localement la courbure des rails, limitant alors l'exploitation de la voie ferrée à une vitesse inférieure à une vitesse prescrite, et impactant la stabilité des trains avec des conséquences importantes sur le confort des utilisateurs voire potentiellement sur la sécurité des circulations.
  • Dans l'optique d'un comportement optimal de roulement sur les rails, il importe donc de corriger ces défauts, par exemple par rectification des rails de la voie ferrée. Par rectification d'un rail ferroviaire, on entend ici le fait de lui imposer localement des déplacements selon des directions privilégiées de l'espace.
  • De manière conventionnelle, de tels déplacements sont déterminés au moyen d'un procédé comportant une étape préliminaire de détermination de paramètres géométriques caractéristiques de la courbure dudit rail, ainsi qu'une étape de calcul de déplacements à imposer en des points dudit rail, dits points de contrôle, dans le but de calibrer ces paramètres géométriques sur des valeurs correspondantes d'un profil de référence, généralement archivées lors de la pose de la voie ferrée.
  • On connait différents procédés selon le principe général décrit ci-avant, notamment ceux visant à déterminer lors de ladite étape préliminaire des distances caractéristiques, dite flèches, mesurant l'écart entre le rail et des cordes réelles ou virtuelles sous-tendant chaque arc reliant deux points de contrôle répartis de part et d'autre, et de manière adjacente, d'un troisième point de contrôle. L'objectif à atteindre étant alors la détermination de distances correctives, dites ripages, à appliquer auxdits points de contrôle afin de modifier lesdites flèches. La méthode dite « Hallade », connue de l'homme du métier, repose sur ce principe général, tout comme par exemple le procédé décrit dans la demande de brevet FR 15 53409 .
  • Il est à noter que lesdites flèches sont déterminées le long du rail au moyen d'un dispositif de mesure, comme par exemple un lorry. Ce dernier est assimilable à un système mécanique comportant un signal d'entrée et un signal de sortie qui sont respectivement une représentation de la géométrie intrinsèque du rail le long duquel il se déplace et lesdites flèches relevées, ces signaux d'entrée et de sortie étant mises en correspondance au moyen d'une fonction de transfert.
  • L'utilisation du formalisme des flèches dans des procédés de ripages est certes avantageuse, et donc largement répandue, car elle permet de décrire des défauts de petites dimensions (quelques millimètres) relativement à certaines grandes distances caractéristiques du rail (plusieurs centaines de mètres pour des rayons de courbure), ce que ne permettent pas des relevés topographiques le long du rail pour ces raisons de rapports de distances, mais aussi de temps et de main d'oeuvre.
  • Il convient néanmoins de tenir compte du fait que ces flèches ne constituent pas une représentation la plus fidèle possible de la géométrie intrinsèque du rail. Elles sont en effet obtenues par application de ladite fonction de transfert à des paramètres de la géométrie intrinsèque du rail, ladite fonction de transfert étant classiquement basée sur un modèle théorique de mesure à trois points, et non pas sur les caractéristiques mécaniques propres dudit dispositif de mesure. Ledit modèle théorique introduit donc un biais entre la géométrie réelle du rail et sa représentation au moyen des flèches. En outre, ce biais peut se propager à la détermination des ripages puisque ces derniers sont calculés à partir desdites flèches, de sorte que la correction des défauts du rail peut être entachée d'erreurs et peut être différente de celle qu'il conviendrait d'appliquer au regard de la géométrie du rail.
  • En outre, l'efficacité des procédés de détermination de ripages connus de l'homme de métier est limitée par le fait que les flèches sont obtenues par une seule fonction de transfert telle que décrite ci-dessus. En effet, certaines caractéristiques fréquentielles de ladite fonction de transfert, telle que sa bande passante, font qu'elle se comporte comme un filtre coupant certaines longueurs d'onde, se révélant ainsi inadaptée à la représentation de défauts du rail associés auxdites longueurs d'onde coupées. Ces défauts ne sont donc pas pris en compte dans lesdits procédés de détermination de ripages et ne peuvent donc pas être corrigés.
  • Le document US 5,012,413 A divulgue un procédé de détermination de ripages similaire.
  • EXPOSÉ DE L'INVENTION
  • La présente invention a pour objectif de remédier à tout ou partie des inconvénients de l'art antérieur, notamment ceux exposés ci-avant, en proposant un procédé qui permette de déterminer les ripages d'un rail d'une voie ferrée en domaine absolu et qui présente des étapes adaptées au calcul des ripages à partir d'une représentation fidèle du rail, ainsi qu'à la prise en compte dans le calcul desdits ripages de défauts du rail répartis sur un spectre étendu (longueurs d'onde de l'ordre du mètre aux longueurs d'onde de l'ordre de la centaine de mètres) grâce à l'utilisation de plusieurs fonctions de transfert.
  • A cet effet, l'invention concerne un procédé de détermination de ripages d'un rail de voie ferrée en domaine absolu, ledit rail étant assimilable à une courbe échantillonnée en des points de contrôle successifs en lesquels ont été obtenues des flèches au moyen d'un nombre N de fonctions de transfert distinctes les unes des autres, N étant supérieur ou égal à 1, et chaque flèche étant associée à une seule desdites fonctions de transfert de sorte à former N ensembles de flèches.
  • En outre, ledit procédé de détermination de ripages d'un rail de voie ferrée en domaine absolu comporte les étapes successives suivantes de :
    • modélisation du rail par un signal dit profil absolu estimé au moyen desdites flèches obtenues,
    • décomposition du profil absolu en un premier sous-signal dit tracé visé du rail et représentatif de la trajectoire du rail telle qu'adaptée à satisfaire des conditions améliorées d'exploitation du rail, ainsi que d'un deuxième sous-signal dit dressage du rail et représentatif des déformations subies par ledit tracé visé du rail au cours du temps,
    • détermination auxdits points de contrôle de ripages de sorte qu'en chaque point de contrôle la somme dudit ripage et dudit dressage du rail est nulle.
  • Dans des modes particuliers de mise en oeuvre, le procédé de détermination de ripages d'un rail de voie ferrée en domaine absolu comporte l'une ou plusieurs des caractéristiques suivantes, prises isolément ou selon toutes les combinaisons techniquement possibles.
  • Dans un mode particulier de mise en oeuvre, ladite décomposition du profil absolu estimé selon le tracé visé du rail et le dressage du rail est obtenue à partir d'un filtrage fréquentiel prédéterminé dudit profil absolu estimé selon respectivement une première fenêtre fréquentielle et une deuxième fenêtre fréquentielle, chacune desdites fenêtres fréquentielles comportant une borne inférieure et une borne supérieure.
  • Dans un mode particulier de mise en oeuvre, lesdites première et deuxième fenêtres fréquentielles sont respectivement un premier intervalle de longueurs d'onde et un deuxième intervalle de longueurs d'onde, la borne supérieure et la borne inférieure respectivement dudit premier intervalle et dudit deuxième intervalle étant égales à une borne asservie à l'abscisse curviligne le long du rail ainsi qu'à une valeur prédéterminée de couloir de ripages.
  • Dans un mode particulier, une épure du rail a été obtenue lors de la pose ou de la maintenance dudit rail au moyen d'une fonction de transfert associée à ladite épure, et dans lequel l'étape de décomposition du profil absolu comporte les sous-étapes successives suivantes de :
    • détermination du tracé visé par application à une épure du rail d'une fonction de transfert inverse de la fonction de transfert associée à ladite épure du rail,
    • détermination d'un dressage intermédiaire par différence entre ledit profil absolu estimé et ledit tracé visé,
    • détermination dudit dressage du rail par filtrage fréquentiel dudit dressage intermédiaire.
  • Dans un mode particulier de mise en oeuvre, lorsque le nombre N de fonctions de transfert est strictement supérieur à 1, ladite étape de modélisation du rail comporte les sous-étapes successives suivantes de:
    • estimation de N profils absolus intermédiaires par application à chaque ensemble de flèches d'une fonction de transfert inverse de la fonction de transfert associée audit ensemble de flèches,
    • filtrage fréquentiel desdits N profils absolus intermédiaires selon respectivement N sous-fenêtres fréquentielles de sorte à obtenir N profils absolus intermédiaires filtrés,
    • recomposition dudit profil absolu estimé à partir desdits N profils absolus intermédiaires filtrés.
  • Dans un mode particulier de mise en oeuvre, lorsque le nombre N de fonctions de transfert est strictement supérieur à 1, ladite étape de modélisation du rail comporte les sous-étapes successives suivantes de:
    • filtrage fréquentiel desdits N ensembles de flèches selon respectivement N sous-fenêtres fréquentielles de sorte à obtenir N ensembles de flèches filtrés,
    • estimation de N profils absolus intermédiaires filtrés par application à chaque ensemble de flèches filtré d'une fonction de transfert inverse de la fonction de transfert associée à l'ensemble de flèches dont le filtrage lors de la sous-étape de filtrage fréquentiel fournit ledit ensemble de flèches filtré,
    • recomposition dudit profil absolu estimé à partir desdits N profils absolus intermédiaires filtrés.
  • Dans un mode particulier de mise en oeuvre, une valeur seuil est associée à chaque sous-fenêtre fréquentielle lors de ladite sous-étape de filtrage fréquentiel, chacune desdites sous-fenêtres fréquentielles comportant en outre une borne inférieure et une borne supérieure déterminant une plage fréquentielle au cours de laquelle la variation du module de ladite fonction de transfert inverse associée à ladite sous-fenêtre fréquentielle est bornée autour de 1 par ladite valeur seuil.
  • Dans un mode particulier de mise en oeuvre, les bornes inférieure et supérieure respectives desdites sous-fenêtres fréquentielles sont déterminées de manière itérative de sorte que la réunion desdites sous-fenêtres fréquentielles décrivent un spectre fréquentiel continu similaire au spectre du rail.
  • Dans un mode particulier de mise en oeuvre, lesdits profils absolus intermédiaires filtrés sont additionnés composante fréquentielle par composante fréquentielle au cours de ladite sous-étape de recomposition dudit profil absolu estimé.
  • Dans un mode particulier de mise en oeuvre, lorsque le nombre N de fonctions de transfert est strictement supérieur à 1, ladite étape de modélisation du rail comporte les sous-étapes successives suivantes de:
    • filtrage fréquentiel desdits N ensembles de flèches selon respectivement N sous-fenêtres fréquentielles de sorte à obtenir N ensembles de flèches filtrés,
    • détermination d'un signal de flèches filtré par recomposition desdits N ensembles de flèches filtrés,
    • détermination dudit profil absolu estimé par application audit signal de flèches filtré d'une fonction de transfert inverse d'une fonction de transfert théorique admettant en entrée et sortie respectivement un profil absolu, représentatif de la géométrie réelle du rail, et ledit signal de flèches filtré.
    PRÉSENTATION DES FIGURES
  • Les caractéristiques et avantages de l'invention seront mieux appréciés grâce à la description qui suit, description qui expose les caractéristiques de l'invention au travers de modes de mise en oeuvre préférés, qui n'en sont nullement limitatifs.
  • La description s'appuie sur les figures annexées qui représentent :
    • Figure 1 : une représentation d'un organigramme d'un exemple de mise en oeuvre d'un procédé de détermination de ripages d'un rail d'une voie ferrée en domaine absolu.
    • Figure 2 : une représentation d'un organigramme d'un exemple de mise en oeuvre du procédé de la figure 1 au cours duquel, lors d'une étape 100, le tracé visé est obtenu par déconvolution d'une épure du rail, et le dressage est obtenu par filtrage fréquentiel d'un dressage intermédiaire obtenu quant à lui par différence entre ledit profil absolu estimé et ledit tracé visé.
    • Figure 3 : une représentation d'un organigramme d'un exemple de mise en oeuvre du procédé de la figure 1 au cours duquel, lors de l'étape 50, le profil absolu estimé du rail est modélisé grâce à des flèches relevées le long du rail au moyen d'un nombre N de fonctions de transfert, N étant strictement plus grand que 1.
    • Figure 4 : une représentation d'un organigramme d'un exemple de mise en oeuvre du procédé de la figure 1, lors de l'étape 50, le profil absolu estimé du rail est modélisé grâce à N ensembles de flèches filtrés, N étant strictement plus grand que 1.
    • Figure 5 : une représentation d'un organigramme d'un exemple de mise en oeuvre du procédé de la figure 1, lors de l'étape 50, le profil absolu estimé du rail est modélisé par déconvolution d'un signal de flèches filtré obtenu par recomposition de N ensembles de flèches filtrés, N étant strictement plus grand que 1.
    DESCRIPTION DÉTAILLÉE D'UN MODE DE MISE EN OEUVRE DE L'INVENTION
  • La figure 1 représente un organigramme d'un exemple de mise en oeuvre d'un procédé de détermination de ripages d'un rail d'une voie ferrée en domaine absolu.
  • L'invention, telle que décrite dans le présent exemple de réalisation, vise spécifiquement un rail ferroviaire, mais reste applicable à tous types de rails, notamment ceux de systèmes de transports guidés destinés à circuler sur un réseau monorail ou multirails.
  • En particulier, une voie ferrée est un réseau double rails. Dès lors, ladite invention s'applique à l'un quelconque des deux rails de ladite voie ferrée, sans perte de généralité, et étant entendu qu'une fois des ripages déterminés pour un desdits deux rails (préférentiellement le rail de plus grand rayon pour l'homme du métier dans la mesure où ledit rail de plus grand rayon est en pratique le rail de guidage en courbe), ces ripages sont aussi reliables à l'autre desdits deux rails grâce à la connaissance de l'entre rail desdits deux rails.
  • De manière classique, un rail de ladite voie ferrée est assimilable dans un plan longitudinal relativement au sens de déplacement le long dudit rail à une courbe s'étendant entre une première extrémité et une deuxième extrémité, et comportant des éléments successifs tels que des alignements ainsi que des arcs de cercles reliés entre eux par des raccordements progressifs, chacun desdits éléments comportant une origine et une fin. En particulier, ces différents éléments participent à la caractérisation de la géométrie du rail dans la mesure où ils présentent une courbure, fonction du déplacement suivant leur trajectoire respective, qui est :
    • nulle pour les alignements,
    • constante non nulle pour les arcs de cercles,
    • de variation linéaire pour les raccordements progressifs.
  • De plus, on note que chaque point de ladite courbe du rail est associé, de manière connue en soi, à un repère de Frenet comportant, d'une part, un vecteur unitaire tangent à ladite courbe en ledit point, et d'autre part un vecteur unitaire normal audit vecteur tangent.
  • On définit pour la suite de la description les côtés gauche et droit dudit rail comme étant les côtés situés sur la gauche et la droite de tout système de transport guidé circulant, le long du rail dans le sens du kilométrage, de ladite première extrémité à ladite deuxième extrémité.
  • Il est à noter que dans la suite de la description, la distance entre deux points de ladite courbe est comprise au sens de la longueur du segment reliant lesdits deux points. Ainsi, ladite distance est différenciée de la distance curviligne séparant lesdits deux points, cette dernière étant comptée le long de la courbe du rail au moyen de l'abscisse curviligne s ayant pour origine l'un desdits deux points.
  • En munissant ledit plan longitudinal d'un repère cartésien comportant une origine quelconque fixe, notée O, ainsi que deux axes, un axe des abscisses X et un axe des ordonnées Y, lesdits deux axes passant par O et étant non parallèles entre eux, ladite courbe admet une représentation sous forme d'un signal (autrement dit d'une fonction numérique) donnant une ordonnée de ladite courbe dans ledit repère cartésien selon l'abscisse X, et est noté dans la suite de la description R(X). Il est cependant important de noter que dans ledit repère cartésien, et pour certaines configurations géométriques du rail, ledit signal R(X) ne correspond pas à une représentation explicite de la courbe du rail, ce qui peut encore être formulé d'un point de vue mathématique en constatant que ladite fonction R(X) est multivaluée, ou encore non univoque (dit encore autrement, il existe au mois une droite parallèle à l'axe des ordonnées Y et qui intersecte la courbe du rail en au moins 2 points). Il est connu de l'homme du métier qu'une telle représentation est désavantageuse lorsqu'il importe d'étudier d'un point de vue numérique la géométrie du rail.
  • Aussi, on définit pour la suite de la description une paramétrisation de la courbe du rail sous la forme d'un signal (autrement dit d'une fonction numérique) noté Yt dépendant de l'abscisse curviligne s le long de la courbe du rail. En notant C la courbure locale de la courbe du rail, qui est également une fonction de s de manière connue en soi, ledit signal Yt(s) est défini au moyen de la formule suivante : Y t s = C s d s 2
    Figure imgb0001
  • Une telle paramétrisation de la courbe du rail est bien connue de l'homme du métier et présente l'avantage d'être explicite, c'est-à-dire que la fonction Yt(s) est monovaluée, ou encore univoque. De cette manière, la fonction Yt(s) est définie uniquement à partir d'invariants euclidiens qui sont respectivement l'abscisse curviligne le long de la courbe du rail ainsi que la courbure de ladite courbe. Cette paramétrisation du rail est par conséquent qualifiée d'absolue intrinsèque car elle constitue une représentation de la géométrie du rail qui est indépendante de son orientation dans ledit plan longitudinal. Dans la suite de la description, ladite paramétrisation absolue intrinsèque est nommée profil absolu du rail. Par ailleurs, ledit profil absolu est classiquement représenté dans un repère orthogonal ayant pour abscisse s et pour ordonnée Yt(s), et dit repère orthogonal absolu.
  • On définit également pour la suite de la description une flèche en un point déterminé du rail, dit point de la flèche, comme étant la distance entre ledit point de la flèche et un point auxiliaire d'une corde, de longueur prédéfinie, sous-tendant un sous-ensemble dudit rail contenant ledit point de la flèche. En outre, ledit point auxiliaire de la corde est la projection, sur ladite corde, du point de la flèche selon la direction du vecteur normal du repère de Frenet associé audit point de la flèche (qui est encore la direction selon l'ordonnée du repère de Frenet).
  • Ladite flèche est une mesure algébrique, c'est-à-dire une longueur affectée d'un signe positif ou négatif selon que le point de flèche est, conventionnellement, situé sur un élément de la courbe du rail dont le centre de courbure est positionné respectivement à droite ou à gauche du rail. De plus ladite flèche participe aussi à la caractérisation de la géométrie du rail dans la mesure où son sens de variation, eu égard à la position du point de la flèche le long du rail, est identique à celui de la courbure, tel que décrit ci-avant.
  • On comprend ainsi que la caractérisation d'une flèche relève non seulement du point du rail en lequel celle-ci est mesurée, mais aussi tout autant de la longueur de la corde associée. Il est donc virtuellement possible de définir une infinité de flèches en un même point d'un rail, en fonction de la longueur de la corde employée et de la position dudit point du rail relativement aux deux points d'intersection de la corde et dudit rail.
  • Par exemple, et à titre d'exemple nullement limitatif, la corde utilisée pour déterminer une flèche est de longueur adaptée de telle sorte que les deux points du rail positionnés à l'intersection de la corde et dudit rail, dits points adjacents, sont équidistants dudit point de la flèche de 10m (en conséquence la corde a une longueur inférieure à 20m). Une telle configuration est classiquement utilisée en ingénierie ferroviaire, c'est-à-dire connue de l'homme du métier, et est avantageuse lorsque, par exemple, un cordeau est utilisé pour mesurer les flèches sur le terrain, ce dernier ne subissant dès lors pas de déformations sensibles sous l'effet de son poids. Rien n'exclut cependant d'avoir d'autres configurations de longueur de corde, par exemple avec une longueur de corde de 10m et des points adjacents positionnés, de manière nullement limitative, à équidistance du point de la flèche. Dans d'autres exemples, ledit point de la flèche n'est pas équidistant desdits points adjacents.
  • Une flèche est mesurée (ou, de manière équivalente, relevée) par des moyens de mesure connus de l'homme du métier, dont le fonctionnement peut être automatique, tel que par exemple des capteurs de déplacement linéaire, ou bien nécessitant une intervention humaine. La mesure de la flèche est par conséquent dépendante de la géométrie du rail, mais aussi surtout desdits moyens de mesure mis en oeuvre et de leurs déplacements le long du rail de sorte qu'elle fournit une paramétrisation du rail qualifiée de relative, par opposition au profil absolu décrit ci-avant.
  • Ledit profil absolu et les flèches relevées le long du rail sont mis en correspondance dans cet ordre au moyen d'une fonction de transfert caractéristique desdits moyens de mesure mis en oeuvre. De cette manière, le profil absolu du rail constitue un signal d'entrée des moyens de mesure, ladite fonction de transfert réalisant des opérations mathématiques à partir de quantités géométriques caractéristiques dudit profil absolu de sorte à fournir un signal de sortie desdits moyens de mesure, à savoir lesdites flèches. Il est à noter que, comme c'est déjà le cas pour le profil absolu, ledit signal de flèches est fonction de l'abscisse curviligne le long de la courbe du rail. Ledit signal de flèches est classiquement représenté dans un repère orthogonal ayant pour abscisse s et pour ordonnée la valeur de la flèche, et dit repère orthogonal des flèches.
  • Dans la suite de la description, on adopte la convention qu'un signal est une fonction numérique qui satisfait les conditions nécessaires et suffisantes (régularité, périodicité) à un développement en harmoniques sinusoïdales, lesdites harmoniques allant des basses aux hautes fréquences (respectivement, de manière équivalente, des grandes aux petites longueurs d'onde) selon un pas prédéterminé, chacune desdites harmoniques sinusoïdales étant en outre caractérisée par son amplitude et sa phase. Dans la suite de la description, on emploie aussi l'expression « composante fréquentielle » au lieu de celle d'harmonique sinusoïdale sans que cela n'entraîne de confusion. En outre, le fait de se placer en régime harmonique en ce qui concerne la description de signaux équivaut de façon connue de l'homme du métier à décrire ces derniers en régime fréquentiel.
  • Dans le présent mode de mis en oeuvre de l'invention, le profil absolu du rail (respectivement le signal de flèches) est une fonction continue s'étendant entre lesdites première et deuxième extrémités du rail, cette fonction pouvant en outre être définie comme la restriction d'une fonction périodique de période égale à la longueur du rail. A cet effet, le profil absolu du rail (respectivement le signal de flèches) satisfait bien les conditions nécessaires et suffisantes à un développement en harmoniques sinusoïdales.
  • On adopte aussi la convention selon laquelle lesdits moyens de mesure, mettant en correspondance le profil absolu du rail avec le signal de flèches, sont des systèmes linéaires, continus et invariants de sorte qu'un signal de sortie est une fonction linéaire du signal d'entrée qui lui est associé, et que le principe de superposition s'applique. Les correspondances qui en découlent entre amplitudes des signaux d'entrée et de sortie d'une part, et phases des signaux d'entrée et de sortie d'autre part, sont connues de l'homme du métier. Il est à noter qu'une telle convention est mise en défaut, et donc non applicable, dans l'hypothèse où lesdits moyens de mesure ont pour entrée non pas ledit profil absolu, mais ledit signal R(X). En effet, le caractère non univoque du signal R(X) implique que les moyens de mesure ne peuvent être modélisés au moyen de systèmes linéaires. Dès lors, le fait d'utiliser ledit profil absolu est clairement un avantage dans la perspective de développements numériques.
  • Il est par ailleurs avantageux d'utiliser le profil absolu Yt(s) pour décrire le rail à la place de la fonction R(X) puisqu'il permet d'estimer une flèche en un point du rail de manière simple à partir d'une mesure de distance directement dans ledit repère orthogonal absolu. En effet, une flèche en un point du rail d'abscisse curviligne sf , définie à partir d'une corde de longueur L prédéterminée ainsi que par les distances respectives des points adjacents au point de la flèche, s'obtient dans ledit repère orthogonal absolu en mesurant directement un écart d'ordonnée, à l'abscisse sf , entre la courbe Yt(s) et une corde annexe. Ladite corde annexe étant définie de sorte que sa projection sur l'axe des abscisses s dans le repère orthogonal absolu est de longueur L, et que ses deux extrémités se projettent sur l'axe des abscisses s en deux points d'abscisses respectives égales aux abscisses curvilignes desdits deux points adjacents. De telles caractéristiques sont bien connues de l'homme du métier.
  • On comprend donc ainsi que le profil absolu Yt(s) est une représentation avantageuse du rail et de sa géométrie. Plus particulièrement, le lien entre ledit profil absolu et les mesures de flèches le long du rail, tel qu'exposé ci-avant, fait qu'on entend par « détermination de ripages en domaine absolu » la détermination de ripages dudit profil absolu Yt(s), ces ripages étant destinés à faire varier l'amplitude dudit profil absolu Yt(s) de sorte à corriger les défauts du rail.
  • Enfin, on introduit pour la suite de la description la notion de couloir de ripages qui correspond à l'écart maximal admissible entre des ripages effectués respectivement à gauche et à droite du rail. Ledit couloir de ripages est typiquement une contrainte technique fixée par un opérateur de maintenance du rail dans le but de respecter, par exemple, des contraintes de gabarits ainsi que d'obstacles imposées audit rail.
  • Le procédé de détermination de ripages en domaine absolu se décompose en plusieurs étapes successives.
  • La courbe du rail est échantillonnée en des points de contrôle successifs en lesquels ont été préalablement obtenues des flèches au moyen d'un nombre N de fonctions de transfert distinctes les unes des autres, N étant supérieur ou égal à 1, et chaque flèche étant associée à une seule desdites fonctions de transfert de sorte à former N ensembles de flèches.
  • Lesdites fonctions de transfert sont respectivement associées à des moyens de mesure distincts les uns des autres. Par exemple, et à titre d'exemple de mise en oeuvre nullement limitatif, les flèches sont obtenues au moyen de deux fonctions de transfert :
    • une première fonction de transfert, préférentiellement associée à des premiers moyens de mesure configurés de sorte que la corde a une longueur de 4,4m, et chaque flèche obtenue par ladite première fonction de transfert est associée à un point de flèche équidistant desdits deux points adjacents,
    • une deuxième fonction de transfert, préférentiellement associée à des seconds moyens de mesure configurés de sorte que la corde a une longueur de 7,4m, et chaque flèche obtenue par ladite deuxième fonction de transfert est associée à un point de flèche distant de 2,2m d'un desdits deux points adjacents.
  • Préférentiellement, ces deux fonctions de transfert sont respectivement associées à deux moyens de mesure différents, comme par exemple, et à titre nullement limitatif, deux lorries.
  • Rien n'exclut cependant d'avoir d'autres fonctions de transfert respectivement associées à des moyens de mesure configurés différemment, comme par exemple des moyens de mesure dits configurés en base 10, c'est-à-dire de sorte que chaque flèche est relevée en un point de flèche équidistant de 5m desdits points adjacents.
  • De manière plus générale, tout moyen de mesure de flèches est caractérisé par la longueur de corde utilisée ainsi que par la position relative d'un point de flèche par rapport aux points adjacents audit point de flèche. Aussi, la présente invention est décrite avec pour seule limitation théorique le fait qu'un point de flèche ne peut être confondu avec un point adjacent qui lui est associé. Cela revient à dire que, pour chaque moyen de mesure, ladite longueur de corde ainsi que ladite position relative du point de flèche, hors points adjacents, sont a priori laissées libres. Toutefois, il convient de noter qu'il existe des configurations des moyens de mesure présentant des caractéristiques connues et, dans certains cas, recherchées par l'homme du métier. Par exemple, pour des moyens de mesure configurés de sorte qu'un point de flèche est positionné à équidistance des points adjacents audit point de la flèche, la fonction de transfert desdits moyens de mesure présente un déphasage nul. A contrario, si ledit point de la flèche n'est pas équidistant desdits points adjacents, la fonction de transfert desdits moyens de mesure présente un déphasage non nul. De préférence, lorsqu'il s'agit d'effectuer des opérations de convolution/déconvolution sur un signal, comme cela est décrit ultérieurement, il est connu de l'homme du métier que l'utilisation de fonctions de transfert présentant un déphasage nul ou constant, de module présentant des zéros isolés ainsi que décroissant lentement aux grandes longueurs d'onde, est avantageux.
  • De cette manière, on comprend que lesdits points de contrôle en lesquels les flèches sont préalablement relevées sont répartis le long du rail sans nécessairement être équidistant les uns des autres. On comprend aussi que si lesdites fonctions de transfert ont toutes en commun le même signal d'entrée, à savoir le profil absolu de la voie, chaque flèche relevée correspond à une sortie d'une unique fonction de transfert. Ainsi, les flèches se répartissent en N ensembles distincts, dits ensembles de flèches, selon les fonctions de transfert dont elles dépendent respectivement.
  • Dans un mode de mise en oeuvre, lesdits points de contrôle sont matérialisés physiquement sur la voie ferrée au moyen de marques identifiables visuellement, comme par exemple des bornes mises en place lors de la pose de ladite voie ferrée, et sont destinés à repérer les emplacements où le rail sera rectifié au moyen d'une modification de flèche, donc d'un ripage, lesdits emplacements étant dès lors confondus avec les points de flèches tels que définis ci-avant.
  • Dans un autre mode de mise en oeuvre, lorsque des dégradations sont subies par le rail, que ce soit pour causes naturelles ou bien d'origine humaine, il arrive que tout ou partie des points de contrôle disparaisse de sorte que la pérennisation de la maintenance de la voie ferrée nécessite l'implantation de nouveaux points de contrôle. A cet effet, et en référence à la demande de brevet FR 14 61093 de la demanderesse, il est connu un procédé de mesure de la courbe du rail adapté à réaliser des mesures justes et fidèles de la géométrie de la courbe du rail, ainsi qu'à mettre en forme lesdites mesures de sorte à fournir lesdits points de contrôle.
  • Dans un exemple de mise en oeuvre dudit procédé de mesure de la courbe du rail, la localisation des points de contrôle est réalisée au moyen du dispositif de mesure décrit dans la demande de brevet FR 14 50897 de la demanderesse.
  • Le procédé de détermination de ripages en domaine absolu comporte dans un premier temps une étape 50, consécutive à l'obtention des flèches en lesdits points de contrôle du rail, de modélisation du rail par un signal dit profil absolu estimé au moyen desdites flèches obtenues. A cet effet, et en référence à la demande de brevet EP 2 806 065 de la demanderesse, il est connu un procédé d'estimation du profil absolu du rail adapté à la détermination, pour une fonction de transfert associée à des moyens de mesure, d'une fonction de transfert inverse de sorte que l'application de ladite fonction de transfert inverse auxdites flèches obtenues fournie une estimation juste et fidèle du profil absolu du rail.
  • On comprend donc que lorsque le nombre N de fonctions de transfert est strictement plus grand que 1, le procédé décrit dans la demande de brevet EP 2 806 065 est adapté à l'estimation de N profils absolus grâce en premier lieu à la détermination de N fonctions de transfert inverses respectivement associées auxdites N fonctions de transfert. Aussi, on adopte dans la suite de la description la convention que ledit profil absolu estimé lors de l'étape 50 correspond à l'une desdites N estimations du profil absolu ou bien encore à une combinaison desdites N estimations du profil absolu. Dans ce dernier cas, ladite combinaison sera détaillée ultérieurement dans un mode préféré de mise en oeuvre de l'invention.
  • Ainsi, à l'issue de l'étape 50 de modélisation du rail par un profil absolu estimé, est supposée disponible à tout instant une base de données comportant les points de contrôle le long de la courbe du rail, des flèches relevées en lesdits points de contrôle et réparties selon N ensembles de flèches, N fonctions de transfert ainsi que le profil absolu estimé du rail. En outre, lorsque le nombre N est strictement supérieur à 1, lesdites N fonctions de transfert sont toutes distinctes puisqu'associées à des moyens de mesure distincts, de sorte que parmi lesdites N fonctions de transfert, il en existe une et une seule, dite fonction de transfert passante, ayant une fréquence de coupure supérieure aux fréquences de coupure respectives des autres fonctions de transfert. Lorsque N est égal à 1, l'unique fonction de transfert utilisée au cours de l'étape 50 est aussi ladite fonction de transfert passante.
  • Le procédé de détermination de ripages d'un rail en domaine absolu comporte ensuite une étape 100 de décomposition du profil absolu estimé en :
    • un premier sous-signal dit tracé visé du rail et représentatif de la trajectoire du rail telle qu'adaptée à satisfaire des conditions améliorées d'exploitation du rail,
    • un deuxième sous-signal dit dressage du rail et représentatif des déformations subies par ledit tracé visé du rail au cours du temps
  • Les conditions améliorées d'exploitation du rail ont trait à la sécurité de déplacement de transports guidés le long du rail à vitesse prescrite ainsi qu'au confort de voyageurs desdits systèmes de transports guidés selon un ensemble de normes, dont notamment la norme NF EN 13848-1. Le dressage du rail, quant à lui, est représentatif de défauts du rail qui s'opposent au respect des conditions de sécurité, et dans une moindre mesure au respect du confort des utilisateurs.
  • Dans un mode de mise en oeuvre de l'étape 100, ladite décomposition du profil absolu estimé selon le tracé visé du rail et le dressage du rail est obtenue à partir d'un filtrage fréquentiel prédéterminé dudit profil absolu estimé selon respectivement une première fenêtre fréquentielle et une deuxième fenêtre fréquentielle, chacune desdites première et deuxième fenêtres fréquentielles comportant une borne inférieure et une borne supérieure. Ledit filtrage fréquentiel selon lesdites première et deuxième fenêtres fréquentielles correspond ainsi à respectivement un premier et deuxième filtrages passe-bande du profil absolu estimé, de sorte à isoler les contributions de longueurs d'onde, respectivement contenues dans lesdites fenêtres fréquentielles, à l'amplitude et à la phase dudit profil absolu estimé.
  • Par exemple, et à titre nullement limitatif, lesdites première et deuxième fenêtres fréquentielles sont respectivement un premier intervalle de longueurs d'onde [3m, 70m] et un deuxième intervalle de longueurs d'onde [70m, 150m], de sorte que la borne supérieure dudit premier intervalle de longueurs d'ondes est égale à la borne inférieure dudit deuxième intervalle de longueurs d'ondes. Cela signifie que le tracé visé du rail (respectivement le dressage du rail) se décompose selon une somme d'harmoniques sinusoïdales, ces dernières étant de longueurs d'onde respectives comprises entre 3m et 70m (respectivement comprises entre 70m et 150m), ainsi que d'amplitudes et de phases identiques respectivement aux amplitudes et phases des harmoniques sinusoïdales du profil absolu estimé ayant mêmes longueurs d'onde. On comprend donc que la somme algébrique des harmoniques sinusoïdales du tracé visé et du dressage correspond au développement en harmoniques sinusoïdales du profil absolu estimé.
  • En outre, un tel choix des bornes desdits premier et deuxième intervalles de longueurs d'onde est avantageux car représentatif des harmoniques du profil absolu estimé contribuant majoritairement, d'une part, à s'opposer aux normes de sécurité, et d'autre part à représenter la trajectoire du rail en conditions améliorées. En effet, il est connu de l'homme du métier que la délimitation entre sécurité et confort se fonde sur le spectre des défauts dans la mesure où se sont les harmoniques des paramètres mesurés, donc ici les flèches, mais donc de manière équivalente le profil absolu estimé par linéarité des fonctions de transfert, qui excitent des modes propres des véhicules circulant sur le rail. Ainsi, il est notamment connu que :
    • ledit premier intervalle de longueurs d'onde correspond à de courtes longueurs d'onde au regard de la longueur totale du rail (plusieurs kilomètres), et sont ainsi les plus contraignantes en termes de sécurité, et dans une moindre mesure en terme de confort,
    • ledit deuxième intervalle de longueurs d'onde correspond à des longueurs d'onde suffisamment grandes pour être représentatives de la trajectoire du rail en conditions améliorées
  • Rien n'exclut, suivant d'autres exemples non détaillés ici, d'avoir d'autres intervalles de longueurs d'onde au moyen desquelles sont obtenus par filtrage fréquentiel du profil absolu estimé ledit tracé visé et ledit dressage.
  • Par exemple, le premier intervalle de longueurs d'onde peut être séparé en deux sous-intervalles qui sont respectivement [3m, 25m] et [25m, 70m], de sorte à isoler avec précision, si la maintenance du rail le nécessite, les harmoniques sinusoïdales de longueurs d'onde comprises entre 3m et 25m représentatives de défauts particulièrement problématiques du point de vue des normes de sécurité.
  • Dans un mode préférentiel de mise en oeuvre de l'étape 100, lesdites première et deuxième fenêtres fréquentielles sont respectivement un premier intervalle de longueurs d'onde et un deuxième intervalle de longueurs d'onde dont respectivement la borne supérieure et la borne inférieure sont égales à une borne asservie à l'abscisse curviligne le long du rail ainsi qu'à une valeur prédéterminée de couloir de ripages. A cet effet, il convient dans un premier temps de déterminer l'abscisse curviligne des origines et fins respectives des éléments composant le rail, ces éléments étant de type alignements, arcs de cercle et raccordements progressifs comme évoqué ci-avant. Les origines et fins desdits éléments, ainsi que leurs abscisses curvilignes respectives, s'obtiennent par exemple visuellement au moyen d'un signal de flèches représenté dans ledit repère orthogonal des flèches, ledit signal de flèche étant alors avantageusement, et de manière connue de l'homme du métier, celui qui est obtenu au moyen de ladite fonction de transfert passante. En effet, ladite fonction de transfert passante est adaptée à fournir un signal de flèches dont la représentation dans le repère orthogonal des flèches laisse apparaître, de manière suffisante nette pour être notée visuellement, l'enchaînement des différents éléments composant le rail.
  • Dès lors, la borne asservie est obtenue selon une règle d'asservissement stipulant que lorsque l'abscisse curviligne le long du rail est comprise dans un intervalle dont les bornes sont les abscisses curvilignes d'un élément du rail (alignement, cercle, raccordement progressif), ladite borne asservie est égale à :
    • la longueur le long du rail dudit élément (soit donc la longueur curviligne dudit élément) si l'amplitude crête à crête du dressage du rail, obtenu par filtrage fréquentiel du profil absolu estimé selon ledit premier intervalle de longueurs d'onde dont la borne supérieure est égale à ladite borne asservie, reste inférieure à ladite valeur prédéterminée du couloir de ripages,
    • sinon une valeur de longueur d'onde obtenue par descente itérative par rapport à la longueur curviligne dudit élément jusqu'à ce que l'amplitude crête à crête du dressage du rail, obtenu par filtrage fréquentiel du profil absolu estimé selon ledit premier intervalle de longueurs d'onde dont la borne supérieure est égale à ladite borne asservie, reste inférieure à ladite valeur prédéterminée du couloir de ripages.
  • Aussi, il convient d'observer que le filtrage fréquentiel du profil absolu estimé, selon ledit mode préférentiel de mise en oeuvre de l'étape 100, se fait par éléments successifs du rail. Or, comme lesdits éléments sont a priori de nature géométrique différente, et que de plus ils comportent chacun des défauts qui leur sont propres, ladite borne asservie est a fortiori différente selon qu'elle est associée à tel ou tel élément du rail. Autrement dit, dans ledit mode préférentiel de mise en oeuvre de l'étape 100, le dressage du rail est un signal fonction de l'abscisse curviligne et défini par morceaux, chaque morceau étant associé à un élément, si bien que ledit signal comporte des discontinuités entre deux morceaux consécutifs. Dès lors, il est connu de l'homme du métier que ledit dressage peut être rendu continu au moyen d'une interpolation linéaire deux morceaux consécutifs présentant une discontinuité, ce qui est effectivement réalisé.
  • Il convent de noter que ledit mode préférentiel de mise en oeuvre de l'étape 100 est particulièrement avantageux dans la mesure où il constitue un filtrage passe-bande adaptatif, c'est-à-dire qu'il permet de sélectionner pour chaque élément du rail des bornes de fenêtres fréquentielles adaptées à isoler au mieux les défauts du rail.
  • La figure 2 représente un autre mode particulier de mise en oeuvre du procédé de la figure 1 au cours duquel, lors de l'étape 100, le tracé visé est obtenu par déconvolution d'une épure du rail, et le dressage est obtenu par filtrage fréquentiel d'un dressage intermédiaire obtenu quant à lui par différence entre ledit profil absolu estimé et ledit tracé visé.
  • Par épure, on entend un ensemble de flèches (et donc de manière équivalente un signal de flèches) mesurées en des points de contrôle, ledit ensemble ayant été relevé lors de la pose initiale ou la maintenance du rail. Une telle épure est représentative de la géométrie du rail telle que dimensionnée à son origine, c'est-à-dire avant que ledit rail ne subisse des déformations en raison de son exploitation, et est donc représentative d'une géométrie optimale dudit rail. En outre, il est connu de l'homme du métier qu'une épure du rail est typiquement consignée dans des cahiers d'exploitation et de maintenance du rail lors de la pose dudit rail, lesdits cahiers étant rendus accessibles à des personnes habilitées à les consulter.
  • Par ailleurs, il convient de noter qu'en tant que signal de flèches, une épure du rail est mise en correspondance avec le profil absolu dudit rail au moyen d'une fonction de transfert, ladite fonction de transfert étant associée à des moyens de mesure utilisés lors de la pose du rail pour relever ladite épure. Typiquement, les moyens de mesure utilisés pour fournir une épure du rail sont de type 3 points purs, avec une corde de longueur 20m et le point de la flèche positionné à équidistance desdits points adjacents, de sorte que la fonction de transfert et la fonction de transfert inverse associées auxdits moyens de mesure sont connues analytiquement. Rien n'exclut cependant qu'une épure du rail puisse être obtenue avec des moyens de mesure configurés selon des caractéristiques techniques différentes. Néanmoins, dans ledit mode particulier de mise en oeuvre de l'étape 100, sont supposées connues la fonction de transfert et la fonction de transfert inverse associées aux moyens de mesure d'une épure du rail.
  • Tel qu'illustré par la figure 2, ladite étape 100 se décompose en plusieurs sous-étapes successives.
  • Dans un premier temps, l'étape 100 comporte une sous-étape 101 de détermination du tracé visé par application à une épure du rail d'une fonction de transfert inverse de la fonction de transfert associée à ladite épure du rail.
  • A l'issue de ladite sous-étape 101, ledit tracé visé est donc différent du profil absolu estimé, et est représentatif de la géométrie intrinsèque du rail lors de la pose de la voie ferrée.
  • L'étape 100 comporte ensuite une sous-étape 102 de détermination d'un dressage intermédiaire par différence entre ledit profil absolu estimé et ledit tracé visé.
  • Ledit dressage intermédiaire est un signal représentatif des variations subies par le rail au cours du temps du fait de son exploitation. Autrement dit, il représente la variation de géométrie intrinsèque du rail entre la date de pose du rail et une date ultérieure à ladite date de pose du rail. Il convient toutefois de noter que le spectre fréquentiel du dressage intermédiaire est aussi étendu que celui du tracé visé, ou bien encore du profil absolu estimé. Cela signifie notamment que ledit dressage intermédiaire comporte des harmoniques sinusoïdales dont la fréquence n'est pas représentative de fréquences caractéristiques de défauts du rail tels que recherchés par l'homme du métier.
  • Aussi, l'étape 100 comporte une sous-étape 103 de détermination dudit dressage du rail par filtrage fréquentiel dudit dressage intermédiaire.
  • Dans un mode de mise en oeuvre de ladite sous-étape 103, ledit filtrage fréquentiel s'effectue au moyen d'une fenêtre fréquentiel de sorte à isoler des défauts du rail considérés comme étant d'intérêt dans l'optique de corriger le rail au moyen de ripages. Autrement dit, ledit filtrage fréquentiel correspond à un filtrage passe-bande. Par exemple, et à titre nullement limitatif, ladite fenêtre fréquentielle correspond à un intervalle de longueurs d'onde [3m, 70m] ou bien encore [3m, 25m] comme exposé ci-avant. Dans un autre exemple, la borne supérieure de ladite fenêtre fréquentielle est déterminée selon des caractéristiques techniques identiques à celles exposées ci-avant dans ledit mode préférentiel de mise en oeuvre de l'étape 100.
  • A l'issue de l'étape 100, est ainsi disponible le tracé visé du rail représentatif de la trajectoire du rail telle que désirée dans des conditions améliorées d'exploitation. Est aussi disponible le dressage du rail qui traduit les défauts du rail, c'est-à-dire les déformations subies par ce dernier au cours du temps et qui, une fois ajoutées audit tracé visé dans le sens algébrique du terme, fournissent le profil absolu estimé du rail tel que calculé au cours de ladite étape 50.
  • Le procédé de détermination de ripages en domaine absolu comporte enfin une étape 200 de détermination auxdits points de contrôle de ripages de sorte qu'en chaque point de contrôle la somme dudit ripage et du dressage du rail est nulle.
  • Dans un mode de mise en oeuvre de l'étape 200, la détermination desdits ripages consiste dans un premier temps à créer une fonction dite fonction de ripages comme étant l'opposé, au sens algébrique du terme, du dressage du rail. Ledit dressage du rail étant une fonction de l'abscisse curviligne, la fonction de ripages est obtenue en multipliant ledit dressage par - 1 de sorte qu'en chaque point de l'abscisse curviligne le long du rail, la somme du dressage et de la fonction de ripages est nulle. Puis, dans un deuxième temps, la détermination des ripages auxdits points de contrôle est donnée par la valeur de ladite fonction de ripages évaluée aux abscisses curvilignes desdits points de contrôle.
  • Une telle façon de déterminer les ripages auxdits points de contrôle est avantageuse car elle est effectuée à partir d'une véritable description géométrique du rail, à savoir le profil absolu estimé de ce dernier. Il s'agit donc d'un avantage par rapport aux méthodes traditionnelles de détermination de ripages, telle que par exemple la méthode Hallade, au cours desquelles des ripages sont calculés à partir de valeurs de flèches qui sont une paramétrisation relative et non absolue du rail, lesdites flèches étant en outre obtenues par application d'au moins une fonction de transfert de moyens de mesure au profil absolu du rail. Or, il convient de noter que l'expression de ladite au moins une fonction de transfert utilisée au cours desdites méthodes traditionnelles de détermination de ripages se base sur le principe de mesure tel que théoriquement mis en oeuvre par lesdits moyens de mesure, s'éloignant ainsi de la réalité empirique découlant des contraintes mécaniques subies par lesdits moyens de mesure. Par exemple, et à titre nullement limitatif, ladite au moins une fonction de transfert est celle, connue de l'homme du métier, d'un dispositif de type 3 points purs.
  • Ainsi, on comprend que les ripages obtenus à l'issue de l'étape 200 ont été avantageusement calculés à partir du profil absolu estimé décrivant avec précision le rail tel que posé sur la voie, et non pas à partir de sa représentation au moyen des flèches. Une telle manière de procéder permet de gagner en précision (de l'ordre du millimètre) dans la détermination des valeurs des ripages.
  • Par ailleurs, l'opération mathématique permettant de déterminer les ripages à partir du profil absolu estimé est simple à effectuer dans la mesure où elle consiste à prendre l'opposé du dressage du rail. Il s'agit là encore d'un avantage au regard des calculs itératifs coûteux en temps réalisés au sein des méthodes traditionnelles (un ripage en un point de contrôle fixé du rail étant dépendant des valeurs des flèches en tous les autres points de contrôle situés de par et d'autre dudit point de contrôle fixé).
  • Enfin, la détermination des ripages telle que décrite dans l'étape 200 est avantageuse car elle permet de corriger les défauts du rail de sorte que la position de ce dernier avant et après correction des défauts n'est pas modifiée. Ce qui permet notamment de respecter les contraintes de gabarits ainsi que d'obstacles imposées à la voie ferrée. En effet, dans la mesure où le dressage du rail consiste en une somme d'harmoniques sinusoïdales, la périodicité desdites harmoniques sinusoïdales impliquent que pour que le dressage du rail soit annulé, les ripages sont répartis de manière uniforme à gauche et à droite du rail. Comme les ripages sont des distances algébriques, on en déduit que la somme des ripages à droite du rail est égale à l'inverse de la somme des ripages à gauche du rail. Dit autrement, le procédé est adapté à avoir autant de ripages à droite du rail qu'à gauche du rail.
  • Alternativement, dans un autre mode de mise en oeuvre de l'étape 200, et lorsqu'il n'est pas nécessaire d'obtenir une fonction de ripages permettant de déterminer des ripages en des points autres que lesdits points de contrôle par continuité de ladite fonction de ripages, les ripages sont déterminés directement à partir des valeurs du dressage du rail en lesdits points de contrôle multipliées par -1. Cette manière de procéder se révèle avantageuse dans l'optique d'un gain de temps de calcul puisque la fonction de ripages toute entière n'est pas calculée. Toutefois, on note qu'une estimation de la fonction de ripages peut être obtenue ultérieurement, sans faire usage du dressage du rail, en interpolant lesdits ripages déterminés en les points de contrôle du rail (interpolation linéaire, splines, etc.). La qualité de cette estimation est classiquement dépendante du pas d'échantillonnage du rail en lesdits points de contrôle.
  • La figure 3 représente un mode préféré de mise en oeuvre du procédé de la figure 1 au cours duquel, lors de l'étape 50, le profil absolu estimé du rail est modélisé grâce à des flèches relevées le long du rail au moyen d'un nombre N de fonctions de transfert, N étant strictement plus grand que 1.
  • Dans ledit mode préféré de mise en oeuvre, et tel qu'illustré par la figure 3, ladite étape 50 se décompose en plusieurs sous-étapes successives.
  • Dans un premier temps, l'étape 50 comporte une sous-étape 51a d'estimation de N profils absolus intermédiaires par application à chaque ensemble de flèches d'une fonction de transfert inverse de la fonction de transfert associée audit ensemble de flèches. Chaque fonction de transfert inverse est par exemple déterminée au moyen du procédé d'estimation de profil absolu décrit dans la demande de brevet EP 2 806 065 .
  • Chacun desdits N profils absolus intermédiaires fournit une approximation du profil absolu réel du rail, les différences entre ces N profils absolus intermédiaires résultant des différences entre les N fonctions de transfert, et donc N fonctions de transfert inverses, auxquelles ils sont respectivement associés. Que lesdites N fonctions de transfert, et donc aussi lesdites N fonctions de transfert inverses, soient distinctes les unes par rapport aux autres est par exemple dû au fait qu'elles sont associées à des moyens de mesure ayant des configurations distinctes, comme évoqué ci-avant dans les exemples où la longueur de corde utilisée et/ou la position du point de la flèche par rapport auxdits points adjacents varient.
  • L'étape 50 comporte ensuite une sous-étape 52a de filtrage fréquentiel desdits N profils absolus intermédiaires selon respectivement N sous-fenêtres fréquentielles de sorte à obtenir N profils absolus intermédiaires filtrés.
  • Dans un mode principal de mise en oeuvre de ladite sous-étape 52a, ledit filtrage fréquentiel selon lesdites N sous-fenêtres fréquentielles correspond à N filtrages passe-bande des N profils absolus intermédiaires. Il s'agit là, pour chaque profil absolu intermédiaire, d'isoler les contributions de longueurs d'onde contenues dans la sous-fenêtre fréquentielle associée audit profil absolu intermédiaire.
  • De manière classique, le profil absolu du rail ne peut être déterminé de manière parfaite à partir d'un ensemble de flèches associées à une fonction de transfert, et donc aussi à une fonction de transfert inverse, dans la mesure où ladite fonction de transfert inverse comporte une bande passante de sorte qu'elle se comporte comme un filtre atténuant, voire coupant, certaines longueurs d'onde. La fonction de transfert inverse est ainsi inadaptée à la restitution de paramètres géométriques et de défauts du rail associés auxdites longueurs d'onde trop atténuées et coupées, et donc au final du profil absolu du rail dans son entièreté. En revanche, il est connu de l'homme du métier que ladite fonction de transfert inverse est adaptée à restituer correctement les harmoniques sinusoïdales (amplitudes et phases) du profil absolu dans un intervalle de longueurs d'onde au cours duquel son module est à variation bornée autour de 1afin de minimiser tout phénomène de distorsion fréquentielle. Il est en outre connu de l'homme du métier qu'il est avantageux d'imposer une contrainte de monotonie de la variation du module de la fonction de transfert inverse sur ledit intervalle de longueurs d'onde. Cela n'est pas envisagé dans la suite de la description.
  • Par exemple, une valeur seuil V s i ,
    Figure imgb0002
    pour i appartenant à l'ensemble discret [1,N], est associée à chaque sous-fenêtre fréquentielle lors de ladite sous-étape de filtrage fréquentiel de sorte que la largeur de chaque sous-fenêtre fréquentielle détermine une plage fréquentielle (et donc de manière équivalente un intervalle de longueurs d'onde) au cours de laquelle la variation du module de la fonction de transfert inverse associée à ladite sous-fenêtre fréquentielle est bornée autour de 1 par ladite valeur seuil de sorte que l'amplitude crête à crête dudit module reste inférieure ou égale au double de ladite valeur seuil. Autrement dit, le module de la fonction de transfert inverse est compris dans l'intervalle 1 V s i ,1 + V s i .
    Figure imgb0003
    Une telle mise en oeuvre revient donc à appliquer, pour chaque valeur seuil, un filtrage passe-bande selon ledit intervalle de longueurs d'onde, les bornes dudit intervalle étant déterminées par ladite valeur seuil. Dans un autre exemple, le module de la fonction de transfert inverse est compris dans un intervalle [0,8 , 1,3].
  • On comprend dès lors que choisir de façon adaptée la valeur seuil de chaque sous-fenêtre fréquentielle permet avantageusement lors de ladite sous-étape 52a de déterminer N profils absolus intermédiaires filtrés, chaque profil absolu intermédiaire filtré se décomposant selon des harmoniques sinusoïdales de fréquences (ou de manière équivalente de longueurs d'onde) comprises dans la sous-fenêtre fréquentielle associée audit profil absolu intermédiaire filtré, ces harmoniques sinusoïdales restituant de manière fidèle les amplitudes et phases des harmoniques de mêmes fréquences appartenant au profil absolu du rail.
  • Il est à noter que par « choisir de façon adaptée la valeur seuil de chaque sous-fenêtre fréquentielle », on entend choisir la valeur seuil de chaque sous-fenêtre fréquentielle de sorte que, d'une part, le module de chacune desdites fonctions de transfert inverses soit borné autour de 1 comme évoqué ci-avant, et aussi que, d'autre part, lesdites sous-fenêtres fréquentielles ne se chevauchent pas et permettent de décrire de manière continue un spectre fréquentiel similaire au spectre du profil absolu, dit spectre du rail.
  • Ainsi, dans un mode préférentiel de mise en oeuvre dudit mode principal de mise en oeuvre de la sous-étape 52a, les bornes inférieure et supérieure respectives desdites sous-fenêtres fréquentielles sont déterminées de manière itérative au moyen de valeurs seuils associées respectivement auxdites sous-fenêtres fréquentielles et de sorte que la réunion desdites sous-fenêtres fréquentielles décrivent un spectre fréquentiel continu similaire au spectre du rail.
  • Par exemple, et dans un premier temps, lesdites N fonctions de transfert sont triées par ordre d'intérêt, ledit intérêt étant fonction du moyen de mesure le plus approprié au type de rail dont on souhaite obtenir un profil absolu estimé. Ce choix est laissé à l'appréciation de l'homme du métier, et peut par exemple consister à trier lesdites fonctions de transfert inverses par longueur de corde décroissante des moyens de mesure qui leur sont associés, ou inversement. On remarque par ailleurs que trier lesdites N fonctions de transfert est équivalent à trier selon le même ordre lesdites fonctions de transfert inverses ainsi que lesdits N profils absolus intermédiaires obtenus à l'issue de l'étape 51a. Dans la suite de la description, lesdits N profils absolus intermédiaires sont triés par ordre d'intérêt décroissant, et sont nommés dans cet ordre : premier profil absolu intermédiaire, deuxième profil absolu intermédiaire, etc. Dès lors, et dans un deuxième temps, un premier profil absolu intermédiaire filtré est déterminé par filtrage dudit premier profil absolu intermédiaire au moyen d'une première sous-fenêtre fréquentielle. Les bornes inférieure et supérieure de ladite première sous-fenêtre fréquentielle sont quant à elles déterminées au moyen d'une première valeur seuil selon des caractéristiques techniques identiques à celles décrites ci-avant pour que le module de la fonction de transfert inverse associée audit premier profil absolu intermédiaire soit borné autour de 1. Puis, dans un troisième temps, un deuxième profil absolu intermédiaire filtré est déterminé par filtrage dudit deuxième profil absolu intermédiaire au moyen d'une deuxième sous-fenêtre fréquentielle. Les bornes inférieure et supérieure de ladite deuxième sous-fenêtre fréquentielle sont quant à elles déterminées au moyen d'une deuxième valeur seuil selon des caractéristiques techniques identiques à celles décrites ci-avant pour que le module de la fonction de transfert inverse associée audit deuxième profil absolu intermédiaire soit borné autour de 1, et par ailleurs sous la contrainte que la borne inférieure de ladite deuxième sous-fenêtre fréquentielle soit égale à la borne supérieure de ladite première sous-fenêtre fréquentielle. Les bornes inférieure et supérieure des autres sous-fenêtres fréquentielles sont ensuite déterminées itérativement, en suivant l'ordre d'intérêt choisi, selon le même principe de sorte que la réunion desdites sous-fenêtres fréquentielles décrivent un spectre fréquentiel continu similaire au spectre du rail.
  • On comprend en effet que lorsque par exemple au moins deux sous-fenêtres fréquentielles se chevauchent de sorte à définir une plage fréquentielle commune auxdites au moins deux sous-fenêtres fréquentielles, les profils absolus intermédiaires filtrés respectivement associés auxdites au moins deux sous-fenêtres fréquentielles comportent des harmoniques sinusoïdales de fréquences respectives identiques puisque comprises dans ladite plage fréquentielle commune. Dès lors, les amplitudes desdites harmoniques sinusoïdales de fréquences respectives communes s'ajoutent lors de la recomposition du profil estimé, comme décrit ci-après, ce qui conduit à une contribution énergétique, aux dites fréquences communes, surévaluée pour ledit profil absolu estimé.
  • Par ailleurs, on comprend qu'il est aussi avantageux que lesdites sous-fenêtres fréquentielles soient contigües de sorte qu'aucune harmonique sinusoïdale, potentiellement porteuse d'énergie à une fréquence appartenant au spectre du rail, ne soit omise lors de la recomposition du profil absolu estimé, comme décrit ci-après.
  • L'étape 50 comporte ensuite une sous-étape 53a de recomposition dudit profil absolu estimé à partir desdits N profils absolus intermédiaires filtrés.
  • Dans un mode de mise en oeuvre de ladite sous-étape 53a, lesdits profils absolus intermédiaires filtrés sont additionnés composante fréquentielle par composante fréquentielle au cours de ladite sous-étape de recomposition dudit profil absolu estimé.
  • Dans un exemple de mise en oeuvre, lorsque les N profils absolus intermédiaires filtrés se décomposent selon des harmoniques sinusoïdales de fréquences (ou de manière équivalente de longueurs d'onde) comprises dans les sous-fenêtres fréquentielles associées respectivement auxdits N profils absolus intermédiaires filtrés, de sorte que lesdites sous-fenêtres fréquentielles ne se chevauchent pas et permettent de décrire de manière continu un spectre fréquentiel similaire au spectre du rail, le profil absolu estimé du rail est obtenu en sommant toutes les harmoniques sinusoïdales desdits N profils absolus intermédiaires filtrés. De cette manière, le profil absolu estimé obtenu par sommation est avantageusement adapté à décrire la géométrie réelle du rail.
  • On comprend donc que l'utilisation d'une pluralité de fonctions de transfert, telle que décrite ci-avant dans l'étape 50 et illustré par la figure 2, permet de contourner des difficultés inhérentes à l'utilisation d'une seule fonction de transfert, dont notamment le fait que ladite une seule fonction de transfert est inadaptée à la représentation de l'ensemble des défauts du rail, et donc a fortiori à la détermination de ripages ajustés à la géométrie réelle du rail tout entier.
  • La figure 4 représente un autre mode particulier de mise en oeuvre du procédé de la figure 1 au cours duquel, lors de l'étape 50, le profil absolu estimé du rail est modélisé grâce à N ensembles de flèches filtrés, N étant strictement plus grand que 1.
  • Tel qu'illustré par la figure 4, ladite étape 50 se décompose en plusieurs sous-étapes successives.
  • Dans un premier temps, l'étape 50 comporte une sous-étape 51b de filtrage fréquentiel desdits N ensembles de flèches selon respectivement N sous-fenêtres fréquentielles de sorte à obtenir N ensembles de flèches filtrés.
  • Dans un mode de mise en oeuvre de ladite sous-étape 51b, ledit filtrage fréquentiel selon lesdites N sous-fenêtres fréquentielles correspond à N filtrages passe-bande des N ensembles de flèches. Il s'agit là, pour chaque ensemble de flèches, d'isoler les contributions de longueurs d'onde contenues dans la sous-fenêtre fréquentielle associée audit ensemble de flèches. Un tel filtrage fréquentiel est réalisé selon des caractéristiques techniques identiques à celles décrites ci-avant dans la sous-étape 52a, à savoir que lesdits ensembles de flèches de ladite sous-étape 51b remplacent lesdits profils absolus intermédiaires de ladite sous-étape 52a.
  • L'étape 50 comporte ensuite une sous-étape 52b d'estimation de N profils absolus intermédiaires filtrés par application à chaque ensemble de flèches filtré d'une fonction de transfert inverse de la fonction de transfert associée à l'ensemble de flèches dont le filtrage lors de l'étape 51b fournit ledit ensemble de flèches filtré. Une telle estimation desdits N profils absolus intermédiaires filtrés est réalisée selon des caractéristiques techniques identiques à celles décrites ci-avant dans la sous-étape 51a, à savoir que lesdits profils absolus intermédiaires filtrés et ensembles de flèches filtrés de ladite sous-étape 51b remplacent respectivement lesdits profils absolus intermédiaires et ensembles de flèches de ladite sous-étape 51b.
  • Enfin, l'étape 50 comporte une sous-étape 53b de recomposition dudit profil absolu estimé à partir desdits N profils absolus intermédiaires filtrés. Une telle recomposition dudit profil absolu estimé est réalisée selon des caractéristiques techniques identiques à celles décrites ci-avant dans la sous-étape 53a.
  • La figure 5 représente un autre mode particulier de mise en oeuvre du procédé de la figure 1 au cours duquel, lors de l'étape 50, le profil absolu estimé du rail est modélisé par déconvolution d'un signal de flèches filtré obtenu par recomposition de N ensembles de flèches filtrés, N étant strictement plus grand que 1.
  • Tel qu'illustré par la figure 5, ladite étape 50 se décompose en plusieurs sous-étapes successives.
  • Dans un premier temps, l'étape 50 comporte une sous-étape 51c de filtrage fréquentiel desdits N ensembles de flèches selon respectivement N sous-fenêtres fréquentielles de sorte à obtenir N ensembles de flèches filtrés. Ladite sous-étape 51c est strictement identique à ladite sous-étape 51b décrite ci-avant. En conséquence, chaque ensemble de flèches filtré, et donc aussi chaque fonction de transfert associée audit ensemble de flèches filtrés, est associé à une plage fréquentielle comportant une borne inférieure et une borne supérieure.
  • L'étape 50 comporte ensuite une sous-étape 52c de détermination d'un signal de flèches filtré par recomposition desdits N ensembles de flèches filtrés. Une telle détermination dudit signal de flèches filtré est réalisée selon des caractéristiques techniques identiques à celles décrites ci-avant dans la sous-étape 53a, à savoir que ledit signal de flèches filtré et lesdits N ensembles de flèches filtrés de ladite sous-étape 52c remplacent respectivement ledit profil absolu estimé et lesdits N profils absolus intermédiaires filtrés de ladite sous-étape 53a.
  • Il convient de noter que le spectre fréquentiel dudit signal de flèches filtré couvre l'ensemble du spectre du rail. Il s'agit néanmoins d'un signal de flèches artificiel dans la mesure où il consiste en un assemblage de signaux de flèches obtenus respectivement à partir de moyens de mesure distincts. Toutefois, dans le domaine fréquentiel, il est théoriquement possible de considérer ledit signal de flèches filtré comme un signal obtenu à partir de moyens de mesure théoriques dont l'entrée est le profil absolu du rail. De tels moyens de mesure théoriques établissent la correspondance entre ledit profil absolu du rail et ledit signal de flèches filtré au moyen d'une fonction de transfert théorique.
  • Ainsi, l'étape 50 comporte une sous-étape 53c de détermination dudit profil absolu estimé par application audit signal de flèches filtré d'une fonction de transfert inverse d'une fonction de transfert théorique admettant en entrée et sortie respectivement ledit profil absolu du rail et ledit signal de flèches filtré.
  • Dans un mode préféré de mise en oeuvre de ladite sous-étape 53c, et dans un premier temps, les modules des fonctions de transfert associées auxdits ensembles de flèches filtrés sont restreints aux plages fréquentielles déterminées lors de ladite sous-étape 51c de sorte à former respectivement des modules restreints. Lesdits modules restreints sont ensuite concaténés de sorte à former le module de ladite fonction de transfert théorique. On comprend donc qu'ainsi le module de la fonction de transfert théorique est défini morceau par morceau, chaque morceau étant associé à une seule desdites plages fréquentielles, et qu'il couvre bien l'ensemble du spectre du rail. Dès lors, et dans un deuxième temps, le module de la fonction de transfert inverse est déterminé par inversion desdits modules restreints sur leurs plages fréquentielles respectives. Enfin, dans un troisième temps, ledit profil absolu estimé est déterminé selon des caractéristiques techniques identiques à celles décrites la demande de brevet EP 2 806 065 .
  • De manière plus générale, il est à noter que les exemples de mise en oeuvre considérés ci-dessus ont été décrits à titre d'exemples non limitatifs, et que d'autres variantes sont par conséquent envisageables.
  • Par exemple, quelques modifications mineures dudit procédé permettent la détermination de relevages d'un rail d'une voie ferrée, c'est-à-dire de déplacements correctifs verticaux adaptés à corriger des défauts relevés sur le profil en long dudit rail.
  • Par ailleurs, l'invention a été décrite en considérant un procédé de détermination de ripages d'un rail d'une voie ferrée. Rien n'exclut, suivant d'autres exemples, d'avoir un procédé de rectification, selon des caractéristiques sensiblement similaires, appliqué à la rectification de poutres métalliques utilisées dans des structures de bâtiment ou bien encore de rails guidés d'ascenseur.

Claims (10)

  1. Procédé de détermination de ripages d'un rail de voie ferrée en domaine absolu, ledit rail étant :
    - assimilable à une courbe échantillonnée en des points de contrôle successifs en lesquels ont été obtenues des flèches au moyen d'un nombre N de fonctions de transfert distinctes les unes des autres, N étant supérieur ou égal à 1, et chaque flèche étant associée à une seule desdites fonctions de transfert de sorte à former N ensembles de flèches le procédé comprenant une étape de modélisation (50) du rail par un signal dit profil absolu estimé au moyen desdites flèches obtenues,
    caractérisé en ce qu'il comporte les étapes successives suivantes :
    - une étape de décomposition (100) du profil absolu estimé en une somme d'un premier sous-signal dit tracé visé du rail et représentatif de la trajectoire du rail telle qu'adaptée à satisfaire des conditions améliorées d'exploitation du rail, ainsi que d'un deuxième sous-signal dit dressage du rail et représentatif des déformations subies par ledit tracé visé du rail au cours du temps,
    - une étape de détermination (200) auxdits points de contrôle de ripages de sorte qu'en chaque point de contrôle la somme dudit ripage et dudit dressage du rail est nulle.
  2. Procédé selon la revendication 1, dans lequel ladite décomposition du profil absolu estimé selon le tracé visé du rail et le dressage du rail est obtenue à partir d'un filtrage fréquentiel prédéterminé dudit profil absolu estimé selon respectivement une première fenêtre fréquentielle et une deuxième fenêtre fréquentielle, chacune desdites fenêtres fréquentielles comportant une borne inférieure et une borne supérieure.
  3. Procédé selon la revendication 2, dans lequel lesdites première et deuxième fenêtres fréquentielles sont respectivement un premier intervalle de longueurs d'onde et un deuxième intervalle de longueurs d'onde, la borne supérieure et la borne inférieure respectivement dudit premier intervalle et dudit deuxième intervalle étant égales à une borne asservie à l'abscisse curviligne le long du rail ainsi qu'à une valeur prédéterminée de couloir de ripages.
  4. Procédé selon la revendication 1, dans lequel une épure du rail a été obtenue lors de la pose ou de la maintenance dudit rail au moyen d'une fonction de transfert associée à ladite épure, et dans lequel l'étape de décomposition (100) comporte les sous-étapes successives suivantes :
    - une sous-étape (101) de détermination du tracé visé par application à une épure du rail d'une fonction de transfert inverse de la fonction de transfert associée à ladite épure du rail,
    - une sous-étape (102) de détermination d'un dressage intermédiaire par différence entre ledit profil absolu estimé et ledit tracé visé,
    - une sous-étape (103) de détermination dudit dressage du rail par filtrage fréquentiel dudit dressage intermédiaire.
  5. Procédé selon l'une des revendications 1 à 4, dans lequel, lorsque le nombre N de fonctions de transfert est strictement supérieur à 1, ladite étape de modélisation (50) du rail comporte :
    - une sous-étape (51a) d'estimation de N profils absolus intermédiaires par application à chaque ensemble de flèches d'une fonction de transfert inverse de la fonction de transfert associée audit ensemble de flèches,
    - une sous-étape (52a) de filtrage fréquentiel desdits N profils absolus intermédiaires selon respectivement N sous-fenêtres fréquentielles de sorte à obtenir N profils absolus intermédiaires filtrés,
    - une sous-étape (53a) de recomposition dudit profil absolu estimé à partir desdits N profils absolus intermédiaires filtrés.
  6. Procédé selon l'une des revendications 1 à 4, dans lequel, lorsque le nombre N de fonctions de transfert est strictement supérieur à 1, ladite étape de modélisation (50) du rail comporte :
    - une sous-étape (51b) de filtrage fréquentiel desdits N ensembles de flèches selon respectivement N sous-fenêtres fréquentielles de sorte à obtenir N ensembles de flèches filtrés,
    - une sous-étape (52b) d'estimation de N profils absolus intermédiaires filtrés par application à chaque ensemble de flèches filtré d'une fonction de transfert inverse de la fonction de transfert associée à l'ensemble de flèches dont le filtrage lors de la sous-étape (51b) fournit ledit ensemble de flèches filtré,
    - une sous-étape (53b) de recomposition dudit profil absolu estimé à partir desdits N profils absolus intermédiaires filtrés.
  7. Procédé selon l'une des revendications 5 à 6, dans lequel une valeur seuil est associée à chaque sous-fenêtre fréquentielle lors de ladite sous-étape de filtrage fréquentiel, chacune desdites sous-fenêtres fréquentielles comportant en outre une borne inférieure et une borne supérieure déterminant une plage fréquentielle au cours de laquelle la variation du module de ladite fonction de transfert inverse associée à ladite sous-fenêtre fréquentielle est bornée autour de 1 par ladite valeur seuil.
  8. Procédé selon la revendication 7, dans lequel les bornes inférieure et supérieure respectives desdites sous-fenêtres fréquentielles sont déterminées de manière itérative de sorte que la réunion desdites sous-fenêtres fréquentielles décrivent un spectre fréquentiel continu similaire au spectre du rail.
  9. Procédé selon l'une des revendications 5 à 8, dans lequel lesdits profils absolus intermédiaires filtrés sont additionnés composante fréquentielle par composante fréquentielle au cours de ladite sous-étape de recomposition dudit profil absolu estimé.
  10. Procédé selon l'une des revendications 1 à 4, dans lequel, lorsque le nombre N de fonctions de transfert est strictement supérieur à 1, ladite étape de modélisation (50) du rail comporte :
    - une sous-étape (51c) de filtrage fréquentiel desdits N ensembles de flèches selon respectivement N sous-fenêtres fréquentielles de sorte à obtenir N ensembles de flèches filtrés,
    - une sous-étape (52c) de détermination d'un signal de flèches filtré par recomposition desdits N ensembles de flèches filtrés,
    - une sous-étape (53c) de détermination dudit profil absolu estimé par application audit signal de flèches filtré d'une fonction de transfert inverse d'une fonction de transfert théorique admettant en entrée et sortie respectivement un profil absolu du rail, représentatif de la géométrie réelle du rail, et ledit signal de flèches filtré.
EP17155780.4A 2016-02-12 2017-02-13 Procédé de détermination de ripages d'un rail d'une voie ferrée en domaine absolu Active EP3205771B1 (fr)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
FR1651177A FR3047814B1 (fr) 2016-02-12 2016-02-12 Procede de determination de ripages d'un rail d'une voie ferree en domaine absolu

Publications (2)

Publication Number Publication Date
EP3205771A1 EP3205771A1 (fr) 2017-08-16
EP3205771B1 true EP3205771B1 (fr) 2019-01-30

Family

ID=56372941

Family Applications (1)

Application Number Title Priority Date Filing Date
EP17155780.4A Active EP3205771B1 (fr) 2016-02-12 2017-02-13 Procédé de détermination de ripages d'un rail d'une voie ferrée en domaine absolu

Country Status (2)

Country Link
EP (1) EP3205771B1 (fr)
FR (1) FR3047814B1 (fr)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114777794B (zh) * 2022-03-28 2024-04-30 中国人民解放军国防科技大学 一种航天器轨道机动逆向移动滑窗检测方法、装置和设备

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3302428A (en) 1965-08-09 1967-02-07 Aldco Inc Device for cooling or keeping cool a beverage container
FR1450897A (fr) 1965-10-04 1966-06-24 Vestiaire pliant
US3493762A (en) 1967-02-27 1970-02-03 Andrew Eng Co Tracing machine having tool offset independent of the scanning head
EP0207197B1 (fr) * 1985-07-02 1989-03-29 Les Fils D'auguste Scheuchzer S.A. Procédé pour la réfection ou la pose d'une voie de chemin de fer
AT389133B (de) * 1987-01-14 1989-10-25 Vnii Transport Stroitelstvu Automatisches steuersystem fuer das richtorgan einer gleisrichtmaschine
US5012413A (en) * 1988-07-27 1991-04-30 Pandrol Jackson, Inc. Railroad track curve lining apparatus and method
FR3005321B1 (fr) * 2013-05-05 2015-10-02 Leyfa Measurement Dispositif de mesure de la geometrie d'une voie ferree et procede d'estimation des profils de nivellement et de dressage de ladite voie ferree

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
None *

Also Published As

Publication number Publication date
EP3205771A1 (fr) 2017-08-16
FR3047814A1 (fr) 2017-08-18
FR3047814B1 (fr) 2019-07-26

Similar Documents

Publication Publication Date Title
EP3162958B1 (fr) Procédé de mesure et d'estimation de géométrie d'une voie ferrée
CN103995288B (zh) 一种高斯束叠前深度偏移方法及装置
FR2872584A1 (fr) Methode pour simuler le depot sedimentaire dans un bassin respectant les epaisseurs des sequences sedimentaires
FR3031817A1 (fr) Methode de correction du calcul d'une caracteristique de vol d'un avion par prise en compte du vent vertical, procede de calcul du coefficient de trainee
Hashim et al. Toward an operating speed profile model for rural two-lane roads in Egypt
FR2984526A1 (fr) Controleur et procede pour diriger des sources
EP3454010B1 (fr) Procédé de compensation d'un champ magnétique, dispositif et programme d'ordinateur associés
EP1211152B1 (fr) Procédé et dispositif de localisation d'un véhicule sur une voie
EP3205771B1 (fr) Procédé de détermination de ripages d'un rail d'une voie ferrée en domaine absolu
EP2573515B1 (fr) Étalonnage de gyroscope vibrant
EP2459982B1 (fr) Procède de prévision d'un bruit de roulement d'un pneumatique
EP3650687A1 (fr) Procede de determination d'un facteur d'induction pour une eolienne equipee d'un capteur de teledetection par laser
EP3712621A1 (fr) Procede de prediction de la vitesse du vent dans le plan du rotor pour une eolienne equipee d'un capteur de teledetection par laser
FR2969753A1 (fr) Procede pour planifier des trajectoires aeroportees sous contrainte de performances plateforme et capteur
EP3781965B1 (fr) Dispositif, procede et produit programme d'ordinateur de validation de donnees fournies par une sonde pluviometrique
WO2022117765A1 (fr) Procédé et dispositif de génération de trajectoire d'un appareil mobile respectant une contrainte temporelle prédéterminée
FR2922321A1 (fr) Methode de simulation de positionnement d'au moins un streamer comprenant une etape de prediction d'un courant en au moins deux points distincts du streamer
WO2015011067A1 (fr) Procede de caracterisation de parametres mecaniques d'une chaussee
CA2464799A1 (fr) Methode pour determiner un modele de vitesse d'ondes sismiques dans une formation souterraine heterogene
EP4172648B1 (fr) Procede de determination des composantes de la vitesse du vent au moyen d'un capteur de teledetection par laser et au moyen d'une coherence temporelle
WO2020002814A1 (fr) Procédé d'établissement d'une cartographie d'une quantité d'un analyte dans l'environnement
EP2491424B1 (fr) Procede de localisation et de cartographie simultanees par filtrage non lineaire elastique
EP1826319B1 (fr) Procédé et système d'insertion d'éléments dans le sol, support d'enregistrement d'informations pour ce procédé
Terefenko et al. CCMORPH v2. 0–New coastal cliffs morphology analysis and visualization toolbox
EP2128731A1 (fr) Procédé de détermination de trajectoire en présence de courants variables dans le temps

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION HAS BEEN PUBLISHED

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20180216

RBV Designated contracting states (corrected)

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

INTG Intention to grant announced

Effective date: 20180822

RIN1 Information on inventor provided before grant (corrected)

Inventor name: LEYMARIE, ALBAN

Inventor name: FAURE, JULIEN

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: LEYFA MEASUREMENT

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

Free format text: NOT ENGLISH

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: AT

Ref legal event code: REF

Ref document number: 1093375

Country of ref document: AT

Kind code of ref document: T

Effective date: 20190215

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

Free format text: LANGUAGE OF EP DOCUMENT: FRENCH

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602017001936

Country of ref document: DE

REG Reference to a national code

Ref country code: NL

Ref legal event code: FP

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG4D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190530

Ref country code: NO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190430

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190501

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190530

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190430

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602017001936

Country of ref document: DE

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20191031

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

REG Reference to a national code

Ref country code: AT

Ref legal event code: UEP

Ref document number: 1093375

Country of ref document: AT

Kind code of ref document: T

Effective date: 20190130

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20170213

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20190130

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: IE

Payment date: 20240223

Year of fee payment: 8

Ref country code: NL

Payment date: 20240226

Year of fee payment: 8

Ref country code: LU

Payment date: 20240226

Year of fee payment: 8

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: AT

Payment date: 20240223

Year of fee payment: 8

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DE

Payment date: 20240226

Year of fee payment: 8

Ref country code: GB

Payment date: 20240229

Year of fee payment: 8

Ref country code: CH

Payment date: 20240301

Year of fee payment: 8

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 20240228

Year of fee payment: 8

Ref country code: BE

Payment date: 20240226

Year of fee payment: 8