Для уточнения оценок статистических характеристик можно проводить сглаживание по ансамблю последовательно вычисленных временных реализаций: после вычисления очередной реализации ансамбля (или пары реализаций для взаимных характеристик) пользователь получает на экране результат текущего осреднения.
Принципы построения. Пакет DISSAN имеет 7 базовых состояний. Это динамическое и статическое состояния, состояния установки, манипуляций с графиком, вычисления статистических характеристик, отображения сечений поверхности и установки параметров внешней нагрузки.
В состоянии установки задаются тип исследуемой системы, ряд ее основных параметров и характеристики, связанные с представлением результатов моделирования. При выборе типа исследуемой системы определяется объект анализа: это либо непосредственно функции, относительно которых составлена система уравнений, либо распределения перемещений, напряжений или ускорений по исследуемому участку поверхности. Если предполагается анализ таких распределений, то соответствующие собственные формы должны быть вычислены заранее и записаны в специальный файл, что можно сделать, например, методом коррекции, наложением упруго-инерционных связей (КНУИС) [4, 1] или методом конечных элементов.
В динамическом состоянии производится интегрирование системы уравнений (решается задача Коши). В режиме "шаг за шагом" каждый шаг интегрирования и его отображение происходят после нажатия соответствующей клавиши. В автоматическом режиме все наоборот: процесс интегрирования продолжается до нажатия любой из клавиш. Текущий результат интегрирования динамически отображается на экране в виде фазовых портретов, отображений Пуанкаре или изменяющихся изображений распределений исследуемых величин. В данном состоянии возможна установка некоторых параметров уравнений и процесса интегрирования, выбор анализируемых точек поверхности, вычисление и отображение распределений средне-квадратических значений.
В статическом состоянии возможна установка параметров дифференциальных уравнений (в том числе для внешней нагрузки) и процесса интегрирования, а также выбор состояния для анализа поведения конструкции. Переход в статическое состояние вызывает масштабирование текущего графика на экране по максимальным и минимальным значениям. Все другие способы отображения доступны в состоянии манипуляций с графиком.
Для численного интегрирования используется модифицированный предиктор-корректор — метод Хемминга. Поскольку указанный метод на рассматриваемом классе задач по совокупности критериев наиболее предпочтителен, то, имея возможность влиять на точность и скорость вычислений, изменяя шаг интегрирования в интерактивном режиме, признано нецелесообразным предлагать пользователю другие методы интегрирования.
Исследуемую модель необходимо представлять в нормальном виде, то есть как систему обыкновенных дифференциальных сравнений первого порядка, разрешенных относительно производных (ее можно получить, если подходящие производные ввести в качестве новых переменных). Далее эта система кодируется по простым правилам на языке ФОРТРАН-77 в виде модуля с заранее оговоренным списком входных и выходных параметров, который отдельно транслируется и включается в исполняемый модуль пакета DYSSAN.
Для вычисления внешней нагрузки используется процесс авторегрессии, параметры которого могут быть скорректированы с помощью специального меню. При этом на экран выводится график спектральной плотности нагрузки с сопутствующей информацией, по которым судят о том, насколько результаты коррекции далеки от желаемых.
Пользователь может запоминать текущее состояние процесса моделирования и возобновлять работу во время следующих сеансов.
Все полученные на экране графики и изображения могут быть выведены на стандартный принтер.
Интерфейс. Пакет DYSSAN имеет простой и удобный интерфейс, построенный на использовании хорошо отработанной системы меню. Каждому базовому состоянию соответствует свое меню, постоянно находящееся в нижней части экрана. Из базовых меню доступны так называемые вызываемые меню (с параметрами или без них) (рис. 3 и 4).
Для выбора необходимых при анализе точек графиков и поверхности конструкции, а также для выбора в базовых меню используются стрелка-указатель, перемещаемая по экрану устройством типа "мышь", либо функциональные клавиши.
На экране постоянно находится строка с подсказкой. Кроме того, из каждого базового меню можно выйти в режим help для просмотра справочника для пользователя.
Опыт использования. Пакет DYSSAN, как и предшествовавшая ему разработка - пакет ЭЛЕРОН для моделирования на базе ЭВМ серии ЕС [2], - создавался и использовался в рамках технологии сопровождающего моделирования при исследованиях акустического нагружения конструкции самолета [3].
Общая схема этой технологии включает в себя два этапа. Первый этап заключается в выборе расчетной схемы конструкции после анализа соответствующей документации и в формировании модели движения, которая представляется динамической системой, полученной разложением движения по собственным формам колебаний. Собственные формы и частоты вычисляются методом КНУИС.
На втором этапе проводится вычислительный эксперимент - моделирование динамического поведения исследуемых конструкций с помощью пакета DYSSAN.
С целью определения степени адекватности модели и реальной конструкции результаты обоих этапов моделирования сравниваются с экспериментальными данными. При недостаточной адекватности корректируется расчетная схема, после чего делается следующая попытка.
По окончании анализа результатов моделирования вырабатываются предложения по планированию дальнейших экспериментов и оценивается степень совершенства конструкции. При необходимости конструкция модифицируется, и для ее нового варианта процесс моделирования повторяется заново.
Мы рассмотрели одно из возможных применений пакета DYSSAN. Он удобен при решении и других задач: как сложных, связанных, например, с исследованиями флаттера, так и более простых - в частности, для обучения студентов некоторым разделам нелинейной динамики.
Список литературы
1. Куравский Л.С. Метод коррекции наложением упругоинерционных связей в задаче о свободных колебаниях тонкостенных авиационных конструкций. В сб. // Вопросы авиационной техники - Жуковский: ЛИИ, 1993- вып. 239.- С.56-72
2. Куравский Л.С. Принципы работы диалоговой системы для моделирования виброакустического нагружения тонкостенных авиационных конструкций // Там же. – С.48-55
3. Куравский Л.С, Арнаутов Е.В. Расчетные оценки параметров колебаний тонкостенных конструкций при акустическом возбуждении // X Юбилейная науч.-техн.конф. по аэроакустике (сент.1992, Суздаль). Москва: ЦАГИ, 1992.- Ч.2.- С.120-123.
4 Kuravsky l.S., Arnautov E.V. On the Approach to Computing Stiffened Structure Natural Modes. - Journal of Sound and Vibration, 1991, v.150, No.1, p.161-166.