CC BY
142
УДК 620.91
ОПРЕДЕЛЕНИЕ РЕЗОНАНСНЫХ ДИАГРАММ ЛОПАСТЕЙ МАЛОГАБАРИТНЫХ ВЕТРОЭНЕРГЕТИЧЕСКИХ УСТАНОВОК А.С. ГУРИНОВ, В.В. ДУДНИК
(Донской государственный технический университет)
Представлена методика определения резонансных характеристик лопасти малогабаритной ветроэнергетической установки. Для вегротурбины диаметром 4,5 м были построены расчетные формы первых трех тонов колебаний лопастей и построена резонансная диаграмма.
Ключевые слова: ветроэнергетическая установка, резонансная диаграмма, формы колебаний.
Введение. Источники возобновляемой энергии - один из современных способов получения электроэнергии с использованием природных явлений. Ветроэнергетические установки среди прочих нетрадиционных источников получили самое широкое распространение. Ветроэнергетика вызывает повышенный интерес по целому ряду причин. Главной из них является независимость ветро-генератора от общей системы электроснабжения, отсутствие выбросов. При правильном выборе ветрогенератора его использование приносит существенную экономию средств. Это особенно актуально при долгосрочной эксплуатации.
Наилучшие характеристики показывают ветроэнергоустановки (ВЭУ), обладающие системами автоматического управления. В случае использования системы автоматического управления с нагрузкой необходимо добиваться изменения частоты вращения при частотах, близких к резонансным. Только в этой ситуации можно добиться длительной работы ВЭУ. Наиболее опасны резонансные колебания лопастей в направлении перпендикулярном плоскости вращения. Частоты таких колебаний могут быть определены аналитически.
Методика построения резонансных диаграмм. Ветротурбина, создавая крутящий момент на валу генератора, создает и переменные нагрузки, обусловленные как наличием косой обдувки, так и дисбалансом самого ветроколеса. Амплитуда переменных нагрузок зависит от конструктивных особенностей лопасти и режимов работы ветроэнергетической установки. Кроме того, амплитуда колебаний лопастей зависит от того, насколько близки частоты внешних сил, возбуждающих колебания к собственным частотам лопастей. При совпадении частот происходит резкое увеличение амплитуды, обусловленное резонансом. Для определения резонансов ветротурбины необходимо построить резонансную диаграмму колебаний в плоскости поперечной плоскости вращения.
В настоящее время существует ряд методов расчета резонансных характеристик лопастей. Однако для малых ВЭУ, самой нижней ценовой категории, необходим наиболее простой метод, расчет которого представлен ниже.
В общем случае уравнение свободных колебаний лопасти в этой плоскости имеет вид [2]
(Ео:)"- (ад)' + т у = о, (1)
где Е: - модуль упругости материала лонжерона лопасти; I х - момент инерции сечения лонжерона лопасти относительно горизонтальной оси; у: - перемещение лопасти в вертикальном направлении; Ыцб - центробежная сила; тл - погонная масса лопасти.
Штрихами обозначено дифференцирование по радиусу, а точками - дифференцирование по времени.
При использовании жесткой заделки лопастей можно принять следующие граничные условия
у ( 0) = 0, у"( 0) = 0, (Е I у") в = 0, (Е I у")' в = 0.
У :(г=0) • у :(г=0) • V : лхУл'ґ =К • V : лхУл'ґ =К
Решение уравнения (1) можно представить в виде произведения двух функций
:ил (г, ^ = £(г) .5(0, (2)
где £(г) - функция только радиуса; 5(0 - функция только времени.
Соответственно выражение распадается на уравнения
5(') + Рлх25(() = 0, (3)
где рлх - угловая частота собственных изгибных колебаний лопасти.
[Е1 £"(г)] " - ["(г)] ' - рлх2т£(г) = 0. (4)
Для решения этих уравнений наиболее распространенным является способ Б.Г. Галеркина [3], согласно которому решение ищется в виде ряда
£(г) = £СФ.(г), (5)
i=0
где Ci - постоянные коэффициенты; ф.(г) - заранее заданные функции радиуса, удовлетворяющие граничным условиям; п - количество членов ряда, принимаемых к расчету.
Функции ф. (г) выбираются так, чтобы каждая из них удовлетворяла всем граничным условиям, а вся совокупность функций ф. (г) сходилась к какому-то значению, а не расходилась. Можно воспользоваться функциями Крылова при соответствующем выборе коэффициентов.
Согласно методу Галеркина ряд (5) подставляется в равенство (4), каждый член умножается на одну из функций ф. (г) и интегрируется по радиусу лопасти
±С } (Ел1тф ") 'ф^г -£ С. | (Мцб ф.) " ф^г -рлх £ |тл ф. Ф^г = 0, (6)
.=0 0 .=0 0 .=0 0
где . - номер функции в ряде (5); k - номер функции на которую умножается уравнение (4).
Для вычисляемых интегралов можно ввести обозначение
Я
| (ЕАЖ) ''фidr = сй
0
1 (N6 ф ;) ' ф^г = Ьік , (7)
| тлфгф^г = тік
0
Таким образом, вместо уравнения (4) можно перейти к системе из п алгебраических уравнений вида (6). Система состоит из бесконечного количества уравнений, но на практике целесообразно ограничиться 4-5 членами и, соответственно, 4-5 уравнениями, которые выглядят следующим образом
£с (с* -ь1к -р1т1к) = 0. (8)
.=0
Форма колебаний в этом случае определяется с точностью до неопределенного множителя с0. Этот множитель может быть найден из условий нормирования. Наиболее распространенное условие нормирования - равенство единице наибольшего отклонения от начального положения. Форма поперечного колебания - это кривая изгиба оси балки в момент наибольшего отклонения от нейтрального положения. При собственных колебаниях внешняя сила отсутствует, поэтому амплитуду колебаний определить невозможно, что и выражается в виде множителя с0 .
При расчетах ряд (5) искусственно обрывают [1]. Ввиду этого, величины р и являются
приближенными, зависящими от количества членов ряда. С повышением использования количества членов ряда точность расчетов растет.
Преобразованный вид формул для частот собственных колебаний дает следующее решение
Р* = — + ®1 —, (9)
где
т»
К
mkk = | mAldr,
0
-*к ^ ($к )2 dr;
0
Скк = 1^ (% )2 dr,
0
К
Sл =| тл Мг .
Г
В случае если юнв = 0, то вращение не происходит, т.е. квадрат частоты собственных колебаний невращающейся лопасти определяется по формуле
р2 =— = р! . (10)
mkk
Для упрощения можно ввести коэффициент
А*
кк=-кк-. (11)
ткк
Соответственно формула (9) преобразуется к виду
рлк = Рквт ^ кк®ет . (12)
Таким образом, собственные частоты колебаний лопасти можно определить, суммируя частоту собственных колебаний невращющейся лопасти и изменение частоты, вызванное увеличением жесткости вследствие центробежной силы. Такой подход пригоден для определения собственных частот колебаний лопасти. Коэффициенты ткк, -кк, ск вычисляются по формам собственных колебаний невращающейся лопасти, также как Ркет. Согласно [1], данный метод может обеспечить погрешность 1-2%.
Выражение (12) показывает, что с ростом угловой скорости ветротурбины собственные частоты изгибных колебаний лопасти повышаются. Физически это объясняется действием растущей центробежной силы, стремящейся распрямить лопасть, т.е. увеличить ее жесткость на изгиб. Как известно, при повышении жесткости собственные частоты колебаний растут.
Уравнение собственных колебаний невращающейся лопасти вытекает из уравнения (1) с учетом отсутствия центробежной силы
(Ел1х $ " ) " - рЛтД = 0. (13)
При незначительном изменении удельной массы и жесткости по длине лопасти уравнение принимает вид
-«Л^ = 0 , (14)
л •
2,
где а; = ЬЪ..
Ел лх
Общее решение данного уравнения, выраженное через функции Крылова, имеет вид
£(г) = AS(аг) + ВТ(аг) + Си(аг) + DU(аг), (15)
где А, В, С, D - постоянные, определенные граничными условиями; S, Т, и, V- функции Крылова
Технические науки
S (ar) = 0,5(ch ar + cos ar)
T(ar) = 0,5(sh ar + sin ar)
. (16)
U(ar) = 0,5(ch ar - cos ar)
V(ar) = 0,5 (sh ar - sin ar)
При r = 0, S (0) = 1, а T (0) = V (0) = U (0) = 0. Производные от одной из функция Крылова выражаются через другие:
— = 0,51 d ch ar + d cos a ) = 0,5a(sh ar - sin ar) = aV,
dr ^ dr dr )
dT = 0,51 d sh ar + dsin a ) = 0,5a(ch ar + cos ar) = aS, (17)
dr ^ dr dr )
d2 S dV 2tt
—— = a = a U.
dr dr
Таким образом, любая производная функции может быть найдена путем круговой замены.
Для определения постоянных A, B, C, D используют четыре граничные условия [4]. Из
(14) поочередно находятся производные и подставляются в соотношения при r=0 и r=R. Это дает
систему уравнений относительно неизвестных A, B, C, D. Такая система имеет нетривиальное решение. Тривиальное решение A=B=C=D=0 отвечает только состоянию покоя лопасти и не представляет интереса. Для нахождения нетривиального решения определитель, составленный из коэффициентов при неизвестных A, B, C, D, приравнивается к нулю. Так получают уравнение для нахождения a;, каждому корню которого соответствует свое нетривиальное решение исходной системы алгебраических уравнений. Подставляя последовательно каждое значение a; в систему алгебраических уравнений, определяются соотношения между коэффициентами A, B, C, D, т.е. соответствующую форму колебаний. Число собственных колебаний балок бесконечно, однако практический интерес представляет собой три-четыре значения a;.
Таким образом, частоты равны
P =а2
E I
. (18)
1
Приравняв определитель к нулю, получаем выражение для нахождения а;
Л(аК) • cos(аК) +1 = 0 . (19)
Его решения равны
1,571
а, =------,
1 К
Л,712
а=
R (20)
7,854
R
10,996
4 R
Значение a; может быть определено приближенно по формуле
j - 0,5
а « л---------------. (21)
’ R
Эта формула тем точнее, чем выше номер тона i.
Определение форм колебаний производится с учетом граничных условий винта
S (а/^ (а1К)
и (а/) --
(22)
Т (а,Л)
Для построения форм колебаний целесообразно воспользоваться условиями нормирова-
ния, в частности, равенства величины перемещения конца лопасти единице:
л.
^( К) = 1, |^2(г ^г = 1.
Соответственно производные равны
£,' (г) = BаS + СаТ + DаU
£, '' (г) = Bа2V + Са2 S + Dа2T £, '' (0) = Ва^ (0) + Са2 S (0) + Dа2T (0) = С а2 -1 = 0 Приведенная масса может быть выражена
| m^2jHdr = т |Ъ^^г = тА ,
(23)
(24)
(25)
где А =
А = —г.
Численное интегрирование можно выполнить методом Симпсона:
Лг (С + 4С2 + 2С3 + 4С4 + 2С5 +... + С11)
А = -
3,
(26)
Для вычисления Ь*. предварительно необходимо найти квадрат первой производной соб-
ственной формы колебаний:
ЬЛ- =| (| т^г)(£,'н )2dr = | т(|Ыг)(£,'н )2dr ,
(27)
где С- =^2, Лг - расстояние между расчетными сечениями лопасти. Резонансные параметры лопасти. На основании представленной методики был запрограммирован расчет колебаний лопастей малогабаритной ветроэнергетической установки, разрабатываемой ООО «Персональные энергосистемы» (г. Ростов-на-Дону). Внешний вид установки представлен на рис.1. В расчете рассматривалась лопасть радиусом 2,25 м со стальным трубчатым лонжероном с жесткой заделкой в комле.
Трубчатый лонжерон выполнял роль несущего элемента. К лонжерону крепился пенопластовый заполнитель с обшивкой. Согласно аэродинамическому расчету оптимальная частота вращения составила 380 об/мин.
Расчет велся до частоты вращения 420 об/мин. Интегрирование проводилось вдоль лопасти с шагом 0,1 II.
На рис.2 показаны расчетные формы первых трех тонов колебаний лопасти, приведенные по условию нормирования к 1.
Рис.1. Внешний вид малогабаритной ветроэнергетической установки
Рис.3. Резонансная диаграмма лопастей малогабаритной ветроэнергетической установки
Как видно из резонансной диаграммы, представленной на рис.3, лопасти имеют хорошие резонансные характеристики, т.е. наступление резонанса возможно только на одной частоте, далекой от оптимальной. Для более наглядного представления значения по оси абсцисс представлены в относительных значениях, приведенных к оптимальной частоте вращения:
где ропт - оптимальная частота вращения ветроколеса.
На данном рисунке кривые тонов колебаний представляют собой графики собственных частот колебаний лопасти при увеличении центробежной силы. Наклонные прямые - графики частот вынужденных колебаний при различной частоте вращения ветроколеса при разных скоростях ветра. На рисунке видно, что первый тон собственных колебаний имеет низкие резонансные частоты, что характерно для всех ветротурбин, у которых лопасти имеют большой удельный вес и низкую жесткость. Пересечение первого тона с гармониками возбуждающей силы на частоте вращения менее 200 об/мин не представляет проблем, так как на этой частоте управление нагрузкой не осуществляется.
Существует проблема с пересечением второго тона с пятой гармоникой на частоте 297 об/мин. Для защиты от возможных колебаний в этой зоне необходимо добавить в алгоритмы управления требования или ослабления нагрузки, с целью выхода ветротурбины на частоты выше 305 об/мин или увеличения, с целью понижения частот вращения ниже 290 об/мин.
Проверить правильность выполнения расчета экспериментально можно, как правило, в одной точке - при отсутствии вращения. В этом случае на лопасти устанавливается виброгенератор малой массы и производится замер частоты появления резонанса по резкому росту амплитуды колебаний. Создать распределенное усилие, аналогичное нагрузкам на ветроколесе при вращении и сгенерировать колебания для нахождения резонансов представляется проблематичным. Выводы. Таким образом, представленный метод может быть использован для определения резонансных характеристик лопасти ветроустановок, имеющих малое изменение массы и жесткости вдоль лопасти.
Расчет, выполненный для двухлопастной ветротурбины диаметром 4,5 м, показал, что не допустить появления резонанса на данном типе лопастей можно, использовав в алгоритмах управления нагрузкой ограничения рабочего диапазона частот 200-290 об/мин и 305-420 об/мин. Библиографический список
1. Кирпикин А.А. Приближенный расчет резонансных диаграмм несущих винтов вертолетов / А.А. Кирпикин. - Харьков: ХАИ, 1999. - 64 с.
2. Миль М.Л. Вертолеты. Расчет и проектирование / М.Л. Миль [и др.]. - Т.2. Колебания и динамическая прочность. - М.: Машиностроение, 1967. - 424 с.
3. Михеев Р.А. Расчет вертолета на прочность. Прочность лопастей несущего винта /
Р.А. Михеев. - М.: МАИ, 1973. - 205 с.
4. Скопинцев Б.И. Собственные колебания лопасти вертолета / Б.И. Скопинцев. - Харьков: ХАИ, 1983. - 200 с.
Материал поступил в редакцию 22.06.2011.
References
1. Kirpikin A.A. Priblizhyonny'j raschyot rezonansny'x diagramm nesushhix vintov vertolyotov / A.A. Kirpikin. - Xar' kov: XAI, 1999. - 64 s. - In Russian.
2. Mil' M.L. Vertolyoty'. Raschyot i proektirovanie / M.L. Mil' [i dr.]. - T.2. Kolebaniya i
dinamicheskaya prochnost'. - M.: Mashinostroenie, 1967. - 424 s. - In Russian.
3. Mixeev R.A. Raschyot vertolyota na prochnost'. Prochnost' lopastej nesushhego vinta / R.A. Mixeev. - M.: MAI, 1973. - 205 s. - In Russian.
4. Skopincev B.I. Sobstvenny'e kolebaniya lopasti vertolyota / B.I. Skopincev. - Xar'kov: XAI, 1983. - 200 s. - In Russian.
SPECIFYING RESONANCE DIAGRAMS OF COMPACT AEROGENERATOR BLADES A.S. GURINOV, V.V. DUDNIK
(Don State Technical University)
The resonance characteristics analysis technique of the compact aerogenerator blade is presented. The design forms of the first three blade vibration tones and a resonance diagram have been built for the wind turbine of 4.5 m in diameter.
Keywords: aerogenerator, resonance diagram, vibration modes.
