Fiabilit

I Fiabilit dun matriel

1. Introduction

La fiabilit est ltude de la dure de vie dun matriel.

Ce chapitre tudie la variable alatoire qui chaque matriel, associe son temps de bon fonctionnement. Son
esprance mathmatique et la moyenne des temps de bon fonctionnement ou MTBF.

Lanalyse dun chantillon permet de justifier que cette variable alatoire suit une loi exponentielle ou une loi de
Weibull.

T dsigne la variable alatoire qui, tout dispositif tir au hasard, associe son temps de bon fonctionnement (Time
Between Failures) ou sa dure de vie avant une dfaillance.
Pour simplifier, nous choisissons comme origine des temps l'instant t = 0 o le dispositif choisi est mis en marche, soit pour
la premire fois, soit aprs une rparation qui la remis neuf.
Alors T mesure ainsi l'instant o apparat la premire dfaillance d'un dispositif pris au hasard dans la population considre,
partir de l'instant t = 0.

2. Fonction dfaillance. Fonction de Fiabilit.

T est une variable alatoire continue valeurs positives et possdant une densit de probabilit f.
t
Dfinition 1 : La fonction dfaillance est dfinie par : F (t ) = P(T t ) = f (t )dt pour t 0.
0
F(t) est la probabilit quun dispositif prlev au hasard dans la population considre ait une dfaillance avant linstant t.

Dfinition 2 : La fonction de fiabilit est dfinie par : R(t ) = P(T > t ) = 1 F (t ) pour t 0. R(t) est la probabilit quun
dispositif prlev au hasard dans la population considre nait pas de dfaillance avant linstant t.

F est la fonction de rpartition de la variable alatoire T.

En anglais fiabilit est traduit par reliability.

Estimation de F(t) et de R(t) :

n
2.1 Mthode des rangs bruts F ( t i ) = i
n
ni
2.2 Mthode des rangs moyens F ( t i ) =
n +1
ni 0,3
2.3 Mthode des rangs mdians F (t ) =
n i +0, 4
EXERCICE : Une usine produit des machines. On tudie la fiabilit de 25 de ces machines ds leur mise en service.
Pour cela on relve leur temps de bon fonctionnement (TBF) en labsence de toute rparation.

TBF en [0;100] ]100,200] ]200,300] ]300,400] ]400,500] ]500,600] ]600,700] ]700,800]

jours
Nombre de 6 6 5 3 0 1 2 1
dfaillances
TBF en ]800,900] ]900,1000]
jours
Nombre de 0 1
dfaillances

Notons Ni le nombre de machines non dfaillantes jusqu linstant ti ;

Ni
Ri = (avec la mthode du rang brut) est la proportion de machines non dfaillantes jusqu linstant ti ;
25
Fi =1 Ri est la proportion de machines dfaillantes linstant ti .
On obtient le tableau suivant :

Fiabilit 1 N. Makni
 i 0 1 2 3 4 5 6 7 8 9 10
ti (en 0 100 200 300 400 500 600 700 800 900 1000
jours)
Ni 25 19 13 8 5 5 4 2 1 1 0
19 13 8 5 5 4 2 1 1 0
Ri 1
25 25 25 25 25 25 25 25 25
6 12 17 20 20 21 23 24 24
Fi 0 1
25 25 25 25 25 25 25 25 25
Dfinition 3 : MTBF : Soit T la variable alatoire qui, chaque matriel, associe son temps de bon fonctionnement.
La MTBF ( Moyenne des Temps de Bon Fonctionnement) est lesprance mathmatique de T : E(T) = MTBF.

3. Taux davarie ( ou taux de dfaillance )

Dans lexemple prcdent, le taux davarie de lintervalle [200,300] est la proportion de machines qui tombent en panne dans
cet intervalle de temps, par rapport celles qui taient en tat de bon fonctionnement au dbut de la priode. Cest dire :
N 2 N 3 13 8
= 0.38 .
N2 13
0,38
Lunit du temps est le jour, on trouve le taux davarie moyen dans lintervalle [200,300] : = 0,0038 . On dit
300 200
que 0,38% des machines tombent en panne par jour.

Dfinition 4 Taux davarie moyen entre les instants t et t+h est gal :

nombre d'lments tombs en panne entre les instants t et t + h

nombre d'lments en tat de marche l'instant t

Cest la probabilit pour quun dispositif en tat de marche linstant t, tombe en panne entre les instants t et t+h. Donc
P(t < T < t + h) F (t + h) F (t ) (1 R(t + h)) (1 R(t )) R (t ) R(t + h)
= = = .
P(T > t ) R (t ) R(t ) R(t )

R(t ) R (t + h)
Dfinition 5 Taux davarie moyen par unit de temps : On divise par h :
R(t ) h
Taux davarie instantan (t ) linstants t : On calcule la limite lorsque h tend vers 0
R (t ) R (t + h ) R (t + h ) R (t ) 1 R '(t )
lim = lim = , do :
h0 R (t ) h h0 h R (t ) R (t )

Dfinition 6 Le taux davarie instantan linstant t dun matriel est la fonction finie par : (t ) =
R '(t )
R (t )
o R est la fonction fiabilit de ce matriel.
Consquences : Si la fonction est connue, la rsolution de lquation diffrentielle linaire du premier ordre sur un
intervalle convenable : R' (t ) + R(t ) = 0 donne la fonction de fiabilit R du matriel.
On peut alors en dduire la fonction de dfaillance F, puisque F(t) = 1 - R(t) , qui est la fonction de rpartition de la
variable alatoire T ; puis la densit de probabilit f de la variable alatoire T : F(t) = f(t).

Remarque: Exprimentalement (t) est une courbe en baignoire.

Pannes pannes
prcoces Vie utile d'usure

Fiabilit 2 N. Makni
4. MTBF : Moyenne des temps de bon fonctionnement

C'est lesprance mathmatique de la variable alatoire T dfinie au dpart

+ x
MTBF = E ( t ) = tf ( t )dt = lim tf ( t )dt .
0 x + 0

Elle reprsente l'esprance de vie du dispositif.

II Loi exponentielle
La loi exponentielle est la loi suivie par la variable alatoire T lorsque le taux d'avarie est constant.
Pour tout t 0 on a (t) = constante strictement positive.

Pour tout t 0 :
Fonction de fiabilit: R( t ) = e t

Fonction de dfaillance: F ( t ) = 1 e t
Densit de probabilit: f ( t ) = e t

y = R(t) y = f(t)

1 t
F(t) R(t)

1
MTBF : E(T) = 1 Ecart type : (T ) =

ln R (t ) = ln(e t ) ; lnR(t) = -t.

L galit R( t ) = e t est quivalente
Posons Y = lnR(t), on obtient Y = t .
Donc les points de coordonnes (t ; Y) tracs dans le mme repre orthogonal sont aligns sur une droite passant par O(0 ;0).

On estime que T suit une lois exponentielle si, pour un grand chantillon, les points connus de coordonnes
(ti ; lnRi) construit dans un repre orthogonal sont approximativement aligns avec lorigine.

L'utilisation de papier semi-logarithmique pour reprsenter R(t) permet de dduire si une loi est exponentielle ou non. S'il
1
s'agit d'une loi exponentielle, la MTBF = est l'antcdent de e1 0,368

Utilisation du papier semi-logarithmique
Reprenons l'chantillon des 25 machines observes au premier paragraphe, pour savoir si T suit une loi exponentielle. :
Les la couples (ti ; Ri) obtenus sont les suivants.

ti 0 100 200 300 400 500 600 700 800 900

Ri 1 0,76 0,52 0,32 0,20 0,20 0,16 0,08 0,04 0,04

Pour savoir si les points de coordonnes (ti ; ln Ri) construits dans un repre, orthogonal sont aligns, on utilise du papier
semi-logarithmique. Sur un tel papier, les ordonnes marques y ont en ralit pour valeur Y = lny.
Le papier semi-logarithmique est constitu de bandes ou modules identiques de hauteur ln10 (de ln 1 = ln 10 ) .
Soit a 1 et p un entier quelconque. Puisque ln(a .10 ) = lna + p ln 10, on passe du point de coordonnes (t; ln(a .10 ))
p p

au point de coordonnes (t ; lna) par une translation de p modules.

Fiabilit 3 N. Makni
 p p
De sorte que chaque module contient les points dont les ordonnes varient de ln(1. 10 ) ln(10.10 )
Le tableau prcdent scrit :

ti 0 100 200 300 400 500 600 700 800 900

Ri 100. 10-2 76. 10-2 52. 10-2 32. 10-2 20. 10-2 20. 10-2 16. 10-2 8. 10-2 4. 10-2 4. 10-2

Portons les points de coordonnes marques (ti; Ri) sur le papier semi-logarithmique; deux modules suffisent, car la
plus grande valeur prise par a est 100. Le nuage de 10 points obtenu est presque rectiligne, ce qui autorise estimer
que T suit une loi exponentielle.
La droite D tracer a pour quation Y = - t . Elle passe donc par l'origine 0 du repre. Notons que les
coordonnes de O sont t = 0, Y = 0. C'est donc le point marqu (0; 100. 10- 2), car ln(100 . 10- 2) = ln(1)= 0.

Construisons une droite passant au plus prs des points du nuage de points et par O.
1 1
Pour t= 
, on obtient
 .  1. Donc la MTBF = est labscisse du point D dont lordonne est


100 2 100 2
marque .10 36,8.10 2 , car ln( .10 ) = 1.
e e
1
On obtient 280 , do 0,0036 .

.

Utilisation de la calculatrice
On peut galement calculer une quation de la droite dajustement des 10 points de coordonnes (ti ; lnRi ) et du
point (0 ; 0) , par la mthode des moindres carrs. On vrifie que le coefficient de corrlation est r = - 0,99. Ce qui
confirme que les points sont presque aligns.

Fiabilit 4 N. Makni
 ln R = 0,0037t + 0,056
On trouve : R = e 0 , 0037 t + 0 ,56 R = e 0 , 0037 t e 0 , 056 or e 0 , 056 1

1
Donc R = e 0 , 0037 t , donc = 0,0037 et MTBF = = 270 jours.

III Loi de Weibull (mathmaticien sudois )

Weibull a choisi une loi sous forme de puissance avec 3 paramtres qui permettent d'obtenir les diverses situations :
dcroissante, constante et croissante.
Dfinition : On dit que la variable alatoire T suit la loi de Weibull lorsque son taux davarie est :
1
t
( t ) = pour t > ; , , sont des constantes avec > 0; > 0;

est le paramtre de reprage qui fixe linstant partir duquel on tudie la fiabilit.
Si = 0 : on tudie la fiabilit ds la premire utilisation de la machine.
Si > 0 : on tudie la fiabilit un certain temps aprs la premire utilisation de lappareil.
Consquences :
On retrouve, pour tout t >
t

Fonction de fiabilit: R ( t ) = exp


t
Fonction de dfaillance: F ( t ) = 1 exp


1 t
t
Densit de probabilit: f ( t ) = exp car F(t) = f(t).

On retrouve la MTBF et l'cart type l'aide de tables (voir formulaire)

MTBF : A + et = B

L'utilisation du papier imagin par Weibull pour reprsenter F(t) permet de dceler une loi de Weibull. Les points de
coordonnes (ti ; F(ti) ) sont aligns lorsque = 0. on retrouve alors graphiquement les valeurs de et de

Exercice ( BTS maintenance)

On tudie la dure de vie dun certain type de composants lectriques fabriqus par une usine.
On dsigne par T la variable alatoire qui, chaque composant , prlev au hasard dans la population, associe sa
dure de vie exprime en mois.
Aprs une tude statistique, on admet que T suit la loi de Weibull de paramtres :
= 0, = 2,4 = 50
Dterminer MTBF = E(t) et (T ) .

Fiabilit 5 N. Makni
Dterminer graphiquement les paramtre de la loi de Weibull :

Fiabilit 6 N. Makni
IV) Fiabilit d'un systme en fonction de ses composants
IV.1 Constituants en srie
Tous les constituants sont ncessaires au bon fonctionnement du systme.

Pour que le systme fonctionne il faut que tous les constituants fonctionnent.

La dure de vie T du systme est donc dfinie par la relation :

R(t) = Prob( T > t ) = Prob { T1 > t et T2 > t et . . . et Tn > t }

= Prob { ( T1 >t ) ( T2 > t ) . . . (Tn > t) }

Si les systmes sont indpendants alors :

R(t) = Prob ( T1 >t ) . Prob ( T2 > t ) . . . . . Prob (Tn > t)

IV.2 Constituants en parallles

Le systme ne sera dfaillant que si tous les constituants sont dfaillants.

La dure de vie T du systme est donc dfinie par la relation :

F(t) = Prob( T t ) = Prob { T1 t et T2 t et . . . et Tn t }

= Prob { ( T1 t ) ( T2 t ) . . . (Tn t) }

EXERCICES
Exercice 1 (corrig)

Taux de dfaillance (ou taux d'avarie)

ti 500 1000 1500 2000 2500 3000 4000
nombre d'lments 13 9 6 4 2 1 0
en fonctionnement
Entre t3 et t4, le nombre d'avaries est N3 N4 = 2.
N3 N4 2 1
Le taux d'avarie entre t3 et t4 est = =
N3 6 3
Entre t4 et t5, le nombre d'avaries est N4 N5 = 2.
N4 N5 1
Le taux d'avarie entre t4 et t5 est =
N4 2
Taux moyen d'avaries par unit de temps:
33,3%
t4 t3 = 500 : Le taux d'avarie tant de 33,3%, le taux moyen sera 0,07%
500
50%
t5 t4 = 500 : Le taux moyen sera = 0,1%
500

Exercice 2 ( BTS maintenance)

Lquipe maintenance a relev durant une anne les temps de fonctionnement, en heures, entre deux rglages
conscutifs dune des machines de conditionnement, des bouteilles et a obtenu les temps de bon fonctionnement,
rangs par ordre croissant :30 ; 50 ; 90 ; 130 ; 170 ; 230 ; 300 ; 410 ; 580.
1) A laide de la mthode des rangs moyens, complter le tableau suivant :

TBF ti 30 50 90 130 170 230 300 410 580

F(ti)
R(ti)
Yi=lnR(ti)
2) Dterminer le taux davarie moyen entre 300ime et la 410ime heure.

Fiabilit 7 N. Makni
 Exercice 3 ( BTS maintenance)
Des pices mtalliques de forme paralllpipdique sont fabriques par des machines pour lesquelles ltude des
temps de bon fonctionnement, exprims en mois , conduit la fonction de fiabilit R telle que :
R(t ) = (0, 05t + 1)e 0.1t .
1) Calculer R(t). Etablir le tableau de variation de R sur lintervalle 0;+ [ [
2) Calculer 10-2 prs :
a) La probabilit quune machine fonctionne plus de 10 mois.
b) La probabilit quune machine tombe en panne au cours de la 1re anne.
3) Sachant que la MTBF est donne par la formule :
+
MTBF = R(t )dt = lim R (t )dt
0 + 0

Calculer laide dune intgration par parties I = R (t )dt .
0

En dduire la MTBF. ( On rappelle que lim e 0 ,1

= 0)
+

Exercice 4 ( BTS maintenance)

Lusine Mcanix est spcialise dans la fabrication de pices mtalliques.
Les dures de vie , en jours , de 9 pices utilises dans des conditions semblables ont donn les rsultats
suivants :
53 ; 112 ; 178 ; 255 ; 347 ; 458 ; 600 ; 805 ; 1151 jours.

1) A laide de la mthode des rangs moyens, complter le tableau suivant :

TBF ti
F(ti)
R(ti)
Yi=lnR(ti)

2) a) Tracer le nuage de points M( ti ; R(ti) ) sur du papier semi-logarithmique. En dduire que la variable
alatoire mesurant la dure de vie des pices suit une loi exponentielle. Dterminer graphiquement la
MTBF.
b) En dduire les paramtres de cette loi, ainsi que lexpression de R(t).
3) a) A laide de la calculatrice, dterminer le coefficient de corrlation entre y et t puis une quation y = at +
b ; de la droite dajustement.
b) En dduire lexpression de R(t) en remplaant a et b par leurs valeurs respectives arrondies 10-3 prs.
4) Dterminer graphiquement et par le calcul linstant t0 o la fiabilit dune pice est 75 %.

Exercice 5 ( BTS maintenance)

Une usine fabrique des engrenages. Le service de maintenance a relev leurs dures de vie en usure acclre.
Les rsultats sont consigns dans le tableau ci-dessous :

Dure de vie en heures 250 350 400 510 550 600 750 800
F(ti) en pourcentage 3 11 18 40 50 63 91 96

F(ti) est le pourcentage dengrenages hors service la date ti.

1) A laide du papier de Weibull justifier que la variable alatoire qui prend pour valeurs la dure de vie des
engrenages peut tre ajuste par une loi de Weibull. Dterminer les paramtre de cette loi.
2) Calculer la MTBF et lcart type de cette loi de Weibull.
3) Dterminer graphiquement, puis par les calculs, le temps au bout duquel 30% des engrenages sont
dfaillants.
4) Une transmission mcanique comporte une srie de trois engrenages identiques dont les dures de vie
suivent la loi prcdente et fonctionnent de faons identiques.
Quelle est la probabilit que la dure de vie dun tel systme soit au moins de 300 heures ?

Fiabilit 8 N. Makni
 Exercice 6 ( BTS)
Sur 24 appareils observs pendant une anne, on relve les nombres de jours de bon fonctionnement suivants :
101 ; 144 ; 179 ; 209 ; 236 ;262 ; 287 ; 310 ; 334 ; 357. Les autres appareils continuant fonctionner
normalement lissue de cette anne :
1) En utilisant la mthode des rangs moyens, complter le tableau suivant :

TBF ti
F(ti)

2) Dterminer laide du papier de Weibull les paramtres de la loi de Weibull ajustant cette distribution.
3) Calculer la dure moyenne des temps de bon fonctionnement pour une dfaillance admise de 80%.
4) Dterminer graphiquement et par les calculs la probabilit quun appareil fonctionne plus de 300 jours sans
dfaillance.

Exercice 7 ( BTS)

On a relev durant une priode de 500 heures la dure de vie de 19 lments identiques et on a le tableau
suivant :
Dure de vie [0;100] ]100,200] ]200,300] ]300,400] ]400,500]
Nombre 4 4 3 2 1
dlments
dfaillant

Les autres lments continuant fonctionner au bout de ces 500 heures.

1) A laide de la mthode des rangs moyens, complter le tableau suivant :

TBF ti
F(ti)
R(ti)
Yi=lnR(ti)

2) a) Tracer le nuage de points M( ti ; R(ti) ) sur du papier semi-logarithmique. En dduire que la variable
alatoire mesurant la dure de vie des lments suit une loi exponentielle. Dterminer graphiquement la MTBF.
b) En dduire les paramtre de cette loi, ainsi que lexpression de R(t).
3)a) A laide de la calculatrice, dterminer le coefficient de corrlation entre y et t puis une quation y = at + b ;
de la droite dajustement.
b) En dduire lexpression de R(t) en remplaant a par sa valeur arrondie 10-4 prs et en remplaant eb par
sa valeur approche 0,01 prs par excs.
4) Dterminer graphiquement et par le calcul linstant t0 o la fiabilit dune pice est 50 %.

Fiabilit 9 N. Makni
Fiabilit 10 N. Makni