¡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:
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")