Технические науки/3.Отраслевое машиностроение                                         

                                                     Кондрашов С.Г.,  Кондрашова И.Г., Иванов И.В.

                                               Херсонский национальный технический университет

Теоретическое исследование автоколебаний при резании

     Рассмотрим случай обработки многолезвийным инструментом с механическим креплением  двух пластин,  расположенным  рядом в одном корпусе и реализующим схему деления припуска,  когда наиболее вероятно возникновение первичных автоколебаний и регенеративного эффекта вследствие отрицательной геометрии пластин 1и 2, ограниченной жёсткости державки  и т.д. Расстояние между пластинами, измеренное по оси инструмента, равно Δ. Первая пластина работает с глубиной t1, вторая t2. Заготовка вращается с угловой скоростью ω0, инструмент двигается вдоль оси со скоростью подачи S.

       Предварительные исследования показали, что волнистость поверхности в этом случае обусловлена в основном колебаниями инструмента в радиальном направлении, поэтому для дальнейших рассуждений принимаем во внимание только колебания по оси y.

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

     Установим взаимовлияние возникающих автоколебаний и их влияние на общие вибрации инструмента.   Пусть колебательная система имеет массу  m, жёсткость K, постоянную затухания С.

    Положения режущей кромки инструмента относительно среднего положения при первом проходе принимаем через  ось У1, при втором через У2, приняв направление вверх за положительное.

     Глубина резания в каждый момент времени при снятии первой стружки:

                            t΄1=t1-y1

      При снятии второй стружки:

                              t΄2=t2-y2+y1

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

    В нашем случае  величина r равна:

                            r = cosφ                                       

где φ – главный угол в плане.

Сила резания равна              Р=rt;

составляющая                        Py=rtcosβ;

введя обозначение                 ry=rcosβ;

будем иметь                           Py=ryt.

Уравнение колеблющейся массы для движения в направлении Y при работе первой пластины имеет вид:

                            my1+cy1+ky1=Pycp-ryy1

Постоянная сила Pycp не влияет на характер колебаний, поэтому:

                            my1+ су1 + (гу + k)y1 = 0

Решение линейного однородного управления второго порядка ищем в виде:

                            y=edt, где d=const

Характеристическое уравнение имеет вид:

                            md2+cd+(k+ry)=0

Корни этого уравнения, в общем случае комплексные, входят попарно сопряженными и равными:

                            d1=-λωy,

                            d2=- λωy ,

где                       λ=с/2m,        .

Общий интеграл уравнения  my1+ су1 + (гу + k)y1 = 0  в случае комплексных корней характеристического уравнения имеет вид:

                            y1=e-λt(A0sinωyt+B0cosωyt)

Подставив начальные условия t=0,  y1=v, получим: 

                            y1=Aye-λtsin(ωyt+φy1)

Поскольку λ>0, то имеет место затухающее движение.

Таким образом, уравнение :

                            y1=Aye-λtsin(ωyt+φy1)

описывает форму волнистой поверхности после первого прохода.

Переменная  глубина резания на втором проходе    

                            2=y1-y2

Сила резания будет изменятся в зависимости от координаты согласно выражению:

                            Py2=ry(y1-y2)

Уравнение движения колебательной системы при втором проходе будет иметь вид:

                            my2+cy2+ky2=(y1-y2)ry ;

                            my2+cy2+(k+ry)y2=ryy1=ry Aye-λtsin(ωyt+φy1).

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

Введём обозначения:

                            f=ryAye-λt/m;

                            P=c/m;

                            q=(k+ry)/m/

Тогда уравнение

                            my2+cy2+(k+ry)y2=ryy1=ry Aye-λtsin(ωyt+φy1)

примет вид:

                            y2+py2+qy2=fsin(ωyt+φy1)

Как известно , общее решение неоднородного уравнения представляется как сумма какого–нибудь частного решения  этого уравнения и общего решения  соответствующего однородного уравнения. Как было показано выше, общее решение однородного уравнения имеет вид:

                            2= Aye-λtsin(ωyt+φy2)

Частное решение неоднородного уравнения ищем в форме:

                            =Mcosβt+Nsinβt

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

                            my1+ су1 + (rу + k)y1 = 0

частота этих колебаний равна ωy. Таким образом, сила Py2 периодически

модулируется частотой ωy, равной частоте колебаний первого резца. Иначе говоря ,    β=ωy .

С учётом последнего выражения уравнения =Mcosβt+Nsinβt примет вид

                            =Mcos ωyt +Nsin ωyt

Подставляя это выражение в исходное дифференциальное уравнение, находим значения М и N:

                                М=-

                           

Прежде чем подставлять найденные значения М и N в равенство:

                            =Mcos ωyt +Nsin ωyt

введём новые постоянные А* и φ, приняв:

                            М=А*sinφ2*, N=A*cosφy*

т.е.

                          

                            tgφ*=M/N

Тогда частное решение неоднородного уравнения можно записать в форме:

Общий интеграл уравнения   y2+py2+qy2=fsin[ωyt+φy1]  равен y2= + ,т.е.

                           

Обозначив

                           

будем иметь

                            .

Уравнения колебательного движения для первого y1=Aye-λtsin(ωyt+φy1) и второго  резцов выведены при условии, что они работают независимо друг от друга. Однако в рассматриваемом случае оба резца жестко закреплены в одном корпусе и, следовательно, результирующее колебательное движение У корпуса инструмента с резцами подчиняется принципу суперпозиции и будет равно сумме колебаний, возбуждаемых первой и второй пластинами, т.е.

                            y=y1+y2=Aye-λt (sinωyt+sin(ωyt+ψy)+Dysin(ωyt+φ*y))

где ψy=φy2-φy1 - сдвиг фаз между колебаниями первой и второй пластин.

       С другой стороны, фазовый угол ψу можно представить в виде:

                            ψy0t0

где t0-время работы первой пластины до момента врезания второй,      ω0-угловая скорость вращения .В свою очередь    t0=L/V,

где  L- путь резани,  V- линейная скорость.

       Эти величины выражаем следующим образом

                            L=nx2πD/2

                            V=ω0D/2

где nx – число оборотов, которое совершит заготовка до момента врезания второго резца, D-диаметр обработки.

Число оборотов:

                            nx=Δ/S

где Δ- расстояние между пластинами,   S- подача в мм/об.

Подставив эти величины в уравнение ψy0t0, получим:

c учётом выражения Ψy= общие колебания инструмента будут описываться уравнением:

                           

     Рассмотрим условия, при которых выражение обращается в нуль, т.е. отсутствуют вибрации инструмента. Очевидно, при наличии первых колебаний величина Ау будет всегда отлично от нуля. Поэтому для выяснения условий виброгашения приравняем к нулю выражение в скобках и режим его относительно ψу. Предварительные расчёты показывают, что величина Dy на два порядка меньше амплитуд первых двух синусоид, поэтому в дальнейших расчётах с большой степенью точности величиной Dysinyty*) можно пренебречь.

Таким образом,

                            sin ωyt+sin(ωyt+2πΔ/s)=0

                            sin ωyt=- sin(ωyt+2πΔ/s)

откуда

                           

где n-любое целое число.

Окончательно будем иметь: Δ/S= (2n+1)/2

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

Литература:

1.     Василенко Н.В. Теория колебаний:  учебное пособие для студентов технических вузов. К.: Высш. шк., 1992 – 430 с. 

2.     Родин П.Р. Геометрия режущей части спирального сверла. К, 1971. 136с.

3.     Кондрашов С.Г. Диссертация на соиск. уч. ст. канд. техн. наук. Киев, КПИ – 1990, - 250 с.