CC BY
63
ТЕХНОЛОГИЯ И ОБОРУДОВАНИЕ МЕХ. И ФИЗ-ТЕХН. ОБРАБОТКИ
УДК 621.923
А. А. ИГНАТЬЕВ, В. А. КАРАКОЗОВА, С. А. ИГНАТЬЕВ
ОЦЕНКА ДИНАМИЧЕСКОГО СОСТОЯНИЯ ШЛИФОВАЛЬНОГО СТАНКА ПО АВТОКОРРЕЛЯЦИОННЫМ ФУНКЦИЯМ ВИБРОАКУСТИЧЕСКИХ КОЛЕБАНИЙ ДЛЯ ОПТИМИЗАЦИИ РЕЖИМОВ ОБРАБОТКИ
Предлагается метод оценки динамического состояния шлифовального станка по автокорреляционной функции виброакустических колебаний, регистрируемых при обработке колец подшипников при стационарном режиме резания, причем входное воздействие на систему рассматривается в виде сигнала типа «белый шум». Динамическая система; виб-роакустические колебания; автокорреляционная функция
В рамках системы мониторинга технологического процесса (СМТП) необходимо контролировать динамические характеристики станков, так как они непосредственно влияют на качество обработки, например, на параметры точности и однородность структуры поверхностного слоя дорожек качения колец подшипников [1, 2]. Для оценки динамических характеристик используют измерение виброакустических (ВА) колебаний основных элементов формообразующей подсистемы станка, возбуждаемых в динамической системе (ДС) различными методами [3-7]. Далее осуществляется обработка результатов измерений и оцениваются различные характеристики ДС (см. табл.).
Идентификация методами теории автоматического управления связана с необходимостью подачи на вход ДС типового сигнала. Наиболее часто в известных работах по экспериментальному определению динамических характеристик станков идентификацию передаточных функций элементов и ДС в целом осуществляют с помощью активного эксперимента, возбуждая систему детерминированными воздействиями: ступенчатым, импульсным или гармоническим. Указанные методы обладают рядом недостатков:
• необходимо формировать эталонное входное воздействие калиброванной амплитуды и формы, что не всегда возможно для станков;
• искажение входного и выходного сигналов снижает точность идентификации;
Контактная информация: (8452) 52-78-23
• сложно оценить реальные динамические характеристики станка, поскольку спектры воздействий отличаются от спектра колебаний, возникающих при резании.
При оценке состояния ДС станков применяется и пассивный эксперимент, когда измерения выполняются в процессе обработки и система возбуждается силами резания. При этом выделяется два метода идентификации: первый - по переходным характеристикам при врезании инструмента в заготовку, второй -при стационарном резании. Оба метода разработаны в СГТУ и реализованы для токарных станков [4]. В обоих случаях обязательно учитывается стохастическая составляющая колебаний в ДС, которая по своему спектру близка к сигналу типа «белый шум». Эта составляющая обусловлена тем, что на ДС оказывают воздействие несколько независимых источников ВА колебаний с различными спектрами (электродвигатели, шпиндели, процесс резания и т. п.). Их суммарное воздействие в диапазоне частот до нескольких килогерц определяет стохастическую составляющую силы резания, которую можно рассматривать как «ограниченный белый шум».
Основное преимущество метода, рассматривающего возбуждение ДС силой резания со стохастической компонентой типа «белый шум», является использование реальных воздействий в технологической системе, что существенно повышает адекватность идентифицированной модели ДС станка [8]. Стохастическая компонента ВА колебаний ДС в этом случае рассматривается как стационарный случайный процесс (ССП). Следовательно, возможность применения метода
обусловливается стационарностью режима резания и линейностью ДС, что выполняется как для токарных, так и для шлифовальных станков при финишных операциях обработки.
Для оценки состояния ДС шлифовального станка представляется целесообразным получить выражение для ее передаточной функции [8, 9]. Для этого используется экспериментально-аналитический метод оценки передаточной функции ДС по автокорреляционной функции (АКФ) ВА колебаний, зарегистрированных в процессе шлифования колец подшипников на основных узлах станков, наиболее близких к зоне резания (шпиндели, опора кольца). Из ранее выполненных исследований известно [10], что при воздействии на объект стохастических колебаний АКФ имеют доста-
точную информативность для оценки его динамических свойств.
Измерение ВА колебаний осуществлялось виброизмерительным прибором ВШВ-003М2, датчик которого устанавливался с помощью магнитного основания на узлы станка. Измерения проводились на различных станках моделей SIW-4, SIW-5 при обработке колец различных типоразмеров на операциях предварительного и окончательного шлифования. Запись и обработка результатов осуществлялись компьютером типа Notebook. Вычисление спектров и АКФ производилось с помощью программного продукта Matlab. Выходные данные формируются в виде сдвоенного графика (спектрограмма и график АКФ).
Т аблица 1
Методы идентификации динамического состояния станков
№ п/п Вид входного воздействия Реализация на станках Основные свойства Выходные характеристики Лите- ра- тура
1 Г армоническое Возбуждение колебаний от вибраторов Большие затраты времени, точен АФЧХ, АЧХ, ФЧХ [1-3]
2 Импульсное Возбуждение колебаний ударом специальным молотком Возбуждение колебаний от математического маятника Электромагнитный возбудитель с обратной связью Сложность задания параметров импульса, легок в интерпретации. Выполнено только на лабораторной установке. На станке не опробовано АФЧХ, ИПФ, передаточная функция АЧХ при различных углах удара [1-3, 8]
3 3.1 3.2 Ступенчатое Без резания При резании Переходной процесс при нагружении через батавскую слезку Переходной процесс при врезании инструмента Не соответствует реальным условиям, легок в интерпретации. Чувствителен к возмущениям, легок в интерпретации Переходная и передаточная функции Переходная и передаточная функции, запас устойчивости [1-3]
4 4.1 4.2 Стохастическое Без резания При резании Возбуждение колебаний от специального вибратора Возбуждение колебаний силой резания с регистрацией спектра и автокорреляционной функции (АКФ) Не соответствует реальным условиям, сложен в применении, нечувствителен к возмущениям Прост в применении, нечувствителен к возмущениям, приближен к реальным условиям АФЧХ, АЧХ, передаточная функция Спектр и АКФ при резании, интегральная оценка спектра и АКФ, передаточная функция, степень устойчивости [4-7]
Если принять во внимание, что переходной процесс при врезном шлифовании колец подшипников средних типоразмеров является кратковременным и износом круга при обработке одного кольца можно пренебречь, то ВА колебания узла формообразующей подсистемы (например, опоры кольца), можно рассматривать как ССП. При условии, что по экспериментальным данным о ВА колебаниях получено математическое выражение АКФ, можно вычислить передаточную функцию динамической системы станка при резании [8].
Для этого используется методика, предложенная в работе [11], согласно которой для передаточной функции линейной системы, находящейся под воздействием ССП типа «белый шум» со спектральной плотностью 50, выполняется соотношение
W(p)W(-p) = Куу (р) + Куу (-р). (1)
Предположим, что АКФ идентифицирована в виде
г ^
~(t)=L-'jWp) I=L-l1 к
K (t) = De-
cos w1t + -
a
(2)
причем K,, (0) = D , где D - дисперсия ССП;
Ю1 - круговая частота.
Изображения по Лапласу от АКФ заданного вида легко получаются, если воспользоваться табличными данными из справочника [12]:
K уу ( p ) =
■ +
Da
D (p + a)
(p + a)° + ю° (p + a)° +w2
(3)
После проведения соответствующих алгебраических операций по формулам (1) и (3), связанных с раскрытием скобок и приведением подобных членов, следует, что передаточная функция ДС шлифовального станка при обработке, рассчитанная при входном воздействии типа «белый шум», записывается в виде
к
W (Р) = —---------------------т, (4)
р +ар + ю0
1 1 ^ 0 0 0 где к = —л — , ю0 = к»! +а .
0 V а
Переходной процесс в ДС при врезании круга в заготовку можно рассматривать как сумму двух составляющих: детерминированной и стохастической. Детерминированная составляющая переходной функции к (г) определяется по передаточной функции W(p) при единичном ступенчатом воздействии на ДС шлифовального станка выражением
2 2 \. (5)
p p +ap + ю0
В качестве ступенчатого рассматривается воздействие силы резания при постоянной подаче шлифовального круга. Детерминированная составляющая h (t) соответствует сглаженной переходной функции колебательного звена. Стохастическая составляющая переходного
процесса h(t) , возникающая под действием
сигнала типа «белый шум» X(t), определяется выражением [13]
О 1 ¥
h(t) = — | exp{- a(t - t)}X(t) sin ro(t -t)dt. (6)
w-¥
После завершения переходного процесса этим же выражением описывается стохастический процесс на выходе системы в стационарном режиме. Следовательно, имеет место преобразование стохастического процесса X(t) в
стохастический процесс h(t), причем если исходный процесс является стационарным случайным с нулевым математическим ожиданием, то и преобразованный относится к тому же
типу. Однако явный вид процесса h(t) установить не удается.
Таким образом, реальный переходной процесс является суммой составляющих h (t) и
h(t) , определяемых, соответственно, детерминированной FP (t) и стохастической Fp (t) составляющими силы резания FP(t). После завершения переходного процесса в ДС наблюдается ССП, определяемый составляющей силы резания FP(t) .
Реальные АКФ, зарегистрированные при шлифовании колец подшипников на станках SIW-5, представлены на рис. 1 и 2. Они аппроксимированы с погрешностью 10... 15% выражениями вида (2).
Для станка № 333:
• черновой проход
Куу (t) = 0,82e-136t(cos 120t +1,14 sin 120t);
• чистовой проход
Куу (t) = 0,55е -252t(cos120t + 2,1sin120t).
Для станка 302:
• черновой проход
Куу (t) = 0,02e-744t(cos 120t + 6,2 sin 120t);
sin w1t
• чистовой проход
К (t) = 0,015e-304t(cos 120t + 2,5 sin і20т).
б
Рис. 1. Автокорреляционные функции при обработке на шлифовальном станке № 333: а - черновой проход; б - чистовой проход
Полученные АКФ свидетельствуют о различных динамических характеристиках станков, в том числе при разных технологических режимах. По идентифицированным передаточным функциям можно не только оценивать реальное состояние ДС, но и выбрать целесообразный режим шлифования по таким характеристикам как степень устойчивости и запас устойчивости.
Вычисляя значения коэффициентов передаточной функции по АКФ ВА колебаний при стационарном резании и определяя степень устойчивости ДС станка при различных значениях параметров технологического режима, можно установить режим с наибольшей степенью
устойчивости, а следовательно, и с более высоким качеством обработки. Динамическое состояние ряда станков одной модели целесообразно сравнивать при близких технологических режимах, что позволит сравнить их техническое состояние и выбрать станок с приемлемыми характеристиками. Обе задачи решаются в рамках системы мониторинга технологического процесса [1].
Автокорреляционная функция екзр-11_1 -ОЗ.щач
0.015 , 0.01 0.005 О
-0.005
L. л . ііі 4 Lu ■AAJ
▼ т* *
т
-5000 -4000 -3000 -2000 -1000 0 1000 2000 3000 4000 5000
Время, отн.ед.
Автокорреляционная функция вк8р-11_3-03.теа\
А А ні ІІІ ІІІ in U,
щ ff" fW
I і і
-5000 -4000 -3000 -2000 -1000 0 1000 2000 3000 4000 5000
Время, отн.ед.
б
Рис.2. Автокорреляционные функции при обработке на шлифовальном станке № 300: а - черновой проход; б - чистовой проход
Известно, что построение ДС станка обычно производится при ряде допущений, направленных на упрощение модели, и не может учесть в явном виде всех процессов, протекающих в системе при резании. Следовательно, идентификация позволяет более четко представить реальную структуру модели ДС и выделить ее наиболее важные составляющие, поскольку измерение ВА колебаний выполняется именно в процессе резания. По идентифициро-
а
а
ванной модели можно сформулировать некоторый критерий, оценивающий динамическое состояние станка, а также установить эталонное значение или диапазон изменения этого критерия, соответствующий заданному качеству обработки деталей. Выход критерия за допустимые границы связан либо с существенным нарушением динамического состояния станка, либо с неправильно выбранным технологическим режимом. В обоих случаях качество обработки деталей не удовлетворяет заданному и требуется осуществить корректирующее воздействие, например, выполнить подналадку станка или выбрать целесообразный технологический режим.
СПИСОК ЛИТЕРАТУРЫ
1. Игнатьев, А. А. Мониторинг станков и процессов шлифования в подшипниковом производстве / А. А. Игнатьев, М. В. Виноградов, В. В. Горбунов [и др.]. Саратов: СГТУ, 0004. 104 с.
0. Игнатьев, А. А. Мониторинг технологического процесса на основе автоматизированного контроля динамических характеристик станков /
A. А. Игнатьев, В. А. Добряков, С. А. Игнатьев // СТИН. 0005. № 7. С. 3-7.
3. Кудинов, В. А. Динамика станков /
B. А. Кудинов. М.: Машиностроение, 1967. 360 с.
4. Попов, В. И. Динамика станков / В. И. Попов, В. И. Локтев. Киев: Техника, 1975. 136 с.
5. Бржозовский, Б. М. Точность и надежность автоматизированных прецизионных металлорежущих станков / Б. М. Бржозовский, А. А. Игнатьев, В. А. Добряков, В. В. Мартынов. Саратов: СГТУ, Ч. 1, 1990. 156 с.; Ч. 0, 1994. 160 с.
6. Есипов, Ю. В. Проблемы мониторинга «здоровья» и ресурса конструкций / Ю. В. Есипов, В. М. Мухортов // Тр. VIII Междунар. науч.-техн. конф. по динамике технологических систем. Ростов-на-Дону: ДГТУ, 0007. Т. 1. С. 44-50.
7. Мерзляков, А. А. Способы возбуждения колебаний при исследовании динамики механических систем / А. А. Мерзляков, Н. А. Серков, Р. О. Си-роткин // Тр. VIII Междунар. науч.-техн. конф. по динамике технологических систем. Ростов-на-Дону: ДГТУ, 0007. Т. 0. С. 194-199.
8. Игнатьев, С. А. Обеспечение качества обработки поверхностей качения колец подшипников на основе контроля динамического состояния шлифовальных станков по стохастическим характеристикам виброакустических колебаний: авто-реф. дис. ... канд. техн. наук / С. А. Игнатьев. Саратов, 0001. 16 с.
9. Игнатьев, С. А. Идентификация динамической системы шлифовального станка по автокорреляционным функциям виброакустических колебаний технологической системы для оптимизации
режима обработки / С. А. Игнатьев, В. А. Каракозова, А. А. Игнатьев // Автоматизация и управление в машино- и приборостроении: сб. науч. тр. Саратов: СГТУ, 0008. С. 80-85.
10. Балакирев, В. С. Экспериментальное определение динамических характеристик промышленных объектов управления / В. С. Балакирев, Е. Г. Дудников, А. М. Цирцлин. М.: Энергия, 1967. 030 с.
11. Лебедев, А. Н. Цифровое моделирование и идентификация стационарных случайных процессов / А. Н. Лебедев, Д. Д. Недосекин, Г. А. Стекло-ва. Л.: ЛЭТИ, 1983. 118 с. (Рукопись депонирована в ВИНИТИ, № 6048-83).
10. Корн, Г. Справочник по математике для научных работников и инженеров / Г. Корн, Т. Корн. М.: Наука, 1984. 830 с.
13. Неймарк, Ю. И. Стохастические и хаотические колебания / Ю. И. Неймарк, П. С. Ланда. М.: Наука, 1987. 404 с.
ОБ АВТОРАХ
Игнатьев Александр Анатольевич, доктор технических наук, профессор, заведующий кафедрой «Автоматизация и управление технологическими процессами» Саратовск. гос. техн. унта.
Каракозова Вера Алексеевна,
зав. лаб. каф. «Автоматизация и управление технологическими процессами» Саратовск. гос. техн. ун-та.
Игнатьев Станислав Александрович, канд. техн. наук, докторант каф. «Автоматизация и
управление технологическими процессами» Саратовск. гос. техн. ун-та.
