CN104063566B - Under the influence of a kind of determination Binary Factor in electrical system element significance level method - Google Patents

Under the influence of a kind of determination Binary Factor in electrical system element significance level method Download PDF

Info

Publication number
CN104063566B
CN104063566B CN201310087838.4A CN201310087838A CN104063566B CN 104063566 B CN104063566 B CN 104063566B CN 201310087838 A CN201310087838 A CN 201310087838A CN 104063566 B CN104063566 B CN 104063566B
Authority
CN
China
Prior art keywords
influence
significance level
probability
reliability
electrical system
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
CN201310087838.4A
Other languages
Chinese (zh)
Other versions
CN104063566A (en
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.)
Liaoning Technical University
Original Assignee
Liaoning Technical 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 Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN201310087838.4A priority Critical patent/CN104063566B/en
Publication of CN104063566A publication Critical patent/CN104063566A/en
Application granted granted Critical
Publication of CN104063566B publication Critical patent/CN104063566B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Testing Electric Properties And Detecting Electric Faults (AREA)
  • Tests Of Electronic Circuits (AREA)

Abstract

The invention discloses the method for element significance level in electrical system under the influence of a kind of determination Binary Factor, feature is just to influence two key factors of electrical equipment reliability:Working time(t)And operating temperature(c)Reliability to discrete component is analyzed, when the working time of particularly each element and different appropriate working temperature, it is difficult using traditional method analysis system reliability, the present invention is based on the influence degree to system reliability, the significance level of research wherein element.The present invention can use accident tree representation system architecture, and abbreviation is carried out to system architecture by accident tree, obtain considering the system fault probability distribution under the influence of t and c Binary Factors, significance level sequence is carried out to element based on probability space importance.It can be widely used under multifactor impact, determine significance level of the element to system.

Description

Under the influence of a kind of determination Binary Factor in electrical system element significance level method
Technical field
The present invention relates to electrical system reliability, more particularly to using electrical system under the influence of a kind of determination Binary Factor The method of middle element significance level.
Background technology
Electrical system is most common system in present every field, the entirety of system where its reliability is directly affected Performance.From system perspective analysis, its reliability can be divided into two parts and be studied.One is the primary element of composition system, this The property of a little elements is applied to the reliability of itself, and then influences the reliability of this electrical system.Two be the knot of system in itself Structure, is exactly the building form of primary element, the different effect journeys that will directly determine elements affect system reliability of building form Degree.The reliability of whole system is both combinations.
For the diode element in electrical system, its probability of malfunction length just with the working time, operating temperature Size, by electric current and voltage etc. there is direct relation.Assuming that the system failure is due to caused by component wear, and by changing member Part carries out failture evacuation.So use time of element is by the key factor as influence component reliability, and this factor influences The degree of probability of malfunction obeys exponential expression.Another factor is exactly operating temperature, it will be apparent that, for electrical equipment temperature mistake The rising of the high and too low decline and fault rate that can all cause its reliability, obeys cosine curve substantially.Electrically member is built first The probability of malfunction space based on use time (t) He operating temperature (c) of part, and it is made up of these elements the time of system (t) with the probability of malfunction space of operating temperature (c).Determine that element is reliable to system in system by probabilistic compct spatial distribution Property effect.
The content of the invention
For preferably invention is described, simple electrical system is designed here and is described, the system is by diode Composition, the rated operation of diode is affected by many factors, wherein importantly t and c.Herein for by the two factors The electrical system of influence is used as research object.There are five primary element X in system1、X2、X3、X4、X5, and it is set to had bright by t and c Develop loud element, shown in its classical accident tree Fig. 1.The accident tree abbreviation of the system is obtained:T=X1X2X3+X1X4+X3X5
1. the fail-safe analysis of electrical equipment
5 essential electronic element X in system1、X2、X3、X4、X5Probability of malfunction, be all to be influenceed by t and c, i.e., it is first The probability of malfunction P of parti(t, c), wherein i ∈ { 1,2,3,4,5 } are the function of t and c as independent variable with.As t and the sides of c two Element just breaks down during the failure of one of face, according to logic or concept Pi(t, c) is as shown in formula (1).
Determine Pi(t, c), it is necessary to first determineWithIf can not be repaiied after discrete component breaks down in system, it is System, which is fixed a breakdown, to be realized by changing element.ThenIt may be considered the cell failure probability of not repairable system[8], and If failure reach 0.9999 element should change (this data can be obtained by given system fault rate back analysis, generally ratio This value is much smaller), as shown in Equation 2.
In formula:λ is cell failure rate.
ForThe normal work of electrical equipment will have certain operating temperature range, higher or lower than the temperature Range element just breaks down, and the rule is expressed as into cosine curve herein, as shown in Equation 3.
In formula:A is range of temperature.
The element of actually distinct type has the scope of different use time life-span and appropriate working temperature
2. the fail-safe analysis of electrical system
Obtained by Fig. 1 systematic failures tree abbreviations, formula (4) is as follows:
T=X1X2X3+X1X4+X3X5 (4)
The system failure (top event) probability of happening is obtained by classical accident tree theory, as shown in formula (5):
PT(t, c)=P1P2P3+P1P4+P3P5-P1P2P3P4-P1P3P4P5-P1P2P3P5+P1P2P3P4P5(5)
From formula (5), PT(t, c) is the function for reflecting electrical malfunction probability, and the function is by P1~5(t, c) is determined, Again by formula (1), it is known that PT(t, c) be byWithThat is PT(t, c) is the function by t and c, by PT(t, c), t and c The three-dimensional probability space distribution of composition and its equivalent curve.
3. the spatial distribution and element importance ranking method of probabilistic compct
Probabilistic compct spatial distribution:The change of i-th of elementary event probability of happening causes top event probability of happening to become The degree of change, in the case of n dimension influence factor changes, the spatial distribution showed in n+1 dimension spaces.WithRepresent, the probabilistic compct spatial distribution of such as the 1st element herein is:
Brief description of the drawings
The accident tree of Fig. 1 electrical systems
Fig. 2 X1~5Probability of malfunction spatial distribution and its equivalent curve
Fig. 3 X1Probabilistic compct spatial distribution
Fig. 4 X1~5Maximum probability importance is distributed
Embodiment
Embodiment is the electrical system shown in Fig. 1.
5 essential electronic element X in system1、X2、X3、X4、X5Probability of malfunction, be all to be influenceed by t and c, i.e., it is first The probability of malfunction P of parti(t, c), wherein i ∈ { 1,2,3,4,5 } are the function of t and c as independent variable with.As t and the sides of c two Element just breaks down during the failure of one of face, according to logic or concept Pi(t, c) such as following formula:
Determine Pi(t, c), it is necessary to first determineWithIf can not be repaiied after discrete component breaks down in system, it is System, which is fixed a breakdown, to be realized by changing element.ThenIt may be considered the cell failure probability of not repairable system[8], and If failure reach 0.9999 element should change (this data can be obtained by given system fault rate back analysis, generally ratio This value is much smaller), as shown in Equation 2.ForThe normal work of electrical equipment will have certain operating temperature range, Just broken down higher or lower than the temperature range element, the rule is expressed as cosine curve herein, as shown in Equation 3.
In formula:λ is cell failure rate, and A is range of temperature.
The element of actually distinct type has the scope of different use time life-span and appropriate working temperature, assumes herein Their use scope, working time scope t ∈ [0,100] day of research, operating temperature interval c ∈ [0,50] DEG C.And root Calculate and obtain according to formula (2) and formula (3)WithExpression functional relation in the range of each.WithRespective It is not continuous, but piecewise function in research range.The segmentation of each function represents as shown in table 1.
System element X can be constructed by table 2 and formula (1)1~5Probability of malfunction spatial distribution and its equivalent curve, such as Fig. 2 It is shown.
Table 1WithExpression formula in survey region
In Fig. 2, X1~5Probability of malfunction spatial distribution and its equivalent curve be all different, this is due to its t and c What influence was caused.For working time t in the search time region of each element, there are two in probability of malfunction spatial distribution map It is due to that element reaches that new element is changed during probability of malfunction 0.9999 to be caused or trizonal probability of malfunction is substantially reduced. Probability of malfunction when actually this is changed can use polynary accident tree space to manage by setting the probability of malfunction of whole system Obtained by inverting, it is much smaller actually to calculate obtained probability of malfunction, in view of herein and this method is not used, and this example is only For illustrating, it is not described in detail here.For operating temperature c, due to using cosine curve as representative function, probability of malfunction Minimum position is in the middle of adaptive temperature scope.From image, the less position of element fault probability concentrates on temperature The intermediate region of scope.But, element accident probability acceptable scope is less on the diagram, and this is due to use binary The inevitable outcome of accident tree representation element fault probability.The superposition of two probability adds element total breakdown probability, this Phenomenon can not be analyzed using classical accident tree.Certainly, the reason for also having element replacement excessive cycle.
Obtained by Fig. 1 systematic failures tree abbreviations, formula (4) is as follows:
T=X1X2X3+X1X4+X3X5 (4)
The system failure (top event) probability of happening is obtained by classical accident tree theory, as shown in formula (5):
PT(t, c)=P1P2P3+P1P4+P3P5-P1P2P3P4-P1P3P4P5-P1P2P3P5+P1P2P3P4P5(5)
From formula (5), PT(t, c) is the function for reflecting electrical malfunction probability, and the function is by P1~5(t, c) is determined, Again by formula (1), it is known that PT(t, c) be byWithThat is PT(t, c) is the function by t and c, by PT(t, c), t and c The three-dimensional probability space distribution of composition and its equivalent curve are as shown in Figure 3.
It was found from Fig. 3 left figures, system fault probability is minimum near the t=0 moment, and main cause is all elements in system Enter use state simultaneously at the t=0 moment, the probability of malfunction of each element this period is all very low, makes the failure of whole system Probability is reduced.In terms of temperature in use, the temperature in use of majority element is all at 20 DEG C to 30 DEG C, so system is in this humidity province Between the probability of malfunction that works it is relatively low.But development over time, the probability of malfunction of element constantly increases, begin with element by for Change, while other elements also maintain original probability of malfunction curvilinear trend to continue to develop, make the new element of replacing to the system failure The effect that probability reduces is cancelled.The ability that each replacement of element cycle difference causes new element to improve system reliability is mutually supported Disappear, make the system failure rate in other regions in addition near t=0 very high.Fig. 3 right figures can be seen that each probability of malfunction forms isolated island, In addition to the characteristics of being analyzed above, the center in temperature of each isolated island is not consistent, and this also reflects has changed member at the moment The adaptive temperature scope of part and these elements is all different.
The electrical system according to caused by defining the change of the i-th element fault probability of happening breaks down the intensity of variation of probability, Probabilistic compct spatial distribution is constituted under the influence of t and c factors, it is the important references of analysis element and system change relation One of, MTLαThe optimization of calculating process will investigate X1~5The probabilistic compct spatial distribution of element.
According to definition for the probabilistic compct spatial distribution such as formula (6) of element.
With regard to X1Probabilistic compct spatial distribution it is following (other are not listed):
By Ig(1) I is arrivedg(5), deploy in t and c three dimensions, formation probability importance spatial distribution is illustrated in figure 3 X1Distribution map.
It is distributed according to Fig. 3 probabilistic compct, the X in gamut is studied here1~5Probabilistic compct sequence.It is aobvious The size of right probabilistic compct determines that its sorts, but the importance sorting in whole survey region be not it is consistent, such as Shown in Fig. 4.
Color corresponding 1,2,3,4,5 represents X respectively in Fig. 41、X2、X3、X4、X5Maximum probability is in some region The distribution of importance.So for whole survey region, taking the method being integrated to curved surface, that is, calculate probabilistic compct point Cloth curved surface and the curved surface of probabilistic compct=0, at 20 ° to 30 °, the volume (X in the range of 0 to 100 days1=224.2744, X2 =8.3744, X3=174.4662, X4=120.7763, X5=94.8089), obtained probabilistic compct is ordered as X1> X3> X4> X5> X2.This result and Fig. 4 result are not quite identical, are due to the size of probabilistic compct not in accordance with occupying What how much region determined, but what integrated height and the volume size of region formation were determined.

Claims (7)

1. a kind of method for determining element significance level in electrical system under the influence of Binary Factor, it is characterised in that just influence electricity Two key factors of gas component reliability:Working time t and operating temperature c are analyzed the reliability of discrete component, often The use of traditional method analysis system reliability is difficult when the working time of individual element and all different appropriate working temperature , based on the influence degree to system reliability, the significance level of research wherein element;It comprises the following steps:Electrical equipment Reliability determine, the reliability of electrical system is determined, the spatial distribution of probabilistic compct and element importance ranking method, can With accident tree representation system architecture, abbreviation is carried out to system architecture by accident tree, under the influence of obtaining considering t and c Binary Factors System fault probability distribution, based on probability space importance to element carry out significance level sequence.
2. a kind of method for determining element significance level in electrical system under the influence of Binary Factor according to claim 1, Characterized in that, the probability of malfunction P of elementi(t, c) is the function of t and c as independent variable, when one of t and the aspects of c two failure Element just breaks down, according to logic or concept Pi(t, c) is:
Pi(t, c)=1- (1-Pi t(t))(1-Pi c(c))
Determine Pi(t, c), it is necessary to first determine Pi tAnd P (t)i c(c), if can not be repaiied after discrete component breaks down in system, system Fix a breakdown is realized by changing element.
3. a kind of method for determining element significance level in electrical system under the influence of Binary Factor according to claim 1, Characterized in that, Pi t(t) the cell failure probability of not repairable system is may be considered, and sets failure and reaches that 0.9999 element should Change, be shown below:
Pi t(t)=0.9999=1-e-λt;λ t=9.2103
In formula:λ is cell failure rate.
4. a kind of method for determining element significance level in electrical system under the influence of Binary Factor according to claim 1, Characterized in that, for Pi c(c), the normal work of electrical equipment will have certain operating temperature range, higher or lower than this Temperature range element just breaks down, Pi c(c) cosine curve is expressed as, is shown below:
<mrow> <msubsup> <mi>P</mi> <mi>i</mi> <mi>c</mi> </msubsup> <mrow> <mo>(</mo> <mi>c</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mn>2</mn> <mi>&amp;pi;</mi> <mi>c</mi> <mo>/</mo> <mi>A</mi> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> <mn>2</mn> </mfrac> </mrow>
In formula:A is range of temperature.
5. a kind of method for determining element significance level in electrical system under the influence of Binary Factor according to claim 1, Characterized in that, probabilistic compct spatial distribution:The change of i-th of elementary event probability of happening causes top event probability of happening The degree of change, in the case of n dimension influence factor changes, the spatial distribution showed in n+1 dimension spaces is usedRepresent, the probabilistic compct spatial distribution of such as the 1st element herein is:
Wherein:PT(t, c) represents system fault probability distribution;Pi(t, c) represents element i Failure probability distribution.
6. a kind of method for determining element significance level in electrical system under the influence of Binary Factor according to claim 1, Characterized in that, probabilistic compct is distributed, the X in gamut is studied here1~5Probabilistic compct sequence, it is clear that probability The size of importance determines that it sorts, but the importance sorting in whole survey region is not consistent.
7. a kind of method for determining element significance level in electrical system under the influence of Binary Factor according to claim 1, its It is characterised by, so for whole survey region, taking the method being integrated to curved surface, that is, calculates probabilistic compct distribution surface With the curved surface of probabilistic compct=0, in 20 ° to 30 °, the volume in the range of 0 to 100 days, X1=224.2744, X2=8.3744, X3=174.4662, X4=120.7763, X5=94.8089, obtained probabilistic compct is ordered as X1> X3> X4> X5> X2
CN201310087838.4A 2013-03-19 2013-03-19 Under the influence of a kind of determination Binary Factor in electrical system element significance level method Active CN104063566B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310087838.4A CN104063566B (en) 2013-03-19 2013-03-19 Under the influence of a kind of determination Binary Factor in electrical system element significance level method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310087838.4A CN104063566B (en) 2013-03-19 2013-03-19 Under the influence of a kind of determination Binary Factor in electrical system element significance level method

Publications (2)

Publication Number Publication Date
CN104063566A CN104063566A (en) 2014-09-24
CN104063566B true CN104063566B (en) 2017-10-31

Family

ID=51551278

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310087838.4A Active CN104063566B (en) 2013-03-19 2013-03-19 Under the influence of a kind of determination Binary Factor in electrical system element significance level method

Country Status (1)

Country Link
CN (1) CN104063566B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105243264A (en) * 2015-09-14 2016-01-13 辽宁工程技术大学 Method for determining component failure probability
CN106055856B (en) * 2015-09-15 2018-08-31 辽宁工程技术大学 A kind of element fault probability based on discrete data determines method
CN105718737A (en) * 2016-01-22 2016-06-29 辽宁工程技术大学 Method for determining clouding factor importance degree and factor joint importance degree
US11416326B2 (en) * 2020-08-28 2022-08-16 Sap Se Systems and methods for failure diagnosis using fault tree

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102436519A (en) * 2011-08-23 2012-05-02 戴志辉 Method for synthetically evaluating dynamic reliability of power system automatic device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6794879B2 (en) * 2000-11-08 2004-09-21 General Electric Company Apparatus and method for detecting and calculating ground fault resistance

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102436519A (en) * 2011-08-23 2012-05-02 戴志辉 Method for synthetically evaluating dynamic reliability of power system automatic device

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于故障树和神经网络的飞机电源系统故障诊断;薄海涛;《创意与实践》;20051231;第4卷;第65-67页 *
计算机辅助故障树分析方法研究与应用;杨建强;《中国优秀硕士学位论文全文数据库 信息科技辑》;20060415(第4期);第5-16页 *

Also Published As

Publication number Publication date
CN104063566A (en) 2014-09-24

Similar Documents

Publication Publication Date Title
CN104063566B (en) Under the influence of a kind of determination Binary Factor in electrical system element significance level method
Schüler et al. Optimal Hubbard Models for Materials with Nonlocal Coulomb Interactions:<? format?> Graphene, Silicene, and Benzene
US20050015740A1 (en) Design for manufacturability
CN102436529B (en) Modeling method for metal oxide semiconductor (MOS) transistor reliability statistics models distributed based on Weibull
US20120025403A1 (en) Design apparatus of semiconductor device, design method of semiconductor device, and semiconductor device
CN107491643A (en) A kind of method, apparatus and system for determining synthesis design flow
EP3175401A1 (en) Predicting and optimizing energy storage lifetime performance with adaptive automation control software
CN103258110B (en) Method for determining accident trend of electrical system on basis of states
CN110222313A (en) Reflect the drought early warning method and apparatus of nonuniformity drought character variable
CN110350540A (en) A kind of fine cutting load method based on On-line Estimation frequency character of load
CN110008076A (en) Processor reliability estimation method, device, electronic equipment and storage medium
CN103955551A (en) Scheme for realizing reliability cost minimization of electrical system
CN103838897A (en) Layered antenna checking method of integrated circuit layout graph verification
CN102694395B (en) Method for evaluating reliability of two terminals of power transmission network
CN103077278B (en) A kind of voltage island energization pins distribution method of SOC (system on a chip)
US6925626B2 (en) Method of routing a redistribution layer trace in an integrated circuit die
CN104239590A (en) Graph structure adjusting method in integrated circuit layout verification
CN103969519A (en) Scheme guaranteeing reliability of electrical system
Zhong et al. Whitespace insertion for through-silicon via planning on 3-D SoCs
CN107863771B (en) Multimode electric system redundancy optimization method based on sequence optimization and Monte Carlo
CN107834538B (en) Multimode electric system redundancy optimization method based on sequence optimization and markov chain
CN104573146A (en) Clock signal transmission adjusting method and related integrated circuit structure
CN106815429B (en) Circuit board lamination deployment method and device
CN104579963A (en) Method and device for optimizing routes of nodes
CN104363293A (en) Web service selection method based on variable coefficients

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
CP02 Change in the address of a patent holder
CP02 Change in the address of a patent holder

Address after: 125105 Longwan South Street, Huludao, Huludao, Liaoning

Patentee after: Liaoning Technical University

Address before: No. 47, Zhonghua Road, Fuxin City, Liaoning Province, Liaoning

Patentee before: Liaoning Technical University