欢迎来到知嘟嘟! 联系电话:13095918853 卖家免费入驻,海量在线求购! 卖家免费入驻,海量在线求购!
知嘟嘟
我要发布
联系电话:13095918853
知嘟嘟经纪人
收藏
专利号: 2021110045367
申请人: 西安理工大学
专利类型:发明专利
专利状态:已下证
专利领域: 测量;测试
更新日期:2024-01-05
缴费截止日期: 暂无
价格&联系人
年费信息
委托购买

摘要:

权利要求书:

1.基于多尺度散布熵构造阈值的故障报警方法,其特征在于,具体按照以下步骤实施:

步骤1、导入滚动轴承整个生命周期从正常状态到发生故障到最终完全失效的带工况标记的振动信号;

步骤2、将步骤1所述的振动数据划分为G个数据段,每个数据段包含N个数据,对所有数据段进行多尺度散布熵MDE值计算,将所有数据段的MDE值存入长度为G的一维向量W,选取一个MDE值记为S作为报警阈值设定在系统中,因为故障工况的MDE值大于正常工况的MDE值,所以报警阈值S介于故障工况MDE值与正常工况MDE值之间;

步骤3、使用自回归模型对轴承未来一段时间的振动信号进行预测,当预测出的振动信号段的MDE值大于报警阈值S则报警,预测出的振动信号段的MDE值小于报警阈值S则进行下一阶段的MDE值计算;

步骤4、导入实际的正常工况和故障工况的轴承振动信号段进行步骤2的多尺度散布熵MDE值计算,根据计算出来的正常工况MDE值和故障工况MDE值选取具体的报警阈值S。

2.根据权利要求1所述的基于多尺度散布熵构造阈值的故障报警方法,其特征在于,所述步骤2具体如下:步骤2.1、将步骤1的轴承振动信号划分为长度为N的时间序列X,X={x1,x2,...,xi,...,xN},i=1,2,...,N,xi为时间序列X中的一个元素,将时间序列X重组成尺度因子为τ的多尺度时间序列Y, 的表达式如公式(1):式中, 表示重组后的尺度因子为τ的多尺度时间序列Y中的一个元素,重组后的多尺度时间序列Y的长度为N‑τ+1,τ表示多尺度分析中的尺度因子,为一个正整数,当τ=1时,Y为原始时间序列;

步骤2.2、从相空间嵌入理论出发,利用嵌入维数m,将步骤2.1中的重组后的多尺度时间序列Y重构为一系列时间轨道 t=1,2,...,N‑τ‑m+2,m<N‑τ‑m+

2,重构的相空间包含大量与原始相空间具有相同吸引子的维度,

Z(m)为将步骤2.1中的重组后的多尺度时间序列Y利用嵌入维数m重构后的时间轨道矩阵;

步骤2.3、计算相邻轨道之间的余弦相似度,得到一系列的余弦相似度d=(d1,d2,...,df,...,dN‑τ‑m+1),f=1,2,...,N‑τ‑m+1,df为第f‑1轨道和第f轨道之间的余弦相似度,定义为:余弦相似度d的取值范围是[‑1,1],余弦相似度d的取值的绝对值趋近于1表示两个轨道之间相似、可预测或周期性的动态变化,相反,余弦相似度d的取值的绝对值趋近于0代表多样性、随机或混沌的动态行为;

步骤2.4、将余弦相似度d的取值范围[‑1,1]平均划分为n个区间,表示为(I1,I2,...,Ik,...,In),k=1,2,...,n,Ik为将[‑1,1]区间n均等分后的第k个区间Ik,然后,计算状态概率向量P=(P1,P2,...Pk,...,Pn),Pk表示统计一系列的余弦相似度d的取值落在第k个区间Ik的概率,其中步骤2.5、多尺度散布熵MDE值根据获得的状态概率Pk计算:

步骤2.6、构建一个长度为G的一维向量W,将所有信号段的MDE值计算完成后存入向量W,以横坐标为[1:G]、纵坐标为向量W,绘制出MDE值分布图,理论上正常工况的MDE值小于故障工况的MDE值,对比观察MDE值分布图和带工况标记的原信号时域图,选取介于正常工况的MDE值与故障工况的MDE值之间的MDE值作为报警阈值记为S,将报警阈值S设定到系统中,一旦MDE值达到报警阈值S,系统将发出报警。