Al intentar conseguir un modelo matemático de la extrusora para intentar cerrar el lazo mediante un PID, se ha comprobado que el sistema en lazo abierto no es lo suficientemente estable y el throughput de la extrusora no es constante. Una hipótesis que se tiene es que la alimentación de la granza en el extrusor no es del todo constante. Por ello una solución que se propone es la de intentar que la alimentación sea constante.
Se diseña una tolva que permita meter más granza para comprobar si el peso que ejerce la propia granza en la entrada del alimentador puede ser una posible solución.
La tolva diseñada permite introducir hasta un máximo de $150gr$ frente a los $42gr$ de la tolva de serie. Se va a hacer una producción en la que se va a muestrear los datos obtenidos para su posterior análisis.
In [187]:
import numpy as np
import pandas as pd
import seaborn as sns
%pylab inline
In [188]:
print ("Numpy v{}".format(np.__version__))
print ("Pandas v{}".format(pd.__version__))
print ("Seaborn v{}".format(sns.__version__))
In [219]:
#Abrimos los ficheros con los datos de las pruebas
#Tolva pequeña
data_p1 = pd.read_csv('prueba1.CSV')
#Tolva Graned
data_p2 = pd.read_csv('prueba2.CSV')
In [220]:
#El ensayo se realizo a distintas velocidades, se filtran los datos por
#velocidad rápida= v1 = 4.42mm/s
#velocidad lenta= v2 = 3.816mm/s
#Filtramos los datos con velocidad de tracción rápida
data_p1_v1 = data_p1[(data_p1['RPM TRAC'] >= 4) & (data_p1['RPM TRAC'] <= 4.42)]
data_p2_v1 = data_p2[(data_p2['RPM TRAC'] >= 4) & (data_p2['RPM TRAC'] <= 4.42)]
#Filtramos los datos con velocidad de tracción lenta
data_p1_v2 = data_p1[(data_p1['RPM TRAC'] > 0) & (data_p1['RPM TRAC'] <= 3.816)]
data_p2_v2 = data_p2[(data_p2['RPM TRAC'] > 0) & (data_p2['RPM TRAC'] <= 3.816)]
In [235]:
#d ={'Tolva pequeña V1':data_p1_v1.ix[:,"Diametro X"],'Tolva Grande V1':data_p2_v1.ix[:,"Diametro X"],'Tolva pequeña V2':data_p1_v2.ix[:,"Diametro X"],'Tolva Grande V2':data_p2_v2.ix[:,"Diametro X"]}
d ={'Tolva pequeña':data_p1.ix[:,"Diametro X"],'Tolva Grande':data_p2.ix[:,"Diametro X"]}
da = pd.DataFrame(data=d)
da.boxplot(return_type='axes',figsize=(0,0))
da.describe()
Out[235]:
In [236]:
da.plot(figsize=(20,10))
Out[236]:
Histrogram of the ratio
In [237]:
da.hist(bins=50)
Out[237]: