CC BY
38
Математическое моделирование. Оптимальное управление Вестник Нижегородского университета им. Н .И. Лоба невского, 2013, № 6 (1), с. 180-184
УДК 534.1; 621.9
© 2013 г.
ОБ УСТОЙЧИВОСТИ ОДНОЙ РАСПРЕДЕЛЕННОЙ СИСТЕМЫ
А.В. Грезина
Нижегородский госуниверситет им. Н.И. Лобачевского
aleksandra-grezina@yandex.ru
Поступила в редакцию 24.10.2013
Построена математическая модель, описывающая самовозбуждение крутильных и поперечных колебаний расточной борштанги, замкнутая на динамику процесса резания. Исследование устойчивости состояния равновесия системы проведено методом D-разбиения. Получены области устойчивости и предельные глубины резания при конкретных значениях геометрических параметров борштанги и технологических режимах резания.
Ключевые слова: математическая модель, устойчивость, D-разбиение, самовозбуждение колебаний, борштанга.
Известно, что растачивание глубоких отверстий на токарных станках при наиболее производительных режимах резания сопровождается вибрациями, приводящими к ряду вредных явлений [1, 2]. Возникновению вибраций наиболее подвержена борштанга как менее жесткое звено технологической системы станка. Она в процессе обработки может совершать продольные, крутильные и поперечные как вынужденные колебания, так и автоколебания. Снижение интенсивности колебаний является одним из основных путей повышения производительности и стойкости инструмента. В целях упрощения при исследованиях обычно не учитывают влияние одних колебаний на другие, хотя в действительности влияние имеет место [1]. Это связано с тем, что математические модели, описывающие взаимосвязь колебаний, представляют собой системы дифференциальных уравнений в частных производных со сложными краевыми условиями, замкнутые на динамику процесса резания. Исследование таких систем представляет трудности теоретического характера. В настоящей работе для исследования устойчивости состояния равновесия системы, описывающей взаимосвязь крутильных и поперечных колебаний, перейдем от распределенной модели к модели с сосредоточенными параметрами, что значительно упрощает процедуру исследования.
1. Построение математической модели
Эквивалентная механическая модель борштанги в распределенной идеализации изображена на рис. 1.
В связи с большой глубиной обрабатываемых отверстий (I / d >> 5, где I, d - глубина и диаметр отверстия) борштанга представляется в
виде упругого стержня с непрерывным распределением массы с учетом ее жесткого крепления к суппорту на одном конце и специфики процесса металлообработки на другом. Предполагается, что материал борштанги обладает свойством однородности и изотропности и справедлива гипотеза плоских сечений. В соответствии с приведенными на рисунке обозначениями ф(х, t), у(х, t) и z(x, t) описывают крутильные и поперечные колебания борштанги относительно осей х, у и г соответственно. В
процессе обработки глубоких отверстий при колебаниях борштанги формируются динамические силы резания Fрез и моменты резания
Мрез и трения МТ . Одни из них связаны с силами трения направляющих, а другие - с динамическими силами резания. Имеется ряд причин [1], вызывающих поперечные автоколебания: во-первых, скорость резания обычно находится в зоне, отвечающей падающему участку зависимости силы резания от скорости; во-вторых, колебания толщины срезаемого слоя, приводящие к изменению сил резания, которые изменяются с запаздыванием по отношению к изменению толщины среза и возбуждают и усиливают автоколебания. Возникновение крутильных автоколебаний может происходить под влиянием тех же факторов, что возникновение поперечных автоколебаний.
Так как исследование устойчивости состояния равновесия системы, описывающей взаимосвязь крутильных и поперечных колебаний, затруднено в распределенной идеализации, то для упрощения задачи перейдем от системы с распределенными параметрами к системе с сосредоточенными параметрами, обладающей таким же энергетическим запасом, что и данная сис-
- п об/
мин
Рис. 1
тема. Рассмотрим случай обработки деталей из хрупких материалов, когда сила и момент резания не зависят от скорости, поэтому основной причиной самовозбуждения крутильных и поперечных колебаний при растачивании таких деталей является запаздывание, а роль трения не велика [1].
Так как борштанга представляет собой упругий стержень, один конец которого жестко закреплен в суппорте, а на другом имеется режущее лезвие, на которое в процессе обработки при колебаниях действуют динамические силы и моменты резания, то для упрощения модели можно заменить борштангу материальной точкой с координатой х = I с приведенными массой и моментом инерции. Поперечные и крутильную жесткости борштанги учтем введением
пружин с соответствующими коэффициентами жесткости (рис. 2).
На рис. 2 обозначены: q1(t), q2(t)- линейные перемещения приведенной массы вдоль осей координат у и z, а q3 ^) - угловое перемещение относительно х. Параметры приведенной системы (масса тпр, момент инерции JIsр, коэффициенты жесткости су, с2 и сф) находятся
из равенств кинетических и потенциальных энергий исходной и приведенной систем. Для определения параметров приведенной системы учтем, что самовозбуждение поперечных и крутильных колебаний борштанги при растачивании глубоких отверстий происходит на частотах, близких к первым собственным частотам системы [2, 3]. Коэффициенты рассеивания
\тЯГ(р)
Рис. 3
Ь.ълы.
Рис. 4
энергии приведенной системы hy, к и кф можно получить, используя следующие соотношения:
hy = к = (тпр ®0 А)/ л
и
кф= ( Jпр®0фА)/ Л
где ю0, ш0ф - первые собственные частоты
приведенной системы, А - логарифмический декремент затухания колебаний, полученный по экспериментальным кривым затухающих колебаний в НИИ ПМК ННГУ [2].
Согласно приведенной эквивалентной механической модели (рис. 2) математическая модель, описывающая взаимосвязь крутильных и поперечных колебаний борштанги, запишется в виде [4]
тпр ^ 0) + тпр гЯз 0) + ^1(t) + ^(t) = ~^у,
тпр ^) + hq 2 () + cq1(t) = -АF:,
(Jпр + тпр г2)Ъ(t) + тпр Щ(t) +
+ кфq з(t) + сф qз(t) = -АМ, где обозначено: к = ку = к, с = су = с;
АМх = АFzг + АFyr, г - внешний радиус борштанги, г - статическое смещение приведённой массы борштанги под действием статической силы резания Fz при растачивании [2], АFy и
АFz - приращения проекций динамической силы резания на оси координат.
Приращения проекций динамической силы резания АFy и АFz определяются динамическими силами, вызываемыми относительными колебаниями резца и детали в зоне резания в данный момент, а также силами, связанными со срезанием следа, и могут быть записаны с учетом основной причины самовозбуждения колебаний - запаздывания, следующим образом [2]: АFy = К21^) + Щз^) -
- Ц(#1(t -т) + Щз^ - т))),
^ = К32 (^) + ^3(ї) - (2)
- ц(ql (/ - Т) + Г^(ї -
где т - запаздывание, равное времени одного оборота обрабатываемой детали, К22, К32 - характеристики резания, ц - коэффициент перекрытия (0 < ц < 1).
Проведенные экспериментальные исследования показывают [2], что статическое смещение приведённой массы под действием статической силы приблизительно составляет 2 мм, тогда можно записать, что |К32 |г >> |К22|г . Поэтому для упрощения дальнейших исследований последний член в выражении приращения момента резания ДМх можно не учитывать.
2. Исследование устойчивости
Исследование устойчивости состояния равновесия системы проведем методом Э-разбиения [4]. Для вывода характеристического уравнения запишем частные решения системы в виде ді (ї) = Сіері (і = 1,2,3) и подставим в систему уравнений (1). В результате несложных преобразований получим характеристическое уравнение, представляющее собой квазиполином следующего вида:
Р4( р) + (К32 ггр( р) + К22 Р2( р))(1 - це -РТ) = 0, (3) где
Р4(Р) = Р1(Р)Р2(Р) - т1рГ2Р“,
Р1( Р) = тпр Р2 + ¥ + ^
Р1(Р) = Р1(Р) - тпрР",
Р2 (Р) = (Лр + тпр Г2) Р 2 + Р + Сф ,
Р2(Р) = Р2(Р) - тпрГ2Р".
Из экспериментальных исследований известно [2], что при растачивании острозаточен-ными резцами характеристики резания остают-
V, м/мин
Рис. 5
ы
Рис. 6
ся пропорциональными глубине резания Ь , поэтому их можно записать как К22 = ЬК22 и
К32 = ЬК32. Тогда в характеристическом уравнении можно выделить параметр Ж = -1/ Ь, и уравнение кривой D-разбиения будет иметь вид:
Ж(р) = (р2(Р)К22 + Р(Р)К32гт)(1 ~рт) , (4)
Р4( Р) ’
где р = /ш, юе (-го,+го) .
Для отыскания предельной глубины резания воспользуемся методом D-разбиения по одному комплексному параметру. Условие грубости при данной структуре характеристического уравнения не нарушается и его можно не проверять, так как множество предельных точек кривой D-разбиения при ю ^ го стягивается в точку Ж = 0 . Отсюда следует также, что на плоскости Ж отсутствуют области D(гo), может быть, за исключением одной точки Ж = 0 .
На рис. 3 представлено D-разбиение на плоскости комплексного параметра Ж для ю е (0,+го). Вторая часть кривой D-разбиения для ю е (-го,0) получается зеркальным отображением построенной кривой относительно реальной оси.
Кривая D-разбиения построена при следующих значениях параметров системы: длина борштанги I = 3м, внешний и внутренний радиусы
борштанги г = 0.04м, гвнутр = 0.035м; скорость
резания V = 90м/мин; подача £ = 0.1мм/об .
Учитывая правило штриховки, найдем область устойчивости D(0). Пересечение кривой D-разбиения с отрицательной действительной полуосью определяет предельную глубину резания Ь* .
На рис. 4 приведен график зависимости предельной глубины резания от запаздывания т при подаче £ = 0.1мм/об, скорости резания V = 70м/мин , длине борштанги I = 3 м и радиусах борштанги г = 0.04м, гвнутр = 0.035м . На рис. 5 представлены графики зависимости предельной глубины резания от скорости резания при длине борштанги I = 3м, радиусах борштанги г = 0.04м, гвнутр = 0.035м . На рис. 6 изображены графики зависимости предельной глубины резания от скорости резания при длине борштанги I = 3.5м и радиусах борштанги
г = 0.035м, гвнутр = 0.03м.
Анализ результатов исследований показывает, что
• с ростом глубины резания система становится неустойчивой;
• с увеличением длины борштанги её жесткость уменьшается, что уменьшает область устойчивости;
• с увеличением длины борштанги и запаздывания предельная глубина резания уменьшается.
Результаты численных исследований хорошо согласуются с результатами экспериментальных исследований, приведенными в работе [1] и полученными в лаборатории динамики систем НИИ ПМК ННГУ [2] .
Список литературы
1. Уткин Н.Ф. Обработка глубоких отверстий. Л.: Машиностроение, 1988. 265 с.
2. Городецкий Ю.И. Создание математических моделей сложных автоколебательных систем в станкостроении // В сб.: «Автоматизация проектирования» / Под ред. акад. В.А. Трапезникова. М.: Машиностроение, 1986. Вып. 1. С. 203-221.
3. Грезина А.В. Об устойчивости борштанги при растачивании глубоких отверстий // Тез. докл. XII Ме-ждунар. конф. «Устойчивость и колебания нелинейных систем управления». Москва, ИПУ. 2012. С. 112-113.
4. Неймарк Ю.И. Динамические системы и управляемые процессы. М.: Наука, 1978. 336 с.
ON STABILITY OF ONE DISTRIBUTED SYSTEM A. V. Grezina
A mathematical model describing self-excitation of boring bar torsional and transverse vibrations is presented. The model is predominantly focused on the cutting process dynamics. The stability of the system equilibrium state is studied using the D-partition technique. Stability domains and maximum depths of cut for specific values of the boring bar geometric parameters and various cutting modes are obtained.
Keywords: mathematical model, stability, D-partition, self-excitation of vibrations, boring bar.
