Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Шабалин Н.Н. Оптимизация процесса переработки вагонов на станциях

.pdf
Скачиваний:
9
Добавлен:
23.10.2023
Размер:
6.13 Mб
Скачать

 

 

 

 

Рис.

16.

Гисто­

 

 

 

 

грамма

статисти­

 

 

 

 

ческого

(ступенча­

 

 

 

 

тая линия) и кри­

 

 

 

 

вая

нормального

 

 

 

 

распределения го­

 

 

 

 

рочного

интервала

О

5

Ю

15

20 tr.MUH

 

 

Судя по небольшой величине коэффициента

вариации,

распределение горочных интервалов может быть аппрок­ симировано кривой, нормального распределения. В табл. 5 приведен порядок проверки гипотезы о нормальном

распределении

горочного интервала при помощи

крите­

рия Пирсона х2 .

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 5

Последовательность

вычисления критерия согласия

 

для распределения

горочного

интервала

 

"i

рт=ыцп

л т

= пРт

ni~'4

( « » - " т ) 2

( " £ - Л т ) 2

Л Т

 

 

 

 

 

 

 

5

10

0,009

 

17

7

49

2,9

6

40

0,018

 

35

5

25

0,72

7

70

0,033

 

64

6

36

0,56

8

100

0,055

 

106

6

36

0,36

9

180

0,081

 

156

24

576

3,7

10

200

0,106

204

4

16

0,07

11

240

0,126

242

2

4

0,01

12

260

0,133

256

4

16

0,06

13

250

0,126

 

242

8

64

0,26

14

190

0,106

204

14

196

0,96

15

140

0,081

 

156

16

256

1,64

16

120

0,055

 

106

14

196

1,85

17

50

0,033

 

61

14

196

3,07

18

40

0,018

 

35

5

25

0,72

19

20

0,009

 

17

3

9

0,53

20

10

0,004

 

8

2

4

0,50

Итого

1920

 

Х2 =17,91

50

В первой вертикальной графе представлены целочис­ ленные значения горочного интервала, которые являются серединами разрядов, каждый из которых равен 1. Во второй графе проставлена статистическая частота каж­ дого разряда. В третьей графе определяется теоретиче­ ская вероятность попадания в данный разряд как произ­ ведение величины разряда на плотность распределения вероятностей по предполагаемому нормальному закону распределения, т. е. Pv=Atf(t). Здесь Д £ = 1 , а значение плотности распределения для каждого значения t на­ ходится при помощи таблицы нормированного рас­ пределения по зависимости

/(о=4-/'^),

где / ' ( — н а х о д и т с я по таблице (например, по табл. 3

книги Вентцель Е. С. «Теория вероятностей»). Для полученных данных

/w-J-z-f-^)-

Так, для * = 6 /'(?=!?) = 0,054 и /(0=0,01 8 Р т =0,018

(см. рис. 16). Произведение объема выборки 1920 на тео­ ретическую вероятность дает значение теоретической часто­ ты каждого разряда, что и отражено в графе 4. В гра­ фах 5 и 6 помещены соответственно отклонения и квадраты отклонений теоретической частоты от статистической. В 7-й графе показаны отношения квадратов отклонений к теоре­ тической частоте. Сумма графы 7 дает значения X2 .

Число степеней свободы равно числу разрядов минус число наложенных связей, что в нашем примере опреде­ лится г=1—3=16—3 = 13. По специальной таблице [2] находим . Р=0,163 . Следовательно, выравнивание с по­ мощью нормального распределения согласуется со стати­ стическими данными. Это видно и по рис. 16, где на гисто­ грамме нанесена кривая распределения. Если по табли­ це получим значение Р < 0 , 0 5 , тогда принятая гипотеза отвергается. Как отмечалось выше, анализ интервалов обслуживания различных видов: расформирования, фор­ мирования, обработки в парке прибытия и отправления и т. д. показывает, что коэффициент вариации их чаще всего не превышает 0,4—0,5, т. е. значительно ниже, чем

51

коэффициенты вариации входящих потоков требований. Это объясняется тем, что операции обслуживания явля­ ются результатом определенной технологии и находятся в руках определенных исполнителей или функционирова­ ния определенных устройств и по сравнению с потоком требований более управляемы. Потоки требований под­ вержены большему влиянию различных объективных факторов с различной частотой и дальностью действия,, что и обусловливает значительно большую их неравно­ мерность по сравнению с потоками обслуживания.

3. ПОКАЗАТЕЛИ ВЗАИМОДЕЙСТВИЯ МЕЖДУ ПРОЦЕССАМИ ПРИБЫТИЯ И РАСФОРМИРОВАНИЯ ПОЕЗДОВ

Успешность процесса расформирования поездов, как и любой системы обслуживания, зависит от перераба­ тывающей способности горки (или вытяжек формирова­

ния,

если

станция

безгорочиая), объема предлагаемой

работы и

степени

неравномерности поступления

поез­

дов

и неравномерности их расформирования. При

всех

условиях средний горочный интервал на расформирова­ ние одного состава должен быть меньше среднего интер­ вала прибытия поездов. В зависимости от соотношения интенсивности прибытия и расформирования поездов в

различные

периоды времени образуются простои в ожи­

дании

расформирования. Величины очередей

зависят

от степени

неравномерности как

потока требований (по­

ездов),

так

и их обслуживания

(расформирования), а

скорость устранения очередей зависит от наличия резер­ ва мощности обслуживающего устройства (горки). От­ сюда возникает задача оценки и учета степени неравно­ мерности прибытия и расформирования при определе­ нии как потребной мощности устройств, так и показа­ телей процесса переработки вагонов.

Для получения средних показателей процесса рас­

формирования

составов

необходимо

его рассматривать

в стационарном

(установившемся)

режиме,

которому

соответствует средняя интенсивность

прибытия

поездов

и средняя интенсивность

их расформирования. Интен­

сивность прибытия поездов можно характеризовать па­ раметром потока X, равным среднему числу поездов, поступающих в час. Интенсивность расформирования

52

зависит от мощности горки и может быть охарактери­ зована средним числом поездов v, которое горка рас­ формировывает за час.

Для начала возьмем наиболее простой случай, когда интервалы между требованиями и интервалы обслужи­ вания независимы и имеют показательное распределе­ ние. В этом случае процесс обслуживания рассматрива­ ется как марковский процесс, при котором вероятность состояния системы в будущем зависит только от состоя­ ния ее в настоящий момент и не зависит от способа, ка­

ким система пришла в это состояние. Под

состоянием

системы

здесь подразумевается число поездов, находя­

щихся в

системе расформирования. Процесс Марко­

ва А. А.,

кроме того, обладает

свойством эргодичности,

согласно

которому существует

предельный

установив­

шийся режим вероятностного процесса, не зависящий от начального состояния системы. Это значит, какое бы ни было состояние системы в данный или начальный мо­ мент, через определенное число шагов переходного про­ цесса система достигает определенного стационарного состояния, определяемого соотношением интенсивности прибытия и расформирования.

Потоки прибытия поездов и их расформирования рассматриваются как простейшие со свойствами стаци­ онарности, ординарности и отсутствием последействия.

Для этих условий рассмотрим вероятности измене­ ния состояний однолинейной системы расформирования за элементарный отрезок времени dt под влиянием пото­ ка требований и обслуживания. Сумма вероятностей раз­

личных состояний системы как полной группы

событий

за период dt равна единице

 

 

P0(dt)

+ Pt(dt)

+

Pz(dt) + ... +

Pa(dt) =

1,

где P0(dt),

...,Pn(dt)

вероятность

появления

0, 1, 2,

 

 

 

п поездов за малый промежу­

 

 

ток dt.

 

 

В соответствии со свойством ординарности потока ве­ роятность появления двух и более событий за малый про­ межуток времени есть величина бесконечно малая посравнению с этим промежутком. В связи с этим практи­ ческое значение имеют только два слагаемых, а осталь­ ными можно пренебречь и написать

P0(dt) + P1(dt)=l.

53

Д ля стационарных потоков

 

 

 

 

 

 

 

 

где X

 

 

P i (dt)

=

Mt,

 

 

 

 

 

интенсивность

потока прибытия поездов.

 

С учетом

предыдущего

равенства

вероятность

того,

•что за период

dt

не прибудет

поезд,

определится

 

 

 

 

Pa(dt)

=

1 — Mt.

 

 

 

 

Аналогично получим вероятности

изменения

системы

под влиянием потока обслуживания

(расформирования).

Вероятность того, что за период

dt

один

поезд

будет

расформирован, т. е. в системе

станет

на

один

поезд

меньше,

определится

 

 

 

 

 

 

 

 

 

 

 

 

Р _ ! (dt)

=

vdt.

 

 

 

 

 

Вероятность того, что

за период dt

не

будет

расформи­

рован ни один поезд

 

 

 

 

 

 

 

 

 

 

 

 

Р _ 0

(dt)

=

1 -

vdt.

 

 

 

 

Теперь надо еще определить вероятность того, что под

влиянием

потока

требований

и обслуживания

система

не изменит своего состояния

за

элементарный

отрезок

времени. Это может произойти практически двумя спо­ собами. За это время ни один поезд не прибудет и не будет расформирован или за это время один поезд при­ будет и один будет расформирован, т. е.

Я п п

(dt) =

(1 -

Ш)

(1 -

vdt)

+ (Ш)

(vdt).

 

Вероятность

оставшихся

событий

имеет

порядок

o(dt).

Если раскрыть скобки

и пренебречь квадратами элемен­

тарных отрезков времени (dt)2,

получим вероятность то­

го, что система не

изменит

своего

состояния

 

 

рпп

=

1 \dt

vdt.

 

 

Здесь определены вероятности изменения системы

расформирования

под

влиянием интенсивности

потока

прибытия и интенсивности расформирования за элемен­ тарный отрезок dt.

Вероятности изменения состояний системы за период dt можно представить наглядно в виде графа переходов из одного состояния в другое, показанного на рис. 17. Как видно, под влиянием интенсивности прибытия поез­

дов число их в системе возрастает, а под

влиянием ин­

тенсивности

расформирования — уменьшается. Под

чис­

лом требований в системе подразумевается

как

число

требований,

находящихся в ожидании обслуживания,

54

 

 

 

1-Ш

t-Mt-M

l-Adt-ki

t-м-Ы

H**))ttt

Рис. 17. Граф перехода системы из одного состояния в другое

так и те, которые обслуживаются. В рассматриваемой од­ нолинейной системе расформирования на горке один состав может находиться в процессе роспуска, а ос­ тальные в очереди.

Для нахождения средних значений числа

требований

в системе или в очереди и других показателей

необхо­

димо знать, как распределены

вероятности

различных

состояний системы, т. е. с какой

вероятностью в

систе­

ме будет О, 1, 2,...,п требований.

Аналитически

состоя­

ние системы расформирования в момент t

и изменения"

за период dt можно представить

системой

дифференци­

альных уравнений [ 7 ] . Вероятность того, что в системе расформирования в конце периода t-\-dt не будет ни од­ ного состава

Р0 (t + dt)

= Я 0 ( 0 ( 1 -

МО + Я ; (0 vdt.

 

 

Первое слагаемое

определяет

одновременное

совпаде­

ние вероятности того, что в момент t в системе

было

нуль составов, и вероятности

того, что за период

dt ни

один поезд не прибыл. Второе слагаемое определяет вероятность того, что в момент t был один состав и за период dt один состав был расформирован.

Вероятность того, что в системе расформирования в конце периода t-\-dt будет п составов, может образо­ ваться тремя способами

Pn(t

+ dt) = Pn(t)(l—ldt—vdt)

+ Pn-iWkdt

+

Pn+xiQvdt.

Здесь первое слагаемое — вероятность

того,

что

в мо­

мент t в системе было ровно п поездов

и за

период

dt

не

изменилось состояние

системы; второе слагаемое

вероятность того, что в момент t в системе

было п—1 со­

ставов, но за период dt

поступил один

состав;

третье

слагаемое — вероятность того, что в момент t в систе­ ме было п + 1 составов, но за период dt один состав расформирован.

55.

Если в этих уравнениях раскрыть скобки, перенести соответственно Po(t) и Pn(t) влево и разделить на dt, то получим систему дифференциальных уравнений:

^= - x p 0 ( 0 + VPI(0;

^±р-

= _

(х + v) р„

(t) + Xря_,

(0 + v

(

0 ,

л >

1;

при начальном

условии:

Р о ( 0 ) =

1, Р п ( 0 )

=

0,

л > 1 .

Рас­

сматривая предельный установившийся режим вероятност­ ного процесса, при котором в силу эргодичности процесса

Pn{t)~+PnVi -»-0 (^-»-со), и приравняв производные

нулю, получим для стационарного состояния алгебраичес­ кие уравнения:

ЬРо = v P i ;

(X + v)P„ = X/V-j + v/>„+

Для упрощения разделим оба уравнения на v и получим

Л = Ро^г = Ро?1

Приняв т г = 1 , получим из второго уравнения

 

(1 + Р ) Р 1

= Р 2 + р Р 0 .

 

 

Подставив сюда Рг

из первого уравнения, получим Р 2 = Я 0

р 2

и в общем

случае

 

 

 

 

 

 

 

 

Рп = Р0?п;

 

(19)

здесь

р = — . = _

=

_ =

Х ^ г < 1 ,

 

где W — число поездов, поступающих в

расформирова­

ние за сутки;

 

 

 

 

 

/ —

средний

интервал

прибытия поездов;

 

tv

средняя

величина

горочного

интервала

на

 

расформирование

одного

состава.

 

Величина р определяет уровень загрузки горки, при­ веденную плотность потока требований, определяющую, сколько поступит поездов за время расформирования одного состава. Эта величина определяет степень опере­ жения темпом расформирования темпа прибытия поез-

.56

ДОБ. Естественно, что в условиях неравномерности за­ грузка горки, как и любого обслуживающего устройст­ ва, должна быть меньше единицы. При заданных раз­ мерах прибытия поездов в расформирование уровень загрузки горки можно снижать сокращением горочного интервала.

Полученная зависимость (19) позволяет определять вероятность различных состояний любой однолинейной системы обслуживания с ожиданием и бесконечной оче­ редностью на входе. Сумма вероятностей всех возмож­ ных состояний как полной группы событий

liPn = 21ЯоР л = 1.

Отсюда вероятность того, что система расформирования свободна

Здесь знаменатель' представляет собой сумму членов геометрической убывающей прогрессии ( 1 + р + р 2 + р 3 + -f-...+р'1 ), которая равна 1/(1—р). Отсюда Рп0рп

=( 1 - р ) р »

Значение Ро определяет резерв мощности горки. Сор­ тировочная горка может быть или занята, вероятность

чего равна

р, или

свободна, вероятность чего равна

1—р = Ро-

Составы,

поступившие в период занятости

горки, будут ожидать расформирования, пока горка не освободится, а составы, поступившие в период свободности горки, такого простоя иметь не будут.

Из зависимости (19) следует, что распределение ве­ роятностей различных состояний системы расформиро­ вания определяется уровнем загрузки горки р и ее ре­ зервами мощности Ро. Чем меньше уровень загрузки горки, тем больше вероятность того, что она будет сво­ бодна пли в системе будет небольшое число составов. Чем больше уровень загрузки горки, тем больше веро­ ятности значительного числа составов в системе рас­ формирования, а значит и в очереди. Это особенно на­ глядно видно из рис. 18, где показано распределение вероятностей нахождения числа составов в системе рас­ формирования, как и в любой однолинейной системе-

 

 

 

 

Рис.

18.

Зависи­

 

 

 

 

мость

вероятностей

 

 

 

 

состояния

системы

 

 

 

 

от

уровня

загруз­

 

 

 

 

 

ки

устройства

обслуживания,

при уровнях загрузки

р = 0,6; 0,8 и 0,9.

При р = 0,6 вероятность

того,

что в

системе

нет требо­

ваний, равна

Р = 1 — 0 , 6 = 0,4

и дальше

вероятность

наличия одного, двух и т. д. требований быстро

падает.

Но при уровне загрузки

р = 0,9 вероятность

того, что в

•системе нет требований,

составляет только 0,1 и вероят­

ности последовательно уменьшающегося числа составов снижаются довольно медленно. Ниже в табл. 6 приведе­

ны численные

значения

вероятностей

различных

состоя­

ний системы

обслуживания

при различных уровнях

 

их

загрузки.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

6

 

Распределение вероятностей

числа требований

 

 

 

 

в системе обслуживания

 

 

 

 

 

 

Число требований в системе

обслуживания я

 

 

 

Уровень

 

 

 

 

 

 

 

 

 

 

 

загрузки р

0

1

2

3

4

5

6

7

8

 

9

 

 

0,3

0,70

0,21

0,06

0,02

0,01

 

 

 

 

 

 

0,5

0,50

0,25

0,12

0,06

0,03

0,02

0,01

 

 

 

 

0,7

0,30

0,21

0,15

0,10

0,07

0,05

0,03

0,024 0,02

0,01

0,9

0,10

0,09

0,08

0,07

0,06

0,059 0,053 0,047 0,043 0,039

.58

Знание распределения вероятностей различных со­ стояний системы позволяет определить среднее числа требований в системе обслуживания, среднее число тре­ бований, ожидающих обслуживания, а также время на­ хождения в них и других показателей процесса.

Среднее число требований в системе обслуживания равно сумме произведений каждого возможного числа требований п на их вероятность Рп, т. е.

со со

 

U

= У, пРп = Р0 ^

яр" = j

~

(20)'

 

 

л = 0

л = 0

 

 

 

 

Здесь

2 пРп

=

Р + V + 3 Р 3 +

- = clif2

 

 

п=0

 

 

 

 

Р

 

 

 

 

 

 

 

 

 

 

Следовательно,

среднее

число составов в

системе-

расформирования зависит от уровня загрузки горки.

Среднее число

составов,

ожидающих

расформирова­

ния (средняя длина очереди), определится как сумма

произведений всех возможных вероятностей Рп

на чис­

ло составов, ожидающих расформирования

(п—1).

Ожидать расформирования будут все составы, кроме одного, находящегося в процессе роспуска,

со со

L 0 K = yi(ti-\)Pn

= PoP^in-l)?"-1

= £ 7 .

(21)

л = 2

л = 2

 

 

со

 

 

 

Здесь 2 (, г -

^ Р " - 1 = Р + 2р2 + Зр3 + ... =

,

л = 2

V

? )

как и в первом случае при определении среднего числа составов в системе расформирования.-

Среднее число составов, находящихся в системе рас­ формирования, равно сумме средней длины очереди и среднего числа составов, находящихся в процессе рос­ пуска, которое равно уровню загрузки горки, т. е.

Lc = L0 + p.

Из зависимостей (20) и (21) следует, что среднее число составов в системе расформирования и в очереди

обратно

пропорционально резерву мощности

горки

(1—р).

Если в формуле (21) уровень загрузки

заме­

нить через интенсивности потока поездов и их расфор­ мирования, то она примет вид

59

Соседние файлы в папке книги из ГПНТБ