CN102449445A - 用于分析具有旋转部件的机器振动的设备和方法 - Google Patents
用于分析具有旋转部件的机器振动的设备和方法 Download PDFInfo
- Publication number
- CN102449445A CN102449445A CN2010800197378A CN201080019737A CN102449445A CN 102449445 A CN102449445 A CN 102449445A CN 2010800197378 A CN2010800197378 A CN 2010800197378A CN 201080019737 A CN201080019737 A CN 201080019737A CN 102449445 A CN102449445 A CN 102449445A
- Authority
- CN
- China
- Prior art keywords
- value
- signal
- length
- mdp
- data
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/003—Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
一种分析具有可以某一旋转速度旋转的部件的机器状态的设备和方法,包括:输入端,用于接收振动信号、所述旋转速度和运算器控制信号。A/D转换器、包络器、降采样器、增强器在时域运行以执行离散自相关,而分析器执行状态监视功能,以便分析机器的状态。
Description
技术领域
本发明涉及机器状态的分析方法,并涉及机器状态的分析设备。本发明也涉及包括这样的设备的系统并涉及操作这样的设备的方法。本发明也涉及使计算机执行分析功能的计算机程序。
背景技术
具有运动部件的机器容易随着时间流逝而受到磨损,这往往使机器的状态退化。这样的具有可运动部件的机器实例为电机、泵、发电机、压缩机、车床和计算机数控机器。可运动部件可以包括轴和轴承。
为了防止机器故障,这样的机器应当取决于机器的状态受到维护。所以优选情况下时常评估这样的机器的运行状态。通过测量从轴承发出的振动或者通过测量机器外壳上的温度,能够确定运行状态,机器的温度取决于轴承的运行状态。对具有旋转或其他运动部件的机器进行这样的状态检查对于安全非常重要,对于这样的机器的寿命长度也非常重要。手工地在机器上进行这样的测量已经公知。普通的做法是由操作员在测量仪器的帮助下在一台或几台机器上的测量点进行测量。
许多商用仪器已经面世,它们取决于以下事实:滚珠轴承中的缺陷产生短脉冲,通常称为冲击脉冲。冲击脉冲测量装置可以产生指示轴承或机器状态的信息。
WO 03062766公开的机器具有测量点以及具有一定轴直径的轴,其中机器使用时轴能够旋转。WO 03062766也公开了具有旋转轴的机器的状态分析装置。所公开的装置具有传感器,用于产生指示测量点振动的测量值。WO 03062766中公开的装置具有数据处理器和存储器。该存储器可以存储程序代码,当程序代码在数据处理器上运行时,将使分析装置执行机器状态监视功能。这样的机器状态监视功能可以包括冲击脉冲测量。
UK专利申请GB 2 190 198 A公开了用于分析具有周期地重复的运行循环的机器的运行的振动分析系统。GB 2 190 198 A教导了提供振动检测装置,在使用中它与该机器操作地相关联,而循环监视器产生与该机器周期率相关的信号。模拟振动信号从振动检测装置经由放大器传递到包络级并从那里到A/D转换器,以某速率对模拟振动信号采样,例如每循环36000次模拟到数字转换。A/D转换后的信号被传递到采样器,它用于将模拟到数字转换器输出表示的十个数值的每个连续组平均从而获得3600个振动值“样点”。这些样点经由时基控制器向样点平均器提供,它具有的寄存器中有多个循环的机器运行循环中对应点的累积的振动值。这许多寄存器的每一个中的总振动值除以运行循环数以给出机器运行循环中每个采样点的线性算术平均振动值。
在对其进行分析操作的预定数目的运行循环结束后,平均振动值的序列将表示对于机器运行循环即在时域中在振动检测装置的位置经历的平均振动。
根据GB 2 190 198的振动分析系统也包括事件识别装置,以确定平均振动值揭示的运行循环中某事件的发生。根据GB 2 190 198,以这种方式或许有可能识别轴承故障和其他事件。
发明内容
本发明的一方面涉及的问题是在具有可以某一旋转速度旋转的部件的机器中能够防止由于机械磨损或故障的未预计到的机器破坏。
通过分析具有可以某一旋转速度旋转的部件的机器状态的设备解决这个问题,包括
输入端,用于接收模拟测量信号,该信号指示具有振动频率和重复频率的振动信号特征;
A/D转换器,用于根据所述模拟测量信号产生数字测量信号,所述数字测量信号具有第一采样频率,所述第一采样频率为所述振动频率的至少两倍;
包络器,用于产生指示所述重复频率的包络信号;
降采样器,用于根据所述包络信号产生降采样后第一数字信号,使得所述降采样后第一数字信号具有降低采样频率;
增强器,具有用于接收所述降采样信号的输入端,其中
所述增强器在时域中运行,对所述降采样后第一数字信号进行离散自相关,以便产生增强器输出信号序列;以及
分析器,用于执行状态监视功能,以便根据所述增强器输出信号序列分析机器的状态。
这种解决方案具有的有利效果是在具有可以某一旋转速度旋转的部件的机器中能够更早检测出初始故障,即更早检测出刚刚开始发展的故障。实现所述早期检测的原因是关于所述增强器输出信号序列中重复信号特征,显著降低了随机信号分量的振幅,从而使所述分析器能够根据所述增强器输出信号序列产生数据指示所述初始故障的诊断。根据这样的诊断可以进行机器保养,并且如果所述诊断数据指示了机器破坏的危急风险,那么所述机器或所述被监视的旋转部件可以停止或减速,以便防止机器破坏。
已经知道在某些机械中发生过未预计到的机器破坏,即使在所述机器已经受到基本上连续的状态监视时。这样的未预计到的破坏对机器拥有者损失很大并可能导致其他的负面效应,如在未预计到的机器故障使机器部件跌落时。
本发明人意识到,在某些机械的机械振动中存在着特别高的噪声级别,并且这样的噪声级别妨碍了对机器故障的检测。因此,对于某些类型的机械,保护性状态监视的常规方法已经无法对即将到来的退化状态提供足够早期和/或可靠的警告。本发明人推断,在这样的机械中可能存在着表明已退化状态的机械振动VMD,但是测量振动的常规方法可能尚未完善。本发明人意识到,常规分析系统的这种尚未完善可能是由于所述测量信号的信噪比SNR低到阻止或妨碍了检测到指示初始故障的任何振动信号分量。
本发明人还意识到,具有缓慢旋转的部件的机器是似乎特别易于突发或未预计到的故障或破坏的机械类型。本发明人还意识到,低旋转速度fROT可以引起机械振动VMD的振幅比旋转速度更高的情况下低。当指示初始机器故障的机械振动VMD具有低振幅时,所述测量信号中的噪声内容将变得相对来说更高。当在具有低于50rpm的旋转速度的机器上测量时,由所述降采样器传递的所述包络和降采样后数字测量信号可能噪声高到在所述降采样后数字测量信号被直接馈入所述分析器时阻止成功的状态监视分析。换言之,所述降采样后数字测量信号的信噪比SNR可能低到阻止或妨碍了检测到任何振动信号分量SD。
本发明人在已经意识到某些机械的机械振动中特别高的噪声级别妨碍对机器故障的检测之后,提出了在噪声环境中能够检测微弱机械信号的方法。如上所述,噪声测量信号SEA中振动信号分量SD的重复频率fD取决于指示被监视机器的旋转部件的初始故障的机械振动VMD。本发明人意识到,有可能检测出初始故障,即刚刚开始发展的故障,只要能够辨别出对应的微弱信号。
因此,测量信号可以包括取决于可旋转运动部件振动的至少一个振动信号分量SD,其中所述振动信号分量具有重复频率fD,它取决于可旋转运动部件的旋转速度fROT。所以,取决于可旋转运动部件的振动的振动信号分量的存在,可以提供被监视机器的退化状态或初始故障的早期指示。
权利要求的解决方案有益地能够检测两个不同的故障特征SD1、SD2,它们可以具有互相不同的重复频率fD1、fD2。因此,所述增强器有益地适于增强具有互相不同的重复频率fD1、fD2的不同特征SD1、SD2。在这个语境中,重复频率值fD1、fD2都与被监视旋转部件的旋转速度fROT成正比,同时重复频率fD1不同于重复频率fD2。
根据所述装置的实施例
所述增强器适于处理所述第一数字信号的第一多个输入样点值,以便产生具有第二多个输出样点值的输出信号序列;所述第二多个为正整数。
可以指示退化状态的相对高频振动信号特征的所述检测和收集要求采样频率,根据尼奎斯特标准,它必须为所述振动信号特征的最高频率的至少两倍。因此,这个收集过程不可避免地导致所述数字测量信号中相对大数目的样点值。因此,尽管在测量样点值的所述检测和收集中需要大数目的样点值,但是期望要在所述增强器中处理的数据量保持尽可能低。
根据权利要求1的解决方案有益地解决了实现以下步骤中的冲突需求
检测和收集具有振动频率fSEA的相对高频振动信号特征,在A/D转换后它产生大数目的样点值;同时也
提供足够小数目的样点值作为所述增强器的输入,以使所述增强器信号处理过程中需要的计算量最小化。
就此而论,提供
包络器和降采样器,以便产生指示所述重复频率的降采样后第一数字信号,使得所述降采样后第一数字信号具有降低的采样频率,向所述增强器传递压缩的重复信号特征,它可以指示退化状态。
根据本发明的实施例,提供了分析具有可以某一旋转速度旋转的部件的机器状态的方法,包括以下步骤:
根据从所述部件的旋转发源的机械振动接收第一数字时域信号;所述第一数字信号具有采样频率;其中,所述第一数字信号可以包括噪声和具有第一重复频率的第一信号特征,进一步包括:
处理所述第一数字信号的第一多个P1输入样点值,以便产生具有第二多个P2输出样点值的输出信号序列;所述输出信号序列为时域信号序列;其中,
所述第一多个P1输入样点值的所述处理包括第三多个P3输入样点值的乘法,用于产生一个输出样点值,使得产生所述第二多个P2输出样点值包括至少第四多个P4乘法。
本发明的一方面涉及的装置分析具有可以某一旋转速度旋转的部件的机器状态,包括
第一输入端,用于根据从所述部件的旋转发源的机械振动接收第一数字信号,所述第一数字信号具有采样频率;其中,所述第一数字信号可以包括噪声和具有第一重复频率的第一信号特征;
增强器,用于处理所述第一数字信号的第一多个输入样点值,以便产生具有第二多个输出样点值的输出信号序列;所述第二多个为正整数;
分析器,用于执行状态监视功能,以便根据所述输出信号序列分析机器的状态;其中,增强器包括:
至少一个第二输入端,用于接收指示所述第二多个的数据;以及
传递所述输出信号序列的输出端;其中
所述增强器适于根据第三多个输入信号乘积产生输出样点值;所述第三多个为正整数;其中,通过将第一样点位置处的第一输入样点值乘以第二样点位置处的第二输入样点值,获得输出样点位置t的输入信号乘积,所述第二输入样点值与所述第一输入样点值分开一定数目的样点位置,其中所述一定数目等于或大于一定值;
所述增强器进一步包括:
至少一个第三输入端,用于接收指示信噪比改进器值的数据;以及
根据以下信息产生指示所述第三多个的数据值的装置
所述信噪比改进器值,以及
指示所述第二多个的所述数据。
这种解决方案能够成功地评估机器的状态,即使在所接收的数字测量信号中噪声非常高时。
根据本发明的实施例,产生指示所述第三多个的数据值的所述装置适于接收
指示所述第一多个的数据,以及
指示所述第二多个的数据,以及
指示所述一定值的数据,以及其中
产生指示所述第三多个的数据值的所述装置适于设置第三多个值CLENGTH等于值ILENGTH与SSTART和OLENGTH之和之间的差。
根据本发明的实施例,产生指示所述第三多个的数据值的所述装置适于接收
指示所述第二多个的数据,以及
指示所述信噪比改进器值的数据,以及其中
产生指示所述第三多个的数据值的所述装置适于设置第三多个值CLENGTH等于第二多个值与所述信噪比改进器值的乘积。
这种解决方案允许用户以有益地容易的和用户友好的方式设置要实现的信噪比改进的级别。
本发明的一方面涉及分析具有可以某一旋转速度旋转的部件的机器状态的方法,包括以下步骤:
根据从所述部件的旋转发源的机械振动接收第一数字时域信号;所述第一数字信号具有采样频率;其中,所述第一数字信号可以包括噪声和具有第一重复频率的第一信号特征,进一步包括:
处理所述第一数字信号的第一多个P1输入样点值,以便产生具有第二P2多个输出样点值的输出信号序列;所述输出信号序列为时域信号序列;其中,
所述第一多个P1输入样点值的所述处理包括第三多个P3输入样点值的乘法,用于产生一个输出样点值,使得产生所述第二多个P2输出样点值包括至少第四多个P4乘法。
根据本发明的实施例,所述方法包括以下步骤:
接收指示信噪比改进器值的数据;以及
根据所述信噪比改进器值和指示所述第二多个的所述数据产生指示所述第三多个的数据值。
附图说明
为了使理解本发明变得简单,将利用实例并参考附图进行说明,其中:
图1显示了根据本发明实施例的状态分析系统2实施例的示意框图;
图2A是图1所示状态分析系统2的一部分的实施例的示意框图;
图2B是传感器接口实施例的示意框图;
图2C是来自振动传感器的测量信号的展示;
图2D展示了冲击脉冲传感器产生的测量信号振幅;
图2E展示了振动传感器产生的测量信号振幅;
图3是根据本发明实施例的冲击脉冲测量传感器的简化展示;
图4是存储器60及其内容的实施例的简化展示;
图5是在具有可运动轴的机器6的客户位置的分析设备实施例的示意框图;
图6展示了根据本发明实施例的预处理器实施例的示意框图;
图7展示了鉴别器230的实施例;
图8展示了鉴别器230的另一个实施例;
图9展示了预处理器200的另一个实施例;
图10A是流程图,展示了增强信号中重复信号模式方法的实施例;
图10B是流程图,展示了产生数字输出信号的方法;
图10C展示了增强器的实施例;
图10D展示了根据增强器方法的实施例的若干信号;
图10E展示了操作增强器用户接口的方法实施例;
图10F展示了操作增强器的方法实施例;
图10G展示了增强器的另一个实施例320;
图10H是表格,展示了输出信号值的一部分计算;
图11是具有多个存储器位置的第一存储器的示意展示;
图12是具有多个存储器位置的第二存储器的示意展示;
图13是包括两个重复信号特征的实例输出信号SMDP的示意展示;
图14A展示了传递给降采样器310的输入的信号中的许多样点值;
图14B展示了对应时间周期的输出样点值;
图15A展示了根据本发明实施例的降采样器;
图15B展示了本发明的另一个实施例;
图16展示的本发明的实施例包括如以上介绍的降采样器和增强器,以及分数降采样器;
图17展示了分数降采样器的实施例;
图18展示了分数降采样器的另一个实施例;
图19展示了降采样器和分数降采样器的另一个实施例;
图20是降采样器和分数降采样器又一个实施例的框图;
图21是流程图,展示了运行图20中降采样器和分数降采样器方法的实施例;
图22A、22B和22C介绍了能够以计算机程序实施的方法;
图23是展示行星齿轮系统的正视图;
图24是图23中行星齿轮系统700的示意侧视图,从图23中箭头SW的方向所见;
图25展示了响应由至少一个传感器10依据行星齿轮系统的旋转所检测的信号,由预处理器200产生并由其输出的示范信号的模拟版本(见图5或图16);
图26展示了图25中所示信号的高振幅区域702A的一部分的实例;
图27展示了包括图26所示小周期扰动903的信号的示范频谱;
图28展示了图25所示的示范信号的一部分的实例;
图29展示了根据本发明实施例的状态分析系统的又一个实施例;
图30是框图,展示了连同用户界面和显示器一起的图29的信号处理组合体的部件;
图31是参数控制器的示意展示。
具体实施方式
在以下说明中不同实施例中的相似特征可以由相同的附图标记指示。
图1显示了根据本发明实施例的状态分析系统2实施例的示意框图。附图标记4涉及具有可运动部件8的机器6的客户机位置。可运动部件可以包括轴承7和轴8,当机器在运行状态时,轴8旋转。当轴旋转时,响应从轴和/或轴承发出的振动,能够确定轴8和轴承7的运行状态。客户机位置4,也可以指客户机部件或用户部件,可以是例如风力发电厂的房屋及地基,即在某位置处的一组风力涡轮,或者某造纸厂的房屋及地基,或者具有可运动部件机器的某个其他制造厂。
当传感器10被连接到机器6的机体上或在其上测量点12处时,状态分析系统2的实施例在运行。尽管图1仅仅展示了两个测量点12,但是应当理解位置4可以包括任何数目的测量点12。图1所示的状态分析系统2包括分析装置14,根据传感器10传递的测量值分析机器的状态。
分析装置14具有通信端口16,用于双向数据交换。通信端口16可连接到通信网络18,如经由数据接口19。通信网络18可以是环球网,也称为因特网。通信网络18还可以包括公共交换电话网。
服务器计算机20被连接到通信网络18。服务器20可以包括数据库22、用户输入/输出接口24和数据处理硬件26,以及通信端口29。服务器计算机20位于位置28,在地理上与客户机位置4分离。服务器位置28可以在第一座城市,比如瑞典的首都斯德哥尔摩,而客户机的位置可以在另一座城市,比如德国的斯图加特或美国密歇根州的底特律。作为替代,服务器位置28可以在某城镇的第一部分,而客户机位置可以在同一城镇的另一部分。服务器位置28还可以称为供应商部件28或供应商部件位置28。
根据本发明的实施例,中央控制位置31包括控制计算机33,具有数据处理硬件和软件,用于测量在客户机位置4的多台机器。机器6可以是风力涡轮或风力涡轮中使用的齿轮箱。作为替代,这些机器可以包括如造纸厂中的机械。控制计算机33可以包括数据库22B、用户输入/输出接口24B和数据处理硬件26B,以及通信端口29B。中央控制位置31可以与客户机位置4分开一段地理距离。利用通信端口29B控制计算机33能够被连接,以经由端口16与分析装置14进行通信。分析装置14可以传递正被部分处理的测量数据,以便允许在中央位置31由控制计算机33进一步执行信号处理和/或分析。
供应商公司占据着供应商部件位置28。供应商公司可以出售和递送分析装置14和/或分析装置14所用的软件。供应商公司还可以出售和递送中央控制位置31处的控制计算机所用的分析软件。以下连同图4讨论了这样的分析软件94、105。这样的分析软件94、105可以通过在所述通信网络18上的传输而递送。
根据系统2的一个实施例,装置14是便携式装置,可以时时连接到通信网络18。
根据系统2的另一个实施例,装置14可以基本上连续地从传感器10接收测量信号,从而允许连续或基本上连续地监视机器状态。根据这个实施例的装置14还可以基本上连续地能够与控制位置31处的控制计算机33通信。因此,根据这个实施例的装置14可以基本上一直“在线”可用于与控制位置31处的控制计算机33进行通信。
根据系统2的一个实施例,装置14基本上被连续地连接到通信网络18。因此,根据这个实施例的装置14可以基本上一直“在线”可用于与供应商计算机20和/或与在控制位置31处的控制计算机33进行通信。
图2A是图1所示的状态分析系统2的一部分的实施例的示意框图。图2A中展示的状态分析系统包括用于产生测量值的传感器单元10。所述测量值可以取决于运动,或者更确切地说,取决于轴旋转时由轴承引起的振动或冲击脉冲。
当设备30被稳固地安装在机器6上或其上的测量点处时,状态分析系统2的实施例在运行。安装在测量点处的设备30可以指双头螺栓30。双头螺栓30可以包括连接耦接32,传感器单元10与其可拆卸地连接。连接耦接32可以例如包括双头螺纹,使得传感器单元能够利用1/4转的旋转与双头螺栓机械地啮合。
测量点12可以包括机器外壳上带螺纹的凹处。双头螺栓30可以具有带有螺纹的突出部分,对应于这些凹处,使得双头螺栓能够像螺栓一样通过纳入到凹处中而牢固地连接到测量点。
作为替代,测量点可以包括机器外壳上带螺纹的凹处,而传感器单元10可以包括对应的螺纹,使得它能够被直接地纳入到所述凹处中。作为替代,仅仅以彩色标记在机器外壳上标注了测量点。
图2A例示的机器6可以具有旋转轴,具有一定的轴直径d1。机器6在使用时,机器中的轴24可以以转速V1旋转。
传感器单元10可以连接到分析机器状态的设备14。参考图2A,分析装置14包括传感器接口40,用于接收由传感器10产生的测量信号即测量数据。传感器接口40被连接到数据处理装置50,它能够按照程序代码控制分析装置14的运行。数据处理装置50也被连接到存储所述程序代码的存储器60。
根据本发明的实施例,传感器接口40包括接收模拟信号的输入端42,输入端42被连接到模-数(A/D)转换器44,其数字输出48被连接到数据处理装置50。A/D转换器44以一定的采样频率fS对所接收的模拟信号进行采样,以便递送具有一定的采样频率fS的数字测量数据信号SMD,并且其中每个样点的振幅都取决于在采样时刻所接收的模拟信号的振幅。
根据展示在图2B中本发明的另一个实施例,传感器接口40包括从冲击脉冲测量传感器接收模拟信号SEA的输入端42、被连接以接收模拟信号的调节电路43以及被连接以从调节电路43接收调节后模拟信号的A/D转换器44。A/D转换器44以一定的采样频率fS对所接收的调节后模拟信号进行采样,以便递送具有一定的采样频率fS的数字测量数据信号SMD,并且其中每个样点的振幅都取决于在采样时刻所接收的模拟信号的振幅。
采样定理保证了有限带宽的信号(即具有最高频率的信号)根据其采样后的形式能够完美地被重建,如果采样频率fS大于被监视的模拟信号SEA的最大频率fSEAmax的两倍的话。等于采样频率一半的频率因此是对能够被采样后的信号SMD明确地表示的最高频率的理论限制。这个频率(采样频率的一半)称为采样系统的尼奎斯特频率。在采样后的信号中能够观察到尼奎斯特频率fN之上的频率,但是它们的频率是不明确的。也就是说,具有频率f的频率分量无法与具有其中B为非零整数的频率B*fN+f和B*fN-f的其他分量区分。这种歧义性称为假频,应对方法可以是先用反假频滤波器(通常是具有接近尼奎斯特频率的截止点的低通滤波器)过滤该信号,再转换到采样后离散表达。考虑到允许非理想滤波器在频率响应中具有一定的斜率,为了提供安全余量,采样频率可以选择高于2的值。因此,根据本发明的实施例,可以将采样频率设置为:
fS=k*fSEAmax
其中,k是具有大于2.0的值的因子。
所以可以选择因子k为大于2.0的值。优选情况下,可以选择因子k为2.0与2.9之间的值,以便提供好的安全余量同时避免产生不必要的许多采样值。根据某实施例,作为有益的情况,选择的因子k使得100*k/2呈现为整数。根据某实施例,因子k可以被设置为2.56。选择k为2.56致使100*k=256=2的8次方。
根据某实施例,数字测量数据信号SMD的采样频率fS可以固定到一定的值fS,比如fS=102.4kHz。
因此,当采样频率fS被固定到一定的值fS时,模拟信号SEA的最高频率fSEAmax将为:
fSEAmax=fS/k
其中fSEAmax是采样后信号中要被分析的最高频率。
因此,当采样频率fS被固定到一定的值fS=102400Hz,而因子k被设置为2.56时,模拟信号SEA的最大频率fSEAmax将为:
fSEAmax=fS/k=102400/2.56=40kHz
所以,响应所接收的模拟测量信号SEA而产生具有一定的采样频率fS的数字测量数据信号SMD。A/D转换器44的数字输出48经由传感器接口40的输出49被连接到数据处理装置50,以便能够向数据处理装置50传递数字测量数据信号SMD。
传感器单元10可以包括振动转换器,传感器单元被构造为在物理上啮合测量点的连接耦接,使得在测量点的机器振动能够被传输到振动转换器。根据本发明的实施例,传感器单元包括具有压电元件的转换器。当测量点12振动时,传感器单元10或者至少其一部分也振动,然后转换器产生电信号,其频率和振幅分别取决于测量点12的机械振动频率和振动振幅。根据本发明的实施例,传感器单元10是振动传感器,提供如频率范围为1.00Hz至10000Hz的10mV/g的模拟振幅信号。这样的振动传感器被设计为基本上是为了传递10mV的相同振幅,而不考虑它是在1Hz、3Hz还是在10Hz被施加1g(9.82m/s2)的加速度。因此典型的振动传感器在高达约10kHz的指定频率范围内具有线性响应。从旋转机器部件中发出的在这个频率范围内的机械振动往往由不平衡或未对准导致。不过,当安装在机器上时,典型情况下,线性响应振动传感器也具有几个不同的机械谐振频率,取决于传感器与振动源之间的物理路径。
滚柱轴承的故障导致相对尖锐的弹性波,称为冲击脉冲,在到达传感器之前沿着机器外壳的物理路径传播。这样的冲击脉冲往往具有很宽的频谱。典型情况下,滚柱轴承冲击脉冲的振幅低于由不平衡或未对准所导致振动的振幅。
宽频谱的冲击脉冲特征使得它们能够激活”鸣震响应”或与传感器相关联的谐振频率的谐振。因此,来自振动传感器的典型测量信号可能具有图2C所示的波形,即在主要低频信号上叠加了更高频率更低振幅谐振的“鸣震响应”。
为了能够分析往往从轴承的故障中发出的冲击脉冲特征,必须滤除低频分量。这能够利用高通滤波器或利用带通滤波器实现。不过,必须调整这些滤波器以便在高频信号部分通过的同时,阻塞低频信号部分。典型情况下,单个的振动传感器将具有与从一个冲击脉冲信号源的物理路径相关联的一个谐振频率,以及与从另一个冲击脉冲信号源的物理路径相关联的不同的谐振频率,正如在US 6,053,047中提及。因此,当使用振动传感器时,旨在使高频信号部分通过的滤波器调整需要个别调整。
当正确地调整了这样的滤波器时,结果信号将由冲击脉冲特征组成。不过,分析从振动传感器发出的冲击脉冲特征稍微受以下事实伤害:振幅响应以及谐振频率固有地取决于冲击脉冲信号的个别物理路径而改变。
作为有益的情况,与振动传感器相关联的这些缺点可以利用冲击脉冲测量传感器减轻。冲击脉冲测量传感器被设计为并适应提供预定的机械谐振频率,以下将进一步详细介绍。
作为有益的情况,冲击脉冲测量传感器的这种特征呈现了可重复的测量结果,因为冲击脉冲测量传感器的输出信号具有稳定的谐振频率,基本上不取决于冲击脉冲信号源与冲击脉冲传感器之间的物理路径。不仅如此,相互不同的各个冲击脉冲传感器提供的谐振频率偏差如果有也非常小。
这个的有益效果是简化了信号处理,因为滤波器不再需要个别调整,与以上介绍的使用振动传感器时的情况大不相同。不仅如此,冲击脉冲传感器的振幅响应定义明确,使得按照由S.P.M.InstrumentAB所定义的适宜测量方法进行测量时,单个测量就提供了可靠的信息。
图2D展示了冲击脉冲传感器产生的测量信号的振幅,而图2E展示了振动传感器产生的测量信号的振幅。两个传感器都被施加了相同系列的机械冲击而没有典型的低频信号内容。正如从图2D和图2E清楚可见,来自冲击脉冲测量传感器的冲击脉冲特征的谐振响应持续时间,要短于对应的来自振动传感器冲击脉冲特征的谐振响应。
冲击脉冲测量传感器提供独特的冲击脉冲特征响应的这种特性,具有的有益效果是从它提供的测量信号有可能区分在短时间跨度内发生的不同机械冲击脉冲之间的差别。
根据本发明的实施例,所述传感器是冲击脉冲测量传感器。图3是根据本发明实施例的冲击脉冲测量传感器10的简化展示。根据这个实施例,传感器包括具有一定质量或重量的部件110和压电元件120。压电元件120有些弹性,所以被施加外力时它能够伸缩。压电元件120在对立表面上分别配备了导电层130和140。当压电元件120伸缩时它产生电信号,由导电层130和140拾取。所以,机械振动被转换为模拟测量电信号SEA,它被传递到输出端子145、150。压电元件120可以位于配重110与表面160之间,在运行期间,它被物理地连接到测量点12,如图3展示。
冲击脉冲测量传感器10具有的谐振频率取决于传感器机械特征,比如配重部件110的质量m和压电元件120的弹性。因此,压电元件具有弹性和弹簧常数k。传感器的机械谐振频率fRM所以也取决于质量m和弹簧常数k。
根据本发明的实施例,传感器的机械谐振频率fRM能够由下面的等式确定:
根据另一个实施例,冲击脉冲测量传感器10的实际机械谐振频率也可以取决于其他因素,比如传感器10与机器6机体连接的性质。
因而谐振冲击脉冲测量传感器10对于以机械谐振频率fRM或其附近的频率的振动特别灵敏。冲击脉冲测量传感器10可以被设计为使得机械谐振频率fRM在28kHz到37kHz的范围内某处。根据另一个实施例,机械谐振频率fRM在30kHz到35kHz的范围内某处。
所以模拟测量电信号具有的电振幅可以在频谱上变化。为了说明理论背景的目的,可以假设要是对冲击脉冲测量传感器10施加从如1Hz到如200000kHz的所有频率振幅一致的机械振动,那么来自冲击脉冲测量传感器的模拟信号SEA的振幅在机械谐振频率fRM处将具有最大值,因为传感器以该频率“推动”时将会谐振。
参考图2B,调节电路43接收模拟信号SEA。调节电路43可以被设计为阻抗适应电路,被设计为适应A/D转换器的从传感器端子145、150所见的输入阻抗,以便出现最佳的信号传输。因此,调节电路43可以运行以适应从传感器端子145、150所见的输入阻抗Zin,以便将最大的电功率传递给A/D转换器44。根据调节电路43的实施例,模拟信号SEA被供给了变压器的初级线圈,而调节后模拟信号由变压器的次级线圈传递。初级线圈有n1匝而次级线圈有n2匝,比值n1/n2=n12。因此,A/D转换器44被连接以从调节电路43接收调节后的模拟信号。A/D转换器44具有输入阻抗Z44,而在调节电路43被连接在传感器端子145、150与A/D转换器44的输入端子之间时,从传感器端子145、150所见的A/D转换器的输入阻抗将是(n1/n2)2*Z44。
A/D转换器44以一定的采样频率fS对所接收的调节后模拟信号进行采样,以便传递具有一定的采样频率fS的数字测量数据信号SMD,并且其中每个样点的振幅都取决于在采样时刻所接收的模拟信号的振幅。
根据本发明的实施例,数字测量数据信号SMD被传递到装置180进行数字信号处理(见图5)。
根据本发明的实施例,用于数字信号处理的装置180包括数据处理器50和使数据处理器50执行数字信号处理的程序代码。根据本发明的实施例,处理器50由数字信号处理器实施。数字信号处理器也被称为DSP。
参考图2A,数据处理装置50被连接到存储所述程序代码的存储器60。优选情况下,程序存储器60是非易失性存储器。存储器60可以是读/写存储器,即从存储器中读取数据和把新数据写到存储器60上都能够进行。根据某实施例,程序存储器60由FLASH存储器实施。程序存储器60可以包括第一存储器段70,用于存储可执行的第一组程序代码80,以便控制分析装置14执行基本的操作(图2A和图4)。程序存储器还可以包括第二存储器段90,用于存储第二组程序代码94。在第二存储器段90中的第二组程序代码94可以包括使分析装置处理所检测的信号或若干信号,以便产生预处理的信号或一组预处理的信号的程序代码。存储器60还可以包括第三存储器段100,用于存储第三组程序代码104。在第三存储器段100中的程序代码组104可以包括使分析装置执行所选定分析功能105的程序代码。当执行分析功能时,它可以使分析装置在用户界面106上呈现对应的分析结果,或者在端口16上传递该分析结果(见图1和图2A以及图7和图8)。
数据处理装置50还连接到进行数据存储的读/写存储器52。不仅如此,数据处理装置50还可以连接到分析装置通信接口54。分析装置通信接口54提供与测量点通信接口56的双向通信,接口56可连接在机器的测量点上,或者在其附近。
测量点12可以包括连接耦接32、可读和可写信息载体58以及测量点通信接口56。
可写信息载体58和测量点通信接口56可以在双头螺栓30附近放置的单独设备59中提供,正如图2展示。作为替代,可写信息载体58和测量点通信接口56可以在双头螺栓30内提供。这在WO 98/01831中进行了更详细的介绍,其内容在此引用作为参考。
系统2被安排为允许测量点通信接口56与分析装置通信接口54之间的双向通信。优选情况下,测量点通信接口56和分析装置通信接口54被构造为允许无线通信。根据某实施例,测量点通信接口和分析装置通信接口被构造为彼此通过射频(RF)信号进行通信。这个实施例包括在测量点通信接口56中的天线,以及分析装置通信接口54中的另一条天线。
图4是存储器60及其内容的实施例的简化展示。该简化展示意在传达对存储器60中存储不同程序函数的总体思路的理解,而这不一定是对在真实存储器电路中存储程序的方式的正确的技术教导。第一存储器段70存储的程序代码用于控制分析装置14以执行基本的操作。尽管图4的简化展示显示了伪代码,但是应当理解,程序代码80可以由机器代码构建,或者由数据处理装置50能够执行或解释的任何级别的程序代码构建(图2A)。
图4展示的第二存储器段90存储着第二组程序代码94。在段90中的程序代码94,当运行在数据处理装置50上时,将使分析装置14执行某功能,比如数字信号处理功能。所述功能可以包括对数字测量数据信号SMD的高级数学处理。根据本发明的实施例,程序代码94适于使处理器装置50执行连同在本文档中图5、图6、图9和/或图16所介绍的信号处理功能。
正如以上连同图1所提及,控制分析装置功能的计算机程序可以从服务器计算机20中下载。这意味着要被下载的程序在通信网络18上被传输。通过在通信网络18上调制载波以携带该程序就能够实现这一点。所以下载后的程序可以被加载到数字存储器中,比如存储器60(见图2A和图4)。因此,信号处理程序94和/或分析功能程序104、105可以经由通信端口接收,比如端口16(图1和图2A),以便将其加载到存储器60中。同样,信号处理程序94和/或分析功能程序104、105也可以经由通信端口29B(图1)接收,以便将其加载到计算机26B的程序存储器位置或数据库22B中。
本发明的一方面涉及计算机程序产品,比如可加载到某装置的数字存储器中的程序代码装置94和/或程序代码装置104、105。计算机程序产品包括软件代码部分,当所述产品运行在分析机器状态的设备的数据处理单元50上时,用于执行信号处理方法和/或分析功能。术语“运行在数据处理单元上”意味着计算机程序加上数据处理单元执行本文档中介绍种类的方法。
措词“计算机程序产品,可加载到状态分析装置的数字存储器中”意味着计算机程序能够被引入到状态分析装置的数字存储器中,以便实现编程为能够或适于执行以上介绍种类的方法的状态分析装置。术语“加载到状态分析装置的数字存储器中”意味着以这种方式编程的状态分析装置能够或适于执行以上介绍种类的方法。
上述计算机程序产品也可以加载到计算机可读介质上,比如光盘或DVD。这样的计算机可读介质可以用于向客户机传递所述程序。
根据分析装置14(图2A)的实施例,它包括用户输入接口102,操作员由此可以与分析装置14互动。根据某实施例,用户输入接口102包括一组按钮104。分析装置14的实施例包括用户输出界面106。用户输出界面可以包括显示器单元106。当数据处理装置50运行基本程序代码80中提供的基本程序功能时,它利用用户输入接口102和显示器单元106提供用户交互。按钮组104可限于几个按钮,比如五个按钮,如图2A展示。中央按钮107可以用作ENTER或SELECT功能,而可以使用其他更多的外围按钮可以用于移动显示器106上的光标。应当理解,以这种方式,经由用户接口可以把符号和文本输入到装置14中。例如,显示器单元106可以显示许多符号,比如字母表中的字母,同时光标在显示器上可移动以响应用户输入,所以允许用户输入信息。
图5是分析装置14的实施例的示意框图,处在具有可运动轴8的机器6所在的客户机位置4。传感器10可以是冲击脉冲测量传感器,被显示为连接到机器6的机体,以便拾取机械振动和向传感器接口40传递表明所检测到的机械振动的模拟测量信号SEA。传感器接口40可以被设计为如连同图2A或图2B的介绍。传感器接口40向用于数字信号处理的装置180传递数字测量数据信号SMD。
数字测量数据信号SMD具有采样频率fS,而每个样点的振幅值都取决于采样时刻所接收的模拟测量信号SEA的振幅。根据某实施例,数字测量数据信号SMD的采样频率fS可以被固定到一定的值fS,比如fS=102kHz。采样频率fS可以由时钟190所传递的时钟信号控制,如图5展示。时钟信号也可以传递给用于数字信号处理的装置180。用于数字信号处理的装置180能够响应所接收的数字测量数据信号SMD、时钟信号以及采样频率fS与时钟信号之间关系,产生有关所接收的数字测量数据信号SMD的短暂持续时间的信息,因为两个连续样点值之间的持续时间等于TS=1/fS。
根据本发明的实施例,用于数字信号处理的装置180包括预处理器200,用于执行数字测量数据信号SMD的预处理,以便将预处理后的数字信号SMDP传递到输出210上。输出210被连接到鉴别器230的输入220。鉴别器230适于鉴别预处理后的数字信号SMDP,以便将鉴别结果传递给用户界面106。作为替代,鉴别结果也可以传递给通信端口16,以便使结果能够传输到比如在控制场所的控制计算机33(见图1)。
根据本发明的实施例,连同用于数字信号处理的装置180中功能框所介绍的功能,预处理器200和鉴别器230可以由计算机程序代码94和/或104实施,正如以上连同图4结合存储器块90和100的介绍。
用户可以仅仅需要几项基本的监视功能,检测某机器的状态是正常还是异常。检测出异常状态后,用户可以要求专门的职业维护人员确定该问题的精确性质,并且执行必要的维修工作。职业维护人员经常需要和使用宽范围的鉴别功能,使得有可能确定异常机器状态的性质和/或原因。因此,分析装置14的不同用户对装置的功能可以提出非常不同的要求。在本文档中术语状态监视功能是用于检测某机器的状态是正常还是有点退化还是异常的功能。术语状态监视功能还包括鉴别功能,使得有可能确定异常机器状态的性质和/或原因。
机器状态监视功能的实例
状态监视功能F1、F2...Fn包括若干功能,比如:振动分析、温度分析、冲击脉冲测量、冲击脉冲测量数据的频谱分析、振动测量数据的快速傅氏变换、状态数据在用户界面上的图形表示、在所述机器上的可写信息载体中存储状态数据、在所述装置中的可写信息载体中存储状态数据、转速测量、不平衡检测以及未对准检测。
根据某实施例,装置14包括以下功能:
F1=振动分析;
F2=温度分析,
F3=冲击脉冲测量,
F4=冲击脉冲测量数据的频谱分析,
F5=振动测量数据的快速傅氏变换,
F6=状态数据在用户界面上的图形表示,
F7=在所述机器上的可写信息载体中存储状态数据,
F8=在所述装置中的可写信息载体52中存储状态数据,
F9=转速测量,
F10=不平衡检测,以及
F11=未对准检测。
F12=从所述机器上的可写信息载体58中检索状态数据。
F13=执行振动分析功能F1和执行功能F12“从所述机器上的可写信息载体58中检索状态数据”,以便能够根据当前振动测量数据和历史振动测量数据进行对比或趋势判断。
F14=执行温度分析F2;并且执行功能“从所述机器上的可写信息载体58中检索状态数据”,以便能够根据当前温度测量数据和历史温度测量数据进行对比或趋势判断。
F15=从所述机器上的可写信息载体58中检索标识数据。
功能F7“在所述机器上的可写信息载体中存储状态数据”,以及F13振动分析和状态数据检索的实施例在WO 98/01831中进行了更详细地介绍,其内容在此引用作为参考。
图6展示了根据本发明实施例的预处理器200实施例的示意框图。在这个实施例中,数字测量数据信号SMD被连接到数字带通滤波器240,它具有较低截止频率fLC、较高截止频率fUC以及较高和较低截止频率之间的通频带带宽。
数字带通滤波器240的输出被连接到数字包络器250。根据本发明的实施例,包络器250的信号输出被传递到输出260。预处理器200的输出260被连接到数字信号处理组合体180的输出210以向鉴别器230的输入220传递。
可以选择数字带通滤波器240的较高和较低截止频率,使得在传感器谐振频率fRM处的信号SMD的频率分量在通频带带宽中。如上所述,由在谐振频率fRM处机械谐振的传感器实现机械振动的放大。所以,模拟测量信号SEA反映了在谐振频率fRM处及其周围振动的放大值。因此,作为有益的情况,根据图6实施例的带通滤波器抑制了在谐振频率fRM之下和之上频率的信号,以便进一步增强测量信号在谐振频率fRM处的分量。不仅如此,作为有益的情况,数字带通滤波器240进一步降低了测量信号中固有地包括的噪声,因为在较低截止频率fLC以下和较高截止频率fUC以上的任何噪声分量也被除去或减弱。因此,使用具有在最低谐振频率值fRML到最高谐振频率值fRMU范围内具有机械谐振频率fRM的谐振冲击脉冲测量传感器10时,数字带通滤波器240可以被设计为使较低截止频率fLC=fRML,和较高截止频率fUC=fRMU。根据某实施例,较低截止频率fLC=fRML=28kHz,和较高截止频率fUC=fRMU=37kHz。
根据另一个实施例,机械谐振频率fRM在从30kHz到35kHz范围内某处,然后数字带通滤波器240可以被设计为具有较低截止频率fLC=30kHz,和较高截止频率fUC=35kHz。
根据另一个实施例,数字带通滤波器240可以被设计为使较低截止频率fLC低于最低的谐振频率值fRML,和较高截止频率fUC高于最高的谐振频率fRMU。例如,机械谐振频率fRM可以是从30kHz到35kHz范围内的频率,然后数字带通滤波器240可以被设计为使较低截止频率fLC=17kHz,和较高截止频率fUC=36kHz。
所以,数字带通滤波器240传递的通频带数字测量数据信号SF具有有益的低噪声内容并反映了通频带内的机械振动。通频带数字测量数据信号SF被传递到包络器250。
所以,数字包络器250收到了通频带数字测量数据信号SF,它可以反映具有正的以及负的振幅的信号。参考图6,所接收的信号被数字整流器270整流,并且整流后的信号可以被可选的低通滤波器280滤波以产生数字包络信号SENV。
所以,信号SENV是响应滤波后的测量数据信号SF而产生的包络信号的数字表示。根据本发明的某些实施例,可以除去可选的低通滤波器280。以下连同图9讨论了一个这样的实施例。所以,当以下连同图9讨论的降采样器310包括低通滤波器功能时,可以除去在包络器250中的可选低通滤波器280。
根据本发明的图6实施例,信号SENV被传递到预处理器200的输出260。因此,根据本发明的实施例,在输出210(图5)上传递的预处理后数字信号SMDP是数字包络信号SENV。
尽管用于响应测量信号而产生包络信号的现有技术模拟设备采用了模拟整流器,该模拟整流器会固有地导致系统误差被引入结果信号中,但是作为有益的情况,数字包络器250将产生真正的整流而没有任何系统误差。所以,数字包络信号SENV将具有良好的信噪比,因为在数字带通滤波器240通频带中的谐振频率处进行机械谐振的传感器会引起高的信号振幅,并且在数字域中进行的信号处理不增加噪声,并且不增加系统误差。
参考图5,预处理后的数字信号SMDP被传递到鉴别器230的输入220。
根据另一个实施例,滤波器240是具有截止频率fLC的高通滤波器。这个实施例通过以高通滤波器240代替带通滤波器而简化了设计,从而将低通滤波处理留给了下游另一个低通滤波器,比如低通滤波器280。高通滤波器240的截止频率fLC被选择为接近谐振冲击脉冲测量传感器10的最低预期机械谐振频率值fRML。当机械谐振频率fRM在30kHz到35kHz范围内某处时,高通滤波器240可以被设计为使较低截止频率fLC=30kHz。然后把高通滤波后的信号传给整流器270以及继续到低通滤波器280。根据某实施例,应当有可能使用谐振频率在20kHz到35kHz范围内某处的传感器10。为了实现这个目的,高通滤波器240可以被设计为使较低截止频率fLC=20kHz。
图7展示了鉴别器230(也见图5)的实施例。鉴别器230的图7实施例包括状态分析器290,适于接收指示机器6状态的预处理后的数字信号SMDP。利用在控制输入端300上传递的选择信号,可以控制状态分析器290执行所选定的状态分析功能。利用用户与用户接口102(见图2A)的互动可以产生在控制输入端300上传递的选择信号。当所选定的分析功能包括快速傅氏变换时,选择信号300将设置分析器290以在频率域中对输入信号进行操作。
取决于要执行的分析是什么类型,状态分析器290可以在时域对输入的预处理后数字信号SMDP进行操作,或者在频率域对输入的预处理后数字信号SMDP进行操作。所以,取决于在控制输入端300上传递的选择信号,如图8所示可以包括FFT 294,或者如图7所展示可以把信号SMDP直接传递给分析器290。
图8展示了鉴别器230的另一个实施例。在图8的实施例中,鉴别器230包括可选的快速傅氏变换器294,被连接以从鉴别器230输入220接收信号。来自快速傅氏变换器294的输出可以被传递到分析器290。
为了分析旋转部件的状态,所期望的是在足够长的时间监视所检测的振动,以便能够检测出重复的信号。一定的重复信号特征指示了旋转部件的已退化状态。分析重复信号特征还可以指示已退化状态的类型。这样的分析也可以引起对已退化状态程度的检测。
因此,测量信号可以包括取决于可旋转运动部件8振动的至少一种振动信号分量SD;其中所述振动信号分量具有取决于可旋转运动部件8的转速fROT的重复频率fD。所以,取决于可旋转运动部件8振动的振动信号分量可以指示所监视机器的已退化状态或故障。事实上,振动信号分量SD的重复频率fD与可旋转运动部件8的转速fROT之间的关系可以指示哪个机械部件有故障。因此,在具有多个旋转部件的机器中,有可能通过使用分析功能105(包括频率分析),利用处理测量信号而识别出单个轻微损坏的部件。
这样的频率分析可以包括对包括振动信号分量SD的测量信号的快速傅氏变换。快速傅氏变换(FFT)使用了一定的频率分辨率。该一定的频率分辨率可以按照频率片段表示,决定了辨别不同频率的界限。术语“频率片段”有时指“线条”。如果需要频率分辨率提供高达轴速的Z个频率片段,那么有必要在该轴旋转X周期间记录该信号。
连同旋转部件的分析,分析高于旋转部件的旋转频率fROT的信号频率可能引人关注。旋转部件可以包括轴和轴承。轴旋转频率fROT往往被称为“级别1”。所关注的轴承信号可能发生每轴转约10次(级别10),即故障重复频率fD(以Hz测量)除以转速fROT(以rps测量)等于10Hz/rps,即级别y=fD/fROT=10Hz/rps。不仅如此,分析轴承信号的谐波可能引人关注,所以测量高达级别100的情况可能引人关注。令最大级别为Y,而在FFT在要使用的频率片段总数为Z,以下情况成立:Z=X*Y。相反X=Z/Y,其中,
X是分析数字信号期间所监视轴的转数;以及
Y是最大级别;以及
Z是以许多频率片段表示的频率分辨率。
考虑把降采样后数字测量信号SMDP(见图5)传递到FFT分析器294时的情况,如图8中介绍:在这样的情况下,当FFT分析器294被设置为Z=1600频率片段,而用户关注的分析频率高达级别Y=100,那么X的值变为X=Z/Y=1600/100=16。
因此,当期望Z=1600频率片段并且用户关注的分析频率高达级别Y=100时,有必要在X=16轴转期间进行测量。
使用用户界面102、106(图2A)可以设置FFT分析器294的频率分辨率Z。
因此,使用用户界面102、106(图2A)可以设置用于状态分析功能105和/或信号处理功能94(图4)的频率分辨率值Z。
根据本发明的实施例,通过从一组值中选择一个值Z可以设置频率分辨率Z。频率分辨率Z的可选择值组可以包括:
Z=400
Z=800
Z=1600
Z=3200
Z=6400
如上所述,采样频率fS可以被固定到一定的值,比如fS=102400kHz,而因子k可以设置为2.56,从而致使要被分析的最高频率fSEAmax为:
fSEAmax=fS/k=102400/2.56=40kHz
对于具有转速fROT=1715rpm=28.58rps轴的机器,选定的级别值Y=100致使要分析的最大频率为:
fROT*Y=28.58rps*100=2858Hz。
快速傅氏变换器294可以适于对所接收的具有一定数目样点值的输入信号,执行快速傅氏变换。当将样点值的该一定数目设置为可以被二(2)除而不会出现分数的偶整数时是有益的。
所以,表示从轴的旋转而发出机械振动的数据信号可以包括重复的信号模式。因此在被监视轴的每转,一定的信号模式可以被重复一定量的次数。不仅如此,重复信号也可以以互相不同的重复频率发生。
在Victor Wowk所著“Machinery Vibration Measurements andAnalysis”(ISBN 0-07-071936-5)的书中,在149页上提供了互相不同重复频率的几个实例:
“轴承保持器损坏频率(FTF)
轴承滚动件损坏(BS)频率
外环(OR)
内环(IR)”
该书还在150页上提供了计算这些特定频率的公式。VictorWowk所著“Machinery Vibration Measurements and Analysis”的书中内容在此引用作为参考。确切地说,计算这些特定频率的上述公式在此引用作为参考。同一书中151页上的表指示,这些频率也随轴承制造商而改变,并且:
FTF可以具有0.378的轴承频率因子;
BS可以具有1.928的轴承频率因子;
OR可以具有3.024的轴承频率因子;以及
IR可以具有4.976的轴承频率因子。
频率因子与轴的转速相乘以获得重复频率。该书指出,对于具有转速为1715rpm,即28.58Hz的轴,从标准类型6311的轴承外环(OR)发出脉冲的重复频率可能为大约86Hz;而FTF重复频率可能为10.8Hz。
当被监视轴以转速旋转时,可以按照被监视轴每时间单位的重复或按照每转的重复讨论这样的重复频率,而不需要在两者之间进行区分。不过,如果机器部件以转速旋转,事情就进一步复杂化,如以下连同图16、图17和图20的讨论。
出现突然故障的机械
某些类型的机械可能非常突然地遭受整体故障或破坏。对于某些机器类型,比如风力电站中的旋转部件,已经知道发生过突然的破坏,并且使维护人员和机器拥有者感到十分惊奇。这样的突然破坏对机器拥有者损失很大并可能导致其他的负面效应,如在未预计到的机器故障使机器部件跌落时。
本发明人意识到,在某些机械的机械振动中存在着特别高的噪声级别,并且这样的噪声级别妨碍了对机器故障的检测。因此,对于某些类型的机械,保护性状态监视的常规方法已经无法对即将到来的退化状态提供足够早期和/或可靠的警告。本发明人推断,在这样的机械中可能存在着表明已退化状态的机械振动VMD,但是测量振动的常规方法可能尚未完善。
本发明人还意识到,具有缓慢旋转的部件的机器是似乎特别易于突发故障的机械类型。本发明人还意识到低转速fROT可以导致机械振动VMD的较低振幅。当机械振动VMD表明初始机械故障具有低振幅,测量信号中的噪声内容相对而言变得更高。当在旋转速度小于50rpm的机器上测量时,由降采样器310传送的包络和降采样数字测量信号SRED可以是多噪声的以致如果直接将降采样数字测量信号SRED送入分析器290将阻碍成功的状态监视分析。换句话说,降采样数字测量信号SRED的信噪比SNR可以是如此地低以致阻碍检测任何振动信号分量SD。
本发明人在已经意识到某些机械的机械振动中特别高的噪声级别妨碍对机器故障的检测之后,提出了在噪声环境中能够检测微弱机械信号的方法。如上所述,测量信号SEA中振动信号分量SD的重复频率fD取决于指示被监视机器6的旋转部件8的初始故障的机械振动VMD。本发明人意识到,有可能检测出初始故障,即刚刚开始发展的故障,只要能够辨别出对应的微弱信号。
因此,测量信号可以包括取决于可旋转运动部件8振动的至少一个振动信号分量SD,其中所述振动信号分量具有重复频率fD,它取决于可旋转运动部件8的转速fROT。所以,取决于可旋转运动部件8的振动的振动信号分量的存在,可以提供被监视机器的退化状态或初始故障的早期指示。
在风力涡轮的应用中,其轴承被分析的转轴可以以小于每分120转的速度旋转,即轴的旋转频率fROT小于每秒2转(rps)。有时要被分析的这种轴以小于每分50转的速度旋转(rpm),即小于0.83rps的轴旋转频率fROT。事实上,典型情况下,转速可能小于15rpm。尽管具有1715rpm转速的轴,如在上述书中讨论,在17.5秒中产生了500转;而以每分50转旋转的轴却要用十分钟才产生500转。某些大型风力电站的轴,典型情况下可以以12RPM=0.2rps旋转。
所以,当要被分析的轴承与缓慢旋转轴相关联,并且监视轴承的检测器产生的模拟测量信号SEA使用大约100kHz的采样频率采样时,与该轴一次完整旋转相关联的样点值的数目变得非常大。作为展示性实例,当该轴以50rpm旋转时,为了描述500转以100kHz的采样频率获得了6千万(60000000)个样点值。
不仅如此,当信号包括如此多的样点时,对信号执行高级数学分析需要大量时间。所以,期望先减少每秒的样点数目,再进一步处理信号SENV。
图9展示了预处理器200的另一个实施例。预处理器200的图9实施例包括数字带通滤波器240和数字包络器250,如以上连同图6的介绍。如上所述,信号SENV是响应滤波后测量数据信号SF产生的包络信号的数字表示。
根据预处理器200的图9实施例,数字包络信号SRED被传递到降采样器310,它适于产生具有降低的采样频率SSR1的数字信号SRED。降采样器310运行以产生输出数字信号,其中两个连续样点值之间的短暂持续时间要长于输入信号中两个连续样点值之间的短暂持续时间。以下连同图14更详细地介绍降采样器。根据本发明的实施例,可以除去可选的低通滤波器280,如上所述。在图9实施例中,当由数字整流器270产生的信号被传递到包括低通滤波的降采样器310时,就可以除去低通滤波器280。
降采样器310的输出312把数字信号SRED传递到增强器320的输入端315。增强器320能够接收数字信号SRED并响应它而产生输出信号SMDP。输出信号SMDP被传递到预处理器200的输出端口260。
图10A是流程图,展示了增强信号中重复信号模式方法的实施例。有益的是,这种方法可以用于在表示具有旋转轴的机器状态的信号中,增强重复信号模式。增强器320可以被设计为根据图10A展示的方法运行。
图10A中方法步骤S1000到S1040表示在实际产生输出信号值之前为了进行设置而采取的准备动作。一旦已经执行了这些准备动作,就可以计算输出信号值,如参考步骤S1050的介绍。
图10B是流程图,展示了产生数字输出信号的方法。更确切地说,图10B展示了已经执行了参考图10A中步骤S1000至S1040所介绍的准备动作时,产生数字输出信号的方法实施例。
参考图10A的步骤S1000,确定了输出信号SMDP的期望长度OLENGTH。
图11是具有多个存储器位置i的第一存储器的示意展示。第一存储器的存储器位置i保存着包括数字值序列的实例输入信号。根据本发明的实施例,实例输入信号用于计算输出信号SMDP。图11显示了输入信号I的许多连续数字值的某些。输入信号I中的数字值2080仅仅展示了在输入信号中出现的几个数字值。在图11中,输入信号中的两个相邻数字值被持续时间tdelta分开。值tdelta是增强器320所接收的输入信号的采样频率fSR的倒数(见图9和图16)。
图12是具有多个存储器位置t的第二存储器的示意展示。第二存储器的存储器位置t保存着包括数字值序列的实例输出信号SMDP。因此,图12展示了将数字值3090存储在连续存储器位置中的一部分存储器。图12显示了输出信号SMDP的连续数字值。输出信号SMDP中的数字值3090仅仅展示了在输出信号中出现的几个数字值。在图12中,输出信号中的两个相邻数字值可以被持续时间tdelta在时间上分开。
参考图10的步骤S1000,可以选择输出信号SMDP的期望长度OLENGTH3010,以便有可能使用输出信号SMDP分析输出信号中的一定的频率。例如,如果较低的频率比较高的频率更受关注,就需要更长的输出信号。使用输出信号能够分析的最低频率是1/(OLENGTH*tdelta),其中OLENGTH是输出信号中样点值的数目。如果fSR是输入信号I的采样频率,那么每个数字样点值之间的时间tdelta将是1/fSR。如上所述,重复信号模式可能出现在表示机械振动的数据信号中。所以,测量信号,比如由包络器250传递的信号SENV和传递给增强器320的信号SRED可以包括取决于可旋转运动部件8的振动的至少一个振动信号分量SD,其中所述振动信号分量SD具有重复频率fD,它取决于可旋转运动部件8的转速fROT。因此,在输出信号SMDP中的连续数字值被持续时间tdelta分开时,为了可靠地检测出重复频率为fREP=fD=1/(OLENGTH*tdelta)的重复信号模式的出现,输出信号SMDP必须包括至少OLENGTH个数字值。
根据某实施例,用户可以输入表示要被检测的最低重复频率fREPmin的值,以及要被监视的轴的最低预期转速的有关信息。分析系统2(图1)包括用于响应这些值计算变量OLENGTH的适合值的功能。
作为替代,参考图2A,分析装置14的用户利用经由用户界面102输入对应值,可以设置输出信号SMDP的值OLENGTH3010。
在下一个步骤S1010中,选择了长度因子L。长度因子L确定了输出信号SMDP中随机信号被抑制的程度。较高的L值比较低的L值在输出信号SMDP中给出更少的随机信号。因此,长度因子L可以被称为信噪比改进值。根据本方法的一个实施例,L是1至10之间的整数,但是也可以将L设置为其他值。根据本方法的实施例,L值可以在增强器320中预置。根据本方法的另一个实施例,L值由本方法的用户经由用户接口102(图2A)输入。因子L的值对计算输出信号所需要的计算时间也有影响。较高的L值比较低的L值需要更长的计算时间。
下一步,在步骤S1020,设置了起始位置SSTART。起始位置SSTART是输入信号I中的位置。
设置起始位置SSTART是为了避免或减少输出信号SMDP中非重复模式的出现。如果设置起始位置SSTART以便使得起始位置之前的输入信号部分2070有一定的时间间隔TSTOCHASTIC_MAX对应的长度,那么具有对应频率fSTOCHASTIC_MAX和更高频率的随机信号将在输出信号O、SMDP中被衰减。
在下一个步骤S1030中,计算输入数据信号的所要求长度。输入数据信号的所要求长度在步骤S1030中根据以下公式(1)计算:
(1)ILENGTH=OLENGTH*L+SSTART+OLENGTH
下一步,在步骤S1040中,计算输入数据信号中的长度CLENGTH。长度CLENGTH是在其上执行输出数据信号计算的长度。根据以下的公式(3)计算这个长度CLENGTH:
(3)CLENGTH=ILENGTH-SSTART-OLENGTH
也能够将公式(3)写为:ILENGTH=CLENGTH+SSTART+OLENGTH
然后在步骤S1050中计算输出信号。根据以下公式(5)计算输出信号。在公式(5)中,对输出信号中时间值t计算了输出信号的值。
输出信号SMDP具有长度OLENGTH,如上所述。为了获得整个输出信号SMDP,必须用公式(5)计算从t=1到t=OLENGTH每个时间值的值。在图11中,数字值2081展示了在输出信号计算中使用的一个数字值。数字值2081展示了在输出信号计算(其中i=1)中使用的一个数字值。数字值2082展示了在输出信号计算中使用的另一个数字值。附图标记2082指以上公式(5)中的数字值I(1+SSTART+t),此时i=1和t=1。因此,附图标记2082展示了输入信号中在位置号P处的数字样点值:
P=1+SSTART+1=SSTART+2
在图12中,附图标记3091指在输出信号(其中t=1)中的数字样点值SMDP(t)。
现在将介绍本方法运行增强器320的另一个实施例,用于增强表示具有旋转轴的机器状态的信号中重复模式。根据某实施例,长度OLENGTH可以在增强器320中预置。根据本方法的其他实施例,通过用户经由用户接口102(图2A)输入可以设置长度OLENGTH。根据本方法的优选实施例,变量OLENGTH被设置为可以被二(2)除而不会出现分数的偶整数。有益的是,根据这条规则选择变量调整输出信号中的样点数,使其适于在可选的快速傅氏变换器294中使用。因此,根据本方法的实施例,优选情况下,变量OLENGTH的值被设置为比如1024、2048、4096的数字。
在特别有益的实施例中,在步骤S1020中设置SSTART值,以便使得起始位置之前的输入信号部分2070与输出信号3040的长度相同,即SSTART=OLENGTH。
正如以上连同公式(1)的例示,输入数据信号的所要求长度是:
ILENGTH=OLENGTH*L+SSTART+OLENGTH
因此,在公式(1)中设置SSTART+OLENGTH致使
ILENGTH=OLENGTH*L+OLENGTH+OLENGTH=OLENGTH*L+OLENGTH*2
所以,输入信号的所要求的长度能够根据以下公式(6)按照输出信号的长度表示:
(6)ILENGTH=(L+2)*OLENGTH
其中L是以上讨论的长度因子,而OLENGTH是输出信号中数字值的数目,如以上讨论。
在本发明的这个实施例中,根据以下公式(7)能够计算长度OLENGTH。
(7)CLENGTH=L*OLENGTH
当已经执行了参考图10A中步骤S1000至S1040所介绍的准备动作时,利用参考图10B所介绍的方法就可以产生数字输出信号。根据本发明的实施例,参考图10B介绍的方法利用DSP 50(图2A)执行。
在步骤S1100(图10B)中,增强器320接收数字输入信号I,在输入端315上具有第一ILENGTH个样点值(见图9和/或图16)。正如以上指出,数字输入信号I可以表示从轴旋转发源的机械振动,只要导致出现了具有重复周期为TR的振动。
接收的信号值被存储(步骤S1120)在与增强器320相关联的数据存储器的输入信号存储器部分中。根据本发明的实施例,该数据存储器可以由读/写存储器52实施(图2A)。
在步骤S1130中,在以上等式(5)中使用的变量t被设置为初始值。初始值可以为1(一)。
在步骤S1140中,计算了样点号为t的输出样点值SMDP(t)。计算可以采用下面的等式:
结果样点值SMDP(t)被存储(图10B步骤S1150)在存储器52的输出信号存储器部分中(见图12)。
在步骤S1160中,本过程检查变量t的值,并且如果t的值表示的数目低于输出样点值的所期望数目OLENGTH,则先执行增加变量t的值的步骤S1160,再重复步骤S1140、S1150和S1160。
在步骤S1160中,如果t的值表示的数目等于输出样点值的所期望数目OLENGTH,,便执行步骤S1180。
在步骤S1180中,输出信号O、SMDP被传递到输出260上(见图9和/或图16)。
如上所述,表示从轴旋转发出机械振动的数据信号可以包括重复信号特征,因此一定的信号特征可以在被监视轴的每转中重复一定数目的次数。不仅如此,也可能出现几种互相不同的重复信号特征,其中互相不同的重复信号特征可以具有互相不同的重复频率。正如以上介绍的增强信号中重复信号特征的方法,能够有益地同时检测具有互相不同重复频率的许多重复信号特征。优选情况下,这使得在单一测量和分析周期内,能够同时检测比如轴承内环故障特征和轴承外环故障特征,正如下面的介绍。
图13是包括两个重复信号特征4010和4020的实例输出信号SMDP的示意展示。输出信号SMDP可以比图13中的展示包括更多的重复信号特征,但是为了展示目的仅仅显示了两个重复信号特征。图13仅仅显示了重复信号特征4010和4020的许多数字值的某些。
在图13中,展示了外环(OR)频率信号4020和内环(IR)频率信号4010。正如在图13中可见,外环(OR)频率信号4020具有低于内环(IR)频率信号4010的频率。外环(OR)频率信号4020和内环(IR)频率信号4010的重复频率分别是1/TOR和1/TIR。
在以上介绍的为增强重复信号模式而运行增强器320的方法实施例中,在步骤S1050中计算输出信号时,放大了重复信号模式。如果在步骤S1010中给因子L更高的值就比给因子L更低的值达到了重复信号模式更高的放大。较高的L值意味着在步骤S1030中要求较长的输入信号ILENGTH。所以,较长的输入信号ILENGTH会引起输出信号中重复信号模式更高的放大。因此,关于输出信号中的重复信号模式,较长的输入信号ILENGTH呈现出更好地衰减随机信号的效果。
根据本发明的实施例,整数值ILENGTH可以响应所期望的随机信号衰减量而选择。在这样的实施例中,可以根据所选定的整数值ILENGTH确定长度因子L。
现在考虑为增强重复信号模式而运行增强器320的方法的示范实施例,其中本方法被用于放大具有某一个最低频率的重复信号模式。为了能够分析具有所述某一个最低频率的重复信号模式,需要一定长度的输出信号。
如上所述,在输出信号的计算中使用更长的输入数据信号与使用更短的输入数据信号相比,会引起重复信号模式放大得更多。如果需要重复信号模式的一定的放大,所以有可能使用一定长度的输入信号,以便实现重复信号模式的这种一定的放大。
为了展示上述实施例,考虑以下实例:
具有最低重复频率fI的重复信号模式受到关注。为了确保检测出这样的重复信号,将有必要生产能够表明完整周期的输出信号,也就是说,需要表现的持续时间为TI=1/fI。当连续的输出信号样点值被样点周期tdelta分开时,输出信号中样点值的最小数目将为OLENGTHmin=TI/tdelta。
如上所述,重复信号的放大量将随着输入信号的长度增加。
如上所述,运行以上参考图10至图13介绍的方法是为了增强从旋转轴发出的测量数据序列中的重复信号特征。措词“重复信号特征”应当理解为样点值[x(t)、x(t+T)、x(t+2T)、...、...x(t+nT)]包括具有非随机振幅值的振幅分量,并且其中这些样点值之间的持续时间T是不变的,只要该轴以不变的转速旋转。参考图13,应当理解,数字值4010由增强输入信号I中多个重复信号值而产生(见图11),其中输入信号值在时间上被持续时间TIR分离。因此在这种情况下,能够推断“重复信号特征”与轴承组件内环处的故障有关,此时重复周期TIR对应于滚珠在内环的通过速率。当然这假设知道了轴的直径和转速。同样,当存在着这样的“重复信号特征”信号分量时,就可能有重复信号分量值x,使得x(t)的振幅类似于x(t+T),它的振幅类似于x(t+2T),它的振幅类似于x(t+nT),等等。当在输入信号中存在着这样的“重复信号特征”时,作为有益的情况,即使该重复信号特征微弱到产生的振幅分量小于随机信号分量的振幅分量时,使用以上介绍的方法也可以检测到它。
当处理器50执行以上连同图4讨论的对应程序代码94时,分析装置14可以执行连同图10至图13介绍的方法。数据处理器50可以包括中央处理单元,用于控制分析装置14以及数字信号处理器(DSP)的运行。DSP可以被安排为实际运行程序代码90,使分析装置14执行程序94,使以上连同图10至图13介绍的过程得以执行。数字信号处理器可以是如TMS320C6722类型,由Texas Instruments制造。以这种方式,分析装置14可以运行以便执行全部的信号处理功能94,包括滤波功能240,包络功能250、降采样功能310和470以及增强功能320。
根据本发明的另一个实施例,信号处理可以在装置14与计算机33之间共享,如上所述。因此,装置14可以接收模拟测量信号SEA并产生对应的数字信号SMD,然后将数字信号SMD传递到控制计算机33,允许在控制位置31执行进一步的信号处理功能94。
图10C展示了增强器的实施例320。增强器320具有输入端315,它由此可以接收具有采样频率fSRED的数字信号SRED。
增强器320可以包括信号处理器325,适于在端口326上接收数字信号SRED。信号处理器325还包括端口327,用于接收指示输入信号I的期望长度ILENGTH的控制值。
增强器320还包括参数设置装置330。参数设置装置330用于产生若干相关控制值,用于执行期望的信号增强。因此参数设置装置330具有输出332,用于向信号处理器325传递控制值ILENGTH。
增强器320可以在输入端335上接收设置指令。在输入端335上接收的设置指令可以包括指示级别值的数据Y、指示频率分辨率的数据Z和指示SNR改进器值的数据L。输入端335可以被连接以向参数设置装置330传递所接收的数据(见图10C)。
增强器320可以整合在以上如参考图1介绍的分析装置14中。作为替代,增强器320也可以是中央控制位置31处的控制计算机33(见图1)的一部分。所以具有采样频率fSRED的数字信号SRED可以如从分析装置14经由通信网络18在端口29B上传递到控制计算机33。
图10D展示了根据增强器方法实施例的若干信号。在图10D的顶部示意地展示了具有采样频率fSR的数字输入信号I。数字输入信号I包括至少ILENGTH个样点值,其中ILENGTH是正整数。
如同以上由等式介绍的计算的执行可以如包括第一信号部分S1和第二信号部分S2的运算所展示。
第一信号部分S1包括输入信号I中前S1L个样点值的副本。S1L=ILENGTH-OLENGTH
第二信号部分S2包括输入信号I中后S2L个样点值的副本。S2L=ILENGTH-SSTART
在图10D底部的信号OS1是响应包括第一信号部分S1和第二信号部分S2的计算获得的输出信号OS1的示意展示。
图10E展示了操作增强器320的方法实施例。在步骤S310用户接口24B、102、104提示用户输入若干增强器设置值。根据某实施例,用户接口适于请求用户指示期望的频率分辨率Z;以及期望的要被检测的最高重复频率fDmax,以及指示期望的信噪比改进的信息。指示期望的SNR改进的信息可以以SNR改进器值L的形式输入。期望的最高重复频率可以以级别数OVHigh,Y的形式输入。在这种语境中,级别数OVHigh,Y等于要被检测的最高重复频率(fDmax)与所述旋转速度(fROT)之间的关系(Y,OV,OVHigh):
OVHigh=Y=fDmax/fROT
根据本发明的实施例,用户接口24B、102、104提示用户输入期望的频率分辨率Z(步骤S310),随后它适于等待以指示期望的频率分辨率Z的数据的形式的输入或者以指令增强器320自动地设置频率分辨率Z的数据的形式的输入(步骤S320)。如果用户输入了指示期望的频率分辨率Z的数据,那么输入的数据将被传递到参数设置器330(步骤S330)。如果用户输入的数据指示期望让频率分辨率Z被自动设置,用户接口将向参数设置器330指示把频率分辨率Z设置为默认值(步骤S340)。
随后,用户接口24B、102、104提示用户输入期望的要被检测的最高重复频率(步骤S350)。用户接口24B、102、104然后适于等待指示期望的最高重复频率的输入或者以指令增强器320自动地设置最高重复频率的数据的形式的输入(步骤S360)。如果用户输入了指示期望的最高重复频率的数据,那么输入的数据将被接收和传递到参数设置器330(步骤S370)。如果用户输入的数据指示期望让最高重复频率被自动设置,用户接口将向参数设置器330指示把最高重复频率设置为默认值(步骤S380)。最高重复频率可以以级别数OVHigh,Y的形式输入和/或设置。正如以上讨论,级别数OVHigh,Y是在要产生的输出信号OS中可检测的最高重复频率的设置。当例如所关注轴承信号可以发生大约被监视轴8、801A、801B、801C、803每转y次时,级别数OVHigh,Y应当被设置为至少y。使用数字,这意味着当所关注轴承信号可以发生大约被监视轴8、801A、801B、801C、803每转100次时,级别数OVHigh,Y应当被设置为至少100。
随后,用户接口24B、102、104提示用户输入期望的SNR改进(步骤S390)。用户接口24B、102、104然后适于等待指示期望的SNR改进的输入或者以指令增强器320自动地设置SNR改进的数据的形式的输入(步骤S400)。如果用户输入了指示SNR改进的数据,那么输入的数据将被接收和传递到参数设置器330(步骤S410)。如果用户输入的数据指示期望让SNR改进被自动设置,用户接口将向参数设置器330指示把SNR改进设置为默认值(步骤S420)。SNR改进可以以SNR改进器值L的形式输入。
这些控制值在端口335上传递到参数设置器320(图10C和图10E)。
参数设置器330包括被连接以接收控制值Z和Y的转数计算器340。转数计算器340用于计算值XE。值XE指示输出信号OS1中的样点(见图10D和/或图12和/或图13)应当对应被监视轴8、801A、801B、801C、803的多少“增强转”。例如,当频率分辨率Z被设置为1600而级别值被设置为Y=100时,那么根据本发明的实施例,输出信号OS1中的样点应当对应被监视轴8、801A、801B、801C、803的XE=Z/Y=16“增强转”。因此,转数计算器340具有输入以接收指示频率分辨率设置值Z的数据和指示设置的级别值Y的数据。转数计算器340响应频率分辨率设置值Z和指示设置的级别值Y的数据产生数据值XE。
转数计算器340向输出信号长度计算器345传递数据值XE。
如图10C所示,增强器具有输入端350,用于接收指示在输入端315上收到的信号SRED的采样频率的数据。参考图9、图16和图30,采样频率值fSR可以对应于值fSR1或fSR2。
增强器320还具有输入端360,用于接收指示旋转速度fROT的数据。对于某些机器,旋转速度被预置为常数值,在这样的情况下该速度值可以被提供到输入端360。作为替代,速度检测器420(见图1和图5和图29)可以被提供为传递指示轴8的旋转速度fROT的信号。轴8的旋转速度fROT可以按照转每秒(rps)即赫兹(Hz)提供。
输出信号长度计算器345用于计算指示输出信号OS1中所需样点值数目OLENGTH(见图10D)的值,以响应收到的数据,即响应轴“增强转”值XE、轴旋转速度值fROT和采样频率值fSR。
因此,如果XE=16,采样频率值fSR=30.72Hz以及轴旋转速度值fROT=0.12rps,那么输出信号OS1中所需样点值数目将为OL=XE*fSR/fROT=4096。
因此,根据某实施例,为了能够以频率分辨率Z进行高达级别OVHigh、Y的重复频率后续分析,输出信号OS1、SMDP中所需样点值的最小数目能够根据参数Y、Z、fROT和fSR计算,其中fSR为采样频率值,使得被监视轴8、801A、801B、801C、803每转的样点数目不变。当旋转速度不变时,以及/或者当分数降采样器被用于补偿可变的轴速度时,被监视轴每转的样点数目不变,正如本文档后面更详细的讨论。
所以,要产生的样点值数目OLENGTH必须为OL或更大,其中OL=XE*fSR/fROT。
输出信号长度计算器345用于计算OL值以及设置OLENGTH值。OLENGTH值被设置为等于或大于算出的OL值的值。
输出信号长度计算器345用于向输入信号长度计算器365的输入端364传递OLENGTH值。输出信号长度计算器345也用于向随机信号频率衰减器设置运算器370传递OLENGTH值。随机信号频率衰减器设置运算器370被安排为设置变量SSTART。变量SSTART控制着随机信号衰减的频率限度。
正如从图10D可见,SSTART值除以采样频率fSR对应于时间周期TS:
TS=SSTART/fSR=SSTART*TSR
其中TSR可以表示两个连续样点之间的延续时间。SSTART为延迟的样点数目,即要被相关的信号S1与S2之间的位移,正如在图10D中也可见。
当变量SSTART被设置为与OLENGTH相同的值时,那么具有对应频率fSTOCHASTIC_MAX=1/TS和更高频率的随机信号将在输出信号O、SMDP中被衰减。所以,有益的是把变量SSTART设置为等于OLENGTH的值或大于OLENGTH的值。因此,随机信号频率衰减器设置运算器370被安排为把变量SSTART设置为等于OLENGTH的值或大于OLENGTH的值。
输入信号长度计算器365用于计算IL值和设置ILENGTH值。产生的变量值IL取决于可以在端口366上收到的指示期望的SNR改进的信息、可以在端口367上收到的变量SSTART的值和可以在端口364上收到的OLENGTH值。
因此,为了能够从具有OLENGTH个样点值的增强器320产生输出信号O、SMDP,该增强器必须在端口315上接收至少IL个样点值。根据某实施例,变量IL的值为:
IL=OLENGTH*L+SSTART+OLENGTH
所以,有益的是把变量ILENGTH设置为等于IL的值或大于IL的值。因此,输入信号长度计算器365被安排为把变量ILENGTH设置为等于IL的值或大于IL的值。输入信号长度计算器365被安排为在端口332上把控制值ILENGTH传递到信号处理器325(见图10C)。
输入信号长度计算器365被安排为把控制值ILENGTH传递到求和确定器375的输入端。求和确定器375还具有一个输入用于接收指示变量SSTART的数据。不仅如此,求和确定器375还具有一个输入用于接收指示要以输出信号O、SMDP的形式产生的样点值的数目OLENGTH的数据。所以,随机信号频率衰减器设置运算器370被安排为把指示变量SSTART的数据传递到求和确定器375,而输出信号长度计算器345被安排为把指示OLENGTH值的数据传递到求和确定器375。
求和确定器375被安排为根据数值ILENGTH、OLENGTH和SSTART产生数值CLENGTH。
数值CLENGTH被设置为基本上等于于数值ILENGTH与数值SSTART和OLENGTH之和之间的差的值。因此,求和确定器375可以传递数值CLENGTH=ILENGTH-SSTART-OLENGTH。
求和确定器375可以被安排为在参数设置装置330的输出377上传递数值CLENGTH。
图10G展示了增强器的另一个实施例320,其中求和确定器的实施例375B具有一个输入用于接收指示要产生的样点值的数目OLENGTH的数据,以及另一个输入用于接收SNR改进器值L。求和确定器375B适于根据数值OLENGTH和SNR改进器值L产生数值CLENGTH。求和确定器375B适于设置CLENGTH=L*OLENGTH。
求和确定器375B可以被安排为在参数设置装置330的输出377上传递数值CLENGTH。
输出信号长度计算器345还用于在参数设置装置330的输出379上传递数值OLENGTH。
如上所述,信号处理器325包括端口326用于接收时域数字信号SRED,以及端口327用于接收指示输入信号I的期望长度ILENGTH的控制值。
信号处理器325与具有多个存储器部分的存储器380协作。根据某实施例,该存储器包括存储器部分382用于存储信号SRED的至少ILENGTH个连续样点值。因此,信号处理器325可以响应在激活输入端384上收到激活信号,用于在输入端327上读取ILENGTH值,随后用于在端口326上接收ILENGTH值个连续样点值。信号处理器325还可以响应在激活输入端384上收到激活信号,与存储器380协作用于把这些样点值存储在存储器部分382中。因此,存储器部分382的内容表示输入信号I,正如图10D所描绘。
输出样点值发生器386用于根据输入信号I产生第一信号部分S1和第二信号部分S2。
响应在激活输入端388上收到激活信号,输出样点值发生器386可以用于读取样点I(i0)至I(i0+CLENGTH),并且把这些样点存储在第二存储器部分390中作为信号部分S1;以及输出样点值发生器386可以用于读取样点I(i0+SSTART+1+1)至I(SSTART+CLENGTH+OLENGTH),并且把这些样点存储在第三存储器部分392中,其中i0为不变的正整数。
因此,存储器部分390和392的内容分别可以表示第一信号部分S1和第二信号部分S2,正如图10D所描绘。随后输出样点值发生器386可以用于对信号S1和S2进行互相关。
作为替代,相关包括读取存储器部分382中存储的输入信号I的样点值,正如图10C和图10D的顶部示意地展示。参考图10F,输出样点值发生器386可以用于执行以下步骤:
步骤S500:把变量t设置为第一个值t0。第一个值t0可以是t0=1。
步骤S510:计算输出样点值:
步骤S520:在输出端口394上传递所产生的样点值SMDP(t)。
步骤S530:增加计数器值t,即设置t:=t+1;
步骤S540:检查数值t是否大于OLENGTH+t0-1。如果数值t大于OLENGTH+t0-1那么产生信号指示已经产生了完整的输出信号(步骤S550)。如果数值t不大于OLENGTH+t0-1那么使用增加的t值重复步骤S510。
图10H是表格,展示了步骤S510中一部分计算的实施例。增强器320适于根据CLENGTH个输入信号乘积P(i,t)产生输出样点值SMDP(t)。CLENGTH为正整数。
参考表1(见图10H),通过将第一样点位置的第一输入样点值I(i)乘以第二输入样点值I(i+t+SSTART)获得输出样点位置t的输入信号乘积P(i,t)。第二输入样点值在输入信号向量I中第二样点位置i+t+SSTART找到(见图10D或图11)。因此,第二输入样点值与所述第一输入样点值分开一定数目NC的样点位置。该一定数目的样点位置可以是NC=(i+t+SSTART)-i=t+SSTART。因此,一定数目NC可以等于输出样点位置值t与一定值SSTART之和。
正如以上连同图10D的说明所指示,一定值SSTART是许多样点位置,它可以对应于时间周期TS。当一定值SSTART被设置为与OLENGTH相同的值时,那么具有对应频率fSTOCHASTIC_MAX和更高频率的随机信号将在输出信号SMDP中被衰减。这个限制频率值fSTOCHASTIC_MAX的值为:
fSTOCHASTIC_MAX=1/TS,其中
TS=SSTART/fSR,其中
fSR=输入信号I的采样频率。
因此,根据优选实施例,一定数目NC等于或大于一定值SSTART。因此,根据优选实施例,输入信号乘积P(i,t)中两项的下标值之间的差等于或大于一定值SSTART。在以上实例中,一项是具有下标i的第一输入样点值I(i),而另一项是具有下标值i+t+SSTART的第二输入样点值I(i+t+SSTART)。在这种连接中重要的是下标值i和i+t+SSTART分别是与输入信号向量I内样点值相关联的值。因此,输入样点值的范围ILENGTH以及下标值i和i+t+SSTART分别必须选择得使这些下标值是输入信号向量内的值。参考图10D的上部,它展示了输入信号向量I的实施例,这意味着下标值i和i+t+SSTART分别必须是从iSTART到iSTART+ILENGTH-1范围内的值。因此,如果常数iSTART被设置为iSTART=1,那么下标值i和i+t+SSTART分别必须是从i=1到i=ILENGTH范围内的值。
输入信号I可以包括ILENGTH个样点值,如上所述。增强器接收具有第一多个ILENGTH输入样点值的输入信号向量。这第一多个ILENGTH输入样点值受到处理以产生具有第二多个OLENGTH输出样点值SMDP(t)的输出信号序列SMDP;所述第二多个为正整数。
根据第三多个CLENGTH输入信号乘积P(i,t),计算输出样点值SMDP(t);所述第三多个(CLENGTH)为正整数。
正如以上等式(1)和(3)中指示,根据某实施例可以使用以下关系:
(1)ILENGTH=OLENGTH*L+SSTART+OLENGTH
和
(3)CLENGTH=ILENGTH-SSTART-OLENGTH
作为实例,可以使用以下数值:
L=10
OLENGTH=1024
SSTART=1024
CLENGTH=10240
ILENGTH=12288
因此,如果例如SSTART=1024,并且t从t=tMIN=1到t=tMAX=OLENGTH=1024变化,并且使用以下等式(8):
那么两项的两个下标值之间的差NC将从NC=1025到NC=2048变化。这是由于最大的下标值差将是NCMAX=SSTART+tMAX=SSTART+OLENGTH=1024+1024=2048;而最小的下标值差将是NCMIN=SSTART+tMIN=1024+1=1025。
因此,如果常数i0=1,那么输入信号向量I必须具有从i=i0=1到i=ILENGTH=12288的下标值范围。
参考图10C,在输出端口394上传递的输出样点值SMDP(t)可以被传递到存储器396,并且存储器396可以存储所接收的样点值,所以它们可读为输出样点值OS1、SMDP的序列,正如图10D左下角示意地展示。
作为替代,在输出样点值发生器386的输出端口394上传递的输出样点值SMDP(t)可以直接传递到增强器320的输出端口398。
根据另一个实施例,产生输出样点值SMDP(t)的等式可以修改为:
以上等式(9)将提供输出信号序列O、SMDP,它与以上等式(5)和(8)产生的输出信号序列O、SMDP等效。可以显示出等式(9)是表达等式(5)的替代方式。
因此,也产生取决于各个输出信号样点值SMDP(t)的等式(9)实施例,要被相乘的样点值之间的样点位置的数目NC将是NC=t+SSTART。
采样频率的降采样
正如以上连同图9的讨论,可能期望提供降采样器310以先降低数字信号的采样频率,再传递到增强器320。有益的是,这样的降采样器310减少了要被分析信号中样点的数目,从而减少了存储要被使用信号所需要的存储空间量。降采样也能够加快在后续增强器320中的处理。
图14A展示了传递给降采样器310的输入的信号中的许多样点值,而图14B展示了对应时间周期的输出样点值。被输入到降采样器310的信号可以具有采样频率fS。正如能够看出,输出信号is具有降低的采样频率fSR1。降采样器310适于执行数字化包络信号SENV的降采样,以便传递具有降低采样频率fSR1的数字信号SRED,使得输出采样频率与输入采样频率fS相比被降低了整数因子M倍。
因此,输出信号SRED仅仅包括了在输入信号SRED中出现的每第M个样点值。图14B展示了M是4的实例,但是M可以是任何正整数。根据本发明的实施例,降采样器可以按照US 5,633,811中的介绍运行,其内容在此引用作为参考。
图15A展示了根据本发明实施例的降采样器310。在根据图15A的降采样器310的实施例310A中,梳状滤波器400对输入信号进行滤波和16∶1比率的降采样。即,与输入采样频率相比,输出采样频率以等于16的第一整数因子M1(M1=16)被降低。有限脉冲响应(FIR)滤波器401接收梳状滤波器400的输出,并且以第二整数因子M2提供了对采样频率的又一次降低。如果整数因子M2=4,FIR滤波器401呈现了采样频率4∶1的降低,所以降采样器310A呈现了64∶1总降采样。
图15B展示了本发明的另一个实施例,其中降采样器310的实施例310B包括低通滤波器402,随后是样点选择器403。样点选择器403适于在从低通滤波器402接收的信号中拾取出每第M个样点值。结果信号SRED1具有fSR1=fS/M的采样频率,其中fS是收到信号SENV的采样频率。低通滤波器402的截止频率由M值控制。
根据一个实施例,M值被预置为一定的值。根据另一个实施例,M值是可设置的。降采样器310可以为可设置的,以进行选定的降采样M∶1,其中M是正整数。可以在降采样器310的端口404上接收M值。
低通滤波器402的截止频率是fSR1/(G*M)赫兹。因子G可以选择为二(2.0)的值,或者大于二(2.0)的值。根据某实施例,G值被选择为2.5与3之间的值。优选情况下,这能够避免假频。低通滤波器402可以由FIR滤波器实施。
由低通滤波器402传递的信号被传递到样点选择器403。样点选择器在一个端口上接收M值,而在另一个端口上接收来自低通滤波器402的信号,并且它响应这些输入产生样点值序列。样点选择器适于在从低通滤波器402接收的信号中拾取出每第M个样点。结果信号SRED1具有fSR1=1/M*fS的采样频率,其中fS是在降采样器310的端口405上接收的信号SENV的采样频率。
可变轴速度的补偿方法
如上所述,作为有益的情况,即使重复信号特征微弱到产生的振幅分量小于随机信号分量时,也可以使用以上介绍的方法检测出输入信号中出现的重复信号特征。
不过,在一定的应用中,轴的转速可能变化。当时,使用输入测量序列执行参考图10-13介绍的方法会导致结果输出信号SMDP的品质退化。
所以,根据本发明一方面的目的是在整个测量序列期间当该轴的转速改变时,实现结果块Y具有与轴的转速不变时同等高的品质。
图16展示了本发明的实施例,包括如以上介绍的降采样器310和增强器320,以及分数降采样器470。
根据本发明的实施例,尽管降采样器310运行以便以M∶1对采样频率降采样,其中M是整数,但是图16的实施例包括分数降采样器470,以U/N对采样频率降采样,其中U和N都是正整数。因此,作为有益的情况,分数降采样器470能够以分数对采样频率降采样。根据某实施例,U和N的值可以在从2到2000的范围内选择。根据某实施例,U和N的值可以在从500到1500的范围内选择。根据又一个实施例,U和N的值可以在从900到1100的范围内选择。
在图16的实施例中,来自降采样器310的输出信号被传递到选择器460。选择器能够对即将输入到增强器320的信号进行选择。对具有不变转速的旋转部件进行状态监视时,可以将选择器460设置在将具有采样频率fSR1的信号SRED传递到增强器320的输入端315的位置,并且可以将禁用分数降采样器470。对具有可变转速的旋转部件进行状态监视时,可以启用分数降采样器470,并且将选择器460设置在将具有采样频率fSR2的信号SRED2传递到增强器320的输入端315的位置。
分数降采样器470具有输入端480。输入端480可以被连接以接收从降采样器310输出的信号。分数降采样器470也具有输入端490,用于接收指示轴8的转速的信息。
可以提供速度检测器420(见图5和图1和图29)以传递信号表明轴8的转速fROT。速度信号可以在处理装置180的端口430上接收,从而使处理装置180能够将该速度信号传递到分数降采样器470的输入端490。可以按照每秒的转数即赫兹(Hz)提供轴8的转速fROT。
图17展示了分数降采样器470的实施例,它能够以分数U/N改变采样频率,其中U和N都是正整数。这使得将被传递到增强器320的采样频率fSR2能够被非常准确地控制,从而即使当轴速度变化时也能够非常好地检测微弱重复信号特征。
在分数降采样器470的输入端490上接收的速度信号被传递到分数发生器500。分数发生器500分别在输出510和520上产生整数输出U和N。U输出被传递到升采样器530。升采样器530经由输入端480接收信号SRED(见图16)。升采样器530包括样点引入器540,用于在端口480上接收的每个样点值之间引入U-1个样点值。每个这样添加的样点值都带有振幅值。根据某实施例,每个这样添加的样点值都被提供了零(0)振幅。
结果信号被传递到低通滤波器550,其截止频率由分数发生器500所传递的值U控制。低通滤波器550的截止频率是fSR2/(K*U)赫兹。因子K可以选择为二(2)的值,或者大于二(2)的值。
结果信号被传递到降采样器560。降采样器560包括低通滤波器570,其截止频率由分数发生器500所传递的值N控制。低通滤波器570的截止频率是fSR2/(K*N)赫兹。因子K可以选择为二(2)的值,或者大于二(2)的值。
由低通滤波器570传递的信号被传递到样点选择器580。样点选择器在一个端口上接收值N,以及在另一个端口上接收来自低通滤波器570的信号,并且它响应这些输入而产生样点值序列。样点选择器适于在从低通滤波器570接收的信号中拾取出每第N个样点值。结果信号SRED2具有fSR2=U/N*fSR1的采样频率,其中fSR1是在端口480上接收的信号SRED的采样频率。结果信号SRED2在输出端口590上传递。
低通滤波器550和570可以由FIR滤波器实施。这有益地排除了使用由样点引入器540引入的零振幅值执行乘法的需要。
图18展示了分数降采样器470的另一个实施例。图18的实施例有益地减少了产生信号SRED2所要求的计算量。
在图18的实施例中,除去了低通滤波器570,所以由低通滤波器550传递的信号被直接传递到样点选择器580。当由硬件实施分数降采样器470时,作为有益的情况,图18的实施例减少了硬件量,从而降低了生产成本。
当由软件实施分数降采样器470时,作为有益的情况,图18的实施例减少了要求执行的程序代码量,从而降低了处理器的负载并提高了执行速度。
参考图17和图18,被传递到分数降采样器470的输出端口上的结果信号SRED2具有fSR2=U/N*fSR1的采样频率,其中fSR1是在端口480上接收的信号SRED的采样频率。分数值U/N取决于在输入端口490上所接收的速率控制信号。如上所述,速率控制信号可以是表明轴8的转速的信号,它可以由速度检测器420(见图1和/或图5)传递。速度检测器420可以由编码器实施,提供的脉冲信号具有适当选择的分辨率,使得速度信号具有所期望的准确度。在一个实施例中,编码器420在轴8的每个整转都传递一个整转标记信号。这样的转标记信号可以是电脉冲的形式,具有能够被准确检测到的边缘并且表明了被监视轴8的一定的旋转位置。根据另一个实施例,编码器420在被监视轴的每转中可以传递许多脉冲信号,以便在该轴的一转内也能够检测出速度变化。
根据某实施例,分数发生器500控制U和N的值,使得降低的采样频率具有使提供的信号SRED2能够有如下性质的值:轴8每转的样点数基本上不变,与轴8的任何速度变化无关。所以,U和N的值越大,分数降采样器470在使轴8每转样点值的数目保持基本上不变方面的能力就越好。
参考图17和图18所介绍的分数降采样可以通过执行对应的方法步骤获得,并且这可以利用存储器60中存储的计算机程序94实现,如上所述。计算机程序可以由DSP 50执行。作为替代,计算机程序也可以由现场可编程门阵列电路(FPGA)执行。
当处理器50执行以上连同图4讨论的对应程序代码94时,分析装置14可以执行连同图10-13介绍的方法以及参考图17和图18介绍的降采样。数据处理器50可以包括中央处理单元50,用于控制分析装置14以及数字信号处理器(DSP)50B的运行。DSP 50B可以被安排为实际运行程序代码90,使分析装置14执行程序94,使以上连同图10至图13介绍的过程得以执行。根据另一个实施例,处理器50B是现场可编程门阵列电路(FPGA)。
图19展示了降采样器310和分数降采样器470的另一个实施例。降采样器310在端口405上接收具有采样频率fS的信号SENV,并且在端口404上接收整数M,如上所述。降采样器310在输出312上传递具有采样频率fSR1的信号SRED1,输出312连接到分数降采样器470A的输入端480。输出采样频率fSR1是:
fSR1=fS/M
其中M是整数。
分数降采样器470A接收具有采样频率fSR1的信号SRED1,作为数据值序列S(j),并且它在其输出590上传递输出信号SRED2作为另一个数据值序列R(q)。
分数降采样器470A可以包括存储器604,适于接收和存储数据值S(j),以及表明被监视旋转部件对应转速fROT的信息。因此存储器604可以存储每个数据值S(j),以便在检测数据值S(j)所对应的传感器信号SEA的值时使它与表明被监视轴转速的值相关联。
当产生输出数据值R(q)时,分数降采样器470A适于从存储器604读取数据值S(j),以及表明对应转速fROT的信息。
从存储器604读取的数据值S(j)被传递到样点引入器540,用于在端口480上接收的每个样点值之间引入U-1个样点值。每个这样添加的样点值都带有振幅值。根据某实施例,每个这样添加的样点值都是零(0)振幅。
结果信号被传递到低通滤波器550,其截止频率由分数发生器500所传递的值U控制,如上所述。
结果信号被传递到样点选择器580。样点选择器在一个端口上接收值N,以及在另一个端口上接收来自低通滤波器550的信号,并且它响应这些输入产生样点值序列。样点选择器适于在从低通滤波器550接收的信号中拾取出每第N个样点值。结果信号SRED2具有fSR2=U/N*fSR1的采样频率,其中fSR1是在端口480上接收的信号SRED的采样频率。结果信号SRED2在输出端口590上传递。
因此,输出数据值R(q)的采样频率fSR2低于输入采样频率fSR1因子D倍。D可以被设置为大于1的任意数,并且它可以是分数。根据若干优选实施例,因子D可设置为1.0至20.0之间的值。在某优选实施例中,因子D是分数,可设置为大约1.3至大约3.0之间的值。将整数U和N设置为适合的值可以获得因子D。因子D等于N除以U:
D=N/U
根据本发明的实施例,整数U和N可被设置为大的整数,以便使因子D=N/U能够以最小的误差跟随速度变化。有益的是,选择变量U和N为大于1000的整数,致使高准确度地调整输出采样频率,以跟踪所监视轴的转速改变。因此,例如将N设置为500而U为1001致使D=1/2.002。
在测量开始时变量D被设置为合适的值,并且该值与要被监视旋转部件的一定转速相关联。此后,在状态监视阶段期间,响应要被监视旋转部件的转速而自动调整分数值D,使得端口590上输出的信号对被监视旋转部件每转提供基本上不变数目的样点值。
如上所述,编码器420在轴8的每个整转都传递一个整转标记信号。这样的整转标记信号可以是电脉冲的形式,具有能够被准确检测到的边缘并且表明了被监视轴8的一定的旋转位置。所述整转标记信号也可以称为索引脉冲,能够在编码器420的输出上生成,以响应在被监视轴旋转时旋转的编码盘片上检测到零角度模式。这能够以几种方式实现,正如本领域的技术人员熟知。该编码盘片可以如配备了零角度模式,它将在盘片的每转产生零角信号。检测速度变化可以如通过被监视轴每次通过该一定旋转位置时在存储器604中登记“整转标记”,以及通过将该“整转标记”与在同一时刻接收的样点值S(j)相关联。以这种方式,当轴旋转较慢时,存储器604将存储两个连续的整转标记之间大量的样点,因为A/D转换器传递每秒不变量的样点。
图20是降采样器310和分数降采样器470又一个实施例的框图。这个分数降采样器实施例被标注为470B。分数降采样器470B可以包括存储器604,适于接收和存储数据值S(j)以及表明被监视旋转部件对应转速fROT的信息。因此存储器604可以存储每个数据值S(j),以便在检测数据值S(j)所对应的传感器信号SEA的值时使它与表明被监视轴转速的值相关联。
分数降采样器470B接收具有采样频率fSR1的信号SRED1,作为数据值序列S(j),并且它在其输出590上传递具有采样频率fSR2的输出信号SRED2作为另一个数据值序列R(q)。
分数降采样器470B可以包括存储器604,适于接收和存储数据值S(j)以及表明被监视旋转部件对应转速fROT的信息。存储器604可以将数据值S(j)存储在若干块中,以便每个块与表明被监视轴相关转速的值相关联,正如以下连同图21的介绍。
分数降采样器470B还可以包括分数降采样变量发生器606,它适于产生分数值D。分数值D可以是浮点数。因此,该分数能够被控制成为响应所接收的速度值fROT的浮点数的值,使得该浮点数的值以一定的误差表明速度值fROT。当如上所述由适宜地编程的DSP实施时,浮点数值的误差可以取决于DSP产生浮点数值的能力。
不仅如此,分数降采样器470B还可以包括FIR滤波器608。FIR滤波器608是低通FIR滤波器,具有一定的低通截止频率,适于利用因子DMAX进行降采样。可以将DMAX设置为合适的值,如20.000。不仅如此,分数降采样器470B还可以包括滤波器参数发生器610。
以下参考图21和图22介绍分数降采样器470B的运行。
图21是流程图,展示了运行图20中降采样器310和分数降采样器470B的方法实施例。
在第一个步骤S2000中,要被监视状态部件的转速fROT被记录在存储器604中(图20和图21),并且这可以与开始测量振动或冲击脉冲在基本上同一时间完成。根据另一个实施例,要被监视状态部件的转速被测量一段时间。检测的最高速度fROTmax和检测的最低速度fROTmin可以被记录在如存储器604中(图20和图21)。
在步骤S2010中,分析了记录的速度值,用于确定转速是否变化的目的。如果确定该速度不变,可以将选择器460(图16)自动设置在将具有采样频率fSR1的信号SRED传递到增强器320的输入端315的位置,并且可以将分数降采样器470、470B禁用。如果确定该速度可变,可以自动启用分数降采样器470、470B,并且将选择器460自动设置在将具有采样频率fSR2的信号SRED2传递到增强器320的输入端315的位置。
在步骤S2020中,用户界面102、106显示了所记录的速度值fROT或速度值fROTmin、fROTmax,并且请求用户输入所期望的级别值OV。如上所述,轴旋转频率fROT常常被称为“级别1”。所关注的信号可能每轴转发生约10次(级别10)。不仅如此,分析某些信号的谐波可能令人关注,所以测量高达级别100或级别500甚至更高的情况可能令人关注。因此,用户可以使用用户界面102输入级别数OV。
在步骤S2030中,确定适当的输出采样频率fSR2。根据某实施例,输出采样频率fSR2被设置为fSR2=C*OV*fROTmin,
其中
C是大于2.0的常数,
OV是表明被监视部件的转速与被分析信号的重复频率之间关系的数,
fROTmin是在即将来临的测量阶段期间,所期望的被监视部件的最低转速。根据某实施例,值fROTmin是在步骤S2020中检测的最低转速,如上所述。
考虑到采样定理,可以选择常数C为2.00(二)或更大的值。根据本发明的实施例,可以把常数C预置为2.40与2.70之间的值。其中,k是具有大于2.0的值的因子。
所以可以选择因子k为大于2.0的值。根据某实施例,作为有益的情况,选择因子C使得100*C/2为整数。根据某实施例,可以把因子C设置为2.56。选择C为2.56,致使100*C=256=2的8次方。
在步骤S2040中,根据被监视部件的所检测的转速fROT选择整数值M。M的值根据被监视部件的所检测的转速可以被自动选择,以使得中间降低采样频率fSR1将高于所期望的输出信号采样频率fSR2。降低的采样频率fSR1的值也根据在测量阶段期间所期望的转速变化程度而选择。根据某实施例,A/D转换器的采样频率fS可以是102.4kHz。根据某实施例,可以将整数值M设置为100与512之间的值,以便呈现1024Hz与100Hz之间的中间降低采样频率fSR1值。
在步骤S2050中,确定了分数降采样变量值D。当被监视状态部件的转速变化时,分数降采样变量值D将根据瞬间检测的速度值变化。
根据步骤S2040和S2050的另一个实施例,设置整数值M时使得中间降低采样频率fSR1高于fSR2(如以上在步骤S2030中确定的)的百分比至少与检测的最高速度值fROTmax除以检测的最低速度值fROTmin之间的关系同样多。根据这个实施例,最大的分数降采样变量值DMAX被设置为值DMAX=fROTmax/fROTmin,而最小的分数降采样变量值DMIN被设置为1.0。此后进行了实际速度值fROT的瞬间实时测量,并且相应地设置了瞬间分数值D。
fROT是表明被监视旋转部件的测量转速的值。
在步骤S2060中,开始实际测量,并且可以确定期望的测量总持续时间。这个持续时间可以根据增强器中需要的随机信号衰减的程度确定。因此,期望的测量总持续时间可以设置为使它对应于或使它超过获得输入信号ILENGTH所需要的持续时间,正如以上连同图10A至图13的讨论。正如以上连同图10A至图13的论述,关于输出信号中的重复信号模式,较长的输入信号ILENGTH致使随机信号的衰减效果更好。
也可以根据期望的被监视部件的转数确定测量的总持续时间。
开始测量时,降采样器310以速率fS接收数字信号SENV,并且它以降低的速率fSR1=fS/M向分数降采样器的输入端480传递数字信号SRED1。以下按照具有样点值S(j)的信号讨论信号SRED1,其中j是整数。
在步骤S2070,将数据值S(j)记录在存储器604中,并且将每个数据值都与转速值fROT相关联。根据本发明的实施例,以每秒fRR=1000次的速率读取和记录转速值fROT。读取和记录速率fRR可以被设置为其他值,取决于被监视旋转部件速度fROT的变化程度。
在后续步骤S2080中,分析已记录的转速值,并且将已记录的数据值S(j)分成取决于转速值的若干块。以这种方式,可以产生许多块数据值S(j)的块,每个数据值S(j)块都与转速值相关联。转速值表明了当记录这个特定块的数据值S(j)时被监视部件的转速。各个数据块的尺寸可以互相不同,即各个块可以持有互相不同数目的数据值S(j)。
例如,假若被监视旋转部件在第一时间阶段期间首先以第一速度fROT1旋转,此后它改变了速度在第二个更短的时间阶段期间以第二速度fROT2旋转,则可以将记录的数据值S(j)分成两个数据块,第一块的数据值与第一速度值fROT1相关联,而第二块的数据值与第二速度值fROT2相关联。在这种情况下,由于第二个时间阶段更短,所以第二个数据块含有的数据值将少于第一个数据块。
根据某实施例,如果已经把全部已记录的数据值S(j)都分成数据块,并且已经将全部块与转速值相关联,那么本方法继续执行步骤S2090。
在步骤S2090中,选择第一块的数据值S(j),并且确定与相关联的转速值fROT所对应的分数降采样值D。将这个分数降采样值D与第一块的数据值S(j)相关联。根据某实施例,如果已经将全部块与对应的分数降采样值D相关联,那么本方法继续执行步骤S2100。因此,根据速度fROT调整了分数降采样值D的值。
在步骤S2100,选择一块数据值S(j)和相关联的分数降采样值D,正如以上在步骤S2090中的介绍。
在步骤S2110中,响应所选择块的输入值S和相关联的分数降采样值D而产生一块输出值R。做到这一点可以如同参考图22的介绍。
在步骤S2120中,检查是否存在着要处理的任何剩余输入数据值。如果存在着要处理的另一块输入数据值,那么重复步骤S2100。如果不存在要处理的剩余输入数据值块,那么该测量阶段完成。
图22A、22B和图22C展示了运行图20中分数降采样器470B方法实施例的流程图。
在步骤S2200中,接收了一块输入数据值S(j)和相关联的特定分数降采样值D。根据某实施例,所接收的数据如同以上在图21步骤S2100中的介绍。在所接收的输入数据值S块中的输入数据值S(j)都与特定分数降采样值D相关联。
在步骤S2210到S2390中,对于在步骤S2200中接收的特定分数降采样值D调整FIR滤波器608,并且产生了一组对应的输出信号值R(q)。
在步骤S2210,选择适于特定分数降采样值D的滤波器设置。正如以上连同图20所述,FIR滤波器608是低通FIR滤波器,具有适于用因子DMAX进行降采样的一定的低通截止频率。可以把因子DMAX设置为适合的值,如20。根据因子DMAX和在步骤S2200中接收的特定分数降采样值D,设置滤波器比值FR的值。步骤S2210可以由滤波器参数发生器610(图20)执行。
在步骤S2220中,选择所接收的输入数据块s(j)中的起始位置值x。应当指出,起始位置值x不需要是整数。FIR滤波器608具有长度FLENGTH,然后将根据滤波器长度FLENGTH和滤波器比值FR选择起始位置值x。滤波器比值FR是在以上步骤S2210中设置的值。根据某实施例,起始位置值x可以被设置为x:=FLENGTH/FR。
在步骤S2230中,准备滤波器求和的值SUM,并且设置初始值,比如SUM:=0.0。
在步骤S2240中,选择在所接收的输入数据中位置x邻近和前方的位置j。可以把位置j选择为x的整数部分。
在步骤S2250中,选择FIR滤波器中的位置Fpos,它对应于所接收的输入数据中选定的位置j。位置Fpos可以是分数。关于滤波器的中间位置,可以把滤波器位置Fpos确定为:
Fpos=[(x-j)*FR]
其中FR是滤波器比值。
在步骤S2260中,检查所确定的滤波器位置值Fpos是否在允许的界限值之外,即在滤波器之外某位置处的某些点。如果这种情况发生,那么继续进行下面的步骤S2300。否则继续进行步骤S2270。
在步骤S2270中,利用插值计算滤波器的值。应当指出,FIR低通滤波器中的相邻滤波器系数值一般具有类似的数字值。因此,作为有益的情况,插值的值将是准确的。首先计算了整数位置值IFpos:
IFpos:=Fpos的整数部分
位置Fpos的滤波器值Fval将是:
Fval=A(IFpos)+[A(IFpos+1)-A(IFpos)]*[Fpos-Ifpos]
其中A(IFpos)和A(IFpos+1)是参考滤波器中的值,而滤波器位置Fpos是这些值之间的位置。
在步骤S2280中,响应信号位置j计算滤波器求和值SUM的更新:
SUM:=SUM+Fval*S(j)
在步骤S2290中,移动到另一个信号位置:
设置j:=j-1
此后,到步骤S2250。
在步骤S2300中,选择与位置x邻近和随后的所接收的输入数据中的位置j。这个位置j可以选为x的整数部分加上1(一),即j:=1+x的整数部分。
在步骤S2310中,选择FIR滤波器中的位置,它对应于所接收的输入数据中选定的位置j。位置Fpos可以是分数。关于滤波器的中间位置,可以把滤波器位置Fpos确定为:
FPOS=[(j-x)*FR]
其中FR是滤波器比值。
在步骤S2320中,检查所确定的滤波器位置值Fpos是否在允许的界限值之外,即在滤波器之外某位置处的某些点。如果这种情况发生,那么继续进行下面的步骤S2360。否则继续进行步骤S2330。
在步骤S2330中,利用插值计算滤波器的值。应当指出,FIR低通滤波器中的相邻滤波器系数值一般具有类似的数字值。因此,作为有益的情况,插值的值将是准确的。首先计算了整数位置值IFpos:
IFpos:=Fpos的整数部分
位置Fpos的滤波器值将是:
Fval(Fpos)=A(IFpos)+[A(IFpos+1)-A(IFpos)]*[Fpos-Ifpos]
其中A(IFpos)和A(IFpos+1)是参考滤波器中的值,而滤波器位置FPOS是这些值之间的位置。
在步骤S2340中,响应信号位置j计算滤波器求和值SUM的更新:
SUM:=SUM+Fval*S(j)
在步骤S2350中,移动到另一个信号位置:
设置j:=j+1
此后,到步骤S2310。
在步骤S2360中,传递输出数据值R(j)。可以把输出数据值R(j)传递到存储器,以便在连续的存储器位置中存储连续的输出数据值。
R(j):=SUM
在步骤S2370中,更新位置值x:
x:=x+D
在步骤S2380中,更新位置值j:
j:=j+1
在步骤S2390中,检查是否已经产生了所期望数目的输出数据值。如果尚未产生所期望数目的输出数据值,那么去往步骤S2230。如果已经产生了所期望数目的输出数据值,那么去往关于图21所介绍方法中的步骤S2120。
实际上,设计步骤S2390是为了确保产生输出信号值R(q)的块,与在步骤S2200中接收的输入数据值S的块相对应,并且当已经产生了输入数据值S对应的输出信号值R时,那么应当执行图21中的步骤S2120。
参考图22介绍的方法可以实施为计算机程序的子程序,并且步骤S2100和步骤S2110可以实施为主程序。
根据本发明的又一个实施例,通过控制由时钟190传递的时钟频率,可以实现对可变轴速的补偿。如上所述,可以提供速度检测器420(见图5)以传递表明轴8转速fROT的信号。速度信号可以在处理装置180的端口430上接收,从而使处理装置180能够控制时钟190。所以,处理装置180可以具有传递时钟控制信号的端口440。因此,处理装置180可以适于响应所检测到的转速fROT而控制时钟频率。
正如连同图2B所述,A/D转换器的采样频率取决于时钟频率。因此,装置14适于响应所检测到的转速fROT而控制时钟频率,以便即使转速变化也使被监视旋转部件每转的样点值数目基本上保持不变。
根据本发明的再一个实施例,通过US 7,010,445中介绍的产生自相关数据的方法,可以实现增强器功能320、94,其内容在此引用作为参考。确切地说,数字信号处理器50可以包括功能94,对数字化的信号执行连续的傅氏变换操作以提供自相关数据。
监视齿轮系统的状态
应当指出,本发明的实施例也可以用于测量、监视和检测齿轮系统的状态。监视包括周转圆传动、齿轮和/或齿轮箱的周转圆齿轮系统时,某些实施例提供了特别有益的效果。这将在下面更详细地介绍。周转圆传动、齿轮和/或齿轮箱也可以被称为行星传动、齿轮和/或齿轮箱。
图23是展示行星齿轮系统700的正视图。行星齿轮系统700包括至少一个或多个外部齿轮702、703、704,围绕着中心齿轮701旋转。外部齿轮702、703、704通常被称为行星齿轮,而中心齿轮701通常被称为太阳齿轮。行星齿轮系统700也可以使用外面的内齿圈705,通常也称为内齿环。行星齿轮702、703、704可以包括P个齿707,太阳齿轮701可以包括S个齿708,而内齿环705可以包括A个齿706。内齿环705上的A个齿被安排与行星齿轮702、703、704上的P个齿啮合,它们又被安排与太阳齿轮701上的S个齿啮合。不过应当指出,太阳齿轮701一般大于行星齿轮702、703、704,从而图23所示的展示不应当解释为是这个方面的限制。当太阳齿轮701与行星齿轮702、703、704的尺寸不同时,分析装置14也可以区分行星齿轮系统700的不同轴和齿轮的所检测的状态,从以下的叙述将变得显而易见。
在许多行星齿轮系统中,太阳齿轮701、行星齿轮702、703、704或内齿环705这三种基本组件之一被保持固定。两个剩余组件之一然后可以用作输入并向行星齿轮系统700提供动力。最后剩余的组件然后可以用作输出并从行星齿轮系统700中接收动力。输入转数与输出转数的比值取决于每个齿轮中的齿数,以及哪个组件被保持固定。
图24是图23中行星齿轮系统700的示意侧视图,从图23中箭头SW的方向所见。包括行星齿轮系统700的示范组合体800可以包括至少一个传感器10和至少一台根据以上介绍的本发明的分析装置14。组合体800例如可用于风力涡轮的齿轮箱。
在组合体800的实施例中,内齿环705保持固定。可旋转轴801具有多个可移动的臂或支架801A、801B、801C,被安排与行星齿轮702、703、704连接。向可旋转轴801提供了输入旋转802后,可旋转轴801和可移动臂801A、801B、801C以及行星齿轮702、703、704可以用作输入并向行星齿轮系统700提供动力。可旋转轴801和行星齿轮702、703、704然后可以相对于太阳齿轮701旋转。可以被安装在旋转轴803上的太阳齿轮701然后可以用作输出并从行星齿轮系统700接收动力。这种配置将产生齿轮比G=1+A/S的增加。比如用作风力涡轮中的齿轮箱时,齿轮比G可以安排为使得输出旋转大约是输入旋转的5-6倍。行星齿轮702、703、704可以分别经由轴承7A、7B和7C安装在可移动臂或支架801A、801B、801C上(在图23和图24中都有显示)。可旋转轴801可以安装在轴承7D上。同样,旋转轴803可以安装在轴承7E上,而太阳齿轮701可以经由轴承7F安装在旋转轴803上。
根据本发明的实施例,所述至少一个传感器10可以连接到行星齿轮系统700的固定内齿环705的测量点12上。传感器10也可以被安排为与分析装置14通信。分析装置14可以被安排为根据传感器10传递的测量数据或信号值,分析行星齿轮系统700的状态,如本文档中以上的介绍。分析装置14可以包括如以上的鉴别器230。
图25展示了响应由所述至少一个传感器10在组合体800中行星齿轮系统700旋转后所检测的信号,由预处理器200(见图5或图16)产生并输出的示范信号的模拟版本。显示信号的持续时间为TREV,它表示在可旋转轴801的一转期间所检测出的信号值。应当理解,由预处理器200在端口260(见图5和图16)上传递的信号可以被传递到鉴别器230的输入220(见图8或图7)。
正如从图25可见,当行星齿轮702、703、704的每一个通过组合体800中传感器10的测量点12时,信号的振幅或信号输出增大。信号的这些部分以下被称为高振幅区域702A、703A、704A,它们可以包括高振幅尖峰901。也能够显示出,在可旋转轴801的一转,即时间周期TREV上信号中尖峰901、902的总量,直接与内齿环705上齿数相关。例如,假若内齿环705上齿数A=73,时间周期TREV期间信号中尖峰的总数将是73;或者假若内齿环705上齿数A=75,时间周期TREV期间信号中尖峰的总数将是75,等等。已经显示出这是真实的,只要组合体800的齿轮702、703、704、705中不存在误差或故障。
图26展示了图25所示信号的高振幅区域702A的一部分的实例。当行星齿轮702通过其机械地离测量点12和传感器10最近的位置时(见图23和图24)可能产生这个信号部分。已经注意到,不时地会出现图26中展示的小周期扰动或振动903。注意,已经将这些小周期扰动903联系到图23和图24所示轴承7A中误差、故障或裂痕的出现,轴承7A可以被安装到移动臂之一801A。小周期扰动903从而可以通过行星齿轮系统700的行星齿轮702从轴承7A传播(或转化)到内齿环705,以上连同图1至图24介绍的传感器10在此可以拾取小周期扰动903。同样,安装到可移动臂801B或801C之一的轴承7B或7C中的误差、故障或裂痕也可能产生这样的小周期扰动903,传感器10可以用与以上相同的方式拾取它们。也应当注意,小周期扰动903也可能从可以安装到旋转轴803的轴承7F中的误差、故障或裂痕发出。在信号中检测出这些小周期扰动可以指示轴承7A、7B、7C和/或7F开始退化,或者指示它们处于其现行预期使用寿命的极限。例如,这可能是重要的,因为它可以有助于预测行星齿轮系统700和/或组合体800何时需要维护或替换。
根据本发明的实施例,分析装置14的鉴别器230中的状态分析器290可以安排为检测从传感器10收到信号中的这些小周期扰动903。利用先前介绍的本发明的实施例使得这成为可能。小周期扰动903也可以被称为冲击脉冲903或振动903。根据本发明的实施例,采用以上介绍的增强器320的分析装置14,使用以上介绍的安装在内齿环705上的传感器10能够检测源自轴承7A(或7B、7C或者7F)的这些冲击脉冲903或振动903。即使机械冲击脉冲或振动信号已经经由一个或几个行星齿轮702、703或704传播,尽管由连接到内齿环705上传感器10所拾取的机械冲击脉冲或振动信号可能微弱,但是提供以上介绍的增强器320使得监视轴承7A(或7B、7C或者7F)的状态成为可能。
正如先前论述和图7-9显示,状态分析器290可以安排为通过操作时域中的信号或者频率域中的信号,执行适当的分析。不过,检测从传感器10收到信号中的小周期扰动903最适于在频率域中描述,如图27所示。
图27展示了包括图26所示小周期扰动903的信号的示范频谱。该信号的频谱包括频率尖峰904,它与行星齿轮702、703、704和内齿环705的齿接合或啮合直接相关。事实上,频谱中尖峰904的频率将位于A×Ω,其中
A是内齿环705的总齿数,以及
Ω是可旋转轴801的每秒转数,这时以不变的转速发生旋转802。
除了频谱中的尖峰904,图26展示的小周期扰动903还可能在以大约频谱中尖峰904为中心的频率f1、f2处产生尖峰904、906。在频率f1、f2处的尖峰904、906从而也可以被称为中心尖峰904周围的对称边带。根据本发明的示范实施例,状态分析器290可以安排为检测频谱中的一个或几个尖峰,从而安排为检测从传感器10收到信号中的小周期扰动。也能够显示出,根据等式1-2,频率f1、f2处的尖峰905、906与中心尖峰904有关:
f1=(A×Ω)-(fD×f702) (等式1)
f2=(A×Ω)+(fD×f702) (等式2)
其中
A是内齿环705的总齿数,
Ω是可旋转轴801的每秒转数,以及
fD是可以表示已退化状态的重复信号特征的重复频率,以及
f702是行星702围绕其自身中心的每秒转数。
重复信号特征的重复频率fD指示了重复信号特征起源的旋转部件之一。重复信号特征的重复频率fD还能够用于区分不同类型的退化状态,正如以上如连同图8的讨论。所以,检测出的重复信号特征的重复频率fD可以表示轴承保持器损坏频率(FTF)、轴承滚动件损坏(BS)频率、外环(OR)频率或内环(IR)频率,它们与图24中组合体800的行星齿轮系统700中轴承7A、7B、7C或7F有关。
因此,正如以上介绍,从一个或几个轴比如可旋转轴801和/或旋转轴803(见图23至图24)的旋转发源的表示机械振动的数据信号可以包括几个重复信号特征,从而一定的信号特征可以在被监控轴之一的每转重复一定的次数。不仅如此,还可能出现几个互相不同的重复信号特征,其中互相不同的重复信号特征可以具有互相不同的重复频率。增强信号中重复信号特征的方法,如以上介绍,作为有益的情况,能够同时检测具有互相不同重复频率的许多重复信号特征。优选情况下,这使得能够使用单一传感器10同时监视与不同轴801、803相关联的几个轴承7A、7B、7C、7F。同时监视还可以使用以下事实:太阳齿轮701和行星齿轮702、703、704的尺寸正常情况是不同的尺寸,这进一步可以实现容易地检测出正在产生小周期扰动903的是图23至图24中轴承7A、7B、7C、7F的哪一个,从而图23至图24中轴承7A、7B、7C、7F的哪一个可能需要维护或替换。增强信号中重复信号特征的方法,如以上介绍,作为有益的情况,也使得有可能在单一信号测量和分析周期中,区分如轴承内环故障特征和轴承外环故障特征。
Ω的相关值表示行星齿轮702、703、704的转速,它能够由传感器420(见图24)指示。传感器420可以适于产生表明轴803关于内齿环705旋转的信号,并且在已知内齿环705、行星齿轮702、703、704和太阳齿轮701的齿数时,从这个信号可以计算出Ω的相关值。
图28展示了图25所示的示范信号的一部分的实例。这个示范部分演示了状态分析器290也可以安排为以类似于以上介绍的方式进行检测的误差或故障的另一个实例。如果齿轮701、702、703、704、705的一个或几个中的齿将破裂或实质上已磨损,状态分析器290可以安排为检测齿已破裂或磨损,因为这也将产生周期扰动,即由于缺失或磨损的齿缺乏接合或啮合。状态分析器290在例如从传感器10接收的信号的频谱中可以检测到这种情况。也应当指出,这种类型的误差或故障可以由状态分析器290在任何类型的齿轮和/或齿轮系统中检测出。这种类型的齿接合误差或啮合误差的频率,在齿轮和/或齿轮系统中往往处于显著高于例如图27中频率f1、f2的频率。
图29展示了根据本发明实施例的状态分析系统2的又一个实施例。传感器10在物理上与机器6相关联,机器6可以包括具有多个旋转部件的齿轮系统700(见图1和图29)。图29的齿轮系统可以是图24的行星齿轮系统700。行星齿轮系统700可以例如被用作风力涡轮的齿轮箱。
传感器单元10可以是冲击脉冲测量传感器,适于产生的模拟信号SEA包括取决于齿轮系统700中可旋转运动部件振动的振动信号分量。传感器10将模拟信号SEA传递给信号处理组合体920。
信号处理组合体920可以包括传感器接口40和数据处理装置50。传感器接口40包括产生数字测量信号SMD的A/D转换器44(图2A、图2B)。A/D转换器44连接到数据处理装置50,以便将数字测量数据信号SMD传递给数据处理装置50。
数据处理装置50连接到用户接口102。用户接口102可以包括用户输入装置104,让用户能够提供用户输入。这样的用户输入可以包括选择所期望的分析功能105、290、290T、290F(图4、图7、图8),以及/或者设置信号处理功能94、250、310、470、470A、470B、320、294(见图4,图30)
用户接口102还可以包括显示器单元106,如连同图2A和图5的介绍。
图30是框图,与用户接口102、104和显示器106一起展示了图29的信号处理组合体920的若干部件。
传感器接口40包括从冲击脉冲测量传感器接收模拟信号SEA的输入端42和A/D转换器44。还可以可选地提供信号调节器43(图2B)。A/D转换器44以一定的采样频率fS对接收的模拟信号进行采样,以便传递具有一定的采样频率fS的数字测量数据信号SMD。
采样频率fS可以设置为:
fS=k*fSEAmax
其中
k是具有大于2.0的值的因子。
所以可以选择因子k为大于2.0的值。优选情况下,可以选择因子k为2.0与2.9之间的值,以便避免假频效应。如在本文档以上所述,选择因子k为大于2.2的值提供了关于假频效应的安全余量。可以选择因子k为2.2与2.9之间的值,以便提供所述安全余量,同时避免不必要地产生许多样点值。根据某实施例,作为有益的情况,选择的因子k使得100*k/2呈现为整数。根据某实施例,因子k可以被设置为2.56。选择k为2.56致使100*k=256=2的8次方。
根据某实施例,数字测量数据信号SMD的采样频率fS可以固定到一定的值fS,比如fS=102.4kHz。
因此,当采样频率fS被固定到一定的值fS时,模拟信号SEA的最高频率fSEAmax将为:
fSEAmax=fS/k
其中fSEAmax是采样后信号中要被分析的最高频率。
因此,当采样频率fS被固定到一定的值fS=102.4kHz,而因子k被设置为2.56时,模拟信号SEA的最大频率fSEAmax将为:
fSEAmax=fS/k=102400/2.56=40kHz
滤波器240接收采样频率为fS的数字测量数据信号SMD。根据某实施例,滤波器240是截止频率为fLC的高通滤波器。这个实施例通过用高通滤波器240替换连同图6介绍的带通滤波器,简化了设计。高通滤波器240的截止频率fLC选择为近似于谐振冲击脉冲测量传感器10的期望机械谐振频率值fRMU的最低值。当机械谐振频率fRM在30kHz到35kHz范围内的某处时,高通滤波器240可以被设计为较低截止频率fLC=30kHz。然后高通滤波后的信号被传递给整流器270并继续到低通滤波器280。
根据某实施例,应当有可能使用谐振频率在20kHz到35kHz范围内某处的传感器10。为了实现这个目标,高通滤波器240可以被设计为截止频率fLC=20kHz。
数字滤波器240的输出信号被传递到数字包络器250。
鉴于响应测量信号产生包络信号的现有技术的模拟设备采用模拟整流器,该模拟整流器会固有地导致结果信号中引入的系统误差,而作为有益的情况,数字包络器250将产生真正的整流结果而没有任何系统误差。所以,数字包络信号SENV将具有高信噪比,因为在数字滤波器240的带宽内谐振频率处机械谐振的传感器会引起信号高振幅。不仅如此,在数字域中执行的信号处理不增加噪声,从而不增加系统误差。
根据本发明的实施例,可以除去包络器250中的可选低通滤波器280。实际上,由于降采样器310包括低通滤波器的功能,所以除去了包络器250中的可选低通滤波器280。因此,图30的包络器250有效地包括数字整流器270,并且由数字整流器270产生的信号被传递给整数降采样器310,它包括低通滤波。
整数降采样器310适于进行数字包络信号SENV的降采样,以便传递具有降采样频率fSR1的数字信号SRED,使得输出采样频率与输入采样频率fS相比被降低了整数因子M倍。
值M可以根据检测出的转速fROT设置。降采样器310可设置为使得所选的降采样为M∶1,其中M是正整数。值M可以在降采样器310的端口404上收到。
优选地,使用低通有限脉冲响应滤波器在多个步骤中执行整数降采样,其中每个FIR滤波器都可设置为所期望的降采样程度。与在多个滤波器中进行降采样相关联的优点在于,只有最后的滤波器需要具有陡斜坡。陡斜坡的FIR滤波器固有地必须具有许多阶,即陡的FIR滤波器必然是长的滤波器。FIR阶数指示了以下信息:
1)实施该滤波器所需要的存储器量,
2)需要的计算量,以及
3)滤波器能够进行的“滤波”量;实际上,更高阶意味着抑止频带的衰减更强,波纹更少、滤波器更窄,等等。因此滤波器越短,它就能够由DSP 50执行得越快。FIR滤波器的长度也与可实现的降采样程度成正比。所以,根据整数降采样器的实施例,在不止两个步骤中执行降采样。
根据优选实施例,在M1、M2、M3和M4四个步骤中进行整数降采样。总的降采样等于M1*M2*M3*M4。这可以通过提供一族不同的FIR滤波器实现,它们可以在几种组合中结合以实现所期望的总降采样M。根据某实施例,在所述族中有八种不同的FIR滤波器。
作为有益的情况,在最后第四个步骤中最大的降采样程度是五(M4=5),致使具有仅仅201阶的适当短的滤波器。以这种方式能够允许步骤1、2和3中的FIR滤波器具有甚至更低的阶。事实上这允许在步骤1、2和3中的每个FIR滤波器具有71阶或更低。为了实现总降采样M=4000,有可能选择提供降采样M1=10、M2=10和M3=10的三个FIR滤波器,以及提供降采样M4=4的FIR滤波器。这致使输出采样频率fSR1=25.6,此时fS=102400Hz,而频率范围为10Hz。这四个FIR滤波器将总共具有414阶,而且抑止频带衰减的结果非常好。事实上,要是在仅仅一个单步骤中进行M4=4000的降采样,它可能将需要大约160000阶以实现同样好的抑止频带衰减。
整数降采样器310的输出312被连接到分数降采样器470以及选择器460的输入端。选择器使得选择的信号能够被输入到增强器320。
当对具有不变转速的旋转部件进行状态监视时,可以把选择器460设置在将采样频率为fSR1的信号SRED传递到增强器320的输入端315的位置,并且可以禁用分数降采样器470。当对具有可变转速的旋转部件进行状态监视时,可以启用分数降采样器470,并且可以把选择器460设置在将采样频率为fSR2的信号SRED2传递到增强器320的输入端315的位置。
分数降采样器470可以由包括连同图20、21和22以及图4介绍的可调整FIR滤波器608的分数降采样器470B、94实施。
分数降采样器470被连接以向选择器460传递具有更低采样频率fSR2的降采样后信号SRED2,所以当状态分析器被设置为监视具有可变转速的机器时,来自分数降采样器470B的输出被传递到增强器320。
可以如连同图10A、10B、11、12和13以及图4的介绍,来实施增强器320、94。输入到增强器320的测量信号是信号SRED(见图30),它也在图11中展示为具有ILENGTH个采样值。信号SRED在图11的说明中也被称为I和2060。增强器信号处理涉及对离散输入信号SRED的离散自相关。输出信号O也称为SMDP,被展示在图12和图13中。
将要输入到增强器的测量信号SRED1、SRED可以包括至少一种振动信号分量SD,取决于所述可旋转运动部件的振动;其中所述振动信号分量具有重复频率fD,它取决于所述第一部件的转速fROT。信号分量SD的重复频率fD可以与被监视旋转部件的转速fROT成正比。
两个不同的故障特征SD1、SD2可以具有不同的频率fd1、fd2并且仍然可以由增强器增强,即改进SNR。因此,作为有益的情况,增强器320适于增强具有互相不同的重复频率fd1和fd2的不同特征SD1、SD2。两个重复频率fd1和fd2都与被监视旋转部件的转速fROT成正比,而fd1与fd2不同(fd1<>fd2)。这可以按照以下方式用数学表达:
fD1=k1*fROT,以及
fD2=k2*fROT,其中
k1和k2是正实数值,以及
k1<>k2,以及
k1大于或等于一(1),以及
k2大于或等于一(1)
增强器将输出信号序列传递给时域分析器290T的输入,以便当用户经由用户接口102、104选择执行时域分析时,时域分析器290T、105(图30和图4)将执行所选定的功能105,并且将相关数据传递到显示器106。增强器320的优点在于它在时域传递输出信号。因此,要求在时域输入信号的状态监视功能105、290T可以被设置为对图12和图13展示的信号输出的信号值进行直接操作。
当用户经由用户接口102、104选择执行频率域的分析时,增强器将把输出信号序列传递给快速傅氏变换器294,并且该快速傅氏变换器将把频率域数据的结果传递给频率域分析器290F、105(图30和图4)。频率域分析器290F、105将执行所选定的功能105并且将相关数据传递给显示器106。
在图29和图30所示的实施例中,作为有益的情况,用户很容易执行采用增强器和分数降采样器的分析。如图30所示,用户接口102、104、24B与参数控制器930协作,适于提供组合体920的若干设置。图31是参数控制器930的示意展示。
以下是参数设置的实例:
为了在频率域执行分析,经由用户接口102、104以下数据:
1)表明所关注的最高重复频率fD的信息。重复频率fD是所关注的特征SD的重复频率。以频率的形式或以表明所关注故障特征SD的最高重复频率的次序号OvHigh的形式可以输入这种信息。
2)表明对重复信号特征SD的SNR值期望改进的信息。以SNR改进器值L的形式可以输入这种信息。连同以上的图10A下面还讨论了SNR改进器值L。
3)表明在期望对增强器输出的信号执行FFT时,FFT 294中期望频率分辨率的信息。这种信息可以被设置为频率片段的值Z。根据本发明的实施例,通过从一组值中选择一个值Z,可设置频率分辨率Z。对频率分辨率Z的可选择值组可以包括
Z=400
Z=800
Z=1600
Z=3200
Z=6400
因此,尽管信号处理相当复杂,但是按照用户所要求的信息,组合体920已经被设计为提供优选情况下简单的用户接口。当用户对以上三个参数输入或选择数值时,全部其他数值都被自动设置或在组合体920中预置。
SNR改进器值L
将要输入的增强器的信号可以包括取决于可旋转运动部件振动的振动信号分量;其中所述振动信号分量具有取决于所述第一部件的转速fROT的重复频率fD;所述测量信号包括噪声也包括所述振动信号分量,所以所述测量信号具有关于所述振动信号分量的第一信噪比。增强器产生具有所述至少一个振动信号分量对应的重复信号分量的输出信号序列(O),所以所述输出信号序列(O)具有关于所述振动信号分量的第二信噪比值。本发明人已经通过测量结果确定,当SNR改进器值L被设置为数值一(1)时,第二信噪比值显著地高于第一信噪比。
不仅如此,本发明人还通过测量结果确定,当SNR改进器值L增大到L=4时,那么输出信号中关于所述振动信号分量的结果信噪比与L=1相关联的SNR值相比加倍。SNR改进器值L增大到L=10似乎致使对于输出信号中振动信号分量,与L=1时相同输入信号的SNR值相比,相关联的SNR值以因子3改进。因此,当SNR改进器值L从L1=1增大到L2时,结果SNR值可能以L2的平方根增大。
另外用户可以输入某设置,使组合体920保持重复测量。用户可以设置它以一定的重复周期TPM重复测量,即一直在时间TPM已经流逝时开始新的测量。可以被设置为一星期、或者一小时或十分钟。对这个重复频率选择的值取决于相关的测量条件。
由于增强器方法要求大量的数据输入值,即输入样点值的数目可能很大,并且它适于对缓慢旋转部件进行测量,所以测量的持续时间有时将相当长。因此,存在着测量的重复频率的用户设置与测量的持续时间不兼容的风险。所以,组合体920在收到以上用户输入后立即执行的步骤之一是计算测量的期望持续时间TM的估计。持续时间TM为:
TM=ILength/fSR2,
其中ILength是根据以下定义的选定的用户设置为了实现测量将要输入到增强器的信号中的样点数,而fSR2的定义如下。
组合体920也适于对比测量的持续时间TM和用户选定的重复周期值TPM。如果重复周期值TPM与测量的期望持续时间TM相比更短或大约相同,参数控制器930适于经由用户接口102、106提供警告指示,如通过显示器上的适宜文本。警告也可以包括声音或闪烁光。
根据某实施例,因为重复周期值TPM取决于算出的测量的持续时间TM的估计,所以组合体920适于计算某建议最小值。
根据以上用户设置,信号处理组合体920的参数控制器930能够设置信号处理功能94(图4)的全部参数,即整数降采样器设置和增强器设置。不仅如此,参数控制器930也能够在需要时设置分数降采样器的全部参数。在期望频率分析时,参数控制器930能够设置FFT294的参数。
以下参数可以在组合体920(图30)中预置:
AD转换器40、44的采样频率fS。
以下参数可以被测量:fROT
如上所述,参数值fROT可以被测量并与信号SRED1的对应样点值相关联地存储,SRED1的样点值被送入分数降采样器470B。
以下参数可以在组合体920中自动设置:
从增强器320输出的信号中的采样频率:
fSR2=C*Ov*fROT
其中
C是大于2.0的常数
Ov是用户输入的或者响应用户选定的要监视的最高频率而算出的次序号
fROT是在实际状态监视期间测量出的旋转部件的瞬时转速;
M=降采样器中使用的整数降采样器值选自包括总整数降采样所用的预定值组的表格。为了选择最适宜的值M,参数控制器930(图30)首先计算相当接近值M_calc=fS/fSR2*fROTmin/fROTmax
其中
fS和fSR2按以上定义,以及
fROTmin/fROTmax是表明测量期间所允许的最低与最高转速之间关系的值。根据值M_calc,选择器然后从预置值的列表中选择适宜的值。如通过从上述表格中选择小于M_calc的最接近值M,可以做到这一点。
fSR1=要从整数降采样器310传递的采样频率。fSR1被设置为fSR1=fS/M
D是分数降采样器的分数降采样器值。D可以被设置为D=fSR1/fSR2,其中fSR1和fSR2按以上定义。
OLENGTH=C*Z
其中
C是大于2.0的常数,比如上述2.56
Z是在期望对从增强器输出的信号进行FFT时,频率片段的选定数目,即表明FFT 294中期望频率分辨率的信息。
SSTART=OLENGTH或大于OLENGTH的值,其中OLENGTH按以上刚刚定义。
ILength=OLENGTH*L+SSTART+OLENGTH
CLength=ILENGTH-SSTART-OLENGTH
SMDP(t)=输出信号的样点值,在等式(5)中定义(见图10A)。
因此,参数控制器930适于产生以上定义的对应设置值,并且将它们传递到相关的信号处理功能94(图30和图4)。
一旦由增强器320已经产生了输出信号,便能够利用在控制输入端300上传递的选择信号,控制状态分析器290执行选定的状态分析功能105、290、290T、290F(图30)。在控制输入端300上传递的选择信号可以利用用户与用户接口102(见图2A和图30)的互动产生。当选定的分析功能包括快速傅氏变换时,将由选择信号300设置分析器290F在频率域中对输入信号进行操作。
快速傅氏变换器294可以适于对接收的具有一定数目的样点值的输入信号进行快速傅氏变换。当将样点值的一定数目设置为可以被二(2)除而不会出现分数的偶整数时是有益的。
根据本发明的优选实施例,从增强器输出的信号中的样点数目OLENGTH根据频率分辨率Z设置。频率分辨率Z与从增强器输出的信号中的样点数目OLENGTH之间的关系为:
OLENGTH=k*Z
其中
OLENGTH是从增强器320传递的信号中的样点值的样点数目。
K是大于2.0的因子
作为有益的情况,因子k可以选择为2.2与2.9之间的值,以便提供足够安全余量,同时避免不必要地产生许多样点值。
根据某实施例,作为有益的情况,选择的因子k使得100*k/2呈现为整数。这种选择致使OLENGTH的值适于作为快速傅氏变换器294的输入。根据某实施例,因子k可以被设置为2.56。选择k为2.56致使100*k=256=2的8次方。
表A表明了用户可选择的频率分辨率Z和对应的OLENGTH值的实例。
表A
本发明的实施例涉及分析机器状态的设备,所述机器具有可以某一旋转速度关于第二机器部件可旋转运动的第一部件,所述装置包括:
传感器,用于监视所述可运动部件,以便产生至少一个模拟测量信号,该模拟测量信号包括取决于所述可旋转运动部件的振动运动的至少一个振动信号分量,其中,所述振动信号分量具有取决于所述第一部件的旋转速度(fROT)的重复频率(fD);所述测量信号包括噪声也包括所述振动信号分量,所以所述测量信号具有关于所述振动信号分量的第一信噪比值;
A/D转换器(40、44),用于响应所述测量信号产生数字测量数据序列(SMD);所述数字测量数据序列(SMD)具有第一采样频率(fS);
第一数字滤波器(240),用于进行所述数字测量数据序列(SMD)的数字滤波,以便获得滤波后的测量信号(SF);
包络器,用于响应所述滤波后的测量信号(SF)产生第一数字信号(SENV、SMDP);
降采样器,用于进行所述第一数字信号(SENV、SMDP)的降采样,以实现具有降低采样频率(fSR1、fSR2)的降采样数字信号(SRED);
所述降采样器(470、470A、470B)具有
第一输入端,用于接收所述第一数字信号(SENV、SMDP);以及
第二输入端,用于接收指示所述可变旋转速度(fROT)的信号;
第三输入端,用于接收指示输出采样频率设置信号的信号;
所述降采样器(470、470A、470B)适于根据以下诸项产生所述降采样数字信号(SRED):
所述第一数字信号(SMD、SENV),
所述指示所述旋转速度(fROT)的信号,以及
所述指示输出采样频率设置信号的信号;其中,所述降采样器(470、470A、470B)适于产生所述降采样数字信号(SRED),使得所述旋转速度变化时所述旋转部件每转的样点值数目保持在基本上恒定值;以及
增强器(320),具有接收所述降采样数字信号(SRED)的输入;所述增强器适于产生具有所述至少一个振动信号分量对应的重复信号分量的输出信号序列(O),使得所述输出信号序列(O)具有关于所述振动信号分量的第二信噪比值;所述第二信噪比值高于所述第一信噪比值;以及
分析器(105;290;290T;294、290F),用于响应所述输出信号序列(O)指示取决于所述可旋转运动部件的所述振动运动的机器状态。这种解决方案通过使滤波器的复杂度最小化同时实现显著的性能改进,有利地提供了非常节约的解决方案。
根据某实施例,所述第一数字滤波器(240)是具有某截止频率的高通滤波器。
根据某实施例,所述冲击脉冲测量传感器(10)在某机械谐振频率值(fRMU)谐振;
所述高通滤波器240的截止频率fLC被选定为低于所述谐振冲击脉冲测量传感器(10)的机械谐振频率值(fRMU)的值。
根据另一个实施例,所述包络器(250)包括数字整流器(270),并且其中
所述降采样器(310)被连接以从所述整流器(270)的输出接收将要被降采样的信号;并且其中,
所述降采样器(310)包括低通滤波。在所述滤波器级中和所述包络器中去除低通滤波器组件使设计简化。不仅如此,由于使用了谐振传感器,也不再需要调整滤波器。
根据某实施例,降采样器包括整数降采样器(310),适于按整数M执行降采样;其中
根据所述旋转速度(fROT)自动设置数值M;
根据某实施例,降采样器包括多个低通有限脉冲响应滤波器。
根据某实施例,每个FIR滤波器都可设置为期望的降采样程度。
根据某实施例,整数降采样器包括一族不同的FIR滤波器,使得所述滤波器族中的多个滤波器可以组合以实现期望的总降采样M。
根据某实施例,所述滤波器族包括不止三个滤波器。
根据某实施例,所述整数降采样器(310)的输出端(312)被连接到分数降采样器(470)的输入端和选择器(460)的端子;所述选择器具有另一个端子;所述另一个端子被连接到所述分数降采样器(470)的输出端;以及
所述选择器(460)的第三端子被连接到所述增强器(320)的输入端。
根据某实施例,存在着使计算机分析具有缓慢旋转的部件的机器的状态的计算机程序,所述计算机程序包括:
计算机可读代码装置,在计算机上运行时,使所述计算机响应接收的模拟测量数据,以采样频率(fS)对模拟测量信号进行采样,以便产生数字测量数据序列(SMD);
计算机可读代码装置,在计算机上运行时,使所述计算机执行状态分析功能(F1、F2、Fn)以分析取决于所述数字包络信号(SENV、SMDP)的机器的状态。
根据某实施例,所述计算机程序产品包括
计算机可读介质;以及
如以上定义的计算机程序,
所述计算机程序记录在所述计算机可读介质上。
本发明的一方面A1涉及分析机器状态的设备,所述机器具有可以某一旋转速度关于第二机器部件可旋转运动的第一部件,所述装置包括:
传感器,用于监视所述可运动部件,以便产生至少一个模拟测量信号,该模拟测量信号包括取决于所述可旋转运动部件的振动运动的至少一个振动信号分量,其中,所述振动信号分量具有取决于所述第一部件的旋转速度(fROT)的重复频率(fD);所述测量信号包括噪声也包括所述振动信号分量,所以所述测量信号具有关于所述振动信号分量的第一信噪比值;以及
A/D转换器(40、44),用于响应所述测量信号产生数字测量数据序列(SMD);所述数字测量数据序列(SMD)具有第一采样频率(fS)。
所述装置还包括增强器(320),具有接收所述数字信号(SRED)的输入端。所述增强器适于产生具有所述至少一个振动信号分量对应的重复信号分量的输出信号序列O,使得所述输出信号序列O具有关于所述振动信号分量的第二信噪比值;所述第二信噪比值高于所述第一信噪比值。所述装置还包括分析器105;290;290T;294、290F,用于响应所述输出信号序列(O)指示取决于所述可旋转运动部件的所述振动运动的机器状态。
根据方面A1的装置的实施例A2,进一步包括参数控制器930;所述参数控制器适于与用户接口102、104、106协作,用于设置所述装置的控制参数,以便适应功能94、310、470、320、294,取决于用户经由所述接口的输入;其中所述装置适于请求经由所述用户接口102、104、106设置信息。
根据所述装置的实施例A2的装置的实施例A3,其中
所述请求的信息包括指示要被监视的频率范围的信息,用于具有所述重复频率fD的所述特征(SD)的检测。
根据实施例A2的装置的实施例A4,其中
所述请求的信息包括指示所关注的最高重复频率(fD)的信息。
根据实施例A2的装置的实施例A5,其中
所述请求的信息包括指示所述第二信噪比值的期望改进程度的信息。
根据实施例A2的装置的实施例A6,其中
所述请求的信息包括指示快速傅里叶变换的期望频率分辨率的信息。
根据实施例A2的装置的实施例A7,其中
所述请求的信息包括指示设置的信息,以使所述组合体920以一定的重复周期TPM保持重复测量。
根据实施例A2、A3、A4、A5、A6或A7的装置的实施例A8,其中
所述装置适于在用户输入的重复周期值TPM短于预计的测量持续时间TM或与其大致相同时经由用户接口102、106提供警告指示。
本发明的实施例A9涉及分析具有选择轴的机器状态的方法,所述方法包括:
接收具有第一多个(ILENGTH)样点值的数字输入信号(I),所述数字输入信号(I)表示旋转所述轴甚至使具有重复周期(TR)的振动发生而发出的机械振动;
在所述增强器中执行相关,以便产生输出信号序列O,其中关于随机信号分量,重复信号振幅分量被放大;
产生具有第二多个OLENGTH样点的数字输出信号O,所述第二多个OLENGTH为正整数并且小于所述第一多个ILENGTH。
根据实施例A9的方法的实施例A10,其中
所述数字输入信号(I)包括从退化的状态发出的信号分量,所述退化信号分量以所述重复周期(TR)重复。
根据实施例A9或A10的方法的实施例A11,其中
所述退化信号分量由噪声隐藏在所述数字输入信号I中。
本发明的一方面B1涉及分析具有选择轴的机器状态的方法,包括:
接收具有第一多个ILENGTH样点值的数字输入信号I,所述数字输入信号I表示旋转所述轴甚至使具有重复周期TR的振动发生而发源的机械振动,;
把所接收的数字输入信号I划分为第一信号部分2070和第二信号部分;
产生具有第二多个OLENGTH样点的数字输出信号O,所述第二多个OLENGTH为正整数并且小于所述第一多个ILENGTH;响应所述第一信号部分2070和所述第二信号部分而产生所述数字输出信号O。
根据方面B1的方法的实施例B2,其中
所接收的数字输入信号I包括从退化的状态发出信号分量,所述退化信号分量以所述重复周期TR重复。
根据实施例B2的方法的实施例B3,其中
所述退化信号分量由噪声隐藏在所述数字输入信号I中。
根据实施例B1、B2或B3的方法的实施例B4,其中
以重复信号振幅分量被放大的这样的方式产生所述数字输出信号O。
根据实施例B1、B2、B3或B4的方法的实施例B5,其中
所述轴的旋转量的增加引起所述重复信号振幅分量的重复发生的次数增加;并且所述方法进一步包括
以下列方式产生所述数字输出信号O
所述第一多个ILENGTH样点值的增加实现所述重复信号振幅分量的更高程度的放大。
根据实施例B1、B2、B3、B4或B5的方法的实施例B6,进一步包括
接收指示要被检测的振动的最低重复速率的值;
响应指示最低重复速率的值,确定第二多个OLENGTH样点的的最小值。
根据实施例B6的方法的实施例B7,进一步包括
计算引起所述轴进行所述最小旋转量的估计的测量持续时间,以便能够检测具有所述最低重复速率的任何振动。
根据实施例B1、B2、B3、B4、B5、B6或B7的方法的实施例B8,进一步包括
让用户能够选择退化信号分量的期望放大量;
响应所述期望放大量,计算第一多个ILENGTH数字输入信号样点值的值。
本发明的一方面B9涉及分析具有选择轴的机器状态的方法,包括:
接收具有第一多个ILENGTH样点值的数字输入信号I,所述数字输入信号I表示某持续时间从旋转所述轴发源的机械振动;
把所接收的数字输入信号I划分为第一信号部分2070和第二信号部分;所述第一部分表示所接收的数字输入信号I的第一持续时间,而所述第二信号部分表示所接收的数字输入信号I的第二持续时间,所述第一持续时间不同于所述第二持续时间;
产生具有第二多个OLENGTH样点的数字输出信号O,所述第二多个OLENGTH为正整数并且小于所述第一多个ILENGTH;其中
在时域中执行所述数字输出信号O中样点值O(t)的产生。
本发明的一方面B10涉及分析具有选择轴的机器状态的方法,包括:
接收具有第一多个ILENGTH样点值的数字输入信号I,所述数字输入信号I表示某持续时间从旋转所述轴发源的机械振动;
把所接收的数字输入信号I划分为第一信号部分2070和第二信号部分;
把所述第一信号部分2070与所述第二信号部分执行相关,以便产生输出信号序列O,其中关于随机信号分量,重复信号振幅分量被放大。
本发明的一方面B11涉及分析具有可以某一旋转速度旋转的部件的机器状态的设备,包括:
第一传感器,适于根据从所述部件的旋转发源的机械振动产生模拟电测量信号SEA;
模拟-数字转换器44,响应所接收的模拟测量数据,以采样频率fS对所述模拟测量信号进行采样,以便产生数字测量数据信号SMD;
第一降采样器,用于进行所述数字测量数据信号SMD、SENV的降采样,以实现具有第一降低采样频率fSR1的第一数字信号SMD、SENV;
增强器,具有接收所述第一数字信号SMD、SENV的输入;所述增强器适于接收第一多个ILENGTH样点值,其中,所述第一数字信号SMD、SENV表示某持续时间从所述部件的旋转发源的机械振动;
所述增强器适于把所接收的第一数字信号SMD、SENV划分为第一信号部分2070和第二信号部分;
所述增强器适于把所述第一信号部分2070与所述第二信号部分执行相关,以便产生输出信号序列O,其中关于随机信号分量,重复信号振幅分量被放大;以及
鉴别器230,用于执行状态分析功能F1、F2、Fn以根据所述第二数字信号SRED2分析所述机器的状态。
Claims (43)
1.一种分析具有可以某一旋转速度(fROT)旋转的部件的机器状态的设备,包括:
输入端(42),用于接收模拟测量信号(SEA),该信号指示具有振动频率(fSEA)和重复频率(fD)的振动信号特征;
A/D转换器(40、44),用于根据所述模拟测量信号产生数字测量信号(SMD),所述数字测量信号(SMD)具有第一采样频率(fS),所述第一采样频率为所述振动频率(fSEA)的至少两倍(k);
包络器(250、94),用于产生指示所述重复频率(fD)的包络信号;
降采样器(310、94),用于根据所述包络信号产生降采样后第一数字信号(SRED、SRED1、SRED2),使得所述降采样后第一数字信号具有降低采样频率(fSR);
增强器(320),具有用于接收所述降采样信号的输入端,其中所述增强器(320)在时域中运行,对所述降采样后第一数字信号(SRED、SRED1、SRED2)进行离散自相关,以便产生增强器输出信号序列(O、SMDP);以及
分析器(290;290T、294、290F),用于执行状态监视功能(F1、F2、Fn),以便根据所述增强器输出信号序列(O、SMDP)分析机器的状态。
2.根据权利要求1的设备,其中,
所述增强器(320)用于对所述降采样后第一数字信号(SRED、SRED1、SRED2)以延迟(SSTART、NC)进行离散自相关。
3.根据权利要求1或2的设备,其中,
所述增强器适于处理所述第一数字信号(SRED、SRED1、SRED2)的第一多个(ILENGTH)输入样点值,以便产生具有第二多个(OLENGTH)输出样点值(SMDP(t))的输出信号序列(O、SMDP);所述第二多个为正整数。
4.根据权利要求2的设备,其中,
所述增强器适于处理所述第一数字信号(SRED、SRED1、SRED2)的第一多个(ILENGTH)输入样点值,以便产生具有第二多个(OLENGTH)输出样点值(SMDP(t))的输出信号序列(O、SMDP);所述第二多个为正整数;以及
所述延迟(SSTART、NC)为整数,具有的值等于或大于所述第二多个值(OLENGTH)。
5.根据任何一个前面的权利要求的设备,其中,
所述增强器(320、382、386)适于根据第三多个(CLENGTH)输入信号乘积(P(i,t))产生输出样点值(SMDP(t));所述第三多个(CLENGTH)为正整数;所述增强器(320、382、386)进一步包括
根据以下信息产生指示所述第三多个(CLENGTH)的数据值的装置(365、375;375B):
信噪比改进器值(L),以及
指示所述第二多个(OLENGTH)的数据(Z、Y、fSR、fROT)。
6.一种分析具有可以某一旋转速度(fROT)旋转的部件的机器状态的设备,包括:
第一输入端(315),用于根据从所述部件的旋转发源的机械振动接收第一数字信号(SRED、SMD、SENV),所述第一数字信号(SRED)具有采样频率(fSR);其中,所述第一数字信号(SRED、SMD、SENV)可以包括噪声和具有第一重复频率(fD1)的第一信号特征(SD1);
增强器(320),用于处理所述第一数字信号(SRED、SMD、SENV)的第一多个(ILENGTH)输入样点值,以便产生具有第二多个(OLENGTH)输出样点值(SMDP(t))的输出信号序列(O、SMDP);所述第二多个为正整数;
分析器(290;290T、294、290F),用于执行状态监视功能(F1、F2、Fn),以便根据所述输出信号序列(O、SMDP)分析机器的状态;其中,增强器(320)包括:
至少一个第二输入端(335、350、360、364),用于接收指示所述第二多个(OLENGTH)的数据(Z、Y、fSR、fROT、OLENGTH);以及
传递所述输出信号序列(O、SMDP)的输出端(394、396、398);其中
所述增强器(320、382、386)适于根据第三多个(CLENGTH)输入信号乘积(P(i,t))产生输出样点值(SMDP(t));所述第三多个(CLENGTH)为正整数;其中,通过将第一样点位置(i)处的第一输入样点值(I(i))乘以第二样点位置(i+t+SSTART)处的第二输入样点值(I(i+t+SSTART)),获得输出样点位置t的输入信号乘积(P(i,t)),所述第二输入样点值与所述第一输入样点值分开一定数目(t+SSTART)的样点位置,其中所述一定数目等于或大于一定值(SSTART);
所述增强器进一步包括:
至少一个第三输入端(335、366),用于接收指示信噪比改进器值(L)的数据;以及
根据以下信息产生指示所述第三多个(CLENGTH)的数据值的装置(365、375;375B):
所述信噪比改进器值(L),以及
指示所述第二多个(OLENGTH)的所述数据(Z、Y、fSR、fROT)。
7.根据权利要求6的设备,其中,所述增强器(320、382、386)适于通过增加(S1140)所述第三多个(CLENGTH)输入信号乘积(P(i,t)),为所述输出样点位置t产生所述输出样点值(SMDP(t))。
8.一种分析具有可以某一旋转速度(fROT)旋转的部件的机器状态的设备,包括:
第一输入端(315),用于根据从所述部件的旋转发源的机械振动接收第一数字信号(SRED、SMD、SENV),所述第一数字信号(SRED)具有采样频率(fS);其中,所述第一数字信号(SRED、SMD、SENV)可以包括噪声和具有第一重复频率(fD1)的第一信号特征(SD1);
增强器(320),用于处理所述第一数字信号(SRED、SMD、SENV)的第一多个(ILENGTH)输入样点值,以便产生具有第二多个(OLENGTH)输出样点值(SMDP(t))的输出信号序列(O、SMDP);所述第二多个为正整数;
分析器(290;290T、294、290F),用于执行状态监视功能(F1、F2、Fn),以便根据所述输出信号序列(O、SMDP)分析机器的状态;其中,增强器(320)包括:
至少一个第二输入端(335、350、360、364),用于接收指示所述第二多个(OLENGTH)的数据(Z、Y、fSR、fROT、OLENGTH);以及
传递所述输出信号序列(O、SMDP)的输出端(394、396、398);其中
所述增强器(320、382、386)适于根据表达式
或者根据表达式
产生输出样点值(SMDP(t));其中,
SSTART为整数,具有的值等于或大于所述第二多个值(OLENGTH);
CLENGTH为整数;以及
i0为整数常量;
并且其中所述增强器进一步包括:
至少一个第三输入端(335、366),用于接收指示信噪比改进器值(L)的数据;以及
根据以下信息产生指示所述整数值CLENGTH的数据值的装置(365、375;375B):
所述信噪比改进器值(L),以及
指示所述第二多个(OLENGTH)的所述数据(Z、Y、fSR、fROT)。
9.根据权利要求5、6、7或8中任何一个的设备,其中,
产生指示所述第三多个(CLENGTH)的数据值的所述装置(375、375B)适于接收
指示所述第一多个(ILENGTH)的数据,以及
指示所述第二多个(OLENGTH)的数据,以及
指示所述一定值(SSTART)的数据;以及其中
产生指示所述第三多个(CLENGTH)的数据值的所述装置(375、375B)适于设置第三多个的值CLENGTH等于值ILENGTH与值SSTART和OLENGTH之和之间的差。
10.根据权利要求5、6、7或8中任何一个的设备,其中,
产生指示所述第三多个(CLENGTH)的数据值的所述装置(375、375B)适于接收
指示所述第二多个(OLENGTH)的数据,以及
指示所述信噪比改进器值(L)的数据;以及其中
产生指示所述第三多个(CLENGTH)的数据值的所述装置(375、375B)适于设置第三多个的值CLENGTH等于第二多个的值(OLENGTH)与所述信噪比改进器值(L)的乘积。
11.根据任何一个前面的权利要求的设备,其中所述增强器(320)进一步包括:
随机信号频率衰减器设置运算器(370),适于接收指示所述第二多个(OLENGTH)的数据;其中
所述随机信号频率衰减器设置运算器(370)适于设置指示所述一定值(SSTART)的整数数据值SSTART等于第二多个的值(OLENGTH)或大于第二多个的值(OLENGTH)的值。
12.根据任何一个前面的权利要求的设备,其中所述增强器(320)进一步包括:
输入信号长度计算器(365),被安排为设置第一多个(ILENGTH)为等于整数IL的值或大于所述整数IL的值,其中IL=OLENGTH*L+SSTART+OLENGTH
其中,OLENGTH为等于所述第二多个的值的整数;
L为所述信噪比改进器值;以及
SSTART为所述一定值。
13.根据任何一个前面的权利要求的设备,其中所述增强器(320)进一步包括:
输出信号长度计算器(345、340、345B),被安排为根据借助输出信号序列(O、SMDP)的分析可得到的指示频率分辨率(Z)的数据,设置第二多个的值(OLENGTH)等于某整数值。
14.根据权利要求13的设备,其中
所述输出信号长度计算器(345、340、345B)适于设置第二多个的值(OLENGTH)等于指示所述频率分辨率(Z)的值与常数C的乘积,其中所述常数C具有等于2的值或大于2.0的值。
15.根据权利要求13的设备,其中
所述输出信号长度计算器(345、340、345B)适于根据以下数据设置所述第二多个的值(OLENGTH):
借助输出信号序列(O、SMDP)的分析可得到的指示频率分辨率(Z)的数据;以及
指示所述旋转速度(fROT)的数据;以及
指示所述采样频率(fSR)的数据;以及
指示级别值(Y、OV、OvHigh)的数据,所述级别值为在所述输出信号序列(O、SMDP)中可被检测的最高重复频率(fDmax)与所述旋转速度(fROT)之间的关系。
16.根据权利要求13的设备,其中
所述输出信号长度计算器(345、340、345B)适于根据以下数据设置第二多个的值(OLENGTH):
所述借助输出信号序列(O、SMDP)的分析可得到的指示频率分辨率(Z)的数据;以及
指示所述采样频率(fSR)的数据;以及
指示在所述输出信号序列(O、SMDP)中可被检测的最高重复频率(fDmax)的数据。
17.根据权利要求1-12中任何一个的设备,其中,
所述输出信号长度计算器(345、340、345B)适于设置第二多个的值(OLENGTH)为预定值。
18.根据权利要求1-12中任何一个的设备,进一步包括
用户接口(102、104、24B),适于允许用户输入指示所述第二多个(OLENGTH)的数据(Z、Y、fSR、fROT、OLENGTH);其中
所述用户接口适于向所述至少一个第二输入端(335、350、360、364)传递指示所述第二多个(OLENGTH)的数据(Z、Y、fSR、fROT、OLENGTH)。
19.根据权利要求18的设备,其中,
所述用户接口(102、104、24B)适于允许用户设置第二多个的值(OLENGTH)。
20.根据任何一个前面的权利要求的设备,进一步包括
控制所述采样频率(fSR、fSR1、fSR2)的装置,使得它等于
fSR=C*Y*fROT
其中
C为具有等于2的值或大于2.0的值的常数;以及
Y为级别值,即在所述输出信号序列(O、SMDP)中可被检测的最高重复频率(fDmax)与所述旋转速度(fROT)之间的关系;以及
fROT是所述可旋转部件的所述旋转速度。
21.根据任何一个前面的权利要求的设备,进一步包括
检测器(420),用于产生指示所述可旋转部件的所述旋转速度(fROT)的值;其中
所述增强器(320)具有输入端(360),用于接收所述旋转速度值。
22.根据权利要求1-20中任何一个的设备,其中,
所述用户接口(102、104、24B)适于允许用户设置旋转速度值(fROT)为可选择的值;其中
所述增强器(320)具有输入端(360),用于接收所述旋转速度值。
23.根据任何一个前面的权利要求的设备,其中,
所述第一数字信号(SRED、SMD、SENV)可以包括第二信号特征(SD2);以及其中
所述增强器(320)适于处理所述第一多个(ILENGTH)输入样点值,使得关于所述噪声,第一重复频率(fD1)和第二重复频率(fD2)都被增强,即使在所述第二重复频率(fD2)不同于所述第一重复频率(fD1)时。
24.根据任何一个前面的权利要求的设备,其中,
所述第一数字信号(SRED、SMD、SENV)可以包括具有第二重复频率(fD2)的第二信号特征(SD2);以及其中
所述增强器(320)适于处理所述第一多个(ILENGTH)输入样点值,使得关于所述噪声,第一重复频率(fD1)和第二重复频率(fD2)都被增强,即使在所述第二重复频率(fD2)不同于所述第一重复频率(fD1),以及第一重复频率(fD1)和第二重复频率(fD2)都与可旋转运动部件(8)的所述旋转速度(fROT)成正比时。
25.根据任何一个前面的权利要求的设备,进一步包括
传感器(10),用于根据从所述部件的旋转发源的所述机械振动产生测量信号;以及
装置(40、44、94、250、310、470),用于根据所述测量信号产生具有所述采样频率(fSR)的所述第一数字信号(SRED、SMD、SENV)。
26.根据任何一个前面的权利要求的设备,其中,
所述传感器(10)是冲击脉冲传感器。
27.根据任何一个前面的权利要求的设备,其中,
所述第一数字信号的所述输入样点值被存储为输入样点值的向量,所述向量具有下标值用于寻址各个输入样点值;所述下标值的范围从最低的下标值(i0、iSTART)到最高的下标值(i0+ILENGTH-1;iSTART+ILENGTH-1);其中
输入样点值的范围和下标值被选择为所述下标值是输入信号向量的范围内的值。
28.根据任何一个前面的权利要求的设备,其中,
输入信号乘积(P(i,t))中两项的两个下标值之间的差(NC)等于或大于所述一定值(SSTART)。
29.一种用于分析具有可以某一旋转速度(fROT)旋转的部件的机器状态的方法,包括以下步骤:
根据从所述部件的旋转发源的机械振动接收第一数字信号(SRED、SMD、SENV),所述第一数字信号(SRED)具有采样频率(fS);其中,所述第一数字信号(SRED、SMD、SENV)可以包括噪声和具有第一重复频率(fD1)的第一信号特征(SD1)及具有第二重复频率(fD2)的第二信号特征(SD2),进一步包括:
处理所述第一数字信号(SRED、SMD、SENV)的第一多个(ILENGTH)输入样点值,以便产生具有第二多个(OLENGTH)输出样点值(SMDP(t))的输出信号序列(O、SMDP),其中
各个输入样点值(I(i))可由输入样点位置值(i)寻址,以及
各个输出样点值(SMDP(t)、O(t))可由输出样点位置值(t)寻址;以及其中
所述处理步骤包括:
A)对于输出样点位置t产生输出样点值SMDP(t),所述输出样点值产生步骤包括:
执行第三多个(CLENGTH)乘法以产生输入信号乘积;通过将第一样点位置(i)处的第一输入样点值(I(i))乘以第二样点位置(i2)处的第二输入样点值(I(i2)、I(i+SSTART+t+t0)、I(i-SSTART+t+t0)),获得输入信号乘积,所述第二输入样点值(I(i+SSTART+t0)、I(i-SSTART+t0))与所述第一输入样点值(I(i))分开一定数目(SSTART+t)的样点位置,所述一定数目等于或大于一定值(SSTART);以及
增加所述第三多个(CLENGTH)输入信号乘积,以产生输入信号乘积之和;所述输入信号乘积之和为对于所述样点位置t的所述输出样点值(SMDP(t));以及
传递对于所述样点位置t的所述输出样点值(SMDP(t));
B)对于基本上全部所述输出样点位置值(t)重复步骤A),使得产生所述输出信号序列(O、SMDP);所述方法进一步包括:
执行状态监视功能(F1、F2、Fn),以便根据所述输出信号序列(O、SMDP)分析机器的状态。
30.一种运行分析具有可以某一旋转速度(fROT)旋转的部件的机器状态的设备的方法,包括以下步骤:
根据从所述部件的旋转发源的机械振动接收第一数字信号(SRED、SMD、SENV),所述第一数字信号(SRED)具有采样频率(fSR);
处理所述第一数字信号(SRED、SMD、SENV)的第一多个(ILENGTH)输入样点值,以便产生具有第二多个(OLENGTH)输出样点值(SMDP(t))的输出信号序列(O、SMDP),其中
各个输入样点值(I(i))可由输入样点位置值(i)寻址,以及
各个输出样点值(SMDP(t)、O(t))可由输出样点位置值(t)寻址;以及其中所述处理步骤包括:
1)设置所述输出样点位置值(t)为初始值(t0);
2)对于相关样点位置(t)产生输出样点值SMDP(t),所述输出样点值产生步骤包括:
a)设置所述输入样点位置值(i)为初始值(i0);
b)通过将第一样点位置(i)处的第一输入样点值(I(i))乘以第二样点位置(i2)处的第二输入样点值(I(i2)、I(i+SSTART+t0)、I(i-SSTART+t0)),产生输入信号乘积(P(i,t)),所述第二输入样点值(I(i+SSTART+t0)、I(i-SSTART+t0))与所述第一输入样点值(I(i))分开一定数目(t+SSTART)的样点位置,所述一定数目大于一定值(SSTART);以及
c)通过一个步骤改变输入样点位置值(i);
d)重复步骤b),以便产生另一个输入信号乘积(P(i,t))并且产生所产生的输入信号乘积之和;其中
e)重复步骤d)直到第三多个(CLENGTH)输入信号乘积(P(i,t))已经求和;
f)根据所产生的输入信号乘积(P(i,t))之和,设置与由输出样点位置值(t)指示的输出样点位置值(t)相关的输出样点值(SMDP(t));
g)传递与由输出样点位置值(t)指示的输出样点位置值(t)相关的输出样点值(SMDP(t));
3)通过一个步骤改变输出样点位置值(t);
4)重复步骤2),直到已经传递了至少所述第二多个(OLENGTH)输出样点值(SMDP(t));所述方法进一步包括:
执行状态监视功能(F1、F2、Fn),以便根据所述输出信号序列(O、SMDP)分析机器的状态。
31.根据权利要求29或30的方法,其中,
所述第二多个(OLENGTH)输出样点值可设置。
32.根据权利要求29、30或31的方法,包括以下步骤:
允许用户输入指示所述第二多个的值OLENGTH的数据。
33.根据权利要求29、30或31的方法,包括以下步骤:
所述数字输入信号(I)表示旋转所述轴甚至引起具有重复周期(TR、1/fD、TIR;TOR)的振动信号分量(SD)的重复发生而发出的机械振动;以及其中
所述第二多个(OLENGTH)的值OLENGTH被设置为大到允许所述输出信号序列(O、SMDP)表示振动信号分量(SD)的至少一个周期(TIR;TOR)的值。
34.根据权利要求29-33中任何一个的方法,包括以下步骤:
允许用户输入指示要被检测的最低重复频率(fREPmin)的数据;以及
根据所述指示最低重复频率(fREPmin)的数据,自动设置第二多个(OLENGTH)的值OLENGTH。
35.根据任何一个前面的权利要求的方法,包括以下步骤:
允许用户输入指示可被检测的最高重复频率(Y、OV、OvHIGH)的数据;以及
根据所述指示可被检测的最高重复频率(Y、OV、OvHIGH)的数据,自动设置第二多个(OLENGTH)的值OLENGTH。
36.根据任何一个前面的权利要求的方法,包括以下步骤:
允许用户输入借助输出信号序列(O、SMDP)的分析可得到的指示期望频率分辨率(Z)的数据;以及
根据所述指示期望频率分辨率(Z)的数据,自动设置第二多个(OLENGTH)的值OLENGTH。
37.根据任何一个前面的权利要求的方法,包括以下步骤:
允许用户输入借助输出信号序列(O、SMDP)的分析可得到的指示期望频率分辨率(Z)的数据;以及
设置值OLENGTH为所述期望频率分辨率(Z)与具有大于2.0的值的常数(C、k)的乘积。
38.根据任何一个前面的权利要求的方法,包括以下步骤:
根据以下数据设置第二多个(OLENGTH)的值OLENGTH:
借助输出信号序列(O、SMDP)的分析可得到的指示频率分辨率(Z)的数据;以及
指示可检测的最高重复频率(Y、OV、OvHIGH)的数据;以及
指示所述旋转速度(fROT)的数据;以及
指示所述采样频率(fSR)的数据。
39.一种分析具有可以某一旋转速度(fROT)旋转的部件的机器状态的方法,包括以下步骤:
根据从所述部件的旋转发源的机械振动接收第一数字时域信号(SRED、SMD、SENV),所述第一数字信号(SRED)具有采样频率(fSR);其中,所述第一数字信号(SRED、SMD、SENV)可以包括噪声和具有第一重复频率(fD1)的第一信号特征(SD1),进一步包括:
处理所述第一数字信号(SRED、SMD、SENV)的第一多个P1(ILENGTH)输入样点值,以便产生具有第二多个P2(OLENGTH)输出样点值(SMDP(t))的输出信号序列(O、SMDP);所述输出信号序列(O、SMDP)为时域信号序列;其中
所述第一多个P1输入样点值的所述处理步骤包括输入样点值的第三多个P3(CLENGTH)乘法,用于产生一个输出样点值(SMDP(t)),使得产生所述第二多个P2(OLENGTH)输出样点值(SMDP(t))包括至少第四多个P4乘法(P4=P2*P3)。
40.根据任何一个前面的权利要求的方法,包括以下步骤:
接收指示信噪比改进器值(L)的数据;以及
装置(365、375;375B),根据以下信息产生指示所述第三多个(CLENGTH)的数据值:
所述信噪比改进器值(L),以及
指示所述第二多个(OLENGTH)的所述数据(Z、Y、fSR、fROT)。
41.一种运行分析具有可以某一旋转速度(fROT)旋转的部件的机器状态的设备的方法,包括以下步骤:
接收模拟测量信号(SEA),该信号指示具有振动频率(fSEA)和重复频率(fD)的振动信号特征;
根据模拟测量信号产生数字测量信号(SMD),所述数字测量信号(SMD)具有第一采样频率(fS),所述第一采样频率为所述振动频率(fSEA)的至少两倍(k);
产生指示所述重复频率(fD)的包络信号;
根据所述包络信号产生降采样后第一数字信号(SRED、SRED1、SRED2),使得所述降采样后第一数字信号具有降低采样频率(fSR);
在时域中对所述降采样后第一数字信号(SRED、SRED1、SRED2)进行离散自相关,以便产生增强的输出信号序列(O、SMDP);以及
执行状态监视功能(F1、F2、Fn),以便根据所述增强的输出信号序列(O、SMDP)分析机器的状态。
42.一种计算机程序,当所述程序在计算机上或分析设备上运行时,包括的计算机可读代码装置使计算机执行权利要求29-41中任何一个的所述步骤。
43.一种计算机程序产品,包括
计算机可读介质;以及
根据权利要求42的计算机程序,
所述计算机程序记录在所述计算机可读介质上。
Applications Claiming Priority (9)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US17550509P | 2009-05-05 | 2009-05-05 | |
US61/175,505 | 2009-05-05 | ||
SE0950312 | 2009-05-05 | ||
SE0950312-9 | 2009-05-05 | ||
SEPCT/SE2009/051494 | 2009-12-22 | ||
PCT/SE2009/051494 WO2010074648A1 (en) | 2008-12-22 | 2009-12-22 | An analysis system |
SE1000044 | 2010-01-18 | ||
SE1000044-6 | 2010-01-18 | ||
PCT/SE2010/050257 WO2010128928A1 (en) | 2009-05-05 | 2010-02-26 | An apparatus and a method for analysing the vibration of a machine having a rotating part |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102449445A true CN102449445A (zh) | 2012-05-09 |
CN102449445B CN102449445B (zh) | 2014-12-17 |
Family
ID=43050276
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201080019737.8A Active CN102449445B (zh) | 2009-05-05 | 2010-02-26 | 用于分析具有旋转部件的机器振动的设备和方法 |
Country Status (6)
Country | Link |
---|---|
US (3) | US9964430B2 (zh) |
EP (1) | EP2427745A4 (zh) |
CN (1) | CN102449445B (zh) |
AU (1) | AU2010245354A1 (zh) |
EA (1) | EA022630B1 (zh) |
WO (1) | WO2010128928A1 (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107643182A (zh) * | 2017-09-08 | 2018-01-30 | 安徽大学 | 基于欠采样信号的轴承故障诊断方法及无线网络节点装置 |
CN108318122A (zh) * | 2017-01-16 | 2018-07-24 | 上银科技股份有限公司 | 自适应振动感测方法 |
CN109051618A (zh) * | 2018-09-06 | 2018-12-21 | 广东寰球智能科技有限公司 | 基于振动信号处理传送星轮故障的监测方法及监测系统 |
CN110546477A (zh) * | 2018-03-29 | 2019-12-06 | 三菱重工业株式会社 | 接触振动检测装置及具备该接触振动检测装置的旋转机械、以及接触振动检测方法 |
CN113670428A (zh) * | 2021-07-19 | 2021-11-19 | 国网四川省电力公司电力科学研究院 | 一种变压器振动在线异常检测方法 |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9304033B2 (en) * | 2008-12-22 | 2016-04-05 | S.P.M. Instrument Ab | Analysis system |
CN102449445B (zh) * | 2009-05-05 | 2014-12-17 | S.P.M.仪器公司 | 用于分析具有旋转部件的机器振动的设备和方法 |
US9903090B2 (en) * | 2012-01-18 | 2018-02-27 | Harnischfeger Technologies, Inc. | System and method for vibration monitoring of a mining machine |
US10337957B2 (en) * | 2012-01-30 | 2019-07-02 | S.P.M. Instrument Ab | Apparatus and method for analysing the condition of a machine having a rotating part |
CN102759448B (zh) * | 2012-06-26 | 2014-11-05 | 西安瑞特快速制造工程研究有限公司 | 基于柔性时域平均的齿轮箱故障检测方法 |
WO2014042582A1 (en) * | 2012-09-11 | 2014-03-20 | S.P.M. Instrument Ab | Apparatus for monitoring the condition of a machine |
CN103063237B (zh) * | 2012-12-25 | 2015-05-20 | 北京华信创通科技有限公司 | 一种使编码器抗干扰的方法及装置 |
WO2014188072A1 (en) | 2013-05-21 | 2014-11-27 | Konecranes Plc | Techniques for monitoring gear condition |
JP6659384B2 (ja) * | 2016-02-02 | 2020-03-04 | 株式会社神戸製鋼所 | 回転機の異常検知装置および回転機の異常検知システム |
US20180284746A1 (en) | 2016-05-09 | 2018-10-04 | StrongForce IoT Portfolio 2016, LLC | Methods and systems for data collection optimization in an industrial internet of things environment |
US11327475B2 (en) | 2016-05-09 | 2022-05-10 | Strong Force Iot Portfolio 2016, Llc | Methods and systems for intelligent collection and analysis of vehicle data |
US11774944B2 (en) | 2016-05-09 | 2023-10-03 | Strong Force Iot Portfolio 2016, Llc | Methods and systems for the industrial internet of things |
US10914701B2 (en) | 2016-06-21 | 2021-02-09 | General Electric Company | Systems and methods for determining rotor deterioration in a dynamoelectric machine |
JP6665062B2 (ja) * | 2016-08-31 | 2020-03-13 | Ntn株式会社 | 状態監視装置 |
RU2646207C1 (ru) * | 2016-11-23 | 2018-03-01 | Общество с ограниченной ответственностью НПЦ "Динамика" - Научно-производственный центр "Диагностика", надежность машин и комплексная автоматизация" | Способ диагностического мониторинга роторных механизмов |
US10921801B2 (en) | 2017-08-02 | 2021-02-16 | Strong Force loT Portfolio 2016, LLC | Data collection systems and methods for updating sensed parameter groups based on pattern recognition |
RU2658233C1 (ru) * | 2017-09-12 | 2018-06-19 | Общество с ограниченной ответственностью НПЦ "Динамика" - Научно-производственный центр "Диагностика, надежность машин и комплексная автоматизация" | Способ диагностики технического состояния роторных механизмов |
US20200133254A1 (en) * | 2018-05-07 | 2020-04-30 | Strong Force Iot Portfolio 2016, Llc | Methods and systems for data collection, learning, and streaming of machine signals for part identification and operating characteristics determination using the industrial internet of things |
US11226315B2 (en) * | 2019-04-09 | 2022-01-18 | The Boeing Company | Passive noise-cancellation for ultrasonic evaluation |
CN111289097A (zh) * | 2020-02-25 | 2020-06-16 | 河海大学 | 一种基于振动信号的锂离子电池故障预警方法及系统 |
CN114563181B (zh) * | 2022-01-10 | 2023-06-27 | 浙江工业大学之江学院 | 基于改进变分模态提取的旋转机械故障特征提取方法 |
CN117740381A (zh) * | 2024-01-22 | 2024-03-22 | 中国矿业大学 | 一种低速重载工况下的轴承故障诊断方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS60195426A (ja) * | 1984-03-16 | 1985-10-03 | Komatsu Ltd | 回転機構の故障診断方法 |
GB2190198A (en) * | 1986-04-28 | 1987-11-11 | Vipac Pty Ltd | Vibrational analysis system for a machine |
JPH01178814A (ja) * | 1988-01-08 | 1989-07-17 | Mitsubishi Heavy Ind Ltd | 信号の大局的構造認識装置 |
JPH04279826A (ja) * | 1991-03-08 | 1992-10-05 | Nkk Corp | 可変速回転系異常診断方法及びその装置 |
WO1998036251A1 (en) * | 1997-02-18 | 1998-08-20 | Dme Corporation | Method and apparatus for predictive diagnosis of moving machine parts |
JPH1178814A (ja) * | 1997-09-11 | 1999-03-23 | Honda Motor Co Ltd | 車両用多板式ブレーキ |
WO2002073150A2 (en) * | 2001-03-13 | 2002-09-19 | Ab Skf | System and method for analyzing vibration signals |
US20040101048A1 (en) * | 2002-11-14 | 2004-05-27 | Paris Alan T | Signal processing of multi-channel data |
JP2006113003A (ja) * | 2004-10-18 | 2006-04-27 | Nsk Ltd | 機械設備の異常診断システム |
CN1926413A (zh) * | 2004-10-18 | 2007-03-07 | 日本精工株式会社 | 机械设备的异常诊断系统 |
Family Cites Families (50)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3554012A (en) | 1968-02-29 | 1971-01-12 | Iko Konsultbyron Ab | Method and arrangement for determining the mechanical state of machines |
US3705516A (en) | 1971-09-30 | 1972-12-12 | Northrop Corp | Method and apparatus for testing the condition of a machine |
US4528852A (en) | 1982-10-21 | 1985-07-16 | Spm Instruments U.S. Inc. | Method and instrument for determining the condition of an operating bearing |
DE3424692A1 (de) | 1984-07-05 | 1986-02-06 | Kletek Controllsysteme GmbH & Co KG, 2820 Bremen | Verfahren und messanordnung zur analyse von periodischen oder quasi periodischen signalen, insbesondere von schallsignalen bei maschinen und anlagen |
US5258923A (en) | 1987-07-22 | 1993-11-02 | General Electric Company | System and method for detecting the occurrence, location and depth of cracks in turbine-generator rotors |
JPH01127934A (ja) | 1987-11-12 | 1989-05-19 | Hitachi Ltd | 損傷診断装置 |
US4912661A (en) | 1987-12-01 | 1990-03-27 | Hewlett-Packard Company | Tracking and resampling method and apparatus for monitoring the performance of rotating machines |
GB2228088B (en) | 1988-12-16 | 1992-09-16 | Nippon Seiko Kk | Method and apparatus for detecting cracks in bearings |
US5109700A (en) | 1990-07-13 | 1992-05-05 | Life Systems, Inc. | Method and apparatus for analyzing rotating machines |
US5201292A (en) | 1991-08-30 | 1993-04-13 | Loral Aerospace Corp. | Apparatus and method for detecting vibration patterns |
US5501105A (en) | 1991-10-02 | 1996-03-26 | Monitoring Technology Corp. | Digital signal processing of encoder signals to detect resonances in rotating machines |
US5365787A (en) | 1991-10-02 | 1994-11-22 | Monitoring Technology Corp. | Noninvasive method and apparatus for determining resonance information for rotating machinery components and for anticipating component failure from changes therein |
US5445028A (en) | 1992-09-18 | 1995-08-29 | Ametek Aerospace Products Inc. | Dynamic digital tracking filter |
US5390545A (en) | 1993-01-11 | 1995-02-21 | Caterpillar Inc. | Apparatus for measuring torsional vibrations of rotating machinery |
WO1995013655A1 (en) | 1993-11-09 | 1995-05-18 | Motorola, Inc. | Method and apparatus for detecting an input signal level |
US5548540A (en) * | 1994-06-24 | 1996-08-20 | General Electric Company | Decimation filter having a selectable decimation ratio |
US5633811A (en) | 1994-12-09 | 1997-05-27 | Computational Systems, Inc. | Hand held data collector and analyzer system |
US5895857A (en) | 1995-11-08 | 1999-04-20 | Csi Technology, Inc. | Machine fault detection using vibration signal peak detector |
SE510771C2 (sv) | 1996-07-05 | 1999-06-21 | Spm Instr Ab | Förfarande för utvärdering av konditionen för en maskin jämte analysapparat samt anordning för samverkan med analysapparaten |
US6351714B1 (en) * | 1998-03-03 | 2002-02-26 | Entek Ird International Corporation | Order tracking signal sampling process |
US6041339A (en) * | 1998-03-27 | 2000-03-21 | Ess Technology, Inc. | Efficient decimation filtering |
FI112972B (fi) | 1998-07-15 | 2004-02-13 | Abb Research Ltd | Laakerin kunnon arviointi |
US6053047A (en) | 1998-09-29 | 2000-04-25 | Allen-Bradley Company, Llc | Determining faults in multiple bearings using one vibration sensor |
AUPQ152499A0 (en) | 1999-07-09 | 1999-08-05 | Commonwealth Scientific And Industrial Research Organisation | A system for monitoring acoustic emissions from a moving machine |
US6868348B1 (en) | 1999-10-29 | 2005-03-15 | Entek Ird International Corporation | Adaptive high frequency energy detection |
US6351713B1 (en) | 1999-12-15 | 2002-02-26 | Swantech, L.L.C. | Distributed stress wave analysis system |
US6332116B1 (en) | 2000-04-19 | 2001-12-18 | National Instruments Corporation | System and method for analyzing signals of rotating machines |
US6591682B1 (en) | 2000-08-14 | 2003-07-15 | Pruftechnik Dieter Busch Ag | Device and process for signal analysis |
US7136794B1 (en) | 2001-05-24 | 2006-11-14 | Simmonds Precision Products, Inc. | Method and apparatus for estimating values for condition indicators |
TW579424B (en) | 2001-07-09 | 2004-03-11 | Shell Int Research | Vibration analysis for predictive maintenance in machinery |
WO2003048701A2 (en) | 2001-12-04 | 2003-06-12 | Skf Condition Monitoring, Inc. | Cyclic time averaging for machine monitoring |
EP1474662B1 (en) | 2002-01-18 | 2011-07-27 | SPM Instrument AB | Analysis systems for analysing the condition of a machine |
US6618128B2 (en) | 2002-01-23 | 2003-09-09 | Csi Technology, Inc. | Optical speed sensing system |
US6668234B2 (en) | 2002-03-22 | 2003-12-23 | Abb Inc. | Method and apparatus for calculating the amplitude of a complex waveform associated with a rotating machine shaft after removing the running speed frequency |
JP3880455B2 (ja) | 2002-05-31 | 2007-02-14 | 中国電力株式会社 | 転がり軸受の余寿命診断方法及びこの余寿命診断装置 |
US7133801B2 (en) | 2002-06-07 | 2006-11-07 | Exxon Mobil Research And Engineering Company | System and methodology for vibration analysis and condition monitoring |
EP1513254A1 (en) | 2003-08-26 | 2005-03-09 | Mitsubishi Electric Information Technology Centre Europe B.V. | Filter enabling decimation of digital signals by a rational factor |
US7640139B2 (en) * | 2004-10-18 | 2009-12-29 | Nsk Ltd. | Abnormality diagnosing system for mechanical equipment |
NZ537244A (en) | 2004-12-16 | 2006-10-27 | Commtest Instr Ltd | Improvements in or relating to vibration analysis |
WO2007137132A2 (en) | 2006-05-17 | 2007-11-29 | Curtiss-Wright Flow Control Corporation | Probabilstic stress wave analysis system and method |
FR2913769B1 (fr) | 2007-03-12 | 2009-06-05 | Snecma Sa | Procede de detection d'un endommagement d'un roulement de palier d'un moteur |
GB0714379D0 (en) | 2007-07-21 | 2007-09-05 | Monition Ltd | Tamping bank monitoring apparatus and method |
US7761256B2 (en) | 2007-12-21 | 2010-07-20 | General Electric Company | Method and system for use in analyzing vibrations of a variable speed rotating body |
DK2085902T3 (en) | 2008-02-02 | 2018-03-12 | Siemens Ind Software Nv | Order tracking method and system |
TW200943037A (en) * | 2008-04-14 | 2009-10-16 | Asustek Comp Inc | Fan speed control method and computer system |
NO2370801T3 (zh) | 2008-12-22 | 2018-04-21 | ||
EP2373971B1 (en) | 2008-12-22 | 2022-01-26 | S.P.M. Instrument AB | An analysis system |
EA021908B1 (ru) | 2008-12-22 | 2015-09-30 | С.П.М. Инструмент Аб | Способ, устройство и машиночитаемый носитель для анализа состояния машины, имеющей вращающуюся часть |
CN102449445B (zh) * | 2009-05-05 | 2014-12-17 | S.P.M.仪器公司 | 用于分析具有旋转部件的机器振动的设备和方法 |
SE535559C2 (sv) | 2010-01-18 | 2012-09-25 | Spm Instr Ab | Förfarande och anordning för analys av tillståndet hos maskin med roterande del |
-
2010
- 2010-02-26 CN CN201080019737.8A patent/CN102449445B/zh active Active
- 2010-02-26 US US13/258,729 patent/US9964430B2/en active Active
- 2010-02-26 WO PCT/SE2010/050257 patent/WO2010128928A1/en active Application Filing
- 2010-02-26 AU AU2010245354A patent/AU2010245354A1/en not_active Abandoned
- 2010-02-26 EA EA201171353A patent/EA022630B1/ru not_active IP Right Cessation
- 2010-02-26 EP EP10772330.6A patent/EP2427745A4/en active Pending
-
2017
- 2017-12-13 US US15/840,720 patent/US10852179B2/en active Active
-
2020
- 2020-12-01 US US17/108,515 patent/US20210223094A1/en active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS60195426A (ja) * | 1984-03-16 | 1985-10-03 | Komatsu Ltd | 回転機構の故障診断方法 |
GB2190198A (en) * | 1986-04-28 | 1987-11-11 | Vipac Pty Ltd | Vibrational analysis system for a machine |
JPH01178814A (ja) * | 1988-01-08 | 1989-07-17 | Mitsubishi Heavy Ind Ltd | 信号の大局的構造認識装置 |
JPH04279826A (ja) * | 1991-03-08 | 1992-10-05 | Nkk Corp | 可変速回転系異常診断方法及びその装置 |
WO1998036251A1 (en) * | 1997-02-18 | 1998-08-20 | Dme Corporation | Method and apparatus for predictive diagnosis of moving machine parts |
JPH1178814A (ja) * | 1997-09-11 | 1999-03-23 | Honda Motor Co Ltd | 車両用多板式ブレーキ |
WO2002073150A2 (en) * | 2001-03-13 | 2002-09-19 | Ab Skf | System and method for analyzing vibration signals |
US20040101048A1 (en) * | 2002-11-14 | 2004-05-27 | Paris Alan T | Signal processing of multi-channel data |
JP2006113003A (ja) * | 2004-10-18 | 2006-04-27 | Nsk Ltd | 機械設備の異常診断システム |
CN1926413A (zh) * | 2004-10-18 | 2007-03-07 | 日本精工株式会社 | 机械设备的异常诊断系统 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108318122A (zh) * | 2017-01-16 | 2018-07-24 | 上银科技股份有限公司 | 自适应振动感测方法 |
CN108318122B (zh) * | 2017-01-16 | 2020-03-17 | 上银科技股份有限公司 | 自适应振动感测方法 |
CN107643182A (zh) * | 2017-09-08 | 2018-01-30 | 安徽大学 | 基于欠采样信号的轴承故障诊断方法及无线网络节点装置 |
CN110546477A (zh) * | 2018-03-29 | 2019-12-06 | 三菱重工业株式会社 | 接触振动检测装置及具备该接触振动检测装置的旋转机械、以及接触振动检测方法 |
CN110546477B (zh) * | 2018-03-29 | 2021-12-14 | 三菱重工船用机械株式会社 | 接触振动检测装置及具备该接触振动检测装置的旋转机械、以及接触振动检测方法 |
US11525753B2 (en) | 2018-03-29 | 2022-12-13 | Mitsubishi Heavy Industries Marine Machinery & Equipment Co., Ltd | Contact vibration detection device, rotary machine including the same, and contact vibration detection method |
CN109051618A (zh) * | 2018-09-06 | 2018-12-21 | 广东寰球智能科技有限公司 | 基于振动信号处理传送星轮故障的监测方法及监测系统 |
CN109051618B (zh) * | 2018-09-06 | 2020-03-31 | 广东寰球智能科技有限公司 | 基于振动信号处理传送星轮故障的监测方法及监测系统 |
CN113670428A (zh) * | 2021-07-19 | 2021-11-19 | 国网四川省电力公司电力科学研究院 | 一种变压器振动在线异常检测方法 |
CN113670428B (zh) * | 2021-07-19 | 2023-12-22 | 国网四川省电力公司电力科学研究院 | 一种变压器振动在线异常检测方法 |
Also Published As
Publication number | Publication date |
---|---|
US9964430B2 (en) | 2018-05-08 |
WO2010128928A1 (en) | 2010-11-11 |
AU2010245354A1 (en) | 2011-11-17 |
US20210223094A1 (en) | 2021-07-22 |
US20180164149A1 (en) | 2018-06-14 |
EP2427745A4 (en) | 2017-11-29 |
EP2427745A1 (en) | 2012-03-14 |
EA201171353A1 (ru) | 2012-10-30 |
US20120072136A1 (en) | 2012-03-22 |
EA022630B1 (ru) | 2016-02-29 |
US10852179B2 (en) | 2020-12-01 |
CN102449445B (zh) | 2014-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102449445A (zh) | 用于分析具有旋转部件的机器振动的设备和方法 | |
CN102308194B (zh) | 分析系统 | |
CN102308190A (zh) | 分析具有旋转部件的机器的状态的方法和系统 | |
CN104081173B (zh) | 用于分析具有旋转部件的机器的状态的装置和方法 | |
EP2370801B1 (en) | An analysis system | |
EP2373971B1 (en) | An analysis system | |
AU2015203361B2 (en) | An apparatus and a method for analysing the vibration of a machine having a rotating part | |
AU2016277551B2 (en) | Method and Apparatus for Analysing the Condition of a Machine Having a Rotating Part |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |