大口径非球面铝反射镜三轴联动加工技术研究

肖建国, 姚同, 张万清, 刘尧, 罗宏, 李茂忠, 黄攀, 康杰, 张若寅

肖建国, 姚同, 张万清, 刘尧, 罗宏, 李茂忠, 黄攀, 康杰, 张若寅. 大口径非球面铝反射镜三轴联动加工技术研究[J]. 红外技术, 2021, 43(9): 829-835.
引用本文: 肖建国, 姚同, 张万清, 刘尧, 罗宏, 李茂忠, 黄攀, 康杰, 张若寅. 大口径非球面铝反射镜三轴联动加工技术研究[J]. 红外技术, 2021, 43(9): 829-835.
XIAO Jiaoguo, YAO Tong, ZHANG Wanqing, LIU Yao, LUO Hong, LI Maozhong, HUANG Pan, KANG Jie, ZHANG Ruoyin. Study on Machining Technology of Large Aspheric Aluminum Reflector with Three Axis Linkage[J]. Infrared Technology , 2021, 43(9): 829-835.
Citation: XIAO Jiaoguo, YAO Tong, ZHANG Wanqing, LIU Yao, LUO Hong, LI Maozhong, HUANG Pan, KANG Jie, ZHANG Ruoyin. Study on Machining Technology of Large Aspheric Aluminum Reflector with Three Axis Linkage[J]. Infrared Technology , 2021, 43(9): 829-835.

大口径非球面铝反射镜三轴联动加工技术研究

详细信息
    作者简介:

    肖建国(1971-), 硕士, 高级工程师, 主要从事超精密单点金刚石车床车削加工技术研究。E-mail:xiaojianguoyn@163.com

  • 中图分类号: TH706

Study on Machining Technology of Large Aspheric Aluminum Reflector with Three Axis Linkage

  • 摘要: 为解决单点金刚石车削技术(single-point diamond turning,SPDT)应用于大口径大弦高非球面铝反射镜加工中遇到的车床导轨行程、工作台回转容积受限和加工质量较低的问题,针对φ682 mm口径、弦高220 mm的凹非球面铝反射镜加工,首先,提出基于SPDT的三轴联动加工方法,在两轴加工基础上加入旋转B轴进行协同加工,使得导轨行程和工作台回转容积能够满足加工需求。然后,设计专用笼状夹具,并通过有限元法研究支撑杆数量、杆径、上下连接板厚度参数对夹具-工件形变特征的影响,通过极差和方差分析研究不同因素对夹具-工件最大变形量影响的显著性,得到一组最佳夹具设计参数组合,即夹具支撑杆数量为24,杆径为22 mm,上下连接板厚度为25 mm。最后,将铝反射镜固定在优化后的笼状夹具上,通过三轴联动加工实现了对φ682 mm口径非球面铝反射镜的加工。对调刀件的检测结果表明:调刀件面形精度Pv值为0.6 μm,表面粗糙度Ra约为10.1 nm,可认为φ682 mm口径非球面铝反射镜的面形精度和表面粗糙度均达到使用要求。本研究能够为同类大口径大弦高非球面反射镜的加工提供一定的理论基础和加工工艺技术参考。
    Abstract: To solve the problem of single point diamond turning(SPDT), technology was used for the processing of large diameter and chord height aspheric aluminum mirrors, which have the problems of limited lathe guide stroke, limited rotary volume of worktable, and low processing quality. To process a concave aspheric aluminum mirror with a diameter of φ682 mm and a chord height of 220 mm, first, a three-axis linkage processing method based on SPDT was proposed, which adds rotary b-axis on the basis of two-axis processing, such that the guide rail travel and table rotation volume can meet the processing requirements. Then, the special cage fixture was designed, and the influence of the number of supporting rods, the diameter of supporting rods, and the thickness of upper and lower connecting plates on the jig-work piece deformation characteristics were examined using the finite element method. The influence of different factors on the maximum deformation of jig and work piece was evaluated via range and variance analysis. A set of optimal jig design parameters was obtained, that is, the number of jig support rods was 24, the diameter of the rods was 22 mm, and the thickness of the upper and lower connecting plates was 25 mm. Finally, the aluminum mirror was fixed on the optimized cage clamp, and the processing of the φ682 mm aspheric aluminum mirror was realized through three-axis linkage processing. The test results show that the surface accuracy Pv of the tool adjusting part was 0.6 μm, and the surface roughness Ra was approximately 10.1 nm. It can be considered that the surface accuracy and surface roughness of the φ682 mm aspheric aluminum mirror can meet the requirements. This study can provide a theoretical basis and technological reference for the processing of the same type of large aperture and high chord aspheric mirror.
  • 智能电网的快速发展,巡检机器人已经在电力系统中大量应用,变电站巡检机器人则是其中一种常用的机器人,该类机器人的主要作用是在无人值守的情况下,对变电站的设备实行巡检[1],采集设备的运行状态信息,监测设备的运行状况,保证电力系统的运行安全。该类机器人在使用过程中,可通过两种方式进行控制,一是遥控方法,二是自主方式。其中自主方式则是当下的主要应用方法,该方式是通过巡检目标、路线等相关内容的设定实现巡检目的[2-3]。但是,由于变电站内存在多种磁场的干扰,会对机器人的巡检轨迹造成不同程度的影响,导致其巡检轨迹发生偏差,与实际轨迹之间存在误差,大规模部署高精度传感器在成本与应用上,都存在较大困难。因此,针对机器人轨迹的三维点云配准[4],很难实现目标的配准。针对这一问题,国内外学者对此提出众多研究,唐志荣等人基于因子分析法提出相关的三维点云配准方法[5];采用了一种新的因子分析法将三维点云的定标问题转化成了模型的参量;利用高斯混合模式进行三维点云的拟合,利用最大值法得到了各因素的系数加载矩阵,通过因素负载矩阵来实现点云的配准。刘剑等人基于特征匹配提出相关的三维点云配准方法[6];首先利用快速点特征直方图(fast point feature histograms, FPFH)对特征进行全面的描述,利用三角网构建特征的区域相关性;在特征点组之间匹配,进行初始化取样的一致性转换,完成三维点云配准。Kadam提出一种绿色、准确、无监督的点云配准方法[7];基于迭代最近点算法和特征配准方法的配准算法提供了准确的配准结果,并且使用一个简单的方法有效地删除注册后的点云中的冗余数据,考虑每个对应点的权重值的算法开发了一个三维应用系统以及OpenGL三维图形库,可以解决点云的配准和冗余消除问题。Fan等提出一种基于边缘匹配的点云到CAD模型的配准方法[8]。分别从点云和CAD模型中提取边缘点,点云数据的提取包括两个步骤,即剖面线数据的提取和识别。然后以相同的间隔对两组数据进行均匀采样,并根据方向角序列的偏差进行匹配,获得点云数据与CAD模型的旋转角度。在模板的三条边上进行验证实验。但是,这些方法对面临着大量的干扰数据,只能在高精度外在设备配合下,才能发挥作用。随着模式识别技术的不断发展,以红外图像特征为基础的信息采集方式逐渐成熟,可通过相关特征融合,将图像中的不同类别的特征实行融合,用于更加可靠地实现图像中内容的识别,降低技术成本。

    因此,本文提出红外图像特征融合下变电站机器人巡视轨迹三维点云配准方法,该方法以融合红外图像特征为基础,实现机器人巡检轨迹的三维点云配准。

    图像特征是由具有一定的时空相关性的图像内容,对所需选取对象的巡视轨迹是对静态或动态的图像进行有效的检索和分离;通过目标追踪技术对图像中所抽取出的目标进行追踪,并获得其巡视轨迹,其主要原理是通过对每个图像中所探测到的物体进行定位,并将它们连接在一起,通过时空关联,获得变电站机器人的运动路径,从而达到对变电站机器人进行追踪的目的。

    保证机器人巡检过程中,能够按照设定轨迹完成,本文研究红外特征融合下的变电站机器人巡检轨迹三维点云配准方法,该方法的整体框架如图 1所示。

    图  1  方法的整体框架
    Figure  1.  The overall framework of the method

    该方法首先通过红外激光仪以及标定的摄像头采集变电站的设备图像信息,然后通过特征提取方法提取采集图像中的特征后,完成特征融合;最后通过三维点云配准方法完成巡检轨迹的配准。

    HOG(histogram of oriented gradient)特征也称为方向梯度直方图特征,该特征能够良好地描述红外图像的形状、轮廓等,并且在机器人采集图像发生旋转和光照变化时,依旧能够有效地描述图像的相关特征[7]

    细胞单元(cell)作为该特征的最小单位,对其全部的像素点梯度实行求解,同时完成方向直方图的统计;除此之外,结合附近的cell形成区块,对每一个细胞单元特征向量实行处理,通过归一化拼接方式完成,其相序步骤为:

    1)特征提取过程中,为保证机器人轨迹图像的局部细节不会受到影响,需对图像实行处理[8],通过Gamma校正法完成,其公式为:

    $$ J_{Z}=I(x, y)^{1/2} (1) $$ (1)

    式中:机器人采集的图像用I(x, y)表示,且为输入;Gamma参数用1/2表示,可通过该参数判断变电站现场的光照变化,>1时,表示高光和暗光部分分别被压缩和扩展,<1时则相反。

    从两个方向逐一求解所有的梯度像素,通过一维离散微分模板完成,其中模板为[-1, 0, +1],两个方向分别为水平和垂直,在此基础上求解梯度模值和方向角,属于每一个像素,其公式为:

    $$ {p_x}\left( {x, y} \right) = \left| {G\left( {x + 1, y} \right) - G\left( {x - 1, y} \right)} \right| $$ (2)
    $$ {p_y}\left( {x, y} \right) = \left| {G\left( {x, y + 1} \right) - G\left( {x, y - 1} \right)} \right| $$ (3)
    $$ p\left( {x, y} \right) = \sqrt {{p_x}{{\left( {x, y} \right)}^2} + {p_y}{{\left( {x, y} \right)}^2}} $$ (4)
    $$ \alpha \left( {x, y} \right) = {\tan ^{ - 1}}\left( {\frac{{{p_y}\left( {x, y} \right)}}{{{p_x}\left( {x, y} \right)}}} \right) $$ (5)

    式中:py(x, y)、px(x, y)均表示梯度,前者对应水平方向,后者对应垂直方向,均属于像素点(x, y)处;p(x, y)、α(x, y)分别为梯度模值和方向角,均属于像素点。设置cell的尺寸,并对梯度方向实行划分以及统计梯度直方图,获取轨迹图像中HOG特征向量。

    局部自相似(local self-similarity, LSS)特征也称为局部自相似描述子,其能够描述红外图像局部区域的形状属性,即使机器人采集图像存在明显的颜色差异和尺寸区别[9],依旧能够较好地获取图像中和轨迹相关的几何布局和形状特征。

    以每一个像素为中心,在输入轨迹图像的局部区域内,选择尺寸合适的邻域作为子窗口,并对其与相对应的局域之间的差平方和实行计算,采用归一化方式对计算结果实行处理,将其转换成相关面[10],其归一化公式为:

    $$ {S_q}\left( {x, y} \right) = \exp \left( { - \frac{{{\rm SSD}_q\left( {x, y} \right)}}{{\max \left( {{{{\rm var} }_{\rm noise}}, {{{\rm var} }_{\rm auto}}\left( q \right)} \right)}}} \right) $$ (6)

    式中:varnoise为常数,红外图像的灰度变化通过其描述;varauto(q)表示对比度,属于子窗。

    在轨迹图像LSS特征提取时,为实现其对图像形变的最大容忍程度以及获取差异性的子区域,需对图像实行划分处理,且以区间的方式完成,同时,极角和极径则是划分的实际位置;获取对各个区域内的相关值,并选择该值中的最大值,用其描述为特征值;LSS特征向量的获取则可通过两种方法结合完成[11],即归一化和组合拼接出。

    为提升巡检轨迹的配准效果,采用多特征自适应融合方法对上述获取的HOG特征和LSS特征实行融合,将两个特征联合观测,形成红外图像轨迹统一特征,则联合观测结果为:

    $$ p\left( {\left| y \right., x} \right) = {\theta _{\rm c}}p\left( {{y_{\rm c}}\left| x \right.} \right) + {\theta _{\rm e}}p\left( {{y_{\rm e}}\left| x \right.} \right) $$ (7)

    式中:$ p\left( {\left| {{y_{\rm c}}} \right., x} \right) $、$ p\left( {\left| {{y_{\rm e}}} \right., x} \right) $分别表示两个特征向量;θcθe则是其对应的特征权重值,两者的和为1。

    设特征数量用M表示,第s个特征的观测结果为:

    $$ {\delta _s} = \frac{1}{N}\sum\limits_{i = 1}^N {\left| {p_s^i - \overline {{p_s}} } \right|} $$ (8)

    式中:$ \overline {{p_s}} = \frac{1}{N}\sum\limits_{i = 1}^N {p_s^i} $;δs表示观测峰结果,属于第s个特征,该值越大表示似然程度越高。第s个特征的权重计算公式为:

    $$ {\theta _s} = \frac{{{\delta _s}}}{{\left| {x_s^{\rm peak} - {{\overline x }^{\rm peak}}} \right|}} $$ (9)

    式中:xspeak表示空间位置,属于δs;$ {\overline x ^{\rm peak}} = \frac{1}{M}\sum\limits_{s = 1}^M {x_s^{\rm peak}} $。δs的值与红外图像的真实状态接近程度越高,则公式(9)的分母越小。

    巡检轨迹三维点云配准包含两个部分,一是初始配准,二是精配准。

    变电站机器人巡视轨迹三维点云配准过程中,随机抽取上文融合后的特征,任意对比抽取的特征和轨迹地图特征中的相似程度。

    在点云数据库中,以采集的轨迹图像中,任意一个确定的特征为例,虽然其特征点具备特征值相等的特征,但是,两者之间相互对应的仅会有一个,只有正确地提取融合后特征的关键点[12],才能够得出最佳的目标位姿参数。初步三维点云配准包含离线和在线两个部分,其步骤如下所述:

    离线部分:

    1)提取特征融合后红外图像模板点云M的关键点集pk

    2)求解每个关键轨迹点的融合特征描述子$ FM{1_{p_i^k}} $,其中,i=1, 2, …, k

    3)求解轨迹空间位置描述子$ FM{2_{p_i^k}} $,属于所有关键点。

    在线部分:

    1)分割融合后红外轨迹图像的点云子集O,且O={Oj, j=0, 1, …, n},并获取目标点云关键点。

    2)求解每一个点云子集的两种描述子,分别用F1oF2o表示,前者对应融合特征,即$ F{1_o} = \left\{ {F{1_{{o_j}, }}j = 1, 2, \cdots , m} \right\} $,后者对应空间位置,即$ F{2_o} = \left\{ {F{2_{{o_j}, }}j = 1, 2, \cdots , m} \right\} $。

    3)在模板点云中,搜索与qi相对应的点,且属于Oj中;如果采用pi表示该对应点,需符合下述条件:

    $$ {{\left| {FM{1_{{p_i}}} - F{1_{{o_{{q_i}}}}}} \right|} \mathord{\left/ {\vphantom {{\left| {FM{1_{{p_i}}} - F{1_{{o_{{q_i}}}}}} \right|} {\left[ {FM{1_{{p_i}}} + F{1_{{o_{{q_i}}}}}} \right]}}} \right. } {\left[ {FM{1_{{p_i}}} + F{1_{{o_{{q_i}}}}}} \right]}} < {\varepsilon _1} $$ (10)
    $$ {{\left| {FM{2_{{p_i}}} - F{2_{{o_{{q_i}}}}}} \right|} \mathord{\left/ {\vphantom {{\left| {FM{2_{{p_i}}} - F{2_{{o_{{q_i}}}}}} \right|} {\left[ {FM{2_{{p_i}}} + F{2_{{o_{{q_i}}}}}} \right]}}} \right. } {\left[ {FM{2_{{p_i}}} + F{2_{{o_{{q_i}}}}}} \right]}} < {\varepsilon _2} $$ (11)

    式中:ε1ε2均表示最小平方误差。依据上述两个公式对空间位置描述子实行判断,如果其不在条件的范围内,则将其删除,反之则保留。

    4)RT均表示矩阵,前者对应旋转,后者对应平移,为实现两者的求解,得出初始配准结果[13],通过奇异值分解算法完成,其公式为:

    $$ {\rm Corr} = \left\{ {h_i^1, h_i^2\left| {h_i^1 \in M, h_i^2 \in {O_j}} \right.} \right\} $$ (12)

    式中:hi表示矩阵元素。

    为保证获取更佳的初步配准的结果,需保证RT的最小化,也就是旋转和平移最小,这样定位精度更高。其公式为:

    $$ f\left( {\boldsymbol{R}, \boldsymbol{T}} \right) = {\sum\limits_{i = 1}^n {\left\| {{L_i} - \left( {\boldsymbol{R}{S_i} + \boldsymbol{T}} \right)} \right\|} ^2} $$ (13)

    式中:n表示数量,属于轨迹匹配点对;LiSi均表示点集,前者对应初始,后者对应目标。

    变电站机器人巡视轨迹三维点云的精配准,采用优化的迭代最近点算法完成,精配准时为获取变换矩阵,需完成误差优化,通过最小二乘实现。该算法完成一次迭代后,可提升匹配精准程度[14],属于两种点云之间,即初始和目标。

    通过该方式可减小M和其匹配轨迹点云L之间的偏移角,使下一次实行匹配点搜索的范围减小,可降低误匹配几率。在初始配准的基础上,保证精配准过程中,每一次迭代之后点对之间的法向量夹角θ为减小的依据是动态角迭代因子;并且,通过动态调整阈值,能够实现迭代次数的调整,提升巡视轨迹位姿的更好估计[15]

    构建误差方程,属于匹配点,该构建以θ的正弦值为依据,其公式为:

    $$ \sin \theta = \sqrt {1 - {{\left( {\frac{{\overrightarrow {{n_i}} \times \overrightarrow {{n_j}} }}{{\overrightarrow {\left| {{n_i}} \right|} \times \overrightarrow {\left| {{n_j}} \right|} }}} \right)}^2}} $$ (14)

    式中:ninj均表示近似法向量,前者属于Li,后者属于Si;法向量夹角用θ表示,属于两个匹配点对之间。匹配点的权重计算公式为:

    $$ {Z}_{t}^{i}\left(j\right)=\left\{\begin{array}{cc}1, & \mathrm{sin}\theta < E\\ 0, & 其他\end{array} \right.$$ (15)

    式中:迭代阈值用E表示,对应动态角度,匹配点的数量求解公式为:

    $$ {n_t} = \sum\limits_{j = 0}^N {\sum\limits_i {Z_t^i\left( j \right)} } $$ (16)

    精配准详细步骤如下所述:

    1)获取初始配准得出的初始对应点集,Si0Li0

    2)获取精配准点集Si1Li1,其通过动态角度阈值方法完成,且Si1={$ {s_1}, \cdots , {s_{{n_1}}} $},Li1={$ {l_1}, \cdots , {l_{{n_1}}} $}。

    3)计算R和平移向量t,两者均属于Si1Li1之间。

    4)求解Si1在一次迭代之后的数据点集Si2,且依据下述公式完成:

    $${{S}_{i}}_{2}={{R}_{1}}{{S}_{i}}_{1}{{t}_{i}}$$ (17)

    5)一次迭代结束后,通过公式(18)对E实行缩小:

    $$E=D_{\mathrm{dec}}, D_{\mathrm{dec}} \in(0,1)$$ (18)

    6)循环步骤2)~5)当满足下述条件时则停止:

    $$ \left\{ {\begin{array}{l} {{d_t} - {d_{t + 1}} < \varepsilon \mathop {}\nolimits^{} \mathop {}\nolimits^{} \mathop {}\nolimits^{} \mathop {}\nolimits^{} } \\ {{d_t} = \frac{1}{n}\sum\nolimits_{i = 1}^n {\left\| {{L_{it}} - {S_{it}}} \right\|} } \end{array}} \right. $$ (19)

    式中:dt表示欧氏距离。

    为测试红外图像特征融合下变电站机器人巡视轨迹三维点云配准方法的应用效果,以某变电站为测试对象,对该变电站采用机器人巡检,选用的机器人在满电情况下能够行驶65 km,巡检的变电站平面结构如图 2所示。

    图  2  变电站结构
    Figure  2.  Substation structure

    为测试文本方法红外图像特征融合效果,随机获取变电站特征融合后的红外图像,并将该图像与融合前的图像作对比,衡量红外图像特征融合下变电站机器人巡视轨迹三维点云配准方法的融合效果,结果如图 3所示。

    图  3  图像特征融合效果
    Figure  3.  Image feature fusion effect

    图 3可知:经过特征融合后的图像能够更清晰地呈现巡检画面,细节没有发生影响,并且图像边缘清晰,因此,三维点云配准方法具备良好的特征融合效果。

    为进一步测试红外图像特征融合下变电站机器人巡视轨迹三维点云配准方法的特征融合效果,随机获取变电不同数量的红外图像实行特征融合,采用相似度作为融合效果评价指标,其中相似度包含3个指标即平均梯度、联合熵、偏差指数,3个指标值的范围均在0~1之间,前两者指标值越大,表示融合效果越佳,期望标准在0.86以上;后者的值越小表示表示融合热效果越佳,期望标准在0.21以下。三者的计算公式依次分别为:

    $$ \Delta \bar G = \frac{1}{{MN}}\sum\limits_{m = 1}^M {\sum\limits_{n = 1}^N {\sqrt {\frac{{\Delta {F_x}{{\left( {m, n} \right)}^2} + \Delta {F_y}{{\left( {m, n} \right)}^2}}}{2}} } } $$ (20)
    $$ U{E_{A, F}} = - \sum\limits_{f = 0}^{L - 1} {\sum\limits_{b = 0}^{L - 1} {{p_{A, F}}} } \left( {f, b} \right){\log _2}{p_{A, F}}\left( {f, b} \right) $$ (21)
    $$ \begin{array}{l} {\rm Diff}_{A, F} = \frac{1}{{MN}}\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {\left| {F\left( {i, j} \right) - A\left( {i, j} \right)} \right|} } \hfill \\ {\rm Diff}_{B, F} = \frac{1}{{MN}}\sum\limits_{i = 1}^M {\sum\limits_{j = 1}^N {\left| {F\left( {i, j} \right) - B\left( {i, j} \right)} \right|} } \hfill \\ {\rm Diff} = {\rm Diff}_{F, A} + {\rm Diff}_{F, B} \hfill \\ \end{array} $$ (22)

    式中:F表示特征融合后图像;AB表示源图像;M×N表示图像大小。

    依据上述公式,统计三维点云配准方法在不同图像数量下的指标测试结果,如图 4所示。

    图  4  融合效果评价结果
    Figure  4.  Evaluation results of fusion effect

    图 4可知:随着图像数量的增加,3个指标发生差异性波动曲线,其中平均梯度、联合熵值均在0.86以上;偏差指数则均在0.2以下,因此,三维点云配准方法的特征融合效果良好,融合后的红外图像质量较高。

    在采用文本方法实行巡视轨迹三维配准之前,需确定ε1ε2的取值,以此保证配准的准确性。以配准误差为衡量标准,获取三维点云配准方法在不同取值下的误差结果,如图 5所示。

    图  5  配准误差测试结果
    Figure  5.  Test results of registration error

    图 5可知:随着ε1ε2的取值的逐渐增加,变电站机器人巡视轨迹三维点云配准误差发生差异性变化,配准误差随着ε1取值逐渐增加而发生先下降后上升的变化趋势,初始取值ε1的配准误差为0.7;当取值为1.0,ε1的配准误差为0.72;其中,在取值为0.3时误差值最小;变电站机器人巡视轨迹三维点云配准误差随着ε2取值的增加而发生不规则的变化波动,ε2的配准误差为0.5,ε2的配准误差为0.56;其中,其取值也为0.3时误差最小,因此,确定ε1ε2的取值均为0.3。

    为衡量三维点云配准方法的三维点云配准性能,采用该方法在不同图像大小下实行配准,获取在不同法向量夹角θ下,完成关键点配准结果,并将配准结果与点云中包含的实际关键点数量实行对比,衡量三维点云配准方法的配准性能,结果如图 6所示。

    图  6  配准性能测试结果
    Figure  6.  Registration performance test results

    图 6结果可知:不同图像大小情况下,随着θ的逐渐增加,三维点云配准方法配准的关键点数量与实际数量之间的误差较小,因此三维点云配准方法具备良好配准性能,能够准确完成图像中关键点的配准,能够完成不同大小图像的配准。

    为直观体现三维点云配准方法的巡视轨迹三维点云配准效果,采用该方法对图 2实行轨迹配准,获取配准后的巡视轨迹与期望轨迹的吻合程度,结果如图 7所示。

    图  7  巡视轨迹三维配准结果
    Figure  7.  3D registration results of inspection trajectory

    图 7可知:经过三维点云配准方法配准后的巡视轨迹与期望轨迹吻合程度较高,机器人能够按照设定轨迹完成变电站巡视,保证巡视的全面性。

    本文提出红外图像特征融合下变电站机器人巡视轨迹三维点云配准方法,对机器人巡检轨迹实行配准,保证机器人的巡检效率。该方法通过初步配准和精配准结合完成巡视轨迹配准,能够更好地保证机器人按照期望轨迹完成巡检。经测试:三维点云配准方法的特征融合效果良好,融合后能够提升图像的边缘清晰度,保证巡视图像的质量,并且,能够在误差极小的情况下,完成巡视轨迹配准。

  • 图  1   非球面铝反射镜几何参数图

    Figure  1.   Geometric parameter diagram of an aspheric aluminum mirror

    图  2   机床结构示意图

    Figure  2.   Schematic diagram of machine tool structure

    图  3   三轴联动加工示意图

    Figure  3.   Schematic diagram of three axis linkage machining

    图  4   坐标变换原理

    Figure  4.   Principle of coordinate transformation

    图  5   夹具几何模型

    Figure  5.   Geometric model of fixture

    图  6   夹具-工件有限元模型

    Figure  6.   Finite element model of fixture work piece

    图  7   最大变形量随各参数变化趋势

    Figure  7.   The variation trend of maximum deformation with various parameters

    图  8   夹具和工件应变云图

    Figure  8.   Strain nephogram of fixture and work piece

    图  9   φ682mm非球面铝反射镜加工设备

    Figure  9.   φ682 mm aspheric aluminum mirror processing equipment

    图  10   非球面铝反射镜加工精度检测

    Figure  10.   Machining accuracy test of aspheric aluminum mirror

    表  1   材料参数

    Table  1   Material parameters

    Material Density/(kg·m-3) Young modulus/(GPa) Poisson ratio Coefficient of thermal expansion
    Aluminum alloy 2770 71 0.33 23
    Stainless steel 7750 193 0.31 17
    下载: 导出CSV

    表  2   坐标变换结果

    Table  2   Coordinate transformation results

    Coordinate point X-coordinate/mm Z-coordinate/mm Tool rotation angle-β
    A1 310 219.878377
    A1' 308.75 220.71 56.54
    A2 309.9 219.727043
    A2' 308.65 220.55 56.53
    A3 309.8 219.575771
    A3' 308.55 220.40 56.52
    A4 309.7 219.424562
    A4' 308.45 220.25 56.51
    A5 309.6 219.273416
    A5' 308.35 220.10 56.50
    下载: 导出CSV

    表  3   夹具参数

    Table  3   Fixture parameters

    Factors
    Code
    Level
    1 2 3
    Support rod number N A 12 18 24
    Support rod diameter R/(mm) B 12 22 26
    Plate thickness T/(mm) C 15 20 25
    Spindle speed n/(r/min) 200
    下载: 导出CSV

    表  4   正交实验方案及结果

    Table  4   Orthogonal experimental scheme and results

    NO.
    Support rod
    number
    Support rod
    diameter
    Plate
    thickness
    Max total
    deformation
    N R/mm T/mm D/mm
    1 12 12 15 0.49531
    2 12 22 20 0.31924
    3 12 26 25 0.28334
    4 18 12 20 0.42477
    5 18 22 25 0.30412
    6 18 26 15 0.35454
    7 24 12 25 0.39917
    8 24 22 15 0.41238
    9 24 26 20 0.36372
    下载: 导出CSV

    表  5   表面粗糙度的方差分析

    Table  5   Analysis of variance of surface roughness

    Variation
    source
    SS df MS F Contribution/%
    A 0.002 2 0.001 1.198 0.7
    B 0.020 2 0.010 14.935 52.6
    C 0.013 2 0.006 9.372 31.6
    Error 0.001 2 0.001 15.1
    Sum 0.036 8 0.018 100
    下载: 导出CSV
  • [1] 陈达任, 李占杰, 靳刚, 等. 金属反射镜超精密加工研究进展[J]. 机械研究与应用, 2018, 31(1): 159-163. https://www.cnki.com.cn/Article/CJFDTOTAL-JXYJ201801053.htm

    CHEN Daren, LI Zhanjie, JIN Gang, et al. Review on Ultra- Precision machining of metal reflector[J]. Mechanical Research & Application, 2018, 31(1): 159-163. https://www.cnki.com.cn/Article/CJFDTOTAL-JXYJ201801053.htm

    [2] 谢启明, 杨静, 徐放, 等. 金属非球面反射镜的加工和检测技术[J]. 红外技术, 2015, 37(2): 119-123. http://hwjs.nvir.cn/article/id/hwjs201502007

    XIE Qiming, YANG Jing, XU Fang, et al. Manufacturing and test technology for metal aspherical reflector[J]. Infrared Technology, 2015, 37(2): 119-123. http://hwjs.nvir.cn/article/id/hwjs201502007

    [3] 李全超, 李蕾, 谭淞年, 等. 大口径铝合金主反射镜设计与分析[J]. 应用光学, 2016, 37(3): 337-341. https://www.cnki.com.cn/Article/CJFDTOTAL-YYGX201603003.htm

    LI Quanchao, LI Lei, TAN Songnian, et al. Design and analysis for large aperture primary aluminum mirrors[J]. Journal of Applied Optics, 2016, 37(3): 337-341. https://www.cnki.com.cn/Article/CJFDTOTAL-YYGX201603003.htm

    [4]

    Horst R T, Tromp N, MD Haan, et al. Directly polished lightweight aluminum mirror[C]//Advanced Optical & Mechanical Technologies in Telescopes & Instrumentation. International Society for Optics and Photonics, 2008: 701808-701808-10.

    [5] 姜伟, 申德新, 王立波, 等. 大口径光学反射镜化学镀镍工艺及质量管理[J]. 电镀与涂饰, 2008(10): 13-16. https://www.cnki.com.cn/Article/CJFDTOTAL-DDTL200810008.htm

    JIANG Wei, SHEN Dexin, WANG Libo, et al. Process and quality management of electroless nickel–phosphorus alloy plating on large-aperture optical reflection mirror[J]. Electroplating & Finishing, 2008(10): 13-16. https://www.cnki.com.cn/Article/CJFDTOTAL-DDTL200810008.htm

    [6] 孟晓辉, 王永刚, 李文卿, 等. 应用旋转法实现大口径非球面反射镜零重力面形加工[J]. 光学精密工程, 2019, 27(12): 2517-2524. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201912003.htm

    MENG Xiaohui, WANG Yonggang, LI Wenqing, et al. Fabrication of zero-gravity surface for large-aperture aspherical mirror by using rotationally method[J]. Optics and Precision Engineering, 2019, 27(12): 2517-2524. https://www.cnki.com.cn/Article/CJFDTOTAL-GXJM201912003.htm

    [7] 孙熠璇, 罗世魁, 高超, 等. 大口径空间反射镜的重力卸载优化方法[J]. 红外与激光工程, 2021, 50(1): 279-287. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202101031.htm

    SUN Yixuan, LUO Shikui, GAO Chao, et al. Optimization method for large-aperture space mirror's gravity unload[J]. Infrared and Laser Engineering, 2021, 50(1): 279-287. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202101031.htm

    [8] 马文静, 徐振源, 曹庭分, 等. 大口径反射镜低应力夹持优化设计[J]. 中国激光, 2020, 47(11): 187-192. https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ202011021.htm

    MA Wenjing, XU Zhenyaun, CAO Tingfen, et al. Optimal design of low-stress mounting for large aperture mirror[J]. Chinese Journal of Lasers, 2020, 47(11): 187-192. https://www.cnki.com.cn/Article/CJFDTOTAL-JJZZ202011021.htm

    [9] 肖建国, 张万清, 刘尧, 等. 金属反射镜的三轴联动加工方法[P]. 中国, 20171460886.8, 2020.

    XIAO Jianguo, ZHANG Wanqing, LIU Yao, et al. Three-axis linkage processing method of metal reflector[P]. China, 20171460886.8, 2020.

  • 期刊类型引用(2)

    1. 邢志坤. 基于LabVIEW的变电站移动机器人轨迹跟踪虚拟仿真系统设计. 自动化与仪表. 2024(07): 67-71 . 百度学术
    2. 李辉,余大成,陈耀. 基于OWA算子和CWAA算子的变电站巡视周期优化. 广西电力. 2024(05): 50-54 . 百度学术

    其他类型引用(1)

图(10)  /  表(5)
计量
  • 文章访问数:  157
  • HTML全文浏览量:  51
  • PDF下载量:  31
  • 被引次数: 3
出版历程
  • 收稿日期:  2021-04-11
  • 修回日期:  2021-04-20
  • 刊出日期:  2021-09-19

目录

/

返回文章
返回