ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ ВЕТРОЭНЕРГЕТИЧЕСКОЙ УСТАНОВКОЙ КОЛЕБАТЕЛЬНОГО ТИПА
ПРИБОРЫ И СИСТЕМЫ АВТОМАТИЧЕСКОГО УПРАВЛЕНИЯ
УДК 62-50
С. А. КАБАНОВ, Е. С. ПЕТУХОВА
ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ ВЕТРОЭНЕРГЕТИЧЕСКОЙ УСТАНОВКОЙ КОЛЕБАТЕЛЬНОГО ТИПА
Рассматривается схема ветроэнергетической установки с колебательным рабочим движением. Представлена математическая модель динамики поворотных аэродинамических поверхностей (планов). Для обеспечения минимума выбранного критерия Красовского разработан алгоритм оптимального управления поворотом рабочих поверхностей с использованием модели прогнозирования движения.
Ключевые слова: ветроэнергетическая установка, оптимальное управление, модель прогнозирования движения.
В настоящее время такие источники энергии, как солнечная, ветровая, океаническая занимают в мировом производстве энергии лишь малую долю (менее 1 %), хотя являются возобновляемыми и относятся к экологически чистым. Ухудшение экологической обстановки обусловливает необходимость разработки устройств, позволяющих увеличить процент использования этих источников энергии.
До недавнего времени исследования в области ветроэнергетики сводились к изучению ветроэнергетических установок (ВЭУ) с вращательным рабочим режимом. Импеллерами (преобразователями энергии ветрового потока) в таких установках служат ветротурбины с горизонтальными или вертикальными осями [1]. ВЭУ с вращательным рабочим режимом имеют достаточно ограниченное применение. Это обусловлено большой массой, приходящейся на единицу полезной мощности; невысоким коэффициентом использования мощности ветрового потока, проходящего через рабочую площадь турбины; ограниченными рабочим диапазоном скоростей ветра и способностью адаптации к изменению ветрового режима; значительной стоимостью и др. [2]. В целях устранения указанных недостатков А. А. Красовским был предложен новый класс ВЭУ — ветроэнергетические установки с колебательным рабочим режимом (ВЭУК) [1, 2]. В них энергия ветрового потока преобразуется в механическую колебательную энергию посредством импеллеров, совершающих автоколебания. Возбуждение автоколебаний производится как за счет аэродинамических, упругих и инерционных сил, так и с помощью автоматического управления поворотом аэродинамических поверхностей — планов импеллера [1].
В настоящей статье рассматривается задача оптимального управления наземной ветроэнергетической установкой с автоколебательным рабочим режимом.
Установка, схема которой представлена на рис. 1, а [3, 4], имеет внешнюю раму 1, не участвующую в колебательном вертикальном движении, и внутреннюю раму 2, подвешенную
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3
Оптимальное управление ветроэнергетической установкой колебательного типа
33
во внешней на пружинах 3. Внутренняя рама (секция импеллера) содержит поворотные аэродинамические поверхности 4 (приводы поворота поверхностей на рисунке не показаны). В рабочем режиме рама 2 совершает возвратно-поступательные колебания, а рабочие поверхности 4 — синхронные колебания вокруг своей фиксированной оси при смене направления поступательного движения секции. Преобразование механической энергии возвратнопоступательного движения в электрическую может осуществляться различными способами.
а) 1
3
4
б)
2 α
ϑv –wв αв
Y ϑ–α Yа
Fа
y.
ϑ–α Х Fтр Fу Ха Fт 4
Рис. 1
Математическая модель ВЭУК [1, 3], приведенная на рис. 1, б, может быть представле-
на системой уравнений
x1 = x2 ;
⎫
x2
=
Fа (x2 , x3 ) − rx2 m
− cпр x1
−
mg
,⎪⎪⎬ ⎪
x3 = u.
⎪⎭
(1)
Здесь x = ( x1, x2 , x3 )T — вектор состояния системы со следующими компонентами:
x1 = y — величина перемещения секции импеллера вдоль вертикальной оси; x2 — скорость перемещения секции импеллера ( y на рис. 1, б); x3 = α — угол атаки рабочей поверхности;
Fa — сумма проекций аэродинамических сил ( Xa , Ya ) на вертикальную ось; rx2 = Fтр , где
Fтр — сила трения, r — коэффициент сопротивления; спр x1 = Fу , где Fу — сила упругости
пружин, спр — суммарная жесткость; mg = Fт , где Fт — сила тяжести, т — масса конструк-
ции, g — ускорение свободного падения.
Величина Fa определяется следующим образом:
Fa = 0,5ρSпΣWv ⎣⎡CY (x3 ) cos αв − CX (x3 ) (sin αв + x2 W )⎦⎤ ,
где ρ — плотность воздуха; SпΣ — суммарная площадь рабочих поверхностей секции; W —
модуль вектора скорости ветра wв ; v = x22 + W 2 − 2x2W sin αв — скорость движения рабочей поверхности относительно невозмущенного потока; αв — угол, образуемый вектором скорости ветра wв и нормалью к плоскости секции; CY (x3 ) , CX (x3 ) — аэродинамические характеристики, для аппроксимации которых используются следующие аналитические выражения:
CY (x3 ) = AY sin 2x3 ; СX (x3 ) = AX (1 − cos x3 ) + CX 0 ,
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3
34 С. А. Кабанов, Е. С. Петухова
здесь AY = 3 , AX = 0,5AY — постоянные коэффициенты; CX 0 = 0,1. Проекции (см. рис. 1, б)
Ya = 0, 5ρv2CY ( x3 ) Sп , Xa = 0, 5ρv2CX ( x3 ) Sп
— подъемная сила и сила лобового сопротивления рабочей поверхности соответственно. Начальные условия считаются заданными:
t0 = 0 , x1(0) = 0 , x2 (0) = 0 , x3 (0) = 0 .
Сигналом управления является скорость изменения угла α: x3 = α = u .
Угол поворота ϑ рабочей поверхности (см. рис. 1, б) связан с углом α следующим со-
отношением:
ϑ
=
α
+
arctg
⎛ ⎜
x2
⎝
W + sin αв cos αв
⎞ ⎟
.
⎠
Ограничений на управление не накладывается.
Требуется найти оптимальное управление uопт , обеспечивающее минимум функционала
Красовского [5] вида
∫I
=
Iт
+
t + ∆t
0, 5
t
u2
+ uо2пт k2
dτ,
Iт = ±[СY (x3 ) cos αв − CX (x3 )(sin αв + x2 W )].
(2)
Здесь Iт — терминальная часть функционала, отражающая полезную составляющую аэроди-
намической силы; t — момент начала оптимизации движения (принимает различные значения); ∆t — интервал оптимизации; k — постоянный коэффициент (его значение выбирается по результатам моделирования управляемого движения исходя из требований к переходному процессу).
В выражении (2) знак минус соответствует повороту рабочей поверхности для получения положительной аэродинамической силы, знак плюс — повороту рабочей поверхности для получения отрицательной аэродинамической силы. Смена знака происходит при изменении на-
правления поступательного движения секции (при смене знака переменной состояния x2 ).
Для минимизации функционала Красовского используем следующую модель прогнозирования движения при u(t) = 0 [5, 6] :
x1 = x2 ;
⎫
x2
=
Fа
(x2 , x3 )
−
rx2 m
−
cпр x1
−
mg
,⎪⎪⎬ ⎪
x3 = 0.
⎭⎪
Гамильтониан системы имеет вид
H
=
Hм
+
p3u + 0,5 u2
+ uо2пт k2
.
Обозначим
f
( x1 ,
x2 ,
x3 )
=
Fа
(x2 , x3 )
−
rx2 m
− cпр x1
−
mg
,
тогда
Hм = p1x2 + p2 f (x1, x2 , x3 ) .
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3
Оптимальное управление ветроэнергетической установкой колебательного типа
35
Здесь p = ( p1, p2 , p 3 )T — вектор сопряженных переменных (множителей Лагранжа).
Оптимальное
управление определяется
из
условия
∂H ∂u
=0
как
uопт
= −k 2 p3 .
Система
уравнений для сопряженных переменных имеет следующий вид:
при граничных условиях
p1
=
−
∂Hм ∂x1
=
−
p2
∂f
(x1, x2 ∂x1
,
x3
)
;
⎫ ⎪ ⎪
p 2
=
−
∂Hм ∂x2
= − p1 −
p2
∂f
(
x1, x2 ∂x2
,
x3
)
,⎬⎪ ⎪
p3
=
−
∂Hм ∂x3
= − p2
∂f (x1, x2 , x3 ) ∂x3
⎪ ⎪ ⎭
pT (t + ∆t) =
∂I т ∂x
t + ∆t .
Отсюда для отрицательных аэродинамических сил получим
p1(t + ∆t)
=
∂I т ∂x1
t
+ ∆t
=
0;
⎫ ⎪ ⎪
p2 (t + ∆t) =
∂I т ∂x2
t + ∆t
=
−СX (x3 ) ; W
⎪⎪ ⎬ ⎪
p3 (t
+
∆t)
=
∂I т ∂x3
t
+
∆t
=
⎡ ⎢ ⎣
∂CY (x3 ∂x3
)
cos
αв
−
∂CX (x3 ∂x3
)
⎛ ⎝⎜
sin
αв
+
x2 W
⎞⎤ ⎠⎟⎥⎦
⎪ ,⎪ ⎪⎭
для положительных аэродинамических сил —
p1 (t
+ ∆t)
=
∂I т ∂x1
t
+ ∆t
=
0;
⎫ ⎪ ⎪
p2 (t + ∆t)
=
∂I т ∂x2
t + ∆t
=
СX (x3 ) , W
⎪⎪ ⎬ ⎪
p3 (t
+ ∆t)
=
∂I т ∂x3
t
+ ∆t
=
⎡⎢− ⎣
∂CY (x3 ∂x3
)
cos
αв
+
∂CX (x3 ∂x3
)
⎛ ⎝⎜
sin
αв
+
x2 W
⎞⎤ ⎠⎟⎦⎥
⎪ .⎪ ⎪⎭
Поскольку задача является нетерминальной (момент окончания процесса заранее не назначается), используется скользящий интервал оптимизации, величина которого определяется по результатам моделирования.
Математическое моделирование управляемого движения ВЭУК производилось при следующих значениях параметров системы и внешней среды: m = 5 кг, α в=0, спр=100 Н/м, W = 5 м/с, r = 38 Н·с/м, SпΣ = 1 м2, ρ = 1, 29 кг/м3, k = 3,5 , ∆t = 0,1 с, шаг интегрирова-
ния 0, 001 с. Результаты моделирования в виде зависимостей величины перемещения x1
секции импеллера, угла поворота ϑ этой секции и угла атаки α от времени представлены на рис. 2—4.
Итак, рассмотрено решение нетерминальной задачи оптимального управления ветроэнергетической установкой с автоколебательным рабочим режимом на основе алгоритма,
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3
36 С. А. Кабанов, Е. С. Петухова
использующего модель прогнозирования движения. Как показали результаты моделирования, применение этого алгоритма позволяет решить задачу оптимального управления, обеспечивая при этом требуемое качество переходных процессов.
x1, м ϑ, рад 01
–0,2
–0,4 0
–0,6
–0,8
–1 –1
1 2 3 4 5 6 7 t, с
1 2 3 4 5 6 7 t, с
Рис. 2 α, рад
Рис. 3
0
1 2 3 4 5 6 7 t, с Рис. 4
При решении задачи были определены возможные направления дальнейших исследований. Моделирование в среде „Fluent“ процесса обтекания профиля поверхности ВЭУК невязким газом показывает, что предельное значение скорости изменения угла атаки ограничено. Поэтому при синтезе установки следует ввести ограничение на управление. Также в системе (1) следует учесть инерционность динамики рабочей поверхности при ее повороте на оптимальный угол. Кроме того, при разработке конструкции ВЭУК должны быть предусмотрены возможности экстренного торможения, снабжения системой антиобледенения и системой контроля уровня вибраций. Решение этих вопросов позволит создать экономичные и экологически чистые системы энергообеспечения.
Исследования, описанные в настоящей статье, выполнены по гранту Российского фонда фундаментальных исследований, № 09-08-00829.
СПИСОК ЛИТЕРАТУРЫ
1. Красовский А. А. Синтез автоколебательных систем с приложением к ветроэнергетической установке нового класса // Изв. РАН. Сер. Техн. кибернетика. 1994. № 6. С. 5—15.
2. Красовский А. А. Синтез цилиндроэллиптических аттракторов с приложением к ветроэнергетической установке нового класса // Изв. РАН. Сер. Теория и системы управления. 1996. № 2. С. 21—30.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3
Оптимальное управление ветроэнергетической установкой колебательного типа
37
3. Красовский А. А. Модульные ветроэнергетические установки с управляемым колебательным рабочим движением — путь решения энергетических проблем // Там же. 2001. № 6. С. 145—151.
4. Красовский А. А. Синтез полиномиальных автоколебательных систем пятой—седьмой степени // Автоматика и телемеханика. 2002. № 10. С. 44—57.
5. Справочник по теории автоматического управления / Под ред. А. А. Красовского. М.: Наука, 1987.
6. Кабанов С. А. Управление системами на прогнозирующих моделях. СПб.: Изд-во СПбГУ, 1997.
Сергей Александрович Кабанов Елена Сергеевна Петухова
Сведения об авторах — д-р техн. наук, профессор; Балтийский государственный технический
университет „Военмех“ им. Д. Ф. Устинова, кафедра систем обработки информации и управления, Санкт-Петербург; E-mail: kabanov@bstu.spb.su — инженер; Балтийский государственный технический университет „Военмех“ им. Д. Ф. Устинова, кафедра систем обработки информации и управления, Санкт-Петербург; E-mail: helpet@yandex.ru
Рекомендована кафедрой систем обработки информации и управления
Поступила в редакцию 07.12.07 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3
УДК 62-50
С. А. КАБАНОВ, Е. С. ПЕТУХОВА
ОПТИМАЛЬНОЕ УПРАВЛЕНИЕ ВЕТРОЭНЕРГЕТИЧЕСКОЙ УСТАНОВКОЙ КОЛЕБАТЕЛЬНОГО ТИПА
Рассматривается схема ветроэнергетической установки с колебательным рабочим движением. Представлена математическая модель динамики поворотных аэродинамических поверхностей (планов). Для обеспечения минимума выбранного критерия Красовского разработан алгоритм оптимального управления поворотом рабочих поверхностей с использованием модели прогнозирования движения.
Ключевые слова: ветроэнергетическая установка, оптимальное управление, модель прогнозирования движения.
В настоящее время такие источники энергии, как солнечная, ветровая, океаническая занимают в мировом производстве энергии лишь малую долю (менее 1 %), хотя являются возобновляемыми и относятся к экологически чистым. Ухудшение экологической обстановки обусловливает необходимость разработки устройств, позволяющих увеличить процент использования этих источников энергии.
До недавнего времени исследования в области ветроэнергетики сводились к изучению ветроэнергетических установок (ВЭУ) с вращательным рабочим режимом. Импеллерами (преобразователями энергии ветрового потока) в таких установках служат ветротурбины с горизонтальными или вертикальными осями [1]. ВЭУ с вращательным рабочим режимом имеют достаточно ограниченное применение. Это обусловлено большой массой, приходящейся на единицу полезной мощности; невысоким коэффициентом использования мощности ветрового потока, проходящего через рабочую площадь турбины; ограниченными рабочим диапазоном скоростей ветра и способностью адаптации к изменению ветрового режима; значительной стоимостью и др. [2]. В целях устранения указанных недостатков А. А. Красовским был предложен новый класс ВЭУ — ветроэнергетические установки с колебательным рабочим режимом (ВЭУК) [1, 2]. В них энергия ветрового потока преобразуется в механическую колебательную энергию посредством импеллеров, совершающих автоколебания. Возбуждение автоколебаний производится как за счет аэродинамических, упругих и инерционных сил, так и с помощью автоматического управления поворотом аэродинамических поверхностей — планов импеллера [1].
В настоящей статье рассматривается задача оптимального управления наземной ветроэнергетической установкой с автоколебательным рабочим режимом.
Установка, схема которой представлена на рис. 1, а [3, 4], имеет внешнюю раму 1, не участвующую в колебательном вертикальном движении, и внутреннюю раму 2, подвешенную
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3
Оптимальное управление ветроэнергетической установкой колебательного типа
33
во внешней на пружинах 3. Внутренняя рама (секция импеллера) содержит поворотные аэродинамические поверхности 4 (приводы поворота поверхностей на рисунке не показаны). В рабочем режиме рама 2 совершает возвратно-поступательные колебания, а рабочие поверхности 4 — синхронные колебания вокруг своей фиксированной оси при смене направления поступательного движения секции. Преобразование механической энергии возвратнопоступательного движения в электрическую может осуществляться различными способами.
а) 1
3
4
б)
2 α
ϑv –wв αв
Y ϑ–α Yа
Fа
y.
ϑ–α Х Fтр Fу Ха Fт 4
Рис. 1
Математическая модель ВЭУК [1, 3], приведенная на рис. 1, б, может быть представле-
на системой уравнений
x1 = x2 ;
⎫
x2
=
Fа (x2 , x3 ) − rx2 m
− cпр x1
−
mg
,⎪⎪⎬ ⎪
x3 = u.
⎪⎭
(1)
Здесь x = ( x1, x2 , x3 )T — вектор состояния системы со следующими компонентами:
x1 = y — величина перемещения секции импеллера вдоль вертикальной оси; x2 — скорость перемещения секции импеллера ( y на рис. 1, б); x3 = α — угол атаки рабочей поверхности;
Fa — сумма проекций аэродинамических сил ( Xa , Ya ) на вертикальную ось; rx2 = Fтр , где
Fтр — сила трения, r — коэффициент сопротивления; спр x1 = Fу , где Fу — сила упругости
пружин, спр — суммарная жесткость; mg = Fт , где Fт — сила тяжести, т — масса конструк-
ции, g — ускорение свободного падения.
Величина Fa определяется следующим образом:
Fa = 0,5ρSпΣWv ⎣⎡CY (x3 ) cos αв − CX (x3 ) (sin αв + x2 W )⎦⎤ ,
где ρ — плотность воздуха; SпΣ — суммарная площадь рабочих поверхностей секции; W —
модуль вектора скорости ветра wв ; v = x22 + W 2 − 2x2W sin αв — скорость движения рабочей поверхности относительно невозмущенного потока; αв — угол, образуемый вектором скорости ветра wв и нормалью к плоскости секции; CY (x3 ) , CX (x3 ) — аэродинамические характеристики, для аппроксимации которых используются следующие аналитические выражения:
CY (x3 ) = AY sin 2x3 ; СX (x3 ) = AX (1 − cos x3 ) + CX 0 ,
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3
34 С. А. Кабанов, Е. С. Петухова
здесь AY = 3 , AX = 0,5AY — постоянные коэффициенты; CX 0 = 0,1. Проекции (см. рис. 1, б)
Ya = 0, 5ρv2CY ( x3 ) Sп , Xa = 0, 5ρv2CX ( x3 ) Sп
— подъемная сила и сила лобового сопротивления рабочей поверхности соответственно. Начальные условия считаются заданными:
t0 = 0 , x1(0) = 0 , x2 (0) = 0 , x3 (0) = 0 .
Сигналом управления является скорость изменения угла α: x3 = α = u .
Угол поворота ϑ рабочей поверхности (см. рис. 1, б) связан с углом α следующим со-
отношением:
ϑ
=
α
+
arctg
⎛ ⎜
x2
⎝
W + sin αв cos αв
⎞ ⎟
.
⎠
Ограничений на управление не накладывается.
Требуется найти оптимальное управление uопт , обеспечивающее минимум функционала
Красовского [5] вида
∫I
=
Iт
+
t + ∆t
0, 5
t
u2
+ uо2пт k2
dτ,
Iт = ±[СY (x3 ) cos αв − CX (x3 )(sin αв + x2 W )].
(2)
Здесь Iт — терминальная часть функционала, отражающая полезную составляющую аэроди-
намической силы; t — момент начала оптимизации движения (принимает различные значения); ∆t — интервал оптимизации; k — постоянный коэффициент (его значение выбирается по результатам моделирования управляемого движения исходя из требований к переходному процессу).
В выражении (2) знак минус соответствует повороту рабочей поверхности для получения положительной аэродинамической силы, знак плюс — повороту рабочей поверхности для получения отрицательной аэродинамической силы. Смена знака происходит при изменении на-
правления поступательного движения секции (при смене знака переменной состояния x2 ).
Для минимизации функционала Красовского используем следующую модель прогнозирования движения при u(t) = 0 [5, 6] :
x1 = x2 ;
⎫
x2
=
Fа
(x2 , x3 )
−
rx2 m
−
cпр x1
−
mg
,⎪⎪⎬ ⎪
x3 = 0.
⎭⎪
Гамильтониан системы имеет вид
H
=
Hм
+
p3u + 0,5 u2
+ uо2пт k2
.
Обозначим
f
( x1 ,
x2 ,
x3 )
=
Fа
(x2 , x3 )
−
rx2 m
− cпр x1
−
mg
,
тогда
Hм = p1x2 + p2 f (x1, x2 , x3 ) .
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3
Оптимальное управление ветроэнергетической установкой колебательного типа
35
Здесь p = ( p1, p2 , p 3 )T — вектор сопряженных переменных (множителей Лагранжа).
Оптимальное
управление определяется
из
условия
∂H ∂u
=0
как
uопт
= −k 2 p3 .
Система
уравнений для сопряженных переменных имеет следующий вид:
при граничных условиях
p1
=
−
∂Hм ∂x1
=
−
p2
∂f
(x1, x2 ∂x1
,
x3
)
;
⎫ ⎪ ⎪
p 2
=
−
∂Hм ∂x2
= − p1 −
p2
∂f
(
x1, x2 ∂x2
,
x3
)
,⎬⎪ ⎪
p3
=
−
∂Hм ∂x3
= − p2
∂f (x1, x2 , x3 ) ∂x3
⎪ ⎪ ⎭
pT (t + ∆t) =
∂I т ∂x
t + ∆t .
Отсюда для отрицательных аэродинамических сил получим
p1(t + ∆t)
=
∂I т ∂x1
t
+ ∆t
=
0;
⎫ ⎪ ⎪
p2 (t + ∆t) =
∂I т ∂x2
t + ∆t
=
−СX (x3 ) ; W
⎪⎪ ⎬ ⎪
p3 (t
+
∆t)
=
∂I т ∂x3
t
+
∆t
=
⎡ ⎢ ⎣
∂CY (x3 ∂x3
)
cos
αв
−
∂CX (x3 ∂x3
)
⎛ ⎝⎜
sin
αв
+
x2 W
⎞⎤ ⎠⎟⎥⎦
⎪ ,⎪ ⎪⎭
для положительных аэродинамических сил —
p1 (t
+ ∆t)
=
∂I т ∂x1
t
+ ∆t
=
0;
⎫ ⎪ ⎪
p2 (t + ∆t)
=
∂I т ∂x2
t + ∆t
=
СX (x3 ) , W
⎪⎪ ⎬ ⎪
p3 (t
+ ∆t)
=
∂I т ∂x3
t
+ ∆t
=
⎡⎢− ⎣
∂CY (x3 ∂x3
)
cos
αв
+
∂CX (x3 ∂x3
)
⎛ ⎝⎜
sin
αв
+
x2 W
⎞⎤ ⎠⎟⎦⎥
⎪ .⎪ ⎪⎭
Поскольку задача является нетерминальной (момент окончания процесса заранее не назначается), используется скользящий интервал оптимизации, величина которого определяется по результатам моделирования.
Математическое моделирование управляемого движения ВЭУК производилось при следующих значениях параметров системы и внешней среды: m = 5 кг, α в=0, спр=100 Н/м, W = 5 м/с, r = 38 Н·с/м, SпΣ = 1 м2, ρ = 1, 29 кг/м3, k = 3,5 , ∆t = 0,1 с, шаг интегрирова-
ния 0, 001 с. Результаты моделирования в виде зависимостей величины перемещения x1
секции импеллера, угла поворота ϑ этой секции и угла атаки α от времени представлены на рис. 2—4.
Итак, рассмотрено решение нетерминальной задачи оптимального управления ветроэнергетической установкой с автоколебательным рабочим режимом на основе алгоритма,
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3
36 С. А. Кабанов, Е. С. Петухова
использующего модель прогнозирования движения. Как показали результаты моделирования, применение этого алгоритма позволяет решить задачу оптимального управления, обеспечивая при этом требуемое качество переходных процессов.
x1, м ϑ, рад 01
–0,2
–0,4 0
–0,6
–0,8
–1 –1
1 2 3 4 5 6 7 t, с
1 2 3 4 5 6 7 t, с
Рис. 2 α, рад
Рис. 3
0
1 2 3 4 5 6 7 t, с Рис. 4
При решении задачи были определены возможные направления дальнейших исследований. Моделирование в среде „Fluent“ процесса обтекания профиля поверхности ВЭУК невязким газом показывает, что предельное значение скорости изменения угла атаки ограничено. Поэтому при синтезе установки следует ввести ограничение на управление. Также в системе (1) следует учесть инерционность динамики рабочей поверхности при ее повороте на оптимальный угол. Кроме того, при разработке конструкции ВЭУК должны быть предусмотрены возможности экстренного торможения, снабжения системой антиобледенения и системой контроля уровня вибраций. Решение этих вопросов позволит создать экономичные и экологически чистые системы энергообеспечения.
Исследования, описанные в настоящей статье, выполнены по гранту Российского фонда фундаментальных исследований, № 09-08-00829.
СПИСОК ЛИТЕРАТУРЫ
1. Красовский А. А. Синтез автоколебательных систем с приложением к ветроэнергетической установке нового класса // Изв. РАН. Сер. Техн. кибернетика. 1994. № 6. С. 5—15.
2. Красовский А. А. Синтез цилиндроэллиптических аттракторов с приложением к ветроэнергетической установке нового класса // Изв. РАН. Сер. Теория и системы управления. 1996. № 2. С. 21—30.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3
Оптимальное управление ветроэнергетической установкой колебательного типа
37
3. Красовский А. А. Модульные ветроэнергетические установки с управляемым колебательным рабочим движением — путь решения энергетических проблем // Там же. 2001. № 6. С. 145—151.
4. Красовский А. А. Синтез полиномиальных автоколебательных систем пятой—седьмой степени // Автоматика и телемеханика. 2002. № 10. С. 44—57.
5. Справочник по теории автоматического управления / Под ред. А. А. Красовского. М.: Наука, 1987.
6. Кабанов С. А. Управление системами на прогнозирующих моделях. СПб.: Изд-во СПбГУ, 1997.
Сергей Александрович Кабанов Елена Сергеевна Петухова
Сведения об авторах — д-р техн. наук, профессор; Балтийский государственный технический
университет „Военмех“ им. Д. Ф. Устинова, кафедра систем обработки информации и управления, Санкт-Петербург; E-mail: kabanov@bstu.spb.su — инженер; Балтийский государственный технический университет „Военмех“ им. Д. Ф. Устинова, кафедра систем обработки информации и управления, Санкт-Петербург; E-mail: helpet@yandex.ru
Рекомендована кафедрой систем обработки информации и управления
Поступила в редакцию 07.12.07 г.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 3