Научная статья на тему 'Определение сопротивления теплопередаче ограждающих конструкций'

Определение сопротивления теплопередаче ограждающих конструкций Текст научной статьи по специальности «Физика»

CC BY
441
31
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ОГРАЖДАЮЩИЕ КОНСТРУКЦИИ / СОПРОТИВЛЕНИЕ ТЕПЛОПЕРЕДАЧЕ / НЕСТАЦИОНАРНЫЙ РЕЖИМ / АРГУМЕНТЫ FO И BI / НЕОГРАНИЧЕННАЯ ПЛАСТИНА / ГРАНИЧНЫЕ УСЛОВИЯ

Аннотация научной статьи по физике, автор научной работы — Хрусталев Б.М., Сизов В.Д.

проведение мероприятий по увеличению термического сопротивления ограждающих конструкций требует определения их теплофизических характеристик, в основу методов определения которых положено решение задач теплопроводности, устанавливающее связь между пространственно-временными изменениями температуры под действием источника теплоты. в данной работе используется решение задачи при нестационарном нагреве ограждающей конструкции в виде неограниченной пластины при граничных условиях iii рода. по известным соотношениям и графикам определяются изменения температуры поверхности от времени прогрева, термического сопротивления конструкции, аргументов fo и bi, т. е. начальных и граничных условий. полученные графические зависимости показывают, что температура поверхности зависит от термического сопротивления, а температура на противоположной поверхности за время теплового воздействия практически не изменяется за время  = 5 ч. таким образом, если наружная температура воздуха изменилась, то по скорости изменения температуры поверхности или относительной температуры  можно определить теплофизические характеристики, решая обратную задачу теплопроводности, используя преобразованное соотношение для определения r как функцию r = f(, ). если использовать построенные графические зависимости r = f(, ) при различных коэффициентах теплоотдачи, то по измеренным температурам в различные промежутки времени можно определить термическое сопротивление в эти же промежутки времени и по их среднему значению искомое сопротивление теплопередаче r. полученные нами расчетные аналитические соотношения и графические зависимости показывают адекватность проведенным натурным измерениям, если выбирать участки с однородным температурным полем и температурной историей, и их можно использовать при определении сопротивления теплопередаче ограждающей конструкции в виде неограниченной пластины при граничных условиях iii рода.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Fulfillment of the activities aimed to an increase of the thermal resistance of enclosing structures requires the determination of their thermophysical characteristics with the use of the determination method based on the solution of problems of heat conduction, establishing the connection between the spatial and temporal temperature changes under the effect of heat source. This work uses the solution of the problem under nonstationary heating of the enclosing structure in the form of unrestricted plate with boundary conditions of the III kind. According to the known relations and graphs alterations in surface temperature depending on warm-up time, on thermal resistance of constructions and on arguments of Fo and Bi, i. e. initial and boundary conditions are determined. The graphic dependencies that have been obtained show that the surface temperature depends on the thermal resistance, while the temperature at the opposite surface during heat exposure remains practically unchanged during  = 5 h. Thus, if the outside air temperature is altered, then the rate of change of surface temperature or relative temperature  make it possible to determine the thermophysical characteristics by solving the inverse problem of thermal conductivity with the use of the converted ratio to determine R as a function R = f(, ). If the constructed graphic dependencies R = f(, ) are used at different heat transfer coefficients, then according to the measured temperatures at different time intervals it is possible to determine thermal resistance in the same time intervals and, according to their average value, determine the required resistance to heat transfer R. The estimated ratio of analytical and graphic dependencies that we have obtained demonstrate the adequacy of the conducted full-scale measurements, if the areas with homogeneous temperature field and temperature history are chosen, and they can be used in determining the heat resistance of the enclosing structure in the form of unrestricted plate with boundary conditions of the III kind.

Текст научной работы на тему «Определение сопротивления теплопередаче ограждающих конструкций»

Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. Т. 61, № 1 (2018), с. 47-59 Energetika. Proc. CIS Higher Educ. Inst. and Power Eng. Assoc. V. 61, No 1 (2018), pp. 47-59 47

DOI: 10.21122/1029-7448-2018-61-1-47-59 УДК 697.112.3:620.179.11

Определение сопротивления теплопередаче ограждающих конструкций

Б. М. Хрусталев1*, В. Д. Сизов1*

^Белорусский национальный технический университет (Минск, Республика Беларусь)

© Белорусский национальный технический университет, 2018 Belarasian National Technical University, 2018

Реферат. Проведение мероприятий по увеличению термического сопротивления ограждающих конструкций требует определения их теплофизических характеристик, в основу методов определения которых положено решение задач теплопроводности, устанавливающее связь между пространственно-временными изменениями температуры под действием источника теплоты. В данной работе используется решение задачи при нестационарном нагреве ограждающей конструкции в виде неограниченной пластины при граничных условиях III рода. По известным соотношениям и графикам определяются изменения температуры поверхности от времени прогрева, термического сопротивления конструкции, аргументов Fo и Bi, т. е. начальных и граничных условий. Полученные графические зависимости показывают, что температура поверхности зависит от термического сопротивления, а температура на противоположной поверхности за время теплового воздействия практически не изменяется за время т = 5 ч. Таким образом, если наружная температура воздуха изменилась, то по скорости изменения температуры поверхности или относительной температуры 0 можно определить теплофизические характеристики, решая обратную задачу теплопроводности, используя преобразованное соотношение для определения R как функцию R = f(0, т). Если использовать построенные графические зависимости R = f(0, т) при различных коэффициентах теплоотдачи, то по измеренным температурам в различные промежутки времени можно определить термическое сопротивление в эти же промежутки времени и по их среднему значению - искомое сопротивление теплопередаче R. Полученные нами расчетные аналитические соотношения и графические зависимости показывают адекватность проведенным натурным измерениям, если выбирать участки с однородным температурным полем и температурной историей, и их можно использовать при определении сопротивления теплопередаче ограждающей конструкции в виде неограниченной пластины при граничных условиях III рода.

Ключевые слова: ограждающие конструкции, сопротивление теплопередаче, нестационарный режим, аргументы Fo и Bi, неограниченная пластина, граничные условия

Для цитирования: Хрусталев, Б. М. Определение сопротивления теплопередаче ограждающих конструкций / Б. М. Хрусталев, В. Д. Сизов // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2018. Т. 61, № 1. С. 47-59. DOI: 10.21122/1029-7448-201861-1-47-59

Адрес для переписки

Хрусталев Борис Михайлович Белорусский национальный технический просп. Независимости, 150, 220013, г. Минск, Республика Беларусь Тел.: +375 17 265-96-56 tgv_fes@bntu.by

Address for correspondence

Khroustalev Boris M.

Belarusian National Technical University

150 Nezavisimosty Ave.,

220013, Minsk, Republic of Belarus

Tel.: +375 17 265-96-56

tgv_fes@bntu.by

Determining Heat Transmission Resistance of Enclosing Structures

B. M. Khroustalev1), V. D. Sizov1)

'-Belarasian National Technical University (Minsk, Republic of Belarus)

Abstract. Fulfillment of the activities aimed to an increase of the thermal resistance of enclosing structures requires the determination of their thermophysical characteristics with the use of the determination method based on the solution of problems of heat conduction, establishing the connection between the spatial and temporal temperature changes under the effect of heat source. This work uses the solution of the problem under nonstationary heating of the enclosing structure in the form of unrestricted plate with boundary conditions of the III kind. According to the known relations and graphs alterations in surface temperature depending on warm-up time, on thermal resistance of constructions and on arguments of Fo and Bi, i. e. initial and boundary conditions are determined. The graphic dependencies that have been obtained show that the surface temperature depends on the thermal resistance, while the temperature at the opposite surface during heat exposure remains practically unchanged during т = 5 h. Thus, if the outside air temperature is altered, then the rate of change of surface temperature or relative temperature 0 make it possible to determine the thermophysical characteristics by solving the inverse problem of thermal conductivity with the use of the converted ratio to determine R as a function R = f(0, т). If the constructed graphic dependencies R = f(0, т) are used at different heat transfer coefficients, then according to the measured temperatures at different time intervals it is possible to determine thermal resistance in the same time intervals and, according to their average value, determine the required resistance to heat transfer R. The estimated ratio of analytical and graphic dependencies that we have obtained demonstrate the adequacy of the conducted full-scale measurements, if the areas with homogeneous temperature field and temperature history are chosen, and they can be used in determining the heat resistance of the enclosing structure in the form of unrestricted plate with boundary conditions of the III kind.

Keywords: enclosing structures, heat transmission resistance, non-stationary mode, arguments of Fo and Bi, unrestricted plate, boundary conditions

For citation: Khroustalev B. M., Sizov V. D. (2018) Determining Heat Transmission Resistance ofEnclosing Structures. Energetika. Proc. CIS Higher Educ. Inst. and Power Eng. Assoc. 61 (1), 47-59. DOI: 10.21122/1029-7448-2018-61-1-47-59 (in Russian)

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

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

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

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

При анализе, решении задач переноса теплоты (массы) в различных материалах используют методы:

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

- вариационные (Био, Канторовича, Лейбензона, Ритца, Трефтца и др.) [7];

- проекционные: интегральных тепловых балансов;

- приведения краевых задач с нелинейными граничными условиями (ГУ) к эквивалентным и нелинейным интегральным уравнениям.

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

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

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

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

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

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

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

На сегодня методы, определяющие теплофизические характеристики, базирующиеся на закономерностях нестационарного теплового потока, разделяются на следующие группы [8]: регулярного теплового режима; основанные на определении параметров нестационарного температурного поля в первой стадии его развития; квазистационарного теплового режима.

Все перечисленные методы в той или иной мере можно использовать для контроля термического сопротивления.

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

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

При этом флюктуации температур при одностороннем подводе теплоты на противоположной плоскости не наблюдаются.

В то же время, используя данные [9], можно сделать вывод, что наиболее близким к описываемому процессу является решение задачи, в отличие от ранее предлагаемого [10], суть которого состоит в следующем.

Начальная температура неограниченной пластины по толщине одинакова. У одной поверхности температура среды известна (задана), на другой поверхности температура равна начальной.

Начальные граничные условия следующие:

Ц т-0 — t0i ; ^ -

дх

— х(8-tx—о); 8 — const; t\x—h — to

х—0

Тогда безразмерная относительная и средняя температуры могут быть рассчитаны из следующих уравнений: - параметр температуры

t— Bi • (1 - п)

t0 Bi + 1 „—i

где An - параметр амплитуды колебаний температуры,

0 — ТГ" = ^BT-T " £ An sin [Vn (1 - n)exp(-Vn2 Fo)], (1) 9-L Bi + 1~T L J

An —--; (2)

n sin Vn (Vn2 + Bi2 + Bi)

параметр средней температуры t -10 Bi

2(Бх + 1) И=1 где Вп - относительная величина,

-£ Bn exp(-v„2Fo), (3)

Bin = — (1 - cosvn); (4)

V

п - число членов ряда; - аргументы:

ат a3Kh х

Fo = ; Bii = "Т; П^,

где аэк - коэффициент теплоотдачи у внешней поверхности, Вт/(м2К); к - толщина слоев ограждающей конструкции, м; а = Хэк - коэффи-

^экРэк

циент температуропроводности, м2/ч; сэк - массовая изобарная теплоемкость, кДж/(кгК); Хэк - коэффициент теплопроводности тела, Вт/(мК); рэк - средняя плотность конструкции, кг/м3; п - параметр толщины.

На основании указанных выше формул [9] построены графические зависимости, из которых можно определить искомые значения температур по следующим выражениям:

- температура поверхности в осевой плоскости пластины (ограждения) при п = 0 и п = 0,5

ч = ч + 0('п - ч); (5)

- средняя температура

1 = ^ + - (о);

- градиент температуры при п = 1,0

— = -и 1о дх к

г Бх

где и =

1 + Б1

Рассмотрим пример при следующих начальных условиях: - для табл. 1:

к = 0,25 м; р = 1100 кг/м3; с = 840 Дж/(кг-К); а = 10,0 Вт/(м2К); Я = 1,0 м2К/Вт; X = 0,25 Вт/(мК); а = 0,00097 м2/ч; 10 = 0 °С; 1к = -5,0 °С;

- для табл. 2:

Я = 2,0 м2К/Вт; X = 0,125 Вт/(м К); а = 0,00048 м2/ч; Б1 = 20,0; ^ = 0 °С; ^ = -5,0 °С; а = 10,0 Вт/(м2К).

Рассчитанные данные сводим в табл. 1 и 2, по которым можно построить графические зависимости ^ = Дт, Я, п) (рис. 1).

Таблица 1

т, ч Fo Bi 0Т| = 0 tn = 0, °С 0| = 0,5 tn = 0,5, °С Gn = 1,0 dt/dxn = j 0, оС/м '

1 0,015 10 0,65 -3,25 0,01 0,05 0 0

2 0,030 10 0,74 -3,70 0,03 -0,15 0 0

5 0,070 10 0,82 -4,10 0,1 -0,85 0,80 1,6

10 0,155 10 0,86 -4,30 0,29 -1,45 0,38 7,6

Таблица 2

т, ч Fo Bi 0| = 0 tn = 0, °С 0| = 0,5 tn = 0^ °С Gn = 1,0 dt/dxn = 10, оС/м '

1 0,03 20 0,80 -4,00 0,05 -0,25 0 0

2 0,06 20 0,85 -4,25 0,12 -0,60 0,40 0,8

5 0,15 20 0,90 -4,50 0,31 -1,55 0,40 8,0

10 0,30 20 0,92 -4,60 0,42 -2,10 0,80 16,0

Рис. 1. Зависимости t = fir, R, n): —■--при R = 1,0 м2К/Вт;

—•--при R = 2,0 м2К/Вт

Fig. 1. Dependencies t = ft, R, n): —■--when R = 1.0 m2K/W;

—•--when R = 2.0 m2K/W

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Из рис. 1 следует, что изменения температур поверхности tn=0 и в осевой плоскости tn=05 неограниченной пластины зависят от термического сопротивления материала ограждающей конструкции R и времени теплового воздействия т.

Вместе с тем при продолжительности воздействия т до 2 ч градиент температуры dt/dx на противоположной поверхности практически равен 0 при R = 1,0 м2К/Вт и при R = 2,0 м2К/Вт изменяется незначительно, а только при т > 5 ч изменения температур существенны, т. е. в этом случае при т < 3 ч справедливы ГУ III рода.

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

В качестве примера можно привести таблицы полученных результатов:

- начальные условия для табл. 3:

R = 2,0 м2К/Вт; h = 0,25 м; р = 1100 кг/м3; с = 0,840 кДж/(кг-К);

а = 10,0 Вт/(м2К); X = 0,125 Вт/(мК); а = 0,00049 м2/ч;

Bi = 20,0; to = 0 °C; tK0H = -10 °С;

- начальные условия для табл. 4:

R = 1,0 м2 К/Вт; h = 0,25 м; р = 1100 кг/м3; с = 0,840 кДж/(кг-К);

а = 10,0 Вт/(м2К); X = 0,25 Вт/(мК); а = 0,00097 м2/ч;

Bi = 10,0; to = 0 °C; tmB = -10 °С.

Таблица 3

т, ч Fo 0П = 0 ¿п = 0 0П = 0,5 ¿П = 0,5 0П = 1,0 ¿П = 1,0

0,2 0,0016 0,5749 -5,75 0,1629 -1,63 0 0

0,5 0,0039 0,6468 -6,47 0,1668 -1,67 0 0

0,8 0,0062 0,6949 -6,95 0,1667 -1,67 0 0

1,0 0,0780 0,7187 -7,19 0,1654 -1,65 0 0

Таблица 4

т, ч Fo 0Т| = 0 ¿П = 0 0Т| = 0,5 ¿П = 0,5 0П = 1,0 ¿П = 1,0

0,2 0,00312 0,4318 -4,32 -0,0049 0,049 0 0

0,5 0,00779 0,5398 -5,40 -0,0009 0,009 0 0

0,8 0,01247 0,6025 -6,02 0,0003 -0,003 0 0

1,0 0,01558 0,6319 -6,32 0,0015 -0,015 0 0

На основании полученных расчетных данных построены зависимости температур при п = 0; 0,5 и 1,0 в разные промежутки времени и при различных величинах термических сопротивлений, которые представлены на рис. 2.

Графики рис. 2 показывают, что изменения температуры поверхности и в осевой плоскости бесконечной пластины зависят от заданных термических сопротивлений и времени воздействия отрицательных наружных температур, и в тот же интервал воздействия (ттах = 1,0 ч) градиент температуры не достигает противоположной границы пластины (внутренней поверхности наружного ограждения) при заданных начальных условиях.

Рис. 2. Зависимости t = fT, R, n): —■--при R = 1,0 м2К/Вт; —•--при R = 2,0 м2К/Вт;

—А--рассчитанные значения температур tn = 0,5 при R = 1,0 м2-К/Вт малы и графики tn=05

при R = 1,0 м2-К/Вт и tn = 10 при R = 1,0 м2-К/Вт и R = 2,0 м2-К/Вт фактически совпадают

с осью координат

Fig. 2. Dependencies t =fix, R, n): —■--when R = 1.0 m2K/W; —•--when R = 2.0 m2K/W;

—А--calculated values of temperatures tn = 0.5 when R = 1.0 m2-K/W are small

and graphics tn = 0.5 when R = 1.0 m2K/W and tn = 10 when R = 1.0 m2-K/W and R = 2.0 m2K/W actually coincide with the axis of the coordinate

При этом продолжительность нестационарного режима сохраняется до времени

52 0,252

т = — = —--= 64 ч.

a 0,00097

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

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

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

В случае нагревания (охлаждения) конструкции в виде неограниченной пластины при ГУ III рода для определения относительной избыточной температуры можно воспользоваться выражением (1).

Из (1) видно, что относительная температура - функция теплофизи-ческих характеристик ограждения, коэффициента теплоотдачи у наружной поверхности и продолжительности измерений, которые входят в критерии Fo и Bi.

Для определения термического сопротивления теплопередаче Яэк воспользуемся соотношениями:

* эк = ^; а = 1 Дэк = 1*4

Яэк СэкР *эк

В! = аА. = 0^эк; Fo = ^ =-Т-.

* эк 8,. СэкРэАЛк

Подставляя приведенные соотношения в (1), получим уравнение , из которого определим сопротивление теплопередаче ограждения:

е= wi) - Wo) _ aR(1 -n)

^<*>,т i) Wo) aR (

-v2

f -X ^„sin [vo (1 -n) ]>

1 „=1

(6)

xexp

8гСэкРэкR3K У

где а - суммарный коэффициент теплоотдачи а = ак +ал; ак, ал - расчетные конвективный и радиационный коэффициенты теплоотдачи.

Из (6) следует, что, определив относительную температуру 9 в любой момент времени, пока тепловой потенциал не достиг противоположной поверхности пластины, по результатам натурных измерений можно найти и термическое сопротивление конструкции.

^ 8.

Для определения сопротивления теплопередаче Яэк =2—— или эф-

*эк

фективного коэффициента теплопроводности *эк вначале задаются ориентировочные расчетные теплофизические характеристики, ориентировочные величины коэффициентов теплоотдачи у наружной поверхности аэк = = ак +ал в пределах нормируемых. По этим значениям определяются критерии Б1 и Бо при различных *эк, сэк, рэк в заданные кратные промежутки времени между т0 и ттах, за которые температурный потенциал не изменит температуру внутренней поверхности ограждения, т. е. ^х=8 т > ).

По полученным расчетным величинам Б1 и Бо, по соотношению (6) или построенным номограммам [9] вычисляются зависимости 9 = ДБо, Б1), по которым можно опосредованно определить зависимость Я = Д9, т) в графическом виде при различных аэк.

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

По разности температур t(0 т ) и ш т ) находят конкретную величину а

по известным соотношениям.

Величину Яэк в зависимости от полученных а и х;- можно определить по графику Я = Д9, т) в разные промежутки времени, найдя ее среднее зна-

_ Я + я2 + ... + яп

чение Я =- или используя преобразованное относитель-

п

но Я соотношение (6).

Все изложенное можно показать на примере.

По выражению (6) построен график зависимости 9 = ДБо, Б1), по которому можно опосредованно определить зависимость Я = Д9, т), подставив различные значения 9 в определенные промежутки времени. Например, для а = 5 Вт/(м2К) при начальных условиях:

to, оС 9, °С а, Вт/(м2-К) h, м с, кДж/(кгК) р, кг/м3 R, м2К/Вт т, ч

0 -5 5 0,25 0,84 1000 0,6 1,0 1,4 1 2 5

Расчет представлен в табл. 5.

Таблица 5

Результаты расчета 0 при а = 5 Вт/(м2К) The results of the calculation 0 when а = 5 W/(m2-R)

R, м2К/Вт 0,6 1,0 1,4

X, Вт/(мК) 0,42 0,25 0,18

а, м2/ч 0,0016 0,0010 0,0007

т, ч 1 2 5 1 2 5 1 2 5

Fo 0,026 0,052 0,130 0,016 0,031 0,078 0,011 0,022 0,056

Bi 3,0 5,0 7,0

9 0,36 0,46 0,58 0,44 0,53 0,66 0,51 0,57 0,72

По результатам табл. 5 построены зависимости R = f(9, т) при а = 5 Вт/(м2К) (рис. 3).

0,75

9

0,65 0,60 0,55 0,50 0,45 0,40 0,35 0,30

1

2 3 4 т, ч

Рис. 3. Зависимости R = f(9, т) при а = 5 Вт/(м2 К) Fig. 3. The dependencies R = f(9, т) when а = 5 W/( m2K)

5

Рассмотрим случай при a = 10 Вт/(м-К), h = 0,25 м для начальных условий:

to, оС 9, °С a, Вт/(м2-К) h, м с, кДж/(кг К) р, кг/м3 R, м2 К/Вт т, ч

0 -5 10 0,25 0,84 1000 0,6 1,0 1,4 1 2 5

Расчеты представлены в табл. 6.

Таблица 6

Результаты расчета 0 при a = 10 Вт/(м2К) The results of the calculation 0 when a = 10 W/(m2-R)

R, м2К/Вт 0,6 1,0 1,4

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

X, Вт/(мК) 0,42 0,25 0,18

а, м2/ч 0,0016 0,0010 0,0007

т, ч 1 2 5 1 2 5 1 2 5

Fo 0,026 0,052 0,130 0,016 0,031 0,078 0,011 0,022 0,056

Bi 6,0 10,0 14,0

0 0,55 0,66 0,77 0,64 0,72 0,82 0,69 0,76 0,84

По результатам табл. 6 построены зависимости R = ßQ т) при a = 10 Вт/(м2К) (рис. 4).

Fig. 4. The dependencies R = ßQ, т) when a = 10 W/(m2K)

Определение сопротивления теплопередаче по графическим зависимостям, построенным на основании расчетных величин (рис. 3 и 4) с использованием (1) и (6), по которым в [9] построены номограммы для определения опосредованно R = ßQ т), производится следующим образом.

В процессе исследований были определены температуры в разные промежутки времени (табл. 7).

Таблица 7

Время т, ч * ОС (0,т0), С * °С * °С * °С 0 Яэк,,-, м2-К/Вт

1 -5,0 -10,0 -8,10 -5,0 0,62 0,90

2 -5,0 -10,0 -8,56 -5,0 0,73 1,11

5 -5,0 -10,0 -9,05 -5,1 0,81 0,93

По рис. 4 получены значения Яэк1 = 0,90; Яэк2 = 1,11; Яэк5 = 0,93.

0,90 + 1,11 + 0,93

Средняя из измеренных величин Я,

эк.ср

3

■ = 0,98 Вт/(м-К)

определяет искомое сопротивление теплопередаче ограждающей конструкции. При этом время окончания измерений составило 5 ч, когда *(Х=3,ттах) - *(0,т0).

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

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

ВЫВОДЫ

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

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

ЛИТЕРАТУРА

1. Руденко, Н. Н. Влияние массивности на максимальный тепловой поток / Н. Н. Руденко, И. Н. Волошановская // Строительство-2001: материалы Междунар. науч.-практ. конф. Ростов н/Д: Рост. гос. строит. ун-т, 2001. С. 84-87.

2. Об утверждении требований к энергетическому паспорту, составленному по результатам обязательного энергетического обследования, и энергетическому паспорту, составленному на основании проектной документации и правил направления копий энергетического паспорта, составленного по результатам обязательного энергетического обследования: Министерство энергетики Российской Федерации, приказ от 19.04.2010 № 182.

3. Фурсова, И. Н. Алгоритм исследования плотности теплового потока через ограждение при нестационарных тепловых условиях / И. Н. Фурсова, А. А. Капралов // Сроитель-ство-2013: материалы Междунар. науч.-практ. конф. Ростов н/Д: РГСУ, 2013. С. 87-88.

4. Тепло- и массообмен: в 2 ч. / Б. М. Хрусталев [и др.]. Минск: Белорус. нац. техн. ун-т, 2007. Ч. 1. 606 с.

5. Фокин, К. Ф. Строительная теплотехника ограждающих частей зданий / К. Ф. Фокин, под ред. Ю. А. Табунщикова, В. Г. Гагарина. 5-е изд., испр. и пересмотр. М.: АВОК-Пресс, 2006. 256 с.

6. Лыков, А. В. Теория теплопроводности / А. В. Лыков. М.: Высш. шк., 1967. 600 с.

7. Био, М. А. Вариационные принципы в теории теплообмена / М. А. Био. М.: Энергия, 1975. 209 с.

8. Методы определения теплопроводности и температуропроводности / А. Г. Шашков [и др.]; под. ред. А. В. Лыкова. М.: Энергия, 1973. 336 с.

9. Пехович, А. И. Расчеты теплового режима твердых тел / А. И. Пехович, В. М. Жидких. Л.: Энергия, 1968. 304 с.

10. К вопросу диагностики сопротивлений теплопроводности наружных ограждающих конструкций зданий / Б. М. Хрусталев [и др.] // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2010. № 4. С. 36-43.

Поступила 27.06.2017 Подписана в печать 25.08.2017 Опубликована онлайн 29.01.2018

REFERENCES

1. Rudenko N. N., Voloshanovskaya I. N. (2001) The Effects of Massiveness on the Maximum Heat Flux. Stroitel'stvo-2001: Materialy Mezhdunar. Nauch.-Prakt. Konf. [Construction-2001: Materials of the International Scientific-and-Practical Conference]. Rostov-on-Don, Rostov State University of Civil Engineering, 84-87 (in Russian).

2. On Approval of Requirements for Energy Passport, Based on the Results of Mandatory Energy Audit and Energy Passport, Compiled on the Basis of Project Documentation and Rules for Sending Copies of the Energy Passports Compiled in Accordance with the Results of Mandatory Energy Audit. Order of the Ministry of Energy of the Russian Federation, 19.04.2010, No 182. Electronic Fund of Legal Normative and Technical Documentation. Available at: http://docs.cntd.ru/document/902220954 (Accessed: 4 September 2017) (in Russian).

3. Fursova I. N., Kapralov A. A. (2013) The Algorithm to Examine the Density of the Heat Flux through the Fence under Non-Stationary Thermal Conditions. Stroitel'stvo-2013: Materialy Mezhdunar. Nauch.-Prakt. Konf. [Construction-2013: Materials of the International Scienti-fic-and-Practical Conference]. Rostov-on-Don, Rostov State University of Civil Engineering, 87-88 (in Russian).

4. Khroustalev B. M., Nesenchuk A. P., Timoshpol'skii V. I., Akel'ev V. D., Sednin V. A., Kop-ko V. M., Nerez'ko A. V. (2007) Heat and Mass Exchange. Part 1. Minsk: Belarusian National Technical University. 606 (in Russian).

5. Fokin K. F., Tabunshchikov Yu. A., Gagarin V. G. (eds.) (2006) Construction Heat Engineering of Enclosing Parts of Buildings. 5th ed. Moscow, AVOK-Press. 256 (in Russian).

6. Lykov A. V. (1967) Theory of Thermal Conductivity. Moscow, Vysshaya Shkola Publ. 600 (in Russian).

7. Bio M. A. (1975) Variational Principles in the Theory of Heat Transfer. Moscow. Energiya Publ. 209 (in Russian).

8. Shashkov A. G., Volokhov G. M., Abramenko T. N., Kozlov V. P., Lykov A. V. (ed.) (1973) Methods of Determination of Thermal Conductivity and Thermal Diffusivity. Moscow. Energiya Publ. 336 (in Russian).

9. Pekhovich A. I., Zhidkikh V. M. (1968) Calculations of the Thermal Regime of Solids. Leningrad, Energiya Publ. 304 (in Russian).

10. Khroustalev B. M., Akeliev V. D., Sizov V. D., Zolotoreva I. M. (2010) Towards Diagnostics of Heat Conductivity Resistance of Building Exterior/Outside Fencing Constructions. Ener-getika. Izvestiya Vysshikh Uchebnykh Zavedenii i Energeticheskikh Ob 'edinenii SNG = Ener-getika. Proceedings of CIS Higher Education Institutions and Power Engineering Associations, (4), 36-43 (in Russian).

Recеived: 27 June 2017 Accepted: 25 August 2017 Published online: 29 January 2018

i Надоели баннеры? Вы всегда можете отключить рекламу.