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

книги из ГПНТБ / Клемин А.И. Инженерные вероятностные расчеты при проектировании ядерных реакторов

.pdf
Скачиваний:
34
Добавлен:
25.10.2023
Размер:
14.61 Mб
Скачать

и при экспоненциальном законе надежности также равна парамет­ ру закона к.

Выясним физический смысл интенсивности отказов на примере экспоненциального закона (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 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

-

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