Fiabilité et Maintenance Industrielle

QUA192

Semaine 3 (Cours 5 & 6)

Estimation de la fiabilité

Roger SERRA

Département de Génie Mécanique

 ETE2013 – QUA192 par Prof. R. SERRA

Semaine 2: Estimation de la fiabilité des

équipements (Cours 3 & 4)
• La loi exponentielle
• Estimation du taux de défaillance et de la durée de vie des
équipements
• Observation en exploitation
Données complètes
Données incomplètes
Tronquées en temps
Tronquées en défaillances
Taux de défaillance maximum
Composants réparables
Observées par intervalles de temps
Modèle de défaillance nulle
• Applications de la loi de Weibull
2
 ETE2013 – QUA192 par Prof. R. SERRA

Semaine 3: Estimation et amélioration de

la fiabilité des équipements (Cours 5 & 6)
• Méthodes graphiques pour déterminer la loi des
modèles :
 Loi exponentielle
 Tracé de hasard
 Loi de Weibull
 Tracé de hasard
 Fonction de répartition
• Amélioration de la fiabilité par les redondances
 Décomposition d’un système
 Configurations en :
 Série
 Parallèle
 Mixte
• Correction TP1 & INTRA
3
 ETE2013 – QUA192 par Prof. R. SERRA

Estimation de la fiabilité – Méthodes graphiques

Objectifs :
 Calcul de la fonction cumulative de hasard (évolution du taux de
défaillance)
 Calcul du pourcentage cumulatif de défaillances (fonction de
répartition)
L'application de l'une de ces fonctions sur des papiers graphiques à
échelle fonctionnelle, appropriés à chaque loi donne une droite si la
loi est la bonne.
F (x)
f (x)
1.0

P( x ≤ a)

x=a 4
 ETE2013 – QUA192 par Prof. R. SERRA

Estimation de la fiabilité – Loi exponentielle

Méthode pouvant être utilisée quelque
soit le type de données :
1. Classer les échantillons par ordre
croissant de durée de vie (TBF).
2. Attribuer des rangs k à chaque
échantillon correspondant à ces
TBF dans l'ordre inverse.
3. Calculer la valeur de hasard =
100/k pour chaque rang k des
éléments défaillants seulement
4. Calculer les valeurs de hasard
cumulatif.
5. Si l’essai a été arrêté pour
certains échantillons, mais qu’ils
n’ont pas failli, on tiendra compte
de leur rang, mais pas du hasard
cumulatif dans le calcul. 5
 ETE2013 – QUA192 par Prof. R. SERRA

Estimation de la fiabilité – Loi exponentielle

Tracer ces valeurs de hasard cumulatif en fonction des TBF sur un
papier de loi exponentielle.(On utilise un papier Team 13-011)

Si on obtient une droite, la loi est exponentielle.

Le MTTF représente le temps correspondant à un hasard cumulatif

de 100% (ou à un pourcentage % de défaillance de 63.2%).

L'échelle de pourcentage cumulatif de défaillance peut être utilisée

pour déterminer le temps avant qu'un certain % de la population ne
fasse défaut.
6
 ETE2013 – QUA192 par Prof. R. SERRA

Estimation de la fiabilité – Loi exponentielle

7
 ETE2013 – QUA192 par Prof. R. SERRA

Estimation de la fiabilité
 Si une droite ne peut être interpolée sur le graphe exponentiel à partir des données,
alors, il faut considérer une autre loi appropriée.

8
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.1

Huit (8) circuits imprimés sont mis en test et on note leur durée de bon
fonctionnement (TBF), soit :
148 h, 200 h, 288 h, 890 h, 50 h, 100 h, 400 h, 550 h.

Vérifiez que la fonction de défaillance répond à une loi exponentielle et

calculez la durée de vie moyenne.

Réponse:
La loi est exponentielle, MTTF = 320 heures.
9
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.2
L’historique de 20 joints en opération est examiné afin de déterminer
leur durée de vie. On a remarqué les fuites sur 13 joints après les
durées d’utilisation suivantes ( en mois) : 32, 39, 58, 66, 70, 75, 88,
106, 109, 130, 155, 185, 210. Les 7 joints restants n’ont pas montré
de défaillances et les observations ont été stoppées, pour raisons
diverses, après les temps d’utilisation suivants (en mois) : 65, 75,
88, 94, 102, 110, 150 mois. Aucun bris n’a été constaté sur ceux-ci.
Vérifiez si la distribution est exponentielle.
Réponse: La loi n’est pas exponentielle.

10
 ETE2013 – QUA192 par Prof. R. SERRA

Estimation de la fiabilité – Loi de Weibull

La technique du traçage est la même que pour la loi exponentielle.

β est la pente de la courbe. On obtient sa valeur en faisant passer

une droite par le point de référence parallèle à celle obtenue. On lit
sa valeur sur l'échelle de forme.

On peut tirer du graphe d'autres informations telles que le % de

défectuosité après un certain temps (établissement des garanties).

θ correspond au temps pour 63.2% de défaillances ou 100% de

hasard cumulatif. Attention, ce n'est pas le MTBF ici !

11
 ETE2013 – QUA192 par Prof. R. SERRA

Point de référence
θ

12
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.3
En utilisant les données de l’exemple précédent,

Calculez les paramètres en utilisant un modèle Weibull

Déterminez le pourcentage de défectuosités qui seront apparues

au terme de la garantie qui est de 45 mois.

Réponses:
β = 2.12,
θ = 130 mois;
10% de défectuosités au bout de 45 mois.

13
 ETE2013 – QUA192 par Prof. R. SERRA

Calcul de MTBF suivant une loi de Weilbull

Il se calcule en utilisant une loi Gamma.

 1
E ( t )= MTBF = θ × Γ 1 + 
 β

où la loi Gamma s’écrit :

∞
Γ ( x)
= ∫ e dt ,
t x −1 − t
x>0
0

Avec Γ(x) est une loi Gamma définie sur le graphe suivant.

14
 ETE2013 – QUA192 par Prof. R. SERRA

Calcul de MTBF suivant une loi de Weilbull

Loi Gamma
y = 0,0729x - 1,1593x + 7,5457x4 - 25,255x3 + 45,769x - 42,388x + 16,486
6 5 2

100

10
Gamma(X)

0,1
1 1,5 2 2,5 3 3,5 4 4,5 5 5,5 6

Exercice 3.4 ; En considérant les résultats de l’exercice précédent,

calculez l ’espérance de vie. Réponse: MTBF = 115 mois.
15
 ETE2013 – QUA192 par Prof. R. SERRA

Calcul de MTBF suivant une loi de Weilbull

On trouve les papiers graphiques sur le site www.weibull.com

Les données sont rangées par ordre croissant de durée de vie.
La fonction de répartition de défaillances est estimée soit
 par la méthode des rangs médians (échantillons de petites tailles)
 par la méthodes des rangs moyens (plus couramment utilisée).

Méthode des rangs médians:

F(i)= 100 (i-0.3)/(n+0.4) en %

où n est le nombre d ’échantillons considérés et i le rang considéré.

Méthode des rangs moyens

F(i) = 100 i/(n+1) en %
16
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.5
On observe 9 roulements à billes durant leur fonctionnement habituel
pour tester une nouvelle série. On relève des défaillances après les
durées de fonctionnement suivantes (en heures) : 800, 300, 400,
200, 700, 1300, 1000, 500 et 600.
a) Déterminez la loi qui régit ces données, en utilisant les
méthodes de fonction hasard, de fonction de répartition par
rangs moyens, par rangs médians.
b) Calculez le MTBF des roulements
c) Calculez la fiabilité des roulements au bout de 600 h de
fonctionnement.

Réponses:
β = 2, θ = 700 heures, MTBF = 620 h, R(600) = 45%
17
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.5 (suite)

18
 ETE2013 – QUA192 par Prof. R. SERRA

Estimation de la fiabilité - Facteur de position positif

ν représente l’instant initial pour lequel F(t) = 0.
Si ν est différent de 0, les données ne se linéarisent
plus dans les graphiques de Weibull.
Lorsque ν est positif, on identifie facilement celui-ci
d’après l’asymptote temporelle.

19
 ETE2013 – QUA192 par Prof. R. SERRA

Estimation de la fiabilité - Facteur de position positif

Lorsque on a identifié ν, on fait le changement de variable, pour

linéariser :
t’=t-ν
On peut répéter l’opération 2 à 3 fois, jusqu’à obtention d’une
droite.
Le temps de bon fonctionnement doit être révisé :
 1
E ( t )= MTBF = θ × Γ 1 +  + ν
 β
Si cela ne fonctionne pas, cela signifie que la loi ne répond pas à
Weibull, ou bien qu’on a un mélange de lois.

20
 ETE2013 – QUA192 par Prof. R. SERRA

Estimation de la fiabilité - ν positif β

21
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.6
Une usine de rechapage de pneus fait le suivi, à l’aide de contrôle de
qualité chez ses clients, de 18 pneus en fonction du seuil de taux
d’usure prédéterminé. Les seuils d’usure ont été obtenus après les
durées d’utilisation suivants (en heures): 45200, 3800, 30000, 3500,
6000, 7000, 8000, 6500, 10000, 17000, 12000, 14000, 11000, 9000,
7500, 5600, 4300 et 4700. Déterminez le MTTF des pneus et la
fiabilité correspondante.

Réponses:
ν = 2000 h, β = 1.5, θ = 6000h, MTTF = 7400 h, R(MTTF) = 43%.

22
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.6 (suite)

23
 ETE2013 – QUA192 par Prof. R. SERRA

Estimation de la fiabilité - ν négatif

• Lorsque ν est négatif, on identifie facilement celui-ci d’après la
continuation de la droite.

• La courbe dans ce cas ne coupe pas l’échelle des temps.

• On poursuit la droite et prend la valeur à laquelle elle coupe

l’échelle du temps.

• On fait le changement de variable

t ’ = t - (-| ν |) = t + ν

• On itère le processus jusqu’à trouver une droite acceptable.

24
 ETE2013 – QUA192 par Prof. R. SERRA

Estimation de la fiabilité - ν négatif

β

25
 ETE2013 – QUA192 par Prof. R. SERRA

Exercie 3.7
Des retours de clients montrent qu’un produit a les durées de vie
suivantes :
200, 900, 1500, 1800, 1900, 2000, 2100, 2200, 2300, 2400, 2500,
2600, 2800 et 3000 heures.
Calculez les paramètres de la loi de Weibull et déterminez l’équation
de fiabilité.

Réponses:
β = 4, ν= -3000 h, θ = 8 000h et R(t)=e-((t+3000)/8000)^4
On ajoute 875 h aux données de durée de vie et on porte l’asymptote. A la
deuxième itération, on trouve 1500 h, et à la troisième itération, on trouve une valeur
de υ égale à 3000 heures. On peut alors espérer obtenir une droite à la quatrième
itération.

26
 ETE2013 – QUA192 par Prof. R. SERRA

Exercie 3.7 (suite)

27
 ETE2013 – QUA192 par Prof. R. SERRA

Amélioration de la fiabilité par les redondances

Pour améliorer la fiabilité, on peut :

 Jouer sur la technologie du composant, ou
bien
 Agencer les sous-systèmes de manière à
les rendre plus fiable.
La plupart des systèmes peuvent être modélisés:
 En série
 En parallèle ou majoritaire
 En attente
 Ou une combinaison de ces configurations.

28
 ETE2013 – QUA192 par Prof. R. SERRA

Amélioration de la fiabilité par les redondances

Système série
 La fiabilité de chaque sous système a été déterminée pour une
mission donnée.
 Un système est monté en série si le système tombe en panne
lorsque un seul de ses éléments tombe en panne.
 La fiabilité d’un système en série est égale au produit de la fiabilité
de chaque sous-système si les défaillances sont indépendantes.
n
Rsystème = ∏ Ri
i =1
 La fiabilité du système est plus petite que la plus petite des
fiabilités des composants.

R1 R2 R3 Rn

Rsys 29
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.8 – Système en série

Soit un système composé de 4 composants agencés en série dont
les composants opèrent dans leur période de maturité. Le
composant 1 suit une loi de weibull avec β = 3, θ = 20 000, ν=0. Le
composant 2 suit une loi linéaire avec a= 0.000000001 et b =
0.00005.
a) Déterminez la fiabilité minimale des composants 3 et 4 du système
si la fiabilité du système est supérieure à 0.5 après 10 000 heures
d’utilisation et que les composants 3 et 4 sont identiques et suivent
une loi exponentielle.
Réponse: 99.10%

b) Déterminez le taux de défaillance du système complet.

Réponse: 6.9 * 10-5 défaillances / heure

30
 ETE2013 – QUA192 par Prof. R. SERRA

Amélioration de la fiabilité par les redondances

Système parallèle
 Dans un système monté en parallèle, le R1
système faillit lorsque toutes les
composantes tombent en panne.
 Si les défaillances sont indépendantes R2
entre elles, la fiabilité du système est
égale à:
n
1 − ∏ Fi ou Fi =
Rsystème = 1 − Ri Ri
i =1

Fi : Fonction de défaillance cumulée du composant i.

 La fiabilité du système est plus grande Rn

que la plus grande des fiabilités des
composants.
 l'étude de l'agencement des systèmes Rsys
permet d'améliorer leur fiabilité.
31
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.9 – Système en parallèle

Déterminez la fiabilité d’un système parallèle composés de 3

composants au bout de 1000 h d’utilisation. Le système 1 suit une loi
log normale avec une moyenne de 7 et un écart-type de 2. Le système
2 suit une loi de Weibull avec β = 0.5, θ = 30 000, ν=0. Le système 3
est dans sa période de maturité avec un temps de bon fonctionnement
de 20 000 h.

Réponse: 99.61%

32
 ETE2013 – QUA192 par Prof. R. SERRA

Amélioration de la fiabilité par les redondances

Méthodologie pour l’élaboration des diagrammes de fiabilité
• Pour analyser les diagrammes de fiabilité où on retrouve des sous-systèmes avec
des configurations en série et en parallèle, on doit suivre une méthodologie
adéquate afin d’évaluer la fiabilité du système au complet à partir de la fiabilité des
sous-systèmes et des composantes.

a1
b1
a2

c
b2
a3

a4

• Les composantes a1, a2, a3 et a4 ont une fiabilité Ra et les composantes b1 et b2 ont
une fiabilité Rb. On assume comme hypothèse de base que les pannes des
différentes composantes sont indépendantes. 33
 ETE2013 – QUA192 par Prof. R. SERRA

Amélioration de la fiabilité par les redondances

Système en file d’attente
 Un seul élément fonctionne, les autres sont en attente.
 A l'avantage de supprimer le vieillissement des éléments qui ne travaillent pas.
 A l'inconvénient de nécessiter un organe de détection de panne et de
commutation.
 Le taux de défaillance de chaque composant doit être connu.

Ractif

R2

Rn
34
 ETE2013 – QUA192 par Prof. R. SERRA

Amélioration de la fiabilité par les redondances

Système en file d’attente
 Si un système a « r » composants identiques indépendants dont la distribution de
vie est exponentielle, alors la fiabilité du système (R(t)) se définit comme suit:

( λt )
r −1 j

Rsys ( t ) = e ∑
− λt

j =0 j!

 La fiabilité R(t) se calcule selon une loi de poisson de moyenne λt . Elle est égale
à la probabilité que tous les systèmes défaillent sauf un.
 La moyenne de la loi de poisson λt est égale au nombre moyen de défaillants au
temps t.
R

R
35
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.10 – Système en file d’attente

3 pompes identiques qui suivent une loi exponentielle ont une fiabilité
de 0.8 après 25 000 h et sont connectées en mode d’attente.
Déterminez la fiabilité du système après 100 000 heures.

Réponse:
99.84%

36
 ETE2013 – QUA192 par Prof. R. SERRA

Amélioration de la fiabilité par les redondances

Système r/n (défaillance partielle ou redondance majoritaire)

 Système à n composants identiques, montés en parallèle,
dont au moins r d'entre eux doivent fonctionner R1
r > (n+1)/2.
 La fiabilité du système Rr/n se calcule selon une loi
binomiale ayant pour fonction de probabilité la fiabilité du R2
composant, qui donne la probabilité que r ou plus r/n
composants tirés des n composants fonctionnent.
Rr
 La fiabilité de chaque élément doit être supérieure à 0.5
pour obtenir un gain de fiabilité.
Rn
n  n! 
Rr ∑ 
j = r  j !( n − j ) !
j n− j
R (1 − R) 
n

37
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.11 – Système majoritaire

Un train est tiré par 5 locomotives et 3 d’entre elles au moins doivent

fonctionner. Les moteurs de locomotives suivent une loi de Weibull
avec β = 5, θ = 30 000, ν=5000). Quelle est la probabilité pour que le
train fonctionne après 20 000 h d’opération.

Réponse :
99.97%

38
 ETE2013 – QUA192 par Prof. R. SERRA

Amélioration de la fiabilité par les redondances

Systèmes liés
On parle de système lié lorsqu’un sous-système est commun à plusieurs autres.
En pratique, un système est lié lorsqu’un de ses composants a plusieurs entrées.
La fiabilité Rsys se calcule par le théorème des probabilités totales.
Rsys P ( syst. fonct./ Comp. liée fonct.) × Rliée + P ( syst. fonct./ Comp. liée ne fonct. pas ) × Fliée

La méthode de résolution consiste à :

 Trouver la composante de liaison sur le parcours,
 Déterminer la fiabilité conditionnelle du système avec la composante de
liaison en fonction (court-circuit),
 Déterminer la fiabilité conditionnelle du système sans la composante de
liaison en fonction (circuit ouvert).
R1 R3 R5

R2 R4 R6
39
 ETE2013 – QUA192 par Prof. R. SERRA

Amélioration de la fiabilité par les redondances

Systèmes liés
Déterminer la fiabilité conditionnelle (P1) du système avec la composante de
liaison en fonction (court-circuit):
R1 R5

P1 = P ( syst. fonct./ Comp. liée fonct.)

R2 R6

Déterminer la fiabilité conditionnelle (P2) du système sans la composante de

liaison en fonction (circuit ouvert):
R1 R5

P2 = P ( syst. fonct./ Comp. liée ne fonct. pas )

R2 R4 R6

Fiabilité totale du système: Rsys = P1 × R3 + P2 × (1 − R3 )

40
 ETE2013 – QUA192 par Prof. R. SERRA

Amélioration de la fiabilité par les redondances

Système complexe
Combinaison de systèmes montés en série et en parallèle.
A D G C D
0.8 0.7 0.9 A 0.7 0.8
0.9

B E H
0.7 0.9 0.8

E F
0.9 0.9
B
C F I 0.95
0.9 0.8 0.7

Réponse: 87.80%
Réponse: 90.17%
41
 ETE2013 – QUA192 par Prof. R. SERRA

Exercice 3.12 – Système complexe

Dans un processus de production, on doit effectuer dans l'ordre les
opérations suivantes:

On choisit l’élément E comme composant critique. L'objectif de fiabilité du

processus est de 0.90.
Déterminez la fiabilité du système.
42
 ETE2013 – QUA192 par Prof. R. SERRA

Questions?

Commentaires?

Suggestions?

43