Next Article in Journal
Composition of Corroded Reinforcing Steel Surface in Solutions Simulating the Electrolytic Environments in the Micropores of Concrete in the Propagation Period
Next Article in Special Issue
Mesoscopic Fluid-Particle Flow and Vortex Structural Transmission in a Submerged Entry Nozzle of Continuous Caster
Previous Article in Journal
Impact Abrasive Wear Property of CrAlN/TiSiN Multilayer Coating at Elevated Temperatures
Previous Article in Special Issue
Exploratory Data Analysis for the Evaluation of Tribological Properties of Wear-Resistant Surface Layers Modified with Rare-Earth Metals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of Heat Source Path on Thermal Evolution during Electro-Gas Welding of Thick Steel Plates

1
School of Engineering, University of Leicester, Leicester LE1 7RH, UK
2
Nanjing Iron & Steel United Co., Ltd., Nanjing 210035, China
3
School of Materials and Physics, China University of Mining and Technology, Xuzhou 221116, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(6), 2215; https://doi.org/10.3390/ma15062215
Submission received: 23 January 2022 / Revised: 3 March 2022 / Accepted: 15 March 2022 / Published: 17 March 2022
(This article belongs to the Special Issue Metallurgical Process Simulation and Optimization)

Abstract

:
In recent years, the shipbuilding industry has experienced a growing demand for tighter control and higher strength requirements in thick steel plate welding. Electro-gas welding (EGW) is a high heat input welding method, widely used to improve the welding efficiency of thick plates. Modelling the EGW process of thick steel plates has been challenging due to difficulties in accurately depicting the heat source path movement. An EGW experiment on 30 mm thickness E36 steel plates was conducted in this study. A semi-ellipsoid heat source model was implemented, and its movement was mathematically expressed using linear, sinusoidal, or oscillate-stop paths. The geometry of welding joints, process variables, and steel composition are taken from industrial scale experiments. The resulting thermal evolutions across all heat source-path approaches were verified against experimental observations. Practical industrial recommendations are provided and discussed in terms of the fusion quality for E36 steel plates with a heat input of 157 kJ/cm. It was found that the oscillate-stop heat path predicts thermal profile more accurately than the sinusoidal function and linear heat path for EGW welding of 30 mm thickness and above. The linear heat path approach is recommended for E36 steel plate thickness up to 20 mm, whereas maximum thickness up to 30 mm is appropriate for sinusoidal path, and maximum thickness up to 35 mm is appropriate for oscillate-stop path in EGW welding, assuming constant heat input.

1. Introduction

In the shipbuilding industry and offshore engineering, the size and structures of ships are increasing rapidly [1,2], thus giving rise to a growing demand for higher control in welding thick steel plates [3,4]. Electro-gas welding (EGW) has become an indispensable welding method for shipbuilding enterprises due to its outstanding advantages, such as large heat input and high welding efficiency [5]. As illustrated in Figure 1, EGW is an automatic welding process using a special flux-cored wire with CO2 gas protection, and used for the welding of vertical position of steel plates [6]. During the welding process, the torch moves along the chosen weld path from bottom to top. A water-cooled copper slider is placed on the front of the weld, and a ceramic “backing” plate is positioned at the rear. The welding torch is also used to feed the welding metal into the groove. In thick plates welding, the welding pool is restrained by the weld pieces, backing plate and sliding copper shoe, so that single pass welding can be accomplished.
In EGW, heat transfer from the heat source to the plates occurs mainly through radiation and convection [7,8]. In thick steel plate welding, the efficiency of heat transfer in the “y” direction (with depth, Figure 1) of single pass EGW is a key factor affecting the quality of the joint [9,10]. In 2016, Hwang, Kim and Lee [11] introduced double ellipsoidal moving heat sources to model the temperature profile and residual stress distribution in EGW of marine steel. However, the ellipsoidal method oversimplifies the oscillate-stop heat source path of industrial welding process to 1 dimensional linear case, whereby the movement path in the “y” direction is not considered. Therefore, the heat source paths are inconsistent with the real EGW process, hence the predicted fusion lines do not agree with experimental results.
To include an oscillate-stop heat source path within weld modelling, in 2017, Xu, Pan and Wan [12] used a sinusoidal function to simulate metal active gas arc welding (GAW). In 2018, Yuan et al. [13] applied a piecewise function to model the real oscillate movement path of gas metal arc welding. In 2019, da Silva Pereira et al. [14] implemented weave patterns by path parameterization, which improved the accuracy and has been further extended to predict defects in welding [15]. However, the effect of different heat source paths on thermal evolution in modelling EGW of thick steel plates and further optimize heat source paths according to thickness of steel plate has not yet been investigated systematically.
This study investigates the effect of three different heat source paths (linear, sinusoidal and oscillate-stop) through FEA on the resulting thermal evolution in thick steel plate welding. The results for EGW welding of E36 marine steel plates with thicknesses ranging from 20–40 mm and a heat input of 157 kJ/cm are compared based on fusion efficiency. The FEA models are verified by comparing against experimental observations. Finally, practical industrial recommendations are provided for each heat source path’s approach, according to variation in plate thickness.

2. Materials and Model Description

2.1. Electro-Gas Welding Experiments

This study investigates the effect of heat source path on thermal evolution during EGW of E36 marine steel plates. For the FEA modelling, experimental parameters, such as the geometry of the welding joints, thermo-mechanical and EGW process variables, and steel composition were considered as inputs, listed in Table 1. The temperature variation of material properties such as density and thermal conductivity of the E36 steel plate was also taken into account. JMatPro was implemented to calculate the material properties of the steel [16,17] for the E36 steel composition, shown in Table 2. In addition to the five main elements: carbon, manganese, phosphorus, and sulphur, microalloying elements such as niobium, titanium, and aluminium were added to improve the mechanical properties of steel. The dimensions of the steel plate model are shown in Figure 2. The thickness of the steel plate is 30 mm, while the length and width of the steel plate are 800 mm and 250 mm, respectively.
In the electro-gas welding process of 30 mm thick steel E36, flux-cored wire JIS Z3319 DW-S60G with a diameter of 1.6 mm was used. The chemical composition and the properties of the weld deposit are presented in Table 3 and Table 4.

2.2. Three-Dimensional EGW Model for the Numerical Simulation

Heat transfer during welding is a complex process dependent on multiple thermo-physical and geometrical factors. In EGW, heat from the heat source is dissipated to the welding plates mainly by conduction, radiation, and convection. Thus, a balance must be found between the rate of heat generation (speed and energy of EGW) and the rate of dissipation, determined by the case specific geometry and chemistry.
The heat transfer equation is given as follows:
ρ C p T t + ( k T ) = Q Q r Q c
where ρ   —density of the material; C p —specific heat capacity; T—absolute temperature; t—time; k—thermal conductivity; Q —heat source; Q r —radiation heat loss; Q c —convection heat loss. During the EGW process, the peak temperature of the heat-affected zone can reach more than 2000 °C, thus the temperature gradient between the steel plate and the environment is massive, causing significant radiation. The heat dissipation terms for radiation and convection are:
Q r = ε σ ( T a m b 4 T 4 )
Q c = h f ( T a m b T )
where ε —thermal emissivity; σ —Stephan Postman’s constant: 5.67 × 10−8 [W/(m2·K4)]; T a m b —ambient temperature [K]; T—steel plate surface temperature [K]; q r —radiant heat transfer flux [W/m2];   q c —convective heat transfer flux [W/m2]; h f —convection heat transfer coefficient between welding parts and the environment.
A three-dimensional EGW model was implemented, as shown in Figure 3a,b. This model uses a free tetrahedral mesh with the following principles designed for improving the computational efficiency:
  • The weld region close to the weld heat source has a larger temperature gradient and, hence, was divided into a finer grid (minimum mesh tetrahedral edge length 7.67 mm).
  • A coarser mesh (maximum mesh tetrahedral edge length 24.24 mm) was applied to the base steel plate areas far away from the heat source, where a smaller change in thermal gradients occurred.
Figure 3b illustrates the thermal profile and fusion line of a steel plate during the EGW process at 520 s at the midsection. The thermal distribution at surface and side of the steel plate is available from the diagram. Experimental fusion lines are used to verify the simulation results. A schematic diagram of the pre-welding bevel and micrograph of real joint after welding is shown in Figure 3c,d, respectively.
The filling of weld metal gradually with the movement of heat source is considered. The study simulates the ‘activate’ and ‘inactive’ of weld metal. Geometry of weld metal is pre-drawn in the work piece gap and properties of weld metal are activate point by point with the heat source movement.
For the welded area of the weld metal, the material properties of the geometry are activated.
For the unwelded area of the weld metal, the material properties of the geometry are inactivated, as shown in Figure 4.

2.3. Heat Source Movement Path

2.3.1. Heat Source Model

The earliest and simplest heat source model is the point heat source model proposed by Rosenthal [18] and it has been widely used in welding simulation. Rosenthal’s heat source model applies a quasi-steady state 3D semi-infinite geometry for point source. A velocity term was added to Rosenthal’s model to simulate heat source movement by Lecoanet et al. [19,20,21]. To describe the heat source distribution, a double ellipsoidal heat source model [22,23] was proposed by Rouquette et al., combining two different ellipses, one in the front quadrant and the other in the rear quadrant [24,25]. Laser Welding Processes have been simulated by a double ellipsoidal heat source model in recent years [26,27,28]. The two heat source models, i.e., Gaussian, semi-ellipsoid and double ellipsoidal heat source model, predict similar temperature distribution and distortion [29,30]. A simplified Gaussian heat source model was to improve calculation efficiency by Cai and Norman [31,32].
In this study, the semi-ellipsoid heat source model was selected. The equation of semi-ellipsoid heat source model is shown below:
q ( x ,   y ,   z ) = 6 Q π r 3 π exp ( 3 ( x 2 + y 2 + z 2 ) r 2 )
where q ( x , y , z ) represents the heat flow density distribution of ( x , y , z ) , Q is the effective power of the arc,   r is the radius of the semi-ellipsoid heating source. The schematic diagram of semi-ellipsoid heat source is shown in Figure 5.

2.3.2. Heat Source Movement Path

The movement of the heat source is expressed using a linear, sinusoidal and oscillate-stop path. In this study, the movement speed in the z -axis is defined as v 1 and in the y-axis as v 2 . As shown in Figure 6, three welding paths are modelled. Their exact coordinate heat source locations at each time step are expressed in terms of v 1 , v 2 and then substituted in Equation (4).
1.
Linear path heat source
The linear path heat source considers the welding heat source moving in the welding direction with a welding speed of v 1 ; the speed of the y -axis direction is zero. The expression of welding heat flux for its movement along the welding direction is shown in Equation (5):
q ( x , y , z ) = 6 Q π r 3 π exp ( 3 ( x 2 + y 2 + ( z v 1 t ) 2 ) r 2 )
The simulation of the welding heat source moves from bottom to top along the welding direction with the welding speed during the electro gas welding process. The heat source moves along the centre of the steel plates for different thicknesses, with a welding speed of 6.9 cm/min, as given in Table 1.
2.
Sinusoidal path heat source
The heat source, moving both in the welding direction with speed v 1 and in the thickness direction with speed v 2 , is considered in the sinusoidal path heat source model. The schematic diagram of the path is shown in Figure 6b.
The defined oscillate range is d , the length of the weld beam is L , the welding speed is v 1 , and the oscillate speed of the heat source in the thickness direction is v 2 ; the equation of the welding period is:
T = 2 d v 2
Total welding time:
t = L v 1
The heat source position of the y -axis direction is considered as:
y = d 4 s i n 2 π t T
Steel plates with three different thicknesses of 30 mm, 35 mm and 40 mm were used in this study to investigate the applicable thickness for different heat source movement paths. The parameters of the sinusoidal for different thickness steel plates are shown in Table 5.
According to the parameters of the sinusoidal heat source movement path, the position of the heat source changes with time during the welding process, as described in Equation (10). The amplitude of sinusoidal function is 10 mm, 12.5 mm and 15 mm for three thick plates, the period of the sinusoidal cycle is 5.5 s, respectively. The oscillate range is 20 mm, 25 mm and 30 mm. Due to the heat source needing to be close to the surface while leaving a certain distance for technological factors, the distance of 5 mm to the top surface and the bottom surface is appropriate.
P 1 ( t ) = A sin ( 2 π t T ) + k
The welding heat flux of the sinusoidal path is described as:
q ( x , y , z ) = 6 Q π r 3 π exp ( 3 ( x 2 + ( y P 1 ( t ) ) 2 + ( z v 1 t ) 2 ) r 2 )
According to Equation (10), the position of heat source change with time for three thickness plates are shows in Figure 7. The y -axis is the position along the thickness of steel plates and the z -axis is the position along the welding direction from the bottom to the top of steel plates. The movement path of the heat source under three thicknesses is drawn using solid lines with red, blue and green colour and the three thickness of 30 mm, 35 mm and 40 mm steel plates are indicated as dotted lines with red, blue and green colour, respectively. The positions of the heat sources for 1 s, 5.5 s, 6 s, 7 s and 11 s are marked in Figure 7.
3.
Oscillate-stop heat source
According to the characteristics of the oscillate-stop path, the heat source not only moves from bottom to top during the welding process but also oscillates in the direction of the depth of the melt pool, which helps the welding bevels on both sides to obtain the same melting depth. According to the welding practice, the oscillate-stop parameters of the heat source for different thickness steel plates are shown in Table 6.
Assume that the oscillate range is R, the distance of the oscillate centre to the origin of the coordinates is C, the period of the cycle is T, cycles is n. according to parameters of the oscillate-stop path; the position of heat source relative to the steel plate with different thickness in the welding period is described as Equation (12):
P 2 ( t ) =   { R 1.5 t + 10 5.5 ( n 1 ) R 1.5 ,                                   T ( n 1 ) < t T ( n 1 ) + 1.5 C + R 2 ,                                               T ( n 1 ) + 1.5 < t T ( n 1 ) + 2.5 R 1.5 t + 30 [ 5.5 ( n 1 ) + 2.5 ] ( R 1.5 ) ,   T ( n 1 ) + 2.5 < t T ( n 1 ) + 4 C R 2 ,                                               T ( n 1 ) + 4 < t T ( n 1 ) + 5.5
Regarding their thickness, the oscillate centre C = 20 mm and period T = 5.5 s are constant, so the equation can be simplified to:
P 2 ( t ) =   { R 1.5 ( t 5.5 n + 5.5 ) + 10 ,                         5.5 n 5.5 < t 5.5 n 4 20 + R 2 ,                                                           5.5 n 4 < t 5.5 n 3 R 1.5 ( t 5.5 n + 30 ) + 30 ,                     5.5 n 3 < t 5.5 n 1.5 20 R 2 ,                                                           5.5 n 1.5 < t 5.5 n
Define a function for the cycle of the equation, P ( m o d ( t , 5.5 ) ) means the equation cycle one time every 5.5 s. The welding heat flux of oscillate stop path be described as:
q ( x , y , z ) = 6 Q π r 3 π exp ( 3 ( x 2 + ( y P 2 ( m o d ( t , 5.5 ) ) 2 + ( z v 1 t ) 2 ) r 2 )
The position of the heat source changes with time for three thickness plates, as shown in Figure 8.

3. Model Verification via Experiment

To verify the accuracy of the model, a cross-section is taken along the thickness (y-axis) of the weld joint. As shown in Figure 9, the calculated melting pool and fusion line using different heat source path modes are compared against the experimental welding joints. The model fusion line is taken at the material melting point, 2055 K.
(1)
As shown in Figure 9b, the application of the linear heat source path did not cause sufficiently high temperatures to melt the top and bottom of the joint; only the centre part of joint is melted, which differs to the experimental observation.
(2)
The application of the sinusoidal path (see Figure 9c) leads to a fully melted weld joint in the thickness direction. The fusion line extended to the outside of the welding groove, but the area covered by the fusion line is smaller than that observed in the experiments (see Figure 9a).
(3)
The application of the oscillate-stop path leads to a fully melted weld joint in the thickness direction. Additionally, the fusion line extended to the outside of the welding groove base metal near the fusion is melted to form a solid joint, and the area covered by the fusion line is similar to that observed in the experiments.
Quantitative analysis is carried out to measure the coordinates of the points of intersection between the fusion lines and top and bottom surface (marked as A, B, C and D in Figure 9). The coordinates for points A, B, C, and D are listed in Table 7. The difference between simulated and experimental coordinates is calculated, defined as “error” as shown in Table 7, where Error = (simulated results-experiment result)/experiment result. The “×” sign refers to the lack of fusion at the given location, i.e., no fusion line in the area.
In the case of the linear heat source path, there are no fusion lines running through the top to the bottom of the welded joint. For the sinusoidal heat source path and the oscillate-stop heat source path, the welding arc not only moves from bottom to top during the welding process but also along with the weld oscillates in the direction of the depth of the melt pool, which helps the welding bevels on both sides to obtain the same melting depth. The “errors” of four selected points with sinusoidal heat source are: 21.4%, 15.5%, 28.6% and 25.0%, respectively, while the error of four selected points with oscillate-stop heat source is 8.8%, 2.0%, 4.8% and 0%, respectively. The “error” of the oscillate-stop heat source is greatly reduced. The heat source with oscillate-stop path holding on at near surface of steel plate for more times than sinusoidal path, which helps to transfer more heat to the surface of the joint during the welding process.
In summary, the path of the heat source can significantly affect the thermal profile of the weld joint, hence FEA models are essential in optimizing and predicting accurately the thermal profile, including the fusion line and shape of the melt zone in different thickness and path conditions.

4. Evolution of Thermal Profile in Heat Affected Zone

4.1. Simulated Heat Source Paths

Heat source paths calculated using the linear heat source path, sinusoidal path and oscillate-stop path are shown in Figure 10. In the figure, thermal profiles indicating heat sources at 340 s, 349 s, and 357 s are shown as examples. The linear heat source moves from the bottom to the top along the centre of the thickness of the steel plate, the sinusoidal heat source moves along with the weld direction and the thickness of the weld pool at the same time, and the oscillate-stop heat source path stays near the surface for a period of time when oscillating to the internal and external position. The numbers in brackets on the right side of every picture are the coordinates of the heat source centre of every picture. The x axis alone the welding line and the y axis alone is the plate thickness direction. When the welding time is 0 s, the heat source centre is in the starting point of the welding line (z) and the centre of plate thickness (y); the coordinates of the heat source centre are the origin coordinates (0, 0).

4.2. Thermal Cycles in Heat Affected Zone

Figure 11a is a schematic diagram indicating the positions of the selected points in the heat affected zone (HAZ). Point 1, point 6 and point 8 are located in the groove which belongs to the molten pool, point 2, point 3, point 4, point 5 and point 7 are located in the base metal, the distance to the heat source of the five points is from 11 mm to 50 mm. Figure 11b–d show the simulated thermal cycle curve of the five selected points during the welding process with the linear path, sinusoidal heat source path and oscillate-stop heat source path, respectively.
As illustrated in Figure 11, the temperature quickly rises when the welding heat source approaches the selected point, reaching a peak value, and then it decreases gradually as the welding heat source moves away from the point. The peak temperature varies according to the distance from the welding central line to the points. As listed in Table 8, the peak temperature of point 1 with a linear heat source path is 1932 K, which does not reach the melting point of 2055 K. The peak temperature of point 1 with sinusoidal heat source path is 2085 K, which exceeds the melting point by 30 K, the point 1 will be melted. The peak temperature of point 1 with the oscillate-stop heat source path is 2156 K, which exceeds the melting point by 106 K, compared with the sinusoidal heat source path; its fusion line advances to a position further away from the groove to form a joint with better quality. The peak temperature values of the selected point are listed in Table 8.
Using the sinusoidal path and the oscillate-stop path, the distance between the heat source and the selected points becomes closer when the heat source moves to the surface. While the oscillate-stop path stays near the internal and external position for a period of time in the cycle of movement. The welding arc can transfer more heat to the surface of steel plates and advance the fusion line in the position. So, the temperature in HAZ with the oscillate-stop path is higher than that with the sinusoidal path. EGW is a single pass method for welding thick steel plates, so the stop-over of the heat source at a near-surface position has great significance on the quality of the weld joint.

5. Recommended Heat Source Path Model for Simulating EGW Welding of Thick Steel Plates

Three different heat source paths have been used for simulating EGW welding of thick steel plates. In this section, the study will examine the effect of the heat source model on the shape of weld pool and the fusion line, so that a different applicable heat source model for simulating EGW thick plates can be defined.
The size of the cross section of the work piece gap is shown in Figure 12. The heat source movement path for 20 mm, 25, 30 mm, 35 mm and 40 mm thick steel plate were designed and listed in Table 5 and Table 6 and Figure 7 and Figure 8.
Simulations were carried out on steel plates with different thicknesses and different heat source models, as listed below. Other modelling parameters are listed in Table 3 and Table 4, and the heat input is 157 kJ/cm for all trials.
  • 20 mm, 25 mm thickness steel plates—using the linear path (Figure 12a,b);
  • 30 mm, 35 mm thickness steel plates—using sinusoidal path (Figure 12c,d);
  • 35 mm and 40 mm thick steel plates—using oscillate-stop path (Figure 12e,f).
For the linear path model: Figure 12a shows a fully fusion welding joint for 20 mm thick steel plates, while Figure 12b shows the lack of fusion in the bottom of the joint for 25 mm thick steel plates. So, the maximum applicable thickness for the linear heat source path is estimated to be 20 mm.
For the sinusoidal heat source path model: Fully fusion of the weld pool was obtained for 30 mm thick steel plates, as shown in Figure 12c, while Figure 12d shows a lack of fusion in the bottom of the joint for 35 mm steel plate. So, the maximum applicable thickness for the sinusoidal function path is estimated to be under 30 mm.
For the oscillate-stop path: Figure 12e shows that a fully fusion of the weld pool is obtained for 35 mm thick steel plates, but the lack of fusion in the bottom of the joint for 40 mm thick steel plates. So, the maximum applicable thickness for the oscillate-stop heat source path is estimated to be around 35 mm.

6. Conclusions

  • Three different types of heat source path models (linear, sinusoidal function and oscillate-stop) were implemented to simulate the EGW process of marine steel with a heat input of 157 kJ/cm.
  • For EGW welding of 30 mm thickness steel plates, the model using the oscillate heat source path predicted a more accurate thermal profile (the shape of melt pool and the fusion line) than those using the sinusoidal and linear heat source paths.
  • The applicable heat source paths for modelling the EGW process of steel plates with different thicknesses were investigated. The linear path model can be used for simulating steel plate thickness up to 20 mm, a maximum thickness of 30 mm is appropriate for the sinusoidal path, and a maximum thickness of 35 mm is appropriate for the oscillate-stop path in EGW welding, with a heat input of 157 kJ/cm.

Author Contributions

Conceptualization, J.F. and H.D.; methodology, Q.T.; software, J.F.; validation, J.F., B.N. and M.L.; formal analysis, Q.T.; investigation, J.F.; resources, H.D.; data curation, X.Y.; writing—original draft preparation, J.F.; writing—review and editing, B.N., J.F.; supervision, H.D.; project administration, B.T.; funding acquisition, B.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

For further data sets, please contact the corresponding author.

Acknowledgments

All authors gratefully acknowledge the support from Nanjing Iron & Steel United Co., Ltd. (NISCO), in particular for providing a PhD scholarship for this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ichimiya, K.; Sumi, H.; Hirai, T. 460 MPa Yield Strength Class Steel Plate with JFE EWEL® Technology for Large Heat Input Welding. JFE Tech. Rep. 2008, 5, 7–12. [Google Scholar]
  2. Kaplan, M.B.; Solomon, S. A coming boom in commercial shipping? The potential for rapid growth of noise from commercial ships by 2030. Mar. Policy 2016, 73, 119–121. [Google Scholar] [CrossRef]
  3. Sasaki, K.; Suda, K.; Motomatsu, R.I.; Hashiba, Y.; Ohkita, S.; Imai, S. Development of Two-electrode Electro-gas Arc Welding Process. Shinnittetsu Giho 2004, 90, 57–64. [Google Scholar]
  4. Fritzsche, A.; Avilov, V.; Gumenyuk, A.; Hilgenberg, K.; Rethmeier, M. High power laser beam welding of thick-walled ferromagnetic steels with electromagnetic weld pool support. Phys. Procedia 2016, 83, 362–372. [Google Scholar] [CrossRef] [Green Version]
  5. Wang, G.; Li, J.B.; Sha, Y.Z.; Ma, J. Study on the penetration of electrogas welding. Weld. Technol. 2002, 5, 8–9. [Google Scholar]
  6. Eakkachai, W.; Koei, H.; Motomichi, Y.; Kenji, S.; Kota, K.; Tadakazu, T.; Hiroshi, Y.; Tsutomu, F.; Shin, N.; Tetsuro, N.; et al. Welding phenomena during vertical welding on thick steel plate using hot-wire laser welding method. Proc. Weld. Soc. 2015, 33, 143s–147s. [Google Scholar]
  7. Haelsig, A.; Mayr, P. Energy balance study of gas-shielded arc welding processes. Weld. World 2013, 57, 727–734. [Google Scholar] [CrossRef]
  8. Bai, X.; Colegrove, P.; Ding, J.; Zhou, X.; Diao, C.; Bridgeman, P.; Hönnige, J.R.; Zhang, H.; Williams, S. Numerical analysis of heat transfer and fluid flow in multilayer deposition of PAW-based wire and arc additive manufacturing. Int. J. Heat Mass Transf. 2018, 124, 504–516. [Google Scholar] [CrossRef] [Green Version]
  9. Aucott, L.; Huang, D.; Dong, H.B.; Wen, S.W.; Marsden, J.A.; Rack, A.; Cocks, A.C.F. Initiation and growth kinetics of solidification cracking during welding of steel. Sci. Rep. 2017, 7, 40255. [Google Scholar] [CrossRef] [Green Version]
  10. Aucott, L.; Huang, D.; Dong, H.B.; Wen, S.W.; Marsden, J.; Rack, A.; Cocks, A.C.F. A Three-stage mechanistic model for solidification cracking during welding of steel. Met. Mater. Trans. A 2018, 49, 1674–1682. [Google Scholar] [CrossRef] [Green Version]
  11. Hwang, S.-Y.; Kim, Y.; Lee, J.-H. Finite element analysis of residual stress distribution in a thick plate joined using two-pole tandem electro-gas welding. J. Mater. Process. Technol. 2016, 229, 349–360. [Google Scholar] [CrossRef]
  12. Xu, G.; Pan, H.; Wang, J. Numerical analysis model of temperature profile in swing-arc narrow gap GMAW welding. Thansactions China Weld. Inst. 2017, 38, 55–60. (In Chinese) [Google Scholar]
  13. Yuan, S.; Liu, W.; Li, L.; Jiang, X. Heat source model of temperature profile in swing-arc narrow gap MAG welding considering sidewall fusion. Thansactions China Weld. Inst. 2018, 39, 95–99. (In Chinese) [Google Scholar]
  14. Pereira, H.A.D.S.; Rodrigues, M.C.; Firmino, J.V.L.D.C. Implementation of weave patterns by path parameterization in the simulation of welding processes by the finite element method. Int. J. Adv. Manuf. Technol. 2019, 104, 477–487. [Google Scholar] [CrossRef]
  15. Wang, Z.; Oliveira, J.P.; Zeng, Z.; Bu, X.; Peng, B.; Shao, X. Laser beam oscillating welding of 5A06 aluminum alloys: Microstructure, porosity and mechanical properties. Opt. Laser Technol. 2019, 111, 58–65. [Google Scholar] [CrossRef]
  16. Guo, Z.; Saunders, N.; Schillé, J.; Miodownik, A. Material properties for process simulation. Mater. Sci. Eng. A 2009, 499, 7–13. [Google Scholar] [CrossRef]
  17. Saunders, N.; Guo, U.K.Z.; Li, X.; Miodownik, A.P.; Schillé, J.-P. Using JMatPro to model materials properties and behavior. JOM 2003, 55, 60–65. [Google Scholar] [CrossRef]
  18. Rosenthal, D. The theory of moving sources of heat and its application of metal treatments. Trans. ASME 1946, 68, 849–866. [Google Scholar]
  19. Lecoanet, A.; Ivey, D.; Henein, H. Simulation of the temperature profile during welding with COMSOL multiphysics® software using Rosenthal’s approach. In Proceedings of the 2014 COMSOL, Cambridge, UK, 17 September 2014. [Google Scholar]
  20. Wang, H.X.; Sun, J.S.; Wei, Y.H.; Zheng, Y.Y. Simulation of GMAW thermal process based on string heat source model. Sci. Technol. Weld. Join. 2005, 10, 511–520. [Google Scholar] [CrossRef]
  21. Eagar, T.; Tsai, N. Temperature fields produced by traveling distributed heat sources. Weld. J. 1983, 62, 346–355. [Google Scholar]
  22. Rouquette, S.; Guo, J.; Le Masson, P. Estimation of the parameters of a Gaussian heat source by the Levenberg–Marquardt method: Application to the electron beam welding. Int. J. Therm. Sci. 2007, 46, 128–138. [Google Scholar] [CrossRef]
  23. Fu, G.; Gu, J.; Lourenco, M.I.; Duan, M.; Estefen, S.F. Parameter determination of double-ellipsoidal heat source model and its application in the multi-pass welding process. Ships Offshore Struct. 2014, 10, 204–217. [Google Scholar] [CrossRef]
  24. Goldak, J.; Chakravarti, A.; Bibby, M. A new finite element model for welding heat sources. Metall. Trans. B 1984, 15, 299–305. [Google Scholar] [CrossRef]
  25. Podder, D.; Mandal, N.; Das, S. Heat source modeling and analysis of submerged arc welding. Weld. J. 2014, 93, 183–191. [Google Scholar]
  26. Kik, T. Computational techniques in numerical simulations of arc and laser welding processes. Materials 2020, 13, 608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Kik, T.; Górka, J. Numerical simulations of laser and hybrid S700MC T-joint welding. Materials 2019, 12, 516. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Kik, T. Heat source models in numerical simulations of laser welding. Materials 2020, 13, 2653. [Google Scholar] [CrossRef] [PubMed]
  29. Hashemzadeh, M.; Chen, B. Comparison between different heat sources types in thin-plate welding simulation. In Developments in Maritime Transportation and Exploitation of Sea Resources; Taylor & Francis: Abingdon, UK, 2013; pp. 329–335. [Google Scholar]
  30. Tong, M.; Duggan, G.; Liu, J.; Xie, Y.; Dodge, M.; Aucott, L.; Dong, H.; Davidchack, R.L.; Dantzig, J.; Barrera, O.; et al. Multiscale, multiphysics numerical modeling of fusion welding with experimental characterization and validation. JOM 2012, 65, 99–106. [Google Scholar] [CrossRef] [Green Version]
  31. Cai, Z.; Wu, S.; Lu, A.; Zhao, H.; Shi, Q. Line Gauss heat source model: An efficient approach for numerical welding simulation. Sci. Technol. Weld. Join. 2001, 6, 84–88. [Google Scholar] [CrossRef]
  32. Norman, A.F.; Ducharme, R.; Mackwood, A.; Kapadia, P.; Prangnell, P.B. Application of thermal modelling to laser beam welding of aluminium alloys. Sci. Technol. Weld. Join. 1998, 3, 260–266. [Google Scholar] [CrossRef]
Figure 1. The electro-gas welding (EGW). (a) The 3D schematic diagram of EGW equipment, (b) The 2D schematic diagram of EGW equipment, (c) Photographs of EGW equipment, (d) Photographs of experimentation.
Figure 1. The electro-gas welding (EGW). (a) The 3D schematic diagram of EGW equipment, (b) The 2D schematic diagram of EGW equipment, (c) Photographs of EGW equipment, (d) Photographs of experimentation.
Materials 15 02215 g001
Figure 2. Geometry of welded steel plate.
Figure 2. Geometry of welded steel plate.
Materials 15 02215 g002
Figure 3. (a) Geometry and mesh of EGW model implemented in this study, (b) the transient temperature distribution simulated of steel plate at 520 s, (c) schematic diagram welding joint of pre-welding bevel, (d) welding joint of after-welding showing the fusion line.
Figure 3. (a) Geometry and mesh of EGW model implemented in this study, (b) the transient temperature distribution simulated of steel plate at 520 s, (c) schematic diagram welding joint of pre-welding bevel, (d) welding joint of after-welding showing the fusion line.
Materials 15 02215 g003
Figure 4. Schematic diagram of the activation process of the weld metal.
Figure 4. Schematic diagram of the activation process of the weld metal.
Materials 15 02215 g004
Figure 5. Schematic diagram of heat source (semi-ellipsoid heat source).
Figure 5. Schematic diagram of heat source (semi-ellipsoid heat source).
Materials 15 02215 g005
Figure 6. Schematic diagram of heat source path. (a) linear heat source path, (b) sinusoidal path heat source, (c) oscillate-stop heat source path.
Figure 6. Schematic diagram of heat source path. (a) linear heat source path, (b) sinusoidal path heat source, (c) oscillate-stop heat source path.
Materials 15 02215 g006
Figure 7. Sinusoidal heat source movement vs. time (t) for 30 mm, 35 mm and 40 mm thick steel plates.
Figure 7. Sinusoidal heat source movement vs. time (t) for 30 mm, 35 mm and 40 mm thick steel plates.
Materials 15 02215 g007
Figure 8. Movement path of oscillate-stop for thickness of 30 mm, 35 mm and 40 mm.
Figure 8. Movement path of oscillate-stop for thickness of 30 mm, 35 mm and 40 mm.
Materials 15 02215 g008
Figure 9. Melt pool and fusion line of experiments and simulated welded joints of 30 mm thickness steel plates, (a) experimental weld joint, (b) modelling using linear path heat source, (c) modelling using sinusoidal path heat source, (d) modelling using oscillate-stop path heat source.
Figure 9. Melt pool and fusion line of experiments and simulated welded joints of 30 mm thickness steel plates, (a) experimental weld joint, (b) modelling using linear path heat source, (c) modelling using sinusoidal path heat source, (d) modelling using oscillate-stop path heat source.
Materials 15 02215 g009
Figure 10. Calculated thermal profiles using three different heat source paths at different welding times of 340 s, 349 s and 357 s. (a) linear path (b) sinusoidal path, (c) oscillate-stop path.
Figure 10. Calculated thermal profiles using three different heat source paths at different welding times of 340 s, 349 s and 357 s. (a) linear path (b) sinusoidal path, (c) oscillate-stop path.
Materials 15 02215 g010
Figure 11. Calculated thermal cycles at selected points in heat affect zones (HAZ), (a) position of selected points in HAZ, (b) thermal cycle curves using linear path, (c) thermal cycle curves using sinusoidal path, (d) thermal cycle curves using oscillate-stop path.
Figure 11. Calculated thermal cycles at selected points in heat affect zones (HAZ), (a) position of selected points in HAZ, (b) thermal cycle curves using linear path, (c) thermal cycle curves using sinusoidal path, (d) thermal cycle curves using oscillate-stop path.
Materials 15 02215 g011
Figure 12. The simulated welding molten pool with three heat sauce paths. Linear heat source path model for 20 mm thickness (a), for 25 mm thickness (b); Sinusoidal path for 30 mm thickness (c), for 35 mm thickness (d); Oscillate-stop path for 35 mm thickness (e) and for 40 mm thickness (f).
Figure 12. The simulated welding molten pool with three heat sauce paths. Linear heat source path model for 20 mm thickness (a), for 25 mm thickness (b); Sinusoidal path for 30 mm thickness (c), for 35 mm thickness (d); Oscillate-stop path for 35 mm thickness (e) and for 40 mm thickness (f).
Materials 15 02215 g012
Table 1. Parameters of welding experiments.
Table 1. Parameters of welding experiments.
Heat Input
(kJ/cm)
Welding Speed (cm/min)Current
(A)
Voltage
(V)
Wire Diameter
(mm)
Wire Feed Rate
(m/min)
1576.9420431.613
Parameters of the torch movement of horizontal oscillations
Internal stop
(s)
External stop
(s)
Amplitude
(mm)
Frequency
(s−1)
11.5102/11
Table 2. Chemical composition of E36 steel (wt %).
Table 2. Chemical composition of E36 steel (wt %).
CSiMnSPNbTiAl
0.080.271.450.0020.010.0140.0140.034
Table 3. Chemical composition of weld deposit.
Table 3. Chemical composition of weld deposit.
GradeCMnSiSPCrNiMoBTi
DW-S60G0.071.680.330.0060.0110.020.770.26/0.02
Table 4. Mechanical properties of weld deposit.
Table 4. Mechanical properties of weld deposit.
Yield Strength Rp0.2
(MPa)
Tensile Strength Rm
(MPa)
Elongation A
(%)
Akv at −20 °C
(J)
53466226124, 139, 120
Table 5. Sinusoidal path parameters of welding heat source.
Table 5. Sinusoidal path parameters of welding heat source.
No.Thickness of Plate
(mm)
Amplitude (A)
(mm)
Period (T)
(s)
Oscillate Range
(mm)
Offset (k)
(mm)
13010.05.52020
23512.55.52520
34015.05.53020
Table 6. Oscillate-stop parameters of welding heat source.
Table 6. Oscillate-stop parameters of welding heat source.
No.Thickness
(mm)
Internal Stop
(s)
External Stop
(s)
Oscillate Range (R)
(mm)
Oscillate Center (C)
(mm)
Period (T)
(s)
1301.01.520205.5
2351.01.525205.5
3401.01.530205.5
Table 7. Comparison of coordinates of selected points between simulated and experimental welding melt pools.
Table 7. Comparison of coordinates of selected points between simulated and experimental welding melt pools.
CoordinateExperiment
(mm)
Linear
(mm)
ErrorSinusoidal
(mm)
ErrorOscillate-Stop
(mm)
Error
A(−15.9, 30)××(−12.5, 30)21.4%(−14.5, 30)8.8%
B(14.8, 30)××(12.5, 30)15.5%(14.5, 30)2.0%
C(6.3, 0)××(4.5, 0)28.6%(6.0, 0)4.8%
D(−6.0, 0)××(−4.5, 0)25.0%(−6.0, 0)0
Note: Error = (simulated result—experiment result)/experiment result.
Table 8. The peak temperature of selected point as shown in Figure 8.
Table 8. The peak temperature of selected point as shown in Figure 8.
Point1 (K)2 (K)3 (K)4 (K)5 (K)6 (K)7 (K)8 (K)
Path
Linear path19321673151314101065181016402420
Sinusoidal path20851848167015511153202118112295
Oscillate-stop path21561919173516091189241019602223
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fu, J.; Tao, Q.; Yang, X.; Nenchev, B.; Li, M.; Tao, B.; Dong, H. The Effect of Heat Source Path on Thermal Evolution during Electro-Gas Welding of Thick Steel Plates. Materials 2022, 15, 2215. https://doi.org/10.3390/ma15062215

AMA Style

Fu J, Tao Q, Yang X, Nenchev B, Li M, Tao B, Dong H. The Effect of Heat Source Path on Thermal Evolution during Electro-Gas Welding of Thick Steel Plates. Materials. 2022; 15(6):2215. https://doi.org/10.3390/ma15062215

Chicago/Turabian Style

Fu, Jun, Qing Tao, Xiaoan Yang, Bogdan Nenchev, Ming Li, Biao Tao, and Hongbiao Dong. 2022. "The Effect of Heat Source Path on Thermal Evolution during Electro-Gas Welding of Thick Steel Plates" Materials 15, no. 6: 2215. https://doi.org/10.3390/ma15062215

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop