CC BY
59
УДК 621.791
ОПРЕДЕЛЕНИЕ ВЛИЯНИЯ ТЕПЛООТВОДА В МЕДНУЮ ПОДКЛАДКУ НА ПРОПЛАВЛЕНИЕ ТОНКОЛИСТОВЫХ МАТЕРИАЛОВ ПРИ АРГОННО-ДУГОВОЙ СВАРКЕ
© 2012 В.М. Бычков, Р.В. Никифоров, А.Н. Паутов, А.А. Логоня
Уфимский государственный авиационный технический университет
Поступила в редакцию 23.03.2012
Статья посвящена определению влияния теплоотвода в медную подкладку на размеры ширины шва при аргонно-дуговой сварке тонколистовых материалов и разработке методики построение численных моделей расчетов теплового поля с применением экспериментально полученных поправочных коэффициентов.
Ключевые слова: дуговая сварка, численное моделирование, теплопередача в подкладку, контактное термическое сопротивление
При освоении технологии сварки новых изделий приходиться оптимизировать геометрию проплавления экспериментальным путем, поэтому все большее применение находят количественные модели сварки для расчета глубины проплавления и ширины шва. Моделирование сварки встык тонколистовых материалов на медной подкладке сопряжено с некоторыми трудностями, поскольку имеет место значительный теплоотвод в подкладку, достаточно точно определить величину которого без экспериментов невозможно. Для получения расчетной модели с экспериментальными коэффициентами, описывающей реальное влияние теплоотвода на ширину стыкового шва при аргонно-дуговой сварке на медной подкладке проводилось сравнение расчетного поля температуры, полученного при численном моделировании, с экспериментально измеренными значениями температуры.
Для описания температурного поля T(x, у, z, ^ в плоском образце размером 50x100x2 мм использовали дифференциальное нелинейное уравнение теплопроводности:
с(Т )р(Т ) = Л(Т ) dt
d2T d2T d2T н--т- +
dx2
dy2
dz2
(1)
где р(Т) - плотность, кг/м3; А(Т) - коэффициент теплопроводности, Вт/(мхК); с(Т) - удельная теплоемкость, Дж/(кг°С). Зависимость те-плофизических свойств стали от температуры р(Т), А(Т), с(Т) принимали согласно литературным данным [4].
В использованной нами расчетной схеме сквозного проплавления тепловая мощность сварочной дуги для одной из свариваемых пластин частично вводится через ее боковую поверхность в плоскости стыка, а частично через верхнюю поверхность пластины. Изотермы температуры плавления от движущихся распределенного и сосредоточенного источников имитируют поверхность солидуса сварочной ванны Та (х, у, z)= Та. Линейный тепловой источник был задан как тепловой поток через плоскость, которая разделяет образец на 2 части (рис. 1) по формуле [1]:
xt ) = тгУехр
тг
( ^ )2 a
(2)
Бычков Владимир Михайлович, кандидат технических наук, доцент кафедры «Оборудование и технологии сварочного производства». E-mail: otsp@ugatu.ac.ru Никифоров Роман Валентинович, аспирант. E-mail: roman. nikiforov. svarka@mail.ru
Паутов Анатолий Николаевич, студент. E-mail: Njkzytrol@inbox.ru
Логоня Артем Алексеевич, студент. E-mail: artemugatu @list. ru
где а - размер конечного элемента, м; И - толщина стали, м; ДХ=Х0-^ТГМБ - перемещение по оси Х.
Нормально-круговой источник нагрева на верхней поверхности рассчитывался по формуле[1]:
q(х, у, г) = - q2 ехр(--(АУ2 +г2)), (3) л
где к - коэффициент сосредоточенности дуги, м-2.
Граничные условия теплопередачи с внешних поверхностей определяли удельный
тепловой поток при конвекции и излучении. Радиационное излучение учитывалось только на поверхности А и рассчитывалось по формуле:
Ягт
где T - температура тела; Тс - температура окружающей среды (20°С); Co - постоянная Сте-фана-Больцмана C0=5,67х10-8 Вт/(м2-К4); в -излучательная способность стали, которую принимали согласно данным работы [2].
Удельный тепловой поток при конвективном теплообмене принимали равным q2k=25 Вт/(м2-К) по данным работы [2] для свободных поверхностей кроме А, B и D, как показано на рис. 1.
2 Л
Кмедъ + Кс
= 1,02 • 104
Вт
м
град'
(5)
где Хвоз - теплопроводность воздуха, Вт/(м^°С); Ьмедъ - микронеровность поверхности медной подкладки, мкм; ксталъ - микронеровность поверхности стали, мкм.
Величина микронеровностей поверхности свариваемого образца из листовой стали
12Х18Н10Т в состоянии поставки и медной подкладки определяли с помощью портативного измерителя шероховатости ТЯ 100. Величина микронеровностей Яа на поверхности стали И31ее1=0,54-0,72 мкм (для расчетов приняли 0,66 мкм), а на поверхности меди ИСи=2,1-3,61 мкм (приняли 2,90 мкм). Зависимость Х(Т) воздуха принимали по данным работы [3]. Для вычисления коэффициента, учитывающего площадь фактического контакта стального образца и медной подкладки, воспользовались выражением, приведенным в справочнике [3]:
8ф=Р/И
(6)
Рис. 1. Расчетная схема модели
Граничные условия контакта сталь-медь рассматривались как сумма термических про-водимостей воздушной прослойки авоз между микровыступами контактируемых материалов и металлического контакта аМ, для расчета которых использовали методику из справочника [3]. Термическая проводимость воздушной прослойки описывается выражением:
где Р - давление на образец, кг/см2; Н - твердость материала по Бринеллю, кг/см2 .
Для расчета этого коэффициента использовали экспериментально измеренное давление прижимов Р=9,85 кг/см2 и данные по твердости меди М1 из марочника [4]. В нашем случае получили для металлического контакта 8ф=0,12, и поскольку эта величина относительно мала, то в соответствие с методикой [3], можно заключить, что термическая проводимость контакта приближенно может быть рассчитана только по значению авоз. Поскольку теплопроводность воздуха зависит от температуры, для расчета была использована линейная зависимость, полученная аппроксимацией табличных данных в программе ЗТАТКТГСА в виде выражения:
а(Т)=10200+32Т
(7)
Эффективный КПД сварочной дуги определяли методом калориметрирования в стоячей воде по формуле:
Лэф
сшАТ VII.....
(8)
где с - теплоемкость воды, Дж/(кг°С); т -масса воды, кг; ДТ - перепад температур при калориметрировании, °С; и - напряжение на дуге, В; I - сила тока, А; 1:св - время сварки, с. Режимы сварки образцов и полученные при этом значения эффективного КПД приведены в таблице 1.
Таблица 1. Данные по режимам сварки и коэффициенты модели
Режим Сила тока, А Скорость сварки, м/ч Напряжение на дуге, В КПД, % Коэффициент сосредоточенности 1 -2 дуги к, см
(а) 115 12 10,10 61 5,5
(б) 115 15 10,27 64 8
(в) 115 18 10,31 65 13
Сварка производилась на образцах из нержавеющей стали 12Х18Н10Т размерами 100х50х2 мм на стенде, приведенном на рис. 2. При сварке на всех режимах использовали вольфрамовые электроды WL-20 диаметром 3 мм с углом заточки 30°, длина дуги 3 мм.
Рис. 2.Общий вид экспериментальной установки
Температурное поле определялось численным решением задачи нелинейной теплопроводности, пример расчета которой приведен на рис. 3.
Рис. 3. Пример расчета в пакете ANSYS/ Multiphysics.
Мощность источников нагрева и коэффициент сосредоточенности дуги подбирались по следующему алгоритму - вначале по полуширине обратного валика подбиралась эффективная мощность линейного источника нагрева, после чего - по полуширине сварного шва выбирался коэффициент сосредоточенности дуги для нормально-кругового источника нагрева. В результате моделирования было установлено, что мощность линейного источника нагрева должна составлять 1/5 от общей мощности для того, чтобы ширина шва с лицевой и обратной стороны соответствовала экспериментальным данным.
Проверка расчетной методики определения проводимости теплового контакта осуществлялась путем записи термических циклов с термопар в медной подкладке, установленных в соответствии с рис. 4 и последующим их сравнением с термическими циклами узлов конечно-элементной сетки медной подкладки при расчете контактной тепловой задачи. При этом термопары типа «К» зачеканивали в медную пластину на глубину 1,5±0,5 мм и заливали термостойким герметиком.
Мед ная подкладка
1
10
77 Т2 ТЗ П
Рис. 4. Схема установки термопар в медную подкладку
Термические циклы снимались с помощью цифрового 4-х канального измерителя температуры СЕМТЕЯ 309. На рис. 5 представлены термические циклы режимов (а)-(в).
Время сварки, сек а)
125 НО
<75
а
^ 80
та
О-
ё 65 2
35 20
0 4 8 12 16 20 24 28 32 36 Время сварки, сек
б)
в)
Рис. 5. Термические циклы режимов а)-в) в соответствии с таблицей 1.
В результате сопоставления экспериментальных термических циклов с расчетными, было установлено, что термическая проводимость, рассчитанная по теплопроводности воздушной прослойки между микронеровностями шероховатостей поверхностей, не отражает реальных условий сварки, по-видимому, в связи с отсутствием равномерного прилегания свариваемых деталей к медной подкладке. Поэтому в расчетной модели для определения значения авоз необходимо использовать
экспериментальный коэффициент, который в наших условиях составил 0,8.
Выводы:
1. Для получения расчетной модели, адекватно описывающей влияние теплоотвода на ширину стыкового шва при аргонно-дуговой сварке на медной подкладке необходимо введение экспериментально полученных коэффициентов контактной теплопроводности между свариваемыми листами и подкладкой.
2. Для аргонно-дуговой сварки стальных образцов толщиной 2 мм поправочный коэффициент к расчетному значению термической проводимости составил 0,8.
СПИСОК ЛИТЕРАТУРЫ:
1. Рыкалин, Н.Н. Расчеты тепловых процессов при сварке. - М: Машиностроительная литература, 1951. 291 с.
2. Шейндлин, А.Е. Излучательные свойства твердых материалов. - М.: Энергия, 1974. 472 с.
3. Шлыков, Ю.П. Контактное термическое сопротивление / Ю.П. Шлыков, Е.А. Ганин, С.Н. Ца-ревский. - М.: Энергия, 1977. 328 с.
4. Сорокин, В.Г. Стали и сплавы. Марочник / В.Г. Сорокин, М.А. Гервасьева. - М.: Интермет Инжиниринг, 2001. 608 с.
DEFINITION THE INFLUENCE OF HEAT TRANSFER IN COPPER PLATE ON THIN SHEET MATERIALS MELTING
DURING TIG-WELDING
© 2012 V.M. Bychkov, R.V. Nikiforov, A.N. Pautov, A.A. Logonya
Ufa State Aviation Technical University
Article is devoted to definition the influence of heat transfer in copper plate on sizes of weld width at TIG-welding of thin sheet materials and working out the technique of creating the numerical models of calculations of a thermal field with application the experimentally received correction factors.
Key words: TIG-welding, numerical modeling, heat transfer in a plate, contact thermal resistance
Vladimir Bychkov, Candidate of Technical Sciences, Associate
Professor at the Department "Equipment and Technologies of
Welding Enterprises". E-mail: otsp@ugatu.ac.ru
Roman Nikiforov, Post-graduate Student. E-mail:
roman. nikiforov. svarka@mail.ru
Anatoliy Pautov, Student. E-mail: Njkzytrol@inbox.ru
Artem Logonya, Student. E-mail: artem_ugatu@list.ru
