CN106383960B - More electric aircraft electrical power system analysis method for reliability based on minimal cut set analytic approach - Google Patents
More electric aircraft electrical power system analysis method for reliability based on minimal cut set analytic approach Download PDFInfo
- Publication number
- CN106383960B CN106383960B CN201610861765.3A CN201610861765A CN106383960B CN 106383960 B CN106383960 B CN 106383960B CN 201610861765 A CN201610861765 A CN 201610861765A CN 106383960 B CN106383960 B CN 106383960B
- Authority
- CN
- China
- Prior art keywords
- node
- arc
- cut set
- minimal
- equivalent
- 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.)
- Expired - Fee Related
Links
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
More electric aircraft electrical power system analysis method for reliability that the present invention relates to a kind of based on minimal cut set analytic approach, comprising: topology diagram is equivalent to network node figure；Source node is sought using Adjacent Matrix Method and reaches specified node minimal path sets；Seek the minimal cut set of failure generation；Non cross link processing is carried out to minimal cut set using " leave out and leave " method；The equivalent fault rate that component physical fault rate is converted to node or arc, further acquires equivalent reliability.The present invention can more explicitly show the weak spot during the flowing of the aircraft energy, while provide improvement idea for aircraft energy operational version, also provide targeted scheme to the maintenance under aircraft difference run the period.
Description
Technical field
The present invention relates to the fields such as aerospace, power system security analysis and graph theory parsing.It can be generally applicable
In aircraft, the underway steamer, high-speed motion of the fail-safe analysis of the poor electric system of self-healing ability, such as mission phase
Electric car etc..It is suitable for reliability when analysis micro-grid system independent operating.
Background technique
More electricity aircrafts replace originally hydraulic, air pressure and mechanical drive system with electric system as much as possible, try hard to make
The secondary energy of aircraft uses electric energy more as far as possible.Electric system occupies leading role, reliable and stable fortune in more electric aircrafts
Row is related to the normal work of entire aircraft system.Aircraft safety accident caused by aircraft electrical power system failure happens occasionally, this
So that the reliability service of electric system caused it is more including aircraft manufacturing producer, airline, major research institution
Aspect concern.
Now both at home and abroad the open source literature for aircraft electrical power system reliability assessment simultaneously it is few, wherein most research at
Fruit all uses Monte Carlo Method or Fault Tree, and obtains relatively good achievement.
Monte Carlo Method is a kind of probability analysis method, has the characteristics that flexible in application, realization is simple.It can effectively be solved
Certainly fail-safe analysis problem.But when utilizing Monte Carlo Method to the fail-safe analysis of electric system, computation burden is big, efficiency
Low is always the inevitable disadvantage of Monte Carlo Method." the aviation based on model importance written by Lv Hong, Yuan Haiwen et al.
The reliability of power-supply system " introduces the concept of model importance, solves on partial extent merely using Monte Carlo method
The lower problem of efficiency.Model importance is significance level of the finger in subsystem, it is by calculating a certain component failure
Caused subsystem failure number accounts for the total Failure count ratio of subsystem and obtains.Model importance has distinguished normal elements and weakness
Impact effect of the component to system reliability.This also inspires us: in systems reliability analysis, treating different components and wants root
Corresponding adjustment is made to the reliability that it is converted according to physical unit action time and significance level.But lack in the method due to each
Interact amendment of the influence to reliability of generation between a subsystem.
Fault Tree be using the undesirable event of system as analysis purpose, by successively push down on trace back it is all
Possible reason, to enumerate component failure, environment influence, human error and various factors that may be present and thrashing
Logical relation, and formed one handstand fault tree.Zhou Suying master thesis " electrical power system of more electric aircraft research " with
It is analyzed and researched for advanced tactical fighter distribution system by reliability of the Fault Tree to its distribution system.In reality
The problem of huge calculation amount during fault tree can bring multiple shot array is established in, and fault tree is once established only
The reliability at some time point of system can be obtained, can not reflect the relationship of reliability and runing time well.But such as
Fruit will retrodict thought in fault tree, and using the knowledge of graph theory, available analysis well is imitated under computer aided calculation
Fruit.The logical source of this fail-safe analysis method namely of the invention.
It is general with modern aircraft in written by Cai Lin, Zhang Ling et al. " large aircraft power supply-distribution system reliability assessment and analysis "
Time power supply-distribution system structure give power supply-distribution system analysis method for reliability process, while by graph theory knowledge utilization
Computer aided calculation obtains aircraft components for the reliability of runing time, and analytical mathematics are perfectly in harmony with the present invention.
But this thinking also promotes improved space at several aspects: being on the one hand in the collection finding process of road using deep
First search is spent, causes road collection matrix that cannot reflect energy flow sequence very well, so that excessive system information is hidden；It is another
Aspect deeply excavate component in different times under dimension reliability variation.
Summary of the invention
The purpose of the present invention is improving to the prior art, a kind of more electric aircraft electrical power system fail-safe analysis sides are provided
Method shows the specific path of energy flow by adjacency Matrix Method, and different in length to aircraft overall operation reliability
Further mining analysis is made under time dimension.This improved method more explicitly show the aircraft energy flowing during it is thin
Weakness, while improvement idea is provided for aircraft energy operational version, needle also is provided with to the maintenance under aircraft difference run the period
To the scheme of property.Technical scheme is as follows:
A kind of more electric aircraft electrical power system analysis method for reliability based on minimal cut set analytic approach, include the following steps:
Step 1: according to following four principle, topology diagram is equivalent to network node figure:
1) network node figure needs to accurately reflect component relationship and working condition；
2) components such as breaker, conducting wire, transformer, rectifier and generator in former topological diagram are abstracted as network node
Arc in figure, wherein control energy flow direction is directed arc, it is otherwise undirected arc；
3) busbar in electric system and load are abstracted as the node in network node figure, while a virtual source is set
Node is equivalent to the source node of whole system power input；
4) network node figure is succinct as far as possible, and the concatenated component of institute is equivalent to an arc when no branch, type
The identical and load derived from a busbar is equivalent to a node；
Step 2: source node is sought using Adjacent Matrix Method and reaches specified node minimal path sets, the method is as follows:
1) meshed network N=(V) is set, wherein V={ v1, v2, v3, vnIt is node collection；
2) adjacency matrix A is defined1=[a1 ij], the topology information in network node figure is described, wherein the subscript of element name
1, which is expressed as path length, is 1 adjacency matrix, while also illustrating that first time Iterative Matrix in calculating process, and i, j respectively indicate square
The row and column of battle array:
3) terminal matrix R=[r is definedjk], reflect the terminal of every arc in meshed network N:
4) A is defined1Operation " ※ " with R acquires A2=[a2 ij], algorithm is as follows, a2 ij={ a1 ij※rjk|i,j,k
=1,2,3, n }:
5) virtual source node is sought to the minimal path sets of specified node using method as above, while minimal path sets being transformed into
Enter minimal path sets matrix and facilitates subsequent computer aided calculation；
Step 3: the minimal cut set of failure generation is sought:
Minimal path sets are observed, if each paths all pass through the same node or the same arc in minimal path sets, that
This node or arc break down, and mean that specified node or arc centainly will appear the situation of power loss, then this node or arc be just
Belong to the single order minimal cut set in grid；Node failure represented by single order minimal cut set is removed, other paths all include
Two components or wherein any one component, mean specified node if two components break down simultaneously
Or arc centainly will appear the situation of power loss, then this two components just belong to second order cut set；And so on found using the method
Source node is to each rank minimal cut set of specified node；That is any one ranks vector in minimal path sets matrix is all 1, then this
The corresponding component failure of column vector is exactly single order minimal cut set；Again by any two column vector in addition to single order minimal cut set
Inclusive-OR operation is carried out, if result all 1, obtains the minimum second order cut set that the corresponding component of the two column vectors is constituted, with
It is such to inquire into out all minimum second order cut sets, similarly obtain three rank cut sets, quadravalence cut set；
Step 4: non cross link processing is carried out to minimal cut set using " leave out and leave " method:
1) a n-dimensional vector E is defined for each minimal cut seti(x1, x2, x3···xn), xi0 or 1 is taken, x is worked asiIt takes
It indicates in minimal cut set to include this node or arc when obtaining 1, works as xiIndicate that minimal cut set does not include this section when obtaining 0
Point or arc；
2) companion matrix T is definedi(t1, t2, tn) and meet the following conditions:
3) E is set1Compare E when i=2 for the minimal cut set of first non cross linkiAnd TiFind TiMiddle all elements are
1 and correspond to EiThe position that middle corresponding position is 0 obtains k according to size order record position serial number and number1, k2, kr；
4) by EiTo kiIt is successively decomposed, is decomposed into Ei(k1) andWherein EiMiddle kiOn position 0 by 1 or-
1 replaces respectively；
5) compare Ei(k1) and Ej, wherein j < i, if all Ei(k1) in there is -1 position to correspond to EjPosition have one
A 1 or several 1 corresponding, and these 1 number summations are equal to line number j, obtain Ei(k1) and EjIt is incompatible, become non cross link
Minimal cut set, therefore by Ei(k1) storage leave, simultaneously forIdentical mode is done to handle；
6) E left is stored for noi(k1) and Ti-1Compare, if Ti-1In do not correspond to E for 0 positioni(k1) in phase
The position answered is not 0, then it is assumed that Ei(k1) by EjIt is absorbed, wherein j < i；Simultaneously forIt does at identical mode
Reason；
If 7) Ei(k1) do not leave and do not delete, it is continued to decompose, by Ei(k1) it is decomposed into Ei(k1)(k2) andWith respectively with 1, -1 instead of k20 on position, according to step 4,5,6 replicate analysis, until all EiIt is leaving
Or it makes a choice in deleting；
Step 5: component physical fault rate is converted to the equivalent fault rate of node or arc based on following principle, into one
Step acquires equivalent reliability；
1) for single component default its qualification factory after, service life should obey exponential law distribution, i.e., this
Component failure rate λ is steady state value, and assert that the faulty whole of institute is on component, and conducting wire connection is good, and tie point is not
Situations such as there are rosin joints；
2) for, by the concatenated arc of several components, equivalent process is carried out using following formula in network node figure
Conversion:
In formula, λeThe failure rate of the arc, λ after equivalentiSuccessively indicate the failure rate of each component in series connection, n indicates network
In node diagram, concatenated component counts in this arc；
3) individual node failure rate consists of three parts, and first is node or the failure rate λ of arc itselfb；Second not
In the case of protection device, expansion property failure rate caused by the failure that adjacent node or arc generate, by all in arc
Non-protected device, quantity are denoted as w, and wherein the expansion failure rate of some component k is denoted as λek；Third is by protecting device situation
When tripping occurs for lower and protection device, adjacent node or arc generate a possibility that being transmitted to the nodes break down after failure,
It is denoted as m by protective devices all in arc, and records before failure is isolated and have passed through several isolating devices or protective switch
Quantity is denoted as u, then the reliability of a node is calculated according to the following formula:
λ in formulaebFor the equivalent rear failure rate of individual node, PfiFor protective device failure transport=relay fail number/guarantor
Shield should move number；
4) failure rate of a node and arc is brought into non cross link treated in minimal cut set, using obtaining respective objects section
Point changes the failure rate of time:
P (S) indicates the destination node failure rate in formula, is equivalent to the reliability of corresponding system, p (Ei) indicate m cut set
In each cut set failure rate.
Detailed description of the invention
Fig. 1 Model in Reliability Evaluation of Power Systems algorithm flow chart
Fig. 2 cut set non cross link flow chart
The total electrical structure diagram of Fig. 3 A380
Fig. 4 A380 network node figure
Specific embodiment
The present invention will be described with reference to the accompanying drawings and examples.
The specific targets of fail-safe analysis are mainly reflected in the unit time on the trouble-proof probability of its aircraft.By general
Failure rate known to rate opinion knowledge and probability (abbreviation reliability) summation of reliability service are 1.
And the mathematical principle of path matrix will be sought using adjacency matrix in the graph theory being mainly concerned in patent and deleted
Going or staying laxative remedy is described below
Path matrix mathematical principle is sought using adjacency matrix:
Meshed network N=(V), wherein V={ v1, v2, v3, vnIt is node collection.Arc can there are two sections
Point indicates that direction is that previous node is directed toward the latter node.
Adjacency matrix A1=[a1 ij], describe the topology information in network.Wherein 1 it is expressed as the adjoining that path length is 1
Matrix, it is also assumed that being first time Iterative Matrix in calculating process.I, j respectively indicate the row and column of matrix.
Terminal matrix R=[rjk], matrix R reflects the terminal of every arc in meshed network N.
Define A1A kind of operation " ※ " with R acquires A2, algorithm such as a2 ij={ a1 ij※rjk| j=1,2,
3, n }.
Leave out the method for leaving:
1) a n-dimensional vector Ei (x1, x2, xn) is defined for each minimal cut set, xi can obtain 0,1.When
Xi indicates to include this node or arc in minimal cut set when obtaining 1, indicate that minimal cut set does not include when xi obtains 0
This node or arc.
2) it defines companion matrix Ti (t1t2, tn) and meets the following conditions
3) assume that E1 is the minimal cut set of first non cross link.When i=2, compares Ei and Ti and find all members in Ti
Element corresponds to the position that corresponding position in Ei is 0 for 1, obtains k1 according to size order record position serial number and number,
K2, kr.
4) successively decomposing according to K1 by Ei, be decomposed into Ei (k1) andThe position ki in Ei is replaced with 1, -1 respectively
On 0 obtain.
5) compare Ei (k1) and Ej (j < i), if the position for having -1 position to correspond to Ej in all Ei (k1) has one
A 1 or several 1 corresponding, and these 1 number summations are equal to line number j, available Ei (k1) is incompatible with Ej, becomes not
The minimal cut set of friendshipization.Therefore leave Ei (k1) storage, simultaneously forIdentical mode is done to handle.
6) Ei (k1) that leaves is stored and Ti-1 compares for no, if not being that 0 position corresponds to Ei (k1) in Ti-1
In corresponding position be not 0, it may be considered that Ei (k1) is absorbed by Ej (j < i).Simultaneously forIt is identical side
Formula processing.
7) it is continued to decompose if Ei (k1) does not leave and do not delete, by Ei (k1) be decomposed into Ei (k1) (k2) andWith 1, -1 replace on the position k2 0, according to step 4,5,6 replicate analysis, until all Ei leaving or
It makes a choice in deletion.
From an aircraft system structural topology figure is obtained to reliability is calculated in its operation to the relationship of time, need with
Lower five steps.Using more electricity aircraft A380 as the exploitativeness of case verification fail-safe analysis method.If Fig. 3 is that A380 is always electrically
System structure chart, and as five steps of example correspondence analysis, verify the exploitativeness of fail-safe analysis method.
1, according to following four principle, topological network figure is abstracted as network node figure:
1) network needs to accurately reflect the component relationship and working condition of former topological diagram.
2) switch, conducting wire, transformer, rectifier, generator etc. are abstracted as the arc in figure in electric system, wherein can control
Energy flow direction processed is directed arc.
3) busbar, load are abstracted as the node in figure in electric system, and a virtual source node is arranged, and are equivalent to entire
The source node of system power output.
4) network is succinct as far as possible, and the concatenated component of institute can be equivalent to an arc when no branch, type
The identical and load derived from a busbar can be equivalent to a node.
Change network node primitive according to structure chart and Fig. 3 is then changed into Fig. 4.
2, it seeks source node and reaches given load node minimal path sets
Using a kind of particular algorithm of adjacency matrix and terminal matrix, the minimal path sets of network are directly sought.Through excessive
A is sought in secondary calculating3, A4,···An-1, All Paths in network, i.e. required minimal path sets can be obtained.And by calculating process
In character matrix be converted to basic road collection matrix, to facilitate subsequent calculating.
3, the minimal cut set of failure generation is sought
Minimal path sets are observed, if road concentrates each paths all to pass through a node, this nodes break down
Then mean destination node centainly and will appear the situation of power loss, then this node element failure is exactly the single order minimal cut in network
Collection.Remove node failure represented by minimal cut set, two elements that other paths all include or wherein any one element,
Then the two element faults are then exactly second order cut set.And so on can find source node to each rank of destination node using the method
Minimal cut set.That is any one ranks vector in minimal path matrix is all 1, then this column for element be exactly single order most
Small cut set.Again by the corresponding progress inclusive-OR operation of any two column vector in addition to single order minimal cut set, if result is all
1, the minimum second order cut set of the two column vectors composition can be obtained, and so on find out all minimum second order cut sets, can similarly obtain most
Small three ranks cut set, quadravalence cut set.
4, the non cross link processing of minimal cut set
Using " leave out and the leave " method of by the agency of carries out non cross link processing to minimal cut set above.This method is mainly thought
Road is gradually melted into disjoint minimal cut set by increasing argument, facilitates the calculating of last reliability.As Fig. 2 process use is left out
Method is left to handle minimal cut set.
5, the conversion of reliability index
Since different node and arc are equivalent by different devices, wherein each device has the corresponding failure rate λ of oneself.And
And when breaking down, each device has interaction, will lead to corresponding failure rate λ and changes.Then being based on following principle will
Component physical fault rate is converted to the equivalent fault rate of node or arc, can further acquire equivalent reliability.
1) after defaulting its qualification factory for discrete component, service life should obey exponential law distribution, i.e. this yuan
Part failure rate λ is steady state value.And assert that the faulty whole of institute is on component, conducting wire connection is good, and tie point is not present
Situations such as rosin joint.
2) for, by the concatenated arc of several components, following formula progress can be used in equivalent process in network
Conversion.Wherein λiSuccessively indicate the failure rate of each component in series connection, n is indicated in network, concatenated member device in this arc
Part number.
3) individual node failure rate is made of three parts, and first is node or the failure rate λ of arc itselfb；Second not
In the case of protection device, expansion property failure rate caused by the failure that adjacent node or arc generate, by all in arc
Non-protected device, quantity are denoted as w, and wherein the expansion failure rate of some element k is denoted as λek；Third is by protecting device
And protection device is when occurring tripping, adjacent node or arc generate a possibility that being transmitted to the nodes break down after failure, warp
It crosses all protective devices in arc and is denoted as m, and record before failure is isolated and have passed through several isolating devices or protective switch number
Amount is denoted as u.Then the reliability of a node can be calculated according to formula
P in formulafiNumber should be moved for protective device failure transport=relay fail number/protection.
Finally the failure rate of each node and arc is brought into the minimal cut set of non cross link processing, passes through probability theory knowledge
Available respective objects node changes the failure rate of time.
Crucial airplane load was required in two hours in " plan of advanced system aviation electronics " that United States Air Force is write
Interior uninterrupted power supply reliability reaches 0.999 975 25.The important load in two hours in the calculating using A380 as example
Reliability reach 0.999 997 88, reach related request.
By instance analysis, the route redundancy for improving critical loads can significantly improve the reliability of long-play；
Reducing the number of components passed through before electric energy arrival critical loads can be improved short time reliability.This is to more electric electricity aircraft early period
Power design provides reference well and helps, and more feasible to one to the weak spot in the more electric aircraft electrical power systems of prediction
Method, while there is corresponding directive function to the optimization and upgrading of later period entire electric system.
Based on following precondition: disconnecting switch, channel paralleling switch have actuator control according to grid weight
Structure requires automatic closure, then power loss situation occur is the case where can not still making service restoration after failure reconfiguration.
Achieved reliability is finally switched into node reliability, calculated result such as table 1 by step 5
The reliability of table 1 specific node and arc
Obtained result is calculated and analyzes available relevant conclusion.Its method is equally applicable to the warship under operating status
The Calculation of Reliability of the not strong independent microgrid system of ship, electric car and self-healing ability.Need the not homology according to practical application
System flexibly holds the principle of network conversion；The physical significance of reasonable definition different rank minimal cut set；Appropriate transformation coefficient
Component physical fault rate is changed into node failure rate.And corresponding conclusion is obtained to interpretation of result, to microgrid operation and maintenance meter
It draws and advises, adjustable strategies appropriate can also be provided to microgrid structure updating.
Claims (1)
1. a kind of more electric aircraft electrical power system analysis method for reliability based on minimal cut set analytic approach, include the following steps:
Step 1: according to following four principle, topology diagram is equivalent to network node figure:
1) network node figure needs to accurately reflect component relationship and working condition；
2) breaker, conducting wire, transformer, rectifier and generator component in former topological diagram are abstracted as in network node figure
Arc, wherein control energy flow direction is directed arc, it is otherwise undirected arc；
3) busbar in electric system and load are abstracted as the node in network node figure, while a virtual source section is set
Point is equivalent to the source node of whole system power input；
4) network node figure is succinct as far as possible, and the concatenated component of institute is equivalent to an arc when no branch, and type is identical
And the load derived from a busbar is equivalent to a node；
Step 2: source node is sought using Adjacent Matrix Method and reaches specified node minimal path sets, the method is as follows:
1) meshed network N=(V) is set, wherein V={ v1, v2, v3, vnIt is node collection；
2) adjacency matrix A is defined1=[a1 ij], the topology information in network node figure is described, wherein the subscript 1 of element name indicates
The adjacency matrix for being 1 for path length, while also illustrating that first time Iterative Matrix in calculating process, i, j respectively indicate matrix
Row and column:
3) terminal matrix R=[r is definedjk], reflect the terminal of every arc in meshed network N:
4) A is defined1Operation " ※ " with R acquires A2=[a2 ij], algorithm is as follows, a2 ij={ a1 ij※rjk| i, j, k=1,
2,3, n }:
5) virtual source node is sought to the minimal path sets for specifying node using method as above, while minimal path sets conversion being entered most
Path collection matrix facilitates subsequent computer aided calculation；
Step 3: the minimal cut set of failure generation is sought:
Minimal path sets are observed, if each paths all pass through the same node or the same arc in minimal path sets, this
A node or arc break down, and mean that specified node or arc centainly will appear the situation of power loss, then this node or arc just belong to
Single order minimal cut set in grid；Remove single order minimal cut set represented by node failure, other paths all include two
A component or wherein any one component mean specified node or arc if two components break down simultaneously
It centainly will appear the situation of power loss, then this two components just belong to second order cut set；And so on using the method find source save
It puts to each rank minimal cut set of specified node；That is any one ranks vector in minimal path sets matrix is all 1, then this arrange to
Measuring corresponding component failure is exactly single order minimal cut set；Any two column vector in addition to single order minimal cut set is carried out again
Inclusive-OR operation obtains the minimum second order cut set that the corresponding component of the two column vectors is constituted, if result all 1 with such
Inquire into whole minimum second order cut sets out, similarly obtains three rank cut sets, quadravalence cut set；
Step 4: non cross link processing is carried out to minimal cut set using " leave out and leave " method:
1) a n-dimensional vector E is defined for each minimal cut seti(x1, x2, x3···xn), xi0 or 1 is taken, x is worked asiObtain 1
When indicate minimal cut set in include this node or arc, work as xiWhen obtaining 0 indicate minimal cut set do not include this node or
Person's arc；
2) companion matrix T is definedi(t1, t2, tn) and meet the following conditions:
3) E is set1Compare E when i=2 for the minimal cut set of first non cross linkiAnd TiFind TiMiddle all elements are right for 1
Answer EiThe position that middle corresponding position is 0 obtains k according to size order record position serial number and number1, k2, kr；
4) by EiTo kiIt is successively decomposed, is decomposed into Ei(k1) andWherein EiMiddle kiOn position 0 by 1 or -1 point
Not instead of；
5) compare Ei(k1) and Ej, wherein j < i, if all Ei(k1) in there is -1 position to correspond to EjPosition have one 1
Or it is several 1 corresponding, and these 1 number summations are equal to line number j, obtain Ei(k1) and EjIt is incompatible, become non cross link most
Small cut set, therefore by Ei(k1) storage leave, simultaneously forIdentical mode is done to handle；
6) E left is stored for noi(k1) and Ti-1Compare, if Ti-1In do not correspond to E for 0 positioni(k1) in it is corresponding
Position is not 0, then it is assumed that Ei(k1) by EjIt is absorbed, wherein j < i；Simultaneously forIdentical mode is done to handle；
If 7) Ei(k1) do not leave and do not delete, it is continued to decompose, by Ei(k1) it is decomposed into Ei(k1)(k2) andK is replaced with 1, -1 respectively20 on position, foundation step 4), 5), 6) and replicate analysis, until all EiIt is staying
It makes a choice in lower or deletion；
Step 5: component physical fault rate is converted to the equivalent fault rate of node or arc based on following principle, is further asked
Obtain equivalent reliability:
1) after defaulting its qualification factory for single component, service life should obey exponential law distribution, i.e. this yuan of device
Part failure rate λ is steady state value, and assert that the faulty whole of institute is on component, and conducting wire connection is good, and tie point is not present
Situations such as rosin joint；
2) for, by the concatenated arc of several components, equivalent process is rolled over using following formula in network node figure
It calculates:
In formula, λeThe failure rate of the arc, λ after equivalentiSuccessively indicate the failure rate of each component in series connection, n indicates network node
In figure, concatenated component counts in this arc；
3) individual node failure rate consists of three parts, and first is node or the failure rate λ of arc itselfb；Second without guarantor
It protects in the case of device, expansion property failure rate caused by the failure that adjacent node or arc generate, by all non-protected in arc
Device, quantity are denoted as w, and wherein the expansion failure rate of some component k is denoted as λek；Third is by protecting device and protects
When protecting device tripping occurs, adjacent node or arc generate a possibility that being transmitted to the nodes break down after failure, by arc
In all protective devices be denoted as m, and record and have passed through several isolating devices or protective switch quantity note before failure is isolated
For u, then the reliability of a node is calculated according to the following formula:
λ in formulaebFor the equivalent rear failure rate of individual node, PfiIt is answered for protective device failure transport=relay fail number/protection
Dynamic number；
4) failure rate of each node and arc is brought into non cross link treated in minimal cut set, using obtaining respective objects section
Point changes the failure rate of time:
P (S) indicates the destination node failure rate in formula, is equivalent to the reliability of corresponding system, p (Ei) indicate each in m cut set
The failure rate of cut set.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610861765.3A CN106383960B (en) | 2016-09-28 | 2016-09-28 | More electric aircraft electrical power system analysis method for reliability based on minimal cut set analytic approach |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610861765.3A CN106383960B (en) | 2016-09-28 | 2016-09-28 | More electric aircraft electrical power system analysis method for reliability based on minimal cut set analytic approach |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106383960A CN106383960A (en) | 2017-02-08 |
| CN106383960B true CN106383960B (en) | 2019-07-02 |
Family
ID=57937068
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610861765.3A Expired - Fee Related CN106383960B (en) | 2016-09-28 | 2016-09-28 | More electric aircraft electrical power system analysis method for reliability based on minimal cut set analytic approach |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106383960B (en) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106980704B (en) * | 2017-02-13 | 2020-07-14 | 上海交通大学 | Multi-target transfer strategy flexible optimization method applied to power failure load of multi-electric aircraft |
| CN107944705B (en) * | 2017-11-27 | 2022-03-11 | 国网河南省电力公司信息通信公司 | Full-end reliability calculation method for dividing communication communities based on modularity |
| CN110555223A (en) * | 2018-06-04 | 2019-12-10 | 北京交通大学 | High-speed train system reliability assessment method based on seepage theory |
| CN110009240B (en) * | 2019-04-12 | 2022-03-29 | 西南交通大学 | Power system reliability assessment method |
| CN111726254B (en) * | 2020-06-22 | 2022-05-24 | 中国人民解放军国防科技大学 | All-demand-oriented minimum path set vector repetition-free search algorithm for multi-state network |
| CN113467984B (en) * | 2021-07-20 | 2024-03-29 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | Method, device, computer equipment and storage medium for diagnosing faults of aircraft system |
| CN115204761A (en) * | 2022-09-19 | 2022-10-18 | 中国兵器科学研究院 | Failure evaluation method, system and medium for multi-node delivery network |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6125453A (en) * | 1998-06-30 | 2000-09-26 | Sandia Corporation | Cut set-based risk and reliability analysis for arbitrarily interconnected networks |
| CN101651343A (en) * | 2009-09-09 | 2010-02-17 | 国家电网公司 | Method and system for checking electric power system model based on hybrid dynamic simulation |
| CN102751723A (en) * | 2012-06-19 | 2012-10-24 | 江苏科技大学 | Ship alternative-current zone power distribution network system and reliability analysis method thereof |
| CN102999788A (en) * | 2012-11-06 | 2013-03-27 | 天津大学 | Fast search method for active splitting surfaces for power system on the basis of multilayer graph partitioning |
| CN103633647A (en) * | 2013-12-05 | 2014-03-12 | 国家电网公司 | Power grid topology based power system reliability calculation method |
| CN103824234A (en) * | 2014-03-18 | 2014-05-28 | 国家电网公司 | Blocking and hierarchical structure based power distribution system reliability evaluation method |
| CN105470949A (en) * | 2015-12-01 | 2016-04-06 | 江苏省电力公司淮安供电公司 | Power distribution network reliability evaluation method considering complex transfer supply |
| CN106355510A (en) * | 2015-07-15 | 2017-01-25 | 中国电力科学研究院 | Power system safety determination method |
-
2016
- 2016-09-28 CN CN201610861765.3A patent/CN106383960B/en not_active Expired - Fee Related
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6125453A (en) * | 1998-06-30 | 2000-09-26 | Sandia Corporation | Cut set-based risk and reliability analysis for arbitrarily interconnected networks |
| CN101651343A (en) * | 2009-09-09 | 2010-02-17 | 国家电网公司 | Method and system for checking electric power system model based on hybrid dynamic simulation |
| CN102751723A (en) * | 2012-06-19 | 2012-10-24 | 江苏科技大学 | Ship alternative-current zone power distribution network system and reliability analysis method thereof |
| CN102999788A (en) * | 2012-11-06 | 2013-03-27 | 天津大学 | Fast search method for active splitting surfaces for power system on the basis of multilayer graph partitioning |
| CN103633647A (en) * | 2013-12-05 | 2014-03-12 | 国家电网公司 | Power grid topology based power system reliability calculation method |
| CN103824234A (en) * | 2014-03-18 | 2014-05-28 | 国家电网公司 | Blocking and hierarchical structure based power distribution system reliability evaluation method |
| CN106355510A (en) * | 2015-07-15 | 2017-01-25 | 中国电力科学研究院 | Power system safety determination method |
| CN105470949A (en) * | 2015-12-01 | 2016-04-06 | 江苏省电力公司淮安供电公司 | Power distribution network reliability evaluation method considering complex transfer supply |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106383960A (en) | 2017-02-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106383960B (en) | More electric aircraft electrical power system analysis method for reliability based on minimal cut set analytic approach | |
| CN107609725B (en) | Power transmission network risk assessment method considering transformer substation influence | |
| CN106295956B (en) | Reliability assessment method considering repairable multi-state complex characteristic of power supply system of nuclear power plant | |
| Schneider et al. | Assessment of interactions between power and telecommunications infrastructures | |
| Srivastava et al. | Expert-system method for automatic reconfiguration for restoration of shipboard power systems | |
| Azeroual et al. | Fault location and detection techniques in power distribution systems with distributed generation: Kenitra City (Morocco) as a case study | |
| CN109977541A (en) | A kind of power information physics system union simulation platform and its synchronous method | |
| Wang et al. | An efficient algorithm for assessing reliability indexes of general distribution systems | |
| Su et al. | Performance evaluation of multiagent systems for navy shipboard power system restoration | |
| CN109064071A (en) | A kind of secondary system of intelligent substation methods of risk assessment based on shielding failure | |
| Zhao et al. | Minimal Cut Sets‐Based Reliability Evaluation of the More Electric Aircraft Power System | |
| Ling et al. | Reliability analysis of DC power distribution network based on minimal cut sets | |
| Zin et al. | Protection system analysis using fault signatures in Malaysia | |
| Belkacemi et al. | Experimental implementation of multi-agent system algorithm for distributed restoration of a smart grid system | |
| Zhang et al. | A novel fault isolation scheme in power system with dynamic topology using wide-area information | |
| Qin et al. | A multiple model filtering approach to transmission line fault diagnosis | |
| Vega et al. | Algorithm to evaluate substation reliability with cut and path sets | |
| Khare et al. | Optimal PMU placement using matrix modification based integer linear programming | |
| Singh et al. | Reliability measures and sensitivity analysis of a complex matrix system including power failure | |
| Billinton et al. | Generalized n+ 2 state system Markov model for station-oriented reliability evaluation | |
| Xiao et al. | Reliability evaluation of protection system in smart substation based on process layer network | |
| Krishna et al. | Complete and incomplete observability analysis by optimal PMU placement techniques of a network | |
| Kang et al. | Unavailability and spurious operation probability of k-out-of-n reactor protection systems in consideration of CCF | |
| Miao et al. | Power supply reliability indices computation with consideration of generation systems, transmission systems and sub-transmission systems' load transfer capabilities | |
| Sobolevskis et al. | Vulnerability assessment of electric power system for the case of Latvian EPS |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190702 Termination date: 20200928 |