CN111319266A - 一种基于等值线的功能梯度材料3d打印路径规划方法 - Google Patents
一种基于等值线的功能梯度材料3d打印路径规划方法 Download PDFInfo
- Publication number
- CN111319266A CN111319266A CN202010161140.2A CN202010161140A CN111319266A CN 111319266 A CN111319266 A CN 111319266A CN 202010161140 A CN202010161140 A CN 202010161140A CN 111319266 A CN111319266 A CN 111319266A
- Authority
- CN
- China
- Prior art keywords
- gradient
- contour
- projection
- contour line
- printing
- 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
Images
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
Landscapes
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Materials Engineering (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
- Optics & Photonics (AREA)
Abstract
本发明涉及一种基于等值线的功能梯度材料3D打印路径规划方法,其包括以下步骤:1)对功能梯度材料模型进行分层切片处理;2)建立模型的材料分布函数;3)建立材料分布函数在层切平面上的梯度投影;4)根据材料梯度投影特性建立基于材料等值线的打印路径。本发明对功能梯度材料零件建立了统一的材料梯度数学模型,能够根据材料梯度投影特性生成基于等值线的扫描路径。该方法只需存储模型的几何信息和材料函数表达式,而无需复杂的内部结构体素计算和距离变换,避免了后续加工路径的二次离散化,从而实现几何建模、材料表达、路径规划的有效融合。
Description
技术领域
本发明涉及3D打印技术领域,具体涉及一种基于等值线的功能梯度材料3D打印路径规划方法。
背景技术
功能梯度材料(Functionally Gradient Materials,FGM)是应现代航天航空工业等高技术领域的需要,通过设计材料的构成要素从表面到内部呈连续梯度变化,从而使构件的整体性能得以满足而发展起来的一种新型功能材料。在制造FGM零件过程中,由于要考虑零件的形状、材料的成分与分布等信息,FGM零件的制造工艺规划远比匀质材料零件要复杂的多。
现阶段对功能梯度材料的研究主要集中在计算机三维建模上,如体素法、有限元法、B样条法等等,但这些方法均未考虑后续的加工工艺规划问题,将会导致加工过程中材料信息的二次离散化,增加梯度材料信息误差。目前利用3D打印技术来制造FGM零件越来越得到广泛地应用,由于3D打印以数字化的形式逐点逐层累加材料,可以通过控制原材料中不同组分的比例变化方便地实现复杂梯度材料的制备,因此利用3D打印技术制造FGM零件具有其他传统工艺无法比拟的优越性。然而,如何将FGM零件的几何表达、材料设计分布、打印路径规划等方面相互关联成一个完整的体系,以提高打印效率和零件质量,对于推动3D打印技术的应用具有重要意义。
发明内容
针对现有技术存在的不足,本发明的目的在于提供一种基于等值线的功能梯度材料的3D打印路径规划方法,通过建立统一的材料梯度数学模型,利用提取等值线的方法直接生成带有材料属性值的打印路径。
为实现上述目的,本发明提供了如下技术方案:
一种基于等值线的功能梯度材料3D打印路径规划方法,其包括以下步骤:
步骤一、构建FGM零件三维几何模型,将模型所定义的几何空域V划分多个子空间域Vm,且每个子空间域Vm的材料分布不同;然后对模型进行分层切片处理以获取各层切片轮廓数据;
步骤二、建立子空间域Vm中任意一点P(x,y,z)关于其中一种材料的分布函数f(x,y,z),并根据f(x,y,z)获取各个材料在该子空间域Vm中的成分百分比变化公式;
步骤二中:对于子空间域Vm由A、B、C三种材料组成的模型,建立Vm中的任意一点P(x,y,z)关于材料A的分布函数:f(x,y,z):Vm→[s,t],然后根据各材料成分百分比之和为1的关系,获得A、B、C三种材料的成分百分比mA、mB、mC分别为:其中,s、t为函数f(x,y,z)在Vm内的最小值和最大值,mAs、mAt为材料A在f(x,y,z)取得最小值、最大值处的成分百分比,H为常数,0≤H≤1,当H=0,子空间域Vm由两种材料构成。
所述步骤四中梯度投影包括零梯度投影、一维梯度投影和二维梯度投影。
所述二维梯度投影包括等距梯度以及复杂梯度。
当二维梯度投影为等距梯度时,即当等值线上梯度投影的模w为常数时,利用轮廓偏置算法从区域边界开始以Δd为间隔计算等值线扫描路径,然后将材料变化函数f(x,y,z)转化为距离di的函数f(di),其中di为第i条等值线路径到边界的距离,di=(i-1)*Δd,i=1,2,...,M,Δd为相邻路径间距,由打印线宽决定。
当二维梯度投影为复杂梯度时,利用等值线对打印区域进行划分,再在每个子区域内生成平行线扫描路径,其包括以下步骤:
1)将每层切片打印区域进行网格化离散,网格间距由打印线宽决定;
2)按照材料分辨率对材料值f(x,y,z)等间隔采样得到序列:c0,c1c2,...,cN,然后从区域边界开始,在离散网格上线性插值追踪获取等值线:L1,L2,...,LN,其中Li={x,y|f(x,y,zp)=ci},i=1,2,...,N,zp为当前切片z坐标;
3)利用等值线L1,L2,...,LN对切片轮廓区域进行连通区域划分,得到子区域S0,S1,S2,...,SN,然后在每个连通子区域Sj内,j=0,1,2,...,N,生成平行线扫描路径;
4)将每个子区域Sj内路径点的材料值f(Sj)设为相同,对于Sj内任意一点f(Sj)按照如下规则赋值:①若f(x0,y0,z0)<c1,则使f(Sj)=c0;②若ci≤f(x0,y0,z0)<ci+1,i=1,2,...,N-1,则使f(Sj)=ci,③若f(x0,y0,z0)≥cN,则使f(Sj)=cN。
本发明的有益效果:对功能梯度材料零件建立了统一的材料梯度数学模型,能够根据模型的材料梯度特性生成基于等值线的扫描路径。该方法只需存储模型的几何信息和材料梯度函数表达式,而无需复杂的内部结构体素计算和距离变换,避免了后续加工路径的二次离散化,从而实现几何建模、材料表达、路径规划的有效融合。
附图说明
图1是本发明的技术流程示意图。
图2a、b、c是本发明的零梯度投影FGM模型路径规划图。
图3a、b、c是本发明的一维梯度投影FGM模型路径规划图。
图4a、b、c是本发明的等距梯度投影FGM模型路径规划图。
图5a、b、c是本发明的复杂梯度FGM模型路径规划图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明公开了一种基于等值线的功能梯度材料3D打印路径规划方法,如图1所示,其包括以下步骤:
步骤一、对FGM零件模型进行分层切片处理
3D打印工艺采用逐层累加方式制备零件,因此首先构建FGM零件三维几何模型,将模型所定义的几何空域V划分多个子空间域Vm,且每个子空间域Vm的材料分布不同,然后对模型进行分层切片处理,以获取各层切片轮廓数据;
步骤二、建立子空间域Vm中任意一点P(x,y,z)关于其中一种材料的分布函数f(x,y,z),并根据f(x,y,z)获取各个材料在该子空间域Vm中的成分百分比变化公式;
步骤二中:对于子空间域Vm由A、B、C三种材料组成的模型,建立Vm中的任意一点P(x,y,z)关于材料A的分布函数:f(x,y,z):Vm→[s,t],然后根据各材料成分百分比之和为1的关系,获得A、B、C三种材料的成分百分比mA、mB、mC分别为:其中,s、t为函数f(x,y,z)在Vm内的最小值和最大值,mAs、mAt为材料A在f(x,y,z)取得最小值、最大值处的成分百分比,H为常数,0≤H≤1,当H=0,子空间域Vm由两种材料构成。
步骤三、获取f(x,y,z)在点P(x,y,z)的梯度向量:
其中,分别代表f(x,y,z)对x、y、z的偏导数。设3D打印的层切平面垂直与z轴,则材料值f(x,y,z)沿z轴的梯度分量只影响不同层切片的材料值,而不影响同一层内的材料值变化,因此将梯度向量投影在层切平面上,得到梯度投影:
其中梯度投影具体为:
(1)零梯度投影
(2)一维梯度投影
(3)二维梯度投影
除(1)(2)以外为二维梯度,具体又分为:
a.等距梯度
当等值线上梯度投影的模w为常数时,相邻等值线的间距相等。利用轮廓偏置算法从区域边界开始以Δd为间隔计算等值线扫描路径,然后将材料变化函数f(x,y,z)转化为距离di的函数f(di),其中di为第i条等值线路径到边界的距离,di=(i-1)*Δd,i=1,2,...,M,Δd为相邻路径间距,由打印线宽决定。
b.复杂梯度
复杂梯度模型的材料等值线呈不规则变化,采用组合路径:先利用等值线对打印区域进行划分,然后在每个子区域内生成平行线扫描路径。具体为:
1)将每层切片打印区域进行网格化离散,网格间距由打印线宽决定;
2)按照材料分辨率对材料值f(x,y,z)等间隔采样得到序列:c0,c1c2,...,cN,然后从区域边界开始,在离散网格上线性插值追踪获取等值线:L1,L2,...,LN,其中Li={x,y|f(x,y,zp)=ci},i=1,2,...,N,zp为当前切片z坐标;
3)利用等值线L1,L2,...,LN对切片轮廓区域进行连通区域划分,得到子区域S0,S1,S2,...,SN,然后在每个连通子区域Sj内,j=0,1,2,...,N,生成平行线扫描路径;
4)将每个子区域Sj内路径点的材料值f(Sj)设为相同,对于Sj内任意一点f(Sj)按照如下规则赋值:①若f(x0,y0,z0)<c1,则使f(Sj)=c0;②若ci≤f(x0,y0,z0)<ci+1,i=1,2,...,N-1,则使f(Sj)=ci,③若f(x0,y0,z0)≥cN,则使f(Sj)=cN。
具体实施例
实验模型为边长10mm、高度为20mm的正六棱柱。
由于等值线位于x2+y2=c(c为常数)的曲线上,等值线上梯度投影的模也为常数,因此区域V1为等距梯度,采用偏置扫描路径;对于区域V2:f(x,y,z)=1,为零梯度投影,区域V2采用平行线扫描路径,如图4b所示,所有切片层的路径效果如图4c所示。
④当材料分布函数f(x,y,z)=sin(x)+cos(y),模型如图5a所示,
实施例不应视为对本发明的限制,但任何基于本发明的精神所作的改进,都应在本发明的保护范围之内。
Claims (8)
1.一种基于等值线的功能梯度材料3D打印路径规划方法,其特征在于,其包括以下步骤:
步骤一、构建功能梯度材料零件三维几何模型,将模型所定义的几何空域V划分为多个子空间域Vm,且每个子空间域Vm的材料分布不同,然后对模型进行分层切片处理以获取各层切片轮廓数据;
步骤二、建立子空间域Vm中任意一点P(x,y,z)关于其中一种材料的分布函数f(x,y,z),并根据f(x,y,z)获取各材料在该子空间域Vm中的成分百分比变化公式;
3.根据权利要求1所述的一种基于等值线的功能梯度材料3D打印路径规划方法,其特征在于,所述步骤四中梯度投影包括零梯度投影、一维梯度投影和二维梯度投影。
6.根据权利要求3所述的一种基于等值线的功能梯度材料3D打印路径规划方法,其特征在于,所述二维梯度投影包括等距梯度以及复杂梯度。
8.根据权利要求6所述的一种基于等值线的功能梯度材料3D打印路径规划方法,其特征在于,当二维梯度投影为复杂梯度时,利用等值线对打印区域进行划分,再在每个子区域内生成平行线扫描路径,其包括以下步骤:
1)将每层切片打印区域进行网格化离散,网格间距由打印线宽决定;
2)按照材料分辨率对材料值f(x,y,z)等间隔采样得到序列:c0,c1c2,...,cN,然后从区域边界开始,在离散网格上线性插值追踪获取等值线:L1,L2,...,LN,其中Li={x,y|f(x,y,zp)=ci},i=1,2,...,N,zp为当前切片z坐标;
3)利用等值线L1,L2,...,LN对切片轮廓区域进行连通区域划分,得到子区域S0,S1,S2,...,SN,然后在每个连通子区域Sj内,j=0,1,2,...,N,生成平行线扫描路径;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010161140.2A CN111319266B (zh) | 2020-03-10 | 2020-03-10 | 一种基于等值线的功能梯度材料3d打印路径规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010161140.2A CN111319266B (zh) | 2020-03-10 | 2020-03-10 | 一种基于等值线的功能梯度材料3d打印路径规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111319266A true CN111319266A (zh) | 2020-06-23 |
CN111319266B CN111319266B (zh) | 2021-06-22 |
Family
ID=71165733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010161140.2A Active CN111319266B (zh) | 2020-03-10 | 2020-03-10 | 一种基于等值线的功能梯度材料3d打印路径规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111319266B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114274511A (zh) * | 2021-12-02 | 2022-04-05 | 乐金显示光电科技(中国)有限公司 | 掩膜的制备方法 |
CN115891174A (zh) * | 2023-01-06 | 2023-04-04 | 江苏长脖鹿智造科技有限公司 | 三维可视化温控3d打印方法及3d打印机 |
WO2024051010A1 (zh) * | 2022-09-06 | 2024-03-14 | 深圳先进技术研究院 | 组织工程支架的打印方法、打印系统、终端及存储介质 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102446229A (zh) * | 2010-10-11 | 2012-05-09 | 王素 | 基于b样条的功能梯度材料零件建模方法 |
CN104361246A (zh) * | 2014-11-24 | 2015-02-18 | 河南理工大学 | 基于距离场的功能梯度材料隐式建模方法 |
CN106682286A (zh) * | 2016-12-08 | 2017-05-17 | 西北工业大学 | 基于等几何分析法的功能梯度材料零件建模方法 |
CN107901423A (zh) * | 2017-12-11 | 2018-04-13 | 杭州捷诺飞生物科技股份有限公司 | 非均质填充的3d打印方法 |
CN107958111A (zh) * | 2017-11-23 | 2018-04-24 | 西安电子科技大学 | 一种功能梯度材料与结构的参数统计识别方法 |
CN108629833A (zh) * | 2018-05-07 | 2018-10-09 | 四川省有色冶金研究院有限公司 | 一种3d打印模型的结构优化方法 |
CN108803481A (zh) * | 2018-07-30 | 2018-11-13 | 山东大学 | 基于费马螺旋线的自由曲面球头刀具路径规划方法和系统 |
JP2019034487A (ja) * | 2017-08-18 | 2019-03-07 | 富士ゼロックス株式会社 | 画像処理装置及び画像処理プログラム |
CN109508489A (zh) * | 2018-11-07 | 2019-03-22 | 山东大学 | 一种各向异性多孔结构的建模方法及系统 |
KR101949790B1 (ko) * | 2018-10-31 | 2019-05-21 | 새한항업(주) | 수치지형도 속성값을 기반으로 점형, 선형, 면형 타입별 패턴 변화값을 적용한 3d 프린팅용 stl포맷 자동 변환시스템 |
CN110001067A (zh) * | 2019-03-27 | 2019-07-12 | 北京机科国创轻量化科学研究院有限公司 | 一种连续纤维增强复合材料3d打印路径规划方法 |
CN110692084A (zh) * | 2017-05-31 | 2020-01-14 | 惠普发展公司,有限责任合伙企业 | 导出场景的拓扑信息 |
-
2020
- 2020-03-10 CN CN202010161140.2A patent/CN111319266B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102446229A (zh) * | 2010-10-11 | 2012-05-09 | 王素 | 基于b样条的功能梯度材料零件建模方法 |
CN104361246A (zh) * | 2014-11-24 | 2015-02-18 | 河南理工大学 | 基于距离场的功能梯度材料隐式建模方法 |
CN106682286A (zh) * | 2016-12-08 | 2017-05-17 | 西北工业大学 | 基于等几何分析法的功能梯度材料零件建模方法 |
CN110692084A (zh) * | 2017-05-31 | 2020-01-14 | 惠普发展公司,有限责任合伙企业 | 导出场景的拓扑信息 |
JP2019034487A (ja) * | 2017-08-18 | 2019-03-07 | 富士ゼロックス株式会社 | 画像処理装置及び画像処理プログラム |
CN107958111A (zh) * | 2017-11-23 | 2018-04-24 | 西安电子科技大学 | 一种功能梯度材料与结构的参数统计识别方法 |
CN107901423A (zh) * | 2017-12-11 | 2018-04-13 | 杭州捷诺飞生物科技股份有限公司 | 非均质填充的3d打印方法 |
CN108629833A (zh) * | 2018-05-07 | 2018-10-09 | 四川省有色冶金研究院有限公司 | 一种3d打印模型的结构优化方法 |
CN108803481A (zh) * | 2018-07-30 | 2018-11-13 | 山东大学 | 基于费马螺旋线的自由曲面球头刀具路径规划方法和系统 |
KR101949790B1 (ko) * | 2018-10-31 | 2019-05-21 | 새한항업(주) | 수치지형도 속성값을 기반으로 점형, 선형, 면형 타입별 패턴 변화값을 적용한 3d 프린팅용 stl포맷 자동 변환시스템 |
CN109508489A (zh) * | 2018-11-07 | 2019-03-22 | 山东大学 | 一种各向异性多孔结构的建模方法及系统 |
CN110001067A (zh) * | 2019-03-27 | 2019-07-12 | 北京机科国创轻量化科学研究院有限公司 | 一种连续纤维增强复合材料3d打印路径规划方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114274511A (zh) * | 2021-12-02 | 2022-04-05 | 乐金显示光电科技(中国)有限公司 | 掩膜的制备方法 |
WO2024051010A1 (zh) * | 2022-09-06 | 2024-03-14 | 深圳先进技术研究院 | 组织工程支架的打印方法、打印系统、终端及存储介质 |
CN115891174A (zh) * | 2023-01-06 | 2023-04-04 | 江苏长脖鹿智造科技有限公司 | 三维可视化温控3d打印方法及3d打印机 |
Also Published As
Publication number | Publication date |
---|---|
CN111319266B (zh) | 2021-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111319266B (zh) | 一种基于等值线的功能梯度材料3d打印路径规划方法 | |
Nelaturi et al. | Manufacturability feedback and model correction for additive manufacturing | |
CN107067471B (zh) | 一种提高悬垂体模型成型质量的自适应扫描速度方法 | |
CN113650301B (zh) | 一种基于水平集的3d打印填充路径规划方法 | |
CN109101741B (zh) | 一种基于三角网格简化的复杂曲面检测自适应采样方法 | |
Hornus et al. | Variable-width contouring for additive manufacturing | |
CN103366069A (zh) | 一种选择性激光烧结的分层算法 | |
CN115455588A (zh) | 涡轮叶片精铸模具型面反变形设计方法 | |
Zhu et al. | Shape transformation perspective for geometric deviation modeling in additive manufacturing | |
Guan et al. | Development of 3D printing entity slicing software | |
Zhong et al. | Data slicing processing method for RE/RP system based on spatial point cloud data | |
Kumbhar et al. | Improved intermediate point curve model for integrating reverse engineering and rapid prototyping | |
Li et al. | Trajectory generation of spray painting robot using point cloud slicing | |
Kong et al. | An improved method for nurbs free-form surface based on discrete stationary wavelet transform | |
CN114571732A (zh) | 一种适用于复杂曲面类零件的3d打印分层方法 | |
Marzoog et al. | Surface segmentation and reconstruction in reverse engineering | |
CN114119628A (zh) | 一种基于特征模板的点云精确分割方法 | |
Tian et al. | Rapid Manufacturing of Turbine Blades Based on Reverse Engineering and 3D Printing Technology | |
Liao et al. | Subdivision Surface Modeling Technology | |
Yuan et al. | Direct 3d printing system: from point cloud to additive manufacturing | |
Lu et al. | An approximation method combined with conformal transformation for the fast measurement of blade section curves | |
Fu et al. | An adaptive sampling method for STL free-form surfaces based on the quasi-Gauss curvature grid | |
Yang et al. | Spline R-function and applications in FEM | |
Liu et al. | Triangular mesh model reconstruction from scan point clouds based on template | |
Guan et al. | Comparison and Error Analysis of Blade Model Reconstruction |
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 |