г;г:ггг£ r? я яим/гсгшг
(35). 2005 -
Ш
ИТЕЙНОЕ ПРОИЗВОДСТВО;
77z<? worked out mathematical model confirms the presence of non-symmetrical front edge of crystallization by cut of casting and allows to estimate qualitativly the dynamics of the casting shell growth in crystallizer in course of time.
ишшшт
_Ж
Е. Б. ДЕМЧЕНКО, БИТУ
КИНЕТИКА ЗАТВЕРДЕВАНИЯ ОТЛИВКИ ПРИ ГОРИЗОНТАЛЬНОМ НЕПРЕРЫВНОМ ЛИТЬЕ
УДК 621.74.047
Известно, что при горизонтальном непрерывном литье в процессе формирования корки отливки в верхней зоне кристаллизатора образуется более тонкая корка, чем в нижней зоне [1—3]. Такое различие обусловлено несколькими факторами, определяющими дальнейший ход процесса затвердевания отливки. Прежде всего это усадка отливки, протекающая при ее формировании и последующем охлажден™. В результате усадки уменьшается размер отливки, за счет сил тяжести она прижимается к нижней зоне кристаллизатора, образуя в верхней зоне значительный по величине газовый зазор. Теплоотвод в верхней зоне ухудшается, корка перегревается и ее рост замедляется. Другой, не менее значительный фактор — это термоконвективные потоки горячего расплава, размывающие верхние участки корки и снижающие скорость затвердевания [4, 5]. Термоконвективные потоки обусловливают неравномерное распределение теплоты перегрева в полости кристаллизатора. В результате фронт кристаллизации смещается вверх и становится несимметричен относительно оси вытяжки отливки, верхние участки затвердевшей корки на 5,0-20% тоньше нижних, кроме того, имеет место запаздывание затвердевания в верхней зоне кристаллизатора (рис. 1) [2]. Все это является причинами нарушения стабильности литья и возникновения прорывов металла на выходе из кристаллизатора.
Если влияние термоконвективных потоков можно снизить, используя разделительные устройства [4], то уменьшить влияние усадки отливки на несимметричность фронта кристаллизации практически нереально. В таких случаях необходимо учитывать неравномерность фронта затвердевания в математических моделях при расчетах кинетики кристаллизации отливки и выборе оптимальных параметров литья.
Покажем это на примере расчета процесса затвердевания сплошной цилиндрической отливки радиусом К в комбинированном кристаллизаторе длиной стальной втулкой и рабочим графитовым вкладышем.
Рис. 1. Запаздывание затвердевания ДА в верхней зоне кристаллизатора сплошной цилиндрической отливки диаметром 0.08 м. чугун С420
Составим для участка длиной Д/.( (рис. 2). ограниченного сечениями / и /+1, уравнение теплового баланса системы «жидкий расплав— твердая корка—кристаллизатор»:
И,ж ■ QLm Й'1,1
Рис.2. Схема к расчету кинетики кристаллизации сплошной цилиндрической отливки при горизонтальном лшье
Раскроем значения каждого слагаемого кроме отведенной теплоты перегрева 0, пер. неравномерность распределения которой учтем отдельно.
Количество тепла Огл, внесенное жидким расплавом через сечение г. запишем как:
/;Г7Г.Т:Г ГГШМЧРТГМ
- 3 (35), 2005
/29
6:.Ж =л\\лУ1(лл-л/)л('иэф+С1лкр). чере з сечение /+1;
а+1.ж =7лИлфУ 1 (л?-лш)л("эф +С17'кр) ,
где -передняя скорость литья: " $ф - эффек-
тивная удельная теплота кристаллизации: "
удельный вес материала отливки: с. - удельная теплоемкость твердого расплава; Ткр температу ра кристаллизации расплава.
Количество тепла 0,л, внесенное в единиц) времени твердой коркой через сечение /:
д|,т ~ ~'ЛсрУ|Л| (2-11 ~ )*ЭД|',ср ,
через сечение /+1:
где Т. ср — средняя температу р>а затвердевшей корки в сечении /; Т1+] ср — средняя температура корки в сечении г+1.'
Количество тепла отведенного через
кристаллизатор к охлаждающей среде:
2; охл -л'"А'"ср4/.с(Л.
где д1ср средний тепловой поток на участке А/:е Л г. - время формирования твердой корки отливки на участке Д/.. (
Для учета неравномерности распределения теплоты псрсгревсО/11с(роспользуемся выражением для интенсивности отвода теплоты перегрева и [5]:
пер
^С пер
</V„
Конкретизиру ем значение г. димеж линейной зависимостью а) вида:
г'пер %+а(\г01.
(21
1ер для чего зада-и(рис.3.
/ =0 - окончание
«0+01Л = О или «1='—.
Со
У^у Уш
Рис. 3. Схема для расчет киношки агверлевания отливки при горизонтальном лшье
Из площади треугольника ЛОВ находим, что
или
«о = лУ1лпер. Подставляя значения и а. в выражение для гпео> получаем
*пер =2Т|Гпер(1-$™
(3)
полагая при этом, что ооъем затвердевшего расплава Иго отношению к полному объему отливки V изменяется пропорционально изменению угловой координаты ф. характеризующей неравномерность отвода теплоты перегрев;! относительно верхней и нижней юны кристаллизатора (рис. 3. б)\
У„=У5Ш|.
Коэффициенты щ определим по координатам точек А 17 В м площади треу гольника ЛОВ, которая представляет собой полную теплоту перегрева О^ли координаты точки Л лл=0, /пср=а0 характеризуют момент начала затвердевания, то координаты точки В I' Г'ж затвердевания
Решая совместно (2) и (3) и интегрируя полученное выражение относителы^?1|ср пределах от ф до 71 (фронт кристаллизации симметричен относительно оси у) для /-го участка, находим полную теплоту перегрева, отводимую с фронта кристаллизации:
а. лер : ""I, VI*-прр(1-51П > V
= «лУГ*перО-5ш»£)л"ллллсри
После подстановки всех слагаемых в уравнение теплового баланса (I) и решения его относительно .и. ркончательно получим
М.=
' 2/?<7,ф
{/пер«' + Аэф(2^,-Л>1
+ ^кч (ЗЛ - ) - 9Л?<ЗЯ - )] I
30 ,
гтт^ггг^гтплглтгггп
(35), 2005 -
где
'пер = гпер Щ<Х-
• Ф\2
-БШ —) 2
19В
теплота пере-
грева, отводимая пропорционально удельному тепловому потоку по периметру кристаллизатора и объему затвердевшего рас-
плава;
т; =-
41
,ср
А/;
1Ср'
коэффициент
Ж пояс Я пояс /пояс
пропорциональности, учитывающий долю отведенной теплоты перегрева для /-й зоны
кристаллизатора; а — средний удельный
ср
тепловой поток в кристаллизаторе. Расчет кинетики затвердевания отливки по уравнению (4) необходимо проводить методом последовательных приближений.
Адекватность математической модели эксперименту проверяли путем сравнения расчетов и результатов исследования процесса формирования сплошной цилиндрической отливки из чугуна методом термического анализа и выливания жидкого остатка. Отливку диаметром 0,10 м получали в комбинированном кристаллизаторе с графитовым вкладышем толщиной 11,3-8,2-Ю-3 м и стальной втулкой 5,2Т0~3-8,3 Т0~3 м. Длина охлаждаемой части кристаллизатора — 0,20 м, неохлаждаемой — 0,14 м. Технологические параметры: скорость литья м>ср=0,3-60~ 1 м/с, температура чугуна в металлоприемнике Т= 1230-1270 °С, время формирования отливки в кристаллизаторе /=40 с.
Схема установки термопар и температурный режим 0 комбинированного кристаллизатора показаны на рис. 4, 5.
По мгновенному распределению температуры по толщине стальной втулки (рис. 6)
определяли удельные тепловые потоки, отводимые с поверхности отливки к охлаждающей среде по толщине стальной втулки в различных зонах кристаллизатора 0 (рис. 7, а), а также
Г----- +
/1 тЛ Ф ~ * 200 ' \ \
220
340
Рис. 4. Схема установки термопар: 1—15 — номера термопар
полный удельный тепловой поток в кристаллизаторе во времени (рис. 7, 5, аппроксимация по методу наименьших квадратов).
Анализ зависимостей показал, что наблюдается существенное различие в значениях теплового потока в верхней и нижней областях кристаллизатора. В условиях плотного контакта корки отливки и нижней поверхности графитового вкладыша в нижней области кристаллизатора 0 пояс термопар) тепловой поток имел максимальное значение #=2,34 МВт/м2. В то же время в верхней области тепловой поток значительнс ниже #=1,42 МВт/м2, что свидетельствует с наличии газового зазора на границе «корка-графитовый вкладыш» при усадке отливки.
Аналогичный характер изменения удельногс теплового потока наблюдался в зоне периодического контакта корки и кристаллизатора (II и II пояса термопар), но при более низких значения* теплового потока вследствие дальнейшей усадк* отливки и увеличения зазора. По мере движенш
Пшшс.
Жххожс
5 / 7 /
6 8
_ „Д, III, ...л, 1Л,,1 < 1 1 и-
- 9 * 11 /
¿С Л- -
10 12
- ■■■* 1 1 » ■ 1.1 1 1
25 ЯО 35 40 45 25 30 35 40 45 25 30 35 40 45
Рис. 5. Температурный режим стальной втулки комбинированного кристаллизатора: 1-12- номера термопар
ГГГТТгГ: ГГ ГГШС№РПГ.Р,
/31
- 3 (35). 2005 i
10 12 14 JS^l©"* .як
16
ж 10 XxlOr* ш
15
10 20 30 40
Рис. 6. Температурный режим стальной втулки комбинированного кристаллизатора
а б
Рис. 7. Значения удельного теплового потока в верхней и нижней зонах кристаллизатора: а — по толщине комбинированной стенки в различных поясах; б — в зависимости от времени фиксирования
отливки
отливки в нижней области кристаллизатора тепловой поток резко снижался до #=1,19 МВт/м2 и на выходе из кристаллизатора составлял #=0,90 МВт/м2. В верхней области расчетные значения удельного теплового потока составляли соответственно 0,81 и 0,61 МВт/м2.
Используя математическую модель, по значениям удельного теплового потока рассчитали зависимости толщины затвердевшей корки отливки от времени ее формирования для верхней и нижней областей кристаллизатора при ф=0—к. Вы-
числения проводили методом последовательных приближений с помощью алгоритма и программы расчета, написанной на языке программирования Pascal. Термофизические коэффициенты выбраны из работы [6]: с = 753,6 ДжДкг-град), с' = 837,4 ДжДкг-град), X = 35 ВтДм-град), г = 215 кДж/кг, Гсо= 1140-С, Тш= 1210-С, гпер= 92,1 кДж/кг, гЭф= 270,7 кДж/кг. Временные интервалы /=4: первый - 0-2 с, второй - 2-18, третий - 1838, четвертый - 38-40 с.
Результаты расчета кинетики затвердевания отливки в верхней и нижней областях комбинированного кристаллизатора показаны на рис. 8. Кривая 1 является расчетной зависимостью fy=f(t) толщины корки в верхней области кристаллизатора, кривая 2 — в нижней (полиномиальная аппроксимация). Области 3, показанные точками, представляют собой толщины корки в верхней и нижней областях кристаллизатора при различных температурах расплава в металлоприемнике, полученных экспериментально по методике выливания жидкого остатка [7]. Темплеты отливки после выливания приведены на рис. 9..
Анализ показал достаточно хорошее совпадение экспериментальных и расчетных значений толщин корки в верхней области кристаллизатора. Расчетная толщина корки на выходе из кристал-
лизатора здесь составила ^расч=9,010_3 м (кривая 1, рис. 8), эксперименты по выливанию показали значения £эксп=9,3-12,2-10~3 м (область 3, верх, температура расплава 7М230— 1270°С). Расхождение в значениях толщин корки незначительно, экспериментальная толщина корки несколько выше, что и должно быть, учитывая условия выливания. Расчетная динамика роста корки отливки во времени соответствует экспериментальной.
В нижней области, на выходе из кристаллизатора, расчетная толщина корки составила 2^асч=20,2-10~3 м (/=40 с), а по результатам экспериментов — ^эксп=28,0-10~3 м (/=36 с) (рис. 9, г). Здесь расхождение в результатах достаточно большое, порядка 10,0-10~3 м. Однако, учитывая значительную погрешность при определении толщины корки в нижней области кристаллизатора
30
о
-1—1-!—1—I
О 10
20 4с
30 40
Рис. 8. Кинетика затвердевания корки отливки в кристаллизаторе: 1 — верхняя область; 2 — нижняя область; 3 — выливание (верх, низ) при Т= 1230-1270 °С
CTTtn ГГ Г/ГГТШТГГПТ;
(35). 2005 -
f=0 с
=7*10"" м
К
=2мо" м
ДЯ=0,007 м
Рис. 9. Толщина корки в верхней и нижней ^
f=28 с =8*10-3 м =24* Ю-3 м AR=0,008 м
верх ,экел ' - низ -эксп
f=36 с верх -9*10"3 М
!-3
К
=28-10'" м AR=0,009 м
частях, а также величина смещения теплового центра горизонтальной отливки после вымывания жидкого остатка в зависимости от времени формирования X в кристаллизаторе
методом выливания, можно считать расчетные значения удовлетворительными, тем более, что определение толщины корки в нижней области кристаллизатора не имеет существенного значения с точки зрения влияния на стабильность процесса литья.
Таким образом, разработанная математическая модель подтверждает наличие несимметричности фронта кристаллизации по сечению отливки и позволяет качественно оценить динамику роста корки отливки в кристаллизаторе во времени.
Основной вывод, который можно сделать из полученных результатов — это удовлетворительное совпадение экспериментальных и расчетных значений толщин корки в верхней области кристаллизатора, что позволяет в дальнейшем использовать разработанный математический аппарат для расчетов кинетики затвердевания заготовок более сложных профилей и определения тепловых и технологических параметров при проектировании технологического процесса и оборудования горизонтального непрерывного литья.
Следует отметить соответствие полученных результатов результатам экспериментальных исследований других авторов [8].
Литература
1. Баранов O.A., Ветров Б.Г., Поль В.Б. и др. Непрерывное литье чугуна. М.: Металлургия, 1968.
2. Шатагин O.A., Сладкоштеев В.Т. и др. Горизонтальное непрерывное литье цветных металлов и сплавов. М.: Металлургия, 1974.
3. Тепловые режимы затвердевания стальных цилиндрических слитков при горизонтальном непрерывном литье / В.В. Соболев, П.М. Трефилов // Промышленная теплотехника. 1988. Т. 10. №2. С. 75-78.
4. Влияние термоконвекции на стабильность процесса горизонтального непрерывного литья / Е.Б. Демченко и др. //Литейное производство. 1987. №1. С. 20-21.
5. Демченко Е.Б. Столярова Г.И. Крутилин А.Н. Аналитический учет влияния конвекции на кинетику затвердевания непрерывнолитой заготовки // Металлургия. Мн.: Выш. шк., 1990. Вып. 24. С. 78-81.
6. Баландин Г.Ф. Основы теории формирования отливок. М.: Машиностроение, 1976.
7. Исследование процесса непрерывного горизонтального литья чугуна и внедрение в производство мероприятий по улучшению технологии литья и качества заготовок / В.А. Гринберг, И.В. Земсков, Е.Б. Демченко и др. Мн., 1979.
8. Исследование тепловых условий затвердевания чугунной заготовки в кристаллизаторе при непрерывном горизонтальном литье / В.И. Самсонов, B.C. Ан // Новые конструкционные материалы и методика обеспечения эксплуатационных характеристик деталей станков. М., 1984. С. 70—84.