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 PDFInfo
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 36
- 238000001514 detection method Methods 0.000 title claims abstract description 28
- 239000000523 sample Substances 0.000 claims abstract description 72
- 230000007547 defect Effects 0.000 claims abstract description 53
- 230000004044 response Effects 0.000 claims abstract description 29
- 239000000284 extract Substances 0.000 claims abstract description 5
- 230000007246 mechanism Effects 0.000 claims description 79
- 238000006073 displacement reaction Methods 0.000 claims description 43
- 238000007689 inspection Methods 0.000 claims description 33
- 230000001364 causal effect Effects 0.000 claims description 18
- 230000000875 corresponding effect Effects 0.000 claims description 18
- 238000001228 spectrum Methods 0.000 claims description 18
- 230000003044 adaptive effect Effects 0.000 claims description 12
- 238000004613 tight binding model Methods 0.000 claims description 12
- 230000003595 spectral effect Effects 0.000 claims description 9
- 230000002950 deficient Effects 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 6
- 238000009792 diffusion process Methods 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000002596 correlated effect Effects 0.000 claims description 3
- 238000005314 correlation function Methods 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 2
- 230000002269 spontaneous effect Effects 0.000 claims description 2
- 238000012423 maintenance Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000000034 method Methods 0.000 description 4
- 230000001066 destructive effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000011897 real-time detection Methods 0.000 description 3
- 239000004567 concrete Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000004570 mortar (masonry) Substances 0.000 description 2
- 244000025254 Cannabis sativa Species 0.000 description 1
- 235000012766 Cannabis sativa ssp. sativa var. sativa Nutrition 0.000 description 1
- 235000012765 Cannabis sativa ssp. sativa var. spontanea Nutrition 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 235000009120 camo Nutrition 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 235000005607 chanvre indien Nutrition 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000011487 hemp Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61K—AUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
- B61K9/00—Railway 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/08—Measuring installations for surveying permanent way
- B61K9/10—Measuring installations for surveying permanent way for detecting cracks in rails or welds thereof
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection 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软件对所接收到的信号被动提取格林函数,对扩散场信号进行互相关,重建阵元之间的格林函数,获取阵元之间未延时的响应,恢复被噪声淹没的早期缺陷信息,然后依据频域合成孔径聚焦成像算法进行轨道板裂缝成像。本发明可以高效、无损、实时的检测无砟轨道板裂缝表面及裂缝内部缺陷。
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.
Description
技术领域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(γ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 in the 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 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:
公式(5)中,G(γ1,γ2,t)代表时域中γ1和γ2之间的格林函数,G(γ1,γ2,-t)代表G(γ1,γ2,t)的时间反转,即对应频域中的共轭操作,*代表卷积运算,代表互相关运算,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, 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)中,P(ω)是发射信号的频谱,ω为发射信号的角频率;f(x,z)代表缺陷目标的点扩散函数,g(ρin,ω)代表发射阵元到缺陷目标的空间格林函数频域响应,g(ρout,ω)代表接收阵元到缺陷目标的空间格林函数频域响应;In formula (6), 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:
公式(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:
公式(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:
公式(9)中,代表缺陷目标等间隔采样的数据域模型,S-1{}代表Stolt映射,完成了从数据域到图像域的转换,p*(ω)代表发射信号的共轭频谱;In formula (9), 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):
然后通过二维傅里叶反变换求得缺陷目标分布的图像Then the image of the defect target distribution is obtained by the two-dimensional inverse Fourier transform
公式(12)中,代表缺陷目标分布的图像域模型。In formula (12), 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
本实施例中,所述超声线性阵列探头4采用市售产品即可,例如,可采用接触式的超声线性阵列探头,具有能量高,衰减小的优势,超声能量几乎不会泄露,更好地入射到待测材料的内部并与缺陷相互作用,形成对检测有利的散射信号,所述超声线性阵列探头4是由若干超声探头组成,例如,可以由8行3列共24个超声探头组成的阵列探头,相应的,所述超声线性阵列探头4即具有24个阵元,超声探头的中心频率可以为1.0MHz,每个探头的阵元可以单独进行激发,也可以组合激发,能够实现聚焦和偏转功能。此外,从上述可见,本实施例中的超声线性阵列探头4具有发射端和接收端,具备自发自收的功能。由于超声波在密度不同的介质中传播速度不同,当其通过两种不同介质的分界面时,会发生反射散射等现象,因无砟轨道属于混凝土构件,是由砂、水泥、石子等混合组成的多孔质非均匀复合材料,超声波发射接收器3发射的超声脉冲波在混凝土中传播时遇到缺陷会发生绕射,在缺陷界面会发生散射和反射,导致到达接收探头时声波能量得幅值显著减小,本申请中,采用超声线性阵列探头4,不仅能够更加有效率的采集扩散场信号,还可以使得其接收的扩散场信号叠加,使得反射扩散场信号在某一个方向的辐射能量最大,而在其它方向的总辐射能量较小,从而实现对扩散场信号的聚焦作用,使得有用的扩散场信号得到增强,进而使得干扰信号得到抑制,从而保证了信号采集及检测结果的精度。In this embodiment, the ultrasonic
参见图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
此外,所述安装板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
此外,在超声线性阵列探头4的至少一侧设有激光测距仪55,通过激光测距仪55所检测的位移,即可判断超声线性阵列探头4的发射接收端面与无砟轨道板7的表面之间的垂直距离,以保持垂直距离在预设的范围内,进而防止沉降距离过大而损伤探头。本实施例中,在超声线性阵列探头的左右两侧分别对称设有激光测距仪55,两个激光测距仪55可以设定为一用一备。In addition, a
参见图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
具体的,所述纵向位移调节机构62包括纵向支架621和纵向电动丝杆调节机构622,所述横向位移调节机构63包括横向支架631和横向电动丝杆调节机构632,所述安装板C61与纵向支架621上下滑动连接,所述纵向支架621与横向支架631横向滑动连接,所述横向支架631固定连接在轨检小车1的前端。本实施例中,所述纵向电动丝杆调节机构622和横向电动丝杆调节机构632均是由驱动电机、一端固定在驱动电机输出端的丝杆和与丝杆螺纹连接的滑块连接件组成,由于此组成结构为已知技术,故在图中未详细示出。Specifically, the longitudinal displacement adjustment mechanism 62 includes a
此外,参见图1和图2所示,在轨检小车1的前端还设有照明灯8,以避免在夜间或较暗环境下摄像仪2拍摄模糊无法识别裂缝的问题,避免漏检,增加了检测效率。照明灯8的数量根据需要灵活设定,例如,本实施例中,是在轨检小车1的前端左右对称设有两个照明灯8,In addition, as shown in FIG. 1 and FIG. 2 , the front end of the on-
此外,在轨检小车1的前端还设有导向机构9,构成导向机构9的导向轮91与钢轨10滚动连接,以对轨检小车1的行驶方向起到导向作用。In addition, a
此外,在轨检小车1上还设有报警模块11,所述报警模块11与计算机14信号连接,以便发现轨道上有裂缝存在时,可以通过计算机14发出指令以驱动报警模块11进行报警。所述报警模块11采用市售报警器即可。In addition, an
此外,在轨检小车1上还设有移动电源12,以实现对所述装置进行移动供电。In addition, a
此外,在轨检小车1上设置座椅13,以便由检测员坐在轨检小车1上进行人工信号采集及检测操作。In addition, a
采用本发明所述成像装置实现轨道板裂缝自动检测的成像方法如下: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
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
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(γ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:
公式(5)中,G(γ1,γ2,t)代表时域中γ1和γ2之间的格林函数,G(γ1,γ2,-t)代表G(γ1,γ2,t)的时间反转,即对应频域中的共轭操作,*代表卷积运算,代表互相关运算,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, 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
首先,由于超声波在无限的均匀介质传播满足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)中,P(ω)是发射信号的频谱,ω为发射信号的角频率;f(x,z)代表缺陷目标的点扩散函数,g(ρin,ω)代表发射阵元到缺陷目标的空间格林函数频域响应,g(ρout,ω)代表接收阵元到缺陷目标的空间格林函数频域响应;In formula (6), 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:
公式(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:
公式(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:
公式(9)中,代表缺陷目标等间隔采样的数据域模型,S-1{}代表Stolt映射,完成了从数据域到图像域的转换,p*(ω)代表发射信号的共轭频谱;In formula (9), 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);
然后通过二维傅里叶反变换求得缺陷目标分布的图像Then the image of the defect target distribution is obtained by the two-dimensional inverse Fourier transform
公式(12)中,代表缺陷目标分布的图像域模型。In formula (12), 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)
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)
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)
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 |
-
2020
- 2020-04-09 CN CN202010274433.1A patent/CN111452830B/en active Active
Patent Citations (11)
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)
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 |