Математическое моделирование. Оптимальное управление Вестник Нижегородского университета им. Н.И. Лобачевского, 2013, № 6 (1), с. 185-188
УДК 534.1; 621.9
О ГАШЕНИИ КРУТИЛЬНЫХ КОЛЕБАНИЙ В ОДНОЙ МЕХАНИЧЕСКОЙ СИСТЕМЕ
© 2013 г. А.В. Грезина, В.Н. Комаров
Нижегородский госуниверситет им. Н.И. Лобачевского
Поступила в редакцию 24.10.2013
Рассматривается задача гашения крутильных колебаний расточной борштанги при обработке глубоких отверстий на токарном станке. Для исследования возникающих при работе вибраций построена математическая модель, описывающая вынужденные и самовозбуждающиеся крутильные колебания борштанги. Проведен анализ результатов численного моделирования.
Ключевые слова: борштанга, крутильные колебания, динамический гаситель.
Гашение паразитных вибраций в машинах и механизмах является одной из основных проблем современного машиностроения, разрешение которой едва ли возможно без предварительного теоретического исследования. В ряду таких проблем стоит и задача о гашении вынужденных колебаний и автоколебаний расточных борштанг, применяемых для обработки глубоких отверстий. В силу специфики режущего инструмента и процесса обработки при растачивании глубоких отверстий при наиболее производительных режимах резания проявляются, в основном, крутильные колебания, которые могут быть как вынужденными, так и самовозбу-ждающимися [1].
Для построения математических моделей, описывающих эти колебания, используется эквивалентная механическая модель токарного станка [1], предназначенного для растачивания глубоких отверстий. В дальнейшем предполагается, что глубокое отверстие растачивается длинными полыми расточными борштангами, для которых отношение глубины (Г) к диаметру отверстия (й) отвечает условию Г / й >> 5 . Наличие внутреннего канала в системе обеспечивает либо подвод смазочно-охлаждающей жидкости (СОЖ), либо отвод стружки в зависимости от принятой схемы обработки. При этом доминирующим колебательным элементом является стебель борштанги с расточной головкой, а обрабатываемая деталь и неподвижные части станка предполагаются недеформируе-мыми телами. Деталь крепится к патрону, который вращается вместе со шпинделем со скоростью п об/мин, а борштанга - к суппорту, который перемещается с подачей £ мм/об. Инструмент, применяемый для растачивания глубоких отверстий, называется расточной головкой. Она
состоит из стального корпуса, в котором крепятся резцы и направляющие.
Из специализированной литературы [1,2] известно, что:
1. При растачивании глубоких отверстий возникают сложные динамические явления, связанные с процессом резания металла и трением направляющих элементов о поверхность вращающейся детали. Возникающие при этом силы и моменты сил зависят от угловой скорости вращения детали, технологических параметров резания, геометрии инструмента, физических свойств смазочно-охлаждающей жидкости и др.
2. Растачивание глубоких отверстий всегда сопровождается вынужденными крутильными и поперечными колебаниями борштанги. Поэтому представляют интерес исследования влияния динамических гасителей на эти колебания.
3. В процессе обработки глубоких отверстий из-за пониженной крутильной жесткости стебля возникают крутильные автоколебания. Основной причиной самовозбуждения колебаний является отрицательное трение. Эффект регенерации колебаний при резании по следу в данном случае не проявляется из-за отсутствия возможности поперечных колебаний расточной головки.
Наличие малого количества работ по гашению колебаний расточных борштанг, на наш взгляд, не в последнюю очередь связано с необходимостью рассмотрения математических моделей, описывающих колебания борштанги в виде систем дифференциальных уравнений в частных производных с краевыми условиями. В связи с этим представляет интерес построение простых и информативных дискретных моделей, для исследования которых можно применять известные и хорошо зарекомендовавшие себя методы.
п л? чКДО
Рис. 1
Для построения математической модели и исследования вынужденных и самовозбуждаю-щихся крутильных колебаний рассмотрим эквивалентную механическую модель [3], представленную на рис. 1.
В общем случае для математического моделирования колебательных процессов записываются дифференциальные уравнения в частных производных с краевыми условиями, описывающие вынужденные крутильные колебания и возбуждение автоколебаний стебля борштанги.
В настоящей работе для решения поставленной задачи строится модальное уравнение, которое несмотря на его простоту отражает специфику распределённого объекта и содержит информацию обо всех его параметрах: погонной массе стебля, моменте инерции его поперечного сечения, длине борштанги и др. Такая модель позволяет получить условия, при которых могут возникнуть вынужденные или само-возбуждающиеся крутильные колебания.
При написании математических моделей предполагается следующее [1]:
• обработка отверстия ведётся по схеме, когда деталь вращается с постоянной скоростью п об/мин, а борштанга перемещается вдоль оси х с постоянной подачей £ мм/об;
• подача £ суппорта с борштангой всегда остаётся много меньше длины борштанги Г = 1350 мм;
• длина стебля борштанги много больше его диаметра;
• плоские сечения стержня (стебля) при колебаниях остаются плоскими;
• расточная головка может совершать крутильные и продольные перемещения как твёрдое тело и не может перемещаться в поперечном направлении;
• угловые перемещения сечения относительно оси х обозначим через ф( х, t).
Математическая модель с учетом сделанных выше предположений представляет собой дифференциальное уравнение в частных производных следующего вида [3,4]:
рТ д2ф(хt) 2 ОТ Э3ф(хt)
рТР о1Отр
дt2
- от,
д 2ф(х, t)
дх2
дtдx
= М (Г, t),
где М(Г,t) = М0е'ш - периодический возмущающий момент; р, ОТ- плотность и крутильная жесткость поперечного сечения стебля борштанги; 81 - коэффициент внутреннего трения материала.
Для задания краевых условий учтём, что борштанга с одного конца жестко закреплена в суппорте в точке х = 0, а на втором конце при х = Г закреплена жестко расточная головка, на которую в процессе растачивания отверстия действуют два момента: момент сил резания МР , зависящий от угловой скорости вращения детали и технологических режимов резания, и момент сил трения Мс между направляющими расточной головки и обработанной поверхностью [1, 5]. Линеаризованные выражения моментов Мс и МР , записанные в окрестности стационарного значения угловой скорости вращения детали, имеют вид:
АМс = Ьсф (Г, t), АМР = ЬР ф (Г, t),
дМг,
ь. ==дМ'
где Ьс = , Р
дю дю
При этом краевые условия запишутся следующим образом:
ф(0, t) = 0,
т
д 2ф(Г, t)
Г дх2
= ОТр 81
д2ф(Г, t)
дх2
- ОТ,
дф(Г, t) дх
(2)
- Ьс ф (Г, t) - ЬР ф (Г, t),
(1)
где ТГ - момент инерции расточной головки относительно оси х.
Для построения модального уравнения и выявления условий возбуждения определенного типа колебаний применим проекционный метод, представляющий собой проецирование операторного уравнения на координатные функции. В качестве координатной функции воспользуемся первой собственной формой крутильных колебаний борштанги. Согласно проекционному методу приближённое решение линейного дифференциального уравнения в частных производных Г1(ф(х, t)) = 0 представляется в виде ф(х, t) = Ф(x)q(t), где Ф(х)- первая собственная форма крутильных колебаний. Условие ортогональности имеет вид
^ Г1 (Ф(x)q(t ))Ф( х)йх = 0. (3)
0
В результате несложных преобразований получим модальное уравнение
) + с^^) + ) = М(Г, t )Ф(Г), (4)
где
L
т1 = рУр ^ Ф2(х)сХ + У Г Ф2(£),
0
L
с1 = GJP 811 Ф'2(х)сХ + (Ис + ИР )Ф2(£),
0
L
к1 = GJP ^ Ф'2 (x)dx.
0
Параметры простого на вид модального уравнения (4) содержат достаточно ёмкую информацию об исследуемом распределенном объекте.
Согласно критерию Рауса-Гурвица состояние равновесия системы q = 0 будет устойчиво, если все коэффициенты характеристического уравнения, соответствующего уравнению (4), положительны. Если в системе результирующее трение отрицательно (кс + кР < 0), то в процессе растачивания возникают крутильные автоколебания.
Таким образом, при растачивании глубоких отверстий могут возникнуть вынужденные крутильные колебания стебля борштанги с режущей головкой (если состояние равновесия устойчиво в малом по Ляпунову), или крутильные автоколебания в противном случае.
Методы динамического гашения колебаний состоят в присоединении к исследуемому объекту динамических гасителей с целью изменения его вибрационного состояния. При действии вибрационных нагрузок более широкого частотного диапазона применяют динамические гасители с трением. Динамические гасители применимы для всех видов колебаний: продольных, поперечных и крутильных. При этом вид колебаний динамического гасителя соответствует виду подавляемых колебаний [1, 2].
Для исследования гашения вынужденных колебаний борштанги запишем математическую модель системы с динамическим гасителем в следующем виде:
т1 41 + С1 41 + С20?1 - 9 2 ) + (5)
+ *191 + к 2 (41 - 42) = ЩЬ, t )Ф(£),
т2 92 + с2(92- 91) + к2(92 -91) = 0
где т1, с1, к1 — параметры модального уравнения, описывающего колебания основной системы; т2, с2, к2 — момент инерции, коэффициенты диссипации и жесткости динамического гасителя.
С помощью методики, описанной в работе С.П. Тимошенко [2], проведем подбор оптимальных параметров гасителя для гашения резонансных амплитуд крутильных колебаний основной системы. Для этого после несложных
преобразований получим выражение, определяющее отношение амплитуды вынужденных колебаний основной системы к ее статическому отклонению, в следующем виде:
-9- = (4ц282у2 + (у2 -82)2)х
Дст
х (у2(^(82 -у2) + ц28(1 -у2)- (6)
-ц 2Ру 28)2 + [у 2(^ц 28 + Р82) -
- (у2 - 1)(у2 -82)]2)-1.
В выражении (6) величины приведены к безразмерному виду следующим образом: ДСТ = = М0/ к1 — статическое отклонение основной системы, вызванное воздействием статического момента М0; р0 = ^к1 / т1 — собственная частота колебаний основной системы; р2 = ^к2 / т2 — собственная частота гасителя; Р = т2/ т1 — отношение момента инерции гасителя колебаний и приведенного момента инерции основной системы; 8 = р2 / рО — отношение собственных частот гасителя и основной системы; у = га / рО — отношение частоты возмущающего момента к собственной частоте основной системы; ц1 = = с1/2т1 рО, ц2 = с2/2т2рО — параметры, определяющие диссипацию основной системы и динамического гасителя.
На рис. 21 представлены кривые АЧХ основной системы при условии, что в ней есть диссипация (кривые, соответствующие значениям
0 2 ’ ’ 0 4 ' ' 0.6 " 0.8" " Г ' 1.2 ' ' 1.4 У
(1;= 0 Ц> = 0.129
Ц:= 0.086
д>= 0.0433
Рис. 2
= 0.0433
Ц!= 0.086
Ц!= 0.129
Рис. 3
Рис. 4
ц1 = 0.0433; 0.086; 0.129), а диссипация гасителя отсутствует (ц2 = 0).
На рис. 3 представлены АЧХ основной системы при подборе параметров гасителя для случая, когда диссипация в основной системе от-сутствует(ц1 = 0). Такая ситуация реализуется тогда, когда «отрицательное трение», вызванное процессом резания, компенсирует силы трения между направляющими головки и обработанной поверхностью. Расчет АЧХ выполнен при оптимальном значении параметра диссипации гасителя (ц 2 = 0.1912).
Проведённый расчёт показал, что в том случае, когда рассеивание энергии в основной системе равно нулю (случай, когда Ъс + ЪР = 0), относительная амплитуда вынужденных колебаний основной системы не превышает 8 единиц.
На рис. 4 представлены АЧХ, когда в основной системе есть диссипация (коэффициент диссипации положительный). В данном случае невозможно осуществить расчёт оптимальных параметров по методике Тимошенко. Поэтому воспользуемся в расчётах оптимальными параметрами динамического гасителя (ц 2 = 0.1912).
Анализ АЧХ показал, что при использовании динамического гасителя с оптимально подобранными параметрами вынужденные колебания ос-
новной системы могут быть уменьшены в 3—4 раза.
В случае, когда в основной системе есть диссипация, а «отрицательное трение» велико, основная система становится неустойчивой и в процессе растачивания возникают крутильные автоколебания.
Примечание
1. На рис. 2—4 по оси абсцисс отложено отношение частоты возмущающего момента к собственной частоте основной системы, по оси ординат — отношение амплитуды вынужденных колебаний к статическому отклонению основной системы.
Список литературы
1. Уткин Н.Ф. Обработка глубоких отверстий. Л.: Машиностроение, 1988. 269 с.
2. Тимошенко С.П. , Яхт Д.Х., Уивер У. Колебания в инженерном деле. М.: Машиностроение, 1985. 471 с.
3. Грезина А.В., Комаров В.Н. О гашении вибраций при растачивании глубоких отверстий. Труды IX Всероссийской научной конференции им. Ю.И. Неймарка «Нелинейные колебания механических систем», Н. Новгород, ННГУ. 2012. С. 299—305.
4. Пановко Я.Г. Основы прикладной теории упругих колебаний. М.: Машиностроение, 1980. 335 с.
5. Зорев Н.Н. Расчёт проекций силы резания. М.: Машиностроение, 1958. 53 с.
ON TORSIONAL VIBRATION ABSORPTION IN ONE MECHANICAL SYSTEM
A. V. Grezina, V.N. Komarov
The problem of boring bar torsional vibration absorption is considered for deep hole drilling on a lathe. A mathematical model is built to describe forced and self-excited torsional vibrations of a boring bar. An analysis of numerical simulation results is presented.
Keywords: boring bar, torsional vibrations, dynamic absorber.