基于点云特征对比的曲面翘曲变形检测方法
Surface warpage detection method based on point cloud feature comparison
通讯作者:
收稿日期: 2020-07-30
Received: 2020-07-30
作者简介 About authors
陈杨波(1997—),男,硕士生,从事缺陷检测的研究.orcid.org/0000-0001-6228-9899.E-mail:
针对目前曲面翘曲变形描述不充分、检测效率低的问题,定义自由曲面的翘曲变形描述方式,提出基于点云特征对比的曲面翘曲变形检测方法. 研究一系列算法,开展曲面单层点云提取,通过对比曲面实测点云与模板点云的空间位置获取翘曲变形区域;通过计算获得翘曲距离和翘曲张角,描述翘曲的变形程度和变化趋势. 实例分析结果表明,提出的方法不需要进行曲面重建,直接通过点云特征对比进行曲面翘曲变形分析,在保证精度的同时,在效率上比三维重建方法检测翘曲有了较大的提高.
关键词:
A description method of free-form surface warpage was defined aiming at the problems of insufficient description and low detection efficiency of surface warpage, and a detection method of surface warpage based on point cloud feature comparison was proposed. A series of algorithms were analyzed. The surface single-layer point cloud was extracted, and the warpage deformation area was obtained by comparing the spatial position of the measured point cloud and the template point cloud of the surface. Then the warpage distance and the warpage angle were calculated to describe the deformation degree and variation trend of the surface warpage. The results of the example analysis show that the proposed method directly analyzes the surface warpage through the comparison of point cloud features without performing surface reconstruction, which ensures the great accuracy and is more efficient than the three-dimensional reconstruction method in detecting the warpage.
Keywords:
本文引用格式
陈杨波, 伊国栋, 张树有.
CHEN Yang-bo, YI Guo-dong, ZHANG Shu-you.
随着产品轻量化要求越来越高,聚合物注塑件由于重量轻、强度高等特点,逐渐替代了部分金属零件,特别是一些薄壁零件,在航空航天、高铁汽车等领域得到广泛应用. 在注塑成型过程中,由于壁厚较薄,聚合物的熔体在模具型腔内的冷却过程难以控制,导致注塑件很容易产生各种缺陷,翘曲变形是其中最复杂的一种缺陷[1].
除了这些针对平整表面的缺陷检测的方法,也有一些针对几何形变缺陷的研究. 吴梅君[13]基于三维点云数据,结合Delaunay三角剖分、曲面拟合、误差分析等方法,实现了对壳体表面的三维重建和外观形变检测;王颖等[14]在对散乱三维点云进行轴向分层和移位规则化处理的基础上,通过寻找相邻点间的几何特征突变进行缺陷检测;Simler等[15]采用光度立体视觉(表面法向量)和立体视觉(三维点云),检测自由形状零件的表面缺陷;Wei等[16]提出快速、鲁棒、自动的缺陷轮廓重构方法,利用随机样本一致性和迭代最近点算法建立缺陷轮廓模型,开展轮廓表面重建. 这些方法针对性较强,不能直接用于翘曲的检测,且基本都需要进行三维重建,在一定程度上制约了检测效率的提升.
本文定义自由曲面的翘曲变形描述方法,提出基于点云特征对比的曲面翘曲变形检测方法. 通过对比曲面实测点云与模板点云的空间位置,检测曲面的翘曲变形。特点在于通过点云直接进行翘曲变形描述指标的计算,不需要进行曲面重建,有助于提高检测的效率.
1. 基于点云的曲面翘曲描述与定义
为了描述曲面翘曲变形的特征,采用翘曲距离和翘曲张角2个指标,对曲面翘曲的变形程度和变化趋势进行描述.
对变形曲面上任意一点的翘曲指标进行描述. 任意一点的翘曲距离用模板曲面的切平面垂直测得的表面翘曲距离表示. 以模板曲面某点处的切平面作为基准面,将该点沿切平面法向量方向投影至变形曲面,投影点与切平面之间的距离即为该投影点的翘曲距离. 任意一点的翘曲张角用变形曲面某点处切平面与对应的模板曲面同一点处切平面之间的夹角来表示. 在上述翘曲距离描述的基础上,变形曲面上投影点的切平面与模板曲面上点的基准切平面之间的夹角为翘曲张角.
任意一点的翘曲距离和翘曲张角描述如图1所示. 取上方模板曲面的某处切平面Sq作为基准面,从切点q处作垂线与下方变形曲面产生交点p,点p到切平面的距离d为该切平面垂直测得的表面翘曲距离,即点p处的翘曲距离;过p点作下方变形曲面的切平面Sp,平面Sp和Sq法向量
图 1
在曲面任意一点翘曲指标描述的基础上,结合曲面翘曲变形的特点,对曲面翘曲变形进行基于点云的描述. 由于薄壁件的曲面翘曲变形一般由固定侧向自由侧单调递增,自由侧边界点一般为翘曲变形较大处。采用若干自由侧边界点的翘曲距离和张角来描述曲面点云的翘曲变形,避免了单个值的偶然性和特殊性,能够得到翘曲变形程度和变形趋势.
典型的曲面翘曲变形的点云描述如图2所示. 不失一般性,对基于点云的曲面翘曲的定义如下.
图 2
1)基于点云的曲面翘曲距离定义.
设变形曲面点云簇上点
从大到小的K个边界点的翘曲距离集合
2)基于点云的翘曲张角的定义.
设上述变形曲面点云簇上的点
上述翘曲距离集合中K个边界点的翘曲张角集合
2. 基于点云对比的翘曲检测关键算法
2.1. 曲面单层点云提取算法
由于点云为物体的表面轮廓,对于薄壁件,提取到的点云一般为形状相似、距离相近的两层,重复表达了物体的特征,而且两层点云距离一般较近,容易产生干扰,不利于曲面翘曲的描述. 采用曲面单层点云提取算法获取双层曲面点云中的一层,具体包括以下步骤.
1)采用表面法线估计算法[17]计算点云P =
2)将点云
式中:
3)遍历移动后点云中的每一点
由于算法针对的是双层薄壁表面,法线估计的邻域半径
2.2. 变形区域的点云簇提取与匹配算法
1)将实测点云、模板点云分别构造为八叉树.
2)遍历提取实测点云中的点作为搜索点,在模板点云中进行体素搜索,即在以搜索点对应相同位置为中心、边长为
3)将模板点云中的点作为搜索点,采用上述相同的方法,在实测点云中进行体素搜索,得到危险区域模板点云.
4)采用欧几里得聚类分割方法[20],对危险区域实测点云进行分割. 将点云中相邻距离小于距离阈值
5)若存在多处变形区域,开展每处变形区域实测点云簇与模板点云簇的对应匹配. 计算实测点云簇与每个模板点云簇之间的距离,选择其中距离最小的模板点云簇作为匹配的点云簇,建立对应关系. 其中点云簇之间的距离计算方法如下:在实测点云簇中随机选择
2.3. 翘曲距离和翘曲张角计算方法
基于上述曲面翘曲变形的点云描述,翘曲距离和翘曲张角的计算可以分为以下步骤.
1)翘曲距离的计算.
采用基于法线估计的边界提取方法[20],提取得到变形区域实测点云簇的边界点. 根据1章的定义,以变形区域实测点云簇上的边界点
根据式(1)计算点
2)翘曲张角的计算.
提取上述翘曲距离集合中对应的K个边界点,根据1章的定义,取这K个边界点中的点
根据式(2)计算点
3. 翘曲变形检测流程
提出的翘曲变形检测方法的流程如图3所示.
图 3
1)点云数据获取. 获取三维点云数据,包括曲面的实测点云数据和模板点云数据.
3)单层点云提取. 采用上述曲面单层点云提取方法,对叶轮实测点云和模板点云分别进行单层点云提取,获取叶轮叶片曲面点云的其中一层.
4)变形区域点云提取. 通过实测点云与模板点云的空间位置对比检测,获取可能出现翘曲的危险区域,提取该危险区域的实测点云及对应的模板点云;采用欧几里得聚类分割方法,分别对危险区域的实测点云和模板点云进行分割,提取得到变形区域的实测点云簇与对应的模板点云簇.
5)变形判断. 根据变形区域实测点云簇的数量进行判断,若得到的点云簇数量为0,则说明该曲面没有翘曲变形;若不为0,则实测点云簇数量即为该曲面出现翘曲变形的数量.
6)对应点云簇匹配. 若有翘曲变形存在,则进行对应点云簇匹配,通过点云簇之间的距离计算方法,得到变形区域实测点云簇与每个模板点云簇之间的距离;选择其中距离最小的模板点云簇作为匹配的点云簇,建立对应关系.
7)翘曲距离和翘曲张角计算. 提取得到变形区域实测点云簇的边界点,计算边界点的翘曲距离,并从大到小选择K个翘曲距离作为翘曲距离集合;选择这K个翘曲距离所对应的边界点,计算这些点的翘曲张角,得到翘曲张角集合. 将这2个集合作为该处曲面翘曲变形的描述.
4. 实例分析
采用典型的薄壁注塑件叶轮作为研究对象,验证所提出的方法. 由于注塑叶轮的叶片厚度一般较小,通常有着复杂的曲面造型,导致在制造过程中很容易出现翘曲变形,尤其是在叶片末端;叶片叶形对叶片的性能有着很大的影响,因此对叶轮叶片翘曲变形的检测显得十分重要.
选取出现翘曲变形的注塑件叶轮为研究对象,注塑叶轮实物图和标准模型图如图4所示.
图 4
图 5
在得到叶轮点云数据后,对原始点云进行预处理. 测量得到叶轮实测点云中两点之间的最小间距,作为之后参数选择的衡量标准. 开展点云滤波和点云配准. 其中,半径滤波的筛选半径根据点云的密度确定,本文选用点云中两点最小距离的10倍,将叶轮点云模型的噪点剔除.
对处理后的叶轮实测点云和模板点云分别进行曲面单层点云提取. 其中,移动距离
图 6
图 7
图 8
将每处翘曲的实测点云簇与对应的模板点云簇进行匹配,其中
图 9
提取变形区域实测点云簇的边界点,计算得到边界点的翘曲距离和翘曲张角,选取其中翘曲距离最大的
表 1 曲面翘曲变形描述结果
Tab.1
翘曲变形区域 | 边界点坐标 | | |
区域1 | (228.4,21.6,61.7) | 2.021 | 8.95 |
区域1 | (227.4,21.7,62.5) | 1.983 | 10.75 |
区域1 | (227.9,21.7,61.8) | 1.978 | 5.06 |
区域1 | (227.5,21.7,61.9) | 1.971 | 5.60 |
区域1 | (244.0,20.7,92.6) | 1.958 | 12.78 |
区域1 | (244.5,20.7,92.4) | 1.957 | 12.83 |
区域1 | (227.5,21.7,61.3) | 1.944 | 3.77 |
区域1 | (229.0,21.6,63.7) | 1.943 | 4.89 |
区域2 | (85.9,17.9,58.1) | 1.628 | 5.43 |
区域2 | (86.8,17.9,57.0) | 1.592 | 3.47 |
区域2 | (86.8,17.9,57.5) | 1.573 | 3.54 |
区域2 | (90.9,17.6,54.0) | 1.570 | 4.35 |
区域2 | (86.8,17.9,57.5) | 1.570 | 3.54 |
区域2 | (87.1,17.9,57.4) | 1.567 | 3.54 |
区域2 | (86.6,17.9,57.9) | 1.562 | 3.53 |
区域2 | (85.3,18.0,58.9) | 1.559 | 3.54 |
图 10
图 11
将本文注塑叶轮用智能三坐标测量仪进行翘曲的测量,得到2处翘曲值分别为2.00和1.60 mm. 将这2个值作为标准值,K个边界点的翘曲距离均值作为各个区域的实测翘曲值,计算得到本文测量方法的相对误差为
上述实例的编译开发在Visual Studio 2017中完成,并借助于PCL点云库,在CPU主频为2.40 GHz、内存为8 GB的电脑设备上运行,测得所需的时间为1 166 ms. 将本文检测方法与三维重建测量方法进行对比. 三维重建测量方法的种类和步骤众多,选取其中最关键的三维表面重建步骤,采用PCL点云库中的贪婪投影三角算法,对叶轮实测点云进行三维表面重建. 在同样配置的运行设备下,测得重建所需时间为12 652 ms. 三维表面重建结果如图12所示.
图 12
由本文方法与三维模型重建方法的对比可知,采用三维模型重建的方法进行翘曲检测,仅三维表面重建这一步骤所需时间就超过本文检测方法所需时间的10倍,三维表面重建后的翘曲检测需要耗费一定的时间. 可以看出,利用本文检测方法显著提高了检测效率. 三维重建一般是需要对点云进行三角化,随着点云数量和密度的增加,重建所需要的时间大大增加. 在提出的检测方法中,大多数步骤只需针对变形区域的点云簇,对检测时间的影响不显著,因此本文检测方法在检测数量与密度大的点云时优势更明显.
5. 结 语
本文定义自由曲面的翘曲变形描述方法,采用翘曲距离和翘曲张角2个指标,分别对曲面翘曲的变形程度和变化趋势进行描述,使得曲面翘曲变形更加直观和全面. 提出基于点云特征对比的曲面翘曲变形检测方法,通过对比曲面实测点云与模板点云的空间位置,开展曲面翘曲变形的检测和描述指标的计算. 该方法的特点在于通过点云直接进行翘曲变形描述指标的计算,不需要进行曲面重建,在保证精度的同时,有助于提高检测的效率,为薄壁注塑件的翘曲变形检测提供了一种便利的方法. 本文方法针对的是自由曲面的翘曲,可以用于电子仪器外壳、汽车覆盖件这类薄壁曲面翘曲变形的描述与检测.
以薄壁注塑件叶轮作为研究对象,进行实例分析. 结果表明,本文检测方法通过翘曲变形程度和变化趋势2个指标,表达了翘曲变形的特征;本文测量精度较好,相对于三维重建方法,在保证精度的同时,检测效率有了较大的提高.
参考文献
塑件翘曲度及其计算方法
[J].DOI:10.3321/j.issn:1000-7555.2008.06.001 [本文引用: 1]
The warpage degree of injection-molded parts and evaluation methods
[J].DOI:10.3321/j.issn:1000-7555.2008.06.001 [本文引用: 1]
Research of three dimensional laser scanning coordinate measuring machine
[J].DOI:10.1051/matecconf/201823202015 [本文引用: 1]
3D shape modeling using a self-developed hand-held 3D laser scanner and an efficient HT-ICP point cloud registration algorithm
[J].
3D point cloud surface reconstruction based on divide-and-conquer method in laser scanner
[J].
3D point cloud analysis for detection and characterization of defects on airplane exterior surface
[J].DOI:10.1007/s10921-017-0453-1 [本文引用: 1]
基于3D扫描技术的木材缺陷定量化分析
[J].DOI:10.11975/j.issn.1002-6819.2017.07.022 [本文引用: 2]
Quantitative analysis of wood defect based on 3D scanning technique
[J].DOI:10.11975/j.issn.1002-6819.2017.07.022 [本文引用: 2]
基于散乱三维点云的缺陷检测和三维重构方法
[J].DOI:10.5768/JAO201637.0303001 [本文引用: 2]
Method for defects detection and 3D reconstruction based on dispersed points cloud
[J].DOI:10.5768/JAO201637.0303001 [本文引用: 2]
Geometric modeling using octree encoding
[J].
A method for registration of 3-D shapes
[J].DOI:10.1109/34.121791 [本文引用: 1]
/
〈 |
|
〉 |
