CN110687206B - Ballastless track functional layer defect imaging method - Google Patents

Ballastless track functional layer defect imaging method Download PDF

Info

Publication number
CN110687206B
CN110687206B CN201911074423.7A CN201911074423A CN110687206B CN 110687206 B CN110687206 B CN 110687206B CN 201911074423 A CN201911074423 A CN 201911074423A CN 110687206 B CN110687206 B CN 110687206B
Authority
CN
China
Prior art keywords
defect
ballastless track
echo
functional layer
model
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
CN201911074423.7A
Other languages
Chinese (zh)
Other versions
CN110687206A (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.)
Shijiazhuang Tiedao University
Beijing Shanghai High Speed Railway Co Ltd
Original Assignee
Shijiazhuang Tiedao University
Beijing Shanghai High Speed Railway Co Ltd
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 Shijiazhuang Tiedao University, Beijing Shanghai High Speed Railway Co Ltd filed Critical Shijiazhuang Tiedao University
Priority to CN201911074423.7A priority Critical patent/CN110687206B/en
Publication of CN110687206A publication Critical patent/CN110687206A/en
Application granted granted Critical
Publication of CN110687206B publication Critical patent/CN110687206B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
    • G01N29/0654—Imaging
    • G01N29/069—Defect imaging, localisation and sizing using, e.g. time of flight diffraction [TOFD], synthetic aperture focusing technique [SAFT], Amplituden-Laufzeit-Ortskurven [ALOK] technique
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/4409—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
    • G01N29/4418—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with a model, e.g. best-fit, regression analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/023—Solids
    • G01N2291/0234—Metals, e.g. steel
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/028—Material parameters
    • G01N2291/0289—Internal structure, e.g. defects, grain size, texture

Landscapes

  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Acoustics & Sound (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

The invention relates to the technical field of nondestructive testing of the pouring quality of a ballastless track functional layer of a high-speed railway, in particular to a ballastless track functional layer defect imaging method, which solves the defects of inaccurate geometric dimension estimation and defect imaging of the existing identification method in the prior art, and comprises the following steps: a. acquiring a Burg power spectrum of an echo signal, wherein k is 100; b. calculate [ f2, f1]The ratio eta of the energy spectrum to the total echo energy; c. sorting η in descending order, with the first 10% as the initial value, indicating the possible existence of defect point PiI-1 … q, and initializing i-1; d. searching edge echo frequency distribution mode with Pi as center, if there is echo frequency distribution mode, using the point to edge region as R (P)i) Otherwise, not recording the point; e. i ═ i +1, e.g. i>And q, ending the search, and otherwise, returning to d to continue the search. The method has the characteristics of no damage and high precision, and is suitable for precise imaging of the internal defects of the ballastless track functional layer.

Description

Ballastless track functional layer defect imaging method
Technical Field
The invention relates to the technical field of nondestructive testing of the pouring quality of a functional layer of a ballastless track of a high-speed railway, in particular to a ballastless track functional layer defect imaging method.
Background
The ballastless track of the high-speed railway is used as a high-speed train carrier, the quality is good, whether defects exist in the track or not is judged, and the development degree of the defects is directly related to the operation safety of the high-speed train. In the process of construction and operation of high-speed railways in China, strict regulations are provided for the quality and the state of ballastless tracks, for example, for the pouring quality of self-compacting concrete of CRTS-III slab ballastless tracks, the surface non-area of the self-compacting concrete is definitely required to be more than 50cm2Air bubbles. At present, the concrete pouring quality of the functional layer mainly passes a uncovering test, the method is time-consuming and labor-consuming, and sampling detection is adopted, so that the quality of the functional layer of the whole line is difficult to represent. Therefore, research on nondestructive and rapid functional layer defect detection, size estimation and imaging methods is an urgent problem to be solved for maintaining the health state of the ballastless track of the high-speed railway and ensuring the normal operation of the high-speed railway.
The ballastless track structure of the high-speed railway is a layered concrete structure. The main detection methods for detecting defects inside concrete structures include electromagnetic wave and elastic wave methods. The electromagnetic wave method is limited by the density of the steel bars of the track slab and the narrow space between the track slab surfaces, particularly between the track bearing platforms, and the accurate detection of the defect size is difficult to realize. The concrete parameter detection by using the elastic wave method starts in the 30 th of the 20 th century, and is gradually applied to the detection of the internal defects of the ballastless tracks along with the requirement of the internal defect detection of the high-speed railway in recent years.
The Fourier transform-based method is a common method for detecting and identifying elastic wave defects. Let the echo signal x (N) of length N, whose discrete fourier transform x (m) is expressed as:
Figure BDA0002261965330000021
meanwhile, according to the wave velocity estimated value, the echo frequency range [ f1l, f1h ] of the track slab is calculated](discrete Fourier frequency point correspondence is [ m1, m2 ]]) And obtaining the relation eta of the echo spectrum of the track slab in the total echo energy：
Figure BDA0002261965330000022
And when the eta is larger than a certain threshold value, the position is considered as the cavity disease.
The fourier transform based approach has two drawbacks: the threshold value is high in subjectivity and lacks of defect boundary echo mode analysis; since the defect edge multi-echo frequency f4 is very close to the hole multi-echo frequency f1, it is difficult to satisfactorily distinguish between edge echoes and defect echoes due to insufficient resolution.
The current research mainly analyzes whether a defect exists or not and a defect type identification method, and the defect boundary echo mode research is lacked, so that the boundary identification method is lacked, and the geometric dimension estimation and the defect imaging are inaccurate.
Disclosure of Invention
The invention aims to solve the defects of inaccurate geometric dimension estimation and defect imaging of the existing identification method in the prior art, and provides a ballastless track functional layer defect imaging method.
In order to achieve the purpose, the invention adopts the following technical scheme:
a ballastless track functional layer defect imaging method comprises the following steps:
s1, constructing a flawless ballastless track layered medium model:
reflection coefficient Ri of elastic wave at different interlayer interfaces:
Figure BDA0002261965330000023
multiple reflection frequency fi of fi＝vi/2hi
In the formula, viIndicates the propagation path h of the elastic waveiAverage wave velocity of (1); f. of1Representing the track slab reflection frequency; f. of2Indicating the functional layer reflection frequency; f. of3Representing the support layer reflection frequency;
s2, constructing a ballastless track layered medium model with defects:
a. when the detection unit is far away from the defect position, only three kinds of layered medium interface reflected waves of the ballastless track exist, which are the same as the defect-free layered medium model;
b. when the detection unit is located above the defect and far from the boundary, the echo signal f1 is dominant due to the large difference between the wave impedances Z4 and Z1 of the hole defect;
c. when the detecting unit is located above the defect and is close to the defect boundary, the distance between the defect edge and the sensor is l, and the frequency f4Expressed as: f. of4＝vi/(l+Z1)；
S3, analyzing the mapping relation between the elastic echo and the defect:
a. taking an elastic echo signal x (N) of the ballastless track as an AR model, wherein the length of the echo signal is N, and the echo signal x (N) is formed by linearly combining the first k time signals and N time disturbances a (N):
Figure BDA0002261965330000031
in the formula (I), the compound is shown in the specification,
Figure BDA0002261965330000032
representing the degree of dependence of the k-order model of the signal x (n) on x (n-i), a (n) obeying a normal distribution
Figure BDA0002261965330000033
And independently of x (n-i) (i ═ 1,2, …, k)
The power spectrum of the echo signal is represented as:
Figure BDA0002261965330000034
b. the elastic echo AR model autocorrelation function is represented as:
Figure BDA0002261965330000041
c. the power spectrum of the echo signal is expressed in a matrix form as:
Figure BDA0002261965330000042
d. model parameters
Figure BDA0002261965330000043
The acquisition of (A) is divided into two parts:
Figure BDA0002261965330000044
is obtained and
Figure BDA0002261965330000045
obtaining;
e. the power spectrum of the echo signal defines a k-th order forward prediction error fk, n:
Figure BDA0002261965330000046
f. using X as the value at Xn-kn-k+1，Xn-k+2，…，XnLinear combination representation
Figure BDA0002261965330000047
Figure BDA0002261965330000048
g. Backward prediction error bk, n:
Figure BDA0002261965330000049
h. obtaining a recurrence formula from k-1 order to k order according to the relation between the forward prediction error and the backward prediction error:
Figure BDA00022619653300000410
in the formula (I), the compound is shown in the specification,
Figure BDA00022619653300000411
is the reflection coefficient;
i. extrapolating k-order model parameters;
s4, establishing an elastic echo frequency distribution mode of the detection position-defect position association:
let the echo signal be composed of sinusoidal signals of two frequencies, i.e. the echo signal x (t) is expressed as: x (t) is A1 sin(2πf1t)+A4sin(2πf4t)
In the formula, A1And A4Are respectively f1And f4A frequency signal amplitude;
s5, building ballastless track functional layer defect imaging method based on Burg power spectrum estimation
a. Acquiring a Burg power spectrum of an echo signal, wherein k is 100;
b. calculate [ f ]2,f1]The ratio eta of the energy spectrum to the total echo energy;
c. sorting eta in a descending order, taking the previous 10% as an initial value to indicate that a defect point possibly exists, and initializing i to 1;
d. searching for edge echo frequency distribution mode with Pi as center, and taking the point to edge region as R (P) if there is echo frequency distribution modei) Otherwise, not recording the point;
e. if i is i +1, if i > q, ending the search, otherwise, returning to the step d to continue the search;
f. merging R (P)i) As the final defective area.
Preferably, the ballastless track structure is divided into a track plate, a functional layer and a supporting layer from top to bottom.
Preferably, the excitation source and the sensor together form a detection unit, and the excitation source and the sensor are spaced at a fixed distance.
Preferably, the elastic wave is excited downwards perpendicular to the track slab, reflected waves are generated on the surface of each layer, and continuously reciprocate in the ballastless track to form multiple waves, so that the multiple waves are received by the sensor.
Preferably, in the step of S3,
Figure BDA0002261965330000051
the value meets the requirements of the mathematical expectation of forward error and backward error tending to 0 and the recursion stability
Figure BDA0002261965330000052
According to
Figure BDA0002261965330000053
The calculating method comprises the following steps:
according to
Figure BDA0002261965330000054
It is obtained that,
Figure BDA0002261965330000055
Figure BDA0002261965330000061
according to the initial values f of the forward error, the backward error and the disturbance variance of the echo signals0,n＝b0,n＝x(n)，
Figure BDA0002261965330000062
And sequentially extrapolating k-order model parameters.
Preferably, in step S4, the model order k satisfies the following condition: there are peaks at f1 and f 4; the peak is the largest proportion of the total spectral space at f1 and f 4.
Preferably, the model order k satisfies:
Figure BDA0002261965330000063
in the formula (I), the compound is shown in the specification,
Figure BDA0002261965330000064
representing the amplitude corresponding to the model frequency f of order k.
The invention has the beneficial effects that:
1. a ballastless track layered medium model with defects is built, the mapping relation between elastic echoes and the defects is analyzed, an elastic echo frequency distribution mode of detecting position-defect position association is built, and a ballastless track functional layer defect imaging method based on Burg power spectrum estimation is provided.
2. The method has the characteristics of no damage and high precision, and is suitable for accurate imaging of the internal defects of the ballastless track functional layer.
Drawings
FIG. 1 is a model diagram of a flawless ballastless track layered medium of a ballastless track functional layer defect imaging method provided by the invention;
FIG. 2 is a ballastless track layered medium model with defects of the ballastless track functional layer defect imaging method provided by the invention;
fig. 3 is a simulated echo signal and a Burg power spectrum of the ballastless track functional layer defect imaging method provided by the invention.
In fig. 1 and 2: zi (i ═ 1,2 and 3) represents the wave impedance of each layer of elastic waves, hi represents the thickness from the bottom of each layer to the surface of a track slab in the ballastless track, Reciter represents an excitation source, and Sensor represents a Sensor;
in FIG. 3: graph (a) is sample data containing two sinusoidal signals; graph (b) shows the Burg power spectrum signal when k is 100.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments.
A ballastless track functional layer defect imaging method comprises the following steps:
s1, constructing a flawless ballastless track layered medium model:
the ballastless track structure can be divided into a track slab, a functional layer and a supporting layer from top to bottom, and because the concrete functions and the concrete strengths of all the layers are different, the wave impedance difference affecting the propagation speed of the elastic wave is further caused, so that the elastic wave is reflected on different interfaces, and the ballastless track structure is often expressed as a layered medium model when the defect detection of the ballastless track is researched by using the elastic wave.
The excitation source and the sensor together form a detection unit, the excitation source and the sensor are fixed in distance relative to hiIn other words, the distance is small and is ignored in the calculation process.
And (2) setting an elastic wave vertical to the track slab to be excited downwards, generating reflected waves on the surface of each layer, continuously reciprocating in the ballastless track to form multiple waves, and further receiving the multiple waves by a sensor, wherein the elastic waves have reflection coefficients Ri at interfaces between different layers:
Figure BDA0002261965330000081
multiple reflection frequency fi of fi＝vi/2hi
In the formula, viIndicates the propagation path h of the elastic waveiAverage wave velocity of (1); f. of1Representing the track slab reflection frequency; f. of2Indicating the functional layer reflection frequency; f. of3Indicating the support layer reflection frequency.
S2, constructing a ballastless track layered medium model with defects:
due to the existence of defects of the functional layer, when the relative position of the detection unit and the cavity changes, the propagation characteristic of the elastic wave changes, so that signals received by the sensor have different characteristics.
a. When the detection unit is far away from the defect position, only three kinds of layered medium interface reflected waves of the ballastless track exist, which are the same as the defect-free layered medium model;
b. when the detection unit is located above the defect and far from the boundary, the echo signal f1 is dominant due to the large difference between the wave impedances Z4 and Z1 of the hole defect;
c. when the detection unit is positioned above the defect and is close to the boundary of the defect, besides the frequencies f1, f2 and f3 formed by interface multiple reflection, the defect edge multiple reflection echo also exists, and the frequency is f 4. The edge reflection echo is formed by exciting an elastic wave, and propagating to the edge of the defect according to the Huygens principle to form a new shockA source propagating back and forth between the inside of the track slab and the edge of the defect and the sensor, the distance between the edge of the defect and the sensor being l, the frequency f4Expressed as: f. of4＝vi/(l+Z1)；
Therefore, the echo parameters are changed due to different relative positions of the defect and the detection unit, and the echo frequency distribution mode when the relative positions of the measuring point and the defect are different is listed in table 1.
TABLE 1 Defect echo frequency distribution Pattern
Position of f1 f2 f3 f4
Far from defect Is provided with Is provided with Is provided with Is free of
Over the defect High strength Weak (weak) Weak (weak) Is free of
Close to the boundary Is provided with Is provided with Is provided with Is provided with
S3, analyzing the mapping relation between the elastic echo and the defect:
a. taking an elastic echo signal x (N) of the ballastless track as an AR model, wherein the length of the echo signal is N, and the echo signal x (N) is formed by linearly combining the first k time signals and N time disturbances a (N):
Figure BDA0002261965330000091
in the formula (I), the compound is shown in the specification,
Figure BDA0002261965330000092
representing the degree of dependence of the k-order model of the signal x (n) on x (n-i), a (n) obeying a normal distribution
Figure BDA0002261965330000093
And independently of x (n-i) (i ═ 1,2, …, k)
The power spectrum of the echo signal is represented as:
Figure BDA0002261965330000094
for the AR model, the model parameters, and the model order k are the main technical indicators of the model design.
b. The elastic echo AR model autocorrelation function is represented as:
Figure BDA0002261965330000095
c. the power spectrum of the echo signal is expressed in a matrix form as:
Figure BDA0002261965330000101
d. model parameters
Figure BDA0002261965330000102
The acquisition of (2) is divided into two parts:
Figure BDA0002261965330000103
is obtained and
Figure BDA0002261965330000104
obtaining;
e. the power spectrum of the echo signal defines a k-th order forward prediction error fk, n:
Figure BDA0002261965330000105
f. using X as the value at Xn-kn-k+1，Xn-k+2，…，XnLinear combination representation
Figure BDA0002261965330000106
Figure BDA0002261965330000107
g. Backward prediction error bk, n:
Figure BDA0002261965330000108
h. obtaining a recurrence formula from k-1 order to k order according to the relation between the forward prediction error and the backward prediction error:
Figure BDA0002261965330000109
in the formula (I), the compound is shown in the specification,
Figure BDA00022619653300001010
is the reflection coefficient;
i、
Figure BDA00022619653300001011
the value meets the requirements of the mathematical expectation of forward error and backward error tending to 0 and the recursion stability
Figure BDA00022619653300001012
According to when
Figure BDA00022619653300001013
The calculation method comprises the following steps:
according to
Figure BDA00022619653300001014
It is obtained that,
Figure BDA00022619653300001015
Figure BDA00022619653300001016
according to the initial values f of the forward error, the backward error and the disturbance variance of the echo signals0,n＝b0,n＝x(n)，
Figure BDA0002261965330000111
Sequentially extrapolating k-order model parameters;
s4, establishing a detection position-elastic echo frequency distribution mode associated with the defect position:
determining model order by simulation test method, wherein four echo signals exist near edge defect edge, wherein f1And f2、f3The interval is large, so that the separation is convenient; and a defect echo f1Defect edge echo f4Are relatively close and difficult to distinguish.
Therefore, the echo signal is composed of sinusoidal signals of two frequencies, i.e. the echo signal x (t) is expressed as: x (t) ═ A1sin(2πf1t)+A4sin(2πf4t)
In the formula, A1And A4Are respectively f1And f4A frequency signal amplitude;
the model order k satisfies the following condition: at f1And f4A peak exists; at f1And f4The proportion of the peak value in the whole spectrum space is the largest, namely the model order k satisfies:
Figure BDA0002261965330000112
in the formula (I), the compound is shown in the specification,
Figure BDA0002261965330000113
representing the amplitude corresponding to the model frequency f of order k.
S5, establishing ballastless track functional layer defect imaging method based on Burg power spectrum estimation
a. Acquiring a Burg power spectrum of an echo signal, wherein k is 100;
b. calculate [ f ]2,f1]The ratio eta of the energy spectrum to the total echo energy;
c. sorting eta in a descending order, taking the former 10% as an initial value to represent that a defect point possibly exists, and initializing i to be 1;
d. searching edge echo frequency distribution mode with Pi as center, if there is echo frequency distribution mode, using the point to edge region as R (P)i) Otherwise, not recording the point;
e. if i is i +1, if i > q, ending the search, otherwise, returning to the step d to continue the search;
f. merging R (P)i) As the final defective area.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art should be considered to be within the technical scope of the present invention, and the technical solutions and the inventive concepts thereof according to the present invention should be equivalent or changed within the scope of the present invention.

Claims (5)

1. A ballastless track functional layer defect imaging method is characterized by comprising the following steps:
s1, constructing a flawless ballastless track layered medium model, wherein the ballastless track structure is divided into a track slab, a functional layer and a supporting layer from top to bottom, an elastic wave is positioned at the top of the ballastless track structure and is downwards excited perpendicular to the track slab, a reflected wave is generated on the surface of each layer and continuously reciprocates in the ballastless track to form a multiple wave, and then the multiple wave is received by a sensor:
reflection coefficient Ri of elastic wave at different interlayer interfaces:
Figure FDA0003552029480000011
multiple reflection frequency fi of fi＝vi/2hi
In the formula, viIndicates the propagation path h of the elastic waveiAverage wave velocity of (1); f. of1Representing the track slab reflection frequency; f. of2Indicating the functional layer reflection frequency; f. of3Representing the support layer reflection frequency; zi represents wave impedance of each layer of elastic waves;
s2, constructing a ballastless track layered medium model with defects:
a. when the detection unit is far away from the defect position, the detection unit is the same as a flawless ballastless track layered medium model, and only three layered medium interface reflected waves of the ballastless track exist;
b. when the detection unit is located above the defect and far from the boundary, the echo signal f1 is dominant due to the large difference between the wave impedances Z4 and Z1 of the hole defect;
c. when the detecting unit is located above the defect and is close to the defect boundary, the distance between the defect edge and the sensor is l, and the frequency f4Expressed as: f. of4＝vi/(l+Z1)；
S3, analyzing the mapping relation between the elastic echo and the defect:
a. taking an elastic echo signal x (N) of the ballastless track as an AR model, wherein the length of the echo signal is N, and the echo signal x (N) is formed by linearly combining first k time signals and N time disturbances a (N):
Figure FDA0003552029480000021
in the formula (I), the compound is shown in the specification,
Figure FDA0003552029480000022
representing the degree of dependence of the k-order model of the signal x (n) on x (n-i), a (n) obeying a normal distribution
Figure FDA0003552029480000023
And independently of x (n-i) (i ═ 1,2, …, k);
the power spectrum of the echo signal is represented as:
Figure FDA0003552029480000024
b. the elastic echo AR model autocorrelation function is expressed as:
Figure FDA0003552029480000025
c. the power spectrum of the echo signal is expressed in a matrix form as:
Figure FDA0003552029480000026
d. model parameters
Figure FDA0003552029480000027
The acquisition of (A) is divided into two parts:
Figure FDA0003552029480000028
is obtained and
Figure FDA0003552029480000029
obtaining;
e. the power spectrum of the echo signal defines a k-th order forward prediction error fk, n:
Figure FDA00035520294800000210
f. using X as the value at Xn-kn-k+1，Xn-k+2，…，XnLinear combination representation
Figure FDA00035520294800000211
Figure FDA00035520294800000212
g. Backward prediction error bk, n:
Figure FDA00035520294800000213
h. obtaining a recurrence formula from k-1 order to k order according to the relation between the forward prediction error and the backward prediction error:
Figure FDA0003552029480000031
in the formula (I), the compound is shown in the specification,
Figure FDA0003552029480000032
is the reflection coefficient;
i. extrapolating k-order model parameters;
s4, establishing an elastic echo frequency distribution mode of the detection position-defect position association:
let the echo signal be composed of sinusoidal signals of two frequencies, i.e. the echo signal x (t) is expressed as: x (t) ═ A1sin(2πf1t)+A4sin(2πf4t)
In the formula, A1And A4Are respectively f1And f4A frequency signal amplitude;
s5, establishing ballastless track functional layer defect imaging method based on Burg power spectrum estimation
a. Acquiring a Burg power spectrum of an echo signal, wherein k is 100;
b. calculate [ f ]2,f1]The ratio eta of the energy spectrum to the total echo energy;
c. sorting eta in a descending order, taking the former 10% as an initial value to represent that a defect point possibly exists, and initializing i to be 1;
d. searching edge echo frequency distribution mode with Pi as center, if there is echo frequency distribution mode, using the point to edge region as R (P)i) Otherwise, not recording the point;
e. if i is i +1, if i > q, ending the search, otherwise, returning to the step d to continue the search;
f. merging R (P)i) As the final defective area.
2. The ballastless track functional layer defect imaging method of claim 1, wherein the excitation source and the sensor together form a detection unit, and the distance between the excitation source and the sensor is fixed.
3. The ballastless track functional layer defect imaging method of claim 1, wherein in the step S3,
Figure FDA0003552029480000041
the value meets the requirements of the mathematical expectation of forward error and backward error tending to 0 and the recursion stability
Figure FDA0003552029480000042
According to
Figure FDA0003552029480000043
The calculation method comprises the following steps:
according to
Figure FDA0003552029480000044
So as to obtain the result that,
Figure FDA0003552029480000045
Figure FDA0003552029480000046
according to the initial values f of the forward error, the backward error and the disturbance variance of the echo signals0，n＝b0，n＝x(n)，
Figure FDA0003552029480000047
And sequentially extrapolating k-order model parameters.
4. The ballastless track functional layer defect imaging method of claim 1, wherein in the step S4, the model order k satisfies the following condition: there are peaks at f1 and f 4; the peak is the largest proportion of the total spectral space at f1 and f 4.
5. The ballastless track functional layer defect imaging method of claim 4, wherein model order k satisfies:
Figure FDA0003552029480000048
in the formula (I), the compound is shown in the specification,
Figure FDA0003552029480000049
representing the amplitude corresponding to the model frequency f of order k.
CN201911074423.7A 2019-11-06 2019-11-06 Ballastless track functional layer defect imaging method Active CN110687206B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911074423.7A CN110687206B (en) 2019-11-06 2019-11-06 Ballastless track functional layer defect imaging method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911074423.7A CN110687206B (en) 2019-11-06 2019-11-06 Ballastless track functional layer defect imaging method

Publications (2)

Publication Number Publication Date
CN110687206A CN110687206A (en) 2020-01-14
CN110687206B true CN110687206B (en) 2022-06-03

Family

ID=69115423

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911074423.7A Active CN110687206B (en) 2019-11-06 2019-11-06 Ballastless track functional layer defect imaging method

Country Status (1)

Country Link
CN (1) CN110687206B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116203131B (en) * 2023-04-28 2023-09-15 中国铁建高新装备股份有限公司 Method and device for detecting tunnel void, electronic equipment and storage medium

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09304365A (en) * 1996-05-20 1997-11-28 Tokimec Inc Method and device for determining flaw in subject for examination
JP2005274444A (en) * 2004-03-25 2005-10-06 Toshiba Corp Ultrasonic flaw detection image processor, and processing method therefor
CN104765062A (en) * 2015-04-13 2015-07-08 四川升拓检测技术有限责任公司 Ballastless track board disengaging nondestructive detection method based on elastic waves
CN105783799A (en) * 2016-03-03 2016-07-20 四川升拓检测技术股份有限公司 Ballastless track plate seam depth non-destructive detection method and equipment based on vibration
CN106546982A (en) * 2016-11-10 2017-03-29 桂林电子科技大学 A kind of waveform inversion method of the layered medium based on simulated annealing
CN109033618A (en) * 2018-07-24 2018-12-18 中南大学 The appraisal procedure that non-fragment orbit typical case hurt influences bullet train safety in operation
CN109164173A (en) * 2018-10-08 2019-01-08 上海工程技术大学 A kind of method and device of multichannel Dynamic Non-Destruction Measurement non-fragment orbit defect
CN109239191A (en) * 2018-09-29 2019-01-18 中国特种设备检测研究院 A kind of supersonic guide-wave defect location imaging method and system
CN109738525A (en) * 2019-03-12 2019-05-10 石家庄铁道大学 The high-speed railway bearing layer concrete velocity of wave detection device and estimation method of gridding
CN109813809A (en) * 2019-04-01 2019-05-28 石家庄铁道大学 Non-fragment orbit defect non-contact non-destructive testing method, terminal device and system

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09304365A (en) * 1996-05-20 1997-11-28 Tokimec Inc Method and device for determining flaw in subject for examination
JP2005274444A (en) * 2004-03-25 2005-10-06 Toshiba Corp Ultrasonic flaw detection image processor, and processing method therefor
CN104765062A (en) * 2015-04-13 2015-07-08 四川升拓检测技术有限责任公司 Ballastless track board disengaging nondestructive detection method based on elastic waves
CN105783799A (en) * 2016-03-03 2016-07-20 四川升拓检测技术股份有限公司 Ballastless track plate seam depth non-destructive detection method and equipment based on vibration
CN106546982A (en) * 2016-11-10 2017-03-29 桂林电子科技大学 A kind of waveform inversion method of the layered medium based on simulated annealing
CN109033618A (en) * 2018-07-24 2018-12-18 中南大学 The appraisal procedure that non-fragment orbit typical case hurt influences bullet train safety in operation
CN109239191A (en) * 2018-09-29 2019-01-18 中国特种设备检测研究院 A kind of supersonic guide-wave defect location imaging method and system
CN109164173A (en) * 2018-10-08 2019-01-08 上海工程技术大学 A kind of method and device of multichannel Dynamic Non-Destruction Measurement non-fragment orbit defect
CN109738525A (en) * 2019-03-12 2019-05-10 石家庄铁道大学 The high-speed railway bearing layer concrete velocity of wave detection device and estimation method of gridding
CN109813809A (en) * 2019-04-01 2019-05-28 石家庄铁道大学 Non-fragment orbit defect non-contact non-destructive testing method, terminal device and system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"AR模型功率谱估计的Burg算法研究及验证";李长远 等;《铁道通信信号》;20180103;全文 *
"一种改进型AR-EWT算法及其在振动模态分解中的应用研究";罗治军 等;《机械设计》;20190731;全文 *
"基于AR模型的Yule-Walker法和Burg法功率谱估计性能分析";*** 等;《计算机与数字工程》;20160531;第44卷(第5期);全文 *

Also Published As

Publication number Publication date
CN110687206A (en) 2020-01-14

Similar Documents

Publication Publication Date Title
Hou et al. Automatic multi-mode Lamb wave arrival time extraction for improved tomographic reconstruction
CN102853791B (en) Method for scanning ultrasonic microscope and measuring thickness, sound velocity, density and attenuation of thin material simultaneously
CN103245311B (en) With the thickness measuring method of Ultrasonic Detection Multilayer Absorbing Material Coating measuring thickness device
CN110243320B (en) Tunnel lining crack depth non-contact measurement method and device
CN106198739A (en) A kind of TOFD near surface blind region defect location detection method based on shape transformation
CN112098526A (en) Near-surface defect feature extraction method for additive product based on laser ultrasonic technology
CN101109732B (en) Ultrasound nondestructive detecting echo signal classificating method based on vague plane characteristic
CN103616439A (en) Method for simultaneously measuring multiple parameters of linear visco-elastic thin layer material by employing ultrasonic flat probe
CN104698089A (en) Ultrasonic relative time propagation technology suitable for inclined crack quantifying and imaging
CN110687206B (en) Ballastless track functional layer defect imaging method
CN104897777A (en) Method for improving longitudinal resolution of TOFD (time of flight diffraction) detection with Burg algorithm based autoregressive spectrum extrapolation technology
CN101571519A (en) Ultrasonic guided wave detection technology for quantifying defects of composite laminated plate
CN113533510B (en) Rail fatigue micro-crack identification method and device
CN106770668B (en) Method for detecting quality of single-hole foundation pile by acoustic transmission method
JP5562118B2 (en) Ultrasonic nondestructive measuring method, ultrasonic nondestructive measuring device, and program
CN113533504B (en) Subsurface crack quantitative measurement method based on laser ultrasonic surface wave frequency domain parameters
CN106680821A (en) Ultrasonic damage-free method for detecting thickness of NiCoCrAlYTa hexabasic coating plasma spraying coating
CN109541689B (en) Method for evaluating compactness of medium based on reflected wave energy characteristics
CN115639157A (en) Surface wave-based surface crack position, length and angle measurement method
CN210180982U (en) Non-contact nondestructive testing system for ballastless track defects
JP2001343365A (en) Thickness resonance spectrum measuring method for metal sheet and electromagnetic ultrasonic measuring method for metal sheet
CN114755300A (en) Defect positioning quantitative detection method based on ultrasonic nondestructive detection
RU2714868C1 (en) Method of detecting pitting corrosion
JP2003149214A (en) Nondestructive inspecting method and its apparatus using ultrasonic sensor
Marin-Cortes Ultrasonic coda wave comparison for quality control of manufactured parts: Proof of feasibility

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant