УДК 630*377
ОПТИМИЗАЦИЯ ТЕРРИТОРИАЛЬНОГО РАСПРЕДЕЛЕНИЯ ТРАНСПОРТНЫХ
ПОТОКОВ В ЛЕСНОМ КОМПЛЕКСЕ
ассистент кафедры промышленного транспорта, строительства и геодезии
А. С. Сушков1
заведующий кафедрой технологии и машин лесозаготовок, доктор технических наук,
профессор О. Н. Бурмистрова2 1 - ФГБОУ ВПО «Воронежская государственная лесотехническая академия» 2 - ФГБОУ ВПО «Ухтинский государственный технический университет»
[email protected], о1аа.Ьигт@таП.ги
Рассматриваемую проблему можно ного программирования, представленных в
решить с помощью класса задач нелиней- виде
™П К (У)Х1 + С 2 (У)Х2 + - + Ср (У)Хр + С0 (У) } ^ А (У) Х > Ь(у)
А2(у)> Ъг(У)
Ар (у)Хр > Ър (у)
Dy > 0
(1)
где, коэффициенты целевой функции а , технологические матрицы А{ и векторы Ъi являются функциями вектора связывающих переменных у. Если зафиксировать значения этих переменных, то задача (1) становится задачей линейного программирования с переменными х. Эта задача имеет блочно-диагональную структуру технологической матрицы и поэтому распадается на р независимых подзадач. Задачи подобного типа являются нелинейным обобщением блочных линейных задач со связывающими переменными. Большинство методов решения таких задач основано на идеях декомпозиции и сводится к решению последовательности линейных подзадач для различных фиксированных значений вектора у.
Задачи типа (1) часто возникают при
составлении графиков производства и распределения продукции лесоперерабатывающих комплексов, когда оптимальный план составлен для целого комплекса предприятий.
При планировании производства возникают следующие задачи:
1. Наилучшим образом распределить различные виды лесопродукции между ЛПК, а внутри - между технологическими установками.
2. Рационализировать материальные потоки между различными установками. (Решение задачи усложняется наличием промежуточных видов продукции, которые могут быть направлены на другие ЛПК для дальнейшей переработки.)
3. Определить оптимальные варианты, режимы использования технологических
установок с учетом различных факторов и условий. Зависимость свойств потоков от этих переменных, как правило, носит ярко выраженный нелинейный характер.
4. Обосновать наилучшие схемы транспортировки лесопродукции при опре-
деленных ограничениях на объемы потребления.
5. Найти рациональный вариант закрепления потребителей за лесопромышленным комплексом ЛПК.
Рассмотрим теперь возможность ис-
Процесс Смешивание Транспортировка Потребитель
Рис. 1. Система производства и распределения лесопродукции
пользования нелинейной модели (1) для решения последней задачи. Пусть предприятия Rl и R2 производят продукты Р1 и Р2 и отправляют их в пункты потребления Т1 и Т2. Суммарный спрос (за весь период) на продукцию в этих пунктах предполагается известным и пусть D(Pi , Т) - величина спроса на продукт Рi в пункте Т. На предприятии Rl три вида лесоматериалов (сырья) М1, М2, М3 могут быть использованы для получения продуктов Р1 и Р2. Кроме того, лесо-
(
материал М1 может быть использован на предприятии R2 для смешения с лесоматериалами М4 и М5. Эту ситуацию схематически можно представить в виде схемы, представленной на рис. 1. Исходя из общей постановки задачи, модель оптимального размещения выпуска однородной продукции с учётом минимизации совокупных затрат на производство и транспорт может быть записана в следующем виде
тт З1 =
т
I I, Шг )
+
г=1 1=1
\
^ тт
(2)
при ограничениях
п
Шг = I Шу < М\ 1=1
доп
1,т;
т
QJ = 1 QlJ > Р7, у = 1,п; 7=1
QlJ > 0,
(4)
(5)
где 1,7 - индексы возможных пунктов производства и соответственно пунктов потребления продукции;
Qi, Мдоп - объёмы производства и соответственно предельно допустимые их величины в 1-ом пункте;
Qlj - интенсивность транспортно -
экономической связи между пунктами I и пунктами у;
fl )- функция производственных затрат в 1-ом пункте;
Су - удельные транспортные затраты на перевозки между пунктами I и у;
Р7
- необходимый объём потребления и соответственно объём доставляемой в у -й пункт продукции.
Из анализа модели (2)-(5) можно установить, что определяющее влияние на размещение производства оказывает не абсолютный уровень транспортных затрат, а дифференциация их по различным связям: чем меньше эта дифференциация, тем слабее влияние транспорта. При малых или относительно постоянных затратах, оптимальный план определяется из условия минимума только производственных затрат.
На практике задача размещения и развития производства решается при следующих условиях.
1. Пункты и объёмы потребления рассматриваемой продукции фиксированы,
критерий - минимум приведённых затрат на производство и транспортировку лесо-продукции.
2. Оптимизация размещения производств выполняется по всему предполагаемому объёму выпуска лесопродукции (а не только по его приросту).
3. В экономико-математических моделях транспортные блоки имеют линейную структуру с постоянными коэффициентами целевой функции - значениями удельных показателей транспортных затрат, не зависящими от объёма перевозок.
4. Расчёты выполняются изолировано по отдельным вариантам без непрерывного обмена информацией, поэтому суммарные загрузки всех видов транспортных потоков лесопродукции ещё неизвестны.
Показатели транспортных затрат должны отражать лишь дополнительные составляющие объёмов и направлений перевозок лесопродукции.
При решении транспортных оптимизационных задач предполагается, что найденные транспортно - экономические связи на этом этапе и транспортные потоки будут осуществляться в реальных текущих грузопотоках. Возникает проблема модификации исходной информации и процедуры решения с целью снижения отрицательного влияния хозяйственного механизма при реализации решения. Здесь можно выделить три следующих этапа:
а) подготовка информации для решения задачи о размещении;
б) выбор решения по минимуму приведенных затрат;
в) реализация принятого решения (определение рациональных потоков сырья и готовой продукции по минимуму тарифных затрат).
Если тарифы совпадают с соответствующими дифференциальными удельными приведёнными транспортными затратами, то потери качества решения не происходит.
Таким образом, проблема учёта транспортного фактора в рассматриваемых задачах требует не только оптимизационного подхода на основе критериев экономической эффективности, но и учёта особенностей, которые вносит в реализацию оптимизационных планов специфика действующего производственного процесса.
Рассмотрим модель размещения пунктов потребления (хранения, переработки) лесопродукции и задачи функционирования лесопромышленных предприятий в разных условиях хозяйствования.
Все они преобразуются в блочно -диагональные линейные задачи, если зафиксировать значения некоторых переменных. Рассмотрим проблему в общем сетевом ракурсе.
Минимизировать
{С 'х + / (у )}^тш (6)
При ограничениях
Ах + F (у )> Ь. (7)
х > 0,у е 5, (8)
где х — п - мерный вектор;
А- постоянная диагональная матрица;
С — п - мерный вектор констант.
Если зафиксировать значение у, то (1)-(3) переходят в задачу линейного про-
граммирования. Величина /(у) - функция (возможно нелинейная) р- мерного вектора у, F(у)- вектор - функция,
5 - произвольное подмножество Ер.
Многие важные в практическом отношении задачи могут быть приведены к виду (1)- (3). Если 5' - множество р - векторов с неотрицательными целочисленными компонентами, а F,f - линейны, то (1)-(3) окажется задачей линейного программирования.
В соответствии с этим алгоритм задачи (1)-(3) решают в следующей последовательности: вначале фиксируется некото*
рое значение у е 5 и решается задача: минимизировать С' х при ограниче-
ниях
Ах > Ь — F
V)
(9)
Затем на основе решения двойственной задачи декомпозиции определяется возможность улучшения решения, полученного на первом шаге и находится новое значение вектора у. Таким образом, полностью используются преимущества частичной линейной задачи, что особенно важно, когда матрица А имеет специальную структуру (например, блочно- диагональную или транспортного типа). В этом случае задача (4) сравнительно легко решается. Эти преимущества не реализуются при использовании алгоритмов, в которых х и у изменяются одновременно.
При формулировке задачи используются как дискретные, так и непрерывные переменные (дискретные - строить или не строить предприятие или склад, непрерыв-
ные соответствуют объёмам производства и поставок лесопродукции).
Функция текущих затрат, представленная на рис. 2, учитывает хорошо известный факт: с увеличением мощности предприятия возрастают условно - постоянные затраты, однако их увеличение компенсируется снижением условно - переменных затрат, что отражает преимущества крупных предприятий (то есть на таких предприятиях затраты на единицу продукции меньше, чем на мелких предприятиях). Сведение задачи размещения предприятий к задаче линейного программирования позволяет, как правило, сделать процедуру
ьО Й
1
СП
3 р , 3 ^ , 3 с - множества номеров
пунктов производства, складов и потребителей соответственно.
Переменные т^ ограничены сверху и снизу и должны удовлетворять дополнительным условиям из-за наличия взаимоисключающих вариантов мощности к.
перебора вариантов более эффективной.
Введём обозначения:
Г у - спрос у-го потребителя;
Су - удельные транспортные затраты на перевозку единицы лесопродукции из пункта I потребителю у;
Ху - объём перевозок из пункта I в
пункт у;
т- количество лесопродукции, привозимой (накапливаемой) предприятием I (складом I) при работе по к-ому варианту (с затратами на прирост единицы продукта, равными Slk );
Чтобы учесть условия, введём переменные У1Л , принимающие значения нуль или единица (0 или 1), так что
Угк = 0 ^ т^ = 0, (10)
У1к =1 ^ Ъг,к-1 ^ тг,к ^ Ъг,к (11)
и не более чем одна переменная Угк = 1 для каждого пункта производства
Объем производства Рис. 2. Кусочно-линейная разрывная функция затрат
Z mik = Z -ij k
ikmik
г. Соответствующие ограничения можно представить в виде
Уik ьгл—1 < mik < у г(12)
I Уik =1 Уik = 0 или 1
k
Так как спрос должен быть удовлетворен полностью, возникают ограничения Iху = ^у ,1 е Jc. (13)
г
Переменные ху и mik связаны между собой ограничением
ху (14)
С ^ см>
2 = ! С«х„ +15.
г, ] г к
Первый член этого выражения - транспортные издержки, а второй и третий соответствует затратам на строительство предприятия и производство продукции. Выражения (12)-(17) могут быть приведены к виду (6)-(8), в которой F и f - линейные функции, а 5 - множество возможных значений переменных ук . Выражения (12)-(17) решаются в два приёма: вначале фиксируются некоторые значения переменных Уiк , и решается линейная задача. Затем в результате решения целочисленной задачи определяются новые значения Уik [1].
Выводы. Оптимальные решения задачи, двойственной к линейной, используется для построения дополнительных ог-
или вида xlj =
Z xij = Z xjk,j e Jcw, (15)
ieJ
keJ,
(то есть с каждого склада входящие и выходящие потоки лесопродукции равны между собой). Ёмкость склада ограничена, поэтому должно соблюдаться ограничение
I ху < Л у, ] е Jw . (16)
iеJp
Общие затраты должны быть минимальными, то есть
+ Zyik (fi,k-1 - Sikbik-1Ц min
(17)
ik
раничении в целочисленном задаче, что уменьшает множество допустимых альтернатив, проверяемых на оптимальность. Таким образом, эта процедура может быть рассмотрена как схема, «автоматизирующая» перебор вариантов, в котороИ информация о ранее рассмотренных случаях используется для выделения случаев, которые следует рассматривать далее.
Библиографический список
1. Бухтояров В.Н., Быков А.В., Суш-ков А.С. Оптимизация транспортно-технологических связеИ предприятии лесопромышленного комплекса // Строительные и дорожные машины. 2011. Вып. 12. С. 20-23.