一种基于分解法建模的可靠性监测方法技术

本发明专利技术公开了一种基于分解法建模的可靠性监测方法。在状态空间分解法的基础上，建立单个部件的SPN模型，然后根据子系统的组成结构建立单个子系统的SPN模型；将子系统SPN模型集成为系统模型，得到系统模型的可达图及缩小的等效半马尔科夫率图，按图建立系统模型的状态方程并求解，根据求解结果得出机械系统的可靠度，并通过计算机显示，实现对系统可靠性的实时监测。分解法使得建模的难度和计算量减小，系统的可靠性监测也更加简单易行。在用马尔科夫模型对机械系统进行可靠性分析时，该方法能解决马尔科夫模型的状态空间爆炸问题，不仅能模拟系统中各个部分间的依赖关系及系统的动态特性，也使得设计者能够合理地选择部件以及了解部件的动态性能。

【技术实现步骤摘要】
一种基于分解法建模的可靠性监测方法
本专利技术涉及一种基于分解法建模的可靠性监测方法，它主要用于机械系统的可靠性监测，属于可靠性

技术介绍
对于可靠性的研究，国内外学者已做了大量工作，并提出了许多模型和方法，例如，可靠性框图、故障树分析法、故障模式及影响分析法、马尔科夫模型和Petri网模型等。随机Petri网(stochasticPetrinets，SPN)是一种用于离散异步并发系统的建模工具，它通过对实际问题构造Petri网并对Petri网进行分析，从而揭示出系统的动态特性等重要信息。随机Petri网为系统的性能监测提供了一条新的、有效的途径。库所中令牌的分布决定变迁的使能和激发，变迁的激发又将改变令牌的分布。以变迁激发导致令牌在库所间的流动，Petri网能够用于模拟系统的动态运行过程，反映系统实时的动态特性。“状态空间法”是一种基于解答空间的问题表示和求解方法，它是以状态和操作符为基础的。在利用状态空间图表示时，从某个初始状态开始，每次加一个操作符，递增地建立起操作符的试验序列，直到达到目标状态为止。由于状态空间法需要扩展过多的节点，容易出现“组合爆炸”，结合分解法有助于减少状态空间爆炸问题。状态空间的分解定理，齐次马尔科夫链的状态空间S可唯一地分解成有限或无限多个互不相交的状态子集的并，即：S＝D∪C1∪C2∪…其中，D是所有非常返状态构成的状态子集；Cn(n＝1,2,…)是所有常返状态构成的不可约闭集。马尔科夫模型广泛应用于相关系统及动态系统的模拟，然而，由于状态空间爆炸，当系统成规模扩大时，马尔科夫模型就显得低效，整个系统难以用单一的模式进行建模。
技术实现思路
本专利技术的目的在于提供一种基于分解法建模的可靠性监测方法。将状态空间分解法应用于马尔科夫模型的建立过程中，用得到的马尔科夫模型建立系统的状态方程，并通过计算机快速求解，根据状态方程的解及计算机终端来对系统的可靠性进行监测。采用状态空间分解法使得建模的难度和计算量减小，而且更易于操作。本专利技术采取的技术方案如下：步骤一：分析机械系统及其层级。根据机械系统的组成结构，利用分解法将系统分解成三个层次。第一层为系统级；第二层将系统分解为若干个子系统；第三层将每个子系统分解为若干个部件。步骤二：单个部件SPN模型的建立。只考虑部件“工作”或“失效”两种状态，根据组件的威布尔分布参数(形状参数β和尺度参数θ)能够得出模型的变迁率：其中，表示第i个部件由工作状态转为失效状态的变迁率。步骤三：将部件模型集成为子系统模型。将子系统模型转换为等效单网模型。根据子系统的组成结构及其部件的威布尔参数得出每个单网模型的等效转换率。步骤四：系统模型的建立。结合系统的组成结构及其动态特征如维修，冗余备份等，通过集成步骤三中所得的子系统单网模型得到系统级模型。步骤五：系统模型的分析与求解。根据系统的SPN模型能够得出所有的可达标识和相应的可达图，并将其转换为等效半马尔科夫率模型。根据得到的等效半马尔科夫率图列出系统的状态方程并通过Matlab软件在计算机上求解。步骤五1：可达标识的推导。将SPN模型当做系统模型以完成从初始状态到所有系统状态的推导。初始标识M0，是一个反映SPN模型中所有库所状态值的矩阵，当库所中无令牌时矩阵中元素为0，当库所中有令牌时元素为1。根据激发的顺序和方式，利用下式能够完成系统从初始标识到可达标识/状态的推导。Mj＝M0+CFT其中，M0为初始标识；Mj为可达标识，j＝1到2n-1是由于变迁激发顺序的不同而导致的不同的可达状态；C为邻接矩阵；FT为计数向量的激发变换。邻接矩阵C由系统SPN模型推导而来，其中的元素值代表模型中的联接方式，“-1”表示联接库所到变迁的弧，“0”表示库所与变迁之间没有弧。“1”表示联接变迁到库所的弧。在半马尔科夫模型中，系统的状态由部件(正常或失效)的组合方式来决定。矩阵Mj的最后一列表示该特定标识下系统是否正常或失效，并用于监测系统的可靠性。步骤五2：可达图及等效半马尔科夫率模型。根据步骤五1中得到的初始标识及所有可达标识能够画出可达图。依据变迁的激发顺序能够得知可达状态互相之间的联系。可达图转换为等效半马尔科夫率模型。从等效半马尔科夫率图中能够清楚地看到从前一状态转到后一状态的路线及变迁率。步骤五3：系统模型的状态方程：在马尔科夫模型中，状态方程为一阶微分方程；状态方程中的变迁率为常数；系统的状态方程如下所示：其中，为状态ii随时间推移的变化率；M为从状态ii转出的变迁率之和；N为从状态j转入状态ii的变迁率；pii(t)为系统在时间t时为状态ii的概率；步骤六：系统可靠性监测。根据上一步对系统模型的分析与求解结果，对系统的可靠性进行监测。根据删除吸收状态的修复率后的等效半马尔科夫率图及上面所列出的微分方程组的解，能够得出系统可靠性R(t)：其中pii(t)表示时间t时系统处于第ii个状态的概率。根据上式得到的系统的可靠性数据得出系统可靠度随时间的变化图。相对于现有技术，本专利技术的有益效果如下：现有技术中，考虑到大多数机械系统的复杂性，对机械系统的可靠性进行实时监测是非常困难的，而本专利技术创造性地将SPN模型、马尔科夫模型以及状态空间分解法结合起来应用于可靠性监测领域，通过该方法，减少了状态空间的产生，从而避免了状态空间爆炸问题，使得建模的难度和计算量减小。同时，在用马尔科夫模型对机械系统进行可靠性分析时，该方法能解决马尔科夫模型的状态空间爆炸问题，此方法不仅能模拟系统中各个部分间的依赖关系及系统的动态特性(常规RBD及FTA技术是无法做到的)，也使得设计者能够合理地选择部件以及了解部件的动态性能。附图说明图1是机械系统可靠性监测的实施步骤。图2是某型号泵送系统的SPN模型。图3是某型号泵送系统的可达图。图4是某型号泵送系统删除吸收状态的修复率后的等效半马尔科夫率图。具体实施方式下面结合附图对本专利技术作进一步的说明：见图1，本专利技术一种基于分解法建模的可靠性监测方法，该方法的具体实施步骤如下：步骤一：将系统分解为三个层次。在图2中，系统(层次1)被分解为若干个子系统(层次2)，每个子系统又被分解为若干个部件(层次3)。分解的过程有助于减少状态空间爆炸问题。例如：一个系统共有九个部件(每个部件的状态都用二进制来表示)，则系统的状态将有29(＝512)种，但在使用分解法的模型中，将系统分为两个子系统，其中一个有四个部件，另一个有五个部件，则分解模型只有24+25(＝48)种状态。可见，分解模型的状态数量远远低于单个模型的状态数量。步骤二：单个部件SPN模型的建立。只考虑部件“工作”或“失效”两种状态。失效变迁的激发按照变迁失效率进行，部件由“工作”转为“失效”的状态变化由库所中令牌的变化来显示，即由变为步骤三：单个子系统模型的建立。层次2中的子系统由若干个部件组成。按照子系统的组成结构(如串联、并联、串并联和混联等)来建立子系统模型，将SPN模型转换为等效的单网模型。机械系统的子系统的基本结构有串联、并联、串并联和混联。不同的结构对应的子系统模型也不同。步骤三中，串联子系统的SPN模型集成过程如下：假设一个串联子系统由三个部件组成，在其SPN模型中，有令牌的库所表示与其对应的状态，如表示第j个层次中第i个本文档来自技高网...

【技术保护点】
一种基于分解法建模的可靠性监测方法，其特征在于：该方法包括如下步骤：步骤一：分析机械系统及其层级：根据机械系统的组成结构，利用分解法将系统分解成三个层次；第一层为系统级；第二层将系统分解为若干个子系统；第三层将每个子系统分解为若干个部件；步骤二：单个部件SPN模型的建立：只考虑部件处于“工作”或“失效”两种状态，根据部件的威布尔分布参数得出模型的变迁率：λi3=(βiθi)(tθi)βi-1]]>其中，表示第i个部件由工作状态转为失效状态的变迁率；威布尔分布参数包括形状参数β和尺度参数θ；步骤三：将部件SPN模型集成为子系统模型：将子系统SPN模型转换为等效单网模型，根据子系统的组成结构及其部件的威布尔参数得出每个单网模型的等效转换率；子系统由若干个部件组成；按照子系统的组成结构来建立子系统模型，将SPN模型转换为等效的单网模型；子系统的组成结构包括串联、并联、串并联和混联；不同的结构对应的子系统模型也不同；步骤四：系统模型的建立：结合系统的组成结构及动态特征，通过集成步骤三中所得的子系统单网模型得到系统级模型；所述动态特征包括维修，冗余备份特征；SPN模型中增加显示系统状态的库所；通过相互之间引入虚拟变迁，使得状态库所与部件库所之间可达。步骤五：系统模型的分析与求解：按照步骤四得到的系统模型，推导出系统的所有可达标识和相应的可达图，并将其转换为等效半马尔科夫率模型；根据得到的等效半马尔科夫率图列出系统的状态方程并通过Matlab软件在计算机上求解；步骤六：系统可靠性监测：根据删除吸收状态的修复率后的等效半马尔科夫率图及系统模型的状态方程得出系统可靠性，对机械系统可靠性进行监测。...

【技术特征摘要】
1.一种基于分解法建模的可靠性监测方法，其特征在于：该方法包括如下步骤：步骤一：分析机械系统及其层级：根据机械系统的组成结构，利用分解法将系统分解成三个层次；第一层为系统级；第二层将系统分解为若干个子系统；第三层将每个子系统分解为若干个部件；步骤二：单个部件SPN模型的建立：只考虑部件处于“工作”或“失效”两种状态，根据部件的威布尔分布参数得出模型的变迁率：其中，表示第i个部件由工作状态转为失效状态的变迁率；威布尔分布参数包括形状参数β和尺度参数θ；步骤三：将部件SPN模型集成为子系统模型：将子系统SPN模型转换为等效单网模型，根据子系统的组成结构及其部件的威布尔参数得出每个单网模型的等效转换率；子系统由若干个部件组成；按照子系统的组成结构来建立子系统模型，将SPN模型转换为等效的单网模型；子系统的组成结构包括串联、并联、串并联和混联；不同的结构对应的子系统模型也不同；步骤四：系统模型的建立：结合系统的组成结构及动态特征，通过集成步骤三中所得的子系统单网模型得到系统级模型；所述动态特征包括维修，冗余备份特征；SPN模型中增加显示系统状态的库所；通过相互之间引入虚拟变迁，使得状态库所与部件库所之间可达；步骤五：系统模型的分析与求解：按照步骤四得到的系统模型，推导出系统的所有可达标识和相应的可达图，并将其转换为等效半马尔科夫率模型；根据得到的等效半马尔科夫率图列出系统的状态方程并通过Matlab软件在计算机上求解；步骤五具体包括如下步骤：步骤五1：可达标识的推导：将SPN模型当做系统模型以完成从初始状态到所有系统状态的推导；初始标识M0，是一个反映SPN模型中所有库所状态值的矩阵，当库所中无令牌时矩阵中元素为0，当库所中有令牌时元素为1；根据激发的顺序和方式，利用下式能够完成系统从初始标识到可达标识/状态的推导：Mj＝M0+CFT其中，M0为初始标识；Mj为可达标识，j＝1到2n-1是由于变迁激发顺序的不同而导致的不同的可达状态；C为邻接矩阵；FT为计数向量的激发变换；邻接矩阵C由系统SPN模型推导而来，其中的元素值代表模型中的联接方式，“-1”表示联接库所到变迁的弧，“0”表示库所与变迁之间没有弧；“1”表示联接变迁到库所的弧；在半马尔科夫模型中，系统的状态由部件正常或失效的组合方式来决定；矩阵Mj的最后一列表示该特定标识下系统是否正常或失效，并用于监测系统的可靠性；步骤五2：可达图及等效半马尔科夫率模型：根据步骤五1中得到的初始标识及所有可达标识能够画出可达图；依据变迁的激发顺序能够得知可达状态互相之间的联系；可达图转换为等效半马尔科夫率模型；从等效半马尔科夫率图中能够清楚地看到从前一状态转到后一状态的路线及变迁率；步骤五3：系统模型的状态方程：在马尔科夫模型中，状态方程为一阶微分方程；状态方程中的变迁率为常数；系统的状态方程如下所示：其中，为状态ii随时间推移的变化率；M为从状态ii转出的变迁率之和；N为从状态j转入状态ii的变迁率；pii(t)为系统在时间t时为状态ii的概率；步骤六：系统可靠性监测：根据删除吸收状态的修复率后的等效半马尔科夫率图及系统模型的状态方程得出系统可靠性，对机械系统可靠性进行监测。2.根据权利要求1所述的基于分解法建模的可靠性监测方法，其特征在于：步骤三中，串联子系统的SPN模型集成过程如下：假设一个串联子系统由三个部件组成，在其SPN模型中，有令牌的库所表示与其对应的状态，表示第j个层次中第i个部件的正常状态；表示第i个部件的失效变迁表示第i个部件失效变迁的变迁率按照变迁率来激发；符号中的参数“f”表示对应部件为失效状态；符号中参数“o”表示对应部件为正常状态；在该模型中，失效变迁率按威布尔分布；如果串联子系统中的任何一个部件的失效变迁被激发则整个系统将失效；然而，根据具体情况，子系统的失效变迁在瞬间发生或延迟发生；矩形代表延迟变迁，单实线代表直接或瞬发变迁；如果部件失效后系统立即随之失效，则称之为瞬发变迁；变迁率或失效率的威布尔参数用公式：

【专利技术属性】
技术研发人员：李俊，王贤琳，邱爽，王飞，
申请(专利权)人：武汉科技大学，
类型：发明
国别省市：湖北;42