УДК 621.316.935
ИССЛЕДОВАНИЕ ВЛИЯНИЯ НАЧАЛЬНЫХ УСЛОВИЙ НА УСТАНОВИВШИЕСЯ ПРОЦЕССЫ В ФЕРРОРЕЗОНАНСНОМ КОНТУРЕ НА ОСНОВЕ РЕШЕНИЯ УРАВНЕНИЯ ДУФФИНГА
И.Е. Кажекин, кандидат технических наук, ФГБОУ ВО «Калининградский государственный технический
университет» e-mail: kazhekin@mail.ru
В статье рассмотрено описание феррорезонансного контура на основе уравнения Дуффинга. Приведены основные выражения, определяющие возможность возникновения феррорезонанса, также амплитуды колебаний в различных состояниях контура. Изложен подход к оценке влияния начальных условий на возможность возбуждения феррорезонансных колебаний в контуре.
феррорезонанс напряжений, нелинейные колебания, уравнение Дуффинга, перенапряжения, однофазные замыкания
Введение
Образование феррорезонансного контура в электросетях, вне зависимости от их напряжения и конфигурации, представляет большую опасность. Такие контуры могут формироваться при различных видах несимметрии в электрических сетях, сочетающих емкостные и нелинейные индуктивные элементы [1]. Один из примеров возникновения феррорезонанса описан в [2]. Нелинейный колебательный контур в рассматриваемом случае возникает при однофазных замыканиях в низковольтных судовых электроэнергетических системах с нейтралью, заземленной через реактор. Особенностью процессов является то, что контур образуется при контакте фазы с корпусом судна [3]. Поэтому на феррорезонансные процессы накладываются еще и дуговые, что значительно усложняет анализ процесса. Дуговые процессы во многом определяют начальные условия при очередном замыкании фазы на корпус и формировании феррорезонансного контура. Это приводит к тому, что перенапряжения в таких сетях могут значительно превосходить перенапряжения в сетях с изолированной нейтралью [4]. Учитывая, что способ заземления нейтрали является основой при построении защиты от однофазных замыканий и при его выборе необходимо учитывать все последствия применения того или иного варианта [5], весьма важным является исследование феррорезонансного контура, формирующегося при заземлении нейтрали через реактор.
Несмотря на то, что исследованию колебательных контуров с нелинейной индуктивностью посвящено достаточно много работ, в настоящее время остаются трудности при его математическом описании. Эти трудности обусловлены тем, что использование более точного описания усложняет анализ процессов и оценку влияния на них различных факторов. Упрощение описания, как правило, ведет к снижению точности расчетов процессов в контуре. В статье предлагается подход к оценке возможности возникновения феррорезонанса на основе решения уравнения Дуффинга, что с одной стороны позволяет несколько повысить точность расчетов, по сравнению, например с [1], с другой - получить выражения, удобные для исследования влияния различных факторов на характер колебаний в нелинейном контуре.
Описание феррорезонансного контура
Схема феррорезонансного контура, образующегося в рассматриваемом случае, приведена на рисунке.
Рисунок - Схема феррорезонансного контура Изображенная схема может быть описана следующим уравнением [6]:
Жщ + Я ■ г + — ГгЛ = Еот $>т(Ш + ф) С
Жг
(1)
При этом для описания нелинейного индуктивного элемента довольно часто предлагается выражение на основе кубического полинома:
г = ащ + р щ3
(2)
Из (1) и (2) можно получить выражение:
Ж 2щ
+ Я ■«■
9 \
1 3 Вщ Жщ а Р 3 / ч
1 + --Ж~ + ~СЩ+С.Щ = к ^ Ет*соЩог + ф)
а
(3)
Используя обозначения а0 =-^а/ С , щ0=а^ Ет0/а1, а = ю/ю0, х = щ/щ0, т = а0 ■ /,
Л = Р1а , щ0={р/а\{а^ Ет02 5 = Я ■ а/®0, выражение (3) можно привести к следующему виду:
Ж х
■ +
2 ■ 5 ■ (1 + 3 ■ Л ■ х2) — + х + Л ■ х3 = к ■ со$>(а ■т + ф)
Жт
(4)
Если в выражении (4) не учитывать 3 ■ Л ■ х2, то оно может быть записано в виде уравнения Дуффинга:
+ 2 ■ 5 — + х + Л ■ х3 = к ■ со$>(а ■т + ф) Жт Жт
(5)
Уравнение (5) достаточно подробно исследовано и может быть использовано для определения условий возникновения феррорезонанса.
Решение уравнения Дуффинга
В литературе предлагаются различные подходы к решению уравнения Дуффинга. Наиболее простым и точным является использование численных методов, например, Булирша-Штера, Розенброка и др., встроенных в ряд программных сред [7]. Однако такое решение значительно затрудняет дальнейшее исследование влияния различных факторов на установившиеся колебания.
Использование же упрощенных методов расчета приводит к снижению точности, однако одновременно с этим упрощается процесс исследования влияния начальных условий и параметров феррорезонансного контура на величину колебаний. Кроме того, при использовании упрощенных методов появляется возможность оценки амплитуды, соответствующей неустойчивому состоянию равновесия.
Один из таких упрощенных подходов к решению уравнения Дуффинга изложен в [8]. Согласно этому подходу число корней уравнения может быть определено по следующему выражению:
(1, при с < 1;
п = 1. , (6) [3, при с > 1,
где с - коэффициент, определяемый по выражению (7):
с = 4 ■(„' - 0 (7)
3 ■ 3/12-А-к
Выражение (6) формулирует условия формирования феррорезонансного контура. При сочетании параметров цепи, соответствующим с < 1, можно сделать вывод, что феррорезонанс ислючен, при с > 1 появляется возможность для его возникновения.
При этом амплитуда колебаний в феррорезонансном контуре может быть найдена следующим образом:
а„ = (8)
где ё - коэффициент, определяемый в результате решения кубического уравнения:
а3 - з ■ с ■ а+2 = о (9)
Уравнение (9) при с > 1 имеет три корня, два из которых соответствуют большим и малым колебаниям или резонансному и нерезонансному состоянию рассматриваемого контура.
Влияние начальных условий на установившиеся колебания в феррорезонансном контуре
Состояние равновесия, в которое перейдет система при с> 1 после образования феррорезонансного контура, определяется начальными условиями. Устойчивые состояния равновесия образуют области притяжения. Любые начальные условия, соответствующие точкам одной из областей притяжения, приводят к определенному виду колебаний. Области притяжения различных видов колебаний могут быть найдены следующим образом [9]:
1. Определяются все решения уравнения.
2. Оценивается устойчивость найденных решений.
3. Границей между областями притяжения является кривая, представляющая собой геометрическое место образов, приближающихся к неустойчивому состоянию системы.
Таким образом, оценка достаточности начальных условий для возбуждения феррорезонансных колебаний в контуре может заключаться в их сопоставлении с амплитудой колебаний, найденной в соответствии с (8) для неустойчивого состояния равновесия. Если начальные условия превышают амплитуду неустойчивых колебаний, то в установившемся режиме контур окажется в состоянии феррорезонанса напряжений.
3
Предложенное решение имеет ряд допущений. В частности, не учитывается возможность возникновения колебаний на субгармониках и на гармониках, кратных трем. Точное определение амплитуды колебаний возможно только для нерезонансного режима. Для неустойчивых и резонансных колебаний необходимо учитывать весь возможный спектр гармоник. В то время как для цепи, изображенной на рисунке, в зависимости от ее параметров и начальных условий, могут быть получены следующие виды феррорезонанса [10]:
- феррорезонанс на основной частоте, сопровождающийся появлением высших гармонических составляющих;
- субгармонический феррорезонанс;
- квазипериодический феррорезонанс, характеризующийся широким спектром гармонических составляющих тока контура;
- хаотический феррорезонанс; режим, при котором возникает непрерывный спектр гармоник.
Однако для оценки возможности возбуждения ферррорезонанса достаточно иметь представления о границах области притяжения нерезонансных колебаний. Если до образования феррорезонансного контура на элементах сформировались начальные значения переменных, выходящие за эти границы, то включение контура приведет к переходу системы в один из видов феррорезонанса.
Заключение
Установленные соответствия между описанием колебательного контура, содержащего нелинейный индуктивный элемент, и решением уравнения Дуффинга позволяют провести оценку возможности возникновения феррорезонанса напряжений исходя из параметров элементов контура. Помимо этого, при установлении возможности возникновения феррорезонансных колебаний в нелинейной цепи на основе приведенных выражений можно оценить величины и соотношения начальных условий, при которых система окажется в состоянии феррорезонанса.
Описанные в работе исследования выполнены при финансовой поддержке РФФИ в рамках научного проекта № 18-38-00622.
СПИСОК ЛИТЕРАТУРЫ
1. Гиндуллин, Ф.А. Перенапряжения в сетях 6-35 кВ / В.Г. Гольдштейн, А.А. Дульзон, Ф.Х. Халилов. - Москва: Энергоатомиздат, 1989. - 192 с.
2. Кажекин, И.Е. Феррорезонансные процессы при однофазных замыканиях в судовых электроустановках с компенсированной нейтралью / И.Е. Кажекин // Известия КГТУ. - 2019. - № 52. - С. 145-153.
3. Благинин, В.А. Описание феррорезонансных процессов при однофазных замыканиях в судовых электроэнергетических системах с компенсированной нейтралью / В.А. Благинин, В.И. Лозовенко, И.Е. Кажекин // Инновации в науке и образовании - 2011: IX Междунар. науч. конф.: в 2-х частях: труды. - 2011. - С. 427-430.
4. Благинин, В.А. Дуговые перенапряжения в судовых электросистемах с различными режимами нейтрали/ В.А. Благинин, В.И. Лозовенко, И.Е. Кажекин // Известия КГТУ. - 2008. - № 13. - С. 118-121.
5. Кажекин, И.Е. Определение требований к защите судовых электроэнергетических систем от однофазных замыканий / И.Е. Кажекин // Электрика. - 2011. - № 3. - С. 35-38.
6. Филиппов, Е. Нелинейная электротехника / Перевод [с нем.] и доп. А. З. Кулебякина. - Москва: Энергия, 1968. - 503 с.
7. Гурский, Д.А. Вычисления в MathCad 12 / Д.А. Гурский, Е.С. Турбина. - Санкт-Петербург: Питер. - 2006. - 544 с.
8. Баев, С.В. Построение аналитического решения уравнения дуффинга / С.В. Баев, Ю.С. Пикуш // Вюник Придншровсько'1 державно'1 академп будiвництва та архггектури. -2009. - № 6-7 (137). - С. 43-48.
9. Хаяси, Т. Нелинейные колебания в физических системах / Т. Хаяси; под ред. В. Е. Боголюбова; пер. с англ. Б. А. Болдова, Г. Г. Гусева. - Москва: Мир, 1968. - 432 с.
10. Ferracci P. Ferroresonance / P. Ferracci // Cahier technique. - Paris: Schneider Electric, 1998. - No. 190. - P. 1 - 28.
RESEARCH THE EFFECT OF PRIMARY CONDITIONS ON STATIONARY PROCESSES IN A FERROSONANCE CIRCUIT BASED ON SOLUTION DUFFING EQUATION
I.E. Kazhekin, Associate Professor, Kaliningrad State Technical University e-mail: kazhekin@mail.ru
The article describes the description of the ferroresonance circuit based on the Duffing equation. The main expressions that determine the possibility of ferroresonance, as well as the amplitudes of oscillations in different states of the circuit are given. The approach to estimating the influence of initial conditions on the possibility of exciting ferroresonant oscillations in the circuit is described.
ferroresonance, nonlinear oscillations, Duffing equation, overvoltage, single-phase circuit