CN103744295A - Reliability analysis system based on GO-FLOW methodology - Google Patents

Abstract

The invention relates to the technical field of industrial system reliability analysis, and particularly relates to a reliability analysis system based on a GO-FLOW methodology. The system comprises an OR gate operator, a AND gate operator, a NOT gate operator, a difference component operator, a signal generator operator, a delayer operator, a stage task operator, a two-stage component operator, a normally-closed valve operator, an normally-open valve operator, a constant probability fault device operator, a open-state valve fault operator, a close-state valve fault operator, and a valve motion fault operator. The reliability analysis system based on the GO-FLOW methodology can be used for security assessment of an industrial system and calculation of fault probability; factors, such as common signals and the like, are taken into account in the model; and thus application range of the GO-FLOW methodology is widened, and dynamic fault probability of a large-scale complex system can be efficiently calculated. The reliability analysis system based on the GO-FLOW methodology is fine and clear in structure; model parameters are clear in meaning and easy to obtain; and modeling process is fast, highly-efficient, precise and strong in practicability.

Description

A kind of reliability analysis system based on GO-FLOW method

Technical field

The present invention relates to industrial system reliability analysis technology field, be specifically related to a kind of reliability analysis system based on GO-FLOW method.

Background technology

Along with industrial system development, the loss that system scale is day by day huge, structure is increasingly sophisticated, the system failure or inefficacy bring is also huger.Therefore, the analytical approach of system reliability also more and more receives people's concern.This technology has become estimated risk, managing risk, has fallen low-risk important tool.In systems reliability analysis method, fault tree analysis (FTA, Failure Tree Analysis) is a kind of the most frequently used effective ways, and is applied in reliability, safety analysis and the risk assessment of various systems.Yet fault tree analysis (FTA, Failure Tree Analysis) has certain limitation, and the modeling of complication system is had to suitable difficulty, particularly for having multiple case and having signal feedback and the system of timing variations.

GO-FLOW sources of law are in GO method (GO methodology), it is a kind of successfully system reliability Safety Analysis Method as leading of take, it is applicable to the safety analysis of the Modernization System of multimode, multiple timings, to there being the safety analysis of production run of actual logistics especially applicable, the method has been applied to the systems such as oil, chemical industry, nuclear industry.GO-FLOW method mainly adopts the mode of figure deduction, describes the function of parts/subsystem by corresponding GO-FLOW operational symbol.Thereby can directly the schematic diagram of system be converted to illustraton of model, then the probability occurring according to the various states of illustraton of model computing system.The method is mainly used in having complex time sequence or state time varying system.

The application that GO-FLOW method relates to is not extensive, and its main cause is that traditional GO-FLOW model analysis method computation process is complicated, workload is huge, is difficult to meet efficiently, requirement accurately.Therefore, a kind of reliability analysis system based on GO-FLOW method of the present invention, the feature such as its extendability is strong, counting yield is high, simulation accuracy is high, modeling is convenient, has important application value.

Summary of the invention

Above-mentioned technical matters of the present invention is mainly solved by following technical proposals:

A reliability analysis system based on GO-FLOW method, is characterized in that, comprises

One operational symbol analysis module: for describing the logical relation of each signal (flow, voltage, electric current etc.) of actual industrial system, comprise AOI logic, differential logic, time delay logic and phased mission logic;

One parts analysis module: comprise two state models and unit failure model for describing the real system parts corresponding with real system, the control signal of parts is produced by signal generator model;

One Switch Analysis module: for describing switch, the valve element of actual industrial system, comprise normal closed gate model, normally open valve door model, close failsafe valve model, open failsafe valve model and valve opening and closing action model.

At above-mentioned a kind of reliability analysis system based on GO-FLOW method, described operational symbol analysis module comprises or door analytic unit, with door an analytic unit, not gate analytic unit, differentiation element, delay unit and phased mission unit; Described parts analysis module comprises two state analysis unit, signal generator unit and unit failure unit; Described Switch Analysis module comprises normal closed gate analytic unit, normally open valve door analytic unit, closes failsafe valve analytic unit, opens failsafe valve analytic unit and valve switch motion analysis unit.

At above-mentioned a kind of reliability analysis system based on GO-FLOW method, described or door analytic unit comprise common or door analyze subelement and with signal correction or door analyze subelement; Described and door analytic unit comprise common and door and analyze subelement and analyze subelement with signal correction with door.

At above-mentioned a kind of reliability analysis system based on GO-FLOW method, the course of work of described or door analytic unit is based on following formula:

$R\left(t\right)=1-\underset{t=1}{\overset{n}{\mathrm{\Π}}}(1-S_i\left(t\right))$

Wherein: R (t) is output signal S _i(t) be input signal;

Described a kind of reliability analysis system based on GO-FLOW method, is characterized in that, described with the course of work of door analytic unit based on following formula:

$R\left(t\right)=\underset{t=1}{\overset{n}{\mathrm{\Π}}}{S}_i\left(t\right)$

Wherein: R (t) is output signal S _i(t) be input signal;

Described with signal correction or door analyze subelement comprise two inputs or door and three inputs or analysis, based on following formula::

The output of two inputs or door is characterized by:

$R\left(t\right)=S_1\left(t\right)+S_2\left(t\right)-\frac{S_1\left(t\right)S_2\left(t\right)}{C\left(t\right)}$

Wherein: R (t) is output signal S ₁(t), S ₂(t) be input signal, C (t) is shared signal;

The output of three inputs or door is characterized by:

$R\left(t\right)=S_1\left(t\right)+S_2\left(t\right)+S_3\left(t\right)-\frac{S_1\left(t\right)S_2\left(t\right)+S_2\left(t\right)S_3\left(t\right)+S_1\left(t\right)S_3\left(t\right)}{C\left(t\right)}+\frac{S_1\left(t\right)S_2\left(t\right)S_3\left(t\right)}{C^2\left(t\right)}$

Wherein: R (t) is output signal S ₁(t), S ₂(t), S ₃(t) be input signal, C (t) is shared signal;

Described with signal correction analyze with door the analysis that subelement comprises two inputs and door and three value and gate, based on following formula::

Two inputs are characterized by with the output of door:

$R\left(t\right)=\frac{S_1\left(t\right)S_2\left(t\right)}{C\left(t\right)}$

Wherein: R (t) is output signal S ₁(t), S ₂(t) be input signal, C (t) is shared signal;

The output of three value and gate is characterized by:

$R\left(t\right)=\frac{S_1\left(t\right)S_2\left(t\right)S_3\left(t\right)}{C^2\left(t\right)}$

Wherein: R (t) is output signal S ₁(t), S ₂(t), S ₃(t) be input signal, C (t) is shared signal.At above-mentioned a kind of reliability analysis system based on GO-FLOW method, the course of work of described not gate analytic unit is based on following formula:

R(t)=1-S(t)

Wherein: R (t) is output signal S _i(t) be input signal.

At above-mentioned a kind of reliability analysis system based on GO-FLOW method, the course of work of described difference unit is based on following formula:

$\left\{\begin{array}{c}R\left(t\right)=S\left(t\right)-O\left(t^{\′}\right)\\ O\left(t\right)=S\left(t\right)\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t') was the signal of O (t) in a upper moment;

The course of work of described delay unit is based on following formula:

$\left\{\begin{array}{c}R\left(t\right)=O\left(t^{\′\left(n\right)}\right)\\ O\left(t\right)=S\left(t\right)\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t' ⁽ⁿ⁾) be O (t) at front n signal constantly, n is delay adjustments value;

The course of work of described phased mission unit is based on following formula:

$\left\{\begin{array}{c}R\left(t\right)=S\left(t\right)-O\left(t^{\′}\right)\\ O\left(t\right)=S\left(t\right)\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t') was the signal of O (t) in a upper moment.

At above-mentioned a kind of reliability analysis system based on GO-FLOW method, the course of work of described two state analysis unit is based on following formula:

R(t)=P _g·S(t)

Wherein: R (t) is that output signal S (t) is input signal, P _gfor element failure rate;

The course of work of described signal generator unit is based on following formula:

Wherein: R (t) is output signal, t ₁, t ₂the time that starts and finish for prearranged signal;

The course of work of described unit failure unit is based on following formula:

$R\left(t\right)=S\left(t\right)\·\exp \{-\mathrm{\λ}\·\underset{i}{\mathrm{\Σ}}\underset{t_k\≤t}{\mathrm{\Σ}}{P}_i\left(t_k\right)\×\min [1,S\left(t_k\right)/S\left(t\right)\left]\right\}$

Wherein: R (t) is that output signal S (t) is input signal, and λ is element failure rate, P _i(t _k) be the value of i control signal when time t (k), S (t _k) be the value of input signal S (t) when time t (k).

At above-mentioned a kind of reliability analysis system based on GO-FLOW method,

The course of work of described normal closed gate analytic unit is based on following formula:

$\left\{\begin{array}{c}R\left(t\right)=S\left(t\right)\·O\left(t\right)\\ O\left(t_1\right)=P_p\\ O\left(t\right)=O\left(t^{\′}\right)+[1-O\left(t^{\′}\right)]\·P\left(t\right)\·P_g\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t') was the signal of O (t) in a upper moment; t ₁for initial time, P _pfor the probability that signal is not opened to valve, P by mistake _gfor the signal probability that then valve is successfully opened;

The course of work of described normally open valve door analytic unit is based on following formula:

$\left\{\begin{array}{c}R\left(t\right)=S\left(t\right)\·O\left(t\right)\\ O\left(t_1\right)={1-P}_p\\ O\left(t\right)=O\left(t^{\′}\right)[1-P\left(t\right)\·P_g]\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t') was the signal of O (t) in a upper moment; t ₁for initial time, P _pfor the probability that signal is not opened to valve, P by mistake _gfor the signal probability that then valve is successfully opened;

The described course of work of closing failsafe valve analytic unit is based on following formula:

$R\left(t\right)=S\left(t\right)\·\{1-\exp [-\mathrm{\λ}\·\underset{i}{\mathrm{\Σ}}\underset{t_k\≤t}{\mathrm{\Σ}}{P}_i\left(t_k\right)\left]\right\}$

Wherein: R (t) is that output signal S (t) is input signal, and λ is element failure rate, P _i(t _k) be the value of i control signal when time t (k);

The described course of work of opening failsafe valve analytic unit is based on following formula:

$R\left(t\right)=S\left(t\right)\·\exp [-\mathrm{\λ}\·\underset{i}{\mathrm{\Σ}}\underset{t_k\≤t}{\mathrm{\Σ}}{P}_i\left(t_k\right)\}$

Wherein: R (t) is that output signal S (t) is input signal, and λ is element failure rate, P _i(t _k) be the value of i control signal when time t (k)

The course of work of described valve switch motion analysis unit is based on following formula:

$R\left(t\right)=\left\{\begin{array}{c}S\left(t\right)\&CenterDot;\{O\left(t^{\&prime;}\right)+[1-O\left(t^{\&prime;}\right)]\&CenterDot;P_1\left(t\right)\&CenterDot;P_g\},{\mathrm{<figure-callout\; id="1"\; label="P"\; filenames="HDA0000433323520000012.png,HDA0000433323520000013.png"\; state="\{\{state\}\}">P</figure-callout>}}_1=1\\ S\left(t\right)\&CenterDot;O\left(t^{\&prime;}\right)[1-P_2\left(t\right)\&CenterDot;P_g],{\mathrm{<figure-callout\; id="2"\; label="P"\; filenames="HDA0000433323520000012.png,HDA0000433323520000013.png"\; state="\{\{state\}\}">P</figure-callout>}}_2=1\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t') is O (t) at upper one signal constantly, P ₁for start signal, P ₂for shutdown signal, P _gfor opening or the shutdown signal probability that then valve is successfully opened or closed.

Therefore, tool of the present invention has the following advantages: 1. meet GO-FLOW basic modeling rule, model structure is meticulous, clear, and each ingredient intrinsic parameter meaning of model clearly, easily obtain; 2. the simulation algorithm carrying in conjunction with simulink, can automatically carry out analytical calculation to GO-FLOW modular concept figure, and calculating process is quick, efficient, accurate; 3. considered in GO-FLOW system modelling process the correction problem about shared signal amount, closer to reality, can accurately reflect the situation of real system; 4. be applicable to the modeling analysis of large-scale complicated system, and modeling required time is short, expense is low, practical.

Accompanying drawing explanation

Fig. 1 structural principle schematic diagram of the present invention.

Fig. 2 a is for being two inputs or gating element model schematic diagram.

Fig. 2 b is two input or the gating element model schematic diagram with shared signal correction.

Fig. 2 c is three input or the gating element model schematic diagram with shared signal correction.

Fig. 3 a is two inputs and gating element model schematic diagram.

Fig. 3 b is two inputs and gating element model schematic diagram with shared signal correction.

Fig. 3 c is the three value and gate component models schematic diagram with shared signal correction.

Fig. 4 is not gate component models schematic diagram.

Fig. 5 a is difference engine component models schematic diagram.

Fig. 5 b is chronotron component models schematic diagram.

Fig. 5 c phased mission event model schematic diagram.

Fig. 6 a is two state element model schematic diagram.

Fig. 6 b is signal generator component model schematic diagram.

Fig. 6 c is unit failure component models schematic diagram.

Fig. 7 a is normal closed gate component models.

Fig. 7 b is normally open valve gating element model.

Fig. 8 a is the fault model that closes state valve.

Fig. 8 b is out the fault model of state valve.

Fig. 8 c is for opening and closing the fault model of valve.

Fig. 9 a is for implementing the circuit diagram in example 1.

Fig. 9 b is for implementing the GO-FLOW system diagram in example 1.

Fig. 9 c is for implementing the result of calculation of example 1.

Figure 10 is example 2GO-FLOW system diagram.

Embodiment

Below by embodiment, and by reference to the accompanying drawings, technical scheme of the present invention is described in further detail.

Embodiment:

Native system embeds in matlab/simulink so that card format to be installed, and the signaling interface with other element of simulink is provided, and can build flexibly the probability model of complication system.System hierarchy as shown in Figure 1.

Wherein, each operational symbol inner structure is described as follows:

1, operational symbol model

(1) operational symbol 22: or door

Or door have a plurality of input and output, when having at least an input signal to exist, output S _iexist.That is:

$R\left(t\right)=1-\underset{t=1}{\overset{n}{\mathrm{\&Pi;}}}(1-S_i\left(t\right))$

Be input as example with two, its simulink model is as Fig. 2 a.In model, utilize tri-kinds of modules of Add, Product and Constant to realize transport function expression formula.For the parts of a plurality of inputs, can be by a plurality of two input block series connection equivalences.

In the GO-FLOW figure of engineering technology system, often there is shared signal, so also inevitably will be in the face of comprising the logic gate computational problem of shared signal input.Shared signal mainly produces in the system that has stand-by equipment.Due to each spare unit phase simple crosscorrelation, while making to calculate, can not calculate by independent probability, must consider modifying factor.Therefore, set up and consider under shared signal or the computation model of door.With shared signal two input or door model and three input or door model respectively as shown in Fig. 2 b, Fig. 2 c.

Operational symbol 30: with door

With the logical and relation of a plurality of signals of door operator representation, its output valve is the probability that all inputs all exist,

$R\left(t\right)=\underset{t=1}{\overset{n}{\mathrm{\&Pi;}}}{S}_i\left(t\right)$

Its simulink model as shown in Figure 3 a.

Consider under shared signal the computation model with door, and or class seemingly.When there is shared signal, need to adopt the shared signal model with signal correction.With shared signal two input with door and three value and gate respectively as Fig. 3 b, bc

(3) operational symbol 23 not gates

Not gate operational symbol has been mainly used in the relation of logic NOT in probability calculation, and its arithmetic unit is as Fig. 4

(4) operational symbol 24: difference engine

Difference engine is the difference relational calculus for analog input amount, and Simulink model is as Fig. 5 a.

(5) operational symbol 28 chronotrons

Chronotron operational symbol is mainly applicable to exist the model construction of the dynamic system of time delay process, and for completing the computing to signal lag operation, its system architecture as shown in Figure 5 b.

(6) operational symbol 40: phased mission

The output input relationship expression of phased mission is as follows:

$R\left(t\right)=\left\{\begin{array}{c}1,t{<t}_i\\ S\left(t\right),t_i\&le;t\&le;t_j\\ S\left(t_j\right),t{>t}_j\end{array}\right.$

Arithmetic logic is as Fig. 5 c.Wherein, the function of simulink custom block is by obtaining the system time judgement current task stage to determine output valve.

2, partial model

(1) 21: two state elements of operational symbol

This operational symbol is usually used in representing only to have the element of normal work and fault two states, as resistance, and circuit, diode etc.Assumed fault rate is P _g, input signal and output signal can be expressed as

R(t)=P _g·S(t)

From mathematic(al) representation, output probability equals input probability and is multiplied by failure rate, can build by gain module in simulink.Corresponding operational symbol model is as Fig. 6 a.

(2) operational symbol 25: signal generator

Signal generator is mainly used to produce the control signal of modules, discrete also continuous, the multiplex step signal that produces of signal generator in the practical application of GO-FLOW, herein for this operational symbol is provided with two parameters, step start time and step duration, and the time reference that provides of simulation clock produces respective waveforms, its simulink model is as Fig. 6 b.Matlab Function custom block wherein for the comparison system time whether in user's setup times interval, thereby produce corresponding output signal

(3) operational symbol 35: unit failure

Unit failure is the most frequently used operational symbol in GO-FLOW fundamental operation symbol, and it is mainly used in describing the parts that lost efficacy when in running order.Its model has a primary input signal, a plurality of branches input signal and an output signal.Suppose that part failure rate is constant λ, output valve is:

$R\left(t\right)=S\left(t\right)\&CenterDot;\exp \{-\mathrm{\&lambda;}\&CenterDot;\underset{i}{\mathrm{\&Sigma;}}\underset{t_k\&le;t}{\mathrm{\&Sigma;}}{P}_i\left(t_k\right)\&times;\min [1,S\left(t_k\right)/S\left(t\right)\left]\right\}$

As can be seen from the above equation, calculating each R (t) constantly needs to consider from initial t ₁each P (t constantly constantly _k) and S (t _k) functional value.Also just mean and need the passing all P (t of storage of array _k) and S (t _k) value.By common simulink module, be difficult to realize, we utilize the self-defining function function in simulink, by M file, realize self-defined emulation module, complete building of summation module in model.Add the modules such as constant, product, MathFunction, structure simulink subsystem is as Fig. 6 c.In Matlab Function module, by recording the value of each moment s (t) and p (t), be used for calculating next output constantly.

3. switch model

(1) operational symbol 26: normal closed gate

Normal closed gate, in the normal state in closed condition, only just can be opened when branch's input signal exists.This operational symbol consists of a primary input, a control inputs and an output.The probability that valve is successfully opened when opening signal arrives is P _g, opening signal not to and the probability of valve in open mode is P _p, the output valve of output signal is:

$\left\{\begin{array}{c}R\left(t\right)=S\left(t\right)\&CenterDot;O\left(t\right)\\ O\left(t_1\right)=P_p\\ O\left(t\right)=O\left(t^{\&prime;}\right)+[1-O\left(t^{\&prime;}\right)]\&CenterDot;P\left(t\right)\&CenterDot;P_g\end{array}\right.$

Wherein: O (t') represents that valve engraves the probability of opening while carving for the moment, t when t ₁for initial time.When building simulink arithmetical unit, utilize O (t) value in a moment in memory component stores.Build transport function expression formula utilizing other common mathematical computing module.Constructed simulink model is as Fig. 7 a.In addition, this operational symbol also can represent normally opened contact

(2) operational symbol 27: normally open valve door

Similar with normal closed gate, when normally working, this valve in open mode, under control signal effect, closes.Therefore it also has a primary input and a control inputs, an output.When control signal arrives, valve is with P _gprobability successfully close, shutdown signal not to and the probability of valve in closed condition is P _p, module output valve is:

$\left\{\begin{array}{c}R\left(t\right)=S\left(t\right)\&CenterDot;O\left(t\right)\\ O\left(t_1\right)=1-P_p\\ O\left(t\right)=O\left(t^{\&prime;}\right)[1-P\left(t\right)\&CenterDot;P_g]\end{array}\right.$

Can construct thus the simulink parts of normal closed gate as Fig. 3 .c.In addition, this operational symbol also can represent normally closed contact.

(3) operational symbol 37: the fault of opening state valve

The valve that this operational symbol model lost efficacy while mainly representing out state.Model hypothesis part failure rate is constant λ, and this operational symbol is by a primary input, and the input of a plurality of branches and an output form.The mathematic(al) representation of output valve is:

$R\left(t\right)=S\left(t\right)\&CenterDot;\exp [-\mathrm{\&lambda;}\&CenterDot;\underset{i}{\mathrm{\&Sigma;}}\underset{t_k\&le;t}{\mathrm{\&Sigma;}}{P}_i\left(t_k\right)\}$

In the building process of simulink model, utilize Memory and Add assembly to realize cumulative summation function.Open the fault model of state valve as Fig. 8 a.

(4) operational symbol 38: the fault of closing state valve

With open status valve class seemingly, crash rate is that the fault operational symbol of the state that the closes valve of constant λ has a primary input, the input of a plurality of branches and an output signal, its mathematic(al) representation is:

$R\left(t\right)=S\left(t\right)\&CenterDot;\{1-\exp [-\mathrm{\&lambda;}\&CenterDot;\underset{i}{\mathrm{\&Sigma;}}\underset{t_k\&le;t}{\mathrm{\&Sigma;}}{P}_i\left(t_k\right)\left]\right\}$

Can build thus simulink model as Fig. 8 b.

(5) operational symbol 39 can open and close valve

The comprehensive operational symbol 26 of this operational symbol and operational symbol 27, be that the valve of frequent switching is carried out to modeling, and this model has a primary input, the input of Liang Ge branch and an output signal, and the probability that valve is successfully opened is P _o, the probability of successfully closing is P _c.

As opening signal P ₁during arrival, output valve is

R(t)=S(t)·{O(t′)+[1-O(t′)]·P ₁(t)·P _g}

As shutdown signal P ₂during arrival, output valve is

R(t)=S(t)·O(t′)[1-P ₂(t)·P _g]

In expression formula, O (t') for valve t constantly before in opening shape probability of state.

Structure simulink operator system as shown in Figure 8 c.By Constant, Add, Product first to P ₁, P ₂process, prevent that the error-logic that two signals arrive simultaneously from occurring.Utilize Memory assembly to realize the storage of signal O (t) is used for producing O (t') signal.In logical diagram, Switch module realizes P ₁, P ₂the selectivity output of two signal expressions.Work as P ₁, P ₂signal does not all arrive, i.e. P ₁=P ₂=0 o'clock, two expression equivalences, were all degenerated to R (t)=S (t) O (t ') so select for a post Yi Yi road output and all can.

GO-FLOW method is mainly applied the reliability safety analysis with industrial circle, the application in conjunction with two example demonstration this patent instruments in real system fail-safe analysis

Implement example 1: circuit dynamic fault analysis modeling

Fig. 9 a is a ball bearing made using system, S during t=0 ₁closure, S during t=10 ₂closed.Wherein the mean lifetime of battery is 1000 hours, analyzes the probability that has the normal work of a lamp at least

Wherein battery failures probability is that 0.1 line fault probability is 0.2, and switch normally cuts out and just cuts out probability 0.1 bulb normal working probability under probability 0.7 original state is by mistake 0.8.

System GO-FLOW arithmograph is as Fig. 9 b:

Adopt fixed step size emulation, step-length is 1 hour, and simulation time 200 hours obtains system reliability temporal evolution result as Fig. 9 c:

Implement the modeling of example 2:GO-FLOW Power System Analysis

Every unit normal condition of certain nuclear power station is powered to power plant by two working bus bars (LGB, LGB) after class, after main electrical breakdown, by two 6.6kV emergency bus (LHA, LHB), to power plant, is powered.When power plant loses after whole external powers, can continue power supply by two diesel-driven generators (LHP, LHQ).After diesel generator set continues to lose efficacy, enter the whole audience loss of power accident, can come recovery system to power by recovering primary power or emergency power pack or diesel generator set.Its electric power system reliability structure is as Figure 10, failure probability that can direct solution whole system after operation simulink.And can directly read result, and this nuclear power station electric power system reliability under normal operation mode is 0.999440, failure probability is 5.599557e-04.

Specific embodiment described herein is only to the explanation for example of the present invention's spirit.Those skilled in the art can make various modifications or supplement or adopt similar mode to substitute described specific embodiment, but can't depart from spirit of the present invention or surmount the defined scope of appended claims.

Claims (8)

1. the reliability analysis system based on GO-FLOW method, is characterized in that, comprises

One operational symbol analysis module: for describing the logical relation of each signal (flow, voltage, electric current etc.) of actual industrial system, comprise AOI logic, differential logic, time delay logic and phased mission logic;

One parts analysis module: comprise two state models and unit failure model for describing the real system parts corresponding with real system, the control signal of parts is produced by signal generator model;

One Switch Analysis module: for describing switch, the valve element of actual industrial system, comprise normal closed gate model, normally open valve door model, close failsafe valve model, open failsafe valve model and valve opening and closing action model.

2. a kind of reliability analysis system based on GO-FLOW method according to claim 1, it is characterized in that, described operational symbol analysis module comprises or door analytic unit, with door an analytic unit, not gate analytic unit, differentiation element, delay unit and phased mission unit; Described parts analysis module comprises two state analysis unit, signal generator unit and unit failure unit; Described Switch Analysis module comprises normal closed gate analytic unit, normally open valve door analytic unit, closes failsafe valve analytic unit, opens failsafe valve analytic unit and valve switch motion analysis unit.

3. a kind of reliability analysis system based on GO-FLOW method according to claim 1, is characterized in that, described or door analytic unit comprise common or door analyze subelement and with signal correction or door analyze subelement; Described and door analytic unit comprise common and door and analyze subelement and analyze subelement with signal correction with door.

4. a kind of reliability analysis system based on GO-FLOW method according to claim 1, is characterized in that, the course of work of described or door analytic unit is based on following formula:

$R\left(t\right)=1-\underset{t=1}{\overset{n}{\mathrm{\&Pi;}}}(1-S_i\left(t\right))$

Wherein: R (t) is output signal S _i(t) be input signal;

Described a kind of reliability analysis system based on GO-FLOW method, is characterized in that, described with the course of work of door analytic unit based on following formula:

$R\left(t\right)=\underset{t=1}{\overset{n}{\mathrm{\&Pi;}}}{S}_i\left(t\right)$

Wherein: R (t) is output signal S _i(t) be input signal;

Described with signal correction or door analyze subelement comprise two inputs or door and three inputs or analysis, based on following formula::

The output of two inputs or door is characterized by:

$R\left(t\right)=S_1\left(t\right)+S_2\left(t\right)-\frac{S_1\left(t\right)S_2\left(t\right)}{C\left(t\right)}$

Wherein: R (t) is output signal S ₁(t), S ₂(t) be input signal, C (t) is shared signal;

The output of three inputs or door is characterized by:

$R\left(t\right)=S_1\left(t\right)+S_2\left(t\right)+S_3\left(t\right)-\frac{S_1\left(t\right)S_2\left(t\right)+S_2\left(t\right)S_3\left(t\right)+S_1\left(t\right)S_3\left(t\right)}{C\left(t\right)}+\frac{S_1\left(t\right)S_2\left(t\right)S_3\left(t\right)}{C^2\left(t\right)}$

Wherein: R (t) is output signal S ₁(t), S ₂(t), S ₃(t) be input signal, C (t) is shared signal;

Described with signal correction analyze with door the analysis that subelement comprises two inputs and door and three value and gate, based on following formula::

Two inputs are characterized by with the output of door:

$R\left(t\right)=\frac{S_1\left(t\right)S_2\left(t\right)}{C\left(t\right)}$

Wherein: R (t) is output signal S ₁(t), S ₂(t) be input signal, C (t) is shared signal;

The output of three value and gate is characterized by:

$R\left(t\right)=\frac{S_1\left(t\right)S_2\left(t\right)S_3\left(t\right)}{C^2\left(t\right)}$

Wherein: R (t) is output signal S ₁(t), S ₂(t), S ₃(t) be input signal, C (t) is shared signal.

5. a kind of reliability analysis system based on GO-FLOW method according to claim 1, is characterized in that, the course of work of described not gate analytic unit is based on following formula:

R(t)=1-S(t)

Wherein: R (t) is output signal S _i(t) be input signal.

6. a kind of reliability analysis system based on GO-FLOW method according to claim 1, is characterized in that, the course of work of described difference unit is based on following formula:

$\left\{\begin{array}{c}R\left(t\right)=S\left(t\right)-O\left(t^{\&prime;}\right)\\ O\left(t\right)=S\left(t\right)\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t') was the signal of O (t) in a upper moment;

The course of work of described delay unit is based on following formula:

$\left\{\begin{array}{c}R\left(t\right)=O\left(t^{\&prime;\left(n\right)}\right)\\ O\left(t\right)=S\left(t\right)\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t' ⁽ⁿ⁾) be O (t) at front n signal constantly, n is delay adjustments value;

The course of work of described phased mission unit is based on following formula:

$\left\{\begin{array}{c}R\left(t\right)=S\left(t\right)-O\left(t^{\&prime;}\right)\\ O\left(t\right)=S\left(t\right)\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t') was the signal of O (t) in a upper moment.

7. a kind of reliability analysis system based on GO-FLOW method according to claim 1, is characterized in that, the course of work of described two state analysis unit is based on following formula:

R(t)=P _g·S(t)

Wherein: R (t) is that output signal S (t) is input signal, P _gfor element failure rate;

The course of work of described signal generator unit is based on following formula:

Wherein: R (t) is output signal, t ₁, t ₂the time that starts and finish for prearranged signal;

The course of work of described unit failure unit is based on following formula:

$R\left(t\right)=S\left(t\right)\&CenterDot;\exp \{-\mathrm{\&lambda;}\&CenterDot;\underset{i}{\mathrm{\&Sigma;}}\underset{t_k\&le;t}{\mathrm{\&Sigma;}}{P}_i\left(t_k\right)\&times;\min [1,S\left(t_k\right)/S\left(t\right)\left]\right\}$

Wherein: R (t) is that output signal S (t) is input signal, and λ is element failure rate, P _i(t _k) be the value of i control signal when time t (k), S (t _k) be the value of input signal S (t) when time t (k).

8. a kind of reliability analysis system based on GO-FLOW method according to claim 1, is characterized in that,

The course of work of described normal closed gate analytic unit is based on following formula:

$\left\{\begin{array}{c}R\left(t\right)=S\left(t\right)\&CenterDot;O\left(t\right)\\ O\left(t_1\right)=P_p\\ O\left(t\right)=O\left(t^{\&prime;}\right)+[1-O\left(t^{\&prime;}\right)]\&CenterDot;P\left(t\right)\&CenterDot;P_g\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t') was the signal of O (t) in a upper moment; t ₁for initial time, P _pfor the probability that signal is not opened to valve, P by mistake _gfor the signal probability that then valve is successfully opened;

The course of work of described normally open valve door analytic unit is based on following formula:

$\left\{\begin{array}{c}R\left(t\right)=S\left(t\right)\&CenterDot;O\left(t\right)\\ O\left(t_1\right)=1-P_p\\ O\left(t\right)=O\left(t^{\&prime;}\right)[1-P\left(t\right)\&CenterDot;P_g]\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t') was the signal of O (t) in a upper moment; t ₁for initial time, P _pfor the probability that signal is not opened to valve, P by mistake _gfor the signal probability that then valve is successfully opened;

The described course of work of closing failsafe valve analytic unit is based on following formula:

$R\left(t\right)=S\left(t\right)\&CenterDot;\{1-\exp [-\mathrm{\&lambda;}\&CenterDot;\underset{i}{\mathrm{\&Sigma;}}\underset{t_k\&le;t}{\mathrm{\&Sigma;}}{P}_i\left(t_k\right)\left]\right\}$

Wherein: R (t) is that output signal S (t) is input signal, and λ is element failure rate, P _i(t _k) be the value of i control signal when time t (k);

The described course of work of opening failsafe valve analytic unit is based on following formula:

$R\left(t\right)=S\left(t\right)\&CenterDot;\exp [-\mathrm{\&lambda;}\&CenterDot;\underset{i}{\mathrm{\&Sigma;}}\underset{t_k\&le;t}{\mathrm{\&Sigma;}}{P}_i\left(t_k\right)\}$

Wherein: R (t) is that output signal S (t) is input signal, and λ is element failure rate, P _i(t _k) be the value of i control signal when time t (k)

The course of work of described valve switch motion analysis unit is based on following formula:

$R\left(t\right)=\left\{\begin{array}{c}S\left(t\right)\&CenterDot;\{O\left(t^{\&prime;}\right)+[1-O\left(t^{\&prime;}\right)]\&CenterDot;P_1\left(t\right)\&CenterDot;P_g\},P_1=1\\ S\left(t\right)\&CenterDot;O\left(t^{\&prime;}\right)[1-P_2\left(t\right)\&CenterDot;P_g],P_2=1\end{array}\right.$

Wherein: R (t) is that output signal S (t) is input signal, and O (t) is the signal of element internal, O (t') is O (t) at upper one signal constantly, P ₁for start signal, P ₂for shutdown signal, P _gfor opening or the shutdown signal probability that then valve is successfully opened or closed.

Images