Calcul d'incertitude

lié au système
d'étalonnage
Note d'application
de thermocouple

Série de notes d'application

relatives aux thermocouples
Il s'agit de la troisième de quatre notes
d'application relatives aux thermocouples :
1. L'essentiel des thermocouples
2. Guide de sélection d'équipement
d'étalonnage de thermocouple
3. Calcul d'incertitude lié au système
d'étalonnage de thermocouple
Scanner de
4. Comment étalonner un thermocouple température
de précision
Super-DAQ 1586A
avec multiplexeur
Budget d'incertitude d'échantillons DAQ-STAQ

Plusieurs incertitudes d'étalonnage doivent être

évaluées pour déterminer l'incertitude totale
L'analyse d'incertitude suivante répertorie
d'un système d'étalonnage de thermocouple.
chaque incertitude à la fois sous sa forme
Cette section illustre comment préparer le
standard et étendue. L'incertitude standard a
budget d'incertitude d'un système incluant les
été multipliée par un facteur de couverture,
équipements suivants lors de l'étalonnage d'un
généralement formulé par la valeur k, pour
thermocouple à 1 000 °C :
obtenir l'incertitude étendue. Autrement,
•• Sonde de référence : thermocouple de type S l'incertitude est divisée par un diviseur de
5650 avec raccord de référence conversion pour estimer la valeur de l'incertitude
•• Afficheur : scanner de température de standard correspondante (k=1). La valeur
précision Super-DAQ 1586A avec multiplexeur k correspond à un intervalle de confiance.
DAQ-STAQ
•• Source thermique : four d'étalonnage pour
thermocouple 9118A
•• Unité testée : thermocouple de type S sans
raccord de référence
La figure 1 présente un exemple de budget
d'incertitude. Elle propose une méthode pour
regrouper certains des plus importants facteurs
d'incertitude du système. Un code est attribué
à chacune des incertitudes et est accompagné
d'un commentaire. Cette analyse d'incertitude
est complète, mais simplifiée pour être destinée
à un public susceptible de comprendre des
concepts techniques généraux. Il est donc Four d'étalonnage
important de consulter l'institut national de pour thermocouple
9118A
métrologie le plus proche ou votre laboratoire
d'étalonnage accrédité pour obtenir davantage
d'informations.
Figure 1. Évaluation de l'incertitude d'un système de calibration par thermocouple à 1 000 °C.
Facteur de
Incertitude standard de type B Code Incertitude standard (°C) Incertitude étendue (°C)
couverture
Bruit de mesure / stabilité de référence (n≥30) A1 0,01 2 0,02

Bruit de mesure / stabilité de l'unité testée (n≥30) A2 0,01 2 0,02

Mesure de l'étalon de contrôle A3 0,01 2 0,02

Précision de l'afficheur 1586A (pour la sonde de référence) B1 0,19 2 0,38

Précision du raccord de référence (pour la sonde de référence) B2 0,01 √3 0,01

Précision de l'afficheur 1586A (pour l'unité testée – Type S) B3 0,28 2 0,56

Précision du raccord de référence (pour l'unité testée – Type S) B4 0,07 2 0,13

Incertitude d'étalonnage de la sonde de référence B5 0,13 2 0,26

Limite de dérive à long terme de la sonde de référence B6 0,06 √3 0,10

Uniformité axiale du four 9118A B7 0,10 2 0,20

Uniformité radiale du four 9118A B8 0,12 2 0,23

Incertitude standard totale (k=1) 0,52

Incertitude étendue totale (k=2) 1,04

Remarque : Le résultat de chaque incertitude est arrondi au centième près.

Par exemple, dans la plupart des cas, k=2 L'équation suivante indique comment effectuer
correspond à un intervalle de confiance 2-sigma ce calcul :
ou 95 %. Plus d'informations sont disponibles s ± 0,05 °C
dans la section « Combinaison d'incertitudes » de UA1 = = = ± 0,009 °C
√n √30
ce document.
•• UA1 : incertitude liée au bruit de mesure (°C)
A1 et A2. Bruit de mesure •• s : écart type de la mesure (°C)
Le bruit de mesure constitue l'incertitude •• n : nombre d'échantillons (30 ou plus)
causée par le bruit et l'instabilité des relevés de Conseil : Si les relevés sont en unité de tension
mesure. Chaque point de mesure repose sur la (mV), divisez UA1 par la sensibilité du thermocouple
moyenne d'au moins 30 relevés ou échantillons. pour la convertir en température (°C).
L'incertitude du bruit ou de l'instabilité des
30 relevés est calculée en divisant l'écart type
des 30 échantillons par la racine carré de n A3. Mesure de l'étalon de contrôle
(n=30). Ce calcul est souvent appelé erreur type La mesure d'un étalon de contrôle est un indicateur
de la moyenne. de la stabilité du processus d'étalonnage. Il est
Généralement, le bruit de mesure couvre à recommandé de mesurer une sonde de contrôle
la fois le bruit de l'afficheur et l'instabilité de la avec l'unité testée. Les résultats de la mesure
source thermique. Les spécifications du fabricant de l'étalon de contrôle sont représentés dans un
quant à la stabilité thermique devraient être graphique de contrôle et analysés pour s'assurer
similaires, mais il est préférable d'effectuer ses qu'aucune erreur inattendue ne se produise.
propres mesures et de contrôler l'écart type réel L'étalon de contrôle doit être similaire à l'unité
pour s'assurer que les exigences de l'analyse testée, afin que les variations du processus aient
d'incertitude sont respectées. des effets similaires sur l'étalon de contrôle et
L'écart type des mesures réelles ou l'erreur l'unité testée. Un étalon de contrôle conçu avec
type de la moyenne peut être insérée dans des matériaux métalliques non nobles risque de
l'analyse d'incertitude à chaque fois qu'une ne pas garantir une répétabilité suffisante pour
mesure est relevée. Sinon, une limite de bruit établir la stabilité du processus d'étalonnage sur
de mesure peut être affectée à chaque point de le long terme. Il est recommandé d'utiliser un
mesure. Dans le cas d'une limite, les mesures thermocouple haute qualité de type S ou R comme
sont contrôlées afin de vérifier que la stabilité étalon de contrôle.
(écart type) répond à la limite établie dans La valeur entrée dans l'analyse d'incertitude
le cadre de l'analyse d'incertitude. Dans cet représente l'écart type des mesures de l'étalon de
exemple, nous utiliserons une limite d'écart type. contrôle sur la durée. Par exemple, la table de la

2 Fluke Calibration Calcul d'incertitude lié au système d'étalonnage de thermocouple

 figure 2 présente les relevés enregistrés lors de les spécifications d'un thermocouple standard
plusieurs sessions d'un processus d'étalonnage. en valeurs susceptibles d'intégrer l'analyse
L'écart type calculé à partir des sept relevés de d'incertitude.
l'étalon de contrôle est de 0,011 °C. Cette valeur Comme la précision de la plupart des
peut intégrer l'analyse d'incertitude, mais il est afficheurs de thermocouples est spécifiée en
préférable de poursuivre et d'ajouter davantage tension, l'équation suivante permet de convertir
de relevés pour mettre à jour le composant de l'incertitude de tension en température :
mesure d'étalon de contrôle jusqu'à ce qu'il U
atteigne au moins 30 relevés. Ce chiffre sert UT(°C) = v
également à établir des lignes de contrôle sur le SMJ
graphique de contrôle de l'étalon de contrôle. La •• UT(°C) : incertitude thermique équivalente
valeur de l'écart type (σ) est multipliée par deux •• Uv : incertitude de mesure de tension
ou trois pour établir des lignes de contrôle 2-σ
•• SMJ : sensibilité thermique du thermocouple au
ou 3-σ, selon vos conditions d'analyse. Le
niveau du raccord de mesure
graphique de contrôle d'échantillon de la
figure 3 présente des lignes de contrôle 3-σ. Conseil : Calculez SMJ avec une table de tension
de thermocouple. Il s'agit de la variation de
Figure 2 : Relevés de contrôle tension (mV) pour une variation de 1 °C (ΔV/ΔT).
d'échantillons
Précision de l'afficheur 1586A pour la sonde
Relevés de l'étalon de contrôle (°C) de référence :
Session 1 1 000,025 •• La spécification de précision de tension du
Session 2 1 000,014
1586A est de 0,004 % + 4 μV (% de relevé,
k=2). Remarque : Ajoutez 2 μV de plus pour les
Session 3 1 000,042 unités testées dont les mesures sont effectuées
Session 4 1 000,022 à partir d'un canal de scanner DAQ-STAQ.
Session 5 1 000,008 •• Le thermocouple est de type S.
Session 6 1 000,031 •• La température est égale à 1 000 °C.
Session 7 1 000,021 •• La tension du thermocouple de type S à
Écart type = ±0,011 °C
1 000 °C s'élève à 9,6 mV.
•• La précision de tension spécifiée pour cet
Figure 3 : Graphique de contrôle intervalle est de ±0,0022 mV (k=1). Voir
d'échantillons l'équation suivante :
9,6 mV
Vérifier le graphique de contrôle standard – 1 000 °C 0,004 × 100 + 0,004 mV
UV(k=1) = = ± 0,0022 mV
1000.060 2
1000.050 •• Le thermocouple de type S a une sensibilité de
±0,0115 mV/°C à 1 000 °C.
Vérifier le relevé standard (°C)

1000.040
•• La précision de l'afficheur est convertie en
1000.030
unités thermiques comme suit :
1000.020
± 0,0022 mV
UT(°C) = = ± 0,191 °C
1000.010
± 0,0115 mV/°C
999.990

999.980 B2. Précision du raccord de référence

(pour la sonde de référence)
31/05/2014 10/06/2014 20/06/2014 30/06/2014 10/07/2014 20/07/2014
Date
Le raccord de référence peut être externe ou
interne. Dans cet exemple, un raccord externe
B1 et B3. Calcul de précision de est utilisé. L'incertitude du raccord de référence
l'afficheur 1586A (pour sonde de externe dépend du type d'équipement utilisé.
La plupart des laboratoires utilisent des sources
référence et l'unité testée) thermiques à 0 °C, telles que des bains de glace
La précision de l'afficheur s'applique à la fois ou des puits secs, tels que le puits sec point zéro
à la sonde de référence et à l'unité testée. Les Fluke 9101. Les incertitudes liées aux bains
calculs suivants montrent comment convertir de glace sont aussi faibles que quelques mK,

3 Fluke Calibration Calcul d'incertitude lié au système d'étalonnage de thermocouple

mais requièrent des techniques méticuleuses B4. Précision du raccord de référence
et souvent complexes. Le principal laboratoire (pour l'unité testée)
d'étalonnage thermique de Fluke emploie une
approche simplifiée en utilisant un bain de glace Dans cet exemple, l'incertitude correspondante
pilée et d'eau du robinet, tout en surveillant la au raccord de référence interne est calculée à
température avec une sonde à thermistance et partir de l'équation suivante.
URJ × SRJ
un afficheur. L'incertitude liée à cette approche URJ_Effective =
est basée sur la plage thermique acceptable du SMJ
bain de glace, qui est de 0 °C ± 0,025 °C. Pour
préparer un bain de glace offrant une incertitude •• URJ_Effective : incertitude thermique effective du
moins élevée, consultez le document E563-11 raccord de référence
de l'ASTM « Standard Practice for Preparation •• URJ : incertitude du raccord de référence
and Use of an Ice-Point Bath as a Reference •• SMJ : sensibilité thermique du thermocouple
Temperature ». au niveau du raccord de mesure
L'équation suivante sert à calculer l'incertitude
•• SMJ : sensibilité thermique du thermocouple
thermique effective du raccord de référence.
au niveau du raccord de référence
Lorsque vous évaluez un thermocouple,
les deux raccords contribuent à la mesure. Calcul de l'incertitude thermique du raccord de
Ainsi, l'incertitude du raccord de référence référence interne de l'exemple suivant :
est multipliée par le rapport des valeurs de •• Le thermocouple est de type S.
sensibilité du thermocouple (Seebeck) au niveau •• La température est égale à 1 000 °C.
de la température du raccord de référence et
du raccord de mesure pour calculer l'incertitude •• La précision du raccord de référence interne
effective du raccord de référence. du 1586A est de 0,25 °C (k=2) ou de 0,13 °C
(k=1).
U ×S
U =
RJ_Effective
RJ RJ
•• Le thermocouple a une sensibilité de
S MJ
0,0115 mV/°C à 1 000 °C.
•• URJ_Effective : incertitude thermique effective du •• Le thermocouple a une sensibilité de
raccord de référence 0,006 mV/°C à 25 °C.
•• URJ : incertitude du raccord de référence •• L'incertitude thermique équivalente est
•• SMJ : sensibilité thermique du thermocouple au calculée comme suit :
niveau du raccord de mesure 0,25 °C × 0,006 mV/°C
U = = 0,13 °C
•• SMJ : sensibilité thermique du thermocouple au
RJ_Effective
0,0115 mV/°C
niveau du raccord de référence
Calcul de l'incertitude thermique du raccord de •• La formule suivante convertit l'incertitude en
référence de l'exemple suivant : incertitude standard :
0,13 °C
•• Le thermocouple est de type S. U = B4 (k=1) = 0,07 °C
2
•• La température est égale à 1 000 °C.
•• La précision du raccord de référence est de B5. Incertitude d'étalonnage de la
0,025 °C. sonde de référence
•• Le thermocouple a une sensibilité de Ce niveau d'incertitude est relevé directement
0,0115 mV/°C à 1 000 °C. sur le certificat d'étalonnage de la sonde
•• Le thermocouple a une sensibilité de de référence. Si l'incertitude à un point de
0,006 mV/°C à 25 °C. température spécifique n'est pas répertoriée sur
•• L'incertitude thermique équivalente est égale le certificat, elle devra être calculée. Consultez
à 0,013 °C. le fournisseur d'étalonnage pour obtenir des
conseils quant à la manière d'effectuer ce
0,025 °C × 0,006 mV/°C
U =
RJ_Effective = 0,013 °C calcul. Une approximation sûre, mais souvent
0,0115 mV/°C surestimée, est de choisir l'incertitude la
plus élevée des deux points de température
•• La formule suivante convertit l'incertitude en
rapportés.
incertitude standard :
0,013 °C
Si la sonde de référence est un thermocouple,
U B2 (k=1) = = 0,007 °C assurez-vous de prendre en compte
√3 l'inhomogénéité. Si l'incertitude d'étalonnage ne
prend pas en compte l'inhomogénéité, il

4 Fluke Calibration Calcul d'incertitude lié au système d'étalonnage de thermocouple

est alors important d'inclure un composant Figure 4 : Spécifications du fabricant
d'incertitude supplémentaire. Consultez un du 9118A
laboratoire d'étalonnage réputé ou l'institut
national de métrologie pour obtenir plus Homogénéité radiale
d'informations sur ce sujet. Température 9118A (14 mm (0,6 po) du 9118A-ITB de cavité
point central géométrique) à cavité
300 °C ±0,5 °C ±0,1 °C
B6. Dérive à long terme de la sonde 700 °C ±0,5 °C ±0,20 °C
de référence 1 200 °C ±0,5 °C ±0,25 °C
La dérive à long terme de la sonde de référence Homogénéité axiale
est difficile à déterminer. Les fabricants sont Température 9118A (longueur axiale 9118A-ITB (à 60 mm
susceptibles de spécifier la dérive à long terme de ±30 mm (1,2 po) à (2,4 po) à immersion
en utilisant une valeur basée sur un certain partir du point central complète)
géométrique)
nombre d'heures d'utilisation à une tempéra-
Plage totale ±0,25 °C ±0,2 °C
ture maximale. Cette méthode prend seulement
en compte le vieillissement du matériau et ne Stabilité de la température
couvre pas d'autres sources de dérive, telles que Spécifications 9118A 9118A-ITB
la manipulation, les chocs mécaniques et les Stabilité ±0,2 °C ±0,1 °C
cycles thermiques. Il est pertinent de prendre Durée de 2 heures, plage totale 2 heures à 700 °C
en compte la spécification du fabricant relative stabilisation maximum, 2 heures au
delà de 700 °C
à la dérive à long terme, mais il est tout aussi Remarque : Stabilité de la température mesurée à 2 sigmas sur
important de surveiller la sonde afin de détecter 30 minutes
toute dérive.
•• La dérive à long terme autorisée est de Combinaison d'incertitudes
±0,1 °C (divisez par la racine carrée de 3 Afin de combiner tous les composants
pour convertir à ±0,06 °C, k=1). d'incertitude, ils doivent être d'abord convertis
en incertitudes standard. Pour ce faire, il est
B7 et B8. Uniformité axiale et radiale important de comprendre le type (la distribution)
du four de chaque incertitude, afin d'effectuer le bon
calcul de conversion. Dans cette analyse
Les principaux facteurs qui affectent l'analyse d'incertitude, le tableau comporte deux types
d'incertitude et qui proviennent du four sont la d'incertitude : la distribution normale et la
stabilité et l'uniformité. La stabilité est déjà prise distribution rectangulaire.
en compte dans les mesures de bruit A1 et A2, Une incertitude basée sur une distribution
car l'incertitude correspondante décrit le bruit normale tend à concentrer la plupart des
général du système d'étalonnage, y compris résultats au centre. La distribution a donc
la stabilité du four. L'uniformité est spécifiée la forme d'une cloche. Dans le cas d'une
de deux façons : radiale (de cavité à cavité) et distribution rectangulaire, le résultat a la
axiale (gradient sur plusieurs profondeurs du même probabilité de se situer n'importe où sur
bloc d'équilibre). Ces valeurs sont spécifiées la limite. La distribution a donc la forme d'un
par le fabricant et sont répertoriées dans la rectangle. Le plus souvent, une distribution
figure 4. Vous pouvez vous fier de ces valeurs rectangulaire est affectée aux limites de
pour établir un budget d'incertitude préliminaire, processus, tels que la dérive autorisée à
mais il est préférable de les vérifier. Dans long terme. Si la distribution est inconnue, il
certains cas, l'incertitude peut être réduite, car est sage de supposer que la distribution est
les performances de l'unité sont meilleures que rectangulaire.
celles spécifiées par le fabricant. Pour estimer Pour convertir une distribution normale
les valeurs d'uniformité à 1 000 °C, utilisez les étendue (k=2 ou k=3) en incertitude standard
valeurs exprimées pour 700 °C et 1 200 °C, puis (k=1), divisez-la simplement par la valeur de
interpolez de façon linéaire. k. Pour convertir une incertitude rectangulaire
•• Uniformité spécifiée par le fabricant à en forme standard, divisez-la par la racine
1 000 °C : axiale ±0,2 °C (k=2) ou ±0,1 °C carrée de 3. Pour convertir d'autres types de
(k=1) et radiale ±0,23 °C (k=2) ou ±0,12 °C distribution d'incertitude, consultez le GUM
(k=1) (Guide to the Uncertainty in Measurement, guide
d'incertitude de mesures) ou d'autres sources
fiables d'informations quant à la façon de
calculer les incertitudes.

5 Fluke Calibration Calcul d'incertitude lié au système d'étalonnage de thermocouple

 Enfin, pour combiner les incertitudes
standard et calculer l'incertitude totale, il est
important de savoir si chaque incertitude
est indépendante ou liée. Une incertitude
indépendante n'a rien en commun avec les
autres incertitudes. Toutes les incertitudes
indépendantes peuvent être combinées ou
regroupées par la somme des carrés. Un
exemple de somme des carrés est indiqué
ci-dessous. Les incertitudes liées ont un certain
type de corrélation ou de relation avec une ou
plusieurs autres incertitudes. Les incertitudes
liées peuvent souvent être regroupées par une
simple addition, puis la somme est ajoutée par
somme des carrées aux autres incertitudes.
Parfois, des techniques plus sophistiquées de
calcul de corrélations sont utilisées, mais ce
n'est pas le cas dans cette analyse d'incertitude.
Dans cette analyse d'incertitude, les
incertitudes de précision de l'afficheur 1586A
sont liées, car le même instrument sert à
mesurer à la fois la sonde de référence et l'unité
testée. C'est pourquoi elles sont regroupées par
une simple addition, puis combinées aux autres
par une somme des carrés. Une fois l'incertitude
standard totale calculée, elle est multipliée par
le facteur de couverture souhaité, k=2 (95 %),
pour obtenir l'incertitude étendue totale.

Somme des carrés

U= √B1 2
+ B22 + B32 + B42 + B52 + B62...

Fluke Calibration. Precision, performance, confidence.™

Electrique RF Température Pression Débit Logiciel

Fluke Calibration Fluke Europe B.V.

PO Box 9090, PO Box 1186, 5602 BD
Everett, WA 98206, États-Unis. Eindhoven, Pays-Bas

Pour plus d’informations, contactez-nous :

Depuis les États-Unis : tél. (877) 355-3225 ou fax (425) 446-5116
Depuis l’Europe/le Moyen-Orient/l’Afrique : tél. +31 (0) 40 2675 200 ou
fax +31 (0) 40 2675 222
Depuis le Canada : tél. (800)-36-FLUKE ou fax (905) 890-6866
Depuis un autre pays : +1 (425) 446-5500 ou fax +1 (425) 446-5116
Site Internet : http://www.flukecal.com
©2015 Fluke Calibration.
Les caractéristiques sont susceptibles d’être modifiées sans préavis.
Imprimé aux États-Unis 7/2015
Pub-ID 13348-fre

La modification de ce document n’est pas permise sans l’autorisation écrite

de Fluke Corporation.

6 Fluke Calibration Calcul d'incertitude lié au système d'étalonnage de thermocouple