CN102501886B  Method and system for detecting states of steel rail welding joints of highspeed railway  Google Patents
Method and system for detecting states of steel rail welding joints of highspeed railway Download PDFInfo
 Publication number
 CN102501886B CN102501886B CN201110344671.6A CN201110344671A CN102501886B CN 102501886 B CN102501886 B CN 102501886B CN 201110344671 A CN201110344671 A CN 201110344671A CN 102501886 B CN102501886 B CN 102501886B
 Authority
 CN
 China
 Prior art keywords
 axle box
 vertical acceleration
 box vertical
 welded joint
 effective value
 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
Links
Abstract
The invention discloses a method and system for detecting the states of steel rail welding joints of a highspeed railway. The method comprises the following steps of: sampling the time of the vertical acceleration of an axle box to obtain a time sampling signal of the vertical acceleration of the axle box; carrying out bandpass filtration on the time sampling signal of the vertical acceleration of the axle box; calculating an effective value of the filtered vertical acceleration of the axle box; determining the positions of the welding joints according to the effective value of the vertical acceleration of the axle box, and extracting the effective value of the vertical acceleration of the axle box at each welding joint; and if the effective value of the vertical acceleration of the axle box at each welding joint is larger than a threshold value, judging that the welding joint is in a poor state. The states of the welding joints on the whole railway are detected according to data of the vertical acceleration of the axle box in the technical scheme of the embodiment of the invention, therefore, the current using states of the welding joints of the highspeed railway can be effectively detected, and a powerful theoretical basis is provided for reasonably and economically guiding the maintenance and the repair of the welding joints.
Description
Technical field
The present invention relates to high speed railway safe operation technical field, particularly a kind of method of inspection and system of rail in high speed railway welded joint state.
Background technology
Steel rail weld joint is the weak link of high speed railway, has become the key equipment of restriction road speed and decision maintenance quality.The judge of steel rail weld joint is comprised to static judge and two kinds of dynamic judges, and main employing static state judge at present, dynamically passes judgment on and lacks.
Static judge operation is fairly simple, first utilizes the straight instrument of electronics to measure the flatness of steel rail weld joint, then according to certain criterion, passes judgment on its state.As, within the scope of the 1m centered by weld seam, if flatness height fluctuation difference in any 200mm section is not more than 0.2mm, thinks that state is qualified, otherwise need repairing.Fig. 1 is the static measurement result schematic diagram of the flatness of steel rail weld joint, and wherein upper and lower bound forms optimum waveform, and the ultimate range between them is called optimum peak value, conventionally value 0.2.
Static optimum waveform and the optimum peak value of passing judgment on rule of thumb determined, and lacks theoretical foundation.In addition, the impact of the flatness of welded joint on vehiclerail system dynamic characteristics, except relevant with its optimum peak value, simultaneously with the contact condition of wheel track, the suspension parameter of vehicle is closely related.Therefore, in order dynamically to grasp the state of steel rail weld joint, instruct maintenance, be necessary very much to carry out the dynamic judge of welded joint.
According to wavelength division, steel rail weld joint irregularity belongs to Short wave irregularity, wavelength is generally between 0.1～1m, and the detection wavelength of track geometry checking system is conventionally between 3～120m, the Short wave irregularity that is difficult to pass judgment on steel rail weld joint, on the safety of high speed railway trackVehicular system operation and the impact of traveling comfort, need to be studied other auxiliary judge means.
Axle box directly with wheel to being connected, track irregularity is by taking turns being directly delivered on axle box.Suppose wheel pairing approximation regard rigid construction as, the impact of the excitation power that axle box vertical acceleration can directly reflect the caused external world of track Short wave irregularity on vehicle dynamics.Therefore, select to adopt the auxiliary track Short wave irregularity of passing judgment on of axle box vertical acceleration on the impact of vehicle dynamics performance [list of references 17] both at home and abroad.Fig. 2 is the axle box vertical acceleration signal schematic representation of surveying on ripple abrasion rail, and Fig. 3 is the axle box vertical acceleration signal schematic representation on normal rail.
Although both at home and abroad carried out a large amount of research [list of references 810] utilizing aspect axle box vertical acceleration diagnosis rail wear, but how to utilize axle box vertical acceleration to pass judgment on welded joint state and yet there are no pertinent literature report and patent description, and the state of welded joint is significant for transporting safely of high speed railway.
Summary of the invention
For the defect of prior art, the embodiment of the present invention provides a kind of method of inspection and system of high speed railway welded joint state, and this scheme utilizes axle box vertical acceleration to detect the state of welded joint, has filled up the domestic and international technological gap in this area.
To achieve these goals, the embodiment of the present invention provides a kind of method of inspection of rail in high speed railway welded joint state, and described method comprises: axle box vertical acceleration is carried out to time sampling, obtain the time sampling signal of axle box vertical acceleration; Time sampling signal to described axle box vertical acceleration carries out bandpass filtering; The effective value of the axle box vertical acceleration after calculation of filtered; According to the effective value of described axle box vertical acceleration, determine the position of welded joint, and extract the effective value of the axle box vertical acceleration at each welded joint place; Whether the axle box vertical acceleration effective value that judges welded joint place is greater than threshold value, if be greater than threshold value, judges that this welded joint state is bad.
Before the time sampling signal to described axle box vertical acceleration carries out bandpass filtering, described method also comprises: the lowerfrequency limit F that presets bandpass filtering
_{l}, and the upper cutoff frequency F that dynamically determines bandpass filtering
_{h}; Wherein, dynamically determine the upper cutoff frequency F of bandpass filtering
_{h}specifically comprise: in frequency limit corresponding to the power spectrum of axle box vertical acceleration, axle box vertical acceleration is carried out to resonance and demodulation, obtain respectively the power envelope spectrum of the axle box vertical acceleration corresponding with each frequency in described frequency limit; Calculate the percentum that the first dominant frequency energy in each power envelope spectrum accounts for this power envelope spectrum gross energy, draw the relation curve of described percentum and its corresponding frequency; The frequency corresponding to peak value of percentum in described relation curve is defined as to the upper cutoff frequency of bandpass filtering.
Described in frequency limit corresponding to the power spectrum of axle box vertical acceleration, axle box vertical acceleration is carried out to resonance and demodulation, and the power envelope spectrum that obtains respectively the axle box vertical acceleration corresponding with each frequency in described frequency limit comprises: the energy of supposing axle box vertical acceleration is mainly distributed in [F
_{a}F
_{b}] between Hz, to frequency from F
_{a}hz is to F
_{b}hz scans, and supposes that frequency increment is F
_{c}hz, the sweep rate of the j time is F
_{j}=Fa+ (j1) * FcHz; Repeat the following step Astep C, until complete [F
_{a}F
_{b}] scanning of Hz frequency limit; A, axle box vertical acceleration is carried out to [F
_{l}, F
_{j}] Hz bandpass filtering, generate filtered signal x
_{j}(t); F wherein
_{l}for lowerfrequency limit; B, utilize Hilbert conversion to extract x
_{j}(t) envelope, obtains a
_{j}(t); C, to a
_{j}(t) carrying out refinement Fourier analysis obtains and F
_{j}the power envelope spectrum of corresponding axle box vertical acceleration.
The effective value of the axle box vertical acceleration after described calculation of filtered comprises: suppose that the axle box vertical acceleration signal after bandpass filtering is { x
_{i}, i=1,2 ... N}, N represents sampling number, the window width of supposing effective value is K; The sum of squares of K point before calculating
obtain the 1st effective value:
to j circulation, j=2,3 ... NK+1, calculates all the other effective values:
${\mathrm{RMS}}_{j}=\sqrt{(S+{x}_{j+K1}^{2}{x}_{j1}^{2})/K}.$
According to the effective value of described axle box vertical acceleration, determine the position of welded joint, and the effective value that extracts the axle box vertical acceleration at each welded joint place comprises: a plurality of effective values are divided into one group, calculate aviation value m and the variances sigma of every group of effective value, according to described aviation value m and variances sigma, draw the reference value R of axle box vertical acceleration effective value
_{t}; Axle box vertical acceleration effective value is greater than to reference value R
_{t}the transfinite location positioning at a place of large value be the large value welded joint that transfinites; According to welded joint periodic spacings, determine the transfinite position of all welded joints of welded joint centre of large value; Get near a plurality of large values same welded joint and transfinite the maxim of an effective value as the effective value of the axle box vertical acceleration at this welded joint place.
Before whether the axle box vertical acceleration effective value that judges welded joint place is greater than threshold value, described method also comprises: according to the effective value of the axle box vertical acceleration at described welded joint place, draw the distribution curve of the axle box vertical acceleration effective value at welded joint place; By in the distribution curve of the axle box vertical acceleration effective value at welded joint place, corresponding probability is greater than sets the effective value of threshold value as the threshold value of detection welded joint state.
To achieve these goals, the embodiment of the present invention also provides a kind of checking system of rail in high speed railway welded joint state, described system comprises: acceleration/accel sampling setup, for axle box vertical acceleration is carried out to time sampling, obtains the time sampling signal of axle box vertical acceleration; Band pass filter means, carries out bandpass filtering for the time sampling signal to described axle box vertical acceleration; Effective value computer device, for the effective value of the axle box vertical acceleration after calculation of filtered; Welded joint registration device, for determine the position of welded joint according to the effective value of described axle box vertical acceleration, and extracts the effective value of the axle box vertical acceleration at each welded joint place; Connector status detecting device, for judging whether the axle box vertical acceleration effective value at welded joint place is greater than threshold value, if be greater than threshold value, judges that this welded joint state is bad.
Described system also comprises the logical upper cutoff frequency generating apparatus of band, the logical upper cutoff frequency generating apparatus of described band comprises: resonance and demodulation unit, for in frequency limit corresponding to the power spectrum of axle box vertical acceleration, axle box vertical acceleration is carried out to resonance and demodulation, obtain respectively the power envelope spectrum of the axle box vertical acceleration corresponding with each frequency in described frequency limit; Percentum Drawing of Curve unit, accounts for the percentum of this power envelope spectrum gross energy for calculating the first dominant frequency energy of each power envelope spectrum, draw the relation curve of described percentum and its corresponding frequency; Upper cutoff frequency determining unit, for being defined as frequency corresponding to the peak value of described relation curve percentum the upper cutoff frequency of bandpass filtering.
Described welded joint registration device comprises: reference value generation unit, for a plurality of effective values are divided into one group, calculate aviation value m and the variances sigma of every group of effective value, and according to described aviation value m and variances sigma, draw the reference value R of axle box vertical acceleration effective value
_{t}; The first welded joint determining unit, for being greater than reference value R by axle box vertical acceleration effective value
_{t}the transfinite location positioning at a place of large value be the large value welded joint that transfinites; The second welded joint determining unit, for according to welded joint periodic spacings, determines the transfinite position of all welded joints of welded joint centre of large value; Welded joint effective value determining unit, transfinites the maxim of an effective value as the effective value of the axle box vertical acceleration at this welded joint place for getting near a plurality of large value same welded joint.
Described system also comprises threshold process unit, described threshold process unit specifically comprises: distribution curve drawing unit, for according to the effective value of the axle box vertical acceleration at described welded joint place, draw the distribution curve of the axle box vertical acceleration effective value at welded joint place; Threshold value determining unit, for by the distribution curve of the axle box vertical acceleration effective value at welded joint place, corresponding probability be greater than set threshold value effective value as the threshold value that detects welded joint state.
The technical scheme of the embodiment of the present invention detects welded joint state completely according to the data of axle box vertical acceleration, can effectively detect the current use state of high speed railway welded joint, for reasonable economy instructs the maintenance of welded joint that strong theoretical foundation is provided.
Accompanying drawing explanation
Fig. 1 is the static measurement result schematic diagram of the flatness of prior art steel rail weld joint;
Fig. 2 is the axle box vertical acceleration signal schematic representation of surveying on ripple abrasion rail;
Fig. 3 is the axle box vertical acceleration signal schematic representation on normal rail;
Fig. 4 is the detailed schematic diagram of the method for inspection of embodiment of the present invention high speed railway welded joint state;
Fig. 4 a is the diagram of circuit of the high speed railway welded joint stateevaluation method of the embodiment of the present invention;
Fig. 5 is the power envelope spectrum schematic diagram of the corresponding axle box vertical acceleration of certain frequency of the embodiment of the present invention;
Fig. 6 is that the embodiment of the present invention the first dominant frequency power percentum is with frequency filtering change curve;
Fig. 7 is the axle box acceleration waveform schematic diagram of embodiment of the present invention actual measurement;
Fig. 8 a is one of effective value schematic diagram of the axle box vertical acceleration that calculates of the embodiment of the present invention;
Fig. 8 b be the axle box vertical acceleration that calculates of the embodiment of the present invention effective value schematic diagram two;
Fig. 8 c be the axle box vertical acceleration that calculates of the embodiment of the present invention effective value schematic diagram three;
Fig. 9 is the distribution curve schematic diagram of embodiment of the present invention actual measurement axle box vertical acceleration effective value;
Figure 10 is embodiment of the present invention high speed railway welded joint state final detection result schematic diagram;
Figure 11 is the schematic diagram of the checking system of embodiment of the present invention rail in high speed railway welded joint state;
Figure 11 a is the refinement functional block diagram of the logical upper cutoff frequency generating apparatus 115 of band in Figure 11;
Figure 11 b is the refinement functional block diagram of welded joint registration device 113 in Figure 11;
Figure 11 c is the refinement functional block diagram of threshold process device 116 in Figure 11.
The specific embodiment
The embodiment of the present invention provides a kind of method of inspection and system of high speed railway welded joint state.Wear away and compare with irregular wear, wave and ripple, the Short wave irregularity at welded joint place had both had highfrequency percussion characteristic, had lowfrequency cycle characteristic simultaneously.Utilize axle box vertical acceleration to pass judgment on steel rail weld joint and will run into two crucial problems, i.e. the selection of bandpass filtering frequency and the location of welded joint.Different in concentrating on narrower highfrequency range from undulatory wear, the highfrequency percussion ripple that welded joint causes is distributed in very wide frequency band, and changes.
According to Fourier transform principle, steel rail weld joint momentary impulse is an infinite bandwidth signal from frequency domain, there are in theory all frequency contents, and uniform along whole frequency axis energy, thereby must contain the paxle box system eigentone composition of wheel, also just must cause the characteristic oscillation of system, produce generalized resonance.The embodiment of the present invention, based on abovementioned principle, utilizes resonance and demodulation to propose new high speed railway welded joint condition detection method and system.
The groundwork of axle box vertical acceleration being carried out to resonance and demodulation is as follows:
1, axle box vertical acceleration signal is carried out to bandpass filtering, remember that filtered axle box vertical acceleration continuous signal is x (t), t is sampling instant.
2, utilize Hilbert conversion to extract the envelope of x (t), obtain a (t).
3, a (t) is carried out to refinement Fourier analysis, obtain the refinement power envelope spectrum of axle box vertical acceleration.
The dynamic testing method of embodiment of the present invention steel rail weld joint mainly comprises following algorithm: the axle box vertical acceleration based on resonance and demodulation moves the fast algorithm of effective value; Axle box vertical acceleration bandpass filtering upper cutoff frequency self adaptation Algorithms of Selecting based on resonance and demodulation; And the accurate location of the welded joint based on resonance and demodulation and statedetection algorithm.
For making object, technical scheme and the advantage of the embodiment of the present invention clearer, below in conjunction with the accompanying drawing in the embodiment of the present invention, technical scheme in the embodiment of the present invention is clearly and completely described, obviously, described embodiment is the present invention's part embodiment, rather than whole embodiment.Embodiment based in the present invention, those of ordinary skills are not making the every other embodiment obtaining under creative work prerequisite, all belong to the scope of the claims in the present invention protection.
Fig. 4 is the detailed schematic diagram of the method for inspection of embodiment of the present invention high speed railway welded joint state, as shown in Figure 4, the method comprises two steps: the 1st step (lefthalf of Fig. 4) utilizes a large amount of actual measurement axle box vertical acceleration data to calculate the effective value at welded joint place, and according to effective value distribution character, calculate threshold value; Particularly, Fig. 7 is the axle box acceleration waveform schematic diagram of actual measurement; Wherein lefthalf is axle box vertical acceleration signal, and right half part is axle box transverse acceleration signal, can find out, when transverse acceleration variable is large, also phase strain is large for vertical acceleration; Fig. 8 aFig. 8 c is the effective value schematic diagram of the axle box vertical acceleration that calculates; Fig. 9 is the distribution curve schematic diagram of axle box vertical acceleration effective value; The 2nd step (right half part of Fig. 4), according to vertical acceleration effective value and the threshold value at welded joint place, detects the state of welded joint online.
The key algorithm of the 1st step and the 2nd step is identical.Key algorithm comprises: axle box vertical acceleration moves fast algorithm, the axle box vertical acceleration bandpass filtering upper cutoff frequency self adaptation Algorithms of Selecting of effective value, the accurate location algorithm of welded joint.In addition, the 1st step also comprises the algorithm of the axle box vertical acceleration effective value distribution of calculating welded joint place, the embodiment of the present invention adopts histogram table to show the probability distribution function of effective value, then it is normalized to the distribution function that obtains effective value, concrete minute 3 steps complete: the 1st step calculates the axle box vertical acceleration effective value at welded joint place, the 2nd step calculates the axle box vertical acceleration effective value histogram at welded joint place, i.e. probability distribution function; To histogram, normalization method obtains its distribution function to the 3rd step.
Fig. 4 a is the diagram of circuit of the high speed railway welded joint stateevaluation method of the embodiment of the present invention, and as shown in Fig. 4 a, the method comprises:
S401, axle box vertical acceleration is carried out to time sampling, obtain the time sampling signal of axle box vertical acceleration.
The collecting device of axle box vertical acceleration is arranged on the axle box of comprehensive detection train, and train, in operational process, gathers the axle box vertical acceleration of a plurality of positions on rail, to the sampling of the vertical acceleration of axle box, is therefore a time sampling.
S402, the time sampling signal of axle box vertical acceleration is carried out to bandpass filtering.
The frequency limit upper limit of bandpass filtering will obtain by described bandpass filtering upper cutoff frequency self adaptation Algorithms of Selecting below, because axle box vertical acceleration is mainly radiofrequency component, lowfrequency component proportion seldom, therefore the selection of lower limit on the result impact of filtering seldom, can preset a lower frequency limit F according to actual needs
_{l}, as 20Hz.
Suppose the upper cutoff frequency F that the embodiment of the present invention adopts bandpass filtering upper cutoff frequency self adaptation Algorithms of Selecting to obtain
_{h}450Hz, with F
_{l}=20HZ is example, and the frequency limit of bandpass filtering is [20,450] Hz.It should be noted that; this concrete frequency limit is only for explaining the embodiment of the present invention; and be not used in, the protection domain of claim is limited, what this case claim needed protection is the method that how to obtain this frequency limit, but not concrete frequency limit numerical value.
The axle box vertical acceleration bandpass filtering upper cutoff frequency self adaptation Algorithms of Selecting that the embodiment of the present invention adopts comprises: first, in frequency limit corresponding to the power spectrum of axle box vertical acceleration, axle box vertical acceleration is carried out to resonance and demodulation, obtain respectively the power envelope spectrum of the axle box vertical acceleration corresponding with each frequency in described frequency limit; Then, calculate the percentum that the first dominant frequency energy in each power envelope spectrum accounts for this power envelope spectrum gross energy, draw the relation curve of described percentum and its corresponding frequency; Finally, the frequency corresponding to peak value of percentum in described relation curve is defined as to bandpass filtering cutoff frequency.
Concrete algorithm is as follows:
(1) in frequency limit corresponding to the power spectrum of axle box vertical acceleration, axle box vertical acceleration is carried out to resonance and demodulation, obtain respectively the power envelope spectrum of the axle box vertical acceleration corresponding with each frequency in described frequency limit.
The energy of supposing axle box vertical acceleration is mainly distributed in [F
_{a}F
_{b}] between Hz, to frequency from F
_{a}hz is to F
_{b}hz scans, and supposes that frequency increment is F
_{c}hz, the sweep rate of the j time is F
_{j}=Fa+ (j1) * FcHz; Repeat the following step Astep C, until complete [F
_{a}F
_{b}] scanning of Hz frequency limit.
A, axle box vertical acceleration is carried out to [F
_{l}, F
_{j}] Hz bandpass filtering, generate filtered signal x
_{j}(t); F wherein
_{l}for lowerfrequency limit, as 20Hz;
B, utilize Hilbert conversion to extract x
_{j}(t) envelope, obtains a
_{j}(t);
C, to a
_{j}(t) carry out refinement Fourier analysis, obtain and F
_{j}the power envelope spectrum of the axle box vertical acceleration that Hz is corresponding.
In the present embodiment, by known to the power spectral analysis of axle box vertical acceleration, the energy of axle box vertical acceleration is mainly distributed between [100600] Hz, therefore frequency is scanned from 100Hz to 600Hz, suppose that frequency increment is 10Hz, the sweep rate of the j time is F
_{j}=100+ (j1) * 10Hz; This frequency increment can also change according to realistic accuracy demand.
Steps A, B, C are referred to as the process of resonance and demodulation, by this resonance and demodulation process just can access and frequency limit [100600] Hz of the power spectrum of axle box vertical acceleration between the corresponding a plurality of power envelopes spectrums of each frequency, Fig. 5 be that the power envelope of the corresponding axle box vertical acceleration of certain frequency is composed schematic diagram.
(2) calculate the percentum that the first dominant frequency energy in each power envelope spectrum accounts for this power envelope spectrum gross energy, draw the relation curve of described percentum and its corresponding frequency.
D, from Fig. 5, therefore the energy of front 5Hz can represent gross energy, calculates respectively the ENERGY E 1 of the first dominant frequency of envelope spectrum and the ENERGY E in front 5Hz;
E, calculating the first dominant frequency energy account for gross energy percentum
Step D, E are also the processes of a circulation, the corresponding power envelope of each frequency spectrum need to be carried out one time, draw the percentum of each power envelope spectrum, draw the relation curve of described percentum and corresponding frequency, as shown in Figure 6.
(3) frequency corresponding to peak value of percentum in described relation curve is defined as to bandpass filtering upper cutoff frequency.
Peak value in Fig. 6 is 450HZ, so obtain the upper limit of frequency filtering scope, is 450Hz.
The effective value of the axle box vertical acceleration after S403, calculation of filtered.
The embodiment of the present invention is when calculating axle box vertical acceleration moves effective value, by time sampling.Suppose that the axle box vertical acceleration signal after bandpass filtering is { x
_{i}, i=1,2 ... N}, wherein N represents sampling number, and the window width of supposing effective value is K, and the algorithm of the mobile effective value of axle box vertical acceleration is as follows:
(1) calculate the sum of squares of front K point
obtain the 1st effective value:
(2) to j circulation, j=2,3 ... NK+1, calculates all the other effective values:
Casehistory is as follows: suppose sampling number N=11, and the long K=7 of window, the axle box vertical acceleration signal after bandpass filtering is { x
_{i}, i=1,2 ... 11}, the 1st definition to NK+1=5 effective value and being calculated as follows:
By calculated amount, analyze knownly, be linear the computing time of the algorithm of the mobile effective value of abovementioned calculating axle box vertical acceleration, and computation complexity is O (N).
S404, according to the effective value of axle box vertical acceleration, determine the position of welded joint, and extract the effective value of the axle box vertical acceleration at each welded joint place.
Specifically comprise: a plurality of effective values are divided into one group, calculate aviation value m and the variances sigma of every group of effective value, according to described aviation value m and variances sigma, draw the reference value R of axle box vertical acceleration effective value
_{t}; Axle box vertical acceleration effective value is greater than to reference value R
_{t}the transfinite location positioning at a place of large value be the large value welded joint that transfinites; According to welded joint periodic spacings, determine the transfinite position of all welded joints of welded joint centre of large value; Get near a plurality of large values same welded joint and transfinite the maxim of an effective value as the effective value of the axle box vertical acceleration at this welded joint place.
Welded joint Precision Orientation Algorithm is as follows:
(1) a plurality of effective values are divided into one group, calculate aviation value m and the variances sigma of the effective value of each group, as, the reference value R of the effective value of axle box vertical acceleration can be calculated by formula (3)
_{t}, the reference value obtaining, as shown in the dotted line in Fig. 8 aFig. 8 c, has only demonstrated the reference value of a group in Fig. 8 aFig. 8 c;
R
_{T}＝m+2.5σ (3)
(2) find out effective value and be greater than reference value R
_{t}point, being defined as large value transfinites a little, the transfinite position at a place of these large values is exactly welded joint, is designated as
result is as being greater than in Fig. 8 a as shown in zero ' o ' point; Because same welded joint has a plurality of large values to transfinite a little around, so there is a plurality of ' o ' point;
(3) right
carry out polymerization, the large value at same welded joint place is transfinited and to retain a maxim, with ' o ' point that guarantees that same welded joint only has 1 large value to transfinite, the welded joint that transfinites of the large value after polymerization is designated as
result is as being greater than in Fig. 8 b as shown in zero ' o ' point;
(4) according to welded joint periodic spacings, be approximately 100m, find out the welded joint that large value transfinites in the middle of welded joint, finally obtain all welded joints, and axle box vertical acceleration effective value corresponding to all welded joints, as being greater than in Fig. 8 c as shown in zero ' o ' point.
If the axle box vertical acceleration effective value at S405 welded joint place is greater than threshold value, judge that this welded joint state is bad.
The mode of definite threshold is as follows: according to the effective value of the axle box vertical acceleration at described welded joint place, draw the distribution curve of the axle box vertical acceleration effective value at welded joint place; By in the distribution curve of the axle box vertical acceleration effective value at welded joint place, corresponding probability is greater than sets the effective value of threshold value as the threshold value of detection welded joint state.
Take Fig. 9 as example, and in the distribution curve of the axle box vertical acceleration effective value at the welded joint place shown in Fig. 9, hypothetical probabilities threshold value is 99.5%, and corresponding probability is that 99.5% effective value 50m/s/s is exactly the threshold value of passing judgment on welded joint state.As shown in figure 10, the effective value of the welded joint that state is bad is apparently higher than the effective value at other welded joint place for final detection result.
Corresponding to the method for Fig. 4 a, the embodiment of the present invention also provides a kind of high speed railway welded joint stateevaluation system.Figure 11 is the schematic diagram of the checking system of embodiment of the present invention rail in high speed railway welded joint state, and as shown in figure 11, this system comprises:
Acceleration/accel sampling setup 110, for axle box vertical acceleration is carried out to time sampling, obtains the time sampling signal of axle box vertical acceleration; Band pass filter means 111, carries out bandpass filtering for the time sampling signal to described axle box vertical acceleration; Effective value computer device 112, for the effective value of the axle box vertical acceleration after calculation of filtered; Welded joint registration device 113, for determine the position of welded joint according to the effective value of described axle box vertical acceleration, and extracts the effective value of the axle box vertical acceleration at each welded joint place; Connector status detecting device 114, for judging whether the axle box vertical acceleration effective value at welded joint place is greater than threshold value, if be greater than threshold value, judges that this welded joint state is bad.
As shown in figure 11, this system also comprises the logical upper cutoff frequency generating apparatus 115 of band, as shown in Figure 11 a, the logical upper cutoff frequency generating apparatus of this band comprises: resonance and demodulation unit 1150, for in frequency limit corresponding to the power spectrum of axle box vertical acceleration, axle box vertical acceleration is carried out to resonance and demodulation, obtain respectively the power envelope spectrum of the axle box vertical acceleration corresponding with each frequency in described frequency limit; Percentum Drawing of Curve unit 1151, accounts for the percentum of this power envelope spectrum gross energy for calculating the first dominant frequency energy of each power envelope spectrum, draw the relation curve of described percentum and its corresponding frequency; Upper cutoff frequency determining unit 1152, for being defined as frequency corresponding to the peak value of described relation curve percentum the upper cutoff frequency of bandpass filtering.
Figure 11 b is the refinement functional block diagram of welded joint registration device 113, as shown in Figure 11 b, welded joint registration device 113 specifically comprises: reference value generation unit 1130, for a plurality of effective values are divided into one group, calculate aviation value m and the variances sigma of every group of effective value, according to described aviation value m and variances sigma, draw the reference value R of axle box vertical acceleration effective value
_{t}; The first welded joint determining unit 1131, for being greater than reference value R by axle box vertical acceleration effective value
_{t}the transfinite location positioning at a place of large value be the large value welded joint that transfinites; The second welded joint determining unit 1132, for according to welded joint periodic spacings, determines the transfinite position of all welded joints of welded joint centre of large value; Welded joint effective value determining unit 1133, transfinites the maxim of an effective value as the effective value of the axle box vertical acceleration at this welded joint place for getting near a plurality of large value same welded joint.
As shown in figure 11, this system also comprises threshold process device 116; As shown in Figure 11 c, threshold process device 116 specifically comprises: distribution curve drawing unit 1160, for according to the effective value of the axle box vertical acceleration at described welded joint place, draw the distribution curve of the axle box vertical acceleration effective value at welded joint place; Threshold value determining unit 1161, for by the distribution curve of the axle box vertical acceleration effective value at welded joint place, corresponding probability be greater than set threshold value effective value as the threshold value that detects welded joint state.
The beneficial effect that technical solution of the present invention is brought: by the indepth analysis of the axle box vertical acceleration data to more than 1000 kilometer of actual measurement of certain line, binding resonant demodulation techniques calculate the distribution function curve of welded joint place effective value can be found out, the distribution curve that under identical operating mode, twice collection calculates coincide consistent.Therefore, new welded joint Precision Orientation Algorithm is stable, and the distribution curve obtaining has good repeatability.And, utilize new method of inspection to pass judgment on welded joint state completely, can effectively detect high speed railway welded joint state, for reasonable economy instructs the maintenance of welded joint that strong theoretical foundation is provided.
Above embodiment only, in order to the technical scheme of the embodiment of the present invention to be described, is not intended to limit; Although the embodiment of the present invention is had been described in detail with reference to previous embodiment, those of ordinary skill in the art is to be understood that: its technical scheme that still can record aforementioned each embodiment is modified, or part technical characterictic is wherein equal to replacement; And these modifications or replacement do not make the essence of appropriate technical solution depart from the spirit and scope of each embodiment technical scheme of the embodiment of the present invention.
List of references:
[1]C.Esweld.Modern Railway Track[M].Duisburg：MRTproductions，2001.
[2]M.Bocciolone，A.
，A.Cigada，A.Collina.A Measurement System for Quick Rail Inspection and Effective Track Maintenance Strategy[J].Mechanical Systems and Signal Processing，2007(21)：12421254.
[3]Sunaga Y.，T.Ide.and M.Kanao.A practical use of axlebox acceleration to control the short wave track irregularities on SHINKANSEN[J].RTRI quarterly Report 1995，9(2)：3540.
[4]Sunaga Y.，Snal，I.and Ide，T..A method to control the short wave track irregularities utilizing axlebox acceleration[J].RTRI quarterly Report，1997，38(4)：176181.
[5]Yoshida，M.，Shinowaki，S.and Sunaga，Y..A study on lateral axlebox acceleration for detecting irregularities of rail joints on sharp curves[J].Jointed Railway Technology Symposium(JRail2004)，2004，11：9154.
[6]Hirofumi TANAKA，Atsushi FURUKAWA.Management Method for Large Wheel Loads and Lateral Forces Using Axlebox Acceleration and Track Irregularity[J].RTRI quarterly Report，2009，50(2)：116122.
[7]Eric G.Berggren，Martin X.D.Li，Jan Spannar.A new approach to the analysis and presentation of vertical track geometry quality and rail roughness[J].Wear 2008(265)：14881496.
[8]A Bracciali and G Cascini.Detection of corrugation and wheelflats of railway wheels using energy and cepstrum analysis of rail acceleration[J].Proc Instn Mech Engrs，1997，211，(Part E)：109116.
[9]Grassie，S.L.Rail corrugation：advances in measurement，understanding and treatment[J].Wear，2005，258：12241234.
[10]Hempelmann，K.and Knothe，K.An extended linear model for the prediction of short pitch corrugation[J].Wear，1996，191：161169.
Claims (9)
1. a method of inspection for rail in high speed railway welded joint state, is characterized in that, described method comprises:
Axle box vertical acceleration is carried out to time sampling, obtain the time sampling signal of axle box vertical acceleration;
Time sampling signal to described axle box vertical acceleration carries out bandpass filtering;
The effective value of the axle box vertical acceleration after calculation of filtered;
According to the effective value of described axle box vertical acceleration, determine the position of welded joint, and extract the effective value of the axle box vertical acceleration at each welded joint place;
Whether the axle box vertical acceleration effective value that judges welded joint place is greater than threshold value, if be greater than threshold value, judges that this welded joint state is bad;
The effective value of the axle box vertical acceleration after described calculation of filtered comprises:
Suppose that the axle box vertical acceleration signal after bandpass filtering is { x
_{i}, i=1,2 ... N}, N represents sampling number, the window width of supposing effective value is K; The sum of squares of K point before calculating
obtain the 1st effective value:
${\mathrm{RMS}}_{1}=\sqrt{S/K};$
To j circulation, j=2,3 ... NK+1, calculates all the other effective values:
2. method according to claim 1, is characterized in that, before the time sampling signal to described axle box vertical acceleration carries out bandpass filtering, described method also comprises: the lowerfrequency limit F that presets bandpass filtering
_{l}, and the upper cutoff frequency F that dynamically determines bandpass filtering
_{h}; Wherein, dynamically determine the upper cutoff frequency F of bandpass filtering
_{h}specifically comprise:
In frequency limit corresponding to the power spectrum of axle box vertical acceleration, axle box vertical acceleration is carried out to resonance and demodulation, obtain respectively the power envelope spectrum of the axle box vertical acceleration corresponding with each frequency in described frequency limit;
Calculate the percentum that the first dominant frequency energy in each power envelope spectrum accounts for this power envelope spectrum gross energy, draw the relation curve of described percentum and its corresponding frequency;
The frequency corresponding to peak value of percentum in described relation curve is defined as to the upper cutoff frequency of bandpass filtering.
3. method according to claim 2, it is characterized in that, described in frequency limit corresponding to the power spectrum of axle box vertical acceleration, axle box vertical acceleration is carried out to resonance and demodulation, and the power envelope spectrum that obtains respectively the axle box vertical acceleration corresponding with each frequency in described frequency limit comprises:
The energy of supposing axle box vertical acceleration is mainly distributed in [F
_{a}F
_{b}] between Hz, to frequency from F
_{a}hz is to F
_{b}hz scans, and supposes that frequency increment is F
_{c}hz, the sweep rate of the j time is F
_{j}=Fa+ (j1) * FcHz; Repeat the following step Astep C, until complete [F
_{a}F
_{b}] scanning of Hz frequency limit;
A, axle box vertical acceleration is carried out to [FL, F
_{j}] Hz bandpass filtering, generate filtered signal x
_{j}(t); F wherein
_{l}for lowerfrequency limit;
B, utilize Hilbert conversion to extract x
_{j}(t) envelope, obtains a
_{j}(t);
C, to a
_{j}(t) carry out refinement Fourier analysis, obtain and F
_{j}the power envelope spectrum of the axle box vertical acceleration that Hz is corresponding.
4. method according to claim 1, is characterized in that, determines the position of welded joint according to the effective value of described axle box vertical acceleration, and the effective value that extracts the axle box vertical acceleration at each welded joint place comprises:
A plurality of effective values are divided into one group, calculate aviation value m and the variances sigma of every group of effective value, according to described aviation value m and variances sigma, draw the reference value R of axle box vertical acceleration effective value
_{t};
Axle box vertical acceleration effective value is greater than to reference value R
_{t}the transfinite location positioning at a place of large value be the large value welded joint that transfinites;
According to welded joint periodic spacings, determine the transfinite position of all welded joints of welded joint centre of large value;
Get near a plurality of large values same welded joint and transfinite the maxim of an effective value as the effective value of the axle box vertical acceleration at this welded joint place.
5. method according to claim 1, is characterized in that, before whether the axle box vertical acceleration effective value that judges welded joint place is greater than threshold value, described method also comprises:
According to the effective value of the axle box vertical acceleration at described welded joint place, draw the distribution curve of the axle box vertical acceleration effective value at welded joint place;
By in the distribution curve of the axle box vertical acceleration effective value at welded joint place, corresponding probability is greater than sets the effective value of threshold value as the threshold value of detection welded joint state.
6. a checking system for rail in high speed railway welded joint state, is characterized in that, described system comprises:
Acceleration/accel sampling setup, for axle box vertical acceleration is carried out to time sampling, obtains the time sampling signal of axle box vertical acceleration;
Band pass filter means, carries out bandpass filtering for the time sampling signal to described axle box vertical acceleration;
Effective value computer device, the effective value for the axle box vertical acceleration after calculation of filtered, specifically comprises: suppose that the axle box vertical acceleration signal after bandpass filtering is { x
_{i}, i=1,2 ... N}, N represents sampling number, the window width of supposing effective value is K; The sum of squares of K point before calculating
obtain the 1st effective value:
to j circulation, j=2,3 ... NK+1, calculates all the other effective values:
${\mathrm{RMS}}_{j}=\sqrt{(S+{x}_{j+K1}^{2}{x}_{j1}^{2})/K};$
Welded joint registration device, for determine the position of welded joint according to the effective value of described axle box vertical acceleration, and extracts the effective value of the axle box vertical acceleration at each welded joint place;
Connector status detecting device, for judging whether the axle box vertical acceleration effective value at welded joint place is greater than threshold value, if be greater than threshold value, judges that this welded joint state is bad.
7. system according to claim 6, is characterized in that, described system also comprises the logical upper cutoff frequency generating apparatus of band, and the logical upper cutoff frequency generating apparatus of described band comprises:
Resonance and demodulation unit, in frequency limit corresponding to the power spectrum in axle box vertical acceleration, carries out resonance and demodulation to axle box vertical acceleration, obtains respectively the power envelope spectrum of the axle box vertical acceleration corresponding with each frequency in described frequency limit;
Percentum Drawing of Curve unit, accounts for the percentum of this power envelope spectrum gross energy for calculating the first dominant frequency energy of each power envelope spectrum, draw the relation curve of described percentum and its corresponding frequency;
Upper cutoff frequency determining unit, for being defined as frequency corresponding to the peak value of described relation curve percentum the upper cutoff frequency of bandpass filtering.
8. system according to claim 6, is characterized in that, described welded joint registration device comprises:
Reference value generation unit, for a plurality of effective values are divided into one group, calculates aviation value m and the variances sigma of every group of effective value, draws the reference value R of axle box vertical acceleration effective value according to described aviation value m and variances sigma
_{t};
The first welded joint determining unit, for being greater than reference value R by axle box vertical acceleration effective value
_{t}the transfinite location positioning at a place of large value be the large value welded joint that transfinites;
The second welded joint determining unit, for according to welded joint periodic spacings, determines the transfinite position of all welded joints of welded joint centre of large value;
Welded joint effective value determining unit, transfinites the maxim of an effective value as the effective value of the axle box vertical acceleration at this welded joint place for getting near a plurality of large value same welded joint.
9. system according to claim 6, is characterized in that, described system also comprises threshold process unit, and described threshold process unit specifically comprises:
Distribution curve drawing unit, for according to the effective value of the axle box vertical acceleration at described welded joint place, draws the distribution curve of the axle box vertical acceleration effective value at welded joint place;
Threshold value determining unit, for by the distribution curve of the axle box vertical acceleration effective value at welded joint place, corresponding probability be greater than set threshold value effective value as the threshold value that detects welded joint state.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201110344671.6A CN102501886B (en)  20111104  20111104  Method and system for detecting states of steel rail welding joints of highspeed railway 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201110344671.6A CN102501886B (en)  20111104  20111104  Method and system for detecting states of steel rail welding joints of highspeed railway 
Publications (2)
Publication Number  Publication Date 

CN102501886A CN102501886A (en)  20120620 
CN102501886B true CN102501886B (en)  20140813 
Family
ID=46214042
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201110344671.6A Active CN102501886B (en)  20111104  20111104  Method and system for detecting states of steel rail welding joints of highspeed railway 
Country Status (1)
Country  Link 

CN (1)  CN102501886B (en) 
Families Citing this family (5)
Publication number  Priority date  Publication date  Assignee  Title 

CN102890143B (en) *  20121019  20150715  西南交通大学  Rail local defect vehiclemounting detection method merging with vehicle speed information and front and rear axle box acceleration information 
CN105509668B (en) *  20151127  20181120  北京锦鸿希电信息技术股份有限公司  Track wave grinds detection system 
CN105346561A (en) *  20151202  20160224  北京交通大学  Rail turnout disease detection system based on operating vehicle and rail turnout disease detection method based on operating vehicle 
CN108318169A (en) *  20180112  20180724  西南交通大学  Maximum dynamic force appraisal procedure and steel rail welding line maintenance system at a kind of steel rail welding line 
CN113501028B (en) *  20210707  20220809  中国铁道科学研究院集团有限公司  Method and device for diagnosing poor welded joint of heavyduty railway steel rail 
Citations (1)
Publication number  Priority date  Publication date  Assignee  Title 

CN1369595A (en) *  20020315  20020918  西南交通大学  Tester of rail state 
Family Cites Families (1)
Publication number  Priority date  Publication date  Assignee  Title 

US6763312B1 (en) *  20030111  20040713  Dynamic Measurement Consultants, Llc  Multiple discriminate analysis and data integration of vibration in rotation machinery 

2011
 20111104 CN CN201110344671.6A patent/CN102501886B/en active Active
Patent Citations (1)
Publication number  Priority date  Publication date  Assignee  Title 

CN1369595A (en) *  20020315  20020918  西南交通大学  Tester of rail state 
NonPatent Citations (1)
Title 

M．Molodova et al.An investigation of the possibility to use axle box acceleration for condition monitoring of welds.《international conference of noise and vibration engineering》.2008,图8、文章第4.1及第5部分. * 
Also Published As
Publication number  Publication date 

CN102501886A (en)  20120620 
Similar Documents
Publication  Publication Date  Title 

CN108845028B (en)  Method and device for dynamically detecting highspeed railway rail corrugation  
CN102501886B (en)  Method and system for detecting states of steel rail welding joints of highspeed railway  
CN108515984B (en)  Wheel damage detection method and device  
CN102114855B (en)  Track detection method and device  
CN106394606B (en)  A kind of railway car wheel loses circle detection method and detection device  
CN202115547U (en)  Track geometry status measuring cart and measuring apparatus  
CN104359444B (en)  Rail outline detection method and system  
CN102874278B (en)  Vehiclemounted detection method fusing vehicle speed information and axle box vertical acceleration information for wheel flats  
CN105292177A (en)  Method for measuring track corrugation by utilizing axle box vibration and impact information  
CN110789566B (en)  Track defect monitoring method and monitoring equipment based on axle box acceleration signal  
CN102706572A (en)  Fault diagnosis and rehabilitation center for train wheel sets  
CN202644332U (en)  Rail top surface irregularity detecting device  
CN106018557A (en)  Track surface indentation detection method based on wavelet analysis  
CN102890143B (en)  Rail local defect vehiclemounting detection method merging with vehicle speed information and front and rear axle box acceleration information  
CN206307052U (en)  A kind of railway car wheel loses loop truss device  
CN103031787A (en)  Method and device for fully automatically and accurately monitoring changes of geometric parameters of track  
CN106092015A (en)  A kind of raceway surface depression length detecting method  
CN111332338A (en)  Track disease detection system  
CN104792275A (en)  Tunnel deformation detection unit and method  
CN104949848A (en)  Hunting instability detecting algorithm based on frequency domain kurtosis  
CN101701976A (en)  Third rail detection device and method  
CN105128888B (en)  Track waveform abrasion method for quick based on many midpoint string synchro measures  
CN204775291U (en)  Device is exerted to track geometry not smooth going detecting system's load based on initiative excitation  
CN106052606A (en)  Track surface recess detection method based on scale average wavelet energy spectrum  
CN108016990A (en)  The detection device and detection method that a kind of portal crane trolley track fastener loosens 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
C14  Grant of patent or utility model  
GR01  Patent grant 