2. Teil
Prof.ErikaHausenblas
MontanuniversitätLeoben,Österreih
10.Jänner2017
1
Zuverlässigkeit gröÿerer Systeme
Komponenten in Seriegeshaltet
Redundanzoder Komponenten parallel geshaltet
Komplexe Kongurationen
MinimalCutSets undPaths
2
Vershiedenes
Standby Shaltungen (ReserveShaltungen)
3
Modellierung von Abhängigkeiten
Betrahten wirein Systemdasaus
n
KomponentenE i
,i =
1, . . . , n
inSeriezusammengeshaltetist.
R i (t)
:Zuverlässigkeitsfunktion einesUntersystemsAnnahme: dieKomponenten funktionierenunabhängig voneinander;
Zuverlässigkeitsfunktion
R S
desGesamtsystemsR S (t) =
n
Y
i=
1R i (t ).
Angenommen:
λ i (t )
ist Ausfallratederi
-ten Komponente zurZeitt >
0Damit gilt
R i (t) =
exp(−λ i (t))
,i =
1, . . . , n
.und somitR S (t ) =
n
Y
i=
1R i (t) =
exp−
n
X
i =
1λ i (t)
!
.
Example
Einuss derFunktionstühtigkeiteinzelnerKomponentenaufdasGesamtsystem:
Zuverl.
derKomp. 0.8 0.85 0.9 0.95 0.98 0.99
Anzahl
derKomp.
1 0.8 0.85 0.9 0.95 0.98 0.99
5 0.32768 0.44370 0.59049 0.77378 0.90392 0.95099
10 0.10737 0.19687 0.348678 0.598737 0.817073 0.904382
20 0.0115292 0.0387595 0.12158 0.35849 0.66761 0.81791
50 1.47
·
10−
5 2.96·
10−
4 5.15·
10−
3 0.07694 0.36417 0.60501Betrahten wirein Systemdasaus
n
KomponentenE i
,i =
1, . . . , n
inSeriezusammengeshaltetist.
R i (t)
:Zuverlässigkeitsfunktion einesUntersystemsAnnahme: dieKomponenten funktionierenunabhängig voneinander;
Hier giltfür
F S
F S (t) =
n
Y
i=
1F i (t),
woraus folgt
R S (t) =
1− F S (t) =
1−
n
Y
i=
1(
1− R i (t)) =
1−
n
Y
i=
1(
1− e − λ i (t) )
Example
Einuss derFunktionstühtigkeiteinzelnerKomponentenaufdasGesamtsystem:
Anzahl Zuverl.des Anstiegder
derKomp. Systems Zuverl. inProzent
1 0.8
· · ·
2 0.960000 0.160000 20.00
3 0.992000 0.032000 24.00
4 0.998400 0.006400 24.80
5 0.999680 0.001280 24.96
6 0.999936 0.000256 24.99
Example
EinBauteilhateinekonstanteAusfallrate
λ
.UmdieZuverlässigkeitzuerhöhen,werden zwei,bzw.dreiBauteileineineMashineeingebaut.ImfolgendenBildistdieR S (t)
füreinBauteil,zweiunddreiBauteileparallelgeshalteteingezeihnet.Mansieht,dass
wenndieZeitwähst,derprozentuelleUntershiedkleinerwird.Dasheiÿt,die
VerbesserungderAusfallratemahtsihvoralleminkurzenZeitspannenbemerkbar,für
grosses
t
gehtR S (t) = e −λ
1t + e −λ
2t − e − ( λ
1+ λ
2)t .
0 10 20 30 40 50 60 70 80
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
R(t)
t
ein Bauteil zwei Bauteile parallel drei Bauteile parallel vier Bauteile parallel
0 10 20 30 40 50 60 70 80
0 50 100 150 200 250 300
Verbesserung in Prozent
t
zwei Bauteile parallel drei Bauteile parallel vier Bauteile parallel
Prof.ErikaHausenblas (Leoben ) Zuverlässigkeit 10.Jänner2017 10/35
Gegeben:
n
Bauteileparallel geshaltetBauteil
j
hateinekonstanteAusfallrateλ j
.Fragen:WielautetdieAusfallratedesGesamtsystems?BerehnensiedieMTFund
dieVerteilungvon
T f
.Sei
p j
dieWahrsheinlihkeitdasBausteinej
funktioniert.DamitgiltP =
1− Q n
i =
1(
1− p j ).
SetztmandieRealibility Funktionein,erhältman
R(t) =
1− Q n
i =
1(
1− R j (t))
undsomit
R(t) =
1− Q n
i=
1(
1− e −t λ j ).
RehnetmansihdieMTFaus,erhältman
MTTF
= P n j=
11
λ j − P n−
1j=
1P n i=j+
11
λ i + λ j
+ P n−
2j=
1P n−
1i =j+
1P n k=i +
11
λ i + λ j + λ k − · · · + (−
1) n−
1P
1λ j .
DadasSystemausfällt,wennalleKomponentenausfallen, gilt
T f =
max{T
1, . . . , T n }
wobeiT i
dieLebensdauerderi
-tenKomponenteist.Fürgroÿe
n
kannmandieVerteilungvonT f
mittelsExtremalverteilungenberehnen.High level redundany
⇔
low level redundanyAngenommen, wirhabeneinSystemzuentwerfen,dasszweiOperationenmittels
Baustein
A
undB
(mitZuverlässigkeitr A
undr B
)inSerieausführenmuss:UmdieZuverlässigkeitdesSystemszuerhöhen,willmanjezweiBauteile
A
undB
verwenden.EskommendabeidiezweimöglihenAnordnungeninBetraht:
WelheobigeAnordnunghatdiehöhereZuverlässigkeit?
EsgiltfürdieAnordnungseriell-parallel
r sp =
1− (
1− r B )
21
− (
1− r A )
2undfürdieAnordnungparallel-seriell
r ps =
1− (
1− r A r B )
2.
BerehnetmandieDierenz
r sp − r ps =
2r A (
1− r A )r B (
1− r B ),
sosiehtman,dassdieseDierenzimmerpositivist. DamitistdieAnordnung
seriell-parallel,besseralsparallel-seriell.
High level redundany
⇔
low level redundanyAngenommen, wirhabeneinSystemzuentwerfen,dasszweiOperationenmittels
Baustein
A
undB
(mitZuverlässigkeitr A
undr B
)inSerieausführenmuss:UmdieZuverlässigkeitdesSystemszuerhöhen,willmanjezweiBauteile
A
undB
verwenden.EskommendabeidiezweimöglihenAnordnungeninBetraht:
WelheobigeAnordnunghatdiehöhereZuverlässigkeit?
SetztmandieWerte
r A =
0.
95undr B =
0.
9ein,bekommtmanfürdieAnordnungseriell-parallelundparallel-seriell
r sp = (
1−
0.
052)(
1−
0.
12) =
0.
987525, r ps = (
1−
0.
05·
0.
1) =
0.
978975.
FürdieWerte
r A = r B =
0.
9,sobekommtmanfürdieAnordnungseriell-parallelundparallel-seriell
r sp = (
1− (
1−
0.
9)
2)
2=
0.
9801, r ps = (
1−
0.
92)
2=
0.
9639.
Gegeben: vier Propeller
mitBuhstaben
A
,B
,C
und
D
.DieAussage dasFlugzeug iegt, wenn
mindestens ein Propel-
ler auf jedem Flügel
funktioniert, führt da-
zu, dass diebeiden Flü-
gelSeriellzubetrahten
- Ausfall der Antriebs-
FunktionaufbeidenFlü-
gelnentsprihtdemSys-
temausfall.
Parallel und Seriellmussniht immermitder geometrishenAnordnung
übereinstimmen. Angenommen manhateinen Stromkreismiteinen
Shalter. Man möhtedie Zuverlässigkeit desShalters(beim unterbrehen
des Stromkreises)erhöhen, indemmaneinen Shalterhinzufügt.Diesen
mussman dannin Serie(physikalish)zushalten,aberin denBerehnung
modelliertmandiese Systemparallel.Umgekehrt,angenommen der
Stromkreisistnormalerweiseunterbrohen undman möhteihnnur im
FalleeinesUnfalls zushalten,mussman die Shalterparallel anbringen,
um die Zuverlässigkeitzu erhöhen.
Ein Systembestehtaus mehrerenSubsystemen,die jeweils vershiedene
Ausfallsraten besitzen. Berehnen Siedie Zuverlässigkeit
R(t)
.Abbildung:ReliabilityBlokDiagramm
figure17
SiehabendieAuswahlzwisheneinenteuren,abersehrZuverlässigenBauteil(B1),oder
parallelgeshaltetbilligenBauteilen(B2).
B1kostet779,dieZuverlässigkeitistexponentialverteiltmitParameter0
.
001;B1kostet249,dieZuverlässigkeitistexponentialverteiltmitParameter0
.
006;NurBauteil1:
R
1(t) =
1− e −
0.
001t
;ZweiBauteile2:
R
2(t) =
1−
1− e −
0.
005t
2=
2e −
0.
005t − e −
20.
005t
;DreiBauteile2:
R
3(t) =
1−
1− e −
0.
005t
3;
VierBauteile2:
R
4(t) =
1−
1− e −
0.
005t
4;
0 10 20 30 40 50
0.95 0.955 0.96 0.965 0.97 0.975 0.98 0.985 0.99 0.995 1
Realibility
Zeit t
kurzes Zeitinterval Bauteil 1 2 Bauteile 2 3 Bauteile 2 4 Bauteile 2
0 50 100 150 200 250 300 350
0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1
Realibility
Zeit t
langes Zeitinterval
Bauteil 1
2 Bauteile 2
3 Bauteile 2
4 Bauteile 2
Abbildung:ReliabilityBlokDiagramm
ebling-105
DasSystem wirdinzwei Systemezerlegt:
Abbildung:ZerlegtesDiagramm
ebling-105-2
Nah derFormelvon Bayesgiltsomit
R ges = R E R zweitesSys + (
1− R E ) R drittesSys . .
DasSystem wirdinzwei Systemezerlegt:
Abbildung:ZerlegtesDiagramm
ebling-105-2
Nah derFormelvon Bayesgiltsomit
R ges = R E R zweitesSys + (
1− R E ) R drittesSys . .
Angenommen
R A = R B =
0.
9,R C = R D =
0.
95 undR E =
0.
80.Berehnen Sie dieZuverlässigkeit desGesamtsystems.
DasBremssystemin einemAutobesteht auseinemhydraulishen System
(Bremspedal) undaus einemmehanishen System(Handbremse). Beide
Untersysteme müssen ausfallen,damitdas Bremssystemversagt.Das
hydraulishe Systemfälltaus sobaldderHauptbremszylinderausfällt
(Ereignis
M
), einerder BlokierkraftregleramRad ausfällt(EreignisseWC
1, . . . , WC
4),oderder Bremsbelag fälltaus (Ereignisse
BP
1, . . . , BP
4 ).Dasmehanishe Bremssystemfälltaus fallsdas Kabelsystemausfällt
(Ereignis
C
), oder beibeidenHinterradbremsenderBremsbelagausfällt.Denition
Ein minimalerWeg istein Weg im Blokdiagramm.FallsalleKomponenten
auf diesenWeg funktionieren,so funktioniertdas Gesamtsystem.
Eine minimaleShnittmenge istein ShnittimBlokdiagramm.Fallen alle
Komponenten diesesShnittes aus,so fälltdasGesamtsystem aus.
MittelssolherMengenkann manuntere und obereShranken der
Gesamtzuverlässigkeit einesSystemsangeben. Hierist eswihtigdassdie
Mengen aus disjunktenElementenbestehen,und dieseunabhängig sind.
Denition
Ein minimalerWeg istein Weg im Blokdiagramm.FallsalleKomponenten
auf diesenWeg funktionieren,so funktioniertdas Gesamtsystem.
Eine minimaleShnittmenge istein ShnittimBlokdiagramm.Fallen alle
Komponenten diesesShnittes aus,so fälltdasGesamtsystem aus.
1
Angenommen die MinimalenShnittmengensind
{S k cut : k =
1, . . . , K }
.Dann giltR ges ≤
K
Y
k=
1n
1
− Y
i∈S k cut
1
− R i o .
2
Umgekehrt, seiendie MinimalenWegegegeben durh
{S k path : k =
1, . . . , K }
.Dann giltR ges ≥
1−
k ′
Y
k=
h
1
− Y
path
R i i
.
1
AngenommendieMinimalenShnittmengensind
{S k cut : k =
1, . . . , K}
.DanngiltR ges ≤
K
Y
k=
1n
1
− Y
i ∈S k cut
1
− R i
o .
2
Umgekehrt,seiendieMinimalenWegegegebendurh
{S k path : k =
1, . . . , K}
.Danngilt
R ges ≥
1−
k ′
Y
k=
1h
1
− Y
i∈S path k
R i
i .
Aufgabe
Angenommen
R D = R A =
0.
95,undR B = R c = R E = R f =
0.
92.Geben Sieeineuntereund obere Abshätzungder ZuverlässigkeitdesGesamtsystemsan.
MinimaleWege:
{A, B} , {A, C } , {D, E } , {D, F }
.MinimaleShnittmengen:
{A, D}
,{B, C, E , F}
.In diesemFallsind die
n
Bauteileparallel angeordnet, abernur einesist in Betrieb. Fälltdieses aus,kommtdas nähste zum Einsatzundwirdangeshaltet. Fälltdieses aus,kommt dasnähste dran. DasSystemfällt
aus, sobaldauh das
n
-teBauteil funktionsuntühtigist.SystememitkalterReserve (Cold standby)
Systememitheiÿer Reserve(Hot standby)
Der Bauteilist während der Reserve ausgeshaltet
⇒
er kann alsonihtkaputtgehen.
Für die Ausfallzeit
T f S
für das Gesamtsystemgilt(Ausfallzeitdes bi
tenBauteils ist
T f i
)T f S = P n i=
1T f i .
In diesemFallgilt
f T S
f (t) = (f T i
f ) (n) ⋆ (t).
Für die MTTFgilt,
MTTF
= − d
f ˆ T S
f
ds = λ n .
EinSystembestehtaus
n
Bauteilen.AngenommendieLebensdauerdern
BauteileistexponentialverteiltmitParameter
λ
.Dannistf T i
f (t) = λe −λ t ,
undf ˆ T i
f (s ) = s+ λ λ , i =
1, . . . n.
DamitgiltfürdieLaplaetransformiertedesGesamtsystems
f ˆ T S f
folgendes
f ˆ T S
f = (s+ λ λ n ) n .
DiesistabergenaudieLaplaetransfomiertederGammaVerteilungundes giltdamit
f T S
f (t) = λ (n− n t n −
11
)! e −λ t .
Da
R (t) =
1− R t
0
f T S
f (s) ds
gilt,ist
R(t)
PoissonverteiltmitParameterλt
.D.h.R(t) = e −λ t P n−
1k=
0( λ t) k
k! .
EinSystembestehtaus2Bauteilen.AngenommendieLebensdauerderbeidenBauteile
istexponentialverteiltmitParameter
λ
1 undλ
2.Dannistf T
1f (t) = λ
1e −λ
1t ,
undf ˆ T
1f (s) = s+ λ
1λ
1
,
und
f T
2f (t) = λ
2e −λ
2t ,
undf ˆ T
2f (s) = s+ λ
2λ
2
,
DamitgiltfürdieLaplaetransformiertedesGesamtsystems
f ˆ T S f
f ˆ T S f = (s+ λ
1λ
1
) λ
2
(s+ λ
2
) .
Somit
f T S
f (t) = λ λ
1λ
21
−λ
2
e −λ
1t − e −λ
2t .
und
R(t) = λ
1e − λ λ
2t −λ
2e − λ
1t
1
−λ
2
und
MTTF
= λ
11
+ λ
12
.
BauteilesindwährendsiewartenimStandbymodus,
⇒
siekönnen ausfallen.BezeihnenwirdieDihteder AusfallwahrsheinlihkeitimStand byModusmit
f ¯ T i
f
,soerhaltenwirfür dieAusfallzeitzeit
T f S
fürdasGesamtsystemT f S =
n
X
i=
1T f i ,
wobeidieAusfallzeitendes
i
tenBauteilsmitT f i
bezeihnetwerden.Eswirdangenommen,dassdieAusfallzeitendereinzelnenBauteileunabhängig
voneinandersind.Indiesem Fallgilt
f T S
f (t ) = (f T i
f ) (n) ⋆ (t ).
EinSystembestehtaus
n
Bauteilenundfunktioniertnur,fallsmindestensk
Bauteilefunktiontühtigsind.
DieWahrsheinlihkeitdasseinBauteilzur Zeit
t
funktioniertistgleihp
Gefragt istdieWahrsheinlihkeit,dasszueinemZeitpunkt
t
genaukBauteilefunktionieren.
P (
genaukBauteilefunktionieren) = n
k
p k (
1− p) n−k .
DamitdasGesamtsystemfunktionstühtigist,müssenmindestens
k
ausn
Bauteilefunktionieren.Esgiltsomit
P (
mindestensk Bauteilefunktionieren) = P k i=
1n i
p i (
1− p) n − i .
Ist
k =
1 erhaltenwir eineParallelshaltung, istk = n
erhalten wireineSerienshaltung.
Commenausefailure
◮
Designormaterialdeieny
◮
Installationerror
◮
Maintenaneerror
◮
harshenvironment(vibration,ontamination,radiation,
. . .
);Casading failure(propagationfailure- Dominoeet)
Negative Abhängigkeit: Angenommen manhat zwei Komponenten.
Sei
A i
das Ereignis,dass Komponentei
ausfällt.◮
DieKomponentenhaben einepositive Anhängigkeit,falls
P(A
1| A
2) > P(A
1)
undP(A
2| A
1) > P(A
2)
gilt.◮
DieKomponentenhaben einenegativeAnhängigkeit,falls
P(A
1| A
2) < P(A
1)
undP(A
2| A
1) < P(A
2)
gilt.Thesquare rootmodel
The
β
fatormodelThebinomialfailurerate model;
The
β
fator modelEs gibt zwei Typen von Fehlern:
FehlerTyp1:Fehler verursaht durhdie Komponente;
FehlerTyp2:Fehler verursaht durhommonlaw;
Es gibt zwei Ausfallraten:
λ (i )
die Ausfallratederi
ten Komponenteaufgrund vonFehlernTyp 1;λ (c)
die AusfallrateaufgrundvonFehlernTyp2.Annahme: EinEintretendesFehlersvon Typ1 undein Eintreten des
FehlersvonTyp2 sindunabhängig.
Commonause fator:
β = λ (i λ ) (c)
+ λ (c) = λ λ (c) .
The
β
fator modelGegeben Bauteile
A
1, · · · A n
.Modellierungeines CommenCausesbei einen in Serieanghängten extraBauteilC
.Annahme:
n
Komponenten mitindividueller Ausfallrateλ (i)
sindgegeben;
EinShok Eregnisverursaht zugleiheinenShaden an mehreren
Komponenten.
DieWartezeitauf das Shokereignisist exponentiellverteilt mit
Parameter
ν
.Trittso ein Shokeregnisein,so gilt