WO2005118366A1 - Processing of railway track data - Google Patents

Processing of railway track data Download PDF

Info

Publication number
WO2005118366A1
WO2005118366A1 PCT/GB2005/001600 GB2005001600W WO2005118366A1 WO 2005118366 A1 WO2005118366 A1 WO 2005118366A1 GB 2005001600 W GB2005001600 W GB 2005001600W WO 2005118366 A1 WO2005118366 A1 WO 2005118366A1
Authority
WO
WIPO (PCT)
Prior art keywords
data
sample
transfer function
track
stored
Prior art date
Application number
PCT/GB2005/001600
Other languages
French (fr)
Inventor
Sandor Matyas Patko
Original Assignee
Deltarail Group Limited
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 Deltarail Group Limited filed Critical Deltarail Group Limited
Priority to EP05738223A priority Critical patent/EP1771327A1/en
Priority to CA002573435A priority patent/CA2573435A1/en
Priority to AU2005249761A priority patent/AU2005249761A1/en
Priority to US11/628,311 priority patent/US20070246612A1/en
Publication of WO2005118366A1 publication Critical patent/WO2005118366A1/en

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B61RAILWAYS
    • B61KAUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
    • B61K9/00Railway vehicle profile gauges; Detecting or indicating overheating of components; Apparatus on locomotives or cars to indicate bad track sections; General design of track recording vehicles
    • B61K9/08Measuring installations for surveying permanent way

Definitions

  • This invention relates to an apparatus and a method for processing data, in particular data obtained by monitoring a railway track, such data for example being used for assessing the quality of the track.
  • Track recording vehicles are known, which are used in surveying a railway track to provide data representing the undulations of the rails in the vertical and horizontal planes, and their curvature.
  • Software packages are also available, for example a software product under the trade mark VAMPIRE (from AEA Technology pic) , for predicting how a particular vehicle will respond when travelling at a particular speed along a track; such software packages, which may be referred to as vehicle dynamics simulations, require input data providing an undistorted representation of the track.
  • VAMPIRE from AEA Technology pic
  • vehicle dynamics simulations require input data providing an undistorted representation of the track.
  • the raw data obtained by the sensors on a track recording vehicle provide information about train movement, and can be processed to determine track data, in particular being filtered to distinguish between short wavelength data and long wavelength data. This filtration process may introduce phase differences.
  • a method of obtaining data on the quality of a railway track comprising: a) receiving from a track recording vehicle data concerning variations of a parameter, the data comprising samples, obtained in either the spatial or the temporal domain, which have been subjected to a filtration process having an associated transfer function (H) ; b) selecting a transfer function H ⁇ which inverts at least the phase differences of the transfer function H of the filter; c) temporarily storing a multiplicity (N) of sequentially-received samples in a memory, each said sample being stored with an indication of the corresponding position or time; d) generating an output data sample by calculating the integral of the product of the stored data samples with an impulse function (F) , wherein the impulse function is deduced from the selected transfer function H ⁇ according to the equation:
  • time (t) is the appropriate variable, or, if expressed in terms of distance (s) :
  • the multiplicity (N) is an odd number; and preferably the impulse function is centred on the middle sample of those stored, that is ((N+l)/2) th sample if N is odd.
  • the impulse function need not be a symmetrical function; it is ⁇ centred' in the sense that it is a function not of absolute time (or distance) but of the time (or distance) relative to that of a specific stored sample. The method described above enables a series of output data samples to be generated substantially in real-time, the only delay being that taken for the receipt of ((N+l)/2) samples.
  • This method may be performed within a track recording vehicle. For example it can enable amplitude and phase distortions of track geometry signals to be removed, so that the corrected signals can be used as input for a vehicle dynamics simulation. Another application is that, once amplitude and phase distortions of track geometry signals have been removed, the signals correctly represent the shape of track features such as dipped rail joints, and so can be used to guide track maintenance.
  • the method of the invention can also remove distortions due to anti-aliasing filters.
  • the present invention also provides an apparatus for performing this method.
  • the method of the invention may be used to provide input data to a vehicle dynamics simulator carried in a track recording vehicle, so that the simulator can deduce the risk of derailment of a particular type of vehicle in substantially real-time.
  • the vehicle dynamics simulator could give a warning signal if the corresponding simulated vehicle would be derailed.
  • the track survey vehicle can, substantially in real-time, provide warnings of track sections that would give high derailment risk for a particular type of vehicle at a particular speed.
  • Warnings might also be given if the simulated vehicle would subject passengers to unacceptable jolts, or if the simulated vehicle would subject the portion of track to unacceptable track forces, and such information could also be reported as soon as the vehicle has passed over that section of the track. This enables track maintenance to be targeted at those sections of track most in need of improvement.
  • an apparatus incorporating the present invention is installed in a track recording vehicle 10, that is to say a rail vehicle incorporating transducers monitoring displacements and accelerations of the bogie and/or the body as the vehicle 10 moves along the track 11.
  • a track recording vehicle 10 that is to say a rail vehicle incorporating transducers monitoring displacements and accelerations of the bogie and/or the body as the vehicle 10 moves along the track 11.
  • transducers monitoring displacements and accelerations of the bogie and/or the body as the vehicle 10 moves along the track 11.
  • it might incorporate an accelerometer monitoring vertical accelerations of the bogie, and a displacement transducer monitoring vertical displacement of the axle relative to the bogie; data from such transducers would enable undulations in the vertical plane of each rail of the track to be monitored.
  • accelerometers measuring horizontal accelerations, along with a displacement transducer to monitor the wheel relative to the bogie enable undulations of the track in the horizontal plane to be monitored.
  • Track recording vehicles normally incorporate several different transducers, data from the transducers being sampled every 1/8 m and digitized, and the output data may involve calculations that combine data from several such transducers.
  • the data is subjected to signal processing (represented diagrammatically by box 12) that includes filtration so as to generate track data, which would typically be displayed to an operator, for example using a graphical interface, and stored for subsequent processing.
  • signal processing represented diagrammatically by box 12
  • the data may also be stored in conjunction with data from other sensors, for example positional data from a GPS sensor.
  • the data typically would represent alignment (a measure of the offset of the rails from the required smooth curve, measured in mm) , and curvature (indicating the reciprocal of the radius of the curve followed by the track, measured in km -1 ) .
  • the cutoff wavelength is set at 70 m, horizontal displacements of shorter wavelength than this being treated as alignment, and horizontal displacements of longer wavelength being treated as curvature.
  • the data typically would represent "top” (a measure of the displacement of the rails from the required smooth curve, measured in mm) , and gradient (indicating the slope of the track, in mm/mm) .
  • the cutoff wavelength in this case is typically also set to 70 m.
  • the track data streams from the processor 12 representing alignment, curvature, and top (and possibly also gradient) , and possibly other data streams such as positional information are transmitted to a data post-processing server 14, and thence to a reporting server 16, and so to various display interfaces 18 and to a data store 20.
  • Data streams representing alignment, curvature, and top (and possibly also gradient) are also supplied by the post-processing server 14 to several different vehicle dynamics modules 22 (three such modules are represented) .
  • Each such module 22 consists of a microprocessor arranged to model the dynamics of a particular vehicle travelling along the track 11 at a particular speed. The output of these vehicle dynamics modules 22 is fed back to the data post-processing server 14, and is supplied to the reporting server 16 along with the corresponding track data (processed as described below) .
  • the processor 12 is used to separate high frequency (short wavelength) components from low frequency (long wavelength) components.
  • Analogue filters or digital infinite impulse response (IIR) filters can perform these tasks efficiently, but they introduce distortion. Methods are known to eliminate this phase distortion, either avoiding it by using finite impulse response (FIR) filters instead of IIR filters, or by back filtering the already distorted data with an identical IIR filter to restore the original phase content.
  • FIR finite impulse response
  • the server 14 performs signal shaping of the incoming data, and forwards it to the rest of the system for storage and/or further processing.
  • the signal processing method can deal with both spatially and temporally sampled data streams. It can also perform 'cross-domain' operations, as well, that is to say to perform temporally defined operations in spatially sampled (taken at equal distances) data, and vice versa.
  • the server 14 consists of: - Digital input and output interfaces - A buffer memory to store N samples of the data stream, including the measured value and a time or distance stamp, indicating the time or distance the measurement was taken.
  • the type of the stamp data depends on the actual operation: if temporal operation is needed, then time stamp, if spatial operation is needed the distance stamp has to be attached to each measured value.
  • the actual sampling method does not affect the operation of the filter. For example, usually the measurements are taken at equal distances, so if the vehicle speed is increasing, then the differences between the consecutive time stamps will decrease, but the system operation will not change.
  • - Memory to store the parameters of the calculations.
  • Arithmetic processing capability is possible to store the parameters of the calculations.
  • the details such as the data transfer protocols, memory type etc. must be adjusted to the system in which the server 14 is used. In certain cases it may be a separate instrument connected to the data bus of the measurement system, in other cases it may be fully integrated into the measurement system.
  • the operation of the server 14 is as follows: 1.
  • the samples of the incoming data are stored in an N- element first-in-first-out (FIFO) buffer, which is initialized with zeros as measured values .
  • FIFO first-in-first-out
  • Each new sample enters the first slot of the buffer, moving the previous measurements one slot forward.
  • the data that had been in the N th slot is deleted, since it is replaced by the one coming from the (N-l) th slot.
  • N is preferably odd.
  • Y(T 0 ) is the output data, time stamped as taken at T 0 .
  • T 0 is the actual time stamp of the ((N+l)/2) th data in the buffer. In a certain sense, the calculation above is centred on T 0 , and the output data stream is always delayed by (N+l)/2 samples.
  • Ti is the time stamp of the oldest (N th ) data in the buffer.
  • T 2 is the time stamp of the latest (1 st ) data in the buffer. It is also true, that ! ⁇ To ⁇ T 2 .
  • X(t) is the data stream stored in the buffer.
  • F(t) is the finite impulse function, derived from the desired restoration.
  • F(t) is integratable between any possible t values.
  • Equation 1 which is expressed above as an integral (implying continuous functions), must in practice be performed as a summation, by a suitable discrete calculation method. Since each sample is processed separately, and has an associated time stamp, if the time intervals or spatial distances between successive samples vary, or there are randomly missing samples, overall operation is not affected. This is a significant advantage. Eq. 1 is shown in the temporal domain. The formula is still valid in the spatial domain, where the time values have to be replaced with distance values:
  • the calculated output is forwarded for further processing.
  • the operation clearly depends on correctly determining the impulse function, F(t) or F(s).
  • the impulse function is defined from the desired system behaviour, described by a transfer function. Transfer functions are complex equations that describe the system behaviour as a function of the cyclic frequency, ⁇ . If H(j ⁇ ) is the transfer function of a filter, then:
  • H(j ⁇ ) I is the ratio of the output to the input amplitude
  • ⁇ (H(j ⁇ )) is angle of the phase delay, where j is the square root of -1.
  • Example 1 Restoring the original phase content of an anti-aliased signal
  • a track recording vehicle 10 will include various transducers which measure aspects of the vehicle movement, such as an accelerometer, gyroscope etc.
  • the signal from such a transducer which is an analogue signal
  • an anti-aliasing filter is low frequency pass analogue filters, eliminating the undesired frequency content.
  • the data processor 12 would then produce digital output signals by sampling the analogue signal at equal distances along the track.
  • Anti-aliasing is essential, but it introduces a non-linear phase delay of the incoming signal. This phase delay will distort the shape of the signal, until now back-filtering was only the way to restore the original phase content. However, back- filtering changes the amplitudes in the transition band and cannot be used if the results are needed in real time.
  • the transfer function H of the analogue antialiasing filter can be given by the amplitude and phase responses as a function of the cyclic frequency:
  • Example 2 Restoring broadband curvature signal from asymmetric versine input
  • curvature It is difficult to measure curvature directly, so different indirect methods are used.
  • One of them is asymmetric versine; the asymmetric versine, v, is measured by considering a fixed length chord between two points on the rail. The chord is divided by a point Y into two unequal parts, Li and L 2 , and v is the distance of the rail from the point Y measured along a line perpendicular to the chord.
  • Asymmetric versine is easy to measure both manually and automatically. It gives a broadband description of the lateral track geometry, recording both short and long wavelengths components in the same output.
  • a complicated transfer function is required, which also introduces phase distortion.
  • Previously-known methods were unable to give a proper reconstruction of curvature from versine in real time.
  • the server 14 can be configured to reproduce broadband curvature from digital asymmetric versine input in real time.
  • a track recording vehicle 10 might include several such vehicle dynamics modules 22 operating in parallel, for example twelve rather than the three modules 22 shown here. Operation of this one vehicle 10 is therefore equivalent to running a fleet of a dozen different vehicles that may use this particular route, each at their own speed, and each of the virtual vehicles is effectively instrumented for assessing the risk of derailment, and also other parameters such as passenger comfort, track forces, vehicle kinematic movements etc..
  • This information is obtained in real-time, and is reported as part of the data provided to the display interfaces 18 as soon as the track recording vehicle 10 has passed over a portion of the track 11.
  • the information is embedded in the same stream of data as the information on track geometry. Hence it can be readily interfaced to track management software .

Abstract

The quality of a railway track (11) may be assessed with transducers on a track recording vehicle (10). The received data are filtered in a way that introduces phase shifts, the filtration process having an associated transfer function (H). An inverting transfer function H<The is therefore selected which inverts at least the phasedifferences of the transfer function (H) of the filter. A multiplicity (N) of successive data samples are stored in a memory, each with an indication of the corresponding position or time, and an output data sample is calculated as the integral of the product of the stored data samples with an impulse function (F) centred on the middle stored sample. The impulse function F (T) is related to the inverting transfer function H<The . As each data sample is moved into the memory the oldest such sample is deleted, and on each occasion an output data sample is calculated. The resulting output data stream represents the original data, without the phase shifts that were caused by the filtration process.

Description

Processing of Railway Track Data
This invention relates to an apparatus and a method for processing data, in particular data obtained by monitoring a railway track, such data for example being used for assessing the quality of the track.
Track recording vehicles are known, which are used in surveying a railway track to provide data representing the undulations of the rails in the vertical and horizontal planes, and their curvature. Software packages are also available, for example a software product under the trade mark VAMPIRE (from AEA Technology pic) , for predicting how a particular vehicle will respond when travelling at a particular speed along a track; such software packages, which may be referred to as vehicle dynamics simulations, require input data providing an undistorted representation of the track. The raw data obtained by the sensors on a track recording vehicle provide information about train movement, and can be processed to determine track data, in particular being filtered to distinguish between short wavelength data and long wavelength data. This filtration process may introduce phase differences. Data from such a track recording vehicle can be subjected to a subsequent filtration process, referred to as "back filtering", to obtain accurate data about the track. However, this process requires all the data about an entire section of track (which might be say 200 km long) , and this entire data stream is then processed in reverse; clearly this can't be done in real-time.
According to the present invention there is provided a method of obtaining data on the quality of a railway track, the method comprising: a) receiving from a track recording vehicle data concerning variations of a parameter, the data comprising samples, obtained in either the spatial or the temporal domain, which have been subjected to a filtration process having an associated transfer function (H) ; b) selecting a transfer function Hτ which inverts at least the phase differences of the transfer function H of the filter; c) temporarily storing a multiplicity (N) of sequentially-received samples in a memory, each said sample being stored with an indication of the corresponding position or time; d) generating an output data sample by calculating the integral of the product of the stored data samples with an impulse function (F) , wherein the impulse function is deduced from the selected transfer function Hτ according to the equation:
Hτ(jα>)= JF(t jωtdt
if time (t) is the appropriate variable, or, if expressed in terms of distance (s) :
Figure imgf000004_0001
e) storing the next successive sample of data in the memory and deleting the oldest sample stored in the memory, and repeating the step of generating an output data sample; and
f) repeatedly performing the preceding step. Preferably the multiplicity (N) is an odd number; and preferably the impulse function is centred on the middle sample of those stored, that is ((N+l)/2)th sample if N is odd. It should be appreciated that the impulse function need not be a symmetrical function; it is λcentred' in the sense that it is a function not of absolute time (or distance) but of the time (or distance) relative to that of a specific stored sample. The method described above enables a series of output data samples to be generated substantially in real-time, the only delay being that taken for the receipt of ((N+l)/2) samples. By appropriately selecting the impulse function, F, the effect of the filtration process on phase, or indeed on both amplitude and phase of the data, can be eliminated.
This method may be performed within a track recording vehicle. For example it can enable amplitude and phase distortions of track geometry signals to be removed, so that the corrected signals can be used as input for a vehicle dynamics simulation. Another application is that, once amplitude and phase distortions of track geometry signals have been removed, the signals correctly represent the shape of track features such as dipped rail joints, and so can be used to guide track maintenance. The method of the invention can also remove distortions due to anti-aliasing filters. The present invention also provides an apparatus for performing this method.
For example, the method of the invention may be used to provide input data to a vehicle dynamics simulator carried in a track recording vehicle, so that the simulator can deduce the risk of derailment of a particular type of vehicle in substantially real-time. The vehicle dynamics simulator could give a warning signal if the corresponding simulated vehicle would be derailed. Hence the track survey vehicle can, substantially in real-time, provide warnings of track sections that would give high derailment risk for a particular type of vehicle at a particular speed.
Warnings might also be given if the simulated vehicle would subject passengers to unacceptable jolts, or if the simulated vehicle would subject the portion of track to unacceptable track forces, and such information could also be reported as soon as the vehicle has passed over that section of the track. This enables track maintenance to be targeted at those sections of track most in need of improvement.
The invention will now be further and more particularly described, by way of example only, and with reference to the accompanying drawings which represents as a block diagram apparatus incorporating the present invention.
In this example, an apparatus incorporating the present invention is installed in a track recording vehicle 10, that is to say a rail vehicle incorporating transducers monitoring displacements and accelerations of the bogie and/or the body as the vehicle 10 moves along the track 11. For example it might incorporate an accelerometer monitoring vertical accelerations of the bogie, and a displacement transducer monitoring vertical displacement of the axle relative to the bogie; data from such transducers would enable undulations in the vertical plane of each rail of the track to be monitored. Similarly accelerometers measuring horizontal accelerations, along with a displacement transducer to monitor the wheel relative to the bogie, enable undulations of the track in the horizontal plane to be monitored. Track recording vehicles normally incorporate several different transducers, data from the transducers being sampled every 1/8 m and digitized, and the output data may involve calculations that combine data from several such transducers. In any event the data is subjected to signal processing (represented diagrammatically by box 12) that includes filtration so as to generate track data, which would typically be displayed to an operator, for example using a graphical interface, and stored for subsequent processing. The data may also be stored in conjunction with data from other sensors, for example positional data from a GPS sensor.
As regards the lateral plane, the data typically would represent alignment (a measure of the offset of the rails from the required smooth curve, measured in mm) , and curvature (indicating the reciprocal of the radius of the curve followed by the track, measured in km-1) . Typically the cutoff wavelength is set at 70 m, horizontal displacements of shorter wavelength than this being treated as alignment, and horizontal displacements of longer wavelength being treated as curvature. As regards the vertical plane, the data typically would represent "top" (a measure of the displacement of the rails from the required smooth curve, measured in mm) , and gradient (indicating the slope of the track, in mm/mm) . The cutoff wavelength in this case is typically also set to 70 m.
In the apparatus shown, the track data streams from the processor 12 representing alignment, curvature, and top (and possibly also gradient) , and possibly other data streams such as positional information, are transmitted to a data post-processing server 14, and thence to a reporting server 16, and so to various display interfaces 18 and to a data store 20. Data streams representing alignment, curvature, and top (and possibly also gradient) are also supplied by the post-processing server 14 to several different vehicle dynamics modules 22 (three such modules are represented) . Each such module 22 consists of a microprocessor arranged to model the dynamics of a particular vehicle travelling along the track 11 at a particular speed. The output of these vehicle dynamics modules 22 is fed back to the data post-processing server 14, and is supplied to the reporting server 16 along with the corresponding track data (processed as described below) .
The data post-processing server 14 is programmed to subject the track data streams from the processor 12 to the filtration process of the invention.
As mentioned above, the processor 12 is used to separate high frequency (short wavelength) components from low frequency (long wavelength) components. Analogue filters or digital infinite impulse response (IIR) filters can perform these tasks efficiently, but they introduce distortion. Methods are known to eliminate this phase distortion, either avoiding it by using finite impulse response (FIR) filters instead of IIR filters, or by back filtering the already distorted data with an identical IIR filter to restore the original phase content. However, there are cases where a signal has already been distorted by an analogue or IIR filter, and an undistorted signal is required. This is taken to be the case here . The server 14 performs signal shaping of the incoming data, and forwards it to the rest of the system for storage and/or further processing. The signal processing method can deal with both spatially and temporally sampled data streams. It can also perform 'cross-domain' operations, as well, that is to say to perform temporally defined operations in spatially sampled (taken at equal distances) data, and vice versa. The server 14 consists of: - Digital input and output interfaces - A buffer memory to store N samples of the data stream, including the measured value and a time or distance stamp, indicating the time or distance the measurement was taken. The type of the stamp data depends on the actual operation: if temporal operation is needed, then time stamp, if spatial operation is needed the distance stamp has to be attached to each measured value. The actual sampling method (equal time or equal distances) does not affect the operation of the filter. For example, usually the measurements are taken at equal distances, so if the vehicle speed is increasing, then the differences between the consecutive time stamps will decrease, but the system operation will not change. - Memory to store the parameters of the calculations. - Arithmetic processing capability.
It will be appreciated that the details such as the data transfer protocols, memory type etc. must be adjusted to the system in which the server 14 is used. In certain cases it may be a separate instrument connected to the data bus of the measurement system, in other cases it may be fully integrated into the measurement system.
The operation of the server 14 is as follows: 1. The samples of the incoming data are stored in an N- element first-in-first-out (FIFO) buffer, which is initialized with zeros as measured values . Each new sample enters the first slot of the buffer, moving the previous measurements one slot forward. The data that had been in the Nth slot is deleted, since it is replaced by the one coming from the (N-l)th slot. N is preferably odd.
2. After the new data sample is inserted into the buffer, the following calculation is performed:
Figure imgf000010_0001
where : Y(T0) is the output data, time stamped as taken at T0. T0 is the actual time stamp of the ((N+l)/2)th data in the buffer. In a certain sense, the calculation above is centred on T0, and the output data stream is always delayed by (N+l)/2 samples. Ti is the time stamp of the oldest (Nth) data in the buffer. T2 is the time stamp of the latest (1st) data in the buffer. It is also true, that !<To<T2. X(t) is the data stream stored in the buffer. F(t) is the finite impulse function, derived from the desired restoration. F(t) is integratable between any possible t values.
It will be appreciated that Equation 1, which is expressed above as an integral (implying continuous functions), must in practice be performed as a summation, by a suitable discrete calculation method. Since each sample is processed separately, and has an associated time stamp, if the time intervals or spatial distances between successive samples vary, or there are randomly missing samples, overall operation is not affected. This is a significant advantage. Eq. 1 is shown in the temporal domain. The formula is still valid in the spatial domain, where the time values have to be replaced with distance values:
Y(S0)= }F(s-S0)x(s)ds s> Eq. lb
3. The calculated output is forwarded for further processing. The operation clearly depends on correctly determining the impulse function, F(t) or F(s). The impulse function is defined from the desired system behaviour, described by a transfer function. Transfer functions are complex equations that describe the system behaviour as a function of the cyclic frequency, ω. If H(jω) is the transfer function of a filter, then:
|H(jω) I is the ratio of the output to the input amplitude, φ(H(jω)) is angle of the phase delay, where j is the square root of -1.
The selected transfer function H is one that reverses at least the phase change, and may also be selected so as to return the amplitude to its original value. The relationship between the transfer function and F(t) is:
Figure imgf000011_0001
The equation above has to be solved for F(t) . Analytical and numerical solutions are both suitable to get a functional F(t), and some examples are discussed below. The final step is to define the size of the buffer memory. First we calculate 11 and Th, such that the following approximation will be true:
Hτ(jω)= |F(t)ejωtdt « JF(t)ejωtdt τ> Eq . 3
Once Ti and T are found, the size of the buffer (N) can be calculated as follows :
1. The temporal window (time period) over which integration is performed is Tw=T- ι. 2. The number of samples in this time period will change as the vehicle changes speed, but if the maximum speed of the vehicle is known, then the number of samples will not exceed N = Tw (Top speed) • (Samples per metre) 3. If the vehicle is going slower than the top speed, some of the stored samples will fall outside this specified time period. However, Eq. 3 shows that we can take F(t)=0 for such samples.
This derivation assumes operation in the temporal domain. If spatial domain operation is needed, F(s) can be generated by replacing the temporal terms with spatial terms, as in Eq. lb.
Example 1 : Restoring the original phase content of an anti-aliased signal
This describes an operation in the temporal domain..
As mentioned earlier, a track recording vehicle 10 will include various transducers which measure aspects of the vehicle movement, such as an accelerometer, gyroscope etc. Typically the signal from such a transducer, which is an analogue signal, is first fed into an anti-aliasing filter, in order to avoid interference of high frequency signals with the digital sampling rate, called aliasing. Anti-aliasing filters are low frequency pass analogue filters, eliminating the undesired frequency content. The data processor 12 would then produce digital output signals by sampling the analogue signal at equal distances along the track. Anti-aliasing is essential, but it introduces a non-linear phase delay of the incoming signal. This phase delay will distort the shape of the signal, until now back-filtering was only the way to restore the original phase content. However, back- filtering changes the amplitudes in the transition band and cannot be used if the results are needed in real time.
The transfer function H of the analogue antialiasing filter can be given by the amplitude and phase responses as a function of the cyclic frequency:
A(ω) = |H(jα>) and φ(ω) = φ(H (jω) Eq.4
where ω is in radians per second. These two functions can be analytically derived, or measured. We must select or create a transfer function Hτ which leaves the amplitude intact, but reverses the phase delay. Hence the amplitude and phase responses of the selected transfer function Hτ should be as follows: . Aτ(ω)= 1, φτ(ω)= -φ(ω). Eqs . 5
This is satisfied by the transfer function:
Hτ(jω) = cos (φτ (ω) ) +jsin (φτ (ω) ) = cos (-φ(ω) ) +jsin (-φ(ω) ) Eq.6 Knowing the target transfer function, F(t) and N can be calculated as described above. Once these have been calculated, the server 14 can restore the original phase content of the incoming signal.
Example 2 : Restoring broadband curvature signal from asymmetric versine input
This describes an operation in the spatial domain.
If we model the railway track as a planar curve, it may be described by its curvature. Curvature for any planar curve is defined as:
Figure imgf000014_0001
where p is the vector pointing to a location on the track, s is the path taken on the track. Usually, the track curvature is split into long and short wavelength parts: the long wavelength part describes the track design, all the bends and straight sections needed to lead the train from A to B, while the short wavelength part describes the local deviations from the design, affecting the ride quality along the track.
It is difficult to measure curvature directly, so different indirect methods are used. One of them is asymmetric versine; the asymmetric versine, v, is measured by considering a fixed length chord between two points on the rail. The chord is divided by a point Y into two unequal parts, Li and L2, and v is the distance of the rail from the point Y measured along a line perpendicular to the chord. Asymmetric versine is easy to measure both manually and automatically. It gives a broadband description of the lateral track geometry, recording both short and long wavelengths components in the same output. Unfortunately, to determine curvature from asymmetric versine a complicated transfer function is required, which also introduces phase distortion. Previously-known methods were unable to give a proper reconstruction of curvature from versine in real time.
The server 14 can be configured to reproduce broadband curvature from digital asymmetric versine input in real time.
The transfer function from curvature to versine is:
Figure imgf000015_0001
where ω is in radians per metre, and Lx and L2 are in metres .
The inverse transfer function:
Hvc(jω)=— is the required transfer function (i.e. the selected 'transfer function Hτ) , and hence F(s) and N can be calculated as described above. Once these are calculated, the server 14 can restore the original curvature .
It will be appreciated that a track recording vehicle 10 might include several such vehicle dynamics modules 22 operating in parallel, for example twelve rather than the three modules 22 shown here. Operation of this one vehicle 10 is therefore equivalent to running a fleet of a dozen different vehicles that may use this particular route, each at their own speed, and each of the virtual vehicles is effectively instrumented for assessing the risk of derailment, and also other parameters such as passenger comfort, track forces, vehicle kinematic movements etc.. This information is obtained in real-time, and is reported as part of the data provided to the display interfaces 18 as soon as the track recording vehicle 10 has passed over a portion of the track 11. The information is embedded in the same stream of data as the information on track geometry. Hence it can be readily interfaced to track management software . Although the method has been described as being performed within a track recording vehicle 10, and so giving information in real-time, it will also be appreciated that data previously obtained using a track recording vehicle 10 may be supplied later to such a phase and amplitude correction microprocessor (equivalent to the post processing server 14), and hence if desired to a plurality of vehicle dynamics modules 22.

Claims

Claims
1. A method of obtaining data on the quality of a railway track, the method comprising: a) receiving from a track recording vehicle data concerning variations of a parameter, the data comprising samples, obtained in either the spatial or the temporal domain, which have been subjected to a filtration process having an associated transfer function (H) ;
b) selecting a transfer function Hτ which inverts at least the phase differences of the transfer function H of the filter;
c) temporarily storing a multiplicity (N) of sequentially-received samples in a memory, each said sample being stored with an indication of the corresponding position or time;
d) generating an output data sample by calculating the integral of the product of the stored data samples with an impulse function (F), wherein the impulse function is deduced from the selected transfer function Hτ according to the equation: Hτ(jω)= JF(t>jωtdt
if time (t) is the appropriate variable, or, if expressed in terms of distance (s) :
Hτ(jω)= jF(s jωsds e) storing the next successive sample of data in the memory and deleting the oldest sample stored in the memory, and repeating the step of generating an output data sample; and
f) repeatedly performing the preceding step.
2. A method as claimed in claim 1 wherein the selected transfer function Hr is such as to reverse both the changes in phase and the changes in amplitude due to the filtration process.
3. A method is claimed in claim 1 or claim 2 wherein the impulse function is centred on the ( (N+l) /2) th stored sample.
4. A method as claimed in claim 3 wherein the multiplicity is an odd number.
5. A method as claimed in any one of the preceding claims wherein the method is performed within a track recording vehicle.
6. An apparatus for performing a method as claimed in any one of the preceding claims.
PCT/GB2005/001600 2004-06-02 2005-04-28 Processing of railway track data WO2005118366A1 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
EP05738223A EP1771327A1 (en) 2004-06-02 2005-04-28 Processing of railway track data
CA002573435A CA2573435A1 (en) 2004-06-02 2005-04-28 Processing of railway track data
AU2005249761A AU2005249761A1 (en) 2004-06-02 2005-04-28 Processing of railway track data
US11/628,311 US20070246612A1 (en) 2004-06-02 2005-04-28 Processing of Railway Track Data

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
GB0412215.6 2004-06-02
GBGB0412215.6A GB0412215D0 (en) 2004-06-02 2004-06-02 Processing of railway track data

Publications (1)

Publication Number Publication Date
WO2005118366A1 true WO2005118366A1 (en) 2005-12-15

Family

ID=32696471

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/GB2005/001600 WO2005118366A1 (en) 2004-06-02 2005-04-28 Processing of railway track data

Country Status (6)

Country Link
US (1) US20070246612A1 (en)
EP (1) EP1771327A1 (en)
AU (1) AU2005249761A1 (en)
CA (1) CA2573435A1 (en)
GB (1) GB0412215D0 (en)
WO (1) WO2005118366A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3042823A1 (en) * 2015-01-08 2016-07-13 SmartDrive Systems, Inc. System and method for aggregation display and analysis of rail vehicle event information
US9663127B2 (en) 2014-10-28 2017-05-30 Smartdrive Systems, Inc. Rail vehicle event detection and recording system
US9902410B2 (en) 2015-01-08 2018-02-27 Smartdrive Systems, Inc. System and method for synthesizing rail vehicle event information
US9908546B2 (en) 2015-01-12 2018-03-06 Smartdrive Systems, Inc. Rail vehicle event triggering system and method

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0601819D0 (en) * 2006-01-31 2006-03-08 Aea Technology Plc Track twist monitoring
CN102353717B (en) * 2011-06-28 2013-03-13 哈尔滨工业大学 Detection apparatus and method for steel rail failures based on characteristics of non-negative tensor resolution
AU2018246236B2 (en) * 2017-03-27 2023-11-30 Harsco Technologies LLC Track geometry measurement system with inertial measurement
AT521420A1 (en) * 2018-07-11 2020-01-15 Plasser & Theurer Export Von Bahnbaumaschinen Gmbh Method and system for monitoring a track

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0697672A2 (en) * 1994-08-19 1996-02-21 MAN Technologie Aktiengesellschaft Procedure to correct a measurement curve or a signal variation, system for carrying out the method, and use for reconstruction of the error position in railway tracks by means of geometrical relative measurements
WO2004009422A1 (en) * 2002-07-19 2004-01-29 Aea Technology Plc Assessment of railway track geometry

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4760797A (en) * 1985-02-20 1988-08-02 Southern Railway Company Method and apparatus for automated tie detection and tamping
US5358202A (en) * 1992-07-21 1994-10-25 Consolidated Rail Corporation Cab signal track code analyzer system
US5768317A (en) * 1995-05-08 1998-06-16 National Semiconductor Corporation Equalization filter compensating for distortion in a surface acoustic wave device
US6424150B2 (en) * 1999-03-17 2002-07-23 Southwest Research Institute Magnetostrictive sensor rail inspection system
DE19704598C1 (en) * 1997-02-07 1998-06-18 Bruker Analytische Messtechnik Process for obtaining an optical FT spectrum
US6715354B2 (en) * 1998-02-24 2004-04-06 Massachusetts Institute Of Technology Flaw detection system using acoustic doppler effect
EP1134945B1 (en) * 2000-03-07 2008-05-14 Alcatel Lucent Method to determine a channel characteristic and discrete wavelet receiver to perform the method
TW577207B (en) * 2002-09-17 2004-02-21 Via Tech Inc Method and circuit adapted for blind equalizer

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0697672A2 (en) * 1994-08-19 1996-02-21 MAN Technologie Aktiengesellschaft Procedure to correct a measurement curve or a signal variation, system for carrying out the method, and use for reconstruction of the error position in railway tracks by means of geometrical relative measurements
WO2004009422A1 (en) * 2002-07-19 2004-01-29 Aea Technology Plc Assessment of railway track geometry

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9663127B2 (en) 2014-10-28 2017-05-30 Smartdrive Systems, Inc. Rail vehicle event detection and recording system
EP3042823A1 (en) * 2015-01-08 2016-07-13 SmartDrive Systems, Inc. System and method for aggregation display and analysis of rail vehicle event information
US9487222B2 (en) 2015-01-08 2016-11-08 Smartdrive Systems, Inc. System and method for aggregation display and analysis of rail vehicle event information
US9902410B2 (en) 2015-01-08 2018-02-27 Smartdrive Systems, Inc. System and method for synthesizing rail vehicle event information
US9981674B1 (en) 2015-01-08 2018-05-29 Smartdrive Systems, Inc. System and method for aggregation display and analysis of rail vehicle event information
US9908546B2 (en) 2015-01-12 2018-03-06 Smartdrive Systems, Inc. Rail vehicle event triggering system and method

Also Published As

Publication number Publication date
AU2005249761A1 (en) 2005-12-15
CA2573435A1 (en) 2005-12-15
US20070246612A1 (en) 2007-10-25
EP1771327A1 (en) 2007-04-11
GB0412215D0 (en) 2004-07-07

Similar Documents

Publication Publication Date Title
US20070246612A1 (en) Processing of Railway Track Data
Lee et al. A mixed filtering approach for track condition monitoring using accelerometers on the axle box and bogie
JP5226437B2 (en) Road surface flatness measuring device
CN109655055B (en) Positioning method and device of rail inspection robot
CN110874450B (en) Railway bridge track irregularity calculation method based on vehicle-mounted monitoring
Seraj et al. RoVi: Continuous transport infrastructure monitoring framework for preventive maintenance
CN107200040B (en) Method and system for determining the vertical profile of a track surface
AU2015264781B2 (en) System for monitoring the operating conditions of a train
CN113358053A (en) Track irregularity detection and evaluation system and method, electronic equipment and railway vehicle
KR20110134547A (en) Rail irregularity measurement system using accelerometers and method thereof
JP6592827B2 (en) Apparatus, method, program, and recording medium for identifying weight of vehicle traveling on traffic road
JP2008013153A (en) Abnormality detection device for railway vehicle
JP3018907B2 (en) Measurement method of ride comfort and vehicle vibration of railway vehicles
DE102008062143B3 (en) Method for determining vertical track bed of rail-road traffic, involves interlinking vertical axle bearing path vectors, three-point longitudinal height vectors, and equally spaced stretching vectors, respectively
Naganuma et al. Development of an inertial track geometry measuring trolley and utilization of its high-precision data
CN105807084A (en) Track geometric parameter detection synchronization and calculation method
Baasch et al. In-service condition monitoring of rail tracks: On an on-board low-cost multi-sensor system for condition based maintenance of railway tracks
Mei et al. Measurement of absolute vehicle speed with a simplified inverse model
CN113415308A (en) Rail corrugation detection method and device
CN112722003B (en) Method and equipment for monitoring train derailment risk
WO2004009422A1 (en) Assessment of railway track geometry
CN105404708A (en) Actually measured acceleration-based outside period motivation and evaluation method for long-span bridge structures
WO2007088321A1 (en) Track twist monitoring
JP2007050805A (en) Data acquisition method for vehicle drive route management
CN110579781B (en) Vehicle positioning method and device and electronic equipment

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KM KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SM SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): BW GH GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2005738223

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2005249761

Country of ref document: AU

WWE Wipo information: entry into national phase

Ref document number: 2573435

Country of ref document: CA

WWE Wipo information: entry into national phase

Ref document number: 11628311

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

WWW Wipo information: withdrawn in national office

Country of ref document: DE

ENP Entry into the national phase

Ref document number: 2005249761

Country of ref document: AU

Date of ref document: 20050428

Kind code of ref document: A

WWP Wipo information: published in national office

Ref document number: 2005249761

Country of ref document: AU

WWP Wipo information: published in national office

Ref document number: 2005738223

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 11628311

Country of ref document: US