Программа модели. ; Модель изготовления деталей

; Модель изготовления деталей. Обратная задача

; Задание исходных данных

Tn_ EQU 35; Среднее время между поступлениями заготовок

T1 EQU 30; Среднее время выполнения 1-й операции, мин

T2 EQU 25; Среднее время выполнения 2-й операции, мин

T3 EQU 35; Среднее время выполнения 3-й операции, мин

To3 EQU 6; Среднеквадратическое отклонение времени выполнения 3-й операции, мин

Tk1 EQU 4; Среднее время контроля после 1-й операции, мин

Tk2 EQU 5; Среднее время контроля после 2-й операции, мин

Tk3 EQU 15; Среднее время контроля после 3-й операции, мин

Tok3 EQU 2; Среднеквадратическое отклонение времени контроля после 3-й операции, мин

Tk EQU 8; Среднее время окончательного контроля, мин

q1_ EQU.12; Доля брака после 1-й операции

q2_ EQU.15; Доля брака после 2-й операции

q3_ EQU.10; Доля брака после 3-й операции

q4_ EQU.80; Доля окончательного брака

Det EQU 4; Количество деталей, которые нужно изготовить

; Описание функции времени подготовки заготовок для детали

Pod FUNCTION RN10,D6

.05,10/.18,14/.34,21/.56,22/.85,28/1,25

; Сегмент имитации изготовления деталей

GENERATE (Exponential(1,0,Tn_)); Источник заготовок

; Подготовка заготовок для деталей

QUEUE Pod; Встать в очередь

SEIZE Pod; Начать подготовку заготовки

DEPART Pod; Покинуть очередь

ADVANCE (Exponential(34,0,FN$Pod)); Имитация подготовки

RELEASE Pod; Закончить подготовку заготовки

; Имитация выполнения 1-й операции

DCount ASSIGN 1,1; Код 1 в Р1-деталь проходит первый раз

ASSIGN 2,1; Код 1 в Р2-признак 1-й операции

Oper1 QUEUE P2; Встать в очередь

SEIZE Konveer1; Начать первую операцию

DEPART P2; Покинуть очередь

ADVANCE (Exponential(23,0,T1)); Имитация 1-й операции

RELEASE Konveer1; Закончить 1-ю операцию

ADVANCE (Exponential(23,0,Tk1)); Контроль 1-й операции

TRANSFER q1_,,Sboi; Отправить брак на пункт контроля

; Имитация выполнения 2-й операции

ASSIGN 2,2; Код 2 в Р2-признак 2-й операции

Oper2 QUEUE P2; Встать в очередь

SEIZE Konveer2; Начать вторую операцию

DEPART P2; Покинуть очередь

ADVANCE (Exponential(23,0,T2)); Имитация 2-й операции

RELEASE Konveer2; Закончить 2-ю операцию

ADVANCE (Exponential(23,0,Tk2)); Контроль 2-й операции

TRANSFER q2_,,Sboi; Отправить брак на пункт контроля

; Имитация выполнения 3-й операции

ASSIGN 2,3; Код 3 в Р2-признак 3-й операции

Oper3 QUEUE P2; Встать в очередь

SEIZE Konveer3; Начать третью операцию

DEPART P2; Покинуть очередь

ADVANCE (Normal(20,(T3#K1),(To3#K1))); Имитация 3-й операции

RELEASE Konveer3; Закончить 3-ю операцию

ADVANCE (Normal(20,Tk3,Tok3)); Контроль 3-й операции

TRANSFER q3_,,Sboi; Отправить брак на пункт контроля

TRANSFER,Met2; Готовые детали

; Сегмент имитации работы пункта контроля

Sboi QUEUE Kont; Встать в очередь на пункт контроля

SEIZE Kontr; Занять пункт контроля

DEPART Kont; Покинуть очередь на пункт контроля

ADVANCE (Exponential(23,0,Tk)); Окончательный контроль

RELEASE Kontr; Освободить пункт контроля

TRANSFER q4_,,EndOper; Отправить в окончательный брак

TEST E P1,1,EndOper; Если второй раз, то в окончательный брак

ASSIGN 1,2; Код 2 в Р1-деталь проходит второй раз

Met1 TRANSFER,(Met1+P2)

TRANSFER,Oper1; На повторное выполнение 1-й операции

TRANSFER,Oper2; На повторное выполнение 2-й операции

TRANSFER,Oper3; На повторное выполнение 3-й операции

EndOper TERMINATE; Счет брака

; Сегмент завершения моделирования и расчета результатов

Met2 TEST L X$Prog,TG1,Met3; Если условие выполняется, то

SAVEVALUE Prog,TG1; X$Prog=TG1 счетчику завершений

SAVEVALUE NDet,0; Обнуление X$NDet

Met3 SAVEVALUE NDet+,1; Счет количества готовых деталей

TEST E X$NDet,Det,Ter1; Если готово Det деталей, зафиксировать один прогон

TEST E TG1,1,Met4; Если содержимое счетчика завершений равно 1,

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

SAVEVALUE Brak,(INT(N$EndOper/X$Prog)); Количество забракованных деталей, шт.

SAVEVALUE DoljaBrak,(X$Brak/(X$Brak+Det)); Общая доля брака

SAVEVALUE DoljaDet,(Det/(X$Brak+Det)); Доля готовых деталей

SAVEVALUE AC3,(AC1-X$AC2)

SAVEVALUE TDet,((X$AC3/X$Prog)/60); Среднее время изготовления Det деталей, час

SAVEVALUE SDet,((X$TDet/N_)#60); Среднее время изготовления одной детали, мин

SAVEVALUE AC2,AC1

SAVEVALUE X$Prog,0 Met4 SAVEVALUE NDet,0; Обнуление X$NDet

TERMINATE 1; Вычитание из счетчика завершений 1

Ter1 TERMINATE; Вывод вспомогательных транзактов

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

Поскольку результатом моделирования является оценка математического ожидания времени TDet изготовления Det деталей, то в ее вычислении используется абсолютное модельное время АС1 (системный числовой атрибут). При проведении дисперсионного анализа встроенный генератор эксперимента имеет две команды START, а между ними - команда RESET. Команда RESET не влияет на абсолютное модельное время. Поэтому АС1 будет суммой абсолютного модельного времени предварительных прогонов до установившегося режима, обозначим его АС2, и абсолютного модельного времени, пусть АС3, основных прогонов, в ходе которых собирается интересующая нас статистика. Нам для расчетов нужно АС3. Для его получения в программу введены строки:

SAVEVALUE AC3,(AC1-X$AC2)

SAVEVALUE AC2,AC1

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

SAVEVALUE X$Prog,0