CN103761353A - 基于三维边界提取的模具气压砂轮进动加工轨迹优化方法 - Google Patents

基于三维边界提取的模具气压砂轮进动加工轨迹优化方法 Download PDF

Info

Publication number
CN103761353A
CN103761353A CN201310571504.4A CN201310571504A CN103761353A CN 103761353 A CN103761353 A CN 103761353A CN 201310571504 A CN201310571504 A CN 201310571504A CN 103761353 A CN103761353 A CN 103761353A
Authority
CN
China
Prior art keywords
array
normal vector
boundary line
mould
grinding wheel
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
Application number
CN201310571504.4A
Other languages
English (en)
Other versions
CN103761353B (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.)
Langxi Pinxu Technology Development Co ltd
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201310571504.4A priority Critical patent/CN103761353B/zh
Publication of CN103761353A publication Critical patent/CN103761353A/zh
Application granted granted Critical
Publication of CN103761353B publication Critical patent/CN103761353B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Processing Or Creating Images (AREA)
  • Numerical Control (AREA)

Abstract

基于三维边界提取的模具气压砂轮进动加工轨迹优化方法，所述轨迹优化方法包括以下步骤：（1）建立模型，保存轨迹规划的CLSF文件和模型的VRML文件；（2）读取VRML文件，采用正则表达式识别文本数据，然后将数据存储到数组中；（3）用户设定“法向量突变阀值”，其取值依据三维模型上的最大曲率值，曲率越大，阀值应取越大；（4）将单位法向量数组内的法向量进行两两作差后取向量模，同阀值进行比较，根据比较结果删除对角线和重复的边界线；（5）将剩余线段进行整理，得到边界线数组；（6）计算轨迹规划上任意加工点的接触面中心与每条边界的距离，对边界线与加工点的相对位置进行姿态优化。该方法考虑气压砂轮进动加工轨迹规划中的边界问题，效果显著。

Description

基于三维边界提取的模具气压砂轮进动加工轨迹优化方法
技术领域
本发明涉及计算机图形领域，是一种基于VRML技术的三维模型边界线提取方法。 
背景技术
目前，对于三维模型边界线提取的研究主要有以下两方面：(1)从二位图像角度出发，通过一些图像处理方法提取边界线，提取结果是另一张图片，这种方法显然不可能包含每条边界线在三维中的解析式；(2)从实物出发，通过扫描或测量得到点云，再经过搜索算法得到属于边界线的点，最后在计算机上逆向建模。第二种方法可以获得边界线在三维中的解析式，能处理密集的点云数据，一般用于不适合人工建模的具有极度不规则形状的对象，且需要额外的三维扫描测量设备，成本较高，操作难度大，计算过程复杂。三维建模软件虽然能显示线框图，但是无法直接将线框图导出为可用于再开发的数据格式，这给应用三维模型边界线数据进行信息再处理造成不便。 
对于任意的复杂三维模型，三维软件在导出VRML文件时，会将所有表面分割成不重复也不遗漏的一系列三角面，各三角面的位置和相互关系可由组成三角面的顶点及顶点外法向计算得到。因此，本发明提出一种边界线提取方法，用于从VRML文件中提取数据，通过一定的算法获得对应三维模型的边界线数据，很有必要。 
发明内容
为了解决无法从三维建模软件中直接获得模型边界线数据的问题，本发明从三维模型出发，提供一种基于VRML技术的三维模型边界线提取方法，该方法适用范围广，可处理任意复杂形貌的三维模型。 
本发明解决其技术问题所采用的技术方案是： 
一种基于VRML技术的三维模型边界线提取方法，所述边界线提取方法包括以下步骤： 
(1)读取VRML文件，文件扩展名为wrl，并采用正则表达式识别文本数据，然后将数据存储到数组中，数组包括单位法向量数组、线段数组和坐标数组； 
(2)用户设定“法向量突变阀值T”，简称“阀值T”，一般取阀值T=0.2，阀值T的取值依据三维模型上最大曲率值，曲率越大，阀值T应取越大； 
(3)将单位法向量数组内每一条法向量进行两两作差，同阀值T进行比较，若差向量的模小于等于阀值T，则认为两个法向量所对应的三角面是连续的面，面内不存在边界线，并且，如果这两个三角面内存在相同顶点的线段，则一定为面内不存在的边界线，称为“对角线”，删除这两条对角线；如果不存在，则表示两个三角面不属于同一面，不作处理；若差向量的模大于阀值T，则认为这两个三角面是不连续的面，并且，如果这两个三角面内存在相同顶点的线段，则一定为边界线，删除其中一条；如果不存在，则表示两个三角面不相交，不作处理； 
(4)将剩余线段进行整理，得到线段数组，并且和坐标数组一起写入数据表格文件，文件扩展名为csv，保存到计算机磁盘中。 
进一步，所述的基于VRML技术的三维模型边界线提取方法，其特征在于：所述步骤(1)中，所述的单位法向量数组是由三角面的三个顶点的法向量求和再取单位向量得到，所述的线段数组是由三角面三顶点两两相连得到。 
再进一步，所述的基于VRML技术的三维模型边界线提取方法，其特征在于：所述步骤(2)中，所述阀值T的物理意义为：当法向量之差的模小于阀值T时，法向量夹角较小，面夹角接近0度，曲面连续性好。 
本发明的技术构思为：本发明涉及的三维模型边界线提取方法，可将VRML格式文本中的图形数据通过识别和计算，得到模型所有边界线的坐标信息。该方法适用于具有复杂形貌的任意三维模型。 
首先使用三维建模软件，如PRO/E、UG、SolidWorks、CATIA等，建立三维模型；然后将模型另存为VRML(扩展名为wrl)文件格式，选择输出版本为VRML97；然后利用本发明提出的三维模型边界线提取方法，从wrl格式文件中提取出三维模型的边界线数据。 
边界线数据由两个数据表格组成，其中一个表格存储了所有的端点坐标值，每行三个元素表示端点空间直角坐标系中的坐标；另一个表格存储了每一条边界线的两个端点坐标的索引，每行两个元素表示一条边界线的两个端点的索引编号。由这两个数组可以确定模型中的所有边界线。 
本发明的有益效果在于：该提取方法无需额外的设备或测量工具，数据量小，辅助数据充分，计算过程简单可靠，易编程，操作难度低，能精确地得到三维模型边界线数据；只要调整阀值即可识别具有不同曲率的曲面的三维模型；合理设置阀值，可以将大区率的圆角特征识别为边界线，而将小曲率的曲面识别为无边界的连续的面。 
附图说明
图1是本发明提取边界线与优化轨迹的整体流程图 
图2是对角线和重复线段去除流程图 
图3是法向量突变阀值的物理含义示意图 
图4是六棱柱处理前所有线段示意图 
图5是六棱柱处理后所有线段示意图 
图6弧形抛光面、抛光路径和边界线 
图7无约束条件切削速度方向分布(俯视) 
图8有约束条件切削速度方向分布(俯视) 
具体实施方式
下面结合附图对本发明作进一步描述。 
参照图1至图8，并以六棱柱提取边界线为例，提供一种基于VRML技术的三维模型边界线提取方法。所述边界线提取方法包括以下步骤： 
第一步，读取VRML文件，文件扩展名为wrl，文件所描述的三维图形信息主要由以下四部分组成：坐标点（Coordinate point）、点法向量（Normal vector）、三角面顶点序列（CoordIndex）、点法向序列（NormalIndex）。六棱柱模型对应的VRML文件主要内容如下： 
Figure BDA0000415553600000041
Figure BDA0000415553600000051
坐标点记录了模型上所有的端点和曲面经过多面体近似后得到的点的坐标值，按出现顺序指定编号0,1,2,3…，组成“点数组”，如表1所示。 
表1点数组 
Figure BDA0000415553600000052
点法向量记录的是每个点在原三维模型中的外法向三坐标值，且相同的外法向向量只记录一次，并按出现顺序指定编号0,1,2,3…，组成“点法向量数组”，如表2所示。 
表2点法向量数组 
Figure BDA0000415553600000053
三角面顶点序列记录了组成整个三维图形的所有三角形小平面对应顶点编号。每个三角面由三个顶点的编号组成，这里称为“三角面”，可表示为i-j-k形式，-1表示一个面的结束；同理，线段可表示为m-n形式。例如，三角面1-3-7表示由1、3、7三个点，即包含1-3、3-7、7-1三条线段。三角面三个顶点两两相连得到“线段数组”，如表3所示。 
表3线段数组 
Figure BDA0000415553600000061
法向序列和三角面序列一一对应，法向序列的一组包含三个对应三角面三个顶点的点法向量编号。将同一个三角面的三个法向量求和后取其单位向量，作为该三角面的“面单位法向量”，如公式
Figure BDA0000415553600000062
其中，
Figure BDA0000415553600000063
表示三角面5-9-4的面单位法向量，
Figure BDA0000415553600000064
分别为顶点4、5、9的点法向量。求出所有三角面的面单位法向量，组成面“面单位法向量数组”，如表4所示。 
表4面单位法向量数组 
Figure BDA0000415553600000065
第二步，用户设定“法向量突变阀值T”，简称“阀值T”，其物理含义为：当面i-j-k21的面单位法向量Ni-j-k22与面l-m-n26的面单位法向量Nl-m-n25之差P24的模小于阀值T时，法向量夹角θ23较小，夹角α27接近180度，曲面连续性好。所以T的取值依据为工件上最大曲率值，曲率越大，T值应取越大，一般取0.2。 
第三步，参照图4，其中存在很多不必要的线段，如：面4-5-9-8由三角面5-9-4和4-9-8组成，面8-9-11-10由三角面9-11-8和8-11-10组成，其中线段4-9、9-4、11-8、8-11属于不真实存在的对角线，而线段9-8和8-9表示的是同一条边界线，需去除其中一条。 
去除对角线和重复边界线的方法如下：(1)将面单位法向量两两作差，取差向量的模：
Figure BDA0000415553600000066
差向量的模的大小代表了三角面夹角大小。(2)将差 向量的模与所述阀值进行比较，若差向量的模小于等于阀值，则两个法向量所对应的三角面是连续的面，面内不存在边界线，并且，如果这两个三角面内存在相同顶点的线段，则一定为面内不存在的边界线，删除这两条线段；如果不存在，则表示两个三角面不属于同一面，不作处理。若差向量的模大于阀值，则这两个三角面是不连续的面，并且，如果这两个三角面内存在相同顶点的线段，则一定为边界线，删除其中一条；如果不存在，则表示两个三角面不相交，不作处理。 
在六棱柱里具体实施如下： 
表3和表4分别列出了三角面的三个线段组和面单位法向量，然后将任意两个面的面单位法向量作差后取差向量的模与阀值比较。在P≤T时是去除对角线的过程，对于该六棱柱模型，由于其没有曲面，可取T=0。编号为0和1这两个面，这两个三角面为1-5-0和0-5-4，根据表4的面单位法向量，求得其差向量的模为
Figure BDA0000415553600000071
将模P0与阀值比较，P0≤T，并且，通过表3看出，编号0和1的线段数组存在相同线段0-5，所以面内存在不存在的边界线，为面内对角线，删除这两条线段0-5和5-0。六棱柱有20个三角面，要除所有对角线，两两比较需要计算
Figure BDA0000415553600000072
次，计算结果列于表5。 
表5对角线去除后线段 
Figure BDA0000415553600000073
当不存在对角线的情况下，可以直接通过表5所有线段两两比较，若顶点相同，则删除其中一条重复线段，也可以通过所述阀值比较的方法删除两个不连续面的重复线段。如三角面1-5-0和三角面2-1-0这两个面，根据表4的面单位法向量，求得其差向量模为
Figure BDA0000415553600000074
将模P1与阀值T进行比较，当P1≥T，可以确定这两个面为不连续面，并且，这两个三角面内存在相同顶点线 段0-1和1-0，则它一定为边界线段，删除其中一条0-1。对于P1≥T的情况，为了删除所有相同顶点的线段，对去除对角线后剩余的36条线段进行两两比较，若顶点相同，则删除其中一条，需要计算
Figure BDA0000415553600000083
次，得到去除重复线段的结果，如下表6。 
第四步，为方便读取，将表6结果线段进行整理，得到线段数组，如表7所示，整理方法为：对该线段数组内容进行拆分，如原来的线段数组是通过字符串格式“0-5”的方式存储，排列成一列，拆分后变成数字格式“0”、“5”的方式存储。最后将表7的数组和坐标数组一起写入数据表格文件，文件扩展名为csv，保存到计算机磁盘中。 
表6结果线段 
表7结果线段数组 
Figure BDA0000415553600000082
本发明的技术构思为：本发明涉及的三维模型边界线提取方法，可将VRML格式文本中的图形数据通过识别和计算，得到模型所有边界线的坐标信息。该方法适用于具有复杂形貌的任意三维模型。 
首先使用三维建模软件，如PRO/E、UG、SolidWorks、CATIA等，建立三维模型；然后将模型另存为VRML(扩展名为wrl)文件格式，选择输出版本为VRML97；然后利用本发明提出的三维模型边界线提取方法，从wrl格式文件中提取出三维模型的边界线数据。 
边界线数据由两个数据表格组成，其中一个表格存储了所有的端点坐标值，每行三个元素表示端点空间直角坐标系中的坐标；另一个表格存储了每一条边界线的两个端点坐标的索引，每行两个元素表示一条边界线的两个端点的索引编号。由这两个数组可以确定模型中的所有边界线。 
本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围的不应当被视为仅限于实施例所陈述的具体形式，本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。 

Claims (4)

1.基于三维边界提取的模具气压砂轮进动加工轨迹优化方法，其特征在于：所述轨迹优化方法包括以下步骤：
(1)通过三维软件建立模具模型，对模型加工表面进行轨迹规划，以CLSF文件格式保存，文件扩展名为cls；同时，将模型以VRWL文件格式保存，文件扩展名为wrl；
(2)读取所述VRML文件，采用正则表达式识别文本数据，然后将数据存储到数组中，数组包括单位法向量数组、线段数组和坐标数组；
(3)用户设定“法向量突变阀值T”，简称“阀值T”，一般取0.2，阀值T的取值依据三维模型上最大曲率值，曲率越大，阈值T应取越大；
(4)将单位法向量数组内每一条法向量进行两两作差，同阀值T进行比较，若差向量的模小于等于阀值T，则认为两个法向量所对应的三角面是连续的面，面内不存在边界线，并且，如果这两个三角面内存在相同顶点的线段，则一定为面内不存在的边界线，称为“对角线”，删除这两条对角线；如果不存在，则表示两个三角面不属于同一面，不作处理；若差向量的模大于阀值T，则认为这两个三角面是不连续的面，并且，如果这两个三角面内存在相同顶点的线段，则一定为边界线，删除其中一条；如果不存在，则表示两个三角面不相交，不作处理；
(5)将剩余线段进行整理，得到线段数组；
(6)提取所述线段数组，即边界线数据，将其用于气压砂轮进动抛光轨迹规划过程中，对轨迹上的任意一个加工点，计算其接触面中心与每条边界的距离，若其距离值小于设定的距离R，则对边界线与加工点的相对位置进行姿态优化。
2.如权利要求1所述的基于三维模型边界线提取的模具气压砂轮进动加工轨迹优化方法，其特征在于：所述步骤(2)中，所述的单位法向量数组是由三角面的三个顶点的法向量求和再取单位向量得到，所述的线段数组是由三角面三顶点两两相连得到。
3.如权利要求1所述的基于三维模型边界线提取的模具气压砂轮进动加工轨迹优化方法，其特征在于：所述的步骤(3)中，所述的阈值T的物理意义为：当法向量之差的模小于阈值T时，法向量夹角较小，面夹角接近0度，曲面连续性好。
4.如权利要求1所述的基于三维模型边界线提取的模具气压砂轮进动加工轨迹优化方法，其特征在于：所述的步骤(6)中，姿态优化的具体实施方式为：添加一个约束条件，即气压砂轮的切削速度方向不能含有沿边界外法向相反方向的分量。
CN201310571504.4A 2013-11-13 2013-11-13 基于三维边界提取的模具气压砂轮进动加工轨迹优化方法 Active CN103761353B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310571504.4A CN103761353B (zh) 2013-11-13 2013-11-13 基于三维边界提取的模具气压砂轮进动加工轨迹优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310571504.4A CN103761353B (zh) 2013-11-13 2013-11-13 基于三维边界提取的模具气压砂轮进动加工轨迹优化方法

Publications (2)

Publication Number Publication Date
CN103761353A true CN103761353A (zh) 2014-04-30
CN103761353B CN103761353B (zh) 2017-04-26

Family

ID=50528589

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310571504.4A Active CN103761353B (zh) 2013-11-13 2013-11-13 基于三维边界提取的模具气压砂轮进动加工轨迹优化方法

Country Status (1)

Country Link
CN (1) CN103761353B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104741994A (zh) * 2015-03-25 2015-07-01 华南理工大学 一种任意曲面砂轮用于曲面精密磨削的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6167159A (en) * 1998-04-30 2000-12-26 Virtue Ltd. Triangle mesh compression
US20030187625A1 (en) * 2002-03-28 2003-10-02 Fujitsu Limited CAD analysis result data processing apparatus, analytic simulation apparatus, CAD analysis result data processing method and CAD analysis result data processing program
CN101706830A (zh) * 2009-11-12 2010-05-12 中国人民解放军国防科学技术大学 对刚性材质物体表面网格模型进行钻孔后模型重建的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6167159A (en) * 1998-04-30 2000-12-26 Virtue Ltd. Triangle mesh compression
US20030187625A1 (en) * 2002-03-28 2003-10-02 Fujitsu Limited CAD analysis result data processing apparatus, analytic simulation apparatus, CAD analysis result data processing method and CAD analysis result data processing program
CN101706830A (zh) * 2009-11-12 2010-05-12 中国人民解放军国防科学技术大学 对刚性材质物体表面网格模型进行钻孔后模型重建的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
唐世伟 等: "一种复杂二值图像边界跟踪与提取算法", 《一种复杂二值图像边界跟踪与提取算法》 *
张献颖 等: "空间三角网格曲线的边界提取方法", 《中国图象图形学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104741994A (zh) * 2015-03-25 2015-07-01 华南理工大学 一种任意曲面砂轮用于曲面精密磨削的方法
CN104741994B (zh) * 2015-03-25 2017-04-19 华南理工大学 一种任意曲面砂轮用于曲面精密磨削的方法

Also Published As

Publication number Publication date
CN103761353B (zh) 2017-04-26

Similar Documents

Publication Publication Date Title
Giblin et al. A formal classification of 3D medial axis points and their local geometry
WO2024077812A1 (zh) 基于点云语义分割与结构拟合的单体建筑三维重建方法
Lopes et al. Improving the robustness and accuracy of the marching cubes algorithm for isosurfacing
WO2019157924A1 (zh) 三维物体的实时检测方法及系统
Xu et al. Unsupervised trajectory clustering via adaptive multi-kernel-based shrinkage
Lien et al. Approximate convex decomposition of polyhedra
US10719549B2 (en) Querying a database based on a parametric view function
CN104200523A (zh) 一种融合附加信息的大场景三维重建方法
CN103729885A (zh) 多视角投影与三维注册联合的手绘场景三维建模方法
CN110176079B (zh) 一种基于准共形映射的三维模型变形算法
Zaharia et al. 3D versus 2D/3D shape descriptors: A comparative study
Olver et al. Affine invariant detection: edge maps, anisotropic diffusion, and active contours
CN113961738A (zh) 一种多特征铸件三维模型检索方法及装置
CN103631993B (zh) 基于vrml技术的三维模型边界线提取方法
CN103761353A (zh) 基于三维边界提取的模具气压砂轮进动加工轨迹优化方法
CN111275805B (zh) 一种基于3d模型纹理的模型网格优化方法和装置
Wang et al. Layoutmp3d: Layout annotation of matterport3d
Ji et al. Point cloud segmentation for complex microsurfaces based on feature line fitting
Setty et al. Region of interest-based 3D inpainting of cultural heritage artifacts
Liu et al. An approach to 3D building model retrieval based on topology structure and view feature
Kang et al. Fast object recognition using dynamic programming from combination of salient line groups
Wang et al. Machined sharp edge restoration for triangle mesh workpiece models derived from grid-based machining simulation
Sintunata et al. Skewness map: estimating object orientation for high speed 3D object retrieval system
Olver et al. Affine invariant edge maps and active contours
Meenakshisundaram Theory and practice of sampling and reconstruction for manifolds with boundaries

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
TR01 Transfer of patent right

Effective date of registration: 20191204

Address after: 313000 No.5-8, Changhe Road, Huaxi street, Changxing County, Huzhou City, Zhejiang Province

Patentee after: Zhejiang Changxing green battery technology Co.,Ltd.

Address before: 313000 Room 1403, 14th Floor, Building B, Freeport, Headquarters 1188 District Fulu Road, Wuxing District, Huzhou City, Zhejiang Province

Patentee before: Zhejiang creation Intellectual Property Service Co.,Ltd.

Effective date of registration: 20191204

Address after: 313000 Room 1403, 14th Floor, Building B, Freeport, Headquarters 1188 District Fulu Road, Wuxing District, Huzhou City, Zhejiang Province

Patentee after: Zhejiang creation Intellectual Property Service Co.,Ltd.

Address before: Hangzhou City, Zhejiang province 310014 City Zhaohui District Six

Patentee before: Zhejiang University of Technology

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201208

Address after: Room 501, Office Building of Market Supervision Bureau of Langchuan Avenue, Jianping Town, Langxi County, Xuancheng City, Anhui Province

Patentee after: Langxi pinxu Technology Development Co.,Ltd.

Address before: 313000 No.5-8, Changhe Road, Huaxi street, Changxing County, Huzhou City, Zhejiang Province

Patentee before: Zhejiang Changxing green battery technology Co.,Ltd.

TR01 Transfer of patent right