CC BY
13
УДК 621.74:669.714
В. Н. Лаптинский, А. А. Романенко, В. Ю. Стеценко, К. Н. Баранов
МЕТОД РАСЧЕТА ПАРАМЕТРОВ ЗАТВЕРДЕВАНИЯ ОТЛИВКИ ПРИ ЛИТЬЕ НАПРАВЛЕННЫМ ЗАТВЕРДЕВАНИЕМ ОТ ВОДООХЛАЖДАЕМОГО СТЕРЖНЯ
UDC 621.74:669.714
V. N. Laptinskii, A. A. Romanenko, V. Yu. Stetsenko, K. N. Baranov
THE METHOD FOR CALCULATING SOLIDIFICATION PARAMETERS OF A CASTING DURING CASTING BY DIRECTIONAL SOLIDIFICATION FROM WATER-COOLED ROD
Аннотация
Созданы опытная установка для литья полых заготовок из силумина методом направленного затвердевания от водоохлаждаемого стержня и методика определения температурного поля в полученных отливках. Разработан полуэмпирический метод численно-аналитического расчета затвердевания отливок по показаниям термопар, установленных в них, позволяющий учитывать геометрию системы струйного охлаждения.
Ключевые слова:
направленное затвердевание, отливка, стержень, силумин.
Abstract
A pilot plant has been created to cast hollow billets of silumin using the method of directional solidification from a water-cooled rod, as well as a technique for determining the temperature field in the castings obtained. A semi-empirical method has been developed for the numerical-analytical calculation of solidification of castings according to the readings of thermocouples installed in them, which makes it possible to take into account the geometry of the jet cooling system.
Keywords:
directional solidification, casting, rod, silumin.
Введение
Расчет затвердевания отливки с помощью численного моделирования на ПК в последние десятилетия стал основным методом в расчетах затвердевания отливок. Тем не менее этот метод, наряду с определенными преимуществами, имеет и ряд недостатков, главными из которых являются большая трудоемкость и отсутствие наглядности. Поэтому, несмотря на расширение применения численных методов и ПК в расчетах затвердевания отливки, полуэмпирические численно-аналитиче-
ские методы не утратили своей актуальности. Эти методы значительно проще математического моделирования на ПК и дают более наглядное представление о влиянии различных факторов на исследуемый процесс. Весьма ценными являются численно-аналитические методы получения соответствующих температурных распределений, которые содержат структурные параметры, достаточно эффективно определяемые из эксперимента. Однако такие методы применительно к расчету затвердевания отливки при интенсивном охлаждении к настоящему времени пока еще
© Лаптинский В. Н., Романенко А. А., Стеценко В. Ю., Баранов К. Н., 2019
не разработаны [1, 2].
Известно, что главная задача теории затвердевания отливки состоит в выявлении связи между комплексом параметров условий литья и свойствами отливки из различных литейных сплавов. Решение этой задачи конструктивными математическими методами позволяет получить соответствующие инженерные формулы, что даст возможность эффективно управлять процессом затвердевания отливки. В результате экспериментальных исследований установлено, что свойства реальных отливок определяются динамикой затвердевания расплава при определенных условиях технологии литья.
Формирование слитка при литье -сложный процесс, сопровождающийся теплофизическими, физико-химическими, гидродинамическими, механическими и другими явлениями. Среди них наиболее важными являются теплофи-зические процессы. В настоящее время в научной литературе отсутствуют эффективные методы расчета затвердевания отливки, позволяющие получить
соответствующие температурные распределения, которые необходимо иметь для расчета затвердевания отливки.
В данной работе описаны опытная установка для литья силуминовой отливки методом направленного затвердевания от водоохлаждаемого стержня и методика для определения температурного поля в отливке [3].
Цель исследования. Разработка полуэмпирического метода численно-аналитического расчета затвердевания силуминовой отливки, полученной методом направленного затвердевания от водоохлаждаемого стержня.
Методика проведения эксперимента
Принципиальная схема процесса литья тонкостенных полых заготовок из силумина методом направленного затвердевания от водоохлаждаемого стержня со струйной системой охлаждения и стальной гильзой представлена на рис. 1 [4].
Рис. 1. Принципиальная схема процесса литья направленным затвердеванием от водоохлаждаемого стержня
Для осуществления процесса литья заготовок из силумина спроектирован соответствующий водоохлаждаемый стержень. В качестве материала этого стержня использовали стальную гильзу с толщиной стенки 4 мм. Наружная поверхность гильзы выполнена с конусностью 0,6° к его продольной оси и имеет длину 2...2,5 диаметра гильзы. Принцип функционирования водоохлаждаемого
стержня представлен на рис. 2. Охлаждение стержня происходит следующим образом. При включении насоса охлаждающая жидкость из подводящего патрубка 3 продавливается через экран с отверстиями 2, после чего в виде затопленных струй равномерно и интенсивно охлаждает внутреннюю поверхность гильзы 1 и отводится через патрубок 4.
Рис. 2. Схема функционирования водоохлаждаемого стержня со струйной системой охлаждения:
1 - гильза; 2 - отверстия; 3 - подводящий патрубок; 4 - отводящий патрубок
Затопленно-струйный метод
охлаждения обеспечивает наиболее интенсивный теплоотвод от внутренней стенки кристаллизующейся отливки [5].
Распределение температуры по времени в отливке из сплава АК12 определяли экспериментально на опытной установке, общий вид которой представлен на рис. 3 [3]. Для получения опытных образцов использовали стальную гильзу наружным диаметром 65 мм, высотой 120 мм и толщиной стенки 4 мм.
В качестве шихты использовали
чушку из сплава АК12, ее химический состав: А1 - основа; - 10,5.11,2 %; Си - 0,51 %; Мп - 0,3 %; Бе - 0,5 %. Плавку вели в термической печи сопротивления марки «8по1-1300» в шамото-графитовом тигле. Разливку жидкого металла проводили при температуре 850 °С в предварительно разогретый до 500 оС стальной съемный захват. Модификаторы не применялись. Расход воды в водоохлаждаемом стержне был равен 12 м3/ч. Для регистрирования температуры применяли термопреобразователи, изготовленные из
термостойкого термоэлектродного провода ПТН/ХА/1/2х0,2(0,3)/1 диаметром 1 мм. В средней части рабочей полости съемного захвата располагали две термопары. Схема расположения термо-
пар в отливке представлена на рис. 4 [3]. Показания термопар фиксировали при помощи блока регистрации температуры фирмы ЬО и ноутбука.
Рис. 3. Общий вид экспериментальной установки с водоохлаждаемым стержнем диаметром 65 мм
Рис. 4. Схема экспериментальной установки с термопарами: 1 - съемный захват; 2 - графитовое кольцо; 3, 4 - фланец; 5 - гильза; 6, 7 - кольцо; 8 - подводящий патрубок; 9 - отводящий патрубок; 10 - заглушка; 11 - экран; 12 - болт; 13 - отливка
Из сплава АК12 были получены опытные образцы диаметром 105 мм, высотой 100 мм и толщиной стен-
ки 20 мм. Графики экспериментальных показаний термопар во время затвердевания отливки представлены на рис. 5 [3].
О
800
600
га
ср >
га 400 ^
ф с
ф
200
10 15 20
Время затвердевания, с
25
30
Рис. 5. Графики показаний термопар во время затвердевания отливки из сплава АК12: 1 - показания
термопары Гь 2 - показания термопары г2
Методика расчета температурного
поля в отливке и получение зависимостей, описывающих движение фронта затвердевания в ней
Расчет затвердевания отливки связан с решением задачи, которую можно отнести к обратным задачам теплопроводности. При этом методика расчета температурных полей является основой для изучения процессов затвердевания отливки. В данной работе ее математическую основу составляют методы интерполяции и метод наименьших квадратов. Методы интерполяции обеспечивают хорошую точность аппроксимации функций на небольших отрезках; метод наименьших квадратов позволяет строить приближенные формулы, пригодные на достаточно больших отрезках.
Методика расчета температурного поля в отливке водоохлаждаемого стержня является полуэмпирической, при этом используются показания тер-
мопар, количество которых определяется из конструктивных особенностей и технологических условий литья. Схема расположения термопар в установке приведена на рис. 6.
Термопара Г3 технологическая,
и она используется для слежения за уровнем температуры, постоянно поддерживаемой в графитовой форме на уровне Т (г) = 550 оС .
Основным структурным элементом этой методики является определение температурного поля Т (г, г) в точках отливки, в которых установлены
термопары Г1
(< г < г < ^2): тей [3, 6, 7]:
и г2 (см. рис. 6) в виде зависимос-
Г2 - Г
Т (Г, г) - р (Г, г ) =
1Т2 (г)+^ т1 (г);
Г1 - Г2
(1)
0
5
1п *
1п ^
Т (г г) = Т2 (г) + -Т- Т1 (Г) ;
1п г 1п г
г <Г <^ <г<
(2)
где Т (г), Т2 (г) - зависимости температуры от времени г в точках расположе-
ния термопар г1 , г2. При этом начальный момент времени г0 = 4 с соответствует началу затвердевания отливки (см. рис. 5), ^ - окончанию затвердевания отливки, т. е. условию Т (Я2, ) = Тк
(здесь Тк - температура кристаллизации силумина).
Рис. 6. Схема расположения термопар: 1 - гильза; 2 - отливка; 3 - графитовая форма; Я0 - внутренний радиус гильзы; Я1 - внешний радиус гильзы (внутренний радиус отливки); Я2 - внешний радиус отливки (внутренний радиус графитовой формы); Я3 - внешний радиус графитовой формы; гь г2 - расстояния от оси до термопар в отливке; г3 - расстояние от оси до технологической термопары в графитовой форме
Значения Т (г, г0) при Я < г < Я2
получаем на основе (1), полагая Т (г, г0)- р (г, г0), т. е. экстраполяцией
на промежуток [Я1, Я2 ]. Аналогично поступаем с формулой (2). В [3] формула (2) получена на основе подхода к интерполированию функций, предложенного в [8].
Очевидно, для описания Т(г, г) в отливке следует использовать именно эти формулы, поскольку для радиальной аппроксимации удобен метод двухточечной аппроксимации из-за небольшого количества термопар в отливке. Метод наименьших квадратов для аппроксимации по времени является более эффективным (количество экспериментальных точек теоретически не ограничено). Правая часть в (2) соответствует решению стационарной задачи об определении распределения температуры в круговой цилиндрической стенке (обо-
лочке) (см., например [9, с. 38]). В силу известного свойства логарифмической функции приближения (1), (2) обладают практически одной точностью для достаточно тонких цилиндрических стенок. К таким стенкам относится стенка данной установки. Очевидно, для толстых стенок формула (1) может значительно уступать формуле (2) по точности. Таким образом, при расчете Т (г, г)
можно учитывать геометрию системы охлаждения.
Для расчета температурного поля в отливке воспользуемся изложенной методикой и формулами (1), (2). Эти аппроксимации получены на основе результатов измерений температуры в отливке с помощью двух термопар.
Простейшие аппроксимации по времени, полученные по методу наименьших квадратов, имеют вид [3, 10]:
Т (г )=Т ( г , г) - р2 ( Г , г )=С() + ^+сС^г2; (3)
Т (г) = Т (г, г) - Р3 (г, г) =
= Ьг) + Ь(г ]г + Ьг V2 + Ьг )г3; (4) т (г ) = т (Г', г) - Р4 (Г', г ) =
с ) + с{' )г + С 2г )г2 + С 3(г )г3 + с 4г )г4,
г = 1,2.
(5)
Числовые значения коэффициентов для (3)-(5) приведены в табл. 1-3, в которых через г обозначены номера термопар в отливке (см. рис. 6).
Табл. 1. Расчетные значения коэффициентов для Р2 (г, г)
г °о а1 а2
1 0,67040735-103 -0,12007794-102 0,28013429
2 0,61236635-103 -0,78197418-101 0,20148279
Табл. 2. Расчетные значения коэффициентов для Р3 (г, г)
г ¿0 ¿2 ь3
1 0,80278625 -103 -0,53942418-102 0,41564778-101 -0,10825945
2 0,67778615-103 -0,28543247•102 0,21171176 •Ю1 -0,53500306 •Ю-1
Табл. 3. Расчетные значения коэффициентов для Р4 (г, г)
г С0 С1 С2 С3 С4
1 0,92772364 •Ю3 -0,10812164 • 103 0,12100994 •Ю2 -0,58358781 0,995202 •Ю-2
2 0,73985664 -103 -0,55460174 -102 0,60640542 -101 -0,28964949 0,494429 •Ю-2
Качественный анализ полученных аналитических аппроксимаций различных порядков может быть выполнен на основе изучения соответствующих графиков. На рис. 7 представлены графики экспериментальных и аппроксимацион-ных зависимостей температуры от времени затвердевания отливки (для термопар Г , Г2), описываемых полиномами 4-го порядка.
Сравнительный анализ полученных зависимостей показывает, что применение формул (1)-(5) представляется вполне оправданным для приближенного описания температурного поля в отливке.
На основе (3)-(5) по форму-
лам (1), (2) получим явные приближенные распределения температурного поля Т (Я2, г) на границе отливка - графитовая форма:
Т (Я2, г)-^Т2 {г) + ^^ {г); (6)
1п 1п -^
т(^= Т2(г)+ Т(г), (7) 1п Г 1п Г
т. е. экстраполяцией на промежуток [ Я2 ] (см. рис. 6).
Т, °С 660
640
620
600
580
560
540
520
6 8 10 12 14 16 1
20 с
Рис. 7. Графики экспериментальных и аппроксимационных по 75 точкам зависимостей температуры от времени для термопар г1 и г2 (точками обозначены экспериментальные данные, сплошными линиями - полиномы четвертого порядка)
Полагая в (6), (7) Т (Я2, г) = Тк (где Тк - температура кристаллизации силумина) и используя (3)-(5), для определения ^ получим алгебраические уравнения степеней соответственно 2, 3, 4. В частности, для полинома 2-й степени на основе приближений (6) получим простую формулу [3, 11]
=
2(- Тк)
- 442 К - Тк) -
(8)
где ^ - коэффициенты, которые находятся приведением подобных членов при одинаковых степенях г в (6) при
Т(Я2, г) = Тк.
Аналогичная формула может быть получена на основе приближения (7).
На рис. 8 представлены графики аппроксимационных зависимостей температуры от времени г, построенные по формулам (6), (7). Значения температуры для термопар г1 и г2 вычислены на основе аппроксимационных зависимо-
стей от времени, полученных на основе полиномов 4-го порядка.
Формулы (1), (2) на основе (3)-(5) использованы для получения приближенных соотношений для закона движения фронта затвердевания г (г) в отливке и расчета времени ее затвердевания г у-.
Полагая в (1), (2) Т (г (г) = Тк,
в результате несложных выкладок получим приближенные формулы для закона движения фронта затвердевания [3, 12] (приближенные решения задачи стефановского типа [13]):
•(0 <
Тк - Т (г) г ^ Тк - Т (г)
Т2(г)-Т(г) 2 Т(г)-Т2(г)
(9)
Тк-Т, (г)
•(О - г/1
Тк-Т (г)
Т1(г )-Т7 (г ) Т (г )-Т1(г
г0 < г < ^ .(10)
На рис. 9 представлены графики, построенные на основе формул (9), (10) с использованием полиномов 4-го порядка для Т (г) и Т2 (г).
4 6 8 10 12 14 16 18 С
Рис. 8. Графики зависимостей Т (Я2,г) (сплошная линия - зависимость на основе формулы (6), точками обозначена зависимость на основе формулы (7))
Рис. 9. Графики зависимостей г (г) (сплошная линия - зависимость на основе формулы (9), точками обозначена зависимость на основе формулы (10))
В результате расчетов на основе (8)-(10) установлено, что при температуре кристаллизации Тк = (578 ± 2) °С
фронт кристаллизации в отливке достигает ее границы Я2 = 52,5 мм при г^- —16 с .
Согласно эксперименту, г^ — 18 с.
Таким образом, формулы (1)-(5) вполне приемлемы для приближенного описания температурного поля в отливке, закона движения фронта затвердевания и расчета времени ее затвердевания.
Замечание. Погрешность экспериментальной технологической установки определяется:
1) погрешностью показаний термопар (3 % - паспортные данные);
2) суммарной погрешностью, определяемой неточностью размещения термопар, неточностью размеров графитовой втулки, гильзы и отливки заданных толщин (не более 7 %);
3) неточностью соблюдения времени заливки (не более 7 %).
Аппроксимация температурных полей является удачной, поскольку отклонение экспериментального времени затвердевания отливки от расчетного составляет около 15 %. Это допустимо условиями эксплуатации технологических установок для литья металлов и сплавов методом направленного затвердевания. Ввиду трудностей изготовления экспериментальной установки необходимое количество экспериментов для полной статистической обработки результатов не проводилось. Погрешность, приведенная в п. 1, определяется паспортными данными термопар; погрешности, приведенные в пп. 2 и 3, вычисляют раздельно, т. е. не прибегая к созданию всей установки. Поскольку толщина отливки достаточно мала, то погрешностью аппроксимации температурного поля в ней можно пренебречь.
Из анализа результатов эксперимента и расчетных данных видно, что эмпирические формулы для описания закона движения фронта кристаллизации, полученные на основе аппроксимации температурного поля в отливке, могут быть использованы в соответствующих расчетах при литье в производственных установках направленного затвердевания (методом намораживания).
Заключение
На основании результатов выполненных исследований можно сформулировать следующие выводы:
- разработана экспериментальная установка литья полых заготовок из силумина методом направленного затвердевания от водоохлаждаемого стержня;
- выполнены экспериментальные измерения температурного поля в отливке;
- разработана полуэмпирическая методика численно-аналитического расчета затвердевания отливки из силумина АК12 при литье направленным затвердеванием от водоохлаждае-мого стержня;
- получены простые приближенные формулы для закона движения фронта и времени затвердевания отливки.
Перечисленные результаты в целом представляют собой простой полуэмпирический метод расчета параметров затвердевания отливок в устройствах литья методом направленного затвердевания от водоохлаждаемого стержня.
Исследования авторов по данной тематике частично поддержаны Белорусским республиканским фондом фундаментальных исследований, грант № Т13-099.
СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ
1. Беленький, А. А. Математическое моделирование и оптимизация процессов литья и прокатки цветных металлов / А. А. Беленький. - Москва: Металлургия, 1983. - 160 с.
2. Кац, А. М. Теплофизические основы непрерывного литья слитков цветных металлов и сплавов / А. М. Кац, Е. Г. Шадек. - Москва: Металлургия, 1983. - 208 с.
3. Разработка численно-аналитического метода расчета скорости затвердевания отливки, полученной методом намораживания на водоохлаждаемом стержне: отчет о НИР (закл.) / ИТМ НАН Беларуси; науч. рук. В. Н. Лаптинский. - Могилев, 2015. - 120 с. - № ГР 20131537.
4. Способ получения полой отливки: пат. БУ 18286 / В. Ю. Стеценко, Е. И. Марукович. - Опубл. 30.06.2014.
5. Марукович, Е. И. Модифицирование сплавов / Е. И. Марукович, В. Ю. Стеценко. - Минск: Беларуская навука, 2009. - 192 с.
6. Баранов, К. Н. Расчет температурного поля в отливке, полученной методом направленного затвердевания / К. Н. Баранов, В. Н. Лаптинский, А. А. Романенко // Материалы, оборудование и ресурсосберегающие технологии: материалы Междунар. науч.-техн. конф. - Могилев: Белорус.-Рос. ун-т, 2015. - С. 21-22.
7. Баранов, К. Н. К методике расчета температурного поля в отливке, полученной методом направленного затвердевания / К. Н. Баранов, В. Н. Лаптинский, А. А. Романенко // Материалы, оборудование и ресурсосберегающие технологии: материалы Междунар. науч.-техн. конф. - Могилев: Белорус.-Рос. ун-т, 2017. - С. 16.
8. Латинский, В. Н. К методике решения линейных интерполяционных задач / В. Н. Лаптинский // Еругинские чтения - XI: тез. докл. Междунар. мат. конф., Гомель, 24-26 мая 2006 г. - Минск: Ин-т мат. НАН Беларуси, 2006. - С. 161.
9. Основы теплопередачи в авиационной и ракетно-космической технике / В. С. Авдуевский [и др.]. -Москва: Машиностроение, 1975. - 624 с.
10. Лаптинский, В. Н. Аналитический расчет времени затвердевания отливки, полученной методом направленного затвердевания / В. Н. Лаптинский, А. А. Романенко // Материалы, оборудование и ресурсосберегающие технологии: материалы Междунар. науч.-техн. конф. - Могилев: Белорус.-Рос. ун-т, 2018. - С. 38.
11. Лаптинский, В. Н. Формула для расчета времени затвердевания отливки, полученной методом направленного затвердевания / В. Н. Лаптинский, А. А. Романенко // Материалы, оборудование и ресурсосберегающие технологии: материалы Междунар. науч.-техн. конф. - Могилев: Белорус.-Рос. ун-т, 2019. - С. 53.
12. Баранов, К. Н. К расчету времени затвердевания отливки, полученной методом направленного затвердевания / К. Н. Баранов, В. Н. Лаптинский, А. А. Романенко // Материалы, оборудование и ресурсосберегающие технологии: материалы Междунар. науч.-техн. конф. - Могилев: Белорус.-Рос. ун-т, 2016. - С. 26.
13. Карташов, Э. М. Аналитические методы в теории теплопроводности твердых тел / Э. М. Кар-ташов. - Москва: Высшая школа, 1985. - 480 с.
Статья сдана в редакцию 20 июня 2019 года
Валерий Николаевич Лаптинский, д-р физ.-мат. наук, проф., ИТМ НАН Беларуси. E-mail: [email protected].
Алексей Андреевич Романенко, канд. физ.-мат. наук, доц., Белорусско-Российский университет. E-mail: [email protected].
Владимир Юзефович Стеценко, канд. техн. наук, ИТМ НАН Беларуси. E-mail: [email protected]. Константин Николаевич Баранов, науч. сотр., ИТМ НАН Беларуси. E-mail: [email protected].
Valery Nikolaevich Laptinskii, DSc (Phys.-Mat.), Prof., Institute of Metal Technology of the National Academy of Sciences of Belarus. E-mail: [email protected].
Alexey Andreyevich Romanenko, PhD. (Phys.-Mat.), Associate Prof., Belarusian-Russian University. E-mail: [email protected].
Vladimir Yuzefovich Stetsenko, PhD (Engineering), Institute of Metal Technology of the National Academy of Sciences of Belarus. E-mail: [email protected].
Konstantin Nikolayevich Baranov, research fellow, Institute of Metal Technology, National Academy of Sciences of Belarus. E-mail: [email protected].
