| G | L \
|___| C /
- |___________| |_________|
4 | | 3
"""
import bms
from bms.physical.electrical import Generator, Resistor, ElectricalNode, Capacitor, Ground, Inductor
from bms.signals.functions import Sinus
U = Sinus('Generator', 2, 5) # Voltage of generator
R = 10 # Resistance in Ohm
C = 0.01 # Capacitance in Fahrads
L = 0.025 # Inductance in Henry
n1 = ElectricalNode('1')
n2 = ElectricalNode('2')
n3 = ElectricalNode('3')
n4 = ElectricalNode('4')
gen = Generator(n4, n1, U)
res = Resistor(n1, n2, R)
cap = Capacitor(n2, n3, C)
ind = Inductor(n3, n4, L)
gnd = Ground(n4)
ps = bms.PhysicalSystem(4, 300, [gen, res, cap, gnd, ind], [])
ds = ps.dynamic_system
#ds._ResolutionOrder3()
d = ds.Simulate()
2
"""
import bms
from bms.physical.electrical import Battery, Resistor, ElectricalNode, Ground
from bms.signals.functions import Sinus
# U=Sinus('Generator',2,5)# Voltage of generator
R = 10 # Resistance in Ohm
Rint = 10
Umax = 12.5
Umin = 6
C = 3600 * 10
soc = 0.8
n1 = ElectricalNode('1')
n2 = ElectricalNode('2')
# n3=ElectricalNode('3')
Bat = Battery(n1, n2, Umin, Umax, C, soc, Rint)
Res = Resistor(n2, n1, R)
G = Ground(n1)
ps = bms.PhysicalSystem(5000, 50, [Bat, Res, G], [])
ds = ps.dynamic_system
# ds._ResolutionOrder3()
d = ds.Simulate()
ds.PlotVariables([[n1.variable, n2.variable], [Bat.soc], [Bat.variables[0]]])
# Validation: analytical solutions