DOI: 10.12737/article_5968b45239a2e4.40971667
Бешевли О.Б., аспирант, Бушуев Д.А., канд. техн. наук, Дуюн Т.А., д-р техн. наук, доц., Рубанов В.Г., д-р техн. наук, проф. Белгородский государственный технологический университет им. В.Г. Шухова
ИМИТАЦИОННОЕ МОДЕЛИРОВАНИЕ ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК МОДУЛЯ ДЛЯ ОБРАБОТКИ КРУПНОГАБАРИТНЫХ ПОДШИПНИКОВ СКОЛЬЖЕНИЯ
tanduun@mail.ru
Предложена технология ремонтной обработки антифрикционной опорной поверхности крупногабаритных подшипников скольжения с использованием разработанного специального приставного модуля, которая позволяет существенно сократить время обработки, повысить качество обрабатываемой поверхности и увеличить срок межремонтного цикла работы крупногабаритных вращающихся агрегатов. Представлен кинематический анализ работы кривошипно-кулисного механизма предложенной конструкции модуля, позволяющий конструктивно предусмотреть необходимые технологические параметры механической обработки. Построена математическая модель и виртуальный прототип приставного модуля. Проведены исследования динамики виртуального прототипа приставного станка, позволяющие оценить сил моментные и скоростные режимы фрезерования, а также скорости перемещения на этапе рабочего хода, когда постоянство скорости обеспечивает стабильность процесса резания, соответствующего величине требуемой подаче режущего инструмента, и повышенной скорости перемещения на этапе холостого хода, способствующее сокращению общего времени обработки изделия.
Ключевые слова: крупногабаритные подшипники скольжения, ремонтная обработка, фрезерование, математическая модель, кинематический анализ, виртуальный прототип, технологические параметры, режимы резания._
Введение. Крупногабаритные подшипники скольжения нашли широкое применение в качестве опорных узлов мельниц самоизмельчения, сушильных барабанов и другого крупногабаритного вращающегося оборудования, используемого в строительной и горнорудной промышленности.
Для повышения эффективности ремонтной обработки антифрикционного слоя крупногабаритных подшипников скольжения предлагается использовать фрезерование, а для его реализации специальный переносной станок (рис. 1), конструкция которого защищена патентами [1, 2]. Замена традиционно используемого ручного шабрения на фрезерование обеспечивает существенное повышение производительности обработки при возможности получения стабильных параметров точности выдерживаемого размера и качества обрабатываемой поверхности [3-5].
Основная часть. Станочный модуль содержит корпус, выполненный в виде двух рам 1 и 2, опирающихся в нижней части опорными элементами, выполненными в виде специальных проточек 3, на соответствующие им проточки сектора подшипника скольжения 4. В верхней части рамы 1 и 2 соединены приводным валом 5, который опирается на их подшипниковые узлы 6.
Приводной вал 5 жестко соединен с конечным звеном кривошипно-коленного механизма 7 привода маятниковой круговой подачи. В привод входят редуктор 8 и электродвигатель 9, которые жестко закреплены болтовым соединением на полке 10, установленной на раме 1. На приводном валу 5 жестко закреплена рама 11, в нижней части которой расположен рабочий узел 12, имеющий направляющие 13 с установленным на них линейным модулем 14. Последний через направляющие 13 связан с шаговым двигателем привода продольного перемещения 15. На линейном модуле 14 закреплена фрезерная бабка 16 с закрепленной фрезой 17. На фрезерной бабке 16 расположен привод главного движения 18, который передает вращение фрезе 17. На концах специальных проточек рам 1 и 2 установлены лазерные датчики 19. На полке 10 закреплен управляющий блок 20 на базе микроконтроллеров с устройством ввода-вывода, входы которого соединены с лазерными датчиками 19, а выходы связаны с шаговым двигателем привода продольного перемещения 15 и линейным модулем 14 через панель управления 21.
В процессе проведения ремонтной обработки важной технологической задачей является обеспечение качества опорной поверхности
скольжения, к которой предъявляют весьма жесткие требования. Качество поверхности скольжения обуславливает коэффициент и условия трения-скольжения, удержание смази, износостойкость и, как следствие, прямолинейность
оси эксплуатируемого агрегата посредством минимального изменения величины зазора в подшипнике в процессе работы, а, следовательно, надежность работы подшипника скольжения и эксплуатируемого агрегата в целом.
Рис. 1. Конструктивная Добиться требуемого качества обрабатываемой поверхности скольжения возможно, осуществляя процесс механической обработки баббита в определенном заданном интервале технологических параметров (табл. 1). Скорость и глубина резания, величина подачи режущего инструмента при фрезеровании баббита формируют микрорельеф поверхности, определяют температурно-силовые характеристики процесса резания, обуславливают износ режущего инструмента.
схема и общий вид станка
В отличие от традиционного фрезерования, где продольная подача - это прямолинейное поступательное перемещение инструмента, особенностью предлагаемого специального оборудования является перемещение инструмента по радиальной поверхности сегмента крупногабаритного подшипника скольжения. Кинематически движение продольной подачи реализовано с использованием четырехзвенного кривошипно-ку-лисного механизма, расчетная схема которого представлена на рис. 2.
Рис. 2. Четырехзвенный шарнирный кривошипно-кулисный механизм
Вращательное движение кривошипу передается от электродвигателя через редуктор. Кулиса и коромысло преобразуют вращательное движение кривошипа в качательное движение стойки коромысла ВЕ, на которой в точке Е расположен режущий инструмент. Таким образом, формируется движение продольной подачи режущего инструмента. За один полный оборот кривошипа режущий инструмент совершает один рабочий и один холостой ход, вернувшись в исходное поло-
Технологические парам«
жение. Движение поперечной подачи осуществляется в конце каждого полного цикла работы механизма - рабочего и холостого хода. Для отвода инструмента и выполнения поперечной подачи использованы шаговые электродвигатели, которые включаются при пересечении фрезы с лучом лазерных датчиков. Главное движение -вращение фрезы вокруг собственной оси реализовано с использованием в качестве привода асинхронного электродвигателя с частотно-импульсным управлением.
Таблица 1
ы фрезерования баббита
Диа- Кол-во Подача Глубина Шероховатость Стойкость Мощность Скорость
метр фрезы, зубьев фрезы на зуб фрезы, резания, мм поверхности Яа, мкм фрезы, мин резания, кВт резания, м/мин
мм мм/зуб
черновая обработка
40 6 0,12 3,7 (4) 10 94 0,95 126
45 6 0,12 3,6 (4) 10 83 0,99(1,06) 141
50 6 0,12 3,5 10(7,6) 73 1,1 157
получистовая обработка
40 6 0,12 3,1 (3,3) 6,3 94 0,76 126
45 6 0,12 3,0 (3,3) 6,3 83 0,84(0,91) 141(150)
50 6 0,12 2,8 (2,3) 6,3(3,4) 73 0,9(0,76) 157
чистовая обработка
40 6 0,12 2,2 (2,5) 3,2 94 0,56(0,62) 126
45 6 0,12 2,0 (2,4) 3,2 83 0,59(0,69) 141(166)
50 6 0,12 1,6 (2,2) 3,2 73 0,56(0,73) 157(183)
Примечание: значения даны для баббита Б16, в скобках - для баббита Б83
Для достижения требуемых условий процесса фрезерования кинематические параметры используемого специального станочного модуля должны обеспечивать заданные технологические режимы. С этой целью разработан алгоритм кинематического расчета модуля, позволяющий определить выходные рабочие кинематические характеристики исполнительных элементов в зависимости от соотношений конструктивных элементов.
Входными параметрами для расчета являются: а - угловая скорость вращения кривошипа ОуА ; 11,12,1з, 14 - длины кривошипа 0\А , связующего звена АВ , плеч коромысла О2В и О2С , соответственно;« - угол при вершине коромысла; I и И - смещение оси коромысла 02 относительно опоры О1 по горизонтали и вертикали, соответственно.
Рис. 3. Расчетная схема четырехзвенного шарнирного кривошипно-кулисного механизма
Начало системы координат расположено в точке О2, направления осей показаны на рис.3. Тогда координаты опоры О1 будут следующими:
а шарнира A -
где
Х01 =1 У01 — -h,
\Xa = Х01 + 11 sinp IУА = У01 -l1co(
( — at.
(1)
(2)
(3)
Координаты шарнира В могут быть определены из условий, что он находится на расстоянии 12 от шарнира А и на расстоянии /3 от оси коромысла О2 :
Jx| + yB = I2
l(xB - xA )2 + (Ув - УА )2 = l2 .
(4)
Полученная система уравнений имеет следующее решение:
CiXA - дСХА^ХАГУАнС^УА^2)
2 , 2
xА + УА
Ув
С1 - XaXB УА
(5)
где
С=1 (3 -122+xA+уА
(6)
Определяем координаты хР и уР мгновенного центра скоростей звена (МЦС) АВ - точки Р (рис.3). Для этого запишем уравнение прямой, проходящей через точки А и О1:
У = уА+1 (x - xa ) - h xA -l
(7)
и уравнение прямой, проходящей через точки O2 и B :
y = yBx. xB
(8)
Решая совместно уравнения (7) и (8), находим точку их пересечения, являющуюся МЦС звена АВ
Хр —
lyA + hXA
УА + h + Ув (l - ха ) Хв
(9)
ур = Ув хр .
Хв
Определяем расстояния от шарниров А и B до МЦС звена AB - точки P :
вр = ^(хв - xP)2 +(Ув - Ур)2
(10)
Располагая найденными расстояниями шарниров А и В от мгновенного центра скоростей, получим:
угловую скорость звена АВ:
cAB =
vA AP'
(11)
где vA —c 11 ; скорость шарнира B
VB — CABBP ,
угловую скорость коромысла
V
B
CO2 =-j~, l3
(12)
(13)
и, наконец, скорость относительного перемещения режущего инструмента по обрабатываемой поверхности
vC = c02 l4 .
(14)
Угловая скорость a 02 меняет знак при тех значениях угла поворота кривошипа p , при которых величина Ув принимает экстремальные (минимальное, максимальное) значения.
Результаты кинематического расчета исполнительного механизма использованы для построения имитационной модели и виртуального прототипа в программном комплексе MSC Adams с целью кинематического и динамического анализа [6-10]. Модель механизма в среде MSC Adams представлена на рис. 4. Она состоит из четырех твердых тел в соответствии с расчетной схемой (рис.3).
Звенья 0iA, AB, B02, 02С заданы стандартными телами системы MSC Adams типа link, массы и моменты инерции которых определены по детализированной модели, созданной в системе автоматизированного проектирования KOMnAC-3D. Звено B02 связано с Ground при помощи RevoluteJoint, а со звеном 02С шарниром типа FixedJoint. Остальные звенья связаны между
x
в
собой посредством Re volute Joint. К звену OiA посредством RotationalJointMotion приложена постоянная скорость вращения п=0,092 мин1 в соответствии с табл. 1. В качестве внешних сил
В результате выполнения кинематического анализа механизма определена траектория движения и скорость перемещения режущего инструмента относительно обрабатываемой поверхности. Рабочие длины кривошипа, кулисы и коромысла подобраны таким образом, чтобы максимальный угол качания фрезы был больше угла сегмента на определенную величину, необходимую для врезания и перебега инструмента, а также для осуществления продольной подачи (рис. 6). На графике рис. 5 отчетливо видны этапы рабочего хода и холостого перемещения
были заданы силы сопротивления в поворотных шарнирах, а в точке С приложена сила резания К направленная по нормали к звену О2С.
инструмента в зависимости от угла поворота кривошипа: на участке угла поворота 0°-180° выполняется быстрое перемещение при холостом ходе, на участке угла поворота 180°-360° - движение с рабочей подачей. Постоянное значение скорости перемещения на этапе рабочего хода обеспечивает стабильность процесса резания и соответствует необходимому значению продольной подачи режущего инструмента на зуб фрезы &=0Д2 мм/зуб (табл. 1). Увеличение скорости перемещения на этапе холостого хода обеспечит сокращение общего времени обработки.
ь
с
SF0RCEJ
Рис. 4. Компьютерная модель механизма станка в среде MSC Adams с наложенными связями и движением
Рис. 5. Результаты кинематического анализа в среде MSC Adams
и.ь-
0 0
S -
и* -0.5-
-1.0-
-1
-1.5
-1.0
-0.5
0.0
Cx, м
0.5
10
1.5
Рис. 6. График движения режущего инструмента
В результате динамического анализа механизма были найдены сил моментные характеристики: момент сопротивления M, который необходимо преодолеть для того, чтобы поддерживать постоянную скорость вращения n на валу кривошипа (рис. 7), и динамические реакции в шарнирах опор (рис. 8). Как видно, пиковые значения реакций наблюдаются при угле поворота кривошипа 180° и 360°, а их уровень свидетельствует, скорее всего, об ударном характере воз-
действия. Полученные результаты расчета момента сопротивления позволяют осуществить выбор электропривода и редуктора с учетом необходимого запаса по мощности, а динамические реакции в шарнирах - выбор подшипников.
Для уточнения влияния этого явления на износ подшипника и его опоры необходимо выполнить моделирование динамики механизма в пространственном представлении, т.е. с интеграцией CAD модели станка в среду MSC.Adams и использовании модуля Adams.Machinery.
300.0 - S ¡X
250 0 -
к es s
s s СП dj 200.0 -150.0 x О с^ со s É-
о &
СО е^ S 100.0 - с о о
О É-
х о
50.0 s о
s
О 0 -
0.0 -250.0
-1000.0 -
90.0
\ I I
2 /А i I /Д
\ / / * i / / i
) / / / / /
/ / \ / /
/ / \ / /
\
■К.___ / \ ¿и
ï 3 'i \ ï |L
270.0 360.0 450.0
Угол поворота, (р
со
H
is se
&
н о
se
s s
5
4
о о
630.0
с -û É-
X
о
Рис. 7. Результаты моделирования механизма в среде MSC Adams: 1 - сила резания; 2 - скорость подачи инструмента; 3 - момент сопротивления
CD
Л
К О
CD
О
¿о
100000.0 80000.0 60000 О 4С000.0 20000.0 О.О -20000.0 -40000.0 -60000 О -80000 О 1.0Е+-005
-ROIyiPhi — R01x:Phi
I I
I ■A— —t ■i-
/ -i-
—
*
90.0 180.0 270.0 360.0 450.0 540.0 630.0
Угол поворота, ф
К
CD
50000.0 О.О -50000.0
о
С
1.0Е+005
—RAx:Ph — RAyiPfli
_
Л Л
\
90.0 180.0 270.0 360.0 450.0 540.0
Угол поворота, ф
со
CD
Л
к
eq
CD
о
Л
С
£р -50000.0 со
-1.0Е+005
-RBxiPhi RBy:Phi
и 1.DE + 005
Э , 50000.0 g
СО CN CD ^
Л CD 00
90 0 180 0 270 0 360 0 450.0 540 0 630 0
Угол поворота, ф
CD
о
-1.0Е-К305
—R02x: Phi — -RQ2y:Phl
90 0 180 0 270 0 360 0 450.0 540 0
Угол поворота, ф
Рис. 8. Графики проекций динамических реакций в соответствующих шарнирах механизма станка
Выводы. Представленная методика кинематического и динамического анализа механизма специального модуля для ремонтной обработки опорной поверхности крупногабаритных подшипников скольжения позволяет конструктивно предусмотреть необходимые технологические параметры механической обработки, обеспечивающие заданное качество обрабатываемой поверхности.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Станок для обработки вкладышей крупногабаритных подшипников скольжения: пат. 132012 Рос. федерация: МПК7 B23D 1/20, B23D 9/00 / Дуганов В.Я., Бешевли О.Б. и др.; заявитель и патентообладатель ФГБОУ ВПО Белгородский государственныйй технологический ун-т им. В.Г. Шухова. № 2013100863/02; заяв. 09.01.13; опубл. 10.09.2013 Бюл. №25.
2. Станок для фрезерования баббитовых вкладышей крупногабаритных подшипников скольжения: пат. 169070 Рос. Федерация: МПК7 В23С 1/20 / Дуганов В.Я., Бешевли О.Б., Дуюн Т.А.; заявитель и патентообладатель ФГБОУ ВПО Белгородский государственныйй технологический ун-т им. В.Г. Шухова. № 2016136002; заяв. 06.09.2016; опубл. 02.03.2017 Бюл. №7.
3. Бешевли О.Б., Дуюн Т.А. Оптимизация технологических параметров при фрезеровании баббита в условиях ремонтной обработки опорной поверхности скольжения крупногабаритных подшипников // Вестник БГТУ им. В.Г. Шухова. 2016. №8. С. 142-148.
4. Бешевли О.Б., Дуюн Т.А. Эмпирические модели шероховатости поверхности при фрезеровании баббита // Вестник БГТУ им. В.Г. Шухова. 2016. №7. С. 122-128.
5. Лозовая С.Ю., Бешевли О.Б., Дуюн Т.А., Воробьев Н.Д. Повышение эффективности ремонтной обработки опорных узлов крупногабаритных вращающихся агрегатов // Вестник Иркутского государственного технического университета. 2015. №7. С. 60-66.
6. Qingni Y., Jian L., Weiji P., Qingyun, Y. Dynamics analysis and simulation of roll grinder components // ManufacturingTechnology. 2014. № 14
(4). Pp. 600-607.
7. Cui Z., Wen G., Jiang C. Real time stuctural analysis and optimization of high speed grinder spindle system based on ADAMS ZhongguoJixieGong-cheng // China Mechanical Engineering. 2009. № 20
(5). Pp. 518-522.
8. Li B., Hu X., Wang H.Analysis and simulation for a parallel drill point grinder - Part 2: Grinding kinematic modeling and simulation // International Journal of Advanced Manufacturing Technology, 2006. №30 (3-4). Pp. 221-226.
9. Сметанников О.Ю., Могильников К.А. Моделирование динамического поведения сложных механизмов в среде ADAMS-EASY5 // Научно-технический вестник Поволжья. 2013. № 5. С. 54-58.
10. Бушуев Д.А., Воробьев Н.Д., Рубанов В.Г. Анализ динамических нагрузок в подшипниках помольно-смесительного агрегата как объекта автоматизации // Вестник Белгородского государственного технологического университета им. В.Г. Шухова. 2014. № 2. С. 143-148.
Beshevli O.B., Bushuev D.A., Duyun T.A., Rubanov V.G.
IMITATION MODELING OF DYNAMIC CHARACTERISTICS OF THE MODULE FOR PROCESSING OF LARGE-SIZED SLIDE BEARINGS
The technology of repair processing of the antifriction abutment surface of large-sized sliding bearings with using of the developed special wing-base unit, which allows to significantly reduce the processing time, improve the quality of the in-process surface and increase the period of the overhaul cycle of large rotating aggregates is proposed. The kinematic analysis of the crank-rocking mechanism which used for the proposed unit design, which allows providing constructively the necessary technological parameters of machining, is presented. Mathematical and computer models of the attachable unit in the MSC.Adams software are constructed. The dynamics of the attachable unit computer model is studied. The results allows to make it possible to evaluate the force-torque characteristics, high-speed milling modes and speed of movement on the stage of the working stroke, when the constancy of the speed ensures the stability of the cutting process that corresponds to the required cutting tool feed rate and the increased speed of movement on the stage of idling, which leads to reducing the overall processing time of the product.
Key words: large-sized sliding bearings, repair processing, milling, mathematical model, kinematic analysis, computer model, MSC.Adams, technological parameters, cutting modes.
Бешевли Олег Борисович, аспирант кафедры технологии машиностроения. Белгородский государственный технологический университет им. В.Г. Шухова. Адрес: Россия, 308012, Белгород, ул. Костюкова, д. 46. E-mail: oleg-beshevli@yandex.ru
Бушуев Дмитрий Александрович, кандидат технических наук, старший преподаватель кафедры технической кибернетики.
Белгородский государственный технологический университет им. В.Г. Шухова. Адрес: Россия, 308012, Белгород, ул. Костюкова, д. 46. E-mail: untame@list.ru
Дуюн Татьяна Александровна, доктор технических наук, доцент, заведующий кафедрой технологии машиностроения.
Белгородский государственный технологический университет им. В.Г. Шухова. Адрес: Россия, 308012, Белгород, ул. Костюкова, д. 46. E-mail: tanduun@mail.ru
Рубанов Василий Григорьевич, доктор технических наук, профессор, заведующий кафедрой технической кибернетики.
Белгородский государственный технологический университет им. В.Г. Шухова. Адрес: Россия, 308012, Белгород, ул. Костюкова, д. 46. E-mail: untame@list.ru