close

Вход

Забыли?

вход по аккаунту

?

Методы и средства таможенного контроля плотности древесины.

код для вставкиСкачать
На правах рукописи
СИМОНЕНКО Антон Анатольевич
МЕТОДЫ И СРЕДСТВА ТАМОЖЕННОГО КОНТРОЛЯ
ПЛОТНОСТИ ДРЕВЕСИНЫ
Специальность 05.11.13 - Приборы и методы контроля
природной среды, веществ, материалов и изделий
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
кандидата технических наук
Санкт-Петербург – 2014
Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Национальный минерально-сырьевой университет
«Горный».
Научный руководитель:
кандидат технических наук, доцент
Павлов Игорь Валерьевич
Официальные оппоненты:
Афонин Пётр Николаевич – доктор технических наук, доцент, Санкт-Петербургский имени В.Б. Бобкова филиал ГКОУ ВПО
«Российской таможенной академии», кафедра технических средств
таможенного контроля и криминалистики, профессор, и.о. заведующего кафедрой
Коновалов Сергей Ильич – кандидат технических наук, доцент, ФГБОУ ВПО «Санкт-Петербургский государственный электротехнический университет «ЛЭТИ» им. В.И. Ульянова (Ленина)»,
кафедра электроакустики и ультразвуковой техники, доцент
Ведущая организация – ФГБОУ ВПО «Санкт-Петербургский
государственный архитектурно-строительный университет»
Защита состоится 08 октября 2014 года в 16 часов 30 мин. на
заседании диссертационного совета Д 212.224.07 в Национальном
минерально-сырьевом университете «Горный» по адресу: 199106,
Санкт-Петербург, Васильевский остров, 21 линия д.2, ауд. 7212.
С диссертацией можно ознакомиться в библиотеке Национального минерально-сырьевого университета «Горный» и на сайте
www.spmi.ru
Автореферат разослан 11 июля 2014 года
УЧЕНЫЙ СЕКРЕТАРЬ
диссертационного совета
Фокин
Андрей Сергеевич
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность работы. Российская Федерация является одной из
крупнейших стран заготовителей и экспортёров лесопродукции. Активная
интеграция России в мировую торговлю требует совершенствования таможенного администрирования. Утверждаются новые федеральные законы,
программы и распоряжения, направленные на развитие внешнеэкономической деятельности. Неуклонно растут объемы контролируемых на границе
грузов и при этом сокращаются сроки проведения таможенных операций и
таможенных процедур.
При экспорте лесоматериалов в документах на груз указываются значения основных параметров древесины, таких как наименование пород,
объем, масса, габариты, влажность и, зачастую, плотность как одно из важных для целлюлозно-бумажной промышленности значений. Для контроля
всех параметров лесопродукции, за исключением плотности, таможенные
органы располагают необходимым набором технических средств таможенного контроля и методик измерений. При этом определение плотности древесины осуществляется на основе табличных методов или в соответствии с
ГОСТ 16483.1-84 «Древесина. Метод определения плотности». В свою очередь табличные методы не гарантируют точность измерений, а измерения
по ГОСТу требуют соответствующих условий и больших временных затрат. Необходимо отметить, что в некоторых случаях, измерение плотности
древесины может быть единственным способом для определения объема
груза при условии наличия весов. И наоборот, в случае их отсутствия, но
при возможности определения объема, значение плотности позволяет рассчитать массу груза, что немаловажно, особенно, при железнодорожных
перевозках. Одним из важнейших достоинств плотности, как показателя
качества древесины, является ее универсальность. Так, например, плотность
древесного сырья определяет важнейший экономический показатель – выход целлюлозы, также она дает представление о механических свойствах
древесины, что можно использовать для прочностной сортировки пиломатериалов. Таким образом, для объективной оценки задекларированных данных и качества лесопродукции, наряду с определением породы, объема и
влажности древесины необходимо определять её плотность.
Опыт показывает, что внедрение современных технологий способствует повышению эффективности таможенного контроля и предупреждает
3
нарушения таможенного законодательства. В связи с этим, актуальность
работы обусловлена необходимостью решения практических задач по оперативной, в полевых и промышленных условиях, оценке плотности древесины, с целью выявления нарушений таможенного законодательства, что
поспособствует повышению эффективности таможенного контроля лесопродукции, экспортируемой с территории Российской Федерации.
Степень разработанности. Работы по исследованию основных физико-механических характеристик древесины проводили Ю.М.Иванов,
В.А.Баженов, Ф.Ф.Садовский, А.А.Рабинович, Н.Ф.Гусев, Н.И.Миронов,
А.Н.Митинский, Л.М.Перелыгин, О.И.Полубояринов, Б.Н.Уголев,
Е.К.Ашкенази,
В.Д.Никишов,
А.М.Боровиков,
А.С.Сапожников,
В.В.Тулузаков, А.Н.Кириллов, В.Н.Волынский и др. Работы в области применения ультразвуковых методов для контроля древесины и изделий из неё
проводили Б.К.Лакатош, Ю.К.Сергиенко, В.И.Яковлев, А.И.Горбунов,
В.И.Федюков, Ф.Ф.Легуша, А.А.Ерофеев, С.В.Скрипец, А.А.Ефимов,
П.М.Мазуркин,
Б.А.Староверов,
Е.В.Саликова,
Е.Б.Темнова,
А.Н.Чубинский, А.А.Федяев, И.Я.Лиманов и др.
Проанализировав работы вышеперечисленных авторов можно сделать вывод, что пристальное внимание в основном уделяется вопросам разработки эффективных методов неразрушающих испытаний применительно
к сортировке пиломатериалов по прочности и к выявлению пороков и дефектов древесины на поточной линии. При этом решений по достоверной
оценке плотности лесопродукции, подходящих для применения в условиях
проведения таможенного контроля, не реализовано.
Цель и задачи исследования. Разработать методы и средства таможенного контроля плотности древесины в широком диапазоне влажности. В
рамках реализации поставленной цели сформулированы и решены следующие задачи:
1. Исследовать зависимости параметров акустических колебаний,
распространяющихся в древесине, от её плотности, породы и влажности, с
учётом времени и места заготовки древесины, и установить возможность
контроля плотности древесины по акустическим параметрам;
2. Теоретически обосновать и экспериментально подтвердить возможность определения плотности древесины в условиях таможенного кон-
4
троля с применением низкочастотного ультразвукового метода и метода
свободных колебаний;
3. На основании проведённых экспериментальных исследований
разработать методы и средства таможенного контроля плотности различных
пород древесины в широком диапазоне влажности;
4. Разработать прикладное программное обеспечение в основу
функционирования которого должны быть положены результаты практического применения разработанных методов и средств, а также дать рекомендации по их применению в таможенной практике.
Идея работы. Установление корреляционных зависимостей между
плотностью древесины различных пород в широком диапазоне влажности,
вычисленной классическим методом по результатам измерения массы и
объема, и приведенной скоростью распространения акустических волн, а
также скоростью распространения ультразвуковых колебаний позволит разработать новые методы контроля плотности древесины. Новые методы и
средства повысят эффективность таможенного контроля лесо- и пиломатериалов, а их применение с использованием персонального компьютера автоматизирует оформление результатов, что усовершенствует систему таможенных операций при комплексной оценке основных параметров лесопродукции.
Научная новизна:
1. Установлены зависимости, учитывающие влияние влажности
различных пород древесины на результаты определения плотности древесины акустическими методами контроля, на основе которых построены градуировочные зависимости «плотность древесины – скорость ультразвуковых колебаний», а также «плотность древесины – приведённая скорость
распространения акустических волн» в широком диапазоне значений массовой доли влаги для различных пород Северо-Западного, Южного и Сибирского регионов.
2. Разработаны методы и средства неразрушающего контроля плотности древесины в широком диапазоне влажности с применением низкочастотного ультразвукового метода и метода свободных колебаний, которые
учитывают условия проведения таможенного контроля лесо- и пиломатериалов.
5
Теоретическое и практическое значение работы. Установленные
на основании обработки и регрессионного анализа статистических данных
по определению плотности и акустических характеристик образцов древесины зависимости между плотностью древесины и приведенной скоростью
распространения акустических волн, а также между плотностью древесины
и скоростью распространения ультразвуковых колебаний расширяют область применения низкочастотных акустических методов контроля.
На основе экспериментальных исследований разработаны и аттестованы методы определения плотности лесо- и пиломатериалов в широком
диапазоне влажности, включающие в себя градуировочные зависимости для
наиболее часто предъявляемых к таможенному контролю пород древесины.
Результаты диссертационной работы используются в таможенных
органах Российской Федерации при таможенном контроле лесо- и пиломатериалов лиственных и хвойных пород древесины, а также могут быть использованы в лесозаготовительной и лесоперерабатывающей отраслях.
Методология и методы исследований. Исследования проводились
с использованием современных методов и средств научного поиска, опираясь на базовые положения науки о древесине, методов акустического контроля и статистических методов обработки экспериментальных данных. В
работе использован системный подход с применением методов планирования, моделирования, анализа и обработки информации.
Положения, выносимые на защиту:
1. Значения скорости распространения ультразвуковых колебаний и
приведённой скорости распространения акустических волн, полученные
при последовательных измерениях влажности и плотности древесины, с
достаточной точностью описываются математической моделью в виде регрессионного уравнения, что позволяет построить градуировочные зависимости «плотность-скорость» для различных пород древесины в широком
диапазоне влажности;
2. Градуировочные зависимости «плотность-скорость», установленные с применением метода свободных колебаний и низкочастотного ультразвукового метода продольно-поперечным способом прозвучивания,
обеспечивают достаточно точное определение физически усреднённых значений плотности древесины различных пород в широком диапазоне влажности;
6
3. Разработанные методы и средства определения плотности древесины учитывают условия применения и принципиально отличаются
от классической реализации существующих методов возможностью
в короткий промежуток времени объективно оценить значения плотности
лесо- и пиломатериалов в широком диапазоне влажности без необходимости их перемещения или изготовления контрольных образцов.
Достоверность научных результатов (защищаемых положений)
подтверждена значительным объемом экспериментальных данных, хорошей сходимостью и воспроизводимостью результатов измерений, а также
положительными результатами опытной эксплуатации и внедрением в таможенную практику.
Реализация работы. На основе проведенных теоретических и экспериментальных исследований создан модуль определения плотности древесины с аккредитованными методиками измерений, который входит в состав
многофункциональных приборов таможенного контроля лесопродукции,
принятых на снабжение таможенных органов России.
Личный вклад автора: анализ разработанности темы, обзор средств
измерений, разработку общей и экспериментальной методик исследований
с последующей оценкой и обоснованием полученных результатов автор
выполнил самостоятельно. Проведение экспериментов, полевых работ, аттестационных мероприятий и опытно-промышленных испытаний осуществлялось с непосредственным участием автора.
Апробация работы. Основные положения диссертационной работы
докладывались и обсуждались на XIII международной научно-практической
конференции молодых ученых, студентов и аспирантов «Анализ и прогнозирование систем управления» (Санкт-Петербург, 2012), 2-й международной научно-практической конференции «Современное машиностроение,
наука и образование» (Санкт-Петербург, 2012), XIХ международной научно-практической конференции студентов, аспирантов и молодых ученых
«Современные техника и технологии» (Томск, 2013) и на II международной
научно-практической конференции «Инновационные системы планирования и управления на транспорте и в машиностроении» (Санкт-Петербург,
2014).
Получен патент № 2 449 265 Российская Федерация, МПК G01N
29/07 «Способ и устройство определения плотности древесины».
7
По теме диссертации опубликовано 9 печатных работ, в том числе 3 в
журналах, рекомендованных ВАК при Минобрнауки России.
Структура и объем работы. Диссертация состоит из введения, пяти
глав, заключения, списка сокращений и условных обозначений, списка терминов, списка литературы и приложений. Содержание работы изложено на
194 страницах машинописного текста, включая 35 рисунков, 42 таблицы и
список литературы из 122 наименований.
Во введении приведено обоснование актуальности выбранной темы
диссертации, сформулированы цель работы и задачи исследования, отражена научная новизна, изложены сведения об апробации и показана практическая ценность работы.
В первой главе выполнен анализ современных методов определения
плотности древесины и степень разработанности темы в области применения акустических методов контроля лесопродукции. Рассмотрены особенности таможенного контроля лесоматериалов с применением технических
средств.
Во второй главе определены требования к методам и средствам таможенного контроля плотности древесины, представлены результаты предварительной оценки методов с выбором оптимальных параметров и аппаратуры.
В третей главе приведены результаты экспериментальных исследований по выявлению корреляционной зависимости между плотностью древесины и скоростью распространения акустических колебаний.
В четвёртой главе уточняются условия выполнения измерений и
оценивается соответствие разработанных методик требованиям к показателям точности с оценкой суммарной погрешности измерений. Приводятся
результаты реализации методов в виде программного обеспечения.
В пятой главе рассмотрены основные правонарушения при экспорте
лесоматериалов и приведены результаты применения разработанных методов и средств в таможенных органах.
В заключении представлены основные научные и практические результаты исследований в соответствии с целью и поставленными задачами.
ЗАЩИЩАЕМЫЕ НАУЧНЫЕ ПОЛОЖЕНИЯ
1. Значения скорости распространения ультразвуковых колебаний и приведённой скорости распространения акустических волн, полу-
8
ченные при последовательных измерениях влажности и плотности древесины, с достаточной точностью описываются математической моделью в
виде регрессионного уравнения, что позволяет построить градуировочные
зависимости «плотность-скорость» для различных пород древесины в
широком диапазоне влажности
Поскольку определение плотности бревен классическим методом по
результатам измерения массы и объема в реальных условиях процесс достаточно трудоемкий и длительный, перед нами была поставлена задача по
разработке методов неразрушающего контроля лесопродукции в реальных
условиях, в том числе без специальных перемещений отдельных бревен для
процедуры измерений, т.е. контроль в штабеле с доступом к верхнему ряду
и/или к торцам. Наиболее подходящими для решения поставленной задачи
являются акустические методы контроля, а именно низкочастотный ультразвуковой метод и метод свободных колебаний. Аппаратура, реализующая
эти методы безопасна для пользователя и обладает небольшими массогабаритными характеристиками, что позволяет их использовать в оперативных
полевых условиях.
При использовании импульсного ультразвукового метода производится измерение времени t, мкс, прохождения импульсом ультразвуковых
колебаний (УЗК) расстояния между излучателем и приемником и определение скорости распространения продольных упругих колебаний в сплошной
среде C, м/с, которая связана с константами упругости материала. Выражение справедливо для изотропного материала, однако может быть также использовано для оценки действующих (эффективных) модулей материалов с
анизотропией свойств.
При использовании метода свободных колебаний измеряется частота
собственных колебаний (ЧСК) изделия, соответствующая определенному
виду колебаний, затем рассчитывается приведенная скорость распространения акустических волн (стержневая скорость звука, Cl). По предварительно
установленным корреляционным зависимостям определяются необходимые физико-механические свойства. Параметр Сl является весьма информативным, и в ряде случаев достаточно полно характеризует физикомеханические свойства изделий, определяющие их поведение при эксплуатации.
Как показал анализ литературы применительно к древесине, параметры C и Cl связаны с такими важными физико-механическими свойствами как плотность, прочность и упругие свойства. Кроме того, на параметры существенное влияние оказывает влажность контролируемых образцов. Таким образом, эти параметры могут быть использованы в качестве
комплексных параметров, характеризующих физико-механические свой-
9
ства контролируемых изделий, либо – в сочетании с контролем влажности для интегральной оценки плотности.
Для исследований были взяты от 4 до 7 образцов в исходном состоянии свежеспиленных неокоренных бревен (берёза, бук, граб, дуб, ель, кедр,
ольха, осина, пихта, сосна, ясень) каждого вида древесины длиной порядка 5
и 6 м. На всех образцах на различных переделах определялись следующие
характеристики: размеры; масса; плотность; скорость распространения продольных УЗК C в различных направлениях и различных областях образцов
(прибор «Пульсар-1.2»); интегральные значения скорости акустических
волн Cl по результатам измерений частот собственных изгибных и продольных колебаний (прибор «Звук-203М»); влажность образцов с применением
приборов, реализующих кондуктометрический метод контроля (прибор
GMH 3850), а также в соответствии с ГОСТ 16483.7-71и с ГОСТ 16588-91.
Как известно, значения C являются величиной локальноинтегральной и могут служить оценкой анизотропии физико-механических
свойств древесины в различных направлениях: по оси образца, т. е. вдоль
волокон - C║, а также перпендикулярно направлению волокон по диаметру
ствола - C. Для оценки возможной неоднородности свойств образцов измерения проводились в нескольких областях бревен по их длине и в нескольких направлениях по диаметру в указанных плоскостях. При сквозном
прозвучивании измерения проводились в двух взаимно перпендикулярных
направлениях по диаметру в каждой области по длине бревна.
Определение скорости C и C║ распространения УЗК по результатам измерения времени распространения ультразвука выполнялось по одной из стандартных методик с расположением датчиков с конусными
насадками, учитывающими особенности и анизотропию объекта контроля
(Рисунок 1).
При этом схема Рисунок 1а реализовывалась во всех областях по
длине объектов контроля в двух взаимно перпендикулярных направлениях
C, м/с, по формуле
 = ср ⁄ ,
(1)
где Dср – средний диаметр бревна, мм;
Ti – время распространения ультразвука, мкс.
Схема Рисунок 1в реализовывалась аналогично со схемой а. Схема
Рисунок 1б реализовывалась при измерениях на расстоянии 0,25 ÷ 1,0 м от
торцов объектов контроля с расположением датчиков по центру торца и по
диаметру как иллюстрирует Рисунок 2 а, б.
Определение скорости распространения УЗК, Cугл ц  и Cугл д , м/с, при
угловом прозвучивании производилось по формулам
10
угл ц = √(т2 + 2 )/ ,
(2)
где Rт – радиус торца бревна, м;
Li – расстояние от торца, м;
угл д = √(т2 + 2 )/ ,
(3)
где Dт – диаметр торца бревна, м.
При анализе неоднородности свойств образцов учитывалось наличие
у некоторых из них дефектов типа трещин или сучков. Наличие таких дефектов может приводить к изменению скоростей C и Cl как в отдельных
областях, так и по всему образцу в целом, что вносит дополнительные погрешности в результаты измерений.
При исследовании отдельно расположенных образцов (вне штабеля)
акустическим методом по результатам измерений ЧСК и приведенной скорости распространения акустических волн Cl их устанавливали на опорах
или складировали без специальной установки на ровной жесткой поверхности (на поддонах, деревянном настиле и т.п.). Схема расположения объектов
контроля при измерении ЧСК продольных колебаний fпр иллюстрирует Рисунок 3а, при измерении ЧСК изгибных колебаний fизг Рисунок 3б.
При измерениях продольных колебаний использовались только две
нижние частоты fпр1 и fпр2, для которых значение коэффициента формы в
формуле
(4)
 =   ,
где fi – ЧСК определенного вида i;
Fi – коэффициент формы, зависящий от формы и размеров изделия, вида
возбуждаемых колебаний и коэффициента Пуассона;
 = √⁄ ,
(5)
где Е – модуль нормальной упругости;
 - плотность;
определяется выражениями
(6)
пр1 = 1⁄2 ,
1
(7)
пр2 = ⁄ ,
где L длина изделия, м.
Расчет Cl по результатам измерения ЧСК продольных колебаний
длинных стержней, к которым можно отнести стволы дерева, с достаточной
точностью производился по формулам
11
 = 2пр1 ,
 = пр2 .
(8)
(9)
Следует отметить, что при установке отдельно расположенных
объектов контроля fизг5 измерялась устойчиво, а в штабеле измерения fизгi
провести не удалось. Полученные значения соотношения Сl(fизг5) / Сl(fпр2)
подтвердили предположение, о том, что ЧСК изгибных колебаний больше
зависят от свойств в радиальном направлении, а ЧСК продольных колебаний - от свойств древесины вдоль ствола.
В результате обработки и корреляционного анализа статистических
данных по определению плотности классическим методом и акустических
характеристик образцов древесины были построены уравнения регрессии
«Плотность – Скорость» («ρ – C║» и «ρ – Cl») для двух интервалов характеристик, соответствующих характеру изменения плотности в зависимости от
влажности.
На основе установленных регрессионных уравнений были построены градуировочные зависимости, которые легли в основу новых разработанных и аттестованных методик измерений плотности древесины. Построенные градуировочные зависимости, представляют собой кусочнолинейные зависимости. Кусочный характер связан с различным характером
изменения свойств древесины в процессе хранения и сушки и условно соответствует сырой древесине (лесоматериалы), имеющей влажность от 100 до
30 %, и высушенной древесине (пиломатериалы) с влажностью менее 30 %.
2. Градуировочные зависимости «плотность-скорость», установленные с применением метода свободных колебаний и низкочастотного
ультразвукового метода продольно-поперечным способом прозвучивания, обеспечивают достаточно точное определение физически усреднённых значений плотности древесины различных пород в широком диапазоне влажности
В рамках аттестации разработанных методов были рассчитаны коэффициенты корреляции и произведена оценка значимости уравнений регрессии в целом, а также отдельных параметров регрессии и коэффициента корреляции с помощью F-критерия Фишера и t-критерия Стьюдента. Обработка выполнялась в программе Excel в пакете «Анализ данных». Для всех зависимостей установлены доверительные границы суммарной погрешности
результатов измерений. При проверке соответствия показателей точности
результатов измерений и способов обеспечения достоверности измерений
показатели точности определения плотности оценивались, как доверительные границы суммарной погрешности результатов измерений по РМГ 29-
12
99. Принимается, что результаты единичных наблюдений плотности и скорости звука подчиняются нормальному закону распределения; для определения доверительных границ погрешности результатов принята доверительная вероятность Р = 0,95. Доверительные границы случайной погрешности результата измерений вычисляются как t(Р,n)∙Sср, где t (Р,n) - коэффициент Стьюдента для числа степеней свободы (n-1) и доверительной вероятности Р (из таблицы распределения Стьюдента), Sср среднеквадратическое значение случайной составляющей погрешности. В нашем случае
была принята вероятность Р=0,95, число измерений n=10, t = 2,262.
Для наглядности следует частично представить результаты измерений, например, полученные методом сквозного прозвучивания для дуба,
которые отражает Таблица 1.
Таблица 1 – Построение градуировочной зависимости «ρ – C║» для дуба
Номер
Скорость, Плотность,
Вид образца
п/п
м/с
кг/м3
бревно
неокор
957
3740
852
4410
770
4
3560
935
5
3620
935
2910
970
3160
970
700
4630
724
9
5050
724
600
2400
10
4830
724
2
3
6
7
8
доска
Плотность, кг/м3
3430
малогабаритные
образцы
сушка
1
Зависимость плотности древесины (порода дуб) от скорости распространения ультразвуковых колебаний
1200
1100
y = -0,1434x + 1420,1
R² = 0,9366
1000
900
800
3400
4400
5400
Скорость Cǁ∞, м/с
Показатели точности определялись на основании серии результатов
измерений при условиях: одним методом в одинаковых (неизменных) условиях – одинаковое время, один оператор, по одному экземпляру использованных средств измерений, на одном и том же образце объекта исследований (сходимость); одним методом, в разных условиях - разное время, разные операторы, разные экземпляры используемых средств измерений, на
одном образце объекта исследований (воспроизводимость).
Пример установленных уравнений регрессии иллюстрируют Таблица
2 и Таблица 3, где результат измерений плотности в виде абсолютного отклонения записывается как
(10)
ρизм = ρ ± ρ ,
где ρ - ширина доверительного интервала для уравнения регрессии.
13
Относительные границы интервала, в котором с вероятностью 0,95
находится погрешность определения плотности рассчитывались по формуле
 = ±100ρ/ρ % .
(11)
С учетом того, что ширина доверительного интервала является постоянной величиной, а значение  изменяется в среднем в 1,5 раза относительные границы интервала оказываются переменными.
Так, например, для определения плотности дуба по результатам измерения C║ методом сквозного прозвучивания получены следующие зависимости:
для бревен неокоренных, окоренных и пиломатериалов в диапазоне
плотности ρ от 1000 до 820 кг/м3 при абсолютной влажности Wабс от 70 до
40 % (соответственно относительной влажности Wотн от 43 до 28 %)
ρ= -0,143∙C║ +1420 .
Неисключенная систематическая погрешность результата измерений  при (Р = 0,95) с учетом того, что воспроизводимость (близость результатов измерений, полученных разными средствами и разными операторами) в одних и тех же условиях измерения на одном объекте не превышает 0,02 % составила 2,6 %.
Среднеквадратическая погрешность суммы неисключенных систематических погрешностей  при равномерном распределении с учетом
того, что воспроизводимость не превышает 0,02 % равна 1,4 %.
Интервал для среднеквадратического значения случайной составляющей погрешности Sср с учетом того, что воспроизводимость не превышает 0,02% составил (3,5…4,4) %. Ширина доверительного интервала
для уравнения регрессии "ρ – C║" ∆ρдов = 80 кг/м3.
Таким образом, Sср/S ≥ 2,6 % и, соответственно, доверительные
границы суммарной погрешности результатов измерений определения
плотности всегда определяются случайной составляющей погрешности,
то есть шириной доверительного интервала уравнения регрессии ∆ρ = ± 80
кг/м3, что составляет ∆ρ/ρ = 8,0 % при Wабс = 70 % и ∆ρ/ρ = 9,8 % при Wабс =
40 %.
Аналогичные расчеты произведены для каждой исследуемой породы древесины. Для наглядности Таблица 2 отражает установленные градуировочные зависимости «ρ – C║» для берёзы, а Таблица 3 «ρ – Cl» для
14
ели сибирской. При этом C║ вычисляется по формуле (12) в соответствии
со значением Кп из таблицы, установленным для соотношения KD = 80 мм
и KL = 500 мм, а Cl вычисляется с учетом формы, геометрических размеров изделия и выбранного рабочего типа колебаний путем расчета по
формулам (8) и (9).
Таблица 2 – Градуировочная зависимость с уравнением для расчёта плотности лесо- и пиломатериалов по значениям C║
Береза
от 10
до 100
Кп
Уравнение для расчета
плотности - брёвна с корой
и без коры, пиломатериалы, кг/м3
ρ = (-0,554 C║ +3629)
± 40
1
Градуировочная зависимость
Плотность, кг/м3
Порода Интервал
древесины Wабс,%
1100
1050
y = -0,554x + 3629
1000
950
900
850
800
750
700
650
600
4650 4750 4850 4950 5050 5150 5250 5350
Значение Cǁ∞, м/с
Таблица 3 – Градуировочная зависимость с уравнением для расчёта плотности лесо- и пиломатериалов по значениям Cl
от 30
до 100
Уравнение для расчета плотности, кг/м3
ρ = (-0,124 Cl +1025)
± 20
ρ = (-0,267 Cl +1767)
± 35
Градуировочная зависимость
Плотность, кг/м3
Ель сибирская
от 5
до 30
Вид
Брёвна окорённые
Пиломатериалы
Интервал
Wабс,%
Брёвна неокоренные
Порода
древесины
850
y = -0,267x + 1767
800
750
700
650
600
550
500
450
400 y = -0,124x + 1025
350
3700 3950 4200 4450 4700 4950
Скорость Сl, м/с
15
3. Разработанные методы и средства определения плотности
древесины учитывают условия применения и принципиально
отличаются от классической реализации существующих методов возможностью в короткий промежуток времени объективно оценить значения плотности лесо- и пиломатериалов в широком диапазоне влажности без необходимости их перемещения или изготовления
контрольных образцов
По итогам исследования зависимости «плотность древесины – скорость распространения УЗК» был разработан метод определения плотности
древесины по результатам измерений времени распространения УЗК, а также были запатентованы способ и устройство определения плотности древесины продольно-поперечным прозвучиванием.
Согласно разработанному методу вначале производятся измерения
относительной влажности образцов древесины влагомером по ГОСТ 1658891. Затем при помощи графиков зависимости «Wабс-Wотн» определяется значение абсолютной влажности Wабс. Затем измеряется скорость распространения УЗК в лесо- или пиломатериалах с использованием схемы углового
(продольно-поперечного) прозвучивания (в условиях складирования и
транспортирования. Контроль лесоматериалов производится с использованием специальных насадок (концентраторов) на преобразователях, пиломатериалы могут контролироваться как с применением насадок, так и без них.
Преобразователи устанавливаются на контролируемый образец с
углублением насадок в древесину на 1-3 мм в соответствии со схемой которую иллюстрирует Рисунок 4. Максимальная постоянная база (KL = 50 см)
расположена вдоль волокон древесины от торца контролируемого образца.
Минимальная база (KD = 8 см) выставляется на торцевой поверхности в
зависимости от вида материала по направлению диаметра бревна или ширины пиломатериалов. Преобразователи должны быть ориентированы перпендикулярно поверхности контролируемого образца. Размещение преобразователей может фиксироваться с использованием специальных приспособлений. Преобразователи прибора «Пульсар» для измерений устанавливаются в места крепления приспособления жестко и не подлежат перемещению.
Длина базы между преобразователями Li для расчета скорости УЗК
определяется как гипотенуза прямоугольного треугольника с катетами KL и
KD по результатам измерения расстояния между преобразователями. Расчет
скорости УЗК производится по результатам измерения времени и длины
базы измерений. После получения устойчивых показаний прибора производится отсчет времени распространения УЗК и повторяются измерения времени (скорости УЗК) в выбранной области контроля 2-3 раза. Показания
16
при последовательных измерениях не должны отличаться на значение, превышающее 3 мкс при измерении времени распространения УЗК или 50 м/с
при измерении C (контроль сходимости результатов измерений). Область
контроля меняется путем изменения расположения преобразователей поворотом относительно центра торца примерно на 90°. За результат измерения
скорости УЗК принимается среднее значение.
Для получения нормированного значения скорости УЗК (C║) вводится поправочный коэффициент Kп, зависящий от величины размера KD.
Значение C║ рассчитывается по формуле
(12)
C║ = Сугл.изм / Кп ,
где Сугл.изм – значение скорости УЗК,
Кп – коэффициент, учитывающий анизотропию свойств в продольном и
радиальном направлении (C / C║), зависящий от соотношения KD/KL.
Значения Кп экспериментально устанавливаются для конкретных пород древесины. После получения нормированного значения скорости УЗК
определяется значение плотности образца по результатам измерения скорости УЗК при помощи градуировочной зависимости «ρ – C║» для соответствующей породы древесины с учетом результатов измерения влажности.
По итогам исследования зависимости «плотность древесины – скорость распространения акустических волн» был разработан метод определения плотности древесины по результатам измерений частот собственных
колебаний. Определение значения плотности древесины производится при
помощи градуировочной зависимости «ρ – Сl» по результатам измерения
ЧСК и расчета скорости Сl для соответствующей породы древесины с учетом результатов измерения влажности.
Измерения производят на образцах без видимых глазом дефектов.
Отбор образцов лесоматериалов из неразобранного штабеля производится
произвольно исходя из удобства доступа для выполнения измерений или
поштучно при раскатывании из штабеля. Отбор образцов бруса и досок
производится произвольно, поштучно из штабеля. При измерениях отдельно расположенных образцов (вне штабеля) их устанавливают на опорах
(Рисунок 3), или складируют без специальной установки на ровной жесткой
поверхности (на поддонах, деревянном настиле и т.п.). Далее измеряются
длины образцов и определяется их относительная влажность по ГОСТ
16588-91 с последующим установлением значений абсолютной влажности
Wабс при помощи графиков зависимости «Wабс-Wотн». Затем определяется
скорость распространения акустических волн Сl с помощью прибора «Звук203М», для чего в режиме ввода исходных данных устанавливаются время
задержки перед измерением – 5 мс и номер поддиапазона в зависимости от
измеренной длины контролируемого образца.
17
Затем необходимо прибор перевести в режим измерения, поднести
микрофон к торцовой поверхности образца на расстояние порядка 10 мм и
не касаясь ее, в области (Рисунок 3), произвести измерение скорости Сl (частоты типа fпр1 или fпр2), для чего нанести ударником (молоток весом не менее 0,5 кг) удар в центре торцовой грани бревна (образца пиломатериала) в
направлении, перпендикулярном поверхности (Рисунок 5), и зафиксировать
показание прибора. Допускается наносить удар и производить прием колебаний микрофоном как с одной, так и с противоположной стороны (торца).
Далее необходимо повторить удар.
Показания при двух последовательных ударах должны отличаться на
значения, не превышающие в режиме индикации скорости Сl ± 50 м/с (контроль сходимости результатов измерений). При этом реальное время (число
импульсов) измерения в процентах от заданного должно быть не менее 80
% (контроль достоверности результатов измерений). В случае выполнения
условий, за результат измерения принимается величина отсчета (наблюдения) по прибору, полученная при втором ударе. Если различия превышают
указанные значения, или реальное время (число импульсов) измерения в
процентах от заданного менее 80 %, необходимо изменить номер поддиапазона и повторить операции. При повторном получении различий, превышающих указанные, необходимо проверить контролируемый образец на
наличие дефектов, например, трещин, гнили и т.д., а также произвести автотестирование прибора путем его выключения и повторного включения.
Частотный коэффициент формы Fi в зависимости от размера (длины)
и рабочего вида колебаний рассчитывается по формулам (6) и (7). Следовательно, определение приведенной скорости распространения акустических
волн Сl с учетом формы, геометрических размеров изделия и выбранного
рабочего типа колебаний производится путем расчета по формулам (8) и (9).
Далее производится запись результатов измерений в протокол (память прибора) и определяется значение плотности образца при помощи градуировочной зависимости «ρ – Сl». Корреляционная связь (градуировочная зависимость) «ρ – Сl» устанавливается предварительно, на основании экспериментальных данных для каждого вида древесины в определенном диапазоне
значений массовой доли влаги. Результаты определения плотности всех образцов из выборки, усредняются и присваиваются всей представленной партии при условии, если размах полученных на выборке значений плотности
не превышает 20 % (контроль однородности партий). В случае превышения
указанного размаха принимается решение либо о неоднородности партии,
либо о необходимости повторных измерений с увеличением выборки.
18
ЗАКЛЮЧЕНИЕ
В диссертации, представляющей собой законченную научноквалификационную работу, на базе выполненных теоретических и экспериментальных исследований была решена актуальная научно-практическая
задача. Обоснована техническая возможность, применения низкочастотных
акустических методов контроля для достаточно точной оценки плотности
лесопродукции, путем применения низкочастотного измерителя времени
ультразвуковых колебаний и измерителя частот собственных колебаний.
Основные научные и практические результаты исследований:
1. Исследованы зависимости параметров акустических колебаний,
распространяющихся в древесине, от её плотности, породы и влажности, с
учётом времени и места заготовки древесины;
2. Установлены градуировочные зависимости «плотность-скорость»
для различных пород древесины (Береза, Бук, Граб, Дуб, Ель европейская,
Ель сибирская, Осина, Ольха, Орех, Пихта, Сосна, Сосна кедровая, Ясень) в
широком диапазоне влажности;
3. Теоретически обоснована и экспериментально подтверждена возможность определения плотности древесины в условиях таможенного контроля с применением низкочастотного ультразвукового метода и метода
свободных колебаний;
4. Разработан и запатентован способ и устройство определения
плотности древесины продольно-поперечным прозвучиванием – патент РФ
2 449 265, МПК G01N 29/07. «Способ и устройство определения плотности
древесины»;
5. Разработаны методы и средства таможенного контроля плотности
различных пород древесины в широком диапазоне влажности;
6. Разработано прикладное программное обеспечение в основу
функционирования которого положены результаты практического применения разработанных методов и средств, а также даны рекомендации по их
применению в таможенной практике.
По теме диссертации опубликованы следующие работы в изданиях, рекомендованных ВАК:
1. Симоненко, А.А. Определение плотности древесины с применением
низкочастотного ультразвукового метода [Текст] / А.А. Симоненко, Е.З. Коварская // Научно-практический журнал «Наука и бизнес: пути развития», №11(29).
– М., 2013. – С. 44-50.
2. Симоненко, А.А. Построение градуировочных зависимостей для
различных пород древесины с применением низкочастотного акустического
метода свободных колебаний [Текст] // Научно-технический вестник Поволжья,
19
№6 2013г. – Казань: Научно-технический вестник Поволжья, 2013. – С. 422426.
3. Симоненко, А.А. Определение плотности древесины с применением
низкочастотного акустического метода свободных колебаний [Текст] / А.А.
Симоненко, Е.З. Коварская // Научно-технический журнал РОНКТД «Контроль. Диагностика», №2 (188). – М. – 2014. – С. 43-49.
4. Пат. 2 449 265 РФ, МПК G01N 29/07. Способ и устройство определения плотности древесины [Текст] / Коварская Е.З., Черноглазов В.С., Воронин
А.А., Смирнова Е.В., Симоненко А.А., Московенко И.Б.; заявитель и патентообладатель Санкт-Петербург, ЗАО «Научно-производственный центр «Инновационная техника и технологии». – № 2010147336/28; заявл. 11.11.2010; опубл.
27.04.2012, Бюл. № 12. – 8 с.: ил.
В других изданиях:
5. Симоненко, А.А. Проблемы таможенного контроля при экспорте леса и возможности их решения [Текст] / А.А. Симоненко, И.В. Павлов // Труды
Х Международной научно-практической конференции молодых ученых, студентов и аспирантов: «Анализ и прогнозирование систем управления» II ч.,
РАН Дом ученых в СПб.: СЗТУ, 2009. – С. 289-293.
6. Симоненко, А.А. Актуальные вопросы определения плотности древесины при таможенном контроле [Текст] / А.А. Симоненко, И.В. Павлов //
Анализ и прогнозирование систем управления: труды XIII Международной
научно-практической конференции молодых ученых, студентов и аспирантов –
СПб.: Петербургский гос. ун-т путей сообщения, 2012. – С. 566-572.
7. Симоненко, А.А. К вопросу о таможенном контроле плотности древесины [Текст] / А.А. Симоненко, Е.З. Коварская // Вестник БелорусскоРоссийского университета. – 2013. – № 3 (40) – С. 67-73.
8. Симоненко, А.А. Плотность древесины - обзор методов и новые способы определения [Текст] // Современные техника и технологии: сборник трудов XIХ Международной научно-практической конференции студентов, аспирантов и молодых ученых. В 3 т. Т. 1 / Томский политехнический университет.
– Томск: Изд-во Томского политехнического университета, 2013. – С. 194-195.
9. Симоненко, А.А. Метрологические трудности использования современных приборов и методов таможенного контроля лесоматериалов [Текст] //
Инновационные системы планирования и управления на транспорте и в машиностроении: сборник трудов II международной научно-практической конференции. Том II // Под ред. Е.И. Пряхина. Гл. ред. И.В. Павлов. Отв. Ред. Д.Ю.
Тимофеев, Т.А. Менухова, А.В. Терентьев – СПб: Национальный минеральносырьевой университет «Горный», 2014. – С. 188-191.
20
а)
б)
в)
Рисунок 1 – Расположение преобразователей прибора «Пульсар-1.2» относительно
объекта контроля при различных стандартных способах прозвучивания – а) сквозное; б) угловое; в) поверхностное
а)
б)
Рисунок 2 – Расположение преобразователей прибора «Пульсар-1.2» при продольно-поперечном прозвучивании
а) fпр
Рисунок 4 – Схема расположения преобразователей при продольно-поперечном прозвучивании
б) fизг
Обозначение: 1- отдельно расположенный объект контроля, 2 – микрофон,
3 – специальные опоры шириной не более 100 мм и длиной не менее
100 мм (деревянные или покрытые слоем резины толщиной 100-120 мм),
расположенные на расстоянии ≈ 0,25L от торцов, - стрелкой показано
место и направление удара.
Рисунок 3 – Схема расположения объекта контроля при измерении ЧСК продольных и изгибных колебаний
Рисунок 5 – Измерение ЧСК бревна прибором «Звук-203М»
Документ
Категория
Без категории
Просмотров
18
Размер файла
1 335 Кб
Теги
древесины, метод, контроля, плотности, средств, таможенного
1/--страниц
Пожаловаться на содержимое документа