Numerical Simulation and Experimental Investigation of Temperature and Residual Stress Distributions in a Circular Patch Welded Structure

: In this study, we performed a numerical simulation and experimental measurements on a steel circular patch welded structure to investigate the temperature and residual stress ﬁeld distributions caused by the application of buried-arc welding technology. The temperature histories during the welding and subsequent cooling process were recorded for two locations, with the thermocouples mounted inside the plate close to the weld bead. On the upper surface of the welded model, the temperature-time changes during the cooling process were monitored using an infrared camera. The numerically calculated temperature values correlated well with the experimentally measured ones, while the maximum deviation of the measured and calculated temperatures was within 9%. Based on the numerical result analysis regarding circumferential and radial stresses after the completion of the welding process, it is concluded that both stresses are primarily tensile within the circular disk. Outside the disk, the circumferential stresses turn from tensile to compressive, while on the other hand the radial stresses disappear towards the ends of the plate.


Introduction
Metal active gas (MAG) welding is one the most reliable procedures for joining various structural components. MAG is applied in many industries due to its low manufacturing cost and fast yet simple performance. While conducting this technique, a large amount of molten wire metal fills the groove between the structural components of the base material. As the heat input of the filler material is very high and its cooling rate is rapid afterwards, permanent plastic deformations and residual stresses occur in the weld area and its vicinity [1][2][3][4][5]. The plastic deformations disturb the external appearance of the welded structure making its assembly difficult. Removing such deformations with mechanical or thermal procedures [6][7][8][9] is often a time consuming task, which significantly impacts the final price of the welded structure. On the other hand, the residual stresses in the weld and the heat affected zone can promote crack initiation, especially in dynamically loaded structures.
Since the occurrence of residual stresses and deformations due to welding is an inevitable phenomena, it is necessary to know their distribution as accurately as possible in order to take adequate measures already at the design stage to reduce them. For this purpose, it is reasonable to apply numerical methods in order to avoid expensive experimental measurements [10][11][12][13][14][15][16]. Regarding the above mentioned residual stresses and deformations caused by MAG welding, there are numerous numerical and experimental studies in the scientific literature for all important forms of welded structures, some of which are elaborated below.
Based on the experimental measurements, Adamczuk et al. [17] presented a methodology for angular deformation predictions for multi pass butt-welded steel plates. Costa et al. [18] investigated residual stresses numerically and experimentally and compared them using the X-ray procedure due to cold-wire and conventional MAG welding. In their work, Heinze et al. [19] studied the influence of thermal material properties on the weld pool size and metallurgical properties of welded joints. Furthermore, Kung et al. [20] numerically simulated and experimentally measured the influence of a welding jig position on the residual stresses and deformations for a multi pass butt-welded sample. Deng et al. [21] took into account geometrical and material nonlinearities and investigated the influence of a weld reinforcement shape on the deformations in the MAG welded plates. In their work, Deng et al. [22] numerically simulated and experimentally measured the deformations of middle-thick T-joint steel plates welded with the MAG procedure. Some possibilities of speeding-up numerical simulations with the application of a combination of shell and solid finite elements on a T-joint welded sample were suggested by Perić et al. [23]. Furthermore, to additionally shorten the welding process simulation time, Perić et al. [24] presented a simplified numerical T-joint model, where only the finite element reactivation or element birth/death technique was applied in the thermal numerical analysis, while the mechanical analysis was performed in one step. Prajadhiana et al. [25] investigated T-joint welding deformations using virtual manufacturing tools and simplified numerical procedures.
Numerical simulations and experimental investigations of deformations and residual stresses are not only limited to the butt-welded and T-joint welded structures described above but are widely performed on large complex MAG welded panel structures that find widespread application in many industries. Deng et al. [26] calculated the welding deformations on a large panel structure using the inherent strain method (ISM). The needed data for ISM calculations were obtained from a small-scale model where the full thermal-elastic-plastic (TEP) simulation procedure was applied. A similar approach was used by Azad et al. [27], Zhang et al. [28], Podder et al. [29], and many other authors. On the other side, Perić et al. [30] presented a simplified TEP procedure for the welding of residual stresses and deformation calculations in large structures. Here, heat input was applied by using the prescribed thermal boundary conditions instead of applying the heat flux to simulate the moving of the electrode in order to speed up the calculation process. Furthermore, Wu et al. [31] and Zhao et al. [32] investigated the residual stresses in MAG butt-welded pipes, while the residual stresses and deformations in MAG welded lap joints are studied by Lin et al. [33].
It is worth noting that all of the above cited studies of residual stresses and deformations were performed on conventional MAG welded structures. One of the main disadvantages of the conventional MAG process is the relatively low deposit rate and depth of penetration, thus the welding of thicker structures requires a large amount of welding passes, which prolongs the product manufacturing process and increases financial costs. With the development of high current MAG technology (also known as buried-arc welding) [34], the number of welding passes can be significantly reduced due to the increased deposit of molten metal. Unlike conventional MAG welded structures, where the residual stress and deformation fields and their magnitudes are well investigated, such studies in the scientific literature data regarding buried-arc welded structures are limited. In their previous studies, Perić et al. [35,36] numerically and experimentally investigated temperature fields, residual stresses, and deformations on butt-weld and T-joint fillet weld samples. In this paper, the numerical and experimental investigations of temperature and residual stress distributions are further extended to circular patch welded structures often used in repair welding applications.

Model Geometry, Welding Conditions and Material Properties
A cylindrical disk of a diameter of 165 mm was welded to a steel plate with 300 mm × 300 mm × 20 mm dimensions, using the buried-arc welding procedure as shown in Figure 1. Such models are common in repair welding [37][38][39], where a small part of a previously damaged material, such as a crack, is replaced with a patch in the form of a circular disk or a rectangular plate. The model geometry and its dimensions are provided in Figure 1a. Before the start of welding, the cylindrical disk and plate with a previously cut hole with a diameter of 171 mm were tack-welded at the bottom side of the sample at positions φ = 0 • , 90 • , 180 • , and 270 • , respectively. The welding was performed in a single pass using a seven-axis robot (from OTC Daihen Corporation, Japan) with a neutral gun inclination and a welding wire in accordance with ISO 14341-A:G 42 4 C/M G3Si1. The welding start/end was at position φ = 270 • , whereby the direction of the welding torch is indicated in Figure 1a. During the welding, the plate and circular disk were kept free without any mechanical fixtures. After welding, a visual testing was conducted according to ISO 17637 norm; no imperfections were found. The material of the welded model was S355J2+N low carbon steel which elemental composition, welding parameters, thermal and mechanical properties are given in Tables 1 and 2

Model Geometry, Welding Conditions and Material Properties
A cylindrical disk of a diameter of 165 mm was welded to a steel plate with 300 mm × 300 mm × 20 mm dimensions, using the buried-arc welding procedure as shown in Figure 1. Such models are common in repair welding [37][38][39], where a small part of a previously damaged material, such as a crack, is replaced with a patch in the form of a circular disk or a rectangular plate. The model geometry and its dimensions are provided in Figure 1a. Before the start of welding, the cylindrical disk and plate with a previously cut hole with a diameter of 171 mm were tack-welded at the bottom side of the sample at positions ϕ = 0°, 90°, 180°, and 270°, respectively. The welding was performed in a single pass using a seven-axis robot (from OTC Daihen Corporation, Japan) with a neutral gun inclination and a welding wire in accordance with ISO 14341-A:G 42 4 C/M G3Si1. The welding start/end was at position ϕ = 270°, whereby the direction of the welding torch is indicated in Figure  1a. During the welding, the plate and circular disk were kept free without any mechanical fixtures. After welding, a visual testing was conducted according to ISO 17637 norm; no imperfections were found. The material of the welded model was S355J2+N low carbon steel which elemental composition, welding parameters, thermal and mechanical properties are given in Tables 1 and 2

Temperature Measurements
Temperature measurement during the welding and subsequent cooling process to the ambient temperature (24 °C) was performed using two NiCr-Ni thermocouples of K-type (TC-101 and TC-102, Figure 1a) mounted in the middle plane of the plate at a depth of 10 mm. The distances of the thermocouples from the center of the circular disk to the thermocouple positions were 99 and 113 mm (Figure 1a, line A-B, ϕ = 90°), respectively. The data transfer from the measuring spots to the computer was enabled via a data logger (from PICO Technology Ltd., UK). The applied thermocouples were able to measure temperatures in an interval of −220 to 1150 °C.
In addition to the thermocouples that measured the temperature responses at particular locations inside the plate, a non-contact temperature measurement was performed on the upper model surface using an infrared (IR) camera. Having completed the welding process, the temperature changes were monitored during the cooling time to ambient temperature. During the measurement, the IR camera was positioned approximately 1.5 m from the upper model surface. The IR camera used in this experiment was Flir SC2000 (from FLIR Systems AB; 320 × 240 pixels, spectral range of 7.5-13 μm) and had a sensitivity of 80 mK at 30 °C, field of view 24° × 18°, minimum focal distance of 0.5 m. and spatial resolution of 1.3 mrad. The thermograms after the weld completion were analyzed using ThermaCAM Researcher and GIMP software.

Temperature Measurements
Temperature measurement during the welding and subsequent cooling process to the ambient temperature (24 °C) was performed using two NiCr-Ni thermocouples of K-type (TC-101 and TC-102, Figure 1a) mounted in the middle plane of the plate at a depth of 10 mm. The distances of the thermocouples from the center of the circular disk to the thermocouple positions were 99 and 113 mm (Figure 1a, line A-B, ϕ = 90°), respectively. The data transfer from the measuring spots to the computer was enabled via a data logger (from PICO Technology Ltd., UK). The applied thermocouples were able to measure temperatures in an interval of −220 to 1150 °C.
In addition to the thermocouples that measured the temperature responses at particular locations inside the plate, a non-contact temperature measurement was performed on the upper model surface using an infrared (IR) camera. Having completed the welding process, the temperature changes were monitored during the cooling time to ambient temperature. During the measurement, the IR camera was positioned approximately 1.5 m from the upper model surface. The IR camera used in this experiment was Flir SC2000 (from FLIR Systems AB; 320 × 240 pixels, spectral range of 7.5-13 μm) and had a sensitivity of 80 mK at 30 °C, field of view 24° × 18°, minimum focal distance of 0.5 m. and spatial resolution of 1.3 mrad. The thermograms after the weld completion were analyzed using ThermaCAM Researcher and GIMP software.

Temperature Measurements
Temperature measurement during the welding and subsequent cooling process to the ambient temperature (24 • C) was performed using two NiCr-Ni thermocouples of K-type (TC-101 and TC-102, Figure 1a) mounted in the middle plane of the plate at a depth of 10 mm. The distances of the thermocouples from the center of the circular disk to the thermocouple positions were 99 and 113 mm (Figure 1a, line A-B, φ = 90 • ), respectively. The data transfer from the measuring spots to the computer was enabled via a data logger (from PICO Technology Ltd., UK). The applied thermocouples were able to measure temperatures in an interval of −220 to 1150 • C.
In addition to the thermocouples that measured the temperature responses at particular locations inside the plate, a non-contact temperature measurement was performed on the upper model surface using an infrared (IR) camera. Having completed the welding process, the temperature changes were monitored during the cooling time to ambient temperature. During the measurement, the IR camera was positioned approximately 1.5 m from the upper model surface. The IR camera used in this experiment was Flir SC2000 (from FLIR Systems AB; 320 × 240 pixels, spectral range of 7.5-13 µm) and had a sensitivity of 80 mK at 30 • C, field of view 24 • × 18 • , minimum focal distance of 0.5 m. and spatial resolution of 1.3 mrad. The thermograms after the weld completion were analyzed using ThermaCAM Researcher and GIMP software.

Thermal Analysis
The numerical simulation was performed following the sequential approach in two separate numerical analyses: thermal and mechanical, using the Abaqus/Standard software. In the thermal analysis, the time-dependent temperature field distribution was first determined, which was afterwards applied as a thermal load in the mechanical analysis. Here, it was assumed that the total heat input from the welding device to the workpiece takes place via melting droplets. The imposed heat flux Q was uniformly distributed over the weld volume and was obtained from the Equation (1): In Equation (1) η, U and I are heat input efficiency, welding voltage, and welding current, respectively. Taking into account the values from Table 2, the total heat flux applied to the weld volume V H calculated according to Equation (1) is Q = 6.18 × 10 10 Jm −3 s −1 . The convective heat transfer coefficient (h c = 10 Wm 2 K −1 ) and emissivity factor (ε = 0.9) were applied as boundary conditions in the thermal analysis on all the outer surfaces of the model. The moving of the heat source was simulated by applying the model change option which is contained as a standard tool in the Abaqus/Standard software. In the numerical simulation, the entire weld volume was divided into 104 element sets, which were deactivated in the first analysis step. In the next 104 steps, the element sets were activated set by set simulating the moving of the heat source and weld filler deposition. In each of the 104 steps, each element set was first added and, after that, the heat flux was imposed, whereby the time duration of each individual step was 1 s. After adding the last set of elements, the moving of the electrode stopped and the cooling process of 7500 s began. Since the tack-welds were very short and their heat input to the model was minor, they were neglected in the numerical simulation. In the thermal numerical simulation first order DC3D8 finite elements were used.

Mechanical Analysis
The calculated thermal field from the numerical analysis was imposed as a thermal load in the mechanical analysis. To shorten the simulation time, this analysis was performed in one step only, without applying the model change option. As mentioned earlier in the text, the circular disk and plate were free welded without any mechanical fixtures, but they were included in the mechanical analysis to prevent the motion of the model as a rigid body, as shown in Figure 1a. The weld filler and base material were considered as homogeneous elastic-perfectly plastic materials that yield according to the von Mises criterion and the associated flow rule [41,42]. Due to the lack of temperature dependent thermal and mechanical data for the weld filler, they were assumed equal to the base metal ones. Deng et al. [43] reported that the influence of the phase transformations on the final residual and deformation field is small in low-carbon steel welded structures; thus, the steel phase transformations were not considered in this study. Furthermore, the creep behavior of the material was not taken into account as the exposure of the material to high temperature cycles was very short. Finally, after introducing the mentioned material behavior simplifications, the total strain increment can be summarized in the form: In Equation (2) dε e , dε p , and dε th represents elastic, plastic, and thermal strain increments, respectively. In the mechanical analysis C3D8I finite elements with incompatible nodes were applied. A more detailed numerical procedure is available in the authors' previous work [30].

Finite Element Mesh
The finite element mesh was modeled with a total of 36,435 finite elements and is shown in Figure 4. The circular disk, plate, and weld filler were modelled as a single part. All the dimensions of the disk, plate, and gap are given in Figure 1. The mesh sensitivity was controlled by applying the submodeling technique [44,45], whereby the areas of high temperature gradients in the weld and its vicinity were carefully controlled. In both the thermal and mechanical analyses, the same mesh of finite elements was used with a conversion of DC3D8 elements from the thermal analysis to C3D8I elements in the mechanical analysis.
Energies 2020, 13, x FOR PEER REVIEW 6 of 11 The finite element mesh was modeled with a total of 36,435 finite elements and is shown in Figure 4. The circular disk, plate, and weld filler were modelled as a single part. All the dimensions of the disk, plate, and gap are given in Figure 1. The mesh sensitivity was controlled by applying the submodeling technique [44,45], whereby the areas of high temperature gradients in the weld and its vicinity were carefully controlled. In both the thermal and mechanical analyses, the same mesh of finite elements was used with a conversion of DC3D8 elements from the thermal analysis to C3D8I elements in the mechanical analysis.

Temperature Field Distribution
The temperature histories recorded on the thermocouple positions TC-101 and TC-102 for the first 500 s after the beginning of the welding process are shown in Figure 5a. Here, it can be seen that the experimental measured peak temperatures reach values of 525 °C and 326 °C, while the numerical calculated ones are slightly lower: 519 °C and 298 °C, respectively. The slope of the temperature history curves show a sharp rise tendency when the electrode passes the measuring spots and were almost identical for the numerical simulation and experimental measurement, which clearly indicates that the actual thermal property values of S355J2+N steel are very close to the data given in Figure 2. Having reached peak temperatures, there is a temperature decrease due to the cooling process. An acceptable match of numerical and experimental results can be seen here, indicating that the temperature independent convective heat transfer coefficient used in the numerical model was a reasonable choice in this case.
The cooling curves for the time intervals of 112 s and 127 s after the beginning of the welding process are plotted in Figure 5b. The results are given for the upper side of the welded model ( Figure  1a, line A-B, ϕ = 90). In Figure 5b, coordinate y = −150 mm denotes point A, while y = 150 mm represents point B (Figure 1a). In Figure 5b, a large difference in peak temperatures for the time of 112 s after the start of welding is noticeable (920 °C vs. 591 °C). The time of 112 s after the beginning of the welding process is in fact only 8 s after the electrode stopped, which was too short to equalize the temperatures throughout the model. On the other hand, the time of 127 s after the start of welding (in this case the time after the electrode stops is 23 s) was long enough for temperature equalization during the cooling process, so the differences in peak temperatures are significantly smaller (635 °C vs. 531 °C). In addition to the cooling curves obtained by the numerical simulation, in Figure 5b

Temperature Field Distribution
The temperature histories recorded on the thermocouple positions TC-101 and TC-102 for the first 500 s after the beginning of the welding process are shown in Figure 5a. Here, it can be seen that the experimental measured peak temperatures reach values of 525 • C and 326 • C, while the numerical calculated ones are slightly lower: 519 • C and 298 • C, respectively. The slope of the temperature history curves show a sharp rise tendency when the electrode passes the measuring spots and were almost identical for the numerical simulation and experimental measurement, which clearly indicates that the actual thermal property values of S355J2+N steel are very close to the data given in Figure 2. Having reached peak temperatures, there is a temperature decrease due to the cooling process. An acceptable match of numerical and experimental results can be seen here, indicating that the temperature independent convective heat transfer coefficient used in the numerical model was a reasonable choice in this case.
The cooling curves for the time intervals of 112 s and 127 s after the beginning of the welding process are plotted in Figure 5b. The results are given for the upper side of the welded model (Figure 1a, line A-B, φ = 90). In Figure 5b, coordinate y = −150 mm denotes point A, while y = 150 mm represents point B (Figure 1a). In Figure 5b, a large difference in peak temperatures for the time of 112 s after the start of welding is noticeable (920 • C vs. 591 • C). The time of 112 s after the beginning of the welding process is in fact only 8 s after the electrode stopped, which was too short to equalize the temperatures throughout the model. On the other hand, the time of 127 s after the start of welding (in this case the time after the electrode stops is 23 s) was long enough for temperature equalization during the cooling process, so the differences in peak temperatures are significantly smaller (635 • C vs. 531 • C). In addition to the cooling curves obtained by the numerical simulation, in Figure 5b  A comparison of numerically calculated and experimentally recorded temperature fields with the use of an IR camera is shown in Figure 6. It is obvious that the numerical results follow the results of the temperature measurements very well, whereas in both cases the temperature range is between approximately 30 °C and 860 °C. It can be seen that the zone of high temperatures is slightly wider in the thermographic measurement than in the numerical simulation. As stated earlier in the text, it was assumed that all the heat from the welding device was transferred to the workpiece via melting droplets, but, in reality, a part of the heat was transferred by radiation to the surface of the workpiece. Since the zone in which the heat was transferred by radiation from the apparatus to the model surface is wider than the weld gap, the electrode trace obtained by the thermographic measurement is somewhat wider than obtained by the numerical simulation.

Residual Stress Field Distribution
The circumferential (tangent to the welding path) residual stress distribution along line C-D (Figure 1a) on the upper side of the welded model is plotted in Figure 7a. It can be seen that these stresses are tensile throughout the entire disk. The tensile circumferential stress reach a peak value of 327 MPa inside the disk on the coordinate x = −59 mm. In the rest of the model, towards the plate edges, the tensile circumferential residual stresses turn into compressive ones. Similar to the circumferential stresses, the radial stresses (perpendicular to the welding path) are primarily tensile throughout the disk area (Figure 7b). Outside the disk, they reach a tensile value of 160 MPa at the  A comparison of numerically calculated and experimentally recorded temperature fields with the use of an IR camera is shown in Figure 6. It is obvious that the numerical results follow the results of the temperature measurements very well, whereas in both cases the temperature range is between approximately 30 • C and 860 • C. It can be seen that the zone of high temperatures is slightly wider in the thermographic measurement than in the numerical simulation. As stated earlier in the text, it was assumed that all the heat from the welding device was transferred to the workpiece via melting droplets, but, in reality, a part of the heat was transferred by radiation to the surface of the workpiece. Since the zone in which the heat was transferred by radiation from the apparatus to the model surface is wider than the weld gap, the electrode trace obtained by the thermographic measurement is somewhat wider than obtained by the numerical simulation. A comparison of numerically calculated and experimentally recorded temperature fields with the use of an IR camera is shown in Figure 6. It is obvious that the numerical results follow the results of the temperature measurements very well, whereas in both cases the temperature range is between approximately 30 °C and 860 °C. It can be seen that the zone of high temperatures is slightly wider in the thermographic measurement than in the numerical simulation. As stated earlier in the text, it was assumed that all the heat from the welding device was transferred to the workpiece via melting droplets, but, in reality, a part of the heat was transferred by radiation to the surface of the workpiece. Since the zone in which the heat was transferred by radiation from the apparatus to the model surface is wider than the weld gap, the electrode trace obtained by the thermographic measurement is somewhat wider than obtained by the numerical simulation.

Residual Stress Field Distribution
The circumferential (tangent to the welding path) residual stress distribution along line C-D (Figure 1a) on the upper side of the welded model is plotted in Figure 7a. It can be seen that these stresses are tensile throughout the entire disk. The tensile circumferential stress reach a peak value of 327 MPa inside the disk on the coordinate x = −59 mm. In the rest of the model, towards the plate edges, the tensile circumferential residual stresses turn into compressive ones. Similar to the circumferential stresses, the radial stresses (perpendicular to the welding path) are primarily tensile throughout the disk area (Figure 7b). Outside the disk, they reach a tensile value of 160 MPa at the

Residual Stress Field Distribution
The circumferential (tangent to the welding path) residual stress distribution along line C-D (Figure 1a) on the upper side of the welded model is plotted in Figure 7a. It can be seen that these stresses are tensile throughout the entire disk. The tensile circumferential stress reach a peak value of 327 MPa inside the disk on the coordinate x = −59 mm. In the rest of the model, towards the plate edges, the tensile circumferential residual stresses turn into compressive ones. Similar to the circumferential stresses, the radial stresses (perpendicular to the welding path) are primarily tensile throughout the disk area (Figure 7b). Outside the disk, they reach a tensile value of 160 MPa at the coordinate x = −99 mm and decrease towards the edges of the plate and decreased where they disappear completely. The compressive radial residual stresses have a peak value of −82 MPa at the coordinate x = 84 mm. In Figure 7a,b coordinate x = −150 mm denotes point C, while x = 150 mm represents point D (Figure 1a). Looking at the full field circumferential and radial residual stress distributions (Figure 8) on the upper side of the welded model, it is noticeable that they are almost symmetric with respect to the x and y axis (Figure 1a). The residual stress symmetries also show a steady state tendency following the path of the electrode, only becoming lost in the welding start/end area, which was also noticed in previous research by other authors [46]. Upon analyzing both circumferential and radial residual stresses, it can be noted that both stresses are primarily tensile inside the circular disk. Outside the disk, the circumferential stresses turn from tensile to compressive, while on the other hand the radial stresses disappear towards the ends of the plate.

Conclusions
In this study, we investigated temperature and residual stress field distributions in a circular patch weld model. The main conclusions can be summarized as follows: Looking at the full field circumferential and radial residual stress distributions (Figure 8) on the upper side of the welded model, it is noticeable that they are almost symmetric with respect to the x and y axis (Figure 1a). The residual stress symmetries also show a steady state tendency following the path of the electrode, only becoming lost in the welding start/end area, which was also noticed in previous research by other authors [46]. Upon analyzing both circumferential and radial residual stresses, it can be noted that both stresses are primarily tensile inside the circular disk. Outside the disk, the circumferential stresses turn from tensile to compressive, while on the other hand the radial stresses disappear towards the ends of the plate. Looking at the full field circumferential and radial residual stress distributions (Figure 8) on the upper side of the welded model, it is noticeable that they are almost symmetric with respect to the x and y axis (Figure 1a). The residual stress symmetries also show a steady state tendency following the path of the electrode, only becoming lost in the welding start/end area, which was also noticed in previous research by other authors [46]. Upon analyzing both circumferential and radial residual stresses, it can be noted that both stresses are primarily tensile inside the circular disk. Outside the disk, the circumferential stresses turn from tensile to compressive, while on the other hand the radial stresses disappear towards the ends of the plate.

Conclusions
In this study, we investigated temperature and residual stress field distributions in a circular patch weld model. The main conclusions can be summarized as follows:

Conclusions
In this study, we investigated temperature and residual stress field distributions in a circular patch weld model. The main conclusions can be summarized as follows: • Temperature histories recorded by thermocouples match well with the numerical calculated ones. The effect of simplified heat flux affects the temperature distribution in the weld bead and its close vicinity, while this effect on the temperature distribution away from the weld bead is significantly smaller. • By measuring temperatures with the application of an IR camera on the upper surface of the model, a good correlation with the numerical simulation results was confirmed, whereas the largest deviations of the results were found in the weld bead and its vicinity.

•
The numerical simulation residual stress results showed that the circumferential and radial residual stresses within the circular disk are primarily tensile.

•
The circumferential residual stresses turn from tensile to compressive outside the disk. The radial residual stresses disappear towards the ends of the plate.

•
The residual stress symmetries also show a steady state tendency following the path of the electrode, except in the welding start/end area.
Finally, it can be concluded that the presented numerical model is an acceptable solution for the analysis of temperature fields and residual stresses in the case of thick plate repair welding. Future research should focus on investigating the parameters of fracture mechanics and fatigue caused by buried-arc welding.