Исследование аэродинамических характеристик самолета "Цикада" с помощью программы Tornado

Коэффициенты отвала поляры и геометрия самолета. Создание графиков зависимости коэффициента подъемной силы от угла атаки с использованием программы Tornado. Получение коэффициентов отвала поляры методом наименьших квадратов в математическом пакете MATLAB.

Рубрика Производство и технологии
Вид курсовая работа
Язык русский
Дата добавления 27.06.2011
Размер файла 3,3 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

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

Размещено на http://www.allbest.ru/

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РФ

САМАРСКИЙ ГОСУДАРСТВЕННЫЙ АЭРОКОСМИЧЕСКИЙ УНИВЕРСИТЕТ имени академика С. П. Королева

Кафедра аэрогидродинамики

ПОЯСНИТЕЛЬНАЯ ЗАПИСКА

к курсовой работе по механике сплошной среды на тему:

Исследование аэродинамических характеристик самолета «Цикада» с помощью программы Tornado

Выполнил:

Группа: 1405

Консультант: Ляскин Антон Сергеевич

САМАРА 2009

РЕФЕРАТ

Пояснительная записка 23 с, 22 рисунка, 3 таблицы, 5 источников

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

В данной курсовой работе необходимо исследовать аэродинамические характеристики самолета «Цикада» с помощью программы Tornado.

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

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

отвал поляр самолет tornado matlab

СОДЕРЖАНИЕ

ВВЕДЕНИЕ

1. МЕТОД ДИСКРЕТНЫХ ВИХРЕЙ

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

3. ГЕОМЕТРИЯ САМОЛЕТА

4. РЕЗУЛЬТАТЫ

ЗАКЛЮЧЕНИЕ

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

ВВЕДЕНИЕ

В данной курсовой работе необходимо исследовать аэродинамические характеристики самолета «Цикада» с помощью программы Tornado.

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

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

1. МЕТОД ДИСКРЕТНЫХ ВИХРЕЙ

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

Рисунок 1.1- К замене непрерывного распределения циркуляции дискретным (1- суммарные вихри, 2- свободные, 3- контрольные точки)

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

Рисунок 1.2- К переходу от непрерывных изменений во времени к ступенчатым

Полагается, что граничные условия и нагрузки скачкообразно изменяются в некоторые расчетные моменты времени ф=0,ф_1,…,ф_r (r=0,1,…), а в промежутках между этими моментами остаются неизменными и равными значениям этих величин в начале каждого промежутка.

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

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

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

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

Метод дискретных вихрей в данной работе:

Рисунок1.3- Разбивка кртыла на 0

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

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

Программа Tornado основана на методе дискретных вихрей. Особенности реализации метода дискретных вихрей в программе Tornado:

Вместо обычных П-образных вихрей, в программе используются семисегментные П-образные вихри, у которых «хвосты» направлены по крылу, по закрылку и потом по потоку (Рисунок 2.3);

Рисунок 2.3- Семисегментный П-образный вихрь

Профиль учитывается с помощью принципа косых нормалей (Рисунок 2.4);

Геометрия, типа угла закрутки конца крыла, учитывается явно (Рисунок 2.5);

Угол атаки учитывается явно;

Отклонение закрылка учитывается явно;

Хвосты вихрей направлены по скорости набегающего потока.

Рисунок 2.4- Замена профиля крыла косыми нормалями

Рисунок 2.5- Закрутка конца крыла

3. ГЕОМЕТРИЯ САМОЛЕТА «ЦИКАДА». ЗАДАНИЕ ГЕОМЕТРИИ В ПРОГРАММЕ TORNADO

Геометрические характеристики самолета «Цикада» (Таблица №1) получены из документации «Спецификация образца легкого самолета «Цикада»» и сняты с чертежей, представленных на рисунке 3.1.

Рисунок 3.1- Три вида самолета «Цикада»

Таблица №1

Параметр

Значение

Размерность

Первое крыло

Число разделов

3

Первый раздел

Координата центра тяжести х

0

Координата центра тяжести у

0

Координата центра тяжести z

0

Координата опорной точки х

0.35

Координата опорной точки у

0

Координата опорной точки z

0

Симметрия относительно плоскости Oxz

1

Средняя аэродинамическая хорда

1.4

м

Профиль в начале раздела

PIII

Число панелей вдоль хорды

10

Угол поперечного V

0.5

град.

Число панелей по полуразмаху

20

Полуразмах

5.46

м

Сужение

1

Профиль на конце раздела

PIII

Стреловидность по 1/4 хорд

0

град.

Угол закрутки конца крыла

0

град.

Тип сетки

1

Наличие закрылков

1

Хорда закрылка в долях к корневой хорде

0.29

Число панелей вдоль хорды закрылка

5

Симметрично ли отклоняются закрылки

1

Второй раздел

Угол поперечного V

0.5

град.

Число панелей по полуразмаху

2

Полуразмах

0.19

м

Сужение

0.84 (1.99)

Профиль на конце раздела

PIII

Стреловидность по линии 1/4 хорд

45

град.

Угол закрутки конца крыла

0

град.

Тип сетки

1

Наличие закрылков

0

Третий раздел

Угол поперечного V

90.5

град.

Число панелей по полуразмаху

2

Полуразмах

0.1

м

Сужение

1

Профиль на конце раздела

0

Стреловидность по 1/4 хорд

0

град.

Угол закрутки конца крыла

0

град.

Тип сетки

1

Наличие закрылков

0

Второе крыло (горизонтальное оперение)

Число разделов

2

Первый раздел

Симметрия относительно плоскости Oxz

1

Координата апекса х

3.85

Координата апекса у

0

Координата апекса z

-0.13

Корневая хорда

0.43

м

Профиль в начале раздела

0

Число панелей вдоль хорды

5

Поворот корневой хорды

-5

град.

Угол поперечного V

0

град.

Число панелей по полуразмаху

3

Полуразмах

0.24

м

Сужение

1.93

Профиль на конце раздела

0

Стреловидность по 1/4 хорд

25.1

град.

Угол закрутки конца крыла

-5

град.

Тип сетки

1

Наличие закрылков

0

Второй раздел

Угол поперечного V

0

град.

Число панелей по полуразмаху

10

Полуразмах

1.09

м

Сужение

1

Профиль на конце раздела

0

Стреловидность по 1/4 хорд

0

град.

Угол закрутки конца крыла

-5

град.

Тип сетки

1

Наличие закрылков

0

Построение геометрии в программе Tornado:

Профиль учитывается с помощью принципа косых нормалей (Рисунок 2.4);

Угол закрутки конца крыла учитывается явно (Рисунок 2.5);

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

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

Применяя линейную разбивку, получим геометрию, представленную на рисунках 3.2-3.7.

Рисунок 3.2- Расположение косых нормалей при линейной разбивке ()

Рисунок 3.3- Геометрия самолета «Цикада» при линейной разбивке ()

Рисунок 3.4- Расположение косых нормалей при линейной разбивке ()

Рисунок 3.5- Геометрия самолета «Цикада» при линейной разбивке ( )

Рисунок 3.6- Расположение косых нормалей при линейной разбивке ()

Рисунок 3.7- Геометрия самолета «Цикада» при линейной разбивке ()

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

Рисунок 3.8- Расположение косых нормалей при втором виде разбивки ()

Рисунок 3.9- Геометрия самолета «Цикада» при втором виде разбивки ()

Рисунок 3.10- Расположение косых нормалей при втором виде разбивки ( )

Рисунок 3.11- Геометрия самолета «Цикада» при втором виде разбивки ()

Рисунок 3.12-Расположение косых нормалей при втором виде разбивки ()

Рисунок 3.13- Геометрия самолета «Цикада» при втором виде разбивки ()

4. РЕЗУЛЬТАТЫ

Результаты получаем в шестом пункте меню программы Tornado в виде графиков, представленных на рисунках 4.1- 4.4.

Рисунок 4.1- Зависимость коэффициента подъемной силы от угла атаки б при линейной разбивке

Рисунок 4.2- Зависимость коэффициента подъемной силы от угла атаки б при разбивке по косинусам- полукосинусам

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

Рисунок 4.3- Поляра при линейной разбивке

Рисунок 4.4- Поляра при разбивке по косинусам- полукосинусам

Таблица №2

, град

л при линейной разбивке

л при разбивке по косинусам- полукосинусам

-515.6071

-528.0517

-43.6265

-41.9238

-14.2480

-13.1831

В таблице №2 приведены коэффициенты отвала поляры, полученные методом наименьших квадратов в математическом пакете MATLAB при разном разбиении.

Таблица №3

, град

при линейной разбивке

при разбивке пол косинусам- полукосинусам

1.043

1.0618

1.925

1.9575

2.5247

2.5691

В таблице №3 приведены коэффициенты подъемной силы, полученные при разном разбиении.

ЗАКЛЮЧЕНИЕ

В данной курсовой работе исследованы аэродинамические характеристики самолета «Цикада» с помощью программы Tornado.

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

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

Белоцерковский О.М. Численное моделирование в механике сплошных сред.--М.:Физматлит,1994.

Tomas Melin A Vortex Lattice MATLAB Implementation for Linear Aerodynamic Wing Applications

Tomas Melin, User's guide and reference manual for Tornado, Royal Institute of Technology (KTH), Department of aeronautics, 2000-12

Размещено на Allbest.ru


Подобные документы

  • Получение путем расчета аэродинамических характеристик самолета Ту-214 в диапазоне изменения высот и чисел Маха полета. Вычисление геометрических характеристик самолета. Подбор аэродинамического профиля крыла и оперения. Полетная докритическая поляра.

    курсовая работа [1,5 M], добавлен 15.02.2014

  • Конструктивно-аэродинамическая компоновка самолета-высокоплана АН-24. Определение аэродинамических характеристик самолета. Подъемная сила и сила сопротивления, их распределение по поверхности. Механизмы возникновения подъемной силы и силы сопротивления.

    контрольная работа [1,2 M], добавлен 29.05.2013

  • Статистическое проектирование облика самолета. Назначение, тактико-технические требования к самолету, условия его производства и эксплуатации, определение аэродинамических и технических характеристик. Разработка технологии изготовления детали самолета.

    дипломная работа [2,6 M], добавлен 21.11.2011

  • Техническое описание самолета. Обоснование проектных параметров. Расчет взлетной массы. Компоновка и расчет геометрических параметров основных частей самолета. Коэффициент максимальной подъемной силы. Определение летно-эксплуатационных характеристик.

    курсовая работа [891,2 K], добавлен 27.06.2011

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

    курсовая работа [1,0 M], добавлен 26.11.2010

  • Схема металлоконструкции неповоротного отвала бульдозера. Описание металлоконструкции, ее элементов, материалов, технологии изготовления и сборки. Определение свойств этой детали в APM Structure 3D. Граничные условия и нагрузки, статический расчет.

    курсовая работа [403,1 K], добавлен 31.01.2016

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

    методичка [500,7 K], добавлен 29.01.2010

  • Описание и анализ надежности шасси самолета Ту-154. Конструктивные усовершенствования тормозного цилиндра и дисков колес, расчет энергоемкости тормоза. Механизмы технического сервиса и разработка передвижной установки обслуживания шасси самолета.

    дипломная работа [1,2 M], добавлен 15.08.2010

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

    дипломная работа [1,9 M], добавлен 18.07.2011

  • Изучение метода Симою и реализация программы в среде MatLab. Определение коэффициентов передаточной функции методом площадей Симою по заданным значениям переходного процесса на определенном интервале времени. Расчет системы автоматического регулирования.

    лабораторная работа [152,5 K], добавлен 22.09.2016

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.