Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Игнатов, В. А. Статистическая оптимизация качества функционирования электронных систем

.pdf
Скачиваний:
4
Добавлен:
21.10.2023
Размер:
9.05 Mб
Скачать

 

 

 

 

Таблица 4-1

к

V 4

Ат.я , ч

Р У 6 -

"<CW - Р«6-

1

87,0

63

31

31

2

150

56

17,12

48,12

3

206

50

16,22

64,34

4

256

48

15,78

80,12

5

304

46

23,98

104,1

6

350

44

18,93

123,03

7

394

43

19,29

142,32

Таким образом, при монотонно возрастающей огра­

ниченной

интенсивности отказов

оптимальная

частота

ѵй.=Літ_1/1 KP также возрастает,

установившееся

значе­

ние

limvft

определяется

установившимся

значением

 

&->оо

 

 

 

 

 

 

Л(/) и параметром х.

 

 

 

 

 

В инженерных расчетах достаточно оценить Лгі по

формуле (4-19), найти

 

 

 

 

 

 

 

lim Д-Cfc = ]/2х мин (а0, а,)-1

 

(4-22)

 

 

 

fe->00

 

 

 

 

и время

Т д установления

Л (0 .

например,

по

уровню

0,9

мин (ао,

аі). Тогда общее число т .проверок работо­

способности до наступления установившегося режима определяет неравенство

ТА

< т <

т А

1

(4-23)

[ Atj

 

Ига ДтЛ

 

 

 

&-»оо

 

в котором используется целая часть отношений. Полагая ТАс^Т0— ІО3 ч, получим 11 < m<31, т. е. примерно че­

рез 20 проверок можно пользоваться (4-22). Расчеты переходного режима оптимального KP значительно упро­ щаются, если воспользоваться квазиоптимальной перио­ дичностью.

В этом случае

Д т ^

Дт, -f Hm Дт„

----- --------- (4-24)

Итак, в данном параграфе получены приближенные аналитические соотношения для расчета переходного и

140

установившегося режимов оптимального KP. Отличи­ тельной особенностью является то, что в них учтены взаимосвязанные внезапные и постепенные отказы, а также дисперсия погрешностей измерений.

4-3. Оптимальная периодичность профилактической замены элементов

Вработах [Л. 17, 18, 54—56] для отдельных случаев показано, что оптимальная периодичность ПЗ элементов (электровакуумные и полупроводниковые приборы, маг­ нетроны, клистроны, модули и т. п.) существует не всегда.

Вэтом параграфе формулируются достаточно общие условия существования и единственности решения транс­ цендентных уравнений, определяющих оптимальную пе­ риодичность ПЗ, о учетом, влияния внезапных отказов, полного и неполного восстановления свойстз элементов после АР и ПЗ. Простой вид этих условий [Л. 53], в ко­ торые входят параметры ФР времени безотказной рабо­ ты элементов и показатели затрат на АР и ПЗ, позволя­ ет легко убедиться в целесообразности поиска оптималь­

ных решений.

Обозначим: F{t) — ФР времени безотказной работы элемента, ср— средние затраты на один АР; сп — сред­ ние затраты на одну ПЗ; %— периодичность ПЗ; с(т) — удельные эксплуатационные расходы (УЭР). Обычно ср и Си измеряют в рублях или часах, а с(т) — руб-ч~1 или в относительных единицах. В последнем случае по физи­ ческому смыслу с (г) является коэффициентом простоя.

Условия существования и характеристики оптималь­ ной ПЗ, итерационные алгоритмы отыскания т0Пт найдем из решения задачи в. постановке, ставшей уже традици­ онной [Л. 15, 17, 18]: F(t), сѵ « са известны, требуется определить т0Пт, минимизирующую математическое ожи­ дание с (т).

Математическое ожидание УЭР

 

A l[g(T)]^Cpf (; ) + Cn[1~ f (т)1.

(4-25)

5

 

141

Дифференцирование выражения (4-25) по х и при­ равнивание производной нулю дает:

Л (т) J Р (t) dt + P (х) = Ср (Ср - сиу ».

(4-26)

о

 

Если решение уравнения (4-26) существует и являет­ ся единственным, то существует и оптимальная пери­ одичность ПЗ, минимизирующая выражение (4-25).

Докажем условия существования и единственности решения (4-26) для ВФИ- и ОВФИ-распределений.

Теорема 4-1. Если F(t) есть ВФИ-распределение с не­ прерывной плотностью и

1<Ср(ср—Сп)-1< ° ° ,

(4-27)

то на интервале [0, оо] решение уравнения (4-26) суще­ ствует.

Доказательство. Обозначим левую часть (4-26) через Ь(т) и исследуем предельное поведение этой функции. Из условия теоремы следует, что она непрерывна, по­ этому .

lim L (х) =

1, lim L (х) = lim А (х) Г Р (t) dt -f- lim P (x) = oo. ]

T->-0

 

1->CO

T-»00

Z~>00 ]

Так как при изменении х от

0 до оо L(t) изменяется

от 1

до

оо, то

при выполнении

условия (4-27) функция

L{т)

пересечет

уровень

ср(ср—сп) _1,

что и требовалось

доказать.

Из теоремы следует, что при ср^ с п оптимальной ПЗ не существует, что вполне соответствует физическим представлениям. Для случая ср = сп этот результат, полу­

ченный различными другими способами,

в настоящее

время широко известен.

 

 

Теорема 4-2. Если

 

 

FXtj= 1 -

,

(4-28)

есть ОВФИ-распределение (2-87) и

 

 

1<ср(ср—сп) - 1< 7 0 мин

(ао, аі),

(4-29)

где То — среднее время безотказной

работы, то решение

(4-26) существует.

 

 

142

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

lira L (т) == lim Л (т) lim [P{t)dt = мин (а0, аt)T 0. t->co т~>ао «J

Пример 4-2. Пусть время безотказной работы элемен­ та имеет гамма-распределение с параметрами п и а, тог­

да 7'о = «а_1 и аТ0=п.

условие

Следовательно, для гамма-распределения

(4-29) принимает особенно простой вид:

 

1<Ср(Ср—Сп)_1< я .

(4-30)

При п = 1 и Ср^Сд это условие не выполняется. Теорема 4-3. Если Л(/) ВФИили ОВФИ-распределе-

ний есть дифференцируемая функция и решение (4-26)

существует, то оно является единственным.

по т и по­

Доказательство. Продифференцируем Ь{т)

сле несложных преобразований получим:

 

L ' (т) = Л' (т) I Р (t) dt > 0.

(4-31)

Следовательно, L (т) является монотонно возрастаю­ щей функцией, которая при изменении %от 0 до оо толь­ ко один раз пересечет уровень ср(ср—сп)_1, что и требо­ валось доказать.

Когда существование и единственность решения (4-26) установлены, можно приступить к отысканию т0птОбычно эту процедуру рекомендуют выполнять графи­ ческим методом. Сократить объем вычислений и гаран­ тировать требуемую точность результатов можно за счет использования итерационного алгоритма отыскания Топт, построенного по методу Ньютона [Л. 57]:

ті+1= Т г+[С р ( С р —

С п ) _І— L (Ті) ] X

 

X [ L ' ( X i ) ] - K

t= 0, 1 ....

(4-32)

где для ВФИ-распределений Xo=kcaT0{cpСо)-1; для

ОВФИ-распределений тo = kcaT0[{cp—сп) (ссГ0—1]_1,

а =

= 1ітЛ (т); £ = 2-4-4 — коэффициент, учитывающий

по-

1-» С О

грешности линеаризации L ( t ) . Так как метод Ньютона обеспечивает сходимость более быструю, чем сходимость геометрической прогрессии с любым знаменателем, мень-

143

шим единицы, то требуемая на практике точность дости­ гается за две-три итерации [Л. 57].

После отыскания т0ПТ целесообразно оценить мини­ мальные УЭР по формуле (4-25) и выигрыш W от опти­

мизации ПЗ. При отсутствии ПЗ lim М[с(т)] = сѵТ ^ 1, еле- т-»со

довательно,

 

W — С ^ ^ [ С р Г *0-— СМин (Т опт)]1'00% —

 

 

 

 

='І1-смин(ТопТ)7’оС-1р]100%.

(4-33)

Пример 4-3, Пусть время безотказной работы элемен­

та имеет гамма-раопределение

с параметрами п = 2; а =

—1,2-10—3 ч_1; показатели ср= 100 руб.; сп=10 руб.

 

Так

как

ср(ср—сп) = 1,111

удовлетворяет

условию

(4-30)

и условие теоремы 4-3

выполняется, то

решение

(4-26) существует и

является

единственным. При

k = \,

 

 

_

ПО-1 667

184 я.

 

 

 

 

 

90(2 -і)

 

 

 

 

 

 

 

 

Применяя

(4-32),

после трех итераций получим

Ті=

= 654 ч, тг=584 ч , тз=574 ч . Вычислительная процедура прекращена, так как последние два значения отличаются менее чем на 2%. Следовательно, тОПт~574. Отметим, что, несмотря на «неудачный» выбор т0 (&=1), итераци­ онный процесс сошелся за три итерации.

Найдем ТЭХ оптимальной ПЗ: М [с5ШН (574)] = = 0,0441 руб -1;

№= 26,5%; ТсштГ'1=0,344.

Целесообразность оптимизации ПЗ очевидна, так как она сокращает средние УЭР примерно на 26,5%.

Перейдем к условиям, характеризующим влияние вне­ запных отказов на решение (4-26). Выберем ОВФИ-рас- пределение типа (4-28), в котором для простоты поло­ жим Хо-=Хі=к, ро=\. Используя неравенства (4-29) и (4-30), получим следующие условия существования опти­ мальных ПЗ при взаимосвязанных внезапных и посте­ пенных отказах:

х< Ср

ъСи

— < • Ср -

Ѵі

— <■ V

Сц

% .

ао <

а6

 

Сп

1^

 

1

аі <

а0;

(4_34)

 

 

1.

ао — аг

 

 

19

 

144

 

Неравенства

 

(4-34) оп­

 

------1

 

 

 

 

Л

ределяют условия, при ко­

 

исто'

 

V

 

 

торых

наличие

 

внезапных

 

\\

 

г

отказов не мешает опреде­

 

Л

/

 

лению

 

Том-

Невыполнение

5,5

 

\

I

 

 

 

неравенства

(4-34)

говорит

 

 

 

.

I

 

1 /

 

о том, что из-за большой ве­

 

 

 

 

 

роятности

внезапных

отка­

 

 

 

Ч

I

 

 

 

зов

 

определение т0Пт с уче­

¥

 

 

 

 

 

 

том

совместного

появления

 

 

Wcftonm

 

внезапных

и

 

постепенных

 

 

 

 

 

 

 

0,2

 

___ U ____ Ѵо*

отказов невозможно.

 

 

 

0,Н

0,6

0,8

 

 

Если

неравенство

(4-34)

 

Рис. 4-1. Графики зависимости

не выполняется, то, так как

 

УЭР от нормированного време­

внезапные отказы не влияют

 

ни при іЯр_1о=0 (1)\ 0,2 (2);

на величину т0ПТ, а изме­

 

0,4 (3).

он-

 

 

г|*___

няют ЛИШЬ Смин (Топт)

15],

 

0 .2

0 ,6

в

уравнении

 

(4-28) целе­

 

V ;

 

V

 

 

 

Рр

сообразно

все

интенсивно­

0,8

 

 

 

 

 

сти внезапных отказов при­

 

 

 

 

 

 

 

 

равнять

нулю.

К

 

этому

0,6

 

 

 

 

 

 

же приему можно прибе­

 

V

 

 

 

 

гать

и

 

тогда,

когда

из-за

Oft

 

 

 

 

 

влияния

внезапных

отказов

 

 

 

 

 

 

итерационная

 

 

процедура

0,2

 

 

 

 

 

 

(4-32) обладает медленной

 

 

 

 

 

 

сходимостью.

 

Не

следует,

 

 

 

 

 

Од 0п

 

однако,

 

забывать,

что в этом

 

 

10

15

20

 

случае

все

характеристики

 

Рис. 4-2.

Графики

зависимо-

оптимальных

ПЗ

относятся

 

 

стей тТ-'о от ро

(1,

2) и

только к постепенным отка­

 

СрС ‘п (3).

 

 

 

 

 

 

зам.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4-25)

при

 

На рис. 4-1 приведены графики функции

СрС,

1

 

10,

Т),

 

10' !

у - 1, Tj1=

l,5-10_3

я

1

и различных

отношениях

Ятір1

 

Кривая

1

построена

при

 

Ят]~’ = 0

(внезапные отказы отсутствуют или не учитываются), кривая 2 —при Хщ1= 0 ,2 , кривая 3 — при Ят^"1= 0 ,4 . Из условия (4-34) следует, что минимум функции (4-25) неразличим в данном случае при Яі)^1> 7,5. Как видим,

влияние внезапных отказов проявляется в том, что мини­ мум функции (4-25) как бы «размывается», однако, что очень важно, величина т0ПТне изменяется. Сростом числа внезапных отказов растет лишь величина ЛГІСмщ^Топт)]-

Ю—385

145

На рис. 4-2 приведена типичная зависимость относи­

тельной величины оптимальной периодичности тД“ 1 от

величины СрС_1п (кривая 3). Как и следовало ожидать, рост затрат на АР влечет за собой уменьшение тТ~^, т. е. в такой ситуации целесообразна более высокая частота проведения ПЗ.

Для случая роф і использование неравенств (4-29) и (4-30) дает следующие условия существования опти­ мальной ПЗ для вероятности полного обновления свойств элемента после АР или ПЗ:

Р о >

CpCt| ttp (Ср Сд) '

а » < а.;

 

(Ср— Си) (а, — «о)

 

Ро>-

_____ ^дао_______ .

аі < а0;

(4-35)

(Ср — с„) (а, — а 0) ’

 

Ро> ~ ъ г ^ 7 ’ а» = а>-

Они говорят о том, что проведение ПЗ нецелесооб­ разно, если вероятность полного обновления элемента

падает ниже

критической,

определяемой затратами на

АР, ПЗ и параметрами

процесса ухудшения

качества

элемента.

 

 

что уменьшение

ро

приводит

Анализ показывает,

к увеличению

tT

т.

е.

если отсутствует

возможность

полностью обновить элемент при ПЗ, то ПЗ следует про­

 

 

 

водить реже.

Если ро мень­

 

 

 

ше

критического

значения,

 

 

 

определяемого

условиями

 

 

 

(4-35), то проведение ПЗ не­

 

 

 

целесообразно

и

приводит

 

 

 

к убыткам. При

увеличении

 

 

 

СрС~1а критическое значение

 

 

 

ро падает. Таким образом,

 

 

 

при

невысокой относитель­

 

 

 

ной

средней

стоимости ПЗ

 

 

 

можно проводить чаще даже

 

 

 

в случае слабого обновле­

Рис. 4-3.

Графики зависимо­

ния элементов. Выигрыш от

оптимизации ПЗ

тем

боль­

стей выигрыша от оптимиза­

ший, чем больше

р0 и СрС-1п.

ции

от

ро при срс-'п = 5 (/);

10

(2) .

 

При

больших

СрС-1д

значи-

146

тельно расширяется диапазон изменения р0, в котором оптимизация целесообразна.

На рис 4-2 (кривые 1, 2) и рис. 4-3 показаны типич­

ные зависимости -сТ ~1и W от ро. построенные при гіо=

= г)і = 2 • 10-3 ч~1, Яо = Лі = 0. Кривые 1 рассчитаны при СрС_1п= 5, кривые 2 — при срс_1п=10. Выигрыш от опти­ мизации существенно зависит от степени обновления элементов после АР или ПЗ.

На рис. 4-4 приведены графики

функции М[с {хТ~1)]

при ро= 1; 0,8; 0,6;

0,4

(соответственно кривые 14),

СрС-1п=10. Они наглядно по­

 

 

 

 

казывают, как влияет умень­

 

 

 

 

шение ро на УЭР, пунктир­

 

 

 

 

ная

кривая

характеризует

 

 

 

 

смещение

экстремума.

От­

 

 

 

 

метим, что (в области т<'т0Пт

 

 

 

 

производная функция значи­

 

 

 

 

тельно больше, чем в обла­

 

 

 

 

сти т>ТоптПоэтому, если

 

 

 

 

по организационным

причи­

 

 

 

 

нам

ПЗ

нельзя проводить

 

 

 

 

своевременно, ее лучше про­

 

 

 

 

вести

позднее — это

приво­

 

 

 

 

дит к меньшим УЭР.

 

 

Рис.

4-4. Графики зависимости

Итак, в этом параграфе

УЭР

от

нормированного

вре­

доказаны условия существо­

мени

при

ро=1 (/); 0,8

(2);

вания и единственности ре­

0,6 (3); 0,4 (4).

 

шения (4-26),

позволяющие

 

 

 

 

оптимизировать периодичность ПЗ элементов с целью обеспечения минимума математического ожидания УЭР,

построен итерационный алгоритм определения х и

проиллюстрировано влияние внезапных отказов, непол­

ного

обновления элементов и стоимостей АР и ПЗ на

на и

и W .

ОГт

Алгоритм расчета ТЭХ оптимальной ПЗ следующий: выбирают целевую функцию, по которой возможна и це­ лесообразна оптимизация; статистическими методами определяют интенсивности ухудшения параметров, ин­ тенсивности внезапных отказов, показатели затрат сѵ, са, вероятность р0; проверяют выполнение необходимых условий (4-27), (4-29), (4-30), (4-34), (4-35), с помощью

147

рекуррентных соотношений

(4-32) и формул (4-25),

(4-33) рассчитывают искомые характеристики.

Достоинство описанного

алгоритма

заключается

в том, что оптимизация имеет наглядный

физический

смысл, использует относительно простые исходные дан­ ные и рекуррентные соотношения, обеспечивающие за­ данную точность. Выигрыш от оптимизации ПЗ оцени­ вается в целом — с учетом взаимосвязанных внезапных и постепенных отказов.

4-4. Оптимальная периодичность профилактического обслуживания систем

При обслуживании сложных устройств (например, приемников и передатчиков РЛС, блоков РНС и т. п.) АР уже не приводит к полному обновлению, поэтому ха­ рактер интенсивности отказов после АР не изменяется. Это обусловлено тем, что из-за большого числа элемен­ тов в устройстве замена одного из них при АР слабо сказывается на суммарной интенсивности отказов. В то же время при ПО таких устройств производится или полная замена всех износившихся элементов, или такая их регулировка, что интенсивность отказов из-за износа элементов падает до нуля — происходит как бы полное обновление устройства после ПО. Ясно, что для слож­ ных устройств, как правило, ср< с п.

Для устройств [Л. 18, 53, 56]

(4-36)

поэтому для минимизации УЭР необходимо решать трансцендентное уравнение

Л(т)т + 1пР(т) = С п С_ 1р.

(4-37)

В работе [Л. 18] доказывается, что необходимыми условиями существования и единственности решения (4-37) являются непрерывность, дифференцируемость, строгое возрастание и неограниченность интенсивности отказов на интервале [0, оо]. Последнее условие на прак­ тике для сложных устройств не выполняется, но, как оказалось, оно и не требуется.

148

Теорема 4-4. [Л. 5'3]. Если распределение (2-87) есть ОВФИ-распределение и

О

СдСр < ln k0,

я0 я,;

 

■ 0

< с пс~'< Л п klt

а0> а ,;

(4-38)

. 0 < с ис“ 1< оо, а0 = а,,

 

то в интервале (0, оо] решение уравнения (4-37) сущест­ вует.

Доказательство. Обозначим левую часть уравнения (4-37) через ЕДт) и исследуем поведение этой функции. Она является непрерывной, так как представляет сумму непрерывных функций. При т-И) Еі(т)-И). Отыскание предела ЕДт) при т-^-оо приводит к неопределенности вида «оо—оо». Применим подстановку Ф(т) = е Ефт).

Из

непрерывности показательной

функции

следует,

что lim Е, (т) — In Е0, если

1ітФ(т) = £0. Подставим Е,(т)

1-»0О

 

t->00

 

 

 

в выражение для Ф(т), получим:

 

 

 

 

lim Ф (т) =

lim \ек w т Р (т)].

(4-39)

 

 

Т-ИЭО

Т -* С С

 

 

 

Теперь

нетрудно показать,

что при а0О ,

1ітФ(т) =

= А0;

при

а0 > а , lim Ф (т) = &,;

при

 

1->00

а0 — а, lim Ф (т) = сю.

Ясно,

 

Т ->00

 

 

Т ->-00

что при выполнении условий теоремы функция Li (т)

пересечет уровень спс~\ что и требовалось доказать.

Теорема 4-5. [Л. 53]. Если решение уравнения (4-37) для ОВФИ-распределения существует, то оно является единственным.

Доказательство. Дифференцируя ЕДт) по т, получим

L'і(т) = Л / (т)т>0, следовательно, ЕДт)

является

моно­

тонно возрастающей функцией, которая

при изменении

т от 0 до оо только один раз пересечет

уровень

спс~‘р,

что и требовалось доказать.

 

(2-9)

Для обобщенного ОВФИ-распределения типа

функция ЕДт) ограничена величиной lnDj„, где

индекс

/ определяет условие

 

 

a.j — мин (а0, <*„_,).

 

(4-40)

Интересно отметить, что т0Пт определяется из реше­ ния (4-37) и в том более общем случае, когда средние

149

Соседние файлы в папке книги из ГПНТБ