К построению теории групповой окорки древесного сырья А. С. Васильев
Петрозаводский государственный университет
Анализ результатов исследований профессоров С. П. Бойкова, С. Б. Васильева, А. В. Житкова, Б. Г. Залегаллера, В. В. Коробова, В. Д. Никишова, Л. М. Перелыгина, Г. И. Торговникова, Б. Н. Уголева, М. Н. Симонова, Н. Л. Леонтьева, И. Р. Шегельмана позволил сделать вывод о недостаточной изученности особенностей процессов, происходящих в окорочных барабанах при групповой обработке древесины. В связи с этим была поставлена задача сформировать научный задел для построения теории групповой обработки древесины
[1], [2], [3], [4], [5], [6], [7] с целью отделения от нее коры.
Процесс функционирования окорочного барабана включает: загрузку в барабан древесного сырья, которое в его поперечном сечении образует сегмент заполнения; поворот сегмента до критического угла наклона его поверхности - угла динамического откоса, при котором древесина достигает наивысшей точки подъема и обрушивается вниз; обрушение древесины вниз и удар поленьев о стенку барабана или о расположенные на его внутренней поверхности рабочие органы; отделение от древесины коры и примесей, выгрузку древесины из барабана [4]. Таким образом, для решения поставленной задачи следует детально изучить движение древесины в окорочном барабане при его вращении, исследовать ударные процессы, происходящие в окорочном барабане, и сам процесс отделения коры от древесины.
Возникающие в результате взаимодействия окариваемых лесоматериалов ударные силы характеризуются кратковременностью интервала времени их действия и значительной величиной. Поэтому координаты перемещающихся отрезков древесины за время удара не изменяются, а скорости получают конечные приращения. Удобной мерой ударной силы (р) служит ее интегральная характеристика - ударный импульс (I), который определим из соотношения:
_ гп0+^_
I = | Р(г )ёг,
гп0
где , ¿0 +т] _ промежуток времени, в течение которого действует сила.
Величина ударной силы будет определять ускорения, с которыми движутся взаимодействующие бревна, а задача определения координат и скоростей сведется к интегрированию дифференциальных уравнений.
Рассмотрим движение отрезка древесины массой ( т), который движется под действием силы
(2), зависящей от радиуса вектора перемещения (г) и скорости:
V =Р = — ёг ’
где г - время перемещения, с.
В соответствии со вторым законом Ньютона:
т=2 г,г) ^
[ г = V
Каждое векторное равенство (1) эквивалентно системе трех скалярных, представляющих собой систему обыкновенных дифференциальных уравнений шестого порядка, решение которого задает параметры перемещения отрезков древесины. В дополнение к конечной силе (2) на отрезок древесины действует ударная сила (р ). На интервале времени (г0, ^ + т) запишем систему (1 ) в интегральной форме:
Г- - 1 - ч> +—
т -V(/0 + т)-V(/0)] = I + |д(г)йГ
Г0 . (2)
Го +т_ 4 '
Г (го + т)— Г (го )= ¡V (г)^Г
Го
Переходя в системе выражений (2) к пределу при т ^ 0, получим разностные уравнения импульсного движения в виде
mДV = т(^ — V2) = I, (3)
где V1 и VI - скорости в начальный и конечный момент действия импульса.
Выражение (3) целесообразно обобщать для случая, когда на бревна действует сразу несколько ударных сил р1,Р2,•••,Рь), такой эффект наблюдается при групповой обработке. Исходя из закона независимости действия сил в механике:
Ь
1 = 1^ .
1=1
Рассмотрим систему из нескольких взаимодействующих отрезков древесины с массами ту (у = 1,2, ...,к). На каждый из них действует некоторая внешняя сила (д), а так же внутренние
силы (Яу), отражающие попарные взаимодействия. В соответствии с третьим законом Ньютона
будут справедливы соотношения
Щ = —Яу.
Рассмотрим импульсное движение системы взаимодействующих отрезков древесины с учетом того, что Ру - ударная сила, действующая на отрезок древесины с номером у. По
аналогии с (1) уравнение движения у отрезка можно представить в виде:
к
т/у = + Ру +Е
1=1
(4)
Изменяя индекс у в пределах от 1 до к получим систему дифференциальных уравнений, описывающих взаимодействие отрезков древесины. Для получения уравнений импульсного движения достаточно проинтегрировать (4) на интервале действия ударных сил. Пренебрегая в полученных соотношениях членами, исчезающими вместе с т ^ 0, запишем эти уравнения в виде:
к
ту ДV] = 1у +Е 1 А 1=1
Г0 +т ,
1}г = ¡ ЯР (Г
Г0
где 1р - ударные импульсы, обусловленные взаимодействием отрезков древесины друг с другом.
Физическая сущность процесса окорки ударным способом заключается в том, что при ударе происходит ослабление силы сцепления коры с древесиной и с каждым ударом эта связь коры и древесины ослабевает, ослабление связи происходит до тех пор, пока не произойдет окончательная ее утрата, и вследствие чего, произойдет отделение участка коры.
Для рассмотрения сущности процесса отделения коры от древесины рассмотрим два типа взаимодействия бревен: при положительной температуре - ''мягкое'' взаимодействие; при отрицательной температуре - "жесткое", ''мягкое'' - допускает относительное перемещение тел, что сопровождается плавным изменением внутренних сил Я у , которые могут быть отнесены к
конечным; ''жесткое'' - характеризуется возникновением значительных усилий, которые приводят к деформации. Им отвечают ударные импульсы Iу . В качестве модели ''мягкого''
взаимодействия рассмотрим систему двух тел соединенных упругими волокнами пренебрежимо малой массы (рис. 1, а)
а) б)
Рис. 1. Расчетная схема
Внутренняя сила (r12), возникающая в системе направлена вдоль прямой соединяющей тела, а ее величина является линейной функцией расстояния между ними:
R12 = C(d + x1 - x2) ,
где C - жесткость соединительных волокон; d - длина волокон в ненапряженном состоянии; X и х - координаты коры и непосредственно древесины.
Допустим, что в начальный момент времени система находится в равновесии, а затем к телу массой m1 прикладывается ударная сила F, равная у в интервале времени t е (о, т) и нулю вне этого интервала.
В этом случае уравнения движения тел в соответствии со вторым законом Ньютона будут иметь вид
I т
42Х2 = Ñ(d + Х1 - Х2 )
В силу допущений начальные значения переменных задаются равенствами Х1 (о) = 0; X2 (о) = d; X1 (0) = 0; X2 (о) = 0 . (6)
Система выражений (5), (6) линейна и допускает явное аналитическое решение, которое в интервале (0; т) имеет вид:
m^1 — Ñ(d + Х1 — Х2 ) (5)
(
X1 =
т(т1 + m2 I • d
1
œ
m2
m1œ2
cos œt
1 2 1 1
—,--------r, — t----— +—- cos œt
т(т1 + m2 )V 2 —2 —2
где
+ т2 Л2 и2 га2
2 = ^{т\ + т“) т1 • т“
К моменту окончания удара г = т и координаты получат приращение порядка т, а для скоростей будут справедливы следующие выражения
Х1(т) = -
I
(
Х2 (т) =
m1 + m2 I
m
m 1 +
v m1
sin œт
œт
1 -
sin ют
(7)
т1 + т“ V гат
Выражения (7) есть решение прямой задачи - воздействие удара на систему тел. Анализируя его можно прийти к выводу, что конечный результат определяется не малостью промежуточного времени т, а зависит от величины произведения ит . В предельном случае при гат ^ 0, который соответствует ''мягкому'' взаимодействию древесины при положительной температуре, имеем
)
Хо =
х1(т)= —
т-1 •
Х2 (т)= 0
Скорости коры и древесины различны, то есть происходит процесс отделения коры от древесины - окорка.
Другой случай ют соответствующий "жесткому" взаимодействию мороженой древесины (см. рис. 1, б)
х1(т) = х2 (т) = —-— • (8)
т*1 + т2
Выражение (8) показывает, что скорости коры и древесины одинаковы, то есть отделение коры от древесины не наступило. Физический смысл выражения (4) - "жесткое"
взаимодействие при помощи жесткой связи длиной ё накладывает ограничения на скорости коры и древесины, что и затрудняет окорку мороженой и засохшей древесины.
При исследовании процесса перемещения древесины в окорочном барабане было установлено, что она перемещается послойного по криволинейным траекториям, подобным внутренней поверхности барабана (рис. 2). Внешний слой древесины перемещается по траектории, задаваемой внутренней поверхностью барабана и поверхностью обрушения древесины. Все остальные слои движутся по подобным траекториям с меньшими радиусами кривизны. Хаотическое перемещение бревен из одного слоя в другой не нарушает общего баланса древесины в слоях, благодаря чему суммарное количество движения слоя и его энергия остаются постоянными. Установлено, что определяющее влияние на процесс смещения одного слоя относительно другого оказывает вес вышерасположенной древесины, если ее проекция на направление перемещения слоя древесины отлична от нуля. Замедление перемещения одного слоя древесины относительно другого характеризуется коэффициентом смещения (к х+1), величина которого определяется скоростями соседних слоев древесины
к и-и+1
1+1 ,
где и и и+1 - скорости перемещения соседних слоев древесины.
Рис. 2. Расчетная схема послойного перемещения древесины
Так как слой состоит из определенного количества отрезков древесины, то скорость слоя можно рассматривать как среднюю арифметическую скоростей составляющих его отрезков древесины:
(°и
+ и2,І + ••• + и 2,і І+І \ + и3 \ і + ••• + \
z + 2
б
и =
(и у+1 + и2,г+1 + - + ^+1)+2\,/+1 + и3\,г+1 + - + \\,г+1) иг+1 =---------------------
1)1 1)о и /- ио' и \
1,г 2,г 2,г — скпппсти блевен 1-го слоя плавее пелвого отлезка' 2 >1 3 >1 2 ,
где ^, 2>г', ... , г’г' — сколости бревен 1-го слоя правее первого отрезка; скорости бревен 1-го слоя левее первого отрезка; и1’г'+1, Ь2,‘+1, ... , Ьк,‘+1 — скорости бревен 1+1-
го слоя правее первого отрезка; 2’г'+1, 3',‘+1, ... , ™'',‘+1 — скорости бревен 1+1-го слоя левее
первого отрезка; 2 — число бревен в 1-ом слое правее первого отрезка; 2' — число бревен в 1-ом
слое левее первого отрезка; ™ — число бревен в 1+1-ом слое правее первого отрезка; ™ — число бревен в 1+1-ом слое левее первого отрезка.
Так как скорость отрезков древесины первого слоя известна и равна скорости внутренней поверхности барабана, то расчет скоростей целесообразно начать со второго слоя.
и2,1 = иа — а1—2,1 ' ^1—2,1;
и3,1 = и2,1 — а2—3,1 ' ?2—3,1; и7,1 = и(;—1) —Л — а(] —1) —Л • ^—1) —Л,
где /1—21, /2—31, ... , ^—1—л — время перемещения бревна из положения 1 в 2, из 2 в 3, ..., из j-1
в j; ] - количество бревен во внешнем слое справа от первого отрезка; °1—2,1, °2—3,1, ... , а(—1)—7,1 — ускорения бревен между положениями 1-2, 2-3, ..., (]-1)-].
Найденные значения скоростей слоев древесины позволили рассчитать коэффициенты смещения соседних слоев древесины и смещения первого слоя относительно внутренней поверхности барабана. Было установлено, что они находятся в интервале 0,85-0,95. Эффективность окорки древесины, зависящая от смещения слоев древесины, будет определяться величиной работы сил трения, которая прямо пропорциональна величине этого смещения.
Изучение процесса перемещения отрезков древесины в окорочном барабане, их ударного взаимодействия и физической сущности процесса отделения коры от древесины, позволили построить математическую модель, позволяющую определять оптимальное время окорки в зависимости от требуемого качества, породы древесины, размеров обрабатываемых бревен, параметров и режима работы окорочного оборудования, которая имеет вид:
1п
Г Ь—1 агЛ
100 • 2 — г=0 г!
100 — К
п • Д2 • 1 • Кс • 0,85 • К, •К^ •-ва с 100 2•л
Т = ---------------------------, (9)
4 •а • Ща • ^уё. о- • Кид •Д
где Ь - число ударов, необходимое для полной окорки элементарного участка, расположенного на поверхности окариваемого бревна; а - среднее количество ударов, получаемых одним элементарным участком за все время обработки; К - требуемое качество окорки, %; Да -внутренний диаметр барабана, м; L - средняя длина обрабатываемых бревен, м; Кс -коэффициент заполнения барабана; К, - коэффициент полнодревесности укладки бревен в барабане; Кё - доля коры в общей биомассе ствола дерева, %; р - угол динамического обрушения, рад; ща - частота вращения барабана, мин-1; Буё — площадь элементарного участка, м2 г; кш — коэффициент, характеризующий породу древесины; Д — средний диаметр
обрабатываемых бревен, м.
Численные значения величин а, Ь, 8уё0- были установлены в результате экспериментальных
исследований [6].
Результаты исследований, представленные в работах [ФТА, патенты] показали, что одним из направлений совершенствования окорочного оборудования является выбор оптимального режима его работы. Для подбора оптимального режима работы окорочного оборудования с учетом различных факторов, касающихся как параметров обрабатываемого сырья, так и параметров окорочного оборудования может быть использована математическая модель (9) которая позволяет определить оптимальное время окорки в зависимости от требуемого качества, породы древесины, геометрических размеров обрабатываемых бревен, параметров и режима работы окорочного оборудования.
Работа выполняется при финансовой поддержке Программы стратегического развития ПетрГУ в рамках реализации комплекса мероприятий по развитию научно-исследовательской деятельности.
Литература
1. Шегельман И. Р. Анализ процесса групповой окорки при положительной и отрицательной температурах / И.Р. Шегельман, А.С. Васильев, А.Ю. Лапатин // Известия высших учебных заведений. Лесной журнал. - 2012. - Вып. 2. С.65-69.
2. Шегельман И. Р. Исследование направлений модернизации техники и технологии лесозаготовок / И. Р. Шегельман // Инженерный вестник Дона [Электронный журнал]. - 2012. -№ 2. URL: http://www.ivdon.ru/magazine/latest/n2y2012/866/
3. Шегельман И. Р. Моделирование технологического процесса очистки древесины в корообдирочном барабане с применением метода дискретных элементов / Шегельман И.Р., Васильев А.С., Колесников Г.Н. - 2008. - № 184. - С. 257.
4. Шегельман И.Р. Функционально-технологический анализ: Методология и приложения. - М: ИПиИ, 2000. - 96 с.
5. Васильев А.С. Технические решения, защищенные правоохранными документами Российской Федерации в отношении оборудования для окорки лесоматериалов / Инженерный вестник Дона. Вып. 2. Ростов-на-Дону, 2012. - [Электронный ресурс] - Режим доступа: http://ivdon.ru/magazine/latest/n2y2012/831/
6. Васильев А.С. Обоснование технических решений, повышающих эффективность режимов групповой окорки древесного сырья: автореф. дис. ... канд. техн. наук : 05.21.01: защищена 18.12.2004 - Петрозаводск: ПетрГУ, 2004. - 19 с.
7. Васильев А. С. Пути повышения эффективности процесса окорки лесоматериалов / Повышение эффективности лесного комплекса Республики Карелия: материалы третьей республиканской научно-практической конференции молодых ученых, аспирантов, докторантов. - Петрозаводск: Изд-во ПетрГУ, 2012. С. 5-7.