Математическая модель единой системы конструкторско-технологической подготовки и производства на промышленном предприятии
В.В. Соседко, А.Г. Янишевская
Для успешного выполнения на промышленном предприятии производственной деятельности, необходимо наличие единой системы конструкторско-технологической подготовки и производства (КТПП). Она охватывает процессы проектирования, конструкторской подготовки, технологической проработки, подготовки документации (конструкторской, технологической, программной), необходимой для производства и непосредственно производства.
Математическая модель единой системы КТПП строится на базе постулатов теории вероятности, в частности, теории массового обслуживания и математической статистики.
В процессе создания моделируются потоки информации, отражающие состояния производства на предприятии, по схеме марковских случайных процессов.
Многие разновидности производственной деятельности на предприятии, которые приходится анализировать при выборе оптимального решения, развиваются как случайные процессы.
Для того чтобы вычислить числовые параметры, характеризующие эффективность в принятии решений при конструировании, технологической подготовке производства приборов и управлении производством, нужно построить вероятностную модель процессов, учитывающую сопровождающие эти процессы случайные факторы.
Для математического описания КТПП, а также для автоматизации управления производством, как правило, развивающихся в форме случайных процессов, может быть применен математический аппарат, разработанный в теории вероятности для так называемых марковских случайных процессов [1, 2].
Определим марковский случайный процесс. Пусть имеется некоторая физическая система S на предприятии, состояние которой меняется с течением времени. Если состояние системы S меняется во времени случайным, заранее непредсказуемым образом, то в системе S протекает случайный процесс [3].
Случайный процесс, протекающий в системе S, можно будет отнести к марковскому процессу, когда он будет протекать без последействия (т.е. в нем будущее состояние не зависит от прошлого, и соответственно прошлое состояние не оказывает влияние на будущее).
Если рассматриваемый процесс на предприятии обладает такими свойствами, когда для каждого момента времени t0 вероятность любого состояния процесса (системы) в
будущем (при t > t0) зависит только от его (ее) состояния в настоящем (при t = t0) и не
зависит от того, когда и каким образом процесс пришел в это состояние, то он будет также марковским.
При анализе случайных процессов с дискретными состояниями, протекающих при КТПП и в управлении производством, удобно пользоваться геометрической схемой - так называемым графом состояний. Графом состояний можно описать потоки информации в управлении сложным производством, случайную последовательность событий в котором следует назвать марковской цепью. Марковская цепь - последовательность случайных событий с конечным или известным числом исходов, характеризующаяся тем свойством, что, говоря нестрого, при фиксированном настоящем будущее независимо от прошлого. Размеченный граф состояний дает возможность определить вероятности состояний
Pl(t X P2O X ■■■, Pn (t)
как функции времени. Эти вероятности удовлетворяют определенного вида дифференциальным уравнениям, так называемым уравнениям Колмогорова - особого вида дифференциальным уравнениям, в которых неизвестными функциями являются вероятности
состояний. Решение этих уравнений дает возможность определить вероятности состояний производства как функции времени.
Поток первичной информации, отражающий ход КТПП и управления производством на предприятии, может поступать в виде документов и сообщений в объединенный комплекс (далее система), в состав которого входят система управления данными об изделии (PDM-система) - организационно-техническая система, обеспечивающая управление всей информацией об изделии и автоматизированная система управления предприятием (АСУП).
Поток информации о ходе разработки изделия или о состоянии производства, как правило, будет регулярным, так как пакет данных или сообщений, получаемый системой, следуют один за другим через строго определенное время. Именно такой поток информации рекомендуется организовывать на предприятии независимо от того, будет ли он образовываться из первичных документов (например, чертежей) или при автоматизированном потоке информации.
Тогда такой поток информации будет стационарным, без последствия и ординарным.
Стационарность потока информации о состоянии производственного процесса на предприятии будет определяться его однородностью по времени.
Отсутствие последействия потока информации в процессах проектирования и управления на предприятии определяется тем, что первичные документы или сообщения появляются в последовательные моменты времени независимо друг от друга.
Ординарность потока информации заключается в том, что первичные документы или сообщения обрабатываются на предприятии поодиночке, а не организованными парами, тройками и т.д.
В идеале поток информации выглядит так, он стационарный, без последствия и ординарный.
Если организовать поток информации на предприятии так, как было предложено выше, то он будет обладать всеми тремя свойствами, и, следовательно, такой поток будет простейшим или стационарным пуассоновским потоком. Поток событий называется простейшим (или стационарным пуассоновским), если он обладает сразу тремя свойствами: стационарен, ординарен и не имеет последействия. Название «простейший» связано с тем, что процессы, связанные с простейшими потоками, имеют наиболее простое математическое описание.
Тогда для этого простейшего потока среднее число событий в единицу времени -поток интенсивности X будет величиной постоянной: Л = const.
Однако если организовать поток сообщений о состоянии производства, который будет передаваться в реальном масштабе времени, то он будет нестационарным пуассоновским потоком. В этом случае поток событий не имеет последействия, ординарен, но нестационарен. В таком потоке интенсивность X зависит от времени
Л = Л($).
Пуассоновский поток событий (как стационарный, так и нестационарный) тесно связан с известным распределением Пуассона. Распределение Пуассона моделирует случайную величину, представляющую собой число событий, произошедших за фиксированное время, при условии, что данные события происходят с некоторой фиксированной средней интенсивностью и независимо друг от друга. Распределение Пуассона играет ключевую роль в теории массового обслуживания.
Число событий потока (от оформления первичных документов до сдачи деталей и узлов в производстве), попадающих на любой участок времени в смене, распределено по закону Пуассона. Вероятность попадания на этот участок времени равна m событий и выражается формулой (1):
а
т
т
т!
е а при т = 0,1,...,
(1)
где а - среднее число событий, приходящееся на конкретный участок времени т. Величина а выражается формулами:
-для стационарного пуассоновского потока а = Я -т, то есть не зависит от того, где на оси 0t взят участок т, но нас интересует нестационарный пуассоновский поток;
того, в какой точке начинается участок т (например, начало обработки документации
цехами и службами течение прошедших суток после 24 ч).
При исследовании состояний КТПП приходится неизбежно сталкиваться с анализом своеобразных систем на предприятии, которые могут быть описаны системой массового обслуживания (СМО). Каждая СМО состоит из какого-то числа обслуживающих единиц, которые предлагается называть "каналами" обслуживания. В данном случае в качестве "каналов" могут фигурировать: потоки информации между отделами и цехами, потоки первичных документов о разработке продукции, потоки информации из цехов в систему управления производством об изготовлении заготовок и деталей и, наконец, потоки сообщений в PDM-систему о процессе обработки данных, поступивших из цехов и служб предприятия.
Законы образования очереди на обслуживание при использовании теории массового обслуживания бывают различными. Они определяются процедурой обслуживания очереди прохождения заявок.
Первый вариант - поступившие на обслуживание сообщения в любой сети становятся в конец очереди и обрабатываются в порядке прихода. Такой вариант обслуживания очереди называется естественным порядком, то есть первое поступившее сообщение обслуживается первым.
Второй вариант - это когда сообщение, ожидающее обработки (очереди), может обслуживаться не в порядке его поступления, а в зависимости от значимости поступившей информации на обслуживание, то есть требуется определить порядок прохождения заявок через приборы обслуживания. Таким образом, получается структура СМО со многими приборами обслуживания.
В реальных условиях на предприятии обслуживание ведется по второму варианту.
При автоматизированном управлении производством система массового обслуживания потоков информации в виде первичных документов или сообщений может быть организована тремя способами: с отказами, с ожиданием и смешанным.
При обработке информации система с отказом может образовываться, если документы или сообщения поступили после 24 часов. Они автоматически переходят в поток информации следующих суток и войдут как информация, отражающая состояние производства, в следующие сутки. Поэтому предлагается в PDM-системе потоки поступающей документации организовывать как поток заявок в системе с ожиданием. В системе с ожиданием заявка не покидает канал, а становится в очередь и ждет, пока не освободится очередной пункт обслуживания.
В качестве универсального метода организации потока информации о состоянии производства на приборостроительном предприятии предлагается поток первичных документов и сообщений организовывать с ограниченным последействием.
Если принять Т как интервал времени между соседними событиями в этом потоке, то плотность распределения / (^) случайной величины Т можно выразить зависимостью:
для нестационарного пуассоновского потока
и, значит, зависит от
/(?) = Яе Я при (? > 0)
Закон распределения с плотностью (1 ) будет показательным или
экспоненциальным.
Найдем числовые характеристики случайной величины Т:
-математическое ожидание (среднее значение) Mt = -1;
Я
-дисперсию Dt = ;
Я
Таким образом, интервалы времени между поступлениями данных о ходе производства
T T T
±1, 12, Li, ...
представляют собой независимые, одинаково распределенные случайные величины, то есть образуют поток Пальма. Потоком Пальма называется поток, обладающий свойствами стационарности, ординарности и независимости интервалов времени Т между событиями.
Как было сказано ранее, на предприятии поток данных, отражающий состояние производства, можно строить в виде стационарного пуассоновского. Простейший поток есть частный случай потока Пальма.
В нем интервалы времени Тх, 72, ., Ti,... представляют собой случайные величины,
распределенные по одному и тому же показательному закону, имеющему экспоненциальное распределение (следовательно, мы имеем дело с потоком Пальма). Их независимость следует из того, что простейший поток без последствия, и расстояние по времени между любыми двумя событиями не зависит от расстояния между другими.
Однако математическое описание потока первичной информации на приборостроительном предприятии по принципу потока Пальма будет очень громоздким.
При декомпозиции пунктов сбора информации на приборостроительном предприятии в разрезе отделов и цехов поток документов или сообщений следует рассматривать в каждой отдельной службе "прореженным", как поток Эрланга.
Поток, образованный каждым п - м событием простейшего потока, называется потоком Эрланга порядка п. Обобщенный закон распределения Эрланга (п -1)- го порядка ((п -1) - го цеха) в данном случае будет иметь плотность в виде выражения:
£„-і(Х) =
~Я;Х
Є
(-1)”-1,ПЯХ-п-----------при х > 0
-■ '=■ п (Я-я )
к ^1
0 при х<0
где к = 1,2, ., у -1, у +1, ., п, кроме Я - Я'
Таким образом, использование закона распределения Эрланга позволит найти усредненную величину потока информации каждой службы, и затем определить плотность потока в целом по предприятию, а также оценить переход производства от состояния к состоянию цехов и производства в целом.
Поток информации, деталей и узлов на заводе в соответствии с маршрутной технологией отображается первичными документами (конструкторской документацией, заявками, накладными, счетами-фактурами, актами отгрузки готовой продукции и т.п.), что переводит производство из состояния в состояние в реальном масштабе времени. Документы фиксируют изменение состояния в производстве с дискретностью 24 часа (сутки). Сдача готовых заготовок, деталей и сборочных единиц - товарной продукции цеха, как правило, осуществляется по маршрутным листам. На автоматизированных рабочих местах приема и сдачи деталей организован запрос информации с формированием маршрутного листа со штрих-кодом, содержащим всю основную
информацию о детали (наименование, материал, партия, сроки) и технологическим маршрутом. Данный документ выводится на принтер, а также может быть прочитан сканером штрих-кодов. Номенклатура реквизитов и объем записей в документе при автоматизированном управлении одинаковы для всех видов продукции, то есть документ унифицирован под идеологию учета. Поэтому маршрутный лист на товарную продукцию цеха можно принять как единицу плотности потока событий X, переводящего службы и предприятие из состояния х. в состояние х . .
В конце суток уведомления из каждого подразделения цеха завода поступают в систему для их обработки и решения задач планирования, учета и регулирования производства. На основании решения комплекса задач оперативного управления производством и выдачи машинограмм службам цехов и завода (к началу нового рабочего дня) они переходят в очередное состояние.
Структуру потоков документов от цеха к цеху, от подразделений в систему и обратно предлагается построить в строгом соответствии с маршрутной технологией (и маршрутным листом в частности), расписанной для каждой детали и сборочной единицы.
Размеченный граф изменения состояний КТПП и производства в подразделениях, и в целом на предприятии, с дискретностью 24 часа (сутки) представлен на рис. 1. Граф состояний охватывает все иерархические уровни структуры завода от инженерных служб и цехов разных профилей до склада готовой продукции.
В связи с особенностями маршрутных технологий при организации производства по изготовлению и выпуску сложных машин на приборостроительном предприятии, а также для более эффективного управления было принято решение подразделения завода группировать по принципу родственных специальностей в одной группе передела.
Предлагается следующая классификация групп, сложившаяся на предприятиях (см. на рис. 1 в направлении сверху вниз):
\т31_4,Я
11и3142Я
31^ [^42,31^ Хр^сгплр^
Рис. 1. - Размеченный граф состояний КТПП и производства с описанием переходов
^4 рст ^
1. Конструкторские отделы. Выполняют конструкторскую подготовку изделий. Разрабатывают сами или осуществляют сопровождение производства изделий. В зависимости от специфики разрабатываемых изделий их количество имеет конечное значение.
2. Технологический отдел. Подразделение представлено в данной группе в единственном числе. Осуществляет технологическую проработку изделий разрабатываемых и подготовку производства разработанных.
3. Цеха заготовительного производства. Данные подразделения изготавливают
детали, используемые в дальнейшей сборке основных изделий. Профили заготовительных цехов: гальвано-малярный, механический, деревообрабатывающий, пластмассовый,
инструментальный и другого направления. В зависимости от выполняемой работы цеха передают изготовленную продукцию либо в сборочные цеха, либо непосредственно на склад готовой продукции, либо при работе по кооперации, в другой цех в соответствии с технологическим маршрутом деталей.
4. Цеха сборочные. Данные цеха обеспечивают сборку, запуск, регулировку, испытания и сдачу блоков и изделий. В зависимости от специфики предприятия могут быть выделены цеха по сборке гражданской продукции и изделий спецтехники либо в соответствии с загрузкой и профилем цеха по сборке крупных изделий, преимущественно механики и цеха по сборке сложных радиоэлектронных изделий.
5. Склад готовой продукции. Здесь решаются задачи сбыта готовой продукции потребителю. Поток документов на данном уровне определяет этап производства оформления товара предприятия - очень важная составляющая в его жизнедеятельности, в то время как уровни 1-4 формируют товарную продукцию цехов.
Первичные информация и документы, сопровождающие поток готовых деталей и сборочных единиц по заводу, образуют потоки информации, переводящие состояние проектирования, технологической подготовки и производства изделий с одного уровня на другой (из одного статуса в другой).
Если представить среднее число событий, описанных в одном типовом первичном документе или партии деталей или узлов (с соответствующим маршрутным листом, тоже являющимся документом), через X, то полный поток документов, поступающих для обработки в систему, будет выражаться зависимостью х • X, где х принимает значения к, t, m в графе состояний (см. рисунок 1). В нашем случае X = const, то есть мы имеем дело с простейшим потоком.
При этом:
к - первичные документы на запуск процессов конструкторской подготовки и технологической проработки - бизнес-процессы, запускаемые в PDM-системе;
t - первичные документы на запуск процессов технологической подготовки -бизнес-процессы, запускаемые в PDM-системе;
m - первичные документы, сопровождающие детали, узлы, блоки и изделия на производстве - маршрутные листы;
Между подразделениями в процессе КТПП и непосредственно производства маршрутной технологией предусмотрены согласование документов и кооперативные услуги по доработке деталей. Документы возвращаются на доработку или на следующий этап работы с ними, а детали после доработки в другом цехе возвращаются в свой цех. Их сопровождают маршрутные листы и уведомления в системе, образуя плотность “потока обслуживания” л, то есть поток, непрерывно обслуживающий заявки других цехов. Суммарный поток таких документов выражается зависимостью х • /л .
Анализ потоков первичных документов, сопровождающих потоки деталей и узлов на приборостроительном предприятии, или сообщений, полученных после обработки документов, показал, что полученный размеченный граф отображает многоканальную СМО.
Для каждого из 14 состояний системы S графа определим предельные вероятности
состояний и составим уравнения Колмогорова:
^11 : _ —к\\,ТО^Ри + ктО,и№тО
ш
*12 : ТГ = -^'12,ТО^Р12 + ктО,12№РтО
ш
*13 : ТГ = —^13,ТО^Р13 + ктО,13№РтО
м
*14 : _ -^14,тО^Р14 + ктО,14№РтО
м
*то : ~7 _ ^11,то^р11 — ^то,иМрто + ^'12,тоДр12 — кто,12№рто + к1зтоЛр13 — кто,1з№то +
ш
+ ^14то^р14 — кто ,14№то — ^то,з1Лрто + ^з1;то№з1 — ^то,з2Лрто + ^з2,то№з2 — ^то,ззЛрто +
+ !зз,то№>зз — Iто ,34^ртО + ^З4,то№з4 ^тО,3а^ртО + ^3а,тоМРза
Шр31
*31 : _ —т31,32^Р31 + т32,31№Р32 + ^тО,З1^РтО — ^З1,тоМРЗ1 — т31,41^Р31 + т41,31№Р41 —
ш
— т31,42ЛРз1 + т42,31^42 — - т31,4рЛРз1 +
т4р,31^р4р т31,ст^р31 + тСГП ,31№РСГП
Мр32
*32 : _ -т32,31№Р32 + т31,32^Р32 + ^тО,З2^РтО — ^З2,то№з2 — т32,41^Р32 + т 41,32№Р 41 —
М
— т32,42Лр32 + т42,32№Р42 т32,4рЛр32 + т4р,32Ур>4р т32,СГПЛр32 + тСГП ,32№РСГП
Мр33
*33 : _ ^тО,ЗЗ^РтО — ^33,тМрз3 — т33,41^Р33 + т41,32№Р41 — тз3,42^Р33 + т42,33№>42 "'
М
— т33,4рЛРз3 + т4р,ззМр4р тзз,СГПЛрзз + тсгп ,ззМрсгп
Шр34
*34 : _ ^тО,З4^РтО — ^З4,то№з4 — т34,41 ДРз4 + т 41,34№> 41 — т34,42 ДРз4 + т42,34№Р42 " '
ш
т34,4рЛр>34 + т4р,34№)4р — т34,СГПЛр>34 + тСГП ,34№>СГП
шр
*3а : ~7~ = ^то,3аЛртО — ^3а,то№за — т3а,41Лр)3а + т41,3а№>41 — т3а,42Лр>3а + т42,3а№>42 "
т
— тз а ,4р^рз а + т4р,3 а №>4р — тз а ,СГПЛр>3 а + тсГП ,3 а №сГП
*41 : _ т31,41^р31 — т41,31^р>41 + т32,41^р32 — т41,32^р>41 + т33,41^р33 — т41,3зМр41 +
ш
+ т34,41^р34 — т41,34^р41... + т3а,41^р3а ~ т41,3аМр41 — т41,СГП^р41 + тСГП,41^р>СГП *42 : _ т31,42^р31 — т42,31^р>42 + т32,42^р32 — т42,32^р>42 + тз3,42^р33 — т42,33^р>42 +
ш
+ тз4,42^рз4 — т42,34^р42... + т3 а ,42^р3 а — т42,3 а №42 — т42,СГП^р42 + тСГП ,42№)СГП
_ шр4р ~ ~ ~
*4в : ^ ~ т31,4р^Р31 — т4р,31№р>4р + т32,4р^р32 — т4р,32№р>4р + т33,4р^р33 — т4р,33№р>4р +
+ тз4,4РЛРз4 т4р,34№Р4р ’" + тз а ,4рЛрз а т4р,3 а №Р4р т4р ,СГП
Яр4в +
тСГП ,4Р№РСГП
*СГП : С _ т31,СГПЛр31 — тСГП,31№РСГП + т32,СГПЛр32 — тСГП,32№РСГП + т33,СГПЛр33 —
ш
— тСГП ,33№РСГП + т34,СГПЛр34 — тСГП ,34№РСГП + т3 а ,СГП^Р3 а — тСГП ,3 а №РСГП +
+ т41,СГП
Лр41 тСГП,41 МрСГП + т42,СГП Хр42 тСГП ,42№РСГП + т4р,СГП^Р4р тСГП ,4в№РСГП
Составим для каждого передела (конструкторская и технологическая подготовка, заготовительные и сборочные цеха) системы уравнений Колмогорова.
Состояние конструкторской подготовки:
Фп
шр
12
йр1з
шР14
— —к11,тоЛр11 + кто,11№рто
— к12,тОЛр12 + кто,12№>то
— —к1зтоЛр13 + кто,13№рто
— -к14,тОЛр14 + ктО,14№РтО
Состояние технологической подготовки:
шр
то
(кцто р11 + к12то р12 + к13,то р13 + к14,то р14 (^то,31 + ^то,32 + ^тО,33 + ^то,34
+ ^тО,3а )РтО + (^31,тОР31 + ^32,тОР32 + ^33,тОР33 + ^34,тОР34 + ^3а,тОР3а (ктО,11 +
+ ктО,12 + ктО,13 + ктО,14 )РтО
Состояние заготовительного производства:
шр3
(то,31 рто (т31,32 + т31,41 + т
131,42 + т31,4в
ш
— ^31,тор31 + т41,31Р41 + т42,31Р42 + т4в,31 Р4в + тСГП,31РСГП )М
ёр3
+ тз1,СГП ) р31)Л + (т32,31 р32
(т31,32Р31 + ^тО,З2РтО (т32,41 + т32,42 + т32,4в + т32,СГП )Р32^Л + (т41,32Р41 +
+ т42,32Р42 " ' + т4в,32Р4в + тСГП,32РСГП (т32,31 + ^32,тОР32)№
шр3
(то,33 рто (т33,41 + т33,42 " ' + т33,4в + т33,СГП ) рзз)Л + (т 41,33Р 41 + т42,33Р 42
+ т4в,зз р4в + тСГП ,33 рСГП — ^зз,то рзз)М
шрз.
(?то,34рто (т34,41 + т34,42 + т34,4р + т34,СГП )Р34)Л + (т41,34Р41 + т42,34Р42
ш
+ т4в,34 р4в + тСГП ,34 РСГП — ^34,тО Р34)М
Щ,р
= (^то,3арто — (т3а,41 + т3а,42 + т3а,4в + т3а,СГП ) р3а)Л + (т41,3аР41 + т42,3аР42 '"
+ т4р,3ар4р + тСГП ,3аРСГП — ^3а,тО Р3а)№
Состояние сборочного производства:
ф
(т31,41 Р31 + т32,41Р32 + т33,41 Р33 + т34,41 Р34 + т3а,41Р3а т41,СГП Р41~)Л +
+ (тСГп ,41РСГП — (т41,31 + т 41,32 + т41,33 + т
141,34 + т41,3а>У 41
,, — (т31,42 Р31 + т32,42 р32 + т33,42 Р33 + т34,42 Р34 + т3а,42 Р3а — т42,СГП р42)Л +
ш
+ (тСГп ,42 РСГП (т 42,31 + т42,32 + т42,33 + т■
) Р42)М
шр.
4в
(т31,4вр31 + т32,4вр32 + т33,4вр33 + т34,4вР34 + т3а,4вР3а т4в,СГП Р42^Л +
+ (тСГП ,4вР СГП — (т4в,31 + т4в,32 + т4в,33 + т4в,34
■ + т4в,3а) Р4в)М
Состояние склада готовой продукции:
шРс
(т31,СГП р31 + т32,СГПР 32 + т33,СГП Р33 + т34,СГП Р34 + т3а,СГП Р3а + т 41,СГПР 41 +
+ т 42,СГП Р 42 + т 4 в,СГП Р4р)Л (тСГП ,31 + ШСГП ,32 + ШСГП ,33 +
+ тсГП ,41 + тСГП ,42 + тСГП ,4 в )№> СГП
----+ т +
СГП,34 ^"1СГП,3а ^
Интегрирование данных систем уравнений дает возможность определить вероятности состояний конструкторско-технологической деятельности, а также производства как функции времени в каждый момент времени.
Литература
1. Суворов М.Д. Интеграция моделей проектирования летательных аппаратов и их систем: учебное пособие / М.Д. Суворов. - Омск: Изд-во ОмГТУ, 1999. - 168 с.
2. Инженерный анализ в САПР. Конспект лекций / Одинец М.Н., Романова И.В., Суворов М.Д., Янишевская А.Г. - Омск: Изд-во ОмГТУ, 2009. - 96 с.
3. Вентцель Е.С. Исследование операций: задачи, принципы, методология: учебное пособие / Е.С. Вентцель . - 5-е изд., стер . - М.: КноРус, 2010. - 192 с.