Abstract

This paper presents the results obtained and the deductions made from an analytical modeling involving friction stir welding of Al 6061-T6. A new database was developed to simulate the contact temperature between the tool and the workpiece. A second-order equation is proposed for simulating the temperature in the contact boundary and the thermal history during the plunge phase. The effect of the preheating temperature on the process time was investigated with the proposed model. The results show that an increase of the preheating time leads to a decrease in the process time up to the plunge and the preheating stage. Results of modeling the thermal history showed acceptable agreement with corresponding experimental data from the literature.

1. Introduction

Friction stir welding (FSW) is a quite innovative process to join metals such as aluminum, magnesium, and steel alloys in the solid state. It was invented by The Welding Institute (TWI) in 1991 [1]. The process consists of a rotating tool which is plunged into two sheets or plates tightly abutted in a line, along which the rotating tool is traversed. During the process, heat is generated by plastic deformation as well as by the friction between the tool and the sample. The sample is eventually joined by the stirring action of the softened (but always solid) material. The advantages of FSW include low residual stress, low energy input, and fine grain size compared to more conventional fusing welding methods. One way to optimize the process is to utilize the thermomechanical modeling of FSW.

A number of academic and industrial institutions have made efforts to develop numerical codes for FSW. Although FSW is simple in concept, the physics behind the process is complex, which includes mechanical heat generation and heat and mass transport. The large strains and strain rates make observing the details of the process difficult, which makes process modeling attractive or essential for understanding it. Colligan [2], Li et al. [3], Guerra et al. [4], and Schmidt et al. [5] studied the material flow in FSW, which is useful for the investigation of the mechanism of FSW. Xu et al. [6, 7] developed a numerical model to capture the main thermomechanical features occurring in the FSW process for the further investigation of the mechanism of FSW.

The microstructure of the friction stir welded 7050-T651 was investigated by Su et al. [8]. Metallurgical, hardness, and quantitative disperse X-ray measurements were performed by Sutton et al. [9] to show that a segregated, banded microstructure consisting of alternating hard particle rich and hard particle poor region is developed in friction stir welding. The superplastic behavior of the welded sections and the one of the base metal were compared by Salem et al. [10]. The feasibility of FSW for joining copper was demonstrated by Lee and Jung [11]. It was shown that the transverse tensile stress reaches about 87% of that of the base metal, which is slightly higher than that of EBW copper joints. It should be noted that the maximum temperature created by the FSW process ranges from 80 to 90% of the melting temperature [12, 13], so that welding defects and large distortions commonly associated with fusion welding are minimized or avoided. The heat efficiency in FSW is 95%, which is higher than that in the traditional fusion welding [14]. FSW contains three stages: plunging, preheating, and traversing. The plunge stage in the FSW process is extremely critical since most of the initial thermomechanical conditions are generated and the material undergoes significant transformation due to the high temperatures and stresses involved in the process. The highly dynamic nature of this phase makes it a challenging research area. A thorough understanding of the plunge stage is also important in the development of tools and processes for successfully stir welding of high strength alloys like steel and titanium-based alloys, as most of the tool wear occurs during this phase [15, 16]. Although there are few experimental studies that focus on the plunge stage [17, 18], there are not many numerical models that focus on the thermomechanical conditions developed during the plunge phase. Many of the FEM-based numerical modeling of FSW simulates either the welding stage or the dwell and the welding stages [1922]. One of the principal difficulties in simulating the plunge stage is the severe mesh distortion of the finite element model which often results in sudden termination of the program. This problem has been successfully dealt with in this research.

A great effort was spent on the analytical models of FSW in order to formulate the heat generation. Russell and Shercliff [23] based the heat generation on a constant friction stress at the tool-matrix interface, equal to the shear yield stress at elevated temperature, which is set to 5 pct of the yield stress at room temperature. The heat input is applied as a point source or line source, and the solution is modified to account for the limited extent of the plate width. Schmidt et al. [19] formulated an analytical model for the heat generation based on different assumptions at the tool-matrix interface. Using this analytical model in comparison with experimental data, the authors suggested that a sticking or close to sticking contact condition is present at the tool/matrix interface. The same authors presented in [24] a new thermal pseudomechanical model where the heat generation is described as a surface flux governed by the material flow stress which in turn depends on the nonuniform temperature at the contact tool matrix interface. In their work, the solution was obtained numerically, and a good correlation with the experimental data was obtained by adjusting only the heat-transfer coefficients’ contact resistance. This was the first attempt to develop a thermal model where the total heat generation is not an input parameter but is actually a result of the model itself. Vilaca et al. [25] developed an analytical model for the FSW process, which simulates the asymmetric temperature field developed below the tool shoulder due to the combination of the rotation and linear velocities and hot-to-cold welding conditions. The main limitation of such a model is the need to calculate, for every different welding condition, the thermal input by an iterative procedure, which refers to experimental data. However, it can be considered a good research tool which allowed for the establishment of a systematic relation between the welding parameters and the final mechanical and metallurgical characteristics [26].

In the present study, analytical modeling of the plunge and the preheating stages has been simulated by a newly developed equation. To validate our simulation the results of the thermal history in the plunge section from modeling are compared with the experimental data which were obtained by Hwang et al. [27].

2. Governing Equation

The heat flux in FSW is primarily generated by the friction and the deformation process. Thus, the governing equation for heat transfer in FSW is given by (1): where is the volumetric heat source term arising from plastic dissipation (W/m3), is the density of materials, is the mass-specific heat capacity, is the coefficient of thermal conductivity, is the temperature, and (=) is the gradient operator.

Ferro and Bonollo [28] suggested an equation to calculate the temperature distribution in the workpiece. Their equation was based on thermal analytical modeling, and both thermal and mechanical parameters were used; that is, where is friction coefficient, is rotational speed of the tool, is the function of the traverse speed for tool [28], is tool force, is contact temperature, is melting point, and the rest are tool parameters which are illustrated in Figure 1.

From this equation, it can be concluded that the contact temperature plays a major role in the determination of the temperature distribution in the whole sample. Hence, in the present work, a second-order equation is proposed for the temperature in the contact boundary relating to process times as well as mechanical parameters. The proposed equation is written as follows: where is contact temperature, is constant and here equal to workpiece initial temperature, is proportional to , and is proportional to . is a tool force, and is a thermal diffusivity. Furthermore, is friction coefficient, is rotational speed, is traverse speed, and is thermal conductivity. This equation is proposed based on the far-field temperature profiles in the joint welds [24], which the temperature raise shows almost second-order behavior. The perfect thermal contact is assumed for the simulation, in which the temperature for the tool interface and the work piece are equal in their contact. Calculated constants for the proposed model from the corresponding experimental condition which was conducted by Hwang et al. [27] are listed in Table 1.

3. Model Validation

For numerical modeling, thermal code was developed and written in visual FORTRAN. Workpiece and the contact boundary were discretized by the FDM method. The 3D domain was discretized with relatively fine and structured mesh, and to reduce the computational time, the half of the domain is used in the simulation (by this simplification, the temperature difference in the advancing and the retrieving side of the tool was assumed to be constant). The heat conduction equation for the current simulation was solved by alternating direction implicit (ADI) method, which is most notably used to solve the problem of heat conduction in two or more dimensions and is an example of an operator splitting method. It is assumed that the material properties (i.e., the density and the thermal conductivity) are constant, and they are not depending on the temperature. The dimensions of the workpiece and the tool were read to program on the basis of Hwang et al. [27] experimental set up (see Figure 2). Thermal properties of Al 6061-T6 were read to the code and initialized the simulation. According to their work [27], the total time for the two stages of plunging and preheating are 83 and 51 seconds, respectively. So the model was set to run for 134 seconds. Results of time temperature for modeling and experimental data in the point near the pin (TC1) are illustrated in Figure 3.

It can be seen that although there are differences between the curves, the trend of the thermal behavior for the modeling and the corresponding experimental data are the same. Maximum deviation between the results occurs around the time 50–100 seconds, where there is transition from stage one (plunge) to the former (preheat). The next reason for the difference is the nature of modeling; there is always intrinsic deviation between results of modeling and experiments. Another reason is the effect of the fact that temperature-dependent parameters are not used in the modeling. Contour of thermal history is shown in Figure 4 with the initial temperature of 30°C and after 400 seconds.

4. Preheating

For materials with high melting point such as steel and titanium or high conductivity such as copper, the heat produced by friction and stirring may be not sufficient to soften and plasticize the material around the rotating tool. Therefore, it is difficult to produce continuous defect-free welds. In these cases, preheating or additional external heating sources can help the material flow and increase the process window. On the other hand, in materials with lower melting point such as aluminum and magnesium, cooling can be used to reduce extensive growth of recrystallized grains and dissolution of strengthening precipitates in and around the stirred zone.

Preheating of the workpiece before welding should be beneficial for improving welding speed and minimizing tool wear. It may be more simple and practical to preheat the initial plunge region of the workpiece before plunging the pin into the workpiece because the tool wear mainly occurs during the initial plunge period at the start of each weld. Moreover, preheating results in decreasing the process time in FSW. By preheating, the initial temperature for plunge increases, and the tool can easily penetrate the metal.

In the present study, the effect of the preheating temperature on the process time was evaluated using the proposed and validated model from the previous part. Five different cases used for the investigation are listed in Table 2.

In fact, as the preheating temperature increases, the initial temperature in the workpiece is increased. By this fact and our proposed model for contact temperature, five different preheating temperatures were read to the database. By assuming that the maximum temperature at the end of the preheating stage in the FSW process is 384°C (according to Hwang et al. experiments [27]), the total time to reach this temperature is evaluated. Results of modeling for five cases are illustrated in Figure 5.

As the preheating temperature increases, the total process time up to the end of the preheating stage is decreased. This results in low-cycle (lower time of the weld) high-quality weld. Moreover, since the temperature of the process in the beginning is high, the deformation of sample becomes easy and voids, and other defects can be avoided. The effect of preheating temperature on process time is illustrated in Figure 6.

From Figure 6, a simple linear behavior for relating preheating time on process time was deducted. This behavior is written as follows: where is the process time and is the preheating temperature. It should be noted that this equation is valid for the two stages of plunging and preheating.

5. Conclusion

An analytical model was developed to simulate the contact temperature in the friction stir welding (FSW). This second-order equation which contains thermal characteristics and welding parameters was compared and validated by experimental data in the literature. The effect of the preheating temperature on the process time was investigated. The results show that the increase of the preheating temperature not only develops the weld quality, but it also decreases the total process time.