CN113704933A - 一种基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法 - Google Patents
一种基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法 Download PDFInfo
- Publication number
- CN113704933A CN113704933A CN202110810291.0A CN202110810291A CN113704933A CN 113704933 A CN113704933 A CN 113704933A CN 202110810291 A CN202110810291 A CN 202110810291A CN 113704933 A CN113704933 A CN 113704933A
- Authority
- CN
- China
- Prior art keywords
- error
- differential motion
- numerical control
- coordinate system
- 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.)
- Pending
Links
- 238000003754 machining Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 19
- 239000011159 matrix material Substances 0.000 claims abstract description 57
- 239000013598 vector Substances 0.000 claims abstract description 41
- 230000009466 transformation Effects 0.000 claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 230000008859 change Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Computer Networks & Wireless Communication (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法,该方法以相邻体间坐标系间的微分运动理论为基础,结合机器人正向运动学理论和以往多体系统的低序体理论,搭建数控外圆磨床正向运动拓扑结构;得到各个运动部件相邻体间的齐次变换矩阵;通过此相邻体间的齐次矩阵得到刀具相对于各个运动部件坐标系的齐次变换矩阵;通过对数控外圆磨床进行几何误差项分析,得到各个运动轴的几何误差的微分运动矢量,再由此得到平动轴与旋转轴的几何误差微分运动矢量表达式。本发明解决了以往的空间误差模型计算复杂问题,具有计算简便快速、应用灵活和创新的优点。
Description
技术领域
本申请涉及机床几何误差建模的领域,尤其是一种基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法。
背景技术
对数控机床进行精度分析,机床误差建模是基础工作,基于多体系统理论搭建几何误差数学模型需要建立多个坐标系,得到的空间加工误差矩阵计算量过程较复杂,且由于二阶及以上的项式对综合几何误差影响较小,经过MATLAB软件计算后,通常需要进行手动消除二阶及以上的项式,才能得到机床空间加工误差模型,较为复杂和麻烦。况且所建立的空间几何误差模型无法反映单个运动部件对于刀具所产生的误差影响,也不能与机床加工代码坐标系一致,方便直观地观察机床加工误差。为了满足要求,一种基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法被提出。
发明内容
为了解决现存的几何误差建模复杂且描述其他运动部件不全面的问题,本发明提出一种基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法。该方法以坐标系间的微分运动关系理论为基础,首先定义各运动轴的坐标系,结合机器人正向运动学理论和以往多体系统的低序体理论,搭建数控外圆磨床正向运动拓扑结构;得到各个运动部件相邻体间的齐次变换矩阵;通过此相邻体间的齐次矩阵得到刀具相对于各个运动部件坐标系的齐次变换矩阵;通过对数控外圆磨床进行几何误差项分析,得到各个运动轴的几何误差的微分运动矢量,再由此得到平动轴与旋转轴的几何误差微分运动矢量表达式;经MATLAB R2018a计算得到各个运动部件相对于刀具坐标系的微分运动矩阵,此过程的计算反映出各个运动部件对数控外圆磨床加工精度的影响;将这些运动部件对磨床的影响叠加,得到刀具坐标系下机床的综合空间加工误差模型;为与机床加工代码坐标系一致方便观察机床加工误差,综合空间加工误差模型应建立在工件坐标系下,得到刀具在工件坐标系下的位置误差和姿态误差。本发明通过微分运动矩阵得到各部件对刀具的影响、刀具坐标系下综合加工误差模型和工件坐标系下的刀具位置误差与姿态误差,解决了以往的空间误差模型计算复杂问题,无法反映各部件对刀具的影响程度问题及可方便的观察到机床加工误差(刀具相对于工件的位置误差与姿态误差),具有计算简便快速、应用灵活和创新的优点。
本发明解决其技术问题所采用的技术方案是:一种基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法,该方法包括以下步骤:
步骤一:依据国际惯例,确定机床各运动轴和坐标系;结合机器人正向运动学理论和多体系统-低序体理论,建立数控外圆磨床工件—砂轮刀具的正向运动拓扑结构。
步骤二:基于坐标系间的微分运动关系理论,得到数控外圆磨床的各个运动部件相邻体间的齐次变换矩阵。
步骤三:依据工件—砂轮刀具的正向运动拓扑结构和各相邻体间的齐次变换矩阵,分别得到刀具相对于各运动部件的齐次变换矩阵。
步骤四:依据刚体六个自由度运动学理论,分析数控外圆磨床各轴的几何误差项,并确定平动轴和回转轴的各项几何误差的微分运动矢量。
步骤五:结合上述的各运动轴几何误差项的微分运动矢量及各部件在刀具坐标系下的齐次变换矩阵,经MATLAB计算得到各运动部件相对于刀具的微分运动矩阵,此矩阵即反映各个运动部件对于数控外圆磨床的刀具加工精度的影响。
步骤六:将各轴的几何误差项的微分运动矢量与对应的微分运动矩阵相乘,便得到各轴的几何误差项在刀具坐标系下的矢量形式,将这些矢量叠加后得到刀具坐标系下的机床综合空间加工误差模型。
步骤七:将此综合空间加工误差模型转换到工件坐标系上。首先需结合刀具坐标系下的综合几何误差矢量,得到它的误差变化矩阵,再根据数控外圆磨床的正向运动拓扑结构得到相对于工件的综合空间误差模型。
步骤八:分别将工件坐标系下的综合空间误差模型与[0 0 0 1]T、[0 0 1 0]T相乘,便得到工件坐标系下的刀具位置误差和姿态误差。
本发明的有益效果是:该方法以坐标系间的微分运动理论为基础,依据机器人正向运动学理论和多体系统-低序体理论搭建数控外圆磨床的工件—砂轮刀具的正向运动拓扑结构;再通过坐标系间的微分运动关系理论,得到数控外圆磨床的各个运动部件相邻体间的齐次变换矩阵;再依据工件—砂轮刀具的正向运动拓扑结构和各相邻体间的齐次变换矩阵,分别得到刀具相对于各运动部件的齐次变换矩阵;然后,由刚体的六个自由度运动学理论分析数控外圆磨床各轴的几何误差项,并确定平动轴和回转轴的各项几何误差的微分运动矢量;之后结合每一个运动轴几何误差项的微分运动矢量和各自在刀具坐标系下的齐次变换矩阵,经MATLAB计算得到各运动部件相对于刀具的微分运动矩阵,将各轴的几何误差项的微分运动矢量与对应的微分运动矩阵相乘,便得到各轴的几何误差项在刀具坐标系下的矢量形式,这些矢量叠加后就是刀具坐标系下的机床综合空间加工误差模型。若要将此综合空间加工误差模型转换到工件坐标系上,直观反映机床加工精度。首先需结合刀具坐标系下的综合几何误差矢量,得到它的误差变化矩阵,再根据数控外圆磨床的正向运动拓扑结构得到相对于工件的综合空间误差模型。分别将工件坐标系下的综合空间误差模型与[0 0 0 1]T、[0 0 1 0]T相乘,便得到工件坐标系下的刀具位置误差和姿态误差。本发明基于微分运动矩阵得到各部件对刀具的影响、刀具坐标系下综合加工误差模型和工件坐标系下的刀具位置误差与姿态误差,解决了以往的空间误差模型计算复杂问题,无法反映各部件对刀具的影响程度问题及可方便的观察到机床加工误差,具有计算简便快速、应用灵活和创新的优点。
下面结合具体实施方式对本发明作详细说明。
具体实施方式
本实例采用某公司的数控外圆磨床。主轴C和转台B的旋转角度分别为c、b,滑座Z上的导轨和滑座X上的导轨移动距离分别为z、x。
基于微分运动矩阵的数控外圆磨床几何误差综合空间加工误差建模方法步骤如下:
步骤一:确定机床各运动轴和坐标系;结合机器人正向运动学理论和多体系统-低序体理论,建立此数控外圆磨床正向运动拓扑结构。该拓扑结构为:工件(初始端)—主轴—Z向滑座—床身—X向滑座—砂轮转台—砂轮刀具(末端执行器)。
步骤三:依据工件—砂轮刀具的正向运动拓扑结构和各相邻体间的齐次变换矩阵,分别得到刀具相对于的齐次变换矩阵。
刀具相对于工件的齐次变换矩阵为:
刀具相对于主轴(回转轴C)的齐次变换矩阵为:
刀具相对于滑座(平动轴Z)的齐次变换矩阵为:
步骤四:依据刚体六个自由度运动学理论,分析数控外圆磨床各轴的几何误差项,并确定平动轴和回转轴的各项几何误差的微分运动矢量。
运动误差也被称为与位置相关的误差,随着运动轴移动和旋转的变换而变化,这里我们将它们称为基本几何误差项,以平动轴Z为例,六项基本几何误差项为:δx(z)、δy(z)、δz(z)、εx(z)、εy(z)、εz(z),这六项基本误差是与z(平动轴Z的移动量/Z方向的加工代码指令)有关的非线性函数。
位置误差是恒定值,与运动轴的运动量无关,包含垂直度误差和运动轴位置偏差。
各运动轴的六项基本几何误差项的微分运动矢量为:
Berror_i=[δx(i) δy(i) δz(i) εx(i) εy(i) εz(i)]T。前三项表示坐标轴i方向上的线位移误差,后三项表示坐标轴i方向上的角位移误差,i表示各运动轴X、Z、B、C。
垂直度误差的微分运动矢量为:Verror_z=[0 0 0 0 vxz 0]T;旋转轴的位置偏差的微分运动矢量为:Perror_b=[pbx 0 pbz 0 0 0]T、Perror_c=[pcx 0 0 0 0 0]T
步骤五:结合上述的各运动轴几何误差项的微分运动矢量及刀具相对于各部件的齐次变换矩阵,得到各运动部件相对于刀具的微分运动矩阵,此矩阵即反映各个运动部件对于数控外圆磨床的刀具加工精度的影响。
结合上述的微分运动矢量,得到平动轴和回转轴的几何误差微分运动矢量分别为:
error_x=Berror_x
=[δx(x) δy(x) δz(x) εx(x) εy(x) εz(x)]T
error_z=Berror_z+Verror_z
=[δx(z) δy(z) δz(z) εx(z) εy(z) εz(z)]T+[01×3 vxz 0]T
=[δx(z) δy(z) δz(z) εx(z) εy(z)+vxz εz(z)]T
error_c=Berror_B+Perror_B
=[δx(b) δy(b) δz(b) εx(b) εy(b) εz(b)]T+[0 pbz 01×3]T
=[δx(b)+pbx δy(b) δz(b)+pbz εx(b) εy(b) εz(b)]T
error_c=Berror_c+Perror_c
=[δx(c) δy(c) δz(c) εx(c) εy(c) εz(c)]T+[pcx 01×5]T
=[δx(c)+pcx δy(c) δz(c) εx(c) εy(c) εz(c)]T
而工件相对于刀具坐标系的微分运动矩阵为:
回转轴C相对于刀具坐标系的微分运动矩阵分别为:
其中,
步骤六:将各轴的几何误差项的微分运动矢量与对应的微分运动矩阵相乘,便得到各运动部件的几何误差项在刀具坐标系下的矢量形式,将这些矢量叠加后得到刀具坐标系下的机床综合空间加工误差模型。
平动轴X的几何误差项在刀具坐标系的矢量形式为:
平动轴Z的几何误差项在刀具坐标系的矢量形式为:
回转轴B的几何误差项在刀具坐标系的矢量形式为:
回转轴C的几何误差项在刀具坐标系的矢量形式为:
床身在刀具坐标系的矢量形式为:
工件在刀具坐标系的矢量形式为:
则机床综合空间加工误差模型为:
步骤七:将此综合空间加工误差模型转换到工件坐标系上。首先需结合刀具坐标系下的综合几何误差矢量,得到它的误差变化矩阵,再根据数控外圆磨床的正向运动拓扑结构得到相对于工件的综合空间误差模型。
刀具的综合误差矢量表示为:
error_t=[δx(t) δy(t) δz(t) εx(t) εy(t) εz(t)]T
可得到刀具坐标系下,刀具的误差矩阵为:
机床综合几何误差矩阵是实际齐次变换矩阵与理想齐次变换矩阵的差值,该数控外圆磨床基于微分运动矩阵的空间加工几何误差模型在工件坐标系下表示为:
步骤八:分别将工件坐标系下的综合空间误差模型与[0 0 0 1]T、[0 0 1 0]T相乘,便得到工件坐标系下的刀具位置误差和姿态误差。
则工件坐标系下的刀具位置误差为:
其中,刀具位置误差的三个方向分量分别为:
工件坐标系下的刀具姿态误差为:
其中,刀具姿态误差的分量分别为:
Tattz=-sin(b)*(εy(b)+εy(x)+εy(z)+vxz+δy(c)*cos(c)+εx(z)*sin(c))
通过本实施例,该数控精密外圆磨床的模型搭建过程迅速、简便且计算量小、各微分运动矩阵物理意义明确,可反映各运动部件对刀具的影响。
Claims (2)
1.一种搭建基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法,其特征在于,包括以下步骤:
步骤一:确定机床各运动轴和坐标系;结合机器人正向运动学理论和多体系统-低序体理论,建立数控外圆磨床工件—砂轮刀具的正向运动拓扑结构;
步骤二:基于坐标系间的微分运动关系理论,得到数控外圆磨床的各个运动部件相邻体间的齐次变换矩阵;
步骤三:依据工件—砂轮刀具的正向运动拓扑结构和各相邻体间的齐次变换矩阵,分别得到刀具相对于各运动部件的齐次变换矩阵;
步骤四:依据刚体六个自由度运动学理论,分析数控外圆磨床各轴的几何误差项,并确定平动轴和回转轴的各项几何误差的微分运动矢量;
步骤五:结合上述的各运动轴几何误差项的微分运动矢量及各部件在刀具坐标系下的齐次变换矩阵,经MATLAB计算得到各运动部件相对于刀具的微分运动矩阵;
步骤六:将各轴的几何误差项的微分运动矢量与对应的微分运动矩阵相乘,便得到各轴的几何误差项在刀具坐标系下的矢量形式,将这些矢量叠加后得到刀具坐标系下的机床综合空间加工误差模型;
步骤七:将此综合空间加工误差模型转换到工件坐标系上;首先需结合刀具坐标系下的综合几何误差矢量,得到它的误差变化矩阵,再根据数控外圆磨床的正向运动拓扑结构得到相对于工件的综合空间误差模型;
步骤八:分别将工件坐标系下的综合空间误差模型与[0 0 0 1]T、[0 0 1 0]T相乘,便得到工件坐标系下的刀具位置误差和姿态误差。
2.根据权利要求1所述的一种搭建基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法,其特征在于,步骤五中,微分运动矩阵反映各个运动部件对于数控外圆磨床的刀具加工精度的影响。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110810291.0A CN113704933A (zh) | 2021-07-19 | 2021-07-19 | 一种基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110810291.0A CN113704933A (zh) | 2021-07-19 | 2021-07-19 | 一种基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113704933A true CN113704933A (zh) | 2021-11-26 |
Family
ID=78648838
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110810291.0A Pending CN113704933A (zh) | 2021-07-19 | 2021-07-19 | 一种基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113704933A (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106502203A (zh) * | 2016-10-08 | 2017-03-15 | 西南交通大学 | 一种数控机床几何误差建模方法 |
CN109613889A (zh) * | 2019-01-07 | 2019-04-12 | 安徽理工大学 | 基于微分变换的数控机床在机测量系统综合误差补偿方法 |
CN110532667A (zh) * | 2019-08-26 | 2019-12-03 | 大连理工大学 | 一种用于提高精密机床空间几何误差模型建模精度的方法 |
-
2021
- 2021-07-19 CN CN202110810291.0A patent/CN113704933A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106502203A (zh) * | 2016-10-08 | 2017-03-15 | 西南交通大学 | 一种数控机床几何误差建模方法 |
CN109613889A (zh) * | 2019-01-07 | 2019-04-12 | 安徽理工大学 | 基于微分变换的数控机床在机测量系统综合误差补偿方法 |
CN110532667A (zh) * | 2019-08-26 | 2019-12-03 | 大连理工大学 | 一种用于提高精密机床空间几何误差模型建模精度的方法 |
Non-Patent Citations (1)
Title |
---|
付国强: "基于指数积理论和坐标系微分运动关系的数控机床几何误差建模与补偿方法研究", 中国博士学位论文全文数据库(工程科技Ⅰ辑), no. 12, pages 022 - 87 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fu et al. | Product-of-exponential formulas for precision enhancement of five-axis machine tools via geometric error modeling and compensation | |
Abbaszadeh-Mir et al. | Theory and simulation for the identification of the link geometric errors for a five-axis machine tool using a telescoping magnetic ball-bar | |
KR100271854B1 (ko) | 다축기계의 3차원입체오차측정방법 | |
CN107450473A (zh) | 一种cfxyza型五轴数控机床旋转轴几何误差计算、补偿及其验证方法 | |
CN108829038B (zh) | 一种刀尖运动轨迹控制算法 | |
Yang et al. | A generalized online estimation algorithm of multi-axis contouring errors for CNC machine tools with rotary axes | |
Li et al. | A sensitivity method to analyze the volumetric error of five-axis machine tool | |
CN102091967A (zh) | 一种多轴数控加工的进给速度平滑方法 | |
CN109375579A (zh) | 一种基于运动学的五轴数控加工刀具姿态规划方法 | |
CN109933920B (zh) | 一种旋转轴位置偏差的误差矢量建模方法 | |
CN113325802B (zh) | 一种五轴机床几何误差补偿方法 | |
Fan et al. | Research on geometric error modeling and compensation method of CNC precision cylindrical grinding machine based on differential motion theory and Jacobian matrix | |
CN116909209A (zh) | 一种考虑动态热误差的数控机床误差建模和预测方法 | |
CN113704933A (zh) | 一种基于坐标系间微分运动关系的数控外圆磨床综合空间加工误差建模方法 | |
Yun et al. | A geometric postprocessing method for 5-axis machine tools using locations of joint points | |
Holub et al. | Application on-the-fly measurement of CNC machine tools | |
She et al. | Development of multi-axis numerical control program for mill—turn machine | |
Sung et al. | Modeling/analysis of four-half axis machine tool via modified denavit-hartenberg notation | |
CN108776459B (zh) | 一种提升五轴数控机床加工精度的工艺方法 | |
CN109933918A (zh) | 旋转轴垂直度误差的误差矢量建模方法 | |
CN113156888B (zh) | 回转摆头型五轴可重构机床敏感几何误差高效计算方法 | |
CN109521729A (zh) | 一种三轴数控机床可补偿及不可补偿几何误差分离方法 | |
Fan et al. | Error modeling and sensitivity analysis of turning-milling compound CNC machine tools | |
Xing et al. | Accuracy enhancement in manufacture of spiral bevel gear with multi-axis CNC machine tools by a new compensation method | |
Chun et al. | Study and analysis of the 3-PRS parallel mechanism |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20211126 |