WO2019144255A1 - 一种数控机床直线运动轴几何精度快速检测方法 - Google Patents

一种数控机床直线运动轴几何精度快速检测方法 Download PDF

Info

Publication number
WO2019144255A1
WO2019144255A1 PCT/CN2018/073695 CN2018073695W WO2019144255A1 WO 2019144255 A1 WO2019144255 A1 WO 2019144255A1 CN 2018073695 W CN2018073695 W CN 2018073695W WO 2019144255 A1 WO2019144255 A1 WO 2019144255A1
Authority
WO
WIPO (PCT)
Prior art keywords
speed
measurement
linear motion
measuring
measuring point
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.)
Ceased
Application number
PCT/CN2018/073695
Other languages
English (en)
French (fr)
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.)
Dalian University of Technology
Original Assignee
Dalian 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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to PCT/CN2018/073695 priority Critical patent/WO2019144255A1/zh
Priority to US16/462,234 priority patent/US11125552B2/en
Publication of WO2019144255A1 publication Critical patent/WO2019144255A1/zh
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/26Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes
    • G01B11/27Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes for testing the alignment of axes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/22Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring angles or tapers; for testing the alignment of axes
    • G01B21/24Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring angles or tapers; for testing the alignment of axes for testing alignment of axes

Definitions

  • the invention belongs to the technical field of numerical control machine tool detection, and particularly relates to a rapid detection method for geometric precision of a linear motion axis of a numerical control machine tool.
  • Linear motion axis is the basic component of CNC machine tools, and its precision performance is the key to affect the quality of CNC machine tools.
  • Geometric accuracy is an important indicator of the linear motion axis and will inevitably deteriorate with the service of the machine tool. As the degree of geometric accuracy declines, CNC machine tools can not only meet the processing requirements of the product, but even cause the friction pair of the key parts of the machine tool to wear and tear, resulting in the scrapping of the machine tool. Therefore, the acquisition of the geometric accuracy state of CNC machine tools is of great significance for the accuracy compensation and machine tool maintenance of the machine tool in the service process. In actual production, machine tool manufacturers and user units currently use laser interferometers, ballbars and other testing instruments to regularly measure the geometric accuracy of machine tools.
  • the research shows that the geometrical accuracy of the linear motion axis of CNC machine tools is fast, and it needs to meet the basic conditions of high measurement efficiency, strong data processing capability, small size, low price and convenient debugging. It has great engineering challenges.
  • the accelerometer is used to measure the acceleration information perpendicular to the direction of motion when the linear motion axis moves.
  • the quadratic time domain integral obtains the motion axis deformation, and the geometric accuracy is estimated, which provides the possibility for the rapid detection of the linear motion axis geometric precision of the CNC machine tool.
  • the object of the invention is to overcome the deficiencies of the prior methods, and to solve the problem of rapid detection of the linear motion axis linear precision of the numerical control machine tool, and to invent a rapid detection method for the linear motion axis geometric precision of the numerical control machine tool.
  • the method adopts the device as an acceleration sensor, has the advantages of strong anti-interference ability and high resolution; the number of device components is small, the cost is low, the volume is small, the structure is compact and simple, and the installation is convenient.
  • the method utilizes the same geometric error to the noise distribution law at different speeds, preserves the high SNR component by filtering, and provides the data basis for the accurate calculation of the geometric accuracy of the linear axis.
  • the high signal-to-noise ratio in each speed is integrated.
  • the measurement data of the bandwidth increases the measurement bandwidth and lays a foundation for the accurate calculation of the geometric accuracy of the machine tool.
  • the measurement system is installed on the linear motion axis of the machine tool, which can realize the rapid measurement of the geometric precision of the linear motion axis of the machine tool, the operation is simple, and the equipment integration is high. Easy to automate.
  • a rapid detection method for linear motion axis linearity of numerical control machine tool adopts an acceleration sensor to measure the acceleration of the linear motion axis when moving in a uniform motion, and obtains the linear motion axis deformation through the second integral, and then calculates the straightness of the linear motion axis.
  • the measuring device is mounted on the linear motion axis, and the linear motion axis moves at three different speeds at a constant speed.
  • the data acquisition card automatically performs multi-channel acquisition and storage of motion point acceleration data; then, based on the same geometric error signal It can be decomposed into different frequency components to filter the acceleration data at different speeds.
  • the displacement data is obtained by quadratic integration of the filtered acceleration data in time domain, and the data of the displacement data at three speeds is superimposed to complete the linear motion. Calculation of the axis deformation; calculating the straightness of the linear motion axis by the end point connection method, and completing the rapid measurement of the straightness of the linear motion axis of the machine tool;
  • the first step is to assemble a linear measuring device for linear motion axes of CNC machine tools.
  • Two uniaxial acceleration sensors are used: 1 # uniaxial acceleration sensor 1 and 2 # uniaxial acceleration sensor 2; two-axis acceleration sensor group; 1 # uniaxial acceleration sensor 1 is measured in the n direction, 2 # uniaxial acceleration
  • the measuring direction of the sensor 2 is t direction, and the two single-axis acceleration sensors are respectively fixed in the detecting box 3 by nuts;
  • the setting direction I is the direction of the linear motion axis, and the detecting box 3 is mounted on the linear motion axis of the numerical control machine tool through the side flange 4 Upper, complete the clamping of the linear motion axis of the measuring device and the numerical control machine tool;
  • the data is continuously moved from one side of the linear motion axis to the other side, and the data acquisition system of the upper computer stores the coordinates of each measuring point and outputs digital signal data; wherein the moving speed is:
  • the maximum feed rate for the linear motion axis in the I direction; the sampling frequencies corresponding to the three speeds are:
  • the sampling frequency of the sensor for high speed motion For the sampling frequency of the sensor under medium speed motion, The sampling frequency of the sensor for low speed motion;
  • Measurement obtained the measured acceleration subset of the I-direction linear motion axis at three speeds k is the uniaxial acceleration sensor label, val is the measurement speed,
  • val is the measurement speed
  • the output signal of the i-th measuring point of the linear motion axis in the I direction r is the number of measuring points in the linear trajectory of the linear axis, The measurement time of the i-th measuring point of the linear motion axis in the I direction of the k # uniaxial acceleration sensor under the condition that the measuring speed is val;
  • the third step is based on the filtering of spatial frequency continuous conditions.
  • the linear motion axes of the obtained I direction are respectively at high speed.
  • Filtering the upper limit frequency for high speed motion Filtering the lower limit frequency for high speed motion; Filter the upper limit frequency for medium speed motion, Filtering the lower limit frequency for medium speed motion; Filtering the upper limit frequency for low speed motion, Filtering the lower limit frequency for low speed motion;
  • the filtered measurement acceleration subset is obtained.
  • the filtered signal K # a uniaxial acceleration sensor is at a measuring speed val conditions, the I signal axis direction of the i-th measuring point filtering;
  • the sampling frequency of the sensor under high-speed motion The maximum feedrate for the linear motion axis in the I direction.
  • the set is the total set of measured displacements
  • the invention has the beneficial effects that the acceleration sensor is used to measure the acceleration of the linear motion axis when moving in a uniform direction, and the deformation of the linear motion axis is obtained by the second integral, thereby realizing the rapid measurement of the geometric precision of the linear motion axis of the numerical control machine tool. .
  • Figure 1 is a block diagram of a geometric accuracy rapid detection device.
  • Figure 2 is a diagram of the geometric accuracy rapid detection device detection.
  • FIG. 4 is a schematic diagram showing the distribution of the filtered measurement acceleration subsets in the n direction at the three speeds of the I-direction motion axis.
  • Fig. 5 is a schematic diagram showing the distribution of displacement subsets in the n direction under low speed and medium speed high speed measurement conditions.
  • Fig. 6 is a schematic diagram showing the distribution of the displacement total set in the n direction.
  • I, n, t are I direction, n direction and t direction.
  • the CNC machine tool to be tested is the vertical machining center bed, and the guide rail length is 1500mm.
  • the first step is to assemble a linear motion axis linear precision detection device for CNC machine tools.
  • Two single-axis accelerometers 1 # single-axis accelerometer 1, 2 # uniaxial accelerometer 2 consisting of 2-way accelerometers; guaranteed 1 # uniaxial accelerometer 1 measuring direction is n-direction, 2 # uniaxial accelerometer 2
  • the measuring direction is the t direction, and the uniaxial acceleration sensor is fixed in the detecting box 3 by the nut respectively;
  • the direction I is the direction of the linear motion axis, and the detecting box 3 is fixed on the linear motion axis through the side flange to complete the measuring device and the numerical control machine tool The clamping of the linear motion axis.
  • the third step is based on the filtering of spatial frequency continuous conditions.
  • the measured acceleration subset under measurement conditions is filtered.
  • the measurement upper limit of the measurement acceleration subset obtained under the high-speed measurement condition is The lower filter limit is The upper limit of the measurement of the acceleration subset is measured under medium speed measurement conditions.
  • the lower filter limit is The upper limit of the measurement of the acceleration subset is measured under low speed measurement conditions.
  • the lower filter limit is Obtained as described above: Figure 4-I direction motion axis is filtered in the n direction at three speeds to measure the acceleration subset distribution.
  • Measurement of acceleration subsets by filtering using discrete-point time domain integration Perform double integration and add the elements of the measured displacement subset under the low speed and medium speed high speed measurement conditions to obtain:
  • Figure 5 Measurement of the displacement subset distribution in the n direction and the measurement of the total displacement distribution in the 6-n direction under low speed and medium speed high speed measurement conditions.
  • the straightness of the 1# sensor measurement direction obtained by the endpoint connection method is 2.65 ⁇ m.
  • the invention relates to a rapid detection method for geometric precision of a linear motion axis of a machine tool, which measures an acceleration of a linear motion axis when moving in a uniform direction by an acceleration sensor, and obtains a linear motion axis deformation by a second integral, and then calculates a linear motion axis.
  • the straightness provides the possibility for rapid detection of geometric accuracy of CNC machine tools.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Machine Tool Sensing Apparatuses (AREA)

Abstract

一种数控机床直线运动轴几何精度快速检测方法,属于数控机床检测技术领域。该方法首先将测量装置安装于直线运动轴,并随着直线运动轴以三种不同速度匀速运动,上层测量系统自动进行多通道采集、存储运动测点加速度数据;然后,基于相同几何误差信号可分解为不同频率分量,对不同测速下加速度信号进行滤波;最后,对滤波后的各加速度数据时域二次积分得到位移数据,并对三种测速下的位移数据进行数据叠加,完成直线运动轴变形计算;通过端点连线法计算直线运动轴的直线度,完成机床直线运动轴的直线度快速测量。本发明具有调试方便、测量效率高、数据处理能力强的优点,可实现机床直线运动轴几何精度的快速测量,设备集成度高、便于实现自动化。

Description

一种数控机床直线运动轴几何精度快速检测方法 技术领域
本发明属于数控机床检测技术领域,特别涉及一种数控机床直线运动轴几何精度快速检测方法。
背景技术
直线运动轴是数控机床的基本组成部分,其精度性能是影响数控机床整机加工质量的关键。几何精度是直线运动轴的重要精度指标,会随着机床的服役不可避免地恶化。随着几何精度衰退程度的增加,数控机床不仅无法满足产品的加工要求,甚至致使机床关键零件的摩擦副过渡磨损,导致机床报废。因此,数控机床几何精度状态的获取,对于机床在服役过程中精度补偿与机床维护具有重要意义。实际生产中,机床厂家与用户单位目前采用激光干涉仪、球杆仪等检测仪器对机床几何精度进行定期测量。然而,这种检测仪器的体积大、造价高、调试复杂、测量耗时长,严重制约机床精度检测效率,势必影响正常生产进程。为此,寻求一种数控机床直线运动轴几何精度快速检测方法,迫在眉睫。
研究表明,数控机床直线运动轴几何精度快速检测,需满足测量效率高、数据处理能力强、装置体积小、价格低廉、调试方便等基本条件,具有极大的工程挑战性。通过加速度传感器测量直线运动轴运动时垂直于运动方向的加速度信息,二次时域积分得到运动轴变形,完成几何精度的预估,为数控机床直线运动轴几何精度快速检测提供了可能性。
2010年西安交通大学郭俊杰、王金栋等在专利发明CN102062575A中公开了一种基于激光多路分时测量的数控机床几何精度检测方法,该方法采用一台激光跟踪仪先后在不同的基点位置对机床相同的3D空间进给运动进行测量,具 有较高的精度。2008年海克斯康测量技术有限公司与国家测量技术研究院亚历山德罗·巴尔萨莫、米凯莱·韦尔迪在专利发明CN101349556中公开了一种确定机床或测量机中的几何误差的方法,该方法通过测量空间内移动标靶的移动单元实现几何误差测量。然而,上述测量方法采用的测量装置体积大、造价高、测量步骤繁杂,无法满足几何精度快速测量的需求。
发明内容
本发明目的在于克服现有方法不足,针对数控机床直线运动轴几何精度快速检测问题,发明了一种数控机床直线运动轴几何精度快速检测方法。该方法采用装置为加速度传感器,具有抗干扰能力强,分辨率高等优势;装置组件数量少,成本低、体积小、结构紧凑简单,便于安装。该方法利用相同几何误差于不同测速下的噪声分布规律,通过滤波保留高信噪比成分,为线性轴几何精度的精确计算提供数据依据;基于空间频率连续分布原则,融合各测速中高信噪比带宽的测量数据,增加了测量带宽,为机床几何精度的精确计算奠定基础;将测量系统安装于机床直线运动轴,可实现机床直线运动轴几何精度的快速测量,操作简单,设备集成度高、便于实现自动化。
本发明的技术方案:
一种数控机床直线运动轴几何精度快速检测方法,采用加速度传感器测量直线运动轴匀速运动时垂于运动方向的加速度,通过二次积分得到直线运动轴的变形,进而计算得到直线运动轴的直线度;首先,将测量装置安装于直线运动轴上,并随着直线运动轴以三种不同速度匀速运动,数据采集卡自动进行多通道采集和存储运动测点加速度数据;然后,基于相同几何误差信号可分解为不同频率分量,对不同速度下加速度数据进行滤波;最后,对滤波后的各加速度数据时域二次积分得到位移数据,并对三种速度下的位移数据进行数据叠加, 完成直线运动轴变形计算;通过端点连线法计算直线运动轴的直线度,完成机床直线运动轴的直线度快速测量;
具体步骤如下:
第一步,组装数控机床直线运动轴几何精度快速测量装置
采用2个单轴加速度传感器:1 #单轴加速度传感器1和2 #单轴加速度传感器2组成的两向加速度传感器组;保证1 #单轴加速度传感器1测量方向为n方向,2 #单轴加速度传感器2测量方向为t方向,两个单轴加速度传感器分别通过螺母固定在检测盒3中;设定方向I为直线运动轴方向,检测盒3通过侧边法兰安装在数控机床直线运动轴4上,完成测量装置与数控机床直线运动轴的装夹;
第二步,直线运动轴几何精度快速测量
将测量装置沿方向I以高速
Figure PCTCN2018073695-appb-000001
中速
Figure PCTCN2018073695-appb-000002
低速
Figure PCTCN2018073695-appb-000003
分别从直线运动轴的一侧匀速运动到另一侧,同时上位机的数据采集系统高频存储各测点坐标及输出数字信号数据;其中运动速度为:
Figure PCTCN2018073695-appb-000004
其中,
Figure PCTCN2018073695-appb-000005
为I方向直线运动轴的最大进给速度;三种速度对应的采样频率分别为:
Figure PCTCN2018073695-appb-000006
其中,
Figure PCTCN2018073695-appb-000007
为高速运动下传感器的采样频率, 为中速运动下传感器的采样频率,
Figure PCTCN2018073695-appb-000009
为低速运动下传感器的采样频率;
测量获得I方向直线运动轴在三种速度下的测量加速度子集
Figure PCTCN2018073695-appb-000010
k为单轴加速度传感器标号,val为测量速度,
Figure PCTCN2018073695-appb-000011
为k #单轴加速度传感器在测量速度val条件下、于I方向的直线运动轴第i测点输出信号,r为直线轴运动轨迹内的测点数,
Figure PCTCN2018073695-appb-000012
为k #单轴加速度传感器在测量速度为val条件下、于I方向的直线运动轴第i测点的测量时间;
Figure PCTCN2018073695-appb-000013
其中,
Figure PCTCN2018073695-appb-000014
为高速运动下第i测点的测量时间,
Figure PCTCN2018073695-appb-000015
为中速运动下第i测点的测量时间,
Figure PCTCN2018073695-appb-000016
为低速运动下第i测点的测量时间;
第三步,基于空间频率连续条件的滤波
首先,分别对获得的I方向的直线运动轴在高速
Figure PCTCN2018073695-appb-000017
中速
Figure PCTCN2018073695-appb-000018
低速
Figure PCTCN2018073695-appb-000019
测量条件下的测量加速度子集
Figure PCTCN2018073695-appb-000020
Figure PCTCN2018073695-appb-000021
进行滤波,在高速、中速和低速的测量速度val条件下,滤波上限频率
Figure PCTCN2018073695-appb-000022
与下限频率
Figure PCTCN2018073695-appb-000023
Figure PCTCN2018073695-appb-000024
Figure PCTCN2018073695-appb-000025
Figure PCTCN2018073695-appb-000026
其中,
Figure PCTCN2018073695-appb-000027
为高速运动下滤波上限频率,
Figure PCTCN2018073695-appb-000028
为高速运动下滤波下限频率;
Figure PCTCN2018073695-appb-000029
为中速运动下滤波上限频率,
Figure PCTCN2018073695-appb-000030
为中速运动下滤波下限频率;
Figure PCTCN2018073695-appb-000031
为低速运动下滤波上限频率,
Figure PCTCN2018073695-appb-000032
为低速运动下滤波下限频率;
按上述操作,得到滤波测量加速度子集
Figure PCTCN2018073695-appb-000033
其中,滤波信号
Figure PCTCN2018073695-appb-000034
为k #单轴加速度传感器在测量速度val条件下、于I方向运动轴第i测点滤波后的信号;
第四步,数据融合与直线度计算
采用时域积分法对滤波测量加速度子集
Figure PCTCN2018073695-appb-000035
积分,得到测量速度子集
Figure PCTCN2018073695-appb-000036
其中
Figure PCTCN2018073695-appb-000037
其中,
Figure PCTCN2018073695-appb-000038
为第i测点的测量速度,
Figure PCTCN2018073695-appb-000039
为第i+1测点滤波后的信号,
Figure PCTCN2018073695-appb-000040
为第m测点滤波后的信号,
Figure PCTCN2018073695-appb-000041
为第1测点滤波后的信号,
Figure PCTCN2018073695-appb-000042
为测量速度val条件下滤波频率;
采用时域积分法对测量速度子集
Figure PCTCN2018073695-appb-000043
进行积分,得到测量位移子集
Figure PCTCN2018073695-appb-000044
其中
Figure PCTCN2018073695-appb-000045
其中,
Figure PCTCN2018073695-appb-000046
为第i测点的测量位移,
Figure PCTCN2018073695-appb-000047
第i+1测点的测量速度,
Figure PCTCN2018073695-appb-000048
为第i测点的测量速度,
Figure PCTCN2018073695-appb-000049
为第1测点的测量速度;
Figure PCTCN2018073695-appb-000050
为i测点沿I方向的测量距离,
Figure PCTCN2018073695-appb-000051
其中,
Figure PCTCN2018073695-appb-000052
高速运动下传感器的采样频率,
Figure PCTCN2018073695-appb-000053
为I方向直线运动轴的最大进给速度。
将低速、中速和高速测量条件下测量得到的测量位移子集中各元素相加,完成数据融合,得到
Figure PCTCN2018073695-appb-000054
元素
Figure PCTCN2018073695-appb-000055
的集合为测量位移总集,
Figure PCTCN2018073695-appb-000056
其中,
Figure PCTCN2018073695-appb-000057
为低速、中速和高速测量条件下的测量位移之和,
Figure PCTCN2018073695-appb-000058
为低速测量条件下的测量位移,
Figure PCTCN2018073695-appb-000059
为中速测量条件下的测量位移,
Figure PCTCN2018073695-appb-000060
为高速测量条件下的测量位移;
最后,计算沿k #单轴加速度传感器测量方向的直线度;直线度计算采用端点连线法,方法如下:
测量位移总集各元素到端点连线的距离集合为
Figure PCTCN2018073695-appb-000061
Figure PCTCN2018073695-appb-000062
其中
Figure PCTCN2018073695-appb-000063
Figure PCTCN2018073695-appb-000064
为第i测点的测量位移之和,
Figure PCTCN2018073695-appb-000065
为第1测点的测量位移之和,
Figure PCTCN2018073695-appb-000066
为第r测点的测量位移之和;
Figure PCTCN2018073695-appb-000067
中,最大值为Δb max,最小值为Δb min; 则直线度f为:
①当Δb min×Δb max≤0时,
f=Δb max-Δb min
②当Δb min≥0且Δb max≥0时,
f=Δb max
③当Δb min≤0且Δb max≤0时,
f=-Δb min
本发明的有益效果:采用了加速度传感器测量直线运动轴匀速运动时垂于运动方向的加速度,并通过二次积分得到直线运动轴的变形的方法,实现了对数控机床直线运动轴几何精度快速测量。
附图说明
图1是几何精度快速检测装置组成图。
图2是几何精度快速检测装置检测图。
图4是I方向运动轴在三种速度下于n方向的滤波测量加速度子集分布示意图。
图5是低速、中速高速测量条件下于n方向测量位移子集分布示意图。
图6是n方向测量位移总集分布示意图。
图中:1-1#加速度传感器;2-2#加速度传感器;3检测盒;
4数控机床直线运动轴;I、n、t为I方向、n方向与t方向。
具体实施方式
下面将结合附图和技术方案详细说明本发明的具体实施方式。
被测数控机床为立式加工中心床身,导轨行程长1500mm。
第一步,组装数控机床直线运动轴几何精度快速检测装置
两个单轴加速度传感器:1 #单轴加速度传感器1、2 #单轴加速度传感器2组成的2向加速度传感器组;保证1 #单轴加速度传感器1测量方向为n方向,2 #单轴加速度传感器2测量方向为t方向,单轴加速度传感器分别通过螺母固定在检测盒3中;方向I为直线运动轴方向,检测盒3通过侧边法兰固定在直线运动轴上,完成测量装置与数控机床直线运动轴的装夹。
第二步,直线运动轴几何精度快速测量
I方向运动轴的最大进给速度
Figure PCTCN2018073695-appb-000068
为531.25mm/s。将测量装置沿I方向以高速
Figure PCTCN2018073695-appb-000069
中速
Figure PCTCN2018073695-appb-000070
5低=速1
Figure PCTCN2018073695-appb-000071
分别从直线运动轴的一侧匀速运动到另一侧,同时上位机的数据采集系统高频存储各测点坐标及输出数字信号数据。高速运动下传感器的采样频率为
Figure PCTCN2018073695-appb-000072
中速运动下传感器的采样频率为
Figure PCTCN2018073695-appb-000073
低速运动下传感器的采样频率为
Figure PCTCN2018073695-appb-000074
测量获得:附图3-I方向运动轴在三种速度下于n方向的测量加速度子集分布;
第三步,基于空间频率连续条件的滤波
分别对测量获得I方向运动轴在高速
Figure PCTCN2018073695-appb-000075
中速
Figure PCTCN2018073695-appb-000076
低速
Figure PCTCN2018073695-appb-000077
测量条件下的测量加速度子集进行滤波。其中,测量获得高速测量条件下测量加速度子集的滤波上限为
Figure PCTCN2018073695-appb-000078
滤波下限为
Figure PCTCN2018073695-appb-000079
在中速测量条件下测量加速度子集的滤波上限为
Figure PCTCN2018073695-appb-000080
滤波下限为
Figure PCTCN2018073695-appb-000081
在低速测量条件下测量加速度子集的滤波上限为
Figure PCTCN2018073695-appb-000082
滤波下限为
Figure PCTCN2018073695-appb-000083
按上述操作获得:附图4-I方向运动轴在三种速度下于n方向的滤波测量加速度子集分布。
第四步,数据融合与直线度计算
采用离散点的时域积分法对滤波测量加速度子集
Figure PCTCN2018073695-appb-000084
进行二重积分,并将得到的低速、中速高速测量条件下测量位移子集中各元素相加获得:
附图5-低速、中速高速测量条件下于n方向测量位移子集分布与附图6-n方向测量位移总集分布。通过端点连线法获得1#传感器测量方向的直线度为2.65μm。
本发明所述的机床直线运动轴几何精度快速检测方法,它通过加速度传感器测量直线运动轴匀速运动时垂于运动方向的加速度,通过二次积分得到直线运动轴的变形,进而计算得到直线运动轴的直线度,为数控机床几何精度快速检测提供了可能性。
以上所述一种数控机床装配变形检测与评估方法仅本发明的较佳方法,故凡依本发明专利申请范围所述的特征及原理所做的等效变化或修饰,均包括本发明专利申请范围内。

Claims (1)

  1. 一种数控机床直线运动轴几何精度快速检测方法,采用加速度传感器测量直线运动轴匀速运动时垂于运动方向的加速度,通过二次积分得到直线运动轴的变形,进而计算得到直线运动轴的直线度;首先,将测量装置安装于直线运动轴上,并随着直线运动轴以三种不同速度匀速运动,数据采集卡自动进行多通道采集和存储运动测点加速度数据;然后,基于相同几何误差信号可分解为不同频率分量,对不同速度下加速度数据进行滤波;最后,对滤波后的各加速度数据时域二次积分得到位移数据,并对三种速度下的位移数据进行数据叠加,完成直线运动轴变形计算;通过端点连线法计算直线运动轴的直线度,完成机床直线运动轴的直线度快速测量;
    其特征在于,具体步骤如下:
    第一步,组装数控机床直线运动轴几何精度快速测量装置
    采用2个单轴加速度传感器:1 #单轴加速度传感器1和2 #单轴加速度传感器2组成的两向加速度传感器组;保证1 #单轴加速度传感器1测量方向为n方向,2 #单轴加速度传感器2测量方向为t方向,两个单轴加速度传感器分别通过螺母固定在检测盒3中;设定方向I为直线运动轴方向,检测盒3通过侧边法兰安装在数控机床直线运动轴4上,完成测量装置与数控机床直线运动轴的装夹;
    第二步,直线运动轴几何精度快速测量
    将测量装置沿方向I以高速
    Figure PCTCN2018073695-appb-100001
    中速
    Figure PCTCN2018073695-appb-100002
    低速
    Figure PCTCN2018073695-appb-100003
    分别从直线运动轴的一侧匀速运动到另一侧,同时上位机的数据采集系统高频存储各测点坐标及输出数字信号数据;其中运动速度为:
    Figure PCTCN2018073695-appb-100004
    其中,
    Figure PCTCN2018073695-appb-100005
    为I方向直线运动轴的最大进给速度;三种速度对应的采样频率分别为:
    Figure PCTCN2018073695-appb-100006
    其中,
    Figure PCTCN2018073695-appb-100007
    为高速运动下传感器的采样频率,
    Figure PCTCN2018073695-appb-100008
    为中速运动下传感器的采样频率,
    Figure PCTCN2018073695-appb-100009
    为低速运动下传感器的采样频率;
    测量获得I方向直线运动轴在三种速度下的测量加速度子集
    Figure PCTCN2018073695-appb-100010
    k为单轴加速度传感器标号,val为测量速度,
    Figure PCTCN2018073695-appb-100011
    为k #单轴加速度传感器在测量速度val条件下、于I方向的直线运动轴第i测点输出信号,r为直线轴运动轨迹内的测点数,
    Figure PCTCN2018073695-appb-100012
    为k #单轴加速度传感器在测量速度为val条件下、于I方向的直线运动轴第i测点的测量时间;
    Figure PCTCN2018073695-appb-100013
    其中,
    Figure PCTCN2018073695-appb-100014
    为高速运动下第i测点的测量时间,
    Figure PCTCN2018073695-appb-100015
    为中速运动下第i测点的测量时间,
    Figure PCTCN2018073695-appb-100016
    为低速运动下第i测点的测量时间;
    第三步,基于空间频率连续条件的滤波
    首先,分别对获得的I方向的直线运动轴在高速
    Figure PCTCN2018073695-appb-100017
    中速
    Figure PCTCN2018073695-appb-100018
    低速
    Figure PCTCN2018073695-appb-100019
    测量条件下的测量加速度子集
    Figure PCTCN2018073695-appb-100020
    Figure PCTCN2018073695-appb-100021
    进行滤波,在高速、中速和低速的测量速度val条件下,滤波上限频率
    Figure PCTCN2018073695-appb-100022
    与下限频率
    Figure PCTCN2018073695-appb-100023
    Figure PCTCN2018073695-appb-100024
    Figure PCTCN2018073695-appb-100025
    Figure PCTCN2018073695-appb-100026
    其中,
    Figure PCTCN2018073695-appb-100027
    为高速运动下滤波上限频率,
    Figure PCTCN2018073695-appb-100028
    为高速运动下滤波下限频率;
    Figure PCTCN2018073695-appb-100029
    为中速运动下滤波上限频率,
    Figure PCTCN2018073695-appb-100030
    为中速运动下滤波下限频率;
    Figure PCTCN2018073695-appb-100031
    为低速运动下滤波上限频率,
    Figure PCTCN2018073695-appb-100032
    为低速运动下滤波下限频率;
    按上述操作,得到滤波测量加速度子集
    Figure PCTCN2018073695-appb-100033
    其中,滤波信号
    Figure PCTCN2018073695-appb-100034
    为k #单轴加速度传感器在测量速度val条件下、于I方向运动轴第i测点滤波后的信号;
    第四步,数据融合与直线度计算
    采用时域积分法对滤波测量加速度子集
    Figure PCTCN2018073695-appb-100035
    积分,得到测量速度子集
    Figure PCTCN2018073695-appb-100036
    其中
    Figure PCTCN2018073695-appb-100037
    其中,
    Figure PCTCN2018073695-appb-100038
    为第i测点的测量速度,
    Figure PCTCN2018073695-appb-100039
    为第i+1测点滤波后的信号,
    Figure PCTCN2018073695-appb-100040
    为第m测点滤波后的信号,
    Figure PCTCN2018073695-appb-100041
    为第1测点滤波后的信号,
    Figure PCTCN2018073695-appb-100042
    为测量速度val条件下滤波频率;
    采用时域积分法对测量速度子集
    Figure PCTCN2018073695-appb-100043
    进行积分,得到测量位移子集
    Figure PCTCN2018073695-appb-100044
    其中
    Figure PCTCN2018073695-appb-100045
    其中,
    Figure PCTCN2018073695-appb-100046
    为第i测点的测量位移,
    Figure PCTCN2018073695-appb-100047
    第i+1测点的测量速度,
    Figure PCTCN2018073695-appb-100048
    为第i测点的测量速度,
    Figure PCTCN2018073695-appb-100049
    为第1测点的测量速度;
    Figure PCTCN2018073695-appb-100050
    为i测点沿I方向的测量距离,
    Figure PCTCN2018073695-appb-100051
    其中,
    Figure PCTCN2018073695-appb-100052
    高速运动下传感器的采样频率,
    Figure PCTCN2018073695-appb-100053
    为I方向直线运动轴的最大进给速度;
    将低速、中速和高速测量条件下测量得到的测量位移子集中各元素相加,完成数据融合,得到
    Figure PCTCN2018073695-appb-100054
    元素
    Figure PCTCN2018073695-appb-100055
    的集合为测量位移总集,
    Figure PCTCN2018073695-appb-100056
    其中,
    Figure PCTCN2018073695-appb-100057
    为低速、中速和高速测量条件下的测量位移之和,
    Figure PCTCN2018073695-appb-100058
    为低速测量条 件下的测量位移,
    Figure PCTCN2018073695-appb-100059
    为中速测量条件下的测量位移,
    Figure PCTCN2018073695-appb-100060
    为高速测量条件下的测量位移;
    最后,计算沿k #单轴加速度传感器测量方向的直线度;直线度计算采用端点连线法,方法如下:
    测量位移总集各元素到端点连线的距离集合为
    Figure PCTCN2018073695-appb-100061
    Figure PCTCN2018073695-appb-100062
    其中
    Figure PCTCN2018073695-appb-100063
    为第i测点的测量位移之和,
    Figure PCTCN2018073695-appb-100064
    为第1测点的测量位移之和,
    Figure PCTCN2018073695-appb-100065
    为第r测点的测量位移之和;
    Figure PCTCN2018073695-appb-100066
    中,最大值为Δb max,最小值为Δb min;则直线度f为:
    ①当Δb min×Δb max≤0时,
    f=Δb max-Δb min
    ②当Δb min≥0且Δb max≥0时,
    f=Δb max
    ③当Δb min≤0且Δb max≤0时,
    f=-Δb min
PCT/CN2018/073695 2018-01-23 2018-01-23 一种数控机床直线运动轴几何精度快速检测方法 Ceased WO2019144255A1 (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/CN2018/073695 WO2019144255A1 (zh) 2018-01-23 2018-01-23 一种数控机床直线运动轴几何精度快速检测方法
US16/462,234 US11125552B2 (en) 2018-01-23 2018-01-23 Method for the rapid detection of the geometric accuracy of the linear motion axis of an NC machine tool

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2018/073695 WO2019144255A1 (zh) 2018-01-23 2018-01-23 一种数控机床直线运动轴几何精度快速检测方法

Publications (1)

Publication Number Publication Date
WO2019144255A1 true WO2019144255A1 (zh) 2019-08-01

Family

ID=67395121

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2018/073695 Ceased WO2019144255A1 (zh) 2018-01-23 2018-01-23 一种数控机床直线运动轴几何精度快速检测方法

Country Status (2)

Country Link
US (1) US11125552B2 (zh)
WO (1) WO2019144255A1 (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10900779B2 (en) * 2018-01-23 2021-01-26 Dalian University Of Technology Method for the rapid detection of the linear axis angular error of an NC machine tool
CN111412866B (zh) * 2020-05-13 2024-11-15 河北省计量监督检测研究院廊坊分院 一种检定刀口形直尺工作棱边直线度的移动测量装置
WO2022192155A1 (en) * 2021-03-10 2022-09-15 Rosendin Electric, Inc. Downlight laser jig
CN116255893A (zh) * 2021-12-10 2023-06-13 财团法人精密机械研究发展中心 工具机智能数字几何精度检测系统及方法
CN115839657B (zh) * 2022-09-15 2025-09-12 中国计量科学研究院 一种基于时空同步的导轨直线度误差惯性测量方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5418611A (en) * 1992-03-25 1995-05-23 Huang; Peisen Multi-degree-of-freedom geometric error measurement system
CN202393362U (zh) * 2011-12-31 2012-08-22 大连民族学院 一种长轴直线度检测系统
CN106152976A (zh) * 2016-06-28 2016-11-23 天津工业大学 一种基于角度测量的电梯导轨垂直度检测装置
CN106863014A (zh) * 2017-02-24 2017-06-20 大连理工大学 一种五轴数控机床直线轴几何误差检测方法
CN107538273A (zh) * 2016-06-29 2018-01-05 长春设备工艺研究所 大型筒形件圆度误差及圆心跳动、直线度在线检测方法及其检测装置

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1987002468A1 (en) * 1985-10-21 1987-04-23 Sundstrand Data Control, Inc. Synchronous fm digital detector
US4896268A (en) * 1987-11-25 1990-01-23 Sundstrand Data Control, Inc. Apparatus and method for processing the output signals of a coriolis rate sensor
KR101123338B1 (ko) * 2007-04-13 2012-03-27 키네틱 인코포레이티드 물체의 운동 감지 장치
FR2917175B1 (fr) * 2007-06-08 2010-04-16 Eurocopter France Procede et systeme d'estimation de la vitesse angulaire d'un mobile
US8025607B2 (en) * 2009-09-16 2011-09-27 Northeastern University Instrumented handle and pedal systems for use in rehabilitation, exercise and training equipment
US20130201316A1 (en) * 2012-01-09 2013-08-08 May Patents Ltd. System and method for server based control
US9814426B2 (en) * 2012-06-14 2017-11-14 Medibotics Llc Mobile wearable electromagnetic brain activity monitor
US20160232811A9 (en) * 2012-06-14 2016-08-11 Robert A. Connor Eyewear System for Monitoring and Modifying Nutritional Intake
US20170164878A1 (en) * 2012-06-14 2017-06-15 Medibotics Llc Wearable Technology for Non-Invasive Glucose Monitoring
US9274136B2 (en) * 2013-01-28 2016-03-01 The Regents Of The University Of California Multi-axis chip-scale MEMS inertial measurement unit (IMU) based on frequency modulation
US9519076B2 (en) * 2014-02-20 2016-12-13 Lockheed Martin Corporation De-centralized control architecture for improved sensitivity of accelerometer-based gravity gradiometers
FR3048781B1 (fr) * 2016-03-11 2018-03-16 Commissariat A L'energie Atomique Et Aux Energies Alternatives Accelerometre lineaire amorti
US11085793B2 (en) * 2016-10-03 2021-08-10 Government Of The United States Of America, As Represented By The Secretary Of Commerce Inertial measurement unit and diagnostic system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5418611A (en) * 1992-03-25 1995-05-23 Huang; Peisen Multi-degree-of-freedom geometric error measurement system
CN202393362U (zh) * 2011-12-31 2012-08-22 大连民族学院 一种长轴直线度检测系统
CN106152976A (zh) * 2016-06-28 2016-11-23 天津工业大学 一种基于角度测量的电梯导轨垂直度检测装置
CN107538273A (zh) * 2016-06-29 2018-01-05 长春设备工艺研究所 大型筒形件圆度误差及圆心跳动、直线度在线检测方法及其检测装置
CN106863014A (zh) * 2017-02-24 2017-06-20 大连理工大学 一种五轴数控机床直线轴几何误差检测方法

Also Published As

Publication number Publication date
US20200064125A1 (en) 2020-02-27
US11125552B2 (en) 2021-09-21

Similar Documents

Publication Publication Date Title
WO2019144255A1 (zh) 一种数控机床直线运动轴几何精度快速检测方法
CN201221938Y (zh) 大型圆柱工件非接触智能离线检测仪器
CN107244335B (zh) 一种波浪磨耗检测装置及方法
CN202278449U (zh) 机床结合部对刀尖点位移量的测试装置
CN105698670B (zh) 一种机床导轨安装平面平行度的快速测量方法
CN202974629U (zh) 精密弹簧片性能测试仪
CN206161290U (zh) 滚动直线导轨副性能测试装置
CN104400560A (zh) 一种数控机床切削工况下主轴轴心轨迹在线测量方法
CN102944417A (zh) 一种机床主轴静刚度的测试平台及方法
CN201397107Y (zh) 一种用于检测减振器连杆的量具
CN102998043B (zh) 用于监测悬浮状态下直线电机动态推力的装置及方法
CN101187538A (zh) 螺纹量规的检测方法
CN104005324B (zh) 一种路面构造信息的检测系统
CN108362493B (zh) 一种数控机床直线轴转角误差快速检测方法
CN108415368B (zh) 一种数控机床直线运动轴几何精度快速检测方法
CN105423977A (zh) 一种液压机滑块检测装置及其检测方法
CN113758545A (zh) 四阀门活塞式动态流量计量标准装置和计量方法
CN203337100U (zh) 结晶器振动台的三维位移检测装置
CN216206874U (zh) 四阀门活塞式动态流量计量标准装置
CN105606012A (zh) 面向精密超精密机床的单测头光学曲面在位测量方法
US10900779B2 (en) Method for the rapid detection of the linear axis angular error of an NC machine tool
CN107957254A (zh) 测量数控机床工作台瞬时运动中心的实验装置与确定方法
CN209028918U (zh) 科氏加速度测量及加速度合成定理检验实验装置
CN115808118B (zh) 粗糙度和轮廓度两用结合协调测量装置及测量方法
CN102401627B (zh) 机床结合部对刀尖点位移量的测试装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18902907

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18902907

Country of ref document: EP

Kind code of ref document: A1