УДК 621.3.062.8
МОДЕЛИРОВАНИЕ РАБОТЫ РЕЗОНАНСНОГО КОММУТАТОРА В ЦЕПИ ПОСТОЯННОГО ТОКА
С.В. Пустынников, Т.Е. Хохлова
Томский политехнический университет E-mail: [email protected]; [email protected]
Исследована работа резонансного коммутатора в сильноточной цепи постоянного тока, позволяющего при помощи резонанса напряжений осуществлять бесконтактное отключение индуктивной нагрузки. Получено расчетное выражение для резонансной емкости. Разработана математическая модель расчета переходного процесса методом переменных состояния, позволяющая определить ток в цепи нагрузки по параметрам резонансного коммутатора. Результаты расчета подтверждены виртуальным моделированием.
Ключевые слова:
Сильноточные и слаботочные цепи, ток нагрузки, резонанс напряжений, резонансный коммутатор, переходный процесс, метод переменных состояния.
Key words:
Low-current and high-current circuits, the load current, the resonance in the series circuit, resonant switch, the transient process, the method of state variables.
Индуктивная нагрузка является инерционным элементом, а ток через индуктивность не может меняться скачком, что приводит к плавному включению индуктивных нагрузок. Коммутация сильноточных цепей постоянного тока, содержащих индуктивную нагрузку или имеющих внутреннюю индуктивность - линий электропередач, линий связи, цепей с генераторами и двигателями постоянного тока и т. д., осуществляется при помощи электромеханических устройств - пускателей, контакторов [1].
Для контактов реле включение индуктивной нагрузки - простая задача, потому что ток через контакты растет медленно. При скачкообразном выключении напряжения на нагрузке запасенная в индуктивности энергия ищет выход, порождая напряжение самоиндукции, по полярности обратное напряжению источника питания. В [2, 3] показано, что размыкание цепей постоянного тока с индуктивностью за время Л^0 приводит к изменению потокосцепления индуктивности от начального значения ^=Ы0 до нуля. При этом теоретически в индуктивности возникает импульс перенапряжения Ы1=№/& бесконечной величины.
На практике Л1>0, что сопровождается возникновением дуги на размыкающих контактах, а также скачком напряжения на индуктивности в 5...7 раз превышающем напряжение источника питания, что приводит к быстрому износу коммутирующего оборудования. При отключении индуктивностей выделяется запасенная в катушке энергия, создавая напряжение самоиндукции, равное рабочему напряжению, умноженному на добротность катушки. Добротность индуктивной нагрузки бывает на практике от 0,5 (катушки с большим внутренним сопротивлением) до 50 (типичные соленоиды электромеханических замков, катушки контакторов и мощных реле, электромоторы и т. д.). Напряжение самоиндукции катушки обычного промышленного реле с рабочим напряжением 24 В может превышать 1 кВ.
Авторами была разработана и исследована модель резонансного коммутатора сильноточной цепи постоянного тока с индуктивной нагрузкой, позволяющая при помощи резонанса напряжений осуществлять бесконтактное отключение индуктивной нагрузки от источника постоянного напряжения. Схема предложенной модели показана на рис. 1.
Рис. 1. Схема резонансного коммутатора
Модель состоит из:
• сильноточной цепи, в которой последовательно включены постоянный источник ЭДС Е1, активно-индуктивное сопротивление нагрузки ЯН, ЬН, тиристор VS;
• резонансного коммутатора, содержащего две индуктивно-связанных катушки индуктивности - Яь Ь1 и Яъ Ь2, включенных встречно (Ь1<Ь2 и Л1<Л2), подключенную последовательно с катушкой индуктивности Ь2 емкость С, ключ £ и источник синусоидальной ЭДС е2. Предложенная модель позволяет осуществлять
бесконтактное размыкание сильноточной цепи путем замыкания цепи резонансного коммутатора.
Подключение нагрузки ЯН, ЬНкисточнику ЭДС Е1 осуществляется путем подачи управляющего сигнала на тиристор VS; по силовой цепи протекает ток ц.
Для отключения нагрузки ЯН, ЬН от ЭДС Е; замыкают ключ Я, и резонансный коммутатор подключается к источнику синусоидальной ЭДС е2, под действием которой в цепи протекает ток /2, по величине в 3...5 раз меньше тока г\. Контур с катушкой индуктивности Л2, Ь2 и емкостью С настроен на резонанс напряжений, при заданной частоте. Тогда на индуктивности Ь2 возникает напряжение, превышающее величину ЭДС е2 в число раз равное величине добротности резонансного контура.
За счет взаимной индуктивности М, при встречном включении катушек, в первой катушке индуктивности Яъ Ь формируется отрицательный импульс тока, проходящий через диод УБ, ЭДС Е; и тиристор УЯ, минуя нагрузку ЛН, ЬН. Ток переходит через нулевое значение, что приводит к запиранию тиристора УЯ и отключению силовой цепи от источника ЭДС Е;. Энергия, запасенная в индуктивности нагрузки ЬН, шунтируется диодом УБ, поэтому перенапряжение на тиристоре УЯ не возникает.
Расчет резонансной емкости С можно произвести при помощи схемы на рис 2, учитывающей индуктивную развязку трансформатора и шунтирование нагрузки ЛН, ЬН.
Я, ]ХЫ -]ХМ
Ч 7ХМ
Ч
Я 2 7Хь2 -7ХМ -]Хс~-
7
^вх
Рис. 2. Схема для расчета резонансной емкости
Комплексное входное сопротивление схемы на рис. 2 равно:
7вх = Я2 + 7 (ХЬ2 - ХМ ) - }ХС +
+ 7хм (Я + 7 (Ху - хм )) = Я Х
ТУ • V вх J вх '
Я, + 7Х 11
По условию резонанса 1вх=В^ех, Хех=0, отсюда получаем:
Хс = Х12 2 - Хм +
Хм%+Х^1 Хм - Хь 1 ХМ
Я + Х^1
Выражение для резонансной ёмкости:
= Я,2 + о21\
рю о2Ь2Я,2 +ю4(Ь2Ь\ -ь,М2)’
где о (рад/с) - угловая частота ЭДС е2.
При расчете и проектировании резонансного коммутатора необходимо знать условия, при которых обеспечивается переход тока нагрузки через нулевое значение, что является обязательным для запирания тиристора У5. Поскольку в первый мо-
мент времени ¡=0+ после замыкания ключа Я нагрузка шунтируется диодом УБ, то параметры нагрузки ЛН, ЬН не оказывают влияния на переходный процесс в остальной цепи. С момента времени ¡=0+ переходный процесс описывается системой уравнений, составленной по законам Кирхгофа для мгновенных значений:
. С12
Л
¿2 %2
&
&
и с е2 .
Си = ¿2 Л
Полученную систему дифференциальных уравнений первого порядка можно решить непосредственно, не преобразуя ее в дифференциальную систему порядка п. Для этого используем метод переменных состояния, согласно которому в матричной форме дифференциальные уравнения переменных состояния записываются в виде [4]:
(Сх / С) = (А)(х) + (В )(Г), где матрица-столбец:
• переменных состояния (величины, подчиняющиеся законам коммутации) - (х);
• источников энергии - (У);
• первых производных по времени от переменных состояния - (аХ/с?0;
матрица:
• состояния (коэффициенты при переменных состояния) - (А);
• коэффициентов при источниках - (В).
Таким образом, получим:
' Сг, / Л ^ ( г, ^ (Е, ^
сИ2 / С
сис / с J
= (А)
Уис
+ (В)
(1)
где
(А) =
¿2 Я - ь2 ь, + М2
МЯ2
М
МЯ,
- ь2 ь + м 2 0
(
(В) =
-ь2
- ь2 ь, + М2
¿1 • я,
- ь2 ь + м 2
_1_
с
-М
-ь2 ь, + М2
ь1
-ь2 ь, + М2 0
-ь2 ь, + М2
- М
-ь2 ь + м 2
- ь2 ь, + М2
-ь1
- ь2 ь, + М2
Кроме того, для решения необходимо задать матрицу-столбец начальных условий, значения которой равны установившимся значениям токов /ь /2 и ?7С:
' Ei / Rh + R ^ '90,909N
(x0) = 0 • (x0)= 0
. 0 ; V 0 ,
Аналитическое решение уравнений состояния записывается в общем виде:
(^ = (+ (xпp)>
где (xсв)=(x0)e(A)' - определяет реакцию цепи, обусловленную ненулевыми начальными условиями (при отсутствии внешнего воздействия),
.( AXt -г)
(B)(V)dr - реакция цепи от внеш-
них воздействий при нулевых начальных условиях, e(A)' - матричная экспоненциальная функция, представляющая собой медленно сходящийся ряд.
Определение коэффициентов этого ряда - трудоемкий и громоздкий процесс, поэтому наиболее часто решение уравнений состояния выполняют численно, используя известные алгоритмы численного интегрирования.
Исследуем работу модели резонансного коммутатора с параметрами: ^=100 В, Д=0,001 Гн, Д=0,1 Ом, Ж=0,01Гн, Х2=0,1 Гн, Л2=0,Юм, 6^=101,5 мкФ, ^=1 Гн, ЛН=1 Ом, e2(t)=100V2sin(314t+180o) В.
Получили дифференциальные уравнения переменных состояния (1), матрицы состояния и коэффициентов при источниках:
(A) =
227.273
27.273 0
27.273
2.273 908,265
' -2273 (B) = -272,727
V 0
и матрицу начальных условий
-272,727
22,727
0
272.727 ^
22.727 0
0 ^
0 0
Для решения системы используем численный метод Рунге-Кутта, основанный на разложении переменных состояния в ряд Тейлора и учете его первых двух членов:
(xk+i) =(xk) + h(dxld \ где (xk) - матрица-столбец переменных состояния для k-го шага; h - шаг интегрирования по времени.
Методическая погрешность данного алгоритма пропорциональна ^и, естественно, уменьшение шага приводит к увеличению точности. Однако с ростом числа шагов погрешность может увеличиваться за счет погрешности округления [4].
На рис. 3 показаны результаты математического моделирования, проведенные в пакете программ MathCad 13.0.
По результатам моделирования получили минимальное значение тока нагрузки /min«-22,574 A.
На рис. 4 представлены результаты виртуального исследования работы резонансного коммутатора с помощью программы Electronic Workbench 5.12. Для моделирования индуктивно-связанных катушек использовался воздушный трансформатор.
Анализ осциллограммы показывает, что минимальное значение тока нагрузки составляет imin~—22,942 A, что практически совпадает с результатами математического моделирования.
Выводы
1. Разработан резонансный коммутатор, который состоит из емкости и индуктивно-связанных катушек, имеющих встречное включение. При этом катушка с меньшей индуктивностью включена в цепь с током нагрузки, а емкость и катушка с большей индуктивностью - к до-
Transformer Model 'ideal'
OK
Sheet 1 I
Primary-to-secondary turns ratio IN): I10
Leakage inductance (LE): |0.01 H
Magnetizing inductance (LM): I0-1 H
Primary winding resistance (RP): |0
Secondary winding resistance 1RS j: |o
Отмена
Oscilloscope
VA1
VB1
59 . 5548 9 . 1049
VA2
VB2
62 . 1 298 -2.2942
T2-T1
VA2-VAI
VB2-VB1
2.5750 - 1 1 .399 1
|l .OOms/div Bdae .f M I 5 V/Div I SkV/Div Reduce
X position j 0.00 - Level I 0.00 Y position I 0.00 - Y position I 0.00 - Reverse
ИЯ W № AC I 0 ПИ AC |Q DC I Save j
Рис. 4. Виртуальное моделирование работы резонансного коммутатора: а) схема; б) окно настройки воздушного трансформатора; в) осциллограмма тока в нагрузке
полнительному источнику синусоидальной ЭДС, образующие последовательный резонансный контур. Это позволяет путем замыкания цепи резонансного коммутатора формировать отрицательный импульс тока в нагрузке, приводящий к запиранию тиристора и отключению силовой цепи от источника ЭДС.
2. Разработана математическая модель расчета переходного процесса методом переменных состояния, позволяющая определить ток в цепи
нагрузки по параметрам резонансного коммутатора.
3. Проведено исследование работы модели резонансного коммутатора с заданными параметрами. Установлено, что подключение резонансного коммутатора приводит к переходу тока нагрузки через нулевое значение. Проведено виртуальное исследование работы резонансного коммутатора, которое подтвердило результаты математического моделирования.
СПИСОК ЛИТЕРАТУРЫ
1. Галинский П.П., Кононов А.П. и др. Теоретические основы электротехники в задачах. - Томск: Изд-во Том. ун-та, 1975. -294 с.
2. Резонансный коммутатор цепей постоянного тока: пат. на ПМ 67799 Рос. Федерация. 2007119003/22; заявл. 21.05.07; опубл. 27.10.07, Бюл. № 30. - 3 с.
3. Пустынников С.В., Хохлова Т.Е., Макенова Н.А. Использование индуктивного размыкателя для коммутации сильноточных цепей постоянного тока // Известия Томского политехнического университета. - 2009. - Т. 314. - № 4. - С. 88-91.
4. Зевеке Г.В., Ионкин П.А. и др. Основы теории цепей. - М.: Энергоатомиздат, 1989. - 528 с.
Поступила 28.10.2010 г.