CN111710022B - 一种避免轮廓相交的快速成型切片处理方法 - Google Patents
一种避免轮廓相交的快速成型切片处理方法 Download PDFInfo
- Publication number
- CN111710022B CN111710022B CN202010515616.8A CN202010515616A CN111710022B CN 111710022 B CN111710022 B CN 111710022B CN 202010515616 A CN202010515616 A CN 202010515616A CN 111710022 B CN111710022 B CN 111710022B
- Authority
- CN
- China
- Prior art keywords
- node
- layer
- slice
- graph
- area
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/005—General purpose rendering architectures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种避免轮廓相交的快速成型切片处理方法,其步骤如下:步骤一、根据成型精度确定层高和切片总数,并获取各层切片平面与STL模型上每个三角面片的交点;步骤二、将每层切片内的交点集合映射为一个图G,并确定图G的极大连通子图Gi;步骤三、获取每个连通子图Gi中的节点度数,并根据节点度数特性识别轮廓相交情况,创建有序封闭轮廓;步骤四、利用每层切片得到的轮廓建立打印封闭区域,并与相邻切片层的区域进行布尔运算以对打印区域进行修正。本发明利用图的节点度数特性,智能识别切片轮廓相交情况,避免了网格拓扑关系的复杂计算,能够有效去除多余轮廓分支,建立正确的打印区域,适用于各种规则和非规则模型的切片处理。
Description
技术领域
本发明涉及快速成型领域,具体涉及一种避免轮廓相交的快速成型切片处理方法。
背景技术
快速成型,又称为3D打印或增材制造,是一种以数字模型文件为基础,通过逐层累积堆叠材料来构造物体的技术。在打印加工前,首先需沿三维模型的成型方向建立一系列切片,以提取模型的截面轮廓,然后根据截面轮廓生成打印填充区域。目前已有的针对STL模型进行切片处理的算法主要分为两大类:
1)基于点的处理方法.首先计算层切平面与模型上所有三角片的交点,然后对交点进行排序以构成有序封闭轮廓,例如文献1(田仁强,张义飞.快速成型中STL模型直接切片新算法研究[J].机床与液压,2019,47(16):55-59);
2)基于三角片的处理方法.首先求得层切平面与第一个三角片的交点,然后根据拓扑关系找到与之相邻的三角片,依次求交后即可得到首尾相连的有向封闭轮廓,例如文献2(徐敬华,盛红升,张树有,等.基于邻接拓扑的流形网格模型层切多连通域构建方法[J].计算机辅助设计与图形学学报.2018,30(1):180-190)。
无论是传统的基于点的方法还是基于三角片的方法,都是在构建成一条有序封闭轮廓后,进行下一条轮廓的创建;但当同一层切片内的多条轮廓相交时,会造成排序结果混乱。因此,改善分层切片算法的有效性和准确性,对于提高快速成形***的可靠性和制作精度具有重要意义。
发明内容
针对现有技术存在的不足,本发明的目的在于提供一种避免轮廓相交的快速成型切片处理方法。
为实现上述目的,本发明提供了如下技术方案:
一种避免轮廓相交的快速成型切片处理方法,其步骤如下:
步骤一、根据成型精度确定层高和切片总数,并获取各层切片平面与STL模型上每个三角面片的交点;
步骤二、将每层切片内的交点集合映射为一个图G,并确定图G的极大连通子图Gi;
步骤三、获取每个连通子图Gi中的节点度数,根据节点度数特性识别轮廓相交情况,并根据轮廓相交情况采用深度优先搜索算法或三角剖分方法创建有序封闭轮廓;
步骤四、利用每层切片得到的轮廓建立打印封闭区域,并与相邻切片层的区域进行布尔运算以对打印区域进行修正。
所述步骤二中,对于每层切片求得的交点集合,首先进行冗余点去除,获得节点集合V,然后根据每个三角面片上的两个交点构成一条边的特性将节点集合V映射为一个图G=(V,E),其中E为节点集合V中的节点构成的所有边集,对图G进行深度优先遍历,得到图G的极大连通子图Gi,i=1,2,3,...,n,n为连通分量总数。
所述步骤三中节点度数特性包括:①连通子图Gi中所有节点的度数都≤2;②连通子图Gi中存在一个节点v的度数deg(v)>2,其中节点v的度数deg(v)定义为与v关联的所有边的关联数目总和。
当连通子图Gi中所有节点的度数都≤2时,采用深度优先搜索算法创建有序封闭轮廓,首先从子图Gi中的某一节点v0出发,搜索与它关联的邻接节点v1,然后从节点v1出发,搜索与v1邻接且未被访问过的节点v2,依次进行搜索,直到回到起点v0,即创建一条有序封闭轮廓。
当连通子图Gi中存在一个节点v的度数deg(v)>2时,采用三角剖分方法创建有序封闭轮廓。
所述三角剖分方法如下:
Step1,将连通子图Gi内的所有节点利用Delaunay三角剖分方法构建一个三角网格DT;
Step2,搜索三角网格DT的边界集合Boundary,边界边按照如下规则定义:Boundary={<vi,vj>|Dep<vi,vj>=1},其中<vi,vj>代表三角网格DT中由节点vi、vj构成的边,Dep<vi,vj>为依附于边<vi,vj>的三角形的个数,当Dep<vi,vj>=1时,<vi,vj>为边界边,vi、vj为边界节点;
Step3,判断边界集合Boundary中的每一条边<vi,vj>是否属于连通子图Gi的边集Ei,①若存在将与该边<vi,vj>相连的三角形删除,更新三角网格DT后,转Step2;②若不存在直接转Step4;
Step4,依次连接三角网格DT的边界节点即为连通子图Gi内的有序封闭轮廓。
所述步骤四中包括以下步骤:
Stepl,在每层切片内,根据各个封闭轮廓之间的内外包含关系建立多边形封闭区域;
Step2,获取第i层切片的封闭区域Ai与第i+j层切片的封闭区域Ai+j的差集d1j=Ai-Ai+j,j=1,2,...,N,并求出它们的并集:获取第i层切片的封闭区域Ai与第i-j层切片的封闭区域Ai-j的差集d2j=Ai-Ai-j,j=1,2,...,N,并求出它们的并集:
Step3,将区域Si=d1∪d2作为物体第i层切片的表层打印区域,区域Ti=Ai-Si作为物体第i层切片的内部打印区域。
本发明的有益效果:利用图的节点度数特性,智能识别切片轮廓相交情况,避免了网格拓扑关系的复杂计算,能够有效去除多余轮廓分支,建立正确的打印区域,适用于各种规则和非规则模型的切片处理。
附图说明
图1为本发明的技术流程示意图。
图2a、b、c、d为本发明的三角剖分交点排序示意图。
图3a、b、c、d、e、f、g为本发明和现有方法处理切片轮廓相交的实例效果图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明公开了一种避免轮廓相交的快速成型切片处理方法,其步骤如下:
步骤一、根据成型精度确定层高和切片总数,并获取各层切片平面与STL模型上每个三角面片的交点;
步骤二、将每层切片内的交点集合映射为一个图G,并确定图G的极大连通子图Gi;
步骤二中,对于每层切片求得的交点集合,首先进行冗余点去除,获得节点集合V,然后根据每个三角面片上的两个交点构成一条边的特性将节点集合V映射为一个图G=(V,E),其中E为节点集合V中的节点构成的所有边集,对图G进行深度优先遍历,得到图G的极大连通子图Gi,i=1,2,3,...,n,n为连通分量总数。
步骤三、获取每个连通子图Gi中的节点度数,根据节点度数特性识别轮廓相交情况,并根据轮廓相交情况采用深度优先搜索算法或三角剖分方法创建有序封闭轮廓;
步骤三中节点度数特性包括:①连通子图Gi中所有节点的度数都≤2;②连通子图Gi中存在一个节点v的度数deg(v)>2,其中节点v的度数deg(v)定义为与v关联的所有边的数目总和。
当连通子图Gi中所有节点的度数都≤2,表示该连通子图内不存在轮廓相交情况。采用深度优先搜索算法创建有序封闭轮廓,首先从子图Gi中的某一节点v0出发,搜索与它关联的邻接节点v1,然后从节点v1出发,搜索与v1邻接且未被访问过的节点v2,依次进行搜索,直到回到起点v0,即创建一条有序封闭轮廓。
当连通子图Gi中存在一个节点v的度数deg(v)>2,表示该连通子图内存在轮廓相交情况。由于快速成型会只在轮廓构成的封闭区域内部进行材料填充,因此对于此种类型的节点拓扑排序,只需找到该连通子图的外边界节点,因此采用三角剖分方法创建有序封闭轮廓,具体为:
Step1,将连通子图Gi内的所有节点利用Delaunay三角剖分方法构建一个三角网格DT;
Step2,搜索三角网格DT的边界集合Boundary,边界边按照如下规则定义:Boundary={<vi,vj>|Dep<vi,vj>=1},其中<vi,vj>代表三角网格DT中由节点vi、vj构成的边,Dep<vi,vj>为依附于边<vi,vj>的三角形的个数,当Dep<vi,vj>=1时,<vi,vj>为边界边,vi、vj为边界节点;
Step3,判断边界集合Boundary中的每一条边<vi,vj>是否属于连通子图Gi的边集Ei,①若存在将与该边<vi,vj>相连的三角形删除,更新三角网格DT后,转Step2;②若不存在直接转Step4;
Step4,依次连接三角网格DT的边界节点即为连通子图Gi内的有序封闭轮廓。
例如,图2a所示的连通图Gi=(Vi,Ei),其节点集合Vi={1,2,...,9},边集Ei={<1,2>,<2,3>,<3,4>,<4,5>,<5,9>,<9,8>,<9,6>,<6,8>,<8,7>,<7,1>},其中节点8和9的度数为3。首先将节点集合Vi进行三角剖分得到三角网格DT,如图2b所示,然后搜索网格DT的边界集合F={<1,2>,<2,3>,<3,4>,<4,5>,<5,7>,<7,1>},由于集合F中的边界这因此将与边<5,7>相连的三角形△578删除,得到新的三角网格DT后(图2c),重新搜索其边界集合F,发现边界边因此删除与之相连的△589,此时新的三角网格DT的边界集合F完全属于图Gi中的边集Ei(图2d),故依次连接边界节点,得到有序封闭轮廓为1→2→3→4→5→9→8→7→1。
步骤四、利用每层切片得到的轮廓建立打印封闭区域,并与相邻切片层的区域进行布尔运算以对打印区域进行修正。具体如下:
Step1,在每层切片内,根据各个封闭轮廓之间的内外包含关系建立多边形封闭区域;
Step2,获取第i层切片的封闭区域Ai与第i+j层切片的封闭区域Ai+j的差集d1j=Ai-Ai+j,j=1,2,...,N,并求出它们的并集:获取第i层切片的封闭区域Ai与第i-j层切片的封闭区域Ai-j的差集d2j=Ai-Ai-j,j=1,2,...,N,并求出它们的并集:
Step3,将区域Si=d1∪d2作为物体第i层切片的表层打印区域,区域Ti=Ai-Si作为物体第i层切片的内部打印区域。
实施例:
如图3a所示的曲柄零件模型,在切片P内产生的交点如图3b所示,若采用传统方法会产生如图3c所示的交点排序错乱。本发明首先将切片P内的交点集合映射为图,并利用Delaunay三角剖分方法构建三角网格,如图3d所示;然后利用图的边集对剖分网格进行裁剪,得到最终的三角网格如图3e所示;顺序连接网格的边界节点即得到有序封闭轮廓,如图3f所示;将该封闭轮廓构建成封闭区域,并和相邻切片层的区域进行布尔运算,得到该层切片的有效打印区域如图3g所示,其中直线填充区域为该零件在切片P的表层打印区域,网格填充区域为该零件在切片P的内部打印区域。
由此,本发明利用图的节点度数特性,智能识别切片轮廓相交情况,避免了网格拓扑关系的复杂计算,能够有效去除多余轮廓分支,建立正确的打印区域,适用于各种规则和非规则模型的切片处理。
实施例不应视为对本发明的限制,但任何基于本发明的精神所作的改进,都应在本发明的保护范围之内。
Claims (4)
1.一种避免轮廓相交的快速成型切片处理方法,其特征在于,其步骤如下:
步骤一、根据成型精度确定层高和切片总数,并获取各层切片平面与STL模型上每个三角面片的交点;
步骤二、将每层切片内的交点集合映射为一个图G,并确定图G的极大连通子图Gi;
步骤三、获取每个连通子图Gi中的节点度数,根据节点度数特性识别轮廓相交情况,并根据轮廓相交情况采用深度优先搜索算法或三角剖分方法创建有序封闭轮廓;
步骤四、利用每层切片得到的轮廓建立打印封闭区域,并与相邻切片层的区域进行布尔运算以对打印区域进行修正,
所述步骤三中节点度数特性包括:①连通子图Gi中所有节点的度数都≤2;②连通子图Gi中存在一个节点v的度数deg(v)>2,其中节点v的度数deg(v)定义为与v关联的所有边的数目总和,当连通子图Gi中存在一个节点v的度数deg(v)>2时,采用三角剖分方法创建有序封闭轮廓,
所述三角剖分方法如下:
Step1,将连通子图Gi内的所有节点利用Delaunay三角剖分方法构建一个三角网格DT;
Step2,搜索三角网格DT的边界集合Boundary,边界边按照如下规则定义:Boundary={<vi,vj>|Dep<vi,vj>=1},其中<vi,vj>代表三角网格DT中由节点vi、vj构成的边,Dep<vi,vj>为依附于边<vi,vj>的三角形的个数,当Dep<vi,vj>=1时,<vi,vj>为边界边,vi、vj为边界节点;
Step3,判断边界集合Boundary中的每一条边<vi,vj>是否属于连通子图Gi的边集Ei,①若存在将与该边<vi,vj>相连的三角形删除,更新三角网格DT后,转Step2;②若不存在直接转Step4;
Step4,依次连接三角网格DT的边界节点即为连通子图Gi内的有序封闭轮廓。
2.根据权利要求1所述的一种避免轮廓相交的快速成型切片处理方法,其特征在于,步骤二中,对于每层切片求得的交点集合,首先进行冗余点去除,获得节点集合V,然后根据每个三角面片上的两个交点构成一条边的特性将节点集合V映射为一个图G=(V,E),其中E为节点集合V中的节点构成的所有边集,对图G进行深度优先遍历,得到图G的极大连通子图Gi,i=1,2,3,...,n,n为连通分量总数。
3.根据权利要求1所述的一种避免轮廓相交的快速成型切片处理方法,其特征在于,当连通子图Gi中所有节点的度数都≤2时,采用深度优先搜索算法创建有序封闭轮廓,首先从子图Gi中的某一节点v0出发,搜索与它关联的邻接节点v1,然后从节点v1出发,搜索与v1邻接且未被访问过的节点v2,依次进行搜索,直到回到起点v0,即创建一条有序封闭轮廓。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010515616.8A CN111710022B (zh) | 2020-06-08 | 2020-06-08 | 一种避免轮廓相交的快速成型切片处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010515616.8A CN111710022B (zh) | 2020-06-08 | 2020-06-08 | 一种避免轮廓相交的快速成型切片处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111710022A CN111710022A (zh) | 2020-09-25 |
CN111710022B true CN111710022B (zh) | 2021-03-23 |
Family
ID=72539869
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010515616.8A Active CN111710022B (zh) | 2020-06-08 | 2020-06-08 | 一种避免轮廓相交的快速成型切片处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111710022B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112184911B (zh) * | 2020-09-28 | 2023-11-14 | 嘉兴学院 | 一种非流行网格模型切片处理方法 |
CN112214811B (zh) * | 2020-10-09 | 2021-12-31 | 中建八局第三建设有限公司 | 基于bfs算法的地质剖面图封闭区域查找填充方法 |
CN112348946A (zh) * | 2020-11-04 | 2021-02-09 | 深圳市纵维立方科技有限公司 | 打印处理方法、装置、电子设备以及存储介质 |
CN113650301B (zh) * | 2021-08-02 | 2023-04-18 | 嘉兴学院 | 一种基于水平集的3d打印填充路径规划方法 |
CN113597119B (zh) * | 2021-09-30 | 2021-12-10 | 杭州捷配信息科技有限公司 | 一种适用于pcb板生产过程的废料处理方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4775946A (en) * | 1985-03-29 | 1988-10-04 | Hitachi, Ltd. | Method for constructing three-dimensional polyhedron model |
US5387844A (en) * | 1993-06-15 | 1995-02-07 | Micron Display Technology, Inc. | Flat panel display drive circuit with switched drive current |
CN101937579A (zh) * | 2010-09-20 | 2011-01-05 | 南京大学 | 一种利用透视草图创建三维曲面模型的方法 |
CN106228549A (zh) * | 2016-07-14 | 2016-12-14 | 嘉兴学院 | 一种基于路径规划的三角网格牙齿分割方法 |
CN109783943A (zh) * | 2019-01-18 | 2019-05-21 | 安徽建筑大学 | 一种建筑平台智能节点拓扑生成与匹配的方法及*** |
CN110236741A (zh) * | 2019-04-15 | 2019-09-17 | 浙江工业大学 | 一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及个性化髁突假体 |
CN110442994A (zh) * | 2019-08-13 | 2019-11-12 | 嘉兴学院 | 一种基于图论的3d打印切片处理方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101510086B (zh) * | 2009-03-26 | 2011-03-16 | 山东理工大学 | 产品三角网格模型的多轴数控加工刀轨生成方法 |
CN104331555B (zh) * | 2014-10-31 | 2017-07-14 | 河海大学常州校区 | 一种针对带有边界的非封闭stl模型的切片处理方法 |
-
2020
- 2020-06-08 CN CN202010515616.8A patent/CN111710022B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4775946A (en) * | 1985-03-29 | 1988-10-04 | Hitachi, Ltd. | Method for constructing three-dimensional polyhedron model |
US5387844A (en) * | 1993-06-15 | 1995-02-07 | Micron Display Technology, Inc. | Flat panel display drive circuit with switched drive current |
CN101937579A (zh) * | 2010-09-20 | 2011-01-05 | 南京大学 | 一种利用透视草图创建三维曲面模型的方法 |
CN106228549A (zh) * | 2016-07-14 | 2016-12-14 | 嘉兴学院 | 一种基于路径规划的三角网格牙齿分割方法 |
CN109783943A (zh) * | 2019-01-18 | 2019-05-21 | 安徽建筑大学 | 一种建筑平台智能节点拓扑生成与匹配的方法及*** |
CN110236741A (zh) * | 2019-04-15 | 2019-09-17 | 浙江工业大学 | 一种具有拓扑优化固定单元和多孔髁突头单元的个性化髁突假体设计方法及个性化髁突假体 |
CN110442994A (zh) * | 2019-08-13 | 2019-11-12 | 嘉兴学院 | 一种基于图论的3d打印切片处理方法 |
Non-Patent Citations (3)
Title |
---|
Improved slicing algorithm employing nearest distance method;M Vatani等;《Proceedings of the Institution of Mechanical Engineers》;20101231;第224卷(第5期);第745-752页 * |
快速成型中 STL 模型直接切片新算法研究;田仁强等;《机床与液压》;20190831;第47卷(第16期);第3-5节 * |
水平集活动轮廓模型的 3 维牙齿重建;吴婷等;《中国图象图形学报》;20160831;第21卷(第8期);第56-59页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111710022A (zh) | 2020-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111710022B (zh) | 一种避免轮廓相交的快速成型切片处理方法 | |
CN108595858B (zh) | 基于bim的轻量化深度处理工作方法 | |
CN112699623B (zh) | 基于非结构网格规则化重构技术的高精度热流计算方法 | |
CN106373184B (zh) | 一种三维打印模型摆放所需支撑量快速估算方法 | |
CN101853525B (zh) | 基于网格分割的带纹理模型细节保持简化方法 | |
CN106504331A (zh) | 基于三维模型检索的牙齿建模方法 | |
CN110956699B (zh) | 一种三角形网格模型gpu并行切片方法 | |
CN109325316B (zh) | 基于共点焊接排序的stl模型高效并行切层方法 | |
CN113650301B (zh) | 一种基于水平集的3d打印填充路径规划方法 | |
CN100585638C (zh) | 基于线框的曲面体三维边界表示模型重建方法及其装置 | |
CN105793848A (zh) | 包括建模步骤的计算机辅助设计方法 | |
CN108229081A (zh) | 一种组合曲面的几何修复的方法 | |
CN109584371B (zh) | 空间曲线覆盖三角网格曲面的方法 | |
CN106445981B (zh) | 一种基于小波变换的stl网格模型切片数据自适应压缩方法 | |
CN109551768A (zh) | 一种基于stl的3d打印文件的数据处理方法 | |
CN108595631B (zh) | 基于图谱理论的三维cad模型双层检索方法 | |
CN107818338B (zh) | 一种面向地图综合的建筑物群组模式识别的方法及*** | |
CN112734934B (zh) | 一种基于相交边映射的stl模型3d打印切片方法 | |
CN107908880A (zh) | 基于虚拟分割的保拓扑的复杂薄壁件模型中面抽取方法 | |
CN114119628B (zh) | 一种基于特征模板的点云精确分割方法 | |
CN116798028A (zh) | 三维零件的尺寸自动标注方法 | |
CN113360716B (zh) | 一种燃气管网结构逻辑化处理方法及*** | |
CN116416377A (zh) | 针对薄壁管激光切割的加工特征的识别方法、装置及*** | |
Zhu et al. | Direct editing on hexahedral mesh through dual operations | |
CN115797601A (zh) | 一种Brep三维模型的中面自动提取方法及*** |
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 |