№4
2008
Таким образом, сложная шестизвенная группа Ассура с четырехугольным замкнутым изменяемым контуром (6ГА6) кинематически и кинетостатически оказывается вполне разрешимой.
СПИСОК ЛИТЕРАТУРЫ
1. Э.Е. Пейсах. Атлас структурных схем восьмизвенных плоских шарнирных механизмов // Теория Механизмов и Машин, №1, т. 4, 2006. - с. 3-17.
2. Дворников Л, Т. О кинематической разрешимости плоской четырехзвенной группы Ассура четвертого класса графо - аналитическим методом, Известия ВУЗов, «Машиностроение», №12, 2004. - с 9-15.
534.11 ; 621.923.5
ДИНАМИКА ПРОЦЕССА ХОНИНГОВАНИЯ И ЕЕ ВЛИЯНИЕ НА ПОГРЕШНОСТИ ФОРМЫ ПРИ ОБРАБОТКЕ ОТВЕРСТИЙ
Канд.техн.паук, доц. СА.ВОРОНОВ, инж.ОЛБОБРЁНКОВ
Разработана модель процесса хонингования отверстий с заданной первоначальной погрешностью формы поверхности. Численно анализируется процесс исправления погрешности формы отверстий в процессе обработки. Вал инструмента схематизируется как вращающийся безмассовый стержень, передающий вращение оправке хона с упруго закрепленными на нем режущими и направляющими брусками. Силы резания моделируются как силы трения, зависящие от нормального давления на контактной поверхности, снятие припуска при прохоэюдении режущих брусков поверхности моделируется согласно гипотезе Престона. При движении брусков вследствие колебаний вала и оправки силы резания меняются, что может привести как к исправлению погрешностей формы, так и к ее ухудшению. Проводится анализ влияния технологических параметров на интегральные критерии качества обработки.
The sample piece of process of a honing of apertures with the set initial error of the molding box of a surface is developed. Process of correction of an error of the molding box of apertures during treatment is numerically analyzed. The shaft of the tool is schematized as rotating безмассовый a core transferringftrans mittingj rotation to a mandrel хона with is elastic fixed on it{him} cutting and directing брусками. Forces of cutting are modelled as the friction forces depending on normal pressure on a contact surface. Removal of an allowance at passage cutting bars surfaces is modelled according to a hypothesis of Preston. At bars movement owing to fluctuations of a shaft and a mandrel of force of cutting vary, that
can result in correction of errors of the molding box, and to its deterioration. The assayring of influence of technological parameters on integrated criteria of quality of treatment is spent.
В настоящее время в серийном и массовом производстве при изготовлении ответственных деталей предъявляются высокие требования к точности и шероховатости поверхности: некруглость менее 1 мкм, волнистость менее 0,2 мкм, нецилиндричность и непрямолинейность образующей менее 2-5 мкм, параметр шероховатости Ra = 0,02-0,8 мкм, отсутствие дефектного слоя металла (структурно-фазовых изменений, напряжений растяжения, микротрещин), определенные значения параметров формы микронеровностей и опорной поверхности [1].
Обеспечение этих требований достигается с помощью таких процессов абразивной обработки, как хонингование брусками из традиционных и сверхтвердых абразивных материалов. Этот процесс относят к процессу доводки; хонингование производится при одновременно выполняемых вращательном и возвратно-поступательном движениях инструмента (головки с брусками). На рис.1 приведена схема рабочего движения инструмента. Подача (разжим) брусков в радиальном направлении при хонинговании производится либо непрерывно, под воздействием постоянного усилия, либо периодически, на каждый двойной ход хонинговальной головки.
При контакте рабочей поверхности бруска с обрабатываемой поверхностью заготовки происходит царапание металла одновременно большим числом абразивных частиц. Размер таких частиц при хонинговании составляет 20—100 мкм, среднее число частиц на поверхности бруска 20—400 зерен на 1 мм2[1]. Основными видами взаимодействия абразивных зерен с металлом являются микрорезание со снятием тончайших стружек и трение с пластическим оттеснением металла. Для интенсивного резания необходимо, чтобы абразивный брусок самозатачивался путем скалывания и вырывания затупившихся зерен из связки. При использовании брусков из сверхтвердых абразивных материалов (алмаза, эльбора) зерна длительное время сохраняют свою остроту, преобладает микроскалывание зерен, а не вырывание их, что значительно повышает стойкость брусков.
Путем выбора оптимальных характеристик брусков и регулирования параметров обработки (скорости, давления) можно управлять процессом обработки, осуществляя на первой стадии непрекращающееся резание металла в течение достаточно длительного времени, необходимого для исправления погрешностей формы заготовки, удаления исходной шероховатости и дефектного слоя. Скорость съема металла при этом составляет 2-4 мкм/с. Для получения поверхности с малой шероховатостью (Ra = 0,1—0,3 мкм при хонинговании), а также для создания благоприятного микрорельефа поверхности деталей и упрочненного поверхностного слоя металла процесс обработки на заключительной стадии
№4
2008
может быть переведен в режим выхаживания, при котором съем металла резко сокращается, а брусок выглаживает обрабатываемую поверхность. Такой переход можно осуществить, изменяя параметры обработки: повышая окружную скорость заготовки или инструмента, снижая давление бруска и частоту
Уок
Рис. 1. Схема движений бруска и детали при хонннговании: 1 —деталь; 2 — брусок: 3 — перебег;
4 — перекрытие (Уок — окружная скорость, Увп — скорость возвратно-поступательного движения,
а. — угол сетки)
При хонинговании в системе поддержания контакта бруска с деталью контакт замыкается кинематически, с помощью клиновой пары (рис. 2) или гидравлической подачи, жесткость системы прижима брусков высока, сила резания при наличии погрешностей формы заготовки непрерывно изменяется. Имеется ряд современных хонинговальных станков, на которых заготовка или хонинговальная головка также совершают дополнительное колебательное движение; такой процесс назван вибрационным хо-нингованием. Этот процесс особенно эффективен при хонинговании глухих отверстий.
О пряв к н
Хонши-овальиыЯ брусок
Деталь
Разжимающий конус
Рис. 2. Схема инструмента для хошнгования
№4
2008
Хонинговалытый инструмент может иметь как один (рис. 3)5 тате и несколько режущих брусков (рис. 4, 5).
На погрешности формы, имеющиеся в обрабатываемой детали, главным образом воздействует контактная сила между нею и бруском хона. Силы хонингования - технологический параметр, который неуправляем, т.к. он зависит от многих других переменных,: параметров инструмента, вращения обрабатываемой детали и скорости возвратно-поступательного движения оправки. Силы хонингования могут вызывать вибрацию и теплонапряженность в процессе хонингования. Данные факторы влияют на точность размеров — например, на отклонение отверстия от круглости.
В настоящее время не разработано полной модели процесса хонингования глубоких отверстий с учетом влияния динамики инструмента на исправление погрешности формы цилиндрических поверхностей. В имеющихся публикациях рассматриваются в основном технологические аспекты процесса, отдельно затрагиваются вопросы динамики вращающихся валов, а также рассматриваются различные модели формирования новых поверхностей при резании.
Деталь
Деталь
Режущий брусок
Брусок
Направл. бруски
Рис. 3. Схема инструмента с одним режущим и двумя направляющими брусками
Рис. 4. Схема мульгибрускового инструмента
Рис. 5. Фотографии хонинговальной головки
В [1, 2] описываются типовые режимы обработки, конструкции станков и хонинговального инструмента, а также пути повышения эффективности обработки отверстий с использованием подходов, типичных для технологов, то есть путем варьирования технологических параметров рассматриваемого процесса, применения различных материалов для режущего инструмента и т.д. в сочетании с проведением большого числа экспериментов.
Исследованиям различных динамических моделей процесса хонингования посвящены [3] - [7]. В [3] приводятся результаты исследования динамических эффектов, сопровождающих процесс виброабразивной обработки на основе разработанных математических моделей поведения инструмента и детали, с учетом вращения инструмента. Подсчитаны собственные частоты колебаний и определены области динамической неустойчивости при гармонических осевой составляющей силы резания и крутящем моменте.
В работе Грязнова А.Н. и Гуськова A.M. [4] рассматривается плоская математическая модель формообразования поверхности при обработке абразивными брусками. Приводятся формулы для расчета и прогнозирования погрешностей формы при хонинговании и суперфинишировании.
Обобщенная модель формирования поверхности при резании представлена в [5]. Разработан новый метод представления процесса формирования новых поверхностей, включающий в себя параметры движения инструмента и погрешности, вызванные процессом обработки и несовершенством инструмента. Модель представлена в виде комплекса уравнений, описывающих модель динамики резания и модель формирования новых поверхностей. В качестве входных параметров системы вводится геометрия инструмента и детали, режимы резания и параметры движения, на выходе выдается трехмерная топография обработанной поверхности.
Моделирование формирования погрешностей формы цилиндрической поверхности при хонинговании коротких отверстий рассматривается в [6]. Разработана плоская модель, учитывающая начальную геометрию отверстия, динамику вала, инструмента и брусков хона, поджимаемых к обрабатываемой поверхности.
В [7] предложена и проанализирована математическая модель процесса хонингования, включающая в себя динамику вала и инструмента с учетом взаимодействия поверхности обрабатываемой детали и брусков хона. Показано, что силы рассматриваемого взаимодействия неконсервативные. Хонин-говальный инструмент моделируется как вращающийся стержень с головкой, закреплённой в среднем сечении вала. Исследование зависимостей для сил резания от режимов обработки при хонинговании проводится в [8].
В настоящей работе исследована модель процесса хонингования при обработке отверстий с заданной неровностью поверхности с целью прогнозирования исправления погрешности формы. Разработана модель динамики инструмента, который схематизируется как вращающийся безмассовый стержень, передающий вращение оправке хона с упруго закрепленными на нем режущими и направляющими брусками. Разработана модель сил резания на основе гипотезы Престона. На базе разработанной модели исследовано влияние технологических параметров на интегральные критерии качества обработки.
Расчетная схема инструмента
Хонинговальный инструмент схематизируется в виде консольного безмассового вала длиной К1ш/1' вращающегося с угловой скоростью \л/, которому сообщается возвратно-поступательное движение со скоростью Уа . На валу шарнирно закреплена хонинговальная головка массой М(рис. 6).
Хонинговальная головка состоит из оправки с вязко-упруго присоединёнными к ней брусками хона длиной 1Ьот квадратного сечения х ¡¥/мпе. Оправка обладает физическими моментами инерции относительно трех осей, JQ - момент относительно продольной оси, У, - момент относительно поперечных осей
Бруски совершают движение в радиальном направлении за счет механизма разжатия. Предполагается, что бруски находятся в постоянном контакте с обрабатываемой поверхностью, вращаясь вместе с валом и оправкой с постоянной угловой скоростью. При этом инерцией брусков пренебрегают. Усилия, передающиеся от вала через механизм разжатия на бруски и далее на обрабатываемую поверхность, могут быть схематизированы как линейные вязко-упругие с коэффициентом жесткости к и коэффициентом демпфирования с и зависят от относительного положения бруска и вала в радиальном направлении.
Рассматривается схема расположения брусков хона, показанная на рис. 7. Имеется три бруска:
№4
2008
один является режущим, он расположен так, что направление его перемещения совпадает с положительным направлением вертикальной оси связанной вращающейся системы координат, два других бруска отклонены на угол в от вертикальной оси и являются направляющими.
Хонинговальная головка схематизируется как твердое тело имеющие четыре степени свободы: два перемещения центра тяжести хона щ V относительно оси вращения и два угла поворота J \ у относительно соответствуюхцих осей во вращающейся системе координат. Крутильными и осевыми колебаниями вала пренебрегаем.
Уравнения динамики вала инструмента с оправкой
Вращение оправке сообщается через вал, соединенный с оправкой шарнирно. При этом вал и оправка могут совершать колебания в поперечном направлении.
Уравнения динамики вала инструмента получим в системе координат (г/,у, и>) с началом в точке закрепления вала, вращающейся вместе с валом относительно неподвижной системы координат (х, у, г) (рис. 8).
Рис. 8. Системы координат Рис. 9. Силы, действующие на вал
Разрежем шарнир, соединяющий вал и головку хона, и приложим к валу силы реакции Ки и , а к хонинговальной головке - силы реакции от вала Яи и , силы инерции Р/7 и , результирующие сил взаимодействия брусков инструмента и детали: силы резания РиЪ и Ру1, моменты сил резания и Гг1, а также инерционные и гироскопические моменты М", М^?, М[ , М7 [9] (рис. 9,10, 11).
Я,
I ^г УГ V
р +р" ±ЫЕ х и
Рис. 10. Силы, действующие на хонинговальную головку
Рис. 11. Моменты, действующие на хонинговальную головку
Рассматривая уравнения динамического равновесия хонинговальной головки, после преобразований получим систему уравнений колебаний хона в виде:
и _ $ Ьт1. -Ъп{РиЪ- Мы)
-/Д. + + Тл + (РЛ- Ми)!**. = о
(1)
где 5П статическая податливость вала от действия единичной поперечной силы, приложенной на его конце.
Модель сил резания
Силы резания, или силы взаимодействия инструмента с обрабатываемой поверхностью, будем моделировать как силы сухого трения с коэффициентом трения , тогда результирующую силу взаимодействия можно разложить на нормальную и касательную составляющие (рис. 12). Касательная составляющая направлена в противоположную сторону к относительной скорости движения брусков инструмента и поверхности детали и считается пропорциональной нормальной составляющей с
№4
2008
коэффициентом /и,
Для данной расчетной схемы инструмента (рис. 13) в произвольный момент времени нормальная
составляющая контактного давления пропорциональна поперечным смещениям оправки хона, скоростям этих смещений, а также радиальным смещениям брусков при обработке неровной поверхности и определяется для каждого из брусков как:
ру - р0 + к [и + б- ¿(у-)соз б + с(г) + 0 - с (у - ЭЛ.£) соя 0 - кУю ,
где р0 - постоянная составляющая радиального давления, задаваемая при начальном разжатии брусков хона после введения их в отверстие (радиальная подача), МПа;
к - эквивалентная радиальная жёсткость бруска, Н/Л/3;
с - эквивалентное демпфирование бруска, Н • с/мъ \
и, V- поперечные перемещения центра масс головки инструмента, мм;
- углы поворота головки относительно соответствующих осей,рад\
Е; - координата, отсчитываемая от центра масс вдоль оси отклонённой хонинговальной головки, мм;
- отклонения обрабатываемой поверхности в точках, соответствующих положению центра масс хонинговальной головки вдоль оси обрабатываемого отверстия и положению каждого бруска на его поверхности в настоящий момент времени
Рис. 13. Схема взаимодействия брусков с обрабатываемой поверхностью
Касательные составляющие сил резания на каждом из брусков определяются следующим образом:
= (3)
где .//2'Лг коэффициенты сухого трения направляющих брусков, а коэффициент резания для режущего бруска.
№ 4 2008
Равнодействующие контактных давлений найдём, проинтегрировав проекции давлений
по всей длине хонинговальной головки 1,юпе и умножив их на ширину бруска Whom . Выражения для равнодействующих сил резания и моментов примут вид:
Pul = -IkA sin2 в-и- 1сA sin2 в • и + кА cos /?(/.+ 2 ft cos2 в) v- а, Л cos р (2fr cos 0 ■- ¿) +
cos p(fe + 2 / cos2 0) V- M (- ) sin tí cos /?[/t, 7Л| - ft (7Л2 4- 7,3) cos , Pvl = (l + 2 cos2 v- Ы (l + 2 cos2 0) v- cos (3 sin2 0 • w + (2 cos 0 -1) -2cAf( cos p sin2 0 • й + M [ - (Уи + 7ЛЗ) cos 0] - ОД cos p (YR2 - 7,3) sin 9, (4)
^ - (1 + 2 cos2 0) 0V -(1 + 2 cos2 в) 4 + ЫД sin /?(/. + 2/; cos2 0) • v+ +cAR sin p [fc + 2 ft cos2 e)-v-kARsmP [fc Ym - ft (YR2 + УА.Я) cos в] + -2/ costf),
Tvl = -2kJbsin2 0• i9v-2cJhsin2 2kARff sin/3sin26-u--2cARft sin Psin2 в • it - kARft sin /? sin 6(Yn- YR3),
/ 3
ГДе A — W" / I — W hmc ■
^ — YY hone1 hone J b ^ ^
12
Модель формирования новых поверхностей
Примем, что, согласно закону Престона, скорость съёма материала при резании прямо пропорциональна работе, совершенной силами трения, действующими на элемент поверхности в данной точке [8] (рис. 14), то есть
—^-р—, (5)
сН куа({
где УК = Ун (¿) - радиальная координата обрабатываемой поверхности под режущим бруском;
р - давление бруска на поверхность;
кх, - безразмерный коэффициент, который зависит от свойств материала и определяется эмпирически;
а* - предел прочности материала детали.
Таким образом, новая координата поверхности в каждой точке после прохождения по ней режущего бруска определяется следующим образом:
№4
2008
(6)
где А/ - время контакта бруска с какой-то определённой точкой поверхности;
(-Т/ - момент времени, когда эта точка поверхности обрабатывалась на предыдущем обороте или проходе;
кч> - удельный коэффициент, характеризующий толщину срезаемой поверхности.
При разработке модели образования новых поверхностей будем считать, что толщина слоя, снимаемого за один проход, незначительна по сравнению с величинами поперечных перемещений вала. Таким образом, она не влияет на величину нормальной составляющей сил резания, определяемой выражением (2). Будем также считать, что координата поверхности в каждой точке изменяется после прохода бруска через эту точку (изменение поверхности, определяемое выражением (б), не входит в уравнение динамики инструмента).
л» %
Рис. 14. Схема образования новых поверхностей
Описание полной модели
После подстановки (4) в (1) получим полную систему уравнения динамики инструмента в
матричном виде: [м] X + [£>] X + [ЛГ] X = В,
(и 4
V
где X - К ; [М] =
Мб,
1]
О
о о
о мзи о о
У, о
О У.
(V)
(8),(9)
и
[О]
2д] ¡сАвт2 в - <5, ¡сА соз 2 соб" (9) 2 <5, б иг б'соз ¡3 5],с А (1 + 2 сое2 0)
^ 2сАЯ^\пршгв 0 -«/„<»
^0)
2с]ь 8Ш" в ;
(10)
№4
2008
К
1 + <уЬ4(1+2сог0)
I.
О
km 1
I Г
-^+cARm/3(fiA2f!m:e) ^+Uh(l+2cos29
гЩтрт^-Ь-
V Ц\
о
о
houe 2
О О
I2
—+2М,ж9 4«, *
В =
f -S,,кА ( Yr2 - Yrз ) sin 0 - , cos /?[/.УД1 - ft ( Yln + Yn ) cos 9 kA [ YR] ~(YR2 + YR3 ) cos в] - kAft cos fi ( YR2 -YK3) sin 9
-kARf, sinfisin - YR3 )
(il)
(12)
Обозначим Y =
'Хл
A
Система, приведённая к каноническому виду Коши, имеет вид:
<
или ? = [А]-7 + [М\х-Н-
'О 0 0 0 1 0 О 0^ 0 0 0 0 0 10 0
где
И'
0 0 0 0 0 0 10 0 0 0 0 0 0 0 1
и
[-С]
H =
f о Л о о о
[М]-1 в
(13)
(14), (15)
Полученная неоднородная система уравнений не имеет аналитического решения, так как правые части уравнений включают переменные }?;, зависящие от положения брусков в настоящий момент времени и состояния поверхности под бруском, образованной в предшествующий момент времени,
№4
2008
когда режущий брусок проходил по этой точке поверхности. Поэтому будем искать приближенное численное решение данной системы.
Определение собственных частот системы. Рассмотрим случай обработки идеально гладкой цилиндрической поверхности, т.е. Уг. = 0 и коэффициент срезания км- = 0 , тогда уравнение (13) будет однородным линейным дифференциальным уравнением, описывающим колебания хона при движении по гладкой поверхности без резания.
Частоты колебаний определяются как собственные значения матрицы А. Решение однородного уравнения (13) ищем в виде:
К(Г) = 70ехрЛг,
где X - безразмерное собственное значение матрицы А, т - безразмерное время,
г = р0и р1 = ел{Щы).
В общем случае собственные значения будут комплексными, т.к. система нагружена неконсервативными силами имеющими циклический характер, и, при определенных значениях параметров системы, колебания могут иметь как периодический характер, так и быть затухающими при отрицательных значениях действительных частей собственных значений, или с возрастающей амплитудой при положительных значениях действительных частей X.
Если положить с = 0, р0 = 0, то можно определить частоты собственных колебаний как мнимые части собственные значения матрицы [А] , т.к. действительные части в этом случае равны нулю. Изменение собственных частот в зависимости от жесткости показано на рис. 15 при 0 = 30 , в зависимос-
Н
ти от угла взаимного расположения брусков 0 - на рис. 16 при к = 1.00 • 108 —-, частота вращения
м
/7 = 1000 об/мин.
Eigen Frequwicies
Eigen Fr*Kiu*ncws
ЗпН4
ВШпвх» А, ЫМг3
Рис.15. Влияние жёсткости на собственные частоты
Рис.16. Влияние угла взаимного расположения Брусков 0 на собственные частота
№4
2008
Из частотного анализа системы можно заключить, что увеличение эквивалентной радиальной жесткости инструмента приводит к увеличению всех собственных частот, в то время как при увеличении угла взаимного расположения брусков в до 60° две пары частот попарно сходятся и затем начинают обратно расходится при 0 >60°.
Eigsn Freqii©noies
Eigen Frequencies
1 .......................L
К** ; ¿gggggSle!ШЩ® | i -е-Ри , ■9 Ру ! Pia ' Pto
-............;..•
J_L
0.5
älues
2.6 x 10е
:UCTy......-J_1_I__
0 10 20 30 43 SO ffi 70 BD Real parts^if^nvalues
0.1
0,05
-0Q5 kvs^-M tki! 4-М I Ui J4-M \Ш \Ш U-M U\-UU\Mi-
......... ...... .......1 1
Щ ; ® -+P« i py
..... Щвм -vviWef Wk
-------L. -----
SUttmss k, №n*
Рис. 17
10 203340 60 63 70 80 50 feg
Рис. 18
На Рис. 17 и 18 показаны графики изменения действительных и мнимых частей комплексных собственных значений для случая, когда С= I * 104 Н - с/М3, = 0.05, /с = 0.3 . Как видно из рис. 17, для 6 = 30 действительная часть собственных значений отрицательна, т.е. колебания являются затухающими на всем диапазоне изменения радиальной жесткости хона. При изменении угла расположения направляющих брусков частоты попарно сближаются и в диапазоне от 55 до 65 градусов действительная часть одного из собственных значений становится положительной, т.е. колебания становятся неустойчивыми. Если значения = 0.3,/с = 0.3, диапазон углов при которых колебания неустойчивы расширяется.
Дискретизация поверхности отверстия. Численное моделирование
Так как при численном интегрировании системы уравнений необходимо знать координаты каждой точки поверхности при каждом шаге интегрирования и пересчитывать ее с учетом закона резания, необходимо описать поверхность в виде сетки в окружном и осевом направлении. Поверхность отверстия задаётся в виде матрицы Уг отклонений в радиальном направлении от базовой поверхности (в данном случае базовой является поверхность идеального цилиндра диаметром £)), число элементов которой может достигать десятков и сотен тысяч.
Матрица Уг имеет размер ЛГ, х ЛГ2. Значение (число точек разбиения в окружном направлении) совпадает с числом шагов интегрирования в окружном направлении. Оно выбирается из сообра-
№4
2008
жений точности и устойчивости процедуры интегрирования. Значение И2 определяется по формуле
/V = Л/ . И 4-д/ +7
где я/1(,Л7 - число витков винтовой траектории бруска; Л^- число дополнительных строк.
Если положить = 0, то траектория камня после одного двойного хода будет накладываться сама на себя, и большая часть поверхности не будет участвовать в обработке. Чтобы избежать этого, к матрице добавляется Nn столбцов, наличие которых обеспечивает сдвиг траектории камня с каждым двойным проходом. В результате через некоторое количество двойных ходов все отдельные траектории дадут в сумме регулярную "сетку", покрывающую всю поверхность. За счет протяженного в осевом направлении режущего камня обеспечивается обработка всей поверхности отверстия.
В процессе интегрирования из матрицы выбираются те элементы, положение которых соответствуют текущему положению брусков инструмента. Новая координата поверхности после прохождения инструмента определяется в зависимости от текущего значения силы резания.
При численном интегрировании системы по времени с автоматическим выбором шага мы не имеем точных значений Уг в произвольных точках. Так как текущее положение инструмента соответствует какому-то углу поворота вала с оправкой, то шаг сетки дискретизированной поверхности и соответствующий ему угловой шаг вращающегося инструмента могут не соответствовать шагу по времени, то есть при интегрировании системы по времени можно попасть в точки, в которых нам отклонение поверхности не известно. Можно эти значения подсчитать путем аппроксимации по отклонениям поверхности в соседних точках, но это вносит дополнительные погрешности, что нежелательно.
Для того, чтобы при интегрировании в расчетах учитывались только известные координаты точек поверхности, удобно ввести в качестве аргумента угол поворота.
Если закон вращения детали ср = со/, то тогда ¿/(р = , и производные по углу поворота можно выразить через производные по времени.
Для интегрирования полученной системы дифференциальных уравнений был использован математический пакет МАТЬАВ 6.5. При этом использовался метод Рунге-Кутта 4-5 порядка (стандартная процедура МАТЪАВ ос!е45).
При численном моделировании варьировались режимы обработки и радиальная жесткость хона. Для остальных параметров принимались значения, представленные в таблице 1.
Таблица 1
Постоянные параметры
Название параметра Параметр Значение
Длина отверстия Ь оп1 500 мм
Диаметр отверстия О, 48 мм
Длина хона Ьопе 50 мм
Ширина бруска хона попе 6 мм
Длина вала ЧиаП 500 мм
Диаметр вала Оы 18 мм
Момент инерции хонинговальной головки относительно продольной оси К 6.431е-5 кг*м2
Момент инерции хонинговальной головки относительно поперечных осей 1.149е-4 кг*м2
Масса хонинговальной головки М 0.397 кг
Отношение окружной скорости инструмента к скорости осевой подачи Бреес1_гаио 3
Из кинематики процесса установлено, что режущий брусок хона при таких значениях параметров обрабатывает все точки поверхности хотя бы один раз за четыре двойных хода.
Анализ результатов
Все расчеты проводились для поверхности, имеющей начальную погрешность формы типа «эллиптичность» (рис. 17).
90 0.0004
Рис. 17. Рассматриваемый вид погрешности формы
В качестве критерия точности обрабатываемого отверстия была выбрана величина, называемая
№4
2008
эксцентричностью поверхности, определяемая для каждого сечения отверстия плоскостью, перпендикулярной оси, по формуле
/у _ Ы „ max min ft-—-— • (16)
' RL. + RL 5
max
гДе ^шах и ^ соответственно наибольшее и наименьшее отклонения поверхности от идеального цилиндра. Итоговая характеристика для всей поверхности получается в результате осреднения полученных значений для каждого сечения:
(17>
е =
Величина снимаемого припуска, характеризующая производительность процесса, вычисляется по среднему радиусу, определяемому как среднее значение всех элементов матрицы отклонений поверхности. В результате моделирования исследовано влияние на изменение эксцентричности и среднего радиуса при обработке величины жесткости хона к, постоянной составляющей контактного давления р0 и скорости вращения инструмента п.
Влияние жесткости на эксцентричность и средний радиус
На рис. 18-19 представлены результаты численного моделирования процесса обработки отверстия при разлргчных значениях жесткости подвески к инструмента. Показаны графики изменения эксцентричности и среднего радиуса отверстий в зависимости от количества двойных ходов движения хона при разных значениях жесткости подвески к инструмента.
При увеличении жесткости к наблюдаются увеличение скорости изменения эксцентричности и незначительное уменьшение скорости изменения среднего радиуса. Влияние жесткости на изменение эксцентричности тем сильнее, чем больше постоянная составляющая контактного давления ру.
Частота вращения инструмента не меняет характер изменения эксцентричности в зависимости от жесткости крепления брусков хонинговальной головки.
Параметры; р0 = 1 кПа, п = 1000 б Iмин .
A/s 4
2008
х ю'4 Eccentricity vs variable stiffness к
5 10 15 20 25 30 5 10 15 20 25 30
Number of douhh feeds Number of double feeds
Рис. 18. Изменение эксцентричности в процессе обработки в рис< 19. Изменение среднего радиуса в процессе зависимости от жесткости крепления брусков хона обработки в зависимости от жесткости крепления
nil 7 / Ъ брусков хона
1 : А = 2-10 Я/.и , 2: к = 2.2 -10 И ! м п 7 . .
1 : к =2-Ю7 н/м2 2 : к = 2.2-10 И/лГ
3 : /с =2.4-107 H/л? 7 . .
3: А: = 2.4-10 Н]м л
Mean radius vs variable stiffness к
Влияние постоянной составляющей контактного давления на эксцентричность и средний радиус
На рис. 20-21 представлены результаты численного моделирования процесса обработки отверстия при различных значениях постоянной составляющей контактного давления р0 инструмента, представлены изменение эксцентричности и среднего радиуса отверстия в зависимости от количества двойных ходов хонинговальной головки.
Увеличение постоянной составляющей контактного давления уменьшает скорость изменения эксцентричности (рис. 20), увеличивает скорость изменения среднего радиуса (рис. 21). Параметры: к = 2• 1 о7 н/м3, п = 1000 об!мин.
Влияние скорости вращення инструмента на эксцентричность и средний радиус
Установлено, что скорость вращения инструмента не влияет на скорость изменения среднего радиуса. На рис. 22 показано изменение эксцентричности после 40 двойных циклов обработки в зависимости от частоты вращения инструмента. В результате анализа было установлено, что при наименьшем значении частоты вращения (100 об/мин) окончательная эксцентричность наименьшая, далее при
№4
2008
увеличении скорости она увеличивается, достигает максимума при 500 об/мин, затем начинает медленно уменьшаться. Значения параметров при моделировании: к- 2 • 107 н/л? , р0 = 0.5 кПа .
■¡Eccentricity vs variable constant part of contact pressure pQ
x 10
•оMean radius vs variable constant part of contact pressure p
10 15 20 25
Number of double feeds
12 16 20 24 28 Number of double feeds
Рис. 20. Изменение эксцентричности в процессе обработки в зависимости от постоянной составляющей контактного давления
\:р0 =0.5 кПа, 2:р0=0.75 кПа, 3 : =1 кПа.
Рис. 21. Изменение среднего радиуса в процессе обработки в зависимости от постоянной составляющей контактного давления
1: р0 = 0.5 кПаь 2:р{)= 0.75 3 : = 1 кПа,
х io'd Г/ш/ eccentricity vs variable frequency of rotation n
5.8 5.75 5.7 5.65 5.6
0 500 1000 1500
Frequency of rotation n, rpm
Рис. 22. Изменение эксцентричности после 40 двойных циклов в зависимости от частоты
вращения инструмента
№4
2008
Выводы
По результатам исследования в данной работе можно сделать следующие выводы:
1. Разработанная математическая модель позволяет оценить влияние параметров инструмента и режимов обработки на частотные характеристики колебаний инструмента. Установлено, что при осесимметричном расположении брусков (55 < в < 60 ) процесс хонингования может быть неустойчивым из-за колебаний с возрастающей амплитудой.
2. Анализ численных результатов моделирования показывает, что:
• на исправление погрешности формы наибольшее влияние оказывают жесткость крепления брусков хона к оправке и постоянная составляющая контактного давления;
• при увеличении жесткости А'наблюдается ускорение исправления погрешности формы поверхности, при этом скорость съема материала уменьшается;
• при увеличении постоянной составляющей контактного давления р0 скорость исправления погрешности формы уменьшается, но при этом увеличивается съем материала.
Работа выполнена благодаря финансовой поддержке РФФИ грант № 07-08-00592.
СПИСОК ЛИТЕРАТУРЫ
1. З.И. Кремень, Н.Х. Стратневский. Хонинговаиие и суперфиниширование деталей/Под ред. Л.Н Филимонова.— 3-е изд., перераб. и доп.- Л.:, Машиностроение, Ленингр. отд-ние, 1 98 8.— 1 37 стр. с ил л.— (Библиотечка шлифовщика. Вып. 9).
2. Бабичев А.II. Хонинговаиие.- М.: Машиностроение. 1965, 96 с.
3. Бндерман В.Л. Прикладная теории я механических колебаний. Учебное пособие для втузов. М., «Высшая школа», 1972.- 416 стр. с и л Л.
4. Воронов С. А., Фатальчук А. В. Моделирование динамических процессов в ибр о абразив но й обработки глубоких отверстий.- Механика твердого тела, №6, 2000.- стр. 167-173.
5. Грязное А.Н., Гуськов A.M. Модель процесса образования поверхности при обработке абразивными брусками.— Известия вузов. Машиностроение, № 1 1980.—стр. 130-1 39.
6. K.F. Ehmann, M.S. Hong. A Generalized Model of the Surface Generation Process in Metal Cutting. Annals of the CIRP Vol. 43/1/1994
7. Kim, J.-D. and M.-S. Choi, "A Study on Prediction of Roundness Variation in Honing Using Hone Dynamics," Journal of Testing and Evaluation, Vol. 25, No.6, 1997, pp.556-564,
8. A.M.Gouskov, S.A.Voronov, E.A.Butcher, S.C.Sinha. Non-conservative oscillations of a tool for deep hole honing. Communications in Nonlinear Science and Numerical Simulation. Vol.11, issue 6, September 2006, p. 685-708.
9. Zcttel H.D. "Verbesserung der Arbeitsgenangkeit von Hon werkzengerT* Jnd. Anz., 1975, 97, N 14, pp. 269-270.