Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Иваночкин П.Г. Механика подвижного состава. Учеб пособ. 2019

.pdf
Скачиваний:
1
Добавлен:
15.03.2026
Размер:
4.91 Mб
Скачать

Дифференциальное уравнение движения при гармоническом возбуждении имеет вид:

q 2nq k 2q f

0

sin( pt ).

(4.36)

 

 

 

Решение уравнения (4.36) будем искать в виде суммы общего решения однородного уравнения и частного решения неоднородного уравнения.

Как было показано ранее, общее решение однородного уравнения в зависимости от соотношения между n и k может быть представлено в одной из трех форм:

q

e nt C cos k t C

sin k t при n < k ;

oo

1

1

2

 

1

q e nt C C t

при n = k ;

oo

1

2

qoo

e nt C1e t C2e t

при n > k ,

 

 

 

 

 

где C1, C2 − произвольные постоянные.

Частное решение (4.36) получим, воспользовавшись методом комплексных амплитуд.

Известно, что

f0ei( pt ) f0 cos( pt ) if0 sin( pt ),

где i − мнимая единица, и, следовательно,

f0 sin( pt ) Im f0ei( pt )

Введем вспомогательное уравнение

y2ny k 2 y f0ei( pt )

инайдем его частное решение yчн . Воспользовавшись линейностью вве-

денного уравнения, для которого справедлив принцип суперпозиции, получим qчн как I m yчн . Задав yчн в виде

yчн Gei( pt ) ,

где G − комплексная амплитуда, получим:

(k 2 p2 2nip)G f0 .

Отсюда

G

 

f0

 

f0

.

(k 2

p2 2nip)

D*ei

 

 

 

 

 

71

 

 

 

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2np

 

D*

(k 2 p2 )2

4n2 p2 ;

arctg

.

 

 

k 2

p2

 

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G De i ,

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

f0

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p2 )2

 

 

 

 

 

 

(k 2

4n2 p2

 

 

 

и, следовательно,

yчн Dei( pt ) .

Таким образом,

qчн

Im yчн

Dsin( pt ),

 

 

(4.37)

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

Общее решение (4.36) будет иметь вид:

q e nt C cos k t C

sin k t Dsin( pt )

при n < k ;

1

1

2

1

q e nt C

C t

Dsin( pt )

при n = k ;

1

2

 

 

q e nt C1e t C2e t

Dsin( pt )

при n > k ,

где С1, C2 − произвольные постоянные, определяемые из начальных условий с учетом вида решения.

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

множителя e nt они стремятся к нулю, и в решении остается только qчн .

В этом случае говорят об установившихся вынужденных колебаниях.

На основании решения (4.37) можно сформулировать основные свойства установившихся вынужденных колебаний:

1)эти колебания являются незатухающими; они длятся так долго, как долго действует возмущающая сила;

2)эти колебания не зависят от начальных условий;

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

4)установившиеся вынужденные колебания отстают по фазе от воз-

, изменяющуюся, как будет показано ниже,

.

72

d называют отношение

Амплитуда D и фазовое запаздывание установившихся вынужден-

ных колебаний зависят от соотношения между частотами р и k и коэффициента затухания n . Проанализируем эти зависимости, называемые ам- плитудно-частотной и фазочастотной характеристиками.

Для большей общности результатов перейдем к безразмерным параметрам.

Безразмерным коэффициентом затухания d 2nk .

Если n k и, следовательно, T1 T , то безразмерный коэффициент затухания можно связать с логарифмическим декрементом колебаний:

 

2nT1T 2nT1 T T

 

 

d

 

 

 

 

 

 

.

kTT

2 T

T

 

1

 

1

 

1

 

 

Добротностью Д называют величину, обратную к d:

Д 1d k 2n. .

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

мический декремент колебаний: Д .

Разделив числитель и знаменатель амплитуды D на k 2 , получим:

D Dст ,

где Dст Q0 c – статическое смещение системы от положения равновесия;– коэффициент динамичности,

 

 

 

1

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

(1

z2 )2

d 2 z2

 

 

 

 

(4.38)

z pk – коэффициент расстройки.

Исследуем зависимость коэффициента динамичности от z и d , представляющую собой амплитудно-частотную характеристику системы в

безразмерном виде:

 

1

при z 0;

0

при z ;

1 d Д

при z 1.

Видим, что добротность Д представляет собой коэффициент динамич-

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

73

сравнению с частотой свободных колебаний, т. е. p k , то амплитуда

вынужденных колебаний близка к статическому смещению, а коэффициент динамичности близок к единице. Если же частота изменения возмущающей силы значительно больше частоты свободных колебаний, т. е. p k , то колебательная система ведет себя как фильтр, практически не

воспринимая возмущения с частотами, существенно превышающими собственную частоту. Выражение для коэффициента динамичности (4.38) показывает, что при малых значениях d вязкое сопротивление становится существенным лишь в достаточно узкой зоне в окрестности резонанса, когда d 2 z2 соизмеримо с (1 z2 )2 . Это же демонстрирует график (z) при d 0 ,

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

Рис. 4.5. Зависимость коэффициента динамичности (z) от безразмерного коэффициента затухания d в случае силового возбуждения колебаний

Для определения экстремальных значений коэффициента динамичности достаточно исследовать подкоренное выражение в уравнении (4.38), поскольку его максимум в силу структуры этого уравнения будет соответ-

ствовать минимуму , и наоборот.

 

 

 

 

 

 

 

 

 

 

 

 

Вычислим производные

по

 

z

 

 

от подкоренного выражения

y(z) (1 z2 )2 d 2 z2 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

)

2d

2

z

2z(d

2

2 2z

2

);

y (z) 4z(1 z

 

 

 

 

 

2

2

2z

2

) 8z

2

12z

2

4(1 d

2

2).

y (z) 2(d

 

 

 

 

 

74

Приравняв к нулю y (z) , получим два значения z, соответствующие экстремумам:

z1 0; z2 1 d 2 2.

Второе значение имеет место только при d 2 , а при d 2 оно совпадает с первым.

Если d 2 , то вторая производная y (z) отрицательна при z1 и положительна при z2 . Следовательно, z1 соответствует максимуму y(z) и локальному минимуму (z) , a z2 – минимуму y(z) и максимуму (z) .

Максимальное значение коэффициента динамичности при этом будет равно:

 

 

(z2 )

 

 

1

 

 

 

 

Д

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(d 2 2)2

d 2

 

1 d 2 4

 

 

 

 

 

d 4 2

 

 

 

 

 

 

 

 

Если d 2 , остается только одно экстремальное

значение при

z1 0 . В этом случае имеют место минимум y(z) и максимум (z) .

На рис. 4.5 представлены кривые,

определяющие зависимость (z)

от z при различных значениях коэффициента d . При d

 

 

 

 

2 максималь-

ное значение коэффициента динамичности, равное единице, имеет место при z 0 , т. е. амплитуда установившихся вынужденных колебаний оказывается равной или меньшей статического смещения. При значениях

d 2 максимум (z) оказывается всегда сдвинутым влево от резонанса.

Для исследования фазочастотной характеристики в безразмерном виде разделим числитель и знаменатель аргумента арктангенса в выражении

для на k 2 :

 

 

 

 

arctg

(2n k)( p k)

arctg

dz

.

1 p2 k 2

1 z2

 

 

 

Так как производная по z независимо от значения d (кроме d=0)

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

arctg(0) 0

при z = 0;

 

arctg( ) 2

при z 1;

arctg( 0)

при z .

Отметим, что при резонансе сдвиг по фазе равен /2 вне зависимости от значения коэффициента d, характеризующего вязкое сопротивление.

На рис. 4.6 представлены кривые, характеризующие зависимость (z) при различных значениях d. При d = 0 (отсутствие вязкого сопротивления) в соответствии с полученными выше результатами (см. рис 4.4) функция

75

(z) имеет разрыв. Отметим, что с ростом d изменяется характер фазовой кривой – из кривой с двумя перегибами она трансформируется в кривую с одним перегибом.

В случае инерционного возбуждения колебаний

f0 f0 p2 .

Тогда

Df0 ин

икоэффициент динамичности при инерционном возбуждении

ин

 

 

 

p2

 

 

 

 

z2

 

 

.

 

 

 

 

 

 

 

 

 

(k 2

p2 )2 4n2 p2

 

z2 )2

 

 

 

(1

d 2 z2

Рис. 4.6. Зависимость сдвига фазы от безразмерного коэффициента затухания

этом

k

Для исследования амплитудно-частотной характеристики заменим в выражении коэффициент расстройки z на обратную ему величину p :

ин

 

 

 

 

1 2

 

 

 

 

 

1

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

1

2 )2 d 2

2

 

2 )2

 

 

(1

(1

d 2 2

Полученная зависимость ( ) полностью совпадает по структуре с(z) при силовом и кинематическом возбуждении. Следовательно,

76

 

 

 

ин 0

при

 

z 0, ;

 

 

 

 

ин 1

при

 

z , 0 ;

 

 

 

 

ин 1 d Д

при

z 1, 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

ин Д

1 d 2 4

 

При

d 2

максимальное значение

имеет место

 

 

 

 

 

при 1 d 2 2 или z 1 1 d 2 2 , т. е. в отличие от силового и кинематического возбуждения максимальное значение коэффициента динамичности смещено вправо от резонанса.

При d 2 максимальное значение ин 1при z .

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

ной частоты k .

Рис. 4.7. Зависимость ин (z) от безразмерного коэффициента затухания в случае инерционного возбуждения колебаний

Фазо-частотная характеристика не зависит от способа возбуждения колебаний (см. рис. 4.6).

77

4.2 Колебания системы с двумя степенями свободы

4.2.1 Дифференциальные уравнения малых колебаний консервативной системы в обобщенных координатах

В частном случае системы с двумя степенями свободы квадратичные формы Т, П, Ф будут соответственно равны [12]:

T

1

 

a11q12

2a12q1q2

a22q22 ;

 

 

(4.39)

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П

 

 

1

c q

2

2c q q c q2

 

;

 

(4.40)

 

 

 

 

 

 

 

2

11 1

 

12 1

2 22 2

 

 

 

 

 

 

 

 

 

 

 

 

Ф

 

1

b11q12

2b12q1q2

b22q22 ,

 

 

(4.41)

2

 

 

 

 

 

 

 

 

 

 

 

 

 

а дифференциальные уравнения малых колебаний примут вид:

 

a11q1 a12q2 b11q1 b12q2

c11q1 c12q2 Q1(t);

 

(4.42)

a12q1 a22q2 b12q1 b22q2 c12q1 c22q2 Q2 (t).

 

 

 

Рассмотрим свободные колебания консервативной системы. В этом

случае

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q1 t Q2 t

0; b11 b12 b22 0

 

и дифференциальные уравнения принимают вид:

 

 

 

 

a11q1 a12q2 c11q1 c12q2 0;

 

 

 

(4.43)

a12q1 a22q2 c12q1 c22q2 0.

 

 

 

 

 

 

 

Начальные условия для q1 и q2

имеют вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(0)

 

(4.44)

t 0; q1 (0) q10; q1 (0) q10;

q2 (0) q20; q2

q20.

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

a11 0; a11a22 a122 0; a22 0,

а аналогичные соотношения для квазиупругих коэффициентов

с11 0; c11c22 c122 0; c22 0

являются достаточными условиями устойчивости положения равновесия системы.

Коэффициенты а12 и с12, связывающие в уравнениях (4.43) обобщенные координаты q1 и q2, называют соответственно коэффициентами инерционной и упругой связи. Если в колебательной системе коэффициент

78

a12 0 , ее называют системой с упругой связью, а если c12 0 – системой с инерционной связью.

Парциальной системой, соответствующей обобщенной координате qj, называют условную колебательную систему с одной степенью свободы, получаемую из исходной системы, если наложить запрет на изменение всех обобщенных координат кроме qj. Парциальными частотами называют собственные частоты 1 и 2 парциальных систем:

12

 

с11

; 22

 

c22

.

(4.45)

 

 

 

 

a11

 

a22

 

4.2.2 Общее решение системы дифференциальных уравнений

Поскольку уравнения (4.43) содержат только обобщенные координаты и их вторые производные по времени, ищем их решение в виде:

q1 A1 sin( t ); q2 A2 sin( t ),

(4.46)

где A1, A2 , , – пока не определенные величины.

Подставив (4.46) в (4.43) и приравняв коэффициенты при синусах, получим однородную алгебраическую систему относительно A1 и A2:

 

 

 

a11 A1 c12

 

a12 A2 0;

 

c11

 

 

 

2

 

 

 

2

 

 

 

 

 

 

2

 

A c

2

 

A 0.

(4.47)

 

 

a

a

 

c

 

 

12

 

12

1

22

 

 

22

2

 

Для того чтобы однородная алгебраическая система (4.47) имела ненулевое решение, она должна быть вырожденной, т. е. ее определитель должен равняться нулю:

с 2a

c 2a

 

 

 

0

 

 

 

 

11

11

12

12

 

 

 

 

c 2a

c

 

2a

 

 

 

 

 

(4.48)

22

22

 

 

 

12

12

 

 

 

 

 

 

Следовательно, решение (4.46) будет иметь смысл только при тех значениях ω, которые удовлетворяют условию (4.47). Раскрывая (4.48), получаем:

 

 

 

c11 2a11 c22

2a22 c12

2a12 2

 

0

 

(4.49)

или

a a

 

4

c a

 

 

 

 

 

2

 

 

 

 

 

 

 

a2

c a

2c a

 

c c

 

c2

0.

(4.50)

 

11

22

12

 

11

22

22

11

12

12

 

11

22

12

 

 

79

4.2.3 Частотное уравнение. Собственные частоты, главные колебания. Парциальные частоты

Уравнение, представленное в форме (4.48), (4.49) или (4.50), называют частотным. Как видно из (4.50), частотное уравнение – биквадратное уравнение. Найденные из (4.48)–(4.50) значения называют собственны-

ми частотами колебаний системы.

Исследование корней частотного уравнения позволяет сделать следующие выводы:

1) если положение равновесия устойчивое, то оба корня частотного уравнения положительны;

2) первая собственная частота системы ω1 всегда меньше меньшей

парциальной частоты, а вторая ω2 – больше большей парциальной частоты. Для колебательных систем с упругой связью (a12= 0) справедливо ра-

венство

12 12 22 22.

 

(4.51)

Запишем два частных независимых решения, соответствующих ча-

стотам 1 и 2 , в виде:

 

 

q11 A11 sin( 1t 1), q21 A21 sin( 1t 1);

(4.52)

q12 A12 sin( 2t 2 ), q22 A22 sin( 2

2 ),

 

где вторая цифра в индексе соответствует номеру частоты, или номеру то-

на колебаний.

Константы A11, A21, A12, A22 не являются независимыми, т. к. система (4.47) вырожденная. Коэффициенты связаны между собой соотношениями:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A21 1 A11 ,

 

(4.53)

 

 

c 2a

 

c 2a

 

 

 

 

где

11

1

11

 

12

 

1

 

12

 

,

 

 

1

c

2a

 

c

2a

 

 

 

 

 

12

1

12

 

 

22

 

1

 

22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A22 2 A12 ,

 

(4.54)

 

2

 

c 2a

 

c

2a

 

 

 

 

 

 

где

 

11

2

11

 

12

2

12

 

.

 

 

 

c

2

a

c

 

2

a

 

 

 

 

 

 

 

 

 

 

22

 

 

 

 

 

 

 

12

2

12

 

 

22

2

 

 

 

 

 

 

 

С учетом (4.53) и (4.54) частные решения (4.52) будут иметь вид:

 

 

 

 

 

q11

A11 sin( 1t 1), q21 1 A11 sin( 1t 1),

(4.55)

 

 

 

 

 

q12 A12 sin( 2t 2 ), q22 2 A12 sin( 2

2 ).

 

 

 

 

 

 

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

80