9 Краткий обзор главы в данной главе рассматриваются основы теплового анализа с учетом фазового перехода. Обсуждаются следующие вопросы

Глава 9. Фазовый переход 9.1. Краткий обзор главы В данной главе рассматриваются основы теплового анализа с учетом фазового перехода. Обсуждаются следующие вопросы: основные определения; реализация процедуры анализа в программе ANSYS (на основе энтальпии); пример вычисления энтальпии по значению удельной теплоемкости; процедуры и рекомендации; пример расчета процесса затвердевания свинцовой отливки. 9.2. Основные определения Фаза представляет собой физически и химически однородную часть некоторого вещества. Изменение фазы является процессом, при котором часть вещества претерпевает изменение физических или химических свойств. Примеры различных фаз в одной системе: лед - вода - пар (различие физических свойств); феррит - перлит - аустенит в железе (химическое различие). Важными задачами при наличии фазовых переходов, которые решаются средствами конечно-элементного анализа, являются затвердевание жидкостей и плавление твердых тел. В частности, анализ плавления металла является реккурентной процедурой, включающей учет изменений фазового состояния. Скрытая теплота Выделение или поглощение теплоты, сопровождающее фазовый переход, рассматривается как скрытая теплота, поскольку изменение температуры при этом незначительно или не наблюдается вовсе. Фазовый переход в чистом веществе(без примесей) происходит без изменения температуры. Во всех других случаях изменение температуры незначительно. 9.3. Реализация процедуры анализа в программе ANSYS Тепловые задачи, включающие фазовые переходы, исследуются в нелинейном нестационарном тепловом анализе (ANTYPE,TRAN). Типичными расчетными величинами такого анализа являются следующими: изменение узловых температур во времени; продолжительность фазового перехода; тепловые характеристики, например, формовочной опоки при литье металла; положение границы раздела между жидкой и твердой фазами. Основное отличие обычного теплового нестационарного анализа от теплового анализа при наличии фазового перехода состоит в том, что в последнем случае необходимо учитывать скрытую теплоту. 9.3.1. Фазовый переход на основе учета энтальпии В программе ANSYS скрытая теплота не учитывается явным образом посредством конечно-элементной формулировки задачи. При этом, однако, скрытая теплота вводится как свойство материала - энтальпия (ENTH). Энтальпия вещества H связана с плотностью , удельной теплоемкостью c и температурой T следующим образом: . Энтальпия имеет размерность "количество теплоты" / "объем". На изменение фазы вещества указывает резкое изменением энтальпии в некотором диапазоне температур. На приведенном графике разность температур Tl - Ts является достаточно малой величиной. Для чистого материала эта разность равна нулю. 9.3.2. Задание свойств материала Для обычного нестационарного анализа требуются следующие свойства материала: теплопроводность (в трех направлениях KXX, KYY, KZZ); плотность (DENS); удельная теплоемкость (C). При анализе с учетом фазового перехода используются теплопроводность (KXX, KYY, KZZ); энтальпия (ENTH). В общем случае, когда любое из этих свойств зависит от температуры, анализ оказывается нелинейным. Анализ фазовых переходов всегда нелинеен, так как энтальпия существенно зависит от температуры. Если для какого-то материала нет данных по энтальпии, то для определения зависимости энтальпии от температуры рекомендуется использовать значения плотности , удельной теплоемкости c и данные по скрытой теплоте. Пример Вычисление величины энтальпии по значению удельной теплоемкости. Дано: - плотность материала (постоянная); Cs - удельная теплоемкость твердого материала; Cl - удельная теплоемкость жидкой фазы; L - скрытая теплота. Используется соотношение . Допущение: если температуры твердой и жидкой фаз (Ts,Tl) не известны, задается приемлемый диапазон температур для переходной зоны. Энтальпия является более гладкой функцией температуры по сравнению с удельной теплоемкостью. Следовательно, в задачах с учетом фазовых переходов целесообразно использовать энтальпию - в этом случае точность решения ухудшается незначительно, даже если шаг интегрирования по времени оказывается сопоставимым с шириной переходной зоны. Преобразование удельной теплоемкости в энтальпию осуществляется следующим образом: при при при T H(расчетная) T1 0 T2 где 9.4. Рекомендации при проведении анализа Используйте энтальпию в качестве свойства материала при расчетах скрытой теплоты. Задавайте небольшое значение минимального шага по времени (команда DELTIM) и используйте автоматический выбор шага по времени (команда AUTOTS), чтобы программа могла корректировать шаг до начала фазового перехода, в течение перехода и после его окончания. Если возможно, используйте элементы низшего порядка, такие как PLANE55 или SOLID70. Если требуются элементы высокого порядка, выбирайте вариант диагональной матрицы удельной теплоемкости, используя соответствующий ключ KEYOPT. Выбирайте обратный метод Эйлера для процедуры интегрирования по времени. Это достигается установкой THETA = 1 в команде TINTP. Имейте в виду, что использование опции линейного поиска (команда LNSRCH) улучшает сходимость при расчетах фазового перехода. Пример Затвердевание свинцовой отливки. Начальная температура Изоляция Основная цель данного примера заключается в определении границы раздела фаз при использовании энтальпии для учета скрытой теплоты. Кроме того, предполагается получить графики изменения температуры во времени и линии равных температур. 9.4.1. Принятые допущения Предполагается идеальный тепловой контакт на границе отливка-песок, т.е. предполагается непрерывность температур. Слиток и опока имеют достаточную длину в направлении оси z, поэтому можно ограничиться двумерным анализом. Изменением объема в процессе затвердевания пренебрегается. Естественная конвекция расплава не учитывается. Предполагается, что форма заливается мгновенно, и в момент времени t=0 как опока, так и отливка имеют одинаковую начальную температуру 400C. 9.4.2. Свойства материала отливки Теплопроводность материала задается графиком и таблицей значений. Т К (С) (Вт/м-С) 704 14.9 Значения свойств материалов, используемые в этом примере, взяты из работы Хольмана “Теплопередача”, 5-ое изд.,McGraw-Hill, 1981 и справочника Baumeister,Marks “Руководство для инженеров-механиков”, 7-ое изд., McGraw-Hill, 1967. Энтальпия: В данном примере зависимость энтальпии материала отливки от температуры рассчитана на основе данных по удельной теплоемкости, скрытой теплоте, плотности и температуре точки плавления. Т Н (С) (Дж/м3) 700 13.763*108 Удельная теплоемкость Cs = 130, Cl = 159 Дж/кг-C Плотность = 11000 кг/м3 Скрытая теплота L = 23300 Дж/кг, Hl - Hs = 2.563.108 Дж/м3 Температура плавления Tm = 327C Принятые температуры перехода Ts = 326.5, Tl = 327.5C Ширина "переходной" зоны Tl-Ts = 1 C 9.4.3. Свойства формовочного песка Теплопроводность k = 0.325 Вт/м-C Удельная теплоемкость c = 0.816 Дж/кг-C Плотность = 1520 кг/м3 Обратите внимание: фазовых изменений в этом материале не предполагается, и энтальпия не используется. 9.4.4. Команды задания свойств материалов mptemp,1,0,100,200,300,371,704 * Задание значений KXX и ENTH для отливки mpdata,kxx,1,1,35.1,33.4,31.5,29. 8,16.1,14.9 mptemp,1,0,326.5,327.5,700,800,900 mpdata,enth,1,1,0,4.669e8,7.248e8,13.763e8 mp,kxx,2,0.325 * Задание значений KXX, DENS и C для песка mpplot,enth,1,0,780,0,1.6e9 * Построение графика энтальпии Материал Nо.1 - свинец Материал Nо.2 - песок Значения энтальпии для температур вне заданного интервала получены экстраполяцией. 9.4.5. Команды для модели и граничных условий et,1,55 * Тепловые твердотельные элементы PLANE55 rectng,0,0.06,0,0.048 * Задание поверхностей для формы и отливки ptxy,0.014,0.012,0.046,0.012,.048,0.036,0.012,0.036 poly nsel,all Сетка элементов с указанием номеров материалов: 9.4.6. Команды задания опций шага нагружения kbc,l * Ступенчатое изменение нагрузки time,10800 * Длительность всего процесса 10800 секунд solve * Инициализация решения 9.5. Итоговая информация Конечные элементы: 401 элемент PLANE55 (108 для отливки, 293 для формовочного песка). Опции: переходный анализ (0 - 10800 с); полный метод Ньютона-Рафсона (по умолчанию); обратный метод интегрирования Эйлера (THETA = 1 в команде TINTP ). Используются процедуры: автоматического выбора шага по времени (пределы шага по времени: min=0.1 с, max=120 с); линейного поиска. Предполагаемые параметры решения: число возможных шагов решения 108 000; количество равновесных итераций на шаг решения равно 25 (по умолчанию); общее число возможных итераций 2 700 000. Действительные параметры решения: 102 шага решения; общее количество итераций 294; две бисекции шага по времени. 9.5.1. Определение границы раздела фаз Напомним: температура затвердевания/плавления свинца равна 327С. Чтобы определить границы раздела фаз, используется команда /CVAL. /postl * Начало постпроцессорной обработки /plopts,minm,off esel,s,mat,,1 /cva1,1,327,400 * Определение границы раздела фаз в свинце set,first aplot * Совмещение графиков результатов В команде /CVAL число 327 означает температуру начала превращения жидкого вещества в твердое. Предполагается, что при 400С свинец находится в жидком состоянии. /noerase plnsol, temp * Начало графического вывода *get, time, active, set, time Для графического вывода результатов используется постпроцессор POST26. /post26 nso1,2,152,temp,,NODE 152 * Графический вывод для нескольких узлов plvar,2,3,4,5,6 9.6. График зависимости температуры от времени 9.7. Полная распечатка для процесса затвердевания /batch,list * Пакетный режим работы /show,ingot,grph * Запись графиков в файл ingot.grph /title, Solidification of a Lead Ingot et,1,55 * Тепловые твердотельные элементы PLANE55 mpdata,kxx,1,1,35.1,33.4,31.5,29.8,16.1,14.9 mp,kxx,2,0.325 * Задание свойств формовочного песка ptxy,0.014,0.012,0.046,0.012,.048,0.036,0.012,0.036 poly /contour, 1, 16, 25, 25 * Задание контуров одинаковых температур finish Температурное поле модели на 3100 секунде Температурное поле модели на 10660 секунде

Похожие:

	Основные положения и обзор возможностей программы В этой главе кратко обсуждаются основные теоретические аспекты процесса теплопередачи и конкретная их реализация в программе ansys....		В дисциплине рассматриваются следующие вопросы Дисциплина имеет практическую направленность на решение широкого спектра прикладных вопросов на всех уровнях иерархии относительно...
	Экзаменационные билеты 4 3 Вопросы к 1 главе 9 4 Вопросы к 2 главе 12 5 Вопросы к 3 главе 17 6 Вопросы к 4 главе 22 7 Вопросы к 5 главе 27 8 Вопросы к 6 главе 32 Контрольные вопросы курсу «Проектирование процесса оказания услуг» «Проектирование процесса оказания услуг» студентам специальности «Сервис» специализация «Сервис электронных систем безопасности (информации,...		Бюллетень осжд. 2010. №3 С. 19-21. Акционерное общество «Чешские железные дороги» Выш школа, 2010. – 600 с. – Рассматриваются теоретические основы и практические методики бухгалтерского учета, финансового анализа...
	Программа курсов повышения квалификации «современные методы рентгеновского фазового анализа» Цель дисциплины: формирование навыков самостоятельного проведения рентгеновского фазового анализа для работников исследовательских...		Universitatea libera internationala din moldova факультет экономических знаний данилова Анастасия, Мастерант Тема: Процессы обучения (Лабораторная работа) Некоторые из этих проблем рассматриваются в других главах настоящего тома, а краткий перечень работ, посвященных этим проблемам,...
	Пирамида общественного мнения в электоральном «зеркале» «данной» в материалах исследований общественного мнения, а кроме того — оценки эффективности теоретических средств его анализа. (Вопросы...		Вестн. Самар. Гос. Техн. Ун-та. Сер. Технические науки. 2011. №4 (32) Краткие сообщения Рассмотрен упрощенный способ решения тепловой задачи нагрева контактной системы выключателя с учетом фазового перехода
	3 экспертные методы в главе обсуждаются основные вопросы теории и практики экспертных оценок, в том числе связанные с типовыми стадиями экспертного опроса Рассмотрены основные идеи современной теории измерений, метода согласования кластеризованных ранжировок и ряда других математических...		Основы защиты населения В настоящей главе рассматриваются основные способы защиты населения от оружия массового поражения, другие вопросы защиты будут рассмотрены...