Desarrollo de software de calculo para determinacion de la

presion y area efectiva, en calibraciones de balanzas de peso

muerto.
Laboratorio custodio de patrones nacionales de presion.
ENAER
(02 Mayo 2016 - 02 septiembre 2016)

Carolina Martnez Guajardo

carolina.martinezg@usach.cl
Profesor patrocinante: Francisco Vivanco Avaria
Representante empresa: Angel Flores Torres

24 de abril de 2017
Indice
1. Objetivo general de la practica. 2
1.1. Objetivos especficos. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2

2. Introduccion 3
2.1. Metrologa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2.2. Incertidumbre de la medicion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.3. Red nacional de metrologa y laboratorios asociados . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.4. Laboratorio nacional de presion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

3. Estado del arte. 7

3.1. Principio de la balanza de presion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
3.2. Montaje de las balanzas de presion y metodo de flotacion cruzada . . . . . . . . . . . . . . . . . 10
3.3. Software y tratamiento de datos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

4. Descripcion resumida del trabajo llevado a cabo 16

5. Procedimiento de trabajo 17
5.1. Investigacion y reconocimiento del area . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
5.2. Exploracion de los procedimientos utilizados y deteccion de fuentes de error . . . . . . . . . . . . 17
5.3. Diseno del nuevo sistema a implementar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
5.4. Desarrollo de la planilla . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
5.5. Comprobacion del programa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
5.6. Mejoras y desarrollo de planillas para calibraciones de instrumentos de menor exactitud . . . . . 23

6. Resultados 25

7. Conclusiones 27

8. Referencias y Bibliografa 28

9. Anexo 29

1
1. Objetivo general de la practica.

Desarrollar e implementar un software de calculo, en entorno Excel, capaz de determinar el valor de la

presion aplicada en una balanza de peso muerto utilizada como patron en una calibracion por medio del metodo
de flotacion cruzada, ademas de informar el area efectiva del ensamble piston-cilindro de la balanza de peso
muerto bajo calibracion. El software debe entregar estos resultados con su respectiva incertidumbre en unidades
SI.

1.1. Objetivos especficos.

1. Estudiar y entender topicos de metrologa: terminologa, parametros de medida, instrumentos de medicion,

metodos de calibracion por variable e instrumento, incertidumbre asociada, etc.
2. Entender y aplicar normas de calibracion para variable presion.
3. Investigar principio de funcionamiento y partes de balanza de peso muerto, identificando los elementos
fundamentales y materiales que lo componen.

4. Conocer y profundizar el modelo de presion utilizado en calibraciones de balanzas de peso muerto.

2
2. Introduccion
2.1. Metrologa
El proceso de medir es un acto que se realiza cotidianamente, casi de manera inconsciente, y es en el comercio
donde mas influencia tienen las magnitudes medidas, desde la comida que se tranza en unidades de masa, los
combustibles fosiles comercializados en unidades de volumen, hasta el consumo electrico expresado en unida-
des de potencia electrica, todo se encuentra cuantificado y regularizado para conservar un sistema ordenado y
universal. A partir de esta necesidad de normalizar las mediciones y establecer un sistema de entendimiento
masivo de las cantidades utilizadas y/o tranzadas se impulsa la metrologa. Se entiende a la metrologa como
la ciencia de las mediciones, comprendiendo las determinaciones experimentales y teoricas a cualquier nivel de
incertidumbre en cualquier campo de la ciencia y la tecnologa[1] .

En la figura se tiene un diagrama que indica las interacciones entre los conceptos y propiedades que influyen
en la metrologa.

Relacion logica entre conceptos fundamentales de metrologa.

En primer lugar se tiene el concepto de cantidad, el que se refiere a un atributo de un fenomeno, cuerpo o
sustancia que puede ser distinguido cualitativamente y determinado cuantitativamente. Una propiedad es una
cantidad si esta permite un ordenamiento lineal del sistema de acuerdo con esa propiedad. En otras palabras,
una propiedad p es una cantidad si siempre podemos decir que dos sistemas que tienen la propiedad p, son
iguales en p o que un sistema es menor que el otro en p.

Se tiene luego la unidad, que corresponde a una cantidad estandarizada de una determinada magnitud fsica,
definida y adoptada por convencion o por ley. Cualquier valor de una cantidad fsica puede expresarse como un

3
multiplo de la unidad de medida. Entre las unidades fundamentales se tienen: kilogramo (masa), metro (lon-
gitud), candela (intensidad luminosa), segundo (tiempo), amperio (corriente electrica), kelvin (temperatura) y
mol (cantidad de sustancia). Para toda unidad fundamental existe un patron primario designado y ampliamente
reconocido que presenta las mas altas cualidades metrologicas y cuyo valor es establecido sin referirse a otros
patrones de la misma magnitud. Las unidades de medida no fundamentales o derivadas pueden escribirse como
combinacion de unidades fundamentales.

El objeto en medicion se compara respecto a una realizacion fsica de unidad correspondiente, por una serie
de operaciones con los miembros de un conjunto de estandares o sus equivalentes, entendiendo un estandar
o patron como una medida materializada, material de referencia o sistema de medicion destinado a definir,
realizar, conservar o reproducir una unidad, o uno o varios valores de una magnitud para servir de referencia.
La medicion que se realice del objeto debe ser trazable, es decir, debe poder ser relacionada con referencias
determinadas, generalmente patrones nacionales o internacionales, por medio de una cadena ininterrumpida de
comparaciones teniendo todas las incertidumbres determinadas, conocida como cadena de trazabilidad.[1]

La trazabilidad esta caracterizada por seis elementos esenciales:

1. Una cadena ininterrumpida de calibraciones.

2. Incertidumbre de la medicion.
3. Reporte de calibracion o de ensayo.

4. Competencia tecnica de laboratorios de calibracion.

5. Referencia a unidades SI.
6. Recalibraciones a intervalos apropiados de tiempo.

2.2. Incertidumbre de la medicion

En general, una medicion tiene factores que estan fuera del control del experimentador que dan origen a error
en el resultado de la medicion. Tradicionalmente, se considera que un error tiene dos componentes: componente
aleatoria y componente sistematica.

Un error aleatorio es aquel que se asocia a variaciones impredecibles o estocasticas, temporales y espaciales
de las magnitudes de influencia. Los efectos de estas variaciones dan origen a las variaciones en observaciones
repetidas del mensurando. Usualmente puede reducirse incrementando el numero de observaciones y su valor
esperado es cero.

Por otra parte, el error sistematico es aquel que se produce de igual modo en todas las mediciones que
se realizan de una magnitud. Puede estar originado en un defecto del instrumento, en una particularidad del
operador o del proceso de medicion, etc. Si la fuente del error es detectada es posible minimizarlo aplicando
una correccion apropiada.

Se tiene entonces que la incertidumbre de una medicion corresponde a una estimacion del posible error en
el valor estimado del mensurando y que caracteriza el intervalo de valores dentro de los cuales se halla el valor
verdadero de la cantidad medida.

La incertidumbre asociada a una medicion puede determinarse por medio de dos formas de evaluacion: tipo
A y tipo B. Ambos metodos estan basados en distribuciones de probabilidad y las componentes de incertidumbre
resultantes de cualquier tipo son cuantificadas por varianzas y desviaciones estandar.

La varianza estimada u2 que caracteriza a una componente de incertidumbre obtenida de una evaluacion tipo
A se calcula mediante series de observaciones repetidas y es la varianza estimada estadstica s2 . La desviacion

4

estandar estimada u = u2 , es entonces u = s y por conveniencia algunas veces se denomina incertidumbre
estandar tipo A. Por otra arte, para una componente de incertidumbre obtenida de una evaluacion tipo B, la
varianza estimada u2 se evalua mediante el uso de la informacion disponible en literatura, documentacion o
cualquier fuente informativa de confianza, y la desviacion estandar estimada u algunas veces es llamada in-
certidumbre estandar tipo B. Generalmente al tipo de incertidumbre tipo A se le asocia una distribucion de
probabilidad de tipo normal, mientras que para aquellas tipo B se considera una distribucion de probabilidad
rectangular.

Para conocer la incertidumbre de la medicion de la manera mas precisa posible es importante identificar
las variables de influencia y sus errores asociados, con lo que se establece un presupuesto de incertidumbre,
consistente en un listado de los contribuyentes a la incertidumbre del proceso de medicion. Una vez identificadas
se procede a la normalizacion de cada una de las incertidumbres (ya que no se pueden mezclar incertidumbres
con distribuciones normales con aquellas con distribuciones rectangulares) y posteriormente se determina la
incertidumbre combinada, la que se obtiene haciendo:
v
uN
uX 2
uc (y) = t [ci u(xi )]
i=1

Donde u(xi ) es la incertidumbre normalizada o estandar de relacionada con la variable de influencia xi , y ci el

coeficiente de sensibilidad de la variable xi , el que esta dado por el efecto de la variable en la medicion y, es decir

y
ci =
xi
Para satisfacer las necesidades de algunas aplicaciones industriales, as como los requerimientos en areas
de la salud y seguridad, se obtiene una incertidumbre expandida up multiplicando la incertidumbre estandar
combinada uc por un factor de cobertura k. El proposito de obtener up es proveer de un intervalo alrededor del
resultado de una medicion en el que puede esperarse que se incluya una fraccion grande de la distribucion de
valores que pueden razonablemente ser atribuidos al mensurando[1] . La eleccion del factor k esta basada en la
probabilidad de cobertura o nivel de confianza requerido para el intervalo y sus posibles valores se listan en la
tabla 1 del anexo.

up (y) = kuc (y)

2.3. Red nacional de metrologa y laboratorios asociados

La Red Nacional de Metrologa (RNM) es una instancia reconocida por el Estado de Chile para articular
y administrar el sistema de aseguramiento metrologico, que garantiza las mediciones realizadas en Chile de
modo que sean comparables, trazables y aceptadas en otros pases. Esta conformada por una serie de institutos
designados, acreditados para la medicion y calibracion de patrones a laboratorios de calibracion, y proveen
material de referencia autentificando a los laboratorios de ensayos. Ademas, actuan como laboratorios pilotos
en los ensayos de aptitud (Proficiency Testing) que se organizan a nivel nacional[2] .

A continuacion se presentan los institutos designados y las variables trabajadas en cada uno de ellos:

CESMEC S.A.: masa, temperatura,

DICTUC S.A.: longitud

IDIC: fuerza
ENAER: presion, humedad
Universidad de Concepcion: magnitudes electricas

5
 CISA S.A.: flujo lquido
Codelco Division Chuquicamata: qumica (metales y sus aleaciones)
Instituto de Salud Publica: qumica (agua y alimentos)

CCHEN: radiaciones ionizantes

2.4. Laboratorio nacional de presion

El Instituto Designado en Presion se encuentra a cargo de la Empresa Nacional de Aeronautica de Chile,
ENAER, quienes desde el ano 1980 se dedican a la metrologa en esta area. En el ano 1999 se designo a ENAER
para implementar un laboratorio especializado en las mediciones de presion, proceso que culmino el ano 2003
con la acreditacion del laboratorio a cargo de esta implementacion.

ENAER, como administradora del Instituto Designado en Presion, tiene la responsabilidad de custodiar y
administrar los Patrones Nacionales de Presion segun consta en el Decreto Supremo N 76 promulgado el 21 de
Abril de 2003 por el Ministerio de Economa, Fomento y Turismo (antes Reconstruccion)[2] .

Actualmente este laboratorio puede emitir certificados de calibracion, segun lo establecido en la norma alema-
na DKD-5 y acreditado ISO/IEC 17025: 2005 por el Deutsche Akkreditierungsstelle (DAkkS) que corresponde
al cuerpo de acreditacion de Alemania.[3]

6
3. Estado del arte.

La variable presion se encuentra presente en variados campos teniendo aplicaciones en ciencia, tecnologa,
industria, medicina, etc. El constante progreso de estas areas ha estimulado el desarrollo de tecnicas mas sofisti-
cadas de medicion de este parametro, creandose as la necesidad de generar mejoras en la precision al momento
de medir, el establecimiento de patrones y la calibracion de instrumentos.

Entre los dispositivos utilizados como patron de medida en calibraciones de presion se encuentran manome-
tros (digitales o analogos), calibradores de presion y balanzas de peso muerto. Los primeros dos tipos de patrones
mencionados son instrumentos cuyo principio de funcionamiento se basa en la deformacion de una membrana o
una pieza elastica, siendo esta deformacion proporcional a la presion aplicada. Son utilizados para la calibracion
de manometros de baja exactitud, en rangos de presion que pueden ir desde -0.1 a 70[MPa], aproximadamente,
usando como medio de transmision de presion un gas, generalmente aire.

(a) Manometro digital. (b) Calibrador de presion. (c) Balanza de peso muerto.

Figura 1: Patrones de referencia para la variable presion.

Por otra parte las balanzas de peso muerto son utilizadas como patron en calibraciones de manometros de
mayor exactitud (del orden o menor a 0,02 % f.s.) y de otras balanzas de peso muerto, por medio del metodo de
flotacion simultanea o flotacion cruzada, cubriendo rangos de presion que van desde 1,5[kPa] y llegando incluso
a 500[MPa].

3.1. Principio de la balanza de presion

Una balanza de presion o balanza de peso muerto consiste en un piston vertical que rota libremente dentro
de un cilindro, donde ambos elementos, de buena factura y calidad mecanica, definen una superficie llamada
area efectiva [4] . En estos dispositivos el principio de funcionamiento se basa en la definicion de presion que se
presenta en la ecuacion 1, donde en particular la presion medida en la balanza es producida por una fuerza
vertical hacia arriba, realizada por efecto del fluido de transmision de presion actuando en la base del piston.

Existen balanzas neumaticas cuyo medio de transmision de presion puede ser nitrogeno o aire y cubren
rangos de medicion desde 1,5[kPa] a 70 [MPa], y balanzas hidraulicas para rangos de presion que van desde
0,1 [MPa] a 500[MPa] aproximadamente, utilizando como fluido de trasmision agua, alcohol o aceite (Sebacate,
Tellus 32 o de similares caractersticas).

F
Pe = (1)
A

7
 Figura 2: Esquema del sistema piston-cilindro.

Figura 3: Conjuntos de masas para balanza de pesos muertos.

La fuerza F resultante es compensada por el peso de una carga conocida puesta sobre el piston. Esta carga
consiste en un conjunto de masas previamente calibrado. Las masas tienen forma de disco y cada uno se identi-
fica con un ndice numeral, como las que se muestran en la figura 3; estas se apilan una sobre otra en la cabeza
del piston, donde distintas combinaciones de los discos generan distintos valores de presion. Luego, la presion
generada por la balanza es obtenida analizando cada una de las distintas componentes de las fuerzas que actuan
en el sistema y se tiene entonces que la ecuacion 1 toma la forma:

P   
aire
gl i mi 1 mi V (f aire ) + 2r
Pe = + (f aire )gl h (2)
A(0,t0 ) (1 + ( + )(t t0 ))(1 + Pn )

 
Donde gl es la aceleracion de gravedad en el sitio donde se lleva a cabo la medicion, i mi 1 aire
P
mi
es la
masa total de las cargas apiladas sobre el piston-cilindro considerando un factor de correccion por la flotabilidad
de las masas. El termino V gl (f aire ) representa el peso ejercido por el volumen libre o adicional del piston

8
sometido al empuje del fluido de transmision de la presion. V puede ser positivo si el piston tiene una forma
como las que se muestran en las figuras 4.a y 4.b, negativo si el piston acaba en un tope como se muestra en la
figura 4.c o nulo cuando es completamente cilndrico.

(a) (b) (c) (d)

Figura 4: Ejemplos de pistones con volumen adicional.

En la ecuacion 2 se tiene tambien la contribucion por la resistencia que ofrece el fluido de transmision dada
por la expresion 2r, donde r es el radio del piston-cilindro y la tension superficial del fluido.

En el denominador se tiene el area efectiva del sistema A(0,t0 ) y sus afecciones por temperatura y presion a
la que se somete. El factor ( + ) corresponde a los coeficientes de dilatacion termica de los materiales tanto
del piston y como del cilindro, mientras que (t t0 ) es la diferencia entre la temperatura del ensamble t y una
temperatura de referencia, tpicamente t0 = 20 C. Por otro lado, (1 + Pn ) representa el cambio en el area
efectiva por efecto de la presion ejercida Pn , con el coeficiente de deformacion con la presion o modulo de
compresibilidad asociado al ensamble piston-cilindro, informado por el fabricante.

Se considera tambien la contribucion a la presion por la diferencia de altura h que exista entre la balanza
patron y la balanza a calibrar, en el ultimo termino de la ecuacion 2, donde influyen ademas las densidades
tanto del aire aire como del fluido de transmision f . La densidad del aire se considera constante, con un valor
aproximado de 1, 2[kg/m3 ], mientras que la densidad del fluido es determinada segun la ecuacion 3a si se trata
de un medio lquido o bien segun la ecuacion 3b si se trata de un gas.

f = 0 (1 + Pn ) (3a)

Pn t 0
f = 0 (3b)
P0 t

Con 0 la densidad del fluido a temperatura y presion ambientales estandar, la compresibilidad del flui-
do, Pn es el valor nominal de la presion aplicada y P0 un valor de presion de referencia, en general P0 = 1[bar][5] .

Luego, considerando que se tiene un sistema cerrado y sin fugas, la presion en la balanza patron debe ser la
misma que en la balanza puesta a prueba, por lo que igualando y despejando se obtiene una expresion para el
area efectiva del sistema piston-cilindro de la balanza bajo calibracion, que depende de la presion como:

9
    
0 aire
V 0 (f aire ) + 2r0
P
gl i m i 1 m0
A0(Pe ,t0 ) =
i
(4)
Pe (1 + (0 + 0 )(t0 t0 ))

Donde las variables con apostrofo se relacionan con el instrumento bajo ensayo.

Luego sera posible determinar el area efectiva del ensamble a calibrar, conociendo los valores de masa, den-
sidad de cada masa, densidad del aire, gravedad local, volumen adicional, densidad del fluido de transmision,
radio del piston-cilindro, coeficiente de dilatacion termica, coeficiente de deformacion con la presion y tempera-
tura en el sistema, tanto del patron como del calibrando.

3.2. Montaje de las balanzas de presion y metodo de flotacion cruzada

La calibracion se lleva a cabo en un recinto controlado y acondicionado, donde la temperatura y humedad

del ambiente deben estar dentro de los margenes de utilizacion tanto de la balanza bajo calibracion, como del
patron. Idealmente la temperatura en el laboratorio se debe encontrar entre 18 C y 25 C, estabilizada en 1 C,
mientras que la humedad relativa del ambiente no debe exceder el 80 % [4] .

La superficie donde se situan las balanzas debe soportar el peso de ambas sin pandear y se debe comprobar
la verticalidad de ambos sistemas piston-cilindro por medio de un nivel de burbuja puesto sobre el piston.

Se comprueba el tiempo de rotacion y cada del piston de la balanza a calibrar para la maxima carga admi-
tida por el equipo. El periodo para ambos casos, en general, es informado por el fabricante, pero si no se cuenta
con este dato se debe confirmar que tanto la rotacion como la cada deben superar los 3 minutos, ya que este
tiempo posibilita llevar a cabo la toma de datos. Esta prueba permite tambien confirmar la verticalidad del
piston, ya que si este se encuentra inclinado el tiempo de rotacion y cada se ven disminuidos considerablemente.

Los equipos se conectan por medio de accesorios como valvulas, racores y tuberas, capaces de soportar las
presiones generadas y compatibles con el fluido de transmision a utilizar[6] . Es recomendable que el circuito de
conexion sea lo mas corto y con el menor numero de uniones posibles, procurando as una medicion mas estable,
reproducible y con menores posibilidades de fugas.

La diferencia de altura entre la balanzas montadas debe ser la mnima posible para as disminuir las contribu-
ciones a la presion y es registrada antes de comenzar la medicion. Para determinar correctamente esta diferencia
se utiliza la ecuacion 5, por lo que sera necesario conocer la altura de la punta de la balanza H, la que se mide
respecto de la mesa donde se encuentra montada, L el largo de cada piston y el recorrido que hace el piston
cuando se encuentra en operacion C, los ultimos dos datos son informados por el fabricante de cada equipo.
Si esta informacion no es suministrada o no se encuentra disponible, el usuario puede estimar la diferencia por
medio de una regla. Cabe mencionar que la diferencia de altura sera negativa cuando el instrumento bajo ensayo
se encuentre a mas altura que el patron y sera positiva en caso contrario.

h = ((Hp Lp ) + Cp ) ((Ht Lt ) + Ct ) (5)

La calibracion se realiza por el metodo de flotacion simultanea o flotacion cruzada. Este consiste en hacer

10
Figura 5: Esquema de montaje para balanzas de presion, con variables que determinan la diferencia de altura entre los instrumentos.

que el fluido empleado en la calibracion como medio de transmision de presion lleve a los dos conjuntos piston-
cilindro, el de la balanza a calibrar y el de la balanza patron, a su posicion de flotacion. El objetivo del metodo
radica en equilibrar la presion en ambas balanzas adicionando o quitando masa en el patron, mientras se mo-
nitorea la cada del piston. Las masas adicionales puestas sobre el patron para lograr el equilibrio entre ambas
balanzas son conocidas como Trim, en general consisten en pequenos pesos de acero inoxidable que suman algu-
nos gramos, con lo que es posible determinar la sensibilidad de la balanza, es decir el valor mnimo de masa para
el cual el equilibrio se rompe y los pesos en la balanza patron comienzan a descender desde su posicion de flota-
cion. Cuando se alcanza el equilibrio en el sistema los pistones se encuentran en flotacion en su nivel de referencia.

(b) Trim sobre campana de la balanza

(a) Conjunto de masas adicionales o Trim.

Figura 6: Masas adicionales

Cuando se encuentran en equilibrio la presion generada por el compresor y la presion debida al peso de las
masas sobre el area del piston, la friccion se reduce girando el piston. Si el conjunto de masas en flotacion se
mantiene girando en su posicion sin descender, en ambas balanzas, es posible medir. Se deberan registrar los
siguientes parametros:

11
 Temperatura en el ensamble piston-cilindro de la balanza patron
Temperatura en el ensamble piston-cilindro de la balanza bajo calibracion
Discos de masa y masas adicionales (o trim) utilizados para el valor de presion en medicion.
Presion ambiental
Temperatura ambiental

Figura 7: Sistema de balanzas de peso muerto conectadas para calibracion por metodo de flotacion cruzada.

En la calibracion se realizan al menos 3 series de medida con un mnimo de 5 puntos de presion[4] incluyendo
el 10 % del rango de medida como valor mnimo y llegando al fondo de escala de medicion de la balanza bajo
ensayo, mientras que el resto de los puntos se distribuyen en el intervalo de manera que los valores representen
de mejor manera la capacidad del equipo. Tpicamente se consideran el 10 %, 20 %, 40 %, 60 %, 80 % y 100 %
de la escala.

La toma de datos se lleva a cabo intercalando series de medida ascendentes y descendentes, cubriendo todo
el rango de medida.

3.3. Software y tratamiento de datos

El laboratorio custodio de patrones nacionales de presion cuenta con una serie de planillas disenadas en la
aplicacion Excel, de Microsoft Office, para calcular las cantidades de interes en calibraciones de balanzas de
peso muerto por metodo de flotacion simultanea: la presion en la balanza patron, area efectiva de la balanza
puesta a prueba e incertidumbres asociadas a estos parametros.

Las planillas de calculo se diferencian por numero de series de medida, numeros de puntos de presion medidos
y tipo de balanza en uso (hidraulica o neumatica), sin embargo son equivalentes al momento de la evaluacion
del modelo de presion, determinacion del area efectiva e incertidumbres asociadas. Cada archivo Excel cuenta
con al menos 5 hojas fundamentales para el tratamiento de los datos recogidos:

12
 (i) Mass: Esta seccion se dispone para la calibracion de los discos de masa asociados al instrumento bajo
ensayo.
(ii) Diff altura: Aqu se determina la diferencia de altura entre las balanzas.

(iii) Lambda: En esta hoja se determina el coeficiente asociado al piston del calibrando.
(iv) Pe y U: En este apartado se realiza la parte fundamental del tratamiento de datos, se evaluan las
ecuaciones 2 y 4, y se evalua la incertidumbre asociada a las cantidades calculadas.
(v) Certificado: Contiene el formato del documento que ha de ser despachado al cliente al termino del proceso
de calibracion.

Ademas, se tiene una planilla aparte donde se determina la combinacion de los discos de masa a utilizar en
cada balanza para los distintos puntos de presion a medir.

En primer lugar, conociendo el rango de medicion de la balanza a calibrar y habiendo definido los puntos de
presion a medir, se selecciona el patron mas adecuado a utilizar en el proceso de calibracion, segun rango de me-
dicion y fluido de transmision. Habiendo definido este, se ingresa en la planilla para determinar las combinaciones
de discos la informacion de masa, densidad de masa, area efectiva y el conjunto de masas asociadas al piston
patron. En esta planilla se han creado 10 rutinas de calculo automatico o macros de Excel que al ser ejecutadas
determinan la combinacion de discos de masas mas optima, entregando como resultado una cadena de texto con
los ndices de cada disco a usar precedidos por la sigla Grque representa la masa base del piston, que debe
considerarse en todos los casos, al igual que el numero 1que representa la masa de la campana de la balanza,
que corresponde a una estructura que se dispone sobre el piston para soportar el conjunto de masas apiladas,
como se muestra en el diagrama de la figura 8. Este tratamiento debe repetirse con el piston y masas de la balan-
za bajo ensayo. Finalmente, conociendo las aleaciones para cada balanza, es posible proceder a la toma de datos.

Figura 8: Diagrama de balanza con campana.

Habiendo realizado la medicion, el operador selecciona la planilla para el tratamiento de datos mas apropia-
da, segun el procedimiento realizado, e ingresa la informacion recogida. Los parametros de posicion y dimension
de los pistones se apuntan en la pagina Diff altura para determinar la diferencia de altura entre los equipos
del montaje; el valor resultante es trasladado a la hoja Pe y MU para ser utilizado en la evaluacion del modelo
descrito en la seccion 3.1. A continuacion, en la hoja Lambda, se debe indicar el modulo de compresibilidad,
coeficiente de Poisson y radio, tanto del piston como del cilindro perteneciente a la balanza en calibracion,

13
cantidades con las que se calcula el coeficiente 0 segun la ecuacion 6.

3c 1 (1 + c )rc2 + (1 0 )rp2
0 = + (6)
2Ep 2Ec (rc2 rp2 )

Por otro lado, los datos de temperatura en el patron, temperatura en el instrumento puesto a prueba, tempe-
ratura ambiental, presion ambiental y masas adicionales, son ordenados en la hoja Pe y MU del archivo Excel.
Se debe ingresar manualmente tambien la suma de masas para cada punto de presion y las caractersticas del
piston patron y del calibrando, esto es, rango de medicion, masa, densidad, largo, area efectiva (informada
por el cliente en el caso de la balanza bajo ensayo), volumen adicional, nivel de referencia, coeficiente de di-
latacion termica y coeficiente de deformacion con la presion; asimismo debe indicar marca, numero de serie y
conjunto de masas asociado. Es responsabilidad del usuario tambien ingresar las constantes asociadas al fluido
de transmision (densidad del medio, compresibilidad y tension superficial), la densidad del aire y gravedad local.

Si la informacion fue ingresada de manera correcta, es decir, en las celdas correspondientes y respetando las
unidades (temperatura en C, masa en Kg y presion en bar), la hoja Pe y U evaluara la ecuacion 2 y 4 para
cada punto de presion medido, de esta forma, si la toma de datos se llevo a cabo con 5 series de medida, cada
una con 6 puntos, la planilla entregara 30 resultados de presion en el patron y 30 resultados de area efectiva
para el piston-cilindro bajo estudio.

Por otra parte, en la pagina Pe y U se determina la incertidumbre asociada al area efectiva calculada, la
que se obtiene con la ecuacion 7:

v !2
u U (A(0,t ) ) 2 s(A0(Pe ,t0 ) )
u
0 0 0
U (A(Pe ,t0 ) ) = 2 A (Pe ,t0 ) t + (7)
2 A0 (Pe ,t0 )

Donde A0 (Pe ,t0 ) es el promedio de los valores de area efectiva entregados por la planilla tras la evaluacion
de la ecuacion 4 para cada punto de presion; U (A(0,t0 ) ) es la incertidumbre del area efectiva del piston-cilindro
patron informada tras su calibracion mas reciente, y s(A0(Pe ,t0 ) ) es la desviacion estandar de los valores de area
efectiva obtenidos en la calibracion.

Se estima tambien la incertidumbre asociada a la capacidad de medicion de presion de la balanza bajo

ensayo, la que se presenta como la ecuacion 8

q
U (Pe ) = 2 s(Pe )2 + s(A0(0,t0 ) )2 Pe2 + s(0 )2 Pe4 (8)

Donde s(Pe ), s(A0(0,t0 ) ) y s(0 ) corresponden a las desviaciones de presion, area efectiva y coeficiente de
compresibilidad, respectivamente, y son determinadas como:

v !2 !2
u
u mgl m0 g
s(Pe ) = t + 0 l (9a)
3A(0,t0 ) 3A(0,t0 )

14
 v !2
u U (A(0,t ) ) 2 (A0(Pe ,t0 ) )
u
0 0
s(A(0,t0 ) ) = t + (9b)
2 A0 (Pe ,t0 )

s 2 2
U (0 )

0 U ()
s( ) = + (9c)
2 2

Finalmente, los resultados obtenidos son trasladados a las secciones dispuestas en la hoja Certificado para
ser incluidos en el reporte que sera despachado, y as informar al cliente, dueno del instrumento estudiado, el
estado y condiciones de medicion de este.

15
4. Descripcion resumida del trabajo llevado a cabo
Para alcanzar los objetivos propuestos por la ENAER y en particular por el laboratorio nacional de pre-
sion, se diseno e implemento un sistema de planillas de calculo, en formato Excel, para el tratamiento de los
datos recogidos en calibraciones de balanzas de peso muerto por medio del metodo de flotacion cruzada. Para
ello se estudiaron los procedimientos utilizados hasta ese entonces en el laboratorio, detectando sus falencias y
proponiendo un metodo mas amigable, cuyo nucleo de calculo se encontrara conforme a la gua de calibracion
internacional vigente, tanto al momento de determinar los parametros de interes, como las incertidumbres aso-
ciadas a ellos.

El sistema desarrollado fue puesto a prueba con cantidades registradas en calibraciones pasadas realizadas
en el laboratorio, ejemplos de guas de calibracion y finalmente con la informacion recogida en una calibracion
real, arrojando resultados favorables en todos los casos.

Si bien el principal objetivo se relacionaba con las calibraciones de balanzas de peso muerto, debido a las
necesidades del laboratorio nacional, se desarrollaron dos sistemas adicionales para el procesamiento de datos
recogidos en calibraciones de distinta instrumentacion asociada a la variable presion, de acuerdo a las guas
de calibracion internacionales correspondientes, y realizando la respectiva comprobacion de su funcionalidad.
Finalmente se apoyo en la preparacion de material a presentar en la auditora del laboratorio nacional de hu-
medad, realizando un analisis estadstico de datos de calibracion de instrumentacion de humedad recogidos a
traves de los anos en el laboratorio.

16
5. Procedimiento de trabajo
Para llevar a cabo el trabajo requerido, se cumplieron 6 etapas que se detallan a continuacion

5.1. Investigacion y reconocimiento del area

El primer acercamiento al entorno de trabajo que se realiza en el laboratorio nacional de presion fue por
medio de la investigacion teorica a traves de la lectura de guas, protocolos, normas y publicaciones relacionadas.

En particular se estudio el libro Metrologa y sus aplicaciones con el fin de profundizar en el area de la me-
trologa conceptos, definiciones, instrumentos y procedimientos, que constituyen nociones fundamentales para
comprender la labor que se realiza al momento de calibrar instrumentacion de diversas areas. Se estudiaron
tambien los protocolos de calibracion validados y normas de calibracion regentes utilizados en el laboratorio
custodio de patrones nacionales de presion, destacando la Gua de calibracion de balanzas de la asociacion eu-
ropea de institutos nacionales de metrologa (Euramet por sus siglas en ingles) cg-3 version 1.0 (03/2011), el
Procedimiento ME-016 para la calibracion de balanzas de presion del Centro Espanol de Metrologa (CEM), la
Recomendacion internacional para distintos tipos de procedimientos de medicion y calibracion de la Organiza-
cion internacional de metrologa legal (OIML por sus siglas en frances), la Gua de usuario para mediciones de
presion de alta exactitud de AMETEK y el procedimiento para la calibracion de instrumentos de presion gauge
DKD-R 6-1, en sus diversas ediciones.

5.2. Exploracion de los procedimientos utilizados y deteccion de fuentes de error

Como se menciono con anterioridad, en la seccion 3.3 del presente reporte, el laboratorio nacional de presion
contaba con un conjunto de planillas disenadas en entorno Excel para el tratamiento de los datos recogidos en las
calibraciones, sistema que fue implementado en el laboratorio en el ano 2013 por el Organismo de acreditacion
Aleman (DAkkS) encargado de auditar las tareas del laboratorio nacional de presion.

Si bien las planillas satisfacen las necesidades de un tratamiento de datos rapido, entregando como resultados
la presion en la balanza patron, area efectiva de la balanza bajo ensayo e incertidumbres asociadas a ambas
cantidades, el procedimiento que realizan no se condice con la norma de calibracion de Euramet cg-3 version
1.0 con la que el laboratorio acredita trabajar, no se encuentran actualizados los modelos y ecuaciones para
la determinacion de la densidad del aire, fluido y area efectiva, establecidos por la norma internacional OIML
R-111-1 y la amplia interaccion que tiene el usuario con el programa en sus distintas secciones, ademas de la
dispersion de las areas a completar en una misma pagina del archivo, lleva a errores de calculo cuya fuente
resulta difcil de rastrear y, por tanto, de corregir rapidamente.

5.3. Diseno del nuevo sistema a implementar

Para comodidad del operario se trabaja en la misma plataforma para el diseno del nuevo sistema, es decir,
se conservara la funcionalidad que ofrece una planilla de calculo en la aplicacion Excel.

En primer lugar se busca disminuir la cantidad de archivos con el fin de desarrollar un metodo mas opti-
mizado, ordenado y por consecuencia mas amigable para el operario. De esta forma se hace necesario incluir la
opcion para distintos protocolos de calibracion, con la posibilidad de variar la cantidad de puntos de medida
y series de medicion. El nuevo entorno debera admitir calibraciones tanto con balanzas hidraulicas, como con
equipos neumaticos, realizando las modificaciones necesarias al momento de evaluar la ecuaciones.

Se debera corregir la forma en que las incertidumbres de los parametros son obtenidas, utilizando ahora el
metodo establecido en la gua para calibracion de balanzas de Euramet cg-3 version 1.0.

17
 Sera necesario tambien disminuir las zonas de interaccion entre la planilla y el usuario, limitando las areas
para el ingreso de datos, con el fin de tener mayor control de posibles errores en el procesamiento. De esta
manera, al restringir la interaccion usuario-planilla a una sola hoja en el archivo, se hace mas sencillo rastrear
malas practicas del operario que desencadenan fallas en los calculos.

Para hacer el procedimiento mas diligente se pretende reestructurar los macros que determinan las combi-
naciones de discos de masas, reduciendo las 10 rutinas existentes a 1, aumentando la eficiencia del proceso.

5.4. Desarrollo de la planilla

En primera instancia se actualiza el modelo para la determinacion de la densidad del fluido de transmision
f , la que se ve afectada por la presion ejercida y temperatura en el ensamble segun la ecuacion 10[2]

0 E
f = (10)
(E Pn )(1 + f (t t0 ))

Donde 0 representa la densidad del fluido bajo condiciones ambientales estandar, E es el modulo de com-
presibilidad, Pn es el valor nominal de la presion aplicada, f el coeficiente de dilatacion con la temperatura
del fluido y (t t0 ) es la diferencia entre la temperatura del ensamble y la temperatura de referencia t0 = 20 C.

Por otra parte, la densidad del aire circundante, aire , se ve afectada por factores ambientales segun la
relacion determinada por R.S. Davis[4] , pero para efectos de calibraciones de balanzas de peso muerto se utiliza
la aproximacion que se tiene en la ecuacion 11[5] donde se involucran la humedad relativa en el lugar H.R.,
temperatura ambiental tamb (en C) y presion del ambiente Pamb , esta ultima expresada en [mbar] o [hPa].

0, 34848Pamb 0, 00902 H.R. exp(0, 061 tamb )

aire = (11)
273, 15 + tamb

Ademas, se tiene que el area efectiva del piston-cilindro bajo calibracion se relaciona con la presion ejercida
en la balanza segun la ecuacion 12, donde A0(0,t0 ) es el area efectiva a presion nula y temperatura de referencia,
0 es el coeficiente de compresibilidad del ensamble y Pe es la presion aplicada.

A0(Pe ,t0 ) = A0(0,t0 ) (1 + 0 Pe ) (12)

As, teniendo los valores de presion Pe y de area A0(Pe ,t0 ) a partir de las ecuaciones 2 y 4 es posible obtener
el valor de area efectiva a presion nula y el coeficiente de compresibilidad del ensamble por medio de un ajuste
lineal a los datos mencionados.

Por otra parte se determinan los coeficientes de sensibilidad para el calculo de incertidumbre de acuerdo
a la gua de calibracion de EURAMET, para ello se evaluan las derivadas parciales respecto a cada variable
influyente en el modelo de presion y area efectiva. La gua de EURAMET informa tambien una distribucion de
probabilidad asociada al error informado por el fabricante o deducido a partir de la sensibilidad del instrumento
de medicion empleado para la determinacion de cada variable de interes. Estos datos son incluidos en el nuevo
sistema de computo.

18
 Cuadro 1: Coeficientes de sensibilidad

Variable Distribucion de probabilidad Coeficiente de sensibilidad

P  aire

gl i 1 m
Pe
mi Normal (k=2) mi = A(0,t0 ) [1+(+)(tt0 )](1+Pn )
i

 

mi 1 aire
P
i +V (f aire )
Pe m
gl Normal (k=3) gl = i
A(0,t0 ) [1+(+)(tt0 )](1+Pn ) + h(f aire )
 P mi 
i m +V
Pe
aire Normal (k=3) aire = gl i
A(0,t0 ) [1+(+)(tt0 )](1+Pn ) + h

P mi
gl aire i 2
Pe mi
mi Normal (k=2) mi = A(0,t0 ) [1+(+)(tt0 )](1+Pn )

Pe gl (f aire )
V Normal (k=2) V = A(0,t0 ) [1+(+)(tt0 )](1+Pn )

 
Pe V
f Normal (k=2) f = gl A(0,t0 ) [1+(+)(tt0 )](1+Pn ) + h

P   

gl i mi 1 aire +V (f aire ) +2r
Pe m
A(0,t0 ) Normal (k=2) A(0,t0 ) = A2(0,t
i
[1+(+)(tt0 )](1+Pn )
0)

 hP   i 
aire
gl i mi 1 m +V (f aire ) +2r Pn
Pe
Normal (k=2) = i
A(0,t0 ) [1+(+)(tt0 )](1+Pn )2

 hP   i 
aire
gl i mi 1 m +V (f aire ) +2r
Pe
Pn Normal (k=2) Pn = i
A(0,t0 ) [1+(+)(tt0 )](1+Pn )2

 hP   i 
aire
gl i mi 1 m +V (f aire ) +2r (tt0 )
Pe
+ Normal (k=2) (+) = i
A(0,t0 ) [1+(+)(tt0 )]2 (1+Pn )

 hP   i 
aire
gl i mi 1 m +V (f aire ) +2r (+)
Pe
t Normal (k=2) t = i
A(0,t0 ) [1+(+)(tt0 )]2 (1+Pn )

Pe
h Normal (k=2) h = gl (f aire )

Pe 2r
Normal (k=2) = A(0,t0 ) [1+(+)(tt0 )](1+Pn )

Pe
2r Normal (k=2) (2r) = A(0,t0 ) [1+(+)(tt0 )](1+Pn )

19
 Variable Distribucion de probabilidad Coeficiente de sensibilidad

F 0
P  aire

m0i Normal (k=2) m0i = gl i 1 0m
i

F 0
P m0i
m0i Normal (k=2) rhom0 = gl aire i rho2 0
i m
i

F 0
V0 Normal (k=2) V 0 = gl (f aire )
 
F 0
P m0i 0
aire Normal (k=2) aire = gl i m0 +V
i

F 0
f Normal (k=2) f = gl V 0

F 0
Normal (k=2) = 2r0

F 0
2r0 Normal (k=2) (2r 0 ) =

Estas formulaciones son consideradas al momento de definir los calculos en el entorno de procesamiento en
desarrollo, el que contara de 8 hojas:

(i) Masas
(ii) Medicion

(iii) Pe y A0
(iv) Incertidumbre
(v) Certificado (1 piston)

(vi) Certificado (2 pistones)

(vii) Comportamiento historico
(viii) Data

Para que el sistema que se quiere implementar cumpla con los objetivos establecidos en la fase de diseno se
crean 14 rutinas, en el editor de Visual Basic que incluye Excel, para agilizar el trabajo de calculo. Las modulos
definidos se listan a continuacion:

Modulo 1: funcion Calculo

Funcion que convierte los caracteres de una cadena de texto en operadores y/o referencias a celdas, lo que
permite evaluar funciones escritas en formato texto.

Modulo 2: macro Piston

El macro esta disenado para filtrar y trasladar las caractersticas del piston patron seleccionado para la
calibracion y su set de masas asociado, desde el listado de pistones y masas ubicados en la pagina Data,

20
a la hoja Pe y A0. Los datos son trasladados a rangos de celdas preestablecidos que son usados en el
calculo de la presion en la balanza patron, Pe .

Modulo 3: macro Indices Masas Patron

Este modulo contiene el macro encargado de determinar y escribir en la hoja Medicion los ndices que
hacen referencia a la posicion de las masas asociadas al patron listadas en la hoja Pe y A0, cuya combina-
cion de pesos sobre el piston seleccionado genera los distintos puntos de presion establecidos previamente
en relacion con el rango total del instrumento bajo calibracion. La rutina entrega cadenas de texto con una
secuencia numerica que representa la combinacion de pesos para cada punto de la calibracion. El texto
resultante consiste en una serie de numeros precedida por Gr; cada uno de los elementos de la serie se
encuentra separado uno de otro por un espacio simple. Si la secuencia de numeros (comenzando en 1 o
2) toma una forma ascendente, teniendose que la diferencia entre dos numeros contiguos es 1, y, ademas,
la secuencia supera 5 numeros consecutivos, el macro entrega la secuencia abreviada en forma de rango,
formado este por dos parentesis que encierran el numero inicial y final de la progresion, separados por dos
puntos verticales (ver figura 8).

Figura 9: Ejemplo de cadena de texto entregada por el modulo 3.

Modulo 4: macro Indices Masas Cal

De manera analoga al macro Indices Masas Patron, este modulo determina la combinacion de masas aso-
ciado a la balanza bajo calibracion que generan una presion aproximada a los valores nominales calculados
en base al rango maximo de medicion del instrumento, entregando una cadena de texto conservando el
formato antes explicado y ejemplificado en la figura 8.

Modulo 5: funcion MASAS

Esta funcion evalua la suma de pesos a partir de los ndices del conjunto de masas del patron, correspon-
dientes a cada punto de presion medido en la calibracion. La entrada de esta funcion es la cadena de texto
entregada por el macro Indices Masas Patron. En el texto, Gr se asocia a la masa del piston patron, el
primer ndice corresponde a la campana de la balanza y el resto de la serie se relaciona con los discos de
masa que son ordenados en la hoja Pe y A0 entre las celdas N2 y N30.

Modulo 6: funcion MASAS CALIBRACION

La funcion MASAS CALIBRACION tiene como entrada la cadena de texto que indica la combinacion
de masas que determina cada punto de presion en la balanza bajo ensayo y, de manera equivalente al
modulo 5, evalua y entrega la suma de los pesos, asignando al caracter Gr el valor de masa del piston, el
ndice 1 correspondera a la campana de la balanza y la serie restante seran los pesos ordenados en la hoja
Medicion entre las celdas U2 y U30.

Modulo 7: funcion MASAS INC

De igual forma que en la funcion MASAS su entrada es la cadena de texto de ndices de masa que carac-
teriza cada punto de presion. De manera analoga a la funcion MASAS, esta herramienta de calculo toma
los ndices y los evalua en su respectivo valor de incertidumbre. En este caso Gr corresponde al valor de
incertidumbre de la masa del piston, que se relaciona con la celda I31 de la hoja Pe y A0, mientras que

21
 1 sera la incertidumbre de la masa de la campana de la balanza de peso muerto patron y se encuentra en
la celda I1 de la misma hoja. El resto de los ndices hace referencia a las posiciones en la columna I2:I30
donde se encuentran ordenados los valores de incertidumbre con su respectivo valor de masa.

Modulo 8: funcion MASAS CALIBRACION INC

Funcion equivalente a MASAS INC para los valores de incertidumbre de las masas de la balanza a cali-
brar, las que se encuentran enlistados en la hoja Medicion en el rango P1:P31 y que son obtenidos de la
calibracion de las masas previamente hecha

Modulo 9: funcion M Ro

Evalua los ndices de combinacion de pesos del patron en el valor de la razon mmii correspondiente. Los
valores de esta razon se encuentran calculados en la columna BB entre las filas 1 y 31 de la hoja In-
certidumbre. Se conservan las relaciones de Gr y 1 con el piston y la campana de la balanza patron,
respectivamente.

Modulo 10: funcion M RoM2

Analogo a las funciones anteriores, tomando como entrada la serie de numeros que indica las masas que

generan la presion en cada punto, evalua dichos indicadores en la razon mm2i , cuyo calculo, para cada valor
i
de densidad de masa y masa, se encuentra en la columna BC de la hoja Incertidumbre entre las mismas
filas antes mencionadas, 1 y 31.

Modulo 11: funcion RoAM

Al igual que las funciones anteriores
 usa
 como entrada los ndices de combinacion de masas para evaluar
en cada punto la relacion 1 aire
mi
necesaria para el calculo de incertidumbre de la presion patron.
Estos valores se encuentran enlistados en la hoja Incertidumbre en el rango que va de BD1 a BD31,
conservando la relacion de posicion que se establecio para las masas con Gr asociado al piston y 1 a la
campana.

Los modulos 12, 13 y 14 (M Ro Cal, M RoM2 Cal, RoAM Cal respectivamente) son equivalentes a los modu-
los 9, 10 y 11, pero en este caso los valores considerados para los calculos a realizar son los asociados a la balanza
bajo ensayo, cuya estimacion se encuentra ordenada en la hoja Incertidumbre entre las celdas BH1 y BK31.

5.5. Comprobacion del programa

La planilla es puesta a prueba para confirmar la eficiencia de calculo que se buscaba al definir las rutinas
creadas, la precision en los resultados y su coherencia con el metodo de calibracion de EURAMET. Para ello se
considera el ejemplo de calibracion presente en el procedimiento ME-016 de CEM, ya que en este documento se
incluyen tanto los datos de entrada como los resultados tras el tratamiento de datos entregando los valores de
presion y area, ademas de los errores relativos a cada parametro, obtenidos de acuerdo al metodo propuesto por
el Centro Espanol de Mediciones, el que consiste en una aproximacion al procedimiento de EURAMET. Para
hacer mas fiel la comparacion, en el archivo Excel, se modifica la hoja asociada a la incertidumbre agregando
secciones donde se realizan los calculos de acuerdo a la iniciativa de CEM. Hecho esto se ingresan los datos del
ejemplo en la seccion Medicion de la planilla de calculo.

El nuevo sistema se evalua tambien con datos tomados en calibraciones anteriores realizadas en el laboratorio
nacional de presion y se comparan los resultados conseguidos con los informados en el perodo en que fue hecha
la calibracion y con los informados por el fabricante.

22
 En las comprobaciones llevadas a cabo se observo un desempeno positivo del sistema desarrollado. En el
caso del ejemplo de CEM los valores resultantes diferan de los informados en el documento entre el sexto y
septimo decimal para los resultados de presion y area, mientras que en los errores la diferencia se perciba entre
la decima y undecima cifra decimal. Por otra parte con los datos de calibraciones anteriores se obtuvo mayor
diferencia entre los resultados, lo que se condice con las modificaciones realizadas al modelo, sin embargo el
resultado entregado por la planilla se aproximaba mas al valor informado por el fabricante. Finalmente en la
verificacion realizada con los datos de entrada de la comparacion internacional realizada por el laboratorio en el
ano 2011, en un rango de 7[MPa] se obtuvieron valores mayores en presion que se aproximaban mas al patron
designado en la medicion, mientras que la incertidumbre se mantuvo en el rango.

5.6. Mejoras y desarrollo de planillas para calibraciones de instrumentos de menor

exactitud
Se complementa el entorno de calculo incluyendo en la seccion Medicion una columna para indicar la
sensibilidad de masa de la balanza patron y a partir de este valor, que debe ser ingresado por el operario, se
determina el equivalente en presion de esta sensibilidad. Se informa tambien la contribucion a la presion por
la diferencia de altura entre los instrumentos y la contribucion a la incertidumbre por el error asociado a la
presion ambiental. Estos valores son entregados para que el usuario tenga una nocion inicial de los aportes a la
incertidumbre por las condiciones de montaje y ambientales.

Se establece ademas que la masa de la campana se asocie ahora al ndice 0 de la secuencia numeral, de esta
forma las masas se ordenan desde 1 a 30 y se evitan mayores confusiones en el listado de masas. Los macros y
funciones relacionados son actualizados para que reconozcan esta asociacion. Tras las modificaciones la planilla
es verificada nuevamente, de manera tal de asegurar una optima ejecucion de sus funciones.

Habiendo comprobado la planilla tras las modificaciones, se procede a la implementacion de archivos de

calculo adicionales para ser utilizados en la calibracion de elementos de menor exactitud como son calibradores
y manometros, para ello se toma como gua de calculo el documento DKD-R 6-1, version 2014, que contiene el
protocolo para calibracion de estos elementos. Estos instrumentos pueden ser calibrados utilizando como patron
una balanza de peso muerto, en el caso de mejor precision de medida (entre 0.015 % y 0.05 % de la lectura),
o bien un manometro de mejor exactitud respecto del instrumento bajo ensayo. De esta forma se disenan dos
archivos de Excel para el calculo de la presion y su incertidumbre en calibraciones de manometros donde el
patron sea una balanza de peso muerto y para aquellos montajes donde el patron sea otro manometro.

Figura 10: Montajes utilizados para la calibracion de otros elementos de medicion de presion como calibradores, manometros o
transductores.

23
 Para el sistema balanza-manometro se conservan 2 macros y 6 funciones de los disenados en la planilla
balanza-balanza, ajustando los parametros necesarios. En este caso se conserva el modelo para calcular la pre-
sion de la ecuacion 2 y se determina su error por medio de los coeficientes de sensibilidad de las variables de
influencia, sin embargo se considera tambien la contribucion por repetibilidad de la medicion, histeresis, resolu-
cion del dispositivo bajo calibracion, desviacion de cero y estabilidad de la medicion.

Por otra parte para el sistema de calibracion manometro-manometro se disena una nueva planilla basandose
en el protocolo detallado en la gua de calibracion DKD-R 6-1 version 2014. En este caso la presion es obtenida
directamente de la lectura de los instrumentos, por lo que ya no se hace necesario su calculo, pero igualmente se
debe hacer una correccion de los mensurando por su desviacion de cero y presion hidrostatica por la diferencia
de altura entre los dispositivos.

La incertidumbre en este caso esta dada por la suma de los cuadrados de los errores relacionados a la diferen-
cia de altura, repetibilidad, histeresis, desviacion de cero, resolucion del calibrando e incertidumbre certificada
del patron.

24
6. Resultados
Los resultados de la validacion realizada con el ejemplo de la gua CEM se muestran en la tabla XX, don-
de los calculos fueron realizados tambien con el software MATLAB para confirmar el buen desarrollo del sistema.

Cuadro 2: Resultados de validacion con ejemplo de gua de calibracion de CEM

Ejemplo CEM Planilla de calculo MATLAB

1002031 1002031,29750872 1002031,2975088

2502043 2502042,80317243 2502042,8031725

4002053 4002053,20647842 4002053,2064785

5002052 5002052,16396107 5002052,1639611

6002027 6002026,70266258 6002026,7205467

6001984 6001984,00394180 6001984,0218258

5002036 5002036,40954544 5002036,4095455

4002030 4002030,15841109 4002030,1584112

2502070 2502069,51892379 2502069,5189239

1002045 1002044,53723134 1002044,5372314

Mientras que las incertidumbres para el primer punto de presion fueron:

Cuadro 3: Resultados de incertidumbres en validacion con ejemplo de CEM

Incertidumbre de la presion en el patron [Pa]
Ejemplo CEM 46
Planilla de calculo 45,53370437
MATLAB 45,533704367
Incertidumbre de la fuerza [N]
Ejemplo CEM 1,1 103
Planilla de calculo 1,11218569 103
MATLAB 1,112185689 103
Incertidumbre del area efectiva [m2 ]
Ejemplo CEM 3,92 109
Planilla de calculo 3,91566607 109
MATLAB 3,915665999 109

Por otra parte los resultados del ajuste lineal para el ejemplo de CEM fueron:

25
 Cuadro 4: Resultados ajuste lineal en validacion con ejemplo de CEM
Intercepto con eje vertical A0(0,t0 ) [m2 ]
Ejemplo CEM 8,06435 105
Planilla de calculo 8,0643516 105
MATLAB 8,064350797 105
Pendiente de la recta 0 A0(0,t0 ) [m2 P a1 ]
Ejemplo CEM 2,9 1016
Planilla de calculo 2,858184 1016
MATLAB 2,8581842625 1016
Modulo de compresibilidad 0 [P a1 ]
Ejemplo CEM 3,57 1012
Planilla de calculo 3,57238579 1012
MATLAB 3,5723857898 1012

Los resultados obtenidos por la planilla se diferencian con los informados en el ejemplo en menos de un 1 %
en todas las cantidades, por lo que se consideran validos, mientras que al comparar los valores arrojados por
la planilla con aquellos conseguidos con Matlab se tiene una congruencia en la mayora de las cifras decimales,
por lo que se concluye que los procedimientos que se realizan en la planilla son correctos y sus resultados son
aptos para ser utilizados en calibraciones de balanzas de peso muerto.

Por otra parte, en la grafica de la figura 11 se tienen los resultados de la evaluacion realizada a la planilla con
los datos de la comparacion internacional donde el laboratorio nacional de presion participo en el ano 2011,
para un rango de medicion de 0-7[MPa].

Figura 11: Grafica comparativa.

En el primer intervalo, se observa que los resultados obtenidos por la planilla (color azul en el grafico) son
mayores a los informados en el ano 2011 por el laboratorio nacional (triangulos de color rojo en el grafico) y
mas cercanos al patron (crculos negros en el grafico), mientras que las incertidumbres se mantienen sin grandes
variaciones, interceptando los datos del patron.

26
7. Conclusiones
Se creo e implemento un sistema de computo en la aplicacion Excel para ser utilizado en calibraciones de
balanzas de presion por medio del metodo de flotacion simultanea o cruzada. El programa calcula los valores de
presion en el instrumento patron, area efectiva e incertidumbres asociadas a cada cantidad. El valor de presion
es corregido al verse afectado por las variables de influencia ambientales y mecanicas, las que se relacionan
segun la ecuacion 2. Por otra parte, el area efectiva del piston-cilindro bajo calibracion es entregado siguiendo
dos formas de calculo, por medio de la ecuacion 4 y por medio de una ajuste lineal a los datos segun la relacion
que se muestra en la ecuacion 12, de donde se obtiene el valor de area efectiva a presion nula. El sistema es
puesto a prueba mostrando un buen desempeno y entregando cifras validas para ser utilizadas en calibraciones
de balanzas de peso muerto.

De los resultados conseguidos por este medio de calculo se observa que el area efectiva de un piston-cilindro
de acero inoxidable, empleado en calibraciones con balanzas neumaticas e hidraulicas, se puede considerar no
afecta por la presion ejercida sobre esta, ya que la diferencia porcentual entre los resultados obtenidos por
el metodo del promedio y de regresion lineal difieren en un margen muy inferior al 1 %, en particular en las
validaciones realizadas se obtuvieron diferencias de 0.0013 % para el ejemplo de CEM, 0.0005 % utilizando
los datos de la comparacion internacional del ano 2011 en rango 7[MPa] y 0.0001 % con informacion de la
comparacion internacional del ano 2011 en rango 120[kPa].

Es importante mencionar tambien, que esta experiencia me permitio conocer un entorno laboral antes desco-
nocido, aproximandome a la metrologa desde un punto de vista mas practico, conociendo normas, procesos,
instrumentacion y las relaciones entre laboratorios nacionales de la red de metrologa e institutos internacionales
que hacen posible el desarrollo de esta disciplina. Cabe destacar la labor del laboratorio custodio de presion
con los laboratorios de calibracion del pas, prestando asistencia y orientacion en intercomparaciones nacionales
de la variable, ademas de encontrarse en constante mejoramiento de sus procesos, lo que se ve reflejado en los
buenos resultados obtenidos en comparaciones internacionales realizadas en el ultimo tiempo.

27
8. Referencias y Bibliografa
[1] Adolfo Escamilla Esquivel, Metrologa y sus aplicaciones
[2] https://www.metrologia.cl
[3] https://www.enaer.cl

[4] EURAMET, Calibration of Pressure Balances, version 1.0 (03/2011)

[5] Physikalisch-Technische Bundesanstalt DKD, Guideline DKD-R 6-1, Calibration of Pressure Gauges
[6] Centro espanol de metrologa, Procedimiento ME-016 para la calibracion de balanzas de presion.
[7] Scott S. Bair, Pressure Rheology for Quantitative Elastohydrodynamics

[8] R. S. Davis, Equation for the Determination of the Density of moist air
[9] Organisation Internationale de Metrologie Legale, International organization of legal metrology

28
9. Anexo
Tabla para determinar coeficiente de expansion.

29