Relaxation oscillations in the Darie wind power plant model

Full Text

1. Предварительные сведения
В данной статье проводится математическое моделирование колебаний ветроэнергетической установ-
ки Дарье (ВЭУ). В статьях [1–6] рассматривалась математическая модель ветрогенератора, состоящая
из двух уравнений. В [1; 2] были предложены новые типы конструкции ветроэнергетической установки
(рис 1.1.). Исследуются движения электромеханической системы, моделирующей стационарные режимы
генератора с ветротурбинным движителем с помощью качественных методов теоретической механики в
32
Кирсанова А.С. Релаксационные колебания в модели ветроэнергетической установки Дарье
Kirsanova A.S Relaxation oscillations in the Darie wind power plant model
статьях [3–5]. Рассмотрены вопросы существования и устойчивости установившихся режимов и описа-
ны области их притяжения [3–5]. Анализ критических режимов, основанный на использовании техники
траекторий уток, был проведен в [7]. В работе [6] исследовано уравнение стационарных режимов при раз-
личных значениях внешнего сопротивления динамической модели. Также описаны различные аспекты
конструктивной теории малых ветроэлектростанций [8; 9]. Математическая модель малой ветроэнерге-
тической установки Дарье представляет собой трехтемповую дифференциальную систему:
J ˙Ω = −kI +M(Ω);
LI˙ = kΩ − (R + r)I;
˙R
= "F (Ω; I;R);
(1.1)
где J — момент инерции турбины, L — индуктивность якоря, I — сила тока в обмотке якоря, Ω —
коэффициент быстроходности, R — внешнее сопротивление, M — относительный момент аэродинами-
ческих сил. Параметр k — коэффициент электромеханического взаимодействия, r — малое внутреннее
сопротивление якоря. Точка означает производную по времени t. График экспериментальных данных
на рис. 2.1 задает функцию M(Ω) [3]. Величины L и " рассматриваются как малые параметры. Это
означает, что в данной системе переменная I является самой быстрой, R — самой медленной из трех
переменных.
Рис. 1.1. Ветрогенератор Дарье
Fig. 1.1. Wind turbine Darrieus
2. Основные результаты
2.1. Аппроксимация функции M(Ω)
Рассмотрим график экспериментальных данных M(Ω), изображенный на рис. 2.1. Относительный мо-
мент аэродинамических сил и коэффициент быстроходности на рис. 2.1 являются безразмерными (M =
= Ma
0;5SbV 2 , Ω = b!
V , где ! — угловая скорость, b — расстояние от эффективного давления лопастей для
оси вращения, V — скорость воздуха, Ma — момент аэродинамических сил, — плотность воздуха,
S — площадь лопастей) [3; 5; 7]. Для построения функции M(Ω) зададим ее как кусочно-линейную
функцию по промежуткам:
Ω ∈ [0; 0:86]; [0:86; 1:8]; [1:8; 4]; [4; 8:64]
и аппроксимируем на каждом из участков отдельно.
Вестник Самарского университета. Естественнонаучная серия 2024. Том 30, № 1. С. 31–39
Vestnik of Samara University. Natural Science Series 2024, vol. 30, no. 1, pp. 31–39 33
Рис. 2.1. Экспериментальные данные
Fig. 2.1. Experimental data
В результате получается следующее задание функции M(Ω):
M(Ω) =
8>>>>>>>>>><
>>>>>>>>>>:
0:02906976745 Ω + 0:12;
где Ω ∈ [0; 0:86];
−0:0638297872 Ω + 0:199893617;
где Ω ∈ [0:86; 1:8];
0:643181818 Ω − 1:0727272725;
где Ω ∈ [1:8; 4];
−0:28556034485 Ω + 2:6422413795;
где Ω ∈ [4; 8:64]:
(2.1)
Рисунок 2.2 показывает график функции M(Ω), заданной (2.1).
Рис. 2.2. График функции M(Ω), заданной (2.1)
Fig. 2.2. Graph of the function M(Ω), given by (2.1)
В данном случае применяется достаточно грубая аппроксимация, однако важно отметить, что более
точная аппроксимация относительного момента аэродинамических сил не влияет существенным образом
на качественное поведение решений дифференциальной системы.
34
Кирсанова А.С. Релаксационные колебания в модели ветроэнергетической установки Дарье
Kirsanova A.S Relaxation oscillations in the Darie wind power plant model
2.2. Медленная кривая
Система (1.1) является сингулярно возмущенной, так как малый параметр L умножается на произ-
водную силы тока. Система (1.1) имеет двумерное инвариантное многообразие [7; 10; 11]
I = h(Ω; I;R): (2.2)
Функция h может быть найдена из уравнения инвариантности
L
@h
@R
"F (Ω; I;R) + L
@h
@Ω
1
J
(−kh +M(Ω)) = kΩ − (R + r)h: (2.3)
Пренебрегая членами порядка O(L), получается
h =
kΩ
R + r
: (2.4)
Движение по инвариантному многообразию описывается дифференциальной системой, где k = 0:5; r =
= 0:1; J = 1:
J ˙Ω = − k2Ω
R+r +M(Ω);
˙R
= "F (Ω; I;R):
(2.5)
Полученная система является быстро-медленной, так как правая часть уравнения для медленной
переменной R умножается на малый параметр ". В связи с этим мы можем использовать аппарат теории
релаксационных колебаний для анализа этой системы [12; 13]. Медленная кривая системы (2.5) задается
уравнением:
− k2Ω
R + r
+M(Ω) = 0: (2.6)
Из (2.6) получим уравнение медленной кривой в явной форме, для этого выразим R(Ω):
R = R(Ω) =
k2Ω
M(Ω)
− r: (2.7)
График медленной кривой представлен на рис. 2.3, при k = 0:5; r = 0:1; J = 1.
Рис. 2.3. Графики медленной кривой
Fig. 2.3. Slow curve charts
Для исследования устойчивости медленной кривой (2.6) необходимо найти производную функции
R(Ω) на каждом из промежутков отдельно:
@R(Ω)
@Ω
=
k2M(Ω) −M′(Ω)k2Ω
M2(Ω)
: (2.8)
На первом участке медленной кривой при Ω ∈ [0; 0:86]: пусть Ω = 0:5;R = 0:8291270529, тогда
@R(Ω)
@Ω (0:5; 0:8291270529) = −0:2400000001 < 0, следовательно, участок является устойчивым.
На втором участке медленной кривой при Ω ∈ [0:86; 1:8]: пусть Ω = 1;R = 1:737373042, тогда
@R(Ω)
@Ω (1; 1:737373042) = −0:1998936100 < 0, следовательно, участок является устойчивым.
Вестник Самарского университета. Естественнонаучная серия 2024. Том 30, № 1. С. 31–39
Vestnik of Samara University. Natural Science Series 2024, vol. 30, no. 1, pp. 31–39 35
На третьем участке медленной кривой при Ω ∈ [1:8; 4]: пусть Ω = 3;R = 0:77533156489, тогда
@R(Ω)
@Ω (3; 0:7753315648) = 0:3575757572 > 0, следовательно, участок является неустойчивым.
На четвертом участке медленной кривой при Ω ∈ [4; 8:64]: пусть Ω = 6;R = 1:514849187, тогда
@R(Ω)
@Ω (6; 1:514849187) = −0:4403735633 < 0, следовательно, участок является устойчивым.
В результате устойчивыми участками медленной кривой являются промежутки:
Ω ∈ [0; 0:86]; [0:86; 1:8]; [4; 8:64];
а неустойчивым участком:
Ω ∈ [1:8; 4]:
Точками срыва являются точки экстремумов медленной кривой, следовательно, исходя из рис. 2.3
имеются две точки срыва с координатами: A(1:8; 5:194117653) и B(4; 0:5666666669). Точки A и B явля-
ются точками смены устойчивости (рис. 2.3).
2.3. Релаксационные колебания
Пусть F(Ω; I;R) = "(  − Ω). Тогда система уравнений (2.5) будет выглядеть следующим образом:
J ˙Ω = − k2Ω
R+r +M(Ω);
˙R
= "(  − Ω):
(2.9)
Зададим точку (Ω∗;R∗) как положение равновесия системы уравнений (2.9):
− k2Ω
R+r +M(Ω∗) = 0;
"(  − Ω∗) = 0:
(2.10)
Выполнив преобразования, получим координаты положения равновесия:
Ω∗ =  ;
R∗ = k2Ω
M(Ω)
− r:
(2.11)
Исследуем на устойчивость положения равновесия при различных значениях параметра  :
Матрица Якоби системы (2.9) в общем виде будет выглядеть следующим образом:
A =
0
BBB@
@f
@Ω
@f
@R
@g
@Ω
@g
@R
1
CCCA
: (2.12)
Рассмотрим 3 случая, когда положение равновесия находится на:
1) на первом устойчивом участке медленной кривой при Ω ∈ [0; 1:8];
2) на втором устойчивом участке медленной кривой при Ω ∈ [4; 8:64];
3) на неустойчивом участке медленной кривой при Ω ∈ [1:8; 4].
Рассмотрим первый случай. Так как момент аэродинамических сил задается кусочно-линейным об-
разом и Ω ∈ [0; 1:8], то возьмем уравнение момента аэродинамических сил для второго участка.
Матрица Якоби в точке ( ; k2 
 
− r), при Ω ∈ [0:86; 1:8]; k = 0:5; r = 0:1; J = 1:
A|
( ; k2 
 
−r) = =
0
B@
−0:199893617 − 0:0638297872 
 
− 0:0638297872
(0:199893617 − 0:0638297872 )2
k2 
−1 0
1
CA
:
(2.13)
Характеристическое уравнение будет равно:
2 + ( 0:19989361700
  ) + 0:016296966936  − 0:10207333629 + 0:15982983247
  = 0: (2.14)
Найдем дискриминант и корни уравнения (2.14) при   = 1:
1 = −0:099946808500 − 0:25310886707i;
2 = −0:099946808500 + 0:25310886707i:
(2.15)
При полученных собственных значениях положение равновесия является асимптотически устойчивым
фокусом (рис. 2.4).
Рассмотрим второй случай. Так как момент аэродинамических сил задается кусочно-линейным об-
разом и Ω ∈ [4; 8:64], то возьмем уравнение момента аэродинамических сил для четвертого участка.
36
Кирсанова А.С. Релаксационные колебания в модели ветроэнергетической установки Дарье
Kirsanova A.S Relaxation oscillations in the Darie wind power plant model
Рис. 2.4. Медленная кривая (тонкая линия) и траектория системы (жирная линия)
при " = 0:01; k = 0:5; r = 0:1; J = 1;   = 1; начальная точка: Ω(0) = 7:5;R(0) = 2:3
Fig. 2.4. Slow curve (thin line) and system trajectory (thick line) at " = 0:01; k = 0:5; r = 0:1; J = 1;   = 1;
starting point: Ω(0) = 7:5;R(0) = 2:3
Действуя аналогичным образом, как и в предыдущем случае, находим корни характеристического урав-
нения при   = 4:5:
1 = −0:29358237550 − 1:2454650976i;
2 = −0:29358237550 + 1:2454650976i:
(2.16)
При полученных собственных значениях положение равновесия является асимптотически устойчивым
фокусом (рис. 2.5).
Рис. 2.5. Медленная кривая (тонкая линия) и траектория системы (жирная линия)
при " = 0:01; k = 0:5; r = 0:1; J = 1;   = 4:5; начальная точка: Ω(0) = 7:5;R(0) = 2:3
Fig. 2.5. Slow curve (thin line) and system trajectory (thick line) at " = 0:01; k = 0:5; r = 0:1; J = 1;   = 4:5;
starting point: Ω(0) = 7:5;R(0) = 2:3
Наконец, рассмотрим третий случай. Так как момент аэродинамических сил задается кусочно-линей-
ным образом и Ω ∈ [4; 8:64], то возьмем уравнение момента аэродинамических сил для третьего участка.
Действуя аналогичным образом, находим корни характеристического уравнения при   = 3:
1 = 0:17878787875 − 0:97308003575i;
2 = 0:17878787875 + 0:97308003575i:
(2.17)
При полученных собственных значениях положение равновесия находится на неустойчивом участке
медленной кривой и является неустойчивым фокусом. В этом случае в системе наблюдается релаксаци-
онный цикл (рис. 2.6). При этом траектория сначала движется вдоль устойчивого участка медленной
кривой, далее особая точка проходит точку экстремума (точку срыва) и при этом теряет устойчивость,
что соответствует бифуркации Андронова — Хопфа, т. е. при прохождении параметра   через бифур-
кационное значение устойчивое положение равновесия дифференциальной системы становится неустой-
чивым. Далее срывается и переходит в "-окрестность другого устойчивого участка, следует вдоль него
до второй точки срыва, после которого возвращается на первый устойчивый участок (рис. 2.6) [12].
Вестник Самарского университета. Естественнонаучная серия 2024. Том 30, № 1. С. 31–39
Vestnik of Samara University. Natural Science Series 2024, vol. 30, no. 1, pp. 31–39 37
Рис. 2.6. Медленная кривая (тонкая линия) и предельный цикл (жирная линия)
при " = 0:01; k = 0:5; r = 0:1; J = 1;   = 3; начальная точка: Ω(0) = 7:5;R(0) = 2:3
Fig. 2.6. Slow curve (thin line) and limit cycle (thick line) at " = 0:01; k = 0:5; r = 0:1; J = 1;   = 3; starting point:
Ω(0) = 7:5;R(0) = 2:3
Суммируем вышесказанное. В зависимости от изменении значения параметра особая точка меняет свое положение на медленной кривой. Если положение равновесия лежит на устойчивой части медленной кривой, то оно является асимптотически устойчивым. С физической точки зрения такая ситуация представляется предпочтительной, поскольку соответствует стационарному режиму работы системы. Случай, когда положение равновесия находится на неустойчивом участке и возникает релаксационный цикл, представляется неприемлемым с физической точки зрения, так как в системе возникают колебания с достаточно большой амплитудой.
Выводы
Рассмотрена трехмерная математическая модель ветрогенератора Дарье. При кусочно-линейной аппроксимации аэродинамического момента, найдены условия возникновения режимов работы ветрогенератора, соответствующие релаксационным колебаниям.

About the authors

A. S. Kirsanova

Author for correspondence.
Email: askirsanova99@gmail.com
ORCID iD: 0009-0005-1799-1423

postgraduate student of the Department of Differential Equations and Control Theory

References

Supplementary files