CN105069312A - 自由曲面测量点分布的规划方法 - Google Patents
自由曲面测量点分布的规划方法 Download PDFInfo
- Publication number
- CN105069312A CN105069312A CN201510523476.8A CN201510523476A CN105069312A CN 105069312 A CN105069312 A CN 105069312A CN 201510523476 A CN201510523476 A CN 201510523476A CN 105069312 A CN105069312 A CN 105069312A
- Authority
- CN
- China
- Prior art keywords
- quadrilateral area
- quadrilateral
- segmentation
- area
- curved surface
- 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.)
- Granted
Links
Landscapes
- Numerical Control (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明涉及一种自由曲面测量点分布的规划方法,属于零件加工制造检测领域。本发明通过设计合理的数值计算方法,并提出了合理的自由曲面局部形状特征的判定参数和判定方法,不但能够达成自适应地根据曲面的局部特性来调整测量点的疏密,实现基于形状特征自适应自由曲面测量点分布规划,而且自由曲面在处理过程中保持一致的四边形拓扑结构,这样有利于实现曲面后续依据测量点数据的重构过程。
Description
技术领域
本发明涉及零件加工制造检测领域,具体涉及一种自由曲面测量点分布的规划方法。
背景技术
自由曲面是指不能由初等解析函数或二次函数进行描述,而必须通过三次或更高次函数进行描述才能得到的曲面。航空、航天等现代装备制造业中涉及到大量的自由曲面零件的检测工作,所涉及到的零件造型复杂、加工精度要求高,对零件的精密检测技术提出了极高的要求。
曲面检测的实质是用一系列离散的采样点提取曲面尽可能准确的实际形状信息,因此曲面测量点的数量及分布将直接决定检测过程的精度和效率,需要运用一定的计算方法进行科学地规划。对于自由曲面进行等间距均匀采样规划是最为简单易行的测量点规划方法,可以通过缩小或增大测量点间距从而在零件测量精度和效率之间取得平衡;但是等间距采样方法冰没有考虑曲面不同位置的形状特征,它在曲面的平顺区域和复杂区域规划相同密度的测量点,这一方面可能使平顺区域规划过多的测量点从而降低曲面数字化的效率,另一方面可能使复杂区域规划过少的测量点从而降低曲面数字化的精度,因此这种测点规划方法不够科学。
实际上,曲面测量点的理想规划方案应当是使测点分布的疏密跟随曲面的形状特征而发生变化,即曲率越大测量点应越密集,反之亦然。然而自由曲面不同位置的形状特征分布非常复杂,曲面的曲率变化没有明显的规律,因此需要设计合理的数值计算方法,实现基于其形状特征自适应测点分布规划。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:如何实现基于形状特征自适应自由曲面测量点分布规划。
(二)技术方案
为了解决上述技术问题,本发明提供了一种自由曲面测量点分布的规划方法,包括以下步骤:
S1、将自由曲面按照u,v方向分解为M×N个四边形区域,并存储每个四边形区域顶点的三维坐标,M、N为正整数,初始化分割次数i为0;
S2、根据所述四边形区域的四个顶点高度坐标和四边形区域的中点高度坐标,确定四边形区域是否需要进一步分解;
S3、根据所述四边形区域的四个顶点高度坐标和四边形区域的中点高度坐标,确定测量点集合。
优选地,步骤S2具体为:
S21、设定分割误差εc作为采样参数,依次读取每一个四边形区域4个顶点的Z向坐标ZV1、ZV2、ZV3、ZV4,并得到四边形中点的Z向坐标ZM,定义Z′M=(ZV1+ZV2+ZV3+ZV4)/4,定义t=|ZM-Z′M|;
S22、将每个四边形区域的t值与εc进行比较,如果全部四边形区域的t值都满足t≤εc,则表明采样条件满足,执行步骤S3;如果有任意四边形区域的t>εc,则将整个自由曲面在该四边形区域进一步分割,得到(2M-1)×(2N-1)个四边形区域,将分割次数i加1,返回步骤S21。
优选地,步骤S3具体为:
S31、预设细分误差εl作为滤波参数,所述滤波参数εl表示区域分割的疏密程度,读取第i次分割后的每个四边形区域4个顶点Z向坐标ZV1 (i)、ZV2 (i)、ZV3 (i)、ZV4 (i)及其中点Z向坐标Z’(i)M,定义Z’(i)M=(ZV1 (i)+ZV2 (i)+ZV3 (i)+ZV4 (i))/4,并得到参数t(i)=|Z(i)M-Z’(i)M|,将t’与εl进行比较,如果存在四边形区域的t(i)≤ε1,则将该四边形区域形状信息过滤,删除该四边形区域本次分割的中点Z向坐标Z(i)M;如果该四边形区域的t(i)>ε1,则该四边形区域形状信息被保留,储存该四边形区域本次分割的中点Z向坐标Z(i)M;直到第i次分割的所有四边形区域都完成判定;然后进行第i-1次分割后各四边形区域的判定,直到完成第1次分割前,M×N多个四边形区域的判定;
S32、输出所有储存的中点Z向坐标信息,形成测量点集合。
(三)有益效果
本发明通过设计合理的数值计算方法,并提出了合理的自由曲面局部形状特征的判定参数和判定方法,不但能够达成自适应地根据曲面的局部特性来调整测量点的疏密,实现基于形状特征自适应自由曲面测量点分布规划,而且自由曲面在处理过程中保持一致的四边形拓扑结构,这样有利于实现曲面后续依据测量点数据的重构过程。
附图说明
图1为本发明实施例的方法流程图;
图2为本发明实施例的方法中判定参数几何含义示意图;
图3为采用B样条及贝塞尔造型方法获得的自由曲面规划得到测点的实例图。
具体实施方式
为使本发明的目的、内容、和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
如图1所示,本发明实施例提供了一种自由曲面测量点分布的规划方法,包括以下步骤:
S1、将自由曲面按照u,v方向分解为M×N个矩形区域(此矩形区域为第1层矩形区域),并存储每个矩形区域顶点的三维坐标,M、N为正整数,初始化分割次数i为0;
S2、根据所述矩形区域的四个顶点高度坐标和矩形区域的中点高度坐标,确定矩形区域是否需要进一步分解,具体步骤包括:
S21、设定分割误差εc作为采样参数,决定了分割区域的大面数量,依次读取每一个矩形区域4个顶点的Z向坐标ZV1、ZV2、ZV3、ZV4,并得到矩形中点的Z向坐标ZM,定义Z′M=(ZV1+ZV2+ZV3+ZV4)/4,定义t=|ZM-Z′M|,Z′M及t的具体几何含义如下图2所示,图2中的“矩形区域”为对图2中的“自由曲面部分”分割得到的对应的矩形,分割的过程将曲面近似为矩形;
S22、将每个矩形区域的t值与εc进行比较,如果全部矩形区域的t值都满足t≤εc,则表明采样条件满足,执行步骤S3;如果有任意矩形区域的t>εc,则将整个自由曲面在该矩形区域进一步分割,得到(2M-1)×(2N-1)个矩形区域,将分割次数i加1,返回步骤S21。至此完成了自由曲面的“采样”阶段。
S3、根据所述矩形区域的四个顶点高度坐标和矩形区域的中点高度坐标,确定测量点集合。步骤S3具体为:
S31、预设细分误差εl作为滤波参数,所述滤波参数εl表示区域分割的疏密程度,读取第i次分割后的每个矩形区域4个顶点Z向坐标ZV1 (i)、ZV2 (i)、ZV3 (i)、ZV4 (i)及其中点Z向坐标Z’(i)M,定义Z’(i)M=(ZV1 (i)+ZV2 (i)+ZV3 (i)+ZV4 (i))/4,并得到参数t(i)=|Z(i)M-Z’(i)M|,将t’与εl进行比较,如果存在矩形区域的t(i)≤ε1,则将该矩形区域形状信息过滤,删除该矩形区域本次分割的中点Z向坐标Z(i)M;如果该矩形区域的t(i)≤ε1,则该矩形区域形状信息被保留,储存该矩形区域本次分割的中点Z向坐标Z(i)M;直到第i次分割的所有矩形区域都完成判定;然后按照上述方式进行第i-1次分割后各矩形区域的判定,再判定第i-2次分割后的各矩形区域的判定,…,按此规律反向地依次判定,直到完成步骤S2第1次分割前,步骤S1分解得到的第1层的M×N多个矩形区域的判定;此过程为曲面的“滤波阶段”。
S32、输出所有储存的中点Z向坐标信息,形成测量点集合。
图3为采用B样条及贝塞尔造型方法获得的自由曲面规划得到测点的实例图。
由以上实施例可以看出,本发明通过设计合理的数值计算方法,并提出了合理的自由曲面局部形状特征的判定参数和判定方法,不但能够达成自适应地根据曲面的局部特性来调整测量点的疏密,实现基于形状特征自适应自由曲面测量点分布规划,而且自由曲面在处理过程中保持一致的四边形拓扑结构,这样有利于实现曲面后续依据测量点数据的重构过程。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (3)
1.一种自由曲面测量点分布的规划方法,其特征在于,包括以下步骤:
S1、将自由曲面按照u,v方向分解为M×N个四边形区域,并存储每个四边形区域顶点的三维坐标,M、N为正整数,初始化分割次数i为0;
S2、根据所述四边形区域的四个顶点高度坐标和四边形区域的中点高度坐标,确定四边形区域是否需要进一步分解;
S3、根据所述四边形区域的四个顶点高度坐标和四边形区域的中点高度坐标,确定测量点集合。
2.如权利要求1所述的方法,其特征在于,步骤S2具体为:
S21、设定分割误差εc作为采样参数,依次读取每一个四边形区域4个顶点的Z向坐标ZV1、ZV2、ZV3、ZV4,并得到四边形中点的Z向坐标ZM,定义Z′M=(ZV1+ZV2+ZV3+ZV4)/4,定义t=|ZM-Z′M|;
S22、将每个四边形区域的t值与εc进行比较,如果全部四边形区域的t值都满足t≤εc,则表明采样条件满足,执行步骤S3;如果有任意四边形区域的t>εc,则将整个自由曲面在该四边形区域进一步分割,得到(2M-1)×(2N-1)个四边形区域,将分割次数i加1,返回步骤S21。
3.如权利要求2所述的方法,其特征在于,步骤S3具体为:
S31、预设细分误差εl作为滤波参数,所述滤波参数εl表示区域分割的疏密程度,读取第i次分割后的每个四边形区域4个顶点Z向坐标及其中点Z向坐标Z’(i)M,定义并得到参数t(i)=|Z(i)M-Z’(i)M|,将t(i)与εl进行比较,如果存在四边形区域的t(i)≤εl,则将该四边形区域形状信息过滤,删除该四边形区域本次分割的中点Z向坐标Z(i)M;如果该四边形区域的t(i)≤εl,则该四边形区域形状信息被保留,储存该四边形区域本次分割的中点Z向坐标Z(i)M;直到第i次分割的所有四边形区域都完成判定;然后进行第i-1次分割后各四边形区域的判定,直到完成第1次分割前,M×N多个四边形区域的判定;
S32、输出所有储存的中点Z向坐标信息,形成测量点集合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510523476.8A CN105069312B (zh) | 2015-08-24 | 2015-08-24 | 自由曲面测量点分布的规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510523476.8A CN105069312B (zh) | 2015-08-24 | 2015-08-24 | 自由曲面测量点分布的规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105069312A true CN105069312A (zh) | 2015-11-18 |
CN105069312B CN105069312B (zh) | 2018-02-09 |
Family
ID=54498677
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510523476.8A Active CN105069312B (zh) | 2015-08-24 | 2015-08-24 | 自由曲面测量点分布的规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105069312B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109387938A (zh) * | 2017-08-14 | 2019-02-26 | 清华大学 | 光学系统自由曲面面形公差分布的分析方法 |
CN112665523A (zh) * | 2020-11-24 | 2021-04-16 | 北京星航机电装备有限公司 | 一种复杂型面的联合测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101968344A (zh) * | 2009-09-22 | 2011-02-09 | 大连海事大学 | 数控加工中心在机三维形面检测系统 |
CN103411574A (zh) * | 2013-08-14 | 2013-11-27 | 西北工业大学 | 航空发动机叶片型面三坐标测量方法 |
CN103486996A (zh) * | 2013-08-14 | 2014-01-01 | 西北工业大学 | 未知cad模型的航空发动机叶片型面测量方法 |
CN103777570A (zh) * | 2014-01-07 | 2014-05-07 | 浙江大学 | 基于nurbs曲面的加工误差快速检测补偿方法 |
-
2015
- 2015-08-24 CN CN201510523476.8A patent/CN105069312B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101968344A (zh) * | 2009-09-22 | 2011-02-09 | 大连海事大学 | 数控加工中心在机三维形面检测系统 |
CN103411574A (zh) * | 2013-08-14 | 2013-11-27 | 西北工业大学 | 航空发动机叶片型面三坐标测量方法 |
CN103486996A (zh) * | 2013-08-14 | 2014-01-01 | 西北工业大学 | 未知cad模型的航空发动机叶片型面测量方法 |
CN103777570A (zh) * | 2014-01-07 | 2014-05-07 | 浙江大学 | 基于nurbs曲面的加工误差快速检测补偿方法 |
Non-Patent Citations (2)
Title |
---|
张海等: "自由曲面测量点分布规划研究", 《机床与液压》 * |
樊喜刚等: "梯度蜂窝面外动态压缩力学行为与吸能特性研究", 《固体力学学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109387938A (zh) * | 2017-08-14 | 2019-02-26 | 清华大学 | 光学系统自由曲面面形公差分布的分析方法 |
CN109387938B (zh) * | 2017-08-14 | 2020-08-11 | 清华大学 | 光学系统自由曲面面形公差分布的分析方法 |
CN112665523A (zh) * | 2020-11-24 | 2021-04-16 | 北京星航机电装备有限公司 | 一种复杂型面的联合测量方法 |
CN112665523B (zh) * | 2020-11-24 | 2022-04-19 | 北京星航机电装备有限公司 | 一种复杂型面的联合测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105069312B (zh) | 2018-02-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110516388B (zh) | 基于调和映射的曲面离散点云模型环切刀轨生成方法 | |
CN109377561A (zh) | 一种基于共形几何的数模表面网格生成方法 | |
CN112784403B (zh) | 基于点云数据建立节理岩体离散元模型的数值模拟方法 | |
CN107392875A (zh) | 一种基于k近邻域划分的点云数据去噪方法 | |
CN111696210A (zh) | 一种基于三维点云数据特征轻量化的点云重构方法及系统 | |
CN103870845A (zh) | 点云聚类去噪过程中新的k值优化方法 | |
CN102306396A (zh) | 一种三维实体模型表面有限元网格自动生成方法 | |
CN103942837B (zh) | 基于序列线性规划的叶片点云模型截面曲线直接构造方法 | |
CN108230452B (zh) | 一种基于纹理合成的模型补洞方法 | |
Lu et al. | Uncertainty-guided intelligent sampling strategy for high-efficiency surface measurement via free-knot B-spline regression modelling | |
CN104616349A (zh) | 基于局部曲面变化因子的散乱点云数据精简处理方法 | |
CN108180918A (zh) | 一种点云测地路径正向跟踪生成方法及装置 | |
CN103106632A (zh) | 一种基于均值漂移的不同精度三维点云数据的融合方法 | |
CN103617336A (zh) | 一种航空噪声等值线图的绘制方法 | |
CN111340723A (zh) | 一种地形自适应的机载LiDAR点云正则化薄板样条插值滤波方法 | |
CN103236043A (zh) | 一种植物器官点云修复方法 | |
CN106815447A (zh) | 基于历史数据的复杂结构件加工特征智能定义与分类方法 | |
CN105069312A (zh) | 自由曲面测量点分布的规划方法 | |
CN105045958A (zh) | 基于bp神经网络的gps高程拟合的实现系统及方法 | |
CN107886573B (zh) | 一种复杂地质条件下边坡三维有限元网格生成方法 | |
CN110310322A (zh) | 一种10微米级高精度器件装配表面检测方法 | |
CN113409332A (zh) | 一种基于三维点云的建筑物平面分割方法 | |
CN106803266A (zh) | 一种船体复杂外板点云肋骨线提取方法及装置 | |
CN116186864B (zh) | 一种基于bim技术的深基坑模型快速建模方法及系统 | |
CN104484900A (zh) | 基于Delaunay三角网的气象要素插值评估方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |