Diferencias

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

--- acemu:ensayos:jornadas:2013_04_13 [2013/04/14 05:46]
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.1 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 55: / Línea 59: @@
 curva01<-read.csv("2013_04_13_curva_lab001_prueba_01.csv")
 head(curva01)
+impulso01<-sum((curva01$ch1+0.50781)*20.1*9.8)/60
 curva02<-read.csv("2013_04_13_curva_lab001_prueba_02.csv")
-# si tenemos en cuenta que la curva de empuje corresponde a 20.1 kg/volt
+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(curva02$time-1.783,20.1*curva02$ch1+0.50781*20.1,main="Curva Empuje LAB01 - 13/04/2013 - Ensayos 1 y 2",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")
-lines(curva01$time-2.89,20.1*curva01$ch1+0.50781*20.1,lwd=3,col=2)
+lines(curva01$time-2.89,(curva01$ch1+0.50781)*20.1*9.8,lwd=3,col=2)
-text(1.1,45,"Ensayo 1:  202g KNDX",col=2)
+text(0.8,480,"Ensayo 1: 202g KNDX",col=2,adj=c(0,0))
-text(1.1,40,"Ensayo 2:  214g KNDX",col=4)
+text(0.8,450,paste("Impulso: ",as.character(format(impulso01,digits=3)),"Ns"),col=2,adj=c(0,0))
+text(0.8,400,"Ensayo 2: 214g KNDX",col=4,adj=c(0,0))
+text(0.8,370,paste("Impulso: ",as.character(format(impulso02,digits=3)),"Ns"),col=4,adj=c(0,0))
 # copiamos a device x11 y salvamos a disco la copia generada como archivo PNG
 dev.copy(device=x11)
@@ Línea 69: / 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