2.1. Background
In tube welding, understanding and manipulating the parameters is critical to achieving a precise and strong joint. According to Table 1, our welding process was performed at a moderate speed of 0.003 m/s to ensure an accurate joint between the pipe sections. In addition, there is an electric current of 100 A, combined with a voltage of 12.5 V. An efficiency of 0.8 indicates the efficiency of energy conversion during the welding process and indicates a well-optimized system.
Table 1
Parameter | Speed | Current | Voltage | Efficiency |
Value | 0.003 m/s | 100 A | 12.5 V | 0.8 |
During the welding process, a non-uniform thermal field was created surrounding the weld, which might cause residual stress to develop in welded structures. To determine the thermal field, transient thermal analysis with the required boundary conditions was used. The 3D transient state heat conduction Equation, which is provided below, governs the transient thermal analysis of the welding process.
$$\frac{\partial }{{\partial x}}\left( {{K_x}\frac{{\partial T}}{{\partial x}}} \right)+\frac{\partial }{{\partial y}}\left( {{K_y}\frac{{\partial T}}{{\partial y}}} \right)+\frac{\partial }{{\partial z}}\left( {Kz\frac{{\partial T}}{{\partial z}}} \right)+Q=\rho {C_P}\left( {\frac{{\partial T}}{{\partial t}} - v\frac{{\partial T}}{{\partial x}}} \right)$$
1
Where, \({C_P}\) = Specific heat, \(J.k{g^{ - 1}}.K\); \(\rho =\) Mass density, \(kg.m{m^{ - 3}}\); \(Q=\) Internal heat generation, \(W.m{m^{ - 3}}\); \(v=\) Relative speed of heat source \(mm.{s^{ - 1}}\). The rate of heat rejection resulting from convection and radiation has been taken into account when solving the transient thermal analysis. Eq. (2) defines the convection and radiation boundary conditions.
$${K_x}\frac{{\partial T}}{{\partial x}}{n_x}+{K_y}\frac{{\partial T}}{{\partial y}}{n_y}+{K_z}\frac{{\partial T}}{{\partial z}}{n_z}+{q_r}+{h_c}\left( {T - {T_\infty }} \right)+\sigma \varepsilon F\left( {{T^4} - T_{r}^{4}} \right)=0$$
2
Where, \({h_c}=\)heat transfer coefficient (convection), \(W.{m^{ - 2}}.K\); \(k=\) Thermal conductivity, \(W.{m^{ - 1}}.K\); \(q=\) heatflux, W.\({m^{ - 1}}\); \(\sigma =\)Stefan-Boltzmann constant, \(5.68 \times {10^{ - 8}}~W.{m^{ - 2}}.{K^4}\); \(F=\) configuration factor; and \({n_x}\), \({n_y}\)and \({n_z}\)are the direction cosines of the boundary.
The Goldak warm temperature delivery model is important for proper welding simulations, in particular, withinside the context of the finite element Method (FEM). Unlike traditional models, Goldak's double ellipsoid design closely resembles actual heat distribution during welding (see Fig. 1). This three-dimensional model can be adapted to different processes and helps to predict metallurgical transformations and optimize weld quality. Its precision ensures structural integrity and improves the overall safety and efficiency of welding. Looking at the dimensions of the heat source in the figure, we see a front length of 5 mm and a rear length that is 15 mm longer, with a uniform width of 5 mm. The marked 3 mm depth demonstrates the penetrating ability of our heat source and strikes a balance between surface fusion and depth integrity. Taken together, these parameters ensure an efficient and precise welding process that promotes high-quality tubular joints that are essential for the intended applications. The equations below state the double ellipsoid heat source model, which best describes the heat source for arc welding:
$$Q\left( {x,~y,~z,~t} \right)=\frac{{6\sqrt 3 {f_f}{Q_w}}}{{{a_f}bc\pi \sqrt \pi }}{e^{ - 3{x^{'2}}/{a_f}^{2}}}{e^{ - 3{y^{'2}}/{b^2}}}{e^{ - 3{z^{'2}}/{c^2}}}$$
3
$$Q\left( {x,~y,~z,~t} \right)=\frac{{6\sqrt 3 {f_r}{Q_w}}}{{{a_r}bc\pi \sqrt \pi }}{e^{ - 3{x^{'2}}/{a_r}^{2}}}{e^{ - 3{y^{'2}}/{b^2}}}{e^{ - 3{z^{'2}}/{c^2}}}$$
4
Where x, y, and z are the welded pipe's local coordinates, \({f_f}\) and \({f_r}\) are parameters that indicate the proportion of heat deposited in the front and rear portions, respectively. Note that \({f_f}\) + \({f_r}\) = 2.0.
The welding heat source's power is \({Q_w}\). Calculations can be made based on the welding current, arc voltage, and arc efficiency. It is assumed that the TIG welding technique has an arc efficiency η of 80%. Parameters that can affect the characteristics of the welding heat source are \({a_f}\), \({a_r}\), b, and c. Under the welding requirements, the heat source's specifications can be changed to produce the necessary melted zone.
As illustrated in Fig. 2, our simulation incorporated one cylinder with a wallthickness of 3 mm, a length of 300 mm, and an outer diameter of 306 mm.
Incorporated within Fig. 3 are the boundary conditions for welding simulations that inherently encompass both thermal and mechanical factors. The thermal component of the solution encompasses free convection and radiation emanating from the material. To predict deformation and stress with precision, mechanical boundary conditions are a necessary prerequisite.
In the context of the static mechanical analysis, thermal loads are modeled as equivalent to body force. To determine structural analysis inputs, a specific course of action derived from the previous transient thermal study was utilized. Subsequently, the resulting stresses and strains were computed based on incremental stress-strain relationships, as described in Eq. (5).
$$d\varepsilon =d{\varepsilon ^l}+d{\varepsilon ^p}+d{\varepsilon ^T}$$
5
Where, \(d\varepsilon =\) total strain; \(d{\varepsilon ^l},d{\varepsilon ^p}{\text{~}}and{\text{~}}d{\varepsilon ^T}\) are the elastic strain, plastic strain, and thermal stress, respectively. In the field of welding, the selection of steel employed can exert a considerable influence on both the welding outcome and the ultimate longevity of the manufactured item. Within the realm of our simulation, we have opted to employ stainless steel (AISI 304) in the context of welding applications, owing to its distinctive attributes and attendant advantages. The AISI 304 alloy is a widely employed stainless steel alloy that encompasses both chromium and nickel, granting it substantial resistance towards rust, corrosion, and staining. Such a combination of traits renders it an exceptional choice for welding applications that necessitate elevated strength, durability, and resilience against harsh environments. Furthermore, the low carbon content of AISI 304 enhances its resistance to sensitization and ensures its exceptional weldability. The material properties of AISI 304 were sourced from [9], and Fig. 4 and Table 2 provide the material properties and their temperature dependence.
Table 2
Material Properties and Their Dependence on Temperature [9].
Temperature [°C] | Thermal conductivity [N/(s. K)] | Specific heatcapacity\(\left[ {m{m^2}/\left( {{s^2}.^\circ C} \right)} \right]\) | Expansion coefficient [1/°C] | Young’s modulus [GPa] |
20 | 15.7 | \(5.1 \times {10^8}\) | \(1.6 \times {10^{ - 5}}\) | 200 |
100 | 16.8 | \(5.25 \times {10^8}\) | \(1.7 \times {10^{ - 5}}\) | 195 |
200 | 17 | \(5.41 \times {10^8}\) | \(1.8 \times {10^{ - 5}}\) | 190 |
400 | 21 | \(5.72 \times {10^8}\) | \(1.9 \times {10^{ - 5}}\) | 180 |
600 | 23.5 | \(6.04 \times {10^8}\) | \(2 \times {10^{ - 5}}\) | 150 |
800 | 26.5 | \(6.3 \times {10^8}\) | \(2.05 \times {10^{ - 5}}\) | 128 |
1000 | 29.2 | \(6.48 \times {10^8}\) | \(2.12 \times {10^{ - 5}}\) | 70 |
1200 | 32.2 | \(6.73 \times {10^8}\) | \(2.2 \times {10^{ - 5}}\) | 15 |
1400 | 35.1 | \(6.91 \times {10^8}\) | \(7 \times {10^{ - 6}}\) | 10 |
1500 | 36.2 | \(7 \times {10^8}\) | \(0.1 \times {10^{ - 6}}\) | 2 |