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

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

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

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

Добавлен: 11.04.2024

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

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

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

Следовательно, искомый средний зазор кладки по формуле (2.10)

оооо

 

Д = Л Ї ( Д ) = J д

/ ( A ) d A =

<J (A al^ix) х

 

 

x e x p f

[ --Y]dH=*2alY

п= 1,13a.

(12.3)

 

L

4 V о / J

 

 

Если известно максимально возможное отклонение размера

кирпича

от номинала

б х м а к о ,

то, согласно условию (3.12),

можно

принять

a = 6 х м а к с / 3 .

Тогда

 

 

 

А = 2 б х м а к с / 3 У л =

0,38бх м а к с .

 

Зная закон (12.2), можно определить и максимально возможную величину зазора кирпичной кладки А м а к 0 , например такую, что вероятность получить больший, чем А м а к с , зазор равна 0,001. По фор­ муле (2.8) вероятность, дополняющая упомянутую до единицы, равна

 

 

мани

мсти

л

X

 

 

 

 

 

X ехр

. _ L j ^ A ^ 2 " ^ д= 2

ф ^ - ^ | с ) = 0,999.

Из этого уравнения

находим Ф (а) = 0,4995

 

или и = 3,29

(см. табл. П.1). Тогда

 

 

 

 

А м а к с =

о-]/2

и = 3,29аУ2

= 4,64а =

1,55бх м а к с .

§ 12.2: Оценка вероятности касания элементов системы «труба в трубе»

Подобные системы получили большое распространение в реакто­ ростроении. Примерами их могут служить рабочий канал активной зоны в цилиндрической ячейке графитовой кладки реактора, кольце­ вая сборка твэлов, элементы парогенераторов и др. Обозначим вну­ тренний диаметрнаружной трубы dHap, внешний диаметр внутрен­ ней трубы dBH, а эксцентриситет их осей є (см. рис. 7). Минималь­ ный зазор между трубами

Д = min -J- ( d H a p — О = - і - (duap—d) — є.

(12.4)

(поф) z

z

 

Величины duaP, dBK и e являются случайными, поскольку на ар> ^вн накладываются погрешности изготовления, а є полностью обусловлен погрешностями сборки. В качестве законов их распре-


деления с достаточной точностью можно принять соответственно нормальный закон (3.9) и закон эксцентриситета (3.29):

Jnap w e

 

 

"нар —"„ар

 

x

p

-•нар

(12.5)

Жвп) = о-вп У 2л.

 

1

» \21

 

ЄХР

 

 

7 =

 

 

 

 

2 \

а в н

 

/ 2 (є) = — j -

є ехр

 

(12.6)

2 I 0-Е

где d"ap, с?вн — номинальные значения диаметров; их можно рас­ сматривать как математические ожидания:

 

 

М (d n a p ) = d"ap,

М (dB H ) =

d"„.

 

 

Чтобы оценить вероятность касания рассматриваемых труб,

необходимо знать закон распределения случайной

величины А.

Это не сложно сделать, так как закон распределения

случайной ве­

личины z = (dnav

— dBH)!2

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

законов, т. е. тоже нормальным законом:

 

 

 

 

 

/ і

(2) =

 

ехр

 

 

 

 

(12.7)

 

 

 

, Т / 2 я

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

° z =

V D

("2"d "a v )+ D

і ~

\ d ^

=

Т ^ 0 " а р + о | , , :

 

гп

= М {y(dHap)-f-

 

7W ( -

у dB 1 I )

=

- у № Р - ^ ™ ) -

І 1 2 - 8 )

• Следовательно, искомый закон / (А) — это закон распределения разности А = z — є двух случайных величин, одна из которых подчиняется нормальному закону (12.11), а другая — закону экс­ центриситета (12.9). По формуле (2.31) интегральный закон распре- - деления случайной величины А имеет вид

F (А) = Р {г—є < А} — Р {г < є + А} =

оо Г в + А

0 Є ,

О 1 - 0

где S — область координатной плоскости (z, е), лежащая ниже пря­ мой z = є + А; А — фиксированное значение рассматриваемой


случайной

величины;

є ^

0. Тогда искомый

дифференциальный

закон, согласно формуле

(2.3), примет вид

 

 

 

 

 

 

/ ( Д ) = / " ( Л ) = 5 M e + A)/a (e)de =

 

 

 

 

 

exp

 

8+ Д-

•є exp

 

 

 

de =

.)

az-\/2n

2

 

 

 

 

 

о

 

 

 

 

 

Оє

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о2 "|72я

 

2

 

1

_ »ї

2

\

a

a z

X

 

 

 

 

о

) .

 

 

X (zH —A) exp

 

- Д \2

 

 

 

(12.9)

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

а 2 = а |

+ о| .

 

 

 

 

 

 

 

 

 

Математически

(формально) условие

касания

можно

записать

в виде А ^

0. Поэтому искомая вероятность касания труб по фор-

'муле

(2.8) равна

 

 

 

 

 

 

 

 

 

 

 

р=

\ /(Д)4Д = 0 , 5 - ф ( - ^ - )

+

 

 

 

 

 

exp

 

2

2

0,5+

Ф

 

 

 

(12.10)

 

 

 

 

CFz + OV

 

 

 

 

 

 

Применение дистанционирующих

устройств для

предотвращения

касания элементов системы «труба в трубе» приводит к снижению

только

среднего

квадратического

отклонения для эксцентриситета

о-8.

Из формулы

(12.10)

видно,

что при а 8

0 р

уменьшается,

стремясь

к р =

0,5—Ф (2H /az ).

Однако, поскольку сами дистан-

ционирующие

устройства

имеют

погрешности

изготовления,

а 8

хотя

и снижается, но практически всегда остается

отличной

от

нуля. Величины

az и оЕ можно оценить по формулам:

 

 

 

 

 

 

az

= (1/2) 1/онор + о1нМ1/6)Кб|,р + б!„,

[(12.11).

где би а г)

и

б в н

— половины допусков

соответственно

ДЛЯ fifjнар

и dBa

[использованы формулы (12.8) и

(3.12)],

и

 

 

 

где е ы

 

•— максимально

к°£ ^

е м а к с / 3 , 7 ,

 

 

(12.12)

а к с

возможный эксцентриситет для конкрет-

 

 

 

 

 

 

ВОЗМОЖНЫЙ ЭКСЦ1

 

 

 

 

ной системы [вероятность превысить его равна 0,001].


Раздел V. ИСПОЛЬЗОВАНИЕ СОВРЕМЕННЫХ СТАТИСТИЧЕСКИХ МЕТОДОВ ДЛЯ ПРИНЯТИЯ РЕШЕНИЙ В УСЛОВИЯХ РЕАКТОРОСТРОЕНИЯ

Г л а в а 13.

В Ы Б О Р Р Е Ш Е Н И Я НА О С Н О В Е Э К С П Е Р И М Е Н Т А Л Ь Н Ы Х Д А Н Н Ы Х ( Э Л Е М Е Н Т Ы Т Е О Р И И С Т А Т И С Т И Ч Е С К И Х Р Е Ш Е Н И Й )

§ 13.1. Выбор решения при помощи проверки простой гипотезы против простой альтернативы

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

сотсутствием необходимых опытных данных вообще.

В§ 4.3 рассмотрена проверка статистических гипотез без учета потерь, к которым приводят неправильные решения. В классиче­ ской статистике потери учитываются лишь косвенно при помощи таких понятий, как уровень значимости, мощность критерия, дове­ рительные вероятности и т. п. Сравнительно недавно разработана теория статистических решений, в которой делается попытка учесть конкретные экономические, психологические и т. п. убытки (выигры­

ши) [115]. Ниже рассмотрены

некоторые положения этой теории

и даны примеры использования

ее в задачах реакторостроения.

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

кое состояние из

S0 , Slt

Sm

имеет

место. Предположив, что

объект находится

в состоянии

S.;,

можно

проверить эту гипотезу,

пользуясь опытными данными. В качестве исходного материала для такой проверки необходимо иметь ряд возможных состояний объек­

та S0 ,

Sj,

Sm, набор соответствующих им

решений уа, у

ут

-относительно

истинности этих состояний и

выборку с элементами

хг, х2,

хп,

представляющую собой п повторных эксперименталь­

ных замеров параметра объекта, по величине которого можно судить

о состоянии объекта Sj. Ситуация наиболее проста при двух

взаимо­

исключающих состояниях 50

и Si.

В этих условиях можно выдвинуть

две гипотезы: Я 0

— о том,

что

объект находится в состоянии S0

и ее альтернативу

Я 3 —• объект

находится в состоянии Sx.

Таким

образом, должна быть проверена простая гипотеза против простой альтернативы (см. § 4.3).


Проверка проводится по следующему правилу

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

(4.46) ]. Если конкретная выборка 2 = и

х2<,

хп)

попадает в об­

ласть приемки G0

= 0 1 ф , то принимается

решение у0

о справедливо­

сти гипотезы Н0

(объект находится в состоянии S0):

х2

хп)

€ Go^-yo-

Если

же точка

г = г,

х2,

 

хп)

 

попадает в критиче­

скую область Gx

=

0,.р , то принимается решение ух

о справедливости

гипотезы

Нх:

( х ь

х2,

 

 

хп)

£G1-+yl.

 

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

При

такой проверке, как уже отмечалось

(4.47)], возможны ошибки первого и второго

рода. Обозначим П 1 0

потери,

например экономические, при ошибке первого-рода (примем

гипотезу Hlt

когда справедлива Я 0 ), а П 0 і

при ошибке второго рода

(принимаем гипотезу Н0,

когда справедлива

Нг).

Обозначая далее

П 0 0

и

П п

выигрыши

(«отрица­

 

 

 

 

 

Т а б л и ц а

13.1

тельные

потери»

 

П 0 0

< О,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П . і < 0 ) ,

получаемые

при

пра­

 

 

П л а т е ж н а я

м а т р и ц а

 

вильных решениях, можно запи­

 

 

 

 

 

 

Гипотезы

 

сать

так

называемую

таблицу

Решения

 

 

я„

 

 

я ,

потерь

(или

платежную

матри­

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

цу)

(табл.

13.1).

 

 

в

том,

 

То

 

 

П 0 0

 

 

Поі

Задача

заключается

 

 

 

 

 

 

 

 

 

 

чтобы найти способ выбора об­

 

Ті

 

 

Пю

 

 

 

ластей G0

и Glt учитывающий ве­

 

 

 

 

 

 

 

 

 

личины этих

потерь. При выборе областей О п р

и 0 к р

в § 4.3 вопрос

о потерях

полностью

игнорировался.

Если

найдем

указанный

выше способ, то

получим

искомое

правило

для

выбора

решения.

Бейесовское правило выбора. Предположим, имеются состояния

Sj> І =

0> 1» •••> т

1 1

соответствующие им решения yh,

k — 0, 1,

т.

Если

наряду

с

таблицей

потерь

Ukj

известны

и априорные ве­

роятности

состояний pj

=

P{Sj),

то для

выбора

решения

можно

воспользоваться бейесовским правилом, которое основывается на

минимизации так называемой

средней функции

риска

R. Послед­

няя

представляет собой

среднее

по

всем возможным

состояниям

Sj

значение (математическое ожидание)

риска:

 

 

 

 

т

 

 

 

m

m

 

 

 

 

R = 2

Р]Г}= 2

 

РІ 2

n f t i

х

 

 

 

і = 0

/ = 0

 

= 0

 

 

(13.1)

 

хР{(хг2,

...,хп)

 

Є

Gk/Sj},

 

 

 

где Г} — условный риск для состояния Sj

(математическое ожидание

потерь); Р {(хъ х2,

xn)£Gk/Sj]

 

—условная

вероятность по­

падания выборки ъ х2,

 

хп)

в область Gh,

если в действитель­

ности имеет место состояние

Sj.

 

 

 

 

 

 

 

Для двух взаимоисключающих событий средняя функция риска

(13.1)

 

 

 

 

 

 

 

 

 

 

 

R =

р < / о

+

р л -

 

 

(13.2)