Материал: Нестационарные аэродинамические силы и аэроупругие колебания самолета. Пентюхов В.И., Будник А.П

где Jm - момент инерции масс сверхзвукового самолета как балки.

Из (2.56) следует, что

Согласно (2.55), величина y1n определяет смещение

центра масс сверхзвукового самолета под действием аэродинамических сил, вызванных аэроупругими колебаниями,

определяет угол поворота продольной оси сверхзвукового самолета под действием аэродинамических сил, вызванных аэроупругими колебаниями сверхзвукового самолета,

соответствующая n -ной форме собственных аэроупругих колебаний сверхзвукового самолета как свободной балки, находящейся под действием массовых и аэродинамических сил, вызванных аэроупругими колебаниями сверхзвукового самолета, определяется следующим равенством:

61

m x y2xn x dx

0

n 1,2,3...

(2.58)

С учетом (2.37), равенство (2.58) принимает вид:

формулу Релея для определения собственной частоты свободной балки, колеблющейся в пустоте.

2.5. Определение форм и частот собственных аэроупругих колебаний

Как уже отмечалось, уравнение (2.32), описывающее формы собственных аэроупругих колебаний сверхзвукового самолета как балки, не имеет общих методов решения, даже если пренебречь изменением массы и жесткости по длине сверхзвукового самолета.

В силу этого решение уравнения (2.32), удовлетворявшее граничным условиям (2.33), (2.34), может быть получено только численными методами.

62

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

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

Полученные при произвольном значении частоты как параметра численные решения уравнения (2.32) не будут, в общем случае, удовлетворять граничным условиям в конечной точке.

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

Именно эти значения частоты как параметра и являются собственными частотами аэроупругих колебаний сверхзвукового самолета как балки.

Соответствующие им численные решения уравнения (2.32) описывают формы собственных аэроупругих колебаний сверхзвукового самолета как балки.

Рассмотрим алгоритм решения уравнения (2.32) методом начальных параметров.

Введем дополнительные функции, определяемые следующими соотношениями:

63

Используя равенства (2.60)-(2.63), сведем уравнение (2.13) к следующей системе из четырех дифференциальных уравнений первого порядка:

Используя векторно-матричную форму записи, уравнения (2.64)-(2.67) можно представить в виде

С учетом соотношений (2.62), (2.63) граничные условия

(2.33) (2.34) принимают вид:

64

Как известно, матричное уравнение (2.68), состоящее из четырех дифференциальных уравнений первого порядка, имеет четыре линейно независимых частных решения, так что его решение можно представить в виде

интегрирования уравнения ( 2. 68) при фиксированном значении частоты как параметра и при следующих начальных условиях:

65