CN109466065B - 面向3d打印的免支撑非均匀三周期极小曲面生成方法及系统 - Google Patents
面向3d打印的免支撑非均匀三周期极小曲面生成方法及系统 Download PDFInfo
- Publication number
- CN109466065B CN109466065B CN201811308405.6A CN201811308405A CN109466065B CN 109466065 B CN109466065 B CN 109466065B CN 201811308405 A CN201811308405 A CN 201811308405A CN 109466065 B CN109466065 B CN 109466065B
- Authority
- CN
- China
- Prior art keywords
- uniform
- support
- printing
- triangular
- supported
- 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
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B29—WORKING OF PLASTICS; WORKING OF SUBSTANCES IN A PLASTIC STATE IN GENERAL
- B29C—SHAPING OR JOINING OF PLASTICS; SHAPING OF MATERIAL IN A PLASTIC STATE, NOT OTHERWISE PROVIDED FOR; AFTER-TREATMENT OF THE SHAPED PRODUCTS, e.g. REPAIRING
- B29C64/00—Additive manufacturing, i.e. manufacturing of three-dimensional [3D] objects by additive deposition, additive agglomeration or additive layering, e.g. by 3D printing, stereolithography or selective laser sintering
- B29C64/30—Auxiliary operations or equipment
- B29C64/386—Data acquisition or data processing for additive manufacturing
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B33—ADDITIVE MANUFACTURING TECHNOLOGY
- B33Y—ADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
- B33Y50/00—Data acquisition or data processing for additive manufacturing
Abstract
本公开公开了一种面向3D打印的免支撑非均匀三周期极小曲面生成方法及系统,根据输入的三维模型、周期函数以及振幅函数计算出模型内部标量场并生成三周期极小曲面相应的非均匀零等值面;对零等值面的隐式表达式中引入一个振幅函数,通过各个三角面片法向量与打印方向,在需要支撑的三角面片区域周围更新振幅函数的值,不断迭代直到获得不需要支撑的模型。
Description
技术领域
本公开涉及一种面向3D打印的免支撑非均匀三周期极小曲面生成方法及系统。
背景技术
本部分的陈述仅仅是提供了与本公开相关的背景技术信息,不必然构成在先技术。
3D打印是一种可以将存储于计算机中的三维数字模型转化为实际物体的新型制造技术,近年来,3D打印技术受到工业和学术界的广泛关注。对于不同类型的打印机,支撑结构的解决方案是不一样的,熔融沉淀式打印机需要大量的支撑结构并且需要在打印后期人工去除。对于粉末打印机,由于他的打印方式是通过粉末凝固形成的因此不需要打印支撑结构处理。由于打印支撑结构耗材无法人工控制以及后期需要人工去除支撑的结构,因此,免支撑在3D打印领域成为一个重要的课题。
三周期极小曲面(Triply Periodic Minimal Surfaces,TPMS)是一种在三维空间中三个独立方向均呈周期性的极小曲面,具有几何形状多样,并可构建参数化数学模型对其进行描述的优点。有一些工作在模型内部加入非均匀的三周期极小曲面使得模型满足特定的物理性质。由于三周期极小曲面是一种准免支撑的结构,在周期尺寸较小时不需要支撑,当尺寸增大时,就必须要考虑支撑的问题。
发明内容
本公开为了解决上述问题,提出了一种面向3D打印的免支撑非均匀三周期极小曲面生成方法及系统,本公开在TPMS的隐式表达式中引入一个振幅函数,根据模型内部的标量场生成相应的三维模型,在需要支撑的三角面片区域周围更新振幅函数的值,从而局部的对三周期极小曲面进行变形,进而不断迭代获得不需要支撑的模型,达到在3D打印中免支撑的效果。
为了实现上述目的,本公开采用如下技术方案:
一种面向3D打印的免支撑非均匀三周期极小曲面生成方法,包括以下步骤:
根据输入的三维模型、周期函数以及振幅函数计算出模型内部标量场并生成三周期极小曲面相应的非均匀零等值面;
对零等值面的隐式表达式中引入一个振幅函数,通过各个三角面片法向量与打印方向,在需要支撑的三角面片区域周围更新振幅函数的值,不断迭代直到获得不需要支撑的模型。
作为进一步的限定,生成非均匀三周期极小曲面的具体过程包括:
定义标量场,利用marching cubes算法重建得到非均匀的三周期极小曲面的mesh表达式。
作为进一步的限定,对于重建出的非均匀的三周期极小曲面表面f,通过各个三角面片法向量与打印方向p的夹角来找到需要支撑的三角面片区域。
作为进一步的限定,当模型表面法向量与打印方向的夹角大于支撑临界角时不需要支撑;同时,对于一块法向量与打印方向的夹角小于支撑临界角的表面区域,如果该区域的表面积小于设定的临界值时,该区域看成是不需要支撑的。
作为更进一步的限定,筛选出需要支撑的三角面片区域的具体过程包括:
(1)计算每个三角面片的法向量与打印方向p的夹角,对于模型夹角大于支撑临界角的三角面片进行标记;
(2)根据被标记的三角面片的连接关系确定不同的联通分支,筛选出需要支撑的表面区域fi,该区域的表面积大于设定的临界值。
作为进一步的限定,如果发现非均匀的三周期极小曲面表面不需要支撑,则输出相应的mesh表达f并退出。
作为进一步的限定,对于根据需要支撑的区域,局部的更新振幅函数。
作为更进一步的限定,对于表面区域fi,局部更新其周围区域其中Pv表示点v的周期,最终得到需要更新的区域X=UXi;更新k′(X)=η·k(X),其中η<1。
一种面向3D打印的免支撑的非均匀三周期极小曲面生成系统,运行于处理器或存储器上,被配置为执行以下步骤:
根据输入的三维模型、周期函数以及振幅函数计算出模型内部标量场并生成三周期极小曲面相应的非均匀零等值面;
对零等值面的隐式表达式中引入一个振幅函数,通过各个三角面片法向量与打印方向,在需要支撑的三角面片区域周围更新振幅函数的值,不断迭代直到获得不需要支撑的模型。
一种3D打印方法,包括上述面向3D打印的免支撑非均匀三周期极小曲面生成方法。
与现有技术相比,本公开的有益效果为:
本公开对于非均匀三周期极小曲面,通过在其隐式表达式中引入一个振幅函数,并且局部的改变振幅函数的值,从而局部的对三周期极小曲面进行变形,达到在3D打印中免支撑的效果。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1为本实施例的方法流程图;
图2a为给定的三维模型kitten;
图2b为Schwarz P-Surface类型的三周极小曲面;
图3a为kitten模型内部非均匀三周期极小曲面,其中需要支撑的三角面片为深色;
图3b筛选掉表面面积小于临界可免支撑最大面积;
图4为最终生成的免支撑的非均匀三周期极小曲面;
具体实施方式:
下面结合附图与实施例对本公开作进一步说明。
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在本公开中,术语如“上”、“下”、“左”、“右”、“前”、“后”、“竖直”、“水平”、“侧”、“底”等指示的方位或位置关系为基于附图所示的方位或位置关系,只是为了便于叙述本公开各部件或元件结构关系而确定的关系词,并非特指本公开中任一部件或元件,不能理解为对本公开的限制。
本公开中,术语如“固接”、“相连”、“连接”等应做广义理解,表示可以是固定连接,也可以是一体地连接或可拆卸连接;可以是直接相连,也可以通过中间媒介间接相连。对于本领域的相关科研或技术人员,可以根据具体情况确定上述术语在本公开中的具体含义,不能理解为对本公开的限制。
首先结合图1对本实施例的一种免支撑的非均匀三周期极小曲面的生成方法进行说明。对于一个初始给定的三维模型S以及周期函数P(x),本实施例在 TPMS隐式表达式中引入一个振幅函数,通过不断迭代的改变振幅函数来对模型内部非均匀三周期极小曲面进行局部变形,从而使极小曲面在打印时不需要支撑。图1为本实施例的流程图,其步骤包括:
步骤一:用户输入三维模型kitten(如图2a)以及非均匀TPMS周期函数P(x) 并初始化振幅函数k(x)=1;
步骤二:根据周期函数P(x)以及振幅函数k(x)计算出模型内部标量场F(x) 并生成相应的零等值面f;
步骤三:利用3D打印自支撑的性质筛选出步骤二重建出的非均匀的TPMS 表面f中需要支撑的三角面片区域;对于重建出的非均匀的TPMS表面f,通过各个三角面片法向量与打印方向p(本实施例中默认打印方向p=(0,0,1))的夹角来找到需要支撑的三角面片区域。
步骤四:如果步骤三发现非均匀的TPMS表面f不需要支撑,则输出三维模型退出最后结果如图4;否则转到步骤五;
步骤五:根据需要支撑的区域,局部的更新振幅函数k(x),转到步骤二。
步骤二中,非均匀三周期极小曲面的生成方法,具体包括以下步骤:
(2-1)在Schwarz P-Surface(如图2b)中引入振幅函数k(x),然后根据周期函数P(x)以及振幅函数k(x)生成一个新的标量场;
(2-2)利用marching cubes算法重建得到非均匀的三周期极小曲面。
具体的,在此以Schwarz P-Surface为例进行说明,其一般形式是由函数的零等值面所定义
其中P(x)表示函数的周期,本实施例通过局部调节函数的振幅来对三周期极小曲面进行调整达到免支撑的要求,因此根据周期函数P(x)以及振幅函数 k(x)生成一个新的标量场
(2-2)利用marching cubes算法重建得到非均匀的三周期极小曲面的mesh 表达f。
步骤三中,需要利用3打印自支撑的性质筛选出需要支撑的三角面片区域,主要步骤包括:
(3-1)标记需要支撑的三角面片,如图3a;
(3-2)根据被标记的三角面片的连接关系确定不同的连通分支,筛选出需要支撑的表面区域,如图3b。
筛选出需要支撑的三角面片区域。为了使模型在打印过程中不需要支撑,对于模型来说有较高的要求。对于3D打印来说,当模型表面法向量与打印方向的夹角大于支撑临界角θthreshold(一般)时不需要支撑;同时,对于一块法向量与打印方向的夹角小于支撑临界角的表面区域f1,如果该区域的表面积 Sf1小于Sthreshold(通过实验得到临界值)时,该区域也可以看成是不需要支撑的。因此筛选需要支撑的三角面片区域的具体包括以下步骤:
(3-1)计算每个三角面片faceti∈f的法向量与打印方向p的夹角θi,对于θi<θthreshold的三角面片进行标记;
(3-2)根据被标记的三角面片的连接关系确定不同的联通分支f1,f2,…,fn,筛选出需要支撑的表面区域fi,其中Sf1>Sthreshold。
步骤五中,通过步骤三获得一些需要支撑的表面区域,局部更新这些区域附近的振幅函数k(x),具体包括以下步骤:
(5-1)对于表面区域fi,本实施例局部更新其周围区域其中Pv表示点v的周期,M表示输入三维模型的内部体空间,最终可以得到需要更新的区域X=UXi;
(5-2)更新k′(X)=η·k(X),其中η<1,本实施例设置η=0.85。转到步骤(2)。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
上述虽然结合附图对本公开的具体实施方式进行了描述,但并非对本公开保护范围的限制,所属领域技术人员应该明白,在本公开的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本公开的保护范围以内。
Claims (5)
1.一种面向3D打印的免支撑非均匀三周期极小曲面生成方法,其特征是:包括以下步骤:
根据输入的三维模型、周期函数以及振幅函数计算出模型内部标量场并生成三周期极小曲面相应的非均匀零等值面;
对非均匀零等值面的隐式表达式中引入一个振幅函数,通过各个三角面片法向量与打印方向p的夹角来找到需要支撑的三角面片区域,在需要支撑的三角面片区域周围更新振幅函数的值,对于表面区域fi,局部更新其周围区域其中Pv表示点v的周期,M表示输入三维模型的内部体空间,最终得到需要更新的区域X=∪Xi;更新k′(X)=η·k(X),其中η<1,不断迭代直到获得不需要支撑的模型;
筛选出需要支撑的三角面片区域的具体过程包括:
(1)计算每个三角面片的法向量与打印方向p的夹角,对于模型夹角大于支撑临界角的三角面片进行标记;
(2)根据被标记的三角面片的连接关系确定不同的联通分支,筛选出需要支撑的表面区域fi,该区域的表面积大于设定的临界值。
2.如权利要求1所述的一种面向3D打印的免支撑非均匀三周期极小曲面生成方法,其特征是:生成非均匀三周期极小曲面的具体过程包括:
定义标量场,利用marching cubes算法重建得到非均匀的三周期极小曲面的mesh表达式。
3.如权利要求1所述的一种面向3D打印的免支撑非均匀三周期极小曲面生成方法,其特征是:如果发现非均匀的三周期极小曲面表面不需要支撑,则输出相应的mesh表达f并退出。
4.一种免支撑的非均匀三周期极小曲面生成系统,其特征是:运行于处理器或存储器上,被配置为执行以下步骤:
根据输入的三维模型、周期函数以及振幅函数计算出模型内部标量场并生成三周期极小曲面相应的非均匀零等值面;
对零等值面的隐式表达式中引入一个振幅函数,通过各个三角面片法向量与打印方向p的夹角来找到需要支撑的三角面片区域,在需要支撑的三角面片区域周围更新振幅函数的值,对于表面区域fi,局部更新其周围区域其中Pv表示点v的周期,M表示输入三维模型的内部体空间,最终得到需要更新的区域X=∪Xi;更新k′(X)=η·k(X),其中η<1,不断迭代直到获得不需要支撑的模型;
筛选出需要支撑的三角面片区域的具体过程包括:
(1)计算每个三角面片的法向量与打印方向p的夹角,对于模型夹角大于支撑临界角的三角面片进行标记;
(2)根据被标记的三角面片的连接关系确定不同的联通分支,筛选出需要支撑的表面区域fi,该区域的表面积大于设定的临界值。
5.一种3D打印方法,其特征是:包括如权利要求1-3中任一项所述的面向3D打印的免支撑非均匀三周期极小曲面生成方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811308405.6A CN109466065B (zh) | 2018-11-05 | 2018-11-05 | 面向3d打印的免支撑非均匀三周期极小曲面生成方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811308405.6A CN109466065B (zh) | 2018-11-05 | 2018-11-05 | 面向3d打印的免支撑非均匀三周期极小曲面生成方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109466065A CN109466065A (zh) | 2019-03-15 |
CN109466065B true CN109466065B (zh) | 2019-12-27 |
Family
ID=65666686
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811308405.6A Active CN109466065B (zh) | 2018-11-05 | 2018-11-05 | 面向3d打印的免支撑非均匀三周期极小曲面生成方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109466065B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111169016B (zh) * | 2019-12-18 | 2021-11-16 | 西北工业大学 | 叶片类零件的3+2轴无支撑3d打印制造方法 |
CN111361145B (zh) * | 2020-03-17 | 2021-06-15 | 华南理工大学 | 一种基于面曝光式的多自由度3d打印方法、装置及系统 |
CN111695259B (zh) * | 2020-06-12 | 2024-04-30 | 重庆大学 | 一种基于3d打印的连续梯度壁厚的tpms结构的加工方法 |
CN111941587B (zh) * | 2020-08-05 | 2021-11-12 | 山东大学 | 一种面向陶瓷3d打印的放样曲面生成方法及系统 |
CN112191849B (zh) * | 2020-10-10 | 2023-03-28 | 重庆邮电大学 | 基于温度分布的梯度多孔散热装置设计及增材制造方法 |
CN112846228B (zh) * | 2020-12-31 | 2023-05-30 | 中核建中核燃料元件有限公司 | 核燃料组件局部下管座免支撑的选区激光熔化成型方法 |
CN113681896B (zh) * | 2021-08-20 | 2023-07-07 | 上海轮廓科技有限公司 | 用于3d打印的方法、装置、设备、存储介质和程序产品 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105904729A (zh) * | 2016-04-22 | 2016-08-31 | 浙江大学 | 一种基于倾斜分层的无支撑三维打印方法 |
WO2016181378A3 (en) * | 2015-05-13 | 2017-11-09 | Massivit 3D Printing Technologies Ltd | A method and apparatus for manufacture of 3d objects |
CN107825701A (zh) * | 2017-09-12 | 2018-03-23 | 上海大学 | 基于五轴打印平台的无支撑3d打印方法 |
CN108327287A (zh) * | 2018-01-16 | 2018-07-27 | 浙江大学 | 一种三周期极小曲面三维打印切片轮廓的快速生成方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10525628B2 (en) * | 2016-04-28 | 2020-01-07 | Wisconsin Alumni Research Foundation | Systems and methods for controlling support structures in manufacturing |
-
2018
- 2018-11-05 CN CN201811308405.6A patent/CN109466065B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016181378A3 (en) * | 2015-05-13 | 2017-11-09 | Massivit 3D Printing Technologies Ltd | A method and apparatus for manufacture of 3d objects |
CN105904729A (zh) * | 2016-04-22 | 2016-08-31 | 浙江大学 | 一种基于倾斜分层的无支撑三维打印方法 |
CN107825701A (zh) * | 2017-09-12 | 2018-03-23 | 上海大学 | 基于五轴打印平台的无支撑3d打印方法 |
CN108327287A (zh) * | 2018-01-16 | 2018-07-27 | 浙江大学 | 一种三周期极小曲面三维打印切片轮廓的快速生成方法 |
Non-Patent Citations (3)
Title |
---|
基于三周期极小曲面的模型轻量化方法;冯紫鑫;《中国优秀硕士学位论文全文数据库》;20170315;第1-56页 * |
复杂曲面约束的薄壁结构优化及其3D打印性研究;江辉;《中国优秀硕士学位论文全文数据库》;20170401;第1-50页 * |
多孔结构的建模方法研究;杨辉;《中国优秀硕士学位论文全文数据库》;20170612;第1-62页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109466065A (zh) | 2019-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109466065B (zh) | 面向3d打印的免支撑非均匀三周期极小曲面生成方法及系统 | |
CN108846148B (zh) | 一种面向3d打印的基于tpms的模型结构优化方法和装置 | |
CN110263384B (zh) | 基于Laplacian微分域变形的三维网格曲面变厚度偏置造型方法 | |
WO2009148157A1 (ja) | 美的曲線を生成する図形情報処理装置、図形情報処理方法及び図形情報処理プログラム | |
US20150127301A1 (en) | Updating A CAD Model To Reflect Global Or Local Shape Changes | |
Fu et al. | Anisotropic simplicial meshing using local convex functions | |
KR20230007513A (ko) | 최소 표면 구조를 적층 제조하는 방법 | |
CN107016214B (zh) | 一种基于有限状态机的参数依赖模型的生成方法 | |
US20210299962A1 (en) | Systems, methods, and media for controlling support structures and build orientation in manufacturing | |
Gupta et al. | Adaptive isogeometric topology optimization using PHT splines | |
CN114077771A (zh) | 气动外形反设计方法 | |
US10573070B1 (en) | Systems and methods for generating a surface that approximates one or more CAD surfaces | |
CN109299499B (zh) | 考虑修正因子的多步骤结构优化设计方法及飞行器 | |
CN109940894B (zh) | 一种基于有限支撑半径控制的卷积曲面混合建模方法 | |
Liu et al. | Automatic sizing functions for unstructured mesh generation revisited | |
Mavriplis | Viscous flow analysis using a parallel unstructured multigrid solver | |
Xie et al. | Geometry-based adaptive mesh generation for continuous and discrete parametric surfaces | |
CN108492372A (zh) | 一种b样条曲面的形状编辑方法 | |
CN110142970B (zh) | 一种用于3d打印技术的壳模型构造方法 | |
CN113536491A (zh) | 点阵结构拓扑优化设计方法、装置及计算机可读存储介质 | |
Kegl et al. | Shape optimal design of elastic planar frames with non‐linear response | |
CN104133923A (zh) | 一种车身快速改型设计方法 | |
CN111274732A (zh) | 一种基于“连接关系-位置”迭代优化的网格修复方法 | |
Lu et al. | Parallel curved mesh adaptation for large scale high-order finite element simulations | |
Zavattieri et al. | Finite element mesh optimization in three dimensions |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20201210 Address after: 311258 No. 1398, Wenbin Road, Wenyan street, Xiaoshan District, Hangzhou, Zhejiang. Patentee after: SHINING 3D TECH Co.,Ltd. Address before: 250101 no.1500 Shunhua Road, high tech Zone, Jinan City, Shandong Province Patentee before: SHANDONG University |