Цель работы: изучение принципов постановки тепловой задачи на распределенном уровне и ее решения с помощью метода конечных элементов в среде ANSYS.
Теоретическая часть
Процессы теплопередачи происходят в пространстве и времени. Поэтому исследование теплопроводности сводится к изучению пространственно-временного изменения температуры, т. е. к нахождению зависимости , где , , — пространственные координаты в декартовой системе, — время.
Совокупность значений температуры во всех точках изучаемого пространства называется температурным полем. Различают стационарные и нестационарные температурные поля. Нестационарным температурным полем называется такое поле, температура которого изменяется не только в пространстве, но и с течением времени. Стационарным температурным полем называется такое поле, температура которого в любой его точке не изменяется во времени.
Для решения задач, связанных с нахождением температурного поля, необходимо иметь дифференциальное уравнение теплопроводности, которое дает зависимость между температурой, временем и координатами элементарного объема. В декартовой системе координат уравнение теплопроводности для изотропных материалов имеет вид:

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

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

    где — температура поверхности тела. В частном случае , то есть температура на поверхности постоянна на протяжении всего процесса теплообмена. Это может быть осуществлено при искусственном поддержании постоянной температуры или при особых условиях теплообмена (см. граничное условие третьего рода).
  2. Граничное условие второго рода состоит в задании плотности теплового потока для каждой точки поверхности тела как функции времени:

    где — нормаль к граничной поверхности. Простейший случай граничного условия второго рода состоит в постоянстве плотности теплового потока: . Такой случай теплообмена имеет место при нагревании тел в высокотемпературных печах, где передача тепла в основном происходит при помощи излучения по закону Стефана-Больцмана, когда температура тела значительно меньше температуры излучающих поверхностей.
  3. Граничное условие третьего рода характеризует закон конвективного теплообмена между поверхностью тела и окружающей средой (закон Ньютона). В этом случае количество тепла, передаваемого в единицу времени с единицы площади поверхности тела в окружающую среду с температурой , прямо пропорционально разности температур между поверхностью тела и окружающей средой:

Чтобы функция была решением тепловой задачи необходимо, чтобы она удовлетворяла дифференциальному уравнению, начальным и граничным условиям. В случае стационарной задачи начальные условия отсутствует.
По теореме единственности решения [1], если некоторая функция удовлетворяет выше указанным условиям, то она является единственным решением данной задачи.
Методы математической физики [1], в том числе метод разделения переменных Фурье, метод источников, операционные методы, методы интегральных преобразований, позволяют эффективно решать узкий круг задач с большим количеством допущений. В инженерных тепловых задачах имеют место конечные объемы сложной геометрической формы, произвольные граничные условия и материалы с нелинейными характеристиками. Поэтому для решения тепловых задач в инженерной практике применяются численные методы.
В пакете ANSYS, предлагаемом для дальнейшего изучения, реализован метод конечных элементов [2].
Список литературы
1. Лыков А. В. Теория теплопроводности. — М.: Высшая школа, 1967. – 598 с.
2. О.Зенкевич, К.Морган Конечные элементы и аппроксимация. Пер. с англ. – М.:Мир,1986 -318с.
3. Системы автоматизированного проектирования: В 9-ти кн. Кн. 4. Математические модели технических объектов: Учеб. пособие для втузов/В. А. Трудоношин, Н. В. Пивоварова; Под ред. И. П. Норенкова. — М.: Высш. шк., 1986. – 160 с.: ил.
4. Чигарев А. В., Кравчук А. С., Смалюк А. Ф. ANSYS для инженеров: Справ. пособие. М.: Машиностроение-1, 2004. 512 с.
5. http://www.fea.ru
Задание 1
Задание 2
Задание 3
Задание 4
Задание 5
Задание 6