Зубарев Д.В.
Московская Академия Приборостроения и
Информатики,  г. Москва, Россия
Ветровой поток, кинетическую
энергию которого ВЭУ преобразует в электрическую, изменяется в крайне широких
пределах, при этом мощность его пропорциональна кубу скорости ветра, а силовое
воздействие на конструкцию ВЭУ пропорционально квадрату скорости ветра. По этим
причинам создание работоспособной, надежной и приемлемой по стоимости ВЭУ без
мехатронных элементов  управления невозможно.  
В
качестве движительной установки ВЭУ выбрано быстроходное ветроколесо с тремя
поворотными лопастями, диаметром 10 метров и номинальной быстроходностью Z ном = 8. Данное ветроколесо обеспечивает мощность 16
кВт при скорости ветра 7.5- 10 м/с. Расположение ветроколеса по направления к
ветру за башней. Управление вращением ветроколеса производится за счет поворота
лопастей. Рабочий режим ветроэнергетической установки 100 об/мин. Диапазон
поворота лопастей 0..90 градусов.
Для
управления ВЭУ используется микропроцессорная САУ, которая позволяет за счет
поддержания оптимальных оборотов и потребляемой мощности достичь  максимального преобразования энергии ветра в
электрическую энергию. В этой статье приведена модель функционирования ВЭУ.
.doc_files\image002.jpg)
Рис. 1.  Модель
управления ВЭУ  в рабочем режиме
          Основное дифференциальное уравнение,
описывающее динамику вращательного движения ветроколеса представляется в виде:
.doc_files\image004.gif)
,где Мвала = Ma – Mc
Где, Ma – аэродинамический момент ВК;
        Mc – момент сопротивления;
        Jвк  
–  момент инерции системы ВК + МЧВ
+ Генератор , приведенный к оси ВК;
          ε – угловое ускорение ВК;
аэродинамический
момент ВК равен:
.doc_files\image006.gif)
Момент
сопротивления на рабочем режиме:
.doc_files\image008.gif)
Момент
сопротивления при разгоне или торможении определяется потерями на трение в МЧВ
и генераторе.
Mc = 0.05*Mном
где, ω – частота вращения ВК
       P – активная электрическая мощность нагрузки;
       Mном – номинальный крутящий момент ВК.
В
результате решения дифференциального уравнения получаем частоту вращения ВК – w
Вычисление
углового ускорения ВК по 7-и точкам основывается на аппроксимации угловой
скорости вращения ВК по линейному закону и вычислении коэффициентов в законе по
методу наименьших квадратов. Принимая, что на интервале  .doc_files\image010.gif) вращение ВК является
равноускоренным, тогда закон вращения ВК запишется:
 вращение ВК является
равноускоренным, тогда закон вращения ВК запишется:
                   .doc_files\image012.gif)
          Период опроса .doc_files\image014.gif) датчика оборотов ВК
равен удвоенному периоду следования окон диска датчика:
 датчика оборотов ВК
равен удвоенному периоду следования окон диска датчика:
                   .doc_files\image016.gif)
Коэффициенты
.doc_files\image018.gif) и
и .doc_files\image020.gif) определяются по методу
наименьших квадратов:
 определяются по методу
наименьших квадратов:
                   .doc_files\image022.gif)
                   .doc_files\image024.gif) ,
,
где              .doc_files\image026.gif)
          Принимая в .doc_files\image014.gif) 
  .doc_files\image028.gif) получим следующую
формулу для вычисления текущего значения углового ускорения ВК:
 получим следующую
формулу для вычисления текущего значения углового ускорения ВК:
          .doc_files\image030.gif)
На
рабочем режиме, когда осуществляется стабилизация заданной частоты вращения ВК,
управляющий сигнал на ЭД формируется следующим образом:
                   .doc_files\image032.gif)
                   .doc_files\image034.gif)
где    S -
управляющий сигнал на ЭД, рад/с;
          Кw  – коэффициент усиления по частоте вращения
ВК;
          Кe  – коэффициент усиления по угловому ускорению
ВК, с;
          wном=10.47рад/с – заданная номинальная
частота вращения ВК;
          w - текущее значение частоты вращения ВК с учетом фильтрации и
ограничения, рад/с;
          e - фильтрованное и ограниченное угловое ускорение ВК, рад/с2;
                   .doc_files\image036.gif)
где    .doc_files\image038.gif) =314рад/с – максимальная частота вращения ЭД соответствующая
максимальному значению управляющего сигнала на ЭД.
=314рад/с – максимальная частота вращения ЭД соответствующая
максимальному значению управляющего сигнала на ЭД.
         В соответствии с данным законом
управления осуществляется формирование, ограничение и выдача управляющего
напряжения на вход ЭД, которое задает направление и скорость вращения ЭД
привода поворота лопастей ВК.
         Далее исходя из уравнения движения
электропривода поворота лопастей имеем:
Ток
управления  .doc_files\image040.gif) где S  - значение
кода, вычисляемое микропроцессором САУ и задающее угловую скорость вращения
электродвигателя.
 где S  - значение
кода, вычисляемое микропроцессором САУ и задающее угловую скорость вращения
электродвигателя.
заданное
напряжение на входе электродвигателя:
.doc_files\image042.gif) где,
 где,
.doc_files\image044.gif) - коэффициент делителя на напряжения, который ограничивает
максимальную угловую скорость вращения электродвигателя.
- коэффициент делителя на напряжения, который ограничивает
максимальную угловую скорость вращения электродвигателя.
.doc_files\image046.gif) где,
 где,
.doc_files\image048.gif) -
коэффициент тахогенератора (обратной связи по угловой скорости вращения
электродвигателя).
 -
коэффициент тахогенератора (обратной связи по угловой скорости вращения
электродвигателя).
.doc_files\image050.gif) - угловая
скорость вращения электродвигателя.
 - угловая
скорость вращения электродвигателя.
Управляемый
электромагнитный момент двигателя:
.doc_files\image052.gif) где,
 где,
.doc_files\image054.gif) -  коэффициент двигателя
 -  коэффициент двигателя
.doc_files\image056.gif) - ток в якорной цепи
электродвигателя
 - ток в якорной цепи
электродвигателя
момент
двигателя:
.doc_files\image058.gif)
.doc_files\image060.gif)
двигатель
преодолевает шарнирный момент создаваемый лопастью:
.doc_files\image062.gif) - шарнирный момент
одной лопасти где,
 - шарнирный момент
одной лопасти где,
.doc_files\image064.gif) - быстроходность ВК
 - быстроходность ВК
.doc_files\image066.gif) - скоростной напор
ветрового потока отсюда шарнирный момент всего ВК:
 - скоростной напор
ветрового потока отсюда шарнирный момент всего ВК:
.doc_files\image068.gif) где,
 где,
.doc_files\image070.gif) - передаточный
коэффициент МПЛ и моментом трения:
 - передаточный
коэффициент МПЛ и моментом трения:
.doc_files\image072.gif)
.doc_files\image074.gif) где,
 где,
.doc_files\image076.gif)
В
установившемся режиме из условия  .doc_files\image078.gif) получаем:
  получаем:
.doc_files\image080.gif) где
где  .doc_files\image082.gif) - в град/с.
 - в град/с.
Литература:
1. Е.М
Фатеев. Ветродвигатели и ветроустановки. Москва 1948 г.
2 .Подураев Ю.В., Кулешов В.С. Принципы построения и
современные тенденции
развития мехатронных систем.
Мехатроника. 2000. №1. С.5-10.
         3. Аршанский М.М., Шалобаев Е.В. Мехатроника: основы
глоссария. 
Мехатроника. 2001. № 4.
С.47-48.
4. ГОСТ
Р51237-98