УДК 674.047
DOI: 10.17238/issn0536-1036.2016.2.114
ИССЛЕДОВАНИЕ КИНЕТИКИ ПРОЦЕССА СУШКИ ПРОФИЛЬНЫХ ЗАГОТОВОК ИЗ ДРЕВЕСИНЫ
П.В. Билей, д-р техн. наук, проф. Р.А. Рокунь, асп.
Национальный лесотехнический институт Украины, ул. Генерала Чупринки, д. 103, г. Львов, Украина, 79057; е-mail: [email protected], [email protected]
Сушка профильных заготовок из массивной древесины имеет ряд особенностей. В первую очередь - это переменное сечение по длине заготовки, что ставит ряд сложных задач по формированию штабелей или пакетов, выбору режимов, технологий и сушильного оборудования. В качестве экспериментального материала выбрана древесина бука. Процесс сушки является сложным тепломассообменным процессом, включающим теплообмен между поверхностью древесины и средой, тепло- и влагопро-водность внутри материала и испарение влаги на его поверхности. Для изучения тепломассообмена необходимо знать распределение влаги по сечению материала. В основном принимают для исследований или расчетов, что влага внутри материала распределяется по параболической или косинусоидальной зависимости. Косинусоидаль-ная зависимость описывается уравнением скорости сушки в виде баланса влаги, которая перемешается внутри материала и испаряется с его поверхности. Полученная физико-математическая модель является исходным положением для составления методики экспериментальных исследований кинетики процесса сушки. Экспериментально необходимо определять влажность центральных и поверхностных слоев древесины, изменение средней влажности во время сушки, линейные размеры заготовок, режимные параметры процесса сушки (температура среды и равновесная влажность древесины). Представлены формулы, по которым определены коэффициенты сушки, вла-гопроводности, влагоотдачи, массообменные критерии Нуссельта и Фурье. По экспериментальным данным построены кривые сушки, отражающие изменение влажности древесины во время сушки, и кривые скорости сушки, на основании которых определяют кинетические характеристики процесса. Найдено, что скорость сушки в ходе эксперимента уменьшается более чем в шесть раз. Полученные значения коэффициентов сушки и влагопроводности можно использовать в теоретических уравнениях для определения продолжительности процесса сушки. Обобщающими характеристиками кинетики процесса сушки являются безразмерные критерии Нуссельта (характе-ризирует соотношение коэффициентов влагоотдачи материала с характерным размером и влагопроводности) и Фурье (характеризирует интенсивность перемещения влаги внутри древесных сортиментов).
Ключевые слова: древесина, влажность, коэффициенты сушки, влагопроводность, влагоотдача, критерии, кинетика.
Введение
В лесной, деревообрабатывающей и мебельной промышленности часто используют профильные заготовки из древесины. Это могут быть столярные изделия, декоративные элементы мебели, деревяные элементы стрелкового оружия, элементы деревяных музыкальных инструментов, деревяные протезы, обувные колодки и дрова. В последние время имеется большой спрос на сухие дрова, идущие на экспорт и внутренний рынок. Режимам и технологии сушки профильных заготовок из древесины уделено очень мало внимания в специальной технической литературе. Профильные заготовки имеют переменное сечение по длине, поэтому возникают различные сложности по способам укладки в пакеты или штабеля, выбору режимов, технологий и сушильного оборудования.
В данной работе рассмотрен конвективный процесс сушки буковых сортиментов. Бук европейский (Fagus sylvatica L.) является наиболее распространенной породообразующей древесиной Карпатского региона Украины [1].
Теоретические положения
Внутренний влагоперенос и внешний влагообмен поверхности древесины с окружающей средой являются переменными величинами во времени и пространстве. Исследование конвективной сушки древесины твердых пород (бук, дуб, граб) показали, что более точно распределение влагосодержания в материале можно описать косинусоидой [1]. При трехмерном потоке влаги (например, для заготовок или дров) закон распределения влажности можно описать следующим уравнением:
КХ KV KZ
U = UnoB + (уц -UnoB)cos—■ cos■ cos—,
где £/пов, иц - соответственно влагосодержание поверхностных и центральных слоев;
R1, R2, R3 - соответственно толщина, ширина и длина заготовки, см.
Среднюю по объему влажность древесины определим по формуле
R RR.>
1 1Н 1 у2 1 у3
W =-J J \udx-dy-dz.
RiR2R3 0 0 0
--RAA = oj—R,^ +—R2R3+—RA) = P(f/n0B-UMR2 +R2R3 +R1R3),
Скорость сушки в данном случае найдем, используя уравнение
7_123 т 12 23 -л
d х ах ду ог
где ат - коэффициент влагопроводности; Р - коэффициент влагоотдачи;
ир - равновесное влагосодержание древесины [1, 3, 4, 6-8]. В конечном варианте запишем уравнение скорости сушки для периода убывающей скорости:
с +
п
ям
К1 К3
1 1 1 : — + — + —;
Я1 Я2 Я3
где Яс - обобщенный характерный размер заготовок, Яс = Я / Яа
С, П - постоянные коэффициенты, С = 4/п2 = 0,405; П = 1,0;
Н - соотношение коэффициентов влагоотдачи рс и влагопроводности ат [2, 6];
Ж и Жр - соответственно средняя и равновесная влажность древесины.
Таким образом, для экспериментальных исследований кинетики процесса сушки необходимо иметь данные о средней влажности древесины (Ж), влажности поверхностных (Жпов) и центральных (Жц) слоев древесины, линейных размерах заготовок (Я1, Я2, Я3), режимных параметрах процесса сушки (температура и равновесная влажность древесины) [2-9].
Методика экспериментального исследования
Как правило, поперечное сечение дров имеет вид треугольника или трапеции. Геометрический центр этих фигур находят по пересечению биссектрис или диагоналей. Для измерения влажности центральных слоев древесины датчики электроконтактного влагомера необходимо внедрить на глубину х = Я. Для измерения влажности поверхностных слоев датчики достаточно внедрить на глубину х = 0,15Я.
Перед процессом сушки (и в процессе сушки) контрольные образцы взвешивают - находят начальную (М„) и среднюю текущую ) массу образца (заготовки). Потом в середине образца выпиливают секции послойной влажности и секцию влажности, которые взвешивают и помещают в сушильный шкаф, где высушивают до абсолютно сухого состояния. Таким образом весовым способом находят распределение влажности по толщине заготовки в конце процесса сушки, а также массу (М0, кг), которую бы имела заготовка в абсолютно сухом состоянии:
100М
М0 =-=4
100 +Ж/
где Wн - средняя начальная влажность образца.
По массе заготовки в абсолютно сухом состоянии (М0) можно воссоздать кинетику изменения средней влажности древесины во время сушки:
м
1
За основу были приняты известные режимы сушки [3, 4].Фиксировались продолжительность процесса сушки на каждой из ступеней режима (т, с), интервал времени (Дт, с) и изменение средней влажности древесины (А Ш, %). Средняя начальная влажность заготовок Шн = 62,4 %. Эквивалентный радиус заготовок Я = 5,2 см.
Результаты экспериментальных исследований
В табл. 1 приведены данные экспериментов по изменению влажности древесины в процессе сушки, а также режимные параметры (ступень режима, равновесная влажность древесины (Щ,), изменение влажности древесины (среднее ее значение (Ш), влажность поверхностных (Шпов) и центральных (Шц) слоев древесины, продолжительность процесса (т)).
Таблица 1
Изменение влажности древесины в процессе сушки
Параметры Продолжительность Влажность древесины Скорость
процесса, % сушки в процессе сушки, % сушки
Ш Шр т-104, с Дт 104, с Ш Шц Ш гг пов ДШ ^-КГ* % Дт ' с
0 0 62,4 68,2 45,8 - -
Шп...40 13,5 5,76 5,76 39,9 49,0 28,6 22,5 3,91
40...30 9,0 10,80 4,32 30,3 36,6 18,9 9,6 2,22
30...20 7,0 16,20 6,12 20,0 25,4 12,7 10,3 1,68
20.15 4,5 25,20 9,0 14,8 15,6 8,2 5,2 0,58
Кинетические характеристики процесса сушки
Таблица 2
Параметры процесса, % Коэффициенты Критерии
Ш Шр К,1/с ат, см2/с Рс, см/с ^ит РОт
Шн...40 13,5 10,710-6 4,3810-4 5,4110-5 0,64 0,93
40.30 9,0 8.61-10"6 2,87-1(Г4 4,42-10-5 0,80 0,46
30.20 7,0 9.53-10"6 3.0210"4 4,80-10-5 0,83 0,68
20.15 4,5 4,54-Ю"6 1,79-Ю"4 2,94-10-5 0,81 0,60
Средние значения 8.35-10"6 3.0210"4 4,39КГ' 0,77 0,68
В табл. 2 приведены кинетические характеристики процесса сушки, определенные по следующим формулам: а) коэффициент сушки
г щ-к
где т - продолжительность 7-й ступени режима сушки, с;
Шн, - соответственно средние начальная и конечная влажности древесины на данный ступени режима, %; Шр - равновесная влажность древесины на данной ступени режима, %;
ßc= — R/(Wt-Wp);
б) коэффициент влагопроводности
„
т &z Дг
г
где A W /Атг- - скорость сушки на данной ступени режима, %/с;
W - средняя влажность древесины на данной ступени режима, %; R - характерний размер (эквивалентный радиус), см; AW/Ax - перепад влажности по толщине материала, %/см;
в) коэффициент влагоотдачи
A W „ . — Ах
г) массообменный критерий Нуссельта, характеризирующий соотношение влагопроводности древесины и влагообмена поверхности древесины со средой,
Num=M/am;
д) массообменный критерий Фурье, характеризирющий интенсивность
перемещения влаги внутри древесных сортиментов
Fom =ax/R2.
т т
Таким образом, построенные по экспериментальным данным кривые сушки (изменения влажности древесины во время сушки) и кривые скорости сушки дают возможность определить кинетические характеристики процесса в целом. Скорость сушки по ходу процесса уменьшается в 6,7 раз. Полученные значения коэффициентов сушки и влагопроводности можно использовать для определения длительности процесса. Обобщающими характеристиками процесса сушки являются безразмерные комплексы - критерии Нуссельта и Фурье.
СПИСОК ЛИТЕРАТУРЫ
1. Билей П.В. Сушка древесины твердых лиственных пород. М.: Экология, 1992. 224 с.
2. Кейс В.М. Конвективный тепломасообмен: пер. с англ. М.: Энергия, 1972.
448 с.
3. Лыков А.В., Михайлов Ю.А. Теория тепло- и массопереноса. Л.: Госэнергоиз-дат, 1963. 535 с.
4. Лыков А.В. Тепломасообмен: справ. 2-е изд.; пераб. и доп. М.: Энергия, 1978. 480 с.
5. Сушка древесины: справ. Серия «Деревообработка» / Сост. С.Ж. Алюшин, В.Ю. Березовский, П.В. Билей, С.А. Бурышев, Н.В. Приходько, С.А. Стрелков. К.: ТМА «Тристан», 2004. 440 с.
6. Шубин Г.С. Теория тепловой обработки и расчет процессов сушки древесины. М.: Лесн. пром-сть, 1973. 248 с.
7. Чудинов Б.С. Вода в древесине. Новосибирск: Наука, 1984. 270 с.
8. Вшей П.В. 1.В. Петришак, I.A. Соколовський, Л.Я. Сорока. Тепломасообмшш процеси деревообробки деревини (Шдручник). Л^в. 2013. 376 с.
9. Brunner - Hildebrand. Die schnitholztrocknung (S. Auflag): Buchdruckwerkstat-ten. Hannover. Gm.bH, 1987. 322 s.
Поступила 08.02.15
UDC 674.047
DOI: 10.17238/issn0536-1036.2016.2.114
Kinetics of Drying of Wood Profile Blanks
P. V. Biley, Doctor of Engineering Sciences, Professor R.A. RokunPostgraduate Student
National Institute of Forestry of Ukraine, General Chuprynka str., 103, Lviv, 79057, Ukraine; e-mail: [email protected], [email protected]
The drying process of profile blanks of solid wood has a number of features. First of all it is a variable cross-section along the length of a blank, which poses a number of challenges to form the staples or packages, to select the modes, technologies and drying equipment. The beechwood is selected as an experimental material, which is often used for firewood. The drying process is a complex process of heat and mass transfer which includes the heat exchange between the wood surface and the medium, thermal conductivity and moisture conduction inside the material, evaporation of moisture on its surface. To study the heat and mass transfer it is necessary to know the distribution of moisture in the cross section of material. Basically, the nature of moisture distribution inside the material is considered for the research in the form of a parabolic or cosine dependence. The cosine dependence is described by the drying rate law in the form of moisture balance. Moisture evaporates from the material surface. The resulting physical and mathematical model is the starting position for the methodology of the kinetics experimental studies of drying process. The moisture content of the central and surface layers of wood, change in the average moisture content at the drying process, linear dimensions of blanks, operational parameters of the drying process (the environmental temperature and the equilibrium moisture content of wood) are determined experimentally. The paper presents formulas, determining the coefficients of drying, moisture conduction, moisture-yielding ability, mass exchange Nusselt and Fourier numbers. According to the experimental data the drying curves reflecting the change in the moisture content during the drying process and the drying rate curves, which help to determine the kinetic characteristics of the process, should be constructed. The drying rate during the test is reduced by more than six times. The obtained values of drying and moisture coefficients can be used in the theoretical equations to determine the duration of the drying process. Kinetics general characteristics of drying process are dimensionless Nusselt and Fourier numbers.
Keywords: wood, moisture content, drying coefficient, moisture conduction, moisture-yielding ability, criteria, kinetics.
REFERENCES
1. Biley P.V. Sushka drevesiny tverdykh listvennykh porod [Wood Drying of Solid Hardwood]. Moscow, 1992. 224 p.
2. Kays W.M. Convective Heat and Mass Transfer. New York, 1966.
3. Lykov A.V., Mikhaylov Yu.A. Teoriya teplo- i massoperenosa [The Theory of Heat and Mass Transfer]. Leningrad, 1963. 535 p.
4. Lykov A.V. Teplomassoobmen [Heat and Mass Transfer]. Moscow, 1978. 480 p.
5. Alyushin S.Zh., Berezovskiy V.Yu., Biley P.V., Buryshev S.A., Prikhod'ko N.V., Strelkov S.A. Sushka drevesiny [Drying of Wood]. Kiev, 2004. 440 p.
6. Shubin G.S. Teoriya teplovoy obrabotki i raschet protsessov sushki drevesiny [Theory of Heat Treatment and Calculation Processes of Wood Drying]. Moscow, 1973. 248 p.
7. Chudinov B.S. Voda v drevesine [The Water in the Wood]. Novosibirsk, 1984.
270 p.
8. Biley P.V., Petrishak I.V., Sokolovs'kiy I.A., Soroka L.Ya. Teplomassoobminni protsesi derevoobrobki derevini (Pidruchnik) [Heat and Mass Transfer Processes of Wood Processing]. Lviv, 2013. 376 p.
9. Brunner-Hildebrand. Die Schnittholztrocknung. Hannover, 1987. 322 p.
Received on February 08, 2015