Математическая модель процесса формирования древесных гранул

МЕХАНИЧЕСКАЯ ОБРАБОТКА ДРЕВЕСИНЫ И ДРЕВЕСИНОВЕДЕНИЕ

УДК 674.08

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ

ПРОЦЕССА ФОРМИРОВАНИЯ ДРЕВЕСНЫХ ГРАНУЛ

© О.Д. Мюллер, канд. техн. наук, доц. В.И. Мелехов, д-р техн. наук, проф. В.К. Любое, д-р техн. наук, проф. Т.В. Тюрикоеа, асп.

Северный (Арктический) федеральный университет имени М.В. Ломоносова, наб. Северной Двины, 17, г. Архангельск, Россия, 163002; е-шаП: ockar@mail.ru

Производство биотоплива из отходов переработки древесины за последние годы приобретает все больший размах. Несмотря на расширение производства древесных гранул, остаются нерешенными вопросы о влиянии отдельных параметров процесса прессования древесных гранул на их качество. Одним из основных параметров, определяющих качество древесных гранул, является давление прессования. Имеющиеся теоретические модели процессов гранулирования показывают, что развиваемое давление прессования является экспоненциальной зависимостью от относительной длины фильеры в матрице. Однако проведенные экспериментальные исследования показали, что эта зависимость носит не экспоненциальную, а линейную зависимость от относительной длины фильеры. Целью настоящей работы была разработка такой математической модели прессования древесной гранулы в цилиндрических фильерах матрицы, которая бы соответствовала экспериментальным исследованиям. Были проанализированы все этапы физического процесса прессования. На основании этого анализа было выявлено, что в процессе прессования древесная шихта претерпевает упругие и пластические деформации. Исходя из этого при разработке новой математической модели процесса прессования древесных гранул за основу была принята теория упруго-пластического деформирования сплошных тел. На базе этой теории рассмотрены два этапа прессования древесной гранулы - выдавливание гранулы через коническое отверстие в цилиндрический канал матрицы и ее движение по цилиндрическому каналу матрицы. Выполненные математические исследования позволили получить зависимость перепада давления по длине гранулы в коническом и цилиндрическом каналах фильеры матрицы и на протяжении всей фильеры в целом. Полученная новая математическая модель процесса прессования древесных гранул показала, что перепад давления в каналах матрицы прямо пропорционален их относительной длине, что полностью соответствует экспериментальным данным.

Ключевые слова: древесные отходы, опилки, древесная шихта, пресс-гранулятор, матрица, фильера, древесные топливные гранулы, пеллеты, биотопливо, давление.

Производство биотоплива из отходов переработки древесины в последние годы приобретает все больший размах. В настоящее время в России уже действуют более 100 заводов по производству древесных топливных гранул, выпускающих около 900 тыс. т древесных топливных гранул (пеллет) в год. При этом около 260 тыс.т произведенных гранул потребляется на внутреннем рынке, остальное количество экспортируется в страны Западной Европы.

Несмотря на все расширяющееся производство древесных гранул, остается нерешенным вопрос о влиянии отдельных параметров процесса прессования древесных гранул на их качество. Информация о теоретических исследованиях процесса прессования древесной шихты весьма ограничена, а разработанные при этом математические модели прессования не всегда соответствуют экспериментальным данным. Например, в работе [1] получена зависимость распределения давления по длине фильеры матрицы, которая имеет экспоненциальный характер от относительной длины цилиндрического канала матрицы:

где р0 - относительное противодавление;

/ - коэффициент трения древесной гранулы о стенки канала матрицы;

Е, - коэффициент бокового давления;

а - угол раскрытия конического участка фильеры матрицы; и (: 3 - длина цилиндрического участка фильеры матрицы до и после конического участка;

12- длина конического участка фильеры матрицы; й и йъ - диаметры цилиндрических участков фильеры матрицы.

В основу данной математической модели были положены условия равновесия выделенного элемента древесной гранулы в канале матрицы под действием сил трения и давления.

Проведенные в Северном (Арктическом) федеральном университете экспериментальные исследования по определению давления выталкивания спрессованной гранулы из цилиндрического канала матрицы показали [2], что распределение носит не экспоненциальный, а линейный характер (рис. 1).

Целью настоящего исследования является разработка математической модели процесса формирования древесных гранул, учитывающей результаты экспериментальных исследований.

Для этого проведены измерения изготовленных древесных гранул [2], анализ результатов которых показал, что наружный диаметр готовых гранул всегда больше диаметра внутреннего канала матрицы.

На основании полученных результатов можно утверждать, что исходная древесная шихта при прессовании древесной гранулы претерпевает упруго-пластические деформации, следовательно, в основу математической модели

Рис. 1. Зависимость давления выталкивания ро гранулы из сосновой шихты (исходная фракция размером 5 < 1 мм) различной влажности V

процесса прессования древесных гранул должна быть положена теория упру-гопластического деформирования.

Для разработки математической модели поделим технологический процесс формирования гранулы на четыре зоны (рис. 2).

Рис. 2. Процесс формирования древесной гранулы (1-1У - зоны): 1- прессовочный ролик; 2 - матрица; 3 - цилиндрический канал фильеры матрицы; 4 - входной конический канал фильеры матрицы; 5 - выходной конический канал фильеры матрицы

В зоне I происходит уплотнение поступившей древесной шихты. На поверхности матрицы образуется слой спрессованной шихты в виде тонкостенного цилиндра (для прессов с цилиндрической матрицей) или плоского диска (для прессов с плоской матрицей). Часть этого слоя при прохождении ролика над входным отверстием фильеры матрицы выдавливается через входной конический канал в калибровочный цилиндрический канал фильеры матрицы. Удаленное количество спрессованной древесной шихты восполняется за счет

прессования поступившей свежей порции древесной шихты. В этой зоне давление прессовочного ролика достигает значения, достаточного для проталкивания спрессованной древесной массы через коническую фаску и калибровочный цилиндрический канал матрицы. Древесная шихта пластически деформируется до максимальной для этого давления плотности. После прохождения ролика оставшаяся часть спрессованной шихты упруго расширяется.

В зоне II (в конической фаске цилиндрического отверстия матрицы) формируется цилиндрическая гранула и происходит истечение спрессованной древесной шихты в цилиндрический канал матрицы. Давление в этой области меньше давления прессования под роликом, вследствие чего находящийся здесь объем спрессованной древесной шихты подвергается упругопластиче-ским деформациям без изменения ее плотности. В этой зоне за счет пластических деформаций спрессованная древесная шихта переходит от формы усеченного конуса к форме цилиндра.

Зона III представляет собой цилиндрический канал матрицы, в котором стабилизируется и закрепляется цилиндрическая форма древесной гранулы, а спрессованная древесная шихта испытывает только упругие деформации. Протяженность цилиндрического канала матрицы во многом определяет давление прессования (давление выталкивания) p0 и, следовательно, плотность спрессованной древесной гранулы.

Зона IV - это зона выхода спрессованной гранулы из цилиндрического канала матрицы. Так как спрессованная древесная цилиндрическая гранула находится под действием давления упругих деформаций, то в момент выхода древесной гранулы из цилиндрического канала матрицы боковое давление со стороны матрицы на нее снижается до атмосферного, происходит увеличение наружного диаметра древесной гранулы на величину упругих деформаций. Резкое увеличение наружного диаметра гранулы может привести не только к образованию трещин, но и даже к разрушению гранулы. Во избежание этого на выходе из цилиндрического канала матрицы выполнен выходной конический канал (5), обеспечивающий плавное снижение упругих деформаций по мере продвижения гранулы.

Из сказанного выше следует, что давления прессование p0 определяется перепадом давлений на входном коническом и цилиндрическом участках фильеры матрицы. Так как основными процессами формирования древесной гранулы являются упругопластические деформации, то для определения перепада давлений на указанных участках воспользуемся основными положениями теории упругих и пластических деформаций [3].

Определение перепада давления в коническом канале матрицы

На рис. 3 представлена расчетная схема для определения усилий при проталкивании древесной шихты через конический канал в цилиндрический канал матрицы. Процесс выдавливания древесной шихты происходит в обла-

сти пластических деформаций, при этом соблюдаются условия постоянства объема.

Учитывая тот факт, что канал представляет собой усеченный конус, за расчетную принята сферическая система координат г, 0, ф.

За верхнюю границу объемно-пластической области примем часть поверхности сферы радиусом Я1 с углом при вершине конуса 2у, опирающуюся на входное отверстие конической фаски диаметром В, за нижнюю - часть поверхности сферы радиусом Я2 с углом при вершине конуса 2у, опирающуюся на выходное отверстие конического канала диаметром й0. Боковая поверхность объемно-пластической области образуется конической поверхностью канала.

Рис. 3. Расчетная схема для конического канала матрицы: О — диаметр сечения конического канала на входе; — диаметр сечения конического канала на выходе; й — диаметр текущего сечения конического канала на радиусе г; Я — радиус телесного угла, опирающегося на дугу диаметром О; Я — радиус телесного угла, опирающегося на дугу диаметром ; 2у — угол при вершине конического канала; р — давление на входе в конический канал; р2 — давление на выходе из конического канала

Давление выталкивания р во входном сечении выделенного элемента определим на основании энергетического баланса сил, действующих на данный элемент.

Выберем промежуток времени А?, за который давление р совершит работу выталкивания АЫт. Эта работа будет затрачена на совершение работы сил давления на выходе из конического канала Аых, работу сил трения о боковую поверхность Аи работу пластической деформации древесной шихты

в коническом канале Ал. Запишем уравнение энергетического баланса для выделенного элемента:

Ал + Аых + А = 0.

Разделив правую и левую части уравнения (2) на время А?, получим уравнение баланса мощностей для выделенного элемента, которое и будет основным при определении давления на входе в конический участок р :

N + N + N + N = 0

1 'выт пл ~ 1 'вых тр V~v

Определение мощности для пластического деформирования выделенного элемента

Уравнения равновесия элементарного объема в сферических координа-

тах:

да

- + -

1 Этг0 1 ЭХ;

+--

Гф

1

я  • яа  я +" [2аг-(ае+аф)+т е^ф]=0; дг гsmф Эе г Эф г 1 4 ф ф *

Эх,

ег

дг

+

дг,

фг

+

1 дае 1 дтеф 1 гsmф де г дф г1 ф *

1 дТфе 1 дат 1

+

дг гsinф де г дф г

71^+7 [3тфг-(oe-oф)ctgФ]=0.

Здесь аг, ае и сф - главные напряжения, действующие в направлении осей г, е и ф; т^,т и тЭф - касательные напряжения в плоскостях, нормальных к осям г, е и ф.

Уравнения связи между скоростями течения и скоростями деформаций в сферических координатах:

дvr

Ке =

1 дvQ VфCtgф

гапф дв

+

+ —; г

г = + ф г дф г

1 д^ дvt

Лге =

Леф =

rsinф дв дг vеctgф

г дф

1

г

+

1

ф

Зv(

ЛфГ =

ф

г rsinф дв 1 дvr

дг

г + г дф

Здесь ^еи £ - скорости относительных деформаций, действующие

в направлении осей г, е и ф; т^и п9ф - скорости относительных угловых

деформаций в плоскостях, нормальных к осям г, е и ф; уф, уе и уг - скорости течения элементарного объема в направлении осей г, е, ф.

Запишем условие постоянства объема (условие несжимаемости):

+ ^ + £ = 0.

Полагая, что движение древесной шихты через конический канал осе-симметричное, в силу симметрии потока получим

э е

= 0;

э ф

= 0;Уе = 0; Уф = 0.

V

V

Г

V

С учетом (7) системы уравнений (4) и (5) принимают следующий вид:

( Заг 1

+ - [2ar - (ае + + тгфс^Ф] = 0;

Зг г Зтег 1

ЗТфГ 1

Зг г

+ - [3т0г + 2тeфctgф] = 0;

V ar + - [3тфг - (ае - аф)с£ф] = 0

dvr

= з7;

vr

^е = ^ф = —; . Лге = Леф = Лфг =

Определим выражение для скорости движения древесной шихты в коническом канале в радиальном направлении. По условию постоянства объема секундный расход является величиной постоянной в любом радиальном сечении конического отверстия.

Пусть радиальная скорость движения шихты на выходе в цилиндрический канал (г = Я{) равна у1 . Площадь поверхности шарового сегмента с углом при вершине 2 у и радиусом Я [4]:

S = пЯ2 (2 + sin2у - 2cosy). Тогда уравнение постоянства расхода

q = vS = vпЯ2 (2 + sin2у - 2cosy) = const.

Для произвольного сечения на радиусе г уравнение (11) примет следующий вид:

q = V

\пЯ?(2

+ sin2y - 2cosy) = v пг2 (2 + sin2 у

(2 + sin2у - 2cosy).

Отсюда определим радиальную скорость (знак «минус» учитывает встречное движение по оси г):

Я2

К = -V

1 2 • г

После подстановки (13) в (9) и преобразований получим

R¡

^Г = 2ví-¿3

R¡

5e = -vi7j;

^Ф = -vi^s".

Следовательно, условие несжимаемости (6) выполняется.

Дальнейший расчет пластической деформации ведем методом энергетического баланса. С одной стороны, приложенные к выделенному объему внешние силы совершают работу Авыт по выталкиванию его через коническую фаску в цилиндрический канал, с другой - работу его пластического деформирования Апл от формы усеченного конуса до формы цилиндра.

Для определения работы пластического деформирования на радиусе r выделим элементарный элемент сферической формы с углом при вершине 2у и толщиной dr. Элементарный объем выделенного элемента с точностью до бесконечно малых величин

dV = nr2

Элементарная мощность пластической деформации на выделенном элементе определяется как

Nm = o^dV,

где o - среднее по объему значение напряжения течения спрессованной древесной шихты; - интенсивность скоростей деформации,

^ =^2 - ^е)2 +

После подстановки (9) и (14) в (17) получим

R2

^ = 2vi —.

r

Работа пластической деформации

Nm = {dNnjI.

V

После подстановки (15), (16) и (18) в (19) будем иметь

Nm = {2oxVj —dV = -2oxVj —L[2(1 - cosy) + sin2y]dr

или

Nm = -2яоЛ—2 [2

—2

Определение мощности сил выталкивания во входном сечении конического канала

Мощность сил выталкивания, действующих на выделенную сферическую поверхность во входном сечении, определим по зависимости

NBbn. = М =

Определение мощности сил давления на выходе из конического канала

Мощность сил давления в выходном сечении запишем по аналогии с входным сечением:

^ых = F2V2 = Р2S2

Определение мощности сил трения древесной шихты о боковую поверхность конического канала

Элементарную силу трения на боковой поверхности конического отверстия определим по зависимости

dFф = -TdS = -т2пгйг.

Приняв т = Рах (где Р = kv; k - коэффициент трения; v - коэффициент Пуассона), получим

dFф = -TdS = -2лРох rdr.

Мощность сил трения определим путем интегрирования элементарной силы трения dF^,:

Я2

Nw =\vdFw = J 2nPaTvrdr

S R1

или с учетом (13):

NTp = J - 2пРа^г = - 2nPoTV1R12ln R.

Я2 R

Определение давления на входе в конический канал

Для определения давления в древесной шихте на входе в конический канал подставим найденные значения мощностей (20) - (22) и (26) в (3):

npvR2 [2 (1 - cosy) + sin2 Y]-2noTvR2 [2 (l - COSY) + sin 2Y] ln — -

Я2

np2vR [2

Я2

Отсюда

V2R-2 {1 + ^ \ ,2 1 o> RL.

Pl = Р2^± + 21 1 + -ZT,-V • 2 I ^т1^

vR ^ 2(1 - cosy) + sin y) R

Р = Р2 + 21 1 + ^-V • 2 I .

^ 2

По уравнению неразрывности (11) имеем vR = v2R2, тогда

L, ß V, nR1

R2

При этом

R2 d0

После подстановки (30) в (29) и преобразований получим уравнение для расчета перепада давления в коническом канале матрицы:

= Pl - Р2 = 2 \ 1 +^-V  2 ^ С> D.

^ 2(1 - cosyj + sin у) d0

Расчет цилиндрического канала матрицы

При движении по цилиндрическому каналу спрессованная древесная гранула имеет форму цилиндра и практически не изменяется. Таким образом, на этом этапе прессования древесной гранулы отсутствуют пластические деформации. В то же время древесная гранула после выхода из конического канала и входа в цилиндрический канал находится в упругодеформированном состоянии (в состоянии «натяга»). Диаметр упругодеформированной гранулы равен диаметру цилиндрического канала do (рис. 4). После выхода из цилиндрического канала за счет упругого расширения диаметр гранулы увеличивается до Djp.

Рис. 4. Схема упругодеформированного состояния гранулы в матрице

Расчетная схема для определения упругодеформированного состояния гранулы представлена на рис. 5. Выберем цилиндрическую систему координат г, в иг.

Рис. 5. Расчетная схема упругодеформированного состояния гранулы в цилиндрическом канале (ог — нормальное напряжение вдоль оси г; о0 — нормальное напряжение вдоль оси 0; Д0 — центральный угол выделенного элемента; г — радиус древесной гранулы в фильере; — радиус древесной гранулы после

выхода из фильеры; иг — смещение конца радиуса наружной поверхности древесной гранулы после выхода из канала фильеры матрицы)

Выделим элементарный объем гранулы и для этого объема запишем основные уравнения теории упругости [5].Так как древесная гранула имеет форму цилиндра, то в качестве расчетной выберем цилиндрическую систему координат г, 0 и г. Исходные уравнения теории упругости в цилиндрической системе координат имеют следующий вид: вектор перемещений

й =

уравнения равновесия

( даг

дг дтгв

дг

+

1дъ

гв

Г

г дв

+

дтг

дг

+ Фг = 0;

2тгв 1 дав дтвг + —^ + —ттт + -Г1 + Фв = 0;

г

г дв

дг

^ дг Г

1 дтв7 да7 г дв дг

закон Гука

аг = Аг +2 ргг; ав = Аг + 2^гв; аг = Аг + 2 у.гг;

тгв = 2№ггв; тгг 2 №ггг;

чТвг = 2№вг> где X и ц - параметры Ламе, А =

vE

V - коэффициент Пуассона; Е - модуль упругости;

; & =

Е

2(1+ч) '

уравнение деформаций

е = е.

Выразим деформации через смещения:

Зиг

£г = "д7;

1иг

~1г+ ;

г ог г

£е

£Г0

£02

оиг ои8 и8\

2 \г 90 Зг г / ; 1 /Зиг Зи2\

= 2\з7 + ~з7/;

1 /Зи8 1 Зи2\ = 2\37 + Г~307'

Для условий нашей задачи, когда цилиндрическую гранулу можно считать бесконечно длинной, частная производная — = 0 . Смещения вдоль оси г

&

равны 0, иг = 0, тогда поставленная задача сводится к решению плоской задачи.

Так как гранула в поперечном сечении имеет форму круга, то в силу осевой симметрии выполняются следующие условия:

смешения не зависят от полярного угла 0, т. е. — = 0 ;

30

тангенциальное смещение щ = 0.

Из массовых сил на гранулу действуют только силы тяжести, которыми ввиду их малости можно пренебречь, тогда компоненты массовых сил равны 0, т. е. Фг ,Ф0 и Фг. Системы уравнений (33)-(36) принимают следующий вид: уравнения равновесия

С ЗОг

От

а8

+ ---= 0;

Зг Зхг8

+

г

2тг8

= 0;

закон Гука

Зг г

+ _ д

Зг г

ог = Аг + 2 цгг; о8 = Аг + 2цг8; о2 = Аг; тГ8 = 0;

ТГ2 0;

-Т82 = 0;

г

гг

уравнение деформаций

е = ег + е0.

Выразим деформации через смещения:

бы

г бг '

е = ^ е0 • г

После подстановки (40) в (38) и (39) и преобразований получим

О =

бг г

О0 = X ^ +

бг г

бы Хы .„„.

е = —^ + —^

бг г

Подставив (41) и (42) в (37), получим дифференциальное уравнение равновесия выделенного элемента гранулы:

^ +1-- "Ыг = 0.

бг г бг г

Для нахождения решения дифференциального уравнения (44) необходимо определиться с начальными и граничными условиями. Так как прессование древесной гранулы идет с некоторыми допущениями фактически непрерывно, то процесс можно считать стационарным. В этом случае для решения поставленной задачи необходимо задаться только граничными условиями. Здесь их два: одно на поверхности цилиндрической оболочки недефор-мированной гранулы при г = Ятр, другое - на поверхности цилиндрической оболочки гранулы после ее упругой деформации до диаметра цилиндрического канала матрицы при г = г0. Граничные условия:

при г = г0 = имеем ог = -р,

Ар

при г = = имеем ог = 0.

Найдем решение дифференциального уравнения (44) в виде

Ыг = г».

Дважды продифференцируем:

бЫ (»"О /-/104

= пг

г

бг

б2ы

бг2

После подстановки (48) и (49) в (44) получим

1 Vй

г г

или

Отсюда

п2 -1 = 0.

Решением данного уравнения являются корни ^ = +1 и п2 = -1, а решением дифференциального уравнения

иг = С1 г + С2 г-1,

где С1 и С2 - константы интегрирования, которые определяются из граничных условий.

Для нахождения констант интегрирования продифференцируем уравнение (53):

ди ^ С,

= С1 - с2.

дг г

Подставляя (53) и (54) в (41) и (42), получим после преобразований:

С

Ог =

г

С

О =

г

Для нахождения констант интегрирования воспользуемся граничными условиями: из условия

С

откуда

С =ъгЧ §-;

(х + ц) ^р

из граничного условия (45) имеем

С

г0

откуда с учетом (58) после преобразований

С2 = -Р.

(^ - го2) V После подстановки (60) в (58):

С =_г__р__

С1 =

Подставив найденные значения констант интегрирования С и С в уравнения (55) и (56), получим выражения для определения радиальных и тангенциальных напряжений при упругой деформации круглого длинного стержня:

г2 (Яр ^

Перемещение иг после подстановки (60) и (61) в уравнение (53) может быть определено по следующей зависимости:

Ы = р , гг02 , | — +1 4|.

г Р 2

Выразив коэффициенты Ламе через модуль упругости Е и коэффициент Пуансона V и подставив в (64), получим

(1 + V) гг2 ( Яр ^

ыг=pl-r

Уравнение (65) является основным для экспериментального определения давления р со стороны древесной гранулы на боковую стенку в цилиндрическом канале матрицы. Положив г = г0, преобразуем уравнение (65):

Радиальное смещение для гранулы при запрессовке ее в фильеру

ыг = Яр - г.

После подстановки (67) в (66) получим выражение для определения бокового давления спрессованной гранулы на боковую цилиндрическую стенку фильеры:

- (*- "0<^ ' ^ г \ 1 - 2v +

или через диаметры

/ / Д2, ^

Рб =

"0 (1 - *+а J

Сила трения спрессованной древесной гранулы о цилиндрическую стенку канала матрицы может быть определена как

Р1р = ш10£крб.

Сила выталкивания гранулы из цилиндрического канала матрицы

^выт = ^ АРц,

где Арц - перепад давления между входным и выходным сечениями цилиндрической части древесной гранулы. Приравняв силы трения (70) и выталкивания (71), определим перепад давления на цилиндрической грануле:

кР

Арц=4-Рб.

"0

После подстановки (69) в (72) получим

4Р, = 4—т-г-----.

ц "0 (7 + V) К J

(Ар -do)l J2

del1" +

Определение давления прессования

На выходе из цилиндрического канала матрицы спрессованная древесная гранула находится под атмосферным давлением, которое можно принять за ноль. В этом случае давление прессования, развиваемое прессовочным роликом,

^рол = АР + АРц.

После подстановки в (74) выражений для Арц из (73) и Арк из (31) и преобразований получим

( d 1

/ Ч (Dгр - do-11

J, р 1 , D k£ E v \d0 j

= + ^- V • ^ 1°>Т + 4Т7ТТЛ—7-

р ^ 2 (1 - cosy) + sin y) d0 d0 (1 + v) ( D2 1

del 1 - 2v + -f

или с учетом, что Р = kv:

(^

/ ч К - ^)|-гр - 11

"о [,1 - 2v + ^ J

Из полученной зависимости следует, что давление, развиваемое роликом для проталкивания спрессованной древесной гранулы, имеет линейный характер. В то же время на давление оказывают влияние такие параметры спрессованной древесной шихты, как модуль Юнга , коэффициент Пуассона и коэффициент трения древесной гранулы о цилиндрические стенки матрицы.

Выводы

1. Полученная на основании теории упруго-пластического деформирования математическая модель прессования древесных гранул адекватна результатам экспериментальных исследований.

2. Для практического использования полученной математической модели при расчете технологических процессов и проектировании оборудования для производства древесных гранул необходимы масштабные экспериментальные исследования, чтобы создать базу данных для спрессованной древесной шихты по модулю Юнга , коэффициенту Пуассона и коэффициенту трения древесной гранулы о металл матрицы, а также по влиянию на них фракционного состава древесной шихты и породы древесины, влажности шихты, плотности спрессованных гранул и температуры процесса гранулирования.

СПИСОК ЛИТЕРАТУРЫ

1. Булатов А.С. Разработка процесса прессового гранулирования мелкодисперсных сред на примере минеральных порошков и древесных отходов: автореф. ...канд. техн. наук. М., 2012. 25 с.

2. Исследование и разработка методов совершенствования технологического процесса получения гранулированного топлива из древесного сырья как возобновляемого источника энергии: отчет о науч.-исслед. работе по гос. контракту № 14.740.11.0097 от 08 сент. 2010 г. (промежуточный, этап № 4) / САФУ; рук. В.К. Любов. Архангельск, 2012. 150 с.

3. Овчинников А.Г. Основы теории штамповки выдавливанием на прессах. М.: Машиностроение, 1983. 200 с.

4. Бронштейн И.Н., Семедяев К.А. Справочник по математике для инженеров и учащихся втузов. М.: Наука, 1981.723 с.

5. Самуль В.И. Основы теории упругости и пластичности. М.: Высш. шк., 1982.

264 с.

Поступила 23.01.14

UDS 674.08

Mathematical Model of Wood Granules Formation

O.D. Myuller, Candidate of Engineering, Associate Professor V.I. Melekhov, Doctor of Engineering, Professor V.K. Lyubov, Doctor of Engineering, Professor T. V. Tyurikova, Postgvaduate Student

Northern (Arctic) Federal University named after M.V. Lomonosov, Naberezhnaya Severnoy Dviny, 17, Arkhangelsk, 163002, Russia; e-mail: ockar@mail.ru

Production of biofuel from waste of processing of wood in recent years gains the increasing scope. Despite all extending production of wood granules, there are unresolved questions of influence of separate parameters of process of wood granules pressing on their quality. One of the key parameters, defining the quality of wood granules pressing pressure is. Available theoretical models of processes of a granulation show that the developed pressure of pressing is exponential dependence on relative length of die hole in a matrix. The conducted pilot researches have shown that this dependence carries not exponential dependence on relative length of die hole, but linear. Development of such mathematical model of a wood granule pressing in cylindrical die holes of a matrix which would correspond to pilot studies, was the purpose of the real work. All stages of physical process of pressing were analyzed. On the basis of this analysis it was revealed that in the course of pressing wood furnace charge undergoes elastic and plastic deformations. Proceeding from it when developing new mathematical model of process of wood granules pressing the theory of elasto-plastic deformation of continuous solid was assumed as a basis. On the basis of this theory two stages of a wood granule pressing - process of expression of a wood granule through a conic opening to the cylindrical canal of a matrix and process of movement of a granule on the cylindrical canal of a matrix are considered. The executed mathematical researches allowed to receive dependence of pressure difference on granule length in conic and cylindrical canals of a die hole in matrix and on all die hole as a whole. The received new mathematical model of process of wood granules pressing has shown that pressure difference in channels of a matrix is directly proportional to their relative length, that completely corresponds to experimental data.

Keywords: wood wastes, sawdust, a wood charge, a press-granulator, a matrix, a die hole, wood combustive granules, pellets, biofuel, pressure.

REFERENCES

1. Bulatov A.S. Razrabotka protsessa pressovogo granulirovaniya melkodispersnykh sred na primere mineral'nykh poroshkov i drevesnykh otkhodov: Avtoref.kand. diss. [The Development of the Press-Granulation Process of the Fine-Dispersed Fluid by the Example of the Mineral Powders and the Wood Waste: Cand. Diss. Abs]. 2012. 25 p.

2. Lyubov V.K. Issledovanie i razrabotka metodov sovershenstvovaniya tekhnolog-icheskogo protsessa polucheniya granulirovannogo topliva iz drevesnogo syr'ya, kak vozobnovlyaemogo istochnika energii [The Research and the Development of Methods of Improving the Process of Granular Fuel from the Wood as the Renewable Energy Source]. Arkhangelsk, 2012. 150 p.

3. Ovchinnikov A.G. Osnovy teorii shtampovki vydavlivaniem na pressakh [The Fundamentals of the Theory of the Extruding Presses]. Moscow, 1983. 200 p.

4. Bronshteyn I.N., Semedyaev K.A. Spravochnik po matematike dlya inzhenerov i uchashchikhsya vtuzov [The Guide of Mathematics for Engineers and Students of Technical Institutes]. Moscow, 1981. 721 p.

5. Samul' V.I. Osnovy teorii uprugosti i plastichnosti [The Fundamentals of the Elasticity and Plasticity Theory]. Moscow, 1982. 264 p.

Received on January, 23, 2014