CN103777911B - 3d打印中的自适应分层方法 - Google Patents

3d打印中的自适应分层方法 Download PDF

Info

Publication number
CN103777911B
CN103777911B CN201410016504.2A CN201410016504A CN103777911B CN 103777911 B CN103777911 B CN 103777911B CN 201410016504 A CN201410016504 A CN 201410016504A CN 103777911 B CN103777911 B CN 103777911B
Authority
CN
China
Prior art keywords
profile
layer
parameter
thickness
precision
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
CN201410016504.2A
Other languages
English (en)
Other versions
CN103777911A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 ZJU filed Critical Zhejiang University ZJU
Priority to CN201410016504.2A priority Critical patent/CN103777911B/zh
Publication of CN103777911A publication Critical patent/CN103777911A/zh
Application granted granted Critical
Publication of CN103777911B publication Critical patent/CN103777911B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种3D打印中的自适应分层方法,包括:根据相应3D打印设备得到最大允许层厚和最小允许层厚;根据最小允许厚度得到需加工实体的SLC文件;读取SLC文件得到实体的每层轮廓信息;保留第一层轮廓,然后设i初始值为2;计算第i层轮廓和保留下来的第一层轮廓之间的相交多边形集;计算面积比率参数和平均距离参数,根据这两个参数和最大允许层厚判断第i层轮廓是否保留;i增加1,如果i的值不大于实体总的分层数目n,则返回步骤5;否则,计算结束,输出实体自适应分层结果。该发明鲁棒性强,适应于各种形状的模型,尤其是对大型的、结构复杂的模型有效果,同时简单,通用性好,更容易嵌入到3D打印设备中,实现商业化。

Description

3D打印中的自适应分层方法
技术领域
本发明属于3D打印技术分层方法领域,尤其涉及一种3D打印中的自适应分层方法。
背景技术
随着3D打印技术的高速发展,其可以为航空航天、机械、生物、医学等领域的发展提供关键的技术支撑。3D打印技术最大优点就是无需机械加工和任何模具,可以大大的缩短产品的研发周期和降低生产成本。而对于3D打印技术,分层是关键技术之一。分层技术涉及到3D打印的加工时间和产品的精度的平衡问题。当产品精度越高时,分层数目就会越大,这样加工时间就会大大增加,而当分层数目减少时,产品的精度又无法保证。那么就需要合适的分层方法可以在保证精度的前提下尽可能的减少分层数目,以提高加工效率。
目前常用分层方法可以分为三种。一种是等层高分层方法,这种方法在现实加工设备上普遍使用,但是其很难保证产品的表面精度,如果想提高零件精度就会大大增大分层数目;一种是针对STL模型的自适应分层方法,这种方法是根据STL模型中面片的法向量等信息进行分层。但是由于STL模型是由三角面片来表示模型的,本身就存在误差,这样精度就很难保证;最后一种是直接针对CAD实体模型进行自适应分层,用精确地三维模型来进行自适应分层,这样来消除STL模型的误差。例如王卫辰等(基于3D CAD模型表面Z向特征曲线的自适应分层方法;王卫辰,朱真才,方亮,邢邦圣,王繁生;《机械科学与技术》;2010年第29卷第5期;607-611页;共5页)提出一种基于模型表面多处兴趣特征的快速自适应分层算法,该算法利用Z向特征曲线提取CAD模型表面的兴趣特征,通过特征曲线几何特征到层厚的直接映射关系快速确定分层点处的适应性层厚,整个层厚计算过程不存在冗余计算。现有的常用的分层误差计算标准有以下几种:一是尖顶高度准则,根据模型的拓扑信息以及所允许的最大尖顶高度计算得到层厚,这种准则普遍使用,但是当模型法向量接近竖直方向时就会失效;二是面积差比率准则,根据相邻两层轮廓的面积差来判断分层高度是否合适,当面积差很大时就减小分层高度,否则增大分层高度,但是这种方法只能处理简单的模型;其他有体积差比率、表面粗糙度、轮廓边界距离等准则。这些分层方法和误差计算准则都有自己的优势和适应范围。
但是没有一种自适应分层方法可以对于所有的模型进行分层,尤其是大型的、结构复杂的模型,这样也导致现有的自适应分层方法也很难通用化,很难实现商业化。
发明内容
本发明的目的在于提出一种3D打印中的自适应分层方法,借助于模型不同轮廓之间的多边形布尔运算达到自适应分层的目的,本发明提供了一种适合于各种3D打印技术的通用自适应分层方法,鲁棒性强,可以对任何形状的模型进行自适应分层,尤其是对大型的、结构复杂的模型也适用。
一种3D打印中的自适应分层方法,包括如下步骤:
步骤1、根据相应3D打印设备得到最大允许层厚和最小允许层厚;
步骤2、根据得到的最小允许厚度,利用切片软件处理需加工实体的三维模型,得到需加工实体的SLC文件;
步骤3、读取得到的SLC文件,得到实体轮廓信息,包括总的分层数目n、每层轮廓的多边形集以及相应的层高,其中:n是大于1的自然数,层高为每层轮廓在Z轴方向的坐标值;
步骤4、保留第一层轮廓作为初始保留下来的最顶层轮廓,然后设i初始值为2;
步骤5、计算第i层轮廓和保留下来的最顶层轮廓之间的相交多边形集;
步骤6、根据相交多边形集以及第i层轮廓的多边形集,计算第i层轮廓与保留下来的最顶层轮廓之间不重合部分的面积比率参数以及两层轮廓边界之间的平均距离参数,根据这两个参数和层厚来判断第i层轮廓是保留还是删除;
步骤7、i增加1,判断i的值是否大于实体总的分层数目n:如果i不大于n,则返回步骤5,进行下一轮计算;否则,计算结束,输出保留下来的轮廓层作为实体自适应分层结果。
步骤1中每个3D打印设备的结构配置等不同,其加工的最大和最大厚度也不同,需要根据具体的设备去获得这两个参数。
步骤2中SLC文件可以通过各种切片软件获得,如MagicRP、Slic3R、Skeinforge等,该文件是用最小层厚等厚分层得到的。
步骤5中,所述的计算第i层轮廓和保留下来的第一层轮廓之间的相交多边形集时采用多边形布尔运算算法。多边形布尔运算算法能够处理带孔洞的多边形以及包含多个多边形的多边形集,布尔运算得到也是多边形集。
步骤6中,计算第i层与保留下来的最顶层轮廓之间不重合部分的面积比率参数以及两层片轮廓边界之间的平均距离参数,根据这两个参数和层高来判断第i层是保留还是删除,包括步骤:
步骤6.1、计算第i层轮廓与保留下来的最顶层轮廓之间不重合部分的面积比率参数,面积比率参数εi的表达式为:
式中Si是第i层轮廓的面积,Sl是保留的最顶层轮廓的面积,Sins是相交多边形集的面积;
步骤6.2、如果得到的面积比率参数εi大于精度所允许的最大距离参数,则保留第i层轮廓,然后进入步骤7;步骤6.3、如果得到的面积比率参数εi不大于精度所允许的最大距离参数,计算第i层轮廓边界与保留下来的最顶层轮廓之间的平均距离参数,平均距离参数ri的表达式为:
式中louter是相交多边形集中表示外边界的多边形的周长和;
步骤6.4、如果平均距离参数ri大于精度所允许的最大距离参数,则保留第i层轮廓,然后进入步骤7;
步骤6.5、如果平均距离参数ri不大于精度所允许的最大距离参数,计算第i层轮廓与保留下来的最顶层轮廓之间的高度差即第i层轮廓的层厚,并判断该层厚是否大于最大允许层厚:如果大于最大允许层厚,则保留第i层轮廓,然后进入步骤7;如果不大于最大允许层厚,则删除第i层轮廓,然后进入步骤7;
所述的精度所允许的最大比率参数和精度所允许的最大距离参数是根据产品精度要求得到。
与现有技术相比,本发明是3D打印中的自适应分层方法,具体的有益效果是:
本发明通过对用最小允许层厚等厚分层得到的SLC文件进行读取,对两个层片轮廓之间采取多边形的布尔运算去计算各个层片的两个参数值,根据参数值去删除或者保留层片,从而实现自适应分层的目的。该方法鲁棒性强,适应于各种形状的模型,尤其是对大型的、结构复杂的模型有效果,同时简单,易实现,通用性好,更容易嵌入到3D打印设备中,实现商业化。
附图说明
图1为本发明3D打印中的自适应分层方法的流程图;
图2为本发明两层片轮廓边界之间的平均距离参数的示意图;
图3(a)为实施例部分是回用汽车后轮的SLC文件图;
图3(b)和图3(c)为利用本发明3D打印中的自适应分层方法得到分层模型与等厚分层模型对比效果图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步说明。
附图1所示为本发明3D打印技术通用自适应分层方法的流程图。以汽车后轮模型为例,来阐明自适应分层方法。
步骤1、获得某3D打印设备(某SLA设备)的最大允许层厚为0.2mm,最小允许层厚为0.05mm;
步骤2、用切片软件MagicRP获得图3(a)所示的汽车后轮的SLC文件,分层厚度为0.05mm,图(b)为得到的后轮的SLC模型;
步骤3、读取后轮SLC文件,得到后轮轮廓信息,包括总的分层数目n=400、每层轮廓的多边形集以及相应的层高,层高为每层轮廓在Z轴方向的坐标值;
步骤4、保留第一层轮廓(第一层轮廓根据实际打印顺序有关,根据不同打印顺序,选择需要打印的第一层作为第一层轮廓)作为初始保留下来的最顶层轮廓,然后设i初始值为2;
步骤5、对第i层轮廓和保留下来的最顶层轮廓之间进行多边形求交布尔运算得到相交多边形集。轮胎每层轮廓包含有带孔洞的多边形以及多个多边形。选择的多边形布尔运算算法能够处理这些情况,得到相交的多边形集。
步骤6、根据相交多边形集以及轮廓多边形,计算第i层轮廓与保留下来的最顶层轮廓之间不重合部分面积比率参数以及两层片轮廓边界之间的平均距离参数,根据这两个参数和层厚来判断第i层是保留还是删除,包含步骤:
步骤6.1、多边形面积Sp的表达式为:
S p = 1 2 [ ( x 1 y 2 - x 2 y 1 ) + ( x 2 y 3 - x 3 y 2 ) + . . . + ( x i - 1 y i - x i y i - 1 ) . . . + ( x m - 1 y m - x m y m - 1 ) + ( x m y 1 - x 1 y m ) ]
式中xi和yi表示多边形顶点坐标,m表示多边形顶点数。
根据此式得到每层轮廓的面积Slay表达式为
S lqy = Σ j = 1 k S pj
上式中,k表示每层轮廓的多边形个数;Spj为每层轮廓的多边形集中第j个多边形的面积。
那么第i层轮廓与保留下来的最顶层轮廓之间不重合部分的面积比率参数εi的表达式为:
上式中Si是第i层轮廓的面积,Sl是保留的最顶层轮廓的面积,Sins是相交多边形集的面积。计算得到第i层面积比率参数后,判断是否大于精度所允许的最大比率参数,精度所允许的最大面积比率参数是根据产品精度要求得到的,可根据不同产品精度要求确定。
步骤6.2、如果第i层面积比率参数大于精度所允许的最大比率参数,则保留第i层轮廓,然后进入步骤7;
步骤6.3、如果第i层面积比率参数不大于精度所允许的最大比率参数,则计算第i层轮廓与保留下来的最顶层轮廓边界之间的平均距离参数。
图2为两层轮廓边界A之间的平均距离参数示意图,表示了尖顶高度、层高和平均距离参数之间的关系,根据图2(a)得到,尖顶高度表达式为
δ = h · r h 2 + r 2 = h · 1 ( h / r ) 2 + 1
上式中,δ是尖顶高度,h是层厚(通过相邻两层轮廓的Z轴方向的坐标值计算得到),r表示轮廓边界之间的距离。那么尖顶高度是关于r的一个升函数,所以用轮廓边界之间的距离参数可以表示尖顶高度。根据图2(b)得到轮廓边界之间的平均距离参数ri的表达式为:
上式中louter是相交多边形集中表示外边界的多边形的周长和。
计算得到平均距离参数后,判断是否大于精度所允许的最大距离参数。精度所允许的最大距离参数是根据产品精度要求得到的,可根据不同产品精度要求确定。
步骤6.4、如果得到的平均距离参数大于精度所允许的最大距离参数,则保留第i层轮廓,然后进入步骤7;
步骤6.5、如果得到的平均距离参数不大于精度所允许的最大距离参数,计算第i层轮廓与保留下来的最顶层轮廓之间的高度差(z轴坐标差),用两个层片的高度值相减即得到层厚。然后判断该层厚是否大于最大允许层厚:如果大于,则保留第i层轮廓,然后进入步骤7;如果不大于,则删除第i层轮廓,然后进入步骤7。
步骤7、i增加1,判断i的值是否大于实体最大分层数n。如果不大于n,则返回步骤5,进行下一轮计算;否则,计算结束,输出保留下来的轮廓层作为后轮自适应分层结果,图3(c)为后轮自适应分层模型,分层数目为232,最大层厚为0.2mm,最小层厚为0.05mm,其即保证了后轮的表面精度,又减少了分层数目,提高加工效率。

Claims (5)

1.一种3D打印中的自适应分层方法,其特征在于,包括如下步骤:
步骤1、根据相应3D打印设备得到最大允许层厚和最小允许层厚;
步骤2、根据得到的最小允许厚度,利用切片软件处理需加工实体的三维模型,得到需加工实体的SLC文件;
步骤3、读取得到的SLC文件,得到实体轮廓信息,包括总的分层数目n、每层轮廓的多边形集以及相应的层高,其中:n是大于1的自然数,层高为每层轮廓在Z轴方向的坐标值;
步骤4、保留第一层轮廓作为初始保留下来的最顶层轮廓,然后设i初始值为2;i=1,2,…,n;
步骤5、计算第i层轮廓和保留下来的最顶层轮廓之间的相交多边形集;
步骤6、根据相交多边形集以及第i层轮廓的多边形集,计算第i层轮廓与保留下来的最顶层轮廓之间不重合部分的面积比率参数以及两层轮廓边界之间的平均距离参数,根据这两个参数和层厚来判断第i层轮廓是保留还是删除;
具体方法为:
步骤6.1、计算第i层轮廓与保留下来的最顶层轮廓之间不重合部分的面积比率参数,面积比率参数εi的表达式为:
式中Si是第i层轮廓的面积,Sl是保留的最顶层轮廓的面积,Sins是相交多边形集的面积;
步骤6.2、如果得到的面积比率参数εi大于精度所允许的最大距离参数,则保留第i层轮廓,然后进入步骤7;步骤6.3、如果得到的面积比率参数εi不大于精度所允许的最大距离参数,计算第i层轮廓边界与保留下来的最顶层轮廓之间的平均距离参数,平均距离参数ri的表达式为:
式中louter是相交多边形集中表示外边界的多边形的周长和;
步骤6.4、如果平均距离参数ri大于精度所允许的最大距离参数,则保留第i层轮廓,然后进入步骤7;
步骤6.5、如果平均距离参数ri不大于精度所允许的最大距离参数,计算第i层轮廓与保留下来的最顶层轮廓之间的高度差即第i层轮廓的层厚,并判断该层厚是否大于最大允许层厚:如果大于最大允许层厚,则保留第i层轮廓,然后进入步骤7;如果不大于最大允许层厚,则删除第i层轮廓,然后进入步骤7;
所述的精度所允许的最大比率参数和精度所允许的最大距离参数是根据产品精度要求得到;
步骤7、i增加1,判断i的值是否大于实体总的分层数目n:如果i不大于n,则返回步骤5,进行下一轮计算;否则,计算结束,输出保留下来的轮廓层作为实体自适应分层结果。
2.根据权利要求1所述的3D打印中的自适应分层方法,其特征在于,所述的切片软件为MagicRP。
3.根据权利要求1所述的3D打印中的自适应分层方法,其特征在于,所述的切片软件为Skeinforge。
4.根据权利要求1所述的3D打印中的自适应分层方法,其特征在于,所述的切片软件为Slic3R。
5.根据权利要求1所述的3D打印中的自适应分层方法,其特征在于,所述的计算第i层轮廓和保留下来的第一层轮廓之间的相交多边形集时采用多边形布尔运算算法。
CN201410016504.2A 2014-01-14 2014-01-14 3d打印中的自适应分层方法 Active CN103777911B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410016504.2A CN103777911B (zh) 2014-01-14 2014-01-14 3d打印中的自适应分层方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410016504.2A CN103777911B (zh) 2014-01-14 2014-01-14 3d打印中的自适应分层方法

Publications (2)

Publication Number Publication Date
CN103777911A CN103777911A (zh) 2014-05-07
CN103777911B true CN103777911B (zh) 2016-08-17

Family

ID=50570203

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410016504.2A Active CN103777911B (zh) 2014-01-14 2014-01-14 3d打印中的自适应分层方法

Country Status (1)

Country Link
CN (1) CN103777911B (zh)

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104325638B (zh) * 2014-05-14 2016-05-11 浙江大学 基于分区的3d打印填充路径生成方法
US10452053B2 (en) 2014-08-29 2019-10-22 Microsoft Technology Licensing, Llc Fabricating three dimensional objects
CN104331933B (zh) * 2014-10-22 2017-03-15 西安电子科技大学 一种分层方向自适应快速选取方法
CN104527078B (zh) * 2015-01-22 2017-05-24 合肥阿巴赛信息科技有限公司 一种三维打印的可打印性获取方法及系统
CN104881891B (zh) * 2015-04-30 2018-07-27 北京敏速自动控制设备有限公司 面向三维打印的三维模型处理方法及装置
EP3109034B1 (en) * 2015-06-24 2020-07-15 British Telecommunications public limited company Printed logic gate
CN113524689B (zh) * 2016-05-31 2023-05-23 耐克创新有限合伙公司 梯度打印三维结构部件
CN106202687A (zh) * 2016-07-05 2016-12-07 河海大学常州校区 一种基于模型面积变化率的自适应分层处理方法
CN106671422B (zh) * 2016-12-20 2019-05-17 华南理工大学 一种制备生物支架的自适应直接切片方法
CN106898050B (zh) * 2017-02-07 2019-08-23 浙江大学 一种基于环形邻域参照轮廓线的网格模型自适应分层方法
CN107336433B (zh) * 2017-07-13 2019-09-20 武汉工程大学 一种3d打印自适应厚度分层方法
CN107803987B (zh) * 2017-10-18 2020-06-30 湖南华曙高科技有限责任公司 用于增材制造的自适应分层处理方法、系统及增材制造设备
US11865787B2 (en) 2018-02-19 2024-01-09 Agency For Science, Technology And Research Method and system for additive manufacturing
WO2019209262A1 (en) 2018-04-24 2019-10-31 Hewlett-Packard Development Company, L.P. Additive manufacturing system
CN109712240B (zh) * 2018-12-13 2020-10-02 华中科技大学 一种适用于增材制造的曲面分层方法
CN112743040B (zh) * 2019-10-29 2023-09-08 共享智能装备有限公司 3d打印方法、打印机及储存介质
CN113414987A (zh) * 2021-06-23 2021-09-21 哈尔滨理工大学 一种3d打印自适应分层厚度方法
CN113619122B (zh) * 2021-08-25 2023-06-27 珠海赛纳三维科技有限公司 三维物体打印方法、装置、设备及存储介质
CN114030181A (zh) * 2021-09-26 2022-02-11 江苏乾度智造高科技有限公司 一种变层厚3d打印切片方法
CN113844035B (zh) * 2021-10-19 2022-05-13 郑州市中心医院 一种应用于颧骨3d打印中三维图像处理方法及系统
CN115430844B (zh) * 2022-11-09 2023-04-07 四川工程职业技术学院 一种变层厚金属零件激光选区熔化成形方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW200718552A (en) * 2005-11-11 2007-05-16 Wei-Hsiang Lai Apparatus for rapid prototyping using computer-printer aided to object realization
CN103025506A (zh) * 2010-04-25 2013-04-03 奥布吉特有限公司 带外壳物体的实体无模制造

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TW200718552A (en) * 2005-11-11 2007-05-16 Wei-Hsiang Lai Apparatus for rapid prototyping using computer-printer aided to object realization
CN103025506A (zh) * 2010-04-25 2013-04-03 奥布吉特有限公司 带外壳物体的实体无模制造

Also Published As

Publication number Publication date
CN103777911A (zh) 2014-05-07

Similar Documents

Publication Publication Date Title
CN103777911B (zh) 3d打印中的自适应分层方法
CN104331023B (zh) 一种用于五轴数控加工的等残留高度刀触点轨迹的生成及优化处理方法
CN104331933B (zh) 一种分层方向自适应快速选取方法
CN106202687A (zh) 一种基于模型面积变化率的自适应分层处理方法
CA2846084C (en) Geometric modeling of a composite part including a ply-stack up and resin
CN106547969A (zh) 一种基于t样条曲面的三维打印切片方法
CN106021854A (zh) 一种天然河道断面水位流量关系的确定方法
CN109101741B (zh) 一种基于三角网格简化的复杂曲面检测自适应采样方法
CN107067471A (zh) 一种提高悬垂体模型成型质量的自适应扫描速度方法
CN104143209B (zh) 基于线画图案的三维模型雕刻方法
CN105912763A (zh) 基于热传导路径拓扑优化的水下滑翔机巡航路径规划方法
CN102799763A (zh) 一种基于点云姿态标准化的点云线特征提取方法
CN103480534B (zh) 机器人喷涂工件曲面造型方法
CN104268317A (zh) 机械零件圆角结构形状优化方法
CN104750892A (zh) 一种变厚度曲面零件内形面的三维建模方法
CN106293547A (zh) 一种用于3d打印的支撑自动生成方法
CN105107901B (zh) 船舶双向曲率板自动成型加工路径确定及成型效果检测方法
Xu et al. Boundary-conformed machining for trimmed free-form surfaces based on mesh mapping
CN103810313A (zh) 一种stl模型到空间分割模型的转换方法
CN104392030B (zh) 一种基于stl三维模型的曲面加工方法
CN105398056A (zh) 三维打印自适应分层算法
Wang et al. Recognition of freeform surface machining features
Huang et al. Automatic CAD model reconstruction from multiple point clouds for reverse engineering
CN107526859B (zh) 设计复合零件的外表面以用于生产模具的方法、装置及系统
CN103971372A (zh) 人工骨快速成型方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant