База имитационных модолей СППР-УДП (стр. 1 из 3)
Лекция 9. БАЗА ИМИТАЦИОННЫХ МОДЕЛЕЙ СППР-УДП
9.1. Математические модели движения поездов на участке железной дороги
Предложенная в [8, 123] математическая схема положена в основу разработки базы имитационных моделей (БИМ) СППР для всех классов объектов (P, М, Т, А) и типов элементов, входящих в функциональные схемы СППР. Это дало возможность рассмотреть вопросы описания структуры технологического комплекса управления движением поездов, позволяющих моделировать функциональные связи между отдельными элементами исследуемой системы.
Так как основной целевой функцией системы является управление движением поездов, то разработка имитационных моделей отдельных элементов технологического комплекса начинается с рассмотрения вопросов моделирования процессов движения (модели класса М).
Для управления процессом движения поездов уже разработаны различные варианты построения математических моделей движения [8,75-80] на участках железных дорог. Применив вышепредложенную математическую схему, основные результаты исследований в рассматриваемой области можно обобщить следующими основными положениями.
1. При моделировании процессов движения на участке длиной L в качестве обслуживающего прибора рассматривается этот участок. В общем случае он характеризуется набором параметров, существенных для расчетов времени его занятия и освобождения. К таким параметрам относятся характеристики плана и профиля участка, величина постоянных и временных ограничений скорости и др.
2. При моделировании процессов движения в качестве заявок рассматриваются поезда, характеризующиеся множествами статистических СЗi и динамических dЗi параметров. Множество СЗi в общем случае включает такие параметры, как сила тяги Fi и торможения fтi состава, вес состава Рi и другие характеристики подвижного состава.
В качестве основных динамических параметров dЗi принимаются скорость Vi и время движения поезда Dti по дискретным отрезкам пути (дискрет-участкам) длиной DI. Обслуживающий прибор DI рассматривается в виде многоканального устройства с динамическим числом каналов обслуживания Уi, которое определяется выражением
где lпi – длина i-го поезда.
В процессе движения поезда по дискрет-участкам Dl (рис. 9.1) складывается следующая технологическая ситуация по времени: (t1, t3, t2 и t4), где
t1 обозначает время занятия поездом данного участка Dii; t2 – освобождение предыдущего Dii-1; t3 – занятие следующего Dii+1; t4 – освобождение занятого участка Dii.
Рис. 9.1. Схема продвижения состава модели
Значения скоростей движения поезда в моменты времени t1, t3, t2, t4
соотносятся в зависимости от режима движения следующим образом:
(V1 = V3 = V2 = V4) – при равномерном движении поезда;
(V1 < V3 < V2 < V4) – в режиме ускорения движения поезда;
(V1 > V3 > V2 > V4) – в режиме торможения поезда.
Предложенный подход позволяет свести моделирование непрерывного процесса движения к расчету значений скоростей в заданных точках пути, где эти значения существенны для моделирования работы всего технологического комплекса.
3. Функция обслуживания заявки Зi по динамическому параметру f i (определение скоростей v1 – v4) рассчитывается по заданному уравнению движения V = f(S) в двух возможных режимах – без прерывания (уравнение в пределах участка не изменяется) и с прерыванием (уравнение движения изменяется). Методы расчета, оптимизации и аппроксимации уравнений движения исследованы в работах [11, 73-77, 83-116].
4. Функция обслуживания заявки Зi по времени fi определяется по вычисленным значениям скоростей v1 – v4 с учетом принятой аппроксимации уравнения движения в пределах участка железной дороги L.
При оценке точности аппроксимации установлено, что во всех случаях достаточная точность моделирования достигается при аппроксимации по средней скорости движения с расчетом времени занятия обслуживающего прибора p i длиной Dl по следующим уравнениям:
,где v1 – средняя скорость движения первой оси; t1 – время движения первой оси; v2 – средняя скорость движения последней оси; t2 – время движения последней оси; t – время занятия обслуживающего прибора pi длиной Dl.
В случае прерывания процесса обслуживания заявки во время движения первой оси (момент t', t1
t' tЗ) время занятия 
.При прерывании процесса обслуживания заявки во время выхода последней оси (момент t", t2
t" t4) время занятия Dl будет равно 
Приведенная математическая модель обеспечивает довольно малое время счета, но точность вычислений определяется принятой длиной обслуживающего прибора Dl и адекватностью зависимостей V = f(S) реальному процессу движения поездов, которые аппроксимируются. При этом эти зависимости не учитывают взаимодействие соседних поездов при их движении.
9.2. моделирование движения поездов с применением тяговых расчетов
При имитационном моделировании уравнение движения определяется блоком выбора машинистом режима движения и решается для каждого поезда отдельно с учетом координат и скоростей впереди и сзади следующих поездов. В общем виде решаемое уравнение [63] выглядит следующим образом:
(9.1)где (fК – fТ – w0 – wВ + iК) – удельные силы, действующие на поезд при его движении (при этом некоторые параметры могут иметь нулевое значение в зависимости от условий движения); fК – удельная сила тяги поезда; fТ – удельная сила торможения; w0 – удельное основное сопротивление движению; ± iК – дополнительное сопротивление движению от уклона и кривизны пути; wВ – удельная сила сопротивления от воздушной среды; e – коэффициент ускорения поезда.
Представим уравнение (9.1) в виде, удобном для интегрирования:
проинтегрировав его, получим:
или 
, (9.2)где wк = fт + w0 + wв + iк.
Умножив левую и правую часть уравнения (9.2) на скорость, определенную в любой точке пути (V = f(S)), получим:
.Проинтегрировав это выражение, получим:
или 
.Обозначив
, имеем 
.Так как
, можно записать:
(9.3)Если принять, что на интервале скорости (V1 – V2) сила, действующая на поезд (fк – wк), постоянна, то из уравнения (9.1) получим:
откуда
или 
, (9.4)или
.Из (9.3) при постоянной силе (fк – wк) получим:
,или
.Зависимость t = f(S) получается путем исключения V из найденных зависимостей V = f(t) и V = f(S). Из (9.4) имеем
. (9.5)Подставляя это выражение в (9.5), находим: