CN111044289B - 基于闭环动态测量的大型高速回转装备对准误差测量方法 - Google Patents

基于闭环动态测量的大型高速回转装备对准误差测量方法 Download PDF

Info

Publication number
CN111044289B
CN111044289B CN201911370022.6A CN201911370022A CN111044289B CN 111044289 B CN111044289 B CN 111044289B CN 201911370022 A CN201911370022 A CN 201911370022A CN 111044289 B CN111044289 B CN 111044289B
Authority
CN
China
Prior art keywords
error
representing
matrix
total
coordinate
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
CN201911370022.6A
Other languages
English (en)
Other versions
CN111044289A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute 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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201911370022.6A priority Critical patent/CN111044289B/zh
Publication of CN111044289A publication Critical patent/CN111044289A/zh
Application granted granted Critical
Publication of CN111044289B publication Critical patent/CN111044289B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M15/00Testing of engines
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明提出了基于闭环动态测量的大型高速回转装备对准误差测量方法,所述测量方法包括以下步骤:将每一个部件认为是一个理想刚体,三个位移自由度δx,δy,δz和三个角度自由度θx,θy,θz,计算出各个部件的几何偏差;将各个部件的几何偏差带入到误差累计求和公式中,借助多刚体系统理论和坐标的矩阵变换计算出多个相互连接的组件的误差累积量,得到总的误差Etotal;结合计算得到的总的误差Etotal对传感器的安装姿态误差进行补偿。本发明的基于闭环动态测量的大型高速回转装备对准误差测量方法,可以计算得到终端的对准误差值,进而通过软件算法将该误差补偿掉,以提高航空发动机单级转子的测量精度。

Description

基于闭环动态测量的大型高速回转装备对准误差测量方法
技术领域
本发明涉及基于闭环动态测量的大型高速回转装备对准误差测量方法,属于测量技术领域。
背景技术
在三维空间测量大尺寸物体时,引入了各种辅助设备以完成全面型测量,例如:升降机构、旋转机构、平移机构等。但是引入的辅助测量设备也会引入测量误差,相比于不同的辅助设备所引入的误差类型以及具体误差值是不同的,这就给最终的测量精度带来了困难,难以建立精确的误差分离模型以消除这些误差所带来的影响。因此从封闭环连接拓扑结构出发,建立各组件的误差传递机理模型,进而实现对误差传递的控制很有必要。
目前主要的误差抑制方法有以下几种:(1)采用精度更高的零部件产品,但各零部件产品精度的提高在一定程度会提高整个装配体的精度,因此会在一定程度上提高装配体的性能指标等,但是这种方法作用有限,无法一味依靠其来提高装配体精度,特别是当机械加工精度达到现有物理极限时候,提高精度所带来的费用是非常高额的;(2)通过外部更高一级的测量设备精确测量出单个组件误差的量值,这种情况存在一个问题,就是容易忽略重要误差和部分关键误差测量困难的问题,此外,不能够考虑各误差之间的耦合关系以及累积效应。因此该方法有一定缺陷。
发明内容
本发明提出基于闭环动态测量的大型高速回转装备对准误差测量方法,所提出的方法可以针对具体的测量设备建立起各误差的传递模型,通过在终端软件算法补偿该误差值,实现最终误差的抑制。这种方法的针对性更强,能够快速准确地抑制测量过程中由于各运动组件对最终测量的影响。同时避免了盲目一味地提高各组件的制造加工精度所带来的高额费用,具有一定的经济意义。实现快速准确的测量误差抑制。
基于闭环动态测量的大型高速回转装备对准误差测量方法,所述测量方法包括以下步骤:
步骤一、将每一个部件认为是一个理想刚体,则它在三维空间拥有六个自由度,三个位移自由度δx,δy,δz和三个角度自由度θx,θy,θz,计算出各个部件的几何偏差;
步骤二、将各个部件的几何偏差带入到误差累计求和公式中,通过借助多刚体系统理论和坐标的矩阵变换计算出多个相互连接的组件的误差累积量,得到总的误差Etotal
步骤三、结合计算得到的总的误差Etotal对传感器的安装姿态误差进行补偿,即通过软件编程算法调整传感器与被测工件之间的相对位姿。
进一步的,在步骤二中,所述矩阵为空间位姿坐标变换矩阵表达为:
Figure BDA0002339425790000021
其中,l=[lx ly lz]T,m=[mx my mz]T和n=[nx ny nz]T表示x,y,z方向向量坐标;p=[px py pz]T表示坐标点平移位置坐标;R为3×3旋转矩阵;P为3×1平移矩阵。
进一步的,将空间位姿坐标变换矩阵改写为如下表达式:
Figure BDA0002339425790000022
Tij表示经过旋转变换和平移变换后的位置坐标矩阵,(a,b,c)表示水平位移平移量,
Figure BDA0002339425790000023
Figure BDA0002339425790000024
表示旋转角度变换量,
一共两条尺寸链条,一条是传感器尺寸链条,另一条是待测工件尺寸链条,根据多体系统理论建立起传递模型,进而根据闭环尺寸链条计算出误差:
Figure BDA0002339425790000025
E表示测量起点和测量终端的偏置误差,Pw表示工件末端位置坐标,Ps表示传感器测量末端坐标位置,P0表示同一起点坐标位置矩阵,Tideal表示理想坐标变换矩阵,Terror表示含有的误差变换矩阵,n表示工件尺寸链的个数,m表示传感器尺寸链中部件的数目,
沿着坐标轴误差分量的计算:
Figure BDA0002339425790000026
本发明具有以下有益效果:本发明的基于闭环动态测量的大型高速回转装备对准误差测量方法,可以计算得到终端的对准误差值,进而通过软件算法将该误差补偿掉,以提高航空发动机单级转子的测量精度,为精准堆叠装配奠定坚实的基础。
附图说明
图1为本发明的基于闭环动态测量的大型高速回转装备对准误差测量方法的方法流程图;
图2为被测部件的装配堆叠过程示意图;
图3为测量设备各组件连接关系图。
具体实施方式
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是建立起针对航空发动机多级转子闭环对准测量误差的计算模型,不单单的通过测量误差进而校正误差,重点对于耦合的误差也进行分析计算。通过终端误差的量级和方向进而采取一定的措施进行误差补偿。典型的测量系统是由多个零部件,大型设备的装配应在同一空间坐标系统,这样可以保证装配精度,但难以测量和测量成本变得更高的体积增加由于工件进行测试。为了保证产品在整个空间的装配精度,将引入精密转盘、精密导轨等实现全方位的辅助装配。不幸的是,这种操作将引入更多的错误,并导致较低的精度。如果不能准确地预测和控制误差,将对后续的装配产生巨大的影响。
参照图1-图3所示,本发明提出了基于闭环动态测量的大型高速回转装备对准误差测量方法的一实施例,所述测量方法包括以下步骤:
步骤一、将每一个部件认为是一个理想刚体,则它在三维空间拥有六个自由度,三个位移自由度δx,δy,δz和三个角度自由度θx,θy,θz,计算出各个部件的几何偏差;
步骤二、将各个部件的几何偏差带入到误差累计求和公式中,通过借助多刚体系统理论和坐标的矩阵变换计算出多个相互连接的组件的误差累积量,得到总的误差Etotal
步骤三、结合计算得到的总的误差Etotal对传感器的安装姿态误差进行补偿,即通过软件编程算法调整传感器与被测工件之间的相对位姿。
在本部分优选实施例中,在步骤二中,所述矩阵为空间位姿坐标变换矩阵表达为:
Figure BDA0002339425790000031
其中,l=[lx ly lz]T,m=[mx my mz]T和n=[nx ny nz]T表示x,y,z方向向量坐标;p=[px py pz]T表示坐标点平移位置坐标;R为3×3旋转矩阵;P为3×1平移矩阵。
在本部分优选实施例中,将空间位姿坐标变换矩阵改写为如下表达式:
Figure BDA0002339425790000041
Tij表示经过旋转变换和平移变换后的位置坐标矩阵,(a,b,c)表示水平位移平移量,
Figure BDA0002339425790000042
Figure BDA0002339425790000043
表示旋转角度变换量,
一共两条尺寸链条,一条是传感器尺寸链条,另一条是待测工件尺寸链条,根据多体系统理论建立起传递模型,进而根据闭环尺寸链条计算出误差:
Figure BDA0002339425790000044
E表示测量起点和测量终端的偏置误差,Pw表示工件末端位置坐标,Ps表示传感器测量末端坐标位置,P0表示同一起点坐标位置矩阵,Tideal表示理想坐标变换矩阵,Terror表示含有的误差变换矩阵,n表示工件尺寸链的个数,m表示传感器尺寸链中部件的数目,
沿着坐标轴误差分量的计算:
Figure BDA0002339425790000045
本发明通过减小了测量系统的误差使得单级转子的测量精度提高,航空发动机多级转子是由单级转子堆叠而成的,因此,该方法可以直接提高后续转子系统整体的装配精度。

Claims (1)

1.基于闭环动态测量的大型高速回转装备对准误差测量方法,其特征在于,所述测量方法包括以下步骤:
步骤一、将每一个部件认为是一个理想刚体,则它在三维空间拥有六个自由度,三个位移自由度δx,δy,δz和三个角度自由度θx,θy,θz,计算出各个部件的几何偏差;
步骤二、将各个部件的几何偏差带入到误差累计求和公式中,通过借助多刚体系统理论和坐标的矩阵变换计算出多个相互连接的组件的误差累积量,得到总的误差Etotal
步骤三、结合计算得到的总的误差Etotal对传感器的安装姿态误差进行补偿,即通过软件编程算法调整传感器与被测工件之间的相对位姿,
在步骤二中,所述矩阵为空间位姿坐标变换矩阵表达为:
Figure FDA0003182059500000011
其中,l=[lx ly lz]T,m=[mx my mz]T和n=[nx ny nz]T表示x,y,z方向向量坐标;p=[pxpypz]T表示坐标点平移位置坐标;R为3×3旋转矩阵;P为3×1平移矩阵,
将空间位姿坐标变换矩阵改写为如下表达式:
Figure FDA0003182059500000012
Tij表示经过旋转变换和平移变换后的位置坐标矩阵,(a,b,c)表示水平位移平移量,
Figure FDA0003182059500000013
Figure FDA0003182059500000014
表示旋转角度变换量,
一共两条尺寸链条,一条是传感器尺寸链条,另一条是待测工件尺寸链条,根据多体系统理论建立起传递模型,进而根据闭环尺寸链条计算出误差:
Figure FDA0003182059500000015
E表示测量起点和测量终端的偏置误差,Pw表示工件末端位置坐标,Ps表示传感器测量末端坐标位置,P0表示同一起点坐标位置矩阵,Tideal表示理想坐标变换矩阵,Terror表示含有的误差变换矩阵,n表示工件尺寸链的个数,m表示传感器尺寸链中部件的数目,
沿着坐标轴误差分量的计算:
Figure FDA0003182059500000016
CN201911370022.6A 2019-12-26 2019-12-26 基于闭环动态测量的大型高速回转装备对准误差测量方法 Active CN111044289B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911370022.6A CN111044289B (zh) 2019-12-26 2019-12-26 基于闭环动态测量的大型高速回转装备对准误差测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911370022.6A CN111044289B (zh) 2019-12-26 2019-12-26 基于闭环动态测量的大型高速回转装备对准误差测量方法

Publications (2)

Publication Number Publication Date
CN111044289A CN111044289A (zh) 2020-04-21
CN111044289B true CN111044289B (zh) 2021-09-03

Family

ID=70239108

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911370022.6A Active CN111044289B (zh) 2019-12-26 2019-12-26 基于闭环动态测量的大型高速回转装备对准误差测量方法

Country Status (1)

Country Link
CN (1) CN111044289B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115265912B (zh) * 2022-07-13 2023-08-01 哈尔滨工业大学 基于双偏置误差传递的大型高速回转装备重心测量方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100901520B1 (ko) * 2008-10-10 2009-06-08 버추얼모션(주) 다물체 동역학 시스템의 엔터티별 자코비안계산 프로세스의 소요시간의 차이로 발생하는 멀티스레드간의 동작시간의 갭을 줄이기 위한 병렬처리 프로세싱 방법
CN104596466A (zh) * 2015-02-04 2015-05-06 厦门大学 大口径光学非球面元件的两段轮廓拼接测量方法
CN107845114A (zh) * 2017-11-10 2018-03-27 北京三快在线科技有限公司 地图的构建方法、装置及电子设备
CN108681634A (zh) * 2018-05-15 2018-10-19 华中科技大学 一种考虑实际工况下的机械产品的运行精度预测方法
EP3467842A1 (en) * 2017-10-09 2019-04-10 Digital Orthopaedics Method for simulating foot and/or ankle
CN109815562A (zh) * 2019-01-07 2019-05-28 南京航空航天大学 一种基于容差空间的装配位姿优化方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100901520B1 (ko) * 2008-10-10 2009-06-08 버추얼모션(주) 다물체 동역학 시스템의 엔터티별 자코비안계산 프로세스의 소요시간의 차이로 발생하는 멀티스레드간의 동작시간의 갭을 줄이기 위한 병렬처리 프로세싱 방법
CN104596466A (zh) * 2015-02-04 2015-05-06 厦门大学 大口径光学非球面元件的两段轮廓拼接测量方法
EP3467842A1 (en) * 2017-10-09 2019-04-10 Digital Orthopaedics Method for simulating foot and/or ankle
CN107845114A (zh) * 2017-11-10 2018-03-27 北京三快在线科技有限公司 地图的构建方法、装置及电子设备
CN108681634A (zh) * 2018-05-15 2018-10-19 华中科技大学 一种考虑实际工况下的机械产品的运行精度预测方法
CN109815562A (zh) * 2019-01-07 2019-05-28 南京航空航天大学 一种基于容差空间的装配位姿优化方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Intelligent measurement and compensation of linear motor force ripple: A projection-based learning approach in the presence of noise;Yang Liu等;《Measurement Science and Technology》;20181231;全文 *
大型试件的精确定位方法与调控系统的研究;谭久彬等;《计量学报》;19950422;全文 *
树形多刚体系统动力学方程;汪恩松;《应用数学和力学》;19851227;全文 *

Also Published As

Publication number Publication date
CN111044289A (zh) 2020-04-21

Similar Documents

Publication Publication Date Title
CN110802585B (zh) 机械臂末端传感器补偿方法及接触力/力矩测量方法
CN111367236B (zh) 一种面向加工过程的移动机器人系统标定方法及系统
CN109822577B (zh) 一种基于视觉伺服的移动式机器人高精度加工方法
EP1274546B1 (en) Pathcorrection for an industrial robot
CN108818536A (zh) 一种机器人手眼标定的在线偏移修正方法及装置
CN109015110B (zh) 一种机床回转轴位置误差建模与辨识方法
CN104977901B (zh) 三轴运动平台改进型交叉耦合控制装置及方法
CN111037542B (zh) 一种逆动力学控制机器人直线加工的轨迹误差补偿方法
CN112327623A (zh) 一种基于负载摆动状态观测的双摆型吊车滑模控制方法
CN102692873A (zh) 工业机器人定位精度标定方法
CN109571477B (zh) 一种改进的机器人视觉与传送带综合标定方法
CN112318498B (zh) 一种考虑参数耦合的工业机器人标定方法
CN111044289B (zh) 基于闭环动态测量的大型高速回转装备对准误差测量方法
CN109291051B (zh) 一种基于角度传感器的串并联机器人末端姿态闭环控制方法
CN115290289B (zh) 一种提高大迎角尾撑试验系统控制精度的优化方法
Saund et al. High accuracy articulated robots with CNC control systems
CN110549333B (zh) 一种用于TriMule卧式混联机器人的重力补偿方法
CN110370271A (zh) 工业串联机器人的关节传动比误差校准方法
CN111489399A (zh) 一种视觉跟踪组件安装参数标定装置和方法
CN110722562A (zh) 一种用于机器人参数辨识的空间雅克比矩阵构造方法
CN114942593A (zh) 一种基于干扰观测器补偿的机械臂自适应滑模控制方法
Liu et al. Laser tracker-based control for peg-in-hole assembly robot
CN116330287A (zh) 在役精度智能诊断校准的工业机器人系统及诊断校准方法
CN114571465B (zh) 一种基于模拟退火算法的四轴并联机器人校准方法
CN110569582B (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