книги из ГПНТБ / Клемин А.И. Инженерные вероятностные расчеты при проектировании ядерных реакторов
.pdfи при экспоненциальном законе надежности также равна парамет ру закона к.
Выясним физический смысл интенсивности отказов на примере экспоненциального закона (6.16). Обозначим Р у (t, t + At) услов ную вероятность, что канал реактора не откажет на интервале вре
мени от t |
до t + At при условии, что до момента t отказов не было. |
||
Согласно теореме умножения вероятностей (см. табл. 1.1), |
|||
|
Р {i + At) = |
Р (/) • Р у (t, |
t + At), |
где P(t + |
At) — безусловная |
вероятность для канала проработать |
|
безотказно период времени |
от 0 до t + |
At; P(t) —• аналогично от |
|
Одо t. Следовательно, условная вероятность отказа канала в интер вале от t до t -f- At
Р у |
(/, t + At) = 1 — P y |
(/, t + At) = |
l — p (/ + |
At)/P (/). |
(6.26) |
||||
При |
достаточно |
малом |
At, |
используя |
выражение |
(6.16), |
получаем |
||
|
Р у (t, t + |
At) == 1 — ехр [~к (t + |
AOl/exp (—kt) |
= |
|
||||
|
|
= |
1 — ехр (—kAt) « |
kAt. |
|
|
(6.27) |
||
Таким образом, при малом |
= 1 величина к равна условной ве |
||||||||
роятности для перемонтируемого изделия отказать в единичном
интервале |
времени |
[t, t + At] около момента |
t при условии, |
что |
|||||||||
до момента t |
оно работало исправно. |
|
|
|
|
|
|||||||
Заметим, |
что безусловная |
вероятность отказа |
изделия |
в ин |
|||||||||
тервале от t до t + |
At определяется по формуле (2.8): |
|
|
||||||||||
Р (i, t + |
At) = |
Р {t < т < |
t + |
At} = |
\ |
f(f)dtg*f |
(t)-At. |
||||||
|
|
|
|
|
|
|
|
|
і |
|
|
|
|
Указанный смысл |
интенсивность |
отказов сохраняет |
и для ремон |
||||||||||
тируемых |
(восстанавливаемых) |
изделий, |
когда |
рассматривается |
|||||||||
их поведение |
в период до первого (или очередного) |
ближайшего |
|||||||||||
отказа. Обозначим |
такую интенсивность отказов для |
восстанавли |
|||||||||||
ваемых изделий через А-ютк- Важно помнить, |
что в общем случае |
||||||||||||
^іотк Д л я |
ремонтируемого изделия |
не равна |
параметру потока от |
||||||||||
казов к этого изделия. И лишь в одном-единственном случае, |
когда |
||||||||||||
к— const, |
т. е. в период |
нормальной эксплуатации, |
интенсивность |
||||||||||
отказов к10ТК |
= |
^- |
|
из |
формул (6.16), (6.23) |
и (6.26) следует, что |
|||||||
Непосредственно |
|||||||||||||
|
|
|
|
к |
= |
— Р ' (/)/Р (t) = / |
(/)/Р (t). |
|
|
(6.28) |
|||
Эта формула имеет общий характер, она справедлива не только для экспоненциального (6.16), но и для любых других законов надеж ности, для которых к = k(t) Ф const.
Для опытного определения интенсивности отказов k(t) необ ходимо проследить за эксплуатацией многих однотипных изделий
в течение определенного периода времени, начиная от момента пуска і = 0. Опытное значение (точечная оценка) интенсивности отказов
М О |
т (i, t+At) _n(i) |
— n |
+ At) |
(6.29) |
||
|
Atn (t) |
|
Atn |
(t) |
|
|
где n(t) — число работоспособных |
изделий |
в момент /; |
m(t, |
t + |
||
+ At) — число отказавших изделий |
в интервале времени |
(/, / + |
At). |
|||
Основные законы надежности и связанные с ними характери |
||||||
стики. Закон (или функция) надежности P(i) наиболее |
полно |
ха |
||||
рактеризует надежность изделия. |
Зная |
закон, можно |
вычислить |
|||
|
Р и с . 2 1 . З а в и с и м о с т ь и н т е н с и в н о с т и о т к а з о в от в р е |
|
||||||
|
|
мени д л я р а з л и ч н ы х |
з а к о н о в н а д е ж н о с т и |
( ц и ф р ы |
|
|||
|
|
у к р и в ы х — н о м е р а в т а б л . 6.1). |
|
|
||||
интенсивность отказов |
Щ) по формуле (6.28), среднюю |
наработку |
||||||
на отказ по формуле (6.24) и у — процентную |
наработку Ту, ко |
|||||||
торая будет |
реализована в у |
• 100% |
случаях: |
t |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
со |
|
|
|
|
X (t) |
= _ Р ' (0/Р (/); |
Г 0 = |
|' Р (0 dt; |
Р (7%) = |
у, |
||
|
|
|
|
|
о |
|
|
|
где Ту |
— период времени-, |
вероятность безотказно |
проработать |
|||||
который |
равна у. |
надежности P(t), как правило, |
|
|
||||
На практике закон |
не известен. |
|||||||
Приходится по статистическим данным эксплуатации или. испы таний устанавливать форму закона, пользуясь методами матема тической статистики, см. § 4.3. Основные законы надежности и
соответствующие им |
характеристики f(t), |
K(t). и Т 0 ; приведены в |
|
табл. 6.1. Характер |
изменения Ці) показан |
на рис. 2І. |
I |
О с н о в н ы е з а к о н ы н а д е ж н о с т и и с в я з а н н ы е с ними характеристики
с
с
%
1
|
Закон надежности Р (і) |
||
Экспоненциальный |
|
||
|
Р ( / ) |
= е х р ( — |
XI), |
см . ф о р м у л у (3.5) |
|
||
Н о р м а л ь н ы й |
|
|
|
Р ( 0 |
|
t-T„ |
(Те>3а). |
= 0 , 5 - Ф |
|||
см . |
ф о р м у л у |
(3.10) |
|
Плотность распределения времени |
Интенсивность отказов |
безотказной работы [ (/) |
М О |
X ехр ( — Xt). с м . ф о р м у л у (3.4) |
X |
|
|
/ - Г 0 \ * |
|
|
ехр |
при |
і < 0, / ( 0 = 0 , так как |
о V 2л 0,5 — Ф |
Т0 |
> З а , см . ф о р м у л у (3.9) |
|
Усеченный нормальный |
Г |
1 |
(t—T\2l |
|
|
|||
|
t — T |
ехр |
|
|
|
|
|
|
0,5 — |
|
|
|
|
|
|
||
Ф |
о~[/2я |
[0, 5 + |
Ф (Т/а)\ |
ехр |
2 |
|||
Р ( 0 = |
|
|||||||
|
|
|
|
|
|
|||
|
при t |
< 0, |
/ (/) = |
0 |
|
|
||
0,5 + |
Ф |
а У 2т 0,5 - - ф ( - н |
||||||
|
|
|
|
|||||
|
|
[отрицательная |
ветвь |
/ (/) |
||||
(любое |
Г > 0) |
|
отсечена] |
|
|
|
||
|
|
|
|
|
|
|||
Т а б л и ц а б.(
Средняя наработка на отказ Го
1
X
а е х р
2 ^ а 1
У " 2 л
с
сЗакон надежности Р (О
4 Л о г а р и ф м и ч е с к и нормальный
(Ы—а\
Р ( 0 = 0 , 5 - Ф ^ — — J,
см . ф о р м у л у (3.16)
5 Гамма |
|
|
|
|
а—1 2 |
с м . |
ф о р м у л у |
(3.21) |
6 В е й б у л л а |
|
|
|
Р(0 = |
е х р ( — W ) , |
см . |
ф о р м у л у |
(3.51) |
Плотность распределения времени |
Интенсивность отказов |
безотказной работы / (<) |
Я (О |
1 |
Г 1 / Ы—а V I |
Г |
1 / |
Ы—аул |
е х Р |
— 2 V |
о ) J |
||
< о У 5 " Ч ~ 2 1 о ) ] ' |
|
|
|
|
см. ф о р м у л у |
(3.14) |
|
|
|
Г ( а ) |
м |
Я, ( А , / ) 0 - 1 |
exp ( — A i ) |
|
|
см . ф о р м у л у (3.17) |
Г ( а ) [ 1 - ї в ( Я - 0 ] |
|
А/уЛ"1 е х р ( - Я Л ) , см . ф о р м у л у (3.52)
Продолжение
Средняя наработка на отказ Г0
exp ( а + у )
а
Х7
>
>
§ 6.3. Учет остановок на ремонты и профилактику при оценке надежности реакторной установки и ее оборудования
Среднее время между восстановлениями. В предыдущем пара графе ничего не говорили о периодах простоя реакторной установ ки для устранения полных отказов — аварийных ремонтов (вос становлений установки) и для перегрузки горючего, плановопредупредительных ремонтов (ППР) и т. п., т. е. для предусмот ренных остановок. Предусмотренные остановки планируются за ранее, и пх длительность, как правило, известна, следовательно, она является не случайной, а постоянной детерминированной ве личиной. Обозначим долю календарного времени, затрачиваемую на остановки для производства плановых ремонтов реакторной
установки, |
б п , а |
для |
перегрузки горючего — бг : |
|
||||
|
|
|
б п = |
пп |
(0 8„//; |
бг = nv (t) |
Q'/t, |
(6.30) |
где nn(t) |
и |
nr(t) |
— число |
плановых |
остановок |
соответственно для |
||
ремонтов |
и перегрузки горючего за |
время t {t достаточно |
велико); |
|||||
бп и 0 Г —средняя длительность соответственно одного планового ремонта и одной перегрузки. Если на практике перегрузку горю
чего всякий |
раз совмещают, например, с остановкой |
для ППР, |
то при таком |
режиме эксплуатации пг(1) = 0 и бг = |
0. Продол |
жительность аварийного ремонта 0Р является случайной величи ной, так как она существенно зависит от характера отказа уста новки, который заранее неизвестен. Среднее календарное время от пуска установки после одного аварийного ремонта до другого
пуска после |
следующего |
аварийного |
ремонта |
можно представить |
|||||||
в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
Т=Т0 |
+ |
8аТ |
+ |
6ГТ |
+ |
Тв, |
|
||
или |
Т = |
(Т0 + |
|
Тв |
|
|
|
П — бР ). |
(6.31) |
||
Эту величину |
иногда называют |
средним |
временем между |
восста |
|||||||
|
)/(1 |
- |
6 |
|
|
|
|||||
новлениями. |
|
|
большой период эксплуатации t |
|
|||||||
Если за сравнительно |
просум |
||||||||||
мировать все время, затраченное на аварийные ремонты установки
0р(^)=2^рУ> т 0 Д л я этого периода среднее время восстановления
і
установки после отказа .
7 V = е р ( 0 / / я ( 0 ,
где m{t) — число полных отказов установки за время t. Обратная величина называется интенсивностью восстановления: Хв = 1 / Г в . Вычислим величину
. |
. J _ = |
i _ 6 n - 6 r = = |
l - 6 „ - 6 f |
( 6 3 2 ) |
Это среднее число восстановлений, а стало быть, и полных отка зов изделия в единицу календарного времени эксплуатации, т. е.
можно |
считать, что и, = ХК{ХК — календарный |
параметр |
потока |
полных |
отказов). |
|
|
По |
данным работы [45], полученным в результате обработки |
||
статистики отказов на 43 АЭС (со средней |
мощностью |
блока |
|
J50 Мет), проработавших в общей сложности (с 1960 по 1969 г.) около 200 лет, календарный параметр потока отказов АЭС (точнее
повреждений, исправления |
которых стоили свыше 10 тыс. долла |
ров или вызвали остановку |
АЭС на срок свыше 4 суток) составляет |
Хк 1 год-1 = 1,14 - КГ*ч_ 1 .
Вычисление коэффициентов технического использования и го товности. Если рассмотреть достаточно большой интервал эксплу атации реакторной установки t, то можно записать t = Tm(t), %{t) = T0m(t). Тогда по формуле (6.4) получаем, что коэффициент технического использования для реакторной установки за этот период эксплуатации
* т . И = |
= 1 ~ 6 Д ~ 6 Г = 1 ~ 6 Д Г 6 Г • |
( 6 - 3 3 ) |
Если период t выбрать в интервале эксплуатации, когда' остановки на плановые ремонты и перегрузку горючего отсутствуют (б п = = бг = 0), то, согласно формуле (6.6),
Я |
г |
т.п |б |
п =б |
г =о |
Т о+ Т в |
%+ |
• |
' |
|
|
= Я т . « | й |
л |
п = |
Т " |
= |
К |
v (6-34) |
||
Формулы (6.33) и (6.34) имеют общий характер, они |
справедливы |
||||||||
не только для реакторной |
установки, но и для любого ремонтиру |
||||||||
емого элемента |
ее оборудования. |
|
|
|
|
|
|||
Оптимальное |
|
число остановок .для |
ППР. Резонно |
считать, что |
|||||
средняя наработка на полный отказ Т0 является монотонно воз
растающей |
функцией доли |
времени б п , затрачиваемой на ППР* |
||
реакторной |
установки: |
|
|
|
|
Т0 = Т0 |
(бп ); |
0 < б п |
< 1 ; |
|
Т 0 (0) = const; |
Г 0 (бп ) |
(6.35) |
|
|
> Г 0 (0). |
|||
В отсутствие ППР Т0 минимальна и равна Г0 (0). При б п Ф 0 Т0(&п) возрастает, однако доля чистого времени работы в общем кален
дарном времени эксплуатации убывает и Кт.ц |
0. |
|
|
Рассмотрим случай линейной зависимости |
Т0(&п) |
— T0(Q) -f- |
|
+ 6бп , где Ь — коэффициент, |
определяемый из опыта |
эксплуата |
|
ции. Найдем то значение б п , |
при котором КТИ максимален. По- |
||
* П о д П П Р п о д р а з у м е в а е м в с е п л а н о в ы е р е м о н т ы , в к л ю ч а я и к а п и т а л ь н ы е , а т а к ж е л ю б ы е п р о ф и л а к т и ч е с к и е м е р о п р и я т и я , о с т а н о в к и р е а к т о р н о й у с т а н о в к и .
те к у щ и е
тр е б у ю щ и е
скольку d2Kr. |
< |
0, |
искомое значение б п равно |
корню |
урав |
нения й/С.т, Jd8n |
— 0: |
|
|
|
|
б°п т = j |
{ У Тв |
[Т0 |
(0) + Ть + b (1 - б Р ) ] - [ Т 0 (0) + |
Гв ]I. |
(6.36) |
Это оптимальная доля календарного времени эксплуатации реак торной установки, которую необходимо расходовать на ППР.
Видно, |
что |
ППР |
в рамках рассматриваемой модели целесообраз |
||
ны (т. е. б°п т > |
0), если |
|
|
||
|
|
|
і |
, Г , ( 0 ) |
(6.37) |
|
|
|
1 — О г |
|
|
|
|
|
|
|
|
Отсюда, |
в |
частности, следует, что |
при Тв <^ Т0(0) ППР произво |
||
дить невыгодно. Согласно формуле (6.37) это целесообразно лишь при очень большой, не достижимой на практике величине Ь. В усло
виях, когда |
Тв |
<С TQ(0), проведение ППР реакторной установки |
||
привело бы к снижению К?, я, т. е. к сокращению |
(по отношению |
|||
к случаю |
б п |
= |
0) доли чистого времени работы |
установки, что, |
например, |
для |
АЭС означало бы недовыработку |
электроэнергии. |
|
Соотношения (6.30)—(6.37) справедливы не только для реак торной установки. Они носят гораздо более общий характер. Ими
можно пользоваться при |
анализе надежности любых устройств |
и элементов оборудования |
АЭС, которые восстанавливаются (ре |
монтируются) в процессе эксплуатации. Для этого вместо доли бг
надо подставить 8 р е з |
— долю календарного времени |
эксплуатации |
|||
изделия, |
в |
течение |
которого оно находится |
в |
резерве (не |
эксплуатируется). |
|
|
|
||
Выбор |
б п |
и б р е з |
для элементов оборудования |
в |
общем случае |
должен базироваться на экономическом анализе. В конечном ито ге они должны быть такими, при которых, например, стоимость электроэнергии, вырабатываемой АЭС, минимальна. Выбор опти мальных (по экономическому критерию) долей времени для ППР и резерва — это большой и сложный вопрос, требующий отдель ного рассмотрения. Он выходит за рамки книги.
Однако, если не учитывать затраты на ППР*, то 6пПТ можно выбирать по формуле (6.36). Для этого необходимо хотя бы при
близительно знать |
коэффициент |
Ь и величину |
Г0 (0). Они |
могут |
|||||||
быть |
найдены |
опытным путем: |
|
|
|
|
|
|
|||
|
|
|
ь s |
[т0 |
(бй |
- |
г 0 |
до;;, |
|
(6.38) |
|
где |
Т0(0) — среднее |
чистое |
время |
работы |
реакторной установки |
||||||
между двумя |
полными |
отказами |
|
(аварийными остановками) при _ |
|||||||
отсутствии ППР (бп |
= 0); Т0(8„) — то же |
при б п = |
бп ^=0. |
Ины |
|||||||
ми словами, |
необходимо |
провести |
эксперимент, |
на |
первой |
фазе |
|||||
* О н и часто с у щ е с т в е н н о м е н ь ш е в ы и г р ы ш а за счет д о п о л н и т е л ь н о й в ы р а б о т к и э л е к т р о э н е р г и и .
которого эксплуатировать |
установку в течение некоторого |
перио |
||||
да времени tx без ППР, а на второй |
фазе работать в течение вре |
|||||
мени t2 уже с ППР. Длины периодов tt и t2 должны |
выбираться |
|||||
из расчета, чтобы на каждом из них произошло |
не |
менее |
5—10 |
|||
аварийных остановок. Тогда бп |
= |
6п 2 ) //2 , где 6П |
2 ) — полное вре |
|||
мя, затраченное на ППР установки за период t2. |
В свою очередь, |
|||||
Т0 (0) = V m i ; |
Т0 |
(бп) = ToJtn2, |
|
|
(6.39) |
|
где хх и mi — чистое время |
работы |
установки и число ее полных |
||||
отказов в период tx; х2 и ш2 |
— то же за период t2. . |
|
време |
|||
Если установку не удастся эксплуатировать в течение |
||||||
ни ti без ППР, то величины b и Т^О) можно оценить другим, |
более |
|||||
грубым способом. Для этого в течение времени tx |
надо эксплуати |
|||||
ровать установку с одной частотой остановок для ППР (при доле
времени бп1 5 ), а в течение |
t, |
— с другой |
|
(при бп |
2 ) ), причем |
бп 1 } ф |
|||||||||
Ф 6<2>. Тогда |
|
|
|
|
|
.Т; |
|
|
|
|
|
|
|
|
|
fii |
,)=e!i,,/fi; |
б п |
|
= 9 |
|
' д |
0 |
(бп2 > )-Г0 (б1 ( |
г |
1 ) ). |
;] |
(6.40) |
|||
|
2 ) |
п 2 |
|
2 |
|
Г 0 |
(6< , , )=т 1 |
//п 1 |
|
||||||
|
|
|
|
|
|
|
|
|
<;(2) |
• -бп1' |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Оп |
|
|
|
|
|
|
Г 0 ( 0 ) = Г 0 ( б ^ ) - 6 б » ' .
Найдя по формуле (6.36) оптимальную долю времени на ППР бпП Т , из выражения (6.30) легко получить оптимальное среднее число остановок реакторной установки для профилактических ремонтов за период эксплуатации
л 2 п т (*)=вг-*/е;. |
(6.41) |
§ 6.4. Определение характеристик надежности
по данным эксплуатации
Определение показателей надежности ремонтируемых элемен тов реакторной установки. Из предыдущих параграфов видно, что все основные показатели надежности можно найти, если из
вестны параметр |
потока |
отказов К, интенсивность восстановления |
|
Яв |
[см. формулы |
(6.15), |
(6.16), (6.33)], а также доли времени б п , |
^г, |
брез- Чтобы |
получить эти характеристики для любого ремон-" |
|
тируемого элемента реакторной установки, необходимо иметь ин формацию, приведенную в табл. 6.2, по каждой единице оборудо вания.
Т а б л и ц а 6.2
с
1
2
о
О
4
5
6
Эксплуатационные |
д а н н ы е |
д л я изделия |
|
|||
|
|
|
Е |
|
|
Причина, |
Характерные моменты |
|
|
а> |
|
д / р е з |
|
|
|
а |
Д'ППР |
характер, |
||
эксплуатации |
Дата А |
' р а б |
и прочие |
последствие |
||
|
< |
я |
|
остановки |
отказа |
|
|
|
|
|
|
|
|
П у щ е н о в |
э к с п л у а т а |
|
ц и ю |
|
|
Н а ч а т |
П П Р |
|
К о н е ц |
П П Р , |
начата эк |
с п л у а т а ц и я |
|
|
О т к а з |
(аварийная оста |
|
н о в к а ) , начат |
ремонт |
|
О к о н ч а н и е |
аварийного |
|
ремонта, |
и з д е л и е переве |
|
дено в резерв |
|
|
П у щ е н о в |
э к с п л у а т а |
|
ц и ю и т . |
д . |
|
и |
/ 2 ^1 |
— |
—. |
— |
tl |
— |
— |
' 3 — ^2 |
— |
и |
|
— |
— |
— И з н о с рабо |
чей ч а с т и ,
|
остановка |
|
реактора |
— |
на t ч |
— |
|
— |
— |
Если просуммировать за некоторый календарный период вре
мени t интервалы времени в столбцах |
4—7 (см. табл. 6.2), |
то соот |
||||||||
ветственно |
получим |
|
|
|
|
|
|
|||
|
|
т (t) = |
2 |
(ti+1 |
— ti) |
по 4-му столбцу; |
|
|||
|
|
Bp (t) |
= |
2 |
(ti+1 |
— |
t() |
по 5-му столбцу; |
|
|
|
|
0П (t) |
— |
2 |
(/; + 1 — tt) |
по 6-му столбцу; |
|
|||
|
|
врез (0 |
— ^ |
(^г-н — ti) |
по 7-му столбцу. |
|
||||
Отсюда для |
выбранного интервала |
времени t |
находим |
все ос |
||||||
новные |
характеристики надежности |
|
|
|
||||||
Т0 |
= |
т (t)/m |
(t); |
X = |
1/Т„; |
Тв = |
0Р (t)/m (t); |
6П = 0П |
(t)/t; |
|
|
|
|
|
К |
= |
УТВ; |
б р е з = |
0 р е з (t)/t, |
|
(6.42) |
гдет(^) — число отказов изделия за время t. Если нас интересуют характеристики надежности по отношению к отказам какого-то конкретного типа, то суммировать по столбцу 5 табл. 6.2 надо вы борочно, отбирая только продолжительности ремонтов после от казов заданного типа. В свою очередь, в качестве tn{t) в формулы (6.42) следует подставить число не всех отказов за время t, а толь ко рассматриваемого типа m.j(t).
Определение показателей надежности неремонтируемых эле ментов реакторной установки. Для получения характеристик на дежности таких элементов необходимо иметь следующую инфор мацию по каждому отдельному типу 'рассматриваемых изделий (табл. 6.3).
|
|
|
|
|
Т а б л и ц а |
6.3 |
|
|
Э к с п л у а т а ц и о н н ы е |
д а н н ы е д л я |
изделий типа |
А |
|
|
|
Номер |
|
Дата |
|
Время |
|
|
|
пуска |
снятия |
|
Причина, |
характер |
|||
изделия |
отказа |
замены |
н последствие |
отказа |
|||
в эксплуа |
с эксплуа - |
на новое |
|||||
|
тацию |
тацин |
|
|
|
|
|
1 |
h |
|
|
А*! |
течь, д е ф е к т |
и з г о |
|
2 |
|
|
h |
|
|||
|
|
|
т о в л е н и я , |
оста |
|||
|
|
|
|
|
|||
|
|
|
— |
|
н о в к а |
реактора |
|
3 |
(з |
' з |
|
на / суток |
|
||
п |
tn |
— |
— |
— |
|
|
|
Определенная часть неремонтируемых элементов реакторных установок отказывает в процессе эксплуатации, остальные могут быть сняты с эксплуатации по истечении некоторого периода вре мени, прежде чем наступит их отказ. Например, каналы реактора по достижении определенного выгорания извлекаются и заменя ются на свежие. Это и отображено в табл. 6.3, см. изделия № 1 и 3.
Для получения средней наработки на отказ Т0 в данном случае нельзя просто просуммировать наработки всех отказавших изде лий и поделить сумму на число отказов, так как часть изделий данного типа проработала большее время и обычно неизвестно, сколько бы эти изделия еще отработали, не будь они сняты с эк
сплуатации. Иными словами, |
информация, содержащаяся в |
табл. 6.3, не отражает надежности |
всех п изделий, она в определен |
ном смысле неполная (усеченная). Тем не менее необходимо найти показатели надежности, которые характеризовали бы все изделия данного типа.
Подобного рода задача в математической статистике называет ся сценкой неизвестных параметров по усеченным наблюдениям.
Решим |
ее |
методом наибольшего правдоподобия [см. формулу |
|
(4.15)]. |
Прежде всего надо задаться законом распределения |
f(t) |
|
случайной |
величины времени работы изделия до отказа т |
(см. |
|
табл. 6.1). |
Рассмотрим три случая. |
|
|
1. Экспоненциалный закон f(t) = (1/Г0 ) exp (—t/T0 ), где Т0 — параметр (численно равный среднему времени работы изделия до отказа), который необходимо определить по данным эксплуатации.
В результате наблюдений |
получено, что у m изделий наработки |
на отказ составили Тг, Т2, |
Тт; остальные п — т изделий про-, |
работали без отказа время |
t3, после чего были сняты с эксплуата |
ции. Полученные результаты можно рассматривать как выборку объемом п значений случайной величины т
Тъ Т2 , . . . , f m , T V C ^ V T . |
(6.43) |
n — m |
- |
