CN105058163B - 磨齿机工作状态时热误差测量装置及测量方法 - Google Patents
磨齿机工作状态时热误差测量装置及测量方法 Download PDFInfo
- Publication number
- CN105058163B CN105058163B CN201510236553.1A CN201510236553A CN105058163B CN 105058163 B CN105058163 B CN 105058163B CN 201510236553 A CN201510236553 A CN 201510236553A CN 105058163 B CN105058163 B CN 105058163B
- Authority
- CN
- China
- Prior art keywords
- mtd
- mrow
- msub
- mtr
- displacement transducer
- 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.)
- Expired - Fee Related
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q17/00—Arrangements for observing, indicating or measuring on machine tools
- B23Q17/006—Arrangements for observing, indicating or measuring on machine tools for indicating the presence of a work or tool in its holder
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Automatic Control Of Machine Tools (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种磨齿机工作状态时热误差测量装置,包括温度测量部分和位移测量部分;温度测量部分为在磨齿机的主要发热部位上设置温度传感器;位移测量部分包括通过支架固定在磨齿机主轴基座上的X向位移传感器和Y向位移传感器,X向位移传感器和Y向位移传感器与检棒接触,检棒固定在转台上。还公开了利用该装置进行热误差测量的方法,包括位移测量、计算检棒圆心坐标序列、排除***误差、计算转台轴线表达式以及热误差建模。该装置在现有磨齿机上设置X向位移传感器、Y向位移传感器以及在磨齿机的主要发热部位布置温度传感器,用该装置进行热误差测量,测量精度高,补偿效果更好。该装置结构简单,测量方法容易掌握,易于推广。
Description
技术领域
本发明属于精密机床加工技术领域,涉及一种磨齿机工作状态时热误差测量装置,本发明还涉及利用上述装置测量磨齿机热误差的方法。
背景技术
在磨齿加工中,冷却液在冷却工件的同时吸收大量磨削热,集聚在工作台的T型槽和导流槽中,造成机床关键部件和加工环境的明显升温,这种情况称为冷却液升温效应。冷却液升温效应是导致磨齿机热误差的主要因素。
目前国内外针对机床热误差的研究主要在非加工状态下进行,没有充分考虑冷却液升温效应对于机床热误差的影响。真实加工条件下,飞溅的冷却液会对激光干涉仪和电涡流传感器造成严重干扰;而不考虑冷却液升温效应的非加工状态热误差预测模型对实际加工状态下的热误差预测效果很差。因此,进行加工状态下热误差测量对提高精密磨齿机加工精度具有重要意义。
发明内容
本发明的目的是提供一种磨齿机工作状态时热误差测量装置,解决了现有装置因无法测量磨齿机工作状态下的热误差而造成的磨齿机加工精度低的问题。
本发明的另一目的是提供利用上述测量装置测量磨齿机热误差的方法。
本发明所采用的第一种技术方案是,一种磨齿机工作状态时热误差测量装置,包括温度测量部分和位移测量部分;温度测量部分包括设置在磨齿机的主轴电机、主轴前后端盖、进给驱动电机、转台面、冷却液导流槽、立柱、床身、转台基座、主轴基座上的温度传感器;位移测量部分包括通过支架固定在磨齿机主轴基座上的X向位移传感器和Y向位移传感器,X向位移传感器和Y向位移传感器与检棒接触,检棒固定在转台上。
本发明第一种技术方案的特点还在于,
检棒长度为680mm-720mm,直径为48mm-52mm。
支架为碳纤维支架。
本发明所采用的第二种技术方案是,磨齿机工作状态时热误差测量方法,采用一种测量装置,其结构为:
包括温度测量部分和位移测量部分;温度测量部分包括设置在磨齿机的主轴电机、主轴前后端盖、进给驱动电机、转台面、冷却液导流槽、立柱、床身、转台基座、主轴基座上的温度传感器;位移测量部分包括通过支架固定在磨齿机主轴基座上的X向位移传感器和Y向位移传感器,X向位移传感器和Y向位移传感器与检棒接触,检棒固定在转台上;
检棒长度为680mm-720mm,直径为48mm-52mm;
支架为碳纤维支架;
测量方法具体包括以下步骤:
步骤1:位移测量:
首先设X向位移传感器轴线与Y向位移传感器轴线交点为测量坐标系原点O,测量坐标系X、Y、Z轴平行于机床坐标系;然后旋转转台,旋转的过程中记录X向位移传感器的数据Dxij和Y向位移传感器的数据Dyij;
步骤2:计算检棒圆心坐标序列:
依据步骤1得到的X向位移传感器的数据Dxij和Y向位移传感器的数据Dyij以及检棒半径rb,利用平面圆方程(x-x0)2+(y-y0)2=rb 2计算出当前测量平面下检棒圆心坐标序列rj,j=(1,2,3,4,5,6);
步骤3:排除***误差:
将步骤2得到的检棒圆心坐标序列rj进行最小二乘圆拟合,得到转台轴心坐标序列Ri,i=(1,2,3,4,5);
步骤4:求转台轴线表达式:
将步骤3得到的转台轴心坐标序列Ri进行空间直线最小二乘拟合,求出转台轴线空间表达式为:
式(1)中,M(Rx0,Ry0,Rz0)为转台轴线上一点,a=(X,Y,Z)为转台轴线方向向量;
步骤5:热误差建模,具体方法为:
步骤5.1:根据步骤3得出的转台轴心坐标序列Ri,求出测量坐标系X、Y方向分量Xi和Yi,再根据所述步骤4得出的转台轴线表达式导出转台轴线角度热误差序列:θp=(θ1,θ2,...θn),其中p为测量数据组数,p=(1,2,……,n);
步骤5.2:提取温度传感器数据Tpq,q为温度传感器编号,q=(1,2,3,4,5,6,7,8,9,10),p为测量数据组数,p=(1,2,……,n);
步骤5.3:根据相关性判别公式:
式(2)中,Rp为某一界面内检棒圆心坐标随着温度变化的热误差序列,
令:其中,λ为相关性系数阀值,
利用多元线性回归理论建立热误差模型:
上式中,A、B、C分别代表热误差序列X、Y、θ的热误差模型中的温度变量系数向量;
步骤5.4:根据最小二乘原理求解系数向量A=(TTT)-1TTX,B=(TTT)-1TTY,C=(TTT)-1TTθ,求得系数向量A、B、C;
步骤5.5:将所述步骤5.4得出的系数向量A、B、C代入热误差模型中即得到磨齿机的热误差。
本发明第二种技术方案的特点还在于,
步骤1具体按照以下步骤实施:
步骤1.1:沿坐标系Z方向每隔90mm-110mm标记5个与坐标系Z轴垂直的测量平面,记为Si,其中,i=(1,2,3,4,5);
步骤1.2:将X向位移传感器和Y向位移传感器置于步骤1.1所述的测量平面Si内,令转台旋转360°,期间每旋转60°记录X向位移传感器的数据Dxij和Y向位移传感器的数据Dyij,其中,i=(1,2,3,4,5),j=(1,2,3,4,5,6);
步骤1.3:间隔10-15分钟执行步骤1.2一次,完成p组测量数据采集。
步骤4求转台轴线表达式的具体方法为:
设di为转台轴心坐标序列Ri到轴线距离的平方和,
则:
令:
求解后得到转台轴线上一点M(Rx0,Ry0,Rz0)以及转台轴线方向向量a=(X,Y,Z);
将M(Rx0,Ry0,Rz0)和a=(X,Y,Z)代入中得出转台轴线表达式。
本发明的有益效果是:磨齿机工作状态时热误差测量装置,在现有磨齿机上设置X向位移传感器和Y向位移传感器,可以对加工状态的磨齿机进行位移测量,在磨齿机的主要发热部位布置温度传感器,可以采集磨齿机工作状态时的温度,将位移与温度结合便可得出磨齿机工作状态时的热误差;用该装置进行热误差测量,测量精度高,补偿效果更好;并且利用磨齿加工Z方向往复运动进行测量,避免引入X、Y方向重复定位误差;该装置结构简单,便于现场安装,测量方法容易掌握,易于推广。
附图说明
图1是本发明一种磨齿机工作状态时热误差测量装置的结构示意图;
图2是利用一种磨齿机工作状态时热误差测量装置测量热误差方法中排除***误差的示意图。
图中,1.支架,2.X向位移传感器,3.Y向位移传感器,4.检棒,5.转台,6.主轴基座。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明一种磨齿机工作状态时热误差测量装置的结构如图1所示,包括温度测量部分和位移测量部分;温度测量部分包括设置在磨齿机的主轴电机、主轴前后端盖、进给驱动电机、转台面、冷却液导流槽、立柱、床身、转台基座、主轴基座上的温度传感器;位移测量部分包括通过支架1固定在磨齿机主轴基座6上的X向位移传感器2和Y向位移传感器3,X向位移传感器2和Y向位移传感器3与检棒4接触,检棒4固定在转台5上,检棒4长度为680mm-720mm,直径为48mm-52mm,支架1为碳纤维支架。
磨齿机工作状态时热误差测量方法,采用一种测量装置,其结构为:
包括温度测量部分和位移测量部分;温度测量部分包括设置在磨齿机的主轴电机、主轴前后端盖、进给驱动电机、转台面、冷却液导流槽、立柱、床身、转台基座、主轴基座上的温度传感器;位移测量部分包括通过支架1固定在磨齿机主轴基座6上的X向位移传感器2和Y向位移传感器3,X向位移传感器2和Y向位移传感器3与检棒4接触,检棒4固定在转台5上,检棒4长度为680mm-720mm,直径为48mm-52mm,支架1为碳纤维支架。
测量方法具体包括以下步骤:
步骤1:位移测量:
首先设X向位移传感器2轴线与Y向位移传感器3轴线交点为测量坐标系原点O,测量坐标系X、Y、Z轴平行于机床坐标系;然后旋转转台5,旋转的过程中记录X向位移传感器2的数据Dxij和Y向位移传感器3的数据Dyij;
步骤1具体按照以下步骤实施,如图2所示:
步骤1.1:沿坐标系Z方向每隔90mm-110mm标记5个与坐标系Z轴垂直的测量平面,记为Si,其中,i=(1,2,3,4,5);
步骤1.2:将X向位移传感器2和Y向位移传感器3置于步骤1.1所述的测量平面Si内,令转台5旋转360°,期间每旋转60°记录X向位移传感器2的数据Dxij和Y向位移传感器3的数据Dyij,其中,i=(1,2,3,4,5),j=(1,2,3,4,5,6);
步骤1.3:间隔10-15分钟执行步骤1.2一次,完成p组测量数据采集;
步骤2:计算检棒4圆心坐标序列:
依据步骤1得到的X向位移传感器2的数据Dxij和Y向位移传感器3的数据Dyij以及检棒4半径rb,利用平面圆方程(x-x0)2+(y-y0)2=rb 2计算出当前测量平面下检棒4圆心坐标序列rj,j=(1,2,3,4,5,6);
步骤3:排除***误差:
将步骤2得到的检棒4圆心坐标序列rj进行最小二乘圆拟合,得到转台5轴心坐标序列Ri,i=(1,2,3,4,5);
步骤4:求转台5轴线表达式:
将步骤3得到的转台5轴心坐标序列Ri进行空间直线最小二乘拟合,求出转台5轴线空间表达式为:
式(1)中,M(Rx0,Ry0,Rz0)为转台轴线上一点,a=(X,Y,Z)为转台轴线方向向量;
具体方法为:
设di为转台轴心坐标序列Ri到轴线距离的平方和,
则:
令:
求解后得到转台轴线上一点M(Rx0,Ry0,Rz0)以及转台轴线方向向量a=(X,Y,Z);
将M(Rx0,Ry0,Rz0)和a=(X,Y,Z)代入中得出转台轴线表达式;
步骤5:热误差建模,具体方法为:
步骤5.1:根据步骤3得出的转台5轴心坐标序列Ri,求出测量坐标系X、Y方向分量Xi和Yi,再根据所述步骤4得出的转台轴线表达式导出转台5轴线角度热误差序列:θp=(θ1,θ2,...θn),其中p为测量数据组数,p=(1,2,……,n);
步骤5.2:提取温度传感器数据Tpq,q为温度传感器编号,q=(1,2,3,4,5,6,7,8,9,10),p为测量数据组数,p=(1,2,……,n);
步骤5.3:根据相关性判别公式:
式(2)中,Rp为某一界面内检棒圆心坐标随着温度变化的热误差序列,
令:其中,λ为相关性系数阀值,λ的范围为0.86到0.9,
利用多元线性回归理论建立热误差模型:
上式中,A、B、C分别代表热误差序列X、Y、θ的热误差模型中的温度变量系数向量;
步骤5.4:根据最小二乘原理求解系数向量A=(TTT)-1TTX,B=(TTT)-1TTY,C=(TTT)-1TTθ,求得系数向量A、B、C;
步骤5.5:将步骤5.4得出的系数向量A、B、C代入热误差模型中即得到磨齿机的热误差。
本发明磨齿机工作状态时热误差测量装置,在现有磨齿机上设置X向位移传感器和Y向位移传感器,可以对加工状态的磨齿机进行位移测量,在磨齿机的主要发热部位布置温度传感器,可以采集磨齿机工作状态时的温度,将位移与温度结合便可得出磨齿机工作状态时的热误差;用该装置进行热误差测量,测量精度高,补偿效果更好;并且利用磨齿加工Z方向往复运动进行测量,避免引入X、Y方向重复定位误差;该装置结构简单,便于现场安装,测量方法容易掌握,易于推广。
Claims (3)
1.磨齿机工作状态时热误差测量方法,其特征在于,采用一种测量装置,其结构为:
包括温度测量部分和位移测量部分;
所述温度测量部分包括分别设置在磨齿机的主轴电机、主轴前后端盖、进给驱动电机、转台面、冷却液导流槽、立柱、床身、转台基座、主轴基座上的温度传感器;
所述位移测量部分包括通过支架(1)固定在磨齿机主轴基座(6)上的X向位移传感器(2)和Y向位移传感器(3),所述X向位移传感器(2)和Y向位移传感器(3)与检棒(4)接触,所述检棒(4)固定在转台(5)上;
所述检棒(4)长度为680mm-720mm,直径为48mm-52mm;
所述支架(1)为碳纤维支架;
测量方法具体包括以下步骤:
步骤1:位移测量:
首先设X向位移传感器(2)轴线与Y向位移传感器(3)轴线交点为测量坐标系原点O,测量坐标系X、Y、Z轴平行于机床坐标系;然后旋转转台(5),旋转的过程中记录X向位移传感器(2)的数据Dxij和Y向位移传感器(3)的数据Dyij;
步骤2:计算检棒(4)圆心坐标序列:
依据所述步骤1得到的X向位移传感器(2)的数据Dxij和Y向位移传感器(3)的数据Dyij以及检棒(4)半径rb,利用平面圆方程(x-x0)2+(y-y0)2=rb 2计算出当前测量平面下检棒(4)圆心坐标序列rj,j=(1,2,3,4,5,6);
步骤3:排除***误差:
将所述步骤2得到的检棒(4)圆心坐标序列rj进行最小二乘圆拟合,得到转台(5)轴心坐标序列Ri,i=(1,2,3,4,5);
步骤4:求转台(5)轴线表达式:
将所述步骤3得到的转台(5)轴心坐标序列Ri进行空间直线最小二乘拟合,求出转台(5)轴线空间表达式为:
<mrow>
<mfrac>
<mrow>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>x</mi>
<mn>0</mn>
</mrow>
</msub>
</mrow>
<mi>X</mi>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>y</mi>
<mn>0</mn>
</mrow>
</msub>
</mrow>
<mi>Y</mi>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mi>z</mi>
<mo>-</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>z</mi>
<mn>0</mn>
</mrow>
</msub>
</mrow>
<mi>Z</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式(1)中,M(Rx0,Ry0,Rz0)为转台轴线上一点,a=(X,Y,Z)为转台轴线方向向量;
步骤5:热误差建模,具体方法为:
步骤5.1:根据所述步骤3得出的转台(5)轴心坐标序列Ri,求出测量坐标系X、Y方向分量Xi和Yi,再根据所述步骤4得出的转台轴线表达式导出转台(5)轴线角度热误差序列:θp=(θ1,θ2,...θn),其中p为测量数据组数,p=(1,2,……,n);
步骤5.2:提取温度传感器数据,记为Tpq,q为温度传感器编号,q=(1,2,3,4,5,6,7,8,9,10),p为测量数据组数,p=(1,2,……,n);
步骤5.3:根据相关性判别公式:
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>R</mi>
<mi>q</mi>
</msub>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>-</mo>
<mover>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>p</mi>
</msub>
<mo>-</mo>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>p</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>-</mo>
<mover>
<msub>
<mi>T</mi>
<mi>q</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>p</mi>
</msub>
<mo>-</mo>
<mover>
<mi>R</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
1
式(2)中,Rp为某一界面内检棒圆心坐标随着温度变化的热误差序列,令:其中,λ为相关性系数阀值,
利用多元线性回归理论建立热误差模型:
<mrow>
<mi>X</mi>
<mo>=</mo>
<mi>T</mi>
<mo>*</mo>
<mi>A</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mn>1</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>2</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>Y</mi>
<mo>=</mo>
<mi>T</mi>
<mo>*</mo>
<mi>B</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mn>1</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>2</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mi>T</mi>
<mo>*</mo>
<mi>C</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mn>1</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>2</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mn>2</mn>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>m</mi>
<mi>k</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mi>k</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
上式中,A、B、C分别代表热误差序列X、Y、θ的热误差模型中的温度变量系数向量;
步骤5.4:根据最小二乘原理求解系数向量A=(TTT)-1TTX,B=(TTT)-1TTY,C=(TTT)-1TTθ,求得系数向量A、B、C;
步骤5.5:将所述步骤5.4得出的系数向量A、B、C代入热误差模型中即得到磨齿机的热误差。
2.根据权利要求1所述的测量热误差的方法,其特征在于,所述步骤1具体按照以下步骤实施:
步骤1.1:沿坐标系Z方向每隔90mm-110mm标记5个与坐标系Z轴垂直的测量平面,记为Si,其中,i=(1,2,3,4,5);
步骤1.2:将X向位移传感器(2)和Y向位移传感器(3)置于所述步骤1.1所述的测量平面Si内,令转台(5)旋转360°,期间每旋转60°记录X向位移传感器(2)的数据Dxij和Y向位移传感器(3)的数据Dyij,其中,i=(1,2,3,4,5),j=(1,2,3,4,5,6);
步骤1.3:间隔10-15分钟执行所述步骤1.2一次,完成p组测量数据采集。
3.根据权利要求1所述的测量热误差的方法,其特征在于,所述步骤4求转台(5)轴线表达式的具体方法为:
设di为转台轴心坐标序列Ri到轴线距离的平方和,
则:令:
求解后得到转台轴线上一点M(Rx0,Ry0,Rz0)以及转台轴线方向向量a=(X,Y,Z);
将M(Rx0,Ry0,Rz0)和a=(X,Y,Z)代入中得出转台轴线表达式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510236553.1A CN105058163B (zh) | 2015-05-11 | 2015-05-11 | 磨齿机工作状态时热误差测量装置及测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510236553.1A CN105058163B (zh) | 2015-05-11 | 2015-05-11 | 磨齿机工作状态时热误差测量装置及测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105058163A CN105058163A (zh) | 2015-11-18 |
CN105058163B true CN105058163B (zh) | 2017-08-25 |
Family
ID=54487788
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510236553.1A Expired - Fee Related CN105058163B (zh) | 2015-05-11 | 2015-05-11 | 磨齿机工作状态时热误差测量装置及测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105058163B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108838627B (zh) * | 2018-07-05 | 2020-04-28 | 大连理工大学 | 一种用于直观检测立式铣床热误差变化的样件加工方法 |
CN108857574B (zh) * | 2018-07-05 | 2020-01-17 | 大连理工大学 | 数控机床主轴的“热误差-温度”环的应用方法 |
CN112476057A (zh) * | 2020-11-25 | 2021-03-12 | 华中科技大学 | 一种数控机床丝杆热误差测量装置及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101972948A (zh) * | 2010-09-26 | 2011-02-16 | 天津大学 | 模拟工况载荷条件下机床主轴热误差试验装置 |
CN102176135A (zh) * | 2011-01-30 | 2011-09-07 | 合肥工业大学 | 数控机床热误差测量集成*** |
CN102759900A (zh) * | 2012-06-29 | 2012-10-31 | 上海三一精机有限公司 | 一种机床主轴热误差的测试***及测试方法 |
CN101797704B (zh) * | 2009-12-31 | 2013-10-30 | 重庆机床(集团)有限责任公司 | 数控滚齿机热变形误差补偿方法 |
CN103801987A (zh) * | 2014-02-25 | 2014-05-21 | 合肥工业大学 | 数控机床主轴旋转热误差测量数据精度的提升方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102011113486A1 (de) * | 2011-09-15 | 2013-03-21 | Liebherr-Verzahntechnik Gmbh | Innenverzahnungsschleifmaschine |
-
2015
- 2015-05-11 CN CN201510236553.1A patent/CN105058163B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101797704B (zh) * | 2009-12-31 | 2013-10-30 | 重庆机床(集团)有限责任公司 | 数控滚齿机热变形误差补偿方法 |
CN101972948A (zh) * | 2010-09-26 | 2011-02-16 | 天津大学 | 模拟工况载荷条件下机床主轴热误差试验装置 |
CN102176135A (zh) * | 2011-01-30 | 2011-09-07 | 合肥工业大学 | 数控机床热误差测量集成*** |
CN102759900A (zh) * | 2012-06-29 | 2012-10-31 | 上海三一精机有限公司 | 一种机床主轴热误差的测试***及测试方法 |
CN103801987A (zh) * | 2014-02-25 | 2014-05-21 | 合肥工业大学 | 数控机床主轴旋转热误差测量数据精度的提升方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105058163A (zh) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108227622B (zh) | 基于一维测头与标准球的机床几何误差测量方法及*** | |
CN103878641B (zh) | 一种五轴数控机床通用的旋转轴几何误差辨识方法 | |
CN111215967B (zh) | 数控机床动态精度检测装置及方法 | |
US10422619B2 (en) | Identification of geometric deviations of a motion guide in a coordinate-measuring machine or in a machine tool | |
CN106736848A (zh) | 数控车床热误差测量补偿***及补偿方法 | |
CN108176897A (zh) | 毛刺去除装置 | |
CN105269404B (zh) | 数控机床刀尖动态特性精度检测装置及其方法 | |
CN105058163B (zh) | 磨齿机工作状态时热误差测量装置及测量方法 | |
CN102768028A (zh) | 单关节臂在线原位测量方法及装置 | |
CN102689234B (zh) | 一种滚齿机床热误差的检测装置及检测方法 | |
CN103372666B (zh) | 柴油机供油单元壳体同轴孔加工方法及同轴度测量工具 | |
CN108723894B (zh) | 一种数控机床分离主轴末端热变形的测量***及其方法 | |
Chen et al. | Prediction and identification of rotary axes error of non-orthogonal five-axis machine tool | |
CN104308657A (zh) | 一种基于球杆仪测量的机床旋转轴几何误差六圈辨识方法 | |
CN102937409A (zh) | 一种极坐标齿轮测量中心及其零点标定方法 | |
CN105783845B (zh) | 一种数控磨齿机在机测量***的齿廓测量方法 | |
CN108620952A (zh) | 一种孔内径在线测量方法 | |
CN107234460A (zh) | 大型原位可重构数字加工方法及球体加工装置及方法 | |
CN109781778A (zh) | 中空冷却结构电主轴内冷工况下的热特性测量装置与方法 | |
CN106736863B (zh) | 一种测量镗铣床空间热变形误差的快速测量规及其方法 | |
Zhang et al. | Three-point method for measuring the geometric error components of linear and rotary axes based on sequential multilateration | |
CN107984300A (zh) | 一种数控机床转台位置无关几何误差的测量方法及*** | |
CN102788550A (zh) | 多传感器形位误差可视化综合测量仪 | |
CN108332642B (zh) | 一种直角头精度检测方法 | |
CN204188144U (zh) | 一种主轴跳动测试平台 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170825 |
|
CF01 | Termination of patent right due to non-payment of annual fee |