CN112687010B - 一种端框钻模的数字化计量方法 - Google Patents

一种端框钻模的数字化计量方法 Download PDF

Info

Publication number
CN112687010B
CN112687010B CN202110022282.5A CN202110022282A CN112687010B CN 112687010 B CN112687010 B CN 112687010B CN 202110022282 A CN202110022282 A CN 202110022282A CN 112687010 B CN112687010 B CN 112687010B
Authority
CN
China
Prior art keywords
point cloud
jig
cloud data
end frame
drilling jig
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110022282.5A
Other languages
English (en)
Other versions
CN112687010A (zh
Inventor
禹静
刘爽
蔡晋辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Jiliang University
Original Assignee
China Jiliang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Jiliang University filed Critical China Jiliang University
Priority to CN202110022282.5A priority Critical patent/CN112687010B/zh
Publication of CN112687010A publication Critical patent/CN112687010A/zh
Application granted granted Critical
Publication of CN112687010B publication Critical patent/CN112687010B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Numerical Control (AREA)

Abstract

本发明公开了一种端框钻模的数字化计量方法。其包括以下步骤:基于端框钻模的设计模型提取虚拟点云数据,并对虚拟点云数据进行点云处理构建数字化标准钻模;基于端框钻模实体获取实测点云数据,并对实测点云数据进行点云处理构建虚拟模型;基于步骤S1获取的数字化标准钻模和步骤S2的虚拟模型进行相应模型处理,并对二者模型进行尺寸参数计算,点云偏差分析,并输出检测结果报告,实现端框钻模的数字化计量。本发明创造性提出了一种新型的计量方法,有限解决传统量值传递加工成本高,传递路线长,检定过程耗时耗力等弊端。本发明可以快速高效获取钻模所需检测的参数,输出钻模质量检测结果,完成计量检定过程。

Description

一种端框钻模的数字化计量方法
技术领域
本发明涉及计量检定领域,具体的展示了一种端框钻模的数字化计量方法。
背景技术
宇航用端框钻模是用来引导火箭导弹等舱体对接的夹具,根据端框钻模量值传递的需要,分别设计并制造计量过程中的标准钻模、协调钻模和工作钻模。
目前宇航端框钻模的量值传递采用的是实物标准逐级传递的方式。从能复现装配单位的产品要求开始,通过标准钻模、协调钻模逐级传递,最后传递给工作钻模。根据目前端框钻模的传递流程以及钻模本身的特点,发现传统的计量保障模式的矛盾愈发突出,存在的问题主要表现在:检定器具繁多,实用性不强;各类标准器数量繁多,结构复杂,在传递到工作器具时,传递路线长,累积误差大,造成精度损失较大,追溯检查误差来源也存在困难;现有标准器种类繁多,加工成本高,对于通用的标准器建立了相应的量传和溯源体系,但是对于大量的专用标准器,缺乏有效的量传体系。而要建立全面的量传体系,工作量巨大费时。
现有端框钻模传统的量值传递方式不适于如今的工业生产发展,所以急需建立一种新型的数字化量值传递方法来满足端框钻模计量工作的需要,以适于航空航天领域的计量保障。
发明内容
为解决上述技术问题,本发明提供了一种端框钻模的数字化计量方法。所述方法包括如下步骤:
S1:基于端框钻模的设计模型提取虚拟点云数据,并对虚拟点云数据进行点云处理构建数字化标准钻模;
S2:基于端框钻模实体获取实测点云数据,并对实测点云数据进行点云处理构建虚拟模型;
S3:基于步骤S1获取的数字化标准钻模和步骤S2的虚拟模型进行相应模型处理,并对二者模型进行尺寸参数计算,点云偏差分析,并输出检测结果报告,实现端框钻模的数字化计量。
为优化上述技术方案,采取的具体措施还包括:
上述步骤S1中的提取设计模型的虚拟点云数据,包括以下步骤:
S11:将设计模型生成以三角形面片或者多边形面片数据形式的STL、PLY或OBJ格式文件;
S12:基于S11步骤生成的文件,提取三角形面片或者多边形面片数据的顶点;
S13:基于S12步骤生成的点云数据,在各个点云之间进行迭代均匀增密;
S14:将三角面片顶点数据和迭代均匀增密数据保存,完成虚拟点云数据提取。
上述步骤S2中的点云数据处理,包括以下步骤:
S21:搜索与采样点P距离最近的k个近邻点,构建kd-tree来划分空间中的点云数据,建立点云数据间拓扑几何关系;
S22:计算各点的单位法向量、高斯曲率、平均曲率、以及全局平均曲率;
S23:基于K近邻域统计分析方法去除点云数据中的离群点;
S24:基于双边滤波方法平滑点云数据中的小噪声点;
S25:基于体素化下采样方法对大量点云数据进行精简;
S26:基于Delaunay三角化方法对处理好的点云数据进行三维重建,获取数字化标准模型和虚拟模型。
上述步骤S3中的数字化标准钻模和虚拟模型处理,包括以下步骤:
S31:将数字化标准钻模和设计模型进行初始拟合对齐,然后进行最佳拟合对齐点,生成三维彩色点云偏差云图;
S32:对数字化标准钻模和进行内环特征提取,重复10次,获取内环半径参数。
S33:基于S1步骤的拟合对齐后,再进行3D比较,然后均匀随机设置6个比较点,获取各点点云偏差情况,输出数字化标准钻模质量检测报告。
S34:将虚拟模型和数字化标准钻模进行模型处理,重复上述步骤S31~S33,输出虚拟模型质量检测报告,完成端框钻模的数字化计量。
附图说明
图1是本发明方法的流程图。
图2是本发明设计模型提取点云数据的流程图。
图3是本发明设计模型提取点云数据的结果图。
图4是本发明点云去躁的流程图。
图5是本发明点云数据精简后的结果图。
图6是本发明点云数据曲面重构后的结果图。
图7是本发明数字化标准钻模和设计模型拟合对齐后结果图。
图8是本发明三维彩色点云偏差色谱图。
具体实施方式
下面结合附图和实施例对本发明的发明内容作进一步地说明。
参阅图1所示,本发明实施例提供了一种端框钻模的数字化计量方法,其具体包括以下实现步骤:
步骤一、基于端框钻模的设计模型提取虚拟点云数据,并对虚拟点云数据进行点云处理构建数字化标准钻模。本实施例中,参阅图2和图3所示,选择PCL(point cloudlibrary)库来提取设计模型的虚拟点云数据,包括以下步骤:
1)利用PCL点云库的指针读取STL文件;
2)并将读取的模型数据传递给该指针;
3)用VTK(Visualization Toolkit)视觉库转换函数将多边形数据转换成三角网格数据;
4)然后根据需要的点云数量在每个三角网格中根据公式(1)~(3)进行迭代均匀采样;
ai′=ai×M1 (1)
bi′=bi×M2 (2)
式中,(ai,bi,ci)表示三角面片顶点中的任意一点,M1和M2表示为采样随机因子,M2=1-R2,/>rand()是随机数生成函数,生成数的范围是0~PCMAX1,2,PCMAX1、PCMAX2是随机产生一个不小于32767的整数。
5)将生成后的点云坐标保存为PCD格式的点云文件。
步骤二、基于端框钻模实体获取实测点云数据,并对实测点云数据进行点云处理构建虚拟模型;
1)选择K近邻域搜索,搜索与采样点距离最近的K个近邻点。其中,K值为15-30之间较为适宜。
2)构建kd-tree来建立离散点间的拓扑几何关系,其中,本实施例中所获取的点云数据皆是三维数据点,因此k=3。
3)根据公式(6)~(9)计算各个采样点的单位法向量、高斯曲率、平均曲率、以及全局平均曲率。
假设曲面方程如公式(4)所示:
Z(x,y)=ax2+bxy+cy2 (4)
曲面拟合方程取最小值如公式(5)所示:
分别对a、b、c求偏导并取0,可求得常系数a、b、c的值。
点Pi的单位法向量可近似为如公式(6)所示:
点Pi的高斯曲率E如公式(7)所示:
点Pi的平均曲率F如公式(8)所示:
点Pi的全局平均曲率G如公式(9)所示:
式中,Zx、Zy、Zxx、Zxy、Zyy分别是
Zx=2ax+by,Zy=bx+2cy,Zxx=2a,Zxy=b,Zyy=2c;
4)对点云数据中的离群点进行去噪,本实施例中,参阅图4所示,首先根据公式(11)计算采样点Pi的k个近邻点的平均距离Dmean(Pi),然后针对整个点云P,根据公式(12)~(13)估计它的全局距离平均值和标准偏差,平均距离在标准范围之外的点,判定为离群点并从点云数据集中滤除。
根据公式(10),计算点Pi与Pj(xj,yj,zj)(j=1,2,..,k)的欧式距离:
Dj=||Pi-Pj|| (10)
根据公式(11),计算点Pi到它的所有k个近邻点的平均距离Dmean(Pi):
根据公式(12),计算点Pi的均值μ:
式中,N是点云集P的所有数据点。
根据公式(13),计算点Pi的标准偏差σ:
根据公式(14),计算点云集P标准范围Span
Span=μ±g·σ (14)
式中,g=1,2...。
5)对点云数据中的小噪声点进行平滑处理,本实施例中,点云平滑包括以下步骤:
a、基于步骤二的k邻域搜索,计算采样点Pi的k个最近邻域点N(Pi);
b、基于步骤二的的点云法向量估计求三维点云模型中的每个数据点Pi的法失;
c、求出每个邻域点的光顺平滑权函数参数和特征保持权函数参数,分别代表点Pi到邻域点Pj的距离x=||pi-pj||和点Pi的法矢与邻域点Pi的法矢两者的内积y=ni·nj
d、将c的结果代入式和/>式计算出光顺平滑权函数wc(x)和特征保持权函数ws(y);
e、将d的结果代入公式(15),计算出改进双边滤波因子的值α;
式中N(Pi)为数据点Pi的邻域点,ni、nj分别为数据点Pi、Pi的法向量。
f、调整数据点Pi的几何位置,通过(16)式计算使用改进双边滤波后的新数据点pi′的位置,将采样点Pi移动到新的几何位置pi′处;
p′i=pi+α·ni (16)
g、当遍历完所有数据点后,得到光顺去噪后的新点云模型,算法终止。
5)对点云进行数据精简处理,本实施例包括以下步骤:
a、根据公式(17)~(22)确定每个三维体素栅格的边长L、小栅格内包含的点云个数n,各个栅格的一维编码
根据公式(17),计算每个三维体素栅格的边长L:
式中,α为根据点云数据大小调节栅格大小的调解系数因子,s为比例系数,n表示点云数据小栅格内的点云数量。
根据公式(18),计算每个小栅格内包含的点云个数:
n=N/V (18)
式中,N代表这个点云数量的总个数,V代表三维体素栅格长方体的体积,V=LxLyLz,Lx为x轴方向最大长度,Ly为y轴方向最大长度,Lz为z轴方向最大长度。
根据公式(19)~(21),对于任一点Pi(xi,yi,zi),计算其小立方体栅格的编号:
mp=ceil(xi-xmin)/L (19)
np=ceil(yi-ymin)/L (20)
lp=ceil(zi-zmin)/L (21)
式中,ceil(x)为取整函数,不小于x的最大整数。
根据公式(22),将pi的栅格编码为(mp,np,zp),转换为一维编码:
b、将编码放入哈希表中,建立点云数据间的拓扑关系;
c、根据公式(23)~(25)计算体素的重心;若体素内包含点云数据,选择体素重心来替换体素中其他的点,若体素内没有点云数据,则直接删除该体素,遍历完所有体素使原始的点云数据得到了良好的点云精简效果,端框钻模的精简结果参阅图5所示。
6)将处理好的点云进行三维重建,将有向点云投影到某一局部二维平面上,然后再坐标平面内进行平面内的三角化,得到各个点的拓扑关系,再依据Delaunay三角化方法,选择一个三角片样本做出初始曲面,不断的扩张曲面边界得到一个完整的三角网格曲面,最后根据投影点的连接关系确定三维点云数据点的拓扑关系得到最终物体的曲面模型,端框钻模重构模型参阅图6所示。
步骤三、基于步骤一获取的数字化标准钻模和步骤二获取的虚拟模型进行相应模型处理,并对二者模型进行尺寸参数计算,点云偏差分析,并输出检测结果报告,实现端框钻模的数字化计量。本实施中,选择在Geomagic Control环境下进行相应处理,包括以下步骤:
1)状态栏的“图形区域”打开点云处理重建生成的数字化标准模型做为测试模型,再导入设计图形SolidWorks模型作为参考模型。
2)在初始菜单下,选择构造几何圆,采取选择多个点拟合内环,均匀提取三个方向的坐标点用来拟合内环轮廓重复上述操作10次,得到10个拟合内环。
3)在尺寸状态栏下,选取半径尺寸计算各个拟合内环半径。计算10次拟合内环半径的平均值,该值即为数字化标准模型内环半径大小,参阅表1所示,钻模内环半径参考值是300mm,经过数字化标准钻模内环特征提取实测值是300.50279mm。半径误差为0.50279mm,有效满足误差在±1mm的范围。
表1内环半径数据表(mm)
4)参阅图7所示,本实例的钻模工件以基于特征对齐的方式和最佳拟合对齐使数字化标准模型与设计图形对齐。选取定位环作为特征约束来匹配模型,保证对齐精度。
5)参阅图8所示,最后进行3D分析,选取3D比较命令,计算方式设置为默认。确定后可得到3D比较偏差色谱图。
6)在初始状态栏下,设置3D比较点,均匀随机选取了6个测试点,可得到数字化标准模型和设计图形的各点偏差数据,参阅表2所示。
表2测试点偏差数据表(mm)
7)最后生成数字化标准模型质量检测报告。
8)再将虚拟钻模和数字化标准钻模导入到Geomagic Control软件重复上述1)~7)过程,进行三维比对分析处理,获取工件结构尺寸参数,与数字标准钻模参数比对,判断工作钻模产品是否合格,完成钻模数字化计量检定过程。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何限制,凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化,均仍属于本发明技术方案的保护范围内。

Claims (3)

1.一种端框钻模的数字化计量方法,其特征在于,包括以下具体步骤:
S1:基于端框钻模的设计模型提取虚拟点云数据,并对虚拟点云数据进行点云处理构建数字化标准钻模;
S2:基于端框钻模实体获取实测点云数据,并对实测点云数据进行点云处理构建虚拟模型;
S3:基于步骤S1获取的数字化标准钻模和步骤S2的虚拟模型进行相应模型处理,并对二者模型进行尺寸参数计算,点云偏差分析,并输出检测结果报告,实现端框钻模的数字化计量;
所述的步骤S1中的提取设计模型的虚拟点云数据,包括以下步骤:
S11:将设计模型生成以三角形面片或者多边形面片数据形式的STL、PLY或OBJ格式文件;
S12:基于S11步骤生成的文件,提取三角形面片或者多边形面片数据的顶点;
S13:基于S12步骤生成的点云数据,在各个点云之间根据公式(1)~(3)进行迭代均匀增密;
ai′=ai×M1 (1)
bi′=bi×M2 (2)
其中,(ai,bi,ci)表示三角面片顶点中的任意一点,M1和M2表示为采样随机因子,M2=1-R2,/>rand()是随机数生成函数,生成数的范围是0~PCMAX1,2,PCMAX1、PCMAX2是随机产生一个不小于32767的整数;
S14:将三角面片顶点数据和迭代均匀增密数据保存,完成虚拟点云数据提取。
2.根据权利要求1所述的一种端框钻模的数字化计量方法,其特征在于,所述的步骤S1和步骤S2中的点云数据处理,包括以下步骤:
S21:搜索与采样点P距离最近的k个近邻点,构建kd-tree来划分空间中的点云数据,建立点云数据间拓扑几何关系;
S22:计算各点的单位法向量、高斯曲率、平均曲率、以及全局平均曲率;
S23:基于K近邻域统计分析方法去除点云数据中的离群点;
S24:基于双边滤波方法平滑点云数据中的小噪声点;
S25:基于体素化下采样方法对大量点云数据进行精简;
S26:基于Delaunay三角化方法对处理好的点云数据进行三维重建,获取数字化标准模型和虚拟模型。
3.根据权利要求1所述的一种端框钻模的数字化计量方法,其特征在于,所述的步骤S3中的数字化标准钻模和虚拟模型处理,包括以下步骤:
S31:将数字化标准钻模和设计模型进行初始拟合对齐,然后进行最佳拟合对齐点,生成三维彩色点云偏差云图;
S32:对数字化标准钻模进行内环特征提取,重复10次,获取内环半径参数;
S33:基于S31步骤的拟合对齐后,再进行3D比较,然后均匀随机设置6个比较点,获取各点点云偏差情况,输出数字化标准钻模质量检测报告;
S34:将虚拟模型和数字化标准钻模进行模型处理,重复上述步骤S31~S33,输出虚拟模型质量检测报告,完成端框钻模的数字化计量。
CN202110022282.5A 2021-01-08 2021-01-08 一种端框钻模的数字化计量方法 Active CN112687010B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110022282.5A CN112687010B (zh) 2021-01-08 2021-01-08 一种端框钻模的数字化计量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110022282.5A CN112687010B (zh) 2021-01-08 2021-01-08 一种端框钻模的数字化计量方法

Publications (2)

Publication Number Publication Date
CN112687010A CN112687010A (zh) 2021-04-20
CN112687010B true CN112687010B (zh) 2024-03-15

Family

ID=75456393

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110022282.5A Active CN112687010B (zh) 2021-01-08 2021-01-08 一种端框钻模的数字化计量方法

Country Status (1)

Country Link
CN (1) CN112687010B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113538373A (zh) * 2021-07-14 2021-10-22 中国交通信息科技集团有限公司 一种基于三维点云的施工进度自动检测方法
CN113779330B (zh) * 2021-09-13 2024-01-26 河海大学常州校区 一种基于相似性度量的超声检测工艺参数智能计算系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104050712A (zh) * 2013-03-15 2014-09-17 索尼公司 三维模型的建立方法和装置
CN109685778A (zh) * 2018-12-12 2019-04-26 重庆大学 一种基于ct切片的机械零件的几种常见几何量检测方法
WO2019242174A1 (zh) * 2018-06-21 2019-12-26 华南理工大学 基于激光雷达的建筑结构自动测量及3d模型生成方法
CN111047542A (zh) * 2019-12-31 2020-04-21 成都奥伦达科技有限公司 一种电力线补点方法
CN111696210A (zh) * 2020-04-22 2020-09-22 北京航天控制仪器研究所 一种基于三维点云数据特征轻量化的点云重构方法及系统
CN111862171A (zh) * 2020-08-04 2020-10-30 万申(北京)科技有限公司 基于多视图融合的cbct与激光扫描点云数据牙齿配准方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104050712A (zh) * 2013-03-15 2014-09-17 索尼公司 三维模型的建立方法和装置
WO2019242174A1 (zh) * 2018-06-21 2019-12-26 华南理工大学 基于激光雷达的建筑结构自动测量及3d模型生成方法
CN109685778A (zh) * 2018-12-12 2019-04-26 重庆大学 一种基于ct切片的机械零件的几种常见几何量检测方法
CN111047542A (zh) * 2019-12-31 2020-04-21 成都奥伦达科技有限公司 一种电力线补点方法
CN111696210A (zh) * 2020-04-22 2020-09-22 北京航天控制仪器研究所 一种基于三维点云数据特征轻量化的点云重构方法及系统
CN111862171A (zh) * 2020-08-04 2020-10-30 万申(北京)科技有限公司 基于多视图融合的cbct与激光扫描点云数据牙齿配准方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"光栅投影技术及逆向工程在重建和测量牙颌数字化模型中的应用";王珍珍等;《华西口腔医学杂志》;第71-74页 *
"基于逆向工程和快速成型技术的髋骨三维实体模型个性化重建";唐通鸣等;《中国组织工程研究》;第4253-4257页 *

Also Published As

Publication number Publication date
CN112687010A (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
CN113348459B (zh) 将网格几何结构转换为水密边界表示的方法、系统及介质
CN112581457B (zh) 一种基于三维点云的管道内表面检测方法及装置
CN107123164B (zh) 保持锐利特征的三维重建方法及系统
Lee et al. Point data reduction using 3D grids
Woo et al. A new segmentation method for point cloud data
JP5703396B2 (ja) 減数された測定点による公差評価
Lee et al. Data reduction methods for reverse engineering
CN108303037B (zh) 一种基于点云分析的工件表面形状差异的检测方法及装置
KR100810294B1 (ko) 3차원 메쉬 데이터의 특징-유지 간략화 방법
CN112687010B (zh) 一种端框钻模的数字化计量方法
CN110176071B (zh) 一种基于特征模板的三维点云重构方法
CN111145240A (zh) 一种基于3d相机的活体西门塔尔牛体尺在线测量方法
Ravishankar et al. Automated inspection of aircraft parts using a modified ICP algorithm
CN116402866A (zh) 基于点云的零件数字孪生几何建模与误差评定方法及系统
CN113111741A (zh) 一种基于三维特征点的装配状态识别方法
CN113971718A (zh) 一种对三维点云模型进行布尔运算的方法
Nieser et al. Patch layout from feature graphs
CN111563861B (zh) 一种基于三维测量点云数据的工件余量快速求取方法
Wu et al. Bottom-up clustering and merging strategy for irregular curvature aero-blade surface extraction
CN113592976B (zh) 地图数据的处理方法、装置、家用电器和可读存储介质
CN112365534B (zh) 一种基于单目相机三维重建的大型煤堆体积测量方法
CN114970247A (zh) 一种面向叶盘结构的高保真有限元模型自动建模方法
Jiang et al. Lifting wavelet algorithm for freeform surface filtering using a Gaussian prediction operator
CN118070434B (zh) 一种汽车零部件的工艺信息模型构建方法及系统
CN113111548B (zh) 一种基于周角差值的产品三维特征点提取方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant