可靠性工程概论03

1、回顾复习回顾复习p维修度M()n对可修产品在发生故障或失效后，在规定的条件下和规定的时间（0, ）内完成修复的概率。p修复率()n修理时间已达到某个时刻但尚未修复的产品，在该时刻后的单位时间内完成修复的概率。p有效度A(t)n可维修产品在某时刻t具有或维持其功能的概率。第三章第三章p可修复系统的可靠性可修复系统的可靠性第三章第三章 可修复系统的可靠性可修复系统的可靠性3.1 马尔可夫过程3.2 状态转移图3.3 n步转移后系统各状态概率3.4 单部件可修系统3.5 串联可修系统3.6 并联可修系统引言引言 可修复系统的组成单元发生故障后，经过修理可以使系统恢复至正常工作状态，如下图所示。如果工

2、作时间和修复时间都服从指数分布，就可以借助马尔可夫过程来描述。 3.1 马尔可夫过程马尔可夫过程p马尔可夫过程定义马尔可夫过程定义 马尔可夫过程是一类“后效性”的随机过程。简单地说，在这种过程中系统将来的状态只与现在的状态有关，而与过去的状态无关。或者说，若已知系统在t0时刻所处的状态，那么t t0时的状态仅与时刻t0的状态有关。3.1 马尔可夫过程马尔可夫过程p马尔可夫过程的数学描述马尔可夫过程的数学描述 设x(t),t0是取值在E=0,1,2,或E=0,1,2,N上的一个随机过程。若对任意n个时刻点0t1t2tn 均有: Px(tn)=in|x(t1)=i1,x(t2)=i2,x(tn-1

3、)=in-1 =Px(tn)=in|x(tn-1)=in-1 i1,i2,inE 则称x(t),t0为离散状态空间E上连续时间马尔可夫过程。 3.1 马尔可夫过程马尔可夫过程p齐次马尔可夫过程齐次马尔可夫过程如果对任意t,u0,均有 Px(t+u)=j|x(u)=i=Pij(t) i,jE 与始点u 无关，则称该马尔可夫过程是齐次的。3.1 马尔可夫过程马尔可夫过程p转移矩阵转移矩阵 Pij(t)称为从状态i到状态j的转移函数，由转移函数的全体组成的矩阵称为转移矩阵。如对n个状态系统的转移矩阵为nn阶方阵，可写为： PPPPPPPPPPnnnnnn2122221112113.1 马尔可夫过程马

4、尔可夫过程p齐次马氏过程的性质齐次马氏过程的性质可以证明，对系统寿命以及故障后的修复时间均服从指数分布时，则系统状态变化的随机过程x(t),t0是一个齐次马尔可夫过程。 njijijPP11; 103.1 马尔可夫过程马尔可夫过程n三条假设三条假设a)，为常数(即寿命和维修时间服从指数分布)b)部件和系统取正常和故障两种状态。c)在相当小的t内，发生两个或两个以上部件同时进行状态转移的概率是t的高阶无穷小，此概率可以忽略不计。3.1 马尔可夫过程马尔可夫过程n可修复系统的可靠性特征量可修复系统的可靠性特征量a)瞬态可用度A(t)、不可用度Q(t);b)稳态可用度A、不可用度Q；c)MTBF、M

5、TTFF(首次故障前平均时间)、MTTR(平均修复时间) 。3.2 状态转移图状态转移图n例例1如一台机器，运行到某一时刻t时，可能的状态为： e1正常； e2故障。如机器处于e1状态的概率P11=4/5,则e1向e2转移的概率P12=1P11=1/5；反过程，如机器处于e2状态，经过一定时间的修复返回e1 状态的概率是3/5，P21=3/5(维修度M();则修不好仍处于e2状态的概率是P22=1P21=2/5.3.2 状态转移图状态转移图p由此可写出系统的转移矩阵为:n转移矩阵Pij也表示事件ei 发生的条件下，事件ej发生的条件概率：Pij=P(ej|ei) ； n矩阵 P:行是起始状态，

6、由小到大；列是到达状态，由小到大排列，建立P时应与转移图联系起来。3.2 状态转移图状态转移图p例例2对于一可修系统，失效率和修复率、为常数，试画出状态转移图： e1正常； e2故障。3.2 状态转移图状态转移图p由此可写出： 通常令t=1,则有p 由此可知，状态转移图是求解(写出)转移矩阵的基础。11P此时转移矩阵P也称为微系数矩阵3.3 n步转移后系统各状态概率步转移后系统各状态概率n设系统初始状态是 的概率 ，由切普曼柯尔莫哥洛夫方程， 可表示为： 式中n=k+l,vE(状态空间)n 此式为由状态i经n步转移到状态j的概率，等于由状态i先经k步转移到状态v，然后由状态v经l步转移到状态j

7、的概率(此处v也可理解为从i到j的通道)。jin步转移Pnij)(Pnij)(PPPPlvjvkivlkijnij)()()(3.3 n步转移后系统各状态概率步转移后系统各状态概率p上式中，若令k=1,l=1,由 可决定 ，即由全部一步转移概率可确定全部两步转移概率。若重复上述方法，就可由全部一步转移概率决定所有的转移概率。p若用矩阵表示n步转移概率，即 ，则有： 转移矩阵 )()(PPnijn),(EjiPijPij)2( )nklPP PP3.3 n步转移后系统各状态概率步转移后系统各状态概率p 一般地，可利用转移概率和系统的初始状态，求出任意转移后系统各状态的概率。公式如下： 式中 P1

8、步转移概率； n步转移概率；n转移步数(次数)；P(0)系统初始状态向量， P(0)= P1(0), P2(0) Pi(0)初始t=0时刻系统处于i状态的概率P(n)n步转移后系统所处状态向量，P(n)= P1(n), P2(n),Pi(n) n步转移后系统处于i状态的概率PPnPn)0()(Pn3.3 n步转移后系统各状态概率步转移后系统各状态概率n例：如下图，已知P(0)=P1(0), P2(0)=1, 0,求n=1,2,等各步(次)转移后系统各状态的概率。 图中e1正常； e2故障。3.3 n步转移后系统各状态概率步转移后系统各状态概率解：依次求得 n=1，n=2， n=3，n=5时的状

9、态矩阵 由此可知，随着n的递增，P1(n)、 P2(n)逐渐趋于稳定。稳定状态概率称为极限概率。 21/21/2(2)(0)(0)(1)1/21/20.450.552/53/5PPPP PPPP555.0445.0)3(P55555. 044445. 0)5(P5.05.05/35/22/12/101)0()1(PPP3.3 n步转移后系统各状态概率步转移后系统各状态概率p本例n时的极限概率为P1()=4/9, P2()=5/9,即n时， 将收敛于一个定概率矩阵，即(本例为)： p在实践中常会遇到这样的情况，不管系统的初始状态如何，在经历了一段工作时间后，便会处于相对稳定状态，在数学上称之为各

10、态历经或遍历性。所谓遍历过程就是系统处于稳定状态的概率与初始状态无关的随机过程。具有这种性质的状态转移矩阵称为遍历矩阵。Pn9/59/49/59/4Pn3.3 n步转移后系统各状态概率步转移后系统各状态概率p如果转移矩阵P经过n次相乘后，所得矩阵的全部元素都大于0,即 (i,jE),(注：常以此为判断马尔可夫链是否为各态历经的或是否存在极限概率)，则这样的转移矩阵都是遍历矩阵。遍历矩阵一定存在极限概率(或稳定状态)。p经过n步转移后的极限状态，就是过程的平稳状态，既然如此，即使再多转移一步，状态概率也不会有变化，这样可以求出平稳状态。0)(Pnij3.3 n步转移后系统各状态概率步转移后系统各

11、状态概率p设平稳状态概率为P(n)= P1, P2Pn, P为一步转移概率矩阵，则求平稳状态概率，只需求解以下方程： 或写成：)()(nPPnPPPPPPPPPPPPPPPPnnnnnnnn21212222111211213.3 n步转移后系统各状态概率步转移后系统各状态概率p 展开后得： (j=1,2,n) （n个方程只有n-1个是独立的，因此必须再加另一个独立方程。） 由此即可求出n个平稳状态概率。 njjijniijPPPP1113.3 n步转移后系统各状态概率步转移后系统各状态概率p 例：求如图所示系统的平稳状态概率。3.3 n步转移后系统各状态概率步转移后系统各状态概率p解：一步转移

12、矩阵为: 设P(n)= P0 P1，则 PPttttPP101011001010PPPP101010PPPPP0P13.4 单部件可修系统单部件可修系统p单部件系统是指一个单元组成的系统(或把整个系统当作一个单元来研究)，部件故障，则系统故障；部件正常，则系统正常。()0( )1tx tt系统状态时刻 系统正常时刻 系统故障3.4 单部件可修系统单部件可修系统p部件的失效率、修复率分别是常数、，则：nt时刻系统处于工作(正常工作)状态，在tt+t之间内发生故障的条件概率为t (即为 )nt时刻系统处于故障状态，在tt+t之间即t时间内修复好的条件概率为t (即为 )01P10P3.4 单部件可

13、修系统单部件可修系统p 单部件可修系统状态转移图3.4 单部件可修系统单部件可修系统p上图中:p同理:ttxttxPtPP0)(| ) 1()(0101条件概率ttxttxPtPP10)(| )0()(0000ttxttxPtPP1)(| )0()(1010ttxttxPtPP11)(| ) 1()(11113.4 单部件可修系统单部件可修系统p上图的转移概率矩阵为:)(111010tPttttP3.4 单部件可修系统单部件可修系统p令p下面研究如何求解 和p首先，利用全概率公式可求出 和 的表达式 1)()(0)()(10txPtPtxPtP、)(0tP)(1tP)(0ttP)(1ttP3.

14、4 单部件可修系统单部件可修系统)()()1 ()()()()(1)(1)(|0)(0)(0)(|0)(0)()(101100000ttPtPttPtPtPtPtxPtxttxPtxPtxttxPttxPttP 此即为 的计算公式 )(0ttP3.4 单部件可修系统单部件可修系统p由上式展开、移项、两边除以 p若令 取极限有： (1) t0t)()(t)() t(lim10000ttPtPtPtP)()()(100tPtPtP3.4 单部件可修系统单部件可修系统p同理可得： (2) p(1)、(2)联立即可求出 和 。 p(1)、(2)的联立方程称为状态方程)()1 ()()(101tPttt

15、PttP)()()(101tPtPtP)(0tP)(1tP3.4 单部件可修系统单部件可修系统p下边求解状态方程p对上述(1)、(2)两边取拉氏变换：)()()()()()(101100tPLtPLtPLtPLtPLtPL)()()0()()()()0()(10111000sPsPPssPsPsPPssP拉氏变换引入的变量引用ssPtPLPssPtPL)()()0()()(3.4 单部件可修系统单部件可修系统p假设t=0时系统为正常状态，即 ， 。代入上式1)0(0P0)0(1P)()()()()(1)(101100sPsPssPsPsPssP)()()(11)(10sssPssssssP3.

16、4 单部件可修系统单部件可修系统p拉氏反变换： ttetPetP)(1)(0)()(3.4 单部件可修系统单部件可修系统p由此瞬态有效度(可用度)：p稳态有效度：p平均有效度：(0 , t)()(0tPtA)(lim)(tAAttetdttAttAtt2)(20)()()(1)(3.4 单部件可修系统单部件可修系统p由上述可归纳出解可修系统有效度的方法步骤如下：p(1)画出系统的状态转移图p(2)写出转移矩阵 p(3)令 ，求出P(也称为转移矩阵)p(4)求状态方程系数矩阵ApA=P-I I为与P同阶的单位矩阵，A又称为转移率矩阵)( tP 1t3.4 单部件可修系统单部件可修系统p(5)写出

17、状态方程式式中 为各状态概率向量 为各状态概率导数向量p (6)求解状态方程通常要给定初始状态 ，且常用拉氏变换及反变换求解法。AtPtP)()()(tP)(,),(),()(10tPtPtPtPn)(tP)(,),(),()(10tPtPtPtPn)0(,),0(),0()0(10nPPPP3.4 单部件可修系统单部件可修系统p如上例： tttttP11)(11PIPAAtPtP)()()(),()(),(1010tPtPtPtP3.4 单部件可修系统单部件可修系统p得状态方程 与前述一致p以下即可用拉氏变换法等求解方程)()()()()()(101100tPtPtPtPtPtP3.5 串联

18、可修系统串联可修系统nn个相同单元组成的串联系统n每个单元： 、 为常数n两种状态：p状态0：n个单元全正常，系统正常状态p状态1：任一单元故障，系统故障状态n因为任一单元故障，系统即停止工作(不会出现两个及以上单元同时故障的情况)3.5 串联可修系统串联可修系统n个相同单元组成的串联系统状态转移图3.5 串联可修系统串联可修系统p用前述方法：tttntntP11)(11nnPnnIPA3.5 串联可修系统串联可修系统p状态方程：p初始条件：AtPtP)()(nntPtPtPtP)(),()(),(10101)0(0P0)0(1P3.5 串联可修系统串联可修系统p用拉氏变换与反变换可解出：tn

19、ennntPtA)(0)()(nA)(3.5 串联可修系统串联可修系统pn个不同单元组成的串联系统系统有n+1个状态：p状态0：n个单元均正常，系统正常状态p状态1：单元1故障，其余正常，系统故障p状态2：单元2故障，其余正常，系统故障 p状态n：单元n故障，其余正常，系统故障ii3.5 串联可修系统串联可修系统3.5 串联可修系统串联可修系统 )(Pt121110120111001001nninnnitPn 令3.5 串联可修系统串联可修系统pA=P-IAtPtP)()(nnnniinntPtPtPtPtPtP0000)(,),(),()(,),(),(1121110103.5 串联可修系统串联可修系统p给定初始条件，(用拉氏正、反变换)解此方程组即可求得： p(瞬态)有效度： p稳态有效度： )(,),(),()(10tPt