Diferencias

Muestra las diferencias entre dos versiones de la página.

--- acemu:ensayos:jornadas:2013_04_13 [2013/04/13 17:16]
kenneth [Prueba de Motor LAB001 - EMA - 13/04/2013]
+++ acemu:ensayos:jornadas:2013_04_13 [2013/04/25 09:40] (actual)
tabare.perez [Prueba de Motor LAB001 - EMA - 13/04/2013]
@@ Línea 1: / Línea 1: @@
+[[acemu:ensayos:jornadas|Retorno a página anterior]]
+----
 ====== Prueba de Motor LAB001 - EMA - 13/04/2013 ======
@@ Línea 5: / Línea 9: @@
 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 kg/volt
+Se hicieron dos ensayos de disparo de motor. Previamente se calibró la balanza para obtener un valor de 20.1 kg/volt (196.98 N/volt)
 **Motor:** LAB001
@@ Línea 15: / Línea 19: @@
 Datos de empuje: {{:acemu:ensayos:jornadas:2013_04_13:2013_04_13_curva_lab001_prueba_01.csv|2013_04_13:2013_04_13_curva_lab001_prueba_01.csv}}
-Video del ensayo: {{:acemu:ensayos:jornadas:2013_04_13_Motor_LAB001_ensayo1.ogg|Ensayo 1}}
+Video del ensayo: [[http://acemu.org/wiki/lib/exe/fetch.php/acemu:ensayos:jornadas:2013_04_13:ensayo1.ogg|Ensayo 1]]
 **Ensayo 2:**
@@ Línea 23: / Línea 28: @@
 Datos de empuje: {{:acemu:ensayos:jornadas:2013_04_13:2013_04_13_curva_lab001_prueba_02.csv|2013_04_13:2013_04_13_curva_lab001_prueba_02.csv}}
-Video del ensayo: {{:acemu:ensayos:jornadas:2013_04_13_Motor_LAB001_ensayo2.ogg|Ensayo 2}}
+Video del ensayo: [[http://acemu.org/wiki/lib/exe/fetch.php/acemu:ensayos:jornadas:2013_04_13:ensayo2.ogg|Ensayo 2]]
 ===== Resultados =====
@@ Línea 31: / Línea 36: @@
 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.
 ==== Apéndice ====
+=== Archivos CSV del DataQ ===
+El formato de los archivos de la toma de datos de empuje con la balanza es el siguiente (se le agregó una header con los nombres de los campos de cada registro, para que R tome esos nombres, o sea que lo que muestro es ligeramente diferente que el header del archivo del DataQ):
+<code>
+time,ch1,ch2,ch3,ch4,date,hour
+.933,-5.0781E-01, 0.0000E+00, 2.4414E+00, 2.4414E+00,04-13-13,14:28:13
+.950,-5.0781E-01, 0.0000E+00, 2.4414E+00, 2.4414E+00,04-13-13,14:28:13
+.967,-5.0781E-01, 0.0000E+00, 2.4414E+00, 2.4414E+00,04-13-13,14:28:13
+.983,-5.0781E-01, 0.0000E+00, 2.4414E+00, 2.4414E+00,04-13-13,14:28:13
+.000,-5.0781E-01, 0.0000E+00, 2.4414E+00, 2.4414E+00,04-13-13,14:28:13
+. . .
+</code>
+Los datos motrados corresponden a las primeras lineas del archivo CSV del ensayo 1
 === Script en R para procesar los datos de los archivos CSV ===
 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.
@@ Línea 39: / Línea 59: @@
 curva01<-read.csv("2013_04_13_curva_lab001_prueba_01.csv")
 head(curva01)
-plot(curva01$time,curva01$ch1)
+impulso01<-sum((curva01$ch1+0.50781)*20.1*9.8)/60
-# si tenemos en cuenta que la curva de empuje corresponde a 20.2 kg/volt
+curva02<-read.csv("2013_04_13_curva_lab001_prueba_02.csv")
+head(curva02)
+impulso02<-sum((curva02$ch1+0.50781)*20.1*9.8)/60
+# si tenemos en cuenta que la curva de empuje corresponde a 20.1 kg/volt (196.98N)
 # 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")
+plot(curva02$time-1.783,(curva02$ch1+0.50781)*20.1*9.8,main="Curva Empuje LAB01 - 13/04/2013 - Ensayos 1 y 2",xlab="tiempo(s)",ylab="Empuje (N)",col=4,lwd=3,type="l")
-# Hacemos lo mismo para la segunda curva de empuje
+lines(curva01$time-2.89,(curva01$ch1+0.50781)*20.1*9.8,lwd=3,col=2)
-curva02<-read.csv("2013_04_13_curva_lab001_prueba_02.csv")
+text(0.8,480,"Ensayo 1: 202g KNDX",col=2,adj=c(0,0))
-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")
+text(0.8,450,paste("Impulso: ",as.character(format(impulso01,digits=3)),"Ns"),col=2,adj=c(0,0))
-# Ahora ponemos los dos ensayos juntos
+text(0.8,400,"Ensayo 2: 214g KNDX",col=4,adj=c(0,0))
-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")
+text(0.8,370,paste("Impulso: ",as.character(format(impulso02,digits=3)),"Ns"),col=4,adj=c(0,0))
-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)
@@ Línea 58: / Línea 78: @@
+=== Calibración de la Balanza ===
+Script de R para procesar datos de calibración de balanza
+<code bash>
+# Calibración de la balanza
+# Ingreso de datos tomados: voltaje y peso
+volt=c(-0.547,-0.449,-0.352,-0.254,-0.156,0.039,0.293)
+weight=c(0,1980,3895,5910,7893,9871,11843,16843)
+# convertimos gramos en kilos
+weight<-weight/1000
+# Ajuste por mínimos cuadrados
+calibration <- lm(weight~volt)
+# Resumen del ajuste
+summary(calibration)
+# Ploteo de puntos y linea de regresión
+plot(volt,weight,lwd=3,main="Calibración de balanza - 13/04/2013",xlab="Voltaje celda de carga",ylab="Peso en la balanza (kg)",pch=5)
+lines(volt,calibration$fitted,col="red",lwd=3)
+# copiamos a device x11 y salvamos a disco como bitmap
+dev.copy(device=x11)
+savePlot("2013_04_13_calibracion_balanza.png",type="png")
+</code>
+Resumen del ajuste por mínimos cuadrados:
+<code>
+Call:
+lm(formula = weight ~ volt)
+Residuals:
+      Min        1Q    Median        3Q       Max
+-0.060084 -0.019677  0.000872  0.024861  0.047092
+Coefficients:
+            Estimate Std. Error t value Pr(>|t|)
+(Intercept) 11.01174    0.01778   619.4 1.20e-15 ***
+volt        20.10699    0.05606   358.6 3.17e-14 ***
+---
+Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
+Residual standard error: 0.04077 on 6 degrees of freedom
+Multiple R-squared:     1,	Adjusted R-squared: 0.9999
+F-statistic: 1.286e+05 on 1 and 6 DF,  p-value: 3.172e-14
+</code>
+Tal como se desprende del ajuste, la pendiente de la recta peso vs. voltaje nos da 20.1 kg/volt y el valor de los diferentes tests estadísticos dan una muy buena correlación lineal.
+El gráfico obtenido al procesar los datos se muestra en la figura siguiente:
+{{ :acemu:ensayos:jornadas:2013_04_13:2013_04_13_calibracion_balanza.png?400 | Gráfico de Calibración de Balanza}}
+----
+[[acemu:ensayos:jornadas|Retorno a página anterior]]

A.C.E.M.U.

Herramientas de usuario

Herramientas del sitio

Diferencias

Herramientas de la página