Сайт Даниила Денисова
В рубрике собраны материалы по методу конечных элементов.
Обновлена программа МКЭ-Ф. Наверное самым большим нововведением стал единый веб-интерфейс. Больше не нужно вводить команды. Хотя при желании — консоль Octave в вашем распоряжении.
Онлайн версия (ограничена по вычислительной мощности): MKE-F & BGModeller
Есть ещё один способ создать сетку для экспериментов в области численных методов. Сетка в 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. Получилось довольно интересно, есть простор для экспериментов по изменению тестовой и весовой функции.
Рассматривается способ вывода на экран матрицы масс и матрицы жесткости. Работаем в МКЭ-Ф и Matlab, поэтому, остановив решение отладчиком в нужном месте уже можно посмотреть, как выглядит матрица масс и матрица жесткости конечного элемента.
Рассматривается статический расчет рамы в МКЭ-Ф. Также приводятся результаты расчета в студенческой версии ANSYS. Следует отметить, что в реализации МКЭ-Ф балочный элемент отличается от Beam188. Вероятно поэтому очень сильного сходства результатов для рамы добиться не удалось (разница в 5%). Необходимо также проверить программу на наличие багов и проанализировать задачу о консольной балке, хоть там и получились замечательные результаты в сравнении с формулами сопротивления материалов (рассмотрена тут).
Я периодически дорабатываю программу МКЭ-Ф. При необходимости верификации результатов, создаются модели стержневых конструкций в ANSYS. Чтобы в разных статьях не повторять описание тривиальных задач по заданию стержневых систем, решил выделить в отдельную статью сам процесс моделирования стержневой системы. Потом буду просто ссылаться, тем более, процесс не сложный и одинаков для любой стержневой системы. Читать далее «Модели стержневых конструкций в ANSYS»