CN114001683B - 一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法 - Google Patents

一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法 Download PDF

Info

Publication number
CN114001683B
CN114001683B CN202111243839.4A CN202111243839A CN114001683B CN 114001683 B CN114001683 B CN 114001683B CN 202111243839 A CN202111243839 A CN 202111243839A CN 114001683 B CN114001683 B CN 114001683B
Authority
CN
China
Prior art keywords
wall thickness
coordinate
target
point cloud
curved surface
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
CN202111243839.4A
Other languages
English (en)
Other versions
CN114001683A (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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN202111243839.4A priority Critical patent/CN114001683B/zh
Publication of CN114001683A publication Critical patent/CN114001683A/zh
Application granted granted Critical
Publication of CN114001683B publication Critical patent/CN114001683B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B15/00Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
    • G01B15/02Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons for measuring thickness

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • General Physics & Mathematics (AREA)
  • Numerical Control (AREA)

Abstract

本发明公开了一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法,属于测量‑加工一体化技术领域。该方法是在具有开放式的数控平台上,利用超声测厚装置测得的构件壁厚数据,重新建立剩余壁厚关联的加工目标曲面,完成包容性分析,进行坐标原点平移量解算,通过加工代码原点修正完成加工补偿。本发明克服了现有剩余壁厚补偿方法的不足,满足复杂构件剩余壁厚、轮廓度同时约束的加工要求,提高了加工精度,减小劳动强度,实现复杂构件加工剩余壁厚自适应补偿。

Description

一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿 方法
技术领域
本发明属于测量-加工一体化技术领域,具体涉及一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法。
背景技术
复杂构件是航空航天等领域重大装备上的典型件,如舱段构件等,其加工后的零件剩余壁厚须满足等壁厚或按一定规律变壁厚的加工要求,且剩余壁厚往往是此类构件加工的核心要求。然而在实际机械加工过程中,某些复杂构件受结构等因素影响,半精加工后只能对外表面进行机械加工,其外表面是逼近设计模型的理想形面;而工件的内表面是铸造形面,未经过机械加工,实际的内轮廓与理想形面偏差较大。若按照工件设计模型的坐标原点进行加工,工件剩余壁厚无法满足精度要求。为此,复杂构件加工剩余壁厚补偿方法已成为制约复杂构件加工并亟待突破的关键问题之一。
天津航天长征火箭制造有限公司在发明专利“一种大型薄壁蒙皮自适应等壁厚铣削系统及其加工方法”,CN104289748A中公开一种大型薄壁蒙皮自适应等壁厚铣削系统及其加工方法;大连理工大学在发明专利“一种大型薄壁筒件数字化减薄加工方法”,CN111195830A中公开一种大型薄壁筒件数字化减薄加工方法。然而,上述这两种方法均无法满足零件轮廓和剩余壁厚同时保证的双重约束加工要求。
上述研究均未提及一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法。
发明内容
本发明主要解决的技术难题是克服现有剩余壁厚补偿方法的不足,面向复杂构件剩余壁厚、轮廓度同时约束的加工要求,提供了一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法。该方法是在具有开放式的数控平台上,利用超声测厚装置测得的构件壁厚数据,重新建立剩余壁厚关联的加工目标曲面,完成包容性分析,进行坐标原点平移量解算,通过加工代码原点修正完成加工补偿。本发明可提高加工精度,减小劳动强度,实现复杂构件加工剩余壁厚自适应补偿。
本发明所采用技术方案为:
一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法,该方法首先将超声测厚装置固定在机床主轴上,获取构件外廓测点坐标和对应点壁厚数据;其次,生成工件壁厚模型,根据要求的壁厚计算实际目标曲面,将实际目标曲面与待加工目标曲面进行点云匹配,建立剩余壁厚关联的加工目标曲面;再次,通过计算待加工目标曲面中心与原始坐标中心的位置平移量完成修调量解算;最后,根据平移量进行加工代码补偿修正。具体步骤如下:
第一步,构件毛坯壁厚在机获取
超声测厚装置1由刀柄2进行夹持;刀柄2装于机床主轴3上;所述超声测厚装置1通过超声线缆4与计算机5进行连接,沿轨迹a扫描测量;计算机5实时完成对超声测厚装置1所得信号的处理,获取刀柄2所在机床坐标系测点坐标数据P={Pj(xj,yj,zj),j∈[1,n]}及测点处毛坯壁厚数据H={hj,j∈[1,n]}。其中,Pj表示第j个测点坐标;hj表示Pj对应点处毛坯壁厚,n为刀柄2所在机床坐标系所有测点个数。
第二步,加工余量包容性分析
通过Shepard局部插值函数将测点坐标数据P构造得到工件外表面轮廓S1。基于毛坯壁厚数据H,在工件外表面轮廓S1的法线方向向内构建壁厚,得到工件实际内轮廓面S4。已知设计壁厚数据H’,在实际内轮廓面S4的基础上沿法方向向外构建壁厚,得到实际目标曲面S2
待加工目标曲面S3具有几何不变性,通过与实际目标曲面S2进行点云匹配确定位置。令实际目标曲面S2的坐标数据Q为目标点云,待加工目标曲面S3在机床坐标系下的坐标数据L为源点云。源点云经旋转矩阵R和平移矩阵T变换后与目标点云重合。点qi(xi,yi,zi),li(ai,bi,ci)为分别属于目标点云和源点云的两个坐标点所构成的第i对匹配点坐标,i∈[1,t],t为目标点云和源点云全部匹配点对数。
由于目标点云和源点云经空间变换后无法完全重合,最小化一个目标函数F判断是否满足匹配精度要求。计算目标函数F:
Figure BDA0003320286190000031
其中,F为目标点云Q与源点云L经点云匹配后的平均距离。设
Figure BDA0003320286190000032
Figure BDA0003320286190000033
分别为源点云L与目标点云Q的中心点坐标。中心点的x,y,z轴坐标通过对一组坐标数据内所有坐标点x,y,z轴坐标分别求均值所得。
计算满足匹配精度的阈值k:
Figure BDA0003320286190000034
对源点云L经空间变换后的坐标数据G构造得到待加工目标曲面S3。此时待加工目标曲面S3与工件外表面轮廓S1在法向上的对应点距离即为对应点处加工余量。
第三步,坐标原点平移量解算
计算实际内轮廓面S4中心点O1(x1,y1,z1),以此基准进行壁厚的修调,实现工件整体壁厚的最优,最终保证加工的精度与质量。实际内轮廓面S4中心点O1(x1,y1,z1)到待加工目标曲面中心点O2(x2,y2,z2)之间的平移量即为坐标原点平移量。其中,待加工目标曲面中心点O2(x2,y2,z2)坐标值在X、Y坐标系下取得。只考虑坐标点在同一水平面内的平移,平移量Δx=x2-x1,Δy=y2-y1
第四步,加工代码补偿修正
利用数控系统提供的6个工件坐标系G54~G59进行加工代码补偿修正。若工件半精加工在G54坐标系里面对刀,则精加工时首先修改刀路程序,采用其他不同的坐标系,将平移量输入从而设定平移后的坐标系:
G5·X(Δx) Y(Δy) (3)
其中,G5·表示G55~G59坐标系,Δx,Δy为G5·坐标系原点在G54坐标系中的x轴和y轴坐标。
将式(3)所示命令增加在加工刀路轨迹代码起始位置;在此坐标系下继续精加工刀路完成加工。
本发明的有益效果:本发明利用在机测量的壁厚数据状态,通过加工曲面包容性分析,同时保证了复杂构件轮廓和剩余壁厚精度。本发明提高了加工精度,减小劳动强度,实现了复杂构件加工剩余壁厚自适应补偿。
附图说明
图1为本发明方法全貌示意图。
图2为本发明中所涉及的曲面截面关系图。
图3为基于坐标原点平移的复杂构件自适应补偿方法流程图。
图中:1超声测厚装置;2刀柄;3机床主轴;4超声线缆;5计算机;6待加工构件;7回转工作台;8辅助压板;S1工件外表面轮廓;S2实际目标曲面;S3待加工目标曲面;S4实际内轮廓面;a测量轨迹。
具体实施方式
下面结合附图和技术方案详细说明本发明的实施方式。
待加工构件6,直径500~550mm,毛坯壁厚7~8mm,剩余壁厚3~4mm。超声测厚装置传感器中心频率10MHz,采用点聚焦方式。自适应补偿方法流程图如图3所示,具体步骤如下:
第一步,构件毛坯壁厚在机获取
超声测厚装置1由刀柄2进行夹持;刀柄2装于机床主轴3上;超声测厚装置1通过超声线缆4与计算机5进行连接,沿轨迹a扫描测量;计算机5实时完成对超声测厚装置1所得信号的处理,获取刀柄2所在机床坐标系测点坐标数据P={Pj(xj,yj,zj),j∈[1,n]}及测点处毛坯壁厚数据H={hj,j∈[1,n]}。
第二步,加工余量包容性分析
通过Shepard局部插值函数将测点坐标数据P构造得到工件外表面轮廓S1。基于毛坯壁厚数据H,在工件外表面轮廓S1的法线方向向内构建壁厚,得到工件实际内轮廓面S4。其中计算机床坐标系测点坐标Pj处的法向量
Figure BDA0003320286190000051
时,以Pj为中心点取k个近邻点,坐标表示为Tm(xm,ym,zm),m∈[1,k],Tm表示第m个近邻点。则
Figure BDA0003320286190000052
与由Tm指向Pj的向量
Figure BDA0003320286190000053
垂直。由此最小化一个目标函数估计点云法向量
Figure BDA0003320286190000054
根据视线方向统一法向量方向,解得法向量
Figure BDA0003320286190000055
则工件实际内轮廓面S4对应的坐标数据E={ej(x'j,y'j,z'j),j∈[1,n]},ej表示实际内轮廓面s4的第j个测点坐标,可由下式得到:
Figure BDA0003320286190000061
已知设计壁厚数据H’,在实际内轮廓面S4的基础上沿法方向向外构建壁厚,得到实际目标曲面S2,见图2。
待加工目标曲面S3具有几何不变性,通过与实际目标曲面S2进行点云匹配确定位置。令实际目标曲面S2的坐标数据Q为目标点云,待加工目标曲面S3在机床坐标系下的坐标数据L为源点云。源点云经旋转矩阵R和平移矩阵T变换后与目标点云重合。点qi(xi,yi,zi),li(ai,bi,ci)为分别属于目标点云和源点云的两个坐标点所构成的第i对匹配点坐标,i∈[1,t],t为目标点云和源点云全部匹配点对数。
Figure BDA0003320286190000062
其中,θ为矩阵旋转角度,tx,ty,tz分别为源点云在x,y,z轴方向上的平移量。
由于目标点云和源点云经空间变换后无法完全重合,最小化一个目标函数F判断是否满足匹配精度要求。利用公式(1)计算目标函数F。利用公式(2)计算F满足匹配精度的阈值k。
对源点云L经空间变换后的坐标数据G构造得到待加工目标曲面S3。此时待加工目标曲面S3与工件外表面轮廓S1在法向上的对应点距离即为对应点处加工余量。
第三步,坐标原点平移量解算
计算实际内轮廓面S4中心点O1(x1,y1,z1),以此基准进行壁厚的修调,实现工件整体壁厚的最优,最终保证加工的精度与质量。实际内轮廓面S4中心点O1(x1,y1,z1)到待加工目标曲面中心点O2(x2,y2,z2)之间的平移量即为坐标原点平移量。其中,待加工目标曲面中心点O2(x2,y2,z2)坐标值在X、Y坐标系下取得。如图2所示。只考虑坐标点在同一水平面内的平移,解得平移量Δx=1.2mm,Δy=0.3mm。
第四步,加工代码自适应补偿
加工代码补偿利用数控系统提供的6个工件坐标系G54~G59。若工件半精加工在G54里面对刀,则精加工时首先修改刀路程序,采用其他不同的坐标系如G55,将平移量输入从而设定平移后的坐标系:
G55 X-1.2 Y-0.3 (6)
其中,-1.2、-0.3为G55坐标系原点在G54坐标系中的x轴和y轴坐标,将该命令增加在加工刀路轨迹代码起始位置。在此坐标系下继续精加工刀路完成加工。

Claims (2)

1.一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法,其特征在于,该方法首先将超声测厚装置固定在机床主轴上,获取构件外廓测点坐标和对应点壁厚数据;其次,生成构件壁厚模型,根据要求的壁厚计算实际目标曲面,将实际目标曲面与待加工目标曲面进行点云匹配,建立剩余壁厚关联的待加工目标曲面;再次,通过计算待加工目标曲面中心与原始坐标中心的位置平移量完成修调量解算;最后,根据平移量进行加工代码补偿修正;具体步骤如下:
第一步,构件毛坯壁厚在机获取
刀柄(2)装于机床主轴(3)上;超声测厚装置(1)由刀柄(2)进行夹持,沿轨迹a扫描测量;超声测厚装置(1)通过超声线缆(4)与计算机(5)连接,计算机(5)实时完成对超声测厚装置(1)所得信号的处理,获取刀柄(2)所在机床坐标系测点坐标数据P={Pj(xj,yj,zj),j∈[1,n]}及测点处毛坯壁厚数据H={hj,j∈[1,n]};其中,Pj表示第j个测点坐标;hj表示Pj对应点处毛坯壁厚,n为刀柄所在机床坐标系所有测点个数;
第二步,加工余量包容性分析
通过Shepard局部插值函数将测点坐标数据P构造得到构件外表面轮廓S1;基于毛坯壁厚数据H,在构件外表面轮廓S1的法线方向向内构建壁厚,得到构件实际内轮廓面S4;已知设计壁厚数据H’,在实际内轮廓面S4的基础上沿法方向向外构建壁厚,得到实际目标曲面S2
待加工目标曲面S3具有几何不变性,通过与实际目标曲面S2进行点云匹配确定位置;令实际目标曲面S2的坐标数据Q为目标点云,待加工目标曲面S3在机床坐标系下的坐标数据L为源点云;源点云经旋转矩阵R和平移矩阵T变换后与目标点云重合;
对源点云L经空间变换后的坐标数据G构造得到待加工目标曲面S3;此时待加工目标曲面S3与构件外表面轮廓S1在法向上的对应点距离即为对应点处加工余量;
第三步,坐标原点平移量解算
计算实际内轮廓面S4中心点O1(x1,y1,z1),以此基准进行壁厚的修调,实现构件整体壁厚的最优,最终保证加工的精度与质量;实际内轮廓面S4中心点O1(x1,y1,z1)到待加工目标曲面S3中心点O2(x2,y2,z2)之间的平移量即为坐标原点平移量;其中,待加工目标曲面S3中心点O2(x2,y2,z2)坐标值在X、Y坐标系下取得;只考虑坐标点在同一水平面内的平移,平移量Δx=x2-x1,Δy=y2-y1
第四步,加工代码补偿修正
利用数控系统提供的6个构件坐标系G54~G59进行加工代码补偿修正;若构件半精加工在G54坐标系里面对刀,则精加工时首先修改刀路程序,采用不同的坐标系,输入平移量从而设定平移后的坐标系:
G5·X(Δx) Y(Δy) (3)
其中,G5·表示G55~G59坐标系;Δx,Δy为G5·坐标系原点在G54坐标系中的x轴和y轴坐标;
将式(3)所示命令增加在加工刀路轨迹代码起始位置,在此坐标系下继续精加工刀路完成加工。
2.根据权利要求1所述的方法,其特征在于,所述第二步中,由于目标点云和源点云经空间变换后无法完全重合,最小化一个目标函数F判断是否满足匹配精度要求;计算目标函数F:
Figure FDA0003732113500000021
其中,F为目标点云Q与源点云L经点云匹配后的平均距离;点qi(xi,yi,zi),li(ai,bi,ci)为分别属于目标点云Q和源点云L的两个坐标点所构成的第i对匹配点坐标,i∈[1,t],t为目标点云Q和源点云L全部匹配点对数;设
Figure FDA0003732113500000022
Figure FDA0003732113500000023
分别为目标点云Q与源点云L中心点坐标;中心点的x,y,z轴坐标通过对一组坐标数据内所有坐标点x,y,z轴坐标分别求均值所得;
计算满足匹配精度的阈值k:
Figure FDA0003732113500000024
CN202111243839.4A 2021-10-26 2021-10-26 一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法 Active CN114001683B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111243839.4A CN114001683B (zh) 2021-10-26 2021-10-26 一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111243839.4A CN114001683B (zh) 2021-10-26 2021-10-26 一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法

Publications (2)

Publication Number Publication Date
CN114001683A CN114001683A (zh) 2022-02-01
CN114001683B true CN114001683B (zh) 2022-09-20

Family

ID=79924039

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111243839.4A Active CN114001683B (zh) 2021-10-26 2021-10-26 一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法

Country Status (1)

Country Link
CN (1) CN114001683B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117282995B (zh) * 2023-11-24 2024-02-02 天津航天长征火箭制造有限公司 一种航天旋压箱底加工余量检测、加工方法及电子设备

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111581874A (zh) * 2020-06-28 2020-08-25 大连理工大学 一种薄壁曲面表层图案激光加工轨迹求解方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10357063B3 (de) * 2003-12-04 2005-04-21 Heraeus Tenevo Ag Vertikalziehverfahren zur Herstellung eines zylinderförmigen Glaskörpers und Vorrichtung zur Durchführung des Verfahrens
CN103176428B (zh) * 2013-04-16 2015-01-21 易兴 基于球坐标的cnc系统插补算法及实现该算法的装置
CN107932215B (zh) * 2017-11-01 2019-05-24 中国航发南方工业有限公司 一种薄壁细长轴类零件顶尖孔的加工方法
CN108672719B (zh) * 2018-05-18 2020-06-23 武汉船用机械有限责任公司 一种轴系加工方法
DE102018214002A1 (de) * 2018-08-20 2020-02-20 Sms Group Gmbh Verfahren und Vorrichtung zum Steuern eines Streckreduzierwalzwerks zwecks Wanddickenkompensation
CN109048296B (zh) * 2018-08-29 2019-10-18 航天材料及工艺研究所 一种复合材料舱段被装配型面变形补偿调整方法
CN109489580B (zh) * 2018-12-10 2021-09-24 华东理工大学 一种航空发动机叶片表面加工的在机点云检测及补偿方法
CN111195830B (zh) * 2020-02-03 2021-01-19 大连理工大学 一种大型薄壁筒件数字化减薄加工方法
CN112589133B (zh) * 2020-12-03 2022-07-01 四川航天长征装备制造有限公司 一种大型不规则母线形面壳体均匀减薄加工方法及装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111581874A (zh) * 2020-06-28 2020-08-25 大连理工大学 一种薄壁曲面表层图案激光加工轨迹求解方法

Also Published As

Publication number Publication date
CN114001683A (zh) 2022-02-01

Similar Documents

Publication Publication Date Title
US9817389B2 (en) Adaptively machining component surfaces and hole drilling
Liao et al. Optimization of robot posture and workpiece setup in robotic milling with stiffness threshold
Kosler et al. Adaptive robotic deburring of die-cast parts with position and orientation measurements using a 3D laser-triangulation sensor
CN111590594A (zh) 基于视觉引导的机器人轨迹跟踪控制方法
Zhao et al. Accuracy analysis in mobile robot machining of large-scale workpiece
CN110716497B (zh) 一种基于平面基准约束与余量约束的配准方法
ElMaraghy et al. Integrated inspection and machining for maximum conformance to design tolerances
CN114001683B (zh) 一种基于坐标原点平移的复杂构件加工剩余壁厚自适应补偿方法
WO2020173111A1 (zh) 一种 3d 微涂覆机器人及其涂覆方法
Mikó et al. Experimental verification of cusp heights when 3D milling rounded surfaces
CN113536488A (zh) 基于配准算法的毛坯质量包容性分析和余量优化方法
Zych Programming of welding robots in shipbuilding
Leco et al. A two-step machining and active learning approach for right-first-time robotic countersinking through in-process error compensation and prediction of depth of cuts
Wu et al. Research on machining error transmission mechanism and compensation method for near-net-shaped jet engine blades CNC machining process
CN116551048A (zh) 一种几何自适应加工设备及方法
Xu et al. Tool orientation and redundancy integrated planning method constrained by stiffness for robotic machining of freeform surfaces
Zhu et al. A new calibration method for a dynamic coordinate system in a robotic blade grinding and polishing system based on the six-point limit principle
Chatelain A level-based optimization algorithm for complex part localization
CN116587268B (zh) 一种空间大面域机器人铣削加工精度提升方法
CN109035238B (zh) 一种面向自由曲面零件的加工余量离线分析方法
CN114833848B (zh) 一种刚度约束的机器人铣削加工刀轴矢量和冗余度集成规划方法
CN114131605B (zh) 一种船体零件打磨轨迹自动配准调整装置及方法
Cserteg et al. Multi-operation optimal blank localization for near net shape machining
CN113799126A (zh) 一种基于变形量和灵巧度的机器人加工路径优化方法
Wu et al. Minimum-acceleration local modification method for planning tool orientations in 5-axis ball-end machining

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