1. Introduction
Mild steel is a type of low-carbon steel that is malleable and ductile, making it ideal for applications that require a low strength-to-weight ratio, such as cages, frames, and fencing. According to Silva et al. (2018), mild steel is preferred over other steels because it endures high stresses and results in longer service life [
1]. Due to its malleability, it can be quickly shaped, drilled, welded, and cut. Welding is the fabrication and repair of metal products and an essential process in every industry. It is a simple, cost-effective, and dependable method for joining metals. In welding processes, weld joint quality is usually affected by the weld material’s attributes, mechanical properties, and the highly concentrated localized zone. Saha et al. (2017) stated that the weldability of a material depends upon the properties changing during the welding process and the strength of the centralized weld zone [
2]. Datta et al. (2008) noted that the weld functionality features also depend on weld joint geometry, which is directly influenced by the input parameters [
3].
Several other welding processes have been proposed to join the specific types and dimensions of work parts: laser welding, TIG welding, electron beam welding, etc. [
4]. Cary et al. noted that TIG welding is extensively used in industries because of high weld quality, protection of the weld pool by using an inert gas shield, its versatility, repairability, and flexibility to adapt. Research showed that an optimal welding speed could produce a joint with high mechanical strength, excellent microstructure, and good morphology [
5]. Eshwar et al. (2014) found that the input parameter has a substantial role in determining the weld quality, which primarily includes bead width and height, and penetrating the area, thermal, and residual stresses, resulting in minimizing the life of a weld joint [
6]. Researchers have tried to work on different welding processes regarding weld bead geometry, temperature, and stress fields. For instance, Dhas et al. (2007) have proposed an adaptive neuro-fuzzy inference system (ANFIS) for determining the optimum input parameters for bead width [
7]. A multi-objective optimization of process parameters was carried out to find the optimal bead geometry using the Taguchi method [
8]. Fernando et al. (2021) investigated the impact of metal transfer modes on the symmetry of bead geometry. They measured bead geometry using a 3D scanner and found that the metal transfer mode affected the bead geometry symmetry. It provides insight into the effect of metal transfer modes on the quality of WAAM aluminum parts [
9]. Improvements in weld quality have been made by combining grey relational analysis with the Taguchi method. Chandrasekhar et al. (2015) have developed an artificial neural network (ANN) to estimate the weld bead width and depth of penetration by using IR images during TIG welding of 6 mm thick 316 LN stainless steel by changing the values of the welding current [
10]. Chokkalingham et al. (2012) have also developed a similar model for better estimation of bead width and depth of penetration [
11]. Kolahan et al. (2010) established a simulated annealing methodology to analyze and improve the geometry of the weld joint (bead width, bead height, and penetration) in gas metal arc welding [
12]. Zhang et al. (2017) have developed an online system for the quality of aluminum alloy welding, which automatically evaluates the quality parameters in real-time using online arc sound, voltage, and spectrum signals [
13]. The microstructure, hardness, and corrosion behavior of gas tungsten arc welding (GTAW) Inconel 625 super alloy over A517 carbon steel was also investigated using ERNiCrMo3 filler metal [
14]. Bandhu et al. (2022) have experimentally studied the low alloy steel after heat treatment [
15]. It was concluded that regulated metal deposition weldments have better mechanical properties than gas welding for the same steel alloy. Jiang et al. (2017) have addressed the quality issues (defects) in the welding process [
16]. They proposed a novel approach for classifying weld defects, such as porosity, slag inclusion, lack of penetration, etc., based on the Dumpster–Shafer evidence theory.
In recent years, the main problem faced by researchers is optimizing the factors that affect the welding joint [
17] and the geometry of the welding area to predict the temperature distribution, minimum residual stresses, and distortion [
18]. Understanding the thermal flow in welding is theoretically and experimentally important. The impact of stress relief heat treatment on the microstructure and mechanical properties of dissimilar GTAW weld joints between Inconel 625 and A106 carbon steel was also investigated [
19]. The results provide insights into the effects of heat treatment on the properties of such weld joints. Sepe et al. (2014) demonstrated the residual stresses and temperature distribution in metal arc welding using FEM analysis, finding higher thermal stress in heat-affected and fusion zones [
20]. The stresses decrease with the increase in temperature and pre-heating and post-heating methods. Katherasan et al. (2014) performed a simulation study for bead geometry and demonstrated the effects of process parameters such as feed rate, voltage, speed, and torch angle on bead width, reinforcement, and penetration depth using particle swarm optimization (PSO) for flux-cored arc welding [
21]. The effect of different process parameters such as voltage, welding current, and wire feed rate on the bead geometry and microstructure of WAAM parts were studied and the results showed that the process parameters significantly impacted the bead geometry and microstructure. It provides insights into optimizing WAAM process parameters for improved part quality and can help guide the WAAM process design [
22]. El-Sayed et al. (2017) also predicted the residual stresses and temperature-induced in friction stir welding using finite element analysis [
23].
A new calibration coefficient was computed for aluminum plates subjected to uniaxial load. An appropriate coefficient and residual stress compensation value could result in in-situ thermal stresses with expected accuracy within acceptable ranges of the industry’s specification limits [
24]. Kulkarni et al. (2019) demonstrated TIG welding joining dissimilar metals for stainless steel (P91 steel and AISI 3161) with Inconel 800 and Inconel 600 [
25]. It was observed that an Inconel 600 interlayer and a fully austenitic structure could be achieved without any mechanical loss. Based on the above literature, it has been realized that very little or no work has been reported to study the thermal and structural analysis and variations in temperature with the TIG welding parameters for the optimized process.
This research primarily focuses on analyzing multi-objective optimization of the TIG welding process for an optimum factor selection affecting bead geometry, temperature, and thermal stresses. Taguchi-based grey relational analysis for multiple objectives was employed to determine significant factors affecting the performance measures. The temperature variations over time and thermal stresses for the optimum geometry were studied. ANSYS APDL was used to generate the model, and finite element analysis (FEA) was employed to predict the responses. Heat flux was also determined by conducting the experiments with the selected levels of factors, which was further used in FEA in the V-shaped butt joint. As per the required dimensions, a 3D solid element with DOF, temperature (SOLID-70), and a 2D solid element with four nodes (PLANE-55) was used. A transient analysis followed by the coupled field analysis was performed to predict temperature variations and thermal stresses. The variation of temperature by welding parameters was also investigated.
2. Materials and Methods
TIG welding was achieved on the mild steel (MS) plates, and its chemical composition is presented in
Table 1 [
26]. A CNC milling machine was used to prepare test samples of dimensions 80 × 20 × 3 mm
3 for the V-shaped butt joint. The samples were appropriately clamped at the ends during welding. In experiments, argon gas was used because of its inert behavior, protecting the molten metal from impurities in the atmosphere and being heavier than air. The joining of metals was performed by melting the work partly with an arc fixed between the electrode (tungsten) diameter of 3.2 mm and the base part. A high-frequency spark delivered a conductive pathway through the air.
Due to the spark produced, the surface of the base work partly melted, and a molten pool was formed.
Figure 1 shows a schematic diagram and setup of TIG welding [
27].
Based on an extensive literature review, the selection of factors and their levels have been provided in
Table 2. The bead geometry was optimized for further analysis of temperature and thermal stresses, as shown in
Figure 2.
The bead geometry was divided into bead height, width, and penetration, which were measured using linear measuring instruments.
3. Experimental Design and Analysis
Taguchi design is an efficient, straightforward approach and saves more time by minimizing the number of experiments, as stated by Choudhury et al. [
28]. Hence it is used in the present study.
Table 3 shows the results of output responses concerning different sets of input factors. Nine experiments were performed by Taguchi orthogonal Array (L
93
3(OA)). A grey relational analysis (GRA) was employed to find the optimal setting for a final geometrical model for further investigation by FEA. GRA is a multi-objective optimization technique in which various responses can be optimized simultaneously. Three input factors, i.e., rate of flow (F), welding current (I), and gap distance (G), were selected to study their influence on bead height and width (H and W) and area of penetration (P) of the weld. Grey’s relational coefficient was determined to find the correspondence in actual and theoretical results.
Two criteria correspond to the output, i.e., the smaller- and larger-the-better cases. In the present study, weld geometry comprises bead height (H) and bead width (W) that corresponds to the lower case that can be represented by Equation (1):
And similarly, the response variable, penetration (P), follows the larger-the-better case, as shown in Equation (2):
In the above expressions (Equations (1) and (2)),
represents the grey relation generated value, where
and
shows the lowest and highest values for all kth responses, respectively, whereas
represents the initial (ideal) sequence for output responses (k = 1–9). The final grey grades show the DOF between sequences
and
. The grey relational coefficient,
was calculated using Equation (3):
In Equation (3),
shows the deviation sequence of the reference,
and
i.e., the ideal and comparability sequences. The symbol
represents the identification coefficient, and its value ranges from 0–1. Usually, it takes a value of 0.5 to ensure equal importance of all factors. The deviation sequence
is determined using Equation (4):
whereas in Equation (5),
show the final grey grades and response variables for corresponding experiments, respectively. Higher values of grey grades show better agreement between the ideal and given reference, i.e.,
The reference sequence shows the better sequence; thus, the high value of grade corresponds to combinations of input parameters closer to the optimum for each experiment using L
93
3(OA). Plots of the main effect and grand mean of grey relational grades are essential for efficient analysis.
4. Grey Relational Analysis
Responses are measured after performing experiments according to the design plan, as provided in
Table 3. After experimentation and data collection, the sequences
are calculated using Equations (2) and (3) by keeping both cases in mind, i.e., smaller and larger-the-better cases, respectively. Similarly, the deviation sequence
was determined using Equation (4) and,
and
are the higher and lower values in deviation sequences in each response, respectively, hence
and
. As discussed earlier, the
value is 0.5, which means all factors have equal importance. In Equation (3),
i.e., the grey relational coefficient for each run was determined, as shown in
Table 4. Grey grades,
were determined by taking the average of grey coefficients using Equation (5).
Table 4 shows that the second experiment gives the highest relational grade, which means it agrees with the theoretical results.
It has been concluded that multiple performance characteristics have been transformed into optimized grey grades. As it is an orthogonal experimental design, it is easy to determine the main effects of input parameters with different grey grades. For instance, the average grey grade value for gas flow rate (F) at the first, second, and third levels were computed by taking the mean value of grey grades from the 1st to 3rd, 4th to 6th, and 7th to 9th experiments, as presented in
Table 5.
Similarly, the mean grey grades can be determined for the other factors, namely welding current and gap distance. The results for each factor at each level have been summarized in
Table 5. As discussed earlier, product quality will be nearer to the optimum as the value of grey grades increases. Therefore, for adequate performance measures, higher grey grades are desirable. The optimal combination of the factors for optimum bead geometry is F1, A2, and G3, as presented in
Table 5. The signal-to-noise ratios for each factorial combination relative to optimization use Equation (6), as provided in
Table 6.
In the above equation,
n is the number of experiments and
shows the measured output response. Mean S/N ratios for the three factors at each level are computed and presented in
Table 6. MINITAB 16 software was used to plot the parameters’ main effects, as shown in
Figure 3. This figure illustrates the optimal condition of factors, i.e., F1, A2, and G3 (i.e., gas flow rate (F) = 10 L/min, welding current (I) = 80 A, and gap distance (G) = 2.5 mm), which is the same as obtained by grey relational analysis.
5. Finite Element Analysis (FEA)
The primary purpose of the finite element analysis method is to determine the solutions to different engineering problems with complex geometries to reduce the calculations and saves time. Experimentation can be used to find unknown factors and optimal combinations, but it requires additional testing and samples to study each behavior. This will has a high cost and time is needed to prepare each model and procure the equipment. The FEA technique predicts the approximate solutions to save time and cost within an acceptable range. Depending upon the accuracy level of the model, the optimized structure or shape is integrated into the finite element analysis. Different boundary conditions must be satisfied according to the structure and field variables. The steps are generally categorized into three main stages: (a) pre-processing, (b) solution, and (c) post-processing. Pre-processing includes the geometric and material properties of the elements, e.g., dimensions, element types, meshing the model, loadings, and boundary constraints. In the solution, the model is built to find the unknown variables. In post-processing, results are evaluated and interpreted.
5.1. Implementation of FEA (Material, Modeling, and Properties)
The present study analyzed the solid model using the finite element method (ANSYS APDL 2020R1) of the butt welded joint (mild steel). The model comprised two mild steel plates with dimensions 80 mm × 20 mm × 3 mm, welded at the ends using a V-shaped butt joint. The welded parts are modeled and analyzed to observe the temperature distribution and stresses induced during welding. The study was conducted by taking a single pass.
Figure 4a,b shows the solid and meshed geometric models. The temperature at each node was adapted in the welding environment for the element types, PLANE 55 (2D), containing four nodes, and SOLID70 (3D) with a single DOF. Meshing was performed in four areas (mm) as A1 (0.0012), A2 (0.0025), A3 (0.005), and A4 (0.0065) for better estimation. Transient thermal analysis was first performed by providing heat flux concerning time to predict the temperature change. Then, stresses and distortion were obtained by coupling the thermal and static structural analyses. The mesh size used in both steps was the same. It was assumed that the base metals and welded regions have the same thermal properties. There is no penetration, and overfilled welds are considered for simplicity.
The temperature-dependent properties, such as Poisson’s ratio, Young’s modulus, conductivity, density, and coefficient of thermal expansion, were characterized through differential scanning calorimetry and dynamic mechanical analysis, as illustrated in
Figure 5a–f. The melting temperature of mild steel is 1450 °C.
5.2. Thermal Loading
The heat flux varies with time after the thermal load is applied during thermal analysis.
Figure 6 shows that the load inputs are provided in three steps (t1–t3), and the total time is up to 1100 s. For the first ten seconds, a ramp input is given, followed by a step input up to 100 s, and then again, a ramp input is given up to the end. When welding begins, the heat steadily increases until it reaches a particular value. For a short period, it remains stable (no change) and then decreases steadily until the temperature of the welded plates reaches room temperature.
The input heat was determined by the product of arc efficiency, provided voltage, and welding current, which are 0.7 V, 15 V, and 80 A, respectively, in the present study. The maximum heat flux determined is 0.62 × 106 W/m2. The node solutions for temperature in this step were regarded as inputs into the stress analysis. This is because the residual stresses produced in heating and cooling periods result from temperature fields.
5.3. Boundary Conditions
The boundary conditions should be set for the thermal and static analysis, which can be applied to the plates for different parameters. The boundary conditions for other parameters used in this study are the following.
5.4. Thermal Analysis
The convection heat was applied to the surface area of the plates, which are meshed in identical elements in the symmetric model, as illustrated in
Figure 7. The Gaussian heat source was applied on surface areas A5 and A10.
The welded region where the plates are in contact is supposed to be insulated, i.e., no heat flow across this region. The value of the coefficient of convection is 20 W/m2 °C. The uniform temperature, 220 °C, was given to other remaining areas.
7. Heat Transfer and Mechanical Analysis
Heat transfer and mechanical analysis is the first step in the FEA while running the TIG welding simulation using ANSYS. The following equation can describe the formulation in FEA for the elements:
The heat conduction equation must be integrated with regard to time for this study. The Crank–Nicolson/Euler theta integration method is used to solve these system equations. This element type can conduct heat in three dimensions. The heat conduction versus time equation must primarily be integrated into this kind of analysis. The heat is provided by the electrode as the heat flux and is considered input for heat transfer from the electrode to the plates based on setting factors and arc efficiency. The heat flux density at any point from the heat source center can be computed using the following equation:
Since, at the y-axis, the value of
x is zero, Equation (8) can be written as:
where
is the maximum heat at the center of the heating source, sub-degree
R is the effective radius of the heating source, and
r is the distance between the center of the arc heating source and point A. The heat source model estimates temperature with finite element analysis in TIG welding. For stress analysis, heat transfer analysis was carried out to compute the temperature at the nodes as the time function, and then the structural analysis was performed using the temperatures attained from the heat analysis. The input values of the selected factors and constant factors used in the simulation are given in
Table 7.
The constant voltage of 15 V was supplied to each experiment, and the selected value for welding current achieved from the grey relational analysis was 80 A. The value of total heat flux at the welded region was calculated using the following equation:
where
Q is heating flux density, and
A is the area of the welded region, which is the same for all the experiments. By computing the value of
Q using the equation and the area of the welded triangular region (
and
), the total heat flux was determined (
q = 0.62 × 10
6 W/m
2). The FEA for the thermal and structured analysis of TIG welding using ANSYS Mechanical APDL works according to the following Equations:
The notations used in the above equations “SDPBS.
http://sdpbs.math.uwm.edu/about.php (accessed on 22 December 2022)” are explained in
Table 8. The total time to arrive at the solution for temperature was 1100 s, illustrated in
Figure 8, and the number of sub-steps in each iteration was 50. The x-axis and y-axis are labeled as cumulative iteration numbers and absolute convergence norms.
The absolute convergence norm uses normalized values to quickly solve the variables using the Newton–Raphson method. A solution within the range means the function value should be higher than the value of the vector norm of the function for each sub-step. It is observed that the HEAT CRIT values are higher than the HEAT L2 (vector norm) at every sub-step, so it can be concluded that the solution is within tolerance and is correctly converged.
9. Conclusions
In this study, we utilized the Taguchi-grey relational analysis approach to optimize the bead geometry in TIG welding. We also conducted numerical simulations using finite element analysis to analyze temperature and stress variations. Our results showed that the welding current was the most influential factor affecting the performance measures, followed by the flow rate. We found that the optimal values to optimize the bead geometry were a heat flux of 0.62 × 106 W/m2, a gas flow rate of 10 L/min, a welding current of 80 A, and a gap distance of 2.5 mm. The maximum deformations induced by the thermal stress along the x-axis, y-axis, and z-axis were 0.044 mm, 0.0047 mm, and 0.556 mm, respectively. Less deformation occurred along the y-direction as there was no motion of the heat source. By changing the voltage, efficiency, and welding speed, we observed that the temperature field increased with voltage and efficiency but decreased by increasing the welding speed of the torch. Our models were validated through additional experiments.
Overall, this study demonstrates the potential of using Taguchi-grey relational and finite element analyses to optimize the bead geometry in TIG welding. By achieving uniform temperature and stress distribution, we can improve the mechanical properties and quality of the weld while saving resources and reducing costs. In future studies, we can consider heat dissipation or losses and multiple passes for even more precise modeling.