автореферат диссертации по энергетике, 05.14.08, диссертация на тему:Методология разработки и создания вертикально-осевых ветроэнергетических установок

Российская академия сельскохозяйственных наук Государственное научное учреждение Всероссийский научно-исследовательский институт электрификации сельского хозяйства

005010431

Соломин Евгений Викторович

На правах рукописи

С,

МЕТОДОЛОГИЯ РАЗРАБОТКИ И СОЗДАНИЯ ВЕРТИКАЛЬНО-ОСЕВЫХ ВЕТРОЭНЕРГЕТИЧЕСКИХ УСТАНОВОК

Специальность 05.14.08 - энергоустановки на основе возобновляемых видов энергии

Автореферат диссертации на соискание ученой степени доктора технических наук

9 ОЕЗ 20М

Москва

2011

005010431

Работа выполнена в ФГБОУ ВПО «Южно-Уральский государственный университет» (Национальный исследовательский университет) на кафедре «Электротехника и возобновляемые источники энергии».

Научный консультант

доктор технических наук, профессор И.М. Кирпичникова.

Официальные оппоненты:

доктор технических наук, профессор П.П. Безруких;

доктор технических наук, профессор

В.В. Елистратов;

доктор технических наук, профессор Ю.Г. Шакарян.

Ведущее предприятие -

Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Уральский федеральный университет имени первого Президента России Б.Н.Ельцина» (УрФУ), г. Екатеринбург.

Защита состоится «-2&> 20(Ь_ г., в /£_ Р° часов, на заседании

диссертационного совета Д 006.037.01 по присуждению учёной степени доктора технических наук при ГНУ «Всероссийский научно-исследовательский институт электрификации сельского хозяйства» (ГНУ ВИЭСХ РАСХН) по адресу: 109456, Москва, 1-й Вешняковский проезд, 2. Факс: (499)170-51-01. E-mail: viesh@dol.ru. Отзывы на автореферат просим высылать по указанному адресу.

С диссертацией можно ознакомиться в библиотеке ГНУ ВИЭСХ РАСХН.

Автореферат разослан « 9 И 20f£ г.

Учёный секретарь диссертационного совета, доктор технических наук, профессор

А.И. Некрасов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы. Внедрение автономных ветроэнергетических установок (ВЭУ) в условиях континентального климата может обеспечить полным или частичным энергоснабжением небольшие поселки, жилые объекты, фермы, офисные помещения или малые цеха с целью развития инфраструктуры на удаленных территориях страны. В настоящее время разработка ВЭУ, перспективных для агропромышленного комплекса РФ, основана в целом на ряде инженерных подходов к проектированию отдельных компонентов, что свидетельствует об отсутствии научно обоснованной методологии разработки конструкции ВЭУ как единого целого. Оптимизация параметров установок затруднена в связи с наличием возмущений, вызванных взаимовлиянием компонентов друг на друга. Создание компонентов осуществляется разными производителями, что осложняет их сопряжение и монтаж. Эти факты являются крупной научной проблемой еще и в связи с тем, что вопросы разработки новых недорогих конструкций, а значит, и ускоренного освоения новых территорий решаются чрезвычайно медленно. Кроме того, с ростом энергодефицита и износом генерирующих мощностей централизованное энергообеспечение существующих и вновь возникающих автономных удаленных объектов становится еще более затруднительным.

Решение проблем проектирования оптимальных технологичных конструкций заключается в создании методологии разработки на основе системного научного подхода с применением совокупности методов и внедрении научно обоснованных технических, технологических и экономических решений. Методология позволит структурировать и логически организовать методы и средства оптимального проектирования. Такой подход имеет важное хозяйственное значение, а внедрение обоснованных решений вносит значительный вклад в развитие экономики страны и повышение ее энергонезависимости и обороноспособности.

Ряд уникальных исследований, проведенных известными Российскими учеными в ХХ-ХХ1 вв. (В.Н. Андриянов, П.П. Безруких, Д.Н. Быстрицкий, В.Р. Вашкевич, В.П. Ветчинкин, В.И. Виссарионов, Я.Б. Данилевич, Г.И. Денисенко, В.В. Елистратов, Н.Е. Жуковский, Е.И. Куклин, М.В. Кузнецов, В.А. Минин, Г.Х. Сабинин, Е.М. Фатеев, В.В. Харитонов, Я.И. Шефтер, Ю.Г. Шакарян и др.) позволяет утверждать, что данная проблема может быть с успехом решена в России.

Исследования, приведенные в диссертационной работе, трижды поддержаны международными грантами по программам финансирования гражданских научноисследовательских работ оборонного комплекса России, осуществляемой Международным научно-техническим центром, г. Москва (2004-2010 гг.), трижды - Минобрнауки РФ, г. Москва (2009-2011 гг.), грантом РФФИ (2011 г.), четыре раза -грантами Администрации г. Челябинска (2009-2011 гг.), а также средствами инвестиционных, инновационных и венчурных проектов в области ветроэнергетики.

Цель работы: создание и реализация методологии разработки вертикальноосевых ветроэнергетических установок с дальнейшей оптимизацией, производством компонентов и установок в целом.

Задачи исследований:

- обосновать необходимость методологического подхода к разработке и созданию вертикально-осевых ветроэнергетических установок (ВО ВЭУ);

- смоделировать и оптимизировать конструкцию ВО ВЭУ покомпонентно с целью получения максимальной выработки энергии;

- разработать технологические процессы изготовления компонентов ВО ВЭУ и методику вибробалансировки;

- провести экспериментальные исследования разработанных ВЭУ, в том числе измерения шума и вибраций;

- обосновать экономические и социальные перспективы массового использования малых ветроэнергетических установок на территории Российской Федерации;

- создать и внедрить семейство малых ВО ВЭУ автономного применения.

Объект исследования: вертикально-осевые ветроэнергетические установки с вновь разработанным модифицированным ротором Дарье.

Предмет исследования: эксплуатационные и конструктивные параметры вертикально-осевой ветроустановки и ее компонентов с учетом закономерностей их функционирования и наличия разнородных возмущающих факторов.

Обоснованность и достоверность научных положений, выводов и рекомендаций подтверждается аргументированностью исходных посылок, вытекающих из основ физики, аэромеханики и электротехники, удовлетворительным совпадением результатов теоретических исследований с результатами экспериментов и опытом эксплуатации.

Научная новизна: разработаны методологические основы создания вертикально-осевых ветроэнергетических установок. Научно обоснован системный подход при проектировании и оптимизации ВО ВЭУ с применением совокупности методов. Разработаны математические и физические модели ВО ВЭУ с целью оптимизации их параметров и режимов использования. Установлены закономерности влияния внешних и внутренних силовых воздействий, приводящих к возникновению шума и вибраций, разработана методика снижения уровня вибраций на основе оптимизации параметров ВЭУ в период разработки. Введено аэродинамическое регулирование частоты вращения ротора ВЭУ. Оптимизированы параметры генерирующей системы. Реализован алгоритм регулирования мощности ВЭУ.

Практическая ценность и реализация результатов исследований.

Разработан и внедрен в производство ряд технологических процессов изготовления комплектующих, разработана и реализована методика предэксплуатаци-онной статической и динамической балансировки роторов вертикально-осевых ВЭУ с целью снижения шума и вибраций.

Создано, испытано, внедрено и находится в опытной и коммерческой эксплуатации более 50 вертикально-осевых ветроэнергетических установок (ВО ВЭУ) с

3, 4 и 6-лопастными роторами мощностью от 0,1 до 30 кВт, предназначенных для работы в широком диапазоне скоростей ветра от 2 м/с до 45 м/с.

Результаты исследований используются в учебном процессе ФГБОУ ВПО «Южно-Уральский государственный университет» (г. Челябинск), ФГАОУВПО «УрФУ» (г. Екатеринбург), Государственном университете Северной Калифорнии (г. Беркли, США), в исследовательских программах Национальной Лаборатории Лоуренс Беркли (г. Беркли, США), в государственных проектах мэрии г. Сан-Франциско (США), мэрии г. Прайя (Кабо Верде) и т.д.

Рекомендации, сформулированные на основе научных положений, результатов и выводов, а также разработанные технологические процессы и методики ис-

пользуются рядом предприятий - производителей комплектующих (ООО «ГРЦ-Вертикаль», ОАО «Государственный Ракетный Центр, КБ им. Макеева», ОАО «КумАПП», Empire Magnetics. Inc. (США), WindSail, Inc. (США) и др.).

На базе ООО «ГРЦ-Вертикаль» создано опытное производство ВЭУ-3 (номинальная мощность 3 кВт). За 2009-2011 гг. реализовано свыше 20 коммерческих установок, показывающих прогнозируемые результаты.

Внедрение результатов исследований.

В настоящее время в России и за рубежом на различных испытаниях и в опытной эксплуатации находятся следующие модификации разработанных ветроэнергетических установок: ВЭУ-0.1 (0,1 кВт, 2 шт), ВЭУ-1 (1 кВт, 4 шт), ВЭУ-3(4) (3 кВт, 1 шт), ВЭУ-З(б) (3 кВт, 5 шт), ВЭУ-5(6) (5 кВт, 1 шт), ВЭУ-30(6) (30 кВт, 1 шт).

Рекомендации и проекты изменений ряда ГОСТ и стандартов организаций по ветроэнергетике переданы в ЗАО «НПЦ Малой Энергетики», г. Москва. Положения, выносимые на защиту.

- Методология разработки и создания ВО ВЭУ на основе поэтапной итерационной оптимизации параметров компонентов.

- Общие принципы построения математических моделей компонентов ВЭУ с учетом внешних и внутренних возмущающих факторов.

- Система аэродинамического регулирования ВО ВЭУ.

- Система регулирования мощности ВО ВЭУ.

- Технологические процессы изготовления лопастей и генератора ВО ВЭУ.

- Методика вибробалансировки.

- Семейство ВО ВЭУ с модифицированными роторами «Дарье».

- Обоснование экономических и социальных перспектив массового внедрения ВО ВЭУ в экономику Российской Федерации.

Апробация работы: разработки награждены двумя дипломами конкурсов разработчиков РФ, четырьмя дипломами лучших проектов РФ, двумя медалями за участие в региональных международных выставках. Разработки получили положительную оценку на четырех НТС различного уровня, включая областной, региональный и федеральный (Минэнерго РФ, Минобразования РФ), пяти межведомственных и международных совещаний в России, Финляндии, США.

Результаты работы были доложены, рассмотрены и одобрены на 16 научнопрактических конференциях регионального, российского и международного уровня, в том числе: на совещании в Министерстве Энергетики РФ о внедрении ВЭУ с водородными накопителями в ряде регионов Российской Федерации, 2009 г.; 21-м научном форуме Международного научно-технологического центра, г. Пусан, Ю.Корея, 2009 г.; Международном совещании комитета ООН по защите окружающей среды, посвященном развитию Африканских стран, г. Прайя, Кабо Верде, 2010 г.; научной сессии Академии электротехнических наук РФ по проблемам «Нетрадиционные и возобновляемые источники электроэнергии», г. Москва, 2010 г.; Международной энергетической конференции «Технологическая основа формирования новой энергетики России» в Московской школе управления Сколково, г. Москва, 2010 г.; Международном форуме инноваций России, Лаппеенран-та, Финляндия, 2010 г.; ежегодных (2007-2011 гг.), научных конференциях Южно-Уральского государственного университета, г. Челябинск и др.

Соискатель является членом Комитета по проблемам использования ВИЭ РосСНИО (Россия), членом научного совета международного научного журнала «Альтернативная энергетика и экология» (Россия), консультантом по альтернативным источникам энергии Уэйнского Государственного Университета (Детройт, США), экспертом по возобновляемым источникам энергии Национальной Лаборатории Лоуренс Беркли и Университета Северной Калифорнии (Беркли, США), экспертом Центра Стратегии республики Кабо Верде, экспертом Инновационного центра «Сколково», экспертом ОАО «НИИЭС» (Русгидро), членом редколлегии научно-технической редакции «АСТА Press» (Калгари, Канада), экспертом 7-й Рамочной Программы Евросоюза и т.д.

Личный вклад: выносимые на защиту результаты получены соискателем лично. В опубликованных совместных работах постановка и исследование задач осуществлялись совместными усилиями соавторов при непосредственном активном участии соискателя.

Публикации: по теме диссертации опубликовано 52 научных работы в центральных, региональных и международных журналах, в том числе 29, рекомендуемых ВАК Минобрнауки РФ. Выпущена 1 монография. Издано 2 учебных пособия. По теме исследований на 7 разработок получены патенты, 1 ноу-хау.

Структура и объем диссертации: диссертация состоит из введения, пяти глав, заключения и приложений, изложенных на 428 страницах машинописного текста, содержит 100 страниц приложений, 159 рисунков, 29 таблиц, список используемой литературы из 239 наименований.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во Введении обоснована актуальность работы, отмечена научная новизна и практическая ценность работы, приведены результаты внедрения.

В главе 1 «Состояние вопроса и задачи исследования» описан ряд технических особенностей основных конструкций ветроэнергетических установок, проведен анализ преимуществ и недостатков вертикально- и горизонтально-осевых ветроэнергетических установок, дан обзор ряда факторов, характеризующих ВЭУ, приведены пути повышения качества работы установки и снижения влияния на окружающую среду. Указаны основные критерии оптимизации. На основе выводов сформулированы цели и задачи исследования.

Исследования вертикально-осевых (ВО) ВЭУ получили интенсивное развитие с начала 80-х годов XX века. Сегодня практически все страны эксплуатируют вертикально-осевые ВЭУ с ротором Дарье. Диапазон мощностей ВО ВЭУ непрерывно расширяется и достиг 4 МВт. Однако интерес к малым ВЭУ также растет в связи с рядом причин (стремление потребителей к автономному энергоснабжению, доступность и простота монтажа, обслуживания, ремонта и т.д.). При этом стоимость малых ВЭУ, используемых для автономного энергоснабжения, неуклонно снижается за счет развития технологий производства комплектующих, снижения установочных расходов и т.д. Наблюдается тенденция приближения малых ВЭУ к местам проживания человека (стремление установить ВЭУ на кровле загородного коттеджа или многоэтажного дома, строении или сооружении). Такие решения весьма привлекательны как с технической, так и экономической точки зрения, и реализованы в США, Канаде, Европе, Индии, Китае. В России эта тенденция так-

же прослеживается, что подтверждается растущим объемом заявок от потенциальных потребителей на размещение ВЭУ до 30 кВт непосредственно на жилых домах. Кроме этого, в России растет активность агропромышленного сектора в части внедрения автономного энергоснабжения удаленных автономных объектов (энергопитание фермерских хозяйств, рыболовецких артелей, тепличных комплексов и т.д.). В связи с этим развитие малой ветроэнергетики, по мнению ряда российских и зарубежных экспертов, является одной из приоритетных задач ветроэнергетики и требует детального изучения, систематизации и совершенствования ВЭУ.

Традиционный процесс проектирования конструкций ВО ВЭУ вплоть до настоящего времени не имел четкого научного обоснования, представляя собой ряд обособленных теорий и методик с генерацией инженерных решений по каждому компоненту конструкции обособленно. В результате такого подхода почти все производимые в мире изделия имели ряд проблем, связанных с наличием шума, вибраций, дезинтегрированного функционирования и т.д. Таким образом, разработка единой методологии проектирования и создания ВО ВЭУ является крупной научной проблемой, имеющей стратегическое хозяйственное значение.

В данной диссертационной работе вертикально-осевая конструкция малой мощности была выбрана для изучения, разработки и оптимизации в связи с рядом факторов: независимость работы ВЭУ от направления ветра, высокий коэффициент использования энергии ветра (до 38%), самостоятельный старт ВЭУ на малых (до 2 м/с) скоростях ветра, выработка энергии со скорости ветра 3 м/с, возможность аэродинамического регулирования частоты вращения, низкий уровень шума и вибраций, отсутствие мультипликатора, простота изготовления компонентов, а также относительно низкая себестоимость (от 60 руб. за 1 Вт мощности).

Критериями оптимизации являлись технические условия (оптимальные характеристики ВЭУ, параметры и технологичность компонентов, уровень влияния автоколебаний и взаимовлияния модулей ВЭУ, самостоятельный старт при малых ветрах, ограничение частоты вращения ротора) и экономические параметры (низкая себестоимость, ускоренная пусконаладка).

В главе 2 «Методология разработки ВО ВЭУ» дано описание системного методологического подхода к разработке ВО ВЭУ на основе последовательнопоэтапной (итерационной) покомпонентной оптимизации конструкции с помощью применения последовательности ряда методов (рис. 1), исследования созданных математических трехмерных и функциональных моделей ротора и ВЭУ в целом. Вычислены внешние и внутренние возмущающие воздействия и их влияние на ВЭУ, обоснован выбор конструкций компонентов, описано аэродинамическое регулирование и оптимальное электронное регулирование мощности ВЭУ.

Процесс исследования на основе эксперимента является чрезвычайно дорогостоящим. С учетом мощного развития компьютерной техники и основанных на ней математических аппаратов, одним из быстрых и менее дорогостоящих путей изучения параметров ветроэнергетических установок является компьютерное моделирование и всесторонний анализ физико-математических моделей.

Методология разработки и создания ВО ВЭУ включает ряд теоретических и экспериментальных методов, применяемых на этапах разработки по итерационной (последовательно-поэтапной) схеме с оптимизацией каждого компонента и конструкции в целом, когда по завершении каждого этапа и ряда итераций проводится

сравнение расчетных и опытных характеристик с проверкой на удовлетворительность по соответствующим критериям.

Рис. 1. Структурная схема методологии разработки и создания ВО ВЭУ

Метод разработки и оптимизации профиля лопасти основан на патентном поиске, теоретических и экспериментальных исследованиях различных профилей лопастей, анализе результатов и выборе наиболее эффективного профиля. Критерий оптимальности профиля - обеспечение максимума интегрированной кинетической энергии, передаваемой профилю потоком ветра силовыми импульсами, при отсутствии автоколебаний профиля во всем диапазоне углов атаки.

Расчет профиля производился в программном комплексе АтуБ СРХ. который основан на конечно-объемном методе (МКО) решения уравнений гидродинамики, таких как уравнение неразрывности, уравнение сохранения энергии и уравнение количества движения. Основная идея МКО заключается в разбиении расчетной области на N-6 число непересекающихся контрольных объемов таким образом, что каждая узловая точка содержится в одном контрольном объеме. Дифференциальное уравнение интегрируется по каждому контрольному объему. Для вычисления интегралов используются кусочно-непрерывные функции, которые описывают из-

менение зависимой переменной (одной из составляющих скорости) между сеточными узлами. В результате находится дискретный аналог дифференциального уравнения. Дискретные уравнения вычисляются с помощью метода Algebraic Coupled Multigrid (AMG), разработанного M.Raw и G.Schneider. Данный метод использует неявную связанную схему решения системы линейных алгебраических уравнений. Вычислительные затраты метода линейно зависят от числа узловых точек. Одним из важных свойств МКО является то, что в нем заложено точное интегральное сохранение таких величин, как масса, количество движения и энергия на любой группе контрольных объемов, а следовательно, и на всей расчетной области. Это свойство проявляется при любом числе узловых точек. Таким образом, даже решение на грубой сетке удовлетворяет точным интегральным балансам. Ansys CFX позволяет проводить расчеты на смешанных сетках, состоящих из различных типов элементов: тетраэдров, призм, клиновидных элементов и гексаэдров. При расчете стационарных вариантов процесс итерации по времени завершается при достижении уровня сходимости.

При моделировании использовалась модель вязкого газа с осреднением турбулентных характеристик (осредненные по Рейнольдсу уравнения Навье-Стокса несжимаемой жидкости):

^V[_q dvi | d(VjVj) _ 1 дР д

()xj . dt dXj p (}Xj r?Xj

где Ху, Ц=1,2 - декартовы координаты (х,у) частицы; I - время; V,- - декартовы составляющие вектора средней скорости потока (у,а); Р - давление; р - плотность; и и и, - молекулярный и турбулентный коэффициенты кинематической вязкости. В расчетах учитывались нестационарные явления (динамический срыв потока, образование сложной системы вихрей, повышение уровня турбулентности в затененной области, взаимодействие вихрей различных размеров, скоростей движения и интенсивности с твердыми поверхностями ротора ВО ВЭУ).

Модель исходного профиля, подлежащего оптимизации, построена в программе Ашуз СБХ и приведена на рис. 2а. Расчетные значения переменных - коэффициентов подъемной силы Су и силы сопротивления Сх для разных углов атаки с поворотом лопасти к моделируемому набегающему потоку ветра со скоростью в диапазоне 2-12 м/с и шагом угла атаки в 1° в диапазоне 0°- 360° получены в процессе расчета. Формулы коэффициентов сил выражаются уравнениями (2).

Су = -4-; Сх = -4- (2)

■> Р* .с .С ’

2 2 ’

где У- подъемная сила, X- сила лобового сопротивления, р - плотность среды, V - скорость потока, 5 - площадь крыла (профиля). Сила лобового сопротивления складывается из профильного Хр (сопротивление давления и трения), индуктивного X; и волнового сопротивлений X,,.:

Х = Хр+Х;+Ху = Хр. (3)

На рис. 2 приведен профиль лопасти и схема сил, действующих на него в

стационарном потоке газа, где V - скорость потока; ХЛ - положение центра давле-

ния; а - угол атаки; Хр; У - сила сопротивления профиля и подъемная сила в поточной системе координат; Мг - момент сил относительно носка профиля; Т, N -тангенциальная и нормальная силы в связанной системе координат; Ь - длина хор-

ды; ц=рУ2/2 -скоростной напор; С^Х^/Ь и Су=У/(ц-8) - коэффициенты сил давления и подъемной; Я.=Ь/Ь - удлинение лопасти; Ь - длина лопасти; 8=Ь-Ь - площадь лопасти; СП1=Мг/^-8-Ь) - коэффициент момента; Сп=1Ч/(ц-8); С,=Т/(ц-8) - коэффициенты нормальной (поперечной) и тангенциальной (продольной) сил; Сх=Схр+ Су-2/(л:-Х) - коэффициент сопротивления лопасти; с - максимальная толщина профиля; хс — положение максимальной толщины от носка профиля.

Рис. 2. Профиль лопасти (а) и схема сил, действующих на профиль в стационарном потоке газа (б)

Оптимизация профиля производится с учетом граничных условий 0,5>с/Ь>0,1 и 0,5>хс/Ь>0,1. После оптимизации по критерию

(2i=6o°°Cyi = MAX (4)

lz?=o° Cxi = MIN

на промежуточных этапах оптимизации было получено несколько профилей, характеристики которых приведены в Приложении к диссертационной работе.

В результате программной оптимизации получен профиль SRCV-2035 с относительной толщиной с/Ь=0,20 (20%) и с хс/Ь=0,35 (35%), относительные координаты которого приведены в табл. 1.

Таблица 1

Координаты профиля 8КСУ2035 с относительной толщиной с/Ь=20% и расположением максимальной толщины на 35% от носка профиля согласно рис. 2а.

х/Ъ 0 0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009

у/с 0 0.03507 0.04960 0.06075 0.07015 0.07843 0.08593 0.09280 0.09921 0.10523

х/Ь 0.01 0.02 0.03 0.04 0.05 0.05192 0.06 0.07 0.08 0.09

_ У/с 0.11092 0.15687 0.19212 0.22185 0.24803 0 25277 0.27117 0.29140 0.30956 0.32633

х/Ь 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55

у/с 0.34223 0.40872 0.454353 0.48235 0.49625 0.5 0.49524 0.48091 0.45740 0.42709

х/Ь 0.6 0.65 0.7 0.75 0.8 0 85 0.9 0.95 0.98681 1

у/с 0.39211 0.35308 0.31076 0.26637 0.22063 0.17386 0.12721 0.08124 0.04529 0

После проведения теоретических вычислений профилей изготовлены модели лопастей, имеющие хорду длиной Ь=0,12 м и относительное удлинение >.=5 (рис. 36). Модели исследованы на предмет получения динамических характеристик профилей в гидродинамической трубе 44К200 Государственного Ракетного Центра (КБ им. В.П.Макеева, Челябинская обл.), г. Миасс (рис. За). Число Рейнольдса в экспериментах составило 630000, что обеспечило режим автомодельности по этому параметру. Скорость потока воды в экспериментах составила 5 м/с, угол атаки менялся с шагом в 1° в диапазоне 0° - 360°. Конструкция экспериментально-

измерительной установки обеспечивает получение сразу профильного сопротивления за счет исключения влияния концов лопасти, дающих вклад в индуктивное сопротивление, с помощью отсекателей. Калибровка проведена на основе сравнения справочных и вновь полученных экспериментальных данных по профилю ЫАСА0015.

а) б)

Рис. 3. Большая скоростная гидродинамическая труба с горизонтальным рабочим участком 44к200 (а) и экспериментальные модели профилей (б)

Испытания показали, что для профилей с острой задней кромкой при углах атаки а, близких к 180° (152°-172°) и (192°-208°) наблюдались сильные вибрации профиля, обусловленные периодическим процессом перехода от безотрывного к отрывному течению и обратно. Этот процесс вызывается обтеканием острого хвоста данного профиля в указанном диапазоне углов атаки. Измеренные при этих углах атаки параметры могут быть интерпретированы как средние величины. В то же время возможно и искажение измеряемых параметров. Так при 152° значение коэффициента центра давления выходит за пределы хорды профиля. У профиля 8ЯСУ2035, имеющего утолщенный скругленный хвост, никаких подобных колебательных процессов не наблюдалось.

Рис. 4. Сравнение коэффициентов расчетной подъемной силы Сус, силы сопротивления Схс, условного центра давления Са с соответствующими экспериментальными значениями Суе, Схе, Сае в поточной системе координат для 8ЯСУ2035

Результаты расчетов и экспериментальные данные приведены на рис. 4. В связи с проведенным анализом, а также на основании патентных исследований был выбран оптимальный профиль 8КСУ2035, на который получен патент на изобретение №2347104 (2006 г.). Полученные экспериментальные данные с учетом результатов численных расчетов гидродинамических характеристик профилей были использованы для расчета аэродинамических характеристик роторов и в математических моделях исследуемых ВЭУ.

Данный метод расчета параметров профиля имеет очевидное преимущество по сравнению с традиционным методом, т.е. построением поляр каждого профиля с предварительной продувкой в аэродинамической трубе и является более быстрым, удобным и менее дорогостоящим.

Разработка и оптимизация ротора ВО ВЭУ основана на последовательности ряда методов и включает в себя патентный поиск, трехмерное компьютерное моделирование компонентов ротора, функциональное математическое моделирование ротора, исследование взаимовлияния компонентов ротора, оптимизацию параметров ротора на основе итерационного метода, разработку агрегатов аэродинамического регулирования и т.д. Критерии оптимальности — самостоятельный старт ротора во всем диапазоне быстроходности, ограничение максимальной частоты вращения, достижение максимального коэффициента мощности (коэффициента использования энергии ветра, КИЭВ), снижение влияния концевых эффектов на лопастях, эргономика. Патентный поиск проведен в Российских и зарубежных базах данных согласно ГОСТ Р 15.011-96. На конструкцию ротора получен патент № 2347104 (2006 г.). Конструкция ротора выбрана двухъярусной 6-лопастной на основании опыта разработчиков и патентных исследований. Второй ярус введен для улучшения конструкции ротора и дополнительного увеличения мощности.

ВЭУ представлена рядом основных модулей (рис. 5) и является интегрированной системой превращения кинетической энергии набегающего потока ветра в электрическую энергию за счет непосредственной передачи момента вращения ротора на вал генератора.

Ротор ВЭУ с лопастями

Г енератор

Ступица

Мачта

Рис. 5. Конструкция малой ветроэнергетической установки ВЭУ-3 с вертикальной осью вращения мощностью 3 кВт

Аэродинамический регулятор

Анализ литературы и предварительных расчетов показывает, что при коэффициенте быстроходности Z < 2 динамический срыв потока с лопастей ротора Дарье доминирует на большей части траектории лопасти. Увеличение числа модулей быстроходности приводит к подавлению процесса формирования динамического срыва потока с вертикальных лопастей, но также снижается возможность самостоятельного старта. Тогда с учетом графика идеального Дарье граничные условия для Ъ\ 6>г>2. В зонах умеренных ветров ВЭУ должна начинать работу с 4 м/с. Ротор получает 100% энергии с периметра. Расчет производится для радиуса II, который вычисляется из соотношения 1,5-11 > Ь > 4Ь, где Ь - длина лопасти. Такие данные взяты из соображений прочности компонентов и минимизации влияния индуктивного сопротивления. Оптимальная быстроходность исследуемого ветро-колеса при наибольшем коэффициенте мощности Ср составляет Z=4■л/m=2, где т

- количество лопастей (6 в исследуемом варианте). Или, с соответствующей поправкой для тщательно спрофилированных лопастей будет на 34% выше, т.е. 2,68. При скорости ветра 4 м/с окружная скорость лопастей \=10,7 м/с. Ветер взаимодействует с лопастями с окружной скоростью. Для устойчивой работы ВЭУ необходимо число Рейнольдса Яе>300000. Тогда хорда лопасти должна быть Ь=11е/(68500-у)=0,4 м. Граничные значения при оптимизации хорды лопасти каждого яруса составили 0,5м>Ь>0,Зм. Приняв длину лопасти Ь в соответствии с вышеприведенными доводами за 2 м, радиус Я получается 1,7 м (с учетом добавки радиуса ступицы 0,4 м). Тогда соотношение хорды Ь к диаметру О составит 11%, что совпадает с соотношениями горизонтально-осевых ВЭУ (5-15% по среднему диаметру лопасти). При этом коэффициент геометрического заполнения ротора получается 0,35. Длина лопасти 2 м выбрана из соотношения Ьлоп - 6Ь для исключения влияния концевых эффектов, а также из соображений эргономики.

Таким образом, получены основные параметры ВЭУ:

Б = 3,4 м - диаметр ротора; 2ЬЛ0П = Н = 4,1 м - высота ротора; 8 = Э Н = 13,94 м2- ометаемая площадь ветроколеса; Ь = 0,4 м - длина хорды симметричного профиля лопасти; с = 35% - относительная толщина профиля; т = 6 - количество вертикальных лопастей.

Расчет главных аэродинамических характеристик ротора выполнен по аналогии с импульсным методом, использующим два активных диска. В этом методе индуктивные скорости, учитывающие влияние на воздушный поток вращающегося ветроколеса, задаются как с наветренной, так и с подветренной стороны ротора. Отличие проведенного расчета заключается в том, что в импульсном методе используют закон сохранения импульса и на основе этого получаются системы уравнений и их решения. В предлагаемом методе с использованием системы параметров индуктивные скорости задаются так, чтобы максимально учесть влияние быстроходности и коэффициента заполнения, равного отношению суммарной площади лопастей в плане к фронтальной площади ветроколеса, а также частично учесть интерференцию ярусов ветроколеса для многоярусных компоновок.

На рис. 6 показано, как при некоторых значениях азимутального угла \|/ и соответствующего положения вектора возмущенной скорости ветра Уп по отношению к лопасти возникает крутящий момент. Вектор скорости \№, относительной по отношению к лопасти, складывается из вектора возмущенной скорости ветра Уп и

вектора -соЯ, противоположного вектору линейной окружной скорости лопасти. Здесь Уп = V + ДУх- і + ДУу- j , где і,- орты системы координат, ДУх, ДУу -компоненты размерной индуктивной скорости, со — угловая скорость вращения, V

- вектор скорости ветра, Я - радиус ветроколеса. Вектор относительной скорости образует угол р с касательной к окружности с центром на оси вращения и прямой, проведенной через точку условного аэродинамического центра профиля лопасти, расположенную на хорде профиля в 35% от его носка. В эту же точку проведен также радиус ветроколеса Я. Если угол атаки лопасти а = р - ср, где ср -установочный угол лопасти, является небольшим, то в этом случае сила сопротивления лопасти X мала, а подъемная сила У, наоборот, большая, и вектор результирующей аэродинамической силы Яа направлен так, что создает положительный момент в сторону вращения ветроколеса. Однако такая картина наблюдается не на всей круговой траектории лопасти.

Рис. 6. Система координат и схема образования крутящего момента

Относительные безразмерные компоненты индуктивной скорости определяются выражениями: и = ДУх/У и V = ДУу/У, где V - модуль скорости ветра. Базовые компоненты индуктивной скорости вычисляются по формулам:

111 = к1-282/к2;

V] = Ц1-(1+5у), 5у = -0.3;

у2 = -и1(1-бу); (5)

и2 = шах (1,1 - (НиОО+г-и!)),

Ъ — ш11/У - быстроходность;

к1= асг'(а/асг)м, если а>асг, иначе к]= ссг ( а/<тсг)м;

к! - параметр, 81 = 0,9 ; Б2 = 0,4 - параметры;

а - 1Ь/(2-К) - коэффициент заполнения; (6)

Ь - хорда лопасти, 1 - количество лопастей в ярусе; осг = (!/2)я-аП1 - переходный коэффициент заполнения;

0Ш = 0,28; Б = 0,57 - параметры; вг = В-(А/С)а;

А = 2,97; В = 0,71; С = 2,8; с!= 0,85 - параметры.

Расчет составляющих индуктивной скорости в зависимости от углового положения лопасти \|/ проводится с помощью следующих соотношений:

При 0 <ц>< л:

и = 111; V = -и1-(со8(\|/) +5у);

При Я < V)/ < 7С +

и = -и! + (и! - и2)(зт(тг((у - ге)/ёрз1 - 1/2)) + 1)/2; у = -у2 + (у„-у2) (5т(я-(('|/- тг)Мр51 - 1/2))+ 1)/2, где у0 = 0 - параметр.

При к + с1рз! < у < 2л- с1рзк

и = -и2; у = у0; _ (7)

При 2л - (1рз1 < у < 2л: '

и = -и2+(и2-и1)-(5ш(7г((у-2л)/с1рз1+ 1)- 1/2))+ 1)/2;

У = Уо+^1-Уо)-(5Ш(71-((\1/-2я)Мр81+1)- 1/2)+ 1)/2, где ёрэ! = л/9 - угол перехода с наветренной к подветренной стороне.

Составляющие относительной скорости набегающего под углом атаки а = Р - ф на лопасть потока определяются формулами:

\ух = 1 + и + г-соз(\|/); \уу = у + г-зт(у);

™х = Ых/У;\уу = Ыу/У. (8)

Вектор \У набегающего на лопасть потока изображен на рис. 6. Графики индуктивных скоростей для разных значений быстроходности приведены на рис. 7.

Компоненты индуктивной скорости при 2 в 1.3 Компоненты индуктивной скорости при 2 = 2.3

Угол поворота лопасти в градусах Угол поворота лопасти 8 градусах

Рис. 7. Компоненты индуктивной скорости для разных Z

Коэффициент силы, действующей на лопасть в направлении оси вращения параллельно радиусу, проведенному в аэродинамический центр, и коэффициент крутящего момента находятся по формулам:

Сг = (Сх(а)этр + Су(а)-со5Р)-(\¥х2 + \уу2)-Ь/(2-11);

С,, = (Су(а)зіпР - Сх(а)-созР)'(\\'х2 + луу2)-Ь/(2-К); (9)

Сч = Сц1 • (1 -(Са-Сао) Ь/(К-эикр)) + Сг'(Са-Сао)'Ь/(К-соз(р),

где:

Сг = Ргасі / (р-У2-11-На) - коэффициент радиальной силы Рга(і;

Сч = Мг / (р-У2'Я2-Н5,) - коэффициент крутящего момента Мг;

р - плотность воздуха;

Н„ - высота яруса ротора (длина лопасти в ярусе);

Счі - значение коэффициента крутящего момента при

совпадении центра давления с аэродинамическим центром;

Сі = Хі / Ь - коэффициент центра давления лопасти;

Ха - положение центра давления от носка профиля;

Сда = 0,35 - относительная координата условного аэродинамического центра профиля лопасти;

Сх(а), Су(а) - коэффициенты силы сопротивления и подъемной силы лопасти в скоростной системе координат, причем:

Сх(а) = Срх(а) + Су2(а)/(зг->.),

где:

Срх(а) - профильный коэффициент сопротивления;

л - удлинение лопасти.

В силу того, что центр давления меняет положение по мере движения лопасти по окружности, в формулах (9) вводится соответствующая поправка на коэффициент крутящего момента.

В расчетной программе, написанной на языке С++ и реализующей описываемый метод, определяются также составляющие сил по осям Ох и Оу (рис. 6) для каждого элемента ветроколеса и вычисляются интегральные силы, действующие на ротор в целом в зависимости от угла его поворота у и быстроходности 2. Для отдельной лопасти составляющие коэффициентов сил будут:

Сх = -СцСОвГу) + Сг8Іп(ц/);

Су = -Сч*5ш(\|/) - Сг-соз(\|/). (10)

Расчет главных аэродинамических характеристик ротора выполнен по следующему алгоритму: найденные коэффициенты С* и Су в каждый дискретный момент (в каждом положении ротора) для каждой скорости ветра при заданной быстроходности Ъ однозначно определяют тангенциальную и радиальную составляющие сил, действующих на ротор (рис. 26). Из них определяющий вклад в значение крутящего момента М* дает тангенциальная составляющая. Интегрируя М2 всех лопастей за полный оборот от 0° до 360° находится значение Сш при заданной

быстроходности. Далее для значений быстроходности Ъ от 0 до 5 вычисляются графики Ст(г) и Ср(2):

Ст = 2'М2/(р-У"'8-К)- коэффициент крутящего момента;

Ср = 2Р/(р-У3 8) = Ст 2 - коэффициент использования энергии ветра, где:

в = 2-Я-Н - фронтальная площадь ветроколеса;

Н - высота ветроколеса;

Р - аэродинамическая мощность ветроколеса;

Мг - крутящий момент.

На рис. 8 представлена главная аэродинамическая характеристика Ср(2) и зависимость коэффициента момента Сш(г) ротора с шестью лопастями.

Эффективность Ср(2) и коэффициент Спі(2)

Рис. 8. Коэффициент использования энергии ветра Ср и коэффициента крутящего момента Ст от быстроходности 2

В процессе расчета аэродинамических характеристик ротора была проведена программная оптимизация установочного угла лопасти с дискретизацией 0.1° при граничных условиях от 0° до 10°. В результате максимальное значение коэффициента использования энергии ветра без учета потерь на трение и люфты было получено при установочном угле лопасти 4,1°.

Можно утверждать, что задача оптимизации параметров ротора ветроэнергетической установки номинальной мощностью 3 кВт выполнена. Малые модули 2, т.е. тихоходность ротора при достаточном значении эффективности не является недостатком. Снижаются нагрузки на ротор и его элементы, снижается уровень вибраций и шума, что положительно сказывается на работе ВЭУ в целом.

На основе главной аэродинамической характеристики можно построить поле мощности, т.е. изолинии аэродинамической мощности в координатах «скорость

прочностных расчетах и при разработке системы регулирования мощности.

Рис. 9. Поле аэродинамической мощности ротора ВЭУ (Вт)

Таблица 2

Экспериментальные данные мощности ВЭУ РЭКСПер (в Ваттах, по вертикали частота вращения ротора, об/мин, по горизонтали скорость ветра, м/с)

1 V, мЗё Т5 1 1 4 6 § Г 8 Тб 11

6 б 0 0 е 0 0 0 0 б 0 б

30 0 0 '6.2 1бъ 6 218 4 261 8 2. її 2 30? 9 346 3 391.8 444 3

35 0 0 Шз 272 § Ш'4 38? У Т99ЇЇ 438.о " 486.5 <42,6

40 0 0 тл 3232 409 3 469 4 611 7 552 2 601.5 660.2

45 6 и 9.9 Ї 73.4 Зб2 0 486 8 о 77 4 641 7 690 5 і 40 4 600,5

40 00 133.2 312.3 558.5 685.8 ?803 849.3 3067 967 2

ь, 55 ’б.б 35 & €Г? і Ш.5 іТзТЗ ШТ?

Я 60 об 327 3 Ш.З 886,9 1063.3 1201 5 1306.3 1389 0

» ?0 0.0 536 4 96 § 0 1197.8 13814 1625.0 1636 2

<5 80 0 0 509 6 1015 3 -333 9 155? 4 1747.1 1896.8

вГ 30 0 0 ш^т-~ Ї7ЖГІ 1М.Ї 2163 ^

В юо 0 0 466,4 896 7 1506 1 1884 9 2133.6 24328

Г? 110 С () 547.8 тьз 2350 3 2693 ?

(§5 120 0 0 859 9 1508 4 2127 5 2535 6 2957.3

130 об 14Й8.2 іті: 27*368 3265 8

140 0 0 1387.6 2163 4 28§6.1 34 4 І 5

150 об 21260 2339 1 3655.2

160 об " 382/3"

1/0 0 0 2346.1 3934.8

180 0 0 2Й54 9

[Коэффициент і!СП5ПЬ50Ва«ИЯ энергии зетрэ Ср = Рзісспер-'Ро = Рзкспвр їр ЗУ V У/2)

і 6 36104 | 0 36500 1 0.36499 1 0.36425 С36531 ( 0.364# | 0 35500 І 0.3653*

, 4 Г~ Я 6 і ? 8 9 10 11

133 8 1 3721 іШІШі Ш9 Б- 1528.2 2171,2 2978,4 Ш 9

После проведения теоретических вычислений изготовлен ряд роторов ВЭУ и проведены испытания на полигонах предприятия «ГРЦ-Вертикаль» в Челябин-

ской области. Результаты испытаний всех роторов различаются не более, чем на 3-7%. Данные выходной мощности ВЭУ-З(б). №03-09-008, полученные с помощью записи на флэш-карту на оконечном выходе электронного регулятора, приведены в табл. 2 с пересчетом и приведением к аэродинамической мощности с учетом КИЭВ. Экспериментальный график мощности приведен на рис. 10.

4500,0 .................. ...............

Мощность Р. Вт

О г 4 6 8 10 12

Рис. 10. График мощности ВЭУ по экспериментальным данным

Экспериментальные измерения проведены на скорости ветра от 3 до 11 м/с. До 3 м/с выработка энергии минимальна и не представляет интерес. На скорости ветра 11 м/с активируются аэродинамические регуляторы и частота вращения стабилизируется на 180 об/мин во всем диапазоне ветров. Мощность ВЭУ нормируется на 11 м/с и составляет минимум 3 кВт с учетом потерь на преобразование.

Анализ данных говорит о том, что реальный КИЭВ составляет 0,365 (в сравнении с теоретическим 0,38-0,40). Этот факт объясняется наличием неучтенных при расчетах потерь (трение, аэродинамические потери в местах болтовых соединений, погрешности сборки, люфт и т.д.). Это позволяет утверждать, что можно получить больший КИЭВ с учетом ряда конструктивных доработок. Однако очевидно, что задача регулирования решена, реализована и отработана. Непосредственно электронный регулятор рассматривается ниже. На основании приведенных данных можно утверждать, что параметры ротора оптимизированы.

В процессе разработки ротора ВЭУ создана рабочая конструкторская документация по ГОСТ ЕСКД. эволюционирующая по мере совершенствования конструкции и ее оптимизации. Проведены поверочные расчеты конструкции с учетом внешних возмущающих воздействий, взаимовлияния компонентов друг на друга и на окружающую среду.

Данный метод расчета характеристик ротора имеет ряд преимуществ по сравнению с инженерными и эмпирические методиками (импульсным и вихревым моделированием), методом моделирования невязкого газа (метод дискретных вихрей, уравнения потенциала и Эйлера) и методом моделирования вязкого газа с учетом нестационарных турбулентных характеристик, и является наиболее применимым на сегодняшнем уровне развития науки, программных средств и процессоров.

Поэтапная оптимизация конструкций механических компонентов ВЭУ на основе метода итераций (повторяемых поэтапных прогонов) двух последовательных математических моделей, трехмерной и функциональной, является менее затрат-

ным способом улучшения изделия по сравнению с экспериментом. Процесс оптимизации параметров конструкции ротора распадается на ряд задач, заключающихся в последовательном исследовании свойств его компонентов (рис. 11).

Разработка и анализ конструкции компонентов ротора производится с точки зрения нанлучшей аэродинамики с учетом требований по безопасности, устойчивости к внешним воздействиям, надежности, эргономики и технической эстетики, с учетом ГОСТ Р 51237-98, ГОСТ Р 51990-2002, ГОСТ Р 51991-2002. Исследования проведены с рядом допущений (граничных условий), не оказывающих существенного влияния на конечный результат (ламинарность потока, устойчивость давления и направления ветра, однородность материалов и т.д.).

Рис. 11. Поэтапная итерационная оптимизация параметров компонентов ВЭУ

Характеристики внешних воздействий на конструкцию ВЭУ определены с учетом технических условий безопасности ветроустановок (Требования Ллойда, стандарты РФ). По результатам проведенного анализа для расчета нагрузок и напряжений в конструкции при создании модели выбран ряд расчетных случаев, максимальный из которых соответствует скорости ветра 45 м/с при континентальных условиях окружающей среды. Коэффициенты безопасности заложены: от гравитационных сил - 1,2; от инерционных сил - 1,3; от аэродинамических сил - 1,5.

С целью вычисления влияния внешнего механического воздействия на внутренние силовые факторы и напряжения в сечениях компонентов ВЭУ расчет проведен сначала отдельно по воздействиям, затем по группам воздействий, и далее по совокупности всех воздействий. В группы объединены воздействия статического характера (статические нагрузки) и воздействия динамического характера (динамические нагрузки). Статическому анализу подвергались отдельные элементы конструкции, прочность которых является определяющей для ВЭУ. Цель статического анализа заключалась в определении прочностных характеристик элементов

конструкции, которые служили основой для динамического анализа. Динамический анализ проведен для конструкции ВЭУ в целом с целью определения запаса прочности при предельных условиях эксплуатации ВЭУ. Совокупность статического и динамического анализа составляет одну итерацию в процессе оптимизации параметров элементов конструкции. Конечная цель оптимизации - минимизация массо-габаритных характеристик и стоимости изделия.

Статический анализ проводится с учетом силового воздействия, учитывающего гравитационные, инерционные и аэродинамические воздействия.

Гравитационные силы представлены произведением ускорения свободного падения на распределение масс, принятое в расчетной схеме конструкции.

Моделирование центробежных сил обеспечивалось заданием в качестве исходных данных частоты вращения ротора. Нагрузки, связанные с дисбалансом ротора, определены из условия, что он имеет статическую несбалансированность, обусловленную поперечным смещением центра масс относительно оси вращения на величину с=0,5% от радиуса ротора. Указанная величина эксцентриситета задана в соответствии с Требованиями Ллойда для ротора, изготавливаемого с выполнением операции предэксплуатационной балансировки.

Аэродинамические силы заданы согласно результатам аэродинамического расчета вращающегося в стационарном ветровом потоке ротора ВЭУ. Указанный расчет включал определение аэродинамических внешних нагрузок (сил и моментов сил), действующих на ротор в целом, а также отдельно на лопасть.

Аэродинамические внешние нагрузки на ротор задавались компонентами Рх=СхЧ8; РУ=СУ-Ч'5; (И)

где Рх, Ру - проекции вектора результирующей аэродинамической силы на оси X, и поточной системы координат ХуУ,^. (ось Ху направлена по вектору горизонтальной скорости ветра, ось Уу нормальна к нему и к оси вращения ротора, ось 7,у совпадает с осью вращения ротора); М7 - крутящий момент вокруг оси ротора; Сх, Су , Ст - аэродинамические коэффициенты, зависящие от углового положения ротора по отношению к вектору скорости ветра (угла ф между плоскостью и плоскостью, проходящей через ось Ъх и продольную ось одной из лопастей), а также от параметра быстроходности 2=(оК/У; q= рУ2/2 - скоростной напор ветра; р=1,225 кг/м3 - плотность воздуха; V - горизонтальная скорость ветрового потока; 8=2-11-Ь - площадь прямоугольника, являющегося проекцией ротора на плоскость У%7^ поверхности, ометаемой ротором; Я =1,7 м - радиус ротора; Ь =4,08 м - высота ротора с учетом толщины среднего кольца; а - угловая скорость вращения.

Вычисленные значения сил Рх, Ру и момента М7 являются периодическими. За один полный оборот ротора реализуются шесть периодов, причем переменная составляющая нагрузок близка к гармонической. Исходя из этого, в расчетной модели было принято, что каждый из действующих на ротор силовых факторов Рх, Еу, Мг представляет собой сумму некоторого среднего за оборот значения (называемого далее статической составляющей) и переменного слагаемого (динамической составляющей). Динамическая составляющая является гармонической, а ее частота равна увеличенной в шесть раз частоте вращения ротора.

В статическом анализе применялся метод конечных элементов (МКЭ), с помощью которого была произведена аппроксимация непрерывных величин, определяемых бесконечным числом значений на рассматриваемой области, т.е. дискретной моделью, которая строится на множестве кусочно-непрерывных функций, определенных на конечном числе подобластей.

При построении учитывалось допущение, что находящееся в равновесии упругое тело под действием внешних нагрузок имеет полную потенциальную энергию:

П(у) = и(у)-\У(у) , (12)

где П - полная потенциальная энергия, и - потенциальная энергия деформации, V/ - потенциал внешних нагрузок, V - произвольное поле возможных перемещений, удовлетворяющее граничным условиям. Из принципа возможной работы следует, что в состоянии равновесия полная потенциальная энергия системы минимальна. Для нахождения действительного поля перемещений необходимо минимизировать приведенное выражение на множестве всех функций V, удовлетворяющих граничным условиям. Тогда функция, представляющая минимум, будет являться искомым полем перемещений \у. С другой стороны, нахождение минимума П(у) заключается в решении дифференциальных уравнений теории упругости, являющихся бесконечномерными задачами. В комплексе М8С 1ЧА8ТКАЫ замена бесконечномерной задачи п-мерной осуществляется с помощью последовательности построения дискретной модели:

- В рассматриваемой области упругого тела фиксируется конечное число точек (узлов). Функция в узле определяется одним значением (скалярная или одномерная векторная функция) и считается переменной, подлежащей определению.

- Область разбивается на конечное число элементарных подобластей, имеющих общие узлы. Подобласти в своей совокупности аппроксимируют форму упругого тела за счет аппроксимации на каждом элементе полиномиальными функциями формы, значения которой внутри элемента и на его границах определяются через значения функции в узлах.

Приведенный алгоритм МКЭ использует метод Релея-Ритца, основанный на введении п пробных функций вместо бесконечномерной задачи, с нахождением частной комбинации действительного поля перемещений V/, которая минимизирует П(у). Вводится ряд неизвестных коэффициентов и;, являющихся узловыми значениями функции Л, заменяющей V,' при минимизации упомянутого выражения по всем возможным функциям п по подпространству V. Тогда приближенное решение для одномерной задачи будет выглядеть в виде кусочно-линейных функций и их комбинаций, сшиваемых в 1-х узлах:

П = (13)

1 = 1

Энергия деформации вычисляется интегрированием по объему произведения функций деформации и напряжения, для чего определяются соответствующие деформации и напряжения. Для перевода перемещений в деформации принят оператор дифференцирования Ь, а для напряжений - оператор 8, определяющийся добавлением к оператору Ь соответствующей матрицы Гука.

a-s|^ujVj|

Тогда выражение для энергии деформации:

(14)

(15)

где интеграл определен на всей конечной области.

В результате расчета определяются эквивалентные напряжения в каждом элементе. При задании свойств трехмерного ортотропного материала использовались следующие соотношения зависимости деформации от напряжений и температуры:

ell

£22 £33 у 12 у23 у13

1

ёТ

-у12

Е1

-у13

Е1

-у21

Е2

Ё2

-VZ3

Е2

-у31 ЕЗ -у32 ' ЕЗ 1

ЁЗ

1

G12

1

С23

1

С13 J

all al'

о22 u2

оЗЗ . +(T-Tief) ■ . a3

о12 0

а23 0

_о13 0

(16)

где Е, - модули упругости, Gij - модули сдвига, v,, - коэффициенты Пуассона, бу -деформации в соответствующих направлениях. Индексы i и j обозначают номера ортогональных осей материала. Изменения температуры Т во внимание не принимались.

Пример ниже приведен для расчета напряженно-деформированного состояния лопасти. Исходными данными для расчета напряжений в программном комплексе MSC NASTRAN являлась ее трехмерная модель, построенная на основе эскизной конструкторской документации. Для моделирования обшивки использовалась модель ортотропного материала и плоские оболочечные элементы типа «PLATE» (пластина). Конечно-элементная модель лопасти приведена на рис. 12. Закрепления модели в пространстве: один конец лопасти зафиксирован, другой конец имеет шарнирное закрепление. Приложенная нагрузка характеризуется усилием от центростремительного ускорения, радиальной аэродинамической силы и динамической составляющей аэродинамических сил. Основной нагрузкой на лопасть является центростремительное ускорение как постоянная нагрузка. Дополнительно на лопасть действуют переменные аэродинамические нагрузки.

Центростремительное ускорение равно:

,2

'ц

-со*- R

где

R - радиус расположения лопасти; Угловая скорость равна:

к ■ II

со-

30

где п=300 оборотов в минуту (максимальное расчетное значение).

(17)

Соответственно, усилие, определяемое центростремительным ускорением, распределено по поверхности модели.

Радиальная аэродинамическая сила равна:

р-У2

рл=сг.£-— гш

2 , (19)

где Сг ~ 0,54-аэродинамический коэффициент при Т = “~ут~; 2 - параметр

быстроходности; р- плотность воздуха; V - скорость ветра; Ь - длина лопасти.

С учетом динамического характера силы для расчетных оценок принята величина, составляющая 30% от инерционной силы. После задания граничных условий выполнен программный расчет. Результаты расчета деформаций лопасти приведены на рис. 126. Оценка прочности лопасти проводится по расчетному коэффициенту запаса прочности, показывающему уровень превышения допускаемого напряжения материала над максимальными действующими эквивалентными напряжениями. Запас прочности равен:

ц = -И- > 1,0

<?эи, . (20)

Минимальный коэффициент запаса прочности равен:

7 = — = 2,0. (21)

а экв

где [о]=280 МПа - допускаемые напряжения для стеклоткани Т-Ю(ЮО) ГОСТ 19170-73; 0экв=14О МПа - максимальные эквивалентные напряжения.

Вывод: на основе трехмерной (ЗБ) компьютерной модели, построенной в программном комплексе М8С NASTRAN, определены силовые факторы и напряженно-деформированное состояние, вычислены частоты собственных колебаний компонентов, проведен подбор удовлетворительных материалов, конструкция проанализирована на удовлетворительность всех параметров. Вычисления приведены для предельно допустимой скорости ветра 45 м/с.

а) б)

Рис. 12. Конечно-элементные модели ступицы ротора (а) и лопасти (б) ВЭУ-3 в МАБТКАЫ (картина распределения эквивалентных напряжений)

Максимальные напряжения в основных элементах силовой конструкции ВЭУ сведены в таблицы. Наиболее нагруженным элементом конструкция является основание ступицы (рис. 12а), в котором максимальные напряжения достигают 1250 МПа (в зонах крепления траверс).

Динамический анализ конструкции выполнен на основе расчетной модели силовой конструкции ВЭУ-З(б) также с учетом применения программновычислительного комплекса М8С КА8Т11АМ, реализующего метод конечных элементов. Число степеней свободы в модели составляло примерно 20000. Общий вид расчетной модели показан на рис. 13. Анализ полученных результатов позволяет считать возможным обеспечение прочности силовой конструкции при принятой конструктивной схеме ВЭУ-3 (6). Результаты расчетов приведены в Приложении к диссертации.

Оценка усталостной прочности конструкции ВЭУ выполнялась по методу допускаемых напряжений. Расчеты на усталостную прочность оптимизированной конструкции показали удовлетворительность предъявляемым требованиям и позволяют прогнозировать срок службы изделия не менее 20 лет.

На основе математической функциональной модели, построенной в программе - визуальном симуляторе функций У188т, проведен анализ воздействия возмущающих аэродинамических и центробежных сил и моментов на конструкцию ВЭУ, выявлены возможные резонансы, получены численные значения виброперемещений и воздействия на фундамент, а затем проведена оптимизация конструкции за счет регулирования жесткости и формы элементов, при необходимости с повторением предыдущих этапов (см. рис. 1) с целью смещения резонансных явлений в безопасные частоты.

1 - лопасть вертикальная

2 - кольцо верхнее

3 - кольпо среднее

4 - кольпо нижнее

5 - труба лопасти горизонтальной 6- основание ротора

7 - ось ступицы ротора

8 - мачта

9 - растяжка

- 9

Рис. 13. Расчетная модель ВЭУ-3

Основными факторами, вызывающими вибрации ВЭУ. передающиеся на фундамент, являются аэродинамические силы и моменты, действующие на ротор, а также периодические инерционные силы и моменты. В созданной модели «Ротор

- гибкая мачта» формульные зависимости представлены в виде структурных схем (примеры интерфейса приведены на рис. 14а и 146). Моделирование силового воздействия ВЭУ на фундамент проведено при различных скоростях вращения ротора при имитации воздействии ветра от 4 до 16 м/с. В процессе моделирования предполагалось, что ротор ВЭУ на выбеге проходит все возможные частоты вращения под аэродинамической нагрузкой.

V.,

Ї ;! —г-

Г VI V"-. ™ ~

і-, * й! г '* : їц Г :

1 • — -у**£ т& -

- г “ - ^ - * *- * Ш,1

І Г

Я*';:*-

а)

б)

Рис. 14. Пример интерфейса (а) и содержание блоков функций (б)

Поведение системы определяется системой дифференциальных уравнений сил и моментов с учетом силового воздействия, перемещения, трения, упругости:

М-х + 11'А' + К-х = Б: 3-ш +- Бй + В-ш = М

(22)

где М, Я. К, 1 О. В - матрицы коэффициентов, Р - матрица сил, М - матрица моментов.

Ввод данных в программу осуществлялся поблочно, с включением обратных связей и выводом графиков. Пульсационная составляющая скорости ветра считалась случайным квазистационарным процессом. Величины горизонтальной и вертикальной сил, действующих на фундамент мачты, можно оценить по формулам:

Р.<р= Руґор-Сі+Мх-С2; Р**пр= Лгг„-С, + Му С2; = Г2 , (23)

где Рх гор, Ру гор - составляющие горизонтальных сил, действующих на вершину мачты со стороны ротора, Р2 - сумма вертикальных сил, действующих на вершину мачты со стороны растяжек, Мх, Му - составляющие аэродинамического изгибающего момента, действующего на мачту со стороны ротора ветроколеса. Для принятых исходных данных по конструкции ВЭУ величины вычисленных коэффициентов Сі и С2 составляют: Сі =0,1; С2 = 0,12.

Пример на рис. 15 приведен для одной из сил Ру, действующей в горизонтальном направлении на мачту со стороны вала ротора. Постоянная составляющая

силы Ру = 10,8 Н. Гармоника с номером 1 имеет частоту вращения ротора. - частота вращения ротора; Т - время, с. По оси абсцисс на спектре Фурье отложена

безразмерная величина V кол № я, равная отношению частоты колебаний силы, приложенной к мачте ВЭУ Vкол , к частоте вращения ротора ук .

600

400

0

-400

-ВОО

-1200

-1600

-2000

-І-А>

т/

1- ¥-і........I

..ІД4

і..

п.

т. С

О. 04......0. И... О, 12.....0.16.....0,20.......0,24 0,28 И...

«мн

; і 20» Й .І 1 кол1-1 'к

1 ...

а) б)

Рис. 15. Зависимость силы Ру (а), действующей на мачту со стороны вала ротора и направленной перпендикулярно ветровому потоку, от времени Т. Спектр Фурье колебательной составляющей силы Ру (б)

Определение виброперемещений компонентов ротора и мачты наглядно демонстрируется графиками колебательных процессов, рис. 16.

а) *'*' б)

Рис. 16. Смещение Ьмнв нижней (а) и ЬИВЕ верхней (б) точек мачты в месте крепления соответственно к фундаменту и ступице в вертикальном направлении (мм) в зависимости от времени Т (с)

На основании исследований функциональной модели ВЭУ получены зависимости колебаний сил и моментов, частотные спектры и величина виброперемещений компонентов. Анализ спектра силовых воздействий проведен способом разложения действующих негармонических колебаний сил, моментов и виброперемещений в ряд Фурье с получением соответствующих спектров частот. На основании спектрального анализа сделаны выводы, на каких частотах вращения ротора ВЭУ ожидается резонанс, получена его амплитуда и преобладающие частоты, составлены таблицы резонансов.

Из анализа полученных данных следует, что резонанс колебаний элементов ВЭУ может вызываться как дисбалансом массы ротора, так и действием переменных аэродинамических сил. Резонансная частота может быть изменена за счет вариаций жесткости компонентов, то есть изменения частот собственных колебаний.

После многократных прогонов модели и итерационной оптимизации компонентов ВЭУ с корректировкой конструкторской документации создана таблица ожидаемых резонансов (табл. 3). Анализ таблицы показывает, что основные резонансы могут возникнуть на частотах 1 и 21 Гц. Частоту 1 Гц (60 об/мин) ротор проходит за считанные секунды, поэтому ее в целом можно считать несущественной.

Таблица 3

Ожидаемые резонансы ВЭУ в диапазоне частот вращения 0-180 об/мин

Номер интервала Начало малого интервала Инерци- онные Частоты возмушаюшнх сил Аэрсдпнамщес кие Примечание

г, с 12 15 у* 24 30 36

* л

1 156 3.525 21.1 42.26 63.40 84.53 105.66 126.8 1-Й резонанс узла подшипников

24: 29~2 17.8 35.66 53.50 71.33 89.16 Ют.О простые колебания

502 1.76’ 10.6 21.20 З1.*0 42.40 53.00 63.6 2-й резонанс узла полшитшкев

4 “53 1.007 6.0 і 2.8-І 13.13 24.17 30.21 36.2 1-Й резонанс мачты

5 1250 0.396 2.3" 4.75 7.12 9.50 11.87 14.2 простые

6 1682 0.167 1.02 2.00 3.01 4.01 5.01 6.01 2-Й резонанс

1900 0.108 0.65 1.30 1.94 2.59 3.24 3.89 простые колебания

2030 0083 0,50 1.00 1.50 2.00 2.50 3.00 3-й резонанс мачты

9 2230 0.056 0.34 0.67 1.04 1.34 1.67 2.01 4-П резонанс мачты

Частота 21 Гц проявляется в форме гармоник малой амплитуды и в целом не оказывает влияния на работу ВЭУ. Эксперименты и опыт эксплуатации ВЭУ подтверждают эти положения. Однако если ВЭУ установлена на инженерном сооружении (здании, мостовой конструкции и т.п.), резонанс на частоте собственных колебаний данного сооружения может привести к нежелательным последствиям. В связи с этим решение проблемы может быть осуществлено с помощью:

- метода особого регулирования мощности и частоты вращения ВЭУ;

- применения виброгасителей;

- проведения расчетов и изменения формы и/или материала мачты.

Ограничение максимальной частоты вращения ротора исследуемой ВЭУ, а

значит, минимизацию генератора и упрощение конструкции электронного регулятора, обеспечивают введенные аэродинамические регуляторы, заключенные в обтекатели аэродинамического профиля (рис. 17), активируемые на определенной частоте вращения ротора и ограничивающие дальнейшее увеличение частоты. Аэродинамические регуляторы имеют щитки с обтекателями, которые при возрастании центробежной силы совершают поворот вокруг оси траверсы на червячной передаче, тем самым снижая скорость вращения. Данное изобретение, защищенное патентом № 2347104 (2006 г.), позволило эксплуатировать ВЭУ на скорости

ветра более номинальной 11 м/с без применения балластной нагрузки и остановки ротора на буревых ветрах. Конструкция также включена в общую функциональную математическую модель, которая в данной работе не приводится в связи с большим объемом.

Поворотные щитки аэродинамического регулятора

Рис. 17. Аэродинамический регулятор вертикально-осевой ВЭУ

Разработка генератора включает в себя методы и стадии традиционного проектирования электромашины с проведением оптимизационных расчетов. В исследуемой ВЭУ генератор представляет собой вентильную электрическую машину переменного тока на постоянных магнитах с осевым (аксиальным) зазором с индуктором на постоянных магнитах и якорной обмоткой (рис. 18).

Индуктор на постоянных магнитах

Якорная обмотка

Индуктор на постоянных магнитах

Рис. 18. Генератор с осевым (аксиальным) зазором

Критериями оптимальности при разработке генератора являлись минимальная масса активных частей, минимальный объем активных частей, минимальная масса постоянного магнита, минимальная масса обмотки, минимальная суммарная стоимость постоянных магнитов и обмотки.

Программа оптимального проектирования, реализованная на языке программирования Ое1рЫ, представлена в виде блок-схемы на рис. 19 и положена в основу пакета программ оптимального проектирования генераторов исследуемых ВО ВЭУ. Постоянными (константами) в программе являются Р2 - выходная мощность; и - выходное напряжение; п - частота вращения; и т.д. Математическая модель имеет ряд ограничений. Оптимизация осуществляется перебором вариантов по методу Фибоначчи в сочетании с методикой координатного спуска, т.е. оптимальным делением интервалов значений независимых переменных с последующим анализом по критериям оптимальности.

После получения оптимальных значений производится поверочный расчет магнитного поля и потока в программном комплексе Ашув ЕМАв, а затем вычисляются параметры генератора в МаШСАО.

Метод разработки качественного генератора в общем случае разбивается на две взаимосвязанные задачи:

- анализ, то есть разработка методики поверочного расчета, при которой полностью известна геометрия и необходимо определить рабочие характеристики;

- синтез, то есть разработка методики оптимального проектирования, при которой для определенных техническим заданием параметров необходимо рассчитать геометрические параметры.

К

о

н

с

т

А

Н

т

ы

пых

ыа-гш-

Дании* ТЗ Р2.и,п,КИД.С

Парам калеггое к

...ПРС-РСДОР

Р2,и.п,С

Независимы*

П«*|1РМР1ШЫ^

Фиксированные размеры________

Ограничения

МиК Допус

Мин. Мкни к тнмал

я мдан ятти

иш:. ал. Р&ЗМ« схп»

тс>лш ТС»ЯЩ рк тогл

кк* кт

тэ

• ' 1 ' ' 1 г

Математн-

чегкяя

МОДСЛЬ

й«ау А ^ Число *

мае* йгмл Ш|>

г кщ.у ЛОЛЭО

•згл со г

Мин масса

Мнн объем

Мгек м*:с* «*гюгг»

I. и.*«гс*1гр'с*'акккж4

4«н.стокм»хт*

Блок-оитими1ато|>

К

р

и

т

Е

Р

И

И

Рис. 19. Схема алгоритма оптимального проектирования генератора

Несмотря на различные конструктивные исполнения, для анализа и синтеза этого типа генератора можно применить единый подход. Для расчета электромагнитной мощности активная часть генератора представляется в виде диска с равномерным токовым слоем, который пронизывается магнитными силовыми линиями. Элементарный участок этого диска создает электромагнитный момент:

с!М = Bi-Aj-R.i-dR.-da, (24)

где с!М - элементарный электромагнитный момент; В; - магнитная индукция элементарного участка; А; - линейная токовая нагрузка элементарного участка; Я; -радиус удаления элементарного участка от оси вращения; а - угол поворота.

Средний электромагнитный момент будет представлять собой поверхностный интеграл по площади активного кольца:

а/5г=|йГог| В1*а/*м*с1к=о.5*х*(°-+?А *{—^у\*в„-*Ахг (25)

О О 5£Х’ \ 2 ) \ 2 ] ’V/

где В, - средний электромагнитный момент; Бп - наружный диаметр активной части; Д- - внутренний диаметр активной части; В5Г - среднее значение индукции в воздушном зазоре; Аа - среднее значение линейной токовой нагрузки (определяется, как количество ампер-витков на единицу длины на среднем диаметре активной части генератора); сШ. - радиус удаления элементарного участка от оси вращения.

Электромагнитная мощность определяется умножением (25) на угловую частоту вращения машины

»!^)*в5г*А5г*±

2 60

(26)

где п - угловая частота вращения ротора, измеренная в об/мин.

Полученное уравнение представляет собой аналог классического уравнения Арнольда для электромашин цилиндрического типа с радиальным воздушным зазором. Средний диаметр генератора эквивалентен диаметру расточки якоря, а толщина кольца активной длине радиальной машины. Из уравнения Арнольда видно, что точность расчета при анализе зависит от точности определения средней линейной нагрузки и средней индукции. Определение средней линейной нагрузки связано с расчетом активного и индуктивного сопротивления фазы. Расчет активного сопротивления, как правило, не представляет сложности. Для расчета индуктивного сопротивления, а также для определения средней индукции необходим качественный расчет магнитного поля. Результаты расчетов приведены на рис. 20 п в табл. 4.

а)

Рис. 20. Графики зависимостей выходной мощности Р2(1) и электрических потерь в фазах Ре1(1) от тока I (а) и график выходного напряжения Щ) от тока I (б)

Таблица 4

Габаритная оптимизация (оптимизация по всем размерам)

На- Внут- Осе Чис- Высо- Тол- То.іш. Чис- Число Диа- Воз- Соп- Масса Мае

руж- рен- вая ло та щи- меди ло вит- меір душ- роїл- ге- са

ный ний дли- маг- маг- на (мм) виг- ков про- ный вле* не- маї-

диа- диа- на ни- ни- дис- ков в сек- вода зазор нме ра- нн-

метр диска (мм) метр диска (мм) (мм) тов на диске тов (ММ) ка (мм) в фазе ции (мм) (мм) фазы (Ом 10- ра (кт) тов (кі)

440 250 47 36 10 6,0 10 338 28 1,5 1,0 0,97 22 17,3

Номинальная точка:

Номинальный ток (А) 16,66

Номинальное напряжение (В) 180

Номинальная мощность (Вт) 3000

Номинальная частота вращения (об/мин) 180

КПДв номинальной точке 0,847

Плотность тока(А/мм") 8,9

Использование данного метода в разработке генератора является перспективным для ВО ВЭУ в связи с введением аэродинамического регулирования (или ограничения мощности генератора сверху).

Разработка электронного регулятора мощности исследуемой ветроэнергетической установки имеет следующие особенности:

- регулирование мощности ВЭУ в диапазоне частот вращения ротора от 40 до 180 об/мин с поддержанием максимального КИЭВ;

- частота вращения ротора ВЭУ ограничивается на 180 об/мин за счет аэродинамических регуляторов и далее не поднимается;

- отсутствие балластной нагрузки.

В связи с этим задачей разработки и оптимизации являлось снижение потерь при преобразовании и максимально полное использование энергии ветра. Таким образом, для получения максимума мощности требуется придерживаться определенной линии, называемой «линией оптимальной мощности» на графике изолиний мощности ротора ВЭУ (рис. 10), которая и выступает в качестве исходных данных при разработке алгоритма электронного преобразователя энергии. Критерием оптимальности при разработке регулятора являлся максимум выходной мощности на основе регулирования параметров выхода с отслеживанием располагаемой мощности ротора. Регулирование электрической мощности ВЭУ осуществлялось с помощью импульсного преобразователя с выбором оптимальных рабочих режимов на основе метода, состоящего из следующих параллельно выполняемых действий:

- измерение внешних параметров (располагаемой мощности, мощности нагрузки, ток заряда аккумуляторной батареи, частота тока и напряжение в фазе генератора и др.);

- определение выхода параметров за пределы безопасных режимов и отключение генератора с переводом преобразователя в режим защиты;

- обеспечение заряда аккумулятора в режиме постоянного напряжения с ограничением тока заряда при условии наличия достаточной генерируемой мощности с предотвращением перезаряда и переразряда аккумулятора.

Схема регулятора построена на основе программируемого микроконтроллера. Для настройки, диагностики и программирования регулятор имеет интерфейс КБ-232. Регулятор построен по схеме понижающего импульсного преобразователя, с несколькими контурами обратной связи.

Алгоритм регулирования мощности, отбираемой от генератора ветроэнергетической установки:

1. Измерение напряжения и тока в цепи нагрузки, расчет мгновенной потребляемой мощности.

2. Измерение частоты тока в фазе генератора и расчет частоты вращения ротора ветроколеса.

3. Определение мощности ротора на основе сравнения частоты вращения ротора с таблицей мощностей ротора.

4. Определение максимально допустимой мощности зарядки аккумуляторной батареи с использованием максимально допустимого тока зарядки батареи и напряжения на ней.

5. Расчет суммарной требуемой мощности для данного момента времени.

6. Сравнение допустимой и требуемой мощностей, определение минимальной из них с регулированием выходной мощности на основе ШИМ-модуляции.

7. С помощью цифровых потенциометров в цепях обратных связей регулятора задание действующей располагаемой мощности на основе предыдущих расчетов. Общая схема регулятора представлена на рис. 21.

Применение данного метода регулирования мощности ВЭУ позволяет значительно повысить уровень удельной выработки энергии ВО ВЭУ в связи с поддержанием максимального КИЭВ во всем диапазоне скоростей ветра до 45 м/с.

В Главе 3 «Разработка технологии изготовления ВЭУ» даны описания технологических процессов изготовления лопастей и генератора, а также приведена методика вибробалансировки ротора.

С целью повышения надежности и увеличения срока службы лопастей ВЭУ разработан технологический процесс производства изделий на основе композитных материалов. Цель разработки - получение новых недорогих конструкционных материалов, облегчение важнейших узлов ВЭУ, повышение эксплуатационных свойств изделия в условиях воздействия больших температурных амплитуд и высоких механических нагрузок. Основными задачами технологического процесса являются получение заготовок и изделий с использованием объемных материалов, монофазное упрочнение поверхностного слоя, структурирование функционального покрытия, сверхпластичная формовка и диффузия. Одна из ключевых задач -получение ультрамелкозернистой структуры с минимумом шероховатости. Изготовление лопасти происходит в оснастке особым способом намотки стеклоткани на распределенную заготовку за одну формовку, с использованием ряда ноу-хау. Изделие монолитно, готово к сборке в составе ВЭУ и включает в себя непосредственно лопасть, основание лопасти (крепление к несущему кольцу) и заделку лопасти (крепление связующего кольца или стяжек). Причем установочный угол выставлен в оснастке. Оснастка для изготовления лопасти приведена на рис. 22а.

Прочностные характеристики материала 8КСУ-ПЭ-300 в среднем на 25% лучше традиционных стеклопластиков:

ов1=620 МПа - предел прочности при растяжении в направлении основы; аВ2=390 МПа- предел прочности при растяжении в направлении утка; о_В1=370 МПа - предел прочности при сжатии в направлении основы;

о_В2=310 МПа- предел прочности при сжатии в направлении утка;

Е!= 31000 МПа - модуль упругости в направлении основы.

Изделия имеют шестикратный запас по прочности в составе ВЭУ, причем вес и себестоимость снижены в 2 раза, а изделие не содержит металлических частей, т.е. не экранирует радиоволны.

Технологический процесс изготовления генератора разработан в части монолитного отверждения якорной обмотки (рис. 226) и снижения общих вибраций генератора.

Рис. 22. Оснастка для изготовления лопасти (а) и обмотка генератора (б)

Технологический процесс заливки обмотки разработан совместно с корпорацией Empire Magnetics, Inc. (США) и использует центрифугу с установленной внутри обмоткой в технологическом кожухе, из центра которого при вращении под действием центробежных сил материал заливки растекается и заполняет пространство внутри кожуха, включая обмотку. Материал заливки изготавливается на основе цементирующего материала или полиэфирных смол.

Методика вибробалансировки ротора разработана для снижения общих вибраций системы «ротор-мачта» и представляет собой трехступенчатую последовательность. На первом этапе балансировки проводится взвешивание компонентов ротора (лопастей, колец, траверс с аэродинамическими регуляторами и стяжек) на предмет формирования комплектов с одинаковой массой. Затем проводится статическая балансировка с контрольным взвешиванием. На втором этапе проводится цеховая сборка ротора на стенде и проверяется перпендикулярность установки лопастей плоскости среднего кольца. При наличии отклонения лопасти от перпендикулярности положение лопасти выравнивается с помощью добавления регулирующих пластин в месте крепления кронштейна торца лопасти к несущему среднему кольцу. На третьем этапе производится подключение генератора ВЭУ к источнику питания, и генератор закручивается в режиме двигателя (вторым вариантом может быть использование привода для раскрутки ротора). Дисбаланс ротора регистрируется с помощью модифицированного виброанализатора К-4102М, после чего проводится последующая динамическая весовая балансировка. Для достижения лучшего результата необходимо принять состояние ротора, отбалансированного в одной плоскости, за исходное состояние и провести вибробалансировку в следующей плоскости. Полный цикл включает в себя балансировку в плоскостях среднего, верхнего и нижнего колец или стяжек согласно конструкции.

Вибрации генератора снижаются за счет виброгасителей, расположенных между якорной обмоткой и валом, между фланцами мачты, а также дополнительных традиционных и разработанных устройств - демпферов растяжек и мачты.

(а)

(б)

В Главе 4 «Экспериментальные исследования» описаны результаты экспериментальной эксплуатации исследуемых ВЭУ, подтверждающие теоретические расчеты, приведены данные исследований шума и вибраций ВЭУ.

Семейство малых вертикально-осевых ветроэнергетических установок мощностью 0.1 - 30 кВт разработано для автономных потребителей электроэнергии. В общей сложности полевые испытания прошло 9 различных по конструкции и мощности модификаций ВЭУ. Примеры приведены на рис. 23. Почти все установки продемонстрировали КИЭВ до 38%.

1 кВт (Челябинск) 3 кВт (США) 3 кВт (Кумертау) 30 кВт (Челябинск)

Рис. 23. Семейство вертикально-осевых ВЭУ

Замеры шума и вибраций нескольких образцов ВЭУ-З(б) (3 кВт, 6 лопастей), проведенные в России и США на расстоянии 10 м от мачты ВЭУ при скоростях ветра 6-8 м/с, во всех случаях показали, что уровень шума окружающего фона выше, чем шум работающей под нагрузкой ВЭУ. Имеются акты соответствующих организаций.

Измерение вибраций проведено на нескольких ВЭУ по разработанной методике. Результат показывает, что данные несколько различаются, однако в целом удовлетворяют нормам. Пример измерений ВЭУ на скорости ветра 8 м/с при нагрузке 1 кВт (прожектор) приведен на рис. 24. Частота вращения составляет 124 об/мин (или 2,06 Гц). Акселерометр установлен на нижний сегмент мачты.

Из рис. 24 видно, что в спектре присутствуют резонансные частоты 2; 6; 12; 21; НО Гц. Гармоника с частотой 2 Гц характеризует дисбаланс непосредственно ротора, т.е. частоту вращения ротора; частота 6 Гц - частоты колебаний трех лопастей каждого яруса ВЭУ; 12 Гц - воздействия аэродинамических сил поочередно на 6 лопастей, 21 Гц - частота, близкая к собственной частоте колебаний подшипникового узла. Частота 110 Гц объясняется наличием вибраций генератора. Поскольку генератор имеет 36 магнитов на верхнем и столько же на нижнем основании (табл. 3), т.е. количество пар магнитных полюсов р=18, с учетом наличия трехфазной обмотки, т.е. т=3, частота пульсаций момента генератора Мгсн будет больше частоты вращения ротора Ырот (Гц) в (р-т = 18-3 = 54) раза. Тогда:

Н-ен = ^0Т • 54 = 2,06 • 54 = 111 Гц. Это значение близко к частоте 110 Гц и говорит о присутствии данной гармоники в связи с наличием нагруженного генератора.

О»»**! Км-Ле»

Снггплг*. ^8ЙЙ|«е«0)М1СГЬ Мроитм. .*-« ».«г исг-г: * я *■- - 5^,4»^ вж

1

I 1 вш %>

« 10 гв у А- згь

вЛОГЗЕЗ . _ , агав

: ? 1 <I §<£ ф

«1 1.: «ч *•» гл Ш'як'ГУ' *»*•»- «« езсгжг ЯШ

М«МЙМШ. КК X «йя/п { 8.Ж Срвяиме Ов0<Х»М. 0 1 0.15- ф

* А} 2* 10 -Х> Й? Ф & !>Э А1Я ж я

Рис. 24. Результаты измерений ВЭУ-3 под нагрузкой прибором К-4102М

Таким образом, в целом теоретические расчеты полностью сходятся с экспериментальными данными. Дальнейшее увеличение скорости ветра не влияет на рост виброскоростей при отсутствии резонансов.

Шум в данном эксперименте не измерялся. В помещении под мачтой некоторый шум присутствовал постоянно, однако на расстоянии 1 метра от основания мачты шум фона перекрывал шум работающей ВЭУ.

Сравнение максимальных значений амплитуды виброскорости (0,25 мм/с или

0,025Т0~2 м/с на частоте 6 Гц, а также 0,019-10~2 м/с на частоте 3 Гц) с табл. 11 ГОСТ 12.1.012-90 («Вибрационная безопасность» Санитарные нормы спектральных показателей вибрационной нагрузки на оператора. Общая вибрация, категория 3, тип "в") дает следующий результат (табл. 5):

Таблица 5

Нормативные значения виброскорости по ГОСТ 12.1.012-90 и эксперимент

Табличные данные ГОСТ 12.1.012-90 Экспериментальные данные

Среднегеометрические частоты полос, Гц Виброскорость, м/с, в 1/1 окт. Частота. Гц Виброскорость, м/с, в 1/1 окт.

3,15 0.063-10“2 3,0 0.019-10~2

6,3 0,032-10“2 6.0 0,025-10~2

Из таблицы 5 следует, что ВЭУ можно размещать в непосредственной близости к производственным зданиям, или даже на них. Аналогичное сравнение экспериментальных данных с СН 2.2.4/2.1.8.566-96 («Производственная вибрация, вибрация в помещениях жилых и общественных зданий», таблица 9 «Допустимые значения вибрации в жилых помещениях, палатах больниц, санаториев») позволяет размещать ВЭУ рядом с жилыми зданиями и непосредственно на них.

В Главе 5 «Экономическая оценка рынка ВЭУ» описывается расчет потенциального рынка малых ветроэнергетических установок, даются рыночные прогнозы, экологические и социальные перспективы внедрения ВЭУ до 30 кВт.

Согласно современным требованиям проведения НИОКР в части коммерциализации вопрос внедрения результатов исследовании является в Российской Федерации одним из стратегических. В соответствии с этим проведение анализа рынка разрабатываемых и создаваемых изделий является побуждающей экономической причиной для создания методологии разработки ветроэнергетических установок на основе научных методов и подходов.

Расчет потенциального рынка малых ветроэнергетических установок произведен исходя из анализа ряда положений, обосновывающих рост износа генерирующих мощностей и необходимости внедрения малой энергетики на основе возобновляемых источников энергии с максимальным приближением электрических и тепловых генераторов к местам потребления. Исходя их этих положений потенциальный рынок для малой ветроэнергетики со стоимостью оборудования 60.000 руб. за кВт мощности ВЭУ прогнозируется следующим образом:

1. Экономически выгодная замена малых дизельных станций агропромышленного комплекса Севера на ВЭУ, по меньшей мере в размере 50%, создаст рынок в объеме 250 тыс. кВт, что эквивалентно нескольким сотням тысяч малых ВЭУ или ориентировочно 250 млрд. руб.

2. При условии планируемой Правительством РФ до 2020 г. замены хотя бы 1% традиционных источников энергии на ветроэнергетику, объем установленной мощности ветроэнергетических установок может составить 2,17 млн. кВт. Из них, по опыту внедрения развитых стран, малые ВЭУ до 30 кВт могут составлять по меньшей мере 5% всего парка, или 100 тыс. кВт установленной мощности, что эквивалентно несколько сотням тысяч малых н сверхмалых ВЭУ, т.е. более 100 млрд. руб.

Таким образом, даже не принимая во внимание объем международного рынка, рост спроса на автономные источники электроэнергии, развитие научнотехнических решений в ветроэнергетике, снижение себестоимости ВЭУ, тенденции развитых стран в области развития малой энергетики и, наконец, пренебрегая данными прогнозов по грядущему энергодефициту, можно смело утверждать, что потенциальный рынок малых ВЭУ в России составляет уже в настоящее время свыше 350 млрд. руб. Решение данной проблемы является перспективной народно-хозяйственной задачей.

ОСНОВНЫЕ ВЫВОДЫ И РЕЗУЛЬТАТЫ РАБОТЫ

1. Научно обоснован методологический подход к разработке ВО ВЭУ.

2. Разработанная конструкция ВО ВЭУ смоделирована в трехмерной и функциональной моделях, оптимизирована покомпонентно по ряду критериев с целью получения максимальной выработки энергии. В результате применения последовательности итерационных методов покомпонентной оптимизации экспериментально полученный КИЭВ ротора составил 36,5% во всем диапазоне быстроходности с сохранением устойчивой работоспособности на скоростях ветра до 45 м/с с прогнозируемой буревой скоростью 60 м/с.

3. Аэродинамическое регулирование частоты вращения ротора и электронное регулирование мощности ВЭУ позволило поддерживать максимальный КИЭВ во всем диапазоне частот вращения, ограниченным до 180 об/мин.

4. Разработаны технологические процессы изготовления лопасти, генератора, аэродинамического регулятора, ротора ВО ВЭУ. Прочность материала лопастей увеличена на 25%, себестоимость снижена в 2 раза.

5. Разработанная методика вибробалансировки позволила снизить вибрации до уровня 0,019-1(Г2 м/с на частоте 3 Гц, что с запасом удовлетворяет требованиям ГОСТ 12.1.012-90 и СН 2.2.4/2.1.8.566-96. Разработанная методика может быть использована для большинства вертикально-осевых ветроэнергетических установок.

6. Проведенные лабораторные, стендовые и полевые испытания разработанных ВО ВЭУ подтвердили теоретические расчеты мощности, шума и вибраций. Выходная мощность ВЭУ-3 соответствует расчетным характеристикам и составляет 3,5 кВт на скорости ветра 11 м/с. Уровень шума работающей ВЭУ не превышает 58 дБ(А) на расстоянии 10 метров и с запасом удовлетворяет ГОСТ 23337-78 и СН 2.2.4/2.1.8.562-96.

7. Результаты теоретических и экспериментальных испытаний вертикальноосевых ВЭУ показали необходимость корректировки действующих стандартов в части разработки и эксплуатации вертикально-осевых ветроустановок, с уточнением требований к вибробезопасности.

8. На основе совокупности теоретических и экспериментальных результатов исследований создана методология разработки вертикально-осевых ветроэнергетических установок, позволяющая структурировать и логически организовать методы и средства оптимального проектирования.

9. Потенциальный рынок малых ВЭУ в Российской Федерации оценен в общую сумму свыше 350 млрд. рублей. За счет применения методологии, основанной на методах оптимизации, получена конструкция, которую можно эксплуатировать на урбанизированных территориях, в том числе и АПК с одновременным снижением установочных и эксплуатационных расходов соответственно на 25% и 15%. Себестоимость ВЭУ без введения автоматизированного производства составила от 60 руб. за Ватт мощности оборудования.

10. Создано семейство ВО ВЭУ мощностью 0,1; 1; 3; 5; 30 кВт автономного применения. ВЭУ мощностью 0,1 и 3 кВт внедрены в систему энергоснабжения фермерских хозяйств.

11. В работе решена крупная научная проблема создания универсальной методологии разработки и создания вертикально-осевых ветроэнергетических установок, что имеет важное народно-хозяйственное значение для повышения энергонезависимости и обороноспособности страны.

Основное содержание диссертации опубликовано в следующих работах:

А. В изданиях, рекомендованных Высшей Аттестационной Комиссией:

1. Соломин, Е.В. Ветроэнергетические установки «ГРЦ-Вертикаль» / И.М. Кирпичникова, В.П. Кривоспицкий, Е.В. Соломин // Вестник МАНЭБ. Приложение «По Материалам I Международной научно-практической конференции «Ресурсосбережение и возобновляемые источники энергии: экономика, экология, опыт применения». - СПб.-Чита, 2008. - Т. 13. - № 3. - С. 129-134.

2. Соломин, Е.В. Ветроэнергетическая установка с вертикальной осью вращения / И.М. Кирпичиикова, Е.В. Соломин // Вестник ЮУрГУ. Серия «Энергетика».

- 2008. - Вып. 10. - № 26. - С. 15-16.

3. Соломин, Е.В. Сопряжение ВЭУ малой мощности с пленочным электронагревателем для обогрева помещений / И.М. Кирпичиикова, И.Н. Панасюк, Е.В. Соломин // Вестник ЮУрГУ. Серия «Энергетика». - 2009. - Вып. 12. - № 34. -С. 78-81.

4. Соломин, Е.В. Сделано в России! Ветроэнергетические установки с вертикальной осью вращения мощностью 1...100 кВт для потребителей электроэнергии / Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2009. -№ 10. - С. 74-78.

5. Соломин, Е.В. Виброгасителн мачт сверхмалых вертикально-осевых ветроэнергетических установок / И.М. Кирпичиикова, Е.В. Соломин // Вестник ЮУрГУ. Серия «Энергетика». - 2010. - Вып. 13. - № 14(190). - С. 78-81.

6. Соломин, Е.В. Инженерный метод и математическое моделирование в проектировании ветроэнергетических установок / Ю.В. Грахов, О.В. Матвеенко, Е.В. Соломин // Вестник ЮУрГУ. Серия «Математика. Механика. Физика». -2010. - Вып. 2. - № 9(185). - С. 45-52.

7. Соломин, Е.В. Экономические аспекты внедрения ветроэнергетических установок / Е.В. Соломин // Вестник ЮУрГУ. Серия «Экономика и менеджмент».

- 2010. - Вып. 14. - № 20(196). - С. 32-36.

8. Соломин, Е.В. Предложение о сотрудничестве при производстве ветроэнергетических установок / Дж.Куль, Е.В. Соломин // Альтернативная энергетика и экология. -М.: НИИЭС, 2010. -№ 1.-С. 115-125.

9. Соломин, Е.В. Технические особенности и преимущества ветроэнергетических установок / Е.В.Соломин, Р.Л.Холстед // Альтернативная энергетика и экология. - М.: НИИЭС, 2010. - № 1. - С. 36-41.

10. Соломин, Е.В. Ветроэнергетические установки ГРЦ-Вертикаль / Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2010. - № 1. -С. 10-15.

11. Соломин, Е.В. Система освещения, основанная на ветроэнергетической установке / А.С. Мартьянов, Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2010. - № 1. - С. 101-105.

12. Соломин, Е.В. Преобразование энергии в ветроэнергетических установках / И.М. Кирпичиикова, А.С. Мартьянов, Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2010. - № 1. - С. 93-97.

13. Соломин, Е.В. Контроллер заряда ветроэнергетической установки / А.С. Мартьянов, Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2010. -№ 1.-С. 106-109.

14. Соломин, Е.В. Отопительная система на базе ВЭУ / Н.В. Пронин, Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2010. - № 1. -

С. 110-114.

15. Соломин, Е.В. Сравнительные характеристики вертикально-осевых ветроэнергетических установок / Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2010. - № 1. - С. 48-53.

16. Соломин, Е.В. Ветроэнергетическая экономика / Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2010. - № 2. - С. 28-30.

17. Соломин, Е.В. Ветроэнергетические установки ГРЦ-Вертикаль для экологических курортов / Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2010. - № 2. - С. 60-64.

18. Соломин, Е.В. Размещение ветроэнергетических установок с вертикальной осью вращения на ответственных инженерных объектах / Е.В. Соломин // Вестник ЮУрГУ. Серия «Строительство и архитектура». - 2010. - Вып. 11. -№33(209).-С. 47-51.

19. Соломин, Е.В. Основы методологии разработки вертикально-осевых ветроэнергетических установок / Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2011.-№ 1.-С. 18-28.

20. Соломин, Е.В. Система отопления на основе ветроэнергетической установки и теплового аккумулятора / А.С. Мартьянов, Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2011. - № 2. - С. 30-33.

21. Solomin, E.V. Development and optimizing of vertical axis wind turbines / E.V. Solomin // Альтернативная энергетика и экология. - М.: НИИЭС, 2011. - № 1.

- С. 29-39.

22. Соломин, Е.В. Итерационная оптимизация параметров и режимов работы вертикально-осевых ветроэнергетических установок / Е.В. Соломин // Вестник ЮУрГУ. Серия «Энергетика». - 2011. - Вып. 15(232). - С. 73-81.

23. Соломин, Е.В. Разработка математической модели ветроэнергетической установки мощностью 3 кВт производства ООО «ГРЦ-Вертикаль в пакете Matlab / Н.В. Пронин, А.С. Мартьянов, Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2011 - №5. - С.41-44.

24. Соломин, Е.В. Методы торможения и регулирования частоты ветроэнергетических установок / А.Н. Киндряшев, Е.В. Соломин // Альтернативная энергетика и экология. - М.: НИИЭС, 2011. - № 5. - С. 38-40.

25. Соломин, Е.В. Разработка математической модели ветроэнергетической установки мощностью 3 кВт производства ООО «ГРЦ-Вертикаль в пакете Matlab / Н.В. Пронин, А.С. Мартьянов, Е.В. Соломин // Вестник ЮУрГУ. Серия «Энергетика». - 2011. - Вып. 7. -№ 35(211). - С. 49-52.

26. Соломин, Е.В. Перспективы использования малых ветроэнергетических установок в агропромышленном комплексе / Е.В. Соломин // Механизация и электрификация сельского хозяйства. - М., 2011. - Вып. 7. - С. 12-15.

27. Соломин, Е.В. Эффективность автономного электроснабжения фермерского хозяйства / А.С. Мартьянов, Е.В. Соломин // Механизация и электрификация сельского хозяйства. - М., 2011. - Вып. 9. - С. 29-30.

28. Solomin, E.V. Joint scientific research of Russian and German scientists in renewable energy / E.V. Solomin // Альтернативная энергетика и экология. - М.: НИИЭС, 2011 -№ 10.-С. 71-82.

29. Solomin, E.V. Photovoltaic or wind solution / E.V. Solomin // Альтернативная энергетика и экология. - М.: НИИЭС, 2011. - № 11. - С. 83-89.

Б. В монографиях, брошюрах, учебных нособннх

1. Соломин, Е.В. Нетрадиционные и возобновляемые источники энергии: учебное пособие / И.М. Кирпичникова, Соломин, Е.В. - Челябинск: Изд-во /ЮУрГУ, 2009,- 167 с.

2. Соломин, Е.В. Возобновляемые источники энергии: учебное пособие к практическим занятиям / И.М. Кирпичникова, Соломин, Е.В. - Челябинск: Изд-во ЮУрГУ.-2009.-50 с.

3. Соломин, Е.В. Методология разработки и создания вертикалыю-осевых ветроэнергетических установок: монография / Соломин, Е.В. - Челябинск: Изд-во ЮУрГУ, 2011.-324 с.

В. В патентах, ноу-хау и свидетельствах ФИПС

1. Пат. 2244996 Российская Федерация, МПК7 Н 02 К 19/16 1/06. Генератор переменного тока / С.А.Ганджа, Е.В. Соломин и др. - № 2003124088/09; заявл. 31.07.03; опубл. 20.01.2005, Бюл. №2.-3 с.

2. Пат. 2347104 Российская Федерация, МПК РОЗО 3/06 (2006.1). Ротор ветряной установки с вертикальной осью вращения / Ю.В.Грахов, Е.В.Соломин и др. -№ 2006117014/06; заявл. 12.05.2006; опубл. 20.02.2009, Бюл. № 5. - 12 с.

3. Пат. 87767 Российская Федерация, МПК Р16Р 15/133 (2006.01). Амортизатор растяжки ветряной установки / И.М.Кирпичникова, Е.В.Соломин и др. -№ 2008138451/22; заявл. 26.09.2008; опубл. 20.10.2009, Бюл. № 29. - 3 с.

4. Пат. 100851 Российская Федерация, МПК Н01Ь 33/00 (2010.01). Комбинированная ветросолнечная энергетическая установка со светодиодным излучателем для освещения общественных мест / Е.В.Соломин - № 2010137760/28; заявл. 10.09.2010; опубл. 27.12.2010, Бюл. № 36. - 2 с.

5. Пат. 101105 Российская Федерация, МПК Р03Б 9/02 (2006.01). Комбинированная ветросолнечная энергетическая установка со светодиодным излучателем для освещения общественных мест / Е.В.Соломин - № 2010137417/06; заявл. 08.09.2010; опубл. 10.01.2011, Бюл. № 1. -2 с.

6. Секрет производства (ноу-хау) «Способ регулирования отбора мощности генератора ветроэнергетической установки» / И.М.Кирпичникова, А.С.Мартьянов, Е.В.Соломин. - Приказ № 60 от 18.02.2011 г. - ГОУ ВПО «Южно-Уральский государственный университет». - 7 с.

7. Пат. 103168 Российская Федерация, МПК Р2189/02 (2006.01). Автономное осветительное устройство / Е.В.Соломин - № 2010138386/07; заявл. 17.09.2010; опубл. 27.03.2011, Бюл. №9.-2 с.

8. Пат. 108504 Российская Федерация, МПК РОЗБ 3/06 (2006.01). Монолитная лопасть ветроэнергетической установки / Е.В.Соломин - № 2011111177/06; заявл. 24.03.2011; опубл. 20.09.2011, Бюл. № 26. -2 с.

Г. В международных, российских, межрегиональных изданиях и материалах конференций

1. Соломин, Е.В. Описание ВЭУ ООО «ГРЦ-Вертикаль» / Е.В. Соломин // Деловой Север. - Екатеринбург, 2008. Вып. 5. - №5. - С.8-9.

2. Соломин, Е.В. Ветроэнергетика для дома и офиса. Ветроэнергетические установки с вертикальной осью вращения мощностью 1... 100 кВт для дома и про-

мышленности / Ю.В. Грахов, Е.В. Соломин и др. // Материалы V Международной научно-практической конференции «Возобновляемые источники энергии. Ресурсы. Системы энергогенерирования на возобновляемых источниках энергии». - М.,

2008.-С. 37-39.

3. Соломин, Е.В. Программно-математическая модель ветроэнергетической установки с вертикальной осью вращения / Ю.В. Грахов, И.М. Кирпичникова, Е.В. Соломин // Материалы V Международной научно-практической конференции «Возобновляемые источники энергии. Ресурсы. Системы энергогенернрования на возобновляемых источниках энергии / М., 2008. - С. 41-42.

4. Соломин, Е.В. Вертикально-осевые ветроэнергетические установки: «Сделано в России» / Е.В. Соломин // Материалы Межрегионального совета по науке и технологиям. «Механика и процессы управления» и «Проблемы машиностроения». Екатеринбург-Миасс: УрО РАН, 2008. - С. 17-19.

5. Соломин, Е.В. Ветроэнергетические установки «ГРЦ-Вертикаль» / Е.В. Соломин // Малая энергетика. М.: Изд-во НИИЭС, 2008. - №3(8). - С. 57-61.

6. Соломин, Е.В. Ветроэнергетические установки ГРЦ-Вертикаль. Перспективы развития // 60-я юбилейная научная конференция, посвященная 65-летию Южно-Уральского государственного университета. Секция "Технические науки". Материалы конференции. - Челябинск: Изд-во ЮУрГУ, 2008. - Т. 1. - С. 23-26.

7. Соломин, Е.В. Социальные, экономические и правовые аспекты размещения ветроустановок на зданиях и сооружениях / Е.В. Соломин, И.М. Кирпичникова // Малая энергетика, М.: Изд-во НИИЭС, 2009. - №1. - С. 56-61.

8. Соломин, Е.В. Методика балансировки ротора ветроэнергетической установки с вертикальной осью вращения / Е.В. Соломин, И.М. Кирпичникова // Малая энергетика, М.: Изд-во НИИЭС, 2009. - №1. - С. 48-52.

9. Соломин, Е.В. Вертикально-осевые установки российской ветроэнергетики / Соломин Е.В. // Материалы Регионального научно-практического семинара «Современное состояние, проблемы и перспективы использования возобновляемых источников энергии» 8-9 октября 2009 г. - Элиста: Изд-во КГУ, 2009. - С. 71-75.

10. Соломин, Е.В. Размещение ветроэнергетических установок мощностью до

10 кВт на кровле зданий / Е.В.Соломин // Материалы IV Международной научнопрактической конференции «Безопасность жизнедеятельности в третьем тысячелетии» 1-3 октября 2009 г. - Челябинск: Изд-во ЮУрГУ, 2009. - С. 284-290.

11. Соломин, Е.В. Социальные, экономические и правовые аспекты размещения ветроустановок на зданиях и сооружениях / Е.В.Соломин, И.М.Кирпичникова // Малая энергетика. - М.: НИИЭС, 2009. - №1-2. - С. 29-39.

12. Соломин, Е.В. Методика балансировки ротора ветроэнергетической установки с вертикальной осью вращения / И.М.Кирпичникова, Е.В.Соломин // «Малая энергетика». - М.:НИИЭС, 2009. - №1-2. - С. 39-43.

13. Solomin, Е. High Efficient Low Noise Vertical Axis Wind Turbines / E.Solomin // 21st ISTC Korea Workshop “Renewable Energy”. Pusan National University, Korea, 2009. - №1-2. - P. 38^14.

14. Соломин, Е.В. Методика вибробалансировки ветроэнергетических установок с вертикальной осью вращения / И.М. Кирпичникова, Е.В. Соломин // Энергетическая политика. - М., 2009 - №5. - С. 72-78.

15. Соломин, Е.В. Энергосберегающая технология отопления с использованием ветроэнергетических установок / Н.В.Пронин, Е.В.Соломин // Энерго- и ресурсосбережение. Нетрадиционные и возобновляемые источники энергии. Сборник материалов Всероссийской студенческой олимпиады 16-19 ноября 2009 г., научно-практической конференции и выставки студентов, аспирантов и молодых ученых 14-18 декабря 2009 г. Екатеринбург. - Екатеринбург: Изд-во УГТУ-УПИ, 2009. - С. 492-495.

16. Соломин, Е.В. Гашение вибраций ветроэнергетической установки с вертикальной осью вращения / Е.В.Соломин // Материалы 61 научно-практической конференции ЮУрГУ, секции технических наук. - Челябинск: Изд-во ЮУрГУ,

2009.-Т. 2. - С. 276-280.

17. Соломин, Е.В. Обеспечение вибробезопасности ветроэнергетических установок с вертикальной осью вращения / И.М. Кирпичникова, Е.В. Соломин // Материалы конференции по ВИЭ 9-10 июня 2009 г. Москва. - М.,2009. - С. 23-27.

18. Соломин, Е.В. Использование ветроэнергетических установок для обеспечения энергосберегающей системы отопления жилых помещений / И.М. Кирпичникова, Е.В. Соломин, И.Н. Панасюк. Н.В. Пронин // Материалы научнопрактического семинара 14-16 апреля 2010 г. Астрахань. - Астрахань: Издат. Дом «Астраханский университет», 2010. - С. 65-68.

19. Соломин, Е.В. Ветроэнергетические установки с вертикальной осью вращения / И.М. Кирпичникова, В.П. Кривоспицкий, Е.В. Соломин // Известия Академии Электротехнических наук РФ. - М., 2010. - №2. - С. 68-72.

20. Соломин, Е.В. Инновационное будущее за регионами / Е.В.Соломин // Debts and financing economic journal. - М., 2010. - №7(44). - С. 68-69.

21. Соломин, Е.В. Ветроэнергетические установки широкого применения / Е.В.Соломин II Материалы Форума «Изменение климата и экология промышленного города». - Челябинск, 2010. - С. 65-66.

22. Соломин, Е.В. Экокультура и фитобиотехнологии улучшения качества жизни на Каспии / И.М. Кирпичникова, Е.В. Соломин, И.Н. Панасюк, Н.В. Пронин // сборник трудов Международной конференции с элементами научной школы для молодежи 7-10 декабря 2010 г. Астрахань. - Астрахань: Издат. Дом ((Астраханский университет», 2010. - С. 25-28.

23. Соломин, Е.В. Контроллер заряда АКБ ветроэнергетической установки / А.С. Мартьянов, Е.В.Соломин // Материалы 62-й научной конференции секции технических наук. - Челябинск: Изд-во ЮУрГУ, 2010. - Т. 3. - С. 123-126.

Издательский центр Южно-Уральского государственного университета

Подписано в печать 22.11.2011. Формат 60x84 1/16. Печать цифровая. Уел. печ. л. 2,32. Уч.-изд. л. 2. Тираж 100 экз. Заказ 392/669.

Отпечатано в типографии Издательского центра ЮУрГУ. 454080, г. Челябинск, пр. им. В.И. Леннна, 76.