Zulauf einer Zisterne in Liter/h berechnen

Aiouh

@paul53
Die Anzeige des Tankvolumens klappt nun gut. Nachdem es heute endlich mal wieder etwas Regen gab, habe ich festgestellt, dass das mit der Durchflussmessung noch nicht klappt.
Die kleinste Volumen Auflösung die ich systembedingt hinbekomme (+/-10mV) sind 29 Liter, was für mich völlig OK ist.
Die Durchflussmenge pendelt jedoch zwischen 0 und 1000, offenbar so wie das letzte digit springt.
Hab ich da noch was übersehen?

Spannung vom Geber (408 =4,08V)
408 	true 	modbus.2 	2020-06-04 19:42:06.895 	
409 	true 	modbus.2 	2020-06-04 19:42:17.889 	
408 	true 	modbus.2 	2020-06-04 19:42:56.992 	
409 	true 	modbus.2 	2020-06-04 19:43:08.092 	
408 	true 	modbus.2 	2020-06-04 19:43:18.889 	
409 	true 	modbus.2 	2020-06-04 19:43:29.889 	
408 	true 	modbus.2 	2020-06-04 19:43:58.990 	
409 	true 	modbus.2 	2020-06-04 19:44:19.888 	
408 	true 	modbus.2 	2020-06-04 19:44:49.889 	
409 	true 	modbus.2 	2020-06-04 19:45:00.955 	
408 	true 	modbus.2 	2020-06-04 19:48:46.886 	
409 	true 	modbus.2 	2020-06-04 19:48:57.886

Volumen in Liter
8055 	true 	javascript.0 	2020-06-04 19:42:06.905 	
8084 	true 	javascript.0 	2020-06-04 19:42:17.897 	
8055 	true 	javascript.0 	2020-06-04 19:42:57.005 	
8084 	true 	javascript.0 	2020-06-04 19:43:08.101 	
8055 	true 	javascript.0 	2020-06-04 19:43:18.896 	
8084 	true 	javascript.0 	2020-06-04 19:43:29.896 	
8055 	true 	javascript.0 	2020-06-04 19:43:58.999 	
8084 	true 	javascript.0 	2020-06-04 19:44:19.894 	
8055 	true 	javascript.0 	2020-06-04 19:44:49.896 	
8084 	true 	javascript.0 	2020-06-04 19:45:01.026 	
8055 	true 	javascript.0 	2020-06-04 19:48:46.894 	
8084 	true 	javascript.0 	2020-06-04 19:48:57.894

Durchfluss 
0 	true 	javascript.0 	2020-06-04 19:42:06.905 	
1000 	true 	javascript.0 	2020-06-04 19:42:17.897 	
0 	true 	javascript.0 	2020-06-04 19:42:57.005 	
1000 	true 	javascript.0 	2020-06-04 19:43:08.101 	
0 	true 	javascript.0 	2020-06-04 19:43:18.896 	
1000 	true 	javascript.0 	2020-06-04 19:43:29.896 	
0 	true 	javascript.0 	2020-06-04 19:43:58.999 	
1000 	true 	javascript.0 	2020-06-04 19:44:19.894 	
0 	true 	javascript.0 	2020-06-04 19:44:49.896 	
1000 	true 	javascript.0 	2020-06-04 19:45:01.027 	
0 	true 	javascript.0 	2020-06-04 19:48:46.894 	
1000 	true 	javascript.0 	2020-06-04 19:48:57.895

paul53

@Aiouh
Der Wert schwankt offenbar, wobei er auch kleiner wird. Ergänze mal um eine Abfrage

    let V = h2V(h);
    if(V > lastV) {
        let fluss = 3600000 * (V - lastV) / (dp.state.lc - lc); // l/h
        lastV = V;
        lc = dp.state.lc;
        setState(idFluss, Math.round(fluss * 100) / 100, true); // 2 Nachkommastellen
    }

und vor dem Trigger:

var lc = getState(idSpannung).lc;

Die 10 Minuten, nach denen der Zufluss auf 0 gesetzt wird, könnten zu kurz sein.

Aiouh

@paul53
Nachdem es wieder angefangen hat zu regnen tut sich was, das sieht gut aus.

Was mir noch aufgefallen ist: Der letzte Durchflusswert nach Ende Zufluss scheint jetzt so lange stehen zu bleiben bis erneut Wasser zufliesst. Der Wert geht nicht auf Null.
Bespiel:
Den letzten Zufluss hatte ich am 5.6. um 14:50 Uhr, der Stand von 109.91 blieb bis heute Nacht 00:06 Uhr stehen, als es zum ersten mal wieder geregnet hat.

109.91 	true 	influxdb.0 	2020-06-05 14:50:13.837 	
0.77 	true 	javascript.0 	2020-06-07 00:06:39.974

Hab ich vielleicht was mit dem Code verbaselt? Sieht jetzt so aus.

const idSpannung = 'modbus.2.holdingRegisters.9_Zisterne_Füllhöhe_U';
const idHoehe = 'javascript.0.Eigene_Datenpunkte.Zisterne_Füllhöhe';
const idVol = 'javascript.0.Eigene_Datenpunkte.Zisterne_Volumen';
const idFluss = 'javascript.0.Eigene_Datenpunkte.Zisterne_Zufluss';

const L = 49.8; // Länge in dm
const R = 8; // Radius in dm
 
function u2h(U) {
   let h = (U-120) /24; // a und b durch konstante Werte ersetzen
   setState(idHoehe, Math.round(h * 100) / 10, true); // in cm mit 1 Nachkommastelle
   return h; // in dm
}

function h2V(h) {
   let V = L * (R * R * Math.acos(1 - h / R) - (R - h) * Math.sqrt(h * (2 * R - h)));
   setState(idVol, Math.round(V), true);
   return V; // in l
}

var h = u2h(getState(idSpannung).val);
var lastV = h2V(h);
var timer = null;
var lc = getState(idSpannung).lc; 

on(idSpannung, function(dp) {
    h = u2h(dp.state.val);
    let V = h2V(h);
    if(V > lastV) {
      let fluss = 3600000 * (V - lastV) / (dp.state.lc - lc); // l/h
    lastV = V;
    lc = dp.state.lc;
    setState(idFluss, Math.round(fluss * 100) / 100, true); // 2 Nachkommastellen
    }
});

paul53

@Aiouh sagte:

Sieht jetzt so aus.

Es fehlt der letzte Teil mit dem Timeout:

on(idSpannung, function(dp) {
    h = u2h(dp.state.val);
    let V = h2V(h);
    if(V > lastV) {
       let fluss = 3600000 * (V - lastV) / (dp.state.lc - lc); // l/h
       lastV = V;
       lc = dp.state.lc;
       setState(idFluss, Math.round(fluss), true); 
       if(timer) clearTimeout(timer);
       timer = setTimeout(function() {
          setState(idFluss, 0, true);
          timer = null;
       }, 1800000); // 30 Minuten
    }
});

Aiouh

@paul53
Hallo, jetzt scheint alles zu klappen.
Vielen Dank für Deine Hilfe UND Geduld :blush:

Aiouh

@paul53
Hallo, ich hab da noch eine Bitte :grin: ...
Mein liegender Tank hat zwar nur einen Durchmesser von 160cm, durch das Domschacht Volumen und das Entnahmerohr erhöht sich jedoch die gemessene maximale Höhe des Wassers, wodurch das angezeigte Tankvolumen daneben liegt. (Der Überlauf liegt höher als der Domschacht)
Wie bekomme ich das hin, dass jede Höhe über "2 x R" einfach ignoriert wird und darüber einfach das maximale Volumen angezeigt wird ? Falls es keine Mühe macht wäre ich dankbar wenn das im gleichen Skript eingebaut wäre, ansonsten bastle ich mal mit Blockly an dem Wert Füllhöhe_U herum.

paul53

@Aiouh sagte in:

jede Höhe über "2 x R" einfach ignoriert wird

function u2h(U) {
   let h = (U-120) /24;
   if(h > 2 * R) h = 2 * R;
   setState(idHoehe, Math.round(h * 100) / 10, true); // in cm mit 1 Nachkommastelle
   return h; // in dm
}

Aiouh

@paul53 :+1:

Aiouh

@paul53
Das Ganze läuft mittlerweile ganz gut, mit einer Ausnahme...
Ich bekomme immer wieder mal falsche Werte, irgendwo zwischen 0 bis ca 500 unterhalb des aktuellen Wertes. Bisher lösche ich die händisch aus der InfluxDB.
Würde meine zusätzliche Zeile " if ((V < (lastV - 500)) V = lastV; " funktionieren?
Ich möchte damit einfach Werte ignorieren die plötzlich mindestens 500 l niedriger liegen.

function u2h(U) {
   let h = (U-120) /24;
   if(h > 2 * R) h = 2 * R;
   if ((V < (lastV - 500)) V = lastV;
   setState(idHoehe, Math.round(h * 100) / 10, true); // in cm mit 1 Nachkommastelle
   return h; // in dm
}

paul53

@Aiouh sagte:

Würde meine zusätzliche Zeile " if ((V < (lastV - 500)) V = lastV; " funktionieren?

Nein, das ist die falsche Stelle. Wird auch die Höhe aufgezeichnet ?

Aiouh

@paul53 Ja, die Höhe wird auch aufgezeichnet

paul53

@Aiouh sagte:

Ja, die Höhe wird auch aufgezeichnet

Da die Höhe zuerst ermittel wird, müsste man dort schon eine entsprechende Beschränkung einbauen. Das neu berechnete Volumen liegt zu dem Zeitpunkt noch nicht vor.
Einfacher wäre eine Auswertung der Spannungsänderung.

on(idSpannung, function(dp) {
    if(dp.state.val > dp.oldState.val - 20) {

Das setzt voraus, dass die starke Abweichung jeweils nur bei einem Messwert auftritt und danach der Messwert wieder stimmt.

Aiouh

@paul53 Danke, das werde ich mal versuchen. Bisher war immer nur ein Messwert falsch, dann kamen wieder viele korrekte.

Aiouh

@paul53 Hallo, ich habe mich nochmal hingesetzt und in der Datenbank nachgeschaut.
Die Grafik hatte zwar Zwischenwerte angezeigt, lag wohl am Zoom factor, in der Datenbank waren die falschen Werte immer "0". Nachdem ich händisch alle Nullwerte entfernt habe ist die Grafik wieder perfekt.
Folgende Überlegung:
Wenn die Spannung unter 120 liegt, das ist quasi die Null-Höhe, soll einfach der letzte Spannungswert verwendet werden.
Würde das mit dem Syntax "if(V <120) V = lastV;" klappen?

on(idSpannung, function(dp) {
    h = u2h(dp.state.val);
    let V = h2V(h);
    if(V <120) V = lastV;
    if(V > lastV) {
       let fluss = 3600000 * (V - lastV) / (dp.state.lc - lc); // l/h
       lastV = V;
       lc = dp.state.lc;

paul53

@Aiouh sagte:

Würde das mit dem Syntax "if(V <120) V = lastV;" klappen?

Teste es.

Aiouh

@paul53 hmmm, bis jetzt bekomme ich zumindest mal normale Werte, habe soweit also nichts "kaputt gemacht". Bin gespannt wie es die nächsten Tage aussieht

Aiouh

@paul53 Funktioniert mit einer Ausnahme. Jede Nacht um 0:00Uhr bekomme ich "0".

paul53

@Aiouh sagte:

Jede Nacht um 0:00Uhr bekomme ich "0".

Da es im Script keinen Bezug zur Uhrzeit gibt, musst Du wohl die Ursache an der Quelle suchen.

Aiouh

@paul53 korrekt, werde auch nach der Ursache suchen. Bin aber schon mal happy dass "meine Zeile" zu funktionieren scheint. :grinning:

Burg_Io

@paul53

Bei meiner Suche nach der Berechnung des Inhaltes meiner Zisterne (liegender Zylinder) bin ich auf diesen Beitrag gestoßen.
Meine Vorstellung ist es mit den Datenpunkten Radius, Länge und Füllhöhe (in dm), das Volumen meiner Zisterne in l zu berechnen und in einen Datenpunkt zu schreiben:

const idR = 0_userdata.0.Test.Zisterne.Zisterne_Radius_dm       // Radius in dm
const idL = 0_userdata.0.Test.Zisterne.Zisterne_Laenge_dm       // Länge in dm
const idVol = 0_userdata.0.Test.Zisterne.Zisterne_Volumen         // zu schreibender Wert in l
const idh = 0_userdata.0.Test.Zisterne.Zisterne_Fuellhoehe_dm   // Höhes in dm

function h2V(h) {
   let V = L * (R * R * Math.acos(1 - h / R) - (R - h) * Math.sqrt(h * (2 * R - h)));
   setState(idVol, Math.round(V), true);
   return V; // in l
}

Danke im Voraus für jeglichen Hinweis

Community Forum

NEWS

Zulauf einer Zisterne in Liter/h berechnen

Support us

743

32.5k

81.7k

1.3m