Tema: Control exhaustivo ciclo de carga/SOC/Cargas/etc

Iniciado por Gabriel 2015

El apantallado podría solucionar el problema. Si tienes por ahí un cacho de algo, prueba y vemos.

Sigue saltando, aún habiendo conectado la malla a Tierra.

- - - - - - - - - - - - - - - - - - -

Manuel, una pregunta sobre el código:

Código PHP:

                if diaok == 1 and horaok == 1:
                        Rele_Out_H[GR][Puerto]+=1      <------?
                        #print "registro activo"
                        #print 


El + (de +=1) ¿es un error al pulsar el teclado? ¿o es que sumas +1 al resultado a cada paso del FOR si se cumplen las 2 condiciones?

Sumo +1

La idea es que puedes tener varias condiciones de horario para un unico rele...

Por ejm uno todos los dias de 10:00 a 12:00 y otro solo los Domingos de 10:00 a 13:00 y de 17:00 a 19:00

En principio por cada rele seria un OR de todas sus condiciones....y por tanto sumo 1 si se cumple una condicion.....y despues al acabar el bucle de condiciones veo por cada rele si tiene un valor > cero......en mi ejm el Domingo a las 11:00 el algoritmo dara un 2

OK, tengo el código listo.

A ver si puedes probar estos cambios para el tema 'hora' ????

Código PHP:

#al empezar
import time

# en bucle principal

        diasemana = time.strftime("%w")        #dia semana, OK
        hora = time.strftime("%H:%M:%S")     #el que uso yo


##### esto en el bucle de condiciones horarias ######

                if (hora >= str(H[I][OPH['valor_h_ON']]) and
                    hora <= str(H[I][OPH['valor_h_OFF']])):
                        horaok=1 


Lo tengo en marcha y no me ha dado problemas, pero como corre en segundo plano no veo nada. Eso sí, sin errores.

EDITO: fijate que en mi caso tengo valor_h_ON y OFF porque así lo definí. Tanto en el diccionario como en la tabla.

Tal como lo tienes te da error lógico por eso de que "2" es mayor que "11"

Habria que poner la hora con dos digitos y pensar que hacer para el caso de ON=23:00:00 y OFF= 04:00:00

Código PHP:

>>> print hora
23:40:25
>>> print str(H[I][OPH['valor_h_on']])
7:00:00
>>> print str(H[I][OPH['valor_h_off']])
23:50:00
print hora >= str(H[I][OPH['valor_h_on']])
False 


Quizás el error está en valor_h_on, debería ser 07:00:00. ¿no?

Iniciado por nikitto

Quizás el error está en valor_h_on, debería ser 07:00:00. ¿no?

Claro...supongo que hay una forma facil de que te ponga en el texto el "0" en las horas < 10:00

Si no lo encuentras,..... se puede hacer por programa viendo la longitud de la cadena, etc.......o usar los segundos como hice yo


...... hecho

Código PHP:

>>> print str(H[I][OPH['valor_h_on']]).zfill(8)
07:00:00 


Ambos codigos dan el mismo resultado con unas pruebas que he hecho

Código PHP:

if (hora >= str(H[I][OPH['valor_h_on']]).zfill(8) and 
                    hora <= str(H[I][OPH['valor_h_off']]).zfill(8)): 
                        horaok=1 


Código PHP:

if (segundosdia >= H[I][OPH['valor_h_on']].seconds and 
                    segundosdia <= H[I][OPH['valor_h_off']].seconds):
                        horaok2=1 


Queda los del cambio de dia dentro del periodo.. me lo pienso

Buena idea, rellenar el valor con ceros al principio sino tiene una longitud de 8. Estas hecho un crack.

Ya tienes el código modificado en el directorio.

Pues, o estoy hoy un poco espeso, o en primera aproximación, no es evidente codificar que se pueda poner periodos que pasen de las 23:59 por ejm:

- ON a las 04:00 y OFF a las 02:00 del día siguiente
- ON a las 22:00 y OFF a las 00:30 del día siguiente (típico de luces nocturnas)
- .......

Mientras lo sacamos toca dar de alta dos registros......por ejm para el caso de las luces

- ON a las 22:00 y OFF a las 23:59:59
- ON a las 00:00 y OFF a las 00:30

...... edito ......

Pues nada mejor que llevar a los hijos temprano al aeropuerto para despejar mente...la solución muy simple...hago que el algoritmo desdoble la condición en los dos registros comentados en el anterior ejemplo

El algoritmo de control horario seria......

Primero sacamos la hora actual teniendo en cuenta el tema de "0" a la izquierda cuando la hora es menor de las 10:00

Código PHP:

hora = time.strftime("%H:%M:%S").zfill(8) #No necesario .zfill(8) dado que la funcion ya pone los ceros a la izquierda 


la parte del control horario

Código PHP:

# -------------------- Bucle de condiciones de horario --------------------------
 
        Rele_Out_H = [[0] * 8 for i in range(8)] # Inicializo variable a 0
        
        for I in range(hcon):  # Bucle por cada registro de condiciones de horario                 
                print "Procesando ",H[I]
                Puerto=(H[I][0]%10)-1   
                GR=int((H[I][0]-Puerto)/10-32) 

                diaok=0
                horaok=0
                                
                if H[I][OPH['parametro_h']]=='T': #Todos los dias de la semana
                        diaok=1
                elif str(NDIA[H[I][OPH['parametro_h']]])==diasemana:
                        diaok=1

                if str(H[I][OPH['valor_h_on']]).zfill(8)>str(H[I][OPH['valor_h_off']]).zfill(8): #True si periodo pasa por 0:00
                        if (hora >= str(H[I][OPH['valor_h_on']]).zfill(8) and 
                            hora <= "23:59:59"): 
                                horaok=1                                                       
                        if (hora >= "00:00:00" and 
                            hora <= str(H[I][OPH['valor_h_off']]).zfill(8)): 
                                horaok=1  

                elif (hora >= str(H[I][OPH['valor_h_on']]).zfill(8) and  # Periodo normal en el mismo dia
                    hora <= str(H[I][OPH['valor_h_off']]).zfill(8)): 
                        horaok=1  
                                       
                if diaok == 1 and horaok == 1:
                        Rele_Out_H[GR][Puerto]+=1
                        print "registro activo"
                        print
                
        print 'Rele_Out_H=',Rele_Out_H  

        for I in range(hcon):
                Puerto=(H[I][0]%10)-1   
                GR=int((H[I][0]-Puerto)/10-32) 
                
                if Rele_Out_H[GR][Puerto]==0:
                        Rele_Out[GR][Puerto]=0 #apago rele 


Con esto debería funcionar ya el control por hora y día de la semana, permitiendo que el periodo de encendido cruce la medianoche



Faltaria como poco implementar en web el estado de los reles y la posibilidad de actualizar/crear registros de reles, y condiciones

Mientras se implementa via web, se puede funcionar sin problemas remotamente, conectándose a la BD y modificando/creando los registros

Luego,..... lo que es la funcionalidad mínima que se especificó inicialmente.......ya estaría acabada

En el camino, al menos a mi, me han salido algún temilla adicional a considerar

A) Conocer la potencia máxima disponible en placas en cada momento
B) Controlar reles via wifi/inhalambrico
C) Controlar local o remotamente dispositivos con control de potencia por angulo de fase (asociado a la A) para minimizar ciclado baterias)

D) Sacar datos (Vbat, Ibat, SOC,...) de dispositivos que se tengan......por ejm en mi caso que tengo:
- Un FM80.....Vbat, Vplacas, Iplacas, Status carga (bulk, abs, flot, equ),.......
- Monitor de baterias......Vbat, Ibat, SOC,...

Podría hacer una instalación mas simple.....


No obstante, este BRICO nació con la idea de ser genérico, y por tanto, es bueno tener uno mismo la posibilidad de sacar por su cuenta los parámetros que quiera del sistema FV (Vbat, Ibat, SOC, temperatura, etc)

Una cosa más sobre el tema del horario que a mi no me funcionó ayer.
Por lo que vi si sólo existe la condición horaria no funciona. Verás en Menú-Relés que puse LUCES_EXT, pues no hubo forma de que se conectara dicho relé.

Código actualizado en dir.

- - - - - - - - - - - - - - - - - - -

Otro problemilla más, al llegar las 11h, no se puso en marcha la depuradora (las condiciones Fv eran favorables). Por ahora he vuelvo al programa anterior.

He realizado una bateria de pruebas y, por ahora, creo que el algoritmo va bien

Los datos dados de alta en la BD son:

Resumen Tablas Reles:



En mi caso como no lo tengo en un sistema real, simulo ,los valores de Vbat, SOC y hora para ver si la salida de los reles es acorde a la programación y los valores simulados de acuerdo a esta tabla



Corriendo el programa ...me sale coherente

Para una mejor logica he puesto que las condiciones se cumplen siempre que sean >= o <= en lugar de > o <.... es decir que si la condicion es OFF si Vbat<24, el programa apaga el Rele si Vbat=24 en lugar de apagarlo si Vbat=23,999

Código PHP:

-----------------------------NUEVO REGISTRO--------------------------------
hora= 12:00:00
Vbat =  25
SOC= 100 %
Situacion Reles antes de algoritmo= [0, 0, 0, 0, 0, 0, 0, 0]

Situacion Reles final= [0, 0, 0, 0, 0, 0, 0, 0]

-----------------------------NUEVO REGISTRO--------------------------------
hora= 13:30:00
Vbat =  25
SOC= 100 %
Situacion Reles antes de algoritmo= [0, 0, 0, 0, 0, 0, 0, 0]

Situacion Reles final= [0, 0, 1, 0, 0, 0, 0, 0]

-----------------------------NUEVO REGISTRO--------------------------------
hora= 15:00:00
Vbat =  26
SOC= 100 %
Situacion Reles antes de algoritmo= [0, 0, 1, 0, 0, 0, 0, 0]

Situacion Reles final= [1, 0, 1, 0, 0, 0, 0, 0]

-----------------------------NUEVO REGISTRO--------------------------------
hora= 15:00:00
Vbat =  26
SOC= 100 %
Situacion Reles antes de algoritmo= [1, 0, 1, 0, 0, 0, 0, 0]

Situacion Reles final= [1, 0, 1, 1, 0, 0, 0, 0]

-----------------------------NUEVO REGISTRO--------------------------------
hora= 15:00:00
Vbat =  24
SOC= 100 %
Situacion Reles antes de algoritmo= [1, 0, 1, 1, 0, 0, 0, 0]

Situacion Reles final= [0, 0, 1, 1, 0, 0, 0, 0]

-----------------------------NUEVO REGISTRO--------------------------------
hora= 14:45:00
Vbat =  28
SOC= 100 %
Situacion Reles antes de algoritmo= [0, 0, 1, 1, 0, 0, 0, 0]

Situacion Reles final= [1, 0, 1, 1, 0, 0, 0, 0]

-----------------------------NUEVO REGISTRO--------------------------------
hora= 15:00:00
Vbat =  28
SOC= 100 %
Situacion Reles antes de algoritmo= [1, 0, 1, 1, 0, 0, 0, 0]

Situacion Reles final= [1, 1, 1, 1, 0, 0, 0, 0]

-----------------------------NUEVO REGISTRO--------------------------------
hora= 19:00:00
Vbat =  28
SOC= 100 %
Situacion Reles antes de algoritmo= [1, 1, 1, 1, 0, 0, 0, 0]

Situacion Reles final= [0, 1, 0, 1, 0, 0, 0, 0]

-----------------------------NUEVO REGISTRO--------------------------------
hora= 19:00:00
Vbat =  28
SOC= 85 %
Situacion Reles antes de algoritmo= [0, 1, 0, 1, 0, 0, 0, 0]

Situacion Reles final= [0, 0, 0, 1, 0, 0, 0, 0]

-----------------------------NUEVO REGISTRO--------------------------------
hora= 19:00:00
Vbat =  28
SOC= 83 %
Situacion Reles antes de algoritmo= [0, 0, 0, 1, 0, 0, 0, 0]

Situacion Reles final= [0, 0, 0, 0, 0, 0, 0, 0]
>>> 


Me surge una duda con este hilo. A la hora de calcular el SOC, se tienen en cuenta como de rápidas son las descargas. Me explico, si la baterías descargan 100 A durante una hora se obtiene el mismo resultado que cuando descargan 10 horas 10 Amperios?

Iniciado por Gabriel 2015

Me surge una duda con este hilo. A la hora de calcular el SOC, se tienen en cuenta como de rápidas son las descargas. Me explico, si la baterías descargan 100 A durante una hora se obtiene el mismo resultado que cuando descargan 10 horas 10 Amperios?

Claro que lo tenemos en cuenta ....usamos Peukert para adaptar el SOC según el régimen de descarga en cada muestra

El calculo/algoritmo y consideraciones esta explicado en este brico

Casi al principio:

http://www.solarweb.net/forosolar/fotovoltaica-sistemas-aislados-la-red/39930-control-exhaustivo-ciclo-carga-soc-cargas-etc-2.html#post313293

Versión Python:

Código PHP:

        ### CALCULO SOC% A C20
        if ibat < 0 :
                Ip1 = -ibat
                Ip1 = Ip1**CP
                Ip1 = AH*Ip1

                Ip2 = AH / 20
                Ip2 = (Ip2**CP)*20
                Ip= -Ip1/Ip2
        else :
                Ip = ibat * EC

        ### esta parte es para corregir poco a poco la desviación real de carga ###
        ### cuando se cumplen las siguientes condiciones es que la batería está prácticamente llena###
        ### algunos reguladores pasan p.e. de 95% a 100% en un segundo, nosotros le damos más tiempo 

        if (ibat>0 and ibat<1 and abs(vbat-vflotacion)<0.2) :
                DS = DS + (AH-DS)/50
        else :
                DS = DS + (Ip * Tiempo_horas)
        ####    

        if DS > AH :
                DS = AH
        if DS < 0 :
                DS = 0

        #print('AH='),AH
        #print('DS='),DS
        #print('DS1='),DS1
        #print('DS2='),DS2
        #print('Ip='),Ip

        if DS > DS1 :
                DS1 = DS
        if DS < DS2 :
                DS2 = 0

        soc = DS/AH*100
        #print('SOC='),soc,('%')
        #print(DS1/AH*100),('%')
        #print(DS2/AH*100),('%')

        ### FIN CALCULO SOC% 


He pillado el fallo que mencionaba Nikitto en la lectura del PCF y copia del valor en la variable Rele_Out_Ant

En realidad son dos fallos


El primero es, de nuevo , por no considerar los 0 a la izquierda en cadenas de texto

Como habia hecho pruebas simulando con los valores de 255 y 254,... obviamente siempre tenia un byte completo de 8 bits

Cuando la lectura da un valor de 4 por ejm, la cadena binaria que se obtiene es "100" que hay que pasar a "00000100" aplicando lo mismo que hicimos con la hora......".zfill(8)"



El segundo es que el valor a introducir en Rele_Out_Ant es con logica positiva mientras que el valor leido del PCF es con logica negativa....luego hay que intercambiar 0 por 1 (yo lo hago con un xor de 255)


Por tanto el código actualizado para copiar el valor de la situación de los PCF en la variable Rele_Out_Ant seria:

Código PHP:

for I in range(nreles):
           Puerto=(TR[I][0]%10)-1
           GR=int((TR[I][0]-Puerto)/10-32)
           addr=GR+32
           estado=bus.read_byte(addr)  #devuelve el valor en decimal del PCF en direccion addr
           estado=bin(estado ^ 255)[2:10].zfill(8)
           Rele_Out_Ant[GR][Puerto]=int(estado[7-Puerto]) 


Voy a probar el nuevo código.

Dale un vistazo a la parte de Menú - Relés. Añadir relé y añadir horario para relé, ya funcionan. Me falta implementar la parte de las condiciones FV.
A ver si más adelante me ayuda algún conocido para dejar la cosa más bonita y agradable a la vista.

- - - - - - - - - - - - - - - - - - -

He visto algo raro, si el relé 328 (por poner el caso que me he encontrado) está definido en la tabla 'reles' pero no en la tabla 'reles_c' ni 'reles_h', al ejecutar el programa con el núevo código, el relé queda activado.

Siento no controlar el interfaz web...pero seguro que encontramos formas de mejorarla visual y operativamente

Veo que has puesto algo para añadir...pero supongo que se usara mas modificar que añadir.... es igualmente valido?

Lo que si hay que poner es algo para cambiar el modo de PRG a ON o a OFF para poder encender o apagar manualmente

- - - - - - - - - - - - - - - - - - -

Iniciado por nikitto

- - - - - - - - - - - - - - - - - - -

He visto algo raro, si el relé 328 (por poner el caso que me he encontrado) está definido en la tabla 'reles' pero no en la tabla 'reles_c' ni 'reles_h', al ejecutar el programa con el núevo código, el relé queda activado.

Si es facil que pase por la codificación hecha, ...... me miro como evitarlo, pero chequea que hace si ese Rele 328 sin condiciones si le poner el modo en OFF

Es un fallo lógico, porque cogemos los reles de la tabla 'reles' y allí sí que existe. Aunque no tiene lógica que esté definido aquí y no en en rele_c y/o rele_h. Sería un caso 'sin sentido'.
Me he dado cuenta porque estaba jugando con la creación de reles y demás, y he visto que el primer número de la izda. (relé 8) estaba encendido (=0). Al borrarlo de la base de datos, se ha solucionado el problema.

- - - - - - - - - - - - - - - - - - -

Por ahora sólo añadir. Quiero ver y aprender como hacer algo más usable. También tengo la página RELES preparada (falta activar) para usar contraseña (MD5). Si no cualquiera me la lía parda desde casa. Jajaja.

Cambia la parte de forzado ON/OFF de sitio en el codigo poniendo TR en lugar de R segun esto...

Código PHP:


# -----------------------Bucle encendido/apagado reles -----------------

        Flag_Rele_Encendido = 0

        for I in range(nreles):
                
                Puerto=(TR[I][0]%10)-1   # desde 0 a 7
                GR=int((TR[I][0]-Puerto)/10-32) # desde 0 donde GR=0 si id_rele=32X, GR=1 si id_rele=33X

                ### forzado ON/OFF
                if TR[I][OP['modo']]=='ON' :
                        Rele_Out[GR][Puerto]=1
                if TR[I][OP['modo']]=='OFF' :
                        Rele_Out[GR][Puerto]=0

                ............. 


De este tipo de errores seguro saldran unos pocos

Iniciado por Mleon

Claro que lo tenemos en cuenta ....usamos Peukert para adaptar el SOC según el régimen de descarga en cada muestra

El calculo/algoritmo y consideraciones esta explicado en este brico

Sí, como el hilo tiene tantos post, no lo recordaba. El caso es que tengo que hacer algo parecido porque en una instalación no me queda otra que utilizar una batería. Tenía pensado utilizar los datos que da el fabricante de C5, C10, C20 y C100 e interpolar entre los mismos con Peukert, o algo parecido.

Iniciado por Gabriel 2015

Sí, como el hilo tiene tantos post, no lo recordaba. El caso es que tengo que hacer algo parecido porque en una instalación no me queda otra que utilizar una batería. Tenía pensado utilizar los datos que da el fabricante de C5, C10, C20 y C100 e interpolar entre los mismos con Peukert, o algo parecido.

Comprueba con los datos CX del fabricante que te da siempre el mismo coeficiente de Peukert (deberia ser asi) y aplica el algoritmo en cada muestra para calcular la Ip (intensidad de Peukert equivalente) con la I medida en descarga y con el dato calculado vas actualizando el SOC

Iniciado por Mleon

Comprueba con los datos CX del fabricante que te da siempre el mismo coeficiente de Peukert (deberia ser asi) y aplica el algoritmo en cada muestra para calcular la Ip (intensidad de Peukert equivalente) con la I medida en descarga y con el dato calculado vas actualizando el SOC

Voy a hacerlo. A ver qué sale. Tratándose de baterías, seguro que el índice no sale igual...

Pues la primera, en la frente. Viendo los datos del fabricante, me sale que el famoso índice cambia de 1,21 a 1,07...

Intensidad Peukert Horas
79 1,21 2
37 1,21 5


37 1,2 5
20,7 1,2 10


20,7 1,14 10
11,25 1,14 20


11,25 1,07 20
2,5 1,07 100

Aunque estamos en verano el desarrollo del brico para RPi sigue. Ahora andamos liados con el manejo de las tablas de la base de datos.

Recordáis que había un MODO: PRG (por programa=condiciones FV y horario), ON (a piñón) y OFF (desconectado 4ever). Pues ya está implementado. Desde la web con el PC, móvil, tableta, podemos activar/desactivar y programar los relés.

Algunas imágenes:







Quiero aprovechar una vez más para dar las gracias a Manuel (mleon), sin ti esto no sería posible; a el_cobarde, por sus lecciones de FV, de las que sigo tomando notas cada día en una libretita que tengo por aquí; a Gabriel 2015, del que seguro también he aprendido algunas cosillas y seguiré aprendiendo; y a todo el foro en general, por no dejarme a nadie. GRACIAS.

Seguimos trabajando ...

PD: Cambiaré donde se escribe el MODO, por algún desplegable con las tres opciones. Para que no se pueda poner nada distinto.

Iniciado por Gabriel 2015

Pues la primera, en la frente. Viendo los datos del fabricante, me sale que el famoso índice cambia de 1,21 a 1,07...

Intensidad Peukert Horas
79 1,21 2
37 1,21 5


37 1,2 5
20,7 1,2 10


20,7 1,14 10
11,25 1,14 20


11,25 1,07 20
2,5 1,07 100

Estas considerando desde C2 a C100...luego puede ser normal que no sea lineal en todo en intervalo

No se la precision que necesitas y en que I te vas a mover normalmente...pero si pones un termino medio (1,14) calcula cuanto se desvia la autonomia con respecto a los datos fabricante

Si quieres mayor precision..pon que el Coef. Peuket vaya cambiando segun la I medida

- - - - - - - - - - - - - - - - - - -

Iniciado por nikitto

......
Quiero aprovechar una vez más para dar las gracias a Manuel (mleon), sin ti esto no sería posible;

.....

Gracias...lo que es seguro es que si ti esto no seria posible.......estaba yo tan tranquilo con mi implementacion en duinomite (en basic) .... y me has puesto a aprender MySQL, Python, web,....

Creo que esta quedando curiosa la implementación en RPi .....cuando le añadamos el control de excedentes por angulo de fase y conmutación remota via wifi nos queda para nota

Por lo pronto, si te parece bien, vamos a añadir como parámetro de condiciones FV que se pueda poner la condición de flotación de las baterías ( "FLOT" por ejm)y asi poder encender o apagar según este parámetro

Después añadiremos algun que otro parámetro (temperaturas, etc) que pueda ser interesante considerar

Tal y como lo hemos codificado...añadir parámetros es muy muy facil


Lo he comentado varias veces...pero reto a que me enseñen algo parecido en potencia/precio para controlar una instalación FV

Unos pocos cambios más:



-----------

Manuel, estuve mirando lo de usar polar chart para los relés y no me convence.

Luego miré lo de solidgauge y va a ser que voy a probarlo. Lo bueno es que el grafico sólo se define una vez, luego cada uno de los parámetros por separado. En el ejemplo de HighCharts: Speed y RPM.