Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Методы обеспечения надежности изделий машиностроения

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
10.92 Mб
Скачать

 

 

 

 

т

r>m\

f>m2

 

 

а_ <^-

Z

 

ут |

- , т .

 

 

 

 

Р"

 

 

Ufli

Ь/12

^

 

'П|

 

 

(6. 12)

 

 

 

 

£П,1+ Ш2 ^

2

^

 

т,

+ ">2

 

 

 

 

 

П1+ п2

 

 

 

Ш,

 

П1+ п2

 

Если верна

гипотеза # 0, то вероятность попадания Ро в кри­

тическое множество [с учетом (6.9) —(6.12)]

не превышает а:

 

п\ + п2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ро=

У \P[m/R\ = R2]P[(m\\ т2)^ х /т \-j- т2 = т }=

 

 

 

т=О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п \ +

п2

(

 

 

 

 

 

 

 

« | + п2

 

 

 

\

 

= I

 

я [т//г,

=

Я2](Я' +

Р")

 

£

P [m /tf,+ tf2]

= а.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.13)

В

соответствии

с работой

[4]

для

сопоставления гипотезы

HO= (RI = R 2)

с

конкурирующими

гипотезами

Н\ =(R\ < /? 2) и

Я2 = (/?1 > /? 2 ) с уровнем

значимости а

используют

следующую

систему правил:

 

испытаний задаются парой (т\т2)у то для

если

результаты

определения т' и т" используют формулы

(6.11),

(6.12);

если

m 'c m iC m " ,

то

гипотезу

 

Но

принимают;

в

случае,

когда

же

 

гипотезу Но отвергают и принимают гипотезу Н\\ если

 

то Но отвергают и принимают гипотезу Н2.

 

Для изделий, обладающих высокой надежностью, т. е. при

<С <

1, можно

пользоваться

приближенной формулой

 

 

Р[х = ( т и m2)/m i + т 2 =

 

О л Г К

+ т 2)!

 

 

т ] = -------- ---------— =

 

 

 

 

 

 

 

 

 

 

 

 

т ,\т2\(л, + п2)

 

 

 

 

 

= с,т- +

(

 

 

 

т.

 

 

I

 

 

 

 

 

 

 

 

\

с

*■

)

 

 

 

(6.14)

 

 

 

 

 

 

VЯ| + Л2/

\

Л| + я2/

 

 

 

 

Таким образом, гипергеометрическое распределение заменя­

ют биномиальным. В табл.

5 приложения [4]

приведены значе­

ния чисел т' и т" для случая,

когда а =0,1;

т\ + т 2 = (0-=-30).

В том случае,

когда отношения

- ^ - > 0 ,1 ,

можно

воспользо­

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

пользовать преобразование i/= 2arcsin-\/jc. Как показано в работе

[4], случайная величина eY= 2arcsinу и м е е т распределе­ ние, близкое к нормальному, со средним значением, равным 2arcsin-\/l — и дисперсией D = Так как эти величины

171

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

М = 0 и дисперсией D = -J—Н— -р = П^ П2 • Отсюда следует,

что случайная величина

Q(/nj, m2)= l2arcsin д /- р ----2arcsin"^/—

) V

/I,+/I2

(6.15)

и имеет нормальное распределение со средним значением

М=? 0

и дисперсией D — 1.

 

 

 

Таким образом, если /?,= /?2 , то вероятность

 

 

^ p Q K ,m 2)|

= а ’

 

(6.16)

где и _а_— квантиль нормального распределения; берется из

2

табл. 1 приложения.

Рассмотренный критерий состоит в том, что гипотеза Но = = (/?i=/?2) отвергается всякий раз, когда выполняется нера­ венство

iQ K .m *)! > и

(6.17)

 

2

В этом случае совокупности двух групп данных объединять не следует, так как они существенно расходятся.

Пример 6.4. Пусть на первом этапе до доработки изделие испытывалось в объеме 100 циклов функционирования (/гi = 100) и при этом было зафиксировано четыре отказа (т \= 4). После проведения доработки изделие испытывалось на функциониро­ вание в объеме 400 инк/мж (//, = 400) и при этом было зафикси­ ровано два отказа ( т 2 = 2).

Определить эффективность доработки с уровнем значимости

а= 0,1 и вероятность безотказной работы изделия.

Ре ш е н и е . По полученным результатам испытаний найдем величину

ni

_

100

__

100

_ Q2

 

/г, + п 2

 

"1Q0 +

400

500

 

 

^Из табл.^ 5 приложения для

mi = m2 = 4 + l = 5

находим, что

т' = 0, т"* = 4. Так

как

т' = 0<.т\ = 4

и т " = 4,

то гипотеза

Но о равенстве вероятностей отвергается с ошибкой а =0,1. Следовательно, эффективность доработки подтверждается с

вероятностью 0,90, а надежность изделия будет определяться только результатами испытаний после доработки:

* т | + т 2 = / , + ^ ; т ' = / ' ; т " = / "

172

ка адежкос™ Й о ™ 01™ Mpa6oTK" подтверждена; 2) сцен-

Для проверки гипотезы H0 = Rt = R2 = R можно воспользо­

ваться приближенным методом, когда mi + m2> 3 0 В этом случае используем выборочную функцию

т2

и = != = = = = = = = = •

(6.18)

Л Г ' - Ч ^ + i )

При подтверждении гипотезы Я0 функция (6.18) имеет нор­ мальное распределение, а величину R можно определить по результатам испытаний по формуле

R =

1 -

Я,

 

где

m, -f т2

 

 

(6.19)

 

П\ +

п 2

 

 

Если \ и \ ^ и а, то гипотеза

Но отвергается, если же \ и \ с и а,

то она принимается. Значения иа приведены в табл. 1 приложения. Пример 6.5. В процессе проведения испытаний на первом этапе было 25 отказов (mi =25) при объеме испытаний 1000 цик­ лов (л| = 1000). По окончании первого этапа выполнили дора­ ботку изделия. При проведении второго этапа испытаний в том же

объеме (л2 = Ю 0 0 циклов)

было зафиксировано 20

отказов

( т 2 = 20).

расхождения этих групп

данных

Определить вероятность

иоценить надежность изделия.

Ре ш е н и е . Для расчета критерия в качестве оценки найдем

п

т, + т 2 _

25 +

20

45 .

= 0,0225.

И — П[+п2 ~

1000 +

1000

2000 ‘

 

Подставляя известные значения в формулу

(6.18), получим

 

 

25 _

20

 

 

и —

 

1000

1000

 

0,76.

 

 

 

 

 

д/о,0255(1 -

0,0225) ( д а

+ - д а )

 

Если выбрать а = 0,10, то иа= 1.281

(из табл 1 приложения).

Поскольку ы = 0 7 6 < Ы а = 1.281,

то гипотеза

Н0 принимается.

Вероятности R i и /?2 близки друг

к другу, следовательно, их мож­

но объединить, т. е. надежность

изделия

R = - Р =

т, + тг

п| + п2

 

25

+

20

0,9775.

1 — 1000

+

Ю00

 

173

О т в е т : 1) эффективность доработки не подтверждена с уровнем значимости а = 0,1; 2) надежность изделия Р =0,9775.

6.4. СРАВНЕНИЕ ВЕРОЯТНОСТЕЙ ОТКАЗА ПО КРИТЕРИЮ СОГЛАСИЯ (НЕПАРАМЕТРИЧЕСКИЙ СЛУЧАЙ)

Допустим, что в процессе испытаний фиксируются отказы, вероятности которых в двух группах данных соответ­ ственно равны Р\ и РгСтавится задача сравнить эти вероятности и определить принадлежность их к одной совокупности. Если в результате п\ независимых испытаний отказы появились т\ раз, а в результате лг независимых испытаний т2 раз, то статисти­ ка

(6.20)

где

1=1

т\ + т2

Р, = Р < = Р 2

П\ + п2

Такая статистика при больших значениях пь имеет распре­ деление %2 с (г — 1)-й степенью свободы.

Это означает, что

P [ T < X i - A r - 1)]* 1 - а,

(6.21)

где xi-a(r—1) находят из условия

где

Г ( Г 2

) — гамма-функция

(выбирают по таблицам).

а =

В табл.

4 приложения даны

значения y}_J<n = r — 1), когда

0,001; 0,01; 0,05; 0,1; 0,9;

п = {г— 1)= 1-f-30.

 

Таким образом, исходная

гипотеза Я0, состоящая в том, что

вероятности отказов равны заданным значениям Л, принимается при статистике T(m|, тг)<Х ?-а(г— 1) и отвергается в противном случае. Вероятность ошибки первого рода приближенно равна а.

Критерий х2 позволяет сравнить вероятности отказов не толь­ ко двух групп данных, но и совокупности групп. В этом случае статистика выглядит так:

* Статистика — функция результатов наблюдений, используемая для оценки па^а^ет^ов^ распределения и (или) для проверки статистических гипотез

174

 

 

T(m„ma, ... ,m r) = J^ ^

,

(6.23)

 

 

 

/=1

 

 

где

Pt =

Я, = P2 =

= P r = P\

 

 

 

P _

m, + m2+

+ mf

 

 

 

 

Л1“Ь Л2 "1“

+

 

 

Ошибка

принятия

гипотезы Н о о равенстве

P i = P

также рав­

на

ос.

 

 

 

 

Критерий согласия %2 может быть использован при доработ­ ках изделий по оценке вероятностей как двух групп данных, так и совокупности нескольких групп.

Пример 6.6. По данным примера 6.4 оценить эффективность

доработки по критерию согласия

%2 с уровнем

значимости (или

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

а = 0,1.

 

 

 

Р е ш е н и е : mi = 4 ,

/11 =

100, m2 = 2, л2 = 400.

Воспользуемся формулой (6.20), предварительно определив

величину

 

 

 

 

 

 

т\ +

т2 _

4 +

2

6

0,012.

П| +

п2

 

100 +

400

500

 

 

Подставим исходные

данные

в формулу

(6.20) и получим

статистику

 

 

 

 

 

 

Т(ти щ ) = (OTi — «Ipf

('Щ n2pf

_ (4 — 100-0,012)г

 

 

 

 

 

 

100- 0,012

(2 - 400-0,012)г _ Q ,g

'400-0,012

Число степеней свободы г = 2 — 1= 1.

Для уровня значимости а = 0,1 по табл. 4 приложения находим

XL»,. ( '= ! ) = 2,71.

Полученное значение Т(т\, т 2)= 8,16>2,71 свидетельствует о том, что гипотеза Н о о равенстве вероятностей Р \ = Р 2 отвер­ гается, т. е. проведенная доработка является эффективной. Для оценки надежности изделия необходимо использовать статисти­ ческие данные, полученные после доработки. Таким образом, критерий согласия %2 подтверждает справедливость критерия сравнения вероятностей, основанного на гипергеометрическом распределении (см. параграф 6.3).

Пример 6.7. Изделие подвергалось испытаниям в три этапа. На первом этапе испытаний было проведено 100 циклов испы­

таний (п,\ = 100) и зафиксировано два отказа (mi = 2 );

на вто­

ром этапе проведено 200 циклов испытаний (п2 = 200)

при трех

175

отказах (m2 = 3) и на третьем — 200 циклов испытаний (яз = 200) при одном отказе ( т з = 1 ) . Оценить принадлежность получен­ ных результатов испытаний к одной совокупности с уровнем значимости а = 0,1. Найти надежность изделия.

Р е ш е н и е . Определим вероятность отказа по формуле

 

 

 

Щ + гп2+ тъ

2 + 3+1

0, 012.

 

 

 

 

п\ + п2+ П3

100 + 200 + 200

 

 

 

 

 

 

 

Подставим исходные данные в формулу (6.20) и найдем ста­

тистику:

 

 

 

 

 

 

 

 

Т(т,.т„ щ ) =

У {" ‘ ~ f f

= <г ~ '0°-„0;0‘2)’

+

 

v

1

2

Li

n.tp

100-0,0.12

'

 

 

(3 - 200-0,012)2 .

(1 -

200-0,012)2

1 м

 

 

 

 

200-0,012

200-0,012

 

 

 

По табл. 4 приложения находим

значение

%2 для

а = 0,1

и

г= 3— 1=2;

xo,i(2)= 4,61.

m2, m3)= 1,84<xo,i(2) = 4,61,

то

Так как

статистика

Т(тi,

следует объединить результаты всех трех этапов испытаний в одну совокупность. Для оценки надежности изделия следует взять

результаты

всех этапов испытаний:

 

 

 

*з =

1 - ^ =

1 - + Г =

° ’995>

 

 

щ + щ + щ

2 + 3+1

0,988.

 

п\ Н " п2 + п3

100 + 200 + 200

 

 

От в е т .

1) результаты

испытаний

следует объединить в

одну совокупность; 2)

надежность изделия равна R = 0,988.

6.5. КРИТЕРИЙ КОЛМОГОРОВА

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

= +

(6-24)

где Xi<x<Lxu+\\ здесь x = ti/T\ — момент остановки испытаний из-за отказа; Т — суммарная наработка изделия; т , — число отка­ зов на /-м интервале (/= 1 ,2 , , /); п — общее число испытаний в циклах.

176

После построения эмпирических функций Ft{x)* и /^(х)* (их можно представить в виде таблиц или графиков) в соответству­ ющих точках этих функций находят максимальное отклонение

A, = sup

1

| Fu (х) — Fa (x)\ =

max

(6.25)

 

 

n

 

Проведенные A. H. Колмогоровым исследования показали, что

{

 

О

при у < 0;

 

( - 1

при у ^ 0.

т -

i

 

k =

- o o

(6.26)

 

 

 

Значения функции К(у) даны в табл. 7 приложения [4]. Вы­ числены также некоторые квантили случайных величин Dn, задава­ емых формулой (6.25). В табл. 7 приложения приведены значе­ ния Dna при а = 0,2; 0,1; 0,05; 0,02; 0,01 и л = 1-М00. Числа Dna удовлетворяют условию

P (D > D na) = а.

(6.27)

Если /г > 100, то проводят нормировку /г, умножая

Dn на -yfn.

Гипотезу о принадлежности двух групп данных одной совокупно­ сти отвергают с уровнем значимости а при Dn-\fn> К\~а [здесь К{КР)= р, т. е. КР— квантиль уровня р для распределения Колмо­ горова]. Если п ^ Ю 0, то Dn сравнивают с числами Dna (табл. 7 приложения). Гипотезу о принадлежности двух групп одной совокупности отвергают с уровнем значимости 1—а, если Dn > Dna- В противном случае гипотезу принимают.

Пример 6.8. Период испытаний отрабатываемого изделия был разбит на два этапа, причем после завершения первого этапа была проведена доработка. Сведения о результатах испытаний сведены в табл. 6.2. Найти эффективность доработки. Оценить надежность изделия.

Р е ш е н и е . По формуле (6.25) с использованием данных табл. 6.2 найдем максимальное отклонение эмпирических функ­ ций:

Dn = 0,25 - 0,15 = 0,1.

Из табл. 7 приложения для величин п = 20 и а = 0,1 находим Dna= 0,265. Так как Dn< D na, то гипотеза Но о принадлежности двух групп данных одной совокупности подтверждается, т. е. для оценки надежности необходимо обе группы данных объеди­ нить:

* f ,(*) = /=■„; F2(X)= F>2.

177

 

6.2. Результаты испытаний (к примеру 6.8)

 

Номер

Число

 

 

f . w

Номер

Число

f.W

испытания

отказов

 

 

испытания

отказов

П\

ТП\1

 

 

 

П2

m2i

 

1

 

 

 

1

 

2

 

 

 

 

2

1

 

3

 

 

 

3

 

4

 

 

0,05

4

0,05

5

1

 

 

5

6

_

 

 

 

6

_

 

7

 

 

 

7

 

8

 

 

 

8

 

9

 

 

 

9

0,05.

10

1

 

 

0,1

10

11

_

 

 

 

11

_

0,10

12

 

 

 

12

1

13

 

 

 

13

 

14

 

 

 

14

 

15

 

 

0,1

15

1

0,15

16

1

 

 

0,15

16

 

 

17

 

 

0,20

17

 

18

1

 

 

18

 

19

 

 

0,25

19

0,15

20

1

 

1

20

1.

 

 

 

-L--------------- 1

£ =

т,

+

т2

-

JL±JL=

1 - - £ г =

0,8.

п |

 

пг

 

+

 

20+ 20

40

 

О т в е т. 1) гипотеза Н0 о принадлежности двух групп данных одной совокупности подтверждается с вероятностью 1—а =0,9; 2) оценка надежности равна /?=0,8.

6.6. ПРОВЕРКА ГИПОТЕЗЫ О РАВЕНСТВЕ ЗНАЧЕНИЙ ДВУХСРЕДНИХ ИЗ НОРМАЛЬНО РАСПРЕДЕЛЕННЫХ

СОВОКУПНОСТЕЙ

На практике эта гипотеза может быть использована при оценке эффективности доработок изделий. Наиболее прием­ лемым является случай, когда выборки получены по данным независимых испытаний, т. е. результаты испытании до прове­ дения доработки не зависят от испытаний после проведения доработки. Суть этого критерия состоит в следующем. Предпола­ гается, что измерения являются выборками из двух нормально распределенных совокупностей со средними значениями pi и \i2 и дисперсиями о2 и Для проверки гипотезы_по = р-1 = Ц2 оп­ ределяют средние выборочные значения х\ и хг параметров и средние выборочные их отклонения S\ и S2 по результатам двух этапов испытаний.

178

Критерием оценки является выборочная функция (критерий Стьюдента)

t =

X, — Х 2

(6.28)

 

/*1 ~ 0 + 52 (П2— О 4

п\ + ^2 *

л/^

(П, + п 2 - 2)

 

где

л - * , ) 2;

(6.29)

 

 

 

( X i - x 2) 2

(6.30)

Величина /, используемая при проверке гипотезы Но, имеет

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

Стьюдента с

пг = п\-\-П2 —2

степенями свободы.

Для заданной

вероятности

а ошибки при

т = п\-\-П2

2 сте­

пенях свободы можно по табл. 8 приложения на^ти значения

ta,m- В случае, если рассчитанная по

формуле

(6.28) величина

\ t \ '^ t a,m, то гипотеза Но

о

равенстве

средних

значении

отвер­

гается; в противном случае, т. е. при |/|<!^а.т,

гипотеза

Но при­

нимается. Формулу (6.28)

можно видоизменить так.

 

 

Х\

%2

4

 

 

(6.31)

 

 

 

 

 

 

где

^ т а в и т с и з а д а ч а с р а в н и т ь upwM-»- эффективность доработки

отказами этапов испытаний и ^ ен^ходные данные испытаний с уровнем значимости а = 0,05.

приведены в табл. 6.3.

179

6.3. Результаты испытаний (к примеру 6.9)

 

П е р в ы й э т а п и с п ы т а н и й

 

В т о р о й э т а п и с п ы т а н и й

Н о м е р

Н а р а б о т к а

 

Н о м е р

Н а р а б о т к а J C 2(

 

Хм

(х, Хм)2

( х 2 X2/)2

о т к а з а

о т к а з а

м е ж д у о т к а з а ­

м е ж д у

о т к а ­

ТП\1

 

П 2/

ми, циклы

 

з а м и ,

циклы

 

 

 

 

 

 

 

1

15

75 625

1

10

130 321

2

25

70 225

2

50

10 272

3

40

62 500

3

100

73 441

4

100

36 100

4

150

48 841

5

150

19 600

5

200

29 241

6

200

8 100

6

250

14 641

7

250

1600

7

350

441

8

280

100

8

450

6 241

9

320

900

9

500

16 641

10

400

12 100

10

600

52 441

11

450

25 600

11

800

184 041

12

480

36 100

12

1000

395 641

13

500

44 100

 

 

 

14

550

67 600

 

 

 

15

600

96 100

 

 

 

Р е ш е н и е . Найдем средние значения наработок и их средние квадратические отклонения по формулам (6.29) и (6.30):

J __ 15+ 25 + 40+100+150 + 200 + 250 + 280 + 320 + 400 + 450 +

Л, ----------------------------------------------------------- --------------------------------------------------------------------------------------------------------------------------► -

 

- + 480+500+550 + 600 =

^

циклов;

~ __ 10+50+100+150+200 + 250+350+450 + 500+600 + 800+1000

2

 

12

 

 

 

 

 

— 371

цикл;

 

 

 

 

5,

, - X,)2 =

199,34; S? = 39736,43;

S2

(x2i -

x2f

=

295,75;

S2 = 87473,0.

_3атем подставим полученные значения

в

формулу (6.28) и

найдем

 

 

 

 

 

 

_______ (290 - 371)______

15-12

0,847.

_

/39736.43-14 + 87473 -7Т

15 +

12

 

 

 

 

 

“ V

----------- Г5_нГ Т 2 _=Г'2

 

 

 

 

 

180