INTRODUCCIÓN
La estadística experimental -en el mundo moderno- es una necesidad real, presente en todas las áreas del conocimiento humano, como una herramienta para auxiliar las decisiones a ser tomadas, permitiendo mayor seguridad a la investigación a ser transformada en tecnología y utilizada por la humanidad o por parte de ella (López & González, 2014).
En la actualidad, los investigadores realizan experimentos virtualmente en todos los campos del saber, por lo general, para descubrir algo acerca de un proceso o sistema en particular; literalmente, un experimento es una prueba o ensayo. Un experimento diseñado es una prueba o serie de pruebas, en las cuales, se inducen cambios deliberados en las variables de entrada de un proceso o sistema, de manera que sea posible observar e identificar las causas de los cambios en la respuesta de salida.
El proceso o sistema bajo estudio, se puede representar por medio del modelo, que permita visualizar el proceso, como una combinación de máquinas, métodos, personas y otros recursos, que transforman alguna entrada (a menudo, un material) en una salida, que tiene una o más respuestas observables. Mediante un diseño adecuado del experimento, se pretende obtener la máxima información de un proceso de la forma más rápida, económica y simple posible.
En sentido general, cuando se habla de diseñar un experimento se debe entender todo aquello relacionado con la planificación, lo más racional posible del mismo, con vistas a su ejecución y un posterior análisis adecuado de sus resultados, que conduzcan a deducciones aceptables del problema (Montgomery, 1991).
Teniendo en cuenta la importancia del uso de las herramientas estadísticas para el análisis de resultados experimentales, se pretende realizar un estudio sobre la velocidad de corrosión del acero, puesto que existen multitudes de construcciones y de equipos metálicos expuestos a la atmósfera, tales como: marcos de ventanas, verjas, postes de alumbrado, maquinaria agrícola, tanques de almacenamiento, puentes, complejas instalaciones fabriles, entre otras estructuras. Se calcula, que un país presenta pérdidas entre un 2 y 5% del PIB, debido a la corrosión (Cadena et al. 2014).
La corrosión puede ser definida como una reacción irreversible de un material con el medio ambiente, lo que resulta en la degradación del material y la pérdida de sus propiedades (Cabrera, 2016).
Son variadas las formas existentes de expresar la velocidad de los procesos corrosivos, debido a la diversidad de formas en que se presenta la corrosión, lo que obliga a buscar la manera más adecuada de reflejar la intensidad del ataque. En este caso, se determina el índice gravimétrico negativo, más comúnmente conocido, como pérdida de peso (DP) (Domínguez & Castro, 1987).
Se debe tener en cuenta que Cuba presenta un clima caracterizado por tener más de la mitad del año temperaturas medias superiores a 25ºC y humedad relativa media de alrededor del 80%. Dada la configuración y la ubicación geográfica del país, la influencia del aerosol marino llega a casi todo el territorio nacional. Estas condiciones favorecen notablemente el deterioro de diferentes materiales, en especial, los metálicos, específicamente, en zonas de agresividad corrosiva alta y extrema. La tendencia actual de los cambios climáticos hace que la atmósfera se torne cada vez más agresiva en Cuba, en cuanto a las penetraciones del mar y el incremento en la fuerza de los vientos y las lluvias (Castañeda et al. 2012).
El Laboratorio de Ensayos de Tropicalización (LABET), se encuentra situado al este de La Habana, en la zona costera de Cojímar, que se caracteriza por poseer uno de los microclimas de mayor agresividad de la República de Cuba y del mundo, dado los altos valores de salinidad, humedad relativa, radiación solar y temperatura (Marrero & González, 2005).
Los ensayos climáticos naturales necesitan tiempos excesivamente largos de exposición, a menudo, varios años, lo que justifica el uso de métodos rápidos de valoración. La utilización de cámaras climáticas de niebla salina neutra con una salinidad de 5% y temperatura de 35ºC, permite someter las muestras, a la acción de una niebla de sal, atomizada en condiciones prefijadas de precisión, de concentración y de temperatura (NC ISO 9227, 2014).
En este trabajo, se evalúa el comportamiento de la velocidad de corrosión del acero, a partir de software estadísticos y modelos matemáticos. Se estima y se compara los efectos de los tratamientos asignados a la variable independiente sobre la variable rendimiento, sin los efectos de los bloques.
MATERIALES Y MÉTODOS
El desarrollo experimental de este trabajo, se efectuó en el Laboratorio de Ensayos de Corrosión del LABET. En este capítulo, se describen los materiales, la metodología seguida en el trabajo experimental y las condiciones en las que se realizaron las pruebas de los sistemas estudiados.
Diseño de experimentos: El diseño de los experimentos realizados en este trabajo pretende evaluar el comportamiento del factor grado de oxidación de las probetas de acero sobre la velocidad de la corrosión. Para ello, se ejecutó la planificación de la actividad de la forma más racional posible, con vistas a su ejecución y un posterior análisis adecuado de sus resultados, que conduzcan a deducciones aceptables del problema.
El trabajo permitirá el uso del diseño en bloques al azar, puesto que distribuye las unidades experimentales en grupos o bloques, de tal forma que, cada una de las unidades experimentales que lo conforman, sean lo más homogéneas posibles y su número igual a la cantidad de tratamientos que se comparan, distribuyéndose estos tratamientos al azar a las unidades experimentales, dentro de cada bloque.
El diseño en bloques completos al azar toma en cuenta los tres principios básicos de la experimentación: repetición, aleatorización y control local. En este diseño, las unidades experimentales se distribuyen en grupos homogéneos. Cada uno de estos grupos es llamado bloque, siendo el tiempo de exposición en la cámara de niebla salina neutra. El número de unidades experimentales dentro de cada bloque es igual al número de tratamientos incluidos en el experimento. Un caso particular de diseño de bloques es el que aparece relacionado con la prueba de t de student, para muestras pareadas, aunque el número de tratamientos es solo dos. Los tratamientos son distribuidos en las unidades experimentales dentro de cada bloque aleatoriamente; así, cada bloque irá a constituir una repetición (Montgomery, 1991).
Variables de investigación:
Variable dependiente: velocidad de corrosión
Donde:
pi y pf: Peso inicial y final de la muestra metálica. (g)
S: Área de muestra expuesta al medio corrosivo. (m2)
t: Tiempo de exposición de la misma al medio corrosivo. (h)
Factores:
-Tiempo de exposición (h): Variable bloque 250 h, 500 horas, 750 h, 1000 h.
-Grado de oxidación de la probeta de acero (Grado A, B, C, D)
Grado A: Superficie de acero con la chapa de laminación intacta en toda la superficie y prácticamente sin corrosión.
Grado B: Superficie de acero con principio de corrosión y, de la cual, la chapa de laminación solo se encuentra parcialmente adherida.
Grado C: Superficie de acero, donde la chapa de laminación se ha perdido por efecto de la corrosión o es fácilmente eliminable por raspado, al encontrarse suelta. La corrosión es generalizada, pero no se han formado aún cavidades visibles.
Grado D: Superficie de acero con corrosión generalizada, exenta de chapa de laminación y gran cantidad de cavidades profundas (NORMA TÉCNICA ECUATORIANA. 2014, ISO 8501-1, 2007).
Teniendo en cuenta los factores y la variable dependiente considerados, se diseñó un plan experimental, con base en bloques al azar, lo que redujo el número total de probetas de ensayos.
Durante la realización del experimento, la temperatura, la salinidad, el pH y la humedad son consideradas constantes, puesto que la cámara está programada para no variar en ningún momento los valores de estos parámetros. El laboratorio se encuentra climatizado, lo que garantiza la temperatura de 23 ± 2°C y humedad relativa de 50 ± 5%.
Descripción de la prueba: La metodología a seguir en la prueba de corrosión, se presenta en el siguiente diagrama, que muestra los pasos a seguir para la preparación, la limpieza y la evaluación de la velocidad de corrosión, por pérdida de masa del acero.
Pasos a seguir:
Selección de las probetas
Preparación de las probetas según el grado de oxidación
Determinación de dimensiones de las probetas
Pesado inicial de las probetas
Exposición de las probetas en cámara de niebla salina neutra
Tratamiento de las probetas posterior a las pruebas
Pesado final de las probetas
Cálculo de la velocidad de la corrosión
Preparación de las probetas según el grado de oxidación: Para la preparación de las probetas de ensayo, se ha partido de chapas de acero de construcción, laminado en caliente y de bajo contenido de carbono, contenido de cobre de 0,03 a 0,10% y de fósforo menor que 0,07%. Las placas son rectangulares, con dimensiones de 100 x 150 mm y espesor de 1mm.
Las probetas fueron clasificadas (NORMA TÉCNICA ECUATORIANA, 2014, según ISO 8501-1, 2007) teniendo en cuenta el grado de oxidación (la cascarilla de laminación y el óxido, que se encuentran comúnmente en las superficies de acero, erigido sin recubrimiento y el acero en almacenamiento) (Figura 1).
Antes de la exposición en la cámara de ensayo, las probetas fueron desengrasadas con disolvente orgánico. No presentaban manchas de óxido visible o productos de corrosión en su superficie.
Determinación de las dimensiones de las probetas: Las probetas representan un ortoedro o paralelepípedo, puesto que todas sus caras son rectángulos y perpendiculares entre sí, por lo tanto, la siguiente ecuación permite determinar el área total de las mismas. Ancho: a (a) = 0,1 m, largo: L(b) = 0,15 m y espesor: e(c) = 0,001m.
Atotal = 2(ab+ac+bc) (Ecuación 2)
Atotal = 2(0,1*0,15+0,1*0,001+0,15*0,001)
Atotal = 0.0305 m2
Tratamiento de las probetas posterior a las pruebas: Al finalizar el período de ensayo, se sacaron las probetas de la cámara de niebla salina neutra y se dejaron secar por un tiempo, entre 0,5 h y 1 h antes de enjuagar, para reducir el riesgo de eliminar los productos de corrosión; después fueron examinadas, removiendo cuidadosamente los residuos de la solución de rociado de sus superficies. Las muestras se enjuagaron con agua corriente limpia, a una temperatura de 23 ± 5°C.
Teniendo en cuenta el tiempo de exposición establecido en el diseño de experimentos, los productos de corrosión formados en las muestras fueron eliminados, de acuerdo con el procedimiento de limpieza química y se determinó el peso final de las probetas.
Descripción del equipamiento: La velocidad de corrosión en las probetas de acero, se determinó por el método gravimétrico de pérdida de peso, empleando la cámara climática de niebla salina neutra, que permite someter las muestras a la acción de una niebla de sal atomizada, en condiciones prefijadas de precisión, de concentración y de temperatura.
La niebla no incide directamente sobre los paneles (que se encuentran colocados en posición vertical o con un ángulo entre 10° y 30° con respecto a la vertical), sino que es recibida por un deflector acrílico. El condensado es reciclado, no así el que escurre de los paneles, que es eliminado por la parte inferior de la cámara.
La cámara climática presenta una capacidad de 1,1m3, empleándose una solución de cloruro de sodio p.a. 50±5g/L, con pH entre 6,5 y 7,2, densidad de la solución entre 1,0255 y 1,0400g/cm3 y una temperatura en el interior de la cámara de 35 ± 2°C (NC ISO 9227, 2014). Las mismas, se ubicaron en portamuestras resistentes a la acción de la solución de rociado, lo que descarta la influencia de la presencia de los portamuestras en los ensayos.
Materiales: Los materiales y equipos empleados en el desarrollo experimental de las pruebas de corrosión, se listan a continuación:
- Cámara de niebla salina neutra
- Balanza analítica digital
- Desecador
- Portamuestras
- Agitadores
- Vasos de precipitado de 500 mL
- Matraz Erlenmeyer de 1000 mL
- Reactivos químicos para limpieza y desengrase
Reactivos: Según la norma (NC ISO 9226, 2015), se realiza el proceso de decapado para el acero por inmersión durante 10 min, con una temperatura entre 20 a 25°C, utilizando las siguientes proporciones de los reactivos: 500mL de ácido clorhídrico (HCl, ρ=1,19 g/mL), 35g de hexametilentetramina y agua destilada para completar a 1000mL.
RESULTADOS Y DISCUSIÓN
Se introducen los datos en el STATGRAPHICS, según lo logrado en todas las combinaciones posibles.
Resumen Multifactorial ANOVA - Velocidad de Corrosión
- Variable dependiente: velocidad de Corrosión
- Factores:
Grado de corrosión: Factor principal
Tiempo de exposición en cámara climática: Variable bloque
- Número de ensayos: 16
El análisis unifactorial de la varianza para la variable dependiente (velocidad de corrosión) permitirá:
- Realizar varias pruebas y gráficas, con el fin de determinar qué factores tienen un efecto estadísticamente significativo en la velocidad de corrosión.
- Las pruebas F de Fisher de la tabla ANOVA permitirá analizar si los diferentes niveles del factor A producen el mismo efecto en la característica observada o rendimiento en las diferentes unidades experimentales.
- El gráfico de medias permitirá la interpretación de los efectos significativos.
- Los gráficos de residuales ayudarán a juzgar si el método matemático es el idóneo.
Análisis de varianza para la Velocidad de Corrosión - Tipo III Suma de cuadrados: La tabla 1 de ANOVA, descompone la variabilidad de la velocidad de corrosión en función de los factores evaluados. Las sumas de cuadrados Tipo III (el valor por defecto) mide la influencia de cada factor y elimina los efectos del resto. Los valores Pvalue indican la significación estadística de cada uno de los factores.
Este experimento evidencia que los factores (grado de corrosión y de tiempo de exposición en la cámara) poseen los valores de Pvalue inferiores a 0,05, por lo que se rechaza la hipótesis nula, es decir, se acepta la alternativa. Los factores tienen un efecto estadísticamente significativo en la velocidad de corrosión, al nivel de confianza de 95,0%.
Prueba de Rango Múltiple para la Velocidad de Corrosión (Grado de Corrosión). Método LSD: 95,0%: La tabla 2 evidencia el procedimiento de comparación múltiple, para determinar las medias significativamente diferentes. En la parte superior de la misma, se identifican 2 grupos homogéneos, usando columnas de X. Dentro de cada columna, los niveles que contienen X forman un grupo de medias, dentro del cual, no hay diferencias estadísticamente significativas. El método actualmente utilizado para discriminar entre las medias es el procedimiento de diferencia mínima (LSD) de Fisher. Con este método, existe un riesgo del 5,0%, de llamar a cada par de medias significativamente diferentes, cuando la diferencia real es igual a 0.
Modelos matemáticos: A continuación, se presenta los modelos matemáticos que describen el comportamiento de la velocidad de corrosión en función del tiempo de exposición en la cámara de niebla salina neutra, para cada gado de oxidación.
Grado A: Velocidad de corrosión = 9222,25-29,4608 *tiempo + 0,0485219 * tiempo2 - 0,0000253615 * tiempo3 R2=100%
Grado B: Velocidad de corrosión= 9509,81- 28,5942 *tiempo + 0,046805 * tiempo2- 0,0,0000246464 * tiempo3 R2=100%
Grado C: Velocidad de corrosión= 5760,08 + 3,84268 *tiempo + 0,0176292 * tiempo2- 0,0000111907 * tiempo3 R2=100%
Grado D: Velocidad de corrosión= 18598,0 - 46,7994*tiempo + 0,0576862 * tiempo2- 0,0000260247 * tiempo3 R2=100%
Demostración de la idoneidad del modelo matemático escogido: El modelo matemático debe cumplir con dos premisas, para que se demuestre su idoneidad:
- Todas las observaciones Yij fueron descritas, según la ecuación del modelo matemático.
- Los errores eij son variables aleatorias independientes.
- eij presenta una distribución normal con media cero y varianza constante para todos los niveles del factor docimado.
- No existe interacción entre el factor analizado y el factor bloque.
Errores aleatorios: El gráfico 1 muestra la distribución de los errores de manera aleatoria, no presentando ningún patrón en su comportamiento.
Distribución normal: El gráfico 2 ilustra el comportamiento de la variable velocidad de corrosión en función de la distribución de probabilidad normal.
Análisis gráfico del comportamiento de los factores y la variable respuesta: En el gráfico 3a, se observa una disminución de la velocidad de corrosión, a medida que transcurre el tiempo de exposición en un mismo grado de oxidación en las muestras.
En el gráfico 3b, se observa una dispersión en los valores de la velocidad de la corrosión en las primeras 250 horas de ensayo, mientras que a las 1.000 presentaron un comportamiento inverso, es decir, se agruparon casi de manera puntual los valores.
Como conclusiones del estudio se destacan:
Se utilizó un diseño de bloques al azar, para la determinación de la influencia del grado de oxidación de las muestras de acero y el factor bloque: tiempo de exposición en la niebla salina neutra, sobre la velocidad de la corrosión, a partir del método gravimétrico.
Se evidencia que los factores (Grado de corrosión y tiempo de exposición en la cámara) poseen los valores de Pvalue inferiores a 0,05, por lo que se rechaza la Hipótesis nula, es decir, existen diferencias significativas entre los tratamientos, lo que es equivalente a plantear que los efectos de los niveles del factor A difieren, con un nivel de significación de 0.05.
Se demostró la idoneidad del modelo matemático, puesto que los residuos cumplen con la condición de la linealidad; no existe una tendencia respecto a los factores en estudio. Se comprobó la condición de homoscedasticidad: varianza constante y la independencia de los errores, puesto que se encuentran normalmente distribuidos y que no presenten ninguna tendencia.
Se observó una disminución de la velocidad de corrosión a medida que transcurrió el tiempo de exposición, para un mismo grado de oxidación en las muestras evaluadas.
Se obtuvo una dispersión en los valores de la velocidad de la corrosión en las primeras horas de ensayo, mientras que a las 1.000 horas presentaron un comportamiento inverso, es decir, se agruparon casi de manera puntual.