Научная статья на тему 'О новом классе динамических моделей поточных линий производственных систем'

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

CC BY
299
78
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
PDE-МОДЕЛЬ / ПРОИЗВОДСТВЕННАЯ ЛИНИЯ / МАССОВОЕ ПРОИЗВОДСТВО / НЕЗАВЕРШЕННОЕ ПРОИЗВОДСТВО / СИСТЕМА УПРАВЛЕНИЯ / БАЛАНСОВЫЕ УРАВНЕНИЯ / ПРОИЗВОДСТВЕННОЙ ЛИНИИ / УРАВНЕНИЕ СОСТОЯНИЯ / ДИСКРЕТНО-СОБЫТИЙНАЯ / МОДЕЛЬ / ТЕОРИЯ МАССОВОГО ОБСЛУЖИВАНИЯ / МОДЕЛЬ ЖИДКОСТИ / CLEARINGФУНКЦИЯ / КВАЗИСТАТИЧЕСКИЙ ПРОЦЕСС / ПЕРЕХОДНЫЙ ПРОЦЕСС / СТОХАСТИЧЕСКИЙ ПРОЦЕСС

Аннотация научной статьи по математике, автор научной работы — Пигнастый О. М.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

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

УДК 658.51.012

О НОВОМ КЛАССЕ ДИНАМИЧЕСКИХ МОДЕЛЕЙ ПОТОЧНЫХ ЛИНИЙ ПРОИЗВОДСТВЕННЫХ СИСТЕМ

О.М. ПИГНАСТЫИ

Национальный Технический Университет "ХПИ"

e-mail: pom7@bk.ru

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

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

Введение

Мировые тенденции развития производственных систем в разных отраслях промышленности требуют повышения гибкости производственных линий. Функционирование промышленного предприятия в условиях неопределенности и неустойчивости внешней среды и потребность в регулярных изменениях номенклатуры выпускаемой продукции обуславливают необходимость построения эффективной системы управления предприятием, разработки новых математических моделей производственных линий [1-6]. Производство в условиях современного рынка достаточно быстро меняет ассортимент, качество и количество выпускаемой продукции [6, с. 12], функционирует в неустановившихся переходных режимах. Снижение себестоимости продукции достигается не традиционным увеличением выпуска продукции, как в случае установившегося производственного процесса, а в результате синхронизации материальных потоков вдоль технологического маршрута. Характерной чертой развития современного промышленного производства является то, что длительность цикла производства занимает существенную часть жизненного цикла изделия. При длительности производственного цикла в несколько месяцев для предприятий, например, с поточным методом организации производства полупроводниковой продукции жизненный цикл изделий не превышает одного года [5], [6, с. 51]. Это приводит к тому, что производственные линии значительную часть времени функционирует в переходном неустойчивом режиме. В связи с этим при проектировании систем управления производственными линиями в последнее десятилетие особо актуальной стала задача построения моделей управления для переходных неустановившихся режимов производства. Актуальность задачи настолько высока, что ведущие мировые предприятия (Intel [5]) создали экспериментальные лаборатории по исследованию переходных неустановившихся режимов.

Таблица 1

Основные параметры производственных линий

Наименование параметра Значение Источник

Кол-во техн. операций, единиц техн. оборудования M=150..1500 шт. [5, с. 4589]

Кол-во изделий в незавершенном производстве N=104.. 105шт. [7, с. 242]

Кол-во производимых номенклатур изделий K=1..10 шт. [8, с. 940]

Продолжительность жизненного цикла изделий 12 мес. [5, с. 4589]

Длительность производственного цикла Td=3..6 мес. [5, с. 4589]

Время принятия управленческого решения 0,5..1.0 час [9, с. 7]

Рис. 1. Структура производственной линии по изготовлению полупроводников [7]

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

Дальнейшее развития теории систем управления потоковыми параметрами производственных линий для неустановившихся переходных режимов привело к тому, что в последнее десятилетие были разработаны качественно новые типы моделей, сочетающие самосогласованные предметно-технологическое на микроуровне и потоковое на макроуровне описания производственной линии [1]. Опыт построения моделей физических систем, состоящих из множества объектов, подсказывает, что дальнейшее развитие теории моделей производственных систем может быть основано на статистическом подходе [11, 12], фундаментом которого является механизм стохастического воздействия оборудования на предмет труда и коллективного взаимодействия предметов труда между собой, и феноменологическом, позволяющем установить основные закономерности технологических явлений [2-5] на основе наблюдений и статистики без использования модельных представлений [1, 12] о характере взаимодействии предметов труда и оборудования. Общие закономерности большинства установившихся технологических процессов известны. Разным технологическим процессам соответствуют разные уравнения состояния [2-5]. Определенные с помощью обработки статистических данных уравнения состояния являются фундаментом феноменологического подхода построения квазистатических моделей производственных линий. Однако, представление переходного неустановившегося производственного процесса последовательной цепочкой квазистатических равновесных уравнений состояния [2-5] приводит к значительной погрешности. В общем виде уравнения состояний могут быть получены в рамках статистического похода, который позволяет однозначно связать макроскопические потоковые характеристики технологического процесса с микроскопическими предметно-технологическими параметрами достаточно большого количества предметов труда [1, 12].

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

Построение единой теории производственных линий

Наличие однообразных задач планирования и управления производством привело исследователей к идее создания единой теории производственных систем. Graves S.C. (1986) [13], Karmarkar U.S. (1989) [14] предложили использовать в качестве основных параметров состояния производственных линий объем незавершенного производства

(WIP: work in process) W, пропускную способность (capacity, throughput) \j\ и длительность производственного цикла (lead time) Td . Корректно определенное Karmarkar U.S. (1989) [14] уравнение состояния производственной линии = Ф(W), получившее название clearing-функции, задавало для устойчивых установившихся производственных процессов мгновенную связь между пропускной способностью , производственным

циклом Td и объемом незавершенного производства W .

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

\j\cl = °(W) = const (1)

ставит фиксированную границу выпуска продукции, предполагая мгновенное наращивание мощности производства. Constant Proportion Clearing-функция (Graves S.C., 1986, [13])

\%\CL = O(W) = a ■ W, a = const (2)

предполагает фиксированное время выхода производства на полную мощность, а при наличие ограничения выпуска продукции для промежутка планирования Д/г-, известна как комбинированная dearing-функция (combined, Karmarkar U.S., 1989) [14]:

[z\cl = ®<Ю = a ■W ^Omax , a = const. (3)

Важный класс составляют нелинейные dearing-функции, используемые для построения однопродуктовых моделей производственных систем, функциональная форма которых для простых случаев описания поточной линии производственной системы определена как:

аУЩ-модель М/М/1 очереди для стационарного состояния (рис. 2) [21] (2009)

\x\cl=Ow )=fW; (4)

б)модель фундаментальной диаграммы трафика для стационарного состояния

\z\cl =o(w) = v w - W2 , (5)

в)двухмоментная модель G/M/i очереди для стационарного состояния (Mehdi J., 1991) [22, с. 350]:

2 2 2 ТТЛ С2 + cs р X л

W = -!~ + Р, Р = -< 1, (6)

2 1-Р М

где ca2 и cs2 характеризуют среднеквадратичное отклонение поступления требований на

изготовления изделий и времени их обработки, м - темп обработки предметов труда, X -

интенсивность поступление предметов труда на первую технологическую операцию.

стесхтес- зватес-

--|1 |J. |J. — |1 Ц —' Ц (1 --Ц AREA NODE

Рис. 2. Схема М/М/1-модели очереди поточной линии производственной системы

Модель G/M/l очереди для стационарного состояния (6) является развитием модели (4) для описания движения предметов труда по технологическому маршруту с

последовательным расположением технологического оборудования. Устойчивое состояние для модели (1.10) обеспечивается при выполнении условия р = (//U)< 1 [23, c. 6]. При Х^и межоперационные заделы становятся бесконечно большими, W ^ да, а при /> и уравнение (6) не может быть использовано, т.к. при его выводе предполагалось р = (// и) < 1. Решение (6) относительно р = (//и) < 1 может быть представлено в виде:

V(W +1)2 + 4W(сг -1) - (W +1) c2a + cS Р =-^-' при С = — > 1' (6)

При с ^ 1 модель M/M/1 очереди (4) является предельным случаем модели (6) G/M/i:

w п U-W

lim р =-, lim= --(7)

с(W +1)' с^ТС 1 + W

Asmundsson J.M [21] предложил распределенные АС-функции (allocated clearing function) для моделирования многопродуктовых производственных систем. AC-функции предполагают модели, в который производится совокупный продукт, агрегирующий в себе используемые технологические ресурсы для производства отдельных продуктов. Альтернативный подход заключается в представлении dearing-функции производственной системы в виде суммы dearing-функций, описывающих выпуск отдельного продукта. Экспериментальные данные свидетельствуют об удовлетворительных результатах использования AC-функций для продуктов, схожих по характеру потребления ресурсов. Интенсивности потребления ресурсов при этом выражаются в единицах времени обработки. Однако, если модели переноса технологических ресурсов для каждой номенклатуры предметов труда являются сложными, использование AC-функций не позволяет удовлетворительно описать производственные процессы. Selguk B., Fransoo J.C., Gok A.G. (2007) аппроксимировали dearing-функцию кусочно-линейной функцией, что позволило использовать для оптимизации параметров производственной системы аппарат линейного программирования.

Clearing-функция может быть получена как аналитически так и численно c использованием стационарных или переходных TQ-моделей, DES-моделей, Fluid-моделей производственных систем или определена эмпирически. Selguk, B., Fransoo J.C.,Gok A.G. (2007) представили методику построения переходных dearing-функций аналитически. В связи с тем, что оперативная информация как о выпуске продукции, так и о состоянии незавершенного производства на предприятии является закрытой, то в большинстве исследований для построения dearing-функции вместо эмпирических данных используется TQ-,DES- и Fluid-модели. Исключением являются работы Haeussler S., Missbauer H. (2012), в которых для построения dearing-функции применены полученные с производственной линии цифровые носители. Kacar N. (2012) использовал для построении dearing-функции оптимизированные параметры производственной линии. Обширные вычислительные эксперименты, связанных с расчетом параметров поточных линий по производству полупроводниковых изделий компании Intel, показали хорошее приближение расчетных и экспериментальных данных для установившихся процессов [22].

Для вывода уравнение состояния использовано множество моделей, из которых следует выделить основные типы моделей и их комбинации: модели массового обслуживания (TQ-model) [15], дискретно-событийные модели (DES-model) [16], модели жидкости (fluid-model) [17]. Каждый тип имеет преимущества, но ни один не подходил в полной мере для моделирования как установившихся, так и переходных процессов работы производственной системы [2-4, 9].

Модели теории массового обслуживания (TQ-model) (рис. 3) используются для описания производственных линии в установившемся режиме [15]. Применение TQ-моделей для переходных режимов приводит к чрезмерному усложнению задачи. Существенным ограничением их применения является то, что современные производственные линии состоят из большого числа технологических операций. Как

правило, используются одномоментные модели теории массового обслуживания. Двухмоментные модели теории массового обслуживания, учитывающие дисперсию параметров обработки предметов труда, для описания производственных линий практически не встречаются в связи их сложностью построения. Следующим ограничением для применения TQ-model является то обстоятельство, что основные формулы получены для установившихся режимов (г ^ ж), что автоматически подразумевает выполнение неравенства Л</ (рис. 4.), где Л - интенсивность поступления предметов труда на обработку; / - интенсивность обработки предметов труда. В действительности, большинства производственных линий с конечным производственным циклом справедливо неравенство Л > /, что не позволяет использовать хорошо разработанный аппарат теории массового обслуживания для установившихся процессов. Решение указанных проблем приводит к высокой размерности задачи, громоздким вычислениям, требует значительных вычислительных ресурсов.

Njt\ W_(Y)[lots] в„

Рис. 3. Базовые модели теории массового обслуживания

Дискретно-событийные модели (DES-модели) удачно применены для описания работы поточных линий в переходном и стационарном режиме, содержащих малое количество технологических операций и обрабатывающих партии предметов труда небольших размеров. Требуемое процессорное время для расчета параметров производственной линии для партии деталей N=104.. 105шт. и 100 технологических операций при миллионе имитаций [8, с. 242] составляет несколько дней, что делает невозможным использования DES-модели при проектировании систем управления производственными линиями. Существенным ограничением использования DES-моделей является сложность построения аналитических зависимостей между входными и выходными данными производственной линии для указанного количества технологических операций и размера партии предметов труда. Значительные вычислительные ресурсы необходимы для реализации алгоритмов генерации случайных чисел с заданным законом распределения параметров обработки предмета труда. Особое внимание следует уделить вопросам устойчивости вычислительных схем. Объем испытаний, необходимый для получения оценок наблюдаемой величины с заданной степенью точности зависит от вида закона распределения случайной величины, корреляции элементов выборки, наличия и длительности переходного режима функционирования моделируемой системы. Недостатком DES-модели является медленная сходимость выборочных средних значений к истинным. На достоверность результатов влияет выбор генератора псевдослучайных чисел и зависимость результатов численного эксперимента от плана проведения эксперимента [20, с. 16]. Большое внимание при анализе результатов численных экспериментов уделяется оценке соответствия дискретно-событийной модели реальной производственной линии с использованием статистических критериев.

Рис. 4. Базовые дискретно-событийные модели

накопитель

накопитель

накопитель

Входящий

т-1 m m=N

S V.

оборудование

Исходящий

7 №1ч

оборудование

оборудование

Рис. 5. Концептуальное представление модели жидкости (Форрестер Дж., 1961 [18])

Дискретные и непрерывные жидкостные модели производственных линий (Fluid-модели) ориентированы на малое количество технологических операций и линейные стационарные решения [16, 17]. Основное применение нашли при изучении переходных режимов работы производственной линии. Повышение количества технологических операций приводит к увеличению размерности системы дифференциальных уравнений, что значительно усложняет вычисления. Основная сложность построения Fluid-моделей заключается в выводе уравнений состояния для параметров производственной линии [17]. В подавляющем большинстве случаев Fluid-модели представлены одномоментным уравнением переноса, дополненным квазистатическим уравнением состояния. В отличие от DES-моделей они менее пригодны для описания стохастических производственных процессов [7].

Предложенный Кармаркаром (Karmarkar U.S) механизм dearing-функций позволил упорядочить процесс описания производственных линий, выделив основные потоковые параметры, которые являются общими для разных типов моделей. Однако, dearing-функции унаследовали недостатки:

• малая эффективность при описании производственных линий, состоящих из большого количества разнообразного технологического оборудования (случай М>>1);

• сложность построения многомоментных моделей (2-х моментных);

• сложность построения многоресурсных и многопродуктовых моделей;

• сложность описания переходных режимов работы линии. Используются для описания квазистатичных производственных систем;

• сложность моделирования межоперационного состояния линии,

которыми обладали существующие TQ -, DES-, Fluid-модели, положенные в основу их построения.

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

Развитие теории производственных линий В последнее десятилетие для проектирования производственных линий разработаны модели, содержащие уравнения в частных производных (PDE-model). Введенный в [1] (2003) новый в класс моделей объединил преимущества TQ-моделей, DES-моделей и Fluid- моделей, что значительно расширило возможности проектирования систем управления производственными линиями. PDE-модели позволили учесть влияние внутренних факторов производства на пропускную способность поточной линии с учетом ограничений, которые определяются максимальной производительностью оборудования и емкостью накопителя [1-5]. Существенным преимуществом PDE-моделей является то, что они подробно позволили описать движение предметов труда от операции к операции, допускают решение в аналитическом виде и не требуют значительных вычислительных ресурсов. При построении PDE-моделей выделилось два подхода: статистический (Пигнастый О.М, Азаренков Н.А., Ходусов В.Д., Демуцкий В.П. [1, 12]) и феноменологический (Armbruster D., Berg R.A., Lefeber E., Ringhofer C., Rooda J., [2-5]). Основным уравнением PDE-модели при феноменологическом описании является уравнение сохранения количества предметов труда, находящих в процессе технологической обработки

шшотшшш

Балансовые уравнения

О

Рис. 6. Использование DES-,TQ- и Fluid-моделей в качестве уравнений состояния для замыкания балансовых уравнений PDE-моделей

ôp(t, x) ÔF (t, x)

dt

+ -

ôx

■ 0, F(t, x) = pt, x) ■ v(t, x) W(t) = jp(t, x)dx, x e[0,l]

(8)

где x - технологическая позиция предмета труда в технологическом маршруте; p(t, x) -плотность расположения предметов труда по позициям в тех. маршруте; F (t, x) - поток предметов труда вдоль тех. маршрута в зависимости от позиции; W(t) - объем незавершенного производства; v(t, x) - скорость движения предметов труда вдоль технологического маршрута. Уравнение (8), записанное в обозначениях, введенных в работах [3-5], является незамкнутым. Замыкание уравнения (8) многими авторами обеспечено с помощью уравнения состояния (dearing-функций). Ниже представлен анализ основных PDE-моделей производственных линий.

0

Для описания синхронизированных конвейерных линий довольно часто используется PDE-модель, содержащая Graves-уравнение состояния v(t, х) = с = const :

Х) + Х) = 0, Fit, х) = fit, х) ■ vit, х) vit, х) = с = const х е [0,l] (9)

Модель (9) позволяет получить основные потоковые параметры конвейера, движение предметов труда на котором происходит с постоянной скоростью v(t, х).

Для описания производственных линий с ограниченной пропускной способность удачно зарекомендовала себя нелинейная PDE-модель Lighthill-Whitham:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

+д1Ш = 0, Fift, х)) = ft, х) ■ vft, х)), vit, х) = vo(1 , vo = const (10)

dt дх R

в которой пропускная способность ограничена максимальной плотностью межоперационных заделов R. При достижении плотности предметов труда в межоперационном заделе критического значения движение предметов труда по технологическому маршруту останавливается. LW PDE-модель (10) нашла широкое применение в задачах передачи трафика и при описании движущего потока автомобилей .

Заслуживает внимания квазистатическая PDE-модель, использующая нелинейное Кагтагкаг-уравнение состояния

+Цтх1 = 0, F it, х) = ft, х) ■ х), vit) ==-/-, (11)

dt дх M + W

где M-количество единиц технологического оборудования вдоль технологического маршрута, / -интенсивность обработки предмета труда, W-объем незавершенного производства. Скорость движения предметов труда vit) в отличие от модели (9) является зависимой от времени. Модель используется для проектирования конвейерных линий с переменной скоростью движения. При увеличении количества предметов труда W в незавершенном производстве скорость конвейера снижается. PDE-модель (11) предложена Д.Армбрустером (Armbruster), получила дальнейшее развитие в работах Э.Лефебера (Lefeber):

+ = 0, F it, х) = Ht, х) ■ vit, х) vit, х) =-л--(12)

dt дх M + pit, х)

Модель (12) учитывает в уравнении состояния для установившегося процесса зависимость распределения предметов труда по технологическим позициям. Увеличение количества предметов труда в межоперационном заделе приводит к снижению скорости обработки изделий. Шагом в развитии двухмоментных моделей стала предложенная Э.Лефебером (Lefeber) двухмоментная PDE-модель, в которой уравнение состояния заменено уравнением Бюргерса.

dpit, х) OF {t, х) п ^ N N / N dvit, х) 1 dvit, х)2 n , л

У +—= 0, Fit,х) = ft,х) ■ vit,х), У + T \ = 0 (13)

t х t 2 х

Скорость движения предметов труда в двухмоментной PDE-модели (13) так же, как и в одномоментной PDE-модели (12) зависит технологической позиции обработки. Однако в отличие от (12) скорость движения предметов труда в модели (13) определяется не из уравнения состояния для установившегося процесса, а путем решения балансового уравнением Бюргерса. Для описания поведения производственных линий при наличие экспериментальных данных Д.Армбрустером (Armbruster) предложена диффузионная PDE-модель [5]

ФМ + = 0^ F it, х) = Cit ) ■ fit, х) - Dit ) дР^, (14)

t х х

в которой для расчета коэффициентов Cit) и Dit) использованы статистические данные,

определяющие технологические траектории предметов труда. Точность вычислительных результатов, полученных с использованием диффузионной PDE-модели определяется точностью имеющихся статистических данных. Одномоментная модель (12) успешно применяется для описания производственных линий с накопленной статистикой данных о

процессе обработки предметов труда. Однако использование диффузионной PDE-модели для проектирования новых производственных линий не представляется возможным. Перспективы развития диффузионной модели тесно связаны с теорией подобия производственных систем [12]. Описанные выше модели используют при замыкании уравнение состояния для установившихся режимов работы производственной линии. Построение уравнений состояния для переходных неустановившихся режимов связано со значительными трудностями [2-5].

Использование автором статистического подхода позволило выразить параметры производственной линии выразить через состояния большого количества предметов труда и построить замкнутую многомоментную PDE-модель для переходных неустановившихся режимов (2003) [1]

ММ) +ММ) ^к > 0, = (15)

а ôS о о

дЫк(t, S) +1(t,S) \х\у (t,S) ô {\x\v(t,S)^ œ

" +--+ k "

У

Ык-i(t, S ) = \nkG(t, S, j)dj

0

ôt ôS \ы]о (t, S) ôS ^ \ы]о (t, S)

где S - технологическая позиция предмета труда в технологическом маршруте, S е \о, Sd ] ; \ы]о (t, S ) = p(t, S ) - плотность расположения предметов труда по позициям в тех. маршруте; \j]i(t, S) = F (t, S) -поток предметов труда вдоль тех. маршрута в зависимости от позиции; %(t,S,p) функции распределения предметов труда по состояниям; \%]\у (t, S) - темп обработки предметов труда вдоль тех. маршрута в зависимости от позиции S ; G(t, S, ju) -

функция, определяющая процесс переноса оборудованием технологических ресурсов на предмет труда [1]. Для производственной линии с количеством предметов труда в незавершенном производстве N>>1 и большим количеством технологических операций (таблица№1) система уравнений (15) в двухмоментном приближении принимает вид [12]

«м+Мм=о, Мм(,.s)-Ы1((,. (16)

ôt as ôt Ыо«,S) ôs v-]1f(' ' ы11 • 4ôS ыо«,s) ( )

Замкнутая двухмоментная система уравнений (16) определяет поведение параметров производственной линии для переходных неустановившихся режимов. Для синхронизированной конвейерной производственной линии, работающей в установившемся режиме, система уравнений (16) упрощается [1]:

â^M + ôWM = о ; + v(t,S) = о ; . v(,,S) (17)

ôt ôs ôt ôs Ы (t, S )

принимает вид (13), а для одномоментного приближения (12):

^ЪМ) +ММ = о ; \x\(t, S ) = \X1r(!.S) (18)

Для синхронизированной линии j = v(t, S) • \ы]о (t, S), откуда скорость движения

предметов труда по технологическому маршруту выразим через интенсивность обработки предметов труда ju

v(t, x) = , (19)

P(t, x)

откуда получим систему уравнений (19) в форме (12). Если скорость движения предметов труда по технологическому маршруту представить усредненной величиной то система уравнений (19) может записана в форме (11).

J v(t, S ) S )dS

V(t) = —S-= , S e[0, Sd ] (20)

Sd W(t)

J\z\o(t, S )dS

0

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

Выводы

При построении нового типа моделей производственных линий использованы два подхода - феноменологический подход (Armbruster D., Berg R.A., Lefeber E., Ringhofer C., Rooda J.) и статистический (Пигнастый О.М., Азаренков Н.А, Ходусов В.Д., Демуцкий В.П.). Феноменологический подход дал возможность построить ряд моделей производственных линий, дополнив уравнения переноса, уравнением состояния в форме clearing-функции. Это позволило построить PDE-модели производственных линий для наиболее простых случаев функционирования. Обоснованность применения определялась сравнительным анализом результатов, полученных с помощью DES-модели и исследуемой PDE-моделью. Однако, хотя построенные PDE-модели с помощью феноменологического подхода и являются прогрессивными, дальнейшее их развитие ограничено. Ограничение обусловлено тем, что уравнения переноса и дополняющие их уравнения состояния для более сложных случаем не могут быть выведены из феноменологического представления. Для построения таких моделей необходимо использовать статистический подход. Показательно то, что общий характер статистических закономерностей не зависит от того, каким способом описывается поведение отдельного предмета труда. При весьма большом числе предметов труда появляются новые своеобразные статистические закономерности. Статистические закономерности обусловлены именно наличием большого числа предметов труда, ни в коем случае нельзя свести к закономерностям поведения отдельного элемента. Использование статистического подхода позволяет получить замкнутые многомоментные балансовые уравнения (уравнения переноса) не из феноменологических соображений, а исходя из законов движения отдельных предметов труда по технологическому маршруту, определенных технологией производства.

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

Список литературы

1. Демуцкий В.П. Теория предприятия: Устойчивость функционирования массового производства и продвижения продукции на рынок / Демуцкий В. П., Пигнастая В. С., Пигнастый О.М. Х.: ХНУ, 2003. - C. 272.

2. Лысенко Д.Э. Автоматизированный синтез моделей технологических процессов на основе прецедентного подхода / Д.Э. Лысенко // Научные ведомости БелГУ. Сер. История. Политология. Экономика. Информатика. - 2013. - № 1 (172). - Вып. 29/1. - С. 121-129

3. Berg R.A. Modelling and Control of a Manufacturing Flow Line using Partial Differential Equations. IEEE Transaction son Control Systems Technology / Berg R.A., Lefeber E., Rooda J.E. -Boston, 2008. - P. 130-136.

4. Lefeber E. Modeling, Validation and Control of Manufacturing Systems / E.Lefeber, R.A.Berg, J.E. Rooda - (Proceeding of the 2004 American Control Conference). - Massachusetts, 2004. -P. 4583-4588.

5. Armbruster D. Continuous models for production flows. In Proceedings of the 2004 American Control Conference / Armbruster D., Ringhofer C., Jo T-J. - Boston, MA, USA, 2004. - P. 4589-4594.

6. Лысенко Ю.Г. Моделирование технологической гибкости производственно-экономических систем / Ю.Г.Лысенко, Н.В. Руменцев. - Донецк: ДонДУ, 2007. - C. 238.

7. Bramson M. Stability of queueing networks, lecture notes in mathematics, Journal of Probability Surveys / M. Bramson. - Netherlands, 2008. - Vol. 5. - P. 169-345.

8. Schmitz J.P. Chaos in Discrete Production Systems / Schmitz J.P., Beek D.A., Rooda J.E. -Journal of Manufacturing Systems, 2002. - Vol. 21, no. 3. - P. 236-246.

9. Kempf К. A Continuum Model for a Re-entrant Factory. Operations research / K. Kempf, D.Marthaler, C.Ringhofer, D.Armbruster, J.Tae-Chang. - 2006. - Vol 54. - № 5. - Р. 933-950 .

10. Vollmann T.E. Manufacturing Planning and Control for Supply Chain Management / T.E.Vollmann, L.Berry, D.C.Whybark, F.R.Jacobs. McGraw-Hill, New York, 2005. - P. 520.

11. Красовский А.А. Фазовое пространство и статистическая теория динамических систем / А. А. Красовский. - М.: Наука, 1974. - C. 232.

12. Азаренков Н.А. К вопросу подобия технологических процессов производственно-технических систем / Н.А.Азаренков, О.М.Пигнастый, В.Д.Ходусов - Доповвд Нащонально! академи наук Украши, 2011. - N02. - С. 29-35, http://dspace.nbuv.gov.ua/handle/123456789/37227/

13. Graves S.C. A tactical planning model for a job shop. Operations Research 34 (4). /S.C. Graves. New York, 1986. - P. 522-533.

14. 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.

15. Gross D. Fundamentals of Queueing Theory / D.Gross, C.M.Harris. - New York, 1974. -

P. 490.

16. Ramadge P.J.The control of discrete event systems. Proceedings of IEEE / P.J. Ramadge, W.M. Wonham. - 1989. - Vol. 77(1). - P. 81-98.

17. Perturbation analysis for on-line control and optimization of stochastic fluid models / C.Cassandras ,Y.Wardi ,B.Melamed ,G.Sun, C.Panayiotou IEEE Trans.Autom.Control. - Netherlands, 2002. - Vol. 47, - № 8. - P. 1234-1248.

18. Форрестер Дж. Основы кибернетики предприятия / Дж. Форрестер. - М.: Прогресс, 1961.

C. 341.

19. Коробецкий Ю.П. Имитационные модели в гибких системах / Ю.П. Коробецкий, С.К. Рамазанов. - Луганск: Изд. ВНУ, 2003. - C. 280.

20. Asmundsson J.M. Production planning models with resources subject to congestion. / Asmundsson J. M., Rardin R. L., [et al]. Naval Res Logist 56, 2009. - P. 142-179.

21. Mehdi J. Stochastic Models in Queuing Theory / J. Mehdi. - New York, 1991. - P. 482.

22. Berg R. Partial differential equations in modelling and control of manufacturing systems / R. Berg. Netherlands, Eindhoven Univ. Technol., 2004. - P. 157.

A NEW CLASS OF DYNAMIC MODELS FLOW LINES OF PRODUCTION SYSTEMS

O.M. PIGNASTY

National Technical University "KhPI"

e-mail: pom7@bk.ru

This paper discusses a class of models of production systems [1]. Models are widely used today for building effective systems of production lines [2, 3]. Class model parameters determine the behavior of the production line with the help of partial differential equations (PDE-model) [2-4], in the last decade successfully used to describe the quasi-static and transient transients The article analyses the shortcomings, appeared in Russia on the whole and Altai Region in particular, of the municipal self-government model and defines some perspectives of the municipal system development^].

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.

i Надоели баннеры? Вы всегда можете отключить рекламу.