Файл: Лукьянов, П. И. Аппараты с движущимся зернистым слоем. Теория и расчет.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 24.10.2024
Просмотров: 72
Скачиваний: 0
|
|
X |
|
|
|
3r)g m |
|
(2 - | ) ( / - ätö* |
|
|
|
|||
|
|
«ft,i [iW + (t - |
kr))2}b/*8y |
|
|
|
||||||||
|
|
|
|
|
|
|
||||||||
|
|
, |
J_ /2.,2_c |
(/ — fer|)' |
|
2T]2 (1 — kr\)2 |
|
|
||||||
|
|
|
/ |
^ |
ta,n2 + |
(^ —кц)' |
+ (* +WГ)2 + (/ — ÄT])2_ |
|
|
|||||
где |
t, |
k ■— |
индексы, |
определяющие соответственно |
столбец |
и |
||||||||
строку |
узловой |
точки; |
I = hlD\ |
т] — шаг интегрирования; |
і |
= |
||||||||
= |
0, |
1,2, |
. . ., |
D/2r|; |
k — |
|
|
|
|
|
||||
= |
0, |
1, |
2, |
. . ., |
kmm — D — |
|
1 |
|
|
|
||||
— d j 2 tg a 3T]. |
|
|
|
|
|
|
|
|
|
|||||
= |
Граничные условия: & = |
|
|
|
|
|
||||||||
|
v0<i = |
0; |
u0,1 — tiQ. |
|
|
|
|
|
||||||
|
Система |
|
уравнений |
(94) |
|
|
|
|
|
|||||
положена в основу оператор |
|
|
|
|
|
|||||||||
ного |
алгоритма |
нахождения |
|
|
|
|
|
|||||||
поля скоростей |
при |
движе |
|
|
|
|
|
|||||||
нии |
сыпучего |
материала в |
|
|
|
|
|
|||||||
емкости. |
|
|
|
|
реали |
|
|
|
|
|
||||
|
Идея |
алгоритма, |
|
|
|
|
|
|||||||
зующего |
|
вычислительную |
|
|
|
|
|
|||||||
процедуру (рис. 91), состоит |
|
|
|
|
|
|||||||||
в следующем: |
|
|
|
|
|
|
|
|
|
|||||
|
по заданным физико-меха |
|
|
|
|
|
||||||||
ническим свойствам сыпучей |
|
|
|
|
|
|||||||||
среды и геометрическим раз |
|
|
|
|
|
|||||||||
мерам аппарата определяются |
|
|
|
|
|
|||||||||
основные размерызоны дви |
|
|
|
|
|
|||||||||
жения; |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
производится подстановка |
|
|
|
|
|
||||||||
граничных |
|
значений |
|
ско |
|
|
|
|
|
|||||
рости; |
|
|
|
|
горизон |
|
|
|
|
|
||||
|
рассчитываются |
|
|
|
|
|
||||||||
тальная |
и |
вертикальная со |
|
|
|
|
|
|||||||
ставляющие |
скорости |
для |
Рис. 91. Блок-схема операторного |
алго |
||||||||||
k-то шага интегрирующей |
|
ритма |
|
|
|
|||||||||
сетки; |
|
|
|
условие |
выдачи |
результатов счета |
vkt { и |
ukt t |
||||||
|
проверяется |
на печать и анализируется условие окончания расчета (сравне ние k с &шах).
Вслучае невыполнения условия k = ßmax производится уве личение k на единицу (переход к следующей строке интегрирую щей сетки), и расчет повторяется. При выполнении этого условия вычисления прекращаются.
Всоответствии с изложенным алгоритм решения уравнений
(94)записан в виде
t^i2>
171
где |
А у |
— оператор вычисления а 3; |
R 2— проверка условия |
|
a t |
а 3; |
a t — текущий |
угол, изменяющийся пропорционально |
|
шагу |
интегрирования; |
Л3, Л4 — блоки |
подстановки граничных |
условий для горизонтальной и вертикальной составляющих ско рости:
z*. i = h — kt}D, xkt i = (h — kt]D) tg a (.; v0,i = 0, u0ti — u0;
Kb— блок присвоения; просчет значений вертикальной и гори зонтальной составляющих скорости для й-го номера сетки интегри
рования; |
Л |
Л7 — операторы |
|
|
|
|
|
|||||
вычисления |
горизонтальной |
и |
О |
0,1 |
0,2 . 0,3 O ßW x,M » |
|||||||
вертикальной составляющих |
ско |
|
|
|
|
|
||||||
рости; Кв — счетчик шагов |
|
|
|
|
|
|
||||||
|
|
|
|
|
ум |
|
|
|
|
|
|
|
|
|
|
|
|
0.9-200 |
|
|
|
|
|
||
|
|
|
|
|
0,8 |
|
|
|
|
|
|
|
|
|
|
|
|
0,7 |
|
|
|
|
|
|
|
|
|
|
|
|
0,6 |
|
|
|
|
|
|
|
|
|
|
|
|
0,5 |
|
|
|
|
|
|
|
Y«A |
г |
|
і |
о |
0,1-200 |
Рис. 93. Графики распределения |
||||||
|
|
|
||||||||||
и>,і |
' |
|
|
|
|
|
вертикальной составляющей скорос |
|||||
Рис. |
92. Графики |
распределения |
го |
ти в горизонтальном сечении: |
||||||||
1 ~ z = |
D; J2 — z = 0,9D; 3 —z = |
0,8D', |
||||||||||
ризонтальной |
составляющей |
скорости |
||||||||||
4 ~ г = |
0,7D; |
5 — z = |
0,60; 6 |
— г = |
||||||||
|
по высоте |
аппарата |
|
|
= 0 ,5 0 ; 7 — z = |
0,40 |
|
|||||
сетки интегрирования; |
Р9— проверка условия вывода на печать |
|||||||||||
результатов |
счета |
vkti |
и uk<i\ Я іа — выдача |
результатов счета |
на печать; Р 1г-— проверка условия окончания счета k < М (где
М = kmiX + 1); |
Я а — конец вычислений. |
Результаты расчета по приведенному алгоритму представлены |
|
на рис. 92 и 93. |
Точность вычислений vk, t и uk, t во многом опре |
деляется величиной шага интегрирования г]. Последний выбран из условия экономии машинного времени и удовлетворительной сходимости результатов расчета со значениями скорости, вы численными по формуле (93). С достаточной для инженерных расчетов точностью принята величина максимальной ошибки 5%. Вычисления показали, что такая точность достигается при r\D 0,5 мм.
172
Качественная картина распределения скорости движения ча стиц сыпучей среды в нижней зоне аппарата, полученная расчет ным путем, подтверждается опытными данными многих исследо вателей.
Для количественного сравнения результатов аналитического расчета поля скоростей с экспериментальными данными исполь зована работа В. А. Битюкова, в которой определено время дви-
/ — по экспериментальным данным; 2 — по результатам расчета
жения меченых частиц от уровня EF (место сужения потока) до
выхода из модели (рис. 94).
Экспериментальные данные представлены графически в без
размерных координатах в виде зависимости — ф ('5‘)> которая справедлива как для плоских, так и для осесимметрич
ных |
аппаратов. |
|
||
Аналитический расчет поля скоростей проведен для сыпучей |
||||
среды |
со |
следующими физико-механическими |
параметрами^: |
|
V = 3 (упругая изотропная среда); / = 0,577 (хорошо сыпучий |
||||
материал); |
у = |
4200 кг/м3 (металлические шарики). |
||
На рис. |
92 |
и 93 приведены расчетные данные, |
характеризую |
щие распределения вертикальной и горизонтальной составляющих скорости в нижней части аппарата. Кривые распределения вер тикальной составляющей скорости имеют отчетливо выраженные точки перегиба, которые принадлежат изолинии сг2/сгг, 0, ограни чивающей зону стока.
173
Для удобства сопоставления результатов расчета с эксперимен
тальными данными использован безразмерный параметр — , где
Тц
|
hN |
т ц — N |
’ |
Уі uk, о
k=0
hN
N
cos
k=0
Соответственно расчетная формула для вычисления —
вид
N
S Uk, О k=0
N
cosa{ £ ( u l , i + vl, г)1/2 k=0
имеет
(95)
Результаты расчета по формуле (95) представлены на рис. 95
(кривая 2) в виде графика зависимости |
х£/тІІ = ср |
. |
Сопо |
||
ставление |
расчетных данных с экспериментальными показывает |
||||
их удовлетворительную |
сходимость. Максимальное |
отклонение |
|||
т(./тц при |
относительном |
расстоянии |
= 0,45ß |
не |
превы |
шает 15%. |
|
|
|
|
ЗАКЛЮЧЕНИЕ
В данной работе рассмотрены основные результаты теорети ческого и экспериментального исследования механики движения сыпучих материалов в аппаратах, опубликованные в литературе за последние годы. В работе менее подробно освещены вопросы движения газов и жидкостей в зернистых насадках и совершенно не затронуты вопросы теплообмена в аппаратах с неподвижным и движущимся слоем сыпучего материала.
На основе предложенной автором модифицированной модели сыпучей среды выявлены закономерности формирования поля напряжений и деформаций, используя которые удалось объяснить особенности движения зернистых материалов в аппаратах и бун керных устройствах.
Особое внимание уделено обоснованию вывода, что изолинии aJ°z, о и aJ°r, о> построенные с помощью уравнений (38) и (39) при постоянном значении коэффициента прочности среды, позво ляют выявить характерные зоны в сферических аппаратах с воз растающей по радиусу плотностью веществ. При рассмотрении стесненных движений компактных масс квазидискретных твердых тел показана возможность объяснения описанных в литературе геотектонических процессов на основе единого механизма обра зования первичных и вторичных самовозбуждающихся полостей пониженного давления (см. стр. 137).
Вывод о единстве механизмов формирования поля напряжений под действием локальной зоны с повышенной скоростью грави тационного погружения веществ подтверждается известными опыт ными данными о деформациях компактных масс дискретных и квазисплошных твердых тел в различных условиях. К их числу относятся процессы образования:
разгружающего свода высотой С0Х при выпуске сыпучего ма териала из цилиндрической емкости (см. рис. 61);
оболочки толщиной С0СХв зернистом слое с полюсом гравита ционного погружения веществ в центре сферического аппарата
(см. рис. 82);
175
кольцевых элементов шириной С0СХ (см. рис. 82) в квазисплошном диске при сгущении двухфазной смеси «газ—твердое» в компактные глобулы (см. стр. 156).
Относительные размеры названных зон деформируемой среды соответственно равны
§§- = 0,36 (рис. 61), |
= 0,41 (рис. 82). |
Это означает, что изменяя характерный размер ООъ можно определить в соответствии с законом Тициуса—Боде величину С0С1 для системы больших размеров.