Cours de Fiabilite
Cours de Fiabilite
Cours de Fiabilite
f(t)
(t)=
R(t)
Synthèse :
− (u ).du = ln(1 − F (t )) e 0
t − ( u ).du
= 1 − F (t ) = R (t )
0
R(t ) = e 0
( u ).du
−
F (t ) = 1 − R(t ) = 1 − e 0
( u ).du
−
= (t ).e 0
dF (t ) − ( u ).du
f (t ) =
dt
MTBF = E (T ) = t .f (t ).dt
0
E1 E2 Ei En
n
Rs = P (S ) = P (S1 S 2 ... Si ...Sn ) = P (S1).P (S 2)....P (Si )....P (Sn ) ➔ Rs = Ri
i =1
En
Rs = 1 − (1 − Ri )
i =1
Dans un système //, la fiabilité du système est plus grande que la plus grande des fiabilités des éléments composant le
système. On utilise ce fait pour améliorer la fiabilité ; cela réalise une redondance active.
Si on désire effectuer un calcul en fonction du temps, on doit introduire la fonction R(t).
n
Si R(t ) = e
− t
, alors Rs = 1 − (1 − e− t ) .
i =1
E1 Hypothèse : le taux de défaillance des éléments E1 et E2 est constant et est égal à e1 et e 2 .
Cette hypothèse a pour conséquence que les lois de fiabilité sont de type exponentiel :
Re1(t ) = e−e1t et Re 2 (t ) = e−e 2t
DC
E2 On fait aussi l’hypothèse que la fiabilité de l’organe DC est égale à 1.
Il sera facile par la suite de la prendre en compte par la suite dans le calcul, cet organe étant en série avec le système
{E1, E2}.
Re1(t ) = e − e1t et Re 2 (t ) = e − e 2t
fe1(t ) = e1e − e1t et fe 2 (t ) = e 2e − e 2t
Le système fonctionnera avec E1 ou E2, ces événements étant mutuellement exclusifs (E1 sans E2 ou E2 sans E1,
mais jamais les 2 en même temps).
Rs ( t ) = 1xe − e1t + ( 0
t
e1e − T .dT xe −
e1
) e 2 ( t −T )
0
t
= e − e1t + e1. e − e1T .dT xe − e 2 .t xe e 2 .T
t t
Rs ( t ) = e − e1t + e1.e − e 2 .t . e − e1T .e e 2 .T .dT = e − e1t + e1.e − e 2 .t . e −( e1 −e 2 )T .dT
0 0
t
− e 2 .t e
− ( e 1 − e 2 )T
t
Rs ( t ) = e − e1t + e1.e − e 2 .t . e −( e1 −e 2 )T .dT = e − e1t + e1.e .
−(e1 − e 2 ) 0
0
− e 1t e −( e1 −e 2 )t
− e 2 .t 1
Rs ( t ) = e + e1.e . −
−(e1 − e 2 ) −(e1 − e 2 )
− e 2 .t 1 − e
− ( e 1 − e 2 )t
Rs ( t ) = e − e1t + e1.e .
(e1 − e 2 )
e1.e − t − e 2 .e − t + e1.e − .t − e1.e − .t .e −(
e1 e1 e2 e2 e 1 − e 2 )t
Rs ( t ) =
e1 − e 2
e1.e − t − e 2 .e − t + e1.e − .t − e1.e −
e1 e1 e2 e 2 .t − e 1 .t + e 2 .t
Rs ( t ) =
e1 − e 2
e1.e − t − e 2 .e − t + e1.e −
e1 e1 e 2 .t
− e1.e − e1.t e1.e − .t − e 2 .e −
e2 e 1t
Rs ( t ) = =
e1 − e 2 e1 − e 2
Si on prend en compte l’élément de détection et de commutation DC, on obtient alors :
Rs ( t ) = e .
e1 − e 2
Remarque : si on considère que tous les éléments ont le même taux de défaillance λ, on obtient alors l’expression
suivante : Rs ( t ) = e− DC .t .e − .t .(1 + .t )
T1
M2 M3 M4 M5 T2 T3
M1 T1
0,99 0,99 0,99 0,99 0,99 0,99
M1 T1
R(t ) = e 0 R(t ) = e − .t
− ( u ).du
et comme (u ) = cte =
t ➔
R(t ) = e 0
.du
− t
= e − .u = e − .t
0
− .t
• Densité de probabilité : f (t ) = .e
− .t
• Fonction de répartition : F (t ) = 1 − R(t ) = 1 − e
1
• Espérance mathématique : E (t ) = MTBF =
5.2 – Durée de vie associée à un seuil de fiabilité :
Il est intéressant de savoir à quel instant la fiabilité atteindra un seuil déterminé.
1 1 1
R(t ) = e − .t ln R(t ) = −.t t = − .ln R(t ) t = .ln
R (t )
Ex : un composant a une MTBF de 2000 heures. A quelle date « tj » ce composant aura une fiabilité de 90% ?
1 1 1 1
tj = .ln = MTBF .ln = 2000 x ln = 211 heures
R (t ) R (t ) 0,9
Au bout de 211 heures, on estime donc que 90% des composants survivront.
5.3 – Représentation graphique de la loi exponentielle :
− .t -.t
Si R(t ) = e , alors ln R (t ) = −.t en logarithmes népériens et logR(t)= en logarithmes décimaux.
2,3
Loi exponentielle sur échelle linéaire Loi exponentielle sur papier semi logarithmique
R(t) logR(t)
Dr
oit
ed
ep
en
0,1 te
/2
,3
1/e=0,368
m=1/ t ou t
t
2,3/
5.4 – Estimation du taux de défaillance :
• Porter sur papier semi logarithmique les N points formés des couples (ti, Ri)
• Tracer la courbe de régression des N points
Dr. M. GUEMANA, FT, Université de Médéa 8
• Si les N points sont sensiblement alignés, alors la loi de fiabilité est exponentielle
• Déterminer λ par la pente de la courbe
• En déduire MTBF = 1/ λ
• En déduire R(t ) = e − .t
VI – ANALYSE DE LA FIABILITE PAR LA LOI DE WEIBULL :
61 – Définition de la loi de Weibüll :
C’est une loi de fiabilité à 3 paramètres qui permet de prendre en compte les périodes où le taux de défaillance n’est
pas constant (jeunesse et vieillesse). Cette loi permet :
• Une estimation de la MTBF
• Les calculs de λ(t) et de R(t) et leurs représentations graphiques
• Grâce au paramètre de forme β d’orienter un diagnostic, car β peut être caractéristique de certains modes de
défaillance
Les 3 paramètres de la loi sont :
β ➔ Paramètre de forme >0 sans dimension:
• Si β>1, le taux de défaillance est croissant, caractéristique de la zone de vieillesse
o 1,5 < < 2,5 : fatigue
o 3 < < 4 : usure, corrosion
• Si β=1, le taux de défaillance est constant, caractéristique de la zone de maturité
• Si β<1, le taux de défaillance est décroissant, caractéristique de la zone de jeunesse
f(t) R(t) (t)
=3
1 =3 1
=1
0,5
=0,5 =0,5
=1
t =3 =1 t t
γ ➔ paramètre de position, - < γ < +, qui s’exprime dans l’unité de temps :
• γ>0 : survie totale sur l’intervalle de temps f(t)
[0, γ]
• γ=0 : les défaillances débutent à l’origine
des temps
• γ<0 : les défaillances ont débuté avant
l’origine des temps ; ce qui montre que la
mise en service de l’équipement étudié a <0 =0 >0
précédé la mise en historique des TBF
t
Relations fondamentales :
−1 t −
t − −
• f (t ) = avec t
Densité de probabilité : . .e
t −
−
• Fonction de répartition : F (t ) = 1 − e
1
E (t ) = MTBF = + .(1 + ) = A +
= B
Γ est une fonction mathématique complexe. A et B sont des paramètres issus de tables et sont calculés à partir de la loi
Γ.
AXE a
AXE A
AXE b
AXE B
AXE A
• Axe A : axe des temps sur lequel on porte les valeurs ti des TBF
• Axe B : valeurs des probabilités de défaillance Fi calculées par la méthode des rangs moyens ou des rangs
médians. On estime R(t) par R(t) = 1 – F(t)
• Axe a : axe des temps en logarithmes népériens : ln(t)
• Axe b : axe qui permet l’évaluation de β
6.4 – Détermination graphique des paramètres de la loi :
1. Préparation des données : détermination des couples (ti, Fi) par les rangs moyens ou les rangs médians
2. Tracé du nuage de points
3. Tracé de la droite de Weibüll
4. Détermination de β, η, γ
5. Détermination des équations de la loi de Weibüll
6. Calcul de la MTBF
7. Exploitation des données issues de la loi
Exemple d’application :
Préparation des données :
Ordre i TBF Fi
η=770 heures
β=1,4
D2 D1