Файл: Волков Е.Б. Основы теории надежности ракетных двигателей.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 220
Скачиваний: 0
mK., = 2n_1arcsin q ,-j n, как следует из выражения (2.85), также
приближенно
|
/?Лг = — V |
arcsin Qn . |
(2. 88) |
|||
|
|
ПС |
|
|
|
|
|
|
‘<i |
|
|
|
|
Это позволяет, используя равенство |
(2.88), предположить, |
что |
||||
и в общем случае в выражениях |
(2.86) и (2.87) для определе |
|||||
ния /(.у может |
быть |
использовано |
соотношение |
(2.88), |
где |
|
В том случае, |
когда |
оценка вида |
Р ~ Р „ = Pi, 2, ..., |
jv + ejv, |
где |
едг — приближенное значение |
e,v удовлетворяет условию Р ^ Р П, |
точность приближения Р « Р Пможет быть оценена так: |
|
Р - Л , |
если Р„>- 0,5, |
0 = |
|
min [Р, (1 — Р)] |
|
где Р„ находится из формулы (2. 78) |
|
Соотношения вида (2.78), |
(2.87) удобны на этапе проекти |
рования потому, что они не связаны с допущениями о виде зако на распределения рассматриваемых случайных величин или век торов.
Как отмечалось, в ряде случаев вместо вероятностей Р*= = Р(Л,-) уже при проектировании приходится иметь дело с их
оценками Р,- |
и доверительными интервалами |
[Р,-, Р,-] для |
Р*. |
|||
Тогда |
|
N |
N |
|
|
|
|
/V |
|
|
|
||
,дг |
П А, |
= П Р , + |
Р ,„ - П Р / |
|
(2.89) |
|
|
|
! =1 |
|
|
|
|
и согласно |
соотношениям (2.65) — (2.66) |
нижняя |
граница |
|||
P i,2,---, .v доверительного интервала для Р|,2,..., к = |
Р ( Г) |
А -) |
i=i
приближенно находится из соотношения
_1,2,..,,ЛГ
V ( |
N |
\ 2 |
(2. 90) |
X |
ря - п р / |
2 ( ‘ - & ) |
|
V |
<“1 |
/ i<j |
|
где все Р* находятся при одной и той же односторонней довери тельной вероятности у. При написании выражения (2. 90) учиты-
98
валось, что а~гы(1 — q2) | / //, где п — минимальное из чисел ис
пытаний, по которым найдены оценки по формуле (2.89). Исследуем важный частный случай, когда рассматри
вается случайный вектор uN= (ии и2, .. ., n.Y) с нормальной функцией распределения (1. 100), имеющий вектор средних цлг,
вектор средних квадратических ст.у и корреляционную матрицу
II Qijll •
Пусть вначале требуется вычислить вероятность Р],2..... n =
= Р (Ui>0, V / = 1, /V), |
выражаемую с помошыо многомерного |
нормального интеграла |
(1. 100), где /г,-= p i Решение даже та |
кой задачи вызывает определенные трудности и получено лишь для двумерного и трехмерного случаев, причем для вычисления Pi, 2 п Pi.2,з требуется применение достаточно сложных специ альных таблиц [63]. Задача упрощается, если все q,j одинаковы
по величине и отношения |
/?г = щ/ст,- |
равны [см. соотношения |
||
(1. 101) |
и (1. 102)]. Однако и в этом частном случае для расчета |
|||
Р |,2, |
дг требуется применить соответствующие таблицы, |
а при |
||
р,-,-> 0,5 — их разработать. |
Попытки |
определить Pi,2, .... .у |
с по |
мощью разложения плотности многомерного нормального рас пределения [40] или использования других методов, например, метода приведения матрицы ||д,-ф к диагональному виду, еще не привели при N>2 к получению аналитических соотношений, до статочно простых для применения на практике. Из выражений (2.86) и (2.87) следует, что для вычисления многомерных нор мальных интегралов вида (1. 100) может быть использовано со отношение
Р (й,- > 0, V £= 1, vV) — Р |
< |
VL— 1, |
N 1 = |
N |
N |
KN+°-K, |
(2.91) |
--=П^(А/) + |
/-1 |
||
1-1 |
|
|
где //.,.=дг,./з,.;
да — коэффициент корреляции Щ и и-\
щ, Oi — среднее значение и среднее квадратическое откло нение иг, F(h) — интеграл Лапласа (1. 103).
Фунция Димерного нормального распределения с учетом вы ражения (2. 87) имеет вид (2. 91), где
= (-*; — М/Л-
Кроме того, из выражения (2.87) может быть найдена и веро ятность Р(6,<Нг<а,-, V /= 1 , N), если положить
(2.92)
где |
|
|
//.w = («; — &/V3,-; |
= |
(«;--«,0/3/. |
В табл. 2. 1 помещены некоторые |
из результатов расчетов |
численных значений многомерных нормальных интегралов [73]
при yV^8.
В графе 1 приведены точные значения Pi, 2, . л' = Р(«г>0, V i= l, N), заимствованные из работы [63], в графах 2 и 3—при
ближенные |
значения |
Р ь 2, . |
д- |
и |
Рп= Р ь 2, • |
• •, л' + еЛ', |
где |
|
еN— приближенное значение |
ел'. |
При |
этом принималось, |
что |
||||
ejV= 0,l |
(1—Р 1, о, - • ., n). |
причем поправка для всех /г,- бралась со |
||||||
знаком минус при N=2 , 3, 4, 5 н с различным знаком при N=6, |
||||||||
7 и 8. |
Из |
таблицы |
видно, что уже первое приближение |
|||||
Р,,2__;Л- « Р 1,2......n для |
ряда |
практических задач |
можно |
счи |
||||
тать удовлетворительным. Погрешность формулы |
(2.91) не име |
ет тенденции к возрастанию при увеличении кратности многомер ного интеграла. К тому же выводу приводит сравнение прибли женных и точных расчетных данных при различных Л,- и q,-; Введение поправок и использование оценки Pi,2, •••, jy—Рп по зволяет при необходимости улучшить результат.
Приведенные соотношения намечают некоторые границы воз можного применения приближенных аналитических соотно шении.
Таблица 2. I
Сравнение данных расчета jV-мерных нормальных интегралов, полученных с помощью ЭВМ [63] и по приближенным аналитическим
соотношениям Л1=2,/1|= /ы = Л
|
|
|
h |
|
|
Номер |
Q |
|
|
|
|
|
|
1,0 |
2,0 |
2,5 |
3,0 |
3,5 |
графы |
|
|
0,74487 |
0,95851 |
0,98824 |
0,99738 |
0,99954 |
1 |
0,50 |
0,75236 |
0,96244 |
0,98968 |
0,99775 |
0,99961 |
2 |
|
0,73760 |
0,95868 |
0,98865 |
0,99753 |
0,99957 |
3 |
|
0,77273 |
0,96294 |
0,98936 |
0,99759 |
0,99957 |
1 |
0,75 |
0,77996 |
0,96703 |
0,99095 |
0,99803 |
0,99966 |
2 |
|
0,75992 |
0,96374 |
0,99004 |
0,99783 |
0,99963 |
3 |
|
0,81084- |
0,97053 |
0,99163 |
0,99811 |
0 99966 |
1 |
0,95 |
0,81439 |
0,97723 |
0,99379 |
0,99865 |
0,99977 |
2 |
|
0,79583 |
0,97495 |
0,99317 |
0,99852 |
0,99975 |
3 |
100
|
|
jV |
> 2; q,•j = |
q = 0,5; !u — h |
|
|
||
|
|
|
|
N |
|
|
Номер |
|
Л |
|
|
|
|
|
|
||
3 |
4 |
5 |
6 |
7 |
8 |
графы |
||
|
||||||||
|
0,912 |
0,895 |
0,874 |
0,859 |
0,8444 |
0,831 |
I |
|
1,8 |
0,918 |
0,897 |
0,876 |
0,856 |
0,837 |
0,819 |
2 |
|
|
0,910 |
0,887 |
0,864 |
0,860 |
0,853 |
0,837 |
3 |
|
|
0,929 |
0,912 |
0,897 |
0,883 |
0,871 |
0,860 |
1 |
|
1,9 |
0,93-1 |
0,917 |
0,899 |
0,884 |
0,866 |
0,851 |
2 |
|
|
0,927 |
0,909 |
0,898 |
0,872 |
0,879 |
0,866 |
3 |
|
|
0,943 |
0,928 |
0,916 |
0,904 |
0,894 |
0,884 |
1 |
|
2,0 |
0,948 |
0,934 |
0,919 |
0,906 |
0,893 |
0,880 |
2 |
|
|
0,943 |
0,927 |
0,911 |
0,897 |
0,904 |
0,892 |
3 |
|
|
0,983 |
0,979 |
0,976 |
0,970 |
0,965 |
0,963 |
1 |
|
2,5 |
0,985 |
0,981 |
0,977 |
0,973 |
0,969 |
0,965 |
2 |
|
|
0,983 |
0,979 |
0,975 |
0,970 |
0,966 |
0,962 |
3 |
|
|
0,996 |
0,995 |
0,994 |
0,993 |
0,992 |
0,991 |
1 |
|
3,0 |
0,997 |
0,996 |
0,995 |
0,994 |
0,993 |
0,992 |
2 |
|
|
0,997 |
0,996 |
0,995 |
0,993 |
0,992 |
0,991 |
3 |
|
Однако вопросы оценки точности в |
определении Pi, 2..... n с |
помощью предлагаемых соотношений остаются недостаточно ис следованными и затруднены тем, что с увеличением N вычисле ние Pi, 2..... .V даже с привлечением самых быстродействующих ЭВМ представляет (особенно при q,-j— М) задачу значительной сложности. Вместе с тем сам факт вычисления с помощью этих соотношений многомерных нормальных интегралов с приемле мой для ряда практических задач точностью вызывает необхо димость не только первоначального рассмотрения, но и дальней ших более углубленных исследований в данном направлении.
Пример 2.9. |
Определение показателя надежности |
конструкции |
по не |
|||
скольким зависимым условиям функционирования. |
критериям, |
согласно |
||||
Конструкция |
рассчитывается на прочность по N = 8 |
|||||
которым для |
се успешной работы должны выполняться условия: |
|
||||
|
.Vi < |
i/i; |
*2<//2;...; |
Х х < у я или /H>0; и2>0; . ..; Hjv>0, |
|
|
где |
щ = !/; — Xi\ i= l,8 ; |
М = 8; |
|
|
||
Xi |
и у 1 — действующее |
и допустимое значение осевой перегрузки; |
101