CN102721824B - 针对具有低速且高加速度特征的质点速度测量方法及装置 - Google Patents
针对具有低速且高加速度特征的质点速度测量方法及装置 Download PDFInfo
- Publication number
- CN102721824B CN102721824B CN201210216938.8A CN201210216938A CN102721824B CN 102721824 B CN102721824 B CN 102721824B CN 201210216938 A CN201210216938 A CN 201210216938A CN 102721824 B CN102721824 B CN 102721824B
- Authority
- CN
- China
- Prior art keywords
- optical fiber
- velocity
- beat interference
- laser
- speed
- 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
- 239000002245 particle Substances 0.000 title claims abstract description 74
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000001133 acceleration Effects 0.000 title claims abstract description 26
- 239000013307 optical fiber Substances 0.000 claims abstract description 67
- 238000012545 processing Methods 0.000 claims abstract description 26
- 238000012360 testing method Methods 0.000 claims abstract description 23
- 238000005259 measurement Methods 0.000 claims description 47
- 230000009466 transformation Effects 0.000 claims description 27
- 238000004891 communication Methods 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 21
- 230000003287 optical effect Effects 0.000 claims description 18
- 239000000835 fiber Substances 0.000 claims description 13
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 238000011430 maximum method Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 6
- 230000035939 shock Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 238000003672 processing method Methods 0.000 description 4
- 239000013077 target material Substances 0.000 description 4
- 229910001250 2024 aluminium alloy Inorganic materials 0.000 description 3
- 238000004880 explosion Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000000630 rising effect Effects 0.000 description 3
- 238000005305 interferometry Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005693 optoelectronics Effects 0.000 description 2
- 239000002243 precursor Substances 0.000 description 2
- 238000005057 refrigeration Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000000827 velocimetry Methods 0.000 description 2
- 229910000530 Gallium indium arsenide Inorganic materials 0.000 description 1
- 238000001069 Raman spectroscopy Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000035559 beat frequency Effects 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000001499 laser induced fluorescence spectroscopy Methods 0.000 description 1
- 238000004556 laser interferometry Methods 0.000 description 1
- 238000001307 laser spectroscopy Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004441 surface measurement Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Landscapes
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开一种针对具有低速且高加速度特征的质点速度测量方法及装置,对于速度为100~101m/s量级,并且加速度在109m/s2以上量级的质点表速度,本发明采用光纤激光多普勒干涉测速硬件***与采用特定流程的短时傅里叶变换或连续小波变换时频变换数据处理相结合的方法,实现具有上述特征运动靶面质点速度的高时间分辨率和速度分辨率测试,其中时间分辨率为0.2ns,速度精度为误差不超过0.1%。而且,本发明可以同时测试多路质点速度,测试效率高,成本较低,易于实现。
Description
技术领域
本发明涉及质点速度测量领域,尤其涉及一种针对具有低速且高加速度特征的质点速度测量方法。
背景技术
***与冲击产生的强动载荷具有高强度、短历时和小尺度的特点,对测量技术和测试装备提出了更高的要求。最近,高速摄影技术、X闪光照相技术、高速数字示波器以及相关的传感器、光源和快电子学设备已经开始应用于***与冲击实验。诊断技术发展的主流是以激光和光电子学技术为基础的快响应、高分辨率、高灵敏度的非接触测量技术,例如各种激光干涉测量和光谱技术(速度和位移干涉仪、VISAR、FPI等)、散斑技术、激光诱导荧光和拉曼光谱技术、光学窗口技术等。
自由表面速度测量是***与冲击实验研究的另一重点与难点。VISAR(Velocity Interferometer System for Any Reflector,任意反射面测速)***是目前应用最广泛的一种测量工具。随着科学技术的进步,VISAR***也不断改进与更新,现代光纤技术和相关光电子器件已开始引入到该测试技术中,在干涉测速原理上也实现了较大的突破,大大降低了干涉***对照明光源相干性和光强的要求。
但是,在冲击波物理中,往往要面对具有低速且高加速度特征的质点速度测量问题。对于速度为100~101m/s量级,并且加速度在109m/s2以上量级的自由表面速度,VISAR等测量手段却难以满足精度要求。
发明内容
本发明针对现有技术中无法精确测量具有低速且高加速度特征的质点速度的问题,提出一种针对具有低速且高加速度特征的质点速度测量方法及装置,依据信号特征有针对性地采取特定步骤的数据处理方法,能够较准确地测量具有低速且高加速度特征的质点速度。
为了解决上述问题,本发明提供一种针对具有低速且高加速度特征的质点速度测量方法,包括:
步骤1,光纤激光多普勒干涉测速***对速度为100~101m/s量级,并且加速度在109m/s2以上量级的运动靶面质点速度进行测量,输出差拍干涉信号至速度计算***;
步骤2,所述速度计算***采用短时傅里叶变换或连续小波变换对所述差拍干涉信号进行时频变换数据处理,以实现具有上述特征的运动靶面质点速度的高时间分辨率和速度分辨率测试。
优选地,上述方法具有以下特点:
所述光纤激光多普勒干涉测速***包括:激光器、光纤环形器、聚焦透镜、探测器和示波器,所述步骤1包括:
激光器发射通讯激光至光纤环形器的第一端口;
所述光纤环形器将接收到的通讯激光从其第二端口发送至聚焦透镜;
所述聚焦透镜将接收到的通讯激光的一部分作为参考光从其端面反射回所述光纤环形器的第二端口,另一部分输出至运动靶面;并将从运动靶面反射回的信号光发送至光纤环形器的第二端口;
所述光纤环形器将接收到的参考光和信号光从其第三端口发送至探测器;
所述探测器探测到所述参考光和信号光产生的差拍干涉信号,并发送至示波器;
所述示波器记录所述差拍干涉信号,并将所记录的差拍干涉信号发送至速度计算***。
优选地,上述方法具有以下特点:
在所述步骤2中,所述速度计算***针对运动靶面的质点速度特征以及差拍干涉信号的特征选择短时傅里叶变换或连续小波变换进行时频变换数据处理;其中,对于质点速度波动剧烈,信噪比大于5dB的差拍干涉信号,选择连续小波变换进行数据处理;对于质点速度波动较为平稳,信噪比小于5dB的干涉信号,选择短时傅里叶变换进行数据处理。
优选地,上述方法具有以下特点:
所述质点速度波动剧烈是指:质点速度波动的每个周期中,vmin≤0.2vmax;
所述质点速度波动较为平稳是指:质点速度波动的每个周期中,vmin≥0.8vmax;
其中,vmin为波谷对应的质点速度,vmax为波峰对应的质点速度。
优选地,上述方法具有以下特点:
采用连续小波变换对所述差拍干涉信号进行时频变换数据处理包括:
对所述差拍干涉信号进行滤波处理;
分别对预处理后的差拍干涉信号和原始的差拍干涉信号进行小波变换,得到它们各自的小波系数,再把对应的小波系数相乘;
用局部模极大值法进行小波脊线的提取,得到质点速度曲线。
优选地,上述方法具有以下特点:
在对所述差拍干涉信号进行滤波处理的步骤之前,所述方法还包括:
对所述差拍干涉信号进行短时傅里叶变换处理,并根据得到的结果估算质点速度的上限和下限,并计算小波基的中心频率及带宽,然后根据短时傅里叶变换的时频图,选择最优的小波基;
在对所述差拍干涉信号进行滤波处理的步骤中,根据估算的质点信号速度的上限和下限,设定滤波速度的上限和下限;
在分别对预处理后的差拍干涉信号和原始的差拍干涉信号进行小波变换的步骤中,采用所述最优的小波基进行小波变换。
优选地,上述方法具有以下特点:
在用局部模极大值法进行小波脊线的提取,得到质点速度曲线的步骤中,找到时频谱的模极大值,作为可靠的小波脊线的时间起点和差频频率起点,再根据加速度的变化范围,设置差频频率的搜寻范围,求出局部极大值,找到时间范围内各时刻的差频频率,从而实现小波脊线的提取。
为了解决上述问题,本发明提供一种针对具有低速且高加速度特征的质点速度测量装置,包括光纤激光多普勒干涉测速***和速度计算***,
所述光纤激光多普勒干涉测速***用于对速度为100~101m/s量级,并且加速度在109m/s2以上量级的运动靶面质点速度进行测量,输出差拍干涉信号至速度计算***;
所述速度计算***用于采用短时傅里叶变换或连续小波变换对所述差拍干涉信号进行时频变换数据处理,以实现具有上述特征的运动靶面质点速度的高时间分辨率和速度分辨率测试。
优选地,上述装置具有以下特点:
所述光纤激光多普勒干涉测速***包括:激光器、光纤环形器、聚焦透镜、探测器和示波器,其中,激光器、光纤环形器、探测器和示波器依次相连,聚焦透镜与光纤环形器相连;
所述激光器用于发射通讯激光至光纤环形器的第一端口;
所述光纤环形器用于将接收到的通讯激光从其第二端口发送至聚焦透镜,以及,将接收到的参考光和信号光从其第三端口发送至探测器;
所述聚焦透镜用于将接收到的通讯激光的一部分作为参考光从其端面反射回所述光纤环形器的第二端口,另一部分输出至运动靶面;并将从运动靶面反射回的信号光发送至光纤环形器的第二端口;
所述探测器用于探测所述参考光和信号光产生的差拍干涉信号,并发送至示波器;
所述示波器用于记录所述差拍干涉信号,并将所记录的差拍干涉信号发送至速度计算***。
优选地,上述装置具有以下特点:
所述光纤激光多普勒干涉测速***还包括光电转换模块、显示模块、第一分光模块和第二分光模块;所述聚焦透镜、光纤环形器、探测器和第二分光模块一一对应,均为N个,N大于等于1;
所述第一分光模块的入端与所述激光器相连,所述第一分光模块的N个出端分别与N个光纤环形器的第一端口相连,将所述激光器发射的通讯激光均分给N个光纤环形器;
每个第二分光模块的入端分别与对应的光纤环形器的第三端口相连,所述第二分光模块的一个出端与对应的探测器相连,另一个出端与光电转换模块相连;
所述光电转换模块用于将接收到的光信号转换成电信号,并发送给显示模块进行显示。
对于速度为100~101m/s量级,并且加速度在109m/s2以上量级的质点速度,本发明实现了具有上述特征运动靶面质点速度的高时间分辨率和速度分辨率测试,其中时间分辨率为0.2ns,速度精度为误差不超过0.1%。而且,本发明可以同时测试多路质点速度,测试效率高,成本较低,易于实现。
附图说明
图1为本发明实施例的针对具有低速且高加速度特征的质点速度测量装置示意图;
图2为本发明实施例的光纤激光多普勒干涉测速***示意图;
图3为本发明实施例的多路测试的光纤激光多普勒干涉测速***示意图;
图4为本发明应用实例的多路测试的光纤激光多普勒干涉测速***示意图;
图5为本发明应用实例的主机箱的示意图;
图6为本发明应用实例的LSP自由表面速度测量示意图;
图7为本发明应用实例的三组LSP自由表面速度测量结果;
图8为本发明应用实例的SHPB自由表面速度测量示意图,其中(a)为测量入射杆端面纵向速度;(b)为测量试样径向速度;
图9为本发明应用实例的SHPB测量不同处理方法比较,其中(a)为测量入射杆端面纵向速度得到的测量结果;(b)测量试样径向速度测得的测量结果。
具体实施方式
下文中将结合附图对本发明的实施例进行详细说明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互任意组合。
对于速度为100~101m/s量级,并且加速度在109m/s2以上量级的质点表速度,本发明采用光纤激光多普勒干涉测速硬件***与短时傅里叶变换或连续小波变换等时频变换数据处理软件相结合的方法,实现具有上述特征运动靶面质点速度的高时间分辨率和速度分辨率测试,其中时间分辨率为0.2ns,速度精度为误差不超过0.1%。
如图1所示,本发明的针对具有低速且高加速度特征的质点速度测量装置,包括光纤激光多普勒干涉测速***和速度计算***,其中,
所述光纤激光多普勒干涉测速***对速度为100~101m/s量级,并且加速度在109m/s2以上量级的运动靶面质点速度进行测量,输出差拍干涉信号至速度计算***;
所述速度计算***采用短时傅里叶变换或连续小波变换对所述差拍干涉信号进行时频变换数据处理,以实现具有上述特征的运动靶面质点速度的高时间分辨率和速度分辨率测试。
如图2所示,本发明实施例的光纤激光多普勒干涉测速***包括:激光器、光纤环形器、聚焦透镜、探测器和示波器,其中,激光器、光纤环形器、探测器和示波器依次相连,聚焦透镜与光纤环形器相连;其中,
激光器发射通讯激光(其初始频率为f0)至光纤环形器的第一端口;
所述光纤环形器将接收到的通讯激光从其第二端口发送至聚焦透镜;
所述聚焦透镜将接收到的通讯激光的一部分作为参考光(其频率仍为f0)从其端面反射回所述光纤环形器的第二端口,另一部分输出至运动靶面;并将从运动靶面反射回的信号光(因多普勒频移效应而具有频率f0)发送至光纤环形器的第二端口;
所述光纤环形器将接收到的参考光和信号光从其第三端口发送至探测器;
所述探测器探测到所述参考光和信号光产生的差拍干涉信号,并发送至示波器;
所述示波器记录所述差拍干涉信号,并将所记录的差拍干涉信号发送至速度计算***。
参考光和信号光产生差拍干涉信号,被探测器探测和示波器记录,其瞬时光强为
其中,I0(t)和Ib(t)分别为参考光和信号光光强幅值,u(t)为运动靶面瞬时质点速度,λ为激光波长,为初始相位差。运动靶面的质点速度为
Δfb(t)为拍频,由瞬时干涉条纹的疏密决定,可取相邻两个波峰或波谷时间间隔的倒数。可见,速度历程隐藏在干涉信号的频域内,而与光强变化无关,因此***的抗干扰能力较强。主要的***误差来源于激光器波长的不稳定性和基频噪声。
在一优选实例中,该光纤激光多普勒干涉测速***采用美国JDSU公司CQF938/400型窄线宽通讯激光器,输出光的波长为1550nm,线宽为200kHz,最大输出功率100mW。探测器采用New Focus公司1544-A型InGaAs PIN产品,带宽为12GHz;示波器采用Lecoy公司WaveMaster 808Zi型产品,带宽为8GHz,单通道最大写入速度40Gs/s。
如图3所示,可以增加分光器件,对图2实施例进行扩展,实现多路(N路)测量。其中,增加一个第一分光模块,聚焦透镜、光纤环形器和探测器一一对应,均为N个,N大于等于1。
所述第一分光模块的入端与所述激光器相连,所述第一分光模块的N个出端分别与N个光纤环形器的第一端口相连,将所述激光器发射的通讯激光均分给N个光纤环形器。N个探测器均接到同一个示波器上进行记录。
另外,可以增加N个第二分光模块,实现对N路光路的监测显示,如图3所示,每个第二分光模块的入端分别与对应的光纤环形器的第三端口相连,所述第二分光模块的一个出端与对应的探测器相连,另一个出端与光电转换模块相连;光电转换模块将接收到的光信号转换成电信号,并发送给显示模块进行显示。
如图4~图5所示,为本发明的一优选实例,采用两路光路测量,分光器件采用耦合器。
如图4所示,实验对象可放在第一微调架或第二微调架上,聚焦透镜采用自聚焦棒(包括第一自聚焦棒和第二自聚焦棒),第一探测器和第二探测器均连接示波器。测量时,采用实验测量通过激光脉宽上升沿同步触发。
如图5所示,在图4中的主机箱内部,包含有:激光源驱动电路板(含激光器)、为激光器制冷的制冷电路板、显示电路板、耦合器(第一耦合器、第二耦合器和第三耦合器)、环形器(第一环形器和第二环形器)、光电转换电路板和触发处理电路板。
其中,显示电路板可用于四个显示——驱动电流、制冷电流、两个探测点的光电流。
光电转换电路板将输入的光信号转换为电信号,并将光电流放大。
触发处理电路板用于触发信号的处理。
第一耦合器的分光比为50%:50%,第二耦合器和第三耦合器的分光比为90%:10%。
速度计算***针对运动靶面的质点速度特征以及差拍干涉信号的特征选择短时傅里叶变换或连续小波变换进行时频变换数据处理;其中,对于质点速度波动剧烈,信噪比大于5dB的差拍干涉信号,选择连续小波变换进行数据处理;对于质点速度波动较为平稳,信噪比小于5dB的干涉信号,选择短时傅里叶变换进行数据处理。
所述质点速度波动剧烈是指:质点速度波动的每个周期中,vmin≤0.2vmax;
所述质点速度波动较为平稳是指:质点速度波动的每个周期中,vmin≥0.8vmax;
其中,vmin为波谷对应的质点速度,vmax为波峰对应的质点速度。
(1)基于短时傅里叶进行数据处理时:
通过选取一个固定宽度的滑动窗口,将非平稳信号逐段截取为近似平稳的信号,然后分别对截得的信号进行傅里叶变换,得到每段中的频谱,进而得到运动靶面的质点速度。
其基本定义式为
STFT(τ,f)=∫x(t)g(t-τ)exp(-i2πft)dt(3)
x(t)表示信号,g(t-τ)表示中心在τ处的窗函数。该STFT(τ,f)即指上文Δfb(t)拍频。在处理时,取窗口函数为Hamming窗,窗函数受到Heisenberg不确定准则的限制,其时间分辨率Δt和频率分辨率Δf满足
Δf·Δt=1(4)
可见,STFT窗函数的时间与频率分辨率不能同时达到最优。
(2)基于连续小波变换进行数据处理时:
连续小波变换具有时频窗口可调的性质,在低频时使用较宽的窗口,在高频时使用较窄的窗口,因而在时域和频域上都有较好的分辨率,能更好地分析非平稳信号的时频特性。定义小波变换
其中,x(t)表示干涉信号,a∈R(实数)且a≠0,a为尺度因子,表示与频率相关的伸缩,b为时间平移因子。在(a,b)空间内得到max|Wψ(a,b)|的路径称为脊线,由此可得瞬时频响。时频变换常用的小波基为修正的Morlet小波或者修正的高斯小波,其时频窗口满足
Δf·Δt=1/4π(6)
通过与式(4)对比可见,CWT除了具有可变的窗口外,其时频分辨率也高于STFT。
采用连续小波变换对差拍干涉信号进行时频变换数据处理包括如下步骤:
(1)对所述差拍干涉信号进行滤波处理;
(2)分别对预处理后的差拍干涉信号和原始的差拍干涉信号进行小波变换,得到它们各自的小波系数,再把对应的小波系数相乘;
(3)用局部模极大值法进行小波脊线的提取,得到质点速度曲线。
其中,步骤(1)之前,还可对所述差拍干涉信号进行短时傅里叶变换处理,并根据得到的结果估算质点速度的上限和下限,并计算小波基的中心频率及带宽,然后根据短时傅里叶变换的时频图,选择最优的小波基;
在步骤(1)中,根据估算的质点信号速度的上限和下限,设定滤波速度的上限和下限;
在步骤(2)中,采用所述最优的小波基进行小波变换。
在步骤(3)中,找到时频谱的模极大值,作为可靠的小波脊线的时间起点和差频频率起点,再根据加速度的变化范围,设置差频频率的搜寻范围,求出局部极大值,找到时间范围内各时刻的差频频率,从而实现小波脊线的提取。
我们设计了两类试验,即激光冲击强化(LSP)和霍普金森杆(SHPB)试验,分别对LSP自由表面速度和SHPB端面及试样径向速度进行了测量。在LSP研究中,通常激光脉宽为几个纳秒,压力持续作用时间为十几个纳秒,弹性前驱波到达毫米量级厚度的靶材背表面所需的时间约为一两百个纳秒,在十几个纳秒内自由表面速度就由几米每秒上升至数百米每秒。因此,LSP自由表面测速对测试设备的灵敏度、时间分辨率和频率响应范围都提出了非常高的要求。Berther等人曾经尝试采用VISAR进行测量,但未能观察到数值计算能够模拟出的弹性前驱波。Arrigoni等人采用Fabry-Perot干涉仪初步观测到了LSP背表面弹性前驱波,但Fabry-Perot干涉仪价格过于昂贵。在SHPB试验中,端面速度和试样径向速度分别仅有十几米每秒和几米每秒,特别是在对后者的测量中出现较大的背景噪声。根据上述分析,需要通过有效的数据处理方法来提高信号的分辨率,以获取上升前沿速度变化历程的细节。利用我们自行研制的光纤激光多普勒干涉测速***,以及基于短时傅里叶(STFT)和连续小波变换(CWT)数据处理方法,能够对上述两类试验实现较好的测量。
实验一,LSP自由表面速度测量:
图6为LSP自由表面测量实验原理图。短脉冲激光经聚焦透镜汇聚到厚度为1.0mm的2024铝合金靶材上,光斑直径约为2mm。在强激光作用面依次紧密覆有0.1mm的铝膜吸收层和4.0mm厚的BK7玻璃约束层,共同装卡在冲击环境箱内。聚焦透镜置于靶体背表面并指向冲击中心位置,进行自由表面速度的测量。高功率密度短脉冲激光辐照在靶材上,产生等离子体并形成冲击波向靶材内部传播,到达靶材背表面时,引起质点运动。冲击波在靶材两个表面来回多次反射,其强度不断衰减,形成自由表面质点速度历程。实验测量通过激光脉宽上升沿同步触发。
图7为采用CWT处理的三个典型LSP试验结果,具有较好的一致性和规律性。由于应力波在靶材的两个表面来回多次反射,自由表面速度出现周期性的峰值。根据靶材厚度及弹塑性波达到背表面时间,可估算2024铝合金弹性波速约为6.11×103m/s,塑性波速约为5.14×103m/s,与文献中弹性波速Ce=6.20×103m/s,塑性波速D=5.52×103m/s相近。
根据Rankin-Hugoniot关系,弹性前驱波引起的自由表面速度对应着Hugoniot弹性极限
其中ρ0、Ce分别为靶材的初始密度和一维应变弹性波速。三个试验测量中最大值为59.7m/s,可得2024铝合金的HEL约为513MPa,与文献中的522MPa接近。在Hugoniot阶段采用的频率分辨率为Δf≈4.9MHz,对应的速度分辨率为3.9m/s。
实验二,SHPB端面及试样径向速度测量:
如图8所示,分别对SHPB端面及试样径向速度进行测量。图9为分别采用STFT和CWT处理结果比较。对于入射杆端面速度历程,STFT和CWT均能较好地解析,两种方法得到了较为一致的结果。对于试样径向速度历程,由于测得干涉条纹具有较大的背景噪声,而CWT对噪声较为敏感,因此不能准确地处理初始段的速度信息。STFT则显示了较好的抗噪声干扰能力,能够基本准确地解析出速度历史。通常试样径向速度在2m/s的低速范围内。就本文作者所知,目前尚未有文献报道出对SHPB试样径向速度的测量结果。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种针对具有低速且高加速度特征的质点速度测量方法,包括:
步骤1,光纤激光多普勒干涉测速***对速度为100~101m/s量级,并且加速度在109m/s2以上量级的运动靶面质点速度进行测量,输出差拍干涉信号至速度计算***;
步骤2,所述速度计算***采用短时傅里叶变换或连续小波变换对所述差拍干涉信号进行时频变换数据处理,以实现具有上述特征的运动靶面质点速度的高时间分辨率和速度分辨率测试;
所述光纤激光多普勒干涉测速***包括:激光器、光纤环形器、聚焦透镜、探测器和示波器,所述步骤1包括:
激光器发射通讯激光至光纤环形器的第一端口;
所述光纤环形器将接收到的通讯激光从其第二端口发送至聚焦透镜;
所述聚焦透镜将接收到的通讯激光的一部分作为参考光从其端面反射回所述光纤环形器的第二端口,另一部分输出至运动靶面;并将从运动靶面反射回的信号光发送至光纤环形器的第二端口;
所述光纤环形器将接收到的参考光和信号光从其第三端口发送至探测器;
所述探测器探测到所述参考光和信号光产生的差拍干涉信号,并发送至示波器;
所述示波器记录所述差拍干涉信号,并将所记录的差拍干涉信号发送至速度计算***。
2.如权利要求1所述的方法,其特征在于,
在所述步骤2中,所述速度计算***针对运动靶面的质点速度特征以及差拍干涉信号的特征选择短时傅里叶变换或连续小波变换进行时频变换数据处理;其中,对于质点速度波动剧烈,信噪比大于5dB的差拍干涉信号,选择连续小波变换进行数据处理;对于质点速度波动较为平稳,信噪比小于5dB的干涉信号,选择短时傅里叶变换进行数据处理;
所述质点速度波动剧烈是指:质点速度波动的每个周期中,vmin≤0.2vmax;
所述质点速度波动较为平稳是指:质点速度波动的每个周期中,vmin≥0.8vmax;
其中,vmin为波谷对应的质点速度,vmax为波峰对应的质点速度。
3.如权利要求1或2所述的方法,其特征在于,
采用连续小波变换对所述差拍干涉信号进行时频变换数据处理包括:
对所述差拍干涉信号进行滤波处理;
分别对预处理后的差拍干涉信号和原始的差拍干涉信号进行小波变换,得到它们各自的小波系数,再把对应的小波系数相乘;
用局部模极大值法进行小波脊线的提取,得到质点速度曲线。
4.如权利要求3所述的方法,其特征在于,
在对所述差拍干涉信号进行滤波处理的步骤之前,所述方法还包括:
对所述差拍干涉信号进行短时傅里叶变换处理,并根据得到的结果估算质点速度的上限和下限,并计算小波基的中心频率及带宽,然后根据短时傅里叶变换的时频图,选择最优的小波基;
在对所述差拍干涉信号进行滤波处理的步骤中,根据估算的质点信号速度的上限和下限,设定滤波速度的上限和下限;
在分别对预处理后的差拍干涉信号和原始的差拍干涉信号进行小波变换的步骤中,采用所述最优的小波基进行小波变换。
5.如权利要求4所述的方法,其特征在于,
在用局部模极大值法进行小波脊线的提取,得到质点速度曲线的步骤中,找到时频谱的模极大值,作为可靠的小波脊线的时间起点和差频频率起点,再根据加速度的变化范围,设置差频频率的搜寻范围,求出局部极大值,找到时间范围内各时刻的差频频率,从而实现小波脊线的提取。
6.一种针对具有低速且高加速度特征的质点速度测量装置,其特征在于,包括光纤激光多普勒干涉测速***和速度计算***,
所述光纤激光多普勒干涉测速***用于对速度为100~101m/s量级,并且加速度在109m/s2以上量级的运动靶面质点速度进行测量,输出差拍干涉信号至速度计算***;
所述速度计算***用于采用短时傅里叶变换或连续小波变换对所述差拍干涉信号进行时频变换数据处理,以实现具有上述特征的运动靶面质点速度的高时间分辨率和速度分辨率测试;
所述光纤激光多普勒干涉测速***包括:激光器、光纤环形器、聚焦透镜、探测器和示波器,其中,激光器、光纤环形器、探测器和示波器依次相连,聚焦透镜与光纤环形器相连;
所述激光器用于发射通讯激光至光纤环形器的第一端口;
所述光纤环形器用于将接收到的通讯激光从其第二端口发送至聚焦透镜,以及,将接收到的参考光和信号光从其第三端口发送至探测器;
所述聚焦透镜用于将接收到的通讯激光的一部分作为参考光从其端面反射回所述光纤环形器的第二端口,另一部分输出至运动靶面;并将从运动靶面反射回的信号光发送至光纤环形器的第二端口;
所述探测器用于探测所述参考光和信号光产生的差拍干涉信号,并发送至示波器;
所述示波器用于记录所述差拍干涉信号,并将所记录的差拍干涉信号发送至速度计算***。
7.如权利要求6所述的装置,其特征在于,
所述光纤激光多普勒干涉测速***还包括光电转换模块、显示模块、第一分光模块和第二分光模块;所述聚焦透镜、光纤环形器、探测器和第二分光模块一一对应,均为N个,N大于等于1;
所述第一分光模块的入端与所述激光器相连,所述第一分光模块的N个出端分别与N个光纤环形器的第一端口相连,将所述激光器发射的通讯激光均分给N个光纤环形器;
每个第二分光模块的入端分别与对应的光纤环形器的第三端口相连,所述第二分光模块的一个出端与对应的探测器相连,另一个出端与光电转换模块相连;
所述光电转换模块用于将接收到的光信号转换成电信号,并发送给显示模块进行显示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210216938.8A CN102721824B (zh) | 2012-06-27 | 2012-06-27 | 针对具有低速且高加速度特征的质点速度测量方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210216938.8A CN102721824B (zh) | 2012-06-27 | 2012-06-27 | 针对具有低速且高加速度特征的质点速度测量方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102721824A CN102721824A (zh) | 2012-10-10 |
CN102721824B true CN102721824B (zh) | 2014-09-24 |
Family
ID=46947643
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210216938.8A Active CN102721824B (zh) | 2012-06-27 | 2012-06-27 | 针对具有低速且高加速度特征的质点速度测量方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102721824B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085601A (zh) * | 2018-09-18 | 2018-12-25 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于弹道靶的高速模型速度连续测量装置及方法 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104316160B (zh) * | 2014-11-20 | 2017-03-01 | 哈尔滨工业大学 | 基于小波脊的水下声信号瞬时频率解调方法 |
CN107064539A (zh) * | 2017-04-27 | 2017-08-18 | 中国工程物理研究院流体物理研究所 | 一种大视场光子多普勒测速装置及方法 |
CN107328952B (zh) * | 2017-05-05 | 2019-06-18 | 四川大学 | 一种无真值时航迹速度稳定度评价方案 |
CN109387853A (zh) * | 2017-08-10 | 2019-02-26 | 中国石油化工股份有限公司 | 一种实验质点运动参数跟踪提取***及方法 |
CN107806933B (zh) * | 2017-12-08 | 2023-06-13 | 中国工程物理研究院激光聚变研究中心 | 光学材料激光诱导冲击波波速的测量装置及其方法 |
CN114545018B (zh) * | 2022-02-22 | 2024-01-30 | 中国工程物理研究院总体工程研究所 | 一种光纤破片测速装置及测速方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101726354A (zh) * | 2009-12-16 | 2010-06-09 | 中国科学院半导体研究所 | 光纤激光矢量水听器 |
CN102095450A (zh) * | 2010-12-09 | 2011-06-15 | 中国工程物理研究院流体物理研究所 | 一种单发次多信息量诊断靶 |
CN102419442A (zh) * | 2011-09-06 | 2012-04-18 | 中国工程物理研究院流体物理研究所 | 一种双源激光干涉测速*** |
-
2012
- 2012-06-27 CN CN201210216938.8A patent/CN102721824B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101726354A (zh) * | 2009-12-16 | 2010-06-09 | 中国科学院半导体研究所 | 光纤激光矢量水听器 |
CN102095450A (zh) * | 2010-12-09 | 2011-06-15 | 中国工程物理研究院流体物理研究所 | 一种单发次多信息量诊断靶 |
CN102419442A (zh) * | 2011-09-06 | 2012-04-18 | 中国工程物理研究院流体物理研究所 | 一种双源激光干涉测速*** |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085601A (zh) * | 2018-09-18 | 2018-12-25 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于弹道靶的高速模型速度连续测量装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102721824A (zh) | 2012-10-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102721824B (zh) | 针对具有低速且高加速度特征的质点速度测量方法及装置 | |
EP3014222B1 (en) | Improvements in fibre optic distributed sensing | |
CN106052842B (zh) | 可消衰落噪声的分布式光纤振动传感***及其解调方法 | |
JP6552983B2 (ja) | ブリルアン散乱測定方法およびブリルアン散乱測定装置 | |
US9335305B2 (en) | Ultrasonic measurement | |
US6504612B2 (en) | Electromagnetic wave analyzer | |
JP6277147B2 (ja) | 光ファイバ振動測定方法及びシステム | |
CN102778256B (zh) | 一种针对强激光驱动的冲击效应试验的多物理场测量*** | |
CN104390935A (zh) | 太赫兹波段测试非线性极化系数和吸收系数的装置及方法 | |
WO2008097288A2 (en) | Laser velocimetry system | |
JP2017533433A (ja) | 弾性および非弾性光散乱を利用したドップラー効果測定のためのシステムと方法 | |
CN106093962B (zh) | 一种干涉测速***及方法 | |
Song et al. | Measurement of fast-changing low velocities by photonic Doppler velocimetry | |
CN113311074A (zh) | 一种薄膜杨氏模量的确定方法、装置及*** | |
WO2014067184A1 (zh) | 基于四象限探测器的脉冲气体激光器腔内流场测量装置 | |
Culshaw et al. | Non-contact measurement of the mechanical properties of materials using an all-optical technique | |
Yang | Temporal coherence of sound transmissions in deep water revisited | |
CN111896787A (zh) | 一种用于太赫兹脉冲辐射器辐射波形的测量***及测量方法 | |
CN109031341B (zh) | 一种使用连续调频激光雷达装置的物体运动速度测量方法 | |
Kuroyama et al. | Characterization as measurement sound source of acoustic cavitation noise from bubble clusters under ultrasonic horn | |
CN114112132B (zh) | 一种激光超声测量梯度残余应力的***和方法 | |
Totems et al. | Advanced signal processing methods for pulsed laser vibrometry | |
Mokrushin et al. | An interferometer with time-and-frequency signal compression for studying properties of materials in shock wave experiments | |
Sorazu et al. | Obtaining complementary Lamb wave dispersion information by two signal processing methods on an all-optical non-contact configuration | |
US10247659B2 (en) | Device for characterizing an interface of a structure and corresponding device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |