Математические г<
УДК 51-72
стрУктуры и^модмирование DOI 10.24147/2222-8772.2023.1.4-14
zuzö. №1 (65). С. 4-14
СКОЛЬЖЕНИЕ ЧАСТИЦ ГРУНТА ПО РАБОЧИМ ОРГАНАМ ВИБРОМАШИН
А.В. Архипенко
канд. техн. наук, доцент, e-mail: iten-miu@mail.ru
Международный инновационный университет, Сочи, Россия
Аннотация. Показано, что в режимах движения «без подбрасывания» коэффициент трения материала о грузонесущую поверхность рабочего органа и угол вибраций оказывают существенное влияние на эффективность вибротранспортирования. Коэффициент передачи горизонтальной составляющей скорости вибраций при постоянном значении коэффициента режима работы сначала возрастает с увеличением коэффициента трения, достигая максимума при определённом соотношении коэффициента трения и угла вибраций, а затем плавно снижается.
Ключевые слова: скольжение частиц, грунт, вибромашины, вибротранспортирование.
На технологических и санитарно-гигиенических показателях производственных процессов, связанных с внутрицеховым перемещением сыпучих материалов, негативным образом сказывается сильный износ движущихся частей исполнительных механизмов, значительные потери транспортируемых продуктов в виде пыли и газов, выделение агрессивных газов и др.
В связи с этим в последнее время всё более широко проводятся изыскания в области создания эффективных транспортирующих, транспортно-подающих и транспортно-технологических машин, использующих принцип вибрационного транспортирования [1-3]. Это обусловлено: возможностью использования герметичных грузонесущих органов для транспортирования газо- и пылеоб-разующих продуктов; возможностью транспортирования горячих (до 800 °С) материалов [4, 5]; простотой включения в автоматизированные технологические линии [4]; сравнительно невысокими затратами на содержание и ремонт; низкой энергоёмкостью и возможностью совмещения транспортирования материала с различными технологическими операциями, такими как охлаждение, нагрев, сушка, возгонка, грохочение, обезвоживание и т. п. [5-11]. Кроме того, вибрационное транспортирование является более «чистым» по экологическим требованиям, так как позволяет транспортировать мелкодисперсные пылевидные материалы по герметичной вибротрубе.
Вибротранспортирующие машины способны заменить такие традиционные виды непрерывного транспорта, как пластинчатые, скребковые, винтовые, ковшевые конвейеры, а при транспортировании горячих сыпучих материалов - и
ленточные. Ограничиваться это может только налипанием пылевидных материалов на рабочие органы машин при угле подъёма рабочего органа более 5°.
Вибротранспортирование отличается простотой конструкций машин и сложностью процессов, протекающих в перемещаемом материале. При выборе кинематических параметров виброконвейеров и вибропитателей учитываются условия их эксплуатации, необходимая производительность и свойства груза (гранулометрический состав, объёмный вес, высота слоя).
В общем случае, при вибротранспортировании реализуется либо режим движения с подбрасыванием, когда груз периодически теряет контакт с грузоне-сущим органом, либо безотрывный режим движения в виде скольжения по рабочему органу.
На характер процесса значительное влияние оказывает соотношение сил трения и сил инерции в слое транспортируемого материала. При движении виброповерхности из нейтрального положения вверх и вперёд силы трения уменьшаются, а силы инерции увеличиваются, и, как только силы инерции превысят силы трения, частица материала начинает перемещаться по виброповерхности. Когда амплитудное значение вертикальной составляющей ускорения виброповерхности превысит ускорение силы тяжести, частица материала от неё отрывается и совершает полёт по параболе.
Виброповерхность с находящейся на ней частицей с массой m колеблется под углом вибрации /3, с амплитудой А и частотой ш по закону S = A sin t.
На частицу действует сила тяжести mg, сила сухого трения ^тр, нормальная реакция N и сила инерции FEE = —m№ sint (рис. 1).
Уравнения движения частицы относительно виброповерхности будут иметь вид:
mx = тА cos (a + [3) sin wt + FTp — mg sin a
ту = mA2 sin [3 sin ut — mg cos a + N.
(1)
При движении частицы по виброповерхности координата y равна нулю, а сила трения Fp определяется соотношениями
—fN при х > 0 ^Тр = I +fN при х < 0 (2)
fiN при х = 0,
где f - коэффициент трения скольжения; fi - коэффициент трения покоя.
При относительном покое частицы на вибрирующей плоскости (совместное движение при х = 0 и у = 0) нормальная реакция N определяется из второго уравнения (1):
N = mg cos a — тA2 sin ¡3 sin ut.
Для обеспечения безотрывного режима работы необходимо, чтобы частица материала оставалась на вибрирующей поверхности, т. е., чтобы соблюдалось условие N > 0, или
9cosa ^ ■ , (3)
> sin шъ. (3)
№ sin /3
Рис. 1. Схема действия сил на частицу насыпного груза: а - угол наклона грузонесущей поверхности; р - угол вибрации; тд - сила тяжести; N - нормальная реакция грузонесущей поверхности; ^ин - сила инерции; ^ин.нач - сила инерции в начальный момент времени; ХОУ - система отсчёта, связанная с грузонесущей поверхностью; X'О'У' - абсолютная
система отсчёта
Частица остаётся в состоянии относительного покоя на плоскости при условии, что сила трения F-ф не превосходит по абсолютной величине предельного значения силы статического трения f1N. Это условие может быть представлено в виде неравенства
—fiN (t) <Fo (t) < f\N(t). (4)
Подставляя выражения для силы трения (2) в уравнение (1) для тх, получим уравнение скольжения частицы по колеблющейся плоскости
х = Аш2 cos (а + Р) (1 ±-¿ЛР- | sin uí ^ g (j cos a — sjn a). (5)
\ cos а — sin а tgp )
Верхние знаки в этом уравнении соответствуют «скольжению частицы вперёд» (обгон плоскости, х > 0), нижние - «скольжению частицы назад» (запаздывание от плоскости, х < 0).
Для горизонтальной (а = 0°) плоскости уравнение (5) примет вид
х = Аш2 cos Р (1 ± Kf) sin ut ^ f д, (6)
где Kf = ftg.
Из соотношений (3) и (4) следует, что нормальная реакция горизонтальной поверхности в любой момент времени положительна при выполнении условия
к = ^ i
9
и частица, попавшая на виброповерхность с нулевой поперечной составляющей скорости (у = 0), в дальнейшем остаётся на ней.
Безразмерный параметр Кр, представляющий собой отношение амплитудного значения вертикальной составляющей ускорения колеблющейся плоскости к ускорению силы тяжести, называется коэффициентом режима работы вибро-транспортирующей установки.
Если это неравенство не выполняется, то нормальная реакция N(t) становится равной нулю и происходит отрыв частицы от грузонесущей поверхности (начало полёта). Момент отрыва t0 определяется уравнением
9 1
sin^0 = л 2 ■ п = ТГ".
Aw2 sin р Кр
Дифференциальные уравнения полёта частицы над плоскостью в подвижной системе координат XOY получаются из уравнения (1) при условии, что N = 0 и FVp = 0.
{
х = Aw2 cos ¡3 sin wt у = Аш2 sin ¡3 sin wt — g.
Так как движение частицы по колеблющейся плоскости без отрыва от неё может происходить лишь при Kp<1, возможны следующие режимы её скольжения:
1. «Скольжение вперёд» и «скольжение назад» с двумя длительными остановками в каждом периоде.
2. «Скольжение вперёд» и «скольжение назад» с двумя мгновенными остановками в каждом периоде.
3. «Скольжение вперёд» и «скольжение назад» с одной длительной и с одной мгновенной остановками в каждом периоде.
4. Скольжение в одном направлении с одной длительной остановкой в каждом периоде.
При движении частицы по плоскости без подбрасывания скорость перемещения её и координата находятся интегрированием уравнения (6) при начальных условиях: х (t*) = х*; х (t*) = 0.
х (t) = — Аш cos (1 ± Kf) (cos wt± — cos wt*) T f g (t± — t*) + x*,
x (t) = Aw cos ¡3 (1 ± Kf) cos wt * (t± — t*) — A cos ¡3 (1 ± Kf) (sin w t± — sin w t *) ^
T f 9*(t ± — t*). Принимая x* = 0 (начало скольжения), имеем
х (t) = —Aw cos (1 ± Kf) (cos wt± — cos wt*) T fg (t± — t*), (7)
ж (t) = Aw cos [3 (1 ± Kf) cos шt* (t± — t*) —
—A cos /3 (1 ± Kf) (sinwt± — sin wt*) T f g~ ±
(t± — t*)2
Момент смены режима движения определяются приравниванием уравнения (7) к нулю:
К/
cos ut — cos ut± ^
K (1 ±Kf)
(ut± — ut *) = 0.
Обозначив через фазовый угол, соответствующий моменту начала скольжения частицы вперёд или назад, а через <р± фазовый угол момента остановки после скольжения частицы и вводя параметр
= ±
K
f
KP(1 ±Kf)'
получим трансцендентное уравнение, определяющее момент смены этапа дви-
жения:
cos p± — cos + z± (p± — í±) = 0.
(8)
Решая уравнение (8) графоаналитическим методом, определяется фазовый угол p± как абсцисса, ближайшая к точке пересечения прямой k = ^z±, проведённой из точки (5±, cos5±) и косинусоиды z = cosí (рис. 2).
Рис. 2. Графическое определение фазовых углов смены этапа движения (Р = 8°; / = 0,6)
Применяя аппарат линейного программирования (рис. 3), можно получить аналитическое решение уравнения (8), и определить граничные параметры режима скольжения. Для 1-го режима скольжения моменты перехода от одного этапа к другому можно обозначить как ¿0+, t+0, ¿о- и ¿_0.
Углы начала скольжения вперёд и скольжения назад определяются из выражений:
= агсвт
{
ж — arcsin 2— при Р + р < 2 arcsin 2— при Р + р > 2.
Фазовые углы смены режимов движения у+0 и у-0 определяются графоаналитическим способом.
Перемещение частицы на этапе скольжения вперёд (у+0 —
= А cos р (1 + Kf)
+
( p+o —
+ ( p+o — 5+) cos 5+ — (sin p+o — sin &+)
Рис. 3. Блок-схема расчёта фазовых углов
Перемещение частицы на этапе скольжения назад (<р-0 — 8-) будет равно
5- = А 008,5(1 — Kf)
-zJ ^ — SJ) +
+ ( p-o — 8-) cos 8- — (sin p-o — sin (L) Средняя скорость движения за время одного периода
S+ + 5-
уср = "^Т"w
или
i
Уср = Аш cos(3 —\ (1 + K/)
— (sinp+o — sin #+)
Z+-
( P+o — 8+У
+ ( p+o — í+)cos 8+—
+ (1 — Kf)
—z +
+ ( p-o — -) cos - (sin p-o — sin -)
(9)
Вводя безразмерный коэффициент К^ (I = 1,2,3,4) передачи горизонтальной составляющей скорости, равный для 1-го режима работы
cos В { .
K*i = !Л (1+K'
2ж
+ (1 — Kf)
( p+o — 8+У
+ ( p+o — ¿+)cos 8+—(sin p+o — sin ¿+)
+
( p-o — -—)2
+ ( p-o — -) cos - — (sin p-o — sin -)
(10)
получим выражение для средней скорости частицы при г-м режиме скольжения
Кр. = АшК^,
где г - режим скольжения частицы.
Для 2-го режима, переход от одного этапа движения к другому будет происходить в моменты времени ¿+-и ¿-+. Учитывая, что при 2-м режиме скольжения = 5-+ и (р-+ = 5+- фазовые углы смены этапов будут равны :
{
cos у+_ = cos у_+ — (у+_ — У_+) 2+ cos у_+ = cos у+_ — (у_+ — у+_ + 2ж) z_.
Решение этих уравнений для определения фазовых углов моментов перехода, будет иметь вид:
у+_ = ж + 2 (1 + Kf) — arcsin — у_+ = 2 (1 — Kf) — arcsin
!
где
ж Kf
- =п
1
2 Kp cos (2Kf
Перемещение частицы за этап скольжения вперёд ( кр+- — и назад
+ 2ж — определится как
S+ = A cos [(1 + Kf)
(У+_ — у_+)
S_ = A cos [ (1 — Kf)
— ( (fi+_ — у_+) cos у_+ — (sin у+_ — sin у_+)
(у_+ +2ж — v+_y
— (у_+ + 2ж — <p+_) cos у+_ — (sin у_+ — sin у+_)
(11)
Тогда коэффициент передачи горизонтальной (полезной) составляющей средней скорости Кг2 в выражении (10) при 2-м режиме скольжения будет равен
cos [
K (1+K/)
— (sin у+_ — sin у_+)
( У + _ — У_+)2 / N
z+----( у+_ — ^_+)cos^_+—
+ (1 — Kf)
(у_+ +2ж — у+_)2
— (у_+ + 2ж — <р+_) cosy+_ — (siny_+ — siny+_)
(12)
2
2
2
При 3-м режиме скольжения, моментами перехода от одного этапа к другому будут ¿о+, t+- и ¿-0. Угол начала скольжения вперёд после длительной остановки ¿о+ = шЬ0+ и фазовые углы, соответствующие моментам перехода от скольжения вперёд к скольжению назад у+- = и от скольжения назад к
длительной остановке у-0 = шЬ-0, определяется также, как и в 1-м режиме.
Выражение для суммарного перемещения частицы за этапы скольжения вперёд ( — £0+) и назад ( у-0 — у+-) будет иметь вид
5 = 5+ + 5- = А 008 , (1 + К1)
о Г
+ ( — ¿0+) 008 ¿0+ — (вт — 8Ш ^0+)
+А сов ,(1—К/)
( <Р-0 — У+-)2
+ ( (р-0 — ) 008 — (вт У-0 — 8Ш (р+-)
(13)
Тогда коэффициент передачи горизонтальной составляющей скорости в уравнении (10) запишется как
К
008 , 2тт
{
(1 + К)
+
( — 50+У
— ( — ^0+) 008 80+ —
— (81п — 8Ш #0+)
+ (1 — К/ )я
(<Р-0 — у+-)'
— (81п у-0 — Б1тр+-)
— (у-0 — 008 —
(14)
Моментами перехода в 4-м режиме будут ¿0+ и ¿+0. Угол начала скольжения вперёд после длительной остановки и фазовый угол, соответствующий моменту перехода от скольжения вперёд к длительной остановке у+0, определяется так же как в 1-м и 3-м режиме.
Перемещение частицы на этапе скольжения вперёд (у+0 — #0+):
5+ = А 008, (1 + К)
+
(у+0 — ¿0+)2
+ (у+0 — 80+) 008 ^0+ — (81п у+0 — 8Ш ^0+)
(15)
Для коэффициента передачи горизонтальной составляющей скорости в уравнении (10) будем иметь:
008 , . К, 4 = (1 + К)
+
( у+0 — ^0+)2
— ( у+0 — 8+0) 008 ¿0+ — (81п у+0 — 8Ш ^0+)
(16)
2
В конкретных случаях при решении практических задач для определения степени влияния коэффициента трения и угла вибрации на эффективность вибротранспортирования целесообразно вычислить все возможные значения Kz и построить графики изменения Kz в зависимости от коэффициента трения при заданных значениях коэффициента режима работы и угла вибрации [.
Для построения графиков зависимости коэффициента передачи горизонтальной составляющей Kz при ограничениях, задаваемых коэффициентом режима работы КР, при определённых значениях угла вибрации и при изменении коэффициента трения f нужно задаться последовательно величиной Кр = [0,2; 1,0] c шагом 0,1 и, при постоянном значении угла вибрации [, выяснить режим движения. Затем по формулам (9), (11), (13) и (15) вычисляются значения Kz при различных значениях коэффициента трения.
Пользуясь этими графиками, легко можно оценить эффективность вибротранспортирования и вычислить скорость виброперемещения при заданных значениях A, ш, [, Кр и известном коэффициенте трения скольжения транспортируемого груза по формуле (10), где значение Kz берётся из соответствующего графика.
При построении графиков зависимости К = F(/) при постоянном значении коэффициента режима работы вычисления достаточно произвести для одного угла вибрации. Для этого находятся условия, при которых коэффициент передачи скорости, определяемый равенствами (10), (12), (14) и (16), будет постоянным при неизменном значении коэффициента режима работы и переменном угле вибраций.
Анализ выражений (10), (12), (14) и (16) показывает, что коэффициент передачи скорости будет постоянным при неизменной величине коэффициента режима работы и изменяющемся угле вибраций лишь при условии
ft д [ = con s t. (17)
Таким образом, зная значение коэффициента передачи скорости для одного угла вибрации [, легко найти по равенству (17) значение коэффициента трения f при другом угле вибрации, которому соответствует такое же значение Kz при неизменном значении коэффициента режима работы.
Анализ полученных данных для значений коэффициента передачи горизонтальной составляющей скорости вибраций показывает, что в режимах движения «без подбрасывания» коэффициент трения материала о грузонесущую поверхность рабочего органа и угол вибраций оказывают существенное влияние на эффективность вибротранспортирования. Коэффициент передачи горизонтальной составляющей скорости вибраций при постоянном значении коэффициента режима работы сначала возрастает с увеличением коэффициента трения, достигая максимума при определённом соотношении коэффициента трения и угла вибраций, а затем плавно снижается.
Это обстоятельство необходимо учитывать при выборе параметров вибро-транспортирующей установки и подбирать такие соотношения между коэффициентом режима работы и параметром [, которые обеспечат максимальное значение коэффициента передачи скорости.
Литература
1. Быховский И.И. Основы теории вибрационной техники. М. : Машиностроение, 1968. 362 с.
2. Бауман В.А., Быховский И.И. Вибрационные машины и процессы в строительстве. М. : Высшая школа, 1977. 255 с.
3. Гончаревич И.Ф., Гудушаури Э.Г., Гаврилина Л.В. Модернизация транспортно-перегрузочной техники методами вибротехнологий // Машиностроение и инженерное образование. 2014. № 2. С. 9-14.
4. Тропман А.Г., Бельков Н.И., Макеева Ю.Н. Вибрационные конвейеры для транспортирования горячих материалов. М. : Машиностроение, 1972. 327 с.
5. Гончаревич И.Ф. Динамика вибрационного транспортирования. М. : Наука, 1972. 244 с.
6. Бельков Н.И. Исследование закономерностей и особенностей вибрационного перемещения тонкоизмельченных сыпучих материалов: дис. канд. техн. наук. Л., 1974. 178 с.
7. Локтионова О.Г. Динамика вибрационных технологических процессов и машин для переработки неоднородных гранулированных сред: дисс. д-ра техн. наук. Курск, 2008. 238 с.
8. Леонтьева А.И., Брянкин К.В., Дегтярев А.А., Орехов В.С. Моделирование и расчёт процесса сушки термолабильных материалов в виброаэрокипящем слое: монография. Тамбов : Изд-во ГОУ ВПО ТГТУ, 2011. 80 с.
9. Бабичев А.П. Технологическое применение колебаний или . . . вибрационные технологии // Вестник ДГТУ. 2005. Т. 5, № 3. С. 1-11.
10. Гончаревич И.Ф., Фролов К.В. Теория вибрационной техники и технологии. М. : 1981. 320 с.
11. Членов В.А., Михайлов Н.В. Виброкипящий слой. М. : Наука, 1972. 340 с.
12. Архипенко В.П. Механико-технологические основы создания вибрационных машин для бункеризации, выпуска и контейнерной обработки: дисс. д-ра техн. наук. Днепропетровск, 1996. 346 с.
13. Гончаревич И.Ф. Вибрация - нестандартный путь: вибрация в природе и технике. М. : Наука, 1986. 209 с.
14. Силин А.Н. Обоснование параметров и разработка вибрационного выпускающщего устройства с маятниковой подвеской для сыпучих материалов: дисс. канд. техн. наук. М., 1993. 182 с.
15. Архипенко А.В. Активные бункерные вибростенки: дис. канд. техн. наук. Белгород, 1999. 143 с.
16. Лютов В.Н. Обоснование параметров и разработка вибропитателей с криволинейным профилем рабочего органа для выпуска сыпучих материалов из бункеров: дисс. канд. техн. наук. М., 1990. 170 с.
CONSERVATION OF THE SOIL PARTICLES BY THE WORKING ORGANS
OF THE VIBROMACHINE
A.V. Arhipenko
Ph.D. (Eng.), e-mail: iten-miu@mail.ru
International Innovation University, Sochi, Russia
Abstract. It is shown that in the modes of motion "without tossing" the coefficient of friction of the material on the load-carrying surface of the working body and the angle of vibrations have a significant impact on the efficiency of vibrotransportation. The transfer coefficient of the horizontal component of the vibration velocity at a constant the value of the operating mode coefficient first increases with an increase in the friction coefficient, reaching a maximum at a certain ratio of the coefficient of friction and the angle of vibration, and then gradually decreases.
Keywords: sliding particles, soil, vibrators, vibrotransportation.
Дата поступления в редакцию: 25.01.2023