CN106886667A - 一种基于事件调度的复杂***可用性分析方法 - Google Patents
一种基于事件调度的复杂***可用性分析方法 Download PDFInfo
- Publication number
- CN106886667A CN106886667A CN201710242847.4A CN201710242847A CN106886667A CN 106886667 A CN106886667 A CN 106886667A CN 201710242847 A CN201710242847 A CN 201710242847A CN 106886667 A CN106886667 A CN 106886667A
- Authority
- CN
- China
- Prior art keywords
- time
- event
- state
- failure
- availability
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Development Economics (AREA)
- Geometry (AREA)
- Educational Administration (AREA)
- Evolutionary Computation (AREA)
- Game Theory and Decision Science (AREA)
- Computer Hardware Design (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于事件调度的复杂***可用性分析方法。该方法在考虑部件和***故障/修理时间概率分布特性以及修理策略的基础上,以事件调度方法为仿真框架,通过建立部件解析模型和概率模型,建立***的可用仿真模型,由此来模拟具有多状态复杂***的可用性,该方法可以实现3种功能:a)不同目标区间内***可用性计算;b)确定关键部件对***可靠性的影响程度;c)准确统计***平均故障时间、***平均维修时间、故障次数/维修次数等参数。本发明可以实现多状态多部件复杂***的可用性计算,其计算过程逼近于***的真实过程,反映出***内各部件之间的关联和影响影响因素和概率特性,有很强的工程应用价值和商业推广价值。
Description
技术领域
本发明涉及一种可用性分析方法,具体涉及一种基于事件调度的复杂***可用性分析方法,属于复杂***可靠性技术领域。
背景技术
目前,复杂***的可靠性和可用性研究是工程领域、试验学以及社会科学研究的重点。研究方法主要有解析方法和仿真方法。解析方法主要采用随机过程理论构建***的解析模型,很难实现复杂***的可用度计算;仿真方法可以研究多部件多状态的***可用性,比较常用的是通过蒙特卡洛仿真方法,但是该方法在计算***可用度时,当***部件和状态较多时,容易出现空间***的问题。
另一方面,随着可用度的研究越来越受到重视,其商业价值也越来越大,国外研究人员开发了比较成熟的商业软件:Avsim(www.isograph-software.com)、Raptor(www.arinc.com/products/raptor/index.html)以及Blocksim(www.reliasoft.com/Blocksim),但是商业软件由于其昂贵的成本限制了其使用。
发明内容
本发明采用基于事件调度的方法,针对复杂***在运行过程中其状态变量是离散变量的实际情况,通过分析***结构的拓扑性和部件的故障特性,采用事件调度方法构建复杂***的可用度仿真模型,从而对***在目标时间内的可用性进行准确的分析,为实现上述目标,其技术方案如下:
步骤一,定义***和部件的状态变量
假设***由n(n≥2)个部件组成,具有k(k≥2)种状态;
假定初始时刻为t0,***状态向量为X(t0),第i个部件的状态向量为Xi(t0),如果第i个部件可用则Xi(t0)=1,反之Xi(t0)=0;
φ(X(t0))为***在t0时刻的状态函数,当***正常运行,所有部件完好时***的状态φ(X(t0))=2;当***正常运行,部分部件不完好时φ(X(t0))=1;***不能运行时φ(X(t0))=0;用ps(t0)表示***在t0时刻处于s状态的概率,即ps(t0)=Pr{φ(X(t0))=s},s=0,1,2;s=0表示***处于不能运行状态,s=1表示***处于正常运行,但部分部件不完好的状态,s=2表示***处于正常运行,所有部件完好的状态;
确定***状态为s的概率ps(t0)和置信区间;
步骤二,分析部件、***结构的拓扑特性和故障行为特性
(1)根据***的运行原理和功能特性,绘制***的可靠性框图;
(2)根据***的可靠性框图结构的拓扑特性,分析部件的独立性和关联性;
(3)根据故障树原理***的最小割集,确定可实现***功能的最短路;
步骤三,设计基于事件调度的仿真流程
该方法以事件调度方法为总体框架,在仿真过程中,考虑两种随机事件,分别为部件故障事件和修理事件,基本步骤如下:
(1)定义模拟变量:仿真时钟、初始化部件状态向量、初始化***状态向量、***的行为参数;
(2)利用已知可用信息,确定部件故障时间分布类型、参数和部件的修理时间;
(3)确定下一个发生的事件:比如部件失效或者部件修理,并且采集与该事件对应的相关信息:事件类型、事件发生的时间、会受到影响的其它部件;
(4)根据上述步骤提供的信息,更新全局时间变量和***状态变量;
(5)如果在仿真时间和下一个时间发生之前没有要观测的事件,则认为***对应的状态函数φ(X(t0))>0,否则为零;
(6)更新下一个事件的部件状态,如果为部件失效,则产生一个随机的部件修理时间,更新其余部件的故障和维修时间,分解详细的与之有关的故障事件;相反,如果下一个事件为部件修理,则产生一个随机的部件故障时间,更新其余部件的故障和维修时间,并且计算与之相关的修理时间;
(7)根据***的当前状态和结构框图的拓扑特征判断下一个事件是否会影响***的状态即根据判断***是否可用的结果,更新***状态变量和统计变量;
(8)更新仿真时钟;
(9)重复子步骤(3)~(8)直至满足截止条件,使之完成当前步骤的迭代,一次迭代可以使仿真时钟完成所有的目标时间并达到最终时间;
(10)记录当前的迭代变量,***在每个时间区域内的可用情况,***的MTTF和MTTR,目标时间区域内***的失效和修理次数;
(11)重复子步骤(1)~(10)多次,直至达到满意的观测效果。
优选地,步骤一中确定***状态为s的概率ps(t0)和置信区间的具体方法为:
(1)用随机变量Y(t0)=φ(X(t0))表示***在t0时刻的观测值,假设进行m次观测,其值分别为Y1(t0),Y2(t0),…,Ym(t0),根据大数定理
(2)根据中心极限定理,当m取值足够大时,置信度为1-α的ps(t0)(s=0,1,2),可表示为其中zα是在标准正态分布中置信水平为1-α/2的百分数,qs(t0)=1-ps(t0);
本发明能够实现多部件、多状态复杂***的可用性计算,其计算过程逼近于***的真实过程,反映出***内各部件之间的关联和影响影响因素和概率特性,主要能够实现以下功能:a)不同目标区间内***可用性;b)确定关键部件以及其对于提高***可靠性的影响程度;c)统计***平均故障时间、***平均维修时间、失效次数、修理次数和***利用率。
附图说明
图1为本发明技术方案示意图;
图2为***在生命周期内的状态变化图;
图3为机载通信网络的***结构图;
图4为基于事件调度的仿真流程图。
具体实施方式
以机载通信网络的可用性分析过程说明本方法的实施过程。
步骤(1):构建发明研究技术方案,如图1所示;
步骤(2):分析***的状态变化过程,如图2所示;
步骤(3):分析***的结构和故障行为特性;
机载通信由一个传输机、一个接收机和6个调制解调机组成,其组成结构如图3所示。
***中,传输机由SPC1-TRC1-TRC2三部分组成,接收机由SPS1-RCR1-RCR2三部分组成,每个调制解调机由SPS1和RLYC1组成;
步骤(4):统计***的故障信息和修理信息,其相关信息如表1所示:
表1部件故障信息和修理信息
步骤(5):根据***结构,确定***的最小路,如表2所示:
表2***的最小路
路径编号 | 包含的部件 | 包含的子部件 |
1 | Tr-Re1-Re3-Re5-Rc | 1-2-3-4-5-8-9-12-13-16-17-18 |
2 | Tr-Re1-Re3-Re4-Re6-Rc | 1-2-3-4-5-8-9-10-11-14-15-16-17-18 |
3 | Tr-Re2-Re3-Re5-Rc | 1-2-3-6-7-8-9-12-13-16-17-18 |
4 | Tr-Re2-Re4-Re5-Rc | 1-2-3-6-7-8-10-11-12-13-16-17-18 |
5 | Tr-Re1-Re3-Re6-Rc | 1-2-3-6-7-8-10-11-14-15-16-17-18 |
步骤(6):设计仿真基本流程,如图4所示,设定算法迭代步骤次数;
步骤(7):根据以上步骤,对***可用度进行评估;
在配置为4.0GRAM,interi3,3.10hz安装windowns7.0的计算机上,采用本发明设计算法,模拟***运行1000h,***不同可用度的置信区间间范围如表3所示:
表3***可用性评估
由结果分析可知,1)***的平均故障时间MTTF为8764h以99%的置信度在(8759,8769)区间内;2)关键部件识别,子部件1,2和16是三个最重要的部件,其在迭代过程引起***故障的概率分别为15.54%、15.51%和15.54%,前两个部件位于传输机中,第三个部件位于接收机中,可见提高其可靠性对提高***可靠性有重要意义。
Claims (2)
1.一种基于事件调度的复杂***可用性分析方法,其特征在于,包括以下步骤:
步骤一,定义***和部件的状态变量:
假设***由n(n≥2)个部件组成,具有k(k≥2)种状态;
假定初始时刻为t0,部件状态向量为X(t0),第i个部件的状态向量为Xi(t0),如果第i个部件可用,则Xi(t0)=1,反之Xi(t0)=0;
φ(X(t0))为***在t0时刻的状态函数,当***正常运行,所有部件完好时***的状态φ(X(t0))=2;当***正常运行,部分部件不完好时φ(X(t0))=1;***不能运行时φ(X(t0))=0;对于确定的目标时间t0>0,用ps(t0)表示***在t0时刻处于s状态的概率,即ps(t0)=Pr{φ(X(t0))=s},s=0,1,2;
确定***状态为s的概率ps(t0)和置信区间;
步骤二,分析部件、***结构的拓扑特性和故障行为特性:
(1)根据***的运行原理和功能特性,绘制***的可靠性框图;
(2)根据***可靠性框图结构的拓扑特性,分析部件的独立性和关联性;
(3)根据故障树原理确定***的最小割集,确定可实现***功能的最短路;
步骤三,设计基于事件调度的仿真流程:
(1)定义模拟变量:仿真时钟、初始化部件状态向量、初始化***状态向量、***的行为参数;
(2)利用已知可用信息,确定部件故障时间分布类型、参数和部件的修理时间;
(3)确定下一个发生的事件,部件失效或者部件修理,并且采集与该事件对应的相关信息:事件类型、事件发生的时间、会受到影响的其它部件;
(4)根据上述步骤提供的信息,更新全局时间变量和***可用状态;
(5)如果在仿真时间和下一个时间发生之前没有要观测的事件,则认为***对应的状态函数φ(X(t0))>0,否则为零;
(6)更新下一个事件的部件状态,如果为部件失效,则产生一个随机的部件修理时间,更新其余部件的故障和维修时间,计算与之有关的故障时间;相反,如果为部件修理,则产生一个随机的部件故障时间,更新其余部件的故障和维修时间,并且计算与之相关的修理时间;
(7)根据***的当前状态和可靠性框图的拓扑特征判断下一个事件是否会影响***的状态,即根据判断***是否可用的结果更新***状态变量和统计变量;
(8)更新仿真时钟;
(9)重复子步骤(3)~(8)直至满足截止条件,使之完成当前步骤的迭代,一次迭代可以使仿真时钟完全所有的目标时间并达到最终时间;
(10)记录当前的迭代变量,***在每个时间区域内的可用情况,***的MTTF和MTTR,在每个目标时间区域内的部件可用性,***的失效和修理次数;
(11)重复子步骤(1)~(10)多次直至达到满意的观测效果。
2.如权利要求1所述的一种基于事件调度的复杂***可用性分析方法,其特征在于,权利要求1步骤一中确定***状态为s的概率ps(t0)和置信区间的具体方法为:
(1)用随机变量Y(t0)=φ(X(t0))表示***在t0时刻的观测值,假设进行m次观测,其值分别为Y1(t0),Y2(t0),…,Ym(t0),根据大数定理
(2)根据中心极限定理,当m取值足够大时,置信度为1-α的ps(t0)(s=0,1,2)的置信区间为其中zα是在标准正态分布中置信水平为1-α/2的百分数,qs(t0)=1-ps(t0)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710242847.4A CN106886667A (zh) | 2017-04-14 | 2017-04-14 | 一种基于事件调度的复杂***可用性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710242847.4A CN106886667A (zh) | 2017-04-14 | 2017-04-14 | 一种基于事件调度的复杂***可用性分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106886667A true CN106886667A (zh) | 2017-06-23 |
Family
ID=59183867
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710242847.4A Pending CN106886667A (zh) | 2017-04-14 | 2017-04-14 | 一种基于事件调度的复杂***可用性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106886667A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109104304A (zh) * | 2018-07-24 | 2018-12-28 | 国网山东省电力公司电力科学研究院 | 一种分布式实时故障处理方法 |
CN109377039A (zh) * | 2018-10-12 | 2019-02-22 | 中国人民解放军92942部队 | 一种***任务可靠性关键故障因素分析方法 |
CN110263446A (zh) * | 2019-06-24 | 2019-09-20 | 广东工业大学 | 一种基于改进fmea的生产线可靠性改进分析方法、*** |
CN113919204A (zh) * | 2021-10-29 | 2022-01-11 | 长春工业大学 | 一种多状态制造***可用性综合重要度分析方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030182599A1 (en) * | 2002-03-21 | 2003-09-25 | Gray William M. | Method and system for assessing availability of complex electronic systems, including computer systems |
CN103592141A (zh) * | 2012-08-13 | 2014-02-19 | Uptime工程有限公司 | 用于测试复杂***的可靠性的方法 |
-
2017
- 2017-04-14 CN CN201710242847.4A patent/CN106886667A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030182599A1 (en) * | 2002-03-21 | 2003-09-25 | Gray William M. | Method and system for assessing availability of complex electronic systems, including computer systems |
CN103592141A (zh) * | 2012-08-13 | 2014-02-19 | Uptime工程有限公司 | 用于测试复杂***的可靠性的方法 |
Non-Patent Citations (2)
Title |
---|
JAVIER FAULIN 等: "Predicting availability functions in time-dependent complex systems with SAEDES simulation algorithms", 《ELSEVIER》 * |
陈志诚 等: "基于可靠性框图的可靠性建模研究", 《工程设计学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109104304A (zh) * | 2018-07-24 | 2018-12-28 | 国网山东省电力公司电力科学研究院 | 一种分布式实时故障处理方法 |
CN109104304B (zh) * | 2018-07-24 | 2021-06-01 | 国网山东省电力公司电力科学研究院 | 一种分布式实时故障处理方法 |
CN109377039A (zh) * | 2018-10-12 | 2019-02-22 | 中国人民解放军92942部队 | 一种***任务可靠性关键故障因素分析方法 |
CN110263446A (zh) * | 2019-06-24 | 2019-09-20 | 广东工业大学 | 一种基于改进fmea的生产线可靠性改进分析方法、*** |
CN113919204A (zh) * | 2021-10-29 | 2022-01-11 | 长春工业大学 | 一种多状态制造***可用性综合重要度分析方法 |
CN113919204B (zh) * | 2021-10-29 | 2024-04-30 | 长春工业大学 | 一种多状态制造***可用性综合重要度分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112202736B (zh) | 基于统计学习和深度学习的通信网络异常分类方法 | |
CN106886667A (zh) | 一种基于事件调度的复杂***可用性分析方法 | |
Agathokleous et al. | Topological robustness and vulnerability assessment of water distribution networks | |
US8577663B2 (en) | System and methods for fault-isolation and fault-mitigation based on network modeling | |
ATE450830T1 (de) | Sensor-fehlerdiagnose und -prognose unter verwendung eines komponentenmodells und zeitskalenorthogonaler entwicklungen | |
Bayer et al. | Beta regression control chart for monitoring fractions and proportions | |
CN105278460B (zh) | 基于级联故障分析的数控机床***组件可靠性评价方法 | |
CN114493049A (zh) | 基于数字孪生的产线优化方法、装置、电子设备及介质 | |
CN105808366B (zh) | 一种基于四变量模型的***安全分析方法 | |
CN102270302B (zh) | 一种基于灰色支持向量机的多应力加速寿命试验预测方法 | |
CN117196159A (zh) | 基于互联网大数据分析的智慧水务分区计量*** | |
Feizollahi et al. | The robust redundancy allocation problem in series-parallel systems with budgeted uncertainty | |
CN104657613A (zh) | 一种复杂机电***使用寿命评估方法 | |
EP3979120A2 (en) | Probabilistic fatigue and blend limit assessment and visualization methods for airfoils | |
CN108614443A (zh) | Phm***模型开发和验证平台设计方法 | |
CN103678858A (zh) | 一种存在竞争失效条件下的设备剩余寿命预测方法 | |
CN103729693A (zh) | 基于确定性库存退化模型的维修及备件供给联合优化方法 | |
CN109471698A (zh) | 云环境下虚拟机异常行为检测***和方法 | |
Qin et al. | Association analysis-based cybersecurity risk assessment for industrial control systems | |
Amari et al. | A fast and robust reliability evaluation algorithm for generalized multi-state $ k $-out-of-$ n $ systems | |
Zhu et al. | Early warning scheme of COVID-19 related internet public opinion based on RVM-L model | |
Hu et al. | Maintenance policy optimization for buffered serial systems considering energy-saving based on dual time windows | |
CN105490836A (zh) | 一种复杂网络可靠度的蒙特卡罗评估方法 | |
Friederich et al. | A Framework for Validating Data-Driven Discrete-Event Simulation Models of Cyber-Physical Production Systems | |
CN112560223B (zh) | 一种航空发动机全寿命周期维修概率建模与成本预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170623 |
|
WD01 | Invention patent application deemed withdrawn after publication |