0, если Sm+1J > SmaXe+i
» если QPmjM>Sm,i иSm+l./<S^m+i ‘
®т,і ’ -ели Qnt'i’At^.Sra f и Sm+ij< S№(tl,
где Smj — значение длины очереди (заполнения накопителя) на входе /п-й ТС на j-м промежутке времени:
$т,і~ $т,і—1 +(Qm—1,/—1 > (37)
Q‘m — возможна^ производительность m-й ТС на /-м промежутке времени:
Q°mJ= Z zm,nj ; (38)
л=1
[ 1, если tj < t'm_ и t'm „+ C
Zm’a-^ 0, если (39)
I . пг, я
с начальными условиями:
Qo(/_1 = Q.1 для / = 1,2, . . . , nt. (40)
В начальный момент времени (при / = 0) длину очереди на входе каждой m-й ТС операции в общем случае выбирают из условий обеспечения номинального режима работы:
$
(41)
т,о=Зн-А<, для т=1,2, . . . , ЛІили задают в виде исходных данных:
Sm,0 =Sm ’ Ддя « = 1.2, . . , Л1 , (42)
где Sm— исходное заполнение m-го накопителя.
В случае необходимости метод параметрической рандомизации можно использовать для определения других показателей надежности. Так например, средний объем продукции V, изготовляемой ТС за время t0 и его среднее квадратическое отклонение оа определяют по выражениям:
Пример. Определить вероятность выполнения сменного задания ТС. технологического процесса с поточной организацией производства подшипниковых щитков электродвигателей переменного тока.
Технологический процесс состоит из трех операций:
фрезерной, выполняемой на продольно-фрезерном станке модели 6Г610;
токарной, выполняемой на токарном полуавтомате модели '.286;
сверлильной, выполняемой на вертикально сверлильном полуавтомате модели 2170 М.
Исходные данные:;
время выполнения задания ta— 8 ч;,
объем задания Vo = 7?5O шт.;
номинальная производительность ТС технологического процесса (?,,= =4’00 шт./ч;
количество ТС операций по технологическому процессу Л1 = 3, причем объем накопителя на входе каждой из них не ограничен;
количество единиц оборудования в каждой ТС операции JV1 = ^=W3=1;
номинальная производительность ТС операций: Qi=150 шт..'ч, Q2= 130шт./ч, Q3=200i шт./ч;>
вероятность отказа элементов (операций) ТС: . <?| = 0,1, <?2=О,0О1, 93=0.01;
для ТС первой и второй операций имеет место единственный способ восстановления работоспособности (nB1=ns„=l) с временами восстановления
= 1 ч и tB. =6 ч соответственно;
для ТС третьей операции имеет место два способа восстановления (пВз = 2) с временем восстановления Гв =2 ч и =-1 ч и вероятностями применения Рв31 = 0,» и Рв3 2 =0,2.
Требуемая точность оценки е=О,0О1..
Вероятность выполнения задания определяем в следующей последовательности.
ГО. 1.1. Подразделяем множество рассматриваемых реализаций на подмножества Av по п. 2.2,.6, при этом число подмножеств вычисляем по формуле (213).
3
NA = 2m=1 =23=8 .
. Определяем в соответствии с п. 2.?.5 векторы, характеризующие подмножества Л, , соответствующие им вероятности по выражению (21) и ранжируем подмножества в порядке убывания величин Р,, ■ Результаты вычислений сведены в табл. L
Таблица 1
Индекс подмножества . V |
Индикатор отказа элемента (операции) ТС |
р |
||
X, |
Л'3 |
|
||
1 |
0 |
0 |
0 |
0,890109 |
2 |
1 |
0 |
0 |
0,09890.1 |
3 |
0 |
0 |
1 |
0,008991 |
4 |
1 |
0 |
1 |
0,0.00999 |
5 |
0 |
г |
п |
0,0'00.891 |
6 |
1 |
г |
0 |
0,000099 |
7 |
0 |
11 |
1 |
0,'000009 |
8 |
1 |
1 |
1 |
0,000001 |
Примечание. В табл. 1 и далее по тексту примера в обрзначениях величин xmn, /тл, tBm n и индекс п, соответствующий номеру рассмат
риваемой единицы оборудования в т-к ТС операции опущен, так как по условию задачи каждая операция содержит одну единицу оборудования.
Исключаем из рассмотрения в соответствии с п. 2.2.6 последние подмножества А.. .Так как
8 8
У Pv= 0,001 = £< 2 Р. =0,00199 ,
то количество исключаемых подмножеств
Nn = 4.
Задаемся дискретом времени в соответствии с п. 2.2.2.3.
* Определяем количество точек разбиения интервала (0, t0) по формуле (30)
Подразделяем подмножества Av (из числа оставшихся) на подмножества Д по п. 2.2.7!. Число таких подмножеств в каждом подмножестве A., (v=l, .... 4) вычисляем по формуле (25);
(«B1)0 («B2)0 (v)0'- ы і=і ;
NA = ("в )1(«в )’(«в )°=Ы-1 = 1 ; rig '01' ' t>2' 4 03 • ’
ХА, - (пВі)«-(Пві)0-(пВі)і=Ы-2=2 ;
КА= (%)1-('Івї),,'('1в3)1 = і-1'2=2 •
Определяем количество реализаций /? в подмножествах Д„,ч по выражениям (28) и (29) и вероятность воспроизведения реализации по выражению (26).
Результаты сведены в табл. 2.
Таблица 2
Индекс подмножества |
Индекс отказа элемента (операции) ТС |
Время восстановления элемента (операции) ТС |
Количество отказывающих элементов (операций) k |
Количество реализации R АР- |
Вероятность воспроизведения реализации Р V. р. |
|||||||||
|
* |
хі |
Х2 |
Хз |
ч |
|
t* в.з |
|
|
|
||||
1 |
1 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
1 |
0,890109 |
||||
2 |
1 |
1 |
0 |
0 |
і |
0 |
0 |
1 |
8 |
0,008901 |
||||
3 |
I |
0 |
0 |
1 |
0 |
0 |
2 |
1 |
8 |
0,0071928 |
||||
3 |
2 |
0 |
0 |
1 |
0 |
0 |
4 |
1 |
8 |
0,0017982 |
||||
4 |
1 |
1 |
0 |
1 |
і |
0 |
2 |
2 |
64 |
0,’0007992 |
||||
4 |
2 |
1 |
0 |
1 |
і |
0 |
4 |
2 |
64 |
0,0001998 |
2.2.10.1.6. Определяем индикаторные функции для каждой реализации Результаты вычислений по выражениям (44), (40), (39), (38), (ЗО), (21) и
(35) для реализации из подмножества Д4ДС параметрами (по п. 2.2.8.1).Номер шага І |
Момент времени |
Длина очереди на входе операции |
zl./ |
Z2,i |
|||||
|
S2,/ |
S3,/ |
|||||||
1 |
0 |
100 |
100 |
100 |
1 |
1 |
|||
2 |
1 |
100 |
100 |
100 |
1 |
1 |
|||
3 |
2 |
100 |
100 |
100 |
1 |
1 |
|||
4 |
3 |
100 |
100 |
100 |
0 |
1 |
|||
5 |
4 |
200 |
0 |
100 |
1 |
1 |
|||
6 |
5 |
150 |
150 |
100 |
1 |
1 |
|||
7 |
6 |
100 |
170 |
230 |
1 |
1 |
|||
8 |
7 |
100 |
140 |
160 |
1 |
1 |
Z3.i |
Возможная производительность операции |
Фактическая производительность операции |
Объем изготовленной продукции V (tp |
||||||
|
0^ |
0° у3,/ |
4/ |
«2,/ |
$з,/ |
|
|||
1 |
150 |
130 |
200 |
100 |
100 |
100 |
100 |
||
1 |
150 |
130 |
- 200 |
100 |
100 |
100 |
200 |
||
1 |
150 |
130 |
200 |
100 |
100 |
100 |
300 |
||
1 |
0 |
130 |
200 |
0 |
100 |
100 |
400 |
||
0 |
150 |
130 |
0 |
150 |
0 |
0 |
400 |
||
0 |
150 |
130 |
0 |
150 |
130 |
0 |
400 |
||
1 |
150 |
130 |
200 |
100 |
130 |
200 |
600 |
||
1 |
150 |
130 |
200 |
100 |
130 |
160 |
760 |
Таблица З
Стр. 28 ГОСТ 27.204—83
xt=9, х8=1, 4=3, 4=4, 4 = 1, 4t=o, 4s=2
сведены в табл. 3.
В соответствии с условием (34)
V(<e)=76O>Vo=75O,
индикаторная функция для данной реализации
<р=Д.
Аналогичным образом вычисляем индикаторные функции других реализаций
Для подмножества А},} значение индикаторной функции <р='1.
Вычисленные значения индикаторных функций для остальных рассматрива емых подмножеств сведены в табл. 4—<8.
Таблица 4
Значение индикаторной функции подмножества Д21 для моментов отказа 1-го элемента (операции)
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
1 |
1 |
1 |
1 |
С |
0 |
0 |
0 |
Таблица 5
Значение индикаторной функции подмножества А3д для моментов отказа 3-го элемента (операции)
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
1 |
1 |
1 |
1 |
I |
0 |
0 |
0 |
Таблица 6
Значения индикаторной функции подмножества А 32
для моментов отказа 3-го элемента (операции)
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
1 |
8 |
0 |
0 |
0 |
0 |
0 |
0 |
Таблица 7
Значение индикаторной функции подмножества Л 41
Момент отказа 1-го элемента (операции) |
Момент отказа 3-го элемента (операции) |
|||||||
1 |
2 |
3 |
4 |
5 |
G |
7 |
8 |
|
1 |
1 |
1 |
1 |
1 |
0 |
0 |
0 |
0 |
2 |
1 |
1 |
1 |
1 |
0 |
0 |
0 |
0 |
3 |
1 |
1 |
1 |
1 |
0 |
0 |
0 |
0 |
4 |
1 |
1 |
1 |
1 |
0 |
0 |
0 |
0 ■ |
5 |
1 |
1 |
1 |
1 |
0 |
0 |
0 |
0 |
6 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
7 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
8 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
|
|
|
|
|
|
|
Та б |
лица 8 |