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

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

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

где 1, ст

 

m2 m1 mп g

– статическая деформация пружины 1.

 

 

 

 

 

c1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Деформация пружины 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l2 2, ст zк zт .

 

 

 

 

 

 

 

 

Статическая деформация пружины 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2, ст

m2 g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c2

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно выражение потенциальной энергии системы принима-

ет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П m g z

 

m g z

 

m g

z

 

 

с2

 

 

z

 

z

 

2

 

 

 

к

т

т

 

2, ст

к

т

 

 

 

 

2

 

 

1

 

п

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.37)

 

 

 

с1

 

 

 

 

 

2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1, ст

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Частные производные от потенциальной энергии по обобщенным координатам равны:

П с2 zк zт ,zк

 

П

с с

z

 

c z

 

с

 

 

 

т

к

 

1

 

2

 

2

 

1 .

 

 

zт

 

 

 

 

 

 

 

 

 

Диссипативная функция Рэлея в данном случае имеет вид

 

 

 

Ф

1

v1 2

2

v2 2 ,

(5.38)

 

 

2

2

 

 

 

 

 

 

 

 

 

где v1 zп zт – разность скоростей верхней и нижней точек крепле-

ния демпфера 1, и соответственно

v2 zк zт .

Окончательно выражение диссипативной функции принимает вид

Ф

1

2

 

2

2

(5.39)

2

zт

2

zк zт .

 

 

 

 

 

Частные производные от диссипативной функции Рэлея по обобщенным скоростям равны:

Ф 2 zк zт ,zк

Ф 1 2 zт 2 zк 1 .zт

101

С учетом выражений (5.37) и (5.39) уравнения Лагранжа второго рода рассматриваемой системы будут иметь вид:

m z c z z

 

 

z z

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

к т

 

 

 

2

 

к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m m z

m

с с

 

z

т

с z

к

с

2

z

т

 

2

z

к

.

 

1

п т

п

 

 

 

1

 

 

2

 

 

 

2

 

1

 

 

1

 

 

 

 

 

1

 

Перенося члены, содержащие неизвестные из правой части в левую,

получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a11q1 a12q2 b11q1 b12q2 c11q1

c12q2 Q1 (t)

 

 

 

 

 

 

 

 

(5.40)

 

 

 

a12q1 a22q2 b12q1 b22q2 c12q1 c22q2 Q2 (t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, для нашей системы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a11 m2 ,

 

 

a12 a21 0 ,

 

 

 

 

a22 m1 mп ,

 

 

 

 

 

b11 2 ,

 

 

b12 b21 2 ,

 

 

 

 

b22 1 2 ,

 

 

 

 

 

с с ,

 

 

c c

 

c

,

 

 

 

c c c .

 

 

 

 

 

 

 

 

11

2

 

 

12

 

 

 

21

 

 

 

 

 

2

 

 

 

22

1

 

2

 

 

 

 

 

 

 

 

Q1(t) 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q2 (t) mп 1 с1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Парциальные частоты 1 ,

2

для рассматриваемой системы будут

равны:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

c11

 

 

c2

 

 

 

 

 

2

 

с22

 

c1 c2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

;

 

 

2

 

 

 

 

 

.

 

 

 

(5.41)

 

 

 

 

a

 

 

m

 

 

 

 

 

a

m m

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

2

 

 

 

 

 

22

 

1

 

 

п

 

 

 

 

 

 

Для железнодорожных экипажей обычно ν21, т. к. m2>m1+mn. Для определения собственных частот колебаний экипажа составим частотное уравнение

 

 

 

 

 

 

 

c

k 2a

 

c

k 2a

0 .

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

11

12

12

 

 

 

 

 

 

 

 

 

 

 

c

k 2a

 

c

k 2a

 

 

 

 

 

 

 

 

 

 

 

 

 

21

 

 

 

21

22

22

 

 

 

 

 

 

 

 

В нашем случае частотное уравнение имеет вид

 

 

 

 

 

 

 

 

 

c c k 2

m m c k 2m

с2

0

 

 

 

 

 

 

 

1

 

2

 

 

1

п

2

2

2

 

 

 

или

 

 

 

 

 

 

k 4 12

22 k 2 12 22 12 32 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

(5.42)

где

2

 

 

c2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

m

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

п

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

2

2

2

2

2 2

 

 

 

 

 

 

 

 

 

k1,2

 

 

 

1

2

1

2

 

4 1 3

 

.

(5.43)

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В соответствии с тем, что ν21, из формулы видно, что k1<ν1, k2>ν2. Меньшую из частот k1 называют основной частотой.

Рассмотрим собственные колебания этой системы при отсутствии сил вязкого сопротивления гасителей. Положим, bij=0 (i, j = 1,2), а правые части уравнений (5.40) приравняем нулю:

102

̈+ − = 0;

2 ̈

2

2

( + ) − + ( + ) = 0.

1

2

1 2

Частными решениями однородной системы (5.44) будут:

zк1

А sin k1t 1 ,

zт1 1А sin k1t 1 ,

 

 

B sin k t

 

,

 

 

 

B sin k t

 

.

z

 

 

z

 

 

 

к2

2

2

 

 

т2

2

2

2

 

(5.44)

(5.45)

Колебания, совершаемые системой с одной из собственных частот k1 или k2, называются главными колебаниями. Главное колебание, соответствующее меньшей собственной частоте k1, называют основным колебанием, т. к. это колебание является основным в результирующем движении системы.

Коэффициенты 1 и 2

 

c k2а

 

2

k2

1

 

k2

0,

 

 

11

1

11

 

 

 

1

1

 

 

 

 

1

 

(5.46)

с k2а

 

 

 

2

 

 

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12

1

12

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

c k2а

 

2 k

2

1

 

k2

0 .

 

 

11

2

11

 

1

2

 

2

(5.47)

 

с

k2а

 

 

2

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

12

2

12

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

называют коэффициентами распределения амплитуд. Они характеризуют отношение амплитуд в главных колебаниях или форму главных колебаний

(коэффициент формы).

Из (5.46)–(5.47) следует, что в первом главном виде колебаний системы (k k1 ) перемещения всегда совпадают по знаку, т. е. обе массы совершают колебание в одной фазе, а во втором главном виде (k k2 ) –

впротивофазах (навстречу друг другу).

Вобщем случае собственные колебания кузова и тележек вагона определяются зависимостями:

zк

Asin(k1t 1) Bsin(k2t 2 );

(5.48)

zт

A 1 sin(k1t 1) B 2 sin(k2t 2 ),

 

где амплитуды A1, A2 и начальные фазы α1, α2 определяются из начальных условий. Периоды (в секундах) колебаний основного (с наименьшей частотой k1) тона Т1 и наложенного Т2 определяются из выражений:

T1 2 ; T2 2 . k1 k2

Частоты колебаний n выражены числом колебаний в секунду (Гц) в отличие от угловых частот k (рад/с):

103

n1

1

 

2

;

n2

1

 

2

.

(5.49)

T1

 

T2

 

 

 

k1

 

 

k2

 

Для современного груженого четырехосного вагона и пути с деревянными шпалами имеем:

mк = 80 т; 2 m т = 14 т; 4 cz = 16 МН/м; 8 cп = 400 МН/м.

Парциальные частоты 1 14,1 с-1, 2 = 66,5 с-1;

главные частоты k1 = 13,8 с-1 (2,2 Гц), k2 = 172 с-1 (27 Гц). Следовательно, частота второго вида колебаний k2 более чем на

порядок выше частоты k1 первого вида.

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

Вынужденные колебания масс рассматриваемой системы характеризует система дифференциальных уравнений (5.40) при условии, что0 . Как и ранее, будем пренебрегать силами сопротивления гасителей

( bij 0 i, j 1, 2 ):

m2 zк с2 zк с2 zт 0;

 

(5.50)

(m1 mп )zт с2 zк (с1 с2 )zт mп с1 .

 

Для упрощения вычислений положим, что неровность рельсовых

нитей может быть выражена как

 

0 sin t,

(5.51)

где 2L v – частота чередования неровностей пути с длиной волны L

при движении вагона со скоростью v.

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

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

zк

C sin t ;

(5.52)

zт

D sin t .

 

104

Подставив выражения (5.52) в уравнения (5.50) и приравнивая коэффициенты при синусах, получим:

 

C( 2 2 ) D 2 0;

 

 

с

 

c c

 

c m 2

 

1

1

 

12

2

, 22

 

1 2

,

u

1 п

.

 

C 12 D( 2 22 ) u 0 ;

 

 

 

 

 

 

 

m2

 

m1 mп

 

m1 mп

Из последних выражений найдем:

 

 

 

 

 

 

 

 

 

 

C

u 2

u ( 2 2 )

 

22 .

 

 

(5.53)

( 2 12 )( 2

22 ) 12 22 ; D

( 2 12 )( 2

22 ) 12

 

 

 

0

1

 

0

1

 

 

 

 

 

 

 

 

 

Следовательно, частные решения уравнений (5.50) получат вид:

 

zк 0 к sin t;

zт 0 т sin t,

 

 

 

 

 

(5.54)

где к , т – коэффициенты нарастания амплитуд (коэффициенты динамичности),

к

 

 

 

 

u 2

 

 

; т

 

 

u( 2

2 )

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

.

(

2

2

)(

2

 

2

2 2

(

2

2

)(

2

 

2

2 2

 

 

 

1

 

2 ) 1 2

 

1

 

2 ) 1 2

Полное решение системы уравнений (5.40) с правой частью получим суммированием общих решений (5.48) для однородной системы с полученным частным решением (5.54) для неоднородной системы.

Учитывая, что собственные колебания исследуемой системы затухнут при условии, когда движение по пути с непрерывной неровностью будет длительным, практический интерес представляет изучение только вынужденных колебаний системы, характеризуемой выражениями (5.54). Амплитуды этих колебаний 0 к и 0 т изменяются по закону изменения ко-

эффициентов нарастания к , т в зависимости от частоты или от скоро-

сти v

При весьма малых скоростях и больших длинах волн неровностей величина мала по сравнению с 1 и 2 . Тогда из уравнения (5.53) с учетом выражения (5.41) получим:

С D

 

u 0

0 .

(5.55)

2

2

 

 

 

 

1

2

 

 

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

Сравнивая знаменатели выражений (5.53) с частотным уравнением (5.42), нетрудно видеть, что значения С и D будут весьма большими приk1 и k2 . Значит, существуют две резонансные частоты 1 и

2 и соответственно две критические скорости v1 и v 2 , при которых амплитуды вынужденных колебаний масс mk и mm неограниченно увеличиваются.

105

Отношение амплитуд колебаний тележек и кузова вагона будет:

D

 

2

2

1

2

 

 

 

 

1

 

 

.

(5.56)

C

 

 

2

2

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

Из выражения (5.56) видно, что при скоростях движения, для которых 1, D C 0, т. е. величина D мала по сравнению с величиной С. Для

грузовых вагонов в груженом состоянии это будет иметь место при движении с эксплуатационными скоростями v = 20−25 м/с по пути с длиной волн неровностей:

L

2

v

2

 

 

 

(20−25) = 9−11 м.

 

1

14.1

 

 

 

 

 

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

5.6. Колебания кузова на рессорах с сухим трением и упругофрикционными связям

Дифференциальное уравнение подпрыгивания кузова на рессорах с постоянным сухим трением (рис. 5.7) имеет вид:

mк z cz Fтр 0

(5.57)

Здесь знак плюс соответствует этапу движения, на котором скорость движения положительна, а знак минус − этапу движения, на котором скорость отрицательна.

Часто уравнение (5.57) записывается в виде

mк z cz Fтр sign z 0.

(5.58)

106

Z

m

с

Рис. 5.7. Схема для исследования колебаний кузова на рессорах с постоянным сухим трением

В последней записи использована функция сигнум, которая определяется формулой:

 

1

 

при z 0;

 

0

при z 0;

sign z

 

 

при z 0.

1

Уравнение (5.58) содержит нелинейное слагаемое. Тем не менее решение этого уравнения можно получить, если рассмотреть последовательные интервалы движения, на каждом из которых скорость z имеет постоянный знак (метод припасовывания). Пусть начальное условие задачи следующее:

при t=0 z z0 ,

z 0 .

(5.59)

Под действием силы упругости пружины кузов будет двигаться внизz 0 и уравнение движения имеет вид:

 

 

 

 

mк z cz Fтр

0

(5.60)

или

 

 

 

 

 

 

 

 

 

 

z 2 z 2 f

тр

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

с

 

 

 

 

 

 

 

– собственная частота колебаний;

 

m

 

 

 

 

 

 

 

 

107

fтр

Fтр

 

 

– стрела трения (прогиб рессор, который соответствует силе

c

 

 

Fтр ).

 

 

При отклонении кузова z0 на величину, меньшую или равную f тр ,

движение не начинается, т. к. силы упругости пружины недостаточно для преодоления силы трения. Полоса fтр z fтр называется зоной застоя

или порогом чувствительности рессорного подвешивания. Общее решение уравнения (5.60) имеет вид:

z fтр c1 cos t c2 sin t .

Определяя постоянные интегрирования из начальных условий, полу-

чим:

z fтр z0 fтр cos t .

(5.61)

Закон движения (5.61) справедлив до тех пор, пока

z 0 . Так как

 

 

z z0 fст sin t , то скорость движения будет отрицательной до мо-

мента времени t1, определенного из условия t1 . В этот момент кузов остановится, а смещение z при этом равно:

z fтр z0 fтр cos z0 2 fтр .

После остановки кузов начнет двигаться вверх. Повторяя приведен-

ные выше расчеты, можно сказать, что движение вверх также будет про-

должаться в течение времени . Максимальное отклонение вверх равно

z0 4 fтр . Процесс движения будет продолжаться до тех пор, пока кузов не остановится в зоне застоя. Зависимость смещения z t (рис. 5.8) на каждом этапе движения представляет собой косинусоиду, смещенную по оси z

на величину fст или fст , с амплитудой,

уменьшающейся по закону

арифметической прогрессии:

 

 

ai 1 i z0 fтрi

,

(5.62)

где i = 0, 1, 2, 3,… – номера точек касания графика колебаний огибающими прямыми;

T 2 – время между двумя соседними максимумами отклонения, ко-

торое условно можно назвать периодом колебаний. Наличие сухого трения не меняет частоту колебаний.

108

Рис. 5.8. График затухания собственных колебаний кузова на рессорах с постоянным сухим трением

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

F c 1 sign ,

где с – жесткость рессорного комплекта;

φ – коэффициент относительного трения, равный отношению силы Fтр

внутреннего трения упругого элемента к силе Р, создающий упругую деформацию;

fст z – деформация упругого элемента;

fст Pcст – статический прогиб, вызванный нагрузкой экипажа брутто

Pст ;

z – дополнительный прогиб, отсчитываемый от положения статического равновесия;

̇– скорость деформации.

109

а

б

Рис. 5.9. Упруго-фрикционная связь с силой трения, пропорциональной перемещению:

а– график силовой характеристики;

б– схематическое изображение этой связи

Жесткость рессор при сжатии c1 c 1 больше, чем при отпуске

c2 c 1 .

 

Дифференциальное уравнение колебаний подпрыгивания кузова на

таких рессорах будут иметь вид:

z 1 sign z 0 .

mк z c fст

 

 

Легко видеть, что это уравнение также может быть решено методом «припасовывания». На первом интервале движения уравнение можно записать в виде:

 

 

 

 

 

 

 

 

 

 

 

 

z 12 z 12 fст

, 1

 

c 1

.

 

 

 

 

 

 

 

 

 

 

 

mк

Если использовать те же начальные условия (5.59), то получаем ре-

шение уравнения:

 

z fст z0

fст cos 1t .

 

 

В момент t1

 

скорость кузова z 0, а смещение a1 z0 zfст .

 

 

1

 

 

 

 

 

 

 

 

Для следующего интервала движения уравнение имеет вид:

 

 

 

 

 

 

 

 

z 22 z 22 fст

, 2

c 1

,

 

 

 

 

 

 

 

 

mк

а его решение:

z fст z1 fст

cos 2t .

 

110