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

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

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

 

 

 

 

 

 

 

 

 

 

Проекция на ось Оz силы упругости F2

2-й пружины, действующей

на тело 2, равна:

упр.2 / z2

 

 

 

 

 

F2z

c2 ( ст.2

z2 z1 ) .

На тело 1 со стороны 2-й пружины действует такая же по величине,

 

 

 

 

 

 

 

 

 

 

но противоположная по направлению сила F *

F , ее проекция на ось

 

 

 

 

 

 

 

 

2

2

Оz равна:

 

 

 

 

 

 

 

 

 

F *

F

c (

ст.2

z

2

z ) .

 

 

2z

2z

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

Проекция на ось Оz силы упругости F1

1-й пружины, действующей

на тело 1, равна:

 

 

 

 

 

 

 

 

 

F1z

упр.1 / z1

c1 ( ст.1

z1 h(x)) .

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

 

 

 

 

 

1 2 1 vz1 2 / 2 2 vz 2 2 / 2 ,

(5.9)

где

v

 

z

v

 

z

 

dh(x)

v

 

– разность скоростей верхней и нижней

z1

Kz

 

0

 

 

1

 

1

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

точек крепления демпфера 1, установленного между телом 1 и колесом K,

взятая в проекции на вертикальную ось;

 

 

vz 2

z2

z1

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

ния демпфера 2, установленного между телами 1 и 2, в проекции на вертикальную ось.

После подстановки vz1 и vz 2 получаем выражение для Ф:

 

 

 

 

 

 

 

 

 

dh(x)

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

z

 

 

v

 

 

 

/ 2

 

z

 

z / 2 .

(5.10)

 

 

 

 

 

 

 

 

 

1

 

2

1

 

1

 

 

dx

 

 

0

 

 

 

 

 

2

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проекция на ось Оz силы вязкого сопротивления R2 2-го демпфера,

действующей на тело 2, равна:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2 z 2 2 (z2 z1 ) .

 

 

 

 

 

 

 

 

 

 

 

 

z2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На тело 1 со стороны второго демпфера действует такая же по вели-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

чине, но противоположная по направлению сила

R*

 

R , ее проекция на

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

ось Оz равна:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

*

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2 z

R2 z 2 (z2

z1 ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проекция на ось Оz силы сопротивления

R1

1-го демпфера, дей-

ствующей на тело 1, равна:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

dh(x)

 

 

 

 

 

 

 

 

 

 

 

R1z

 

 

1

z1

 

 

 

v0

.

 

 

 

 

 

 

 

z1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

91

Используя выражения (5.7), (5.8), (5.9), получаем уравнения Лагранжа 2-го рода для рассматриваемой системы:

m1 z1

(c1 c2 )z1 ( 1 2 )z1 c2 z2 2 z2

c1h(x) 1

dh(x)

v0 ,

(5.11)

 

 

 

 

dx

 

m2 z2 c2 (z2 z1 ) 2 (z2 z1 ) .

Начальные условия движения:

при t = t0, z1(t0 ) z1(0) , z1(t0 ) z1(0) , z2 (t0 ) z2(0) , z2 (t0 ) z2(0) .

Найдем положение равновесия данной механической системы (см. рис. 5.2), и исследуем его устойчивость.

Потенциальная энергия системы (5.8) имеет вид:

Π Π(0)+

1

(c1 c2 )z12

2 ( c2 ) z1z2 c2 z22 .

2

 

 

 

 

 

Вычисляем частные производные от потенциальной энергии по

обобщенным координатам:

 

 

 

 

 

Π (c1 c2 )z1 c2 z2 ,

 

Π

c2 z1 c2 z2 .

 

 

z1

 

 

 

z2

Условия равновесия (3.5) принимают вид:

(c1 c2 )z1 c2 z2 0 ,c2 z1 c2 z2 0 .

Положение равновесия получаем, решая эту систему линейных однородных алгебраических уравнений: z1 = 0, z2 = 0.

Исследуем найденное положение равновесия на устойчивость, применяя критерий Сильвестра.

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

 

с11 =

 

c1

c2 ,

 

z

 

 

 

1

 

 

 

с12

= z z

2

c2 ,

с21 = с12,

 

1

 

 

 

 

 

с22

= z

c2 .

 

 

 

 

2

 

 

Согласно (3.12) для устойчивости должны выполняться неравенства:

1 c11 = с1 + с2 > 0,

 

 

c12

 

 

 

c1 c2

c2

 

= с1с2 > 0,

2

 

c11

 

 

 

 

 

c21

c22

 

 

 

c2

c2

 

 

откуда следует, что положение равновесия устойчиво.

92

5.3 Исследование продольных колебаний системы «локомотив – состав»

Рассмотрим показанную на рис. 5.3 механическую систему, состоящую из тел А и В, массы которых равны соответственно m1 и m2, расположенных на горизонтальном участке пути и соединенных между собой

сцепкой, которая имеет коэффициент жесткости с и коэффициент вязкого

сопротивления . Будем считать, что к первому телу приложена сила

 

, а

F1

ко второму

 

 

 

 

 

 

 

 

 

 

 

F2 . Эти силы, направленные по горизонтали, в общем случае

являются переменными величинами, представим

 

 

 

 

и

их как F1 F1 (t)e1

 

 

 

 

 

 

 

 

 

 

 

 

 

F2

F2 (t)e1

, где

 

e1

– орт горизонтальной оси Ох. Требуется определить

движение тел системы.

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r1

rЦ

 

 

m1

 

c

 

 

m2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e1

A

 

F1

Ц

 

 

B

 

F2

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

O

Рис. 5.3. Двухмассовая модель системы «локомотив – состав»

Положение тел, размерами которых пренебрегаем, определяется их

радиус-векторами r1 и r2 (рис. 5.3). Рассматриваемая система имеет две

степени свободы. Абсциссу тела А обозначим как х1, тела В х2. Обобщенные координаты q1 и q2 выберем таким образом, чтобы

 

 

 

 

 

 

 

x1 q1 ,

 

 

x2 q1 q2 ,

 

т. е. q1 определяет положение тела А, а q2

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

А и В.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(q1

 

 

 

 

 

r1

x1e1 q1e1

,

 

r2

x2e1

 

q2 )e1 .

(5.12)

Перейдем к составлению уравнений Лагранжа. Кинетическая энер-

гия системы запишется как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T T

 

T

 

 

1

m q 2

 

1

m

 

(q

q

 

)2

 

 

A

B

 

 

2

2

 

 

 

 

 

 

 

 

2

 

1 1

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5.13)

 

 

1

(m m

 

)q 2

2m q

q

 

 

m

q 2

.

 

 

 

 

 

 

 

2

2

 

 

 

2

 

 

1

 

 

1

 

2

1

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

93

Остановимся на вопросе классификации сил, действующих на систе-

му. К внешним силам относятся: силы тяжести и силы реакции гладкой го-

 

 

 

 

 

 

 

 

 

 

 

 

ризонтальной плоскости, направленные по вертикали, а также силы

F1 и

 

 

 

 

 

 

 

 

 

 

 

 

F2 . К внутренним силам относятся силы упругости пружины и силы вязко-

го сопротивления демпфера, приложенные к телам А и В.

 

Длина пружины, соединяющей тела системы, равна в текущем поло-

жении l = q2, деформация пружины

 

 

 

 

 

 

 

 

l = l l0 = q2 l0,

 

 

где l0 – длина пружины в недеформированном состоянии.

 

Потенциальная энергия пружины

 

 

 

 

 

 

 

 

 

 

1

c( l)2

 

1

c(q l )2 .

(5.14)

упр

 

 

 

2

2

2

0

 

 

 

 

 

 

 

 

Обобщенные силы, соответствующие Пупр, находим как

 

Q1П упр / q1

 

0,

 

 

Q

упр / q2 c(q2

l0 ) .

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

 

 

 

1

v 2 ,

 

 

 

 

(5.15)

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где разность абсолютных скоростей тел А и В

v q2 .

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

Q1 / q1 0 ,

Q2 / q2 q2 .

(5.16)

С учетом вышеизложенного уравнения Лагранжа 2-го рода (2.22) для изучаемой системы «локомотив – состав» принимают вид:

(m1 m2 )q1 m2q2

F1 (t) F2

(t),

(5.17)

m2q1 m2q2 c(q2 l0 ) q2 F2 (t).

 

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

x

x ,

x v ,

при t = 0, 1

 

 

10

1 10

x2 x20 ,

x2 v20 ,

q q x ,

q q v ,

 

 

 

 

 

 

 

 

 

1

10

10

1

10

10

(5.18)

q2 q20 x20 x10 ,

q2

q20 v20 v10.

 

94

Интегрируя дифференциальные уравнения движения (5.17) с учетом начальных условий (5.18), можно получить искомые зависимости q1(t) и q2(t). Однако в данной задаче целесообразно поступить иначе: сначала определить движение центра масс системы, а затем изучить колебания тел А и В в подвижных осях, связанных с центром масс, т. е. разбить абсолютное движение на переносное и относительное.

Положение центра масс системы Ц задается его радиус-вектором:

 

 

 

 

 

m r

m r

 

 

 

 

r

1 1

 

2

2

.

 

 

 

 

m1 m2

 

 

 

 

Ц

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С учетом (5.12) получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

m2

 

 

 

 

r

q

 

 

 

 

q

e .

(5.19)

 

 

 

Ц

 

1

 

m1 m2

 

2

1

 

 

 

 

 

 

 

 

 

 

 

Осуществим переход от координат q1, q2 к координатам qЦ, q2. Пола-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

гая

rЦ

qЦe1 , имеем выражение для обобщенной координаты qЦ

центра

масс системы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qЦ q1

m2

 

 

 

q2 .

 

 

 

 

 

 

 

 

 

(5.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m1 m2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Видим, что с учетом (5.20) первое из уравнений (5.17) дает диффе-

ренциальное уравнение движения центра масс связки двух тел:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F2

.

 

 

 

 

 

 

 

(5.21)

 

 

 

 

 

 

(m1 m2 )qЦ F1

 

 

 

 

 

 

 

 

 

Начальные условия движения имеют вид:

 

 

 

 

 

 

 

 

 

при t = t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qЦ qЦ0 q10

 

 

m2

q20 ;

 

 

 

 

 

 

 

 

 

 

m1 m2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m2

 

 

 

 

 

(5.22)

 

 

 

 

 

 

 

qЦ

qЦ0

q10

 

 

 

 

q20 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m1 m2

 

 

 

 

 

 

 

Перейдем к изучению относительного движения тел A и B. Выражая

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q1 из первого уравнения (5.17) и подставляя во второе, получим:

 

 

 

 

 

 

 

 

 

k

2

(q2

l0 ) f2 (t) ,

 

 

 

 

(5.23)

 

 

 

 

 

 

q2

2bq2

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 2 c

m1 m2

,

2b

m1 m2

.

f

(t)

F1 (t)

 

F2 (t)

 

(5.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m1m2

 

 

 

m1m2

 

 

 

 

2

 

 

 

 

m1

m2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим деформацию пружины в текущий момент как = l = q2 l0, тогда уравнение (5.23) примет вид:

 

 

2

(t) ,

(5.25)

 

2b k f2

95

оно представляет собой дифференциальное уравнение затухающих колебаний тел A и B относительно центра масс системы.

Начальные условия для (5.25) имеют вид:

при t = t0, = 0 = q20 l0,

 

 

q20

.

(5.26)

0

Найдем положение равновесия данной механической системы (см.

рис. 5.3), и исследуем его устойчивость.

 

 

 

 

 

Потенциальная энергия (5.14) имеет вид:

 

 

 

 

 

1

c( l)2

1

c(q l )2

,

 

 

2

 

 

 

2

 

2

0

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

Π

0 ,

Π

c(q l ).

 

 

 

 

 

 

q1

q2

 

 

2

0

 

 

 

 

 

 

 

Делаем вывод о том, что положение равновесия является безразличным по отношению к q1. Что касается второй координаты, то равновесие

будет иметь место при q2 = l0, и оно будет устойчивым, т. к. c 0 .

q2

5.4 Исследование колебаний подпрыгивания и галопирования кузова на рессорах

На рис. 5.4 показана схема двухосной тележки. Кузов представляет собой твердое тело массы m, момент инерции которого относительно оси, проходящей через центр масс Ц перпендикулярно плоскости чертежа, равен JЦ. Массами колес и подвески пренебрегаем.

Подвеска каждой оси состоит из пружины жесткости с и гидравлического демпфера с коэффициентом вязкого сопротивления . Точки L1 и L2 крепления подвески к кузову расположены симметрично относительно центра масс на расстоянии b от него (т. е. база тележки равна 2b).

z

 

 

 

 

 

 

 

 

 

 

 

Fупр.1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fупр.2

v0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R2

 

 

r

 

 

 

 

 

 

 

 

L1

 

 

 

 

 

 

 

 

 

 

 

 

 

Ц

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

L2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R1

 

 

 

 

 

 

 

 

 

 

 

 

O

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

e1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

96

Предполагая, что центр масс движется по горизонтали с заданной постоянной скоростью V0, требуется составить дифференциальные уравнения колебаний подпрыгивания и галопирования кузова при движении по пути с неровностью h = h(x).

На рис. 5.4 горизонтальной пунктирной линией отмечен уровень, соответствующий положению статического равновесия центра масс кузова при h = 0. Ось Ох направлена вправо вдоль этой линии, ось Оz – вертикально вверх. Кузов совершает плоскопараллельное движение, его положение задается горизонтальным перемещением центра масс хЦ = V0t, вертикальным отклонением центра масс от его уровня статического равновесия z (подпрыгивание) и углом (галопирование).

При изучении колебаний подпрыгивания и галопирования примем в качестве обобщенных координат q1 = z и q2 = . Число степеней свободы равно двум.

Кинетическая энергия кузова согласно (2.9) запишется в виде:

 

1

mz

2

 

1

 

 

2

 

 

T

 

 

 

 

J Ц

 

 

.

(5.27)

2

 

2

 

 

 

 

 

 

 

 

 

 

 

Потенциальную энергию силы тяжести и сил упругости пружин

находим согласно (2.19):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

Π mgz

c ( l j )2 .

(5.28)

 

 

2

 

 

 

 

 

 

 

 

 

j 1

 

Длина пружины задней подвески равна

 

l1 l0

ст z bsin h1 ,

 

где ст = mg/2c – статическая осадка пружины;

h1 = h(x1) – вертикальное перемещение оси заднего колеса K1 за счет наезда на неровность; координата х1 = хЦ b = v0t b.

Предполагая угол малым, т.е. полагая sin , cos 1, получаем деформацию пружины задней подвески в виде:

l1 l1 l0 ст z b h(x1) .

Аналогичным образом для пружины передней подвески находим:

l2 l2 l0 ст z b h(x2 ) ,

где координата x2 = хЦ + b = v0t + b.

Выражения для l1 и l2 подставляем в (5.28) и находим потенциальную энергию П, затем вычисляем согласно (2.15) обобщенные силы:

Q1

/ z 2cz c(h(x1 ) h(x2 )),

 

(5.29)

Q

/ 2cb2 cb(h(x ) h(x )).

 

2

1

2

 

97

Диссипативная функция (2.21) в нашем случае принимает вид:

 

1

2

 

 

vj 2 ,

(5.30)

 

 

 

 

2 j 1

 

где vj – разность скоростей точек Lj и Kj ( j =1, 2) крепления демпферов к кузову и к осям.

Нахождение vj проводим аналогично тому, как это было проделано в п. 5.1, получаем:

 

dh

 

 

 

 

v1 z b v0

 

,

v2 z b v0

dh

x x2 .

dx

 

 

 

 

 

x x1

dx

 

 

 

 

 

Находим обобщенные силы:

Q

/ z 2 z v

dh

 

 

 

dh

 

 

 

 

 

 

 

 

 

1

0

dx

 

x x1

dx

 

x x2

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

/ 2 b2 v b dh

 

dh

 

 

 

2

 

0

dx

 

x x1

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Записываем уравнения Лагранжа согласно (2.22):

,

.

xx2

 

2

 

2c

 

c

 

v

 

 

dh

 

 

 

 

 

 

 

 

 

 

 

 

z

z

z

h(x1) h(x2 )

dh

 

 

 

 

 

 

 

 

0

 

 

 

 

,

(5.31)

m

m

m

m

 

 

 

 

 

 

dx

x x

dx

x x

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

2 b2

 

2cb2

 

cb

h(x1 ) h(x2 )

v b dh

 

 

 

dh

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

. (5.32)

 

 

 

 

 

 

 

J Ц

 

J Ц

 

J Ц

 

J Ц

dx

 

x x1

 

dx

 

x x2

 

 

 

 

 

 

 

 

 

Как следует из полученных уравнений, колебания подпрыгивания и галопирования кузова разделены, т.к. в уравнение (5.31) не входит угол , а в уравнение (5.32) – координата z.

Выражение для неровности пути принимаем в виде:

 

 

 

2 x

 

h(x ) h0

 

cos

 

 

 

 

 

 

 

1

L p

 

.

 

Начальные условия движения:

при t

=

t0, z(t0 ) z0 ,

z(t0 ) z0 ,

(t0 ) 0 , (t0 ) 0 .

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

Потенциальная энергия дается выражением (5.28):

Π Π(0)+cz2 cb2 2 .

98

Вычисляя производные и применяя критерий Сильвестра (3.12), получаем, что положение равновесия z = 0, φ = 0 будет устойчивым.

5.5 Исследование колебаний подпрыгивания пути и железнодорожного экипажа с одноступенчатым рессорным подвешиванием

Составим уравнения колебаний для модели рельсового экипажа в виде механической системы, включающей в себя две массы: масса m2, моделирующая кузов, и масса m1, моделирующая тележку. Эти массы связаны между собой через кузовную подвеску, состоящую из пружины жесткостью с2 и гасителя колебаний с коэффициентом демпфирования 2

(рис. 5.5). Экипаж движется по деформируемому пути с неровностью t .

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

зонтали система перемещается с постоянной скоростью.

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

Из сравнения расчетных схем на рис. 5.5 и 5.6 легко видеть, что c1 8c , c2 4сz , m2 mк , m1 2m .

Рассматриваемая система имеет две степени свободы. За обобщенные координаты примем q1 zт и q2 zк , отсчитываемые от равновесного положения и определяющие положение центров масс тележки и кузо-

99

ва. При условии безотрывного качения колеса по рельсу перемещение точки пути под колесом zп zт .

Рис. 5.6. Расчетная схема четырехосного вагона с одинарным рессорным подвешиванием

Уравнения Лагранжа второго рода для рассматриваемой системы запишем в виде:

 

d

 

T

 

 

 

T

 

 

 

 

 

 

 

 

z

 

dt

z

 

 

 

 

 

к

 

 

 

 

к

 

 

 

 

 

 

T

 

 

 

 

T

 

d

 

 

 

 

 

 

 

 

 

 

z

 

dt

z

 

 

 

 

 

т

 

 

 

 

т

 

 

 

П Ф ;zк zк

(5.33)

П Ф .zт zт

Кинетическая энергия системы может быть представлена в виде:

2

2

 

2

T

m1zт

 

m2 zк

 

mп zт

.

 

 

 

2

2

2

 

Потенциальная энергия сил тяжести запишется как

Птяж m2 gzк m1gzт mп gzп.

Потенциальная энергия сил упругости имеет вид:

Пупр с22 l2 2 с21 l1 2.

Деформация пружины, моделирующей упругость пути

ся как

l1 1, ст zт ,

(5.34)

(5.35)

(5.36)

l1 , находит-

100