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

Транспортная логистика

.pdf
Скачиваний:
0
Добавлен:
30.11.2025
Размер:
10.85 Mб
Скачать

ны в матричные. С помощью моделей производится оптимизация маршрутов транспортных потоков. С этой целью используется извест­ ная «транспортная задана линейного программирования»:

целевая функция

тп

;=1 7=1

при условиях

т

п

X ХЧ = B j

X ХУ = А‘ ХУ - 0 ’

,=1

м

где - стоимость перевозки груза из пункта / в пункт j на единицу перевозимого груза;

Ху - искомая величина грузопотока из пункта / в пунктj; Bj - количество получаемых грузов в пунктахj;

А, - количество отправления грузов из пунктов /;

1

, 2

пункты отправления груза;

1

, 2„..J,...,п - пункты назначения груза.

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

8.2.1.Метод потенциалов для решения транспортной задачи

вматричной форме с ограничениями пропускной способности

Задача 1

1. Построить оптимальный план перевозок каменного угля с пя­ ти станций Aj (i = 1.. .5), обслуживающих шахты, до девяти крупных потребителей, имеющих подъездные пути В) (/' = 1,..., 9).

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

2. Определить объем тонно-километровой работы для всех пла­ нов начального и оптимального планов перевозки груза.

280

Исходные данные

Данные о наличии ресурсов на станциях отправления А; приве­ дены в табл. 8.3; данные о размерах прибытия груза Bj на станциях назначения - в табл. 8.4. Расстояние перевозки от каждой ;-й стан­ ции отправления до каждой у'-й станции назначения указано в форме 1 в правом верхнем углу каждой клетки. В левом верхнем углу формы 1 указаны ограничения пропускной способности, если тако­ вые имеются. Матрица расстояний и ограничений пропускной спо­ собности принимается одинаковой для всех вариантов.

Таблица 8.3 Ресурсы станции отправления Аь тыс. т____________

Номер стан­

 

 

 

 

Варианты

 

 

 

 

ции отправ­

 

2

 

 

5

6

7

8

9

0

ления («ро­

1

3

4

ка матрицы)

150

150

150

145

155

400

400

150

400

145

А\

а 2

150

145

150

155

150

155

150

155

155

400

As

145

155

155

150

400

145

150

150

150

150

л 4

155

400

145

150

150

150

155

145

150

150

А5

400

150

400

400

145

150

145

400

145

155

ИТОГО

1000

1000

1000

1000

1000

1000

1000

1000

1000

1000

отправлено

 

 

 

 

 

 

 

 

 

 

Таблица 8.4 Объем прибытия груза В, на станции назначения, тыс. т

Номер стан­

 

 

 

 

Варианты

 

 

 

 

ции отправ­

 

 

 

 

 

 

7

 

9

0

ления (стро­

1

2

3

4

5

6

8

ка матрицы)

100

90

100

100

100

85

85

115

95

80

в,

Вг

100

90

100

100

100

85

100

120

80

115

в3

І 0 0 І

110

95

100

100

80

120

80

100

80

В4

100

110

95

150

100

80

120

85

115

95

Bs

100

100

105

90

100

100

80

100

80

100

в.

150

150

155

90

І00

100

110

50

170

50

В-,

100

100

100

110

90

120

80

165

85

155

В*

100

100

100

110

90

120

80

120

50

120

Д,

150

150

150

150

220

230

225

165

225

205

ИТОГО

1000

1000

1000

1000

1000

1000

1000

1000

1000

1000

отправлено

 

 

 

 

 

 

 

 

 

 

281

Методические указания к выполнению задачи 1

1. Решение задачи начинают с составления исходной матрицы (форма 1). Для решения задачи в столбец Л/ заносят объемы ресур­ сов станций отправления, а в строку Bj - объемы прибытия грузов на станции назначения.

После записи ресурсов и потребностей груза по своему варианту на исходную матрицу для решения задачи строят начальный план любым известным способом (к примеру, методом минимального элемента).

После составления начального плана необходимо проверить ба­ ланс по строкам и столбцам матрицы. Число базисных клеток (ба­ зисной клеткой называется клетка, имеющая корреспонденцию, т. е. объем перевозки из i-й станции отправления на j станцию назна­ чения), должно быть:

К - т + п ~ \ ,

где т - число строк матрицы; п - число столбцов.

Для рассматриваемой матрицы АТ = 9 + 5 - 1 = 13. Если число ба­ зисных клеток больше числа К, то начальный план составлен не­ правильно. Клетки, где величина перевозки равна пропускной спо­ собности, называются небазисными.

В форме 1 в базисных клетках корреспонденция перевозки отме­ чена жирным шрифтом.

Форма 1

Пример построения начального плана для задачи № 1 на матрице

X

Л ,~ « 5

Вг ~ 100

В , - 120

Д ,”

120

S i - 8 0

 

«SO

B*>~ SO

& »«225

 

 

80

40

90

 

103

150

30 50

?S

80

 

90

 

«4,-iSO

 

 

 

*000

 

 

 

 

1 0 *

 

30

 

 

 

3 0

 

 

 

 

 

 

 

 

 

 

 

 

Ю

30

43

 

40

25

65

 

*0

10

so

ІШ5

 

 

т

 

 

 

 

3S

 

3C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 » -1 5 5

10

20

35

75

 

160

90

80

70

40

 

60

1030

 

 

 

 

 

 

 

 

 

SB

 

75

 

 

 

 

 

 

 

 

 

 

 

 

Л 4 - И 5

 

45

5

15

 

30

IW

40

75

JO

 

20

1070

 

 

t m

 

 

 

 

 

 

 

 

43

 

 

 

 

 

 

 

 

 

 

 

 

Л . -»400

 

І5

15 25

10

20

*5

25

 

iso

70

 

«5

, « , f

 

 

 

130

20

4&

e a

5 8

 

 

U

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

*5

101S

1075

I05S

1105

І030

lo s s

1025

1070

1

0 )

 

 

 

 

 

....

 

 

 

 

 

 

 

 

 

282

Объем тонно-километровой работы равен 42725 т-км.

2. Оптимальный план перевозок на заданной матрице составим по методу потенциалов.

Любой допустимый план является оптимальным тогда, когда каждой строке и каждому столбцу могут быть присвоены некото­ рые числа U, и У„ называемые потенциалами и отвечающие сле­ дующим условиям:

 

Уг -и ,< с ц

д л я

лг„ = 0;

(8.1)

 

Vі - Ui = Су

д л я

0 < Хц < di{.

(8.2)

 

Vj -U j> с,,

д л я

Хц = dih

(8.3)

где Vj -

потенциал j -го столбца;

 

 

 

U, -

потенциал /-й строки;

 

 

 

c,/ -

расстояние перевозки от і-го поставщика до /-го потребителя;

х,у -

корреспонденция (размер перевозки) от г-го поставщика до

у'-го потребителя;

 

 

 

d,j -

величина пропускной способности.

 

3.

Присвоение потенциалов начинаем со строки, в которой среди

базисных клеток имеется максимальное расстояние. Этой строке присваиваем потенциал, равный 1000. Затем, используя условие оп­ тимальности (8.2), находим потенциалы остальных строк и столб­ цов следующим образом:

для у-го столбца

для /-Й строки

4. После присвоения всем строкам и столбцам потенциалов определяем, имеются ли нарушения неравенств (8.1) и (8.3), по форійуле

Ht/ = Vj - Ui - Cij.

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

283

5. Улучшение допустимого плана начинаем с клетки, имеющей

максимальное (по модулю) нарушение Z/™3*. Для этой клетки строим

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

Следует заметить, что для каждой клетки с нарушением сущест­ вует только один контур улучшения плана. Нумерация клеток кон­ тура начинается с клетки с нарушением. Если клетка с нарушением свободная, то ей присваиваем № 1. Для клеток с поставками, рав­ ными пропускной способности, нумерация начинается с нуля. Да­ лее номера присваиваются по ходу контура. Число клеток в контуре всегда четное.

В найденном замкнутом контуре определяем корреспонденцию улучшения допустимого плана на данном этапе решения. Коррес­ понденция улучшения плана находится следующим образом:

* У = m in [х0четн, ( d ij - * „ ) _ , ] .

На величину хул изменяются все корреспонденции контура, на­ чиная с клегки с нарушением: уменьшаются корреспонденции, за­ писанные в четных клетках, и увеличиваются корреспонденции, записанные в нечетных клетках контура.

6. После пересмотра корреспонденции необходимо пересоста­ вить систему потенциалов всей матрицы и проверить соблюдение условия оптимальности (8.1) и (8.3). Если небазисные клетки удов­ летворяют этим условиям, то найдено оптимальное решение. Если имеются нарушения условий оптимальности, то расчет на матрице следует продолжить до тех пор, пока все клетки матрицы не будут удовлетворять условиям (8.1)-(8.3). Число нарушений и их величина всегда стремятся к нулю.

Используя вышеизложенный алгоритм и произведя необходимое количество итераций, получим оптимальный план перевозок (в ле­ вом верхнем углу каждой клетки находится величина нарушения, если таковая имеется), (табл. 8.5—8.12).

284

Оттенком I

I отмечены четные клетки замкнутого контура;

1 \

I - нечетные;

[

I - клетки с максимальным (по модулю) на­

рушением;

■небазисные клетки.

Таблица 8.5 Построение оптимального плана перевозок. Шаг 1

 

f t -

*5

f t

-100

f t * 120

f t -

120

By-9 0

f t - 110

X

 

 

 

 

 

 

 

 

 

 

J50

 

90

35

40

90

 

 

 

150

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

l

 

 

І Ш

Лгш150

 

10

 

30

45

 

 

 

 

65

ts

 

 

 

 

I T

 

 

 

 

 

 

 

 

 

 

ш

і і

 

 

 

 

 

 

 

 

 

 

Aj * 155

25

20

)0

35

75

 

l«0

 

90

SO

 

 

 

 

 

 

 

 

 

 

Л -М 5

 

45

 

5

35

5

30

 

no

40

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

400

55

15

45

25

10

 

 

 

 

SO

 

 

 

 

m

 

 

w

m

* •

 

 

 

 

 

few rS S ?

 

 

 

 

 

 

& ш .< д

 

К

1075

1075

1015

1105

 

10Э0

1015

* i 8

75

И

TO

75

20

1025

ft-to

f t -225

Ц

SO

■ y ? : M

1000

 

30

30

(065

 

 

40

«0

юзо

M

75

 

30

20

1070

 

45

 

 

70

 

1005

 

 

1070

1090

 

Таблица 8.6 Построение оптимального плана перевозок. Шаг 2

 

S.-85

B-i« 100

5, = 120

84»120

Bs-m

tf6~!I0

 

«,-eo

tf,-S0

Й ,-225

 

 

,4,“ t50

<0

35

40

90

 

 

 

150

 

 

 

5

75

80

, / / / / / /

1000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

30

 

f

 

P

25

H

65

 

 

 

30

30

 

 

 

 

 

45

 

10

 

 

 

 

 

Ai=150

 

 

 

 

 

 

 

 

1065

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25 20 i 10

35

75

 

 

160

90

30

80

Я

70

40

«о

 

 

155

 

 

 

 

1030

 

 

 

 

 

 

 

 

 

 

 

 

 

 

80

IS

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A -145

45

 

5

35

 

5

30

110

30

40

 

 

75

30

20

1070

 

 

106

 

 

 

 

 

1

 

 

 

 

 

. 45 .

 

 

 

 

 

 

 

 

 

 

 

i .

.

.

* =400

j j j j j j l

 

25

10

1

 

Щ

25

 

80

 

 

20 !

70

 

1005

 

 

I

120

 

Ш

W

(0

 

!

 

я

 

ш

 

 

 

 

 

 

 

 

 

 

 

Kr

1075

j 1075

1015

j

1105

10»

1 10*5

1

 

г

 

 

 

 

1025 1

1070

1090

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

285

Таблица 8.7 Построение оптимального плана перевозок. Шаг 3

Таблица 8.8

286

Таблица 8.9 Построение оптимального плана перевозок. Шаг 5

Таблица 8.10 Построение оптимального плана перевозок. Шаг 6

287

Таблица 8.11 Построение оптимального плана перевозок. Шаг 7

Таблица 8.12

Итоговый оптимальный план перевозок

288

8.2.2.Метод потенциалов для решения транспортной задачи

всетевой форме

Задача 2

Построить оптимальный план перевозки груза на сети (рис. 8.4) от трех станций отправления до девяти станций назначения.

Рис. 8.4. Сеть транспортных узлов

Для данной задачи студент приводит две сети: с начальным пла­ ном перевозок и с оптимальным планом перевозок.

Исходные данные

Размеры груза со станции отправления приведены в табл. 8.13, а прибытия на станцию назначения - в табл. 8.14. Для всех вариантов задан один и тот же полигон (см. рис. 8.4), на котором производится решение задачи. В знаменателе для звеньев 1-7, 1-11, 2-9, 3-10 зада­ ны ограничения пропускной способности.

Методические указания к выполнению задачи 2

1. Решение задачи в сетевой форме начинают с составления на­ чального плана, который не допускает встречных перевозок на уча­ стках заданного полигона. Начальный (или любой допустимый)

289