Доработал программу МКЭ-Ф до версии v0.5. Статья о том, как моделировать частотный отклик фермы в МКЭ-Ф. Приводятся графики результатов для разных частот.
Частотный отклик фермы в МКЭ-Ф: Работа программы
Рассматривается ферма из задачи А. В первых рассматриваемых запусках, до резонанса, гармоническая (синусоида) сила прикладывается к 5-му узлу вниз.
Входные файлы изменились! В том месте, где задается величина нагрузки, теперь добавилась ещё частота, а тип ГУ сменился с 10 на 11. Для типа 10 и 11 длина строки разная. А надо одинаковую, поскольку dlmread читает прямоугольный массив. Придется добавлять ноль справа в строку статической нагрузки (тип 10). Это «пустая ячейка», препроцессор прочитал строку, и увидев тип ГУ=10, ничего с этим нулем не делает, он только для размера строки.
В данной статье рассматриваем силу приложенную к узлу 5. Амплитуда 1000Н, частоты: 50 Гц, 100Гц, 120Гц, 140Гц и 160Гц. Для каждой частоты придется создавать отдельный файл т.к. сила одна и та же.
Кроме того, в конце статьи пример вывода программы для трех гармонических нагрузок с разными направлениями, частотами и амплитудами. Если читатель будет экспериментировать, необходимо помнить, что третья компонента силы не работает (оси Z нет — ферма плоская).
Запуск выглядит так:
P1 = StructFEProblem('CaseAFreqResponse120Hz.txt'); P1.RunTransient(0.0000005,0.2,5,2); clear all; clc;
Так для каждой задачи. Clear all (удалить все переменные) и clc (очистить экран) я делаю всегда после работы больших программ, на всякий случай, чтобы текущие переменные не повлияли на работу.