ПРОЧИЕ ВОПРОСЫ СТРОИТЕЛЬСТВА И АРХИТЕКТУРЫ / OTHER QUESTIONS RELATED TO
CONSTRUCTION AND ARCHITECTURE
DOI: https://doi.Org/10.18454/mca.2022.27.2
ПОСТРОЕНИЕ И ПРОВЕРКА АДЕКВАТНОСТИ УРАВНЕНИЯ РЕГРЕССИИ КОЭФФИЦИЕНТА ТЕПЛОПРОВОДНОСТИ ОТ СООТНОШЕНИЯ ИСХОДНЫХ КОМПОНЕНТОВ КОМПОЗИТНОГО МАТЕРИАЛА НА ОСНОВЕ ГИПСОВОГО ВЯЖУЩЕГО И ДРЕВЕСНЫХ ОПИЛОК
Научная статья
Ковалев Н.С.1, Горохов Т.И.2, Ерофеев А.В.3 *
1 2 3 Тамбовский государственный технический университет, Тамбов, Россия
* Корреспондирующий автор (av.erofeev[at]yandex.ru)
Аннотация
В настоящее время в строительстве используется большой спектр материалов, большинство из которых являются композитными, т. е. состоящими из двух и более компонентов. В связи с этим встает вопрос о подборе оптимального состава композитного материала, который может быть выполнен с использованием математического планирования эксперимента. В работе для определения оптимального состава композитного материала на основе гипсового вяжущего с добавлением древесных опилок по критерию теплопроводности было выполнено математическое планирование эксперимента. Для измерения коэффициента теплопроводности и определения теплового сопротивления при стационарном тепловом режиме использовать измеритель теплопроводности ИТП-МГ4. Также в работе получено уравнение регрессии, построена поверхность отклика и выполнена проверка адекватности полученных коэффициентов регрессии.
Ключевые слова: теплопроводность, композит, экспериментальные испытания, математическое планирование, уравнения регрессии, адекватность, величина отклика.
CONSTRUCTION AND VERIFICATION OF THE ADEQUACY OF REGRESSION EQUATION OF THE THERMAL CONDUCTION COEFFICIENT FROM THE RATIO OF THE INITIAL COMPONENTS OF A COMPOSITE MATERIAL BASED ON GYPSUM BINDER AND SAWDUST
Research article
Kovalev N.S.1, Gorokhov T.I.2, Erofeev A.V.3 *
1 2 3 Tambov State Technical University, Tambov, Russia
* Corresponding author (av.erofeev[at]yandex.ru)
Abstract
Today, a wide range of materials are used in construction, most of which are composite, i.e., they consist of two or more components. In this regard, the question arises about the selection of the optimal composition of composite material, which can be performed with mathematical planning of the experiment. In order to determine the optimal composition of a composite material based on gypsum binder with a supplement of sawdust in accordance with the criterion of thermal conduction, mathematical planning of the experiment was conducted. In order to measure the coefficient of thermal conduction and to determine thermal resistance in a steady temperature condition, an ITP-MG4 thermal conductivity meter must be used. The regression equation is also obtained, the response surface is constructed and the adequacy of the obtained regression coefficients is verified.
Keywords: thermal conduction, composite, experimental tests, mathematical planning, regression equations, adequacy, response value.
Введение
В современном строительстве с каждым годом все чаще используют различные виды отделочных материалов [1]. Как правило, они представляют собой композитные материалы, под которыми понимаются материалы, изготавливаемые из двух и более компонентов, имеющих различные физические и химические свойства, но при соединении дающие новый материал с новыми свойствами. Одним из новых отделочных материалов является материал на основе гипсового вяжущего с добавлением древесных опилок (рис.1). Материал предназначен для внутренней отделки стен помещений, в том числе и наружних.
Рис. 1 - Отделочный материал на основе гипсового вяжущего с добавлением древесных опилок
Разработка любого композитного материала сопряжена с определением оптимального состава, который дает наилучшие характеристики. Немаловажную роль для рассматриваемого материала, в случае, когда он будет применяться для отделки наружных стен, является такой показатель как коэффициент теплопроводности [2]. Использование отделочного материала с более высоким коэффициентом термического сопротивления повышает энергетическую эффективность здания. Кроме того, применение данного материала при отделке помещений объектов капитального ремонта позволит сократить расходы на утепление наружных стен. Именно поэтому представляется актуальной задачей определение зависимости изменения коэффициента теплопроводности отделочного материала на основе гипсового вяжущего с добавлением древесных опилок от соотношения исходных компонентов.
Вопросы по рациональному использованию древесных отходов для получения материалов и изделий теплоизоляционного назначения, характеристике этих материалов и применяемому технологическому оборудованию рассматривались в работах ученых Абраменкова Н.И., Бужевича Г.А., Валуевой Е.Ф., Егорова А.Д., Запруднова В.И., Косимова О.Б., Кучерявого В.И., Мельниковой Л.В., Мурзина B.C., Прусса Б.Н., Разинькова Е.М., Рахимова Р.З., Савина В.И., Удербаева С.С., Уитнея Ч., Федорова C.B., Хакимова Ш.А., Ходжаева Ш.А., Чепелева С.Р., Шамаева В.А., Щербакова A.C. и др. Однако, исследования наиболее близкого по своей структуре материала было отражено при разработке гипсостружечных плит. Однако, некоторые особенности технологических режимов изготовления древесно-гипсовых плит и их эксплуатационные показатели сейчас изучены недостаточно.
Методы и принципы исследования
Для измерения коэффициента теплопроводности и определения теплового сопротивления при стационарном тепловом режиме по ГОСТ 7076-99 [3] использовался измеритель теплопроводности ИТП-МГ4 (рис. 2), принцип работы которого основан на создании стационарного теплового потока, проходящего через плоский образец определенной толщины и направленного перпендикулярно к лицевым граням образца. Конструктивно прибор состоит из нагревателя - преобразователя, выполненного в виде стационарной установки и электронного блока [4]. Стационарная установка прибора в свою очередь состоит из блока управления холодильником и нагревателем, выполненного на элементах Пельтье (охлаждение осуществляется вентилятором), устройства преобразования первичных сигналов датчиков, платиновых датчиков температуры, тепломера, а также источника питания. Установка в верхней части имеет прижимной винт, снабженный динамометрическим устройством с трещоткой для создания постоянного усилия прижатия испытываемого образца и отсчетным устройством для измерения толщины образца.
Рис. 2 - Измеритель теплопроводности ИТП-МГ4
Определение коэффициента теплопроводности и теплового сопротивления проходило следующим образом. Измеритель теплопроводности ИТП-МГ4 (рис. 2) подключается к сети переменного тока 220 В через сетевой шнур. После подключения ИТП-МГ4 к сети установка открывается и в нее помещается образец, толщина которого предварительно измеряется штангенциркулем. Открытие установки производится путем ослабления прижимного микрометрического винта с последующим поворотом против часовой стрелки эксцентрикового замка (происходит освобождение поводка коромысла) и поднятием поводка. Для окончательного открытия установки подвижная Г-образная стенка установки отводится на 90°, а коромысло с закрепленной на нем плитой нагревателя поднимается. После помещения образца в установку коромысло опускается. При этом зазор между образцом и плитой нагревателя устанавливается путем вращения микрометрического винта в интервале от 2 до 5 мм. Далее поводок поднимается, а Г-образная стенка закрывается. Закрепление поводка осуществляется путем его опускания в паз эксцентрикового замка с последующим его поворотом по часовой стрелке.
После включения питания электронного блока на дисплей вводятся фактическая толщина образца и устанавливаются требуемые температуры холодильника и нагревателя, которые выбираются в соответствии с рекомендациями руководства по эксплуатации измерителя теплопроводности ИТП-МГ4 в зависимости от прогнозируемой теплопроводности материала. В дальнейшем программное устройство прибора устанавливает на поверхностях образца заданные температуры и поддерживает их до тех пор, пока тепловой поток, проходящий через образец, не стабилизируется. Наблюдение за тепловым потоком осуществляется автоматически. Вычисление значений коэффициента теплопроводности и теплового сопротивления также происходит автоматически через 30 минут после начала испытания.
Основные результаты
Для определения оптимального состава композитного материала по критерию теплопроводности было выполнено математическое планирование эксперимента [5], целью которого является нахождение уравнения регрессии зависимости коэффициента теплопроводности от соотношения исходных компонентов, а также проверка ее адекватности. В ходе проведения математического планирования эксперимента [6] были определены четыре точки плана эксперимента, каждой из которых соответствует определенные граничные соотношение компонентов: гипсового вяжущего, воды и древесных опилок. Процентные содержания исходных компонентов с учетом их кодировок (+1 и -1) были занесены в таблицу 1, в которую также были занесены данные экспериментальных исследований теплопроводности опытных образцов материала.
Для повышения точности данных измерения проводились согласно техническим рекомендациям, которые шли с прибором: каждый образец испытывался не менее трёх раз для того, чтобы проверить адекватность эксперимента [7], которая показывает способность предсказывать результаты эксперимента в некой области с требуемой точностью. После чего полученные результаты усреднялись.
Таблица 1 - Матрица планирования
№ опыта План в кодированных переменных План в отношениях План в исходных переменных Отклик у, теплопроводность X, Вт/м2К Теплопроводность по формуле Некрасова Хн, Вт/м2К
51 52 21 22 Х], % Х2, % Хз, % У1 У2 Уз У
1 -1 -1 0,78 0,22 38,89 50,00 11,11 0,215 0,210 0,212 0,212 0,40
2 -1 +1 0,78 0,44 35,00 45,00 20,00 0,187 0,185 0,172 0,181 0,32
3 +1 -1 1,00 0,22 45,00 45,00 10,00 0,229 0,224 0,227 0,227 0,40
4 +1 +1 1,00 0,44 40,91 40,91 18,18 0,236 0,235 0,237 0,236 0,35
Полученные экспериментальные данные сопоставлялись с теоретическими результатами теплопроводности, полученные по формуле Некрасова:
Я = 1,163 -70,0196 + 0,22^2 - 0,16 (1)
^Рт Рш , 3
= — =- - относительная плотность материала, кг/м3, здесь
рш 1000 р
рт - плотность материала, кг/м3; = 1000 кг/м3 - плотность воды.
Анализируя результаты, представленные в таблице 2, стоит отметить, что коэффициент теплопроводности,
полученный по формуле Некрасова, почти в 2 раза выше, чем экспериментально полученный. По-видимому, это связано
с тем, что формула Некрасова является эмпирической формулой и позволяет определять только ориентировочные
значения коэффициентов теплопроводности, что в свою очередь только доказывает необходимости нахождения
уравнения регрессии, позволяющей на основе экспериментальных данных определять истинные значения
коэффициента теплопроводности материала в зависимости от соотношения исходных компонентов.
Уравнение регрессии в общем виде может быть представлено следующим образом:
Хз
у = Ьп + й-,--+ Ь2 —
0 1Х2 2Х2 (2)
где Ьоо, Ь], Ь2 - коэффициенты уравнения регрессии.
В качестве аргументов Х1, Х2, Хз приняты соответственно вода, гипс марки Г-16 (вяжущее вещество), опилки (мелкий заполнитель).
Тогда уравнение регрессии с учетом условий кодирования может быть представлено в виде:
у = Й0 + ¿1^1 + Й2^2, (3)
где
(0) (0) ^-Ж0
^ = - ^ ) = Х2 УХ2У д = ^2 - ^2 ) = Х2 УХ2У 1 Р21 Р21 , 2 Рг2 Рг2 (4)
^А0 „ „
где ^—^ , ^—^ - отношения относительных содержаний компонентов в исходной точке плана, а р21, р22 -
интервалы варьирования новой переменных г1 и г2 , равные [г1(0) — 21(н)], ^(в) — ^2(0)] соответственно; /хЛ(в) (хЛ(в) „ /хЛ(н) /хЛ(н) „
^—^ , ^—^ - верхний уровень отношения компонентов, ^—^ , ^—^ - нижний уровень отношения компонентов,
На основании данных таблицы 1 были найдены коэффициенты уравнения регрессии Ьо, Ь1 и Ь2: г, У1+У2+У3+У4 0,212+0,181+0,227+0,236 ппЛЛ
Ь0=-=-= 0,214 /с\
0 4 4 (5)
и -У1-У2+У3+У4 -0,212-0,181+0,227+0,236
=- = - = 0,017 (6)
1 4 4 (6)
и -У1+У2-У3+У4 -0,212+0,181-0,227+0,236 -
От =- = - = —0,005 пл
2 4 4 (7)
у = 0,214 + 0,0175: - 0,00552
После чего были найдены и р22:
Р21 = (0) - г^н) = 0,89 - 0,78 = 0,11
Рг2 = ^(в) - ^2(0) = 0,44 - 0,33 = 0,11
где р21, р22 - интервалы варьирования новой переменных. Подставив полученные значения в формулу (4), получим:
(8)
(9) (10)
^ - 0,89
*1 = Х2
0,11
ХЛ
= 9,09 — - 8,09,52 =^
г 0,зз х
Хп лз
2-= 9,09 — - 3,00
0,11 , х2 ,
Следовательно, уравнение регрессии будет иметь следующий вид:
у = 0,091 + 0,155^2 - 0,045*2
(11)
(12)
Далее проверяем воспроизводимость эксперимента. Результаты определения однородности оценок дисперсий представлены в таблице 2.
2
Таблица 2 - Однородность оценок дисперсий
^ =( У - У1)2 ^ = ( У - У 2 ! ^ =(У - Уз )2 О
0,000007 0,000005 0,000000 0,561
0,000032 0,000013 0,000087 0,657
0,000005 0,000007 0,000000 0,561
0,000000 0,000001 0,000001 0,500
Критическое значение Окр, для У1вос = т-1 = 2 и У2вос= N = 4 и уровня значимости ^вос = 5 %, составляет 0,7679. В Расчетная величина О для экспериментальных данных не превышает критического значения Окр, следовательно, гипотеза об однородности выборочных дисперсий отвечает результатам наблюдений.
*в2ос{у} = 0,00002
Проверим гипотезу о значимости коэффициентов уравнения регрессии Ь с помощью критерия Стьюдента:
5{Ь} =
(13)
N
0,00002
3 • 4
= 0,00129
£ = ■
!Ы
(14)
(15)
0,00129
По критерию Стьюдента при числе степеней свободы Узн= Щт-1) = 8 и при заданном уровне значимости qзн (5 %), 1кр = 2,3060. Коэффициенты регрессии представлены в таблице 3.
Таблица^ 3 - Коэффициенты регрессии
Ь t
Ьо 0,214 165,956
Ьг 0,017 13,372
Ь2 -0,005 4,199
Найденная величина параметра t превышает значение ^ для оценок Ьо, Ьг и Ь2, поэтому оценка всех коэффициентов является значимыми.
Проверим гипотезу об адекватности математического описания при Vlад=N-d = 1, V2ад=Vзн=N(m-1) = 8 и при заданном уровне значимости qад (5 %).
Таблица 4 - Оценка дисперсии воспроизводимости отклика
У ¿1 ¿2 y = 0,214 + 0,01751 - 0,00552 ly -yl |у-у|2
0,212 -1 -1 0,202 0,010 0,00011
0,181 -1 1 0,192 -0,011 0,00011
0,227 1 -1 0,236 -0,009 0,00009
0,236 1 1 0,226 0,010 0,00010
Среднее значение 0,00010
F = 0,00129/0,0001 = 12,9000 (16)
Вычисленное по результатам наблюдений значение критерия Фишера
¥ = 12,90. Критическое значение ¥кр = 239. Эмпирическое значение критерия ¥ меньше критического Fкр, следовательно, гипотеза об адекватности не отвергается, и полученное нами математическое описание является адекватным. Из гипотезы о значимости следует что все коэффициенты Ъ1, Ъ2, Ъ3 значимы. 52
После проверки, было выполнено построение поверхности отклика (рисунок 3) по уравнению отклика (формула 12) и по данным таблицы 1.
■ 0,19-0,20 ■ 0,20-0,21 ■ 0,21-0,22 "0,22-0,23 ■ 0,23-0,24
Рис. 3 - Поверхности отклика теплопроводности
Анализ полученного уравнения регрессии и поверхности отклика (рис. 3) позволяет сделать вывод о том, что оптимальный состав рассматриваемого композита с точки зрения теплозащитных характеристик материала соответствует кодировку (-1;+1), т. е. вода - 35%; гипсовое вяжущее - 45%; древесные опилки - 20%.
Заключение
На основе использования математического планирования эксперимента была получена уравнение регрессии, позволяющее установить математическую зависимость между коэффициентом теплопроводности отделочного материала на основе гипсового вяжущего с добавлением древесных опилок от соотношения исходных компонентов материала. В кодированных значениях оно описывается формулой 4. Проведение данного экспериментального испытания позволило вывести уравнения регрессии, используя математическое планирование. Проверка коэффициентов уравнения регрессии показала их адекватность. Анализ уравнения регрессии и поверхности отклика позволил получить оптимальный состав рассматриваемого композита с точки зрения теплозащитных характеристик.
Conflict of Interest Конфликт интересов
None declared. Не указан.
Список литературы / References
1. Власов О.Е. Строительная теплофизика. Состояние и перспективы развития / О.Е. Власов. - М.: Гостстроиздат, 1961. - 290 с.
2. Горчаков Г.И. Прогнозирование теплопроводности композиционных материалов различного строения / Г.И. Горчаков, И.И. Лифанов, А.А. Багаутдинов и др. // Строительные материалы. - 1992. - № 4. - С. 27-28.
3. ГОСТ 7076-99 Межгосударственный стандарт «Материалы и изделия строительные» - Метод определения теплопроводности и термического сопротивления при стационарном тепловом режиме ОКС 27.220 ОКСТУ 5709
4. Общество с ограниченной ответственностью «Специальное конструкторское бюро Стройприбор» ОКП 42 7646 «Измеритель теплопроводности ИТП-МГ4» Руководство по эксплуатации Э 12.102.010 РЭ - Челябинск
5. Макаричев Ю.А. М 30 Методы планирование эксперимента и обработки данных: учеб. пособие / Макаричев Ю.А., Иванников Ю.Н. - Самара: Самар. гос. техн. ун-т, 2016. - 131 с.: ил.
6. Кузнецова Е.В. Э413 Математическое планирование эксперимента: Учебно-методическое пособие для студентов очного и заочного обучения специальностей «Технология обработки металлов давлением», «Динамика и прочность машин», «Компьютерная механика», «Компьютерная биомеханика» / Е.В. Кузнецова - Пермь: Перм. гос. техн. ун-т, 2011. - 35 с.
7. Соколовская И.Ю. Полный факторный эксперимент / И.Ю. Соколовская // Методические указания для самостоятельной работы студентов. - Новосибирск: НГАВТ, 2010. - 36 с.
Список литературы на английском языке / References in English
1. Vlasov O.E. Stroitel'naya teplofizika. Sostoyanie i perspektivy razvitiya [Construction thermophysics. State and prospects of development] / O.E. Vlasov. - M.: Goststroizdat, 1961. - 290 p. [in Russian]
2. Gorchakov G.I. Prognozirovanie teploprovodnosti kompozicionnyh materialov razlichnogo stroeniya [Prediction of thermal conduction of composite materials of various structures] / G.I. Gorchakov, I.I. Lifanov, A.A. Bagautdinov et al. // Stroitel'nye materialy [Building materials]. - 1992. - No. 4. - pp. 27-28. [in Russian]
3. GOST 7076-99 Mezhgosudarstvennyj standart "Materialy i izdeliya stroitel'nye" - Metod opredeleniya teploprovodnosti i termicheskogo soprotivleniya pri stacionarnom teplovom rezhime OKS [Interstate standard "Materials and construction products" - Method for determining thermal conduction and thermal resistance in steady temperature condition CCO] 27.220 OKSTU 5709 [in Russian]
4. Limited Liability Company "Special'noe konstruktorskoe byuro Strojpribor" RCP 42 7646 "Thermal conductivity meter ITP-MG4" Operating Manual E 12.102.010 RE - Chelyabinsk [in Russian]
5. Makarichev Yu.A. M30 Metody planirovaniya eksperimenta i obrabotki dannyh [Methods of experiment planning and data processing]: study guide / Yu.A. Makarichev, Yu.N. Ivannikov - Samara: Samara State Technical University, 2016. -131 p.: ill. [in Russian]
6. Kuznetsova E.V. E413 Matematicheskoe planirovanie eksperimenta: Uchebno-metodicheskoe posobie dlya studentov ochnogo i zaochnogo obucheniya special'nostej "Tekhnologiya obrabotki metallov davleniem", "Dinamika i prochnost' mashin", "Komp'yuternaya mekhanika", "Komp'yuternaya biomekhanika" [Mathematical planning of the experiment: An educational and methodical manual for full-time and part-time students of the specialties "Technology of metal processing by pressure", "Dynamics and strength of machines", "Computer mechanics", "Computer biomechanics"] / E.V. Kuznetsova. - Perm: Perm State Technical University. 2011. - 35 p. [in Russian]
7. Sokolovskaya I.Yu. Polnyj faktornyj eksperiment [Complete factorial experiment] / I.Yu. Sokolovskaya // Metodicheskie ukazaniya dlya samostoyatel'noj raboty studentov [Methodological guidelines for independent work of students]. - Novosibirsk: NGAVT, 2010. - 36 p. [in Russian]_