Fiabilité Cours
Fiabilité Cours
1. Introduction
2. Rappels mathématiques
3. Mesure de la fiabilité
4. Maintenabilité
5. Disponibilité
6. Sécurité
03/03/2018 2
1 - Introduction
03/03/2018 4
Exemple : télévision
Ma télé doit
fonctionner dès
que je la branche
Mais si elle est en
panne, elle doit être
réparée tout de suite
Et surtout, le tube
ne doit pas
imploser!
Disponibilité
Fiabilité
03/03/2018 + Maintenabilité + Sécurité = 5
24h/24
Définitions (norme NF EN 13306)
03/03/2018 8
Définitions (norme NF EN 13306)
03/03/2018 10
2 – Rappels mathématiques
P
Tmoy
t (heures)
500 1000 1500
1
0,8
t (heures)
500 1000 1500
Ex : proba pour que 80% des équipements soient en panne avant 1500h
2 – Rappels mathématiques
4. Densité de probabilité
C’est la dérivée de la fonction de répartition : c’est
encore la probabilité pour qu’un événement arrive
dans un intervalle de temps donné
Je teste sa 1. La fiabilité
fiabilité, 2. Modélisation mathématique
Ah, Ah,
Ah.. 3. Fiabilité des équipements
complexes
4. Paramètres fiabilistes
5. Cas typiques de calcul de
fiabilité
6. Amélioration de la fiabilité
03/03/2018 15
3.1 – La fiabilité
03/03/2018 17
B – Le temps
La base :
temps jusqu’à défaillance pour un
système non réparable
ou temps entre 2 défaillances
consécutives pour un système réparable
s’exprime en secondes ou bien en « unités
d’usage » (heures, distances, cycles ou
toutes grandeurs appropriées : 1 heure, 3
mois, 1000 km, etc...)
définit une variable aléatoire T que l’on
traite selon les méthodes usuelles du
calcul des probabilités
03/03/2018 18
Dispositifs non réparables
TTF
Arrêt
t
TBF
Fonctionnement
UT
Arrêt
TTR TTR
DT DT
03/03/2018 20
TBF = Time Between Failure (temps s’écoulant entre deux
défaillances)
03/03/2018 21
C – Fiabilité = probabilité
90 chances sur
La fiabilité est une
100 pour qu’elle
caractéristique d’un
me raconte des
équipement qui s’exprime
salades !..
sous forme probabiliste
Son estimation:
- à partir de résultats
obtenus sur une période
donnée
- à partir de résultats
obtenus sur un échantillon
d’équipements
03/03/2018 22
C – Fiabilité = probabilité
1. Expression mathématique
2. Estimation statistique
3. Application
x 4. Taux de défaillance
5. Modèle mathématique de la
fiabilité
∫sin3udu
03/03/2018 24
1 – Expression mathématique
03/03/2018 27
2 – Estimation statistique
28
2 – Estimation statistique
N(i)-N(i+1) N(i)-N(i+1)
i < t < i+1 f(i) = λ(i) =
N(0) N(i)
F(i+1)=1–R(i+1)
N(i+1)
R(i+1) = i
t = i+1
03/03/2018
N(i+1) N(0) =∑f(k)
k =0
31
3 - Application
Période Nb de
L’observation des défaillances défaillance
s
sur un ensemble de 40
0 – 1000h 1
moteurs asynchrones de
1000 – 2000h 4
même type, utilisés par la 2000 – 3000h 8
société MECAMALO, a donné 3000 – 4000h 12
les résultats suivants 4000 – 5000h 11
5000 – 6000h 3
03/03/2018 32
Cumul
Date ou Nb Défail- Nb survi- défaillants
période lants vants n R(t) F(t) f(t) λ(t)
0 0 40 0 1 0
[0,1000[
1000
[1000,2000[
2000
[2000, 3000[
3000
[3000, 4000[
4000
[4000, 5000[
5000
[5000, 6000[
6000
Date ou Cumul défail
période Nb Défaillants Nb survivants lants n R(t) F(t) f(t) λ(t)
0 40 0 100,0% 0,0%
[0,1000[ 1 2,5% 2,50%
1000 39 1 97,5% 2,5%
[1000,2000[ 4 10,0% 10,26%
2000 35 5 87,5% 12,5%
[2000, 3000[ 8 20,0% 22,86%
3000 27 13 67,5% 32,5%
[3000, 4000[ 12 30,0% 44,44%
4000 15 25 37,5% 62,5%
[4000, 5000[ 11 27,5% 73,33%
5000 4 36 10,0% 90,0%
[5000, 6000[ 3 7,5% 75,00%
6000 1 39 2,5% 97,5%
03/03/2018 34
Allure de la fiabilité R(t)
R(t)
%
100
50
10
t
03/03/2018 1000 2000 3000 4000 5000 6000 35
Allure de la fonction de défaillance F(t)
F(t)
%
100
50
10 3500 h
t
03/03/2018 1000 2000 3000 4000 5000 6000 36
Allure de la fréquence d’apparition des
défaillances f(t)
p(t)
%
E(t)
30
20
10
t
1000 2000 3000 4000 5000 6000
03/03/2018 37
Allure du taux de défaillance λ(t)
λ (t)
%
100
50
λ(t) Max
10
t
03/03/2018 1000 2000 3000 4000 5000 6000 38
Remarques sur le calcul de F(t)
quand l’effectif de
l’échantillon est inférieur à F(t) = n − 0,3
N(0) + 0,4
20, on utilise la méthode
des rangs médians
03/03/2018 40
4 – Taux de défaillance
03/03/2018 41
A – Taux de défaillance moyen
03/03/2018 42
A1 - Caractérisation d’un matériel unique
03/03/2018 43
Exemple
une ligne de production a fonctionné pendant un cycle de
200 heures
03/03/2018 45
Exemple
10 dispositifs identiques ont été testés sur une durée de 1000 h.
6
λ= = 7,71.10 − 4 def/h
260 + 425 + 624 + 712 + 803 + 955 + 4 x1000
03/03/2018 46
Exemples de taux moyen de défaillances
Composants électroniques
03/03/2018 47
Exemples de taux moyen de défaillances
Composants mécaniques
λ = 35 - 30 = 0,014 def/h
35x10
03/03/2018 50
B – Taux de défaillance instantané
Composant
mécanique
Composant
électronique
Composant
électrique
t
03/03/2018 53
C - Evolution du taux de défaillance dans le temps
(« déverminage »)mécanique
pour des composants électroniques
Composant
électronique
Composant
électrique
t
03/03/2018 54
C - Evolution du taux de défaillance dans le temps
Zone 2
c’est la période, dans la vie d’un matériel, pendant laquelle le
tauxλ instantané de défaillance est pratiquement constant pour le
1 2 3
composant électronique, moins pour le composant mécanique
elle peut être de durée plus ou moins importante (plus pour
l’électronique que pour la mécanique)
les défaillances sont aléatoires et liées le plus souvent à la
Composant
dérive des composants
mécanique
Composant
on appelle aussi cette période « période de maturité »
électronique
Composant
en maintenance, c’est la période où l’on metélectrique
en place du
t
préventif,
03/03/2018
même si le correctif reste nécessaire 55
C - Evolution du taux de défaillance dans le temps
Zone 3
c’est la période de défaillance par vieillissement (ou période
λ
d’usure ou1 « fin de vie ») 2 3
pendant cette période, le taux de défaillance croît rapidement
à un certain point de λ(t), le matériel est « mort »
on peut éliminer cette période par des politiques appropriées de
déclassement ou de remplacement systématique
Composant
mécanique
la tendance est toutefois d’effectuer une maintenance
Composant
conditionnelle, ce qui permet de prévoir les défaillances et
électronique
Composant
d’exploiter le matériel au maximum de ses possibilités
électrique
t
03/03/2018 56
C - Evolution du taux de défaillance dans le temps
λ
Sur cet exemple, les pannes de
jeunesse ont été éliminées
Les défaillances
n’apparaissent que lorsque
l’équipement est âgé
t
Exemple : équipements
mécaniques et
Modèle B électromécaniques simples
03/03/2018 58
C - Evolution du taux de défaillance dans le temps
λ
Cet exemple indique qu’on
ne sait pas dire à quel
moment vont commencer
les ennuis, puis on a un
t taux croissant lié à l’usure
Modèle C
03/03/2018 59
C - Evolution du taux de défaillance dans le temps
λ
Cette courbe montre des
systèmes sans défaut au
départ
t Cas des moteurs
électriques achetés neufs
Modèle D
03/03/2018 60
C - Evolution du taux de défaillance dans le temps
λ
Cette courbe montre des
systèmes souffrant de
défaillances précoces
t
C’est le signe, le plus
souvent, d’une mauvaise
conception
Modèle F
03/03/2018 62
C - Evolution du taux de défaillance dans le temps
03/03/2018 63
Exemple 1 : roulement à billes
03/03/2018
t
∫0 λ ( x ).dx =
t dF(x)
∫0 1 − F(x) = [− ln(1 − F( x ))]t
0 66
5 - Modèle mathématique de la fiabilité
Soit encore :
∫0 λ( x ).dx = − ln[R( x )]
t t
0 = − ln R( t ) + ln R(0)
t
Mais comme R(0) = 1, il vient : ln R( t ) = − ∫ λ( x ).dx
0
t
− ∫ λ( x ).dx
ou encore : R( t ) = e 0
03/03/2018 68
3.4 - Paramètres fiabilistes de type
« moyenne »
On va obtenir, avec M = Mean (moyenne) :
le MTBF (Mean Time Between Failures)
le MTTR (Mean Time To Restore)
le MUT (Mean of Up Times)
le MDT (Mean of Down Times)
le MTTF (Mean Time To Failure)
Parmi ces indicateurs, le MTBF a un très grand intérêt
03/03/2018 69
MTBF
C’est l’espérance mathématique du temps de bon
fonctionnement entre deux défaillances
Son expression est donc :
+∞ +∞
MTBF = ∫ t.f ( t ).dt = ∫ t.dF( t )
−∞ −∞
+∞
MTBF=∫ R(t).dt
0
03/03/2018 70
Remarques sur le MTBF
MTBF ≅
∑ UT
Nombre de défaillanc es
03/03/2018 71
Exemple
R(t) = e ∫0 = e ∫0
− λ(t).dt − λ.dt
t
Ici λ est constant donc : ∫0 λ.dt = λ.t
Dans ces conditions : R( t ) = e − λt
+∞
On sait que MTBF=
∫0 R(t).dt
donc MTBF=∫
+∞
0
[
e−λt.dt = − 1.e−λt
λ
]+∞
R1 R2 Rn
Calcul de MTBF
1
MTBF =
λ1 + λ2 + ...λn
Calcul de la Fiabilité d’un Système
(λ=constante)
Systèmes constitués de plusieurs composants en parallèle:
Le système est défaillant, si et seulement
R1
si tous les éléments du système sont
simultanément défaillants ,
R2
La probabilité de panne de l' ensemble :
F (t ) = F1 (t ).F2 (t ).F3 (t )....Fn (t )
La Fiabilité du Système :
n n
R (t ) = 1 − ∏ [1 − Ri (t )] = 1 − ∏ (1 − e
− λi t
Rn )
i =1 i =1
79
2 - Cas d’un taux de défaillance non
constant
Si le taux de défaillance est de la forme
λ(t) = at + b
at 2
alors t −( + bt )
− ∫ ( at +b )dt
R( t ) = e =e
2
0
03/03/2018 80
2 - Cas d’un taux de défaillance non
constant
Si l’on connaît les paramètres a et b, on n’aura aucun
mal à calculer R
Par contre, le calcul du MTBF ne sera pas simple
at 2
+∞ −( +bt )
MTBF = ∫ e 2
.dt
0
03/03/2018 83
Modèle de Weibull
β
t−γ
− t>γ
η
R( t ) = e
Pour γ = 0 et β = 1, on retrouve la loi exponentielle
Celle-ci est donc un cas particulier de la loi de
Weibull
03/03/2018 85
Paramètre β
C’est un nombre sans dimension qui détermine la
forme de la distribution f(t) des défaillances :
si β = 1, le taux de défaillance est constant
si 0 < β < 1, le taux de défaillance est décroissant
(période de jeunesse)
si β > 1, le taux de défaillance est croissant
03/03/2018 88
2 – Expressions mathématiques
A - Fonction de défaillance
Par définition, F(t) = 1 – R(t)
β
t−γ
−
η
si t > γ
F( t ) = 1 − e
03/03/2018 89
2 – Expressions mathématiques
B – Densité de probabilité
C’est la dérivée de F(t)
C’est aussi la fréquence d’apparition
d’une défaillance entre 2 instants
β
β−1 − t − γ
β t− γ η
si t > γ
f ( t ) = . e
η η
03/03/2018 90
B - Densité de probabilité
f(t)
β = 3
β = 0 ,5
β = 1
t
Quand β > 1, le fréquence de défaillance passe
03/03/2018 par un maximum 91
Application à la maintenance
A - Taux de défaillance
f(t)
Par définition, λ(t) =
1 - F(t)
β−1
β t−γ
λ( t ) = . si t > γ
η η
03/03/2018 92
A – Taux de défaillance
β−1
λ(t) β t−γ
λ( t ) = .
β=2 η η
β=1
β = 0,5
On montre que :
2 − .Γ1 + 1
2
V(T) = σ2(T) = η2.Γ1 +
β β
03/03/2018 95
β A B β A B β A B
03/03/2018 98
4 - Détermination graphique des
paramètres de la loi de Weibull
1
On pose Y = lnln[ ] X = ln( t − γ )
03/03/2018
1 − F( t ) 100
A – Justification théorique
6) Détermination du MTBF
P1
On les choisit aussi de
manière que les
projections de P1P2 et
P2P3 sur l’axe b soient t 22 − t1 t 3
égales γ=
03/03/2018 2t 2 − t1 − t 3 106
E – Redressement de la concavité
A
03/03/2018 107
3.6 – Amélioration de la fiabilité
1. Caractéristique de la maintenabilité
2. Taux de réparation
3. Temps moyen de réparation après
défaillance
4. Temps moyen d’indisponibilité après
défaillance
5. Amélioration de la maintenabilité
03/03/2018 110
4.1 - Caractéristique de la maintenabilité
Ca risque
d’être
Pour un matériel travaillant dans
long!
des conditions données
d’utilisation, c’est la probabilité
pour qu’une opération de
maintenance active, avec
l’utilisation de procédures et
moyens prescrits, puisse être
effectuée dans un intervalle de
temps donné
03/03/2018 111
4.1 - Caractéristique de la maintenabilité
−∞ 0
03/03/2018 116
Exemple de calcul de taux de
maintenabilité
Une pompe industrielle a fonctionné 10000 heures
pendant lesquelles on a relevé 7 pannes de durées:
4h - 2,5h - 6h - 12h - 1,5h - 36h et 3,5h
Son MTTR est donc
µ= 1 = 1 = 0,11 pannes/h
MTTR 9,357
03/03/2018 117
4.4 - Temps moyen d’indisponibilité
après défaillance
Temps de non
Temps de réparation (TTR)
détection de la
défaillance
Temps d'appel
Temps annexes Temps de remise
maintenance
(administratif, techniques, ..) en condition
03/03/2018 120
4.5 - Amélioration de la maintenabilité
Temps Améliorations possibles
Vérification, • voyants, capteurs intégrés à la conception
localisation • appareils de mesure avec ou sans
enregistrement
• supervision
Diagnostic • documentation opérationnelle complète
(plans mis à jour, notices d’entretien, etc..)
• repérage et accessibilité des points de
mesure
• facilité de diagnostic : diagramme causes-
effets, arbre de diagnostic, logigramme de
03/03/2018 dépannage, système expert 121
4.5 - Amélioration de la maintenabilité
1. Définition
2. Les temps en maintenance
3. Disponibilités théoriques
4. Disponibilité opérationnelle
5. Autres disponibilités
6. Amélioration de la disponibilité opérationnelle
03/03/2018 124
5.1 - Définition
Temps requis
Temps effectif de disponibilité Temps effectif d’indisponibilité
03/03/2018 131
5.2 – Les temps en maintenance
t4
03/03/2018
t5 134
5.3 – Disponibilités théoriques
1. Disponibilité instantanée
1. Disponibilité instantanée
Selon les lois suivis par les MTBF et MTTR,
D(t) s’exprime par des fonctions plus ou
moins complexes de λ(t) et de µ(t)
Si λ et µ constants, on obtient :
µ λ
D( t ) = + .e −( λ +µ )t
µ+λ λ+µ
2. Disponibilité asymptotique
Encore appelée disponibilité
intrinsèque
C’est la limite, lorsqu’elle existe,
de la disponibilité instantanée,
lorsqu’on fait tendre le temps
vers l’infini
03/03/2018 137
5.3 – Disponibilités théoriques
2. Disponibilité asymptotique
Si λ et µ sont constants alors :
µ
Da = limt →∞ D( t ) =
µ+λ
Or λ = 1/MTBF et µ = 1/MTTR donc :
MTBF
Da =
MTBF + MTTR
03/03/2018 138
5.3 – Disponibilités théoriques
2. Disponibilité asymptotique
Par rapport aux temps exprimés précédemment, on
peut encore exprimer la disponibilité asymptotique
sous la forme :
t1
Da =
t1 + t 4
La disponibilité asymptotique est donc celle
fournie par le constructeur d’un équipement.
Celui-ci ne peut pas en effet prendre en compte
les causes extérieures d’indisponibilité et les
03/03/2018 carences du service maintenance 139
5.4 – Disponibilité opérationnelle
MUT MUT
Do = =
MUT + MDT MTBF
03/03/2018 140
5.4 – Disponibilité opérationnelle
2. Disponibilité globale
Elle caractérise le taux global
d’utilisation de l’équipement
temps de disponibilité
DG =
temps total
Autre formule
MTBF
DG =
MTBF + MDT
03/03/2018 143