CC BY
150
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015. №7 (204). Выпуск 34/1
137
УДК 658.51.012
ОБЗОР МОДЕЛЕЙ УПРАВЛЯЕМЫХ ПРОИЗВОДСТВЕННЫХ ПРОЦЕССОВ ПОТОЧНЫХ ЛИНИЙ ПРОИЗВОДСТВЕННЫХ СИСТЕМ
REVIEW OF GOVERNANCE MODELS PRODUCTION LINES MANUFACTURING SYSTEMS
О.М. Пигнастый O.M. Pignasty
Национальный Технический Университет "ХПИ”, Украина, 61002, Харьков, ул. Фрунзе, 21 National Technical University “KhPI”, 21 Frunze St, Kharkov, 61002, Ukraine
e-mail: pom7@bk.ru
Аннотация. В статье приведен обзор основных моделей управляемого производственного процесса, используемых для проектирования систем управления современными поточными линиям. Выполнен сравнительный анализ разных типов моделей и показаны области применения моделей. Дана краткая характеристика основных параметров моделей. Детально рассмотрены кинетические модели и потоковые модели с использованием уравнений в частных производных.
Resume. The article provides an overview of the main models managed manufacturing process used for the design of control systems with modern flow lines. A comparative analysis of different types of models and shows the application of models. A brief description of the main parameters of the models. A detailed analysis of kinetic models and flow models using partial differential equations.
Ключевые слова: PDE-модель, производственная линия, массовое производство, незавершенное производство, система управления, балансовые уравнения, производственной линии, уравнение состояния, дискретнособытийная, модель, теория массового обслуживания, модель жидкости, Clearing- функция, квазистатический процесс, переходный процесс , стохастический процесс.
Keywords: PDE-model, production line, mass production, work in progress, management system, balance equations of the production line, equation of state, discrete-event model, queuing theory, model fluid, Clearing-function, quasi-static process, transient process, stochastic process.
Введение
Вопросы планирования и управления производством промышленной продукции определяют ключевые задачи для многих теоретических дисциплин. Большое количество работ, посвященных проектированию систем управления производственными процессами, используют математический аппарат теории исследования операций (Arrow K.J.,Karlin S., Бир С., Первозванский А.А. Шкур-ба В.В., Бессонов В.А., Бигель Дж., Прыткин Б.В.), теории оптимизации (Интрилигатор М., Kempf K., Якимович С.Б., системной динамики (Форрестер Дж.), теории массового обслуживания (Gross D, Harris С.) , теории управления запасами (Рыжиков Ю.И.), теории планирования и управления производством (Modigliani F., Holn C., Бусленко Н.П., Дудорин В.И., Емельянов С.В., Митрофанов С.П., Соколицын С.А., Разумов И.М., Балашевич В.А.), статистической теории динамических систем (Власов А.А., Казаков И.Е., Красовский А.А.), статистической теории производственных систем (Азарен-ков Н.А., Демуцкий В.П., Пигнастый О.М., Локтев И.И., Петров Б.Н., Тихомиров И.А., Ходусов В.Д., Armbruster D., Ringhofer C.). Проблемам общей методологии моделирования сложных производственно-технических систем посвящены труды отечественных и зарубежных ученых Arrow K.J., Solow R.M., Вильсона А., Власова В.А., Глушкова В.М., Голдрата Э., Занга В.Г., Канторовича Л.В., Ле-тенко В.А., Разумова И.М., Соколицына С.А., Хакена Г., Шананина А.А., Шкурбы В.В.. Вопросам применения DES-,TQ- и Fluid-моделей при проектирования высокоэффективных систем управления поточными линиями уделено внимание в трудах Cassandras C., Heymann M., Hopp W.J., Eekelen J., Ramadge P., Roset B., Wardi Y., Wonham W., Бусленко Н.П., Брагинского К.А., Лысенко Ю.Г., Коро-бецкого Ю.П., Рамазанова С.К., Руменцева Н.В. Использованы математические модели для принятия управленческих и проектно-технологических решений в условиях неопределенности посвящены работы в работах Buzacott J.A., Harrison J., Shanthikumar J.G., Бессонова В.А., Дудорина В.И., Ермольева Ю.М., Зарубы В.Я., Зарубина B.C., Иванилова И.П., Леонтьева В.В., Лотова А.В., Новикова ДА., Прангишвили И.В., Редькина А.К., Тихомирова А.К., Тихонова А.Н.
Модель стратегического развития предприятия определяется спросом на выпускаемую продукцию. Промышленное предприятие имеет доступ в течение производственного цикла к ограни-
138
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015 № 7(204). Выпуск 34/1
ченному набору технологических ресурсов, производит в единицу времени ограниченное количество продукта, определяемое производственной мощностью (эффективностью использования технологического оборудования) [1]. Конкурентноспособность предприятия в значительной степени характеризует объем выпуска продукции в единицу времени, длительность производственного цикл и объем незавершенного производства (WIP: work in process or in-process inventory) [1, 2]. В большинстве работ они приняты в качестве основных параметров модели. Предметы труда, находящиеся в незавершенном производстве, распределены вдоль технологического маршрута. Это распределение определено динамикой поступления предметов труда на первую технологическую операцию и выходом готовой продукции с последней, неравномерностью времен обработки и ограничениями на использовании ресурсов для каждой технологической операции [1,3-7]. Оптимизация параметров систем управления производственными линиями привела к появлению двух взаимосвязанных задач планирования и управления производством. Прямая задача заключается в оценке длительности производственного цикла и пропускной способности производственной системы в зависимости от количества распределенных вдоль технологического маршрута предметов труда [8, c.4591]. Обратная задача заключается в определении необходимого количества предметов труда в незавершенном производстве и его распределения вдоль технологического маршрута для обеспечения выпуска готовой продукции с заданной во времени интенсивностью. Для решения указанных задач использованы разнообразные модели управляемых производственных процессов, основные характеристики некоторых из них приведены ниже.
Дискретно-событийные модели производственного процесса (DES-модели) [1, 9].
В настоящее время при проектировании систем управления поточными линиями широко используются дискретно-событийные модели управляемого процесса (Ankenman B. E., Bekki J. M., Fowler J., 2010) [9], основанные на детальной имитации процесса переноса технологических ресурсов на предмет труда для каждой единицы оборудования поточной линии. Uzsoy R. (2011) [10], Kacar N. (2012) использовали DES-модели описания поточной линии для уточнения информации о состоянии предметов труда вдоль технологического маршрута. Погрешность модели определена факторами, связанными с выбором масштаба времени планирования. Применение DES-моделей планирования производства с разработкой детального графика загрузки производственных мощностей во времени позволило получить зависимость выпуска готовой продукции от вида функции распределения предметов труда по технологическим операциям (Lu S., [11]). Диспетчерские правила (Шкурба В.В., [7]), описывающие стратегию управления, регулируют возникновение дискретных событий начала и окончания обработки предмета труда на каждой технологической операции. Более сложные алгоритмы планирования при построении диспетчерских правил оптимизируются длительность производственного цикла и темп движения предметов труда вдоль технологического маршрута (Lasserre J. [12]). Точность соответствия DES-модели реальному производственному процессу, может быть сколь угодно высокой за счет детализации описания технологической операции и увеличения количества повторений численного эксперимента [13], теоретически ограничивается погрешностью вычислительной схемы. Затраты машинного времени, необходимого для вычисления (частота процессора 4 ГГц) миллион итераций DES-модели (Berg R., Lefeber E.,[13]) для партии 10 тыс. деталей, движущейся по технологическому маршруту, состоящему из 20 технологических операций, составляют от нескольких часов до нескольких дней [14]. В то же время, требуемая длительность интервала планирования для систем управления параметрами поточной линии современных производственных предприятий не должна превосходить несколько минут [15, с.7], (рис.1). Следующая проблема дискретного моделирования связана с тем, что процесс переноса технологических ресурсов на предмет труда в результате выполнения технологической операции является стохастическим [8,c.4589]. Случайными являются как время выполнения технологической операции, так и количество перенесенных ресурсов на предмет труда. Продолжительность вычислений, превышающая требуемый промежуток планирования, стохастичность процесса переноса технологических ресурсов на предмет труда и отсутствие конечной функциональной связи между параметрами поточной линии не позволяют эффективно использовать DES-модели в качестве инструмента управления параметрами современной производственной поточной линии. DES-модели, применяемые для описания производственных систем, например, для моделирования поточных линий по производству полупроводниковой продукции [16], являются «относительно медленными» [8,c.4589], что не позволяет их эффективно использовать для построения систем управления и планирования производством.
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015. №7 (204). Выпуск 34/1
139
проект
( MRP
{ точно-е-срок
f
поточное
________1__________I_________I__________I_______I_____
О секунды минуты дни недели месяцы
Рис.1. Длительность интервала принятия решения [15, p.7] Fig. 1. The duration of the interval decision making [15, p.7]
Модели теории массового обслуживания (TQ-модели) [1,5,17,18] широко используются для описания производственных линии в установившемся режиме. Повышение производительности вычислений достигается переходом от дискретного описания поведения предметов труда к непрерывному описанию с использованием усредненных числовых характеристик. С точки зрения использования вычислительных ресурсов модели теории массового обслуживания [17] для стационарных режимов функционирования производственной линии являются эффективнее DES-моделей, позволят представить аналитическую функциональную зависимость между потоковыми параметрами модели. Достаточно хорошо представлены TQ-модели производственных линий в работах Buzacott J.A., Shanthikumar J.G., Chao X., Hopp W.J. [18], Lefeber E.[19], которые сосредоточили свое внимание на взаимной связи между длительностью производственного цикла изготовления партии деталей и количеством предметов труда в межоперационных накопителях (рис.2). Очередь изделий, ожидающих обработку, представляет незавершенное производство (WIP) производственной системы. Достаточно точную оценку выходной траектории параметров удалось получить только для установившегося режима работы производственной линии, состоящей из небольшого количества технологических операций. Постоянное наличие переходных процессов в ходе функционирования производственной линии требует использование более совершенных нестационарных моделей, что приводит к значительному увеличению времени счета и к отсутствию наглядной функциональной связи между параметрами производственной линии [14]. При описании переходных нестационарных процессов TQ-модели теряют свои преимущества перед DES-моделями. Применение TQ-моделей для переходных режимов приводит к чрезмерному усложнению задачи. Существенным ограничением их применения является то, что современные производственные линии состоят из большого числа технологических операций. Как правило, используются одномоментные модели теории массового обслуживания. Двухмоментные модели, учитывающие дисперсию параметров обработки предметов труда, для описания производственных линий практически не встречаются в связи со сложностью их построения. Следующим ограничением для применения TQ-model является то обстоятельство, что основные формулы получены для установившихся режимов (t ), что автоматически подразумевает выполнение неравенства Я< л (рис.2), где Я - интенсивность поступления предметов труда на обработку; /л - интенсивность обработки предметов труда.
Рис. 2. Основные параметры TQ-модели [19, p.10] Fig. 2. The basic parameters TQ-model [19, p.10]
В действительности, для большинства производственных линий с конечным производственным циклом справедливо неравенство Я > /л, что не позволяет использовать хорошо разработанный аппарат теории массового обслуживания для установившихся процессов. Решение указанных про-
140
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015 № 7(204). Выпуск 34/1
блем приводит к высокой размерности задачи, громоздким вычислениям, требует значительных вычислительных ресурсов.
Модели жидкости производственных процессов (Fluid-модели) [6,20,21]. Широкий класс описания управляемых производственных процессов содержат уравнения переноса. Данный класс эффективно используется для описания переходных нестационарных, в основном детерминированных производственных процессов. Традиционно различают два подхода к представлению потока предметов труда вдоль технологического маршрута. Первый подход заключается в агрегации предметов труда по состояниям [21] и с последующим построением уравнений переноса, второй, более распространенный, в представлении потока изделий в виде непрерывной жидкости (Fuid-модель, Форрестер Дж., 1961, [6]), для описания которой используются уравнения переноса. Концептуально дискретные Fluid-модели (рис.3) представляют m-ое рабочее место в сети массового обслуживания в виде m-ого накопителя емкостью и задвижки между (т-1)-ым и m-ым накопителем, регулирующей с интенсивностью um-i движение предметов труда между ними (Kefeli, Uzsoy , Fathi, Kay, 2011) [1]:
dqm
dt
Mm—1 Mm ,
m = l.M,
M =Л,
(1)
где qm - незавершенное производство (WIP) в межоперационном заделе перед m-ым технологическом оборудованием (длина очереди), jum - темп обработки предметов труда на m-ом технологиче-
ском оборудовании, Mo = ^ - интенсивность поступление предметов труда на первую технологическую операцию. Дискретизация модели обусловлена удобством взаимосвязи уравнений модели с конкретными потоковыми параметрами для участка технологического маршрута производственной линии: межоперационными заделами в накопителе qm, темпом обработки предметом труда на те-
кущем ит и предыдущем um_i участке. Уравнения (1) Fluid-модели являются непрерывными во времени и дискретными в пространстве, описывают N пространственных частей потока предметов труда. Набор обыкновенных дифференциальных уравнений (1) представляет временную эволюцию длины очереди предметов труда в межоперационном заделе, может быть представлено в форме уравнений Форрестера (1961) [6], моделируют поток предметов труда по технологическому маршруту с дискретно расположенным технологическим оборудованием. Этапы производства рассматриваются как изменение непрерывных переменных, описывающее состояние предметов труда в технологическом пространстве состояний. В отличие от реальной жидкости, переменные состояния не описывают физическое пространство, а обозначают степень завершенности части или стадии производства. В качестве переменной, определяющей расположение оборудования, при расчете поточных
линий часто используют переменную С (метр), представляющую путь, который проходит предмет труда в результате обработки с начальной технологической операции до текущей технологической операции при общей протяженности технологического маршрута D (метр) [22]. При этом в качестве переменной, определяющей темп обработки предметов труда вдоль технологического маршрута, принимается скорость (метр/час) движения предметов труда по технологическому маршруту поточной линии (для конвейерных линий -скорость движения конвейера). Система уравнений модели (1) представляется в виде:
dq{tfm m =M(t,Cm—l) — M(tCm), m = l.M, M(0) = Л. (2)
dt
При количестве технологических операций в технологическом маршруте M>>1
(M * 102 [14], M * 250 [8], M * 300 [23,с.445]) для построения систем управления переходными производственными процессами использование Fluid-модели (1), (2) становится малоэффективным. Fluid-модели утрачивают преимущества перед TQ-моделями. Однако, факт наличия большого количества технологических операций M>>1 позволяет осуществить предельный переход при (Cm —^Cm—1 )=A << D к представлению системы уравнений (2) в виде нестационарной континуальной (непрерывной по координате пространства) модели поточной линии [8, c.4591]:
M(t,Cm ) = M(t,Cm—1) +
5u(t,С)
дС
ACm+0(л£ж2),
5p(t,C)
C=Cm
5t
5u(t, С)
С=С.
5С
-0(ACm), m = l.M, (3)
С=С
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015. №7 (204). Выпуск 34/1
141
где p(t,£m)
q(t,Cm)/ 7^m
[1,21]- усредненная плотность предметов труда в межоперационном заделе
q(t,Cm) m-ой технологической операции для участка технологического маршрута ]Cm_i,Cm],
М(0) = Л. Функция /u(t,0 задана, определяет темп работы отдельных участков поточной линии. Тот факт, что производственные процессы являются не стационарными и стохастическими, затрудняет использование Fluid-моделей в системах управления поточными линиями, требует их дальнейшего совершенствования. Другой особенностью модели является то, что уравнения (3) необходимо дополнить уравнением состояния, определяющим функцию /u(t,0. Фундаментом для постро-
ения уравнения состояния является механизм детального взаимодействия предметов труда с технологическим оборудованием. Это привело к необходимости использования для построения уравнения состояния DES-моделей со всеми описанными выше их недостатками. Применение эмпирических зависимостей, определяющих функцию /u(t,0, для нестационарных переходных процессов оказались бесперспективными из-за своей сложности построения и недостаточной точности между численными и практическими результатами [6].
Рис. 3. Концептуальное представление Fluid-модели [6] Fig. 3. Conceptual view of Fluid-model [6]
Наряду с расчетом количественного распределения предметов труда вдоль технологического маршрута при заданной во времени интенсивности поступления заказов м(0) = Л и темпе выпуска
готовой продукции m(D) = Мм представляет интерес задача определения для каждой технологиче-
ской операции уровня запаса разных видов сырья и материалов, необходимого для обеспечения бесперебойной работы поточной линии. Данная задача планирования производства в классической постановке, сформулированной Modigliani F. и Hohn F. (1955) [24], определяет шаблон для интенсивности поступления ресурсов, необходимых для выполнения технологической операции. Модель производства, предложенная Modigliani F. и Holn F., обеспечивала в дискретные моменты времени координацию действий различных частей производства с поставщиками сырья и потребителями готовой продукции. Общее потребление ресурсов на технологических операциях поточной линии связывалось с дискретными моментами времени, разделенными периодами планирования. Это позволило использовать упрощенные ограничения, которые накладывают отпечаток на поведение ресурсов в агрегированной модели оптимизации. Предложенный подход являлся доминирующим в научной литературе на протяжении десятилетий [15]. Однако, наличие проигнорированных ограничений при агрегировании сужает рамки применения моделей (Johnson L.A.,Montgomery D.C.) [25], (VoB S.,Woodruff D.L.).
Модели с использованием функции ожидания [26]. При постановке задач планирования производства возникает вопрос о масштабах времени планирования. Несмотря на то, что планирование производства осуществляется в дискретные моменты времени, планы фактически генерируются в непрерывном режиме времени в соответствии с поступившими заказами. В связи с этим Schneeweiss С. (2003) [27] предложил в качестве уравнения состояния производственной системы использовать функцию ожидания, обеспечивающую периодами времени поступления заказов и параметрами в модели планирования производства в дискретные моменты времени. В стационарных условиях ожидаемая длительность производственного цикла является нелинейной функцией использования ресурсов (Buzacott J.A., Shanthikumar J.G. [26], Hopp W.J. [18]).
Для построения нестационарных функций ожидания предложено использовать дискретные Fluid-модели. Дискретная Fluid-модель управляемого процесса «одно изделие - один технологический ресурс» представлена в виде:
I = I- + R, - Dt = I-
+X - D, X = R,,
(4)
142
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015 № 7(204). Выпуск 34/1
где I -количество готовой продукции на период планирования А/г, ti = Ц..АТ, R -количество материала поступившего в периоде Att, Dt и X, - спрос на продукцию и объем производства за период Att. Балансовое соотношение (4) соответствует простейшему типу функции ожидания (время цикла много меньше периода планирования). Предполагается, что поступающий материал для обработки в
период Ati доступен для использования в конце периода. Из-за небольшой продолжительности производственного цикла Td, Td << Ati материал перерабатывается за период Ati и незавершенным производством (WIP) можно пренебречь. Когда длительность Td превышает несколько периодов планирования, количество продукции на период планирования I i становится зависимым от количества незавершенного производства. Объем производства X за период Att. связан с количеством материала R , поступившего в период At , L = 1,2,3
Ii = Ii-\ + Xi - Di = Ii-\ + R i-L - Di
Xt = R
i i-L
(5)
Балансовое уравнение (5) широко использовалось при планирования потребностей в материалах промышленных MRP-систем (Vollmann Т.Е., 2005) [15], (Vob S., Woodruff D.L., 2003) [28]. Хотя существуют модели (Hackman S.T., Leachman R.C., 1989) [29], в которой сроки поставки технологических ресурсов равны дробному числу периодов планирования, общим, как для теории, так и для производственной практики, является предположение о том, что время поставки соответствует це-
лому числу периодов планирования. Большинство моделей вводят ограничения вида X — C на мак-
симальный выпуск X за период планирования, где C - предельная мощность производства. В кон-
це периода t производственная система имеет уровень незавершенного производства:
i i+L
W = I Rn -I xn (6)
n=i-L+\ n=i+\
Материал, поступивший в период t, остается в системе в течение интервала времени
X = (ti - ti-L). Общепринятая точка зрения на использование в дискретных Fluid-моделях MRP-систем ограничений производственных мощностей привела к задаче линейного программирования:
N
I(hi • 1г +°1 •R Ц min ,
Ii = Ii-\ + Ri-L - Di,
R-L — C, Ri > 0, It > 0,
(7)
i=\
где R -стоимость единицы незавершенного производства, стоимость единицы производственных ресурсов, используемых в момент времени t .
Hopp W.J., Spearman M.L. [18] , основываясь на детальной DES-модели взаимодействия отдельных предметов труда с оборудованием, представили в течение интервала планирования Att зависимость производительности поточной линии от интенсивности поступления технологических ресурсов, необходимых для обработки изделий. При этом Liu J., Li C., Yang F., Wang H., Uzsoy R. [10] указывают на необходимость использование для решения задачи больших ресурсов машинного времени. Для глубокого детального изучения рассмотренных проблем требуются более производительные процессоры (Kacar N.). Применение данных моделей давало хорошее совпадение теоретических и практических данных при описании квазистационарных производственных процессов. Однако их возможность описывать нелинейные зависимости между темпом движения предметов труда по технологическому маршруту и длительностью производственного цикла при интенсивно используемых технологических ресурсах вызывала сомнения [14].
Основная проблема при определении масштабов времени Ati = (ti - ti-L) заключается в том,
что система планирования и управления производством при Td >> Ati требует оценки последствия влияния решений на состояние параметров производственной системы. При использовании фиксированных промежутков At = (ti - ti-L) игнорируются эффекты, возникающие в пределах промежутка планирования. Достижения максимальной емкости накопителя межоперационных заделов в преде-
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015. №7 (204). Выпуск 34/1
143
лах промежутка планирования приводило к остановке поточной линии, что ограничивало объем производства Xt за период At,-.
Модели управляемых производственных процессов с использованием clearing-функции [30,31]. Наличие однообразных задач планирования и управления производством, для решения которых применялись различные модели или их комбинации, привело к идее создания единой оптимизационной теории производственных систем с поточным способом организации производства, для построения которой Graves S.C. (1986) [30], Karmarkar U.S. (1989) [31] предложили
использовать в качестве основных параметров состояния пропускную способность \z\cl производственной системы, объем незавершенного производства W и длительность производственного цикла T. Для описания поведение параметров системы Karmarkar U.S. ввел уравнение состоя-
ния \x\cl = Ф(^) (рис.4), задающее связь между пропускной способностью и объемом незавершенного производства, получившее название clearing-функции [31]. Clearing-функция может быть определена для группы машин, оборудования, поточных линий, одного или нескольких заводов,
включенных в единый производственный процесс. Clearing-функция \x\cl = ®(W) = const ставит фиксированную границу выпуска продукции, предполагая мгновенное наращивание мощности
производства, а \zC = 0(W) = a -W, a = const [30] предполагает фиксированное время выхода производства на полную мощность, которая при наличие ограничения выпуска продукции для промежутка планирования At,, известна как комбинированная dearing-функция (Karmarkar U.S.,1989) [31]. Важный класс составляют нелинейные dearing-функции, используемые для построения однопродуктовых моделей:
а) TQ-модель М/М/1 очереди для стационарного состояния (рис.5) [1]
\x\cl =®(W) = ^W , (8)
б) модель фундаментальной диаграммы трафика для стационарного состояния
\x\cl =®(W) = P-W - W2, (9)
в)модель G/M/1 очереди для стационарного состояния (Mehdi J., 1991) [1], (Berg R.,2004) [13,c.7]:
P1 Я
~ + P, P = ~<1 (10)
1-P P
2 2 w = Ca + Cs
2
где c 2 и c2 характеризуют среднеквадратичное отклонение поступления требований на изготовления изделий и времени их обработки, р - темп обработки предметов труда, Я - интенсивность поступление предметов труда на первую технологическую операцию [1]. Модель G/M/1 очереди для стационарного состояния (10) является развитием модели (8), движение предметов труда по технологическому маршруту с последовательным расположением технологического оборудования. Berg R. указал на то, что устойчивое состояние для модели (10) обеспечивается при выполнении условия
Р< 1 [13,с.6]. При Я^р межоперационные заделы бесконечно большие (W ), а при Я> р уравнение (10) не может быть использовано, т.к. при его выводе предполагаось р = (Я/р)< 1. Решение
(10) относительно р< 1 описано в [1]:
Р = -
yl(W +1)2 + 4W(c2 -1) - (W +1)
при
c =
C 2 + c 2
2(c2 -1)
При c ^ 1 модель M/M/1 очереди (8) является предельным случаем модели (10):
limp = W lim[^\CL = P'W
c^i (W +1) ’ cJrC 1 + W
(11)
(12)
2
144
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015 № 7(204). Выпуск 34/1
Рис. 4. Clearing-функция производственной системы [31] Fig. 4. Clearing-function production system [31]
л __ л_________ л
— 11 Г га- 1 Г д 4 1 ]
ц ^ р fi ^ ц J р
Рис. 5. Схема М/М/1-модели очереди поточной линии
. Схема М/1
Fig. 5. Model M/M/1 queue production line
Asmundsson J.M (2006) [32] предложил распределенные АС-функции (allocated clearing function) для моделирования многопродуктовых линий. AC-функции предполагают, что производится совокупный продукт, агрегирующий в себе технологические ресурсы для производства отдельных продуктов. Альтернативном подходом является представление dearing-функции в виде суммы dearing-функций отдельных продуктов. Экспериментальные данные свидетельствуют об удовлетворительных результатах использования AC-функций для продуктов, схожих по характеру потребления ресурсов. Интенсивности потребления ресурсов при этом выражаются в единицах времени обработки [1]. Однако, если модели переноса ресурсов для каждой номенклатуры предметов труда являются сложными, использование AC-функций не позволяет описать производственные процессы. Selguk B., Fransoo J.C., Gok A.G. (2007) аппроксимировали dearing-функцию кусочнолинейной функцией [1], что позволило использовать для оптимизации параметров поточной линии аппарат линейного программирования. Clearing-функция может быть получена как аналитически так и численно c использованием TQ-моделей, DES-моделей, Fluid-моделей производственных систем или определена эмпирически. Selguk, B., Fransoo J.C.,Gok A.G. (2007) [1] представили методику построения переходных dearing-функций аналитически. В связи с тем, что оперативная информация о выпуске продукции и о состоянии незавершенного производства на предприятии является закрытой, то в большинстве исследований для построения dearing-функции вместо эмпирических данных используется TQ- и DES-модели. Исключением являются работы Haeussler S., Missbauer H. (2012) [1], в которых для построения dearing-функции применены полученные с производственной линии цифровые носители. Kacar N. (2012) использовал для построении dearing-функции оптимизированные параметры поточной линии. Вычислительные эксперименты, связанные с расчетом параметров поточных линий компании Intel, показали хорошее приближение расчетных и экспериментальных данных для установившихся процессов [1].
Несмотря на то, что dearing-функция являются удачным инструментом для определения мгновенной связи между пропускной способностью производственной системы и объемом незавершенного производства, наличие ограниченного количества параметров, входящих в уравнение состояния, не позволяет эффективно моделировать временное изменение параметров производственного процесса, обусловленное факторами обработкой изделий (Armbruster D, Kempf K., 2012) [1]. Попытки создания нестационарных dearing-функций ограничены специальными теоретическими уточнениями и экспериментальными исследованиями (Fonteijn J., Wienke M., 2012) [1]. В качестве уточнения Lefeber Е. (2008) [20] вводит в dearing-функцию \x\cl = ®(W(t -г0)) эффективное время м
обработки т0 = - (основное технологическое время) предметов труда на технологических
т=0
операциях [13]. Missbauer H (2009) предпринял попытки расширить применение dearing-функции на переходные производственные процессы. При этом обращалось внимание на существенную зависимость пропускной способности производственной системы от начального распределения предметов труда по технологическому маршруту и необходимость обеспечения условий перехода производственной системы из одного стационарного устойчивого состояния в другое. Производственные процессы являются стохастичными [21,33], но несмотря на это при построении dearing-функций практически не уделено внимание исследованию устойчивости потоковых параметров производ-
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015. №7 (204). Выпуск 34/1
145
ственных линии, отсутствуют оценки времени затухания случайных возмущении потоковых параметров и оценки их абсолютных величин. Исследования поточной линии Intel по производству полупроводниковых изделий, проведенные Armbruster D., Kempf K.G. (2012) [1] показали, что повседневные стохастичные факторы возмущения параметров производства имеют время затухания от 12 дней до недели, что требует наличие страховых запасов до 20% от нормативного количества. Предположение о квазистационарности переходного процесса является значительным ограничением для широкого применения уравнения состояния (dearing-функции) в аналитическом виде, построенного в основном с использованием TQ-моделей.
Модели производственных процессов, использующие уравнения в частных производных. PDE-модели [8]. В современной литературе выделились основные три типа моделей и их комбинации для вывода уравнения состояния, определяющего связь потоковых параметров производственных линий. Это модели массового обслуживания (TQ-model), дискретно-событийные модели (DES-model) и Fluid-модели [1]. Каждый тип моделей имеет свои преимущества, но ни один из них не подходит в полной мере для моделирования как установившихся, так и переходных процессов функционирования производственной системы [13, с.2]. Существующие TQ-модели описывают поточные линии в стационарном режиме [17]. Использование их при описании переходных процессов приводит к чрезмерному усложнению и большим затратам машинного времени. DES-модели используются для описания поточных линий в переходном и стационарном режиме, но являются дискретными и требует больших затрат машинного времени. Fluid-модели ориентированы на малое количество интервалов разбиения технологического маршрута и линейные стационарные решения в рамках заданного интервала. Требование повышения точности модели приводит к увеличению количества интервалов обобщения и к усложнению модели в связи с увеличением размерности системы дифференциальных уравнений (1) (Kefeli A., Uzsoy R., 2011).
В последнее десятилетие при проектировании производственных поточных линий используются модели, описывающие поведение производственной системы с помощью уравнений в частных производных (PDE-model) [1,9,13,20,21]. Введенный класс моделей объединил преимущества TQ-моделей, DES-моделей и Fluid-моделей, значительно расширил возможности проектирования систем управления поточными линиями. PDE-модели в общем случае являются непрерывными, могут быть успешно использованы при описании стационарных и переходных режимов работы поточной линии, не требуют больших затрат машинного времени [1].
Ключевым вопросом при построении PDE-модели поточных линии является выбор системы координат. Распространенным подходом является использование в качестве переменной, определяющей место обработки предмета труда на технологическом маршруте, стоимости S (грн.) перенесенных технологических ресурсов на предмет труда (Дабагян А.В.,2008) [1], (Федюкин В.К., 2004) S е[0, S ] (Sd (грн)- себестоимость изготовления продукции), эффективного времени обработки
M
предмета труда тт (час), тт е[0,гм] (Eekelen J., 2006), (Ramadge P., Wonham W.) [1] (тм = ^Лтт
m=1
(час)-общее эффективное время обработки предмета труда, Лтт -среднее время обработки предмета труда на m-ой технологической операции ) или степени незавершенности изготовления изделия X (Armbruster D.,Ringhofer С., Berg V., Lefeber E.,2004) [1], X e [0.1] . Под степенью незавершенности изготовления изделия X понимается позиция предмета труда в технологическом маршруте, которая может быть представлена как отношение среднего времени Лтт обработки предмета труда к его общему времени обработки тм [13,с.1б]. Для предмета труда, прошедшего обработку на m-ой опера-
Т f m / M \ m
ции, можно записать х = —m = 1 "^Лтк Г^Лтк I. Каждому времени обработки тт=^Лтк одно-
Тм V *=1 / *=1 ) k=1
значно соответствует стоимость ресурсов Sm = S(тт) , перенесенных на предмет труда, а общему времени тм - себестоимость S = S(гм). Таким образом, степень незавершенности изготовления изделия X может быть определена через время обработки т или стоимость перенесенных затрат S = S(т) на предмет труда. Целесообразно при моделировании производственной поточной линии с
одним обобщенным технологическим ресурсом использовать безразмерную переменную х =
S
определяющую позицию предмета труда в технологическом маршруте [3], х e [0.1] . Безразмерная т
переменная х =
используется в случае, когда модель поточной линии не рассматривает при
146
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015 № 7(204). Выпуск 34/1
переноса на предмет труда структуру ресурсов (Armbruster D., Ringhofer С.) [1]. Для модели поточной линии, которая учитывает в результате выполнения технологической операции потребление несколько взаимосвязанных между собой технологических ресурсов, на использование которых наложены ограничения, применение безразмерной переменной X является затруднительным. Если ввести функцию плотности предметов труда p(t, x) в состоянии X в момент времени t, то общее количество предметов труда, находящихся в различных стадиях готовности есть величина (Armbruster D., Ringhofer С.) [1,20,21]:
W(t) = J p(t, x)dx, x e [0,l] (13)
0
Так как процесс обработки предметов труда является стохастическим, то в результате выполнения операции предмет труда может находиться в том или ином состоянии [8,c.4544]. Это позволило записать усредненную плотность предметов труда p(t, x) (шт.) и поток предметов труда
F(t, x) (шт./час) по технологическому маршруту через функцию распределения предметов труда
f (t, r, x) по состояниям r = Лтт1 (Armbruster D., Ringhofer С., 2005) [1]:
p(t, x) = J f (t,r, x)dr, F(t, x) = J- f (t,r, x)dr (14)
0 0r
Положение предмета труда в пространстве состояний характеризуется точкой с координатами (q,q,..,q ,..,qn), которые определяют количественное значение параметров предмета труда.
Пространство состояний используется при построения моделей многопоточных линий, потребляющих в ходе производства несколько ресурсов. Использование в одномерном описании безразмерной переменной x [1, 3] для углубленного исследования изменения состояния предмета труда является затруднительным. Напротив, применение в качестве переменной модели параметра S, характеризующего технологическую позицию обработки предмета труда через стоимость позволяет применить детально разработанный математический аппарат производственных функций [21], что дает возможность суммировать ресурсы путем суммирования их стоимостей. Однако, несмотря на раскрывающиеся перспективы в связи с использованием стоимостного представления изменения состояния предмета труда, подавляющее большинство авторов (Armbruster D., Ringhofer С. (2005), Berg R., Lefeber E., Rooda J. (2008) [13], Wienke M., Fonteijn J., (2012), Kempf K., (2012) [16]), применяют PDE-модели поточных линий без ограничений на потребление технологических ресурсов, вводят для описании состояния предмета труда переменные состояния (r, x) .
В PDE-модели поток предметов труда F(t, x) = pt, x) • v(t, x) (пропускная способность поточной линии) представляется в виде произведения плотности pt, x) предметов труда и скорости v(t, x) их движения [33]. Предполагая, что дефектные предметы труда отсутствуют (нет источников и стоков), движение предметов труда по технологическому маршруту удовлетворяет закону сохранения:
dQ = Я-р, Q(t) = W(t) = Jp(t, x)dx (15)
dt 0
где Я - интенсивность поступление предметов труда на первую технологическую операцию, р -
выпуск продукции [1]. Уравнение (15) является интегральной формой закона сохранения количества предметов труда, находящихся в процессе обработки, может быть представлено в дифференциальном виде [1]:
dp(t, x) dF (t, x)
dt + dS
F(t, x) = pt, x) • v(t, x).
(16)
Граничное условие F (t,0) = Я() задает поток заготовок на первую операцию. Профиль не-
завершенного производства в начальный момент времени определяется условием p(0, S) = p0 (S), которое характеризует распределение заготовок по технологическим операциям поточной линии. Для свободной поточной линии p(0, S) =0. Уравнение (16) обеспечивает взаимосвязь во времени плотности pt, x) распределения предметов труда и темпа их движения F(t, x) (пропускной способности поточной линии) для каждой точки x технологического маршрута [21]. Неравномерность распреде-
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015. №7 (204). Выпуск 34/1
147
ления незавершенного производства вдоль маршрута обусловлена разными эффективными временами обработки предметов труда на каждой операции. Неравномерность в производительности оборудования вдоль технологического маршрута определяет динамику изменения плотности предметов труда, что существенным образом влияет на пропускную способность линии. Модели производственного процесса, содержащие уравнение (16) учитывают влияние внутренних факторов на пропускную способность, а также ограничения, определяемые максимальной производительностью оборудования и емкостью межоперационных накопителей. Это позволило PDE-моделям составить значительную конкуренцию DES-моделям, преимущество которых в том, что они допускают решение в аналитическом виде и не требуют значительных вычислительных ресурсов. Трудность построения PDE-моделей определяется тем, что балансовые уравнения вида (16) незамкнуты [13, 21, 33]. Для замыкания уравнение (16) дополняется уравнением состояния. Если уравнение состояния зада-
i i
но линейной dearing-функцией Ф(Щ = JF(t,x)dx = cjp(t,x)dx = с-W, (Constant Proportion, Graves S.C.,
0 0
1986, [30]) (рис.4), то система уравнений PDE-модели:
а^х) + dF(t х) =0, ^(t, x) = с = const, F(t, x) = p(t, x) - c (17)
ot OS
допускает аналитическое решение. Если скорость движения предметов труда вдоль технологического маршрута постоянна v(t, x) = c , то для единичного потока предметов труда, поступающего в момент времени t=0 на первую технологическую операцию поточной линии, уравнения (17) имеет решение:
f 0, t< V,
c-pit, x) = H(c -1 - x) = \ /c 2(t) = H(t). (18)
1 1 e %,
Встроенная в модель величина т = c1 определяет длительность задержки между временем поступления сырья на первую операцию и временем выхода готового продукта. Постоянной скорости движения предметов труда по технологическому маршруту соответствует постоянное время задержки т. Связь потока предметов труда представлена линейной зависимостью Следующая PDE-модель (Lighthill-Whitham) [1] использует уравнения состояния в виде нелинейной зависимость потока F (p) предметов труда от плотности p :
0p(t, x) OF (t, x)
dt + =0, F(P) = P-v(P)
p-v0(1 -p) ,
R
R<-1. a«
Широкое распространение получили PDE-модели поточных линий, содержащие стационарные уравнения состояния. Для завода, уравнение состояния поточной линии которого представлено М/М/1-моделью очереди размером W(t), параметры 2, р, W (15) и длительность производственного цикла T = (l+ W)/р для стационарного состояния определяются зависимостью (Gross D., Harris C.) [17]:
2 =
p-W
1 + W
где
Td =-
W = -
2
W (t) =Jp(t, x)dx. (20)
р-2 p-2 0
Скорость движения предметов труда v(t, x) и локальный поток F(t, x) для М/М/1- модели выражаются через интенсивность выпуска продукции р [13]
0p(t, x) OF (t, x)
v(t, x)'
P
F(t, x) = p(t, x) - v(t, x):
P-P(t, x)
(21)
Ot dS 1 + W 1 + W
Уравнение состояния F(t, x) в интегральной форме для устойчивого режима работы поточ ной линии (Armbruster D., Fonteijn J., Wienke M., 2012) имеет вид [1]
F (t,W) = j F (t, x)dx = P
p-W
+ W
(22)
М/М/1-модели dearing-функции (6). Предприняты попытки построить совершенные PDE-модели поточных линий, уточняющие уравнения состояния путем использования экспериментальных измерений на заводе, детальной имитационной модели или TQ-моделей [1]. Широко применяются при исследовании поточных линий модели:
148
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015 № 7(204). Выпуск 34/1
F(t, x) = Ф(р)
Fo
F (t, x) =
1 + p(t, x)
F (t, x) = k -(l - ehw) 1 - x
pt, x), F (t, x) =
Fo
1+p(t, x)+k - p(t, x)(l - x) ’
, k-w
F (t, x) = P
k + w
(i
p(t,x),
r(t, x) =
T(t, x)
I p(t, z)dz - x -т0
1
F
(23)
(24)
(25)
V 0
которые предложили D.Armbruster, K.Kempf (2012) (23), J.Asmundsson., R.Uzsoy (24), Ringhofer С (2012) (25). r(t, x) - время, требуемое для завершения производства предметов труда, находящихся
на стадии x; k , k, k2 -технологические коэффициенты; M-характеристика максимальной емкость
накопителя.
Детальный анализ функционирования полупроводниковой поточной линии с применением PDE-модели (сплошная линия) и DES-модели (затенение) выполнили Perdaen D. (2008), Lefeber E. (2010) [1]. Характерное поведение выходного потока предметов труда, рассчитанного с помощью PDE-модели и DES-модели, представлено на рис.6 (Lefeber E., 2010). Принципиальным недостатком методов, использующих DES-модели, является чрезвычайно большое количество вычислений для сложных производственных систем, таких как полупроводниковые поточные линии [14]. Именно в этой области применение PDE-моделей предлагают самые высокие перспективы для проектирования систем управления поточными линиями. Эти методы способны оптимально сочетать точность DES-моделей при применении гораздо менее производительных процессоров [13]. Использование в PDE-моделях dearing-функции проявляло значительные перспективы на ранней стадии их развития. Растущий объем публикаций, связанных с развитием и уточнением уравнения состояния для PDE-моделей свидетельствует о том, что дальнейшее развитие подхода с использованием clearing-функций является не эффективным. Необходимо развивать статистические методы, позволяющие строить многомоментные PDE-модели для переходных режимов, в которых для замыкания использовано нестационарное уравнение, в основанное на детальном механизме взаимодействия предметов труда между собой и оборудованием.
Рис. 6. Сравнение решений PDE-модели и DES-модели, Lefeber E. 2010 Fig. 6. Comparison of PDE-model and DES-model, Lefeber E. 2010
Кинетические модели управляемых производственных процессов [3, 21, 34]. В серии работ Armbruster D , Ringhofer С., Lefeber E., Kempf K. [1, 9, 16] представили кинетические модели поточных линий. Armbruster D., Ringhofer С. (2004) вводят функции распределения f (x,v, t) предметов труда по состояниям, характеризующую количество деталей в состоянии x в момент времени t. Типичный подход, с помощью которого определяется эволюция функции распределения предметов труда по состояниям, заключается в выводе замкнутых уравнений для моментов функции распределения. Построение PDE-моделей с использованием кинетической теории, содержит иерархический набор уравнений. Это позволяет выйти за пределы применимости квазистатистических моделей. Предложенный Н.Н.Боголюбовым метод, основанный на выборе малого параметра, позволяет обрезать количество уравнений на нужном уровне. С использование кинетического подхода записаны уравнения для первых моментов функции распре-деления предметов труда f (x, v,t) (Armbruster D., Ringhofer С.) [35, с.819]:
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015. №7 (204). Выпуск 34/1
149
dp(t, x) dp(t, x) • v(t, x)
+ o,
dt dx
p(t,0) • v(t,0) = X(t),
dv(t, x) . dv(t, x)
v ■+ v(t, x) • . 7 =o,
v(t,1) =
dt
M
1 + W (t)
dv(t,0) dt
(
= —a
v(t,0) -
M
1 + W (t)
, при X<p , v(t ,0) =
M
dx
, при Л> m .
(26)
(27)
(28)
0,5 + W (t)
с граничными условиями (27) M/M/1-модели, полученной для установившегося устойчивого режима, где а -экспериментальная величина. Интегрирование первого уравнения (26) по x позволяет записать уравнение:
dW(t) = p(t,0) • v(t,0) - p(t,1) • v(t,1) = X(t) - X\ t - -dt ( c
(1.29)
которое получил Lefeber Е. (2008) [1] из весьма общих соображений. Известно из практических исследований, что выход с обработки первого изделия партии происходит через некоторое время задержки относительно времени поступления партии на обработку [21]. При переходе от p(t, x) к агрегированной переменной Fluid-модели W (t) (15) влияние неравномерного распределения предметов труда вдоль технологического маршрута [5] и наличие ограничений емкости накопителей не учитывается. Результаты расчетов потоковых параметров поточной линии, полученные с использованием модели (26)-(28), более приближены к экспериментальным данным, чем результаты расчетов с использованием М/М/1-модели (8), G/M/1-модели (1.10), непрерывной Fluid-модели (15) (Fonteijn J., Missbauer H.), хотя, как утверждает Armbruster D. (2012), детальное изучение приближений, связанных уравнениями (26)-(28), на текущий момент не сделано. Не ясным остался вопрос, для каких производственных систем та или иная модель является наиболее удачной. Наиболее существенные проблемы, возникающие при построении кинетических моделей производственных систем, заключаются в том, что кинетическое уравнение является интегро-дифференциальным, решение которого представляется сложной математической задачей [34]. В виду сложности закона воздействия оборудования на предмет труда, кинетическое уравнение не может быть записано в точном виде для конкретных производственных процессов. Даже при простых предположениях о характере воздействия оборудования на предмет труда не удается получить решения в точном аналитическом виде. В связи с этим приобретают особое значение построение эффективных методов приближенного решения кинетического уравнения производственной линии.
Выводы
Представлен обзор основных моделей, которые использованы при проектировании систем управления поточными производственными линиями. Рассмотрены области применения моделей и ограничения, которые не позволяют их эффективно использовать для проектирования систем управления. Уделено значительное внимание новым типам моделей - кинетическим и моделям, содержащим уравнения в частных производных (PDE-моделям). Проведен анализ использования моделей для наиболее простых случаев функционирования производственных поточных линий. Обоснованность применения определялась сравнительным анализом результатов, полученных с помощью DES-модели и исследуемой PDE-моделью. Показано, как для построения PDE-моделей использованы статистические методы описания больших систем. При этом общий характер статистических закономерностей не зависит от того, каким способом описывается поведение отдельного предмета труда. Использование статистического подхода позволяет получить замкнутые многомо-ментные балансовые уравнения (уравнения переноса) не из феноменологических соображений, а исходя из законов движения отдельных предметов труда по технологическому маршруту, определенных технологией производства. Развитие и использование PDE-моделей требует решения вопросов: 1.Вывод нестационарных уравнений состояния, основанных на детальной технологии обработки предмета труда с учетом схемы оборудования. 2.Построение многомоментных замкнутых балансовых моделей для установившихся и переходных нестационарных режимов функционирования производственной линии. 3.Построение двухуровневых моделей управления параметрами производственной линии для установившихся и переходных режимов с учетом параметров оборудования, схемы его расстановки и приоритетов движения предметов труда
150
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015 № 7(204). Выпуск 34/1
Список литературы
References
1. Пигнастый О. М. О новом классе динамических моделей поточных линий производственных систем / О. М. Пигнастый // Научные ведомости Белгородского государственного университета. Белгород: БГУ. -2014. - № 31/1. - С. 147-157.
Pignastyj O. M. O novom klasse dinamicheskih modelej potochnyh linij proizvodstvennyh si-stem /
O. M. Pignastyj // Nauchnye vedomosti Belgorodskogo gosudarstvennogo universiteta. Belgorod: BGU. - 2014. -№ 31/1. - S. 147-157.
2. Демуцкий В.П. Теория предприятия: Устойчивость функционирования массового производства и продвижения продукции на рынок. / Демуцкий В. П., Пигнастая В. С., Пигнастый О. М. - Х.: ХНУ, 2003. -C. 272.
Demuckij V.P. Teorija predprijatija: Ustojchivost' funkcionirovanija massovogo proizvodstva i prodvizhenija produkcii na rynok. / Demuckij V. P., Pignastaja V. S., Pignastyj O. M. - H.: HNU, 2003. - C. 272.
3. Armbruster D. Kinetic and fluid model hierarchies for supply chains supporting policy attributes / D. Armbruster., D. Marthaler, C. Ringhofer // Bulletin of the Institute of Mathematics. - Academica Sinica. -2006. - P. 496 - 521.
Armbruster D. Kinetic and fluid model hierarchies for supply chains supporting policy attributes / D. Arm-bruster., D. Marthaler, C. Ringhofer // Bulletin of the Institute of Mathematics. - Academica Sinica. - 2006. -
P. 496 - 521.
4. Бусленко Н. П. Математическое моделирование производственных процессов / Н. П. Бусленко. - М.: Наука, 1964. - 363 с.
Buslenko N. P. Matematicheskoe modelirovanie proizvodstvennyh processov / N. P. Buslenko. - M.: Nauka, 1964. - 363 s.
5. Первозванский А. А. Математические методы в управлении производством / А. А. Первозванский. -М.: Наука, 1975. - 616 с.
Pervozvanskij A. A. Matematicheskie metody v upravlenii proizvodstvom / A. A. Pervozvanskij. - M.: Nauka, 1975. - 616 s.
6. Форрестер Дж. Основы кибернетики предприятия / Дж. Форрестер. - М.: Прогресс, 1961. - 341 с. Forrester Dzh. Osnovy kibernetiki predprijatija / Dzh. Forrester. - M.: Progress, 1961. - 341 s.
7. Шкурба В. В. Планирование дискретного производства в условиях АСУ / В .В. Шкурба, В .А. Болдырева, A. A. Вьюн и др. / под ред. В. М. Глушкова. - К.: Техника, 1975. - 296 с.
Shkurba V. V. Planirovanie diskretnogo proizvodstva v uslovijah ASU / V .V. Shkurba, V .A. Boldyreva,
A. A. V'jun i dr. / pod red. V. M. Glushkova. - K.: Tehnika, 1975. - 296 s.
8. Armbruster D. Continuous models for production flows.. / D. Armbruster, C. Ringhofer., Jo T- J. // In Proceedings of the 2004 American Control Conference. - Boston, MA, USA. - 2004. - P. 4589 - 4594.
Armbruster D. Continuous models for production flows.. / D. Armbruster, C. Ringhofer., Jo T- J. // In Proceedings of the 2004 American Control Conference. - Boston, MA, USA. - 2004. - P. 4589 - 4594.
9. Ankenman B. E. Simulation in planning: Planning Production and Inventories in the Extended Enterprise. /
B. E. Ankenman, J. M. Bekki, J. Fowler et al. // A State of the Art Handbook New York: Springer-Verlag. - 2010. -Vol. 151. - P. 565-592.
Ankenman B. E. Simulation in planning: Planning Production and Inventories in the Extended Enter-prise. / B. E. Ankenman, J. M. Bekki, J. Fowler et al. // A State of the Art Handbook New York: Springer-Verlag. - 2010. -Vol. 151. - P. 565-592.
10. Liu J. Production planning for semiconductor manufacturing via simulation optimization / J.Liu, C.Li, F.Yang, R.Uzsoy, S.Jain, R .Creasey/ / Simulation Conference (WSC) - New York: IEEE. -2011. - P. 256 - 322.
Liu J. Production planning for semiconductor manufacturing via simulation optimization / J.Liu, C.Li, F.Yang,
R. Uzsoy, S.Jain, R .Creasey/ / Simulation Conference (WSC) - New York: IEEE. -2011. - P. 256 - 322.
11. Lu S.R. Effi cient scheduling policies to reduce mean and variance of cycle time in semiconductor plants / Lu S., Ramaswamy D., Kumar P.R. // IEEE Transactions on Semiconductor Manufacturing - 1994. №7(3). -P. 374 - 388.
Lu S.R. Efficient scheduling policies to reduce mean and variance of cycle time in semiconductor plants / Lu
S. , Ramaswamy D., Kumar P.R. // IEEE Transactions on Semiconductor Manufacturing - 1994. №7(3). -P. 374 - 388.
12. Dauzere-Peres S. On the importance of sequencing decisions in production planning and scheduling. / S.Dauzere-Peres, J. Lasserre // Journal International Transactions in Operational Research. - 2006. - №9. -P. 779 - 793.
Dauzere-Peres S. On the importance of sequencing decisions in production planning and scheduling. / S.Dauzere-Peres, J. Lasserre / / Journal International Transactions in Operational Research. - 2006. - №9. -P. 779 - 793.
13. Berg R. Partial differential equations in modelling and control of manufacturing systems / R. Berg. - Netherlands, Eindhoven Univ. Technol., 2004. - 157 p.
Berg R. Partial differential equations in modelling and control of manufacturing systems / R. Berg. - Netherlands, Eindhoven Univ. Technol., 2004. - 157 p.
14. Lefeber E. Modeling, Validation and Control of Manufacturing Systems / E.Lefeber, R.A.Berg, J.E. Rooda // - Proceeding of the 2004 American Control Conference. -Massa- chusetts. - 2004. - P. 4583 - 4588.
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015. №7 (204). Выпуск 34/1
151
Lefeber E. Modeling, Validation and Control of Manufacturing Systems / E.Lefeber, R.A.Berg, J.E. Rooda // -Proceeding of the 2004 American Control Conference. -Massa- chusetts. - 2004. - P. 4583 - 4588.
15. Vollmann T.E. Manufacturing Planning and Control for Supply Chain Management. / T.E.Vollmann, L.Berry, F.RJacobs // New York: McGraw-HilL-2005. - P. 520.
Vollmann T.E. Manufacturing Planning and Control for Supply Chain Management. / T.E.Vollmann, L.Berry, F.R.Jacobs // New York: McGraw-HilL-2005. - P. 520.
16. Kempf K.G. Simulating Semiconductor Manufacturing Systems / K.G. Kempf // Winter Simulation Conference, Institute of Electrical and Electronics Engineers - Piscataway, New Jersey, - 1996. - P. 3 - 11.
Kempf K.G. Simulating Semiconductor Manufacturing Systems / K.G. Kempf // Winter Simulation Conference, Institute of Electrical and Electronics Engineers - Piscataway, New Jersey, - 1996. - P. 3 - 11.
17. Gross D. Fundamentals of Queueing Theory. / D.Gross, C.M.Harris. - New York: Wiley, 1985. - 587 p.
Gross D. Fundamentals of Queueing Theory. / D.Gross, C.M.Harris. - New York: Wiley, 1985. - 587 p.
18. Hopp W.J. Factory Physics: Foundations of Manufacturing Management. / W. J. Hopp, M.L. Spearman. -Boston: Irwin. McGraw-Hill, 2001. - P. 698.
Hopp W.J. Factory Physics: Foundations of Manufacturing Management / W. J. Hopp, M.L. Spearman. -Boston: Irwin. McGraw-Hill, 2001. - P. 698.
19. Lefeber E. Modeling and Control of Manufacturing Systems / E.Lefeber // Decision Policies for Production Networks. - 2012.- Р.9-30.
Lefeber E. Modeling and Control of Manufacturing Systems / E.Lefeber // Decision Policies for Produc-tion Networks. - 2012.- R.9-30.
20. Пигнастый О. М. Континуальш моделi прогнозування виробничного функщонування поточних лшш / О. М. Пигнастый, В. Я. Заруба // Сучаст концепци прогнозування розвитку складних сощально-eкономiчних систем. - Бердянск: Видавець Ткачук О.В. - 2013. - С. 74-89.
Pignastyj O. M. Kontinual'ni modeli prognozuvannja virobnichnogo funkcionuvannja potochnih linij / O. M. Pignastyj, V. Ja. Zaruba // Suchasni koncepcii prognozuvannja rozvitku skladnih social'no-ekonomichnih sis-tem. - Berdjansk: Vidavec' Tkachuk O.V. - 2013. - S. 74-89.
21. Пигнастый О. М. Статистическая теория производственных систем / О. М Пигнастый. - Харюв: ХНУ, 2007. - 388 с.
Pignastyj O. M. Statisticheskaja teorija proizvodstvennyh sistem / O. M Pignastyj. - Harkiv: HNU, 2007. -
388 s.
22. Соколицын С. А. Применение математических методов в экономике и организации машиностроительного производства / С. А. Соколицын. - Л.: Машиностроение, 1970. - 345 с.
Sokolicyn S. A. Primenenie matematicheskih metodov v jekonomike i organizacii mashinostro-itel'nogo pro-izvodstva / S. A. Sokolicyn. - L.: Mashinostroenie, 1970. - 345 s.
23. Tian F. An iterative approach to item-level tactical production and inventory planning. / F.Tian, S.P.Willems, K.G.Kempf // International Journal of Production Economics. - 2011. - №133. - P. 439 - 450.
Tian F. An iterative approach to item-level tactical production and inventory planning. / F.Tian, S.P.Willems, K.G.Kempf / / International Journal of Production Economics. - 2011. - №133. - P. 439 - 450.
24. Modigliani F. Production planning over time and the nature of the expectation and planning / F.Modigliani, F.Hohn// Econometrica №23. - New York, 1955. - P. 46 - 66.
Modigliani F. Production planning over time and the nature of the expectation and planning / F.Modigliani, F.Hohn// Econometrica №23. - New York, 1955. - P. 46 - 66.
25. Johnson L.A. Operations Research in Production Planning, Scheduling and Inventory Control / L.A. Johnson, D.C. Montgomery. - New York: Wiley, 1974. - 544 р.
Johnson L.A. Operations Research in Production Planning, Scheduling and Inventory Control / L.A. Johnson, D.C. Montgomery. - New York: Wiley, 1974. - 544 p.
26. Buzacott J. A. Stochastic Models of Manufacturing Systems. / J. A. Buzacott, J. G. Shanthikumar. - Prentice-Hall, Englewood Cliffs, NJ, 1993. - 403 p.
Buzacott J. A. Stochastic Models of Manufacturing Systems. / J. A. Buzacott, J. G. Shanthikumar. - Pren-tice-Hall, Englewood Cliffs, NJ, 1993. - 403 p.
27. Schneeweiss C. Distributed Decision Making. / C. Schneeweiss. - Berlin: Springer-Verlag, 2003. -508 P.
Schneeweiss C. Distributed Decision Making. / C. Schneeweiss. - Berlin: Springer-Verlag, 2003. -508 P.
28. Vob S. Introduction to Computational Optimization Models for Production Planning in a Supply Chain / S. Vob, D. L. Woodruff // Berlin: Springer - Verlag. -2003. - 261 р.
Vob S. Introduction to Computational Optimization Models for Production Planning in a Supply Chain / S. Vob, D. L. Woodruff // Berlin: Springer - Verlag. -2003. - 261 р.
29. Hackman S.T. A general framework for modeling production / S.T.Hackman, R.C. Leachman // Management Science. - 1989. - №4. - P. 478 - 495.
Hackman S.T. A general framework for modeling production / S.T.Hackman, R.C. Leachman // Management Science. - 1989. - №4. - P. 478 - 495.
30. Graves S.C. A tactical planning model for a job shop. Operations Research 34 (4). /S.C. Graves - New York,1986. - P. 522 - 533.
Graves S.C. A tactical planning model for a job shop. Operations Research 34 (4). /S.C. Graves - New York,1986. - P. 522 - 533.
31. Karmarkar U.S. Capacity Loading and Release Planning with Work-in-Progress (WIP) and Leadtimes / U.S. Karmarkar / / Journal of Manufacturing and Operations Management. -№2. -1989. - P.105-123.
152
НАУЧНЫЕ ВЕДОМОСТИ
Серия Экономика. Информатика. 2015 № 7(204). Выпуск 34/1
Karmarkar U.S. Capacity Loading and Release Planning with Work-in-Progress (WIP) and Leadtimes /
U. S. Karmarkar / / Journal of Manufacturing and Operations Management. -№2. -1989. - P.105- 123.
32. Asmundsson J. M. Production planning models with resources subject to congestion. / J. M Asmundsson, R. L Rardin, C. H. Turkseven, R. Uzsoy // Naval Res Logist. - 2009. - N 56 - P. 142 - 179.
Asmundsson J. M. Production planning models with resources subject to congestion. / J. M Asmundsson, R. L Rardin, C. H. Turkseven, R. Uzsoy // Naval Res Logist. - 2009. - N 56 - P. 142 - 179.
33. Пигнастый О. М. Стохастическое описание экономико-производственных систем с массовым выпуском продукции / В. П. Демуцкий, В. С. Пигнастая, О. М. Пигнастый // Доповда Нащонально1 академп наук Украши. - Ки1в: Видавничий д1м "Академперюдика''. - 2005. - № 7. - С. 66-71.
Pignastyj O. M. Stohasticheskoe opisanie jekonomiko-proizvodstvennyh sistem s massovym vypuskom produkcii / V. P. Demuckij, V. S. Pignastaja, O. M. Pignastyj // Dopovidi Nacional'noi akademii nauk Ukraini. - Kiiv: Vidavnichij dim "Akademperiodika". - 2005. - № 7. - S. 66-71.
34. Азаренков Н.А. Кинетическая теория колебаний параметров поточной линии / Н. А. Азаренков, О. М. Пигнастый, В. Д. Ходусов // Доповда Нащонально1 академп наук Украши. 2014. № 12. - С. 36 -43.
Azarenkov N.A. Kineticheskaja teorija kolebanij parametrov potochnoj linii / N. A. Azarenkov, O. M. Pignastyj,
V. D. Hodusov // Dopovidi Nacional'noi akademii nauk Ukraini. - 2014. - № 12. - - S. 36 -43.
35. Armbruster D. Modeling production planning and transient clearing functions, / D. Armbruster, J. Fonteijn, M. Wienke // Logistics Research. - 2012. - VOL 87 - №3. - Р.815-822.
Armbruster D. Modeling production planning and transient clearing functions, / D. Armbruster, J. Fonteijn, M. Wienke // Logistics Research. - 2012. - VOL 87 - №3. - Р.815-822.
