Simposio de Metrología 25 al 27 de Octubre de 2006

TENDENCIAS ACTUALES DE LA CALIBRACION DE INSTRUMENTOS

PARA PESAR DE FUNCIONAMIENTO NO AUTOMATICO
1
Pablo Canalejo, 2 Jorge Mendoza
1
IBSEI S. A. de C. V.; 2Laboratorio SIMMA S. A. de C. V.
1
Rayas No. 111 Col. V. Gómez 15620 México D.F.; 2 Zona B, Casa 75 A, Los Reyes Iztacala, Edo. Mex.
1
Tel (55) 57 59 08 58; Fax: (55) 55 37 45 66; 2Tel y Fax: (55) 58 83 35 18
1
[email protected]; [email protected]

RESUMEN

La ponencia presenta las tendencias actuales para la calibración de instrumentos para pesar de
funcionamiento no automático y su importancia en los laboratorios de calibración de México. Se mencionan
algunas magnitudes de influencia que pueden producir mayor impacto en la calibración, que se consideran
actualmente en la Comunidad Económica Europea, su mejor estimado y su incertidumbre para instrumentos
de bajo y mediano alcance de medición. Se presenta un ejemplo comparativo del presupuesto de
incertidumbre obtenido siguiendo los criterios actuales usados comúnmente por los laboratorios de
calibración en México y siguiendo las tendencias actuales que se presentan y analizan (Anexo A).

1 INTRODUCCIÓN 2 EL PROCEDIMIENTO DE CALIBRACIÓN

Existen pocos documentos que establecen los 2.1 Generalidades

métodos de calibración de estos instrumentos,
algunos de ellos fueron publicados hace algunos El procedimiento debe considerar la determinación del
años [1]. En general solo se consideraba como error de indicación (prueba de linealidad) y su
contribuyentes de la incertidumbre, la repetibilidad
del instrumento, la resolución del mismo y la
incertidumbre del patrón utilizado. incertidumbre de medición para un número limitado
de cargas de
Existen otras magnitudes de influencia que
contribuyen a la incertidumbre del instrumento. Se prueba y debe reproducir las operaciones de pesada
considera necesario tomar en consideración estos que realiza el usuario de manera rutinaria tanto como
factores a fin de que la incertidumbre asociada a sea posible.
su calibración sea representativa de las
características del instrumento. Se puede considerar que en términos resumidos, una
calibración consiste en aplicar al instrumento cargas
Existen ya publicaciones que tratan este tema con de prueba bajo condiciones especificadas, determinar
la finalidad de que se utilicen como guía de el error de la indicación y estimar la incertidumbre de
requisitos mínimos para la normalización medición asociada al error de indicación.
correspondiente [2], [3].
Una calibración se realiza normalmente en el lugar
En este documento solo se mencionan algunas de donde el instrumento está instalado y debe abarcar
las magnitudes que se pueden considerar de todo el intervalo de pesada desde cero hasta el
mayor impacto en la incertidumbre de la alcance máximo (Max). El cliente puede especificar
calibración. cierta parte del intervalo de pesada, o puede
especificar ciertas cargas nominales individuales para
las que requiere la calibración.

Se debe también tomar en cuenta ciertos requisitos

para que un instrumento pueda ser calibrado, entre
estas se encuentran que: el instrumento pueda ser

1
 Simposio de Metrología 25 al 27 de Octubre de 2006

identificado claramente; todas sus funciones estén Sustituyendo (3) en (1) obtenemos:
libres de contaminación o daño y operen
correctamente; las indicaciones puedan leerse con ε = IL + δIdigL + δIrep + δIexc − (I0 + δIdig0 ) − Mc (4)
facilidad sin ambigüedades; haya sido conectado
a la fuente de alimentación durante el tiempo Las tendencias actuales en la calibración de
apropiado; sea nivelado cuando aplique y haya instrumentos para pesar incluyen la expresión de los
sido excitado con carga de valor cercano a Max al resultados de la calibración mediante una ecuación de
menos una vez. aproximación lineal [2]. No obstante, para fines de
simplificación, no es objetivo de los autores hacer uso
de este recurso.
3 RESULTADO DE LA CALIBRACIÓN
En adelante abordaremos solo el caso en que la
El resultado de la calibración de un instrumento carga de referencia está conformada por pesas.
para pesar es el error de indicación (ε) que se
La masa convencional (Mc) de las pesas se puede
determina por la expresión:
expresar como:
ε = I − Mc (1) Mc = Vn + δ ± U (5)

donde: I es la indicación del instrumento y Mc es donde: Vn es el valor nominal de la carga; δ es la

el valor de masa convencional de la carga de desviación con respecto al valor nominal determinada
referencia. en condiciones convencionales y U es la
incertidumbre certificada en el Informe de calibración.
Cualquier indicación I correspondiente a una carga
de prueba es básicamente la diferencia de la Sustituyendo (5) en (4) se puede escribir:
indicación con carga IL y la indicación sin carga I0:
ε = IL + δIdigL + δIrep + δIexc − (I0 + δIdig0 ) − Vn + δ (6)
I = IL - I0 (2)
El valor de masa convencional de la carga de prueba
Las indicaciones de un instrumento están
puede estar afectado por la probable deriva de su
limitadas por su resolución, y a su vez están
valor en el tiempo (Cp), los efectos del empuje del aire
afectadas por diversos factores entre los que se
(Cemp); los efectos de convección (Cconv) y los efectos
encuentran la excentricidad; las condiciones
magnéticos (Cmag).
ambientales; las condiciones de instalación y
operación (nivelación, corrientes de aire, fuentes
Considerando lo anterior en la ecuación (6),
de calor, radiación, alimentación eléctrica,
obtenemos para el error la expresión siguiente:
vibraciones, influencias magnéticas) y la
variabilidad en los resultados de medición.
ε = IL + δI digL + δIrep + δI exc − (I0 + δIdig 0 )
Por eso, se puede afirmar que: − Vn + δ + C p + C emp + C conv + C mag (7)
I = IL + δI digL + δIrep + δI exc − (I 0 + δI dig0 ) (3) Se debe considerar que aun cuando el mejor
estimado de una corrección sea cero se debe tener
Donde IL es la lectura del el instrumento con carga; I0 en cuenta su incertidumbre estándar siempre que
es la lectura del instrumento sin carga; δIdig0 es la ésta pueda afectar de manera significativa a la
corrección por el redondeo de la indicación sin carga; incertidumbre de la calibración.
δIdigL es la corrección por el redondeo de la indicación
con carga; δIrep es la corrección debida a la repetibilidad Los errores de indicación se determinan durante la
imperfecta y δIexc es la corrección debida a la carga prueba que conocemos como prueba de exactitud,
excéntrica. que debe realizarse usando al menos 5 cargas de
prueba diferentes distribuidas simétricamente a lo
Todas estas correcciones pueden tener largo del alcance normal de pesada ajustando la
esperanzas cero [2]. indicación a cero antes de la prueba.

2
 Simposio de Metrología 25 al 27 de Octubre de 2006

4 EVALUACIÓN DE LAS MAGNITUDES de la prueba la indicación sin carga debe ser cero. La
DE INFLUENCIA Y SU carga debe aplicarse al menos 5 veces o 3 veces
INCERTIDUMBRE cuando LT ≥ 1 000 kg.

4.1 Redondeo de la indicación o resolución Para la evaluación de la incertidumbre debida a la

(δIdig0 ); (δIdigL ) repetibilidad (δIrep ) se asume una distribución normal,
estimada como:
Las indicaciones se obtienen normalmente como
múltiplos enteros de la división de escala d. En u(δIrep ) = s(I j ) (10)
este caso es recomendable asumir el mejor
estimado de ambas correcciones como cero,
considerando que la probabilidad de encontrar su donde j es el número de mediciones realizadas en la
valor entre los límites ± d0/2 o ± dL/2 prueba de repetibilidad.
respectivamente es del 100 %,
El mejor estimado de la incertidumbre estándar
d0 es la división de escala del instrumento para la asociada a la indicación del instrumento es:
carga cero y dL la división de escala del
instrumento para la carga L. u(δIrep ) = s(I j ) / n (11)

La incertidumbre estándar para δIdig0 (resolución

donde n es el numero de mediciones realizadas en la
sin carga) se estima, asumiendo una distribución
prueba de exactitud con carga.
rectangular con límites ± d0/2, mediante la
ecuación siguiente:.
Cuando se realiza solo una prueba de repetibilidad,
esta desviación estándar puede considerarse como
d0 representativa para todas las indicaciones del
u(δI dig0 ) = (8)
(2 3 ) instrumento en todo el intervalo de pesar considerado.
Cuando se determinan varias con diferentes cargas
La incertidumbre estándar para δIdigL (resolución de prueba, debe usarse el mayor de los valores.
con carga). se estima, asumiendo una distribución
rectangular con límites ± dI/2; mediante la 4.3 Excentricidad (δIexc )
ecuación siguiente:
La prueba consiste en colocar una carga de prueba
en diferentes posiciones del receptor de carga de tal
u(δIdigL ) = dL (9) modo que el centro de gravedad de la carga coincida
(2 3 ) con las posiciones indicadas en el esquema que sigue
o en posiciones equivalentes tanto como sea posible.
4.2 Repetibilidad (δIrep )
Posiciones de la carga para la prueba de
Normalmente el mejor estimado de esta excentricidad [3]
corrección es cero y se determina durante la
prueba de repetibilidad.
1 1
La prueba consiste en la colocación repetida de 3 4 3 4
una misma carga sobre el receptor de carga, bajo 2 5 5
2
condiciones idénticas de manipulación de las
cargas y el instrumento, tanto como sea posible.
La carga de prueba no tiene que ser calibrada y
debe estar concentrada en un solo punto. Las cargas de prueba deben tener al menos un valor
de masa aproximado a Max/3 y no tiene que estar
La prueba se realiza al menos con una carga de calibrada.
prueba LT. Para los instrumentos con división de
escala d constante es apropiado usar una carga El efecto de carga excéntrica puede determinarse
0,5 Max ≤ LT ≤ Max; esta carga puede reducirse si calculando las ∆Iexc,i: (diferencias entre las
LT > 0,5 Max asciende a varias toneladas. Antes indicaciones Ii obtenidas en diferentes posiciones de

3
 Simposio de Metrología 25 al 27 de Octubre de 2006

la carga con respecto a la indicación I1 en el centro Si no esta disponible la información, D puede

del receptor de carga): estimarse a partir de la calidad de las pesas, la
frecuencia y el cuidado en su uso, como un múltiplo
∆Iexc,i = Ii – I1 (12) de su incertidumbre expandida U(Mc):

El efecto de carga excéntrica (δIexc) aplicable a D = k DU(Mc ) (17)

cualquier indicación I del instrumento se basa en
el modulo de la mayor de las diferencias I∆Iexc,iI:max donde el valor de kD puede elegirse entre 1 y 3.
determinadas en la prueba de carga excéntrica, la
carga de prueba Lexc y la indicación I, y su valor se Para la evaluación de la incertidumbre de Cp, es
estima mediante la ecuación siguiente [2]: recomendable asumir una distribución simétrica con

{ }
límites ± D (distribución rectangular). La incertidumbre
δIexc = ∆Iexc,i max
/(2L exc ) I (13) es entonces:

se asume una distribución rectangular de modo u(C p ) = D (18)

que la incertidumbre estándar es: 3

4.6 Empuje del aire (Cemp)

u(δI exc ) = I ∆I exc , i max /(2L exc 3 ) (14)
El valor de esta corrección depende de la densidad
real ρ del material de la pesa de calibración, de la
4.4 Incertidumbre estándar de la masa de densidad del aire ρa en el lugar de calibración y del
referencia lugar en que se ajusta el instrumento.

La corrección δ es la desviación de la pesa dada Se pueden considerar tres situaciones típicas.

en el Informe de Calibración conjuntamente con la
incertidumbre de la calibración U y el factor de Cuando se realiza el ajuste del instrumento en el
cobertura k. La incertidumbre estándar es: mismo lugar en que se calibra, inmediatamente antes
de la calibración, la corrección se obtiene de la
siguiente expresión [2]:
u(δ) = U (15)
k
Cemp = - Mc [(ρa – ρ0)(1/ρ - 1/ρc)] (19)
Cuando la carga de prueba esta conformada por
más de una pesa patrón las incertidumbres se donde: ρ0 es la densidad convencional del aire; ρc es
suman aritméticamente para considerar las la densidad convencional de las pesas y . ρ es
correlaciones asumidas. En tal caso la expresión densidad de las pesas;
(15) se transforma en:
Si el instrumento fue ajustado en un lugar diferente
u(δ) = ∑ Ui (16)
del lugar de la calibración, en un ambiente con
k densidad del aire ρas desconocida, la corrección se
i
obtiene de la siguiente expresión:
4.5 Deriva del valor de masa convencional
Cemp = - Mc [(ρa – ρ0)(1/ρ - 1/ρc) + ∆ρas /ρc] (20)
de la carga de prueba, Cp

Para la corrección por la posible deriva de Mc Donde ∆ρas es la diferencia entre las densidades del
desde la última calibración es recomendable aire en el lugar del ajuste y en el lugar de la
asumir un valor límite D basándose en la calibración respectivamente.
diferencia en Mc obtenida a partir de dos
calibraciones consecutivas. Es aceptable asumir Si se asume que la densidad del aire en el lugar del
que el valor de esta corrección es cero y que la ajuste es la densidad convencional, lo cual es una
probabilidad de encontrar su valor entre los límites solución posible, entonces la ecuación (20) se
simétricos ± D es del 100 % transforma en la siguiente:

4
 Simposio de Metrología 25 al 27 de Octubre de 2006

Cemp = - Mc [(ρa – ρ0)/ρ] (21)

Aplicando la Ley de propagación de las
Si el instrumento se ajusta antes de la calibración, Incertidumbres [5] [6]., la varianza de la densidad del
la varianza se obtiene de la siguiente expresión: aire se determina por:

u 2 C emp = mN [u 2 (ρ a )(1/ ρ − 1/ ρ c ) + (ρ a − ρ 0 ) u 2 (ρ) / ρ 4 ]

2 2 2 2 2 2
⎛ ∂ρ ⎞ 2 ⎛ ∂ρ ⎞ ⎛ ∂ρ ⎞ (27)
u 2 ρ a = ⎜⎜ a ⎟⎟ u p + ⎜ a ⎟ u 2h + ⎜ a ⎟ u 2 t
(22) ⎝ δp ⎠ ⎝ δ h ⎠ ⎝ ∂t ⎠

Si el instrumento fue ajustado en forma 4.9 Efectos de convección (Cconv)

independiente a la calibración en un ambiente con
densidad del aire desconocida la varianza se Si la temperatura de las pesas (tp) es diferente a la del
obtiene de: ambiente (t), se manifiesta una diferencia de
temperatura ∆T = tp - t. que se evalúa usando un
⎡ ⎛1 1 ⎞
2
u2(∆ρas)⎤
termómetro junto a las pesas en el interior de su
2 u (ρ)
2
u Cemp = m u (ρa )⎜ − ⎟ + (ρa − ρ0 ) 4 +
⎢ ⎜ ⎟ ⎥
2 2
2
N estuche.
⎢⎣ ⎝ ρ ρc ⎠ ρ ρc2 ⎥
⎦
(23) Cuando se coloca la pesa sobre el receptor de carga,
∆T > 0 genera un flujo de aire y fuerzas parásitas que
Si la densidad del aire en el lugar del ajuste es provocan un cambio aparente en su masa
igual a la densidad convencional del aire, la convencional ∆Mc. El signo de ∆Mc es opuesto al
varianza se obtiene de: signo de ∆T.

u 2 C emp = m N [u 2 ( ρ a ) / ρ 2 + (ρ a − ρ 0 ) u 2 (ρ ) / ρ 4 ] (24) El mejor estimado de esta corrección puede

2 2

considerarse cero siempre que su efecto sobre la

4.7 Incertidumbre estándar asociada a la masa convencional sea despreciable: En este caso la
corrección puede ser asumida como cero y su
densidad del material de las pesas u(ρ)
incertidumbre se puede estimar asumiendo una
distribución rectangular con límites ± ∆Mc .
La incertidumbre asociada a la densidad del
material de fabricación de las pesas, se determina
La corrección se determina por la expresión [3]:
en el laboratorio o puede estimarse a partir de la
aleación del material con que se fabrica la pesa.
(Ver Tabla B.7 correspondiente al Método F2 3/4 ∆Tconv
C conv = −k v Mc − k hMc ∆Tconv (28)
recomendado en la del Anexo B de la OIML 1/ 4
R111:2004) [5].
∆Tconv

De modo que: Si ∆Tconv >0: ; kv = 215 x 10 -9 y kh = 75,4 x 10 -9.

Si ∆Tconv < 0 ; kv = 119 x 10 -9 y kh = 20,2 x 10 -9.
uρ = UTabla / 2 (25)
Y la incertidumbre, por la expresión:
4.8 Incertidumbre estándar asociada a la
densidad del aire u(ρa) C conv
u(C conv ) = (29)
3
La densidad del aire se puede evaluar usando la
siguiente expresión aproximada siguiente [5].: 4.10 Efectos magnéticos (Cmag.)

0,34848p − 0,009(hr )(e 0,062 t ) Para considerar los efectos de interacción magnética
ρa = (26) entre las pesas y los instrumentos de muy alta
273,15 + t
resolución es recomendable detectar si esta es
apreciable. Se debe pesar la pesa patrón
donde: ρa es la densidad del aire (kg/m3); h es la conjuntamente con un separador no metálico, por
humedad relativa (%); p es la presión barométrica ejemplo de plástico o de madera. El separador se
del ambiente (hPa); t es la temperatura del debe colocar primero encima de la pesa y luego
ambiente (º C). debajo de ella para obtener dos indicaciones. Si la

5
 Simposio de Metrología 25 al 27 de Octubre de 2006

diferencia de las indicaciones es desigual de cero, tendencias actuales aceptadas por la Comunidad
esto puede considerarse como una advertencia en Europea.
el informe de calibración.
Lo anterior significa que en las calibraciones de
En principio, la corrección aplicable debería instrumentos para pesar que se realizan en México no
corresponder a la diferencia de las indicaciones de toman en cuenta algunas magnitudes de influencia
la balanza. Sin embargo lo masa apropiado seria que tienen aportaciones significativas a la
eliminar el efecto indeseable. incertidumbre del resultado de la calibración por lo
que está presente una subestimación del valor
La corrección aplicable y su incertidumbre deben correcto.
ser objeto de discusiones posteriores.
Las insuficiencias mencionadas afectan la mejor
NOTA: capacidad de medición manifestada por los
laboratorios, así como la elaboración de los
No es objetivo de este trabajo considerar los presupuestos de incertidumbres que la soportan.
posibles efectos debidos a: histéresis, deriva del
instrumento debido a los cambios de temperatura REFERENCIAS
durante la calibración, no evaluación de los
dispositivos de tara y no evaluación de la [1] Becerra La calibración de instrumentos para
movilidad ni considerar las particularidades de los pesar:- Publicación CENAM
instrumentos con más de un valor de división en el [2] EA-10/18 Guidelines on the calibration of
alcance normal de pesada. nonautomatic weighing instruments (June 2005).
[3] Paschalis Tsimitras.- Calibration of electronic non
5 INCERTIDUMBRE ESTÁNDAR DEL automatic weighing instruments. Error analysis. OIML
ERROR Bulletin Number 2. April 2003.
[4] OIML D 28 Conventional value of the result of
5.1 Agrupando las incertidumbres weighing in air. 2004.
consideradas, la incertidumbre estándar [5] OIMLI R 111-1 Weights of classes E1, E2, F1, F2,
combinada asociada al error de indicación se M1, M2, M2-3. Edition 2004 (E).
calcula por: [6] Guide for the expression of the uncertainty (GUM)
ISO, 1993.
u 2 (ε ) = d 02 / 12 + d I2 / 12 + s 2 ( I j ) / n + u 2 (δI exc )
(30)
+ u 2 (δmc ) + u 2 (C emp ) + u 2 (C p ) + u 2 (C conv )

Todas las magnitudes de entrada se asumen

independientes, por lo que no se consideran las
covarianzas.

5.2 Incertidumbre expandida de la

calibración

La incertidumbre expandida del error es:

U(ε) = k u(ε) (31)

Donde k es el factor de cobertura.

6 CONCLUSIONES

En México no se cuenta en la actualidad con una

Guía de Calibración de instrumentos para pesar.
No obstante, los laboratorios acreditados siguen
un procedimiento más o menos uniforme, sin
embargo es insuficiente comparado con las

6
 Simposio de Metrología 25 al 27 de Octubre de 2006

ANEXO A Ejemplo de aplicación para una

calibración particular. Efectos de convección

Instrumento: Balanza electrónica ∆t (conv) = -4º C

Alcance máximo ( Max ): 300 g Kv = 1,19E-07
División mínima ( d ): 0,000 1g Kh = 2,02E-08

Condiciones de referencia Ver los resultados en la tabla 2

Temperatura de las pesas 20 º C
Temperatura ambiente promedio 22,5º C Incertidumbre estándar por empuje del aire
Humedad relativa promedio: 40 %
Presión barométrica promedio: 777 hPa Densidad del aire (ρa) = 0,914 kg/m3
Densidad de las pesas ( ρ ): 7 850 kg/m3 u(ρa) = 0,003 kg/m3
Densidad de las pesas ( ρ ) = 7 850 kg/m3
Prueba de carga excéntrica: u( ρ ) = 157 kg/m3
Carga de prueba: 100 g
Resultados obtenidos en gramos Ver los resultados en la tabla 2
99,999 4; 99,999 7; 99,999 2; 99,999 6; 99,999 5;
99,999 2

Prueba de repetibilidad
Carga de prueba: 100 g
Resultados obtenidos en gramos:
99,999 3; 99,999 4; 99,999 3; 99,999 4; 99,999 5;
99,999 4; 99,999 4; 99,999 5; 99,999 4; 99,999 4

Prueba de linealidad: Ver datos en la tabla 1

para 5 cargas de prueba.

Incertidumbre estándar por los patrones: Ver

resultados para 3 cargas de prueba en la tabla 2.

Incertidumbre estándar por deriva de los

patrones: Ver resultados para 3 cargas de prueba
en la tabla 2.

Incertidumbre estándar por resolución

u(d0) = 2,89E-05 g
u(dL) = 2,89E-05 g

Incertidumbre estándar por carga excéntrica

δΙexc = 0,000 4 g
u(δΙexc) = 1,154 6 x 10 –6 (este valor es relativo).
La contribución en unidades de masa se obtiene
multiplicando por el valor de la carga. Por ejemplo,
si se calcula para la carga de prueba usada (100g)

u(δΙexc) = 1,154 6 x 10 –4 g

Incertidumbre estándar por repetibilidad

s = 6,67 x 10 -5 g

7
 Simposio de Metrología 25 al 27 de Octubre de 2006

Tabla 1 Resultados de la prueba de linealidad (determinación del error de indicación del instrumento)
Carga (g) mc δmc Cemp I ε
M g G g g g
60 60,000022 0,022 4,096E-05 59,9998 -0,00018
120 119,999954 -0,046 8,192E-05 119,9998 -0,00007
180 179,999978 -0,022 1,229E-04 179,9997 -0,00016
240 240,000134 0,134 1,638E-04 239,9996 -0,00037
300 300,000068 0,068 2,048E-04 299,9993 -0,00056

Tabla 2 Presupuestos de incertidumbre (expandida con k =1) para 3 valores de carga

60 g 180 g 300 g
Magnitud Distribución grados de libertad (ν) U (k =1) U (k =1) U (k =1)
g g g
Resolución Rectangular ∞ 2,89E-05 2,89E-05 2,9E-05
Resolución Rectangular ∞ 2,89E-05 2,89E-05 2,9E-05
Repetibilidad Normal 9 6,67E-05 6,67E-05 6,7E-05
Excentricidad Rectangular ∞ 6,93E-05 2,08E-04 3,5E-04
Masa convencional Normal 100 2,70E-05 6,55E-05 7,5E-05
Empuje Normal 100 4,92E-05 1,48E-04 2,5E-04
Deriva de patrones Rectangular ∞ 1,56E-05 3,78E-05 4,3E-05
Convección Rectangular 100 1,21E-05 3,11E-05 4,9E-05
Error ε (g ) -1,2E-4 -1,6E-4 -5,6E-4
Uc (k = 2) g 2,4E-4 5,6E-4 8,9E-4
ν eff 21 10 8

Tabla 3 Comparación de resultados (incertidumbres expandidas con k=2)

Carga aplicada U tradicional U Propuesta
g mg mg
60 0,2 0,2
120 0,2 0,4
180 0,2 0,5
240 0,2 0,7
300 0,2 0,9