CN102914364B - Dynamic identifying method for instant contact loss of wheel rail under impact effect of wheel rail - Google Patents

Dynamic identifying method for instant contact loss of wheel rail under impact effect of wheel rail Download PDF

Info

Publication number
CN102914364B
CN102914364B CN201210399620.8A CN201210399620A CN102914364B CN 102914364 B CN102914364 B CN 102914364B CN 201210399620 A CN201210399620 A CN 201210399620A CN 102914364 B CN102914364 B CN 102914364B
Authority
CN
China
Prior art keywords
frequency
time
vibration acceleration
wheel
wheel rail
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201210399620.8A
Other languages
Chinese (zh)
Other versions
CN102914364A (en
Inventor
林建辉
丁建明
伍川辉
刘璐
易彩
陈春俊
张兵
苏燕辰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201210399620.8A priority Critical patent/CN102914364B/en
Publication of CN102914364A publication Critical patent/CN102914364A/en
Application granted granted Critical
Publication of CN102914364B publication Critical patent/CN102914364B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

The invention discloses a dynamic identifying method for instant contact loss of a wheel rail under an impact effect of the wheel rail and belongs to a traffic information safety technology. The dynamic identifying method comprises the following steps: firstly, mounting a vibration acceleration sensor above a wheel axle box; collecting axle box vibration acceleration information measured by the vibration acceleration sensor through a background collector; extracting a time-frequency feature for the axle box vibration acceleration, thereby obtaining a time-frequency distribution feature; and drawing a contour map for the time-frequency distribution feature and judging if the wheel rail instantly loses the contact. According to the scheme, the purpose of immediately detecting the ultra-large hidden danger of the instant contact loss of the wheel rail under the impact effect of the wheel rail is achieved according to the dynamic identifying method provided by the invention; the possibility of derailing of a vehicle is eliminated; the technology is mature; the reliability is high; the cost is low; the quantity of the channels is less; the system is simple; the anti-interference property is high; and the social and economical benefits are very obvious.

Description

A kind ofly because causing wheel track moment, wheel-rail impact effect loses the dynamic identifying method of contact
Technical field
The invention belongs to traffic and transport field, specifically, relate to and a kind ofly because wheel-rail impact effect causes wheel track moment, lose the dynamic identifying method of contact.
background technology
Rolling stock is under arms in process, and during wheel brake hard locking, if wheel exists at Rail Surface while moving on the soft objects such as ice and snow, leaf, rotation of wheel motion is become and slided or the motion of trackslipping by rolling, and wheel track will produce sliding friction.For relative rolling friction, sliding friction is violent, and violent sliding friction produces that high temperature makes that wheel tread occurs abrading, peels off, rhegma and fall the wheel faults such as piece, and these defects are brought strike and the injury of intermittence to wheel rail system.Existing theoretical calculation and test test shows, having the wheel-rail impact power of the wheel generation of defect is the several times of wheel track static load, wherein HFS directly has a negative impact to track and wheel, thereby reduce serviceable life, the increase energy consumption of vehicle and circuit and improve maintenance cost, serious is to cause axle temperature to raise to occur that axletree is fervent; Accelerate on the other hand axletree, tread Crack Extension causes axletree cold cut.What pay particular attention to is simultaneously, when rolling stock travelling speed reaches certain value, when flaw size reaches certain value, wheel-rail impact can cause that later wheel track moment loses contact, some meetings make wheel track repeatedly intermittence lose contact, wheel track loses contact and means that the vertical acting force of wheel track is zero, moment derailment coefficients be infinity, as there is the section of large transverse force at bend circuit, track Ultrahigh Bypass etc. in fast running, the derailing that vehicle occurs may be very large.
summary of the invention
The object of the invention is to lose the especially big hidden danger of contact because wheel-rail impact causes wheel track moment, provide a kind of and lose the dynamic identifying method of contact because wheel-rail impact effect causes wheel track moment, the possibility of avoiding vehicle that derailing occurs occurs.
To achieve these goals, the technical solution used in the present invention is as follows:
Because wheel-rail impact effect causes wheel track moment, lose a dynamic identifying method for contact, realize by the following method:
(1) vibration acceleration sensor is installed on to wheel axle box top, and gathers by the collector on backstage the axle box vibration acceleration information that vibration acceleration sensor measures;
(2) axle box vibration acceleration information is carried out to time-frequency characteristics extraction, obtain time-frequency distributions feature;
(3) draw the contour map of time-frequency distributions feature, judge whether moment loses contact to wheel track.
In order to guarantee the accurate measurement of vibration acceleration sensor, in described step (1), the range of vibration acceleration sensor is ± 50g that frequency response is not less than 2kHz, and the sample frequency of collector is not less than 10kHz.
Further obtain video distribution feature, the specific implementation step of described step (2) is as follows:
(2a) select frequency slice function, wherein frequency coordinate , the centre frequency that represents conversion,
Figure DEST_PATH_387190DEST_PATH_IMAGE003
represent frequency resolution,
Figure DEST_PATH_80339DEST_PATH_IMAGE004
expression departs from the number of the resolution of centre frequency, so the basis function of frequency choice function is:
Figure DEST_PATH_735312DEST_PATH_IMAGE005
According to
Figure DEST_PATH_19662DEST_PATH_IMAGE006
size, frequency choice function calculates the following three kinds of situations that comprise:
Situation one: when
Figure DEST_PATH_901031DEST_PATH_IMAGE007
and
Figure DEST_PATH_522505DEST_PATH_IMAGE008
time, frequency choice function for:
Figure DEST_PATH_362602DEST_PATH_IMAGE010
Figure DEST_PATH_262425DEST_PATH_IMAGE011
and
Figure DEST_PATH_422011DEST_PATH_IMAGE012
;
Situation two: when
Figure DEST_PATH_723679DEST_PATH_IMAGE008
and time, frequency choice function for:
Figure DEST_PATH_903491DEST_PATH_IMAGE014
Figure DEST_PATH_262928DEST_PATH_IMAGE015
and ;
Situation three: when
Figure DEST_PATH_262294DEST_PATH_IMAGE017
and
Figure DEST_PATH_904628DEST_PATH_IMAGE018
time, frequency choice function
Figure DEST_PATH_915309DEST_PATH_IMAGE009
for:
Figure DEST_PATH_476740DEST_PATH_IMAGE019
Figure DEST_PATH_572872DEST_PATH_IMAGE020
and ;
Wherein,
Figure DEST_PATH_618506DEST_PATH_IMAGE022
for the frequency discrimination factor,
Figure DEST_PATH_350839DEST_PATH_IMAGE006
centered by coordinate corresponding to frequency, meet
Figure DEST_PATH_403108DEST_PATH_IMAGE023
;
(2b) axle box vibration acceleration is carried out to Fourier transform, transformation for mula is as follows:
Figure DEST_PATH_652824DEST_PATH_IMAGE024
(2c) frequency slice function is selected the frequency of Fourier, selects formula as follows:
Figure DEST_PATH_231573DEST_PATH_IMAGE025
(2d) for various centre frequencies, carry out the time-frequency distributions feature that inverse-Fourier transform obtains vibration acceleration:
Figure DEST_PATH_541331DEST_PATH_IMAGE026
In order better to realize the present invention, the concrete methods of realizing of described step (3) is as follows:
(3a) draw the contour map of the time-frequency distributions feature of vibration acceleration;
(3b) find time period corresponding to this contour map maximum energy peak value;
(3c) at ceiling capacity, in the time period, find successively second largest energy peak, if there is the decay that local energy increases with frequency in section at one time, exist wheel track that the phenomenon departing from occurs.
The present invention compares compared with prior art, has the following advantages and beneficial effect:
(1) the invention solves a most important problem of present stage, how to judge because wheel-rail impact effect causes its moment and lose contact, avoided a kind of the most severe factor that endangers vehicle operating safety of bringing because of wheel fault, due to the method that can not address the aforementioned drawbacks at present, therefore the present invention is as an important breakthrough in this area, thoroughly changed this sky window, because wheel track loses and contacts with wheel suddenly, the harm bringing is immeasurable, therefore the present invention is a new technology, has thoroughly solved this harm.
(2) recognition methods of the present invention is very simple, only need vibration acceleration sensor and the collector of a passage to collect information needed, and by calculating, feature extraction, relatively can recognize immediately the phenomenon that this instantaneous wheel and wheel track depart from, there is real-time, and process is very simple; Simultaneously the present invention does not need to transform circuit, can carry out system-wide and cover and detect, and reliability is high, easily fill easy care, can meet the needs that wheel moment loses the Dynamic Recognition of contact completely.
(3) the present invention, owing to being a brand-new technology, has advance and superiority, and technology maturation, reliability are high, have good prospect and application market.
accompanying drawing explanation
Fig. 1 is realization flow figure of the present invention.
Fig. 2 is vibration acceleration sensor installation site schematic diagram in the present invention-embodiment.
Fig. 3 is acceleration and the time-frequency characteristics figure thereof when picking rate is 150km/h in the present invention-embodiment.
Fig. 4 is acceleration and the time-frequency characteristics figure thereof when picking rate is 130km/h in the present invention-embodiment.
embodiment
Below in conjunction with accompanying drawing and embodiment, the invention will be further described, and embodiments of the present invention include but not limited to the following example.
Embodiment 1
As shown in Figure 1, acceleration transducer has been installed on axle box in the present invention, the dynamic vibration information for Dynamic Acquisition locomotive vehicle-mounted under wheel fault, and on backstage, do real-time detection computations.
As shown in Figure 1, performing step of the present invention is as follows:
(1) vibration acceleration sensor is arranged on axle box
Installation of sensors orientation as shown in Figure 2, be selected in ± 50g of the range of vibration transducer wherein, frequency response is not less than 2000Hz, the sample frequency of pick-up transducers is not less than 10kHz, and obtains acceleration signal
Figure DEST_PATH_80897DEST_PATH_IMAGE027
.
Wherein as shown in part on Fig. 3 (a) and Fig. 4, collect one group of axle box Vertical Acceleration signal that contains wheel fault, the speed of collection signal is 150km/h and 130km/h, and corresponding frequecy characteristic and same time peak energy are respectively as shown in part under Fig. 3 (b) and Fig. 4.
(2) carry out temporal frequency signature analysis
Applying frequency section wavelet transformation carries out time-frequency characteristics extraction to axle box vibration acceleration, and the detailed step of feature extraction is as follows:
The first step: the selection of frequency slice function and realization
Frequency coordinate is set , wherein
Figure DEST_PATH_567559DEST_PATH_IMAGE002
the centre frequency that represents conversion,
Figure DEST_PATH_517061DEST_PATH_IMAGE003
represent frequency resolution,
Figure DEST_PATH_809502DEST_PATH_IMAGE004
expression departs from the number of the resolution of centre frequency.
The basis function of frequency choice function is described as:
Then calculate predominant frequency choice function.
(1a) when and time, frequency is selected choice function
Figure DEST_PATH_713556DEST_PATH_IMAGE009
for:
Figure DEST_PATH_108765DEST_PATH_IMAGE028
Wherein
Figure DEST_PATH_126399DEST_PATH_IMAGE011
and
Figure DEST_PATH_214441DEST_PATH_IMAGE012
.
(1b) when
Figure DEST_PATH_74950DEST_PATH_IMAGE008
and
Figure DEST_PATH_273850DEST_PATH_IMAGE013
time, frequency is selected choice function
Figure DEST_PATH_145991DEST_PATH_IMAGE009
for:
Figure DEST_PATH_998409DEST_PATH_IMAGE014
Wherein
Figure DEST_PATH_752739DEST_PATH_IMAGE015
and
Figure DEST_PATH_958592DEST_PATH_IMAGE016
.
(1c) when and
Figure DEST_PATH_505297DEST_PATH_IMAGE018
time, frequency choice function
Figure DEST_PATH_950185DEST_PATH_IMAGE009
for:
Figure DEST_PATH_490888DEST_PATH_IMAGE019
Wherein
Figure DEST_PATH_739553DEST_PATH_IMAGE020
and
Figure DEST_PATH_74719DEST_PATH_IMAGE021
.
In above-mentioned formula,
Figure DEST_PATH_272482DEST_PATH_IMAGE022
for the frequency discrimination factor, centered by coordinate corresponding to frequency, meet
Figure DEST_PATH_442749DEST_PATH_IMAGE029
.
Second step: the fast fourier transform that axle box vibration is accelerated:
Figure DEST_PATH_417659DEST_PATH_IMAGE030
The 3rd step: frequency slice function is selected the frequency of Fourier:
Figure DEST_PATH_633876DEST_PATH_IMAGE031
The 4th step: carry out the time-frequency distributions feature that inverse-Fourier transform obtains vibration acceleration for different center frequency:
(3) judgement wheel track moment loses the identification of contact
The first step: the contour map of drawing the time-frequency distributions of vibration signal;
Second step: the corresponding time period of finding maximum energy peak value;
The 3rd section: at ceiling capacity, in the time period, find successively other time maximum energy peak value, if there is the decay that local energy increases with frequency in section at one time, exist wheel track that the phenomenon departing from occurs.
According to above-described embodiment, just can realize well the present invention.

Claims (2)

1. because wheel-rail impact effect causes wheel track moment, lose a dynamic identifying method for contact, it is characterized in that, realize by the following method:
(1) vibration acceleration sensor is installed on to wheel axle box top, and gathers by the collector on backstage the axle box vibration acceleration information that vibration acceleration sensor measures;
(2) axle box vibration acceleration information is carried out to time-frequency characteristics extraction, obtain time-frequency distributions feature;
(3) draw the contour map of time-frequency distributions feature, judge whether moment loses contact to wheel track;
Wherein, in described step (1), the range of vibration acceleration sensor is ± 50g that frequency response is not less than 2kHz, and the sample frequency of collector is not less than 10kHz;
The specific implementation step of described step (2) is as follows:
(2a) select frequency slice function, wherein frequency coordinate
Figure 2012103996208100001DEST_PATH_IMAGE002
,
Figure 2012103996208100001DEST_PATH_IMAGE004
the centre frequency that represents conversion, represent frequency resolution,
Figure 2012103996208100001DEST_PATH_IMAGE008
expression departs from the number of the resolution of centre frequency, so the basis function of frequency choice function is:
Figure 2012103996208100001DEST_PATH_IMAGE010
According to
Figure 2012103996208100001DEST_PATH_IMAGE012
size, frequency choice function calculates the following three kinds of situations that comprise:
Situation one: when
Figure 2012103996208100001DEST_PATH_IMAGE014
and
Figure 2012103996208100001DEST_PATH_IMAGE016
time, frequency choice function for:
Figure DEST_PATH_IMAGE020
Figure DEST_PATH_IMAGE022
and
Figure DEST_PATH_IMAGE024
;
Situation two: when and
Figure DEST_PATH_IMAGE026
time, frequency choice function
Figure DEST_PATH_IMAGE028
for:
Figure DEST_PATH_IMAGE030
Figure DEST_PATH_IMAGE032
and
Figure DEST_PATH_IMAGE034
;
Situation three: when
Figure DEST_PATH_IMAGE036
and
Figure DEST_PATH_IMAGE038
time, frequency choice function
Figure 730507DEST_PATH_IMAGE028
for:
Figure DEST_PATH_IMAGE040
and
Figure DEST_PATH_IMAGE044
;
Wherein,
Figure DEST_PATH_IMAGE046
for the frequency discrimination factor,
Figure 884145DEST_PATH_IMAGE012
centered by coordinate corresponding to frequency, meet
Figure DEST_PATH_IMAGE048
;
(2b) axle box vibration acceleration is carried out to Fourier transform, transformation for mula is as follows:
Figure DEST_PATH_IMAGE050
(2c) frequency slice function is selected the frequency of Fourier, selects formula as follows:
Figure DEST_PATH_IMAGE052
(2d) for various centre frequencies, carry out the time-frequency distributions feature that inverse-Fourier transform obtains vibration acceleration:
2. according to claim 1ly a kind ofly because wheel-rail impact effect causes wheel track moment, lose the dynamic identifying method of contact, it is characterized in that, the concrete methods of realizing of described step (3) is as follows:
(3a) draw the contour map of the time-frequency distributions feature of vibration acceleration;
(3b) find time period corresponding to this contour map maximum energy peak value;
(3c) at ceiling capacity, in the time period, find successively second largest energy peak, if there is the decay that local energy increases with frequency in section at one time, exist wheel track that the phenomenon departing from occurs.
CN201210399620.8A 2012-10-19 2012-10-19 Dynamic identifying method for instant contact loss of wheel rail under impact effect of wheel rail Active CN102914364B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210399620.8A CN102914364B (en) 2012-10-19 2012-10-19 Dynamic identifying method for instant contact loss of wheel rail under impact effect of wheel rail

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210399620.8A CN102914364B (en) 2012-10-19 2012-10-19 Dynamic identifying method for instant contact loss of wheel rail under impact effect of wheel rail

Publications (2)

Publication Number Publication Date
CN102914364A CN102914364A (en) 2013-02-06
CN102914364B true CN102914364B (en) 2014-02-12

Family

ID=47612835

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210399620.8A Active CN102914364B (en) 2012-10-19 2012-10-19 Dynamic identifying method for instant contact loss of wheel rail under impact effect of wheel rail

Country Status (1)

Country Link
CN (1) CN102914364B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104777054B (en) * 2015-05-18 2017-07-25 浙江工业大学 A kind of parameter identification method of the resonant mode fatigue crack propagation test vibrational system based on soft-measuring technique
DE102015119392A1 (en) * 2015-11-11 2017-05-11 Knorr-Bremse Systeme für Schienenfahrzeuge GmbH Method and device for comparison-controlled derailment detection
JP6884343B2 (en) * 2017-06-02 2021-06-09 曙ブレーキ工業株式会社 Derailment sign detection system, control device, derailment sign detection method, and derailment sign detection program
CN107521511B (en) * 2017-08-18 2019-06-21 唐智科技湖南发展有限公司 A kind of method and system of wheel derailment alarm
CN107521520B (en) * 2017-08-18 2019-04-05 唐智科技湖南发展有限公司 A kind of method and system of wheel derailment early warning
CN110207810B (en) * 2018-02-28 2021-11-30 常州路航轨道交通科技有限公司 Axle box vibration detection system
CN108749849B (en) * 2018-06-05 2020-01-10 唐智科技湖南发展有限公司 Wheel derailment early warning method, device and system

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AT413974B (en) * 2003-05-15 2006-07-15 Siemens Transportation Systems DISCHARGE DETECTION BY FALL SPEED DETERMINATION
JP5191163B2 (en) * 2007-04-27 2013-04-24 株式会社ブリヂストン Tire contact state estimation method and tire contact state estimation device
CN102173299A (en) * 2011-01-21 2011-09-07 合肥安迅铁道应用技术有限公司 Positioning monitoring and managing system for mobile derailer
CN102706568A (en) * 2012-05-22 2012-10-03 西南交通大学 Portable train vibration tester

Also Published As

Publication number Publication date
CN102914364A (en) 2013-02-06

Similar Documents

Publication Publication Date Title
CN102914364B (en) Dynamic identifying method for instant contact loss of wheel rail under impact effect of wheel rail
CN110658005B (en) Method for identifying rail corrugation diseases based on vehicle body acceleration
CN106394606B (en) A kind of railway car wheel loses circle detection method and detection device
CN108845028B (en) Method and device for dynamically detecting high-speed railway rail corrugation
CN108515984B (en) Wheel damage detection method and device
CN106018557B (en) A kind of raceway surface recess detection method based on wavelet analysis
CN105092467B (en) A kind of quick monitoring device of high-speed railway wheel tread flat and method
CN102636577B (en) Method for extracting rim ultrasonic flaw detection signal on basis of wavelet transformation
CN102890143B (en) Rail local defect vehicle-mounting detection method merging with vehicle speed information and front and rear axle box acceleration information
CN105806604A (en) Locomotive vehicle running gear bearing holder fault pre-alarm method
CN103175897B (en) A kind of high-speed switch hurt recognition methods based on vibration signal end-point detection
CN113343928B (en) Method and device for detecting rail wave mill of speed change section high-speed railway and computer equipment
CN107403139A (en) A kind of municipal rail train wheel flat fault detection method
CN108229254B (en) Impact diagnosis method for polygonal out-of-round fault of tread of rail transit wheel
CN103115668A (en) Fault feature extraction method and fault feature recognition method of vibration signal of rolling mill chatter mark
CN109696609A (en) A kind of acquisition methods and its acquisition system of novel bow net current carrying quality
CN111055881A (en) Wheel-rail interface damage evolution monitoring method based on noise signals
CN109059839B (en) Method, device and system for diagnosing out-of-roundness fault of wheel tread
CN107521511B (en) A kind of method and system of wheel derailment alarm
CN103679849B (en) A kind of method that accident doubtful point is judged based on on-board driving recorder
CN104512433A (en) Data collecting and processing method for linear array type infrared axle temperature probes
CN105501250A (en) Fault diagnosis method of scuffed train wheel treads based on vehicle-mounted detection device
CN110514861B (en) High-speed rail running speed estimation method utilizing amplitude spectrum autocorrelation
CN110956164A (en) Train wheel damage identification method based on noise signals
CN106052606A (en) Track surface recess detection method based on scale average wavelet energy spectrum

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