发明名称 一种机械故障程度的评价方法
摘要 本发明涉及一种机械故障程度的评价方法,其步骤:(1)利用转子实验台模拟机械设备正常运行状态,采用振动信号;(2)利用转子实验台模拟机械设备故障的轻度故障程度、中度故障程度和重度故障程度,采集振动信号;(3)对所有振动信号进行双谱操作获得双谱值矩阵;(4)将双谱值矩阵定义为劣化演化矩阵;计算正常运行状态下的劣化演化矩阵以及故障下的劣化演化矩阵;(5)计算劣化演化矩阵的特征值、均值、方差;(6)计算机械设备正常运行状态下的均值最大值、方差最大值;(7)根据以上获得劣化区间划分图;(8)采集实际设备的振动信号得到劣化均值、劣化方差,将其绘制在劣化区间划分图中,根据坐标点所处的区域判断设备运行状态。
申请公布号 CN104155133A 申请公布日期 2014.11.19
申请号 CN201410384028.X 申请日期 2014.08.06
申请人 北京信息科技大学 发明人 蒋章雷;左云波;徐小力;吴国新;谷玉海
分类号 G01M99/00(2011.01)I 主分类号 G01M99/00(2011.01)I
代理机构 北京远大卓悦知识产权代理事务所(普通合伙) 11369 代理人 贺持缓
主权项 一种机械故障程度的评价方法,其包括以下步骤:(1)利用转子实验台模拟机械设备正常运行状态,并采用现有数据采集设备采集转子实验台在正常运行状态下的振动信号x<sub>w</sub>(n)={x<sub>1</sub>,…x<sub>N</sub>},其中,N为每组数据个数;w为数据组别,w=1,…,9,即1至9组为正常运行状态;(2)利用转子实验台模拟机械设备故障的轻度故障程度、中度故障程度和重度故障程度三种故障程度,并采用现有数据采集设备采集转子实验台在三种故障下的振动信号x<sub>w</sub>(n)={x<sub>1</sub>,…x<sub>N</sub>},其中,N代表每组数据个数;w=10,11,12代表数据组别,即第10组数据代表轻度故障程度状态、第11组数据代表中度故障程度状态、第12组数据代表重度故障程度状态;(3)对所有振动信号x<sub>w</sub>(n)进行双谱操作,获得不同运行状态下的双谱值矩阵F<sub>w</sub>(ω<sub>1</sub>,ω<sub>2</sub>);(4)将描述各个频率点双谱值增减幅度的双谱值矩阵定义为劣化演化矩阵F<sub>pq</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>F</mi><mi>pq</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>F</mi><mi>q</mi></msub><mo>-</mo><msub><mi>F</mi><mi>p</mi></msub></mrow><msub><mi>F</mi><mi>p</mi></msub></mfrac><mo>)</mo></mrow><mo>*</mo></msup><mo>,</mo></mrow>]]></math><img file="FDA0000549977720000011.GIF" wi="362" he="152" /></maths>式中,F<sub>q</sub>、F<sub>p</sub>分别表示状态q、状态p对应的双谱值矩阵,p∈w,q∈w;“*”表示分子位置处矩阵与分母位置处矩阵的相除操作为两个矩阵各个对应元素之间的相除;计算机械设备正常运行状态下的劣化演化矩阵:F<sub>1‑2</sub>、F<sub>1‑3</sub>、F<sub>1‑4</sub>、F<sub>1‑5</sub>、F<sub>1‑6</sub>、F<sub>1‑7</sub>、F<sub>1‑8</sub>、F<sub>1‑9</sub>,以及故障下的劣化演化矩阵:F<sub>1‑10</sub>、F<sub>1‑11</sub>、F<sub>1‑12</sub>;(5)计算劣化演化矩阵F<sub>pq</sub>的特征值,得到12组特征值,再计算每组特征值的均值、方差:<img file="FDA0000549977720000012.GIF" wi="1730" he="84" /><img file="FDA0000549977720000013.GIF" wi="1308" he="84" />其中a表示劣化均值,b表示劣化方差;(6)计算机械设备正常运行状态下对应的均值最大值a<sub>max</sub>为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>a</mi><mi>max</mi></msub><mo>=</mo><mo>{</mo><msub><mi>a</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>2</mn></mrow></msub></msub><mo>,</mo><msub><mi>a</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>3</mn></mrow></msub></msub><mo>,</mo><msub><mi>a</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>4</mn><mo>,</mo></mrow></msub></msub><msub><mi>a</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>5</mn></mrow></msub></msub><mo>,</mo><msub><mi>a</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>6</mn></mrow></msub></msub><msub><mrow><mo>,</mo><mi>a</mi></mrow><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>7</mn></mrow></msub><mtext></mtext></msub><mo>,</mo><msub><mi>a</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>8</mn></mrow></msub></msub><mo>,</mo><msub><mi>a</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>9</mn></mrow></msub></msub><mo>}</mo><mo>;</mo></mrow>]]></math><img file="FDA0000549977720000014.GIF" wi="1096" he="78" /></maths>计算机械设备正常运行状态下对应的方差最大值b<sub>max</sub>为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>b</mi><mi>max</mi></msub><mo>=</mo><mo>{</mo><msub><mi>b</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>2</mn></mrow></msub></msub><mo>,</mo><msub><mi>b</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>3</mn></mrow></msub></msub><mo>,</mo><msub><mi>b</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>4</mn></mrow></msub></msub><mo>,</mo><msub><mi>b</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>5</mn></mrow></msub></msub><mo>,</mo><msub><mi>b</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>6</mn></mrow></msub></msub><msub><mrow><mo>,</mo><mi>b</mi></mrow><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>7</mn></mrow></msub><mtext></mtext></msub><mo>,</mo><msub><mi>b</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>8</mn></mrow></msub></msub><mo>,</mo><msub><mi>b</mi><msub><mi>F</mi><mrow><mn>1</mn><mo>-</mo><mn>9</mn></mrow></msub></msub><mo>}</mo><mo>;</mo></mrow>]]></math><img file="FDA0000549977720000015.GIF" wi="1109" he="73" /></maths>(7)以劣化均值为横坐标,以劣化方差为纵坐标,将均值最大值a<sub>max</sub>、方差最大值b<sub>max</sub>;故障状态下的劣化均值、劣化方差<img file="FDA0000549977720000016.GIF" wi="834" he="79" />的数值绘制在同一个图形中,得到劣化区间划分图;(8)采集实际设备的振动信号,对该振动信号进行步骤(3)~(6)的操作,得到对应的劣化均值、劣化方差,将其绘制在步骤(7)中的劣化区间划分图中,根据坐标点所处的区域判断设备运行状态。
地址 100192 北京市海淀区清河小营东路12号166信箱
您可能感兴趣的专利