CC BY
15
ISSN 2304-120X
ниепт
научно-методический электронный журнал
Елаховский Д. В. Оценка температурных полей элементов строительных конструкций в программе физического образования студентов строительной специальности // Научно-методический электронный журнал «Концепт». - 2017. - № V7. - 0,3 п. л. - URL: http://e-kon-cept.ru/2017/170154. htm.
ART 170154 УДК 378.147
Елаховский Дмитрий Вячеславович,
кандидат физико-математических наук, доцент кафедры общей физики ФГБОУ ВО «Петрозаводский государственный университет», г. Петрозаводск [email protected]
Оценка температурных полей элементов строительных конструкций в программе физического образования студентов строительной специальности
Аннотация. Внедрение в программу курса общей физики элементов строительной физики с целью повышения мотивационной составляющей процесса обучения предполагает более глубокое рассмотрение физических явлений, играющих важную роль в строительстве, но выходящих за рамки традиционного курса физики. Рассмотренный нами ранее механизм теплопроводности, ответственный за теплообмен, соответствовал стационарности его реализации. Реальные условия эксплуатации зданий осуществляются при изменяющейся температуре наружного воздуха и работе отопительных систем, с учетом этого температурное поле ограждения уже не является стационарным и его расчет основан на теории нестационарной теплопроводности. Учет начальных и граничных условий значительно осложняет процесс получения аналитических решений для наиболее важных ситуаций, поэтому используются альтернативные математические методы, рассмотренные на примере одномерной модели.
Ключевые слова: стационарная и нестационарная теплопроводность, температурное поле, уравнение теплопроводности, метод конечных разностей. Раздел: (01) отдельные вопросы сферы образования.
Использование элементов строительной физики в базовом курсе физического образования студентов строительной специальности университетов способствует повышению мотивационной составляющей процесса обучения [1]. Включение раздела строительной физики под названием «Физические основы архитектурной климатологии» вызвано тем обстоятельством, что эксплуатация зданий различного предназначения предполагает наличие соответствующего теплового режима, обеспечивающего ощущение теплового комфорта людей, реализацию производственных процессов в соответствии с техническим заданием и долговечность элементов конструкций самого здания. Тепловая обстановка помещения формируется в зависимости от температуры наружного воздуха и температуры внутри помещения, наличия отопительного оборудования и эксплуатационных характеристик. Требуемые условия санитарных норм помещений вызывают необходимость при проектировании ограждающих конструкций зданий учитывать как наружные условия их эксплуатации, так и полный спектр процессов, протекающих в помещении. В работе [2] достаточно подробно (по сравнению с традиционным курсом общей физики) рассмотрен процесс теплопроводности, существенно влияющий на формирование климатической среды помещения. Учитывая рамки журнальной статьи, данный процесс рассматривали в условиях, обеспечивающих его стационарность. На практике обеспечить такие условия достаточно сложно. Среди факторов, нарушающих режим стационарности, следует отметить, например, суточные колебания температуры наружного воздуха, достигающие в некоторых регионах 300. Система отопления помещений (особенно печи, использующие дрова) также не обеспечивает постоянство температуры воздуха в течение суток,
ниегп
issn 2304-120X Елаховский Д. В. Оценка температурных полей элементов строительных конструкций в программе физического образования студентов строительной специальности // Научно-методический электронный журнал «Концепт». - 2017. - № V7. - 0,3 п. л. - URL: http://e-kon-cept.ru/2017/170154. htm.
научно-методический электронный журнал
и это не способствует тепловому комфорту. Ослабление влияния колебательных процессов воздуха внутри помещений связано со способностью ограждения аккумулировать тепловую энергию за счет их теплоемкостных свойств. Все сказанное вызывает необходимость оценки нестационарных полей температуры, что будет способствовать более детальному знакомству студентов с проблемами формирования тепловой среды помещений и улучшению их профессиональных навыков.
Как отмечалось в [3], явление теплопроводности играет важную роль в проектировании конструктивных элементов зданий, но в традиционном курсе физики это явление представлено достаточно поверхностно в виде уравнения Фурье и зависимости коэффициента теплопроводности от различных параметров. Поэтому необходимо знакомить студентов с наиболее важными особенностями теплопроводности, особенно твердых тел, а также с некоторой терминологией, не встречающейся в обычном учебнике физики. Общее уравнение теплопроводности имеет вид:
CT г
Р■Cv ¥
Cjx j Cj , м,
— + — + — | + q. (1)
dx dy dz
где р - плотность материала ограждения, су - удельная теплоемкость, ^ =-ж(ёт/ёх)
- плотность теплового потока вдоль оси х (аналогично для других проекций вектора ]). Стоящая в правой части уравнения (1) величина я учитывает мощностные характеристики внешних и внутренних источников тепловой энергии, но так как они находятся вне строительного ограждения, то их наличие при решении дифференциальных уравнений учитывается заданием граничных условий на поверхностях ограждения. С учетом этого уравнение теплопроводности принимает вид:
CT d ( CT 1 d Р-cV — = —I ж— 1 + —
dt dx ^ dx) dy
r CT 1 d ( dT
ж— + -|ж— |. (2)
dy) dz ^ dz
Коэффициент теплопроводности ® можно считать постоянной величиной благодаря его слабой зависимости от температуры. Тогда уравнение (2) принимает вид:
CT
р •Cv-= ®
dt
{ d 2Т d 2Т d 2ТЛ ■ + —^ + -
vdx2 dy2 dz2 )
(3)
Отношение к = ® /(сур) называется температуропроводностью данного материала и характеризует интенсивность нестационарных процессов. Тогда в окончательном виде уравнение нестационарной теплопроводности принимает вид:
— = кдт. (4)
где Л - оператор Лапласа.
Данное уравнение - это дифференциальное уравнение в частных производных второго порядка, а постоянство коэффициента теплопроводности обеспечивает его линейность. Это обстоятельство облегчает процедуру его решения за счет использования определенных математических методик, однако явный или удобный для конкретных расчетов аналитический вид решения данного уравнения может быть получен далеко не для всех важных с точки зрения практических приложений случаев. Наиболее используемое математическое обеспечение - это операционный метод и
ниегп
issn 2304-120X Елаховский Д. В. Оценка температурных полей элементов строительных конструкций в программе физического образования студентов строительной специальности // Научно-методический электронный журнал «Концепт». - 2017. - № V7. - 0,3 п. л. - URL: http://e-kon-cept.ru/2017/170154. htm.
научно-методическии электронный журнал
метод частных производных. Аналитическое решение данного уравнения, как правило, представляется в виде быстро сходящихся бесконечных рядов в определенных интервалах изменения переменных или в виде протабулированных специальных функций. Главная проблема решения такого уравнения связана с необходимостью использования начальных и граничных условий, которым должно удовлетворять решение данного уравнения для конкретной задачи:
1. Искомое нестационарное распределение температур Т (х, у, z, ^ должно при t = 0 соответствовать начальному распределению:
Т (х, у, ^ 0) = f (х, у, ^ (5)
2. На поверхностях ограждения для любой конкретной задачи должны быть заданы поддерживаемые на этих поверхностях температуры или подводимые потоки. К примеру, если поверхность граничит с воздушной средой заданной температуры Твозд, то должно выполняться условие:
ст
-х—
Сх
= «( Тпов - Твозд ) , (6)
3. Для неоднородных ограждений - непостоянная величина) решения для соседних областей Т1 (х, у, z, ^ ®1) и Т2 (х, у, z, ^ ®2) на границе зон должны удовлетворять условию непрерывности температуры
Т1 гран — Т2 гран (7)
и непрерывности тепловых потоков
сТ,
-®i —
сХ
сТ7
= -х2—i сХ
(8)
Перечисленные начальные и граничные условия весьма затрудняют решение уравнения нестационарной теплопроводности (особенно в случаях сложной конфигурации ограждений). В связи с этим аналитические решения удается получать для плоских, цилиндрических и сферических поверхностей. Даже для стационарных температурных полей при решении уравнения Лапласа с учетом граничных условий (5)-(8) аналитический вид удается получить только для простых конфигураций.
Это обстоятельство вызывает необходимость замены реального решения уравнения его приближенным численным интегрированием методом конечных разностей.
Для простоты (не умаляющей суть процедуры) рассмотрим одномерную задачу.
8Т - к 81!. (9)
д 8х
Разобьем ось х на ряд интервалов шириной d (см. рисунок) и пронумеруем значения температур на границах интервалов: Т1, Т2, ... Tм, На данном 1-м интервале
Гат ^ т, - т..
производную 8Т дх заменим на I — I - ——1-1 и т. д.
I Ах I d
ISSN 2304-120X
ниепт
научно-метолический электронный журнал
Елаховский Д. В. Оценка температурных полей элементов строительных конструкций в программе физического образования студентов строительной специальности // Научно-методический электронный журнал «Концепт». - 2017. - № V7. - 0,3 п. л. - URL: http://e-kon-cept.ru/2017/170154.htm.
Метод конечных разностей [4]
Аналогично поступаем со вторыми производными:
А (CT
dx l dx
А
CT dx
Ax
AT Ax
AT Ax
= Ti+1 - 2Tj + TM d2
Тогда одно уравнение в частных производных (9) заменится системой обыкновенных связанных дифференциальных уравнений:
dTi_kT+i - 2Ti + Ti-i dt d2
(10)
где i = 1, 2, ... с заданными начальными условиями Т|(0) и определенными граничными условиями для крайних значений Т|.
Если еще и время разбить на малые конечные интервалы h и повторить вышеприведенные действия, то получим алгоритм для вычисления значений температур спустя интервал времени К
T - 2T + T T = T +h ■ k^1 2Ti + Ti-1
d2
(11)
Найдя отсюда все Т|, то есть распределение температур в момент времени t + ^ можно аналогичным путем вычислить ТГ в момент времени t + 2h и т. д. Для большей точности необходимо уменьшать интервалы d и 1л, что резко увеличивает объем необходимых вычислений. Современные ЭВМ легко справляются с этой задачей, поэтому метод конечных разностей с успехом используется при анализе нестационарных температурных полей.
Рассмотрим пример аналитического решения уравнения (9). При этом используется такая модель.
Плоская стенка толщиной d и температуропроводностью k имеет при t < 0 температуру, равную температуре окружающей среды Тн. В момент времени t = 0 за счет внутреннего источника тепла температура левой поверхности мгновенно приняла значение Тв и в дальнейшем не изменялась. Требуется вычислить нестационарное температурное поле в стенке Т (х, ^ в любой ее точке (0 < х < d) в любой последующий момент времени (0 < t < °°). Для этого требуется решить одномерное уравнение теплопроводности (9) с начальными и граничными условиями:
Т (х, 0) = Тн при х > 0; (12)
T (0, 0 = Tв, Ж 0 = Тн а > 0). (13)
ниегп
issn 2304-120X Елаховский Д. В. Оценка температурных полей элементов строительных конструкций в программе физического образования студентов строительной специальности // Научно-методический электронный журнал «Концепт». - 2017. - № V7. - 0,3 п. л. - URL: http://e-kon-cept.ru/2017/170154. htm.
научно-методическии электронный журнал
Для простоты будем считать, что коэффициенты теплоотдачи на граничных поверхностях ав и ан бесконечно большие.
В подобной формулировке при решении такой задачи значения параметров Тн, Тв, d, к варьируются в широких пределах. Но некоторые преобразования позволяют унифицировать данную задачу.
Исходное уравнение (9) запишем в эквивалентном виде (для разности температур Т(х, ^ - Тн) с учетом постоянства Тн
8(Т - Тн) к 82(Т - Тн)
8г дх2
с начальными и граничными условиями:
Т (х, 0) - Тн = 0, Т 0 - Тн = 0, Т (0, 0 - Тп — Тв - Тн.
Таким образом, при одинаковой величине перепада температур Тв - Тн температурный ход для разности температур Т - Тн не зависим от абсолютного значения температуры среды Тн.
Если ввести в рассмотрение относительный разогрев
6(х,1) - Т(ТМ)т Т, (14)
в н
(при этом 0 < 0 < 1), то искомое уравнение принимает вид:
= (15)
5t 5x2
с начальными и граничными условиями:
0 (х, 0) — 0, 0 0 = 0, 0 (0, 0 = 1.
Это уравнение не содержит температурных параметров Тн и Тв, тем самым для данной стенки при самых различных значениях Тн и Тв относительный разогрев 0 (х, ^ в тех же точках и в те же моменты времени t будет одинаков.
Дальнейшее упрощение задачи связано с введением относительных координат, т. е. для стенок различной толщины d следует анализировать относительные разо-гревы не на одинаковых расстояниях х от края пластины, а на одинаковых относительных расстояниях а — х / d, где 0 < а < 1. Переходя от абсолютных координат к относительным, имеем:
8 1 8 82 1 82 ах - а • аа, — ---, —- -——-.
8х d 8а 8х а 8а
Тогда искомое уравнение переходит в
8?-^ (16)
81 а2 8а2
Из этого уравнения следует, что ход относительного разогрева в соответственных точках для стенок разной толщины d одинаков только в том случае, если температуропроводность стенок к пропорциональна d2.
ниегп
issn 2304-120X Елаховский Д. В. Оценка температурных полей элементов строительных конструкций в программе физического образования студентов строительной специальности // Научно-методический электронный журнал «Концепт». - 2017. - № V7. - 0,3 п. л. - URL: http://e-kon-cept.ru/2017/170154. htm.
научно-методический электронный журнал
Если разделить обе части дифференциального уравнения (16) на k АЗ2, внести эту
к 1
величину под знак дифференциала и ввести безразмерную величину у = — г = —, то по-
<г г
лучим уравнение:
(17)
5у 5а
при этом 0 (а, 0) = 0, 0 (1, у) = 0, 0 (0, у) = 1.
Величина Г - естественный масштаб времени для всей данной группы задач:
[ t* ] =
м2
= с
м с
Таким образом, для всей совокупности значений параметров Тн, Тв, k и d получим единое уравнение и единое его решение - некоторую функцию 0 (а, у). Это значит, что для всех стенок на одинаковых относительных расстояниях а = х ^ в сходственные моменты времени у = t /Г будут одинаковые значения относительных разогревов 0 = (Т - Тн) / (Тв - Тп).
В том случае, когда у = t /Г ^ те, получается стационарное состояние, при этом
0 (а, у) уже не зависит от времени (^5 = 0) и уравнение (17) принимает вид:
^ = 0. (18)
аа
Решением этого уравнения, удовлетворяющего граничным условиям при а = 0 и а = 1, будет:
0 (а, «) = 0стац(а) = 1 - а. (19)
Для конечных значений у при аналитическом решении уравнения (17) это решение представляется в виде бесконечного ряда по различным функциям. В частности, при у > 1 (то есть t > Г) используются тригонометрические и показательные функции, которые быстро сходятся для больших времен. Кривые 0 (а, у) протабулированы для ряда значений у = t /Г.
Из решения уравнения (17), представленного в виде ряда по тригонометрическим функциям, можно получить оценку значения уk « 0,1, при котором нестационарное распределение 0 (а, у) практически переходит в стационарное. Отсюда можно оценить время тепловой релаксации:
а2
tk~ykt*~0.1 —. (20)
к
Если стенка имеет толщину d = 0,1 мм и температуропроводность k = 10-6 м2/с, то время тепловой релаксации tk « 20 мин.
Если наружная температура Тн изменяется так, что 6Т << Тв - Тн (1+2° в течение часа), то распределение температуры внутри стенки все время практически «подстраивается» под изменение Тн, что обеспечивает линейное снижение температуры, как и в стационарном случае, но наклон (Тв - Тн) ^ постепенно будет изменяться в соответствии с изменением Тн В этом случае формула (19) описывает квазистационарное распределение температуры и, перейдя к размерным переменным, получим:
ISSN 2304-120X
ниепт
научно-методический электронный журнал
Елаховский Д. В. Оценка температурных полей элементов строительных конструкций в программе физического образования студентов строительной специальности // Научно-методический электронный журнал «Концепт». - 2017. - № V7. - 0,3 п. л. - URL: http://e-kon-cept.ru/2017/170154. htm.
T_(x,t) = TH(t) + [Тв -Tn(t)]
d - x d ■
(21)
Таким образом, рассмотрение нестационарной теплопроводности расширяет кругозор будущих строителей, что позволяет оценить возможности теоретического рассмотрения теплофизических свойств ограждающих конструкций, обеспечивающих режим комфортного проживания людей.
Ссылки на источники
1. Елаховский Д. В. Повышение мотивационной составляющей физического образования студентов строительной специальности вуза // Физическое образование в вузах. - 2013. - Т. 19. - № 4. - С. 74. Елаховский Д. В. Элементы архитектурной климатологии как фрагменты физического образования студентов строительной специальности университетов // Научно-методический электронный журнал «Концепт». - 2015. - № 5 (май). - С. 86-90. - и^: http://e-koncept.ru/2015/15146.htm. Там же.
Тодес О. М. Курс лекций и заданий по строительной физике. - Л., 1972.
2.
Dmitry Elakhovskiy,
Candidate of Physical and Mathematical Sciences, Associate Professor, General Physics Chair, Petrozavodsk
State University, Petrozavodsk
Building components temperature field estimation in the program of physics for civil engineering students
Abstract. Introduction of building physics elements in the course of general physics with the purpose of learning process motivational component increasing involves a deeper studying of the physical phenomena that play an important role in construction, but beyond the scope of traditional physics course. The mechanism of thermal conductivity considered earlier is responsible for heat transfer in steady state conditions. But the temperature field of retaining walls in real environment is not stationary because of varying outdoor air temperature and heating system operation. That is why its estimation should be based on the theory of unsteady-state conductivity. Calculation of initial and boundary conditions significantly complicates the process of analytical solutions obtaining for the most important cases so we suggest the use of alternative mathematical methods presented here on example of one-dimensional model.
Key words: steady-state and unsteady-state thermal conductivity; temperature field; heat conductivity equation; finite-difference method. References
1. Elahovskij, D. V. (2013). "Povyshenie motivacionnoj sostavljajushhej fizicheskogo obrazovanija studentov stroitel'noj special'nosti vuza", Fizicheskoe obrazovanie v vuzah, t. 19, № 4, p. 74.
2. Elahovskij, D. V. (2015). "Jelementy arhitekturnoj klimatologii kak fragmenty fizicheskogo obrazovanija studentov stroitel'noj special'nosti universitetov", Nauchno-metodicheskij jelektronnyj zhurnal "Koncept", № 5 (maj), pp. 86-90. Available at: http://e-koncept.ru/2015/15146.htm.
3. Ibid.
4. Todes, O. M. (1972). Kurs lekcij i zadanij po stroitel'noj fizike, Leningrad.
Рекомендовано к публикации:
Горевым П. М., кандидатом педагогических наук, главным редактором журнала «Концепт»
Поступила в редакцию Received 15.05.17 Получена положительная рецензия Received a positive review 22.05.17
Принята к публикации Accepted for publication 22.05.17 Опубликована Published 31.07.17
www.e-koncept.ru
© Концепт, научно-методический электронный журнал, 2017 © Елаховский Д. В., 2017
977230412017307
