Россия
сотрудник
Минск, Беларусь
УДК 69.002.5 Строительные машины, механизмы и оборудование
В статье рассмотрены вопросы, связанные с исследованием планетарных машин. Показано, что они способны значительно интенсифицировать большинство применимых к ним процессов. Это происходит за счет усложнения траектории движения рабочего органа и увеличения инерционных сил. Кроме того, агрегаты планетарного типа являются энергоэффективными. Представлено влияние взаимного воздействия мелющих тел на их силовые характеристики в планетарной мельнице с внешней и внутренней обкатками барабанов. Показана зависимость разрушающих усилий от угловой скорости, степени заполнения помольных барабанов, геометрических соотношений самой мельницы. Получены достаточно простые зависимости для расчета полной скорости в любой точке планетарного рабочего органа заглаживающей машины. Их использование дает возможность проанализировать изменение скорости в широких пределах в зависимости от ряда конструктивных и технологических параметров. Предложенный метод позволяет выполнять расчет скоростей рабочих органов, совершающих планетарное движение. Методами классической механики с использованием мгновенного центра скоростей получена формула для расчета скорости в любой точке перемешивающей лопасти. Установлено, что скорость лопасти во времени изменяется по синусоидальному закону. Общая мощность с увеличением угловой скорости плавно возрастает. Полученные результаты в данной статье могут быть применены при конструировании и эксплуатации планетарных агрегатов.
планетарная мельница, рабочий орган, заглаживание, смеситель, перемешивание, вращение, водило
Введение. Интенсификация работы технологических машин возможна за счет усложнения траектории движения рабочих органов и обрабатываемого материала, увеличения скорости их движения и перевода стационарных процессов в нестационарные. Все это может быть реализовано в машинах планетарного типа [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. Вайтехович П.Е. Интенсификация и моделирование процессов диспергирования в поле инерционных сил. Минск: БГТУ. 2008. 220 с.
2. Бушуев Л.П. Многорежимная планетарная мельницы // Изв. вузов. Горный журнал. 1965. № 10. С. 148-154.
3. Звавич Л.И., Рязановский А.Р. Геометрия. Справочное пособие. Москва. 2012.128 с.
4. Марков В.Г. Планетарные механизмы: учеб. пособие по курсу «Теория механизмов и машин». Л.: ЛТХП. 1969. 28 с.
5. Болотный А.В. Теория и процессы заглаживания бетонных поверхностей. Л.: Строиздат. 1975. 127 с.
6. Ким Б.Г. Относительное движение мелющих тел в барабане планетарной мельницы при водопадном режиме ее работы // Изв. вузов. Горный журнал. 1975. № 9. С. 81-86.
7. Бушуев Л.П. Экспериментальные исследования и вопросы теории планетарных центробежных мельниц // Науч. докл. высшей школы. Горное дело. 1959. № 2. С. 220-226.
8. Доброборский Г.А., Лянсберг Л.М., Рабин А.Н. Определение основных режимов движения загрузки в барабанах многобарабанной планетарно-центробежной мельницы с вертикальными осями // Изв. вузов. Горный журнал. 1993. № 1. С. 85-89.
9. Семененко Д.В. Влияние конструктивных и технологических параметров горизонтальной планетарной мельницы на эффективность процесса измельчения: автореф. дис. … канд. техн. наук; Белорус. гос. технол. ун-т. Минск. 2014. 20 с.
10. Зубко И.Ю., Пермяков П.В., Трусов П.В. Динамика многих соударяющихся тел: приложение к механическому легированию // Известия вузов. Северо-кавказский регион. Естественные науки. 2008. № 2. С. 29-35.
11. Зубко И.Ю., Зайцев А.В. Моделирование движения мелющих тел в процессе механического легирования // Вестник ПНИПУ. Аэрокосмические техника. 2017. № 2. С. 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. Р. 815-819.
13. Вайтехович П.Е., Боровский Д.Н. Разрушающее воздействие мелющих тел в горизонтальной планетарной мельнице при различных способах обкатки помольных барабанов // Химическое и нефтегазовое машиностроение. 2019.№ 2. С. 3-6.
14. Рысс-Березарк С.А. Определение параметров заглаживающих машин с дисковым рабочим органом с простым и сложным движением. Л.: Строиздат. 1999. 167 с.
15. Вайтехович П.Е., Сидоров Н.Н. Определение скорости движения технологичных машин планетарного типа // Труды БГТУ. Сер. химии и технологии неорган. Вып. 12. 2004. С. 28−35.
16. Вайцяховіч П.Я., Францкевіч В.С., Грэбянчук П.С. Машыны і абсталяванне прадпрыемстваў будаўнічых матэрыялаў: вучэб. дапам. Мінск: БДТУ. 2018. 297 с.
17. Вайтехович П.Е., Хвесько Г.М., Боровский Д.Н., Семененко Д.В. Определение основных параметров эффективного планетарного смесителя / Труды БГТУ. Серия 2. Химические технологии, биотехнология и геоэкология. 2019. №2(223). С. 114-119.
