Дозирование
 Скачать 5.11 Mb.
|
|
Рис. 5.22. Перемещение центра тяжести порции сыпучего материала вдоль наклонного вибрирующего лотка при углах наклона его к горизонту: ♦ 1 – α = 9°; ■ 2 – α = 11°; ▲ 3 – α = 13° Для устранения недостатков данной физической модели автор работы [31] воспользовался подходом к описанию процесса движения тел под действием вибрации, который основан на рассмотрении процесса движения как дискретного, т.е. движение тела складывается из множества независимых перемещений, причем время одного перемещения равно периоду колебания. Для детального описания данного процесса был проведен ряд исследований, в которых на первом этапе вместо сыпучего материала использовалось твердое тело [32]. На рисунке 5.23 показана схема действующих сил на твердое тело, расположенное на наклонном лотке. При отсутствии вибрации и угле наклона лотка меньше, чем угол трения покоя, тело неподвижно. Вибрационное воздействие, перпендикулярное к поверхности лотка, позволяет на определенном полупериоде колебания (лоток движется вниз) уменьшить силу нормального давления, а следовательно, и силу трения до такого значения, что тело начнет свое движение под действием скатывающей силы mg sin α, причем будет двигаться с переменным положительным ускорением. Затем (лоток движется вверх) тело будет двигаться с переменным отрицательным ускорением вплоть до остановки. Графически силовое воздействие силы инерции можно представить следующим образом (рис. 5.23). Рис. 5.23. Схема сил, действующих на твердое тело, находящееся на наклонном вибрирующем лотке mg sin α f тр (mg cos α – P ц ) α mg cos α mg Рис. 5.24. Изменение ускорения движения тела в направлении, перпендикулярном к поверхности лотка Если положительного ускорения K p y K a p ′′ > достаточно для того, чтобы преодолеть силу трения и обеспечить движение тела вниз по наклонному лотку, то на протяжении ∆ϕ д тело будет двигаться с переменным положительным ускорением (от точки А до точки В). На участке от точки В до точки. С тело будет двигаться с отрицательным ускорением и в точке С остановится. На участке от точки С до точки А тело будет неподвижно, а затем цикл движения повторится. Ускорение движения тела от точки А до точки С можно определить, используя следующее соотношение )) ( cos ( sin тр t y m mg f mg ma ′′ + α − α = (5.37) Поскольку для начала движения необходимо преодолеть силу трения, ускорение на участке разгона определяется следующим образом: 0 0 2 тр ) ( sin cos sin A K t А g f g a p − ω ω + α − α = . (5.38) Скорость тела в конце участка разгона будет равна 1 0 0 1 тр 1 1 0 0 | ) cos( ) cos sin ( t A K t А t g f g V p t t t ∆ − ω ω − ∆ α − α = ∆ + . (5.39) Скорость тела в конце участка торможения можно определить, используя зависимость 2 0 0 2 тр 1 2 2 1 0 1 0 | ) cos( ) cos sin ( t A K t А t g f g V V p t t t t t ∆ − ω ω − ∆ α − α + = ∆ + ∆ + ∆ + . (5.40) Если V 2 = 0, то 2 0 0 2 тр 1 2 1 0 1 0 | ) cos( ) cos sin ( t A K t А t g f g V p t t t t t ∆ + ω ω + ∆ α − α − = ∆ + ∆ + ∆ + . (5.41) Приравняем (5.39) и (5.41), получим ) ( cos ) ( cos ) cos sin ( cos ) ( cos ) cos sin ( 2 0 1 0 0 2 1 0 0 2 тр 1 0 0 0 1 0 0 1 тр t A K t t А t t t А t g f g t A K t А t t А t g f g p p ∆ − ∆ + ω ω − − ∆ + ∆ + ω ω + ∆ α − α − = = ∆ − ω ω + ∆ + ω ω − ∆ α − α (5.42) Анализ показывает, что данное равенство справедливо при выполнении двух условий: 0 cos sin тр = α − α g f g и 0 0 = A K p . Учитывая постановку задачи, данные условия невыполнимы, следовательно, ∆t 2 всегда меньше ∆t 1 . Следует отметить, что участок торможения всегда меньше чем участок разгона. Численное решение уравнения (5.42) в диапазоне изменения углов наклона лотка к горизонту и амплитуд колебаний по длине лотка показало, что участок торможения составляет от 0,54 до 0,43 от времени разгона. Учитывая незначительное влияние на общее время движения, будем оценивать время движения по времени разгона. Поскольку амплитуда колебаний лотка уменьшается по его длине, то также будет уменьшаться время движения тела за одно колебание пластины. Для определения времени движения за одно колебание рассмотрим участок разгона, который представлен на рис. 5.25. Функция, описывающая ускорение колебаний лотка в точке загрузки порции, может быть записана в виде 1 0 1 sin α = L A y , а на расстоянии S как 2 0 2 sin ) ( α − = S k A y L L . Тогда для определения ∆t(S) необходимо решить следующую пропорцию: Рис. 5.25. Схема к расчету времени движения за одно колебание лотка 2 1 0 2 2 ) ( α − π α − π = ∆ ∆ t S t , где ( ) р a L р K S k a y K arcsin ; arcsin 2 0 0 1 = α − ′′ = α Окончательно получим: 0 0 0 arcsin 2 arcsin 2 ) ( t K S k a y K S t р a р ∆ − π − ′′ − π = ∆ , (5.43) где K р , k a – коэффициенты пропорциональности; 0 y ′′ – ускорение, перпендикулярное к поверхности лотка в точке загрузки порции сыпучего материала, м/с 2 ; S – путь, пройденный центром тяжести твердого тела, м; ∆t 0 – дискретная величина промежутка времени, с. Для описания процесса движения твердого тела по наклонному лотку в дискретной форме, исходя из теоретического анализа процесса, были сделаны [31] следующие допущения: 1. За период одного колебания лотка тело движется с некоторой постоянной скоростью, т.е. тело за отрезок времени ∆t пройдет путь t V S ∆ > =< ∆ , (5.44) где <V > – средняя скорость, м/с; ∆t – время движения за интервал времени ∆t, с. 2. Средняя скорость, пропорциональная пройденному пути, т.е. расстоянию от загрузочного края до центра тяжести тела, на последующем периоде колебания уменьшается: S K V V v − > >=< < 0 , (5.45) где <V 0 > – средняя начальная скорость, мм/с; K v – коэффициент пропорциональности. 3. Время движения ∆t уменьшается по длине лотка в соответствии с (5.43). Исходя из данных допущений, были разработаны алгоритм и программа расчета на ПЭВМ. Для проверки адекватности предложенной модели провели серию экспериментов. Методика эксперимента заключалась в следующем. В наклонный лоток (рис. 5.26) помещалось твердое тело (стальной параллелепипед размером 20 × 15 × 10 мм), как показано на рис. 5.26, а. Включали вибратор и проводили видеосъемку процесса движения (видеокамера Panasonic NV-DS65). Полученное изображение раскладывалось на отдельные кадры. По полученным видеоматериалам определяли положение центра тяжести A 0L A L (S ) K p y 0 ″ ∆t(S) y ″ t ∆t 0 твердого тела в каждый конкретный момент времени. Затем процесс повторяли снова, чтобы получить среднее значение координаты центра тяжести твердого тела при фиксированных параметрах процесса (угол наклона лотка, амплитуда вибрации загрузочного края лотка, частоты колебаний лотка). По полученным средним экспериментальным значениям проводили идентификацию параметров (V 0 , K p ). На рисунке 5.27 сплошной линией показана функция S(t), полученная по модели, а точками – экспериментальные данные. а) 0 с д) 8 с б) 2 с и) 10 с в) 4 с к) 12 с г) 5 с л) 14 с Рис. 5.26. Разложение по кадрам видеоизображения процесса движения твердого тела по наклонному вибрирующему лотку Рис. 5.27. Экспериментальные (точки) и расчетные (линия) данные движения центра тяжести твердого тела Анализ результатов расчета по данному алгоритму и их сравнение с экспериментальными данными показали, что значение скорости 0 > прямо пропорционально синусу угла наклона лотка к горизонту, а значение K p прямо пропорционально косинусу данного угла. Эти закономерности не противоречат физической сущности рассматриваемого процесса. Действительно, ускорение движения тела, а следовательно, и скорость прямо пропорциональны синусу угла наклона лотка к горизонту, сила трения прямо пропорциональна косинусу данного угла. При переходе от модели движения центра тяжести твердого тела к модели движения центра тяжести порции сыпучего материала вдоль наклонного вибрирующего лотка были приняты следующие допущения [31]: – на протяжении каждого цикла движения центр тяжести ускоряется, затем затормаживается и в конечном итоге скорость его перемещения становится равной нулю; – средняя скорость движения, в пределах одного цикла, изменяется линейно по длине лотка, уменьшаясь от загрузочного к разгрузочному краю; – время движения центра тяжести с указанной средней скоростью, в пределах одного цикла, уменьшается при перемещении его от загрузочного к разгрузочному краю согласно уравнению (5.43); – после начала ссыпания материала с лотка центр тяжести порции сыпучего материала движется с постоянной скоростью, равной скорости его движения на момент начала ссыпания. Как и в модели движения отдельного твердого тела, в математической модели движения отдельной порции сыпучего материала по наклонному вибрирующему лотку идентификации подлежат два параметра V 0 и K p Экспериментально было подтверждено [33], что средняя скорость движения центра тяжести порции сыпучего материала прямо пропорциональна синусу угла наклона лотка к горизонту, а значение K P прямо пропорционально косинусу того же угла. 5.3.4. МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ПРОЦЕССА ПРЕОБРАЗОВАНИЯ ОТДЕЛЬНЫХ ПОРЦИЙ В НЕПРЕРЫВНЫЙ ПОТОК Процесс преобразования отдельных порций в непрерывный поток, помимо изменения формы и движения отдельной порции сыпучего материала, включает в себя объединение двух соседних порций. Данный процесс может быть описан с использованием математического аппарата случайных марковских процессов, дискретных в пространстве и времени. Детальное рассмотрение данного процесса выполнено в работе [31]. Предварительное исследование процесса объединения порций проводилось с использованием прозрачного короба, показанного на рис. 5.28. а) а) б) б) в) в) Рис. 5.28. Соединение порций (одинаковый объем порций) Рис. 5.29. Соединение порций (разный объем порций) В короб засыпались две порции равного объема одного и того же сыпучего материала, отличающиеся по цвету, как показано на рис. 5.28, а. Затем включали вибратор и проводили видеосъемку, результаты которой представлены на рис. 5.28, б, в. Как видно из снимков, при соединении порций образуется строгая граница. Иная картина складывается при загрузке порций с разными объемами, как это показано на рис. 5.29, а, б, в. Как видно из рисунка, материал из порции с большим объемом накладывается на материал из порции меньшего объема. Аналогичная картина складывается при соединении двух порций в районе разгрузочного края лотка (рис. 5.30). По результатам экспериментов было сделано допущение о том, что в процессе преобразования порций сыпучего материала на наклонном вибрирующем лотке поведение отдельных порций можно считать независимым, а результирующее количество сыпучего материала на интервале объединения порций рассматривать как сумму двух независимых функций. 0 с 6 с 2 с 8 с 4 с 10 с Рис. 5.30. Соединение порций на наклонном лотке С учетом сделанного допущения процесс преобразования отдельных порций в непрерывный поток был представлен в следующем виде. Через равные промежутки времени ∆Т в загрузочную часть лотка подается i-я порция материала. Для конкретного сыпучего материала форма и объем порции (при единичной ширине лотка – площадь продольного сечения) однозначно определены амплитудой A 0 (5.31) и полупериодом T 0 / 2 (5.30). С течением времени t центр тяжести i-й порции сыпучего материала перемещается вдоль лотка. Поскольку процесс рассматривается как дискретный (время t изменяется дискретно с шагом ∆t), то на каждом переходе, т.е. для каждого нового значения времени t, по зависимостям (5.32) и (5.33) рассчитывается амплитуда A(t) и полупериод T(t) / 2. Координаты центра тяжести и амплитуды предыдущей порции i – 1 и последующей i + 1 рассчитываются по тому же алгоритму при значениях времени t – ∆T и t + ∆T, соответственно. При начале ссыпания материала с лотка, т.е. при выполнении условия S i (t)+ T(t) / 4 ≥ L L , скорость движения центра тяжести считаем постоянной. Можно для дискретных моментов времени t = i ∆t по зависимостям (5.32), (5.33), (5.44) определить численные значения A(t), Т(t) и S(t). Условие начала ссыпания материала с лотка имеет следующий вид: 0 4 ) ( ) ( S L k T k S L − ≥ + , (5.46) где L L – длина лотка; S 0 – координата центра тяжести порции сыпучего материала после загрузки ее на лоток. Если для порции с номером z ссыпание начнется по истечению времени TF с момента ее загрузки на лоток, то предыдущая порция с номером z – 1 будет находиться в движении в течение времени TF + ∆Т. На рисунке 5.31 схематично показана ситуация к моменту начала ссыпания материала порции z с лотка, при этом часть материала порции z – 1 к этому моменту времени будет еще оставаться на лотке. Рис. 5.31. Схема к расчету производительности вибрационного дозатора в момент времени t Количественно эту часть материала можно определить, используя следующее выражение: [ ] π ϕ − = ∆ ∆ + − ∆ + ∆ + − T TF z T TF T TF z Т A A F 0 , 1 0 , 1 cos , (5.47) где A TF+ ∆T – амплитуда, рассчитанная по зависимости (5.32) в момент времени t = TF + ∆T; F 0 – площадь продольного сечения порции после ее загрузки на лоток. Численное значение ϕ z–1, 0 определяется следующим образом: [ ] TF F T TF T TF z T TV T T 25 , 0 25 , 0 5 , 0 0 , 1 + ∆ − π = ϕ ∆ + ∆ + − , (5.48) где Т TF+ ∆T – период, рассчитанный по зависимости (5.33) в момент времени t = TF + ∆T; V F – скорость центра тяжести порции сыпучего материала на момент начала ссыпания; Т TF – период, рассчитанный по зависимости (5.33) в момент времени t = TF. Состояние системы по истечению интервала времени ∆t схематично показано на рис. 5.32. Центры тяжести порции z и z – 1 переместятся на расстояние S = V F ∆t, соответственно изменятся амплитуда и полупериод. Часть материала порции z сыпется с лотка, а от порции z – 1 на лотке еще останется часть материала, но меньше чем в предыдущем случае. Количество материала, ссыпавшегося с лотка из порции z, можно определить, используя следующее выражение: [ ] π ϕ − π − − = ∆ ∆ + ∆ + ∆ + T TF z T TF T TF z Т A A F F ) ( cos 1 , 0 1 , , (5.49) где [ ] TF t TF F T TF z T T tV T 25 , 0 25 , 0 5 , 0 1 , − + ∆ π = ϕ ∆ + ∆ + . (5.50) |