Сайт Даниила Денисова
Для дальнейшей разработки программы Aero-PM понадобилась генерация профиля NACA в Matlab. То есть процедура, которая могла бы генерировать координаты четырехзначных профилей NACA с разной хордой и выдавала бы раздельно верхнюю поверхность, нижнюю и среднюю линию. В релиз пока не добавил, но с исходников на GitHub скачать уже можно. Файл NACA4Digit.m.
В вычислительной аэродинамике есть семейство методов, называемых потенциальными. Это например, метод дискретных вихрей, метод вихревых панелей или вихревых решеток. По этому семейству методов есть несколько книг. Наиболее детально методы изложены в книге J. Katz, A. Plotkin — Low-speed Aerodynamics.
Программа Aero-PM v.0.1.0 реализует самый простой метод дискретных вихрей, в соответствии с книгой. В данной заметке опишу полученные результаты для пластинки при малых углах атаки. Программа написана на Matlab и размещается в репозитории Aero-PM на GitHub.
Есть ещё один способ создать сетку для экспериментов в области численных методов. Сетка в Salome строится, как правило в графическом режиме. Рассмотрим сегодня вот это:
Программа Salome бесплатна. Скачать можно тут. Есть версия под винду, она очень долго распаковывается и проверяется антивирусом (20 минут у меня вышло). Итак, приступим.
Предположим, что некоторая программа получила скалярное поле в результате своих расчетов. Поле сохранено в текстовом виде, в файле CSV (Comma-separated values). Рассмотрим, как происходит визуализация CSV в ParaView на несложном примере.
Теперь рассмотрим чтение сетки Gmsh в Matlab. Сетка Gmsh должна быть сохранена в msh-файл. Версия формата «2.2 0 8». Она указывается в начале msh-файла. Не исключено что работать будут и другие версии. Пишем на Matlab класс, конструктором которого является функция с аргументом в виде имени текстового msh-файла. Поля и методы сделаем доступными (public). Читать далее «Чтение сетки Gmsh в Matlab»
Продолжаю тему популярного открытого сеточного генератора. Теперь рассматривается структурированная 2D сетка в Gmsh. Вот такая:
С наступившим 2018-м годом! Для численных экспериментов в области сеточных методов, полезно иметь возможность строить простые сетки. Рассмотрим, как создается неструктурированная 2D сетка в Gmsh 3.0.
Фундаментальная тема 🙂 Написал программу-блокнот, реализующую метод взвешенных невязок в Mathematica. Получилось довольно интересно, есть простор для экспериментов по изменению тестовой и весовой функции.
В справочнике авиаконструктора под редакцией Лурье и Горяинова (ЦАГИ, 1939 год, том 3, стр. 623) приводится формула, позволяющая выполнить экспериментальное определение момента инерции. Вот такая:
Формула позволяет определить момент инерции относительно оси, проходящей через центр масс тела (центр тяжести).
Рассмотрим, что такое физический маятник и как получилась такая формула.
Читать далее «Экспериментальное определение момента инерции»
Рассматривается способ вывода на экран матрицы масс и матрицы жесткости. Работаем в МКЭ-Ф и Matlab, поэтому, остановив решение отладчиком в нужном месте уже можно посмотреть, как выглядит матрица масс и матрица жесткости конечного элемента.