Иваночкин П.Г. Механика подвижного состава. Учеб пособ. 2019
.pdf
|
|
|
|
|
|
|
|
|
|
Проекция на ось О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, |
|
|
Q2П |
упр / 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
