CN109255201B - 一种基于som-mqe的滚珠丝杠副健康评估方法 - Google Patents
一种基于som-mqe的滚珠丝杠副健康评估方法 Download PDFInfo
- Publication number
- CN109255201B CN109255201B CN201811253853.0A CN201811253853A CN109255201B CN 109255201 B CN109255201 B CN 109255201B CN 201811253853 A CN201811253853 A CN 201811253853A CN 109255201 B CN109255201 B CN 109255201B
- Authority
- CN
- China
- Prior art keywords
- value
- health
- point
- degradation
- screw pair
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Computer Hardware Design (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种基于SOM‑MQE的滚珠丝杠副健康评估方法,涉及滚珠丝杠副健康评估技术。该方法首先采集丝杠副***两端的轴承支座以及丝母上的优选通道振动信号,对该信号进行预处理后进行特征提取,然后将提取到的特征值通过SOM‑MQE映射为健康值;通过自适应阈值检测算法获得健康值曲线的报警阈值,当健康值超过报警阈值时,自动报警,并将出现问题的位置定位到丝杠副***具体的部件上面。本发明所述的方法能够评估丝杠副中具体部件的健康状态,为丝杠副的精准维修提供依据。
Description
技术领域
本发明涉及滚珠丝杠副健康评估技术。
背景技术
随着装备智能化水平的不断提高,以及人工智能技术的飞速发展,设备的智能化维护逐渐成为研究热点。设备的有效维护能够保障其安全高效的运行。如何在设备运行时对于其健康状况进行监测与评估,如何对健康状况较差的部件进行定位,以及如何安排维修等都是急需解决的问题。滚珠丝杠作为重要的传动部件,在数控机床等行业应用广泛,其安全运行对于整个工业设备来说至关重要。因此,有必要对其进行健康评估。
方法1:专利“一种滚珠丝杠副健康状态的评估方法”将数控机床的转速与扭矩信号作为输入,结合拉普拉斯降维与马氏距离分析模型建立了传感器信号样本点在特征空间中与健康值之间的关系,从而获得滚珠丝杠副性能衰退程度的量化评估。但是,该方法没有提供自动预警***,当丝杠副的健康状况较差时,并没有定位到具体的部件位置,这样无法对部件提供精准的维修。
方法2:数控机床滚珠丝杠副性能退化评估技术采用量子遗传算法与灰色神经网络对滚珠丝杠副的性能进行评估。但是其用于模型训练的数据为不同运行时间下的数据,仅仅将运行时间作为性能退化的量化指标,不符合实际情况。比如,运行时间为0.5年的丝杠有可能比如运行时间为1年的丝杠退化严重,因为中间的工况有可能不同。因此,通过现场经验对丝杠副退化模型的输出进行设定也不符合实际情况。所通过这些数据训练出来的模型所得到的评估结果缺乏一定的可信度。
方法3:基于性能退化模型的滚珠丝杠副寿命预测研究通过监测电流获得丝杠副的轴向力再得知转速可以获得丝杠副的期望剩余寿命,将振动信号与期望的剩余寿命一起输入到动态模糊神经网络中进行训练输出丝杠的剩余寿命,该方法侧重于寿命预测,并且期望剩余寿命通过计算获得,预测结果缺乏一定的可信度。
发明内容
本发明的目的是为了解决现有技术存在上述问题,提供一种基于SOM-MQE的滚珠丝杠副健康评估方法。
本发明所述的基于SOM-MQE的滚珠丝杠副健康评估方法包括以下步骤:
步骤一、采集丝杠副***两端的轴承支座以及丝母上的振动信号,所述振动信号通过安装在丝杠副***两端的轴承支座以及丝母上面的三通道振动传感器进行采集;
步骤二:从每个三通道振动传感器输出的三个通道信号中选择最能代表对应部件的退化趋势的通道信号作为健康评估的原始信号;
步骤三:分别将每个部件上面优选的通道信号进行预处理、并对预处理后的信号进行特征提取;
步骤四:将步骤三最终得到的特征值通过SOM-MQE映射为健康值;
步骤五:通过自适应阈值检测算法获得健康值曲线的报警阈值,所述健康值曲线为步骤四获得的健康值随时间变化的曲线;
步骤六:当健康值超过报警阈值时,自动报警,并将出现问题的位置定位到丝杠副***具体的部件上面。
进一步地,步骤三中的预处理包括小波降噪。
进一步地,步骤三中对预处理后的信号提取的特征包括:均方根值、方差、标准差、最大值、最小值、平均幅值、峭度因子、波形系数、峰值、峰值因子、脉冲指标、方根幅值、裕度系数以及偏度。
进一步地,步骤三还包括特征选择,所述特征选择是指对提取到的特征进行降维处理,以选择出对于丝杠退化贡献大的特征。
进一步地,采用Fisher准则对提取到的特征进行降维处理。
进一步地,降维处理后得到特征为:均方根值、标准差、最大值以及方差。
进一步地,步骤四的具体方法为:
步骤四一、通过正常数据来训练自组织映射神经网络;
步骤四二、输入测试数据,并找到每个测试数据的BMU;
步骤四三、计算测试数据与BMU之间的距离,该距离即为健康值。
进一步地,步骤五中的报警包括一级报警及二级报警,一级报警阈值及二级报警阈值分别对应一级临界退化点与二级临界退化点,所述一级临界退化点与二级临界退化点采用3σ与6σ准则进行检测。
本法明通过特征提取以及特征选择的方法获得振动信号的特征值,并通过自组织映射-最小量化误差法(SOM-MQE)将这些特征值映射为健康值。通过自适应阈值检测算法获得丝杠的报警阈值。在两端轴承支座以及丝母上面分别布置振动传感器,可以定位到具体的部件。本发明所述的方法能够评估丝杠副中具体部件的健康状态,为丝杠副的精准维修提供依据。
附图说明
图1为本发明实施方式中滚珠丝杠副健康评估的装置结构示意图,图中3-A指代三通道振动传感器;
图2为本发明实施方式中滚珠丝杠副***健康评估的原理示意图;
图3为本发明实施方式中滚珠丝杠副的健康值曲线,其中两个实心正方形分别代表一级退化临界点与二级退化临界点。
具体实施方式
具体实施方式一:如图1所示,本实施方式所述的一种基于SOM-MQE的滚珠丝杠副健康评估方法包括以下步骤:
步骤一、采集丝杠副***两端的轴承支座以及丝母上的振动信号,所述振动信号通过安装在丝杠副***两端的轴承支座以及丝母上面的三通道振动传感器进行采集,如图2所示。
步骤二:从每个三通道振动传感器输出的三个通道信号中选择最能代表对应部件的退化趋势的通道信号作为健康评估的原始信号。
步骤三:将三个优选的通道信号通过小波降噪进行预处理,处理后的信号作为特征提取的输入。
特征提取是健康评估中的关键环节,提取特征的好坏将影响着健康评估的准确性,为了全面考虑信号的信息,本实施方式从时域与频域空间提取信号的14个特征,分别为均方根值、方差、标准差、最大值、最小值、平均幅值、峭度因子、波形系数、峰值、峰值因子、脉冲指标、方根幅值、裕度系数以及偏度,这14个特征为丝杠副的敏感特征。
由上述14个时频特征所构成的特征向量维数较高,并且各个特征之间存在一定的耦合关系,这样会导致所观测的数据存在一定的重叠,如果把这14个特征均作为滚珠丝杠副的敏感特征,则会造成信息处理量大,不利于后续实时预测算法的实现。因此需要对提取到的特征做降维处理,以选择出最优的特征进行后续的处理。本实施方式采用Fisher准则对提取到的丝杠副的敏感特征进行降维处理,以筛选出对于丝杠退化贡献大的特征。通过降维处理所筛选出的重要特征为:均方根值、标准差、最大值以及方差。
步骤四:将步骤三最终得到的特征值通过SOM-MQE映射为健康值。健康评估是为了获得丝杠副的性能衰退曲线,并根据衰退曲线来评估其健康状态。本实施方式采用自组织映射-最小量化误差法(SOM-MQE)对丝杠健康状况进行评估。
自组织映射属于神经网络的一种算法,它可以将高维的数据映射到低维。在训练阶段,计算每个输入样本的特征与映射层神经元之间的距离,找到在映射层中距离输入样本最近的神经元,并将此神经元定义为最佳匹配单元(BMU)。
最小量化误差法是应用SOM(自组织映射)并通过测量失效数据与正常基线之间的距离来衡量丝杠的退化水平。
应用SOM-MQE的算法流程如下所示:
第一步:通过正常数据来训练自组织映射神经网络;
丝杠在运行过程中,每隔一段时间采集一定时间的振动数据,比如每隔10min采集15s的振动数据,并将这段数据记为一个样本,将该样本数据通过步骤三处理得到该样本的特征值,分别为均方根值、标准差、最大值以及方差。在丝杠正常运行状态下获得一部分样本数据,将这些数据记为正常数据。
第二步:输入测试数据,并找到每个测试数据的BMU;
测试数据是指丝杠从正常运行到完全失效在这个过程中的所有样本数据,其中每个样本数据是由步骤三中的均方根值、标准差、最大值以及方差四个特征值组成。
第三步:计算测试数据与BMU之间的距离,该距离为丝杠健康评估的健康值。
步骤五:通过自适应阈值检测算法获得健康值曲线的报警阈值,所述健康值曲线为步骤四获得的健康值随时间变化的曲线。
步骤六:当健康值超过报警阈值时,自动报警,并将出现问题的位置定位到丝杠副***的具体部件上面。
自动报警部分采用自适应阈值的检测算法自动检测健康值曲线中的临界点,根据退化情况,将临界点分为一级退化临界点与二级退化临界点。相应的在软件部分对应一级报警与二级报警。
如图3所示,一级报警对应刚开始的微弱退化阶段;二级报警对应丝杠的退化阶段。当一级报警触发时,维修人员应***件,对零部件进行维修或者更换。当一级报警触发时,如果维修人员认为设备还可以正常运行不需要维修,那么出现这种情况可能有两种原因,第一,由于数据的不确定性,出现了误报警;第二,处于微弱退化阶段,退化不明显,与正常阶段差别较小。如果一级报警没有采取措施,这时可以采用二级报警作为辅助,这种方式的好处是能够增加预警的成功率。
由于在两端轴承与中间丝母上面分别安装了振动传感器,因此当***出现报警时可以自动定位到具体的部件上面,这样可以精准地对具体部件进行维修或者更换。
自动报警部分将采用3σ与6σ准则对一级退化临界点与二级退化临界点进行检测。
3σ准则是统计学中常用的异常值处理方法,当选择3σ准则时,根据随机变量符合正态分布的情况,测量值落在(μ-3σ,μ+3σ)内的概率为99.74%,出现在此范围之外的概率为0.26%,认为是小概率事件,偏离了原状态。然而对于丝杠这种运行中不确定性较多的情况,仅仅采用一个点落在区间之外,不能够正确反映丝杠是否发生了退化,因此采用连续m个点位于区间之外来检测丝杠退化的临界点。在这里m的取值为3。
当选择6σ准则时,根据随机变量符合正态分布的情况,测量值落在(μ-6σ,μ+6σ)的概率几乎为100%,落在区间之外的概率极小。同样,对于丝杠这种运行工况较为复杂的情况,如果只采用一个点不足以确定其性能是否发生了退化,这里采用连续n个点位于区间之外用来检测丝杠退化的临界点。在这里n取值为3。
自适应阈值一级退化临界点检测与二级退化临界点检测的算法流程如下:
定义由SOM-MQE生成的健康值序列为{HV},其中,正常状态下的健康值序列为{HVNormal},新引入的健康值为HVNew。
第一步:在设备运行中,健康值是不断生成的。记当前点前面点的均值为μ,当前点前面点的方差为σ。
第二步:记当前点与前面点均值μ之间的差值为d,判断差值d与方差σ之间的大小。如果连续3个点满足d>3σ,将此点定义为一级退化临界点;如果3个点满足d>6σ,则将此点定义为二级退化临界点。
第三步:当检测到一级退化临界点与二级退化临界点中的任意一个时,会触发相应的报警,并定位到相应的部件上面。
Claims (3)
1.一种基于SOM-MQE的滚珠丝杠副健康评估方法,其特征在于,包括以下步骤:
步骤一、采集丝杠副***两端的轴承支座以及丝母上的振动信号,所述振动信号通过安装在丝杠副***两端的轴承支座以及丝母上面的三通道振动传感器进行采集;
步骤二:从每个三通道振动传感器输出的三个通道信号中选择最能代表对应部件的退化趋势的通道信号作为健康评估的原始信号;
步骤三:分别将每个部件上面优选的通道信号进行预处理、并对预处理后的信号进行特征提取;
预处理包括小波降噪;
对预处理后的信号提取的特征包括:均方根值、方差、标准差、最大值、最小值、平均幅值、峭度因子、波形系数、峰值、峰值因子、脉冲指标、方根幅值、裕度系数以及偏度;
特征选择是指对提取到的特征进行降维处理,以选择出对于丝杠退化贡献大的特征;
步骤四:将步骤三最终得到的特征值通过SOM-MQE映射为健康值,步骤四的具体方法为:
步骤四一、通过正常数据来训练自组织映射神经网络;
步骤四二、输入测试数据,并找到每个测试数据的BMU;
步骤四三、计算测试数据与BMU之间的距离,该距离即为健康值;
步骤五:通过自适应阈值检测算法获得健康值曲线的报警阈值,所述健康值曲线为步骤四获得的健康值随时间变化的曲线;报警包括一级报警及二级报警,一级报警阈值及二级报警阈值分别对应一级临界退化点与二级临界退化点,所述一级临界退化点与二级临界退化点采用3σ与6σ准则进行检测;自适应阈值一级退化临界点检测与二级退化临界点检测的算法流程如下:
定义由SOM-MQE生成的健康值序列为{HV},其中,正常状态下的健康值序列为{HVNormal},新引入的健康值为HVNew;
第一步:在设备运行中,健康值是不断生成的;记当前点前面点的均值为μ,当前点前面点的方差为σ;
第二步:记当前点与前面点均值μ之间的差值为d,判断差值d与方差σ之间的大小;如果连续3个点满足d>3σ,将此点定义为一级退化临界点;如果3个点满足d>6σ,则将此点定义为二级退化临界点;
第三步:当检测到一级退化临界点与二级退化临界点中的任意一个时,会触发相应的报警,并定位到相应的部件上面;
步骤六:当健康值超过报警阈值时,自动报警,并将出现问题的位置定位到丝杠副***具体的部件上面。
2.根据权利要求1所述的方法,其特征在于,采用Fisher准则对提取到的特征进行降维处理。
3.根据权利要求2所述的方法,其特征在于,降维处理后得到特征为:均方根值、标准差、最大值以及方差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811253853.0A CN109255201B (zh) | 2018-10-24 | 2018-10-24 | 一种基于som-mqe的滚珠丝杠副健康评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811253853.0A CN109255201B (zh) | 2018-10-24 | 2018-10-24 | 一种基于som-mqe的滚珠丝杠副健康评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109255201A CN109255201A (zh) | 2019-01-22 |
CN109255201B true CN109255201B (zh) | 2023-07-14 |
Family
ID=65046586
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811253853.0A Active CN109255201B (zh) | 2018-10-24 | 2018-10-24 | 一种基于som-mqe的滚珠丝杠副健康评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109255201B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111618655B (zh) * | 2019-11-29 | 2021-12-28 | 东莞先知大数据有限公司 | 一种数控机床滚珠丝杆健康程度量化评价方法 |
CN112101142A (zh) * | 2020-08-27 | 2020-12-18 | 深圳市行健自动化股份有限公司 | 泥浆泵运行状态评估方法、监控终端和计算机可读存储介质 |
CN112621381B (zh) * | 2020-12-25 | 2022-07-26 | 上海交通大学 | 机床进给***健康状态智能评估方法及装置 |
CN115014759A (zh) * | 2022-04-29 | 2022-09-06 | 南通市恒瑞精密机械制造有限公司 | 一种滚珠丝杠副的实时健康状态检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101870076A (zh) * | 2010-07-02 | 2010-10-27 | 西南交通大学 | 一种基于性能退化模型的数控机床导轨副寿命预测方法 |
CN101870075A (zh) * | 2010-07-02 | 2010-10-27 | 西南交通大学 | 一种基于性能退化模型的数控机床丝杠副寿命预测方法 |
CN105718876A (zh) * | 2016-01-18 | 2016-06-29 | 上海交通大学 | 一种滚珠丝杠健康状态的评估方法 |
CN108062586A (zh) * | 2017-11-30 | 2018-05-22 | 中国船舶工业***工程研究院 | 基于衰退贡献度的船舶主机关联部件状态监测方法与*** |
-
2018
- 2018-10-24 CN CN201811253853.0A patent/CN109255201B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101870076A (zh) * | 2010-07-02 | 2010-10-27 | 西南交通大学 | 一种基于性能退化模型的数控机床导轨副寿命预测方法 |
CN101870075A (zh) * | 2010-07-02 | 2010-10-27 | 西南交通大学 | 一种基于性能退化模型的数控机床丝杠副寿命预测方法 |
CN105718876A (zh) * | 2016-01-18 | 2016-06-29 | 上海交通大学 | 一种滚珠丝杠健康状态的评估方法 |
CN108062586A (zh) * | 2017-11-30 | 2018-05-22 | 中国船舶工业***工程研究院 | 基于衰退贡献度的船舶主机关联部件状态监测方法与*** |
Non-Patent Citations (2)
Title |
---|
一种新型深度自编码网络的滚动轴承健康评估方法;佘道明等;《东南大学学报(自然科学版)》;20180920(第05期);全文 * |
基于自组织映射与受限玻尔兹曼机的滚动轴承健康评估;刘浩等;《机械传动》;20170615(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109255201A (zh) | 2019-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109255201B (zh) | 一种基于som-mqe的滚珠丝杠副健康评估方法 | |
CN109298697B (zh) | 基于动态基线模型的火电厂***各部件工作状态评估方法 | |
JP7282184B2 (ja) | 工業プロセスで使用されるコンポーネントから発生する信号の異常を検出及び測定するためのシステムと方法 | |
CN113255848B (zh) | 基于大数据学习的水轮机空化声信号辨识方法 | |
US20140365179A1 (en) | Method and Apparatus for Detecting and Identifying Faults in a Process | |
CN111767930A (zh) | 物联网时序数据异常检测方法及其相关设备 | |
CN112414694B (zh) | 基于多元状态估计技术的设备多级异常状态识别方法及装置 | |
CN111538311B (zh) | 一种基于数据挖掘的机械设备柔性多状态自适应预警方法及装置 | |
US9009001B2 (en) | Monitoring for invalid data from field instruments | |
US11544554B2 (en) | Additional learning method for deterioration diagnosis system | |
CN110736968A (zh) | 基于深度学习的雷达异常状态诊断方法 | |
CN111191855B (zh) | 一种基于管网多元水质时序数据的水质异常事件识别预警方法 | |
CN115277464A (zh) | 基于多维时间序列分析的云网络变更流量异常检测方法 | |
CN113574480A (zh) | 用于预测设备损坏的装置 | |
CN116625683A (zh) | 一种风电机组轴承故障识别方法、***、装置及电子设备 | |
CN109211564B (zh) | 一种用于滚珠丝杠副健康评估的自适应阈值检测方法 | |
EP4206838A1 (en) | Forecasting and anomaly detection method for low density polyethylene autoclave reactor | |
CN112528227A (zh) | 一种基于数理统计的传感器异常数据识别方法 | |
CN117115637A (zh) | 一种基于大数据技术的水质监测预警方法及*** | |
CN116720073A (zh) | 一种基于分类器的异常检测提取方法与*** | |
CN116306262A (zh) | 一种基于集成学习算法的镀锌钢腐蚀预测方法 | |
CN110869859A (zh) | 用于确定机床主轴状态的装置和方法 | |
Aye et al. | Fault detection of slow speed bearings using an integrated approach | |
CN116956089A (zh) | 电器设备温度异常检测模型训练方法以及检测方法 | |
WO2020183781A1 (ja) | 異常診断装置 |
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 |