Optical System Design of Suspended Infrared Night Vision Based on Low Light Level Helmet Observation
-
摘要: 头盔夜视仪由单波段向多波段图像融合的方向发展。本文对基于微光头盔观察、悬挂式红外夜视仪的技术方案、图像配准精度进行分析并进行光学仿真。首先分析悬挂式红外夜视仪与微光头盔组合使用的工作模式以及图像旋转、圆形视场的设计方案；其次根据悬挂式红外夜视仪的设计指标，对其红外物镜及投影物镜进行光学仿真；第三从悬挂精度、光轴一致性及畸变等三方面分析图像配准精度；最后根据仿真结果及图像配准精度分析说明基于微光头盔观察、悬挂式红外夜视仪的技术方案可行，能达到预期的效果。Abstract: Helmet night vision systems are developed from single-band to multi-band image fusion. In this study, we analyzed the technical program and image registration accuracy based on low-light-level helmet observation and a hanging infrared night vision device. Optical simulation analysis was also conducted. First, we analyzed the working mode of the combination of hanging infrared night vision and low-light-level helmet, as well as the design scheme of image rotation and circular field of view. Second, according to the design index of hanging infrared night vision, optical simulation of an infrared lens and projection lens was carried out. Third, the image registration accuracy was analyzed from three viewpoints: suspension accuracy, optical axis consistency, and distortion. Finally, according to the simulation results and image registration accuracy analysis, a technical scheme based on low-light-level helmet observation and suspended infrared night vision is feasible and can achieve the targeted effect.
-
Keywords:
- projection lens /
- infrared lens /
- image fusion /
- image registration /
- suspended /
- infrared night vision
-
0. 引言
机载三轴稳定平台广泛应用于遥感成像领域，其主要作用就是隔离载机运动等因素带来的扰动，减小飞行器姿态角带来的变化，保证成像质量[1]。现需研制一台具有高精度的三轴稳定平台，使搭载在平台上的垂视相机和高光谱相机两种线阵相机工作在稳定的成像环境中，获得高质量的图像数据。
传统的机载三轴平台机械传动往往采用电机直驱或带有减速箱的传动机构。这些方式均不能满足成像载荷质量大，稳定平台质量小于59.8 kg的指标要求。钢丝绳传动具有高精度、低空回、重量轻等特点，所以采用了钢丝绳传动机构。但是钢丝绳传动会带来一定的柔性传动，在线阵相机成像这种高精度的场合往往不能忽略，会带来图像的抖动和畸变。传统的PID（proportion integration differentiation）控制很难满足高成像质量的要求，张洪文针对某遥感扫描组件消旋组件同步提出双电机驱动，采用独立的位置环和速度环控制钢丝绳传动，同步精度达到了0.5 mrad[3]；谢宏伟在实验中确定了预紧力、负载、绳组数、转速对传动精度的影响[4]。但是以上稳定平台的成像载荷均比较小，没有明显暴露出钢丝绳传动的缺点，在这次实际工程中，两个成像载荷的重量高达21.3 kg，稳定精度的技术指标为0.35mrad，对控制算法的设计提出更高的要求。
本文主要是为了改善钢丝绳传动的三轴转台非线性干扰，设计出一种模糊自适应和前馈复合控制的策略，最后通过仿真和实验验证该控制策略的有效性。
1. 控制系统的数学模型
1.1 钢丝绳传动系统的建模
本实验的稳定平台都是小幅有限角运动，钢丝绳机构如图 1所示，结构图如图 2所示，采用ROT-LOK型结构[3-4]。在往复运动中，钢丝绳上不断变化的张力会对传动系统引入传动空回。
根据结构的空间和刚性选择4组绳组数，被动轮上每组钢丝绳的一端都串联预紧弹簧用以调节钢丝绳的预紧力，可以在实验中精确地调节弹簧的长度来选取合适的预紧力[5-6]。钢丝绳上的预紧力为Fp，负载力为Fl，钢丝绳的变形量为δ，被动轮与主动轮的半径是Ro, Ri，AE是刚度即横截面积与弹性模量的乘积，μ为摩擦系数，钢丝绳的模型为：
$$ {\delta _{\text{o}}} = \int {\frac{{{R_{\text{o}}}({F_{\text{p}}} + {F_{\text{l}}}){{\text{e}}^{ - \mu \theta }}}}{{{\text{AE}}}}} {\text{d}}\theta - \int {\frac{{{R_{\text{o}}}({F_{\text{p}}} - {F_{\text{l}}}){{\text{e}}^{\mu \theta }}}}{{{\text{AE}}}}{\text{d}}\theta } $$ (1) $$ {\delta _{\text{i}}} = \int {\frac{{{R_{\text{i}}}({F_{\text{p}}} + {F_{\text{l}}}){{\text{e}}^{ - \mu \theta }}}}{{{\text{AE}}}}} {\text{d}}\theta - \int {\frac{{{R_{\text{i}}}({F_{\text{p}}} - {F_{\text{l}}}){{\text{e}}^{\mu \theta }}}}{{{\text{AE}}}}{\text{d}}\theta } $$ (2) 则钢丝绳的传动空回e为：
$$ e = \frac{{{\delta _{\text{o}}}}}{{{R_{\text{o}}}}} + \frac{{{\delta _{\text{i}}}}}{{{R_{\text{i}}}}} $$ (3) 1.2 直流电机模型
机载三轴稳定平台会涉及到多框架耦合等问题，但本文研究的是钢丝绳传动对伺服造成的影响，所以等效为3个相同的控制系统进行设计，由电机的工作原理可以推出电压平衡方程[7]：
$$ {U_{\text{d}}} - {K_{\text{i}}}{I_{\text{d}}} = L{\dot I_{\text{d}}} + R{I_{\text{d}}} + {C_{\text{e}}}{\dot \theta _m} $$ (4) 式中：Ud为电枢两端电压；L、R、Id分别为电枢回路等效的总电感、电阻、电流；Ki为电流反馈系数；Ce是伺服电机电动势系数；$ {\dot \theta _{\text{m}}} $和θm是电机输出轴的转速和转角。
直流伺服电机的动力学方程可以写为：
$$ {J_{\text{m}}}{\ddot \theta _{\text{m}}} = {k_{\text{t}}}i - {b_\text{m}}{\dot \theta _{m} } - {K_{\text{L}}}\left( {{\theta _{\text{m}}} - {\theta _{\text{L}}}} \right) $$ (5) 式中：kti为电磁转矩；kt为电机的转矩系数；电动机输出轴的转动惯量为Jm；电机的粘性阻尼系数为bm；KL为电机和框架的耦合刚度系数；$ {\dot \theta _{\text{L}}} $和θL为负载的转速和转角。此系统的负载的转动惯量大以及精度要求高，所以将电机和惯性负载视为二质量系统进行建模分析[8]。设负载的黏性阻尼系数为bL，TmL为负载输出力矩，可以得知负载的数学模型如下：
$$ {J_{\text{L}}}{\ddot \theta _{\text{L}}} = {T_{{\text{mL}}}} - {b_{\text{L}}}{\dot \theta _{\text{L}}} $$ (6) $$ {K_{\text{L}}}\left( {{\theta _{\text{m}}} - {\theta _{\text{L}}}} \right) - {T_{{\text{mL}}}} = 0 $$ (7) 在仿真中将脉冲宽度调制（Pulse Width Modulation，PWM）环节近似的等效为一个比例环节，置为ku；速度环放大系数为ks；反馈系数为Kv。
1.3 建立摩擦模型
在伺服系统中，摩擦模型的建立相当重要，根据相关文献，仿真建立了较为经典的Lugre模型[9]。设状态变量z代表接触面鬓毛的平均变形，Lugre的数学模型可由下式表示：
$$ F = {\sigma _0}z + {\sigma _1}\dot z + \alpha \dot \theta $$ (8) $$ \dot z = \dot \theta - ({\sigma _0}\left| {\dot \theta } \right|/g\left( {\dot \theta } \right))z $$ (9) $$ g\left( {\dot \theta } \right) = {F_{\text{c}}} + \left( {{F_{\text{s}}} - {F_{\text{c}}}} \right){{\text{e}}^{ - {{\left( {\frac{{\dot \theta }}{{{v_{\text{s}}}}}} \right)}^2}}} + \alpha \dot \theta $$ (10) 式中：σ0与σ1为动态摩擦参数，其中σ0为刚毛的刚性系数，σ1为滑动阻尼系数。Fc、Fs、α、vs为静态摩擦参数，其中Fc为库仑摩擦；Fs为静摩擦；α为粘性摩擦系数；vs为stribeck特性阶段的动静切换速度，而本文针对项目情形进行了实测，得出相关参数的具体数值。
2. 控制器设计
2.1 模糊自适应控制设计
由于机载平台的特点，传统PID控制难以满足其精度要求，实验设计了模糊自适应PID控制器。其结构如图 3所示。
从结构图可以看出模糊控制采用两输入三输出二维模糊控制模块，以摆扫速度误差e和误差变化率de为输入语言变量，以Δkp, Δki, Δkd为输出语言变量。对于每个状态的3个参数，有：
$$ k_{\mathrm{p}}(t)=k_{\mathrm{p}}(t-1)+\Delta k_{\mathrm{p}}(t) $$ (11) $$ k_{\mathrm{i}}(t)=k_{\mathrm{i}}(t-1)+\Delta k_{\mathrm{i}}(t) $$ (12) $$ k_{\mathrm{d}}(t)=k_{\mathrm{d}}(t-1)+\Delta k_{\mathrm{d}}(t) $$ (13) 输入变量e和de的模糊集合可以划分为{NB, NM, NS, Z, PS, PM, PB}，输出变量Δkp, Δki, Δkd的论域划分和输入变量的划分是一致的。输入输出变量边缘采用高斯隶属函数曲线，其余部分采用三角形隶属函数曲线如式(14)：
$$ f(x, a, b, c) = \left\{ {\begin{array}{*{20}{c}} 0&{x \leqslant 0} \\ {\frac{{x - a}}{{b - a}}}&{a \leqslant x \geqslant 0} \\ {\frac{{c - x}}{{c - b}}}&{b \leqslant x \leqslant c} \\ 0&{x \geqslant c} \end{array}} \right. $$ (14) e和de分别有7个论域，故一共有49条控制规则，可以根据相应领域专家的技术和经验编写控制规则表[10-12]。
2.2 前馈补偿控制设计
前馈控制的引入就可以直接将干扰引入到控制装置，在控制系统中同时使用前馈控制加反馈控制的方法称为复合控制[12-14]。
实际钢丝绳传动会引起一定的滞后，前馈补偿基于不变性原理，通过前馈补偿环节把可以测量的扰动误差引入到控制器的设计中，从而能很好地抑制干扰，减小系统的误差，控制模型如图 4所示。
采用基于输入的自适应前馈补偿来提高伺服系统的性能，可以求得：
$$ {G_1}(s) = \frac{{Ls + R + {K_i}}}{{{G_2}(s)}} $$ (15) $$ \begin{array}{l} {G_2}(s) = \{ {J_{\text{L}}}{J_{\text{m}}}L{s^4} + [{J_{\text{L}}}{J_{\text{m}}}(R + {K_{\text{i}}}) + \hfill \\ \quad ({J_{\text{L}}}{b_{\text{m}}} + {J_{\text{m}}}{b_{\text{L}}})L]{s^3} + [({J_{\text{L}}}{b_{\text{m}}} + {J_{\text{m}}}{b_{\text{L}}})(R + {K_{\text{i}}}) \hfill \\ \quad + {b_{\text{m}}}{b_{\text{L}}}L]{s^2} + {b_{\text{m}}}{b_{\text{L}}}(R + {K_{\text{i}}})s\} {k_{\text{u}}}{k_{\text{t}}}{k_{\text{s}}}{K_{\text{v}}} \hfill \\ \end{array} $$ (16) $$ {G_{\text{r}}}(s) = \frac{{{G_2}(s)}}{{Ls + R + {K_{\text{i}}}}} $$ (17) 式中：G1为电机机电部分和负载的等效传递函数；Gr为前馈补偿装置的传递函数。
3. 三轴转台模糊自适应前馈控制仿真
3.1 三轴转台的参数分析
在本次实验中，所用电机参数如下：L＝2.2 mH，R＝10 Ω，kt＝2.1 N·m/A，Ce＝2.1 V/(rad/s)，Jm＝0.0037 kg·m2，JL＝7.91 kg·m2，bm＝0.01，bL＝15.0，Ki＝0.001，KL＝6.0；摩擦模型参数：σ0＝2600，σ1＝25，α＝0.2，Fc＝5.2，Fs＝6.4，vs＝0.01；钢丝绳选用的是SUS304材质，直径1 mm，它的实际参数是μ＝0.0015，AE＝8×103；在调试平台时，测得预紧力Fp＝21.3 N效果最好；PWM环节放大倍数ku＝11，速度环放大倍数ks＝6，反馈系数Kv＝2；PID控制器3个参数kp＝100，ki＝1.0，kd＝50；根据以上相关参数求出前馈部分传递函数：
$$ {G_{\text{1}}}(s) = \frac{{0.041{s^4} + 81.210{s^3} + 373.203{s^2} + 415.842s}}{{0.0022s + 10.001}} $$ (18) 前馈控制的输入是根据飞机上的主惯导获得数据，经过控制板处理之后用来补偿飞机姿态，由于条件限制，仿真的时候所用的是以前在飞机下载的惯导数据。基于传统PID控制仿真模型和模糊自适应PID前馈补偿控制仿真模型框图如图 5、图 6所示。
![]()
图 6 模糊自适应PID前馈补偿控制仿真模型Figure 6. Fuzzy adaptive PID feedforward compensation control simulation model3.2 三轴转台仿真结果分析
运行仿真得到的仿真曲线如图 7，采用传统PID控制，在速度过零点时，波形明显发生畸变，出现位置跟踪“平定”现象和速度跟踪“死区”现象，位置跟踪出现很大的滞后。而采用了模糊控制PID控制器和前馈控制符合的控制算法，消除了“平定”和“死区”现象，对系统滞后明显改善，响应速度更快。
4. 实验
4.1 摇摆台抗干扰实验
钢丝绳传动误差是制约伺服性能的关键，传动误差是指当输入轴单向转动时，输出轴的实际值对理想值的误差，在摇摆台实验室，模拟飞机环境对稳定平台进行实验测量，图 8是平台测试情况，在稳定平台横滚框架安装平面反射镜，横滚框架位于零位，反射镜与地面平行，在反射镜正下方放置自准直仪，自准直仪精度设置为1 μrad，俯仰轴以0.3°，0.5 Hz摆动。用传统PID控制和模糊自适应前馈控制算法分别测试曲线如图 9。从输出曲线可以看出，传统PID存在很大的抖动，传动精度在1mrad，而模糊自适应前馈控制传动精度达到了0.2 mrad，很大地改善了钢丝绳对系统带来的抖动，能达到指标要求。
4.2 速度稳定性实验
用平行光管模拟无穷远处，设备摆放如图 10，稳定平台俯仰框架以固定速度摆扫，相机通过反射镜对靶标进行成像，平行光管采用模块化光电测试系统METS L19型号，分辨率靶标选择使用美军标USAF1951分辨率靶。拍出的靶标如图 11所示。图 11(a)采用传统PID控制拍出来的靶标，存在抖动弯曲的现象，用Matlab查看每行相差两个像素，相机角分辨率为10 μrad，而图 11(b)采用模糊自适应前馈控制算法拍出来的靶标用Matlab查看靶标每行像素是基本一致的，满足指标要求。
4.3 实验结果分析
在摇摆台模拟飞机环境做的抗干扰实验，将控制系统的稳定精度从1 mrad改善到0.2 mrad，在速度平稳性实验中，采用平行光管测靶标成像的方式，采用复合控制的算法能较好地改善钢丝绳带来的速度抖动性，以上实验与仿真结果也相吻合，仿真和实验都说明了采用复合控制策略较传统PID对控制精度带来了大幅的改善，可以为以后类似的工程提供参考。
5. 结论
本文针对采用钢丝绳传动的大质量线阵相机的伺服控制精度，提出模糊自适应PID前馈补偿的复合控制策略，将按输入设计的前馈控制方法与模糊自适应控制结合起来，通过仿真和实验证明，该策略能够有效地消减钢丝绳传动误差，平台稳定精度从1 mrad提高到0.2 mrad，满足了成像的指标要求。
-
![]()
图 1 微光头盔与悬挂式红外夜视仪组合方式
Figure 1. Combination of low light level helmet and suspended infrared night vision
![]()
图 2 微光头盔与悬挂式红外夜视仪融合图像
Figure 2. Fusion image of low light level helmet and hanging infrared night vision
![]()
图 3 微光头盔与悬挂式红外夜视仪工作示意图
Figure 3. Working diagram of low light level helmet and suspended infrared night vision
![]()
图 10 红外夜视仪倾斜悬挂时成像关系
Figure 10. Imaging relationship of infrared night vision with tilted suspension
![]()
图 11 光轴无偏差时融合成像示意图
Figure 11. Schematic diagram of fusion imaging without optical axis deviation
![]()
图 12 光轴有偏差时融合成像示意图
Figure 12. Schematic diagram of fusion imaging with optical axis deviation
表 1 悬挂式红外夜视仪光学参数
Table 1 Optical parameters of suspended infrared night vision
Infrared lens Focal length
Field
F/#
Band
Detector type13.88 mm
20°(circular)
1
8～12 μm
UFPA 384×288, 17 μmProjection lens Focal length
Field
F/#
Band20.58 mm
20°(circular)
4
0.486～0.656 μmDetector type OLED 800×600, 12.6 μm Exit pupil distance 3.7 mm Suspended infrared night vision Field
Magnification
Temperature20°(circular)
1×
−40℃～60℃
下载: 导出CSV
表 2 红外物镜、投影物镜在相同视场点处的畸变
Table 2 Distortion of infrared lens and projection lens at the same field of view
Field of view Infrared lens Projection lens 0.1ω −0.01999022% 0.02068782% 0.2ω −0.07997603% 0.08263791% 0.3ω −0.18000125% 0.18550392% 0.4ω −0.32013343% 0.32868725% 0.5ω −0.50045517% 0.51130229% 0.6ω −0.72105160% 0.73212184% 0.7ω −0.98199306% 0.98949535% 0.8ω −1.28331225% 1.28122817% 0.9ω −1.62497467% 1.60440259% 1ω −2.00684062% 1.95510906%
下载: 导出CSV
表 3 红外物镜公差
Table 3 Tolerance of infrared objective lens parts
Parameter Tolerance N ±3 aperture ΔN ±0.8 aperture Aspheric error ±0.00006 mm Thickness of optical parts ±0.02 mm Surface tilt ±0.006 mm Air distance ±0.02 mm Element tilt ±0.02 mm Element eccentricity 0.025 mm
下载: 导出CSV
表 4 红外物镜公差分析结果
Table 4 Tolerance analysis results of infrared objective lens
Lens percentage/% MTF minimum
(Nyquist frequency)90 0.219 80 0.246 50 0.291 20 0.348 10 0.365
下载: 导出CSV
表 5 投影物镜零件公差
Table 5 Tolerance of projection lens parts
Parameter Tolerance N ±4 aperture ΔN ±0.5 aperture Thickness of optical part ±0.02 mm Air distance ±0.04 mm Surface tilt ±6′ Element tilt ±6′ Element eccentricity ±0.052 mm nd ±0.0009 vd ±0.95%
下载: 导出CSV
表 6 投影物镜公差分析结果
Table 6 Tolerance analysis results of projection lens
Lens percentage/% MTF minimum(40 lp/mm) 90 0.603 80 0.626 50 0.669 20 0.687 10 0.697
下载: 导出CSV
-
[1] 任桃桃, 邱亚峰. 头盔式单目低照度CMOS夜视仪结构设计与分析[J]. 红外技术, 2016, 38(8): 653-658. http://hwjs.nvir.cn/article/id/hwjs201608005 REN Taotao, QIU Yafeng. Structure design and analysis of monocular low illumination CMOS night vision device for helmet[J]. Infrared Technology, 2016, 38(8): 653-658. http://hwjs.nvir.cn/article/id/hwjs201608005
[2] 白瑜, 程习敏, 冯成, 等. 拼接式全景头盔红外夜视仪系统设计[J]. 光电技术应用, 2014, 29(1): 26-31. https://www.cnki.com.cn/Article/CJFDTOTAL-GDYG201401008.htm BAI Yu, CHENG Ximin, FENG Cheng, et al. Design of segmented panoramic helmet-mounted infrared night vision[J]. Electro-Optic Technology Application, 2014, 29(1): 26-31. https://www.cnki.com.cn/Article/CJFDTOTAL-GDYG201401008.htm
[3] 王思博. 三通道微光夜视仪光学系统设计[D]. 长春: 长春理工大学, 2014. WANG Sib. Optical of Three Channels of Night Vision Instrument[D]. Changchun: Changchun University of Science Technology, 2014.
[4] 史训豪. 夜视仪在单兵观瞄系统中的应用[J]. 科技风, 2019(10): 18-19. https://www.cnki.com.cn/Article/CJFDTOTAL-KJFT201928018.htm SHI Xunhao. Application of night vision device in single soldier observation and aiming system[J]. Technology Wind, 2019(10): 18-19. https://www.cnki.com.cn/Article/CJFDTOTAL-KJFT201928018.htm
[5] 李力, 李训牛. 直升机飞行员头盔夜视系统视觉及舒适性优化技术分析[J]. 红外技术, 2017, 39(10): 890-896. http://hwjs.nvir.cn/article/id/hwjs201710004 LI Li, LI Xunniu. Analysis of optimization technology for a night vision system for helicopter pilots with special focus on the visual characteristics and comfort[J]. Infrared Technology, 2017, 39(10): 890-896. http://hwjs.nvir.cn/article/id/hwjs201710004
[6] 许为. 航空夜视镜的人机工效学问题[J]. 国际航空, 1991(12): 52-53. https://www.cnki.com.cn/Article/CJFDTOTAL-ZHHK200504027.htm XU Wei. Ergonomic problems of aerial night vision goggles[J]. International Aviation, 1991(12): 52-53. https://www.cnki.com.cn/Article/CJFDTOTAL-ZHHK200504027.htm
-
期刊类型引用(0)
其他类型引用(2)
计量
- 文章访问数: 194
- HTML全文浏览量: 44
- PDF下载量: 83
- 被引次数: 2
