Изгибные колебания. Получение системы дифференциальных уравнений
Допущения:
- лопатку рассматриваем как закрученный стержень переменного сечения;
- отсутствует связь изгибных и крутильных колебаний;
- ось лопатки проходит через центры тяжести всех сечений и нерастяжима;
- при колебаниях все сечения остаются плоскими и перпендикулярными упругой линии (гипотеза Кирхгофа – Лява).
Примем систему координат OXYZ, начало которой O расположено в центре тяжести корневого сечения, оси X и Y совпадают с главными центральными осями инерции корневого сечения (пусть ), ось Z направлена по оси лопатки. Оси вспомогательной системы координат параллельны осям главной системы координат, начало находится в центре тяжести текущего сечения. Выделим в лопатке бесконечно малый элемент высотой dz (рис. 14.1).
Рис.14.1. Бесконечно малый элемент лопатки
и динамические силы, приложенные к нему при изгибных колебаниях
Рассмотрим условия равновесия элемента с учетом сил инерции
и
,
где и - упругие смещения в направлении осей X и Y. Знак волны над переменной означает то, что она зависит от времени.
В направлении оси
, откуда
(14.1)
Из условия равновесия в направлении оси аналогично получаем
(14.2)
Вокруг оси :
.
Приводя подобные и пренебрегая бесконечно малыми второго порядка, получаем
. (14.3)
Вокруг оси :
.
Приводя подобные и пренебрегая бесконечно малыми второго порядка, получаем
. (14.4)
Установим связь моментов и углов поворота сечений. В плоскости OXZ (рис. 14.2)
Рис.14.2. Изгиб оси лопатки в плоскости OXZ
1 – положение сечения в равновесном состоянии; 2 - положение изогнутой оси лопатки;
3 - касательная к изогнутой оси лопатки в месте расположения сечения.
Из рис. 14.2 видно, что
С учетом знаков момента и угла поворота
(14.5)
Рис.14.3. Изгиб оси лопатки в плоскости OYZ
1 – положение сечения в равновесном состоянии; 2 - положение изогнутой оси лопатки;
3 - касательная к изогнутой оси лопатки в месте расположения сечения.
В другой плоскости из рис. 14.3 видно
(14.6) .
Сечение лопатки поворачивается вокруг осей X и Y. Из рис. 14.4 видно, что (с учетом знаков углов поворота) перемещение произвольной точки сечения в направлении оси Z вследствие поворота сечения вокруг оси X составляет
,
Рис.14.4. Поворот сечения лопатки в плоскостях OXZ и OYZ
вследствие поворота сечения вокруг оси Y составляет
.
Есть также перемещение в направлении оси Z сечения в целом, связанное с перемещением его центра, точки О1 . Полное перемещение точки в направлении оси Z
.
Согласно формулам Коши
.
Тогда напряжение в окрестности точки
.
Это напряжение создает бесконечно малые моменты вокруг осей (рис. 14.5)
Рис. 14.5. Определение знаков моментов сил от напряжения на бесконечно малой площади dF
,
.
Полный момент
В первое слагаемое входит статический момент вокруг оси . Поскольку – главная ось инерции, этот момент равен нулю. Во второе слагаемое входит момент инерции относительно оси , в третье – центробежный момент инерции. Таким образом,
Аналогично:
Поскольку деформации лопатки при колебаниях малы, будем считать, что оси главной и вспомогательной систем координат приблизительно параллельны. Тогда
; , моменты инерции тоже приблизительно одинаковы.
Запишем два последних уравнения в матричной форме.
.
Обозначив матрицу как [A] , получим
.
Пусть матрица имеет вид . Тогда
(14.7)
(14.8)
Уравнения (14.1)…(14.8) образуют систему дифференциальных уравнений, описывающих изгибные колебания лопатки:
Параметры в первых четырех уравнениях, а именно: qx ; βy ; My ; Qx имеют отношение к плоскости OXZ, поэтому первые четыре уравнения описывают колебания в плоскости OXZ, а последние четыре – определяют колебания в плоскости OYZ. Связь между колебаниями в этих плоскостях осуществляется через коэффициент axy при Mx и My. Если лопатка не слишком закрученная, то оси и приблизительно совпадают с главными осями инерции во всех сечениях, в таком случае центробежные моменты инерции близки к нулю, коэффициент axy = 0 и связь между колебаниями в плоскостях OXZ и OYZ отсутствует. В этом случае
; .
Рассмотрим колебания в плоскости OXZ, т. е. только первые четыре уравнения. Их решение будем искать в гармоническом виде (где p – собственная частота колебаний):
После подстановки в дифференциальные уравнения получаем:
Из первого уравнения . Подставим это во второе уравнение, получим . Подставив в третье уравнение, получим . Подставив в четвертое уравнение, получим . Окончательно получили дифференциальное уравнение четвертого порядка
.
Обозначим
. Тогда
,
.
Решение уравнения четвертого порядка - сумма четырех линейно независимых частных решений.
Частными решениями уравнения являются , , а также гиперболические косинус и синус
В решение будем подставлять их суммы, которые называются функциями Крылова.
,
,
,
.
Преимущество этих функций состоит в том, что при дифференцировании они переходят друг в друга: dS/dz = V ; dT/dz = S ; dU/dz = T; dV/dz = U.
При этом S(0)=1; T(0) = U(0) = V(0) = 0.
Подставим решение вида в полученные уравнения.
,
,
.
Неизвестные постоянные С1, С2, С3, С4 определяются из граничных условий.