![](/user_photo/_userpic.png)
книги из ГПНТБ / Игнатов, В. А. Статистическая оптимизация качества функционирования электронных систем
.pdf
|
|
|
|
Таблица 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 |
(соответственно кривые 1—4), |
|||||||
СрС-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