книги / Типовые задачи оперативного управления непрерывным производством
..pdf2. Определяются не только материальные потоки и запасы, но и график ремонта отдельных агрегатов. При определении графика учитываются плановые ограниче ния на выпуск продукции, переработку сырья и техноло гические связи между отдельными агрегатами. Кроме того, учитываются факторы, связанные с рациональным планированием занятости рабочей силы при производст ве ремонтных работ.
Вторая постановка задач рассматривается в [96]. Эти задачи носят комбинаторный характер с большим чис лом переменных, что существенно затрудняет их практи ческое решение. Кроме того, ряд факторов, учитывае мых при планировании загрузки рабочей силы, трудно формализуется.
Исследование подобных задач требует предваритель ного изучения более простых задач с фиксированным графиком ремонтов. Поэтому в настоящей главе основ ное внимание будет уделяться этим задачам.
Задача календарного планирования может быть сформулирована для всего технологического комплекса, например предприятия в целом или для отдельного про изводства. В реальном производстве с непрерывными технологическими процессами между отдельными агре гатами имеются буферные емкости. В зависимости от их величины и характера использования задача решается для всего предприятия или отдельного производства.
Принимая во внимание типовые операции по преоб разованию материальных потоков, введенные ранее, сле дует отметить, что такие независимые производства до статочно часто представляют собой либо отдельные ти повые операции, либо типовые комплексы этих операций. Таким образом, типовые задачи календарного планиро вания можно связать со структурами комплексов типо вых операций.
В дальнейшем будем рассматривать задачу календар ного планирования применительно к комплексу, состоя щему из параллельно соединенных смесительных опера ций. Перенесение полученных результатов и процедур на другие параллельно-последовательные комплексы типо вых операций не представляет трудностей.
Критериями оптимальности для календарного плана работы комплекса могут быть приняты следующие [11]:
1. Максимум прибыли от производства готовых про дуктов на конец планового периода.
'Z. минимум суммарных издержек от производства заданного количества готовых продуктов в плановом периоде, т. е. минимум себестоимости товарной про дукции.
3. Максимум выпуска определенных готовых продук тов (ранжированных по значимости для производства) при выполнении плановых заданий на остальные (или все) готовые продукты.
Кроме того, могут использоваться и другие критерии оптимальности, например равномерность выпуска про дукции в целом или по отдельным видам при заданной прибыли.
Следует отметить, что введенные критерии рассматри ваются при полной определенности графиков поставки и качественных характеристик поступающих сырьевых компонентов. При значительных неопределенностях в указанных параметрах могут быть введены другие кри терии, например максимизация вероятности выполнения плана [96].
Как и ранее, будут рассмотрены детерминированные постановки. Перейдем теперь к строгой формулировке задач.
а) Задача технико-экономического планирования
Рассмотрим задачу оптимального технико-экономиче ского планирования промышленного комплекса, состоя щего из параллельно соединенных смесительных опера ций (рис. 6-1). На вход комплекса поступают потоки исходных сырьевых компонентов. После переработки на агрегатах комплекса на выходе получается ряд готовых продуктов с гарантированными показателями качества. С целью «развязки» работы комплекса с работой постав щиков и потребителей на входе и выходе комплекса име ется парк емкостей для хранения сырья и готовой про
дукции. |
|
_ |
|
|
поставка |
|
Введем обозначения: |
Xj — агрегированная |
|||||
/-го |
сырьевого |
компонента |
за |
планируемый |
период; |
|
Xjs— количество /-го |
сырьевого |
компонента, |
идущее |
|||
на производство s-ro готового продукта; |
— коли |
|||||
чество производимого 5-го готового продукта; |
п — чис |
|||||
ло |
сырьевых |
компонентов; |
m — число готовых про |
дуктов.
Уравнение модели каждой 5-й параллельной смеси тельной операции можно записать:
Xjs—UjJs, /= 1, . . я, 5=1, т,
где Ujs—подлежащий выбору рецепт смешения,
eV s; Vs — множество допустимых рецептов смешения.
Рис. 6-1. Пример схемы про мышленного комплекса для формулировки задачи тех- нико-экономического плаиироваиия.
а — входные |
емкости; |
б — па |
раллельные |
смесительные опе |
|
рации; о — выходные |
емкости. |
а
Рассмотрим случай, когда все Vs заданы с помощью систем равенств и неравенств вида
—s = ] .........т ’
7=1 _
V i s < v i s < v is, i = \ , . . . , / s;
П
vIS |
|
|
где Vis — i - я качественная |
характеристика |
(показатель |
качества) для 5-го готового |
продукта; v ^ |
- i - я качест |
венная характеристика (показатель качества) /-го ис ходного компонента; u,-s, v u — соответственно верхнее и
нижнее гарантированные значения /-го показателя каче
ства для 5-го готового продукта;
В этом случае модель комплекса может быть записа на следующим образом:
Ys = X,st 5— 1.......т; |
( - ) |
/=1 |
6 1 |
|
|
> s х п < ^ X j; < y la у x rs, 5 = 1 |
/,. (6-2) |
r=l |
/=> |
1 |
Естественно также, что для всех /, s |
(6-3) |
Xie> 0 t Ys^O. |
Потребуем, чтобы неиспользованные исходные ком поненты по объему не превышали отведенного под них объема парка входных емкостей, т. е.
__ |
т |
__ |
о < ^ |
- 2 |
х /5+ / ; х< / ; \ у= 1 ...... (6-4) |
|
S=1 ' |
где / “ —запасы /-го сырьевого компонента к началу пла
нируемого периода; / “*— ограничения объема резервуар
ного парка по /-му сырьевому компоненту. Аналогично произведенные, но не отгруженные потребителю готовые продукты по объему не должны превышать объема пар ка выходных емкостей, т. е.
0 < - ? s + 2 Х,а + Г * < Т " , s = |
1,... ,m, |
(6-5) |
/=i |
|
|
где — плановая поставка 5-го готового |
продукта |
по |
требителям; /®ых—запасы 5-го продукта к началу плано ного периода; /“ых — ограничения объема резервуарного парка по 5-му продукту.
Необходимо найти план работы комплекса на теку щий период, т. е. выпуск готовых продуктов Ys, 5=1, ...
..., т и затраты компонентов на выпуск этих продуктов XjS, /=1, ..., п, 5=1, ..., ш, удовлетворяющие условиям (6-1) —(6-5) и обеспечивающие максимум (минимум) функционала
п т |
|
F = % 2 cisXjs, |
(6-6) |
/= ! s=I
выражающего один .из выбранных ранее критериев опти мальности.
Для случая максимизации прибыли
Cjs==Cs
где cs—отпускная цена s-ro готового продукта; Cj — се бестоимость /-го исходного компонента.
Для случая минимизации себестоимости CjS=Cj.
Для случая максимизаций выпуска s-ro готового про дукта функционал (6-6) запишется в виде
/-1
В дальнейшем без потери общности будем рассматри вать задачу максимизации прибыли, т. е. с функцио налом
ит
/= 1 5 — 1
Первоначально можно искать решение задачи в тер
минах XjS, а на втором этапе с помощью (6-1) |
получить |
|||||||
значения переменных Ys. |
|
|
|
|
||||
Выпишем полученную задачу. |
|
т, обеспе |
||||||
Требуется найти Xjs, j= 1, |
п, s= 1, |
|||||||
чивающее |
максимум |
функционалу |
(6-7) |
при |
наличии |
|||
ограничений |
|
|
|
|
|
|
|
|
3 * 5 2 |
Х Г8 < |
2 |
Vkn X i* ** ViS 2 |
X r S ' |
5 = 1 |
...........ТП' |
|
|
Г=1 |
/=1 |
|
/=I |
|
|
|
|
|
|
|
|
|
ls\ |
|
|
|
|
0<:- у*+2^ / 5+ Сх< / Г - |
|
»l> (6-8) |
||||||
|
|
/ - 1 |
|
|
|
|
|
|
0 < * / |
m |
|
/ = I......«; |
|
||||
~ |
2 |
+ |
|
|||||
|
|
s—I |
|
|
|
|
|
X js ^ 0,
Сформулированная задача (6-8) является задачей линейного программирования с размерностью, определя емой размерностью модели промышленного комплекса. В результате решения находится оптимальный техникоэкономический план рассматриваемого комплекса, т. е
6) Задача календарного планирования
При внедрении полученного оптимального плана не обходимо детализировать этот план применительно к пе риодам времени, меньшим, чем период технико-экономи
ческого планирования. Например, следует разделить го довой план на квартальные или месячные, в свою оче редь квартальные или месячные на декадные или пяти дневные и т. д.
Предположим, что весь рассматриваемый период пла нирования разделен на Т интервалов. В этом случае мо дель комплекса можно записать следующим образом:
y8t = '%xist* |
t — 1...... T; |
(6-9) |
||
/=i |
|
|
|
|
v1st = 2 |
( XJst / 2 *rst |
< * = |
1. ••• . la* |
(6-10) |
/-> |
r=l |
|
|
|
где приняты те же обозначения, что и в |
(6-1), (6-2), но |
все переменные отнесены к соответствующему t-му ин тервалу, t= 1, ..., Т.
Выпускаемая продукция должна удовлетворять за данным показателям качества, которые будем считать неизменными во времени
v |
is |
Ü |
<т>. |
- |
|||
|
1st |
|
Для каждого временного интервала заданы объемы сырьевых компонентов хц, j= 1, ..., п, t=\, ..., Т, посту пающих на переработку в этот интервал. При этом
т_ |
_ |
2 |
/ — !»•••> |
t=\ |
|
Если интервалы выбрать так, что можно было прене бречь динамикой запасов внутри интервала, то условия ограниченности емкостей можно проверять лишь в нача ле и в конце каждого интервала. Для входных емкостей запишем:
/ = 1...... п.
(6- 11)
Ограничения, связанные |
с |
Парком |
Выходных |
ем |
||||
костей, запишутся в виде |
|
|
|
|
|
|
||
° - |
2 л , + |
i l |
+ с |
|
|
s = |
1 •••••т ' |
|
|
/=1 |
/=1 |
|
|
|
|
|
|
|
|
/ = 1 , . . . , Г - 1 , |
|
(6- 12) |
||||
где г/я/ ■— заданный |
режим отгрузки |
5-го |
готового |
про |
||||
дукта. |
|
|
|
|
|
|
|
|
Если заданы лишь предельные значения отгрузки, то, |
||||||||
вводя новую |
переменную |
у ° т( |
|
(отгрузку |
5-го готового |
|||
продукта в l-й период), вместо |
(6-12) |
будем иметь: |
|
|||||
О |
|
|
вых ^ |
J вых |
5 = 1 , . . . , т , |
|
||
1= |
1=1 |
s |
|
s |
|
|
|
|
|
|
|
|
|
|
t= 1,... ,Т — 1;
от
ysi^ y
рт
s i '
где y s l , у — заданные |
пределы отгрузки. |
|
|
Требуется при условии выполнения оптимального те |
|||
кущего плана |
|
|
|
т |
к* 5 = 1 , |
, m |
(6-13) |
2 ^ = |
|||
/—1 |
|
|
|
найти показатели работы комплекса на каждом времен
ном |
интервале, т. е. |
значения |
y st |
и X jSt> s = l , .... m, / = |
||||
= 1 , |
..., |
п, t= 1, .... |
Т, удовлетворяющие |
|
ограничениям |
|||
(6-11), |
(6-12) и обеспечивающие максимум прибыли |
|||||||
|
|
T |
/ m __ |
п |
ni |
\ |
|
|
|
|
* = 2 |
2 « * - - 2 |
e/ 2 * / - |
J |
• |
(6' 14) |
|
|
|
1=1 \ i = I |
/ = I |
s=l |
|
|
Как и в случае задачи технико-экономического пла нирования, можно выразить переменные задачи через X jstt а после решения найти с помощью (6-9) значения переменных y st.
Выпишем полученную |
задачу. Требуется |
найти Xjisit |
||||
/ = 1, |
п, s = l , |
т , |
/= 1 , |
Т, доставляющие мини |
||
мум функционалу |
|
|
|
|
||
|
|
Т |
п |
!П |
_ |
|
|
|
? = 2 |
2 |
2 |
-«#)■*/- |
|
|
|
/ = 1 |
/= 1 5 = 1 |
|
|
|
при наличии ограничений |
|
|
|
|||
|
и |
« |
|
|
/г |
1 |
|
^ » 5 |
jt-*-jut |
is |
2 ^ js t > 7 1 ’ • • • |
> |
|
|
/=1 |
/=1 |
|
|
J=I |
|
|
о<-2Р*+2 2^,+'Г<7Г' |
|
||||
|
/=i |
<=i/= |
1 |
|
||
|
|
/г= 1 |
, , |
T — 1; |
|
|
|
|
i £ |
x |
lsl=Y*s, |
|
|
|
|
I- 1j-i |
|
__ |
|
|
|
t |
^ |
t m |
|
o^S-Sï-SS-^+C ^
/ = ] |
/= 1 5 = 1 |
Xjst^O, /= 1 ,... ,n, |
s = l,...,m , f = l,... ,7\ |
Полученная задача является задачей линейного про граммирования. В отличие от' задачи (6-8) ее можно назвать задачей детализации текущего плана или про сто детализированной задачей. Размерность детализи рованной задачи (6-15) зависит как от размерности мо дели комплекса, так и от числа интервалов времени Т.
Задача (6-15) для реальных непрерывных комплек сов даже при небольшом числе интервалов планирова ния приводит к задачам линейного программирования достаточно большой размерности (включающих сотни переменных и ограничений). Для решения подобных за дач могут быть использованы методы разложения (де композиции) .
6-2. Применение общих процедур декомпозиции для решения детализированной задачи
Предварительно введем переменные г ' и |
харак |
теризующие остатки исходных и готовых продуктов во входных и выходных емкостях на конец /-го временного интервала:
t t т
|
г“ = |
2 * у « - 2 |
2 * / « + 7Г |
- ' = |
1...... п; |
|||||||
|
|
1= |
|
1=] |
5=1 |
|
|
|
|
|
|
|
С = |
- 2 |
*- + 2 |
2 |
|
•*« + С 1- * = 1 |
m |
||||||
|
|
/=1 |
/=1 |
/= |
|
|
|
|
|
|||
И Л И |
|
|
|
//2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
г";ч = г -/х_1+ |
л> —2 - v * |
|
|
|
|
/= !» • |
||||||
|
|
|
|
5=1 |
|
|
|
|
|
|
|
|
г5 /'Х = С - . — |
“1“2 |
"'7s*’ |
|
^= |
Ь • • • > |
^ = |
1,...> т, |
|||||
причем |
|
|
|
/=| |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В Ч ____ |
IBX . |
^ В Ы Х ____ ГВЫХ |
|
|
|||||
|
|
|
6 . |
“ i , |
I ^ |
“1* |
« |
|
|
|||
|
|
|
JO |
|
I |
|
|
SO |
S |
|
|
|
Запишем детализированную задачу (6-15) в следую |
||||||||||||
щей форме: |
|
|
|
|
|
|
|
|
|
|
|
|
Y <Ра< |
ViJxiS(<0, |
i = \ , . . . , l s, 5 = |
1,... ,т, t— 1,...,7\" |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
П
пх
V
2 (0|, — «III) XI4<0;
/=--1 т
1 |
1 —X |
+ 2 x ist |
x i 1' |
n \ |
|
|
|
i = i , . . . , |
i= l
уВЫХ _1_ >ШЛХ |
и |
_ |
|
5 = 1 , . . . , m; |
|
~5/ 1 s<- , + 2 |
|
|
|
у—1 |
/ = = 1..... я; |
гпх <7nx, |
||
it |
/ |
|
В1ЛХ< |
Г \ |
5 = 1 , . . . |
z st |
5 |
|
2 |
^ x Jst= Y *a, 5==1, — m; |
|
|
||||
|
1/=I |
|
|
|
|
|
|
|
|
|
В Ы Х |
0, |
|
|
|
|
T |
n |
2sf/ |
|
|
|
|
|
m |
|
шах. |
|
|
||
F = |
2 |
2 |
2 (C.-C,)-*/.! |
|
|
||
|
/=1 /=1 s=l |
|
|
|
|
||
На рис. 6-2 представлена таблица со структурой де |
|||||||
тализированной |
задачи. Вектор х* |
является |
вектором |
||||
|
|
|
|
1 |
|
В .Ч В Ы Х |
— |
критериальных потоков для г.-го интервала; |
zf , z/ |
векторы запасов входных и выходных продуктов для то
го же интервала; |
Е;г — единичная матрица; |
B/t — матри |
|
ца |
с элементами |
0 и 1; V/. — матрицы с |
элементами |
Vi't |
V%s ИUis V%jt- |
|
Структура задачи представлена набором независи мых блоков Vt и общей системой ограничений, связы вающих эти блоки (эта система обведена пунктиром). Для решения задач с подобной структурой может при меняться стандартная процедура метода декомпозиции Данцига—Вульфа [20, 97]. Существенно, что общая система имеет сильно разреженную матрицу коэффи циентов 0 и 1. Эта особенность может быть использова на для облегчения реализации процедуры декомпо зиции.
Детализированная задача может быть представлена и как задача со ступенчатой структурой. Для этого вве дем переменные ysi\
П
Уst '— Уsi - 1 I Xjsi » ^ 1 » • • • » ^ S 1, ... , Ш,
/=1
Уso ”
Согласно ограничениям детализированной задачи
ysT==Y*s, s = l, .... m.
С учетом этих переменных ступенчатая структура де тализированной задачи изображена на рис. 6-3,а.
Для ступенчатой структуры, |
изображенной на |
|
рис. 6-3,6, |
D/, В/ — матрицы, а х* — вектор переменных; |
|
сt — вектор |
коэффициентов целевой |
функции; d* — век |
тор постоянных в правых частях для |
/-го интервала. |