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

книги / Типовые задачи оперативного управления непрерывным производством

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
19.07 Mб
Скачать

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

Вторая постановка задач рассматривается в [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* — век­

тор постоянных в правых частях для

/-го интервала.

Соседние файлы в папке книги