Файл: Лукьянов, П. И. Аппараты с движущимся зернистым слоем. Теория и расчет.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 для системы больших размеров.