CC BY
21
УДК 658.516.3:621
Динамическая диагностика процесса обработки методами поверхностного пластического деформирования
В. Л. Заковоротный, В. С. Ми на ков, М. М. Ханукаев
(Донской государственный технический университет)
Рассматривается процесс выглаживания. В зависимости от внедрения выглаживателя в обрабатываемую поверхность изменяется уравнение в вариациях относительно точки равновесия процесса выглаживания. Для обеспечения рациональных параметров настройки процесса предлагается метод, основанный на анализе сигнала виброакустической эмиссии, в котором отображаются параметры градиента физико-механических свойств поверхностного слоя.
Ключевые слова: процесс выглаживания, градиентные свойства, динамика, виброакустическая эмиссия.
Введение. Методы поверхностного пластического деформирования деталей для повышения их эксплуатационных свойств являются достаточно распространёнными в технологии машиностроения [1]. Один из примеров таких технологий — процесс выглаживания алмазными инструментами [2]. В результате алмазного выглаживания изменяются свойства микрорельефа поверхности, распределение остаточных напряжений в поверхностном слое, микротвёрдость и её распределение по толщине приповерхностного слоя и др. Таким образом, в приповерхностном слое формируются функционально градиентные структуры, характеризуемые, например, изменением параметров матрицы жёсткости в зависимости от сближения контактирующих поверхностей. При варьировании технологических режимов, геометрии наконечника алмазного выглаживающего инструмента, физико-механических свойств материала детали, шероховатости обрабатываемой детали и пр. обработанные поверхности могут обладать повышенной коррозионной и износостойкостью, контактной выносливостью и пр. Однако улучшение эксплуатационных характеристик деталей в этом случае зависит от динамики контактного взаимодействия и других факторов.
В статье ставится задача и предлагаются решения для выбора режимов выглаживания алмазными инструментами на основе анализа сигнала временной последовательности колебаний инструмента, сопровождающей процесс выглаживания. Сигнал временной вибрационной последовательности формируется в результате преобразования силовой эмиссии в зоне обработки динамической системой, которая состоит из подсистем со стороны инструмента и заготовки, взаимодействующих через динамическую связь, формируемую процессом выглаживания. В свою очередь, формируемая динамическая связь несёт информацию о свойствах поверхности детали, формируемой в процессе выглаживания.
Математическое моделирование динамики процесса выглаживания. Схема процесса выглаживания приведена на рис. 1. Она соответствует наиболее распространённому способу выглаживания с упругим закреплением инструмента. Эта схема соответствует процессу трения инден-тора с образцом, который достаточно подробно рассмотрен в работах [3, 4]. Там же показано, что упругие и диссипативные свойства подсистемы со стороны инструмента после предварительного нагружения в вариациях относительно точки равновесия могут быть схематизированы в линейном приближении в виде матриц скоростных коэффициентов и упругости. Инструмент для процесса выглаживания состоит из корпуса, в котором через упругие элементы устанавливается алмазный наконечник, осуществляющий процесс выглаживания. При этом силы контактного взаимодействия формируются в результате деформаций алмазного наконечника относительно корпуса через упругие элементы. Одновременно имеют место деформационные свойства корпуса относительно несущей системы станка, которую можно считать абсолютно жёсткой.
Рис. 1. Схема динамической системы процесса выглаживания: а — схема деформационных смещений вершины алмазного наконечника от статической установки Х01 до точки равновесия Х^ ; 6 — функции сближения сил контактного взаимодействия; в — схема контакта наконечника выглаживающего инструмента с заготовкой; г — схематизация деформационных смещений
Если вершина инструмента устанавливается в точке О на рис. 1, а, которой соответствует координата Х01, то в результате конечной жёсткости в установившемся состоянии инструмент
смещается в точку О1. Этой точке соответствует координата точки положения равновесия Х[‘] в отжимающем направлении. В процессе выглаживания в установившемся состоянии наблюдается баланс сил, формируемых в зоне контакта наконечника с образцом, и сил упругости со стороны инструмента. Причём силы зависят от упругих смещений наконечника в направлении Х1 и предварительного его смещения относительно поверхности заготовки Х01. Качественная характеристика зависимости сил от деформационных смещений, названная функцией сближения, приведена на рис. 1, б. Упругие деформационные смещения при заданных силах определяются матрицей жёсткости, и в динамическом режиме они зависят от матрицы скоростных коэффициентов. Поэтому справедлива схематизация, представленная на рис. 1, г. Схема ориентации эллипса жёсткости в контакте приведена на рис. 1, в.
Величина натяга ґр = (Х01 -Х^, интерпретируемого как глубина резания, характеризует
основную координату настройки динамической системы выглаживания, влияющую на силы контактного взаимодействия (рис. 2). На рис. 2 приведены примеры зависимости установившихся сил контактного взаимодействия от величины натяга для заготовки диаметром 80 мм обработанной чистовым точением с = 2,0 мКм.
а) б)
Рис. 2. Зависимость сил от величины натяга при выглаживании (обработка с СОЖ — индустриальное масло 20; скорость выглаживания — 40 м/мин): 1 — Ст45; 2 — 25ХН1М: а — изменение силы в направлении Л!; б — изменение силы в направлении Х2
Приведённые зависимости сил от натяга в пределах высоты микронеровностей заготовки хорошо аппроксимируются параболической кривой, то есть представимы в виде
/1(Х1) = к1(Х0Л-Х1)2-
МХ') = кЛХод-^)2'
где к1Г к2 — экспериментально определяемые коэффициенты, имеющие смысл давления, кг/мм2. Они зависят от геометрии выглаживателя, физико-механических свойств заготовки, её шероховатости и скорости резания.
В неустановившемся состоянии уравнение динамики процесса выглаживания в соответствии с рис. 1 можно представить в виде
(1)
Л2
(2)
где ^ (*1) = {*! (Ход -X,)2,к2 (Ход - ^ )2} — вектор-функция зависимости двух составляющих
сил от текущего значения глубины внедрения наконечника ^ =(Х01 -Х^, причём, Х01 — расстояние между поверхностью заготовки и вершиной наконечника без учёта его упругих деформационных смещений; X = {ХиХ2}г — величина упругих деформационных смещений,
отсчитываемая от точки О; т =
т 0" Хі кГ3" 1 "с1,1 ^2,1
0 , ь = / С =
т Рі,2 2,2 _^1,2 1 и1
— соответственно
матрицы инерционных, диссипативных и упругих коэффициентов подсистемы инструмента; /■(£) = {/^(£),/2(£)}г — силовой шум процесса выглаживания, характеризующийся составляющими сил, не объяснимыми в координатах упругих деформационных смещений инструмента относительно заготовки. Он имеет широкий спектральный состав.
Если в (2) Х01 = const , то в динамической системе устанавливается равновесие X* = {X* ,Х2|г, координаты которого определяются из системы
^1,1^1 + ^2,1^2 = ^1 (*0,1 “ ) '
^1,2^1 + ^2,2^2 = ^2 (*0,1 “ ) 1
(3)
или
Х[ =
1 + 2x^,1 - д/1 + 4хЛд _
где Xi =
^1^2,2 k2C2\
х2 =
2Xi
^=x2Ki-^*)2, ^2^1,1 “ ^1^2,1
(4)
. При анализе (4) необходимо иметь в виду, что мат-
^1,1^2,2 (^2,1 ) ^1,1^2,2 (^2,1 )
рица с в (2) является симметричной и положительно определённой. Кроме этого из двух корней Х[ выбран один, который соответствует области определения нелинейных функций (1). Будем полагать, что в процессе выглаживания установилось стационарное состояние, которому соответствуют точки равновесия X* = {Х*,Х2|г, являющиеся асимптотически устойчивыми. Рассмотрим уравнение в вариациях относительно X* = {х*,Х2}т. При этом учтём, что силы
контактного взаимодействия зависят не только от величины внедрения, но и от колебательных скоростей в вариациях относительно точки равновесия. Нелинейную функцию зависимости сил от деформационных смещений в общем случае можно представить в виде Р {Х^^Х^сИ,дХ2/сИ).
Скорость обработки будем считать неизменной. После замены X (£) = X* +-*'(£) имеем уравнение в вариациях относительно точки равновесия
d2x ,dx / „л * u\
m^-^r + h — + сх = <$[xl,X ] + f(f),
Л2 - Л — -V'-!"- / ■■ Г/' (5)
где ф(л-1,Х*) = ^(х* +Л-,<&■/<#)-Р(X*) — новая нелинейная функция, обладающая свойством ф(0,Х*) = 0. Причём для фиксированной точки равновесия {х = 0,<&■/<# = 0} ф(л-1) = {ф: (л'1),ф2 (л-!)}7 = {0,0}г . Тогда для малых вариаций в окрестности равновесия имеем
(6)
d х , dx ./.ч
m—— + hy — + сух = f(t), 1 dt 1 w
dt
где cy =
ftPi С. . - —— С. , h,. /7, , _д%
С1,1Д ^1,2,1 , hy = 1,1 дхг 1,2 дх2
Зф2 с ——— с 2Д дх1 2,2 _^2,1,1 ^2,2,1 _ • L ^2,1 Зф2 дХ1 ^2,2 дф2 дх2
^1,1,1 ^1,2,1 ^2,1,1 ^2,2,1
суммарные
матрицы упругих и скоростных коэффициентов с учётом линеаризованной реакции со стороны процесса обработки. Матрицы с1 и в отличие от с и /7 не являются симметричными. Более того, их симметричные составляющие могут за счёт реакции со стороны процесса обработки стать отрицательно определёнными. В этом случае равновесие системы будет неустойчивым. Вопрос устойчивости равновесия будет рассмотрен в наших следующих публикациях. Здесь же будем считать, что равновесие устойчиво.
Пусть задан спектральный состав /■(£) = {/![(£),/2 (£)}г в виде спектральной матрицы 5Г г (со) 5Г г (/о)]
11 2 . Тогда из (6) получаем наблюдаемые автоспектры колебаний
5'Г'Л^ ^/2(» (/о)
инструмента, например в направлении Х1
5х1ЛП=-
А (со) + В (со) - С (со)
Кд)2-®2] +4(^і“од)2 [К,2)2-“2] + 4(^2со0,2 )
(7)
где
характеристиче-
А N = V, И [(С2,2Д - яхо2 )2 + (/72,2Дсо)2 ]; В (со) = 5Гг/г (со) [(с2ДД )2 + (/72ДДсо)2 С (со) = 2ке5/г ^ (/о)[(с2 2Д - яхо + У^2 2дш|(с2дд - Л2дд(о/)^, т2 [р2 + 2^соодр + (со0д )2 ] [р2 + 2^2со0 2р + (со0 2 )2 =
= (тр2 + ЬхгР + С1Х1) (я?р2 + М2 21Р + С2 2Д ) - (/712 2р + С12Д ) (/72ддр + С2ДД )2
ский полином системы (6), приведённый к мультипликативной форме.
Матрица 5^г (/о), как показывают экспериментальные исследования, имеет широкий
спектральный состав. Поэтому при заданных и неизменных параметрах подсистемы инструмента в зависимости от градиентов функции ср[х^х/Л) по направлению сближения контактируемых
поверхностей наблюдается смещение собственных частот системы. Другими словами, при изменениях текущих значений градиентов свойств системы в поверхностной области наблюдается смещение собственных частот и соответствующих им всплесков в автоспектре колебательных смещений (со). Поэтому по девиации собственных частот можно оценивать градиенты функций ф(л-1) = (-х^Фг (-х-1)}Г ■ Важно подчеркнуть, что девиация частот наблюдается при изме-
рении колебаний инструмента в любом направлении, так как характеристический полином системы (5) остаётся неизменным.
Алгоритм настройки. При заданной геометрии и физико-механических свойствах заготовки и выглаживателя качество формируемой поверхности зависит от величины припуска, скорости резания и шероховатости поверхности обрабатываемой заготовки. Если двигаться от точки ^ = О на рис. 2, которой соответствует поверхность заготовки, в сторону увеличения заглубления выглаживателя, то изменяются выходные параметры процесса. При больших глубинах формируемая поверхность становится нестационарной, при малых — процесс выглаживания не сглаживает все микронеровности исходной поверхности. Существует оптимальное значение внедрения, при котором выглаженная поверхность имеет минимальное значение микронеровностей. Причём на оптимальное значение внедрения оказывает влияние не столько само внедрение, сколько градиенты функций зависимости сил от внедрения по самому внедрению. Этому условию соответствует и равномерное распределение микротвёрдости по глубине на всех участках поверхности. Стационарным является и распределение остаточных сжимающих напряжений в приповерхностной области [5].
Важно подчеркнуть, что на эффективность процесса выглаживания влияет не столько глубина внедрения, сколько градиенты рассматриваемых силовых функций по глубине. В свою очередь, зависимости, приведённые на рис. 2, обусловлены при прочих неизменных условиях шероховатостью обрабатываемой заготовки. Поэтому при настройке технологической системы и в процессе её функционирования важно обеспечить требуемое значение градиента функции дР1/дХ1, которому для данной заготовки соответствует дР2 /дХ1. Как показано выше, изменение градиентов дР1/дХ1 и дР2/дХ1 отображается в вариациях частоты основных осцилляторов вибрационных последовательностей.
Рис. 3. Смещения корней характеристического полинома АРС-модели по мере увеличения внедрения инструмента в заготовку
Для оценивания девиации частоты основных осцилляторов удобно воспользоваться методами авторегрессионного спектрального анализа (АРС) [6]. Эти методы позволяют в ходе процесса оценивать распределение корней характеристического полинома АРС-моделей в комплексной плоскости и тем самым — отклонения градиентов рассматриваемых функций. В качестве управления или параметра настройки выступает смещение Х01 в (2). Приведём характерный пример смещения корней характеристического полинома АРС моделей вибрационных последовательностей по мере увеличения Х01, которому соответствует для данной упругой системы возрастание внедрения инструмента в заготовку (рис. 3). Значение деформационного смещения инструмента по направлению Х2 не имеет принципиального значения, так как величины деформационных смещений малы по сравнению с диаметром обрабатываемой поверхности. Характеристики смещения корней приведены для случая выглаживания стали 25ХН1М алмазным выглаживателем диаметром 4,0 мм. Заготовка имеет шероховатость Яа = 2,0 мКм. Скорость выглаживания — 120 м/мин, величина подачи — 0,01 мм/об. Диаграммы приведены для подсистемы инструмента, имеющей следующие матрицы жёсткости, диссипации и инерционных коэффициентов:
ст = 100^ ММ 60,0 — мм , = 6,7.10- КГ'С2 ММ 1,1-10 = кГ'с2 ММ
60,0 — мм 800 — ММ _ 1,1-ю-3 кГс2 ММ 1б,о.ю-! кГс2 ММ
ш'1' =
1,4-10-
0,064-10
кГ-с2
мм
кГ -с2
0,064-10 2,0 -10~4
4 кГ • с мм кГ -с2
мм мм
На рис. 3 даны диаграммы двух основных значимых осцилляторов для следующих значений величины внедрения выглаживателя в образец: 1 — 0,5 мкм, 2 — 1,0 мкм, 3 — 1,5 мкм, 4 — 1,8 мкм, 5 — 2,2 мкм. Кроме этого на рисунках выделены центры группирования для каждого значения внедрения и области 95 % разброса их значений.
Выполненные исследования, в том числе проиллюстрированные на рис. 3, дополненные данными о шероховатости обработанной поверхности, позволяют сделать следующие выводы о свойствах частотного состава наблюдаемых в процессе обработки вибрационных последовательностях.
1) По мере внедрения выглаживателя в обрабатываемую поверхность наблюдается смещение частот основных осцилляторов вибрационной последовательности, что естественно, так как возрастает динамическая жёсткость контакта инструмента с заготовкой. При этом частоты основных осцилляторов имеют тенденцию смещаться в высокочастотную область.
2) При увеличении внедрения инструмента в заготовку частоты основных осцилляторов также смещаются в сторону единичной окружности и при этом эллипсы рассеяния расположения корней характеристического полинома АРС моделей возрастают. Это свидетельствует об увеличении коэффициентов затухания в единой динамической системе, состоящей из инструмента и динамической связи, формируемой в области сопряжения инструмента с заготовкой.
3) Замечено, что увеличение радиального биения шпинделя станка и (или) величины припуска приводит к увеличению оси эллипса рассеяния частот (наблюдается уширение спектральной линии колебаний основных осцилляторов). В свою очередь, кинематические возмущения от станка, вариации припуска и свойств заготовки влияют на предельно достижимую шероховатость формируемой поверхности.
4) По мере увеличения внедрения инструмента в заготовку изменяется шероховатость поверхности обработанной детали (рис. 4). Зависимость величины микронеровностей от внедрения носит экстремальный характер. При этом минимальное значение шероховатости наблюдается при величинах внедрения примерно равных 0,6—0,7 от величины микронеровности заготовки. Для конкретного примера эта величина равна (1,2—1,5) мкм. На приведённой иллюстрации жирная линия проведена через математические ожидания. Кроме этого точками показан разброс экспериментальных данных. Характерно, что дисперсия существенно возрастает при приближении глубины внедрения к величине шероховатости исходной поверхности.
5) При увеличении скорости резания некоторые осцилляторы смещаются в высокочастотную область примерно пропорционально росту скорости. Их частоты отличаются от частот собственных колебаний системы, и при вариации внедрения эти частоты, лежащие в низкочастотной области, практически не смещаются в частотной области. Следовательно, существуют пространственно зависимые возмущения, действующие со стороны процесса обработки. На рис. 3 эти частоты показаны едиными кластерами, обозначенными звёздами.
Рис. 4. Изменение шероховатости поверхности обработанной детали по мере увеличения внедрения инструмента в заготовку
Кроме этого, необходимо отметить, что предельно достижимая шероховатость поверхности зависит от точности станка (прежде всего, от радиальных биений шпинделя) и вариаций поверхности по отношению к оси вращения заготовки, обусловленных точностью установки детали и исходной точностью изготовления её геометрии в поперечном сечении.
Заключение. Наблюдение за вибрационной последовательностью колебаний инструмента в процессе выглаживания поверхностей представляет собой надёжный источник информации о состоянии процесса и параметрах качества формируемой поверхности.
Библиографический список
1. Шнейдер, Ю. Г. Технология финишной обработки давлением / Ю. Г. Шнейдер. — Санкт-Петербург: Политехника, 1998. — 414 с.
2. Торбило, В. М. Алмазное выглаживание / В. М. Торбило. — Москва: Машиностроение, 1972. — 104 с.
3. Заковоротный, В. Л. Динамика трибосистем. Самоорганизация, эволюция / В. Л. Заково-ротный. — Ростов-на-Дону: Изд. центр ДГТУ, 2003. — 502 с.
4. Заковоротный, В. Л. Введение в динамику трибосистем / В. Л. Заковоротный, В. П. Блохин, М. И. Алексейчик. — Ростов-на-Дону: ИнфоСервис, 2004. — 680 с.
5. Папшев, Д. Д. Отделочно-упрочняющая обработка пластическим деформированием / Д. Д. Папшев. — Москва: Машиностроение, 1978. — 152 с.
6. Марпл, С. Л. Цифровой спектральный анализ и его приложения / С. Л. Марпл. — Москва: Мир, 1990. — 584 с.
Материал поступил 19.10.2011.
References
1. Shnejder, Yu. G. Texnologiya finishnoj obrabotki davleniem / Yu. G. Shnejder. — Sankt-Peterburg: Politexnika, 1998. — 414 s. — In Russian.
2. Torbilo, V. M. Almaznoe vy'glazhivanie / V. M. Torbilo. — Moskva: Mashinostroenie, 1972. — 104 s. — In Russian.
3. Zakovorotny'j, V. L. Dinamika tribosistem. Samoorganizaciya, e'volyuciya / V. L. Zakovorotny'j. — Rostov-na-Donu: Izd. centr DGTU, 2003. — 502 s. — In Russian.
4. Zakovorotny'j, V. L. Vvedenie v dinamiku tribosistem / V. L. Zakovorotny'j, V. P. Bloxin, М. I. Aleksejchik. — Rostov-na-Donu: InfoServis, 2004. — 680 s. — In Russian.
5. Papshev, D. D. Otdelochno-uprochnyayushhaya obrabotka plasticheskim deformirovaniem / D. D. Papshev. — Moskva: Mashinostroenie, 1978. — 152 s. — In Russian.
6. Marpl, S. L. Cifrovoj spektral'ny'j analiz i ego prilozheniya / S. L. Marpl. — Moskva: Mir, 1990. — 584 s. — In Russian.
DYNAMIC DIAGNOSTICS OF SURFACE YIELDING TREATMENT TECHNIQUES V. L. Zakovorotniy, V. S. Minakov, М. M. Khanukayev
(Don State Technical University)
The smoothing process is considered. Depending on the smoother digging into the work surface, the equation in variations about the smoothing balance point varies. For providing the rational setting of the process, the method based on the vibroacoustic emission signal analysis, in which the gradient parameters of the mechanical-and-physical properties of the surface coating are represented, is offered.
Keywords: smoothing process, gradient properties, dynamics, vibroacoustic emission.
