УДК 621.926
П. Е. Вайтехович, Д. Н. Боровский, Е. О. Ильяшенко, С. А. Герасимчик
Белорусский государственный технологический университет
МЕТОД РАСЧЕТА РАДИАЛЬНЫХ И ОКРУЖНЫХ НАГРУЗОК ПРИ ЗУБЧАТОМ ПРИВОДЕ ПЛАНЕТАРНОЙ МЕЛЬНИЦЫ
Рассмотрены вопросы, связанные с определением динамических нагрузок, возникающих в подвижных узлах планетарной мельницы с шестеренчатым приводом. Выполнен анализ влияния конструктивных и технологических параметров на их величину. Для проведения исследования была выбрана горизонтальная планетарная мельница с внешней обкаткой помольных барабанов. Так как на характер движения загрузки в помольных барабанах планетарной мельницы большое влияние оказывают геометрические параметры отдельных узлов агрегата, были введены геометрические критерии, отражающие соотношение между этими параметрами. В ходе силового расчета привода планетарной мельницы определены инерционные силы, оказывающие существенное влияние на движение. С учетом принципа Даламбера составлены уравнения кинетостатики. Из этих уравнений получены выражения для определения окружной составляющей силы зацепления и реакции на оси барабана. Апробация методики расчета этих сил проводилась на планетарной мельнице с размерами реального промышленного объекта. В ходе расчета варьировались такие параметры, как угловая скорость водила, соотношение геометрических параметров и степень загрузки барабанов. Изменения искомых параметров анализировались на протяжении одного полного оборота водила. Произведя оценку полученных результатов, можно отметить, что изменение опорной реакции и окружной силы происходит по косинусоидальному закону в небольшом диапазоне. Также установлено, что увеличение угловой скорости и коэффициента загрузки приводит к значительному повышению динамических нагрузок. Общим итогом аналитических исследований является разработка метода расчета динамических нагрузок на опорные узлы планетарной мельницы с учетом влияния максимального количества параметров.
Ключевые слова: планетарная мельница, частота вращения, загрузка мельницы, инерционные силы, зубчатый привод, динамические нагрузки, реакция опоры, окружная сила.
P. Ye. Vaitekhovich, D. N. Borowskiy, E. O. Ilyashenko, S. A. Gerasimchik
Belarusian State Technological University
THE METHOD OF CALCULATION OF RADIAL AND CIRCUMFERENTIAL LOADS IN A GEAR DRIVE OF PLANETARY MILL
The issues associated with the determination of dynamic loads encountered in moving the nodes of the planetary mill with gear drive. The analysis of influence of constructive and technological parameters on their size. To conduct the study was selected horizontal planetary mill with external running of the grinding drums. As traffic load in the grinding drums, planetary ball mills are largely influenced by geometric parameters of individual units was introduced geometric criteria reflecting the correlation between these parameters. Ho de the power calculation drive planetary mill was determined as the inertial force of an eye-Sabaudia a significant impact on the movement. Taking into account the principle of d'alembert the equations terms of kinetostatical. From these equations the expressions for determining the circumferential force component of engagement and reaction on the axis of the drum. Approbation of the methods of calculation of these forces was carried out in a planetary mill with dimensions of a real industrial object. In the course of computation varied parameters such as the angular velocity of the carrier, the ratio of geometric parameters and the degree of loading of the reels. Change the desired parameters of the analysis were zeroality for one full revolution of the carrier. Having evaluated the received re-of results, it can be noted that the change in support reactions and the circumferential force occurs at the cosine law in a small range. It is also established that the increase of angular velocity and load factor leads to a significant increase of dynamic loads. The overall outcome of the analytical studies is to develop a method of calculating the dynamic loads on the bearings of the planetary mill with the influence of the maximum number of parameters.
Key words: planetary mill speed, mill load, inertial forces, gear drive, dynamic loads, reaction support, ring power.
Введение. Более десятка лет авторы данной работы занимаются исследованием планетарных мельниц, что подтверждается монографией [1]
одного из них, а также многочисленными публикациями [2-5]. Вначале исследовалось движение одиночного мелющего тела [1-2]. Были
определены режимы движения, установлено влияние конструктивных и технологических параметров на величину инерционных сил, действующих на тело. В дальнейшем учтено влияние конструктивных и технологических параметров на характер перемещения элементов загрузки [3-4]. Исследование кинематики и динамики движения мелющих тел позволило определить весь комплекс нагрузок, действующих на измельчаемый материал с их стороны. А это в свою очередь обеспечило возможность приступить к составлению модели разрушения материала в планетарной мельнице [5]. Такая модель была разработана, апробирована и показала хорошую сходимость расчетных и экспериментальных данных.
После проведенных аналитических и экспериментальных исследований встал вопрос о создании промышленного образца планетарной мельницы. На этом этапе наиболее важным является выбор оптимальной схемы привода и разработка методов его расчета. Работа в этом направлении уже начата. Был предложен метод определения динамических нагрузок при ременном приводе [6] и дан анализ влияния различных факторов на их величину. Очевидно, что невысокий коэффициент полезного действия ременных передач всегда ограничивает возможность их использования. Кроме того, у планетарных мельниц нагрузки изменяются в широком диапазоне, даже на протяжении одного оборота. Это в еще большей степени снижает надежность и так не долговечного ременного привода.
Основная часть. Указанные выше недостатки переориентировали авторов на исследование зубчатого привода планетарных мельниц. Да, он более сложен в изготовлении и эксплуатации, но в нем полностью устранены проблемы ременного привода.
Рис. 1. Расчетная схема
Для исследования принята зубчатая передача, представляющая собой шестерню, установ-
ленную соосно с размольным барабаном и обкатываемую по наружной поверхности неподвижного зубчатого колеса (рис. 1). Обкатка осуществляется с помощью водила, соединяющего оси барабана и неподвижного колеса.
Согласно схеме на рис. 1, с использованием принципа Даламбера составлены уравнения кинетостатики:
Ж = 0,
k=1
ФА + Фс + ФГ • sina-(GA + GC)• sinф-
-FV t§aW - XA = 0
(1)
Ж = o.
k=1
-Фг • cosa-(GA + GC )• cosф-FT + YA = 0,
¿m (F )=0,
k=1
ФС • 1ас • COs a - gc • lac • sin Ф+ Fx • Га = 0.
где 0А и ОС - вес барабана и загрузки соответственно; Ф^ - инерционная сила от массы барабана; ФС и Фг - переносная и относительная силы инерции загрузки соответственно; ¥г и ¥х -окружная и радиальная силы соответственно; ХА и УА - проекции реакций связи на соответствующие оси координат; а - угол, на который сегмент загрузки отстает от водила; ф - угол поворота водила; ю и & - угловые скорости барабана и водила соответственно; аш - угол зацепления; 1АС - расстояние от центра барабана до центра тяжести загрузки; гА - радиус шестерни.
Стандартный угол для эвольвентного зацепления аш = 20° [7].
Расстояние 1АС рассчитывается по формуле [8]:
lac =
4 • гб • (sin у)3 3 •( 2 •у- sin2y)'
(2)
где у - половина угла сектора, который ограничивает сегмент загрузки.
Угол а, определяющий положение загрузки в барабане мельницы, составляет 40-60° [9].
Из системы уравнений (1) можно определить силы в зацеплении: - окружную
Ft = — •(Ge • Iас • sin ф-Фc • 1ас • cos а); (3)
радиальную:
F = F • tga
w-
(4)
а
Эти силы используются при расчете зубчатых колес на усталость.
По уравнениям (1) определяются составляющие нормальной реакции связи:
Xa = Ф л
ФС + Ф г
• sin а -
-(ga + gc )•sin ф-t§aw;
(5)
Ya = Фг • cos a + ( Ga + Gc )• cos ф + FT. (6)
Искомой величиной будет опорная реакция на оси барабана:
ra =4xa+ya.
(7)
По величине этой реакции вычисляется динамическая грузоподъемность и подбираются подшипники для опорных узлов.
Из приведенных формул (3)-(7) видно, что их использование возможно при предварительном расчете сил тяжести Ол и ОС и инерционных сил ФА, ФС и Фг. Силы инерции определяются по следующим формулам:
фл = ma • aa = ma
• п2 •( + гл);
Фс = тз • п2
•( гп + гл);
Ф = т • ю2 • l
AC-
(8) (9)
(10)
Причем масса барабана рассчитывалась как суммарная масса боковых стенок с крышками и шестернями. При расчете массы загрузки тз учитывалась доля стальных шаров и измельчаемого материала. Соотношение угловых скоростей барабана ю и водила О находилось по зависимости для внешней обкатки, полученной нами ранее [1]:
1 + к
(11)
Описанный выше алгоритм был реализован при расчете промышленного образца планетарной мельницы с диаметром барабана 0,2 м (гб = 0,1 м) и длиной 1 м. В качестве варьируемых параметров принимались угловая скорость водила О = 30-80 рад/с и угол его поворота ф = 0-360°. Основным расчетным параметром являлась опорная реакция ЯА, дополнительным - окружная сила ¥х. Влияние степени загрузки помольных барабанов на характер изменения опорной реакции оценивали коэффициентом загрузки Кз:
К =
S= 0,5 • г2 • (2 •у-sin у)
о 2 :
п • г
(12)
где & и ^б - площадь сегмента и поперечного сечения барабана соответственно.
Возможные изменения геометрических размеров привода учитывались геометрическими критериями. Причем в данной работе при их назначении мы придерживались того же принципа, что и при изучении механики мелющей загрузки [1-4].
В результате первый из критериев, характеризующий соотношение размеров зубчатых колес, представлен в виде
к = -
loa га
а Г
(13)
Второй геометрический критерий устанавливает соотношение между размерами барабана и приводной шестерни:
b = —
(14)
Результаты расчетов показали, что при фиксированных значениях геометрических критериев к = 0,4, Ь = 1,0, постоянных коэффициенте загрузки Кз = 0,5 и скорости водила О = 30 рад/с реакция на протяжении одного оборота водила изменяется по косинусоидальному закону, но амплитуда ее изменения, как и при ременном приводе [6], не очень велика. Максимальное значение ЯА = 614 кН достигается в диапазоне изменения угла поворота водила ф = 270-300° (рис. 2).
а
е р
н р
о п О
620 615 610 605 600 595 590
585
0 30 60 90 120 150 180 210 240 270 300 330 360 Угол поворота водила ф, град
Рис. 2. Зависимость изменения величины опорной реакции от угла поворота водила
Аналогичным образом изменяется реакция связи и при других значениях угловой скорости. Но ее увеличение приводит к значительному повышению величины опорной реакции (рис. 3).
^ 5000 с^ 4000 | 3000
I 2000 « 1000 0
н р
о п О
30 50 60 80
Угловая скорость водила О, рад/с
Рис. 3. Зависимость изменения величины опорной реакции от угловой скорости водила
г
а
К такому же повышению приводит и увеличение степени загрузки барабана до Кз < 0,5 (рис. 4). При дальнейшем ее увеличении значение опорной реакции стабилизируется прежде всего за счет уменьшения расстояния 1АС.
л <и а
к а о а О
3000 2500 2000 1500 1000 500 0
0,0288 0,196 0,500 0,804 Коэффициент заполнения К,
Рис. 4. Зависимость изменения величины опорной реакции от коэффициента заполнения
Следует отметить, что оба предыдущих графика построены по максимальным значениям нормальной реакции при угле поворота водила ф = 300°.
Повышение величины опорной реакции в зависимости от двух проанализированных параметров вполне объяснимо. Во-первых, увеличивается масса загрузки, а соответственно, сила тяжести ОС и инерционные силы ФС и Фг, во-вторых, все инерционные силы значительно возрастают при повышении угловой скорости.
Что касается влияния геометрических критериев на величину опорной реакции, то здесь закономерности несколько иные.
При фиксированном значении критерия Ь = 1,0, что соответствует равенству диаметров барабана и приводной шестерни, повышение критерия к возможно за счет уменьшения диаметра подвижного колеса, т. е. радиуса гп. При этом уменьшается длина водила (гп + гА). Естественно, что снижаются значения инерционных сил ФС и ФА, которые существенно влияют на величину опорной реакции, уменьшая ее. Правда, одновременно с этим, согласно формуле (11), повышается угловая скорость барабана ю и вместе с ней инерционная сила Фг. Но при малых значениях критерия к < 0,2 его влияние на угловую скорость барабана и указанную инерционную силу не очень существенно. По этой причине опорная реакция в указанном диапазоне резко снижается (рис. 5).
В дальнейшем это влияние усиливается и величина опорной реакции стабилизируется. Увеличение геометрического критерия Ь при фиксированных диаметре барабана (гб) и критерии к = 0,4 равнозначно уменьшению диаметра приводной шестерни (гА), а соответст-
венно, и длины водила (гп + гА). А все это в совокупности приводит к почти линейному уменьшению опорной реакции (рис. 6).
3 5000
<=сГ 4000
ция 3000
цак 2000
еа
р 1000
ая 0
на
р
о п
о
Л
к:
0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0
Геометрический критерий к
Рис. 5. Зависимость изменения величины опорной реакции от величины
геометрического критерия к
900 800 700 600 500 400 300 200 100
н р
рпо 0
0,6
0,7 0,8 0,9 1,0 1,1 1,2 1,3 1,4 Геометрический критерий Ь
Рис. 6. Зависимость изменения величины опорной реакции от величины геометрического критерия Ь
Кроме того, была рассчитана тангенциальная сила и построена ее зависимость от угла поворота водила на протяжении полного оборота (рис. 7).
24,0 23,5 23,0 22,5 22,0 21,5 21,0
и О
10 30 50 7 0 90 110 130 150 170 190 210 230 250 270 290 310 330 350
Угол поворота водила ф, град
Рис. 7. Зависимость изменения величины окружной силы от угла поворота водила
По полученной зависимости можно сказать, что характер изменения тангенциальной силы аналогичен изменению опорной реакции. Однако ее величина значительно меньше.
Заключение. В данной работе предложен алгоритм расчета зубчатого привода планетарной мельницы на основе расчетной схемы, в которой учтены силы тяжести размольного барабана, приводной шестерни и загрузки, а также силы инерции. В качестве определяемых параметров приняты опорная реакция на оси барабана и окружная сила в зубчатом зацеплении. Для определения этих сил составлена система уравнений кинетостатики, в результате решения которой получены расчетные зависимости.
Реализация составленного алгоритма позволила провести анализ изменения опорной реак-
ции и окружной составляющей силы зацепления от угла поворота водила, его угловой скорости, степени заполнения размольного барабана и ряда геометрических параметров. Опорная реакция является основным параметром для расчета динамической грузоподъемности при выборе подшипников, а окружная сила - при расчете на выносливость зубчатого зацепления. Методика расчета прошла апробацию на реальном промышленном объекте. Ее можно рекомендовать для оптимизации конструктивных элементов шестеренчатого привода планетарных мельниц.
Литература
1. Вайтехович П. Е. Интенсификация и моделирование процессов диспергирования в поле инерционных сил. Минск: БГТУ, 2008. 220 с.
2. Вайтехович П. Е., Семененко Д. В., Юхневич Д. В. Специфика движения мелющих тел в вертикальной планетарной мельнице // Химическое и нефтегазовое машиностроение. 2009. С. 7-10.
3. Вайтехович П. Е., Семененко Д. В., Хвесько Г. М. Анализ моделей относительного движения мелющих тел в планетарной мельнице // Труды БГТУ. Сер. III, Химия и технология неорган. в-в. 2010. Вып. XVIII. С. 167-171.
4. Вайтехович П. Е.,Семененко Д. В. Влияние взаимодействия между мелющими телами на характер их движения в планетарной мельнице // Химическое и нефтегазовое машиностроение. 2010. № 9.С. 13-15.
5. Вайтехович П. Е., Хвесько Г. М., Козлов П. С. Эффективность разрушающего воздействия мелющих тел в горизонтальной планетарной мельнице // Труды БГТУ. 2013. № 3: Химия и технология неорган. в-в. С. 182-187.
6. Вайтехович П. Е., Хвесько Г. М., Боровский Д. Н. Определение динамических нагрузок на опорные узлы при ременном приводе планетарной мельницы // Труды БГТУ. 2016. № 3: Химия и технология неорган. в-в. С. 121-126.
7. Детали машин и основы конструирования / А. Ф. Дулевич [и др.]. Минск: БГТУ, 2006. 219 с.
8. Звавич Л. И., Рязановский А. Р. Геометрия: справ. пособие. М.: Дрофа, 2002. 128 с.
9. Ким Б. Г. Определение положения помольной среды в планетарной мельнице при каскадном режиме // Изв. вузов. Горный журнал. 1976. № 1. С. 149-153.
References
1. Vaytekhovich P. Ye. Intensifikatsiya i modelirovaniye protsessov dispergirovaniya v pole iner-tsionnykh sil [Intensification and modeling of dispersion in the field of inertial forces]. Minsk, BGTU Publ., 2008. 220 p.
2. Vaytekhovich P. Ye., Semenenko D. V., Yukhnevich D. V. Specifics of the motion of grinding bodies in vertical planetary stranded-nice. Khimicheskoye i neftegazovoye mashinostroeniye [Chemical and petroleum engineering], 2009, pp. 7-10 (In Russian).
3. Vaytekhovich P. Ye., Semenenko D. V., Khves'ko G. M. The analysis models the relative motion of grinding bodies in a planetary mill. Trudy BGTU [Proceedings of BSTU], serias III, Chemistry and Technology of Inorganic Substances, 2010, issue XVIII, pp. 167-171 (In Russian).
4. Vaytekhovich P. Ye., Semenenko D. V. Effect of interaction between grinding bodies the character of their motion in a planetary mill. Khimicheskoye i neftegazovoye mashinostroeniye [Chemical and petroleum engineering], 2010, no. 9, pp. 13-15 (In Russian).
5. Vaytekhovich P. Ye., Khves'ko G. M., Kozlov P. S. The effectiveness of the damaging effects of grinding bodies in a horizontal rest a planetary mill. Trudy BGTU [Proceedings of BSTU], 2013, no. 3: Chemistry and Technology of Inorganic Substances, pp. 182-187 (In Russian).
6. Vaytekhovich P. Ye., Khves'ko G. M., Borovskiy D. N. Determination of dynamic loads on the bearings in the belt drive planetary mill. Trudy BGTU [Proceedings of BSTU], 2016, no. 3: Chemistry and Technology of Inorganic Substances, pp. 121-126 (In Russian).
7. Dulevich A. F., Osoko S. A., Los' A. M., Tsaruk F. F. Detali mashin i osnovy konstruirovaniya [Machine parts and design principles]. Minsk, BGTU Publ., 2006. 219 p.
8. Zvavich L. I., Ryazanovskiy A. R. Geometriya [Geometry]. Moscow, Drofa Publ., 2012. 128 p.
9. Kim B. G. Determination of the position of the grinding environment in a planetary mill at CAS-cascade mode. Izvestiya vuzov. Gornyy zhurnal [Higher education. Mining journal], 1976, no. 1, pp. 149— 153 (In Russian).
Информация об авторах
Вайтехович Петр Евгеньевич - доктор технических наук, профессор, заведующий кафедрой машин и аппаратов химических и силикатных производств. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: vpe51@mail.ru
Боровский Денис Николаевич - кандидат технических наук, ассистент кафедры машин и аппаратов химических и силикатных производств. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: dzianis23@rambler.ru
Ильяшенко Егор Олегович - магистр технических наук, аспирант кафедры машин и аппаратов химических и силикатных производств. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: ilyashenko@tut.by
Герасимчик Сергей Анатольевич - студент 5-го курса факультета химической технологии и техники. Белорусский государственный технологический университет (220006, г. Минск, ул. Свердлова, 13а, Республика Беларусь). E-mail: elf01@tut.by
Information about the authors
Vaytekhovich Petr Yevgen'yevich - DSc (Engineering), Professor, Head of the Department of Machines and Apparatus for Chemical and Silicate Production. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: vpe51@mail.ru
Borovskiy Denis Nikolaevich - PhD (Engineering), assistant lecturer, the Department of Machines and Apparatus for Chemical and Silicate Production. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: dzi-anis23@rambler.ru
Ilyashenko Egor Olegovich - Master of Engineering, PhD student, the Department of Machines and Apparatus for Chemical and Silicate Production. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: ilyashenko@tut.by
Gerasimchik Sergey Anatol'yevich - student of 5th course of Chemical Technology and Engineering Faculty. Belarusian State Technological University (13a, Sverdlova str., 220006, Minsk, Republic of Belarus). E-mail: elf01@tut.by
Поступила 15.11.2016