飞机蒙皮机器人化制造曲面分片算法研究

Civil Aircraft Design & Research

2019 No. 3

Sum No. 134

http： //myfj. cnjoumals. com myfj_sadri@comae, cc (021)20866796

DOI： 10. 19416/j. enki. 1674 -9804.2019.03.002

飞机蒙皮机器人化制造曲面分片算法研究

王乾杰赵兴炜陶波*

(华中科技大学数字制造装备与技术国家重点实验室，武汉430074)

摘要：

为实现飞机曲面机器人加工的自动轨迹规划，研究了基于模型几何特征的飞机曲面分片算法。通过分析飞机构件加工的工 艺参数信息，以表面加工允许的最大曲率或最大偏角为阈值，提出基于近似曲率的三角网格模型分片方法。以三角网格中任 意面片为起点，沿网格邻接关系向外拓展联接，形成若干近似平面化的模型分片结构，当超过阈值时则停止拓展，从而得到一 整片曲面。通过对飞机CAD模型进行分片，发现采用本文提出方法可以保证分片结果的均匀性，为后期加工轨迹自动规划规 划提供保障。

关键词：曲面分片；表面加工；网格模型

中图分类号:V262.2 + 3 文献标识码：A OSID：

〇引言

随着机器人自动化加工的逐渐普及，自动化打 磨与自动化喷涂技术已广泛应用于汽车及家具加工 等工业领域，且相较传统的手动加工,机器人自动化 加工表现出了良好的加工精度和极高的加工效率。 飞机加工领域对比汽车与家电加工行业有着更高的 精度要求,为保障飞机飞行时的安全性，加工厚度与 一致性必须满足设计需要，因此将机器人自动化加 工引人飞机加工领域势在必行。

洛克希德•马丁公司针对F-35战机开发机器 人飞机精整系统,实现对尾翼和小零件的自动喷涂， 涂层厚度变化范围较手工喷涂提升80%[1];国内清 华大学[2]就飞机涂装自动喷涂技术进行解析和系 统设计，喷涂误差在10%范围以内；中国民航大 学[3]针对飞机维护设计的涂层打磨系统，对飞机大 机翼打磨的人力需求从6人降低至2人，极大地提 升了加工效率。

一方面，大型飞机构件可以根据其实际功能分

解为多个功能结构，经人工分割后分别进行加工。

另一方面，大型飞机构件属于大型复杂构件，需分解 为多个功能结构,分别进行加工。传统加工方式采 用手动分片，在规划加工轨迹的过程中便已完成各 功能区域划分[4]，或是根据机械臂可达范围对模型 进行分片处理[>3]。文献[5-7]在规划轨迹前对 工件三角网格模型进行分片处理，对每一个分片区 域分别进行轨迹规划，实现分操作的自动运行，拓展 了机器人自动化加工的涵盖范围，但从实现角度来 看依然存有一些问题：（1)分片结果不能保证任意 两三角面片间法向矢量偏角或近似曲率小于设定阈 值，而这正是保证加工质量的关键所在；（2)分片结 果面积小数量多，极为琐碎，不利于加工效率的提 升。随后文献[8]提出基于分水岭法的曲面分片方 式，寻找并连接曲率发生突变的特征点形成曲面分 片边界;文献[9-10]提出结合区域增长法的曲面 分片方式，连接起点邻接具有相同曲面特征的面片 形成加工分片。文献介绍方法可以一定程度上增大 分片结果的面积，减少加工轨迹规划的次数,同时让

[基金项目]国家重点研发计划（批准号：2017YFB1301504)资助项目。

* 通信作者• E-mail: taobo@hust.edu.cn

引用格式：王乾杰，赵兴炜，陶波.飞机蒙皮机器人化制造曲面分片算法研究[J].民用飞机设计与研究，2019(3):9-14.

WANG Q J, ZHAO X W, TAO B. Investigation on Surface Segmentation Algorithm for Robotic Manufacturing of Aircraft Skin[J]. Civil Aircraft Design and Research，2019(3) :9-14( in Chinese).

9

专栏：第三届COMAC国际科技创新周总第134期

每一个区域具有相同的曲率特征，但依然无法从数 工区域的最大曲率限制为：

字上保证区域中的偏角和曲率小于阈值。文献 [11]采用区域增长法，以三角面片法向矢量偏角为 k ^2-cUl--71 - -52)--

(4)

Lc

依据进行加工区域划分，保证分片区域的最大偏角 本文使用飞机三角网格模型，数据格式为一 小于偏角阈值。但文献针对复杂曲面喷涂加工，仅 系列三角形片体构成，保存三角面片的三个顶点 考虑偏角对加工质量的影响，不便于在打磨加工轨 坐标和法向矢量坐标。在计算近似曲率的过程 迹规划中使用，仍具有提升空间。

中，计算两三角面片的法向矢量偏角0及质心间 本文将根据表面加工的工艺模型对加工时的厚

距L近似认为质心间距为一段曲率圆弧长，偏角 度及力学信息进行分析，提出以三角网格模型近似 为圆弧对应的圆心角。由此可计算得到曲率圆弧 曲率为判据的曲面分片方法，提升分片区域加工性 的半径：

能。该方法可以生产相对比较平滑的曲面，有利于 (5)

后期机器人加工轨迹的生成。并在Visual Studio 2017平台编写算法实现代码，验证大型飞机构件分 由此计算得到近似曲率计算方法：

片算法的效率和可靠性。

=e ^.d{ \\ -s2)K appro -

i —

2

(6)

1分片判据阈值的计算

在近似过程中将两三角面片质心间距视为曲率 大型飞机曲面加工，无论是涂层喷涂或是表 圆的一段圆弧，则作为对比的这两个三角面片必须 面打磨，因其复杂的表面特征，都需进行曲面分 足够接近，则在使用判据过程中，对比的三角面片应 片操作。引言中介绍了几种曲面加工分片方法， 为相邻两三角面片。

但都或多或少的存在不足，需对喷涂加工与打磨 加工模型进行分析，获得贴合加工模型的分片判 断依据。

通过分别分析曲面喷涂加工与曲面打磨加工的 厚度模型与加工力学模型，分析曲面分片划分区域 应满足的曲面特征，从而实现对大型飞机曲面分片 阈值判据的推导与计算。

文献[12]根据赫兹接触理论，推导复杂曲面砂 轮打磨的力学模型，得到截面线的压强分布函数：

nx) --

(i)

图1接触面压强分布

如图1所示，其中P。为打磨区域最大压强j

为压力分布区域求解数据：

文献[13 ]分析喷涂加工的涂层厚度累积模型， 令

分

别

为

最

小

、平均、最大漆膜厚度，

a为任意位置s漆膜厚度,为该位置曲面法向矢

R = y/R • d (3)量与与之相切平面法向矢量之间的偏角，则任意点 其中，ft,、/?2分别为砂轮变形后中心点曲率半 $漆膜厚度应满足式(7):

径和工件接触中心点曲率半径，为压痕深度。

Q^coia.h < Q, < Q„,„ (7)

为保证分片后划分范围内加工质量，加工区域 若令4匕为涂层厚度最大允许偏差，则对喷涂 的压强分布应尽量均匀，至少不能出现砂轮边界脱 区域整体喷涂厚度应满足：

离工件表面的情况。若无脱离情况，则可认为A = Q,-_

(

Qd _ AQd)

(9)

(Z/2)^SP。，其中4为磨削接触宽度，可以求得加

Q/rain

10

2019年第3期王乾杰，等飞机蒙皮机器人化制造曲面分片算法研究

在计算近似曲率阈值的过程中，将质心间距Z 默认为1，使近似曲率数值阈值上与法向矢量偏角相等。

分片操作开始后，以三角网格模型中的随机 三角面片7>,为起点，搜索与之相邻的三角面片

计算两相邻三角面片法向矢量偏角若偏角 7满足曲面分片判据，即M < /W。或是y <〜，则

将三角面片b与三角面片k相连形成组合面片 心，并更新组合面片边界三角面片7>^。随后搜 索与边界三角面片相邻的三角面片，并对相邻两 个三角面片的偏角进行计算与分片判据进行对 比，将满足者连接入组合面片。重复此类操作，直 至组合面片心中所有边界三角面片均不存在满足 分片判断的相邻三角面片，则将组合面片心视为 一个三角面片几何并输出为三角网格模型，模型 所描述便是一个曲面分片结果。随后再以任意零 散三角面片为起点，重复上述操作，完成下一个分 片区域的划分和储存，直至三角网格模型S中所 有的三角面片均不存在满足分片判据的邻接三角 面片，如图2所示。2. 2对分片区域近似平面化

在简单根据第1节介绍的判据对模型进行分 片操作过程中，由于判据针对对象为三角网格模 型中的邻接三角形，分片结果可能不能直接使用 于机械臂自动化加工。在对飞机模型进行分片的 过程中，飞机的机身可近似视为一个圆柱，在曲面 分片过程中可能会被完整提取出来如图3所示， 对比后续加工中机械臂的可达范围，机身曲面的 面积过大且不具备近似平面性的特征，需在后续 进行优化。

2曲面分片及优化原理

2. 1基于区域增长法的曲面分片实现

本文对大型飞机模型进行曲面分片基于区 域增长法，根据后续加工工艺的不同（曲面打磨 或曲面喷涂）分别选取曲面分片判据，对模型进 行处理。

获得飞机三角网格模型数据，将模型数据记录 为一组三角面片集：

S = {.Trt\\ i = 1,2,... (1〇)

其中，7>,表示三角网格模型中的第；个三角面 片，模型中三角面片总数为〜每个三角面片7>,.数 据中包含面片的三个顶点坐标z,}，以及面片 的法向矢量若进行打磨操作的曲面分片，则 加上面片的质心坐标。获取后续加工工艺参

数，根据式(6)和式(9)计算得到曲面分片判断阈值

Kappro^ a,h〇图3未优化分片效果

为使分片结果可近似为平面，需对上文中的分 片结果进行优化处理。将每个代表分片区域的组合

面片定义为一个三角面片集：

Scb = {Tr,. I i = 1,2,... ,n> (11)

定义中各参数含义于式（10)中相同。从三角

面片集&中边界三角面片为起点，构件优化面片

具体操作如下：

11

专栏：第三届COMAC国际科技创新周总第134期

1) 设定近似平面化区域最大偏角九ax，表示在

分片结果优化过程中可以允许的优化结果区域中最 大的法向矢量偏角。

2) 计算优化面片&的面积平均法向矢量no™%:

V Ai • /iorm-n〇rmavg = —~~n--------------

(12)

S A

其中，n为优化面片中成员三角面片的数量，/l, 和norm,.分别为第i个成员三角面片的面积和法向 矢量。初始优化面片中仅有起点三角面片,平均法 矢即三角面片法矢。

3)

搜索与优化面片中边界三角面片邻接的三

角面片，并计算三角面片法向矢量与优化面片平均

矢量的偏角T,若y <0max，将三角面片连接进人优

化面片，并更新优化面片的平均法向矢量。4) 重复上述操作，直至优化面片中边界三角面

片不存在与之相邻且满足角度阈值限制的三角面

片，则优化面片构建完成，保存优化面片集合，完成

—组曲面分片优化结构。

5) 重复上述操作，直至组合面片中所有的三角

面片均加人优化面片，且所有优化面片构建完成，则 曲面区域近似平面化优化完成。

近似平面化处理后的模型表面将所有的法向矢 量偏角限制在了最大偏角以内，符合近似平面 化的思想，更便于机器人自动化加工的实现。

2. 3小型区域与狭长区域的识别与剔除

随着分片的进行，组合面片或优化面片可能出 现包含三角面片数量极少便构造完成的情况。此类 特征区域对应曲面模型中曲率较大的区域,不满足 保证加工质量的曲面几何限制，在实际加工中不适 用于机器人自动加工,或者不适用于作为参数输入 的表面加工工艺及工具。仅对曲面分片算法辅助的 表面加工规划，为保证加工表面厚度均匀性及维持 自动化加工的高效性，不宜对此类特征区域进行加 工，需进行识别与剔除。

完成近似平面化的加工区域优化后，计算每个 加工区域的面积'，并计算飞机模型表面积4,设 置曲面分片存贮面积最小占比若加工区域面积 与飞机模型面积的比值小于I即

'

. 4,则将

该区域从分片结果中剔除，在后续加工中有工人手 工加工或更换加工参数对此区域进行规划。

12

3算法验证

大型飞机曲面分片算法的验证采用式（9)描述 的曲面喷涂加工的角度阈值判据对飞机模型进行曲 面分割,飞机模型下载自网络，如图4所示。

图4

飞机模型

程序代码运行于Windows 10操作系统，Visual

Studio 2017平台下使用C ++语言编写，结合0pen-

CasCade库实现算法显示模块实现。程序的开发和 运行平台：CPU:Inter® Core™ i7 -7700 @ 3.60GHz, 内存:8GB,显卡:3GB。

导人原始模型（* • igs)，使用商用软件 SOLIDWORKS将模型转换为三角网格模型（*. STL),飞机三角网格中包含18016个三角面片，模 型分片前设置法向矢量阈值％为〇. 35mrf,近似平 面化区域最大偏角为lrad，存贮面积最小占比

s设置为10% ,输入飞机三角网格模型，运行曲面 分片算法。

分片算法输出三角网格模型置于同一视图

中，分片效果图如图5所示，从图中可以看出，飞 机曲面分片算法成功将飞机模型划分成为机翼上 下表面（3、5)、尾翼上下表面（4、6),以及被拆分 为若干近似平面的飞机机身（1、2、7)共计7个区 域。所有分片区域均具有近似平面的特性，三角 网格模型中显示的锯齿状边界在曲面拟合后可化 为平滑边界。

U)俯视图

2019年第3期王乾杰，等飞机蒙皮机器人化制造曲面分片算法研究

[3]祝世兴，蔡高存.自适应飞机大翼打磨设备的设计

与分析[J].机械设计与研究，2018, 34(6): 115- 118.

[4 ][5 ]

林焰，衣正尧，李玉平，等.大型船用坞内外板涂装 机器人[J].机器人，2018, 40(1): 115-128.Wei Chen, Dean Zhao. Path planning for spray painting robot of workpiece surfaces [ J ]. Mathematical Problems in Engineering, 2013： 1-6.

0>)仰视图

图5

分片算法结果视图

[6 ]

Heping Chen, Ning Xi, Weihua Sheng. A general framework for automatic CAD-guided tool planning for surface manufacturing [ C ]. Proceedings of the 2003 IEEE International Conference on Robotics and Automa­tion. Taipei： IEEE： 2003.

[7 ]

Mayur V. Andulkar, Shital S. Chiddarwar. Incremental approach for trajectory generation of spray painting robot [J]. Industrial Robot： an International Journal, 2015, 42(3)： 228-241.

[8 ]

韩丽，高小山，楚秉智.离散曲率约束的三角网格模 型拓扑分割算法[J].计算机辅助设计与图形学学 报，2009, 21(6): 831-835.

[9] 刘胜兰，周儒荣，安鲁陵.三角网格模型的数据分块

算法[J].南京航空航天大学学报，2003, 35(6): 653-658.

文件大小4.86Mb

[10]

Zezhong C. Chen, Zuomin Dong, Geoffrey W. Vickers. Automated surface subdivision and tool path generationfor 3 -

-

算法验证过程中的模型数据与程序运行内存与 时间消耗如表1所示，从表中可以看出，分块整体时

间平均在5. 82s左右，对比实际模型加工耗时，时间 占比小，拥有很高的处理效率;曲面分片算法处理后 的曲面模型表面特征均符合喷涂加工的参数需求，便 于加工规划，保证了加工厚度的均匀性;算法简单，实 现了大型飞机表面加工的全面自动化，提高模型生成 效率，从而实现了加工轨迹规划流程整体效率。由此 看来，本算法在加工中使用可带来显著的效果。

表1

解决方案Release

时间消耗5.82s

算法验证参数内存消耗2.60Mb

文件格式ASCII

4结论

分析曲面喷涂和曲面打磨的工艺模型，对曲面 分片算法的阈值判据进行推导与计算，采用区域增 长法对目标模型进行曲面分片，使分片区域拥有一 致的表面特征，保证后续加工在分片区域内的均匀 性;通过对分片结果进行优化处理，使分片曲面具有 近似平面化的特征，符合轨迹规划习惯，且不会过于 琐碎，便于对模型进行加工轨迹规划以提升规划效 率。最后通过实际模型对算法进行验证，确保算法 的可实现性。参考文献：

[1 ]

刘亚威.机器人喷涂在F-35的应用[J]•航空科学 技术，2011(5): 19-18.

[2]缪东晶，吴聊，徐静，等.飞机表面自动喷涂机器人

系统与喷涂作业规划[J].吉林大学学报（工学版）， 2015, 45(2)： 547-553.

[12]

-^--axis CNC machining of sculptured parts

[J]. Computers in Industry, 2003 , 50 ： 319-331.

[11] 冯浩，龙迁瑶，潘海鹏.适用于喷涂的点云模型分块

处理[J]•机械设计与制造，2016(3): 257-260.丁小天.基于三维点云的大型复杂构件机器人自主 加工轨迹规划[D].武汉：华中科技大学，2019: 14- 18.

[13] 刘亚举，面向曲率组合曲面的喷涂机器人喷枪轨

迹优化研究[D].兰州：兰州理工大学，2016: 21-23.

作者简介

王乾杰男，硕士研究生。主要研究方向：机器人自主加工

分片规划。E-mail:381235095@qq.com

赵兴炜男，博士，讲师。主要研究方向：机器人动力学与控制 0 E-mail ：zhaoxingwei@ hust. edu. com

陶波男，博士，教授，博士生导师。主要研究方向：机器人学。E-mail: taobo@ hust. edu. cn

13

专栏：第三届COMAC国际科技创新周总第134期

Investigation on Surface Segmentation Algorithm for

Robotic Manufacturing of Aircraft Skin

WANG Qianjie ZHAO Xingwei TAO Bo *

(State Key Lab of Digital Manufacturing Equipment and Technology, Huazhong University of

Science and Technology, Wuhan 430074, China)

Abstract ： In order to realize the automatic trajectory planning of robotic machining, surface segmentation algorithm

for aircraft based on the geometric characteristics of the machining model was studied. Based on the analysis of the processing parameter, the segmentation algorithm for triangle mesh based on approximate curvature was proposed, and the maximum curvature or maximum deviation angle of the processing was used as the threshold. After that, adjacent triangle meshes were connected together, which started from an arbitrary patch in the mesh and end when the threshold was reached. As a result, a approximated planar segmentation surfaces was formed. The proposed al­gorithm was utilized for surface segmentation of a CAD model of aircraft. Its results showed the uniformity of the segmented surface, which could ensure the smoothness of the generated robotic machining trajectory.Keywords： surface segmentation； surface processing； triangular mesh

* Corresponding author. E-mail： taobo@hust.edu.cn14