Научная статья на тему 'Разработка математической модели многопроходной дуговой сварки'

Разработка математической модели многопроходной дуговой сварки Текст научной статьи по специальности «Физика»

CC BY
675
144
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ / МНОГОПРОХОДНАЯ СВАРКА / ТЕРМИЧЕСКИЙ ЦИКЛ / ЗОНА ТЕРМИЧЕСКОГО ВЛИЯНИЯ / МИКРОСТРУКТУРА / MATHEMATICAL MODEL / MULTI-PASS ARC WELDING / THERMAL CYCLE / HEAT-AFFECTED ZONE / MICROSTRUCTURE

Аннотация научной статьи по физике, автор научной работы — Макарчук Александра Васильевна, Макарчук Наталия Васильевна, Старцев Василий Николаевич

В статье рассматриваются подходы к созданию модели прогнозирования микроструктуры металла и механических свойств сварного соединения при многопроходной сварке. Основной составляющей этой модели является математический расчет термодинамических циклов проходов сварки для конкретной точки сварного шва и зоны термического влияния. На основе физических процессов, сопровождающих сварку, разработана математическая модель многопроходной дуговой сварки с учетом фазовых переходов при плавлении, затвердевании и испарении металла. Введение в разработанную модель функции состояния металла(доля металла в твердом состоянии) позволила учесть выделение и поглощение скрытой теплоты фазового перехода при плавлении и затвердевании в источниковом члене уравнения теплопроводности. Разработанная квазистационарная модель позволяет установить влияние режима сварки, последовательности укладки сварочных валиков, формы и размеров изделия и условий подогрева на термический цикл, от которого зависят структура и свойства основного металла в зоне термического влияния. Кроме того, решение уравнения теплопроводности в данной форме позволяет легко переходить от температурных зависимостей по координате к зависимостям по времени. Разработанная математическая модель применена к расчету термических циклов при многопроходной сварке листового материала малоуглеродистой низколегированной стали толщиной 20 мм с V-образной разделкой кромок стыка. Результаты расчета позволили определить поперечную форму сварного соединения и положение основных температурных линий в зоне шва и зоны термического влияния. Определены основные количественные параметры термических циклов при многопроходной сварке, такие как максимальные температуры, время пребывания и скорости нагрева и охлаждения металла шва и зоны термического влияния в интервалах температур превращения аустенита. При многопроходной сварке каждая точка шва и зоны термического влияния проходит большое количество разных термических циклов. Такие повторные термические циклы воздействуют на металл ранее выполненных проходов шва и зоны термического влияния. В результате сварное соединение представляет собой агрегат слоев с неоднородной структурой и механическими свойствами. На основании рассчитанных термических циклов прогнозируется микроструктура зоны шва и зоны термического влияния.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Макарчук Александра Васильевна, Макарчук Наталия Васильевна, Старцев Василий Николаевич

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

MATHEMATICAL MODEL ELABORATION OF MULTI-PASS ARC WELDING

The article discusses approaches to the creation of predictive models of metal microstructure and mechanical properties of welded joints in multi-pass welding. The main component of this model is the calculation of thermodynamic cycles of weld passes for a particular point of the weld and heat-affected zone (HAZ). The mathematical model of multi-pass arc welding, based on physical processes, taking into account phase transitions at melting, solidification and evaporation of metal was elaborated. Introduction in the developed model, the function state of the metal (metal part in the solid state) has allowed to take into consideration the release and absorption of latent heat of phase transition during melting and solidification in the source term of the conductivity equation. The developed model allows to evaluate the influence of welding parameters, sequence of welding passes, design and dimensions of the structure and heating conditions on the thermal cycle which determines mechanical properties and structure of the metal in the weld zone and in HAZ. Additionally, calculation the equation of conductivity in this form, moving seamlessly from the temperature’s dependencies by the coordinate into dependencies by time. The developed mathematical model was applied to the calculation of thermal cycles in multi-pass welding sheet material, low carbon and low alloy steel with a thickness of 20 mm with a V-shaped groove of the joint. The calculation results allowed to determine the transverse shape of welded joint and the position of the main temperature lines in the area of the weld and heat-affected zone. The following main parameters of the thermal cycle at multi-pass welding were evaluated: maximal temperatures, holding time and rate of heating and cooling of the metal in the weld zone and in HAZ in temperature transformation range of the austenite. In multi-pass welding each point of the weld and heat-affected zone undergoes a large number of different thermal cycles. Such repeated thermal cycles affect the metal previously completed passes of weld and HAZ. In the result, the welded joint is an assembly of layers with heterogeneous structure and mechanical properties. Microstructure in the weld zone and in HAZ is predicted based on the calculated thermal cycles.

Текст научной работы на тему «Разработка математической модели многопроходной дуговой сварки»

ВЕСТНИК*]

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ВЯцР^

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА^,

DOI: 10.21821/2309-5180-2017-9-1-121-131

MATHEMATICAL MODEL ELABORATION OF MULTI-PASS ARC WELDING

A. V. Makarchuk1, N. V. Makarchuk1, V. N. Startsev2

1 — Admiral Makarov State University of Maritime and Inland Shipping, St. Petersburg, Russian Federation.

2 — Academician I. V. Gorynin Central Research Institute of Structural Materials "Prometey", St. Petersburg, Russian Federation.

The article discusses approaches to the creation of predictive models of metal microstructure and mechanical properties of welded joints in multi-pass welding. The main component of this model is the calculation of thermodynamic cycles of weld passes for a particular point of the weld and heat-affected zone (HAZ). The mathematical model of multi-pass arc welding, based on physical processes, taking into account phase transitions at melting, solidification and evaporation of metal was elaborated. Introduction in the developed model, the function state of the metal (metal part in the solid state) has allowed to take into consideration the release and absorption of latent heat of phase transition during melting and solidification in the source term of the conductivity equation. The developed model allows to evaluate the influence of welding parameters, sequence of welding passes, design and dimensions of the structure and heating conditions on the thermal cycle which determines mechanical properties and structure of the metal in the weld zone and in HAZ. Additionally, calculation the equation of conductivity in this form, moving seamlessly from the temperature's dependencies by the coordinate into dependencies by time. The developed mathematical model was applied to the calculation of thermal cycles in multi-pass welding sheet material, low carbon and low alloy steel with a thickness of 20 mm with a V-shaped groove of the joint. The calculation results allowed to determine the transverse shape of welded joint and the position of the main temperature lines in the area of the weld and heat-affected zone. The following main parameters of the thermal cycle at multi-pass welding were evaluated: maximal temperatures, holding time and rate of heating and cooling of the metal in the weld zone and in HAZ in temperature transformation range of the austenite. In multi-pass welding each point of the weld and heat-affected zone undergoes a large number of different thermal cycles. Such repeated thermal cycles affect the metal previously completed passes of weld and HAZ. In the result, the welded joint is an assembly of layers with heterogeneous structure and mechanical properties. Microstructure in the weld zone and in HAZ is predicted based on the calculated thermal cycles.

Keywords: mathematical model, multi-pass arc welding, thermal cycle, heat-affected zone, microstructure.

For citation:

Makarchuk, Aleksandra V., Nataliya V. Makarchuk, and Vasilii N. Startsev. "Mathematical model elaboration of multi-pass arc welding." Vestnik Gosudarstvennogo universiteta morskogo i rechnogo flota imeni admirala S.O. Makarova 9.1 (2017): 121-131. DOI: 10.21821/2309-5180-2017-9-1-121-131.

УДК 621.791

РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ МНОГОПРОХОДНОЙ ДУГОВОЙ СВАРКИ

А. В. Макарчук1, Н. В. Макарчук1, В. Н. Старцев2

1 — ФГБОУ ВО «ГУМРФ имени адмирала С. О. Макарова», Санкт-Петербург, Российская Федерация.

2 — ГНЦ ФГУП «ЦНИИ КМ «Прометей» НИЦ «Курчатовский институт», Санкт-Петербург, Российская Федерация.

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

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

(доля металла в твердом состоянии) позволила учесть выделение и поглощение скрытой теплоты фазового перехода при плавлении и затвердевании в источниковом члене уравнения теплопроводности. Разработанная квазистационарная модель позволяет установить влияние режима сварки, последовательности укладки сварочных валиков, формы и размеров изделия и условий подогрева на термический цикл, от которого зависят структура и свойства основного металла в зоне термического влияния. Кроме того, решение уравнения теплопроводности в данной форме позволяет легко переходить от температурных зависимостей по координате к зависимостям по времени. Разработанная математическая модель применена к расчету термических циклов при многопроходной сварке листового материала малоуглеродистой низколегированной стали толщиной 20 мм с V-образной разделкой кромок стыка. Результаты расчета позволили определить поперечную форму сварного соединения и положение основных температурных линий в зоне шва и зоны термического влияния. Определены основные количественные параметры термических циклов при многопроходной сварке, такие как максимальные температуры, время пребывания и скорости нагрева и охлаждения металла шва и зоны термического влияния в интервалах температур превращения аустенита. При многопроходной сварке каждая точка шва и зоны термического влияния проходит большое количество разных термических циклов. Такие повторные термические циклы воздействуют на металл ранее выполненных проходов шва и зоны термического влияния. В результате сварное соединение представляет собой агрегат слоев с неоднородной структурой и механическими свойствами. На основании рассчитанных термических циклов прогнозируется микроструктура зоны шва и зоны термического влияния.

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

Для цитирования:

Макарчук А. В. Разработка математической модели многопроходной дуговой сварки / А. В. Макарчук, Н. В. Макарчук, В. Н. Старцев // Вестник Государственного университета морского и речного флота имени адмирала С. О. Макарова. — 2017. — Т. 9. — № 1. — С. 121-131. DOI: 10.21821/2309-51802017-9-1-121-131.

Введение

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

[1] - [3].

В связи с этим большой интерес представляет моделирование физических процессов, происходящих при сварке. При этом моделирование теплового состояния объекта сварки является наиболее востребованным, так как температурное поле определяет форму осей кристаллитов, а, следовательно, и микроструктуру металла шва, что влияет на характер распределения деформа-^ ций и деформационную способность металла в температурном интервале хрупкости. Теоретиче-г ское исследование процессов теплообмена в настоящее время в значительной степени базируется ^ на их численном моделировании с использованием ЭВМ [3] - [6].

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

Математическая модель многопроходной дуговой сварки

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

флюса к поверхности металла в разделке подается сварочная проволока, между ней и металлом горит электрическая дуга. Сварочная проволока движется вдоль кромок по разделке основного металла. Расплавленный металл проволоки и основного металла формирует валик наплавленного металла и нагревает листы основного металла. Валики накладываются последовательно друг за другом в определенном порядке, заполняя разделку листов.

Сформулируем математическую модель распространения тепла в металле в соответствии с описанной физической картиной. Так как предполагается произвольная сборка листов металла, будем искать решение в трехмерной области ХУ2 (X — продольное направление, У — ширина, 2 — высота). Движением жидкого металла в ванне расплава пренебрегаем. Физические параметры металла, такие как плотность, теплоемкость, теплопроводность, являются функциями температуры. Тогда уравнение нестационарной теплопроводности и движения электрода можно записать в следующем виде:

р(Т) • с, (Т). ^ - ¿IV ((Т) • grad(T)) = 5; (1)

- = V , (2)

Ж св 4 7

где р(Г) — плотность; Ср (Т) — теплоемкость; ЦТ) — теплопроводность; £ — объемные источники тепла; Т — температура; Vсв — скорость сварки; t — время; х — продольное направление сварки.

Введем в рассмотрение функцию фазового состояния металла (долю металла в твердом состоянии):

Solid = 0,5 • erfc

4 (T - Tl ) + 2

Tl - TS

(3)

где TS — температура солидуса; TL — температура плавления.

Эта функция равна единице, если T < TS, и равна нулю, если T > TL, а в интервале температур TS < T < TL принимает значение от 1 до 0. При фазовом переходе от твердого состояния к жидкому или от жидкого к твердому поглощается или выделяется скрытая теплота плавления. Если предположить, что количество тепла, выделившееся при фазовом переходе, пропорционально количеству расплавленного металла, то

1

Q = \L • d(Solid), (4)

0

где L — скрытая теплота плавления; Q — теплота.

Тогда объемный источник тепла в уравнении (1) будет иметь вид

5 = p(T).L (5)

dt

Перейдем теперь к системе координат, связанной с электродом. Тогда вместо системы уравнений (1) - (2) получим одно стационарное уравнение:

Р(Т) • с, (T) • Усв • <dT = div(X(T) • grad(T)) + p(T) • L ■ VCT ■ d(SOid). (6)

dx dx

Для задания граничных условий построим методику расчета наплавки одного валика. Разделка кромок металла по высоте делится в соответствии с чертежами на сварку на n слоев. Причем высота предыдущего слоя является верхней поверхностью внутри разделки. В точке с координатами Х0, Y0, Z0 на площади, ограниченной радиусом r, действует поверхностный источник тепла с заданным распределением плотности мощности q. Высота металла валика шириной £вал увеличивается от этой координаты на длине L от высоты предыдущего слоя валиков до высоты следующего слоя, и далее до конца разделки по длине высота валика остается постоянной. Считается, что металл непрерывно поступает в плоскости Y - Z в начале координат и на участке роста высоты

2 О

7

валика в направлении оси Х со скоростью 'св. После расчета температурного поля переходят к рас -чету следующего валика и так далее до полного заполнения разделки.

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

В начале расчетной области в плоскости У - 2 и на участке роста высоты валика задается начальная температура и скорость сварки:

Т = Т

V = V .

св

'0' V 'св. (7)

На верхней поверхности листов основного металла в полосе шириной h в обе стороны от оси разделки, которая покрыта флюсом, задается условие теплоизоляции

йТ

йп

= 0,

(8)

где п — направление нормали.

На остальной части верхней поверхности, боковой и нижней части листов задается условие конвективного теплообмена

Х- — = -а-к-(Т- Т ),

йп У воз;'

(9)

где а — коэффициент теплоотдачи; k — коэффициент, учитывающий направление нормали к направлению силы тяжести.

В области воздействия теплового источника задан тепловой поток и учитывается его изменение, связанное с фазовыми переходами плавления и испарения металла:

X • — = q - р • V • L - т • LL,

(10)

где q — плотность теплового потока сварочного источника; L — скрытая теплота плавления; LL — скрытая теплота испарения; т — поток пара металла с единицы поверхности.

Следует отметить, что учет испарения металла с поверхности ванны расплава является важной составляющей баланса тепла, так как при плотности теплового потока, превышающей 104 Вт/см2, температура поверхности сварочной ванны в области воздействия дуги будет превышать температуру кипения от нескольких десятков до сотен градусов. Расчет потока испаряющегося металла с поверхности ванны расплава проводился в рамках одномерной модели испарения Найта [7].

Таким образом, решение уравнений (3), (6) совместно с граничными условиями (7) - (10) определяют распределение температуры при наложении одного валика металла при многопроходной сварке.

Особенностью данного расчета является то, что температурные поля в различных сечениях по координате Х соответствуют температурным полям для одного сечения, но при разном времени. При этом время определяется как отношение координаты к скорости сварки, а скорость нагрева и охлаждения металла — как произведение скорости сварки на производную температуры по координате Х (11):

Г, =

X

V

ат_

Ж

ат

(11)

ЕЭ

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

ВЕСТН1

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

Результаты численного моделирования многопроходной сварки

Разработанная математическая модель была применена к расчетной оценке параметров термического цикла сварного шва и зоны термического влияния при сварке листовой малоуглеродистой низколегированной стали толщиной 20 мм. Рассматривается Г-образная разделка кромок с углом 45° и притуплением 1 мм. Раскладка валиков внутри разделки листов металла осуществлялась в шесть слоев. Расчет проводился только в разделке, области усиления сварных швов не рассчитывались. Теплофизические характеристики в расчетах были следующие: теплоемкость С = 690 Дж/кгК, теплопроводность X = 0,268 Вт/смК, плотность р = 0,0078 кг/см3, температура плавления Тпл = 1500 °С. Скорость сварки в расчетах варьировалась в диапазоне Гсв = 10 - 20 мм/с. В пятне нагрева радиусом г = 6,5 мм плотность мощности принималась q = 18 кВт/см2. Начальная температура листов металла для проведенных расчетов принималась 20 °С.

Работоспособность сварного соединения при действии внешних нагрузок зависит от локальной микроструктуры и локальных механических свойств сварного шва, зоны термического влияния и основного металла. Классическая схема структурных изменений металла шва и зоны термического влияния основана на схеме однопроходной сварки, в которой каждая зона проходит один раз свой термический цикл. Максимальная температура цикла, скорость нагрева в интервале температур от температуры критической точки Ас3 до Ттах, время пребывания металла выше критической точки Ас3 соответственно при нагреве, охлаждении и суммарное, скорости охлаждения в интервалах температур превращения аустенита 800 °С - 500 °С и 600 °С - 500 °С, а также соответствующие этим интервалам температур длительности охлаждения; мгновенная скорость охлаждения при температуре наименьшей устойчивости аустенита ТтП являются важнейшими параметрами термических циклов сварки [3], [4], [8], [9]. Характер изменений структуры зависит от химического состава основного металла и его предшествующей термической и механической обработки. Здесь необходимо отметить, что при многопроходной сварке каждая точка шва и зоны термического влияния проходит большое количество разных термических циклов. Такие повторные термические циклы воздействуют на металл ранее выполненных проходов шва и зоны термического влияния. В результате сварное соединение представляет собой агрегат слоев с неоднородной структурой и механическими свойствами (рис. 1). На фотографии шлифа хорошо видны зоны термического влияния от валиков, наложенные на зоны термического влияния ранее наплавленных валиков.

Рис. 1. Макроструктура сварного шва и зоны термического влияния с Х-образной разделкой кромок листа, выполненного MIG/MAG сваркой

,ВЕСТНИК

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

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

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

Рис. 2. Расчетная форма линий сплавления при многопроходной сварке с К-образной разделкой кромок листа толщиной 20 мм (номера валиков не связаны с очередностью наплавки валиков)

Рассмотрим термические циклы небольшой области металла вокруг точки с координатами У = 1,0 см, 2 = 0,8 см вблизи от границы сплавления валика 2 и основного металла (рис. 3).

e>j

В26

30. L (см)

0. 5. 10. 15. 20. t (с)

Рис. 3. Термические циклы при многопроходной сварке для малой области с координатами центра Y = 1,0 см, Z = 0,8 см со скоростью V = 15 мм/с (циклы приведены к единому началу)

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

при наложении различных сварочных валиков в разные зоны термического влияния. Рассматриваемая область металла нагревается до максимальной температуры от валика 2 и достигает величины чуть больше Ттах = 1400 К (~1130 °С), которая несколько превышает температуру Ас3 а^у превращения, т. е. находится на участке перекристаллизации (нормализации). Время пребывания металла при температуре выше критической точки Ас3 составляет А= 2,8 с, из них нагрев осуществляется в течение А^ = 0,44 с со средней скоростью дТ/д t = 743 К/с, охлаждение метала до точки Ас3 проходит со средней скоростью дТ/дtохл = -138 К/с. Длительности охлаждения в интервалах температур превращения аустенита 1073 - 773 К (800 °С - 500 °С) и 873 - 773 К (600 °С - 500 °С), а также соответствующие этим интервалам температур скорости охлаждения, соответственно, равны Аt8/5 = 4,5 с, дT/дt8/5 = -67 К/с и Аtб/5 =1,9 с, дТ/д t6/5 = -52 К/с. Мгновенная скорость охлаждения в точке минимальной устойчивости аустенита равна дТ/д = -97 К/с. При выполнении этого прохода металл переходит в состояние аустенита, который при последующем охлаждении распадается, и можно ожидать образования мартенситно-бейнитной микроструктуры. В малоуглеродистой низколегированной стали на участке нормализации образуется мелкозернистая структура, характеризующаяся в целом достаточно высоким комплексом механических свойств (прочностью, пластичностью, вязкостью). На закаливающихся сталях после данного прохода металл в данной точке будет соответствовать микроструктуре закалки.

Вторым по значимости термическим циклом для данной области является наплавка валика 5. Максимальная температура, достигаемая при проходе 5, чуть ниже температуры Ас3, Ттах = 1060 К (~790 °С), т. е. наша точка оказывается в области неполной перекристаллизации (см. рис. 3 и 4). Его основными параметрами являются: время пребывания металла в межкритическом интервале температур Аt = 3,55 с, Аt = 1,18 с со средней скоростью дТ/д t = 74 К/с

г г ^ г преб 5 5 нагр 5 г ^ г- нагр

и охлаждение после достижения Ттах со средней скоростью дТ/д ^хл = -40 К/с. В интервале температур распада аустенита параметры термического цикла: А^/5 = 5,9 с, дТ/д ts/5 = -51 К/с и А^/5 = 2,0 с, дТ/д ^/5 = -50 К/с. Нагрев и охлаждение малоуглеродистой низколегированной стали по термическому циклу 5 не изменяет его структуру. Она остается мартенситно-бейнитной, но размер зерна может измениться. Качество сварного соединения значительно меньше зависит от изменений структуры при наплавке валика 5, чем это было при наложении валика 2.

2 О

7

[127 I

Рис. 4. Графики термического цикла и скорости нагрева, охлаждения металла в точке У = 1,0 см, 7 = 0,8 см при наплавке валика 5

При наплавке валика 8 максимальная температура в рассматриваемой точке достигает Ттах = 900 К (~630 °С). Это участок рекристаллизации, на котором металл, нагревается от 800 К (~500 °С) до температуры несколько ниже Ас1. Время достижения максимальной температуры от начала нагрева составляет около 5 с. После начала охлаждения, рассматриваемая область сразу попадает в зону интервала температур 873 - 773 К (600 °С - 500 °С), параметры которой Д* = 3,53 с, дТ/д * = -28 К/с. Структура и свойства металла при нагреве в этом диапазоне могут быть различными в зависимости от исходного его состояния перед сваркой. Холодная пластическая деформация металла перед сваркой при нагреве до температур ниже Ас1 вызывает рекристаллизацию к значительному росту зерна.

Термический цикл от валика 9 проходит по нижней границе участка рекристаллизации Ттах = 770 К (~500 °С). Максимальная температура достигается за 6,6 с. Скорость охлаждения на начальном участке составляет дТ/д * = -20 К/с. Поэтому структура металла, полученная после наплавки валиков 2, 5 и 8, после термического цикла валика 9, как и от всех остальных валиков, термические циклы которых располагаются ниже его, в процессе сварки не будет претерпевать видимых изменений.

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

Аналогичным образом можно получить термическую картину для любой точки сварного шва и зоны термического влияния — например для области разделки листов, где материалом заполнения в основном является сварочная проволока. Рассмотрим область с центром в точке с координатами Y = 1,2 см, 2 = 1,2 см. Термические циклы, которые испытывает металл в зоне шва, носят особый характер, так как часть из них имеет максимальную температуру выше температуры плавления, а другая часть меньше нее (рис. 5). Время пребывания металла выше точки Ас3 достигает от 4 - 5 с при скорости сварки 20 мм/с и до 8 - 10 с при скорости 10 мм/с. Максимальные температуры нагрева в этой области 2500 - 2800 К свидетельствуют о большом перегреве металла. Повторные термические циклы могут иметь максимальные температуры в диапазоне 500 - 1800 К. После того как металл стал твердым, скорость охлаждения изменяется от 200 - 250 К/с в начале до 10 - 20 К/с в конце. Для повторных термических циклов скорость охлаждения находится в диапазоне 80 - 10 К/с.

Т|

2500

2000

1500

ЕЗ

1000

500

20. |_(СМ)

0. 2.5 5. 7.5 10. Мсек.)

Рис. 5. Графики термических циклов при наплавке валиков 5 и 8 в области зоны шва с координатами центра Y = 1,2 см, 2 = 1,2 см со скоростью сварки V = 20 мм/с

ВЕСТНИК«

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ......^

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА,

Как уже отмечалось, химический состав стали и термический цикл сварки определяет локальные микроструктурные и механические свойства шва и зоны термического влияния. Вместе с тем, взаимозависимость этих связей изучена недостаточно, поэтому на практике используются регрессивные модели образования фаз и механических свойств (твердость HV, временное сопротивление оВ, предел текучести а02, относительное удлинение 5, относительное сужение у, энергию удара KVC) [10] - [12], которые учитывают влияние максимальной температуры термических циклов, времени пребывания и скорости охлаждения в интервалах температур превращения аустенита At8/5, дТ/дt .

Выводы

1. На основе физических процессов, сопровождающих сварку, разработана математическая модель многопроходной дуговой сварки. Введение в разработанную модель функции состояния металла (доля металла в твердом состоянии) позволила учесть выделение и поглощение скрытой теплоты фазового перехода при плавлении и затвердевании металла свариваемой конструкции. Кроме того, в граничных условиях принято во внимание изменение теплового потока дугового источника, связанное с плавлением и испарением наплавляемого металла. Учет скрытой теплоты плавления и затвердевания, которая оказывает влияние на форму и размеры сварочной ванны, позволяет в целом более точно рассчитать температурные поля при сварке, а, следовательно, и термические циклы для конкретных точек сварного соединения.

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

3. Определены важнейшие параметры термических циклов многопроходной сварки: максимальная температура цикла; скорость нагрева в интервале температур от критической точки Ас3 до максимальной температуры; время пребывания металла выше критической точки Ас3 соответственно при нагреве, охлаждении и их сумма; скорости охлаждения в интервалах температур превращения аустенита 800 °С - 500 °С и 600 °С - 500 °С, а также соответствующщая этим интервалам температур длительность охлаждения; мгновенная скорость охлаждения при температуре наименьшей устойчивости аустенита.

СПИСОК ЛИТЕРАТУРЫ

1. Березовский Б. М. Математические модели дуговой сварки: монография / Б. М. Березовский. — Челябинск: Изд-во «Челябинский ЦНТИ», 2006. — Т. 4. Основы тепловых процессов в свариваемых изделиях. — 547 с.

2. Березовский Б. М. Математические модели дуговой сварки: монография / Б. М. Березовский. — Че-лябинск: Изд-во «Челябинский ЦНТИ», 2008. — Т. 5. Температурные поля свариваемых изделий и термические циклы сварных соединений. — 656 с.

3. Кархин В. А. Тепловые процессы при сварке: монография / В. А. Кархин. — СПб.: Изд-во СПбГПУ, 2013. — 646 с.

4. Ерофеев В. А. Физико-математическая модель процесса многопроходной дуговой сварки / В. А. Еро-феев, А. В. Масленников // Известия ТулГУ. Сер. Компьютерные технологии в соединении материалов. — 2005. — № 3. — С. 246-255.

5. Собко С. А. Моделирование тепловых процессов при аргонодуговой сварке / С. А. Собко, И. С. Га-реев, М. С. Писарев, Е. В. Брунеткина // Сварочное производство. — 2013. — № 4. — С. 6-10.

6. Rao R. V. Multi-pass turning process parameter optimization using teaching-learning-based optimization algorithm / R. V. Rao, V. D. Kalyankar // Scientia Iranica. — 2013. — Vol. 20. — Is. 3. — Pp. 967-974. DOI: 10.1016/j.scient.2013.01.002.

7. Кривцун И. В. Численный анализ процессов нагрева и конвективного испарения металла при обработке излучением импульсного лазера / И. В. Кривцун, И. Л. Семенов, В. Ф. Демченко // Автоматическая сварка. — 2010. — № 1. — С. 5-10.

2 О

7

'129 I

ЛВЕСТНИК

............ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

Х^ОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

8. Донченко Б. А. Расчет термодинамических циклов точек при автоматической сварке и наплавке с учетом особенностей плавления основного металла / Б. А. Донченко // Сварочное производство. — 2011. — № 9. — С. 3-9.

9. Донченко Б. А. Особенности расчета максимальных температур при автоматической сварке и наплавке / Б. А. Донченко // Сварочное производство. — 2011. — № 2. — С. 3-7.

10. Karkhin V. A. Prediction of microstructure and mechanical properties of weld metal with consideration for real geometry / V. A. Karkhin, P. N. Khomich, V. G. Michailov // Proceedings of Joint International Conference "Computer Technology in Welding and Manufacturing and Information Technologies in Welding and Related Processes". — Kiev, 2006. — Pp. 162-166.

11. Касаткин О. Г. Расчетная оценка ударной вязкости низколегированного металла шва / О. Г. Касаткин // Автоматическая сварка. — 2005. — № 1. — С. 57-58.

12. Krasovskyy A. On the residual stresses in multi-pass welds: coupling of welding simulation and fatigue analysis / A. Krasovskyy, S. Sonnichsen, D. Bachmann // Procedia Engineering. — 2011. — Vol. 10. — Pp. 506-511. DOI: 10.1016/j.proeng.2011.04.085.

REFERENCES

1. Berezovskii, B. M. Osnovy teplovykh protsessov v svarivaemykh izdeliyakh. Chelyabinsk: Izd-vo «Chelyabinskii TsNTI», 2006.

2. Berezovskii, B. M. Temperaturnye polya svarivaemykh izdelii i termicheskie tsikly svarnykh soedinenii. Chelyabinsk: Izd-vo «Chelyabinskii TsNTI», 2008.

3. Karkhin, V. A. Teplovyeprotsessy pri svarke: monografiya. SPb.: Izd-vo SPbGPU, 2013.

4. Erofeev, V. A., and A. V. Maslennikov. "Fiziko-matematicheskaya model' protsessa mnogoprokhodnoi dugovoi svarki." Izvestiya TulGU. Ser. Komp'yuternye tekhnologii v soedinenii materialov 3 (2005): 246-255.

5. Sobko, S. A., I. S. Gareyev, M. S. Pisarev, and Ye. V. Brunetkina. "Thermal process modelling at argon-arc welding." Svarochnoe Proizvodstvo 4 (2013): 6-10.

6. Rao, R. Venkata, and V. D. Kalyankar. "Multi-pass turning process parameter optimization using teaching-learning-based optimization algorithm." Scientia Iranica 20.3 (2013): 967-974. DOI: 10.1016/j.scient.2013.01.002.

7. Krivtsun, I. V., I. L. Semenov, and V. F. Demchenko. "Chislennyi analiz protsessov nagreva i konvektivnogo ispareniya metalla pri obrabotke izlucheniem impul'snogo lazera." Avtomaticheskaya svarka 1 (2010): 5-10.

8. Donchenko, Ye. A. "Calculation of a dot thermal cycle at machine welding and surfacing taking into consideration melting features of the parent metal." Svarochnoe Proizvodstvo 9 (2011): 3-9.

9. Donchenko, Ye. A. "Calculation features of peak temperatures at automatic welding and deposition." Svarochnoe Proizvodstvo 2 (2011): 3-7.

10. Karkhin, V. A., P. N. Khomich, and V. G. Michailov. "Prediction of microstructure and mechanical properties of weld metal with consideration for real geometry." Proceedings of Joint International Conference "Computer Technology in Welding and Manufacturing and Information Technologies in Welding and Related Processes". Kiev, 2006: 162-166.

11. Kasatkin, O. G. "Raschetnaya otsenka udarnoi vyazkosti nizkolegirovannogo metalla shva." Avtomaticheskaya svarka 1 (2005): 57-58.

12. Krasovskyy, A., S. Sonnichsen, and D. Bachmann. "On the residual stresses in multi-pass welds: coupling of welding simulation and fatigue analysis." Procedia Engineering 10 (2011): 506-511. DOI: 10.1016/j. proeng.2011.04.085

ИНФОРМАЦИЯ ОБ АВТОРАХ INFORMATION ABOUT THE AUTHORS

Макарчук Александра Васильевна — Makarchuk, Aleksandra V. — PhD

кандидат технических наук Admiral Makarov State University

ФГБОУ ВО «ГУМРФ имени of Maritime and Inland Shipping

адмирала С.О. Макарова» 5/7 Dvinskaya Str., St. Petersburg 198035, 198035, Российская Федерация, Санкт-Петербург, Russian Federation

ул. Двинская 5/7 e-mail: makarchukav@gumrf.ru e-mail: makarchukav@gumrf.ru

ВЕСТНИК«)

ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ^^

МОРСКОГО И РЕЧНОГО ФЛОТА ИМЕНИ АДМИРАЛА С. О. МАКАРОВА

Макарчук Наталия Васильевна —

кандидат технических наук ФГБОУ ВО «ГУМРФ имени адмирала С.О. Макарова»

198035, Российская Федерация, Санкт-Петербург,

ул. Двинская 5/7

e-mail: makarchuknv@gumrf.ru

Старцев Василий Николаевич —

кандидат технических наук

ГНЦ ФГУП «ЦНИИ КМ «Прометей»

НИЦ «Курчатовский институт»

191015, Российская Федерация, Санкт-Петербург,

ул. Шпалерная 49

e-mail: star1v@mail.ru

Makarchuk, Nataliya V. — PhD

Admiral Makarov State University

of Maritime and Inland Shipping

5/7 Dvinskaya Str., St. Petersburg 198035,

Russian Federation

e-mail: makarchuknv@gumrf.ru

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Startsev, Vasilii N. — PhD

Academician I.V. Gorynin Central

Research Institute of Structural

Materials "Prometey"

49 Shpalernaya Str., St. Petersburg, 191015,

Russian Federation

e-mail: star1v@mail.ru

Статья поступила в редакцию 20 декабря 2016 г.

Received: December 20, 2016.

i Надоели баннеры? Вы всегда можете отключить рекламу.