Иваночкин П.Г. Механика подвижного состава. Учеб пособ. 2019
.pdf
Дифференциальное уравнение движения при гармоническом возбуждении имеет вид:
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 и фазовое запаздывание установившихся вынужден-
ных колебаний зависят от соотношения между частотами р и k и коэффициента затухания n . Проанализируем эти зависимости, называемые ам- плитудно-частотной и фазочастотной характеристиками.
Для большей общности результатов перейдем к безразмерным параметрам.
Безразмерным коэффициентом затухания d 2n
k .
Если 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 p
k – коэффициент расстройки.
Исследуем зависимость коэффициента динамичности от 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
