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

книги из ГПНТБ / Кацев, П. Г. Статистические методы исследования режущего инструмента

.pdf
Скачиваний:
16
Добавлен:
20.10.2023
Размер:
8.47 Mб
Скачать

висимость между теми же переменными. В случае линейной за­ висимости эти две характеристики связи совпадают. Это позволяет использовать величину разности (ii„ — г) в качестве меры отклоне­

ния регрессионной зависимости от линейного вида.

Найдем корреляционное отношение для зависимости (102) стойкости сверл диаметром 28 мм от заднего угла (см. табл. 25).

Для расчета составим табл. 29. Получаем

 

 

*а ы

£ и л - у ) 2

93024,82 =

2114,2.

 

 

 

п— 1

45 — 1

 

 

 

 

 

Для расчета s2 сначала по данным табл. 25

рассчитаем

£

у )2

= 148 217 374,

затем

получаем

 

 

 

 

 

 

£ (У -У )а

148 217 374

33685,1.

 

 

 

 

п — 1

 

45 — 1

 

 

 

 

 

 

 

 

 

 

Теперь

найдем

выборочное

корреляционное

отношение

 

 

■Пв

& Ш _

2114,2

0,25.

 

 

 

 

s {(/}

33685,1

 

 

Таблица 29

 

 

 

 

 

 

 

 

К

расчет у

м еж групповой дисперсии средних

 

 

 

 

Задний угол

Л

Л

 

1Л ~ У \

( л - у у

(Ui-У)"- Л

 

а

 

 

6

О

116

178,33

31 801,55

63 603,10

 

7

5

273

21,33

 

454,96

2 274,80

 

8

3

304

9,67

 

93,50

280,50

 

9

8

319

24,67

 

608,61

17 041,10

 

10

1

271

23,33

1

544,28

544,28

 

11

6

255

39,33

546,84

9 281,04

 

 

у = 294,33

 

 

 

 

£ =

93 024,82

Следует заметить, что выборочное корреляционное отношение удовлетворительно отвечает своему назначению только тогда, когда число наблюдений в каждой группе _(серии) достаточно ве­

лико для образования групповых средних yt. В противном случае

следует вычислять выравненное * корреляционное отношение г), которое представляет собой отношение среднего квадратического отклонения предсказанных групповых средних от общей средней к общему среднему квадратическому отклонению, т. е.

s2 (yi) _ £ (Л —У)2

( 1 1 7 )

Л 52 £ ( * - * ) ’

В литературе укоренился неточный термин «теоретическое».

100

Вычислим ii для нашего примера с помощью табл. 30

S (у уУ~

246306,97

5597,88,

'Ш =

и— 1

4 5 — 1

тогда

 

 

 

~

5597,88

0,40.

 

11

33685,10

 

 

 

Таблица 30

К расчету дисперсии предсказанных групповых средних

Задний угол

 

 

1VI ~ v |

( U i - v y -

{ 4 ~ v y n t

а

 

V I

6

2

47,49

246,84

60 929,90

121 859,97

7

5

206,42

87,91

7 728,16

38 640,84

8

3

303,43

9,10

82,81

248,43

9

28

338,52

44,19

1 952,75

54 677,17

10

1

311,69

17,36

301,36

301,36

11

6

222,94

71,39

5 096,53

30 579,19

у

= 294,33

 

 

£

== 246 306,97

Проверка гипотезы о статистической значимости связи

Какую величину выборочного коэффициента корреляции сле­ дует считать достаточной для статистически обоснованного вывода о наличии корреляционной связи между исследуемыми перемен­ ными? Ответить на этот вопрос помогает знание закона вероят­ ностного распределения /'. В случае совместной нормальной рас­ пределенности исследуемых переменных и при достаточно боль­ шом объеме выборки п распределение г можно считать прибли­

женно нормальным со средним, равным своему теоретическому значению р, и дисперсией

s, W = J£ r r -

(118)

В этом случае для оценки значимости г надо вычислить нор­

мированное отклонение

При значении t, большем, чем tHр при заданной вероятности

и числе испытаний, можно считать нулевую гипотезу отвергну­ той, т. е. признать данное значение г значимым.

101

В нашем примере имеем

г = 0,35

и п

95

s {г} =

1 — 0,352 =

0,09; t

0,35

3,89.

0,09

1^95^1

 

 

Очевидно, что коэффициент корреляции обладает высокой до­ стоверностью, уровень значимости ниже, чем 0,01. Нулевая ги­ потеза, что р = 0, опровергается.

Однако следует учитывать, что при малых значениях п и при малых или близких к ± 1 значениях г это приближение оказы­

вается очень грубым. Поэтому используется тот факт, что величина

 

 

 

 

 

 

t =

V l - г5

 

 

(120)

 

 

 

 

 

 

 

 

 

 

v

при условии р =

0 распределена по закону Стыодента с п — 2 сте­

пенями свободы. Поэтому,

если

окажется,

что

 

 

 

 

 

 

\ r \ V n -

2

*кр.

 

 

( 121)

 

 

 

 

 

 

V 1 — Г“

<

 

 

 

 

 

 

 

 

 

 

 

 

 

то

гипотеза

об

отсутствии

корреляционной

связи принимается.

 

В

нашем

случае

 

 

 

 

 

 

 

 

 

 

 

 

 

| г 1V~п — 2

_0,35 V 95

2 __^

 

. j

 

 

 

 

 

V 1 — Г2

 

~

V 1 — 0,352

~

'

 

По

таблице

работы

[И]

для

=

0,025

(р = 0,95) и / =

=

п — 2 =

95 — 2 = 93

 

находим,

что

tKP =

1,98. Так как

3,58 > 1 ,9 8 ,

 

то

гипотеза

 

р = 0 отвергается.

 

 

Значимость корреляционного отношения г| рассчитывается

точно так же,

как и для г.

 

 

 

 

 

 

указанный метод

 

Перевод

г

в число г.

В некоторых случаях

определения ошибки коэффициента корреляции может оказаться недостаточным в связи со значительным отклонением распределе­ ния г от нормального (особенно при высоких значениях г).

Р. А. Фишер [60] предложил преобразование случайной вели­ чины г в величину 2, распределение которой значительно ближе к нормальному, чем распределение г. Преобразование г в 2 за­

дается

формулой

 

 

 

2 =

4 - ! п4 ^ .

(123)

Имеются также готовые таблицы данного преобразования

[11 ].

Средняя

квадратическая ошибка для 2 вычисляется по формуле

 

s { z l

= Vn —2,

(124)

102

Оценка значимости г может производиться с помощью /, при

этом

Для нашего примера г = 0,35; п = 95. Определяем z по фор­ муле (123), находим г = 0,37. Ошибка для г

s z

= 0,104,

J/95 — 3

тогда

 

0,37

3,55.

0,104

 

Этому значению t при / —93 соответствует вероятностьр 0,99,

откуда следует, что величина г

значима.

Д о в е р и т е л ь н ы е г р а н и ц ы д л я

р. Если достоверность выбороч­

ного коэффициента корреляции г доказана, то с помощью средней

ошибки можно установить доверительные границы для коэффи­ циента корреляции р той генеральной совокупности, из которой

взята выборка.

Когда /г достаточно

велико, а г близко

к 0,5, это

можно сделать

обычным

способом

по формуле

 

 

 

/

Р

^ ""Н^

к

р

(1

Однако, учитывая, что распределение величины z более близко к нормальному, чем распределение г, лучше определять довери­ тельные границы с помощью z по формуле

 

z — /Kps{z} < z 0^ z +

/Kps{z|.

(126)

В дальнейшем z

переводятся в г. Для

нашего примера z =

0,37;

s {2} = 0,104.

доверительные границы

(при Р = 0,95; t

= 2)

Определяем

для z0 генеральной совокупности по формуле (125)

 

0,37 — 2 • 0,104 < 20 < 0,37

+ 2-0,104,

 

0,162 < z0 < 0,578.

Теперь в обратном порядке определяем по формуле (123) гра­ ничные значения г по z и получаем доверительные границы для р:

0,16 <С р <С 0,66. Границы для р очень велики. Кроме того, мини­ мальное значение не подтверждает наличия корреляции.

С р а в н е н и е з н а ч и м о с т и р а з л и ч и я д в у х к о э ф ф и ц и е н т о в к о р р е л я ­ ции . С помощью числа z можно определить значимость различия

между двумя коэффициентами корреляции или между фактически полученным коэффициентом корреляции и теоретически ожидае­ мым, а также провести объединение данных по нескольким корреля­

103

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

 

 

= .

Zl~ 2a

(127)

 

 

 

( z ^

+ S 2 )Z2}

 

 

где zh 22

определяется по формуле

(123);

определяется

no

формуле

(124).

 

 

 

 

Эта функция приближенно удовлетворяет нормированному

нормальному

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

Поэтому сравнивая значения

W

со значением

t 0, взятым по таблице интеграла

вероятности

при

данном уровне вероятности Р, можно принять или отвергнуть

нулевую гипотезу.

Установление зависимости стойкости от ряда параметров на основе множественной корреляции

Если исследуется связь между несколькими признаками, то применяется метод математической статистики, называемый мно­ жественной корреляцией. Теоретически можно построить уравне­ ние множественной корреляции для любого числа факторов. Но практически трудоемкость вычислений резко возрастает по мере увеличения числа факторов. Поэтому при составлении корреля­ ционных уравнений необходима механизация вычислительных работ. Простота и однообразие вычислительных операций (сумми­ рование, умножение, возведение в степень) упрощает составление алгоритма для расчета корреляционных уравнений на электронновычислительных машинах. Разработаны схемы расчета корреля­ ционных уравнений для различного числа факторов, а также схемы использования для этих целей электронно-вычислительных машин.

Расчет уравнения множественной корреляции для трех пере­ менных. При изучении корреляции трех переменных одна из них рассматривается как функция, две другие — как аргументы. В простейшем случае связь эта линейная

м

{у} = Ро +

Pi*! + ря*а-

 

 

(128)

Расчет оценок параметров ро, Ри Р2 производим по способу

наименьших квадратов.

относительно

неизвестных

Ь0, Ьг

Решая систему

уравнений

и Ь2 и учитывая,

что

 

 

 

 

 

получаем

 

i>xSi ~Ь b

12S2 = гю^О)

 

 

 

 

 

 

 

 

 

 

Ь хГ12^1

 

— Г2()Sq,

 

 

 

где s(. — среднее

квадратическое

отклонение

фактора (t

= 1,2);

г12— коэффициент

корреляции

между х г

и

лг2; ri0 — коэффи­

циент корреляции

у с х ± и х 2.

 

 

 

 

104

О т с ю д а

_ г 10

Т 2 ( / 1 2 .

S °

и 1

т \ 2

 

1

( 1 2 9 )

 

Г 1 0 Г 1 2 ,

г 2 0

s o

и 2

Г 2

 

1

^ 2

1

Г 12

 

Разделим первое уравнение системы нормальных уравнений

на п, получим

 

 

Ь0у ^1х1

Ь2х2,

 

где Ьх и Ь2 уже определены формулами (129).

Вычисленные зна­

чения Ь0, Ьх и Ь2 подставляем в уравнение (128)

и получаем урав­

нение множественной корреляции у

от х х и х 2.

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

*1.2. 8. 4...... P = V l ± ’

(130)

где

 

 

 

 

 

 

Г12

Г13

• • г 1р

0

Д* = (— 1!)*У

1

Г23

' ' ' Г2р г 21

Г32

1

■• ■ Г3р

Г21

 

Гр 2

ГР з

. . .

1 г р 1

1

Г23

Г24

 

^2р

 

f 32

1

Га*

• • •

г 3р

.

 

 

 

 

 

Гp i

г р з

Гр4

. . .

1

 

 

 

 

Второй определитель равен первому без первой строки и послед­ ней колонки. Здесь индекс 1 соответствует результативному приз­ наку, а индексы 234. . — факториальному признаку. На основе формулы (130) находим, что для двух аргументов х х и х 2 коэффи­ циент множественной корреляции R рассчитывается по формуле

R1. 23----

^ Г21'Г3\Г33

2

 

Г

 

23

При независимых переменных, т. е. когда г12^

множественной корреляции рассчитывается по формуле

(131)

0, коэффициент более простой

R — V r \i - \ - г\\.

( 1 3 2 )

Применим формулу (130) для расчета коэффициента множе­ ственной корреляции при трех факторах. На основе обработки результатов измерений величин биения по наружному диаметру плашек на отдельных операциях была получена корреляцион­ ная зависимость

* ll234 =

0,045 + 0,705х'2 — 0,370*3 + 0,692*,

(133)

и парные коэффициенты

корреляции

 

 

г12 = 0,48;

г13 =

0,38;

г23 — 0,87;

 

 

ги = 0,69;

гз, =

0,53;

/'24 = 0,61,

 

где * х — величина биения

по наружному диаметру готовых пла­

шек (в конце

процесса);

* 2, * 3,

*4 —

величина биения

соответ­

ственно после операций: автоматная обработка, шлифование на­ ружного диаметра, нарезание резьбы.

Определим коэффициент множественной корреляции для за­ висимости (133) с помощью формулы (130)

0,48

0,38

0,69

0

1

0,87

0,61

0,48

0,87

1

0,53

0,38

/

ш

0,61 0,53

1

0,69

0,074

V

ч

1

0,87

 

0,61

0,15

 

 

 

 

 

0,87

1

 

0,53

 

 

 

0,61

0.53

' 1

 

Отсюда

0,7.

Таким образом, зависимость (133) имеет высокую тесноту связи, это означает, что, пользуясь этим уравнением, можно установить допустимые значения * 2, *3, х4, обеспечивающие нахождение * х в допустимых пределах.

Рассмотрим расчет уравнения множественной корреляции на примере обработки результатов испытаний партии сверл диа­ метром 6 мм в количестве 95 шт. Обработка результатов испытаний сверл в работе позволила получить зависимости стойкости от каж­ дого параметра (всего 13 уравнений), в том числе уравнения, ко­ торые имеют достаточно высокую тесноту зависимости (табл. 31). Также установлено уравнение зависимости между величинами толщины сердцевины К и ширины ленточки /:

К = 0,77 + 0,15/; г = 0,30.

J36

Таблица 31

Уравнения парной корреляции стойкости с параметрами сверл диаметром 6 мм

Уравнение связи стойкости

Коэффи­

Среднее

Пределы

 

с параметрами

циент

значение

колебания

 

 

корреляции

параметра

параметров

Т =

1441,0 — 11,82ср

—0,48

117,2

110—122

Т =

180,58 + 26,5К

0,32

0,9

0,76—1,07

Т =

30,9 + 109,97

0,25

0,86

0,5—1,1

Угол при вершине сверла 2ср оказался независимым от К и /.

Найдем коэффициенты и свободный член уравнения на основе

уравнений

(129):

Т = ь0 +

bkK + bfh

 

 

 

 

 

 

и

б/гт

Г[тrkf

S \т\

0,32 — 0,25 — 0.30

60,9

 

1

r2

s ( / г )

 

 

1 — 0,09

0,074

 

 

1

rkf

 

 

 

 

 

 

 

и

ffr

rk i rkf

s{T)

0,25

0,32-0.30

60,9

 

l

. 2

s

1 / 1

 

1 — 0,09

0,141

 

 

1

r M

 

 

 

 

 

 

b0 =

T — bkk b,[ =

58 — 221,37 • 0,9 — 73,09 • 0,86 = — 204,09.

Уравнение

множественной

корреляции стойкости с К и /

 

 

 

Т = —204,09 +

221,37К + 73,09/.

Коэффициент множественной

корреляции

 

 

 

Ят fk' =у^

~4~ Пт

^ гктГ[тГkf

 

 

 

 

 

1

 

 

 

 

 

0,322 +

0,252 — 2 ■0,32 • 0,25 • 0,3

= 0,58.

 

 

= V -

 

1 — 0,32

 

 

Это значение свидетельствует о высокой тесноте зависимости. Включим в полученное уравнение множественной корреляции третий параметр — угол в плане 2ср. Так как этот параметр яв­ ляется независимым по отношению к параметрам К и /, то для его включения достаточно сложить уравнение Т / х (К, /) с уравне­

нием 7\»= / 2 (2ср) (табл. 31)

г =

— 204,09 + 221.37К + 73,09/

Т =

1441,0—

11,8-2ср

________________

2Т =

1236,91 +

221,37К +

73,09/— 11,8-2<р

Т = 618,45 +

110,68К +

36,54 — 5,9-2ф

Коэффициент множественной корреляции для трех параметров

R = V ^ + Г2фт = V 0,582 + 0,482 = 0,75.

107

Значение R свидетельствует о высокой тесноте связи. Под­

ставив в уравнение средние значения параметров, получим рас­ четное среднее значение стойкости Трасч = 58,04 мин против фак­ тической Тф = 58,0 мни.

Корреляция качественных признаков

икорреляция рангов

Вряде случаев исследуемые признаки не могут быть оценены количественно, а только качественно, или им могут быть присвоены условные порядковые оценки — ранги. В этих случаях теснота связи определяется методами корреляции качественных призна­ ков и рангов. Рассмотрим эти методы на примере решения конкрет­ ной задачи.

Установление критерия технологической пластичности быстро­ режущей стали [33]. Применение методов пластической деформа­ ции при изготовлении заготовок режущих инструментов (прокат сверл, штамповка фрез, долбяков и т. д.) выдвинуло новые тре­ бования к инструментальным сталям, в частности к их пластич­ ности в горячем состоянии. Практика показала, что при изготовле­ нии заготовок инструмента из быстрорежущей стали методами пластической деформации может появляться брак в виде трещин. Возникла необходимость в создании метода оперативного кон­ троля быстрорежущей стали по пластичности для определения возможности проката сверл. Для этого надо было найти такой показатель свойств стали, который можно было бы легко опреде­ лять и который был бы тесно связан с поведением стали при про­ кате.

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

 

nDсо

где D — диаметр рабочей

части образца; 10— длина рабочей

части образца; со — число

оборотов одной головки образца отно­

сительно другой при его деформации до разрушения.

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

Для испытаний образцов на горячее кручение во ВНИИ соз­ даны специальные установки, которые были направлены четырем инструментальным заводам. Исследования проводили следующим образом. Поступающую на заводы быстрорежущую сталь подвер­ гали испытанию на кручение при различных температурах 700,

108

800,900,1000, 1100,1150 и 1200° С. Фиксировали химический состав плавки, величину зерна, карбидную неоднородность, центральную рыхлость, диаметр проката, а также результаты проката сверл из данной плавки стали (наличие или отсутствие брака по трещинам при прокате). Всего с четырех заводов получены результаты по 45 плавкам.

Задача заключалась в определении тесноты связи между каж­ дой характеристикой плавок и результатами проката — наличие или отсутствие брака. Характеристика плавки, имеющая наибо­ лее тесную связь с результатами проката, и может быть выбрана как критерий для оценки пригодности плавки для проката. Учи­ тывая, что результат проката (брак или годный) имеет качествен­ ный характер, для определения тесноты связи применена корре­ ляция качественных признаков.

Таблица 32

Количество случаев брака при прокатке сверл в зависимости от пластичности материала

 

Количество плавок

прокатан­

Общее

Число оборотов образца

ных

 

до разрушения

 

 

количество

 

с браком

без брака

плавок

 

 

со «г: 7,57

19

7

26

со > 7,57

1

18

19

Общее количество плавок

20

25

45

Зависимость результатов проката от каждого из свойств стали (теснота связи между ними) характеризуется коэффициентом ассо­ циации А, который вычисляют следующим образом. Определяют

среднее значение данного факториального признака. Затем ча­ стные значения этого признака разбивают на две группы, в пер­ вую включают значения признака меньшие и равные среднему, а во вторую— больше среднего. Для каждой группы подсчиты­ вают количество плавок с различными результативными приз­ наками (наличие или отсутствие брака при прокате сверл). По по­ лученным данным (табл. 32) рассчитывают коэффициент ассоциа­ ции

А=

ad be

(134)

 

+ с) iP + б) (a -J- b) (с -|- d)

где а и с — количество плавок, прокатанных с браком при числе оборотов образца до разрушения соответственно при со < со и

со >• со; b u c l — количество плавок,

прокатанных без брака соот-

------- 2

о,

ветственно при со « со и со > со;

со =

— -------- среднее число

оборотов образца до его разрушения, рассчитанное из всех испы­ таний; п — общее количество испытанных образцов.

109

Соседние файлы в папке книги из ГПНТБ