Compare commits
2 Commits
pv_forecas
...
4727364048
| Author | SHA1 | Date | |
|---|---|---|---|
|
|
4727364048 | ||
|
|
666eb211a3 |
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
Binary file not shown.
87
main.py
87
main.py
@@ -26,57 +26,56 @@ db = DataBaseInflux(
|
|||||||
bucket="allmende_db"
|
bucket="allmende_db"
|
||||||
)
|
)
|
||||||
|
|
||||||
# hp_master = HeatPump(device_name='hp_master', ip_address='10.0.0.10', port=502)
|
hp_master = HeatPump(device_name='hp_master', ip_address='127.0.0.1', port=8111)
|
||||||
# hp_slave = HeatPump(device_name='hp_slave', ip_address='10.0.0.11', port=502)
|
hp_slave = HeatPump(device_name='hp_slave', ip_address='127.0.0.1', port=8111)
|
||||||
# shelly = ShellyPro3m(device_name='wohnung_2_6', ip_address='192.168.1.121')
|
shelly = ShellyPro3m(device_name='wohnung_2_6', ip_address='192.168.1.121')
|
||||||
wr_master = PvInverter(device_name='solaredge_master', ip_address='192.168.1.112', unit=1)
|
wr = PvInverter(device_name='solaredge_master', ip_address='192.168.1.112')
|
||||||
wr_slave = PvInverter(device_name='solaredge_slave', ip_address='192.168.1.112', unit=3)
|
|
||||||
meter = SolaredgeMeter(device_name='solaredge_meter', ip_address='192.168.1.112')
|
meter = SolaredgeMeter(device_name='solaredge_meter', ip_address='192.168.1.112')
|
||||||
|
|
||||||
es.add_components(wr_master, wr_slave)#hp_master, hp_slave, shelly, wr_master, wr_slave, meter)
|
es.add_components(hp_master, hp_slave, shelly, wr, meter)
|
||||||
# controller = SgReadyController(es)
|
controller = SgReadyController(es)
|
||||||
#
|
|
||||||
# # FORECASTING
|
# FORECASTING
|
||||||
# latitude = 48.041
|
latitude = 48.041
|
||||||
# longitude = 7.862
|
longitude = 7.862
|
||||||
# TZ = "Europe/Berlin"
|
TZ = "Europe/Berlin"
|
||||||
# HORIZON_DAYS = 2
|
HORIZON_DAYS = 2
|
||||||
# weather_forecaster = WeatherForecaster(latitude=latitude, longitude=longitude)
|
weather_forecaster = WeatherForecaster(latitude=latitude, longitude=longitude)
|
||||||
# site = Location(latitude=latitude, longitude=longitude, altitude=35, tz=TZ, name="Gundelfingen")
|
site = Location(latitude=latitude, longitude=longitude, altitude=35, tz=TZ, name="Gundelfingen")
|
||||||
#
|
|
||||||
# p_module = 435
|
p_module = 435
|
||||||
# upper_roof_north = PvWattsSubarrayConfig(name="north", pdc0_w=(29+29+21)*p_module, tilt_deg=10, azimuth_deg=20, dc_loss=0.02, ac_loss=0.01)
|
upper_roof_north = PvWattsSubarrayConfig(name="north", pdc0_w=(29+29+21)*p_module, tilt_deg=10, azimuth_deg=20, dc_loss=0.02, ac_loss=0.01)
|
||||||
# upper_roof_south = PvWattsSubarrayConfig(name="south", pdc0_w=(29+21+20)*p_module, tilt_deg=10, azimuth_deg=200, dc_loss=0.02, ac_loss=0.01)
|
upper_roof_south = PvWattsSubarrayConfig(name="south", pdc0_w=(29+21+20)*p_module, tilt_deg=10, azimuth_deg=200, dc_loss=0.02, ac_loss=0.01)
|
||||||
# upper_roof_east = PvWattsSubarrayConfig(name="east", pdc0_w=7*p_module, tilt_deg=10, azimuth_deg=110, dc_loss=0.02, ac_loss=0.01)
|
upper_roof_east = PvWattsSubarrayConfig(name="east", pdc0_w=7*p_module, tilt_deg=10, azimuth_deg=110, dc_loss=0.02, ac_loss=0.01)
|
||||||
# upper_roof_west = PvWattsSubarrayConfig(name="west", pdc0_w=7*p_module, tilt_deg=10, azimuth_deg=290, dc_loss=0.02, ac_loss=0.01)
|
upper_roof_west = PvWattsSubarrayConfig(name="west", pdc0_w=7*p_module, tilt_deg=10, azimuth_deg=290, dc_loss=0.02, ac_loss=0.01)
|
||||||
# cfgs = [upper_roof_north, upper_roof_south, upper_roof_east, upper_roof_west]
|
cfgs = [upper_roof_north, upper_roof_south, upper_roof_east, upper_roof_west]
|
||||||
# pv_plant = PvWattsPlant(site, cfgs)
|
pv_plant = PvWattsPlant(site, cfgs)
|
||||||
#
|
|
||||||
# now = datetime.now()
|
now = datetime.now()
|
||||||
# next_forecast_at = (now + dt.timedelta(hours=1)).replace(minute=0, second=0, microsecond=0)
|
next_forecast_at = (now + dt.timedelta(hours=1)).replace(minute=0, second=0, microsecond=0)
|
||||||
while True:
|
while True:
|
||||||
now = datetime.now()
|
now = datetime.now()
|
||||||
if now.second % interval_seconds == 0 and now.microsecond < 100_000:
|
if now.second % interval_seconds == 0 and now.microsecond < 100_000:
|
||||||
state = es.get_state_and_store_to_database(db)
|
state = es.get_state_and_store_to_database(db)
|
||||||
# mode = controller.perform_action(heat_pump_name='hp_master', meter_name='solaredge_meter', state=state)
|
mode = controller.perform_action(heat_pump_name='hp_master', meter_name='solaredge_meter', state=state)
|
||||||
#
|
|
||||||
# if mode == 'mode1':
|
|
||||||
# mode_as_binary = 0
|
|
||||||
# else:
|
|
||||||
# mode_as_binary = 1
|
|
||||||
# db.store_data('sg_ready', {'mode': mode_as_binary})
|
|
||||||
|
|
||||||
# if now >= next_forecast_at:
|
if mode == 'mode1':
|
||||||
# # Start der Prognose: ab der kommenden vollen Stunde
|
mode_as_binary = 0
|
||||||
# start_hour_local = (now + dt.timedelta(hours=1)).replace(minute=0, second=0, microsecond=0)
|
else:
|
||||||
# weather = weather_forecaster.get_hourly_forecast(start_hour_local, HORIZON_DAYS)
|
mode_as_binary = 1
|
||||||
# total = pv_plant.get_power(weather)
|
db.store_data('sg_ready', {'mode': mode_as_binary})
|
||||||
# db.store_forecasts('pv_forecast', total)
|
|
||||||
#
|
if now >= next_forecast_at:
|
||||||
# # Nächste geplante Ausführung definieren (immer volle Stunde)
|
# Start der Prognose: ab der kommenden vollen Stunde
|
||||||
# # Falls wir durch Delay mehrere Stunden verpasst haben, hole auf:
|
start_hour_local = (now + dt.timedelta(hours=1)).replace(minute=0, second=0, microsecond=0)
|
||||||
# while next_forecast_at <= now:
|
weather = weather_forecaster.get_hourly_forecast(start_hour_local, HORIZON_DAYS)
|
||||||
# next_forecast_at = (next_forecast_at + dt.timedelta(hours=1)).replace(minute=0, second=0, microsecond=0)
|
total = pv_plant.get_power(weather)
|
||||||
|
db.store_forecasts('pv_forecast', total)
|
||||||
|
|
||||||
|
# Nächste geplante Ausführung definieren (immer volle Stunde)
|
||||||
|
# Falls wir durch Delay mehrere Stunden verpasst haben, hole auf:
|
||||||
|
while next_forecast_at <= now:
|
||||||
|
next_forecast_at = (next_forecast_at + dt.timedelta(hours=1)).replace(minute=0, second=0, microsecond=0)
|
||||||
|
|
||||||
|
|
||||||
time.sleep(0.1)
|
time.sleep(0.1)
|
||||||
|
|||||||
238
pv_inverter.py
238
pv_inverter.py
@@ -1,155 +1,139 @@
|
|||||||
# pv_inverter.py
|
|
||||||
# -*- coding: utf-8 -*-
|
|
||||||
from typing import Optional, Dict, Any, List
|
|
||||||
from pymodbus.client import ModbusTcpClient
|
|
||||||
from pymodbus.exceptions import ModbusIOException
|
|
||||||
import struct
|
|
||||||
import time
|
import time
|
||||||
|
import struct
|
||||||
|
import pandas as pd
|
||||||
|
from typing import Dict, Any, List, Tuple, Optional
|
||||||
|
from pymodbus.client import ModbusTcpClient
|
||||||
|
|
||||||
|
EXCEL_PATH = "modbus_registers/pv_inverter_registers.xlsx"
|
||||||
|
|
||||||
|
# Obergrenze: bis EXKLUSIVE 40206 (d.h. max. 40205)
|
||||||
|
MAX_ADDR_EXCLUSIVE = 40121
|
||||||
|
|
||||||
class PvInverter:
|
class PvInverter:
|
||||||
"""
|
def __init__(self, device_name: str, ip_address: str, port: int = 502, unit: int = 1):
|
||||||
Minimaler Reader für einen SolarEdge-Inverter hinter Modbus-TCP→RTU-Gateway.
|
|
||||||
Liest nur die bekannten Register (wie im funktionierenden Skript).
|
|
||||||
Kompatibel mit pymodbus 2.5.x und 3.x – kein retry_on_empty.
|
|
||||||
"""
|
|
||||||
|
|
||||||
def __init__(
|
|
||||||
self,
|
|
||||||
device_name: str,
|
|
||||||
ip_address: str,
|
|
||||||
port: int = 502,
|
|
||||||
unit_id: int = 1,
|
|
||||||
timeout: float = 1.5,
|
|
||||||
silent_interval: float = 0.02,
|
|
||||||
):
|
|
||||||
self.device_name = device_name
|
self.device_name = device_name
|
||||||
self.host = ip_address
|
self.ip = ip_address
|
||||||
self.port = port
|
self.port = port
|
||||||
self.unit = unit_id
|
self.unit = unit
|
||||||
self.timeout = timeout
|
|
||||||
self.silent_interval = silent_interval
|
|
||||||
self.client: Optional[ModbusTcpClient] = None
|
self.client: Optional[ModbusTcpClient] = None
|
||||||
self._connect()
|
self.registers: Dict[int, Dict[str, Any]] = {} # addr -> {"desc":..., "type":...}
|
||||||
|
self.connect_to_modbus()
|
||||||
|
self.load_registers(EXCEL_PATH)
|
||||||
|
|
||||||
# ---------------- Verbindung ----------------
|
# ---------- Verbindung ----------
|
||||||
def _connect(self):
|
def connect_to_modbus(self):
|
||||||
# retries=0: keine internen Mehrfachversuche
|
self.client = ModbusTcpClient(self.ip, port=self.port, timeout=3.0, retries=3)
|
||||||
self.client = ModbusTcpClient(self.host, port=self.port, timeout=self.timeout, retries=0)
|
|
||||||
if not self.client.connect():
|
if not self.client.connect():
|
||||||
raise ConnectionError(f"Verbindung zu {self.device_name} ({self.host}:{self.port}) fehlgeschlagen.")
|
print("❌ Verbindung zu Wechselrichter fehlgeschlagen.")
|
||||||
print(f"✅ Verbindung hergestellt zu {self.device_name} ({self.host}:{self.port}, unit={self.unit})")
|
raise SystemExit(1)
|
||||||
|
print("✅ Verbindung zu Wechselrichter hergestellt.")
|
||||||
|
|
||||||
def close(self):
|
def close(self):
|
||||||
if self.client:
|
if self.client:
|
||||||
self.client.close()
|
self.client.close()
|
||||||
self.client = None
|
self.client = None
|
||||||
|
|
||||||
# ---------------- Low-Level Lesen ----------------
|
# ---------- Register-Liste ----------
|
||||||
def _read_regs(self, addr: int, count: int) -> Optional[List[int]]:
|
def load_registers(self, excel_path: str):
|
||||||
"""Liest 'count' Holding-Register ab base-0 'addr' für die konfigurierte Unit-ID."""
|
xls = pd.ExcelFile(excel_path)
|
||||||
try:
|
df = xls.parse()
|
||||||
rr = self.client.read_holding_registers(address=addr, count=count, slave=self.unit)
|
# Passe Spaltennamen hier an, falls nötig:
|
||||||
except ModbusIOException:
|
cols = ["MB Adresse", "Beschreibung", "Variabel Typ"]
|
||||||
time.sleep(self.silent_interval)
|
df = df[cols].dropna()
|
||||||
return None
|
df["MB Adresse"] = df["MB Adresse"].astype(int)
|
||||||
except Exception:
|
|
||||||
time.sleep(self.silent_interval)
|
|
||||||
return None
|
|
||||||
|
|
||||||
time.sleep(self.silent_interval)
|
# 1) Vorab-Filter: nur Adressen < 40206 übernehmen
|
||||||
if not rr or rr.isError():
|
df = df[df["MB Adresse"] < MAX_ADDR_EXCLUSIVE]
|
||||||
return None
|
|
||||||
return rr.registers
|
|
||||||
|
|
||||||
|
self.registers = {
|
||||||
|
int(row["MB Adresse"]): {
|
||||||
|
"desc": str(row["Beschreibung"]).strip(),
|
||||||
|
"type": str(row["Variabel Typ"]).strip()
|
||||||
|
}
|
||||||
|
for _, row in df.iterrows()
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
# ---------- Low-Level Lesen ----------
|
||||||
|
def _try_read(self, fn_name: str, address: int, count: int) -> Optional[List[int]]:
|
||||||
|
fn = getattr(self.client, fn_name)
|
||||||
|
# pymodbus 3.8.x hat 'slave='; Fallbacks schaden nicht
|
||||||
|
for kwargs in (dict(address=address, count=count, slave=self.unit),
|
||||||
|
dict(address=address, count=count)):
|
||||||
|
try:
|
||||||
|
res = fn(**kwargs)
|
||||||
|
if res is None or (hasattr(res, "isError") and res.isError()):
|
||||||
|
continue
|
||||||
|
return res.registers
|
||||||
|
except TypeError:
|
||||||
|
continue
|
||||||
|
return None
|
||||||
|
|
||||||
|
def _read_any(self, address: int, count: int) -> Optional[List[int]]:
|
||||||
|
regs = self._try_read("read_holding_registers", address, count)
|
||||||
|
if regs is None:
|
||||||
|
regs = self._try_read("read_input_registers", address, count)
|
||||||
|
return regs
|
||||||
|
|
||||||
|
# ---------- Decoding ----------
|
||||||
@staticmethod
|
@staticmethod
|
||||||
def _to_int16(u16: int) -> int:
|
def _to_i16(u16: int) -> int:
|
||||||
return struct.unpack(">h", struct.pack(">H", u16))[0]
|
return struct.unpack(">h", struct.pack(">H", u16))[0]
|
||||||
|
|
||||||
@staticmethod
|
@staticmethod
|
||||||
def _apply_sf(raw: int, sf: int) -> float:
|
def _to_f32_from_two(u16_hi: int, u16_lo: int, msw_first: bool = True) -> float:
|
||||||
return raw * (10 ** sf)
|
b = struct.pack(">HH", u16_hi, u16_lo) if msw_first else struct.pack(">HH", u16_lo, u16_hi)
|
||||||
|
return struct.unpack(">f", b)[0]
|
||||||
|
|
||||||
|
# Hilfsfunktion: wie viele 16-Bit-Register braucht dieser Typ?
|
||||||
@staticmethod
|
@staticmethod
|
||||||
def _read_string_from_regs(regs: List[int]) -> Optional[str]:
|
def _word_count_for_type(rtype: str) -> int:
|
||||||
b = b"".join(struct.pack(">H", r) for r in regs)
|
rt = (rtype or "").lower()
|
||||||
s = b.decode("ascii", errors="ignore").rstrip("\x00 ").strip()
|
# Passe hier an deine Excel-Typen an:
|
||||||
return s or None
|
if "uint32" in rt or "real" in rt or "float" in rt or "string(32)" in rt:
|
||||||
|
return 2
|
||||||
|
# Default: 1 Wort (z.B. int16/uint16)
|
||||||
|
return 1
|
||||||
|
|
||||||
# ---------------- Hilfsfunktionen ----------------
|
def read_one(self, address_excel: int, rtype: str) -> Optional[float]:
|
||||||
def _read_string(self, addr: int, words: int) -> Optional[str]:
|
"""
|
||||||
regs = self._read_regs(addr, words)
|
Liest einen Wert nach Typ ('INT' oder 'REAL' etc.).
|
||||||
if regs is None:
|
Es werden ausschließlich Register < 40206 gelesen.
|
||||||
|
"""
|
||||||
|
addr = int(address_excel)
|
||||||
|
words = self._word_count_for_type(rtype)
|
||||||
|
|
||||||
|
# 2) Harte Grenze prüfen: höchstes angefasstes Register muss < 40206 sein
|
||||||
|
if addr + words - 1 >= MAX_ADDR_EXCLUSIVE:
|
||||||
|
# Überspringen, da der Lesevorgang die Grenze >= 40206 berühren würde
|
||||||
return None
|
return None
|
||||||
return self._read_string_from_regs(regs)
|
|
||||||
|
|
||||||
def _read_scaled(self, value_addr: int, sf_addr: int) -> Optional[float]:
|
if words == 2:
|
||||||
regs = self._read_regs(value_addr, 1)
|
regs = self._read_any(addr, 2)
|
||||||
sf = self._read_regs(sf_addr, 1)
|
if not regs or len(regs) < 2:
|
||||||
if regs is None or sf is None:
|
return None
|
||||||
return None
|
# Deine bisherige Logik interpretiert 2 Worte als Float32:
|
||||||
raw = self._to_int16(regs[0])
|
return self._to_f32_from_two(regs[0], regs[1])
|
||||||
sff = self._to_int16(sf[0])
|
|
||||||
return self._apply_sf(raw, sff)
|
|
||||||
|
|
||||||
def _read_u32_with_sf(self, value_addr: int, sf_addr: int) -> Optional[float]:
|
|
||||||
regs = self._read_regs(value_addr, 2)
|
|
||||||
sf = self._read_regs(sf_addr, 1)
|
|
||||||
if regs is None or sf is None:
|
|
||||||
return None
|
|
||||||
u32 = (regs[0] << 16) | regs[1]
|
|
||||||
sff = self._to_int16(sf[0])
|
|
||||||
return self._apply_sf(u32, sff)
|
|
||||||
|
|
||||||
# ---------------- Öffentliche API ----------------
|
|
||||||
def get_state(self) -> Dict[str, Any]:
|
|
||||||
"""Liest exakt die bekannten Register und gibt ein Dict zurück."""
|
|
||||||
state: Dict[str, Any] = {}
|
|
||||||
|
|
||||||
# --- Common Block ---
|
|
||||||
state["C_Manufacturer"] = self._read_string(40004, 16)
|
|
||||||
state["C_Model"] = self._read_string(40020, 16)
|
|
||||||
state["C_Version"] = self._read_string(40044, 8)
|
|
||||||
state["C_SerialNumber"] = self._read_string(40052, 16)
|
|
||||||
|
|
||||||
# --- Inverter Block ---
|
|
||||||
state["I_AC_Power_W"] = self._read_scaled(40083, 40084)
|
|
||||||
state["I_AC_Voltage_V"] = self._read_scaled(40079, 40082)
|
|
||||||
state["I_AC_Frequency_Hz"] = self._read_scaled(40085, 40086)
|
|
||||||
state["I_DC_Power_W"] = self._read_scaled(40100, 40101)
|
|
||||||
state["I_AC_Energy_Wh_total"] = self._read_u32_with_sf(40093, 40095)
|
|
||||||
|
|
||||||
status_regs = self._read_regs(40107, 2)
|
|
||||||
if status_regs:
|
|
||||||
state["I_Status"] = status_regs[0]
|
|
||||||
state["I_Status_Vendor"] = status_regs[1]
|
|
||||||
else:
|
else:
|
||||||
state["I_Status"] = None
|
regs = self._read_any(addr, 1)
|
||||||
state["I_Status_Vendor"] = None
|
if not regs:
|
||||||
|
return None
|
||||||
|
return float(self._to_i16(regs[0]))
|
||||||
|
|
||||||
return state
|
def get_state(self) -> Dict[str, Any]:
|
||||||
|
"""
|
||||||
|
Liest ALLE Register aus self.registers und gibt dict zurück.
|
||||||
# ---------------- Beispiel ----------------
|
Achtet darauf, dass keine Adresse (inkl. Mehrwort) >= 40206 gelesen wird.
|
||||||
if __name__ == "__main__":
|
"""
|
||||||
MODBUS_IP = "192.168.1.112"
|
data = {"Zeit": time.strftime("%Y-%m-%d %H:%M:%S")}
|
||||||
MODBUS_PORT = 502
|
for address, meta in sorted(self.registers.items()):
|
||||||
|
words = self._word_count_for_type(meta["type"])
|
||||||
master = PvInverter("solaredge_master", MODBUS_IP, port=MODBUS_PORT, unit_id=1)
|
# 3) Nochmals Schutz auf Ebene der Iteration:
|
||||||
slave = PvInverter("solaredge_slave", MODBUS_IP, port=MODBUS_PORT, unit_id=3)
|
if address + words - 1 >= MAX_ADDR_EXCLUSIVE:
|
||||||
|
continue
|
||||||
try:
|
val = self.read_one(address, meta["type"])
|
||||||
sm = master.get_state()
|
if val is None:
|
||||||
ss = slave.get_state()
|
continue
|
||||||
|
key = f"{address} - {meta['desc']}"
|
||||||
print("\n=== MASTER ===")
|
data[key] = val
|
||||||
for k, v in sm.items():
|
return data
|
||||||
print(f"{k:22s}: {v}")
|
|
||||||
|
|
||||||
print("\n=== SLAVE ===")
|
|
||||||
for k, v in ss.items():
|
|
||||||
print(f"{k:22s}: {v}")
|
|
||||||
|
|
||||||
finally:
|
|
||||||
master.close()
|
|
||||||
slave.close()
|
|
||||||
|
|||||||
Binary file not shown.
Reference in New Issue
Block a user