基于聚类状态主控边界点的单调多态关联系统可靠性分析

张永进， 孙有朝， 张燕军

1.安徽工业大学 数理科学与工程学院， 马鞍山 243002 2.南京航空航天大学 民航学院， 南京 211106 3.扬州大学 机械工程学院， 扬州 225127

基于聚类状态主控边界点的单调多态关联系统可靠性分析

张永进1， 孙有朝2,*， 张燕军3

1.安徽工业大学 数理科学与工程学院， 马鞍山 243002 2.南京航空航天大学 民航学院， 南京 211106 3.扬州大学 机械工程学院， 扬州 225127

鉴于组成系统单元的多状态单调关联性特征，将多元离散函数理论引入描述系统状态结构函数，发展了控制状态等价类主导状态向量的状态等价类主控边界点的逻辑方法，推导了多态单调关联系统的状态结构函数、可靠性和期望状态表达式；面向顾客的需求偏好，将负效用函数嵌入系统平均性能效用模型；鉴于元件状态引起的计算复杂性，提出了集合运算的德摩根律法和新型的框图式算法，简化了系统可靠度的表达式。结合某型航空发动机的简化演算，验证了主导等价类向量方法和框图算法的合理性与有效性，为工程系统的可靠性设计和可靠性管理提供理论依据。

主控边界点； 离散函数； 多状态关联系统； 单调性； 可靠性

随着系统设计的结构越来越复杂化，组成单元从“完美工作状态”到“完全失效状态”的演化过程中往往会经历若干中间状态，而系统的状态取决于系统内部若干个元件提供的状态组合，从而系统表现出多状态特征，它是工程机械系统性能的重要表征，如锅炉燃烧的燃油量由若干喷油器喷出，汽轮机所需蒸汽量由若干锅炉提供等。随着多态系统研究的深入，工程上提出系统中任何一个元件都有其特定的功能，而且一个部件质量的改善不会降低系统的可靠性[1]，这就使得单调关联多态系统(Multi-State Coherent System：MSCS)的可靠性研究具有重要的现实意义。

最早Barlow[2]和EI-Neweihi[3]等对多态关联系统进行了研究，对多态关联系统给出了一系列基本定义，随后Natring[4]进一步完善了多态关联系统基本概念。目前，常见关于多态关联系统的可靠性研究方法有多状态故障树[5-6]、系统结构函数[7]、多状态最小路与割[8-13]、齐次马尔可夫模型[14-15]、非齐次马尔可夫模型[16]、随机衰退模型[17-18]、通用生成函数[19]和寿命分布函数[20]等方法，这些方法中多数并不考虑面向使用者的系统效能。进一步，系统可靠性评估通常建立在子单元可靠性数据的基础上，为此Barlow和Wu[2]将二态关联系统的可靠性理论推广到多态关联系统，EI-Neweihi[3]获得了系统性能与部件性能之间的关系，分别给出了一个并联系统和一个串联系统中存在的可靠性上下边界，Kundu[21]和Franko[22]等研究了单调关联系统可靠性问题。考虑部件性能改进对系统性能改进的影响，Griffith[23]给出了一个多态系统结构函数公理化方法，考虑系统效用，得到了系统性能度量方法。基于最小路和最小割原理，Butler[24]提出用模块分解方法来改进多态系统的可靠性边界，并研究了可靠性界限。为研究关联系统结构特性，Block和Savits[25]推导了多态系统结构函数的分解理论，给出了系统元件关联性与重要性基本概念。面向顾客偏好的多态关联系统，Hudson与Kapur[26]首次提出了等价类和边界点基本概念，分析了系统所处固定状态的概率，应用容斥原理和不交集对系统可靠性的状态边界进行了改进。基于顾客对系统元件状态变化如何影响系统状态的陈述，Boedigheimer和Kapur[7]建立了顾客偏好下系统的结构函数，给出了系统处于指定状态的概率计算方法。应用等价类和边界点结构函数方法，Lisnianski与Levitin[27]对系统可靠性性能进行了分析。鉴于结构函数描述系统状态的有效性，Liu[28]和Magana[29]等进一步完善了等价类结构函数评估系统可靠性的方法，建立了多态单调关联系统的结构函数，应用最小路和最小割方法研究了系统可靠性，然而随着单元数以及元件状态数的增加，计算复杂性也快速增加。

针对工程中常见的多状态单调关联系统，本文拟应用结构函数与离散函数相关理论，基于Lisnianski[27]、Liu[30]以及Magana[29]等给出的系统状态等价类及边界点方法，将一般状态等价类推广到可靠的状态聚类，应用状态等价聚类的边界点方法对多态关联系统可靠性进行研究，提出基于逻辑代数的德摩根律方法以及简易框图算法，优化计算复杂程度。最后给出一个算例分析，将航空发动机系统分解成一般性可靠性框图结构，分析该系统所处的各个状态时的可靠度，验证所给出的新方法的有效性，为工程实际的复杂多态关联系统可靠性设计提供理论基础。

1 单调多态关联系统

考虑一个含有n个元件的单调MSCS。设一个多态系统具有M+1个状态，不妨记为状态集{0,1,…,M}，设该系统由n个元件构成，而且这些元件具有多个状态，记为{0,1,…,mi}(i=1,2,…,n)，即元件共有mi+1个状态。若一个系统具有初始性能状态，并非完好状态时，不妨记这些状态为{l,l-1,…,l-k+1}⊂{0,1,…,M}，其中l-k+1为系统的最低可用状态，l为系统初始状态，k为任意状态。

由于一个系统所处状态总是由其组成元件的状态所确定，不妨记xi为第i个元件所处状态，即xi∈{0,1,…,mi}，于是系统所有元件状态构成一个状态向量x=[x1,x2,…,xn]，则一定存在一个系统状态s∈{0,1,…,M}及从部件状态到系统状态的映射φ，满足

φ(x)=s

φ:[0,1,…,m1]·…·[0,1,…,mn]→

[0,1,…,M]

称φ为系统状态的结构函数。

x=[x1,x2,…,xn]

ΩCi∈{0,1,…,mi}，ΩS∈{0,1,…,M}

(0,1,…,mn)

为方便构建多状态单调关联系统的结构函数φ(x)，首先给出状态向量序概念。

定义1设x=[x1,x2,…,xn]，y=[y1,y2,…,yn]是2个元件的状态向量。若单元状态向量满足“xi>yi且至少存在一个元件i满足xi>yi(i=1,2,…,n)”，则状态向量x优于状态向量y，记作x≻y。

下面给出MSCS的数学定义及性质[31]。

定义2对于一个由n个元件构成的多状态系统{C,φ}，其中C表示元件，记xi和yi是系统第i个元件所处的状态，则有

① 规范性(边界条件)。若系统所有部件失效，则系统处于完全失效状态；若系统所有部件处于完美状态，则系统处于完美状态，即有

φ(0,0,…,0)=0，φ(m1,m2,…,mn)=M

② 可达性。对于系统的每一个状态k，至少有一个部件状态向量x=[x1,x2,…,xn]满足

φ(x)=φ(x1,x2,…,xn)=k∈{0,1,…,M}。

③ 关联性。系统中每个元件均具有其特定作用，不存在无关元件，即系统中不存在不关联元件。若标记元件状态向量

(αi,x)=(x1,x2,…,xi-1,αi,xi+1,…,xn)

(βi,x)=(x1,x2,…,xi-1,βi,xi+1,…,xn)

则关联性应用数学方法可描述为：若αi≠βi，则φ(αi,x)≠φ(βi,x)。

④ 单调性。部件性能改进不会导致系统性能恶化，即系统元件的状态函数φ(x)是单调增的。若2个部件状态向量满足序x≻y，则系统状态向量满足φ(x)>φ(y)；若向量序x≤y，则系统状态向量满足φ(x)≤φ(y)。

通常，满足上述性质的系统{C,φ}称为多状态单调关联系统，结构函数φ(x)为单调关联结构函数。

2 可靠状态类与边界

考察系统所处状态的性能水平，需要分析对应系统性能水平下部件的状态组合，由于不同的部件状态组合可能使得系统性能处于同一个状态水平，从而需要找到同一个系统水平下各个不同元件所处状态的集合类，即状态等价类。不妨将部件状态向量空间划分为一系列等价类，每个等价类对应于系统的特定状态。

定义3对于系统的任意状态k，称集合

k∈{0,1,…,M}

为系统状态等价类。

不同顾客对象对于系统性能水平的要求不同，多态关联系统的可靠性可以理解为系统所处的状态不低于系统规定的性能状态。鉴于此，应用定义3从面向顾客要求的性能水平角度给出多态系统的可靠状态等价类定义4。

定义4对于规定的系统最低可靠性状态k，标记系统状态不低于k的元件状态集合表示为等价类

k∈{0,1,…,M}

由于给出的等价类分类中元素随着组成部件的增加会快速增加，为方便计算，下面给出主导每一个等价类中的边界元素定义。

根据状态等价类定义知，任一个可靠等价类可能包含不止一个状态组合，而且可能有些状态组合对系统可靠性评估计算时是冗余的，导致在系统结构函数的简化过程中会出现无关的元件状态向量，为此给出控制状态类的边界。

{φ(y)}max=k-1

基于已给出的假设和定义，多态单调关联系统在水平k时的可靠性可表示为

根据定义5给出的主控边界点定义，表明边界点的作用是用较少的等价类元素控制和主导等价类中其他元素，避免重复计算而使得计算量大大减少。

3 单调MSCS可靠性

要获取单调关联多态系统的结构函数，需要将二元逻辑运算扩展到多元逻辑运算，利用离散函数相关定义，首先给出逻辑运算的析取形式。

3.1 离散函数析取运算式

根据离散函数的析取运算，对于任意实数u和v，二元逻辑运算可以被延伸为下述析取运算形式：

u∪v=max(u,v)

为了将逻辑运算应用于多状态变量，需要将二元逻辑运算扩展到n元逻辑运算，下面给出二元变量集合运算{u1,u2,…,un}为

对于一个二元系统，考虑系统元件只有2个状态，若二元部件“成功”，则xi取值1，否则取值0，从而对应的指示函数I(xi=1)=xi，否则I(xi=0)=1-xi。于是根据枢轴量分解式，元件i的结构函数φ(x)可表示为“析取”范式形式：

φ(x)=xiφ(1i,x)+(1-xi)φ(0i,x)=

I(xi=1)φ(1i,x)+I(xi=0)φ(0i,x)=

(1)

式中：1i和0i表示第i个元素取值为1和0；若xi=α，则指示变量I(xi=α)=1，否则指示变量I(xi=α)=0。

从而根据式(1)，多状态系统的结构函数枢轴量分解可以表示为如下“析取”范式形式[26]：

(2)

式中：元件i的最低状态为0，最高状态为mi。显然，式(2)给出了单个元件的结构函数“析取”形式，于是对于一个由n个子单元构成的多态单调关联系统而言，根据定义5，其结构函数枢轴分解可以表示为“析取”范式形式[26]：

(3)

显然式(3)给出了一个多态系统结构函数表示形式，但随着状态数目的增加，结构函数的复杂程度会快速增加。

3.2 边界点结构函数

从式(3)知，若要计算一个单调关联多态系统结构函数，需要计算所有可能组合表达式，这是极其繁琐的，需要寻找新方法优化结构函数，为此将最小路集基本概念推广到MSCS，给出MSCS的最小路结构函数表示形式。

一个最小路集是一个确保系统正常工作的最少部件集合。一个系统通常有多个最小路集，不妨标记Pj是该系统的第j个最小路，j=1,2,…,p，p为系统的最小路集数，则一个二元关联系统结构函数可表示为

j=1,2,…,p

记最小路集中所有的部件正常取值1，否则取值为0。由于一个系统正常运行的充要条件是至少有一个最小路中所有部件正常运行，从而该系统的二元结构函数可表示为

(4)

(5)

若将各个下边界点的结构函数式(5)代入系统的结构函数式(4)，可以得到MSCS的结构函数形式为

(6)

式(6)为使用边界点的系统结构函数，通过使用第i个子单元状态与边界点中第i个元素进行逻辑比较，获取指示函数值。

为便于逻辑代数运算，可利用德摩根律，将式(6)重新表示和描述。首先不妨引入记号：

(x1x2…xn)

则式(6)可表示为

为求MSCS的可靠性，Lisnianski[27]和Liu[30]等给出了系统状态等价类的边界问题，以下简记为L-L方法,该方法的结构函数形式为

(7)

3.3 MSCS可靠度计算

对于一个MSCS，其可靠性可定义为在指定的时间和环境条件下，系统在一个指定或更高状态下执行任务的概率，从而MSCS以状态k或更高状态执行任务的可靠性可表示为

Rk=Pr{φ(x)≥k}=

X2=s2,…,Xn=sn}·

Pr(X1=s1,X2=s2,…,Xn=sn)=

Pr(X1=s1,X2=s2，…,Xn=sn)=

E{I[φ(x)≥k]}

(8)

式中：k=0,1,…,M，si∈{0,1,…,mi}，指示变量

若部件寿命具有统计独立性，则有

(9)

由式(8)可以观察到，计算Rk需要计算所有向量下的系统结构函数，随着系统状态的增加复杂程度快速增加。若应用式的最小析取范式，MSCS的可靠性计算会有所简化。

若给出运算标记

则上述可靠度计算式可重新分别简记为

(10)

显然，在式(10)的运算过程中，需要用到2个指示函数乘积问题。为计算Rk，首先给出其运算规律：

I(xi≥u)·I(xi≥v)=I(xi≥max(u,v))

(11)

式中：u和v是实数。基于规则式(11)，指示函数连乘的期望运算为

(12)

式中：pi,si为第i个元件处于状态si的概率。将式(11)、式(12)代入式(10)，系统处于不低于状态k下运行的可靠性Rk可通过式(13)进行计算：

(13)

而应用L-L方法，基于边界点计算系统处于状态k运行的可靠性表达式为

(14)

3.4 系统性能效用

对于系统设计者和使用者而言，多态系统处于某一个状态时的概率具有重要指导意义。根据式(13)和式(14)，系统处于某个固定状态k时的可靠性可表示为

Pr{φ(x)=k}=Rk-Rk+1

(15)

进一步，根据系统所处每一个状态的概率，系统所处平均状态可表示为

(16)

事实上，系统开始进入使用或贮存时的可靠性可表示为R0=Pr{φ(x)=M}=RM，从而若此多状态系统处于状态k时的可靠度包含有初始可靠性，那么系统仅仅因为本次使用或贮存而产生的可靠度可表示为

工程上，一个负效用函数可以衡量顾客偏好下的效用损失，经济损失与系统完美性能和当前状态差量密切相关，即若一个系统在最佳状态时，则顾客的负效用度达到最低。通常应用一个指数负效用函数去度量这种关系[29]：

DU[Δ(x)]=ρ[eβ(M-φ(x))-1]

式中：ρ、β为系统退化参数。当系统处于完美工作状态时，负效用值等于0；当系统完全失效时，系统负效用值为ρ[exp(βM)-1]。应用该函数可以计算出该系统的期望负效用值，具体为

Pr{M-φ(x)=k}=ρ·

(17)

事实上，式(17)的负效用能够刻画顾客角度的不满意程度，对于复杂系统可靠性设计具有一定的指导意义。

4 算例分析

下面给出一个某型号航空发动机系统的算例分析，考虑该整机系统有多个状态，其子系统也具有多个状态，应用本文提出的方法计算系统的可靠性。某型号双转子航空发动机系统结构如图1所示。

根据图1(d)中可靠性框图，假设各个单元状态列于表1～表4,x1、x2、x3和x4分别为系统各单元所处状态，根据这些子单元状态，验证所提出方法和原始方法的算法比较。

在航空发动机系统中，表征压气机性能好坏的主要参数有[32]：空气流量、增压比、效率和喘振裕度。由于增压比为压气机出口气流压强与其进口气流的压强之比，总增压比的高低，在设计时根据特定型号发动机的需要来确定，它是影响涡轮喷气发动机工作性能的一个重要的循环参数，对发动机的单位推力和耗油率有较大的影响。工程上，总增压比越高，发动机性能就越好(推重比大、耗油率低)。考虑某先进军用涡轮喷气发动机的总增压比为8～12，给出状态分类如表1所示。

由于发动机的燃烧室的功能是把燃料的化学能给释放出来，转换为热能，使气体的总焓增大，提高燃气的膨胀做功能力，其状态好坏，直接影响发动机的工作性能。而燃烧室主要由主燃烧室和加力燃烧室构成，主燃烧室位于压气机和涡轮之间，而加力燃烧室位于涡轮后喷管前，是为较短时间内增大发动机推力而设置的[32]。下面就燃烧室工作原理，将这一部分分成2个子单元(主燃烧单元B1，加力燃烧单元B2，如图1(b)和图1(d))所示，根据点火可靠、燃烧稳定与燃烧完全等划分为3个状态，具体数据见表2和表3。

图1 双转子航空发动机系统Fig.1 Twin-rotor aero engine system

涡轮的主要目的是将气流的能量转换为机械能，高温高压燃气在涡轮内膨胀，向外输出功带动压气机，气流流过时产生作用力，对转子叶片作功使其转动，而将气流的能量转换成机械能输出[32]。由于涡轮转子叶片是在高温燃气包围下工作，不仅承受自身气动力、离心力、热应力和振动应力等作用，而且要抵制高温燃气的冲刷、氧化与腐蚀作用，从而产生磨损、腐蚀和高温氧化等表面失效，以及蠕变、高低周疲劳断裂等整体失效模式，于是涡轮转子叶片的工作状况对航空发动机的运行性能与可靠性起着决定性作用。根据这些失效模式和机理，考虑涡轮转子具有4个不同状态。

表1 压气机单元A状态概率Table 1 State probability of compressor A

表2 主燃烧单元B1状态概率Table 2 State probability of main combustor B1

表3 加力燃烧单元B2状态概率Table 3 State probability of afterburner B2

表4 涡轮单元C状态概率Table 4 State probability of turbine C

根据航空发动机系统结构(如图1(a)所示)以及工作原理(图1(b)所示)，整个发动机系统的性能状态输出，考虑由图1(d)中框图结构形式输出。考虑整个航空发动机推动力输出符合下列公式

S=φ(a,b1,b2,c)=min{a,b1+b2,c}

(18)

式中：S为系统推动力输出量；a为压气机单元功能贡献量；b1为主燃烧室单元B1功能贡献量；b2为加力燃烧室单元B2功能贡献量；c为涡轮转子单元功能贡献量。贡献量为整体推力贡献的单位分量，主要是用来方便确定系统状态。

4.1 状态类与边界点

根据发动机系统结构状态情况的定量描述，以及式(18)、图1(d)和表1～表4相关数据。为方便，用(x1x2x3x4)表示系统状态输出，则系统状态等价类的计算结果展示于表5。

表5 系统输出的状态等价类Table 5 State equivalence Ω of system output

进一步，不妨将上述9个等价类按照贡献量输出情况进行适当合并，简化为5个等价类。依据等价类定义3，表5中给出了细分9个等价类时的边界点，显然从表5可以看出等价类中边界点个数要比等价类所有等价元素个数少很多，这也说明了边界点方法可以减少计算量。根据定义4、定义5，以及表5中数据，下面给出上等价类如表6 所示。

根据定义5以及表5、表6中相关数据，下面给出状态类下边界点以及上状态类的下边界点相关数据，如表7所示。

表6 系统输出的上等价类Table 6 Upper equivalence Ω of system output

表7 系统输出的下边界点Table 7 Lower boundary points of system output

根据以上计算得到的边界点数据，应用文中给出的分析方法计算多态系统的可靠度。

4.2 MSCS的结构函数运算

应用析取式(6)中的结构函数计算方法，下面给出基于表7中边界点数据的结构函数表达式。

基于表7中边界点及式(6)，给出主控边界下结构函数的框图形式如图2所示。

为便于应用，图2可表示为

(19)

若应用表7中状态类边界点，以及式(7)中的L-L方法，则析取形式的结构函数应为

图2 新方法框图形式Fig.2 Diagram form of new approach

(20)

为验证结构函数的有效性，下面给出单元状态向量[2,2,0,3]时系统所处状态计算。根据式(19)和式(20)有

事实上，根据表5中状态等价类的划分，单元状态向量[2,2,0,3]对应的系统状态为3。这表明基于状态等价类的主控边界点能够完全控制所有其他状态向量。

4.3 MSCS的可靠性计算

为演示新方法的有效性，考虑新框图算法和L-L方法，应用边界点和表7中相关数据，计算各状态下可靠度。下面仅给出状态1的可靠度计算过程，其他计算结果参见表8。

基于表7中上等价类的下边界点数据，应用图2的框图结构算法以及式(13)，系统处于状态1的可靠度应为

表8 系统可靠度计算结果Table 8 Results of system reliability

E[(1101)]+E[(1011)]-E[(1111)]=

(p1,1+p1,2)(p2,1+p2,2)·1·

(p4,1+p4,2+p4,3)+(p1,1+p1,2)·1·

(p3,1+p3,2)(p4,1+p4,2+p4,3)-

(p1,1+p1,2)(p2,1+p2,2)(p3,1+p3,2)·

(p4,1+p4,2+p4,3)=(0.04+0.94)·

(0.11+0.88)(0.06+0.13+0.80)+

(0.04+0.94)(0.08+0.91)·

(0.06+0.13+0.80)-(0.04+0.94)·

(0.11+0.88)(0.08+0.91)·

(0.06+0.13+0.80)=0.970 102 98

若应用表7中一般等价类的边界点数据及式(14)，L-L方法析取形式下系统处于状态1的可靠度应为

1-E{[1-(1101)][1-(1011)]·

[1-(1012)][1-(1202)]·

[1-(2202)][1-(2112)]·

[1-(2022)][1-(2212)]·

[1-(2123)][1-(2213)]}=

E[(1101)]+E[(1011)]-

E[(1111)]=0.970 102 98

显然析取形式下，应用L-L方法的计算工作量较大，而本文提出的方法较为简洁，减少计算复杂度。为便于分析比较，将上述计算结果列于表8。

从表8中数据知，提出的新方法与L-L方法计算结果均为一致。上述计算结果也表明，若整个发动机系统处于状态 4，即次完美状态，则压气机的工作状态必为完美状态，而涡轮至少处于低度磨损与腐蚀、轻度高温氧化状态，而主力和助力燃烧室至少一个子单元处于完美状态。

根据表8中计算结果，系统处于某个特定状态k的概率应为

Pr{φ(x)=0}=R0-R1=1-0.970 1=0.029 9

Pr{φ(x)=1}=R1-R2=0.970 1-0.910 3=

0.059 8

Pr{φ(x)=2}=R2-R3=0.910 3-0.872 5=

0.037 8

Pr{φ(x)=3}=R3-R4=0.872 5-0.857 1=

0.015 4

Pr{φ(x)=4}=R4-R5=0.857 1-0.730 4=

0.126 7

Pr{φ(x)=5}=R5=0.730 4

从而系统状态的期望性能水平为

即系统可靠性性能平均状态水平为4.340 3。

进一步，根据式(17)，不妨取c=100,β=0.5，基于顾客角度，该航空发动机系统的负效用水平应为

Pr{φ(x)=5-k}}=95.68

系统完全失效时的最大负效用为200.593 9，表明系统性能需要改进。

5 结 论

1) 鉴于工程复杂系统处于多个状态运行的现实问题，将多元离散函数理论引入结构函数描述系统状态，给出了单调关联系统的数学概念，状态等价类，控制等价类的主导向量等相关定义与逻辑方法。

2) 基于离散函数理论的析取与合取运算的代数规则，将二元逻辑运算推广到多元逻辑运算，给出了多状态系统的析取与合取范式，推导了系统多状态边界点的结构函数。为便于运算，提出了集合运算的德摩根律方法，简化了边界点下系统的结构函数和可靠度表达式。为评估整个系统所处的平均性能，给出了系统的期望状态计算式，考虑系统面向顾客的偏好与需求，建立了基于负效用函数的性能效用模型。

3) 为验证模型的有效性与合理性，将一个某型航空发动机系统进行框图分解，基于各个子单元对系统的贡献度，给出了系统等价类，获取了系统状态相关的上下边界点数据。为便于运算，算例中构建了新型的框图式算法，给出了析取与合取形式下的结构函数表示，大大减少了计算复杂度，增加了问题描述的清晰度。通过新框图算法和德摩根律方法，降低了可靠度的计算复杂程度。

4) 通过算例分析，基于本文提出的状态等价类的主控边界点数据，以L-L方法作为比较对象，提出的方法能够寻找到完全控制等价类中其它向量，而且计算量简洁许多，并不需要计算所有满足≥0状态的所有边界点，只需要计算控制其他边界点的状态等价类中主控边界点即可，大大减少计算复杂度。

5) 本文给出的状态等价类的主控边界点方法，推广了L-L方法，通过实例验证了提出方法的有效性与合理性。然而随着单元数和子单元状态数的增加，复杂度也快速增加，可以进一步考虑使用计算机，给出系统状态和可靠度的计算程序以便于应用。

[1] CHAHKANDI M, RUGGERI F, SUAREZ L A. A generalized signature of repairable coherent systems[J]. IEEE Transactions on Reliability, 2016, 65(1): 434-445.

[2] BARLOW R E, WU A S. Coherent systems with multi-state components[J]. Mathematics of Operations Research, 1978, 3(4): 275-281.

[3] EI-NEWEIHI E. Multistate coherent systems[J]. Journal of Applied Probability, 1978, 15(4): 675-688.

[4] NATRING B. Two suggestions of how to define a multistate coherent system[J]. Advances in Applied Probability, 1982, 14(2): 434-455.

[5] HUANG H Z, TONG X, ZUO M J. Posbist fault tree analysis of coherent systems[J]. Reliability Engineering & System Safety, 2004, 84(2): 141-148.

[6] AVEN T. On performance measures for multistate monotone systems[J]. Reliability Engineering & System Safety, 1993, 41(3): 259-266.

[7] BOEDIGHEIMER R A, KAPUR K C. Customer-driven reliability models for multistate coherent systems[J]. IEEE Transactions on Reliability, 1994, 43(1): 46-50.

[8] LUO T, TRIVEDI K S. An improved algorithm for coherent system reliability[J]. IEEE Transactions on Reliability, 1998, 47(1): 73-78.

[9] BOUTSIKAS M V, KOUTRAS M V. Generalized reliability bounds for coherent structures[J]. Journal of Applied Probability, 2000, 37(3): 778-794.

[10] ESARY J D, PROSCHAN F. A reliability bound for systems of maintained, interdependent components[J]. Journal of the American Statistical Association, 1970, 65(329): 329-338.

[11] HSIEH Y C. New reliability bounds for coherent systems[J]. Journal of the Operational Research Society, 2003, 54(9): 995-1001.

[12] LI J A, WU Y, LAI K K, et al. Reliability estimation of multi-state components and coherent systems[J]. Reliability Engineering & System Safety, 2005, 88(1): 93-98.

[13] XUE J, YANG K. Dynamic reliability analysis of coherent multistate systems[J]. IEEE Transactions on Reliability, 1995, 44(4): 683-688.

[14] CUI L, LI H. Analytical method for reliability and MTTF assessment of coherent systems with dependent components[J]. Reliability Engineering & System Safety, 2007, 92(3): 300-307.

[15] LIU Y W, KAPUR K C. Reliability measures for dynamic multistate nonrepairable systems and their applications to system performance evaluation[J]. IIE Transactions, 2006, 38(6): 511-520.

[16] NATVIG B. On the deterioration of non-repairable multistate strongly coherent systems[J]. Journal of Applied Probability, 2014, 51(51): 69-81.

[17] KUNDU P, HAZRA N K, NANDA A K. Reliability study of a coherent system with single general standby component[J]. Statistics & Probability Letters, 2016, 110: 25-33.

[18] ERYILMAZ S. A new look at dynamic behavior of binary coherent system from a state-level perspective[J]. Annals of Operations Research, 2014, 212(1): 115-125.

[19] ZHANG X, WILSON A. System reliability and component importance under dependence: A Copula approach[J]. Technometrics, 2016: 1-28.

[20] POLPO A, SINHA D, DE B. Nonparametric Bayesian estimation of reliabilities in a class of coherent systems[J]. IEEE Transactions on Reliability, 2013, 62(2): 455-465.

[21] KUNDU P, HAZRA N K, NANDA A K. Reliability study of a coherent system with single general standby component[J]. Statistics & Probability Letters, 2016, 110: 25-33.

[22] FRANKO C, OZHKUT M, KAN C. Reliability of coherent systems with a single cold standby component[J]. Journal of Computational & Applied Mathematics, 2015, 281(C): 230-238.

[23] GRIFFITH W S. Multistate reliability models[J]. Journal of Applied Probability, 1980, 17(3): 735-744.

[24] BUTLER D A. Bounding the reliability of multi-state systems[J]. Operations Research, 1982, 30(3): 530-544.

[25] BLOCK H W, SAVITS T H. A decomposition for multistate monotone systems[J]. Journal of Applied Probability, 1982, 19(2): 391-402.

[26] HUDSON J C, KAPUR K C. Reliability bounds for multi-state systems with multistate components[J]. Operations Research, 1985, 33(1): 153-160.

[27] LISNIANSKI A, LEVITIN G. Multi-state system reliability assessment, optimization, application[M]. Singapore: World Scientific, 2003: 89-153.

[28] LIU Y W, KAPUR K C. Customer’s cumulative experience measures for reliability of non-repairable aging multistate systems[J]. Quality Technology and Quantitative Management, 2007, 4(2): 225-234.

[29] MAGANA S, ALBERTO C. Dynamic reliability based performance measures for multi-state systems[D]. Seattle, WA: University of Washington, 2010: 116-130.

[30] LIU Y W. Multi-state system reliability: Models, dynamic measures and applications[D]. Seattle, WA: Univesity of Washington, 2006: 115-121.

[31] XUE J, YANG K. Dynamic reliability analysis of coherent multistate systems[J]. IEEE Transactions on Reliability, 1995, 44(4): 683-688.

[32] 刘长富, 邓明. 航空发动机结构分析[M]. 西安: 西北工业大学出版社, 2006: 72-112.

LIU C F, DENG M. Structural analysis of aeroengine[M]. Xi’an: Northwestern Polytechnical University Press, 2006: 72-112 (in Chinese).

(责任编辑： 张晗)

*Correspondingauthor.E-mail:sunyc@nuaa.edu.cn

Reliabilityanalysisformulti-statecoherentsystemwithmonotoniccomponentsbasedonpivotalboundarypointsofclusteringstates

ZHANGYongjin1，SUNYouchao2,*，ZHANGYanjun3

1.SchoolofMathematicsandPhysics,AnhuiUniversityofTechnology,Maanshan243002,China2.CollegeofCivilAviation,NanjingUniversityofAeronauticsandAstronautics,Nanjing211106,China3.CollegeofMechanicalEngineering,YangzhouUniversity,Yangzhou225127,China

Consideringthemonotoneandcoherenceofthemulti-statesystem,themultiplediscretefunctiontheoryisintroducedtodescribethestructurefunctionofsystemstate.Thelogicapproachesfortheequivalenceclassofthecomponentstatewhichcontrolthestatevectorofsystemareproposed,andtheexpressionsforthestatestructurefunction,reliabilityandexpectedstatesarederivedforthemulti-statecoherentsystem.Toavoidthecomplexityofcomputationcausedbythenumberofthestate,theDemogenlawandthenewblockdiagramalgorithmaredevelopedtosimplifytheexpressionforthesystemreliability.Anillustrativeexampleofacertaintypeofaeroengineverifiestheeffectivenessofthelogicvectormeasurecontrollingthestateequivalenceclassandtheblockdiagramalgorithm.Itprovidestheoreticalbasisforreliabilitydesignandreliabilitymanagementofsystemengineering.

pivotalboundarypoint;discretefunction;multi-statecoherentsystem;monotonicity;reliability

2016-10-19;Revised2017-02-26;Accepted2017-04-05;Publishedonline2017-04-191331

URL：www.cnki.net/kcms/detail/11.1929.V.20170419.1331.008.html

s：NationalNaturalScienceFoundationofChina(U1333119,60979019,60572171,51605424,71601002);CAACScienceandTechnologyProject(MHRD201123,MHRD200908,MHRD0722);NaturalScienceFoundationofJiangsuProvince(BK20150455).

2016-10-19；退修日期2017-02-26；录用日期2017-04-05； < class="emphasis_bold">网络出版时间

时间：2017-04-191331

www.cnki.net/kcms/detail/11.1929.V.20170419.1331.008.html

国家自然科学基金 (U1333119,60979019,60572171,51605424,71601002)； 中国民航局科技基金 (MHRD201123,MHRD200908,MHRD0722)； 江苏省自然科学基金 (BK20150455)

.E-mailsunyc@nuaa.edu.cn

张永进， 孙有朝， 张燕军． 基于聚类状态主控边界点的单调多态关联系统可靠性分析J． 航空学报,2017,38(8):220868.ZHANGYJ,SUNYC,ZHANGYJ.Reliabilityanalysisformulti-statecoherentsystemwithmonotoniccomponentsbasedonpivotalboundarypointsofclusteringstatesJ.ActaAeronauticaetAstronauticaSinica,2017,38(8):220868.

http://hkxb.buaa.edu.cnhkxb@buaa.edu.cn

10.7527/S1000-6893.2017.220868

V231; TB114.3； TH112

A

1000-6893(2017)08-220868-12