«Балтийский государственный технический университет «ВОЕНМЕХ» им. Д.Ф. Устинова»

(БГТУ «ВОЕНМЕХ» им. Д.Ф. Устинова»)

 

  Факультет   А   «Ракетно-космическая техника»
    шифр   наименование
Кафедра   А9   «Плазмогазодинамика и теплотехника»
    шифр   наименование
Дисциплина  

«Моделирование с применением CAD/CAM/CAE-систем»

 

 

ОТЧЕТ

на тему

Моделирование движения тела в потоке
 
 

 

Выполнил студент группы

 

А9М22

 

Бутин М. А.

Фамилия И.О.

ПРЕПОДАВАТЕЛЬ

 

Яковчук М. С.

 

 

Фамилия И.О. Подпись

Оценка

 

 

«_____»

 

20____ г.

             

 

 

САНКТ-ПЕТЕРБУРГ

2018г.

Содержание

Введение 3

Цель работы 3

Геометрическая модель 3

Построение сетки 4

Настройка решателя 5

Результаты 9

Заключение 11

Введение

В данной работе проводилось компьютерное моделирование движения тела в потоке. В ходе подготовки к расчётам, были применены настройки «Скользящей сетки» Для проведения численных расчетов использовался программный комплекс ANSYS Workbench 18.0.

 

Цель работы

 

 

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

 

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

 

Построение модели проводилось в программном пакете Ansys Design Modeler.

В разделе «Создание эскиза» (Sketching) был построен эскиз и с помощью операции «Вращение» (Revolve)была получена твердотельная модель. Также была создана пластина стабилизатора. С помощью функции «Шаблон» (Pattern) были созданы еще 3 стабилизатора, аналогичные созданному ранее. Размеры исследуемого тела: длина – 450 мм ,радиус скруглённой части – 50 мм, Результат построения исследуемого тела приведен на рисунке 1.

Рисунок 1 – Исследуемое тело.

 

С помощью операции «Ограждение» (Enclosure) были созданы сферическая и кубическая расчетные области (рисунок 2), так как для придания телу трёх степеней свободы требуется фигура типа «сфера». Диаметр сферы – 500 мм, значение выбрано так, чтобы обеспечить плавный рост ячеек сетки при переходе к неподвижной сеточной области. Размер неподвижной области – 500х500х500 мм.

Рисунок 2 – Разбиение расчётной области.

 

 

Построение сетки

 

Построение сеточной модели проводилось в сеточном построителе Ansys Meshing. На рисунке 3 показана тетраэдральная сетка состоящая из 110180 ячеек.

 

 

Рисунок 3 – Сеточная модель.

 

Качество сетки приведено на рисунке 4.

 

 

Рисунок 4 – Качество сетки.

Настройка решателя

 

 

Рисунок 5 – Положение входной и выходной границы относительно тела.

 

В данном расчете используется k-epsilon realizable модель турбулентности.

 

 

Рисунок 6 – Модель турбулентности.

 

Тип входной границы задан как «velocity-inlet», скорость потока равна 20 м/с.

 

 

Рисунок 8 – Настройка входной границы.

 

Граница выхода имеет тип «Pressure Outlet».

 

 

Рисунок 9 – Настройка выходной границы.

 

На остальных гранях параллелепипеда применено граничное условие «Симметричности» (symmetry) (стенка без прилипания).

 

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

Для решения данной задачи была использована «

В расчете контролируется сила, действующая на ось x (сила лобового сопротивления или Drag force) и момент по оси z (момент тангажа или Mz).

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

 

Рисунок 10 – Настройка отображения графиков лобового сопротивления и момента тангажа.

 

Рисунок 11 – Настройка отображения силы лобового сопротивления для исследуемого тела.

 

 

Результаты

 

Рисунок 12 – График невязок.

 

 

Рисунок 13 – Отображение значений сил на границе «wall».

 

 

Рисунок 14 – Отображение значений моментов на границе «wall».

 

Рисунок 15 – Поле распределения давления по поверхности тела.

 

 

Рисунок 16 – Линии тока.

Заключение

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