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

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

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

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

Добавлен: 11.04.2024

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

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

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

В действительности р и q не являются детерминированными вели­ чинами, равными р яд, а распределены около этих средних значений по вероятностному закону. Как показывают непосредственные ис­ следования, закон распределения р (характеристик прочности ма­ териалов) в абсолютном большинстве практических случаев можно описать одним из следующих четырех законов: нормальным (3.9), логарифмически нормальным (3.14), Вейбулла (3.52) и гамма-рас­ пределением (3.17). Что касается нагрузок, а стало быть, и напря­ жений q, то в литературе пока очень мало данных о их вероятност­ ных законах распределения. Интуиция инженера и теоретические предпосылки говорят о том, то здесь следует ожидать закон, близ­ кий к нормальному. В некоторых случаях, по-видимому, возможны и исключения. Например, в работе [15] принят (правда, бездоказа­ тельно) логарифмически нормальный закон для нагрузки.

В упомянутой работе рассматриваемая задача решена в пред­ положении, что и напряжение q, и характеристики прочности мате­ риала р распределены одинаково логарифмически нормально. Бо­ лее реальным является случай, когда q подчиняется нормальному закону, а р одному из вышеперечисленных, допустим, распределе­ нию Вейбулла. Именно для этого случая и попытаемся вычислить искомую вероятность Fa (t). Согласно формулам (3.9) и (3.52), плот­ ности распределения q и р имеют вид:

'1

1 i q - q

;

/ 2 ( p ) = ^ . Y P Y - i

 

 

ехр(-ХрУ),

(6.103)

где ад — среднее квадратическое отклонение напряжения от сред­ него значения q. Для большей наглядности результатов предполо­ жим, что у = 2 (закон Релея). Тогда / 2 (р) = 2Хр ехр (—Хрй) и, согласно формуле (3.53),

p = M(p)=(l/2)Yn/X

или ^=(я/4)(1/р 2 ) .

(6.104)

Условием отказа рассматриваемого элемента конструкции является неравенство

x = p—q<0.

(6.105)

Поскольку q и р — случайные величины, х — тоже случайная ве­ личина; найдем ее закон распределения. Для этого воспользуемся


формулой (2.31). Предполагая, что q и р — независимые величины т. е. / (<7, р) = h (<?)7а (р), получаем

 

 

 

 

 

со

 

со

 

 

 

F(x)=P{p-q<x}=P{q>p-x}=§[

 

 

 

 

J

 

Ш Ы р ) ^ ] dp\

 

 

 

 

 

О P—x

 

 

f ( X ) =

F> (x) =

f h (p -

x) U (P) dp =

 

0 g V 2 : r t

X

 

 

 

 

 

 

 

 

 

 

X f pexp( —2 )-ехр

 

1 (

P

-

x

- q ^

 

 

 

 

 

 

 

 

 

 

2Xo0

1

 

±

/ \

£

±

!

 

 

 

1-^2X02

'-|/2n

 

 

 

 

 

2

V

 

a „

 

 

 

2A

 

 

 

(x +

 

q)exp

 

 

 

 

 

 

 

( l + 2 X o = ) 3 / 2

2

^ с г д ] Л і + 2 Я а

2

 

 

 

 

+

l + 2 A a 2

9/ J

 

 

 

 

(6.106)

Искомая вероятность отказа конструктивного элемента за время

Интегрируя функцию (6.106)

[в частности, второе

слагаемое—по

гіФ . ,

ехр( — 0 , 5 u 2 )

dul

 

частям, помня, что — (и) = — ' . — -

— ,

получаем

dx к '

У2к

dxj

J

' • » - Ї + Ф ( £ ) - 7 І Т Я Е ; 2

X

Хехр

(6.107)

1 + 2 Ь =

 

Время t в правую часть выражения (6.107) входит неявно через % [см. формулу .(6.104)], так как р — р {t). Выражение (6.107) су­ щественно упрощается, если учесть, что обычно aq ^ Q,lq.

Отсюда

l + 2 W < l + 2

. = ! + _ i L . _ L

100

|200 / с |

Поскольку всегда

К3

> 1,

то с достаточной точностью

можно

считать 1 + 2 Х о | «

1.

Таким

образом, O ( g / a 5 ) « O ( 1 0 ) =

1/2 (см!


табл. П.1), аналогично

и Ф .

.

: — .

Окончательно

 

 

 

 

 

 

2

 

 

 

формулу (6.107) можно

записать в

виде

 

 

 

 

F n ( 0 =

l - e x p [ - ^ ] = l

•ехр

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1 —ехр

п

 

 

(6.108)

 

 

 

 

 

 

 

 

Из выражения

(6.102) следует, чго К3 p[t)/q = K3(t).

Следова­

тельно,

надежность

конструкции

R (t)

[см. выражение

(6.101)],

точнее,

вероятность,

что она проработает безотказно

в

течение

времени t при

принятых законах

распределения

q

и

р и при

 

 

 

R(t)=exp

4/СІ

(t)

 

 

(6.109)

 

 

 

 

 

 

 

 

Например, при К3 (0 =

3/? ( 0 = 0,92; при К3 ( 0 = 2

R it) = 0,82;

при К3

(0 = 1

R (0 =

0,45.

 

 

 

 

 

Воспользовавшись формулой (6.109), легко получить интенсив­

ность

отказов

для

рассматриваемого

конструктивного

элемента

4 0 =

R'(t)

1

dK3

 

Rit)

 

dt

2K3(t)

dt

 

 

K

'

(Здесь, как и выше,

предполагаем,

ято-q

не зависит

от времени.)

Таким образом, зная зависимость р

(t)

или изменение

во времени

коэффициента запаса

К3 (0 = Р (0/<7> п

о

формулам (6.109) и (6.110)

несложно рассчитать основные характеристики надежности эле­ мента конструкции. Некоторые фактические данные по прочности реакторных сталей, необходимые для оценки R (t), можно найти в работе [15].

Расчет надежности трубопроводов и арматуры. Наиболее прос­ той и рациональный подход к расчету надежности трубопроводов и арматуры реакторных установок базируется на статистике отказов соответствующих элементов в эксплуатации. Для арматуры метод расчета полностью идентичен изложенному в § 6.2 и 6.5. Достаточно знать интенсивность отказов Хг для каждого устройства (вентиля, клапана, задвижки и т. д.), тогда вероятность безотказной работы отдельного устройства в течение времени t можно вычислить по формуле (6.16). Величины Хг можно получить путем постановки специальных испытаний элементов арматуры на надежность либо в результате сбора данных об их отказах в эксплуатации (см. § 6.4), а также из литературных источников (см. табл. П. 10 в приложении).

Расчет надежности трубопроводов тоже можно вести по формуле (6.16). Интенсивность отказов в этом случае

%1 Xi дог / j "Т~ Xi г и | з / V r n 6 -4- Xf ш її •

(6.111)


где Я І П О Г параметр потока отказов одного погонного метра пря­ мого участка г-го трубопровода длиной /»; кігпб и кіт — параметры потоков отказов отдельного гиба н одного погонного метра,сварного• шва t'-ro трубопровода; NiTa6 и — число гибов и длина сварных швов на 1-м трубопроводе. Обычно доминирующим слагаемым в пра­ вой части выражения (6.111) является последнее.'Для системы тру­ бопроводов, например АЭС в целом, интенсивность отказов

 

т

 

л

= 2 Х

(6.П2)

где т — полное число ниток

различных

по типу трубопроводов

на АЭС. Зная Л, по формуле (6.16) легко вычислить вероятность, что -за время t не произойдет ни одного отказа системы: Р (t)=exp (—At). Среднее число отказов такой системы за время t, согласно формуле

(6.10), составит т (t) = At. Исходные данные ( і п о г ; Хгяб; кш) для расчета надежности трубопроводов разных типов правильнее всего

определять по результатам специальных испытаний или на основе опыта эксплуатации, также как для элементов арматуры. Некото­

рые данные по этим к

можно найти в таблице П. 10 в приложении.

При оценочных, прикидочных расчетах

надежности допустимо для

системы

всех трубопроводов АЭС (D >

100 мм) задать одинаковые

значения

к с запасом, например

 

 

К«* = Ко? + U

- ~ - + К у =

Ю8 -МО» 4 - і • м-\

(6.113)

Когда никаких (даже ориентировочных) статистических данных об отказах конкретного трубопровода, используемого на АЭС, нет (в си­ лу его уникальных размеров, форм, материала, условий работы и т. п.), оценку надежности можно провести по методике, аналогич­ ной изложенной в начале параграфа, рассматривая только сварные соединения.

§ 6.8. Особенности оценки теплотехнической

надежности реактора

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


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

Однако часто

статистических

данных об отказах

либо вообще

нет (проектируется уникальная

установка), либо

недостаточно.

В этих условиях

прогноз надежности, учитывающий

все факторы

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

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

в том числе уникальных

реакторов на этапах их проектирования

и в период эксплуатации.

Расчет теплотехнической надежности R

позволяет, во-первых, оценить уровень общей надежности реактора і ? о б щ , в частности, найти для него верхнюю границу, так как всегда •^общ ^ Я; во-вторых, оптимизировать параметры активной зоны по критерию теплотехнической надежности (пример такой оптими-' зации приводится в работе [61]).

Формула для расчета показателя R теплотехнической надежно­ сти. Показатель R представляет собой вероятность, что за время t работы реакторной установки в стационарном режиме не произойдет ни одного отказа активной зоны реактора по теплотехническим при­ чинам,.т. ё. пережога, прогорания, расплавления, связанного с на­ рушением нормального теплоотвода от твэлов. Такую вероятность R можно найти, используя подход, изложенный в § 5.3 [см. формулу 5.49)] и работах [6, 61—64]:

Я = WKh=

П [0,5 + Ф(т]Н/аА )]П А .

(6.114)

А=1

й=1

 

Поскольку обычно величины R и тем более Rh (теплотехническая надежность отдельного канала) близки к единице (иначе реактор будет ненадежен в теплотехническом отношении), формулу (6.114)