CC BY
52
РАСЧЕТ ХАРАКТЕРИСТИК ПРОТИВОДЕЙСТВУЮЩЕЙ ПРУЖИНЫ ДЛЯ УСТРОЙСТВА ОГРАНИЧЕНИЯ МОЩНОСТИ ВЭУ
Петько Виктор Гаврилович
д-р техн. наук, профессор, Оренбургский ГАУ, г. Оренбург
Е-mail: VGPetko @mail. ru Пугачёв Владимир Валерьевич аспирант, Оренбургский ГАУ, г. Оренбург Е-mail: olorin777@mai.ru
CALCULATION OF CHARACTERISTICS AGAINST SPRING FOR UNIT LIMITS OF POWER WIND TURBINE
Victor Petko
Professor, Doctor of Technical Sciences, Orenburg State Agrarian University,
Orenburg Vladimir Pugachev
graduate student, Orenburg State Agrarian University, Orenburg АННОТАЦИЯ
В статье дано описание устройства ограничения мощности ветроэнергетической установки и приведены расчеты параметров пружины, противодействующей давлению ветра на боковую лопату.
ABSTRACT
The article describes the device limits the power of wind turbines and the calculations of the parameters of the spring, the wind pressure on the opposing side shovel.
Ключевые слова: ветроэнергетическая установка, ограничение мощности. Keywords: wind power plant, power limit.
В настоящее время в ветроэнергетике используются различные способы ограничения мощности ветроустановок при увеличении скорости ветра выше номинальной. Одними из основных являются: 1) изменение угла атаки лопастей ветротурбины; 2) использование определенного профиля лопастей, изменяющегося по всей длине лопасти; 3) вывод ветроколеса из-под ветра, т. е. изменение его положения относительно направления скорости ветра.
Использование такого устройства как виндроза 1 (небольшое ветроколесо, устанавливаемое за основной ветротурбиной) заключается в основном как один из способов ориентации на ветер. Нами же предлагается придать виндрозе возможность вращения относительно вертикальной 2 оси под действием потока ветра на связанную с ней боковую лопату 3 (рис. 1).
Для предотвращения отклонения боковой лопаты виндрозы при случайных порывах ветра ее необходимо закрепить с помощью пружины 4 или другого устройства к неподвижной части головки ветроагрегата 5. В качестве неподвижной части может выступать диск или дуга постоянного радиуса, к которому через трос закреплена пружина. При скоростях ветра не выше расчётной величины плоскость виндрозы повёрнута относительно плоскости ветроколеса 6 на угол 90о благодаря первоначальному натяжению пружины. Такая конструкция виндрозы будет не только ориентировать ветроустановку по ветру, но и ограничивать ее мощность при превышении скорости ветра выше номинальной.
При увеличении скорости ветра, он начнет действовать на боковую лопату виндрозы и повернет ее на определенный угол у. В результате крыльчатка виндрозы установится под углом к ветру и начнет вращаться, пока не
1
46
Рис. 1. Устройство ограничения мощности ВЭУ
установиться в прежнее состояние. Но при ее вращении начнет изменять свое положении основное ветроколесо, которое повернется на тот же угол у относительно направления ветра. Ограничение мощности в такой установке будет происходить за счет поворота ветроколеса на определенный угол, зависящий от скорости ветра.
Требуемый закон изменения у от скорости ветра, обеспечивающий неизменную мощность ветроколеса при скорости ветра V больше расчётной скорости ветра V>Vр может быть получен из выражения [1]:
^Гопт =7 (1.1)
Тогда
Гопт = (1.2)
Однако под действием боковой лопаты и виндрозы фактический угол поворота ветроколеса будет отличаться от оптимального его значения уфуопт. В результате и мощность ветроколеса при V>Vр не будет постоянной, равной расчётной мощности. С целью выявления, насколько эта разница будет существенной, следует определить, каков будет фактический угол поворота ветроколеса под действием боковой лопаты при той или иной скорости ветра V.
С одной стороны на виндрозу действует момент боковой лопаты, зависящий от её геометрических размеров и скорости:
МЛ = ЬРЛ = 0,5Ь раЬУ2 = КЛУ2 (1.3)
где ^ = 0,5раЬУ2- сила напора ветра, действующего на лопату, Н; Ь - плечо приложения силы боковой лопаты, м; аЬ - площадь лопаты, м2; р - плотность воздуха, кг/м ;
КЛ = 0,5ЬраЬ - коэффициент момента лопаты, 1/кг.
С другой стороны на виндрозу действуют противодействующие моменты системы пружин:
£МП = С0(у - уо) + С1(у - у±) + 62(7 - 72) + - + СП(У - Гп), (1.4)
где /0, /2 и /п - угловое смещение точек упора вторых концов соответственно нулевой (основной), первой, второй, и л-ной пружин, рад. При этом /0 < / <
/2 < - < /п;
С0, Сд С2 и Сл - угловые коэффициенты жёсткости соответственно основной, первой, второй и л-ной корректирующих пружин, численно равные моменту противодействия пружины при закручивании её на угол один радиан, Нм. В установившемся состоянии при У> Ур
Мл _ ТМП
или
КлУ2 _ ^оС/ - /о) + ^С/ - /1) + ^2С/ - /2) + - + СпС/ - /п) (1.5) Отсюда, задавшись коэффициентами жёсткости пружин и начальными углами их установки, можно определить скорость ветра, при которой поворот ветроколеса осуществится на требуемый угол у. При этом, если выражение в скобках будет отрицательно, то составляющая момента от соответствующей пружины и всех последующих принимаются равными нулю. Следовательно, в пределах изменения у от нуля до у1 (работает основная пружина):
У=ГК^ (1.6)
В пределах изменения у от у1 до у2 (работает основная и первая корректирующая пружина):
У _ 2 1со(7~7о) + с1(7~71) (1 7)
у/ кл
В пределах изменения у от у2 до у3 (работает основная, первая и вторая корректирующие пружины):
У _ 22^о(7~7о) + ^1(7~71) + ^2(7~72) (1 8)
у/ кл
В пределах изменения у от уп и выше (работает основная и л корректирующих пружин):
У _ 2|Со(7-7о) + С1(7-71) + С2(7-72) + -+Сп(7-7п) (1.9)
2 Кл
Полученные уравнения позволяют определить скорость ветра необходимую для поворота ветротурбины на фиксированный угол у. Так скорость ветра, при которой только что начнётся поворот виндрозы (у=0),
У° = Ур =
2 СоС-Го).
(1.10)
л! Кл '
Она же будет и расчётной скоростью ветра, а мощность ветротурбины при этой скорости - расчётной (номинальной) скоростью ветроагрегата.
Скорость ветра, при которой у станет равно ух (момент вступления в действие первой корректирующей пружины):
У± =
2 СоСп-Го)
Кл
Скорость ветра, при которой у станет равно у2 (момент вступления в действие второй корректирующей пружины):
У2 =
2 СоСГ2-Уо)+С1СГ2-У1)
Кл
Скорость ветра, при которой у станет равно уп (момент вступления в действие л-ной корректирующей пружины):
Уп =
2 СоСГп-Го) + С1СГп-Г1) + - + Сп-1СУп-Гп-1)
Кл
Тогда определить, как будет изменяться у в зависимости от скорости ветра можно будет по формулам:
КлУ2 + СоГо.
при У0<У<У1
У =
(111)
при У<У<У2
У =
КлУ2 + СоУо + С1Г1. Со + С1
(1.12)
при У2<У<У3
У =
КлУ2 + СоУо + С1У1 + С2У2. Со + С1 + в2
(1.13)
при Ул<У
У =
КлУ2 + СоУо + С1Г1 + С2Г2 + - + СпУп. Со ^ С1 ^ С 2 ^ "' ^ Сп
(1.14)
Так как угловое смещение точки крепления первой пружины с целью создания начального напряжения осуществлено в обратную сторону, то у± < 0.
Выражения позволяют подобрать такие значения входящих в них коэффициентов, при которых поворот виндрозы начнётся при расчётной скорости ветра. Например, можно предварительно задавшись величинами Кл и у0 определить нужную величину 60, или, наоборот, задавшись величинами Кл и 60, определить предварительную закрутку первой пружины у0.
По выведенным законам будет происходить фактическое изменение у при изменении скорости ветра. Подбором коэффициентов жёсткости пружин фактическое изменение у можно приблизить с заданной точностью к оптимальному изменению.
При расчётной скорости ветра и мощность ветроколеса будет считать расчётной Рр:
РР = КРУ?
(1.15)
Р = ККр Vр ,
где Кр - геометрический коэффициент мощности ветроколеса, зависящий от его диаметра, конструкции и количества лопастей, а также от плотности воздуха,
м2/кг.
Тогда мощность при скорости ветра больше расчётной и у=0:
Р = КРУ3,
Однако в связи с тем, что при скорости ветра больше расчётной происходит поворот ветроколеса на угол у, мощность ветроколеса снижается до величины:
Руу =Роу(^у)3, Вт (116)
Для того чтобы рассчитать параметры ветроколеса при различных скоростях ветра зададимся некоторыми постоянными величинами: У0 = 8 м/с; а = 0,2 м; Ь = 0,2 м; L = 1,2 м; у0 = —2,67рад (при угловом коэффициенте жесткости основной пружины 60=0,74 Нм); р = 1,29 кг/м3; КР = 50 м2/кг.
Отклонение мощности фактической РуУ при регулировании от мощности расчетной Р0У будем считать по формуле:
ДР% = ^-^100% (1.17)
Р0 V
Поочередно вводя корректирующие пружины и подбирая их угловые коэффициенты жесткости добьемся отклонения мощности ДР = ±5%. Фактическое отклонение ветроколеса при работе основной и корректирующих пружин рассчитаем по формулам (1.6)—(1.8). Оптимальный угол отклонения найдем по формуле (1.2). Начиная с расчетной скорости ветра и изменяя ее, найдем параметры ветроколеса при работе основной и корректирующих пружин (табл. 1).
Таблица 1
Расчёт характеристик ветроколеса__
V, м/с Мп, Нм Уф, рад Уопт, рад Pov, Вт PYv, Вт АР,% Уф, Град
8 1,98144 0 0 25600 25600 0 0
8,1 2,031286 0,067359 0,157297 26572,05 26391,68 3,092512 3,859382
8,25 2,107215 0,169966 0,246809 28075,78 26879,51 4,998087 9,738347
8,5 2,23686 0,345162 0,344701 30706,25 25587,28 -0,04968 19,77634
8,6 2,289802 0,416705 0,37575 31802,8 24317,19 -5,01099 23,87543
8,7 2,343362 0,430474 0,403887 32925,15 24710,82 -3,47337 24,66433
8,8 2,397542 0,444402 0,4297 34073,6 25077,88 -2,03952 25,46234
8,9 2,452342 0,458489 0,453598 35248,45 25416,4 -0,71721 26,26948
9 2,50776 0,472735 0,475882 36450 25724,38 0,485868 27,08574
9,1 2,563798 0,487141 0,496782 37678,55 25999,9 1,562103 27,91111
9,2 2,620454 0,501706 0,516475 38934,4 26241,04 2,504052 28,74561
9,3 2,67773 0,516429 0,535105 40217,85 26445,95 3,304476 29,58923
9,5 2,79414 0,546355 0,569621 42868,75 26740 4,45311 31,30383
9,6 2,853274 0,561556 0,585686 44236,8 26825,81 4,788325 32,1748
9,7 2,913026 0,576917 0,60105 45633,65 26868,77 4,956134 33,0549
9,9 3,03439 0,608116 0,62991 48514,95 26820,71 4,768387 34,84246
10 3,096 0,623954 0,643501 50000 26727,34 4,403653 35,74992
10,2 3,221078 0,656108 0,669207 53060,4 26397,27 3,114345 37,5922
10,3 3,284546 0,672423 0,681388 54636,35 26159,28 2,184701 38,52702
10,4 3,348634 0,688898 0,69316 56243,2 25872,14 1,063034 39,47096
10,5 3,41334 0,705532 0,704547 57881,25 25535,72 -0,25109 40,42402
10,6 3,478666 0,722325 0,715574 59550,8 25150,17 -1,75714 41,3862
10,7 3,54461 0,739278 0,72626 61252,15 24715,87 -3,45364 42,3575
10,75 3,577815 0,747814 0,731481 62114,84 24480,63 -4,37254 42,84657
10,78 3,597812 0,752954 0,734576 62636,33 24333,75 -4,94627 43,14111
11 3,74616 0,765865 0,756456 66550 24920,76 -2,65329 43,88085
11,1 3,814582 0,77182 0,765955 68381,55 25167,86 -1,68804 44,22204
11,2 3,883622 0,777829 0,775193 70246,4 25401,92 -0,77376 44,56632
11,3 3,953282 0,783892 0,784184 72144,85 25622,42 0,087573 44,91369
11,5 4,09446 0,796179 0,801469 76043,75 26020,77 1,643651 45,61768
12 4,45824 0,827839 0,841069 86400 26745,87 4,476062 47,4317
13 5,23224 0,895202 0,907923 109850 26865,35 4,94278 51,29131
14 6,06816 0,967954 0,962551 137200 25007,54 -2,3143 55,45969
14,2 6,242774 0,983151 0,972325 143164,4 24395,57 -4,70479 56,33042
16 7,92576 1,129625 1,047198 204800 15944,46 -37,717 64,72276
18 10,03104 1,312852 1,110242 291600 4840,427 -81,0921 75,2209
20 12,384 1,517636 1,159279 400000 60,00928 -99,7656 86,95412
При расчете дойдя до значения скорости ветра 8,6 м/с отклонение мощности от расчетной составило порядка 5 % в сторону убывания. При этом отклонение ветротурбины от направления ветра составило почти 24°. В этот момент включается в работу первая корректирующая пружина (Gi=3,15 Нм^ и постепенно происходит стабилизация мощности ветроколеса. При достижении скорости ветра значения 10,8 м/с будет задействована вторая корректирующая пружина (G2=7,6 Нм), посредством которой мощность ветроколеса останется в заданных пределах. Угол отклонения ветротурбины составит при этом порядка 43°. При отклонении ветроколеса на угол порядка 1 радиана и больше коэффициент использования энергии ветра резко падает, что вызывает снижение вырабатываемой мощности. Это доказывают и результаты расчетов: при V>14 м/с значение фактической мощности уменьшается и при V=20 м/с составляет около 60 Вт. Ветроколесо полностью выводится из под ветра о чем говорит его угол отклонения уф=86°.
По полученным результатам расчёта построим зависимости мощности ветроколеса и угла отклонения от скорости ветра.
А
Р, кВт
зо
А 7- 2Рад
Р=№
20 80
60
10 М
20
О
18 К м/с
Рис. 2. Характеристики ветроколеса в результате работы противодействующих пружин
В интервале скоростей ветра от 0 до 8 м/с мощность ветроколеса изменяется
устройство виндрозы, которое поддерживает мощность на заданном уровне
). Пунктирной линией показано изменение мощности при отсутствии устройства ограничения. С достижением скорости ветра значений около 14 м/с мощность ветроколеса падает.
Список литературы:
1. Адрианов, В. Н. Ветроэлектрические станции [Текст] / В. Н. Адрианов, Д. Н. Быстрицкий - М.: Государственное Энергетическое Издательство, 1960. - 310 с.
2. Д. де Рензо. Ветроэнергетика [Текст] / Д. де Рензо - М.: Энергоатомиздат, 1982. - 272 с.
3. Фатеев Е. М. Ветродвигатели и ветроустановки [Текст] / Е. М. Фатеев - М.: ОГПЗ-СЕЛЬХОЗГИЗ, 1956. - 544 с.
по закону
. При скорости ветра больше 8 м/с включается в работу
(
) путем поворота ветроколеса на определенный угол (зависимость
4. Харитонов, В. П. Автономные ветроэлектрические установки [Текст] / В. П. Харитонов. - М.: ГНУ ВИЭСХ, 2006. - 280 с.
