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

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

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

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

Добавлен: 11.04.2024

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

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

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

 

Гипотеза о равенстве

дисперсий

двух

нормально

распределен­

ных

величин. На

основе

выборок (4.52)

вычислены

эмпирические

дисперсии

0и- и а\у. Требуется

проверить гипотезу Я 0 о равенстве

генеральных дисперсий: ох =

а£. Вычисляем

к р и т е р и й

Ф и-

ш е р a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЇЇ =

оЪ/ofy.

 

 

 

 

 

(4.55)

 

Если

вычисленное

по

результатам

наблюдений

значение §

попадает

в область

х

р ^

f

^ хр, то

гипотеза

Н0

принимается

на

 

 

 

1

_ 2

 

 

2

 

 

 

[34]. Величины

хр/2

уровне значимости

р\ иначе — отвергается

и

х^

р находятся по табл. П.5 при следующих

аргументах:

 

 

 

 

 

 

 

( у !

>Ч — \\

/?3 —1);

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

=

 

(4.56)

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- 1 ;

пг-~ 1

 

 

 

 

Гипотеза о корреляции двух величин. Инженера часто интере­ сует вопрос: значима ли корреляционная связь между рассматри­ ваемыми величинами X и Y, или ею можно пренебречь? Для про­ верки этой гипотезы необходимо найти величину

 

 

 

 

 

z=-Linl+P

(4.57)

 

 

 

 

 

 

2

1 - р

 

где р — эмпирический

коэффициент

корреляции, вычисляемый по

формуле (4.11) на основе результатов наблюдений

над величина­

ми X

и Y: х

у

х2, у2,

хп,

уп

(п — объем выборки по каждой

величине). Если | Z | ^

u0 ,s|(i-p)/]/n—3, то гипотеза

об отсутствии

корреляции

(р =

0) между X

и Y принимается, иначе — отвер­

гается

на

уровне

значимости

(3. Значение Uo,5(i-p)

определяется

по табл.

П.1.

 

 

 

 

 


Раздел I I . РЕШЕНИЕ ВЕРОЯТНОСТНЫХ ЗАДАЧ, СТОЯЩИХ ПЕРЕД КОНСТРУКТОРОМ И ИНЖЕНЕРОМ-РАСЧЕТЧИКОМ ЯДЕРНОГО РЕАКТОРА

Г л а в а 5.

О Ц Е Н К А Т О Ч Н О С Т И И Н Ж Е Н Е Р Н Ы Х Р А С Ч Е Т О В

ЯД Е Р Н Ы Х Р Е А К Т О Р О В

§5.1. Постановка задачи

Точность любого инженерного расчета количественно характе­ ризуется погрешностями результатов расчета. Естественно считать, что точность не может быть больше единицы или 100%, поэтому связь между точностью t и погрешностью А записывают в виде

/ = 1 _ | Д | или t% = 100 | А°/о |.

(5.1)

Результатом инженерного расчета может быть одна величина г или целая группа (непрерывная функция или поле) величин Ги Задача оценки точности любого инженерного расчета может быть решена, если известна процедура вычисления погрешности отдель­ ного результата расчета

Лг = ки — / ' ф ,

(5.2)

где г„ — результат расчета (номинальное

значение); г ф — факти­

ческое значение. Основные факторы, влияющие на точность инже­ нерного расчета, следующие:

1) адэкватность принятой математической модели реальному явлению (глубина знания природы, физики конкретного явления или процесса в конечном итоге проявляется в выборе определенной методики расчета, которая всегда является каким-то приближением

кдействительности);

2)погрешности исходных данных (отклонение фактических ве­ личин от принятых в качестве исходных для расчета, связанное, например, с наличием допусков на параметры и т.п.);

3)погрешности эмпирических соотношений, формул и констант,

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

4)погрешности вычислительных операций, численных матема­ тических методов, округлений ит . д.;

5)грубые ошибки, описки, промахи.

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

3*

51


соотношениям; в-третьих, путем проведения повторных расчетов, тщательного анализа промежуточных результатов и всего хода вы­ числения.

При инженерных расчетах обычно приходится задаваться до­ пустимыми погрешностями вычислительных операций е в ы ч . Надо придерживаться следующего простого правила: погрешность вы­ числения должна быть приблизительно на порядок меньше мини­ мальной из погрешностей исходных данных. Использование элек­ тронных цифровых вычислительных машин (ЭЦВМ) позволяет вы­ полнять вычисления с достаточно высокой точностью.

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

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

результатов

расчета.

Оставшиеся две группы факторов (см. п.2 и 3) имеют случайный

характер.

Действительно,

практически никогда не известно зна­

чение фактической погрешности Ах любой

величины х, задаваемой

в исходных

данных или определяемой

по

эмпирической формуле

для конкретных условий рассматриваемой инженерной задачи. Обычно задают номинальное (при котором проводится расчет) зна­

чение х н

и максимально возможную погрешность А (0 <

| Ах | ^ Д),

которая,

например, в случае геометрического размера

является

половиной допуска для этого размера, в случае эмпирической фор­

мулы — максимальной

ошибкой результата. Таким образом, фак­

тическое значение х

заключено внутри интервала х н — Д <! х <1

< ха + Д, т.е. х типичная непрерывная случайная величина,

распределенная 'б области хн — Д

х

<

хн

+

А по какому-то ве­

роятностному закону / (х).

 

 

 

 

 

 

 

 

 

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

погрешностей пп. 2 и 3

на точность инженерного

расчета

возможен

только

методами

те­

ории вероятностей. Правда, иногда, желая получить

верхнюю пре­

дельную

оценку для погрешности

расчета

(грубую,

но гаранти­

рующую,

что в действительности хуже

не будет),

закрывают

гла­

за на случайный характер

факторов и

считают,

что

все они

дей-


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

Учитывая сказанное, задачу оценки точности расчета можно

сформулировать

следующим

образом.

Известна функция

*

 

г =

г

г,

Хо,

 

ХІ,

xk),

(5.3)

которая при х х

= х \ ,

х 2

=

х \ ,

xk

=

х \ представляет

собой ре­

зультат расчета /-„; xt — аргументы функции, которые имеют по­ грешности, упомянутые в пп. 2 и 3 приведенного выше перечня фак­ торов. Величины xt — это, например, геометрические размеры, режимные параметры, свойства материалов, эмпирические коэф­ фициенты в формулах и т.д.

Желательно, чтобы в списке аргументов xt были только неза­ висимые величины, например, если /• = г (х, у), а у = у (х), то надо в качестве аргумента рассматривать только х : г (х).

Если в алгоритме расчета присутствуют эмпирические формулы, содержащие несколько экспериментальных коэффициентов, допус­ тим z — у а — г/Р), то удобнее вместо нескольких погрешностей от­

дельных коэффициентов а , (3, у рассматривать одну общую

погреш­

ность формулы. Для этого, во-первых, к аргументам

х г

надо

до­

бавить одну

новую

величину

Xj с областью изменения 1

Aj

^

<]х;

1 +

Aj, где А; — максимальная относительная

погрешность

рассматриваемой формулы

в

долях

единицы; во-вторых,

исклю­

чить из

аргументов

функции

г все

эмпирические

коэффициенты,

входящие

в данную формулу.* Причем в

алгоритме

расчета

эту

формулу

надо домножить

на

Xj, например

вместо

прежнего z

за­

писать z =

xj

уа"(Уп

/")•

 

 

 

 

 

 

 

 

Итак,

задано г

( х 1 г

х . , ,

xit

x k ) , требуется

оценить

по­

грешность

Аг результата

расчета

rH = г (х", х \ ,

 

xf),

если из­

вестны максимальные погрешности Аг аргументов х п а иногда и законы их распределения / (х-г) в областях

j e ? - A i < * i < * y + Ai.

(5-4)

Как правило, эти законы /(*,•) близки к нормальному (3.9). Ког­ да нет ни практических, ни теоретических предпосылок для приня­ тия в качестве закона / (х[) нормального или любого другого рас­ пределения, обычно в запас можно считать, что xt распределено по равновероятному закону (3.1).

§ 5.2. Два метода оценки точности расчета

 

Поскольку аргументы функции

г — случайные

величины, то

и г — типичная непрерывная случайная величина,

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

по закону / (г), имеющая дисперсию

D (г) = а 2 и

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

* Н е о б я з а т е л ь н о з н а т ь а н а л и т и ч е с к и й в и д ф у н к ц и и г, д о с т а т о ч н о у м е т ь

ВЫЧИСЛИТЬ ЄЄ В ЛЮбоЙ ТОЧКе ъ Х2

A-ft).


ожидание М (л). В этих условиях задача заключается'в следующем: в результате расчета получено некоторое значение случайной вели­

чины /•„ - л (х\, x'i,

.v"), требуется

оценить возможный

разброс

Дг значений

г

около

/•„. Существенно,

что нас интересует

разброс

именно

около

г ш

а

не около М (г),

который характеризует дис­

персия.

Так

что если М (г) совпадает с /-„, то искомая погрешность

расчета

Дг

зависит только от дисперсии

D (г), в противном случае

Дг увеличивается

на

величину |г„ — М

(г)\ .

 

Метод линеаризации. Так как г —функция случайных

аргумен­

тов .г,-, область возможных значений которых достаточно мала (х" —

— Д; < х , .г" + Д,), для вычисления ее дисперсии можно восполь­ зоваться формулой (2.28). Получим

где индекс н

у производных дгідхі означает,

что они вычисляются

в точке (xj,

Л " ,

Л'А). Если аналитический

вид функции г не из­

вестен, то производные следует вычислять по конечно-разностным

формулам,

например

 

 

/ Э М

_

г (х», 4

х? + А;, .... *»)- г{х», х»

xg)

\дХі ) н

~~

 

Д;

 

В формуле (5.5) ри — коэффициент корреляции (2.25) величин хг и Xj, вычисляемый по формуле (4.11). В условиях полного отсутствия данных по ри можно (это будет в запас) положить

Р « = 0 ,

если

dxj

<0,

 

 

(5.7)

 

 

дхі

911= 1'

если

> 0 .

дхі

 

 

 

Для практики инженерных расчетов в реакторостроении интерес­ ны следующие два частных случая.

1. Законы распределения / (хг-) аргументов функции г (5.3) сим­

метричны

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

середины интервала

х" (5.4), т.е. М (х,) =

= х* для

всех і — 1, 2,

k. В этих условиях для линеаризованной

функции г, согласно формуле (2.28), М (г) =

га.

Если законы распределения хг близки к нормальным с математи­ ческими ожиданиями М (хї) ж х", то максимальная погрешность