Chapitre 3

LA FIABILITE DES SYSTEMES DE PRODUCTION

I – APPROCHE DE LA FIABILITE PAR LES PROBABILITES :
Définition selon la NF X 06–501 : la fiabilité est la caractéristique d’un dispositif exprimée par la
probabilité que ce dispositif accomplisse une fonction requise dans des conditions d’utilisation
données et pour une période de temps déterminée.
1. Probabilité : c’est le rapport :
Nb cas favorables
1
Nb cas possibles
On notera R(t) la probabilité de fonctionnement à l’instant t. Le symbole R provient de l’anglais Reliability.
On notera F(t) la fonction définie par F(t)=1-R(t). C’est la probabilité complémentaire. F(t) est la probabilité de
défaillance à l’instant t. F(t)+R(t)=1.
2. Fonction requise : ou accomplir une mission ou rendre le service attendu. La définition de la fonction requise
implique un seuil d’admissibilité en deçà duquel la fonction n’est plus remplie.
3. Conditions d’utilisation : définition des conditions d’usage, c’est à dire l’environnement et ses variations, les
contraintes mécaniques, chimiques, physiques, etc. Il est évident que le même matériel placé dans 2 contextes
de fonctionnement différents n’aura pas la même fiabilité.
4. Période de temps : définition de la durée de mission T en unités d’usage. Ex : on se fixe un minimum R(Tm) =
0,9 pour une durée de mission Tm = 8000 heures ; à tout instant Ti de la mission est associée une fiabilité R(ti).
Ex : moteur de voiture préparé pour les 24 heures du Mans :
 Probabilité : c’est celle de terminer ; fiabilité requise=0,98
 Fonction requise : 200 km/h de moyenne (seuil minimal)
 Conditions d’utilisation : de jour, de nuit, avec de la pluie, n ravitaillements, etc.
 Période de temps : au bout de 24 heures (durée de la mission)
II – EXPRESSIONS MATHEMATIQUES :
21 – Fonctions de distribution et de répartition :
Notion de variable aléatoire : on appelle variable aléatoire X une variable telle qu’à chaque valeur x de la VA X on
puisse associer une probabilité F(x). Une variable aléatoire est donc une fonction qui à chaque évènement d’une
expérience aléatoire associe un nombre réel.
Une VA peut être :
 Continue : intervalle de temps entre 2 défaillances consécutives
 Discrète : nombre de défaillance sur un intervalle de temps
Soit une loi de probabilité relative à une VA continue T.
Cette loi est caractérisée par sa fonction de distribution (appelée aussi densité de probabilité) f(t) et par sa fonction de
répartition F(t) telles que :
dF (t ) P (t  T  t + dt )
f (t ) = = lim
dt dt �0 dt
La fonction F(t) représente la probabilité qu’un évènement (défaillance) survienne à l’instant T dans l’intervalle [0,t].
F (t ) = P (T  t )
ti
Comme f (t ).dt = P (t  T  t + dt ) � F (t ) = �f (t )dt = P (T  ti )
-�
n
Remarque : si la VA est discrète, l’expression devient : F (tn ) = �f (ti ) =P (T  tn )
0

Chapitre 3 LA FIABILITE Page 44

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
22 – Application à la fiabilité :
Un dispositif mis en marche la 1ère fois à t=0 tombera inexorablement en panne à un instant T non connu à priori.
T (date de la panne), est une VA de la fonction de répartition F(t).
 F(t)  probabilité de défaillance avant un instant ti
 R(t)  probabilité de bon fonctionnement à ti
 R(t) + (F(t) = 1
t �
 �f (t )dt + �f (t )dt = 1
0 t

23 – Taux de défaillance :
On définit le taux de défaillance de la manière suivante :
nombre de défaillants sur un intervalle de temps
l (t)=
nombre de survivants au début de la période x intervalle de temps
On définit :
 N0 le nombre initial de dispositifs
 Ns(t) est le nombre de dispositifs survivants à l’instant t
 Ns(t + Δt) est le nombre de dispositifs survivants à l’instant t + Δt
Au niveau d’une défaillance, 2 cas peuvent se produire :
 Les défaillants sont remplacés
 Les défaillants ne sont pas remplacés
Les défaillants sont remplacés : Ns(t) sera toujours égal à N0 :
On nomme C(Δt) le nombre de défaillants durant Δt.
C( Dt )
D’après la formule générale du taux de défaillance, on a : l (t)= .
N 0.Dt
Ns(t ) - Ns(t + Dt )
Les défaillants ne sont pas remplacés : l (t)=
Ns (t ).Dt
Ce taux de défaillance est une valeur moyenne sur une période Δt connue. Or, au même titre que F(t) et R(t), il est
intéressant de connaître l’évolution de λ(t) au cours du temps.
C’est le taux de défaillance instantané :
On fait tendre Δt  dt et (Ns(t) – Ns(t + Δt))  dN. dN sera précédé du signe « - » car il y a moins de survivants à (t +
Δt) qu’à t.
-dN -dN
l (t)=  l (t).dt=
N ( t ).dt N (t )
l (t).dt est appelé probabilité conditionnelle de défaillance sur [t, t+dt].
Applications :
Cas N°1 : les défectueux sont remplacés. Une étude a été
menée sur 70 véhicules pendant une période allant de 80000km C( Dt ) 41
l (t ) = = = 0,585.10-4 panne / km
à 90000km. 41 défaillances ont été réparées. Déterminer le taux No.Dt 70.(90000 - 80000)
de défaillance pour cette période.

Chapitre 3 LA FIABILITE Page 45

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
Cas N°2 : les défectueux ne sont pas remplacés. On teste un lot de 50 électrovannes soumises en continu à 8
impulsions par minute. A la 50ème heure, il en reste 33. A la 60ème heure, il en reste 27. Déterminer le taux de défaillance
sur cette classe, par heure et par impulsion.
Ns(t ) - Ns(t + Dt ) 33 - 27
l (t ) = = = 18.10 -3 def / heure = 3,79.10 -5 def / imp.
Ns(t ).Dt 33.10
Si les électrovannes étaient remplacées, on obtiendrait :
C( Dt ) 33 - 27
l (t ) = = = 12.10 -3 def / heure
No.Dt 50 x10

Liaison entre le taux de défaillance et la fiabilité :

« Probabilité d’avoir une panne entre t et dt » = « probabilité de survivre à l’instant t » x « probabilité conditionnelle de
défaillance entre t et t+dt ».
Cette expression est identique à : f (t ).dt = R(t ).l (t).dt � f(t)=R(t).l (t)
Il vient donc l’expression du taux de défaillance en fonction de la loi de fiabilité et la densité de probabilité :
f(t)
l (t)=
R(t)
Synthèse :

Chapitre 3 LA FIABILITE Page 46

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
III – EXPRESSIONS DES LOIS DE FIABILITE :
dF (u )
f (u ) =
du
f (u ) dF (u ) dF (u ) dF (u )
l (u ) = = = � l (u ).du =
R(u ) R(u ).d (u ) (1 - F (u )).du 1 - F (u )
Intégrons les 2 membres entre 0 et t :

t t dF (u ) t t -dF (u )
�l(u ).du = �1 - F (u ) � - �l(u ).du = �1 - F (u )
0 0 0 0

l (u ).du = [ ln(1 - Fu )] 0 = [ ln(1 - F (t )) - ln(1 - F (0))]

t t
-�
0

A t=0, il n’y a pas de défaillance, donc F(0) = 0, donc ln(1-F(0)) = ln1 = 0

t

l (u ).du = ln(1 - F (t )) � e �
t l ( u ).du
-
-� 0
= 1 - F (t ) = R(t )
0

On obtient donc les expressions générales des lois de fiabilité :

t

R (t ) = e �
l ( u ).du
-
0

F (t ) = 1 - R (t ) = 1 - e �
l ( u ).du
-
0

t
dF (t )
= l (t ).e �
- l ( u ).du
f (t ) = 0

dt
�
MTBF = E (T ) = �t .f (t ).dt
0

La MTBF est définie comme étant l’espérance mathématique de la VA T.

IV – LOIS DE COMPOSITION EN FIABILITE : ASSOCIATIONS DE MATERIELS :
Le problème qui se pose à la maintenance au niveau de la fiabilité est son amélioration constante. Il peut pour cela
intervenir sur la technologie du composant, agencer les composants ou sous-systèmes de manière à les rendre plus
fiables par l’utilisation de redondances dont on distingue 3 grandes catégories :
 Les redondances actives
 Les redondances passives ou « stand-by »
 Les redondances majoritaires
41 – Redondance active :
Une redondance active est réalisée par la mise en parallèle d’éléments assurant les mêmes fonctions et
travaillant en même temps.
On a donc à faire à un système appelé par les fiabilistes « système parallèle ».
Hypothèses de départ :
 Les défaillances sont indépendantes les unes des autres
 La fiabilité de chaque sous-système ou de chaque élément a été déterminée

Chapitre 3 LA FIABILITE Page 47

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
Système série :
On dit qu’un système est un système série d’un point de vue fiabilité si le système tombe en panne lorsqu’un seul de
ses éléments est en panne.

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

Cette association est caractéristique des équipements en ligne de production.

Système // :
On dit qu’un système est un système // d’un point de vue fiabilité si, lorsqu’un ou plusieurs de ses éléments tombent en
panne, le système ne tombe pas en panne.
Pour calculer la fonction fiabilité d’un système // à n éléments, ils est plus aisé de passer par la
E1 fonction défaillance F.
F = 1 - R = 1 - P (S ) = P (S )
E2 F = P (S1).P (S 2)....P (Si )....P (Sn ) = F1.F 2....Fi....Fn
F = (1 - R1).(1 - R 2)....(1 - Ri )....(1 - Rn)
Ei Rs = 1 - (1 - R1).(1 - R 2)....(1 - Ri )....(1 - Rn )
n
Rs = 1 - �(1 - Ri )
En 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
- lt
, alors Rs = 1 - �(1 - e - lt ) .
i =1

42 – Redondance passive :
Dans ce cas, un seul élément fonctionne, les autres sont en attente. Ceci a l’avantage de
diminuer ou de supprimer le vieillissement des éléments ne travaillant pas. En contrepartie,
on a l’inconvénient d’être obligé d’avoir un organe de détection des pannes et de
commutation d’un système sur un autre.
Le calcul d’un système à redondance passive ou « stand-by » se fait en tenant compte de la
variable temps. Il faut donc connaître au préalable, pour chaque composant, son taux de
défaillance λ(t) et sa loi de fiabilité R(t).

Calcul d’un système à redondance passive à 2 éléments en // :

E1 Hypothèse : le taux de défaillance des éléments E1 et E2 est constant et est égal à le1 et le 2 .
Cette hypothèse a pour conséquence que les lois de fiabilité sont de type exponentiel :
Re1(t ) = e - le1t et Re 2 (t ) = e - le 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 - le1t et R e 2 (t ) = e - le 2t
fe1(t ) = le1e - le1t et fe 2 (t ) = le 2e - le 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).

R(S) = [Prob(S marche sachant que E1 marche) x Prob(E1 marche)]

Chapitre 3 LA FIABILITE Page 48

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
+ [Prob(S marche sachant que E1 ne marche pas) x Prob(E1 ne marche pas)]
 Prob(E1 ne marche pas)  probabilité que E1 soit défaillant
 Prob(S marche sachant que E1 marche)  = 1 (tant que E1 marche, S fonctionnera toujours)
 Prob(E1 marche)  probabilité que E1 fonctionne  Re1(t ) = e - le1t
t t
- le 1T
 Probabilité que E1 tombe en panne sur l’intervalle [0, t] à l’instant T = �f (t )dt = �l e
0 e1 0 e1 .dT

 Probabilité que S marche sachant que E1 ne marche plus à partir de T = Re 2 (t > T ) = e - le 2 ( t -T )

Rs ( t ) = 1xe - le1t + ( �
t
l
0e1e - le1T
.dT xe )
- le 2 ( t -T )
= e - le1t
+ le1.�e - le 1T
.dT xe - le 2 .t
xe
t
le 2 .T
0

t t
Rs ( t ) = e - le1t + le1.e - le 2 .t .�e - le1T .e le 2 .T .dT = e - le1t + le1.e - le 2 .t .�e - ( le1 -le 2 )T .dT
0 0
t
- le 1t - le 2 .t t - ( le 1 - le 2 )T - le 1t - le 2 .t �e -( le1 -le 2 )T �
Rs ( t ) = e + le1.e .�e .dT = e + le1.e .� �
0
�-(le1 - le 2 ) � 0

- le 1t - le 2 .t �e -( le1 - le 2 )t 1 �
Rs ( t ) = e + le1.e .� - �
�-(le1 - le 2 ) -(le1 - le 2 ) �

- le 1t - le 2 .t 1 - e -( le1 - le 2 )t �
�
Rs ( t ) = e + le1.e .� �
�(le1 - le 2 ) �
le1.e - le1t - le 2 .e - le1t + le1.e - le 2 .t - le1.e -le 2 .t .e -( le1 -le 2 )t
Rs ( t ) =
le1 - le 2
le1.e - le1t - le 2 .e - le1t + le1.e - le 2 .t - le1.e - le 2 .t -le1.t +le 2 .t
Rs ( t ) =
le1 - le 2
le1.e - le1t - le 2 .e - le1t + le1.e - le 2 .t - le1.e - le1.t le1.e - le 2 .t - le 2 .e - le1t
Rs ( t ) = =
le1 - le 2 le1 - le 2
Si on prend en compte l’élément de détection et de commutation DC, on obtient alors :

- lDC .t le1.e - le 2 .t - le 2 .e - le1t

Rs ( t ) = e .
le1 - le 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 - lDC .t .e - l .t .(1 + l.t )
�i = n -1
(l.t )i �
Pour n éléments de taux de défaillance identiques montés en //, on trouve : Rs ( t ) = e
- ( lDC + l ).t
. �� �
�i =0 i ! �

Chapitre 3 LA FIABILITE Page 49

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
43 – Redondance majoritaire :
La redondance majoritaire est telle que la fonction est assurée si au moins la majorité des éléments est en état de
fonctionnement.
Cette redondance concerne surtout des signaux de grande sécurité, et en particulier les équipements électroniques. Le
signal de sortie est celui de la majorité des composants. Le cas le plus simple comporte 3 éléments.
On considère que l’organe D de décision a une fiabilité égale à 1.
E1
RS=probabilité d’avoir plus de 2 éléments en fonctionnement correct
Si Re1=Re2=Re3=R
E2 D k =3
RS = �C3k .R k .(1 - R )3-k = 3R 2 - 2R 3
k =2
E3
Si on généralise à n (impair obligatoirement pour avoir une majorité) éléments, on obtient :
k =n
n +1
RS = �Cnk .R k .(1 - R )n -k avec c =
k =c 2
La formule de calcul de « c » permet d’obtenir la majorité des éléments.
En tenant compte de la fiabilité du composant de décision :
k =n
n +1
RS = RD .�Cnk .R k .(1 - R )n -k avec c =
k =c 2
44 – Application :
Un processus est M1 M2 M3 M4 M5 T1 T2 T3
représenté par le processus 0,85 0,99 0,99 0,99 0,99 0,8 0,99 0,99
ci-contre :
La fiabilité du système entier est le produit de toutes les fiabilités élémentaires : Rs = 0,64
Pour améliorer cette fiabilité, on peut appliquer des redondances sur les systèmes les moins fiables : M1 et T1.
Une des solutions peut consister à utiliser 3 T1 et 2 M1. Economiquement, il va de soi que cette solution coûterait trop
cher. On se contentera de redonder les éléments faibles des systèmes M1 et T1

T1

M2 M3 M4 M5 T2 T3
M1 T1
0,99 0,99 0,99 0,99 0,99 0,99

M1 T1

Rs = �
1 - (1 - 0,85)2 �
� 1 - (1 - 0,8)3 �
�x 0,99 x �
4
� �x 0,99 = 0,91  Résultat satisfaisant.
2

Chapitre 3 LA FIABILITE Page 50

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
V – ANALYSE DE LA FIABILITE PAR LA LOI EXPONENTIELLE :
51 – Définition de la loi exponentielle :
Rappel sur la durée de vie d’un matériel : On constate que durant la période de maturité d’un équipement,
tx de défaillance
λ(t) est constant ou sensiblement constant. C’est le champ
d’application de la loi exponentielle qui repose sur l’hypothèse λ
= constante.
Maturité Les défaillances émergent sous l’action de causes diverses et
indépendantes.
 Si λ=cte, alors MTBF = 1/ λ en fiabilité
t  Si μ=cte (taux de réparation), alors MTTR = 1/ μ en
maintenabilité
t

R (t ) = e �
l ( u ).du
-
0
et comme l(u ) = cte = l
t  R(t ) = e - l .t
R (t ) = e �
- l .du t
0
e - l .u �
=�
� �= e
- l .t
0

 Densité de probabilité : f (t ) = l.e - l .t

 Fonction de répartition : F (t ) = 1 - R(t ) = 1 - e - l .t
1
 Espérance mathématique : E (t ) = MTBF =
l
52 – 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 - l .t � ln R(t ) = -l.t � t = - .ln R(t ) � t = .ln
l l 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
l R(t ) R (t ) 0,9
Au bout de 211 heures, on estime donc que 90% des composants survivront.
53 – Représentation graphique de la loi exponentielle :
-l.t
Si R(t ) = e - l .t , alors ln R (t ) = -l.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
l/2
,3
1/e=0,368

m=1/l lt ou t
t
2,3/l

Chapitre 3 LA FIABILITE Page 51

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
54 – 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
 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 - l .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) l(t)

=3
1 =3 1
=1
0,5
=0,5 =0,5
=1
t =3 =1 t t

Remarque : pour γ=0 et β=1, on retrouve la distribution

exponentielle, cas particulier de la loi de Weibüll :
f(t)
1 1
l= = h1
h MTBF
avec h2 < h1

η  Paramètre d’échelle >0 qui s’exprime dans l’unité de h2

temps t

γ  paramètre de position,  - < γ < + , qui s’exprime dans l’unité de temps :

Chapitre 3 LA FIABILITE Page 52

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
 γ>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 g<0 g=0 g>0
précédé la mise en historique des TBF
t

Relations fondamentales :

 -1 �t -g �
 �t - g � -� �
 Densité de probabilité : f (t ) = . � .e �h �
avec t �g
h �h � �

�t -g �
-� �
 Fonction de répartition : �h �
F (t ) = 1 - e

�t -g �
-� �
 Loi de fiabilité : �h �
R ( t ) = 1 - F (t ) = e
Taux de défaillance :

 -1 �t -g �  -1
f (t ) f (t )  �t - g � -� �h
� 1  �t - g �
l (t ) = = = .� � .e
�
. � l(t ) = . �
R (t ) 1 - F ( t ) h � h � h �h �

�t - g �
-� � �
�h �
e
MTBF et écart type :

 -1 �t -g �
+� x  �t - g � -� x �
E (t ) = MTBF = �t .f (t ).dt = lim � t .f (t ).dt = lim � t. . � .e �h �
.dt
0 x �� 0 x �� 0 h �h � �
1
E (t ) = MTBF = g + h.G(1 + ) = Ah + g

s = Bh
Γ est une fonction mathématique complexe. A et B sont des paramètres issus de tables et sont calculés à partir de la loi Γ.

Chapitre 3 LA FIABILITE Page 53

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION

Ex : pour β=1,2, γ=0 et η=550 heures  MTBF = 0,9407x550+0≈517 heures.

Chapitre 3 LA FIABILITE Page 54

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
62 – 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é, en particulier les roulements à billes.
 1 1
�t -g �  
-� �
�h �
�t - g � 1 �t - g � t - g � 1 � � 1 �
R (t ) = e � ln R( t ) = - � � � ln =� �� =�
ln � � t = h. �
ln � +g
�h � R (t ) � h � h � R (t ) � � R( t ) �
63 – Papier Weibüll :
C’est un papier log / log qui comporte 4 axes :

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 β
64 – 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

Chapitre 3 LA FIABILITE Page 55

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
Exemple d’application :
Préparation des données :
Ordre i TBF Fi
1 165 0,11
2 330 0,26
3 515 0,42
4 740 0,58
5 915 0,73
6 1320 0,89
Tracé du nuage de points :

η=770 heures

β=1,4

D2 D1

Tracé de la droite de Weibüll D1 : le tracé se fait sans difficulté « au jugé ».

Détermination des paramètres de la loi :
 Le fait d’obtenir directement une droite D1 sans faire de redressements indique que γ=0 (paramètre de
position)
 La droite D2, // à D1, passant par l’origine coupe l’axe « b » en un point β=1,4. C’est la valeur du
paramètre de forme
 La droite D1 coupe l’axe des temps à t=η=770 heures. C’est le paramètre de la loi de Weibüll
Equations de la loi :
0,4
�t �
-� �
R (t ) = e �770 �

Chapitre 3 LA FIABILITE Page 56

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
Détermination de la MTBF :
Les tables annexes donnent les valeurs de A et B pour β=1,4 : A=0,911 et B=0,660. On en déduit
MTBF = Ah + g = 0,911x 770 = 700 heures et s = Bh = 0,660 x 770 = 508 heures.
Remarque sur la forme du nuage de points :
 Si le nuage de points approxime une droite, la détermination de γ est instantanée puisque γ=0.
 Dans le cas où ce n’est pas une droite mais une courbe (concave ou convexe) qui est approximée, il existe des
méthodes de redressement de la courbe pour obtenir une droite et donc γ. Dans ce cas, l’utilisation de logiciels
spécialisés est conseillée.
VII – METHODES D’APPROXIMATION DES VALEURS DE LA FONCTION DE REPARTITION :
On dispose pour nos études de fiabilité d’un certain nombre de données expérimentales ou réelles sur les TBF ; TBF
dont on veut étudier la fonction de répartition.
Ces données représentent un échantillon « n » de la population que l’on veut appréhender. Elles doivent être classées
par ordre croissant de durée (en heures, jours, etc), suivant l’unité la plus adaptée.
i
L’estimation de la fonction de densité pour une durée ti est donnée par : f (ti ) =
n +1
Or, ce n’est pas la fonction de densité qui nous intéresse mais la fonction de répartition F(ti). Cette fonction de
répartition peut être estimée selon plusieurs méthodes dont 2 sont particulièrement applicables pour les lois de fiabilité
(exponentielle et Weibüll) : ce sont les méthodes des rangs médians et des rangs moyens. Le choix entre l’une ou
l’autre des méthodes est fonction de la taille « n » de l’échantillon.
i - 0,3
 Si n �20 , on utilise la méthode des rangs médians et F (ti ) =
n + 0,4
i
 Si n > 20 , on utilise la méthode des rangs moyens et F (ti ) =
n +1
Des tables donnent les valeurs de F(ti) directement en fonction de la taille n de l’échantillon.
Ex : table des rangs médians :

Chapitre 3 LA FIABILITE Page 57

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
VIII – APPLICATION DE WEIBULL : OPTIMISATION D’UNE PERIODE D’INTERVENTION
SYSTEMATIQUE : ABAQUES DE KELLY :
La question qui revient sans cesse dans un service maintenance pour un équipement est : faut-il choisir de garder le
correctif ou de mettre en œuvre un préventif systématique ? Pour répondre à cette question, il existe plusieurs outils
(abaques de Noiret par exemple) dont l’utilisation de la loi de Weibüll.
La mise en pratique de cette loi va permettre de répondre aux 2 questions suivantes :
 Existe-t-il une période d’intervention systématique T telle que la maintenance préventive soit plus économique
que la maintenance corrective ?
 Si oui, quelle est cette période optimisée θ ?
Cet outil d’optimisation sera nommé outil « r, β ».
71 – Mise en œuvre de la méthode :
Sur un système réparable, dont un constituant « fragile » est interchangeable, comment faire pour déterminer la période
θ de remplacement préventif ?
Il faut en 1er lieu connaître :
 La loi comportementale R(t) du constituant
 Le coût « p » du correctif qui, par hypothèse, est égal au coût de l’intervention préventive liée au remplacement
du constituant défectueux
 Le coût indirect « P » des conséquences de la défaillance
On appellera r=P/p le ratio de « criticité économique » de la défaillance. Domaine de validité : 2 < r < 100.
Evaluation du coût C1 de l’intervention corrective :
Le coût moyen d’une intervention corrective est p + P.
p+P
Le coût moyen par unité d’usage devient : C1 =
MTBF
Evaluation du coût C2(θ) d’une intervention préventive systématique :
Si θ est la période de remplacement systématique du composant, le coût aura 2 termes :
 Le coût de l’intervention p
 Le coût du correctif résiduel lié au risque de défaillance avant θ et évalué par sa probabilité F(t) avec t< θ. Ce
coût est égal à : P.F (t ) = P.(1 - R (t )) .
p + P.(1 - R (t ))
Le coût moyen par unité d’usage est donc C 2(q ) = , avec m(θ) la durée de vie moyenne des
m(q )
q
composants ne dépassant pas θ, puisqu’ils ont été changés à cette date. m(q ) = �R(t )dt .
0

Critère de choix :
C 2(q )
Le préventif systématique sera choisi s’il existe une période θ telle que C2(θ)<C1 ou encore  1.
C1
Principe de l’optimisation de θ :
C 2(q )
On étudie les variations de quand θ varie.
C1
C 2(q )
 Si > 1, alors il n’y a pas de solutions
C1
 Si le rapport à minimum inférieur à 1, la valeur t = θ correspondant au minimum est optimisée
C 2(q ) p + P (1 - R(q )) MTBF
On forme le rapport : = .
C1 m(q ) p+P

q�
�
-� �
R(θ) est modélisable par une loi de Weibüll à 2 paramètres (γ=0). h�
R(q ) = e �

Chapitre 3 LA FIABILITE Page 58

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
� 1� �
La MTBF est aussi une fonction de η et β : MTBF = h.G �
1 + �avec G = �e - t .t a -1.dt qui est une fonction
� � 0

mathématique complexe.
� 1�
Si on pose x =
q
et r =
P
, alors
C 2(q ) p + P (1 - R(q )) MTBF
= . devient : C 2( x )
=
1+ r 1- e-x ( 

) .G�
1+ �
� �
h p C1 m(q ) p+P C1 x
-t  1+ r
�e .dt 0

C 2( x )
On constate que le rapport est dépendant de 2
C1
paramètres : β qui caractérise la forme de la distribution et r,
paramètre économique, qui caractérise le rapport des coûts
indirects / directs (criticité des défaillances).
C 2( x )
Exploitation du rapport :
C1
En plus des 2 paramètres cités précédemment, le rapport fait
aussi intervenir le temps. On trace alors sur un graphique une
C 2( x )
série de courbes = f ( x ) pour des valeurs successives
C1
de θ et de r. On obtient alors des abaques (appelées abaques
de KELLY) tels que celle ci-contre :

72 – Méthodes de gestion des matériels :

Gestion individuelle en maintenance préventive
systématique : en cas de défaillance résiduelle, le
remplacement du composant défaillant initialise une nouvelle
période θ pour l’échéancier. C’est la méthode la plus
fréquente.

Défaillance
q q q q q q q

t1 t2 t3 t4 t5 t6 t7
duréeq
Gestion collective en maintenance préventive systématique : en cas de défaillance résiduelle, le remplacement du
composant défaillant ne modifie pas l’échéancier prévu.
Défaillance
q q q q q

t1 t2 t3 t4 t5 t6
q
Cette notion de gestion des équipements nous intéresse dans le cas de l’optimisation d’une période de remplacement,
puisqu’elle nous conduit à l’obtention de 2 abaques.

Chapitre 3 LA FIABILITE Page 59

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
IX – DETERMINATION EXPERIMENTALE DU TAUX DE DEFAILLANCE :
Rappels :
Le taux de défaillance, noté l(t), est un indicateur de la fiabilité. Il représente une proportion de dispositifs
survivants à un instant t.
nombre de défaillances
Sa forme générale est : . Le plus fréquemment, il s’exprime en « pannes / heure ».
durée d' usage
Attention : utilisé en fiabilité, le taux de défaillance devra exclure les défaillances extrinsèques à
l’ensemble analysé, telles que les pannes dues à des fautes de « conduite » (accidents, consignes
non respectées) ou dues à une influence accidentelle du milieu extérieur (inondation, incendie, etc.).

Exposé de la méthode :
Cette détermination nécessite un nombre important de données sur une période relativement longue de la vie des
matériels. Si ce n’est pas le cas, seule une partie de la courbe sera mise en évidence.
La méthode utilisée est celle de l’actuariat qui consiste à faire des calculs de probabilités à partir de renseignements
statistiques.
Il s’agit donc de déterminer expérimentalement le taux de défaillance l (t ) qui correspond à la probabilité d’avoir une
défaillance dans les intervalles de temps constituant la vie du matériel étudié. Une estimation de l (t ) par tranche de
temps est déterminée par le calcul suivant :
ni
lˆ(t i ) =
Ni .Dt i
 ni le nombre de défaillant durant Dt i ,
 Ni le nombre de survivants au début de la tranche t i
 Dt i = t i +1 - t i l’intervalle de temps observé
La détermination du nombre de classes doit être telle que la courbe ne soit pas trop déformée. Ce nombre dépend du
nombre total de défaillants. On peut déterminer le nombre de classes « r » tel que :

r = �n i
ou r = 1 + 3,3.log �ni
Le nombre de classes ainsi déterminé, il reste à construire la courbe en baignoire à partir des données.
La synthèse de la méthode est donnée ci-après :
Phase 1 : choix des classes :
Ex : on a 112 défaillances au bout de 1000 heures de fonctionnement. On choisit k = 10 classes de 100 heures.
Phase 2 : tableau :
Classe Nombre de machines en Cumul des temps de Nombre de Taux moyen de
fonctionnement fonctionnement défaillances défaillance

Phase 3 : exploitation :
Ce tableau permet de tracer l’histogramme des
défaillances (répartition dans le temps) et de tracer la
courbe en baignoire l(t) :

Chapitre 3 LA FIABILITE Page 60

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
Méthodologie d’étude du taux de défaillance moyen :
1. Avoir un nombre significatif d’équipements semblables (N 0), utilisés dans des conditions identiques.
2. Posséder l’historique de chaque matériel.
3. Initialiser l’origine des mesures à un temps t0.
4. Si le nombre d’équipement est supérieur à 50 alors découper en classes : le nombre de classes = N0 et
la durée d’une classe (Dt) = (durée de vie maxi – durée de vie mini) / N0.
5. Etablir un tableau permettant le calcul des l(ti+Dt)
6. Tracer l’histogramme d’évolution de l(ti+Dt)
7. Analyser l’histogramme et conclure.

Application :
Soit un bien non réparable sur lequel ont été réalisés des essais sur 55 matériels, depuis l’instant t0, pendant une
durée totale de 3490 heures.
N0 Ns(t1) Ns(t1+ Dt)
t0 t1 t1+ Dt t

 N0 : nombre initial de matériels

 Ns(t1) : nombre de survivants à l’instant t1
 Ns(t1+∆t) : nombre de survivants à l’instant t1+∆t
 ∆t : intervalle de temps
 Nd∆t : nombre de défaillants dans ∆t
 Ns(t) : nombre de survivants à l’instant t (en début de classe)

l(t+Dt)
Nd
Détermination du t t+Dt NdDt Ns(t)
nombre de
classes : 0- 500 3 55 1,09E-4 l (t + Dt ) = Dt

55 = 7, 4  7 500- 1000 8 52 3,08E-4 Ns (t ) �Dt

classes 1000- 1500 10 44 4,55E-4
La colonne Nd∆t est issue des essais
Détermination de 1500- 2000 12 34 7,06E-4 réalisés : elle correspond au nombre de
l’intervalle de 2000- 2500 7 22 6,36E-4 matériels défaillants dans chaque intervalle
classe : 3490 / 7 = de temps.
2500- 3000 8 15 1,07E-3
498,6 heures  500
heures 3000- 3500 7 7 2,00E-3

Chapitre 3 LA FIABILITE Page 61

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
X – ESTIMATIONS DES DIVERSES FONCTIONS EMPIRIQUES DE FIABILITE ET ETUDES DE
LEURS RELATIONS :
N0 : nombre d’éléments bons à l’instant t 0
Ni : nombre d’éléments bons à l’instant t i
ni : nombre d’éléments défaillants entre t i et t i +1 , noté aussi DNi
Dt i : intervalle de temps observé égal à t i +1 - t i
ni
On estime l̂ (t ) le taux de défaillance par tranche Dt : lˆ(t i ) =
Ni .Dti
n
On estime fˆ(t i ).Dt i la fonction défaillance sur l’intervalle Dt i par : fˆ(t i ).Dt i = i
N0
i

On estime Fˆ (t i ) la fonction de défaillance cumulée par : ˆ i �n i

N0 - N i N
F (t i ) = �f (t i ).Dt i = 0
= = 1- i
0 N0 N0 N0
N
On estime Rˆ (t i ) la fonction de fiabilité par : Rˆ (t i ) = 1 - Fˆ (t i ) = i
N0

ni ni
ni N .Dt N .Dt fˆ(t i )
On peut calculer alors lˆ(t i ) par : lˆ(t i ) = = 0 i = 0 i =
Ni .Dt i Ni .Dt i Ni Rˆ (t i )
N0 .Dt i N0
fˆ(t i ) fˆ(t i )
 lˆ(t i ) = et lˆ(t i ).Dt i = .Dti (relations servant au calcul des lois de fiabilité)
Rˆ (t i ) Rˆ (t i )

� �
ni 1
On peut aussi calculer la MTBF par : MTBF = �ti .f (ti ).Dti =�ti . = ( n1t1 + n2t2 + ... + ni t i + ... + n�t �) car
0 0 N0 N0
en général t0=0

Chapitre 3 LA FIABILITE Page 62

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
XI – LE ROCOF : Rate of Occurrence Of Failure :
Le « ROCOF » (Rate of OCurence Of Failure) ou « INTENSITE DE DEFAILLANCE » est un indicateur de fiabilité
équivalent au taux de défaillance.
Il ne concerne uniquement que les matériels réparables.
Analyse mathématique :
Soient des variables aléatoires (T1, T2, T3, ... Ti ) les instants de défaillance et (X1, X2, X3...Xi ) l’intervalle de temps entre
la défaillance (i-1) et la défaillance i.
N(t) est le nombre de défaillances observées sur [0, t]
On s’intéresse ici aux instants successifs de défaillance d’un matériel.
Le comportement de Xi est intéressant à étudier, en particulier ses tendances qui peuvent mettre en
évidence de mauvaises maintenances, des phénomènes de vieillissement, etc. On s’intéresse donc ici
aux instants successifs de défaillance d’un matériel (alors qu’en fiabilité, on s’intéresse plus souvent à la
première défaillance survenue sur plusieurs matériels différents mais de même conception et exploités dans
des conditions identiques).
Hypothèses :
 Le retour d’expérience débute dés la mise en service du matériel ;
 La maintenance est du type « As bad as old » (Matériel remis dans l’état qui précédait la défaillance) et non du
type « As good as old » (Matériel remis à neuf après défaillance)
Définitions du ROCOF :
Le ROCOF est estimé par le paramètre z(t) :

E[ N (t + dt ) - N (t )] d
z (t ) = lim = E[ N (t )]
dt �0 dt dt
Le ROCOF est la limite (si elle existe), du quotient de l’espérance mathématique du nombre de
défaillances d’une entité réparée, pendant un intervalle de temps [t, t + dt] par la durée de l’intervalle
de temps lorsque cette durée tend vers 0.

Estimateur : zˆ (t ) = (Nombre de défaillances observées dans [ t, t + Dt ]) / Dt

L’estimateur z(t) se construit et s’analyse de la même manière que le taux de défaillance.
Graphe de Nelson – Aalen :
L’évolution des défaillances peut être représentée par un graphe appelée graphe de Nelson-Aalen.
L’allure de l’évolution permet d’apporter un jugement sur l’état du matériel.

Chapitre 3 LA FIABILITE Page 63

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION

Détérioration du matériel
Nombre cumulé de
défaillances
Durée de vie indépendante

Amélioration du matériel

Temps cumulé

Application : changeur d’outil :

On désire analyser les défaillances d’un changeur d’outils.
Pour cela, on dispose d’un extrait de l’historique de 1999 à
2007.
Nb de Heures de
Années
défaillances fonctionnement
1999 8 2956
2000 8 3021
2001 5 2895
2002 2 1800
2003 4 2905
2004 3 2896
2005 1 1600
2006 2 3135
2007 2 3296
Nb de
On construit un tableau ayant le trame Dt
Années Défaillances z(t) ni
suivante : (heures)
ni
 z(t) : intensité de défaillance (estimateur du ROCOF) : z(t) = ni / ∆t
 ∑ni : nombre de défaillances cumulées (Graphe de Nelson Aalen)
On aboutit donc au tableau suivant :
Nb de ∆t
Années z(t) ∑ni
défaillances (heures)
1999 8 2956 2,71E-03 8
2000 8 3021 2,65E-03 16
2001 5 2895 1,73E-03 21
2002 2 1800 1,11E-03 23
2003 4 2905 1,38E-03 27

Chapitre 3 LA FIABILITE Page 64

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
2004 3 2896 1,04E-03 30
2005 1 1600 6,25E-04 31
2006 2 3135 6,38E-04 33
2007 2 3296 6,07E-04 35

Mettre en place une politique d’amélioration du matériel et de la

Matériel en fin de période de jeunesse.
maintenance.

Chapitre 3 LA FIABILITE Page 65

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION
XI – SYNTHESE :

Chapitre 3 LA FIABILITE Page 66

Chapitre 3
LA FIABILITE DES SYSTEMES DE PRODUCTION

Chapitre 3 LA FIABILITE Page 67