МАШИНОСТРОЕНИЕ
УДК 621.9
А. В. Казаков, А. А. Жолобов
ПРОГНОЗИРОВАНИЕ И ОБЕСПЕЧЕНИЕ ТОЧНОСТИ ТОКАРНОЙ ОБРАБОТКИ СТУПЕНЧАТЫХ ВАЛОВ
UDC 621.9
A. V. Kazakov, A. A. Zholobov
FORECASTING AND ENSURING THE ACCURACY OF STEPPED SHAFTS TURNING
Аннотация
Рассмотрены вопросы оценки и прогнозирования погрешностей обработки, возникающих вследствие деформации элементов ТСР, влияния внешних силовых факторов на упругую деформацию элементов ТСР при токарной обработке ступенчатых валов. Разработана математическая модель, которая позволила выявить и оценить влияние деформации заготовки, инструмента, задней и передней бабок станка на точность и форму обрабатываемых поверхностей ступенчатого вала.
Ключевые слова:
деформация, силы резания, технологическая система резания, токарная обработка, повышение точности, оптимизация, математическая модель, ступенчатые валы, прогнозирование, коррекция траектории движения режущего инструмента.
Abstract
The paper considers the problems of evaluating and forecasting processing errors caused by the deformation of elements of a technological cutting system (TCS) and the influence of external forces on elastic deformation of TCS elements when turning stepped shafts. The mathematical model has been developed, which reveals and assesses the impact of the deformation of a workpiece, of a cutting tool, and that of a headstock and a tailstock of the machine, on the accuracy and form of processed surfaces of a stepped shaft.
Key words:
deformation, cutting forces, technological cutting system (TCS), turning, increased accuracy, optimization, mathematical model, stepped shafts, forecasting, adjustment of the contour of the cutting tool movement.
При токарной обработке валов существенное влияние на точность формируемых поверхностей оказывает геометрическая точность технологической системы резания (ТСР) и податливость её элементов под воздействием сил резания [1-3]. Для обеспечения заданных конструктором параметров
© Казаков А. 1В ., Жолобов А. А., 2014
изделия необходим строгий учёт возникающих и изменяющихся в процессе обработки погрешностей.
Многими учеными проведены теоретические и экспериментальные исследования этих погрешностей, в том числе и в Белорусско-Российском университете [3-6]. Однако в этих работах не
уделялось особого внимания специфике формирования погрешностей на поверхностях ступенчатых валов.
Для прогнозирования точности токарной обработки ступенчатых валов авторами сделана попытка создания математической модели, позволяющей выявить и оценить влияние деформации заготовки, инструмента, задней и передней бабок станка на точность размеров и форму обрабатываемых поверхностей ступенчатого вала.
В математической модели были учтены рекомендуемые в ранее выполненных работах [3, 4, 7-12] параметры: пространственные отклонения заготовки, погрешность установки заготовки, смещение центровых отверстий заготовки (при обработке в центрах); погрешность базирования заготовки в патроне; смещение оси задней бабки; податливость шпиндельного узла; податливость задней бабки; жесткость инструмента; упругая деформация заготовки.
Кроме того, при разработке модели было принято допущение, что материал поверхности заготовки имеет однородную структуру и равномерную твердость [3].
Рассматриваемая математическая модель имеет отличие от моделей, предлагаемых в [3, 4, 7-12]: применяется дифференциальное уравнение изогнутой оси балки для расчёта упругих деформаций, возникающих в технологической системе. Подход направлен на более полное раскрытие сущности моделирования процессов механической (токарной) обработки ступенчатых валов на станках с ЧПУ.
Основными причинами возникновения погрешностей цилиндрических поверхностей, связанных с податливостью ТСР, является влияние упругой деформации заготовки и инструмента и упругой деформации шпиндельного узла и задней бабки. Указанные факторы в разной степени сказываются на погрешностях обрабатываемой поверхности [3].
Для установления влияния различных внешних силовых факторов на упругие деформации элементов ТСР и, как следствие, на погрешности формируемых в процессе обработки цилиндрических поверхностей ниже рассмотрено общее решение дифференциального уравнения изогнутой оси балки.
Прогиб (величину упругой деформации) заготовки с учётом упругой деформации шпиндельного узла и задней бабки можно определить с помощью дифференциального уравнения изогнутой оси балки [13]
Д
ЗУ ( 2 )
( *)''=■
Мх (У)( х)
Е1 ( х)
(1)
где ДЗУ(2)(х) - уравнение упругих деформаций заготовки в плоскости УОХ (2ОХ), м; М(х) - уравнение изгибающего момента, вызванного силой резания в рассматриваемой плоскости, Н-м; Е1(х) - уравнение изменения жёсткости при изгибе заготовки, Н-м2.
Уравнение (1) описывает закон изменения упругой деформации по всей длине обрабатываемой заготовки. Значение же упругой деформации, определяющее точность цилиндрической поверхности после механической обработки, будет находиться в точке приложения силы резания. То есть в любой момент времени в точке касания вершины резца и обрабатываемой поверхности отмечается такое значение упругой деформации заготовки, которое впоследствии оказывает влияние на точность обработанной поверхности. Для этого уравнение (1) приводится к виду, позволяющему определять значение упругой деформации в заданной точке.
Необходимо решить дифференциальное уравнение (1). Следует зафиксировать резец в произвольный момент обработки, а его действие на обрабатываемую заготовку заменить соответствующими составляющими силы резания. Таким образом, обрабатываемую заготовку можно представить с определёнными
допущениями как набор участков, жёсткость которых постоянна [14].
В процессе точения на заготовку действуют силовые факторы в виде сосредоточенных сил и изгибающих моментов, а также собственный вес заготовки как равномерно распределённая нагрузка по каждой из ступеней.
В качестве примера на рис. 1 представлен вал, содержащий N участков (ступеней) в зафиксированный момент обработки. Расстояние от правого торца заготовки до текущего положения резца обозначено через Ь0, а вылет заготовки - через ЬЗ. Пусть на заготовку действует силовой фактор £ (изгибающий момент или поперечная сила), при-
ложенный в начале участка с номером пР (см. рис. 1). Необходимо определить, какое влияние £ оказывает на деформацию в конце участка с номером п.
Закон изгибающего момента Мг(у)(х) в пределах участка с номером п от силового фактора £, приложенного в начале участка с номером пР, выразим следующей зависимостью:
М2 (У )( Х)
= т2 (У )(£ ПР , п) '
(2)
где /п - длина участка с номером п, м.
0
Рис. 1. Схема для определения изгибающего момента от силовых факторов
Интегрирование выражения (1) в пределах от 0 до /п даёт значение угла поворота в конце участка с номером п
[13]:
1}МУ (2 )(х)
в2(У) (П) = Д32(У) (П)'= I -+ вг(у)0п , (3)
0 Е1п
где в2(У )(п) - угол поворота в конце п-го участка, рад; в2 (У )0п - угол поворота в начале п-го участка, рад.
Угол поворота в начале п-го участка определяется из условия целостно-
сти ступенчатого вала [13]:
в2 (У )0п =в2 (У )(п -1) . (4)
Исходя из вышеизложенного, угол поворота в2(У)(£,пР,п) в конце участка с номером п от силового фактора £, приложенного в начале участка с номером пР:
в2 (У )(£ пР, п ) =
= ] ш-+ в2(У )(£ пр, п -1). (5)
0 Е1п
После суммирования выражения начального угла поворота в2{Т)(£,пр,п-1') от пр до п - 1
02(I )(£ пр, п -1) =
п-1 ¡1
= Т1
1=пр 0
т1 (2)(£ пР, 1)
ЕГ
+ 02 (I )01. (6)
С учётом (5) выражение (6) принимает следующий вид:
02 (I )(£ пр, п ) =
п '
Т I
1=пр 0
4 т1 (2 )(& пР, 1)
ЕГ
+ 02 (I )01, (7)
где в2(у)01 - угол поворота в начале первого участка, рад.
Значение упругой деформации в конце участка с номером п можно получить интегрированием выражения (3)
[13]:
лз2 (i)( п) =
„п
=11
MI (2)( * )
ЕГ
+ ¡п ' 02(I)0п + Л2(I)0п, (8)
где АЗГ(2)(п) - значение упругой деформации в конце п-го участка, м; АЗТ(2)0п - значение упругой деформации в начале п-го участка, м.
Значение упругой деформации в начале п-го участка определяется из условия целостности ступенчатого вала [13]:
АЗ2(I)0п = АЗ2(I)(п - 1). (9)
Исходя из (2), (8), (9) значение упругой деформации Лг(Т)(£,пр,п) в конце участка с номером п от силового фактора приложенного в начале участка с номером пр, определяется зависимостью
АЗ2(I )(£ пР, п) = I ^ (2)(£ пР, п)
0 Е1п
+ 1п 'в2(I)(£ пР ,п - 1) +
+ АЗ2(I)(£ пр, п -1). (10)
После суммирования выражения упругой деформации ЛЗщ)(£,пр,п - 1) от пр до п - 1
АЗ2(I)(£ пр, п - 1) =
п-1 ¡к
= Т 11
п-1 ¡к ^ (2 )(£ пр, 1)
к=пр 0
Е1
+
+ ¡к -в(£, пр, к -1) + АЗ2 (I )01, (11)
где ЛЗ2(Т)01 - упругая деформация в начале первого участка (при пр = 1), м.
С учётом (7) выражение (11) принимает следующий вид:
АЗ2(I )(£ пр, п - 1) =
п-1
Т
к=пр
¡1 гЩ (2 )(£ пр, к) + 1 11 Е1 + к
( к-1 4
ТI
т7 (2 )(#, пр, 1)
2 (I )01
/
+ А
З2 (I )01
(12)
С учётом (10) окончательно
АЗ2 (I )(£ пр,п) =
к=пр
¡к 11
т (2 пр>к)
Е1к
+¡кх
'Т I mI(2)(£ ^1)
У1="р 0
I (2 'р
Е
2 (I )01
У
+А,
З2 (I )01
.(13)
Таким образом, совместное влияние различных силовых факторов на деформацию заготовки сводится к определению функции, описывающей уравнение изгибающего момента в пределах рассматриваемого участка (14) и двух интегралов (15) и (16):
0
0
0
т (z nP, п);
lk i i
mY(Z)(£ ПР , k) .
EI„
fmY(Z)(£ ПР , 0
i
EL
(14)
(15)
(16)
Внешние силовые факторы, действующие на обрабатываемую заготовку, можно разделить на три группы. Первая -сосредоточенный изгибающий момент. В качестве примера может служить изгибающий момент, вызванный силой Рх.
Вторая - сосредоточенная поперечная сила (силы резания Ру, Рг). Третья - собственный вес детали, который необходимо учитывать при обработке длинных и нежёстких валов (в виде равномерно распределённой по каждой ступени, нагрузки).
Далее выражения (14)...(16) преобразовываются для каждой группы силовых факторов.
Для постоянного момента £ = М (рис. 2) выражение (14) принимает вид:
mM,M
Y (Z )V
Y(Z)'' "Р
, Пр, n) = My (Z). (17)
о
Рис. 2. Закон изменения изгибающего момента при М = const
В (15) и (16) подставлено (17), получены следующие выражения:
\ mYM(Z) (MY(Z) , ПР , 0 = MY(Z)h n
J-™-- ™ ; (18)
EL
EL
h mMz)(My(z),Пр,k) = MY(Z$ .
ii
EIk 2 • EIk
(19)
В (13) подставлены (18) и (19), получена зависимость, выражающая влияние постоянного изгибающего момента на деформацию заготовки:
A MZ (y )( MY ( Z ^ ПР, П ) =
M
Y ( Z )
E
Z 4
k—пр
1 k-1 1 —+ Z -
2I ^ I
k i—nP i
. (2о)
Для поперечной силы (£ = Р) (рис. 3) зависимость (14) принимает вид:
( п-1
тр( Z )(PZ (Y Ъ ПР, п) = PZ
Z (Y )
Л
Z h+*
V J—nP
(21)
о
Рис. 3. Закон изменения изгибающего момента от поперечной силы р
Далее последовательно подставлены (21) в (15) и (16), (22) и (23) в (13), получена зависимость (24), выражающая влияние поперечной силы на деформацию заготовки:
(2)(р2(I), пр , 1) р2(I )11
I (2 ) V 2 (I )>' *р
Е
Е1
1-1 I
т ¡,+2
\]="р
Л
;(22)
к (2)(р (I), пр,к)
Е1ь
р 12 ( к-1
2 (I у к
2 • ЕГ
к
Т Ь + |
V )=«р
(23)
р.
АPI(2)(р7(2),пр ,п) = ^ ) Т ¡к Х
Е
к=пр
(
2 А
к-1 7
Т ¡, + ^
V}=пр 3
Л
А
у
Л"
к-11 1 -1 ; + Т } Т ¡> + £
1 = пр 11 V ] = пр У
(24)
Собственный вес заготовки можно рассматривать как равномерно распределённую нагрузку по каждой ступени [13], приложенную от первого до последнего участка (рис. 4).
0
Рис. 4. Закон изменения изгибающего момента от собственного веса
Поэтому зависимость (14) принимает вид:
(
т
2(У)
(п) = Ц ] + I/
1=1
I
V ]+1
2
. (25)
Далее последовательно подставлены (25) в (15) и (16), (26) и (27) в (13), получена зависимость (28), выражающая влияние собственного веса заготовки на её деформацию:
Н(У)(/) / г^ .
I У = ЕТ I ]
0 Е1 Е1 ^ 1=1
(
/ '-1 / 2 + и +2
V2 *=1+1 2
У
(26)
У
•т2 (у )(к) = /
Л—*
0
А,. , /
Е1к
2 Е1к
к-1
I ] ] +1: ¡.+¡3
V1=1 V 2 *=]+1 3
Л ,2 л + Чкк
12
У
;(27)
у
Д 3 (п) = 1 4
к=1
2ЕА
к-1
Г ' к-1 '
I ] + 11 '« + Т
V]
=1
V №=]+1
^ /2 Л
+ чА 12
к-1 7 Г'-1 Г / '-1 7
^Ег I: Ч]/] к +1 +£
г=1 ^ г V ]=1 V 2 1+1 2
6
у У
.(28)
Таким образом, с учётом вышеизложенных зависимостей и преобразований уравнение (1) принимает следующий вид:
д(п)зу (г) = Д ЗУ (г )01 +в01У (г) I /г +
г=1
+ !Д М (2 )(Мг2 (У), пРг, п ) +
+ Д Р (2 )(Ру (2), прг, п) + Д 3 (п ). (29)
Это уравнение описывает значение упругой деформации в конце участка с номером п.
Влияние податливости шпиндельного узла АШУ(2), задней бабки Абщ), инструмента АИУ(2), а также смещения оси задней бабки Асбщ) и геометрических неточностей ТСР на деформацию заготовки подробно рассмотрено в [14]. Эти факторы могут быть учтены посредством задания соответствующих начальных параметров в выражении (29).
На основе уравнения (29), методик структурно-параметрической оптимизации [15] и коррекции траектории движения режущего инструмента [16] разработана интегрированная система автоматизированного проектирования ступенчатых валов и управляющих программ ЧПУ, позволяющая управлять точностью токарной обработки на этапе технологической подготовки производства.
Для проверки адекватности и значимости математической модели (29) проводилась обработка партии валов на токарном станке модели 16К20Ф3, оснащенном системой ЧПУ N0-201 с последующим изменением на координат-но-измерительной машине БигаМах. Исходные данные эксперимента: диаметры обрабатываемых ступеней А = 20, = 25, А3 = 30 мм; длины обрабатываемых ступеней Ь1 = Ь2 = Ь3 = 40 мм; диаметр заготовки Пзаг = 38 мм; длина заготовки ЬЗ = 240 мм; предел прочности материала заготовки ав = 610 МПа; модуль упругости материала заготовки Е = 2,1 • 1011 Па.
Параметры инструмента: материал режущей пластины - Т15К6; передний угол у = 0°; главный угол в плане Ф = 90°; радиус при вершине г = 2 мм.
Результаты измерений приведены на рис. 5. В табл. 1 представлены сравнительные данные измерения погрешности профиля продольного сечения ступеней вала.
Применение интегрированной системы автоматизированного проектирования ступенчатых валов и программ
п—1
п—1
Ч
X
токарной обработки, основанных на управлении траекторией формообразующих элементов, позволило повысить точность обработки деталей в условиях эксперимента на 20...40 %. Данная сис-
тема может быть использована на промышленных предприятиях для прогнозирования и обеспечения точности изготовления ступенчатых валов.
Рис. 5. Погрешности профиля продольного сечения ступеней вала: а-в - построенные на основе формулы (29); г-е - построенные на основе измерений после обработки без коррекции траектории движения режущего инструмента; ж-и - построенные на основе измерений после обработки с коррекцией траектории движения режущего инструмента
Табл. 1. Результаты измерения
Диаметр, мм Погрешность профиля продольного сечения
Наименование без коррекции, мкм с коррекцией, мкм
20 Конусность 100,4 61,2
25 57,1 41
30 41 32,5
СПИСОК ЛИТЕРАТУРЫ
1. Жолобов, А. А. К вопросу технологического обеспечения качества поверхностей деталей при их формировании / А. А. Жолобов, С. В. Гонорова, В. А. Широченко // Отделочно-упрочняющая технология в машиностроении : сб. тез. докл. Междунар. науч.-техн. конф. - Могилев : Белорус.-Рос. ун-т, 1994. - Ч. 1. - С. 67-69.
2. Жолобов, А. А. Прогнозирование и контроль возможностей технологической системы СПИД / А. А. Жолобов // Теоретические и технологические основы упрочнения и восстановления изделий машиностроения : сб. науч. тр. - Полоцк : Технопринт ; Полоцкий гос. ун-т, 2001. - С. 78-79.
3. Жолобов, А. А. Прогнозирование и обеспечение качества технологических систем на этапах их проектирования и изготовления : монография / А. А. Жолобов. - Могилев : Белорус.-Рос. ун-т, 2005. - 304 с. : ил.
4. Будкевич, А. М. Автоматизация прогнозирования деформации технологической системы при точении / А. М. Будкевич, К. А. Шалыжин ; науч. рук. А. А. Жолобов, канд. техн. наук, доц. // Новые материалы, оборудование и технологии в промышленности : материалы Респ. науч.-техн. конф., Могилев, 29 янв. 2004 г. / М-во образования Респ. Беларусь [и др.] ; редкол. : И. С. Сазонов [и др.]. - Могилев, 2004. - С. 22.
5. Пашкевич, В. М. Повышение точности механической обработки и сборки изделий машиностроительного производства на основе построения и использования компьютерных обучающихся систем : автореф. дис. ... д-ра техн. наук : 05.02.08 / В. М. Пашкевич ; Белорус. нац. техн. ун-т. - Минск, 2007. - 45 с.
6. Пашкевич, В. М. Самообучающиеся системы искусственного интеллекта в машиностроении /
B. М. Пашкевич, Ж. А. Мрочек. - Могилев : МГТУ, 2003. - 434 с.
7. Жолобов, А. А. К вопросу прогнозирования и контроля геометрической точности технологических систем / А. А. Жолобов // Перспективные технологии, материалы и системы : сб. науч. тр. -Могилев : МГТУ, 2001. - С. 23.
8. Жолобов, А. А. Математическая модель для прогнозирования деформации заготовки при точении / А. А. Жолобов, А. М. Будкевич // Машиностроение : сб. науч. тр. - Минск : Технопринт, 2003. - № 19. - С. 310-314.
9. Жолобов, А. А. Прогнозирование точности деталей при обработке на токарных станках / А. А. Жолобов, А. М. Будкевич // Теория и практика машиностроения. - 2004. - № 1. - С. 43-46.
10. Жолобов, А. А. Автоматизация процесса прогнозирования качества поверхностей деталей /
A. А. Жолобов, А. М. Будкевич // Вестн. МГТУ. - 2002. - № 2. - С. 11-14.
11. Жолобов, А. А. Методика контроля геометрической точности технологического оборудования и деталей / А. А. Жолобов, А. М. Будкевич // Современные методы и приборы контроля качества и диагностики состояния объектов : материалы Междунар. науч.-техн. конф. - Могилев, 2004. -
C. 160-161.
12. Жолобов, А. А. Прогнозирование шероховатости поверхности при точении / А. А. Жоло-бов, А. М. Будкевич // Машиностроение : Респ. межведомств. сб. науч. тр. - Минск : Технопринт, 2002. - № 18. - С. 23-25.
13. Писаренко, Г. С. Справочник по сопротивлению материалов / Г. С. Писаренко, А. П. Яковлев,
B. В. Матвеев ; отв. ред. Г. С. Писаренко. - 2-е изд., перераб. и доп. - Киев : Наук, думка, 1988. - 736 с.
14. Жолобов, А. А. Особенности расчёта упругих деформаций технологической системы при токарной обработке ступенчатого вала / А. А. Жолобов, А. В. Казаков // Вестн. Белорус.-Рос. ун-та. -2009. - № 3. - С. 87-95.
15. Казаков, А. В. Оптимизация режимов резания технологических переходов наружной токарной обработки ступенчатых валов на станке с ЧПУ / А. В. Казаков, А. А. Жолобов // Вестн. Белорус.-Рос. унта. - 2012. - № 4. - С. 23-31.
16. Жолобов, А. А. Коррекция траектории движения режущего инструмента на токарных станках с ЧПУ / А. А. Жолобов, А. В. Казаков // Вестн. Полоцкого гос. ун-та. - 2009. - № 9 (серия С). -
C. 124-133.
Статья сдана в редакцию 17 января 2014 года
Александр Алексеевич Жолобов, канд. техн. наук, проф., Белорусско-Российский университет. Тел.: 8-0222-23-04-51.
Алексей Владимирович Казаков, ассистент, Белорусско-Российский университет. Тел.: 8-0222-23-04-51.
Aleksandr Alekseyevich Zholobov, PhD (Engineering), Prof., Belarusian-Russian University. Phone: 8-0222-23-04-51.
Aleksei Vladimirovich Kazakov, assistant lecturer, Belarusian-Russian University. Phone: 8-0222-23-04-51.