УДК 656.212.5 (23.01)
В. І. БОБРОВСЬКИЙ, Є. Б. ДЕМЧЕНКО (Дніпропетровський національний університет залізничного транспорту імені академіка В. Лазаряна)
МОДЕЛЮВАННЯ ПРОЦЕСУ НАСУВУ ТА РОЗПУСКУ СОСТАВІВ НА СОРТУВАЛЬНІЙ ГІРЦІ
Розроблено імітаційну модель процесу насуву та розпуску составів на сортувальній гірці. Запропонована модель передбачає спільне імітування режиму руху маневрового составу та процесу скочування його відче-пів. Це дозволить комплексно оцінювати якість сортувального процесу, що необхідно для визначення раціональних режимів розформування составів.
Ключові слова: сортувальна гірка, режим розформування составів, гірковий локомотив
Разработана имитационная модель процесса надвига и роспуска составов на сортировочной горке. Предложенная модель предусматривает совместное имитирование режима движения маневрового состава и процесса скатывания его отцепов. Это позволит комплексно оценивать качество сортировочного процесса, что необходимо для определения рациональных режимов расформирования составов.
Ключевые слова: сортировочная горка, режим расформирования составов, горочный локомотив
The simulation model of pushing and breaking-up process at the hump yard was developed. The proposed model provides co-simulation of train motion regime and cut rolling down. It allows to evaluate comprehensively the quality of classification process, which is necessary for determination of efficient mode of breaking-up process.
Key words: hump yard, breaking-up mode, hump engine
В сучасних умовах одним з основних напрямків підвищення ефективності функціонування сортувальних станцій є впровадження ресурсозберігаючих технологій переробки ва-гонопотоків. При цьому вирішення проблеми зниження енерговитрат в підсистемі розформування вимагає комплексного розгляду процесів насуву та розпуску составів. Проте, як показав аналіз [1], в теперішній час задачі удосконалення процесу насуву вирішуються окремо від розпуску составів. Так, існуючі моделі насуву [2, 3] дають можливість встановити залежність витрат палива гірковими локомотивами від швидкості розпуску, конструкції поздовжнього профілю сортувальної гірки та параметрів составу, що розформовується. В той же час, не враховуючи необхідність системному підходу до підвищення ефективності сортувального процесу, вказані моделі не дозволяють оцінювати вплив обраного режиму розформування составу на умови інтервального та прицільного гальмування його відчепів. Крім того, наявні моделі насуву побудовані з використанням спрощеного алгоритму керування гірковим локомотивом, що не дозволяє детально імітувати процес руху маневрового составу.
З метою вирішення вказаної задачі було розроблено модель процесу насуву та розпуску составів на сортувальній гірці. Вказана модель враховує особливості роботи гіркових тепловозів і дає можливість визначати точну початкову швидкість відчепів в моменти відриву від сос-
таву на вершині гірки (ВГ), що необхідно для подальшого моделювання процесу їх скочування. Крім того, розроблена модель насуву та розпуску дозволяє визначати витрати палива гірковим локомотивом, величина яких необхідна для визначення раціональних режимів розформування составів.
В моделі состав, що насувається на гірку, представляється сукупністю відчепів з певними параметрами (кількість вагонів, їх тип, довжина, маса та основний питомий опір руху) і розглядається як нерозтяжний гнучкий стержень з рівномірно розподіленою по довжині масою. Така модель составу дозволяє найбільш повно врахувати зміни умов його руху при переходах з одного елемента профілю на інший та після відриву від составу чергового відчепу.
Для вирішення задачі моделювання насуву та розпуску достатньо розглядати керований поступальний рух маневрового составу, тому при його імітуванні необхідно враховувати лише зовнішні сили, що співпадають з напрямком руху або протилежні йому. У відповідності з цим враховувались сили: Fx - дотична сила тяги локомотива; W - сила опору руху составу; Вт - гальмова сила локомотива. В рівняннях руху маневрового составу розглядаються відповідні питомі сили fK, Wu, Ьт; рівнодіюча вказаних сил f залежить від режиму роботи гірково-го локомотива та дорівнює f = f ± wR в режимі тяги, f = ± wR - в режимі вибігу та f = ± wR - Ьт - в режимі гальмування.
© Бобровський В. І., Демченко Є. Б. 2012
13
В роботі розроблено методику розрахунку сил, що діють на маневровий состав в процесі насуву та розпуску. Як показав аналіз [1], загальноприйнята методика визначення дотичної сили тяги локомотива /к, запропонована в ПТР [4], не відповідає реальним умовам виконання маневрової роботи на сортувальній гірці. У зв’язку з цим, дотична сили тяги локомотива / в розробленій моделі визначається за частковими (проміжними) тяговими характеристиками, які можуть бути реалізовані за умовами зчеплення, і далі по наступним проміжним характеристикам аж до виходу на автоматичну (зовнішню) характеристику [5].
Питома сила опору руху wx визначається як
= W0 + W” + WCB + WGR + wi + ^вр,
де w'0 - основний питомий опір руху локомотива;
W0 - основний питомий опір руху вагонів составу;
wсв - додатковий питомий опір від середовища та вітру;
Wck - додатковий питомий опір від стрілок та кривих;
Wj - додатковий питомий опір від ухилу;
w^ - додатковий питомий опір при зрушенні з місця.
Основний питомий опір руху локомотива W0 та додатковий питомий опір при зрушенні составу з місця w^ визначаються за методикою ПТР [4]. Основний питомий опір руху вагонів составу w0 розраховується як середньозважена величина основного питомого опору руху його відчепів і корегується після відриву від составу чергового відчепу; при цьому величина w”” для кожного вагону входить в структуру моделі составу. Величини додаткових опорів руху від середовища та вітру wCI!, стрілок та кривих wCl5 розраховуються за методикою [6].
Додатковий питомий опір від ухилу wі приймається чисельно рівним середньому ухилу і, на якому знаходиться маневровий состав; при цьому на підйомі значення і додатне (і > 0), а на спуску - від'ємне (і < 0). Середній ухил і, на якому знаходиться состав, коли його голова перебуває в точці Sj, визначається за різницею відміток голови h(Sj) та хвоста h(Sj - lC) составу
h(Sj)-h(Sj -lc)
i(Sj) =--^-----------’
(c
де lC - довжина маневрового составу, м.
У зв’язку з тим, що в процесі розформуван-
ня на сортувальній гірці автоматичні гальма составу, як правило, не включаються, зменшення його швидкості здійснюється лише за рахунок гальмівної сили маневрового локомотива, питоме значення Ьт якої розраховується за наступною формулою [4]:
Ьт = 1000 #р^кр, (1)
де Ор - розрахунковий гальмівний коефіцієнт;
фкр - розрахунковий коефіцієнт тертя гальмівних колодок.
Значення величин та фк^ визначаються за
допомогою методики [7] в залежності від кількості гальмівних осей локомотива, а у випадках включення гальм в составі - і від кількості осей вагонів, в яких задіяні автогальма. При цьому слід зауважити, що гальмівна сила Ьт з моменту переведення крана машиніста в гальмове положення зростає поступово до свого максимального значення. Крім того, при моделюванні гальмування маневрового составу слід враховувати час, потрібний на реакцію машиніста. В цьому зв'язку, розрахунковий гальмівний коефіцієнт Ор в (1) розглядається як функція від
тривалості гальмування їгал, значення якої приймається згідно з [7].
Особливістю моделювання розпуску составу є зміна його параметрів при відриві відчепів. В цьому зв’язку необхідно на кожному кроці At контролювати можливість відриву чергового відчепу. Після фіксування відриву відчепу в моделі виконується відповідне зменшення довжини та маси составу; при цьому змінюється координата його першої осі Sj та перераховується основний питомий опір W0. Це, в свою чергу, викликає відповідні зміни в режимі роботи гіркового локомотива, що спрямовані на підтримання заданої швидкості розпуску v0.
Рух составу в моделі описується за допомогою диференційного рівняння другого порядку S' = ft, S, S'), в якому незалежною змінною є час t:
S ” =
d 2 S dt2
g /р10-3
1 + у
(2)
де ---- - прискорення сили тяжіння з ураху-
1 + у
ванням інерції мас, що обертаються.
Використання рівняння руху (2) дозволяє виконувати сумісне моделювання процесів розпуску составу та скочування його відчепів.
14
Відомо, що рівняння (2) має єдиний розв'язок, якщо його права частинаf(S, V) неперервна та диференційована. Проте характер зміни сили f у даному виразі не завжди відповідає вказаній умові. Так, в моменти переключення позицій контролера може відбуватися стрибкоподібна зміна дотичної сили тяги /к. Гальмівна сила Ьт в процесі зменшення швидкості составу також змінюється нерівномірно. Опір руху в кривих являє собою ступінчасту функцію, розриви якої мають місце в точках заняття та звільнення маневровим составом криволінійних ділянок колії. Крім того, в момент відриву на ВГ чергового відчепу змінюються параметри составу. Тому в моделі було прийнято, що в межах кроку інтегрування At режим руху маневрового составу залишається постійним; з цією метою обраний достатньо малий крок At (At = 1 с). В режимі тяги в межах At позиція контролера не змінюється; при цьому в кінці кроку здійснюється аналіз швидкості составу і при необхідності корегується режим роботи локомотива та моделювання повторюється.
Аналогічно, гальмівна сила Ьт на кожному кроці At приймається постійною і, якщо швидкість составу в кінці кроку стала нижче допустимої, то, в залежності від умов руху, відбувається перехід до режиму тяги або вибігу.
Крім того, в межах кроку інтегрування At ні голова маневрового составу, ні його хвіст не повинні переходити границі початку або кінця кривої. Якщо на окремому кроці ця умова не виконується, то цей крок розділяється на окремі частини з відомою довжиною ASi, AS2, ... , ASn При цьому на кроках Н(п-1) моделювання переміщення составу виконується з використанням диференційного рівняння руху першого порядку V =f(S,V) з незалежною змінною S:
V =dV-S-А V > 0 (3)
dS 1 + у V
На останньому п-му кроці моделювання переміщення составу виконується за допомогою рівняння (2) на час At' = At - At*, що залишився до кінця початкового кроку At. Аналогічний алгоритм використовується і для моделювання руху составу при відриві від нього на ВГ чергового відчепу.
Для забезпечення неперервності функцій i(S) та f(V) в моделі використовувалась методика сплайн-апроксимації поздовжнього профілю колії [8] та тягових характеристик локомотива. Інтегрування диференційних рівнянь руху (2) та (3) виконувалось за методом Рунге-Кутта IV порядку [9].
Витрати палива маневровим тепловозом G в процесі насуву та розпуску составу визначаються як
a? n 60
G = G
j=i
де Gj - витрати палива на кроці At, кг/хв.
Величина витрат палива Gj в режимі тяги визначається за навантажувальними характеристиками в залежності від ефективної потужності Ne та частоти обертання колінчастого валу дизеля пд при роботі локомотива на окремій позиції контролера пк Відомо, що регулятори тепловозів забезпечують стале значення частоти пд для кожної позиції контролера [10].
Величина Ne може бути розрахована за наступною формулою [11]
ЛГ N к
Ne =- к
(1 - Р)п
де Р - частка потужності двигуна, що витрачається на привід допоміжного обладнання локомотива;
^ - к. к. д. тягової передачі;
N - дотична потужність локомотива, що визначається як [12]
Nк
Fv. 270
Значення частоти пд для окремої позиції контролера та величин Р, ^ приймаються згідно з паспортними характеристиками локомотива. Криві витрат палива G = f(Ne) для кожної позиції контролера машиніста побудовані на основі паспортних навантажувальних характеристик маневрового тепловоза [14] (рис. 1, а) та представлені відповідними регресійними моделями (рис. 1, б).
Витрати палива в режимі вибігу або гальмування залежать від частоти обертання колінчастого валу дизеля на холостому ході та приймаються як стала величина згідно з даними, наведеними в [4, табл. 29].
Таблиця 1
Режими роботи гіркового локомотива
Режим Положення органів керування Керовані параметри
п к п т F к B т
Тяга 1 2 min 0
8 2 max 0
Вибіг 0 2 0 0
Г альмування 0 3 0 min
0 6 0 max
15
Рис. 1. Витрати палива тепловозом ТЭМ2 в режимі тяги:
а) паспортна навантажувальна характеристика; б) регресійні моделі витрат палива
При моделюванні состав, що розформовується, розглядається як керована система, що функціонує в умовах дії зовнішніх та внутрішніх факторів, а також керуючих впливів [13]. Керований рух составу визначається режимом роботи гіркового локомотива. При цьому основними керованими параметрами є дотична сила тяги Fk та гальмівна сила Вт маневрового тепловоза, що залежать від встановленої позиції контролера пк та положення крана допоміжного гальма (ум. № 254) пт (див. табл. 1).
В процесі моделювання необхідно обрати таке керування гірковим локомотивом, щоб швидкість составу в момент завершення насуву їзн дорівнювала потрібній швидкості розпуску v0, а подальша фазова траєкторія V(t) для всіх tjH < t < tK відповідала встановленому режиму розформування, де tK - кінцевий момент руху составу. Початковий момент часу t0 приймається рівним 0; в якості кінцевого моменту tK розглядається момент відриву від локомотива останнього відчепа.
В роботі було розроблено алгоритм керування гірковим тепловозом, який, окрім вимог з безпечного виконання маневрової роботи та експлуатації локомотивів, враховує й біхевіо-ральні фактори, що пов’язані з керуючими діями машиніста. Даний алгоритм забезпечує плавний розгін та подальший рух составу з швидкістю, близькою до встановленої v,j. При цьому фактична швидкість уц0п на кожному кроці At може відхилятись від заданої швидкості руху Vj(t) на величину похибки реалізації
б: Удоп Є [vmin; vmax], vmm = Vj - б, vmax = v + б.
В початковий момент часу t0 для зрушення составу з місця встановлюється режим тяги, при цьому контролер машиніста переміщується на першу позицію (пк = 1). Подальший набір позицій виконується з мінімально допустимою затримкою контролера в кожному положенні
(3 с). У випадку, коли состав не вдається зрушити з місця протягом 10 с [14], то констатується необхідність попереднього його осаджування, та моделювання насуву припиняється.
При керуванні локомотивом машиніст, по-кладаючись на власний досвід, прогнозує характер руху маневрового составу та заздалегідь вживає заходи щодо дотримання встановленого режиму насуву. Так, при наближенні швидкості руху до встановленої швидкості розпуску машиніст забезпечує плавний вихід на неї, заздалегідь зменшуючи силу тяги локомотива. Аналогічно здійснюються й зниження швидкості руху, при якому машиніст, беручи до уваги інерційність систем локомотива та оцінюючи запас кінетичної енергії составу, почергово застосовує режими гальмування та вибігу.
Враховуючи вказані обставини в моделі було встановлено зону регулювання швидкості [v№n - 0,5; удоп], при досягненні якої процес розгону з максимальним темпом переводу контролера завершується, і модель переходить в режим прогнозного керування рухом составу.
В даному режимі на кожному кроці моделювання At, в залежності від величини питомого опору руху wK та режиму роботи локомотива на попередньому кроці, виконується пошук такої комбінації керуючих впливів (пк, пт), яка б на протязі періоду прогнозування ^р забезпечувала допустимий режим руху составу. При цьому слід зауважити, що робота маневрового тепловоза в режимі постійних розгонів-сповільнень призводить до підвищення витрат палива [15] та негативних динамічних реакцій состава [16]. У зв’язку з цим необхідно обмежити частоту зміни режиму роботи гіркового локомотива, для чого в моделі величина періоду tnp приймається рівною ^р = 2At = 2 с. Нижче наведений алгоритм прогнозного керування
16
рухом составу для різних режимів роботи локомотива.
В режимі тяги виконується пошук максимально можливої позиції контролера машиніста пк, при якій би в межах періоду їпр забезпечувалась допустима швидкість руху удоп. Причому, перевід контролера в бік скидання навантаження дозволяється виконувати одразу на декілька позицій та до завершення мінімально допустимої витримки часу (3 с). Якщо така позиція не знайдена, то область пошуку допустимого керування гірковим локомотивом на даному кроці переміщується на режим вибігу.
В режимі вибігу за допомогою (2) виконується розрахунок швидкості v;+2, з якою рухатиметься состав в кінці поточного періоду прогнозування ^р. Якщо швидкість v;+2 перевищить максимально допустиму (v;+2 > vmax), то для моделювання руху составу в межах даного періоду tnp встановлюється режим гальмування.
У випадку, коли в режимі вибігу прогнозована швидкість Vj+2 опуститься нижче мінімально допустимої (v+2 < vmin), необхідно перевірити, чи не виконувався на поточному кроці tj пошук допустимого керування в режимі тяги. Якщо такий пошук мав місце, то з міркувань безпеки виконання маневрової роботи для мо-
оду tnj, встановлюється режим вибігу. В протилежному випадку відбувається перехід до режиму тяги, причому набір навантаження обов’язково виконується починаючи з першої позиції (пк = 1).
Аналогічно, в режимі гальмування визначається швидкість Vj+2 в кінці поточного періоду прогнозування ^р. Якщо дана швидкість виявиться нижчою за мінімально допустиму (Vj+2 < vmln), а состав знаходиться на прискорюючому ухилі (»к < 0), то необхідно перевірити, чи не виконувався на поточному кроці tj пошук допустимого керування в режимі вибігу. У випадку, коли такий пошук мав місце, для моделювання руху в межах даного періоду ^р встановлюється режим гальмування. В протилежному випадку виконується перехід до режиму вибігу. Інакше, коли v;+2 < vmin та w > 0, для моделювання приймається режим тяги (пк = 1).
Розроблена модель дозволяє імітувати різноманітні режими розформування составів. Для прикладу на рис. 2 наведено результати моделювання насуву та розпуску составу масою 3869 т гірковим локомотивом ТЭМ2. При цьому розглядався режим, в якому розгін составу здійснювався до встановленої швидкості розпуску v0 = 1,7 м/с, після чого його розформування
делювання руху составу в межах даного пері- виконувалось з постійною швидкістю.
а)
vo,
м/с
1,6
1,2
0,8
0,4
0
A B C
1,8 - Перехід А 1,8 - Ділянка B-C
1,7 -
1,5 - / 1,5 -
7
820 825 830 835 840 845 850 2100 2120 2140 2160 2180
800 1000 1200 1400 1600 1800 2000 S', м
Sj ,м V, м/c t, с Реж Им,ткм Ис,ткм Ит,ткм E/G
804.0 0.0000 О о П1 0.00000 0.00000 0.00000 0.00000
804.0 0.0000 1.0 П1 0.00000 0.00000 0.00000 0.00333
804.0 0.0000 2.0 П1 0.00000 0.00000 0.00000 0.00667
804.0 0.0000 о СП П1 0.00000 0.00000 0.00000 0.01000
804.0 0.0055 4.0 П2 0.00003 0.00001 0.00000 0.01317
804.0 0.0108 5.0 П2 0.00011 0.00004 0.00000 0.01633
1407.33 1.6970 376.0 П3 5.39102 2.40666 0.00000 1.11804
Состав после отрыва . отцепа . 1: sj=1392.71 L = 7 8 9.38 Q= 3827.00 os =224
1394.40 1.6926 377.0 П3 5.40268 2.41345 0.00000 1.12111
1405.70 1.5448 832.0 П1 9.29912 4.04906 0.00000 2.53835
Состав после отрыва отцепа 35: sj=1391.78 L = 16.4 5 Q=0 .00 os=0
Рис. 2. Результати моделювання насуву та розпуску составу:
а) графік швидкості составу; б) протокол роботи моделі
17
Як видно з рис. 2, а, розроблений алгоритм моделювання насуву та розпуску забезпечує плавний розгін составу та подальший його рух зі швидкістю Кф, близькою до встановленої швидкості розпуску v0 (Кф є [1,52; 1,74]). Слід зауважити, що на заключному етапі розформування (рис 2, а, ділянка В-С), коли довжина составу не перевищує 10 вагонів, режим його руху різко змінюється. Ці зміни спричинені зменшенням маси составу в кінці розпуску, що викликає завищену інтенсивність розгону навіть на перших позиціях контролера.
Розроблена модель насуву та розпуску дозволяє детально імітувати режим розформування составу на сортувальній гірці та визначати пов'язані з цим витрати палива гірковим локомотивом. При цьому структура побудованої моделі передбачає її інтеграцію в існуючу модель скочування відчепів [17] для створення об’єднаної моделі розформування составів на основі системного підходу. Така модель дозволить комплексно оцінювати якість сортувального процесу та може бути використана в системі підтримки прийняття рішень щодо визначення ефективних режимів функціонування сортувальних комплексів станцій в умовах змінної інтенсивності вхідного потоку поїздів.
БІБЛІОГРАФІЧНИЙ СПИСОК
1. Бобровский, В. И. Совершенствование имитационной модели процесса надвига и роспуска составов на сортировочных горках [Текст] /
B. И. Бобровский, Е. Б. Демченко // Зб. наук. праць ДНУЗТ «Транспортні системи та технології перевезень».- Д.: Вид-во ДНУЗТ, 2012. - Вип. 3. - С. 5-9.
2. Огар, О. М. Науковий підхід до визначення раціональних конструктивно-технологічних параметрів сортувальних гірок [Текст] / О. М. Огар // Зб. наук. праць ДонІЗТ. - 2009. - № 18. - С. 9-16.
3. Овчинников, В. М. Сокращение расхода дизельного топлива на маневрах [Текст] / В. М. Овчинников, С. А. Пожидаев, Н. Г. Швец, В. В. Скре-жендевский // Зб. наук. праць ДНУЗТ «Транспортні системи та технології перевезень».- Д.: Вид-во ДНУЗТ, 2011. - Вип. 1. - С. 62-70.
4. Правила тяговых расчетов для поездной работы [Текст]. - М.: Транспорт, 1985. - 287 с.
5. Назаров, Л. С. Повышение эффективности маневровой работы [Текст] / Л. С. Назаров, С. Л. Назаров // Железнодорожный транспорт. - 2001. - №8. -
C. 56-57.
6. Правила и нормы проектирования сортировочных устройств на железных дорогах Союза ССР. ВСН 207-89. - М.: Транспорт, 1992. - 104 с.
7. Тяговые расчеты [Текст]: справочник / Под ред. П. Т. Гребенюка. - М.: Транспорт, 1987. - 272 с.
8. Бобровский, В. И. Представление продольного профиля сортировочных горок в АСУ расформированием составов // Информационно-управляющие системы на железнодорожном транспорте. - 1996. -№ 1, 2. - С. 19-25.
9. Бобровский, В. И. Дифференциальные уравнения движения отцепа и методы их решения // Ин-формационно-управляющие системы на железнодорожном транспорте. - 1996. - № 6. - С. 34-39.
10. Теорія та конструкція локомотивів. Основи проектування [Текст] : підручник для ВНЗ залізнич. трансп. / Б. Є. Боднар, Є. Г. Нєчаєв, Д. В. Бобирь / Під ред. Б. Є. Боднара. - Д.: ПП «Ліра ЛТД», 2010. -358 с.
11. Теория локомотивной тяги [Текст] : учебник для вузов ж.-д. транспорта / В. Д. Кузьмич, В. С. Руднев, С. Я. Френкель / Под ред. В. Д. Кузьмича. -М.: Маршрут, 2005. - 448 с.
12. Справочник по электроподвижному составу, тепловозам и дизель-поездам [Текст] / Под ред.
А. И. Тищенко. - М.: Транспорт, 1976. - Т. I. - 432 с.
13. Бобровский, В. И. Имитационная модель развязки линий в железнодорожном узле // Концепція підвищення ефективності вантажних перевезень на залізничному транспорті: Міжвуз. зб. наук. праць -Вип. 38. - Харків: ХарДАЗТ, 1999. - С. 35-42.
14. Тепловозы ТЭМ1 и ТЭМ2 [Текст]: Изд. 2-е, испр. и доп. / Под ред. Е. Ф. Сдобникова. - М.: Транспорт, 1978. - 279 с.
15. Механизация и автоматизация формирования поездов [Текст] / Ю. А. Муха, В. А. Король, Н. М. Иванков и др. / Под общ. ред. Ю. А. Мухи. - К.: Техніка, 1987. - 136 с.
16. Бурченков, В. В. Исследование ускорений, развиваемых локомотивами в процессе роспуска состава [Текст] / В. В. Бурченков // Теоретические и экспериментальные исследования систем автоматизации процесса расформирования составов на сортировочных горках: Труды ВНИИЖТ. - М.: Транспорт, 1983. - вып. 661. - С. 81-91.
17. Оптимизация режимов торможения отцепов на сортировочных горках [Текст] : Монография / под. ред. В. И. Бобровского - Д.: Изд-во Маковецкий, 2010. - 260 с.
Надійшла до редколегії 27.11.2012.
Прийнята до друку 28.11.2012.
18