Теория и практические методы плавки в холодном тигле

--PAGE_BREAK--1.1.1     Преимущества индукционной плавки


Индукционная плавка в холодных тиглях отличает­ся от дуговой и плавки электронным пучком следую­щими преимуществами:

1.  При индукционной плавке возможно  интенсивное перемешивание расплава без использования специальных устройств. При индукционном нагреве тепло выде­ляется в значительном объеме расплава или в поверх­ностном слое, объема. Это создает благоприятные усло­вия для конвективного перемешивания и получения рав­номерной температуры по объему. Но особенно важно, что при индукционном нагреве создастся сильное при­нудительное перемешивание расплава вследствие элек­тродинамического взаимодействия токов, индуктирован­ных в расплаве, с током индуктора.

2.  Используя перемешивание, расплав можно пере­греть во всем объеме тигля. Предел перегрева опреде­ляется только мощностью источника тока и  устойчивостью расплава при высоких температурах.

3.  Индукционную плавку в холодных тиглях можно проводить в любой среде, начиная от глубокого ваку­ума, необходимого для рафинирования расплава от ле­тучих примесей, кончая атмосферой повышенного дав­ления, необходимой дли предотвращения термической диссоциации соединений с летучей компонентой.

4.  Поверхность расплава при индукционной плавке не закрыта электродом, как в дуговой печи, поэтому возможна лучшая очистка расплава от летучих приме сей и газов. В сочетании с регулируемым нагревом и интенсивным перемешиванием наличие свободной по­верхности обеспечивает   условия рафинирования при нормальном давлении,  аналогичные электроннолучевой плавке, но без перегревов и испарения расплава.

5.  Процесс индукционного нагрева надежен и устой­чив, подаваемую мощность можно легко регулировать в широких пределах. В установках для индукционного нагрева отсутствуют такие недолговечные элементы, как катод в установках электроннолучевого нагрева.

6.  Вследствие того, что токи, наведенные в распла­ве при наличии гарнисажа, и токи, наведенные в ме­таллическом тигле, не являются продолжением друг друга, как это происходит в дуговых печах, дуги меж­ду расплавом и тиглем не возникает, и нет опасности на­рушения стенки тигля.  

Схема с цилиндрическим индуктором также имеет ряд достоинств:

-                   всесторонний нагрев с перемешиванием расплава позволяет получить глубокую ванну с равномерным распределением температуры по высоте и радиусу

-                   индуктор имеющий высокое напряжение вынесен из зоны термоионизированных газов над расплавом, что снижает вероятность электрических пробоев и повышает надежность работы печи

Особенностью плавки в холодных тиглях является наличие гарнисажа. В металлургии гарнисажем называется особый слой на границе расплава и тигля, предохраняющий тигель от разъедания. Плавка с гарнисажем гарантирует абсолютную чистоту расплава, так как он отделен от возможного источника загрязнения, то есть стенки тигля, кристаллическим слоем. Именно эта особенность привлекает к ней повышенный интерес в металлургии тугоплавких материалов.

Именно за счет схемы с цилиндрическим индуктором и соответственно наличия естественного перемешивания, ванна расплава является не стационарной. Для учета конвекционных потоков в модели решается гидродинамическая задача. Наличие гарнисажа также отразится на графике распределения тепловых полей.




1.2       Коммерческий пакет программ ANSYS


1.2.1     Описание пакета «ANSYS»


Универсальная программа ANSYSпозволяет выполнять практически любой вид анализа (расчета) с помощью метода конечных элементов. МКЭ означает, что задача может решаться в любой геометрически сложной области. Универсальность программы заключается в том, что она может быть использована для решения задач, связанных со многими разделами науки и техники, такими, как динамика и прочность машин, электротехника, электромагнетизм, электроника, теплофизика, гидродинамика, газовая динамика, аэродинамика, биомеханика и т. д.

Программа «ANSYS» предназначена для численного решения уравнений (систем уравнений) с частными производными.

Следовательно, численные решением является таблица чисел или множество таких таблиц, соответствующих дискретным моментам времени.

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

— геометрию расчетной области,

— значения коэффициентов уравнений или теплофизические свойства среды, которая находится в расчетной области,

— начальные условия,

— граничные условия.
    продолжение
--PAGE_BREAK--1.2.2      Задание свойств материалов (среды)


Для большинства элементов необходимо задавать свойства материалов (среды). Свойства материалов могут быть:

— линейными или нелинейными,

— изотропными или анизотропными,

— зависимыми или независимыми от температуры.

Естественно, в модели определяются только те свойства, которые необходимы для решения задачи.

Часто теплофизические свойства материалов зависят от температуры. Эту зависимость можно задавать в виде полинома или таблично.

Полиномы могут быть линейными, квадратичными, кубическими или четвертой степени:

Свойство = С0 + С1Т + С1Т2 + С1Т3 + С1Т4.                                             (2)

Если определено только С0, свойство постоянно, если определены С0и С2, свойство линейно зависит от температуры и т.д.
1.2.3     Задание свойств жидкостей при решении задач гидродинамики


При решении задач гидродинамики в основном задают значения следующих теплофизических свойств:

— плотность,

— динамическая вязкость,

— теплопроводность,

— удельная теплоемкость.

При этом при задании свойств могут быть определены многие дополнительные параметры жидкостей.

Свойства жидкостей также могут задаваться как полиномом, так и табличным методом. При этом задание свойств табличным методом дает меньшую относительную погрешность, чем при расчете полинома.

Кроме выбора расчета плотности, вязкости, теплопроводности и удельной теплоемкости имеется возможность пересчитывать их от одной глобальной итерации к другой при поиске решения, тем самым улучшая сходимость решения.
1.2.4      Создание геометрической модели


Основной целью конечно-элементного анализа является матема­тическое воссоздание поведения реальной технической системы. Дру­гими словами, анализ должен быть точной математической моделью физического прототипа. В широком смысле модель включает в себя все узлы, элементы, свойства материалов,, константы, граничные ус­ловия и другие особенности, которые используются для представле­ния физической системы.

В принятой для программы ANSYSтерминологии понятие генера­ция модели подразумевает создание узлов и элементов, которые ото­бражают пространственную конфигурацию и связи реальной системы. То есть генерация (создание) модели означает задание геометриче­ской конфигурации узлов и элементов модели. Программа ANSYSпред­лагает пользователю несколько различных методов построения моде­ли, допуская при ее создании использование комбинации способов.

На практике задают не только геометрию модели, то есть не только разбивают пространство на области с разными свойствами, но еще и вводят более мелкое разбиение – сетку. Сетка вводится для того, чтобы можно было использовать метод конечных элементов. Данный метод позволяет даже самые сложные процессы приближенно представить, как совокупность множества простых, развивающихся по линейному закону. Пример разбиения на такую сетку представлен ниже на рисунке 1.



Рис 1. Пример разбиения области на сетку

2.       Практическая глава
2.1       Постановка задачи


Пакет «ANSYS» представляет собой программу предназначенную для моделирования различных типов задач от механических до тепловых с последующей их визуализацией и возможной обработкой результатов. В данном пакете можно решать одно-, двух- и трехмерные задачи.

На данном этапе моделирования решалась двухмерная электротепловая задача. Была создана модель, включающая в себя цилиндрический холодный тигель с расплавом высотой 200 мм, дно и двухвитковый индуктор высотой 200 мм. Так как система является цилиндрической и имеет ось симметрии, то рассматривается только половина данной системы. Созданная модель служит для расчета и визуализации температурного распределения, распределения источников теплоты в расплаве. До начала моделирования возможно задание различного времени нагрева, подаваемого напряжения, тока индуктора или плотности тока на индукторе, и изменение параметров расплава. Также созданная модель имеет возможность перерасчитывать в процессе моделирования распределение источников теплоты через заданный промежуток времени для получения более точного результата. В модели отсутствует стартовый элемент, помещаемый в расплав, поэтому время нагрева расплава до необходимой температуры существенно больше реального.

Следующим шагом по модернизации существующей модели будет включение в нее гидравлического расчета. Без этого расплав представляет собой твердое тело. В тугоплавких стеклах движение расплава составляет 2-5 мм/с, поэтому, при учете гидродинамики, температура, самой горячей области расплава, понизится, а в других областях возрастет.

Далее представлены результаты тестового расчета модели. При этом использовались следующие основные исходные данные:

— время нагрева: 10000 сек. Это примерное время достижения установившегося состояния температуры самой горячей зоны расплава;

— напряжение на индукторе: 1535 В;

— время перерасчета источников теплоты: 100 сек. Источники теплоты перерасчитываются каждые 100 секунд;

— начальная температура: 20 ˚С;

— максимальный шаг по времени: 10 сек. Данный шаг по времени обеспечивает максимальную сходимость результатов расчетов

— частота: 1,76·106 Гц;

— плотность: 2700 кг/м3;

— электропроводность: 5,9·10-2 Ом∙м;

— теплоемкость: 1700 Дж/(кг∙К);

— теплопроводность: 1,2 Вт/(м∙К);

Также заданы:

— величина тепловых потерь в боковую стенку: 20 Вт/см2;

— величина тепловых потерь в дно: 5 Вт/см2;

— величина тепловых потерь с зеркала расплава: 10 Вт/см2;

При этом также возможно:

— изменять значения потерь в бок и дно тигля;

— изменять плотность, электропроводность, теплоемкость и теплопроводность расплава;

— выбирать между решением электрической задачи по току индуктора, по напряжению на индукторе и по плотности тока на индукторе;

— выбирать начальное время.

В разработке электрогидродинамической модели не учитывается холодный тигель (ХТ). Для учета влияния ХТ в модели используется так называемый метод приближения, который заключается в том, что берется площадь занимаемая трубками ХТ и данное кольцо полученное из трубок тигля вычитается из площади сечения индуктора. В модели индуктор приближен к расплаву на ту площадь, которую занимали трубки ХТ. Это позволяет провести расчеты, сымитировать ХТ.



    продолжение
--PAGE_BREAK--2.2       Опыт проведенный на кафедре ЭТПТ

2.2.1     Температурные измерения


Проведение теста заключалось в  измерении профиля температур с помощью 8 термопар. Измерение производилось при 3х различных режимах работы генератора. Опускание блока термопар вглубь расплава производилось не плавно, а ступенчато с шагом 10мм при заглублении 0 – 150 мм и с шагом 5 мм при опускании на глубину 150 – 180 мм в первом режиме. Так как термопары имели кварцевые чехлы, то после второго режима некоторые из них вышли из строя. В результате высоких температурных перепадов произошло растрескивание кварцевых чехлов и для получения результатов было решено увеличить шаг: во втором режиме шаг равен 10 мм при 0 – 180 мм и 10 — 20 мм при третьем режиме работы генератора. Измерения производились только при опускании термопар, из-за высоких температур нахождение долгое время в расплаве термопар было невозможно.

Как уже упоминалось ранее, для проверки правильности и точности модели модель разрабатывалась под уже существующую установку, и процессы смоделированы по опыту, который проводился на кафедре ЭТПТ. Целью опыта было по измерение профилей температур в расплаве стекла для трех режимов генератора с помощью высокотемпературных термопар.

Измерение температуры расплава и получение картины теплового поля в расплаве производилось с помощью 8 термопар. Система крепления и схема самой термопары рассмотрены ниже.
Схема термопары и ее геометрические размеры:



Рис 2. Схема термопары и ее геометрические размеры


Расположение термопар:




Рис.3 Рисунок термопар с кварцевым покрытием (чехлом)

1 – подвижная трубка; 2 – соединение трубки с крепление блока термопар; 3 – крепёжный пластина;

4 – термопары; 5 – скобы крепления; 6 – вертикальный алюминиевый крепеж.

На рис № 4 представлено расположения блока термопар. Так как представленная индукционная система является осесимметричной, то отпадает необходимость установки термопар по всему поперечному сечению холодного тигля.



Рис 4. Вид смонтированной системы крепления и перемещения Pt-Pt/Rdтермопар
2.2.2     Результаты эксперимента


В результате проведенного эксперимента были получены данные температурных полей, которые представлены в виде  сводных таблиц результатов измерений, которые показывают значение температуры каждой термопары на определенной глубине в расплаве при различных режимах плавки.

Показания действующих термопар в первом режиме представлены в таблице 1.

Таблица 1. Показания действующих термопар в первом режиме


Показания действующих термопар во втором режиме представлены в таблице 2.

Таблица 2. Показания действующих термопар во втором режиме


Показания действующих термопар в третьем режиме представлены на таблице 3. В таблице отсутствие термопар с номерами 1 и 2 связано с их выходом из строя, следовательно, данные с них в 3 режиме не были получены.

Таблица 3. Показания действующих термопар в третьем режиме



В таблицах 1-3 помечены курсивом значения выпадающие из допустимого диапазона температур это может быть связано с электромагнитными наводками на спае термопар при соприкосновении с расплавом; отсосом тепла в спае термопар и кварцевый чехол с зеркала расплава. Жирным показаны оптимальные значения температур. Прочерки в таблицах это отсутствие данных с термопары на данном этапе.

В ходе эксперимента показания термопар записывались в реальном времени через цифровую система сбора данных.

На рисунках ниже представлены графики измерения термопарами полей температур в холодном тигле в реальном времени. Каждая «полка» (горизонт. фронт) представляет собой один шаг.

 


 


Рис 5. Показания измерений термопар в реальном времени при первом режиме
   


Рис 6. Показания измерений термопар в реальном времени при втором режиме

   


Рис 7. Показания измерений термопар в реальном времени при третьем режиме

На рисунке 6 и 7 всплески (резкие пики) связаны с выходом из стороя термопар и отсутствием получаемого с них сигнала.



    продолжение
--PAGE_BREAK--