CN112346408A - 一种构建二元一次方程组的数控磨床几何误差补偿方法 - Google Patents

一种构建二元一次方程组的数控磨床几何误差补偿方法 Download PDF

Info

Publication number
CN112346408A
CN112346408A CN202011056386.XA CN202011056386A CN112346408A CN 112346408 A CN112346408 A CN 112346408A CN 202011056386 A CN202011056386 A CN 202011056386A CN 112346408 A CN112346408 A CN 112346408A
Authority
CN
China
Prior art keywords
grinding machine
error
coordinate system
numerical control
matrix
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.)
Granted
Application number
CN202011056386.XA
Other languages
English (en)
Other versions
CN112346408B (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.)
Beijing University of Technology
Original Assignee
Beijing 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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202011056386.XA priority Critical patent/CN112346408B/zh
Publication of CN112346408A publication Critical patent/CN112346408A/zh
Application granted granted Critical
Publication of CN112346408B publication Critical patent/CN112346408B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00—Programme-control systems
    • G05B19/02—Programme-control systems electric
    • G05B19/18—Numerical 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/404—Numerical 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 control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00—Program-control systems
    • G05B2219/30—Nc systems
    • G05B2219/33—Director till display
    • G05B2219/33133—For each action define function for compensation, enter parameters

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)
  • Constituent Portions Of Griding Lathes, Driving, Sensing And Control (AREA)

Abstract

本发明公开了一种构建二元一次方程组的数控磨床几何误差补偿方法，该方法以多体系统运动学为基础，首先建立磨床拓扑结构，建立基础坐标系和刀具坐标系、工件坐标系。之后建立磨床运动部件的位置、运动矩阵及其误差矩阵，按照拓扑结构进行对应矩阵的连乘，得到有误差情况下砂轮磨削点在工件坐标系中的位置方程，将其中x方向和z方向的方程进行简化，联立得到一个关于x和z的二元一次方程组等号左侧部分，将方程组等号右侧部分赋值，带入误差参数的数值，进行求解之后得到的x和z的值即为修正后的指令坐标。本发明通过求解二元一次方程组快速得到修正的指令值，解决了以往通过迭代求解时计算繁琐、费时等缺点，便于操作，实用性好。

Description

一种构建二元一次方程组的数控磨床几何误差补偿方法
技术领域
本发明涉及机床几何误差补偿的领域，尤其是一种构建二元一次方程组的数控磨床几何误差补偿方法。
背景技术
随着机械加工领域对加工精度要求的提高，多种提高机床加工精度的方法被提出，大多数要改变机床结构，成本较高或迭代计算，繁琐费时。在控制成本、简化计算的前提下，提高加工设备的精度、减小几何误差的方法越来越重要。为了满足要求，一种构建二元一次方程组的数控磨床几何误差补偿方法被提出。
发明内容
为了解决现有几何误差补偿方法繁琐、费时的问题，本发明提出一种构建二元一次方程组的数控磨床几何误差补偿方法。该方法以多体系统运动学为基础，根据数控磨床各运动部分的结构来建立其拓扑结构，然后建立基础坐标系和刀具坐标系、工件坐标系，找到刀具坐标系、工件坐标系在基础坐标系中的相对位置。之后建立数控磨床每个运动部件的位置、运动变换矩阵及其误差矩阵，按照部件间相对结构顺序和拓扑结构进行对应矩阵的连乘，经过变换，即可得到有误差情况下砂轮磨削点在工件坐标系中的位置方程，将位置方程中x方向和z方向的方程进行简化，消去高阶误差参数，保留一阶误差参数和常量，将简化后的方程进行联立可得到一个关于x和z的二元一次方程组等号左侧部分，将方程组等号右侧部分赋予实际加工中的点位坐标值，并带入对该磨床实际测量得到各个误差参数的数值，进行求解之后得到的x和z的值即为修正磨床几何误差后的可驱动磨床进行加工的指令坐标。本发明通过求解二元一次方程组快速得到修正的指令值，解决了以往通过迭代求解时计算繁琐、费时等缺点，便于操作，实用性好。
本发明解决其技术问题所采用的技术方案是：一种构建二元一次方程组的数控磨床几何误差补偿方法，该方法包括以下步骤：
步骤一：以多体系统运动学为基础，根据数控磨床各运动部分的结构来建立其拓扑结构。
步骤二：建立磨床基础坐标系和刀具坐标系、工件坐标系，找到刀具坐标系、工件坐标系在基础坐标系中的相对位置。
步骤三：建立数控磨床每个运动部件的位置、运动变换矩阵及其误差矩阵，按照部件间相对结构顺序和拓扑结构进行对应矩阵的连乘，得到“床身-工件”分支和“床身-刀具”分支的矩阵表达式。
步骤四：矩阵表达式经过变换，得到有误差情况下砂轮磨削点在工件坐标系中的位置方程，将位置方程中x方向和z方向的方程进行简化，消去高阶误差参数，保留一阶误差参数和常量。
步骤五：将简化后的x方向和z方向的方程进行联立，得到一个关于x和z的二元一次方程组的等号左侧部分，将方程组等号右侧赋予实际加工中的工件坐标系中的点位坐标值(x和z值)。
步骤六：通过对该磨床实际测量得到各个误差参数的数值，并带入到上述二元一次方程组中。
步骤七：将该关于x和z的二元一次方程组进行求解，得到的x和z的值即为修正磨床几何误差后的可驱动磨床进行加工的指令坐标值。
步骤八：将求得的修正后的指令坐标值写入数控代码中，传入数控磨床进行加工，即可实现该数控外圆复合磨床的几何误差补偿。
本发明的有益效果是：该方法以多体系统运动学为基础，根据数控磨床各运动部分的结构来建立其拓扑结构，然后建立基础坐标系和刀具坐标系、工件坐标系，找到刀具坐标系、工件坐标系在基础坐标系中的相对位置。之后建立数控磨床每个运动部件的齐次坐标变换矩阵及其误差矩阵，按照部件间相对结构顺序和拓扑结构进行对应矩阵的连乘，经过变换，即可得到有误差情况下砂轮磨削点在工件坐标系中的位置方程，将位置方程中x方向和z方向的方程进行简化，消去高阶误差参数，保留一阶误差参数和常量，将简化后的方程进行联立可得到一个关于x和z的二元一次方程组等号左侧部分，将方程组等号右侧部分赋予实际加工中的点位坐标值，并带入对该磨床实际测量得到各个误差参数的数值，进行求解之后得到的x和z的值即为修正磨床几何误差后的可驱动磨床进行加工的指令坐标。本发明通过求解二元一次方程组快速得到修正的指令值，解决了以往通过迭代求解时计算繁琐、费时等缺点，便于操作，实用性好。
下面结合具体实施方式对本发明作详细说明。
具体实施方式
本实例采用某机床厂的数控精密外圆复合磨床。砂轮半径为R＝150mm。
构建二元一次方程组的数控外圆复合磨床几何误差补偿方法步骤如下：
步骤一：以多体系统运动学为基础，根据数控磨床各运动部分的结构来建立其拓扑结构。该磨床拓扑结构包括两分支,“床身-刀具”分支依次为：床身-x向滑座-下砖塔-上转塔-砂轮箱(砂轮)；“床身-刀具”分支依次为：床身-z向工作台-头架主轴-工件。
步骤二：首先使机床各部件移动到机床绝对零点，以此为初始条件设置工件和砂轮的体参考坐标系分别位于工件的中心O2和砂轮的中心O3处，其余部件的体参考坐标系位置均为头架主轴端面中心O1处。
步骤三：建立数控磨床每个运动部件的齐次坐标变换矩阵及其误差矩阵，按照部件间相对结构顺序和拓扑结构进行对应矩阵的连乘，得到的“床身-工件”分支和“床身-刀具”分支的矩阵表达式分别为
“床身-工件”分支的矩阵表达式为：
{Pw}0＝[S01]p[S01]pe[S01]s[S01]se[S12]p[S12]pe[S12]s[S12]se*{rw}w
“床身-刀具”分支的矩阵表达式为：
{Pt}0＝[S03]p[S03]pe[S03]s[S03]se[S34]p[S34]pe[S34]s[S34]se[S45]p[S45]pe[S45]s[S45]se*{rt}t
其中，[S01]p、[S01]pe、[S01]s、[S01]se分别为z向工作台的位置矩阵、位置误差矩阵、运动矩阵、运动误差矩阵。
[S12]p、[S12]pe、[S12]s、[S12]se分别为头架主轴的位置矩阵、位置误差矩阵、运动矩阵、运动误差矩阵。
[S03]p、[S03]pe、[S03]s、[S03]se分别为x向滑座的位置矩阵、位置误差矩阵、运动矩阵、运动误差矩阵。
[S34]p、[S34]pe、[S34]s、[S34]se分别为上转塔的位置矩阵、位置误差矩阵、运动矩阵、运动误差矩阵。
[S45]p、[S45]pe、[S45]s、[S45]se分别为砂轮的位置矩阵、位置误差矩阵、运动矩阵、运动误差矩阵。
步骤四：矩阵表达式经过变换，得到有误差情况下砂轮磨削点在工件坐标系中的位置方程,将方程简化结果为[Xw；Yw；Zw；1]，其中
Xw＝δx(x)-δx(z)-q2x+x+εz(z)*q5y+εy(z)*z+δx(b)*cos(b)+q5z*(sin(b)+εy(b)*cos(b))+δz(b)*sin(b)-εy(z)*(R*sin(b)+q5z*cos(b)-q5x*sin(b))+R*(sin(b)*(εxz+εy(x))-cos(b)+εy(b)*sin(b))-q5y*(εxb+εxy+εz(x)+εz(b)*cos(b)-εx(b)*sin(b))-q5x*(sin(b)*(εxz+εy(x))-cos(b)+εy(b)*sin(b))
Yw＝δy(b)+δy(x)-δy(z)-q2y+q5y-εyz*q2z+εxy*x-εx(z)*z-εyz*z-εz(z)*(x-R*cos(b)+q5x*cos(b)+q5z*sin(b))+(εx(z)+εyz)*(R*sin(b)+q5z*cos(b)-q5x*sin(b))-R*(εz(b)+sin(b)*(εbz+εx(x))+cos(b)*(εxb+εxy+εz(x)))-q5z*(εx(b)+cos(b)*(εbz+εx(x))-sin(b)*(εxb+εxy+εz(x)))+q5x*(εz(b)+sin(b)*(εbz+εx(x))+cos(b)*(εxb+εxy+εz(x)))
Zw＝δz(x)-δz(z)-q2z+z+εyz*q2y-εxz*x+εy(z)*x+εy(z)*(-R*cos(b)+q5x*cos(b)+q5z*sin(b))+δz(b)*cos(b)-δx(b)*sin(b)-q5x*(sin(b)+cos(b)*(εxz+εy(x))+εy(b)*cos(b))+q5y*(εbz+εx(x)+εx(b)*cos(b)+εz(b)*sin(b))-q5y*(εx(z)+Eyz)+R*(sin(b)+cos(b)*(εxz+εy(x))+εy(b)*(cos(b)-sin(b)*(εxz+εy(x))))-q5z*(sin(b)*(εxz+εy(x))-cos(b)+εy(b)*sin(b))
其中，εx(x)、εy(x)、εz(x)、δx(x)、δy(x)、δz(x)分别表示表示X向滑座沿导轨运动时产生的滚转角误差、俯仰角误差、偏摆角误差、定位误差、垂直直线度误差、水平直线度误差。
εx(z)、εy(z)、εz(z)、δx(z)、δy(z)、δz(z)分别表示表示z向工作台沿导轨运动时产生的滚转角误差、俯仰角误差、偏摆角误差、定位误差、垂直直线度误差、水平直线度误差。
εx(b)、εy(b)、εz(b)、δx(b)、δy(b)、δz(b)分别表示表示转塔B轴转动时时产生的滚转角误差、俯仰角误差、偏摆角误差、定位误差、垂直直线度误差、水平直线度误差。
q2x、q2y、q2z为工件坐标系原点在基础坐标系中的位置。
q5x、q5y、q5z、为刀具坐标系原点在基础坐标系中的位置。
步骤五：将简化后的x方向和z方向的方程进行联立，得到一个关于x和z的二元一次方程组的等号左侧部分，该部分为
Figure BDA0002710979670000041
将方程组等号右侧赋予实际加工中的工件坐标系中的点位坐标值(x和z值)，选取加工点位坐标为(70,100)，组合成的二元一次方程组为：
Figure BDA0002710979670000042
步骤六：通过对该磨床实际测量得到各个误差参数的数值，各线性误差参数的测量结果均处于0.001mm-0.002mm之间，角误差测量结果处于0.001rad-0.002rad之间，将误差参数带入到上述二元一次方程组中。
步骤七：将该关于x和z的二元一次方程组进行求解，结果x＝70.0488，z＝100.0202，得到的x和z的值即为修正磨床几何误差后的可驱动磨床进行加工的指令坐标值。
步骤八：将求得的修正后的指令坐标值写入数控代码中，传入数控磨床进行加工，磨床的几何误差因素即可被消除，实现了该数控外圆复合磨床的几何误差补偿。
通过本实施例，该数控外圆复合磨床的几何误差降低了91.32％，结合几何误差约占总误差的50％-70％，该磨床的综合精度提高了约45.5％-63.7％。

Claims (1)

1.一种构建二元一次方程组的数控磨床几何误差补偿方法，其特征在于：该方法包括以下步骤：
步骤一：以多体系统运动学为理论基础，根据数控磨床各运动部分的结构来建立其拓扑结构；
步骤二：建立磨床基础坐标系和刀具坐标系、工件坐标系，找到刀具坐标系、工件坐标系在基础坐标系中的相对位置；
步骤三：基于步骤一的拓扑结构，建立数控磨床每个运动部件的位置、运动变换矩阵及其误差矩阵，按照部件间相对结构顺序和拓扑结构进行对应矩阵的连乘，得到“床身-工件”分支和“床身-刀具”分支的矩阵表达式；
步骤四：矩阵表达式经过变换，得到有误差情况下砂轮磨削点在工件坐标系中的位置方程，将位置方程中x方向和z方向的方程进行简化，消去高阶误差参数，保留一阶误差参数和常量；
步骤五：将简化后的x方向和z方向的方程进行联立，得到一个关于x和z的二元一次方程组的等号左侧部分，将方程组等号右侧赋予实际加工中的工件坐标系中的点位坐标值即x和z值；
步骤六：通过对该磨床实际测量得到各个误差参数的数值，并带入到上述二元一次方程组中；
步骤七：将该关于x和z的二元一次方程组进行求解，得到的x和z的值即为修正磨床几何误差后的可驱动磨床进行加工的指令坐标值；
步骤八：将求得的修正后的指令坐标值写入数控代码中，传入数控磨床进行加工，即可实现该数控外圆复合磨床的几何误差补偿。
CN202011056386.XA 2020-09-30 2020-09-30 一种构建二元一次方程组的数控磨床几何误差补偿方法 Active CN112346408B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011056386.XA CN112346408B (zh) 2020-09-30 2020-09-30 一种构建二元一次方程组的数控磨床几何误差补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011056386.XA CN112346408B (zh) 2020-09-30 2020-09-30 一种构建二元一次方程组的数控磨床几何误差补偿方法

Publications (2)

Publication Number Publication Date
CN112346408A true CN112346408A (zh) 2021-02-09
CN112346408B CN112346408B (zh) 2023-11-24

Family

ID=74361431

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011056386.XA Active CN112346408B (zh) 2020-09-30 2020-09-30 一种构建二元一次方程组的数控磨床几何误差补偿方法

Country Status (1)

Country Link
CN (1) CN112346408B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112949098A (zh) * 2021-04-19 2021-06-11 广东工业大学 运动学误差映射矩阵的迭代修正方法及其迭代修正系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107966957A (zh) * 2017-11-27 2018-04-27 重庆大学 一种数控蜗杆砂轮磨齿机空间误差解耦补偿方法
CN109471408A (zh) * 2018-11-01 2019-03-15 西北工业大学 基于nc代码重构的三轴机床几何误差补偿方法
CN109732402A (zh) * 2019-03-14 2019-05-10 西安交通大学 基于激光干涉仪的多线机床空间几何误差测量辨识方法
CN110109418A (zh) * 2019-05-19 2019-08-09 重庆理工大学 一种大型龙门五面加工中心的几何误差快速辨识方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107966957A (zh) * 2017-11-27 2018-04-27 重庆大学 一种数控蜗杆砂轮磨齿机空间误差解耦补偿方法
CN109471408A (zh) * 2018-11-01 2019-03-15 西北工业大学 基于nc代码重构的三轴机床几何误差补偿方法
CN109732402A (zh) * 2019-03-14 2019-05-10 西安交通大学 基于激光干涉仪的多线机床空间几何误差测量辨识方法
CN110109418A (zh) * 2019-05-19 2019-08-09 重庆理工大学 一种大型龙门五面加工中心的几何误差快速辨识方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
熊成: "机床旋转轴几何误差测量与补偿技术研究", 《工程科技I辑》, no. 3 *
范晋伟等: "基于多体系统理论的汽车凸轮轴磨削几何误差建模与辨识技术理论研究", 《中国机械工程》, vol. 24, no. 8, pages 2216 - 2222 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112949098A (zh) * 2021-04-19 2021-06-11 广东工业大学 运动学误差映射矩阵的迭代修正方法及其迭代修正系统

Also Published As

Publication number Publication date
CN112346408B (zh) 2023-11-24

Similar Documents

Publication Publication Date Title
CN109471408B (zh) 基于nc代码重构的三轴机床几何误差补偿方法
CN107336023A (zh) 一种增减材复合加工中心
CN102069222B (zh) 一种双刀加工叶片的卧式机床结构
CN107966957B (zh) 一种数控蜗杆砂轮磨齿机空间误差解耦补偿方法
CN201685166U (zh) 一种复合数控磨削中心
CN108829038B (zh) 一种刀尖运动轨迹控制算法
CN205437784U (zh) 一种自动钻铣系统及钻铣生产线
CN102069223B (zh) 一种双刀加工叶片的立式机床结构
US10146204B2 (en) Machining method and a control device for a machine tool
CN112558549B (zh) 一种面向大部件群孔加工孔位误差最小的基准选取方法
CN111702762A (zh) 一种工业机器人作业姿态优化方法
CN112139513B (zh) 基于机床的b轴单点车削自由曲面的加工方法
CN102621929B (zh) 双转台五轴联动数控机床的夹具高度及加工路径优化方法
CN112346408B (zh) 一种构建二元一次方程组的数控磨床几何误差补偿方法
CN107350748A (zh) 一种v型机机体斜面及斜面孔系的加工方法
CN113579766B (zh) 一种六自由度串并混联数控机床及其后处理方法
CN208068016U (zh) 转动机构及其磨床
CN207087324U (zh) 一种增减材复合加工中心
CN219649210U (zh) 一种六轴数控机床
CN106774163B (zh) 高精度三维轮廓控制方法及装置
CN108549323B (zh) 一种五轴数控工具磨床通用后置处理方法
CN102554379A (zh) 一种超硬刀具电火花加工装置及操作方法
CN115609088A (zh) 齿轮及其齿面加工反调修正方法
CN108453622A (zh) 转动机构及其在机加工中的应用
CN107861467A (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