从傅里叶望远镜接收信号中获取目标速度的方法

李杨; 相里斌; 张文喜

doi:10.3788/HPLPB20132509.2193

从傅里叶望远镜接收信号中获取目标速度的方法

doi: 10.3788/HPLPB20132509.2193

李杨^1,2, ,,
相里斌³,
张文喜^2,3

1.
中国科学院西安光学精密机械研究所, 西安 71 01 1 9;
2.
中国科学院大学, 北京 1 00049;
3.
中国科学院光电研究院, 北京 1 00094

详细信息

通讯作者:
李杨

计量
- 文章访问数: 1524
- HTML全文浏览量: 290
- PDF下载量: 220
- 被引次数: 5
出版历程
- 收稿日期: 2013-01-04
- 修回日期: 2013-04-26
- 刊出日期: 2013-07-17

Target velocity acquisition from signals received by Fourier telescopy

Li Yang^{1,2
, ,},
Xiang Libin³,
Zhang Wenxi^2,3

1.
Xi’an Institute of Optics and Precision Mechanics,Chinese Academy of Sciences,Xi’an 710119,China;
2.
Graduate University of Chinese Academy of Sciences,Beijing 100049,China;
3.
Academy of Opto-Electronics,Chinese Academy of Sciences,Beijing 100094,China

摘要

摘要: 分析了傅里叶望远镜接收信号，结果表明，接收信号中不仅包含了被测目标的图像信息，还包含了目标的二维速度矢量信息。通过提取接收周期信号的频率可以推算出目标各分量的速度，从而得到二维速度矢量。信号频率的提取采用快速傅里叶变换分析进行频率粗估计，频率精确校正采用加Hanning窗的比值校正法。仿真了利用傅里叶望远镜测量导弹飞行速度。结果表明：在无噪声影响下，利用该方法获取的目标速度几乎无误差；存在噪声时，若接收信噪比高于5 dB，测量的目标速度精度高于0.02%。说明从傅里叶望远镜接收信号中获取目标速度的方法精度高、抗噪性好，理论上是可行的。
- 傅里叶望远镜 /
- 测速 /
- 频谱校正
Abstract: The signals received by Fourier telescopy are analyzed and demonstrated to contain not only image information but also two-dimensional velocity information of the measured target. The velocities on each component can be calculated by the periodic signal frequencies that can be extracted from received signals. Then the two-dimension velocity is composited. The frequency of the signal is extracted by fast Fourier transform analysis for coarse frequency estimate, and by ratio frequency correction algorithms using Hanning window for precise frequency correction. Simulation is made to measure the flight velocity of the missile by Fourier telescopy. The result shows that the target velocity acquired by this method has almost no error without noise, while it has relatively high precision with noise when SNR is higher than 5 dB. The velocity acquisition method from receiving signal of Fourier telescopy is demonstrated to have high precision, a good anti-noise property and to be theoretically feasible.
- Fourier telescopy /
- velocity measurement /
- spectrum correction

HTML全文

激光自1960年产生以来，便得到了迅速的发展和广泛的使用。在工程领域，金属材料的焊接^[1-2]、切割、熔覆^[3-8]等都离不开激光。激光与金属材料相互作用过程中，金属材料吸收的光能转化为热能，导致金属材料的温度升高。激光加工具有作用时间短、温度变化速度快等特点，除此之外，激光所产生的温度与金属材料的材质有关。因此，研究激光作用于金属材料产生的瞬态温度场分布和激光的中心温度对激光加工有十分重要的意义。目前，已经有诸多的实验与理论工作，例如：毕娟等^[9]利用二维轴对称的物理模型，采用有限元方法模拟了双层材料的瞬态温度场，通过数值模拟计算了激光的瞬态温度场；周杰等^[10]基于脉冲激光辐照金属材料的热学模型，通过热传导方程实现了对激光辐照下瞬态温度场的预测。激光的中心温度的测量一直是一个重要而困难的问题，到目前为止仍没有得到很好的解决。本文利用XL-800WF光纤传输多功能脉冲激光加工系统，在输入电流为150 A、脉冲宽度为2.5 ms的情况下，照射一块较大面积的厚钢板，拍摄了钢板的凝固显微照片。通过建立无热源的热传导方程，研究理想情况下激光中心的温度场，模拟了激光温度场的分布图，推导了激光中心温度的数值范围，并对可能存在的误差进行了分析。

1. 实验

实验上，采用XL-800WF光纤传输多功能脉冲激光加工系统。其最大输出功率800 W，激光波长1064 nm，有效光斑尺寸约1.0 mm，单脉冲最大激光能量90 J，脉冲宽度0.5~5 ms，脉冲频率1~40 Hz。采用光纤传输，聚焦透镜的焦距为160 mm，输入电流为150 A，脉冲持续时间为2.5 ms，设备配置了CCD监视系统。实验使用的钢板表面光洁、平整，表面积足够大，热扩散效果良好。钢板较厚，在激光照射过程中不会熔穿。激光对钢板的一个固定点进行照射，凝固后形成激光熔斑。清除熔斑周围的残渣，拍摄其显微图片。图 1 (a), (b), (c)分别为激光照射1, 2, 3个脉冲的熔斑显微图。由图可以获得清晰的液态钢流动路径。测量结果表明，单次脉冲照射和多次脉冲照射得到的凝固图极为相似，凝固圆的外直径也相同。由此可知，激光器输出稳定，激光每一次照射均具有相同的能量，实验数据可靠。

图 1 激光脉冲分别照射1, 2, 3次的显微图

Figure 1. Micrographs of laser pulses illuminating 1, 2, 3 times, respectively

下载: 全尺寸图片幻灯片

2. 瞬态温度场模拟

由于照射的钢板较大，不影响激光热量的扩散。在数学上可以认为是自由边界。以激光熔斑中心为原点建立直角坐标系，r表示任意点的位置，激光作用结束的瞬间为初始时刻。理想情况的热传导方程可写为

$\frac{\partial u(r, t)}{\partial t}=a^2 \frac{\partial^2 u(r, t)}{\partial^2 r^2}$

(1)

$u(r, 0)=\phi(r)$

(2)

式中：u(r, t)表示任意位置、任意时刻的温度；t表示时间；a为热扩散系数；ϕ(r)是初始时刻的温度分布。对方程(1)、(2)进行傅里叶变换，得

$\frac{\mathrm{d} U(\omega, t)}{\mathrm{d} t}=-a^2 \omega^2 U(\omega, t)$

(3)

$U(\omega, 0)=\varPhi(\omega)$

(4)

式中：U(ω, t)为u(r, t)作傅里叶变换的对应表达式；ω为r的对应算符。其中方程(4)中的Φ(ω)为

$\varPhi(\omega)=\int_{-\infty}^{\infty} \phi(r) \mathrm{e}^{-\mathrm{i} \omega r} \mathrm{~d} r$

(5)

式中：ϕ(r)为积分核。联合方程(3)~(4)，解得

$U(\omega, t)=\varPhi(\omega) \mathrm{e}^{-a^2 \omega^2 t}$

(6)

将方程(6)作傅里叶逆变换，得

$u(r, t)=\frac{1}{2 {\rm{ \mathsf{ π}}}} \int_{-\infty}^{\infty} \phi(r) \mathrm{e}^{-\omega^2 a^2 t} \mathrm{e}^{\mathrm{i} \omega t} \mathrm{~d} \omega$

(7)

将方程(5)代入方程(7)，得

$u(r, t)=\frac{1}{2 {\rm{ \mathsf{ π}}}} \int_{-\infty}^{\infty}\left(\int_{-\infty}^{\infty} \phi(\zeta) \mathrm{e}^{-\mathrm{i} \omega \zeta} \mathrm{d} \zeta\right) \mathrm{e}^{-\omega^2 a^2 t} \mathrm{e}^{\mathrm{i} \omega r} \mathrm{~d} \omega$

(8)

式中：ζ为积分变量。将方程(8)进一步化简，得

$u(r, t)=\frac{1}{2 a \sqrt{{\rm{ \mathsf{ π}}} t}} \int_{-\infty}^{\infty} \phi(\zeta) \mathrm{e}^{-\frac{(r-\zeta)^2}{4 a^2 t}} \mathrm{~d} \zeta$

(9)

方程(9)是一个卷积形式，利用卷积性质将其写成如下形式

$u(r, t)=\phi(r) * K(r, t)$

(10)

式中: $K(r, t)=\frac{1}{2 a \sqrt{{\rm{ \mathsf{ π}}} t}} \mathrm{e}^{-\frac{r^2}{4 a^2 t}}$ 。实验使用的XL-800WF光纤传输脉冲激光加工系统，其输出光束为高斯型光束。由于光束的线度远小于光热的影响范围，故用δ函数来表示，即ϕ(r)=δ(r)。因此，计算方程(9)得

$u(r, t)=\frac{1}{2 a \sqrt{{\rm{ \mathsf{ π}}} t}} \mathrm{e}^{-\frac{r^2}{4 a^2 t}}$

(11)

方程(11)即为理想情况下激光的中心温度演化函数。结合光与物质相互作用的反应时间和激光熔斑的凝固时间，即可推导、模拟出激光的中心温度。

3. 理论分析

测量光与物质相互作用的反应时间和激光熔斑的凝固时间对瞬态温度场的模拟有重要影响，亦是解决瞬态温度场问题的核心所在。尽管这方面已经有诸多的工作^[11-14], 但仍然没有形成统一的认识和标准。目前，反应时间及凝固时间的获取较为成熟的方法有高速摄像^[15-16]、数值模拟^[17-18]等。其中高速摄像技术获得的数值更为精确，但其设备昂贵，对实验室的要求较高。相对而言，数值模拟更具通用性和实用性。本实验无法实现直接测量，反应时间和凝固时间均通过间接方法得到。同时，本文给出了一个简单且合理的推断，以证明数据的合理性。

反应时间是激光作用于金属材料后，光与物质相互作用需要一个转换时间，使得激光的能量为金属材料吸收并产生热量，转化为金属材料的温度场。通过参考相关的数据，并结合多年的实验经验，反应时间取t=1.0 ms。凝固时间是指金属获得自身的温度场后，熔化为液态，从液态凝固至固态的时间。凝固时间取Δt=0.8 ms。另外，为了说明凝固时间的合理性，凝固时间可以通过白光干涉仪的测量结果及金属的流速来模拟解决。

用白光干涉仪测量熔斑的外半径，测量的结果为r=0.4 mm(图 2)。液态钢的流速^[19]一般可以取为v=1.0 m/s。由于凝固时液态金属的粘滞阻力不断增加，流速会逐渐变缓，是一个减速过程，可以用匀减速过程来模拟液态钢的流动，由此可得平均速度为v=0.5 m/s。因此，可计算激光熔斑的凝固时间为

$\Delta t=\frac{r}{\bar{v}}=0.8 \mathrm{~ms}$

(12)

图 2 白光干涉仪测量激光熔斑的半径

Figure 2. Radius of laser welding spots measured by white light interferometer

下载: 全尺寸图片幻灯片

由以上分析模拟可知，钢板的反应时间为t=1.0 ms。激光熔斑从液态转变为固态的凝固的时间为0.8 ms, 整个过程在1.8 ms时刻完成(如图 3所示)。

图 3 钢板的瞬态温度图。其中t=1.0 ms为熔化时的温度，t=1.8 ms为凝固时的温度

Figure 3. Transient temperature diagram of steel plate. Where t=1.0 ms corresponds to temperature of melting time, t=1.8 ms corresponds to temperature of solidification

下载: 全尺寸图片幻灯片

将时间t=1.0 ms和t=1.8 ms代入方程(11)，其比值即为凝固时间内激光中心温度的衰减率

$\eta=\frac{u(r, 1.8)_{\max }}{u(r, 1.0)_{\max }} \times 100 \%=\frac{6.31}{8.92} \times 100 \%=70.74 \%$

(13)

方程(13)表明，在0.8 ms内，激光的中心温度减为原来的70.74%。最后，通过纯钢熔点(凝固温度)逆推激光的中心温度。钢的熔点(凝固温度)为1811 K^[20]。由此可以推导激光的中心瞬时温度为

$u(r, 1) \approx \frac{1811 \mathrm{~K}}{\eta}=\frac{1811 \mathrm{~K}}{70.74 \%}=2560 \mathrm{~K}$

(14)

激光作用于金属材料瞬态温度场模拟的类似工作有很多。王馨露^[11]利用非线性瞬态热传导微分方程来模拟激光打孔温度场，并利用ABAQUS有限元分析方法进行计算，得到的结果与本文相吻合。本质上，自由热传导方程与非线性瞬态热传导微分方程对温度场的描述具有一致性。相比之下，非线性瞬态热传导微分方程侧重于金属材料的性质，而自由热传导方程的侧重点在激光加工系统；裴旭等人^[21]采用COMSOL系统研究了脉冲激光与金属铝相互作用时的瞬态温度场变化，得到了瞬态温度分布亦在同一范围内。

长期以来，研究者们在激光中心瞬态温度场和激光的中心温度探测方面做了大量工作，根据实验、经验和模拟方法得到了部分数据。将上述结果与已有的工作经过对比，本文得到的数据与以前的结果比较吻合，得到的数据在同一个温度范围。这也充分表明本文模拟方案较为可靠。

4. 结论

激光的中心温度的测量一直是一个重要而困难的问题，到目前为止仍没有得到很好的解决。本文利用XL-800WF型光纤传输多功能脉冲激光加工系统对厚钢板进行多次脉冲照射，钢板熔化后在扩散过程中凝固，形成凝固图。由白光干涉仪测量了熔斑的直径，通过参考相关文献及液态金属的流速推算金属的凝固时间。根据液态钢的凝固时间，建立了理想条件下的热传导方程，求解热传导方程来计算激光的温度，从而模拟了激光瞬态温度场分布。结果表明，激光中心的温度约为2560 K，与以前的实验结果相比较，本文的结果与其吻合得较好。

参考文献(0)

施引文献

期刊类型引用(2)

1.	李小平，赵汝和，刘言生. 数控电解机械激光多向切割温度场模拟分析. 激光杂志. 2020(05): 114-117 . 百度学术
2.	李国会，徐宏来，向汝建，杜应磊，吴晶，向振佼，张越. 金属和玻璃胶合体热膨胀仿真及实验验证. 强激光与粒子束. 2019(12): 10-15 . 本站查看

其他类型引用(3)

资源附件(0)

访问统计

点击查看大图

计量

文章访问数: 1524
HTML全文浏览量: 290
PDF下载量: 220
被引次数: 5

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

从傅里叶望远镜接收信号中获取目标速度的方法

doi: 10.3788/HPLPB20132509.2193

通讯作者:
李杨

计量

Target velocity acquisition from signals received by Fourier telescopy

1. 实验

2. 瞬态温度场模拟

3. 理论分析

4. 结论

期刊类型引用(2)

其他类型引用(3)

计量

目录

1. 实验

2. 瞬态温度场模拟

3. 理论分析

4. 结论

留言板

从傅里叶望远镜接收信号中获取目标速度的方法

doi: 10.3788/HPLPB20132509.2193

通讯作者: 李杨

计量

出版历程

Target velocity acquisition from signals received by Fourier telescopy

1. 实验

2. 瞬态温度场模拟

3. 理论分析

4. 结论

期刊类型引用(2)

其他类型引用(3)

计量

出版历程

目录

1. 实验

2. 瞬态温度场模拟

3. 理论分析

4. 结论

通讯作者:
李杨