Estimation des incertitudes
Etape 4
La loi de propagation.
Remarque : vous avez souvent des craintes avec
cette étape…
…mais ce ne sont que des mathématiques à
appliquer !
137
Estimation des incertitudes
Méthodes d’évaluations des incertitudes
Etape 4 : la loi de propagation
La valeur du mesurande est exprimé par une formule mathématique :
y = f(x , x ,…, x )
1 2 n
Comment une incertitude sur une variable x se propage sur la valeur
de l’incertitude de y ?
A la formule mathématique précédente on
applique la loi de propagation pour
trouver les incertitudes
138
Estimation des incertitudes
Méthodes d’évaluations
Etape 4 : la loi de propagation
Les termes sont
appelés coefficient de
sensibilité
Suivant la pente de la
tangente, une incertitude
uxi peut avoir une influence
grande ou petite dans
l’incertitude uy d’où le nom
de coefficient de sensibilité
Il faut comprendre le coefficient
de sensibilité comme le quotient
de l’accroissement de la réponse
par l’accroissement du
signal d’entrée. 139
Estimation des incertitudes
Méthodes d’évaluations
Etape 4 : la loi de propagation
Pour propager l’incertitude sur x dans
i l’incertitude sur y, on considère que dans la
plage de variation de xi (u xi), la fonction y est assimilable à sa tangente :
u(dû à xi dans y)
Dans le cas général il y a plusieurs facteurs
d’influence xi . Comme un écart-type se
carré les différentes composantes !
Le mieux sera
d’illustrer tout cela
dans un exemple…
Voir diapos suivante 140
Estimation des incertitudes
Méthodes d’évaluations
Etape 4 : la loi de propagation (exemple 1/3)
(ubalance s =0,2g)
(umasse volumique ?)
(upied a coulisse=0,03mm)
141
Estimation des incertitudes
Méthodes d’évaluations
Etape 4 : la loi de propagation (exemple 2/3)
Application Numérique
Masse : m=9,6g±um (um=0,2g)
Longueur d’arête : a=2,15mm±ua (ua=0,03mm)
uc?
Application de la loi de propagation
142
Estimation des incertitudes
Méthodes d’évaluations
Etape 4 : la loi de propagation (exemple 3/3)
Illustration pour l’influence de la propagation de l’incertitude de a sur l’incertitude de ρ
Tangente à la courbe en
a=2,15mm
2,15-0,03 (a -ua) 2,15+0,03 (a +ua)
a=2,15mm
Incertitude sur ρ
résultante de Tangente à la courbe en
a=2,15mm
2,15-0,03 (a -ua) 2,15+0,03 (a +ua)
a=2,15mm 143
Estimation des incertitudes
Attendez, ne partez pas, c’est pas fini !
Il est possible que les paramètres xi soient
corrélés entre eux !
La loi de propagation est un peut plus complexe…
144
Estimation des incertitudes
Méthodes d’évaluations
Etape 4 : la loi de propagation
Si corrélation entre les causes d’incertitudes 1/2
Existence d’une corrélation : formule générale de la loi de propagation
Covariance
r(xi ,xj ) est le coefficient de corrélation
r(xi ,xj ) = 0, si totalement non corrélées
r(xi ,xj ) = 1, si totalement corrélées
145
Estimation des incertitudes
Relation entre le dessin de définition
et les incertitudes de mesures :
la pièce est-elle conforme ?
147
Notion de capabilité
Relation entre IT et Incertitudes – ISO 14253-1 (2013)
U constructeur du micromètre : 5μm 19.85 19.95
IT
Je cherche à mesurer la cote de 20, j’ai une pièce bonne si
la mesure est entre 19,85 et 19,95…
Mais mon moyen de mesure a une incertitude de 5μm…
148
Notion de capabilité
Relation entre IT et Incertitudes – ISO 14253-1 (2013)
19.845 19.855 19.945 19.955
U constructeur du micromètre : 5μm 19.85 19.95
IT
Plage d’incertitude
-U ; + U
Donc si je trouve avec mon micromètre : 19.85, la pièce est peut-être bonne ou
peut-être fausse…le raisonnement est le même sur la valeur de 19.95…on peut
149
Notion de capabilité
Relation entre IT et Incertitudes – ISO 14253-1 (2013)
19.845 19.855 19.945 19.955
U constructeur du micromètre : 5μm 19.85 19.95
Si la valeur lue sur le
IT Si la valeur lue sur le
micromètre est micromètre est
supérieure à 19.955, on
est sûr que la cote n’est est sûr que la cote n’est
pas bonne
Zone de non Zone de Zone de non
conformité conformité conformité
Si la valeur lue
sur le
entre 19.855 et Plage d’incertitude
19.945, on est
sûr que la cote -U ; + U
Zone a risque… Zone a risque…
… où l’on peut déclarer conforme une pièce qui ne l’est pas et inversement 150
Notion de capabilité
Définition Ancienne Norme NF E02-204 et CNOMO et MSA
La capabilité d’un système de mesure est donnée par le coefficient de
Capabilité du Moyen de Contrôle ( CMC ), c’est à dire la valeur du rapport
entre la tolérance affectée à la valeur de la grandeur et l’incertitude de mesure
de cette grandeur. La norme spécifie une incertitude maximale admissible.
IT = valeur de la tolérance
Rem : CMC (CNOMO : CMC >8 ; MSA CMC >15)
Remarque 1 : La valeur de 2 est admise comme limite minimale dans le cas de mesures
délicates (IT faible)
Remarque 2 : Plus l’incertitude de mesure de l’instrument est faible, plus le nombre
de pièces déclarées non-conformes diminue…mais plus la mesure coute chère
Remarque 3 : En production série, la notion de capabilité se remplacera par le Risque
Client (étude statistique), voir nouvelles normes 14253, 151
VOILA, LE COURS EST FINI, MAINTENANT
PASSONS A DES EXEMPLES
152
Estimation des incertitudes
Exemples
153
Exemple pratique de calcul des
incertitudes
Exemple : Mesure d’une cote a l’aide d’une jauge de profondeur
Jauge en acier
Coefficient de dilatation : α1 = 11.5.10 -6
m/°C/m
a : Cote à mesurer dont le nominal vaut 10mm
Pièce en alliage d’aluminium
α2= 24.10 -6 m/°C/m
Processus de mesure
154
Exemple pratique de calcul des
incertitudes
Etape 1 : Faire la liste des composants de l’incertitude avec méthode des 5M
Méthode Milieu
Résultat de
mesure
Matière Main d’oeuvre Moyens
155
Exemple pratique de calcul des
incertitudes
Exemple : Mesure d’une cote a l’aide d’une jauge de profondeur
Méthode Milieu
Nombre de mesure
Température de la pièce
Répétabilité
Défauts géométriques Résultat de
des surfaces
mesure
Défaut géométrique
Efforts de mesure
des becs Résolution
Dilatation Justesse
Habileté
Dilatation Fidélité
Matière Main d’oeuvre Moyens
156
Exemple pratique de calcul des
incertitudes
Etape 2 : Mise en équation de la valeur mesurée ( y=f(x1 ,x2 ,…,xn) )
a= ar + Cf + C j+ C T
ar= résultat de la mesure (moyennée ou non si une seule mesure)
Ici trois mesure :
avec Cqi : correction de la résolution à chaque lecture.
Cf = correction sur défaut de forme des surfaces
Cj = correction sur la justesse de la jauge (donnée certificat d’étalonnage ou norme)
CT1 = correction due à la dilatation de la jauge (CT1= arnominale x α1 x ΔT)
α1 : coefficient de dilatation de l’acier de la jauge, ΔT : écart à 20°c de la température de l’instrument
CT2 = correction due à la dilatation de la pièce mesurée (CT2= arnominale x α2 x ΔT)
α2 : coefficient de dilatation de l’acier de la jauge, ΔT : écart à 20°c de la température de l’instrument 157
Exemple pratique de calcul des
incertitudes
Etape 3 : Propagation des incertitudes
Voir formules en annexe A
Exemple pratique de calcul des
incertitudes
Etape 3 : Propagation des incertitudes
159
Exemple pratique de calcul des
incertitudes
Etape 3 : Propagation des incertitudes
On néglige l’incertitude sur le coefficient de dilatation :
A calculer avec le A calculer avec le
Est estimé coefficient de coefficient de
Sera déterminé Est imposé par la
suivant une loi de dilatation de dilatation de
par une méthode résolution de
distribution l’acier de la jauge l’acier de la pièce
de type A l’instrument
Rem : le fait de diviser par Est donné par un
3 revient a prendre certificat
Estimé suivant une loi de
l’écart-type de la d’étalonnage 160
distribution
moyenne
Exemple pratique de calcul des
incertitudes
Exemple pratique de calcul des
incertitudes
Application numérique :
= 0,083+2,80+196+81+0,005776+0,0256
Des termes sont négligeables par rapport à d’autres (notamment la dilatation)
Uc(a) = 16,73 µm
Avec k = 2 : l’incertitude élargie U = 32 µm
a = 10,00 mm ± 0.03 mm avec k = 2
162