CN112816216A - 一种滚动轴承性能测试台和异常测试样本的识别与修正方法 - Google Patents

一种滚动轴承性能测试台和异常测试样本的识别与修正方法 Download PDF

Info

Publication number
CN112816216A
CN112816216A CN202110008755.6A CN202110008755A CN112816216A CN 112816216 A CN112816216 A CN 112816216A CN 202110008755 A CN202110008755 A CN 202110008755A CN 112816216 A CN112816216 A CN 112816216A
Authority
CN
China
Prior art keywords
data
sample
bearing
abnormal
rolling bearing
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
CN202110008755.6A
Other languages
English (en)
Other versions
CN112816216B (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.)
China Three Gorges University CTGU
Original Assignee
China Three Gorges University CTGU
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 China Three Gorges University CTGU filed Critical China Three Gorges University CTGU
Priority to CN202110008755.6A priority Critical patent/CN112816216B/zh
Publication of CN112816216A publication Critical patent/CN112816216A/zh
Application granted granted Critical
Publication of CN112816216B publication Critical patent/CN112816216B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • G01M13/04Bearings
    • G01M13/045Acoustic or vibration analysis
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Abstract

本发明提供了一种滚动轴承性能测试台和异常测试样本的识别与修正方法,测试台主要由空间环境子平台、轴承性能测试台架子平台、轴承性能数据采集子平台组成,能实现对轴承任意时刻的性能数据进行实时测试采集并存储。异常样本的识别与修正方法,针对滚动轴承全寿命周期全部采样点的性能数据进行检验,根据各个样本点间性能数据计算其相对于其它邻近样本的变化率和偏移率,与预先设定的阈值进行比较,如果变化率或者偏移率大于预先设定的阈值,则判定该轴承对应的此个样本数据异常,需要采用后续手段对其进行修正。由此即完成了对轴承性能测试中异常监测样本的识别与修正,可用于对轴承的性能进行安全测评和预测。

Description

一种滚动轴承性能测试台和异常测试样本的识别与修正方法
技术领域
本发明涉及滚动轴承性能测试,安全测评及运维保障领域,特别涉及一种滚动轴承性能测试台和针对滚动轴承测试试验中异常测试样本的识别与修正方法。
背景技术
滚动轴承是机械设备中使用最为广泛的关键零部件,而受各种因素的影响又是整个机械系统中可靠性较差的零部件,因此其运行的性能直接关乎到整个旋转机械设备的安全运行状态,如果能够搭建一套滚动轴承性能测试台对其性能进行测试,对于其运行数据进行分析识别,检测出其中的异常数据样本并进行修正,即可提高用于后续安全测评和服役性能预测数据集的数据质量,减少异常数据对于结果产生的误差影响。其研究内容主要涉及机械动力学、材料学、电气控制、信号处理技术、数据挖掘、计算机技术等,其主要方法是通过搭建的滚动轴承性能测试试验台采集数据,分析其中存在的异常数据样本,并利用权值修正法对其进行修正。
传统在滚动轴承的异常数据样本修正过程中,有异常数据识别与异常数据修正两个必要环节。在异常数据识别环节,技术人员是利用数据挖掘方法获取数据样本点的变化情况,通过逻辑算法判定其是否为异常样本数据点;在异常数据修正环节,技术人员是利用插值拟合等方法将变化异常的数据修改到正常波动范围内的数据。
发明内容
本发明所要解决的技术问题是提供一种滚动轴承性能测试台和异常测试样本的识别与修正方法,本发明基于设计的滚动轴承性能测试台,对于采集到的样本数据进行监测,结合两个逻辑算法识别其中的畸变数据和缺失数据等异常信息,并采用数据权值修正法对于异常数据进行修正,形成了由“性能测试实验→异常样本识别→异常数据修正”构建的一体化滚动轴承异常样本的识别与修正方法,同时又维持了样本数据的一致性等要求。这样的测试台不仅可以完成滚动轴承的性能测试实验,而且还可以对采集到的滚动轴承性能数据样本进行异常数据识别与修正。在搭建的滚动轴承测试台中,采用温湿度控制器模拟轴承在运行过程中各类环境,采用PLC控制器结合温湿度等输入信号控制电机、磁粉制动器等输出机构,使用振动加速度传感器采集振动信号并借助双积分信号调理器和NI9234采集卡进行数据的调理和传输;在异常监测样本的识别环节,通过样本的变化率和偏移率等指标信息识别并检验出缺失异常数据和畸变异常数据。采用变化率的波动指标识别出畸变的异常数据,同时联合同一试验台架上多个轴承通过偏移率指标判定出是否存在缺失的数据,弥补了传统滚动轴承在进行性能评估之前并未对数据质量进行识别检测的缺陷;在异常数据样本的修正环节,采用数据权值修正法修正识别为异常的数据样本。不仅能够针对单个轴承的异常数据进行修正,而且还可以对同一试验台架上多个轴承样本数据进行联合修正。这样既能够提升整体数据样本的质量,又避免了异常数据对于后续处理所带来的误差。
为了实现上述的技术特征,本发明的目的是这样实现的:一种滚动轴承性能测试台,它包括空间环境罩、温湿度控制器、轴承安装台架X轴、轴承安装台架Y轴、X轴磁粉加载机构、Y轴磁粉加载机构、X轴驱动电机、Y轴驱动电机、振动加速度传感器、双积分信号调理器、数据采集卡NI9234、计算机和PLC控制器;
所述X轴驱动电机信号控制端连接PLC控制器X轴驱动电机信号发送端,Y轴驱动电机信号控制端连接PLC控制器Y轴驱动电机信号发送端;X轴磁粉加载机构信号控制端连接PLC控制器X轴磁粉加载机构信号发送端,Y轴磁粉加载机构信号控制端连接PLC控制器Y轴磁粉加载机构信号发送端;
所述空间环境罩用于隔离外界环境干扰,温湿度控制器模拟轴承运行时的各类环境,温湿度控制器的各类环境信号控制端连接PLC控制器各类环境信号发送端;双积分信号调理器信号输入端连接振动加速度传感器的信号输出端,双积分信号调理器信号输出端连接数据采集卡NI9234的信号输入端,数据采集卡NI9234的信号输出端连接计算机的信号输入端,计算机通过信号采集软件将采集的振动数字信号显示并记录。
所述各类环境包括模拟轴承运行时的温度和湿度环境参数。
所述振动加速度传感器安装时,需要对轴承外圈表面进行清洁:先用砂纸打磨平面,然后用丙酮清洗,吹干;将适量混合好的高温胶涂抹在粘贴处,粘贴振动加速度传感器在轴承外圈上45°,135°,225°,315°位置;室温固化数小时,保持振动加速度传感器相对轴承端面平行以及相对径向垂直,不能因粘结剂的流动性导致传感器旋转;对所用振动加速度传感器进行标定,调节偏置电压使输出为零;振动加速度传感器采集到轴承振动信号后,输入到双通道电荷电压滤波积分放大器中,接端子板和数据采集卡,输入到计算机,即可完成对轴承的振动信号采集。
对任意一项所述滚动轴承性能测试台所采集的振动信号进行异常样本识别与修正的方法,针对滚动轴承全寿命周期全部采样点的性能数据进行检验,根据各个样本点间性能数据计算其相对于其它邻近样本的变化率和偏移率,与预先设定的阈值进行比较,如果变化率或者偏移率大于预先设定的阈值,则判定该轴承对应的此个样本数据异常,需要采用后续手段对其进行修正,修正之后的数据能够直接用于实现对滚动轴承服役性能进行评估和预测。
所述滚动轴承异常样本的识别与修正方法,具体包括如下步骤:
步骤1,异常监测样本的识别:
步骤1-1,针对滚动轴承全寿命周期全部采样点的性能数据进行检验,根据各个样本点间性能数据的变化率来对异常数据进行识别,依据下式求出第i个样本点的数据变化率;
Δli=(li-li-1)/li
式中:li和li-1分别是滚动轴承在第i和i-1个样本点的性能数据的统计量,Δli是滚动轴承第i个样本点相对于第i-1个样本点的数据变化率;
步骤1-2,根据下式算出前n个样本的平均数据变化;
Figure BDA0002884165390000031
式中:Δli,av为样本点i的前n个样本的平均数据变化;j代表实验台架上的轴承编号,n代表实验台架上的轴承总数目,Δlij代表实验台架上第j个轴承的第i个样本的数据变化率;
步骤1-3,如果|Δli|≥kΔli,av,则视为异常监测数据;k为样本数据的突变系数,根据实际的载荷工况特性自设定,均匀载荷时的突变系数k小于变载荷工况;
步骤1-4,选取同一个实验台架上n个轴承的全寿命周期性能测试数据作为样本,将n个轴承的全寿命周期性能测试数据看为横向量为m,纵向量为n的数组,m为轴承全寿命周期性能测试数据的样本数,根据下面式子求出各个样本点的期望均值;
Figure BDA0002884165390000041
式中:E(j)代表实验台架上n个轴承的第j个样本的期望均值,x(i,j)代表第i个轴承的第j个样本的统计量;
步骤1-5,根据下面式子求出各个样本点的方差:
Figure BDA0002884165390000042
式中:σ(i,j)为第i个轴承的第j个样本的方差;
步骤1-6,根据下面式子计算各个样本点的偏移率:
Figure BDA0002884165390000043
式中:ρ(i,j)为第i个轴承的第j个样本的偏移率;
步骤1-7,针对每个轴承全寿命周期上各个性能测试样本分别计算出他们的偏移率,与预先设定的阈值λ进行比较,如果偏移率大于阈值λ,则判定该轴承对应的此个样本数据异常;
轴承试验台架上的所有轴承在全寿命周期上测试所得的所有服役性能样本,均需要经过上步骤1-3~步骤1-7来检验,任意一个逻辑算法检定为异常样本数据,即可判定该样本为异常样本数据,需要采用后续手段对其进行修正或校正;
步骤二,异常监测样本的修正:
步骤2-1,滚动轴承的性能测试数据在时间序列范围内有n个原始样本数据,记为φn,这n个原始样本数据中,#1轴承的振动数据为a1,a2,…,an;#2轴承的振动数据为b1,b2,…,bn;#n轴承的振动数据为c1,c2,…,cn,从φn中取最近的
Figure BDA0002884165390000051
个序列数据,即:
Figure BDA0002884165390000052
步骤2-2,单个轴承的异常样本数据修正,充分利用异常样本数据点周围数据,采用数据权值修正法对插值数据进行修正,设x(t1),x(t2),…,x(tn)为单个滚动轴承在不同时刻点ti的性能测试值,取异常样本点x(ti)两端ti-1和ti+1时刻的n个测试样本数据作为该异常样本点数据x(ti)的基础修正数据,则异常样本数据x(ti)在基础测试数据x(ti+1)和x(ti-1)下的修正公式如下:
Figure BDA0002884165390000053
式中:x′(ti)为异常样本x(ti)在第1次修正之后的结果,x(ti-1,j)和x(ti+1,j)分别为第j个轴承在第i-1和第i+1个时刻点样本统计量数值;
步骤2-3,同一试验台架上多个轴承样本数据的联合修正,设同一试验台架上第j个轴承在ti时刻的样本数据分别为x(ti,j),x(ti)为待修正的异常样本数据;该试验台架上同时测得了n个轴承的实验样本数据,对φn中的样本数据,基于同一试验台架上多个轴承的其他样本数据进行联合修正,得到φ′n,并使得φ′n中的历史数据与当前数据相一致:
Figure BDA0002884165390000054
式中:x″(ti)为异常样本x(ti)在第2次修正之后的结果,x(ti,j)为第j个轴承在第i个时刻点样本统计量数值;
步骤2-4,综合步骤2-2和步骤2-3得到滚动轴承异常数据样本的最终修正结果为:
Figure BDA0002884165390000055
采用相同的步骤可以实施对实验台架上其他轴承在ti时刻的异常样本数据b(ti),c(ti)进行修正,由此即可完成对轴承性能测试中所有异常监测样本的修正。
至此,即完成了滚动轴承性能测试台的设计,性能测试和异常样本的识别与修正的全过程。本发明设计的滚动轴承性能测试台和异常样本的识别与修正方法,采用上述技术方案可以取得如下技术效果:
1、本测试台适用于滚动轴承工作过程中运行情况的性能测试,能够较好的反映出滚动轴承的运行状况,温湿度控制器用于模拟轴承在运行过程中各类环境,采用PLC控制器结合温湿度等输入信号控制电机、磁粉制动器等输出机构,使用振动加速度传感器采集振动信号并借助双积分信号调理器和NI9234采集卡进行数据的调理和传输。
2、异常监测样本的识别能够检测出滚动轴承振动序列数据中存在的异常数据,通过样本的变化率和偏移率等指标信息识别并检验出缺失数据和畸变数据。在滚动轴承的服役性能进行分析和预测之前采用异常检测样本的识别技术,弥补了传统滚动轴承在进行性能评估之前并未对数据质量进行识别检测的缺陷。
3、数据权值修正法是一种对于异常样本数据进行修正的方法,不仅能够针对单个轴承的异常数据进行修正,而且还可以对同一试验台架上多个轴承样本数据进行联合修正,较之于传统的数据插补方法,数据权值修正之后的数据与原始样本数据拟合性更高、误差更小,可以在后续中直接用于对滚动轴承服役性能进行评估和预测。
综上所述,本发明基于设计的滚动轴承性能测试台,对于采集到的样本数据进行监测,结合两个逻辑算法识别其中的畸变数据和缺失数据等异常信息,并采用数据权值修正法对于异常数据进行修正,形成了由“性能测试实验→异常样本识别→异常数据修正”构建的一体化滚动轴承异常样本的识别与修正方法,同时又维持了样本数据的一致性等要求,为滚动轴承样本质量检测与修复提供了一种新的思路方法。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或者通过本发明的实践了解到。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是本发明滚动轴承性能测试台的安装布置及区域划分示意图。
图2是本发明中滚动轴承性能测试台进行测试实验的步骤流程图。
图3是本发明中振动加速度传感器在滚动轴承上的安装角度和安装位置示意图。
图4是本发明滚动轴承性能测试异常样本的识别与修正方法的工作原理示意图。
图5是本发明步骤2-2中基础修正样本与异常样本的权值修正示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
如图1所示,一种滚动轴承性能测试台,包括空间环境罩、温湿度控制器、轴承安装台架X轴、轴承安装台架Y轴、X轴磁粉加载机构、Y轴磁粉加载机构、X轴驱动电机、Y轴驱动电机、PLC控制器、振动加速度传感器、双积分信号调理器、数据采集卡NI9234和计算机组成。
进一步的,所述空间环境罩能够将测试试验台工作环境与外部环境隔离开,避免实验台运行过程中外部环境中的粉尘、水雾、硬沙粒等进入机器中,减少硬质颗粒状异物对测试中的滚动轴承造成的损伤,同时减少外部噪声和微弱振动等所带来的实验误差,其中双积分信号调理器、数据采集卡NI9234以及计算机均位于外部环境中,而其它设备均位于空间环境罩内。
进一步的,温湿度控制器用于在封闭空间中调整温湿度等环境参数信息,模拟滚动轴承在不同的实际工况中的运行环境;轴承安装台架X轴和轴承安装台架Y轴分别用于固定其上分布的滚动轴承,为实验测试提供一个稳定的底座与台架,对于每一个安装台架其上均分布安装有两个滚动轴承,台架上连接轴承的测试轴一端通过联轴器与电机转轴相连接,同时另一端与磁粉加载机构即磁粉制动器相连接;X轴和Y轴磁粉加载机构均与对应轴末端相连接并接通到PLC控制器的输出端口,分别用于对X轴、Y轴台架上测试的滚动轴承进行加载,磁粉制动器在空载时没有转矩的输出,而处于加载状态时负载大小可以通过PLC控制器输出端控制激磁电流大小进行调节;X轴和Y轴驱动电机用于提供动力,其信号控制端接到PLC控制器的信号输出端,通过PLC输出端分别调整每个电机的转速;PLC控制器选用西门子S7-1500系列的1516-3PN/DP,其配备了电源模块、输入模块和输出模块3个扩展模块,其主要作用是接受温湿度等输入信息调整驱动电机和磁粉制动器等输出机构的运行状态;振动加速度传感器通过粘接剂粘贴在滚动轴承的外圈上,同时保持传感器相对轴承端面平行以及相对径向垂直,用于收集滚动轴承运行过程中的振动信号。
由于传感器所输出的信号较为微弱,采用双积分信号调理器用于对信号源和采集设备之间的信号进行非线性补偿等操作提高分辨率和降低噪声,信号调理包括模块主要包括电荷转换,电压放大,滤波三个部分,其输入端连接振动加速度传感器,输出端连接到数据采集卡NI9234;数据采集卡能够将调理之后的信号传输给计算机进行存储分析等,而后续的性能测试数据分析与修复等步骤均是在计算机中进行完成。
如附图1所示,该测试台架既能完成对轴承的全寿命周期的性能测试,也能实现对单一轴承任意时刻的性能测试。当对轴承的全寿命周期性能进行测试时,将待测试轴承安装至轴承台架上,调整好环境参数和载荷工况参数,粘贴好振动加速度传感器,设置采集时间间隔,假定为10分钟,则每间隔10分钟,数据采集设备自动采集1次测试数据,直至轴承磨损失效,由此即构建了该测试轴承的性能测试序列样本数据。当对单一轴承进行性能测试时,将待测试轴承安装至轴承台架上,调整好环境参数和载荷工况参数,粘贴好振动加速度传感器,即可直接实现对单一轴承的性能测试。
下面结合附图2对本发明滚动轴承性能测试的实施例作详细说明。首先,将测试的滚动轴承外圈表面打磨平整并用丙酮清洗吹干,将适量混合好的高温胶涂抹在接触处并将加速度传感器粘贴轴承外圈位置上,粘贴方式如图3所示,同时对所用传感器进行标定,调节偏置电压使输出为零。其次,将测试试验台的线路按上述说明连接好之后关上测试试验台的空间环境罩,由温湿度控制器模拟实际的运行工况,并通过线路将温湿度等信息传输到PLC控制器的输入端口,PLC控制器根据输入信息控制X轴、Y轴驱动电机输出预设的转速,磁粉制动器调节激磁电流大小控制实验的负载参数。最后,振动加速度传感器将采集的振动信号输出到处于外部环境中的双积分调理器,经过信号的处理之后通过数据采集卡NI9234将振动信号结果输出到计算机上。到此,完成了整个滚动轴承性能测试实验中数据的采集与传输部分。
下面结合附图4对本发明轴承性能测试异常样本的识别与修正方法的实施例作详细说明。本实施例的主要目的是样本的变化率和偏移率等指标信息识别并检验出轴承性能测试中的缺失数据和畸变数据,并采用数据权值修正法对单个轴承的异常数据或者同一试验台架上多个轴承样本数据进行联合修正,最终达到轴承在性能测试过程中,所有采集的性能样本数据拟合性更高、误差更小,以便在后续中直接用于对滚动轴承服役性能进行评估和预测。实施例包括如下具体步骤:
步骤1,异常监测样本的识别:
步骤1-1,针对滚动轴承全寿命周期全部采样点的性能数据进行检验,根据各个样本点间性能数据的变化率来对异常数据进行识别,依据下式求出第i个样本点的数据变化率;
检验数据是由采集来的滚动轴承采样点性能数据所组成;其中,异常数据包括畸变异常数据和缺失异常数据;数据变化率能够反映出当前样本与前后样本之间的变化幅度;
Δli=(li-li-1)/li
式中:li和li-1分别是滚动轴承在第i和i-1个样本点的性能数据的统计量,Δli是滚动轴承第i个样本点相对于第i-1个样本点的数据变化率。
步骤1-2,根据下式算出前n个样本的平均数据变化;
前n个样本是指位于测试样本点时间节点之前的n个样本数据,平均数据变化能够反映出正常数据样本的平均数据变化率,可以通过当前测试样本点的数据变化率与平均率之间的关系判定其是否属于畸变异常数据。
Figure BDA0002884165390000091
式中:Δli,av为样本点i的前n个样本的平均数据变化;j代表实验台架上的轴承编号,n代表实验台架上的轴承总数目,Δlij代表实验台架上第j个轴承的第i个样本的数据变化率。
步骤1-3,如果|Δli|≥kΔli,av,则视为异常监测数据;k为样本数据的突变系数,根据实际的载荷工况特性自设定,均匀载荷时的突变系数k小于变载荷工况;
采用当前逻辑算法与判定条件所得到的结果包含正常数据与畸变异常数据,畸变异常数据是在滚动轴承性能测试系统正常时,由于周围环境的随机扰动因素导致滚动轴承的振动序列数据骤然增大或者减小所产生的。突变系数应当根据设定的载荷工况特性进行设定,其值随着载荷工况的变化而变化。
步骤1-4,选取同一个实验台架上n个轴承的全寿命周期性能测试数据作为样本,将n个轴承的全寿命周期性能测试数据看为横向量为m,纵向量为n的数组,m为轴承全寿命周期性能测试数据的样本数,根据下面式子求出各个样本点的期望均值;
这一步骤的目的是通过同一个实验台架上多个滚动轴承判定当前样本数据是否为缺失异常数据,缺失异常数据是测试环节故障导致振动序列数据丢失或者不完整所产生的;期望均值能够反映出数据样本的平均值信息。
Figure BDA0002884165390000101
式中:E(j)代表实验台架上n个轴承的第j个样本的期望均值,x(i,j)代表第i个轴承的第j个样本的统计量。
步骤1-5,根据下面式子求出各个样本点的方差:
方差能够反映出每个样本值与全体样本平均值的变化程度,用于衡量该样本数据的波动大小,方差越大表明当前数据的波动越大,其与正常状态下的样本数据的偏差也就越大。
Figure BDA0002884165390000102
式中:σ(i,j)为第i个轴承的第j个样本的方差。
步骤1-6,根据下面式子计算各个样本点的偏移率:
偏移率是指当前数据样本与正常样本之间的偏移程度,可以通过当前测试样本点的偏移率与预设的阈值进行对比判定其是否属于缺失异常数据。
Figure BDA0002884165390000111
式中:ρ(i,j)为第i个轴承的第j个样本的偏移率。
步骤1-7,针对每个轴承全寿命周期上各个性能测试样本分别计算出他们的偏移率,与预先设定的阈值λ进行比较,如果偏移率大于阈值λ,则判定该轴承对应的此个样本数据异常;采用当前逻辑算法与判定条件所得到的结果包含正常数据与缺失异常数据;偏移率的阈值也应当根据设定的载荷工况特性进行设定,其值随着载荷工况的变化而变化。
轴承试验台架上的所有轴承在全寿命周期上测试所得的所有服役性能样本,均需要经过上步骤1-3~步骤1-7来检验,任意一个逻辑算法检定为异常样本数据,即可判定该样本为异常样本数据,需要采用后续手段对其进行修正或校正;
步骤二,异常监测样本的修正:
步骤2-1,滚动轴承的性能测试数据在时间序列范围内有n个原始样本数据,记为φn,这n个原始样本数据中,#1轴承的振动数据为a1,a2,…,an;#2轴承的振动数据为b1,b2,…,bn;#n轴承的振动数据为c1,c2,…,cn,从φn中取最近的
Figure BDA0002884165390000112
个序列数据,即:
Figure BDA0002884165390000113
对于识别为异常的数据样本点,不仅能够直接使用单个轴承的样本数据进行修正,也可以联合同一试验台上其它滚动轴承的样本数据进行联合修正。
步骤2-2,单个轴承的异常样本数据修正,充分利用异常样本数据点周围数据,采用数据权值修正法对插值数据进行修正,设x(t1),x(t2),…,x(tn)为单个滚动轴承在不同时刻点ti的性能测试值,取异常样本点x(ti)两端ti-1和ti+1时刻的n个测试样本数据作为该异常样本点数据x(ti)的基础修正数据,如附图5所示。则异常样本数据x(ti)在基础测试数据x(ti+1)和x(ti-1)下的修正公式如下:
Figure BDA0002884165390000121
式中:x′(ti)为异常样本x(ti)在第1次修正之后的结果。x(ti-1,j)和x(ti+1,j)分别为第j个轴承在第i-1和第i+1个时刻点样本统计量数值。
步骤2-3,同一试验台架上多个轴承样本数据的联合修正,设同一试验台架上第j个轴承在ti时刻的样本数据分别为x(ti,j),x(ti)为待修正的异常样本数据。该试验台架上同时测得了n个轴承的实验样本数据,对φn中的样本数据,基于同一试验台架上多个轴承的其他样本数据进行联合修正,得到φ′n,并使得φ′n中的历史数据与当前数据相一致:
Figure BDA0002884165390000122
式中:x″(ti)为异常样本x(ti)在第2次修正之后的结果。x(ti,j)为第j个轴承在第i个时刻点样本统计量数值。
步骤2-4,综合步骤2-2和步骤2-3得到滚动轴承异常数据样本的最终修正结果为:
Figure BDA0002884165390000123
采用相同的步骤可以实施对实验台架上其他轴承在ti时刻的异常样本数据b(ti),c(ti)进行修正。
针对于同一试验台架上多个轴承数据样本点,充分利用该异常数据点同一个台架上同一时间节点上其它轴承的数据,采用联合修正法对异常数据进行修正。
由此即完成了对轴承性能测试中异常监测样本的修正。
当检验数据中的每一个样本点均完成识别检验与修正之后,即完成了对于对滚动轴承性能测试中异常监测样本的识别与修正。至此,即完成了对轴承性能测试中异常监测样本的修正。

Claims (5)

1.一种滚动轴承性能测试台,其特征在于:它包括空间环境罩、温湿度控制器、轴承安装台架X轴、轴承安装台架Y轴、X轴磁粉加载机构、Y轴磁粉加载机构、X轴驱动电机、Y轴驱动电机、振动加速度传感器、双积分信号调理器、数据采集卡NI9234、计算机和PLC控制器;
所述X轴驱动电机信号控制端连接PLC控制器X轴驱动电机信号发送端,Y轴驱动电机信号控制端连接PLC控制器Y轴驱动电机信号发送端;X轴磁粉加载机构信号控制端连接PLC控制器X轴磁粉加载机构信号发送端,Y轴磁粉加载机构信号控制端连接PLC控制器Y轴磁粉加载机构信号发送端;
所述空间环境罩用于隔离外界环境干扰,温湿度控制器模拟轴承运行时的各类环境,温湿度控制器的各类环境信号控制端连接PLC控制器各类环境信号发送端;双积分信号调理器信号输入端连接振动加速度传感器的信号输出端,双积分信号调理器信号输出端连接数据采集卡NI9234的信号输入端,数据采集卡NI9234的信号输出端连接计算机的信号输入端,计算机通过信号采集软件将采集的振动数字信号显示并记录。
2.根据权利要求1所述一种滚动轴承性能测试台,其特征在于:所述各类环境包括模拟轴承运行时的温度和湿度环境参数。
3.根据权利要求1所述一种滚动轴承性能测试台,其特征在于:所述振动加速度传感器安装时,需要对轴承外圈表面进行清洁:先用砂纸打磨平面,然后用丙酮清洗,吹干;将适量混合好的高温胶涂抹在粘贴处,粘贴振动加速度传感器在轴承外圈上45°,135°,225°,315°位置;室温固化数小时,保持振动加速度传感器相对轴承端面平行以及相对径向垂直,不能因粘结剂的流动性导致传感器旋转;对所用振动加速度传感器进行标定,调节偏置电压使输出为零;振动加速度传感器采集到轴承振动信号后,输入到双通道电荷电压滤波积分放大器中,接端子板和数据采集卡,输入到计算机,即可完成对轴承的振动信号采集。
4.对权利要求1-3任意一项所述滚动轴承性能测试台所采集的振动信号进行异常样本识别与修正的方法,其特征在于,针对滚动轴承全寿命周期全部采样点的性能数据进行检验,根据各个样本点间性能数据计算其相对于其它邻近样本的变化率和偏移率,与预先设定的阈值进行比较,如果变化率或者偏移率大于预先设定的阈值,则判定该轴承对应的此个样本数据异常,需要采用后续手段对其进行修正,修正之后的数据能够直接用于实现对滚动轴承服役性能进行评估和预测。
5.根据权利要求4所述的滚动轴承异常样本的识别与修正方法,其特征在于,具体包括如下步骤:
步骤1,异常监测样本的识别:
步骤1-1,针对滚动轴承全寿命周期全部采样点的性能数据进行检验,根据各个样本点间性能数据的变化率来对异常数据进行识别,依据下式求出第i个样本点的数据变化率;
Δli=(li-li-1)/li
式中:li和li-1分别是滚动轴承在第i和i-1个样本点的性能数据的统计量,Δli是滚动轴承第i个样本点相对于第i-1个样本点的数据变化率;
步骤1-2,根据下式算出前n个样本的平均数据变化;
Figure FDA0002884165380000021
式中:Δli,av为样本点i的前n个样本的平均数据变化;j代表实验台架上的轴承编号,n代表实验台架上的轴承总数目,Δlij代表实验台架上第j个轴承的第i个样本的数据变化率;
步骤1-3,如果|Δli|≥kΔli,av,则视为异常监测数据;k为样本数据的突变系数,根据实际的载荷工况特性自设定,均匀载荷时的突变系数k小于变载荷工况;
步骤1-4,选取同一个实验台架上n个轴承的全寿命周期性能测试数据作为样本,将n个轴承的全寿命周期性能测试数据看为横向量为m,纵向量为n的数组,m为轴承全寿命周期性能测试数据的样本数,根据下面式子求出各个样本点的期望均值;
Figure FDA0002884165380000022
式中:E(j)代表实验台架上n个轴承的第j个样本的期望均值,x(i,j)代表第i个轴承的第j个样本的统计量;
步骤1-5,根据下面式子求出各个样本点的方差:
Figure FDA0002884165380000031
式中:σ(i,j)为第i个轴承的第j个样本的方差;
步骤1-6,根据下面式子计算各个样本点的偏移率:
Figure FDA0002884165380000034
式中:ρ(i,j)为第i个轴承的第j个样本的偏移率;
步骤1-7,针对每个轴承全寿命周期上各个性能测试样本分别计算出他们的偏移率,与预先设定的阈值λ进行比较,如果偏移率大于阈值λ,则判定该轴承对应的此个样本数据异常;
轴承试验台架上的所有轴承在全寿命周期上测试所得的所有服役性能样本,均需要经过上步骤1-3~步骤1-7来检验,任意一个逻辑算法检定为异常样本数据,即可判定该样本为异常样本数据,需要采用后续手段对其进行修正或校正;
步骤二,异常监测样本的修正:
步骤2-1,滚动轴承的性能测试数据在时间序列范围内有n个原始样本数据,记为φn,这n个原始样本数据中,#1轴承的振动数据为a1,a2,…,an;#2轴承的振动数据为b1,b2,…,bn;#n轴承的振动数据为c1,c2,…,cn,从φn中取最近的
Figure FDA0002884165380000032
个序列数据,即:
Figure FDA0002884165380000033
步骤2-2,单个轴承的异常样本数据修正,充分利用异常样本数据点周围数据,采用数据权值修正法对插值数据进行修正,设x(t1),x(t2),…,x(tn)为单个滚动轴承在不同时刻点ti的性能测试值,取异常样本点x(ti)两端ti-1和ti+1时刻的n个测试样本数据作为该异常样本点数据x(ti)的基础修正数据,则异常样本数据x(ti)在基础测试数据x(ti+1)和x(ti-1)下的修正公式如下:
Figure FDA0002884165380000041
式中:x′(ti)为异常样本x(ti)在第1次修正之后的结果,x(ti-1,j)和x(ti+1,j)分别为第j个轴承在第i-1和第i+1个时刻点样本统计量数值;
步骤2-3,同一试验台架上多个轴承样本数据的联合修正,设同一试验台架上第j个轴承在ti时刻的样本数据分别为x(ti,j),x(ti)为待修正的异常样本数据;该试验台架上同时测得了n个轴承的实验样本数据,对φn中的样本数据,基于同一试验台架上多个轴承的其他样本数据进行联合修正,得到φ'n,并使得φ'n中的历史数据与当前数据相一致:
Figure FDA0002884165380000042
式中:x″(ti)为异常样本x(ti)在第2次修正之后的结果,x(ti,j)为第j个轴承在第i个时刻点样本统计量数值;
步骤2-4,综合步骤2-2和步骤2-3得到滚动轴承异常数据样本的最终修正结果为:
Figure FDA0002884165380000043
采用相同的步骤可以实施对实验台架上其他轴承在ti时刻的异常样本数据b(ti),c(ti)进行修正,由此即可完成对轴承性能测试中所有异常监测样本的修正。
CN202110008755.6A 2021-01-05 2021-01-05 一种滚动轴承性能测试台和异常测试样本的识别与修正方法 Active CN112816216B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110008755.6A CN112816216B (zh) 2021-01-05 2021-01-05 一种滚动轴承性能测试台和异常测试样本的识别与修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110008755.6A CN112816216B (zh) 2021-01-05 2021-01-05 一种滚动轴承性能测试台和异常测试样本的识别与修正方法

Publications (2)

Publication Number Publication Date
CN112816216A true CN112816216A (zh) 2021-05-18
CN112816216B CN112816216B (zh) 2022-08-16

Family

ID=75857454

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110008755.6A Active CN112816216B (zh) 2021-01-05 2021-01-05 一种滚动轴承性能测试台和异常测试样本的识别与修正方法

Country Status (1)

Country Link
CN (1) CN112816216B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117113009A (zh) * 2023-10-23 2023-11-24 合肥亚明汽车部件有限公司 一种数字化工厂设备运行风险预警方法
CN117713914A (zh) * 2024-02-06 2024-03-15 北斗天汇(北京)科技有限公司 一种北斗卫星智能遥测遥控通信一体化终端系统

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060167659A1 (en) * 2002-08-30 2006-07-27 Nsk Ltd. Method and device for monitoring status of mechanical equipment and abnormality diagnosing device
CN102519443A (zh) * 2011-11-26 2012-06-27 东南大学 车载微机械陀螺仪异常测量数据的识别与修正方法
CN104766175A (zh) * 2015-04-16 2015-07-08 东南大学 一种基于时间序列分析的电力系统异常数据辨识与修正方法
CN106548270A (zh) * 2016-09-30 2017-03-29 许昌许继软件技术有限公司 一种光伏电站功率异常数据辨识方法及装置
CN107024353A (zh) * 2017-05-11 2017-08-08 浙江省机电设计研究院有限公司 一种密封轴承综合性能模拟试验机
WO2017145222A1 (ja) * 2016-02-22 2017-08-31 株式会社日立製作所 軸受劣化診断装置、軸受劣化診断方法及び軸受劣化診断システム
CN109084981A (zh) * 2018-10-22 2018-12-25 中国矿业大学 一种轴承冲击摩擦磨损试验机
CN111238815A (zh) * 2020-02-26 2020-06-05 中国民航大学 一种样本不平衡下基于数据增强的轴承故障识别方法
CN111896254A (zh) * 2020-08-10 2020-11-06 山东大学 一种变速变载大型滚动轴承故障预测系统及方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060167659A1 (en) * 2002-08-30 2006-07-27 Nsk Ltd. Method and device for monitoring status of mechanical equipment and abnormality diagnosing device
CN102519443A (zh) * 2011-11-26 2012-06-27 东南大学 车载微机械陀螺仪异常测量数据的识别与修正方法
CN104766175A (zh) * 2015-04-16 2015-07-08 东南大学 一种基于时间序列分析的电力系统异常数据辨识与修正方法
WO2017145222A1 (ja) * 2016-02-22 2017-08-31 株式会社日立製作所 軸受劣化診断装置、軸受劣化診断方法及び軸受劣化診断システム
CN106548270A (zh) * 2016-09-30 2017-03-29 许昌许继软件技术有限公司 一种光伏电站功率异常数据辨识方法及装置
CN107024353A (zh) * 2017-05-11 2017-08-08 浙江省机电设计研究院有限公司 一种密封轴承综合性能模拟试验机
CN109084981A (zh) * 2018-10-22 2018-12-25 中国矿业大学 一种轴承冲击摩擦磨损试验机
CN111238815A (zh) * 2020-02-26 2020-06-05 中国民航大学 一种样本不平衡下基于数据增强的轴承故障识别方法
CN111896254A (zh) * 2020-08-10 2020-11-06 山东大学 一种变速变载大型滚动轴承故障预测系统及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
陈亚红 等: "短期电力负荷预报中几种异常数据的处理", 《东北电力学院学报》 *
陈浩 等: "基于三次样条函数的激光雷达数据可视化插值法", 《仪器仪表学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117113009A (zh) * 2023-10-23 2023-11-24 合肥亚明汽车部件有限公司 一种数字化工厂设备运行风险预警方法
CN117113009B (zh) * 2023-10-23 2024-01-16 合肥亚明汽车部件有限公司 一种数字化工厂设备运行风险预警方法
CN117713914A (zh) * 2024-02-06 2024-03-15 北斗天汇(北京)科技有限公司 一种北斗卫星智能遥测遥控通信一体化终端系统

Also Published As

Publication number Publication date
CN112816216B (zh) 2022-08-16

Similar Documents

Publication Publication Date Title
CN111051852B (zh) 多核传感器系统及其隔离和恢复的方法
CN1103946C (zh) 分析测量值的方法以及实施此方法的测量值分析器
CN112816216B (zh) 一种滚动轴承性能测试台和异常测试样本的识别与修正方法
CN110363339B (zh) 一种基于电机参数进行预测性维护的方法与系统
CN114615702B (zh) 基于5g通信网络故障检测诊断系统
CN111323158B (zh) 一种联轴螺栓的实时应力监测系统及安全系数检测方法
CN115935243B (zh) 一种基于数据处理的故障分析方法
CN112131797B (zh) 一种基于应力分析的主轴轴承寿命预测和可靠性评估方法
CN116295948A (zh) 大温差环境下工业温度传感器的异常检测方法、系统及存储介质
CN110530639A (zh) 一种高速列车轴箱轴承故障诊断方法
CN113358958A (zh) 一种伺服驱动器加速寿命试验装置及测试方法
CN114035086B (zh) 基于信号处理的电池组多故障诊断方法
CN113378425B (zh) 基于雨流计数法的汽车换热器产品疲劳寿命测试方法
CN102129029A (zh) 一种基于低频噪声测量的模拟电路故障诊断方法
CN113405590A (zh) 轨道车辆关键部件状态测试装置、系统及方法
CN117091237A (zh) 一种用于轨道列车的空调设备检测方法及系统
KR20180042897A (ko) 압축센싱 기반 구조물 상태진단 시스템 및 그 방법
CN113532776B (zh) 一种发电机轴瓦瓦枕绝缘垫失效诊断方法及系统
CN116878831A (zh) 一种螺栓松动监测方法及监测装置
CN115455596A (zh) 基于粒子群算法的聚氨酯胶密封件贮存可靠性评估方法
CN106403880B (zh) 一种压缩机转子和定子之间的间隙检测方法和装置
KR20070087139A (ko) 진공펌프의 고장 방지 및 예보유지를 위한 정밀진단 기법과그 구현 시스템.
CN116990014B (zh) 一种齿轮箱性能测试系统
CN109388829B (zh) 一种电子产品寿命测算方法
CN112485009A (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