CC BY
52
УДК 625.1:51.001.57
В. В. ЛАГУТА (ДИИТ)
НЕПРЕРЫВНЫЙ КОНТРОЛЬ СДВИЖЕК В ЗАДАЧЕ АВТОМАТИЗАЦИИ РАСЧЕТА ВЫПРАВКИ ЖЕЛЕЗНОДОРОЖНЫХ КРИВЫХ
Наведено чисельний метод розв'язання диференщального pibHHHHH здвигiв, який дозволяе обчислити ве-личини здвигiв в будь-якш точцi залiзничноl криво! при проведенш и виправки. Здвиги, обчисленнi мiж точками вимiру стрiл вигибу, можуть бути використанi як контрольна
Приводится численный метод решения дифференциального уравнения сдвижек, позволяющий вычислять величины сдвижек в любой точке железнодорожной кривой при проведении ее выправки. Сдвижки, вычисленные между точками измерения стрел изгиба, могут быть использованы как контрольные.
The numerical method of solution of a differential equation of railway shifts is presented allowing to calculate the quantity of shifts in any point of a railway curve during its straightening. Shifts calculated between the points of measurement of bend arrows can be used as control ones.
Ритмичная, безаварийная работа железнодорожного транспорта зависит от технического состояния пути, его положения в плане. Решение задачи проектирования выправки железнодорожного пути имеет важное практическое значение для текущего содержания, в связи с этим проводятся исследования и предлагаются различные методики расчета.
Весомый вклад в развитие методов расчета выправки кривых принадлежит В. И. Ангелей-ко, В. Б. Бредюку, А. В. Гавриленкову, И. В. Гоникбергу, Б. Н. Евдаеву, Д. Г. Ковтуну, П. Г. Козийчуку, И. П. Корженевичу, А. А. Лебедеву, А. И. Проценко, Е. М. Рапопорту, Н. Г. Ренгачу, И. Я. Туровскому, А. П. Шутову из зарубежных ученых и инженеров можно назвать Dobrovolski B., Degorski T., Ellis D., Hull H. A., Kellie A. C., Moriyama T., Rodwan A. E., Rajaram B., Rao M. S., Sitko S. и др. Среди последних изысканий в решении этой проблемы следует отметить работы
А. Я. Когана и С. В. Петуховского [1, 2], Л. Маркса [3], Г. Г. Орлова [4].
Среди методов расчета сдвижек для проектных кривых можно указать такие: методы, использующие сглаживающие функции; усредняющие методы; методы, использующие аппарат приближения функций; методы, использующие аппарат теории вероятностей; методы, использующие аппарат математического программирования; смешанные методы.
Натурная кривая рассматривается для некоторого интервала железнодорожного пути 5 е [0,1], 5 — длина кривой. Предполагается, что проектная кривая получается из натурной кривой с помощью сдвижек и (5) по нормалям натурной.
В [5] приводится вывод точного дифференциального уравнения сдвижек на основе аппарата теории кривых
U "(S) = Kg (5)
Kp
1 + и(s)Kg (s) -
(1 + и(s)Kg(s))2 -u'2(s) 2 -u'(s)(2и'Kg(s) + и(s)Kg'(s))
1 + и (s)Kg (s)
(1)
Из точного дифференциального уравнения циальное уравнение сдвижек, которое имеет сдвижек получено приближенное дифферен- вид
и "(и, и') = к8 - Кр +(К2 - 2К8Кр ) и + 0 • и' +
-2 К\ • Кр • и2
+ (4К8 -3Кр)и'2 + 2К'&и • и'1 +о(и2,и'1 ), (2)
где К (5) - кривизна натурной кривой, Кр (5) - кривизна проектной кривой. Из уравнения (2) получается еще одно приближенное дифференциальное уравнение сдвижек
и " -(К2 - 2 К К
)и-К -КР) = 0. (3)
которое предлагается использовать в случаях, когда традиционное уравнение сдвижек
и " = К - Кр
(4)
не обеспечивает необходимую точность вычислений. Уравнения (1), (2) являются нелинейными, в этом состоит сложность решения поставленной задачи даже численными методами. Уравнения (3), (4) - линейные, для таких типов уравнений существуют весьма эффективные методы.
При проведении выправки железнодорожной кривой в основном руководствуются расчетными величинами сдвижек в местах измерения стрел изгиба. Величины сдвижек между точками измерений обычно не рассчитываются и во время проведения выправки не контролируются. Надо ли проводить этот контроль или не надо? Автор этот вопрос не оспаривает. На этот счет никаких рекомендаций в инструкциях по содержанию пути не говорится.
Последние факты введения скоростного пассажирского движения (Киев - Днепропетровск, Киев - Харьков, Днепропетровск - Донецк) говорят о том, что необходимость в качественном содержании пути актуально существует. Как известно, при высоких скоростях движения путь расстраивается интенсивней, что с необходимостью приводит к повышенному вниманию относительно его содержания. В статье не рассматривается какой-либо вычислительный метод по определению рациональных величин сдвижек. Рассматривается только численный метод, позволяющий решать нелинейное дифференциальное уравнение сдвижек как первую краевую задачу. Решение данной задачи позволяет определить величины сдвижек в любой точке железнодорожной кривой, если известны ее натурная и проектная кривизна. Полностью численный метод описан в статьях [6, 7]. Решение основывается на аппроксимации линейным сплайном второй производной сдвижки на интервале железнодорожной кривой.
Пусть известны кривизна натурной кривой и кривизна проектной кривой. Правую часть уравнения (1), равно как уравнений (2), (3), (4), обозначим через /(и, и', 5), где 5 - длина кривой от начальной точки измерения (выправки), и - величина сдвижки в точке 5 , и' - производная от величины сдвижки в точке 5 .
Для определения сдвижек рассматривается первая краевая задача для обыкновенного дифференциального уравнения 2-го порядка с однородными граничными условиями
и" = / (5, и, и'), и (а) = 0, и (Ь) = 0. (5)
Приближенное решение строится в классе функций С3. Предполагается, что выполнены условия существования решения нелинейной краевой задачи [8], которое может быть получено последовательными приближениями
иР = f(5,ир-ииР-1) , Р = 1,2••• .
(6)
Каждое решение ир должно удовлетворять
краевым условиям (5). Предел последовательности решений ир будет приближенным решением исходной задачи.
Отрезок интегрирования [а, Ь] приводится
заменой переменой
5 - а
t =- к единичному
Ь - а
интервалу 0 < t < 1. Функция и ^) также должна удовлетворять однородным граничным условиям
и (0) = 0, и (1) = 0.
(7)
Вторая производная искомой функции и" = f ^, и, и') аппроксимируется кусочно-линейной функцией ^, звенья которой соединяются в узлах сетки {ti} на отрезке [0, 1]. Приближенное решение в этом случае получается в виде кубического сплайна й^). Аналогичный подход реализован для решения нелинейных интегральных уравнений Вольтерра и описан в
[9, 10].
На отрезке интегрирования вводится регулярная сетка {^ }
ti = i • И, i = 0, п . Вторая производная
и" = f (^ ир), и'р))
заменяется кусочно-линейной функцией с вер- где индекс р обозначает номер итерации (при-шинами в узлах сетки ближения) и /¿р-1 = / Ь, пр-1, и' р-1 ),
РгР (*) = ' /г-1 (* - и-1)
1 , /гР-1 - / --1 и , ) ир-1 = иР-1((,), и'р-1 = и'Р-1((,), р = 1,2,.... Для
отрезка [*;-1, ] получается следующее приближенное решение для узлов сетки [6]
*г - *г-1
Гг-1 < * < 1г, г = 1, п ,
— г
(г - п) -Г Ё[(3з - 2) /А1 + (3 з -!) /р-1 ] -
6 з =Г
—3 п _
-г— £[(3п - 3 з + 2)/зр-1 + (3п - 3 з + 1)/зр-1 ], г = 1, п -1, (8)
6 1=г +1
г
р
и- = -- Ё[(31 - 2)/Зр-11 + (3 з -1) /р-1 ]■
6 з=Г
— 2 п _
-— Ё[(3п-3з + 2))-1 +(3п-3з +1)//-1 ], г = 0,п , (9)
6 з =г+Г
где У0 = / (*з, и0, и'0), и0 = 0, и'0 = 0, з = М, р = 0,1,2,... .
Приближенное решение и(*) для любой ляется следующим образом: точки отрезка интегрирования [0, 1] опреде-
й(з) = (* -1) Ё [(3 з - 2) / -1 + (3 з -1) / ] + (* -1) ф ^+*{£ -1) у ©ё ^ -
6 з=1 I, *
(¿>1) г
—3 п _
— Ё [(3п - 3з + 2)/-1 + (3п - 3з +1)/ ], *г < * < 1г+1, г = 0, п -1. (10)
6 У=г+2
(г< п-1)
Точность данного приближенного решения 4
проанализирована в [7]. Погрешность для лю- |и(*) -„(*)| < max|/"(t)| • —*(1 - *).
бой точки * решения (8, 10) составляет 0<*< 24
2 Данная ошибка достигает наибольшего значе-
Ли (*) = *(1 - *)-• 7шк|/''(*)|,
ния и оценивается как
'„2
учитывая линейное представление функции в |и(0.5) - ¿1(0.5) < max|/"(*)|__.
л -Ж-Л -Ж- -Ж- I I 1 I '
окрестности точки в форме Коши. Наиболь-
0<*<11 ' 96
шая погрешность определяется На реальном примере покажем, как будут
7 отличаться сдвижки, определённые из диффе-
|и(0.5) -и(0.5)| < max|/"(0|---2. ренциальных уравнений (1), (3), (4). Считаем,
1 1 0<* <11 1 24 Т.Г
что заданы натурная кривизна К и проектная
В граничных точках погрешность равна нулю кривизна Кр . Натурная и проектная кривые вследствие выполнения граничных условий для
точного и приближенного решений. При оценке заданы таблицей значений стрел. Исходные
разности точного и приближенного решения с данные и результаты приближённого численно-
помощью остаточного члена интерполяционно- го решения дифференциальных уравнений (1),
го многочлена Лагранжа, получим (3), (4) представлены в табл1
и
Таблица 1
Результаты вычислений сдвижек соответствующие различным дифференциальным уравнениям
Номер точки Стрелы Сдвижки согласно дифференциального уравнения, (мм)
натурной кривой (мм) проектной кривой (мм) классического [уравнение (4)] линейного [уравнение (3)] Точного [уравнение (1)]
0 0.00 0.00 0.00 0.00 0.00
1 5.00 19.74 -4.91 -4.98 -4.98
2 39.00 39.47 -29.63 -29.77 -29.78
3 84.00 59.21 -51.63 -51.86 -51.88
4 54.00 78.95 -49.06 -49.39 -49.41
5 130.00 98.69 -61.04 -61.50 -61.53
6 80.00 99.10 -46.23 -46.79 -46.82
7 120.00 99.10 -39.60 -40.26 -40.29
8 95.00 91.09 -10.29 -11.00 -11.04
9 58.00 82.89 22.84 22.11 22.08
10 81.00 74.70 26.19 25.48 25.44
11 72.00 66.50 31.47 30.80 30.77
12 32.00 66.16 35.00 34.39 34.36
13 90.00 66.16 2.88 2.32 2.29
14 83.00 66.16 -3.23 -3.74 -3.77
15 32.00 66.16 9.67 9.22 9.18
16 82.00 66.16 -12.08 -12.48 -12.51
17 61.00 66.16 -25.83 -26.16 -26.20
18 80.00 66.16 -36.56 -36.83 -36.86
19 69.00 66.16 -29.62 -29.81 -29.84
20 63.00 66.16 -15.34 -15.43 -15.47
21 79.00 66.16 -0.05 -0.05 -0.08
22 40.00 66.16 22.59 22.69 22.66
23 71.00 62.07 17.13 17.33 17.30
24 59.00 57.01 15.22 15.51 15.48
25 48.00 51.96 17.62 18.00 17.97
26 44.00 46.90 14.44 14.87 14.86
27 40.00 45.90 4.58 5.09 5.07
28 55.00 45.90 -10.76 -10.18 -10.21
29 51.00 45.90 -14.24 -13.59 -13.61
30 40.00 45.90 -9.85 -9.13 -9.15
31 47.00 45.90 -11.27 -10.47 -10.49
32 49.00 45.90 -12.16 -11.27 -11.29
33 46.00 45.90 -8.52 -7.55 -7.56
34 42.00 45.90 -5.02 -3.95 -3.96
35 57.00 45.90 -2.99 -1.82 -1.83
36 44.00 45.90 11.91 13.18 13.17
37 34.00 45.25 24.00 25.37 25.36
38 37.00 36.07 20.38 21.85 21.84
39 30.00 26.89 15.27 16.83 16.82
40 12.00 17.71 12.71 14.35 14.35
41 9.00 8.53 3.72 5.44 5.43
42 7.00 0.00 -4.00 -2.22 -2.22
43 0.00 0.00 0.00 0.00 0.00
Приближённое численное решение находилось предложенным методом [6], граничные условия выбирались однородными. Перед подстановкой в соответствующее дифференциальное уравнение кривизна натурной и проектной
кривой определяются по значениям стрел. Между точками измерений проводится линейная интерполяция. При приближенном решении дифференциального уравнения сдвижек необходимо решить еще и задачу аппроксимации
натурной кривизны (5) и проектной кривизны Кр (5). В том случае, если кривизна проектной кривой выбрана геометрически правильной структуры, то для этой кривой задача аппроксимации отпадает (она уже предопределена линейной зависимостью). Аппроксимация как натурной кривизны, так и проектной зависимостью, отличной от линейной, например, приближение К? (5) и Кр (5) параболическими
сплайнами, приводит к появлению эффекта Гиббса в виде характерных колебательных процессов, отсутствующих у исходных функций. Наиболее это свойственно для участков железнодорожной кривой с очень медленно меняющейся кривизной, а также в местах сопряжения участков кривых различной структуры с явно присутствующими особенностями. Эффект Гиббса может привести к существенным ошибкам при определении величин сдвижек, особенно при аппроксимации прямых и круговых участков пути.
ВЫВОДЫ
Для контроля качества проведения выправки железнодорожной кривой можно использовать значения сдвижек в точках лежащих между точками измерения стрел изгиба. Значение сдвижек могут быть получены как решение первой краевой задачи для обыкновенного дифференциального уравнения второго порядка (1), (3), (4), определяемое в соответствии с выражением (10).
В зависимости от вычислительных требований при определении сдвижек, можно пользоваться соответствующим дифференциальным уравнением.
При аппроксимации кривизны натурной и проектной кривых следует избегать эффекта Гиббса.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Коган А. Я. К вопросу о расчете положения пути в плане по стрелам несимметричной измерительной хорды / А. Я. Коган, С. В. Петухов-ский // Вестник ВНИИЖТ. - М. - 2002, № 1 -С. 26-29.
2. Коган А. Я. Расчет оптимальных сдвигов при выправке криволинейных участков пути / А. Я. Коган, С. В. Петуховский. // Вестник ВНИИЖТ. - М. - 2002, № 5.
3. Маркс Л. Измерение параметров пути с привязкой к системе глобального позиционирования // Железные дороги мира. - 2002, № 7. - С. 71-75.
4. Орлов Г. Г. Обоснование требований к точности специальной реперной сети для контроля профиля пути в плане и в профиле / Труды ЦНИИС. - М. -2001. - С. 62-75.
5. Лагута В. В. Совершенствование проектирования кривых железнодорожного пути в плане: Диссертация на соискание ученой степени к.т.н.: 05.22.06. - Д.: ДИИТ. - 2002. - 186 с.
6. Лагута В. В. Метод итерации решения первой краевой задачи с нелинейной правой частью / Д.: Вюник ДНУЗТ. - Вип. 18, 2007.
7. Лагута В. В. Погрешность приближенного решения первой краевой задачи с аппроксимацией второй производной линейным сплайном / Д.: Вгсник ДНУЗТ. - Вип. 19, 2007.
8. Камке Э. Справочник по обыкновенным дифференциальным уравнениям: Пер. с нем. - М.: Наука. Д.: Вюник ДНУЗТ. - Вип. 18, 2007. 1976. - С.61-65.
9. Верлань А. Ф. Интегральные уравнения: Методы, алгоритмы, программы. Справочное пособие / А. Ф. Верлань, В. С. Сизиков - К.: Наук. думка - 1986. - С. 71-79, 84-94, 142-143, 148152.
10. Мейнарович Е. В. О применении интерполяционных сплайнов к решению нелинейных интегральных уравнений Вольтерра / Е. В. Мейна-рович, Р. В. Поляков, Л. Н. Шлепаков - В кн.: Линейные и нелинейные краевые задачи математической физики. - К.: Ин-т кибернетики АН УССР, 1974 - С.204-212.
Надшшла до редколеги 03.01.2008.
