CC BY
24
Л ТЕПЛОСНАБЖЕНИЕ, ВЕНТИЛЯЦИЯ, КОНДИЦИОНИРОВАНИЕ ВОЗДУХА, ГАЗОСНАБЖЕНИЕ И ОСВЕЩЕНИЕ
УДК 532.5:621.9.044 Князева И.А. - аспирант E-mail: [email protected] Вятский государственный университет
Адрес организации: 610000, Россия, г. Киров, ул. Московская, д. 36 Золотоносов Я.Д. - доктор технических наук, профессор E-mail: [email protected]
Багоутдинова А.Г. - кандидат технических наук, доцент E-mail: [email protected]
Казанский государственный архитектурно-строительный университет
Адрес организации: 420043, Россия, г. Казань, ул. Зеленая, д. 1
Выбор математической модели для описания теплового процесса лазерной сварки пружинно-витых каналов
Аннотация
Данная статья посвящена выбору математической модели описания теплового процесса лазерной сварки для производства пружинно-витых каналов теплообменных аппаратов, которые являются базовым элементом высокоэффективных кожухотрубчатых теплообменных аппаратов нового поколения. Предлагаемая математическая модель основывается на решении уравнения теплопроводности Фурье. Приведены выражения и обоснование их выбора для построения термических циклов сварки тонких пластин линейным быстродвижущимся источником тепла.
Ключевые слова: лазерная сварка, термический цикл, математическая модель, быстродвижущийся линейный источник.
Прямой пружинно-витой канал представляет собой пружину круглого сечения, витки которой жестко соединены между собой. Одним из возможных вариантов процесса получения таких каналов является технология плотной намотки проволоки на оправку с последующим объединением витков наматываемой проволоки между собой при помощи лазерной сварки.
Лазерная сварка обладает рядом неоспоримых преимуществ:
- высокая производительность сварочного процесса;
- малые зоны термического влияния;
- снижение времени теплового воздействия на деталь и, как следствие, уменьшение внутренних напряжений;
- возможность управления параметрами лазерного излучения в широком диапазоне значений.
Однако, несмотря на преимущества, лазерная сварка требует повышения стабильности, управляемости процесса. Сложность физических процессов, протекающих в расплаве сварочной ванны, образованной воздействием лазерного излучения, затрудняет опытное изучение данного процесса.
Стабильность процесса лазерной сварки, качество сварного шва напрямую зависят от правильности выбора режимов лазерной обработки. Подбор режимов обработки материалов лазерным излучением, в частности лазерной сварки, является трудоёмким процессом, так как требует проведения большого количества экспериментов.
Параметры режима обработки могут быть достаточно точно и сравнительно недорого получены путем математического моделирования технологического процесса. Моделирование физических процессов лазерной сварки позволяет ускорить решение задачи получения оптимальных режимов обработки за счет использования вычислительной техники. Также математическое моделирование сокращает затраты на выполнение опытных работ.
При воздействии сварочного источника в металле возбуждаются процессы тепло- и массопереноса, которые определяют изменение структуры металла, фазовые переходы и динамику напряженно-деформационного состояния [1].
Основными параметрами теплового процесса при сварке являются температура Т (х, у, г, ^ и плотность мощности концентрированного источника энергии д (х, у, х, ^
(электрической дуги, потока электронов, излучения лазера, низкотемпературной плазмы, газового пламени и т.д.). Параметры Т ид распределены в пространстве, заданном декартовыми координатами х, у, г, и во времени /. поэтому процесс теплопереноса является системой с распределенными параметрами.
В теории сварочных процессов разработаны методы расчета, основанные на решении прямой задачи теплопроводности. Вычислительный эксперимент над решением прямой задачи позволяет построить распределение температуры в пространстве и во времени при заданном источнике нагрева. Модель источника строится в соответствии с выбранной схемой, а выбор режима сварки ограничивается набором расчетных схем источника.
Математической моделью теплового процесса является уравнение теплопроводности с соответствующими начальными и граничными условиями. Начальные условия задаются в виде распределения температуры в заданный момент времени, чаще при t=0. В качестве граничных условий задаются либо температура на границах рассматриваемой области тела (граничные условия первого рода), либо тепловой поток (граничные условия второго рода).
В качестве системы координат наиболее часто используются линейная (декартова), цилиндрическая или сферическая системы координат. Выбор системы координат определяется формой тела, его размерами и характером действия концентрированного источника энергии. Для тел вращения используют цилиндрическую систему координат, для прочих тел - линейную.
Размерность математической модели определяется количеством пространственных переменных, вдоль которых рассматривается процесс теплопереноса и время. Выбор размерности модели определяется требуемой точностью, временем расчета, располагаемыми вычислительными ресурсами и математическим описанием тела.
В данном случае необходимо рассматривать тепловой процесс при стыковой сварке тонкостенных цилиндров быстродвижущимся источником, с мощностью, равномерно распределенной по толщине свариваемых деталей. Начало системы координат находится на пересечении оси луча лазера с верхней точкой стыка двух тонкостенных цилиндров. Луч лазер движется со скоростью сварки по направлению оси х.
В математической модели процессы теплопереноса при сварке описываются квазистационарным уравнением теплопроводности, а теплообмен с окружающей средой и металлом - различными нелинейными граничными условиями.
С целью упрощения решения поставленной задачи в математическом моделировании сварочных процессов применяют схемы нагрева тела: бесконечное или полубесконечное тело, бесконечная или полубесконечная пластина, плоский слой, бесконечный или полубесконечный стержень.
Источники тепла, встречающиеся в практических случаях сварки, также разнообразны. Их классифицируют по следующим признакам [1]:
- по признаку распределенности: сосредоточенные (точечный, плоский, объемный) и распределенные (по определенному закону ввода тепла в изделие) источники тепла;
- по времени действия: мгновенно и непрерывно действующие;
- по расположению относительно рассматриваемой точки во времени: неподвижные, подвижные, быстродвижущиеся источники тепла.
Для тела цилиндрической формы это уравнение теплопроводности необходимо преобразовать в пространстве цилиндрических координат г, в, г.
При сварке кольцевого соединения цилиндрических деталей подвижным сварочным источником, перемещающимся вдоль стыка кольцевого соединения со скоростью V (рис. 1), математической моделью теплового процесса является неоднородное дифференциальное уравнение в частных производных [2]:
(1)
где Т (г, в, г, ^ - температурное поле в цилиндрических координатах, а - коэффициент температуропроводности.
Система цилиндрических координат г, в, неподвижная. Она связана с источником, который перемещается со скоростью V вдоль стыка кольцевого соединения двух витков пружинно-витого канала.
Рис. 1. Схема кольцевого соединения витков пружинно-витого канала:
К, в, 2- система цилиндрических координат; q (г, I) - плотность мощности объемного источника
Для моделирования теплового процесса используем неоднородное уравнение теплопроводности с граничными условиями второго рода.
Так как сварка ведется лазерным излучением, то в качестве источника тепла для математического моделирования процесса необходимо выбрать быстродвижущийся источник тепла. Мощный быстродвижущийся источник является частным случаем непрерывно действующих источников. При перемещении источника с большой скоростью рассматривают случай предельного состояния, когда V—>°о. В этом случае предполагают, что энергия источника выделяется вдоль линии его перемещения, а тепло распространяется в направлении, перпендикулярном оси перемещения источника.
Необходимым условием лазерной сварки является то, что на удельный тепловой поток с/, пространство и время воздействия источника наложены ограничения, характерные для реальных высококонцентрированных источников энергии [2]:
0 <д (х, у, 0 <дтах х, у е (0, (Лу)
т е (о, о,
где с!„ - диаметр пятна нагрева теплового источника;
Т- время воздействия теплового источника.
Применение схемы быстродвижущегося источника определяется условием:
К-_?х=д-Х.о. Р)
дв дб д2~
При математическом моделировании процесса распространения теплоты при лазерной сварке витков пружинно-витого канала, несмотря на кривизну витков, он может быть приравнен к случаю нагрева тонких пластин линейным источником теплоты, так как любую окружность можно разделить на множество прямолинейных участков. Данное обстоятельство существенно упрощает уравнение теплопроводности. Но при этом необходимо учитывать, что при сварке спиральных швов происходит наложение тепловых потоков, так как на замыкающем участке кольцевого сварного шва уже
произошел нагрев проволоки от сварки предыдущего витка. Тогда схема лазерной сварки пружинно-витого канала будет иметь следующий вид:
Известно уравнение, описывающее изменение температуры в пластине при воздействии линейного источника теплоты мощностью с/ с равномерным распределением ее по толщине пластины, движущегося с постоянной скоростью [3]:
АТ =-----4й м, (3)
г?<5Л/4-77:ЛсрЬ
где АТ - приращение температуры, К;
с[ - эффективная мощность источника тепла, Вт;
V - скорость перемещения источника тепла, м/ч; ё - толщина пластины, м;
X - коэффициент теплопроводности, Вт/(м*К); ср - объемная теплоемкость, Дж/(м’К): t - время, с;
а - коэффициент температуропроводности, м2/с;
Ь - коэффициент температуроотдачи, с"1.
Наложение тепловых потоков при моделировании лазерной сварки по винтовой линии учитывается путем суммирования температурных полей. С учетом данного факта формула (3) приобретает следующий вид:
АТ = V-------- 4 рГ^к~Ып: (4)
“ 17(5^/477:ЛсрЬп п - число пересечений источником теплоты плоскости 1-1;
£?г - время, прошедшее с момента пересечения источником теплоты плоскости 1-1.
Исходя из принятой математической модели, на основе параметров, полученных при проведении опытных работ, в системе МаШСАБ были построены термические циклы в точках х=0,1, х=0,5, х=2,0, представленные на рис. 2.
Рис. 4. Термические циклы в точках х=0Д. х=0,5,х=2,0
Как видно из графика, на расстоянии 0,1 мм от оси шва максимальная температура составляет 1250 "С. На основании термических циклов можно судить о физических процессах, происходящих в сварочной ванне и структуре металла сварного шва.
Заключение
Описана математическая модель, применяемая для математического моделирования процесса лазерной сварки пружинно-витых каналов. Приведены зависимости и графики термических циклов, построенные на основе данных проведенных опытных работ. Установлено, что для построения более точной математической модели необходимо учитывать форму свариваемых деталей, что в значительной степени усложнит процесс построения модели, но позволит более качественно подобрать режимы лазерной сварки. Данный факт будет учтен в последующих работах по математическому моделированию процесса лазерной сварки, решению прямой и обратной задачи теплопроводности для определения режимов изготовления змеевиковых теплообменных элементов.
Список литературы
1. Петров Г.Л., Тумарев А.С. Теория сварочных процессов (с основами физической химии). - М.: Высшая школа, 1977. - 392 с.
2. Коновалов А.В., Куркин А.С., Неровный Э.Л., Якушин Б.Ф. Теория сварочных процессов: Учебник для вузов / Под ред. В.М. Неровного. - М.: Изд-во МГТУ им.
Н.Э. Баумана, 2007. - 752 с.
3. Негода Е.Н. Тепловые процессы при сварке: учеб. пособие; Дальневосточный государственный технический университет. - Владивосток: Изд- во ДВГТУ, 2008. - 125 с.
4. Язовских В.М., Беленький В.Я., Трутников Д.Н. Построение математических моделей сварки цилиндрических тел // Вестник ПГТУ. Сварка. - Пермь, 2012. - С. 5-17.
5. Язовских В.М. Математическое моделирование и инженерные методы расчета в сварке: Ч. 2. Тепловые процессы при сварке и моделирование в пакете МаШсас!. -Пермь: Изд-во ПГТУ, 2008. - 119 с.
TennocfiaS^kefiiie, Befimunaijiia, KotiQiiijiioHiipoBaHiie BcsQyxa,
sasocHaS^ceHue u ocBemetitie
Knyazeva I.A. - post-graduate student E-mail: [email protected] Vyatka State University
The organization address: 610000, Russia, Kirov, Moscovskaya St., 36 Zolotonosov Ya.D. - doctor of technical sciences, professor E-mail: [email protected]
Bagoutdinova A.G. - candidate of technical sciences, associate professor E-mail: [email protected]
Kazan State University of Architecture and Engineering
The organization address: 420043, Russia, Kazan, Zelenaya St., 1
Selection of mathematical model to describe the thermal process of the laser welding of spring-twisted channels Resume
Modeling of physical processes of laser welding speeds up the solution of optimum processing conditions with the help of computer technology. The same mathematical simulation reduces the cost of implementation of the experimental work.
Mathematical models of the thermal conductivity of the process with the appropriate boundary conditions. There are several methods for detecting the heat equation. They differ from each other in the accepted system of coordinates, and as a source.
This article deals with the choice of a mathematical model describing the heat of the laser welding for the production of spring-twisted-channel heat exchangers, which are the basic element of high-performance shell and tube heat exchangers of the new generation. Dependence of thermal cycles and graphs based on experimental studies conducted with data. Found that for the construction of more precise mathematical models must take into account the form of blanks, which greatly complicates the process of model building, but will allow higher quality modes to choose laser welding.
Keywords: laser welding, thermal cycle, mathematical model, a fast-moving line source.
References
1. Petrov G.L., Tumarev A.S. The theory of welding processes (with the fundamentals of physical chemistry). - M.: High. School, 1977. - 392 p.
2. Konovalov A.V., Kurkin A.S., Nerovny E.L., Jakushin B.F. The theory of welding processes: Textbook for Universities / ed. V.M. Nerovny. - M.: Publishing House of the MSTU. NE Bauman, 2007. - 752 p.
3. Nehoda E.N. Thermal processes in welding: studies, allowance; Far Eastern State Technical University. - Vladivostok: Publishing house FENTU, 2008. - 125 p.
4. Yazovskih V.M., Belenky V.Y., Trushnikov D.N. Construction of mathematical models of welding cylindrical bodies // Herald PGTU. Welding. - Perm, 2012.-P.5-17.
5. Yazovskih V.M. Mathematical modeling and computational methods in engineering welding: in 2 Part 2. Thermal processes in welding and modeling package Mathcad. -Perm: Perm State Technical University Publishing House, 2008. - 119 p.
