Исследование динамики потока в модели организации грузоперевозок по круговой цепочке станций
Исследование динамики потока в модели организации грузоперевозок по круговой цепочке станций
Аннотация
Код статьи
S042473880013024-5-1
Тип публикации
Статья
Статус публикации
Опубликовано
Авторы
Хачатрян Нерсес Карленович 
Должность: Заместитель директора по научной работе, ЦЭМИ РАН
Аффилиация: Центральный экономико-математический институт РАН
Адрес: Москва, РФ
Бекларян Левон Андреевич
Аффилиация: ЦЭМИ РАН
Адрес: Российская Федерация, Москва
Выпуск
Страницы
83-91
Аннотация

Статья посвящена построению и исследованию модели организации железнодорожных грузоперевозок по круговой цепочке станций. Каждая станция характеризуется определенным числом путей, каждый из которых в произвольный момент времени может быть задействован, а также эффективностью их использования.  Движение грузопотока осуществляется с помощью двух технологий. Первая технология основана на взаимодействии соседних станций и задает интенсивность потока между ними в зависимости от соотношения свободных путей на них и нормативного коэффициента, характеризующего пропускную способность перегонов (участков железнодорожной линии между станциями) и технические характеристики железнодорожного подвижного состава, осуществляющего перевозки.  Ее дополняет вторая технология, задача которой заключается в том, чтобы в полной мере использовать пропускную способность станций, которая выражается количеством перевозимого груза за единицу времени в зависимости от загруженности станций и эффективности использования путей. Такая модель описывается системой дифференциальных уравнений с разрывными правыми частями. Аналитическое решение такой системы крайне затруднительно, поэтому, она была исследована численно. Приведены результаты численного исследования указанной системы, основная цель которых является определение динамики грузопотока, а также изучение ее зависимости от параметров модели, характеризующих пропускную способность станций и перегонов, а также от загруженности станций в начальный момент времени.

Ключевые слова
организация железнодорожных грузоперевозок, интенсивность грузопотока, дифференциальные уравнения, численное решение, стационарный поток.
Источник финансирования
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований в рамках проектов №19-01-00147 и №19-010-00958.
Классификатор
Получено
25.03.2021
Дата публикации
29.03.2021
Всего подписок
24
Всего просмотров
1421
Оценка читателей
0.0 (0 голосов)
Цитировать   Скачать pdf
1

ВВЕДЕНИЕ

2 Транспорт играет важную роль в экономике страны. Он обеспечивает развитие, связь и координацию работы всех отраслей экономики, а также связь между регионами страны. Развитие транспорта является важнейшим условием модернизации экономики. Оно невозможно без оптимального планирования сетей, улучшения организации движения и ряда других задач, решение которых невозможно без математического моделирования транспортных потоков. Здесь можно выделить два крупных направления: моделирование загрузки транспортных сетей городов и моделирование динамики транспортного потока (Швецов, 2003). Первое направление представлено моделями расчета корреспонденций такими как гравитационная модель (Carrothers G. A. P. 1956; Wilson A. G., 1971), энтропийная модель (Harris B., Wilson A. G., 1978; Popkov Yu. S., 1995), модели семейства конкурирующих центров (Fotheringham A. S., 1983; Fotheringham A. S., 1986), а также моделями распределения потоков по сети, включая различные варианты равновесного распределения и алгоритм оптимальных стратегий (Shvetsov V.I., 2009; Leventhal T., Nemhauser G.L., Trotter L., 1973; Lo H.K., Chen A., 2000; Bar-Gera H., 2002; Spiess H., Florian M., 1989). Второе направление представлено основными классами динамических моделей – макроскопическими (гидродинамическими), кинетическими (газодинамическими) и микроскопическими. Макроскопические модели (Daganzo C. F. 1994; Daganzo C. F. 1995; Гасников и др., 2013; Сухинова и др., 2009; Иносэ, Хамада, 1983) описывают движение транспортных средств в усредненных терминах, таких как плотность, средняя скорость, поток и др. При таком подходе транспортный поток уподобляется движению специфической жидкости, поэтому модели этого класса также называют гидродинамическими. Микроскопическими называются модели, в которых явно моделируется движение каждого автомобиля. Такой подход позволяет теоретически достичь более точного описания движения транспортных средств по сравнению с усредненным макроописанием, однако этот подход требует больших вычислительных ресурсов при практических применениях. Примерами таких моделей являются модели следования за лидером (Brackstone M., McDonald M., 2000), модели оптимальной скорости (Bando M., Hasebe K., Nakayama A., Shibata A., Sugiyama Y., 1995), модель Трайбера (Treiber M., Hennecke A., Helbing D., 2000), а также модели клеточных автоматов (Cremer M., Ludwig J., 1986; Chowdhury D., Santen L., Schadschneider A., 2000). Промежуточное место занимает кинетический подход, при котором поток описывается плотностью распределения транспортных средств в фазовом пространстве, т.е. пространстве координат и скоростей автомобилей (Helbing D., Treiber M., 1998; Nelson P., 1995). Динамика фазовой плотности описывается кинетическим уравнением. Это уравнение основано на усредненном описании эффектов взаимодействий индивидуальных транспортных средств, и в этом смысле оно ближе к микроуровню, чем гидродинамические уравнения. Отметим, что модели, в которых движение транспортного средства уподобляется какому-либо физическому потоку (гидро - и газодинамические модели), как правило применяются при описании движения автомобильного транспорта.
3 В данной статье представлена макроскопическая динамическая модель, описывающая движение железнодорожного транспорта по сети, представляющей собой круговую цепочку станций (Beklaryan L.A., Khachatryan N.K., 2006; Khachatryan N.K., Akopov A.S., 2017 ; Khachatryan N.K., Akopov A.S., Belousov F.A., 2018; Beklaryan Levon A., Khachatryan Nerses K., Akopov Andranik S., 2019; Бекларян Л.А., Хачатрян Н.К., 2013; Бекларян Л.А. , Хачатрян Н.К. , 2019). Транспортный поток в железнодорожной сети формируется в результате применения технологий, осуществляющих взаимодействие соседних станций с учетом их пропускной способности, а также пропускной способности перегонов. Такая модель позволяет прогнозировать загруженность станций и динамику потока, возникающего в железнодорожной сети. В отличие от аналогичных моделей, описанных в работах, указанных выше, станции имеют различные характеристики (количество путей и эффективность их использования).
4

1. ПОСТАНОВКА ЗАДАЧИ

5 Рассмотрим движение грузопотока по круговой цепочке, состоящей из n станций. Обозначим через Δi,    i=1,  2,...,n число путей на i-й станции. Каждая станция в произвольный момент времени характеризуется количеством задействованных путей. Пусть zi(t),    i=1,  2,...,n число путей, задействованных на i-й станции в момент времени t. Тогда Δi-zi(t),    i=1,  2,...,n - число свободных путей на i-й станции в момент времени t. Очевидно, что функции zi(t) должны удовлетворять следующим ограничениям:
6 0zi(t)Δi,    i=1,  2,...,n . (1)
7 Движение грузопотока осуществляется с помощью двух технологий. Первая технология основана на процедуре взаимодействия соседних станций. Согласно ей произвольная станция
8
  1. принимает груз с предыдущей станции, если число свободных путей на ней больше чем на предыдущей станции и предыдущая станция не пуста. При этом интенсивность приема пропорциональна как разности чисел свободных путей на этих станциях, так и числу свободных путей на данной станции.
  2. Отправляет груз на следующую станцию, если число свободных путей на ней меньше чем на следующей станции, и она не пуста. При этом интенсивность отправки пропорциональна как разности чисел свободных путей на этих станциях, так и числу свободных путей на следующей станции.
9 Таким образом, каждая станция с номером i ( 2in-1 ) в рамках первой технологии должна принять груз с предыдущей станции с интенсивностью, равной α[ (Δi-zi)-(Δi-1-zi-1)](Δi-zi)sign(Δi-zi)sign(zi-1) , если Δi-zi>Δi-1-zi-1 и отправлять груз на следующую станцию с интенсивностью, равной α[ (Δi+1-zi+1)-(Δi-zi)](Δi+1-zi+1)sign(Δi+1-zi+1)sign(zi) , если Δi+1-zi+1>Δi-zi .
10 Для 1 -й станции предыдущей является n -я станция, соответственно, для n -й станции последующей является 1 -я станция. Поэтому, 1 -я станция в рамках первой технологии должна принять груз с предыдущей станции с интенсивностью, равной α[ (Δ1-z1)-(Δn-zn)](Δ1-z1)sign(Δ1-z1)sign(zn) , если Δ1-z1>Δn-zn и отправлять груз на следующую станцию с интенсивностью, равной α[ (Δ2-z2)-(Δ1-z1)](Δ2-z2)sign(Δ2-z2)sign(z1) , если Δ2-z2>Δ1-z1 .
11 Аналогично, n -я станция должна принять груз с предыдущей станции с интенсивностью, равной α[ (Δn-zn)-(Δn-1-zn-1)](Δn-zn)sign(Δn-zn)sign(zn-1) , если Δn-zn>Δn-1-zn-1 и отправлять груз на следующую станцию с интенсивностью, равной α[ (Δ1-z1)-(Δn-zn)](Δ1-z1)sign(Δ1-z1)sign(zn) , если Δ1-z1>Δn-zn .
12 Функция sign(.) определяется следующим образом:
13 sign(t)=1,     t>00,    t0.
14 Параметр α характеризует потенциал первой технологии и ограничивается пропускной способностью перегонов (участков железнодорожной линии между станциями) и техническими характеристиками железнодорожного подвижного состава, осуществляющего перевозки.
15 Каждая станция обладает определенной пропускной способностью. В данной модели она выражается количеством перевозимого груза за единицу времени в зависимости от загруженности станций, причем в общем случае эта зависимость у каждой станции своя из-за разницы в числе путей и эффективности их использования. Предполагается, что станции стремятся использовать пропускную способность в полной мере. Данная задача решается в рамках второй технологии, согласно которой произвольная станция дополнительно принимают и отправляют грузопоток с интенсивностью, определяемой пропускной способностью, соответственно, на ней и на следующей станции. Она описывается с помощью непрерывных убывающих функций φi(.),    i=1,2,...,n , определенных на отрезке [0, Δi] , причем φi(Δi)=0,    i=1,2,...,n .
16 Таким образом, интенсивность движения грузопотока на станциях задается с помощью следующей системы дифференциальных уравнений:
17 z˙1(t)=α (Δ1-z1-Δn+zn)sign(Δ1-z1-Δn+zn)(Δ1-z1)sign(Δ1-z1)sign(zn)- -α (Δ2-z2-Δ1+z1)sign(Δ2-z2-Δ1+z1)(Δ2-z2)sign(Δ2-z2)sign(z1)+ +φ1(z1)-  φ2(z2),      t[0, +) ; (2)
18 z˙i(t)=α (Δi-zi-Δi-1+zi-1)sign(Δi-zi-Δi-1+zi-1)(Δi-zi)sign(Δi-zi)sign(zi-1)- -α (Δi+1-zi+1-Δi+zi)sign(Δi+1-zi+1-Δi+zi)(Δi+1-zi+1)sign(Δi+1-zi+1)sign(zi)+ +φi(zi)-  φi+1(zi+1),      i=2,  3,...,n-1,    t[0, +) ; (3)
19 z˙n(t)=α (Δn-zn-Δn-1+zn-1)sign(Δn-zn-Δn-1+zn-1)(Δn-zn)sign(Δn-zn)sign(zn-1)- -α (Δ1-z1-Δn+zn)sign(Δ1-z1-Δn+zn)(Δ1-z1)sign(Δ1-z1)sign(zn)+ +φn(zn)-  φ1(z1),      t[0, +). (4)
20 Отметим, что конструкция правых частей системы (2)-(4) гарантирует, что всякое ее решение, удовлетворяющее условию (1) в начальный момент времени, будет удовлетворять ему и в последующие моменты времени.
21 Перепишем систему (2)-(4) в следующем виде:
22 z˙1(t)=α (zn-z1+Δ1-Δn)sign(zn-z1+Δ1-Δn)(Δ1-z1)sign(Δ1-z1)sign(zn)- -α (z1-z2+Δ2-Δ1)sign(z1-z2+Δ2-Δ1)(Δ2-z2)sign(Δ2-z2)sign(z1)+ +φ1(z1)-  φ2(z2),      t[0, +) ; (5)
23 z˙i(t)=α (zi-1-zi+Δi-Δi-1)sign(zi-1-zi+Δi-Δi-1)(Δi-zi)sign(Δi-zi)sign(zi-1)- -α (zi-zi+1+Δi+1-Δi)sign(zi-zi+1+Δi+1-Δi)(Δi+1-zi+1)sign(Δi+1-zi+1)sign(zi)+ +φi(zi)-  φi+1(zi+1),      i=2,  3,...,n-1,    t[0, +) ; (6)
24 z˙n(t)=α (zn-1-zn+Δn-Δn-1)sign(zn-1-zn+Δn-Δn-1)(Δn-zn)sign(Δn-zn)sign(zn-1)- -α (zn-z1+Δ1-Δn)sign(zn-z1+Δ1-Δn)(Δ1-z1)sign(Δ1-z1)sign(zn)+ +φn(zn)-  φ1(z1),      t[0, +). (7)
25 Помимо системы (5)-(7) рассмотрим ее частный случай соответствующий идентичности всех станций. т.е. в предположении, что все станции содержат одинаковое число путей и с одинаковой эффективностью их могут использовать. В этом случае Δi=Δj=Δ,   φi(.)=φj(.)=φ(.) для всех i,j=1,2,...,n и система уравнений (5)-(7) приобретает вид:
26 z˙1(t)=α (zn-z1)sign(zn-z1)(Δ-z1)sign(Δ-z1)sign(zn)- -α (z1-z2)sign(z1-z2)(Δ-z2)sign(Δ-z2)sign(z1)+ +φ(z1)-  φ(z2),      t[0, +) ; (8)
27 z˙i(t)=α (zi-1-zi)sign(zi-1-zi)(Δ-zi)sign(Δ-zi)sign(zi-1)- -α (zi-zi+1)sign(zi-zi+1)(Δ-zi+1)sign(Δ-zi+1)sign(zi)+ +φ(zi)-  φ(zi+1),      i=2,  3,...,n-1,    t[0, +) ; (9)
28 z˙n(t)=α (zn-1-zn)sign(zn-1-zn)(Δ-zn)sign(Δ-zn)sign(zn-1)- -α (zn-z1)sign(zn-z1)(Δ-z1)sign(Δ-z1)sign(zn)+ +φ(zn)-  φ(z1),      t[0, +). (10)
29 Несложно заметить, что система уравнений (8)-(10) имеет бесконечное множество стационарных решений вида:
30 zi=c,    0cΔ,    i=1,  ...,n . (11)
31 Аналитическое исследование других решений системы (8)-(10), а также решений системы (5)-(7), крайне затруднительно так как правые части дифференциальных уравнений в общем случае являются разрывными функциями. В связи с этим, системы (8)-(10) и (5)-(7) были исследованы численно.
32

2. РЕЗУЛЬТАТЫ ЧИСЛЕННЫХ ЭКСПЕРИМЕНТОВ

33 Перейдем к изложению результатов численных экспериментов, в которых функции φi(.) были определены следующим образом:
34 φi(zi)=ai(Δi2-zi2),    i=1,2,...,n . (12)
35 Напомним, что функции φi(.) должны быть убывающими и удовлетворяющими условию φi(Δi)=0 . Параметр ai , участвующий в определении этих функций описывает эффективность использования путей на i -ой станции. Также отметим, что определяя такой класс функций φi(.) , мы предполагаем, что пропускная способность станций зависит как от числа путей, так и от эффективности их использования, причем с ростом эффективности использования путей растет линейно, а с увеличением числа путей больше чем линейно. Основная задача численных экспериментов заключалась в исследовании динамики решений системы (5)-(7) (в частности (8)-(10)), определяющей динамику интенсивности грузопотока, а также ее зависимость от начальных значений и параметров модели α,  ai,  Δi,    i=1,2,...,n .
36 Начнем с численного исследования системы (8)-(10). Отметим, что для системы (8)-(10) φi(zi)=φ(zi)=a(Δ2-zi2),    i=1,2,...,n , т.е. параметрами являются α,  a,  Δ .
37 2.1. Численное решение системы (8)-(10)
38 Как было отмечено выше, система (8)-(10) имеет стационарные решения вида (11). Как показывают численные эксперименты при любых значениях параметров α,  a,  Δ система (8)-(10) имеет непрерывное решение, причем любое ее решение со временем выходит на стационарный режим, т.е. существует t->0 такое, что всякое решение {zi()}1n системы (8)-(10) удовлетворяет условию zi=c-,   t[t-, +) . Более того значение с- зависит только от начальных условий и определяется следующим образом:
39 c-=1ni=1nzi(0) .
40 Момент времени t- зависит как от начальных условий так и от параметров модели. На рис.1 приведен график решений системы (8)-(10) при α=2,  a=0.1,  Δ=10   и следующих начальных условиях: z1(0)=3,  z2(0)=2,  z3(0)=5,  z4(0)=3,  z5(0)=4,  z6(0)=4,  z7(0)=3,  z8(0)=2,  z9(0)=1,  z10(0)=4 .
41

42 Рис.1. График решения системы (8)-(10)
43 Из правых частей уравнений (8)-(10) следует что начиная с момента времени t- первые два слагаемых в них равны нулю, ненулевыми являются только слагаемые определяемые функцией φ(.) . Это означает, что со временем движение грузопотока осуществляется только в рамках второй технологии, каждая станция принимает груз с предыдущей станции и отправляет на следующую станцию с постоянной интенсивностью φ(c-) , т.е. в системе грузоперевозок устанавливается стационарный поток На рис.2. приведены графики интенсивностей потоков, возникающих между станциями.
44

Рис.2 Интенсивность потока в случае идентичных станций

45 Как следует из определения функции φ(.) величина этого потока прямо пропорциональна параметрам a,  Δ и обратно пропорциональна значению c- . Таким образом, в случае идентичности станций в системе грузоперевозок начиная с некоторого момента времени устанавливается стационарный поток, величина которого увеличивается по мере роста пропускной способности станций и уменьшается с увеличением загрузки станций в начальный момент времени ( c- ). Отметим, что величина этого потока не зависит от характеристики первой технологии, т.е. параметра α . От него зависит время выхода на стационарный поток ( t- ). Чем меньше значение параметра α , тем позже наступит это время. Наконец отметим, что все выводы полученные в этом параграфе относительно динамики решений системы (8)-(10) и интенсивности потока остаются справедливыми и при α=0 , т.е. в системе грузоперевозок устанавливается такой же стационарный поток как и при α>0 , но в более отдаленном будущем.
46 2.2. Численное решение системы (5)-(7)
47 Перейдем к исследованию общей системы (5)-(7). В результате численных экспериментов было установлено, что система (5)-(7) имеет два типа решений. Первый тип решений характеризуется непрерывными правыми частями. Существуют t~>0 такое, что zi=c-i,  0<c-i<Δi,   i=1,  2,...,n,    t[t~, +) . (13)
48 Кроме того, выполнятся условие
49 1ni=1nc-i=c- , где c-=1ni=1nzi(0) .
50 Второй тип решения системы (5)-(7) характеризуется разрывными правыми частями. Как показывают численные эксперименты такие решения являются следствием того, что одна или несколько компонент решения начинают колебаться в правой окрестности нуля, принимая периодически нулевые значения (терпят разрывы). Это приводит к разрывам остальных компонент решения. Вследствие такого поведения решений интенсивность грузопотока меняется скачкообразно. Отметим, что такое поведение характерно для решений с малыми начальными значениями и решениями, полученными при значениях параметра α , близкого к нулю. Также отметим, что при α=0 решения системы (5)-(7) являются разрывными независимо от начальных условий. На практике такой режим грузоперевозок крайне сложно реализуем, поэтому, такой тип решений мы в дальнейшем рассматривать не будем.
51 Итак, наши дальнейшие исследования связаны с первым типом решений. На рис.3 приведен пример такого типа решения при α=2 и следующих значениях остальных параметров и начальных условий:
52 Δ1=10,  Δ2=11,  Δ3=12,  Δ4=8,  Δ5=10,  Δ6=11,  Δ7=9,  Δ8=10,  Δ9=12,  Δ10=10 . a1=0.1,  a2=0.2,  a3=0.3,  a4=0.05,  a5=0.4,  a6=0.1,  a7=0.1,  a8=0.2,  a9=0.2,  a10=0.02 . (14) z1(0)=4,  z2(0)=2,  z3(0)=5,  z4(0)=7,  z5(0)=6,  z6(0)=8,  z7(0)=5,  z8(0)=6,  z9(0)=8,  z10(0)=3 .
53

Рис.3. Графики решений системы (5)-(7) первого типа

54 Из за симметричности правил приема и отправки грузов (отправка грузов с i -й станции на i+1 -ю станцию с определенной интенсивностью равносильна приему грузов на i+1 -й станции с той же интенсивностью) и условия (13) следует, что каждая станция принимает груз с предыдущей станции и отправляет на следующую станцию с постоянной и одинаковой интенсивностью т.е. в системе грузоперевозок, как и в случае идентичных станций, устанавливается стационарный поток (рис.4).
55

Рис.4. Интенсивность потока в общем случае

56 Величина данного потока зависит как от параметров модели, так и от начальных значений, а точнее от их среднего арифметического ( c- ). С помощью численных экспериментов исследуем эти зависимости.
57 Начнем с параметра α , который характеризует потенциал первой технологии. На рис.5 приведен график зависимости интенсивности стационарного потока от α при фиксированных значениях остальных параметров и начальных значений (см. (14)).
58

Рис.5. Зависимость интенсивности стационарного потока от параметра α

59 Как видим с увеличением параметра α увеличивается интенсивность стационарного потока, однако скорость ее роста уменьшается. Отметим, что при этом уменьшаются и слагаемые в правой части системы (5)-(7), отвечающие за прием и отправку грузов в рамках первой технологии. Это означает, что, имея большой потенциал первой технологии, его невозможно реализовать при сильно ограниченных пропускных способностях станций.
60 Перейдем к исследованию зависимости интенсивности стационарного потока от характеристик пропускной способности станций, т.е. от ai и Δi . Начнем с параметра ai . На рис. 6 приведена зависимость интенсивности стационарного потока от a1 при фиксированных значениях остальных параметров и начальных условий.
61

Рис.6. Зависимость интенсивности стационарного потока от a1 .

62 Как видно, с увеличением a1 увеличивается поток, причем вначале скорость роста является возрастающей, однако начиная с некоторого значения она становится убывающей. Такая же тенденция наблюдается для других параметров ai . Таким образом, увеличение эффективности использования путей на некоторой выбранной станции до определенного уровня дает неплохой вклад в увеличение интенсивности потока во всей системе грузоперевозок. Однако, в силу ограниченности пропускной способности на остальных станциях, дальнейшее увеличение эффективности использования путей на одной станции уже не может привести к ощутимому вкладу в увеличение интенсивности грузопотока.
63 Перейдем к исследованию зависимости стационарного потока от другой характеристики пропускной способности станций - числа путей на них ( Δi ). На рис. 7 приведена эта зависимость для числа путей на первой станции при фиксированных значениях остальных параметров и начальных условий. Как видно данная зависимость является линейной. Напомним, что пропускная способность станций зависит от количества путей согласно формуле (12). Такая же зависимость наблюдается и отношении других параметров.
64

Рис.7. Зависимость интенсивности стационарного потока от числа путей на 1-ой станции

65 Наконец, исследуем зависимость интенсивности стационарного потока от начальной загруженности станций ( c- ). Как показали эксперименты эта зависимость является убывающей. На рис.8 приведен график этой зависимости при фиксированных значениях параметров.
66

Рис.8. Зависимость интенсивности стационарного потока от начальной загруженности станций

67 Отметим, что максимальная начальная загруженность станций, равная 10.3 соответствует ситуации, когда в начальный момент времени на всех станциях все пути задействованы, т.е. zi(0)=Δi,   i=1,2,...,n . В этом случае организация грузопотока в рамках данной модели невозможна. Она также невозможна в случае, когда в начальный момент времени на всех станциях все пути свободны, т.е. zi(0)=0,   i=1,2,...,n .
68

ЗАКЛЮЧЕНИЕ

69 Исследована модель организации грузоперевозок по круговой цепочке станций. Движение грузопотока осуществляется с помощью двух технологий, задающих интенсивность грузопотока в зависимости от пропускной способности станций и перегонов. Такая модель описывается системой дифференциальных уравнений с ограничениями, описывающими допустимый диапазон изменения числа задействованных путей на станциях. Эти ограничения могут порождать разрывы в правых частях уравнений. В силу сложности получения аналитического решения данная система была исследована численно. Проведенные эксперименты позволили выявить указанные разрывы. Они возникают вследствие того, что со временем некоторые компоненты решения начинают колебаться у левой границы указанных выше диапазонов (нуля) периодически выходя на нее. Такое поведение характерно для решений с малыми начальными значениями и решениями, полученными при значениях параметра α , близкого к нулю, т.е. в случае недостаточной обеспеченности станций грузом в начальный момент времени и низкой пропускной способности перегонов. Отметим, что данное свойство решений не распространяется на частный случай системы, описывающей модель с идентичными станциями (одинаковое число путей на всех станциях и эффективность их использования). В этом случае независимо от характеристик станций и перегонов, а также их начальной загруженности, в системе грузоперевозок устанавливается стационарный поток.
70 Разрывные решения приводят к большим колебаниям в величине потока, возникающего в системе грузоперевозок. На практике такой режим грузоперевозок крайне сложно реализуем, поэтому, более детально изучены режимы грузоперевозок, определяемые непрерывными решениями. Как показали численные эксперименты, такие режимы грузоперевозок приводят к установлению стационарного потока в системе грузоперевозок. Исследована зависимость величины такого потока от параметров модели, определяющих пропускную способность перегонов и станций, а также от загруженности станций в начальный момент времени. С помощью такой модели, предварительно оценив ее параметры на основе статистических данных, можно прогнозировать эффект, полученный от мероприятий, связанных с увеличением пропускной способности той или иной станции.

Библиография

1. Бекларян Л.А., Н.К. Хачатрян. (2013). Об одном классе динамических моделей грузоперевозок // Журнал вычислительной математики и математической физики. Т. 53. № 10. C. 1649–1667.

2. Бекларян Л.А., Хачатрян Н.К. (2019). Динамические модели организации грузопотока на железнодорожном транспорте// Экономика и математические методы. Т. 55. № 3. С. 62-73.

3. Гасников А.В., Кленов С.Л., Нурминский Е.А., Холодов Я.А., Шамрай Н.Б. (2013). Введение в математическое моделирование транспортных потоков. Под ред. Гасникова А.В. М.: МЦНМО.

4. Иносэ Х., Хамада Т. (1983). Управление дорожным движением. М.: Транспорт.

5. Сухинова А.Б., Трапезникова М.А., Четверушкин Б.Н., Чубарова Н.Г. (2009). Двумерная макроскопическая модель транспортных потоков // Математическое моделирование. Т. 21. № 2. С. 118–126.

6. Швецов В.И. (2003). Математическое моделирование транспортных потоков // Автоматика и телемеханика. № 11. С. 3–46.

7. Bando M., Hasebe K., Nakayama A., Shibata A., Sugiyama Y. (1995). Dynamical model of tra?c congestion and numerical simulation // Physical Review. E. Vol. 51. P. 1035–1042.

8. Bar-Gera H. (2002). Origin-based algorithm for the traffic assignment problem// Transportation Science. Vol. 36. No. 4. P. 398–417.

9. Beklaryan L.A., Khachatryan N.K. (2006). Traveling Wave Type Solutions in Dynamic Transport Models // Functional Differential Equations. Vol. 13. No. 12. P. 125–155.

10. Beklaryan Levon A., Khachatryan Nerses K., Akopov Andranik S. (2019). Model for organization cargo transportation at resource restrictions// International Journal of Applied Mathematics. Vol. 32. No. 4. P. 627-640.

11. Brackstone M., McDonald M. (2000). Car following: A historical review // Transportation Research. F. Vol. 2. P. 181–196.

12. Carrothers G. A. P. (1956). An historical review of the gravity and potential concepts of human interaction // J. American Instit. Planners. Vol. 22. P. 94–102.

13. Chowdhury D., Santen L., Schadschneider A. (2000). Statistical physics of vehicular tra?c and some related systems // Physics Reports. Vol. 329. P. 199–329.

14. Cremer M., Ludwig J. (1986). A fast simulation model for tra?c ?ow on the basis of Boolean operations // Mathematics and Computers in Simulation. Vol. 28. P. 297–303.

15. Daganzo C. F. (1994). The cell transmission model: A dynamic representation of highway tra?c consistent with the hydrodynamic theory // Transportation Research. B. Vol. 28. P. 269–287.

16. Daganzo C. F. (1995). The cell transmission model, Part II: Network tra?c // Transportation Research. B. Vol. 29. P. 79–93.

17. Fotheringham A. S. (1983). A new set of spacial-interaction models: the theory of competing destinations // Environment and Planning. A. Vol. 15. P. 15–36.

18. Fotheringham A. S. (1986). Modelling hierarchical destination choice // Environment and Planning. A. Vol. 18. P. 401–418.

19. Harris B., Wilson A. G. (1978). Equilibrium values and dynamics of attractiveness terms in production-constrained spatial-interaction models // Environment and Planning. A. Vol. 10. P. 371–388.

20. Helbing D., Treiber M. (1998). Gas-kinetic-based tra?c model explaining observed hysteretic phase transition // Physical Review Letters. Vol. 81. P. 3042–3045.

21. Khachatryan N.K., Akopov A.S. (2017). Model for Organizing Cargo Transportation with an Initial Station of Departure and a Final Station of Cargo Distribution // Business Informatics. No. 1. P. 25–35.

22. Khachatryan N.K., Akopov A.S., Belousov F.A. (2018). About Quasi-Solutions of Traveling Wave Type in Models for Organizing Cargo Transportation // Business Informatics. No. 1 (43). P. 61–70.

23. Leventhal T., Nemhauser G.L., Trotter L. (1973) A column generation algorithm for optimal traffic assignment // Transportation Science. No 7. P. 168–176.

24. Lo H.K., Chen A. (2000) Traffic equlibrium problem with route-specific costs: formulation and algorithms// Transportation Research. B. Vol. 34. No 6. P. 493–513.

25. Nelson P. (1995). A kinetic model of vehicular tra?c and its associated bimodal equilibrium solutions // Transport Theory and Statistical Physics. Vol. 24. P. 383–409.

26. Popkov Yu. S. (1995). Macrosystems theory and its applications. Berlin: Springer Verlag.

27. Prigogine I., Herman R. (1971). Kinetic Theory of Vehicular Tra?c. N.Y.: Elsevier.

28. Shvetsov V.I. (2009) Algorithms for distributing traffic flows// Automation and Remote Control. Vol. 70. No 10, P. 1728–1736.

29. Spiess H., Florian M. (1989). Optimal strategies: a new assignment model for transit networks // Transportation Research. B. Vol. 23. P. 83–102.

30. Treiber M., Hennecke A., Helbing D. (2000). Congested tra?c states in empirical observations and microscopic simulations // Physical Review. E. Vol. 62. P. 1805–1824.

31. Wilson A. G. (1971). A family of spatial interaction models and associated developments // Environment and Planning. A. Vol. 3. P. 255–282.

Комментарии

Сообщения не найдены

Написать отзыв
Перевести