УДК 621.9
А. А. Жолобов, канд. техн. наук, проф., А. В. Казаков
ОСОБЕННОСТИ РАСЧЁТА УПРУГИХ ДЕФОРМАЦИЙ ТЕХНОЛОГИЧЕСКОЙ СИСТЕМЫ ПРИ ТОКАРНОЙ ОБРАБОТКЕ СТУПЕНЧАТОГО ВАЛА
Рассмотрены вопросы оценки погрешностей обработки от деформации технологической системы при точении ступенчатого вала на станке с ЧПУ. Использован подход, базирующийся на методе профессора В. С. Корсакова. Разработана математическая модель, которая позволила выявить и оценить влияние деформации заготовки, инструмента, задней и передней бабок станка на точность и форму обрабатываемых поверхностей ступенчатого вала.
При токарной обработке ступенчатого вала по мере удаления припуска и продвижения резца от задней бабки к шпиндельному узлу станка изменяются жесткость и деформация заготовки, величина упругой деформации передней и задней бабок. В результате получается деталь с переменным диаметром, а также с погрешностями в продольном сечении. При неблагоприятных условиях величина данной погрешности может превысить допуск на размер детали [1].
Целью работы является исследование погрешностей формы обрабатываемых поверхностей в результате деформаций элементов технологической системы. Для этой цели была разработана математическая модель. В созданной математической модели были учтены следующие параметры: отклонение от круглости
(овальность) и погрешности профиля продольного сечения (конусность) исходной заготовки; смещение центровочных гнезд заготовки (при обработке в центрах); исходная шероховатость поверхности заготовки; погрешность базирования заготовки в патроне; смещение оси задней бабки; податливость шпиндельного узла; податливость задней бабки; жесткость инструмента; упругая деформация заготовки.
При разработке модели были сделаны некоторые допущения: материал поверхности заготовки имеет однородную структуру и равномерную твердость; шпиндельный узел и задняя бабка имеют одинаковую величину жесткости в направлении осей У и 2.
Разработанная математическая модель на основе сведений о материале заготовки, размерах и шероховатости её поверхности, значениях режимов резания, геометрии инструмента, жесткости технологической системы, а также погрешности базирования позволила определить значения упругой деформации в каждой точке её поверхности по методу профессора В. С. Корсакова. Суть метода заключается в следующем. Если в статическом состоянии настроить станок на требуемую (расчетную) глубину резания X, то при резании действительная глубина XДЕЗ окажется меньше расчетной на величину упругой деформации, т. е.
Д =1 - 1 Дез . (1)
При определении действительной глубины резания XДЕЗ в зависимости от
способа базирования заготовки возможны три варианта расчета: при базировании в патроне; при базировании в патроне с поджимом задним центром; при базировании в центрах.
Во всех трех случаях на величину действительной глубины резания X ДЕЗ
оказывают влияние следующие факторы: жесткость технологической системы (шпиндельного узла ]Ш, задней бабки _|Б, инструмента и заготовки); смещение правого и левого центровочных гнезд относительно оси заготовки ДСП и ДСЛ, смещение оси задней бабки ДСБ, погрешность базирования заготовки при обра-
ботке в патроне ДП; величина овальности ДО и конусности ДК поверхности заготовки; её исходная шероховатость Я^з.
В общем случае глубина резания X для каждой точки поверхности заготовки определяется как разность между радиусом заготовки Я3 и радиусом получаемой детали Яд. Эта глубина резания отличается от действительной глубины резания X ДЕЗ, имеющей место в реальных условиях, на величину деформации технологической системы. Данная величина деформации в математической модели была учтена при расчёте действительного радиуса детали ЯД . Величины исходной погрешности заготовки и погрешности базирования также учтены при расчёте действительного радиуса заготовки ЯД . Таким образом,
действительная глубина резания X ДЕЗ может быть найдена из соотношения [1]
1ДЕЗ — КД — КД + А
і Д Д
(2)
где Я Д, ЯД - действительные радиусы
детали и заготовки, рассчитываемые для любой точки поверхности заготовки, соответственно, м; А- величина, на которую изменяется глубина резания вследствие наличия исходной шероховатости поверхности Я^, м.
Учёт упругой деформации заготовки. Величину прогиба в любой точке заготовки, при базировании в центрах, патроне или в патроне с поджимом задним центром [1] можно определить с помощью дифференциального уравнения изогнутой оси балки [2]
А ЗУ® (х) — '
М У(Г) (х)
Е1(х)
(3)
где АЗУ(2)(х) - закон изменения упругих
деформаций заготовки в рассматриваемой плоскости; МУ(г)(х) - изгибающий момент, вызванный соответствующей силой или силами резания в рассматриваемой
плоскости; Е1(х) - закон изменения жёсткости при изгибе заготовки.
Уравнение (3) описывает закон изменения упругой деформации по всей длине обрабатываемой заготовки, но значение упругой деформации, определяющее точность цилиндрической поверхности после механической обработки, будет находиться в точке приложения силы резания. В любой момент времени в точке касания вершины резца и обрабатываемой поверхности находится значение упругой деформации заготовки, которое впоследствии оказывает влияние на точность обработанной поверхности. Для этого необходимо привести уравнение (3) к виду, который позволит определять значение упругой деформации в заданной точке.
Для решения дифференциального уравнения (3) резец был зафиксирован в произвольный момент обработки, а его действие на обрабатываемую заготовку заменено на соответствующие составляющие силы резания. Это состояние представлено на рис. 1.
Для дальнейших пояснений обрабатываемую заготовку удобно разделить на участки с постоянной жёсткостью (рис. 2).
Каждая ступень является участком, который характеризуется следующими параметрами: угол поворота в начале участка; угол поворота в конце участка; упругая деформация в начале участка; упругая деформация в конце участка; изгибающий момент в начале участка МУ^; изгибающий момент в
конце участка МУ^; диаметр участка
Б;; длина участка 1; (рис. 3).
В результате интегрирования выражения (3) получен закон изменения угла поворота [2]:
гМУ(7)(х)
0(х) = Дзу(7)(х)1 = { + 80. (4)
Рис. 1. Схема нагружения вала силами резания
Рис. 2. Разделение обрабатываемого вала на участки
Рис. 3. Параметры участков
Исходя из вышеизложенного для каждого из участков справедливо выражение
1 г
8 ;(х)=ЕГ ■!М у(7);(х)+8 “, (5)
где 91(х) - угол поворота в поперечном сечении заготовки с координатой х (0 < х < 11) 1-го участка, рад; МУ(2)1(х) -
закон изменения изгибающего момента на 1-м участке, Н-м; 9 01 - угол поворота в
начале 1-го участка, рад.
Закон изменения изгибающего момента на 1-м участке можно представить в виде линейной функции
Мк - Мн
МуИ(х) — М" + У(2)1. У(2)1 • x, (6)
1.
где М У(7); - изгибающий момент в начале 1-го участка, Н-м; МУ^ - изгибающий момент в конце 1-го участка, Н-м; ^ -
длина 1-го участка, м.
После подстановки М;(х) в (5) и интегрирования полученного выражения (7) угол поворота в поперечном сечении с координатой х (0 < х < 1;) 1-го участка оп-
ределяется выражением (8).
1
9Дх) — — I МУ(2)1 +
+-
Мк — Мн
±ЇХУ®1 ±У±У®1
1;
1
• х + 9о1; (7)
9і(х)—ЕТ •(М ”(2)‘ •х+
м к — м н
1¥±У®1 1¥±У^> 2
2 • I
х2) + 9 01. (8)
Значение угла поворота в конце 1-го участка (9) определяется выражением (8) при х = 1;:
0,(1,) = ^ ■ (М У,7, ■ 1,+
+-
Мк - Мн
1у±У®1 1у±У®1 \2
2 • 1;
12)+00,. (9)
После упрощения выражения (9) значение угла поворота в конце 1-го участка примет окончательный вид (10):
е1(11)=—■(МУ(7» + МУ(7)1) •11 + 0 (10)
11 2Е I 01
Таким образом, для того чтобы получить значение угла поворота в конце к-го участка необходимо просуммировать выражение (10):
0(к) ^вД); (11)
1=1
0(к)=гТЕ • Е
(Мн + Мк ) ■ 1
V У(7)1 У(7)1/ I
I
.(12)
Закон изменения упругих деформаций можно получить интегрированием выражения (5):
Д ЗУ(7)1 (х) = |81 (х) + Д ЗУ(7)о1 ; (13)
ДЗУ(7)1(х) I
V Е11
(М У(7)1 ■ х +
Мк - М н
У(7)1 У (7)1 2
2 ■ I
х ) + 0 0
+ Д ;(14)
^ ЗУ(7)01 ’ 4 '
ДЗУ(7)1(х) = ■ (
1 - М У(7)1 ■ х
Е1
2
У(7)1 Х,АУ(7)1 3
6 ■:
х ) + 8 01 ■ х + ДЗУ(7)01, (15)
где Д ЗУ(7)1 (х) - значение упругой деформации в сечении с координатой х (0 < х < 11) 1-го участка; ДЗУ(7)01 - значение упругой
деформации в начале 1-го участка.
Таким образом, упругая деформация в конце 1-го участка определяется зависимостью
^ ъцн 12
ДЗУ(7)1 (11 ) т^т
Е1
М
У(7)1 '' 2
У(7)1 У(7)1 ,2
6
+ 801 ■ 1 + ДЗУ(7)01 ; (16)
Д ЗУ(7)1 (11 ) ^
1 (2Мн + Мк ) ■ 12
1 ^1У±У(7)^ У(7)1у Х1
6Е 11
+ ^01 ■ 1 + Д ЗУ(7)01 .
+
(17)
Угол поворота в начале 1-го участка определяется из условия целостности бруса (001 = 81-1) [2]:
801 = 8(1 -1).
(18)
Тогда выражение (17) примет вид:
ДЗУ(7)1 (11)
1 (2МУ + МУ(7)1) ■ 12
6Е
I
+
(
+
— ■ Е
2Е к=1
(Мн + Мк ) ■ 1
У(7)к У(7)к к
I,.
+ 801 ) )
(7)01 •
(19)
Упругая деформация в конце к-го участка определяется как сумма
ДЗУ(7)1 (11 ) , т. е.
Д зу(7) (к) = Е Дзу(7)Л); (20)
1 к
ДЗУ(7)(к) =----------Е
ЗУ(7)Ч/ 6Е 1=1
(2М У(7)1 + М У (7)1) ■ 12
V
к 1
I
+801 ■Е11 + — (Е11 Х
1=1 2Е \1=2
ХЕ
j=1
V
Ij
+Д (21)
Т£АЗУ(7)01- 4 '
))
Для упрощения расчётов удобны следующие замены:
Д
Д
ЗУ(7)1
(2М У(7)1 + МУ(7)1) ■ 12
; (22)
Д6ЗУ(7)1 = (М У(7)1 +1МУ(7)1) ■11 . (23)
Тогда выражение (21) принимает
вид:
1 Л
ДЗУ(7) (к) = ■ Е ДЗУ(7)1 + 801 ■ Е1* +
6— 1=1
1=1
При базировании вала в центрах начальный угол поворота первого участка 801 определяется из выражения (24) и условий ДЗУ(7) (п) =
0, ДЗУ(7)01 0 :
1=1
I
1 П л П
— •Ё АЗУ(Z)1 + 901 •£ 11 +
6Е 1—1
+-
2Е
п 1-1
1—1
л
Ё 1гЕ АЗута — 0; (25)
V1—2 .И1 У
8 - 6- ■ ЕД- - 2-'(&'■ ЕДЧ(26>
801 = п
1=1
Для схемы, приведенной на рис. 5, угол поворота в начале первого участка определяется по формуле (26) при
Д ЗУ(7) (4) “ 0, ДЗУ(7)01 “ 0 :
901 —
6Е
Е • ( АЗУ2)1 + 33
+А3
>У(.)2 ^ЗУ2)3
+ АзУ^ ) '
2е ^ (12 • АЗУ(Z)1 +13 • (АЗУ(Z)1 + АЗУ(Z)2) +
+ АзУ^^
/(І1 +12 + !з +14).
+1 • (А9 + А9
^ Ч lL1ЗУ(Z)Г ЗУ(.2
(27)
ного узла ДШ, а также деформация задней бабки ДБ изменяются в соответствии с изменением силы резания по мере продвижения резца к шпинделю станка. Их величину можно определить исходя из величины действительной жесткости шпиндельного узла ]Ш и задней бабки _|Б соответственно, а также с учетом расстояния Ь0 от торца заготовки до рассматриваемой точки [1]. Принимая во внимание допущение, что их жесткость одинакова в направлении осей У и 7, получаются формулы для расчета величин упругой деформации шпиндельного узла и задней бабки в направлении соответствующих осей У и 2 при обработке заготовки в центрах либо в патроне с поджимом задним центром [1, 3]:
А
ШУ
Р • Т
ГУ О ■
•ш • ТЗ
Р • Т
А — Р. Ь° ; Ц1П7, • у >
•ш • ТЗ
(29)
Значение упругой деформации обрабатываемой заготовки (в конце третьего участка) по формуле (24) при к = 3.
Учет упругой деформации инструмента. В процессе резания на резец воздействует сила резания, которая деформирует его в направлении осей У и 7. Величину деформации инструмента от действия сил резания в направлениях осей У и 7 можно определить по формулам [1]:
А —. РУ • ТИ
ЛИУ
Е и • В • Н
(28)
Р • Т3
А — р Ти ‘—±ТЛГ7
Аиг
ЗЕ • Т
И АИ
где Ьи - вылет резца, м; Еи - модуль упругости для материала державки, Па; В, Н -ширина и высота сечения державки резца соответственно, м; Ти - момент инерции
сечения державки резца, I —
В • Н3
12
Учет деформаций шпиндельного узла и задней бабки. Деформация шпиндель-
А — Ру • (Тз Т°) ;
БУ • т 7
JБ • ТЗ А — Pz•(Lз - Т°)
Б. ‘
(30)
-)б • ТЗ
где _]Ш - жесткость шпиндельного узла, Н/м; _]Б - жесткость задней бабки, Н/м.
В результате совместного действия деформаций шпиндельного узла ДШ и задней бабки ДБ суммарная величина деформации на длине Ь0 может быть определена из рис. 4 по формуле
А Т° - А А - А
ШБУ® ШУ® БУ® ШУ®
Т - Т
З О
А Т° — А
ШБУ ШУ
А Т° — А
шб. тпу.
Т
(31)
(АБУ АШУ) • (ТЗ ТО) ;
Тз
(АБ. - АШ. ) • (ТЗ - ТО )
(32)
и
При обработке заготовки в патроне учитывается только упругая деформация шпиндельного узла, которая не зависит
от положения рассматриваемой точки на поверхности заготовки. Её величина изменяется лишь в зависимости от действующей силы резания. Расчетные формулы имеют следующий вид [1]:
Д =
ШУ • 5
-)ш р.
(33)
ДШ7
— 7
.1ш
Рис. 4. Схема определения суммарного прогиба задней бабки и шпиндельного узла
Учет погрешности формы и смещения центровых отверстий заготовки. Овальность ДО И конусность ДК поверхности заготовки оказывают влияние на исходный радиус заготовки Я3, соответствующий расстоянию от оси заготовки до точки, расположенной на её поверхности, и не зависят от изменения силы резания [1]. Наличие исходной овальности поверхности До оказывает влияние на изменение радиуса заготовки от нуля до ДО во время вращения заготовки вокруг своей оси. Величина конусности ДК изменяет его величину от Я3 до Я3 - ДК. Можно записать формулу для расчета величины исходного радиуса заготовки Я3 для каждой точки её поверхности, выраженного через угол ф полярной системы координат [1, 3]:
где Я3 - исходный радиус заготовки, м; ДК - величина конусности поверхности заготовки, м; ф - полярная координата точки, расположенной на поверхности заготовки, град; ДО - величина овальности поверхности заготовки, м.
Таким образом, величина исходного радиуса поверхности заготовки Я3 для каждого поперечного сечения изменяется по мере перемещения резца от правого торца заготовки на длину ЬО на величину ДК • ЬО / ЬЗ [1].
Кроме этого, наличие смещения центровочных отверстий ДСп, ДСЛ и погрешности базирования заготовки ДП приводит к дополнительному смещению её оси относительно оси центров станка. Эти факторы также не зависят от силы резания и в данной модели учтены при нахождении действительного радиуса заготовки Яд , определяемого как расстояние от оси центров станка до точки на поверхности заготовки [1]. Таким образом, действительный радиус заготовки Яд с учетом смещения центровых отверстий заготовки при обработке в центрах можно определить по формуле [3]
+ АСМ • С^(Ф + Юl),
(35)
где ДСМ - суммарная величина смещения в рассматриваемой точке, расположенной на длине ЬО, м; ю1 - угол между большой полуосью эллипса и прямой, проведенной из центра заготовки через правое центровочное отверстие [1], град.
Суммарная величина смещения ДСМ определяется как геометрическая сумма смещений центровых отверстий правого и левого торцов заготовки на длине ЬО [1]:
АСМ
А2 + А2 - у
^ СМП ^ ^ СМЛ ^^СМП л (36)
Х АСМЛ • СОБ
где ДСМП, ДСМЛ - величины смещения центровых отверстий правого и левого торца заготовки на длине ЬО соответственно [1], м,
АСМП АСП
АСМЛ АСЛ
с и ^ 1 - Т
V ТЗ У
Ґь л _Ч°
V ТЗ У
(37)
(38)
где ДСП, ДСЛ - величины смещения центровых отверстий правого и левого торца заготовки соответственно, м; фцц - угол между прямыми линиями, проведенными от оси заготовки через центровые отверстия, град.
При обработке в патроне с поджимом задним центром расчет также ведется по формуле (35), только суммарная величина смещения ДСМ определяется как геометрическая сумма смещения центровочного отверстия правого торца заготовки ДСП и смещения, вызванного погрешностью базирования заготовки в патроне ДП на длине ЬО [1]:
АСМ
А 2
СМП
2А,
с ь у
А •
П Ьз,
• А, • То
(39)
СМП П
СОБ
(фдц).
При обработке в патроне действительный радиус заготовки можно определить по формуле [1]
ЯД — -у/яЗ - (АП • б1п(ф + ю2))2 + АП • соб(ф + ю2),
(40)
где ю2 - угол между большой полуосью эллипса и вектором смещения заготовки при базировании в патроне станка, град.
Учет смещения оси задней бабки. При обработке заготовки в центрах либо в патроне с поджимом задним центром на действительный радиус детали Я Д будет непосредственно влиять смещение оси задней бабки АСБУ и Асб. [1].
Эту величину смещения оси задней бабки по направлению осей У и . для любой точки заготовки, расположенной на длине Ь° от правого торца, можно определить из соотношений [1]:
А —А •1—
СБУЬ СБУ ’
Асбуь Асбу •
З
1 - и
(41)
О
З
где ДСБУ, ДСБ7 - смещение оси задней бабки станка, измеренное в направлениях осей У и 7 соответственно, м.
Учет исходной шероховатости поверхности. Для определения изменения глубины резания, вызванного исходной шероховатостью поверхности заготовки Д^, введем допущение, что высота неровностей изменяется случайным образом в соответствии с нормальным законом распределения. Тогда, в соответствии с правилом «шести сигм», для каждой точки поверхности заготовки можно принять значение шероховатости поверхности соответственно формулам [4]:
З
Rz
(42)
(43)
где о^ - среднеквадратическое отклонение высоты неровностей поверхности, мкм; К^3 - значение шероховатости поверхности, мкм; ЯапёБ - случайная функция, значения которой соответствуют стандартизованному нормальному распределению.
После определения значений деформации всех элементов технологической системы определяется действительный радиус детали по формуле [1, 5]
R Д =
*-Д -
ЗZ 1 + АИZ
2
Абу ) + (44) AБz)-.
Математическая модель деформаций элементов технологической системы позволяет определить, какое влияние на точность обрабатываемых поверхностей оказывают группа, все или один из следующих параметров: отклонение от круг-лости (овальность) и погрешности профиля продольного сечения (конусность) исходной заготовки; смещение центровочных гнезд заготовки (при обработке в центрах) и смещение оси задней бабки;
исходная шероховатость поверхности заготовки и погрешность её базирования в патроне; податливость шпиндельного узла и задней бабки, жёсткость инструмента; упругая деформация заготовки. А также данная модель может применяться в САПР, направленных на повышение качества механической обработки ступенчатых валов при использовании станков с ЧПУ.
СПИСОК ЛИТЕРАТУРЫ
1. Жолобов, А. А. Прогнозирование и обеспечение качества технологических систем на этапах их проектирования и изготовления : монография / А. А. Жолобов. - Могилев : Бело-рус.-Рос. ун-т, 2005. - 304 с. : ил.
2. Биргер, И. А. Сопротивление материалов : учеб. пособие / И. А. Биргер, Р. Р. Мавлю-тов. - М. : Наука, 1986. - 224 с.
3. Жолобов, А. А. Математическая модель для прогнозирования деформации заготовки при точении / А. А. Жолобов, А. М. Будкевич // Машиностроение : сб. науч. тр. - Минск : Технопринт, 2003. - № 19. - С. 310-314.
4. Жолобов, А. А. Прогнозирование шероховатости поверхности при точении / А. А. Жолобов, А. М. Будкевич // Машиностроение : респ. межведомств. сб. науч. тр. - Минск : Технопринт, 2002. - № 18. - С. 23-25.
5. Жолобов, А. А. Застосування 3Б модел1 ГВМ для тдвещення продуктивност1 та якост1 механ1чно1 обробки / А. А. Жолобов, А. В. Казаков // Весн. ЖДТУ. - 2006. - № 4 (39). -С. 21-25.
°Rz _
6
Белорусско-Российский университет.
Материал поступил 16.10.2008
A. A. Zholobov, A. V. Kazakov Features of calculation of elastic deformations of technological system at turning processing of a step shaft
The problems of estimation of errors of processing from deformation of technological system at turning of a step shaft on the machine tool with numerical program control are considered. The approach based on the method of professor V. S. Korsakov is used. The mathematical model allowing to reveal and estimate the influence of deformation of a half-finished product, a tool, rear and front mandrels of the machine tool on accuracy and form of processing surfaces of a step shaft is developed.