Файл: Клемин А.И. Инженерные вероятностные расчеты при проектировании ядерных реакторов.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: х1г |
у1г |
х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), М (г) = |
га. |
Если законы распределения хг близки к нормальным с математи ческими ожиданиями М (хї) ж х", то максимальная погрешность