Дописал для программы МКЭ-Ф динамический отклик стержневой конструкции. Вообще говоря, результаты неоднозначные. Отличие от того, что выдает ANSYS, около 10%. На правду похоже и программа работает, но есть над чем подумать. В статье приведу результаты расчетов динамического отклика для той фермы, что рассматривалась в ANSYS ранее.
Работа в МКЭ-Ф v0.3
Программа переписана в объектно-ориентированном стиле. Поэтому теперь изменился ход работы. Вначале создается объект класса StructFEProblem:
P1 = StructFEProblem('CaseATransientSite.txt');
Здесь создан объект P1. Единственным аргументом в конструктор объекта передается имя текстового файла расчетного случая. Этот файл также немного поменялся. Об этом чуть позже.
Далее вызывается метод запуска нужного типа решателя:
P1.RunTransient(0.00005,0.1,5,2);
RunTransient — метод анализа на динамический отклик. В аргументах: 0.00005 — шаг по времени, 0.1 — общее время моделирования, 5 — интересующий узел стержневой системы, 2 — номер интересующей степени свободы в узле.
Файл CaseATransientSite.txt находится в архиве программы. Пояснения к файлу даны в самом файле, под данными.
Важная особенность. Программа пока может только прикладывать нагрузку, но не может её убрать. Поэтому, в нашем случае, колебания конструкции происходят под действующей в 5-м узле нагрузкой -1000 Н.
Отклик стержневой конструкции (ферма)
Результаты статического и модального анализа
Результаты моделирования отклика
P1.RunTransient(0.00005,0.5,5,2);
P1.RunTransient(0.00005,0.2,5,2);
Исходники
Программу можно скачать на соответствующей странице моего сайта.