Файл: Клемин А.И. Инженерные вероятностные расчеты при проектировании ядерных реакторов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 236

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Считаем, что распределение каналов по величине отклонения кон­ струкционного параметра х£ от номинального значения подчиняется нормальному закону, симметричному относительно середины поля допуска [х" — Дг ; х" -+- Д; ]. В этом случае можно записать (ин­ декс і для простоты опустим)

где N — полное число технологических каналов в активной зоне, da — количество технологических каналов с отклонениями Ах, ле­ жащими между а и а + da; о — среднее квадратическое отклонение конструкционного параметра.

Из общих геометрических соотношений для активной зоны и из условия б) монотонности изменения конструкционного параметра от центра к периферии можно найти связь между г и a, dr и da. Дей­ ствительно, число каналов активной зоны, расположенных между mr~\~dr,

dn = 2nrdr/S

(5.71)

где S m — площадь элементарной ячейки реактора, т. е. часть попе­ речного сечения активной зоны, приходящаяся на один канал. При­ равняем величины dn, выраженные формулами (5.70) и (5.71). Это и будет означать выполнение условия б), если считать, что при изме­ нении г от 0 до R величина а = х н — х меняется от —А до -f-A. Имеем

или

(5.72)

С другой стороны,

 

яг

 

R

Д

nR2 = Г 2nrdr = Sm

С Гехр

откуда

 

а

 

(5.73)

Подставляя найденное выражение для rdr (5.72) в формулу (5.69), получаем

 

 

 

 

 

 

С

 

 

 

 

Г

\

f а

\ 21

\ d a

 

 

 

 

 

 

 

J

Ч» ['(«)] <* exp

-J[—)

 

 

 

 

 

Ap=2^=4

 

 

 

 

 

J _

 

:

 

-•

(5-74)

 

 

Vі

[r (a)] exp

 

/ a _ \2 -

 

 

 

 

 

 

 

J

 

2

[ a

 

da

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрируя

по частям в числителе,

получаем

 

 

 

J W [r (a)] a exp

~ ^

[ ^ ) 2

] d

a =

— ^ ( a ) a 2 e x p

2 V cr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a2 Г

Ч ^ . - ^ е х р

 

J _

/ a_

 

da.

(5.75)

 

 

 

+

 

2 1

a

 

 

 

 

 

 

J

 

dr

da

 

 

 

 

 

 

 

 

 

 

—д

 

 

 

 

 

 

 

 

 

 

 

 

Из выражения

(5.73)

находим

 

 

 

 

 

 

 

 

 

 

 

dr

 

 

 

 

exp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

|l/2

 

 

 

da

 

20"1/2я

J 2 0 ^ j K t ) + « ( t ) ]

 

 

 

 

 

 

 

 

Поскольку обычно для конструкционных параметров

активной

зоны

Д / а ^ З ,

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф (Л/а) « 0 , 5

 

 

 

 

(5.76)

[см. табл. П.1]. Тогда выражение

для dr/da

упрощается и

j > [ r ( a ) ] a

exp

1

(SL ^'

da=

 

a2

exp

I f

A

 

—д

 

 

 

 

 

2 I a

 

 

 

 

 

 

2 I a

 

 

 

 

 

 

 

 

 

 

 

 

 

exp Г-

/ «

~V" da

 

^

(

«

я

+

^

 

? l r ( . ) ] f x -

.

V a J .

 

 

 

 

 

/

a V 1/2

 

 

 

 

 

 

 

— Д

 

 

 

 

0.5Ч

\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о

 

Согласно выражению (5.73), r(a=

 

—A) = 0 и r (a =

+ A ) = i?;

A/a ^

3,

т.

e.

exp

 

1 / Д

\ 2 <^

1,

поэтому первое

слагаемое

 

 

 

 

 

 

 

 

or

 

 

 

 

 

 

 

 

 

существенно меньше второго, и им можно пренебречь. Оставшийся интеграл можно вычислить, имея в виду, что exp [—(a/a)2] являет­ ся наиболее резко меняющейся функцией в подынтегральном выра­ жении. Для оценки величины Ар допустимо слабо меняющиеся функ-



цим вынести из-под знака интеграла, взяв их в точке максимума функции ехр [—(а/0)2 ], т. е. при а = 0. В итоге находим

| ¥ -

[/• (а)] а ехр

 

 

 

da :

 

 

\V2)

dr

r = R/)R2

—д

 

 

 

 

 

 

 

 

 

так как

г (0)

R/V2.

 

 

 

 

 

 

 

 

 

 

 

Аналогичным образом вычисляем интеграл в знаменателе фор­

мулы (5.74):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

' Ч'2

[г (а)] ехр

 

 

 

а_

\ 2

 

 

 

 

 

 

 

 

 

а

da:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y2 [/-(0)] j

ехр

 

 

 

da « У2noW

[

^

 

 

 

—д

 

 

 

 

 

 

 

 

 

 

 

Воспользовавшись полученными выражениями, формулу (5.74)

можно преобразовать к виду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d4

 

 

 

 

Ар = 2 / г

\ І

і

dr

r = R/Y2

 

R

dr

r=R/Y2

Vі 4 rr

(^i 77\

 

]г = к/ r

 

 

 

'\

"* \г=к/у 'і

В частном

случае,

когда W (г) =

J0

 

2,405

г),

где

#э —-эффек-

тивныи

радиус

активной

зоны,

получаем

 

 

 

 

 

 

 

 

 

j

(

2 ' 4

0 5

 

 

 

 

 

 

 

Ар

R

 

2,405

 

у 2

-

 

 

 

 

 

(5.78)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

УТл

 

Rs

(

2,405

 

R

 

 

 

 

 

 

 

 

 

 

 

У2

'

Rs

 

 

 

 

Здесь J0

и /х —функции

Бесселя

 

1-го рода соответственно нулевого

и первого порядков.

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом,

формулы (5.77)

и

(5.78)

устанавливают

связь

между характеристиками

рассеяния

конструкционных

параметров

(в пределах допусков) и максимальным отклонением запаса реактив­ ности, возможным при специальном размещении технологических каналов в активной зоне. Средние квадратические отклонения а г можно вычислить по результатам замеров конструкционных пара­ метров xt на определенной партии каналов, см. гл. П . Если партия еще не изготовлена или недостаточно исследована, то можно оценить величины ot приближенно, положив at = Aj/З или же (это будет в запас) приняв ot = Аг / |/2я , как для треугольного распределения Симпсона [6].


Г л а в а

6.

О Ц Е Н К А Н А Д Е Ж Н О С Т И Р Е А К Т О Р Н О Й У С Т А Н О В К И *

§ 6.1.

Основные количественные показатели

общей

надежности реакторной установки

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

Основными понятиями современного количественного подхода к проблеме надежности ядерных реакторных установок являются следующие.

Н а д е ж н о с т ь р е а к т о р н о й у с т а н о в к и (в широ­ ком смысле) есть свойство, обусловленное ее безотказностью, ремон­ топригодностью и долговечностью и обеспечивающее нормальное выполнение установкой требуемой задачи в заданном объеме и в за­ данных условиях эксплуатации.

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

В отличие от отказа под н е и с п р а в н о с т ь ю

будем пони­

мать такой дефект оборудования, который не требует

немедленного

прекращения его эксплуатации, т. е. с которым оборудование в тече-

*

П о д

р е а к т о р н о й у с т а н о в к о й

д л я о п р е д е л е н н о с т и т а м , где это не

огово­

р е н о о с о б о ,

б у д е м п о н и м а т ь л и б о

р е а к т о р н у ю п а р о п р о и з в о д и т е л ь н у ю

уста­

н о в к у , л и б о о т д е л ь н ы й б л о к А Э С , л и б о м о н о б л о ч н у ю А Э С в ц е л о м .

 

4 Зак.

1282

 

81


ниє определенного времени еще может выполнять свои основные функции.'

Б е з о т к а з н о с т ь р е а к т о р н о й у с т а н о в к и (или надежность в узком смысле) количественно характеризуется вероят­ ностью Р (і) того, что в заданных режиме и условиях эксплуатации установка проработает без отказов в течение времени і:

 

Р ( 0 = Р { т > 0 .

 

(6-і)

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

реакторной

установки, или нара­

ботка на отказ, — случайная величина.

 

 

Р е м о н т о п р и г о д н о с т ь р е а к т о р н о й у с т а н о в -

к и

количественно характеризуется суммарным

временем аварий­

ных

ремонтов 9р и суммарным временем планово-предупредитель­

ных

ремонтов Ґ1ПР (профилактических работ) 6П

за весь период ее

эксплуатации

 

 

 

 

8р = 2 6р ;,

6П = 2 0D ;.

 

/

/

 

 

Д о л г о в е ч н о с т ь р е а к т о р н о й

у с т а н о в к и ко­

личественно характеризуется ее сроком службы т0 . с

и полным тех­

ническим ресурсом

 

т т . р = X Xj,

(6.2)

т. е. полным временем безотказной работы реакторной установки до момента вывода ее из эксплуатации (т; — интервалы работоспособ­ ности). Срок службы

т с . о = Т т . р + Є = 2 (т, + BPJ- + 0 n J + Q'j)

(6.3)

і

 

есть полная календарная продолжительность эксплуатации реактор­

ной установки (6 — полное

время простоя установки,

включающее

в себя периоды аварийного

ремонта — 0р ;, ППР —

0п ; - и прочие

остановки — Q'j, например, для перегрузки ядерного топлива и т. д.). Срок службы стационарных энергетических реакторных установок обычно составляет 20—30 лет [41—43].

Количественные показатели общей надежности реакторной уста­ новки за период времени t следующие.

1. Коэффициент технического использования (КТИ), называемый иногда реакторостроителями коэффициентом использования кален­ дарного времени (КИКВ), равен отношению «чистого» времени рабо­

ты т (І) установки за некоторый

календарный период эксплуатации

к величине этого периода t:

 

 

/Ст. и О = т {t)lt

= х (0ЛЄ (0 + т (01,

(6.4)

где G (t) — полное время простоя установки за период t. Для установки, снятой с эксплуатации, Кт.,, = тт > р с , с . Для доста­ точно большого периода t КТИ численно равен вероятности застать