Влияние динамики процесса сверления на формирование погрешностей глубоких отверстий Текст научной статьи по специальности «Механика и машиностроение»

УДК 62:531.391.3 В.С. БЫКАДОР

ВЛИЯНИЕ ДИНАМИКИ ПРОЦЕССА СВЕРЛЕНИЯ НА ФОРМИРОВАНИЕ ПОГРЕШНОСТЕЙ ГЛУБОКИХ ОТВЕРСТИЙ

Рассматривается механизм формирования нелинейной динамической характеристики процесса сверления в радиальном направлении на различных этапах сверления. Анализируется влияние на нелинейную динамическую характеристику некоторых факторов. Особое внимание уделено динамике колебательных движений сверла в радиальном направлении с позиции формирования погрешностей отверстия.

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

Введение. Во многих областях машиностроения используются детали, имеющие глубокие отверстия [1]. Для получения глубоких отверстий, имеющих отношение глубины I к диаметру D равное или более 5, используются сверла, обладающие малой жесткостью в радиальном направлении особенно при обработке отверстий малого диаметра. Это зачастую приводит к появлению брака, связанного с геометрической неточностью отверстия, прежде всего с отклонением действительной оси отверстия от его идеальной оси и отклонениями диаметра отверстия в поперечных сечениях.

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

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

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

Рис.1. Эквивалентная схема механической системы

Необходимо проанализировать деформационные смещения вершины инструмента во вращающейся системе координат, для этого представим силы, действующие на сверло в координатах упругих деформаций инструмента, тогда в общем виде [2] уравнения динамики сверла в радиальных направлениях будет иметь вид:

+ Н — + СУ = PY (у,5Р,УР),

где И = т1 0 - матрица приведенных масс, Н-с2/мм; Н = И1 0

0 т2 0 ^2

коэффициентов диссипации, Н-с/мм; С

0 с-

(1)

- матрица приведенных

матрица приведенных коэффициентов

жесткости, Н/мм; Y = ||_у1, у 2|| - вектор отклонений сверла в радиальных направлениях,

мм; Рч (ч, ^р, VP) = ||РУ12 (ч, ^р, VP), РУ22 (Y, ^р, Ур |т - вектор-функций динамических характеристик (сил резания) по радиальным направлениям, Н ; Sр - подача пиноли, мм/рад; Ур - скорость резания, мм/с. Отметим, что подача Sр и скорость Ур задают стационарную траекторию движения инструмента. Матрицы М, Н и С определяют параметры подсистемы инструмента в ортогональных радиальных направлениях.

Следующие обстоятельства раскрывают связь деформационных смещений вершины инструмента с точностью отверстия.

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

- если координаты состояния не только смещаются относительно оси вращения, но при этом наблюдается изгиб инструмента, то возможно формирование нового направления движения инструмента, что приведет к уводу оси отверстия.

Базовые уравнения динамики. Моделирование упругих радиальных деформаций режущего инструмента не представляет сложностей, если рассматриваются колебания в частотном диапазоне, ограниченном первой формой колебаний. Параметры уравнения (1) М, Н и С могут быть идентифицированы на основе анализа колебательных реакций системы на импульсные 5 -образные воздействия и статических измерений жесткости сверла. Особое значение имеет представление радиальных составляющих сил РУ12 (у1, у2) и РУ22 (у1, у2) в координатах упругих деформаций, определяемых относительно стационарной траектории. Для раскрытия динамических характеристик примем базовую гипотезу образования сил, утверждающую, что силы на i -ом режущем лезвии сверла пропорциональны площади fi срезаемого слоя материала заготовки соответствующим режущим лезвием сверла. Данная гипотеза не противоречит положениям теории резания металлов.

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

но точки закрепления 01 на угол а01 точки 02 (рис.2).

и одновременного поворота на угол а02 относительно

0

с

Рис.2. Схема деформирования сверла в радиальном направлении

Рассматриваемая схема деформирования сверла приводит к нелинейному изменению площадей £ (у1) и f2 (у) срезаемых слоев каждым режущим лезвием сверла (рис.3). В начале

упругого деформирования сверла влияние угла а01 на смещение вершины инструмента будет больше (за счет вылета сверла), чем угла а02, что приведет к уменьшению площади £ (у) и увеличению площади £ (у1), т.е. £ (у1) < £ (у1) (рис.3,6). По мере деформирования сверла, влияние угла а02 на смещение вершины инструмента будет увеличиваться, а угла а01 -уменьшаться, что приведет к увеличению площади £ (у1) и уменьшению площади £ (у1), в результате £ (у1 ) = £ (у1) (рис.3,в). При дальнейшем деформировании в направлении у1, площадь £ (у1) будет продолжать увеличиваться, а площадь £ (у1) - уменьшаться, в результате £ (у1 )> А2 (у1) (рис.3,/). Приведенная схема деформирования рассматривается для инструмента, обладающего симметричной геометрией режущей части, например, неравенство углов в плане ф1 и ф2 может привести к дополнительным изменениям площадей срезаемого слоя

А (у1) и £ (у1) и соответственно к изменению сил резания.

Рис.3. Влияние отклонения сверла у1 на площади срезаемых слоев £ (у1) и f2 (у1):

? - а 01 _ а 02 _ 0 ; б- а 01 ^ а 02 ; В - а 01 ^ 0, а 02 ^ 0 • £1 ~ £2 ; Г~ а 01 ^ а 02

Тогда динамическая характеристика Р^ (у1), в силу базовой гипотезы, будет определяться значениями площадей £ (у1 ) и £2 (У1 ) :

Ру1х (У1) = Ру(1)(У)-Ру(2)(У), РУ )(У1)=р у А, (У1) і = 1,2,

(2)

(3)

где Ру)(у1) - радиальная составляющая силы резания на / -ом режущем лезвии в направлении у1, Н; ру - коэффициент резания (удельная сила резания) по радиальному направлению

У1

, МПа; А, (у1) - площадь срезаемого слоя на і -ом режущем лезвии,

мм

Одновременно изменение площадей срезаемых слоев £1 (у1) и £2 (у1) приводит к обра-

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

возникают как составляющая силы PY1Y (у1, у2), так и силы (у1, у2), действие этих сил бу-

дет направлено по ортогональным направлениям у1 и у2 соответственно.

Для дальнейшего анализа важны следующие обстоятельства, которые приведены ниже:

- упругие деформации сверла будем рассматривать как стержень круглого сечения, известно, что в этом случае, независимо от выбранной системы координат, матрицы М, Н и С в системе (1) являются не только диагональными, но и имеют равные значения коэффициентов т1, ^ и с, деформация же сверла осуществляется по направлению действия силы. В этом случае можно выбрать произвольную систему координат, при этом матрицы М, Н и С являются независимыми от выбранной системы координат;

- примем допущение о том, что ориентации сил резания в зависимости от деформации сверла являются постоянными.

Эти обстоятельства позволяют при оценке деформации сверла выбрать новую систему

координат у0 - у2, повернутую относительно системы координат у1 - у2 на угол у, который определяется ориентацией силы в системе координат у1 - у2, тогда деформационные смещения в направлении у0 будут вызывать изменение силы, действующей в этом же направлении

(рис.4).

Рис.4. Преобразование исходной системы координат у1 — у2 к новой системе координат У10 — у 2

Можно проанализировать движение системы в направлении у0 (у):

Уравнение (4) на фазовой плоскости запишем относительно фазовых координат:

■ = г

dt т

— (ртх (у ) — — су )

(5)

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

На рис.5,а показаны нелинейные зависимости изменения площадей £ (у) и А2 (у), а на рис.5,6 соответствующая им нелинейная зависимость динамической характеристики резания Рта (у) в радиальном направлении, которая вычислена по выражениям (2) и (3) при у1 = у, данные для расчета имели следующие значения: диаметр сверла В = 12 мм, вылет сверла Ь = 150 мм, углы в плане ф1 = ф2 = 60°.

Рис. 5. Нелинейные зависимости: а-площадей £1 (у) и £2(у); б- динамической характеристики Рта (у) без учета функции сближения FS ; в) - динамической характеристики Рта (у) с учетом функции сближения FS

В процессе сверления по мере заглубления сверло соприкасается своей боковой поверхностью (ленточками) с цилиндрической поверхностью образуемого отверстия. В результате

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

Fs = g ^(у^в ,ws ), (б)

FS ^ FS max : ^(y) ^ ^mln ’ /зacв ^ /зacв max , (7)

где 5(у) - расстояние между поверхностями сверла и отверстия в радиальном направлении, мм; /^ - глубина засверливания, мм; Ws - ширина ленточки сверла, мм, Ws = const.

Влияние функции сближения Fs на динамическую характеристику Р^ (у) проявляются

в образовании участков {- B; - A} и {A; в} (рис.5,в).

.у, мм

з)

у, мм

6)

Рис.6. Влияние на динамическую характеристику (у) ■

а - глубины засверливания /засв ; б- двойного угла в плане 2ф

На рис.6 показано изменение динамической характеристики РуХ (у) в зависимости от глубины засверливания ^ (рис.6,а) и двойного угла в плане 2ф сверла (рис.6,б).

При увеличении глубины засверливания /3aci5 точки - A и A (рис.5,в) приближаются к значению PYS (у) = 0 , что является непосредственным проявлением влияния функции сближения Fs , значение которой увеличивается с увеличением глубины засверливания /3aci5 . Приведены фазовые портреты численного решения системы (5) для значений динамической характеристики PYS (у) (рис.б,а). Как можно наблюдать при малых значениях глубины засверливания

(рис.7,а), при незначительном влиянии функции сближения Fs на сверло, на фазовой плоскости образуются три точки равновесия Q1, Q2 и Q3 со своими областями притяжения. Эти области разделены седлообразными сепаратрисами.

-0.04 -0.02 0.00 0.02 0.04

у, ММ

з)

у; ММ

Рис. 7. Фазовые портреты решений системы (5) для динамической характеристики PYS (у), соответствующей различным глубинам засверливания / 3ac]5 : а- /3aci5 = 5 мм; б- PYS (у) мм; в - /3£юв = 30 мм; г - /3£юв = 60

мм

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

сверливания /засв , точки равновесия Q2 и Q3 вырождаются, а область притяжения точки Q1

занимает всю фазовую плоскость.

Проанализируем преобразование фазовых портретов в зависимости от величины двойного угла в плане 2ф сверла (рис.8).

у, мм

а)

у, ММ

-0.04 -0.02 0.00 0.02 0.04

V; ММ

у, мм

I)

Рис. 8. Фазовые портреты решений системы (5) для динамической характеристики Ру£ (у), соответствующей различным значениям двойного угла в плане 2ф сверла: а- 2ф = 120о ; б- 2ф = 130о ; в- 2ф = 140о ; г- 2ф = 150о

На фазовых портретах (рис.8) также отчетливо проявляются три стационарные точки б\, Q2 и б3 с областями притяжения соответственно. Но в отличие от ранее рассмотренного случая увеличение угла 2ф способствует расширению областей притяжения точек равновесия б2, б3 и сужению области притяжения стационарной точки , таким образом, вероятность

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

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

например, с уменьшением угла 2ф затрудняется отвод стружки из зоны резания вследствие ее пакетирования. Таким образом, фазовые портреты показывают, что на величину области притяжения точки равновесия Q1 принципиальное влияние оказывает величина двойного угла в плане 2ф (рис.8).

Все рассмотренные выше случаи относились к режущему инструменту, обладающему идеальными и симметричными режущими свойствами. Однако если режущие свойства лезвий сверла отличаются друг от друга, то имеет место смещение точки равновесия Q1.

Изменение режущих свойств лезвий сверла может быть вызвано следующими факторами:

- погрешность заточки лезвий сверла, проявляющаяся в наличие эксцентриситета вершины сверла относительно его идеальной оси;

- неравномерность развития износа;

- анизотропия физических свойств обрабатываемого материала;

- неравномерность поступления СОЖ в зону резания и прочими факторами.

Влияние на смещение точки равновесия Q1 и изменение фазового портрета системы

асимметрии геометрии режущей части сверла проявляется в наличии эксцентриситета вершины сверла относительно его идеальной оси.

На рис.9,6 показано изменение динамической характеристики Рта (у) в результате эксцентриситета вершины сверла (е ф 0) (рис.9,з).

V, мм

а) 6)

Рис. 9. Влияние неравенства углов в плане фі и Ф2: а- на образование эксцентриситета, е ; б- динамическую характеристику, (у)

Эксцентриситет е вызывает не только изменение областей притяжения точек равновесия (область притяжения точки Q3 больше, чем точки Q2), но и смещение точки равновесия Q1 относительно начала координат фазового портрета, т.е. относительно идеальной оси вращения сверла, что приводит к возникновению геометрических погрешностей отверстия (рис.10).

.у

I 100

тз

і

тз

50

0

-50

-100

Рис.10.

Выводы. Динамическая характеристика Р^ (у) зависит от координат системы и является нелинейной. Анализ влияния на Р^ (у) глубины засверливания /засв , двойного угла в плане 2ф и

наличия асимметрии режущей части сверла, а также анализ решений уравнения колебательного движения системы в радиальном направлении под действием различных модификаций динамической характеристики Р^ (у) показал ряд особенностей процесса сверления.

1) В начале сверления при врезании сверла в заготовку система чувствительна к вариациям радиальных отклонений, что способствует стабилизации системы в точках равновесия, находящихся не в начале системы координат (0; 0) фазового портрета, но по мере увеличения

заглубления вероятность стабилизации системы в точке с координатами (0; 0) также увеличивается. С технологической точки зрения важно обеспечить точность позиционирования сверла в точке с координатами (0; 0) в момент его врезания и на начальном этапе сверления. По мере заглубления сверла будет усиливаться влияние функции сближения Fs на динамическую характеристику РуХ (у), что приведет к расширению области притяжения точки Q1 за счет вырождения стационарных точек Q2 и Qз, тем самым будет сформирован естественный механизм стабилизации сверла вблизи его идеальной оси.

2) Увеличение двойного угла в плане 2ф сверла способствует расширению областей

притяжения точек Q2 и Qз, это может спровоцировать отклонение вершины сверла от его идеальной оси, однако чрезмерное уменьшение угла 2ф может вызвать другие негативные факторы процесса сверления, поэтому выбор величины угла 2ф необходимо выполнять исходя из некоторого компромисса.

0.04

0.02

0.00

0.02

0.04

у, мм

Фазовый портрет решений системы (5) для динамической характеристики Ру2 (у) при наличии эксцентриситета (в > 0) вершины сверла

3) Наличие эксцентриситета е вершины сверла относительно его идеальной оси, вызванного неравенством углов в плане ф1 и ф2, приводит к тому, что стабилизация изображающей точки в начале координат фазовой плоскости (0,0) невозможна, т.е. при асимметрии режущих свойств лезвий сверла всегда существует отклонение вершины сверла от его идеальной оси, что приводит к образованию геометрических погрешностей отверстия, которые могут проявляться в виде разбивки диаметра или увода оси отверстия в зависимости от величины упругих изгибных деформаций сверла.

4) Область притяжения точки Q1 является, в общем случае, ограниченной. Характерно, что область притяжения точки Q1 зависит не только от величины упругих смещений, но и от их скорости. Отметим также, что область притяжения точки Q1 не является монотонной при варьировании упругих деформаций сверла.

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

Библиографический список

1. Троицкий Н.Д. Глубокое сверление / Н.Д. Троицкий. - Л.: Машиностроение, 1971. -

176 с.

2. Заковоротный В.Л. Динамика процесса резания. Синергетический подход / В.Л. Зако-воротный, М.Б. Флек. - Ростов н/Д: Терра. 2006. - 876 с.

3. Заковоротный В.Л. Динамика трибосистем. Самоорганизация, эволюция / В.Л. Заковоротный. - Ростов н/Д: Издательский центр ДГТУ, 2003. - 502 с.

Материал поступил в редакцию 24.09.10.

References

1. Troickii N.D. Glubokoe sverlenie / N.D. Troickii. - L.: Mashinostroenie, 1971. - 176 s. - In

Russian.

2. Zakovorotnyi V.L. Dinamika processa rezaniya. Sinergeticheskii podhod / V.L. Zakovorotnyi, M.B. Flek. - Rostov n/D: Terra. 2006. - 876 s. - In Russian.

3. Zakovorotnyi V.L. Dinamika tribosistem. Samoorganizaciya, evolyuciya / V.L. Zakovorotnyi.

- Rostov n/D: Izdatel'skii centr DGTU, 2003. - 502 s. - In Russian.

V.S. BYKADOR EFFECT OF DRILLING OPERATION DYNAMICS ON DEEP HOLE ACCURACY

Radial non-linear dynamic performance of drilling operation on various stages is considered. Effect of some factors on non-linear dynamic performance is analysed. Special attention is paid to dynamics of radial oscillatory

motions of the drill taking into account hole accuracy.

Keywords: non-linear dynamic performance, dynamics of drilling operation, cutting process.

БЫКАДОР Виталий Сергеевич (р. 1980), старший преподаватель кафедры «Автоматизация производственных процессов» Донского государственного технического университета. Окончил Ростовскую-на-Дону государственную академию сельскохозяйственного машиностроения (2002), Донской государственный технический университет (2007).

Область научных интересов: динамика и диагностика технологических процессов, теория управления.

Автор 4 публикаций.

vit-bull@mail.ru

Vitaly S. BYKADOR (1980), Senior Lecturer of the Production Automation Department, Don State Technical University. He graduated from Rostov-on-Don State Agricultural Engineering Academy (2002), Don State Technical University (2007).

Research interests: technological processes dynamics and test operation, control theory.

Author of 4 publications.