Например, Бобцов

Модели расписаний переработки двух видов сырья с различной кратностью использования каскада экстракторов

УДК 641.524.6:664.5
Модели расписаний переработки двух видов сырья с различной кратностью использования каскада экстракторов
Перов А.Г., Косачев В.С., Кошевой Е.П.
Кубанский государственный технологический университет
Использование одного каскада на различных видах сырья увеличивает время ожидания измельченного сырья перед его загрузкой в экстрактор. Снижение этого времени, возможно, при использовании специализированных каскадов и уменьшении кратности их использования. Из этого следует, что при организации краткосрочного оперативного планирования число переходов установки с одного вида сырья на другое должно быть минимально. Определяющим фактором такого расписания становится операция дробления. В этом случае при оптимальной переработке сырья эта стадия для рассматриваемого производства определяет временной график простоев расписания.
Ключевые слова: экстракция, теория расписаний, эфиромасличном сырье.
В реальных производственных условиях экстракторы периодического действия высокого давления (CO2-экстракторы) объединяют в установки (два параллельно работающих экстрактора), имеющие общую систему обеспечения (конденсатор, испаритель, и т.д.). В этом случае процесс переработки сырья проходит в этих экстракторах параллельно (загрузка установка, экстрагирование, слив установка).
Рассмотрим модели бинарных расписаний образующихся при переработке двух видов сырья с различной кратностью использования каскада экстракторов на данном виде сырья. При составлении графика Ганта для моделируемого расписания учитывалось, что при переходе каскада с одного сырья на другое необходима очистка каскада от остатков предыдущего сырья. Основными моделируемыми параметрами процесса являются число каскадов, используемых для переработки данного вида сырья (k) и кратность использования этих каскадов (n1, n2). Функциями отклика в этом случае являются время завершения работ (D) и простой каскадов (∆).
Рассмотрим простейшее расписание, возникающее при использовании одного каскада на двух видах сырья (k=1, n1=1, n2=1). Даже в этом простейшем случае возможно два вида расписаний, имеющих различные функции отклика (Рисунок 1).

Сравнивая эти два графика Ганта, видим, что время завершения работ и простой каскада в первом случае меньше чем во втором. Это соответствует упорядочиванию по правилу Джонсона, т.к. представленные расписания фактически образуют двухстадийную систему и минимум соответствует более короткой первой стадии. Для проверки этого предположения рассмотрим расписания, возникающие при использовании каскада дважды на кориандре и однократно на укропе (k=1, n1=2, n2=1). Как и в предыдущем случае, оптимальное расписание (Рисунок 2) подчиняется правилу Джонсона (D=640, ∆1=110, ∆2=335). Рассмотрим процесс построения модели расписания при многократном использовании каскада на двух видах сырья.
Рис. 1. Диаграмма Ганта экстракции кориандра и укропа однократно на одном каскаде (k=1, n1=1, n2=1)

Рис. 2. Диаграмма Ганта экстракции кориандра и укропа на одном каскаде (k=1, n1=2, n2=1)

В простейшем случае этот вариант расписания соответствует

двукратному использованию экстрактора на обоих видах сырья (k=1, n1=2, n2=2). Очевидно, в этом случае возможно создание четырех расписаний с различными вариантами переработки (Рисунок 3).

Как видно из представленных графиков Ганта оптимальным также

является расписание, удовлетворяющее алгоритму Джонсона.

Следовательно, бинарные расписания образованные на одном каскаде

фактически представляют собой расписания двух машин на нескольких

партиях сырья.

Особенностью данного процесса является тот факт, что при переходе с

одного сырья на другое длительность второй стадии увеличивается на время

очистки каскада. Следовательно, при расчете графика выполнения работ

матричным методом этот факт должен учитываться в матричных

уравнениях.

В общем виде матрица длительностей выполнения работ в этом случае

имеет следующую блочную структуру:

Wb =

d1,1 s2,1

Λ Λ

d1,i Λ s2,i Λ

d1,n1 o2,n1

d1,n1 +1 s2,n1 +1

Λ Λ

d1, j Λ s2, j Λ

d1,n1 + n2 o2,n1 + n2

(1)

где n1 - кратность использования каскада на первом сырье; n2 - кратность использования каскада на втором сырье; d1, i - дробление первого сырья, i = 1, 2, …, n1; d1, j - дробление второго сырья, j = n1+1, n1+2, …, n1+n2; s2, i - Σ (загрузка, экстракция, слив) первого сырья, i = 1, 2, …, n1-1;

s2, j - Σ (загрузка, экстракция, слив) второго сырья, j = n1+1, n1+2, n1+n2-1; o2,n1 - Σ (загрузка, экстракция, слив, очистка) первого сырья; o2,n2 - Σ (загрузка, экстракция, слив, очистка) второго сырья.
Рис. 3. Диаграмма Ганта экстракции кориандра и укропа на одном каскаде (k=1, n1=2, n2=2).

В этом случае матрица завершения работ (Z) может быть рассчитана по

следующей матричной формуле:

Zb =

w1,1

w1,1 + w2,1

∑2 w1,i
{ }i=1
w2,2 + z1,2 ; z2,1 max

Λ Λ

∑j w1,i
{ }i=1
w2, j + z1, j ; z2, j −1 max

Λ Λ

∑n1 +n2 w1,i

{ }i=1

m + z ; z2,n1 +n2

1,n1 + n2 2,n1 + n2 −1 max

(2)

Представленное матричное уравнение позволяет определить не только

завершение всех стадий процесса, но и их начало - вычитанием из матрицы

(Zb) матрицы (Wb). Использование элементов полученной матрицы начала работ (Rb=Zb-Mb) позволяет определить длительность простоев каскадов в ожидании загрузки их сырьем по формуле:

∑ ( )n1+n2

Sb∆ =

r2, j − z1, j .

j =1

(3)

Аналогично возможно построение расписаний для произвольного числа

партий различного сырья и различной кратности использования одного

каскада на этих партиях сырья.

Одним из возможных вариантов реализации матричного алгоритма

является среда инженерных расчетов MathCAD.В качестве исходных данных

могут быть использованы: вектор кратностей использования каскада (myN);

вектор длительностей дробления сырья (myD); вектор длительностей

суммарного времени загрузки экстракции и слива каскада (myS) и вектор

длительностей очистки каскада (myO). В этом случае матрица

длительностей будет содержать столько блоков, сколько элементов

содержится в векторе (myN). Длина каждого блока равна величине каждого

элемента этого вектора. Матрица длительностей (myW) формируется на

основе следующего модуля:

myW := k ← 0

for i ∈ 1 .. rows ( myN )

for j ∈ 1 .. myN i− 1

mw 0 , k ← myD i− 1 mw 1 , k ← myS i− 1

(4)

k← k+ 1

mw 1 , k − 1 ← mw 1 , k − 1 + myO i− 1

mw

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

myZ := k ← 0

for i ∈ 1.. rows(myN)

for j ∈ 1.. myNi−1

mw0, k ← myDi−1

mw1, k ← mySi−1

k←k+ 1

mw1, k−1 ← mw1, k−1 + myOi−1 for j ∈ 0.. k − 1

(5)

j

∑mz0, j ←

mw0, i

i=0

mz1, j ← mw0, 0 + mw1, 0 if j 0
( )mz1, j ← max mz1, j−1, mz0, j + mw1, j otherwise

mz

В качестве примера построения такого расписания можно рассмотреть переработку трех видов сырья (кориандр, укроп и хмель) на основе следующих исходных векторов:

3

myN

:=

 

2

 

5

Вектор кратностей для трех видов сырья

 10 

myD

:=

 

15

 

 25 

Вектор длительностей дробления трех видов сырья

 120 

myS

:=

 

150

 

 180 

Вектор длительностей суммарного времени загрузки экстракции и слива трех видов сырья

(6)

 120 

myO

:=

 

120

 

 120 

Вектор длительностей очистки каскада для трех видов сырья

Использование исходных данных (6) в модуле (4) позволяет получить матрицу длительностей, которая оптимально упорядочена (минимальные длительности вначале первой строки этой матрицы):

myW

=



10 120

10 120

10 240

15 150

15 270

25 180

25 180

25 180

25 180

25 300



.

(7)

Использование этой матрицы в модуле (5) позволяет получить матрицу

завершения работ:

myZ

=



10 130

20 250

30 490

45 640

60 910

85 1090

110 1270

135 1450

160 1630

185 1930



.

(8)

Вычитая из матрицы завершения работ матрицу длительностей,

получаем матрицу начала работ (myR=myZ-myW):

myR

=



0 10

10 130

20 250

30 490

45 640

60 910

85 1090

110 1270

135 1450

160 1630



.

(9)

Матрица (myR): позволяет рассчитать вектор времен простоя каскада

(∆) при выполнении данной производственной программы:

∆ j := myR1, j − myZ0, j ,
где j - текущий элемент вектора ∆, j=1, 2, …, Σ(myN). В данном примере вектор (∆) имеет следующие элементы:

(10)

∆T = ( 0 110 220 445 580 825 980 1135 1290 1445)

(11)

Выводы:
Использование одного каскада на различных видах сырья увеличивает время ожидания измельченного сырья перед его загрузкой в экстрактор. Снижение этого времени возможно при использовании специализированных каскадов и уменьшении кратности их использования. Из этого следует, что при организации краткосрочного оперативного планирования число переходов установки с одного вида сырья на другое должно быть минимально. Определяющим фактором такого расписания становится операция дробления. В этом случае при оптимальной переработке сырья эта стадия для рассматриваемого производства определяет временной график простоев расписания.
Список литературы
1. Танаев B.C., Сотсков Ю.Н., Струсевич В.А. Теория расписаний. Многостадийные системы. М.: Наука, 1989. 322 с.
2. Кафаров, В.В., Мешалкин В.П. - Анализ и синтез химикотехнологических систем. М. : Химия, 1991.-431 с.
3. Колчин А.Ф., Овсянников М.В., Стрекалов А.Ф., Сумароков С.В. Управление жизненным циклом продукции. - М. : Анахарсис, 2002. - 304 с.
4. Конвей Р.В., Максвелл В.Л., Миллер Л.В. Теория расписаний. М., Наука, 1975 г.