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

Решение системы уравнений (3.12.) ищем в виде:

Z1 = a1sinωt;   Z2 = a2sinωt;   Z3 = a3sinωt,(3.13)

где а1, а2, a3- амплитудные коэффициенты; ω - частота колебаний.

Подставляя (3.13) в (3.12), получаем систему однородных алгебраических уравнений.

(r11 - ω2m1) a1 + r12a2 + r13a3 = 0
r21a1 + (r22 - ω2m2)a2 + r23a3 = 0
r31a1 + r32a1 + (r33 - ω2m3)a3 = 0
(3.14)

54

Вычисляя коэффициенты системы канонических уравнений метода перемещений rik, имеем

r1 = r;  r22 = r33 = 2r;   r12 = r23= - r;   r13 = 0,(3.15)

где r = 4i/3=16EJ/9.

Введем параметр λ = ω2т/r.

Тогда с учетом принятых обозначений система уравнений (3.14) может быть представлена в виде:

(1 - λ) a1 - a2   = 0
-a1 + 2 (1 - λ) a2 - a3 = 0
  - a2 + (2 - 3λ) a3 = 0
(3.16)

Как обычно, для вычисления собственных значений системы однородных уравнений (3.16) используем условие равенства нулю определителя, составленного из коэффициентов этой системы [27]:

(1 - λ) -1 0
-1 2(1 - λ) -1
0 -1 (2 - 3λ)
 = 0
(3.17)

Раскрывая этот определитель, получаем кубическое уравнение относительно искомого параметра λ:

3 -16λ2 + 10λ - 1 = 0.(3.18)

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

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

Для нахождения первого корня λ1 уравнения (3.18) применяем итерационный метод. Первое приближение λ1(1) принимаем в виде отношения последних коэффициентов уравнения (3.18)

λ1(1) = 1/10 = 0,1.

Подставляя это значение λ = λ1(1) в уравнение (3.18) и обозначая его левую часть через A= 6λ3-16λ2+ 10λ-1=0, находим А (1) = -0,154.

Для второй итерации принимаем λ1(2) =0,15. В этом случае имеем

А(2) = 0,160. Очевидно, истинное значение корня λ1 заключено в интервале 0,1 ≤ λ1 ≤ 0,15.

55

Продолжая итерационный процесс, в пятом приближении получаем λ1(5)= 0,123; при этом A (5) = -0,001.

Понятно, что продолжая этот итерационный процесс, можно обеспечить сколь угодно высокую точность вычисленияλ1.

Принимая λ = 0,123, разделим уравнение (3.18) на множитель (λ - 0,123). В таком случае получаем квадратное уравнение

2 - 15,262λ + 8,123 = 0,

корни которого легко вычисляются:

λ2 = 0,758;                               λ3=1,786.

Для проверки правильности вычисления корней уравнения (3.18), приводим определитель (3.17) к виду

|D-λE| =
(1 - λ) -1 0
-1
2
(1 - λ)
-1
2
0
-1
3
(2 - λ)
3
 = 0

где Е - единичная диагональная матрица.

Здесь должно выполняться условие λi = SD , где

∑λi = 0,123 + 0,758 + 1,786 = 2,667

- сумма корней уравнения (3.18);

SD = 1 + 1 + 2/3 = 2,667

- след матрицы D, равный сумме диагональных элементов этой матрицы.

Возвращаясь к введенным выше обозначениям параметров λ и r, вычисляем частоты собственных колебаний по формуле

ωi
 λi
m
(3.19)
где E=2,7·107кН/м2; J = bh3 = 0,3·0,33 = 0,000675м4
12 12




Имеем:

ω1 = 6,313 сек-1;   ω2= 15,67 сек-1;   ω3= 24,05 сек-1.

Отсюда находим периоды колебаний: Ti = 2π/ωi:

T1 = 0,995 сек;   Т2 = 0,401 сек;   Т3 = 0,261 сек;

56

Определение форм собственных колебаний

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

Находим:

V1 = { 1 }; V2 = { 1 }; V3 = { 1 },
0,877 0,242 -0,786
0,538 -0,883 0,236
(3.20)

где Zi- векторы, определяющие i-ую форму колебаний.

Проверяем ортогональность форм собственных колебаний по формуле(2.42)

 k=1n  mkVikVjk = 0

Например, для первой и второй форм имеем:

т (1·1·1 + 2·0,877·0,242 – 3·0,538·0,883) = 0,0007≈ 0.

Аналогично проверяется ортогональность первой и третьей, второй и третьей форм.

57

Lib4all.Ru © 2010.
Корпоративная почта для бизнеса Tendence.ru