CN108279643B - 一种基于测量点与自适应差分进化算法的工件姿态调整方法 - Google Patents

一种基于测量点与自适应差分进化算法的工件姿态调整方法 Download PDF

Info

Publication number
CN108279643B
CN108279643B CN201810081803.2A CN201810081803A CN108279643B CN 108279643 B CN108279643 B CN 108279643B CN 201810081803 A CN201810081803 A CN 201810081803A CN 108279643 B CN108279643 B CN 108279643B
Authority
CN
China
Prior art keywords
actual
theoretical
matching
curved surface
machine tool
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
CN201810081803.2A
Other languages
English (en)
Other versions
CN108279643A (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201810081803.2A priority Critical patent/CN108279643B/zh
Publication of CN108279643A publication Critical patent/CN108279643A/zh
Application granted granted Critical
Publication of CN108279643B publication Critical patent/CN108279643B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/408Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by data handling or data format, e.g. reading, buffering or conversion of data
    • G05B19/4086Coordinate conversions; Other special calculations
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/35Nc in input of data, input till input file format
    • G05B2219/35356Data handling

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Numerical Control (AREA)
  • Automatic Control Of Machine Tools (AREA)

Abstract

本发明公开一种基于测量点与自适应差分进化算法的工件姿态调整方法,包括初始匹配与精确匹配两步;首先,在理论工件曲面与实际曲面上分别选取参考点建立对应的局部坐标系,通过局部坐标系的粗对齐来实现初始匹配;在精确匹配中,通过初始匹配后的点集建立了基于最小二乘距离的目标函数,利用自适应差分进化算法搜索使目标函数值最小的最优解,从而获得最佳精确匹配矩阵;经上述两步,可求得最终的空间变换矩阵,并结合机床结构计算出对应机床的调整量,通过调整机床实现工件姿态的调整,保证工件精确定位与余量的均布。本发明实现了ISO五轴数控机床加工测试件的姿态调整与精确定位,验证了该方法的有效性。

Description

一种基于测量点与自适应差分进化算法的工件姿态调整方法
技术领域
本发明涉及数控加工技术领域,具体为一种基于测量点与自适应差分进化算法的工件姿态调整方法。
背景技术
工件精确定位与余量分布均化是保证加工精度的关键,通过实时调整工件姿态,使其与理论姿态保持一致是实现精确定位与余量均布的直接有效的方法。余量分布对保证加工精度至关重要,而工件姿态对余量分布影响显著。在机械加工中,由于多次装夹等原因,工件姿态容易发生变化,因此有必要寻找一种方法来保证工件姿态的准确性。
关于匹配问题,Besl最先提出了迭代最近点(ICP)算法,该算法被广泛应用于点云配准,在特定的初始条件下能够较好地实现曲面匹配;在传统ICP算法基础上,研究者提出了新的匹配算法,Liu提出了一种新颖的自由曲面匹配算法,该方法通过基于共线性和接近性约束条件建立的传统ICP准则来直接处理可能的点匹配,而不需要任何特征提取、图像预处理或异常数据的运动估计;Du提出了一种基于期望最大化(EM)估计的有界尺度的概率迭代最近点算法,用于噪声点集的各向同性尺度配准;针对自由曲面及其他类型曲面的匹配问题,Zhang提出了一种点到最近点的迭代匹配算法,该算法采用基于距离分布的统计方法实现异常点的处理。然而,就ICP算法而言,其鲁棒性和收敛速度有待提高,且采用该算法获得的解很可能是局部最优解。为了提高ICP算法的收敛速度与运算精度,Lena提出了收敛迭代最近点算法,该方法使得局部曲面匹配的精度大幅提升;针对3D点云,Ge提出了基于非线性Gauss-Helmert 模型和非线性最小二乘的点云数据配准方法,并完成相关三维模型的重构,该方法在收敛速度和配准精度上优于ICP算法。在机械加工中,曲面匹配主要用于复杂零件的定位与余量优化,在曲面匹配的实际应用中,众多研究者也进行了相关的研究工作。Yan为解决数控加工中复杂曲面类零件的定位问题,提出了两步匹配法,算法分为粗匹配与精确匹配,在精确匹配中引入了遗传算法(GA);Ko提出了一种针对自由曲面造型三维物体的匹配方法,该方法将匹配问题归结为采用区间投影法(IPP)求解非线性方程系统,并将曲面的高斯曲率和平均曲率作为匹配特征,该方法可实现全局及局部特征的匹配;Mehrad通过粗定位与精确定位,实现了基于曲率和距离相似性的点到设计模型的匹配,该方法稳定性较好,但由于B样条曲面的重构,该方法在运行效率上有改进空间;Shen针对大型的复杂自由曲面类毛坯件的最佳定位,提出了最佳匹配算法,该算法极大地改善了毛坯件余量分布,提高了加工效率;Zhang提出了基于三坐标测量的加工余量优化模型,该方法通过坐标对齐以确保毛坯件具有足够的加工余量;Sun针对复杂曲面加工提出了一种统一定位模型,并通过BFGS算法实现余量优化,该方法也涉及到了匹配过程,以满足用户定义的需优先保证加工余量的特定表面的约束。
上述匹配算法大多是针对点云与设计模型的匹配,且大多算法需要在已知曲面参数方程或其他曲面特性参数的情况下实现。对于不规则曲面,其曲面参数方程的确定较为困难。
发明内容
针对上述问题,本发明的目的在于提供一种基于测量点与自适应差分进化算法(Self-adaptive Differential Evolution Algorithm,后面简称SADE)的工件姿态调整方法,通过对工件姿态的实时调整,辅助完成工件的精确定位,以确保工件每一工序余量分布均匀,提高加工精度;该方法只需少数测量点,因此对规则与非规则曲面都具有良好的适用性。技术方案如下:
一种基于测量点与自适应差分进化算法的工件姿态调整方法,包括:
进行点集匹配:设工件特征面实际测量点集P={Pi|i=1,2,…,n},理论测量点集Q={Qi|i=1,2,…,n},则P到Q的匹配表示为:
Figure BDA0001561224470000021
式中,R为空间旋转变换矩阵,α、β、γ分别表示实际测量点相对于XYZ轴的旋转角度; T为平移变换矩阵,tx、ty、tz分别表示实际测量点相对于XYZ轴的平移量;匹配流程如下:
Step1:获取测量点集Q和P;
Step2:在理论曲面与对应的实际曲面上分别选取3个对应的参考点,建立理论与实际局部坐标系;
Step3:初始匹配:通过局部坐标系的对齐实现初始匹配,并计算获取初始匹配矩阵R0和T0
Step4:点集变换:利用下式对实际点集P进行变换,获得新的实际点集P';
P′=PR0+T0 (2)
Step5:精确匹配:采用SADE算法搜索最佳精确匹配矩阵R1和T1
Step6:按以下公式计算最终的变换矩阵R和T;
Figure BDA0001561224470000022
Step7:计算机床调整参数并对机床进行调整;
Step8:重新测量并获取新的点集P”,计算点集P”与理论点集Q之间的平均距离dist_ad,并与工件加工面容差ε进行比较,如果dist_ad≤|ε|,结束流程,否则令P'=P”,转至Step 5 重新计算。
进一步的,所述局部坐标系建立的方法为:
在理论曲面上选取不共线的Q0,Q1,Q2作为参考点,在实际曲面选取不共线的P0,P1,P2作为参考点;根据上述参考点确定各局部坐标系XYZ轴单位方向矢量,据此建立局部坐标系,具体如下:
设ia,ja,ka为实际局部坐标系XYZ轴单位方向矢量,P0为原点;it,jt,kt为理论局部坐标系XYZ轴单位方向矢量,Q0为原点;则有:
Figure BDA0001561224470000031
Figure BDA0001561224470000032
根据公式(4)、(5)求得坐标系各轴的单位方向矢量,建立实际局部坐标系P0-XaYaZa和理论局部坐标系Q0-XtYtZt
更进一步的,所述初始匹配矩阵的求解方法包括:
初始匹配通过局部坐标系P0-XaYaZa与Q0-XtYtZt的对齐实现,两坐标系与初始匹配矩阵之间存在以下关系:
Figure BDA0001561224470000033
式中R0为初始匹配矩阵,其形式如下:
Figure BDA0001561224470000034
根据R0,实际曲面与理论曲面之间的方向误差,即实际测量点相对于理论坐标轴Xt,Yt和Zt的旋转角度α00和γ0,如下式所示:
Figure BDA0001561224470000041
式中,R0(i,j)为矩阵R0第i行第j列元素(i,j=1,2,3);
坐标系P0-XaYaZa相对于坐标系Q0-XtYtZt的平移量由下式计算得到:
T0=μ-μ0′R0=[tx0 ty0 tz0]( 8) 式中,
Figure BDA0001561224470000042
分别表示实际测量点集与理论测量点集的质心坐标。
更进一步的,所述实际曲面与理论曲面之间的初始姿态误差通过下式所示姿态误差矩阵 pose0描述:
Figure BDA0001561224470000043
通过pose0能够直观地表示实际曲面相对于理论曲面的旋转角度值与平移量,理想情况下,实际曲面与理论曲面完全重合,因此pose0为零矩阵。
更进一步的,所述精确匹配矩阵的求解方法为:
精确匹配旋转变换与平移变换矩阵分别为R1和T1,其中:
Figure BDA0001561224470000044
T1=[tx1 ty1 tz1]
初始匹配后,实际点集由P转换为P′={P′i|i=1,2,…,n},将精确匹配归结为一个最小化问题,依据最小二乘原理确定其目标函数如下:
Figure BDA0001561224470000045
由式(10)知,目标函数包含α111和tx1,ty1,tz1 6个变量,由式(8)可知,初始匹配后的实际点集P'相对于理论点集Q的平移量可表示为:T1=μ-μ1′R1=[tx1 ty1 tz1],其中P'的质心坐标
Figure BDA0001561224470000046
因此最终的目标函数由式(11)表示;经降元处理,目标函数只含有α1113 个变量;
Figure BDA0001561224470000051
更进一步的,在所述精确匹配中,引入SADE算法对目标函数进行求解,其中的自适应变异算子如下式所示:
Figure BDA0001561224470000052
式中,F为变异算子,F0为初始变异算子,Gm为最大进化代数,G表示当前进化代数。
更进一步的,所述Step7具体为:
通过最终的空间变换矩阵R和T,根据五轴数控机床结构类型求解出对应机床的调整量,包括旋转轴的旋转角度以及平移轴的平移量;根据调整量对机床进行相应的调整,以使工件姿态与匹配结果一致,进而保证余量分布均匀;
在CZFXYB结构的五轴数控机床中,设绕C轴转角为ΔC,绕B轴转角ΔB,沿XYZ轴的平移量分别为Tx、Ty、Tz,则有:
Figure BDA0001561224470000053
T=T0+T1=[Tx Ty Tz] (14)
由式(13)、(14)求得机床调整量ΔB、ΔC及Tx、Ty、Tz,实现机床的调整,进而完成工件姿态的调整。
本发明的有益效果是:本发明的方法分初始匹配与精确匹配,通过两次匹配并结合数控机床的结构,可以求得最佳机床调整参数,最后对机床进行调整以实现工件姿态的调整。可成功用于ISO五轴数控机床典型加工测试件(ISO CD 10791-7:2017,Testconditions for machining centres-Part 7:Accuracy of finished test piece)姿态的调整,确保其定位的准确性与余量分布的均匀性;试验证明,本方法不需要已知工件曲面参数方程或其他参数,只需少数测量点即可完成匹配,实现姿态调整,操作简单可靠,具有更好的适用性;且相对于GA和 BFGS算法,SADE算法更加简便、可靠和稳定,其收敛速度较其他进化算法更快,且不受初值的影响;此外,对大空间、非线性和不可求导的连续问题,其求解效率较其他进化算法好。
附图说明
图1是匹配流程图。
图2是局部坐标系示意。
图3是SADE算法流程图。
图4是测试件模型。
图5是初始姿态示意。
图6是SADE优化曲线图。
图7是姿态调整前后点集间的距离。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细说明。本发明提出了一种新的匹配算法以实现对工件姿态的实时调整。所提算法包括初始匹配与精确匹配,在精确匹配环节引入自适应差分进化算法(SADE)来搜索最佳精确匹配矩阵,通过分步匹配求得实际点集相对于理论点集的最佳空间旋转变换矩阵R和平移变换矩阵T。最后,根据最佳变换矩阵结合机床结构求得机床调整参数,并实现工件姿态的调整。
进行点集匹配:
设工件特征面实际测量点集P={Pi|i=1,2,…,n},理论测量点集Q={Qi|i=1,2,…,n},则P到Q的匹配可以表示为:
Figure BDA0001561224470000061
式中,R为空间旋转变换矩阵,α、β、γ分别表示实际测量点相对于XYZ轴的旋转角度; T为平移变换矩阵,tx、ty、tz分别表示实际测量点相对于XYZ轴的平移量。图1所示为匹配流程图,具体流程如下:
Step1.获取测量点集Q和P。
Step 2.在理论曲面与对应的实际曲面上分别选取3个对应的参考点,建立理论与实际局部坐标系。
如图2所示,理论曲面与实际曲面之间存在姿态误差。在理论曲面与实际曲面上分别选取Q0,Q1,Q2和P0,P1,P2作为参考点。通常情况下,参考点由3个不共线的点构成,且尽可能地分布其他测量点的外侧。根据参考点可以确定各局部坐标系XYZ轴单位方向矢量,据此建立局部坐标系,具体方法如下:
设ia,ja,ka为实际局部坐标系XYZ轴单位方向矢量,P0为原点;it,jt,kt为理论局部坐标系XYZ轴单位方向矢量,Q0为原点。则有:
Figure BDA0001561224470000071
Figure BDA0001561224470000072
根据公式(2)、(3)可以求得坐标系各轴的单位方向矢量,建立如图2所示的实际局部坐标系P0-XaYaZa和理论局部坐标系Q0-XtYtZt
Step3.初始匹配:通过局部坐标系的对齐实现初始匹配,并计算获取初始匹配矩阵R0和T0
初始匹配通过局部坐标系P0-XaYaZa与Q0-XtYtZt的对齐实现,经过初始匹配,可近似确定后续精确匹配中变量的边界值。理论上,初始匹配后实际局部坐标系与理论局部坐标系完全重合,因此两局部坐标系与初始匹配矩阵之间存在式(4)所示关系:
Figure BDA0001561224470000073
式中R0为初始匹配矩阵,其形式如下:
Figure BDA0001561224470000074
根据R0,实际曲面与理论曲面之间的方向误差(即实际测量点相对于理论坐标轴Xt,Yt和 Zt的旋转角度α00和γ0可以求得,其结果如下式(5)所示:
Figure BDA0001561224470000075
式中,R0(i,j)为矩阵R0第i行第j列元素(i,j=1,2,3)。
坐标系P0-XaYaZa相对于Q0-XtYtZt的平移量可由式(6)计算得到。
T0=μ-μ′0R0=[tx0 ty0 tz0] (6)
式中,
Figure BDA0001561224470000081
分别表示实际测量点集与理论测量点集的质心坐标。
求得α000,tx0,ty0,tz0后,实际曲面与理论曲面之间的初始姿态误差可以通过式(7)所示姿态误差矩阵pose0描述:
Figure BDA0001561224470000082
通过pose0可以直观地表示实际曲面相对于理论曲面的旋转角度值与平移量,理想情况下,实际曲面与理论曲面完全重合,因此pose0为零矩阵。
Step4.点集变换:利用公式(8)对实际点集P进行变换,获得新的实际点集P';
P′=PR0+T0 (8)
Step5.精确匹配:采用SADE算法搜索最佳精确匹配矩阵R1和T1
设精确匹配旋转变换与平移变换矩阵分别为R1和T1,其中:
Figure BDA0001561224470000083
T1=[tx1ty1tz1](9)
初始匹配后,实际点集由P转换为P′={Pi′i=1,2,…,n},精确匹配是为了寻找最佳匹配矩阵R1和T1,使匹配后的点集到理论点集Q距离的均方误差最小,因此,精确匹配可以归结为一个最小化问题,依据最小二乘原理可确定其目标函数如式(10)所示。
Figure BDA0001561224470000084
由式(10)可知,目标函数包含α111和tx1,ty1,tz1 6个变量,且属于非线性问题。为了提高运算效率,有必要对该函数进行降元处理,由式(6)可知,初始匹配后的实际点集P'相对于理论点集Q的平移量可以表示为:T1=μ-μ′1R1=[tx1 ty1 tz1],其中P'的质心坐标
Figure BDA0001561224470000085
因此最终的目标函数可由式(11)表示。经降元处理,目标函数只含有α1113个变量,使运算简化。
Figure BDA0001561224470000086
在精确匹配中,本发明首次引入SADE算法对目标函数进行求解,以获得最佳矩阵参数α111和tx1,ty1,tz1。SADE是传统差分进化算法(DE)的改进,其操作主要包括交叉、变异、选择和边界条件处理。在SADE中,可以灵活设计自适应变异算子,能够有效避免因变异率过大造成的搜索效率低下和全局最优解精度低,以及因变异率过小导致的种群多样性降低进而出现早熟等现象,自适应变异算子设计如式(12)所示。
Figure BDA0001561224470000091
式中,F为变异算子,F0为初始变异算子,Gm为最大进化代数,G表示当前进化代数。在算法开始时自适应变异算子为2F0,具有较大值,在初期可以保持个体多样性,避免早熟;随着算法的进展,变异算子逐步降低,到后期其值接近F0,保留优良信息,避免最优解遭到破坏,增加搜索到全局最优解的概率。图3为SADE流程图。
Step6.按公式计算最终的变换矩阵R和T;
Figure BDA0001561224470000092
Step7.计算机床调整参数并对机床进行调整。
根据初始匹配矩阵R0和T0及精确匹配矩阵R1和T1,最终的空间变换矩阵R和T可由式(13)计算获得。
通过R和T,根据五轴数控机床结构类型可以求解出对应机床的调整量,包括旋转轴的旋转角度以及平移轴的平移量。根据调整量对机床进行相应的调整,以使工件姿态与匹配结果一致,进而保证余量分布均匀。以CZFXYB(工作台主轴双回转)结构的五轴数控机床为例,设绕C轴转角为ΔC,绕B轴转角ΔB,沿XYZ轴的平移量分别为Tx、Ty、Tz,则有[20]
Figure BDA0001561224470000093
T=T0+T1=[Tx Ty Tz] (21)
由式(14)、(15)可以求得机床调整量ΔB、ΔC及Tx、Ty、Tz,实现机床的调整,进而完成工件姿态的调整。
以ISO五轴数控机床典型加工测试件为例,对本发明所开发的算法进行验证。图4所示为测试件模型,加工曲面容差为ε=±0.035mm,加工曲面初始姿态如图5所示。具体实施步骤及相关测量数据和试验结果如下。
1)初始匹配
采用三坐标测量机对加工曲面进行测量,获取理论测量点集Q与对应的实际测量点集P,选取参考点建立局部坐标系,测量数据如表1所示。点集初始平均距离dist_in=2.02194mm,初始匹配结果如表2所示。
2)精确匹配
首先,根据初始匹配后姿态误差矩阵pose0中α000的值确定精确匹配中α111的边界值,理论上,其边界值不会超过α000,但考虑到匹配中存在误差,可在α000的基础上适当扩大其边界值。则α111的边界值可分别定义为:[-α0-k10+k1], [-β0-k20+k2],[-γ0-k30+k3],其中ki(i=1,2,3)为[0°,3°]的任意常数。在实例中优化变量边界值如下:
1 (L)1 (U)]=[-5°,5°]
1 (L)1 (U)]=[-5.5°,5.5°]
1 (L)1 (U)]=[-3°,3°]
根据本发明确定的目标函数及优化变量,结合文献[13][17]确定SADE算法相关的参数值如表3所示。图6所示为SADE优化曲线,经过200次进化获得全局最优解α1=-0.02994048°,β1= 0.00881839°,γ1=0.00256112°。精确匹配矩阵R1、T1如下所示。
Figure BDA0001561224470000101
T1=[0.337889015 0.037944696 -0.072290709]
3)机床调整
根据公式(3)可求得最终空间变换矩阵R和T,再结合公式(14)、(15)可以求得机床调整参数ΔC,ΔB,Tx,Ty,Tz。最后对机床进行调整,以实现工件姿态的调整。dist_ad为姿态调整后理论点集与实际点集之间的平均距离,用以评估调整效果。调整参数及调整后的测量结果分别见表4、表5。
4)结果分析
图7所示为匹配前、初始匹配后及姿态调整后实际点集与理论点集之间的平均距离,从中可知,经过两次匹配,实际点集与理论点集间的误差大幅降低,最终收敛到加工曲面容差范围内,调整后点集间距离dist_ad=0.028077mm<|ε|。
表1初始测量数据(dist_in表示初始点集间平均距离/mm)
Figure BDA0001561224470000102
Figure BDA0001561224470000111
Figure BDA0001561224470000121
表2初始匹配结果
Figure BDA0001561224470000122
表3SADE算法参数值
Figure BDA0001561224470000123
表4调整参数及计算结果
Figure BDA0001561224470000124
Figure BDA0001561224470000131
表5初始匹配与调整后的测量点集
Figure BDA0001561224470000132
Figure BDA0001561224470000141

Claims (3)

1.一种基于测量点与自适应差分进化算法的工件姿态调整方法,其特征在于,包括:
进行点集匹配:设工件特征面实际测量点集P={Pi|i=1,2,…,n},理论测量点集Q={Qi|i=1,2,…,n},则P到Q的匹配表示为:
Figure FDA0002455634060000011
式中,R为空间旋转变换矩阵,α、β、γ分别表示实际测量点相对于XYZ轴的旋转角度;T为平移变换矩阵,tx、ty、tz分别表示实际测量点相对于XYZ轴的平移量;匹配流程如下:
Step1:获取测量点集Q和P;
Step2:在理论曲面与对应的实际曲面上分别选取3个对应的参考点,建立理论与实际局部坐标系;
Step3:初始匹配:通过局部坐标系的对齐实现初始匹配,并计算获取初始匹配矩阵R0和T0
Step4:点集变换:利用下式对实际点集P进行变换,获得新的实际点集P';
P′=PR0+T0 (2)
Step5:精确匹配:采用SADE算法搜索最佳精确匹配矩阵R1和T1
Step 6:按以下公式计算最终的变换矩阵R和T;
Figure FDA0002455634060000012
Step7:计算机床调整参数并对机床进行调整;
Step8:重新测量并获取新的点集P”,计算点集P”与理论点集Q之间的平均距离dist_ad,并与工件加工面容差ε进行比较,如果dist_ad≤|ε|,结束流程,否则令P'=P”,转至Step 5重新计算;
所述局部坐标系建立的方法为:
在理论曲面上选取不共线的Q0,Q1,Q2作为参考点,在实际曲面选取不共线的P0,P1,P2作为参考点;根据上述参考点确定各局部坐标系XYZ轴单位方向矢量,据此建立局部坐标系,具体如下:
设ia,ja,ka为实际局部坐标系XYZ轴单位方向矢量,P0为原点;it,jt,kt为理论局部坐标系XYZ轴单位方向矢量,Q0为原点;则有:
Figure FDA0002455634060000021
Figure FDA0002455634060000022
根据公式(4)、(5)求得坐标系各轴的单位方向矢量,建立实际局部坐标系P0-XaYaZa和理论局部坐标系Q0-XtYtZt
所述初始匹配矩阵的求解方法包括:
初始匹配通过局部坐标系P0-XaYaZa与Q0-XtYtZt的对齐实现,两坐标系与初始匹配矩阵之间存在以下关系:
Figure FDA0002455634060000023
式中R0为初始匹配矩阵,其形式如下:
Figure FDA0002455634060000024
根据R0,实际曲面与理论曲面之间的方向误差,即实际测量点相对于理论坐标轴Xt,Yt和Zt的旋转角度α00和γ0,如下式所示:
Figure FDA0002455634060000025
式中,R0(i,j)为矩阵R0第i行第j列元素(i,j=1,2,3);
坐标系P0-XaYaZa相对于坐标系Q0-XtYtZt的平移量由下式计算得到:
T0=μ-μ′0R0=[tx0 ty0 tz0] (8)
式中,
Figure FDA0002455634060000026
分别表示实际测量点集与理论测量点集的质心坐标;
所述精确匹配矩阵的求解方法为:
精确匹配旋转变换与平移变换矩阵分别为R1和T1,其中:
Figure FDA0002455634060000031
T1=[tx1 ty1 tz1]
初始匹配后,实际点集由P转换为P′={Pi′|i=1,2,…,n},将精确匹配归结为一个最小化问题,依据最小二乘原理确定其目标函数如下:
Figure FDA0002455634060000032
由式(10)知,目标函数包含α111和tx1,ty1,tz1 6个变量,由式(8)可知,初始匹配后的实际点集P'相对于理论点集Q的平移量表示为:T1=μ-μ′1R1=[tx1 ty1 tz1],其中P'的质心坐标
Figure FDA0002455634060000033
因此最终的目标函数由式(11)表示;经降元处理,目标函数只含有α111 3个变量;
Figure FDA0002455634060000034
在所述精确匹配中,引入SADE算法对目标函数进行求解,其中的自适应变异算子如下式所示:
Figure FDA0002455634060000035
式中,F为变异算子,F0为初始变异算子,Gm为最大进化代数,G表示当前进化代数。
2.根据权利要求1所述的基于测量点与自适应差分进化算法的工件姿态调整方法,其特征在于,所述实际曲面与理论曲面之间的初始姿态误差通过下式所示姿态误差矩阵pose0描述:
Figure FDA0002455634060000036
通过pose0能够直观地表示实际曲面相对于理论曲面的旋转角度值与平移量,理想情况下,实际曲面与理论曲面完全重合,因此pose0为零矩阵。
3.根据权利要求1所述的基于测量点与自适应差分进化算法的工件姿态调整方法,其特征在于,所述Step7具体为:
通过最终的空间变换矩阵R和T,根据五轴数控机床结构类型求解出对应机床的调整量,包括旋转轴的旋转角度以及平移轴的平移量;根据调整量对机床进行相应的调整,以使工件姿态与匹配结果一致,进而保证余量分布均匀;
在CZFXYB结构的五轴数控机床中,设绕C轴转角为ΔC,绕B轴转角ΔB,沿XYZ轴的平移量分别为Tx、Ty、Tz,则有:
Figure FDA0002455634060000041
T=T0+T1=[Tx Ty Tz] (14)
由式(13)、(14)求得机床调整量ΔB、ΔC及Tx、Ty、Tz,实现机床的调整,进而完成工件姿态的调整。
CN201810081803.2A 2018-01-29 2018-01-29 一种基于测量点与自适应差分进化算法的工件姿态调整方法 Active CN108279643B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810081803.2A CN108279643B (zh) 2018-01-29 2018-01-29 一种基于测量点与自适应差分进化算法的工件姿态调整方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810081803.2A CN108279643B (zh) 2018-01-29 2018-01-29 一种基于测量点与自适应差分进化算法的工件姿态调整方法

Publications (2)

Publication Number Publication Date
CN108279643A CN108279643A (zh) 2018-07-13
CN108279643B true CN108279643B (zh) 2020-08-07

Family

ID=62805341

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810081803.2A Active CN108279643B (zh) 2018-01-29 2018-01-29 一种基于测量点与自适应差分进化算法的工件姿态调整方法

Country Status (1)

Country Link
CN (1) CN108279643B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109773593B (zh) * 2019-01-25 2021-01-22 重庆理工大学 一种基于余量约束条件下的磨削方法
CN109598761B (zh) * 2019-02-16 2022-08-09 西南科技大学 一种曲线测量定位方法
CN110202318B (zh) * 2019-06-18 2021-11-05 华东理工大学 一种基于双侧超声滚压加工的航空叶片定位与姿态调节方法
CN110497088B (zh) * 2019-08-31 2020-05-19 大连理工大学 基于曲面映射的柔性共形天线激光加工误差控制方法
CN112184561B (zh) * 2020-09-25 2024-02-20 合肥工业大学 一种基于多元搜索的火箭级段位姿纠偏方法及其装置
CN112184562B (zh) * 2020-09-25 2024-02-23 合肥工业大学 一种基于近似矩阵的火箭级段位姿纠偏方法及装置
CN112529945B (zh) * 2020-11-17 2023-02-21 西安电子科技大学 一种多视角三维isar散射点集配准方法
CN112817271A (zh) * 2020-12-29 2021-05-18 苏州千机智能技术有限公司 基于在机测量的铸造机匣毛坯加工余量优化方法
CN112923887B (zh) * 2021-01-21 2021-08-20 西南交通大学 任意平面测量控制网非线性最小二乘平差数据检验方法
CN112525161B (zh) * 2021-02-09 2021-05-14 南京景曜智能科技有限公司 一种旋转轴标定方法
CN113405497B (zh) * 2021-05-10 2022-07-12 中航西安飞机工业集团股份有限公司 通过优选与构造测量点提高零组件局部位姿精度的方法
CN113405496B (zh) * 2021-05-10 2022-07-12 中航西安飞机工业集团股份有限公司 通过优选与构造测量点提高零组件整体位姿精度的方法
CN113405511B (zh) * 2021-05-24 2022-07-05 杭州电子科技大学 一种基于ipde算法的关节式坐标测量机标定方法
CN114611362B (zh) * 2022-03-22 2023-08-15 中国工程物理研究院流体物理研究所 一种大型器械工作面的安装调试方法、电子装置及介质
CN114626767A (zh) * 2022-05-16 2022-06-14 中国科学院光电技术研究所 一种离轴变斜率非球面光学元件加工匹配度评价方法
CN118131716A (zh) * 2024-04-30 2024-06-04 成都飞机工业(集团)有限责任公司 一种加工余量评估与优化方法、装置、设备及存储介质

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103913172B (zh) * 2013-12-06 2016-09-21 北京航空航天大学 一种适用于复杂低空下飞行器的路径规划方法
CN106325207B (zh) * 2016-10-08 2019-03-12 南京工业大学 一种五轴数控制齿机床几何误差实际逆向运动学补偿方法

Also Published As

Publication number Publication date
CN108279643A (zh) 2018-07-13

Similar Documents

Publication Publication Date Title
CN108279643B (zh) 一种基于测量点与自适应差分进化算法的工件姿态调整方法
CN111590594B (zh) 基于视觉引导的机器人轨迹跟踪控制方法
CN108563186B (zh) 一种五轴球头铣削几何误差补偿方法
CN103777570B (zh) 基于nurbs曲面的加工误差快速检测补偿方法
CN104408299B (zh) 基于距离识别冗余运动学参数的机器人位置误差补偿方法
CN110125455B (zh) 一种用于机器人钻孔中优化钻头位姿的方法
CN109765848B (zh) 五轴机床运动轴几何误差贡献值计算及影响评估方法
CN108253911B (zh) 一种基于测量点几何特征迭代配准的工件位姿调整方法
CN108803487B (zh) 一种零件侧铣表面的点位轮廓误差预测方法
JP7276788B2 (ja) 5軸数値制御工作機械の主要な幾何学的誤差の補正の比率を最適化させる補正方法
CN110686655B (zh) 三维曲面工件的定位方法、装置及系统
CN112731864B (zh) 一种机床加工误差补偿方法、装置及零件加工机床
CN110625600A (zh) 一种机器人末端工件坐标系标定方法
CN110103071B (zh) 一种变形复杂零件的数字化寻位加工方法
Geng et al. A novel welding path planning method based on point cloud for robotic welding of impeller blades
CN105574884B (zh) 一种医疗机器人dh参数标定方法
CN109978991A (zh) 基于视觉快速实现复杂构件装夹位姿误差在线测量的方法
CN110153582B (zh) 焊接方案生成方法、装置以及焊接系统
CN109773593B (zh) 一种基于余量约束条件下的磨削方法
CN111085902A (zh) 一种视觉在线检测及修正的工件打磨系统
CN112526925B (zh) 基于三维凸轮型面实体化模型偏差补偿的型面精加工方法
CN116587268B (zh) 一种空间大面域机器人铣削加工精度提升方法
CN113240753A (zh) 机器人与双轴变位机构基坐标系标定球面拟合法
CN109035238B (zh) 一种面向自由曲面零件的加工余量离线分析方法
CN110727246A (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