CN111008310B - 不考虑维修的间歇性工作逻辑门及其故障树的仿真方法 - Google Patents

不考虑维修的间歇性工作逻辑门及其故障树的仿真方法 Download PDF

Info

Publication number
CN111008310B
CN111008310B CN201911268898.XA CN201911268898A CN111008310B CN 111008310 B CN111008310 B CN 111008310B CN 201911268898 A CN201911268898 A CN 201911268898A CN 111008310 B CN111008310 B CN 111008310B
Authority
CN
China
Prior art keywords
event
time
simulation
state
unit
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.)
Active
Application number
CN201911268898.XA
Other languages
English (en)
Other versions
CN111008310A (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 CN201911268898.XA priority Critical patent/CN111008310B/zh
Publication of CN111008310A publication Critical patent/CN111008310A/zh
Application granted granted Critical
Publication of CN111008310B publication Critical patent/CN111008310B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/90—Details of database functions independent of the retrieved data types
    • G06F16/901—Indexing; Data structures therefor; Storage structures
    • G06F16/9027—Trees

Landscapes

  • Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Test And Diagnosis Of Digital Computers (AREA)

Abstract

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

Description

不考虑维修的间歇性工作逻辑门及其故障树的仿真方法
技术领域
本发明涉及系统可靠性及故障仿真的技术领域，尤其涉及一种间歇性工作逻辑门的仿真方法，以及在不考虑维修情况下含有该间歇性工作逻辑门的系统故障树模型的仿真方法。
背景技术
随着科技水平的高速发展，系统的功能和组成越发庞大复杂，系统的工作环境也越发严酷。很多系统在运行过程中一旦发生故障将会导致巨大的经济损失甚至具有严重的安全隐患，这使得系统在工作过程需要非常可靠。
具体地，对于某核电系统来说，该系统具有较强的辐射性、组成单元众多，故障逻辑也较为复杂。该核电系统与其他系统相比，不仅在组成和规模上较为特殊，还存在系统呈间歇性工作的特点。
故障树分析方法是工程应用中适用性最强的可靠性建模分析方法。
FTA(Fault Tree Analysis，故障树)模型根据系统由上至下的结构功能对应关系，用与、或、非等逻辑门来表示单元与单元、单元与系统之间的故障因果逻辑关系。对于组成单元较少、单调关联系统的可靠性分析，都很容易建立其故障树模型，进而进行定性、定量分析，但对于单元数量较多，单元间耦合关系复杂的系统，FTA模型普遍存在问题是：无法准确描述系统的可靠性逻辑关系，即无法建立系统的FTA模型。
动态故障树模型是在故障树模型的基础上增加了几个动态逻辑门，扩展了故障树的适用范围，动态故障树模型不仅可以用于建立常见的串联、并联和表决等静态逻辑结构的故障树模型，还可以建立单元故障按照一定顺序发生等动态逻辑结构的故障树模型。常用的动态故障树的逻辑门符号包括顺序相关门、优先与门、功能触发门和备份门，备份门包括冷备门、温备门和热备门。
在现有的故障树逻辑门计算中，静态逻辑门如或门、与门以及表决门等可以通过结构函数，计算其输出事件的发生概率，而动态逻辑门如功能触发门、顺序相关门以及备份门需要转化为马尔科夫链计算，其计算方法也存在很大的局限性。
对于间歇性工作逻辑门来说，相对于现有的动态故障树逻辑门更加复杂，其中包含了事件的间歇性工作等特点。
针对这一问题，本文提出基于蒙特卡罗仿真的方法解决，建立间歇性工作逻辑门的仿真方法，来找出间歇性工作系统的薄弱环节，无论对于改进系统的设计，提高系统的可靠性水平，减轻安全问题发生，减少经济损失，都具有极为重要的意义。
发明内容
为克服现有技术的缺陷，本发明要解决的技术问题是提供了一种间歇性工作逻辑门仿真方法，还给出了不考虑维修情况下，含有间歇性工作逻辑门的故障树仿真方法，找出间歇性工作系统的薄弱环节，无论对于改进系统的设计，提高系统的可靠性水平，减轻安全问题发生，减少经济损失，都具有极为重要的意义。
本发明的技术方案是：这种间歇性工作逻辑门仿真方法，间歇性工作逻辑门包括：
输入事件：主单元(E1)和辅助单元(E2)的故障模式；
输出事件(Y)：是输入事件发生后，通过逻辑门的输出结果；
第一时间(T1)：表示主单元(E1)规定的工作时间；
第二时间(T2)：表示辅助单元(E2)规定的工作时间；
事件状态：包括正常、故障、停机三个状态，分别用0、1和2表示；其中，正常状态是指单元处于正常工作状态，即故障事件不发生且单元处于正常状态；故障状态是指单元发生了故障，不具备工作能力，即故障事件发生；停机状态是指单元工作一定时间之后，由于一些原因使得该单元暂时处于不可用状态，但是单元未发生故障，属于无故障停机状态，即故障事件不发生且单元处于无故障停机状态；该间歇性工作逻辑门模型表示：当主单元输入事件在工作期间发生故障之后或者达到规定的工作时间T1之后，输出事件为故障状态或者停机状态；当辅助单元输入事件为故障状态时，若现有输出事件为故障状态，则输出事件状态不变，若现有输出事件为停机状态，则输出事件变为故障状态；
该仿真方法包括以下步骤：
(1)获得输入事件的状态及时间阶段：主事件状态和辅助事件的状态；间歇性工作逻辑门的输入事件均有三个状态：正常、故障和停机，分别用0、1和2表示；
(2)判断当前工作时间阶段：若为主单元工作时间阶段(T1)跳转步骤(3)；若为辅助单元工作时间阶段(T2)跳转步骤(6)；
(3)判断输入事件主单元状态：若为输入事件主单元状态为0，跳转步骤(4)；若为输入事件主单元状态为1，跳转步骤(5)；
(4)若为辅助单元状态为2，输出0，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(5)若为辅助单元状态为2，输出1，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(6)判断主单元状态：若为主单元状态2，跳转步骤(7)；若为主单元状态1，跳转步骤(8)；否则，仿真输入错误，跳转步骤(9)；
(7)判断辅助单元状态：若为辅助单元状态为0，则输出2，跳转步骤(9)；若为辅助单元状态为1，则输出1，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(8)若为辅助单元状态为0，则输出1，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(9)输出事件的状态，间歇性工作逻辑门仿真结束。
本发明以基于蒙特卡罗仿真的方法解决，建立间歇性工作逻辑门的仿真方法，来找出间歇性工作系统的薄弱环节，无论对于改进系统的设计，提高系统的可靠性水平，减轻安全问题发生，减少经济损失，都具有极为重要的意义。
还提供了不考虑维修情况下含间歇性工作逻辑门的故障树仿真方法，该仿真方法包括以下步骤：
(I)仿真初始化：获得仿真输入信息，变量初始化，进行失效时间抽样，并对全部事件进行排序；
(II)获得最先发生事件及时间；
(III)判断最先发生事件类型：若是故障事件，则转到步骤(IV)；若是主单元切换事件，则转到步骤(V)；若是辅助单元切换事件，则转到步骤(VI)；
(IV)故障事件处理：更改因故障事件发生影响的底事件状态，更改仿真中因故障事件发生受影响的事件发生的时间，跳转步骤(VII)；
(V)主单元切换事件处理：更改因主单元切换事件影响的底事件状态，更改仿真中因主单元切换事件影响的事件发生的时间，跳转步骤(VII)；
(VI)辅助单元切换事件处理：更改因辅助单元切换事件受影响的底事件状态，更改仿真中因辅助单元切换事件影响的事件发生的时间，跳转步骤(VII)；
(VII)调用相关逻辑门仿真算法，通扫故障树，输出顶事件状态；
(VIII)判断仿真时间是否大于设定仿真时间，若否，则本次仿真还未结束，事件队列重新排序，转到步骤(II)继续仿真；若是，则本次仿真结束，继续步骤(IX)；
(IX)判断仿真次数是否大于设定仿真次数，若否，则仿真未结束，仿真次数加1，转到步骤(I)，重新进行失效时间抽样及排序，继续仿真；若是，则继续步骤(X)；
(X)统计仿真数据；
(XI)结束。
附图说明
图1是间歇性工作逻辑门的符号。
图2是根据本发明的间歇性工作逻辑门仿真方法的流程图。
图3是根据本发明的含间歇性工作逻辑门在不考虑维修情况下的故障树仿真方法的流程图。
图4是核电系统的功能结构图。
图5是核电系统A的故障树逻辑门。
图6是冷备份门仿真方法的流程图。
图7是核电系统故障时间的频率分布图。
图8是核电系统可靠度曲线。
具体实施方式
间歇性工作逻辑门包括：
输入事件：主单元(E1)和辅助单元(E2)的故障模式；
输出事件(Y)：是输入事件发生后，通过逻辑门的输出结果；
第一时间(T1)：表示主单元(E1)规定的工作时间；
第二时间(T2)：表示辅助单元(E2)规定的工作时间；
事件状态：包括正常、故障、停机三个状态，分别用0、1和2表示；其中，正常状态是指单元处于正常工作状态，即故障事件不发生且单元处于正常状态；故障状态是指单元发生了故障，不具备工作能力，即故障事件发生；停机状态是指单元工作一定时间之后，由于一些原因使得该单元暂时处于不可用状态，但是单元未发生故障，属于无故障停机状态，即故障事件不发生且单元处于无故障停机状态；该间歇性工作逻辑门模型表示：当主单元输入事件在工作期间发生故障之后或者达到规定的工作时间T1之后，输出事件为故障状态或者停机状态；当辅助单元输入事件为故障状态时，若现有输出事件为故障状态，则输出事件状态不变，若现有输出事件为停机状态，则输出事件变为故障状态；
如图2所示，这种间歇性工作逻辑门仿真方法包括以下步骤：
(1)获得输入事件的状态及时间阶段：主事件状态和辅助事件的状态；间歇性工作逻辑门的输入事件均有三个状态：正常、故障和停机，分别用0、1和2表示；
(2)判断当前工作时间阶段：若为主单元工作时间阶段(T1)跳转步骤(3)；若为辅助单元工作时间阶段(T2)跳转步骤(6)；
(3)判断输入事件主单元状态：若为输入事件主单元状态为0，跳转步骤(4)；若为输入事件主单元状态为1，跳转步骤(5)；
(4)若为辅助单元状态为2，输出0，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(5)若为辅助单元状态为2，输出1，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(6)判断主单元状态：若为主单元状态2，跳转步骤(7)；若为主单元状态1，跳转步骤(8)；否则，仿真输入错误，跳转步骤(9)；
(7)判断辅助单元状态：若为辅助单元状态为0，则输出2，跳转步骤(9)；若为辅助单元状态为1，则输出1，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(8)若为辅助单元状态为0，则输出1，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(9)输出事件的状态，间歇性工作逻辑门仿真结束。
本发明以基于蒙特卡罗仿真的方法解决，建立间歇性工作逻辑门的仿真方法，来找出间歇性工作系统的薄弱环节，无论对于改进系统的设计，提高系统的可靠性水平，减轻安全问题发生，减少经济损失，都具有极为重要的意义。
本发明的间歇性工作系统由主单元和辅助单元组成，两个单元均按照规定的时间间歇性工作，通过扩充输入事件的状态，分别为正常、故障和停机，利用蒙特卡罗抽样方法，结合基于事件的仿真思想，编写了间歇性工作逻辑门事件状态之间逻辑关系的仿真方法。该方法能够找出间歇性工作系统的薄弱环节，无论对于改进系统的设计，提高系统的可靠性水平，减轻安全问题发生，减少经济损失，都具有极为重要的意义。
有关含有间歇性工作逻辑特点的系统，以图4中某核电系统加以说明，该核电系统A包括A1子系统和A2子系统两部分。这两个子系统都是间歇性工作子系统，每个子系统中的主单元和辅助单元按照一定时间交替性工作，其中的单元均包含正常、故障和停机三种状态。
如图3所示，本发明还提供了不考虑维修情况下含间歇性工作逻辑门的故障树仿真方法，该仿真方法包括以下步骤：
(I)仿真初始化：获得仿真输入信息，变量初始化，进行失效时间抽样，并对全部事件进行排序；
(II)获得最先发生事件及时间；
(III)判断最先发生事件类型：若是故障事件，则转到步骤(IV)；若是主单元切换事件，则转到步骤(V)；若是辅助单元切换事件，则转到步骤(VI)；
(IV)故障事件处理：更改因故障事件发生影响的底事件状态，更改仿真中因故障事件发生受影响的事件发生的时间，跳转步骤(VII)；
(V)主单元切换事件处理：更改因主单元切换事件影响的底事件状态，更改仿真中因主单元切换事件影响的事件发生的时间，跳转步骤(VII)；
(VI)辅助单元切换事件处理：更改因辅助单元切换事件受影响的底事件状态，更改仿真中因辅助单元切换事件影响的事件发生的时间，跳转步骤(VII)；
(VII)调用相关逻辑门仿真算法，通扫故障树，输出顶事件状态；
(VIII)判断仿真时间是否大于设定仿真时间，若否，则本次仿真还未结束，事件队列重新排序，转到步骤(II)继续仿真；若是，则本次仿真结束，继续步骤(IX)；
(IX)判断仿真次数是否大于设定仿真次数，若否，则仿真未结束，仿真次数加1，转到步骤(I)，重新进行失效时间抽样及排序，继续仿真；若是，则继续步骤(X)；
(X)统计仿真数据；
(XI)结束。
优选地，所述步骤(I)包括以下分步骤：
(I.1)输入仿真所需信息：输入所有底事件的失效分布函数；设定仿真总时间SimAllTime、设定仿真总次数SimAllNumber；仿真时间用SimTime表示，仿真次数用SimNumber表示，事件队列用一个数组MyEvent()表示；设置主单元的规定运行时间为T1，辅助单元的规定运行时间为T2；设置所有底事件的初始状态，底事件状态表示：0代表正常，1代表故障，2代表停机等待；
(I.2)初始仿真次数SimNumber赋值为1：SimNumber＝1；
(I.3)系统仿真时间初始化：仿真时间SimTime赋值为0，
SimTime＝0；主单元切换时间MainSwitchTime赋值为0，
MainSwitchTime＝0；辅助单元切换时间AuxiliarySwitchTime赋值为0，AuxiliarySwitchTime＝0；单元失效时刻FailureTime赋值为0，FailureTime＝0；
(VI.4)失效时间抽样：利用蒙特卡罗抽样方法，抽取所有输入事件的失效时间；令主单元切换时间等于主单元的规定工作时间，即MainSwitchTime＝T1；辅助单元切换时间等于辅助单元的规定工作时间，即AuxiliarySwitchTime＝MainSwitchTime+T2；将抽取的失效时间和主、辅单元的切换时间放入事件队列MyEvent()中，并对事件进行排序。
优选地，所述步骤(II)中，从排序后的事件队列MyEvent中，获得最先发生事件的时间MyEvent(0)，仿真时间步进到该事件发生的时刻，即SimTime＝MyEvent(0)，将该事件从事件队列中移除。
优选地，所述步骤(IV)包括以下分步骤：
(IV.1)输入事件状态＝1；
(IV.2)判断故障发生单元是否为主单元：若是主单元故障，跳转步骤(IX)；若不是主单元故障，跳转步骤(X)；
(IV.3)辅助单元＝0，下次辅助单元切换时间＝仿真时间+规定工作时间T2，并重新放入事件队列进行更新，删除事件队列中主单元切换事件，转到步骤(X)。
优选地，所述步骤(V)包括以下分步骤：
(V.1)令切换事件的主单元＝2，切换事件的辅助单元＝0；
(V.2)切换事件的辅助单元下次切换时刻AuxiliarySwitchTime＝仿真时间SimTime+规定工作时间T2，并放入事件队列进行更新；
(V.3)切换事件的主单元下次失效的时刻＝下次该主单元失效事件发生时间+规定工作时间T2，并放入事件队列进行更新，转到步骤(X)。
优选地，所述步骤(VI)包括以下分步骤：
(VI.1)判断主单元是否发生故障：若主单元发生故障，切换事件的辅助单元＝2，切换事件的主单元＝1，将该事件的辅助单元下次失效的时刻从事件队列中删除，更新仿真队列，跳转步骤(X)；若主单元未发生故障，切换事件的辅助单元＝2，切换事件的主单元＝0，转到步骤(VI.2)；
(VI.2)切换事件的主单元下次切换时刻＝仿真时间+规定工作时间T1，并放入事件队列进行更新，转到步骤(VI.3)；
(VI.3)切换事件的辅助单元下次失效的时刻＝下次该辅助单元失效事件发生时间+规定工作时间T1，并放入事件队列进行更新，转到步骤(VII)。
优选地，所述步骤(VII)包括以下分步骤：
(VII.1)根据建立的故障树模型，将故障事件作为子节点，搜寻连接的故障树逻辑门，并调用与门、或门、间歇性工作逻辑门，获得输出事件的状态；调用间歇性工作逻辑门的算法如图2所示；
(VII.2)判断输出事件的状态是否为正常状态：若是，则转到步骤
(VII.4)；若否，则输出事件为故障状态，转到步骤(VII.3)仿真；
(VII.3)判断步骤(VII.2)中得到的输出事件是否为顶事件：若是，则不再进行故障树搜寻，转到步骤(VII.4)；若否，则转到步骤(VII.1)仿真；
(VII.4)判断故障树模型其他子树中是否存在该故障事件：若是，则转到步骤(VII.1)进行故障树搜寻；若否，转到步骤(VIII)。
对于核电系统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的功能结构图如图4所示：
根据间歇性工作逻辑门以及现有的动态故障树模型，建立核电系统A的故障树模型。故障树模型中，a表示核电系统A的故障模式，a1表示第一子系统A1的故障模式，a11表示第一主单元A11的故障模式，a12表示第一辅助单元A12的故障模式；a2表示第二子系统A2的故障模式，a21表示第二主单元A21的故障模式，a22表示第二辅助单元A22的故障模式。其故障树模型如图5所示。
在上述逻辑门中，存在两个动态逻辑门，冷备份门和间歇性工作逻辑门，其中，含有间歇性工作逻辑门的故障树模型的仿真方法如上所述，冷备份门的仿真方法需要重新建立。
图6冷备份门的仿真方法步骤如下：
步骤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
故障率(/h) 主单元运行时间(h) 辅助单元运行时间(h)
A11 0.005 200 \
A12 0.005 \ 150
A21 0.005 200 \
A22 0.005 \ 150
编写核电系统A的仿真程序，仿真输入为：每个单元的相关参数以及仿真次数3000次、仿真时间8760h。仿真计算核电系统的可靠性指标。
统计核电系统的使用可用度在总仿真次数3000次下，核电系统故障时间的频率分布图如图7、核电系统可靠度曲线如图8所示。
通过仿真程序计算核电系统的可靠度随系统运行时间的结果如表2所示。
表2
系统运行时间/h 50 100 150 200 250 300 350 400
可靠度仿真结果 0.97 0.91 0.82 0.74 0.59 0.46 0.36 0.29
以上所述，仅是本发明的较佳实施例，并非对本发明作任何形式上的限制，凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰，均仍属本发明技术方案的保护范围。

Claims (8)

1.不考虑维修情况下的故障树仿真方法，其特征在于：该方法包括间歇性工作逻辑门仿真方法，
间歇性工作逻辑门包括：
输入事件：主单元E1和辅助单元E2的故障模式；
输出事件Y：是输入事件发生后，通过逻辑门的输出结果；
第一时间T1：表示主单元E1规定的工作时间；
第二时间T2：表示辅助单元E2规定的工作时间；
事件状态：包括正常、故障、停机三个状态，分别用0、1和2表示；其中，正常状态是指单元处于正常工作状态，即故障事件不发生且单元处于正常状态；故障状态是指单元发生了故障，不具备工作能力，即故障事件发生；停机状态是指单元工作一定时间之后，由于一些原因使得该单元暂时处于不可用状态，但是单元未发生故障，属于无故障停机状态，即故障事件不发生且单元处于无故障停机状态；该间歇性工作逻辑门模型表示：当主单元输入事件在工作期间发生故障之后或者达到规定的工作时间T1之后，输出事件为故障状态或者停机状态；当辅助单元输入事件为故障状态时，若现有输出事件为故障状态，则输出事件状态不变，若现有输出事件为停机状态，则输出事件变为故障状态；
间歇性工作逻辑门仿真方法包括以下步骤：
(1)获得输入事件的状态及时间阶段：主事件状态和辅助事件的状态；间歇性工作逻辑门的输入事件均有三个状态：正常、故障和停机，分别用0、1和2表示；
(2)判断当前工作时间阶段：若为主单元工作时间阶段T1跳转步骤(3)；若为辅助单元工作时间阶段T2跳转步骤(6)；
(3)判断输入事件主单元状态：若为输入事件主单元状态为0，跳转步骤(4)；若为输入事件主单元状态为1，跳转步骤(5)；
(4)若为辅助单元状态为2，输出0，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(5)若为辅助单元状态为2，输出1，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(6)判断主单元状态：若为主单元状态2，跳转步骤(7)；若为主单元状态1，跳转步骤(8)；否则，仿真输入错误，跳转步骤(9)；
(7)判断辅助单元状态：若为辅助单元状态为0，则输出2，跳转步骤(9)；若为辅助单元状态为1，则输出1，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(8)若为辅助单元状态为0，则输出1，跳转步骤(9)；否则，仿真输入错误，跳转步骤(9)；
(9)输出事件的状态，间歇性工作逻辑门仿真结束。
2.根据权利要求1所述的不考虑维修情况下的故障树仿真方法，其特征在于：不考虑维修情况下的故障树仿真方法包括以下步骤：
(I)仿真初始化：获得仿真输入信息，变量初始化，进行失效时间抽样，并对全部事件进行排序；
(II)获得最先发生事件及时间；
(III)判断最先发生事件类型：若是故障事件，则转到步骤(IV)；若是主单元切换事件，则转到步骤(V)；若是辅助单元切换事件，则转到步骤(VI)；
(IV)故障事件处理：更改因故障事件发生影响的底事件状态，更改仿真中因故障事件发生受影响的事件发生的时间，跳转步骤(VII)；
(V)主单元切换事件处理：更改因主单元切换事件影响的底事件状态，更改仿真中因主单元切换事件影响的事件发生的时间，跳转步骤(VII)；
(VI)辅助单元切换事件处理：更改因辅助单元切换事件受影响的底事件状态，更改仿真中因辅助单元切换事件影响的事件发生的时间，跳转步骤(VII)；
(VII)调用相关逻辑门仿真算法，通扫故障树，输出顶事件状态；
(VIII)判断仿真时间是否大于设定仿真时间，若否，则本次仿真还未结束，事件队列重新排序，转到步骤(II)继续仿真；若是，则本次仿真结束，继续步骤(IX)；
(IX)判断仿真次数是否大于设定仿真次数，若否，则仿真未结束，仿真次数加1，转到步骤(I)，重新进行失效时间抽样及排序，继续仿真；若是，则继续步骤(X)；
(X)统计仿真数据；
(XI)结束。
3.根据权利要求2所述的不考虑维修情况下的故障树仿真方法，其特征在于：所述步骤(I)包括以下分步骤：
(I.1)输入仿真所需信息：输入所有底事件的失效分布函数；设定仿真总时间SimAllTime、设定仿真总次数SimAllNumber；仿真时间用SimTime表示，仿真次数用SimNumber表示，事件队列用一个数组MyEvent()表示；设置主单元的规定运行时间为T1，辅助单元的规定运行时间为T2；设置所有底事件的初始状态，底事件状态表示：0代表正常，1代表故障，2代表停机等待；
(I.2)初始仿真次数SimNumber赋值为1：SimNumber＝1；
(I.3)系统仿真时间初始化：仿真时间SimTime赋值为0，
SimTime＝0；主单元切换时间MainSwitchTime赋值为0，
MainSwitchTime＝0；辅助单元切换时间AuxiliarySwitchTime赋值为0，AuxiliarySwitchTime＝0；单元失效时刻FailureTime赋值为0，FailureTime＝0；
(I.4)失效时间抽样：利用蒙特卡罗抽样方法，抽取所有输入事件的失效时间；令主单元切换时间等于主单元的规定工作时间，即
MainSwitchTime＝T1；辅助单元切换时间等于辅助单元的规定工作时间，即AuxiliarySwitchTime＝MainSwitchTime+T2；将抽取的失效时间和主、辅单元的切换时间放入事件队列MyEvent()中，并对事件进行排序。
4.根据权利要求3所述的不考虑维修情况下的故障树仿真方法，其特征在于：所述步骤(II)中，从排序后的事件队列MyEvent中，获得最先发生事件的时间MyEvent(0)，仿真时间步进到该事件发生的时刻，即SimTime＝MyEvent(0)，将该事件从事件队列中移除。
5.根据权利要求4所述的不考虑维修情况下的故障树仿真方法，其特征在于：所述步骤(IV)包括以下分步骤：
(IV.1)输入事件状态＝1；
(IV.2)判断故障发生单元是否为主单元：若是主单元故障，跳转步骤(IX)；若不是主单元故障，跳转步骤(VII)；
(IV.3)辅助单元＝0，下次辅助单元切换时间＝仿真时间+规定工作时间T2，并重新放入事件队列进行更新，删除事件队列中主单元切换事件，转到步骤(VII)。
6.根据权利要求5所述的不考虑维修情况下的故障树仿真方法，其特征在于：所述步骤(V)包括以下分步骤：
(V.1)令切换事件的主单元＝2，切换事件的辅助单元＝0；
(V.2)切换事件的辅助单元下次切换时刻AuxiliarySwitchTime＝仿真时间SimTime+规定工作时间T2，并放入事件队列进行更新；
(V.3)切换事件的主单元下次失效的时刻＝下次该主单元失效事件发生时间+规定工作时间T2，并放入事件队列进行更新，转到步骤(VII)。
7.根据权利要求6所述的不考虑维修情况下的故障树仿真方法，其特征在于：所述步骤(VI)包括以下分步骤：
(VI.1)判断主单元是否发生故障：若主单元发生故障，切换事件的辅助单元＝2，切换事件的主单元＝1，将该事件的辅助单元下次失效的时刻从事件队列中删除，更新仿真队列，跳转步骤(VII)；若主单元未发生故障，切换事件的辅助单元＝2，切换事件的主单元＝0，转到步骤(VI.2)；
(VI.2)切换事件的主单元下次切换时刻＝仿真时间+规定工作时间T1，并放入事件队列进行更新，转到步骤(VI.3)；
(VI.3)切换事件的辅助单元下次失效的时刻＝下次该辅助单元失效事件发生时间+规定工作时间T1，并放入事件队列进行更新，转到步骤(VII)。
8.根据权利要求7所述的不考虑维修情况下的故障树仿真方法，其特征在于：所述步骤(VII)包括以下分步骤：
(VII.1)根据建立的故障树模型，将故障事件作为子节点，搜寻连接的故障树逻辑门，并调用与门、或门、间歇性工作逻辑门，获得输出事件的状态；
(VII.2)判断输出事件的状态是否为正常状态：若是，则转到步骤(VII.4)；若否，则输出事件为故障状态，转到步骤(VII.3)仿真；
(VII.3)判断步骤(VII.2)中得到的输出事件是否为顶事件：若是，则不再进行故障树搜寻，转到步骤(VII.4)；若否，则转到步骤(VII.1)仿真；
(VII.4)判断故障树模型其他子树中是否存在该故障事件：若是，则转到步骤(VII.1)进行故障树搜寻；若否，转到步骤(VIII)。
CN201911268898.XA 2019-12-11 2019-12-11 不考虑维修的间歇性工作逻辑门及其故障树的仿真方法 Active CN111008310B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911268898.XA CN111008310B (zh) 2019-12-11 2019-12-11 不考虑维修的间歇性工作逻辑门及其故障树的仿真方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911268898.XA CN111008310B (zh) 2019-12-11 2019-12-11 不考虑维修的间歇性工作逻辑门及其故障树的仿真方法

Publications (2)

Publication Number Publication Date
CN111008310A CN111008310A (zh) 2020-04-14
CN111008310B true CN111008310B (zh) 2023-08-25

Family

ID=70115425

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911268898.XA Active CN111008310B (zh) 2019-12-11 2019-12-11 不考虑维修的间歇性工作逻辑门及其故障树的仿真方法

Country Status (1)

Country Link
CN (1) CN111008310B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06102321A (ja) * 1992-09-22 1994-04-15 Nec Ic Microcomput Syst Ltd 論理回路の故障シミュレーション方法
EP1980964A1 (en) * 2007-04-13 2008-10-15 Yogitech Spa Method and computer program product for performing failure mode and effects analysis of an integrated circuit
EP2169486A2 (en) * 2008-09-30 2010-03-31 Honeywell International Inc. Intermittent fault detection and reasoning
FR3006470A1 (fr) * 2013-05-29 2014-12-05 Dassault Aviat Dispositif et procede informatises d'analyse de panne dans un systeme
CN104778370A (zh) * 2015-04-20 2015-07-15 北京交通大学 基于蒙特卡洛仿真求解动态故障树模型的风险分析方法
CN106021647A (zh) * 2016-05-06 2016-10-12 北京航空航天大学 一种基于割序集的动态故障树蒙特卡洛仿真定量计算方法
CN106055729A (zh) * 2016-04-20 2016-10-26 西北工业大学 一种基于蒙特卡洛仿真的故障树分析方法
CN107844641A (zh) * 2017-10-24 2018-03-27 北京航空航天大学 一种基于故障机理树和故障树相结合的可靠性建模方法
CN109270851A (zh) * 2018-08-17 2019-01-25 北京航空航天大学 人机交互动态故障树认知过载故障逻辑门的设计方法
CN110489932A (zh) * 2019-09-17 2019-11-22 北京航空航天大学 单事件条件约束门模型及仿真方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7681086B2 (en) * 2007-09-20 2010-03-16 Embraer- Empresa Brasileira De Aeronautica S.A. Fault tree map generation
US20160327607A1 (en) * 2015-05-06 2016-11-10 Airbus Operations Limited Methods and systems for transforming fault tree diagrams of engineering systems

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06102321A (ja) * 1992-09-22 1994-04-15 Nec Ic Microcomput Syst Ltd 論理回路の故障シミュレーション方法
EP1980964A1 (en) * 2007-04-13 2008-10-15 Yogitech Spa Method and computer program product for performing failure mode and effects analysis of an integrated circuit
EP2169486A2 (en) * 2008-09-30 2010-03-31 Honeywell International Inc. Intermittent fault detection and reasoning
FR3006470A1 (fr) * 2013-05-29 2014-12-05 Dassault Aviat Dispositif et procede informatises d'analyse de panne dans un systeme
CN104778370A (zh) * 2015-04-20 2015-07-15 北京交通大学 基于蒙特卡洛仿真求解动态故障树模型的风险分析方法
CN106055729A (zh) * 2016-04-20 2016-10-26 西北工业大学 一种基于蒙特卡洛仿真的故障树分析方法
CN106021647A (zh) * 2016-05-06 2016-10-12 北京航空航天大学 一种基于割序集的动态故障树蒙特卡洛仿真定量计算方法
CN107844641A (zh) * 2017-10-24 2018-03-27 北京航空航天大学 一种基于故障机理树和故障树相结合的可靠性建模方法
CN109270851A (zh) * 2018-08-17 2019-01-25 北京航空航天大学 人机交互动态故障树认知过载故障逻辑门的设计方法
CN110489932A (zh) * 2019-09-17 2019-11-22 北京航空航天大学 单事件条件约束门模型及仿真方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
计算机辅助可视化故障树分析系统;吴多胜等;《微计算机信息》(第8期);全文 *

Also Published As

Publication number Publication date
CN111008310A (zh) 2020-04-14

Similar Documents

Publication Publication Date Title
CN110221198B (zh) 基于故障树的煤矿井下组合开关故障诊断方法
EP3675306A1 (en) Quick search method for cascading failure of large ac/dc power grid, and system
US20090083014A1 (en) Method for analyzing the reliability of technical installations with the use of physical models
CN108388740B (zh) 一种基于故障树分析旁联冗余系统可靠性的方法
CN112487592B (zh) 基于贝叶斯网络的任务可靠性建模分析方法
CN103308824A (zh) 一种基于概率Petri网的电力系统故障诊断方法
Hartonas-Garmhausen et al. ProbVerus: Probabilistic symbolic model checking
CN105825045A (zh) 多阶段任务系统可修备件需求预测方法
Yevkin An improved modular approach for dynamic fault tree analysis
CN110474327B (zh) 配电网cps信息-物理组合预想故障生成方法及系统
CN111008310B (zh) 不考虑维修的间歇性工作逻辑门及其故障树的仿真方法
CN104597375A (zh) 一种电网故障诊断方法
CN117909864A (zh) 一种电力故障预测系统及方法
Abdelghany et al. Cause-consequence diagram reliability analysis using formal techniques with application to electrical power networks
CN116562167A (zh) 基于混合备份的发火测试系统可靠性评估及优化方法
CN109522640B (zh) 基于电流再分配的电路系统级联失效传播路径确定方法
Bouissou A simple yet efficient acceleration technique for Monte Carlo simulation
CN112016239A (zh) 一种基于动态贝叶斯网络快速求解系统可靠度的方法
CN110705051B (zh) 时间约束触改逻辑的动态故障树模型及仿真方法
CN111046556B (zh) 考虑维修的含间歇性工作逻辑门的动态故障树仿真方法
CN110175314A (zh) 一种基于ism法的软硬件混合指挥控制系统故障分层方法
Hsiao et al. Special protection system reliability assessment
Zhang et al. Key fault propagation path identification of CNC machine tools based on maximum occurrence probability
CN113673721A (zh) 一种基于深度强化学习的集群系统预防性维修方法
Niloofar et al. Learning Temporal Truth Tables of Dynamic Fault Trees from Time Series Data on Faults

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