Russian Federation
employee
Minsk, Belarus
UDK 69.002.5 Строительные машины, механизмы и оборудование
The article deals with issues related to the study of planetary machines. It is demonstrated that they are able to significantly intensify most of the processes applicable to them. This is due to the complication of the trajectory of the movement of the working body and the increase in inertial forces. In addition, planetary units are energy efficient. The influence of the mutual action of grinding bodies on their power characteristics in a planetary mill with external and internal running of the drums is presented. The dependence of destructive forces on the angular velocity, the degree of filling of the grinding drums, the geometric ratios of the mill itself is shown. Sufficiently simple dependences have been obtained for calculating the total speed at any point of the planetary working body of the smoothing machine. Their use makes it possible to analyze the change in speed over a wide range, depending on a number of design and technological parameters. The proposed method allows to calculate the speeds of the working bodies that make planetary motion. Using the methods of classical mechanics, using the instantaneous center of velocities, a formula is obtained for calculating the speed at any point of the mixing blade. It is established that the speed of the blade changes in time according to a sinusoidal law. The total power gradually increases with increasing angular velocity. The results obtained in this article can be applied in the design and operation of planetary units.
planetary mill, working body, smoothing, mixer, mixing, rotation, carrier
Введение. Интенсификация работы технологических машин возможна за счет усложнения траектории движения рабочих органов и обрабатываемого материала, увеличения скорости их движения и перевода стационарных процессов в нестационарные. Все это может быть реализовано в машинах планетарного типа [1, 2].
Характерной особенностью таких машин является наличие колеса, которое с помощью водила обкатывается по неподвижной поверхности. В зависимости от способа обкатки (внутренней или внешней) траектория отдельных точек колеса может описывать одну из «замечательных» кривых: гипоциклоиду или эпициклоиду [3, 4].
Изменяя геометрические соотношения основных элементов планетарного механизма, можно получать самую разнообразную форму этих кривых, а соответственно и траекторию. В дополнение к этому, варьируя скоростью вращения водила, можно в широком интервале изменять скорости и ускорения отдельных точек вращающегося колеса. Таким образом, планетарные машины и механизмы имеют широчайший диапазон регулирования параметров.
Основная часть. Машины планетарного типа имеют большие перспективы для использования в механических процессах, таких как измельчение, смешение и довольно специфическом процессе – заглаживании бетонных поверхностей [5]. Несмотря на кажущуюся простоту, эти процессы сложны по своей физической сущности и слабо изучены до настоящего времени. Все они довольно энергозатратны и сопровождаются повышенным абразивным износом рабочих органов технологического оборудования. В этой связи интенсификация процесса воздействия рабочего органа на обрабатываемую среду может привести к сокращению времени их контакта, а соответственно, к снижению энергозатрат износа.
Степень интенсификации за счет придания рабочим органам планетарного движения в указанных выше процессах проявляется по-разному. Она зависит прежде всего от максимально допустимой скорости, определяемой частотой вращения водила.
Для смесителей в режиме механического перемешивания эта скорость определяется удержанием твердых частиц на лопасти. В противном случае они начинают отбрасывать частицы от центра на периферию и наступает режим сепарации с уменьшением эффективности перемешивания. По этой причине перемешивающие агрегаты преимущественно тихоходные. Скорость лопастей при этом не превышает 10 м/с.
Подобным образом ограничивается и скорость движения рабочего органа дисковой заглаживающей машины. Она представляет собой вращающийся диск, который с определенным усилием прижимается к не затвердевшей бетонной поверхности и перемещается вдоль нее. Задача диска состоит в выравнивании плоскости поверхности и заполнении мелких раковин. Показателем качества является заглаживающая способность – это длина линии, по которой рабочий орган воздействует на каждую элементарную площадку бетонной поверхности. Естественно, что при усложнении траектории движения ее длина будет больше. Указанный факт подтверждается простым сравнением. Любая точка диска с обычным вращением оставляет на поверхности след в виде окружности, а при его планетарном вращении – гипоциклоиды или эпициклоиды. Длина последних двух кривых при всех равных параметрах в 3–4 раза больше, чем окружности [3, 4]. Однако увеличение скорости вращения диска может привести к отбрасыванию частиц бетона на соседние участки уже обработанной поверхности. Такая скорость является верхним ограничением для дисковых заглаживающих машин. Она обычно выше 10 м/с, в общем случае лучше сказать составляет несколько десятков метров в секунду. По сравнению с быстроходной планетарной мельницей такую машину можно назвать среднеходной.
Следует отметить, что перед разработчиками планетарных мельниц ставятся совершенно иные задачи по сравнению со смесителями и заглаживающими машинами. Основная из них − это повышение эффективности процесса разрушения материала. Очевидно, что скорости движения рабочих органов при этом должны быть как можно больше. Основным элементом планетарной мельницы является размольный барабан, вовлеченный в планетарное движение, в который загружаются измельчающие тела (шары) и материал, подлежащий измельчению.
В конструктивном исполнении планетарные мельниц могут быть вертикальными, горизонтальными и наклонными, с внешней и внутренней обкаткой [1, 6–8]. В процессе изучения движения шаров в планетарной мельнице (рис. 1) удалось установить соотношения угловых скоростей и углов поворота водила и барабана [1, 9]:
(1)
где ω и Ω – угловые скорости барабана и водила, рад/с; φ и ψ – углы поворота барабана и водила; k – геометрический критерий, который равен k = r/R, где r – радиус приводного элемента, м; R – радиус поверхности обкатки, м.
Знак «+» в формулах – внешняя обкатка,
«-» – внутренняя.
Рис. 1. Силовая схема планетарной мельницы с внутренней обкаткой
На измельчающее тело, находящееся внутри барабана, действует сила тяжести и две инерционные силы FА и FВА, связанные с поворотом барабана и водила соответственно:
(2)
где b − второй геометрический критерий, который определяется как b = rб/r, где rб – радиус барабана.
Расчет силы давления шара на стенки барабана по формуле (3), равной нормальной реакции
( ), показал, что она меняется циклически. Причем установлена зависимость цикличности от критерия k.
(3)
Изменение силы давлении за один цикл свидетельствует о нестационарности процесса разрушения материала, обусловленного воздействием измельчающих тел. Поскольку относительная сила давления (отнесенная к единице массы шара) – это ни что иное как его ускорение, то из расчетной зависимости очевидно, что она значительно превышает ускорение свободного падения. Таким образом, основным разрушающим силовым фактором в планетарных мельницах становятся инерционные силы. Экспериментальные исследования процесса помола в планетарной мельнице подтвердили ее высокую эффективность.
Естественно модель движения одиночного мелющего тела является довольно упрощенной. В дальнейшем было изучено движение сегмента загрузки в целом, рис. 2.
Рис. 2. Расчетная схема относительного движения мелющего тела в сегменте загрузки при внешней
обкатке барабанов
Рассмотрим относительное движение мелющих тел помольной камеры планетарной мельницы в подвижной системе координат АXY, которая вращается с угловой скоростью w, равной скорости вращения барабана [1, 9].
Уравнение относительного движения мелющего тела можно записать в следующем виде:
(4)
где m – масса мелющего тела, кг; ad – относительное ускорение мелющего тела, м/с2; G – сила тяжести, Н; FT – сила трения, Н; ,
– соответственно переносная и относительная силы инерции, действующие на мелющее тело, Н; Fс – кориолисова сила инерции, Н; Fр – сила давления, Н.
В то время, когда мелющее тело находится внутри сегмента загрузки, то на него оказывают воздействие другие тела. Последнее предлагается оценить дополнительной силой, названной силой давления Fр.
Для возможности учета воздействия других мелющих тел, выраженного силой давления Fр, обратимся к модели (рис. 3). Рассмотрим послойное движение шаров по дуге окружности. Тангенциальное усилие в момент сдвига слоев будет зависеть от давления, оказываемого в радиальном направлении на каждое мелющее тело анализируемого слоя [9–12].
Рис. 3. Расчетная схема для определения силы давления
Составив уравнения для силы давления на элементарном участке dr и проведя интегрирование, получили расчетную зависимость для нее на любом радиусе ri:
(5)
где rш – радиус шара, м; ri – текущий радиус, м; r0 – расстояние от центра барабана до свободной поверхности сегмента загрузки, м.
С учетом всех силовых факторов, входящих в уравнение (5) удалось в пределах сегмента загрузки определить зоны отрыва мелющих тел, а также зоны безотрывного движения со скольжением и без него [13]. В каждой из этих зон реализуется преобладающий способ разрушающего воздействия: удар, истирание и раздавливание. Разработана методика расчета напряжений при каждом из них. В качестве примера показано изменение ударных напряжений (рис. 4), которые значительно выше сжимающих.
Рис. 4. Изменение ударных напряжений в планетарной мельнице: 1 ‒ ω = 90 рад/с; 2 ‒ ω = 180 рад/с
Таким образом можно оценить влияние на эффективность изменения каждого из них в отдельности и провести суммирование этих разрушающих эффектов.
Специфика процесса заглаживания, геометрические размеры рабочего органа и его скоростные режимы оказывают влияние на показатели эффективности и методы исследования заглаживающей машины.
Показателем эффективности заглаживания бетонных поверхностей является заглаживающая способность [5, 14], определяемая для дискового рабочего органа по формуле:
, (6)
где vд – скорость диска, м/с; vз – скорость поступательного движения машины, м/с; R – радиус диска, м.
Для определения скорости диска, вовлеченного в планетарное движение, использован совершенно иной подход. Она определялась путем дифференцирования параметрических уравнений траекторий движения: эпициклоиды при внешней обкатке и гипоциклоиды при внутренней [15].
Методологию определения скорости рассмотрим на примере эпициклоиды (рис. 5).
Рис. 5. Эпициклоида
Параметрические уравнения эпициклоиды имеют вид:
(7)
где R ‒ радиус обкатки, м; r ‒ радиус обкатываемого круга (диска), м; rd ‒ радиус производящего круга (рабочего органа), м; φ ‒ угол поворота водила, рад.
Первая производная по каждой координате для уравнений эпициклоиды по переменной φ:
(8)
Так как скорость ‒ это производная по времени, то необходимо перейти от переменной φ к новой переменной t. По каждой из координат получим:
(9)
где Ω – угловая скорость водила, рад/с.
Используя геометрические критерии k = r/R и b = rd/r и учитывая соотношения углов (10), получим уравнение для расчета полной скорости диска (11).
(10)
(11)
На рис. 6 показано изменение модуля полной скорости в зависимости от частоты вращения водила n за один полный оборот.
Очевидно, что зависимость линейной скорости от угла поворота носит периодический характер. Причем увеличение частоты вращения приводит только к росту амплитуды изменения скорости. Изменение геометрического критерия k не очень существенно влияет на амплитуду скорости.
Соблюдая последовательность в изложении материала от быстро- и среднеходных технологических агрегатов, переходим к тихоходным. Таковыми в нашем случае являются лопастные смесители с планетарным движением рабочего органа. Задача процесса перемешивания заключается в достижении гомогенности среды по всему объему смесителя и поддержании ее в таком состоянии длительное время. Это особенно важно в высокопроизводительных технологических процессах, например, в шлам-бассейнах в производстве цемента.
Рис. 6. Зависимость скорости от угла поворота при различных n
При наличии в перемешиваемой среде нескольких компонентов эффективность процесса обычно оценивается по одному из них – ключевому. По степени его распределения судят о качестве перемешивания. Для этого в разных местах смесителя отбирают пробы и определяют долю ключевого компонента. Ее равенство во всех точках свидетельствует об идеальном перемешивании. Однако в реальных условиях это равенство не достигается. Его нарушение оценивается среднеквадратичным отклонением по ключевому компоненту. Делением среднеквадратичного отклонения на долю ключевого компонента в исходном состоянии получают коэффициент вариации, который используется в качестве показателя эффективности перемешивания. К сожалению в этом показателе совершенно не учитываются ни геометрические, ни технологические параметры самого перемешивающего агрегата. В этой связи мы были вынуждены обратиться к альтернативным источникам [16], предлагающим другие методы оценки качества перемешивания.
Эффективность перемешивания в смесителе, в том числе и планетарном [17], оценивается критерием λ, с−1, определяемым по следующей формуле:
(12)
где v – средняя скорость лопасти, м/с; S – активная площадь лопасти, м2; V − объем смесителя, м3.
Важные численные параметры, характеризующие работу смесителя, зависят от скорости движения лопастей. Для ее определения рассмотрим расчетную схему (рис. 7).
Рис. 7. Схема расположения лопастей
Если предположить, что один перемешивающий орган содержит три лопасти, то каждая из них жестко скреплена с колесом стержнем АС. При этом за один цикл водила ОА все лопасти совершают одинаковые движения. В этой связи расчет скорости, а затем и мощности можно выполнять для одной лопасти с последующим увеличением на их количество.
При известной угловой скорости водила ω в зафиксированном положении механизма φ = ω · t определим скорость в произвольной точке K лопасти. Примем, что в начальный момент движения φ = 0, стержень АС3 совпадал с осью Ох. После поворота водила на угол φ стержень АС3 повернется по отношению к водилу на угол φ1. Тогда формула примет вид:
(13)
где k – геометрический критерий, равный k = r1/r2.
Угол поворота колеса вместе с лопастями в абсолютном движении равен углу поворота вокруг мгновенного центра скоростей E и находится согласно выражениям:
(14)
Аналогичная связь будет между угловыми скоростями водила и колеса, которая находится по формуле:
(15)
С учетом этого скорость в точке K определяется по формуле:
(16)
где lEK – длина линии EK.
В результате геометрических преобразований получена формула длины lEK:
(17)
Для упрощения дальнейших преобразований переменную lСK обозначим через aK: lСK = aK. По расчетной схеме (рис. 7) видно, что lЕА = r2. Неизменную длину стержня lАС обозначим lАС = l.
После введения новых обозначений и преобразования формулы (17) для lEK, скорость лопасти в произвольной точке примет вид:
(18)
При фиксированных значениях угловой скорости водила ω = 1 рад/с и геометрическом критерии k = 3 по формуле (18) произведем расчет и анализ изменения скорости по ширине лопасти за один полный оборот водила (φ = 360°).
Результаты расчета представлены в виде трехмерной графической зависимости (рис. 8).
Рис. 8. Изменение скорости по ширине лопасти за один полный оборот водила
Из графика видно, что скорость vK равномерно повышается с увеличением ширины лопасти (aK). А вот в процессе вращения (угол φ) изменение скорости носит синусоидальный характер. Причем, экстремальные значения на соседних пиках графика разные. Максимальные значения скорости достигают vK = 2,74 м/с.
В дальнейшем был проведен расчет мощности, затрачиваемой на перемешивание, которая напрямую зависит от скорости перемещения лопастей. Расчет показал, что, несмотря на локальные пульсации скорости, мощность в зависимости от угловой скорости изменяется по линейному закону.
Выводы. Подводя итог, можно отметить, что в работе представлены методы расчета кинематических и динамических характеристик трех разных видов технологических машин планетарного типа. Причем как конструкции агрегатов, так и методики расчета имеют существенные отличия. Это свидетельствует о том, что процесс их исследования находится в постоянном развитии и у него есть хорошие перспективы.
1. Vaytekhovich P.E. Intensification and modeling of dispersion processes in the field of inertial forces [Intensifikaciya i modelirovanie processov dispergirovaniya v pole inercionnyh sil]. Minsk: BSTU. 2008. 220 p. (rus)
2. Bushuev L.P Multimode planetary mill [Mnogorezhimnaya planetarnaya mel'nicy]. Izv. universities. Mining magazine. 1965. No 10. Pp 148-154. (rus)
3. Zvavich L.I., Ryazanovskij A.R. Geometry. Reference manua. [Geometriya. Spravochnoe posobie]. M.: Drofa. 2012. 128 p. (rus)
4. Markov V.G. Planetary mechanisms: textbook. manual for the course "Theory of mechanisms and machines" [Planetarnye mekhanizmy: ucheb. posobie po kursu «Teoriya mekhanizmov i mashin»]. L.: LTHP. 1969. p. 28. (rus)
5. Bolotnyj A.V. Theory and processes of smoothing concrete surfaces [Teoriya i processy zaglazhivaniya betonnyh poverhnostej]. L: Stroizdat. 1975. 127 p. (rus)
6. Kim B.G. Relative motion of grinding bodies in the drum of a planetary mill in the waterfall mode of its operation [Otnositel'noe dvizhenie melyushchih tel v barabane planetarnoj mel'nicy pri vodopadnom rezhime ee raboty]. Izv. universities. Mining magazine. 1975. No. 9. Pp. 81-86. (rus)
7. Bushuev L.P. Experimental research and questions of the theory of planetary centrifugal mills [Eksperimental'nye issledovaniya i voprosy teorii planetarnyh centrobezhnyh mel'nic]. Scientific report high school. Mining. 1959. No 2. Pp. 220-226. (rus)
8. Dobroborskij G.A., Lyansberg L.M., Rabin A.N. Determination of the main modes of movement of the load in the drums of a multi-drum planetary centrifugal mill with vertical axes [Opredelenie osnovnyh rezhimov dvizheniya zagruzki v barabanah mnogobarabannoj planetarno-centrobezhnoj mel'nicy s vertikal'nymi osyami]. Izv. universities. Mining magazine. 1993. No 1. Pp. 85-89. (rus)
9. Semenenko D.V. Influence of design and technological parameters of a horizontal planetary mill on the efficiency of the grinding process: abstract of the thesis. dis. … cand. tech. sciences; BSTU. Minsk. 2014. 20 p. (rus)
10. Zubko I. YU., Permyakov P.V., Trusov P.V. Dynamics of many colliding bodies: application to mechanical alloying [Dinamika mnogih soudaryayushchihsya tel: prilozhenie k mekhanicheskomu legirovaniyu]. Izv. universities. North Caucasian region. Natural Sciences. 2008. No. 2. Pp. 29-35. (rus)
11. Zubko I.YU., Zajcev A.V. Simulation of milling balls motion within mechanical alloying process [Modelirovanie dvizheniya melyushchih tel v processe mekhanicheskogo legirovaniya]. Bulletin of PNRPU. Aerospace engineering. 2017. No. 2. Pp. 34-43.
12. Feng Y.T., Han K., Owen D.R.J. Discrete element simulation of the dynamics of high energy planetary ball milling processes // Materials Science and Engineering A. 2004. Vol. 375-377. Pp. 815-819.
13. Vaytekhovich P.E., Borovskiy D.N. Destructive effect of grinding media in a horizontal planetary mill with different methods of running in grinding drums [Razrushayushchee vozdejstvie melyushchih tel v gorizontal'noj planetarnoj mel'nice pri razlichnyh sposobah obkatki pomol'nyh barabanov]. Chemical and petroleum engineering. 2019. No 2. Pp. 3-6. (rus)
14. Ryss-Berezark S.A. Determination of the parameters of smoothing machines with a disk working body with a simple and complex movement [Opredelenie parametrov zaglazhivayushchih mashin s diskovym rabochim organom s prostym i slozhnym dvizheniem]. L: Stroizdat. 1999. 167 p. (rus)
15. Vajtekhovich P.E., Sidorov N.N. Determination of the speed of movement of technological machines of planetary type [Opredelenie skorosti dvizheniya tekhnologichnyh mashin planetarnogo tipa]. Trudy BGTU. Ser. himii i tekhnologii neoorgan. 2004. No 12. Pp. 28-35.
16. Vajtekhovich P.E., Franckevich V.S., Hrebianchuk P.S. Machinery and equipment for construction materials companies: a textbook [Mashyny і abstalyavanne pradpryemstvaў budaўnіchyh materyyalaў: vucheb. dapam.]. Minsk: BSTU. 2018. 297 p. (rus)
17. Vajtekhovich P.E., Hves'ko G.M., Borovskij D.N., Semenenko D.V. Determination of the main parameters of an efficient planetary mixer [Opredelenie osnovnyh parametrov effektivnogo planetarnogo smesitelya]. Trudy BGTU. Series 2. Chemical technologies, biotechnology and geoecology. 2019. No. 2 (223). Pp. 114-119. (rus)
