Общие принципы
При решении задачи теплопроводности неизвестной является температура в узлах расчетной схемы (аналогично перемещениям в прочностной задаче). Тогда Уравнение равновесия выглядит:
Задача теплопроводности
Прочностная задача
Матрица теплопроводности формируется из коэффициентов теплопроводности К, что задаются в жесткости конечного элемента. В случае пластин, матрица теплопроводности умножается на толщину элемента Н.
При решении нестационарной задачи теплопроводности в уравнение равновесия включается матрица теплопоглощения (демпфирования) С:
Матрица демпфирования зависит от произведения коэффициента демпфирования С на плотность материала RO (и умноженных на толщину элемента Н, если это пластина).
Типы конечных элементов
Конечные элементы теплопроводности
№ 1505
Стержневой КЭ теплопроводности
№ 1508
Треугольный КЭ теплопроводности
№ 1509
Четырехугольный КЭ теплопроводности
№ 1512
КЭ теплопроводности в форме тетраэдра
№ 1514
Объемный шестиузловой КЭ теплопроводности
№ 1516
Объемный восьмиузловой КЭ теплопроводностиКонечные элементы конвективного теплообмена
№ 1551
Одноузловой КЭ конвективного теплообмена
№ 1555
Двухузловой КЭ конвективного теплообмена
№ 1558
Треугольный КЭ конвективного теплообмена
№ 1559
Четырехугольный КЭ конвективного теплообменаВиды нагрузок
Стационарные виды нагрузок
1. Постоянная температура в узле(или граничное условие первого рода) Когда задана температура одного из узлов:
Аналогом в прочностной задаче есть заданное перемещение узла.
|
Рис. 3. Нагрузка на узел
|
, °С
2. Тепловой поток (или граничное условие второго рода)
Когда внешняя нагрузка 
равна:
По определению тепловой поток – это количество теплоты, которое проходит через изотермическую поверхность (условную поверхность с одинаковой температурой) за единицу времени.
Измеряется в Дж/c*м2 (или Вт/ м2, поскольку 1 
= 1 Дж/c).
В ЛИРА-САПР есть возможность задавать тепловой поток на узел, пластину или объемное тело.
2.1. Тепловой поток на узел
При задании теплового потока на узел необходимо вычислить площадь, через которую проходит поток (в случае рассмотрения балки, это будет площадь сечения). К примеру, если величина теплового потока q = 100 Дж/c*м2, а площадь через которую проходит поток равна 0,1 м2, то величина задаваемой нагрузки на узел Р равна 100*0,1 = 10 Дж/c .
2.2. Тепловой поток на ребро пластины
При задании теплового потока на ребро необходимо чистую величину теплового потока умножать на толщину пластины. Например, если величина теплового потока q = 100 Дж/c*м2, а толщина пластины равна 0,1 м, то величина задаваемой нагрузки на узел Р равна 100*0,1 = 10 Дж/c*м.
2.3. Тепловой поток на объемное тело
В случае задания теплового потока на грань объемного тела: не нужно умножать, ни на площадь, ни на длину ☺. Необходимо только указать номер грани, на которую воздействует поток и величину самого потока q = 100 Дж/c*м2.
3. Конвективный теплообмен (или граничное условие третьего рода)
Когда внешняя нагрузка P равна:
Конвективный теплообмен – это процесс потери тепла телом за счет того, что при движении воздуха происходит соприкосновение отдельных частиц, у которых различная температура. При конвективном теплообмене происходит теплоотдача между поверхностью тела и движущейся средой (например, воздухом).
Создание конвективной нагрузки на схему совершается в 2 этапа:
На первом этапе нужно создать конечный элемент конвекции в месте контакт схемы с воздухом.
- Для стержня – контакт моделируется через одноузловой КЭ конвекции (КЭ №1551).
- Для пластины – контакт моделируется через двухузловой КЭ конвекции (КЭ №1555).
- Для объемного тела – через пластинчатые КЭ конвекции (КЭ №1558 и №1559).
В жесткости такого конечного элемента нужно указать коэффициент конвективной теплоотдачи и площадь контакта (для одноузлового КЭ) или высоту контактной поверхности (для 2х-узлового КЭ). Коэффициент конвективной теплоотдачи h измеряется в Дж/c*м2*оС.
На втором этапе нужно приложить нагрузку на контактный элемент. Через задание нагрузки формируется внешняя температура воздуха.
Пример задания конвективной нагрузки смотрите здесь: «Решение стационарной задачи теплопроводности» и «Решение динамической теплопроводности для задачи огнестойкости».
Нестационарные виды нагрузок
Из всех видов нестационарных нагрузок: непосредственно на узел (элемент) можно задать только изменяемый во времени тепловой поток на узел. Все остальные виды нестационарных нагрузок можно задать, только используя диалоговое окно «Формирование динамических загружений из статических».
1. Изменяемый во времени Тепловой поток
1.1. На узел
Есть три вида изменяемой во времени нагрузки на узел (аналогично прочностной задаче):
- ломанная с произвольным шагом;
- синусоидальная;
- ломанная с равномерным шагом.
Величина задаваемой нагрузки определяется аналогично стационарному тепловому потоку, то есть значение чистого теплового потока умножается на площадь, через которую проходит этот поток.
1.2. На ребро пластины или на грань объемного КЭ
Чтобы задать нестационарный тепловой поток на ребро пластины или на грань объемного КЭ нужно задать стационарный тепловой поток в любом Загружении, кроме Загружения 4. Поскольку, Загружение 4 предназначено для задания демпфирующих нагрузок.
После этого, в окне «Динамика во времени из статических загружений» нужно сформировать Нестационарное загружение из соответствующего стационарного, путем задания закона изменения нагрузки во времени.
Для формирования теплового потока можно использовать следующие законы преобразования:
- ломанный с произвольным шагом (1);
- синусоидальный (2);
- ломанный с равномерным шагом (4).
2. Изменяемый во времени Конвективный теплообмен.
Использовать такой тип нагрузки можно, если температура движущейся среды (к примеру, воздуха) может значительно изменяться во времени.
Как было сказано выше, задание нестационарного конвективного теплообмена происходит так же, как и задание теплового потока на ребро или на грань (то есть, через формирование динамических нагрузок из статических).
3. Изменяемый во времени Лучистый теплообмен
(или радиационный теплообмен)
Когда внешняя нагрузка равна:
– угловой коэффициент
– коэффициент поглощения (степень черноты поверхности конструкции)
– коэффициент излучения (степень излучения источника)
– постоянная Стефана Больцмана (равна 5,67*10-8Вт/м2К4)
Лучистый теплообмен происходит за счет того, что часть тепловой (внутренней) энергии тела (которая существует благодаря механическому колебанию элементарных частиц, из которых состоит тело) превращается в энергию излучения. Энергия излучения - это энергия электромагнитных колебаний с волнами различной длины. Возникают электромагнитные волны за счет колебания заряженных частиц (электронов или ионов), которые входят в состав тела. При попадании лучистой энергии на какое-либо тело, часть этой энергии поглощается, часть – проходит сквозь тело, а часть - отражается. Степень поглощения телом лучевой энергии учитывается через коэффициент . К примеру, для бетона этот коэффициент равен 0,75, а для абсолютно чёрного Тела – 1.
Угловой коэффициент Ф обычно принимается равным 1, он учитывает расположение источника излучения по отношению к поглотителю.
Пример задания лучистой нагрузки смотрите здесь «Решение динамической теплопроводности для задачи огнестойкости».
Можно подвести итоги. Тепло - это энергия, которая возникает за счет движения частиц, из которых состоит тело (атомов или молекул). Всего существует три основных вида передачи тепла: теплопроводность, конвекция и лучистый теплообмен.
Пример. Решение стационарной задачи теплопроводности
Рассмотрим пример расчёта теплопроводности внешней стены кирпичного дома.
Создаем задачу в 15м признаке схемы, который существует специально для решения задач теплопроводности. Рассмотрим участок стены, длиной 1 метр.
ШАГ 1. Геометрия
ШАГ 2. Создание элементов КОНВЕКЦИИ
Чтобы задать температуру воздуха, добавляем стержневые элементы на внутренней и внешней поверхности стены, и меняем их тип на КЭ №1555.
Рис. 11. КЭ модель с элементами конвекции
ШАГ 3. Характеристики материалов
Зададим соответствующие коэффициенты теплопроводности K для слоёв стены. Значения коэффициента теплопоглощения C и удельного веса R0 в статическом расчете не учитываются, поэтому можно их задать равными единицам.
Рис. 12. Свойства теплопроводности соответствующих слоев стены
Для элементов конвекции тоже создаём жесткость, и задаем там коэффициенты конвекции внутреннего и внешнего слоя:
|
|
|
|
Рис.13. Жесткости слоёв стены
ШАГ 4. Внешняя нагрузка
Через внешнюю нагрузку мы задаем температуру воздуха для элементов конвекции. Для этого, в разделе Нагрузки открываем конвективный теплообмен 
и задаём температуру внутри и снаружи стены.
ШАГ 5. Просмотр результатов
Пример. Решение динамической теплопроводности для задачи огнестойкости
Решение задачи теплопроводности - это важный этап при определении огнестойкости конструкции
В Eurocode EN 1992-1-2 (2004) (пункт 4.3.1 (3)) написано «Advanced calculation methods should include calculation models for the determination of: the development and distribution of the temperature within structural members (thermal response model)…» ;
В СТО 36554501-006-2006 «Правила по обеспечению огнестойкости и огнесохранности железобетонных конструкций» (пункт 4.13) сказано: «Расчет огнестойкости и огнесохранности рекомендуется производить по приведенному сечению, когда сечение элемента разбивается на малые характерные участки, нагретые до различных температур, и каждый малый участок приводится к ненагретому бетону с учетом соответствующих понижающих характеристик прочности бетона.», и в пункте 5.3. добавлено: «Температуру бетона определяют теплотехническим расчетом …».
В ДСТУ-Н EN 1991-1-1:2010 (пункт 4.3.1) сказано: «Уточнені методи розрахунку включають розрахункові моделі для визначення наступного: зростання та розподілення температури по всіх елементах конструкції (теплотехнічний розрахунок);…».
Поэтому можно смело сказать, что без теплопроводности никуда ☺
Рассмотрим теперь пример того, как определить распространение тепла от пожара по сечению балки, размерами 160х320 мм. Поскольку наша балка не переменного сечения, то можно рассматривать плоскую задачу, где ширина всех элементов будет равна одному погонному метру балки (100 см).
ШАГ 1. Геометрия
Создаем схему сечения балки в 15м признаке схемы (рис. 17).
|
Назначим жесткость для элементов теплопроводности:
Рис. 16. Свойства теплопроводности
|
Рис. 17. КЭ модель
|
В СТО 36554501-006-2006 «Правила по обеспечению огнестойкости и огнесохранности железобетонных конструкций» (пункт 6.3) для теплотехнического расчета рекомендуют принимать теплопроводность и теплоёмкость бетона по формулам (рассмотрим пример тяжелого бетона на карбонатном заполнителе):
- теплопроводность: λ=1,14-0,00055t, [Дж/c*м*оС.]
- коэффициент теплоёмкости: C=0,71+0,00083t, [кДж/кг*оС.]
где 
– это температура, 
.
Поскольку теплопроводность и теплоёмкость зависят от температуры, то для нашей задачи мы используем свойства бетона, что рассчитаны для температуры 300 °С. Именно с такими характеристиками бетона получаются максимально схожие результаты температурных кривых прогрева сечения с теми, что приведены в нормах.
ШАГ 2. Внешняя нагрузка. ПРЕДЫСТОРИЯ
После этого начинаем задавать граничные условия (внешнюю нагрузку). В Загружении 1 выделяем все узлы схемы, и задаем им постоянную температуру, которая равна 20 °С. Это стандартное значение для предыстории, поскольку считается, что средняя температура конструкции до начала пожара является именно таковой. После этого все узлы окрашиваются в зеленый цвет, рис.19.
ШАГ 3. Внешняя нагрузка. КОНВЕКЦИЯ
Как было уже сказано, существует 3 основных вида передачи тепла. Саму теплопроводность внутри тела мы уже учли путем задания жесткости. Остается конвекция и лучистый теплообмен.
Поэтому, для формирования конвективного теплообмена, мы создадим контактные двухузловые конечные элементы с тех сторон балки, которые предположительно будут подвержены воздействию огня: слева, снизу и справа (Рис. 19, б). Создаем «жесткость» для элементов конвекции, и задаем там толщину нашего сечения Н = 100см, а также коэффициент конвекции а=25 Дж/(м2*с*°С). Коэффициент конвекции, равный 25 – это нормативная величина для стандартного графика развития пожара.
|
После этого, необходимо сформировать внешнюю температуру среды. Для этого используем вид нагрузки под названием заданная температура t:
Нагрузку на конвективные элементы формируем в Загружении 5, при этом Загружения 2-4 оставляем свободными. В окне задания температуры (рис.18) записываем единицу. Чуть позже, на стадии формирования Динамической нагрузки из статической, мы будем задавать график изменения температуры при пожаре. Величины температур из этого графика будут автоматически умножаться на число, которое мы введем в этом окне. Поскольку, числа из графика нам изменять не нужно, то в окне заданной температуры мы просто записываем единицу. |
Рис. 18. Нагрузка на элементы конвекции
|
После задания конвективной нагрузки, контактные элементы окрашиваются в оранжевый цвет.
а)
б)
в)
г)
ШАГ 4. Внешняя нагрузка. ЛУЧИСТЫЙ ТЕПЛООБМЕН (или РАДИАЦИЯ)
Как было сказано раньше, тепловое излучение – это вид теплообмена, который совершается путем поглощения электромагнитных волн. Если в случае конвективной нагрузки, площадь контакта с поверхностью мы задаем через жесткость конвективного элемента (мы задавали Н=100 см), то в случае лучистого теплообмена, площадь контактной поверхности нужно задавать через нагрузку вида Тепловой поток на узел или на ребро пластины, (в случае объемного конечного элемента, тепловой поток нужно задавать на грань элемента). При расчете, на каждом шаге интегрирования, это число будет умножаться на произведение: Ф*εm*εf*σ*(Твнеш4-
Твнутр4), где 
- температура из графика пожара, который будет задаваться на следующем шаге, а 
- температура в узле расчетной схемы.
Сразу появляется вопрос, зачем площадь задавать через нагрузку? И он справедлив ☺
Постараемся найти на него ответ. Для этого, рассмотрим отдельно два варианта формирования лучистого теплообмена.
Вариант 1 (если задавать поток тепла на узел)
Тепловой поток на узел определяется как чистый тепловой поток Дж/c*м2, умноженный на площадь, через которую проходит этот поток [м2]. В случае радиации мы принимаем, что чистый тепловой поток =1. Нам нужно только найти площадь поглощения, которая приходит на каждый узел схемы.
Если мы рассмотрим часть расчетной схемы (рис. 20), где размеры каждого КЭ равны 20х20 мм, а толщина = 1 м, то контактная поверхность, которая припадает на 1 узел равна:
(0,01м +0,01м)*1м =0,02 м2.
Значит, в Загружении 6, на каждый контактный узел задаем тепловой поток, который равен 0,02 м2∙1[Дж/c*м2]= 0,02 [Дж/c], и получаем схему, которая изображена на Рис.19, г.
Вариант 2 (если задавать поток тепла на ребро пластины)
Как и в предыдущем варианте, тепловой поток на ребро определяется как величина чистого теплового потока Дж/c*м2, умноженная, в данном случае, не на площадь, а на толщину пластины Н, [м]. На длину этого же ребра нагрузка умножается автоматически при расчёте (аналогично нагрузке на ребро пластины в прочностной задаче).
|
Чтобы задать нагрузку только на определенные грани конечных элементов, нужно выделить все узлы этих граней, а также выделить сами конечные элементы, которым принадлежат эти грани. При задании Нагрузки на пластины (рис. 21) поставьте галочку для автоматического выбора граней, на которые задается нагрузка.
Таким образом, для нашей схемы, в Загружении 6, на все контактные грани нужно задать тепловой поток, который равен:
1
|
Рис. 21. Выбор загруженных граней
|
∙ 1 м = 1 
.
ШАГ 5. ЗАДАНИЕ ПОЖАРА. Формирование динамического загружения
Для задания огневого воздействия используем стандартный график развития пожара, который определяется по формуле:
– время нагрева в секундах,
– начальная температура (предыстория) = 20 °С.
Или можно воспользоваться таблицей 6.1. из СТО
Нажимаем на кнопку Формирование динамических загружений из статических, и задаем изменение конвекции во времени, которая зависит от температуры огня. Формируем динамику из 5го Загружения (конвекции), (Рис.20, а). Задаем там 21 точку из таблицы 6.1, плюс одна точка на момент времени 0сек, значение температуры в котором, делаем таким же как в предыстории – 20 °С. Время задаем в секундах (левая колонка), температуру – в градусах (правая колонка).
а)
б) 
После этого, аналогично формируем радиацию, которая также зависит от температуры огня (Рис.22, б). Только теперь, формируем динамику из Загружения 6 (радиации), и используем Закон преобразования №15 – Тепловое излучение. И задаем коэффициент поглощения = 0.75 (для бетона), излучения = 1 (для огня), и угловой = 1, поскольку считаем, что огонь расположен непосредственно вблизи балки.
И нажимаем на параметр Выполнить расчет динамики во времени:
ШАГ 6. ПРОСМОТР РЕЗУЛЬТАТОВ
Выше приведен пример температурных полей по сечению 7200 секунде (что равно 120 минутам горения).
В завершение можно добавить, что нет необходимости считать простые сечения вручную, каждое из них пересчитывать в 15 признаке и импортировать в 5й для прочностного расчета. Это можно сделать через вкладку Огнестойкость в Расчетном модуле Железобетон, при задании материалов. После формирования учета огнестойкости, в элементе будет подобрана арматура с учетом предела огнестойкости конструкции.
Но благодаря 15му признаку схему, у пользователя есть возможность рассчитать на огнестойкость сечения или пространственные задачи любой конфигурации и сложности.
Статья "Решение задачи теплопроводности" в формате PDF
Пример "Динамическая теплопроводность"
Заметили ошибку? Выделите ее и нажмите Ctrl+Enter, чтобы сообщить нам.
Комментарии 1