Если одно из понятий незнакомо, вернитесь к предыдущему занятию и воспроизведите его численный пример. После этого переходите к новой постановке.
Место занятия в инженерной методике
Расчётная лабораторная работа · Определение теплового сопротивления и установившегося нагрева
Инженерный вопрос. Остаётся ли расчёт внутри заданной зависимости теплопроводности от температуры?
Температура, тепловое сопротивление и постоянная времени связывают потери с охлаждением и длительностью включения. Это позволяет сравнить материалы, отвод тепла и рабочие циклы по количественным критериям.
Последовательность решения, отчёт и критерии результата
Результат получают из наблюдаемых величин и проверяют в условиях данного практикума.
Последовательность решения
Обосновать постановку и способ измерения. Согласуйте источники мощности, теплопроводность и условия отвода; выберите температуру, поток или тепловое сопротивление как искомую характеристику.
Результат этого занятия. Проверка постоянного предела и объяснение изменения сопротивления материала.
Данные именно этого опыта. Сохраните температурный диапазон, поток и сопротивление при постоянной и табличной теплопроводности.
Оформление наблюдений. Температуры в заданных местах, потоки через границы, мощности по областям и свойства; для сопротивления — определённая разность температур и мощность.
Как обосновать результат
Проверка из практикума.h задаётся извне. Решатель теплопроводности не определяет скорость обтекания и не рассчитывает коэффициент теплоотдачи. Прямое численное сравнение возможно только при совпадении всех исходных данных и постановки.
Независимый контроль. Сопоставьте суммарную мощность источников с выходящим тепловым потоком и аналитическим пределом, допустимым для этой геометрии.
Критерий выполнения. Постоянная таблица воспроизводит Q=4{,}5454545\,\mathrm{Вт}; все полученные температуры находятся в заявленном диапазоне применимости материала либо отмечены как выход за него.
Точность и расхождения. Сеточную устойчивость температуры и потока оценивайте отдельно. Неопределённость теплоотдачи и теплопроводности нельзя устранить уменьшением конечного элемента.
Основание сравнения. Опора — закон Фурье и аналитический тепловой баланс из занятия. Теплоотдача и свойства материала являются исходными данными; одномерный предел применим только при согласованном пути тепла.
Вывод для устройства
Определите доминирующий тепловой путь и влияние материала или границы; решение о допустимом нагреве требует заданной температуры среды и предела материала.
Границы вывода. Результат зависит от источников, среды, теплоотдачи и теплоёмкости. Упрощённый эталон не заменяет неизвестные условия охлаждения реального устройства.
Теория и вывод формул
Коэффициенты k_x и k_y описывают ортотропную теплопроводность в координатных направлениях. Зависимость k(T) в рассматриваемой реализации относится к изотропному материалу; возможность совмещения свойств проверяется по поддерживаемой модели.
Соглашения об обозначениях и единицах
\(i(t)\), \(u(t)\) — мгновенные значения; \(\widehat I\), \(\widehat U\) — комплексные амплитуды; \(\underline I_{\text{д}}\), \(\underline U_{\text{д}}\) — комплексные действующие значения. Для синусоиды \(\widehat I=\sqrt{2}\,\underline I_{\text{д}}\); действующее значение тока равно \(|\underline I_{\text{д}}|\). Нормировка обязательно учитывается в формулах мощности.
\(\Phi=\int_S\mathbf B\cdot\mathrm d\mathbf S\) — поток через ориентированную поверхность; \(\Psi=\sum_k\oint_{C_k}\mathbf A\cdot\mathrm d\mathbf l\) — потокосцепление последовательных витков с согласованными замкнутыми контурами. Только при одинаковом сцепленном потоке каждого витка \(\Psi=N\Phi\). Обе величины измеряются в веберах: число витков безразмерно. Запись «Вб·виток» подчёркивает суммирование.
\(\theta_{\text{м}}\), \(\Omega\) (также \(\omega_{\text{мех}}\)) — механические угол и угловая скорость. Для основной гармоники \(\theta_{\text{э}}=p\theta_{\text{м}}\), \(\omega_{\text{э}}=p\Omega\), где \(p\) — число пар полюсов. Геометрический период зубцов задаётся отдельно. В производных по углу используйте радианы.
Статическая индуктивность \(L_{\text{ст}}=\Psi/i\) — наклон секущей из начала координат; это определение используется при \(\Psi(0)=0\) и \(i\ne0\). Динамическая индуктивность \(L_{\text{дин}}=\partial\Psi/\partial i\) — наклон касательной (дифференциальная, инкрементальная); положение и остальные токи фиксированы. Для линейной характеристики без смещения обе равны одной индуктивности \(L\). При постоянных магнитах и взаимном возбуждении сначала учитывают ненулевое смещение потокосцепления.
Положительные направления, оси и специальные индексы определяются в конкретной задаче. Единицу берите из условия или заголовка таблицы; перед подстановкой в формулу СИ переводите миллиметры в метры, миллигенри в генри и градусы в радианы.
Численные результаты обычно показаны с четырьмя значащими цифрами; для сравнения близких значений в таблицах сходимости сохранено больше цифр. Полные значения остаются в исходных расчётных данных. Малые остатки вычисления координат на осях показаны как ноль; малые физические величины сохраняются в показательной записи.
Постоянные k_x и k_y описывают ортотропную проводимость по координатным направлениям. Температурная кривая k(T) в описанной реализации задаёт изотропный материал; произвольное совмещение моделей нужно сверять с поддерживаемыми свойствами.
Между заданными точками зависимости k(T) руководство предусматривает линейную интерполяцию. За пределами диапазона сохраняется ближайшее крайнее значение коэффициента. Такое правило определяет поведение вычислительной модели, но не подтверждает свойства реального материала вне исследованного интервала. Поэтому при анализе решения необходимо контролировать диапазон полученных температур.
Вывод контрольного соотношения
Отделите сопротивление материала от внешнего охлаждения и оцените чувствительность потока к изменению k.
Одномерная стационарная теплопроводность с внешне заданной теплоотдачей. Radiation использует абсолютные температуры в четвёртой степени.
Последовательные тепловые сопротивления складываются: R=11\,\mathrm{К/Вт}.
Полный поток равен 50/11=4{,}5454545 Вт.
Температура охлаждаемой поверхности равна 300+Q\cdot 10=345{,}454545\,\mathrm{К}.
Границы применимости. Одномерная стационарная теплопроводность с внешне заданной теплоотдачей. Radiation использует абсолютные температуры в четвёртой степени.
Связь с опытом в FEMM
Замените постоянный k таблицей, совпадающей с k=10 в контрольном диапазоне; затем измените k и сравните сопротивления.
Сначала вычислите контрольное значение вручную, затем сформулируйте ожидаемый знак и зависимость от одного параметра. Повторите опыт в рабочей копии с изменением только этого параметра. Сохраняйте геометрию, единицы, глубину, материалы, источники и граничные условия вместе с результатом. Успешное завершение решателя означает решение введённой модели, но само по себе не доказывает её соответствие устройству.
Проверка и границы вывода
Постоянная таблица воспроизводит Q=4{,}5454545\,\mathrm{Вт}; все полученные температуры находятся в заявленном диапазоне применимости материала либо отмечены как выход за него.
h задаётся извне. Решатель теплопроводности не определяет скорость обтекания и не рассчитывает коэффициент теплоотдачи.
Идеальный контакт между слоями предполагает непрерывную температуру; контактное тепловое сопротивление нельзя получить простым переименованием материала.
Пояснения к официальному руководству
Подробности, которые следует учитывать при настройке и проверке модели. Все учебные пояснения.
Тепловые поля ввода, кривая \(T\)–\(k\) и суммарный поток проводника
Boundary задаёт Fixed Temperature, Heat Flux, Convection, Radiation, Periodic или Anti-Periodic. Перед выбором отличите температуру в К, поверхностную плотность теплового потока в Вт/м², коэффициент конвекции в Вт/(м²·К) и безразмерную излучательную способность. Излучение использует абсолютную температуру, поэтому градусы Цельсия напрямую в соответствующие поля не подставляют.
Constant Thermal Conductivity активирует kx,ky в Вт/(м·К). Thermal Conductivity Depends on Temperature отключает эти поля и открывает Edit T-k Curve: табличный k(T) изотропен. Точки сортируются по температуре, между ними применяется линейная интерполяция, вне диапазона удерживается ближайшее крайнее значение. Это правило экстраполяции программы не подтверждает материал за пределами эксперимента.
Brick — имя материала в текущем документе, а не тепловое граничное условие. Чтобы область стала кирпичом, её метке назначают эту запись и проверяют теплопроводность, тепловыделение и при необходимости объёмную теплоёмкость.
Обозначения на рисунке
Add
Добавить
Delete
Удалить
Modify
Изменить
Brick
Кирпич
Увеличить рисунок · Источник: Д. Микер, [33]. Подпись и пояснение переведены для курса.
Показан режим постоянной теплопроводности, поэтому таблица T–k неактивна. Тепловыделение задаётся в Вт/м³, а объёмная теплоёмкость — в МДж/(м³·К): массовую теплоёмкость в Дж/(кг·К) нельзя вводить сюда без умножения на плотность и пересчёта масштаба.
Обозначения на рисунке
Constant Thermal Conductivity
Постоянная теплопроводность
Thermal conductivity, x (W/m-K)
Теплопроводность по x, Вт/(м·К)
Thermal conductivity, y (W/m-K)
Теплопроводность по y, Вт/(м·К)
Edit T-k Curve...
Редактировать зависимость k(T)
Volume heat generation (W/m^3)
Объёмное тепловыделение, Вт/м³
Volumetric heat capacity (MJ/(m^3*K))
Объёмная теплоёмкость, МДж/(м³·К)
Увеличить рисунок · Источник: Д. Микер, [33]. Подпись и пояснение переведены для курса.
На самом изображении показана библиотека теплопроводности. В печатной подписи оригинала ошибочно указана Electrostatics Material Library; русская подпись исправлена по заголовку окна и контексту теплового раздела. Перенос из библиотеки создаёт запись материала, после чего её параметры и назначение области проверяют в модели.
Обозначения на рисунке
Library Materials
Материалы библиотеки
Model Materials
Материалы текущей модели
Heat Flow Material Library
Библиотека тепловых материалов
Metallic Solids
Твёрдые металлы
Nonmetallic Solids
Твёрдые неметаллические материалы
Liquid Metals
Жидкие металлы
Saturated Liquids
Насыщенные жидкости
Gases at 1 atm
Газы при давлении 1 атм
Исправлена явная ошибка подписи оригинала: на рисунке тепловая библиотека, а не электростатическая.
Увеличить рисунок · Источник: Д. Микер, [33]. Подпись и пояснение переведены для курса.
Import from file для T–k отдельно спрашивает порядок столбцов и единицы; температура в °C, °F или {}^\circR и проводимость в Btu/(ft·hr·F) либо kcal/(m·hr·K) преобразуются в К и Вт/(м·К). После импорта откройте Plot T-k Curve и проверьте хотя бы одну исходную пару вручную. Volume heat generation вводится в Вт/м³; Volumetric heat capacity — в МДж/(м³·К) и действует только в переходном расчёте с предыдущим решением и шагом времени.
Каждая строка связывает одну температуру с соответствующей теплопроводностью. В показанном примере чистой меди k убывает с температурой; это свойство конкретного набора данных, а не универсальный закон для всех материалов. После импорта проверьте порядок столбцов и единицы хотя бы одной пары.
Обозначения на рисунке
Copper, Pure
Чистая медь
T (K)
Температура, К
k (W/(m*K))
Теплопроводность, Вт/(м·К)
Import from file...
Импорт из файла
Plot T-k Curve
Построить график k(T)
Log Plot T-k Curve
Построить график k(T) в логарифмическом масштабе
Увеличить рисунок · Источник: Д. Микер, [33]. Подпись и пояснение переведены для курса.
Тепловой Conductor объединяет изотермические поверхности. Prescribed temperature задаёт общую T и вычисляет суммарный поток; Prescribed heat flux в этом диалоге означает заданный полный тепловой поток в Вт, в отличие от плотности потока на обычной границе. Вторая величина неактивна и вычисляется. Одинаковое имя на раздельных поверхностях связывает их одной изотермической идеализацией.
Проверьте на своей модели
Импортируйте T–k с одной температурой 20 °C и убедитесь, что она стала 293,15 К.
Сравните единицы Heat Flux в Boundary и Prescribed heat flux в Conductor; подтвердите полный поток через площадь.
Укажите, почему ненулевая теплоёмкость сама по себе не превращает стационарный расчёт в переходный.
Тепловое сопротивление пластины и конвекции — аналитический эталон
Дано.L=10\,\mathrm{мм}, k=10\,\mathrm{Вт/(м\cdot К)}, S=0{,}001\,\mathrm{м^2}; слева T=350\,\mathrm{К}; справа среда T_{\infty }=300\,\mathrm{К} и заданный h=100\,\mathrm{Вт/(м^2\cdot К)}.
Последовательные тепловые сопротивления складываются: R=11\,\mathrm{К/Вт}.
Полный поток равен 50/11=4{,}5454545 Вт.
Температура охлаждаемой поверхности равна 300+Q\cdot 10=345{,}454545\,\mathrm{К}.
Ответ:Q=4{,}5454545\,\mathrm{Вт}; T поверхности=345,5 К.
Проверка:h задаётся извне. Решатель теплопроводности не определяет скорость обтекания и не рассчитывает коэффициент теплоотдачи. Значения получены аналитически и сами по себе не подтверждают выполнения полевого расчёта FEMM.
Практикум в FEMM 5
Замените постоянный k таблицей, совпадающей с k=10 в контрольном диапазоне; затем измените k и сравните сопротивления.
В копии тепловой пластины задайте слева Fixed Temperature=350 К.
Создайте Convection boundary с T_{\infty }=300\,\mathrm{К} и h=100\,\mathrm{Вт/(м^2\cdot К)}; назначьте только правой стороне.
Проверьте, что старое фиксированное значение на правой стороне заменено; не задавайте там одновременно две независимые связи.
После решения измерьте температуру правой поверхности и полный поток; затем удвойте h и проверьте изменение сопротивления.
Ожидаемый результат.Q=4{,}5454545\,\mathrm{Вт}; T поверхности=345,5 К. Для FEMM сравнение действительно только при совпадении идеализации и граничных условий.
Как проверить.h задаётся извне. Решатель теплопроводности не определяет скорость обтекания и не рассчитывает коэффициент теплоотдачи. Прямое численное сравнение возможно только при совпадении всех исходных данных и постановки.
Обработка и анализ результатов
Остаётся ли расчёт внутри заданной зависимости теплопроводности от температуры?
Сохраните температурный диапазон, поток и сопротивление при постоянной и табличной теплопроводности.
Восстановите профиль или тепловое сопротивление из таблицы; сопоставьте суммарные источники и выходящие потоки.
С чем сравнить. Опора — закон Фурье и аналитический тепловой баланс из занятия. Теплоотдача и свойства материала являются исходными данными; одномерный предел применим только при согласованном пути тепла.
Итог работы. Проверка постоянного предела и объяснение изменения сопротивления материала.
Расчётные исследования по теме
Изучите физическую постановку, выполните программу во встроенном Python FEMM 5 и проверьте результаты по сохранённым моделям и независимым критериям.
Многослойный кабель позволяет проверить тепловую модель значительно строже, чем сравнением одной максимальной температуры. Для каждого пассивного слоя есть независимое аналитическое сопротивление, а через все слои проходит одна и та же мощность жилы. Учебная задача учит различать удельное тепловое сопротивление материала, теплопроводность и сопротивление целого цилиндрического слоя. Эти три величины имеют разные определения и не должны смешиваться в таблицах.
Два медных участка длиной по 45 мм соединены прослойкой длиной 10 мм, ширина 10 мм, глубина 1 м. Только прослойка выделяет 10 Вт на метр; внешние торцы поддерживаются при 20 °C.
Пластина 20 × 20 мм, глубина 1 мм; короткий контакт 60 °C на левой нижней стороне и контакт 20 °C на правой верхней стороне. Теплопроводности по осям 20 и 2 Вт/(м К).
Снимки интерфейса FEMM показывают геометрию, настройки и работу с результатами в FEMM 5. Откройте изображение, чтобы рассмотреть подписи и значения в полном размере.
Поперечная геометрия многослойного кабеля в модуле теплопроводности FEMM 5. Текстовые подписи скрыты в отображении, чтобы различались тонкие слои; зелёные точки — метки областей. Открыть в полном размереСохранённое решение «Тепловое сопротивление кабеля». Температура T, К, и изотермы многослойного кабеля; стационарная теплопроводность. Открыть в полном размереРабочая область модели «Тепловое узкое место контактной прослойки». Нативный препроцессор FEMM 5. Открыть в полном размереСохранённое решение «Тепловое узкое место контактной прослойки». Температура T, К, в контактной прослойке между медными теплоотводами. Изотермы скрыты для читаемости узкой нагретой области. Открыть в полном размереРабочая область модели «Ортотропный теплораспределитель аппарата». Нативный препроцессор FEMM 5. Открыть в полном размереСохранённое решение «Ортотропный теплораспределитель аппарата». Температура T, К, и изотермы ортотропного теплораспределителя; стационарная теплопроводность. Открыть в полном размере
Самостоятельная работа
Замените постоянный k таблицей, совпадающей с k=10 в контрольном диапазоне; затем измените k и сравните сопротивления. Сопоставьте исходный и изменённый режимы, сохраните исходные файлы, результаты и объяснение различия.
Критерий выполнения: Постоянная таблица воспроизводит Q=4{,}5454545\,\mathrm{Вт}; все полученные температуры находятся в заявленном диапазоне применимости материала либо отмечены как выход за него.
Пояснения и типичные ошибкиЧто происходит вне табличного диапазона k(T)?
Численно используется крайнее значение по описанию руководства; физическую допустимость такого продолжения автор модели проверяет отдельно.
Какое контрольное число и единица проверяют этот опыт?
Q=4{,}5454545\,\mathrm{Вт}; T поверхности=345,5 К.
При каких предположениях получена приведённая формула?
Одномерная стационарная теплопроводность с внешне заданной теплоотдачей. Radiation использует абсолютные температуры в четвёртой степени.
На что обратить внимание
Идеальный контакт между слоями предполагает непрерывную температуру; контактное тепловое сопротивление нельзя получить простым переименованием материала.
Одномерная стационарная теплопроводность с внешне заданной теплоотдачей. Radiation использует абсолютные температуры в четвёртой степени.
