[go: up one dir, main page]

CN111452830A - An imaging method and device for realizing automatic detection of track slab cracks - Google Patents

An imaging method and device for realizing automatic detection of track slab cracks Download PDF

Info

Publication number
CN111452830A
CN111452830A CN202010274433.1A CN202010274433A CN111452830A CN 111452830 A CN111452830 A CN 111452830A CN 202010274433 A CN202010274433 A CN 202010274433A CN 111452830 A CN111452830 A CN 111452830A
Authority
CN
China
Prior art keywords
receiving
ultrasonic
function
green
signal
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.)
Granted
Application number
CN202010274433.1A
Other languages
Chinese (zh)
Other versions
CN111452830B (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.)
Shanghai University of Engineering Science
Original Assignee
Shanghai University of Engineering Science
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 Shanghai University of Engineering Science filed Critical Shanghai University of Engineering Science
Priority to CN202010274433.1A priority Critical patent/CN111452830B/en
Publication of CN111452830A publication Critical patent/CN111452830A/en
Application granted granted Critical
Publication of CN111452830B publication Critical patent/CN111452830B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B61RAILWAYS
    • B61KAUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
    • B61K9/00Railway vehicle profile gauges; Detecting or indicating overheating of components; Apparatus on locomotives or cars to indicate bad track sections; General design of track recording vehicles
    • B61K9/08Measuring installations for surveying permanent way
    • B61K9/10Measuring installations for surveying permanent way for detecting cracks in rails or welds thereof
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种实现轨道板裂缝自动检测的成像方法及装置,所述成像方法是指,先利用摄像仪检测轨道板裂缝表面缺陷,再利用超声波发射接收器激发并经由超声线性阵列探头发射频率在1~5MHz范围内的超声波信号至轨道板内部,由超声线性阵列探头采用自发自收的模式获取轨道板中的扩散场信号并将扩散场信号传输给计算机,最后由计算机采用MATLAB软件对所接收到的信号被动提取格林函数,对扩散场信号进行互相关,重建阵元之间的格林函数,获取阵元之间未延时的响应,恢复被噪声淹没的早期缺陷信息,然后依据频域合成孔径聚焦成像算法进行轨道板裂缝成像。本发明可以高效、无损、实时的检测无砟轨道板裂缝表面及裂缝内部缺陷。

Figure 202010274433

The invention discloses an imaging method and device for realizing automatic detection of track slab cracks. The imaging method refers to firstly using a camera to detect surface defects of track slab cracks, and then using an ultrasonic transmitter and receiver to excite and transmit through an ultrasonic linear array probe. The ultrasonic signal with a frequency in the range of 1-5MHz is sent to the inside of the track plate. The ultrasonic linear array probe adopts the mode of self-sending and self-receiving to obtain the diffuse field signal in the track plate and transmits the diffuse field signal to the computer. Finally, the computer uses MATLAB software to analyze the signal. The received signal passively extracts the Green's function, performs cross-correlation on the diffuse field signal, reconstructs the Green's function between the array elements, obtains the undelayed response between the array elements, recovers the early defect information submerged by noise, and then according to the frequency Orbital plate crack imaging using domain synthetic aperture focusing imaging algorithm. The invention can efficiently, non-destructively and real-time detect the crack surface of the ballastless track slab and the defects inside the crack.

Figure 202010274433

Description

一种实现轨道板裂缝自动检测的成像方法及装置An imaging method and device for realizing automatic detection of track slab cracks

技术领域technical field

本发明涉及一种实现轨道板裂缝自动检测的成像方法及装置,属于轨道缺陷检测技术领域。The invention relates to an imaging method and device for realizing automatic detection of track slab cracks, belonging to the technical field of track defect detection.

背景技术Background technique

高速铁路对于缩短城市距离、推动区域协调发展有着重要的意义,是现代交通运输方式中不可或缺的一部分。随着高速铁路的迅猛发展,无砟轨道作为一种主要的轨道结构,是由无砟无砟轨道板、CA砂浆层、支撑层和基床构成,其运用范围越来越广。然而近年来,高铁线下结构出现了越来越多的病害,包括线下结构离隙、贯穿裂缝以及CA砂浆层脱空等病害。主要原因一方面是列车高速重载运行过程中会对无砟轨道产生挤压、冲击等作用,导致其内部可能会出现不密实、裂缝或空洞,外部形成损伤层或蜂窝麻层等各种各样的缺陷;另一方面是由于无砟轨道在前期制作中因施工工艺、施工经验可能存在问题,导致本身就存在缺陷;另外,雨雪的侵蚀、环境温度变化等自然灾害下也将导致产生缺陷。因缺陷的存在将严重影响无砟轨道的承载力和耐久性,将会致使无砟轨道结构失效,无法保证高速铁路无砟轨道及线下结构的稳定性和平顺性,而稳定性和平顺性恰恰是保证高铁快速和安全运营的重要前提条件,将直接关系到列车的正常运营和乘客的人身安全。High-speed railways are of great significance for shortening urban distances and promoting coordinated regional development, and are an indispensable part of modern transportation methods. With the rapid development of high-speed railway, ballastless track, as a main track structure, is composed of ballastless track slab, CA mortar layer, support layer and foundation bed, and its application range is getting wider and wider. However, in recent years, more and more diseases have appeared in the structure of high-speed rail lines, including the under-line structure relief, through cracks, and voids in the CA mortar layer. The main reason is that, on the one hand, during the high-speed and heavy-load operation of the train, the ballastless track will be squeezed, impacted, etc., which may cause incompactness, cracks or voids inside, and the formation of damage layers or honeycomb hemp layers on the outside. On the other hand, the ballastless track may have problems in the pre-production due to construction technology and construction experience, resulting in its own defects; in addition, natural disasters such as rain and snow erosion and changes in ambient temperature will also lead to defect. The existence of defects will seriously affect the bearing capacity and durability of the ballastless track, which will lead to the failure of the ballastless track structure. It is precisely the important prerequisite for ensuring the rapid and safe operation of high-speed rail, which will be directly related to the normal operation of the train and the personal safety of passengers.

但目前,我国实现无砟轨道缺陷检测主要依靠人工静态检测技术,由于轨道交通用于可供线路检修维护的有效天窗时间仅为2-3小时,且高速铁路的线程又很长,若采用现有的检测手段不仅耗费大量人力物力,而且效率十分低下,检测维护成本很高,以致不能满足轨道安全预警需求。However, at present, the detection of ballastless track defects in my country mainly relies on manual static detection technology. Since the effective skylight time for rail transit to be used for line maintenance and maintenance is only 2-3 hours, and the thread of high-speed railway is very long, if the existing Some detection methods not only consume a lot of manpower and material resources, but also have very low efficiency and high detection and maintenance costs, so that they cannot meet the needs of rail safety early warning.

发明内容SUMMARY OF THE INVENTION

针对现有技术存在的上述问题,本发明的目的是提供一种实现轨道板裂缝自动检测的成像方法及装置,以实现高效、无损、实时检测出无砟轨道板裂缝缺陷,为高铁的安全运营提供及时预警和有力保障。In view of the above-mentioned problems existing in the prior art, the purpose of the present invention is to provide an imaging method and device for realizing automatic detection of track slab cracks, so as to realize efficient, non-destructive and real-time detection of ballastless track slab crack defects, and for the safe operation of high-speed railways Provide timely warning and strong guarantee.

为实现上述目的,本发明采用如下技术方案:To achieve the above object, the present invention adopts the following technical solutions:

一种实现轨道板裂缝自动检测的成像方法,包括如下步骤:An imaging method for realizing automatic detection of track slab cracks, comprising the following steps:

a)利用摄像仪检测轨道板裂缝表面缺陷,具体是指:使用摄像仪对轨道板表面进行拍摄,并将拍摄的图片传递给计算机,计算机根据预存的正常轨道板图片与之比对,以判断轨道板表面是否存在裂缝;a) Using a camera to detect surface defects of track slab cracks, specifically: using a camera to shoot the surface of the track slab, and transfer the captured image to the computer, and the computer compares it with the pre-stored normal track slab image to judge Whether there are cracks on the surface of the track plate;

b)利用超声波发射接收器激发并经由超声线性阵列探头发射频率在1~5MHz范围内的超声波信号至轨道板内部,由超声线性阵列探头采用自发自收的模式获取轨道板中的扩散场信号并将扩散场信号传输给计算机,具体是指:当计算机判断结果显示存在裂缝时,则调节超声线性阵列探头的位置使其下沉压置于裂缝处,然后启动超声波发射接收器,使超声波发射接收器的发射端发射频率在1~5MHz范围内的超声波信号,并由超声线性阵列探头的发射端发射至轨道板内部,然后由超声线性阵列探头的接收端采集扩散场信号并将其传输给超声波发射接收器的接收端,然后由超声波发射接收器的接收端将接收的扩散场信号传输给计算机;b) Use the ultrasonic transmitter and receiver to excite and transmit ultrasonic signals with a frequency in the range of 1 to 5MHz to the inside of the track plate through the ultrasonic linear array probe. Transmit the diffused field signal to the computer, specifically: when the computer determines that there is a crack, adjust the position of the ultrasonic linear array probe to make it sink and place it at the crack, and then start the ultrasonic transmitter and receiver to transmit and receive the ultrasonic wave. The transmitting end of the ultrasonic sensor transmits ultrasonic signals with a frequency in the range of 1 to 5 MHz, and is transmitted to the inside of the track plate by the transmitting end of the ultrasonic linear array probe, and then the diffused field signal is collected by the receiving end of the ultrasonic linear array probe and transmitted to the ultrasonic wave. The receiving end of the transmitter and receiver, and then the received end of the ultrasonic transmitter receiver transmits the received diffuse field signal to the computer;

c)由计算机采用MATLAB软件对所接收到的信号被动提取格林函数,对扩散场信号进行互相关,重建阵元之间的格林函数,获取阵元之间未延时的响应,恢复被噪声淹没的早期缺陷信息,然后依据频域合成孔径聚焦成像算法进行轨道板裂缝成像。c) The computer uses MATLAB software to passively extract the Green's function from the received signal, cross-correlate the diffuse field signal, reconstruct the Green's function between the array elements, obtain the undelayed response between the array elements, and recover from being submerged by noise. The early defect information is obtained, and then the track plate crack imaging is carried out according to the frequency domain synthetic aperture focusing imaging algorithm.

一种实施方案,所述步骤c)具体包括如下操作:In one embodiment, the step c) specifically comprises the following operations:

1)数据处理:1) Data processing:

首先,由于超声波信号在无砟轨道板内部传播的过程中,超声波信号和缺陷目标之间发生相互作用,经过散射和多次的反射后会形成一个近似均匀的声场,该声场即为扩散场,而如步骤b)所示,通过超声线性阵列探头即可获取无砟轨道中的扩散场信号,设定超声线性阵列探头中任意的两个接收阵元分别为γ1和γ2,两个接收阵元γ1和γ2处于封闭的曲面空间内,则接收阵元γ1和γ2处声场互谱的所有声源积分等于γ1和γ2之间频域因果格林函数与非因果格林函数之差:First of all, due to the interaction between the ultrasonic signal and the defect target during the propagation of the ultrasonic signal inside the ballastless track slab, an approximately uniform sound field will be formed after scattering and multiple reflections, which is the diffuse field. As shown in step b), the diffused field signal in the ballastless track can be obtained through the ultrasonic linear array probe, and any two receiving array elements in the ultrasonic linear array probe are set to be γ 1 and γ 2 respectively, and the two receiving elements are If the array elements γ 1 and γ 2 are in a closed surface space, then the integral of all sound sources of the sound field cross-spectrum at the receiving array elements γ 1 and γ 2 is equal to the frequency domain causal Green function and non-causal Green function between γ 1 and γ 2 Difference:

G(r1,r2ω)-G*(r1,r2,ω)=-2jω∫G(r1,r,ω)G*(r2,r,ω)dV (1):G(r 1 , r 2 ω)-G * (r 1 ,r 2 ,ω)=-2jω∫G(r 1 ,r,ω)G * (r 2 ,r,ω)dV(1):

公式(1)中,等式的左侧是两个接收阵元γ1和γ2之间的频域格林函数的因果和非因果函数之差,等式右侧表示接收阵元γ1和γ2处声场互相关的所有声源积分,具体的,G(γ12,ω)代表频域中γ1和γ2之间的因果格林函数,上角标*代表共轭复数,G*12,ω)代表频域中γ1和γ2之间的非因果格林函数,j代表虚数单位,j2=-1,ω代表发射信号的角频率,γ是处于扩散场积分密闭空间V中的任意一点位置,代表任意的缺陷目标,可视为噪声源,G(γ1,γ,ω)代表γ和γ1之间的频域格林函数传播公式,G(γ2,γ,ω)代表γ和γ2之间的频域格林函数传播公式,dV代表密闭空间V的微分;In formula (1), the left side of the equation is the difference between the causal and non-causal functions of the frequency domain Green function between the two receiving elements γ 1 and γ 2 , and the right side of the equation represents the receiving elements γ 1 and γ Integral of all sound sources in the cross-correlation of sound fields at 2 places, specifically, G(γ 12 ,ω) represents the causal Green function between γ 1 and γ 2 in the frequency domain, the superscript * represents the conjugate complex number, G *12 ,ω) represents the non-causal Green's function between γ 1 and γ 2 in the frequency domain, j represents the imaginary unit, j 2 =-1, ω represents the angular frequency of the transmitted signal, γ is in the diffusion The position of any point in the closed space V of field integration represents any defect target and can be regarded as a noise source, G(γ 1 ,γ,ω) represents the frequency domain Green’s function propagation formula between γ and γ 1 , G(γ 2 , γ, ω) represents the frequency domain Green's function propagation formula between γ and γ 2 , dV represents the differential of the closed space V;

然后,假设扩散场噪声源强度是空间均匀分布且互不相关,则功率谱密度|q(ω)|2同时与噪声源的位置无关:Then, the power spectral density |q(ω)| 2 is also independent of the position of the noise source, assuming that the intensity of the diffuse-field noise source is spatially uniformly distributed and uncorrelated with each other:

(q(r1,ω)q*(r2,ω)>=δ(r1-r2)|q(ω)|2 (2);(q(r 1 , ω)q * (r 2 , ω)>=δ(r 1 -r 2 )|q(ω)| 2 (2);

公式(2)代表频域γ1和γ2处信号的互功率谱的期望值等于γ1和γ2间噪声源的功率谱密度,具体的,公式(2)中,<q(γ1,ω)q*2,ω)>代表频域γ1和γ2处信号的互功率谱的期望值,q(γ1,ω)代表γ1处的场强,q(γ2,ω)代表γ2处的场强,<>代表统计平均的期望值,上角标*代表共轭复数,δ(γ1-γ2)代表γ1和γ2之间信号的互功率谱与噪声源的位置无关;频域中两个接收阵元γ1和γ2位置处声场之间的互相关为:Equation (2) represents the expected value of the cross-power spectrum of the signals at γ 1 and γ 2 in the frequency domain equal to the power spectral density of the noise source between γ 1 and γ 2. Specifically, in formula (2), <q(γ 1 , ω )q *2 , ω)> represents the expected value of the cross-power spectrum of the signals at γ 1 and γ 2 in the frequency domain, q(γ 1 , ω) represents the field strength at γ 1 , q(γ 2 , ω) represents the The field strength at γ 2 , <> represents the expected value of the statistical average, the superscript * represents the conjugate complex number, δ(γ1-γ2) represents the cross-power spectrum of the signal between γ 1 and γ 2 is independent of the position of the noise source; The cross-correlation between the sound fields at the positions of the two receiving array elements γ 1 and γ 2 in the frequency domain is:

<p(r1,ω)p*(r2,ω)>=|q(ω)|2∫G(r1,r,ω)G*(r2,r,ω)dV (3);<p(r 1 , ω)p * (r 2 , ω)>=|q(ω)| 2 ∫G(r 1 , r, ω)G * (r 2 , r, ω)dV (3);

公式(3)中,p(γ1,ω)代表接收阵元γ1处的声场,p(γ2,ω)代表接收阵元γ2处的声场,<p(γ1,ω)p*2,ω)>代表两个接收阵元γ1和γ2位置处声场之间的互相关函数,上角标*代表共轭复数;In formula (3), p(γ 1 , ω) represents the sound field at the receiving element γ 1 , p(γ 2 , ω) represents the sound field at the receiving element γ 2 , <p(γ 1 , ω)p *2 , ω)> represents the cross-correlation function between the sound fields at the positions of the two receiving array elements γ 1 and γ 2 , and the superscript * represents a complex conjugate number;

然后,由公式(1)和(3)得到:Then, from formulas (1) and (3), we get:

(G(r1,r2,ω)-G*(r1,r2,ω))|q(ω)|2=-2jω<p(r1,ω)p*(r2,ω)> (4);(G(r 1 , r 2 , ω)-G * (r 1 , r 2 , ω))|q(ω)| 2 =-2jω<p(r 1 , ω)p * (r 2 , ω) >(4);

公式(4)中,等式左侧为格林函数G(γ1,γ2,ω)以及它对应的时间反转,即对应频域中的共轭操作,把这两项乘以随机噪声的功率谱密度后,其结果等于等式右侧扩散场中两个的接收阵元γ1和γ2互相关结果,上角标*代表共轭复数;In formula (4), the left side of the equation is the Green's function G(γ 1 , γ 2 , ω) and its corresponding time reversal, that is, the conjugate operation in the corresponding frequency domain. Multiply these two terms by the random noise. After the power spectral density, the result is equal to the cross-correlation result of the two receiving array elements γ 1 and γ 2 in the diffuse field on the right side of the equation, and the superscript * represents the conjugate complex number;

然后,求出公式(4)所对应的时域表达式,频域2jω对应时域2d/dt,根据卷积定理可知,频域中的乘积对应时域中的卷积,从而得到:Then, the time domain expression corresponding to formula (4) is obtained. The frequency domain 2jω corresponds to the time domain 2d/dt. According to the convolution theorem, the product in the frequency domain corresponds to the convolution in the time domain, thus obtaining:

Figure BDA0002444267950000031
Figure BDA0002444267950000031

公式(5)中,G(γ1,γ2,t)代表时域中γ1和γ2之间的格林函数,G(γ1,γ2,-t)代表G(γ1,γ2,t)的时间反转,即对应频域中的共轭操作,*代表卷积运算,

Figure BDA0002444267950000041
代表互相关运算,Cq(t)表示扩散场中的噪声q(t)的自相关结果,d/dt代表对t求导,p(γ1,t)代表时域中接收阵元γ1处的声场,p(γ2,t)代表时域中接收阵元γ2处的声场;In formula (5), G(γ 1 , γ 2 , t) represents the Green's function between γ 1 and γ 2 in the time domain, and G(γ 1 , γ 2 , -t) represents G(γ 1 , γ 2 , t) time reversal, which corresponds to the conjugate operation in the frequency domain, * represents the convolution operation,
Figure BDA0002444267950000041
represents the cross-correlation operation, C q (t) represents the autocorrelation result of the noise q(t) in the diffuse field, d/dt represents the derivation of t, p(γ 1 ,t) represents the receiving array element γ 1 in the time domain The sound field at , p(γ 2 ,t) represents the sound field at the receiving array element γ 2 in the time domain;

公式(4)和(5)的结果表明扩散场中接收阵元声场p(γ1,t)和p(γ2,t)互相关并求导,其结果等于这两接收阵元之间的格林函数响应,根据声波的互易性原理,在时间轴上,这两接收阵元之间的格林函数响应具有对称性;The results of formulas (4) and (5) show that the sound fields p(γ 1 ,t) and p(γ 2 ,t) of the receiving array elements in the diffuse field are correlated and derived, and the result is equal to the difference between the two receiving array elements. Green's function response, according to the principle of reciprocity of sound waves, on the time axis, the Green's function response between the two receiving array elements has symmetry;

2)频域合成孔径聚焦(SAFT)成像:2) Frequency Domain Synthetic Aperture Focusing (SAFT) imaging:

首先,由于超声波在无限的均匀介质传播满足Helmholtz波动方程,因此,假设在二维空间x-z坐标轴下,超声波发射阵元的坐标为(u,0),接收阵元的坐标为(v,0),缺陷目标距离发射阵元和接收阵元的距离分别为ρin和ρout,任意一对发射接收阵元对应的时域接收信号为e(v,u,t),那么对应的发射接收阵元所对应的频率响应E(v,u,ω)为:First of all, since the propagation of ultrasonic waves in an infinite homogeneous medium satisfies the Helmholtz wave equation, it is assumed that the coordinates of the ultrasonic transmitting array element are (u, 0) and the coordinates of the receiving array element are (v, 0 under the xz coordinate axis of the two-dimensional space). ), the distances of the defect target from the transmitting array element and the receiving array element are ρ in and ρ out respectively, and the time domain receiving signal corresponding to any pair of transmitting and receiving array elements is e(v, u, t), then the corresponding transmitting and receiving The frequency response E(v, u, ω) corresponding to the array element is:

E(v,u,ω)=P(ω)∫∫f(x,z)g(ρtn,ω)g(ρout,ω)dxdz (6):E(v, u, ω)=P(ω)∫∫f(x, z)g(ρ tn , ω)g(ρ out , ω)dxdz (6):

公式(6)中,

Figure BDA0002444267950000042
P(ω)是发射信号的频谱,ω为发射信号的角频率;f(x,z)代表缺陷目标的点扩散函数,g(ρin,ω)代表发射阵元到缺陷目标的空间格林函数频域响应,g(ρout,ω)代表接收阵元到缺陷目标的空间格林函数频域响应;In formula (6),
Figure BDA0002444267950000042
P(ω) is the spectrum of the transmitted signal, ω is the angular frequency of the transmitted signal; f(x, z) represents the point spread function of the defective target, and g(ρ in , ω) represents the spatial Green's function from the transmitting array element to the defective target Frequency domain response, g(ρ out , ω) represents the spatial Green's function frequency domain response from the receiving array element to the defect target;

此外,格林函数为:Also, the Green's function is:

Figure BDA0002444267950000043
Figure BDA0002444267950000043

公式(7)中,g(x,z,ω)代表二维空间x-z坐标轴下任意点的格林函数,波数k=ω/c,c为声波在介质中的纵波速度,空间变量x和z与波数kx和kz相对应;In formula (7), g(x, z, ω) represents the Green’s function of any point under the xz coordinate axis of the two-dimensional space, the wave number k=ω/c, c is the longitudinal wave velocity of the sound wave in the medium, and the spatial variables x and z Corresponds to the wave numbers k x and k z ;

使用变量x和z代替掉变量ρin和ρout,即g(x,z,ω)代替g(ρin、ρout,ω),将公式(7)带入公式(6),重新获得接收响应E(v,u,ω),然后用波数ku和kv分别表示发射阵元(u,0)和接收阵元(v,0)的空间信息,为了实现将积分转换二维傅里叶变换,对接收响应E(v,u,ω)中变量u和v进行傅里叶变换得到:Use the variables x and z to replace the variables ρ in and ρ out , that is, g(x, z, ω) instead of g(ρ in , ρ out , ω), bring the formula (7) into the formula (6), and regain the reception Response E(v, u, ω), and then use the wave numbers ku and k v to represent the spatial information of the transmitting array element (u, 0) and the receiving array element (v, 0), respectively. In order to realize the transformation of the integral into a two-dimensional Fourier Leaf transform, the Fourier transform of the variables u and v in the received response E(v, u, ω) is obtained:

Figure BDA0002444267950000051
Figure BDA0002444267950000051

公式(8)中,E(ω,ku,Kv)代表频域响应E(v,u,ω)的空间傅里叶变换,F代表点扩散函数f(x,z)的空间傅里叶变换;In formula (8), E(ω, ku , K v ) represents the spatial Fourier transform of the frequency domain response E(v, u , ω), and F represents the spatial Fourier transform of the point spread function f(x, z) leaf transform;

相应的,点扩散函数f(x,z)的空间傅里叶变换是F(kx,kz),为了求得图像域,反变换模型为:Correspondingly, the spatial Fourier transform of the point spread function f(x, z) is F(k x , k z ). In order to obtain the image domain, the inverse transform model is:

Figure BDA0002444267950000052
Figure BDA0002444267950000052

公式(9)中,

Figure BDA0002444267950000053
代表缺陷目标等间隔采样的数据域模型,S-1{}代表Stolt映射,完成了从数据域到图像域的转换,p*(ω)代表发射信号的共轭频谱;In formula (9),
Figure BDA0002444267950000053
represents the data domain model of the defect target sampled at equal intervals, S -1 {} represents the Stolt map, which completes the conversion from the data domain to the image domain, p * (ω) represents the conjugate spectrum of the transmitted signal;

对公式(9)进行非线性坐标变化,将波数ku,kv和k映射到和图像域缺陷目标相关联的波数kx和kz,则:The nonlinear coordinate change of formula (9) is performed to map the wavenumbers ku, kv and k to the wavenumbers kx and kz associated with the defect target in the image domain, then:

kx=ku+kv (10):k x = k u + k v (10):

Figure BDA0002444267950000061
Figure BDA0002444267950000061

然后通过二维傅里叶反变换求得缺陷目标分布的图像Then the image of the defect target distribution is obtained by the two-dimensional inverse Fourier transform

Figure BDA0002444267950000062
Figure BDA0002444267950000062

公式(12)中,

Figure BDA0002444267950000063
代表缺陷目标分布的图像域模型。In formula (12),
Figure BDA0002444267950000063
An image domain model representing the distribution of defect targets.

一种实现轨道板裂缝自动检测的成像装置,包括轨检小车及设置在轨检小车上的摄像仪和计算机,所述摄像仪与计算机信号连接,还包括超声波发射接收器和超声线性阵列探头,所述超声线性阵列探头的发射端与超声波发射接收器的发射接口信号连接,所述超声线性阵列探头的接收端与超声波发射接收器的接收接口信号连接,所述超声波发射接收器的接收接口与计算机信号连接,所述超声线性阵列探头与一自适应伸缩机构相连接,所述自适应伸缩机构与一双轴位移调节机构相连接,且所述双轴位移调节机构固定连接在轨检小车的前端。An imaging device for realizing automatic detection of track slab cracks includes a track inspection trolley, a camera and a computer arranged on the track inspection trolley, the camera is connected with a computer signal, and also includes an ultrasonic transmitter receiver and an ultrasonic linear array probe, The transmitting end of the ultrasonic linear array probe is signally connected to the transmitting interface of the ultrasonic transmitting receiver, the receiving end of the ultrasonic linear array probe is signally connecting to the receiving interface of the ultrasonic transmitting and receiving device, and the receiving interface of the ultrasonic transmitting and receiver is connected to the receiving interface of the ultrasonic transmitting receiver. computer signal connection, the ultrasonic linear array probe is connected with an adaptive telescopic mechanism, the adaptive telescopic mechanism is connected with a biaxial displacement adjustment mechanism, and the biaxial displacement adjustment mechanism is fixedly connected to the front end of the rail inspection trolley .

一种实施方案,所述自适应伸缩机构包括垂向设置的、互相平行的安装板A和安装板B,所述安装板A和安装板B之间连接有可水平伸缩的伸缩支架,所述超声线性阵列探头与安装板A相连接,所述双轴位移调节机构与安装板B相连接。In one embodiment, the self-adaptive telescopic mechanism includes a mounting plate A and a mounting plate B that are arranged vertically and parallel to each other, and a horizontally telescopic telescopic bracket is connected between the mounting plate A and the mounting plate B. The ultrasonic linear array probe is connected with the mounting plate A, and the biaxial displacement adjusting mechanism is connected with the mounting plate B.

一种优选方案,所述安装板A远离安装板B的一侧设有探头连接件,所述超声线性阵列探头固定安装在探头连接件上。In a preferred solution, a probe connector is provided on the side of the mounting plate A away from the mounting plate B, and the ultrasonic linear array probe is fixedly installed on the probe connector.

一种优选方案,在超声线性阵列探头的至少一侧设有激光测距仪。In a preferred solution, a laser rangefinder is provided on at least one side of the ultrasonic linear array probe.

一种实施方案,所述双轴位移调节机构包括安装板C、纵向位移调节机构和横向位移调节机构,所述自适应伸缩机构与安装板C固定连接,所述安装板C与纵向位移调节机构滑动连接,所述纵向位移调节机构与横向位移调节机构滑动连接,所述横向位移调节机构固定连接在轨检小车的前端。In one embodiment, the biaxial displacement adjustment mechanism includes an installation plate C, a longitudinal displacement adjustment mechanism and a lateral displacement adjustment mechanism, the adaptive telescopic mechanism is fixedly connected with the installation plate C, and the installation plate C is connected with the longitudinal displacement adjustment mechanism. Sliding connection, the longitudinal displacement adjusting mechanism is slidably connected with the lateral displacement adjusting mechanism, and the lateral displacement adjusting mechanism is fixedly connected to the front end of the rail inspection trolley.

一种优选方案,所述纵向位移调节机构包括纵向支架和纵向电动丝杆调节机构,所述横向位移调节机构包括横向支架和横向电动丝杆调节机构,所述安装板C与纵向支架上下滑动连接,所述纵向支架与横向支架横向滑动连接,所述横向支架固定连接在轨检小车的前端。In a preferred solution, the longitudinal displacement adjustment mechanism includes a longitudinal bracket and a longitudinal electric screw adjustment mechanism, the lateral displacement adjustment mechanism includes a transverse bracket and a transverse electric screw adjustment mechanism, and the mounting plate C is slidably connected to the longitudinal bracket up and down. , the longitudinal support and the transverse support are connected in a lateral sliding manner, and the transverse support is fixedly connected to the front end of the rail inspection trolley.

一种优选方案,在轨检小车的前端还设有照明灯。In a preferred solution, a lighting lamp is also provided at the front end of the rail inspection trolley.

一种优选方案,在轨检小车的前端还设有导向机构,构成导向机构的导向轮与钢轨滚动连接。In a preferred solution, a guide mechanism is also provided at the front end of the rail inspection trolley, and the guide wheels constituting the guide mechanism are connected with the rail in a rolling manner.

一种优选方案,在轨检小车上还设有报警模块,所述报警模块与计算机信号连接。In a preferred solution, an alarm module is also provided on the rail inspection trolley, and the alarm module is signal-connected to the computer.

一种优选方案,在轨检小车上还设有移动电源。In a preferred solution, a mobile power supply is also provided on the rail inspection trolley.

与现有技术相比,本发明具有如下有益技术效果:Compared with the prior art, the present invention has the following beneficial technical effects:

1)本发明先通过摄像仪可以对轨道板表面的裂缝缺陷进行拍照检测,然后通过超声波发射接收器和超声线性阵列探头对裂缝内部缺陷进一步检测,从而可以实现轨道板裂缝的自动检测;1) In the present invention, the crack defects on the surface of the track plate can be photographed and detected by a camera first, and then the internal defects of the crack can be further detected by the ultrasonic transmitter receiver and the ultrasonic linear array probe, so that the automatic detection of the track plate crack can be realized;

2)本发明采用的超声线性阵列探头具有能量高、衰减小的优势,超声能量几乎不会泄露,更好地入射到待测材料的内部并与缺陷相互作用,形成对检测有利的散射信号,且每个探头的阵元可以单独进行激发,也可以组合激发,能够实现聚焦和偏转功能,可以更好的采集扩散场信号;2) The ultrasonic linear array probe used in the present invention has the advantages of high energy and low attenuation, the ultrasonic energy hardly leaks, and is better incident on the inside of the material to be tested and interacts with the defects, forming a scattering signal favorable for detection, And the array elements of each probe can be excited individually or in combination, which can realize focusing and deflection functions, and can better collect diffuse field signals;

3)本发明对信号进行处理及成像时,计算机先运用接收阵元扩散场信号之间的互相关运算重建格林函数的原理,获得阵元之间未延时的响应,提取出被噪声淹没的早期缺陷散射信息,然后运用频域合成孔径聚焦成像技术,不仅能确保具有高分辨率的能力,而且还有效降低成像的计算时间,速度快、精度高;3) When the present invention processes and image the signal, the computer first uses the principle of reconstructing the Green's function through the cross-correlation operation between the diffused field signals of the receiving array elements, obtains the un-delayed response between the array elements, and extracts the signal submerged by the noise. Early defect scattering information, and then using frequency-domain synthetic aperture focusing imaging technology can not only ensure the ability to have high resolution, but also effectively reduce the computing time of imaging, with high speed and high precision;

综上所述,本发明可实现高效、无损、准确、实时检测无砟轨道裂缝表面及裂缝内部缺陷,能为高铁的安全运营提供及时预警和有力保障,可为后续轨道维修工作提供有力支撑;因此,本发明相对于现有技术,具有显著进步性和应用价值。To sum up, the present invention can realize efficient, non-destructive, accurate and real-time detection of ballastless track crack surface and crack internal defects, can provide timely early warning and strong guarantee for the safe operation of high-speed rail, and can provide strong support for subsequent track maintenance work; Therefore, compared with the prior art, the present invention has significant progress and application value.

附图说明Description of drawings

图1是本发明实施例中提供的轨道板裂缝自动检测的成像方法及装置的结构示意图;1 is a schematic structural diagram of an imaging method and device for automatic detection of track plate cracks provided in an embodiment of the present invention;

图2是本发明实施例中提供的成像装置用于检测工作时的状态图;FIG. 2 is a state diagram of the imaging device provided in the embodiment of the present invention when it is used for detection;

图3是本发明实施例中提供的自适应伸缩机构的结构示意图;3 is a schematic structural diagram of an adaptive telescopic mechanism provided in an embodiment of the present invention;

图中标号示意如下:1、轨检小车;2、摄像仪;3、超声波发射接收器;4、超声线性阵列探头;5、自适应伸缩机构;51、安装板A;52、安装板B;53、伸缩支架;54、探头连接件;55、激光测距仪;6、双轴位移调节机构;61、安装板C;62、纵向位移调节机构;621、纵向支架;622、纵向电动丝杆调节机构;63、横向位移调节机构;631、横向支架;632、横向电动丝杆调节机构;7、轨道板;8、照明灯;9、导向机构;91、导向轮;10、钢轨;11、报警模块;12、移动电源;13、座椅;14、计算机。The symbols in the figure are as follows: 1. Track inspection trolley; 2. Camera; 3. Ultrasonic transmitter and receiver; 4. Ultrasonic linear array probe; 5. Self-adaptive telescopic mechanism; 51. Mounting plate A; 53. Telescopic bracket; 54. Probe connector; 55. Laser distance meter; 6. Biaxial displacement adjustment mechanism; 61. Mounting plate C; 62. Longitudinal displacement adjustment mechanism; 621. Longitudinal support; Adjustment mechanism; 63, lateral displacement adjustment mechanism; 631, lateral support; 632, lateral electric screw adjustment mechanism; 7, track plate; 8, lighting lamp; 9, guide mechanism; 91, guide wheel; 10, steel rail; 11, Alarm module; 12, mobile power supply; 13, seat; 14, computer.

具体实施方式Detailed ways

以下结合附图和实施例对本发明的技术方案做进一步详细描述。The technical solutions of the present invention will be described in further detail below with reference to the accompanying drawings and embodiments.

实施例Example

结合图1至图3所示:本发明提供的一种实现轨道板裂缝自动检测的成像装置,包括轨检小车1及设置在轨检小车1上的摄像仪2和计算机14,所述摄像仪2与计算机14信号连接,还包括超声波发射接收器3和超声线性阵列探头4,所述超声线性阵列探头4的发射端与超声波发射接收器3的发射接口信号连接,所述超声线性阵列探头4的接收端与超声波发射接收器3的接收接口信号连接,所述超声波发射接收器3的接收接口与计算机14信号连接,所述超声线性阵列探头4与一自适应伸缩机构5相连接,所述自适应伸缩机构5与一双轴位移调节机构6相连接,且所述双轴位移调节机构6固定连接在轨检小车1的前端。1 to 3 : an imaging device for realizing automatic detection of track slab cracks provided by the present invention includes a track inspection trolley 1 and a camera 2 and a computer 14 arranged on the track inspection trolley 1. The camera 2 is signal-connected with the computer 14, and also includes an ultrasonic transmitter receiver 3 and an ultrasonic linear array probe 4, the transmitting end of the ultrasonic linear array probe 4 is signal-connected with the transmitter interface of the ultrasonic transmitter receiver 3, and the ultrasonic linear array probe 4 The receiving end of the ultrasonic transmitter and receiver 3 is signal-connected with the receiving interface of the ultrasonic transmitter receiver 3, the receiver interface of the ultrasonic transmitter receiver 3 is signal-connected with the computer 14, and the ultrasonic linear array probe 4 is connected with an adaptive telescopic mechanism 5. The adaptive telescopic mechanism 5 is connected with a biaxial displacement adjusting mechanism 6 , and the biaxial displacement adjusting mechanism 6 is fixedly connected to the front end of the rail inspection trolley 1 .

本实施例中,所述超声线性阵列探头4采用市售产品即可,例如,可采用接触式的超声线性阵列探头,具有能量高,衰减小的优势,超声能量几乎不会泄露,更好地入射到待测材料的内部并与缺陷相互作用,形成对检测有利的散射信号,所述超声线性阵列探头4是由若干超声探头组成,例如,可以由8行3列共24个超声探头组成的阵列探头,相应的,所述超声线性阵列探头4即具有24个阵元,超声探头的中心频率可以为1.0MHz,每个探头的阵元可以单独进行激发,也可以组合激发,能够实现聚焦和偏转功能。此外,从上述可见,本实施例中的超声线性阵列探头4具有发射端和接收端,具备自发自收的功能。由于超声波在密度不同的介质中传播速度不同,当其通过两种不同介质的分界面时,会发生反射散射等现象,因无砟轨道属于混凝土构件,是由砂、水泥、石子等混合组成的多孔质非均匀复合材料,超声波发射接收器3发射的超声脉冲波在混凝土中传播时遇到缺陷会发生绕射,在缺陷界面会发生散射和反射,导致到达接收探头时声波能量得幅值显著减小,本申请中,采用超声线性阵列探头4,不仅能够更加有效率的采集扩散场信号,还可以使得其接收的扩散场信号叠加,使得反射扩散场信号在某一个方向的辐射能量最大,而在其它方向的总辐射能量较小,从而实现对扩散场信号的聚焦作用,使得有用的扩散场信号得到增强,进而使得干扰信号得到抑制,从而保证了信号采集及检测结果的精度。In this embodiment, the ultrasonic linear array probe 4 may be a commercially available product. For example, a contact-type ultrasonic linear array probe can be used, which has the advantages of high energy and low attenuation, and the ultrasonic energy hardly leaks. Incident to the interior of the material to be tested and interact with defects to form scattering signals favorable for detection, the ultrasonic linear array probe 4 is composed of several ultrasonic probes, for example, it can be composed of 24 ultrasonic probes in 8 rows and 3 columns. Array probe, correspondingly, the ultrasonic linear array probe 4 has 24 array elements, the center frequency of the ultrasonic probe can be 1.0MHz, and the array elements of each probe can be excited individually or in combination, which can realize focusing and deflection function. In addition, it can be seen from the above that the ultrasonic linear array probe 4 in this embodiment has a transmitting end and a receiving end, and has the function of self-transmitting and self-receiving. Due to the different propagation speed of ultrasonic waves in media with different densities, when it passes through the interface of two different media, reflection and scattering will occur. Because the ballastless track is a concrete component, it is composed of a mixture of sand, cement, and stones. Porous heterogeneous composite materials, the ultrasonic pulse wave emitted by the ultrasonic transmitter and receiver 3 will be diffracted when it encounters defects in the concrete, and will be scattered and reflected at the defect interface, resulting in a significant amplitude of the sound wave energy when it reaches the receiving probe. In this application, using the ultrasonic linear array probe 4 can not only collect the diffuse field signals more efficiently, but also superimpose the diffuse field signals received by the ultrasonic linear array probe 4, so that the radiated energy of the reflected diffuse field signals in a certain direction is the largest, The total radiation energy in other directions is smaller, so as to realize the focusing effect on the diffused field signal, so that the useful diffused field signal is enhanced, and then the interference signal is suppressed, thereby ensuring the accuracy of signal acquisition and detection results.

参见图1至图3所示,本实施例中,所述自适应伸缩机构5包括垂向设置的、互相平行的安装板A51和安装板B52,所述安装板A51和安装板B52之间连接有可水平伸缩的伸缩支架53,所述超声线性阵列探头4与安装板A51相连接,所述双轴位移调节机构6与安装板B52相连接,自适应伸缩机构5与双轴位移调节机构6相结合,可以实现超声线性阵列探头4三轴位移调节,从而可根据要求灵活调节超声线性阵列探头4的位置,进而灵活调节以超声线性阵列探头4每次扫描检测的范围。所述伸缩支架53可以由若干伸缩杆(图中未标记)交叉连接组成,伸缩支架53连接有伸缩驱动机构(未显示),通过伸缩驱动机构控制伸缩支架53的伸缩,伸缩支架53、伸缩驱动机构及自适应伸缩机构5的整体均属于现有技术,此处就不再一一赘述。Referring to FIG. 1 to FIG. 3 , in this embodiment, the adaptive telescopic mechanism 5 includes a mounting plate A51 and a mounting plate B52 that are vertically arranged and parallel to each other, and the mounting plate A51 and the mounting plate B52 are connected to each other. There is a horizontally retractable telescopic bracket 53, the ultrasonic linear array probe 4 is connected with the mounting plate A51, the biaxial displacement adjusting mechanism 6 is connected with the mounting plate B52, and the adaptive telescopic mechanism 5 is connected with the biaxial displacement adjusting mechanism 6 In combination, the three-axis displacement adjustment of the ultrasonic linear array probe 4 can be realized, so that the position of the ultrasonic linear array probe 4 can be flexibly adjusted according to requirements, and then the detection range of the ultrasonic linear array probe 4 in each scan can be flexibly adjusted. The telescopic support 53 can be composed of a plurality of telescopic rods (not marked in the figure) that are cross-connected. The telescopic support 53 is connected with a telescopic drive mechanism (not shown), and the telescopic drive mechanism controls the telescopic support 53. The entire mechanism and the self-adaptive telescopic mechanism 5 belong to the prior art, and will not be repeated here.

此外,所述安装板A51远离安装板B52的一侧设有探头连接件54,所述超声线性阵列探头4固定安装在探头连接件54上,使用的时候,将超声线性阵列探头4固定安装在探头连接件54上并使超声线性阵列探头4的发射接收端面始终与无砟轨道表面保持垂直。In addition, the side of the mounting plate A51 away from the mounting plate B52 is provided with a probe connector 54, and the ultrasonic linear array probe 4 is fixedly installed on the probe connector 54. When in use, the ultrasonic linear array probe 4 is fixedly installed on the On the probe connector 54, the transmitting and receiving end faces of the ultrasonic linear array probe 4 are always kept perpendicular to the surface of the ballastless track.

此外,在超声线性阵列探头4的至少一侧设有激光测距仪55,通过激光测距仪55所检测的位移,即可判断超声线性阵列探头4的发射接收端面与无砟轨道板7的表面之间的垂直距离,以保持垂直距离在预设的范围内,进而防止沉降距离过大而损伤探头。本实施例中,在超声线性阵列探头的左右两侧分别对称设有激光测距仪55,两个激光测距仪55可以设定为一用一备。In addition, a laser range finder 55 is provided on at least one side of the ultrasonic linear array probe 4, and the displacement detected by the laser range finder 55 can be used to determine the distance between the transmitting and receiving end face of the ultrasonic linear array probe 4 and the ballastless track plate 7. The vertical distance between the surfaces to keep the vertical distance within the preset range, thereby preventing the probe from being damaged due to the excessive settlement distance. In this embodiment, the left and right sides of the ultrasonic linear array probe are respectively provided with laser rangefinders 55 symmetrically, and the two laser rangefinders 55 can be set as one for use and one for backup.

参见图1和图2所示,所述双轴位移调节机构6包括安装板C61、纵向位移调节机构62和横向位移调节机构63,所述自适应伸缩机构5与安装板C61固定连接(具体的,自适应伸缩机构5中的安装板B52与安装板C61固定连接,从而使得自适应伸缩机构5与双轴位移调节机构6相连接),所述安装板C61与纵向位移调节机构62滑动连接,所述纵向位移调节机构62与横向位移调节机构63滑动连接,所述横向位移调节机构63固定连接在轨检小车1的前端。1 and 2, the biaxial displacement adjusting mechanism 6 includes a mounting plate C61, a longitudinal displacement adjusting mechanism 62 and a lateral displacement adjusting mechanism 63, and the adaptive telescopic mechanism 5 is fixedly connected to the mounting plate C61 (specifically , the mounting plate B52 in the self-adaptive telescopic mechanism 5 is fixedly connected with the mounting plate C61, so that the self-adaptive telescopic mechanism 5 is connected with the biaxial displacement adjusting mechanism 6), and the mounting plate C61 is slidably connected with the longitudinal displacement adjusting mechanism 62, The longitudinal displacement adjusting mechanism 62 is slidably connected with the lateral displacement adjusting mechanism 63 , and the lateral displacement adjusting mechanism 63 is fixedly connected to the front end of the rail inspection trolley 1 .

具体的,所述纵向位移调节机构62包括纵向支架621和纵向电动丝杆调节机构622,所述横向位移调节机构63包括横向支架631和横向电动丝杆调节机构632,所述安装板C61与纵向支架621上下滑动连接,所述纵向支架621与横向支架631横向滑动连接,所述横向支架631固定连接在轨检小车1的前端。本实施例中,所述纵向电动丝杆调节机构622和横向电动丝杆调节机构632均是由驱动电机、一端固定在驱动电机输出端的丝杆和与丝杆螺纹连接的滑块连接件组成,由于此组成结构为已知技术,故在图中未详细示出。Specifically, the longitudinal displacement adjustment mechanism 62 includes a longitudinal bracket 621 and a longitudinal electric screw adjustment mechanism 622, the lateral displacement adjustment mechanism 63 includes a lateral bracket 631 and a transverse electric screw adjustment mechanism 632, and the mounting plate C61 is connected to the longitudinal The bracket 621 is slidably connected up and down, and the longitudinal bracket 621 is slidably connected to the transverse bracket 631 , and the transverse bracket 631 is fixedly connected to the front end of the rail inspection trolley 1 . In this embodiment, the longitudinal electric screw adjustment mechanism 622 and the transverse electric screw adjustment mechanism 632 are both composed of a drive motor, a screw with one end fixed at the output end of the drive motor, and a slider connector threadedly connected to the screw. Since this composition is a known technology, it is not shown in detail in the figure.

此外,参见图1和图2所示,在轨检小车1的前端还设有照明灯8,以避免在夜间或较暗环境下摄像仪2拍摄模糊无法识别裂缝的问题,避免漏检,增加了检测效率。照明灯8的数量根据需要灵活设定,例如,本实施例中,是在轨检小车1的前端左右对称设有两个照明灯8,In addition, as shown in FIG. 1 and FIG. 2 , the front end of the on-track inspection trolley 1 is also provided with a light 8 to avoid the problem of blurred and unrecognized cracks when the camera 2 is photographed at night or in a dark environment, avoiding missed inspections and increasing detection efficiency. The number of the lighting lamps 8 can be set flexibly according to the needs. For example, in this embodiment, two lighting lamps 8 are symmetrically arranged at the front end of the rail inspection trolley 1.

此外,在轨检小车1的前端还设有导向机构9,构成导向机构9的导向轮91与钢轨10滚动连接,以对轨检小车1的行驶方向起到导向作用。In addition, a guide mechanism 9 is also provided at the front end of the rail inspection trolley 1 , and the guide wheels 91 constituting the guide mechanism 9 are rollingly connected with the rail 10 to guide the running direction of the rail inspection trolley 1 .

此外,在轨检小车1上还设有报警模块11,所述报警模块11与计算机14信号连接,以便发现轨道上有裂缝存在时,可以通过计算机14发出指令以驱动报警模块11进行报警。所述报警模块11采用市售报警器即可。In addition, an alarm module 11 is also provided on the rail inspection trolley 1, and the alarm module 11 is connected with the computer 14 by signal, so that when a crack is found on the track, the computer 14 can issue an instruction to drive the alarm module 11 to give an alarm. The alarm module 11 can be a commercially available alarm device.

此外,在轨检小车1上还设有移动电源12,以实现对所述装置进行移动供电。In addition, a mobile power supply 12 is also provided on the rail inspection trolley 1 to implement mobile power supply for the device.

此外,在轨检小车1上设置座椅13,以便由检测员坐在轨检小车1上进行人工信号采集及检测操作。In addition, a seat 13 is provided on the rail inspection trolley 1 so that the inspector can sit on the rail inspection trolley 1 to perform manual signal collection and detection operations.

采用本发明所述成像装置实现轨道板裂缝自动检测的成像方法如下:The imaging method for realizing automatic detection of track slab cracks using the imaging device of the present invention is as follows:

a)使轨检小车1沿着轨道板7方向行驶,行驶时,摄像仪2对轨道板7表面进行拍摄,并将拍摄的图片传递给计算机14,计算机14根据预存的正常轨道板图片与之比对,以判断轨道板7表面是否存在裂缝;a) Make the track inspection trolley 1 run along the direction of the track plate 7. When driving, the camera 2 takes a picture of the surface of the track plate 7, and transmits the photographed picture to the computer 14, and the computer 14 matches it with the pre-stored normal track plate picture Compare to judge whether there are cracks on the surface of the track plate 7;

b)当判断结果显示存在裂缝时,则控制器控制报警模块11进行报警,轨检小车1停止行驶,调节超声线性阵列探头4的位置使其下沉压置于裂缝处,然后启动超声波发射接收器3,使超声波发射接收器3的发射端发射频率在1~5MHz范围内的超声波信号,并由超声线性阵列探头4的发射端发射至轨道板7内部,超声波信号在轨道板7内部传播过程中传播在缺陷界面发生散射和反射从而形成扩散场信号,然后由超声线性阵列探头4的接收端采集扩散场信号并将其传输给超声波发射接收器3的接收端,然后由超声波发射接收器3的接收端将接收的扩散场信号传输给计算机14;即:b) When the judgment result shows that there is a crack, the controller controls the alarm module 11 to give an alarm, the rail inspection trolley 1 stops running, adjusts the position of the ultrasonic linear array probe 4 to make it sink and press at the crack, and then starts the ultrasonic transmission and reception 3, so that the transmitting end of the ultrasonic transmitter and receiver 3 transmits ultrasonic signals with a frequency in the range of 1 to 5 MHz, and is transmitted from the transmitting end of the ultrasonic linear array probe 4 to the inside of the track plate 7, and the ultrasonic signal propagates inside the track plate 7. In the medium propagation, scattering and reflection occur at the defect interface to form a diffuse field signal, and then the diffuse field signal is collected by the receiving end of the ultrasonic linear array probe 4 and transmitted to the receiving end of the ultrasonic transmitter receiver 3, and then the ultrasonic transmitter receiver 3. The receiving end of the device transmits the received diffuse field signal to the computer 14; namely:

c)由计算机采用MATLAB软件对所接收到的信号被动提取格林函数,对扩散场信号进行互相关,重建阵元之间的格林函数,获取阵元之间未延时的响应,恢复被噪声淹没的早期缺陷信息,然后依据频域合成孔径聚焦成像算法进行轨道板裂缝成像,具体为:c) The computer uses MATLAB software to passively extract the Green's function from the received signal, cross-correlate the diffuse field signal, reconstruct the Green's function between the array elements, obtain the undelayed response between the array elements, and recover from being submerged by noise. Then, the track plate crack imaging is performed according to the frequency domain synthetic aperture focusing imaging algorithm, which is as follows:

1)数据处理:1) Data processing:

首先,由于超声波信号在无砟轨道板内部传播的过程中,超声波信号和缺陷目标之间发生相互作用,经过散射和多次的反射后会形成一个近似均匀的声场,该声场即为扩散场,而如步骤b)所示,通过超声线性阵列探头即可获取无砟轨道中的扩散场信号,设定超声线性阵列探头中任意的两个接收阵元分别为γ1和γ2,两个接收阵元γ1和γ2处于封闭的曲面空间内,则接收阵元γ1和γ2处声场互谱的所有声源积分等于γ1和γ2之间频域因果格林函数与非因果格林函数之差:First of all, due to the interaction between the ultrasonic signal and the defect target during the propagation of the ultrasonic signal inside the ballastless track slab, an approximately uniform sound field will be formed after scattering and multiple reflections, which is the diffuse field. As shown in step b), the diffused field signal in the ballastless track can be obtained through the ultrasonic linear array probe, and any two receiving array elements in the ultrasonic linear array probe are set to be γ 1 and γ 2 respectively, and the two receiving elements are If the array elements γ 1 and γ 2 are in a closed surface space, then the integral of all sound sources of the sound field cross-spectrum at the receiving array elements γ 1 and γ 2 is equal to the frequency domain causal Green function and non-causal Green function between γ 1 and γ 2 Difference:

G(r1,r2,ω)-G*(r1,r2,ω)=-2jω∫G(r1,r,ω)G*(r2,r,ω)dV (1);G(r 1 , r 2 , ω)-G * (r 1 , r 2 , ω)=-2jω∫G(r 1 , r, ω)G * (r 2 , r, ω)dV (1);

公式(1)中,等式的左侧是两个接收阵元γ1和γ2之间的频域格林函数的因果和非因果函数之差,等式右侧表示接收阵元γ1和γ2处声场互相关的所有声源积分,具体的,G(γ1,γ2,ω)代表频域中γ1和γ2之间的因果格林函数,上角标*代表共轭复数,G*1,γ2,ω)代表频域中γ1和γ2之间的非因果格林函数,j代表虚数单位,j2=-1,ω代表发射信号的角频率,γ是处于扩散场积分密闭空间V中的任意一点位置,代表任意的缺陷目标,可视为噪声源,G(γ1,γ,ω)代表γ和γ1之间的频域格林函数传播公式,G(γ2,γ,ω)代表γ和γ2之间的频域格林函数传播公式,dV代表密闭空间V的微分;In formula (1), the left side of the equation is the difference between the causal and non-causal functions of the frequency domain Green function between the two receiving elements γ 1 and γ 2 , and the right side of the equation represents the receiving elements γ 1 and γ Integral of all sound sources for cross-correlation of sound fields at 2 places, specifically, G(γ 1 , γ 2 , ω) represents the causal Green function between γ 1 and γ 2 in the frequency domain, the superscript * represents the conjugate complex number, G *1 , γ 2 , ω) represents the non-causal Green's function between γ 1 and γ 2 in the frequency domain, j represents the imaginary unit, j 2 =-1, ω represents the angular frequency of the transmitted signal, γ is in the diffusion The position of any point in the closed space V of field integration represents any defect target and can be regarded as a noise source, G(γ 1 , γ, ω) represents the frequency domain Green’s function propagation formula between γ and γ 1 , G(γ 2 , γ, ω) represents the frequency domain Green's function propagation formula between γ and γ 2 , dV represents the differential of the closed space V;

然后,假设扩散场噪声源强度是空间均匀分布且互不相关,则功率谱密度|q(ω)|2同时与噪声源的位置无关:Then, the power spectral density |q(ω)| 2 is also independent of the position of the noise source, assuming that the intensity of the diffuse-field noise source is spatially uniformly distributed and uncorrelated with each other:

<q(r1,ω)q*(r2,ω)>=δ(r1-r2)|q(ω)|2 (2);<q(r 1 , ω)q * (r 2 , ω)>=δ(r 1 -r 2 )|q(ω)| 2 (2);

公式(2)代表频域γ1和γ2处信号的互功率谱的期望值等于γ1和γ2间噪声源的功率谱密度,具体的,公式(2)中,<q(γ1,ω)q*2,ω)>代表频域γ1和γ2处信号的互功率谱的期望值,q(γ1,ω)代表γ1处的场强,q(γ2,ω)代表γ2处的场强,<>代表统计平均的期望值,上角标*代表共轭复数,δ(γ1-γ2)代表γ1和γ2之间信号的互功率谱与噪声源的位置无关;频域中两个接收阵元γ1和γ2位置处声场之间的互相关为:Equation (2) represents the expected value of the cross-power spectrum of the signals at γ 1 and γ 2 in the frequency domain equal to the power spectral density of the noise source between γ 1 and γ 2. Specifically, in formula (2), <q(γ 1 , ω )q *2 , ω)> represents the expected value of the cross-power spectrum of the signals at γ 1 and γ 2 in the frequency domain, q(γ 1 , ω) represents the field strength at γ 1 , and q(γ 2 , ω) represents the The field strength at γ 2 , <> represents the expected value of the statistical average, the superscript * represents the conjugate complex number, δ(γ1-γ2) represents the cross-power spectrum of the signal between γ 1 and γ 2 is independent of the position of the noise source; The cross-correlation between the sound fields at the positions of the two receiving array elements γ 1 and γ 2 in the frequency domain is:

<p(r1,ω)p*(r2,ω)>=|q(ω)|2∫G(r1,r,ω)G*(r2,r,ω)dV (3);<p(r 1 , ω)p * (r 2 , ω)>=|q(ω)| 2 ∫G(r 1 , r, ω)G * (r 2 , r, ω)dV (3);

公式(3)中,p(γ1,ω)代表接收阵元γ1处的声场,p(γ2,ω)代表接收阵元γ2处的声场,<p(γ1,ω)p*2,ω)>代表两个接收阵元γ1和γ2位置处声场之间的互相关函数,上角标*代表共轭复数;In formula (3), p(γ 1 , ω) represents the sound field at the receiving element γ 1 , p(γ 2 , ω) represents the sound field at the receiving element γ 2 , <p(γ 1 , ω)p *2 , ω)> represents the cross-correlation function between the sound fields at the positions of the two receiving array elements γ 1 and γ 2 , and the superscript * represents a complex conjugate number;

然后,由公式(1)和(3)得到:Then, from formulas (1) and (3), we get:

(G(r1,r2,ω)-G*(r1,r2,ω))|q(ω)|2=-2jω<p(r1,ω)p*(r2,ω)> (4):(G(r 1 , r 2 , ω)-G * (r 1 , r 2 , ω))|q(ω)| 2 =-2jω<p(r 1 , ω)p * (r 2 , ω) > (4):

公式(4)中,等式左侧为格林函数G(γ12,ω)以及它对应的时间反转,即对应频域中的共轭操作,把这两项乘以随机噪声的功率谱密度后,其结果等于等式右侧扩散场中两个的接收阵元γ1和γ2互相关结果,上角标*代表共轭复数;In formula (4), the left side of the equation is the Green's function G(γ 1 , γ 2 , ω) and its corresponding time reversal, that is, the conjugate operation in the corresponding frequency domain. Multiply these two terms by the random noise. After the power spectral density, the result is equal to the cross-correlation result of the two receiving array elements γ 1 and γ 2 in the diffuse field on the right side of the equation, and the superscript * represents the conjugate complex number;

然后,求出公式(4)所对应的时域表达式,频域2jω对应时域2d/dt,根据卷积定理可知,频域中的乘积对应时域中的卷积,从而得到:Then, the time domain expression corresponding to formula (4) is obtained. The frequency domain 2jω corresponds to the time domain 2d/dt. According to the convolution theorem, the product in the frequency domain corresponds to the convolution in the time domain, thus obtaining:

Figure BDA0002444267950000121
Figure BDA0002444267950000121

公式(5)中,G(γ12,t)代表时域中γ1和γ2之间的格林函数,G(γ12,-t)代表G(γ12,t)的时间反转,即对应频域中的共轭操作,*代表卷积运算,

Figure BDA0002444267950000122
代表互相关运算,Cq(t)表示扩散场中的噪声q(t)的自相关结果,d/dt代表对t求导,p(γ1,t)代表时域中接收阵元γ1处的声场,p(γ2,t)代表时域中接收阵元γ2处的声场;In formula (5), G(γ 12 ,t) represents the Green’s function between γ 1 and γ 2 in the time domain, and G(γ 12 ,-t) represents G(γ 12 , t) time reversal, which corresponds to the conjugate operation in the frequency domain, * represents the convolution operation,
Figure BDA0002444267950000122
represents the cross-correlation operation, C q (t) represents the autocorrelation result of the noise q(t) in the diffuse field, d/dt represents the derivation of t, p(γ 1 ,t) represents the receiving array element γ 1 in the time domain The sound field at , p(γ 2 ,t) represents the sound field at the receiving array element γ 2 in the time domain;

公式(4)和(5)的结果表明扩散场中接收阵元声场p(γ1,t)和p(γ2,t)互相关并求导,其结果等于这两接收阵元之间的格林函数响应,根据声波的互易性原理,在时间轴上,这两接收阵元之间的格林函数响应具有对称性;The results of formulas (4) and (5) show that the sound fields p(γ 1 ,t) and p(γ 2 ,t) of the receiving array elements in the diffuse field are correlated and derived, and the result is equal to the difference between the two receiving array elements. Green's function response, according to the principle of reciprocity of sound waves, on the time axis, the Green's function response between the two receiving array elements has symmetry;

2)频域合成孔径聚焦(SAFT)成像:2) Frequency Domain Synthetic Aperture Focusing (SAFT) imaging:

频域合成孔径聚焦利用的是相位延迟方法,在波动方程反演理论的基础上对格林函数进行傅立叶变换分解,所以又称为波数算法,之后利用Stolt Mapping和插值实现坐标变换,滤波后处理后再进行二维傅立叶反变换最终生成空间-时间下的检测图像,因此,可以对缺陷目标使用爆炸模式的思想,把一个缺陷目标离散成独立的散射点,这些散射点被视为二次声源,总的散射声场大小是各个独立的散射点声场的叠加,然后根据每个阵元的接收信号即可反推缺陷的具体位置和大小;自发自收的传播路径(声程)正好等于该阵元到缺陷目标距离的2倍;The frequency-domain synthetic aperture focusing uses the phase delay method. The Green's function is decomposed by Fourier transform on the basis of the wave equation inversion theory, so it is also called the wavenumber algorithm. After that, the coordinate transformation is realized by using Stolt Mapping and interpolation. After filtering and post-processing Then the two-dimensional inverse Fourier transform is performed to finally generate the detection image in space-time. Therefore, the idea of explosion mode can be used for the defect target, and a defect target can be discretized into independent scattering points, which are regarded as secondary sound sources. , the total scattered sound field size is the superposition of the sound fields of each independent scattering point, and then the specific position and size of the defect can be reversed according to the received signal of each array element; the propagation path (sound path) of spontaneous and self-receiving is exactly equal to the array 2 times the distance from the element to the defect target;

首先,由于超声波在无限的均匀介质传播满足Helmholtz波动方程,因此,假设在二维空间x-z坐标轴下,超声波发射阵元的坐标为(u,0),接收阵元的坐标为(v,0),缺陷目标距离发射阵元和接收阵元的距离分别为ρin和ρout,任意一对发射接收阵元对应的时域接收信号为e(v,u,t),那么对应的发射接收阵元所对应的频率响应E(v,u,ω)为:First of all, since the propagation of ultrasonic waves in an infinite homogeneous medium satisfies the Helmholtz wave equation, it is assumed that the coordinates of the ultrasonic transmitting array element are (u, 0) and the coordinates of the receiving array element are (v, 0 under the xz coordinate axis of the two-dimensional space). ), the distances of the defect target from the transmitting array element and the receiving array element are ρ in and ρ out respectively, and the time domain receiving signal corresponding to any pair of transmitting and receiving array elements is e(v, u, t), then the corresponding transmitting and receiving The frequency response E(v, u, ω) corresponding to the array element is:

E(v,u,ω)=P(ω)∫∫f(x,z)g(ρin,ω)g(ρout,ω)dxdz (6);E(v, u, ω)=P(ω)∫∫f(x, z)g(ρ in ,ω)g(ρ out ,ω)dxdz (6);

公式(6)中,

Figure BDA0002444267950000131
P(ω)是发射信号的频谱,ω为发射信号的角频率;f(x,z)代表缺陷目标的点扩散函数,g(ρin,ω)代表发射阵元到缺陷目标的空间格林函数频域响应,g(ρout,ω)代表接收阵元到缺陷目标的空间格林函数频域响应;In formula (6),
Figure BDA0002444267950000131
P(ω) is the spectrum of the transmitted signal, ω is the angular frequency of the transmitted signal; f(x, z) represents the point spread function of the defective target, and g(ρ in , ω) represents the spatial Green's function from the transmitting array element to the defective target Frequency domain response, g(ρ out , ω) represents the spatial Green's function frequency domain response from the receiving array element to the defect target;

此外,格林函数为:Also, the Green's function is:

Figure BDA0002444267950000132
Figure BDA0002444267950000132

公式(7)中,g(x,z,ω)代表二维空间x-z坐标轴下任意点的格林函数,波数k=ω/c,c为声波在介质中的纵波速度,空间变量x和z与波数kx和kz相对应;In formula (7), g(x, z, ω) represents the Green’s function of any point under the xz coordinate axis of the two-dimensional space, the wave number k=ω/c, c is the longitudinal wave velocity of the sound wave in the medium, and the spatial variables x and z Corresponds to the wave numbers k x and k z ;

使用变量x和z代替掉变量ρin和ρout,即g(x,z,ω)代替g(ρin、ρout,ω),将公式(7)带入公式(6),重新获得接收响应E(v,u,ω),然后用波数ku和kv分别表示发射阵元(u,0)和接收阵元(v,0)的空间信息,为了实现将积分转换二维傅里叶变换,对接收响应E(v,u,ω)中变量u和v进行傅里叶变换得到:Use the variables x and z to replace the variables ρ in and ρ out , that is, g(x, z, ω) instead of g(ρ in , ρ out , ω), bring the formula (7) into the formula (6), and regain the reception Response E(v, u, ω), and then use the wave numbers ku and k v to represent the spatial information of the transmitting array element (u, 0) and the receiving array element (v, 0), respectively. In order to realize the transformation of the integral into a two-dimensional Fourier Leaf transform, the Fourier transform of the variables u and v in the received response E(v, u, ω) is obtained:

Figure BDA0002444267950000133
Figure BDA0002444267950000133

公式(8)中,E(ω,ku,kv)代表频域响应E(v,u,ω)的空间傅里叶变换,F代表点扩散函数f(x,z)的空间傅里叶变换;In formula (8), E(ω,k u ,k v ) represents the spatial Fourier transform of the frequency domain response E(v, u, ω), and F represents the spatial Fourier transform of the point spread function f(x, z) leaf transform;

相应的,点扩散函数f(x,z)的空间傅里叶变换是F(kx,kz),为了求得图像域,反变换模型为:Correspondingly, the spatial Fourier transform of the point spread function f(x, z) is F(k x , k z ). In order to obtain the image domain, the inverse transform model is:

Figure BDA0002444267950000141
Figure BDA0002444267950000141

公式(9)中,

Figure BDA0002444267950000142
代表缺陷目标等间隔采样的数据域模型,S-1{}代表Stolt映射,完成了从数据域到图像域的转换,p*(ω)代表发射信号的共轭频谱;In formula (9),
Figure BDA0002444267950000142
represents the data domain model of the defect target sampled at equal intervals, S -1 {} represents the Stolt map, which completes the conversion from the data domain to the image domain, p * (ω) represents the conjugate spectrum of the transmitted signal;

对公式(9)进行非线性坐标变化,将波数ku,kv和k映射到和图像域缺陷目标相关联的波数kx和kz,则:The nonlinear coordinate change of formula (9) is performed to map the wave numbers ku , k v and k to the wave numbers k x and k z associated with the defect target in the image domain, then:

kx=ku+kv (10);k x =k u +k v (10);

Figure BDA0002444267950000143
Figure BDA0002444267950000143

然后通过二维傅里叶反变换求得缺陷目标分布的图像Then the image of the defect target distribution is obtained by the two-dimensional inverse Fourier transform

Figure BDA0002444267950000144
Figure BDA0002444267950000144

公式(12)中,

Figure BDA0002444267950000145
代表缺陷目标分布的图像域模型。In formula (12),
Figure BDA0002444267950000145
An image domain model representing the distribution of defect targets.

综上所述,本发明先利用摄像仪检测轨道板裂缝表面缺陷,再利用超声波发射接收器发射频率在1~5MHz范围内的超声波信号至轨道板内部,同时通过超声线性阵列探头的自发自收模式对轨道板裂缝进行等间距扫描,用超声线性阵列探头的全矩阵捕获功能获取轨道板中的扩散场信号并将扩散场信号传输给计算机,最后由计算机MATLAB软件对扩散场信号进行处理并依据频域合成孔径聚焦成像算法进行轨道板裂缝成像,可以直观的获得轨道板裂缝表面缺陷情况及裂缝内部缺陷情况;此外,本发明采用超声线性阵列探头对扩散场信号进行采集,不仅信号采集过程中灵敏度高、分辨率高、速度快、效率高,且不易受外界条件干扰、精度高,使得后续的无损检测具有检测速度快、效率高、检测精度高等优点,可实现高效、无损、准确、实时检测无砟轨道裂缝,能为高铁的安全运营提供及时预警和有力保障,可为后续轨道维修工作提供有力支撑;因此,本发明相对于现有技术,具有显著进步性和应用价值。To sum up, the present invention first uses a camera to detect the surface defects of the track plate cracks, and then uses the ultrasonic transmitter receiver to transmit ultrasonic signals with a frequency in the range of 1 to 5MHz to the inside of the track plate. The mode scans the cracks of the track slab at equal intervals, and uses the full-matrix capture function of the ultrasonic linear array probe to obtain the diffuse field signal in the track slab and transmits the diffuse field signal to the computer. Finally, the diffuse field signal is processed by the computer MATLAB software and based on the The frequency domain synthetic aperture focusing imaging algorithm is used to image the track slab cracks, which can intuitively obtain the surface defects of the track slab cracks and the internal defects of the cracks; High sensitivity, high resolution, high speed, high efficiency, and it is not easily disturbed by external conditions and has high precision, which makes subsequent non-destructive testing have the advantages of fast detection speed, high efficiency and high detection accuracy, and can achieve efficient, non-destructive, accurate and real-time detection. The detection of ballastless track cracks can provide timely early warning and strong guarantee for the safe operation of high-speed rail, and can provide strong support for subsequent track maintenance work; therefore, compared with the prior art, the present invention has significant progress and application value.

最后有必要在此指出的是:以上所述仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。Finally, it is necessary to point out here: the above is only a preferred embodiment of the present invention, but the protection scope of the present invention is not limited to this, and any person skilled in the art is within the technical scope disclosed by the present invention. Any changes or substitutions that can be easily thought of within the scope of the present invention should be covered within the protection scope of the present invention.

Claims (10)

1.一种实现轨道板裂缝自动检测的成像方法,其特征在于,包括如下步骤:1. an imaging method that realizes automatic detection of track plate cracks, is characterized in that, comprises the steps: a)利用摄像仪检测轨道板裂缝表面缺陷;a) Use a camera to detect surface defects of track slab cracks; b)利用超声波发射接收器激发并经由超声线性阵列探头发射频率在广5MHz范围内的超声波信号至轨道板内部,由超声线性阵列探头采用自发自收的模式获取轨道板中的扩散场信号并将扩散场信号传输给计算机;b) Use the ultrasonic transmitter and receiver to excite and transmit ultrasonic signals with a frequency in a wide range of 5MHz to the inside of the track plate through the ultrasonic linear array probe, and the ultrasonic linear array probe adopts the mode of spontaneous and self-receiving to obtain the diffused field signal in the track plate and send it to the track plate. The diffuse field signal is transmitted to the computer; c)由计算机采用MATLAB软件对所接收到的信号被动提取格林函数,对扩散场信号进行互相关,重建阵元之间的格林函数,获取阵元之间未延时的响应,恢复被噪声淹没的早期缺陷信息,然后依据频域合成孔径聚焦成像算法进行轨道板裂缝成像。c) The computer uses MATLAB software to passively extract the Green's function from the received signal, cross-correlate the diffuse field signal, reconstruct the Green's function between the array elements, obtain the undelayed response between the array elements, and recover from being submerged by noise. The early defect information is obtained, and then the track plate crack imaging is carried out according to the frequency domain synthetic aperture focusing imaging algorithm. 2.根据权利要求1所述的成像方法,其特征在于,所述步骤c)具体包括如下操作:2. The imaging method according to claim 1, wherein the step c) specifically comprises the following operations: 1)数据处理:1) Data processing: 首先,由于超声波信号在无砟轨道板内部传播的过程中,超声波信号和缺陷目标之间发生相互作用,经过散射和多次的反射后会形成一个近似均匀的声场,该声场即为扩散场,而如步骤b)所示,通过超声线性阵列探头即可获取无砟轨道中的扩散场信号,设定超声线性阵列探头中任意的两个接收阵元分别为γ1和γ2,两个接收阵元γ1和γ2处于封闭的曲面空间内,则接收阵元γ1和γ2处声场互谱的所有声源积分等于γ1和γ2之间频域因果格林函数与非因果格林函数之差:First of all, due to the interaction between the ultrasonic signal and the defect target during the propagation of the ultrasonic signal inside the ballastless track slab, an approximately uniform sound field will be formed after scattering and multiple reflections, which is the diffusion field. As shown in step b), the diffused field signal in the ballastless track can be obtained through the ultrasonic linear array probe, and any two receiving array elements in the ultrasonic linear array probe are set to be γ 1 and γ 2 respectively, and the two receiving elements are If the array elements γ 1 and γ 2 are in a closed surface space, then the integral of all sound sources of the sound field cross-spectrum at the receiving array elements γ 1 and γ 2 is equal to the frequency domain causal Green function and non-causal Green function between γ 1 and γ 2 Difference: G(r1,r2,ω)-G*(r1,r2,ω)=-2jω∫G(r1,r,ω)G*(r2,r,ω)dV (1);G(r 1 , r 2 , ω)-G * (r 1 , r 2 , ω)=-2jω∫G(r 1 , r, ω)G * (r 2 , r, ω)dV (1); 公式(1)中,等式的左侧是两个接收阵元γ1和γ2之间的频域格林函数的因果和非因果函数之差,等式右侧表示接收阵元γ1和γ2处声场互相关的所有声源积分,具体的,G(γ1,γ2,ω)代表频域中γ1和γ2之间的因果格林函数,上角标*代表共轭复数,G*1,γ2,ω)代表频域中γ1和γ2之间的非因果格林函数,j代表虚数单位,j2=-1,ω代表发射信号的角频率,γ是处于扩散场积分密闭空间V中的任意一点位置,代表任意的缺陷目标,可视为噪声源,G(γ1,γ,ω)代表γ和γ1之间的频域格林函数传播公式,G(γ2,γ,ω)代表γ和γ2之间的频域格林函数传播公式,dV代表密闭空间V的微分;In formula (1), the left side of the equation is the difference between the causal and non-causal functions of the frequency domain Green function between the two receiving elements γ 1 and γ 2 , and the right side of the equation represents the receiving elements γ 1 and γ Integral of all sound sources for cross-correlation of sound fields at 2 places, specifically, G(γ 1 , γ 2 , ω) represents the causal Green function between γ 1 and γ 2 in the frequency domain, the superscript * represents the conjugate complex number, G *1 , γ 2 , ω) represents the non-causal Green's function between γ 1 and γ 2 in the frequency domain, j represents the imaginary unit, j 2 =-1, ω represents the angular frequency of the transmitted signal, γ is in the diffusion The position of any point in the closed space V of field integration represents any defect target and can be regarded as a noise source, G(γ 1 , γ, ω) represents the frequency domain Green’s function propagation formula between γ and γ 1 , G(γ 2 , γ, ω) represents the frequency domain Green's function propagation formula between γ and γ 2 , dV represents the differential of the closed space V; 然后,假设扩散场噪声源强度是空间均匀分布且互不相关,则功率谱密度|q(ω)|2同时与噪声源的位置无关:Then, the power spectral density |q(ω)| 2 is also independent of the position of the noise source, assuming that the intensity of the diffuse-field noise source is spatially uniformly distributed and uncorrelated with each other: <q(r1,ω)q*(r2,ω)>=δ(r1-r2)|q(ω)|2 (2);<q(r 1 , ω)q * (r 2 , ω)>=δ(r 1 -r 2 )|q(ω)| 2 (2); 公式(2)代表频域γ1和γ2处信号的互功率谱的期望值等于γ1和γ2间噪声源的功率谱密度,具体的,公式(2)中,<q(γ1,ω)q*2,ω/)>代表频域γ1和γ2处信号的互功率谱的期望值,q(γ1,ω)代表γ1处的场强,q(γ2,ω)代表γ2处的场强,<>代表统计平均的期望值,上角标*代表共轭复数,δ(γ1-γ2)代表γ1和γ2之间信号的互功率谱与噪声源的位置无关;频域中两个接收阵元γ1和γ2位置处声场之间的互相关为:Equation (2) represents the expected value of the cross-power spectrum of the signals at γ 1 and γ 2 in the frequency domain equal to the power spectral density of the noise source between γ 1 and γ 2. Specifically, in formula (2), <q(γ 1 , ω )q *2 , ω/)> represents the expected value of the cross-power spectrum of the signals at γ 1 and γ 2 in the frequency domain, q(γ 1 , ω) represents the field strength at γ 1 , q(γ 2 , ω) represents the field strength at γ 2 , <> represents the expected value of the statistical average, the superscript * represents the conjugate complex number, δ(γ1-γ2) represents the cross-power spectrum of the signal between γ 1 and γ 2 is independent of the location of the noise source ; The cross-correlation between the sound fields at the positions of the two receiving array elements γ 1 and γ 2 in the frequency domain is: <p(r1,ω)p*(r2,ω)>=|q(ω)|2∫G(r1,r,ω)G*(r2,r,ω)dV (3);<p(r 1 , ω)p * (r 2 , ω)>=|q(ω)| 2 ∫G(r 1 , r, ω)G * (r 2 , r, ω)dV (3); 公式(3)中,p(γ1,ω)代表接收阵元γ1处的声场,p(γ2,ω)代表接收阵元γ2处的声场,<p(γ1,ω)p*2,ω)>代表两个接收阵元γ1和γ2位置处声场之间的互相关函数,上角标*代表共轭复数;In formula (3), p(γ 1 , ω) represents the sound field at the receiving element γ 1 , p(γ 2 , ω) represents the sound field at the receiving element γ 2 , <p(γ 1 , ω)p *2 , ω)> represents the cross-correlation function between the sound fields at the positions of the two receiving array elements γ 1 and γ 2 , and the superscript * represents a complex conjugate number; 然后,由公式(1)和(3)得到:Then, from formulas (1) and (3), we get: (G(r1,r2,ω)-G*(r1,r2,ω))|q(ω)|2=-2jω<p(r1,ω)p*(r2,ω)> (4);(G(r 1 , r 2 , ω)-G * (r 1 , r 2 , ω))|q(ω)| 2 =-2jω<p(r 1 , ω)p * (r 2 , ω) >(4); 公式(4)中,等式左侧为格林函数G(γ1,γ2,ω)以及它对应的时间反转,即对应频域中的共轭操作,把这两项乘以随机噪声的功率谱密度后,其结果等于等式右侧扩散场中两个的接收阵元γ1和γ2互相关结果,上角标*代表共轭复数;In formula (4), the left side of the equation is the Green's function G(γ 1 , γ 2 , ω) and its corresponding time reversal, that is, the conjugate operation in the corresponding frequency domain. Multiply these two terms by the random noise. After the power spectral density, the result is equal to the cross-correlation result of the two receiving array elements γ 1 and γ 2 in the diffuse field on the right side of the equation, and the superscript * represents the conjugate complex number; 然后,求出公式(4)所对应的时域表达式,频域2jω对应时域2d/dt,根据卷积定理可知,频域中的乘积对应时域中的卷积,从而得到:Then, the time domain expression corresponding to formula (4) is obtained. The frequency domain 2jω corresponds to the time domain 2d/dt. According to the convolution theorem, the product in the frequency domain corresponds to the convolution in the time domain, thus obtaining:
Figure FDA0002444267940000021
Figure FDA0002444267940000021
公式(5)中,G(γ1,γ2,t)代表时域中γ1和γ2之间的格林函数,G(γ1,γ2,-t)代表G(γ1,γ2,t)的时间反转,即对应频域中的共轭操作,*代表卷积运算,
Figure FDA0002444267940000022
代表互相关运算,Cq(t)表示扩散场中的噪声q(t)的自相关结果,d/dt代表对t求导,p(γ1,t)代表时域中接收阵元γ1处的声场,p(γ2,t)代表时域中接收阵元γ2处的声场;
In formula (5), G(γ 1 , γ 2 , t) represents the Green's function between γ 1 and γ 2 in the time domain, and G(γ 1 , γ 2 , -t) represents G(γ 1 , γ 2 , t) time reversal, which corresponds to the conjugate operation in the frequency domain, * represents the convolution operation,
Figure FDA0002444267940000022
represents the cross-correlation operation, C q (t) represents the autocorrelation result of the noise q(t) in the diffuse field, d/dt represents the derivation of t, p(γ 1 , t) represents the receiving array element γ 1 in the time domain The sound field at , p(γ 2 , t) represents the sound field at the receiving array element γ 2 in the time domain;
公式(4)和(5)的结果表明扩散场中接收阵元声场p(γ1,t)和p(γ2,t)互相关并求导,其结果等于这两接收阵元之间的格林函数响应,根据声波的互易性原理,在时间轴上,这两接收阵元之间的格林函数响应具有对称性;The results of formulas (4) and (5) show that the sound fields p(γ 1 , t) and p(γ 2 , t) of the receiving array elements in the diffuse field are correlated and derived, and the result is equal to the difference between the two receiving array elements. Green's function response, according to the principle of reciprocity of sound waves, on the time axis, the Green's function response between the two receiving array elements has symmetry; 2)频域合成孔径聚焦(SAFT)成像:2) Frequency Domain Synthetic Aperture Focusing (SAFT) imaging: 首先,由于超声波在无限的均匀介质传播满足Helmholtz波动方程,因此,假设在二维空间x-z坐标轴下,超声波发射阵元的坐标为(u,0),接收阵元的坐标为(v,0),缺陷目标距离发射阵元和接收阵元的距离分别为ρin和ρout,任意一对发射接收阵元对应的时域接收信号为e(v,u,t),那么对应的发射接收阵元所对应的频率响应E(v,u,ω)为:First of all, since the propagation of ultrasonic waves in an infinite homogeneous medium satisfies the Helmholtz wave equation, it is assumed that the coordinates of the ultrasonic transmitting array element are (u, 0) and the coordinates of the receiving array element are (v, 0 under the xz coordinate axis of the two-dimensional space). ), the distances of the defect target from the transmitting array element and the receiving array element are ρ in and ρ out respectively, and the time domain receiving signal corresponding to any pair of transmitting and receiving array elements is e(v, u, t), then the corresponding transmitting and receiving The frequency response E(v, u, ω) corresponding to the array element is: E(v,u,ω)=P(ω)∫∫f(x,z)g(ρin,ω)g(ρout,ω)dxdz (6);E(v, u, ω)=P(ω)∫∫f(x, z)g(ρ in ,ω)g(ρ out ,ω)dxdz (6); 公式(6)中,
Figure FDA0002444267940000031
P(ω)是发射信号的频谱,ω为发射信号的角频率;f(x,z)代表缺陷目标的点扩散函数,g(ρin,ω)代表发射阵元到缺陷目标的空间格林函数频域响应,g(ρout,ω)代表接收阵元到缺陷目标的空间格林函数频域响应;
In formula (6),
Figure FDA0002444267940000031
P(ω) is the spectrum of the transmitted signal, ω is the angular frequency of the transmitted signal; f(x, z) represents the point spread function of the defective target, and g(ρ in , ω) represents the spatial Green's function from the transmitting array element to the defective target Frequency domain response, g(ρ out , ω) represents the spatial Green's function frequency domain response from the receiving array element to the defect target;
此外,格林函数为:Also, the Green's function is:
Figure FDA0002444267940000032
Figure FDA0002444267940000032
公式(7)中,g(x,z,ω)代表二维空间x-z坐标轴下任意点的格林函数,波数k=ω/c,c为声波在介质中的纵波速度,空间变量x和z与波数kx和kz相对应;In formula (7), g(x, z, ω) represents the Green’s function of any point under the xz coordinate axis of the two-dimensional space, the wave number k=ω/c, c is the longitudinal wave velocity of the sound wave in the medium, and the spatial variables x and z Corresponds to the wave numbers k x and k z ; 使用变量x和z代替掉变量ρin和ρout,即g(x,z,ω)代替g(ρin、ρout,ω),将公式(7)带入公式(6),重新获得接收响应E(v,u,ω),然后用波数ku和kv分别表示发射阵元(u,0)和接收阵元(v,0)的空间信息,为了实现将积分转换二维傅里叶变换,对接收响应E(v,u,ω)中变量u和v进行傅里叶变换得到:Use the variables x and z to replace the variables ρ in and ρ out , that is, g(x, z, ω) instead of g(ρ in , ρ out , ω), bring the formula (7) into the formula (6), and regain the reception Response E(v, u, ω), and then use the wave numbers ku and k v to represent the spatial information of the transmitting array element (u, 0) and the receiving array element (v, 0), respectively. In order to realize the transformation of the integral into a two-dimensional Fourier Leaf transform, the Fourier transform of the variables u and v in the received response E(v, u, ω) is obtained:
Figure FDA0002444267940000041
Figure FDA0002444267940000041
公式(8)中,E(ω,ku,kv)代表频域响应E(v,u,ω)的空间傅里叶变换,F代表点扩散函数f(x,z)的空间傅里叶变换;In formula (8), E(ω, k u , k v ) represents the spatial Fourier transform of the frequency domain response E(v, u, ω), and F represents the spatial Fourier transform of the point spread function f(x, z) leaf transform; 相应的,点扩散函数f(x,z)的空间傅里叶变换是F(kx,kz),为了求得图像域,反变换模型为:Correspondingly, the spatial Fourier transform of the point spread function f(x, z) is F(k x , k z ). In order to obtain the image domain, the inverse transform model is:
Figure FDA0002444267940000042
Figure FDA0002444267940000042
公式(9)中,
Figure FDA0002444267940000043
代表缺陷目标等间隔采样的数据域模型,S-1{}代表Stolt映射,完成了从数据域到图像域的转换,p*(ω)代表发射信号的共轭频谱;
In formula (9),
Figure FDA0002444267940000043
represents the data domain model of the defect target sampled at equal intervals, S -1 {} represents the Stolt map, which completes the conversion from the data domain to the image domain, p * (ω) represents the conjugate spectrum of the transmitted signal;
对公式(9)进行非线性坐标变化,将波数ku,kv和k映射到和图像域缺陷目标相关联的波数kx和kz,则:The nonlinear coordinate change of formula (9) is performed to map the wavenumbers ku, kv and k to the wavenumbers kx and kz associated with the defect target in the image domain, then: kx=ku+kv (10);k x =k u +k v (10);
Figure FDA0002444267940000044
Figure FDA0002444267940000044
然后通过二维傅里叶反变换求得缺陷目标分布的图像Then the image of the defect target distribution is obtained by the two-dimensional inverse Fourier transform
Figure FDA0002444267940000045
Figure FDA0002444267940000045
公式(12)中,
Figure FDA0002444267940000046
代表缺陷目标分布的图像域模型。
In formula (12),
Figure FDA0002444267940000046
An image domain model representing the distribution of defect targets.
3.一种实现轨道板裂缝自动检测的成像装置,包括轨检小车及设置在轨检小车上的摄像仪和计算机,所述摄像仪与计算机信号连接,其特征在于:还包括超声波发射接收器和超声线性阵列探头,所述超声线性阵列探头的发射端与超声波发射接收器的发射接口信号连接,所述超声线性阵列探头的接收端与超声波发射接收器的接收接口信号连接,所述超声波发射接收器的接收接口与计算机信号连接,所述超声线性阵列探头与一自适应伸缩机构相连接,所述自适应伸缩机构与一双轴位移调节机构相连接,且所述双轴位移调节机构固定连接在轨检小车的前端。3. an imaging device that realizes automatic detection of track plate cracks, comprising a track inspection trolley and a camera and a computer arranged on the track inspection trolley, the camera is connected with the computer signal, and it is characterized in that: also include an ultrasonic transmitter receiver and an ultrasonic linear array probe, the transmitting end of the ultrasonic linear array probe is signal-connected with the transmitting interface of the ultrasonic transmitting receiver, the receiving end of the ultrasonic linear array probe is signally connected with the receiving interface of the ultrasonic transmitting and receiving device, and the ultrasonic transmitting The receiving interface of the receiver is connected with the computer signal, the ultrasonic linear array probe is connected with an adaptive telescopic mechanism, the adaptive telescopic mechanism is connected with a biaxial displacement adjustment mechanism, and the biaxial displacement adjustment mechanism is fixedly connected At the front end of the rail inspection trolley. 4.根据权利要求3所述的成像装置,其特征在于:所述自适应伸缩机构包括垂向设置的、互相平行的安装板A和安装板B,所述安装板A和安装板B之间连接有可水平伸缩的伸缩支架,所述超声线性阵列探头与安装板A相连接,所述双轴位移调节机构与安装板B相连接。4 . The imaging device according to claim 3 , wherein the self-adaptive telescopic mechanism comprises a mounting plate A and a mounting plate B that are vertically arranged and parallel to each other, and the mounting plate A and the mounting plate B are located between the mounting plate A and the mounting plate B. 5 . A horizontally retractable telescopic bracket is connected, the ultrasonic linear array probe is connected with the mounting plate A, and the biaxial displacement adjusting mechanism is connected with the mounting plate B. 5.根据权利要求4所述的成像装置,其特征在于:所述安装板A远离安装板B的一侧设有探头连接件,所述超声线性阵列探头固定安装在探头连接件上。5 . The imaging device according to claim 4 , wherein a probe connector is provided on the side of the mounting plate A away from the mounting plate B, and the ultrasonic linear array probe is fixedly installed on the probe connector. 6 . 6.根据权利要求4所述的成像装置,其特征在于:在超声线性阵列探头的至少一侧设有激光测距仪。6 . The imaging device according to claim 4 , wherein a laser range finder is provided on at least one side of the ultrasonic linear array probe. 7 . 7.根据权利要求3所述的成像装置,其特征在于:所述双轴位移调节机构包括安装板C、纵向位移调节机构和横向位移调节机构,所述自适应伸缩机构与安装板C固定连接,所述安装板C与纵向位移调节机构滑动连接,所述纵向位移调节机构与横向位移调节机构滑动连接,所述横向位移调节机构固定连接在轨检小车的前端。7 . The imaging device according to claim 3 , wherein the biaxial displacement adjusting mechanism comprises a mounting plate C, a longitudinal displacement adjusting mechanism and a lateral displacement adjusting mechanism, and the adaptive telescopic mechanism is fixedly connected with the mounting plate C 7 . , the mounting plate C is slidably connected with the longitudinal displacement adjustment mechanism, the longitudinal displacement adjustment mechanism is slidably connected with the lateral displacement adjustment mechanism, and the lateral displacement adjustment mechanism is fixedly connected to the front end of the rail inspection trolley. 8.根据权利要求7所述的成像装置,其特征在于:所述纵向位移调节机构包括纵向支架和纵向电动丝杆调节机构,所述横向位移调节机构包括横向支架和横向电动丝杆调节机构,所述安装板C与纵向支架上下滑动连接,所述纵向支架与横向支架横向滑动连接,所述横向支架固定连接在轨检小车的前端。8 . The imaging device according to claim 7 , wherein the longitudinal displacement adjustment mechanism comprises a longitudinal bracket and a longitudinal electric screw adjustment mechanism, and the lateral displacement adjustment mechanism includes a transverse bracket and a transverse electric screw adjustment mechanism, 9 . The mounting plate C is slidably connected up and down with the longitudinal bracket, the longitudinal bracket is slidably connected with the transverse bracket, and the transverse bracket is fixedly connected to the front end of the rail inspection trolley. 9.根据权利要求3所述的成像装置,其特征在于:在轨检小车的前端还设有导向机构,构成导向机构的导向轮与钢轨滚动连接。9 . The imaging device according to claim 3 , wherein a guide mechanism is further provided at the front end of the on-rail inspection trolley, and the guide wheels constituting the guide mechanism are connected with the rail in a rolling manner. 10 . 10.根据权利要求3所述的成像装置,其特征在于:在轨检小车上还设有报警模块,所述报警模块与计算机信号连接。10 . The imaging device according to claim 3 , wherein an alarm module is further provided on the rail inspection trolley, and the alarm module is signal-connected to the computer. 11 .
CN202010274433.1A 2020-04-09 2020-04-09 An imaging method and device for realizing automatic detection of track slab cracks Active CN111452830B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010274433.1A CN111452830B (en) 2020-04-09 2020-04-09 An imaging method and device for realizing automatic detection of track slab cracks

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010274433.1A CN111452830B (en) 2020-04-09 2020-04-09 An imaging method and device for realizing automatic detection of track slab cracks

Publications (2)

Publication Number Publication Date
CN111452830A true CN111452830A (en) 2020-07-28
CN111452830B CN111452830B (en) 2021-04-09

Family

ID=71673653

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010274433.1A Active CN111452830B (en) 2020-04-09 2020-04-09 An imaging method and device for realizing automatic detection of track slab cracks

Country Status (1)

Country Link
CN (1) CN111452830B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111845844A (en) * 2020-07-30 2020-10-30 湖南第一师范学院 A mobile dual-track flaw detection and positioning device and its use method
CN112781654A (en) * 2020-12-31 2021-05-11 西南交通大学 Crack steel rail gap fault detection system
CN114152676A (en) * 2021-11-16 2022-03-08 上海工程技术大学 Method for realizing automatic detection of wind power blade defects based on ultrasonic waves
CN117825512A (en) * 2024-01-04 2024-04-05 浙江大学 High-speed ultrasonic dynamic imaging detection method and device for rail head defects

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003260052A (en) * 2002-03-08 2003-09-16 Hitachi Medical Corp Image diagnosis support device
WO2007011909A1 (en) * 2005-07-18 2007-01-25 Atlantis Components, Inc. Registration of 3-d imaging of 3-d objects
CN105004792A (en) * 2015-07-20 2015-10-28 北京工业大学 Nonlinear ultrasonic phased array imaging method for micro-crack detection
CN108657223A (en) * 2018-07-23 2018-10-16 中国安全生产科学研究院 A kind of urban track traffic automatic tour inspection system and tunnel deformation detecting method
CN108693253A (en) * 2018-05-02 2018-10-23 南昌航空大学 A kind of rapid phase-control battle array ultrasound total focus imaging technique
CN208541707U (en) * 2018-01-04 2019-02-26 吴申田 Multifunctional abdominal radiotherapy precise locator
CN109781849A (en) * 2019-01-25 2019-05-21 江苏大学 A damage imaging method based on frequency-domain sub-modal superposition
CN110057829A (en) * 2019-04-30 2019-07-26 南水北调工程设计管理中心 A kind of comprehensive detection device and detection method in inner wall crack and deformation
CN110758457A (en) * 2019-11-14 2020-02-07 义乌轩久铁路技术有限公司 Ultrasonic rail detection device
CN110823904A (en) * 2019-10-31 2020-02-21 王佩洁 Hydraulic engineering crack extraction method
CN211856387U (en) * 2020-04-09 2020-11-03 上海工程技术大学 A signal acquisition device for crack detection of ballastless track slab

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003260052A (en) * 2002-03-08 2003-09-16 Hitachi Medical Corp Image diagnosis support device
WO2007011909A1 (en) * 2005-07-18 2007-01-25 Atlantis Components, Inc. Registration of 3-d imaging of 3-d objects
CN105004792A (en) * 2015-07-20 2015-10-28 北京工业大学 Nonlinear ultrasonic phased array imaging method for micro-crack detection
CN208541707U (en) * 2018-01-04 2019-02-26 吴申田 Multifunctional abdominal radiotherapy precise locator
CN108693253A (en) * 2018-05-02 2018-10-23 南昌航空大学 A kind of rapid phase-control battle array ultrasound total focus imaging technique
CN108657223A (en) * 2018-07-23 2018-10-16 中国安全生产科学研究院 A kind of urban track traffic automatic tour inspection system and tunnel deformation detecting method
CN109781849A (en) * 2019-01-25 2019-05-21 江苏大学 A damage imaging method based on frequency-domain sub-modal superposition
CN110057829A (en) * 2019-04-30 2019-07-26 南水北调工程设计管理中心 A kind of comprehensive detection device and detection method in inner wall crack and deformation
CN110823904A (en) * 2019-10-31 2020-02-21 王佩洁 Hydraulic engineering crack extraction method
CN110758457A (en) * 2019-11-14 2020-02-07 义乌轩久铁路技术有限公司 Ultrasonic rail detection device
CN211856387U (en) * 2020-04-09 2020-11-03 上海工程技术大学 A signal acquisition device for crack detection of ballastless track slab

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111845844A (en) * 2020-07-30 2020-10-30 湖南第一师范学院 A mobile dual-track flaw detection and positioning device and its use method
CN111845844B (en) * 2020-07-30 2021-07-09 湖南第一师范学院 Movable double-rail flaw detection device and using method thereof
CN112781654A (en) * 2020-12-31 2021-05-11 西南交通大学 Crack steel rail gap fault detection system
CN112781654B (en) * 2020-12-31 2023-03-17 西南交通大学 Crack steel rail gap fault detection system
CN114152676A (en) * 2021-11-16 2022-03-08 上海工程技术大学 Method for realizing automatic detection of wind power blade defects based on ultrasonic waves
CN117825512A (en) * 2024-01-04 2024-04-05 浙江大学 High-speed ultrasonic dynamic imaging detection method and device for rail head defects
CN117825512B (en) * 2024-01-04 2024-06-11 浙江大学 High-speed ultrasonic dynamic imaging detection method and device for rail head defects

Also Published As

Publication number Publication date
CN111452830B (en) 2021-04-09

Similar Documents

Publication Publication Date Title
CN111452830A (en) An imaging method and device for realizing automatic detection of track slab cracks
CN108956761B (en) Steel plate full-coverage ultrasonic detection device and method
Pelivanov et al. High resolution imaging of impacted CFRP composites with a fiber-optic laser-ultrasound scanner
CN102967657B (en) Nondestructive testing device and method based on synthetic aperture ultrasonic imaging technique
Iyer et al. Evaluation of ultrasonic inspection and imaging systems for concrete pipes
CN109164173B (en) A method and device for multi-channel dynamic non-destructive detection of ballastless track defects
CN111323485A (en) Imaging method and device for detecting internal defects of track slab
CN211856387U (en) A signal acquisition device for crack detection of ballastless track slab
CN111337574B (en) Topological imaging method and device based on ultrasonic lamb wave plate type ballastless track cracks
CN111044569A (en) Tunnel concrete structure defect detection method
CN109765295A (en) Laser ultrasonic rapid detection method and device for micro-cracks on concrete surface
CN111307945B (en) An imaging method and device for detecting near-surface defects of ballastless track based on ultrasonic array
Li et al. Crack detection in monocrystalline silicon solar cells using air-coupled ultrasonic lamb waves
He et al. Quantitative detection of surface defect using laser-generated Rayleigh wave with broadband local wavenumber estimation
CN108344801A (en) II plate-type ballastless track mortar bond quality detection devices of CRTS and method
CN106546661B (en) A kind of interference ultrasonic synthetic aperture is to inside concrete imaging method
CN114942268A (en) Tunnel lining crack detection method and inspection system based on acoustic wave change matching
Sun et al. Using Lamb wave TDTE method to realize ultrasonic array super-resolution imaging of multiple asymmetric defects
CN119757366A (en) A multi-modal defect nondestructive detection device and detection method for metal rolling element
CN209103356U (en) A Wind Turbine Basic Quality Data Acquisition System Based on Ultrasonic CT
CN110057915B (en) Underwater full-section three-dimensional foundation imaging detection method for channel engineering
IEM-RM B-scan ultrasonic image analysis for internal rail defect detection
Zhu et al. Time‐Domain Topological Energy Imaging Method of Concrete Cavity Defect by Lamb Wave
CN207366506U (en) The model test apparatus of concrete for hydraulic structure hidden defects detection
CN115015913A (en) A millimeter-wave near-field high-precision three-dimensional imaging method for building curtain wall detection

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