Моделирование катушки индуктивности с гистерезисом.
В данном примере показано, как изменение коэффициентов уравнения магнитного гистерезиса Джайлса-Атертона влияет на результирующую кривую B-H. Параметры моделирования настроены на выполнение четырех полных циклов переменного тока с начальной напряженностью поля (H) и плотностью магнитного потока (B), установленными на ноль.
Схема модели:
Определение функции для загрузки и запуска модели:
function start_model_engee()
try
engee.close("inductor_with_hysteresis", force=true) # закрытие модели
catch err # в случае, если нет модели, которую нужно закрыть и engee.close() не выполняется, то будет выполнена её загрузка после catch
m = engee.load("$(@__DIR__)/inductor_with_hysteresis.engee") # загрузка модели
end;
try
engee.run(m) # запуск модели
catch err # в случае, если модель не загружена и engee.run() не выполняется, то будут выполнены две нижние строки после catch
m = engee.load("$(@__DIR__)/inductor_with_hysteresis.engee") # загрузка модели
engee.run(m) # запуск модели
end
end
Запуск модели:
try
start_model_engee() # запуск симуляции с помощью специальной функции, реализованной выше
catch err
end;
Вывод результатов моделирования из переменной simout:
res = collect(simout)
Запись результатов в переменные:
H = collect(res[4]) # напряжённость поля
B = collect(res[5]) # плотность магнитного потока
using Plots
plot(H[100:end,2], B[100:end,2], linewidth=3, xlabel= "H, А/м", ylabel= "B, Т", legend=:bottomright)
Определение нового коэффициента обратной намагниченности:
Вывод параметров блока Nonlinear Inductor:
engee.get_param("inductor_with_hysteresis/Nonlinear Inductor")
Переопределяем значения коэффициентов уравнения магнитного гистерезиса с помощью функции set_param!, скрытой под маской:
engee.load("$(@__DIR__)/inductor_with_hysteresis.engee")
# @markdown Коэффициент обратимой намагниченности:
c = 0.2 # @param {type:"slider", min:0, max:1, step:0.1}
# @markdown Коэффициент объемной связи:
K = 200.0 # @param {type:"slider", min:0, max:1000, step:1}
# @markdown Коэффициент междоменной связи:
alpha = 0.0001 # @param {type:"slider", min:0, max:0.001, step:0.0001}
engee.set_param!("inductor_with_hysteresis/Nonlinear Inductor", "c" => c)
engee.set_param!("inductor_with_hysteresis/Nonlinear Inductor", "K" => K)
engee.set_param!("inductor_with_hysteresis/Nonlinear Inductor", "alpha" => alpha)
engee.run("inductor_with_hysteresis")
res1 = collect(simout)
H1 = collect(res1[7])
B1 = collect(res1[9])
engee.close("inductor_with_hysteresis", force = true)
plot(H[100:end,2], B[100:end,2], label = "Исходные параметры модели", xlabel= "H, А/м", ylabel= "B, Т", linewidth=3)
plot!(H1[100:end,2], B1[100:end,2], label = "Изменённые параметры модели", legend=:bottomright, linewidth=3)
Выводы:
В данном примере с помощью программного управления была проведена симуляция модели катушки индуктивности с гистерезисом. На графиках показано, как отдельные коэффициенты уравнений Джайлса-Атертона влияют на кривую гистерезиса.