发明名称 一种基于概率神经网络的轴承健康状态识别方法
摘要 一种基于概率神经网络的轴承健康状态识别方法,步骤如下:1)采集原始信号;2)选择和评估参数:运用统计学原理计算震动信号产生的特征集,选取的统计特征,评估选取参数对已知数据分配的能力。3)构造概率神经网络模型;4)将输入参数输入概率神经网络进行状态识别。本发明将轴承的运行状态分为正常,亚健康和故障三种状态,采用概率神经网络(PNN)对轴承的健康状态进行估计,使用样本熵等作为概率神经网络(PNN)的输入特征参数,评估数据的分布能力,通过与传统轴承状态的正常-故障模型进行实验对比,并给出了轴承健康状态的识别结果。
申请公布号 CN103512751A 申请公布日期 2014.01.15
申请号 CN201310278309.2 申请日期 2013.07.03
申请人 辽宁大学 发明人 张利;田立;刘萌萌;陈朋杰;赵中洲
分类号 G01M13/04(2006.01)I;G06N3/02(2006.01)I 主分类号 G01M13/04(2006.01)I
代理机构 沈阳杰克知识产权代理有限公司 21207 代理人 罗莹
主权项 1.一种基于概率神经网络的轴承健康状态识别方法,其特征在于,步骤如下:1)采集原始信号:通过数据采集系统采集正常轴承和单点驱动端缺陷的数据;2)选择和评估参数:运用统计学原理计算震动信号产生的特征集,选取的统计特征如下:幅度:Range=max(x<sub>i</sub>)-min(x<sub>i</sub>)                    (8)有效值:<maths num="0001"><![CDATA[<math><mrow><msub><mi>&mu;</mi><mi>x</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>绝对平均值:<maths num="0002"><![CDATA[<math><mrow><msub><mi>&mu;</mi><mrow><mo>|</mo><mi>x</mi><mo>|</mo></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>|</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>均方值:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>&psi;</mi><mi>x</mi><mn>2</mn></msubsup><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>x</mi><mi>i</mi><mn>2</mn></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>均方根(RMS):<maths num="0004"><![CDATA[<math><mrow><msub><mi>&psi;</mi><mi>x</mi></msub><mo>=</mo><msqrt><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>x</mi><mi>i</mi><mn>2</mn></msubsup></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths>方差:<maths num="0005"><![CDATA[<math><mrow><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mi>x</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow></math>]]></maths>标准差:<maths num="0006"><![CDATA[<math><mrow><mi>SD</mi><mo>=</mo><msqrt><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mi>&mu;</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow></math>]]></maths>偏度:<maths num="0007"><![CDATA[<math><mrow><mi>S</mi><mo>=</mo><mfrac><msub><mi>M</mi><mn>3</mn></msub><msup><mi>&delta;</mi><mn>3</mn></msup></mfrac><mo>=</mo><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>[</mo><msup><mrow><mo>(</mo><mi>X</mi><mo>-</mo><mi>&mu;</mi><mo>)</mo></mrow><mn>3</mn></msup><mo>]</mo></mrow><msup><mi>&delta;</mi><mn>3</mn></msup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow></math>]]></maths>峰态:<maths num="0008"><![CDATA[<math><mrow><mi>K</mi><mo>=</mo><mfrac><msub><mi>M</mi><mn>4</mn></msub><msup><mi>&delta;</mi><mn>4</mn></msup></mfrac><mo>=</mo><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>[</mo><msup><mrow><mo>(</mo><mi>X</mi><mo>-</mo><mi>&mu;</mi><mo>)</mo></mrow><mn>4</mn></msup><mo>]</mo></mrow><msup><mi>&delta;</mi><mn>4</mn></msup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow></math>]]></maths>峰值:peak=max(X<sub>i</sub>)                   (17)波形指标:<maths num="0009"><![CDATA[<math><mrow><mi>WI</mi><mo>=</mo><mfrac><mi>peak</mi><msub><mi>&mu;</mi><mi>x</mi></msub></mfrac><mo>=</mo><mfrac><mrow><mi>max</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mi>N</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>x</mi><mi>i</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>18</mn><mo>)</mo></mrow></mrow></math>]]></maths>脉冲指标:<maths num="0010"><![CDATA[<math><mrow><mi>II</mi><mo>=</mo><mfrac><mi>peak</mi><msub><mi>&mu;</mi><mrow><mo>|</mo><mi>x</mi><mo>|</mo></mrow></msub></mfrac><mo>=</mo><mfrac><mrow><mi>max</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><mrow><mn>1</mn><mo>/</mo><mi>N</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover></mrow><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>|</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>19</mn><mo>)</mo></mrow></mrow></math>]]></maths>样本熵:其计算过程如下步骤一:N点序列,首先计算N*N的距离矩阵D,第i行第j列的元素就是d,并定义成点i,j之间的距离。<maths num="0011"><![CDATA[<math><mrow><msub><mi>d</mi><mi>ij</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>1</mn><mo>,</mo><mo>|</mo><mi>u</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mi>u</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo><mo>&lt;</mo><mi>r</mi></mtd></mtr><mtr><mtd><mn>0</mn><mo>,</mo><mi>u</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><mi>u</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>|</mo><mo>&GreaterEqual;</mo><mi>r</mi></mtd></mtr></mtable></mfenced><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>;</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>,</mo><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>20</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤二:用矩阵D中的元素很容易得到<img file="FDA00003457188700032.GIF" wi="58" he="82" />和<img file="FDA00003457188700033.GIF" wi="89" he="92" />假设窗长度是2;<maths num="0012"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>i</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>2</mn></mrow></munderover><msub><mi>d</mi><mi>ij</mi></msub><mo>&cap;</mo><msub><mi>d</mi><mrow><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>21</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0013"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>i</mi><mn>3</mn></msubsup><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>2</mn></mrow></munderover><msub><mi>d</mi><mi>ij</mi></msub><mo>&cap;</mo><msub><mi>d</mi><mrow><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>&cap;</mo><msub><mi>d</mi><mrow><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mrow><mo>(</mo><mi>j</mi><mo>+</mo><mn>2</mn><mo>)</mo></mrow></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>22</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤三:通过<img file="FDA00003457188700036.GIF" wi="66" he="85" />和<img file="FDA00003457188700037.GIF" wi="56" he="82" />计算得到<img file="FDA00003457188700038.GIF" wi="66" he="101" />和<img file="FDA00003457188700039.GIF" wi="78" he="102" /><maths num="0014"><![CDATA[<math><mrow><msup><mi>&phi;</mi><mi>m</mi></msup><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mo>-</mo><mi>m</mi><mo>+</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mi>m</mi><mo>+</mo><mn>1</mn></mrow></munderover><mi>ln</mi><msubsup><mi>C</mi><mi>i</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>23</mn><mo>)</mo></mrow></mrow></math>]]></maths>步骤四:计算样本熵<maths num="0015"><![CDATA[<math><mrow><mi>SampEn</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>r</mi><mo>,</mo><mi>N</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mi>ln</mi><mfrac><mrow><msup><mi>&phi;</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow></mrow><mrow><msup><mi>&phi;</mi><mi>m</mi></msup><mrow><mo>(</mo><mi>r</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>24</mn><mo>)</mo></mrow></mrow></math>]]></maths>评估上述13个参数对已知数据分配的能力,并另外选取性能比较好的绝对值平均值,峰态,标准差,均方根和样本熵作为输入特征。3)构造概率神经网络模型:输入层有5个节点表示输入特征,求和层和输出层都是3个节点,表示健康的三个状态,模式层有150个节点,表示选取的150个样本;4)将输入参数输入概率神经网络进行状态识别,根据健康度的定义HD=P(healthy)-P(failure),计算健康度的值,根据轴承健康度的阈值来判定轴承所处的运行状态为健康、亚健康或故障。当0.4&lt;HD&lt;1时,轴承为健康状态,当-0.4&lt;HD&lt;0.4时,轴承为亚健康状态,当-1&lt;HD&lt;-0.4时,轴承为故障状态。
地址 110136 辽宁省沈阳市沈北新区道义南大街58号