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

кандидата технических наук
Буюкли, Татьяна Васильевна
город
Москва
год
2011
специальность ВАК РФ
05.07.05
цена
450 рублей
Диссертация по авиационной и ракетно-космической технике на тему «Моделирование динамического поведения лопаток компрессоров авиационных двигателей в нестационарном потоке воздуха»

Автореферат диссертации по теме "Моделирование динамического поведения лопаток компрессоров авиационных двигателей в нестационарном потоке воздуха"

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

Буюкли Татьяна Васильевна

МОДЕЛИРОВАНИЕ ДИНАМИЧЕСКОГО ПОВЕДЕНИЯ ЛОПАТОК КОМПРЕССОРОВ АВИАЦИОННЫХ ДВИГАТЕЛЕЙ В НЕСТАЦИОНАРНОМ

ПОТОКЕ ВОЗДУХА

Специальность 05.07.05 «Тепловые, электроракетные двигатели и энергоустановки летательных аппаратов»

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

2 О ОПТ 2011

Москва-2011

4857594

Работа выполнена на кафедре «Теория воздушно-реактивных двигателей» Московского авиационного института (национального исследовательского университета).

Научный руководитель

доктор технических наук, профессор Агульник Алексей Борисович

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

- доктор технических наук, профессор Колотников Михаил Ефимович

- кандидат технических наук, доцент Белоусов Григорий Геннадьевич

Ведущая организация

ОАО «Авиадвигатель»

Защита состоится 14 ноября 2011 г. в 13 часов на заседании диссертационного совета Д 212.125.08 при Московском авиационном институте (национальном исследовательском университете) по адресу: 125993, г.Москва, А-80, ГСП-3, Волоколамское шоссе, д.4.

С диссертацией можно ознакомиться в библиотеке Московского авиационного института (национального исследовательского университета)

Автореферат разослан «_¿_» ¿^¡¿^Z^C-__ 2011 г.

Ученый секретарь диссертационного

совета Д 212.125.08, профессор, д.т.н.

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

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

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

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

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

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

Флаттер лопаток тесно связан с нестационарными течениями в турбомашинах. Но задачи о неустановившемся обтекании колеблющихся лопаток турбомашин существенно отличаются от аналогичных задач для одиночного крыла. Это объясняется следующими причинами:

• существенно большей жесткостью лопаток по сравнению с жесткостью крыла (это приводит к тому, что частота флаттера незначительно отличается от частот собственных колебаний лопатки, в то время, как у крыла это отличие существенно);

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

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

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

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

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

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

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

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

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

Основные задачи исследования.

1. Исследование влияния вязкости потока на расчет условий возникновения флаттера в компрессорной решетке профилей энергетическим методом;

2. Исследование возможности возникновения изгибно-крутильного и сверхзвукового безотрывного флаттера в рабочем колесе вентилятора с помощью энергетического метода;

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

4. Исследование течения воздуха в системе «воздухозаборник -рабочее колесо вентилятора»;

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

Метод исследования.

При численном моделировании задач нестационарной газовой динамики применялся разработанный в ЦИАМ ПК Cobra, позволяющий определять стационарные и нестационарные аэродинамические характеристики элементов силовых установок и летательных аппаратов. Метод расчета основывается на неявном модифицированном методе С.К.Годунова и позволяет использовать подвижные сетки. Расчет течений в колеблющихся решетках проводился в слое переменной толщины.

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

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

Основные результаты исследования, выносимые на защиту.

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

• определять условия возникновения аэроупругой неустойчивости и флаттера решеток компрессоров в вязкой и невязкой постановках;

• анализировать свободные колебания профилей в вязком нестационарном потоке воздуха;

• исследовать обтекание вязким потоком воздуха систем двигателя, например, «воздухозаборник - рабочее колесо вентилятора», и определять нестационарные аэродинамические нагрузки на лопатки рабочего колеса вентилятора;

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

Научная новизна.

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

2. Исследована решетка свободно колеблющихся профилей в вязком нестационарном потоке воздуха (совместная задача) на двух режимах обтекания:

а) На режиме, близком к границе газодинамической устойчивости, получены условия возникновения вращающегося срыва;

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

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

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

Обоснованность и достоверность результатов исследования.

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

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

Практическая значимость. Разработанная методика позволяет моделировать и анализировать динамическое поведение лопаток компрессоров в нестационарном потоке воздуха с учетом и без учета вязкости на обычном персональном компьютере без больших затрат вычислительных и временных ресурсов.

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

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

Апробация работы. Основные результаты диссертационный работы докладывались и обсуждались на научно-практическом семинаре «Теория, численные методы и математический эксперимент в газовой динамике», посвященная юбилею А.Н. Крайко (Москва, 2009); XII международном семинаре «Супервычисления и математическое моделирование» (Саров, 2010); Международной конференции MSC.Software (Москва, 2010), III Международной Научно-Технической Конференции, посвященной 80-летию ЦИАМ им. П.И.Баранова (Москва, 2010), на научно-технических совещаниях ММПП «Салют» (Москва, 2011), Международном семинаре «Проблемы авиационного двигателестроения» (Москва, 2011), на заседаниях кафедры «Теория ВРД» МАИ.

Публикации. По материалам диссертации опубликовано 8 печатных работ, в том числе 2 статьи в изданиях, рекомендованных ВАК, 3 материала докладов на Всероссийских и Международных научно-технических конференциях. Выпущено 5 отчетов о НИР.

Струюура и объем работы. Диссертация состоит из введения, четырех глав, основных выводов, списка литературы (56 наименований). Основная часть работы содержит 132 страницы, 59 иллюстраций, 8 таблиц.

ОСНОВНОЕ СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

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

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

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

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

3. Расчеты обтекания колеблющихся решеток и рабочих колес без учета взаимного влияния конструкции и потока;

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

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

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

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

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

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

При использовании этого метода исследования аэроупругой неустойчивости решеток турбомашин, движение профилей определяется строго заданным гармоническим законом движения. Вывод о флаттере делается на основе знака коэффициента аэродинамического демпфирования D = -W /ка2,

1/<а„ jа

где W= [ М—di - аэродинамическая работа за цикл колебаний, М -о dt

аэродинамический момент, а и ú)0 - амплитуда и частота колебаний лопатки соответственно. Отрицательный коэффициент D соответствует положительной работе потока над профилем и говорит о флаттере, положительный D — о демпфировании колебаний потоком.

В ходе работы по созданию и отработке методики было выполнено несколько расчетов.

1. Расчет аэроупругих характеристик решетки профилей в вязком дозвуковом потоке.

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

сдвига фаз 6 = 2к -п!Nя (здесь п = 1, 2, З...и т.д., Nп - число лопаток) были получены значения коэффициента аэродемпфирования с учетом и без учета вязкости, представленные на рисунке 2. Сравнение полученных значений коэффициента аэродемпфирования В с экспериментальными данными дает удовлетворительные результаты, а игнорирование вязкости существенно влияет

Рисунок 1. Исследуемая решетка профилей

Рисунок 2. Коэффициент аэродинамического демпфирования в зависимости от межлопаточного угла сдвига фаз 5

2. Расчетная оценка возможности возникновения сверхзвукового безотрывного флаттера в рабочем колесе вентилятора.

Целью данной работы было исследование изгибно-крутильных колебаний лопатки и оценка возможности возникновения флаттера. Объектом исследования была лопатка перспективного КНД. При расчете лопатки на собственные частоты и формы колебаний было получено, что 2-я и 3-я частоты близки. Вторая форма колебаний является крутильной, а третья - изгибной. Поэтому, было принято решение исследовать изгибно-крутильные колебания данной лопатки на 2-й собственной частоте. Явление равенства частот колебаний невозможно в реальных конструкциях. Но именно для такого исследования изгибно-крутильных колебаний лопатки было принято это допущение.

Для исследования поведения решетки выбран профиль, находящийся на среднем радиусе лопатки. Колебания лопатки по изгибной форме для выбранного профиля соответствуют перемещению И. Колебания по крутильной форме - повороту вокруг центра профиля на угол а. Таким образом, профили совершали поступательные перемещения И и поворот на угол а (рисунок 3).

смещенный профиль

профиль О

Рисунок 3. Система обозначений в исследуемой решетке

Расчеты проводились в вязкой (Unsteady Reynolds Averaged Navier-Stocks equations - URANS) и невязкой (Euler) постановках.

В соответствии с применяемым для исследования энергетическим методом анализа флаттера, в задаче моделировались гармонические колебания профилей, т.е. уравнения движения каждого профиля имели вид: h = Hh sin a>0t для перемещения профиля и а = На sin coQt для вращения. Здесь Hh,Ha -начальные амплитуды, которые принимаются равными единице, со0 -частота собственных колебаний лопатки по выбранной форме, t - время.

Параметрами данной задачи являлись:

• сдвиг фаз между перемещениями соседних профилей, соответствующими изгибным колебаниям лопатки - 8и;

• сдвиг фаз между вращательными движениями соседних профилей, соответствующими крутильным колебаниям - 8к (в данной задаче принимаем, что 8и = 8к);

• сдвиг фаз между поступательным и вращательным движениями профиля, соответствующий сдвигу фаз между изгибными и крутильными колебаниями лопатки) - А.

В результате последовательной «переборки» сдвигов фаз между перемещением и вращением соседних профилей 8t=Sk и сдвига фаз между поступательным и вращательным движениями профиля А, было получено, что при значениях 8и =8к =-160° и А = -40° к профилям решетки подводится максимальная работа, т.е. наблюдаются условия для возникновения флаттера (рисунок 4а,б). На рисунке 4а также приведена точка, соответствующая значению работы аэродинамических сил над профилем, полученная в расчете без учета вязкости (Euler).

а) в зависимости от 5и = 8к при б) в зависимости от А при

фиксированном А = -40° фиксированных <5„ = дк = -160°

Рисунок 4. График изменения работы над профилями

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

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

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

В результате расчетов удлиненной решетки только на режиме, соответствующем т.9 на рисунке 5, и только при одном значении угла сдвига фаз между изгибными колебаниями соседних профилей 8 = 20° наблюдается подвод работы. Но полученный набор близких к нулю значений работы дает надежду полагать, что расчеты проводятся вблизи границы флаттера (рисунок 6).

Рисунок 5. Напорная ветка, ппр= 1,0

Рисунок 6. Изменение суммарной работы за период изгибных колебаний в зависимости от угла сдвига фаз между изгибными колебаниями соседних

профилей

3. Решение совместной задачи (аэродинамика-колебания) на двух режимах обтекания.

Исследования проводятся в решетке, полученной сечением цилиндрической поверхностью с относительным (отнесенным к радиусу вентилятора) радиусом Я = 0.95 рабочего колеса модельного вентилятора.

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

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

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

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

На рисунке 8 для наглядности сопоставляются Фурье-разложения момента аэродинамических сил, действующего на 0-ой профиль, полученные при решении нестационарной задачи при неподвижных профилях (рисунок 8а) и при решении совместной задачи (рисунок 86). Здесь 5 - сдвиг фаз между соседними профилями, N - количество профилей в решетке.

Из представленных зависимостей видно, что в совместной задаче амплитудно-частотные характеристики момента изменились таким образом, что его максимальный период колебаний стал кратен (7) периоду собственных колебаний лопатки. Частоты остальных гармоник оказались кратными этой минимальной частоте (7/2, 7/3, 7/4 и т.д.), а амплитуды гармоник 7/1 и 7/2 возросли по сравнению с их аналогами примерно в 2 раза.

и.

245 255 265 275

Л

151

245 255 265 275

Рисунок 7. Напорная ветка, ппр= 1,05

8 = 2ллм

/

5 = 4лЛ \

1(|1 1

а) неподвижные профили б) свободно колеблющиеся профили

(аэродинамическая задача) (совместная задача)

Рисунок 8. Фурье-разложение момента аэродинамических сил, действующего на 0-ой профиль решетки для задач разного типа

Т.е. колебания потока «подстраиваются» под колебания лопаток по первой собственной форме. При этом две зоны срыва, образованные течением в решетке с неподвижными профилями, постепенно превращаются в одну большую зону срыва, движущуюся против вращения решетки со скоростью, составляющей 75% от скорости вращения решетки. Такое явление можно классифицировать как вращающийся срыв.

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

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

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

Исследованная конфигурация, состоящая из упрощенного осесимметричного

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

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

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

Изменение момента при пересечении лопаткой зоны максимальной окружной неравномерности в периферийной зоне установившегося течения (с 23 по 24 оборот) достигает 7% от среднего значения этого параметра, тогда как в неустановившемся течении значение момента достигает 30% от среднего значения этого параметра.

Рисунок 10. Временная зависимость момента аэродинамических сил действующего на лопатки рабочего колеса вентилятора за 23 оборота венца

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

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

Рисунок 9. Исследуемая конфигурация

течения в системе «воздухозаборник - рабочее колесо вентилятора», позволяет сделать выводы о наличии или отсутствии «раскачки» лопатки, что дает возможность судить о степени ее устойчивости. Также, рассчитав напряженно-деформированное состояние (НДС) лопатки, можно сделать выводы о ее статической прочности и ресурсе.

Для определения вынужденных колебаний лопатки используется традиционный подход, основанный на решении уравнения вынужденных колебаний:

[м]{«}+[см+№Мпо}, (1)

где: {u(t)\ - вектор узловых перемещений, [М], [С], [К] - матрицы масс, демпфирования и жесткости, - вектор узловых аэродинамических

сил.

Решением уравнения (1) является искомое движение узловых точек лопатки, которое представляется в виде разложения по ее собственным формам с коэффициентами, зависящими от времени:

{и(0} = 1>у(0{ру->, (2)

7=1

где yj(t)- коэффициенты разложения (главные координаты), {<pj) - векторы

собственных форм колебаний.

С учетом (2), линейной связи между матрицей демпфирования [С] и матрицами массы [М] и жесткости [К], а также свойством ортогональности собственных форм, уравнение (1) можно свести к системе обыкновенных дифференциальных уравнений для коэффициентов разложения y^t):

У jit) + 2 a>jZj?j(t) + ojjyj(t) = Fj(t), (3)

где: Fj(t) = {q>j}r{F(t)} - проекция аэродинамической силы на вектор

перемещений по j-ой собственной форме; coj- частота колебаний лопатки по j-й

собственных форме, - коэффициент механического демпфирования. В

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

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

- проводится расчет трехмерной лопатки на собственные частоты и формы колебаний в любом пакете конечно-элементного анализа (например, ПК Patran/Nastran);

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

- для каждой собственной формы определяется функция времени Fj(t) и

решается обыкновенной дифференциальное уравнение (3) с начальными условиями yj(0) = 0, y'j(O) = 0;

- после определения всего набора коэффициентов разложения у^) для

j=l,..., п по формуле (2) определяются перемещения всех узлов конечно-элементной механической сетки {«(г)};

- по известным коэффициентам разложения у^) для ]=1,..., п по

п

формуле {<т(<)} = X Л(0} определяются напряжения во всех узлах конечно-/=1

элементной механической сетки. Здесь {сг,} - вектор напряжений в узлах лопатки при ее колебаниях по собственным формам, который получается, как и вектор перемещений по собственным формам колебаний, из ПК Ра^ап/^авКап.

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

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

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

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

• определять условия возникновения аэроупругой неустойчивости и флаттера решеток компрессоров в вязкой и невязкой постановках;

• анализировать свободные колебания профилей в вязком нестационарном потоке воздуха;

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

В ходе работ по созданию методики были решены следующие задачи:

1. Задача поиска условий возникновения флаттера в компрессорных решетках профилей. В результате решения выявлена важность учета влияния вязкости потока;

2. Задача исследования решетки свободно колеблющихся профилей в потоке воздуха на двух режимах:

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

периоду собственных колебаний лопатки, что может привести к вращающемуся срыву;

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

3. Задача обтекания воздухом системы «воздухозаборник - рабочее колесо вентилятора». Такое исследование позволяет уточнить картину течения и образования отрыва в воздухозаборнике, что позволяет более точно оценить нестационарные нагрузки на лопатки, возникающие из-за их взаимодействия с отрывом в воздухозаборнике;

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

ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

Публикации в изданиях, рекомендованных ВАК:

1. Буюкли Т.В., Шорстов В.А. К вопросу о постановке задачи расчета колебаний решеток компрессоров для определения аэроупругих характеристик [Электронный ресурс] // Труды МАИ. 2010. №39. URL: http://www.mai.ru/science/trudy/

2. Буюкли Т.В., Шорстов В.А. Расчет обтекания колеблющихся профилей неявным методом для оценки аэроупругой устойчивости решетки в вязком потоке // Известия вузов. Авиационная техника. 2011. №1. С. 9-12.

Прочие публикации:

3. Макаров В.Е., Шорстов В.А., Федорченко Ю.П., Андреев С.П., Буюкли Т.В., Фролов В.Н.. Развитие модифицированного метода С.К. Годунова для решения нестационарных вязких 2D/3D задач аэродинамики // Теория, численные методы и математический эксперимент в газовой динамике: Сборник тезисов научно-практического семинара, ЦИАМ, 25-26 августа 2009 г. - Москва, ЦИАМ, 2009, С.86-87.

4. Макаров В.Е., Андреев С.П., Шорстов В.А., Буюкли Т.В.. Аэроупругость лопаточных машин: 3D моделирование вынужденных колебаний лопаток вентилятора и расчетная оценка условий

возникновения автоколебаний в квази-ЗБ постановке. // Супервычисления и математическое моделирование: Сборник тезисов докладов XII Международного семинара, РФЯЦ-ВНИИЭФ, 11-15 октября 2010г.- Саров, РФЯЦ-ВНИИЭФ, 2010. С.92.

5. Макаров В.Е., Андреев С.П., Шорстов В.А., Буюкли Т.В.. Исследование условий возникновения автоколебаний в вентиляторе ТРДЦ для БСМС // Основные результаты научно-технической деятельности ЦИАМ (2009-2010) / под общ.ред. Скибина В.А., Солонина В.И., Гуревича О.С. - М., ЦИАМ, 2010. С.285-289. - ISBN 978-5-94049-027-2.

6. Макаров В.Е., Андреев С.П., Шорстов В.А., Буюкли Т.В., Федорченко Ю.П., Фролов В.Н.. Оценка нестационарных нагрузок на лопатки вентилятора ТРДД для БСМС при стационарной неравномерности на входе в двигатель // Основные результаты научно-технической деятельности ЦИАМ (2009-2010) / под общ.ред. Скибина В.А., Солонина В.И., Гуревича О.С.. - М„ изд. ЦИАМ, 2010. С.289-293. -ISBN 978-5-94049-027-2.

7. Макаров В.Е., Андреев С.П., Шорстов В.А., Буюкли Т.В. Аэроупругость лопаточных машин: 3D моделирование вынужденных колебаний лопаток вентилятора и расчетная оценка условий возникновения автоколебаний в квази-ЗБ постановке [Электронный ресурс] // Авиадвигатели XXI века: материалы конф. - Электрон.дан. -М.: ЦИАМ, 2010. - 1 электрон.опт.диск (CD-ROM). - Системные требования: IBM PC, Windows 2000 или выше. - Загл.с этикетки диска. - ISBN 978-5-94049-026-5.

8. Шорстов В.А., Буюкли Т.В., Макаров В.Е., Андреев С.П. Аэроупругость лопаток вентилятора: 3D моделирование вынужденных колебаний и расчетная оценка условий возникновения автоколебаний в квази-ЗО постановке // Супервычисления и математическое моделирование. Труды XII международного семинара / Под ред. Р. М. Шагалиева. - Саров: ФГУП «РФЯЦ-ВНИИЭФ», 2011. - 418с.:ил.

Подписано в печать 28.09.2011 Формат А5

Бумага для цифровых печатных машин 80 г/кв.м Печать цифровая. Тираж 100 экз.

Типография «Принт-сервис» 125481, г. Москва, ул. Туристская, д. 27, корп. 1 Тел. 8 (495) 767-7824,496-0737,8 (965) 306-1653 Ьир://принтсервис.сот, e-mail: tci.print@gmail.coni

Оглавление автор диссертации — кандидата технических наук Буюкли, Татьяна Васильевна

Условные обозначения.

ВВЕДЕНИЕ.

Цель работы.

Метод исследования.

Научная новизна.

Обоснованность и достоверность результатов исследования.

Практическая значимость.

Основные результаты исследования, выносимые на защиту.

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

ГЛАВА 1. ДИНАМИЧЕСКОЕ ПОВЕДЕНИЕ ЛОПАТОК В ПОТОКЕ ВОЗДУХА: АЭРОУПРУГАЯ НЕУСТОЙЧИВОСТЬ И ФЛАТТЕР.

1.1. Динамическое поведение лопаток в потоке воздуха.

1.2. Разновидности флаттера.

1.3. Методы анализа флаттера.

Выводы по главе 1.

ГЛАВА 2. ОСНОВНЫЕ УРАВНЕНИЯ.

2.1. Нестационарная аэродинамика.

2.1.1. Осредненные по Рейнольдсу уравнения Навье-Стокса.

2.1.2. Конечноразностная схема.

2.1.3. Особенности метода решения нестационарных задач.

2.1.4. Использование подвижных сеток.

2.1.5. Используемые модели турбулентности.

2.1.5.1. Однопараметрическая модель турбулентности Спаларта-Аллмараса.

2.1.5.2. Двухпараметрическая модель турбулентности k-s.

2.1.5.3. Двухпараметрическая модель турбулентности SST.

2.1.6. Метод расчета обтекания конструкций в ПК Cobra.

2.1.7. Определение нестационарных характеристик обтекания изолированного профиля.

2.2. Механическая задача.

2.3. Алгоритм исследования динамического поведения решеток профилей и лопаток компрессоров в потоке воздуха.

Выводы по главе 2.

ГЛАВА 3. ИССЛЕДОВАНИЕ ДИНАМИЧЕСКГО ПОВЕДЕНИЯ ПРОФИЛЕЙ ЛОПАТОК КОМПРЕССОРОВ.

3.1. Анализ динамического поведения профилей в вязком дозвуковом потоке.

3.1.1. Особенности постановки задачи.

3.1.2. Экспериментальные исследования.

3.1.3. Расчет колебаний профилей.66*

3.1.4. Результаты расчетов колебаний и их обработка.

3.2. Расчетная оценка возможности возникновения сверхзвукового безотрывного флаттера в рабочем колесе вентилятора.

3.2.1. Геометрическая модель лопатки.

3.2.2. Собственные частоты и формы колебаний лопатки.

3.2.3. Расчет стационарного течения.

3.2.4. Расчет нестационарного течения.

3.2.5. Результаты исследования.

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

3.3.1. Получение стационарного решения.

3.3.2. Поиск условий возникновения аэроупругой неустойчивости.

3.3.3. Расчет аэроупругой неустойчивости в решетке профилей с использованием измененной расчетной области.

3.3.3.1. Постановка задачи.

3.3.3.2. Результаты исследования аэроупругой неустойчивости в удлиненной расчетной области.

3.4. Решение задачи о свободных колебаниях профилей в вязком нестационарном потоке.

3.4.1. Объект исследования и постановка задач обтекания и колебаний

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

3.4.2.1. Расчет обтекания без учета механических колебаний.

3.4.2.2. Решение совместной задачи (обтекание с учетом механических колебаний).

3.4.3. Решение совместной задачи на режиме запирания.

Выводы по главе 3.

ГЛАВА 4. ЗАДАЧА О ВЫНУЖДЕННЫХ КОЛЕБАНИЯХ ЛОПАТКИ, НАХОДЯЩЕЙСЯ ПОД ДЕЙСТВИЕМ НЕСТАЦИОНАРНОЙ АЭРОДИНАМИЧЕСКОЙ НАГРУЗКИ.

4.1. Расчетная оценка нестационарных нагрузок на лопатки вентилятора перспективного ТРДД при стационарной неравномерности на входе в двигатель.

4.1.1. Исследуемая геометрия.

4.1.2. Постановка задачи.

4.1.3. Результаты расчетов.

4.2. Решение задачи о вынужденных колебаниях лопатки.

4.2.1. Решение задачи о вынужденных колебаниях лопатки в матричном виде с использованием главных координат.

4.2.1.1. Теоретическое обоснование решаемой задачи.

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

4.2.2. Задача о вынужденных колебаниях в пакете конечно-элементного анализа.

4.2.3. Расчет напряженно-деформированного состояния лопатки.

Выводы по главе 4.

ВЫВОДЫ.

Введение 2011 год, диссертация по авиационной и ракетно-космической технике, Буюкли, Татьяна Васильевна

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

Сложность моделирования динамического поведения лопаток компрессоров »заключается в комплексном характере этой проблемы: в ней I тесно переплетены вопросы нестационарной аэродинамики (т.к. колебания лопаток во множестве случаев связаны с переменными аэродинамическими силами), задачи колебаний тонких закрученных лопаток, находящихся в упругом взаимодействии с другими элементами ротора, вопросы усталостной прочности материала [1];

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

С явлением флаттера конструкторы авиационных ГТД впервые столкнулись в середине 50-х годов при создании реактивных двигателей второго поколения [2]. Возникновение флаттера связано с тем, что при улучшении аэродинамических качеств компрессорных ступеней, получаемом увеличением окружных скоростей и степени повышения полного давления, разработчики уменьшают толщины профилей лопаток, замковых соединений и дисков.

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

В1 литературе отмечены случаи разрушения лопаток от флаттера, которые приводили к серьезным последствиям для двигателей и самолетов. На рисунке 1 в качестве примера представлен характер разрушений в двигателе при усталостном разрушении от флаттера двух лопаток вентилятора [2].

Рисунок 1 - Схема расположения мест разрушения узлов при флаттере лопаток вентилятора [2]: 1 - разрыв защитного фланца корпуса вентилятора; 2 - разрушение по елочному замку 2-х лопаток и трещины еще на 2-х лопатках; 3 -разрушение, перехлест и задиры по антивибрационным полкам; 4 -разрывы корпуса вентилятора; 5 - смещение вперед носка КСД; 6 - срез заклепочного соединения корпуса с силовым корпусом; 7 - срез шпилек крепления корпуса КСД с силовым корпусом; 8 - обрыв 2-х лопаток КСД; 9 - разрушение вала КСД

Флаттер лопаток тесно связан с нестационарными течениями в турбомашинах и может быть объяснен только при изучении неустановившихся течений [3].

Задачи о неустановившемся обтекании колеблющихся лопаток турбомашин существенно отличаются от аналогичных задач для одиночного крыла. Это объясняется следующими причинами:

• существенно большей^ жесткостью лопаток по сравнению с жесткостью крыла (это приводит к тому, что частота флаттера незначительно отличается от частот собственных колебаний лопатки, в то время, как у крыла это отличие существенно);

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

• аэродинамической и механической связанностью лопаток в колесе (на лопатку действует спектр аэродинамических и упругих сил, зависящих от колебаний других лопаток) [2].

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

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

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

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

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

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

Цель работы

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

Метод исследования

При численном моделировании задач нестационарной газовой динамики применялся разработанный в ЦИАМ ПК Cobra, позволяющий определять стационарные и нестационарные1 аэродинамические характеристики элементов силовых установок и летательных аппаратов. Метод расчета основывается на неявном модифицированном методе С.К.Годунова и позволяет использовать подвижные сетки. Расчет течений в колеблющихся решетках проводился в слое переменной толщины.

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

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

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

1. 3D расчет лопатки на собственные частоты и формы колебаний;

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

3. построение расчетной модели решетки профилей и задание граничных условий в ПК Cobra;

4. получение стационарного обтекания решетки;

5. определение закона движения для профиля на основе расчета лопатки на собственные формы и частоты (точный закон движения для расчета флаттера энергетическим методом и «осциллятор» - для совместной задачи);

6. анализ обтекания колеблющихся профилей (для задачи флаттера на всех межлопаточных углах сдвига фаз 8 = Im/Njj, где п = 1,2,3., Ыл - число лопаток);

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

8. для задачи флаттера: вывод о возможности возникновения аэроупругой неустойчивости и флаттера на основе полученных данных о работе аэродинамического потока над профилями.

Научная новизна

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

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

3. Исследована, решетка свободно колеблющихся профилей* в вязком нестационарном потоке воздуха (совместная задача) на двух режимах обтекания.

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

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

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

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

Обоснованность и достоверность результатов исследования

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

Достоверность результатов решения задачи о вынужденных колебаниях подтверждается сравнением с коммерческими ПК, основанными на МКЭ.

Практическая значимость

Разработанная методика позволяет моделировать и анализировать динамическое поведение лопаток компрессоров в нестационарном потоке воздуха с учетом и без учета вязкости на обычном ПК без больших затрат вычислительных и временных ресурсов.

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

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

Основные результаты исследования, выносимые на защиту

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

• определять условия возникновения аэроупругой неустойчивости и флаттера решеток компрессоров в вязкой и невязкой постановках;

• анализировать свободные колебания профилей в вязком нестационарном потоке воздуха;

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

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

Разработанная методика используется в ФГУП ЦИАМ им. П.И. Баранова для исследования динамического поведения профилей и лопаток компрессоров в нестационарном потоке воздуха с учетом и без учета вязкости.

Заключение диссертация на тему "Моделирование динамического поведения лопаток компрессоров авиационных двигателей в нестационарном потоке воздуха"

Выводы по главе 3

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

2. При исследовании изгибно-крутильных колебаний лопатки в рассмотренных условиях обтекания решетки профилей в рамках используемого подхода и расчетной методики, получено, что при значениях сдвигов фаз между перемещением и вращением, соседних профилей», соответствующих сдвигу фаз между изгибными и крутильными* колебаниями соседних лопаток, 6и-8к- -160° и сдвига фаз между поступательным-№ вращательным движениями профиля; А = -40° к лопаткам подводится максимальная работа, т.е. имеют место условия для возникновения сверхзвукового безотрывного > флаттера.

3. При исследовании • чисто изгибных колебаний получено, что большинство изгибных колебаний демпфируются потоком. А возникновение аэроупругой неустойчивости возможно только в сильно отрывном течении.

4. С использованием разработанной методики проведены исследования структуры' нестационарных возмущений- и анализ свободно колеблющихся профилей в решетке, полученной сечением рабочего колеса модельной ступени вентилятора перспективного ТРДД, на 2 режимах работы.

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

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

ГЛАВА 4. ЗАДАЧА О ВЫНУЖДЕННЫХ КОЛЕБАНИЯХ ЛОПАТКИ, НАХОДЯЩЕЙСЯ ПОД ДЕЙСТВИЕМ

НЕСТАЦИОНАРНОЙ АЭРОДИНАМИЧЕСКОЙ НАГРУЗКИ

4.1. Расчетная оценка нестационарных нагрузок на лопатки вентилятора перспективного ТРДД при стационарной неравномерности на входе в двигатель

В" экстремальных условиях полета (например, отказ одного из двигателей при взлете)гили в условиях старта при сильном боковом ветре и максимальном, расходе воздуха через двигатель, самолет вынужден переходить на углы атаки больше допустимых (а > 20°). При этом в воздухозаборнике на внутренней стороне наветренной части обечайки мотогондолы может возникать отрыв потока. Это происходит из-за того, что на стартовом режиме работы двигателя с расходами воздуха, близкими к максимальным, и при малой- скорости набегающего потока (Мн < 0.1), обтекание относительно тонкой (для ТРДД большой степени двухконтурности) обечайки воздухозаборника осуществляется с сильным разгоном потока на внутренней стороне кромки. Эта зона сверхзвукового течения, в зависимости от качества профилирования обечайки, может быть достаточно> велика, и замыкающий ее скачок уплотнения может вызывать отрыв пограничного слоя внутри воздухозаборника. Как известно, заметный боковой ветер ^>10 м/с) в условиях старта приводит к увеличению уровня скоростей на кромке и, как следствие, вероятность возникновения нерасчетного (отрывного) течения значительно возрастает. Этот отрыв может привести к возникновению зон окружной и радиальной неравномерностей потока перед рабочим колесом вентилятора, которые, в свою очередь, могут вызывать большие нестационарные нагрузки на его лопатки [54] .

Особенность постановки задачи заключалась в следующем:

- течение во внутренней области рассчитывалось в системе координат, вращающейся вместе с рабочим колесом;

- течение в остальной части расчетной области рассчитывалось в неподвижной системе координат;

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

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

В качестве граничных ставились следующие условия:

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

- лопатка — вращающаяся вязкая стенка;

- внутренняя поверхности обечайки — вязкая стенка, вращающая в противоположном направлении;

- сечение за рабочим колесом вентилятора - заданное статическое давление с условием радиального равновесия;

- внешняя граница — условия в набегающем потоке.

Используемая расчетная сетка, иллюстрируемая на рисунке 4.5, содержала около 4 млн. ячеек.

110000

105000

100000

95000

90000

85000

80000

75000

70000

65000

60000

55000

50000

45000

40000

35000

30000

25000 Р

Рисунок 4.6. Мгновенное поле статического давления и пристеночные линии тока на обтекаемых поверхностях после 23-го оборота

На рисунке 4.7 и рисунке 4.8 для того же момента времени (после 23-го оборота) показаны поля осевой компоненты скорости с пристеночными линиями тока в плоскости симметрии течения и в поперечном сечении воздухозаборника перед рабочим колесом вентилятора.

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

Рисунок 4.7. Мгновенное поле осевой скорости и пристеночные линии тока в плоскости симметрии течения после 23-го оборота систему «воздухозаборник — рабочее колесо вентилятора» даже при работе на углах атаки, больше допустимых.

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

4.2. Решение задачи о вынужденных колебаниях лопатки

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

Решение совместной задачи в п. 3.4. осуществлялось с учетом простейшей механической связи между колеблющимися профилями. Т.е. набора профилей представлялся, как ряд осцилляторов, колебания которых определяется уравнением т-1 + ск - 7 = Ра. Однако, такой способ задания механической связи'позволяет решать совместную задачу только в квази-ЗО-постановке.

Решение трехмерной совместной задачи сопряжено с множеством трудностей, таких как огромные затраты расчетного времени или затруднения создания алгоритма- расчета движения трехмерных сеток. Поэтому решение совместной задачи в 3-0 постановке в большинстве случаев можно заменить решением задачи о вынужденных колебаниях конструкции под действием газовых сил. Такая постановка задачи позволяет сделать выводы о наличии или отсутствии' «раскачки» конструкции, что позволяет судить о степени ее устойчивости. Также, рассчитав напряженно-деформированное состояние (НДС) конструкции, можно судить о ее статической прочности и ресурсе.

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

4.2.1. Решение задачи о вынужденных колебаниях лопатки в матричном виде с использованием главных координат

4.2.1.1. Теоретическое обоснование решаемой задачи

Рассмотрим уравнение свободных колебаний системы с конечным числом степеней свободы в матричном виде [48] :

МЙ + [1]{М} = {0} (4.1)

Где [М] - матрица масс, [К] - матрица жёсткости, {и} - вектор перемещений.

Решение этого уравнения ищем в виде: м(0} = {^}зт(©/) (4.2)

Подставив [2] в [1] и проведя преобразования, получаем:

-^2ям]+т){^}={0} (4.з)

Решение такого уравнения существует, только если его определитель равен нулю, т.е.

К]-<о][М]\ = 0 (4.4)

Из решения уравнения (4.4) получаем собственные частоты системы сосо]. Частота в Гц вычисляется по формуле: ^ = --

PjY[cm) = см = 1 v

4.14)

Т.к. {<Pj}r[C]{(Pi) = ZgjODj

Матрицу жесткости, согласно свойствам ортогональности, можно представить в виде:

Р;}Г1КШ} = К^ (4.15)

С учетом (4.11) получаем: = со*

Используя полученные преобразования матриц жесткости, масс и демпфирования в векторы, имеем систему уравнений, для каждой собственной формы cpj: у, (О + 2 а&у, (0 + а>)у, (0 = FJ (0, (4.16) где Fj(i) - {(Pj}T{F{t)} - аэродинамическая сила.

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

Задача о вынужденных колебаниях лопатки рабочего колеса под воздействием заданной аэродинамической силы решается в следующей последовательности: проводится расчет трехмерной лопатки на собственные частоты и формы колебаний в любом пакете конечно-элементного анализа (например, ПК Patran/Nastran); осуществляется пересчет распределения давления на поверхности лопатки, полученного из расчета обтекания конструкции) с узлов аэродинамической сетки в узлы механической сетки;

- для каждой собственной формы определяется функция времени Fj(t) и решается обыкновенной дифференциальное уравнение (4.16) с начальными условиями уу(0) = 0, ]йу(0) = 0. Здесь у- - перемещение каждого узла по ]-й собственной форме, а у^ - скорость каждого узла, перемещающегося по ]-й собственной форме.

- после определения всего набора коэффициентов разложения у^ для ]=1,., п по формуле (4.6) определяются перемещения всех узлов конечно-элементной механической сетки (и(/)};

- по известным коэффициентам разложения у^) для ]=1,., п по п формуле {сг(/)} = X! -У/ (0 } > гДе {с,-} - вектор напряжений в узлах лопатки

1=1 при ее колебаниях по собственным формам (получается, как и вектор перемещений по формам колебаний, из ПК Ра^ап/ЫаБ^ап) определяются напряжения во всех узлах конечно-элементной механической сетки (сг(/)}.

Ниже иллюстрируются перечисленные выше этапы решения поставленной задачи.

Для основного исследуемого режима с ппр = 1 (пвр = 17520 об/мин) вычислены 10 собственных частот колебаний лопатки, представленных в таблице 4.2.

Библиография Буюкли, Татьяна Васильевна, диссертация по теме Тепловые, электроракетные двигатели и энергоустановки летательных аппаратов

1. Писаренко Г.С., Ольштейн JI.E. Проблемы аэроупругости лопаток турбомашин // Проблемы прочности, №86, 1974 С. 3-8.

2. Хориков А. А. Обеспечение отсутствия флаттера лопаток компрессоров на различных этапах создания турбомашин / Проблемы прочности, 1976, №3. С. 25-28.

3. Самойлович F.C. Нестационарное обтекание и- аэроупругие колебания.решеток турбомашин. М.:.издательство.«Наука», 1969. 444 с.

4. J.G. Marshall; М: Imregun A- Review of Aeroelasticity Methods with, Emphasis on Turbomachinery Applications // Journal' of 'Fluids and Structures, Vol.10, №.3, 1996. pp. 237-267.

5. A.V. Shrinivasan. Flutter and. Resonant Vibration Characteristics of Engine Blades // Journal, of Engineering- for Gas Turbines and Power, Vol.119, №3, 1997. pp. 742-775.

6. J.M. Verdon Review of Unsteady Aerodynamic Methods for Turbomachinery Aeroelastic and Aeroacoustic Applications // AIAA Journal Vol.31, №2, 1993. pp. 235-249.

7. Фершинг Г. Основы аэроупругости. М.: Машиностроение, 1984.600 с.

8. Ольштейн JI-1E. Новые аспекты проблемы аэроупругости турбомашин // Проблемы прочности, №3, 1976. С. 3-7.

9. Н. Doi. Fluid/Structure Coupled Aeroelastic Computations For Transonic Flows In Turbomachinery. Ph.D. Thesis Stanford University, Department of Aeronautics - and Astronautics, 2002.

10. H. Thermann, R. Niehuis. Unsteady Navier-Stokes Simulation of a Transonic Flutter Cascade Near Stall Conditions Applying Algebraic Transition Models // Proceedings of ASME Turbo Expo 2005, June 6-9, Reno-Tahoe, Nevada, USA.

11. R. Srivastava, M.A. Bakhle, T.G. Keith, G.L. Stefko. Flutter Analysis of a Transonic Fan // Proceedings of ASME Turbo Expo 2002, June 3-6, Amsterdam, The Netherlands.

12. Y.N. Shmotin, R.Y. Starkov, P.V. Chupin, V.I. Gnesin, L.V. Kolodyazhnaya. Numerical Solution of Aeroelastic Behavior of Isolated Fan Bladerow // Turbomachines: Aeroelasticity, Aeroacoustics, and Unsteady

13. Aerodynamics / Edited by V.A. Skibin, V.E. Saren, N.M. Savin, S.M. Frolov. Moscow: Torus Press, 2006. pp. 24-26.

14. Bendiksen Recent Developments in Flutter Suppression Techniques for Turbomachinery Rotors // Journal of Propulsion and Power, Vol.4, №2, 1988. pp. 164-172.

15. P. Petrie-Repar, A. McGhee, P. Jacobs, R. Gollan Analytical Maps of Aerodynamic Damping as a Function* of Operating Condition for a Compressor Profile // Proceedings of ASME Turbo Expo 2006, May 8-11, Barcelona, Spain.

16. Снайдер, Коммерфорд. Сверхзвуковой, безотрывный флаттер в роторах вентиляторов: теория и эксперимент. // Энергетические машины и установки (Пер. Trans ASME. Ser.A/1974. #4)

17. S. Fleeter. Aeroelasticity Research for Turbomachine Applications // Journal'of Aircraft, Vol.16, №5, 1979. pp. 320-342.

18. Борисов Г.А., Локштанов E.A., Олыптейн Л.Е. Вращающийся срыв в осевом компрессоре // Промышленная аэродинамика, вып.24. М.: Оборонгиз, 1962.

19. Кулагина'А.Т. Некоторые особенности колебаний лопаток в условиях вращающегося срыва // Проблемы прочности, №3, 1976: С. 45-48.

20. Хориков А.А. О влиянии близости собственных частот колебаний лопатки, по различным формам, на устойчивость однородного компрессорного венца к флаттеру // Проблемы прочности*№8, 1974. С. 83-87.

21. Хориков А.А. О возможности возникновения «классического» флаттера рабочих лопаток турбомашин // Проблемы прочности, №36 1976. С. 25-28:

22. О: Bendiksen, P. Freidmann. Coupled Bending-Torsion Flutter in Cascades // AIAA Journal, Vol.18, №2, 1980. pp. 194-201.

23. M. Montgomery, M. Tartibi, F. Eulitz, S. Shmitt. Application of Unsteady Aerodynamics and Aeroelasticity in Heavy-Duty Gas Turbines // Proceedings of ASME Turbo Expo 2005, June 6-9, Reno-Tahoe, Nevada, USA

24. Т.Н. Fransson. Analysis of Experimental Time-Dependent Blade Surface Pressures from an Oscillating Turbine Cascade Using the Influence-Coefficient Technique // Journal de Physique III-, Vol.2, №4, 1992. pp. 575-594.

25. Srivastava R., Bakhle M.A., Keith T.G and Stefko G.L.: Flutter analysis of a transonic fan. / National Aeronautics ans Space Administration, Washington, DC, 2002.

26. Reddy T.S.R., Bakhle M.A., Trudell J.J., Mehmed O., Stefko G.L. LINFLUX-AE: A Turbomachinery Aeroelastic Code Based on a 3-D'Linearized Euler Solver. / National Aeronautics ans Space Administration, Washington, DC, 20041

27. Reddy T.S.R., Srivastava-R., Mehmed O. Stefko G.L. Flutter and forced response analysis of cascades usinga two dimensional, linearized Euler Solver. / National Aeronautics ans* Space Administration; Washington, DC, 1999.

28. Srivastava R., Bakhle 1УГ.А., Keith T.G., Stefko G.L. Aeroelastic stability computations for turbomachinery. / National- Aeronautics ans Space Administration, Washington; DC, 200 к

29. M. Sadeghi, F. Liu. Computation of cascade flutter by uncoupled and' coupled methods < // International Journal of. Computational' Fluid Dynamics, Vol.19,' №8, 2005. pp.» 559-569.

30. I. McBean, K. Hourigan, M. Thompson, F. Liu. Prediction of Flutter of Turbine Blades in a Transonic Annular Cascade // Journal of Fluids Engineering; Vol. 1>27, №4'; 2005. pp. 1053-1058.

31. Годунов^ C.K., Забродин» A.B., Прокопов Г.П. Конечноразностная схема для двумерных нестационарных задач газовой, динамики и расчет-течения с отошедшей ударной волной. // Журнал вычисл. матем. и матем. физ. 1961, т.1, №6, стр.1020-1050:

32. Годунов С.К., Забродин А.В., Иванов М.Я., Крайко.А.Н., Прокопов» Г.П. Численное' решение многомерных задач газовой динамики. М., Наука; 1976, 400 стр.

33. Колган В1П. Применение принципа минимальных значений производных к построению конечноразностных схем разрывных решений газовой.динамики. // Уч. зат ЦАГИ, 1972, т.З, №6, стр.68-77.

34. Тилляева, Н.И: Обобщение модифицированной схемы С.К.Годунова на произвольные нерегулярные сетки. // Уч. зап. ЦАГИ, 1986, т. 17, №2, стр.25-33.

35. Родионов А.В. Монотонная схема второго порядка точности для скозного счета неравновесных течений. // Журнал вычисл. матем. и матем. физ., 1987, т. 27, № 4, стр.585-593.

36. Копченов В.И., Топеха Е.А. Неявная релаксационная конечно-разностная схема для системы уравнений Эйлера // Научно-технический отчет ЦИАМ №11543, 1989.

37. Гуляев А.Н., Козлов В.Е., Секундов. А.Н. К созданию универсальной однопараметрической модели» для турбулентной вязкости. // Изв. РАН, МЖГ, 1993, №2, стр.69.

38. Wilcox, D. С. Turbulence Modeling for CFD. DCW Industries, Inc., 5354 Palm Drive, La Cañada, Calif., 1993.

39. Spalart, P. R.; and Allmaras, S. R. A One-Equation Turbulence Model for Aerodynamic Flows. // AIAA Paper 92-0439, 1992.

40. Spalart, P. R.; and Allmaras, S. R. A One-Equation Turbulence Model for Aerodynamic Flows. // La Recherche Aerospatiale, no. 1, 1994, pp. 5-21.

41. Jones, W. P.; and Launder, В. E. The Predictiomof Laminarization with a Two-Equation Model of Turbulence. // International Journal of Heat and Mass Transfer, vol. 15, 1972, pp. 301-314.

42. Launder, В. E.; and Sharma, В. I. Application of the Energy Dissipation Model of Turbulence to the Calculation of Flow Near a Spinning Disc. // Letters in Heat and Mass Transfer, vol. I, no. 2, 1974, pp. 131-138.

43. Martelli, F., Belardini, E., Adami P. Unsteady Flow Simulation: A Numerical Challenge. //RTO-MP-89, France, 2002, pp.37-40.

44. Menter, F. R. Two-Equation Eddy Viscosity Turbulence Models for Engineering Applications. // AIAA J., vol: 32, Nov. 1994, pp. 1299-1310.

45. Menter, F. R. Influence of Freestream Values on k-e Turbulence Model Predictions. // AIAA J., vol. 30, no. 6, June 1992, pp. 1657-1659.

46. Menter, F. R. Zonal Two Equation k- 8 Turbulence Models for Aerodynamic Flows. // AIAA Paper 93-2906, Orlando, Fla., 1993.

47. W.J. McCroskey, L.W. Carr,| K.W. McAlister. Dynamic Stall Experiments on Oscillating Airfoils. // ALAA J., vol.14, no.l, January 1976, pp. 57-63.

48. М.И. Бабаков. Теория колебаний. Издательство «НАУКА», главная редакция*физико-математической литературы. М., 1968 г. 560 с.

49. Boles, Т.Н. Fransson. Aeroelasticity in turbomachines comparison of theoretical and experimental cascade results. // Communication du laboratoire de thermique appliqué et de turbomachines, Lausanne, EPFL, 1986.*

50. Буюкли T.B., Шорстов В.A. Расчет обтекания колеблющихся, профилей неявным методом для оценки аэроупругой устойчивости решетки в вязком потоке // Известия вузов. Авиационная техника. 2011. №1. С. 9-12.

51. Буюкли Т.В., Шорстов В.А. К вопросу о постановке задачи расчета колебаний решеток компрессоров, для определения аэроупругих характеристик Электронный ресурс. // Труды МАИ. 2010. №39. URL: http://www.mai.ru/science/trudy/

52. В.Н. Ершов. Неустойчивые режимы турбомашин. Вращающийся срыв. М.: издательство «Машиностроение», 1966, 180стр.

53. Численное моделирование нестационарных явлений в газотурбинных двигателях: Научное издание / Августинович В.Г., Шмотин Ю.Н. и др. М.: Машиностроение, 2005. - 536 с.