主权项 |
一种马尔可夫故障趋势预测方法,该方法基于1.5维谱频带能量均值实现,其包括以下步骤:(1)利用转子实验台模拟旋转机械设备正常运行状态,利用现有数据采集设备采集转子实验台在正常运行状态下的振动信号x<sub>w</sub>(n)={x<sub>1</sub>,...,x<sub>N</sub>},其中,N代表每组数据个数,w代表数据组别,w=1,表示正常运行状态;(2)利用转子实验台模拟旋转机械设备故障的轻度故障程度、中度故障程度和重度故障程度三种故障程度,并利用现有数据采集设备采集转子实验台在三种故障下的振动信号x<sub>w</sub>(n)={x<sub>1</sub>,...,x<sub>N</sub>},其中,N代表每组数据个数;w=2、3、4代表数据组别,w=2代表轻度故障程度状态,w=3代表中度故障程度状态,w=4代表重度故障程度状态;(3)计算所有振动信号x<sub>w</sub>(n)中每组振动信号的1.5维谱;(4)计算各组振动信号的1.5维谱频带能量均值<img file="FDA0001110907670000011.GIF" wi="115" he="68" /><maths num="0001"><math><![CDATA[<mrow><msub><mover><mi>S</mi><mo>‾</mo></mover><mrow><mi>w</mi><mo>,</mo><mn>3</mn><mi>x</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>S</mi><mrow><mi>w</mi><mo>,</mo><mn>3</mn><mi>x</mi></mrow></msub><mrow><mo>(</mo><msub><mi>ω</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001110907670000012.GIF" wi="469" he="134" /></maths>式中,S<sub>w,3x</sub>(ω<sub>r</sub>)为振动信号的1.5维谱;ω<sub>r</sub>代表频率,r=1,2,…,N,N为正整数;(5)获得频带能量区间:根据1.5维谱频带能量均值<img file="FDA0001110907670000013.GIF" wi="91" he="74" />对旋转机械设备故障的三种故障程度进行量化,获得划分故障状态的频带能量区间:频带能量均值数值在区间<img file="FDA0001110907670000014.GIF" wi="172" he="79" />的状态分类为正常运行状态,记为状态1;频带能量均值数值在区间<img file="FDA0001110907670000015.GIF" wi="238" he="79" />的状态分类为轻度故障状态,记为故障状态2;频带能量均值数值在区间<img file="FDA0001110907670000016.GIF" wi="243" he="79" />的状态分类为中度故障状态,记为故障状态3;频带能量均值数值在区间<img file="FDA0001110907670000017.GIF" wi="258" he="83" />的状态分类为重度故障状态,记为故障状态4;频带能量均值数值在区间<img file="FDA0001110907670000018.GIF" wi="234" he="82" />的状态分类为破坏状态,记为故障状态5;记状态序列为:{S<sub>1</sub>,S<sub>2</sub>,...,S<sub>n</sub>},n为状态数目,为正整数;状态空间E={1,2,...,5};(6)采集实际旋转机械设备的振动信号,对该振动信号进行步骤(3)至步骤(4)的操作,得到各组振动信号的1.5维谱频带能量均值,结合步骤(5)给出的频带能量区间,获得实际旋转机械设备的状态序列:{S<sub>1</sub>,S<sub>2</sub>,...,S<sub>n</sub>};(7)利用马尔可夫链对实际旋转机械设备的状态进行趋势预测。 |