Не отобразилась форма расчета стоимости? Переходи по ссылке

Не отобразилась форма расчета стоимости? Переходи по ссылке

Дипломная работа на тему «Задача маршрутизации транспорта с временными окнами и ограничением на грузоподъемность транспортных средств рассматривается на примере пермской торговой фирмы ООО ‘Фабрика еды'»

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

Написание диплома за 10 дней

Введение

С точки зрения экономической теории, понятие транспортной логистики рассматривается как планирование оптимальных маршрутов грузоперевозок с минимизацией затрат, что является особенно актуальным для торговых компаний, ежедневно доставляющих продукцию покупателям. В связи с растущей загруженностью автомобильных дорог, водители предприятий часто опаздывают в установленные временные интервалы доставки, что приводит к возвратам продукции, росту неудовлетворенности клиентов, и, в конечном счете, к потерям выручки и убыткам. Существует ряд моделей, применив которые, фирма может существенно сократить транспортные расходы и оптимизировать последовательность посещения контрагентов. Одна из них — классическая задача коммивояжера (КЗК), решением которой является кратчайший замкнутый маршрут водителя, проходящий через всех покупателей по одному разу. Задача была поставлена английскими математиками У. Гамильтоном и Т. Киркманом в середине ХIХ века. В течение последующих десятилетий ученые из разных стран изучали КЗК, разрабатывая новые математические постановки задачи и алгоритмы их решения. Так появились точные и приближенные методы достижения оптимальных маршрутов (Бронштейн, Заико, 2010). К точным методам относятся алгоритм полного перебора и метод ветвей и границ, а к приближенным (позволяют найти маршрут близкий к точному решению задачи) — метод имитации отжига, алгоритм муравьиной колонии и др.

Несмотря на популярность классической задачи коммивояжера, она редко решается на практике, поскольку ее ограничения не учитывают, например, маршрутизацию нескольких машин в автопарке компании, пробки на дорогах и временные интервалы доставки покупателям. Желая приблизить теоретическую задачу к реальности, Данциг и Рамсер, предложили разновидность КЗК, включив в нее ограничение на количество автомобилей и назвав ее — задачей маршрутизации транспорта (ЗМТ) (Danzig, Ramser, 1959). Решение задачи представляет собой оптимальные маршруты для двух и более автомобилей, которые проходят через все точки на карте один раз, возвращаясь в исходную. Данную модель успешно применяют в таких отраслях науки как медицина, машиностроение и программирование и т.д. Подчеркнем, что задача маршрутизации транспорта является обобщением классической задачи коммивояжера и решается похожими алгоритмами.

В настоящем исследовании рассматривается одна из областей применения задачи маршрутизации транспорта — логистика. В работе осуществляется поиск решения ЗМТ на примере пермской торговой фирмы ООО «Фабрика еды», занимающейся оптовой торговлей мясопродуктами. Компания имеет двух водителей, отвечающих за доставку продукции контрагентам. Летом 2015 года у предприятия возникли трудности с развозом продукции покупателям: в связи с увеличивающимся количеством заявок, водители часто опаздывали в желаемое время поставки, установленное покупателями, за что были оштрафованы возвратом продукции в полном размере. Общие непредвиденные транспортные расходы по вине двух водителей составляли от 25 до 35 тыс. в месяц, что послужило причиной поиска их минимизации. В исследовании предполагается, что решив задачу маршрутизации транспорта, торговая компания сократит свои логистические расходы. Отметим, что в работе рассматривается модифицированная ЗМТ, которая учитывает специфику бизнеса: в задачу добавляется ограничение на грузоподъемность транспортных средств, а также два условия, устанавливающие временные окна покупателей. Описываемая задача анализируется на реальных данных торговой компании (рассматривается день с 51 покупателем). ЗМТ решается современными математическими алгоритмами: муравьиным методом и алгоритмом поиска с запретами. Применяемые методы являются эвристическими (приближенными) способами нахождения оптимального решения задачи.

На практике любой путь водителей будет оптимальным, если он удовлетворяет временным интервалам доставки. В исследовании рассматривается как детерминированная, так и стохастическая постановка ЗМТ, процесс пересчета маршрута водителей, а также проводится сценарный анализ пробок в г. Пермь c построением сетевых моделей планирования и управления логистикой. Результаты данной работы могут найти применение не только в исследуемой фирме, но и в различных предприятиях связанных с доставкой продукции и объездом контрагентов, например, в службе инкассации, пожарной и скорой помощи и т.д.

Целью выпускной квалификационной работы является постановка и решение задачи поиска оптимального плана доставки продукции торговой компании.

Задачами работы являются:

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

) Постановка задачи маршрутизации транспорта с временными интервалами доставки продукции и ограничением на грузоподъемность транспортных средств;

) Анализ и сравнение алгоритмов решения различных модификаций задачи маршрутизации транспорта;

) Разработка методологии сведения стохастической задачи маршрутизации транспорта к серии детерминированных задач;

) Сбор данных, необходимых для математического моделирования;

) Применение ряда методов в современных программных пакетах (Умные маршруты, Муравьиная логистика) на данных компании;

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

) Формирование маршрутов развоза продукции водителями фирмы ООО «Фабрика еды» с учетом временных ограничений покупателей и грузоподъемности автомобилей;

) Оценка и анализ полученных результатов.

Объект исследования — математические модели оптимизации логистики торговых предприятий. Предметом выпускной квалификационной работы является применение задачи маршрутизации транспорта с временными окнами и ограничением на грузоподъемность автомобилей на примере торговой фирмы ООО «Фабрика еды». В исследовании использовались следующие методы: анализ специальной и научной литературы, методы математического моделирования и оптимизации.

Выпускная квалификационная работа состоит из введения, четырех разделов, заключения и списка литературы. В первом разделе подробно рассматривается ряд существующих работ по исследуемой теме. Теоретическое обоснование структурно разделено на две части: первая часть посвящена классической задаче коммивояжера, а вторая — задаче маршрутизации транспорта с временными окнами и ограничением на грузоподъемность автомобилей (capacitated vehicle routing problem with time windows — CVRPTW). Отметим, что последняя задача является обобщением первой, следовательно, необходимо провести как совместный, так и сепаратный анализ двух задач. Второй раздел, описывает предметную область (компанию), ограничения и предпосылки, возникающие в ходе моделирования. Рассматривается математическая постановка исследуемой задачи. В третьем разделе подробно описываются используемые данные и алгоритмы решения задачи. Также в третьей части проводится сравнительный анализ основных методов решения по ряду критериев. Большое внимание уделяется муравьиному методу и алгоритму «поиска с запретами». Последний раздел подразумевает обсуждение основных результатов исследования. В данной части работы в программных пакетах реализуются математические методы нахождения оптимальных маршрутов водителей торгового предприятия, строятся сетевые модели планирования и управления, проводится пересчет маршрутов в случае непредвиденных дорожных ситуаций, а также решается классическая задача об оптимальной загрузке автомобиля. Полученные результаты обобщаются в виде выводов, а также даются рекомендации по практическому применению моделей и внедрению их в компании.

маршрутизация транспорт развоз грузоподъемность

1. Теоретическое обоснование

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

Начнем с подробного описания и изучения классической задачи коммивояжера — КЗК («задача странствующего торговца»), поскольку не проанализировав данную задачу, понять более сложную ее модификацию не представляется возможным (Макконнелл, 2004). Отметим, что КЗК представляет собой одну из самых известных математических задач XX века. С экономической точки зрения решением задачи является самый выгодный маршрут, проходящий через указанные города (точки на карте) по одному разу, с возвращением в исходный. В условиях данной задачи традиционно указываются какие-либо критерии оптимальности (наименьшее расстояние или время). Говоря об актуальности исследования задачи, которая появилась в середине ХVIII века, подчеркнем, что с 1959 по 2008 годы около 1 000 научных журналов публиковали статьи про задачу коммивояжера на первых страницах. Более того, научная литература на тему КЗК растет экспоненциально каждый год со скоростью 6% (Baldacci, et al., 2012).

Представим классическую задачу коммивояжера в виде модели на графе. Пусть имеется ориентированный конечный граф  в котором N — множество вершин, А — множество ребер (Хайруллин, 2014). Тогда длина ребра  будет являться его весом и будет связывать вершины . Граф должен быть полным, поскольку необходимо учесть все возможные ребра. Таким образом, решение задачи коммивояжера — это такой маршрут, который является ориентированным полным остовным циклом (проходит через все вершины) минимального веса (сумма всех расстояний между вершинами) в графовой модели. Предполагается, что переменная  является бинарной, то есть принимает значение 1, если ребро графа включено в оптимальный маршрут, проделанный коммивояжером, и 0, если не включено.

Математическая постановка классической задачи коммивояжера выглядит следующим образом:

 

Подробно постановка задачи будет описана в следующем разделе. КЗК относится к числу трансвычислительных (Irnich, 2008): уже при относительно небольшом числе вершин графа (> 67) ее невозможно решить методом полного перебора маршрутов современным компьютером менее чем за несколько миллиардов лет. Задача имеет экспоненциальную алгоритмическую сложность: коммивояжер в каждом городе встает перед выбором следующего из еще не посещенных, а значит, существует

маршрутов для асимметричной (расстояние из города i в город j, не равно обратному) задачи, где n — число городов. В настоящем исследовании будут использоваться именно асимметричные матрицы расстояний и времени, что существенно увеличит алгоритмическую сложность. Классическая задача коммивояжера является NP-полной (Костюк, 2013), что означает невозможность решения каким-либо алгоритмом за полиномиальное время. За нахождение такого метода, американское математическое сообщество предлагает приз в миллион долларов, что говорит об огромном интересе к задаче как с теоретической, так и практической точки зрения.

Приведем наиболее значимые работы о КЗК расположив их в хронологическом порядке. Первым упоминанием о классической задаче коммивояжера (traveling salesman problem — TSP) стали труды английских математиков У. Гамильтона и Т. Киркмана в середине ХIХ века

(Laporte, 1992). В то время авторы не имели необходимых методов для решения задачи и могли сформулировать только вид предполагаемого решения — маршрут коммивояжера, проходящий через каждый город только один раз (необходимо объехать все города без исключения). Описываемый путь торговца получил название гамильтонов цикл. С развитием такого направления науки как комбинаторная оптимизация, исследователи стали придавать задаче коммивояжера теоретический смысл. На основе КЗК были разработаны такие известные методы дискретной оптимизации, как «жадный» алгоритм, метод ветвей и границ, а также ряд современных эвристик (Moon, et al., 2012). Спустя столетие после постановки классической задачи коммивояжера, ее решением занялись американские сотрудники RAND Corporation, в том числе, Д. Фалкерсон. Применив метод ветвей и границ, ученый нашел маршрут для 21 города США и доказал его оптимальность (Toth, Vigo, 2002). В середине 70-х годов ХХ века КЗК активно применялась в биологии, информатике и медицине. Через два десятилетия Д. Апплегейт и У. Кук создали «библиотеку» стандартизированных постановок и модификаций классической задачи коммивояжера — TSPLIB, куда вошли экземпляры, содержащие более 10 000 точек на карте (Applegate, Cook, 1993). В середине 2010 года европейские математики решили КЗК, включающую тысячу вершин. Используя методы кластеризации, удалось доказать, что полученное решение всего на 1% отличается от оптимального (Schneider, et al., 2010).

По данным электронных библиотек, в научной литературе существует более 6 000 исследований, связанных с задачей коммивояжера. Рассмотрим наиболее значимые из них. Кроуз в 1958 году в своей работе продемонстрировал несколько возможных способов упрощения вычислительной сложности задачи (Croes, 1958). Один из них — использование неориентированного графа при построении модели (матрица расстояний — симметрична). Математик обосновал, что симметричную КЗК решать значительно легче по сравнению с асимметричной. Автор показал на практике, что задача с неориентированными дугами графа на практике становится менее реалистичной. Более того, ученый считал, что для упрощения процедуры нахождения оптимального решения, необходимо использовать только линейные условия, добавляя в постановку ограничение элиминирующее подциклы в маршруте коммивояжера (Croes, 1958). Заметим, что преследуя цель нахождения оптимального решения КЗК, на первой итерации необходимо получить допустимое решение. Впоследствии, итеративное улучшение решения, путем ряда перестановок в исходной последовательности городов, приведет к образованию наилучшего решения задачи. Рассмотрим другую немаловажную работу, посвященную решению классической задачи коммивояжера. Лин в 1973 году впервые проанализировал простейший эвристический алгоритм нахождения решения КЗК — «жадный» метод. Автор исследовал вычислительное время задачи и продемонстрировал, что «жадный» алгоритм способен его существенно сократить (Lin, Kernighan, 1973). Подчеркнем, что решение, получаемое в ходе использования «жадного» алгоритма не всегда будет оптимальным, а значит, необходимы дополнительные перестановки в полученной последовательности городов. Переходя к современным статьям о КЗК, отметим, что несмотря на развитие эвристических алгоритмов, многие ученые на протяжении XX века продолжали решать задачу точными методами. Многие ученые были убеждены, что задачу коммивояжера необходимо решать методом ветвей и границ. Заметим, что описываемый алгоритм был популярен в научном сообществе и положил начало большому количеству статей, реализующих его на разном количестве городов (вершин на графе). Метод заключался в разделении множества маршрутов на подмножества, в каждом из которых производился поиск «нижней границы»- кратчайшего маршрута в подмножестве. Наилучший маршрут соединялся из всех «нижних границ» подмножеств (Toth, Vigo, 2002). Осознав необходимость применения задачи на практике, только к концу ХХ века ученые стали решать классическую задачу коммивояжера при помощи приближенных алгоритмов, таких как: алгоритм минимального остовного дерева, локального поиска и др.

Классическая задача коммивояжера решается в логистической практике редко, поскольку ее ограничения не включают, например, моделирование дорожных пробок, маршрутизацию нескольких машин в автопарке и, конечно, временные интервалы доставки покупателям. Поэтому Данциг и Рамсер предложили разновидность КЗК, дополнив задачу ограничением на количество автомобилей и назвав ее задачей маршрутизации транспорта (ЗМТ). В своей статье исследователи поставили вопрос оптимального развоза бензина на автозаправочные станции в городе N, предложив математическую постановку задачи и решив ее эвристическим алгоритмом. Данциг и Рамсер придали ЗМТ следующую форму задачи математического программирования, которая будет анализироваться в следующем разделе работы (Danzig, Ramser, 1959):

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

Отметим, что с математической точки зрения ЗМТ (6) — (11) отличается от классической задачи коммивояжера (1) — (6), наличием индекса k над переменной , который отвечает за номер транспортного средства в задаче. Несмотря на незначительное нововведение в математическую постановку, алгоритмическая сложность задачи существенно увеличивается.

Через несколько лет после публикации статьи Данцига и Рамсера (1959), Кларк и Райт (1964) улучшили подход американских исследователей, предоставив вниманию математического сообщества алгоритм «сбережений», который был легко применим к решению задачи маршрутизации транспорта и на тот момент являлся одним из передовых, поскольку находил оптимальные маршруты за кратчайшее время. Метод «сбережений» начинает поиск оптимального маршрута с первоначального решения задачи, когда каждому транспортному средству ставится в соответствие единственная случайная точка на карте (количество автомобилей равно количеству вершин графа). На каждой последующей итерации алгоритм пытается так уменьшить количество машин (путем объединения двух случайных маршрутов), что «сбережения» (сокращение общей длины или времени маршрута) были максимальными. Кларк и

Райт (1964) впервые предложили проводить оптимизацию маршрутов, решая ЗМТ как последовательно, так и параллельно, доказав, что одновременный выезд автомобилей из начальной точки приводит к более эффективному решению задачи (Clarke, Wright, 1964). Вышеупомянутая статья имела огромный успех, что и привлекло ученых из разных стран к анализу как постановок, так и алгоритмов касающихся стандартной задачи маршрутизации транспорта (обобщение задачи коммивояжера без усложнений и дополнительных ограничений).

Существенное внимание в первых работах о ЗМТ уделялось вычислительной сложности задачи. Так, Кан и Ленстра (1975) определили масштаб и характер задачи маршрутизации транспорта с математической точки зрения, описав новые алгоритмы решения ЗМТ, связывающие точные и эвристические методы. Авторы доказали, что задача маршрутизации транспорта является NP-сложной (как и задача коммивояжера), что говорит о ее неразрешимости за полиномиальное время (Kan, Lenstra, 1975). Тем не менее, многие исследователи конца ХХ века продолжали решать задачу маршрутизации транспорта точными методами.

В своей обзорной статье Лапорт (1992) детально проанализировал ряд работ, связанных с точными алгоритмами решения ЗМТ (полный перебор, метод ветвей и границ и др.) (Laporte, 1992). Ученый выяснил, что огромное количество задач остались нерешенными из-за длительности обработки данных описываемыми методами. Более того, в большинстве статей рассматривались задачи малых размерностей. Таким образом, необходимость углубленного исследования эвристических алгоритмов и новых модифицированных постановок задачи росла с каждым годом.

Отметим, что одной из первых и основных разновидностей классической задачи маршрутизации транспорта стало добавление в математическую постановку ограничения на грузоподъемность автомобилей, которое позволило оптимально распределять грузы между машинами, не превышая суммарную нагрузку на каждое транспортное средство. Как и стандартная ЗМТ, ее модификация с грузоподъемностью в литературе в основном решается классическими методами оптимизации (методом имитации отжига, жадным алгоритмом, и др.). Важно заметить, что задача маршрутизации транспорта с ограничением только на грузоподъемность редко используется на практике, поскольку оптимизация реальных бизнес процессов подразумевает введение и анализ таких ограничений как: временные интервалы поставки, пробки на дорогах в конкретном городе в определенное время дня и др. Такие ограничения неизбежно возникают, когда задача решается для конкретных организаций, которые занимаются транспортировкой людей или доставкой товаров (продукции). К примеру, Дерошье в своей статье ввел понятие «жестких» и «мягких» временных интервалов доставки грузов, что существенно изменило взгляд на ЗМТ, образовав новую математическую постановку задачи — с временными окнами и ограничением на грузоподъемность (capacitated vehicle routing problem with time windows — CVRPTW) (Desrochers, et al., 1992). Основным отличием «жестких» временных окон от «мягких» является возможность отгрузки продукции покупателю в случае нарушения временного интервала: если водитель опаздывает к покупателю с «жесткими» временными окнами — продукция будет непринята в полном размере (в случае с «мягкими» интервалами водитель получит штраф, но сможет произвести отгрузку). Более того, Дерошье в 1992 году обосновал, что решение задачи будет оптимальным, если автомобиль не имеет возможности приехать к покупателю позже заданной правой границы «жесткого» временного окна, в отличие от случая с «мягкими» временными интервалами.

В конце ХХ века большинство научных работ посвящались именно задаче маршрутизации транспорта с «жесткими» временными окнами, поскольку она в большей степени применяется на практике, начиная от деятельности инкассаторов и почтовых работников до маршрутизации школьных автобусов (Irnich, 2008). Математик выяснил, что на современных компьютерах, точные методы решают вышеописанную задачу с 50 — 60 точками на карте более чем за 50 часов, что недопустимо для оптимизации реальных логистических процессов. Так, Соломон в 1987 году рассмотрел первый эвристический алгоритм, позволяющий решить ЗМТ с временными интервалами и грузоподъемностью транспортных средств за существенно меньшее время, по сравнению с точными методами (Solomon, 1987). Эвристика была основана на алгоритме сбережений Кларка и Райта (1964) для классической задачи маршрутизации транспорта и улучшала текущее решение на каждой итерации, позволяя решать задачу на ориентированном графе (матрица расстояний между каждой точкой на карте асимметрична). Более того, Соломону удалось ускорить процесс решения задачи почти в 2 раза (Solomon, 1987).

Помимо эвристических методов, в современных статьях появилась новая категория алгоритмов — метаэвристики (обобщенные способы решения задачи коммивояжера) (Пантелеев, Метлицкая, Алешина, 2013). Одним из наиболее популярных метаэвристических методов является алгоритм муравьиной колонии, который будет использоваться в настоящем исследовании в разделе Данные и методология. Рассмотрим его подробно. Впервые в литературе М. Дориго сформулировал и реализовал алгоритм муравьиной колонии, который эмпирически решает ЗМТ быстрее, в сравнении с другими существовавшими и доступными алгоритмами конца ХХ века (Dorigo, 1999). В работах отмечается, что муравьиный алгоритм эффективен для количества вершин больше 20, то есть позволяет учесть временные ограничения в маршруте у наибольшего количества покупателей, в отличие от других методов оптимизации (гнездовой метод Монте-Карло, модифицированный генетический алгоритм и др.) (Курейчик, Мартынов, 2014). Более того, алгоритм муравьиной колонии меньше подвержен неоптимальным начальным решениям по сравнению с обычными эвристическими алгоритмами. Маршрут строится случайным образом и колония при выборе каждой следующей вершины опирается на память всех предшествующих муравьев, а не нескольких, в отличие от алгоритмов, построенных по принципу марковских процессов (текущее состояние системы зависит только от предыдущего). Итак, алгоритм муравьиной колонии основан на поведении множества искусственных муравьев, которые кооперируются для нахождения решения задачи путем обмена феромонами, находящимися на дугах графа (дороги в городе). Как известно, муравьи способны находить наилучший маршрут от источника еды до муравейника, оставляя след из феромона, по которому будут двигаться муравьи-последователи. Поскольку на самом выгодном маршруте феромона будет скапливаться больше, чем на других, через определенный момент времени, вся колония предпочтет его. Данный процесс получил название коллективное обучение (Dorigo, 1999). Необходимо отметить, что алгоритм муравьиной колонии учитывает как пространственные, так и временные характеристики задачи. С его помощью планируется найти решение, близкое к оптимальному.

Говоря о работах, связанных с подробным рассмотрением получаемых решений ЗМТ, отметим, что многие авторы считают недостаточным одного лишь решения задачи без последующего анализа. Экономисты утверждают, что для оптимизации логистики также необходимо обращать внимание на сам процесс формирования заказов, который может быть не оптимальным. Описываемая ситуация моделируется в виде задачи о ранце 0-1 (0-1 knapsack problem) (Toth, Vigo, 2002), которая имеет следующую математическую постановку:

 

где:  — цена j-го заказа;

— вес j-го заказа;

— переменная, отвечающая за нахождение j-го заказа в автомобиле;

— суммарная грузоподъемность транспортного средства.

Таким образом, ограничение (12) является целевой функцией классической задачи о ранце и представляет собой дневную выручку компании, которую необходимо максимизировать. Ограничение (13) показывает, что вес всех заявок в автомобиле не должен превышать его грузоподъемность. Последнее ограничение (14) — говорит о том, что переменная принимает только два значения — 1, если заказ загружен в машину, 0 — в противном случае.

В научных статьях существует огромное количество модификаций и алгоритмов решения данной задачи (Иглин, 2003). Существуют задачи о многомерном ранце, целочисленная задача о рюкзаке, задача о рюкзаке с мультивыбором и др. Заметим, что задача о ранце (как и задача коммивояжера) является NP-сложной и не решается ни одним полиномиальным алгоритмом за мыслимое время (Макконнелл, 2004). В выпускной квалификационной работе частично затрагивается процесс фильтрации покупателей в случае возникновения излишнего спроса на продукцию фирмы. Подчеркнем, что решение задачи о рюкзаке не является целью настоящего исследования и будет рассмотрено в дополнение к анализу решения задачи маршрутизации транспорта с ограничением на грузоподъемность и временными окнами покупателей в разделе Описание результатов.

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

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

2. Постановка исследовательской проблемы

В выпускной квалификационной работе задача маршрутизации транспорта с временными окнами и ограничением на грузоподъемность транспортных средств рассматривается на примере пермской торговой фирмы ООО «Фабрика еды». Компания занимается оптовой торговлей скоропортящейся продукцией: мясопродуктами, курицей, а также полуфабрикатами. В июне 2015 года в фирме появился второй водитель, который также как и первый, отвечает за развоз продукции к покупателям. Фактическим адресом фирмы является ул. Академика Курчатова 9, что находится в 12 минутах езды до центра города Пермь и, следовательно, накладывает определенные сложности на быстроту доставки продукции контрагентам. Прием заявок осуществляется офисом, который находится по вышеупомянутому адресу, за день до отгрузки (доставки). В заявке покупателей указывается желаемое время отгрузки с 11:00 до 20:00 (время работы водителя). Обозначаемый в заявке временной интервал должен быть не меньше 30 минут (для того, чтобы водитель успел доставить продукцию в отдаленные районы). Цех предприятия начинает работу с 6:00 утра каждый день, для того чтобы максимально удовлетворить все текущие заявки покупателей. Среднее количество контрагентов, ожидающих поставку в случайный день, представлено на графике:

Рис. 1. Среднее количество покупателей на одного водителя

В определенный момент времени, а именно, осенью 2015 года, у компании возникли сложности с доставкой продукции покупателям: в связи с увеличивающимся количеством заявок один из двух водителей ездил не оптимально, не попадая в желаемое время отгрузки, за что покупатели штрафовали предприятие, делая возврат в полном объеме. Несмотря на то, что некоторые контрагенты были лояльны к времени отгрузки и опозданиям, ООО «Фабрика еды» каждый месяц несла потери в виде возвращенной (не принятой) продукции на 25 — 35 тыс. рублей по вине как первого, так и второго водителя.

Описывая типичный рабочий день водителей, необходимо отметить, что они развозят продукцию на машинах ИЖ-2717, которые оборудованы специальными рефрижераторами для скоропортящихся продуктов. В среднем один из водителей не успевал к одному контрагенту в неделю (потеря продукции на 4 — 5 тыс. рублей). Из причин, на которые ссылались водители, выделим одну наиболее значимую: «Производя доставку более приоритетным покупателям в центре города, практически невозможно успеть вовремя в отдаленные районы, которые также желают получить продукцию в определенный срок». Таким образом, предприятие хотело, по возможности, уменьшить свои транспортные расходы, определив строгую последовательность посещения заказчиков двумя водителями.

Данный процесс развоза продукции для одного водителя можно смоделировать в виде задачи коммивояжера, представленной системой ограничений (1) — (6). Итак, в один из дней работы исследуемой фирмы имеется N покупателей. Водитель, выезжает один раз в начале рабочего дня с предприятия в 11:00 (предварительно загрузив продукцию). Он доставляет заказы всем контрагентам (соответственно, побывав у каждого один раз) и возвращается после максимального удовлетворения всех заявок на предприятие (не обязательно в 20:00). В задаче будет находиться наименьшее время, а не расстояние, поскольку у водителя основной целью является успеть во временные окна как можно к большому числу покупателей (увеличение или уменьшение расстояния на несколько километров не приведет к существенным потерям для фирмы, по сравнению с возвратом продукции). Если водитель возвращается в фирму раньше окончания рабочего дня, он может заниматься, например, доставкой продукции от поставщиков. Таким образом, задача водителя заключается в определении последовательности объезда покупателей, которая будет минимизировать его путь, а также удовлетворять временному интервалу, установленному каждым покупателем (Lorini, 2011). Заметим, что в исследовании для нас представляет интерес не только задача коммивояжера, но, в большей степени, задача маршрутизации транспорта (6) — (11). Вторая является обобщением первой задачи. В работе используется модифицированная постановка ЗМТ — задача маршрутизации транспорта с временными интервалами и ограничением на грузоподъемность транспортных средств:

 

Переменные задачи:  — бинарная неотрицательная переменная, принимающая значение 1, если водитель под номером k едет от покупателя i к контрагенту j, и значение 0 в остальных случаях;  — показывает фактическое время прибытия водителей к покупателям;  — неотрицательная переменная, предотвращающая распадение гамильтонова цикла на подциклы. Отметим, что две последние переменные не зависят от номера водителя (k), поскольку данные переменные непосредственно связаны с покупателями, а не с водителями фирмы (например, в задаче не задается какой именно водитель приедет к пятому покупателю в 12:00, а учитывается сам факт приезда одного из водителей в конкретное время). ЗМТ с временными окнами и грузоподъемностью автомобилей — (6) — (11), (15) — (17) является обобщением КЗК (1) — (6) и ЗМТ (6) — (11), и именно поэтому будет подробно рассмотрена ниже. (7) — целевая функция задачи — означает, что суммарное время водителей через всех покупателей должно быть наименьшим. (15) — обозначает, что суммарный спрос на продукцию (в кг) не должен превышать максимальную грузоподъемность каждого транспортного средства. (8) и (9) показывают, что оба водителя после отгрузки продукции покидают каждого покупателя только один раз и приезжают к другому контрагенту для поставки только один раз. Поскольку, решением задачи коммивояжера являются замкнутые маршруты водителей, так называемые гамильтоновы циклы, ограничение (10) означает невозможность распадения данного цикла на подциклы (объезд всего нескольких покупателей из всей заявки). Ограничение (11) указывает на то, что переменные задачи  являются бинарными и неотрицательными, а (6) означает, принадлежность переменных  к множеству вещественных чисел. Большой интерес в данной постановке представляют два последних ограничения. В повседневной деятельности, предприятие заключает договоры поставки с новыми контрагентами, которые склонны к оппортунистическому поведению, а именно к сокращению желаемого времени доставки продукции (которое прописано в контракте) в новых ежедневных заявках. Торговому предприятию зачастую приходится подстраиваться под своих покупателей для того, чтобы их сохранить, а следовательно, встает необходимость удовлетворения заявок контрагентов именно в желаемое для них время, что ведет к модификации традиционной логистической задачи ЗМТ (Moon, et al., 2012). С математической точки зрения, в задачу маршрутизации транспорта необходимо добавить ряд дополнительных предпосылок, основной из которых будет ограничение на временные интервалы, в которые покупатели ожидают доставку товара. Важным свойством задачи маршрутизации транспорта с временными интервалами и ограничением на грузоподъемность транспортных средств является то, что минимальные по времени маршруты начинаются из ООО «Фабрика еды» в 11:00, в дальнейшем проходя через все вершины (покупателей) только один раз, в течение соответствующих временных интервалов и заканчиваются в точке начала маршрута. Исходя из этого, каждому покупателюв дневной заявке будет поставлено в соответствие «жесткое» временное окно , где границы интервала будут обозначать начало ожидаемого времени поставки и его окончание соответственно. Временные интервалы покупателей должны устанавливаться во время работы водителей, которое обозначим — [0, ]. Отметим, что как опытному, так и новому водителю, необходимо определенное время  для того, чтобы обслужить i-го покупателя. В модифицированной задаче, у каждого ребра весом будет являться время проезда от одного покупателя к другому (). Таким образом, в задачу (6) — (11), наряду с ограничением на грузоподъемность, добавятся два новых ограничения (16) и (17). (16) — задает временные окна у каждого покупателя, в которые водителю необходимо отгрузить продукцию. Ограничение (17) показывает временную взаимосвязь между последовательными покупателями в маршруте. Рассмотрим двух контрагентов к которым водитель приезжает напрямую (). Очевидно, что время приезда к покупателю j будет больше либо равно сумме времени приезда к покупателю i, времени разгрузки у контрагента i, а также времени проезда от i к j: . Однако такое ограничение справедливо только для тех пар клиентов i и j, которых хотя бы один из водителей посещает последовательно. Поэтому описываемое ограничение (17) становится нелинейным, что накладывает определенные сложности в выборе алгоритма решения данной логистической проблемы. Таким образом, решив детерминированную задачу (6) — (11), (15) — (17), фирма достигнет своей логистической цели — успеть вовремя ко всем покупателям в их временные интервалы, получив минимальные по времени маршруты водителей.

В работе, помимо описываемой детерминированной задачи, затрагивается стохастическая постановка ЗМТ, в которой целевая функция направлена на минимизацию риска возникновения потерь продукции. Предположим, что решив задачу (6) — (11), (15) — (17), водитель приезжает к пятому покупателю в 11:58, а правая граница временного окна — 12:00. Очевидно, что в случае, например, затора на дороге или задержки у предыдущего покупателя, водитель не успеет доставить продукцию в желаемый интервал отгрузки. Для фирмы более осмысленной представляется задача минимизации математического ожидания потерь, связанных с доставкой продукции позже указанного срока:

 

где:  — вероятность недоставки продукции i-му покупателю;

— выручка от i-го контрагента.

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

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

 

Рис. 2. Вероятность недоставки продукции покупателям

В исследовании подразумевается, что вероятность непринятия продукции i-ым покупателем при  принимает значение 0,2. Данное значение было получено на основе статистического анализа неотгруженных заказов.

Таким образом, ограничение (18) можно представить в виде:

 

где: k и c — параметры линейной убывающей функции.

После изменения функционала задачи, возникает стохастическая постановка — (6), (8) — (11), (15) — (17), (18.1), решение которой, с экономической точки зрения, позволит снизить риск недоставки продукции путем наиболее раннего приезда водителей к покупателям. Очевидно, что описываемая модификация ЗМТ с временными окнами и грузоподъемностью автомобилей неединственная. Данную задачу можно преобразовать под любые особенности бизнеса, например, введя ограничения на случайный размер заказа или на обслуживание одного покупателя несколькими транспортными средствами. Оказалось, что исследуемая компания не ставит цель наиболее раннего приезда ко всем контрагентам (есть покупатели с временными окнами длинной в 6 часов и фирме без разницы: приедет водитель за 2 или за 3 часа до конца интервала). В ходе анализа расходных накладных и личных интервью с водителями ООО «Фабрика еды» было установлено, что минимальная разность (позволяющая максимально снизить риск недоставки) между правой границей временного окна у каждого покупателя и приездом водителя  должна составлять 10 минут, что приведет к трансформации линейной убывающей функции в кусочно-линейную невозрастающую:

Рис. 3. Вероятность недоставки продукции ООО «Фабрика еды»

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

 

Таким образом, в выпускной квалификационной работе будет анализироваться детерминированная ЗМТ (6) — (11), (15) — (17), решение которой (кратчайшие по времени маршруты для двух водителей, удовлетворяющие всем временным окнам покупателей) путем изменения параметра надежности (временного резерва) будет приближено к решению стохастической задачи. Оптимальные по времени маршруты водителей будут удовлетворять специфике реального бизнеса и минимизировать транспортные расходы компании при математическом ожидании риска недоставки продукции равном нулю. Отметим, что, несмотря на относительную простоту постановки задачи (6) — (11), (15) — (17), решать ее достаточно сложно, поскольку необходимо не только сформировать определенный набор данных с учетом различных характеристик предприятия, но и подобрать для них алгоритмы нахождения оптимального решения (Rodriguez, Ruiz, 2012).

. Данные и методология

Информационной базой исследования являются данные бухгалтерского и управленческого учета фирмы ООО «Фабрика еды». На первой стадии их сбора анализировались все дневные заказы покупателей начиная с лета 2015 года (когда у фирмы появился второй водитель) по настоящее время. На основании расходных накладных и счет фактур, которые формировались в программе 1С: Предприятие 8.2, было выяснено, что у фирмы имеются ежемесячные «проблемные» поставки клиентам. Количество покупателей в описываемых заявках варьируется от 20 до 32 на одного водителя и зависит от времени года (осенью и зимой количество клиентов больше, чем летом). В настоящем исследовании рассматривается один из «проблемных» дней с 52 контрагентами (включая фирму-поставщика). Заявки данных покупателей выполнялись 18 сентября 2015 года. Анализируемый день был проблемным для предприятия, поскольку компания несла убытки по вине одного из водителей в размере 14 560 рублей. Фактические маршруты водителей в рассматриваемом дне представлены в таблицах 1 и 2, а их реальные адреса доставки в Приложении 13.

Таблица 1

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диплом

Фактический маршрут 1-го водителя в исследуемый день*

* — в таблицах представлены реальные данные компании ООО «Фабрика Еды», поскольку решение ЗМТ рассматривается на примере описываемой фирмы.

Таблица 2

Фактический маршрут 2-го водителя в исследуемый день

 

В таблицах отражен 51 контрагент (с помощью видеорегистраторов производилась фиксация реального времени проезда от одного клиента к другому, а также время обслуживания каждого покупателя), которым водитель должен доставить заказ с 11:00 до 20:00 в установленное временное окно. В последнем столбце таблицы приведено фактическое время, которое водители тратили на разгрузку заказов у каждого клиента, а также на оформление необходимых документов. Отметим, что в рассматриваемом дне деятельности фирмы первый водитель не успел к контрагенту № 19, который находится в г. Добрянка (примерно 70 км от центра города). Сумма продукции для данного покупателя в исследуемой заявке составляла 14 560 рублей. Водитель, не успев к данному покупателю на 16 минут, был оштрафован возвратом продукции в полном размере. Компания пыталась сделать данный заказ приоритетным, что привело к опозданиям к последующим покупателям в центре города. Таким образом, фирма сталкивается со следующими проблемами, которые будут решены в ходе анализа ЗМТ:

а) Формирование последовательности обслуживания покупателей;

б) Сокращение потерь продукции, как в исследуемой заявке, так и в ряде других;

в) Оптимальное распределение загруженности водителей.

Таким образом, для математического решения логистической задачи на следующей итерации сбора данных формировались матрицы расстояний и времени между 52 точками на карте г. Пермь (покупателями и фирмой). Каждое из 2 704 расстояний/времени проезда между контрагентами рассчитывалось в современных онлайн-картах Google Maps, в которых имеется возможность учесть дорожную ситуацию (пробки, заторы), а также специфику автомобильного движения в г. Пермь (большое количество дорог с односторонним движением). С помощью этих сервисов были собраны расстояния (в км) и время проезда (в мин) между покупателями, максимально близкие к реальным. Таким образом, удалось получить асимметричные матрицы расстояний и времени между всеми покупателями.

Итак, после сбора данных для исследования, необходимо решить задачу маршрутизации транспорта с временными окнами покупателей и ограничением на грузоподъемность транспортных средств математическими методами, выбор которых происходил не случайным образом. На первом этапе анализа проверяется необходимость найма второго водителя фирмой — решается задача коммивояжера, в случае, если бы у фирмы имелся один водитель и он бы развозил 51 заявку индивидуально. На следующем этапе, фактические маршруты двух водителей оптимизируются и решение сравнивается с реальным. Далее задача разделяется на две части и каждая решается по отдельности. Более того, в выпускной квалификационной работе предлагается методология сведения стохастической задачи маршрутизации транспорта к серии детерминированных ЗМТ с использованием экспертной информации о вероятности опозданий к покупателям в зависимости от имеющихся временных запасов. В исследовании рассматривается решение ряда детерминированных задач с разными параметрами надежности (. Помимо одинакового параметра напряженности для всех контрагентов (10 минут), рассматривается вариант установления разных временных резервов для покупателей с различным весом заказов (см. Описание результатов, таблица 4). В конечном счете, целая задача маршрутизации решается с учетом трех сценариев дорожной ситуации в г. Пермь: «благоприятная», «умеренная», «неблагоприятная» и полученные результаты обобщаются в виде сетевой модели планирования и управления (PERT). Отметим, что каждый из сценариев математически представляет собой оптимальный план доставки продукции покупателям, скорректированный на время проезда утром, днем и вечером, путем увеличения или уменьшения времени проезда от одного покупателя к другому (. В ходе достижения оптимального решения ЗМТ используются два основных метода: муравьиный алгоритм и метод поиска с запретами. Проведем анализ каждого, предварительно рассмотрев различные классы алгоритмов с точки зрения их практического применения.

Точные методы решения логистических задач (полный перебор, метод ветвей и границ) позволяют найти глобальный оптимум — наилучший из всех возможных маршрутов и широко используются для решения задач маршрутизации без дополнительных ограничений. Доказано, что с увеличением количества вершин в задаче время ее решения точными методами увеличивается экспоненциально (Toth, Vigo, 2002). Вторая группа алгоритмов — эвристические методы. В отличие от точных методов решения задачи, они позволяют найти приближенные маршруты водителей (локальные оптимумы), которые, с теоретической точки зрения, не всегда являются наилучшими. На практике любая последовательность объезда покупателей, улучшающая текущий маршрут и удовлетворяющая временным интервалам, будет приемлемой для предприятия. В работе применяется один метод из данной категории: алгоритм муравьиной колонии, который рассматривался в теоретическом обосновании исследования. В данном разделе работы опишем алгоритм формально. Каждой дуге графа будет поставлено в соответствие два веса: мера привлекательности — , называемая феромоном, и мера актуальности — . В задаче маршрутизации транспорта матрицы расстояния и времени являются асимметричными, а следовательно, . Суть алгоритма заключается в том, что в самом начале процедуры  муравьев находится в точке начала маршрута. Каждый муравей выбирает следующую вершину графа исходя из вероятностно-пропорционального правила (Dorigo, 1999):

 

где:  — вероятность перехода муравья под номером k из вершины i в вершину j;

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диплом

, алгоритм вырождается в «жадный», если , муравьи реагируют только на количество феромона).

Знаменатель формулы позволяет вероятности находиться в пределах единицы ( — множество еще не посещенных муравьями вершин, смежных с i. Таким образом, муравьи предпочитают выбирать каждую последующую вершину, исходя только из меры актуальности и количества феромона. После расчета вероятности перехода из одной вершины в другую, в алгоритме используется генератор случайных чисел от 0 до 100. Если найденная вероятность больше, чем случайное число, муравей передвигается в данную вершину, если нет, то рассматривается следующая вершина из еще не использованных в маршруте. Когда муравей проходит каждую последующую точку, количество феромона уменьшается. После прохождения всего маршрута одним из муравьев используется правило обновления феромона:

 

где: — случайное число, отвечающее за интенсивность испарения феромона (0<<1);- момент времени.

Алгоритм муравьиной колонии находит решение либо через фиксированное количество итераций, либо через определенное количество времени. Важным моментом в решении задачи муравьиным алгоритмом является связь расстояния и времени, которая будет моделироваться через параметр  — актуальность посещения каждой конкретной вершины (разница между правой границей временного окна i-го покупателя и текущим моментом времени) (Лебедев, 2012):

 

где: , является шириной временного интервала до начала обслуживания покупателя j;

— текущий момент прохождения маршрута, — левая и правая граница временного окна у покупателя под номером j;

— время проезда от контрагента i к j;

— коэффициент штрафующий за не включение j-ой вершины в маршрут (с экономической точки зрения показывает приоритет обслуживания j-го покупателя).

Заметим, что муравьиным алгоритмом невозможно решить одновременно или последовательно целую задачу маршрутизации транспорта в современных программных пакетах: необходимо вершины графа разбить каким-то образом между транспортными средствами. Так, в работе будет использован простейший алгоритм разбивки (кластеризации) данных, который реализован в программном пакете «Муравьиная логистика». С точки зрения экономической теории существует огромное количество способов разбить один массив данных более чем на две части. В исследовании будет представлена кластеризация вершин графа (точек на карте г. Пермь) при помощи алгоритма минимального остовного дерева (Schneider, et al., 2010). Принцип метода основан на «жадном» алгоритме: между всеми покупателями (за исключением фирмы-поставщика) будет найдены наименьшие расстояния. Таким образом, удастся получить путь с началом в самое удаленной точке на карте и концом в противоположной части г. Пермь. Затем, установив ограничение на «центры масс» кластеров (в исследовании необходимо разбить ЗМТ на две части) не менее чем 25 точек, удастся распределить массив покупателей между двумя водителями, назначив каждому случайный кластер. Отметим, что кластеры будут сформированы на основе расстояния между покупателями, не учитывая временной критерий. С практической точки зрения, любое торговое предприятие может самостоятельно выбрать необходимый критерий кластеризации. Только после распределения покупателей между двумя водителями каким-либо математическим методом, можно применить алгоритм муравьиной колонии. В зависимости от оптимальности полученного решения, имеет смысл решить задачу без разбивки ее на кластеры, использовав алгоритм поиска с запретами, который может учесть вышеизложенный факт. Проанализируем данный метод решения ЗМТ в деталях.

Используемый в исследовании алгоритм одновременного решения задачи (без деления на отдельные задачи коммивояжера для каждого водителя), является метаэвристическим и основан на алгоритме имитации отжига, а также механизмах функционирования краткосрочной памяти человека. На данный момент, эмпирически доказано, что поиск с запретами является наиболее эффективным методом решения описываемой задачи, а именно, он позволяет найти решение за наименьшее время (в сравнении с другими метаэвристиками) (Wiener, et al., 2009). Поиск с запретами начинает решать задачу с первоначального маршрута (, найденного, к примеру, «жадным» методом. На каждой итерации (t), алгоритм (путем случайной перестановки последовательности вершин графа) улучшает решение (), минимизируя целевую функцию. Итерации продолжаются до тех пор, пока не будет выполнен критерий остановки (конечное количество этапов алгоритма, либо заданное время). Метод поиска с запретами очень похож на алгоритм имитации отжига, следовательно необходимо описать последний для выявления расхождений.

В методе имитации отжига очередная последовательность вершин (точек на карте) также формируется случайным образом, путем незначительного изменения предыдущего маршрута. Один из самых распространенных в литературе вариантов изменения — перестановка двух случайно выбранных вершин графа. Если маршрут получился лучше всех предыдущих итераций, то он рассматривается как следующий для улучшения. Если же маршрут получается хуже, то один из вариантов — не использовать его. Придерживаясь описываемой логики, решение задачи может попасть в один из локальных экстремумов (минимумов), которых достаточно много в задаче маршрутизации транспорта. Именно исходя из данных предположений маршрут не отвергается (Лопатин, 2005). Вероятность его принятия рассчитывается по формуле:

 

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диплом

где:  — неотрицательная разность между тестируемым и ранее полученным решением;- постоянно убывающий параметр, снижение которого происходит, например, по закону полураспада (при коэффициенте «альфа» равном 0,5): . В случае ЗМТ, чем больше параметр T, тем выше вероятность, того, что алгоритм преодолеет локальный минимум. По мере уменьшения данного параметра, алгоритм достигает конечного решения. Метод имитации отжига эффективен для задач с количеством вершин более 15, поскольку он находит приблизительное решение всего с 2 % отклонений от оптимального за несколько минут. Основным его недостатком является невозможность учета дополнительных ограничений в задаче, например, времени, пробок, грузоподъемности и др.

Главным отличием алгоритма поиска с запретами от метода имитации отжига является то, что последним нельзя решить ЗМТ из-за слишком большого количества локальных минимумов, в то время как алгоритм поиска с запретами обладает краткосрочной памятью и имеет свойство «выходить» из локальных оптимумов задачи с помощью «листа табу», в который заносятся маршруты, сформированные на предыдущих итерациях после достижения одного из локальных экстремумов (Чеблоков, Ченцов, 2005). Таким образом, у алгоритма появляется возможность «не возвращаться» в предыдущие локальные оптимумы, а итерационно двигаться к следующему, вне зависимости от того, лучше он или нет. Наличие краткосрочной памяти позволяет алгоритму поиска с запретами запоминать все найденные экстремумы и на последней итерации выбирать тот, в котором значение целевой функции задачи маршрутизации транспорта минимальное. Описываемый алгоритм позволяет решить ЗМТ с различными ограничениями быстро и эффективно, по сравнению с другими методами нахождения оптимального маршрута. Более того, на данный момент, поиск с запретами — один из немногих алгоритмов, который имеет программное обеспечение и привязан к геоинформационным системам (ГИС), что позволяет решить задачу еще эффективнее. Таким образом, после обработки данных и выявления наиболее приемлемых как с теоретической, так и практической/программной точки зрения алгоритмов решения ЗМТ, рассмотрим вышеизложенную методологию в разрезе результатов исследования.

. Описание результатов

Применив алгоритмы на реальных данных предприятия ООО «Фабрика еды», удалось решить задачу маршрутизации транспорта с временными окнами покупателей и ограничением на грузоподъемность автомобилей двумя способами: разделить задачу, распределив контрагентов между водителями кластерным методом, а также решить целую задачу одновременно. Математические алгоритмы реализовывались в современных программных пакетах, таких как «Муравьиная логистика» и «Умные маршруты». Опишем подробно получившиеся результаты.

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

ООО «Фабрика еды» о найме второго водителя, поскольку одно транспортное средство фактически не могло развезти 51 заявку с 11:00 до 20:00 (см. Приложение 1). Муравьиный алгоритм, нашел решение задачи за

3 минуты. Несмотря на быстроту получения маршрута, он включал 19 опозданий к покупателям, что по факту выразилось бы в значительных денежных потерях. Отметим, что водитель доставил бы продукцию только к утру следующего дня после первого выезда. Так, возникла необходимость рассмотреть и оптимизировать логистику не для одного, а для двух водителей фирмы.

На следующем этапе были по отдельности оптимизированы реальные маршруты водителей, полученные при помощи видеорегистраторов транспортных средств и личных интервью (см. Приложения 2, 3). Данный вариант оптимизации логистики компании проводился для того, чтобы понять, действительно ли водители ездили не лучшим образом в исследуемом дне. Оказалось, что муравьиный алгоритм смог найти оптимальные решения для двух задач коммивояжера, улучшив фактические дневные планы водителей. Так, первый водитель избежал опоздания на 16 минут в г. Добрянка, что позволило сократить потери фирмы как в денежном, так и в количественном выражении. Отметим, что у водителя получился достаточно «плотный» график объезда покупателей (без учета различной дорожной ситуации, водитель превысил свое рабочее время на несколько минут). Для второго водителя удалось уменьшить простои (ожидания) перед началом временных окон некоторых покупателей, что также говорит о более рациональном решении задачи математическими методами, чем ее решение, основанное на личном опыте водителей.

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

В исследовании удалось оптимизировать логистику ООО «Фабрика еды» путем одновременного решения задачи маршрутизации транспорта с временными окнами и ограничением на грузоподъемность транспортных средств алгоритмом поиска с запретами. Результаты метода представлены в Приложениях 6, 10. Анализируя данные таблицы, заметим, что алгоритму поиска с запретами было задано конечное время итераций — 10 минут, за которое удалось получить оптимальное решение задачи. У первого водителя аннулировалось опоздание в г. Добрянка, а у второго общий дневной простой сократился до 5 минут, а также увеличилась нагрузка на автомобиль (стало больше покупателей в центре г. Пермь). Таким образом, в работе удалось решить детерминированную ЗМТ методом поиска с запретами, оптимизировав в первом приближении логистику фирмы ООО «Фабрика еды». Более того, после получения наилучших маршрутов, необходимо приблизить их к реальной недетерминированной жизни, снижая риск недоставки продукции (как было отмечено ранее, водителям достаточно иметь временной резерв минимум в 10 минут до конца временного окна у каждого покупателя). В данном случае, одним из способов достижения оптимальных маршрутов водителей является сокращение всех временных интервалов покупателей на 10 минут и повторное решение задачи. Проведя анализ первоначального решения ЗМТ (см. Приложения 6, 10), было выяснено, что только первый водитель имеет четырех покупателей, к которым приезжает меньше чем за 10 минут до конца временного окна. Второй водитель, несмотря на большее количество заявок, удовлетворяет описываемому ограничению, в виду значительного количества широких временных интервалов (минимальный резерв времени — 16 минут).

Подчеркнем, что в работе на примере водителя № 1 решается серия детерминированных задач, потому что только он опаздывал к клиентам в исходном решении (см. Приложение 6). С точки зрения математической теории, локальная оптимизация (решение серии задач только для первого водителя, не меняя оптимальный план для второго) сужает область допустимых значений, что может привести к отсутствию решения. Торговой фирме на первом этапе предлагается решить ЗМТ для двух водителей методом поиска с запретами. Затем, если решение найдено и удовлетворяет всем требованиям бизнеса, оно может считаться наилучшим. Если решение ЗМТ не удовлетворяет требованию временного резерва равного 10 минутам у одного из водителей, то оптимизация будет проводится локально, поскольку второй водитель уже может ехать по кратчайшему маршруту без дополнительных модификаций пути. Если решение для одного из водителей не может быть найдено, фирме необходимо перейти к глобальной оптимизации — сократить временные окна как у первого, так и у второго водителя одновременно. В работе помимо подробного рассмотрения локальной оптимизации задачи был учтен вариант сокращения временных окон покупателей для двух водителей одновременно. В описываемом случае, целевая функция ЗМТ возросла на 34 минуты для первого водителя, по сравнению с локальной оптимизацией маршрута (см. таблицу 3). Заметим, что увеличение произошло из-за изменившейся последовательности покупателей и их количества (у водителя № 1 стало на 3 покупателя больше по сравнению с локальной оптимизацией). Таким образом, фирме сначала необходимо попытаться оптимизировать маршрут одного из водителей (если второй маршрут оптимален) и только потом переходить к глобальной оптимизации, если первый вариант покажет свою несостоятельность. Решение детерминированной задачи без временных резервов алгоритмом поиска с запретами позволит изначально разделить покупателей на две группы между водителями, следовательно, возможна как локальная так и глобальная оптимизация задачи.

Итак, в ходе решения серии детерминированных задач с различными параметрами надежности ( было установлено, что максимальный временной запас у первого водителя может составлять не 10, а 12 минут:

Таблица 3

Оптимальный маршрут водителя № 1 (с учетом временного резерва — 12 минут)

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диплом

 

Отметим, что на первом этапе все временные окна были сокращены на 10 минут и найдено решение задачи. Затем, рассматривался максимально возможный параметр напряженности, который на третьей итерации составил 12 минут (именно данное решение анализируется далее в работе). Алгоритм за разумное время не выдает решения задачи с временным резервом в 13 минут и более, которое бы удовлетворяло всем ограничениям (например, водитель начинал опаздывать к различным покупателям, что является недопустимым для любой торговой фирмы). Подчеркнем, что маршрут водителя будет оптимальным только в том случае, если удовлетворяет всем временным интервалам покупателей. Используя описываемую процедуру, торговая компания сможет снизить вероятность недоставки продукции до минимума. Более того, компания может рассматривать установление разных временных резервов для покупателей (например, чем больше сумма заказа, тем больше резерв). Рассмотрим описываемый вариант решения задачи подробнее. Разделим всех клиентов, которым водитель № 1 доставлял продукцию, на 6 групп в зависимости от веса (в кг) (см. Рис. 4). На основе экспертных оценок водителя было установлено, что  для первой группы клиентов из-за незначительного веса заказа

Рис. 4. Количество покупателей у водителя № 1 в зависимости от веса продукции (в кг)

Предположим, что исследуемые временные резервы увеличивается линейно при переходе к каждой следующей группе покупателей на

параметр  В работе для демонстративных целей предполагается, что

3 минуты. Таким образом, первый водитель будет приезжать ко всем покупателям не за 10 минут до конца временного окна, а за время, увеличивающееся в зависимости от размера заказа (к примеру, к контрагенту с весом продукции в 30 кг, водитель вынужден приехать за 15 минут до конца временного интервала). Отметим, что рассматриваемая эвристика может быть внедрена в любую торговую компанию путем изменения параметра  Обобщенные сведения серии детерминированных задач к стохастической представлены ниже:

Таблица 4

Маршруты водителя № 1 в случае разных параметров надежности **

** — Вариант 1 — решение детерминированной задачи для водителя № 1 (без временных резервов) (см. Приложение 6), вариант 2 — решение детерминированной задачи для водителя № 1 с разными временными резервами в зависимости от размеров заказов (см. Приложение 7), вариант 3 — решение детерминированной задачи для водителя № 1 с максимальным временным резервом в 12 минут у каждого покупателя (см. таблицу 3). Сравнение резервов производится относительно варианта № 1.

Заметим, что оптимальный маршрут второго водителя представленный в Приложении 10, остался без изменений, поскольку априори удовлетворял ограничениям компании. Таким образом, в работе удалось приблизить решение серии детерминированных задач для водителя № 1 к стохастической путем установления различных временных резервов. Отметим, что наименьшее значение целевой функция достигается при параметре надежности равном 12 минутам у всех покупателей. Подчеркнем, что на первой итерации решения, задача рассматривалась для временного резерва в 10 минут. Затем временные окна покупателей сокращались на одну минуту, что привело к максимально возможному резерву времени — 12 минут. Вариант № 2 решения задачи (с временными резервами, зависящими от веса заказов) незначительно отличается от оптимального решения в терминах значений целевой функции. Заметим, что с точки зрения практического применения задачи, данный вариант не будет являться наилучшим, поскольку исследуемая фирма желает приезжать ко всем покупателям минимум за 10 минут до конца отгрузки (в описываемой модификации предполагается наличие меньших временных резервов). Тем не менее, вариант решения детерминированной задачи с различными временными резервами может применяться в других торговых компаниях, задавая приоритетность посещения покупателей в зависимости от параметра  (см. таблицу 4, рис. 4).

После рассмотрения различных модификаций локального решения серии детерминированных задач для водителя №1, в работе предлагается разновидность ЗМТ с учетом дорожной ситуации в конкретный день рабочей недели для того, что максимально приблизить решение детерминированной ЗМТ для двух водителей к стохастической. Воспользуемся современными геоинформационными системами и сценарным анализом, предложив оптимизационную модель сетевого планирования для торговых предприятий. При помощи статистики Яндекс карт, которые фиксируют дорожную ситуацию в каждый момент времени в конкретный день, формируя хранилища информации о пробках на дорогах в разных городах России, удалось оптимизировать маршруты водителей в трех измерениях. С математической точки зрения, были взяты оптимальные планы доставки продукции покупателям у каждого водителя (см. таблицу 3, Приложение 10) и скорректированы на время проезда утром, днем и вечером. Фактически данный этап означал составление двух дополнительных асимметричных матриц расстояний на 2 704 ячейки каждая (первая матрица — утренняя, когда дорожная ситуация «благоприятная», вторая — вечерняя, в случае «неблагоприятных» пробок). «Умеренная» дорожная ситуация уже была получена и описана выше. В итоге, необходимо было решить 2 задачи маршрутизации транспорта с временными окнами покупателей и ограничением на грузоподъемность транспортных средств и получить 4 новых таблицы с маршрутами для двух водителей (см. Приложения 8, 9, 11, 12). Данный процесс автоматизирован в логистической программе «Умные маршруты» и позволяет решить вышеперечисленные задачи за небольшой промежуток времени. Таким образом, рассмотрев три сценария дорожной ситуации в г. Пермь в исследуемый день и скорректировав маршруты на коэффициенты пробок, посчитанные в Яндекс Картах, удалось наглядно представить результаты в виде сетевых графиков (см. Приложение 15, 16). Опишем вкратце две получившиеся сетевые модели, точнее их критические пути, поскольку алгоритм поиска с запретами уже сформировал оптимальные решения, которые получались бы при расчете критических путей в данных объемных сетевых диаграммах (с всеми возможными дорогами между каждым покупателем). На практике, водителям будут выдаваться именно разработанные в исследовании сетевые графики, представленные в Приложениях 15, 16, поскольку они являются не только удобными представлениями 6 решений оптимизационных задач коммивояжера, но и индивидуальными планами каждого водителя с собственными полными и свободными резервами времени. Отметим, что резервы времени рассчитываются каждым водителем индивидуально и самостоятельно по мере фактического посещения каждого контрагента в дневной заявке. Так, полный резерв времени, будет равен позднему времени отбытия от конкретного покупателя за вычетом раннего времени прибытия к данному контрагенту в случае «неблагоприятной» и «благоприятной» дорожной ситуации соответственно. Свободный резерв времени, подразумевает также вычет реального времени разгрузки у каждого покупателя. Таким образом, водитель будет знать и самостоятельно ориентироваться во времени между заказами в зависимости от разной дорожной обстановки утром (когда практически нет дорожных пробок), днем (когда они умеренные) и вечером (когда дорожная ситуация оценивается в 5 — 8 баллов). Характеризуя критические пути сетевых моделей, отметим, что у каждого покупателя имеется номер посещения водителем по порядку, среднее время проезда к нему в случае «умеренной» дорожной обстановки, раннее время прибытия при «благоприятных» пробках и позднее время отбытия при «неблагоприятной» дорожной ситуации. Более того, в сетевых моделях, представленных в Приложениях 15, 16, имеется время разгрузки у отдельно взятого покупателя (поскольку модели строились на примере конкретного дня). На практике, водитель не может знать и даже приблизительно оценивать сколько времени ему потребуется у каждого контрагента, поэтому полные и свободные резервы времени водитель сможет посчитать самостоятельно только после реальной доставки.

Отметим, что помимо сценарного анализа пробок, в работе приводится пример позиционного управления логистикой компании, который в отличие от программного управления, учитывает возникновение непредвиденных дорожных ситуаций. Алгоритм поиска с запретами обладает возможностью улучшать решение на каждой итерации, что на практике можно использовать для перерасчета маршрутов водителей, в случае если они выбиваются из графика. Рассмотрим данную особенность алгоритма на следующем примере. Представим, что водитель № 1, выезжая от контрагента № 16, в 15:07, столкнулся с дорожными работами, либо у него сел аккумулятор, подзарядка которого, как и проезд по объездной дороге, заняли бы в сумме 30 минут. Таким образом, водитель не успеет к покупателю № 17 в г. Добрянка и предприятие понесет потери продукции (возникает отрицательный свободный резерв времени):

Таблица 5

Пересчет маршрута в реальном времени (на примере водителя 1)

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

 

Несмотря на то, что водитель не выехал к контрагенту № 17, алгоритм поиска с запретами пересчитал последующий маршрут, изменив последовательность покупателей с № 21 по № 23 (для сравнения последовательности см. Приложение 6). Более того, фирма, возможно, успеет перепродать неотгруженную продукцию, так как водитель приедет в депо за 1:51 до конца рабочего дня. В итоге, произойдет оптимизация логистики в реальном времени (без данного решения водитель пытался бы доехать до

г. Добрянка и только на подъезде понял бы, что опоздал, к примеру, из-за дорожного затора на мосту). Случай с возникновением излишнего свободного резерва времени (простоя у покупателя) представлен в Приложении 14. Используя механизм пересчета маршрутов водителей, компания ООО «Фабрика еды», как и другая торговая фирма сможет предотвратить риск недоставки продукции контрагентам при помощи позиционного управления логистикой.

Итак, в результате решения модифицированной задачи маршрутизации транспорта с временными окнами покупателей и ограничением на грузоподъемность автомобилей, удалось найти оптимальные планы (маршруты) доставки продукции контрагентам максимально приближенные к реальности при помощи авторских сетевых моделей планирования и управления логистическими процессами. Повторно отметим, что данная оптимизация проводилась при наличии фактических маршрутов водителей ООО «Фабрика еды», а именно было известно реальное время доставки и разгрузки, а также вес заказов. Анализируя литературу, можно прийти к выводу, что для оптимизации логистики торгового предприятия, помимо построения оптимальных маршрутов, необходимо исследовать сам процесс формирования заказов у фирмы-поставщика. Данную процедуру применяют для фильтрации покупателей и упрощения формирования оптимальных маршрутов водителей.

В результатах исследования, будет частично затронут процесс отбора заявок покупателей компанией ООО «Фабрика еды». Рассмотрим следующую ситуацию. В один из дней у фирмы имеется 40 заказов покупателей. Зная, что водитель максимум сможет развести только 30 (достаточно с плотным временным графиком доставки), необходимо каким-то образом сократить количество заявок до 30, максимизировав дневную выручку компании и загрузку автомобиля (по весу или объему). Очевидно, что в случае когда количество заявок у одного водителя не превышает 30, необходимость сокращения числа заказов априори отпадает — водитель с большой вероятностью сможет доставить всю заказанную продукцию (необходимо решить ЗМТ и предоставить водителю сетевые графики). Проблема отбора заявок возникает тогда, когда их больше 30 (например, возникает излишний спрос на мясопродукты ближе к выходным дням) и предприятию необходимо фильтровать покупателей в зависимости от суммы и размеров заказов. Именно после получения результатов ЗМТ, в данном исследовании удалось осознать необходимость решения задачи о ранце в дополнение к задаче коммивояжера и вкратце рассмотреть стандартную задачу о рюкзаке (12 — 14) в сфере логистики. Вернемся к предыдущему примеру, когда у ООО «Фабрика еды» за день до доставки имеется 40 заказов на одного водителя и представим один из вариантов их оптимального сокращения. Введем два основных предположения при решении стандартной задачи о ранце. Во-первых, заполняя автомобиль заказами, будем ориентироваться на цену и вес каждого заказа. Отметим, что в исследовании не случайно был взят именно вес, как один из параметров модели, а не объем заказа, поскольку габариты груза напрямую зависят от его веса (чем больше мясной продукции в заказе, тем больше объем). Во-вторых, в работе, решение задачи о рюкзаке не является основной целью, поэтому анализируется классическая формулировка задачи о ранце 0-1 (любой заказ может быть загружен в автомобиль только один раз) без дополнительных ограничений исходя из исключительно демонстративных целей. Отметим, что данную задачу можно так же подробно анализировать, как и задачу маршрутизации транспорта, применяя различные алгоритмы и рассматривая различные модификации задачи.

Для наглядности решения будем ориентироваться на суммарную грузоподъемность автомобиля равную 300 — 305 кг, поскольку эмпирически выяснено, что 300 кг продукции соответствует максимальному объему рефрижератора у транспортного средства для доставки охлажденного мяса. Более того, на практике было выяснено, что один водитель максимум может увести 30 заявок в день, не опаздывая ни к одному покупателю и не перерабатывая. Таким образом, экономико-математические ограничения являются не случайными, и диктуются реальными условиями работы фирмы ООО «Фабрика еды». Поскольку анализируемая задача является дополнением к общему процессу поиска оптимального плана перевозок торгового предприятия, решим ее достаточно простой эвристикой — «жадным алгоритмом», доказав оптимальность полученного решения. На первом этапе оптимизации, одному из двух водителей компании, был предложен список из 40 заказов, а также их веса и цены. Водителю необходимо было в течении 5 минут, отобрать к себе в автомобиль 30 заявок, которые бы максимизировали выручку фирмы и не превышали 300 — 305 кг. По истечении времени, водитель ООО «Фабрика еды», предоставил решение с выручкой равной

172 770 рублей и суммарным весом 326 кг, что было больше грузоподъемности автомобиля. В ходе решения задачи водитель использовал простейший алгоритм сортировки всех заказов по убыванию цены: отбирались первые 30 заявок в списке (с максимальной ценой), без учета веса (см. Приложение 19).

Следующим этапом решения стала, так называемая, «случайная эвристика» — заявки были расположены по мере их фактического поступления и были взяты первые 30 из списка. В результате выручка стала равна 149 350 кг, а вес стал меньше и составил 284 кг (см. Приложение 17). Заметим, что первые два вида решений задачи о рюкзаке являются неоптимальными поскольку не учитывают одновременно два входных параметра — вес и цену заказа, а учитывают их по отдельности. Рассмотрим решение задачи жадным алгоритмом, в основе которого лежит простая эвристика — необходимо поделить цену на вес заказа (получив стоимость 1 кг заявки в рублях) и отсортировать их по убыванию данного показателя, взяв первые 30. В итоге, выручка получилась 162 480 рублей, а вес 302 кг, что удовлетворяет всем ограничениям задачи (см. Приложение 18). Очевидно, что существуют различные алгоритмы решения задачи о рюкзаке, которые находят более точное и оптимальное решение, чем вышеописанное. Жадный алгоритм был использован для простоты и наглядности достижения лучшего решения задачи. Полученный список заявок является оптимальным, поскольку, он выгодней, чем список сформированный водителем на практике. Решение сотрудника фирмы, имея большую выручку, должно быть сокращено по весу, что приведёт к снижению дохода в среднем на 12 000 — 13 000, что является меньшим значением по сравнению с жадным алгоритмом.

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

Возвращаясь, к решению задачи маршрутизации транспорта с временными окнами и ограничением на грузоподъемность автомобилей, можно заметить, что не все контрагенты в исследуемом дне являются «проблемными» для фирмы. Оказалось, что существует закономерность между длинной временного интервала доставки у i-го покупателя и возвратами продукции (Azi, et al., 2012). Как отмечалось выше, в исследовании рассматривается случайно выбранный день в котором один из водителей опоздал в г. Добрянка к покупателю с временным окном отгрузки в 90 минут. Более того, работая с базой данных предприятия, а именно с реальными дневными отгрузками на протяжении 6 — 7 месяцев (с того момента, когда в фирме появился второй водитель), было выяснено, что возврат продукции осуществляется покупателями у которых временной интервал доставки менее двух часов. Отметим, что неудовлетворенность некоторых контрагентов, связанная с опозданием водителей, не коррелировала с местонахождением покупателей, поскольку у

ООО «Фабрика еды» были «проблемные» заявки как в центре города Пермь, так и в удаленных районах. Таким образом, чтобы избежать дополнительных транспортных расходов по вине водителя, предприятию необходимо не только решить ЗМТ и задачу о рюкзаке, но и сформировать оптимальный дизайн контракта для новых контрагентов и, по возможности, перевести некоторых имеющихся клиентов на модифицированный договор поставки.

В результатах данной работы предлагается способ оптимизации дизайна контракта торговой фирмы, путем создания сигналов для покупателей к расширению временных окон доставки. Таким образом, компания еще до начала формирования оптимальных маршрутов сможет оптимизировать транспортировку продукции. Предприятие может создать такие стимулы тремя способами: сделать скидку за расширение временного интервала доставки, ввести надбавку к цене за сокращение окна получения продукции или мотивировать покупателя не денежным способом. Рассмотрим каждую стратегию, выбрав наиболее состоятельную и эффективную. Представим, что ООО «Фабрика еды» делает минимальную скидку в размере 5 % от общей суммы заказа при установке временного интервала доставки продукции не менее 2 часов. Покупатели (рациональные экономические агенты), узнав о скидке изменят свои временные окна так, чтобы любым способом достичь даже незначительного уменьшения стоимости заказа. С другой стороны, фирме — поставщику, необходимо сравнить свои издержки от возвратов продукции, например, за месяц, с альтернативными издержками (упущенной выручкой при введении скидки). Для ООО «Фабрика еды», упущенные выгоды от использования скидки превышают издержки возвратов продукции покупателей в 4 — 5 раз, что является недопустимым для торгового предприятия. Анализируя стратегию повышения цены (суммы заказа) для покупателей, можно также прийти к ее неэффективности. В условиях современной экономической ситуации, а также наличия большого количество конкурентов на рынке мясной продукции, ООО «Фабрика еды» не сможет постоянно увеличивать цены выше рыночных. Придерживаясь стратегии «штрафования» покупателей за сужение временного окна доставки, фирма-поставщик будет либо «отпугивать» новых контрагентов, либо нести убытки из-за их отсутствия. Таким образом, два варианта оптимизации контракта, апеллирующие к цене продукции предприятия или к обще сумме заказа являются неэффективными, поскольку не приводят к изменению поведения контрагентов. Покупатели либо переключатся на другого поставщика из-за высокой конкуренции на рынке, либо будут продолжать устанавливать те интервалы доставки, которые им необходимы.

Рассмотрим возможность оптимизации контракта торговой фирмы не денежным способом. Приведём простой пример избегания дополнительных юридических издержек фирмой, заключающей неоклассический контракт с покупателем. Известно, что в контракте невозможно прописать все будущие разногласия сторон, поэтому в договор добавляется фраза: «Все конфликтные ситуации регулируются на основании законодательства Российской Федерации и решаются через суд», которая позволяет прояснить как способы действий покупателя, так и сократить юридические издержки фирмы. С точки зрения логистики, торговому предприятию также необходимо подобрать такие слова в контракте, которые бы косвенно повлияли на решение контрагента о сокращении временного окна доставки. Примером такой фразы для ООО «Фабрика еды» могут являться следующие предложения: «Покупатели, ожидающие доставку продукции в интервал длительностью менее двух часов, имеют вероятность 5 — 7 % опоздания водителя. Контрагенты, устанавливающие время отгрузки более двух часов, имеют близкую к максимальной вероятность доставки продукции вовремя». Таким образом, фирма будет сигнализировать покупателям, о возможности не доставки продукции в ожидаемый интервал, что создаст стимул для контрагентов расширить свои временные окна отгрузки. В результате данного анализа выяснилось, что для управления логистикой на предприятии недостаточно одного решения задачи маршрутизации транспорта и задачи об оптимальной загрузке транспортных средств. Фирме необходимо также проводить анализ своих покупателей, выявляя «проблемных» и оптимизировать контракты с ними.

Заключение

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диплом

В исследовании рассматривалось решение задачи маршрутизации транспорта с ограничением на грузоподъемность и временными окнами покупателей на примере компании ООО «Фабрика еды», находящейся в

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

В работе проводится анализ основных научных статей по теме. Выяснено, что первое упоминание о КЗК (TSP) появилось в трудах У. Гамильтона и Т. Киркмана в середине ХIХ века. Математики не обладали необходимыми алгоритмами решения задачи и могли сформулировать лишь предполагаемый результат — маршрут торговца, проходящий через каждый город только один раз с возвращением в исходный (при этом необходимо объехать все города без исключения). Данный маршрут (цикл) получил название гамильтонов. По мере развития математики, а именно такого направления как комбинаторная оптимизация, ученые стали придавать задаче коммивояжера не только практический, но и теоретический смысл. Литература на тему КЗК росла экспоненциально и, к примеру, на 2016 год существует более 6 000 исследований, связанных с задачей коммивояжера.

Классическая задача коммивояжера решается в логистической практике редко, поскольку ее ограничения не учитывают, например, пробки на дорогах, маршрутизацию нескольких машин в автопарке и временные окна доставки клиентам. Приближая КЗК к реальности, была разработана модификация задачи с ограничением на количество автомобилей — задача маршрутизации транспорта — ЗМТ (VRP). Важно заметить, что задача маршрутизации транспорта с ограничением только на количество машин редко используется на практике, поскольку оптимизация реальных бизнес процессов подразумевает введение и анализ различных специфичных ограничений.

Учитывая вышеперечисленные факты, в выпускной квалификационной работе рассматривалось решение КЗК и ЗМТ на примере фирмы, находящейся в г. Пермь. В определенный момент времени у предприятия возникли сложности с доставкой продукции покупателям: из-за увеличивающего количества заявок, водители ездили неоптимально, не попадая в желаемое время отгрузки, за что покупатели штрафовали предприятие, делая возврат в полном объеме. Несмотря на то, что большинство контрагентов были лояльны к времени отгрузки и опозданиям, ООО «Фабрика еды» каждый месяц несла потери продукции на 25 — 35 тыс. рублей по вине то одного, то второго водителя. В среднем один из водителей не успевал к одному контрагенту в неделю и фирма теряла продукцию на

4 — 5 тыс. рублей.

Решить данную проблему для компании удалось, рассмотрев ЗМТ, которая несмотря на относительную простоту постановки задачи, решалась достаточно сложно. В работе рассматривались как детерминированная, так и стохастическая постановки задачи. На основе ряда эвристик в исследовании показывается, что компания может решать серию детерминированных задач, приближая получаемые маршруты к реальности. Безусловно, для того чтобы решить описываемые задачи, необходимо было подобрать репрезентативные данные и выбрать эффективные методы нахождения наилучших маршрутов. Информационной базой исследования стали данные бухгалтерского и управленческого учета фирмы ООО «Фабрика еды». На первом этапе их сбора были изучены все дневные отгрузки продукции покупателям начиная с лета 2015 года (когда у фирмы появился второй водитель) по настоящее время. Для данного исследования был взят один из вышеуказанных дней с 51 покупателем. Заявки именно данных клиентов выполнялись 18 сентября 2015 года. Необходимо отметить, что в анализируемом дне деятельности предприятия первый водитель не успевал к контрагенту № 19, который находился в г. Добрянка (70 км от центра г. Пермь). Более того, предприятие понесло убытки по вине одного из водителей в размере 14 560 рублей (штраф в виде возврата продукции в полном размере). Предприятие пыталось поставить данный заказ в приоритет, что привело к неудовлетворению последующих заявок у нескольких покупателей в центре города. Таким образом, появилась необходимость при помощи математических методов сформировать такую последовательность обслуживания покупателей, чтобы сократить данные потери, как в исследуемой повторяющейся заявке, так и в ряде других.

В итоге, первым этапом анализа стало решение задачи коммивояжера, в имитационном случае: если бы у фирмы имелся один водитель и он бы развозил 51 заявку индивидуально. Данный эксперимент был проведен, поскольку необходимо математически проверить нужность второго водителя в фирме. Затем, фактические маршруты двух водителей были оптимизированы и решение было сравнено с реальным. На следующем этапе работы задача CVRPTW была разделена на 2 части и каждая решена по отдельности. В конечном счете, целая задача маршрутизации была решена с учетом трех сценариев дорожной ситуации в г. Пермь: «благоприятного», «умеренного», «неблагоприятного» и полученные результаты были обобщены в виде сетевой модели планирования и управления (PERT). На всех вышеперечисленных этапах работы использовались два основных метода решения задачи: муравьиный алгоритм и метод поиска с запретами.

Таким образом, в результате исследования удалось математически подтвердить управленческое решение руководства ООО «Фабрика еды» о найме второго водителя, поскольку одно транспортное средство фактически не могло развезти 51 заявку с 11:00 до 20:00. Затем были по отдельности оптимизированы реальные маршруты водителей, полученные при помощи видеорегистраторов транспортных средств и личных интервью. Оказалось, что муравьиный алгоритм смог найти оптимальные решения для двух задач коммивояжера, улучшив фактические дневные планы водителей. В ходе анализа также выяснилось, что торговому предприятию нет необходимости проводить любого вида кластеризацию данных, поскольку не удастся ни в одном из случаев учесть все реальные бизнес ограничения.

В работе рассматривались две модификации ЗМТ: разновидность с пересчетом маршрута, а также с учетом реальной дорожной ситуации. Воспользовавшись пересчетом маршрутов, водители смогут учитывать непредвиденные ситуации, что максимально приблизит решение детерминированной задачи к стохастической. Более того, при помощи современных геоинформационных систем и сценарного анализа была предложена оптимизационная модель сетевого планирования для торговых предприятий. Заметим, что на практике для оптимизации логистики торгового предприятия, помимо решения различных модификаций ЗМТ, необходимо исследовать сам процесс формирования заказов, который с экономической точки зрения моделируется в виде задачи о ранце 0-1

(0-1 knapsack problem). Решив данную задачу, удалось отфильтровать заявки покупателей, что фактически привело к дополнительной оптимизации логистики.

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

Таким образом, применяя описываемые в работе способы оптимизации транспортной логистики, ООО «Фабрика еды», как и другая торговая фирма, сможет сократить свои транспортные расходы на 25 — 35 тыс. рублей в месяц.

Исследование можно развить в следующих направлениях:

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диплом

) Решить более сложные модификации задачи о ранце 0-1, внедрив в нее различные ограничения;

) Рассмотреть подробно оптимизацию дизайна контрактов, описав функцию полезности покупателей;

) Оптимизировать процесс управления скоропортящимися запасами.

Таким образом, несмотря на ряд как достоинств, так и недостатков математических моделей КЗК и ЗМТ, применение их в торговых компаниях для минимизации транспортных затрат может быть достаточно эффективно.

Список литературы

1) Бронштейн, Е.М., Заико, Т.А. (2010), «Детерминированные оптимизационные задачи транспортной логистики», Автоматика и телемеханика, №10, с. 133-147.

) Иглин, С.П. (2003), «Решение некоторых задач теории графов в MATLAB», Математика в приложениях, №4, с. 28-33.

) Костюк, Ю.Л. (2013), «Эффективная реализация алгоритма решения задачи коммивояжера методом ветвей и границ», Прикладная дискретная математика, Т. 20,№ 2, с. 78-90.

) Курейчик, В.М., Мартынов, А.В. (2014), «Об алгоритмах решения задачи коммивояжера с временными ограничениями», Информатика, вычислительная техника и инженерное образование, Т. 16, № 1, с. 1-13.

) Лебедев, Б.K., Лебедев, О.Б. (2012), «Моделирование адаптивного поведения муравьиной колонии при поиске решений, интерпретируемых деревьями», Известия ЮФУ. Технические науки, Т. 132, № 7, с. 27-34.

) Лопатин, А.С. (2005), «Метод отжига», Стохастическая оптимизация в информатике, Т.1. №1, с. 133-139.

) Макконнелл, Дж. (2004), Основы современных алгоритмов, М: Техносфера,с. 233-255.

) Пантелеев, А.В., Метлицкая, Д.В., Алешина, Е.А. (2013), Методы глобальной оптимизации. Метаэвристические стратегии и алгоритмы, Москва, с. 105-130.

) Хайруллин, Р.З. (2014), «Математическое моделирование развоза грузов по разветвленной сети автодорог», Вестник МГСУ, № 7, с. 184-191.

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

) Чеблоков, И.Б., Ченцов, А.Г. (2012), «Об одной задаче маршрутизации с внутренними работами», Вестник Удмуртского университета. Математика. Механика. Компьютерные науки, № 1, с. 96-119.

) Applegate D., Cook J. (1993), “Solving large-scale matching problems”, Network Flows and Matching, pp. 557-576.

) Azi, N., Gendreau, M. and Potvin, J.-Y. (2012), “A dynamic vehicle routing problem with multiple delivery routes”, Annals of Operations Research, pp. 1-10.

) Baldacci, R., Mingozzi, A. and Roberti, R. (2012), “Recent exact algorithms for solving the vehicle routing problem under capacity and time window constraints”, European Journal of Operational Research, Vol. 218, No. 1, pp. 11-26.

) Clarke, G., Wright, J.W. (1964), “Scheduling of vehicles from a central depot to a number of delivery points”, Operations Research, No. 12, 568-581.

) Croes, G.A. (1958), “A method for solving traveling salesman problems”, Operations Research, No. 6, pp. 791-812.

) Dantzig, G. B., Ramser, J. H. (1959), “The truck dispatching problem”, Management Science, Vol. 6, No. 1, pp. 80-91.

) Desrochers, M., Desrosiers, J. and Solomon, M. (1992), “A new optimization algorithm for the vehicle routing problem with time windows”, Operations Research, No. 40,pp. 342-354.

) Dorigo, M. (1999), “Ant algorithms for discrete optimization”, Artificial life, Vol. 5, No. 2, pp. 137-172.

) Irnich, S. (2008), “Solution of real-world postman problems”, European Journal of Operational Research, Vol. 190, No. 1, pp. 52-67.

) Kan, A.R., Lenstra, J.K. (1975), “Some simple applications of the travelling salesman problem”, Operations Research, No. 26, pp. 717-734.

) Laporte, G. (1992), “The Vehicle Routing Problem: An Overview of Exact and Approximative Algorithms”, European Journal of Operational Research, Vol. 59, No. 3,pp. 345-358.

) Lin, S., Kernighan, B.W. (1973), “An effective Heuristic Algorithm for the TSP”. Operations Research, No. 21, pp. 498-516.

) Lorini, S., Potvin, J.-Y. and Zufferey, N. (2011), “Online vehicle routing and scheduling with dynamic travel times”, Computers & Operations Research, Vol. 38, No. 7,pp. 1086-1090.

) Moon, I., Lee, J. and Seong, J. (2012), “Vehicle routing problem with time windows considering overtime and outsourcing vehicles”, Expert Systems with Applications, Vol. 39, No.18, pp. 13202-13213.

) Rodriguez, A., Ruiz, R. (2012), “A study on the effect of the asymmetry on real capacitated vehicle routing problems”, Computers & Operations Research, Vol. 39, No. 9,. 2142-2151.

) Schneider, J.J., Bukur, T. and Krause, A. (2010), “Traveling salesman problem with clustering”, Journal of Statistical Physics, Vol. 141, No. 5, pp. 767-784.

) Solomon, M. M. (1987), “Algorithms for the Vehicle Routing and Scheduling Problems with Time Window Constraints”, Operations Research, Vol. 35, No. 2, pp. 254-265.

) Toth, P. and Vigo D. (2002), The Vehicle routing problem, Society for Industrial and Applied Mathematics, p. 386

) Wiener, J. M., Ehbauer, N. N. and Mallot, H. A. (2009), “Planning paths to multiple targets”, Psychological Research, No. 73, pp. 644-658.

Приложение 1

 

Приложение 2

Оптимизированный маршрут водителя № 1

 

Приложение 3

Оптимизированный маршрут водителя № 2

 

Приложение 4

Маршрут водителя № 1. Кластерный анализ

 

Приложение 5

Маршрут водителя № 2. Кластерный анализ

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диплом

 

Приложение 6

Решение ЗМТ («умеренная» дорожная ситуация) для водителя № 1

без учета резерва времени в 10 минут

 

Приложение 7

Решение ЗМТ («умеренная» дорожная ситуация) для водителя № 1

с различными резервами времени (зависимость от веса)

 

Приложение 8

Решение ЗМТ («благоприятная» дорожная ситуация) для водителя № 1

 

Приложение 9

Решение ЗМТ («неблагоприятная» дорожная ситуация) для водителя № 1

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Подробнее

 

Приложение 10

Решение ЗМТ («умеренная» дорожная ситуация) для водителя № 2

 

Приложение 11

Решение ЗМТ («благоприятная» дорожная ситуация) для водителя № 2

 

Приложение 12

Решение ЗМТ («неблагоприятная» дорожная ситуация) для водителя № 2

 

Приложение 13

Адреса доставки продукции

 

Приложение 14

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Цена диплома

Таблица 1

Маршрут с незапланированным простоем

Таблица 2

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

 

Приложение 15

Задача о рюкзаке (фактическая дневная заявка)

 

Приложение 16

Решение задачи о рюкзаке «жадным алгоритмом»

 

Приложение 17

Решение задачи о рюкзаке. Эксперимент с водителем. Цена по убыванию

 

Нужна помощь в написании диплома?

Мы - биржа профессиональных авторов (преподавателей и доцентов вузов). Сдача работы по главам. Уникальность более 70%. Правки вносим бесплатно.

Заказать диплом

Средняя оценка 0 / 5. Количество оценок: 0

Поставьте оценку первым.

Сожалеем, что вы поставили низкую оценку!

Позвольте нам стать лучше!

Расскажите, как нам стать лучше?

1331

Закажите такую же работу

Не отобразилась форма расчета стоимости? Переходи по ссылке

Не отобразилась форма расчета стоимости? Переходи по ссылке