CN111046556A - 考虑维修的含间歇性工作逻辑门的动态故障树仿真方法 - Google Patents

考虑维修的含间歇性工作逻辑门的动态故障树仿真方法 Download PDF

Info

Publication number
CN111046556A
CN111046556A CN201911268889.0A CN201911268889A CN111046556A CN 111046556 A CN111046556 A CN 111046556A CN 201911268889 A CN201911268889 A CN 201911268889A CN 111046556 A CN111046556 A CN 111046556A
Authority
CN
China
Prior art keywords
event
time
simulation
unit
state
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.)
Granted
Application number
CN201911268889.0A
Other languages
English (en)
Other versions
CN111046556B (zh
Inventor
赵广燕
刘陆杰
孙宇锋
胡薇薇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN201911268889.0A priority Critical patent/CN111046556B/zh
Publication of CN111046556A publication Critical patent/CN111046556A/zh
Application granted granted Critical
Publication of CN111046556B publication Critical patent/CN111046556B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/20—Administration of product repair or maintenance

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Test And Diagnosis Of Digital Computers (AREA)

Abstract

一种考虑维修的含间歇性工作逻辑门的动态故障树仿真方法，即考虑维修情况下，含有间歇性工作逻辑门的故障树的仿真方法，找出间歇性工作系统在可维修情况下的薄弱环节，无论对于改进系统的设计，提高系统的可靠性水平，减轻安全问题发生，减少经济损失，都具有极为重要的意义。

Description

考虑维修的含间歇性工作逻辑门的动态故障树仿真方法
技术领域
本发明涉及系统可靠性及故障仿真的技术领域，尤其涉及一种考虑维修的含间歇性工作逻辑门的系统故障树的仿真方法，即含有该间歇性工作逻辑门的系统故障树模型在考虑维修情况下的仿真方法。
背景技术
随着重大工程及复杂装备的研发，系统的功能和组成越发庞大复杂，系统的工作环境也越发严酷。很多系统在运行过程中一旦发生故障将会导致巨大的经济损失甚至具有严重的安全隐患，系统故障后及时维修可以极大提高系统的可用度。
具体地，对于某核电系统来说，该系统具有较强的辐射性、组成单元众多，故障逻辑也较为复杂。该核电系统与其他系统相比，不仅在组成和规模上较为特殊，还主要体现在以下两个方面：1)系统中存在系统呈间歇性工作特点；2)系统中单元故障需要进行维修的特点。
故障树分析方法是工程应用中适用性最强的可靠性建模分析方法。
FTA(Fault Tree Analysis，故障树)模型根据系统由上至下的结构功能对应关系，用与、或、非等逻辑门来表示单元与单元、单元与系统之间的故障因果逻辑关系。对于组成单元较少、单调关联系统的可靠性分析，都很容易建立其故障树模型，进而进行定性、定量分析，但对于单元数量较多，单元间耦合关系复杂的系统，FTA模型普遍存在问题是：无法准确描述系统的可靠性逻辑关系，即无法建立系统的FTA模型。
动态故障树模型是在故障树模型的基础上增加了几个动态逻辑门，扩展了故障树的适用范围，动态故障树模型不仅可以用于建立常见的串联、并联和表决等静态逻辑结构的故障树模型，还可以建立单元故障按照一定顺序发生等动态逻辑结构的故障树模型。常用的动态故障树的逻辑门符号包括顺序相关门、优先与门、功能触发门和备份门，备份门包括冷备门、温备门和热备门。
在现有的故障树逻辑门计算中，静态逻辑门如或门、与门以及表决门等可以通过结构函数，计算其输出事件的发生概率，而动态逻辑门如功能触发门、顺序相关门以及备份门需要转化为马尔科夫链计算。对于具有间歇性工作逻辑的系统来说，现有的动态故障树逻辑门无法对其间歇性的工作逻辑进行描述，同时，对于含有该逻辑的系统可靠性指标也无法通过解析方法进行求解，而考虑维修的情况更加增加了计算的复杂性。
针对这一问题，本文提出了基于蒙特卡罗仿真的解决方法，提出了间歇性工作逻辑门的仿真方法，并给出了在考虑维修情况下的含间歇性工作逻辑门的系统故障树模型仿真方法，来找出含间歇性工作逻辑可维修系统的薄弱环节，无论对于改进系统的设计，提高系统的可靠性水平，减轻安全问题发生，减少经济损失，都具有极为重要的意义。
发明内容
为克服现有技术的缺陷，本发明要解决的技术问题是提供了一种考虑维修的含间歇性工作逻辑门的动态故障树的仿真方法，即考虑维修情况下，含有间歇性工作逻辑门的故障树的仿真方法，找出间歇性工作系统在可维修情况下的薄弱环节，无论对于改进系统的设计，提高系统的可靠性水平，减轻安全问题发生，减少经济损失，都具有极为重要的意义。
本发明的技术方案是：这种含间歇性工作逻辑门的系统故障树可维修时的仿真方法，
间歇性工作逻辑门包括：
输入事件：主单元(E1)和辅助单元(E2)的故障模式；
输出事件(Y)：是输入事件发生后，通过逻辑门的输出结果；
第一时间(T1)：表示主单元(E1)规定的工作时间；
第二时间(T2)：表示辅助单元(E2)规定的工作时间；
事件状态：包括正常、故障、停机三个状态，分别用0、1和2表示；其中，正常状态是指单元处于正常工作状态，即故障事件不发生且单元处于正常状态；故障状态是指单元发生了故障，不具备工作能力，即故障事件发生；停机状态是指单元工作一定时间之后，由于一些原因使得该单元暂时处于不可用状态，但是单元未发生故障，属于无故障停机状态，即故障事件不发生且单元处于无故障停机状态；该间歇性工作逻辑门模型表示：当主单元输入事件在工作期间发生故障之后或者达到规定的工作时间T1之后，输出事件为故障状态或者停机状态；当辅助单元输入事件为故障状态时，若现有输出事件为故障状态，则输出事件状态不变，若现有输出事件为停机状态，则输出事件变为故障状态；
该考虑维修的含间歇性工作逻辑门的动态故障树仿真方法包括以下步骤：
(1)仿真初始化：获得仿真输入信息，变量初始化，进行失效时间抽样，并对全部事件进行排序；
(2)获得最先发生事件及时间；
(3)判断最先发生事件类型：若是故障发生事件，则转到步骤(4)；若是修复事件，则转到步骤(5)；若是时间切换事件，则转到步骤(6)；
(4)故障事件处理：更改因故障事件发生影响的底事件状态，更改仿真中因故障事件发生受影响的时间切换事件及修复事件发生的时间，跳转步骤(9)；
(5)修复事件处理：更改因修改事件发生影响的底事件状态，更改仿真中因修复事件发生受影响的事件发生的时间，跳转步骤(9)；
(6)判断切换事件的属性：若是主单元切换事件，跳转步骤(7)；辅助单元切换事件，跳转步骤(8)；
(7)主单元切换事件处理：更改因主单元切换事件影响的底事件状态，更改仿真中因主单元切换事件影响的事件发生的时间，跳转步骤(9)；
(8)辅助单元切换事件处理：更改因辅助单元切换事件受影响的底事件状态，更改仿真中因辅助单元切换事件影响的事件发生的时间，跳转步骤(9)；
(9)调用相关逻辑门仿真算法，通扫故障树，输出顶事件状态；
(10)判断仿真时间SimuTime是否大于设定仿真时间SimAllTime，若否，则本次仿真还未结束，事件队列重新排序，转到步骤(2)继续仿真；若是，则本次仿真结束，继续步骤(11)；
(11)判断仿真次数SimNumber是否大于设定仿真次数SimAllNumber，若否，则仿真未结束，仿真次数加1，转到步骤(1)，重新进行失效时间抽样及排序，继续仿真；若是，则继续步骤(12)；
(12)统计仿真数据；
(13)仿真结束。
本发明利用蒙特卡罗抽样方法，结合基于事件的仿真思想，编写了含有间歇性工作逻辑门的动态故障树在可维修情况下的仿真算法，因此能够找出含间歇性工作单元的可维修系统的薄弱环节，无论对于改进系统的设计，提高系统的可靠性水平，减轻安全问题发生，减少经济损失，都具有极为重要的意义。
附图说明
图1是根据本发明的含间歇性工作逻辑门的系统故障树可维修时的仿真方法的流程图。
图2是间歇性工作逻辑门仿真方法的流程图。
图3是核电系统A的功能结构图。
图4是核电系统A的故障树。
图5是冷备份门仿真方法的流程图。
图6是可用度频率-概率图。
图7是MTBCF频率-概率图。
图8是MTTR频率-概率图。
图9是首次故障时间频率-概率图。
具体实施方式
间歇性工作逻辑门包括：
输入事件：主单元(E1)和辅助单元(E2)的故障模式；
输出事件(Y)：是输入事件发生后，通过逻辑门的输出结果；
第一时间(T1)：表示主单元(E1)规定的工作时间；
第二时间(T2)：表示辅助单元(E2)规定的工作时间；
事件状态：包括正常、故障、停机三个状态，分别用0、1和2表示；其中，正常状态是指单元处于正常工作状态，即故障事件不发生且单元处于正常状态；故障状态是指单元发生了故障，不具备工作能力，即故障事件发生；停机状态是指单元工作一定时间之后，由于一些原因使得该单元暂时处于不可用状态，但是单元未发生故障，属于无故障停机状态，即故障事件不发生且单元处于无故障停机状态；该间歇性工作逻辑门模型表示：当主单元输入事件在工作期间发生故障之后或者达到规定的工作时间T1之后，输出事件为故障状态或者停机状态；当辅助单元输入事件为故障状态时，若现有输出事件为故障状态，则输出事件状态不变，若现有输出事件为停机状态，则输出事件变为故障状态；
如图1所示，这种考虑维修情况下含间歇性工作逻辑门的动态故障树的仿真方法，包括以下步骤：
(1)仿真初始化：获得仿真输入信息，变量初始化，进行失效时间抽样，并对全部事件进行排序；
(2)获得最先发生事件及时间；
(3)判断最先发生事件类型：若是故障发生事件，则转到步骤(4)；若是修复事件，则转到步骤(5)；若是时间切换事件，则转到步骤(6)；
(4)故障事件处理：更改因故障事件发生影响的底事件状态，更改仿真中因故障事件发生受影响的时间切换事件及修复事件发生的时间，跳转步骤(9)；
(5)修复事件处理：更改因修改事件发生影响的底事件状态，更改仿真中因修复事件发生受影响的事件发生的时间，跳转步骤(9)；
(6)判断切换事件的属性：若是主单元切换事件，跳转步骤(7)；辅助单元切换事件，跳转步骤(8)；
(7)主单元切换事件处理：更改因主单元切换事件影响的底事件状态，更改仿真中因主单元切换事件影响的事件发生的时间，跳转步骤(9)；
(8)辅助单元切换事件处理：更改因辅助单元切换事件受影响的底事件状态，更改仿真中因辅助单元切换事件影响的事件发生的时间，跳转步骤(9)；
(9)调用相关逻辑门仿真算法，通扫故障树，输出顶事件状态；
(10)判断仿真时间SimuTime是否大于设定仿真时间SimAllTime，若否，则本次仿真还未结束，事件队列重新排序，转到步骤(2)继续仿真；若是，则本次仿真结束，继续步骤(11)；
(11)判断仿真次数SimNumber是否大于设定仿真次数SimAllNumber，若否，则仿真未结束，仿真次数加1，转到步骤(1)，重新进行失效时间抽样及排序，继续仿真；若是，则继续步骤(12)；
(12)统计仿真数据；
(13)仿真结束。
本发明利用蒙特卡罗抽样方法，结合基于事件的仿真思想，编写了含有间歇性工作逻辑门多动态故障树在可维修情况下的仿真算法，因此能够找出含间歇性工作单元的可维修系统的薄弱环节，无论对于改进系统的设计，提高系统的可靠性水平，减轻安全问题发生，减少经济损失，都具有极为重要的意义。
优选地，所述步骤(1)包括以下分步骤：
(1.1)输入仿真所需信息：输入所有底事件的失效分布函数、维修分布函数；设定仿真总时间SimAllTime、设定仿真总次数SimAllNumber；仿真时间用SimTime表示，仿真次数用SimNumber表示，事件队列用一个数组MyEvent()表示；设置主单元的规定运行时间为T1，辅助单元的规定运行时间为T2；设置所有底事件的初始状态，底事件状态表示：0代表正常，1代表故障，2代表停机等待；
(1.2)初始仿真次数SimNumber赋值为1：SimNumber＝1；
(1.3)系统仿真时间初始化：仿真时间SimTime赋值为0，SimTime＝0；主单元切换时间MainSwitchTime赋值为0，MainSwitchTime＝0；辅助单元切换时间AuxiliarySwitchTime赋值为0，AuxiliarySwitchTime＝0；单元失效时刻FailureTime赋值为0，FailureTime＝0；单元修复时刻RepaireTime赋值为0，RepaireTime＝0；
(1.4)失效时间抽样：利用蒙特卡罗抽样方法，抽取所有输入事件的失效时间；令主单元切换时间等于主单元的规定工作时间，即MainSwitchTime＝T1；辅助单元切换时间等于辅助单元的规定工作时间，即AuxiliarySwitchTime＝MainSwitchTime+T2；将抽取的失效时间和主、辅单元的切换时间放入事件队列MyEvent()中，并对事件进行排序。
优选地，所述步骤(2)中，从排序后的事件队列MyEvent中，获得最先发生事件的时间MyEvent(0)，仿真时间步进到该事件发生的时刻，即SimTime＝MyEvent(0)，将该事件从事件队列中移除。
优选地，所述步骤(4)包括以下分步骤：
(4.1)输入事件状态＝1；
(4.2)抽取该输入事件下次维修需要的时间；
(4.3)输入事件的下次修复时刻RepaireTime＝仿真时间SimTime+抽取的下次维修时间，并放入事件队列；
(4.4)判断发生故障单元的类型：若辅助单元故障，跳转步骤(4.5)；若主单元故障，跳转步骤(4.6)；若普通单元故障，跳转步骤(9)；
(4.5)下次主单元切换时间＝下次主单元切换事件发生时间+抽取的维修时间；下次辅助单元切换时间＝下次辅助单元切换事件发生时间+抽取的维修时间；并重新放入事件队列进行更新；转到步骤(9)；
(4.6)辅助单元状态＝0；下次辅助单元切换时间＝仿真时间+规定工作时间T2，并重新放入事件队列进行更新，转到步骤(13)；
(4.7)判断发生故障的主单元维修时间是否小于T2：若是，下次主单元切换时间＝仿真时间+规定工作时间T2+规定工作时间T1；若否，下次主单元切换时间＝仿真时间+抽取的维修时间+规定工作时间T1，转到步骤(9)。
优选地，所述步骤(5)包括以下分步骤：
(5.1)输入事件状态＝0；
(5.2)抽取该输入事件下一次失效时间；
(5.3)该输入事件下一次失效时刻FailureTime＝仿真时间SimTime+抽取的下一次失效时间，并重新放入事件队列进行更新，转到步骤(9)。
优选地，所述步骤(7)包括以下分步骤：
(7.1)令切换事件的主单元＝2，切换事件的辅助单元＝0；
(7.2)切换事件的辅助单元下次切换时刻AuxiliarySwitchTime＝仿真时间SimTime+规定工作时间T2，并放入事件队列进行更新；
(7.3)切换事件的主单元下次失效的时刻＝下次该主单元失效事件发生时间+规定工作时间T2，并放入事件队列进行更新；跳转步骤(9)。
优选地，所述步骤(8)包括以下分步骤：
(8.1)判断主单元是否发生故障：若主单元发生故障，切换事件的主单元＝1，切换事件的辅助单元＝2，切换事件的辅助单元下次失效的时刻＝下次改辅助单元失效事件发生时间+(下次主单元切换事件-当前仿真时间)，并更新仿真队列，跳转步骤(9)；若主单元未发生故障，切换事件的辅助单元＝2，切换事件的主单元＝0，跳转步骤(8.2)；
(8.2)切换事件的主单元下次切换时刻＝仿真时间+规定工作时间T1，并放入事件队列进行更新，跳转步骤(8.3)；
(8.3)切换事件的辅助单元下次失效的时刻＝下次该辅助单元失效事件发生时间+规定工作时间T1，并放入事件队列进行更新，跳转步骤(9)。
优选地，所述步骤(9)包括以下分步骤：
(9.1)根据建立的故障树模型，将故障事件作为子节点，搜寻连接的故障树逻辑门，并调用与门、或门、间歇性工作逻辑门，获得输出事件的状态；
(9.2)判断输出事件的状态是否为正常状态：若是，则转到步骤(9.4)；若否，则输出事件为故障状态，转到步骤(9.3)仿真；
(9.3)判断步骤(9.2)中得到的输出事件是否为顶事件：若是，则不再进行故障树搜寻，转到步骤(9.4)；若否，则转到步骤(9.1)仿真；
(9.4)判断故障树模型其他子树中是否存在该故障事件：若是，则转到步骤(9.1)进行故障树搜寻；若否，转到步骤(10)。
优选地，如图3所示，所述步骤(9.1)包括以下分步骤：
(9.1.1)获得输入事件的状态及时间阶段：主事件状态和辅助事件的状态，间歇性工作逻辑门的输入事件均有三个状态：正常、故障和停机，分别用0、1和2表示；
(9.1.2)判断当前工作时间阶段：属于主单元工作时间阶段T1，跳转步骤(9.1.3)；属于辅助单元工作时间阶段T2，跳转步骤(9.1.6)；
(9.1.3)判断输入事件主单元状态：若输入事件主单元状态为0，跳转步骤(9.1.4)；若输入事件主单元状态为1，跳转步骤(9.1.5)；
(9.1.4)若辅助单元状态为2，输出0，跳转步骤(9.1.9)；否则，仿真输入错误，跳转步骤(9.1.9)
(9.1.5)若辅助单元状态为2，输出1，跳转步骤(9.1.9)；否则，仿真输入错误，跳转步骤(9.1.9)；
(9.1.6)若主单元状态2，跳转步骤(9.1.7)；若主单元状态1，跳转步骤(9.1.8)；否则，仿真输入错误，跳转步骤(9.1.9)；
(9.1.7)判断辅助单元状态：若辅助单元状态为0，则输出2，跳转步骤(9.1.9)；若辅助单元状态为1，则输出1，跳转步骤(9.1.9)；否则，仿真输入错误，跳转步骤(9.1.9)；
(9.1.8)若辅助单元状态为0，则输出1，跳转步骤(9.1.9)；否则，仿真输入错误，跳转步骤(9.1.9)；
(9.1.9)输出事件的状态，间歇性工作逻辑门仿真结束。
优选地，所述步骤(12)中，仿真数据包括：系统或者其他单元的工作时间，故障次数，维修时间，然后计算相应的可靠性指标，评价系统或者单元的可靠性。
有关含有间歇性工作逻辑特点的系统，以图3中某核电系统加以说明，该系统包括第一子系统A1和第二子系统A2两部分。第一子系统A1和第二子系统A2都是间歇性工作子系统，系统中的主单元和辅助单元按照一定时间交替性工作，其中的单元均包含正常、故障和停机三种状态。
对于核电系统A来说，核电系统A包括A1和A2两个完全一致的子系统，A2作为A1的冷备份。以A1为例说明其组成，A2组成与A1一致。A1包括第一主单元A11和第一辅助单元A12两个组成部分，A1的工作模式为A11工作200h，停止工作，切换到A12工作150h。对于A1而言，A11在200h内发生故障或者工作时间达到200h后，A11切换到第一辅助单元A12工作，此时A1状态为故障(前者)或者停机(后者)，第一辅助单元A12无故障工作时间150h之后，A1恢复正常。核电系统A的功能结构图如图3所示。
根据间歇性工作逻辑门以及现有的动态故障树模型，建立核电系统A的故障树模型。故障树模型中，a表示核电系统A的故障模式，a1表示第一子系统A1的故障模式，a11表示第一主单元A11的故障模式，a12表示第一辅助单元A12的故障模式；a2表示第二子系统A2的故障模式，a21表示第二主单元A21的故障模式，a22表示第二辅助单元A22的故障模式。其故障树模型如图4所示。
在上述逻辑门中，存在两个逻辑门，冷备份门和不同单元呈交替性工作逻辑门，其中，含有间歇性工作逻辑门系统模型的仿真方法如上所述，冷备份仿真方法需要重新建立。
图5冷备份门的仿真方法步骤如下：
步骤1：获取输入事件的状态，其中，0表示正常状态，1表示故障状态，2表示停机状态，11表示备份状态，事件的状态用EventState表示。
步骤2：判断工作件是否处于正常状态？若是，则输出事件处于正常状态，转到步骤4；若否，工作件则由正常状态变为了无法正常工作的状态，转到步骤3继续仿真。
步骤3：判断备份件的状态是否处于正常状态？若是，则启动备份件，备份件由正常停机状态变为正常状态，备份件的状态值：EventState(备份件)＝0，转到步骤4；若否，备份件的状态为故障状态，即，EventState(备份件)＝1，输出事件的状态为故障状态，即，EventState(输出事件)＝1。
步骤4：输出事件的状态为在正常状态，故障不会再传递下去，输出事件的状态值EventState(输出事件)＝0。
步骤5：得到输出事件的状态值EventState，退出冷备份门。
已知核电系统A中组成单元A11、A12、A21和A22的失效时间和维修时间均服从指数分布，其对应的仿真输入数据如表1所示。
表1
Figure BDA0002313615070000121
编写核电系统A的仿真程序，仿真输入为：每个单元的相关参数以及仿真次数1000次、仿真时间8760h。仿真计算核电系统的可用性指标。
统计核电系统的相关可靠性指标，在总仿真次数1000次下，核电系统可用度、平均严重故障间隔时间(MTBCF)、平均任务修复时间和首次故障时间的频率-概率分布图分别如图6、7、8、9所示。
通过仿真程序计算核电系统的可靠性相关指标结果如表2所示。
表2
Figure BDA0002313615070000131
以上所述，仅是本发明的较佳实施例，并非对本发明作任何形式上的限制，凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属本发明技术方案的保护范围。

Claims (10)

1.考虑维修的含间歇性工作逻辑门的动态故障树的仿真方法，
间歇性工作逻辑门包括：
输入事件：主单元(E1)和辅助单元(E2)的故障模式；
输出事件(Y)：是输入事件发生后，通过逻辑门的输出结果；
第一时间(T1)：表示主单元(E1)规定的工作时间；
第二时间(T2)：表示辅助单元(E2)规定的工作时间；
事件状态：包括正常、故障、停机三个状态，分别用0、1和2表示；其中，正常状态是指单元处于正常工作状态，即故障事件不发生且单元处于正常状态；故障状态是指单元发生了故障，不具备工作能力，即故障事件发生；停机状态是指单元工作一定时间之后，由于一些原因使得该单元暂时处于不可用状态，但是单元未发生故障，属于无故障停机状态，即故障事件不发生且单元处于无故障停机状态；该间歇性工作逻辑门模型表示：当主单元输入事件在工作期间发生故障之后或者达到规定的工作时间T1之后，输出事件为故障状态或者停机状态；当辅助单元输入事件为故障状态时，若现有输出事件为故障状态，则输出事件状态不变，若现有输出事件为停机状态，则输出事件变为故障状态；
其特征在于：该考虑维修的含间歇性工作逻辑门的动态故障树的仿真方法包括以下步骤：
(1)仿真初始化：获得仿真输入信息，变量初始化，进行失效时间抽样，并对全部事件进行排序；
(2)获得最先发生事件及时间；
(3)判断最先发生事件类型：若是故障发生事件，则转到步骤(4)；若是修复事件，则转到步骤(5)；若是时间切换事件，则转到步骤(6)；
(4)故障事件处理：更改因故障事件发生影响的底事件状态，更改仿真中因故障事件发生受影响的时间切换事件及修复事件发生的时间，跳转步骤(9)；
(5)修复事件处理：更改因修改事件发生影响的底事件状态，更改仿真中因修复事件发生受影响的事件发生的时间，跳转步骤(9)；
(6)判断切换事件的属性：若是主单元切换事件，跳转步骤(7)；辅助单元切换事件，跳转步骤(8)；
(7)主单元切换事件处理：更改因主单元切换事件影响的底事件状态，更改仿真中因主单元切换事件影响的事件发生的时间，跳转步骤(9)；
(8)辅助单元切换事件处理：更改因辅助单元切换事件受影响的底事件状态，更改仿真中因辅助单元切换事件影响的事件发生的时间，跳转步骤(9)；
(9)调用相关逻辑门仿真算法，通扫故障树，输出顶事件状态；
(10)判断仿真时间SimuTime是否大于设定仿真时间SimAllTime，若否，则本次仿真还未结束，事件队列重新排序，转到步骤(2)继续仿真；若是，则本次仿真结束，继续步骤(11)；
(11)判断仿真次数SimNumber是否大于设定仿真次数SimAllNumber，若否，则仿真未结束，仿真次数加1，转到步骤(1)，重新进行失效时间抽样及排序，继续仿真；若是，则继续步骤(12)；
(12)统计仿真数据；
(13)仿真结束。
2.根据权利要求1所述的考虑维修的含间歇性工作逻辑门的动态故障树的仿真方法，其特征在于：所述步骤(1)包括以下分步骤：
(1.1)输入仿真所需信息：输入所有底事件的失效分布函数、维修分布函数；设定仿真总时间SimAllTime、设定仿真总次数SimAllNumber；仿真时间用SimTime表示，仿真次数用SimNumber表示，事件队列用一个数组MyEvent()表示；设置主单元的规定运行时间为T1，辅助单元的规定运行时间为T2；设置所有底事件的初始状态，底事件状态表示：0代表正常，1代表故障，2代表停机等待；
(1.2)初始仿真次数SimNumber赋值为1：SimNumber＝1；
(1.3)系统仿真时间初始化：仿真时间SimTime赋值为0，SimTime＝0；主单元切换时间MainSwitchTime赋值为0，MainSwitchTime＝0；辅助单元切换时间AuxiliarySwitchTime赋值为0，AuxiliarySwitchTime＝0；单元失效时刻FailureTime赋值为0，FailureTime＝0；单元修复时刻RepaireTime赋值为0，RepaireTime＝0；
(1.4)失效时间抽样：利用蒙特卡罗抽样方法，抽取所有输入事件的失效时间；令主单元切换时间等于主单元的规定工作时间，即MainSwitchTime＝T1；辅助单元切换时间等于辅助单元的规定工作时间，即AuxiliarySwitchTime＝MainSwitchTime+T2；将抽取的失效时间和主、辅单元的切换时间放入事件队列MyEvent()中，并对事件进行排序。
3.根据权利要求2所述的考虑维修的含间歇性工作逻辑门的动态故障树的仿真方法，其特征在于：所述步骤(2)中，从排序后的事件队列MyEvent中，获得最先发生事件的时间MyEvent(0)，仿真时间步进到该事件发生的时刻，即SimTime＝MyEvent(0)，将该事件从事件队列中移除。
4.根据权利要求3所述的考虑维修的含间歇性工作逻辑门的动态故障树的仿真方法，其特征在于：所述步骤(4)包括以下分步骤：
(4.1)输入事件状态＝1；
(4.2)抽取该输入事件下次维修需要的时间；
(4.3)输入事件的下次修复时刻RepaireTime＝仿真时间SimTime+抽取的下次维修时间，并放入事件队列；
(4.4)判断发生故障单元的类型：若辅助单元故障，跳转步骤(4.5)；
若主单元故障，跳转步骤(4.6)；若普通单元故障，跳转步骤(9)；
(4.5)下次主单元切换时间＝下次主单元切换事件发生时间+抽取的维修时间；下次辅助单元切换时间＝下次辅助单元切换事件发生时间+抽取的维修时间；并重新放入事件队列进行更新；转到步骤(9)；
(4.6)辅助单元状态＝0；下次辅助单元切换时间＝仿真时间+规定工作时间T2，并重新放入事件队列进行更新，转到步骤(13)；
(4.7)判断发生故障的主单元维修时间是否小于T2：若是，下次主单元切换时间＝仿真时间+规定工作时间T2+规定工作时间T1；若否，下次主单元切换时间＝仿真时间+抽取的维修时间+规定工作时间T1，转到步骤(9)。
5.根据权利要求4所述的考虑维修的含间歇性工作逻辑门的动态故障树的仿真方法，其特征在于：所述步骤(5)包括以下分步骤：
(5.1)输入事件状态＝0；
(5.2)抽取该输入事件下一次失效时间；
(5.3)该输入事件下一次失效时刻FailureTime＝仿真时间SimTime+抽取的下一次失效时间，并重新放入事件队列进行更新，转到步骤(9)。
6.根据权利要求5所述的考虑维修的含间歇性工作逻辑门的动态故障树的仿真方法，其特征在于：所述步骤(7)包括以下分步骤：
(7.1)令切换事件的主单元＝2，切换事件的辅助单元＝0；
(7.2)切换事件的辅助单元下次切换时刻AuxiliarySwitchTime＝仿真时间SimTime+规定工作时间T2，并放入事件队列进行更新；
(7.3)切换事件的主单元下次失效的时刻＝下次该主单元失效事件发生时间+规定工作时间T2，并放入事件队列进行更新；跳转步骤(9)。
7.根据权利要求6所述的考虑维修的含间歇性工作逻辑门的动态故障树仿真方法，其特征在于：所述步骤(8)包括以下分步骤：
(8.1)判断主单元是否发生故障：若主单元发生故障，切换事件的主单元＝1，切换事件的辅助单元＝2，切换事件的辅助单元下次失效的时刻＝下次改辅助单元失效事件发生时间+(下次主单元切换事件-当前仿真时间)，并更新仿真队列，跳转步骤(9)；若主单元未发生故障，切换事件的辅助单元＝2，切换事件的主单元＝0，跳转步骤
(8.2)；
(8.2)切换事件的主单元下次切换时刻＝仿真时间+规定工作时间T1，并放入事件队列进行更新，跳转步骤(8.3)；
(8.3)切换事件的辅助单元下次失效的时刻＝下次该辅助单元失效事件发生时间+规定工作时间T1，并放入事件队列进行更新，跳转步骤(9)。
8.根据权利要求7所述的考虑维修的含间歇性工作逻辑门的动态故障树仿真方法，其特征在于：所述步骤(9)包括以下分步骤：
(9.1)根据建立的故障树模型，将故障事件作为子节点，搜寻连接的故障树逻辑门，并调用与门、或门、间歇性工作逻辑门，获得输出事件的状态；
(9.2)判断输出事件的状态是否为正常状态：若是，则转到步骤(9.4)；若否，则输出事件为故障状态，转到步骤(9.3)仿真；
(9.3)判断步骤(9.2)中得到的输出事件是否为顶事件：若是，则不再进行故障树搜寻，转到步骤(9.4)；若否，则转到步骤(9.1)仿真；
(9.4)判断故障树模型其他子树中是否存在该故障事件：若是，则转到步骤(9.1)进行故障树搜寻；若否，转到步骤(10)。
9.根据权利要求8所述的考虑维修的含间歇性工作逻辑门的动态故障树仿真方法，其特征在于：所述步骤(9.1)包括以下分步骤：
(9.1.1)获得输入事件的状态及时间阶段：主事件状态和辅助事件的状态，间歇性工作逻辑门的输入事件均有三个状态：正常、故障和停机，分别用0、1和2表示；
(9.1.2)判断当前工作时间阶段：属于主单元工作时间阶段T1，跳转步骤(9.1.3)；属于辅助单元工作时间阶段T2，跳转步骤(9.1.6)；
(9.1.3)判断输入事件主单元状态：若输入事件主单元状态为0，跳转步骤(9.1.4)；若输入事件主单元状态为1，跳转步骤(9.1.5)；
(9.1.4)若辅助单元状态为2，输出0，跳转步骤(9.1.9)；否则，仿真输入错误，跳转步骤(9.1.9)
(9.1.5)若辅助单元状态为2，输出1，跳转步骤(9.1.9)；否则，仿真输入错误，跳转步骤(9.1.9)；
(9.1.6)若主单元状态2，跳转步骤(9.1.7)；若主单元状态1，跳转步骤(9.1.8)；否则，仿真输入错误，跳转步骤
(9.1.9)；
(9.1.7)判断辅助单元状态：若辅助单元状态为0，则输出2，跳转步骤(9.1.9)；若辅助单元状态为1，则输出1，跳转步骤(9.1.9)；否则，仿真输入错误，跳转步骤(9.1.9)；
(9.1.8)若辅助单元状态为0，则输出1，跳转步骤(9.1.9)；
否则，仿真输入错误，跳转步骤(9.1.9)；
(9.1.9)输出事件的状态，间歇性工作逻辑门仿真结束。
10.根据权利要求9所述的考虑维修的含间歇性工作逻辑门的动态故障树仿真方法，其特征在于：所述步骤(12)中，仿真数据包括：系统或者其他单元的工作时间，故障次数，维修时间，然后计算相应的可靠性指标，评价系统或者单元的可靠性。
CN201911268889.0A 2019-12-11 2019-12-11 考虑维修的含间歇性工作逻辑门的动态故障树仿真方法 Active CN111046556B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911268889.0A CN111046556B (zh) 2019-12-11 2019-12-11 考虑维修的含间歇性工作逻辑门的动态故障树仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911268889.0A CN111046556B (zh) 2019-12-11 2019-12-11 考虑维修的含间歇性工作逻辑门的动态故障树仿真方法

Publications (2)

Publication Number Publication Date
CN111046556A true CN111046556A (zh) 2020-04-21
CN111046556B CN111046556B (zh) 2021-10-26

Family

ID=70235809

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911268889.0A Active CN111046556B (zh) 2019-12-11 2019-12-11 考虑维修的含间歇性工作逻辑门的动态故障树仿真方法

Country Status (1)

Country Link
CN (1) CN111046556B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663226A (zh) * 2012-03-07 2012-09-12 北京航空航天大学 具有表决功能的可修混合公共贮备系统模型及可靠性仿真算法
CN106055729A (zh) * 2016-04-20 2016-10-26 西北工业大学 一种基于蒙特卡洛仿真的故障树分析方法
US20160327607A1 (en) * 2015-05-06 2016-11-10 Airbus Operations Limited Methods and systems for transforming fault tree diagrams of engineering systems
CN108388740A (zh) * 2018-03-01 2018-08-10 哈尔滨工程大学 一种基于故障树分析旁联冗余系统可靠性的方法
CN109740187A (zh) * 2018-12-11 2019-05-10 南京理工大学 基于fta-mc的起重机电气系统可靠性建模与仿真分析方法
CN110489932A (zh) * 2019-09-17 2019-11-22 北京航空航天大学 单事件条件约束门模型及仿真方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663226A (zh) * 2012-03-07 2012-09-12 北京航空航天大学 具有表决功能的可修混合公共贮备系统模型及可靠性仿真算法
US20160327607A1 (en) * 2015-05-06 2016-11-10 Airbus Operations Limited Methods and systems for transforming fault tree diagrams of engineering systems
CN106055729A (zh) * 2016-04-20 2016-10-26 西北工业大学 一种基于蒙特卡洛仿真的故障树分析方法
CN108388740A (zh) * 2018-03-01 2018-08-10 哈尔滨工程大学 一种基于故障树分析旁联冗余系统可靠性的方法
CN109740187A (zh) * 2018-12-11 2019-05-10 南京理工大学 基于fta-mc的起重机电气系统可靠性建模与仿真分析方法
CN110489932A (zh) * 2019-09-17 2019-11-22 北京航空航天大学 单事件条件约束门模型及仿真方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
LIU L , ZHAO G , SUN Y , ET AL.: "Reliability Analysis Method of Nuclear Cleaning System Considering Performance Degradation", 《019 PROGNOSTICS AND SYSTEM HEALTH MANAGEMENT CONFERENCE (PHM-QINGDAO)》 *
李志强; 徐廷学; 顾钧元; 刘玉东; 李凯: "基于动态贝叶斯网络的某控制单元可靠性分析", 《航空兵器》 *
海宁; 张彩珍; 孙国营; 于洋: "基于动态故障树的新型区间占检系统可用性分析", 《计算机应用与软件》 *
赵广燕; 孙宇锋; 赵刚: "电子产品性能与可靠性一体化设计分析技术研究", 《APPLIED COMPUTING,COMPUTER SCIENCE, AND COMPUTER ENGINEERING（ACC 2011 V4）》 *
钱虹; 骆建波; 金蔚霄; 周金明; 臧旭东: "基于动态故障树的DEH系统优化的可靠性评价", 《中国电力》 *

Also Published As

Publication number Publication date
CN111046556B (zh) 2021-10-26

Similar Documents

Publication Publication Date Title
Klein et al. Attribute-based architecture styles
CN108376221A (zh) 一种基于aadl模型扩展的软件系统安全性验证与评估方法
CN108388740B (zh) 一种基于故障树分析旁联冗余系统可靠性的方法
US20090083014A1 (en) Method for analyzing the reliability of technical installations with the use of physical models
CN106155035A (zh) 基于维修类数据的故障诊断方法及故障诊断系统
Hartonas-Garmhausen et al. ProbVerus: Probabilistic symbolic model checking
CN112560269A (zh) 基于Rhapsody状态机的高容错电子系统任务可靠性仿真分析方法
CN112287624A (zh) 基于生存分析加go法的继电保护系统可靠性评价方法
CN110175359B (zh) 基于业务流程的复杂系统安全性建模方法及装置
Lipaczewski et al. Using tool-supported model based safety analysis--Progress and experiences in SAML development
CN110489932A (zh) 单事件条件约束门模型及仿真方法
CN112765827B (zh) 一种功能相关系统的可靠性分析方法
CN108595959B (zh) 基于确定性随机Petri网的AADL模型安全性评估方法
Khan et al. Synergizing reliability modeling languages: BDMPs without repairs and DFTs
CN111046556B (zh) 考虑维修的含间歇性工作逻辑门的动态故障树仿真方法
Khan et al. Explaining boolean-logic driven markov processes using GSPNs
CN110531608B (zh) 基于冗余设计的高可靠电子设备定量fmeca分析方法和系统
Khan et al. Scalable reliability analysis by lazy verification
Bouissou et al. Boolean approximation for calculating the reliability of a very large repairable system with dependencies among components
Tang et al. Evaluation of software dependability based on stability test data
Fota et al. Dependability evaluation of an air traffic control computing system
CN110705051B (zh) 时间约束触改逻辑的动态故障树模型及仿真方法
CN111008310A (zh) 不考虑维修的间歇性工作逻辑门及其故障树的仿真方法
CN106354930B (zh) 一种空间飞行器的自适应重构方法及系统
Wang et al. Reliability modeling and verification method for dynamic reconfiguration system

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
GR01 Patent grant
GR01 Patent grant