CC BY
166
УДК 620.91
Математическая модель динамики вращения ветротурбины и поворота малой ветроэнергетической установки для фермерских хозяйств посредством автоматической системы ориентации
А. С. Гуринов, В. Л. Гапонов
(Донской государственный технический университет)
Представлено математическое моделирование динамики вращения ветротурбины и поворота ветроэнергетической установки. Показано распределение углов атаки лопасти по азимуту вращения ветротурбины при угле косой обдувки 30°. Видно, что обтекание при наличии ветра несоосного с осью вращения приводит к существенной неравномерности углов атаки и, соответственно, к изменению сил на лопастях в процессе вращения.
Ключевые слова: ветроэнергетическая установка, моделирование динамики поворота, моделирование динамики вращения вегрогенератора, лопасть.
Введение. Эффективность горизонтально-осевой ВЭУ зависит от многих факторов. Одним из важных показателей является способность установки быстро перестраиваться навстречу ветру с наименьшим углом к нему. В то же время, вследствие инерционных сил, которые затрудняют перестройку ветроустановки вслед за изменениями метеорологической обстановки, ВЭУ постоянно находится в состоянии присутствия косой обдувки, т. е. наличия угла между осью вращения и направлением ветра. Свою лепту в наращивание величины вносит и трение в опоре горизонтальноосевой ВЭУ. Для того чтобы оценить влияние различных факторов и определить эффективность ВЭУ уже на стадии разработки, целесообразно иметь математическую модель ветроустановки. Для её создания необходимо рассмотреть основные силы, действующие на ветротурбину. Моделирование динамики поворота ВЭУ. Как известно, вращение ветротурбины обусловлено наличием аэродинамических сил, создающих крутящий момент на валу генератора. Однако этому вращению препятствует момент генератора, появляющийся вследствие наличия сил трения в мультипликаторе, опорах вала и сил сопротивления вследствие электромагнитных сил.
Суммарное уравнение моментов на валу ветрогенератора можно представить следующим образом:
мк-мг-му^ =(1вт+шмп(1)
где Мк — крутящий момент, обусловленный наличием аэродинамических сил; Мг — момент сопротивления на генераторе, обусловленный наличием электрической нагрузки; Ми — мо-
Утр
мент, обусловленный силами трения в мультипликаторе и опорах вала; 1ВТ — момент инерции ветротурбины; п— передаточное отношение мультипликатора генератора; 1МП — момент инерции ротора генератора; соу —угловая скорость вращения ветротурбины.
При этом разворачивающий момент ВЭУ равен:
м -м -м -м =т с^>1 т
"пер "і "гир ^ВЭУ ^ ^
где М1 — момент от аэродинамических сил на ветротурбине, обусловленный косой обдувкой;
Мпер ~ момент от сил, действующих на оперение; Мгир — гироскопический момент ветротур-
1763
бины; М7 — момент от сил трения в опоре ветротурбины; I — момент инерции ветротур-
тр ВЭУ
бины относительно оси поворота; со2 — угловая скорость поворота ВЭУ относительно вертикальной оси.
Момент генератора является функцией электрической нагрузки Мг = /"(£), трение в опоре величина постоянная
где кп — коэффициент трения в подшипнике; твэу — масса поворачивающейся части ВЭУ; с!п — диаметр сепаратора опорного подшипника.
Моделирование динамики вращения горизонтально-осевой ветротурбины. Крутящий момент, обусловленный наличием аэродинамических сил, зависит от характера обтекания ветротурбины и, в первую очередь, наличием косой обдувки. Он может быть рассчитан несколькими методами. Наиболее оптимальным, с точки зрения сложности и точности, является метод расчёта единичного элемента лопасти [1, 2, 3]. Этот метод расчёта основывается на допущении, что каждый элемент лопасти ветротурбины можно рассматривать как отрезок профиля, двигающийся по винтовой линии. Аэродинамические силы рассчитываются, исходя из результирующей скорости потока, обтекающего этот профиль. Причём считается, что прилегающие части лопасти не влияют на характеристики рассматриваемого отрезка профиля. Силы и крутящий момент винта получаются интегрированием элементарных сил и моментов отдельных элементов лопасти. Соответственно задача определения аэродинамических характеристик ветротурбины решается при следующих допущениях:
— лопасть является абсолютно жёсткой на изгиб и кручение;
— аэродинамические характеристики профиля сечения лопасти не зависят от скольжения потока вдоль лопасти и влияния центробежных сил на пограничный слой.
На элементе лопасти при вращении возникают аэродинамические силы, обусловленные наличием как окружной, так и набегающего потока через ветротурбину. Представленное на рис. 1 соотношение скоростей сечения лопасти соответствует лопасти, находящейся в верхнем положении с обдувкой набегающего потока под углом ав к оси вращения. В этом случае вертикальная составляющая скорости потока иу отрицательна и угол атаки ар меньше угла установки.
Скорость продольного движения представляет собой сумму скорости набегающего потока и окружной скорости.
где г — текущее значение радиуса; со — угловая скорость несущего винта; 1/К — скорость ветра; Ф — азимут лопасти.
Вертикальная составляющая в плоскости сечения профиля скорости элемента лопасти будет равна:
(3)
их = СОЛ + IV БІП ф
(4)
иу = 1/К соэ аБ,
(5)
где ав — угол между осью вращения и направлением ветра. Угол притекания имеет вид:
3 = Ф + ар
и
3 = агсЬд^
(6)
(7)
ор = 3-Ф
(8)
1764
Рис. 2. Распределение аэродинамических сил на сечении лопасти
1765
Значение осевой силы элемента лопасти равно:
ОТ = дХпр эт 3 + бУпр соэ 3.
Значение элемента силы, направленной на вращение ветротурбины:
(9)
сі(2 = сіГпр5т$-сІХпрсо5$
сім к = я/<?
(Ю)
(П)
(12)
(13)
где су — коэффициент подъёмной силы профиля на участке элемента лопасти; сх — коэффициент силы сопротивления профиля на участке элемента лопасти; иг — суммарная скорость обтекания элемента лопасти; их, иу — проекции скорости движения элемента лопасти; Ь —
хорда лопасти; р — плотность воздуха; с1г — длина элемента лопасти вдоль радиуса.
Величина коэффициента подъёмной силы в линейной зоне может быть определена по
где а, — производная коэффициента подъёмной силы по углу атаки; ор — текущий угол атаки
сечения относительно точки нулевой подъёмной силы профиля; ф — угол установки профиля. Текущий угол установки профиля представляет собой:
Принимая во внимание тот факт, что помимо вращения ветротурбины, присутствует и второе вращательное движение, а именно вращение всей установки ветрогенератора, тогда скорость продольного движения, а также вертикальная составляющая скорости элемента лопасти принимают несколько иной вид.
Плечо скорости вращения V , действующей на произвольную точку лопасти, равно:
где гт — плечо вращающей силы; л — расстояние от центра лопасти до точки приложения сил; т— расстояние от оси вращения ветрогенератора до его ветротурбины.
Для проведения интегрирования величины аэродинамических сил ветротурбины по радиусу необходимо учесть потерю эффективности ветротурбины на краю диска вследствие перетекания воздуха на концах лопастей. Вследствие разницы давления воздух перетекает с нижней поверхности на верхнюю, уменьшается разница в давлении и уменьшаются и силы в области концов
формуле:
Г и \
Су =орая =а„ ф + -^ ,
V ху
(14)
(15)
(16)
4- 111
у = агйд — г
(17)
Кр =«А,С05У
их = СОЛ + 1/К БІП ав СОБф і!у = 1/К созав + ]/вр біп ф
(18)
(19)
(20)
1766
лопастей. Зона распространения концевых потерь зависит от удельной нагрузки на лопасть и на всю ветротурбину.
Однако мало нагруженные ветротурбины можно рассчитывать при принятии некоторых простых допущений, касающихся концевых потерь. Одним из них является допущение о том, что часть лопасти длиной, равной половине концевой хорды, не создаёт аэродинамических сил. В этом случае относительный радиус рабочей части лопасти равен:
В = 1- —.
2 Л
Осевое усилие на ветротурбине равно:
360 5/?
Т = г [ с/ф | с1Т .
О га
Интегрированное значение момента лопасти определяется так:
360 5/?
М2=г\(^\ гтс!Т,
(21)
(22)
(23)
где г — количество лопастей.
Таким образом, момент ветротурбины представляет собой:
((
Мх =г [ с/ф|
р |(сол + 1/К зтав соб ф)2 + (1/К соб ав + солт соб у эт ф)2)
Л
с„Ьс!г
р |(сол + 1/К зтав соб ф)2 + (1/К соб ав + согт соб у эт ф)2)
суЬс!г
собЗ
гЛг
(24)
Крутящий момент определяется следующим образом:
360 5/?
Мк =г [ с/ф| яй?.
0 га
При численном интегрировании лопасть разбивается на элементы длиной:
", '
где пь — количество отрезков, на которые разбивается лопасть.
Радиус /-го элемента лопасти равен:
г -я + —
1М1 ^0^2'
= ГМ(1-1) + ^Гм •
Выводы. На основании предложенной выше методики было выполнено моделирование экспериментальной ветротурбины диаметром 2 м на разных углах косой обдувки. На рис. 3 показано распределение углов атаки лопасти по азимуту вращения при угле косой обдувки 30°. Видно, что обтекание при наличии ветра, несоосного с осью вращения, приводит к существенной неравномерности углов атаки и, соответственно, распределения сил на лопастях в процессе вращения. Причём если угол между направлением вращения и направлением ветра присутствует в горизонтальной плоскости, то основной перекос аэродинамических усилий присутствует в вертикальной плоскости.
(25)
(26)
(27)
(28)
1767
Рис. 3. Распределение аэродинамических углов атаки профиля ветротурбины диаметром 2 м (ветер 6 м/с, угол косой обдувки 30°)
М
Л 1 І
470,00 д Д
370,00 Л д
А 270,00 Л
д 170,00 д
1 70,00 0 1
-60
40_______________-20 __________-20,00__________________20________________40
Рис. 4. Мощность на валу ВЭУ в зависимости от угла косой обдувки
1768
6С
—* 4-Ш—
ег мма аэр. моментов
Су
60 -40 -20 ) \ 2( 40 момент сил трения
У МО мент сил трения
/4,00
-22.00
Рис. 5. Значения моментов на ВЭУ от аэродинамических сил и сил трения
Расчёт крутящего момента при этом показал, что косая обдувка приводит к существенному снижению эффективности ВЭУ. Так, при угле 15° снижение вырабатываемой мощности составляет 11 %, а при 30° эта величина вырастает до 35 % (рис. 4).
Выполненные далее расчёты поворота всего ветрогенератора при постоянной частоте вращения позволили определить моменты аэродинамических сил на ветроколесе и стабилизаторе (рис. 5). Значение момента трения выглядит как две горизонтальные линии. Это фактически зона нечувствительности ВЭУ к повороту. Видно, что до значения угла косой обдувки 14° ВЭУ не будет совершать поворот относительно своей вертикальной оси.
Таким образом, можно констатировать, что ВЭУ всегда находится в зоне косой обдувки не только вследствие инерционных сил, но сил трения.
Библиографический список
1. Дудник, В. В. Определение оптимального азимута установки лопастей соосного сверхлёгкого вертолёта / В. В. Дудник, В. А. Колот // Вестник ДГТУ. — 2011. — Т. 11, № 5 (56). — С. 667-675.
2. Вертолёты. Расчёт и проектирование. В 3 книгах. Книга 1. Аэродинамика / М. Л. Миль [и др.]; под ред. М. Л. Миля. — Москва: Машиностроение, 1966. — 450 с.
3. Расчёт распределения рабочих и испытательных нагрузок на лопасть крупной ВЭУ /
Н. В. Белан [и др.]. — Харьков: ХАИ, 2006. — 49 с.
Материал поступил в редакцию 01.12.2011.
1769
References
1. Dudnik, V. V. Opredelenie optimal'nogo azimuta ustanovki lopastej soosnogo sverxlyogkogo vertolyota / V. V. Dudnik, V. A. Kolot // Vestnik DGTU. — 2011. — T. 11, # 5 (56). — S. 667—675. — In Russian.
2. Vertolyoty'. Raschyot i proektirovanie. V 3 knigax. Kniga 1. Ae'rodinamika / M. L. Mil' [i dr.]; pod red. M. L. Milya. — Moskva: Mashinostroenie, 1966. — 450 s. — In Russian.
3. Raschyot raspredeleniya rabochix i ispy'tatel'ny'x nagruzok na lopast' krupnoj VE'U / N. V. Belan [i dr.]. — Xar'kov: XAI, 2006. — 49 s. — In Russian.
MATHEMATICAL MODEL OF ROTATIONAL DYNAMICS OF WIND TURBINE AND SMALL AEROGENERATOR TURN FOR FARMS USING AUTOMATIC ORIENTATION SYSTEM
A. S. Gurinov, V. L. Gaponov
(Don State Technical University)
The mathematical simulation of wind turbine rotation and gondola turning is presented. The blade angle-of-attack distribution in azimuth of the wind turbine rotation with 30° angle of lateral wind is shown. The study shows that the nonaxial rotation flow causes critical irregularity of the angles of attack and appropriately force variation on the blades during the rotation.
Keywords: aerogenerator, turn dynamic simulation, wind turbine rotation dynamic simulation, blade.
1770
