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

Тип работы:
Диссертация
Предмет:
Механика
Страниц:
108


Узнать стоимость

Детальная информация о работе

Выдержка из работы

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

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

В настоящее время наибольшее развитие и применение в инженерной практике находят методы первой группы (см., например, [41], [42], [6], [56]). Суть таких методов состоит в последовательном решении прямой задачи с последующей модификацией формы профиля для достижения наилучшего совпадения получаемых свойств с желаемыми.

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

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

Теоретическую основу обратных методов аэродинамического проектирования профилей крыльев и гидродинамических решеток, составляют обратные краевые задачи аэрогидродинамики (ОКЗА) (см., например, [49], [47], [12], [13], [14]), являющиеся частью общей теории обратных краевых задач (ОКЗ). В отличие от прямых краевых задач, в которых требуется найти функцию, удовлетворяющую в заданной области некоторому дифференциальному уравнению в частных производных, а на границе — заданному условию, в ОКЗ граница области и функция в этой области отыскиваются по дополнительному краевому условию на искомой границе.

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

Одно из первых полных исследований обратной задачи по f (s) дано в работе W. Mangler’a [60], где использован так называемый прицнип сопоставления плоскостей. W. Mangler установил зависимость между дуговыми абсциссами контура профиля и окружности в канонической области и свел ОКЗА к решению прямой задачи (задачи Шварца) в канонической области. В этой работе были также получены условия разрешимости задачи и предложен способ их удовлетворения. Г. Г. Тумашев при исследовании ОКЗА в постановке W. Mangler’a ввел в рассмотрение вспомогательную область в плоскости комплексного потенциала течения, что позволило не только дать другой способ построения аналитического решения краевой задачи, но и рассмотреть в дальнешем ряд других ОКЗА (см., например, [48], [49]).

Метод решения ОКЗА существенным образом зависит от выбора параметра, в функции которого задается распределение скорости. Результаты по решению ОКЗА в случае, когда распределение скорости задается не в виде функции дуговой координаты, получены В. М. Шурыгиным [53] (i> = v (x)), Г. Ю. Степановым [47] (v = ^(7), v = v (6)). Л. А. Симонов (см., например, [45], [46]) дал постановку и построил решение ОКЗА по годографу скорости.

Другой существенной особенностью ОКЗА является тот факт, что произвольным исходным данным соответствует, как правило, физически нереализуемое решение задачи, то есть контур получаемого профиля может оказаться незамкнутым и самопересекающимся. Кроме того, заданная скорость Voo набегающего потока может не совпадать с величиной скорости, определяемой в ходе решения задачи. Таким образом, эти задачи являются некорректными. Условия, обеспечивающие замкнутость контура (условия замкнутости) и согласование указанных величин скорости, называются условиями разрешимости. По существу условия разрешимости содержатся в работах А. Бетца [54] и подробно выведены в статьях В. Манглера [60], М. Лайтхилла [57] и Г. Г. Тумашева [51].

Основываясь на общей идее В. К. Иванова (см., например, [17]) построения квазирешения некорректных задач, А. М. Елизаровым [9] введено определение и доказана корректность квазирешения ОКЗА. Далее в работах А. М. Елизарова, Н. Б. Ильинского, А. В. Поташева (см., например, [13]) разработан метод построения квазирешения ОКЗА. Суть его состоит в минимальной коррекции исходного распределения скорости с тем, чтобы удовлетворить условиям разрешимости. Однако, применение квазирешения для замкнутости искомого контура иногда может очень сильно изменить исходное распределение скорости или давления.

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

Первый подход к нахождению оптимальных форм состоит в постановке и решении & quot-модельных"- задач. При их формулировке, как правило, используются простейшие модели. Это позволяет построить аналитические решения и сделать выводы об их оптимальных аэродинамических характеристиках. Именно такой подход использовался в работе М. А. Лаврентьева [22], который доказал, что решением вариационной задачи о дуге максимальной подъемной силы заданной длины и ограни-ченнной кривизны в потоке идеальной несжимаемой жидкости является дуга окружности.

Аналогичные результаты были получены и другими авторами. В. И. Зубов [16] предложил численный метод оптимизации формы профиля, обладающего максимальной подъемной силой, при дополнительных ограничениях, выражающих приближенно условие безотрывности обтекания. Этот метод основывается на интегральном представлении функции, конформно отображающей единичный круг на внешность профиля. Через коэффициенты Фурье этой функции выражена величина подъемной силы, максимум которой достигается за счет варьирования этих коэффициентов при указанных дополнительных ограничениях. А. М. Елизаровым [10], Д. Ф. Абзалиловым, Н. Б. Ильинским, Р. Ф. Мардановым [3] даны постановки и решения задачи об отыскании среди профилей с заданной длиной контура такого профиля, который обладает максимальной подъемной силой при некоторых ограничениях на форму профиля и распределение скорости на его поверхности.

При анализе и использовании результатов решения таких & quot-модельных"- задач, естественно, следует учитывать, что они показывают пути совершенствования аэродинамических форм. Получаемые при этом объекты бывают физически нереализуемы.

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

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

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

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

Более быстродействующие методы базируются на описании оптимизируемых величин и ограничений через решение ОКЗА. В качестве управляющей в данном подходе выбирается некоторая вспомогательная функция, позволяющая с малыми вычислительными затратами отыскивать требуемые величины (см., например, [11], [2], [15]).

Целью настоящей работы является разработка методов оптимизации форм крыловых профилей и профилей гидродинамических решеток. Сначала внимание уделено решению & quot-модельных"- задач. К ним относятся задачи построения профилей с кусочно-постоянным распределением скорости. В данном классе задач удалось аналитически определить максимум коэффициента подъемной силы. Далее были поставлены и решены оптимизационные задачи (бесконечно-тонких профилей, гидродинамической решетки, телесных профилей), решение которых осуществлялось с помощью итерационных процедур с модификацией геометрии контура для оптимизации аэродинамических характеристик на каждом шаге итерации. Основное преимущество разработанных методов оптимизации данной диссертации является то что, выбор математической модели и метода решения, позволили записать оптимизируемый функционал и накладываемые ограничения записываются через одну функцию, которая выбиралась в качестве управляющей. При этом решение задачи отыскивается сразу в физической области, что является существенным отличием от работ [2], [3], [10], [15], где для построения решения использовалась вспомогательная область комплексного потенциала.

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

Заключение

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

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

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

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

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

ПоказатьСвернуть

Содержание

Используемые аббревиатуры и обозначения.

I. Высоконесущие профили постоянной скорости

§ 1. Профили с постоянными скоростями на верхней и нижней поверхностях.

§ 2. Профили со скачком скорости на верхней поверхности

§ 3. Профили с обдувом верхней поверхности.

§ 4. Профили при наличии плоского экрана.

II. Оптимизация бесконечно-тонких профилей

§ 5. Профиль в неограниченном потоке.

§ 6. Профиль над экраном.

§ 7. Оптимизация формы бесконечно-тонких профилей прямой однорядной гидродинамической решетки.

III. Профили конечной толщины

§ 8. Профиль конечной толщины в неограниченном потоке

§ 9. Профиль конечной толщины при наличии экрана.

Список литературы

1. Абзалилов Д. Ф., Ильинский Н. Б., Степанов Г. Ю. Построение крылового профиля с отбором внешнего потока // Известия РАН МЖГ. 1996. № 6. С. 23−28.

2. Абзалилов Д. Ф., Ильинский Н. Б. Построение и оптимизация крыловых профилей с отбором внешнего потока // Ученые записки ЦАГИ. 1998, Т. 29, — № 3−4, — С. 52−61.

3. Абзалилов Д. Ф., Ильинский Н. Б., Марданов Р. Ф. Задача максимизации циркуляции скорости при обтекании гладкого контура с источниками и стоками // Журнал вычислительной математики и математической физики, 2000, том 40, № 1, С. 82−90.

4. Аульченко С. М. Метод оптимизации профилей в дозвуковом потоке идеального газа. Ин-т теор. и прикл. мех. СО АН СССР (Новосибирск) Препр. 1987. № 30−87, 45 с.

5. Боллхауз У. Ф. Некоторые новейшие достижения в численном исследовании трансзвуковых течений // В кн. Численные методы в динамике жидкости. М. :Мир, 1981, С. 152−242.

6. Викторов Г. В. Гидродинамическая теория решеток. — М.: Высшая школа, 1969, 368 с.

7. Гуревич М. И. Теория струй идеальной жидкости. М.: Наука, 1979. — 536 с.

8. Дорфман JI. А. Расчет безвихревого обтекания решеток профилей и построение решеток по заданному распределению скоростей на профилях // Прикладная математика и механика, 1952, 16, № 5, С. 599−612.

9. Елизаров A.M. О квазирешениях внешней обратной краевой задачи // Известия вузов. Математика. 1984. — № 10. — С. 42−50.

10. Елизаров А. М. Некоторые экстремальные задачи теории крыла // Известия вузов. Математика. 1988. — № 10. — С. 71−74.

11. Елизаров А. М., Фёдоров Е. В. Оптимизация формы крыловых профилей методом обратных краевых задач // Труды семинара по краевым задачам. Казанский университет, 1989, № 25, С. 20−25.

12. Елизаров А. М., Ильинский Н. Б., Поташев А. В. Обратные краевые задачи аэрогидродинамики // Итоги науки и техники. Механика жидкости и газа. М.: ВИНИТИ. — 1989. — Т. 23 — С. 3−115.

13. Елизаров А. М., Ильинский Н. Б., Поташев А. В. Обратные краевые задачи аэрогидродинамики. М.: Наука, 1994. 440 с.

14. Елизаров А. М., Ихсанова А. Н., Фокин Д. А. Численная оптимизация формы крыла экраноплана методами вариационных обратных краевых задач // Обозрение прикладной и промышленной математики. Том.8. Москва: Научное издательство & quot-ТВП"-, 2001., С. 165 167.

15. Зубов В. И. К вопросу об оптимальном профиле крыла в потоке идеальной несжимаемой жидкости // Журнал вычислительная математика и математическая физика, 1980, № 1, С. 241−245

16. Иванов В. К., Васин В. В., Таната В. П. Теория линейных некорректных задач и её приложения. М.: Наука, 1978. 206 с.

17. Ильинский А. Н., Ильинский И. Б., Поляков Д. В., Поташев А. В Степанов Г. Ю. Уточнение критерия отрыва турбулентного пограничного слоя с использованием эмпирических данных // Казань. Препринт НИИММ КГУ, № 98−2, 1998.

18. Киселев О. М. Построение прямой однорядной решетки по заданной хордовой диаграмме // Известия вузов. Авиационная техника, I960, № 4, С. 31−39.

19. Кочин Н. Е. Гидродинамическая теория решеток. M. -J1.: Госте-хиздат, 1949. — 103 с.

20. Лаврентьев М. А. Об одной экстремальной задаче в теории крыла аэроплана. // М.- Л.: Гостехтеориздат, 1934. 40 с.

21. Лаврентьев М. А., Шабат Б. В. Методы теории функции комплексного переменного. М.: Наука, 1973. 736 с.

22. Леонтьев В. Г. Высоконесущие модельные крыловые профили с двумя участками постоянной скорости на контуре / В. Г. Леонтьев, А. В. Поташев // Известия РАН. МЖГ. 2001, № 6. С. 15−20.

23. Леонтьев В. Г. К проблеме построения высоконесущих модельных аэродинамических профилей / В. Г. Леонтьев, А. В. Поташев Ц Известия Вузов. Авиационная техника. 2002. № 2. С. 24−28.

24. Леонтьев В. Г. Численное решение задачи Лаврентьева для дужки вблизи горизонтального экрана // VIII Четаевская международная конференция & quot-Аналитическая механика, устойчивость и управление

25. Ф движением& quot-: Тезисы докладов. Казань: Издательство Казанскогогосударственного технического университета, 2002 С. 267.

26. Леонтьев В. Г. Численное решение об оптимизации дужки вблизи горизонтального экрана // Гидродинамика больших скоростей. Тезисы докладов Международной летней научной школы. Казань:

27. Издательство казанского математического общества, 2002 С. 99 100.

28. Лойцянский Л. Г. Механика жидкости и газа. М.: Наука, 1987. — 840 с.

29. Маклаков Д. В. Нелинейная задача о движении профиля произвольной формы вблизи границы раздела сред // Тр. семинара по краевым задачам. Казань: казан, ун-т. — 1984. — Вып. 21. — С. 126 133.

30. Маклаков Д. В. О существовании решения задачи о движении-профиля произвольной формы в потоке двухслойной жидкости // Известия вузов. Математика. 1985 — № 6, С 30−36.

31. Нугманов 3. X., Овчинников В. А., Павлов В. Г. Аэродинамическое проектирование с учетом условия безотрывности. Известия вузов. Авиационная техника, 1985, № 3, С. 47−50.

32. Нугманов 3. X., Овчинников В. А., Павлов В. Г., Романов В. М. Численные методы расчета обтекания профиля идеальным несжимаемым потоком. Казань: Казанский авиационный институт, 1986, 64 с.

33. Поташев А. В. Построение крылового профиля с закрылком конечных размеров // Известия РАН. МЖГ, 1995, № 1. С. 173−180.

34. Рыжик И. М., Градштейн И. С. Таблицы интегралов, сумм, рядов и произведений. Изд. 3-е, перераб. M. -J1.: Гостехиздат, 1951. -464 с.

35. Симонов JI. А. Построение профилей по годографу скоростей. -Прикладная математика и механика, 1940, 4, № 4, С. 97−116.

36. Симонов JI. А. Построение профилей по годографу скоростей. -Прикладная математика и механика, 1941, 5, № 2, С. 193−222.

37. Степанов Г. Ю. Гидродинамика решеток турбомашин. М.: Физ-матгиз, 1962. 512 с.

38. Тумашев Г. Г., Нужин М. Т. Обратные краевые задачи. Уч. зап. Казан, ун-т. 1955, 115, № 6 167 с. — РЖМат, 1956, 7298.

39. Тумашев Г. Г., Нужин М. Т. Обратные краевые задачи и их приложения. Казань: Казанский университет. — 1965. — 333 с.

40. Тумашев Г. Г. Построение профилей по заданному распределению скоростей // Уч. зап. казан, ун-та. 1949. Т. 109. № 1. С. 73−87.

41. Тумашев Г. Г. Определение формы границ потока жидкости по заданному распределению скорости или давления // Уч. зап. казан, ун-та. 1952. Т. 112. № 3. С. 3−41.

42. Черноусько Ф. ji. О движении идеальной жидкости с разрывом давления вдоль границы // ПММ. 1962. Т. XXVI. Вып.2. С. 373 375.

43. Шурыгин В. М. Определение контура профиля по заданному распределению давления // Тр. ЦАГИ. 1948. — № 660. — 20 с.

44. Betz A. Anderung der Profilform zur erzielung einer vorgebenen Anderrung der Druckverteilung // Luftfahrtforschung. 1934. — 11. -№ 6. — S. 158−164. № 978, 11 p.

45. Costello G. R. Method of designing cascade blades with prescribed velocity distributions in compressible potential flows. NACA Rept., 1950, № 978, 11 p.

46. Hirose N., Tkanashi S., Kawai N. Transonic airfoil design based on Navie-Stocks equations // Techn. Report Nat. Aerospace Lab. 1986 № 901. — P. l-16 № 978, 11 p.

47. Lighthill M. J. A new method of two-dimensional aerodynamic design // Aeronaut, res. Counc. Repts. and Mem. 1945. — № 2112 — 53 p.

48. Lighthill M. J. A mathematical method of cascade design Aeronaut res. Counc. Repts. and Mem. — 1945. — № 2104 — 18 p.

49. McBride E. J. Blade profiles. Proc. First U. S. Nat. Congr. Appl. Mech., Publ. Amer. Soc. Mech. Engrs, N. Y., 1952, 699−704 — РЖМех, 1954, 4061.

50. Mangier W. Die Bereshung eines Tragflugelprofiles mit vorgeshriebener Drucverteilung // Jahrb. Dtsch. Luftfarhrtforchung. 1938. Bd.l. S. 46−53

51. Nasyrov S. Robin capacity and lift of infinitely thin airfoils. // Complex Variables, 2002, Vol. 47, № 2, pp. 93−107.

52. Nelder /. A., Mead R. A simplex method for function minimization // Computer Journal. 1965. № 7, pp. 308−313.

53. H. B. Squire, A. D. Young. The calculation of the profile drag of airfoils- ARC Rep. and Mem. 1838 (1938).

Заполнить форму текущей работой