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

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

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

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

за функцию, а объем партии п за независимую переменную,

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

уравнения

 

У = ^ - = 0 , 2 2 +

^

(135)

или

 

 

 

5 = я ( 0,22 +

- ^

)

 

с корреляционным отношением 11 = 0,8.

Зависимость (135) позволяет рассчитывать среднее квадрати­ ческое отклонение s стойкости по величине размаха R и объему выборки п. Расчет на ЭВМ «НАИРИ-2» дал эту зависимость в сле­

дующем виде:

'у = -^- = 0,5882п-°-2Ш,

(136)

или

 

s = 0,5882/г-0’2428/?.

 

В табл. 38 приведены значения отношения

подсчитанные

по зависимостям (135) и (136) и табулированные теоретические значения для нормального закона распределения.

Таблица 38

Сравнение значений

по разным источникам

 

 

 

S

 

 

 

Колнчестпо

испытаний п

 

 

 

 

 

 

 

 

 

 

 

 

 

ПГ

2

4

6

10

14

20

26

30 '

50

100

 

По фор­

0,655

0,437

0,365

0,307

0,282

0,263

0,253

0,247

0,237

0,229

муле

 

 

 

 

 

 

 

 

 

 

(140)

0,694

0,420

0,380

0,336

0,305

0,284

0,267

0,255

0,228

0,19

По фор­

муле

 

 

 

 

 

 

 

 

 

 

(141)

0,705

0,446

0,374

0,314

0,286

0,263

0,248

0,242

0,219

 

Теоре­

 

тиче­

 

 

 

 

 

 

 

 

 

 

ское

 

 

 

 

 

 

 

 

 

 

значе­

 

 

 

 

 

 

 

 

 

 

ние

 

 

 

 

 

 

 

 

 

 

120

Исследование зависимости стойкости режущего инструмента

от балла карбидной неоднородности. Исследование

влияния балла

карбидной неоднородности на механические и

режущие

свой­

ства инструмента выполнялось на

быстрорежущей стали

Р18

с карбидной неоднородностью от 3

до 8 баллов.

Испытания на

износостойкость выполняли в лаборатории изотопов ВНИИ по методике ускоренных испытаний, основанной на использовании радиоактивных изотопов.

Характеристикой режущих свойств инструмента служил удельный весовой износ режущей кромки Ар, выражаемый в мил­

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

В результате корреляционного анализа получены следующие зависимости и коэффициенты корреляции (табл. 39). Увеличение балла карбидной неоднородности инструментального материала способствует уменьшению стойкости инструмента. При этом наиболее сильно эта связь проявляется при подаче s = 0,5 мм/об, когда на всех скоростях резания связь является достоверной.

Таблица 39

Зависимости удельного весового износа Ар и стойкости Т от б а лл а В карбидной неоднородност и д л я однозубой червячной фрезы

(сталь Р18, гл уб и н а t = 5,63 м м , s = 0,5 м м /об)

с ! Я

Ш

 

 

 

 

 

 

 

 

н

Зависимость удельного

5 2

Зависимость

 

 

5

Л 2

с

11

о

износа

от В

 

стойкости от В

 

 

(В)

5

°

о

Др.10* мг/г = /

 

Т мин =

/ (В)

 

 

 

 

 

 

 

а

»

ч

 

 

 

S+5

 

 

о к

и

 

 

 

 

 

б я

У S2

 

 

 

« § v .

 

 

 

 

 

 

 

 

 

 

21

138

5,10+

0,2В

 

0,22

272,1 - - 8,74В

32

140

9,97 +

0,47В

0,21

46,4 +

107,8

 

 

 

 

 

 

 

 

В

40

130

13,87 +

1,92В

0,49

11,96 +

65,27

 

 

 

 

 

 

 

 

В

50

133

253 +

2,92В .

0,42

11,40 +

14,86

В

 

 

 

 

 

 

 

 

60

120

101,93 — 20,86В +

2,85В2

0,33

5,19 +

6,51

В

 

 

 

 

 

 

 

 

Коэффи­

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

г ( Л )

—0,25

0,52

Г

О

СО

0,38

0,30

Зависимости механических свойств стали Р18 от балла кар­ бидной неоднородности: а) ударной вязкости (вдоль волокна) ан кгс-м/см2 от балла карбидной неоднородности В

aH= 0,374 + ^ g - , п = 0,81;

121

пределы изменений параметров: ударная вязкость от 0,5 до 2,2 кгс-м/см2, балл карбидной неоднородности от 3 до 8; б) пре­ дела прочности при изгибе (вдоль волокна) а„ кгс/мм2 от балла

карбидной неоднородности

В

 

а1

327,9

>1 = 0,62;

„0,29

пределы изменения параметров: предел прочности при изгибе от 142 до 278 кгс/мм2, балла карбидной неоднородности от 3 до 8.

Испытания ручных ножовочных полотен. Ручные ножовочные полотна из стали У10А испытывали на ножовочном станке с чис­ лом двойных ходов в минуту, равным 60, и статическим давлением полотна на образец 6 кгс. При испытании разрезалась сталь 45 сечением 20x20, твердость НВ 160— 190. Критерием оценки ре­

жущих свойств полотен служила продолжительность 8-го реза в минутах (у). Обрабатывали результаты испытаний двух партий

полотен, изготовленных на различных заводах.

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

Таблица 40

 

 

 

 

 

 

Р езульт ат ы обработ ки

испы т аний р учн ы х

нож овочных

 

полот ен

 

 

 

 

 

 

 

Л'Ь

Параметр полотна

 

 

Уравнение связи

Коэффи­

Предел ы

 

 

циент

изменения

пар­

x i

 

 

 

</=/ (Х()

корреля-

параметра

тии

 

 

 

 

цин г

в мм

I

Высота зуба х1

 

у

=

5,4 — 2,5а-,

—0,36

0,55—1,0

 

Высота разводки

а.

у

=

0,37 4~ 3,25ха

0,4

0,9—1,4

п

Высота зуба а.

 

У =

6,4 — 4,6а,

—0,4

1,0-1,15

 

Высота разводки

х2

у

=

0.78 + 0,3а2

0,34

1,5-2,7

Дисперсионный,

анализ

В

предыдущем

изложении предполагалось, что наблюдае­

мый

разброс результатов связан лишь со случайными причи­

нами. Именно это предположение лежало в основе всех про­ веряющихся нулевых гипотез. Все подконтрольные факторы поддерживались на одном и том же уровне. Но может возникнуть

122

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

Рассмотрим самый простой случай, когда дисперсия наблюде­ н и й о2 известна заранее и исследуется один переменный фактор Л. Пусть при измерении фактора А получились результаты наблюде-

НИИ X1 , Х2, . . . , Хц.

Найдем выборочную дисперсию

Сравним эту дисперсию, имеющую f = а — 1 степеней свободы,

с генеральной дисперсией а 2. Если s2 от а 2 отличается незначи­ мо, то и влияние фактора А нужно признать незначимым, так

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

тем, что дисперсия суммы двух независимых случайных величин

равна сумме

их дисперсий.

В нашем случае складываются эффект

случайности

(с дисперсией

о 2) и эффект воздействия фактора А

(с дисперсией

ал’), которые

очевидным

образом

независимы.

Поэтому общая

дисперсия

наблюдений

должна

быть равна

а2 + о^. А величина s2 является оценкой этой общей дисперсии. Следовательно,

0

0

9

Oa ^

s~— а“.

Число а;1 называется дисперсией фактора А.

Это название дано по аналогии с обычной дисперсией, но нужно

помнить, что числа х ъ х 2, . . ., хп не являются случайными,

по­

этому 0 л не связана ни с какой случайной величиной. Выбор

с~а

удобен по двум причинам. Во-первых, дисперсия является про­ стейшей мерой рассеяния. Во-вторых, показатель влияния фак­ тора А определен теперь аналогично показателю влияния слу­

чайного фактора (т. е. обычной дисперсии о2), что позволит не­ посредственно сравнивать фактор А с эффектом случайности.

Изучение переменных факторов по их дисперсиям называется дис­ персионным анализом.

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

123

о

факторов А и В. Будем анализировать дисперсию сгц величины

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

по — 0л ~Ь °чз -j~ Gab 0с>

(137)

где cr3i — дисперсия, связанная с влиянием фактора А\ а% — дис­ персия, связанная с влиянием фактора В ; а%в — дисперсия, связанная с влиянием взаимодействия факторов А и В; aj — оста­

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

Рассмотрим этот метод на конкретной задаче. При решении одной из проблем производства режущего инструмента возникла необходимость оценить различие в качестве разных марок инстру­ ментальных материалов и качестве изготовления инструмента двумя изготовителями. Для этой цели изготовители I и II изгото­ вили концевые фрезы из быстрорежущих сталей Р9К5 и ЕМоБСоб. После изготовления были проведены стойкостные испытания. По результатам этих испытаний необходимо было ответить на поставленные выше вопросы. Можно было бы решить задачу пу­ тем оценки существенности различия средних значений испыта­ ний отдельно по каждой группе (по изготовителям и по маркам стали), используя один из статистических критериев (t — крите­

рий, критерий Вилкоксона и др.). Однако такое решение не поз­ волит использовать в каждом из ответов весь объем испытаний, а только его часть и, кроме того, не позволит учесть влияние взаи­ модействия факторов. Если последнее значительно, то можно получить ошибочное решение. Например, если имеется сильное влияние опыта изготовителя в термической обработке какой-то марки стали, то тогда можно ошибиться в оценке качества сравни­ ваемых марок.

Выполним дисперсионный анализ для исследования различия в качестве сталей Р9К.5 и ЕМо5Со5 и качества изготовления инстру­ мента двумя изготовителями на результатах стойкостных испыта­ ний концевых фрез диаметром 16 мм. В результате испытаний по­

лучено

N

= пег значений xijk стойкости

фрез,

где

i — номер

повторения

испытания

(i

= 1,

. . .,

п),

j

— номер

уровня

фактора

А

(/ =

1, . .

.,

с),

k

номер

уровня

фактора

В (k =

1...........г). В

нашем случае

п =

5,

с =

г = 2.

 

 

В табл. 41 представлены данные для дисперсионного анализа,

которому

подвергаются

концевые

фрезы

двух

марок

быстро­

режущей

стали, и качества изготовления двумя изготовите­

лями.

В

графе 1 указаны изготовители,

в последующих столб­

цах указаны

марки

стали.

На пересечении

строк

и

столбцов

указаны результаты испытаний концевых фрез (стойкость в мин) при соответствующей комбинации изготовителей (фактор А) и

124

Таблица 41 Стойкость концевы х фрез диам ет ром 16 м м из двух м а р о к быст рореж ущ ей ст али

 

 

 

 

Марка стали (В)

Изгото-

испыта­

Р9К5

 

витель

 

 

 

 

(А)

ний

 

Расчет 0 (-у и х

*ijk

 

кт

 

 

25,7

0 П =

104,1

19,0

 

19.0

 

 

 

21.3

 

19.0

0 ц =

10 836,8

21.3

 

20,2

 

 

 

15,7

 

20,2

*п =

20,83

23,5

 

37,2

©21 =

187,2

21,4

 

35,7

 

 

 

12,1

II

38.6

 

=

35 043,8

20,1

 

35.7

 

 

 

28,3

 

40,0

*21

=37,5

10,0

©.1 = 291,3

е*

02, = 84 855,69

х Л = 29,1

£ „ * * * =

ЕМо5Со5

Расчет

н х

| ©12 = 100,8

■0J2 = 10 160,6

х12 ' 20

©22 = 86,9

0^2 = 7551,6

*22 = 17,6

© .2= 187,7 0 2, = 35 281,29

Расчет

0 ,. =204,8

А'ц = 20,4

0{. = 41 943,0

0 2. =274,1

*2- =27,4

= 75 130,8

0 = 478,9

х 0 — 18,8

0 2 = 229 345,2

а = 23,9

12 944,74

ю

сл

марки стали (фактор В). Здесь же приведены исходные расчетные данные: сумма значений стойкостей 0 и их квадраты и средние

стойкости х. На основе этих данных рассчитываются суммы квад­

ратов отклонений. Расчеты выполняются по следующим рабочим формулам (исходные формулы для дисперсионного анализа, а также описание его приведено в литературе [13, 48, 54]).

1. Общая сумма квадратов

 

 

 

 

2 х'

0-

 

 

229345.2

 

 

 

 

 

 

 

гсп

2944,74

2-2-5

 

 

 

 

 

 

 

 

 

=

12944,73 — 11467,26 =

1477,47

 

 

с

числом

степеней

свободы

/общ = N

1 =

19.

 

 

 

2.

 

Сумма квадратов для

варьирования

по фактору А (изгот

витель)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

0 ■— ®1 =

J -

117073,8 — 11467,26 =

240,12

 

 

 

 

1

 

ГСП

10

 

 

 

 

 

 

 

с

числом

степеней

свободы

fA ~ с — 1

= 1-

 

 

 

 

 

3.

 

Сумма

квадратов для варьирования по фактору В (мар

стали)

 

V 0 *

 

 

 

 

 

 

 

 

 

 

 

 

-----— —- ^ г

121804,58 — 11467,26 =

713,19

 

 

nr

 

h

 

ГСП

10

 

 

 

 

 

 

с

числом

степеней

свободы

fB = г — 1

= 1.

 

факторов А х В

 

4.

Сумма

квадратов

для

взаимодействия

 

 

 

 

п

 

 

 

 

 

I

 

 

®1

 

 

 

 

 

 

 

 

 

 

пг

 

 

ГСП

 

 

 

=

_1_ 63592,8 — 11707,38— 12180,46+ 11467,26 =

297,98

с

числом

степеней

свободы

fAB ( с —

1)

1)

=

1.

 

5.

Сумма

квадратов

случайных отклонений

 

 

 

 

 

 

 

 

0 -/; =

12944,74— 12718,56 =

226,18

с

числом

степеней

свободы

focr = сг (а — 1) =

16.

 

 

 

Сумма квадратов для случайных отклонений может быть под­

считана также на основе соотношения (137), как разность между

общей суммой квадратов и суммой квадратов по пп. 2,

3,

4,

т. е.

I 5 | = | 1 | — | 2 | + | 3 | + | 4 |.

 

 

 

 

и их

Далее рассчитывается критерий F для каждого фактора

взаимодействия

по формуле

F =

S“

например,

Фа^т°Ра . Так,

 

 

240 12

So ct

 

 

 

для фактора А

получаем г =

16,98. Табличное

значе-

- щ ^ - =

126

мне FKр для факторов А, В и их взаимодействия

А х В находим

по

таблице работы

[11] (для а = 0,05, f 1 =

1 и

= 16), рав­

ное

4,49.

 

 

 

 

Результаты расчетов сведены в табл. 42. Фактические значе­

ния критерия F для факторов А (изготовитель),

В (марка стали) и

для

взаимодействия

факторов (А х В ) больше табличного, следо­

вательно, можно сделать следующие выводы: а)

качество изготов-

Таблица 42

Дисперсионный анализ данных о влиянии изготовителя и марки быстрорежущей стали на стойкость концевых фрез диаметром 16 мм из сталей Р 9 К 5 и ЕМооСоб

 

Сумма

Число

Средний

 

Критерии

 

степеней

 

 

квадратов

свободы

квадрат

 

 

 

Источник варьирования

 

 

 

 

2

 

 

S S

Г

s~

п

ф

а к т о р а

 

 

~

s 2

 

 

 

 

 

 

 

 

 

 

 

О С Т

Изготовитель (фактор Д)

240,12

1

240,12

 

16,98

Сталь (фактор В)

713,19

1

713,19

 

50,44

Взаимодействие А X В

297,98

1

297,98

 

21,07

Случайные отклонения (оста­

226,18

16

14,14

 

 

 

ток)

 

 

 

 

 

 

Сумма

1477,47

19

 

 

 

 

ления концевых фрез более высокое у изготовителя /; б) каче­ ство быстрорежущей стали Р9К.5 выше, чем стали ЕМо5Со5; в) существенность взаимодействия факторов А х В означает, что

качество марки стали по разному проявляется у разных изгото­ вителей, а именно изготовитель II обеспечивает лучшее исполь­

зование возможностей стали Р9К.5, чем изготовитель / (возможно, что изготовитель II лучше, чем изготовитель I, осуществляет терми­

ческую обработку стали Р9К5).

Г л а в а V. Статистическое планирование экспериментов

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

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

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

Лабораторный активный эксперимент предусматривает при­

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

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

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

128

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

Возникновение современных статистических методов планиро­ вания эксперимента связано с именем Р. А. Фишера. В СССР теория планирования эксперимента развивается с 1960 г. под руковод­ ством В. В. Налимова [41—47]. Планирование эксперимента,

vкак и вся математическая теория эксперимента, не является пока еще законченной, стройной теорией и развивается в основном как раздел прикладной математики.

Задача построения математической модели объекта с помощью планирования эксперимента требует формулировки цели иссле­ дования. Эта задача возникает в любом исследовании, но для по­ строения математической модели требуется количественная фор­

 

мулировка цели. Характеристика цели исследования,

заданная ка-

г чественио,

называется

параметром

оптимизации или критерием

[

эффективности. Параметр оптимизации является

результатом

I изучаемого

процесса, его выходом или реакцией (откликом).

 

Изменяя условия процесса, можно получить различные зна­

 

чения того или иного отклика и при необходимости можно опти­

 

мизировать процесс потому или иному отклику, принятому за кри­

 

терий или параметр оптимизации. Параметр оптимизации должен

 

быть эффективным с

 

точки зрения достижения конечной цели

 

данного исследования. Важно, чтобы параметр оптимизации был

 

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

 

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

 

шем случае на двух уровнях (например, годный и брак). Следует

 

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

 

чественный, и ограничивает возможности исследования. Пара­

 

метр оптимизации должен быть однозначным, т. е. заданному на­

 

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

 

оптимизации.

 

 

 

 

 

Следует учитывать, что фактически мы имеем дело с зависи­

 

мостью оценки математического ожидания от факторов. Поэтому

 

важна также статистическая эффективность параметра оптимиза­

 

ции. Это требование

означает, что число различимых состояний

 

в пределах изменения

значения параметра должно

быть макси-

. мально возможным. В

статистике эффективными называются та­

 

кие оценки, которые имеют наименьшую возможную при данных

 

условиях дисперсию. Так, например, при исследовании обрабаты­

 

ваемости металлов их твердость является статистически неэффек­

 

тивным параметром оптимизации, так как ее измерение дает боль­

 

шую ошибку. Лучше принять другой параметр, например, вре­

 

менное сопротивление на разрыв, если известно, что он тесно свя­

 

зан с твердостью.

 

 

 

 

 

Объект исследования может характеризоваться рядом пара­

 

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

I один, если возможно,

то обобщенный. В сокращении числа пара-

j

метров

оптимизации

важную роль

играет оценка

корреляции

 

9 п.

Г. Кацев

 

 

 

129

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