Комплексное оценивание при выборе оптимального плана производства на тактическом уровне с учетом нечетких критериев и ограничений Текст научной статьи по специальности «Математика»

УДК 658.5.012+122.011.56 ББК 22.193

КОМПЛЕКСНОЕ ОЦЕНИВАНИЕ ПРИ ВЫБОРЕ ОПТИМАЛЬНОГО ПЛАНА ПРОИЗВОДСТВА НА ТАКТИЧЕСКОМ УРОВНЕ С УЧЕТОМ НЕЧЕТКИХ КРИТЕРИЕВ И ОГРАНИЧЕНИЙ 1

Вожаков А. В.2, Гитман М. Б.3, Федосеев С. А.4

(Пермский государственный технический университет, Пермь)

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

Ключевые слова: календарное планирование производства, нечеткие множества, математическое моделирование, оптимизация.

1. Введение

На крупных предприятиях с дискретным производством и сложной высокотехнологичной продукцией возникает серьезная проблема при календарном планировании - эффективно распределить множество технологических операций по имеющимся производственным мощностям (тактический уровень планирования [3]). На предприятиях, как правило, используются не оптимальные, а лишь допустимые производственные планы, по-

1 Работа выполнена при поддержке РФФИ (проект № 10-08-00539).

2 Артем Викторович Вожаков, кандидат технических наук ([email protected]).

3 Михаил Борисович Гитман, доктор физико-математических наук, профессор ([email protected]).

4 Сергей Анатольевич Федосеев, кандидат физико-математических наук ([email protected]).

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

С точки зрения оптимизации все решения, получаемые в результате выработки планов, можно разделить на:

- оптимальные по всем критериям;

- оптимальные по одному/нескольким критериям;

- приемлемые по всем ограничениям.

2. Постановка задачи

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

Рассмотрим многокритериальную задачу оптимизации планирования производства на тактическом уровне [4]. Для описания календарного плана производства вводится матрица Р. Ее элементы рІй определяют количество операций 1-го вида, запланированных на день с номером й, І є 1, ..., Ж, й є 1, ..., Т, где Т -количество дней в плановом периоде; Ж - количество видов операций. Введем матрицу баланса номенклатурных единиц В. Ее элементы Ьій определяют количество номенклатурных единиц с номерами і є 1, ..., С (С -количество номенклатурных позиций в справочнике), находящихся на цеховых складах на конец дня с номером й (й є 1, ..., Т); значение индекса й = 0 используется в матрице баланса для определения остатков номенклатурных единиц на начало планируемого периода. Очевидно, что значения матрицы напрямую зависят от плана производства Р и от начального значения остатков на складе г. Мат-

рица В зависит от главного календарного плана производства (ГКПП) [2]. В дни, когда запланирован выпуск продукции, происходит отгрузка готовой продукции со склада, уменьшая тем самым остаток. Элементы матрицы В можно определить следующим образом:

где гі - остаток /-ой номенклатурной единицы на начало периода; g - главный календарный план производства (ГКПП); а/}- -состав изготавливаемых изделий; w — количество операций в /ом технологическом этапе; N - количество изготавливаемых

Очевидно, что при составлении плана производства должны учитываться следующие ограничения:

1. Суммарное количество операций одного вида в плане-графике должно быть равно общему количеству операций данного вида, которые необходимо выполнить в плановом периоде согласно ГКПП:

где О - количество операций, которые необходимо выполнить.

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

ная загрузка оборудования, данная величина может быть нечеткой, т. е. может быть задана в виде нечеткого числа [4]; qk - ко-

Г, і є1, С, й = 0,

N

3 =1

номенклатурных позиций в справочнике; =

(3) Ри0ш £ АЧк,

личество рабочих центров к-го типа; іі - трудоемкость операции і-го вида; Єї - номер типа рабочего центра для совершения операции і-го вида; Е - количество типов рабочих центров, имеющихся на производстве.

3. В любой момент времени остаток номенклатурных единиц не может быть отрицательным. Математически данное ограничение может быть записано как

(4) Ьій > 0 для всех і є 1, С , й є 0, Т .

Введем четыре частных критерия оптимальности плана - графика производства:

1. Критерий комфортности производства - план - график производства должен быть скомпонован таким образом, чтобы операции одного вида запускались в производство как можно большими партиями (в этом случае не потребуется переналадка оборудования при переходе от выполнения одной операции к другой):

(5) Ji = ZZßd ® min,

l=1 d=1

ß IX Pid * o, где ßid = i 0 0 l0, Pu = 0

2. Производство должно быть равномерным (уровень загрузки рабочих центров должен изменяться постепенно, без резких перепадов):

T-1 E

(6) J 2 = ZZ \zkd - zkd+i\ ® min,

d=1 k=1

Z Pld^lßlk

lk rl, e, = k,

- j=i___________________ /? - j

kd

где -------------------, Л = Г г *

tkqk I0, е/ * к

3. Риск срыва плана производства должен быть минимален. Максимальное опережение производства т является нормативом предприятия и измеряется в днях. Будем говорить, что продукция выпущена с опережением, если продукция произведена

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

¿1 * <г/2

1 Ш р І

(7) л = 22 -к,2-•' '

где X, = і

1, 2 8 м * 0,

і=1

-> Ш1П.

1=1

N

0, 2 §• = 0.

1=1

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

(8) /4 = тща(ё) ® шт.

,є1,Т

где со(,) =

Ш Т

• 22 Рій = °’

і=1 и=,

Ш Т

0’ 22 Ріи *0-

і=1 и=,

На основе предложенных частных критериев введен обобщенный критерий оптимальности с использованием расширенного специального нечеткого множества над частными критериями оптимальности

^ = {т\/'^\; т2/'^2 ; ; т4/'^4 } ,

где т £ [0; 1] , 7 £ 1,4, - экспертная оценка значимости 7-го критерия.

Для сравнения двух элементов J и Jr2 множества J будем использовать четкую функцию от нечеткого аргумента (специальный индекс ранжирования [1]):

(9) H(Jr1,Jr2) = sign Ct,

7/r1 Tr1 Ur1 Tr 2 где C> m J - m J ;

u

i доставляет max|(mr1 • Jp - MpJi2 )/u\;

- функция принадлежности (значимость), Jr ; mi2 - функция принадлежности (значимость) Jrt 2;

U = max(Jir1, J[2) , i e 1,4 .

При этом если значение i единственно и если sign C, = «+», то Jr1 > Jr2; если sign C = «-», то Jr1 < Jr2.

Если значение i не единственно, то определим k - количество равных по модулю экстремумов (k < n) и вычислим

1 = £ sign C. i=1

Если 1 = 0, , то Jr1 = J2; если 1 < 0, то Jr1 < J2; если 1 > 0, то Jr1 > J2.

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

Найти значение матрицы плана производства pld, l e 1, W, d e 1, T , приводящее к минимуму выражение

(10) Jr ® min, где Jr = {m1/J1; m2/J 2; m3/J3; m4/J4},

где J1, J2, J3, J4 - частные критерии оптимальности плана производства при ограничениях (1)-(4); ль /л2, ц3, л4 - оценки значимости частных критериев.

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

- показатели должны быть возрастающими, т. е. при улучшении значения критерия показатель возрастает;

- показатели должны быть отнормированы к единице, при этом за 0 берется самое худшее значение критерия, а за 1, соответственно, самое лучшее;

- значения показателей должны в точности соответствовать значениям соответствующего частного критерия оптимальности. Для любых двух значений частного критерия оптимальности (А и В) соответствующие им показатели оптимальности (РА и РВ) должны находиться в отношении: если А лучше В, то РА > РВ, если А = В, то РА = РВ;

- минимальные и максимальные значения показателей должны быть достижимыми в рамках конкретной задачи.

Будем использовать следующие показатели оптимальности плана производства, которые соответствуют всем требованиям.

Показатель комфортности:

Ш Т

1

1=1 ¿=1

1—

Т

Показатель равномерности:

E

W

£ Pid^iPik

гДе zkd =

і=i

hVk

I1, et = k, 0, e, * k.

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

T-1

£ rd

it;

(13) J2 =

d=1

2 T -1

Показатель риска срыва производства:

d*<г/2

T d-d* >0 ^ W P t id-d* l

■y XlV p

Zj Zj j* Zj

(14) J3* = 1 -

=1 d* t! q/e,E

d <t/2 T d-d*>0

£ £

гДе Xd = i

d=1 d =1

N

1 £ gtd *0

d *

i=1

N

0 £ gid = °.

i=1

Показатель минимальности срока производства: (15) J4 = 1 - max ®(d) -1

где w(d) =

de1,T T - 2

WT

d, ££ Pi« = 0

l=1 u=d W T

0 ££Pi« *0.

l=1 u=d

3. Решение задачи, оценка результатов

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

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

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

2. Использование трех критериев вместо четырех позволит на порядок снизить машинное время, требуемое для проведения исследования. В противном случае проведение столь длительных исследований стало бы невозможным.

ГКПП рассматриваемого примера представлен в таблице 1.

Таблица 1. Пример производства крупных партий

Дата выпуска (день от начала планового периода) Число изделий

22 50

44 50

66 50

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

1. Максимально комфортный план производства. В данном случае использовался только критерий комфортности производства, и задача, по сути, стала однокритериальной. Экспертные оценки важности критериев следующее: ^1 = 1, ^2 = 0, ^3 = 0. Результат обозначим как «комфортный план производства».

2. Достаточно комфортное и максимально равномерное про-

изводство. Экспертные оценки важности критериев следующее: [л1 = 0,7, [л2 = 1, = 0,3. Результат обозначим как «экспертный

план производства».

Полученные результаты будут использоваться для сравнения.

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

ны при решении данной задачи. Для исследования были выбраны значения функции принадлежности от 0,1 до 1 с шагом в 0,1

В результате проведенных вычислений было получено множество решений (размер множества совпадает с количеством проведенных вычислений). Из данного множества были исключены решения, заведомо не оптимальные. Заведомо неоптимальными решениями были признаны значения Iх, для которых существует решение I2 такое, что для всех I выполняются условия: J¡11 < 12*; Я < т2. Таким образом, получили набор решений (множество выбираемых решений), представленный в таблице 2.

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

- линейная свертка значений частных критериев без учета важности критериев;

- линейная свертка значений частных критериев с учетом важности критериев в качестве множителей;

- поиск наилучшего значения с использованием индекса ранжирования (9).

Дополнительно стоит указать физический смысл индекса ранжирования (9): наилучшим значением множества I будет признано то значение 1г, у которого больше количество пиковых значений частных критериев, с учетом их значимости. В таблице 2 представлены оценки по всем трем способам выбора наилучшего. Для удобства таблица отсортирована по значению линейной свертки частных критериев без учета их важности. Значения в столбце «Ранг» указывают оценку решения с помощью индекса ранжирования (9), где 1 - наилучшее значение, 2 - наилучшее из оставшихся решений и т.д.

Из таблицы видно, что оценки, полученные с помощью индекса ранжирования, сильно отличаются от оценок, полученных с помощью линейных сверток. Далее каждое решение было рассмотрено группой экспертов, которые проводили оценки качества решений с точки зрения сбалансированности и отражения поставленных целей. Эксперты указали, что значения, полученные с помощью индекса ранжирования (9), являются более адекватными. Таким образом, данный индекс следует рассматривать как альтернативу линейной свертке, однако для каждой задачи вопрос применения той или иной оценки является открытым. Итак, наилучшим признано решение №10 из таблицы 2 -это решение соответствует оценкам важности частных критериев для данной задачи: = 1, ^2 = 0,6, = 0,8.

Таблица 2. Множество выбираемых решений задачи

№ т т т 3 * 3 * 2 3 * X3* /=1 X м* 1=1 Ранг

1 0,4 1 0,3 0,40 0,97 0,36 1,74 1,24

2 0,7 1 0,3 0,74 0,84 0,37 1,9 1,48

3 1 0,8 0,4 0,82 0,75 0,4 1,97 1,58

4 0,4 0,7 1 0,42 0,59 0,96 1,98 1,55

5 0,5 0,4 1 0,53 0,47 1 2,01 1,45

6 0,5 1 0,7 0,49 0,89 0,62 2,01 1,58

7 0,5 1 0,8 0,39 0,85 0,79 2,04 1,68 3

8 1 0,3 1 0,87 0,27 0,97 2,11 1,92 5

9 1 0,4 0,5 0,93 0,43 0,76 2,13 1,49

10 1 0,6 0,8 0,83 0,5 0,8 2,14 1,77 1

11 0,9 0,6 0,8 0,81 0,51 0,83 2,16 1,70 2

12 1 0,4 0,9 0,89 0,31 0,95 2,16 1,87 4

13 0,9 0,3 1 0,83 0,36 0,98 2,18 1,84

Отметим, что решение с экспертной оценкой функции принадлежности (предыдущий пример) не попало в число наилуч-

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

Рассмотрим решение задачи при наилучших значениях функций принадлежности: ¡л1 = 1, ¡л2 = 0,6, ¡лъ = 0,8.

На рис. 1 представлен график изменения показателя равномерности плана производства ра для четырех планов производства: допустимый план до оптимизации (допустимый); план после оптимизации по критерию комфортности (комфортный); план после оптимизации по комплексному критерию с экспертными оценками (экспертный) и план после оптимизации по комплексному критерию с наилучшими оценками важности критериев (далее оптимальный). Точки на графике указывают значение показателя ра между днями d и d + 1.

Раїмомермоєть загрузки оборудования

-До оптимизации - Комфортность -Экспертный -Оптимальный

Рис. 1. Равномерность загрузки оборудования

На рис. 2 представлены значения показателя минимальности риска срыва производства для допустимого, комфортного, экспертного и оптимального производственных планов.

Зініаче%:іиіе гпопсазателпяі рмспса

1,2

1

т

| 0,8

§

а

і о.е

Ш 0,4

о

с

0,2

0

□Допустимый

©Комфортный

□ Экспертный

ШОптимальный

Рис. 2. Риск срыва производства

На рис. 3 представлены значения показателя комфортности У* для допустимого, комфортного и комплексного производственных планов.

1 1

I °>9

0 « «

1 ! #0'7

0 0,6 £ 0,5

1 °’4 I

Е 0,2 0,1 о

Іиачеию© иоссазагшаша тм< прошшадапша

□Допустимый Ы Комфортный □ Экспертный 0 Оптимальный

Рис. 3. Комфортность производства

4. Заключение

Из приведенных выше графиков видно, что после оптимизации допустимого плана производства по комплексному критерию произошли следующие характерные изменения:

1. Значение показателя комфортности улучшено по сравнению с допустимым планом и меньше оптимального значения (комфортный план) менее чем на 10%.

2. Показатель равномерности загрузки оборудования значительно улучшен по отношению к допустимому и комфортному плану производства, однако на 40% меньше, чем в экспертном плане производства.

3. Значение показателя риска срыва производства в два с половиной раза лучше, чем в допустимом плане, и менее чем на 20% меньше оптимального значения;

4. Производственный график стал более равномерным за счет ограничения на досрочный выпуск деталей.

5. В последние дни планового периода нагрузка снижается, чтобы дать возможность оперативно решить возможные проблемы.

Литература

1. БОРИСОВ АН., АЛЕКСЕЕВ А.В., МЕРКУРЬЕВА Г.В. Обработка нечеткой информации в системах принятия решений - М.: Радио и связь, 1989. - 304 с.

2. ГАВРИЛОВ Д.А. Управление производством на базе стандарта MRP II. - СПб.: Питер, 2002.

3. СТОЛБОВ В.Ю., ФЕДОСЕЕВ С.А. Модель интеллектуальной системы управления предприятием // Проблемы управления.- 2006.- №5. - С. 36-39.

4. ФЕДОСЕЕВ С.А., ВОЖАКОВ А.В., ГИТМАН М.Б. Модель календарного планирования производства с нечеткими целями и ограничениями // Системы управления и информационные технологии. - 2009. - №3. - С. 21-24.

COMPLEX EVALUATION FOR OPTIMAL TACTICAL PRODUCTION SCHEDULE SELECTION UNDER FUZZY CRITERIA AND CONSTRAINTS

Artem Vozhakov, Perm State Technical University, Cand. Sc., ([email protected]).

Mikhail Gitman, Perm State Technical University, Doctor of Science, professor ([email protected]).

Sergey Fedoseev, Perm State Technical University, Cand. Sc., assistant professor ([email protected]).

Abstract: The model is considered for evaluation ofperformance of an automated production scheduling and planning system under fuzzy criteria and constraints. Optimal schedules obtained using different models of integrated estimation are calculated and compared.

Keywords: production scheduling, fuzzy sets, mathematical modelling, optimization.

Статья представлена к публикации членом редакционной коллегии Д. А. Новиковым