Menu

Исследование влияния основных технологических параметров плоского шлифования синтетических монокристаллов корунда на величину остаточных внутренних напряжений методами дисперсионного и регрессионного анализов

Добровинская Е. Р., Чукаев В. И., Бабийчук И. П.


Как известно, при обработке синтетических монокристаллов корунда абразивно-алмазным инструментом в прилежащем к обработанной поверхности слое возникают напряжения определенной величины [1, 2]. Эти напряжения в большинстве случаев определяют эксплуатационную надежность деталей, изготовленных из монокристаллов корунда. Для гарантируемого обеспечения требуемой эксплуатационной стойкости изделий необходимо знать влияние технологических параметров обработки и характеристик абразивно-алмазного инструмента на величину остаточных напряжений возникающих при обработке корунда. В связи с этим авторами была предпринята попытка найти такие зависимости применительно к условиям плоского торцевого шлифования монокристаллов корунда абразивно-алмазным инструментом зернистостью 7/5 на связке М1.

В качестве факторов, влияющих на величину остаточных напряжений, т. е, независимых переменных, рассматривались следующие:

h — скорость вращения инструмента, об/мин;

P — нормальная нагрузка на обрабатываемую заготовку, кг;

К — концентрация алмазного порошка в алмазоносном слое инструмента, %;

Δτ — интервал между зачистками, мин;

F — суммарная поверхность заготовок, подвергающихся одновременной обработке, мм2.

Методика проведения экспериментов изложена в работе [1].

Цель проведения исследований состояла в том, чтобы определить:

какие из перечисленных факторов влияют на величину остаточных напряжений;

степень влияния этих факторов на величину остаточных напряжений;

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

Ответ на поставленные вопросы дают дисперсионный и регрессионный анализы и тесно связанное с ними планирование эксперимента. Для определения связи между указанными независимыми переменными и величиной остаточных объемных напряжений необходимо было спланировать и реализовать эксперимент так, чтобы можно было установить главные эффекты и парные эффекты взаимодействия. Изучались пять факторов и их вариации на двух уровнях. Поэтому матрица плана строилась таким образом, чтобы главные эффекты и их взаимодействия входили одинаковое число раз. Такой матрицей является четверть реплика полного факторного эксперимента 25 (2)5-2, т. е. проводилось всего 8 опытов вместо 32.

В табл. 1 представлены натуральные значения и условные обозначения независимых переменных на нижнем и верхнем уровнях.

Таблица 1

Значения факторов K1, % P, кг n, об/мин Δτ, мин P, мм2
ну ву ну ву ну ву ну ву ну ву
Натуральное 25 150 1 6 650 5100 1 3 200 400
Кодированное 0 1 0 1 0 1 0 1 0 1

Прнмечание: ну — нижний уровень; ву — верхний уровень.

В табл. 2 представлена матрица плана 25-2 [3] в условных обозначениях.

Таблица 2

Номер
опыта
K P n Δτ F
1 0 0 0 0 0
2 1 1 0 0 0
3 0 0 1 1 0
4 1 0 1 0 1
5 0 1 1 0 1
6 1 0 0 1 1
7 0 1 0 1 1
8 1 1 1 1 0

Таким образом, каждый эффект встречается 4 раза, а каждое взаимодействие 2 раза.

В табл. 3 представлены результаты реализованного эксперимента. Численные значения являются средними значениями трех повторений опыта.

Таблица 3

Номер опыта 1 2 3 4 5 6 7 8
Величина
напряжений,
кг/мм2
1,77 1,50 1,56 1,40 1,30 1.46 1,70 1,17

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

Для ответа на этот вопрос был использован F-критерий с k - 1 и N - k степенями свободы [4].

где x — значение величины остаточных напряжений;

•• — общее среднее;

j — среднее значение главного эффекта или эффекта взаимодействия на данном уровне;

xij — текущее значение;

N — число опытов;

k — число уровней;

nj — число опытов на данном уровне;

индексы

• — номер опыта;

j — номер уровня.

Значения P — критерия для каждого фактора, вычисленные по данным табл. 3, сравниваются с табличным значением Fτ , которое для случая выявления значимости главного эффекта равно (при достоверности 1 - α = 0,75, N = 8, k = 2) Fτ =1,62, а для эффекта взаимодействия (N = 8, k = 4) Fτ = 2,05 [4].

Если табличное значение Fτ больше вычисленного по опытным данным, то влияние этого фактора или сочетания факторов незначимо. Все расчеты проводились на электронной вычислительной машине БЭСМ-4. Результаты определения значимости главных эффектов и парных эффектов взаимодействия приведены в табл. 4 и табл. 5 соответственно.

Таблица 4

Средние значения величин остаточных объемных напряжений, кг/мм2

Уровень
фактора
Факторы
K P n Δτ F
ну 1,58 1,55 1,61 1,49 1,50
ву 1,38 1,42 1,36 1,47 1,46
Коэффициент корреляции -0,53 -0,35 -0,67 -0,053 -0,094
F-критерий 2,87 0,85 5,85 2,4.10-5 1,5.10-5

Таблица 5

Средние значения остаточных напряжений, кг/мм2

Уровень
сочетания
факторов
kP k•n k•Δτ k•F P•n P•Δτ P•F n•Δτ n•F Δτ•F
00 1,66 1,73 1,53 1,66 1,61 1,58 1,66 1,68 1,63 1,63
01 1,50 1,43 1,65 1,50 1,48 1,51 1,43 1,58 1,58 1,35
10 1,43 1,48 1,45 1,33 1,60 1,40 1,33 1,35 1,36 1,36
11 1.33 1,28 1,31 1,43 1,23 1,43 1,50 1,35 1,35 1,58
F-критерий 1,510-2 3,710-2 17,810-2 23,410-2 17,710-2 3,710-2 68,510-2 1,510-2 0,510-2 1,39

Прнмечание: 0 — нижний уровень; 1 — верхний уровень.

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

Задача регрессионного анализа формулируется следующим образом: нужно получить некоторое представление о поверхности отклика факторов, которое в общем случае можно математически изобразить в виде функции

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

Зависимость (1) удобно представить в логарифмической форме, апроксимирующее уравнение которой имеет вид

Задача сводится к определению значений коэффициентов

Однако, учитывая результаты дисперсионного анализа, выражение (2) можно существенно упростить, так как значимыми оказались только два главных эффекта. Следовательно.

Значения коэффициентов b0, b1, b2 определялись методом наименьших квадратов [5], основное положение которого: сумма квадратов невязок должна быть минимальной, т. е.

После проведения всех необходимых расчетов уравнение регрессии можно записать в виде

В уравнении (5) k выражено в %, а скорость вращения инструмента в об/мин уменьшена в 100 раз. Окончательно уравнение (5) запишется в виде

Для проверки адекватности принятой математической модели строят F — критерий Фишера [3]. Им проверяют гипотезу о том, что дисперсия относительно модели значимо пре¬вышает дисперсию опыта против альтернативы о незначимом различии между этими дисперсиями. Если различие незначимо, то гипотеза об адекватности модели может быть принята.

Значение критерия Фишера вычисляют по формуле

где S20 —дисперсия опыта;

S2ад — дисперсия адекватности.

Дисперсию опыта, в свою очередь, вычисляют по формуле

где S2i — дисперсия i-того опыта;

fi—число степеней свободы в i-том опыте, равном числу параллельных опытов ni минус 1, т. е. fi =ni— 1.

Таблица 6

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

№ опыта Условия проведения опыта Номер параллель
-ного опыта
Численное значение, кг/мм2 ( yi — y̅i )2 Расчетный критерий Кохрена
Gp
Табличное значение
Gт
Численное значение свободного
члена
Дисперсия опыта
S20
К, % n, об/мин Р, кг Δτ, мин
I 25 2500 1,35 3 1 1,70 0,0001 0,535 Gт (1,4,5)=0,906 A=1,0417
I
0,0489
2,97 3 2 1,45 0,0676
2,97 1 3 1,99 0,0784
1,35 1 4 1,73 0,0004
Средние значения 1,71 Σ 0,1465
II 25 1800 1,35 3 1 1,78 0,0016 0,334 Gт (1,4,5)=0,906 A=1,0768
II
0,0025
1,35 1 2 1,85 0,0009
2,97 3 3 1,87 0,0025
2,97 1 4 1,77 0,0025
Средние значения 1,82 Σ 0,0075
III 25 650 1,35 1 1 2,14 0,0256 0,347 Gт (1,4,5)=0,906 A=1,0765
III
0,0138
2,97 1 2 1,93 0,0025
2,97 3 3 1,86 0,0144
1,35 3 4 1,98 0,0000
Средние значения 1,98 Σ 0,0415 A̅ =1,0650 0,0218
IV 150 650 2,97 3 1 1,88 0,0663 0,415 Gт (1,4,5)=0,906 A=1,2854
IV
0,0985
1,35 3 2 1,87 0,0715
1,35 1 3 2,50 0,1314
2,97 1 4 2,30 0,0264
Средние значения 2,14 Σ 0,2956
V 150 1800 2,97 1 1 2,00 0,0001 0,616 Gт (1,4,5)=0,906 A=19,2978
V
0,0772
2,97 3 2 1,80 0,0351
3,50 1 3 1,90 0,0076
3,50 3 4 2,25 0,0689
Средние значения 1,99 Σ 0,1117
V 150 1800 2,97 1 1 2,00 0,0001 0,616 Gт (1,4,5)=0,906 A=19,2978
V
0,0772
2,97 3 2 1,80 0,0351
3,50 1 3 1,90 0,0076
3,50 3 4 2,25 0,0689
Средние значения 1,99 Σ 0,1117
VI 150 2500 2,97 3 1 1,35 0,0025 0,500 Gт (1,2,3)=0,998 A=0,9743
VI
0,005
2,97 3 2 1,45 0,0025
Средние значения 1,40 Σ 0,0050 A̅ =1,1853 0,0469
VII 150 650 6,00 1 1 2,50 0,3549 0,592 Gт (1,7,5)=0,727 A=1,1699
VII
0,0996
6,00 1 2 2,10 0,0383
1,00 1 3 1,90 0,0000
2,00 1 4 1,90 0,0000
3,00 1 5 1,64 0,0698
4,00 1 6 1,63 0,0751
5,00 1 7 1,66 0,0596
Средние значения 1,90 Σ 0,5977
VIII 150 2000 6,00 1 1 2,05 0,0058 0,340 Gт (1,3,5)=0,967 A=1,2988
VIII
0,0086
6,00 1 2 2,00 0,0007
6,00 1 3 1,87 0,0106
Средние значения 1,97 Σ 0,0171
IX 150 2960 6,00 1 1 1,84 0,0000 - - A=1,2616
IX
0,0000
6,00 1 2 1,84 0,0000
Средние значения 1,84 Σ 0,0000
X 150 3480 6,00 1 1 1,66 0,0000 - - A=1,1751
X
0,0000
6,00 1 2 1,66 0,0000
6,00 1 3 1,66 0,0000
Средние значения 1,66 Σ 0,0000
XI 150 5100 6,00 1 1 1,60 0,0025 0,348 Gт (1,7,5)=0,727 A=1,1978
XI
0,0996
1,00 1 2 1,80 0,0225
2,00 1 3 1,77 0,0144
3,00 1 4 1,66 0,0001
4,00 1 5 1,64 0,0001
5,00 1 6 1,56 0,0081
6,00 1 7 1,52 0,0169
Средние значения 1,65 Σ 0,0646 A̅ =1,2206 0,0237

Прнмечание: Образцы опытов I, II и III — вырезаны из стержня 1. Образцы опытов IV, V и VI — вырезаны из стержня 2. Образцы опытов VII, VIII, IX, X и XI — вырезаны из стержня 3.


Дисперсия адекватности определяется по формуле

где N — число различных опытов;

ni — число параллельных опытов;

i — среднее арифметическое из параллельных опытов;

ŷi — предсказанное по уравнению значение в этом опыте;

f — число степеней свободы f = N — k — 1;

k — число определяющих коэффициентов регрессии по данным N опытов.

Для оценки адекватности принятой модели воспользуемся результатами ранее проведенных экспериментов [1, 2]. Эти данные в зависимости от различных параметров приведены в табл. 6. Численные значения остаточных внутренних напряжений, приведенные в табл. 6, получены на образцах с разными значениями первоначальных объемных напряжений. При подсчете усредненной дисперсии определяемой по формуле (8) для каждого опыта проверялось однородность дисперсии по критерию Кохрена

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

где S2i — дисперсия опыта.

Результаты расчетов оценки дисперсии опыта и ее однородности приведены в табл. 6. Для оценки адекватности принятой модели изменения остаточных внутренних напряжений в силу разных значений начальных объемных напряжений необходимо подсчитать значение свободного члена в уравне¬нии (5) для каждого уровня значений первоначальных на¬пряжений. Численные значения свободных членов и их средние значения для каждого значения параметра и сочетания параметров приведены в табл. 6 На основании этого подсчитываются предсказанные уравнением (5) значения внутренних напряжений применительно к конкретным условиям.

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

Таблица 7

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

№ опыта Среднее значение опыта y̅i Предсказанное значение ŷi (y̅i — ŷi)2 Дисперсия адекватности S2ад Расчетное значение
критерия Фишера
Табличное значение
критерия Фишера
I 1,71 1,75 0,0016 0,0012 0,055 Fт (2, 27, 5) = 3,3
II 1,82 1,80 0,0004
III 1,98 1,96 0,0004
Σ(y̅i — ŷi)2=0,0024
IV 2,14 1,93 0,0414 0,0973 2,08 Fт (2, 21, 5) = 3,4
V 1,99 1,78 0,0443
VI 1,40 1,96 0,1089
Σ(y̅i — ŷi)2=0,1946
VII 1,90 2,08 0,0158 0,0127 0,536 Fт (4, 80, 5) = 2,5
VIII 1,97 1,82 0,0219
IX 1,84 1,77 0,0042
X 1,66 1,74 0,0068
XI 1,65 1,69 0,0014
Σ(y̅i — ŷi)2=0,0508

На рис. 1 и 2 представлены графические зависимости величины изменения объемных внутренних напряжений от К к n. На рис. 3 представлены графические зависимости изменения остаточных внутренних напряжений от тех же параметров, но для различных уровней начальных объемных напряжений, и нанесены результаты измерения напряжений в предварительных опытах. Относительное изменение остаточных объемных внутренних напряжений может быть определено по уравнению

где σ0, σx — соответственно начальное и текущее значение остаточных объемных напряжений;

nx , n0 — соответственно текущее и начальное значение скорости вращения инструмента;

kx, k0 — соответственно текущее и начальное значение концентрации алмазного порошка в алмазоносном слое инструмента.

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

Проведенный регрессионный анализ позволил установить математическую зависимость влияния концентрации алмазного порошка и скорости вращения инструмента на величину остаточных внутренних напряжений. Изменение скорости вращения инструмента от 650 до 5100 об/мин уменьшает начальную величину напряжений на 18%, а изменение концентрации алмазного порошка от 25 до 150% — на 14%.

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

ЛИТЕРАТУРА

1. Бабийчук И. П., Добровинская Е. Р , Литвинов Л. А., Пищик В. В, Чукаев В. И. Физика и химия обработки материалов в печати).

2. Добровинская Е. Р., Бабийчук И. П., Литвинов Л. А., Чукаев В. И. Монокристалла и техника. Л., вып. 4, 1971.

3. Адлер Ю. П. Введение в планирование эксперимента. М., изд. «Металлургия, 1969.

4. Хикс Ч. Основные принципы планирования эксперимента. М.-Л.,изд. «Мир», 1967.

5. Гутер Р., Овчинский Б. Элементы численного анализа и математической обработки результатов опыта. М., Физматгиз, 1962.