ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
УДК 629.42.015
С. В. МЯМЛИН1, Н. Е. НАУМЕНКО2, О. М. МАРКОВА3, Е. Н. КОВТУН4,
В. В. МАЛЫЙ5*
*Каф. «Вагоны и вагонное хозяйство», Днепропетровский национальный университет железнодорожного транспорта имени академика В. Лазаряна, ул. Лазаряна, 2, Днепропетровск, Украина, 49010, тел. +38 (056) 776 84 98, эл. почта sergeymyamlin@gmail.com, ORCID 0000-0002-7383-9304
2Отдел «Динамика многомерных механических систем», Институт технической механики Национальной академии наук Украины и Государственного космического агентства Украины, ул. Лешко-Попеля, 15, Днепропетровск, Украина,
49005, тел. +38 (056) 375 68 01, эл. почта dep7@ukr.net, ORCID 0000-0002-6913-2568
3Отдел «Динамика многомерных механических систем», Институт технической механики Национальной академии наук Украины и Государственного космического агентства Украины, ул. Лешко-Попеля, 15, Днепропетровск, Украина,
49005, тел. +38 (056) 375 68 03, эл. почта dep7@ukr.net, ORCID 0000-0002-9390-4046
4Отдел «Динамика многомерных механических систем», Институт технической механики Национальной академии наук Украины и Государственного космического агентства Украины, ул. Лешко-Попеля, 15, Днепропетровск, Украина,
49005, тел. +38 (0562) 46 79 35, эл. почта dep7@ukr.net , ORCID 0000-0002-2961-3682
5 Отдел «Динамика многомерных механических систем», Институт технической механики Национальной академии наук Украины и Государственного космического агентства Украины, ул. Лешко-Попеля, 15, Днепропетровск, Украина,
49005, тел. +38 (0562) 46 79 35, эл. почта dep7@ukr.net, ORCID 0000-0002-6501-1398
ОЦЕНКА ДИНАМИЧЕСКИХ ПОКАЗАТЕЛЕЙ ШЕСТИОСНОГО ЛОКОМОТИВА
Цель. Работа посвящена оценке динамических качеств локомотива перспективной конструкции и определению возможности его использования на железных дорогах Украины. Методика. Для достижения поставленной цели использованы методы математического и компьютерного моделирования динамики железнодорожных транспортных средств, а также методы численного интегрирования систем обыкновенных нелинейных дифференциальных уравнений. Результаты. Для решения поставленной задачи построена расчетная схема локомотива на трехосных тележках, которая представляет собой систему твердых тел, соединенных элементами различной реологии. С использованием этой расчетной схемы разработана математическая модель пространственного движения локомотива, позволяющая исследовать его пространственные колебания при движении по прямолинейным и криволинейным участкам пути со случайными неровностями в плане и профиле. При составлении математической модели принимались во внимание как геометрические (нелинейность профиля поверхности катания колеса), так и физические нелинейности системы (работа сил сухого трения, нелинейность характеристик сил взаимодействия колес и рельсов). С помощью компьютерного моделирования проведены многовариантные расчеты, позволившие оценить динамические качества локомотива при его движении по прямолинейным и криволинейным участкам пути. При этом определялись показатели плавности хода локомотива в горизонтальной и вертикальной плоскостях, рамные силы, коэффициенты вертикальной динамики в первой и второй ступенях подвешивания, коэффициенты запаса устойчивости против схода колеса с рельса. Кроме того, проведено сравнение полученных результатов с аналогичными характеристиками широко используемого на железных дорогах Украины шестиосного тепловоза ТЭ116. Определено влияние скорости движения и технического состояния пути на безопасность движения локомотива. Научная новизна. Разработана математическая модель пространственного движения шестиосного локомотива перспективной конструкции, учитывающая особенности его конструкции, и сделана оценка динамических характеристик при движении экипажа по пути произвольного очертания. Практическая значимость. Разработанная математическая модель движения шестиосного локомотива может быть использована при проведении научных исследований по созданию перспективных конструкций локомотивов.
Ключевые слова: подвижной состав; локомотив; математическая модель; динамические показатели
Введение в полной мере отвечает требованиям, предъяв-
ляемым к современным транспортным систе-В настоящее время техническое состояние мам. Одной из самых острых проблем желез-
тягового подвижного состава Украины не ных дорог Украины является физический и мо-
doi: 10.15802/stp2015/46078
© С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
186
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
ральный износ локомотивного парка. Так, в настоящее время средний износ парка тепловозов составляет 97 %, электровозов - 90 %. В Программе обновления тягового подвижного состава железных дорог Украины на период до 2020 года заложены потребности в закупке подвижного состава нового поколения. Использование таких локомотивов позволит повысить безопасность движения поездов и значительно уменьшить потребление топлива и эксплуатационные расходы на их содержание [2, 4, 6, 10, 12].
Большую работу по созданию современных локомотивов проводят практически все мировые производители, в том числе и компания Siemens. В созданных ими тепловозах реализована значительная часть требований, предъявляемых к перспективным локомотивам будущего: возможность эксплуатации во всех странах Европы, использование как в грузовом, так и в пассажирском движении. Одним из таких перспективных локомотивов является шестиосный тепловоз с электрической передачей серии ER20, который может быть использован и на железных дорогах Украины после соответствующей адаптации [8, 14, 15].
Цель
Цель работы заключается в оценке динамических качеств локомотива перспективной конструкции при его движении по прямолинейным и криволинейным участкам железнодорожного пути и определение возможности его использования на железных дорогах Украины.
Методика
Для решения поставленной задачи была выбрана расчетная схема локомотива и составлена математическая модель, описывающая его пространственные колебания. На основе построенной математической модели разработан алгоритм и составлена компьютерная программа, предназначенная для оценки динамических характеристик движущегося локомотива.
Расчетная оценка динамических показателей локомотива проведена путем анализа результатов решения нелинейных дифференциальных уравнений, представляющих динамическую модель пространственных колебаний экипажа.
Результаты
Кузов тепловоза типа ER20 выполнен как самонесущая модульная конструкция и подразделяется на следующие компоненты: главную раму, боковые стенки, кабины управления и смонтированные лобовые части. Главная рама состоит из двух буферных брусьев, двух наружных и двух центральных продольных балок, четырех поперечных балок. Буферные брусья воспринимают ударно-тяговые усилия, которые через автосцепное устройство передаются от состава поезда на локомотив. Силы тяги и торможения передаются от тележек на кузов через штангу, присоединенную к буферному брусу. Кузов тепловоза опирается на две трехосные тележки. Конструкция тележки бесчелюстная с двухступенчатым рессорным подвешиванием с использованием винтовых пружин и гидравлических гасителей колебаний. Подвешивание тяговых электродвигателей опорно-рамное. Нагрузка от кузова передается на каждую тележку через упругое рессорное подвешивание второй ступени [1, 9].
Расчетная схема тепловоза типа ER20 (рис. 1) построена как механическая система, состоящая из девяти твердых тел (кузова, двух рам тележек и шести колесных пар).
Ш4 , . ШГГ
1ГГ —ПГ
$7 п ^4 І А ] [Ґ и У А I У АА ] [Iі и
И і ц
Н х -А
Рис. 1. Расчетная схема тепловоза типа ER20
Fig. 1. Design scheme of the locomotive, type ER20
doi: 10.15802/stp2015/46078
© С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
187
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
Выбранная конфигурация твердых тел расчетной схемы позволяет учесть следующие особенности тепловоза типа ER20:
- упруго-диссипативную с вязким трением связь кузова и рамы тележки;
- упруго-вязкую связь между рамой тележки и колесными парами: комплекты пружин буксового подвешивания допускают взаимные перемещения этих тел во всех направлениях;
- форму профилей поверхности катания колес.
Железнодорожный путь моделируется в виде приведенной к каждому колесу массой (двенадцать приведенных масс), имеющей только вертикальное и поперечное горизонтальное перемещения и опирающейся в этих направлениях на пружины и демпферы вязкого трения, моделирующие упруго-диссипативные свойства рельсов и подрельсового основания. Продольные и угловые перемещения элементов верхнего строения пути даной моделью не учитываются.
Конфигурация рассматриваемой механической системы описывается 9^6+12^2=78 координатами.
При описании движения тепловоза по участкам пути произвольного очертания выбрана неподвижная система координат , а для
каждого твердого тела - по две подвижные системы координат: естественная Oxyz и связанная с твердым телом Cx'y'z' (Cx', Cy', Cz' -главные центральные оси инерции). Все системы координат приняты правыми, причем считается, что оси ОЕ, , Ox , Cx' направлены слева
направо, а оси О)Z, Oz, Cz' - вниз (рис. 2). Оси естественной системы координат направлены соответственно по касательной, нормали и бинормали - к оси пути. Начало координат О для каждого твердого тела находится на расстоянии s от его положения в начальный момент времени (здесь s - пройденный путь). Положение естественной системы координат относительно неподвижной характеризуется дуговой координатой вдоль оси пути s, углом X между осями
Ox и О£, в плане, а также углами фК и 9h между этими осями в вертикальных плоскостях, которые определяются возвышением наружного рельса hr. Параметры рельса х, К - задан-
ные функции координаты s, для круговой кривой X и hr - постоянные, для прямой - равняются нулю [11].
Fig. 2. Coordinate system
Положение твердого тела относительно естественной системы координат описывается величинами x, y, z, у, ф, 9 . Поступательные перемещения x, y, z и углы поворота у, ф, 9 описывают подергивание, боковой относ, подпрыгивание и виляние, галопирование, боковую качку твердого тела соответственно [3].
При обозначении координат твердых тел системы введены следующие индексы: кузов -f , рама тележки - Ък (к =1, 2), колесная пара
- і (і = 1,6 - номер колесной пары по ходу движения экипажа), колесо - ij (j=1 - левая сторона вагона, j =2 - правая), рельсы в точках контакта - rij.
Для определения числа степеней свободы рассматриваемой механической системы приняты во внимание ограничения, наложенные на перемещения тел с учетом общепринятых допущений и конструктивных особенностей ходовых частей локомотива.
Вертикальные перемещения и боковая качка колесных пар выражаются через вертикальные перемещения колес:
z-i + z^ z-, — z л —
z = -Vі2; 9i = -A—2; (i 1,6), (1)
2 2d1
где 2d1 - расстояние между средними кругами катания колесной пары.
Предполагается, что радиусы r кругов ката-
doi: 10.15802/stp2015/46078
© С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
188
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
ния всех колес равны между собой. Тогда при совпадении продольных плоскостей симметрии пути и тележки при перемещении вагона вдоль оси пути все колеса при отсутствии проскальзывания повернутся на один и тот же угол:
Фі =- , (i = 1,6), (2)
r
tj =1 m [(s+xj_ у і x і )2+(yj+x x j)2 +
+(zі - )2 ] + 21„ (0і +0)2 +
+2 ijv і +% і )2 + - і,уііфі+(ф¥ )2, (6)
где s - текущее значение пройденного пути.
В продольном направлении кузов тепловоза и рамы тележек совершают одинаковые линейные перемещения:
Xf = % = ХЬ2. (3)
С учетом связей (1) - (3) система «локомо-тив-рельсовая колея» имеет 78-20=58 степеней свободы.
Уравнения движения локомотива составлены в форме уравнений Лагранжа II рода [3]:
А (<Г_) _дф дП дф = Q
dt dq, dq, dq, dq, 1 ' ’
(i = 1,..., 58), (4)
где mj - масса j -го твердого тела; Ij с соответствующими индексами обозначают главные центральные моменты инерции j -го тела;
X j = VKj, Kj - кривизна пути под j -ым телом; V - скорость движения; hj - возвышение пути под центром масс j -го тела, вызванное возвышением в кривой наружного рельса на величину hr =0h 2d1.
Кинетическая энергия системы, моделирующей путь, в соответствии с принятыми допущениями записывается следующим образом:
тп = 2mrhЕЕ4 +2mrvЕЕ4, (7)
2 і=1 і=1 2 і=1 і=1
где T и П - кинетическая и потенциальная энергии системы; Ф - функция рассеяния; qt, qt - обобщенные координаты и их скорости; Qi - соответствующие им обобщенные силы, не имеющие потенциала; S, - внешние силы.
Кинетическая энергия системы получена путем сложения кинетических энергий тел, входящих в систему, определенных по теореме Кенига:
т = Tf +ЕТЪк +Ет + Тп , (5)
к=1 1=1
где Tf - кинетическая энергия кузова тепловоза; Tbk - кинетическая энергия к -ой рамы тележки; Ti - кинетическая энергия i -ой колесной пары; Тп - кинетическая энергия пути.
В общем случае выражение для вычисления кинетической энергии для j -го твердого тела имеет следующий вид:
где mrh mrv - массы пути, приведенные к одному колесу в горизонтальном поперечном и вертикальном направлениях.
Потенциальная энергия системы (П) равна сумме потенциала сил тяжести (ПЬ1), энергии сжатия упругих элементов подвешивания (Пь 2) и потенциальной энергии приведенных масс пути (ПП ).
Потенциальная энергия, обусловленная опусканием центров тяжести тел, определяется из выражения:
2
ПL1 = _ mfS (0к/ьУ/ + Zf ) _ mbg Е0>*кУък _
к=1
6
+zbk) _ mkpg Е (0hjy, +z, ), (8)
і=1
где mf, mb, mkp - массы кузова, рамы тележки и колесной пары соответственно; g - ускорение свободного падения; 0h]- - угол между
горизонтальной плоскостью под центром масс j -го тела и плоскостью пути, вызванный возвышением наружного рельса над внутренним
doi: 10.15802/stp2015/46078
© С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
189
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
в криволинейном участке пути.
Составляющая потенциальной энергии nL 2 определяется по теореме Клапейрона как сумма энергий, накапливаемых в упругих элементах системы при их деформациях. В общем случае nL 2 имеет вид:
nL 2 = 2(i kcz A^ + key A 2yi ) +
2 1=1
+i(VV(k A2 .. + k A2 .. + k A2..) +
^ v v ax axij ay ayij az azij'
2 2=1 j=1
+2 kz XAz., (9)
2 . =1
где Aczl, Acy1, - изменение длины l -ых упругих
элементов центрального подвешивания в вертикальном и горизонтальном поперечном направлениях; kcz, key - жесткости пружин центрального подвешивания соответственно в вертикальном и горизонтальном поперечном направлениях; Aaz,j , Ay] , Axj - деформации пружин буксового подвешивания в вертикальном, горизонтальном поперечном и продольном направлениях; kclz, key, kax - соответствующие жесткости пружин буксового подвешивания; kz - жесткость скользуна, установленного на раме тележки, ограничивающего вертикальные перемещения кузова; Azi - взаимное вертикальное перемещение кузова и рамы тележки в зоне i -го скользуна.
Выражение для определения потенциальной энергии пути имеет вид:
Пп =
и 4 ІІ(У A 2
-k A ..
ry rytj
),
i =1 j=1
(10)
где Arzrj = zij ; Aryjj = Угі. (i =1, 6 j =1,2).
Функция рассеивания энергии в экипаже учитывает действие сил вязкого трения в обеих ступенях подвешивания и при взаимных поворотах кузова относительно ходовых частей в плане:
1 4 2
Фт = -(В VA2. +в YA2 i) +
L ~ \rcz / у czi rcy / у cy /
2 i=1 i=1
+1 Pez yA2i + 2 2i, (11)
2 i =1 2 i =1
где Pcz, Pcy - коэффициенты вязкого трения
в центральном подвешивании в вертикальном и горизонтальном поперечном направлениях соответственно; Pez - коэффициент вязкого трения в буксовом подвешивании в вертикальном направлении; Р2 - коэффициент вязкого
трения при взаимных поворотах кузова относительно рамы тележки в плане; A - относительные скорости перемещений тел, связанных диссипативными элементами.
Относительные скорости A перемещений тел, связанных диссипативными элементами, определяются как производные по времени от относительных перемещений соответствующих тел.
Функция рассеивания энергии в пути имеет следующий вид:
Фп = 2 УУСРгу A lij +Prz A2zij ), (12)
2 i =1 j =1
где Pry, Prz - коэффициенты рассеивания энергии в пути в горизонтальном и вертикальном поперечном направлениях.
Обобщенные силы Q. являются силами взаимодействия колес и рельсов. Действие сил в контакте колес и рельсов (силы псевдоскольжения) моделируется по методике, описанной в работе [5].
Составляющие сил крипа определяются по формулам:
'т' __ xV 74 .
1 xij- _ Jij ;
ьї
s sec a —
Taij =—-------j T. (i = 1,6, j = 1,2), (13)
sij
где Tj - проекция полной силы Ту на продольную ось x; Таі. - проекция силы Т. на касательную к профилю поверхности катания ко-
doi: 10.15802/stp2015/46078
© С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
190
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
леса в точке контакта; а. - угол между касательной к профилю поверхности колеса в точке контакта и поперечной осью у; в. - относительная скорость проскальзывания j -го колеса i -ой колесной пары по рельсу (крип); s, вyij. -
продольная и поперечная составляющие относительной скорости проскальзывания j -го колеса i -ой колесной пары по рельсу;
T = -F в..
j j j
Fij в
ij ij \4
ftrNij
)4 +1
(14)
где N. - нормальное давление j -го колеса і -ой колесной пары на рельс; ftr - коэффициент трения скольжения; F. - коэффициент крипа j -го колеса i -ой колесной пары.
Полный крип в. определяется выражением:
1
в] = (в2j +в2j sec2 ац)2. (15)
Обобщенные силы Qi вычисляются как коэффициенты при вариациях обобщенных координат в выражениях возможных работ сил псевдоскольжения T j , Taij .
Внешними силами, действующими на локомотив, являются: сила тяги FT , развиваемая локомотивом, тормозная сила B, возникающая при включении тормозов, и силы сопротивления движению W, к которым относятся все остальные внешние силы.
В общем случае внешние силы, действующие на локомотив, могут быть представлены в виде:
S = W + B + FT . (16)
После подстановки в уравнение Лагранжа (4) выражений кинетической и потенциальной энергий, функции рассеивания, обобщенных и внешних сил получена система нелинейных дифференциальных уравнений, описывающая движение рассматриваемого локомотива по прямолинейным и криволинейным участкам пути.
На основе разработанной математической модели разработан алгоритм и создана компьютерная программа, предназначенная для
оценки динамических характеристик движущегося локомотива. Разработанная компьютерная программа позволяет имитировать широкий спектр различных расчетных случаев, охватывающих всевозможные сочетания параметров ходовых частей, вагона в целом и пути. В результате расчетов может быть получен практически любой набор необходимых выходных величин.
Расчетная оценка динамических показателей локомотива проведена путем анализа результатов решения нелинейных дифференциальных уравнений, представляющих динамическую модель колебаний экипажа. Решения нелинейных дифференциальных уравнений получены численно по методу Адамса-Баш-форта.
Ходовые качества, безопасность движения и плавность хода локомотива оценивались по следующим показателям [7]:
WG, WV - показатели плавности хода в горизонтальной и вертикальной плоскостях; Иг -рамные силы; kdv1, kclv2 - коэффициенты вертикальной динамики в первой и второй ступенях подвешивания; kst - коэффициенты запаса устойчивости против схода колеса с рельса.
Для сравнительного анализа также проведены расчеты по оценке динамических показателей серийного локомотива ТЭ116 [16], который широко эксплуатируется на железных дорогах Украины в настоящее время.
Движение локомотивов моделировалось с учетом действия внешних возмущений, связанных с геометрическим несовершенством рельсовой колеи (просадки, уширения, углы в плане, стыки, неравномерный износ). Возмущения, имеющие случайный характер, моделировались под колесами локомотива с транспортным запаздыванием, которое зависит от расстояния между колесными парами и скорости движения вагона. В качестве возмущений принимались случайные неровности рельсовой колеи в вертикальном и поперечном горизонтальном направлениях, соответствующие удовлетворительному состоянию пути [13].
Рассмотрено движение обоих локомотивов по пути, который включает в себя как прямолинейные, так и криволинейные участки пути (входные переходные кривые, круговые кривые
doi: 10.15802/stp2015/46078
© С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
191
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
разного радиуса, выходные переходные кривые).
При этом начальная скорость движения экипажа принималась равной 40 км/ч, а затем увеличивалась с шагом 10 км/ч до 170 км/ч на прямолинейном участке пути, до 120 км/ч на кривой радиусом 650 м и до 100 км/ч на кривой радиусом 300 м. Участок пути длиной 1 250 м включает в себя прямолинейный отрезок длиной 500 м, входную переходную кривую длиной 150 м, круговую кривую длиной 300 м, выходную переходную кривую длиной 150 м и еще один прямолинейный отрезок длиной 150 м. Радиусы круговых кривых составляли 650 и 300 м. При движении по круговой кривой радиусом 650 м (допустимое непогашенное ускорение согласно Нормам равно 0,7 м/с2 [7]) с возвышением наружного рельса 0,11 м наибольшая допускаемая скорость движения составляет 108 км/ч. Для круговой кривой радиусом 300 м с возвышением наружного рельса 0,15 м наибольшая допускаемая скорость движения составляет 80 км/ч.
При движении экипажа по круговой кривой дополнительным критерием оценки его динамических характеристик является коэффициент устойчивости экипажа от опрокидывания (kopr).
На рис. 3-8 приведены максимальные расчетные значения показателей плавности хода (WG, WV ), рамных сил (Hr), коэффициентов вертикальной динамики (км, kdv2) и коэффициентов запаса устойчивости против схода колес с рельсов (kst), полученные для рассматриваемых локомотивов. Как видно из приведенных рисунков, значения динамических показателей как тепловоза ТЭ116, так и локомотива типа ER20 по большинству нормируемых величин не превышают допустимых значений во всем исследованном диапазоне изменения скоростей движения.
При этом значения некоторых рассматриваемых динамических характеристик для обоих экипажей несущественно отличаются друг от друга. Существенно худшие показатели по коэффициентам вертикальной динамики обеих ступеней рессорного подвешивания имеет тепловоз ТЭ116.
Рис. 3. Зависимости максимальных значений показателя плавности хода в горизонтальном направлении от скорости движения локомотивов ER20 и ТЭ116 по прямым участкам пути
Fig. 3. Dependences of maximum values of smoothness in the horizontal direction on the speed of the locomotive ER20 and TE116 on straight sections of the track
Рис. 4. Зависимости максимальных значений показателя плавности хода в вертикальном направлении от скорости движения локомотивов ER20 и ТЭ116 по прямым участкам пути
Fig. 4. Dependences of maximum values of smoothness in the vertical direction on the speed of the locomotive ER20 and TE116 on straight sections of the track
Рис. 5. Зависимости максимальных значений рамных сил от скорости движения локомотивов ER20 и ТЭ116 по прямым участкам пути
Fig. 5. Dependences of maximum values of frame forces on the speed of the locomotive ER20 and TE116 on straight sections of the track
doi: 10.15802/stp2015/46078 © С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
192
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
Рис. 6. Зависимости минимальных значений коэффициента запаса устойчивости против схода колеса с рельсов от скорости движения локомотивов ER20 и ТЭ116 по прямым участкам пути
Fig. 6. Dependences of the minimum values of safety factor against derailment of the wheel from the rail on the speed of locomotive ER20 and TE116 on straight
sections of the track
Рис. 7. Зависимости максимальных значений коэффициентов вертикальной динамики первой ступени подвешивания от скорости движения локомотивов ER20 и ТЭ116 по прямым участкам пути
Fig. 7. Dependences of maximum values of the coefficients of the vertical dynamics of the first stage hanging on the speed of the locomotive ER20 and TE116 on straight sections of the track
Рис. 8. Зависимости максимальных значений коэффициентов вертикальной динамики второй ступени подвешивания от скорости движения локомотивов ER20 и ТЭ116 по прямым участкам пути
Fig. 8. Dependences of maximum values of the coefficients of the vertical dynamics of the second stage hanging on the speed of the locomotive ER20 and TE116 on straight sections of the track
На рис. 9-18 приведены графики зависимостей максимальных значений рамных сил, коэффициента устойчивости экипажа от опрокидывания, коэффициентов вертикальной динамики, коэффициента запаса устойчивости против схода колеса с рельсов для тепловозов ТЭ116 и ER20, движущихся по круговым кривым радиуса 650 м (рис. 9-13) и радиуса 300 м (рис. 14-18) от скорости движения.
Рис. 9. Зависимости максимальных значений рамных сил от скорости движения локомотивов ER20 и ТЭ116 по круговой кривой радиуса 650 м
Fig. 9. Dependences of maximum values of frame forces on the speed of the locomotive ER20 and TE116 along a circular curve, radius of 650 m
Рис. 10. Зависимости минимальных значений коэффициента устойчивости экипажа от опрокидывания от скорости движения локомотивов ER20 и ТЭ116 по круговой кривой радиуса 650 м
Fig. 10. Dependences of the minimum values of stability factor against overturning of the crew on the speed of locomotive ER20 and TE116 on a circular curve, radius of 650 m
doi: 10.15802/stp2015/46078 © С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
193
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
Рис. 11. Зависимости максимальных значений коэффициентов вертикальной динамики первой ступени подвешивания от скорости движения локомотивов ER20 и ТЭ116 по круговой кривой радиуса 650 м
Fig. 11. Dependences of maximum values of the coefficients of the vertical dynamics of the first stage hanging on the speed of the locomotive ER20 and TE116 along a circular curve, radius of 650 m
Рис. 12. Зависимости максимальных значений коэффициентов вертикальной динамики второй ступени подвешивания от скорости движения локомотивов ER20 и ТЭ116 по круговой кривой радиуса 650 м
Fig. 12. Dependences of maximum values of the coefficients of the vertical dynamics of the second stage hanging on the speed of the locomotive ER20 and TE116 along a circular curve, radius of 650 m
Рис. 13. Зависимости минимальных значений коэффициента запаса устойчивости против схода колеса с рельсов от скорости движения локомотивов ER20 и ТЭ116 по круговой кривой радиуса 650 м
Fig. 13. Dependences of the minimum values of safety factor against derailment of the wheel from the rail on the speed of locomotive ER20 and TE116 on a circular curve, radius of 650 m
Рис. 14. Зависимости максимальных значений рамных сил от скорости движения локомотивов ER20 и ТЭ116 по круговой кривой радиуса 300 м
Fig. 15. Dependences of maximum values of frame forces on the speed of the locomotive ER20 and TE116 along a circular curve, radius 300 m
Рис. 15. Зависимости минимальных значений коэффициента устойчивости экипажа от опрокидывания от скорости движения локомотивов ER20 и ТЭ116 по круговой кривой радиуса 300 м
Fig. 15. Dependences of the minimum values of stability factor against overturning the crew on the speed of locomotive ER20 and TE116 in a circular curve, radius of 300 m
Рис. 16. Зависимости максимальных значений коэффициентов вертикальной динамики первой ступени подвешивания от скорости движения локомотивов ER20 и ТЭ116 по круговой кривой радиуса 300 м
Fig. 16. Dependences of maximum values of the coefficients of the vertical dynamics of the first stage hanging on the speed of the locomotive ER20 and TE116 along a circular curve, radius of 300 m
doi: 10.15802/stp2015/46078
© С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
194
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
Рис. 17. Зависимости максимальных значений коэффициентов вертикальной динамики второй ступени подвешивания от скорости движения локомотивов ER20 и ТЭ116 по круговой кривой радиуса 300 м
Fig. 17. Dependences of maximum values of the coefficients of the vertical dynamics of the second stage hanging on the speed of the locomotive ER20 and TE116 along a circular curve, radius of 300 m
Рис. 18. Зависимости минимальных значений коэффициента запаса устойчивости против схода колеса с рельсов от скорости движения локомотивов ER20 и ТЭ116 по круговой кривой радиуса 300 м
Fig. 18. Dependences of the minimum values of safety factor against derailment of the wheel from the rail on the speed of locomotive ER20 and TE116 in a circular curve, radius of 300 m
Проведенные исследования показали, что при всех рассмотренных скоростях движения, для которых допустимо движение по описанным выше круговым кривым, нормируемые динамические показатели не превышают допустимых значений.
Научная новизна и практическая значимость
Разработана математическая модель пространственного движения шестиосного локомотива перспективной конструкции, учитывающая особенности его конструкции, которая может быть использована при проведении научных исследований по созданию перспективных конструкций локомотивов.
Дана оценка динамических характеристик локомотива перспективной конструкции при его движении по пути произвольного очертания.
Выводы
Оценка динамических характеристик локомотива типа ER20 показала возможность использования его на железнодорожной колее шириной 1520 мм при движении со скоростями до 140 км/ч. Естественно, это требует последующей экспериментальной проверки.
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ
1. Грузовой тепловоз ER20 CF для Литвы // Ж.-д. мира. - 2009. - № 5. - С. 17-24.
2. Короп, Г. В. Разработка и исследование математической модели движения шестиосного локомотива с управляемой радиальной установкой колесных пар / Г. В. Короп, С. А. Клюев, Е. В. Зубарь // Вісн. СНУ ім. В. Даля. - 2013. -№ 4 (193), ч. 2. - C. 65-68.
3. Лазарян, В. А. Динамика транспортных средств / В. А. Лазарян. - Киев : Наук. думка, 1985. - 528 с.
4. Маслиев, В. Г. К 50-летию создания первого отечественного тепловоза на пневматическом рессорном подвешиваннии / В. Г. Маслиев // Вісн. НТУ «ХПІ». Сер. : Трансп. машинобудування. -Харків, 2013. - № 32 (1005). - С. 60-66.
5. Математическое моделирование колебаний рельсовых транспортных средств / В. Ф. Уш-калов, Л. М. Резников, В. С. Иккол [и др.]. -Киев : Наук. думка, 1989. - 240 с.
6. Мямлин, С. В. Совершенствование конструкции ходовых частей локомотивов / С. В. Мям-лин, Л. А. Недужая // Наука та прогрес трансп. Вісн. Дніпропетр. нац. ун-ту. залізн. трансп. -
2013. - № 5 (47). - C. 124-136. doi: 10.15802/-stp2013/17977.
7. Нормы для расчета и оценки прочности несущих элементов, динамических качеств и воздействия на путь экипажной части локомотивов железных дорог МПС РФ колеи 1520 мм. -Москва : МПС РФ : ВНИИВ, 1998. - 145 с.
8. Обновление локомотивного парка Литовских железных дорог / С. Дайлидка, С. В. Мямлин, Л. П. Лингайтис [и др.] : зб. наук. пр. / ДонІЗТ. - Донецк, 2011. - Вип. 28. - С. 174-179.
9. Пономарева, Е. Тепловоз ER20 CF : особенности конструкции / Е. Пономарева // Локомотив-
информ. - 2012. - № 2. - С. 12-17.
doi: 10.15802/stp2015/46078 © С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
195
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
10. Програма оновлення локомотивного парку
залізниць України на 2012-2016 роки : постанова Кабінету Міністрів України від 1 серпня 2011 р. №840 [Электронный ресурс]. - Режим доступа: http://zakon4.rada.gov.ua/laws/show-
/840-2011-%D0%BF. - Загл. с экрана. - Проверено : 31.03.2015.
11. Радченко, Н. А. Криволинейное движение рельсовых транспортных средств / Н. А. Радченко. - Киев : Наук. думка, 1988. - 212 с.
12. Стан локомотивного парку Укрзалізниці // Магистраль. - 2012. - № 20. - С. 2.
13. Ушкалов, В. Ф. Расчетные возмущения для оценки динамических качеств грузовых вагонов / В. Ф. Ушкалов, Л. Г. Лапина, И. А. Ма-щенко // Наука та прогрес трансп. Вісн. Дніп-ропетр. нац. ун-ту. залізн. трансп. - 2013. -
№ 4 (46). - C. 135-144. doi: 10.15802/stp2013-/16600.
14. Myamlin, S. Construction Analysis of Mechanical Parts of Locomotives / S. Myamlin, M. Luchanin, L. Neduzha // TEKA Commission of Motorization and Power Industry in Agriculture. - 2013. -Vol. 13, № 3. - P. 162-169.
15. Myamlin, S. Mathematical Modeling of a Cargo Locomotive / S. Myamlin, S. Dailidka, L. Neduzha // Transport Means : Proc. of 16th Intern. Conf. / Kaunas : Kaunas University of Technology, 2012. - Р. 310-312.
16. Ushkalov, V. F. The influence of locomotives state on their wheelset wear under operation conditions / V. F. Ushkalov, O. M. Markova, E. N. Kovtun // Supplement to Vehicle System Dynamics. - 1998. - Vol. 29. - P. 18-29. doi: 10.1080/00423119808969549.
С. В. МЯМЛІН1, Н. Ю. НАУМЕНКО2, О. М. МАРКОВА3, О. М. КОВТУН4, В. В. МАЛИЙ5*
*Каф. «Вагони і вагонне господарство», Дніпропетровський національний університет залізничного транспорту імені академіка В. Лазаряна, вул. Лазаряна, 2, Дніпропетровськ, Україна, 49010, тел. +38 (056) 776 84 98, ел. пошта sergeymyamlin@gmail.com, OrCID 0000-0002-7383-9304
2Відділ «Динаміка багатовимірних механічних систем», Інститут технічної механіки Національної академії наук України і Державного космічного агентства України, вул. Лешко-Попеля, 15, Дніпропетровськ, Україна, 49005, тел. +38 (056) 375 68 01, ел. пошта dep7@ukr.net, ORCID 0000-0002-6913-2568
3Відділ «Динаміка багатовимірних механічних систем», Інститут технічної механіки Національної академії наук України і Державного космічного агентства України, вул. Лешко-Попеля, 15, Дніпропетровськ, Україна, 49005, тел. +38 (056) 375 68 03, ел. пошта dep7@ukr.net, ORCID 0000-0002-9390-4046
4Відділ «Динаміка багатовимірних механічних систем», Інститут технічної механіки Національної академії наук України і Державного космічного агентства України, вул. Лешко-Попеля, 15, Дніпропетровськ, Україна, 49005, тел. +38 (0562) 46 79 35, ел. пошта dep7@ukr.net, ORCID 0000-0002-2961-3682
5*Відділ «Динаміка багатовимірних механічних систем», Інститут технічної механіки Національної академії наук України і Державного космічного агентства України, вул. Лешко-Попеля, 15, Дніпропетровськ, Україна, 49005, тел. +38 (0562) 46 79 35, ел. пошта dep7@ukr.net, ORCID 0000-0002-6501-1398
ОЦІНКА ДИНАМІЧНИХ ПОКАЗНИКІВ ШЕСТИВІСНОГО ЛОКОМОТИВА
Мета. Робота присвячена оцінці динамічних властивостей локомотива перспективної конструкції та визначенню можливості його використання на залізничних коліях України. Методика. Для досягнення поставленої мети використані методи математичного та комп’ютерного моделювання динаміки залізничних транспортних засобів, а також методи числового інтегрування систем звичайних нелінійних рівнянь. Результати. Для розв’язання поставленої задачі побудована розрахункова схема локомотива на тривісних візках, яка представляє собою систему твердих тіл, з’єднаних элементами різної реології. З використанням цієї розрахункової схеми розроблено математичну модель просторового руху локомотива, яка дозволяє проводити дослідження його просторових коливань при русі по прямолінійним та криволінійним ділянкам колії із випадковими нерівностями в плані та профілі. При складанні математичної моделі приймались до уваги як геометричні (нелінійність профілю кочення колеса), так і фізичні нелінійності системи (робота сил сухого тертя, нелінійність характеристик взаємодії коліс та рейок). За допомогою комп’ ютерного моделювання проведено багатоваріантні розрахунки з метою оцінки динамічних властивостей локомотива при його русі по прямолінійним та криволінійним ділянкам колії. При цьому визначались показники плавності ходу локомотива в горизонтальній та вертикальній площинах, рамні сили, коефіцієнти вертикальної динаміки в першій та другій ступенях підвішування, коефіцієнти запасу стійкості проти сходу коліс із рейок. Крім того, проведено порівняння отриманих результатів із аналогічними характеристиками широко використовуваного
doi: 10.15802/stp2015/46078 © С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
196
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
на залізницях України шестивісного тепловоза ТЭ116. Визначено вплив швидкості руху та технічного стану залізничної колії на безпеку руху локомотива. Наукова новизна. Розроблена математична модель просторового руху шестивісного локомотива перспективної конструкції, що враховує особливості його конструкції, та наведена оцінка динамічних характеристик при русі екіпажа по колії довільного обрису. Практична значимість. Розроблена математична модель руху шестивісного локомотива може бути використана при проведені наукових досліджень по створенню перспективных конструкцій локомотивів.
Ключові слова: рухомий склад; локомотив; математична модель; динамічні показники
S. V. MYAMLIN1, N. YU. NAUMENKO2, O. M. MARKOVA3, H. M. KOVTUN4,
V. V. MALIY5*
*Dep. «Car and Car Facilities», Dnipropetrovsk National University of Railway Transport named after Academician V. Lazaryan, Lazaryan St., 2, Dnipropetrovsk, Ukraine, 49010, tel./fax +38 (056) 776 84 98, e-mail sergeymyamlin@gmail.com,
ORCID 0000-0002-7383-9304
2Dep. «Dynamics of Multidimensional Mechanical Systems», Institute of Technical Mechanics of the National Academy of Sciences of Ukraine and State Space Agency of Ukraine, Leshko-Popelia St., 15, Dnipropetrovsk, Ukraine, 49005, tel. +38 (056) 375 68 01, e-mail dep7@ukr.net, ORCID 0000-0002-6913-2568
3Dep. «Dynamics of Multidimensional Mechanical Systems», Institute of Technical Mechanics of the National Academy of Sciences of Ukraine and State Space Agency of Ukraine, Leshko-Popelia St., 15, Dnipropetrovsk, Ukraine, 49005, tel. +38 (056) 375 68 03, e-mail dep7@ukr.net, ORCID 0000-0002-9390-4046
4Dep. «Dynamics of Multidimensional Mechanical Systems», Institute of Technical Mechanics of the National Academy of Sciences of Ukraine and State Space Agency of Ukraine, Leshko-Popelia St., 15, Dnipropetrovsk, Ukraine, 49005, tel. +38 (0562) 46 79 35, e-mail dep7@ukr.net, ORCID 0000-0002-2961-3682
5 Dep. «Dynamics of Multidimensional Mechanical Systems», Institute of Technical Mechanics of the National Academy of Sciences of Ukraine and State Space Agency of Ukraine, Leshko-Popelia St., 15, Dnipropetrovsk, Ukraine, 49005, tel. +38 (0562) 46 79 35, e-mail dep7@ukr.net, ORCID 0000-0002-6501-1398
EVALUATION OF DYNAMIC INDICATORS OF SIX-AXLE LOCOMOTIVE
Purpose. The paper is devoted to dynamic characteristics evaluation of the locomotive with prospective design and determination the feasibility of its use on the Ukrainian railways. Methodology. The methods of mathematical and computer modeling of the dynamics of railway vehicles, as well as methods for the numerical integration of systems of ordinary nonlinear differential equations were used to solve the problem. Findings. The calculated diagram of a locomotive on three-axle bogies was built to solve the problem, and it is a system of rigid bodies connected by various elements of rheology. The mathematical model of the locomotive movement, allowing studying its spatial vibrations at driving on straight and curved sections of the track with random irregularities in plan and profile was developed with use of this calculated diagram. At compiling the mathematical model took into account both geometric (nonlinearity profile of the wheel roll surface) and physical nonlinearity of the system (the work forces of dry friction, nonlinearity characteristics of interaction forces between wheels and rails). The multivariate calculations, which allowed assessing the dynamic qualities of the locomotive at its movement along straight and curved sections of the track, were realized with the use of computer modeling. The smoothness movement indicators of the locomotive in horizontal and vertical planes, frame strength, coefficients of vertical dynamics in the first and second stages of the suspension, the load factor of resistance against the derailment of the wheel from the rail were determined at the period of research. In addition, a comparison of the obtained results with similar characteristics is widely used on the Ukrainian railways in six-axle locomotive TE 116. The influence of speed and technical state of the track on the locomotive traffic safety was determined. Originality. A mathematical model of the spatial movement of a six-axle locomotive with prospective design, taking into account the peculiarities of its design was developed. Evaluation of dynamic characteristics at the carriage movement along a track of arbitrary shape was made. Practical value. The developed mathematical model of the six-axle locomotive movement can be used at conducting research on development of advanced locomotives designs.
Keywords: rolling stock; locomotive; mathematical model; dynamic indicators
REFERENCES
1. Gruzovoy teplovoz ER20 CF dlya Litvy [Freight locomotive ER20 CF for Lithuania]. Zheleznyye dorogi mira - Railways of the World, 2009, no. 5, pp. 17-24.
doi: 10.15802/stp2015/46078 © С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
197
ISSN 2307-3489 (Print), ISSN 2307-6666 (Online)
Наука та прогрес транспорту. Вісник Дніпропетровського національного університету залізничного транспорту, 2015, № 3 (57)
РУХОМИИ СКЛАД ЗАЛІЗНИЦЬ І ТЯГА ПОЇЗДІВ
2. Korop G.V., Klyuyev S.A., Zubar Ye.V. Razrabotka i issledovaniye matematicheskoy modeli dvizheniya shes-tiosnogo lokomotiva s upravlyayemoy radialnoy ustanovkoy kolesnykh par [Development and research of mathematical model of six-axle locomotive with a controlled radial installation of wheel pairs]. Visnyk Skhidnoukrainskoho natsionalnoho universytetu imeni V. Dalia [Bulletin of the East Ukrainian National University named after V. Dahl], 2013, no. 4 (193), part 2, pp. 65-68.
3. Lazaryan V.A. Dinamika transportnykh sredstv [The dynamics of vehicles]. Kyiv, Naukova dumka Publ.,
1985. 528 p.
4. Masliyev V.G. K 50-letiyu sozdaniya pervogo otechestvennogo teplovoza na pnevmaticheskom ressornom podveshivannii [On the 50th anniversary of the first Russian locomotive on the pneumatic spring suspension]. Visnyk Natsionalnoho Tekhnichnoho Universytetu Kharkivskoho Politekhnichnoho Instytutu [Bulletin of National Technical University of Kharkiv Polytechnic Institute], 2013, no. 32 (1005), pp. 60-66.
5. Ushkalov V.F., Reznikov L.M., Ikkol V.S. Matematicheskoye modelirovaniye kolebaniy relsovykh transportnykh sredstv [Mathematical modeling of vibrations of railway vehicles]. Kyiv, Naukova dumka Publ., 1989. 240 p.
6. Myamlin S.V., Neduzhaya L.A. Sovershenstvovaniye konstruktsii khodovykh chastey lokomotivov [Design
improvement of the locomotive running gears]. Nauka ta prohres transportu. Visnyk Dnipropetrovskoho natsionalnoho universytetu zaliznychnoho transportu - Science and Transport Progress. Bulletin of Dni-propetrovsk National University of Railway Transport, 2013, no. 5 (47), pp. 124-136. doi:
10.15802/stp2013/17977.
7. Normy dlya rascheta i otsenki prochnosti nesushchikh elementov, dinamicheskikh kachestv i vozdeystviya na put ekipazhnoy chasti lokomotivov zheleznykh dorog MPS RF kolei 1520 mm [Standards for the calculation and assessment of the strength of load-bearing elements, dynamic qualities and impact on the track of undercarriage locomotives Ministry of Railways of the Russian 1,520 mm]. Moscow, Ministerstvo putey soob-shcheniya Rossiyskoy Federatsii, Vsesoyuznyy nauchno-issledovatelskiy institut vagonostroeniya Publ., 1998. 145 p.
8. Daylidka S., Myamlin S.V., Lingaytis L.P., Neduzhaya L.A., Yastremskas V. Obnovleniye lokomotivnogo parka Litovskikh zheleznykh dorog [Upgrading the locomotive fleet of Lithuanian Railways]. Zbirnyk naukovykh prats Donetskoho instytutu zaliznychnoho transportu [Proc. of Donetsk Institute of Railway Transport], 2011, vol. 28, pp. 174-179.
9. Ponomareva Ye. Teplovoz ER20 CF: osobennosti konstruktsii [The locomotive ER20 CF: design features].
Lokomotiv-Inform - Lokomotiv-Inform, 2012, no. 2, pp. 12-17.
10. Prohrama onovlennia lokomotyvnoho parku zaliznyts Ukrainy na 2012-2016 roky [The upgrade program for railways locomotive park of Ukraine for 2012-2016]. Available at: http://zakon4.rada.gov.ua/laws/show/840-2011-%D0%BF (Accessed 31. 03.2015).
11. Radchenko N.A. Krivolineynoye dvizheniye relsovykh transportnykh sredstv [Curvilinear motion of railway vehicles]. Kyiv, Naukova dumka Publ., 1988. 212 p.
12. Stan lokomotyvnoho parku Ukrzaliznytsi [The state of railways locomotive fleet]. Magistral - Magistral, 2012, no. 20, p. 2.
13. Ushkalov V.F., Lapina L.G., Mashchenko I.A. Raschetnyye vozmushcheniya dlya otsenki dinamicheskikh kachestv gruzovykh vagonov [Calculated disturbances for evaluation of dynamical properties of freight cars]. Nauka ta prohres transportu. Visnyk Dnipropetrovskoho natsionalnoho universytetu zaliznychnoho transportu - Science and Transport Progress. Bulletin of Dnipropetrovsk National University of Railway Transport, 2013, no. 4 (46), pp. 135-144. doi: 10.15802/stp2013/16600.
14. Myamlin S., Luchanin M., Neduzha L. Construction Analysis of Mechanical Parts of Locomotives. TEKA Commission of Motorization and Power Industry in Agriculture, 2013, vol. 13, no. 3, pp. 162-169.
15. Myamlin S., Dailidka S., Neduzha L. Mathematical Modeling of a Cargo Locomotive. Proc. of 16th Intern. Conf. Transport Means 2012, Kaunas, 2012, pp. 310-312.
16. Ushkalov V.F., Markova O.M., Kovtun E.N. The influence of locomotives state on their wheelset wear under operation conditions. Supplement to Vehicle System Dynamics, 1998, vol. 29, pp. 18-29. doi: 10.1080/00423119808969549.
Статья рекомендована к публикации д.т.н. В. Л. Горобцом (Украина)
Поступила в редколлегию 03.02.2015
Принята к печати 20.04.2015
doi: 10.15802/stp2015/46078
© С. В. Мямлин, Н. Е. Науменко, О. М. Маркова, Е. Н. Ковтун, В. В. Малый, 2015
198