CN110231006B - 空气耦合超声干涉法 - Google Patents
空气耦合超声干涉法 Download PDFInfo
- Publication number
- CN110231006B CN110231006B CN201910494935.2A CN201910494935A CN110231006B CN 110231006 B CN110231006 B CN 110231006B CN 201910494935 A CN201910494935 A CN 201910494935A CN 110231006 B CN110231006 B CN 110231006B
- Authority
- CN
- China
- Prior art keywords
- frequency
- ultrasonic
- air
- difference
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000008878 coupling Effects 0.000 title claims abstract description 28
- 238000010168 coupling process Methods 0.000 title claims abstract description 28
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000005259 measurement Methods 0.000 claims abstract description 37
- 230000008859 change Effects 0.000 claims abstract description 19
- 230000005284 excitation Effects 0.000 claims description 11
- 230000005428 wave function Effects 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 claims description 2
- 238000005305 interferometry Methods 0.000 claims 7
- 238000004458 analytical method Methods 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 9
- 238000002592 echocardiography Methods 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 7
- 239000000463 material Substances 0.000 description 5
- 230000001427 coherent effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000002310 reflectometry Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000001028 reflection method Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 239000011358 absorbing material Substances 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000004556 laser interferometry Methods 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B17/00—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations
- G01B17/06—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations for measuring contours or curvatures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/534—Details of non-pulse systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/06—Systems determining the position data of a target
- G01S15/08—Systems for measuring distance only
- G01S15/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S15/34—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/06—Systems determining the position data of a target
- G01S15/08—Systems for measuring distance only
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/06—Systems determining the position data of a target
- G01S15/08—Systems for measuring distance only
- G01S15/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
- G01S15/36—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated with phase comparison between the received signal and the contemporaneously transmitted signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/93—Sonar systems specially adapted for specific applications for anti-collision purposes
- G01S15/931—Sonar systems specially adapted for specific applications for anti-collision purposes of land vehicles
Abstract
本发明公开了空气耦合超声干涉法，将空气耦合超声换能器作为测头正对工件表面，超声波在工件表面发生来回反射，测得第一次工件表面反射回波达到空气耦合超声换能器的相位差；基于超声波频率及波长的变化，则将距离的测量转变成声波相位相对于声波角频率的变化率；所述声波角频率的变化为空气耦合超声换能器带宽fB内最高频率F2与最低频率F1之差乘以2π；所述声波相位的变化为F2频率下与F1频率下声程除以波长所得整倍数的差值乘以2π。该超声技术基于长时间的持续发射和大跨度的频率扫描，具有测量精度高和抗干扰能力强的优点，其在激光干涉中的对应技术比较难以实现。
Description
技术领域
本发明涉及空气耦合超声测距技术领域,特别是空气耦合超声干涉法。
背景技术
复杂结构高精度纳米加工以后，需要对表面的自由曲面进行高精度测量。但是静态的自由曲面的高精度测量方法和技术一直是一个难题。超精密加工后表面轮廓精度超精密测量难、测量成本高、表面表征理论和测量方法未有应用性突破。
目前对复杂结构纳米加工表面的自由曲面测量的方法中高精度的测量方法行业内普遍采用三轴坐标机技术平台。这个技术平台包括高精度的标尺系统如线纹尺、精密丝杆、感应同步器、光栅尺、磁尺及光波波长等；包括导轨系统如滑动导轨、滚动轴承导轨和气浮导轨等；包括驱动、平衡、转台等技术模块系统，另外再加上技术上与三轴坐标机其它部分相对独立的三维测头。总体来说，对于三坐标机系统，位置和尺寸测量精度达到微米量级的已经是高端产品，关键要看测头测量精确度的技术水平。测头分为接触式与非接触式。接触式的探针类测头是目前主要使用的测头类型，探针类的测量精度一般不会优于微米量级。
非接触式测头的种类主要分为光学、超声和电磁几类。其中光学类测头是目前国内外研究的热点。光学类测头里面又分为线结构光法、激光三角法、激光半焦斑成像法等。目前公开文献见到的国内最好的实验室指标是采用的半光斑成像式激光测头实现的，其重复瞄准不确定度最好的效果达到微米量级，测量灵敏度可达30mV/μm，激光跟踪瞄准被测曲面倾角可达25°。但是该方法受材料类别影响非常大，材料的颜色以及材料表面的吸光性质都会影响测量的灵敏度，所以对不同的材料需要不同的标准试样。有时候材料的标准试样是很难获取的。
激光干涉式测头主要最对微小的相对变化量敏感，对绝对量不敏感，绝对测量精度受制于控制台的三维坐标瞄准。而且激光干涉测头对被测工件表面的激光反射特性要求很高，反射率太高太低都不好，有些特殊的吸波材料还需要要更换不同波长的激光。所以激光干涉法目前没有受到大规模推广应用。
空气超声法测距主要用时差法。一般来说，频率越高，测距精度越高。但是1MHz量级的高频超声波在空气中的衰减很大，以至于无法用于常见距离的测量。另外空气超声时差法受外界干扰较大，所以实际应用中时差法的测距精度大打折扣。举例来说汽车行业广泛应用的空气耦合超声测距倒车雷达采用的频率一般在几十KHz量级，其测距精度通常在一个发射波长的量级。比如一个40KHz的超声波倒车雷达的测距精度大约在10mm左右。
发明内容
本发明的技术方案是：空气耦合超声干涉法，如果空气耦合超声测头正对着工件表面，超声波换能器发射出的超声波在工件表面就会发生来回发射。超声波在空气中传播的波长λ为λ＝c/f (1)
式中c是空气声速，f是超声波频率。
如果频率变化，空气超声波测头声程不变，但是相位会改变。第一次工件表面反射回波到达空耦超声换能器的相位差Φ为：
Φ＝2πL/λ (2)
式中L是超声波从换能器到工件表面再反射回换能器的声程，空气耦合超声换能器到工件的距离则为L/2。(1)式代入(2)式中有
Φ＝2πLf/c (3)
如果超声波频率发生微小变化，相位差也发生微小变化。(3)式对微小频率变化Δf展开得到微小相位差变化ΔΦ：
ΔΦ＝2πLΔf/c＝LΔω/c (4)
其中Δω为微小角频率变化。继续整理后有
L＝c·ΔΦ/Δω (5)
从公式(5)可知，我们把距离的测量转变成了相位相对角频率ω的变化率，这个变化率是一个与被测量距离线性相关的常数。令式(5)中的ΔΦ为一个周期2π，则式(5)变为
L＝c·Δf＝c·(f2-f1) (6)
式中f1和f2是相邻的两个相位周期对应的频率。
优选的是，公式(6)还可以以另一种方法推出。假设在某一频率f1下总声程为波长λ1的整数n倍，则
L＝nλ1＝nc/f1 (7)
如果连续增加频率，可以得到一个频率f2，此频率下总声程为波长λ2的整数(n+1)倍
L＝(n+1)λ2＝(n+1)c/f2 (8)
联合(7)和(8)式消去n求解L也可以得到公式(6)。
定义空气耦合超声换能器带宽fB内最低的一个频率F1，该频率下声程L为波长的整数倍N1；定义空气耦合超声换能器带宽fB内最高的一个频率F2，该频率下声程L也为波长的整数倍N2。我们后面将会说明N1和N2可以用声时差法或者干涉法准确地求出。我们后面也将说明，通常换能器带宽与这两个频率的差值相差极小，即
F2-F1≈fB (9)
为了表述方便，我们以后会用频带宽度fB来代替公式(9)中的频率差来说明该频率差的物理意义。把频率差和相位差代入到公式(5)中有，
L＝c·(N2-N1)/(F2-F1) (10)
因为c，N2,N1都为非常容易确定的常数，由误差分析可知式(10)的总声程相对误差为
δL/L＝δ(F2-F1)/|F2-F1|＝(|δF2|+|δF1|)/|F2-F1| (11)
或者空气耦合超声换能器测头到工件的距离相对测量误差为：
(δL/2)/(L/2)＝(|δF2|+|δF1|)/|F2-F1| (12)
总声程L的测量绝对误差为
δL＝L·(|δF2|+|δF1|)/|F2-F1| (13)
或者空气耦合超声换能器测头到工件的距离测量绝对误差为
δL/2＝L/2·(δF2|+|δF1|)/|F2-F1|≈L/2·(|δF2|+|δF1|)/fB (14)
公式(12)和(14)为我们提供了估算测量相对精度和绝对精度的计算公式。公式(14)说明，要想提高测量的绝对精度，需要增大带宽fB(或者F2-F1)，减小被测距离L/2,减小F1和F2的测量不确定度。对于空气耦合超声换能器，由于单探头在被测工件表面的反射率几乎为100％，信号非常强，所以在空气超声测距应用中可以适当放宽应用带宽的范围，比如可以选用-20dB带宽。通常商业应用级别的压电型空气耦合超声换能器在空气中的插入损耗-6dB带宽在20％～50％左右，-20dB带宽能够达到100％。因此，公式(14)中的带宽fB可以近似用换能器的中心频率fW代替。不失一般性，公式(14)可以整理为：
δL/2≈L/2·(|δF2|+|δF1|)/fw (15)
举例说明，如果我们选用一个非常常见的300KHz的空气耦合超声压电换能器，换能器到工件的距离选为45mm，频率步进选为10Hz或者说频率测量误差为10Hz，根据公式(15)可知距离测量误差为
δL/2≈45mm·(10+10)/(300k)＝3×10-3mm＝3μm (16)
这个测量精度已经远远高于脉冲反射方法。
我们再回头论证一下公式(9)的合理性。由以上所选取的实验参数可知相邻两个周期的频率差为(f2-f1)＝c/L＝340/(90×10-3)＝3800Hz＝3.8KHz (17)
公式(9)中fB与(F2-F1)的最大差值百分比为：
2×3.8KHz/300KHz＝2.5％ (18)
所以公式(9)是有根据近似成立的，在说明物理意义时我们可以用带宽fB来代替F2-F1的频率差。这进一步说明带宽越宽，干涉测量方法越精确。
我们采用一种特殊的超声波激励方式来实现上述干涉法测量。我们用正负变化的正弦波或者方波一个固定频率持续激发空气耦合超声换能器，激发的时间大于一次回波的时间，小于第二次回波的时间。第一次回波时间和第二次回波时间可以通过一个短周期的脉冲反射方法快速获得。然后我们变换激发频率按照上面的激励方法进行频率扫描检测。以扫查角频率为横坐标，以发生干涉的时域信号的绝对值积分为纵坐标作二维图。则该二维图为一个类似正弦波变化的曲线，可以用一个正弦波函数使用回归算法来拟合该曲线。频带越宽，模拟的结果就越准确。模拟出来的正弦波函数周期可以认为是权利要求2中在ΔΦ为2π时的Δω。发生干涉的时域信号选取范围为从第二次回波时间始至超声波持续激发时间加上第一次回波时间为至的时间区间以内。
本发明的优点是：
第一、该超声技术基于超过一个回波时间的长时间的持续发射和相对带宽内大跨度的频率扫描，具有测量精度高和抗干扰能力强的优点，其在激光干涉中的对应技术比较难以实现。
第二、该技术理论相对测距精确度正比于频率的测量误差值与空气耦合超声换能器的带宽比值，其短距离测距精度可以达到亚微米量级，长距离测距精度可以优于1％的波长。
第三、该测距技术短距离可以应用于自由曲面轮廓的高精度扫描，长距离可以应用于汽车自动驾驶时超声波雷达抗干扰测距。
附图说明
下面结合附图及实施例对本发明作进一步描述：
图1是空气耦合超声测头测距示意图；
图2是空气耦合超声测头实物图；
图3是399.8KHz、50Vpp、120周期连续激发产生相干相长的频率图；
图4是402KHz、50Vpp、120周期连续激发产生相干相消的频率图；
图5是399.8KHz、50Vpp、3周期连续激发用于时差法测量的频率图；
图6是550KHz换能器置于工件上方25mm时的干涉部分幅值绝对值积分平均值与频率的变化；
具体实施方式
实施例：
选取一个空气耦合超声换能器垂直对准被测工件自由曲面，如图1所示。选用聚焦式换能器可以缩小焦斑，提高横向分辨率；选用直径小的频率高的空气耦合超声平面换能器也能达到类似效果。空气耦合超声换能器的外径一般在10mm到50mm之间，如图二的实物照片所示，加上适当的夹具后可以直接作为测头安装在三轴坐标测量机上。聚焦换能器焦斑一般在1mm到5mm之间，焦距在5mm到50mm之间；平面换能器晶片直径一般在10mm～50mm，近场区一般在10mm～100mm之间。空耦超声换能器在高频段(1MHz以上)一般都选用平面型。
选用空气耦合超声自动扫描系统作为测试平台，分别选择一个中心频率为400KHz的压电空气耦合超声换能器垂直固定在离工件表面大约42mm的上方，以及一个中心频率为550KHz的压电空气耦合超声换能器垂直固定在离工件表面大约25mm的上方。
一个400KHz换能器自发自收的声压电信号随时间变化的曲线，如图3至5所示。其中最图5中3个周期激发的数据可以用作比较粗略的时差法测距，多次重复测量精度优于1mm。图3和图4是连续波激发，连续激发的时间超过了第一次回波的时间，小于第二次回波时间，于是超过第一次回波时间的连续激发部分与第一次回波发生干涉效应。有的频率相干相长，有的频率相干相消，并且在第二次回波时反应出来(如图中方框所标示处)。将干涉部分幅值取绝对值积分平均，就会得到随频率类似正弦波变化的曲线。作为实例，图6为一个更高频率550KHz换能器放置于工件上方约25mm时得到的干涉部分幅值绝对值积分平均值与频率的变化。实际应用中，该曲线可以进行一次带通滤波变成比较光滑和易于数据处理的正弦波曲线。
从图6可以用肉眼估算大约从530KHz到570kHz干涉相位变化了6个周期，将该数据代入到公式(10)有：
也就是说换能器与工件上被测点的距离为51mm/2＝25.5mm，这与实验设置时放置的大致位置25mm比较接近。精确度达到微米量级的计算可以从图四中提取更加精确的频率数据得到。比如对图6进行滤波后用一个正弦波函数进行回归拟合，得到的精确的拟合的正弦函数周期。该周期乘以2π即为公式(5)中ΔΦ为2π时的Δω，从而声程L可以通过公式(5)计算出。
当换能器降低频率到40KHz与100KHz之间工作在与工件距离1米到20米的范围时，作者在实际应用实验中能够快速得到的绝对误差通常优于一个波长的1％量级。
在频率扫描时，为了提高测距速度，除了带宽两端以外，带宽以内的频率可以按照奈奎斯特最低采样定律快速采样扫频以确定方程式(10)中的N1和N2。按照前述的持续发射得到相干信号之前，可以先按照图三中的下图以低周期脉冲发射得到粗略的距离LC，根据下式可得带宽内连续最大频率扫描步进FS为
为了更进一步提高测距速度，该带宽中间的频率步进扫描也可以省略。图四中的周期可以大致估算为c/LC，这样可以得到式(11)中N1和N2的初值。这样在带宽的两端各精细地扫描一个频率周期足以精确地确定式(10)中的F1和F2。
本发明所提出的空气耦合超声高精度测距干涉技术中，测量的相对精度等于电子设备对频率测量的误差与空气耦合超声具体应用频带的比值。实际工程应用中，对频率的测量精度非常高，加上空气耦合超声换能器在空气中的应用相对频带非常高，所以该技术可以给出极高的距离测量精度。其短距离测距精度可以达到亚微米量级，长距离测距精度可以优于1％的波长。其中，短距离可以应用于自由曲面轮廓的高精度扫描，长距离可以应用于汽车自动驾驶时超声波雷达抗干扰测距等。
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明的。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明的所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。
Claims (8)
1.空气耦合超声干涉法，将空气耦合超声换能器作为测头正对工件表面，其特征在于：超声波在工件表面发生来回反射，测得第一次工件表面反射回波达到空气耦合超声换能器的相位差；基于超声波频率及波长的变化，则将距离的测量转变成声波相位相对于声波角频率的变化率；所述声波角频率的变化为空气耦合超声换能器带宽fB乘以2π；所述声波相位的变化为带宽上边界频率时与带宽下边界频率时声程除以波长所得数值的差值乘以2π；
声波相位相对于声波角频率的变化率，即声波相位差相对于声波角频率差的比值，其与声程的公式关系为：
L＝λ·(ΔΦ/2π)/(Δf/f)＝c·ΔΦ/Δω，
其中ΔΦ为声波相位差，Δω为声波角频率差；
声波相位相对于声波角频率的变化率是一个与被测距离线性相关的常数，其也可以转化为相邻的两个相位周期对应的频率之差，其与声程的公式关系为：
L＝c·Δf＝c·(f2-f1)；
将声波角频率差和声波相位差用空气耦合超声换能器带宽fB内声程为波长整数倍的最高频率F2与声程为波长整数倍的最低频率F1之间的关系来替换，则声程的公式转化为：
L＝c·(N2-N1)/(F2-F1)。
2.根据权利要求1所述的空气耦合超声干涉法，其特征在于：超声波在一个固定频率检测时持续激发的时间大于一次回波的时间，小于第二次回波的时间。
3.根据权利要求2所述的空气耦合超声干涉法，其特征在于：超声波进行扫频检测，以扫查角频率为横坐标，以发生干涉的时域信号绝对值积分为纵坐标作二维图，对该二维图所呈现的正弦波变化曲线用正弦波函数回归算法来拟合；扫频跨度的大小与拟合数据的误差大小反相关；模拟出来的正弦波函数周期为相邻的两个相位周期对应的角频率之差2π(f2-f1)，且此时声波相位差ΔΦ为2π。
4.根据权利要求3所述的空气耦合超声干涉法，其特征在于：发生干涉的时域信号范围为：从第二次回波时间开始，至超声波持续激发时间加上第一次回波的时间为结束的时间区间。
5.根据权利要求4所述的空气耦合超声干涉法，其特征在于：根据误差分析法，声程的相对误差公式为：δL/L＝δ(F2-F1)/|F2-F1|＝(|δF2|+|δF1|)/|F2-F1|。
6.根据权利要求5所述的空气耦合超声干涉法，其特征在于：声程的绝对误差公式为：δL＝L·(|δF2|+|δF1|)/|F2-F1|。
7.根据权利要求6所述的空气耦合超声干涉法，其特征在于：第一次工件表面反射回波达到空气耦合超声换能器的相位差Φ＝2πL/λ；其中λ为波长；L是超声波从换能器到工件表面再反射回换能器的声程，声程等于两倍的被测距离。
8.根据权利要求7所述的空气耦合超声干涉法，其特征在于：超声波频率发生变化时，波长发生变化，相位也发生变化；相位变化的差值对频率变化的差值的关系为：
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910494935.2A CN110231006B (zh) | 2019-06-10 | 2019-06-10 | 空气耦合超声干涉法 |
| PCT/CN2019/119709 WO2020248516A1 (zh) | 2019-06-10 | 2019-11-20 | 空气耦合超声干涉法 |
| JP2021573800A JP7285029B2 (ja) | 2019-06-10 | 2019-11-20 | 空気連成超音波干渉法 |
| EP19932344.5A EP3982080A4 (en) | 2019-06-10 | 2019-11-20 | AIR COUPLED ULTRASONIC INTERFEROMETRIC PROCESS |
| US17/614,538 US11892541B2 (en) | 2019-06-10 | 2019-11-20 | Air-coupled ultrasonic interferometry method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910494935.2A CN110231006B (zh) | 2019-06-10 | 2019-06-10 | 空气耦合超声干涉法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110231006A CN110231006A (zh) | 2019-09-13 |
| CN110231006B true CN110231006B (zh) | 2020-07-17 |
Family
ID=67859360
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910494935.2A Active CN110231006B (zh) | 2019-06-10 | 2019-06-10 | 空气耦合超声干涉法 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US11892541B2 (zh) |
| EP (1) | EP3982080A4 (zh) |
| JP (1) | JP7285029B2 (zh) |
| CN (1) | CN110231006B (zh) |
| WO (1) | WO2020248516A1 (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110231006B (zh) * | 2019-06-10 | 2020-07-17 | 苏州博昇科技有限公司 | 空气耦合超声干涉法 |
| CN110470253A (zh) * | 2019-09-26 | 2019-11-19 | 杭州电力设备制造有限公司 | 基于超声波的厚度测量方法、装置、电子设备及存储介质 |
| CN112255155B (zh) * | 2020-10-28 | 2021-10-15 | 浙江大学 | 一种颗粒浓度和粒径二维分布的旋转测量系统和方法 |
| CN117741602A (zh) * | 2022-09-22 | 2024-03-22 | 华为终端有限公司 | 一种信号处理方法与电子设备 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2745107A1 (en) * | 2011-08-17 | 2014-06-25 | EMPA Eidgenössische Materialprüfungs- und Forschungsanstalt | Air coupled ultrasonic contactless method for non-destructive determination of defects in laminated structures |
| KR101851706B1 (ko) * | 2017-06-19 | 2018-04-25 | 주식회사 디에스비연구소 | 암반 홀 표면 검사시스템 |
| CN108226304A (zh) * | 2017-12-29 | 2018-06-29 | 大连交通大学 | 一种基于测量模型的超声相控阵线扫描灵敏度计算方法 |
Family Cites Families (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5871421A (ja) * | 1981-10-23 | 1983-04-28 | Mitsubishi Heavy Ind Ltd | 振動物体の変位計測装置 |
| JPS60218086A (ja) * | 1984-04-13 | 1985-10-31 | Nec Corp | 距離測定装置 |
| SU1357708A1 (ru) * | 1985-10-09 | 1987-12-07 | Киевский технологический институт легкой промышленности | Способ бесконтактного измерени толщины плоских изделий |
| JPS6453186A (en) * | 1987-08-24 | 1989-03-01 | Suzuki Motor Co | Ultrasonic range finder |
| SU1755047A1 (ru) * | 1990-05-29 | 1992-08-15 | Киевский технологический институт легкой промышленности | Способ определени рассто ни |
| JP3790873B2 (ja) * | 1997-05-01 | 2006-06-28 | 泰 石井 | 音響式表面積計 |
| KR100707103B1 (ko) | 2004-06-09 | 2007-04-13 | 학교법인 포항공과대학교 | 공기중 파라메트릭 어레이를 이용한 초음파 거리 측정장치 및 방법 |
| DE102006001964A1 (de) * | 2006-01-13 | 2007-07-19 | Robert Bosch Gmbh | Messgerät, insbesondere Entfernungsmessgerät |
| JP2012037294A (ja) * | 2010-08-05 | 2012-02-23 | Jtekt Corp | 超音波計測方法および超音波工作物径測定装置 |
| CN102865839B (zh) * | 2012-09-21 | 2015-07-01 | 华南理工大学 | 一种基于宽带调频及接收补偿的超声波测厚方法及装置 |
| CN103616691A (zh) | 2013-11-19 | 2014-03-05 | 郑州中智电子科技有限公司 | 一种超声波测距的装置 |
| CH708881B1 (de) * | 2013-11-20 | 2017-06-15 | Besi Switzerland Ag | Durchlaufofen für Substrate, die mit Bauteilen bestückt werden, und Die Bonder. |
| CN103676240B (zh) * | 2013-12-20 | 2017-08-01 | 合肥京东方光电科技有限公司 | 表面检测装置 |
| JP2015175700A (ja) | 2014-03-14 | 2015-10-05 | 沖電気工業株式会社 | 位置推定装置、位置推定方法およびプログラム |
| WO2016103839A1 (ja) * | 2014-12-22 | 2016-06-30 | オリンパス株式会社 | 超音波診断装置、超音波診断装置の作動方法および超音波診断装置の作動プログラム |
| CN105004792A (zh) * | 2015-07-20 | 2015-10-28 | 北京工业大学 | 一种用于微裂纹检测的非线性超声相控阵成像方法 |
| JP6815786B2 (ja) | 2016-08-03 | 2021-01-20 | 株式会社東京精密 | 超音波変位計測装置及び超音波変位計測方法 |
| GB201618378D0 (en) * | 2016-10-31 | 2016-12-14 | Heriot-Watt Univ | Microwave sensor |
| CN106473777A (zh) * | 2016-12-12 | 2017-03-08 | 广东技术师范学院 | 一种超声波诊断装置、工作方法以及工作程序 |
| CN108445517B (zh) * | 2018-03-20 | 2019-09-06 | 北京邮电大学 | 一种定位信号滤波方法、装置及设备 |
| CN109696680B (zh) * | 2018-12-27 | 2021-07-16 | 北京哈工科教机器人科技有限公司 | 基于相位检测的高精度超声波测距装置及方法 |
| CN110231006B (zh) * | 2019-06-10 | 2020-07-17 | 苏州博昇科技有限公司 | 空气耦合超声干涉法 |
-
2019
- 2019-06-10 CN CN201910494935.2A patent/CN110231006B/zh active Active
- 2019-11-20 WO PCT/CN2019/119709 patent/WO2020248516A1/zh unknown
- 2019-11-20 JP JP2021573800A patent/JP7285029B2/ja active Active
- 2019-11-20 EP EP19932344.5A patent/EP3982080A4/en active Pending
- 2019-11-20 US US17/614,538 patent/US11892541B2/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2745107A1 (en) * | 2011-08-17 | 2014-06-25 | EMPA Eidgenössische Materialprüfungs- und Forschungsanstalt | Air coupled ultrasonic contactless method for non-destructive determination of defects in laminated structures |
| KR101851706B1 (ko) * | 2017-06-19 | 2018-04-25 | 주식회사 디에스비연구소 | 암반 홀 표면 검사시스템 |
| CN108226304A (zh) * | 2017-12-29 | 2018-06-29 | 大连交通大学 | 一种基于测量模型的超声相控阵线扫描灵敏度计算方法 |
Non-Patent Citations (2)
| Title |
|---|
| 《The science case for simultaneous mm-wavelength receivers in radio astronomy》;Richard Dodson等;《New Astronomy Reviews》;20170919;全文 * |
| 《中距离电磁超声导波检测及应用》;赵军辉等;《无损检测》;20151130;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| JP7285029B2 (ja) | 2023-06-01 |
| JP2022538783A (ja) | 2022-09-06 |
| EP3982080A1 (en) | 2022-04-13 |
| US20220229177A1 (en) | 2022-07-21 |
| CN110231006A (zh) | 2019-09-13 |
| US11892541B2 (en) | 2024-02-06 |
| EP3982080A4 (en) | 2023-07-12 |
| WO2020248516A1 (zh) | 2020-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110231006B (zh) | 空气耦合超声干涉法 | |
| CN107688051B (zh) | 一种基于激光超声表面波的亚表面缺陷宽度的测量方法 | |
| US9091664B2 (en) | Pulsed eddy current sensor for precision measurement at-large lift-offs on metallic surfaces | |
| JP2008545123A (ja) | 超音波減衰量を使用して材料特性を決定する方法及びシステム | |
| JP6392331B2 (ja) | 誘電性物質の厚さまたは深さの非破壊的絶対測定 | |
| Hunter et al. | Calibration and validation of ultrasonic reflection methods for thin-film measurement in tribology | |
| Zhang et al. | Calibration of the ultrasonic lubricant-film thickness measurement technique | |
| Haitjema | Uncertainty analysis of roughness standard calibration using stylus instruments | |
| Murugarajan et al. | Measurement, modeling and evaluation of surface parameter using capacitive-sensor-based measurement system | |
| Kuschmierz et al. | Optical, in situ, three-dimensional, absolute shape measurements in CNC metal working lathes | |
| Fei et al. | Simultaneous velocity, thickness and profile imaging by ultrasonic scan | |
| CN102012401A (zh) | 固体材料非均匀性质的无损检测方法 | |
| Cinbis et al. | A noncontacting technique for measuring surface tension of liquids | |
| Fomitchov et al. | Extrinsic and intrinsic fiberoptic Sagnac ultrasound sensors | |
| JP5256450B2 (ja) | 表面粗さ評価方法 | |
| Barakat et al. | Evaluation of new designed reference blocks for calibration and NDT by optical and ultrasonic techniques | |
| Lonardo et al. | Measurement and topography characterisation of surfaces produced by selective laser sintering | |
| Yamauchi et al. | A study of the stability, reliability, and accuracy of Neubrescope-based pipe thinning detection system | |
| KR20190128565A (ko) | 선형 단면 형상을 갖는 대상물의 표면 측정장치 및 측정방법 | |
| Etaix et al. | Thickness or phase velocity measurements using the Green's function comparison method | |
| Kazys et al. | Online profiling of nonplanar objects by high-resolution air-coupled ultrasonic distance measurements | |
| Suparta | Estimation of solid material surface roughness using time-of-flight ultrasound immerse transducer | |
| Brasunas et al. | Thickness measurement | |
| Sasaki et al. | Precise displacement measurements using phase information of 40 kHz ultrasonic waves in pinhole-based air-coupled ultrasonic system | |
| Pant et al. | Paint thickness measurement using acoustic interference |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 215000 room 906, building B, Wusongjiang business district, 169 Changhong North Road, Luzhi Town, Wuzhong District, Suzhou City, Jiangsu Province Applicant after: SUZHOU PHASERISE TECHNOLOGY Co.,Ltd. Address before: 215123 No. 188 Chuangyuan Road, Dushu Lake Science and Education Innovation Zone, Suzhou Industrial Park, Jiangsu Province Applicant before: SUZHOU PHASERISE TECHNOLOGY Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |