一种基于马尔科夫模型电机系统可靠性定量评估方法与流程

本发明涉及一种可靠性定量评估方法，尤其是适用于各种类型、各种拓扑结构的电机系统可靠性马尔科夫模型的改进方法。
背景技术：
：马尔科夫模型(markovmodel)是对随机过程可靠性评估的一个有效的方法，可以清楚的描述出系统从健康状态到失效状态的转变过程。电机系统是由很多个零件组成，每一个零件可能发生故障由多种。当任意零件发生故障后，电机系统的性能都会发生变化，电机系统的状态就随着发生变化。通过判断电机在发生故障后的性能是否满足工作的需求来判断电机系统在发生故障后是否可靠，尤其是对于高容错性能的电机系统而言，具有较强的带故障运行能力。电机系统的某个零件发生一个故障后，系统仍然能够可靠运行，在此基础上系统发生第二个故障。如果还能够继续可靠运行，系统再发生第三个故障。系统分析能够带故障数量的情况越详细可靠性定量评估的结果越准确。传统的可靠性框图方法在考虑了健康状态和失效状态没有考虑故障工作状态。一般的马尔科夫模型可靠性方法考虑最多的状态只有两到三个故障情况，并且其计算过程需要计算常微分方程，计算过程复杂，不能快速直观的得到每个状态和系统的可靠度函数。技术实现要素：本发明的目的是为了更加准确对系统的可靠性进行定量评估，减小计算量，快速计算每个状态发生的概率，提出了一种改进的马尔科夫模型计算方法。为达到上述技术目的，本发明是通过一下技术方案实现的：一种基于马尔科夫模型电机系统可靠性定量评估方法，具体实现步骤如下：s1、确定电机系统在发生一个故障、两个故障、i(i＞2)个故障下的系统所有可能运行情况，通过判断其是否满足可靠性标准来确定系统在每个故障情况下的可靠性状态；s2、假设根据步骤1得到系统一共有n个状态，其中有一个健康状态、一个失效状态和n-2个故障工作状态。x个一个故障工作状态，y个两个故障工作状态，z个三个故障工作状态等等各种故障工作状态。建立电机系统的马尔科夫状态转移图以及状态转移矩阵a；s3、依据状态转移矩阵a中的元素确定每个可靠工作状态之间的影响系数矩阵r，根据公式r(t)＝r·e快速计算得到系统的可靠度函数r(t)，根据平均故障前时间公式计算系统的平均故障前时间数值，完成系统可靠性定量评估。作为优选，所述步骤s2中，状态转移矩阵a的转置表示为：状态转移矩阵a的转置是n行n列的矩阵，矩阵中的每一个元素的含义是两个状态之间的转移概率，状态本身的转移概率等于本身转移出去的所有转移概率之和的相反数。作为优选，所述步骤s3中，影响系数矩阵r为：影响系数矩阵r是n-1行n-1列的下三角矩阵，矩阵的每一行矩阵表示每一个可靠工作状态，每一个元素的含义是先发生的状态对后发生状态的影响系数，状态j对状态i的影响系数ri，j表示为：初始状态条件是状态1发生的概率等于1，其他状态发生的概率等于0，即p(0)＝[10…0]。影响系数矩阵r的每一行依次对应一个可靠工作状态，每一行元素之和等于系统初始时每个状态发生的概率。初始情况下状态1即健康状态发生的概率等于1，所以r1，1＝1，其他状态在初始情况下发生的概率等于0，例如r2，1+r2，2＝0。可靠度函数r(t)计算公式为：式中t表示时间。快速计算得到系统的可靠度函数r(t)。由系统可靠度函数r(t)计算得到系统的平均故障前时间：从而计算得到系统的可靠性定量评估。本发明采用上述技术方案，具有以下有益效果：(1)改进的方法在可靠性定量评估中对电机系统所有可能发生的故障工作状态都进行了分析。(2)改进的方法计算量显著减小。(3)改进的方法能够快速的直观的得到每一个可靠工作状态发生的概率，和状态之间的影响系数。附图说明图1为本发明的n状态电机系统的马尔科夫状态转移图。图2为本发明实施某一台三相永磁电机系统的马尔科夫状态转移图。图3为本发明实施某一台三相永磁电机和三相逆变器组成的三相永磁电机系统示意图。图4为本发明实施的某一台三相永磁电机系统马尔科夫改进方法解得的可靠度随时间变化曲线图。具体实施方式下面结合实例和附图对本发明做进一步的详细说明。系统同时发生多个故障的情况下，系统状态数目为n时的马尔科夫状态转移图如图1所示。图1中状态1为健康状态，状态2到状态x为一个故障工作状态，状态x+1到状态y为两个故障工作工作状态，状态n为失效状态，其他的故障工作状态在图中省略。以一个三相永磁电机系统为例来讨论一种改进的马尔科夫模型计算方法的使用过程。s1、首先确定永磁电机系统可能发生的主要故障，如表1所示。其次确定电机系统发生一个故障时的可靠工作状态情况，如表2所示。在系统发生一个故障的可靠工作的基础上模拟发生第二个故障情况得到系统发生两个故障时的可靠工作状态情况，如表3所示。在系统发生两个故障的基础上，电机系统再发生一个故障，即发生三个故障时发现系统都处于失效状态。也就是说此永磁电机系统最多发生两个故障。表1永磁电机系统可能发生的主要故障表2永磁电机系统发生一个故障时的状态故障模式状态故障模式状态igbt开路失效电流传感器信号偏移可靠igbt短路失效电流传感器信号放大可靠电容开路失效电流传感器损坏失效电容短路失效位置传感器损坏失效相间短路失效位置传感器信号丢失失效匝间短路可靠永磁体性能下降可靠绕组开路失效永磁体损坏可靠表3永磁电机系统先后发生两个故障时的状态表3中不发生表示这种情况不发生。s2、确定此电机系统所有可能发生的可靠工作状态后得到系统的马尔科夫状态转移图，如图2所示。图2中状态1到状态22所对应的故障情况见表2和3。状态1表示永磁电机系统处于健康状态，状态2到状态6表示永磁电机系统处于一个故障工作状态，状态7到状态22表示永磁电机系统处于两个故障工作状态，状态23表示永磁电机系统处于失效状态。表4马尔科夫状态转移图中每个状态的含义图2中状态之间的转移概率用a表示，下标表示对应状态转移情况，例如a2，1表示从状态1转移到状态2的状态转移概率，a3，1表示从状态1转移到状态3的状态转移概率，其他的状态转移概率含义相同。同时状态本身的转移概率的大小等于本身转移出去的所有转移概率之和的相反数，例如，a1，1等于a2，1、a3，1、a4，1、a5，1、a6，1和a23，1的和，其他状态本身的转移概率的大小相同的方法确定。根据故障发生的概率得到图2对应的状态转移矩阵a，它的转置为：表5矩阵a转置中每个非零元素的数值a1，1-1.046×10-5a2，2-9.36×10-6a3，3-1.036×10-5a4，4-1.036×10-5a5，5-9.1×10-6a6，6-8.46×10-6a7，7-9.26×10-6a8，8-8×10-6a9，9-7.36×10-6a10，10-9×10-6a11，11-8.36×10-6a12，12-9.26×10-6a13，13-9×10-6a14，14-8.36×10-6a15，15-8×10-6a16，16-9×10-6a1７，17-9×10-6a18，18-7.1×10-6a19，19-7.36×10-6a20，20-8.36×10-6a21，21-8.36×10-6a22，22-7.1×10-6a2，11.1×10-6a3，10.1×10-6a4，10.1×10-6a5，11.36×10-6a6，12×10-6a23，15.8×10-6a7，20.1×10-6a8，21.36×10-6a9，22×10-6a10，31.36×10-6a11，32×10-6a12，41.1×10-6a13，41.36×10-6a14，42×10-6a15，51.1×10-6a16，50.1×10-6a17，50.1×10-6a18，52×10-6a19，61.1×10-6a20，60.1×10-6a21，60.1×10-6a22，61.36×10-6a23，25.9×10-6a23，37×10-6a23，45.9×10-6a23，55.8×10-6a23，65.8×10-6a23，79.26×10-6a23，88×10-6a23，97.36×10-6a23，109×10-6a23，118.36×10-6a23，129.26×10-6a23，139×10-6a23，148.36×10-6a23，158×10-6a23，169×10-6a23，179×10-6a23，187.1×10-6a23，197.36×10-6a23，208.36×10-6a23，218.36×10-6a23，227.1×10-6s3、根据状态转移矩阵a，并利用公式：同时考虑系统初始情况下各个状态发生的概率，即p(0)＝[10…0]。计算得到每个可靠工作状态之间的影响系数矩阵r，影响系数矩阵r为：表6矩阵r中每个非零元素的数值r1，11r2，21r3，31r4，41r5，51r6，61r7，70.917r8，80.447r9，90.355r10，100.068r11，110.048r12，120.083r13，130.068r14，140.047r15，150.553r16，160.932r17，170.932r18，180.405r19，190.645r20，200.952r21，210.952r22，220.595r2，1-1r3，1-1r4，1-1r5，1-1r6，1-1r7，10.083r8，10.553r9，10.645r10，10.932r11，10.952r12，10.917r13，10.932r14，10.952r15，10.447r16，10.068r17，10.068r18，10.595r19，10.355r20，10.048a21，10.048r22，10.405r7，2-1r8，2-1r9，2-1r10，3-1r11，3-1r12，4-1r13，4-1r14，4-1r15，5-1r16，5-1r17，5-1r18，5-1r19，6-1r20，6-1a21，6-1r22，6-1利用公式：式中t表示时间。计算得到此永磁电机系统的可靠度函数：由可靠度函数r(t)计算得出此永磁电机系统的平均失效前时间：从而实现电机系统可靠性定量评估。例如，对于一台三相永磁电机和三相逆变器组成的永磁电机系统，如图3所示。通过图2所示的永磁电机系统的马尔科夫状态转移图，建立系统的状态转移矩阵a，解得每个状态之间的影响系数，最后计算得到永磁电机系统的可靠度函数，可靠度函数随时间的变化曲线如图4所示。对可靠度函数r(t)积分得到三相永磁同步电机系统的平均失效前时间，从而实现可靠性的定量评估。平均故障前时间对应可靠度函数随时间变化曲线与横轴(时间轴)和纵轴(可靠度轴)组成的闭合平面的面积大小。本发明提出的一种基于马尔科夫模型电机系统可靠性定量评估方法适用与其他不同拓扑结构的电机系统。以上所述仅为本发明的具体实施例，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。当前第1页12