НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 • 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Исследование характеристик токарной обработки на основе
получения и анализа многомерных полиномиальных зависимостей
# 08, август 2013
Б01: 10.7463/0813.0589299
Малькова Л. Д., Чихарева М. А.
УДК 621.9
Россия, МГТУ им. Н.Э. Баумана ldm@hotbox.ru chihareva ma@mail.ru
Введение
Математические модели в резании металлов можно классифицировать как теоретические и экспериментальные (эмпирические) зависимости. Теоретические зависимости, представленные в работах Зорева, Розенберга, Клушина и других ученых, опираются на схему стружкообразования и включают в себя механические, технологические и физические свойства материала. Достоинством такого подхода является его универсальность, недостатком -затруднения при определении на практике требуемых характеристик. Экспериментальные зависимости, широко представленные в справочной литературе, являются эмпирическими моделями, полученными в результате математической обработки экспериментальных данных. Наиболее распространены степенные зависимости с корректирующими коэффициентами, которые получались на основе линейной аппроксимации при обработке данных в логарифмических координатах. А также модели, являющиеся произведением площади срезаемого слоя на дополнительные коэффициенты. Оба типа указанных зависимостей получены на основании серий однофакторных экспериментов и должны сопровождаться таблицами коэффициентов. Их достоинствами являются доступные для варьирования на практике параметры зависимости, простота получения и использования. Недостатками -оценочный характер, отсутствие учета совместного влияния параметров зависимости на результирующее значение, большое количество опытов и правомерность только в области проведенных экспериментов.
Самостоятельным направлением в математическом описании процессов механической обработки являются полиномиальные модели, полученные в результате спланированных
многофакторных экспериментов. Такие зависимости, проходя через экспериментальные точки, с высокой степенью достоверности описывают все пространство эксперимента за счет возможности учета совместного влияния параметров на результат. Но полиномиальные модели требуют подготовки и планирования эксперимента, а также обработки и анализа их на базе специального математического аппарата.
Целью данной работы является формирование принципиального алгоритма получения и анализа многомерных полиномиальных зависимостей в резании металлов на примере токарной обработки и сравнение экспериментально полученных результатов с результатами математического анализа.
Научная новизна работы заключается в подтверждении положений, экспериментально полученных в [1]. В частности, математически доказано, что для минимизации энергозатрат на токарной обработке со снятием припуска за два прохода максимально возможный припуск должен сниматься на черновом проходе, обработка должна вестись на максимально возможных скоростях резания и подачах на оборот, причем подача на черновой обработке должна превышать подачу на чистовой.
Для достижения представленной цели и получения научных результатов поставлены и решены следующие задачи:
- экспериментально или на основе литературных источников определить факторы, наиболее сильно влияющие на динамические характеристики токарной обработки, и характер этого влияния;
- составить план эксперимента для получения полиномиальной модели тангенциальной составляющей силы резания и использования ее для модели энергопотребления при точении;
- оптимизировать параметры токарной обработки по критерию минимального суммарного энергопотребления черновой и чистовой обработки на базе специального математического аппарата.
1. Определение факторов эксперимента
При планировании эксперимента фактором называется независимая переменная величина, влияющая на параметр оптимизации. Каждый фактор имеет область определения, т.е. совокупность всех значений, которые он может принимать. Факторы должны быть управляемыми, т.е. позволяющими экспериментатору установить их требуемые значения и поддерживать постоянными эти значения в течение опыта, непосредственно воздействующими на объект исследования, совместимыми и независимыми друг от друга [2].
Особенностью планирования эксперимента с механической обработкой металлов являются технологические мероприятия, позволяющие сохранить условия обработки в процессе
эксперимента близкими к постоянным. Например, для достоверности полученных результатов при исследованиях динамических характеристик при точении необходимо предельно сократить время непосредственного участия режущего лезвия в процессе обработки, чтобы величина износа, установленная в плане эксперимента, была близка к постоянной; вести обработку в области с минимальным отжатием заготовки, вызывающем вибрации и изменение глубины резания; проводить замеры составляющих силы резания на одном диаметре заготовки, чтобы сохранить идентичные условия в зоне резания и т.д. Перечисленные требования определяют необходимость предельно возможного сокращения количества опытов и их продолжительности для определения составляющих силы резания.
Одним из направлений указанного сокращения опытов является планирование дробного факторного эксперимента, подробно рассмотренного в работе [3] для токарной обработки. Но исследования по такой схеме показали, что достоверность результирующей полиномиальной функции во многом зависит от выбранной схемы смешивания оценок и установки на ее основе номеров кодов факторов при смешанном коэффициенте. Неточности и погрешности на этом этапе моделирования дробного факторного эксперимента приводят к неадекватности результирующей модели.
Другим направлением сокращения опытов является углубленный анализ параметров обработки, вводимых в качестве факторов в полном факторном эксперименте. В частности, по литературным источникам было установлено, что на силу резания при точении оказывают влияние параметры режима обработки (глубина резания I, мм; подача на оборот Sо, мм/об; скорость резания V, м/мин), параметры режущего инструмента (передний угол у, град; главный угол в плане ф, град; величина износа режущей кромки по задней поверхности кз, мм) и твердость материала НВ. При этом значение главного угла в плане ф в первую очередь оказывает влияние на перераспределение долей осевой и радиальной составляющих силы резания [4], влияние на тангенциальную составляющую, которая определяет мощность резания и энергопотребление, не превышает 10%; кроме того, значение угла ф часто определено конструктивными особенностями детали и не может подвергаться варьированию. В свою очередь возможность управления передним углом у при использовании сменных многогранных пластин исключена из-за наличия сложного рельефа передней поверхности.
Экспериментальные исследования, подготовленные и проведенные для анализа влияния глубины и скорости резания, подачи на оборот и износа по задней поверхности, как потенциальных факторов полиномиальной зависимости, дали результаты, представленные ниже и проиллюстрированы рисунками 1 ... 5. Во всех случаях рассмотрена токарная обработка стали 45 ГОСТ 1050-74 резцами из твердого сплава Т15К6 ГОСТ 3882-74.
На рисунке 1 представлены экспериментальные зависимости тангенциальной составляющей силы резания от глубины резания при различных подачах на оборот. Общую зависимость можно характеризовать как монотонно возрастающую, близкую к линейной. При изменении глубины резания 0,5 мм до 3,5 мм величина силы резания возрастает более чем в 3...5 раз.
Pz, Н 3000 -
2500
2000
1500
1000
500
0,5
1 1,5
- So = 0,07 мм/об
2 2,5
■ So = 0,23 мм/об -
3 3,5
■ So = 0,39 мм/об
1 мм
0
0
4
Рисунок 1. Зависимость тангенциальной составляющей силы резания от глубины резания при различных подачах на оборот (у = 91 м/мин, НВ 200, у = 10°, ф = 45°, Из = 0).
Влияние подачи на тангенциальную составляющую силы резания при различных скоростях резания представлено на рисунке 2. Аналогично предыдущей зависимость характеризуется как монотонно возрастающая, близкая к линейной. При изменении подачи от 0,07 мм/об до 0,39 мм/об величина тангенциальной составляющей силы резания возрастает более чем в 3...4 раза.
Рг, Н
2000 и 1600 -1200 -800 -400 -0 -0 0,1 0,2 0,3 0,4
So, мм/об
-♦-V = 30 м/мин -■-V = 61 м/ми^ V = 93,5 м/мин = 118 м/мин -Ж-v = 152 м/мин
Рисунок 2. Зависимость тангенциальной составляющей силы резания от подачи при различных скоростях резания (Г = 2 мм, НВ 200, у = 10°, ф = 45°, Из = 0).
Зависимость силы резания от скорости резания не является монотонной в диапазоне технологически допустимых значений аргумента (рисунок 3). Кривая имеет экстремумы, положение которых определяется сочетанием обрабатываемого и инструментального материалов и другими параметрами обработки.
Pz, Н
0 20 40 60 80 100 120 140 160
V, м/мин
Рисунок 3. Зависимость тангенциальной составляющей силы резания от скорости резания (г = 2 мм, & = 0,3 мм/об, НВ 200, у = 10°, ф = 45°, К = 0).
Влияние износа по задней поверхности получено с помощью однофакторных экспериментов при различных параметрах режима обработки. Замеры проводились при эксплуатации треугольных СМП, работавших в одинаковых производственных условиях, текущий износ которых составил соответственно 0; 0,146; 0,235; 0,35 и 0,581 мм (рисунок 4). Результаты экспериментов представлены на рисунке 5.
Рисунок 4. Вид задней поверхности экспериментальных изношенных пластинок.
Рг, Н 3500
3000
2500
2000
1500
1000
500
0
___к ■
II-- -а--- » — -А- -▲
---Я- у -т- *- -1— У
0 0,1 Режим обработки:
0,2
0,3 ■3
4
0,4 5
0,5 ■6 —I— 7
ь 0,6
Ьз, мм — 8
Рисунок 5. Графики зависимости тангенциальной составляющей силы резания от износа по задней поверхности при различных режимах обработки:
Номер режима обработки: Глубина резания I, мм: Подача на оборот £о, мм/об: Скорость резания V, м/мин:
1 2 3 4
0,5 3 0,5 3
0,097 0,097 0,39 0,39
56 56 56 56
5 6 7 8
0,5 3 0,5 3
0,097 0,097 0,39 0,39
139 139 139 139
Во всех сериях эксперимента, согласно графикам, с увеличением износа изменяется сила резания, суммарное увеличение которой составляет от 15% на режиме 4 до 78% на режиме 5. На графиках видно, что при небольшой величине износа может наблюдаться уменьшение силы резания, а затем ее монотонное увеличение на всех рассмотренных режимах обработки.
Экспериментальные исследования твердости поковок позволили принять, что твердость заготовок в партии из одной плавки является одинаковой с некоторыми отклонениями от математического ожидания.
В результате глубина резания I, подача на оборот £о, скорость резания V и величина износа режущей кромки по задней поверхности Из удовлетворяют всем требованиям к факторам полного факторного эксперимента, причем при назначении уровней двух последних следует учесть нелинейный характер их влияния на объект исследования.
2. Составление плана эксперимента
Число опытов N в полном факторном эксперименте определяется числом всех сочетаний уровней факторов, т.е. N = тп, где т - число уровней каждого фактора, п - число факторов.
Таким образом, при п = 4 минимальное число опытов на двух уровнях составляет 16, а на трех уровнях 81. Поскольку увеличение числа опытов может привести к изменениям условий эксперимента с течением времени, желательно ограничиться двумя уровнями его варьирования, если нет априори известной экстремальной зависимости результата от какого-либо исследуемого фактора. В рассматриваемом примере получения полиномиальной зависимости тангенциальной составляющей силы резания при точении уровни подачи и глубины резания могут быть назначены на всем технологически допустимом диапазоне варьирования. При назначении уровней износа и скорости резания должны быть учтены участки монотонности: для скорости нижний уровень должен превышать 50 м/мин; для износа верхний уровень либо не должен превышать 0,2 мм (что соответствует рекомендациям многих производителей), либо сравнивать новый инструмент со значительно изношенным до 0,6...0,8 мм, линеаризируя экспериментальные кривые.
Математическая обработка результатов может осуществляться при равномерном дублировании опытов, при неравномерном дублировании и при отсутствии дублирования. В управляемом эксперименте желательно обеспечить равномерное дублирование.
Уровни варьируемых факторов для проведения экспериментальных исследований представлены в таблице 1, а матрица планирования эксперимента в таблице 2.
Таблица 1.
Значения уровней варьируемых факторов.
Факторы
1 2 3 4
натуральное значение натуральное значение натуральное значение натуральное значение
Уровни и интервалы значение в кодах значение в кодах значение в кодах значение в кодах
глубина подача скорость величина износа по
резания лх1 на оборот Х2 резания Х3 задней х4
г, мм Sо, мм/об V, м/мин поверхност и Кз, мм
Уровень 1 1 -1 0,1 -1 86 -1 0 -1
Уровень 2 5 +1 0,5 + 1 138 +1 0,17 +1
Основной уровень 3 - 0,3 - 112 - 0,085 -
Интервал варьирования 2 - 0,2 - 26 - 0,085 -
Значение в кодах получено по формуле [2]:
X. X. 0
X = -^ , (1)
где ¡х^ - натуральное значение /-того фактора (например, ^ = 5 мм), 0 - натуральное значение основного уровня /-того фактора (100 = 3 мм), 8г- - интервал варьирования /-того фактора (определяется от основного уровня, е{ = 2 мм).
Представленная матрица сформирована для наиболее простой полиномиальной зависимости, в которой не включены степени выше первой для каждого фактора, но представлен весь спектр возможного совместного влияния факторов эксперимента. В общем виде уравнение выглядит следующим образом:
у = К +
К1Х1 + К2 Х2 + Ь3 Х3 + К4 Х4 + К12 Х1Х2 + К13 Х1Х3 + К14 Х1Х4 + К23 Х2 Х3 + К24 Х2 Х4 +
(2)
I К4Х3Х4 I ¿123Х>1«Х2Х3 I ¿1242I 134ХХ ХХ*Х4 I ¿^234*Х2ХХ*Х4 + К?1234ХХ *Х2ХХ*Х4
Таблица 2.
Матрица планирования полного факторного эксперимента 24.
№ Х0 Х1 Х2 Х3 Х4 Х1Х2 Х1Х3 Х1Х4 Х2Х3 Х2Х4 Х3Х4 Х1Х2Х3 Х1Х2Х4 Х1Х3Х4 Х2Х3Х4 Х1Х2Х3Х4 у
1 + + + + + + + + + + + + + + + + У1
2 + + + + - + + - + - - + - - - - У2
3 + + + - + + - + - + - - + - - - У3
4 + + + - - + - - - - + - - + + + У4
5 + + - + + - + + - - + - - + - - У5
6 + + - + - - + - - + - - + - + + У6
7 + + - - + - - + + - - + - - + + У7
8 + + + + + + + + - - У8
9 + - + + + - - - + + + - - - + - У9
10 + - + + - - - + + - - - + + - + У10
11 + - + - + - + - - + - + - + - + У11
12 + - + - - - + + - - + + + - + - У12
13 + - - + + + - - - - + + + - - + У13
14 + - - + - + - + - + - + - + + - У14
15 + - - - + + + - + - - - + + + - У15
16 + - - - - + + + + + + - - - - + У16
3. Обработка результатов эксперимента
Опыты проведены на станке модели 16К20 сборным токарным резцом с пластиной
Т15К6 без покрытия (а = 8°, у = 12°, а' = 6°, ф = 45°) на прутке диаметром 88 мм из стали 45 ГОСТ 1050-74. Замеры тангенциальной составляющей силы резания Р2 проведены на динамометрической установке, включающей динамометр УДМ600, тензометрический усилитель 8АНЧ-21 и компьютер с платой аналого-цифрового преобразователя (АЦП) LA70M4. Реализованный план полного факторного эксперимента с результатами опытов и статистической обработки представлен в таблице 3.
Для каждого опыта, т.е. строки в плане эксперимента, определены среднее значение тангенциальной составляющей силы резания Р2 сред, дисперсия и среднее квадратическое отклонение. С помощью С-критерия Кохрена при 5%-ном уровне значимости принята гипотеза однородности дисперсий, на основании чего определены коэффициенты регрессии в уравнении (2). Все расчеты выполнены по рекомендациям [2]. Далее по критерию Стьюдента при 5%-ном уровне значимости выявлены незначащие коэффициенты: ¿13, ¿23, ¿123, ¿134, ¿234, ¿1234. В результате получена следующая модель:
На основании ^-критерия Фишера для адекватности при 5%-ном уровне значимости принята гипотеза адекватности полученной модели (3).
После преобразований, обратных кодированию, по формуле (1) полиномиальная зависимость тангенциальной составляющей силы резания при точении в указанных ранее условиях выглядит следующим образом:
у = 920,54 + 609,15х1 + 442,79х2 -12,09х3 + 22,40х4 + +288,18х1х2 + 20,08х1 х4 -32,34х2х4 +14,44х3х4 -30,77х1 х2х
1Л2Л4
(3)
Р2 = 133,76 + 55,32г - 16,49£о -- 1066,53Кз + +797,38^0 + 389,64йз + 812,98£оИз + 6^Из - 905,14^0Из
(4)
Таблица 3.
План полного факторного эксперимента с результатами опытов и статистической обработки.
№
я и н а
з
е р
а н и
б у
л Г
т
о р
об об
а
к 2 а
ача
д
о оП
я и н а
з
е р
ь т с
о р
о к С
о
пи ат сс оо нн
к
ей К
к ^ «
и
рм ем
вм
о
п
л е В
е
К 5 К
X
е ч а н з е е
н
д
е р
С
о К ь
ы л и с
е ещ
ю
я
и н
на р
зр
яе
5 ^
а т с о с
к
я,
и
с р
е п с
и Д
е о к с е
§ к
те аи рн де ан
к л е к ет
но
д
е р
С
й
иа
рн тр
ро к- К
е и н е ч
ев
а н
е р
х о
нК
3 е о н
4 и л
я
и р
е т
и р
бк а
а ри
К
ш и Ф й и ре ед
та и
ря
кл
£ ^
т с о н
нат
в к
^ §
д
е и н е ч а н
3 е о н
4 и л б а Т
раи рт
К
ш и Ф я
и ^
и р
к
с о н
нат
в к е
01
02
03
04
05
06
07
08
09
10
11
12
13
14
15
16
0,5
138
0,17
0,5
138
0,00
0,5
86
0,17
0,5
86
0,00
0,1
138
0,17
0,1
138
0,00
0,1
86
0,17
0,1
86
0,00
0,5
138
0,17
0,5
138
0,00
0,5
86
0,17
0,5
86
0,00
0,1
138
0,17
0,1
138
0,00
0,1
86
0,17
0,1
86
0,00
2266
2258
2214
2305
898
677
910
709
454
411
479
520
162
140
159
165
7581
4950
4734
36
320
712
3670
1556
234
519
27
56
12
149
35
104
87
70
69
18
27
61
39
15
23
12
10
0,307
0,322
0,69
2,4
5
5
5
5
6
5
5
5
5
5
7
3
6
Полученная полиномиальная зависимость с высокой точностью описывает и экспериментальные точки, и промежуточные точки пространства эксперимента. Но полином не дает возможности визуальной оценки влияния каждого из параметров обработки (факторов эксперимента) на результирующее значение, а определение минимального или максимального
значения требует использования дополнительного математического аппарата. В частности, для определения экстремума полиномиальной зависимости (функции) необходимо решить систему уравнений, являющуюся совокупностью приравненных к нулю частных производных по всем аргументам (найти градиент функции и приравнять его к нулю). Но особенностью инженерных задач является то, что их математическое решение должно удовлетворять техническим и технологическим критериям, поэтому при решении на каждый из аргументов накладываются ограничения.
Поскольку на границе заданной области отсутствует полноценная 5-окрестность, необходимая для математических понятий максимума и минимума функции, задача нахождения минимального значения функции с учетом заранее наложенных ограничений делится на два самостоятельных алгоритма: нахождение минимального значения на границе и внутри заданного пространства [6].
Поиск минимального значения на границе ограниченной функции может быть осуществлен по методу множителей Лагранжа - методу нахождения минимума функции / (Х) относительно наложенных к ограничений, где каждая точка п-мерного пространства функции /(Х)описывается вектором Х, т.е. Х е Яп - с наложенным дополнительным условием Каруша-Куна-Такера.
4. Решение задачи оптимизации параметров токарной обработки по критерию минимального суммарного энергопотребления черновой и чистовой обработки
В работе [1] было экспериментально установлено, что степень влияния соотношения припусков чернового и чистового проходов на тенденцию изменения суммарного энергопотребления зависит от соотношения подач при обработке. Решение на основе перечисленных выше математических методов не позволит наглядно отобразить какие-либо тенденции, но позволит решить оптимизационную задачу, рассматривая еще и факторы скорости резания и износа в многомерном пространстве эксперимента. Расчеты произведены для обработки диаметра 90 мм на длине 100 мм с припуском под два прохода 5 мм на сторону аналогично однофакторным экспериментальным исследованиям, представленным в [1].
Энергия Е, затрачиваемая на один технологический проход при токарной обработке:
Р • V • г
2_О
3,6-106
Е = Р2 ■ V • го, Дж или Е = ^ ^ 1 О , кВт ч (5)
где го - основное технологическое время, мин; го =-1-, I - длина обработки, мм, пшп -
п • X
' '»ли 1 ' ,
шп о
частота вращения шпинделя, об/мин.
^ п • ё • пшп
Принимая, что V =-—, где а - диаметр заготовки, мм, имеем:
1000
^ „п^а •п / р -я-а •/
Е = Р2--^--= -, Дж (6)
2 1000 Пшп • 1000 • К)
где Р2 = / (г, £о, V, Кз ) определено уравнением (4).
Суммарное энергопотребление двух технологических проходов записывается как:
р, •па •/ р7 2 • я • (а - 2г1 )•/
Есм = Е1 + Е2 = -+ -Ь-, дж (7)
с>м 1 2 1000 • 5о1 1000 • £о 2 ^
где индексы 1 и 2 принадлежат параметрам чернового и чистового прохода соответственно; Р21 = / (гь ^01, Vl, Из1), Р2 2 = / (гь8о 2, V2, Кз 2 ), учитывая, что припуск чистового прохода
г 2 = г -
Ограничения, для применения их в данном методе, должны быть представлены неравенствами вида gi (х) < 0, где г - порядковый номер от 1 до к. Они представляют собой
технические или технологические ограничения, определяющие область исследования. Причем ограничения, заданные двойными неравенствами, вырождаются в К = 2к ограничений, используемых в методе множителей Лагранжа. Ограничения, накладываемые на аргументы функции (7), определены экспериментальными исследованиями и рекомендациями из справочной литературы и сведены в таблицу 4.
Таблица 4.
Ограничения на аргументы функции (7).
Параметр обработки Нижняя граница включительно Верхняя граница включительно
значен ие ограничение в общем виде ограничение на границе значен ие ограничение в общем виде ограничение на границе
чернивии проход г1 3,5 gi (ti ) = 3,5 - ti < 0 gi (ti )= 0 4,0 g2 (ti ) = ti - 4,0 < 0 g 2 (ti )= 0
0,02 g3 (S0i )=0,02 - S0i < 0 g 3 (Soi )= 0 0,4 g4 (Soi ) = Soi - 0,4 < 0 g 4 (Soi )= 0
П 80 g5 (vi ) = 80 - vi < 0 g5 (Vi )= 0 i40 g6 (vi ) = vi -140 < 0 g 6 (vi )= 0
Кз1 0 g7(h3i)= 0 -h3i < 0 g 7 (h3i )= 0 0,2 g8 (h3i ) = h3i - 0,2 < 0 g8 (h3i )= 0
чисшвии проход ^о2 0,01 g9 (So2 ) = 0,0i - So2 < 0 g9 (So2 )= 0 0,05 gi0 (So2) = So2 - 0,05 < 0 gi0 (So2 )= 0
v2 80 gii(v2 )= 80 - V2 < 0 gii (v2 )= 0 i40 gi2(v2) = v2 -140 < 0 gi2 (v2 )= 0
Кз 2 0 gi3(h32)= 0-h32 < 0 gi3 (h32 )= 0 0,2 gi4 (h32 ) = h32 - 0,2 < 0 gi4 (h3 2 )= 0
Функция (7) удовлетворяет всем требуемым условиям для использования метода множителей Лагранжа: на нее наложены ограничения, она дифференцируемая и гладкая на всей области исследования. Функция Лагранжа примет следующий вид:
1 = / ( Х ) + ( Х ) (8)
;=1
где / ( X ) = / (¿1, 8о1, г1, кз1, 8о2, г2, кз2 ) - исследуемая функция, д , > 0 - множители Лагранжа. Дополнительное условие Каруша-Куна-Такера:
^^ L (, 8*1, V* , И*1, 8*2, V*, И*2, Ц*, , Ц3, Ц4, Ц5, Ц6' Ц7, ' Ц9' Мю' Мц, ^12' ^13' ^14 ) 0
Д, > 0 , (9)
Д, ■ ^ ( ¿8 ' 8*1'v8 ' ' 8*2'v8' К2 ) = 0
( 888888*1 / \ ¿1, 8о1, Vl, И31, 8*2, V2, к32) - точка минимума, а условия gi (¿1, 8о1, Vl, И31, 8*2, V2, к32 ) = 0
представлены в таблице 4.
Для определения экстремумов функции внутри области исследования найдем ее градиент и приравняем его к нулю:
grad Е(¿1'8о1'^8о^v2'кз2) = (10)
Точка экстремума является точкой минимума функции, если ее матрица Гессе - матрица вторых частных производных, описывающая поведение функции во втором порядке -положительно определенная.
Решение систем (9) и (10) с предположением, что все ограничения активы, т.е. д, ^ 0 V,,
и с проверкой матрицы Гессе произведено с помощью программного пакета Ма1ЬаЬ. Результат решения представлен в таблице 5.
Таблица 5.
Координаты минимальных точек на границе и внутри области исследования функции (7).
Параметр обработки ¿1' 8о1' г1' кз1' 8о2' г2' кз 2, Е ^сум,
мм мм/об м/мин мм мм/об м/мин мм Дж
Значение 4,0 0,4 140 0 0,05 140 0 165675
В результате расчетов установлено, что при указанных выше условиях на токарной обработке снятие припуска за два прохода требует минимальных затрат электроэнергии при параметрах обработки, лежащих на границе исследуемой области. Аналитический расчет
подтвердил тенденции, экспериментально полученные в [1]: максимально возможный припуск должен сниматься на черновом проходе, обработка должна вестись на максимально возможных скоростях резания и подачах на оборот, причем подача на черновой обработке должна превышать подачу на чистовой. Расчетная величина затрачиваемой энергии в исследуемой области составляет 165675 Дж или 0,046 кВт-ч.
Выводы.
1. В результате анализа характера влияния параметров обработки на тангенциальную составляющую силы резания при точении с целью введения их в качестве факторов в полный факторный эксперимент получено, что глубина резания t, подача на оборот £о, скорость резания v и величина износа режущей кромки по задней поверхности Из удовлетворяют всем требованиям к факторам полного факторного эксперимента, причем при назначении уровней двух последних следует учесть нелинейный характер их влияния на объект исследования. В полученной экспериментальной зависимости отражены линейные эффекты и эффекты совместного влияния параметров обработки на результат.
2. Математическое исследование полиномиальных многомерных зависимостей с наложенными технологическими ограничениями требует математического аппарата, обеспечивающего решение двух самостоятельных задач: нахождение минимального значения на границе и внутри заданного пространства. Для чего может быть использован метод множителей Лагранжа с наложенным дополнительным условием Каруша-Куна-Такера и проверка матрицы Гессе экстремальных точек внутри пространства.
3. В результате расчетов подтверждены тенденции, экспериментально полученные в [1]. В частности, для минимизации энергозатрат на токарной обработке со снятием припуска за два прохода максимально возможный припуск должен сниматься на черновом проходе, обработка должна вестись на максимально возможных скоростях резания и подачах на оборот, причем подача на черновой обработке должна превышать подачу на чистовой.
Список литературы
1. Малькова Л.Д. Энергосбережение при проектировании технологической операции токарной обработки // Известия ВУЗов. Машиностроение. 2011. Спец. выпуск. С. 43-49.
2. Спиридонов А.А. Планирование эксперимента при исследовании технологических процессов. М.: Машиностроение, 1981. 184 с.
3. Древаль А.Е., Малькова Л.Д. Совместное влияние параметров механической обработки на величину составляющих силы резания // Известия ВУЗов. Машиностроение. 2007. № 8. С. 53-61.
4. Грановский Г.И., Грановский В.Г. Резание металлов: учебник для машиностроительных и приборостроительных специальностей вузов. М.: Высшая школа, 1985. 304 с.
5. Львовский Е.Н. Статистические методы построения эмпирических формул: учеб. пособие для втузов. 2-е изд., перераб. и доп. М.: Высшая школа, 1988. 239 с.
6. Зорич В.А. Математический анализ. Часть 1. 2-е изд., испр. и доп. М.: ФАЗИС, 1997. 554 с.
SCIENTIFIC PERIODICAL OF THE BAUMAN MS TU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
A study of lathing characteristics on the basis of obtaining and analyzing
multidimensional polynomial dependences
# 08, August 2013
DOI: 10.7463/0813.0589299
Mal'kova L.D., Chihareva M.A.
Bauman Moscow State Technical University, 105005, Moscow, Russian Federation
ldm@hotbox.ru chihareva ma@mail.ru
This article considers a possibility and an algorithm of the mathematical study of polynomial dependences which are used for describing characteristics of machining. The authors describe advantages of polynomial dependences obtained on the basis of a complete factorial experiment in multidimensional parameters spaces. For illustrative purposes, an analysis of energy consumption during two-stage lathing was carried out on the basis of the obtained polynomial dependence of the tangential cutting force on processing factors. The influence of the specified factors was established a priori on the basis of conducted experimental studies or literary sources. Conditions and results of all the conducted experiments are presented in this paper along with the algorithm for determining the minimal value of energy consumption.
Publications with keywords: energy consumption, polynomials, a full factorial experiment Publications with words: energy consumption, polynomials, a full factorial experiment
References
1. Mal'kova L.D. Energosberezhenie pri proektirovanii tekhnologicheskoy operatsii tokarnoy obrabotki [Energy saving at the designing technological operation of turning]. Izvestiya VUZov. Mashinostroenie, 2011, spec. iss., pp. 43-49.
2. Spiridonov A.A. Planirovanie eksperimentapri issledovanii tekhnologicheskikhprotsessov [Design of experiment in the study of technological processes]. Moscow, Mashinostroenie, 1981. 184 p.
3. Dreval' A.E., Mal'kova L.D. Sovmestnoe vliyanie parametrov mekhanicheskoy obrabotki na velichinu sostavlyayushchikh sily rezaniya [The joint effect of the parameters of mechanical treatment on the value of the components of the cutting force]. Izvestiya VUZov. Mashinostroenie, 2007, no. 8, pp. 53-61.
4. Granovskiy G.I., Granovskiy V.G. Rezanie metallov [Cutting of metals]. Moscow, Vysshaya shkola, 1985. 304 p.
5. L'vovskiy E.N. Statisticheskie metody postroeniya empiricheskikh formul [Statistical methods of constructing of empirical formulas]. Moscow, Vysshaya shkola, 1988. 239 p.
6. Zorich V.A. Matematicheskiy analiz. Chast' 1 [Mathematical analysis. Part 1]. Moscow, FAZIS, 1997. 554 p.