Higrómetro Digital para Altitudes de 3000m
Higrómetro Digital para Altitudes de 3000m
NACIONAL
DIRECTOR:
ING. ORTEGA LÓPEZ MIGUEL DARÍO, [Link].
[Link]@[Link]
CODIRECTOR:
ING. ZÚÑIGA PUEBLA HUGO FRANCISCO, [Link].
[Link]@[Link]
Certifico que el presente trabajo fue desarrollado por LUCIANO MARTÍN ALBUJA
JARAMILLO y CARLOS ALEXIS ALMEIDA ANSASOY, bajo mi supervisión.
_______________________ ________________________
Ing. Miguel Ortega, [Link]. Ing. Hugo Zúñiga, [Link].
DIRECTOR DE PROYECTO CODIRECTOR DE PROYECTO
i
DECLARACIÓN
__________________________ ____________________________
Luciano Martín Albuja Jaramillo Carlos Alexis Almeida Ansasoy
ii
DEDICATORIA
Luciano Albuja
Este trabajo lo dedico a mi familia que siempre está ahí cuando necesito un consejo y
apoyo.
Carlos Almeida
A mi madre Gladys, ya que con todo su apoyo y cariño he llegado a culminar esta etapa
tan importante de mi vida.
A mi papito Carlos, que desde el cielo sé que me apoya y me da fuerzas para luchar cada
día.
A mis hermanos, Marlon y Jaider, que me han inspirado a ser mejor persona.
A mi hermano mayor Roberth, por ser la persona que me ha ayudado a cumplir este sueño
tan grande.
A mi hijita Valentina, que ha sido mi fuerza y mi voluntad para lograr todo lo que me
proponga con el fin de darle un mejor futuro.
A mi mujercita hermosa, Mayra, por ser la persona que ha compartido en toda mi etapa
universitaria, me ha apoyado y brindado todo su amor.
A mis amigos y familiares que han sido participes con su apoyo y cariño en toda mi vida.
“Todos nuestros sueños pueden hacerse realidad si sólo tenemos el coraje de perseguirlos”
Walt Disney
iii
AGRADECIMIENTO
Luciano Albuja
iv
AGRADECIMIENTO
A Dios por haberme otorgado la vida y la maravillosa familia que tengo, por guiarme
cada día y brindarme la fuerza para salir siempre adelante.
A mi madre Gladys, ya que con sus valores y consejos he llegado a ser la gran persona
que soy hoy en día. Por estar a mi lado en todos los momentos buenos y malos que hemos
tenido, y seguir luchando por mantener siempre nuestra familia unida. Deseo que Diosito
siempre le vida para que me vea cumplir muchas metas más.
A mi hermano Roberth, por ser, más que un hermano mayor, un padre que me ha
brindado todo su apoyo para cumplir esta meta tan importante en mi vida. Sabes que te
quiero mucho y gracias por todas tus palabras de apoyo.
A mi hija Valentina, que a pesar de que llego a mi vida cuando estaba muy joven, me
ha enseñado a luchar cada día por ser un gran padre. Por brindarme todo su amor y
hacerme sonreír cuando los días están más oscuros. Todo lo hago por ti, y para verte feliz
toda la vida. Te amo mi pedacito de cielo.
A mi mujercita Mayra, por estar a mi lado en todo momento, por darme una hermosa
hija y un gran amor. Gracias por apoyarme y ayudarme a culminar este trabajo. Te amo
mucho.
A mis amigos de universidad, con los cuales he compartido buenos y malos momentos,
pero siempre han estado ahí para brindarme todo su cariño. Con ustedes he formado una
pequeña familia lejos de casa.
A todos mis familiares que han sido participes en el transcurso de mi vida. Los quiero
mucho.
Carlos Almeida
v
ÍNDICE
CERTIFICACIÓN……………………………………………………………….….i
DECLARACIÓN………………………………………………...………………...ii
DEDICATORIA ……………………...…………………………………………...iii
AGRADECIMIENTO …………………………………………………………….iv
AGRADECIMIENTO…………………………………………………………..……….v
RESUMEN ....................................................................................... xiv
ABSTRACT....................................................................................... xv
INTRODUCCIÓN................................................................................. 1
1. MARCO TEÓRICO .......................................................................... 3
1.1. Psicrometría ........................................................................................................... 3
1.1.1. Cartas psicrométricas ............................................................................................. 3
1.1.2. Humedad relativa (hr) ............................................................................................. 4
1.1.3. Temperatura de bulbo seco (tbs) ............................................................................ 5
1.1.4. Temperatura de bulbo húmedo (tbh) ...................................................................... 5
1.1.5. Temperatura de punto de rocío .............................................................................. 6
1.2. Higrómetro ............................................................................................................. 7
1.2.1. Higrómetros mecánicos (por deformación) ............................................................. 7
1.2.2. Higrómetros de bulbo seco y bulbo húmedo ........................................................... 8
1.2.3. Higrómetros por condensación ............................................................................. 10
1.2.4. Higrómetros de absorción infrarroja ...................................................................... 10
1.2.5. Higrómetros de sales higroscópicas ..................................................................... 11
1.2.6. Higrómetro digital ................................................................................................. 12
1.3. Arduino y sensores ............................................................................................... 13
1.3.1. ¿Qué es Arduino? ................................................................................................ 13
1.3.2. Software ............................................................................................................... 14
1.3.3. Tipos de placas Arduino ....................................................................................... 14
1.3.4. Sensores de humedad y temperatura ................................................................... 15
1.3.5. Sensores barométricos......................................................................................... 17
1.4. Definiciones de error ............................................................................................ 19
1.5. Ajuste de curvas ................................................................................................... 20
1.5.1. Método numérico directo para ajuste de curvas ................................................... 20
1.5.2. Herramientas matemáticas para ajuste de curvas ................................................ 21
2. METODOLOGÍA ............................................................................ 24
vi
2.1. Verificación del estado de funcionamiento de los termómetros del psicrómetro de
flujo continuo (Ps FC) ...................................................................................................... 24
2.2. Verificación del comportamiento de varios higrómetros digitales comerciales a
diferentes altitudes........................................................................................................... 25
2.3. Verificación del comportamiento del prototipo de higrómetro digital sin
acondicionamiento de señal............................................................................................. 26
2.3.1. Análisis de selección de materiales electrónicos .................................................. 26
2.3.2. Materiales del higrómetro digital sin acondicionamiento de señal ......................... 27
2.3.3. Prototipo electrónico de higrómetro digital sin corrección ..................................... 31
2.3.4. Análisis de las propiedades psicrométricas del aire con el prototipo de higrómetro
digital sin acondicionamiento de señal ............................................................................. 36
2.4. Análisis de alternativas para la corrección de altitud (A) y humedad relativa (HR) del
prototipo sin corrección .................................................................................................... 37
2.4.1. Acondicionamiento y corrección de altitud ............................................................ 37
2.4.2. Acondicionamiento y corrección de valores de humedad relativa ......................... 40
2.4.3. Análisis de sensibilidad de la humedad relativa respecto a la temperatura de bulbo
húmedo............................................................................................................................ 42
2.4.4. Ecuación de corrección para los datos preliminares ............................................. 45
2.5. Implementación de las alternativas de corrección y cálculos de psicrometría en el
prototipo........................................................................................................................... 52
2.5.1. Código de funcionamiento del higrómetro digital con corrección .......................... 53
2.6. Verificación del comportamiento del prototipo con acondicionamiento de señal de 0
a 3040 msnm ................................................................................................................... 56
2.7. Higrómetro digital con acondicionamiento ............................................................ 56
2.7.1. Materiales y componentes .................................................................................... 57
2.7.2. Características físicas y técnicas del higrómetro digital corregido ........................ 63
2.7.3. Uso del higrómetro digital fabricado ..................................................................... 67
2.7.4. Costo de fabricación del higrómetro digital con corrección ................................... 69
3. RESULTADOS Y DISCUSIÓN ...................................................... 71
3.1. Resultados ........................................................................................................... 71
3.1.1. Programa para el cálculo de humedad relativa ..................................................... 71
3.1.2. Pruebas de higrómetros digitales comerciales a diferentes altitudes .................... 71
3.1.3. Pruebas del higrómetro digital sin acondicionamiento de señal ............................ 77
3.1.4. Pruebas del higrómetro digital con acondicionamiento de señal ........................... 85
3.2. Discusión.............................................................................................................. 88
3.2.1. Datos obtenidos con el higrómetro digital sin corrección ...................................... 88
3.2.2. Datos de humedad relativa del aire analizadas en Tulcán (3040 msnm) .............. 89
3.2.3. Selección de la ecuación adecuada para la corrección de temperatura de bulbo
húmedo............................................................................................................................ 91
vii
3.2.4. Datos obtenidos con el higrómetro digital con corrección ..................................... 92
4. CONCLUSIONES Y RECOMENDACIONES ................................. 93
4.1. Conclusiones ........................................................................................................ 93
4.2. Recomendaciones ................................................................................................ 94
BIBLIOGRAFÍA ................................................................................ 95
ANEXOS ........................................................................................... 97
viii
CONTENIDO DE FIGURAS
Capítulo 1
Figura 1.1. Carta psicrométrica a temperaturas normales y presión barométrica de 101.325
kPa (al nivel del mar). ....................................................................................................... 3
Figura 1.2. Líneas de humedad relativa. ........................................................................... 5
Figura 1.3. Líneas de temperatura de bulbo seco. ............................................................ 5
Figura 1.4. Líneas de temperatura de bulbo húmedo. ....................................................... 6
Figura 1.5. Líneas de temperatura de punto de rocío........................................................ 6
Figura 1.6. Esquema de funcionamiento de higrómetros mecánicos. ............................... 7
Figura 1.7. Higrómetro de cabello humano, FIQ-EPN. ...................................................... 8
Figura 1.8. Esquema de funcionamiento de higrómetros de bulbo seco y bulbo húmedo. 9
Figura 1.9. Higrómetro de bulbo seco y bulbo húmedo (psicrómetro de flujo continuo) de la
FIM-EPN. .......................................................................................................................... 9
Figura 1.10. Dispositivo empleado para encontrar el punto de rocío. ...............................10
Figura 1.11. Higrómetro de absorción infrarroja. ..............................................................11
Figura 1.12. Esquema de funcionamiento de un higrómetro de sales higroscópicas........12
Figura 1.13. Higrómetros digitales del laboratorio de termodinámica de la FIM-EPN. ......13
Figura 1.14. Placa Arduino uno básica.............................................................................13
Figura 1.15. Software utilizado por Arduino. ....................................................................14
Figura 1.16. Tipos de placas Arduino...............................................................................15
Figura 1.17. Sensor DHT11 y sus pines de conexión.......................................................16
Figura 1.18. Sensor DHT22 y sus pines de conexión.......................................................17
Figura 1.19. Sensor BMP280 para Arduino. .....................................................................18
Figura 1.20. Sensor BMP180 para Arduino. .....................................................................19
Figura 1.21. Ajuste de curvas, regresión por mínimos cuadrados. ...................................20
Figura 1.22. Ejemplos de interpolación polinomial: a) de primer grado (lineal) que une dos
puntos, b) de segundo grado (cuadrática o parabólica) que une tres puntos y c) de tercer
grado (cúbica) que une cuatro puntos. .............................................................................21
Figura 1.23. Funciones de Excel interconstruidas relacionadas con el ajuste de datos por
regresión..........................................................................................................................22
Figura 1.24. Ajuste de un modelo logarítmico. .................................................................22
Figura 1.25. Algunas funciones de Matlab para implementar interpolación, regresión, ....23
Figura 1.26. Ejemplo de ajuste polinomial realizado en Matlab. a. Se observa los valores
interpolados linealmente. b. Ajuste polinomial de quito grado. .........................................23
Capítulo 2
Figura 2.1. Metodología del proyecto. ..............................................................................24
Figura 2.2. Higrómetros digitales utilizados en la toma de datos de este proyecto...........25
Figura 2.3. Sensor AM2302, estructura física. .................................................................28
ix
Figura 2.4. Sensor BMP180 estado físico. .......................................................................28
Figura 2.5. Pantalla LCD 1602 forma física. ....................................................................29
Figura 2.6. Arduino uno forma física. ...............................................................................29
Figura 2.7. Agua desmineralizada....................................................................................30
Figura 2.8. Ventilador 8x8. ...............................................................................................30
Figura 2.9. Cartón pluma. ................................................................................................31
Figura 2.10. Dimensiones de la caja fabricada de cartón pluma. .....................................32
Figura 2.11. Caja de cartón pluma montada en psicrómetro de flujo continuo. ................32
Figura 2.12. Ubicación de sensores para el prototipo inicial.............................................33
Figura 2.13. Conexión de sensor DTH22 en placa Arduino uno. ......................................33
Figura 2.14. Conexión de sensor BMP180 en placa Arduino uno. ...................................34
Figura 2.15. Conexiones para el prototipo de higrómetro digital. .....................................34
Figura 2.16. Código para definir librería e inicio de encendido. ........................................35
Figura 2.17. Código de lectura de sensores. ....................................................................36
Figura 2.18. Gráfica de altitudes de referencia y experimentales, ecuación de corrección en
altitud. ..............................................................................................................................39
Figura 2.19. Curva de humedad relativa real vs humedad relativa experimental (higrómetro
fabricado).........................................................................................................................41
Figura 2.20. Programa de cálculo de Th. .........................................................................41
Figura 2.21. Curva de temperatura de bulbo húmedo real vs temperatura de bulbo húmedo
experimental (cálculo) ......................................................................................................42
Figura 2.22. Carta psicrométrica para una altitud de 2800 msnm a temperatura normal. .43
Figura 2.23. Ejemplo de cálculo de valores de humedad relativa variando la temperatura
de bulbo húmedo. ............................................................................................................44
Figura 2.24. Curvas de sensibilidad de humedad relativa por cada grado Celsius de
temperatura de bulbo húmedo (%/°C). .............................................................................44
Figura 2.25. Gráfica de temperatura de bulbo húmedo real vs experimental, ecuación de
corrección Excel. .............................................................................................................45
Figura 2.26. Gráfico de ecuación de corrección de temperatura de bulbo húmedo utilizando
Matlab. .............................................................................................................................51
Figura 2.27. Gráfica de temperatura de bulbo húmedo real vs experimental, ecuación de
corrección Matlab. ...........................................................................................................52
Figura 2.28. Código para la corrección y lectura de altitud. ..............................................53
Figura 2.29. Código del valor de referencia de la temperatura de bulbo húmedo.............53
Figura 2.30. Código para el cálculo de la temperatura de bulbo húmedo sin la corrección.
........................................................................................................................................54
Figura 2.31. Código para corrección de la temperatura de bulbo húmedo y cálculo de la
humedad relativa utilizando la temperatura corregida. .....................................................54
Figura 2.32. Código para plasmar los valores corregidos y no corregidos en el higrometro
fabricado. .........................................................................................................................55
x
Figura 2.33. Códigos mostrados en la pantalla lcd del higrómetro fabricado. ...................55
Figura 2.34. Placa Arduino uno utilizada en el higrómetro digital con corrección. ............57
Figura 2.35. a) Sensor DHT22 de humedad y temperatura b) Sensor barométrico BMP180.
........................................................................................................................................58
Figura 2.36. Placas o shields de expansión. ....................................................................58
Figura 2.37. Fuente de alimentación para Arduino uno, batería de 9 voltios. ...................59
Figura 2.38. Modelo 3D de la carcasa del higrómetro digital con corrección. ...................60
Figura 2.39. Pantalla LCD Nokia 5110 de Arduino. ..........................................................62
Figura 2.40. Código elaborado para el encendido de la pantalla Nokia 5110. ..................63
Figura 2.41. Partes del higrómetro digital.........................................................................63
Figura 2.42. Propiedades mostradas en la pantalla Nokia 5110 del higrómetro digital. ....64
Figura 2.43. Parte posterior del higrómetro digital con corrección. ...................................65
Figura 2.44. Parte superior del higrómetro digital con corrección. ....................................65
Figura 2.45. Diagrama de flujo del funcionamiento del higrómetro digital fabricado………68
Figura 2.46. Metodología desarrollada en el proyecto………………………………………. 70
Capítulo 3
Figura 3.1. Hoja de presentación programa para el cálculo de humedad relativa. ..........71
Figura 3.2. Datos de humedad relativa de los higrómetros digitales (H1, H2 y H3), y
psicrómetro de flujo continuo (2806msnm valores reales). ..............................................73
Figura 3.3. Datos de humedad relativa del higrómetro digital H4 y del psicrómetro de flujo
continuo (2806msnm valores reales). ..............................................................................74
Figura 3.4. Datos de humedad relativa de los higrómetros digitales (H5 y H6), y psicrómetro
de flujo continuo (2806msnm valores reales). ..................................................................76
Figura 3.5. Datos de humedad relativa de higrómetros digitales (H1, H2, H3, H4), y las
cartas psicrométricas (72 msnm valores reales). .............................................................76
Figura 3.6. Datos de humedad relativa real y experimental sin corrección en Atacames. 78
Figura 3.7. Datos de humedad relativa real y experimental sin acondicionamiento de señal
en Santo Domingo. ..........................................................................................................79
Figura 3.8. Datos de humedad relativa real y experimental sin acondicionamiento en
Rocafuerte. ......................................................................................................................80
Figura 3.9. Datos de humedad relativa real y experimental sin corrección en Quito Centro.
........................................................................................................................................82
Figura 3.10. Datos de humedad relativa real y experimental sin corrección en Quito Sur.
........................................................................................................................................83
Figura 3.11. Datos de humedad relativa real y experimental sin corrección en Tulcán. ...84
Figura 3.12. Datos de humedad relativa real y experimental con corrección en Quito Centro.
........................................................................................................................................86
Figura 3.13. Datos de humedad relativa real y experimental con corrección en Tulcán. ..87
Figura 3.14. Datos mensuales de humedad relativa para la ciudad de Quito y Tulcán 2012-
2014. ...............................................................................................................................90
xi
CONTENIDO DE TABLAS
Capítulo 1
Tabla 1.1. Características técnicas de sensor DHT11......................................................16
Tabla 1.2. Características técnicas de sensor DHT22......................................................17
Tabla 1.3. Características técnicas del sensor BMP280. .................................................18
Tabla 1.4. Características técnicas del sensor BMP180. .................................................19
Capítulo 2
Tabla 2.1. Análisis de selección del microcontrolador Arduino. ........................................26
Tabla 2.2. Análisis de selección para sensores de humedad y temperatura. ...................27
Tabla 2.3. Análisis de selección para sensores de altitud. ...............................................27
Tabla 2.4. Localidades del Ecuador utilizadas para estudio de propiedades psicrométricas.
........................................................................................................................................37
Tabla 2.5. Cálculo de error porcentual entre datos de altitud con Google Earth (altitud de
referencia) y datos de altitud con GIS en msnm...............................................................38
Tabla 2.6. Tabla de datos corregidos mediante uso de Excel. .........................................40
Tabla 2.7. Ejemplo de cálculo de humedad relativa variando la temperatura de bulbo
húmedo............................................................................................................................43
Tabla 2.8. Factores de la ecuación de corrección (4) .......................................................46
Tabla 2.9. Resultados de los cálculos respectivos para la obtención de la ecuación de
corrección. .......................................................................................................................47
Tabla 2.10. Resultados de los cálculos para la verificación de la ecuación de corrección.
........................................................................................................................................49
Tabla 2.11. Resultados del cálculo del coeficiente de determinación utilizando varios
métodos. ..........................................................................................................................49
Tabla 2.12. Ecuaciones de Matlab con error promedio en humedad relativa ...................50
Tabla 2.13. Factores de ecuación de corrección 16. ........................................................50
Tabla 2.14. Significado de los códigos mostrados junto al valor de humedad relativa en el
higrómetro digital construido. ...........................................................................................56
Tabla 2.15. Análisis de selección para el proceso de fabricación de la carcasa del
higrómetro digital. ............................................................................................................61
Tabla 2.16. Análisis de selección para los materiales de la impresión 3D. .......................62
Tabla 2.17. Funcionamiento de los botones del higrómetro digital con corrección. ..........64
Tabla 2.18. Características técnicas del higrómetro digital con corrección.......................66
Tabla 2.19. Tabla de costos de fabricación del prototipo e higrómetro digital con corrección.
........................................................................................................................................69
Capítulo 3
Tabla 3.1. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital
H1 y el psicrómetro de flujo continuo a 2806 msnm (Quito). ............................................72
xii
Tabla 3.2. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital
H2 y el psicrómetro de flujo continuo a 2806 msnm (Quito). ............................................72
Tabla 3.3. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital
H3 y el psicrómetro de flujo continuo a 2806 msnm (Quito). ............................................73
Tabla 3.4. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital
H4 y el psicrómetro de flujo continuo a 2806 msnm (Quito). ............................................74
Tabla 3.5. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital
H5 y el psicrómetro de flujo continuo a 2806 msnm (Quito). ............................................75
Tabla 3.6. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital
H6 y el psicrómetro de flujo continuo a 2806 msnm (Quito). ............................................75
Tabla 3.7. Errores porcentuales de humedad relativa de los higrómetros digitales
estudiados a 72msnm. .....................................................................................................77
Tabla 3.8. Valores de varias propiedades psicrométricas del aire con el psicrómetro de flujo
continuo y el higrómetro digital construido sin corrección en la ciudad de Atacames (nivel
del mar). ..........................................................................................................................78
Tabla 3.9. Valores de varias propiedades psicrométricas del aire con el psicrómetro de flujo
continuo y el higrómetro digital sin corrección en la ciudad de Santo Domingo (586 msnm).
........................................................................................................................................79
Tabla 3.10. Valores de varias propiedades psicrométricas del aire con el psicrómetro de
flujo continuo y el higrómetro digital sin corrección en la ciudad de Rocafuerte (950 msnm).
........................................................................................................................................80
Tabla 3.11. Valores de varias propiedades psicrométricas del aire con el psicrómetro de
flujo continuo y el higrómetro digital sin corrección en la ciudad de Quito Centro (2804
msnm). ............................................................................................................................81
Tabla 3.12. Valores de varias propiedades psicrométricas del aire con el psicrómetro de
flujo continuo y el higrómetro digital sin corrección en la ciudad de Quito Sur (2951 msnm).
........................................................................................................................................82
Tabla 3.13. Valores de varias propiedades psicrométricas del aire con el psicrómetro de
flujo continuo y el higrómetro digital sin corrección en la ciudad de Tulcán (3046 msnm).
........................................................................................................................................83
Tabla 3.14. Datos de altitudes de referencia y altitudes obtenidas con el higrómetro digital
sin corrección en diferentes localidades del Ecuador. ......................................................84
Tabla 3.15. Valores de varias propiedades psicrométricas del aire con el psicrómetro de
flujo continuo y el higrómetro digital con corrección en la ciudad Quito (2804 msnm). .....85
Tabla 3.16. Valores de varias propiedades psicrométricas del aire con el psicrómetro de
flujo continuo y el higrómetro digital corregido en la ciudad de Tulcán (2804 msnm). ......86
Tabla 3.17. Datos de altitudes de referencia y altitudes obtenidas con el higrómetro digital
con corrección en diferentes localidades del Ecuador......................................................87
Tabla 3.18. Resultados del cálculo de error porcentual promedio de humedad relativa con
las tres ecuaciones de corrección. ...................................................................................88
Tabla 3.19. Valores de varias propiedades psicrométricas del aire con el psicrómetro de
flujo continuo y el higrómetro digital sin corrección la ciudad de Tulcán (3046 msnm). ....90
xiii
RESUMEN
xiv
ABSTRACT
Keywords: Hygrometer, digital hygrometer, relative humidity, wet bulb temperature, dry
bulb temperature.
xv
DISEÑO Y CONSTRUCCIÓN DE UN HIGRÓMETRO DIGITAL PARA
ALTITUDES DE HASTA 3000 METROS SOBRE EL NIVEL DEL
MAR PARA EL LABORATORIO DE TERMODINÁMICA FIM-EPN.
INTRODUCCIÓN
Las mediciones de los valores de humedad relativa en los últimos años han estado
directamente relacionadas con la facilidad de medición, siendo el higrómetro digital el
dispositivo de mayor uso para obtener estos valores [1]. Estos medidores al no tener un
análisis correspondiente al cambio de altitudes, es decir desde el lugar que fueron
fabricados hasta lugares con una altitud superiores a los 1000 msnm, presentan errores en
su respectiva medición. Por lo tanto, confiar únicamente en el valor marcado no es muy
conveniente, si se desea realizar procesos de producción sensibles al cambio de humedad
relativa.
En la actualidad el uso de la humedad relativa en diversos procesos industriales y
científicos es de mucha utilidad, debido a que una medición errónea en estos valores afecta
indirectamente a la salud de las personas. La característica principal de los higrómetros
digitales es obtener los valores de humedad relativa y temperatura ambiente de una
manera rápida reduciendo los tiempos de búsqueda de variables que afecten a cualquier
proceso, industrial o científico.
Uno de los objetivos de este trabajo es obtener un higrómetro digital que relacione la altitud
con la humedad relativa, mediante la utilización de un microcontrolador programable de
software abierto con sensores de humedad, temperatura y altitud. El higrómetro fabricado
debe presentar un margen de error promedio porcentual bajo en humedad relativa en
comparación a otros higrómetros existentes en el mercado ecuatoriano.
Para que el higrómetro digital fabricado presente un bajo margen de error se realiza una
corrección al valor de altitud del sensor BMP180, y otra corrección indirecta a la humedad
relativa (corrección a la temperatura de bulbo húmedo experimental en función de la real).
De esta manera el higrómetro fabricado realiza mediciones con un bajo error porcentual en
los valores de humedad relativa para altitudes de hasta 3000 msnm. El higrómetro digital
fabricado será propiedad del Laboratorio de Termodinámica de la Facultad de Ingeniería
Mecánica de la Escuela Politécnica Nacional (FIM-EPN).
1
Objetivo general
Diseñar y construir un higrómetro digital para altitudes de hasta 3000 msnm para el
laboratorio de Termodinámica FIM-EPN.
Objetivos específicos
2
1. MARCO TEÓRICO
Este capítulo aborda de manera general los conceptos relacionados con el estudio de las
propiedades psicrométricas del aire y los diferentes tipos de higrómetros más utilizados.
También se explica los sensores de humedad relativa, temperatura y altitud disponibles en
el mercado ecuatoriano que trabajan con la plataforma Arduino.
1.1. Psicrometría
La psicrometría se define como la ciencia que estudia las propiedades termodinámicas del
aire real (aire con contenido de humedad), así como también el efecto que causa la
humedad atmosférica sobre los materiales y el confort humano. [2]
Las cartas psicrométricas son diagramas de doble entrada en las cuales se relacionan
varias propiedades del aire, como la humedad relativa, punto de rocío, temperatura,
volumen específico, etc. Se utilizan para estudiar el cambio de dichas propiedades
termodinámicas con variación de humedad. [3]
Figura 1.1. Carta psicrométrica a temperaturas normales y presión barométrica de 101.325 kPa (al
nivel del mar).
(Fuente: [6] )
3
En una carta psicrométrica se encuentran varias propiedades del aire, a una presión o
altitud específica. Las propiedades de mayor importancia para este trabajo son las
siguientes:
a) Humedad específica
b) Humedad absoluta
Por lo tanto, la humedad relativa es la relación entre la presión parcial del vapor de agua
contenido en una muestra de aire a la presión y temperatura del punto de rocío, y la presión
de agua saturada a la presión y temperatura de bulbo seco. En otras palabras, es la relación
entre el contenido de agua de una mezcla respecto al contenido de agua en saturación de
la misma mezcla en iguales condiciones de presión y temperatura, este valor es expresado
en porcentaje (%) [4].
En las cartas psicrométricas la humedad relativa se muestra con líneas curvas con
tendencia creciente de izquierda a derecha, manteniendo siempre un valor constante. Tal
como se muestra en la Figura 1.2.
4
Figura 1.2. Líneas de humedad relativa.
(Fuente: Propia)
5
bulbo se encuentra cubierto por un material higroscópico húmedo, con el fin de medir la
energía que se desprende por la evaporación del agua (proceso adiabático). [7]
En las cartas psicrométricas la temperatura del punto de rocío se muestra con líneas
horizontales de izquierda a derecha, manteniendo siempre un valor constante. Tal como
se muestra en la Figura 1.5.
6
1.2. Higrómetro
Es un instrumento que mide valores de humedad relativa y temperatura de bulbo seco en
el ambiente, su aplicación depende de la altitud en la cual se realice las mediciones [8].
Existen diferentes tipos de higrómetros o sensores de humedad relativa, de los cuales se
presentan los más importantes:
Para encontrar la humedad relativa y la temperatura de bulbo seco, los higrómetros utilizan
elementos mecánicos y escalas graduadas diseñadas para las diferentes altitudes sobre el
nivel del mar. Los mecanismos utilizados no deben presentar mayor complicación en su
funcionamiento y deben resistir ambientes corrosivos. El principio de funcionamiento
utilizado en los higrómetros mecánicos se muestra en la Figura 1.6.
Estos higrómetros tienen un error en su medición de hasta 3% según los datos del
fabricante, trabajan en rangos de temperaturas de entre -20 °C y 70 °C y humedades
relativas entre 15% y 95% [8].
7
En la Figura 1.7 se observa un higrómetro mecánico que utiliza cabello humano de la
Facultad de Ingeniería Química de la Escuela Politécnica Nacional (FIQ-EPN). Este
higrómetro utiliza un papel que presenta escalas de tiempo, humedad relativa y
temperatura de bulbo seco.
Los higrómetros de bulbo seco y bulbo húmedo han ido mejorando ya que no se utiliza
termómetros comunes, sino se ha optado por utilizar termocuplas o termo-resistencias que
calculan la temperatura con alta tasa de confiabilidad. De igual manera, para la temperatura
de bulbo húmedo ya no se utiliza algodones o mechas comunes porque se llenan de polvo
y no absorben adecuadamente el agua, por lo cual, se está utilizando materiales cerámicos
porosos que se encuentran en contacto con el agua y por el fenómeno de capilaridad
siempre humedecen al termopar [8].
8
En la Figura 1.8 se observa el fenómeno empleado por estos higrómetros:
Los higrómetros de bulbo seco y bulbo húmedo son capaces de determinar la humedad
relativa del ambiente en un rango de 10% a 90% y en temperaturas de entre 0 °C a 90 °C.
La humedad calculada tiene un porcentaje de 0.5% de error, por lo cual estos higrómetros
son utilizados para calibración de otros equipos. [8]
Figura 1.9. Higrómetro de bulbo seco y bulbo húmedo (psicrómetro de flujo continuo) de la FIM-
EPN.
(Fuente: Propia)
9
1.2.3. Higrómetros por condensación
Estos tipos de higrómetros utilizan como base fundamental el punto de rocío para calcular
la humedad relativa del ambiente. Para esto se utiliza una cámara compuesta por una
resistencia calefactora, un espejo, un equipo de refrigeración, una fuente luminosa y dos
foto-resistencias, en la cual circula una mezcla gaseosa de aire que humedece el espejo y
con la utilización de los demás componentes se logra temperar su superficie hasta llegar
al punto de rocío. Luego, se mide la temperatura en que se encuentra y utilizando
ecuaciones y cartas psicrométricas se determina el valor de la humedad relativa [8].
Este dispositivo trabaja en rangos de entre -70 °C a 40 °C para encontrar el valor del punto
de rocío, y presenta un error de 0.5% igual que el higrómetro de bulbo seco y bulbo
húmedo, según el fabricante [8].
Para determinar la humedad relativa el higrómetro produce un haz de luz infrarroja sobre
el material sólido y con la ayuda de una fotocélula se mide la energía reflejada en 2
longitudes de onda. La intensidad de estas dos señales se compara utilizando circuitos
10
electrónicos obteniéndose una variación que es directamente proporcional a la cantidad de
agua que presenta el material. [9]
Las longitudes de onda que ayudan a determinar la humedad son diferentes para cada
material, por lo que utilizando este equipo para una sola aplicación ayudaría que los
sensores estén siempre calibrados [9].
Estos dispositivos están compuestos por un tubo metálico, el cual se encuentra enrollado
de tejidos como la fibra de vidrio que a su vez están en contacto con las sales higroscópicas
saturadas (cloruro de litio y agua). Mediante la utilización de dos hilos de oro que se
envuelven en el tejido de fibra de vidrio y una resistencia limitadora de corriente, se genera
una disipación de potencia que evapora el agua salina. Gracias a este fenómeno la
temperatura que se obtiene se dispersa por todo el tubo metálico, en donde utilizando una
termo-resistencia de platino se conoce su valor. Este valor de temperatura es proporcional
a la temperatura de punto de rocío, por lo cual con la ayuda de las cartas psicrométricas
se determina la humedad relativa.
11
Figura 1.12. Esquema de funcionamiento de un higrómetro de sales higroscópicas.
(Fuente: [8] )
12
aplicaciones dentro del laboratorio y principalmente para la realización de prácticas y tesis.
En la Figura 1.13 se observa los higrómetros mencionados.
13
Estas placas electrónicas utilizadas por Arduino son “Placas de Circuito Impreso” o por sus
siglas en ingles PCB (Printed Circuit Board), es decir consta de una serie de circuitos
electrónicos internos, que facilitan al usuario el desarrollo de diferentes aplicaciones sin
tener que preocuparse de las conexiones del microcontrolador.
1.3.2. Software
Existes varios modelos de placas Arduino, todos diseñados para diferentes aplicaciones y
con características variadas (número de entradas y salidas, tamaño físico, modelo del
microcontrolador). A pesar de que tienen diferentes características, comparten cosas en
común por ejemplo librerías, documentación ya que los microcontroladores son de la
misma familia (Microcontroladores AVR marca Atmel) [10]. En la Figura 1.16 se observa
los diferentes tipos de placas Arduino que existen:
14
Figura 1.16. Tipos de placas Arduino.
(Fuente: [10] )
Existen varios sensores de humedad y temperatura que trabajan en conjunto con Arduino,
pero los más utilizados son los sensores DHT11 y DHT22. Estos sensores brindan
resultados óptimos y de manera sencilla, además son fáciles de encontrar en el mercado
ecuatoriano. Están compuestos internamente por un sensor de humedad capacitivo y un
termistor, contiene un circuito integrado el cual tiene la función de convertir los valores
analógicos a digitales, ver ANEXO I.
a) Sensor DHT11
Estos sensores están calibrados y aseguran alta confiabilidad en sus resultados, siempre
y cuando no sean utilizados a altitudes superiores a los 1000 msnm, ya que a estas
altitudes presentan varios errores en sus mediciones. Están compuestos por un
microcontrolador de 8 bits integrado el cual brinda resultados instantáneos. En la Tabla 1.1
se encuentran las características técnicas de funcionamiento según los fabricantes y en la
Figura 1.17 se observa el sensor DHT11 con sus pines de conexión.
15
Tabla 1.1. Características técnicas de sensor DHT11
(Fuente: [11] )
b) Sensor DHT22
Estos sensores al igual que el sensor DHT11 presenta gran confiablidad y estabilidad en
la recolección de datos de humedad y temperatura, pero de igual manera hasta altitudes
no mayores de 1000 msnm. El sensor de humedad capacitivo y el termistor se encuentran
conectados a un microcontrolador de 8 bits y los valores resultantes son digitales. En la
Tabla 1.2 se observa las características técnicas dadas por el fabricante y en la Figura 1.18
el sensor DHT22 con sus respectivos pines de conexión.
16
Tabla 1.2. Características técnicas de sensor DHT22.
(Fuente: [11] )
De acuerdo con estas características, este sensor es más preciso que el sensor DHT11,
además, cuenta con un mayor rango de medición. Principalmente este sensor es utilizado
con mayor frecuencia porque cuenta con una resolución en su medición de 0,1 para
humedad relativa y temperatura, es decir, los resultados finales contienen decimales, lo
que hace que sean más precisos.
17
a) Sensor BMP280
(Fuente: [11] )
b) Sensor BMP180
Es un sensor de presión barométrica digital, el cual se usa principalmente como altímetro
ya que es capaz de medir la presión atmosférica en función de la altura en la que se
encuentre, además de la temperatura ambiente. En la Tabla 1.4 se observan varias las
características técnicas.
18
Tabla 1.4. Características técnicas del sensor BMP180.
(Fuente: [11])
Este sensor también consta de un regulador de 3.3V incorporado que permite un correcto
funcionamiento a partir de una fuente de 5 V, también consta de un sistema de Ultra-bajo
consumo de energía (0.5uA a 1 Hz). En la Figura 1.20 se observa el sensor con sus pines
de conexión.
a) Error Absoluto
19
b) Error Relativo
𝑬𝒓𝒓𝒐𝒓 𝑨𝒃𝒔𝒐𝒍𝒖𝒕𝒐
𝑬𝒓𝒓𝒐𝒓 𝑹𝒆𝒍𝒂𝒕𝒊𝒗𝒐 =
𝑽𝒂𝒍𝒐𝒓 𝑹𝒆𝒂𝒍
Este error relativo se puede multiplicar por 100% para expresarlo como porcentaje y se
denota como error porcentual.
𝑬𝒓𝒓𝒐𝒓 𝑨𝒃𝒔𝒐𝒍𝒖𝒕𝒐
𝑬𝒓𝒓𝒐𝒓 𝑹𝒆𝒍𝒂𝒕𝒊𝒗𝒐 𝑷𝒐𝒓𝒄𝒆𝒏𝒕𝒖𝒂𝒍 = ∗ 𝟏𝟎𝟎%
𝑽𝒂𝒍𝒐𝒓 𝑹𝒆𝒂𝒍
La aproximación por mínimos cuadrados más simple que se utiliza es ajustar a los puntos
obtenidos en una línea recta, para lo cual se tiene la siguiente ecuación:
𝑦 = 𝑎0 + 𝑎1 𝑥 + 𝑒
20
Donde 𝑎0 𝑦 𝑎1 son coeficientes que representan la intersección con el eje 𝑦 y 𝑒 representa
el error entre los valores reales y los experimentales. Para ajustar de mejor manera la línea
se minimiza la suma de los errores de los datos [12].
Con los valores calculados se procede a realizar una validación o cuantificación del error
para la regresión lineal. Se calcula el valor de 𝑟 que representa el coeficiente de correlación,
es decir, si este coeficiente se acerca a 1 garantiza que los datos de la ecuación son
adecuados, mientras que si se acercan a 0 la ecuación encontrada es ineficiente y el ajuste
no sería valido.
b) Interpolación
Este método se utiliza cuando los valores graficados aseguran un porcentaje de error
pequeño, es decir que son muy precisos, con lo cual se grafica una curva que pase por
todos los puntos de forma directa [12]. El procedimiento más adecuado es la interpolación
polinomial, es decir hallar un polinomio de grado n que se ajuste al número de puntos de
la gráfica, en la Figura 1.22 se observa una imagen con varios ejemplos de interpolación
polinomial.
Figura 1.22. Ejemplos de interpolación polinomial: a) de primer grado (lineal) que une dos puntos,
b) de segundo grado (cuadrática o parabólica) que une tres puntos y c) de tercer grado (cúbica)
que une cuatro puntos.
(Fuente: [12])
Estos programas utilizan las ecuaciones del método manual, por lo cual existen varias
bibliotecas y paquetes de software que ayudan con el ajuste de curvas, de los cuales los
más utilizados son el Excel y el Matlab.
a) Excel
Excel es una hoja de cálculo fabricada por Microsoft Inc. Es un tipo de software matemático
computarizado que permiten realizar operaciones en renglones y columnas de datos.
Cuenta con varias funciones numéricas interconstruidas para resolver ecuaciones, ajustes
21
de curvas y optimización. Además, incluye un lenguaje macro (VBA) para realizar cálculos
numéricos y consta de herramientas para la visualización de gráficos, diagramas con los
cuales se complementan los análisis numéricos [12].
Excel brinda varias funciones para el ajuste curvas, de las cuales la función más útil es el
análisis de regresión por mínimos cuadrados, además existen varias funciones
interconstruidas para el método de regresión, las cuales se las observa en la Figura 1.23.
Figura 1.23. Funciones de Excel interconstruidas relacionadas con el ajuste de datos por
regresión.
(Fuente: [12])
A parte de las funciones interconstruidas existen comandos muy utilizados para el ajuste
de curvas, los más utilizados son el Trendline y el Data Analysis Toolpack. Estos comandos
permiten obtener varias tendencias en una gráfica, pueden ser ajustes polinomiales,
logarítmicos, exponenciales, y de potencia [12].
En la Figura 1.24 se observa un ejemplo de ajuste por regresión con un modelo logarítmico.
22
b) Matlab
Matlab es una software fabricado por Mathworks Inc., fundada por los analistas numéricos
Cleve Moler y John N. Little. Cuenta con varias funciones matemáticas de resolución de
matrices, funciones numéricas, cálculos simbólicos y herramientas para visualización de
gráficos [12]. Para el ajuste de curvas Matlab tiene varias funciones preconstruidas las
cuales se las observa en la Figura 1.25.
Figura 1.26. Ejemplo de ajuste polinomial realizado en Matlab. a. Se observa los valores
interpolados linealmente. b. Ajuste polinomial de quito grado.
(Fuente: [12])
Matlab a diferencia del Excel incluye herramientas de ajuste curvas básico y cajas de
herramientas donde se realizan operaciones estadísticas especializadas. La caja de
herramientas para ajuste de curvas contiene una interfaz gráfica que ajusta los datos sin
necesidad de usar solo polinomios, por ejemplo, transformada discreta de Fourier e
interpolación de datos con segmentaria cúbica [12].
23
2. METODOLOGÍA
En el presente capítulo se detalla la evaluación de los higrómetros digitales presentes en
el mercado a nivel del mar y a 2800 msnm, tomando como referencia el psicrómetro de
flujo continuo (Ps FC, ver Figura 1.9) del laboratorio de termodinámica de la FIM-EPN.
También, se muestra el proceso de fabricación del higrómetro digital para altitudes de entre
0 y 3000 (msnm) con corrección, es decir un higrómetro capaz de obtener valores de
humedad relativa, temperatura y altitud con bajo error porcentual. Este proceso de
fabricación implica varias tomas de datos de propiedades psicrométricas (humedad
relativa, temperatura de bulbo seco, temperatura de bulbo húmedo, altitud) en diferentes
altitudes, y posteriormente la corrección de señal del higrómetro digital construido. En la
Figura 2.1 se observa un diagrama de flujo de las actividades desarrolladas en este trabajo,
el diagrama completo se observa en la Figura 2.46, p. 70.
24
termómetros de mercurio del psicrómetro de flujo continuo es de ±0.2°C, valor que se
encuentra dentro del rango de aceptación de la norma ASTM E337-15 [13].
Los dispositivos utilizados para el análisis previo a la realización del prototipo de higrómetro
digital sin corrección se especifican con las siguientes abreviaturas (H1, H2, H3, H4, H5,
H6), donde H1 y H2 pertenecen al laboratorio de termodinámica de la FIM-EPN, H3
pertenece a la Facultad de Ingeniería Química de la EPN, H4 pertenece al laboratorio de
fluidos de la FIM-EPN, y por último H5 y H6 pertenece al LAEV (laboratorio de análisis de
esfuerzos y vibraciones) de la FIM-EPN. Estos higrómetros se los observa en la Figura 2.2.
Las mediciones son realizadas en dos altitudes diferentes, una a nivel del mar en la que se
espera que los higrómetros funcionen adecuadamente, y otra a una altitud de 2806 msnm
25
(Quito) en las que se verifica un funcionamiento con errores considerables. En estas dos
altitudes se comparan los valores de humedad relativa (HR) y la temperatura de bulbo seco
(TBS) de cada uno de los higrómetros digitales, con los resultados obtenidos en las cartas
psicrométricas (CP). Los resultados de esta verificación se observan en la sección 3.1.2,
pp. 71-77.
Cabe señalar que las propiedades psicrométricas son muy importantes para una infinidad
de procesos industriales y científicos, un margen de error alto presentaría una afectación
importante. El error en la medición de la humedad relativa se ha venido cometiendo de
manera no intencionada desde que se introdujeron los psicrómetros digitales al mercado.
Existen varios tipos de placas Arduino y sensores los cuales se describieron en el capítulo
1, por lo cual se realiza un análisis de selección para las diferentes placas Arduino,
sensores de humedad y temperatura y los sensores barométricos presentes en el mercado
ecuatoriano. Con este análisis se define cuáles serían los adecuados para la construcción
del prototipo de higrómetro digital sin acondicionamiento de señal. Para el análisis se
utilizaron valores que varían de 1 a 5 en donde, 5 es el valor óptimo para la característica
evaluada, y 1 es el valor más bajo. En la Tabla 2.1 se observa el análisis realizado para la
selección de la placa Arduino adecuada.
26
Con este análisis se determina que el Arduino uno cumple con los requerimientos
necesarios para el presente trabajo.
En la Tabla 2.2 se observa el análisis para la selección del sensor de humedad relativa y
temperatura.
Con el análisis realizado en la Tabla 2.2 se determina que el sensor DHT22 es la mejor
opción y cumple con las características necesarias para la realización de este trabajo.
Con el análisis de la Tabla 2.3 se determina que el sensor BMP180 cumple con las
especificaciones necesarias para el presente trabajo. Después de realizar esta selección
en la placa Arduino y en los sensores se procede a la construcción del prototipo de
higrómetro digital sin acondicionamiento de señal con la ayuda de otros materiales que se
mencionan a continuación.
En este apartado se detallan los materiales elegidos para el montaje del prototipo de
higrómetro digital sin corrección, y también los elementos electrónicos seleccionados
anteriormente.
a) Sensor de humedad
El sensor AM 2302 (DHT22) permite obtener las propiedades psicométricas del aire
(humedad relativa, temperatura de bulbo seco). En la Figura 2.3 se muestra el sensor AM
2302 en su forma física.
27
Figura 2.3. Sensor AM2302, estructura física.
(Fuente: Propia)
Este sensor presenta 4 pines de conexión, de los cuales uno se anula debido a que no
presenta variación en las mediciones si es o no conectado. Por lo tanto, se utilizará los 3
pines restantes los cuales se especifican en el ANEXO I y se ilustra en el apartado 4.1 de
ese anexo.
b) Sensor de Altitud
El sensor de altitud BMP 180, que se muestra en la Figura 2.4, permite obtener la altitud
del lugar donde se realice las respectivas mediciones tomando en cuenta como referencia
los 0 msnm. Las características del sensor, los rangos de medición y funcionamiento se
muestran en el ANEXO II.
Cabe mencionar que el sensor de altitud no censa presión atmosférica directamente, por
lo cual, para la determinación de la presión atmosférica se trabaja en base a la ecuación
(1) [4]. Esta ecuación al estar en función de la altitud permite obtener la presión atmosférica
propia del lugar donde se encuentre el sensor y este valor es utilizado para la mayoría de
los cálculos que sen explicados en este proyecto.
28
c) Pantalla de cristal líquido 1602
La Figura 2.5 muestra la pantalla en la cual se trabaja para obtener los valores de humedad
relativa, temperatura de bulbo seco (temperatura ambiente), altitud y presión atmosférica
sin corrección.
Esta pantalla presenta un potenciómetro que permite variar el contraste y nitidez de los
resultados, además incluye botones los cuales facilitan su funcionamiento. La pantalla
permite una conexión directa al Arduino, es decir cuenta con pines que están distribuidos
de igual manera que los sockets de la placa Arduino uno, esto facilita su programación. Las
especificaciones de esta pantalla se muestran con detalle en el ANEXO III.
Para obtener datos de manera digital y análoga se utilizará una placa Arduino uno, ya que
es de fácil programación y de un costo relativamente bajo, su forma física se muestra en la
Figura 2.6. Los datos de operación y esquema de la placa se encuentran detallados en el
ANEXO IV.
29
e) Agua desmineralizada
Para obtener los datos de humedad relativa reales es necesario conocer el valor de la
temperatura de bulbo húmedo utilizando un material higroscópico y agua. Como el agua
presenta altas cantidades de solidos se utiliza agua desmineralizada para evitar que los
minerales sedimenten alrededor del bulbo de termómetro de mercurio. El agua utilizada se
muestra en la Figura 2.7.
f) Ventilador 8x8
Para reducir el tiempo de estabilización de los bulbos del psicrómetro de flujo continuo se
usa un ventilador de 8x8 cm2 que produce una velocidad del viento de entre 0.5 a 0.9 m/s.
Este trabaja con un voltaje constante de 5 v y un flujo de corriente de 0.11 [A], la forma
física de este ventilador se muestra en la Figura 2.8. Este flujo es únicamente usado para
reducir el tiempo la estabilización de los bulbos y de los sensores.
30
g) Cartón pluma
El cartón pluma es usado para fabricar la caja que evita que la radiación solar llegue
directamente a los bulbos del psicrómetro de flujo continuo, además facilita la redirección
del flujo de aire generado por el ventilador. Este cartón está formado por 2 materiales, uno
es un polímero y el otro es cartón común. Como presenta una baja conductividad térmica
no afecta directamente en las mediciones de temperaturas. La forma física del cartón pluma
se muestra en la Figura 2.9.
Para la fabricación del prototipo electrónico de higrómetro digital sin corrección se utilizan
los materiales mencionados en la sección 2.3.2, los cuales se colocan en una caja fabricada
de cartón pluma. Esta caja se realiza con el fin de lograr que los sensores se encuentren
en las mismas condiciones de temperatura de los bulbos del psicrómetro de flujo continuo.
El diseño y las dimensiones de la caja se observan en la Figura 2.10.
31
Figura 2.10. Dimensiones de la caja fabricada de cartón pluma.
(Fuente: Propia)
32
b) Instalación de sensores en la caja de cartón pluma y conexión en placa Arduino
uno.
El sensor de humedad se ubica junto al sensor de altitud dentro de la caja de cartón pluma.
Debido a la distribución de los sensores en la caja, el flujo de aire generado por el ventilador
llega antes a los sensores y luego a los bulbos. Esto evita que la humedad del material
higroscópico modifique las lecturas de humedad relativa y temperatura del sensor DHT22.
La ubicación de los sensores en la caja se observa en la Figura 2.12.
El pin bidireccional del sensor DHT22 se conecta a la entrada digital (n°3) del
microcontrolador como se muestra en la Figura 2.13, este puerto permite la lectura de la
temperatura y humedad relativa del aire. Además, el sensor requiere una entrada de voltaje
continuo en un rango de 3.3 a 5.5 Voltios, para una adecuada medición. Si, el voltaje se
encuentra fuera de este rango puede presentar interferencia en las mediciones, por lo
tanto, el valor de alimentación para el sensor es de 5 voltios.
33
El sensor BMP180 trabaja con dos entradas de voltaje, una de 3.3 V y otra de 5 V para su
funcionamiento, por lo tanto, se utiliza 5 V al igual que el sensor DHT22, también tiene una
conexión a tierra para protección de saltos de corriente. Este sensor tiene salidas SDA Y
SDL que van conectados a la placa Arduino en las entradas análogas A4 y A5
respectivamente como se observa en la Figura 2.14.
En la Figura 2.14 se visualiza las conexiones del sensor BMP180 con la placa Arduino uno.
Los valores censados se muestran en una pantalla de cristal líquido 1602, ubicada en la
parte superior de la caja, conectada a la placa Arduino, tal como se muestra en la Figura
2.15. Cuando se enciende la pantalla se muestra la humedad relativa sin corrección junto
con la temperatura ambiental. También, presionando el botón “down” se observa la altitud
y presión atmosférica.
34
c) Código para el prototipo de higrómetro digital sin corrección
En este apartado, luego de realizadas las conexiones en la placa Arduino uno, se explica
el código que ejecuta la placa Arduino para la medición con los sensores DHT22 y BMP180,
este código se realiza en el programa IDE de la empresa Arduino.
35
Al inicio se define p0 como variable flotante, esta variable define la presión de referencia
para el sensor de altitud. Al colocar p0 con un valor de 101365 [Pa] se considera como
valor de referencia 0 [msnm].
Las variables involucradas en este código son la humedad relativa (Hr), temperatura (T),
presión (P) y altitud (A). Estas variables están representadas como punto flotante
(FLOATES) para que los resultados presenten decimales. El código para obtener los
valores de los sensores se observa en la Figura 2.17.
Las propiedades psicrométricas del aire analizadas en este apartado son la humedad
relativa, temperatura de bulbo seco y temperatura de bulbo húmedo, estas propiedades
varían dependiendo del lugar de análisis, también se calcula los valores de altitud.
Utilizando el prototipo electrónico de higrómetro digital construido el cual se menciona en
la sección 2.3.3, pp. 31-36, se realizan mediciones en varias localidades del Ecuador con
respecto a humedad relativa y temperatura ambiente.
36
Tabla 2.4. Localidades del Ecuador utilizadas para estudio de propiedades psicrométricas.
(Fuente: Propia)
Los datos de altitudes de las localidades analizadas se obtienen de Google Earth, estos
datos se aprecian en la Tabla 2.4 y varían desde los 0 msnm hasta los 3040 msnm.
En cada lugar, se toman datos de valores reales (psicrómetro de flujo continuo) y de valores
experimentales (prototipo electrónico de higrómetro digital sin corrección), con el fin de
conocer los errores presentes en altitud y humedad relativa. Estas mediciones se las realiza
cada 20 minutos en el trascurso del día. Las tablas y gráficas obtenidas se observan en la
sección 3.1.3, pp. 77-84, en las cuales se evidencia los errores en valores de altitud y
humedad relativa.
Para realizar la corrección en altitud se considera como altitud cercana a la real a los datos
satelitales de Google Earth (altitud de referencia). Estos valores son comparados con datos
satelitales del GIS (Geographic Information System) para varias localidades del país, con
el objetivo de verificar el error que presenta Google Earth en relación con la altitud. Con
esta verificación se tiene que los datos son confiables para el desarrollo del presente
trabajo, estos valores se muestran en la Tabla 2.5.
37
Tabla 2.5. Cálculo de error porcentual entre datos de altitud con Google Earth (altitud de referencia)
y datos de altitud con GIS en msnm.
-78,826912 -0,561996 2846 2831 0,53% -78,811844 0,213959 1085 1076 0,83%
-78,826912 -0,554474 2444,34 2528 3,42% -78,811844 0,244095 780 772 1,03%
-78,826912 -0,546936 2254 2260 0,27% -78,796776 -0,50927 2995 2973 0,73%
-78,826912 -0,539406 2031 2049 0,89% -78,77417 -0,501732 3023 3005 0,60%
-78,826912 -0,531868 2025 1996 1,43% -78,314621 -0,561996 4168 4172 0,10%
-78,826912 -0,524338 1887 1910 1,22% -78,299553 0,100952 3879 3883 0,10%
-78,826912 -0,5168 1827 1852 1,37% -78,299553 0,131088 3939 3934 0,13%
-78,826912 -0,50927 1973 1960 0,66% -78,826912 0,251625 831 823 0,96%
-78,826912 -0,494202 2018 2006 0,59% -78,826912 0,259163 953 952 0,10%
-78,826912 -0,441467 2090 2079 0,53% -78,826912 0,266685 1133 1130 0,26%
-78,826912 -0,396271 1583 1586 0,19% -78,826912 0,304367 1244 1251 0,56%
-78,826912 -0,373665 1335 1334 0,07% -78,819374 -0,418861 1994 1983 0,55%
-78,826912 -0,335999 1504 1512 0,53% -78,819374 -0,411339 1765 1767 0,11%
-78,826912 -0,23806 1887 1893 0,32% -78,819374 -0,403793 1632 1633 0,06%
-78,826912 0,349564 1767 1769 0,11% -78,819374 -0,396271 1656 1658 0,12%
Promedio 0,61%
(Fuente: Propia)
El error porcentual promedio que presenta Google Earth en la altitud es de 0,61%, por lo
cual se consideran validos los datos que presenta esta plataforma.
En la Tabla 3.14 del capítulo 3, p. 84, se observan los datos de altitud de referencia vs los
calculados, con estos valores se calcula la ecuación de corrección en altitud.
38
Figura 2.18. Gráfica de altitudes de referencia y experimentales, ecuación de corrección en altitud.
(Fuente: Propia)
Con esta gráfica se obtiene la ecuación lineal, la cual permite determinar los valores de
altitud corregida en el higrómetro digital para localidades superiores a los 1000 msnm, se
la observa en la ecuación (2):
Donde:
Donde:
Con las ecuaciones (2) y (3) se realiza una verificación de corrección, para esto se utiliza
Excel como herramienta de cálculo, obteniendo así, resultados óptimos y reduciendo el
error considerablemente.
39
Los valores que se muestran en Tabla 2.6 se obtienen a partir de la ecuación (2) y (3). Los
datos de altitudes superiores a los 1000 msnm son corregidos con la ecuación (2), mientras
que los datos de altitudes menores a los 1000 msnm son corregidos con la ecuación (3).
Con el uso de estas dos ecuaciones se reduce el error porcentual promedio de altitud a 1,2
%.
Altitud Altitud
Altitud Real Diferencia
Lugar obtenida Corregida Error [%]
[msnm] [msnm]
[msnm] [msnm]
Teleférico
3123 2949,2 3121,04 1,96 0,1%
Quito
Quito
2804 2652,54 2795,99 8,01 0,3%
CCICEV
Santo
531 577,72 556,07 25,07 4,7%
Domingo
PROMEDIO 1,2%
(Fuente: Propia)
Para realizar la corrección se parte del gráfico de una línea recta con el universo de datos
experimentales (temperatura de bulbo húmedo y humedad relativa) obtenidos sin
corrección en las diferentes localidades que se muestran en la Tabla 2.4, p. 37. En la Figura
2.19 se observa los valores de humedad relativa experimental y humedad relativa real de
todas las localidades analizadas (ver Tabla 2.11) graficados en una línea recta. La cantidad
de los valores contenidos por la línea recta son reducidos, por lo tanto, la complejidad al
realizar la corrección a los valores de humedad relativa es eleva y no se recomienda
realizar la corrección directamente a los datos de humedad relativa.
40
Figura 2.19. Curva de humedad relativa real vs humedad relativa experimental (higrómetro
fabricado)
(Fuente: Propia)
Con los resultados de la Figura 2.19 se considera la opción de corregir los valores de
temperatura de bubo húmedo debido a que la temperatura de bulbo seco no presenta
errores porcentuales promedios elevados (ver sección 3.1.3, pp. 77-84).
41
El programa es elaborado siguiendo el procedimiento que se especifica en el libro ASRHAE
Fundamental Handbook, estos pasos se especifican más a detalle en el anexo V.
Figura 2.21. Curva de temperatura de bulbo húmedo real vs temperatura de bulbo húmedo
experimental (cálculo)
(Fuente: Propia)
42
Figura 2.22. Carta psicrométrica para una altitud de 2800 msnm a temperatura normal.
(Fuente: [5])
Tabla 2.7. Ejemplo de cálculo de humedad relativa variando la temperatura de bulbo húmedo.
43
Figura 2.23. Ejemplo de cálculo de valores de humedad relativa variando la temperatura de bulbo
húmedo.
(Fuente: [5], realización: Propia)
Las diferencias de los valores de humedad relativa por cada grado Celsius de temperatura
de bulbo húmedo varía entre 1.5 (%/°C) y 14 (%/°C) para esta carta psicrométrica.
Figura 2.24. Curvas de sensibilidad de humedad relativa por cada grado Celsius de temperatura
de bulbo húmedo (%/°C).
(Fuente: [5], realización: Propia)
44
Por el análisis de la sección 2.4.3, p. 42-45, y el resultado de la Figura 2.21 se realiza una
corrección directamente a la temperatura de bulbo húmedo para obtener los valores de
humedad relativa corregidos.
➢ Excel.
➢ Regresión lineal con mínimos cuadrados.
➢ Programación en Matlab.
a) Uso de Excel
Con los datos obtenidos a partir del higrómetro digital sin corrección (ver sección 3.1.3, pp.
77-84), se obtiene la ecuación de corrección utilizando Excel.
La ecuación de corrección que se obtiene con este análisis se muestra en la Figura 2.25.
Para verificar la confiabilidad de esta ecuación se obtiene el coeficiente de determinación
R2.
32,00
27,00
y = 0,0172x2 + 0,2895x + 6,8399
R² = 0,9976
Th real [°C]
22,00
17,00
12,00
7,00
7,00 9,00 11,00 13,00 15,00 17,00 19,00 21,00 23,00 25,00 27,00
Th experimental [°C]
Mediante Excel se obtiene la ecuación de corrección (4) que corrige los valores de
temperatura de bulbo húmedo.
45
𝑇ℎ𝑐𝑜𝑟𝑟 = 𝐷(1) ∗ 𝑇ℎ2 + 𝐷(2) ∗ 𝑇ℎ + 𝐷(3) (Ec. 4)
Donde:
Factor Valor
D(1) 0,0172
D(2) 0,2895
D(3) 6,8399
(Fuente: Propia)
Para realizar este análisis se trabaja con los datos de temperatura de bulbo húmedo que
se obtienen con el higrómetro digital fabricado y con el psicrómetro de flujo continuo, en la
cual se realiza un ajuste lineal. Las ecuaciones utilizadas se obtienen del libro de métodos
numéricos para ingenieros de Steven Chapra [12].
Para obtener la ecuación de corrección utilizando una regresión lineal por mínimos
cuadrados, se calcula la pendiente de la línea recta usando la ecuación (5). Esta ecuación
usa los datos de temperatura de bulbo húmedo reales y experimentales que se muestran
en la sección 3.1.3, pp. 77-84.
Donde:
𝑛: Número de datos.
46
𝑦𝑖 : Valores de temperatura de bulbo húmedo reales
𝑥𝑖 : Valores de temperatura de bulbo húmedo experimentales
𝑎1 : Pendiente de la línea recta
Σ𝑦𝑖 Σ𝑥𝑖
𝑦̅ = , 𝑥̅ = (Ec. 6)
𝑛 𝑛
Donde:
ao = y̅ − a1 ∗ x̅ (Ec. 7)
Donde:
En la Tabla 2.9 se muestran los resultados para obtener la ecuación de corrección (9).
Tabla 2.9. Resultados de los cálculos respectivos para la obtención de la ecuación de corrección.
(Fuente: Propia)
𝑦 = 𝑎𝑜 + 𝑎1 𝑥 (Ec. 8)
Donde:
Al realizar la regresión por mínimos cuadrados, se presenta una línea recta que no
necesariamente pasa por todos los datos de las mediciones de temperatura de bulbo
húmedo [12]. Para realizar la cuantificación del error de la ecuación (9), se parte de los
coeficientes calculados usando las ecuaciones (10) y (11) respectivamente. Con los
coeficientes se procede al calculó de la desviación estándar usando la ecuación (12).
Donde:
𝑡𝑆
𝑆𝑦 = √𝑛−1 (Ec. 12)
Donde:
𝑆𝑦 : Es la desviación estándar.
𝑟 𝑆
𝑆𝑦 = √𝑛−2 (Ec.13)
𝑥
Donde:
𝑆𝑡 −𝑆𝑟
𝑟2 = (Ec. 14)
𝑆𝑡
Donde:
𝑟 2 : Coeficiente de determinación
𝑟: Coeficiente de correlación
En la Tabla 2.10 se observan los resultados del análisis de verificación de error para la
ecuación de corrección (9).
48
Tabla 2.10. Resultados de los cálculos para la verificación de la ecuación de corrección.
𝑆𝑦 𝑆𝑦 𝑆𝑡 𝑆𝑟 𝑟2 r
𝑥
(Fuente: Propia)
Donde:
𝑟 2 : Coeficiente de determinación
Tabla 2.11. Resultados del cálculo del coeficiente de determinación utilizando varios métodos.
(Fuente: Propia)
c) Uso de Matlab
49
Tabla 2.12. Ecuaciones de Matlab con error promedio en humedad relativa
(Fuente: Propia)
Para analizar el error porcentual promedio de humedad relativa con los datos corregidos
se utilizan los coeficientes de cada ecuación, los resultados del cálculo del error se
observan en la Tabla 2.12.
Con el análisis planteado, al observar las curvas de la Figura 2.26, la ecuación que contiene
la mayor cantidad de datos reales y experimentales de temperatura de bulbo húmedo es la
de 5to grado, además esta ecuación presenta un error promedio de humedad relativa
corregida igual a 3,144%. Por estas dos razones se elige la ecuación (16) de 5to grado
para realizar la corrección en temperatura de bulbo húmedo.
Donde:
Constante Valor
Ar(1) -6,8138E-05
Ar(2) 0,0063888
Ar(3) -0,23327
Ar(4) 4,1441
Ar(5) -34,935
Ar(6) 122,81
(Fuente: Propia)
50
Figura 2.26. Gráfico de ecuación de corrección de temperatura de bulbo húmedo utilizando Matlab.
(Fuente: Propia)
51
En la Figura 2.27 se observa el valor del coeficiente de determinación obtenido a partir de
la ecuación (16) mediante Matlab.
Al evaluar los resultados con las ecuaciones de corrección, se tiene que la mejor opción
entre las tres ecuaciones es la calculada con programación en Matlab. Con la ecuación
(16) se reduce el porcentaje de error en la humedad relativa en varias localidades del
ecuador con altitudes mayores a los 2000 msnm. Los resultados de los errores
porcentuales promedios al utilizar estas ecuaciones se encuentran en la Tabla 3.18 del
capítulo 3, p. 88.
52
2.5.1. Código de funcionamiento del higrómetro digital con corrección
A partir de las ecuaciones obtenidas en la sección 2.4.4, pp. 45-52, se realiza un ajuste a
la programación del higrómetro digital sin corrección.
Primero se realiza una corrección en altitud con las ecuaciones (2) y (3). En la Figura 2.28
se muestra el código para la lectura y corrección de altitud.
Para encontrar el valor de humedad relativa corregido, se parte de los valores medidos de
temperatura de bulbo seco y temperatura de bulbo húmedo sin corrección, de los cuales
se calcula una diferencia. Esta diferencia permite determinar la temperatura de bulbo
húmedo de referencia, en la Figura 2.29 y la Figura 2.30 se observa el código de
programación para este cálculo.
53
Figura 2.30. Código para el cálculo de la temperatura de bulbo húmedo sin la corrección.
(Fuente: Propia)
Utilizando la ecuación (16), p. 50, se elabora el código que corrige la temperatura de bulbo
húmedo, y a partir de estos valores de encuentra la humedad relativa corregida para varias
altitudes. En la Figura 2.31 se observa el código de corrección para la temperatura de bulbo
húmedo y humedad relativa.
Figura 2.31. Código para corrección de la temperatura de bulbo húmedo y cálculo de la humedad
relativa utilizando la temperatura corregida.
(Fuente: Propia)
Para demostrar que el higrómetro digital realiza las correcciones a los valores mostrados
en la pantalla, se programa un código que permite visualizar los datos con corrección y sin
corrección, esta programación se observa en la Figura 2.32. Para variar estos dos valores
se debe presionar el pulsador uno del higrómetro digital fabricado.
54
Figura 2.32. Código para plasmar los valores corregidos y no corregidos en el higrometro
fabricado.
(Fuente: Propia)
En la Figura 2.33 se observa los diferentes códigos que se muestran en la pantalla lcd
Nokia 5110 para diferentes ambientes y diferentes localidades.
Estos códigos varían dependiendo del lugar de medición, en la Tabla 2.14 se detallan los
códigos mencionados.
55
Tabla 2.14. Significado de los códigos mostrados junto al valor de humedad relativa en el higrómetro
digital construido.
Las mediciones se las realiza cada 20 minutos en el trascurso del día, para la ciudad de
Quito se tomaron un total de 13 datos y para la ciudad de Tulcán se tomaron un total de 17
datos. Con estos valores (reales y experimentales) se calcula el error porcentual promedio
tanto en humedad relativa como en temperatura de bulbo seco. Las tablas y gráficas
obtenidas se observan en las secciones 3.1.4, pp. 85-88.
56
2.7.1. Materiales y componentes
Para la construcción del higrómetro digital con corrección se utilizan diferentes materiales,
estos se dividen en componentes internos (microcontrolador, sensores, placas de
conexión) y componentes externos (material base de la carcasa, pantalla, botones).
a) Componentes Internos
Los componentes internos utilizados para el higrómetro digital con corrección son similares
a los utilizados en el prototipo electrónico sin corrección. Además, se agregaron dos placas
o shields de expansión para aumentar el número de conexiones de voltaje y tierra, y una
batería de 9 V. para el funcionamiento de todo el higrómetro digital. A continuación, se
detallan los componentes internos.
Figura 2.34. Placa Arduino uno utilizada en el higrómetro digital con corrección.
(Fuente: Propia)
Para obtener los valores de humedad relativa y temperatura de bulbo seco se utiliza el
sensor DTH22. El cual se encuentra ubicado en la parte superior del higrómetro con
entradas de aire. Para obtener los valores de altitud sobre nivel del mar se utiliza el sensor
BMP180 el cual se posiciona dentro del higrómetro junto a la placa Arduino uno. Estos
57
sensores se seleccionan por que cumplen con las especificaciones necesarias para el
correcto funcionamiento del higrómetro digital. Las tablas del análisis de selección de los
sensores se las observa en la sección 2.3.1, pp. 26-27. En la Figura 2.35 se observa los
sensores mencionados.
Para la conexión de la pantalla y los botones se utiliza dos placas o shields de expansión,
los cuales se conectan directamente en la placa Arduino uno. Estas placas se utilizan con
el fin de aumentar las salidas de voltaje y las conexiones a tierra. En la Figura 2.36 se
observan las placas conectadas al Arduino uno.
58
➢ Fuente de alimentación
Como fuente de alimentación para el Arduino uno se utiliza una batería de 9 voltios la cual
se encuentra posicionada en la parte posterior del higrómetro digital. Se encuentra
conectada al interruptor (ON-OFF) y al Arduino uno mediante un socket para baterías. En
la Figura 2.37 se observa la batería utilizada con el socket de conexión.
b) Componentes externos
Los componentes externos del higrómetro digital con corrección son la carcasa,
pulsadores, interruptor y la pantalla lcd, los cuales se mencionan a continuación:
59
Figura 2.38. Modelo 3D de la carcasa del higrómetro digital con corrección.
(Fuente: Propia)
Arranque
de viruta
4 3 3 5 3 18
(maquina
CNC)
Impresión
4 4 5 5 5 23
3d
Moldeo por
inyección 1 5 4 1 3 15
plástica
(Fuente: Propia)
El proceso de impresión 3D utiliza varios filamentos, los más utilizados en el mercado son
PLA, ABS y PETG.
Los filamentos ABS son materiales de bases petroquímicas, es decir no son amigables con
el medio ambiente, ya que emiten gases nocivos al ser procesados. Este tipo de filamento
no se usa en lugares pequeños y sin ventilación, una gran ventaja es sus altas propiedades
térmicas [16].
Los filamentos PLA son de bases naturales, es decir, no emiten gases nocivos al ser
procesados y son amigables con el medio ambiente. Este tipo de filamentos presenta una
mayor variedad de colores, la desventaja son sus bajas propiedades térmicas (a partir de
los 50 °C presenta una reducción en las propiedades mecánicas) [16].
Los filamentos PETG son los materiales plástico más utilizado en el mundo, entre uno de
sus productos se tiene a las botellas plásticas. Es un filamento que presenta propiedades
mecánicas y térmicas altas y al ser un material brillante se observa claramente cualquier
rayón en la superficie [16].
Para evaluar las características de los materiales utilizados en la impresión 3D, se utiliza
una numeración de 1 a 5 donde el valor optimo es 5 y el mínimo valor es 1.
61
Tabla 2.16. Análisis de selección para los materiales de la impresión 3D.
Con el análisis que se encuentra en la Tabla 2.15 y la Tabla 2.16, se fabrica la carcasa
para el higrómetro digital corregido mediante la utilización de impresión 3D y con filamento
PLA.
En la Figura 2.39 se presenta una pantalla de cristal líquido de 48x84 pixeles, utilizada para
mostrar los valores con corrección y sin corrección de humedad relativa, temperatura de
bulbo húmedo, temperatura de bulbo seco, altitud. Para la presentación final de higrómetro
fabricado se utiliza esta pantalla, debido a que presenta mejores prestaciones que la
pantalla de 16x2 pixeles utilizada en el prototipo de higrómetro digital sin corrección. Las
características eléctricas de funcionamiento se detallan en el ANEXO VI, cabe mencionar
que esta pantalla presenta una luz interna que permite observar los valores en la noche.
62
Figura 2.40. Código elaborado para el encendido de la pantalla Nokia 5110.
(Fuente: Propia)
El higrómetro digital con corrección está formado de tres secciones principales, estas se
ubican de tal manera que conforman un equipo ergonómico, ligero y fácil de utilizar, ver
Figura 2.41.
➢ Sección frontal:
En la parte frontal del higrómetro digital con corrección se encuentra ubicado la pantalla lcd
Nokia 5110, en la cual se visualiza los valores de humedad relativa, temperatura de bulbo
seco, temperatura de bulbo húmedo y altitud. Estas propiedades se encuentran distribuidas
63
de tal manera que la humedad relativa y la temperatura ambiente se las observe con más
relevancia. En la Figura 2.42 se observa los valores de les propiedades antes mencionadas
en la pantalla LCD Nokia 5110.
Figura 2.42. Propiedades mostradas en la pantalla Nokia 5110 del higrómetro digital.
(Fuente: Propia)
Tabla 2.17. Funcionamiento de los botones del higrómetro digital con corrección.
(Fuente: Propia)
➢ Sección posterior
En la parte posterior de higrómetro digital con corrección se ubican tres agujeros, en los
cuales se insertan tres tornillos de sujeción de manera que la unión con la parte frontal sea
la adecuada y soporte cargas externas. También se ubica la cavidad donde se inserta la
64
batería de 9 V cubierta por una tapa sujetada con dos tornillos pequeños. En la Figura 2.43
se observa la parte posterior del equipo.
➢ Parte superior
En la parte superior se encuentra ubicada la cubierta del sensor DHT22. Como el sensor
debe estar en contacto con el aire se fabrica una tapa con entradas de aire, las cuales
están distribuidas uniformemente en toda la superficie (ocho entradas laterales y una
superior). En la Figura 2.44 se observa la cubierta del higrómetro digital corregido.
65
➢ Características físicas
𝐼𝐷𝐷𝑙𝑜𝑤 100 𝜇𝐴
𝐼𝐷𝐷ℎ𝑖𝑔ℎ 640 𝑚𝐴
Exactitud ±6 %HR
Resolución 0,1 °C
Temperatura Rango* 10 30 °C
Exactitud ±2 °C
(Fuente: Propia)
66
2.7.3. Uso del higrómetro digital fabricado
Para el correcto funcionamiento del higrómetro digital con corrección se deben seguir los
siguientes pasos:
a) Advertencias de uso:
• El sensor de humedad relativa del higrómetro digital debe estar expuesto a flujos de
aire, por lo que no se recomienda ponerlo en contacto directo e indirecto con agua.
• No se debe presionar los tres pulsadores al mismo tiempo, ya que esto genera cortos
y dañar el equipo.
• Para reducir el tiempo de respuesta del higrómetro digital se recomienda ponerlo en
contacto con flujos de aire de entre 0.5 y 1 m/s.
• Si el higrómetro después de presionar el interruptor de encendido no muestra valores
en la pantalla, quiere decir que la batería esta baja, por lo cual se la debe cambiar.
67
On
PRESIONAR INTERUPTOR
Energizar
microcontrolador
Mostrar pantalla de
inicio Definir código de
LCD Nokia 5110 permite
medición inicial Si NO Cambiar bateria
visualizar texto
“B2”
NO
SI
Corregir altitud con
Definir código de Definir código de Definir código de Definir código de
A >1000msnm NO ecuación 3 “Código
medición “B2” medición “A3” medición“B2” medición “A3”
de Medición X3”
SI
Presionar Presionar Presionar
NO pulsador 1 por 2 pulsador 2 por 2 pulsador 3 por 2
A<1000msnm segundos segundos segundos
SI SI SI
SI
NO
SI
Eliminar
Definir código de
correcciones en Apagar luz de LCD Encender luz de LCD
medición “A0”
A y HR Nokia 5110 Nokia 5110
PRESIONAR INTERRUPTOR
Figura 2.45 Diagrama de flujo del funcionamiento del higrómetro digital fabricado.
(Realización y fuente: Propia)
68
2.7.4. Costo de fabricación del higrómetro digital con corrección
El costo de fabricación del higrómetro digital con corrección tiene un valor de $101,70
únicamente considerando materiales de la Tabla 2.19.
Tabla 2.19. Tabla de costos de fabricación del prototipo e higrómetro digital con corrección.
Sensor BMP180*
15,00
(ver Figura 2.4)
69
Alternativa
Alternativa 2 Alternativa 2 Alternativa Diseño de la
corrección Psicrometría
A HR corrección A carcasa
HR
Termómetros del
Termómetros del Alternativa 1 Alternativa 1 Alternativa 3 Programación
Laboratorio de Disposición
Psicrómetro de A HR H3 en la tarjeta Impresión de
Termodinámica definitiva de
Flujo continúo la carcasa
(TLT) elementos
(Ps FC)
(apreciación 1/20°C)
HR HR HR
Selección Selección de
Tarjeta sensores:
programable HR, Tbs y A
Psicrómetro de Higrómetros digitales Prototipo sin acondicionamiento Higrómetro digital con acondicionamiento
Flujo continuo comerciales
Ps FC
Figura 2.46 Metodología desarrollada en el proyecto.
(Realización y fuente: Propia) 70
3. RESULTADOS Y DISCUSIÓN
En esta sección se encuentran los resultados de todas las mediciones de las propiedades
psicrométricas del aire realizadas en diferentes localidades del Ecuador, con higrómetros
comerciales, con el higrómetro digital fabricado sin corrección y con corrección, así como
también su respectivo análisis. Además, se detalla el desarrollo de un programa para
reducir los errores humanos en mediciones de humedad relativa con cartas psicrométricas
y un estudio de condiciones ambientales entre Tulcán y Quito.
3.1. Resultados
3.1.1. Programa para el cálculo de humedad relativa
Para evitar que los valores reales de humedad relativa obtenidos con el psicrómetro de
flujo continuo presenten errores humanos (lectura de las cartas psicrométricas) se realiza
un programa con macros en Excel. Esta programación utiliza las ecuaciones del ASHRAE
Fundamentals Handbook [4], las cuales se observan en el ANEXO VII.
En la Tabla 3.1, Tabla 3.2 y Tabla 3.3, se observan los resultados de humedad relativa y
temperatura de bulbo seco del higrómetro (H1), higrómetro (H2), e higrómetro (H3)
respectivamente, respecto a los valores encontrados en las cartas psicrométricas a una
altitud de 2806 msnm (EPN-Quito). También, se visualiza los errores calculados entre
71
dichas propiedades. Cabe mencionar que los errores calculados en este trabajo son errores
porcentuales (ver sección 1.4, pp. 19-20).
Tabla 3.1. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital H1 y el
psicrómetro de flujo continuo a 2806 msnm (Quito).
Tabla 3.2. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital H2 y el
psicrómetro de flujo continuo a 2806 msnm (Quito).
72
Tabla 3.3. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital H3 y el
psicrómetro de flujo continuo a 2806 msnm (Quito).
Las mediciones de los tres higrómetros (H1, H2, H3) presentan errores considerables
únicamente en valores de humedad relativa (HR), mientras que en los valores de
temperatura de bulbo seco (TBS) se encuentra dentro de las características técnicas de
cada equipo. En la Figura 3.2 se observan los valores de humedad relativa de los tres
higrómetros digitales y del psicrómetro de flujo continuo a 2806 msnm.
Figura 3.2. Datos de humedad relativa de los higrómetros digitales (H1, H2 y H3), y psicrómetro de
flujo continuo (2806msnm valores reales).
(Fuente: propia)
73
En la Tabla 3.4 se observan los resultados de humedad relativa y temperatura de bulbo
seco del higrómetro digital del laboratorio de fluidos (H4), respecto a los valores
encontrados en las cartas psicrométricas a una altitud de 2806 msnm. De igual manera, se
visualiza los errores porcentuales calculados entre dichas propiedades.
Tabla 3.4. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital H4 y el
psicrómetro de flujo continuo a 2806 msnm (Quito).
Las mediciones de este higrómetro muestran errores en la humedad relativa (HR), mientras
que en la temperatura de bulbo seco (TBS) no tiene mayores complicaciones. En la Figura
3.3 se muestran los valores de humedad relativa del higrómetro digital (H4) y del
psicrómetro de flujo continuo.
Figura 3.3. Datos de humedad relativa del higrómetro digital H4 y del psicrómetro de flujo continuo
(2806msnm valores reales).
(Fuente: propia)
74
Los higrómetros digitales (H5 y H6) pertenecientes al LAEV, fueron analizados en el mismo
instante y bajo las mismas condiciones de temperatura y humedad relativa. En la Tabla 3.5
y en la Tabla 3.6 se observan los resultados de la humedad relativa (HR) y temperatura de
bulbo seco (TBS) de los higrómetros digitales, respecto a los valores tomados con el
psicrómetro de flujo continuo. Además, se visualiza el error porcentual calculado para
ambas propiedades.
Tabla 3.5. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital H5 y el
psicrómetro de flujo continuo a 2806 msnm (Quito).
Tabla 3.6. Datos de humedad relativa y temperatura de bulbo seco del higrómetro digital H6 y el
psicrómetro de flujo continuo a 2806 msnm (Quito).
Las mediciones de los dos higrómetros (H5 y H6) presentan un error considerable en los
valores de humedad relativa, mientras que en los valores de temperatura de bulbo seco se
encuentran dentro del margen de las especificaciones técnicas del fabricante. En la Figura
75
3.4 observa los datos de humedad relativa de ambos higrómetros digitales y del
psicrómetro de flujo continuo medidos a 2806 msnm (EPN-QUITO).
Figura 3.4. Datos de humedad relativa de los higrómetros digitales (H5 y H6), y psicrómetro de
flujo continuo (2806msnm valores reales).
(Fuente: propia)
Figura 3.5. Datos de humedad relativa de higrómetros digitales (H1, H2, H3, H4), y las cartas
psicrométricas (72 msnm valores reales).
(Fuente propia)
76
Tabla 3.7. Errores porcentuales de humedad relativa de los higrómetros digitales estudiados a
72msnm.
1 3% 3% 0% 2%
2 9% 9% 5% 9%
3 6% 6% 4% 8%
4 7% 7% 2% 8%
5 7% 7% 4% 7%
6 8% 8% 4% 7%
7 6% 6% 5% 5%
8 4% 4% 2% 4%
9 3% 3% 2% 2%
10 4% 4% 1% 3%
11 7% 7% 4% 6%
12 3% 3% 0% 2%
13 9% 9% 2% 7%
14 6% 6% 2% 6%
Promedio 6% 6% 3% 5%
(Fuente: propia)
El higrómetro digital sin corrección se evalúa en 3 localidades menores a los 1000 msnm.
Las mediciones a nivel del mar se las realiza en la ciudad de Atacames, los valores de las
propiedades psicrométricas anteriormente mencionadas se las observa en la Tabla 3.8. De
igual manera, se observa el error porcentual calculado entre los valores de humedad
relativa experimental y real (psicrómetro de flujo continuo).
77
Tabla 3.8. Valores de varias propiedades psicrométricas del aire con el psicrómetro de flujo continuo
y el higrómetro digital construido sin corrección en la ciudad de Atacames (nivel del mar).
MEDICIÓN TBS (°C) TBH (°C) HRR [%] HRSC [%] ERROR [%]
1 30 26,5 76,219 79 3,65
2 29,5 26 75,998 78 2,63
3 29,5 26 75,998 77 1,32
4 29,2 25,7 75,863 76 0,18
5 29,2 25,6 75,22 75 0,29
6 29 25,3 74,484 75 0,69
7 29,2 25,7 75,863 76 0,18
8 28,5 25 75,539 75 0,71
9 28,3 25 76,76 75 2,29
10 28,2 25 77,378 78 0,80
11 28,1 25 78,001 77 1,28
12 28,6 25,2 76,239 77 1,00
13 28,5 25,4 78,17 80 2,34
14 28,8 25,5 76,982 80 3,92
15 28,8 25,5 76,982 81 5,22
16 28,8 26 80,284 82 2,14
17 28,2 25 77,378 80 3,39
PROMEDIO 1,89
(Fuente: Propia)
Como se observa en la Tabla 3.8, las mediciones a nivel del mar no presentan mayores
complicaciones, ya que el error porcentual promedio en humedad relativa es de 1,89%.
Este valor se encuentra dentro de las especificaciones técnicas del fabricante. Los valores
de humedad relativa real y humedad relativa calculada sin corrección se observan en la
Figura 3.6.
Figura 3.6. Datos de humedad relativa real y experimental sin corrección en Atacames.
(Fuente: Propia)
78
En la Tabla 3.9 se observa los valores de las propiedades psicrométricas del aire
analizadas en la ciudad de Santo Domingo con el psicrómetro de flujo continuo y con el
higrómetro digital sin corrección.
Tabla 3.9. Valores de varias propiedades psicrométricas del aire con el psicrómetro de flujo continuo
y el higrómetro digital sin corrección en la ciudad de Santo Domingo (586 msnm).
MEDICIÓN TBS (°C) TBH (°C) HRR [%] HRSC [%] ERROR [%]
1 25 22 77,689 82 5,55
2 24,5 22 81,033 83 2,43
3 25 22 77,689 82 5,55
4 25 22 77,689 81 4,26
5 24,8 21,8 77,596 81 4,39
6 24,9 21,8 76,941 81 5,28
7 25 21,8 76,291 81 6,17
8 24,9 22 78,347 83 5,94
9 24,8 22 79,01 86 8,85
10 24,4 22 81,719 87 6,46
11 24,2 22 83,107 88 5,89
12 24 22 84,518 87 2,94
13 24 21,5 80,831 86 6,39
14 24 21,8 83,035 87 4,78
15 24,1 22 83,81 85 1,42
16 24,1 21,9 83,071 87 4,73
17 24 22 84,519 87 2,94
PROMEDIO 4,94
(Fuente: Propia)
Figura 3.7. Datos de humedad relativa real y experimental sin acondicionamiento de señal en
Santo Domingo.
(Fuente: Propia)
79
En la Tabla 3.10 se observa los valores de las propiedades psicrométricas del aire
analizadas en el recinto de Rocafuerte (Esmeraldas) con el psicrómetro de flujo continuo y
con el higrómetro digital sin corrección.
Tabla 3.10. Valores de varias propiedades psicrométricas del aire con el psicrómetro de flujo
continuo y el higrómetro digital sin corrección en la ciudad de Rocafuerte (950 msnm).
MEDICIÓN TBS (°C) TBH (°C) HRR [%] HRSC [%] ERROR [%]
1 24 22 84,71 88 3,88
2 24 22 84,71 88 3,88
3 23,7 21,7 84,61 89 5,19
4 23,5 21,7 86,04 90 4,60
5 23,4 21,7 86,76 90 3,73
6 23,4 21,8 87,52 91 3,98
7 23,2 21,3 85,2 91 6,81
8 23,1 21,3 85,92 91 5,91
9 23,8 21,8 84,65 86 1,59
10 23,3 21,6 86,74 90 3,76
11 23,4 21,4 84,52 90 6,48
12 23,2 21,3 85,2 90 5,63
13 23,5 21,7 86 88 2,33
14 23,2 21,3 85,2 90 5,63
PROMEDIO 4,53
(Fuente: Propia)
Figura 3.8. Datos de humedad relativa real y experimental sin acondicionamiento en Rocafuerte.
(Fuente: Propia)
80
Los valores que se obtienen mediante el prototipo de higrómetro digital sin corrección en
las localidades menores a 1000 msnm presentan errores relativamente bajos (dentro de
las especificaciones técnicas de los proponentes del proyecto). Por lo cual, se considera
anular la corrección en humedad relativa para estas altitudes.
El higrómetro digital sin corrección se avalúa en localidades con altitudes superiores a los
2000 msnm. Las propiedades psicrométricas del aire analizadas son la temperatura de
bulbo seco (TBS), temperatura de bulbo húmedo (TBH), humedad relativa real (HRR) con
el psicrómetro de flujo continuo y la humedad relativa sin corrección del higrómetro digital
(HRSC), de igual manera se calcula el porcentaje de error entre estas humedades.
Tabla 3.11. Valores de varias propiedades psicrométricas del aire con el psicrómetro de flujo
continuo y el higrómetro digital sin corrección en la ciudad de Quito Centro (2804 msnm).
MEDICIÓN TBS (°C) TBH (°C) HRR [%] HRSC [%] ERROR [%]
1 18 11 47,926 40 16,54
2 18,9 11,2 44,621 40 10,36
3 19,5 11,9 45,97 39 15,16
4 21,5 12,8 41,969 38 9,46
5 21,2 12,7 42,73 37 13,41
6 21,9 13 41,346 37 10,51
7 21,5 12,2 38,66 34 12,05
8 21,9 12,2 37,024 33 10,87
9 22,2 13,5 42,799 35 18,22
10 23 13,2 37,896 35 7,64
11 23 13,3 38,416 34 11,50
12 23,1 13,3 38,019 34 10,57
13 23 13,2 37,896 33 12,92
14 23 13 36,865 32 13,20
15 22,8 12,8 36,613 30 18,06
16 22,8 13,2 38,691 31 19,88
PROMEDIO 13,15
(Fuente: Propia)
El error porcentual promedio de humedad relativa medido a esta altitud (2804 msnm) es de
13,15%, este valor se encuentra fuera de las especificaciones técnicas del fabricante. La
variación de la humedad relativa real respecto a la experimental sin corrección se la
observa en la Figura 3.9.
81
Figura 3.9. Datos de humedad relativa real y experimental sin corrección en Quito Centro.
(Fuente: Propia)
En la Tabla 3.12 se observa los valores de las propiedades del aire analizadas con el
psicrómetro de flujo continuo y con el higrómetro digital sin corrección en una altitud de
2951 msnm (Quito-SUR).
Tabla 3.12. Valores de varias propiedades psicrométricas del aire con el psicrómetro de flujo
continuo y el higrómetro digital sin corrección en la ciudad de Quito Sur (2951 msnm).
MEDICIÓN TBS (°C) TBH (°C) HRR [%] HRSC [%] ERROR [%]
82
El error porcentual promedio de humedad relativa medido a esta altitud (2951 msnm) es de
15,87%, este valor se encuentra fuera de las especificaciones técnicas del fabricante. La
variación de la humedad relativa real respecto a la experimental sin corrección se la
observa en la Figura 3.10.
Figura 3.10. Datos de humedad relativa real y experimental sin corrección en Quito Sur.
(Fuente: Propia)
En la Tabla 3.13 se observa los valores de las propiedades del aire analizadas con el
psicrómetro de flujo continuo y con el higrómetro digital sin corrección en una altitud de
3046 msnm (Tulcán).
Tabla 3.13. Valores de varias propiedades psicrométricas del aire con el psicrómetro de flujo
continuo y el higrómetro digital sin corrección en la ciudad de Tulcán (3046 msnm).
ERROR
Medición TBS (°C) TBH (°C) HRR [%] HRSC [%]
(%)
1 15,3 12,8 78,53 72 8,32%
2 15 12,5 78,26 73 6,72%
3 14,9 12,2 76,55 73 4,64%
4 14,8 12,5 79,82 75 6,04%
5 15 12,8 80,76 74 8,37%
6 15 12,8 80,76 78 3,42%
7 15,1 12,7 79,14 78 1,44%
8 14,7 12,5 80,61 77 4,48%
9 14,8 12,7 81,51 80 1,85%
10 15,1 12,9 80,81 80 1,00%
11 15,5 13 78,53 78 0,67%
12 15,3 12,9 79,24 80 0,96%
13 15,8 13,3 78,68 83 5,49%
14 16 13,8 81,22 83 2,19%
15 16,1 14 82,08 81 1,32%
16 15,8 13,1 77,07 80 3,80%
17 16 13,8 81,22 83 2,19%
PROMEDIO 3,70%
(Fuente: Propia)
83
El error porcentual promedio de humedad relativa medido a 3046 msnm es de 3,70%, este
valor se encuentra dentro de las especificaciones técnicas del fabricante. La variación de
la humedad relativa real respecto a la experimental sin corrección se la observa en la Figura
3.11.
Figura 3.11. Datos de humedad relativa real y experimental sin corrección en Tulcán.
(Fuente: Propia)
Tabla 3.14. Datos de altitudes de referencia y altitudes obtenidas con el higrómetro digital sin
corrección en diferentes localidades del Ecuador.
(Fuente: propia)
84
3.1.4. Pruebas del higrómetro digital con acondicionamiento de señal
En esta sección se muestran los resultados de las mediciones de humedad relativa con el
higrómetro digital corregido, estas mediciones se realizan en diferentes altitudes.
Las propiedades analizadas son la temperatura de bulbo seco (TBS), temperatura de bulbo
húmedo (TBH), humedad relativa real (HRR) con el psicrómetro de flujo continuo, y la
humedad relativa con corrección del higrómetro digital (HRCC). Además, se obtiene el
porcentaje de error entre las humedades relativas.
El higrómetro digital con corrección se analiza en altitudes superiores a los 2000 msnm. En
la Tabla 3.15 se observa los valores de las propiedades psicrométricas analizadas a una
altitud de 2804 msnm (Quito-Centro).
Tabla 3.15. Valores de varias propiedades psicrométricas del aire con el psicrómetro de flujo
continuo y el higrómetro digital con corrección en la ciudad Quito (2804 msnm).
MEDICIÓN TBS [°C] TBH [°C] HRR [%] HRCC [%] ERROR [%]
1 28,2 17,5 39,39 39,78 0,99
2 27,9 17 38,18 36,89 3,38
3 27,2 16,5 38,34 36,35 5,19
4 25,2 15 38,5 36,34 5,61
5 24,2 14,8 41,37 40,35 2,47
6 23,2 13,9 40,74 38,53 5,42
7 23 13,8 41,04 37,8 7,89
8 22 13,5 43,67 40,06 8,27
9 21,7 13,2 43,33 42,2 2,61
10 21,2 13,2 45,59 42,63 6,49
11 20,8 12,9 45,71 42,39 7,26
12 20,1 12,7 47,88 45,87 4,20
13 20 13 50,22 47,83 4,76
PROMEDIO 4,96
(Fuente: Propia)
En la Figura 3.12 se observa la variación que existe entre la humedad relativa real y la
humedad relativa experimental con el higrómetro digital corregido.
85
Figura 3.12. Datos de humedad relativa real y experimental con corrección en Quito Centro.
(Fuente: Propia)
En la Tabla 3.16 se observa los valores de las propiedades del aire analizadas con el
psicrómetro de flujo continuo y con el higrómetro digital corregido en una altitud de 3046
msnm (Tulcán).
Tabla 3.16. Valores de varias propiedades psicrométricas del aire con el psicrómetro de flujo
continuo y el higrómetro digital corregido en la ciudad de Tulcán (2804 msnm).
MEDICIÓN TBS [°C] TBH [°C] HRR [%] HRCC [%] ERROR [%]
86
Como se observa en la Tabla 3.16 el error porcentual promedio en humedad relativa
calculado en la ciudad de Tulcán (3046 msnm) es de 1,86%. Este valor es bajo comparado
con los resultados de la Tabla 3.13 el cual es 3,70%.
En la Figura 3.13 se observa la variación que existe entre la humedad relativa real y la
humedad relativa experimental con el higrómetro digital con corrección de señal.
Figura 3.13. Datos de humedad relativa real y experimental con corrección en Tulcán.
(Fuente: Propia)
Tabla 3.17. Datos de altitudes de referencia y altitudes obtenidas con el higrómetro digital con
corrección en diferentes localidades del Ecuador.
(Fuente: Propia
En la Tabla 3.18 se encuentra los resultados finales del porcentaje de error porcentual
promedio de humedad relativa utilizando las tres ecuaciones (4), (9) y (16).
87
Tabla 3.18. Resultados del cálculo de error porcentual promedio de humedad relativa con las tres
ecuaciones de corrección.
Mínimos
Matlab Excel
Cuadrados
Error en Humedad
3,14% 4,06% 2,87%
relativa
(Fuente: Propia)
3.2. Discusión
3.2.1. Datos obtenidos con el higrómetro digital sin corrección
Para obtener los valores de referencia en altitud se utilizan los datos satelitales de Google
Earth en las localidades analizadas. Con el fin de garantizar que estos datos sean
adecuados para la corrección del higrómetro digital se realiza una verificación en donde se
comparan las altitudes de varias localidades del Ecuador obtenidas de Google Earth con
los datos satelitales del GIS (Geographic Information System). De esta verificación se
obtiene que Google Earth presenta un error porcentual promedio de 0,61%. Los valores de
este análisis se observan en la Tabla 2.5, p. 38.
La altitud experimental de los lugares analizados presenta varios errores al comparar con
los valores de Google Earth (altitud de referencia), donde se obtiene un error porcentual
promedio de 15,4%, tal como se observa en la Tabla 3.14, p. 84.
Las mediciones de varias propiedades psicrométricas del aire se las realiza en diferentes
localidades del Ecuador, las cuales van desde los 0 msnm hasta los 3040 msnm. Los
resultados de estas mediciones en localidades menores a los 1000 msnm se observan en
la Tabla 3.8, la Tabla 3.9 y la Tabla 3.10 (ver pp. 78-80). En las cuales se visualiza que la
temperatura de bulbo seco y la humedad relativa de cada lugar se encuentra dentro de las
especificaciones técnicas de los fabricantes del higrómetro digital.
La diferencia entre los errores promedios de humedad relativa de las localidades menores
a los 1000 msnm está directamente relacionada con la temperatura de bulbo húmedo. Esta
temperatura en la ciudad de Atacames se encuentra entre 25 y 26 °C, en cambio para las
88
localidades de Santo Domingo y Rocafuerte el rango se encuentra entre 21 y 22 °C. Esto
se justifica con el análisis de sensibilidad que se realiza en la sección 2.4.3, pp. 42-45.
Los resultados de las propiedades psicrométricas del aire medidos con el higrómetro digital
sin corrección en altitudes superiores a los 2000 msnm se observan en la Tabla 3.11, la
Tabla 3.12 y la Tabla 3.13 (ver pp. 81-83). Estas mediciones muestran un elevado error
porcentual promedio en los valores de humedad relativa y se encuentran fuera las
especificaciones técnicas del higrómetro digital.
El error porcentual promedio que presenta la humedad relativa en la ciudad de Quito centro
(2804 msnm) es de 13,15% y para Quito sur (2951 msnm) es de 15,87%. Estos valores
son elevados en comparación con el error porcentual promedio que se obtiene en la ciudad
de Tulcán (3046 msnm) que es de 3,70%, a pesar de que la cuidad de Tulcán se encuentra
a una altitud más elevada a la de Quito. Esto se debe a que las condiciones ambientales
de estas ciudades son diferentes, con esto se evidencia que la altitud no es la única
causante de errores de humedad relativa, sino también el ambiente mismo del lugar. Para
asegurar que esta aseveración es correcta se realiza un análisis variando las condiciones
ambientales de Tulcán con resistencias eléctricas, estos resultados se observan en la
sección 3.2.2.
3.2.2. Datos de humedad relativa del aire analizadas en Tulcán (3040 msnm)
Para justificar esta variación se grafica la humedad relativa de tres años consecutivos de
estas dos localidades como se observa en la Figura 3.14. Los datos son obtenidos de los
anuarios meteorológicos del INAMHI específicos para cada año. Como se observa en la
Figura 3.14, la diferencia promedio en humedad relativa que se encuentra entre estas dos
localidades varía entre 14% y 18%, siendo los valores en Tulcán más elevados durante
todo el año. Además, los valores de humedad relativa en Tulcán se encuentran fuera del
rango de confort humano (40-70%) según ASRHAE.
89
Figura 3.14. Datos mensuales de humedad relativa para la ciudad de Quito y Tulcán 2012-2014.
(Fuente: [17], realización: Propia)
Con el fin de comparar los errores de humedad relativa en Tulcán se realiza un análisis
variando la cantidad de humedad y temperatura mediante resistencia eléctricas, hasta
alcanzar un ambiente similar al de la ciudad de Quito. En la Tabla 3.19 se observa los datos
que se obtienen con este análisis.
Tabla 3.19. Valores de varias propiedades psicrométricas del aire con el psicrómetro de flujo
continuo y el higrómetro digital sin corrección la ciudad de Tulcán (3046 msnm).
% HR T °C
HR valores ERROR
Medición TBS (°C) TBH (°C) Psicrómetro Psicrómetro
reales [%] (%)
1 1
1 23 14 42,53 50 20,4 17,56%
2 22,5 13,8 43,57 49 20,6 12,46%
3 23 14,1 43,06 50 20,7 16,12%
4 23,5 14,5 43,08 49 21,2 13,74%
5 22,9 14,3 44,55 51 20,7 14,48%
6 22,5 13 44,66 52 20,5 16,44%
7 22,5 14,5 47,41 52 20,5 9,68%
8 22,8 14,5 46,07 51 20,4 10,70%
9 22,5 14 44,66 52 20,2 16,44%
10 23 15 47,93 54 20,3 12,66%
11 23 14,1 43,06 52 20,5 20,76%
12 24 15 43,63 50 21,3 14,60%
13 24,5 16 46,74 52 22 11,25%
14 25 16 44,66 51 22,4 14,20%
15 22 15 52,62 63 18,7 19,73%
16 19 14 62,54 70 17,7 11,93%
PROMEDIO 14,55%
(Fuente: Propia)
90
Al comparar los datos de la Tabla 3.13, p. 83, con los datos de la Tabla 3.19 de la ciudad
de Tulcán, se observa un gran cambio, el error porcentual promedio en humedad relativa
del higrómetro digital sin corrección presenta un incremento del 3,7% al 14,55%. Con esto
se verifica que el incremento de error de humedad relativa no solo depende de la altitud,
sino también del ambiente en el que se encuentre.
Los valores de humedad relativa real y experimental mostrados en la Figura 2.19, p.41,
muestran una distribución de los datos de manera desordenada. La ecuación que
representa esta curva tiene un coeficiente de determinación de 0,9703, este resultado
indica que si es aceptable esta ecuación de corrección. Pero después de demostrar que la
altitud ya no es la única causante del error de humedad relativa mostrado en la sección
3.2.2 y analizando los valores del análisis de sensibilidad de la humedad relativa (ver
sección 2.4.3 pp. 42-45) se procede a corregir los valores de temperatura de bulbo húmedo.
Para encontrar los valores de temperatura de bulbo húmedo a partir del sensor DHT22 se
realiza un programa en macros Excel en cual se ingresa el valor de humedad relativa y la
temperatura ambiente (ver Figura 2.20, p. 41). Con los datos obtenidos de temperaturas
de bulbo húmedo experimentales y reales se realiza una gráfica y se encuentra que los
datos tienen una distribución lineal (ver Figura 2.21 p. 42). Esta gráfica presenta un
coeficiente de determinación de 0,995 el cual es mayor que el obtenido con humedades
relativas.
Para encontrar la ecuación de corrección, se analiza los datos por tres métodos diferentes,
obteniendo tres ecuaciones. Para seleccionar una de las tres ecuaciones, se consideran
dos valores, el primero es el coeficiente de determinación (el más cercano a uno es la
ecuación más eficiente) y el segundo es el cálculo del error en humedad relativa al corregir
los valores obtenidos sin corrección. Este último valor se lo observa en la Tabla 3.18, p.
88.
La ecuación más eficiente con estos resultados es la que se obtiene mediante el Excel
(ecuación 4), para corroborar dicha eficiencia se realizaron mediciones en Quito y en
Tulcán con la ecuación de corrección programada. Los resultados con la ecuación de Excel
no fueron óptimos ya que el porcentaje de error en humedad relativa aumentó. Por lo tanto,
se utiliza la ecuación (16) obtenida mediante Matlab para verificar y comparar ambos
resultados. Los resultados de las dos programaciones permiten corroborar que el uso de
91
la ecuación de Excel no se puede utilizar en localidades de alta humedad, los resultados
de la ecuación de Matlab se utilizan para bajas y altas humedades, obteniendo un resultado
aceptable.
Para la corrección de altitud se utiliza dos ecuaciones, la ecuación (2) se utiliza para
corregir la altitud de lugares superiores a los 1000 msnm, mientras que la ecuación (3) se
utiliza para lugares menores a los 1000 msnm. Mediante estas dos ecuaciones se obtiene
una reducción en el error porcentual promedio de altitud de 15,4% (ver Tabla 3.14) a 0,3%
(ver Tabla 3.17) para localidades superiores a los 2000 msnm.
92
4. CONCLUSIONES Y RECOMENDACIONES
4.1. Conclusiones
En el presente trabajo se cumple con el objetivo de fabricar un higrómetro digital con
prestaciones similares a los disponibles en el mercado y con corrección en humedad
relativa y altitud.
La corrección más apropiada resulta ser a la temperatura de bulbo húmedo a través de una
ecuación totalmente independiente de la altitud (ver análisis de la sección 2.4.2, pp. 40-42,
sección 2.4.3, pp. 42-45 y sección 2.4.4, pp. 45-52). De esta corrección se deriva luego el
valor de humedad relativa final.
Los resultados del coeficiente de determinación entre el Excel, regresión por mínimos
cuadrados (calculo manual) y Matlab, son muy cercanos a la unidad, es decir las
correcciones al usar estas ecuaciones son confiables. A partir del análisis realizado en la
sección 3.2.3, pp. 91-92, se concluye que la ecuación de corrección que brinda los mejores
resultados para el cálculo de humedad relativa es la ecuación (16) obtenida mediante
Matlab.
93
Tabla 3.13, p. 83). Cabe señalar que la humedad relativa en esta localidad se encuentra
sobre el 70% durante todo el año, ver Figura 3.14, p. 90.
4.2. Recomendaciones
Para ampliar el rango de medición del higrómetro fabricado y de otros higrómetros digitales
se requiere fabricar un túnel con acondicionamiento de aire que sea portátil y que permita
modificar la temperatura y humedad relativa del aire. De tal manera que se varíe la altitud
llevándolo a diferentes localidades del Ecuador.
Para reducir las restricciones de trabajo del higrómetro fabricado se requiere utilizar
materiales con altas propiedades térmicas para la fabricación de la carcasa. Esto debido a
que el sensor de humedad relativa trabaja en rangos de -40 °C a 80 °C y el material utilizado
en la impresión 3D de la carcasa (PLA) no puede sobrepasar los 50 °C limitando así el
rango del equipo.
El sensor de humedad y temperatura del higrómetro digital fabricado debe estar expuesto
a flujos de aire, por lo que se recomienda no ponerlo en contacto directo e indirecto con
agua. También no se debe presionar los tres pulsadores al mismo tiempo ya que se
generan cortos dentro del equipo.
94
BIBLIOGRAFÍA
[2] D. P. Gatley, Understanding Psychrometrics, Third ed., Atlanta: Tracy Becker, 2013.
[12] C. R. P. Chapra Steven C., Métodos Numéricos para Ingenieros, Quinta ed., México:
McGraw-Hill, 2006.
[13] ASTM, «Standard Test Method for Measuring Humididy with a Psychrometer (the
Measurement of Wet-and Dry-Bulb Temperatures),» ASTM E337-15, 2018.
[15] G. F. Jesús Belchí, Excel 2000, Herramientas avanzadas, macros y programación con
VBA, FC, 2000.
[16] Tresdpro, «¿Qué material utilizan las impresoras 3d?,» 10 octubre 2019. [En línea].
Available: [Link] [Último
acceso: 15 enero 2020].
95
[17] INAMHI, «Servicios Meteorológicos,» 2 Enero 2020. [En línea]. Available:
[Link]/. [Último acceso: 5 Febrero 2020].
96
ANEXOS
97
ANEXO II. CONDICIONES DE OPERACIÓN DEL SENSOR BMP180.
98
99
ANEXO III. CONDICIONES DE OPERACIÓN DEL LCD 1602.
100
ANEXO IV. CONDICIONES DE OPERACIÓN Y ESQUEMA DEL
MICROCONTROLADOR ATMEGA328P DEL ARDUINO UNO.
101
ANEXO V. ECUACIONES PARA EL DESARROLLO DEL PROGRAMA
“CÁLCULO DE TEMPERATURA DE BULBO HÚMEDO”
Número
Nombre de la ecuación Ecuación de
ecuación
21
35
23
20
102
ANEXO VI. CONDICIONES DE OPERACIÓN DEL LCD NOKIA 5110.
103
104
105
ANEXO VII. ECUACIONES PARA EL DESARROLLO DEL PROGRAMA
“CÁLCULO DE HUMEDAD RELATIVA”
Número
Nombre de la ecuación Ecuación de
ecuación
21
35
23
106
24,94 1
3
19
30
2
36
179,12
4
6
R12,50