CC BY
32
Л.В. Кудюров, Ю.К. Мустафаев
МОДЕЛИРОВАНИЕ УГОНА РЕЛЬСОВОГО ПУТИ ПОД ДЕЙСТВИЕМ ДИНАМИЧЕСКИХ НАГРУЗОК
Рассматривается динамическая модель угона железнодорожного пути в процессе его эксплуатации. Рельсовая нить представлена как механическая система из п элементарных участков, связанных упругими и диссипативными связями друг с другом и с основанием, на котором располагается рельсовая плеть. Применение принципа Даламбера к выделенному элементу рельсового пути дало возможность свести задачу к системе алгебраических уравнений, позволяющих оценить угон пути в зависимости от нагрузки.
Явление угона пути, т. е. продольного перемещения рельсовой плети известно давно. Причинами угона является косой удар, испытываемый рельсом при переходе колеса через стык, температурные деформации, прогибы рельсов под действием динамических сил [1]. Последнее обстоятельство вынуждает колесо все время идти на подъем, заставляя изогнутый рельс, выпрямляясь под действием веса колеса, «выстилаться» впрямую, вызывая тем самым перемещение рельса и упругие деформации рельсовых креплений и подрельсового основания в направлении движения (рис. 1).
Р и с. 1. Процесс «выстилания» рельса перед набегающим колесом
Из всех указанных причин особый интерес представляет последняя, так как она обусловлена только динамическими факторами. Особую актуальность эта гипотеза приобретает при рассмотрении бесстыкового пути.
Существующие методы расчета динамического угона пути достаточно сложны и требуют решения систем дифференциальных уравнений в частных производных. При этом возникают сложности вычислительного характера, связанные с необходимостью учета большого числа движущихся сил, равных количеству колесных пар состава. Поэтому представляет интерес разработка эффективного алгоритма вычислений, отличающегося быстродействием и достоверностью расчетов.
Для реализации поставленной задачи была построена модель, использующая следующие допущения:
рельс рассматривается как система элементов, связанных друг с другом и с основанием упругими связями (рис. 2);
каждый элемент представляет собой материальную точку, объединяющую в себе шпалу и участок рельса над ней, длина которого равна длине междушпального расстояния; масса эле-
мента равна суммарной массе шпалы, рельсовых скреплений и заменяемого участка рельса (рис. 3);
возмущения не передаются мгновенно от элемента к элементу; существует определённая скорость распространения возмущения (для металлов — скорость звука, равная примерно 5000 м/с).
Р и с.2 Модель рельсовой плети, состоящая из отдельных элементов:
1 - элементы массой т, равной массе моделируемого участка; 2 - упругие связи;
3 - демпферы
Р и с. 3. Разбиение рельсового пути на элементарные участки
Очевидно, что данная система склонна к распространению колебательных (волновых) процессов. Наличие конечной скорости распространения возмущающего воздействия позволяет рассматривать выделенный участок рельса на коротком промежутке времени как отдельную изолированную колебательную систему. Интервалы времени должны быть таковыми, чтобы изменение характеристик системы происходило с достаточной плавностью (не более 3 - 5 % от любого параметра). Исходя из заданной точности, нетрудно подсчитать, что при чистом синусоидальном колебательном процессе необходимо делать 200 выборок на период (примерно 3% точность).
Учитывая, что длина междушпального промежутка равна примерно 0,5 м, а скорость распространения волны 5000 м/с, а также то, что нет физического смысла в волне, длина которой меньше междушпального расстояния (нет опоры для колебаний), делаем вывод, что необходимые характеристики должны быть следующими: интервал времени (период волны) -
Т = 500005 () >= 0,0001 (с), (1)
5000 (м/ с) шаг по времени (период волны) -
Т 7
Ы =----= 5 • 10-7 (с). (2)
200
При скорости движения возмущающей силы, равной 20 м/с, её время движения по рельсовой плети длиной 15 м составит:
1Р 15
/дв = — = — = 0,75 с, (3)
дв и 20
где 1р — длина рельса, и- скорость движения возмущающего воздействия.
Пусть необходимо дополнительно рассмотреть свободные колебания на промежутке времени, равном 3/дв. Тогда общее время процесса равно:
Т = 4-/дв = 4-0,75 = 3 с. (4)
Учитывая, что шаг времени Д/ равен 5 • 10-7 с, приходим к выводу, что необходимо сделать следующее количество выборок:
п =— = 6 -106. (5)
Ы 5-10-7
Кроме того, на каждом шаге необходимо просчитать параметры каждого элемента, число которых для стандартной рельсовой плети равно 27, следовательно, необходимо сделать
27п = 27^106 = 1,62 108 (6)
количество одинаковых расчетов, что, очевидно, под силу лишь достаточно производительной ЭВМ.
Упростить задачу можно, применив наиболее эффективный формульный ряд с точки зрения арифметики ЭВМ. Это значит, что уравнения должны по возможности состоять из простых арифметических действий, не используя сложных математических функций, для которых ЭВМ использует отдельные встроенные алгоритмы.
Таким условиям, на наш взгляд, удовлетворяет система уравнений равновесия Д’Аламбера. Запишем ее применительно к нашей модели. Выделим элемент модели вместе с системой действующих на него сил (рис. 4).
"►х Ухф,!)
г
Р и с. 4. Силовая схема
На рис. 4 введены обозначения: [Р] — внешняя вертикальная сила (вес состава); [Т] — горизонтальная внешняя сила, вызванная сопротивлением колес движению; Q(i) — поперечная (срезывающая) сила, возникающая из-за разницы в 7 - координатах (/ - 1) и i -того элементов; Q(i+l) — аналогично, для i и ( i + 1) - го элементов; Л>1) — продольная сила вызванная деформацией рельса между (/ + 1) и i -тым элементами; Л(,+1) — то-же, для ьтых и (/ + 1)-тых элементов;
и ^2С(,) — диссипативные силы, всегда направление противоположно вектору скорости;
Рхпо({) и Гтпоа) — силы упругости подрельсового основания; Jx(i) и — силы инерции.
Для каждого элемента может быть составлена система двух уравнений с двумя неизвестными ах и а7:
-тах(у) + [Т](/) — Рхе(^) — -Ршп(у) + Л(1-1з) — Л(1+1£ = 0;
-та2(1^) + [р](^) — Ргс(ц-) — Ршп.^) + Q(i-1з + Q(i+lj) = 0, (7)
откуда получим ускорения в виде
а _ [Т](Г) - -^хфу) - ^Хпп(у) + Л\i-1J) - Ла+1,Л . х(у) _ >
т
а _ [Р](}) - ^’7с(},3) - ) + ^-1,Л + ^М^") (о)
) _----------------------------------------- . (0)
т
Исходя из поставленных ранее условий, свидетельствующих тому, что сила (возмущение) не могут передаваться от элемента к элементу мгновенно, можем заключить, что на промежутке времени А/ система сил будет кинетостатической.
В этом случае найденные ускорения аХ1 и абудут также постоянны и будут отличаться от предыдущих не более чем на 0,03, что обеспечивает достаточную точность аппроксимации. Используя формулу для равноускоренного движения
и = и 0 + аД/, (9)
определим скорость в конце интервала времени Д/, исходя из следующих граничных условий:
конечная скорость в конце предыдущего интервала будет начальной в начале данного интерва-
ла, т.е.
°х(у) _ их^3) + ах(у) Д/ ;
иг(у) _ иг0(у) + ^(у) Д/ . (10)
Аналогичные граничные условия используем и для перемещений:
если а _ 0, то 8 _ 80 + и0Д/;
V V + и2-и2 (11)
иначе V _ V 0 +
0 2а
Применительно к нашей системе запишем:
2 2 Ux(i,j) Uxo(i,j) .
x(i, j) xo (i, j) _
2ax(i, j)
2 _ 2
о _ о , Uz(i,j) _ Vzo(i,j) /1
Sz(iJ) Szo(i,j) + 2 a ' '
2az(i,j)
Таким образом, на данном шаге были получены все интересующие нас параметры элемента модели, как-то: ускорения ах и az, скорости vx и vz, перемещения Sx и Sz для каждого элемента системы. При переходе к следующему интервалу времени необходимо передать начальные условия. Это делается следующим образом:
uxo(ij+1) = Uz(ij), uzo(y+1) _ Uz(y), Sxo(ij+1) _ Sx(ij), Szo(y+1) _ Sz(ij). (13)
Силы, входящие в систему уравнений (7), определяются из условий
Fxci = ku * Uxoi , Fzci = ku * Uzoi , Fxnoi = Sxoi * Жгю« , Fzпоi = Szoi * жnob ;
N,_1 = Жp« • ([xi _ x(i_1)]_\ [Sxo, _ Sxo(,_1) + x _ x(i_1)f + lSzo, _ Szo(,_1)f ] ; (14)
Ni+1 = жp« • ^[x(i+1) _ xi]_ \ [x i+1) _ xi + Sxo(i+1) _ Sxoi]2 + [Szo(i+1) _ Szoi] j ;
Q,_1 = _Жpb • {Szo, _ Szo(,_1)) , Q,+1 = _жpb • (Szo, _ Szo(,+1)) , Jxi = _OTax0i , Jzi = _maz0i ,
где - коэффициент пропорциональности, больший нуля; жр, жпо - коэффициенты жесткости рельса и подрельсового основания соответственно.
Если координата хр силы Р находиться между Xi и х^), то величина воздействия от силы Р на i - тый элемент определяется как
[р]_ Хр Х-1 р . (15)
X - хы
Если координата хр силы Р находиться между Xi и Х(+1), то величина воздействия от силы Р на i - тый элемент определяется как
[р ]_х+1 - Хр р. (16)
xi+1- х
Если сила Р находится в окрестности i - того элемента (хр ~ xi ), то [Р] = Р.
Сила продольного сопротивления движению равна:
Т _ , (17)
где ю0 = 2,5 Н/кН — основное удельное сопротивление движению.
На основе приведённого алгоритма была составлена программа для ЭВМ и произведены расчёты, анализ которых подтверждает справедливость предложенного метода моделирования. Полученные результаты во многом схожи с результатами аналогичных работ [2]. Наблюдались ярко выраженные колебательные процессы двух основных тонов: высокочастотные — с частотой собственных колебаний отдельных элементов, представляющих собой, в общем, физический пружинный маятник, и низкочастотные — колебания всей рельсовой плети в целом как
единой системы.
Данные расчетов приведены в таблице, а графики на основании этих данных — на рис. 5— 9.
Исходные данные
Вариант 1 Вариант 2 Вариант 3
Вертикальная жесткость подрельсового основания ж„оТ, Н/м 2,5-107 2,5-107 2,5-107
Г оризонтальная жёсткость подрельсового основания ’Жпо«? Н/М 2-107 1,5107 107
Максимальная сила удержания рельсового скрепления Рмах, Н 900 2000 2000
С целью упрощения вычислений (с сохранением наглядности работы алгоритма), все варианты расчетов приведены для десятишпального отрезка рельсовой плети, закрепленной в стыках.
На рис. 5-9 цифрами обозначено: 1 - график движения участка рельса над первой шпалой;
2 - график движения участка рельса над пятой шпалой; 3 - график движения участка рельса над десятой шпалой.
На графиках наглядно показано поведение различных участков рельса при прохождении по нему единичной силы, эквивалентной воздействию одной колёсной пары (для вагона с одноосной тележкой) или одной тележки вагона (с двухосной тележкой).
Анализ графиков указывает на то, что изменение горизонтальной жесткости подрельсового основания оказывает сильное влияние на характер протекающего колебательного процесса в горизонтальной плоскости и почти не оказывает влияния на характер вертикальных колебаний. С увеличением горизонтальной жесткости колебательный процесс приобретает более выраженные экстремумы, что говорит о больших внутренних силах, развиваемых в рельсе, и наоборот, снижение горизонтальной жесткости подрельсового основания делает график колебательного процесса более плавным.
Снижение прочности рельсового скрепления влечет за собой более раннее наступление колебательного процесса с проскальзыванием подошвы рельса по основанию. На графиках (рис. 8 и рис. 9), отражающих процессы при проскальзывании рельса отчетливо видны моменты наступления процесса угона. Данные моменты обуславливают появление «площадок» на графике колебательного процесса, что объясняется превышением внутренних сил в рельсе над силами удержания рельсового скрепления, когда происходит срыв подошвы рельса. При начавшемся
процессе скольжения подошвы рельса по подрельсовому основанию силы трения скольжения уравниваются с вызвавшей этот процесс возмущающей силой («площадка»). Кроме того, повышается частота колебаний, равная при обычном колебательном процессе 40^50 Гц, при проскальзывании повышающейся до 450^500 Гц.
Р и с. 5. Перемещения подошвы рельса по горизонтали для первого варианта расчета (одно деление по вертикали — 5-10"4 м, по горизонтали - 0,025 с ):
1 — перемещения для участка рельса над первой шпалой; 2 - перемещения для участка рельса над пятой шпалой; 3 — перемещения для участка рельса над десятой шпалой
Р и с. 6. Перемещения подошвы рельса по горизонтали для второго варианта расчета (одно деление по
вертикали — 5 • 10-4 м, по горизонтали - 0,025 с):
1 — перемещения для участка рельса над первой шпалой; 2 - перемещения для участка рельса над пятой шпалой; 3 — перемещения для участка рельса над десятой шпалой
Р и с. 7. Перемещения подошвы рельса по горизонтали для третьего варианта расчета (одно деление по вертикали — 5-10"4 м, по горизонтали — 0,025 с):
1 — перемещения для участка рельса над первой шпалой; 2 - перемещения для участка рельса над пятой шпалой; 3 — перемещения для участка рельса над десятой шпалой
Р и с. 8. Перемещения верхней части подрельсового основания по горизонтали при проскальзывании рельса (одно деление по вертикали - 10"4 м, по горизонтали - 0,025 с):
1 - перемещения для участка рельса над первой шпалой; 2 - перемещения для участка рельса над пятой шпалой; 3 - перемещения для участка рельса над десятой шпалой
Р и с. 9. Искусственно смоделированные перемещения верхней части подрельсового основания по горизонтали при проскальзывании рельса при малой жесткости подрельсового основания и малой силе удержания в рельсовых скреплениях (одно деление по вертикали - 10-4 м, по горизонтали - 0,025 с):
1 - перемещения для участка рельса над первой шпалой; 2 - перемещения для участка рельса над пятой шпалой;
3 - перемещения для участка рельса над десятой шпалой
Представленный выше упрощенный алгоритм расчёта и проведенные предварительные исследования причин угона и численной оценки его величины позволяют получить достоверные результаты и дать рекомендации по снижению силового воздействия, влияющего на угон пути.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК:
1. ШахуньянцГ. М., Железнодорожный путь. М.: Транспорт, 1987. 479 с.
2. Меньшикова В. И. Динамические продольные силы и перемещения рельсов железнодорожного пути (угон рельсов) // Динамические исследования пути и корректировка правил расчётов железнодорожного пути на прочность. Труды ЦНИИ. Вып. 466. М.: Транспорт, 1972. С. 83-187.
Поступила 22.07.2004 г.
