книги / Методы обеспечения надежности изделий машиностроения
..pdf
|
|
|
|
т |
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