Расчёт тиристора
alexevsУправляемый тиристор Т123-200 с охладителем О123
В этом интерактивном скрипте показана методика разработки и проверки модели однофазного тиристорного регулятора с силовым тиристором Т123-200 и охладителем О123. Материал полезен инженерам и студентам, которые строят физические модели силовой электроники в Engee и хотят связать параметры блоков с паспортными данными прибора, тепловым расчетом и алгоритмом импульсно-фазового управления. На примере модели VT.engee рассмотрены выбор уровня детализации тиристора, расчет потерь открытого состояния, подбор тепловой цепи, обоснование RL-нагрузки, RC-цепи защиты, развязанной цепи управления gate и подсистемы СИФУ. Кодовые ячейки выполняют контрольные аналитические вычисления, запускают модель без изменения исходного файла и строят графики тока, напряжения, усредненных величин и регулировочной характеристики. Формулы и порядок расчета опираются на datasheet и справочную литературу.
1. Методика построения модели
Разработка модели силового тиристорного регулятора начинается не с выбора блоков, а с выбора физической границы задачи. В учебной и инженерной постановке удобно разделить модель на четыре части: силовую цепь, цепь управления, систему синхронизации и фазового управления, тепловую цепь. Такое разделение повторяет структуру реального устройства: энергия передается через анод-катод тиристора и нагрузку, управляющий импульс вводится через отдельный маломощный канал, а тепловая мощность отводится от полупроводниковой структуры через корпус и радиатор.
В Engee для силового тиристора выбран блок Thyristor, а не идеальный ключ. Причина принципиальна: в расчете нужны ток удержания, ток и напряжение управления, падение напряжения в открытом состоянии, тепловой порт и табличная ВАХ. Идеальный ключ хорош для логической проверки алгоритма, но он не дает корректных потерь и температуры. Документация Engee прямо указывает, что блок Thyristor можно параметризовать либо эквивалентной транзисторной схемой, либо таблицей ВАХ открытого состояния; для паспортной модели промышленного прибора удобнее таблица, потому что она ближе к графикам datasheet и устойчивее в длительном моделировании [8].
Силовая часть строится как управляемое выпрямление на индуктивной нагрузке. Индуктивность не допускает мгновенного обрыва тока, поэтому в модели нужна свободноходовая ветвь с диодом. RC-цепь ставится параллельно тиристору как демпфер коммутационных перенапряжений и как элемент ограничения скорости изменения напряжения . Тепловой порт тиристора соединяется с радиатором и тепловой массой среды.
2. От паспорта тиристора к параметрам блока
Паспортные данные сначала переводятся в набор параметров, которые реально используются моделью. Для Т123-200 в расчете нужны: класс повторяющегося напряжения , средний и действующий токи открытого состояния и , ток удержания , параметры управления и , пороговое напряжение открытого состояния , динамическое сопротивление , предельная температура перехода , тепловые сопротивления и .
Выбор класса напряжения выполняют от максимального напряжения на закрытом тиристоре с запасом. Для сети амплитуда равна
Класс 6 на для такой сети недостаточен, поэтому в модели установлен , соответствующий исполнению Т123-200-16. Это дает запас на сетевой пик и повторяющиеся перенапряжения.
Открытое состояние задается табличной зависимостью . Если паспорт дает график ВАХ, из него выбирают несколько точек и заносят в I_on_vector и V_on_vector. Если доступна только линейная аппроксимация, используют
Такую аппроксимацию применяют в справочных расчетах потерь силовых полупроводниковых приборов; она же приведена в приложении C ГОСТ Р МЭК 61800-4-2012 для потерь в открытом тиристоре [4].
3. Расчет потерь и управляющего импульса
Средняя мощность потерь в открытом состоянии оценивается по среднему и действующему значениям тока:
Первое слагаемое описывает потери на условном пороговом падении напряжения, второе - джоулеву составляющую на наклонном сопротивлении ВАХ. Для первичной оценки берут значения при максимальной температуре перехода, потому что тепловой расчет должен быть консервативным.
Управляющая цепь должна обеспечить ток gate выше паспортного . В практических схемах закладывают запас, потому что зависит от температуры, разброса прибора и формы импульса. В этой методике используется проверка
Для источника , отпирающего напряжения и последовательных сопротивлений оценочный ток равен
Это выше . Поэтому в модели управляющий импульс подается через реальную электрическую ветвь с источником, ключом, трансформатором, диодом и резисторами, а не прямым логическим сигналом в gate: такой способ сохраняет физический смысл тока управления и согласуется с рекомендацией Engee учитывать суммарное сопротивление цепи gate [8].
4. Параметризация модели радиатора тиристора О123
Радиатор в модели параметризуется по паспорту. Важно не подменять радиатор обычным постоянным тепловым сопротивлением, если требуется динамическая тепловая картина: паспорт задает массу, материал и тепловые сопротивления, а блок Радиатор в Engee позволяет задать не только статический отвод тепла, но и тепловую емкость.
Основные паспортные данные для О123:
| Параметр | Обозначение | Значение |
|---|---|---|
| Номинальная рассеиваемая мощность | ||
| Тепловое сопротивление при естественном охлаждении | ||
| Тепловое сопротивление при принудительном охлаждении | ||
| Масса радиатора | ||
| Удельная теплоемкость алюминиевого сплава |
Динамические характеристики
Масса радиатора берется непосредственно из паспорта: . Материал указан как коррозионностойкий алюминиевый сплав, поэтому для учебно-инженерной модели принимается типовое значение
Тепловая емкость радиатора:
Этот параметр нужен для переходного теплового процесса: чем больше , тем медленнее растет температура радиатора при той же рассеиваемой мощности.
Опорная точка установившегося режима
Формула теплового сопротивления:
Отсюда перепад температур при номинальной мощности и естественной конвекции:
Эта точка используется как опорная для построения табличной зависимости блока Радиатор: , .
Степенная модель естественной конвекции
Для естественного воздушного охлаждения тепловой поток удобно аппроксимировать степенной зависимостью:
Для естественной конвекции принимается показатель . Коэффициент находится из опорной точки:
После этого мощность, отводимая радиатором при выбранном векторе температур, рассчитывается как
В модели используется вектор . Такой набор точек удобен для табличной параметризации: он покрывает малый, средний и близкий к номинальному перегрев радиатора, не перегружая модель лишними точками.
Альтернативная линейная оценка
Для контроля можно сравнить степенную модель с линейной:
Линейная модель проходит через начало координат и предполагает постоянное тепловое сопротивление. Степенная модель лучше отражает естественную конвекцию, где интенсивность теплоотдачи меняется с перегревом. Поэтому для блока Радиатор используется табличный вектор , рассчитанный по степенному закону, а линейная оценка остается контрольной.
#=
Параметризация радиатора О123.
Эта ячейка повторяет расчет параметров, которые вводятся в блок `Радиатор`:
масса, удельная теплоемкость, вектор перепадов температур и вектор отводимой мощности.
Степенная зависимость используется для естественной конвекции; линейная модель оставлена
как контрольная оценка по постоянному тепловому сопротивлению.
=#
using Printf
# Паспортные параметры О123
P_nominal = 120.0 # рассеиваемая мощность, Вт
Rth_natural = 0.71 # тепловое сопротивление при естественном охлаждении, °C/Вт
Rth_forced = 0.212 # тепловое сопротивление при принудительном охлаждении, °C/Вт
mass = 2.0 # масса радиатора, кг
cp_aluminum = 900.0 # удельная теплоемкость алюминиевого сплава, Дж/(кг*К)
println("Основные параметры радиатора О123:")
println("Масса: $mass кг")
println("Удельная теплоемкость: $cp_aluminum Дж/(кг*К)")
println("Номинальная мощность: $P_nominal Вт")
println("Тепловое сопротивление, естественное охлаждение: $Rth_natural °C/Вт")
println("Тепловое сопротивление, принудительное охлаждение: $Rth_forced °C/Вт")
# Динамическая часть модели: тепловая емкость задает инерцию нагрева.
heatsink_heat_capacity = mass * cp_aluminum
println()
println("=== ДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ===")
println("Масса радиатора: ", mass, " кг")
println("Удельная теплоемкость: ", cp_aluminum, " Дж/(кг*К)")
println("Тепловая емкость радиатора: ", heatsink_heat_capacity, " Дж/К")
# Опорная точка для естественной конвекции из паспортной мощности и Rth.
delta_T_nominal = P_nominal * Rth_natural
println()
println("=== РАСЧЕТ ОПОРНОЙ ТОЧКИ ===")
println("Номинальная мощность P = ", P_nominal, " Вт")
println("Тепловое сопротивление Rth = ", Rth_natural, " °C/Вт")
println("Перепад температур ΔT = P * Rth = ", round(delta_T_nominal, digits=2), " K")
println("Опорная точка: (ΔT = ", round(delta_T_nominal, digits=2), " K, Q = ", P_nominal, " Вт)")
# Для естественной конвекции используем степенной закон Q = k * ΔT^n.
n = 1.25
k = P_nominal / (delta_T_nominal^n)
println()
println("=== ПАРАМЕТРЫ СТЕПЕННОЙ ЗАВИСИМОСТИ ===")
println("Показатель степени n = ", n)
println("Коэффициент k = ", round(k, digits=4))
println("Формула: Q = ", round(k, digits=4), " * ΔT^", n)
# Вектор ΔT соответствует табличной параметризации радиатора в модели.
delta_T_vector = [10.0, 30.0, 50.0, 70.0, 90.0]
Q_vector = [k * (dT^n) for dT in delta_T_vector]
println()
println("=== ВЕКТОРЫ ДЛЯ ПАРАМЕТРИЗАЦИИ МОДЕЛИ ===")
println("Вектор разности температур ΔT, K:")
println(delta_T_vector)
println("Вектор отводимой мощности Q, Вт:")
println(round.(Q_vector, digits=1))
println()
println("Таблица значений:")
println("│ ΔT, K │ Q, Вт │")
println("├───────┼─────────┤")
for (i, dT) in enumerate(delta_T_vector)
@printf("│ %6.1f │ %7.1f │\n", dT, Q_vector[i])
end
# Проверка: в опорной точке степенная модель должна вернуть P_nominal.
Q_check = k * (delta_T_nominal^n)
error = abs(Q_check - P_nominal)
println()
println("=== ПРОВЕРКА РАСЧЕТА ===")
println("При ΔT = ", round(delta_T_nominal, digits=2), " K:")
println(" Расчетная мощность Q = ", round(Q_check, digits=2), " Вт")
println(" Номинальная мощность = ", P_nominal, " Вт")
println(" Погрешность = ", round(error, digits=4), " Вт")
println()
println("="^60)
println("ИТОГОВЫЕ ПАРАМЕТРЫ ДЛЯ МОДЕЛИ")
println("="^60)
println("Способ параметризации: табличная зависимость Q(ΔT)")
println("Тип конвекции: естественная")
println("Вектор ΔT, K: [", join(round.(delta_T_vector, digits=1), ", "), "]")
println("Вектор Q, Вт: [", join(round.(Q_vector, digits=1), ", "), "]")
println("Масса радиатора: ", mass, " кг")
println("Удельная теплоемкость: ", cp_aluminum, " Дж/(кг*К)")
# Линейная модель нужна только для сравнения с нелинейной естественной конвекцией.
Q_linear = [dT / Rth_natural for dT in delta_T_vector]
println()
println("=== ЛИНЕЙНАЯ МОДЕЛЬ ДЛЯ СРАВНЕНИЯ ===")
println("Тепловое сопротивление Rth = ", Rth_natural, " °C/Вт")
println("Вектор Q линейной модели, Вт:")
println(round.(Q_linear, digits=1))
println()
println("│ ΔT, K │ Q степенная, Вт │ Q линейная, Вт │")
println("├───────┼─────────────────┼────────────────┤")
for (i, dT) in enumerate(delta_T_vector)
@printf("│ %6.1f │ %15.1f │ %14.1f │\n", dT, Q_vector[i], Q_linear[i])
end
string(
"Cth=", round(heatsink_heat_capacity; digits=1),
" J/K; delta_T_nom=", round(delta_T_nominal; digits=2),
" K; Q_vector=", join(round.(Q_vector; digits=1), ","),
)
5. Расчет электрических и управляющих блоков
Источник сети задается амплитудой, а не действующим значением:
При и период равен
RL-нагрузка нужна для демонстрации режима с непрерывностью тока и фазовым сдвигом. Ее параметры оцениваются через индуктивное сопротивление, модуль импеданса и электромагнитную постоянную времени:
Для и получаем , , . Диод свободного хода нужен именно из-за этой постоянной времени: ток индуктивности не должен исчезать в момент закрывания тиристора.
RC-цепь имеет постоянную времени
Она не заменяет полноценный расчет перенапряжений, но задает демпфирующий путь для быстрых процессов. Тепловая масса воздуха рассчитывается по формуле
Этот может понадобиться для моделирования инерции окружающей среды полигона.
6. Система импульсно-фазового управления
СИФУ строят как синхронизированную задержку от перехода сетевого напряжения через ноль. Блок Hit Crossing выделяет начало полупериода, RS-триггер формирует логическое состояние ожидания, а Variable Time Delay переносит момент отпирания на величину, заданную углом управления.
Задание угла в модели:
Ограничитель Saturation задает допустимый диапазон . Для текущей синусоиды фактический диапазон , поэтому нижнее ограничение не искажает сигнал, но защищает модель при другой амплитуде задания. Связь угла и задержки:
При это дает диапазон задержек . Такая структура выбрана вместо генератора произвольных импульсов, потому что она сохраняет главное свойство фазового управления: момент открытия всегда отсчитывается от фазы питающего напряжения, а не от независимых часов.
7. Аналитический расчет перед запуском модели
Следующая кодовая ячейка выполняет контрольный расчет паспортных, электрических и тепловых величин. В ячейке проверяются амплитуда сети, потери тиристора, суммарные тепловые сопротивления, оценка температуры перехода, параметры RL-нагрузки, RC-цепи, диапазон задержки СИФУ и запас управляющего тока.
Эта ячейка не обращается к VT.engee; она нужна как инженерная проверка исходных чисел перед численным моделированием.
#=
Здесь собраны только аналитические оценки. Их цель - быстро проверить,
что паспортные данные и параметры блоков согласованы до запуска модели.
Сами формулы вынесены в текст выше, а ниже оставлены численные подстановки.
=#
# Расчет паспортных и модельных параметров
using Printf
f = 50.0
T = 1 / f
U_rms = 600.0
U_peak = sqrt(2) * U_rms
ITAVM = 200.0
ITRMS = 315.0
UT0 = 1.05
rT = 1.17e-3
P_TAV = UT0 * ITAVM + rT * ITRMS^2
Rth_jc = 0.080
Rth_ch = 0.020
Rth_ha_nat = 0.710
Rth_ha_forced = 0.212
Rth_ja_nat = Rth_jc + Rth_ch + Rth_ha_nat
Rth_ja_forced = Rth_jc + Rth_ch + Rth_ha_forced
T_ambient = 25.0
Tj_est_nat = T_ambient + P_TAV * Rth_ja_nat
Tj_est_forced = T_ambient + P_TAV * Rth_ja_forced
R_load = 6.0
L_load = 5e-3
X_L = 2*pi*f*L_load
Z_load = sqrt(R_load^2 + X_L^2)
tau_load = L_load / R_load
Rs = 100.0
Cs = 1e-4
tau_snubber = Rs * Cs
alpha_min = 90.0 - 82.0
alpha_max = 90.0 + 82.0
delay_min = alpha_min * T / 360
delay_max = alpha_max * T / 360
Vgate_source = 15.0
VGT = 3.5
IGT = 0.20
R2 = 7.0
Rg = 10.0
Ig_est = (Vgate_source - VGT) / (R2 + Rg)
rows = [
("Период сети T", T, "с"),
("Амплитуда сети Vm", U_peak, "В"),
("Потери тиристора при ITAVM/ITRMS", P_TAV, "Вт"),
("Rthja естественное", Rth_ja_nat, "К/Вт"),
("Rthja принудительное", Rth_ja_forced, "К/Вт"),
("Оценка Tj при естественном охлаждении", Tj_est_nat, "°C"),
("Оценка Tj при принудительном охлаждении", Tj_est_forced, "°C"),
("Индуктивное сопротивление нагрузки XL", X_L, "Ом"),
("Модуль импеданса нагрузки |Z|", Z_load, "Ом"),
("Постоянная времени нагрузки L/R", tau_load, "с"),
("Постоянная времени RC-цепи", tau_snubber, "с"),
("Минимальная задержка СИФУ", delay_min, "с"),
("Максимальная задержка СИФУ", delay_max, "с"),
("Оценка тока управления gate", Ig_est, "А"),
("Запас к 3*IGT", Ig_est / (3*IGT), "раз"),
]
for (name, value, unit) in rows
@printf("%-46s = %12.6g %s\n", name, value, unit)
end
string("U_peak=", round(U_peak; digits=2), " V; P_TAV=", round(P_TAV; digits=2), " W; Ig/(3IGT)=", round(Ig_est/(3IGT); digits=3))
8. Программный запуск модели и подготовка данных
Следующая кодовая ячейка загружает существующую модель VT.engee, временно ограничивает время расчета в открытой сессии, включает логирование нужных портов, запускает моделирование и преобразует результаты в массивы для графиков. Исходный файл модели не сохраняется.
Ограничение StopTime до введено из практических соображений. В модели задан фиксированный шаг , поэтому расчет исходных может дать десятки миллионов точек и перегрузить память при построении графиков. Для методической проверки формы сигналов достаточно одной секунды, а при тепловых исследованиях время можно увеличивать осознанно.
#=
Запуск модели.
Эта ячейка выполняет три задачи: готовит окружение,
запускает существующую модель Engee и приводит результаты к массивам.
=#
include("$(@__DIR__)/useful_functions.jl")
# Перед загрузкой закрываем одноименную модель, чтобы не смешивать результаты разных запусков.
try_close_model(model_name)
model_path = resolve_existing_file(model_file)
model = engee.load(model_path; force=true)
engee.open(model)
# Запоминаем исходные настройки только для отчета в статусной строке.
original_stop_time = safe_model_param(model, "StopTime", "")
original_fixed_step = safe_model_param(model, "FixedStep", "")
original_solver = safe_model_param(model, "SolverName", "")
status_message = ""
# Временное ограничение расчета в текущей сессии. Файл модели не сохраняется.
safe_set_model_param!(model, "StopTime", string(simulation_stop_time))
# Временная запись сигналов в текущей сессии. Исходный файл модели не сохраняется.
log_ports = [
"Генератор синусоиды/1",
"Датчик напряжения/1",
"Датчик тока/1",
"Смещение/1", # температура теплового узла после перевода из K в °C
]
for port in log_ports
safe_set_log(port)
end
# Запускаем модель после настройки логирования; результат Engee читается как набор именованных сигналов.
results = engee.run(model)
# Извлекаем три ключевых сигнала: угол открытия, напряжение нагрузки и ток нагрузки.
(angle_xy, angle_key) = get_signal(results, ["alpha", "Генератор синусоиды.1", "Генератор синусоиды.main_out", "Генератор синусоиды"])
(vload_xy, vload_key) = get_signal(results, ["Vload", "Датчик напряжения.Vload", "Датчик напряжения.1", "Датчик напряжения.v_out"])
(iload_xy, iload_key) = get_signal(results, ["Iload", "Датчик тока.Iload", "Датчик тока.1", "Датчик тока.i_out"])
# Прореживаем мгновенные сигналы до 1000 точек перед визуализацией.
t_alpha, alpha = thin(angle_xy[1], angle_xy[2]; max_points=max_plot_points)
t_v, vload = thin(vload_xy[1], vload_xy[2]; max_points=max_plot_points)
t_i, iload = thin(iload_xy[1], iload_xy[2]; max_points=max_plot_points)
# Считаем средние значения за один сетевой период и тоже прореживаем результат.
v_avg_full = moving_average_by_time(vload_xy[1], vload_xy[2], network_period)
i_avg_full = moving_average_by_time(iload_xy[1], iload_xy[2], network_period)
t_vavg, v_avg = thin(vload_xy[1], v_avg_full; max_points=max_plot_points)
t_iavg, i_avg = thin(iload_xy[1], i_avg_full; max_points=max_plot_points)
try_close_model(model_name)
9. Осциллограммы тока и напряжения
Следующая кодовая ячейка строит два базовых графика: мгновенное напряжение на нагрузке и мгновенный ток нагрузки. По этим осциллограммам проверяют, что тиристор открывается с фазовой задержкой, ток индуктивной нагрузки не обрывается мгновенно, а амплитуды находятся в ожидаемом диапазоне.
# Эта ячейка показывает мгновенные формы сигналов после фазового управления.
p1 = plot(t_v, [vload, iload]; label=["Vload" "Iload"], xlabel="t, с", ylabel="U, В; I, А", title="Напряжение и ток нагрузки", legend=:topright, size=(800, 400))
period = 600:800
p2 = plot(t_v[period], [vload[period], iload[period]]; label=["Vload" "Iload"], xlabel="t, с", ylabel="U, В; I, А", title="Напряжение и ток нагрузки", legend=:topright, size=(800, 400))
plot(p1, p2; layout=(2, 1), size=(800, 400))
10. Сравнение усредненных величин с углом открытия
Следующая кодовая ячейка сопоставляет усредненные за период сети напряжение и ток с сигналом угла открытия . Усреднение выполняется в предыдущей ячейке моделирования скользящим окном , то есть за один период сети . Такой график показывает не отдельные коммутационные фронты, а медленную регулировочную реакцию нагрузки на изменение угла.
# Здесь сравниваются средние за период величины с медленно меняющимся углом открытия.
p3 = plot(t_vavg, [v_avg, alpha]; label=["mean(Vload), окно 20 мс" "alpha"], xlabel="t, с", ylabel="Uavg, В; alpha, град", title="Усредненное напряжение", legend=:topright)
p4 = plot(t_iavg, [i_avg, alpha]; label=["mean(Iload), окно 20 мс" "alpha"], xlabel="t, с", ylabel="Iavg, А; alpha, град", title="Усредненный ток", legend=:topright)
plot(p3, p4; layout=(2, 1), size=(800, 400))
11. Регулировочная характеристика
Следующая кодовая ячейка строит зависимость усредненных выходных величин от угла открытия. Начальные точки отбрасываются индексом 40:end, чтобы уменьшить влияние переходного процесса в начале расчета. Такой график удобен для инженерного вывода: при увеличении среднее напряжение и ток управляемого выпрямителя уменьшаются, а при уменьшении угла - возрастают.
# График в координатах alpha показывает регулировочную характеристику нагрузки.
plot(alpha[40:end], [v_avg[40:end], i_avg[40:end]]; label=["Vload" "Iload"], xlabel="alpha, град", ylabel="Uavg, В; Iavg, А", title="Регулировочные характеристики", legend=:topright, size=(800, 400))
12. Температура корпуса тиристора
В тепловой части модели датчик абсолютной температуры подключен к тепловому узлу тиристора и радиатора, а блок Смещение переводит результат из Кельвинов в градусы Цельсия. Поэтому график ниже следует читать как температуру корпуса/контактного теплового узла, а не как температуру кристалла. Для оценки температуры перехода используется связь
Здесь берется из моделирования, - тепловая мощность, рассеиваемая тиристором, - паспортное тепловое сопротивление переход-корпус. Если требуется именно температура перехода, в модель добавляют соответствующий внутренний тепловой вывод или вычисляют ее после расчета через потери и .
Следующая кодовая ячейка извлекает температурный сигнал из результатов моделирования, ограничивает число точек до max_plot_points и строит две кривые: абсолютную температуру корпуса и перегрев относительно начального значения. На коротком интервале график обычно растет медленно, потому что радиатор О123 имеет большую тепловую емкость . Такой результат физически ожидаем: электрические процессы устанавливаются за миллисекунды, а тепловой процесс развивается на существенно большем времени.
# Эта ячейка извлекает и строит температуру корпуса/теплового узла тиристора.
# Список имен оставлен расширенным, потому что Engee может вернуть сигнал по имени линии,
# имени блока или автоматически сформированному имени порта.
(temp_xy, temp_key) = get_signal(results, [
"temp_VT_C",
"tempC",
"Смещение.temp_C",
"Смещение.1",
"Смещение.main_out",
"Датчик абсолютной температуры.1",
])
t_temp, temp_case = thin(temp_xy[1], temp_xy[2]; max_points=max_plot_points)
temp_rise = temp_case .- first(temp_case)
p_temp1 = plot(
t_temp,
temp_case;
label=temp_key,
xlabel="t, с",
ylabel="Tcase, °C",
title="Температура корпуса тиристора",
legend=:topright,
)
p_temp2 = plot(
t_temp,
temp_rise;
label="ΔTcase",
xlabel="t, с",
ylabel="ΔT, °C",
title="Перегрев корпуса относительно начального значения",
legend=:topright,
)
plot(p_temp1, p_temp2; layout=(2, 1), size=(800, 400))
Библиографический список и справочные источники
- Розанов Ю. К., Рябчицкий М. В., Кваснюк А. А. Силовая электроника: учебник для вузов. М.: Издательский дом МЭИ, 2016.
- Чиженко И. М., Руденко В. С., Сенько В. И. Основы преобразовательной техники.
- Герман-Галкин С. Г. Компьютерное моделирование полупроводниковых систем в MATLAB 6.0. СПб.: КОРОНА принт, 2001.
- ГОСТ Р МЭК 61800-4-2012. Системы силовых электроприводов с регулируемой скоростью.
- Гребенников В. В. Исследование характеристик тиристоров.
- SEMIKRON Danfoss. Application Manual Power Semiconductors.
- Engee. Программное управление скриптами.
- Engee. Thyristor.
- Engee. Радиатор.
- Паспорт
Тиристор T123-200 datasheet.pdf. - Паспорт
Охладитель О123 datasheet.pdf.