Engee 文档
Notebook

超级电容器的充放电

在此示例中,对超级电容器进行建模,在其充电和随后的放电期间从中获取电压读数。 为了给超级电容器充电100秒,从市电向其提供100mA的电流。 然后超级电容器一分钟不使用。 在接下来的一个小时内,每50秒连接一个50mA负载一秒钟,以使超级电容器放电。

模型图:

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]: