一种深沟球滚动轴承内圈故障劣化状态识别方法及系统专利检索-轴承内圈轴承套圈滚动轴承轴承轴承类专利检索查询-专利查询网

积极推动地理标志专门立法

2022-03-10 地理标志，立法，知识产权
保护知识产权是对创新最大的激励

2022-03-10 保护知识产权，创新，激励
谢商华：加快制定知识产权基本法

2022-03-10 知识产权基本法
擦亮“双奥之城”品牌

2022-03-10 双奥，知识产权
让冰雪运动“热”力全开

2022-03-10 冰雪运动，知识产权
携手共奋进　走好强国路

2022-03-10 强国，知识产权
坚持创新引领　方能稳中求进

2022-03-10 创新，稳中求进，知识产权
答好“两张卷” 奋进新征程

2022-03-10 知识产权
专家解读政府工作报告中的创新和知识产权相关部署

2022-03-10 政府工作报告，创新，知识产权
今年政府工作报告指出：加强知识产权保护和运用

2022-03-10 政府工作报告，知识产权保护

一种深沟球滚动轴承内圈故障劣化状态识别方法及系统

阅读：59发布：2021-03-02

IPRDB可以提供一种深沟球滚动轴承内圈故障劣化状态识别方法及系统专利检索，专利查询，专利分析的服务。并且本发明公开一种深沟球滚动轴承内圈故障劣化状态识别方法及系统。本发明基于滚动轴承内圈故障劣化过程中采集的振动加速度时间序列信号，通过分别对振动响应序列信号采用非线性和线性自回归各态历经时间序列模型进行辨识，再计算两种模型之间响应的最小二乘积分误差得到非线性估计值，能够对机械系统退化状态特征进行评估，具有计算步骤简单、容易实现、结果准确等优点。，下面是一种深沟球滚动轴承内圈故障劣化状态识别方法及系统专利的具体信息内容。

权利要求

1.一种深沟球滚动轴承内圈故障劣化状态识别方法，其特征在于，包括：在含有内圈故障深沟球滚动轴承稳定转速运行状态采集设定时长的振动加速度时间序列信号；

对所述振动加速度时间序列信号进行去噪处理，得到去噪后振动响应序列信号；

对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

根据所述辨识非线性时间序列模型所用的去噪振动响应时间序列信号，建立对应所述辨识非线性时间序列模型结构的最优线性时间序列模型，所述最优线性时间序列模型为单输入单输出；

根据所述最优线性时间序列模型和所述非线性时间序列模型，确定非线性估计值；

根据所述非线性估计值识别深沟球滚动轴承内圈故障劣化状态。

2.根据权利要求1所述的深沟球滚动轴承内圈故障劣化状态识别方法，其特征在于，所述对所述振动加速度时间序列信号进行去噪处理，得到去噪后振动响应序列信号，具体包括：对所述振动加速度时间序列信号采用小波等滤波算法进行去噪处理，得到去噪后振动响应序列信号。

3.根据权利要求1所述的深沟球滚动轴承内圈故障劣化状态识别方法，其特征在于，所述对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型，具体包括：采用非线性自回归各态历经时间序列模型对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

所述非线性自回归各态历经时间序列模型公式如下：

y(t)＝f(y(t-1),y(t-2),…,y(t-na),u(t),u(t-1),…,u(t-nb+1))其中，f表示非线性函数，y(t)表示去噪后的输出振动时间序列信号，u(t)表示去噪后的输入时间序列信号，na表示输出时间序列回归量的个数，nb表示输入时间序列回归量的个数。

4.根据权利要求1所述的深沟球滚动轴承内圈故障劣化状态识别方法，其特征在于，所述根据所述辨识非线性时间序列模型所用的去噪振动响应时间序列信号，建立对应所述辨识非线性时间序列模型结构的最优线性时间序列模型，具体包括：采用线性自回归各态历经时间序列模型，建立与所述非线性时间序列模型对应结构的线性单输入单输出时间序列模型；

所述线性自回归各态历经时间序列模型公式如下：

y(t)+a1y(t-1)+a2y(t-2)+…+anay(t-na)＝b1u(t)+b2u(t-1)+…+bnbu(t-nb+1)+e(t)其中，y(t)表示去噪后的输出振动时间序列信号，u(t)表示去噪后的输入时间序列信号，na表示输出时间序列回归量的个数，nb表示输入时间序列回归量的个数，a1,a2,…,ana,b1,b2,…,bnb表示常数系数，e(t)表示误差函数。

5.根据权利要求1所述的深沟球滚动轴承内圈故障劣化状态识别方法，其特征在于，所述根据所述最优线性时间序列模型和所述非线性时间序列模型，确定非线性估计值，具体包括：根据所述最优线性时间序列模型和所述非线性时间序列模型采用最小二乘法，计算所述最优线性时间序列模型和所述非线性时间序列模型的最小二乘积分误差；

根据所述最小二乘积分误差，确定非线性估计值。

6.一种深沟球滚动轴承内圈故障劣化状态识别系统，其特征在于，包括：信号采集模块，用于在含有内圈故障深沟球滚动轴承稳定转速运行状态采集设定时长的振动加速度时间序列信号；

去噪处理模块，用于对所述振动加速度时间序列信号进行去噪处理，得到去噪后振动响应序列信号；

非线性模型建立模块，用于对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

线性模型建立模块，用于根据所述辨识非线性时间序列模型所用的去噪振动响应时间序列信号，建立对应所述辨识非线性时间序列模型结构的最优线性时间序列模型，所述最优线性时间序列模型为单输入单输出；

估计值确定模块，用于根据所述最优线性时间序列模型和所述非线性时间序列模型，确定非线性估计值；

状态识别模块，用于根据所述非线性估计值识别深沟球滚动轴承内圈故障劣化状态。

7.根据权利要求6所述的深沟球滚动轴承内圈故障劣化状态识别系统，其特征在于，所述去噪处理模块，具体包括：去噪处理单元，用于对所述振动加速度时间序列信号采用小波等滤波算法进行去噪处理，得到去噪后振动响应序列信号。

8.根据权利要求6所述的深沟球滚动轴承内圈故障劣化状态识别系统，其特征在于，所述辨识非线性模型建立模块，具体包括：非线性模型建立单元，用于采用非线性自回归各态历经时间序列模型对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

所述非线性自回归各态历经时间序列模型公式如下：

9.根据权利要求6所述的深沟球滚动轴承内圈故障劣化状态识别系统，其特征在于，所述线性模型建立模块，具体包括：线性模型建立单元，用于采用线性自回归各态历经时间序列模型，建立与所述非线性时间序列模型对应结构的线性单输入单输出时间序列模型；

所述线性自回归各态历经时间序列模型公式如下：

10.根据权利要求6所述的深沟球滚动轴承内圈故障劣化状态识别系统，其特征在于，所述估计值确定模块，具体包括：估计值确定单元，用于根据所述最优线性时间序列模型和所述非线性时间序列模型采用最小二乘法，计算所述最优线性时间序列模型和所述非线性时间序列模型的最小二乘积分误差；

根据所述最小二乘积分误差，确定非线性估计值。

说明书全文

一种深沟球滚动轴承内圈故障劣化状态识别方法及系统

技术领域

[0001] 本发明涉及机械系统状态监测与评估技术领域，特别是涉及一种深沟球滚动轴承内圈故障劣化状态识别方法及系统。

背景技术

[0002] 滚动轴承是现代机械系统中大量使用的机械零件，其运行状态是否正常往往直接影响到整台机器的性能，例如精度、可靠性及寿命等。为了实现预知维修保证滚动轴承正常运行，需要提早地发现轴承故障劣化程度，所以对滚动轴承内圈故障劣化状态识别研究具有重要的价值。据统计，滚动轴承故障在机械设备故障中占有很大的比例，而各种轴承失效中内圈或外圈裂纹占失效总数的比例很高。在正常工作条件下，由于受到载荷、安装、润滑状态或轴承工作条件恶劣等因素的影响，滚动轴承运转一段时间后将会产生故障，大部分可归结为轴承元件表面的劣化，进而使振动加剧。初始阶段故障程度较为轻微，继续运行中故障会进一步劣化，使其运转状态逐渐恶化以至完全失效。

[0003] 目前滚动轴承故障劣化识别方法中，第一种是经典的参数统计估计方法，这种方法有很大的局限性。第二种方法是智能识别方法，如人工神经网络和支持向量机。该方法克服了传统参数估计方法的困难，能自适应地学习和调整网络的大小。但是它需要大量的典型故障数据样本或经验知识，而在工程实际当中，一般很难获得大量典型的故障样本。目前滚动轴承的运行状态监测与故障诊断技术在诊断信息的处理与分析、故障类型的识别等方面已经取得了很大的进步，但是对滚动轴承的故障程度估计以及劣化状态和趋势预测方面的研究所做的工作还较少。现有的振动法退化状态识别技术一般基于线性模型的假设，或者把非线性模型线性近似化，这种方法在某一些具有弱动力学行为的机械系统上应用是可行的，但对具有强非线性动力学行为的滚动轴承出现较大的误差。而随着机械系统高速化、轻量化等趋势的发展，含故障滚动轴承运行时将出现非线性动力学行为或者其动力学行为非线性程度变得更高，基于非线性理论对滚动轴承劣化状态进行识别是未来的趋势之一。

发明内容

[0004] 本发明的目的是提供一种深沟球滚动轴承内圈故障劣化状态识别方法及系统，能够简便、快捷且精确的实现深沟球滚动轴承内圈故障劣化状态识别。

[0005] 为实现上述目的，本发明提供了如下方案：

[0006] 一种深沟球滚动轴承内圈故障劣化状态识别方法，包括：

[0007] 在含有内圈故障深沟球滚动轴承稳定转速运行状态采集设定时长的振动加速度时间序列信号；

[0008] 对所述振动加速度时间序列信号进行去噪处理，得到去噪后振动响应序列信号；

[0009] 对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

[0010] 根据所述辨识非线性时间序列模型所用的去噪振动响应时间序列信号，建立对应所述辨识非线性时间序列模型结构的最优线性时间序列模型，所述最优线性时间序列模型为单输入单输出；

[0011] 根据所述最优线性时间序列模型和所述非线性时间序列模型，确定非线性估计值；

[0012] 根据所述非线性估计值识别深沟球滚动轴承内圈故障劣化状态。

[0013] 可选的，所述对所述振动加速度时间序列信号进行去噪处理，得到去噪后振动响应序列信号，具体包括：

[0014] 对所述振动加速度时间序列信号采用小波等滤波算法进行去噪处理，得到去噪后振动响应序列信号。

[0015] 可选的，所述对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型，具体包括：

[0016] 采用非线性自回归各态历经时间序列模型对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

[0017] 所述非线性自回归各态历经时间序列模型公式如下：

[0018] y(t)＝f(y(t-1),y(t-2),…,y(t-na),u(t),u(t-1),…,u(t-nb+1))[0019] 其中，f表示非线性函数，y(t)表示去噪后的输出振动时间序列信号，u(t)表示去噪后的输入时间序列信号，na表示输出时间序列回归量的个数，nb表示输入时间序列回归量的个数。

[0020] 可选的，根据所述辨识非线性时间序列模型所用的去噪振动响应时间序列信号，建立对应所述辨识非线性时间序列模型结构的最优线性时间序列模型，具体包括：

[0021] 采用线性自回归各态历经时间序列模型，建立与所述非线性时间序列模型对应结构的线性单输入单输出时间序列模型；

[0022] 所述线性自回归各态历经时间序列模型公式如下：

[0023] y(t)+a1y(t-1)+a2y(t-2)+…+anay(t-na)＝b1u(t)+b2u(t-1)+…+bnbu(t-nb+1)+e(t)

[0024] 其中，y(t)表示去噪后的输出振动时间序列信号，u(t)表示去噪后的输入时间序列信号，na表示输出时间序列回归量的个数，nb表示输入时间序列回归量的个数，a1,a2,…,ana,b1,b2,…,bnb表示常数系数，e(t)表示误差函数。

[0025] 可选的，所述根据所述最优线性时间序列模型和所述非线性时间序列模型，确定非线性估计值，具体包括：

[0026] 根据所述最优线性时间序列模型和所述非线性时间序列模型采用最小二乘法，计算所述最优线性时间序列模型和所述非线性时间序列模型的最小二乘积分误差；

[0027] 根据所述最小二乘积分误差，确定非线性估计值。

[0028] 一种深沟球滚动轴承内圈故障劣化状态识别系统，包括：

[0029] 信号采集模块，用于在含有内圈故障深沟球滚动轴承稳定转速运行状态采集设定时长的振动加速度时间序列信号；

[0030] 去噪处理模块，用于对所述振动加速度时间序列信号进行去噪处理，得到去噪后振动响应序列信号；

[0031] 非线性模型建立模块，用于对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

[0032] 线性模型建立模块，用于根据所述辨识非线性时间序列模型所用的去噪振动响应时间序列信号，建立对应所述辨识非线性时间序列模型结构的最优线性时间序列模型，所述最优线性时间序列模型为单输入单输出；

[0033] 估计值确定模块，用于根据所述最优线性时间序列模型和所述非线性时间序列模型，确定非线性估计值；

[0034] 状态识别模块，用于根据所述非线性估计值识别深沟球滚动轴承内圈故障劣化状态。

[0035] 可选的，所述去噪处理模块，具体包括：

[0036] 去噪处理单元，用于对所述振动加速度时间序列信号采用小波等滤波算法进行去噪处理，得到去噪后振动响应序列信号。

[0037] 可选的，所述非线性模型建立模块，具体包括：

[0038] 非线性模型建立单元，用于采用非线性自回归各态历经时间序列模型对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

[0039] 所述非线性自回归各态历经时间序列模型公式如下：

[0040] y(t)＝f(y(t-1),y(t-2),…,y(t-na),u(t),u(t-1),…,u(t-nb+1))[0041] 其中，f表示非线性函数，y(t)表示去噪后的输出振动时间序列信号，u(t)表示去噪后的输入时间序列信号，na表示输出时间序列回归量的个数，nb表示输入时间序列回归量的个数。

[0042] 可选的，所述线性模型建立模块，具体包括：

[0043] 线性模型建立单元，用于采用线性自回归各态历经时间序列模型，根据所述辨识非线性时间序列模型所用的去噪振动响应时间序列信号，建立与所述非线性时间序列模型对应结构的线性单输入单输出时间序列模型；

[0044] 所述线性自回归各态历经时间序列模型公式如下：

[0045] y(t)+a1y(t-1)+a2y(t-2)+…+anay(t-na)＝b1u(t)+b2u(t-1)+…+bnbu(t-nb+1)+e(t)

[0046] 其中，y(t)表示去噪后的输出振动时间序列信号，u(t)表示去噪后的输入时间序列信号，na表示输出时间序列回归量的个数，nb表示输入时间序列回归量的个数，a1,a2,…,ana,b1,b2,…,bnb表示常数系数，e(t)表示误差函数。

[0047] 可选的，所述估计值确定模块，具体包括：

[0048] 估计值确定单元，用于根据所述最优线性时间序列模型和所述非线性时间序列模型采用最小二乘法，计算所述最优线性时间序列模型和所述非线性时间序列模型的最小二乘积分误差；

[0049] 根据所述最小二乘积分误差，确定非线性估计值。

[0050] 根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明提供一种深沟球滚动轴承内圈故障劣化状态识别方法，基于滚动轴承内圈故障劣化过程中采集的振动加速度时间序列信号，通过分别对振动响应序列信号采用非线性和线性自回归各态历经时间序列模型进行辨识，再计算两种模型之间响应的最小二乘积分误差得到非线性估计值，能够对机械系统退化状态特征进行评估，具有计算步骤简单、容易实现、结果准确等优点。

附图说明

[0051] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

[0052] 图1为本发明实施例深沟球滚动轴承内圈故障劣化状态识别方法流程图；

[0053] 图2为本发明的振动响应序列非线性估计原理图；

[0054] 图3为本发明中实施例的振动响应序列排列熵计算图(信号长度10240)；

[0055] 图4为本发明中实施例的振动响应序列排列熵计算图(信号长度4096)；

[0056] 图5为本发明中实施例的振动响应序列非线性估计结果图(信号长度10240)；

[0057] 图6为本发明中实施例的振动响应序列非线性估计结果图(信号长度4096)；

[0058] 图7为本发明实施例深沟球滚动轴承内圈故障劣化状态识别系统结构图。

具体实施方式

[0059] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0060] 本发明的目的是提供一种深沟球滚动轴承内圈故障劣化状态识别方法及系统，能够简便、快捷且精确的实现深沟球滚动轴承内圈故障劣化状态识别。

[0061] 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

[0062] 图1为本发明实施例深沟球滚动轴承内圈故障劣化状态识别方法流程图。如图1所示，一种深沟球滚动轴承内圈故障劣化状态识别方法，包括：

[0063] 步骤101：在含有内圈故障深沟球滚动轴承稳定转速运行状态采集设定时长的振动加速度时间序列信号；采用接触式或者非接触式测振方式得到滚动轴承外圈可测位置的振动信号；

[0064] 步骤102：对所述振动加速度时间序列信号进行去噪处理，得到去噪后振动响应序列信号；

[0065] 步骤103：对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

[0066] 步骤104：根据所述辨识非线性时间序列模型所用的去噪振动响应时间序列信号，建立对应所述辨识非线性时间序列模型结构的最优线性时间序列模型，所述最优线性时间序列模型为单输入单输出；

[0067] 步骤105：根据所述最优线性时间序列模型和所述非线性时间序列模型，确定非线性估计值；

[0068] 步骤106：根据所述非线性估计值识别深沟球滚动轴承内圈故障劣化状态。

[0069] 步骤102，具体包括：

[0070] 对所述振动加速度时间序列信号采用小波等滤波算法进行去噪处理，得到去噪后振动响应序列信号。

[0071] 步骤103，具体包括：

[0072] 采用非线性自回归各态历经时间序列模型对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

[0073] 所述非线性自回归各态历经时间序列模型公式如下：

[0074] y(t)＝f(y(t-1),y(t-2),…,y(t-na),u(t),u(t-1),…,u(t-nb+1))[0075] 其中，f表示非线性函数，y(t)表示去噪后的输出振动时间序列信号，u(t)表示去噪后的输入时间序列信号，na表示输出时间序列回归量的个数，nb表示输入时间序列回归量的个数。

[0076] 步骤104，具体包括：

[0077] 根据所述辨识非线性时间序列模型所用的去噪振动响应时间序列信号，采用线性自回归各态历经时间序列模型，建立与所述非线性时间序列模型对应结构的线性单输入单输出时间序列模型；

[0078] 所述线性自回归各态历经时间序列模型公式如下：

[0079] y(t)+a1y(t-1)+a2y(t-2)+…+anay(t-na)＝b1u(t)+b2u(t-1)+…+bnbu(t-nb+1)+e(t)

[0080] 其中，y(t)表示去噪后的输出振动时间序列信号，u(t)表示去噪后的输入时间序列信号，na表示输出时间序列回归量的个数，nb表示输入时间序列回归量的个数，a1,a2,…,ana,b1,b2,…,bnb表示常数系数，e(t)表示误差函数。

[0081] 步骤105，具体包括：

[0082] 根据所述最优线性时间序列模型和所述非线性时间序列模型采用最小二乘法，计算所述最优线性时间序列模型和所述非线性时间序列模型的最小二乘积分误差；

[0083] 根据所述最小二乘积分误差，确定非线性估计值。

[0084] 本发明公开的深沟球滚动轴承内圈故障劣化状态识别方法，基于滚动轴承内圈故障劣化过程中采集的振动加速度时间序列信号，通过分别对振动响应序列信号采用非线性和线性自回归各态历经时间序列模型进行辨识，再计算两种模型之间响应的最小二乘积分误差得到非线性估计值，能够对机械系统退化状态特征进行评估，具有计算步骤简单、容易实现、结果准确等优点。

[0085] 时间序列非线性估计的定义：

[0086] 动力学非线性程度估计是将非线性系统N(u)的动力学行为和某个作为标准的线性系统L(u)的动力学行为差别进行度量的一种评价指标(包含相对误差︱N(u)-L(u)/L(u)︱和绝对误差︱N(u)-L(u)︱)。其中N(u)，L(u)分别代表非线性系统和线性系统的动力学行为，可采用振动信号，加速度信号，速度信号等来对其动力学行为进行描述。而振动信号，加速度信号，速度信号均为时间序列信号，非线性系统的动力学行为非线性程度的变化也体现在上述时间序列信号中，因此本发明是通过估计时间序列信号中非线性程度的变化来实现对机械系统退化状态导致系统动力学行为非线性程度变化的估计。

[0087] 振动加速度时间序列信号采集与处理：

[0088] 为能够实现本方法对滚动轴承内圈故障劣化状态的预测及说明本方法的可行性，以滚动轴承内圈故障振动测试实验为基础，在不同的转速和内圈故障长度设置下采集振动加速度时间序列信号，每组时间序列信号的长度为10240，共采集了四种内圈故障长度下的外圈振动加速度信号，将上述信号采集以后形成mat文件，在Matlab里面采用小波方法进行去噪处理。图3为本发明中实施例的振动响应序列排列熵计算图(信号长度10240)。图4为本发明中实施例的振动响应序列排列熵计算图(信号长度4096)。图5为本发明中实施例的振动响应序列非线性估计结果图(信号长度10240)。图6为本发明中实施例的振动响应序列非线性估计结果图(信号长度4096)。

[0089] 非线性自回归各态历经时间序列模型辨识：

[0090] (1)对采集的振动加速度时间序列信号v1,v2,…,vT(T为时间序列长度)采用小波方法进行去噪处理；

[0091] (2)令时间序列信号v4,v5,…,vT为输出时间序列y(t)，时间序列信号v2,v3,…,vT-2为输入时间序列u(t)；

[0092] (3)选择3个输出时间序列回归量，2个输入时间序列回归量；

[0093] (4)采用树分解算法确定非线性函数f的结构，基于输出时间序列y(t)和输入时间序列u(t)，采用Matlab自带函数isnlarx辨识获得非线性自回归各态历经时间序列模型如下：

[0094] y(t)＝f(y(t-1),y(t-2),y(t-3),u(t),u(t-1))

[0095] 线性自回归各态历经时间序列模型辨识：

[0096] (1)对采集的振动加速度时间序列信号v1,v2,…,vT(T为时间序列长度)采用小波方法进行去噪处理；

[0097] (2)令时间序列信号v4,v5,…,vT为输出时间序列y(t)，时间序列信号v2,v3,…,vT-2为输入时间序列u(t)；

[0098] (3)选择3个输出时间序列回归量，2个输入时间序列回归量；

[0099] (4)建立线性自回归各态历经时间序列模型结构如下：

[0100] y(t)+a1y(t-1)+a2y(t-2)+a3y(t-3)＝b1u(t)+b2u(t-1)+e(t) (5)[0101] 其中，y(t)表示去噪后的输出振动时间序列信号，u(t)表示去噪后的输入时间序列信号，a1,a2,a3,b1,b2表示常数系数，e(t)表示误差函数。

[0102] (5)基于输出时间序列y(t)和输入时间序列u(t)采用最小二乘法辨识或者建立状态方程由矩阵分解算法获得线性自回归各态历经时间序列模型参数值。

[0103] 不同内圈故障长度下非线性估计值的计算：

[0104] 将不同长度内圈故障下采集到的振动加速度时间序列信号辨识得到的非线性自回归各态历经时间序列模型与线性自回归各态历经时间序列模型的响应误差进行计算，并在时间序列长度内进行积分求得不同时刻点的非线性估计值，计算示意图如图2所示。图2为本发明的振动响应序列非线性估计原理图。

[0105] 图7为本发明实施例深沟球滚动轴承内圈故障劣化状态识别系统结构图。如图7所示，一种深沟球滚动轴承内圈故障劣化状态识别系统，包括：

[0106] 信号采集模块201，用于在含有内圈故障深沟球滚动轴承稳定转速运行状态采集设定时长的振动加速度时间序列信号；

[0107] 去噪处理模块202，用于对所述振动加速度时间序列信号进行去噪处理，得到去噪后振动响应序列信号；

[0108] 非线性模型建立模块203，用于对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

[0109] 线性模型建立模块204，用于根据所述辨识非线性时间序列模型所用的去噪振动响应时间序列信号，，建立对应所述辨识非线性时间序列模型结构的最优线性时间序列模型，所述最优线性时间序列模型为单输入单输出；

[0110] 估计值确定模块205，用于根据所述最优线性时间序列模型和所述非线性时间序列模型，确定非线性估计值；

[0111] 状态识别模块206，用于根据所述非线性估计值识别深沟球滚动轴承内圈故障劣化状态。

[0112] 所述去噪处理模块202，具体包括：

[0113] 去噪处理单元，用于对所述振动加速度时间序列信号采用小波等滤波算法进行去噪处理，得到去噪后振动响应序列信号。

[0114] 所述非线性模型建立模块203，具体包括：

[0115] 非线性模型建立单元，用于采用非线性自回归各态历经时间序列模型对所述去噪后振动响应序列信号进行辨识，得到描述所述振动信号的辨识非线性时间序列模型；

[0116] 所述非线性自回归各态历经时间序列模型公式如下：

[0117] y(t)＝f(y(t-1),y(t-2),…,y(t-na),u(t),u(t-1),…,u(t-nb+1))[0118] 其中，f表示非线性函数，y(t)表示去噪后的输出振动时间序列信号，u(t)表示去噪后的输入时间序列信号，na表示输出时间序列回归量的个数，nb表示输入时间序列回归量的个数。

[0119] 所述线性模型建立模块204，具体包括：

[0120] 线性模型建立单元，用于采用线性自回归各态历经时间序列模型，根据所述辨识非线性时间序列模型所用的去噪振动响应时间序列信号，建立与所述非线性时间序列模型对应结构的线性单输入单输出时间序列模型；

[0121] 所述线性自回归各态历经时间序列模型公式如下：

[0122] y(t)+a1y(t-1)+a2y(t-2)+…+anay(t-na)＝b1u(t)+b2u(t-1)+…+bnbu(t-nb+1)+e(t)

[0123] 其中，y(t)表示去噪后的输出振动时间序列信号，u(t)表示去噪后的输入时间序列信号，na表示输出时间序列回归量的个数，nb表示输入时间序列回归量的个数，a1,a2,…,ana,b1,b2,…,bnb表示常数系数，e(t)表示误差函数。

[0124] 所述估计值确定模块205，具体包括：

[0125] 估计值确定单元，用于根据所述最优线性时间序列模型和所述非线性时间序列模型采用最小二乘法，计算所述最优线性时间序列模型和所述非线性时间序列模型的最小二乘积分误差；

[0126] 根据所述最小二乘积分误差，确定非线性估计值。

[0127] 本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

[0128] 本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

标题	发布/更新时间	阅读量
滚子轴承的内圈-专利编号CN102239338B	2020-05-12	1013
防开裂轴承内圈-专利编号CN102829085A	2020-05-12	95
一种轴承内圈-专利编号CN102537079A	2020-05-13	752
轴承及其内圈-专利编号CN105370734B	2020-05-11	978
一种轴承内圈-专利编号CN106015345A	2020-05-11	669
轴承内圈冷碾机-专利编号CN102886475B	2020-05-11	352
一种轴承内圈-专利编号CN104696374A	2020-05-11	212
一种轴承内圈-专利编号CN103322061A	2020-05-12	759
一种轴承内圈-专利编号CN102022435A	2020-05-13	675
滚子轴承的内圈-专利编号CN102239338A	2020-05-13	605

一种深沟球滚动轴承内圈故障劣化状态识别方法及系统

一种深沟球滚动轴承内圈故障劣化状态识别方法及系统

技术领域

背景技术

发明内容

附图说明

具体实施方式

IPRDB

热门服务

关于我们

友情链接

联系方式