I.
Erreurs dues à l’appareillage utilisé :
Un instrument de mesure permet d'établir une relation entre la valeur X du mesurande (grandeur
faisant l'objet de la mesure) et la valeur lue X i du résultat de la mesure.
- la fidélité
- la justesse
Figure 01 : Représentation de la fidélité et de la justesse sous forme de cible .
On peut en donner les définitions suivantes :
1- Fidélité :
Une méthode est fidèle lorsqu’elle donne toujours le même résultat ou des résultats voisins si on
la répète sur le même échantillon. Elle caractérise la dispersion des mesures X i d'une même
grandeur. On en définit l'écart type σ ou la variance σ 2 :
√
n
∑ ( X i− X )2
i=1
σ=
n−1
n
Xi
Avec X =∑ (moyenne des X i )
1=1 n
NB : La définition de σ implique un grand nombre de mesures au cours desquelles il convient de
s'assurer que le mesurande n'a pas évolué et que l'ambiance est la même.
2- Justesse (exactitude) :
Étroitesse d'accord entre le résultat d’une mesure et la valeur attendue (CIBLE ou valeur réputée
vraie). Une méthode est réputée juste quand la moyenne X d'un grand nombre de mesures X i est
1
confondue avec la valeur X du mesurande, quelle que soit la dispersion. L'erreur de justesse J est
définie par :
J= X− X
3- Erreur due à l’évolution de la température :
Ce type d'erreur est fréquent et il faut y penser constamment. On retiendra la relation qui lie la
variation dimensionnelle à l'élévation de la température.
∆ L=α . L0 . ∆T
Longueur initiale à température t0 = L0
Longueur à la température t1 = L1
Coefficient de dilatation linéaire du matériau = α en mm/mm.°C
Un solide de 1 mm de long s’allonge de α mm lors d’une élévation de température de 1°C
4- Composantes de l’erreur de mesure :
Il est important de connaître la structure des erreurs de mesure si l’on veut déterminer la valeur de
l’incertitude. Quelle que soit la grandeur d’une erreur de mesure et le nombre des paramètres qui
en seront à l’origine, celle-ci comprendra toujours deux parties distinctes.
4-1Partie systématique :
Cette forme d’erreur se répétera toujours de la même façon et dans le même sens. Elle peut être
constante, quand elle est due par exemple au défaut de dimension d’un étalon, ou évolutive, si
elle provient par exemple de la dilatation thermique de la pièce mesurée. Elle peut être minimisée
lorsque l’on connaît avec précision ses origines en réalisant les corrections appropriées sur les
résultats de la mesure.
4-2Partie aléatoire :
Cette forme d’erreur se reproduira d’une façon et dans un sens totalement imprévisibles, elle
provient de la multiplicité des paramètres indépendants qui interviennent lors de la réalisation de
la mesure. De par sa nature aléatoire, elle est souvent régie par des lois de probabilité
(paragraphe 4-3), dont on peut estimer les paramètres en utilisant des méthodes statistiques afin
de déterminer approximativement son étendue.
4-3Quelques lois de probabilité :
On peut obtenir le tracé de la probabilité P (ou f sur la figures de distribution) en utilisant la relation
suivante :
2
−1 ( x−X )
1 2 σ
2
p ( x )=f ( x )= .e
σ . √2 π
x est comprise dans un intervalle [ Xmin ; Xmax]. Cet intervalle est défini à partir des mesures
obtenues.
2
4-3.1 Loi de distribution rectangulaire :
L’emploi de cette loi suppose que la variable x ait la même probabilité de prendre n’importe quelle
valeur dans l’intervalle¿ – U ,+U ¿ .
Figure 02 : Courbe de distribution d’une loi rectangulaire
Les calculs donnent comme valeur de l’écart type d’une distribution rectangulaire équiprobable
I
dans un intervalle de largeur I ;σ = ce qui appliqué à notre problème d’incertitude donnerait :
2√ 3
2U U
u= = et donc U ≈ 1 ,73 u.
2 √ 3 √3
Dans notre cas l’interprétation que l’on peut faire de ce qui précède est la suivante : dans un
intervalle ¿ X – u , X + u¿ c’est-à-dire pour k = 1, il y aurait une probabilité p = 0,577 de trouver le
résultat de la mesure X i , dans un intervalle¿ X – 2u , X +2 u ¿ c’est-à-dire pour k = 2 cette probabilité
serait bien entendu de 1.
4-3.2 Loi de distribution triangulaire :
Les calculs donnent comme valeur de l’écart type d’une distribution triangulaire symétrique dans
I
un intervalle de largeur I ;σ = ; ce qui, appliqué à notre problème d’incertitude donnerait :
2√ 6
2U U
u= = et donc U ≈ 2 , 45 u.
2 √ 6 √6
Figure 03 : Courbe de distribution d’une loi triangulaire
L’interprétation que l’on peut faire de ce qui précède est la suivante : dans un intervalle
¿ X – u , X + u¿ c’est-à-dire pour k = 1, il y aurait une probabilité p = 0,65 de trouver le résultat de la
mesure X i ; et dans un intervalle¿ X – 2u , X +2 u ¿ c’est-à-dire pour k = 2 cette probabilité serait de
0,965 ; naturellement dans un intervalle¿ X – 3 u , X +3 u ¿ c’est-à-dire pour k = 3 la probabilité serait
de 1.
4-3.3 Loi de distribution normale :
3
Les calculs donnent comme valeur de l’écart type d’une distribution normale dans un intervalle de
I 2U U
largeur I ;σ = ; ce qui, appliqué à notre problème d’incertitude donnerait : u= = et donc
6 6 3
U =3u .
Figure 04 : Courbe de distribution d’une loi normale
L’interprétation que l’on peut faire de ce qui précède est la suivante : dans un intervalle
¿ X – u , X + u¿ c’est-à-dire pour k = 1, il y aurait une probabilité p = 0,683 de trouver le résultat de la
mesure X i ; et dans un intervalle¿ X – 2u , X +2 u ¿ c’est-à-dire pour k = 2 cette probabilité serait de
0,954 ; naturellement dans un intervalle¿ X – 3 u , X +3 u ¿ c’est-à-dire pour k = 3 la probabilité serait
de 0,9973 c'est-à-dire pratiquement 1.
4-4Evaluation d’une incertitude simple :
L’incertitude étant le moyen de prendre en compte les erreurs inévitables que l’on commet lors de
la mesure, erreurs que l’on ne connaît pas (en effet si l’on connaissait ces erreurs il suffirait alors
d’effectuer les corrections nécessaires pour obtenir la valeur vraie). En aucun cas on ne pourra
calculer une valeur exacte de l’étendue de l’incertitude. On ne pourra qu’estimer une valeur plus
ou moins proche de la réalité.
L’estimation de l’incertitude sera calculée à partir d’outils statistiques, c’est-à-dire en considérant
les résultats de plusieurs mesures Xi (échantillon) en faisant des hypothèses sur les lois de
distribution de ces mesures, et en réalisant les calculs correspondants. En général, les résultats
issus de cette méthode seront exprimés par une moyenne (Xi) et un écart type σ (Xi).
Naturellement, dans les calculs d’incertitudes par « une méthode de type A » on admettra que u
(Xi) est égale à σ (Xi). L’incertitude globale sur toutes les mesures effectuées avec un appareil de
mesure sera une fonction de ces incertitudes partielles :
Uglobale = f(Uenvironnement, Uopérateur, Urésolution, Ujustesse, Ufidélité,…)
4-4.1 Application :
Voici un exemple de la détermination de l’incertitude sur des mesures de longueur effectuées
avec un pied à coulisse (Tableau ). Ces mesures ont été réalisées par les étudiants d’une façon
totalement indépendante sur une pièce en PVC (coefficient de dilatation 78 micron/°C/m), dont la
température de la pièce à mesurer fut 20°±10°. Afin de tenir compte de la justesse du pied à
4
coulisse, supposons 10 mesures, d’une cale étalon de 60±0,0004 à 95% (Tableau ), prise avec le
même pied à coulisse.
59,14 59,22 59,7 59,82
59,42 59,72 59,7 59,56
Tableau 01 : Résultats des mesures de la pièce 1 (alésage)
60,06 60,08 60,06 60,06 60,08
60,08 60,08 60,06 60,08 60,1
Tableau 02 : Résultats des mesures d’un étalon Ø60mm
Probabilté
4 2 .3 8 6 4 2 .7 8
.1 .2 59 .3 .4 .5 .6 59 .7
59 59 59 59 59 59 59
Mesures
Figure 05 : Tracé de la courbe de distribution d’une loi normale à partir des mesures sur la
pièce
1. Fidélité (Répétabilité) :
La répétabilité est un essai rapide qui permet de se faire une idée de la fidélité de l’appareil
Répétabilité=±2urep avec urep=σ
√
n
∑ ( X i− X )2
i=1
σ= =urep =0,234254
n−1
2. Justesse :
L’erreur de justesse = moyenne des mesures de l’étalon - valeur étalon
J= X etalon− X etalon=0,074
J n’est pas un écart type.
En utilisant la loi de distribution normale, on a :
2J j
u J= = =0,02466666
6 3
5
3. Incertitude étalon :
L’étalon utilisé est de 60±0,0004 à une probabilité de 95%, donc uet =0,0004/2
uet =0,0002
4. Variation de T° de la pièce :
Une température de 20°±10°C implique une dilatation ∆L :
∆ L=α . ∆ T . X
−6
∆ L=78.10 × 10× 59,535=0,0464 mm
2∆ L ∆L
uT ° = = =0,01546666666
6 3
Nous nous limiterons à ces paramètres principaux, on pourrait en ajouter d’autres comme la
Résolution de l’appareil … etc.
Incertitude globale est la composition des écarts types (somme quadratique) :
√ 2 2 2
uGlobale = ( urep ) + ( u j ) + ( u et ) + ( uT ° )
2
uGlobale ≈ 0,2360
Incertitude de mesure élargie : ± U=± 2 ×uGlobale =± 0,4721 à 95 % (facteur d’élargissement k=2)
Ecriture de l’incertitude de mesure :
X=59,535 mm Mesure
±0,4721 mm Incertitude de la mesure
95% niveau de confiance
La mesure est: 59,535±0,4721mm à 95%
5. Conformité des mesures :
Après avoir vu les matériels de mesure et le calcul d’une incertitude simple, se pose la question
de la conformité de la cote mesurée. Toute mesure est associée à une incertitude. Lorsque la
mesure est au centre de l’intervalle de tolérance et que l’incertitude est faible, on ne se pose pas
de questions. Par contre, si la mesure est à la limite de l’intervalle de tolérance, on peut se
demander s’il n’y a pas un risque d’accepter une pièce non-conforme.
Il y a deux manières de résoudre ce problème :
I. Ancienne méthode :
6
Si l’incertitude de mesure est faible (<=1/4 IT), on accepte toutes les pièces dont la mesure est
contenue dans l’intervalle de tolérance. On prend un risque représenté sur le schéma de la Figure
.
II. Nouvelle méthode :
On doit jeter les pièces à risque. On jette donc plus de pièces mais il n’y a pas de risque d’envoyer
un mauvais produit au client.
Cette méthode est avantageuse si on a une bonne maîtrise de la métrologie. Tout investissement
en métrologie (matériel, formation etc.) sera valorisé car on éliminera moins de pièces.
Figure 06 : Conformité d’une mesure
III. Exemples :
Si on mesure une pièce, cotée 100±1 dans un dessin de définition, avec un pied à coulisse. Un
calcul d’incertitude a donné ±0,25 à 95%
1er cas : Si la mesure est 100,08, il y a 95% de chance que la pièce ait une dimension comprise
entre 99,83 et 100,33
La pièce sera déclarée conforme sans risque (acceptée)
2ème cas : Si la mesure est 100,92, il y a 95% de chance que la pièce ait une dimension comprise
entre 100,67 et 101,17
La pièce sera déclarée conforme avec risque (acceptée avec risque) selon l’ancienne
méthode ou non-conforme (refusée) selon la nouvelle méthode
3ème cas : Si la mesure est 101,02, il y a 95% de chance que la pièce ait une dimension comprise
entre 100,77 et 101,27
7
La pièce sera déclarée non-conforme avec risque (refusée avec risque) selon l’ancienne
méthode ou non-conforme selon la nouvelle méthode
4ème cas : Si la mesure est 101,3, il y a 95% de chance que la pièce ait une dimension comprise
entre 101,05 et 101,55
La pièce sera déclarée non-conforme selon l’ancienne et la nouvelle méthode.
II.Le Travail Demandé :
1-Tableaux de données :
59.82 59.70 59.22 59.14
59.56 59.70 59.72 59.42
59.00 59.01 59.50 59.03
59.30 59.77 59.60 59.44
Tableau 01 : Résultats des mesures de la pièce 1 (alésage).
60.08 60.06 60.06 60.08 60.06
60.10 60.08 60.06 60.08 60.08
60.50 60.70 60.04 60.80 60.10
60.60 60.33 60.44 60.40 60.20
Tableau 02 : Résultats des mesures d’un étalon Ø60mm .
2-Déterminer et tracé de la probabilité P :
√
n
∑ ( X i− X )2
i=1
σ=
n−1
2
−1 ( x−X )
1 2 σ
2
p ( x )=f ( x )= .e
σ . √2 π
Danc les résultats de P(x) selon les mesures est :
P X
0,440 59
0,464 59,01
0,514 59,03
0,825 59,14
1,059 59,22
1,257 59,3
8
1,401 59,42
1,402 59,44
1,364 59,5
1,270 59,56
1,181 59,6
0,903 59,7
0,903 59,7
0,844 59,72
0,696 59,77
0,557 59,82
Tracé de la courbe de distribution d’une loi normale à partir des mesures sur la pièce :
P
1.600
1.400
1.200
1.000
P
0.800
0.600
0.400
0.200
0.000
58.9 59 59.1 59.2 59.3 59.4 59.5 59.6 59.7 59.8 59.9
2-Déterminer les incertitudes de fidélité, justesse, étalant et température :
Fidélité :
Répétabilité=±2urep avec urep=σ
√
n
∑ ( X i− X )2
i=1
σ= =urep =0,28457
n−1
justesse :
J= X etalon− X etalon=60.2425−60=0.2425
En utilisant la loi de distribution normale, on a :
2 J j 0.2425
u J= = = =0,080833
6 3 3
9
étalant :
L’étalon utilisé est de 60±0,0004 à une probabilité de 95%, donc 𝑢𝑒𝑡 = 0,0004/2 , 𝑢𝑒𝑡 = 0,0002
Variation de T° de la pièce :
Une température de 20°±5°C implique une dilatation ΔL : ∆ L=α . ∆ T . X
−6
∆ L=78.10 × 5× 59 , 43=0,0231777 mm
2∆ L ∆L −3
uT ° = = =7,7259. 10
6 3
3- Déterminer l’incertitude globale et l’incertitude élargie :
l’incertitude globale :
√ 2 2 2
uGlobale = ( urep ) + ( u j ) + ( u et ) + ( uT ° )
2
√
uGlobale = ( 0,28457 ) + ( 0,080833 ) + ( 0,0002 ) + ( 7,7259. 10 )
2 2 2 −3 2
uGlobale ≈ 0,2959
l’incertitude élargie :
Incertitude de mesure élargie : ± U=± 2 ×uGlobale =± 0,5918 à 95 % (facteur d’élargissement k=2)
4-Ecrire l’incertitude de mesure :
X=59,44 mm Mesure
±0,5918 mm Incertitude de la mesure
95% niveau de confiance
La mesure est: 59,44±0,5918mm à 95%
5- Est-ce que vous remarquez sur la conformité de la cote mesurée :
OUI,La mesure est: 59,44±0,5918mm à 95% .
10
11