Engee 文档
Notebook

超级电容器的充放电

在本示例中,超级电容器在充电和放电时会产生电压读数。从电源施加 100 mA 的电流,对超级电容器充电 100 秒。然后,超级电容器在一分钟内不使用。在接下来的一小时内,每隔 50 秒施加 50 毫安的负载,持续一秒钟,对超级电容器放电。

模型示意图:

supercapacitor_charging_and_discharging_behavior_1743426433042.png

定义加载和运行模型的函数:

In [ ]:
function start_model_engee()
    try
        engee.close("supercapacitor_charging_and_discharging_behavior", force=true) # закрытие модели 
        catch err # в случае, если нет модели, которую нужно закрыть и engee.close() не выполняется, то будет выполнена её загрузка после catch
            m = engee.load("$(@__DIR__)/supercapacitor_charging_and_discharging_behavior.engee") # загрузка модели
        end;

    try
        engee.run(m, verbose=true) # запуск модели
        catch err # в случае, если модель не загружена и engee.run() не выполняется, то будут выполнены две нижние строки после catch
            m = engee.load("$(@__DIR__)/supercapacitor_charging_and_discharging_behavior.engee") # загрузка модели
            engee.run(m, verbose=true) # запуск модели
        end
end
Out[0]:
start_model_engee (generic function with 1 method)

运行模拟

In [ ]:
start_model_engee();
Building...
Progress 0%
Progress 5%
Progress 10%
Progress 15%
Progress 21%
Progress 26%
Progress 31%
Progress 36%
Progress 41%
Progress 47%
Progress 52%
Progress 57%
Progress 62%
Progress 68%
Progress 73%
Progress 78%
Progress 83%
Progress 88%
Progress 93%
Progress 98%
Progress 100%

将模拟数据写入变量:

In [ ]:
t = simout["supercapacitor_charging_and_discharging_behavior/SysOutput_1"].time[:]
load_current = simout["supercapacitor_charging_and_discharging_behavior/SysOutput_1"].value[:]
supercapacitor_voltage = simout["supercapacitor_charging_and_discharging_behavior/SysOutput_2"].value[:]
Out[0]:
WorkspaceArray{Float64}("supercapacitor_charging_and_discharging_behavior/SysOutput_2").value

数据可视化

In [ ]:
using Plots
In [ ]:
p1 = plot(t, load_current, linewidth=2, label="Нагрузка, мА", ylabel="Ток, мА", lc=2)
p2 = plot(t, supercapacitor_voltage, linewidth=2, label="Напряжение суперконденсатора", xlabel="Время, с", ylabel="Напражение, В")
plot(p1, p2, layout=(2,1))
Out[0]: