Simulation of the mode of electron-beam welding of a thin-wall structure from AD31 alloy

Sergei O. Kurashkin; Курашкин Сергей Олегович; Yuri N. Seregin; Серегин Юрий Николаевич; Vadim S. Tynchenko; Тынченко Вадим Сергеевич; Aleksandr V. Murygin; Мурыгин Александр Владимирович; Svetlana V. Kotelnikova; Котельникова Светлана Владимировна

doi:10.31772/2712-8970-2022-23-2-296-304

Моделирование режима электронно-лучевой сварки тонкостенной кон-струкции из сплава АД31

Авторы: Курашкин С.О.¹, Серегин Ю.Н.¹, Тынченко В.С.¹, Мурыгин А.В.¹, Котельникова С.В.¹
Учреждения:
1. Сибирский государственный университет науки и технологий имени академика М. Ф. Решетнева
Выпуск: Том 23, № 2 (2022)
Страницы: 296-304
Раздел: Раздел 3. Технологические процессы и материалы
URL: https://journals.eco-vector.com/2712-8970/article/view/546134
DOI: https://doi.org/10.31772/2712-8970-2022-23-2-296-304
ID: 546134

Цитировать

Полный текст

Аннотация
Полный текст
Об авторах
Список литературы
Дополнительные файлы
Статистика

Аннотация

В статье содержатся результаты, полученные авторами при исследовании возможности применения электронно-лучевой сварки (ЭЛС) для тонкостенных конструкций из алюминиевого сплава АД31. Сегодня ЭЛС подобных конструкций не используется вследствие отсутствия технологии. В настоящее время на производстве используются другие технологии соединения подобных конструкций, однако они имеют высокую себестоимость, причина которой объясняется высоким процентом брака. Предложенный авторами способ применения ЭЛС позволит существенно повысить качество соединения в тонкостенных конструкциях и воспроизводимость технологического режима сварки изделий. Авторы разработали технологическое решение представленной задачи, опираясь на многолетний опыт по использованию моделей тепловых процессов, сопровождаемых при электронно-лучевой сварке. В качестве предмета исследования предложено моделирование параметров электронно-лучевой сварки тонкостенных труб для волноводных трактов, изготавливаемых из алюминиевого сплава АД31. В статье приведены результаты математического моделирования технологических параметров при нагреве алюминиевого сплава источниками энергии, эквивалентными электронному пучку при ЭЛС. Анализ и оценка результатов моделирования проводилась с применением критерия оптимальности, разработанного авторами. В основе расчетов, выполненных авторами, применен функционал, использующий математические модели нагрева металла комплексным источником тепла, состоящим из подвижных мгновенных точечных и линейного источников энергии. В статье приведены результаты расчетов для пластины толщиной 0,12 см, которая соответствует опытным образцам, применяемых при изготовлении волноводных трактов. В результате путем изменения таких значений, как: ток пучка и скорость сварки, было получено распределение температуры на поверхности изделия в процессе ЭЛС, что показало, применимость моделирования для отработки нового технологического процесса.

Ключевые слова

электронно-лучевая сварка, моделирование, технологические параметры, программа, оптимизация, распределение энергии

Полный текст

Введение

В настоящее^¹ время при производстве конструкций волноводных трактов применяются разнообразные технологии, несущие большие затраты вследствие низкой воспроизводимости качества соединений конструкций [1–3]. Малая толщина и низкая температура плавления соединяемых изделий ставит под сомнение целесообразность применения электронно-лучевой сварки в этом направлении [4–6]. В работах [7–9] авторы моделируют процесс электронно-лучевой сварки (ЭЛС) для титанового сплава с целью нахождения оптимального режима. Полученный опыт моделирования привел авторов к идее исследования возможности применения электронно-лучевой сварки тонкостенных алюминиевых конструкций. Таким образом была поставлена цель смоделировать процесс сварки исследуемого материала с заданной толщиной и сравнить результаты моделирования с возможностями электронно-лучевого оборудования.

За основу в исследованиях была принята классическая теория тепловых процессов [10; 11]. Распределение температур в тонкослойном материале при действии на его поверхности подвижного точечного источника может быть рассчитано с использованием метода суперпозиции (сложения источника) [12–15]. Действительный точечный источник, изображенный на рис. 1, принимают перемещающимся по поверхности полубесконечного тела.

Отражение потока теплоты, отдаваемая источником Q₁, от границы II учитываем введением источника Q₄, симметричного источнику Q₁ относительно границы II, т.е. на расстоянии 2δ от источника Q₁, действующего на границе I. В свою очередь, граница I для источника Q₄ будет учтена, если ввести фиктивный источник Q₃, удаленный от границы I на расстояние 2δ. Для учета толщины материала δ при моделировании теплового процесса был введен линейный движущийся источник Q₂.

Рис. 1. Схема введения фиктивных источников для расчета температур в плоском слое

Fig. 1. Scheme of introducing fictitious sourcesfor calculating temperatures in a flat layer

Математический аппарат

Приняв выбранный комплекс источников нагрева и применив метод наложения [10], получаем итоговую модель температуры нагрева (1):

$T_{\sum} = T_{1} + T_{2} + T_{3} + T_{4}$ , (1)

где $T_{1} = T_{H} + \frac{Q_{1}}{cρ \sqrt{{(4 π a)}^{3}}} e^{- \frac{v x}{2 a}}_{0}^{t} exp (- \frac{v^{2} τ}{4 a} - \frac{x^{2} + y^{2} + z^{2}}{4 a τ}) \frac{d τ}{τ^{3 / 2}}$ – температура нагрева от движущегося по поверхности мгновенного точечного источника;

$T_{2} = T_{H} + \frac{Q_{2}}{4 πλδ} e^{- \frac{v x}{2 a}}_{0}^{t} exp (- \frac{v^{2} τ}{4 a} - \frac{2 λτ}{cρδ} - \frac{(x^{2} + y^{2})}{4 a τ}) \frac{d τ}{τ}$ – температура от нагрева движущегося линейного источника;

$T_{3} = T_{H} + \frac{Q_{3}}{cρ \sqrt{{(4 π a)}^{3}}} e^{- \frac{v x}{2 a}}_{0}^{t} exp (- \frac{v^{2} τ}{4 a} - \frac{x^{2} + y^{2} + {(z - 2 δ)}^{2}}{4 a τ}) \frac{d τ}{τ^{3 / 2}}$ – температура нагрева от движущегося над поверхностью фиктивного мгновенного точечного источника;

$T_{4} = T_{H} + \frac{Q_{4}}{cρ \sqrt{{(4 π a)}^{3}}} e^{- \frac{v x}{2 a}}_{0}^{t} exp (- \frac{v^{2} τ}{4 a} - \frac{x^{2} + y^{2} + {(z + 2 δ)}^{2}}{4 a τ}) \frac{d τ}{τ^{3 / 2}}$ – температура нагрева от движущегося над поверхностью фиктивного мгновенного точечного источника.

При моделировании теплового процесса авторами используется разработанный ими функционал:

$J = \sqrt{\frac{1}{n - 1}_{i = 1}^{n} {(T_{н i} - \bar{T_{н}})}^{2}}$ , (2)

где $T_{н i} = \frac{T_{\sum}}{T_{max}}$ ; $T_{max} = max T_{\sum}$ .

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

Рис. 2. График изменения функционала от скорости перемещения источника нагрева

Fig. 2. Graph of the change in the functional from the speed of movement of the heating source

Расчеты температуры нагрева (2) проводились при значениях времени интегрирования, соответствующих моменту достижения границы зоны термического влияния нижней поверхности нагреваемого материала, пример показан на рис. 3.

В качестве границы зоны термического влияния принято максимальное значение температуры, которой соответствуют справочные [11] свойства исследуемого материала (c, ρ, a, λ).

Рис. 3. Моделирование поперечного сечения зоны термического влияния:h – толщина сварного шва; x – ширина изделия

Fig. 3. Modeling a Cross-Section of a Heat-Affected Zone:h – the thickness of the weld; x – product width

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

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

Результаты

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

Для наглядности получаемой при моделировании информации результаты вычислений сведены в таблицу.

Результаты вычисления модели (2)

1	2	3	4
Q, [калл]	J	v, [см/с]	t, [с]
3	0,0979	0,21	1,107
4	0,0958	0,27	0,869
5	0,0939	0,32	0,721
6	0,0922	0,36	0,619
7	0,0907	0,40	0,545
8	0,0892	0,44	0,489
9	0,0874	0,49	0,446
10	0,0865	0,52	0,409
11	0,0857	0,56	0,379
12	0,0837	0,60	0,355
13	0,0835	0,64	0,333
14	0,0815	0,67	0,314
15	0,0806	0,70	0,298
16	0,0809	0,74	0,284
17	0,0798	0,77	0,271
18	0,0795	0,80	0,259
19	0,0787	0,84	0,249
20	0,0782	0,86	0,239
21	0,0774	0,90	0,231
22	0,0770	0,93	0,223
23	0,0764	0,96	0,216
24	0,0762	0,99	0,209
25	0,0755	1,02	0,203
26	0,0742	1,05	0,198
27	0,0744	1,07	0,192
28	0,0743	1,11	0,187
29	0,0727	1,13	0,183
30	0,0733	1,16	0,178
31	0,0725	1,18	0,174
32	0,0727	1,22	0,170
33	0,0713	1,24	0,167
34	0,0718	1,27	0,163
35	0,0709	1,29	0,160
36	0,0705	1,32	0,157
37	0,0704	1,35	0,154
38	0,0702	1,37	0,151
39	0,0695	1,41	0,149
40	0,0693	1,42	0,146
41	0,0686	1,45	0,144
42	0,0690	1,47	0,141
43	0,0685	1,50	0,139
44	0,0680	1,52	0,137
45	0,0679	1,55	0,135
46	0,0676	1,57	0,133
47	0,0677	1,60	0,131
48	0,0676	1,62	0,129
49	0,0678	1,64	0,127
50	0,0667	1,67	0,126

В таблицу помещены численные значения заданного комплексного источника нагрева Q, согласно уравнению (1). Численные значения скорости движущегося источника нагрева v и времени интегрирования t выбраны по аналогии с рис. 2 и 3 соответственно. Расчеты проводились с применением упрощенной модели (2), то есть варьировались координаты y и z, а координата x выбрана по максимальной площади зоны термического влияния. При таком выборе параметров моделирования согласно результатам, полученным авторами ранее [8–10], форма зоны термического влияния максимально приближена к форме поперечного сечения сварного шва. Численные значения функционала J в таблице определены по уравнению (2) для всего объема материала, ограниченного зоной термического влияния. Авторы допустили наличие корреляции между значениями всех столбцов таблицы.

По результатам, представленным в таблице, авторами получен график зависимости функционала J от значения вводимой энергии комплексного источника нагрева, показанный на рис. 4. Если принять за основу гипотезу о пропорциональности между качеством сварного шва и численного значения функционала (2), то согласно полученной зависимости, изображенной на рис. 4, авторами сделано следующее предположение о том, что при сварке тонкостенных конструкций не требуется оборудование, обеспечивающее высокостабильный очень малый ток сварки.

Рис. 4. График зависимости функционала от энергии источника нагрева

Fig. 4. Graph of the dependence of the functional on the energy of the heating source

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

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

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

Рис. 5. График зависимости энергии источника нагрева от скорости его движения

Fig. 5. Graph of the dependence of the energy of the heating source on the speed of its movement

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

Заключение

В статье разработан новый подход к поиску технологических режимов в электронно-лучевой сварке. В ходе исследования авторы провели моделирование тонкостенной пластины из алюминиевого сплава АД31 толщиной 0,12 см. В результате путем изменения таких значений, как: ток пучка и скорость сварки, было получено распределение температуры на поверхности изделия в процессе электронно-лучевой сварки, что показало, что можно отработать новый технологический процесс с помощью моделирования.

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

___________________________________

^[1] Исследование выполнено при финансовой поддержке РФФИ, Правительства Красноярского края и Краевого фонда науки в рамках научного проекта № 20-48-242917 «Модели и методы управления процессом электронно-лучевой сварки тонкостенных конструкций».

The reported study was funded by Russian Foundation for Basic Research, Government of Krasnoyarsk Territory, Krasnoyarsk Regional Fund of Science, to the research project: “Models and methods for controlling the process of electron beam welding of thin-walled structures”, project No 20-48-242917.