Авторефераты по всем темам >>
Авторефераты по техническим специальностям
ОАО УАвиадвигательФ На правах рукописи
Сипатов Алексей Матвеевич
Методология расчетного анализа нестационарных трехмерных процессов в авиационных двигателях
05.07.05 - Тепловые, электроракетные двигатели и энергоустановки летательных аппаратов
АВТОРЕФЕРАТ
диссертации на соискание ученой степени доктора технических наук
Пермь - 2010
Работа выполнена в ОАО УАвиадвигательФ, г. Пермь.
Научный консультант:
Заслуженный деятель науки РФ, доктор физико-математических наук, профессор, Соколкин Юрий Викторович, ПГТУ, г. Пермь.
Официальные оппоненты:
Академик РАН, Липанов Алексей Матвеевич, ИПМ УрО РАН, г. Ижевск Доктор физико-математических наук, Старший научный сотрудник, Осипов Анатолий Алексеевич, ЦИАМ, г. Москва Доктор физико-математических наук, профессор, Егоров Михаил Юрьевич, ПГТУ, г. Пермь
Ведущая организация: ОАО УСамарский научно-технический комплекс им.
Н.Д. КузнецоваФ
Защита состоится 20 декабря 2010г. в 14 ч. 00 мин. на заседании диссертационного совета Д 212.188.06 в Пермском государственном техническом университете по адресу: 614990 г. Пермь, Комсомольский пр., 29.
С диссертацией можно ознакомиться в библиотеке ПГТУ.
Автореферат разослан _________ 2010 г.
Ученый секретарь Диссертационного совета д.т.н., проф., В.И. Свирщев
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы В настоящее время в области создания и производства газотурбинных двигателей существует жесткая конкуренция между ведущими двигателестроительными фирмами во всем мире. Можно отметить две ключевые проблемы, успешное решение которых позволяет обеспечить необходимый объем заказов для дальнейшего существования фирмы и её развития: повышение конкурентноспособности разрабатываемых двигателей и соблюдение (желательно с существенным запасом) постоянно ужесточающихся международных норм по их экологическим характеристикам. В направлении повышения конкурентноспособности перед разработчиками и производителями авиационных и наземных газотурбинных двигателей ставятся довольно противоречивые требования: повышение надежности, повышение экономичности, увеличение ресурса работы, снижение массы, снижение себестоимости их производства. В рамках соблюдения экологических характеристик постоянно исследуются возможности дальнейшего снижения уровня шума двигателей и уменьшения выбросов вредных веществ, образующихся при сгорании топлива. С целью повышения точности расчетного анализа начинают переходить к моделированию многодисциплинарных процессов (multi-physics). Об этом свидетельствуют все последние научные конференции по двигателестроению. Современное развитие технологий проектирования достигло такого уровня, что дальнейшее их совершенствование невозможно без всестороннего исследования процессов, протекающих при работе газотурбинных двигателей. Такие исследования требуют применения последних научно-технических достижений в области методов численного анализа. Таким образом обязательной составляющей успешного решения отмеченных выше задач является применение в процессе проектирования методов численного анализа с учетом моделирования трехмерных и нестационарных явлений. Не менее важной является также задача постоянного совершенствования расчетных методов, что позволяет сохранить высокий научный уровень применяемых математических моделей.
Следует также отметить стратегический характер задачи повышения надежности и конкурентноспособности отечественных авиационных двигателей. Успешное решение данной задачи является важным фактором в развитии авиационного направления страны в целом.
Цель работы Целью настоящей работы являлась разработка методологии расчетного анализа, позволяющая значительно повысить уровень численного моделирования рассматриваемых в работе явлений, а именно:
- динамического нагружения лопаток турбомашин при статор-ротор взаимодействии;
- теплового состояния охлаждаемых лопаток турбин с учетом трехмерных и нестационарных процессов; применение сопряженной (газ - твердое тело) постановки;
- генерации и распространения акустических возмущений от вентиляторной ступени и реактивного сопла.
Таким образом в работе исследуются эффекты, возникающие при решении многодисциплинарных задач в области авиационного двигателестроения. В частности исследуются задачи по генерации акустических пульсаций турбулентным потоком, по взаимодействию этих пульсаций с элементами двигателя и их распространению в области ближнего и дальнего поля.
Отдельно исследуется вопрос оценки динамических напряжений в лопатках турбомашин. Представленная методология позволяет описать как нагружение исследуемых лопаток за счет действия нестационарных газовых сил, так и их динамический отклик. С точки зрения обеспечения длительного ресурса работы двигателя крайне важным является вопрос оценки теплового состояния лопаток турбины высокого давления (ТВД) с максимально возможным уровнем точности. В данном случае на основе наиболее критического элемента ТВД (рабочей лопатки первой ступени) представлена методология оценки теплового состояния лопаток в сопряженной постановке. Такая постановка позволяет полностью описать процесс теплообмена на границе газ - твердое тело без дополнительного задания коэффициентов теплоотдачи на поверхностях раздела сред.
Достоверность Достоверность разработанной методологии численного анализа подтверждается сопоставлением результатов расчетов с аналитическими решениями при рассмотрении модельных задач и сопоставлением с экспериментальными данными при исследовании натурных объектов.
Научная новизна В соответствии с отмеченными выше направлениями исследований следует отметить научную новизну полученных результатов:
- Оценка динамических напряжений в лопатках турбомашин.
Разработана математическая модель газодинамического статор-ротор взаимодействия (неотражающие граничные условия на входе и выходе из расчетной области, обобщенные условия периодичности, возможность выполнения расчетов с 4-м порядком точности аппроксимации исходных уравнений по пространству), позволяющая повысить точность численного моделирования газодинамической нагрузки на лопатки турбомашин до 50% по сравнению с разработанным ранее газодинамическим программным комплексом Pluton.
Разработана методика оценки коэффициента демпфирования лопаток турбомашин, которая позволяет сузить диапазон его возможного изменения примерно в 2-3 раза.
- Оценка теплового состояния охлаждаемых лопаток турбин.
Разработана методика оценки теплового состояния лопаток турбин на основе решения сопряженной задачи по теплообмену между газом и металлом лопатки в трехмерной постановке и выполнения нестационарного анализа.
Сопряженная постановка позволяет исключить необходимость задания коэффициентов теплоотдачи на обтекаемых поверхностях.
При выполнении нестационарного моделирования статор-ротор взаимодействия в исследуемой ступени впервые выполнена количественная оценка сегрегации газа в межлопаточном канале ротора с учетом трехмерных и вязких эффектов. Выявлен эффект эжекции горячего газа на спинку рабочей лопатки. Проведена его количественная оценка.
- Оценка акустического совершенства вентиляторной ступени авиационного двигателя.
Разработан эффективный расчетный метод оценки акустического и аэродинамического совершенства вентиляторной ступени авиационного двигателя, учитывающий трехмерные, вязкие и нестационарные эффекты газодинамического взаимодействия лопаток ротора и статора.
Разработана математическая модель процесса распространения акустических возмущений на основе решения линеаризованных уравнений Эйлера с применением схемы DRP, высокого порядка точности по пространству и времени с низкими диссипативными и дисперсионными ошибками. Сопоставление результатов расчетов с данными зарубежных авторов, полученных с применением аналогичных подходов, показывает, что при использовании разработанной программы получено более точное соответствие эксперименту. При этом уточнение результатов численного анализа достигает 5 дБ.
- Оценка шума, генерируемого реактивной струей авиационного двигателя.
Разработана методика численного моделирования процесса генерации шума при истечении реактивных струй, позволяющая корректно описывать генерацию шума до чисел Струхаля 3.5 - 4.0 с использованием коммерческого ПО ANSYS-Fluent. Эти результаты соответствует передовому уровню численного моделирования генерации шума от реактивных струй во всем мире.
Таким образом на защиту выносится:
- Математическая модель газодинамического статор-ротор взаимодействия, позволяющая повысить точность численного моделирования газодинамической нагрузки на лопатки турбомашин до 50% по сравнению с газодинамическим ПО ANSYS-CFX, ANSYS-Fluent.
- Методика оценки коэффициента демпфирования лопаток турбомашин, которая позволяет сузить диапазон его возможного изменения примерно в 2-раза.
- Методика оценки теплового состояния лопаток турбин на основе решения сопряженной задачи по теплообмену между газом и металлом лопатки в трехмерной постановке и выполнения нестационарного анализа. С её помощью впервые выполнена количественная оценка сегрегации газа в межлопаточном канале ротора с учетом трехмерных и вязких эффектов.
Выявлен эффект эжекции горячего газа на спинку рабочей лопатки. Проведена его количественная оценка.
- Расчетный метод оценки акустического и аэродинамического совершенства вентиляторной ступени авиационного двигателя, учитывающий трехмерные, вязкие и нестационарные эффекты газодинамического взаимодействия лопаток ротора и статора.
- Математическая модель процесса распространения акустических возмущений на основе решения линеаризованных уравнений Эйлера с применением схемы DRP. Сопоставление результатов расчетов с данными зарубежных авторов показывает, что при использовании разработанной программы получено более точное соответствие эксперименту. При этом уточнение результатов численного анализа достигает 5 дБ.
- Методика численного моделирования процесса генерации шума при истечении реактивных струй, позволяющая корректно описывать генерацию шума до чисел Струхаля 3.5 - 4.0 с использованием коммерческого ПО ANSYS-Fluent.
Практическая значимость В результате проведенных исследований были получены следующие практические результаты:
- Разработан газодинамический программный комплекс Euler_3D, позволяющий моделировать нестационарную газодинамическую нагрузку на лопатки турбомашин.
- Разработан программный комплекс LEEAxi, позволяющий моделировать распространение акустических возмущений в осесимметричной постановке.
- Выполнено исследование различных вариантов геометрии спрямляющего аппарата вентиляторной ступени двигателя ПД-14 с точки зрения минимизации акустических источников статор-ротор взаимодействия.
По результатам расчетов спроектирован профиль лопатки спрямляющего аппарата (СА) вентиляторной ступени, который позволил снизить уровень тонального шума вентиляторной ступени в источнике на 4 дБ.
- Выполнена оценка динамических напряжений для различных вариантов лопаток рабочего колеса и направляющего аппарата первой ступени компрессора высокого давления двигателя ПД-14. По результатам расчетов определена конструкция лопатки рабочего колеса, позволяющая снизить максимальный уровень расчетных динамических напряжений в 6 раз по сравнению с исходной геометрией лопатки.
- Проведена серия исследований теплового состояния ряда рабочих лопаток ТВД для различных двигателей в совместной (газ - твердое тело) постановке. В частности выполнено расчетное исследование теплового состояния рабочей лопатки первой ступени ТВД двигателя ПД-14. По результатам расчетов намечено дальнейшее совершенствование системы перфорации лопатки.
Апробация работы Основные результаты работ были представлены на различных конференциях и семинарах. В частности можно отметить следующие конференции:
2000г.: международная научно-техническая конференция Двигатели века, Москва - декабрь. 2002г.: научно-техническая сессия по проблемам газовых турбин, Москва - сентябрь. 2004г.: международная конференция по газотурбинным двигателям ASME-2004, Вена, Австрия - июнь. 2005г.:
международная конференция по газотурбинным двигателям ASME-2005, Рино, штат Невада, США - июнь; международная конференция по авиационным двигателям ISABE-2005, Мюнхен, Германия - сентябрь; семинар пользователей Fluent, Санкт-Петербург - октябрь; международная научно-техническая конференция Авиадвигатели XXI века, Москва - декабрь. 2006г.:
Международная конференция УВычислительный эксперимент в аэроакустике 2006Ф, Светлогорск - сентябрь; XVIII сессия Российского акустического общества, Таганрог - сентябрь. 2007г.: всероссийская научно-техническая конференция УАэрокосмическая техника, высокие технологии и инновацииФ, ПГТУ, Пермь - апрель; международная конференция по газотурбинным двигателям ASME-2007, Монреаль, Канада - Май; международная конференция по авиационным двигателям ISABE-2007, Пекин, Китай - сентябрь; международная конференция УНовые рубежи авиационной наукиФ, Москва, в рамках МАКС-2007 - август; XIX сессия Российского акустического общества, Нижний Новгород - сентябрь. 2008г.: Международная конференция УВычислительный эксперимент в аэроакустике 2008Ф, Светлогорск - сентябрь;
Присуждение премии Лауреата премии Пермского края в области науки I степени за лучшую работу в области информационных технологий и вычислительных систем за серию работ УВиртуальная среда для проектирования турбомашинФ, Пермь - февраль. 2009г.: Семинары: институт механики сплошных сред УрО РАН, г. Пермь - ноябрь; ООО УК ОДК, г.
Москва - октябрь; ЦИАМ, г. Москва - август, ноябрь; на фирме Pratt&Whitney, Хартфорд, штат Коннектикут, США - апрель. Конференции: международная конференция по газотурбинным двигателям ASME-2009, Орландо, штат Флорида, США - июнь; всероссийская конференция по авиационной акустике, Москва - октябрь; международная конференция по авиационным двигателям ISABE, Канада, Монреаль - сентябрь. 2010г.: Научно-технический конгресс по двигателестроению (НТКД2010), Москва - апрель; международная конференция по аэроакустике (AIAA/CEAS), Стокгольм, Швеция - июнь.
Публикации Основные положения и результаты исследований отражены в 36 научных работах. Из них 10 работ напечатаны в журналах, рекомендованных ВАК для опубликования научных результатов докторских диссертаций, 18 работ - в зарубежных журналах и трудах международных конференций, 2 работы - в трудах сессий Российского акустического общества. Две работы изданы в виде монографий.
Структура диссертации Диссертация состоит из введения, пяти глав, заключения и списка литературы.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обосновывается актуальность рассматриваемых задач, указываются поставленные цели, отмечается научная новизна проведенных исследований и разработанных математических моделей, а также приводится краткая характеристика работы в целом.
Первая глава В этой главе представлен обзор работ по численному моделированию газодинамических процессов применительно к авиационным двигателям.
Данный обзор сделан по научным направлениям, которые будут детально рассматриваться в последующих главах.
Развитие расчетных схем и их применение.
В настоящее время процесс развития расчетных схем во всем мире направлен на повышение их точности и устойчивости. Он включает в себя совершенствование дисперсионных и диссипативных свойств применяемых разностных схем, а также повышение их устойчивости, прежде всего, при описании скачков уплотнения в потоке. Такое развитие схем существенно расширяет диапазон газодинамических задач, которые могут быть решены численно с высоким уровнем точности. Низкая схемная вязкость позволяет описывать процесс распространения акустических возмущений и их генерацию в турбулентном потоке, а устойчивое и корректное описание скачков уплотнения в сочетании с низкой схемной вязкостью позволяет, в частности, более корректно смоделировать сильное взаимодействие скачка уплотнения с пограничным слоем на лопатках турбомашин.
В нашей стране в этом направлении следует отметить работы, выполняемые в институте математического моделирования РАН под руководством член-корр. РАН Б.Н. Четверушкина и к.ф.-м.н. Т.К. Козубской. В данном случае проводится развитие разностных схем высокого порядка точности и их адаптация на неструктурированные сеточные модели. В ЦИАМ (г. Москва) для решения задачи генерации тонального шума при ротор-статор взаимодействии применяется разработанная ранее разностная схема Годунова - Колгана - Родионова второго порядка аппроксимации по времени и пространству. Необходимая точность при этом достигается за счет более мелкой дискретизации расчетной области. Примеры применения такой постановки можно найти в работах д.ф.-м.н. А.А. Осипова. Для моделирования процессов распространения акустических пульсаций широко применяются численные алгоритмы на основе схемы высокого порядка точности DRP (Dispersion Relation Preserving). Применение такого подхода представлено в работах к.т.н. В.И. Милешина. Достаточно интересные результаты были получены с применением схемы Кабаре. Её развитием занимаются в настоящее время д.ф.-м.н. В.М. Головизнин (Московский институт безопасного развития атомной энергетики РАН) и д.ф.-м.н. С.А. Карабасов (технический факультет Кембриджского университета, Великобритания). Сопоставление различных расчетных схем по решению одномерной тестовой задачи (распространение импульса Гаусса малой амплитуды) свидетельствует о её низких диссипативных и дисперсионных свойствах. Следует также отметить применяемую в настоящей работе модифицированную схему ЛаксаВендроффа. Её развитие представлено в работах, выполненных в институте прикладной механики УрО РАН (г. Ижевск) под руководством академика РАН А.М Липанова. В данном случае гашение нефизичных пульсаций, вызванных применением центрально-разностного шаблона, выполняется за счет введения схемной вязкости и дополнительной монотонизации решения. Как показывают одномерные тесты, выполненные при помощи данной схемы с четвертым порядком точности по пространству, она несколько уступает по точности схеме DRP, но вместе с тем, обладает большей устойчивостью при описании скачков уплотнения, что является важным при моделировании течения в турбомашинах.
Помимо указанных выше отмечу научно-исследовательские центры, внесшие значительный вклад в развитие численных методов газовой динамики в нашей стране. Это научные коллективы МГУ (А.А. Самарский, академик РАН, д.ф.-м.н., Г.Г. Черный), институт теоретической и прикладной механики СО РАН (академик РАН, д.ф.-м.н., В.М. Фомин, д.ф.-м.н., проф., М.С. Иванов), институт математики им. С.Л. Соболева СО РАН (академик РАН, д.ф.-м.н., С.К.
Годунов), институт прикладной математики им. М. В. Келдыша, ЦИАМ им.
П.И. Баранова (д.ф.-м.н., профессор, А.Н. Крайко) и ряд других научных центров, указанных в первой главе в диссертации.
В США в 90-х годах XX века по инициативе NASA была создана рабочая группа для разработки нового метода численного решения законов сохранения - метода пространственно-временных элементов (CE/SE). Основное отличие метода заключается в том, что запись интегральных уравнений сохранения выполняется 4-мерном евклидовом пространстве, включающем пространственные и временные координаты. Сопоставление данного метода с коммерческим ПО и методом DRP представлено во второй главе настоящей работы.
Введение специальной системы базисных функций для решения газодинамических задач методом Галеркина позволяет, как свидетельствуют публикации, достаточно эффективно выполнять численное моделирование рассматриваемых процессов. Среди достоинств такого подхода авторами отмечается: устойчивость численного алгоритма, низкие диссипативные свойства, обеспечение заданного порядка аппроксимации, задание неотражающих граничных условий. Учитывая относительную новизну применения данного подхода к решению газодинамических задач, как правило, в открытой печати представлены примеры решения одно- и двумерных задач.
Публикации, в которых бы рассматривались примеры решения трехмерных нелинейных задач газовой динамики в области авиационного двигателестроения с использованием такого подхода, практически отсутствуют.
Достаточно широкое применение в процессе исследования нестационарных газодинамических задач получил метод крупных частиц. Суть метода заключается в разделении этапов численного интегрирования уравнений газовой динамики на каждом шаге по времени на два этапа: Эйлеров и Лагранжев (разработчики метода - академик РАН О.М. Белоцерковский, академик РАЕН Ю.М. Давыдов) В настоящее время под их руководством продолжаются работы по совершенствованию схемы с целью улучшения её диссипативных и дисперсионных свойств. В г. Перми её развитием занимается д.ф.-м.н., профессор, М.Ю. Егоров.
Высокий уровень получаемых результатов, применительно к моделированию генерации шума реактивными струями, показывают работы д.ф.-м.н. М.Х. Стрельца и к.ф.-м.н. М.Л. Шура (СПбГТУ). Применяемая математическая модель базируется на решении нестационарных уравнений Навье-Стокса неявной конечно-объемной схемой повышенного порядка аппроксимации. Благодаря ряду оригинальных приемов, математическая модель позволяет на относительно грубых сетках (2Ц5 миллионов узлов) получать приемлемый уровень точности в частотном спектре до числа Струхаля порядка 2-3, а в ряде случаев, и для более высоких частот.
Новое направление, принципиально отличное от рассмотренных выше направлений развития газодинамического анализа, связанно с применением уравнения Больцмана вместо системы уравнений Навье-Стокса. Прямая аппроксимация данного уравнения в семимерном пространстве (пространство, импульсы молекул и время) приводит к недоступным на сегодняшний день требованиям по вычислительным ресурсам. В связи с этим используются различные аппроксимации исходного уравнения, позволяющие снизить его размерность. В настоящее время наиболее активно данное направление развивается к.ф.-м.н. И. Ибрагимовым (фирма Elegant Mathematics, Германия).
В результате проведенного анализа разностных схем, можно отметить, что, на сегодняшний день, не существует Усамого оптимальногоФ расчетного метода, позволяющего решать весь спектр задач газовой динамики, который наблюдается при моделировании процессов в авиационных двигателях. Для каждого класса рассматриваемых задач существует свой набор расчетных схем, позволяющих наиболее адекватно описать исследуемые явления и процессы.
Именно исходя из этих позиций, были выбраны применяемые в последующих главах математические постановки, которые позволили корректно описать исследуемые процессы.
Моделирование теплового состояния лопаток турбин Как хорошо известно, имеется много различных факторов, влияющих на тепловое состояние лопаток турбин. Прежде всего требуется учесть влияние трехмерных и нестационарных эффектов течения газа. Не менее важным является выбор области, где выполняется численный анализ, а также корректное определение на её границах граничных условий. В целом опыт моделирования свидетельствует, что только совместная (CHT - Conjugate Heat Transfer) постановка позволяет адекватно описать процесс выдува охлаждающего воздуха через отверстия перфорации и связанный с ним конвективный теплообмен газа с поверхностью лопатки.
В случае применения совместной постановки многие исследователи пытаются учесть эффект влияния нестационарности течения на тепловое состояние лопатки. В настоящей работе такой учет выполнен с помощью гибридной постановки. Отдельно моделировалось нестационарное газодинамическое статор-ротор взаимодействие, оценивалась величина влияния нестационарных эффектов на адиабатическую температуру стенки и на этой основе выполнялась оценка её влияния на температуру металла лопатки. В этом направлении следует отметить работу He (Оксфордский университет, 2009). В данной работе проблема различия временных масштабов тепловой и газодинамической задачи (примерно на 4 порядка) решалась с помощью перевода нестационарной тепловой задачи в частотную область (в предположении её линейности). При этом полученное решение не являлось точным, однако данная аппроксимация предлагает путь возможного одновременного решения тепловой и газодинамической задачи в сопряженной постановке с учетом нестационарных эффектов.
Оценка динамических напряжений в лопатках турбомашин В работе для оценки динамических напряжений применяется расчетная методика, состоящая из следующих двух этапов:
1. Расчет газодинамической задачи (например статор-ротор взаимодействия) и определение нестационарной газодинамической нагрузки на лопатки исследуемого венца.
2. Решение нестационарной прочностной задачи, где в качестве граничных условий используется нестационарная газодинамическая нагрузка на лопатки, полученная при выполнении первого этапа. В предположении, что вынужденные колебания лопатки могут быть представлены через сумму её колебаний по собственным формам и, предполагая линейность рассматриваемой задачи, она сводится к определению собственных форм и частот колебаний лопатки, а также нахождению модальных координат через систему несвязанных между собой обыкновенных дифференциальных уравнений для исследуемых форм колебаний. В этом случае демпфирование лопатки описывается с помощью коэффициентов демпфирования по рассматриваемым собственным формам.
Описанные выше этапы представляют собой базовую постановку по расчетной оценке динамических напряжений. Одной из серьезных трудностей, возникающей при применении данного подхода, является выбор коэффициента демпфирования для исследуемой лопатки по заданной форме колебаний. Для его корректного определения требуется учесть много различных факторов: тип крепления лопатки к диску, номер формы, имеющиеся центробежные нагрузки, связанность лопатки и диска, неидентичность лопаток и др. Анализ имеющихся экспериментальных данных позволил сделать некоторые обобщения, которые существенно ограничивают область возможного изменения коэффициентов демпфирования для каждой исследуемой формы колебаний лопатки. В частности в результате анализа имеющихся экспериментальных данных по величине логарифмического декремента, найденных автором в открытой печати, удалось построить упрощенную зависимость (штриховая линия) суммарного уровня демпфирования от номера резонансной формы вида d=10%/n.
В частности, на рисунке 1 можно выделить три группы форм колебаний:
первая группа это формы колебаний с 1-ой до 6Е7-ой (среди них - подгруппа изгибных форм). Для этой группы в суммарный уровень демпфирования лопатки существенный вклад вносит как конструкционное демпфирование, так и аэродемпфирование. Во второй группе форм колебаний (с 7-ой по 13-ю формы - крутильно-изгибные) вклад конструкционного демпфирования незначительный и основную роль играет аэродемпфирование. В третьей группе форм (с 14-ой формы и далее - пластинчатые) основную роль играет материальное демпфирование. Соответственно, с увеличением номера формы уменьшаются демпфирующие свойства лопаток.
- Расчетная оценка Авиадвигатель, - Экспериментальная оценка Авиадвигатель, - Зарубежные публикации Рис. 1 - Упрощенная зависимость логарифм. декремента от номера формы Моделирование генерации и распространения акустических возмущений В целом оценить последние достижения в области вычислительной аэроакустики можно по материалам ежегодных международных конференций AIAA-CEAS, ASME Turbo-EXPO, AIAA ISABE и УВычислительный эксперимент в аэроакустикеФ. Кроме того, следует отметить международный журнал по аэроакустике (International journal of aeroacoustics ISSN 1475 472 X), где публикуются наиболее интересные разработки. По результатам анализа материалов отмеченных выше конференций можно выявить некоторые тенденции и особенности современного состояния по численному моделированию акустических процессов, а также понять эффективность их применения для решения конкретных задач авиационного двигателестроения.
В случае применения разностных схем высокого порядка точности (четвертого и выше) для моделирования процесса генерации акустических возмущений при истечении реактивной струи появляется возможность проводить расчет одного варианта конструкции в течение одной-двух недель.
Это связано с тем, что такие схемы имеют существенно меньшую схемную диссипацию и дисперсию. Необходимая сеточная размерность при этом составляет порядка 6-10 млн. узлов. Отмечается также важность корректного моделирования режима течения в пограничном слое вблизи среза сопла. В работе Eastwood S. (Кембриджский университет, 2008) показано, что в зависимости от режима течения в пограничном слое, расчетные данные по уровню шума в дальнем поле могут отличаться на 5 дБ! Одним из исследуемых в настоящее время направлений по уменьшению шума реактивных струй является применение специально профилированного выходного участка сопла (установка шевронов). В настоящее время эффекты воздействия шевронов на характеристики турбулентного истечения струи всесторонне исследуются. В этой связи интересно отметить ряд работ по моделированию обтекания крылового профиля и взаимодействия следа от стержня с крыловым профилем. При этом по второй задаче проводилось сравнение 5-ти некоммерческих газодинамических пакетов по их возможности описывать генерацию акустических возмущений. Как следует из полученных результатов, ни один из пакетов не дал полного согласования с имеющимися экспериментальными данными. Это свидетельствует о научной привлекательности данной тематики в плане дальнейшего развития применяемых математических моделей.
Генерация шума в процессе газодинамического статор-ротор взаимодействия лопаток турбомашин как правило выполняется на основе гибридного метода. Так для моделирования процесса нелинейного ротор-статор взаимодействия и образования акустических источников применяется трехмерная газодинамическая программа. Далее, с помощью акустического решателя (на основе линеаризованных уравнений Эйлера) проводится расчет распространения акустических возмущений в области ближнего поля и, затем, проводится расчет шума в области дальнего поля на основе интегрального решения волнового уравнения для неограниченного пространства.
Глава Как правило, характерной чертой газодинамического программного обеспечения собственной разработки является его специализация на решение определенного класса задач. Такая специализация позволяет выбрать математическую постановку, наиболее полно отражающую все происходящие процессы. Вместе с тем, при разработке такой математической модели становится возможным ввести определенные допущения, которые не вносят заметный вклад в получаемое решение, но существенно упрощают её программную реализацию. Представленная в настоящей главе разработка специализированного ПО не является исключением. Класс решаемых задач был ограничен моделированием газодинамического нестационарного статор-ротор взаимодействия между соседними венцами лопаток турбомашин. Цель такого моделирования - получение нестационарных динамических нагрузок на лопатки со стороны газового потока для дальнейшего расчета в них динамических напряжений. Для решения данного класса задач была выбрана математическая модель, которая отражает все характерные процессы рассматриваемого явления и не является чрезмерно громоздкой. В качестве исходной системы уравнений была выбана система нестационарных трехмерных уравнений Эйлера. Выбор конечно разностной схемы был обусловлен следующим моментом. Как правило, для решения стандартных газодинамических задач, считается достаточным применение схем второго порядка точности по пространству и времени. Вместе с тем, для корректного описания процесса распространение акустических возмущений, вызванных статор-ротор взаимодействием, вопрос выбора конечно-разностной схемы становится крайне важным. Так при моделировании процесса распространения акустических пульсаций четвертый порядок конечно-разностной схемы является минимальным для получения результатов с приемлемым уровнем точности. В соответствии с этим выбор конечно-разностной схемы был обусловлен возможностью при необходимости увеличить её порядок точности без принципиальной реорганизации численной процедуры решения. Это позволит в дальнейшем решать, в том числе, и класс акустических задач, связанных со статор-ротор взаимодействием. В результате была выбрана явная двухшаговая схема Лакса-Вендроффа варианта Рихтмайера. При этом вся расчетная область дискретизируется ортогональными прямолинейными ячейками. В этом случае интерфейс между статором и ротором отсутствует.
Вместе с тем, возникает проблема описания структуры течения в пограничном слое. Однако цель создания газодинамического программного комплекса - получение нестационарных газодинамических нагрузок на лопатки турбомашин при статор-ротор взаимодействии. В этом случае оказывается возможным пренебречь описанием структуры течения в пограничном слое и моделировать стенки расчетной области как абсолютно скользкие.
Представим теперь применяемую математическую постановку. Исходную систему уравнений можно записать в следующем виде:
Wt + AWx + BWy + CWz = 0 (1) здесьW ={, u, v, w, E} - вектор неизвестных, A, B, C - матрицы Якоби Используя расщепление по пространственным направлениям для одномерного случая применяемую схему можно записать в следующем виде (3D случай рассмотрен в диссертации):
n Wi+1/ 2 = Li+1/ 2(Wi n ) n+1/ n n Wi+1/ 22 = Wi+1/ 2 - t / 2A(Wi+1/ 2)Ri+1/ 2(Win )x-1 (2) n+1/ n+ Win+1 = Wi n - t(Fi+1/ 22 - Fi-1/1/ 2)x-1 - t(Din - Din )x-2 +1/ 2 -1/ Здесь F - потоковый вектор в направлении i, D - член, описывающий искусственную вязкость, а величины L, R, S аппроксимируют вектор неизвестных W требуемым образом. Точность схемы при этом определяется количеством членов ряда в выражениях (3):
N / m+1 n n Li+1/ 2(Win ) = (-1) am (Wi+m +Wi-m+1) m=N / m+1 n n Ri+1/ 2(Wi n ) = (3) (-1) bm (Wi+m -Wi-m+1) m=N / m+1 n n Si+1/ 2(Win ) = (-1) dm (Wi+m -Wi-m+1) m=Di+1/ 2 = i+1/ 2Si+1/ 2(Win ), Коэффициенты am, bm, dm выбираются таким образом, чтобы обеспечить заданный порядок аппроксимации (N) исходной системы уравнений.
Кроме того, при выполнении газодинамических расчетов с 4-м порядком точности схемы дополнительно применялась монотонизация решения с помощью нелинейной коррекции конвективных потоков (FCT) на промежуточном (n+1/2) слое по времени.
Верификация разработанного решателя была выполнена как на одномерных тестах: задача Сода и распространение импульса Гаусса малой амплитуды в однородном потоке, так и на трехмерной тестовой задаче:
обтекание колеблющихся турбинных лопаток газовым потоком (одиннадцатая стандартная конфигурация, Одномерные тесты показали высокий уровень результатов, получаемых с помощью применяемой схемы Лакса-Вендроффа. Остановимся более подробно на полученных результатах при обтекании турбинной лопатки.
Экспериментальные работы по данной лопатке были выполнены в Шведском Рис. 2 - Экспериментальная модель а) корыто вх. кромка спинка б) корыто вх. кромка спинка Рис. 3 - Распределение амплитуды (а) и фазы (б) нестационарного давления вдоль профиля лопатки в среднем сечении (сдвиг фаз - 180 градусов).
Безразмерная амплитуда Фаза, (град.) королевском институте. Вид лопаток представлен на рис.2. Экспериментальные исследования были проведены для заданной формы (первая изгибная) и амплитуды колебаний лопатки. Также были зафиксированы параметры газового потока на входе и выходе. По результатам выполненных исследований определены безразмерная амплитуда и фаза нестационарного давления на первой гармонике вдоль профиля исследуемой лопатки в среднем сечении по высоте для различных сдвигов фаз в колебаниях соседних лопаток. На рис. представлено сопоставление расчетных и экспериментальных данных по распределению амплитуды и фазы нестационарного давления для случая, когда сдвиг фаз между колебаниями соседних лопаток составляет 180 градусов. При этом расчетные данные получены как с использованием разработанного автором программного комплекса Euler_3D, так и с помощью программного комплекса Pluton (совместная разработка ОАО Авиадвигатель и ПГТУ, 1997г.). Как видно, применение Euler_3D позволяет повысить точность моделирования амплитудно-фазовых характеристик газодинамического нагружения лопатки до 50%. Особое внимание при этом следует обратить на повышение точности по распределению фазы нестационарного давления вдоль профиля, поскольку именно фазовые характеристики определяют возможность возникновения флаттера. Следует также отметить, что в рамках разработанной программы Euler_3D реализованы обобщенные условия периодичности, которые позволяют применять условия периодичности при некратном соотношении лопаток статора и ротора, а также неотражающие граничные условия на входе и выходе.
Глава В настоящее время методика оценки вибронапряжений, как правило, включает в себя газодинамическое моделирование статор-ротор взаимодействия в двумерной или трехмерной постановке, вычисление собственных форм колебаний лопатки, вычисление модальной нагрузки, а также выполнение корректной оценки коэффициента демпфирования для заданной формы колебаний. Иногда вместо задания коэффициента демпфирования решается отдельная контактная задача по моделированию характеристик демпфера. Такой подход является более точным.
В данной главе представлены результаты работ, последовательно проводимые по этой тематике. В настоящее время для определения нестационарных аэродинамических нагрузок на рабочих лопатках используются как квази-трехмерная так и полностью трехмерная постановки.
При оценке уровня вибронапряжений учитываются все основные факторы, влияющие на их уровень. Несмотря на это, учитывая сложность рассматриваемой проблемы, основной целью данного исследования ставилась разработка методики оценки максимального уровня возможных вибронапряжений в рассматриваемом венце, а не точное предсказание действительных значений вибронапряжений в анализируемой лопатке. Кроме того, было уделено внимание исследованию влияния выбора газодинамического решателя на уровень получаемых результатов и описанию основных принципов задания коэффициентов демпфирования для заданной формы колебаний лопатки, уровня газодинамической нагрузки, а также наличия или отсутствия демпфера.
Представим ниже математическую модель по оценке динамических напряжений. Уравнение движения, используемое для определения перемещений и напряжений, записанное в конечно-элементной форме, можно записать в виде:
[M ]{u}+[C]{u}+[K]{u}= {F(t)} (4) && & где {u} - вектор узловых перемещений в узлах конечно-элементной модели, [M], [C], [K] - матрицы масс, демпфирования и жесткости, {F(t)} - вектор узловых сил.
Предполагая, что любая вынужденная или свободная динамическая реакция системы может быть представлена в виде взвешенной суммы мод её колебаний, можно записать:
n n {u(t)} = yi (t){i}, { (t)} = y (t){ } (5) j j i=1 j=В этом случае, выполняя ряд преобразований, и, делая некоторые разумные допущения, можно привести уравнение (4) к системе несвязанных между собой обыкновенных дифференциальных уравнений относительно модальных координат y (t) :
j...
y (t) + 2 y (t) + 2 y (t) = Fj (t) (6) j j j j j j где Fj (t) = { }T {F(t)} - обобщенная сила, возбуждающая j-тую собственную j форму, y (t) - отклик по j-той собственной форме; - собственная круговая j j частота j-той формы; { } - собственный вектор, вектор j-той собственной j формы; - коэффициент демпфирования по j-той собственной форме.
j Фактически, Fj (t) является проекцией распределенной нагрузки на вектор перемещений по j-той собственной форме.
Вибронапряжения в точке r определяются по формуле:
v (r) = (t[t1,t1+T ] (r,t) - minT ] (r,t) ) / 2 (7) max t[t1, t1+ где T - период изменения давления; (r,t) - какая-либо из компонент ~(r,t) тензора напряжений .
В рамках представленной методики оценки динамических напряжений в лопатках турбомашин, можно выделить следующие этапы численного анализа:
1. Определение собственных форм и частот лопатки с помощью пакета ANSYS.
2. Нахождение резонансных режимов работы двигателя по диаграмме Кэмпбелла.
3. Определение нестационарных газодинамических нагрузок на резонансных режимах на основе результатов расчета статор-ротор взаимодействия. На данном этапе использовалось как программное обеспечение собственной разработки Euler_3D, так и коммерческое ПО ANSYS-CFX.
4. Вычисление проекций возбуждающей силы на вектор перемещений по собственным формам (с помощью собственных подпрограмм к ПО ANSYS).
5. Определение коэффициентов демпфирования по собственным формам. Для их нахождения применяется методика собственной разработки (для предварительной оценки - рис.1). Более точная оценка требует анализа дополнительных факторов: уровень газодинамической и центробежной нагрузки, тип закрепления в замке, связанность диска и лопаток, неидентичность лопаток, величина аэродемпфирования и др. (методика оценки влияния этих факторов представлена в диссертации).
6. Вычисление уровня резонансных динамических напряжений с учетом полученных данных по проекциям возбуждающей силы и коэффициентов демпфирования (на основе подпрограмм собственной разработки).
7. Определение местных запасов по многоцикловому нагружению с учетом статической нагрузки по диаграмме Гудмэна.
Представлена верификация разработанной методики оценки динамических напряжений на основе рабочей лопатки двигателя ПС-90А2.
Экспериментальные данные по величине динамических напряжений в лопатках исследуемого колеса были получены в ходе испытаний газогенератора данного двигателя.
При оценке вибронапряжений в лопатках турбины большое значение имеет задача получения корректных граничных условий по ее нестационарному газодинамическому нагружению. В соответствии с этим было исследовано три различных газодинамических подхода по статор-ротор моделированию с точки зрения получения газодинамических нагрузок:
1. Трехмерный расчет с осреднением между венцами. На его основе выполнение расчетов в нестационарной квази-трехмерной постановке для трех сечений по высоте лопатки (втулка, среднее сечение, периферия). В данном случае для нестационарного моделирования применялась программа Pluton, разработанная на ОАО Авиадвигатель.
2. Полностью трехмерный нестационарный расчет с использованием коммерческого пакета CFX.
3. Полностью трехмерный нестационарный расчет с использованием собственного решателя Euler_3D.
Рис. 4 - Распределение числа Маха в среднем сечении расчетной области и вблизи лопаток В качестве режима работы двигателя, на котором было выполнено газодинамическое моделирование, выбран режим близкий к таким условиям работы двигателя, при которых имелся резонанс по 12-ой и 34-ой гармоникам.
На входе в расчетную область были заданы распределения полного давления, полной температуры и направление потока. Значение осредненного статического давления было использовано в качестве граничного условия на выходе. Для корректного описания нестационарного нагружения лопаток рабочего колеса по 12-ой гармонике, соответствующей количеству жаровых труб в двигателе, был выполнен расчет камеры сгорания. Распределение температуры в выходном сечении камеры сгорания было использовано для задания граничных условий на входе при моделировании статор-ротор взаимодействия. Распределение числа Маха в расчетной области представлено на рис. 4.
Как и ожидалось, в процессе испытаний был обнаружен резонанс по 1-ой изгибной (1И) и 1-ой крутильной (1К) формам, возбуждаемым 12-ой и 34-ой гармониками (количество жаровых труб и лопаток соплового аппарата соответственно).
На режиме 10100 об./мин. наблюдался резонанс по 1И форме. По результатам записи сигнала в тензодатчике и его Фурье обработки была построена спектральная характеристика сигнала. На основе огибающей линии f0.(красная линия на рис. 5) и следующей формулы: = 100% был 3 fres вычислен логарифмический декремент по 1И форме. Методика, на основе которой была получена данная формула, была разработана в ЦИАМ. Здесь fres - резонансная частота, f0.5 - ширина резонансного пика на уровне половины значения амплитуды мощности на резонансной частоте. Для исследуемого случая было получено, что логарифмический декремент примерно равен 12%.
Для 1К формы не удалось оценить уровень логарифмического декремента на основе аналогичного подхода. Это связано с низким измеренным уровнем динамических напряжений. Расчетная оценка свидетельствует о том, что логарифмический декремент будет находится в данном случае в диапазоне 23%.
Рис. 5 - Динамический отклик по 1И форме В результате применения описанного выше расчетного анализа было получено распределение динамических напряжений по 1К и 1И формам (рис.
6). Количественные результаты расчетов представлены в таблице 1. В данной таблице представлено сопоставление результатов расчета и экспериментальных данных. В частности приведены максимальные измеренные значения вибронапряжений. Данные по 1И форме представлены в соответствии со значением логарифмического декремента в 12%.
Рис. 6 - Распределение вибронапряжений по 1И (слева) и 1К (справа) формам Табл. 1. Сопоставление расчетных и экспериментальных данных CFD Частота Расчетная величина Эксп.
Номер вращ., виброанпряжений данные формы v v подход об./мин. , МПа (=12%) , МПа 1 (1И) 1 23.12-ая 11020 18.2 13.гармон.
=2.5% =2.0% (среднее (консерв.
значение оценка) ) 1 10.5 12.3 (1К) 34-ая 11372 9.2 8.6 10.гармон.
3 8.9 11.Существенно более низкий уровень напряжений был вычислен при использовании CFX решателя (второй CFD подход). Как можно видеть в этом случае уровень вибронапряжений на 42% ниже экспериментального. В то же время первый CFD подход завышает уровень вибронапряжений на 24%. Вместе с тем можно было бы ожидать, что мы получим уровень вибронапряжений в диапазоне значений на 0 - 20% ниже экспериментального. Это связано с тем, что для сравнения были взяты максимальные значения вибронапряжений, полученные в эксперименте. Эффект неидентичности лопаток мог привести к некоторому увеличению вибронапряжений, а данный эффект не был смоделирован в расчетах. Детальный анализ показывает, что наиболее вероятная причина полученных результатов связана с упрощенным моделированием выдува охлаждающего воздуха. В частности, во втором CFD подходе ряды отверстий перфорации были заменены меньшим количеством узких щелей. В то же время в первом CFD подходе выдув охлаждающего воздуха через лопатки соплового аппарата и рабочего колеса не моделировался вообще и заданный перепад давления моделировался за счет изменения высоты проточной части H в квази-3D решателе. Это привело к сохранению более резких градиентов температуры вблизи лопатки ротора идущих от камеры сгорания. Таким образом полученные результаты свидетельствуют о важности точного моделирования выдува охлаждающего воздуха для того, чтобы корректно описать газодинамическое нагружение по 12-ой гармонике. Данные по применению 3-го CFD подхода для получения динамических напряжений по 12-ой гармонике не представлены, поскольку для ускорения в нем расчетов в качестве расчетной области был выбран сектор меньший, чем периодичность по полю температур на входе в исследуемую ступень.
По 1К форме все представленные подходы показали хорошее согласование с экспериментальными данными. Таким образом в данном случае все подходы хорошо описали распространение следов от лопаток статора и их взаимодействие с лопатками ротора.
Интересно отметить, что консервативная оценка уровня вибронапряжений v по 1К форме находится в диапазоне 11.1 - 12.9 МПа. Этот уровень позволяет выполнить оценку вибронапряжений с некоторым запасом, который учитывает возможное неблагоприятное сочетание различных факторов для реальной геометрии турбины.
Представленная методика расчетного анализа была применена в процессе проектирования турбины и компрессора высокого давления перспективного авиационного двигателя ПД-14. В данном случае было показано, что расчетная оценка динамических напряжений в рабочих лопатках первой ступени ТВД находится в пределах допустимых значений. Вместе с тем для лопатки рабочего колеса первой ступени КВД был получен высокий уровень динамических напряжений по 23 форме колебаний, возбуждаемой лопатками направляющего аппарата, располагающимися за рабочим колесом. В результате проведенной совместно с ЦИАМ работы по доводке геометрии исследуемой ступени, удалось снизить максимальный уровень динамических напряжений в исследуемой лопатке в 6 раз.
Глава В данной главе представлены основные этапы разработанной методологии по расчету теплового состояния лопаток турбин с учётом трехмерных и нестационарных эффектов. Прежде всего был выполнен нестационарный расчет статор-ротор взаимодействия и исследованы эффекты сегрегации газа в межлопаточном канале ротора (первая часть исследований), далее была применена совместная (газ - твердое тело) постановка для получения данных по тепловому состоянию исследуемой лопатки (вторая стационарная часть). В процессе первого этапа работ были проанализированы нестационарные газодинамические эффекты. Дополнительно к этому, для оценки влияния температурной неоднородности на выходе из камеры сгорания на тепловое нагружение лопаток ротора, был выполнен расчет течения в камере сгорания данного двигателя с моделированием процессов горения. Во второй части проводилось расчетное исследование течения в ступени турбины с осреднением по углу параметров потока между статором и ротором. В данном случае использовалась существенно более подробная дискретизация расчетной области. Это позволило детально описать все отверстия перфорации, как на лопатках статора, так и на лопатках ротора, а также структуру течения в пограничном слое (Y+ в первом узле от стенки выдерживался в среднем в диапазоне 2-3). Кроме того, было детально смоделировано течение во внутренней полости рабочей лопатки. Одновременно с этим решалась тепловая задача в металле лопатки с учетом теплообмена между поверхностью лопатки и обтекающим её газом (CHT постановка).
В качестве объекта исследования выбрана рабочая лопатка первой ступени турбины высокого давления авиационного двигателя. Проточная часть турбины, геометрическая модель лопатки, а также расположение внутренних полостей и отверстий перфорации представлены на рис. 7. Как видно из приведенной схемы система охлаждения лопатки - конвективно-пленочная.
При этом шесть рядов перфорации расположено со стороны корыта и 3 ряда со стороны спинки. В задней полости лопатки были установлены штырьки, дополнительно турбулизирующие поток и увеличивающие поверхность теплообмена. Следует отметить, что для выбранной ступени ТВД имелись экспериментальные данные по тепловому состоянию рабочих лопаток. Эти данные были получены в ходе испытаний газогенератора данного двигателя.
Измерения выполнялись как с применением термопар, так и с использованием термокристаллов.
Рис. 7 - Проточная часть ТВД и модель исследуемой лопатки В результате была возможность выполнить идентификацию предлагаемой методики на основе сравнения поля температур на поверхности лопатки, полученного в результате численного моделирования с имеющимися экспериментальными данными. Как показывает такое сравнение, изложенная в настоящей главе методика оценки теплового состояния лопаток, позволяет выполнять процесс проектирования лопаток ТВД на качественно новом уровне.
Появляется возможность ещё на стадии проектирования выявлять наиболее опасные участки на лопатке с высокой степенью надежности и так же в режиме проектирования проводить оценку влияния мероприятий, направленных на улучшение её теплового состояния.
По результатам расчетов (рис. 8) видно, что входная температурная неоднородность потока практически полностью транслируется на лопатки рабочего колеса. Кроме того, на рис. 8 хорошо заметен эффект сегрегации газа в рабочем колесе: горячий газ располагается ближе к поверхности корыта лопатки, а холодный, наоборот, устремляется к спинке лопатки. Для того, чтобы данный эффект оценить количественно, было выполнено два расчета на одинаковой сеточной модели. Первый расчет выполнен в стационарной постановке, т.е. с осреднением параметров потока на интерфейсе между статором и ротором, а второй расчет - без осреднения. Результаты нестационарного расчета были осреднены по времени за период статор-ротор взаимодействия. В таком случае разность между осредненными параметрами потока и параметрами потока, полученными в результате стационарного расчета, будет характеризовать эффекты, вызванные нестационарным статорротор взаимодействием. Следует также отметить, что для данных двух расчетов все стенки в расчетной области предполагались адиабатическими. Таким образом анализировались следующие параметры:
T = Ttime _ averaged -Tsteady, V = Vtime _ averaged -Vsteady (8) Их распределения в двух сечениях по высоте лопатки представлено на рис. 9. Как можно видеть, более горячий газ (по сравнению со стационарным расчетом), обладающий более высоким уровнем скорости, находится со стороны корыта лопатки, а более холодная и медленная его часть располагается вблизи спинки. В соответствии с этим, разность вектора скорости вблизи корыта будет направлена по направлению потока, а вблизи спинки, в обратном направлении.
Также важно отметить, что часть горячего Рис. 8 - Распределение газа, располагающегося со стороны корыта осредненной полной проскальзывает на спинку за счет эффекта температуры в среднем эжекции. Это приводит к образованию более сечении по высоте высокой адиабатической температуры поверхности лопатки в нестационарном расчете, как со стороны корыта, так и со стороны спинки. Для корневого и периферийного сечений, вторичные течения несколько отодвигают горячий газ от стенок лопатки. Этот эффект приводит к существенному различию температур непосредственно на поверхности лопатки. Так максимальная разность температур в приведенных сечениях достигает 70-75K. Однако распределение разности адиабатической температуры поверхности лопатки существенно отличается по её высоте. Эти отличия прежде всего связаны с влиянием вторичных течений вблизи концевых поверхностей. Кроме того, наличие радиальной компоненты скорости также оказывает определенное влияние на эффект сегрегации газа. На рис. показано распределение разности температуры и вектора скорости вдоль поверхностей корыта и спинки лопатки.
Рис. 9 - Распределение разности полной температуры и векторов скорости в среднем (слева) и верхнем (справа) сечениях.
Рис. 10 - Распределение разности температуры и скорости на корыте (слева) и на спинке (справа).
Поведение линий тока в межлопаточном канале ротора свидетельствует о том, что максимальное влияние эффекта сегрегации достигается в тех местах на лопатке, где вблизи стенки радиальная компонента скорости незначительна.
Второй этап применяемого подхода основан на моделировании течения в исследуемой ступени турбины с осреднением параметров потока в окружном направлении на интерфейсе между статором и ротором. Такая постановка позволила:
1. Использовать существенно более мелкую дискретизацию расчетной области для того, чтобы детально смоделировать все отверстия перфорации, как на лопатках статора, так и на роторных лопатках и концевых поверхностях.
2. Применить совместную CHT постановку для моделирования теплового состояния лопатки рабочего колеса. В данной постановке одновременно с газодинамической задачей решалась тепловая задача в металле лопатки и моделировалась теплопередача на границе газ - твердое тело.
В процессе расчета использовались различные временные шаги для газодинамической задачи и тепловой задачи. Это связано с различными временными интервалами по теплопроводности в металле и конвективному теплообмену в газе. В данном случае использовался временной шаг порядка 1e 4 для газодинамической задачи и порядка 1e+1 для тепловой задачи в металле.
Возможность использования при численном моделировании различных временных шагов была обусловлена стационарностью применяемой постановки.
Сеточная модель роторной части расчетной области состояла из 6.4 млн.
узлов (24 млн. элементов) и включала в себя межлопаточный канал, внутреннюю полость и стенки лопатки. Для описания течения в пограничном слое использовалось 8 призматических слоев. Модель статорной части расчетной области состояла из 7.2 млн. узлов (21.3 млн. элементов) и не включала в себя внутреннюю полость лопаток статора. На рис. 11 представлена дискретизация внутренней поверхности лопатки рабочего колеса вместе с описанием дискретизации каналов перфорации. Следует отметить, что на поперечное сечение каждого канала перфорации приходилось по тетраэдральных ячейки и по 16 (два раза по 8) призматических слоев. Такая дискретизация позволила корректно описать выдув охлаждающего воздуха.
Кроме того, сгущение ячеек вблизи отверстий перфорации также было подобрано таким образом, чтобы процесс смешения охлаждающего воздуха с основным потоком слабо зависел от дальнейшего измельчения сеточной модели.
Численное моделирование было выполнено со вторым порядком точности по пространству и времени. Для получения сошедшегося решения потребовалось 2-е - 3-е суток процессорного времени с распараллеливанием задачи на 12 процессоров (точнее 12 ядер на вычислительном кластере 20года выпуска).
Рис. 11 - Сеточная дискретизация лопатки рабочего колеса.
Результаты применения такого подхода представлены на рис. 12-14. В частности, на рис. 12 приведено распределение концентрации охлаждающего воздуха вблизи лопатки РК, выдуваемого через неё. Как видно из рисунка со стороны корыта можно выделить два наиболее серьезных участка, где эффективность пленочного охлаждения недостаточна (эти области выделяются, по результатам газодинамического расчета, синим цветом). Первый участок расположен вблизи торца лопатки примерно на расстоянии половины хорды лопатки от входной кромки. С учетом наличия достаточно сильных вторичных течений в этой области присутствие такого участка на лопатке является опасным с точки зрения обеспечения длительного ресурса её работы. Второй участок расположен в нижней части лопатки за последним рядом отверстий перфорации. Отсутствие синих следов за первым, вторым и пятым отверстиями перфорации в данном ряду свидетельствует о недостаточном перепаде давлений на этих отверстиях. При анализе эффективности пленочного охлаждения лопатки со стороны спинки можно отметить, что пленочное охлаждение работает слабо в области вторичных течений, что является опасным особенно для верхней части лопатки. Однако для окончательного вывода по эффективности пленочного охлаждения лопатки в этих областях требуется проанализировать рис. 13.
Рис. 12 - Концентрация охлаждающего воздуха, выдуваемого через лопатки РК Рис. 13 - Концентрация охлаждающего воздуха, выдуваемого через лопатки СА На данном рисунке представлено распределение концентрации охлаждающего воздуха вблизи поверхности исследуемой лопатки, выдуваемого через лопатку спрямляющего аппарата и концевые поверхности статора. Такое представление результатов расчета позволяет судить о суммарном влиянии охлаждающего воздуха на пленочное охлаждение лопатки рабочего колеса.
Можно отметить слабое влияние данного воздуха на охлаждение лопатки со стороны корыта. Вместе с этим сильное воздействие указанный воздух оказывает на защиту лопатки со стороны спинки в областях вторичных течений. Так концентрация данного воздуха в отмеченных выше областях достигает 50%, что оказывает существенное благоприятное влияние на тепловое состояние лопатки.
Итоговым результатом численного анализа является полученное распределение температуры в металле лопатки (рис. 14). Несколько излишний расход охлаждающего воздуха через центральный ряд отверстий перфорации со стороны корыта приводит к низким значениям температуры лопатки в этой области. Также виден высокий уровень температур на торце лопатки как со стороны корыта, так и со стороны спинки. Центральное горячее пятно со стороны спинки вызвано плохим конвективным охлаждением лопатки во внутренней полости на данном участке. Резкая диффузорность средней полости лопатки в районе корневого сечения, вследствие этого возникновение вихревого движения и дальнейшее влияние сил Кориолиса обусловили слабое конвективное охлаждение лопатки на данном участке. Для того, чтобы снизить температуру лопатки в данной области требуется модифицировать геометрию средней полости лопатки с целью добиться безвихревого характера течения и, тем самым, интенсифицировать теплообмен со стороны спинки.
Рис. 14 - Распределение температуры в металле лопатки На данном рисунке черными кривыми представлено также три сечения по высоте лопатки, в которых выполнялись экспериментальные измерения температуры металла на поверхности. Сопоставление расчетных и экспериментальных данных в отмеченных выше трех сечениях представлено на рис. 15. Расчетные данные при этом приведены для трех различных расчетных постановок: расчет ступени с распределением поля температуры на входе по экспериментальным данным (синяя кривая), по расчетным данным (зеленая кривая), а также расчет только межлопаточного канала с расчетным полем температур на входе в РК (красная кривая). Детальный анализ результатов свидетельствует, что полученное различие в расчетном распределении температур вызвано различием в определении граничных условий по температуре на входе в исследуемую ступень или на входе в рабочее колесо.
Вместе с тем, различие в полученных результатах не носит принципиального характера. Это свидетельствует о том, что на основе результатов численного моделирования, полученных с использованием любой из отмеченных выше постановок, можно выявить наиболее проблемные участки на лопатке и проработать мероприятии по улучшению её теплового состояния с высоким уровнем надежности.
спинка вх. кромка корыто Рис. 15 - Распределение расчетных и экспериментальных температур стенки лопатки в трех сечениях по высоте При анализе соответствия результатов численного моделирования и экспериментальных данных важно отметить наличие также определенной погрешности при выполнении экспериментальных измерений на натурном двигателе, в том числе вызванной имеющимися отклонениями в пределах допусков в изготовлений узлов и деталей двигателя. Об уровне такой погрешности можно судить по показаниям термокристаллов, установленных на разных лопатках. Как видно из рисунка 15 данная погрешность достаточно существенна.
Разработанная методика была применена в дальнейшем для оценки теплового состояния данной лопатки на взлетном режиме работы двигателя.
Такой расчет позволил оценить уровень тепловой нагрузки на исследуемую лопатку, определить значения максимальных температур в металле.
Верхнее сечение Среднее сечение Нижнее сечение Глава Прежде всего изложена разработанная методика оценки акустического и аэродинамического совершенства вентиляторной ступени по уровню дипольных акустических источников, генерируемых ротор-статор взаимодействием. Суть данной методики представлена на рис. 16. Она включает в себя следующие этапы: 1) нестационарный расчет ротор-статор взаимодействия; 2) определение нестационарного поля давления на лопатке спрямляющего аппарата (СА) и нахождение разности нестационарного давления на её средней поверхности; 3) Фурье-анализ полученной разности нестационарного давления; 4) осреднение по площади амплитуд разности нестационарного давления на частоте следования и кратных ей гармониках.
Полученные характеристики рассматриваются как интегральные характеристики тонального шума, вызванного ротор-статор взаимодействием.
Выявлено, что снижение осредненной амплитуды можно получить наклоном СА в сторону вращения вентилятора. В данном случае большее количество следов от лопаток вентилятора попадает на лопатку СА одновременно. Это приводит к тому, что на меньшей площади лопатки изменение давления происходит синфазно. Данная методика была апробирована на заранее заданных тестовых конфигурациях вентиляторной ступени двигателя ПС-90А и показала свою эффективность.
1 этап 2 этап 3 этап 4 этап Рис. 16 - Методика оценки акустического совершенства В дальнейшем на её основе было исследовано шесть различных вариантов конфигураций лопаток СА для перспективного авиационного двигателя ПД-14.
По результатам работ выбран такой криволинейный наклон СА, который позволил получить существенно более высокие акустические характеристики при обеспечении заданных газодинамических параметров ступени. Для получения данных по уровню шума вентиляторной ступени в области дальнего поля (где выполняются стандартные акустические измерения) необходимо было корректно смоделировать их распространение от нелинейной области ротор-статор взаимодействия в области ближнего поля. Таким образом требовалось специализированное ПО, позволяющее описать распространение акустических возмущений с учетом особенностей геометрии и неравномерного среднего потока. В результате был создан специализированный акустический решатель, основанный на решении линеаризованных уравнений Эйлера, в осесимметричной постановке. В этом случае решение моделируется во временной области для каждой окружной моды m отдельно, а исходная система уравнений записывается в цилиндрической системе координат в следующем виде:
qm qm qm 1 = - A + B + Cqm + Dqm (9) t x r r Здесь qm = {'m,u'm,v'm, w'm} - вектор неизвестных, записанный относительно пульсационных составляющих параметров потока, A, B, C, D - коэффициенты, зависящие от средних параметров потока.
Для возможности проведения расчетов в областях с криволинейными границами, используется переход от физических координат (x, r) к расчетным (, ).
Аппроксимация пространственных производных осуществляется на основе семиточечной центрально-разностной схемы DRP (Dispersion Relation Preserving) четвертого порядка точности:
N F F(xl + jx) (10) a j x l x j=- N Для обеспечения минимальных дисперсионных ошибок был использован метод наименьших квадратов, т.е. один из коэффициентов aj находился из условия минимума следующего функционала:
E = k*x - kx d(kx) min (11) Для подавления нефизичных высокочастотных пульсаций, в конце каждого шага по времени к рассчитанным газодинамическим полям применялся фильтр вида:
N Q (x0 ) = Q(x0 ) - DQ (x0 ), DQ (x0 ) = Q(x0 + jx) (12) d d j j=- N Для интегрирования по времени использовалась явная двухслойная шестишаговая оптимизированная схема Рунге-Кутты четвертого порядка точности.
Около входной и выходной границ расчетной области были созданы буферные зоны, в которых вид решаемых уравнений изменялся следующим образом:
q = -F(q)- R (q - q0) (13) damp t Коэффициент демпфирования Rdamp является экспоненциальной функцией, плавно изменяющейся от нуля вблизи внутренних границ расчетной области до заданного значения на внешних границах. Введение дополнительного члена в системе (13) позволяет добиться плавного сведения выходящих из расчетной области волн к заданному вектору q0. Таким образом УвносФ возмущений в расчетную область осуществляется посредством задания в качестве вектора qнекоторой исследуемой акустической моды (или суммы нескольких мод с одинаковым окружным числом m и возможно разными частотами и радиальными числами). УВыносФ возмущений осуществляется посредством приравнивания вектора q0 к нулю.
Разработанная программа была верифицирована на модельных тестовых задачах, имеющих аналитическое решение, а также на основе экспериментальных данных по уровню тонального шума, генерируемого вентилятором двигателя JT15D фирмы Pratt&Whitney. Эксперименты проводились в исследовательском центре NASA в начале 80-х гг. В частности в открытой печати были найдены экспериментальные диаграммы направленности тонального шума в дальнем поле для пяти режимов работы двигателя, а также расчетные данные по данному двигателю, полученные в 2004г. специалистами компании Boeing. При проведении расчетов специалистами Boeing использовался подход, аналогичный представленному в настоящей работе. Были также известны и некоторые другие параметры, необходимые для повторения расчетов: геометрия входного устройства, частоты вращения вентилятора, числа Маха на входе в рабочее колесо вентилятора, способ генерации акустических мод.
Геометрия входного устройства приведена на рис 17.
Рис. 17 - Геометрия воздухозаборника Результаты сопоставления расчетных и экспериментальных данных представлены на рис. 18 (собственные расчеты - красная кривая, расчеты Boeing - синяя). Как видно из графиков, получено хорошее соответствие результатов расчетов с экспериментальными данными. Следует отметить, что для режима 5 получено лучшее соответствие с экспериментом по сравнению с результатами Boeing.
На основе разработанной программы была выполнена оценка вклада ротор-статор взаимодействия вентиляторной ступени в суммарный уровень шума двигателя ПС-90А при стендовых испытаниях. Было показано, что данный источник тонального шума, играющий заметную роль в полетных условиях (особенно на режиме захода на посадку), вносит относительно небольшой вклад в суммарный уровень шума при стендовых испытаниях.
Взаимодействие неоднородного набегающего потока с лопатками вентилятора и шум реактивной струи сильно затеняют шум ротор-статор взаимодействия.
Режим №Режим №Режим №Режим №Режим №Рис. 18 - Диаграммы направленности для 5 режимов работы двигателя JT15D В связи с этим в дальнейшем планируется выполнение экспериментальных работ с противотурбулентным входным устройством и специальными разделительными стенками между передней и задней полусферами двигателя.
Эти конструкции позволят в стендовых условиях приблизиться к условиям полета.
Также в данной главе рассмотрены особенности разработанной методологии для расчета шума реактивных струй. В данном случае в качестве газодинамического ПО использовался пакет ANSYS-Fluent. Результаты численного моделирования по уровню акустических пульсаций в дальнем поле были сопоставлены с экспериментальными данными для двух одноконтурных моделей сопел (дозвуковой и сверхзвуковой случай), по которым имелись экспериментальные данные.
Как известно, основная цель моделирования шума реактивных струй - определение таких форм сопел, которые бы позволили максимально снизить уровень излучаемого шума (до 3-5 дБ и более). Следует однако отметить, что перед тем как приступать к решению данной задачи требуется провести исследовательскую работу по разработке методологии выполнения таких расчетов, поскольку стандартное применение ПО ANSYS-Fluent в данном случае приводит к невозможности корректно описать весь спектр турбулентных пульсаций в требуемом диапазоне частот.
В ходе работы были исследованы: различные постановки выполнения расчетов - включение и исключение внутренности сопла из расчетной области;
применение различных типов дискретизации расчетной области - структурированная, неструктурированная, смешанная; применение дополнительной адаптации; применение различных математических постановок - LES, DNS (Direct Numerical Simulation); тип и точность применяемой расчетной схемы; метод аппроксимации потоков на границах контрольных объемов и др.
Рис. 19 - Распределение давления в расчетной области и спектр пульсаций в дальнем поле на угле направленности 120 градусов В частности, в качестве дозвуковой тестовой модели была выбрана модель сопла Boeing с диаметром на срезе в 2 дюйма и числом Маха на выходе 0.9. В результате выполненной работы, удалось разработать такую методологию выполнения расчетов, которая позволяет адекватно разрешить спектр акустических пульсаций до чисел Струхаля 3,5 - 4,0. Результаты численного моделирования и их сопоставление с экспериментом представлены на рис. 19.
Основное содержание результатов диссертации опубликовано в работах:
1. Августинович В.Г. Математическое моделирование нестационарного силового и теплового нагружения лопаток турбины / В.Г. Августинович, Ю.Н.
Шмотин, А.М. Сипатов, И.А. Повышев // Труды международной научнотехнической конференции Двигатели 21 века. - М: ЦИАМ, 2000. - С.55Ц56.
2. Koutsenko I.G. Application of CFD-Based Analysis Technique for design and optimization of gas turbine combustors / I.G. Koutsenko, S.F. Onegin, A.M. Sipatov // Proceeding of ASME TURBO EXPO 2004 / June 14-17, 2004, Vienna, Austria. - ASME, 2004, paper GT2004-53398, 8P.
3. Сипатов А.М. К оценке интенсивности тонального шума вентиляторной ступени / А.М. Сипатов А.М., М.В. Усанин // журнал Авиационная техника, Известия высших учебных заведений. - Казань: КГТУ-КАИ им. А.Н. Туполева, 2004. - Вып. №4. - С.23-27.
4. Усанин М. В. Оценка интенсивности источников тонального шума взаимодействия вентилятор - спрямляющий аппарат / М. В. Усанин, А.М.
Сипатов, П.В. Трусов // Труды всероссийской научно-технической конференции Аэрокосмическая техника и высокие технологии 2004. - Пермь:
ПГТУ, 2004. - С.88Ц94.
5. Sipatov A.M. Computational Analysis of Tonal Noise Generated High-Bypass Ratio Fan Stage / A.M. Sipatov, M.V. Usanin, V.G. Avgustinovich, N.O.
Chuhlantseva // Proceeding of ASME TURBO EXPO / June 6-9, 2005, Reno-Tahoe, Nevada, USA. - ASME, 2005, paper GT2005-68857, 10P.
6. Сипатов А.М. Применение упрощенной методики для оценки различных вариантов вентиляторной ступени с точки зрения акустического совершенства / А.М. Сипатов, М.В. Усанин, В.Г. Августинович, Н.О. Чухланцева // Труды международной научно-технической конференции Авиадвигатели 21 века. - М: ЦИАМ, 2005. - Т.1. - С.89.
7. Sipatov A.M. Hybrid Approach for Computation of Tonal Noise Applied to High Bypass Ratio Fan Stage / A.M. Sipatov, V.G. Avgustinovich, M.V. Usanin, N.O.
Chuhlantseva // Proceeding of ISABE conference / September 4-9, 2005, Munich, Germany. - ISABE, 2005, paper 2005-1187, 11P.
8. Синер А.А. Моделирование тонального шума авиационного двигателя / А.А. Синер, А.М. Сипатов // Cборник научных трудов Вычислительная механика. - Пермь: ПГТУ, 2006. - Вып. №4. - С.43Ц50.
9. Bormalev S.V. Designing Quieter Fans for Turbo-Jet Engines / S.V. Bormalev, V.G. Avgustinovich, A.M. Sipatov, M.V. Usanin // J.: ANSYS Solutions, Volume 7, Issue 2, 2006. pp. 24-26.
10. Синер А.А. Моделирование тонального шума авиационного двигателя / А.А. Синер, А.М. Сипатов, В.А. Чурсин // Сборник трудов XVIII сессии Российского акустического общества / 11-15 сентября, 2006, Таганрог. - М.:
ГЕОС, 2006. - Т. 3. - C.205Ц208.
11. Сипатов А.М. Моделирование акустических характеристик вентиляторной ступени / А.М. Сипатов, Н.О. Чухланцева // журнал Математическое моделирование. - М.: РАН, 2007. - Т.19. - Вып. №8. - С.97 - 104.
12. Сипатов А.М. Применение трехмерного газодинамического анализа с целью проектирования малошумной вентиляторной ступени авиационного двигателя / А.М. Сипатов // Сборник тезисов всероссийской научнопрактической конференции Вычислительный эксперимент в аэроакустике / 27-30 сентября 2006, Светлогорск. - М.: МАКС Пресс, 2006. - С.62Ц63.
13. Сипатов А.М. Результаты численного моделирования акустических характеристик вентиляторной ступени авиационного двигателя / А.М. Сипатов, Н.О. Чухланцева, М.В. Усанин // журнал Вестник ИжГТУ. - Ижевск: ИжГТУ, 2007. - Вып. №2. - С.28Ц31.
14. Иноземцев А.А. Моделирование акустических характеристик перспективного авиационного двигателя / А.А. Иноземцев, А.М. Сипатов, М.В.
Усанин // журнал Известия самарского научного центра РАН. - Самара: СН - РАН, 2008. - Спец. выпуск посвящённый памяти А.П. Меркулова, С.85Ц92.
15. Сипатов А.М. Анализ акустических характеристик вентиляторной ступени авиационного двигателя / А.М, Сипатов, М.В. Усанин // Информационно-технический журнал Пермские авиационные двигатели. - Пермь: ОАО Авиадвигатель, 2008. - Вып. №17. - С.24Ц27.
16. Сипатов А.М. Трехмерный акустический анализ в процессе проектирования авиационного двигателя / А.М. Сипатов, М.В. Усанин, Н.О.
Чухланцева // журнал Ученые записки ЦАГИ. - М.: ЦАГИ, 2010. - Вып. №2.
- С.3Ц10.
17. Русаков С.В. Разработка методики модального анализа шума турбомашин с помощью подвижной решетки микрофонов / С.В. Русаков, А.А. Синер, А.М.
Сипатов, В.А. Чурсин // Труды всероссийской научно-технической конференции Аэрокосмическая техника и высокие технологии 2008. - Пермь:
ПГТУ, 2008. - С.337Ц339.
18. Синер А.А. Об одном подходе к решению обратных задач акустики турбомашин / А.А. Синер, С.В. Русаков, А.М. Сипатов, М.В. Усанин // журнал Ученые записки ЦАГИ. - М.: ЦАГИ, 2010. - Вып. №1. - С.53Ц58.
19. Чухланцева Н.О. Исследование перспективных численных схем газовой динамики на основе решения одномерных тестовых задач / Н.О. Чухланцева, М.В. Усанин, А.М. Сипатов // журнал Авиационная техника, Известия высших учебных заведений. - Казань: КГТУ-КАИ им. А.Н. Туполева, 2010. - Вып. №3.
- С.32-37.
20. Сипатов А.М. Численный расчет звуковых полей во входных устройствах с осевой симметрией / А.М. Сипатов, М.В. Усанин, А.А. Синер, Н.О.
Чухланцева // Сборник тезисов международной конференции Новые рубежи авиационной науки. - М.: ЦАГИ, 2007. - С.133,252Ц253.
21. Сипатов А.М. Трехмерное численное моделирование акустических источников от реактивной сверхзвуковой струи / А.М. Сипатов, М.В. Усанин, Н.О. Чухланцева // Сборник трудов XIX сессии Российского акустического общества / 24-27 сентября, 2007, Нижний Новгород. - М.: ГЕОС, 2007. - Т. 3. - C.293Ц296.
22. Сипатов А.М. Оценка шума реактивной струи на основе трехмерного газодинамического моделирования / А.М. Сипатов, Н.О. Чухланцева, М.В.
Усанин // Труды всероссийской научно-технической конференции Аэрокосмическая техника и высокие технологии 2008. - Пермь: ПГТУ, 2008.
- С.408Ц410.
23. Чухланцева Н.О. Трехмерный анализ шума реактивной струи / Н.О.
Чухланцева, М.В. Усанин, А.М. Сипатов // Материалы VII Международной конференции по неравновесным процессам в соплах и струях (NPNJТ2008) / 2431 мая 2008г., Алушта. - М.: Изд-во МАИ, 2008., стр. 430-433.
24. Sipatov A.M. Applying Fluent Software for Jet Noise Generation Modeling / A.M. Sipatov, M.V. Usanin, N.O. Chuhlantseva // Proceeding of International Acoustics Conference AIAA/CEAS / 7-9 June, 2010, Stockholm, Sweden. - AIAA/CEAS, 2010, paper 2010-3843, 8P.
25. Августинович В.Г. Численное моделирование нестационарных явлений в газотурбинных двигателях / В.Г. Августинович, Ю.Н. Шмотин, А.М. Сипатов, Д.Б. Румянцев, Д.Н. Ташлыков, Р.Ю. Старков, И.А. Повышев, А.И. Полулях // М.: Машиностроение, 2005. - 535с.
26. Загитов Р.А. Исследование статор-ротор взаимодействия в турбине высокого давления / Р.А. Загитов, А.М. Сипатов // Тезисы докладов всероссийской конференции Аэрокосмическая техника и высокие технологии 2005. - Пермь: ПГТУ, 2005 - С.98.
27. Августинович В.Г. Виртуальная среда для проектирования турбомашин / Августинович В.Г., Сипатов А.М. // Труды всероссийской научно-технической конференции Аэрокосмическая техника и высокие технологии 2007, посвященной 90-летию П.А. Соловьева. - Пермь: ПГТУ, 2007, С45.
28. Sipatov A.M. Tools for Estimating Resonant Stresses in Turbine Blades / A.M.
Sipatov, V.G. Avgustinovich, N.V. Gladisheva, I.A. Povishev // Proceeding of ASME TURBO Expo / May 14-17, 2007, Montreal, Canada. - ASME, 2007, paper GT2007-27196, 10P.
29. Иноземцев А.А. Трехмерный анализ вибронапряжений в лопатках турбомашин / А.А. Иноземцев, В.Г. Августинович, А.М. Сипатов // журнал Известия самарского научного центра РАН. - Самара: СН - РАН, 2008. - Спец.
выпуск посвящённый памяти А.П. Меркулова, С.93Ц103.
30. Августинович В.Г. Методология оценки вибронапряжений в лопатках турбин / В.Г. Августинович, А.М. Сипатов, Н.В. Гладышева, И.А. Повышев // журнал Газотурбинные технологии. - Рыбинск: Издательский дом Газотурбинные технологии, 2008. - Вып. №4, С.34Ц37.
31. Сипатов А.М. Трехмерный расчет вибронапряжений в лопатках турбин / Сипатов А.М., Гладышева Н.В., Августинович В.Г., Повышев И.А. // журнал Математическое моделирование. - М.: РАН, 2009. - Т.21. - Вып. №7. - С.106Ц120.
32. Sipatov A.M. Advanced Numerical Heat Analysis of High Pressure Turbine Blade / A.M. Sipatov, V.G. Latyshev, L.Yu. Gomzikov, N.V. Gladysheva // Proceeding of ISABE conference / September 7-11, 2009, Montreal, Canada. - ISABE, 2009, paper 2009-1217, 12P.
33. Sipatov A. Three Dimensional Heat Transfer Analysis of High Pressure Turbine Blade / A. Sipatov, L. Gomzikov, V. Latyshev, N. Gladysheva // Proceeding of ASME TURBO Expo / June 8-12, 2009, Orlando, Florida, USA. - ASME, 2009, paper GT2009-59163, 9P.
34. Гомзиков Л.Ю Трехмерное численное исследование теплового состояния рабочей лопатки турбины высокого давления / Л.Ю. Гомзиков, А.М. Сипатов, Н.В. Гладышева // Сборник трудов научно-технического конгресса по двигателестроению (НТКД2010). - М.: ЦИАМ, 2010. - С.46Ц47.
35. Сипатов А.М. Оценка теплового состояния рабочей лопатки турбины высокого давления на основе трехмерного моделирования, / А.М. Сипатов, Л.Ю. Гомзиков // журнал Известия самарского научного центра РАН. - Самара:
СН - РАН, 2010. - Вып. 4. - С.215-226.
36. Сипатов А.М. Решение многодисциплинарных задач в процессе проектирования авиационных двигателей (газовая динамика, аэроакустика, прочность) / А.М. Сипатов // Екатеринбург: УрО РАН, 2010. - 320с.
Авторефераты по всем темам >>
Авторефераты по техническим специальностям