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

УДК 674.093

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

А.И. АГАПОВ, проф., Вятский государственный университет, д-р техн. наук

kaf_mtd@vyatsu. ги (1) Вятский Государственный Университет, 610046, Киров, ул. Московская, дом 83, кв. 94

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

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

Постановка задачи

Брусово-развальный способ раскроя пиловочника широко используется в промышленности, так как позволяет за два прохода сформировать обрезные доски. При раскрое пиловочника средних и больших размеров (диаметр бревна в вершинном торце 24-60 см) часто используется схема раскроя с выпиливанием одного бруса и трех пар боковых досок (рисунок) [1]. В теории максимальных поставов обычно определяются оптимальные размеры бруса и досок без учета ширины

Рисунок. Схема раскроя пиловочника с выпиливанием одного бруса и трех пар боковых досок с учетом ширины пропила Figure. The scheme of plank log cut with sawing a beam and three pairs of sideboards, taking into account the best opening face

пропила [1]. Затем при расчете поставов учитываются и ширина пропила и припуск на усушку. Такой подход упрощает постановку и решение задачи оптимизации раскроя пиловочника, хотя и вносит определенную погрешность в определение оптимальных размеров пилопродукции. Безусловно, ширина пропила повлияет на оптимальные размеры пиломатериалов, но математическая модель при учете ширины пропила усложняется. Следует отметить, что в теории максимальных поставов ранее не учитывался раскрой бруса на обрезные доски (второй проход). В связи с этим была поставлена задача оптимизации с учетом ширины пропила и количества пропилов при раскрое бруса на обрезные доски. В таком варианте задача оптимизации раскроя пиловочника поставлена впервые.

Для такой схемы раскроя пиловочника пока не определены оптимальные размеры бруса и досок.

Составление математической модели

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

7 = ИЛ + 2Г1-Ь1 + 2Т2Ь2 +2Г3-Ь3, (1) где Н - толщина бруса;

А - ширина пласти бруса;

Т1, Т2, Т3 - толщины первой, второй и третьей пары боковых досок соответственно;

Ь Ь Ь3 - ширины наружной пласти первой, второй и третьей пары боковых досок соответственно. Такое представление целевой функции отвечает физической сущности теории максимальных поставов при раскрое пиловочника. В формуле (1) одновременно суммарно учитываются площади поперечных сечений бруса и боковых досок. Можно предположить, что с увеличением объема бруса объем досок уменьшается и наоборот. Очевидно, имеется такое соотношение размеров бруса и досок, при котором объем этой пилопродукции получается максимальным. Такой системный подход отвечает также требованиям постановки и решения оптимизационных задач.

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

сС2 - Н2 - Л2 = 0,

где С - диаметр пиловочника в вершинном торце.

Уравнение связи для первой пары боковых досок

С - Ь12 - (Н + 2Т1 + 2е)2 = 0,

сС2 - Ь12 - Н2 - 4Т12 - 4е2 -- 4ЯТ1 - 4Не - 8Т1е = 0.

Уравнение связи для второй пары боковых досок

С2 - ь22 - (Н + 2Т1 + 2Т2 + 4е)2 = 0, С2 - Ь22 - Н2 - 4Т12 - 4Т22- 16е2 - 4НТ1 -- 4НТ2 - 8Не - 8Т1Т2 - 16Т1е - 16Т2е = 0. (4) Уравнение связи для третьей пары боковых досок

С2 - Ь32 - (Н + 2Т1 + 2Т2 + 2Т3 + 6е)2 = 0, сС2 - Ь32 - Н2 - 4Т^ - 4Т22 - 4Т32 - 36е2 -- 4НТ1 - 4НТ2 - 4НТ3 -2 8Т1Т2 - 8Т1Т3 -- 8Т2Т3 - 12Не - 24Т1е - 24Т2е - 24Т3е = 0. (5) Полагаем, что математическая модель составлена.

решение математической модели

Для решения задачи воспользуемся методом множителей Лагранжа [4]. Функция Лагранжа будет иметь вид

Ь = НА + 2Т& + 2Т2Ъ2 + 2ТъЪъ + Х(с12 -

-Н2 - А2)+ Х1(с12 -Ъ2 -Н2 -Ае2--АНТХ - АТ2 - АНе - Ще) +Х2 (,й2 - Ъ22 --И2 - АТ2 - АТ2 -1 бе2 - АНТ, - АНТ2 --87\Т2 -Ше-ХЩе-Х6Т2е) +Х3(с/2 -Ь2 --Я2 - АТ2 - АТ2 - АТ2 - 36е2 - АНТХ -

-АНТ2 - АНТЪ - 87\Т2 - 87\ТЪ --8Т2Т3 -12Не-2АТ1е-2АТ2е-2АТ3е), где X, Х1, Х2, Х3 - множители Лагранжа.

Находим частные производные от функции Лагранжа и приравниванием их к нулю

— = Н-2ХА = 0,

ал

— = А - 2ХН - 2Х,Н - АХХ - А\е -

ан 1111

-2Х2Н - АХ2Тх - АХ2Т2 - 8Х2е --2Х3Н - АХ3Тх - АХ3Т2 - АХ3Т3 --\2Х3е = 0,

— = 2Т-2ХЪ, =0,

а\ 1115

— = 2Ь.~ АХ,Н - ЯХТ, - «к.е- АХгН

атх 1 1 1112

-8А,27; -8Х2Т2 АбХ^-ЪХ^ --АХ3Н-ЪХ3Т2-8Х3Т3-2АХ3е = 0,

Л = 2Т,-2Х2Ь2=0,

аь2 аь

аъ

= 2Ъ2 - 8Х2Т2 - АХ2Н - 8Х2Тх -16Х2е -

-8Х3Т2-АХ3Н -8Х37;-8?1зГз-24Х3е = 0, - = 2Т3-2Х3Ь3 = 0,

аъ3 3 3 3 '

^ = 26, - %Х3Т3 - АХ^Н -

аъ

-%Х3Тх-ЪХ3Т2-2АХ3е = 0.

Решаем полученную систему уравнений (7) совместно с уравнением связи. Из предпоследнего уравнения системы (7) находим

Т3 = ^3, Х3 = Т3 / ьъ.

Рассматриваем последнее уравнение системы (7), которое представляем в следующем виде

Ь3 = 4Х3Т3 + 2Х3Н + 4Х3Т1 + 4Х3Т2 + 12Х3е.

В последнее равенство (9) подставим выражение (8), получим

Ь32 = 2Т3

Полученное равенство (10) подставляем в уравнение связи (5), получим

8Т32 + 6

Решая полученное квадратное уравнение (11), получим

Тъ =±(л18с12+(Н + 2Т1+2Т2+6е)2 (12) -Ъ{Н + 2Тх+2Т2 + 6е)).

Рассматриваем пятое уравнение системы (7), из которого находим

Т2 = X = Т2 / Ь2.

Рассматриваем шестое уравнение системы

В последнее равенство (14) подставим выражение (13) и получим

Ь22 = 4Т22 +2НТ2 + 4Т1Т2 + 8Т2е + Ь2Ь3.

Из последнего равенства (15) можно выразить ширину второй пары боковых досок 2Т

Ъъ = Ъ2--ЦН + 2^+ 2Тг + 4е).

Ъ2

Используя уравнение связи (5) и зная ширину третьей пары боковых досок, можно определить толщину третьей пары боковых досок

Т3=^ф2-Ь2-

Рассматриваем третье уравнение системы (7), из которого находим

Т1 = ХД, X = Т / Ь1.

Рассматриваем четвертое уравнение системы

В последнее равенство (19) подставляем выражение (18) , получим

Ь12 = 2НТ1 + 4Т12 + 4Т1е + Ь1Ь2.

Из последнего равенства (20) определяем ширину наружной пласти второй пары боковых досок

2 Т

К = + 2Т1+Ае).

Ъ1

Зная ширину второй пары боковых досок, толщину этих досок можно определить из уравнения связи

Т2=\ф2-Ь22-

Рассматриваем первое уравнение сис-темы(7), которое представляем в виде

Н = 2ХА, X = У / 2Л.

Рассматриваем второе уравнение сис-темы(7), которое с учетом четвертого уравнения этой системы можно записать как

Л = 2ХЯ + Ь1.

В последнее равенство(24) подставляем выражение (23), получим

Л2 = Н2 + АЬ1.

Из уравнения связи

Используя равенство (26), выражение (25) можно представить в виде

с!2 - 2Н2 = АЬ1.

Из последнего равенства (27) можно определить ширину первой пары боковых досок

ъ^А2-Н2 ^2-2Н2

1 Л Л '

Зная ширину первой пары боковых

досок, из уравнения связи (3) можно определить толщину этих досок

Т^ф2-Ь2-

Если в формулу (29) подставить уравнение (28), а затем решить его, то получим

1 2

Н

3d2 -АН2

d2 -Н2

~

Таким образом, рассмотрены все уравнения системы (7) и учтены все уравнения связи. Получены формулы для определения оптимальных размеров бруса и досок, но по которым непосредственно невозможно рассчитать оптимальные размеры бруса.

Алгоритм решения задачи

Для решения задачи воспользуемся численным методом. Задаемся толщиной бруса, а остальные размеры и величину целевой функции определяем по формулам. С целью упрощения расчетов и анализа результатов все полученные ранее формулы представляем в относительных единицах, полагая mH = H/d. Тогда алгоритм решения задачи будет иметь следующий вид [5].

Относительная ширина пласти бруса

а

н •

Относительная ширина первой пары боковых досок

6, _тА-тн _ 1 -2/и <1 у/1 -т2н

та

тТ

Относительная толщина первой пары боковых досок

= 2е)).

Относительная ширина наружной пласти второй пары боковых досок

тъ2=;=тъ1--к

а тъ, 1

Относительная толщина второй пары

боковых досок

Т 1 I-

«г, =7 = 2 Ч1" "++

Относительная ширина наружной

пласти третьей пары боковых досок 2шт

шът ~7 = тЧ---

м

Относительная толщина третьей пары боковых досок

Т 1 /-

'= ~с1 = 2" ^ " ^+2^+2^+66)).

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

7д! = • тТ1.

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

7д2 = 2тЬ2 * тТ2.

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

7д3 = 2ть3 * тТ3.

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

7 = 7, + 72 + 73.

д д1 д2 дЗ V '

Относительная площадь поперечного сечения бруса

7бр = тн • тл.

Суммарная относительная площадь поперечных сечений бруса и досок

т.

7 = 7 +7 .

бр д

При дальнейшем раскрое бруса (второй проход) получаются обрезные доски, а также образуются потери в виде опилок. Объем опилок зависит от числа выпиливаемых досок из бруса. Поэтому количество пропилов может быть различным. Если боковые пропилы находятся (расположены) вне плас-

ти бруса, то количество учитываемых пропилов следует определять по формуле

= п - 1,

где п - число выпиливаемых обрезных досок из бруса.

Если один из пропилов расположен в зоне пласти бруса, то количество учитываемых пропилов необходимо определять по выражению

12 = п.

Если два боковых пропила расположены в зоне пласти бруса, то количество учитываемых пропилов необходимо определять по выражению

13 = п + 1.

Число досок, получаемых из бруса, может находиться, как правило, в пределах 3.. .8 в зависимости от диаметра пиловочника и толщины получаемых обрезных досок из бруса. Тогда потери древесины в опилки (площадь поперечного сечения пропилов) можно определить по формуле

7 = I • т„ • т .

пр Н е ^ '

С учетом потерь древесины в опилки при распиловке бруса площадь поперечного сечения получаемых обрезных досок из него определится по формуле

7 = 2. - 7 .

д бр бр пр

С учетом потерь древесины в опилки, получаемых при раскрое бруса, объем обрезных досок определится по формуле

7 = 7 + 7 „ .

д д бр

результаты расчетов

В расчетах принимаем следующие варианты: I = 0, I = 2, I = 4, I = 6, I = 8.

Расчеты производятся следующим образом. Прежде всего, принимаем число пропилов i и для этого варианта выполняются расчеты. Вначале задаемся относительной шириной пропила, а затем относительной толщиной бруса в предполагаемом диапазоне 0,2.0,5 с градацией 0,05 и определяем все размеры бруса и досок, а также величину целевой функции. В результате расчета находим диапазон тн, в котором целевая функция принимает максимальное значение. Далее для этого диапазона производим расчеты с градацией тН = 0,01 и определяем все разме-

Таблица

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

пропилов при разной относительной ширине пропила The optimal relative sizes of beams and boards of various amounts of cuts at different relative kerf width

I m e тн mA mbl mT ть2 mT mb3 mT Z б д бр Z д Z

0 0 0,359 0,933 0,795 0,123 0,606 0,094 0,359 0,069 0,335 0,360 0,69586

0,005 0,372 0,928 0,779 0,122 0,582 0,088 0,336 0,059 0,343 0,333 0,67671

0 0,01 0,386 0,922 0,761 0,121 0,554 0,081 0,307 0,049 0,352 0,306 0,65813

0,015 0,4 0,916 0,742 0,120 0,524 0,075 0,279 0,039 0,360 0,279 0,64019

0,02 0,415 0,909 0,720 0,119 0,491 0,068 0,246 0,028 0,369 0,254 0,62298

0,025 0,434 0,900 0,692 0,119 0,443 0,062 0,192 0,017 0,380 0,226 0,60667

0,03 0,456 0,889 0,656 0,119 0,382 0,054 0,117 0,004 0,392 0,199 0,59161

2 0 0,359 0,933 0,795 0,123 0,606 0,094 0,359 0,069 0,335 0,360 0,69586

0,005 0,372 0,928 0,779 0,122 0,582 0,088 0,336 0,059 0,341 0,333 0,67485

0,01 0,385 0,922 0,762 0,121 0,556 0,081 0,312 0,049 0,347 0,306 0,65428

0,015 0,399 0,916 0,743 0,119 0,527 0,075 0,284 0,039 0,353 0,280 0,63420

0,02 0,414 0,910 0,721 0,118 0,493 0,068 0,251 0,029 0,360 0,254 0,61469

0,025 0,433 0,901 0,693 0,118 0,446 0,062 0,197 0,017 0,368 0,227 0,59585

0,03 0,45 0,893 0,666 0,117 0,402 0,054 0,152 0,006 0,374 0,203 0,57802

4 0 0,359 0,933 0,795 0,123 0,606 0,094 0,359 0,069 0,335 0,360 0,69586

0,005 0,371 0,928 0,780 0,122 0,584 0,087 0,340 0,059 0,337 0,334 0,67114

0,01 0,383 0,923 0,764 0,120 0,561 0,081 0,321 0,049 0,338 0,308 0,64660

0,015 0,396 0,918 0,747 0,119 0,535 0,075 0,299 0,039 0,339 0,282 0,62227

0,02 0,41 0,912 0,728 0,117 0,505 0,068 0,271 0,029 0,341 0,257 0,59821

0,025 0,425 0,905 0,705 0,116 0,471 0,061 0,239 0,019 0,342 0,232 0,57449

0,03 0,443 0,896 0,677 0,116 0,425 0,054 0,192 0,008 0,343 0,207 0,55124

6 0 0,359 0,933 0,795 0,123 0,606 0,094 0,359 0,069 0,335 0,360 0,69586

0,005 0,37 0,929 0,781 0,121 0,587 0,087 0,345 0,059 0,332 0,334 0,66743

0,01 0,382 0,924 0,766 0,120 0,564 0,081 0,326 0,049 0,330 0,308 0,63895

0,015 0,394 0,919 0,750 0,118 0,541 0,074 0,308 0,040 0,326 0,284 0,61041

0,02 0,407 0,913 0,732 0,117 0,514 0,068 0,286 0,030 0,322 0,259 0,58186

0,025 0,421 0,907 0,711 0,115 0,483 0,061 0,260 0,020 0,318 0,235 0,55333

0,03 0,436 0,899 0,689 0,114 0,447 0,054 0,229 0,009 0,313 0,211 0,52489

8 0 0,359 0,933 0,795 0,123 0,606 0,094 0,359 0,069 0,335 0,360 0,69586

0,005 0,37 0,929 0,782 0,121 0,587 0,087 0,345 0,059 0,328 0,335 0,66373

0,01 0,381 0,924 0,767 0,119 0,567 0,081 0,331 0,050 0,321 0,309 0,63132

0,015 0,392 0,919 0,752 0,118 0,546 0,074 0,317 0,040 0,313 0,285 0,59862

0,02 0,404 0,914 0,736 0,116 0,522 0,067 0,300 0,030 0,304 0,260 0,56564

0,025 0,417 0,908 0,717 0,114 0,494 0,061 0,279 0,020 0,295 0,237 0,53239

0,03 0,43 0,903 0,698 0,113 0,466 0,054 0,259 0,010 0,285 0,214 0,49890

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

Результаты расчетов показывают, что с увеличением ширины пропила при первом проходе раскроя бревна толщина бруса возраста-

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

Анализируя результаты расчетов, следует отметить, что при раскрое двухкантного бруса с числом пропилов г = 4...6 и относительной ширине пропила т = 0,01...0,025 оп-

тимальная относительная толщина этого бруса изменяется в пределах тН = 0,38...0,42 и в среднем составляет 0,4 от диаметра бревна в вершинном торце.

Выводы и рекомендации

Решена задача оптимизации раскроя пиловочника средних размеров с получением бруса и трех пар боковых досок с учетом ширины пропила и количества пропилов при раскрое бруса на обрезные доски. Численным методом получены оптимальные относительные размеры бруса и досок в зависимости от ширины пропила, которые рекомендуется использовать при выборе варианта постава. Подтверждается ранее принятая гипотеза о том, что с увеличением объема бруса объем боковых досок уменьшается и наоборот. Таким образом, имеется такое соотношение размеров бруса и боковых досок, при котором объем этой пилопродукции получается максимальным. С увеличением ширины пропила толщина бруса возрастает, а размеры боковых досок уменьшаются. При раскрое бруса с увеличением количества пропилов оптимальная толщина бруса незначительно уменьшается,

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

Библиографический список

1. Аксенов, П.П. Теоретические основы раскроя пиловочного сырья / П.П. Аксенов. - М: Лесная пром-сть, 1960. - 216 с.

2. Агапов, А.И. Алгоритм определения оптимальных размеров брусьев и досок при раскрое пиловочника бру-сово-развальным способом / А.И. Агапов // Механика технологических процессов в лесном комплексе: между -нар. научно-практическая конференция, Воронеж, 25-27 марта 2014 г. - Воронеж: ВГЛТА, 2014. - С. 287-291.

3. Агапов, А.И. Оптимизация раскроя пиловочника средних размеров при брусово-развальном способе распиловки / А.И. Агапов // Современная наука : актуальные проблемы и пути их решения : 7-я междунар. дистанционная науч. конф., Липецк, 20-21 февраля 2014 г - Липецк : МАКСИМАЛ, 2014. - С. 16-24.

4. Агапов, А.И. Оптимизация раскроя пиловочника с выпиливанием трех брусьев разной толщины и двух пар боковых досок / А.И. Агапов // Актуальные проблемы развития лесного комплекса: междунар. научно-технич. конф., Вологда, 3-4 декабря 2013 г. - Вологда: ВОГУ, 2014. - С. 62-66.

5. Агапов, А.И. Влияние ширины пропила на оптимальные размеры брусьев и досок при раскрое пиловочника с выпиливанием трех брусьев одинаковой толщины и четырех пар боковых досок / А.И. Агапов // Лесотехнический журнал. - 2014. - Т.4, № 2. - С. 128-135.

OPTIMIZATION ALGORITHM OF SOLVING THE TASK OF PLANK MEDIUM-SIZED AND LARGE-SIZED LOG CUT BY LOG SQUATTER BREAK-UP METHOD Agapov A.I., Prof. Vyatka State University, Dr. Sci.

kaf_mtd@vyatsu.ru

(1) Vyatka State University, 610046, Kirov, st. Moscow, 83, Apt. 94 Summary: the mathematical model of plank log cut by cutting out a bar and three pairs of lateral boards taking into account the width has been developed. The mathematical model includes the target function and four equations of link. To solve the mathematical model the method of Lagrange multipliers has been used. The determination of the optimal size of a beam and side boards was carried out by the numerical method. In this regard, the algorithm of solving this problem in relative units has been developed. With increasing width of cut, the beam width increases, and the size of lateral boards decreases. When cutting plank logs with increasing width of the cuts, the beam width is slightly reduced, and the size of lateral boards increases.

Keywords: cutting of logs, a beam and boards, the optimality criterion, mathematical model, the objective function, equations of link, Lagrange method, numerical method, an algorithm of the problem.

References

1. Aksenov P.P. Teoreticheskie osnovy raskroya pilovochnogo syr 'ya [Theoretical Foundations of cutting sawlogs]. Moscow: Forest Industry, 1960. 216 p.

2. Agapov A.I. Algoritm opredeleniya optimal 'nykh razmerov brus 'ev i dosok pri raskroepilovochnika brusovo-razval 'nym sposobom [The algorithm for determining the optimal size of beams and planks when cutting logs by log squatter method]. Mechanics of technological processes in the forest complex: Intern. Scientific Conference, Voronezh, 25-27 March 2014. Voronezh: VGLTA, 2014. pp. 287-291.

3. Agapov A.I. Optimizatsiya raskroya pilovochnika srednikh razmerov pri brusovo-razval'nom sposobe raspilovki [Optimization of cutting logs with a medium-sized lumber sawing break-up method]. Modern science: current problems and their solutions: 7th Intern. Remote scientific. Conf., Lipetsk, 20-21 February 2014. Lipetsk: MAXIMUM, 2014. pp. 16-24.

4. Agapov A.I. Optimizatsiya raskroya pilovochnika s vypilivaniem trekh brus'ev raznoy tolshchiny i dvukh par bokovykh dosok [Optimization of cutting sawing logs with three beams of different thicknesses and two pairs of side boards]. Actual problems of forest sector development: Intern. Scientific-Technical. Conf., Vologda, December 3-4, 2013. Vologda: VOGU, 2014. pp. 62-66.

5. Agapov A.I. Vliyanie shiriny propila na optimal'nye razmery brus'ev i dosok pri raskroe pilovochnika s vypilivaniem trekh brus'ev odinakovoy tolshchiny i chetyrekh par bokovykh dosok [The effect of cutting width on the optimum size of beams and boards when cutting logs into three beams of the same thickness and four pairs of sideboards]. Vyatka State University, FAM, Dept. of Forestry MTD. 2014. T. 4, № 2. pp. 128-135.