回顾复习维修度M(τ)对可修产品在发生故障或失效后,在规定的条件下和规定的时间(0,τ)内完成修复的概率。修复率μ(τ)修理时间已达到某个时刻但尚未修复的产品,在该时刻后的单位时间内完成修复的概率。有效度A(t)可维修产品在某时刻t具有或维持其功能的概率。第三章可修复系统的可靠性第三章可修复系统的可靠性3.1马尔可夫过程3.2状态转移图3.3n步转移后系统各状态概率3.4单部件可修系统3.5串联可修系统3.6并联可修系统引言可修复系统的组成单元发生故障后,经过修理可以使系统恢复至正常工作状态,如下图所示。如果工作时间和修复时间都服从指数分布,就可以借助马尔可夫过程来描述。3.1马尔可夫过程马尔可夫过程定义马尔可夫过程是一类“后效性”的随机过程。简单地说,在这种过程中系统将来的状态只与现在的状态有关,而与过去的状态无关。或者说,若已知系统在t0时刻所处的状态,那么tt0时的状态仅与时刻t0的状态有关。3.1马尔可夫过程马尔可夫过程的数学描述设{x(t),t≥0}是取值在E={0,1,2,…}或E={0,1,2,…,N}上的一个随机过程。若对任意n个时刻点0≤t1t2…tn均有:P{x(tn)=in|x(t1)=i1,x(t2)=i2,…,x(tn-1)=in-1}=P{x(tn)=in|x(tn-1)=in-1}i1,i2,…,in∈E则称{x(t),t≥0}为离散状态空间E上连续时间马尔可夫过程。3.1马尔可夫过程齐次马尔可夫过程如果对任意t,u≥0,均有P{x(t+u)=j|x(u)=i}=Pij(t)i,j∈E与始点u无关,则称该马尔可夫过程是齐次的。3.1马尔可夫过程转移矩阵Pij(t)称为从状态i到状态j的转移函数,由转移函数的全体组成的矩阵称为转移矩阵。如对n个状态系统的转移矩阵为n×n阶方阵,可写为:PPPPPPPPPPnnnnnn2122221112113.1马尔可夫过程齐次马氏过程的性质可以证明,对系统寿命以及故障后的修复时间均服从指数分布时,则系统状态变化的随机过程{x(t),t≥0}是一个齐次马尔可夫过程。njijijPP11;103.1马尔可夫过程三条假设a),为常数(即寿命和维修时间服从指数分布)b)部件和系统取正常和故障两种状态。c)在相当小的t内,发生两个或两个以上部件同时进行状态转移的概率是t的高阶无穷小,此概率可以忽略不计。3.1马尔可夫过程可修复系统的可靠性特征量a)瞬态可用度A(t)、不可用度Q(t);b)稳态可用度A、不可用度Q;c)MTBF、MTTFF(首次故障前平均时间)、MTTR(平均修复时间)。3.2状态转移图例1如一台机器,运行到某一时刻t时,可能的状态为:e1-正常;e2-故障。如机器处于e1状态的概率P11=4/5,则e1向e2转移的概率P12=1-P11=1/5;反过程,如机器处于e2状态,经过一定时间的修复返回e1状态的概率是3/5,P21=3/5(维修度M());则修不好仍处于e2状态的概率是P22=1-P21=2/5.3.2状态转移图由此可写出系统的转移矩阵为:转移矩阵Pij也表示事件ei发生的条件下,事件ej发生的条件概率:Pij=P(ej|ei);矩阵P:行是起始状态,由小到大;列是到达状态,由小到大排列,建立P时应与转移图联系起来。3.2状态转移图例2对于一可修系统,失效率和修复率λ、μ为常数,试画出状态转移图:e1——正常;e2——故障。3.2状态转移图由此可写出:通常令Δt=1,则有由此可知,状态转移图是求解(写出)转移矩阵的基础。11P此时转移矩阵P也称为微系数矩阵3.3n步转移后系统各状态概率设系统初始状态是的概率,由切普曼—柯尔莫哥洛夫方程,可表示为:式中n=k+l,vE(状态空间)此式为由状态i经n步转移到状态j的概率,等于由状态i先经k步转移到状态v,然后由状态v经l步转移到状态j的概率(此处v也可理解为从i到j的通道)。jin步转移Pnij)(Pnij)(PPPPlvjvkivlkijnij)()()(3.3n步转移后系统各状态概率上式中,若令k=1,l=1,由可决定,即由全部一步转移概率可确定全部两步转移概率。若重复上述方法,就可由全部一步转移概率决定所有的转移概率。若用矩阵表示n步转移概率,即,则有:转移矩阵][)()(PPnijn),(EjiPijPij)2(()nklPPPP3.3n步转移后系统各状态概率一般地,可利用转移概率和系统的初始状态,求出任意转移后系统各状态的概率。公式如下:式中P-1步转移概率;-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.3n步转移后系统各状态概率例:如下图,已知P(0)=[P1(0),P2(0)]=[1,0],求n=1,2,…等各步(次)转移后系统各状态的概率。图中e1——正常;e2——故障。3.3n步转移后系统各状态概率解:依次求得n=1,n=2,n=3,n=5时的状态矩阵由此可知,随着n的递增,P1(n)、P2(n)逐渐趋于稳定。稳定状态概率称为极限概率。21/21/2(2)(0)(0)(1)1/21/20.450.552/53/5PPPPPPPP555.0445.0)3(P55555.044445.0)5(P5.05.05/35/22/12/101)0()1(PPP3.3n步转移后系统各状态概率本例n时的极限概率为P1()=4/9,P2()=5/9,即n时,将收敛于一个定概率矩阵,即(本例为):在实践中常会遇到这样的情况,不管系统的初始状态如何,在经历了一段工作时间后,便会处于相对稳定状态,在数学上称之为各态历经或遍历性。所谓遍历过程就是系统处于稳定状态的概率与初始状态无关的随机过程。具有这种性质的状态转移矩阵称为遍历矩阵。Pn9/59/49/59/4Pn3.3n步转移后系统各状态概率如果转移矩阵P经过n次相乘后,所得矩阵的全部元素都大于0,即(i,jE),(注:常以此为判断马尔可夫链是否为各态历经的或是否存在极限概率),则这样的转移矩阵都是遍历矩阵。遍历矩阵一定存在极限概率(或稳定状态)。经过n步转移后的极限状态,就是过程的平稳状态,既然如此,即使再多转移一步,状态概率也不会有变化,这样可以求出平稳状态。0)(Pnij3.3n步转移后系统各状态概率设平稳状态概率为P(n)=[P1,P2…Pn],P为一步转移概率矩阵,则求平稳状态概率,只需求解以下方程:或写成:)()(nPPnPPPPPPPPPPPPPPPPnnnnnnnn21212222111211213.3n步转移后系统各状态概率展开后得:(j=1,2,…n)(n个方程只有n-1个是独立的,因此必须再加另一个独立方程。)由此即可求出n个平稳状态概率。njjijniijPPPP1113.3n步转移后系统各状态概率例:求如图所示系统的平稳状态概率。3.3n步转移后系统各状态概率解:一步转移矩阵为:设P(n)=[P0P1],则PPttttPP101011001010PPPP101010PPPPP0P13.4单部件可修系统单部件系统是指一个单元组成的系统(或把整个系统当作一个单元来研究),部件故障,则系统故障;部件正常,则系统正常。()0()1txtt系统状态时刻系统正常时刻系统故障3.4单部件可修系统部件的失效率、修复率分别是常数λ、μ,则:t时刻系统处于工作(正常工作)状态,在t→t+Δt之间内发生故障的条件概率为λΔt(即为)t时刻系统处于故障状态,在t→t+Δt之间即Δt时间内修复好的条件概率为μΔt(即为)01P10P3.4单部件可修系统单部件可修系统状态转移图3.4单部件可修系统上图中:同理:ttxttxPtPP}0)(|)1({)(0101条件概率ttxttxPtPP1}0)(|)0({)(0000ttxttxPtPP}1)(|)0({)(1010ttxttxPtPP1}1)(|)1({)(11113.4单部件可修系统上图的转移概率矩阵为:))((111010tPttttP3.4单部件可修系统令下面研究如何求解和首先,利用全概率公式可求出和的表达式}1)({)(}0)({)(10txPtPtxPtP、)(0tP)(1tP)(0ttP)(1ttP3.4单部件可修系统)()()1()()()()(}1)({}1)(|0)({}0)({}0)(|0)({}0)({)(101100000ttPtPttPtPtPtPtxPtxttxPtxPtxttxPttxPttP——此即为的计算公式)(0ttP3.4单部件可修系统由上式展开、移项、两边除以若令取极限有:(1)t0t)()(t)()t(lim10000ttPtPtPtP)()()(100tPtPtP3.4单部件可修系统同理可得:(2)(1)、(2)联立即可求出和。(1)、(2)的联立方程称为状态方程)()1()()(101tPtttPttP)()()(101tPtPtP)(0tP)(1tP3.4单部件可修系统下边求解状态方程对上述(1)、(2)两边取拉氏变换:)]([)]([)]([)]([)]([)]([101100tPLtPLtPLtPLtPLtPL)()()0()()()()0()(10111000sPsPPssPsPsPPssP拉氏变换引入的变量—引用ssPtPLPssPtPL)()]([)0()()]([3.4单部件可修系统假设t=0时系统为正常状态,即,。代入上式1)0(0P0)0(1P)()()()()(1)(101100sPsPssPsPsPssP)()()(11)(10sssPssssssP3.4单部件可修系统拉氏反变换:ttetPetP)(1)(0)()(3.4单部件可修系统由此瞬态有效度(可用度):稳态有效度:平均有效度:(0,t))()(0tPtA)(lim)(tAAttetdttAttAtt2)(20)()()(1)(3.4单部件可修系统由上述可归纳出解可修系统有效度的方法步骤如下:(1)画出系统的状态转移图(2)写出转移矩阵(3)令,求出P(也称为转移矩阵)(4)求状态方程系数矩阵AA=P-II为与P同阶的单位矩阵,A又称为转移率矩阵)(tP1t3.4单部件可修系统(5)写出状态方程式式中为各状态概率向量为各状态概率导数向量(6)求解状态方程通常要给定初始状态,且常用拉氏变换及反变换求解法。AtPtP)()()(tP)](,),(),([)(10tPtPtPtPn
本文标题:可靠性工程概论03
链接地址:https://www.777doc.com/doc-174164 .html