CN108062433B - 基于增材再制造点云模型的梯度曲面分层方法 - Google Patents

基于增材再制造点云模型的梯度曲面分层方法 Download PDF

Info

Publication number
CN108062433B
CN108062433B CN201711199464.XA CN201711199464A CN108062433B CN 108062433 B CN108062433 B CN 108062433B CN 201711199464 A CN201711199464 A CN 201711199464A CN 108062433 B CN108062433 B CN 108062433B
Authority
CN
China
Prior art keywords
data
point cloud
distance
gradient
max
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
CN201711199464.XA
Other languages
English (en)
Other versions
CN108062433A (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.)
Academy of Armored Forces of PLA
Original Assignee
Academy of Armored Forces of PLA
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 Academy of Armored Forces of PLA filed Critical Academy of Armored Forces of PLA
Priority to CN201711199464.XA priority Critical patent/CN108062433B/zh
Publication of CN108062433A publication Critical patent/CN108062433A/zh
Application granted granted Critical
Publication of CN108062433B publication Critical patent/CN108062433B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明涉及一种基于增材再制造点云模型的梯度曲面分层方法,包括如下步骤:①成形距离计算;计算对应网格的质心点的欧式距离d:距离d反映了缺损点云在该点处的缺损程度,也就是需要进行增材再制造的成形距离;②梯度距离确定;③分层数据的生成;梯度分层时,真实存在的分层曲面族不需要被确切的描述出来,而只是以各曲面上的分层数据来插值表示就足够了;分层数据包括两部分:轮廓数据和轮廓内部数据;对于第i层的分层数据,其轮廓数据直接记录为分层区间号num=i的数据点;轮廓内部数据P’(x’,y’,z’)则可由低于切层曲面(num<i)的点云采用梯度距离增长得到。

Description

基于增材再制造点云模型的梯度曲面分层方法
技术领域
本发明涉及增材再制造工程技术领域,特别涉及一种基于增材再制造点云模型的分层方法。
背景技术
目前增材再制造技术是一系列用于恢复损伤零(构)件的缺损尺寸和服役性能的先进技术的统称,其对损伤零(构)件尺寸恢复和性能提升的过程是以数字模型驱动为特征的离散―堆积过程。其中,离散过程是以零件缺损部分的三维模型为对象的切片分层过程。零件模型无论是在造型软件中生成还是由反求工程构建,都必须经过分层处理才能输入到增材再制造设备中,因此分层方法是增材再制造过程中的一个关键环节。分层方法不仅影响增材再制造的精度,对其效率也有重要的影响,因此分层方法一直是增材再制造研究的重点和难点问题。
国内外学者提出了基于STL(Stereo Litbography)模型的分层算法,主要有等层厚分层算法、适应性分层算法、斜边分层算法以及曲面分层算法。等层厚分层算法实现简单、程序执行速度快,但台阶效应明显;适应性分层算法采用适应性变化层厚的方法进行分层,有效地减小了台阶效应,但并不能完全消除台阶效应;采用斜边分层或曲面分层等先进分层算法能够完全消除台阶效应,但在系统实现上较为困难。同时,因上述分层算法均是基于STL模型开发的,须经历数据处理中最复杂耗时的曲面重构过程,由此导致了运算过程的复杂化。
发明内容
本发明要解决的技术问题是提供一种基于增材再制造点云模型的精确、高效、台阶效应小的梯度曲面分层方法。
实现本发明目的的技术方案是提供一种基于增材再制造点云模型的梯度曲面分层方法,包括如下步骤:
①成形距离计算;
首先认定增材再制造模型是完整的,由缺损点云和标准点云组成;将增材再制造模型进行空间网格划分,对于任意点P(x,y,z),计算出它的网格号(m,n,q):
Figure RE-GDA0001609557940000021
其中,box_size指网格宽度,xmin,ymin,zmin分别为空间点云在x,y,z方向距离基准点的最小值;
计算网格(m,n,q)的质心点Pp(xx,yy,zz),假设网格中点的个数为sum,则:
Figure RE-GDA0001609557940000022
建立缺损点云所处的网格(m,n,q)与标准点云所处的网格(m’,n’,q’)的对应关系,对应原则是网格号满足:m=m’且n=n’;并计算对应网格的质心点的欧式距离:
Figure RE-GDA0001609557940000023
距离d反映了缺损点云在该点处的缺损程度,也就是需要进行增材再制造的成形距离;选择dmax作为曲面分层的高度范围值;
②梯度距离确定;
2.1)首先确定分层数目n;考虑兼顾到增材再制造的快速性的要求,n值由下式进行估算:
n=int(dmax/hmax)+1;其中hmax为曲面分层梯度距离的最大值;
2.2)然后,按照如下算法确定出点云的分层区间号num和梯度距离hλ
2.2.1)读入网格划分后的增材再制造模型,提取出缺损点云所在网格的质心数据P 以及其对应的成形距离d,构建数据结构D=[P;d;num=NULL;hλ];
2.2.2)令i=0,dimax=dmax+1,输入分层数目n;
2.2.3)令dimin=(n-i)hmin,hmin为曲面分层梯度距离的最小值,遍历num值为NULL的成形距离d,如果dimin≤d<dimax,记录所对应的点P所属的分层区间号num=i,令其分层的梯度距离为:
Figure RE-GDA0001609557940000024
2.2.4)令i=i+1,dimax=dimin,若i<n-1,返回2.2.3);
2.2.5)记录点云中num值为NULL的点,令num=i-1,hλ=d;
③分层数据的生成;
梯度分层时,真实存在的分层曲面族不需要被确切的描述出来,而只是以各曲面上的分层数据来插值表示就足够了;分层数据包括两部分:轮廓数据和轮廓内部数据;对于第i层的分层数据,其轮廓数据直接记录为分层区间号num=i的数据点;轮廓内部数据P’(x’,y’,z’)则可由低于切层曲面(num<i)的点云采用梯度距离增长得到。
进一步的,步骤③中,轮廓内部数据P’(x’,y’,z’)由低于切层曲面(num<i)的点云采用梯度距离增长得到的公式如下:
针对所有数据结构D=[P(x,y,z);d;num;hλ],num=0,1,L,i,
Figure RE-GDA0001609557940000031
进一步的,步骤②中,引入的梯度距离hλ应该处在区间[hminhmax]中,这保证了层与层之间不会出现超出系统再制造能力的厚度;然而,据上计算出的hλ不一定满足该条件,可通过调节[hminhmax]的值来保证;
对公式
Figure RE-GDA0001609557940000032
进行推导,可以计算出第i层区间的hλ满足:
Figure RE-GDA0001609557940000033
可见,随着点云所属的分层区间号i的增加,hλ越界的可能越大;保证hλ不越界的充分必要条件就是:
hmax≥2hmin
所以说,实现梯度曲面分层时的[hminhmax]必须是在满足该条件的前提下进行调节。
本发明具有积极的效果:(1)针对增材再制造过程中数字模型分层算法误差大、效率低、台阶效应明显等问题,本发明提供基于增材再制造点云模型的直接分层方法,具体涉及等距曲面分层方法。通过对增材再制造点云模型直接进行分层,省去了数据处理中复杂耗时的曲面重构过程,既提高了效率、减少了误差来源,又较好的消除了台阶效应。
(2)对于小倾斜度表面的成形,提出本发明的基于增材再制造点云模型的梯度曲面分层法,增材再制造基面点云的增长值可同时作为成形数据,为后续的路径生成提供了一条捷径,保证了增材再制造成形的“快速性”和“精确性”。
附图说明
图1为梯度曲面分层的效果示意图。
具体实施方式
(实施例1)
为了解决小倾斜度表面的成形精度问题,本实施例提出了梯度曲面分层的方法。该方法是指在成形高度范围内构造出若干梯度变化的中间分层面,以实现从再制造基面向成形表面的过渡变化。根据模型特征确定出分层曲面是梯度分层的关键。
对于典型的呈现
Figure RE-GDA0001609557940000043
型或
Figure RE-GDA0001609557940000044
型的再制造模型,为了研究方便,分层方向一般统一定为 Base下的Z轴方向,则此时模型都具备严格的保凸性,如果将再制造基面点云采用梯度距离进行逐层增长就形成了一系列的过渡分层面。同时点云的增长值也构成了分层曲面上的成形点,这还为再制造路径的生成提供了数据基础。
本实施例的基于增材再制造点云模型的梯度曲面分层方法包括如下步骤:
①成形距离计算;
首先认定增材再制造模型是完整的,由缺损点云和标准点云组成。将增材再制造模型进行空间网格划分,对于任意点P(x,y,z),计算出它的网格号(m,n,q):
Figure RE-GDA0001609557940000041
其中,box_size指网格宽度,xmin,ymin,zmin分别为空间点云在x,y,z方向距离基准点的最小值。
计算网格(m,n,q)的质心点Pp(xx,yy,zz),假设网格中点的个数为sum,则:
Figure RE-GDA0001609557940000042
建立缺损点云所处的网格(m,n,q)与标准点云所处的网格(m’,n’,q’)的对应关系,对应原则是网格号满足:m=m’且n=n’。并计算对应网格的质心点的欧式距离:
Figure RE-GDA0001609557940000051
距离d反映了缺损点云在该点处的缺损程度,也就是需要进行增材再制造的成形距离。选择dmax作为曲面分层的高度范围值。
②梯度距离确定;
2.1)首先确定分层数目n。考虑兼顾到增材再制造的快速性的要求,n值由下式进行估算:
n=int(dmax/hmax)+1;其中hmax为曲面分层梯度距离的最大值。
2.2)然后,按照如下算法确定出点云的分层区间号num和梯度距离hλ
2.2.1)读入网格划分后的增材再制造模型,提取出缺损点云所在网格的质心数据P 以及其对应的成形距离d,构建数据结构D=[P;d;num=NULL;hλ];
2.2.2)令i=0,dimax=dmax+1,输入分层数目n;
2.2.3)令dimin=(n-i)hmin,hmin为曲面分层梯度距离的最小值,遍历num值为NULL的成形距离d,如果dimin≤d<dimax,记录所对应的点P所属的分层区间号num=i,令其分层的梯度距离为:
Figure RE-GDA0001609557940000052
2.2.4)令i=i+1,dimax=dimin,若i<n-1,返回2.2.3);
2.2.5)记录点云中num值为NULL的点,令num=i-1,hλ=d。
引入的梯度距离hλ应该处在区间[hminhmax]中,这主要保证了层与层之间不会出现超出系统再制造能力的厚度。然而,据上计算出的hλ不一定满足该条件,可通过调节[hminhmax]的值来保证。
对公式
Figure RE-GDA0001609557940000053
进行推导,可以计算出第i层区间的hλ满足:
Figure RE-GDA0001609557940000054
可见,随着点云所属的分层区间号i的增加,hλ越界的可能越大。保证hλ不越界的充分必要条件就是:
hmax≥2hmin
所以说,实现梯度曲面分层时的[hminhmax]必须是在满足该条件的前提下进行调节,这对于现实的堆焊成形工艺优化具有重要的指导意义。
③分层数据的生成。
梯度分层时,真实存在的分层曲面族不需要被确切的描述出来,而只是以各曲面上的分层数据来插值表示就足够了。分层数据包括两部分:轮廓数据和轮廓内部数据。对于第i层的分层数据,其轮廓数据直接记录为分层区间号num=i的数据点;轮廓内部数据P’(x’,y’,z’)则可由低于切层曲面(num<i)的点云采用梯度距离增长得到,针对所有数据结构D=[P(x,y,z);d;num;hλ],num=0,1,L,i,
Figure RE-GDA0001609557940000061
如果该方法得到的轮廓内部数据,在数量分布上满足要求,就可直接作为成形数据进行保存,省去了后续插值构建的过程,这为再制造路径生成提供了一条捷径。图1即为应用该方法的梯度曲面分层效果。它还可以借鉴到快速成形中去,对于呈现△型或
Figure RE-GDA0001609557940000062
型的曲面分层可以将高于切层面的点云向该层面映射,来生成轮廓内部数据。
显然,上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而这些属于本发明的精神所引伸出的显而易见的变化或变动仍处于本发明的保护范围之中。

Claims (2)

1.一种基于增材再制造点云模型的梯度曲面分层方法,其特征在于包括如下步骤:
①成形距离计算;
首先认定增材再制造模型是完整的,由缺损点云和标准点云组成;将增材再制造模型进行空间网格划分,对于任意点P(x,y,z),计算出它的网格号(m,n,q):
Figure FDA0003040115640000011
其中,box_size指网格宽度,xmin,ymin,zmin分别为空间点云在x,y,z方向距离基准点的最小值;
计算网格(m,n,q)的质心点Pp(xx,yy,zz),假设网格中点的个数为sum,则:
Figure FDA0003040115640000012
建立缺损点云所处的网格(m,n,q)与标准点云所处的网格(m’,n’,q’)的对应关系,对应原则是网格号满足:m=m’且n=n’;并计算对应网格的质心点的欧式距离:
Figure FDA0003040115640000013
距离d反映了缺损点云在该点处的缺损程度,也就是需要进行增材再制造的成形距离;选择dmax作为曲面分层的高度范围值;
②梯度距离确定;
2.1)首先确定分层数目n;考虑兼顾到增材再制造的快速性的要求,n值由下式进行估算:
n=int(dmax/hmax)+1;其中hmax为曲面分层梯度距离的最大值;
2.2)然后,按照如下算法确定出点云的分层区间号num和梯度距离hλ
2.2.1)读入网格划分后的增材再制造模型,提取出缺损点云所在网格的质心数据P以及其对应的距离d,构建数据结构D=[P;d;num=NULL;hλ];
2.2.2)令i=0,dimax=dmax+1,输入分层数目n;
2.2.3)令dimin=(n-i)hmin,hmin为曲面分层梯度距离的最小值,遍历num值为NULL的成形距离d,如果dimin≤d<dimax,记录所对应的点P所属的分层区间号num=i,令其分层的梯度距离为:
Figure FDA0003040115640000021
2.2.4)令i=i+1,dimax=dimin,若i<n-1,返回2.2.3);
2.2.5)记录点云中num值为NULL的点,令num=i-1,hλ=d;
引入的梯度距离hλ处在区间[hmin,hmax]中,这保证了层与层之间不会出现超出系统再制造能力的厚度;然而,计算出的hλ不一定满足该条件:引入的梯度距离hλ处在区间[hmin,hmax]中,可通过调节[hmin hmax]的值来保证;
对公式
Figure FDA0003040115640000022
进行推导,可以计算出第i层区间的hλ满足:
Figure DEST_PATH_FDA0002929943110000024
可见,随着点云所属的分层区间号i的增加,hλ越界的可能越大;保证hλ不越界的充分必要条件就是:
hmax≥2hmin
所以说,实现梯度曲面分层时的[hmin,hmax]必须是在满足该条件:hmax≥2hmin的前提下进行调节;
③分层数据的生成;
梯度分层时,真实存在的分层曲面族不需要被确切的描述出来,而只是以各曲面上的分层数据来插值表示就足够了;分层数据包括两部分:轮廓数据和轮廓内部数据;对于第i层的分层数据,其轮廓数据直接记录为分层区间号num=i的数据点;轮廓内部数据P’(x’,y’,z’)则由低于切层曲面、num<i的点云采用梯度距离增长得到。
2.根据权利要求1所述的基于增材再制造点云模型的梯度曲面分层方法,其特征在于:步骤③中,轮廓内部数据P’(x’,y’,z’)由低于切层曲面、num<i的点云采用梯度距离增长得到的公式如下:
针对所有数据结构D=[P(x,y,z);d;num;hλ],num=0,1,…,i,
Figure FDA0003040115640000031
CN201711199464.XA 2017-11-26 2017-11-26 基于增材再制造点云模型的梯度曲面分层方法 Active CN108062433B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711199464.XA CN108062433B (zh) 2017-11-26 2017-11-26 基于增材再制造点云模型的梯度曲面分层方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711199464.XA CN108062433B (zh) 2017-11-26 2017-11-26 基于增材再制造点云模型的梯度曲面分层方法

Publications (2)

Publication Number Publication Date
CN108062433A CN108062433A (zh) 2018-05-22
CN108062433B true CN108062433B (zh) 2021-06-18

Family

ID=62134940

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711199464.XA Active CN108062433B (zh) 2017-11-26 2017-11-26 基于增材再制造点云模型的梯度曲面分层方法

Country Status (1)

Country Link
CN (1) CN108062433B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108763801B (zh) * 2018-06-04 2022-06-03 辽宁工程技术大学 一种激光增材再制造熔覆层几何特征与稀释率建模方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106547969A (zh) * 2016-11-02 2017-03-29 浙江大学 一种基于t样条曲面的三维打印切片方法
CN106808681A (zh) * 2016-12-30 2017-06-09 华中科技大学 一种提高增材制造零件精度的方法
CN106898050A (zh) * 2017-02-07 2017-06-27 浙江大学 一种基于环形邻域参照轮廓线的网格模型自适应分层方法
CN106985395A (zh) * 2017-03-14 2017-07-28 北京航空航天大学 基于特征的增材制造方法及装置
CN107067469A (zh) * 2016-12-27 2017-08-18 中国人民解放军装甲兵工程学院 获取损伤零件缺损部位模型的系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6584062B2 (ja) * 2014-10-27 2019-10-02 キヤノン株式会社 再生産方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106547969A (zh) * 2016-11-02 2017-03-29 浙江大学 一种基于t样条曲面的三维打印切片方法
CN107067469A (zh) * 2016-12-27 2017-08-18 中国人民解放军装甲兵工程学院 获取损伤零件缺损部位模型的系统
CN106808681A (zh) * 2016-12-30 2017-06-09 华中科技大学 一种提高增材制造零件精度的方法
CN106898050A (zh) * 2017-02-07 2017-06-27 浙江大学 一种基于环形邻域参照轮廓线的网格模型自适应分层方法
CN106985395A (zh) * 2017-03-14 2017-07-28 北京航空航天大学 基于特征的增材制造方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
增材再制造_面向高端在役装备的智能熔敷成形技术;朱胜;《焊接技术》;20170531;第46卷(第5期);70-73 *
复杂曲面正向/逆向快速设计关键技术与增材制造数据处理方法研究;钟山;《中国博士学位论文全文数据库_工程科技Ⅱ辑》;20140615;C029-28 *

Also Published As

Publication number Publication date
CN108062433A (zh) 2018-05-22

Similar Documents

Publication Publication Date Title
CN113348459B (zh) 将网格几何结构转换为水密边界表示的方法、系统及介质
CN108022307B (zh) 基于增材再制造点云模型的自适应平面分层方法
CN106898050B (zh) 一种基于环形邻域参照轮廓线的网格模型自适应分层方法
CN109002639B (zh) 一种基于重位点阵搜索的扭转和倾斜晶界自动建模方法
CN109584357B (zh) 基于多轮廓线的三维建模方法、装置、系统及存储介质
CN109377561A (zh) 一种基于共形几何的数模表面网格生成方法
CN109472870B (zh) 一种基于网格重构和多影响域修正的模型匹配方法
CN110060342B (zh) 一种三维曲面拟合方法
CN107038308B (zh) 一种基于线性内插的规则格网地形建模方法
CN102831283B (zh) 一种基于曲面特征的复杂产品模型构建方法
CN107886573B (zh) 一种复杂地质条件下边坡三维有限元网格生成方法
CN108062433B (zh) 基于增材再制造点云模型的梯度曲面分层方法
CN111079326A (zh) 二维各向异性网格单元度量张量场光滑化方法
Wang et al. An adaptive slicing thickness adjustment method based on cloud point in 3D printing
CN113342999A (zh) 一种基于多层跳序树结构的变分辨率点云简化方法
CN104036552B (zh) 基于最远点优化的蓝噪声网格生成方法
CN105096379A (zh) 一种基于k-邻域的三角网格曲面重建方法
CN116664758A (zh) 一种基于断层面的四面体网格自适应分级细化方法
Remacle et al. Fast and robust mesh generation on the sphere–application to coastal domains
Navangul et al. A vertex translation algorithm for adaptive modification of STL file in layered manufacturing
KR101235506B1 (ko) 초기 클러스터링을 사용한 병합식 계층적 클러스터링 방법 및 장치
Ye et al. Information reuse to accelerate customized product slicing for additive manufacturing
Wang et al. A slicing algorithm to guarantee non-negative error of additive manufactured parts
Li et al. Printing orientation optimization of 3D model
CN111524236A (zh) 一种基于平面分析的自适应三角网简化方法

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