- •2 Математическое и физическое моделирование в электротехнике и электронике
- •2.1 Математическое моделирование электрических цепей и процессов
- •2.1.1 Общие сведения о моделировании электрических цепей
- •2.1.1.1 Краткий обзор и классификация моделирующих программ
- •2.1.1.2 Общие сведения о вычислителях моделирующих программ
- •2.1.1.3 Понятие о структурном моделировании
- •2.1.1.4 Методология математического моделирования цепей
- •2.1.2 Применение основных методов и их машинная реализация
- •2.1.2.1 Применение метода переменных состояния
- •Порядок получения упс в электрических цепях
- •2.1.2.2 Применение узлового анализа
- •2.1.2.3 Модели цепи в базисе расширенных узловых уравнений
- •2.1.2.3.1 Идеальные приемники электрической энергии
- •2.1.2.3.2 Идеальные источники электрической энергии
- •2.1.2.3.3. Управляемые источники электрической энергии
- •2.1.2.3.4 Идеальные электрические ключи
- •2.1.2.3.5 Многополюсник (макромодель)
- •2.1.3 Моделирование динамических цепей активными rc-цепями
- •2.1.3.1 Формирование уравнений для построения модели rc-цепи
- •2.1.3.2 Построение модели rc-цепи
- •2.1.3.3 Исследование модели в частотной и временной области
- •2.1.4 Математическое моделирование режимов силовых цепей эту
- •2.1.4.1 Особенности компенсации реактивной мощности эту
- •2.1.4.2 Построение математических моделей силовых цепей
- •2.1.4.3 Моделирование динамических режимов силовых цепей эту
- •2.1.4.4 Моделирование установившихся режимов силовых цепей эту
2.1.4.3 Моделирование динамических режимов силовых цепей эту
В динамическом режиме дифференциальные уравнения индуктивностей и емкостей при анализе электрических схем аппроксимированы дискретными резистивными моделями, для неявного метода численного интегрирования. При построении программного комплекса выбран алгоритм метода трапеций.
Математическая модель динамической индуктивности реализованная в указанном методе интегрирования показана ниже.
, (2.24)
где h – шаг интегрирования во времени.
Для емкостного элемента в машинном моделировании также применена дискретная резистивная схемная модель и используется соответствующее выражение численного итерационного процесса:
. (2.25)
В приведенных дискретных моделях параметры режима изменяются от одного временного шага tn к другому tn+1. Использование дискретных моделей по сути сводит анализ переходной характеристики динамического режима к анализу резистивной цепи по постоянному току.
Пример моделирования переходных характеристик при пуске индукционной установки периодического действия показан на рис. 2.2.114, 2.2.115. Моделирование выполнено для схемы индукционной установки с детализированными схемными моделями ВДТ и БК. Исследованы переходные режимы при различных начальных условиях, соответствующих установившемуся режиму в начале нагрева заготовки. Переходные характеристики в схеме силовой цепи получены непосредственно для переменных состояния, в качестве которых приняты ток индуктора и напряжения на батарее конденсаторов.
|
|
|
|
|
|
Рис. 2.114 - Фазовая траектория переменной состояния il(t) для индуктора (а) и ее фазовые портреты на плоскости переменных состояния (б) |
|
|
|
При многовариантных исследованиях на ЭВМ пусковых режимов УИН с применением рассмотренных моделей установлено, что при всех значениях резистивных и реактивных элементов, определенных в электромагнитной задаче переходный процесс в УИН моделей ОКБ-84А и ИНМ-75, при мощностях до 295 кВт не приобретает продолжительного колебательного характера и завершается за 2 – 3 периода напряжения промышленной частоты. Коммутационные перенапряжения индукционных нагревателей не превышают значения 1000 В при всех значениях температуры и всех типоразмерах нагреваемых заготовок, а токи находятся в допустимых пределах, предусмотренных режимами эксплуатации коммутационной аппаратуры.
Кроме того, в ходе моделирования режимов силовых цепей однофазного и трехфазного индукционных нагревателей проведено исследование величин тока при коммутациях, связанных с переключением секций конденсаторных батарей и ступеней регулирования вольтодобавочных трансформаторов.
|
|
|
|
|
|
Рис. 2.115 - Фазовая траектория переменной состояния uC(t) для БК (а) и ее фазовые портреты на плоскости переменных состояния (б) |
|
Даже при нулевой начальной фазе коммутации источника для выбранных с применением оптимизационных процедур параметров регуляторов, величина тока в источнике питания не достигает ударного значения.
