16.02.09. № 75

Авторефераты по всем темам >> Авторефераты по разное

На правах рукописи

ГАЗИЗОВ Асгат Мазхатович

СОВЕРШЕНСТВОВАНИЕ ТЕХНОЛОГИИ РОТОРНОЙ ОКОРКИ ЛЕСОМАТЕРИАЛОВ ПУТЕМ ОПТИМИЗАЦИИ ОСНОВНЫХ ПАРАМЕТРОВ ПРОЦЕССА

05.21.01. - Технология и машины лесозаготовок и лесного хозяйства

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

доктора технических наук

Санкт-Петербург - 2010.

Работа выполнена на кафедре технологии лесозаготовительных производств Санкт-Петербургской государственной лесотехнической академии имени С.М. Кирова

Научный консультант - доктор технических наук, профессор

Григорьев Игорь Владиславович

Официальные оппоненты - доктор технических наук, профессор

Минаев Александр Николаевич

доктор технических наук, профессор

Шегельман Илья Романович

доктор технических наук, профессор

Герц Эдуард Федорович

Ведущая организация - Государственное образовательное

учреждение высшего

профессионального образования

Московский государственный

университет леса

Защита диссертации состоится л___ я 2010 г. в _______ на заседании диссертационного Совета Д.212.008.01 в Северном (Арктическом) федеральном университете / 163002, Архангельск, набережная Северной Двины 17/.

С диссертацией можно ознакомиться в библиотеке университета.

Автореферат разослан л_____ сентября 2010 г.

Ученый секретарь

диссертационного Совета

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. В настоящее время перед лесным сектором экономики РФ остро стоит ряд проблем, одной из которых, является повышение эффективности лесоэксплуатации и деревопереработки. Известно, что окорка древесины является одной из наиболее энерго- и трудоемкой операцией первичной лесопереработки, которая осуществляется на всех типах лесопромышленных складов.

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

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

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

Развитие отраслевой науки должно быть направлены на создание лэффективной системы использования природных ресурсов. Именно это требование содержалось в одном из посланий Президента РФ Федеральному Собранию. Оно конкретизировано и развито в Концепциях развития лесного хозяйства и лесопромышленного комплекса, одобренных Правительством России. В Перечень Приоритетных направлений развития науки, технологий и техники в РФ, утвержденного Президентом РФ 21 мая 2006 г. Пр-843 включен пункт Рациональное природопользование. Вышесказанное позволяет утверждать, что повышение эффективности окорки лесоматериалов является весьма актуальной проблемой.

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

Объект исследований. Кора и древесина основных лесорастительных пород.

Предмет исследования. Процесс механической окорки лесоматериалов в роторных окорочных станках.

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

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

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

Апробация работы. Основные положения диссертации и отдельные ее разделы докладывались и обсуждались на Научной конференции профессорско-преподавательского состава БГАУ (Уфа, 1993 г.); Международной научной конференции Лес - 2000 (Брянск, 2000 г.); Научно-практической конференции Интеграция науки, образования и производства для развития лесного хозяйства и лесопромышленного комплекса (Воронеж, 2001 г.); Всероссийской научно-практической конференции Повышение эффективности и устойчивости развития АПК (Уфа, 2005 г.); Международной научно-практической конференции Актуальные проблемы лесного комплекса (Брянск, 2008 г.); Всероссийской научно-практической конференции Интеграция аграрной науки и производства: состояние, проблемы и пути решения (Уфа, 2008 г.); Первой и второй международных научно-практических Интернет конференциях Леса России в XXI веке (Санкт-Петербург, 2009 г.); Международной научно-технической конференции Актуальные проблемы развития лесного комплекса (Вологда, 2009 г.); Межвузовской научной конференции Молодая мысль: Наука. Технологии. Инновации (Братск, 2009 г.); Международной научно-практической конференции Современные проблемы лесозаготовительных производств, производства материалов и изделий из древесины: пиломатериалы, фанера, деревянные дома, заводского изготовления, столярно-строительные изделия (Санкт-Петербург, 2009 г.); Всероссийской научно-практической конференции аспирантов, соискателей иадокторантов (МГТУ, г. Майкоп, 2008 г.); и ежегодных научно-технической конференции Санкт-Петербургской государственной лесотехнической академии имени С.М. Кирова в 20072009 гг. Часть материалов работы вошла в проект Селективная окорка лесоматериалов, который по итогам конкурса на лучшие инновационные проекты в сфере науки и высшего образования в 2009 году, проводимого Комитетом по науке и высшей школе Санкт-Петербурга, признан победителем в номинации Лучшая научно-инновационная идея. Часть материалов работы получена при выполнении НИР по государственному контракту П1209 по ФЦП Научные и научно-педагогические кадры инновационной России на 2009-2013 годы, по направлению Переработка и утилизация техногенных образований и отходов в рамках мероприятия 1.3.1.

Публикации. Материалы диссертации опубликованы в сорока печатных работах, включая одну монографию. Десять статей опубликованы в журналах, рекомендованных ВАК РФ для публикации результатов докторских диссертаций. Результаты исследований отражены в научно-технических отчетах по НИР.

Структура и объем работы. Диссертация состоит из введения, 6 глав, основных выводов и рекомендаций, списка литературы. Общий объем работы 350 страниц. Диссертационная работа содержит 93 рисунка, 26 таблиц. Список литературы содержит 319 наименований.

На защиту выносятся следующие положения:

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

СОДЕРЖАНИЕ РАБОТЫ

Во введении сформулирована актуальность темы диссертационной работы, определена цель исследований, изложены научная новизна, значимость для теории и практики, сформулированы основные положения выносимые на защиту.

1. Состояние проблемы и задачи и исследования.

В разделе произведен обзор и анализ литературных источников по повышению качества и эффективности окорки. Рассмотрено место окорки в современных технологических процессах лесозаготовительного и деревообрабатывающего производства; дан сравнительный анализ способов окорки круглых лесоматериалов; выполнен анализ типоразмерных рядов роторных окорочных станков и их инструментального оснащения. Проанализировано влияние строения и состояния коры на параметры процесса окорки. Сделан обзор авторских свидетельств и патентов по окорочно-зачистным станкам и устройствам к ним. Критически проанализированы современные методы расчета параметров процесса окорки древесины.

Фундаментальный вклад в теорию окорки лесоматериалов внесли С.П. Бойков, С.Б. Васильев, А.В. Житков, Б.Г. Залегаллер, В.Д. Никишов, Л.М. Перелыгин, Г.И. Торговников, Б.Н. Уголев, М.Н. Симонов, Н.Л. Леонтьев, И.Р. Шегельман, Н.Ф. Пигильдин, Ф.И. Копернн и др. ученые СПбГЛТА, МГУЛ, ВГЛТА, ПетрГУ, ЦНИИМЭ, КарНИИЛПКа и др.

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

На основании анализа сформулированы выводы и задачи исследования:

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

2. Методика оценки и управления эффективностью работы окорочного оборудования

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

3. Теоретические исследования процесса разрушения коры при роторной окорке древесины

При окорке круглых лесоматериалов в сложных условиях, характеризуемых факторами влажности (W,%), температуры (T,oC), большими диаметрами бревна (dб, м) и рядом других, эффективность процесса отделения коры при использовании роторных станков зависит от большого числа параметров, в ряде случаев противоположного действия. Основными из этой совокупности параметров являются угловые, силовые и кинематические параметры окорки.

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

Рассмотрим процесс отделения коры от древесины с помощью скребкового короснимателя при следующих угловых параметрах (рис. 1, а). Примем в качестве угла окорки (резания) - угол >/2 между передней гранью кулачка и плоскостью, касательной к поверхности бревна диаметром dб в месте контакта ее с рабочей кромкой, а за установочный угол резания о > /2 примем угол между передней гранью кулачка и плоскостью, проходящей через ось качания и рабочую кромку. Выбор угловых характеристик определяется: диаметром бревен dб и параметрами станка расстоянием от оси качания короснимателей до рабочей кромки (L) и диаметром оси ротора (R).

Отделение коры от древесины происходит путем реализации механизма сдвига под действием давления передней грани кулачка. Рассмотрим нормальную силу Fс, перпендикулярную к плоскости передней грани, и действующую на участок коры толщиной hк с площадью контакта s (рис. 1, б), величина которой зависит от геометрических и угловых параметров рабочей кромки короснимателя и диаметра бревна dб. Указанная сила связана с силой прижима короснимателя F1 соотношением:

. (1).

Нормальной силе Fс препятствуют две силы (рис. 1 б): Т1 сила трения коры о коросниматель

, (2)

где тр коэффициент трения и Т2 сила внутреннего сцепления частиц коры (С):

, (3)

k=tg коэффициент внутреннего трения.

Суммируя силы Т1 и Т2 и, разделив результат на площадь контакта s, определим величину горизонтального давления qг.

Рис. 1. Схема разрушения массива коры:

а) взаимодействие короснимателя с корой; б) слой коры под давлением; в) предельные круги Мора; 1-древесина; 2- кора; 3-коросниматель; 4- окоренная поверхность

Зависимость давления qг от вертикального давления qв установим с помощью обобщенной диаграммы Мора с учетом уменьшения силы внутреннего сцепления вследствие сокращения площади контакта при увеличении сдвига частиц коры:

, МПа, (4)

где j bк сдвиг частиц коры, предельное значение которого равно величине подачи бревна за 1 оборот короснимателя; bк - ширина снимаемой коры. Необходимо отметить, что в пределе отношение , т.е. оно является обратной величиной коэффициента перекрытия Kп.

Таким образом, анализ соотношения (4) показывает, что при деформации сдвига j= bк значение сил внутреннего сцепления снижается до нуля.

Диаграмма Мора для двух кругов предельных напряжений с диаметрами соответственно равными пределу прочности на разрыв р и сж представлена на рис. 1, в. Секущая Sт, проходящая через вершины А и В окружностей, в первом приближении представляет собой огибающую кругов Мора, т.е. угол ее наклона к оси Ох примем за (k =tg), а ординату точки пересечения с осью Оу за сцепление С. В этом случае, когда наряду с сжимающими (отрицательными) напряжениями присутствуют растягивающие (положительные), величины С и следует рассматривать не как истинные физические характеристики, а как параметры аппроксимированной диаграммы сдвига, отражающие влияние сцепления и угла внутреннего трения на развитие процесса разрушения.

Тогда получим соотношения для определения параметров С и k:

. (5)

Вертикальное давление qв, действующее на переменный слой hс коры, определим из решения задачи о вдавливании штампа (рис. 1, б).

Деформация (дифференциал dhс) элементарного слоя в направлении Z действия qв равна:

, (6)

где: Е - модуль общей деформации коры при отсутствии сдвига, равный начальному давлению q0 при hс=hr, где hr величина углубления рабочей кромки короснимателя в массив коры; ; геометрический параметр короснимателя (штампа).

Наряду с абсолютными значениями hс и hк будем рассматривать их относительные величины и .

Интегрируя выражение (6) по переменной Z в пределах деформируемого слоя, получим соотношение для определения вертикального давления qв с учетом коэффициента уплотнения ky, отражающего увеличение общей деформации массива коры при погружении ядра уплотнения:

. (7)

Коэффициент уплотнения может быть учтен через несущую способность среды, через коэффициент k внутреннего трения, и через коэффициент kпр пористости до начала деформации (kпр=1+0, где 0 отношение объема пор в массиве коры к объему твердых частиц) и его изменение в процессе линейного деформирования. Последний подход применительно к деформации коры наиболее целесообразен. Поэтому принимаем, что уплотнение происходит пропорционально величине .

Установленные соотношения (4) и (7) позволяют определить (рис. 1, б) приведенное давление интегральную характеристику нагрузки в произвольной точке массива коры под совокупным сжимающим действием вертикального и горизонтального давлений.

В качестве критерия качественного отделения коры принимается следующее условие: на границе раздела кора-древесина величина приведенного давления достигает величины предельной характеристики прочности на скалывание ск, т.е. должно выполнятся условие:

. (8)

Принимая во внимание, что значение ск коры отличается от аналогичной характеристики для древесины в 210 раз, в зависимости от породы, влажности и температуры, можно заключить, что принятое условие (8) обеспечит качественную окорку на полную глубину без повреждения заболонного слоя древесины.

Для апробации разработанной модели были выбраны несколько пород древесины, отличающихся как по физико-механическим свойствам, так и по условиям деформирования в соответствии с диаграммами Мора. Для семи различных древесных пород был выполнен расчет С и k, который показал, что минимальной характеристикой сцепления обладает кора лиственницы, максимальной - осины и дуба. Кора хвойных деревьев характеризуется достаточно узким диапазоном изменения величины С=1,41,61 МПа. Диапазон изменения величины внутреннего трения достаточно широкий и изменяется на порядок от минимального значения k=0,08 для осины до максимальной величины k=0,8 для коры березы.

Учитывая различия в плотностях, прочностных характеристиках и параметрах диаграмм Мора, для сравнительного анализа были выбраны три породы древесины: сосны, осины и березы. Выбор этих пород обусловлен и тем обстоятельством, что, основываясь на известных регрессионных уравнениях зависимости толщины коры hк от dб, при заданном диаметре бревна, в частности dб=0,4м, толщина их коры практически совпадает и равна соответственно hк=12,8, 12,1 и 12,4 мм. Это позволило на данном этапе исследований исключить при проведении сравнительного анализа фактор влияния диаметра бревна.

На рис. 2 представлен характер изменения вертикального qв, горизонтального qг и приведенного давлений в массиве коры березы с ростом толщины hс снимаемого слоя (% от толщины коры hк). Приведенные данные соответствуют следующим параметрам окорки свежесрубленного бревна диаметром dб=0,4м на станке ОК-63: скорость подачи uп=0,45 м/с; F1=800 Н; число короснимателей 4; число оборотов в минуту 135, частота вращения ротора - 2,7 с-1; L= 0,366 м, его задний и передний углы - соответственно 0,785 и 0,25 рад.; =2,02 рад.; радиус заточки 0,0015 м; hr=0,0016 м; bк=0,05 м; Kп =1, kпр=1,25, тр=0,3. При таких параметрах удельная сила окорки достигает =22,17 кН/м.

Для данных условий окорки березы величина предела прочности ск коры поперек волокон по камбиальному слою при Т=0оС принимается равной 1,42, а древесины- 3,86 МПа. Таким образом, на основании анализа данных рис. 2 можно отметить, что достигнутое значение приведенного давления =0,81 МПа составляет лишь 57% от величины ск, т.е. условие (8) не выполняется и параметры окорки нуждаются в корректировке.

Рис. 2. Изменение давления в массиве коры березы:

1 изменение вертикального qв давления; 2 изменение горизонтального qг давления; 3 изменение приведенного давления

Сравнение, для аналогичных условий окорки, характера изменения приведенного давления в массиве коры для трех пород древесины: сосны, осины и березы показало, что величина ск для первых двух пород составляет, соответственно ск=0,33 и 0,79 МПа. То есть, принятые параметры окорки обеспечат выполнение условия (8) только для массива коры сосны. Причем в этом случае они нуждаются в существенной корректировке в сторону снижения величины . Для качественной окорки осины, в отличие от березы, требуется незначительное увеличение давления.

Рассмотрим три основных способа изменения величины (управления процессом разрушения массива коры) путем изменения: 1) угловых, 2) силовых, 3) кинематических параметров окорки.

Управление угловым параметром. Задаваясь диапазоном изменения угла окорки =1,742,35 рад при фиксированных значениях остальных параметров, была получена функциональная логарифмическая зависимость величины от . Которая, в частности, показала, что для условий окорки березы, при постоянной силе прижима F1=800 Н, за счет изменения углового параметра удельная сила окорки изменялась от 8,6 до 35,87 кН/м. Однако достигнутое значение =1,21 составляет лишь 86% предельной характеристики прочности ск. И только увеличение силы F1 до 1000 Н (соответственно силы до 44,84 кН/м) обеспечивает выполнение условия (8).

Полученные результаты позволили провести исследования степени влияния угла окорки на выполнение критерия (8) для различных пород. Для этого необходимо установить зависимость безразмерной функции от размерного угла , после чего производная будет являться характеристикой искомой интенсивности (степени влияния). На рис. 3 представлены графики зависимости от для сосны, осины и березы. Линейный характер зависимостей позволяет сделать вывод о том, что производные являются угловыми коэффициентами прямых и составляют соответственно 3,74, 1,07 и 3,24. Эти данные показывают, в частности, что при прочих равных условиях окорка сосны и березы характеризуется существенно большим влиянием углового фактора, чем окорка осины.

Управление силовым параметром. На рис. 4 для условий окорки березы представлен функциональный линейный характер изменения приведенного давления от силы F1. Как видно, увеличение F1 до 1600 Н обеспечивает выполнение критерия (8). Удельная сила составила =44,33 кН/м. Этот результат, при его сравнении с соответствующим показателем =44,84 предыдущего (углового) способа управления, свидетельствует о том, что параметр является характеристикой сопротивления массива коры разрушающему действию короснимателя и не зависит от способа управления процессом окорки. Результаты исследований по установлению степени влияния силы прижима F1 на процесс окорки по аналогии с влиянием угла представлены на рис. 5, где отложены значения относительного критерия , и F1, Н.

Рис. 3. Влияние угла окорки на процесс разрушения коры:

1 - сосны; 2 - осины; 3 березы

Рис. 4. Зависимость величины приведенного давления от силы прижима

Сравнивая данные рис. 5 и 3, отметим качественное совпадение влияния параметров управления, при имеющихся двух количественных отличия: во-первых, изменилась в сторону увеличения пропорция степени относительного влияния (отношение угловых коэффициентов прямых) и, во-вторых, диапазон изменения для углового параметра (=12,9) значительно шире соответствующего диапазона для силового (=11,8). Это означает, что результаты процесса разрушения коры в большей степени зависят от изменения угловых параметров, чем силовых.

Управление кинематическим параметром. Зафиксируем угол окорки и силу прижима: =1,74, F1=800Н. Переменной величиной является скорость uп подачи бревна в диапазоне uп=0,250,45м/с (Kп =12).

Рис. 5. Влияние силы прижима на процесс разрушения коры:

1 - сосны; 2 - осины; 3 березы

На рис. 6 в рамках вышеотмеченных (рис. 3 и 5) обозначений для трех пород древесины представлены , и uп. Как видим, кинематический параметр при окорке осины по степени влияния существенно превысил аналогичные результаты для сосны и березы, тем самым, компенсируя слабое влияние углового и силового параметров. Необходимо отметить уменьшение пропорции степени относительного влияния, т.е. кинематический параметр уменьшает дифференциацию в развитии процесса разрушения коры различных пород. В тоже время, диапазон изменения расширился до значения =14,5, что свидетельствует о том, что данный параметр оказывает существенное влияние на процесс разрушения коры.

Рис. 6. Влияние скорости подачи на процесс разрушения коры:

1 - сосны; 2 - осины; 3 - березы

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

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

По своей структуре, как известно, кора деревьев представляет собой многослойный материал, состоящий из корки, луба и камбия. С позиции механики сплошных сред каждый слой и кора в целом представляют собой трехкомпонентную среду, содержащую: 1) твердую (перидерма, рыхлая паренхима, каменистые клетки, волокна); 2) жидкую (вода, при низких температурах - лед); 3) газообразную (защемленный воздух) компоненты.

Известно, что кора легко отделяется при положительных температурах, а также при влажности не менее 4050%. Однако процесс образования окоренной поверхности сухого или мерзлого сырья усложняется и характеризуется проявлением слабо изученных механизмов смятия, уплотнения коры с последующим сдвигом вдоль плоскости раздела кора-древесина.

Таким образом, на процесс отделения коры, помимо соотношения деформационных и прочностных характеристик, оказывают влияние ее влажность (W,%) и температура (То,С), поскольку вода и лед обладают различной сжимаемостью. Наряду с этим агрегатное состояние влаги характеризуются проявлением отличных друг от друга упругопластических и вязких свойств, что существенно изменяет прочность адгезионных связей между твердыми компонентами коры и заболонного слоя древесины.

На данном этапе исследований была дана оценка влияния влажности на развитие процесса разрушения коры. В результате статистической обработки физико-механических свойств коры различных древесных пород установлено, что существует корреляционная связь между величинами сж и k, тогда как между р и k она статистически не значима. В результате расчетов была установлена зависимость увеличения относительной величины предела прочности на сжатие от относительного увеличения плотности по сравнению с начальным состоянием, которую с коэффициентом детерминации R2= 0,63 можно выразить в виде:

= 2,62 - 1,7151. (9)

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

Обобщив известные опытные данные для коры четырех пород - ели, сосны, березы и лиственницы в комлевой и срединной частях хлыста, установлен логарифмический закон связи () для всех пород деревьев. Расчеты для различных пород деревьев показали, что коэффициент при натуральном логарифме (для ели он равен Kw=0,6724) является функцией начальной плотности коры kо.

В итоге получена зависимость () в виде:

(10)

Зависимость (10) отличается от известной теоретической зависимости = 1+W, полученной для оценки влияния влажности на плотность трехкомпонентной среды и более полно отражает протекание этого сложного процесса.

При анализе физико-механических свойств коры различных пород выделяется широкий диапазон изменения характеристики сцепления лиственных деревьев (С=1,185,77 МПа), тогда как кора хвойных деревьев характеризуется более узким диапазоном изменения величины С=1,41,61 МПа. Разрушение коры зависит от величины сцепления С, которая, как показывает взаимосвязь соотношений (9)(10), в свою очередь, зависит от влажности W. Полученные исходные данные о влиянии W на характеристики коры позволили перейти к рассмотрению вопроса оценки влияния влажности среды на механизм развития разрушения массивов коры различных пород. Для сопоставительного анализа были выбраны кора сосны, осины и березы для условий окорки свежесрубленного (W=100-130%) бревна диаметром dб=0,4м на станке ОК-63.

Рассмотрим влияние W на развитие процесса разрушения через механизм зависимости ск(W), учитывая, что ск снижается по мере увеличения относительной влажности . Степень этого снижения была определена при обработке известных опытных данных статистическими методами распознавания образов для различных пород деревьев.

Анализ показал, что рост W сильнее сказывается на снижении прочности, чем на увеличении давления в слое коры, в связи с чем, зависимость критерия S(W) - практически линейная, положительная. Предельное состояние (S1) наблюдается при снижении W до 50%. Полученные результаты показали, что чем меньше исходная плотность коры, тем больше влияние влажности на параметры процесса разрушения.

Оценка интенсивности влияния W на S выполнялась посредством производной dS/dW, поскольку в силу линейной связи S(W) производная является угловым коэффициентом наклона прямых. Установлено, что при плотности коры сосны более, чем в 2 раза ниже плотности коры березы, степень влияния влажности на процесс ее разрушения в 4 раза превосходит этот показатель при разрушении высокоплотной коры березы по мере насыщения ее водой. Таким образом, разработанная на данном этапе математическая модель позволяет аналитически учесть влияния влажности коры различных пород на развитие процесса ее разрушения.

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

С понижением температуры внутреннее сцепление частиц С будет усиливаться. При отрицательных температурах (T<0) вода в порах замерзает, а величина С во многом зависит от изменения сил сцепления льда с твердыми компонентами коры. В результате анализа данных об изменении относительной величины внутреннего сцепления в зависимости от относительной плотности материала при различных отрицательных температурах, для любых значений температуры, получена зависимость:

= 0,33 еxp{ (-0,04Т+1,07) }. (11)

Полученные соотношения являются основой для расчета величины приведенного давления и в случае анализа процесса разрушения коры при понижении температуры Т. Для использования критерия разрушения (3) необходимо уточнить влияние Т на величину предела прочности коры на скалывание ск. Статистическая обработка известных опытных данных позволила установить корреляционные зависимости ск(Т) для широкого ряда свежесрубленных (W 100%) материалов коры различных пород. Необходимо также учитывать то обстоятельство, что величины ск зависит не только от Т, но и диаметра бревна dб, поскольку от последнего зависит толщина коры и степень ее насыщения влагой (льдом). Известные экспериментальные исследования на образцах летних и зимних бревен диаметром dб =0,080,2м для различных пород позволили установить зависимости ск(dб). Полученные результаты позволили сделать ряд выводов: во-первых, установлен линейный характер возрастания относительной прочности коры с увеличением относительного диаметра, т.е. масштабный эффект среды разрушения является весьма существенным; и, во-вторых, что проявление этого масштабного эффекта в зимнее время существенно выше, чем при окорке бревен при положительной температуре.

На основании полученных результатов стало возможным оценить влияние фактора понижения температуры среды на механизм разрушения массива коры различных пород. Апробация модели при изменении Т от 0 до -20оС выполнена для окорки сосны на станке ОК-63. Некоторые характеристики процесса принимались переменными в следующих пределах: dб=0,1-0,63м, uп=0,15-0,6 м/с, F1=7302900 Н, =1,742,35 рад.; =135-180 об/мин, задний и передний углы короснимателя - 0,44 и 0,25 рад., Kп=14.

На рис. 7 показан характер изменения величины по мере увеличения глубины снимаемого слоя коры (в % от hk) для четырех температурных состояний. Данные соответствуют следующим значениям параметров окорки: dб =0,4м; uп=0,3м/с; F1=730 Н; =1,74 рад.; =180 об/мин; Kп=2. Полученные результаты сопоставлены с предельной характеристикой прочности ск сосны для данных условий, которая по мере снижения Т от 0 до -20оС возрастает и составила соответственно ск=0,34; 0,76; 1,12 и 1,66 МПа.

Рис. 7. Изменение приведенного давления в массиве коры сосны:

1-Т=0; 2-Т=-5; 3-Т=-10; 4-Т=-20оС

Сравнительный анализ показывает, что предельное условие (8) качественной окорки на полную глубину выполняется до значения Т=-7оС (удельная сила окорки достигает величины =16,8 кН/м). Дальнейшее снижение Т требует адекватной корректировки параметров окорки в направлении увеличения давления и выполнения условия (8).

В рамках разработанной модели принимаем, что основными параметрами управления величиной являются: угловой, силовой и кинематический, т.е. путем изменения одного параметра при постоянстве остальных стремимся к выполнению условия (8). Установлено, в частности, что при снижении Т до -15оС и ниже для обеспечения качественной окорки необходимо увеличить угол до 2,05рад (в 1,18 раз) и силу прижима F1- до 2000Н (в 2,74 раза), а также снизить скорость подачи uп до 0,15м/с (в 2 раза), увеличив тем самым Kп до 4. В совокупности эти меры приводят к тому, что удельная сила окорки возрастает до =45,9 кН/м. Зависимость (Т) практически функциональная и подчиняется параболическому закону. Фактор температуры в большей степени отражается на силовом и кинематическом параметрах, чем на угловом.

Наиболее важным моментом в процессе разрушения сплошного мерзлого массива при изменении внешней среды является переход температуры Т через нулевое значение. В этой связи поставлена задача изучения процесса разрушения коры в узком диапазоне изменения Т: от 2 до -6 оС, поскольку в этом промежутке происходит фазовое превращение жидкости.

Разработанная модель расчета была реализована на примере трех пород древесины - сосны, осины и березы при следующих параметрах окорки: dб=0,4м; uп=0,45м/с; F1=7302900 Н; =1,741,9 рад.; =135 об/мин; Kп=1, т.е. при постоянном перекрытии путем управления угловым параметром в узком диапазоне и, в основном, за счет силового параметра F1 обеспечивали выполнение критерия полной окорки (8). В качестве интегральной характеристики процесса разрушения принято достигнутое значение удельной силы .

Расчет удельной силы окорки с понижением температуры вблизи нулевого значения для различных древесных пород показал, что, в наибольшей степени изменение температуры сказывается на параметре при окорке осины. Объяснение этому - максимальное значение характеристики внутреннего сцепления С осины и установленное влияние Т на величину С в соответствии с полученным соотношением (11).

Результаты расчета показали, что линейный характер зависимости (Т) удовлетворительно описывает процесс до значения Т=-7 оС. При дальнейшем понижении Т наблюдается нелинейный характер увеличения силы. Кроме этого, что при переходе через нулевую температуру по мере увеличения прочности коры проявляется локальный скачок (максимум) удельной силы . Это наглядно видно, например, (рис. 8) при сравнительном анализе результатов окорки осины и березы. Можно сделать вывод, что диапазон изменения Т при переходе через нулевую отметку влияет на изменение внутреннего состояния массива коры и развитие процесса его разрушения при роторной окорке.

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

Поскольку толщина коры определяется диаметром бревна, механизм разрушения слоев массива зависит от данного фактора и заслуживает более детального внимания.

Рис. 8. Характер изменения силы при переходе через Т=0оС:

1 - осина; 2 береза

Влияние диаметра бревна на развитие механизма роторной окорки. Анализ регрессионных известных зависимостей позволил получить связь между hк и dб для различных пород. Очевидно, что чем толще кора древесины, тем сложнее развитие механизма силового воздействия на нее короснимателя. Геометрические параметры dб и hк бревна предопределяют достижение необходимой удельной силы окорки . Помимо силовых факторов, необходимо учесть влияние диаметра бревна на кинематические параметры окорки. Толстые бревна рекомендовано окаривать на меньших скоростях подачи uп, снижая при необходимости частоту вращения ротора и увеличивая угол окорки. Это связано с тем, что при входе толстых бревен на больших скоростях наблюдаются динамические явления, приводящие к деформациям короснимателей и износу подшипников ротора.

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

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

С целью оптимизации режимов окорки целесообразно сортировать бревна по диаметру и запускать их в роторный станок вершиной вперед с тем, чтобы по мере роста dб плавно увеличивать силу прижима F1 (при постоянной скорости подачи uп) либо постепенно снижать значения uп при постоянной силе F1. Отмеченные обстоятельства диктуют необходимость оценки такого параметра как сбежистость бревна Сб и его влияния на установление допустимых диапазонов вариаций силовых и кинематических параметров окорки. Таким образом, актуальной является задача дополнения положений разработанной математической модели расчета параметров окорки учетом диаметра бревна и его возможных вариаций.

Расчеты выполнены для условий окорки влажных бревен в зимних условиях (T0оС). На рис. 9 представлены зависимости критерия разрушения коры S от диаметра dб, м. Графики соответствуют следующим значениям параметров: uп=0,3 м/с, F1=730 Н, =1,74 рад., =180 об/мин, W=100%, Т=-5оС, Kп=2.

Расчеты свидетельствуют, что при таких параметрах окорки сосны условие качественной окорки (8) выполняется для бревен диаметром до 0,45 м, для осины - до 0,35 м. Удельная сила окорки при этом составляет соответственно 16,42 и 17,21 кН/м. При больших значениях dб необходимо корректировать параметры окорки с целью выполнения условия (8).

В табл. 1 представлены результаты расчетов качественной окорки бревен сосны различного диаметра. Выполнение условия (8) достигалось путем изменения трех параметров uп, F1, Kп при постоянстве остальных.

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

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

Таблица 1. Показатели окорки бревен сосны различного диаметра

dб,м	hк,м	Kп	F1, Н	,кН/м	Ак, МПа/м
0,20	0,003	1	730	10,77	2,56
0,30	0,007	1	1000	13,17	1,32
0,45	0,017	2	730	16,42	0,63
0,50	0,021	2	900	19,37	0,59
0,55	0,026	2	1000	21,07	0,50
0,60	0,031	2	1080	22,28	0,43

Рис. 9. Влияние диаметра бревна на выполнение критерия качественной окорки:

1 - сосна; 2 осина

При постоянной величине Lк, задавшись различными значениями наибольших диаметров комля dк и вершины dв, сбежистость Сб определялась в процентах. Далее, в диапазоне допустимых значений диаметров бревна (dк, dв) с помощью программы генерации случайных чисел по закону равномерного распределения формировалась выборка диаметров di, для которых реализовывалась разработанная математическая модель и, в результате, формировались выборки расчетных значений параметров и показателей окорки (Xi). Статистическая обработка полученных выборок позволила установить значения их характеристик - математического ожидания M(Xi), дисперсии D(Xi) и коэффициента вариации v(Xi) в зависимости от параметра Сб. На рис. 10 для условий окорки сосны представлены результаты имитационного моделирования применительно к оценке вариации силовых характеристик (X1= F1, X2 =).

Полученные результаты позволяют, в частности, задавшись фактической величиной сбега, установить допустимые диапазоны отклонений силовых параметров от своих математических ожиданий. Так, при сбеге Сб = 10% и M() = 19,37 кН/м, среднее квадратичное отклонение () составило:

кН/м,

т.е. допустимый диапазон вариации удельной силы равен:

или =17,9520,79 кН/м. (12)

Рис. 10. Влияние сбега на вариации силовых параметров

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

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

Так, на основании обработки известных опытных данных физико-механических свойств сухой коры, в табл. 2 представлены в соответствующих масштабных единицах (м.е.) значения математических ожиданий (М), средних квадратичных отклонений () и коэффициентов вариаций (v, %) плотности коры k (м.е.- кг/м3), пределов прочности на сжатие сж и разрыв р (м.е.- МПа) для трех пород древесины. С повышением влажности вариация свойств коры будет возрастать.

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

Таблица 2. Статистические показатели свойств коры

Порода древесины	Статистические характеристики показателей свойств коры
Плотность сухой коры, k	Предел прочности на сжатие, сж	Предел прочности на разрыв, р
М		v	М		v	М		v
Сосна	370,5	74,5	20,1	6,2	1,7	27,4	1,5	0,5	33,3
Осина	608,0	67,0	11,0	13,3	4,1	30,6	11,0	1,1	9,6
Береза	770,0	57,5	7,5	19,2	3,5	18,2	2,3	0,3	11,1

Принимая во внимание, что технологически наиболее вариативным является силовой параметр - сила прижима F1, а угловой и кинематический параметры целесообразно выдерживать стабильными при окорке бревен в условиях заданной температуры Т, была принята схема моделирования (рис. 11), в соответствии с которой реализуется следующий алгоритм:

1. С помощью метода статистических испытаний (Монте-Карло) датчиком случайных чисел формируется выборка нормально распределенных в интервале (0,1) чисел i, i=100, которые подставляются в формулу общего вида:

, j=1,Е,4, (13)

где M(Xj) и (Xj) - математические ожидания и средние квадратичные отклонения трех показателей физико-механических свойств коры (табл. 2) и ее влажности W. При расчетах статистические характеристики влажности принимались: M(W)=40100% с шагом 20%; (W)=015% с шагом 5 %;

2. Выработанные числа Хj при заданной температуре Т принимаются переменными параметрами окорки, которые с учетом угловых, кинематических и силовых параметров используются в математической модели процесса разрушения коры. Выходным параметром является величина приведенного давления на границе раздела кора-древесина, которая сравнивается с характеристикой прочности коры на скалывание ск;

3. Если условие /ск 1 не выполняется, осуществляют корректировку силы прижима F1 в сторону увеличения так, чтобы М(F1) являлось величиной, достаточной для выполнения критерия разрушения (8).

Рис. 11. Схема реализации вариационного метода

Рис. 12. Корреляционное поле силы прижима

На рис. 12 представлено корреляционное поле силы прижима F1 в зависимости от вариаций влажности W вблизи заданных М(W)=100% и (W)=5% при следующих условиях окорки мерзлой сосны на станке ОК-63: Т=-5 оС; dб =0,4 м; =1,8 рад.; Kп=2. Коэффициент вариации силы составил v(F1)=11,7%. Обработка и обобщение результатов моделирования для других пород древесины в широком диапазоне изменения влажности позволили установить линейный характер влияния величины (W) на рост величины v(F1).

Полученный результат позволяет решать следующие задачи: задавшись предельно-допустимым значением (W), например 4%, определить допустимую величину v(F1)=8,47% и, в частности, при М(F1)=850Н установить величину (F1)=70 Н. Таким образом, диапазон изменения силы прижима F1 составит:

Н.

Обратная задача заключается в том, что, задавшись и выдерживая в определенном диапазоне значения силы F1, определяем допустимый диапазон изменения влажности бревна W, что позволяет осуществлять предварительную сортировку бревен по критерию стабилизации параметров окорки. Анализируя полученные результаты, необходимо отметить, что разработанный вариационный метод позволяет произвести расчет допустимых диапазонов отклонений силовых параметров окорки от заданных значений и тем самым обеспечить стабилизацию технологических показателей в условиях изменчивости физико-механических свойств и характеристик состояния влажности лесоматериала.

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

Одним из основных факторов, влияющих на качество окорки является влажность коры и, следовательно, показатель влажности ее основных структурных составляющих - корки и луба. При этом необходимо учитывать различия W как для коры в целом, так и для корки и луба, а также влияние этого различия на прочностные характеристики материалов. Разрушение более сухой и хрупкой корки существенно отличается от процесса разрушения более влажной и пластичной среды лубяного слоя. Влажность корки существенно меньше влажности луба для свежесрубленной древесины, тогда как для сплавной древесины указанное различие меньше. Влажность коры свежесрубленной древесины у различных пород колеблется в диапазоне W=71139%, причем для сосны влажность луба (Wл) достигает 180%, а корки (Wк) 40%, у ели соответственно - 120 и 45%.

Принимая влажность коры за 100%, можно оценить, что для хвойных свежесрубленных бревен имеют место следующие соотношения:

Wк = 0,45- 0,55 W, Wл = 1,55-0,75 W (14, а);

для сплавных материалов:

Wк =0,7-0,8 W, Wл =1,1-1,2 W (14, б);

Наиболее трудными представляются условия селективной окорки свежесрубленной березы, поскольку показатели влажности составляют:

Wк = 0,15-0,25W , Wл =1,1-1,15 W (14, в).

Рассмотрим процесс селективного отделения коры от древесины с помощью скребкового короснимателя с углом окорки >/2 (рис. 13) на базе двухроторного станка типа 2ОК63-1.

Рис. 13. Схема разрушения массива коры:

1-древесина; 2- кора (корка, луб); 3-коросниматели; 4- окоренная поверхность;5-сучок

При взаимодействии короснимателя первого ротора с пеньком сучка возникает дополнительная растягивающая сила Fр, связанная с нормальной силой Fс известным соотношением:

, (15)

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

= о + , (16)

о установочный угол, .

При работе с углами окорки >/2 огибание пенька возможно при условии достижения критического угла значения:

кр =с + >arcctg б,

где , dc и hc диаметр и высота сучка.

Тогда предельное значение коэффициента усиления:

В зависимости от значений R, db, и L второе угловое слагаемое в соотношении (16), которое примем в качестве угловой поправки, может принимать как положительные, так и отрицательные значения, что существенно влияет на величину угла окорки . При работе скребковыми короснимателями параметры R, db, и L принимаются таковыми, чтобы угловая поправка >0.

Расчет влияния длины L на величину поправки при различных значениях db показал, что, существует некоторое критическое значение величины Lкр 0,3 м, до и после которого наблюдается прямо противоположное влияние dб на угловую поправку . При использовании короснимателей с длиной L, сопоставимой с Lкр, при любом диаметре бревна угловая поправка стремится к постоянной величине, равной: =const=0,7 рад. (40,1 град.). Анализ показал, что максимальные значения поправки имеют место при малых диаметрах бревна и больших значениях длины короснимателя L.

На практике нашли применение два типа короснимателей: у первых жесткость в направлении подачи превышает жесткость в плоскости ротора, у вторых - наоборот. В первом случае применяемое на практике значение о = 60-70 град. следует считать завышенным и приводящим к частым поломкам короснимателей, поскольку при таких значениях установочного угла огибание сучка сопровождается критическими деформациями самого короснимателя. Во втором случае изгиб происходит в направлении подачи и коросниматель обходит пенек сучка сбоку. Оптимальной представляется такая конструкция короснимателя, когда передняя грань выполнена с переменным значением о. Указанные обстоятельства требуют более детального исследования особенностей развития процесса селективной окорки применительно к работе внешнего ротора, коросниматели которого непосредственно взаимодействуют с пеньками сучков.

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

; ; . (17)

При выполнении условий (17) фиксируются достигнутые угловые ( и о), силовые (F1 ,Fс) и кинематические (скорость uп) параметры. В качестве интегрального показателя силовых затрат принимается величина удельной силы окорки , кН/м. Для реализации разработанной модели выбраны три породы: сосна, ель и береза при условии окорки свежесрубленного (W=100130%) бревна на станке типа 2ОК63-1. Технические характеристики станка: скорость подачи обоих роторов одинакова и изменяется в диапазоне uп=0,21,0 м/с; F1=7302900 Н; число короснимателей -6; число -135, 180 об/мин; частота вращения роторов - 2,5-5с-1; dб=0,1-0,55м; L =0,2-0,38 м, задний и передний углы - соответственно 0,785 и 0,25 рад.; =1,742,35 рад., R=0,47. Толщина корки hкр варьировалась в пределах от 5 до 50% величины hк. Толщина луба принималась как hкл=hк-hкр. Величины пределов прочности коры сосны, ели и березы ск поперек волокон соответственно: ск= 0,33, 0,78 и 1,42 МПа. В качестве базового варианта расчета для станка 2ОК63-1 принимался случай рассмотрения массива коры как единого целого. Отталкиваясь от базового варианта были выполнены расчеты при условии реализации селективных принципов окорки.

В табл. 3 при исходной влажности коры W=100%, и перерасчете влажностей корки и луба в соответствии с соотношениями (14, а)(14, в) представлены результаты расчетов удельной силы окорки , а также соответствующих удельных сил и , необходимых для разрушения слоев корки и луба.

Таблица 3. Результаты расчетов усилий окорки

Порода	, кН/м	, кН/м	, кН/м
Сосна	8,03	16,02	1,41
Ель	13,73	32,51	6,11
Береза	34,83	178,1	27,5

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

Исследуем влияние толщины слоев корки и луба на трансформацию и распределение удельных сил окорки. На рис. 14 представлены результаты выполненных исследований, где по оси абсцисс отложена толщина корки (в % к толщине коры), а на оси ординат - безразмерная величина, равная отношению полученной средневзвешенной удельной силы окорки =hкр+hкл к . Как следует из анализа полученных данных эффективность селективного метода особенно явно проявляется при окорке такого трудоемкого лесоматериала как береза. Лишь при значении hкр=0,05hк имеет место равенство = .

Рис. 14. Влияние толщины корки на соотношение удельных сил при окорке:

1 - сосны; 2 Цели; 3 - березы

При больших значениях толщины корки березы наблюдается существенный рост удельных затрат на ее разрушение, что вызывает необходимость увеличения нормальной силы прижима Fс. Это, в свою очередь, при взаимодействии короснимателя с сучком, обусловливает возникновение растягивающей силы Fp, определяемой с помощью соотношения (15) с учетом коэффициента усиления Ку.

На рис. 15 представлена зависимость коэффициента Ку от диаметра сучка dс=0,005-0,025м и диаметра бревна dб =0,25-0,55м. Данные соответствуют значениям угла кр=0,251,0 рад. (1457 град.), L=0,25м, hc= dс/2. Как видно, при снижении dб с ростом характерного размера dс коэффициент усиления возрастает и при критических величинах может достичь значений 45, что негативно сказывается на физическом состоянии короснимателя и условиях окорки в целом.

На рис. 16 представлено влияние критического угла кр на величину Ку при L=0,25, dс=0,02, hc=0,01 и dб =0,25м.

Полученный результат с учетом данных табл. 3 о кратном увеличении значений при разрушении корки, свидетельствует о целесообразности применения параметров окорки с тенденцией роста угловых значений кр при соблюдении условия: кр>arcctg б. Следует также отметить, основываясь на результатах выполненных исследований, что с ростом параметра L коэффициент усиления Ку снижается до 1, а при малых диаметрах бревна - до 0,5 и менее. Интенсивность и величина этого снижения в зависимости от dб представлены на рис. 17.

Рис. 15. Влияние диаметров бревна и сучка на величину коэффициента усиления Ку

Рис. 16.Зависимость коэффициента Ку от критического угла кр

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

Рис. 17. Зависимость величины Ку от длины L: 1 dб=0,25 м; 2 dб=0,35 м; 3 dб=0,55 м

Разработанные математические модели реализованы в виде макроса в программе Microsoft Excel, созданного в виде автоматизированного рабочего места Расчет параметров процесса роторной окорки.

4. Объект, аппаратура, методика и условия проведения

экспериментальных исследований

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

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

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

5. Результаты теоретических и экспериментальных исследований

В результате исследования математических моделей, которая производилась при помощи прикладных программ MatСad 2005. Professional и Excel 2005, входящей в пакет прикладных программ Office XP Professional для операционной системы Windows ХР, установлено, что значение толщин массива коры в целом, отдельных ее слоев пробки и луба, в зависимости от диаметра бревна, наиболее точно описываются полиномиальной зависимостью с коэффициентом аппроксимации от 0,737 до 0,937. Результаты экспериментальных исследований отличаются от теоретических не более чем на 8%. Сравнение результатов теоретических и экспериментальных исследований, значение коэффициента неравномерности с определением доверительного интервала, позволяет утверждать, что разработанные математические модели адекватны объекту исследования. Автоматизированное рабочее место Расчет параметров процесса роторной окорки, разработанное на основании математических моделей, позволяет автоматизировать и существенно сократить трудоемкость расчета основных параметров роторной окорки круглых лесоматериалов в различных условиях.

6. Оценка технико-экономической эффективности внедрения селективной окорки круглых лесоматериалов на предприятиях лесопромышленного комплекса

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

Общие выводы.

Повышение качества окорки круглых лесоматериалов является одним из основных направлений научно-технического прогресса лесной промышленности. Данное направление научных исследований входит в створ Перечня Приоритетных направлений развития науки, технологий и техники в Российской Федерации, утвержденного Президентом РФ 21 мая 2006 г. Пр-843 (пункт Рациональное природопользование).
Отделенные друг от друга в процессе окорки пробковый и лубяной слои коры основных лесообразующих пород являются ценными полуфабрикатами для производства различных товаров для народного хозяйства. Особенно важно иметь возможность селективной окорки хвойных пород древесины, поскольку до настоящего времени не разработано технологий эффективной утилизации отходов окорки хвойных деревьев, включающих пробку и луб.
Разработанная математическая модель роторной окорки с использованием тупых короснимателей, основанная на механизмах деформаций элементарного объема сплошной среды и построении обобщенных диаграмм Мора, позволяет оценить нормальные и касательные компоненты приведенного давления в толще массива коры и условия его разрушения на границе с древесиной различных пород деревьев.
Установленные на основе реализации математической модели количественные соотношения зависимости величины приведенного давления от угловых, силовых и кинематических параметров управления процессом окорки позволяют произвести сравнительный анализ степени их влияния на достижение заданного качества отделения коры от древесины в конкретных технологических условиях.
Рассматривая массив коры как трехкомпонентную среду, содержащую твердую, жидкую и газообразную компоненты, обоснован детерминированный метод оценки влияния влажности, температуры, диаметра бревна, величины сбега на достижение необходимой и достаточной удельной силы окорки.
Учитывая, что в реальных производственных условиях факторы влияния и параметры управления процессом окорки являются переменными величинами и варьируются в широких пределах, разработанный вариационный метод на основе имитационного моделирования устанавливает допустимые диапазоны их вариаций, которые обеспечивают стабилизацию качественных показателей окорки.
Разработанный метод расчета параметров селективной роторной окорки лесоматериалов при использовании типоразмерного ряда двухроторных окорочных станков, учитывающий отличительные особенности разрушения слоев корки и луба различной толщины, позволяет обосновать такой выбор угловых и геометрических характеристик короснимателя в сочетании с силовыми параметрами, обеспечивающих выполнение процесса селективной окорки.
Установлен линейный характер возрастания относительной прочности коры с увеличением относительного диаметра, т.е. масштабный эффект среды разрушения является существенным. Проявление масштабного эффекта в зимнее время существенно выше, чем при окорке летних бревен.
Запатентованный стенд для измерения и регистрации процессов роторной окорки круглых лесоматериалов показал надежную и точную работу в процессе проведения экспериментальных исследований, что позволяет рекомендовать его для дальнейших исследований.
Значение толщин массива коры в целом, отдельных ее слоев пробки и луба, в зависимости от диаметра бревна, наиболее точно описываются полиномиальной зависимостью с коэффициентом аппроксимации от 0,737 до 0,937.
Сравнение результатов теоретических и экспериментальных исследований, значение коэффициента неравномерности с определением доверительного интервала, позволяет утверждать, что разработанные математические модели адекватны объекту исследования.
Автоматизированное рабочее место Расчет параметров окорки, разработанное на основании математических моделей, позволяет автоматизировать и существенно сократить трудоемкость расчета основных параметров роторной окорки, что позволяет рекомендовать его для включения практику предприятий лесопромышленного комплекса, отраслевых НИИ и КБ, и в учебный процесс вузов, для студентов специальности 250401 Лесоинженерное дело, и направления 250300 Технология и оборудование лесозаготовительных и деревообрабатывающих производств.
Согласно петле качества, возможность эффективной утилизации отходов окорки является неотъемлемой составляющей общей оценки качества работы окорочного оборудования, наряду с площадью и шероховатостью поверхности лесоматериала очищенного от коры.
Решение экономических аспектов, связанных с эффективной утилизацией отходов окорки имеет важные специфические особенности: формирование затрат при комплексном использовании древесной массы и коры, формирование стоимости и потребительской стоимости при оценке различных видов коры, обоснование цен на продукцию, получаемую при утилизации отходов окорки, организация маркетинга, и т.д.
В зависимости от того, какой потенциал отходов окорки используется (материальный, энергетический и др.), отходы окорки целесообразно подразделять на материальные (сырье для дальнейшей переработки по видам продукции) и энергетические.
Финансовый успех предприятий ЛПК, внедряющих селективный метод окорки круглых лесоматериалов, будет определяться спросом отходов окорки, как технологического сырья, так и энергетической составляющей в топливном балансе самого предприятия. При этом рекомендуется предусматривать методический учет нормативных параметров отходов окорки по видам.
Расчеты показывают, что на среднем лесопромышленном складе (бирже сырья потребителя), при минимальных отпускных ценах на луб, и утилизацию пробковой части коры в собственной котельной, при удаленности потребителей лубяного слоя на расстояние 5075 км, срок окупаемости инвестиций в мероприятия, связанных с внедрением селективной окорки, составит около 3,5 лет.

Основное содержание диссертации опубликовано в следующих работах:

Газизов А.М., Бойко В.В. Ваньков A.M. Опыт эксплуатации окорочных станков // Лесная промышленность, 1988, № 9, С. 19.
Газизов А.М. К вопросу обоснования параметров станка для окорки хлыстов // Лесная промышленность, 1997, № 4, С. 2324.
Газизов А.М., Шапиро В.Я., Григорьев И.В. Моделирование процесса разрушения коры при роторной окорке древесины // Вестник Красноярского государственного аграрного университета, № 5. 2008 г. С. 271-279.
Газизов А.М., Шапиро В.Я., Григорьев И.В. Влияние влажности на развитие процесса разрушения коры при роторной окорке // Вестник Московского государственного университета леса Ц Лесной вестник. Вып. № 6 (63). Издательство МГУЛ. 2008. С. 129133.
Газизов А.М., Шапиро В.Я., Григорьев И.В. Исследование процесса разрушения коры при роторной окорке бревен различного диаметра // Справочник. Инженерный журнал. № 2. 2009. С. 4550.
Газизов А.М., Шапиро В.Я., Григорьев И.В. Вариационный метод расчета и стабилизации параметров роторной окорки // Справочник. Инженерный журнал. № 7. 2009. С. 4751.
Шапиро В.Я., Григорьев И.В., Газизов А.М. Анализ теории расчета параметров процесса окорки древесины в роторных окорочных станках //Известия Санкт-Петербургской лесотехнической академии. 2009. № 186, С. 110-120.
Газизов А.М., Григорьев И.В., Гумерова О.М. Повышение качества окорки лесоматериалов // Вестник Красноярского государственного аграрного университета, № 10. 2009 г. С. 132-141.
Шапиро В.Я., Григорьев И.В., Газизов А.М. Обоснование и расчет параметров селективной роторной окорки лесоматериалов//Известия СПбГЛТА. 2010. № 190, С. 104-117.
Газизов А.М., Григорьев И.В., Гумерова О.М. Моделирование процесса разрушения коры при окорке резанием//Известия СПбГЛТА. 2010. № 193, С. 18-24.
Газизов А.М., Григорьев И.В., Кацадзе В.А., и др. Повышение эффективности механической окорки лесоматериалов. Монография. СПб.: ЛТА, 2009. - 240 с.
Газизов А.М., Дмитренков Г.Л. Определение нагрузок, действующих на механизм подачи роторных окорочных станков типа ОК 100-2 и ОК 80-2 / Переработка и энергоиспользование низкокачественной древесины. Труды ЦНИИМЭ. - Химки, 1989, С. 60 - 64.
Симонов М.Н., Минакова Т.В., Газизов А.М. Режимы круглогодовой окорки лесоматериалов с отбором коры на технологическое сырье / Переработка и энергоиспользование низкокачественной древесины. Труды ЦНИИМЭ. - Химки, 1989, С. 10-17.
Газизов А.М. Экспериментальное исследование нагрузок действующих на механизм подачи роторного станка / ВНИПИЭИ и Леспром 15.06.1990 №2670. Химки, 1990.
Газизов А.М. Определение внешних нагрузок, действующих на механизм подачи в момент захода хлыста в ротор окорочного станка / Материалы научной конференции профессорско-преподавательского состава БГАУ. Уфа, БГАУ, 1993, С. 53-54.
Газизов А.М., Посняк В.К. Обоснование параметров окорочного станка / Материалы научной конференции профессорско-преподавательского состава БГАУ Уфа, БГАУ, 1993, С. 54-55.
Газизов А.М. Обоснование параметров оборудования для окорки хлыстов и полухлыстов. Уфа: БГАУ, 1995. 16 с.
Газизов А.М. Эффективность перехода на окорку древесины в хлыстах // Лесной экономический вестник, Москва, 1997, № 2 С. 18-20.
Газизов А.М. Окорка древесины в хлыстах // Сельские узоры. Республиканский аграрный журнал, Уфа, 1997, № 3, С. 2627.
Газизов А.М. Состояние и перспективы механизации процессов лесоэсплуатации в Республике Башкортостан // Актуальные проблемы лесного комплекса: Сборник информационных материалов международной научной конференции Лес - 2000. Вып.2. - Брянск: РИО БГИТА, 2000, С. 2930.
Газизов А.М. Состояние и тенденции развития лесопромышленного комплекса Республики Башкортостан / Принципы формирования высокопродуктивных липняков, Уфа, 2000. С. 3538.
Газизов А.М. Основные направления комплексного использования древесины / Материалы научно-практической конференции Интеграция науки, образования и производства для развития лесного хозяйства и лесопромышленного комплекса, Воронеж, 2001. С. 1922.
Газизов А.М. Пути использования древесины без потерь / Материалы всероссийской научно-практической конференции Повышение эффективности и устойчивости развития АПК, Уфа, 2005. С. 3649.
Газизов А.М., Шапиро В.Я., Григорьев И.В. Анализ современных методик расчета основных параметров окорки режущим инструментом и пути их уточнения / Сборник трудов международной научно-практической конференции Актуальные проблемы лесного комплекса под общей редакцией А.Я. Панфилова. Вып. № 21. Брянск 2008 г. С. 231235.
Газизов А.М., Шапиро В.Я., Григорьев И.В., Гумерова О.М. Пути совершенствования методик расчета технологических параметров механической окорки круглых лесоматериалов / Деп. рукописи: библиогр. указ. ВИНИТИ, 30.09.08. № 547 - В 2008. 25 с.
Газизов А.М., Григорьев И.В., Гумерова О.М. Анализ технических решений по повышению качества поштучной механической окорки круглых лесометриалов / Деп. рукописи: библиогр. указ. ВИНИТИ, 30.09.08. № 546 - В 2008. - 23 с.
Газизов А.М., Григорьев И.В., Гумерова О.М., Ильюшенко Д.А. Аппаратное обеспечение экспериментальных исследований селективной роторной окорки круглых лесоматериалов / Технология и оборудование лесопромышленного комплекса Межвузовский сборник научных трудов. Выпуск 2. СПб.: ЛТА 2008 г. С. 4249.
Газизов А.М., Григорьев И.В., Гумерова О.М., Ильюшенко Д.А. Анализ типоразмерных рядов роторных окорочных станков и их инструментального оснащения / Деп. рукописи: библиогр. указ. ВИНИТИ, 16.02.09. № 75 - В 2009. - 21 с.
Шапиро В.Я., Григорьев И.В., Гумерова О.М., Газизов А.М. Совершенствование метода расчета параметров окорки круглых лесоматериалов трением / Материалы первой международной научно-практической Интернет конференции Леса России в XXI веке. СПб.: ЛТА, 2009. С. 218222.
Газизов А.М. Математическая модель процесса роторной окорки и исследование способов управления его параметрами / Материалы первой международной научно-практической Интернет конференции Леса России в XXI веке. СПб.: ЛТА, 2009. С. 162167.
Газизов А.М. Модель влияния температуры на эффективность роторной окорки круглых лесоматериалов / Технология и оборудование лесопромышленного комплекса Межвузовский сборник научных трудов. Выпуск 3. СПб.: ЛТА 2009 г. С. 4755.
Григорьев И.В., Газизов А.М., Теппоев А.В. Стенд для исследования процесса окорки лесоматериалов. Патент на полезную модель № 76597 опубл. 27.09.2008 Бюлл. № 27.
Газизов А.М., Григорьев И.В., Гумерова О.М. Обоснование показателей качества работы окорочных станков / Технология и оборудование лесопромышленного комплекса Межвузовский сборник научных трудов. Выпуск 3. СПб.: ЛТА 2009 г. С. 8185.
Газизов А.М., Григорьев И.В., Ильюшенко Д.А. Влияние строения и состояния коры на параметры процесса селективной окорки / Материалы международной научно-технической конференции Актуальные проблемы развития лесного комплекса. Вологда: ВоГТУ, 2009. С. 8284.
Газизов А.М., Гумерова О.М., Ильюшенко Д.А. Сравнительный анализ способов окорки круглых лесоматериалов / Молодая мысль: Наука. Технологии. Инновации: Материалы Межвузовской научной конференции студентов, магистрантов, аспирантов и молодых ученых. - Братск: ГОУ ВПО БрГУ, 2009. - С 7478.
Газизов А.М., Гумерова О.М., Ильюшенко Д.А. Пути повышения качества окорки лесоматериалов / Материалы Международной научно-практической конференции Современные проблемы лесозаготовительных производств, производства материалов и изделий из древесины: пиломатериалы, фанера, деревянные дома, заводского изготовления, столярно-строительные изделия. СПб.: НП НО - МТД, 2009. Том. 2. С.179-184.
Газизов А.М. Влияние технологических факторов на формирование напряжений в массиве коры при роторной окорке древесины / Сборник материалов Всероссийской научно-практической конференции докторантов, аспирантов и соискателей. -аМайкоп, 2008. С. 5463.
Мурашкин Н.В., Газизов А.М., Гумерова О.М., Куприянова Е.А. Экономические вопросы использования отходов окорки в условиях рыночных отношения на предприятиях ЛПК / Технология и оборудование лесопромышленного комплекса Межвузовский сборник научных трудов. Выпуск 4. СПб.: ЛТА 2009 г. С. 4653.
Газизов А.М., Григорьев И.В., Ильюшенко Д.А., Гумерова О.М. Статистические характеристики толщин слоев коры на пиловочных бревнах / Материалы второй международной научно-практической Интернет конференции Леса России в XXI веке. СПб.: ЛТА, 2009. С. 96100.
Газизов А.М. Перспективы рационального использования лесных ресурсов путем внедрения окорки древесины / Всероссийская научно-практическая конференция "Интеграция аграрной науки и производства: состояние, проблемы и пути решения". Уфа, 2008.-с.41-43"

Просим принять участие в работе диссертационного Совета Д.212.008.01 или прислать Ваш отзыв на автореферат в двух экземплярах с заверенными подписями по адресу: 163002, Архангельск, набережная Северной Двины 17 Северный (Арктический) федеральный университет.

Авторефераты по всем темам >> Авторефераты по разное

Blog