Метод разложения по собственным формам колебаний применительно к конечноэлементным моделям




В конечноэлементной трактовке задача о нагружении сводится к решению системы обыкновенных дифференциальных уравнений относительно вектора независимых обобщённых координат Х

(7.5.1)
с начальными данными , (7.5.2)

где , , – матрица масс, жёсткостей и диссипативных свойств системы, – вектор внешних сил.

Среди методов определения динамических нагрузок, используемых в РКТ, наиболее широкое распространение получили метод разложения по собственным тонам колебаний и метод прямого интегрирования уравнений движения.

При использовании метода разложения по собственным тонам колебаний задача о динамическом нагружении решается в два этапа. Сначала определяются собственные частоты и формы колебаний. Затем независимо от вида модели (балка, оболочка, упругое тело) решение отыскивается в виде разложения по собственным тонам колебаний. С помощью обобщённого метода Галеркина либо вариационного принципа Лагранжа получаются обыкновенные дифференциальные уравнения для обобщённых координат.

Воспользуемся для исследования вынужденных колебаний конечно-элементного уравнения (7.5.1) методом разложения по собственным формам колебаний. Соответствующая задача о собственных колебаниях формируется из (7.5.1) путём отбрасывания диссипативных сил и правых частей уравнения. При этом решение отыскивается в виде стационарных колебаний

(7.5.3)

В итоге приходим к следующей задаче о собственных значениях

(7.5.4)

В результате решения этой задачи определяются собственные частоты и собственные формы колебаний . При этом собственные формы колебаний определяются с точностью до констант в виде множителей в силу однородности задачи (7.5.4).

Размерность вектора совпадает с размерностью вектора Х (с количеством независимых обобщённых координат задачи N 1). Максимальное количество тонов колебаний при решении этой задачи также совпадает с размерностью задачи N 1 (размерностью вектора Х).

Однако достоинством метода разложения по тонам колебаний является то, что при его использовании можно ограничиться гораздо меньшим количество тонов по сравнению с размерностью задачи. Спектр учитываемых в решении тонов колебаний должен увязываться со спектром внешних воздействий и контролироваться характером сходимости решения.

Выразим вектор кинематических параметров Х в (7.5.1) через вектор обобщённых координат T (главных координат) и матрицу собственных векторов , учитываемых в решении задачи:

(7.5.5)

Столбцами матрицы являются l собственных векторов, учитываемых в решении задачи (7.5.1).

совокупность l векторных равенств, отвечающих различным собственным векторам и собственным частотам (7.5.4), можно записать в виде следующего матричного равенства:

(7.5.6)

Матрица имеет порядок где – порядок вектора Х, а l - количество учитываемых в решении тонов колебаний, которое может быть гораздо меньше размерности задачи .

Матрица имеет диагональную структуру

(7.5.7)

Подставив Х из (7.5.5) в уравнение движения (7.5.1) (предварительно опустив в нем диссипативный член) и учитывая равенство (7.5.6), будем иметь

(7.5.8)

Умножив равенство (7.5.8) слева на матрицу , получим уравнение движения в обобщённых координатах

. (7.5.9)

Здесь учтена ортогональность собственных форм колебаний и их нормировка по единичной приведенной массе

(7.5.10)

В равенстве (7.5.10) – единичная матрица l -го порядка (l - количество учитываемых тонов колебаний).

Чтобы система уравнений была распадающейся, введем диссипативные силы, пропорциональные скоростям обобщённых координат. Тогда уравнение движения примет следующий окончательный вид:

(7.5.11)

где матрица диссипативных сил:

(7.5.12)  

Начальные данные для интегрирования системы (7.5.11) определяются на основе заданных полей перемещений и скоростей (7.5.2) с использованием соотношения (7.5.5):

(7.5.13)

Умножая эти равенства слева на с учётом (7.5.10), придём к начальным значениям для обобщённых координат

, . (7.5.14)

В процессе интегрирования системы (7.5.11) c начальными данными (7.5.14) вычисляются обобщённые координаты (абсолютные перемещения и углы поворота), а затем векторы силовых факторов и напряжения в элементах.



Поделиться:




Поиск по сайту

©2015-2024 poisk-ru.ru
Все права принадлежать их авторам. Данный сайт не претендует на авторства, а предоставляет бесплатное использование.
Дата создания страницы: 2019-06-03 Нарушение авторских прав и Нарушение персональных данных


Поиск по сайту: