Файл: Волков Е.Б. Основы теории надежности ракетных двигателей.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