¡Esta es una revisión vieja del documento!
Presentes: Enrique Madruga, Tabaré Pérez y Kenneth Irving Ausentes: el resto del mundo…
Se hicieron dos ensayos de disparo de motor. Previamente se calibró la balanza para obtener un valor de 20.3 kb/volt
Motor: LAB001
Ensayo 1:
Combustible: 202 g de KNDX (batch 1 del 10/11/2012)
Datos de empuje: 2013_04_13:2013_04_13_curva_lab001_prueba_01.csv
Video del ensayo:
Ensayo 2:
Combustible: 214 g de KNDX (batch 3/11/2012)
Datos de empuje: 2013_04_13:2013_04_13_curva_lab001_prueba_02.csv
Video del ensayo:
Tal como surge de las curvas de empuje, los resultados son muy similares, la diferencia se debe a que en el ensayo 2 la masa de KNDX era un poco mayor. El perfil de las curvas es idéntico y corresponden a dos batches preparados en fechas diferentes, por lo tanto dos cocinadas de candy independientes. Eso demuestra que hemos logrado una buena uniformidad en la preparación de candy.
De los archivos originales levantados con el Dataq, se recortaron y aislaron los segmentos donde están los datos de la curva de empuje. Estos datos en formato CSV se leyeron con el programa R (RStudio) y se procesaron con los siguientes comandos par generar la figura de más arriba.
El script empleado es el siguiente:
# Leemos los datos de la curva de empuje curva01<-read.csv("2013_04_13_curva_lab001_prueba_01.csv") head(curva01) plot(curva01$time,curva01$ch1) # si tenemos en cuenta que la curva de empuje corresponde a 20.2 kg/volt # y que arrancamos con la señal de -0.50781, podemos ya graficar directamente # en kg vs tiempo plot(curva01$time-2.933,20.3*curva01$ch1+0.50781*20.3,main="Curva Empuje LAB01 - 13/04/2013 - Ensayo 1",xlab="tiempo(s)",ylab="Empuje (Kg)",col=4,lwd=3,type="l") # Hacemos lo mismo para la segunda curva de empuje curva02<-read.csv("2013_04_13_curva_lab001_prueba_02.csv") plot(curva02$time-1.783,20.3*curva02$ch1+0.50781*20.3,main="Curva Empuje LAB01 - 13/04/2013 - Ensayo 2",xlab="tiempo(s)",ylab="Empuje (Kg)",col=4,lwd=3,type="l") # Ahora ponemos los dos ensayos juntos plot(curva02$time-1.783,20.3*curva02$ch1+0.50781*20.3,main="Curva Empuje LAB01 - 13/04/2013 - Ensayos 1 y 2",xlab="tiempo(s)",ylab="Empuje (Kg)",col=4,lwd=3,type="l") lines(curva01$time-2.89,20.3*curva01$ch1+0.50781*20.3,lwd=3,col=2) text(1.1,45,"Ensayo 1: 202g KNDX",col=2) text(1.1,40,"Ensayo 2: 214g KNDX",col=4) # copiamos a device x11 y salvamos a disco la copia generada como archivo PNG dev.copy(device=x11) savePlot("2013_04_13_graficos_ensayos.png",type="png")