振动信号的采集与预处理.docx

振动信号的采集与预处理几乎所有的物理现象都可看作是信号，但这里我们特指动态振动信号。振动信号采集与一般性模拟信号采集虽有共同之处，但存在的差异更多，因此，在采集振动信号时应注意以下几点：1. 振动信号采集模式取决于机组当时的工作状态，如稳态、瞬态等；2. 变转速运行设备的振动信号采集在有条件时应采取同步整周期采集；3. 所有工作状态下振动信号采集均应符合采样定理。对信号预处理具有特定要求是振动信号本身的特性所致。信号预处理的功能在一定程度上说是影响后续信号分析的重要因素。预处理方法的选择也要注意以下条件：1. 在涉及相位计算或显示时尽量不采用抗混滤波；2. 在计算频谱时采用低通抗混滤波；3. 在处理瞬态过程中 1X 矢量、2X 矢量的快速处理时采用矢量滤波。上述第 3 条是保障瞬态过程符合采样定理的基本条件。在瞬态振动信号采集时，机组转速变化率较高，若依靠采集动态信号（一般需要若干周期）通过后处理获得1X 和 2X 矢量数据，除了效率低下以外，计算机（服务器）资源利用率也不高，且无法做到高分辨分析数据。机组瞬态特征（以波德图、极坐标图和三维频谱图等型式表示）是固有的，当组成这些图谱的数据间隔过大（分辨率过低）时，除许多微小的变化无法表达出来，也会得出误差很大的分析结论，影响故障诊断的准确度。一般来说，三维频谱图要求数据的组数（rpm 分辨率）较少，太多了反而影响对图形的正确识别；但对前面两种分析图谱，则要求较高的分辨率。目前公认的方式是每采集 10 组静态数据采集 1 组动态数据，可很好地解决不同图谱对数据分辨率的要求差异。影响振动信号采集精度的因素包括采集方式、采样频率、量化精度三个因素，采样方式不同，采集信号的精度不同，其中以同步整周期采集为最佳方式；采样频率受制于信号最高频率；量化精度取决于 A/D 转换的位数，一般采用 12 位，部分系统采用 16 位甚至 24 位。振动信号的采样过程，严格来说应包含几个方面：1. 信号适调由于目前采用的数据采集系统是一种数字化系统，所采用的 A/D 芯片对信号输入量程有严格限制，为了保证信号转换具有较高的信噪比，信号进入 A/D 以前，均需进行信号适调。适调包括大信号的衰减处理和弱信号的放大处理，或者对一些直流信号进行偏置处理，使其满足 A/D 输入量程要求。2. A/D 转换 A/D 转换包括采样、量化和编码三个组成部分。采样（抽样），是利用采样脉冲序列p(t)从模拟信号 x(t)中抽取一系列离散样值，使之成为采样信号 x(nt)（n=0,1,2,）的过程。t 称为采样间隔，其倒数称1/tf 之为采样频s率。采样频率的选择必须符合采样定理要求。由于计算机对数据位数进行了规定，采样信号 x(nt)经舍入的方法变为只有有限个有效数字的数，这个过程称为量化。由于抽样间隔长度是固定的（对当前数据来说），当采样信号落入某一小间隔内，经舍入方法而变为有限值时，则产生量化误差。如 8 位二进制为2 256，即量化增量为所测信号最大电压幅值的 1/256。8编码是把采样数据转变为计算机能识别的数字格式。一、采样定理1. 采样定理采样定理解决的问题是确定合理的采样间隔 t 以及合理的采样长度T，保障采样所得的数字信号能真实地代表原来的连续信号x(t)。衡量采样速度高低的指标称为采样频率 f 。一般来说，采样频率 f 越高，采样点越密，ss所获得的数字信号越逼近原信号。为了兼顾计算机存储量和计算工作量，一般保证信号不丢失或歪曲原信号信息就可以满足实际需要了。这个基本要求就是所谓的采样定理，是由Shannon 提出的，也称为 Shannon 采样定理。Shannon 采样定理规定了带限信号不丢失信息的最低采样频率为f ³ 2 f 或w ³ 2wsmsm式中 f 为原信号中最高频率成分的频率。m采集的数据量大小 N 为TN =Dt因此，当采样长度一定时，采样频率越高，采集的数据量就越大。使用采样频率时有几个问题需要注意。一，正确估计原信号中最高频率成分的频率，对于采用电涡流传感器测振的系统来说，一般确定为最高分析频率为12.5X，采样模式为同步整周期采集，若选择频谱分辨率为400 线，需采集1024 点数据，若每周期采集32 点，采样长度为 32 周期。二，同样的数据量可以通过改变每周期采样点数提高基频分辨率，这对于识别次同步振动信号是必要的，但降低了最高分析频率，如何确定视具体情况而定。2. 采样定理解析采样定理实际上涉及了 3 个主要条件，当确定其中 2 个条件后，第 3 个条件自动形成。这 3 个条件是进行正确数据采集的基础，必须理解深刻。条件 1 采样频率控制最高分析频率采样频率（采样速率）越高，获得的信号频率响应越高，换言之，当需要高频信号时，就需要提高采样频率，采样频率应符合采样定理基本要求。这个条件看起来似乎很简单，但对于一个未知信号，其中所含最高频率信号的频率究竟 有多高，实际上我们是无法知道的。解决这个问题需要2 个步骤，一是指定最高测量频率，二是采用低通滤波器把高于设定最高测量频率的成分全部去掉（这个低通滤波器就是抗混滤波器）。现实的抗混滤波器与理论上的滤波器存在差异，因此信号中仍会存在一定混叠成分，一般在计算频谱后将高频成分去掉，一般频谱线数取时域数据点的 1/2.56，或取频域幅值数据点的 1/1.28，即 128 线频谱取 100 线，256 线频谱取 200 线，512 线频谱取 400 线等等。图 、采样过程示意图抗混滤波器的使用主要是针对频谱分析的，对于涉及相位计算的用途反而会引入相位误 差。几乎所有的滤波器的相位特性远比幅值特性差。为说明该条件，我们举例进行说明。 要想在频谱中看到 500Hz 的成分，其采样频率最少为 1000Hz。 若采样频率为 32 点转，频谱中最高线理论上可达到 16X。条件 2 总采样时间控制分辨率频谱的分辨率（谱线间隔）受控于总采样时间，即1Df =T其中f 为频谱分辨率，T 为总采样时间。 如果采样总时间为 0.5 秒，则频谱分辨率为 2Hz； 若区分 6cpm（0.1Hz）的频谱成分，则总采样时间至少为 10 秒； 对于总采样时间为 8 转的时间信号，频谱分辨率为 1/8X。条件 3 采样点数控制频谱线数解释这个条件，需要对FFT 计算频谱的过程有一个了解。如果对于一个2048 点的时间波形数据，我们可以获得2048 点频域数据1024 线频谱（每条谱线有两个值，直接值和正交值，或者说幅值和相位两个值）。对旋转机械来说，频谱仅仅画出了FFT 复数输出的幅值部分，对于相位部分一般不画，因此频谱中的线数最多为时域点数的一半，考虑到混叠的影响，频谱线数一般会低于时域数据点数。QuadraturePhaseDirect小结采样定理是实现正确采样的基准，上述 3 个条件中，可以根据需要设置其中 2 个条件，第 3 个条件就会自动固定。 如果采样总时间为 0.5 秒，想获得 3200 线频谱，则有11Df = 2Hz条件 2T0.5sec条件 3 3200 线频谱实际需要 4096 点频谱数据（考虑到混叠问题），8192 点时域数据= 8192 / 0.5 =16384 f（Hz）s f =16384 / 2 = 8192HzMaxf= 3200lines ´ 2Hz / line = 6400Hz3200 若在频谱上能区分 0.2Hz 间隔的频率成分，频谱确定为 800 线，则有11= 5（秒）条件 2 TDf0.2条件 3 800 线频谱实际需要 1024 点频域数据，2048 点时域数据= 2048 ¸ 5 = 409.6 f（Hz）sf = 409.6 / 2 = 204.8 （Hz）Maxf = 800lines ´0.2Hz / line =160Hz800 若在频谱上能区分 0.1Hz 间隔的频率成分，且能在频谱上最大看到180Hz，则有条件 1 f条件 2 T³ 2 f ³ 2´180 ³ 360（Hz）Maxs11=0.1=10（秒）Df因此，按不低于 360 点秒的采样速率采集 10 秒钟，可采集时域数据最少 3600 点。为方便 FFT 计算，数据点数应为 2 的整数次幂，与 3600 最接近的数值是 4096，由此可获得 2048 点频域数据，即可获得 1600 线频谱。1600 线、频率间隔为 0.1Hz 的频谱最高分析频率为 160Hz，显然不能满足需要。4096 下一个 2 的整数次幂的数值是 8192，由此可获得 3200 线的频谱，其最高分析频率达到了 320Hz，可以满足要求，可以通过提高采样速率来实现这一要求。采样频率（Hz） ： 分析频率（上限频率）（Hz） * 2.5651210242048HzHz1020405012848161002562482005121245001280256051201280025600512001280000.40.20.10.040.020.010.0040.80.40.20.080.040.020.0081.60.80.40.160.080.040.0161000200050001000020000500001052.510025012562.5 1000002560000.0025000.0042500.008125 在同步整周期采样时，若采集 32 点转，共采集 8 转，则可获得 256 点时域数据和100 线频谱，有T=8 转f1/T1/8 转1/8Xf 32 点转÷216XMaxf 100 线×0.125 转线12.5X100用通用的方式表达为设x （n0，1,2,N-1）为一采样序列若每周期等角度采集m 点，共采集 L 周，则n有mLN设该旋转机械的转动频率为 f，则采样间隔为1Dt =fm变换后的频率分辨率1fm fDf =NDt N L或f = DfL显然，工频分量正好处于第 L 条线上。相应地，kfkfL，即第 k 阶分量也处于整数f 上，这样就保证了特征频率成分在频谱上的准确定位。采用同步整周期采样可获得的最高分析频率为NNfmff =mDf =2.562.56L 2.56【问题 01】在不重采样的前提下，能否提高频谱的分辨率？【问题 02】对于旨在分析齿轮故障和滚动轴承故障的振动数据采集应注意什么事项？理解了采样定理的实质，我们就会对某些仪器系统中列出的技术指标有了正确的认识，频谱分辨率并不是衡量采样质量的唯一指标，即 400 线频谱与 400 线的频谱之间有可能存在差异；在分析齿轮故障时就不会出现没有啮合频率成分这样的尴尬；在分析开停车过程时出现分辨率过低问题频谱的功率泄漏问题数字信号分析需要选择取合理的采样长度，虽然在采样过程中充分考虑了采样定理和分析要求，但毕竟是一个用区间为（TT）的有限长度信号来近似 t信号的过程， |t|T 的 x(t)值为零，因此所得到的频谱和实际频谱存在一定差异，这种现象称之为泄漏现象。 l 影响数据采集过程的几个关键环节l A/D 转换位数（转换精度）l 采样方式（自由采集与整周期采集、同步自由采集与同步整周期采集）l 数据采集的效率l 数据采集中相位信号的作用振动信号的预处理1. 低通抗混滤波抗混滤波器是一种低通滤波器，如广泛采用的 8 阶椭圆滤波器。在线系统采样单元中采用的抗混滤波器，应具备截至频率可跟踪性，即随着机器转速的变化，低通滤波器的截至频率也随之变化。抗混滤波器的使用目的是避免频谱分析时高频分量折叠到低频段，但每一种低通滤波器的相频响应曲线并不像幅频曲线那样平坦。如某型号 4 阶低通滤波器，在其截至频率处相位偏移达180°，在 3kHz 处偏移也有45°，这个相位偏移足以影响幅值的正确计量。图 、某低通滤波器的响应特性曲线一般来说，用于幅值计量、相位计算、轨迹显示等用途时，最好不使用抗混滤波器。2. 数字矢量滤波数字矢量滤波器是一种用特殊数字技术实现的跟踪带通滤波器，可实现 1X、2X 或设定分频矢量值的功能。数字矢量滤波器在低速和高速下的带宽设置不同，低速时带宽设置窄一些，在高速时带宽适当增大。通带越窄，需要的响应时间越长，因此债的通带只适用于低速段。如本特利公司的 DVF3 型数字矢量滤波器设置的带宽及响应时间为：表 、DVF3 滤波参数设置表转速范围1000rpm1000rpm带宽范围120±10cpm12±1cpm7.30 秒 fc图、带通滤波器3. 积分电路当采用壳体振动测量传感器时，用户可能需要使用某种特定的振动单位作为监测参量，以方便与相关标准对应。当采用加速度传感器测量时，积分电路就是必要的预处理手段。积分器位移速度加速度图 、振动信号积分流程图需要说明的是，对于壳体振动信号，我们可以通过积分将加速度信号转换为速度信号，或进一步积分转换为位移，但仅对壳体振动信号有效，不能把壳体振动位移与转子相对振动位移相比较，因为它们之间不存在规律性的关系。另外需要注意的是，加速度积分成的速度与速度传感器直接测量的速度是存在区别的。 图 、绝对位移相对位移4. 泄漏数字信号分析需要选取合理的采样长度，这个长度就是数据采样对原始信号的截断，若截断长度为（TT），则对于|t|T 的 x(t)值均为零，因此得到的频谱与实际频谱存在一定差异，表现为频谱上出现旁瓣，主瓣的幅值与输入的时间信号幅值产生了差异（降低）。为了降低时间波形因截断引起的频域振荡现象，在信号分析时多采用加窗的方式进行修正，窗函数就是截断函数，不同的窗函数对旁瓣的抑制能力也有区别。一般来说，窗函数必须具备以下基本要求：窗谱的主瓣要窄且高，旁瓣要小，正负交替接近相等，以减小泄漏或负谱现象。常用的窗函数及其性能指标如表 所示。表 、典型窗函数的性能特点旁瓣衰减dB/10oct20606020201.28B1.33B1.20B1.23B1.30B1.36B高斯1.55B1.64B55窗函数是数据处理时要考虑的因素，但因与数据采集有关，在此先作提示。5. 其它问题与信号采集有关的问题还涉及到系统的输入阻抗与输入阻抗的问题。对于框架式仪表提供的缓冲输出信号，其输出阻抗多在几百欧姆左右，如本特利公司的7200 系列仪表和 3300 系列仪表的缓冲输出阻抗为 100 欧姆，3500 系列仪表的缓冲输出阻抗为 550 欧姆。数据采集系统的输入阻抗不易过低，理论上讲输入阻抗值越高越好。要保证信号损失率 在 1%以下，输入输出阻抗比应控制在 100:1 以上。信号损失率的定义为：éùZSL(%) = 1-´100%êúinZ + Zoutëûin如某系统信号输出阻抗为 9090 欧姆，接入一个输入阻抗为 100k 欧姆的系统，其信号损失率为：éù100000Zé´100% = 1-êùúSL(%) = 1-´100% = 8.33%êúinZ + Zoutë 9090 +100000ûëûin这样的信号损失率显然是不符合工程要求的。