CN108490879A - 一种基于近似熵的数控机床丝杠健康状态评估方法 - Google Patents
一种基于近似熵的数控机床丝杠健康状态评估方法 Download PDFInfo
- Publication number
- CN108490879A CN108490879A CN201810374717.0A CN201810374717A CN108490879A CN 108490879 A CN108490879 A CN 108490879A CN 201810374717 A CN201810374717 A CN 201810374717A CN 108490879 A CN108490879 A CN 108490879A
- Authority
- CN
- China
- Prior art keywords
- screw
- machine tool
- lead screw
- numerical control
- leading screw
- 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
Links
- 230000036541 health Effects 0.000 title claims abstract description 37
- 238000011156 evaluation Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 37
- 230000003862 health status Effects 0.000 claims abstract description 13
- 238000012360 testing method Methods 0.000 claims abstract description 11
- 230000008569 process Effects 0.000 claims abstract description 10
- 230000001133 acceleration Effects 0.000 claims abstract description 7
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000005520 cutting process Methods 0.000 claims description 2
- 238000002474 experimental method Methods 0.000 abstract 1
- 239000013598 vector Substances 0.000 description 15
- 239000000523 sample Substances 0.000 description 4
- 238000012423 maintenance Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000013074 reference sample Substances 0.000 description 2
- 241000287828 Gallus gallus Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
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/406—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 monitoring or safety
- G05B19/4065—Monitoring tool breakage, life or condition
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于近似熵的数控机床丝杠健康状态评估方法，其包括如下步骤：(a)设置评测参数，生成专用评测G代码；(b)数控机床运行评测G代码，采集丝杠工作时的传感器信号，并进行信号预处理；(c)将丝杠运动过程划分为加速、匀速、减速三个阶段，分阶段计算丝杠工作时信号的近似熵值；(d)与不同健康状态丝杠的标准样本比较，评估丝杠当前的健康状态。本发明具有使用方便、推广成本低的优点，同时不需要搭建实验平台，无需拆装丝杠，不影响丝杠动态特性，可实现丝杠健康状态的快速评估。
Description
技术领域
本发明属于丝杠健康状态评估的技术领域，更具体地，涉及一种基于近似熵的数控机床丝杠健康状态评估方法。
背景技术
丝杠是数控机床的重要组成部分，其性能直接影响到数控机床加工质量、可靠性和稳定性水平，随着数控机床加工朝着高速、高精和重载的方向发展，数控机床丝杠往往不能达到设计寿命便出现性能衰退或损坏。对数控机床丝杠健康状态进行实时评估，可在丝杠性能衰退或损坏发生前及时采取措施，提高数控机床整机的可靠性和稳定性水平，保证加工的质量和精度，提高良品率和生产效率。
目前，丝杠性能或健康状态评估技术领域已经发展了一些方法与技术，但均存在一定的局限性。
如公告号为CN201610032369，名称为“一种滚珠丝杠健康状态的评估方法”，该专利通过建立丝杠不同健康状态下传感器信号样本点在特征空间中与丝杠健康值的映射关系，而得到滚珠丝杠性能衰退程度的量化评估。但该方法未明确提取信号特征的内容，实际应用的效果存在很大不确定性，评估结果可能存在较大的偏差。
公告号为CN201610186983，名称为“一种全程实时数据统计的丝杠健康保障方法”，该专利通过数据的全程实时采集，丝杠位置区间划分，对采集的数据分区间统计，以衡量丝杠的不均衡工作状态，并对丝杠工作状况和各区间相对健康状况作出预测判断。但该方法仅能对丝杠各区间相对健康状况进行评估，无法评估丝杠整体的健康状态情况。
因此，针对现有技术的局限性，本领域的技术人员致力于开发一种适合生产现场应用的丝杠健康状态评估方法。
发明内容
针对现有技术的以上缺陷或改进需求，本发明提供了一种基于近似熵的数控机床丝杠健康状态评估方法，其通过数控机床运行专用评测G代码，计算丝杠工作时传感器信号的近似熵值，与不同健康状态丝杠的标准样本进行比较，实现丝杠当前的健康状态的实时评估。
为实现上述目的，按照本发明的一个方面，提供了一种基于近似熵的数控机床丝杠健康状态评估方法，其特征在于，包括下列步骤：
(a)根据待测数控机床与丝杠的配置，设置有关测试参数，生成专用评测G代码；
(b)数控机床运行评测G代码，采集丝杠工作时的传感器信号，并进行信号预处理；
(c)将丝杠运动过程划分为加速阶段、匀速阶段、减速阶段三个阶段，分阶段分别计算丝杠工作时在各个工作阶段的信号近似熵值ApEn；
(d)与不同健康状态丝杠的标准样本进行比较，用距离评估方法比较基准样本与本次计算结果之间的相关性，评估丝杠当前的健康状态。
作为进一步优选的，步骤(a)所述的设置有关测试参数包括设置丝杠运动行程、丝杠进给速度。
作为进一步优选的，步骤(b)所述的数控机床运行评测G代码时为空载运行，不进行实际切削，丝杠轴向载荷可认定为固定常量。
作为进一步优选的，步骤(b)所述的传感器信号为丝杠转速、丝杠进给轴电流信号。
作为进一步优选的，步骤(b)所述的信号预处理为对传感器信号进行低通滤波，以滤除噪声，提高信噪比。
作为进一步优选的，步骤(c)所述的划分丝杠过程为三个阶段是通过丝杠转速信号实现阶段划分。
作为进一步优选的，步骤(d)所述的不同健康状态丝杠包括健康丝杠、各类非健康和损坏的丝杠。
作为进一步优选的，步骤(d)所述的距离评估方法为欧几里得距离评估法。
作为进一步优选的，步骤(d)所述的丝杠当前的健康状态，若其评估结果为非健康状态，则报警提示维修人员进行处理。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果：
1、本发明通过计算丝杠工作时传感器信号的近似熵值实现丝杠健康状态的评估，与现有技术相比较，不需要额外的搭建实验平台，不需要拆装丝杠，可以在丝杠正常工作的条件下实现，不影响丝杠的动态特性；
2、本发明根据待测数控机床与丝杠的配置，设置有关测试参数，可最大限度模拟实际加工工况，相比于现有技术和方法，更符合丝杠实际工作条件，丝杠健康状态评估结果更准确。
3、通过本发明实现对丝杠健康状态的评估，能实时反馈丝杠健康状态，可在丝杠处于亚健康水平时提示维修人员及时处理，提高丝杠的可靠性水平，降低丝杠故障发生率，进而提高设备整机的可靠性水平。
附图说明
图1是按照本发明的优选实施例所构建的丝杠健康状态评估方法的流程图；
图2是按照本发明的优选实施例所构建的评测G代码示例图；
图3是按照本发明的优选实施例所构建的丝杠运动过程阶段划分示意图；
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
图1是按照本发明的优选实施例所构建的性能评估方法的流程图，如图1所示，本实例选取对象为宝鸡机床厂VMC850L加工中心的X轴丝杠，X轴左右行程800mm，丝杠类型为滚珠丝杠，数控系统为华中数控HNC-818B。本发明实施包括以下步骤：
步骤(a)：根据待测数控机床与丝杠的配置，设置有关测试参数，生成专用评测G代码。本实例中，测试参数为设置丝杠运动行程、丝杠进给速度，根据数控机床与丝杠的配置，将丝杠运动行程设置为800mm，丝杠进给速度设置为800mm/min。生成专用评测G代码如图2所示。
步骤(b)：数控机床运行评测G代码，采集丝杠工作时的传感器信号，并进行信号预处理。
本实例中，直接将丝杠寿命预测模型植入数控系统内部，通过数控系统内部获取丝杠工作时的传感器数据，采集的传感器信号包括丝杠转速、丝杠进给轴电流。然后对采集的传感器信号进行低通滤波，以滤除噪声，提高信噪比。
步骤(c)：将丝杠运动过程划分为加速阶段、匀速阶段、减速阶段三个阶段，分阶段分别计算丝杠工作时在各个工作阶段的信号近似熵值ApEn；
本实例中，通过丝杠转速将丝杠的一次运动过程划分为加速、匀速、减速三个阶段。通过将丝杠转速和丝杠进给轴电流信号对齐，可截取出加速、匀速、减速三个阶段对应的丝杠进给轴电流信号，如图3所示。然后对三个阶段的电流信号分别计算其近似熵值ApEn，计算步骤如下：
步骤(c1)：提取待求阶段的电流信号，记为{u(1),u(2),…,u(N)}，数据个数为N，预先给定模式维数m和相似容限r。
本实例中，模式维数m取值为2，相似容限r由电流信号的标准差求得，计算公式为：
r＝k*std({u(1),u(2),…,u(N)})
其中，k为比例系数，取k＝0.4。
步骤(c2)：根据序列{u(i)}按顺序依次构造m维向量X(i)：
X(i)＝[u(i),u(i+1),…,u(i+m-1)],i＝1,2,…,N-m+1
步骤(c3)：对于任意两个向量X(i)、X(j)，向量间的距离：
d[X(i),X(j)]＝max1≤k≤m(|u(i+k-1)-u(j+k-1)|)
步骤(c4)：对每一个X(i)，向量X(i)与其它所有向量X(j)(j＝1,2,…,N-m+1，j≠i)间的相似度：
其中，
步骤(c5)：定义Φm(r)：
步骤(c6)：将维数m改为m+1，重复步骤(c2)-(c5)，得到Φm+1(r)。
步骤(c7)：此阶段电流信号的近似熵值可由下式计算求得：
ApEn(m,r,N)＝Φm(r)-Φm+1(r)
步骤(d)：与不同健康状态丝杠的标准样本进行比较，用距离评估方法比较基准样本与本次计算结果之间的相关性，评估丝杠当前的健康状态。
本实例中，不同健康状态的丝杠包括健康丝杠、各类非健康或损坏的丝杠，其标准样本数据是执行上述专用评测G代码后采集传感器数据计算得到的特征向量Y：
Y＝[ApEn加速,ApEn匀速,ApEn减速]
其中，ApEn加速、ApEn匀速、ApEn减速分别为加速、匀速、减速阶段电流信号近似熵值。
具体的，在步骤(c)求得各个阶段电流信号的近似熵值后，可得丝杠当前状态的特征向量Y′。用向量间欧几里得距离d(Y′,Y)依次评估丝杠当前状态的特征向量Y′与不同的基准样本特征向量Y之间的相关性，计算公式为：
其中，n为特征向量长度，Y′i、Yi分别为丝杠当前状态的特征向量、不同基准样本特征向量的第i个元素。
设定相似度阈值T＝0.2，当丝杠当前状态的特征向量Y′与某一基准样本特征向量Y之间的欧几里得距离d(Y′,Y)小于T时，则认为该基准样的丝杠状态即为丝杠当前健康状态。若评估结果表明丝杠当前健康状态为非健康或损坏状态，则报警提示维修人员及时处理。
综上，本发明通过设置评测参数，生成专用评测G代码，使数控机床运行专用评测G代码并采集传感器信号，然后计算丝杠工作时传感器信号的近似熵值，与不同健康状态丝杠的标准样本进行比较，实现丝杠当前的健康状态的实时评估。可达到快速准确评估丝杠当前健康状态的目的，具有不需要搭建实验平台，不需要拆装丝杠，不影响数控机床加工和丝杠的动态特性等优点。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于近似熵的数控机床丝杠健康状态评估方法，其特征在于，包括如下步骤：
(a)根据待测数控机床与丝杠的配置，设置有关测试参数，生成专用评测G代码；
(b)数控机床运行评测G代码，采集丝杠工作时的传感器信号，并进行信号预处理；
(c)将丝杠运动过程划分为加速阶段、匀速阶段、减速阶段三个阶段，分阶段分别计算丝杠工作时在各个工作阶段的信号近似熵值ApEn；
(d)与不同健康状态丝杠的标准样本进行比较，用距离评估方法比较基准样本与本次计算结果之间的相关性，评估丝杠当前的健康状态。
2.如权利要求1所述的数控机床丝杠健康状态评估方法，其特征在于，所述设置有关测试参数为设置丝杠运动行程、丝杠进给速度。
3.如权利要求1所述的数控机床丝杠健康状态评估方法，其特征在于，所述评测G代码的功能是使丝杠按照设定的测试参数运动。
4.如权利要求1所述的数控机床丝杠健康状态评估方法，其特征在于，所述数控机床运行评测G代码时为空载运行，不进行实际切削，丝杠轴向载荷可认定为固定常量。
5.如权利要求1所述的数控机床丝杠健康状态评估方法，其特征在于，所述采集的传感器信号为丝杠转速、丝杠进给轴电流信号。
6.如权利要求1所述的数控机床丝杠健康状态评估方法，其特征在于，所述丝杠运动过程三个阶段的划分是通过丝杠转速实现划分。
7.如权利要求1所述的数控机床丝杠健康状态评估方法，其特征在于，所述不同健康状态丝杠包括健康丝杠、各类非健康和损坏的丝杠。
8.如权利要求1所述的数控机床丝杠健康状态评估方法，其特征在于，所述距离评估方法为欧几里得距离评估法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810374717.0A CN108490879B (zh) | 2018-04-24 | 2018-04-24 | 一种基于近似熵的数控机床丝杠健康状态评估方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810374717.0A CN108490879B (zh) | 2018-04-24 | 2018-04-24 | 一种基于近似熵的数控机床丝杠健康状态评估方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108490879A true CN108490879A (zh) | 2018-09-04 |
| CN108490879B CN108490879B (zh) | 2021-01-15 |
Family
ID=63314147
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810374717.0A Active CN108490879B (zh) | 2018-04-24 | 2018-04-24 | 一种基于近似熵的数控机床丝杠健康状态评估方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108490879B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113238528A (zh) * | 2021-05-31 | 2021-08-10 | 华中科技大学 | 一种机床健康状态实时评价方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103760820A (zh) * | 2014-02-15 | 2014-04-30 | 华中科技大学 | 数控铣床加工过程状态信息评价装置 |
| CN104569814A (zh) * | 2014-12-24 | 2015-04-29 | 南京航空航天大学 | 一种基于近似熵的直流牵引电机健康状态实时分析方法 |
| CN104808585A (zh) * | 2015-04-13 | 2015-07-29 | 华中科技大学 | 一种机床健康状态快速检查方法 |
| CN105974886A (zh) * | 2016-06-28 | 2016-09-28 | 华中科技大学 | 一种数控机床的健康监测方法 |
| CN106482639A (zh) * | 2016-10-17 | 2017-03-08 | 南京航空航天大学 | 基于近似熵计算的低速冲击位置辨识方法 |
-
2018
- 2018-04-24 CN CN201810374717.0A patent/CN108490879B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103760820A (zh) * | 2014-02-15 | 2014-04-30 | 华中科技大学 | 数控铣床加工过程状态信息评价装置 |
| CN104569814A (zh) * | 2014-12-24 | 2015-04-29 | 南京航空航天大学 | 一种基于近似熵的直流牵引电机健康状态实时分析方法 |
| CN104808585A (zh) * | 2015-04-13 | 2015-07-29 | 华中科技大学 | 一种机床健康状态快速检查方法 |
| CN105974886A (zh) * | 2016-06-28 | 2016-09-28 | 华中科技大学 | 一种数控机床的健康监测方法 |
| CN106482639A (zh) * | 2016-10-17 | 2017-03-08 | 南京航空航天大学 | 基于近似熵计算的低速冲击位置辨识方法 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113238528A (zh) * | 2021-05-31 | 2021-08-10 | 华中科技大学 | 一种机床健康状态实时评价方法及系统 |
| CN113238528B (zh) * | 2021-05-31 | 2022-08-02 | 华中科技大学 | 一种机床健康状态实时评价方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108490879B (zh) | 2021-01-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112629863B (zh) | 变工况下动态联合分布对齐网络的轴承故障诊断方法 | |
| CN108959778B (zh) | 一种基于退化模式一致性的航空发动机剩余寿命预测方法 | |
| Wang et al. | High-dimensional process monitoring and fault isolation via variable selection | |
| CN108760327B (zh) | 一种航空发动机转子故障的诊断方法 | |
| CN108490880B (zh) | 一种数控机床切削刀具磨损状态实时监测方法 | |
| CN113834657B (zh) | 一种基于改进mset及频谱特征的轴承故障预警及诊断方法 | |
| CN108181105B (zh) | 基于逻辑回归和j散度的滚动轴承故障预诊方法及系统 | |
| WO2016155241A1 (zh) | 基于Kalman滤波器的容量预测方法、系统和计算机设备 | |
| JP2015075821A (ja) | 状態診断方法および状態診断装置 | |
| WO2024065777A1 (en) | Method, apparatus, electronic device, and storage medium for diagnosing industrial fault | |
| CN108444715A (zh) | 轴承状态诊断方法、装置、存储介质及电子设备 | |
| CN115034137A (zh) | 基于rvm和退化模型的轴承剩余寿命两阶段混合预测方法 | |
| CN116861313B (zh) | 基于振动能量趋势的卡尔曼滤波工况识别方法、系统 | |
| CN110308713A (zh) | 一种基于k近邻重构的工业过程故障变量识别方法 | |
| Ma et al. | Application of variational auto-encoder in mechanical fault early warning | |
| CN110837248A (zh) | 机床健康度测评方法、系统、综合系统及数控机床 | |
| CN108490879B (zh) | 一种基于近似熵的数控机床丝杠健康状态评估方法 | |
| CN110108474B (zh) | 一种旋转机械运行稳定性在线监测与评估方法及系统 | |
| CN109612726A (zh) | 一种用于振动信号特征提取的多重超阶分析方法 | |
| CN112067298A (zh) | 一种基于层次全局模糊熵的滚动轴承故障诊断方法 | |
| CN113761650B (zh) | 一种柴电动力机车系统序贯概率比检验故障检测方法 | |
| CN114088389A (zh) | 一种齿轮变速箱的数据处理方法和相关装置 | |
| CN112621381B (zh) | 机床进给系统健康状态智能评估方法及装置 | |
| TW201633025A (zh) | 工具機主軸故障形式的診斷方法及其系統 | |
| Thanagasundram et al. | A fault detection tool using analysis from an autoregressive model pole trajectory |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220106 Address after: 441705 Wudang Road, Shihua Town, Gucheng County, Xiangyang City, Hubei Province Patentee after: HUBEI GUCHENG COUNTY DONGHUA MACHINERY CO.,LTD. Address before: 441053 Luzhong Road, Xiangcheng District, Xiangyang, Hubei Province, No. 296 Patentee before: HUBEI University OF ARTS AND SCIENCE Patentee before: Advanced Manufacturing Engineering Research Institute of Xiangyang Huazhong University of science and technology |
|
| TR01 | Transfer of patent right |