Математическое моделирование воздействия потока энергии на материалы с учетом фазовых переходов

Введение 5
1. Моделирование воздействия потока энергии на материалы 6
1.1. Нестационарное уравнение теплопроводности 10
1.2. Аналитические решения 11
1.3. Фазовые переходы 13
1.4. Метод прогонки 20
1.5. Решение уравнение теплопроводности в среде Matlab. 24
2. Численное исследование нестационарного температурного поля при воздействии потока энергии на материалы с учетом фазовых переходов 27
Заключение 36
Список литературы 37
Приложения 38

Весь текст будет доступен после покупки

Для модификации структуры приповерхностного слоя изделий и материалов с целью повышения их эксплуатационных характеристик применяется электронно-пучковая импульсная обработка. В настоящее время это является одним из наиболее перспективных способов упрочнения и изменения характеристик приповерхностного слоя материалов. При электронно-пучковой обработке достигаются сверхвысокие скорости (108-1010 К/с) нагрева до температур плавления, последующее охлаждения тонкого приповерхностного слоя материала (10-7-10-6 м) происходит за весьма малые (10-6-10-3 с) времена воздействия высоких температур, формирование предельных градиентов температуры (до107 …108 К/м) обеспечивает охлаждение приповерхностного слоя за счет теплоотвода в основной объем материала со скоростью (104 …109 К/с). Все это создает условия для образования в приповерхностном слое аморфной, нано- и субмикрокристаллической структуры [1-3].

Весь текст будет доступен после покупки

1. Моделирование воздействия потока энергии на материалы
Изучение любого физического явления сводится к установлению зависимости между величинами, характеризующими это явление. Для сложных физических процессов, в которых определяющие величины могут существенно изменяться в пространстве и времени, установить зависимость между величинами очень трудно. В этих случаях на помощь приходит метод математической физики, в котором ограничивается промежуток времени и из всего пространства рассматривается лишь элементарный объем. Это позволяет в пределах элементарного объема и выбранного малого отрезка времени пренебречь изменением некоторых величин, характеризующих процесс, и существенно упростить зависимость.
Выбранные таким образом элементарный объем и элементарный промежуток времени , в пределах которых рассматривается изучаемый процесс, с математической точки зрения являются величинами бесконечно малыми, а с физической точки зрения - величинами еще достаточно большими, чтобы в их пределах можно было игнорировать дискретное строенное среды и рассматриваться ее как континуум (сплошную). При решении задач, связанных с нахождением температурного поля, необходимо иметь дифференциальное уравнение теплопроводности.
Для облегчения вывода этого дифференциального уравнения сделаем следующие допущения:
+тело однородно и изотропно;
+физические параметры постоянны;
+деформация рассматриваемого объема, связанная с изменением температуры, является очень малой величиной по сравнению с сами объектом;
+внутренние источники теплоты в теле, которые в общем случае могут быть заданы как , распределены равномерно.
В основу вывода дифференциального уравнения теплопроводности положен закон сохранения энергии, который в рассматриваемом случае может быть сформулирован следующим образом: Количество теплоты dQ, введенное в элементарный объем извне за время вследствие теплопроводности, а также от внутренних источников, равно изменению внутренней энергии или энтальпии вещества (в зависимости от рассмотрения изохорического или изобарического процесса), содержащегося в элементарном объеме:
(1)
где - количество теплоты, Дж, введенное в элементарный объем путем теплопроводности за время ; - количество теплоты, Дж, которое за время выделилось в элементарном объеме ; - изменение внутренней энергии или энтальпии вещества, содержащегося в
Рис. 1.
элементарном объеме , за время .
Для нахождения составляющих уравнения выделим в теле элементарный параллелепипед со сторонами , (рис. 1). Параллелепипед расположен так, чтобы его грани были параллельны соответствующим координатным плоскостям.
Количество теплоты, которое подводится к граням элементарного объема за время в направлении осей , обозначим соответственно .
Количество теплоты, которое будет отводиться через противоположные грани в тех же направлениях, обозначим соответственно . Количество теплоты, подведенное к грани в направлении оси за время , составляет , где
(2)
- проекция плотности теплового потока на направление нормали к указанной грани. Количество теплоты, отведенное через противоположную грань элементарного параллелепипеда в направлении оси , запишется как
(3)
Разница количеств теплоты, подведенных к элементарному параллелепипеду и отведенных от него за время в направлении оси , представляет собой количество теплоты
(4)
Или
(5)
Функция является непрерывной в рассматриваемом интервале и может быть разложена в ряд Тейлора:

Если ограничиться двумя первыми членами ряда, то уравнение (4) запишется в виде
(6)
Аналогичным образом можно найти количество теплоты, подводимое к элементарному объему и в направлениях двух других координатных осей и .
, (7)
. (8)
Полное количество теплоты , подведенное теплопроводностью к рассматриваемому объему, будет равно:
(9)
С другой стороны, согласно закону сохранения энергии:
, (10)
Сопоставив выражения для теплоты и произведя сокращения, получим дифференциальное уравнение теплопроводности:
, (11)
Или в сокращенной записи:
. (12)
Величина

- оператор Лапласа .

Весь текст будет доступен после покупки

1. Иванов IO.Ф., Коваль К.Н., Колубаева Ю.А., Тересов А.Д. и др. Модификация поверхности технически чистого титана электронным пучком с субмиллисекундной длительностью // Изв. вузов. Физика. 2008. № 11/3. С. 17–21.
2. Денисова Ю.А., Иванов Ю.Ф., Иванова О.В., Иконникова И.А., Коваль Р.П., Крысина О.В., Петрикова Е.А., Тересов А.Д., Шугуров В.В. Эволюция структуры поверхностного слоя стали, подвергнутого электронно-ионно-плазменным методам обработки. – Томск: Изд-во НТЛ, 2016. – 302 с.

3. Ю.И.Иванов, С.В.Карпий и др. // Структура, фазовый состав и свойства титана после электровзрывного легировании и электронно-пучковой обработки. Серия «Фундаментальные проблемы современного материаловедения» Новокузнецк, 2010. - С.123-138.

4. А. В. Лыков. Теория теплопроводности – изд «Высшая Школа» Москва 1967.592 с.
5. Исаченко В. П. и др. Теплопередача. Учебник для вузов, изд. 3-е, перераб. и доп. М, «Энергня», 1975. 488с.

Весь текст будет доступен после покупки