Tool Temperature and Process Modeling of Friction Stir Welding ()
1. Introduction
Friction stir welding (FSW) invented by the Welding Institute is a solid state joining technique initially intended to weld aluminum alloys [1] . Presently, it is used to weld many kinds of materials, such as aluminum, magnesium and titanium alloys and even plastics [2] . FSW welds materials in a solid state, so that it has many advantages over conventional liquid-solid welding methods: higher weld strength, fine grain size, low residual stress, no consumable materials and low energy input [3] [4] . This process consists of a rotating tool and plates or sheets to be welded, where the tool is plunged into the materials and traversed along a line to finish joining them. Despite the simple configuration, phenomena inside the materials under welding and at the interface between the rotating tool and materials are too complicated to fully understand because there is thermal softening of the materials and heat transmission between them at high temperatures.
It is of great importance to investigate these phenomena during FSW for sound joining of workpieces and high product qualities. Elucidation of the thermo-mechanical phenomena between tool and workpiece can lead to the optimization of welding parameters, effective tool design, and the application of FSW to new materials and products. For deep understanding of the FSW process, both experimental methods and numerical simulations can be adopted.
Welding parameters, material flows and tool/workpiece interface phenomena have direct influences on heat generation during FSW. Although much research has investigated the important aspect of welding phenomena based on temperatures measured at specific points of workpiece, only a few experimental studies have been reported about tool temperatures because very complicated settings are required for tool temperature measurements [5] [6] [7] [8] [9] . Accurate experimental data about tool surface temperature are necessary to understand exact temperatures at the tool/workpiece interface. They are also needed to precisely validate developed simulation models by comparing measured and calculated temperatures. Therefore, a method and device for effectively measuring tool temperatures are essential to investigate FSW processes from theoretical and practical viewpoints.
Temperature analyses of the tool and workpiece during FSW have been conducted in many research studies. One way to optimize the process is to utilize the thermo-mechanically coupled modeling of FSW, considering thermal boundaries and friction at tool/workpiece interfaces. Hence, numerical models to simulate many phenomena that occur inside the material and tool and at their interfaces have been increasingly required. However, it is very difficult to establish a model that exactly meets the experimental results over a wide range of FSW conditions. FSW models include some thermal boundaries, and temperatures of tool and workpiece are sensitive against heat flows through the boundaries, so that the optimum values of heat transfer coefficients are necessary for these thermal boundaries. However, almost no study has investigated the precise thermal boundary conditions for FSW modeling because relatively low traverse speeds have been adapted for the FSW simulations.
Relationships between process parameters, such as rotational and traverse speed, and inside temperatures are important information for production sites. Although some studies investigated the relationship between the traverse speed and FSW temperatures for aluminum alloys, they analyzed the temperatures at low traverse speed under around 3.5 mm/s [10] [11] . For more practical use, simulation and experimental data at higher traverse speeds should be discussed.
In this study, a thermos-mechanically coupled process model was developed to investigate FSW phenomena inside a tool and workpiece and at their interface using the DEFORM-3D software. As a workpiece, 6061-T6 aluminum alloy was employed. The system of FSW process model includes several thermal boundaries. Among heat flows through these thermal boundaries, the heat transfers into the exterior of the system become more sensitive to tool and workpiece temperatures than the heat transfers within the system. This paper especially focused on a heat transfer coefficient at the workpiece bottom, and optimized it through experiments and finite element (FEM) analyses for a reasonable simulation model to be established. Although other FEM studies have investigated FSW process phenomena under low tool traverse speeds corresponding to laboratory scale experiments, this study developed an analysis model that is able to investigate material internal phenomena at a wide range of traverse speeds for practical high speed welding. Both as fundamental data and reference data for simulation modeling, the tool temperatures during FSW were measured with a special tooling system with imbedded thermocouples within a tool. Then, the heat transfer coefficient at the workpiece bottom was optimized and the accuracy of developed FEM model was validated. Finally, the temperatures and stress distribution around workpiece/tool interfaces were investigated with the developed simulation model.
2. Procedure
2.1. Experiment to Investigate Tool Temperature and Welding Force
Several experiments with 6061-T6 aluminum alloy were conducted to measure tool temperatures and welding forces while butt welding of the aluminum sheets. The obtained results served as the fundamental data to investigate FSW phenomena both inside the tool and workpiece and at their interface, and are used to develop an FEM analysis model which accurately represents internal phenomena around the tool/workpiece interface under a wide range of traverse speeds.
The experiments were conducted on a Hitachi Power Solution 2D-FSW machine and traverse and axial forces were measured with load cells equipped under the machine table. The geometry of each workpiece was 200 mm long, 70 mm wide and 4 mm thick and the sides of two workpieces were butt welded. The joining surfaces were machined previously to meet each other without small gaps. The tool had a shoulder 12 mm in diameter and 3˚ in taper angle, and a threaded probe 6 mm in diameter and 3.8 mm in height.
With respect to welding temperature, tool side measurement was adopted because tool temperatures are under nearly stable conditions during FSW process.
The tool inside temperatures during FSW were measured with MULTI INTELLIGENCE R tooling system for FSW. Figure 1 shows a general view of the temperature measuring tool used for experiments. Figure 2 shows the thermocouple positions imbedded in the tool. The thermo-couples were arranged at three different points to measure the temperatures of the probe-tip, probe-root and shoulder part respectively. The tool system was precisely assembled and thermo-couples were imbedded at positions around 1 mm from the tool surface to measure temperatures as close to tool/workpiece interface temperatures as possible.
Table 1 shows experiment conditions. To obtain data not only for fundamental
Figure 1. General view of temperature measuring tool.
Figure 2. Cross-section view of positions for thermocouples imbedded in tool.
Table 1. Experiment conditions to measure temperatures and forces.
study but also for practical welding, traverse speeds were set for a wide range from 1 mm/s to 20 mm/s, close to the practical traverse speed at production sites. In addition, the tool plunge speed was set to 0.5 mm/s, and then the tool was held for 5 s for the dwelling before traversing.
2.2. FEM Simulation Modeling
To clarify the thermal and mechanical phenomena of workpieces and tools during the FSW process, it is essential to understand temperature distributions, stress fields and strains inside tools and workpieces or at their interfaces. These inside and interface phenomena are difficult to investigate only by experiments. Hence, a numerical analyses model with FEM that is able to cover a wide range of welding speeds considering precise thermal boundary conditions was developed in this study.
The commercial software DEFORM-3D was used for the FEM simulation, and Lagrangian-code was employed for the analyses. The workpiece was set to a continuum model of 70 × 70 × 4 mm with meshing into about 32 thousand tetrahedral elements. The element sizes around the tool/workpiece interfaces were set to be about 0.1 - 0.3 mm in length, finer than those in the surroundings. The tool in this model had the same geometries as that used in experiments except that the probe height was 3.6 mm. That is, the probe height was 0.2 mm smaller in the FEM simulation than in experiments and mesh breaking at model’s bottom while re-meshing was avoided by an increase in the clearance between the tool tip and workpiece bottom. The probe part of the model was a cylinder pin without threads, giving priority to analysis stability during FSW simulations. There were 15 thousand tool elements, and tool/workpiece interface elements were set to be finer like the workpiece modeling. The workpiece model had a rigid plastic body and tool model a rigid body, and thermo-coupled viscous-plastic analyses were conducted.
Figure 3 represents the flow stress of aluminum alloy 6061-T6 used for the FEM analyses. This flow stress model includes not only thermal softening showing a dramatic decrease in flow stress at high temperatures but also strain ratehardening. First, quasi-static compression tests for 6061-T6 aluminum alloy were conducted under three different temperature conditions of room temperature (about 20˚C), 100˚C and 200˚C to obtain the stress-strain relationship. The geometry of a cylinder specimen was 10 mm in diameter and 10 mm in height. The stroke speed was set to 1 mm/min for quasi-static condition. The stress-strain relationships over 200˚C were based on references [12] [13] . Figure 4 represents the relationships between temperature and 0.2% yield strength at a strain rate of 0.002 s−1 rearranged from Figure 3. The strain rate dependence in the model can follow Johnson and Cook constitutive model shown in Equation (1), and its strain rate parameter C was set to 0.02 [14] to design the temperature and strain rate-dependent flow stress model.
(1)
Figure 3. Flow stress model of aluminum alloy 6061-T6 for temperatures of 20˚C, 200˚C, 316˚C and 428˚C and strain rates of 0.002 s−1, 100 s−1 and 100,000 s−1.
Figure 4. Relationship between temperature and 0.2% yield strength of 6061-T6 obtained from compression test results and literature values at strain rate of 0.002 s−1.
A: yield stress, B: strain hardening constant, ε: equivalent plastic strain,
n: hardening exponent, C: strain rate constant,
: dimensionless strain rate,
: homologous temperature, m: temperature exponent.
2.3. Optimization of Thermal Boundary Condition
A system of FSW process model includes several thermal boundaries in itself. A thermal boundary within the system is at the tool/workpiece interface, where the heat flowing between them remains inside the system itself. As for thermal boundaries against the exterior of the system, there are heat convections from the tool surface and the workpiece-top surface to the air, and heat conductions from the workpiece-bottom surface to the back plate. The amount of heat transfer from a solid body to the air is relatively small. In contrast, the heat transfer between two solid bodies in contact is not negligibly large when there are large temperature differences between them. For this reason, the heat transfer coefficient between them has a large influence on the amount of heat removed from the system and corresponding temperature reduction. Therefore, its value could be sensitive to analysis results, so that the heat transfer coefficient between the workpiece bottom and back plate was optimized in this research. Table 2 shows the thermal boundary conditions used for modeling.
The optimization of the workpiece-bottom heat transfer coefficient for the FEM analysis model was carried out by comparing temperatures numerically and experimentally obtained. The traverse speed and the tool rotation speed were set to 5 mm/s and 1000 min−1, respectively. Five different values of the heat transfer coefficient at the workpiece bottom surface were set between 1.0 to 5.0 kW/m2/K for optimization. For modeling a convective heat transfer as the thermal boundary condition, Equation (2) was used.
(2)
q: heat flux, h: heat transfer coefficient, Tw: object1 surface temperature,
Tf :ambient temperature or object 2 temperature.
The constraint condition of the workpiece was that the bottom surface of the workpiece was in contact with the back plate without friction and the side surfaces are fixed in all directions. Workpiece and tool surfaces conduct heat transfer against the air of ambient temperature 20˚C. The friction model between the tool and workpiece was followed by a shear friction model shown in Equation (3), and the frictional coefficient m = 1.0 was employed [15] .
(3)
: shear stress, m: friction coefficient, k: shear yield stress.
After the heat transfer coefficient was optimized for the workpiece bottom surface at a traverse speed of 5 mm/s, the validation analyses were carried out at different traverse speeds to compare the analysis and experiment results of tool temperatures and traverse force.
3. Results
3.1. Tool Temperature and Welding Force Measurement Results
Figures 5-8 represent the data profile of tool temperature histories while welding at traverse speeds of 1, 5, 10 and 20 mm/s, respectively. The tool rotational speeds were 1000 min−1. In each graph, temp1 is the temperature at the probe tip, temp2 at the probe root and temp3 at the tool shoulder. Each graph also indicates the welding forces of the axial and traverse directions. It should be noted that no defect was visually detected for the above four welding conditions.
Table 2. Boundary conditions used for analysis model.
Figure 5. Data profile of tool temperatures and welding forces (1000 min−1, 1 mm/s).
Figure 6. Data profile of tool temperatures and welding forces (1000 min−1, 5 mm/s).
Figure 7. Data profile of tool temperatures and welding forces (1000 min−1, 10 mm/s).
The above results confirmed that temperatures of each portion increased monotonously from the tool plunge stage to the dwelling (rotation holding)
Figure 8. Data profile of tool temperatures and welding forces (1000 min−1, 20 mm/s).
stage. Here, traversing began after a 5-second tool rotation hold. After the traverse began, probe tip temperatures (temp1) became almost constant and this steady state continued until all processes finished. On the other hand, the temperatures of the probe root (temp2) and tool shoulder (temp3) increased continuously, but slightly throughout the traverse process. Generally, the measured temperatures of tool tips were about 100˚C to 150˚C higher than those of tool shoulders.
As for welding forces, they were almost constant during the traverse process at traverse speeds of 1, 5 and 10 mm/s, but at a traverse speed of 20 mm/s they kept increasing until the process finished.
Figure 9 shows the relationship between tool traverse speeds and tool temperatures. The results showed that probe tip temperatures were about 100˚C - 120˚C higher than shoulder temperature, and about 160˚C - 180˚C higher than probe root temperature. It was indicated that tool temperature at each measurement portion decreased by about 30˚C to 50˚C with an increase in tool traverse speed from 1 to 20 mm/s. In contrast, Assidi et al. [5] reported that tool temperatures were almost constant and nearly independent of traverse speed in a range of low tool traverse speeds from 1 to 3 mm/s. The plastic work and friction work consumed near the tool and at the tool/workpiece interface respectively can be regarded as a kind of moving heat source in high-speed welding. Because the convective heat transfer from the tool to moving workpiece increased with traverse speed, the tool temperature would decrease with increasing traverse speed.
Figure 10 shows the relationship between tool rotation speeds and tool temperatures for a traverse speed of 5 mm/s. Tool temperature increased with tool rotational speed. The probe tip temperature increased up to about 500˚C, close to the workpiece solidus temperature 582˚C. On the other hands, shoulder temperature and probe root temperature were likely to converge to about 370˚C and 300˚C, respectively at high tool rotational speeds.
Figure 9. Tool temperatures and tool traverse speeds (rotation: 1000 and 750 min−1).
Figure 10. Tool temperatures and tool rotational speeds (at traverse speed of 5 mm/s).
Figure 11 shows the results of the force measurements. Each value represents the data in the middle of welding path. It was clearly seen that the welding forces increased linearly with traverse speed. These characteristics were the same as those reported for aluminum alloy 2024-T3 [11] [16] [17] [18] .
3.2. FEM Model Optimization Result
Figure 12 shows analysis and experiment results for the tool-tip and shoulder temperatures and the welding forces in the traverse direction. Welding conditions were a traverse speed of 5 mm/s and rotational speed of 1000 min−1 for both analyses and experiment, and data 1.0 second after traversing were compared in the figure.
In the experiment result, the tool-tip temperature was 462˚C, the shoulder temperature 313˚C and the traverse force 1113 N. As for the analysis results calculated for five different heat transfer coefficients, tool-tip temperature decreased by 67˚C with increasing heat transfer coefficient from 1.0 to 5.0 kW/ m2/K. On the other hand, the shoulder temperature changed little with heat transfer coefficient. These suggest that the amount of heat flows from workpiece
Figure 11. Welding force measurement results.
Figure 12. Experiment and analysis results for optimizing heat transfer coefficient.
bottoms to the back plate affects tool-tip temperatures sensitively since the tool-tip is close to the workpiece bottom. In contrast, the traverse force increased with heat transfer coefficient. This is because the material flow stress increased with decreasing tool tip temperature.
Figure 13 shows the error ratios between the analysis and experiment results for five values of heat transfer coefficient. The error ratio was defined by Equation (4),
(4)
where Tt,a and Tt,m are tool tip temperatures analyzed and measured, respectively, Ts,a and Ts,m are shoulder temperatures analyzed and measured, respectively, and Ft,a and Ft,m are traverse forces analyzed and measured, respectively. The error ratio was the minimum at 2.0 kW/m2/K, the optimum value of heat transfer coefficient. It was also confirmed from Figure 12 that analysis results showed
Figure 13. Difference between analysis and experiment results.
the best agreement with the experiment ones for a heat transfer coefficient of 2.0 kW/m2/K.
This optimization process suggests that workpiece-bottom thermal boundary conditions remarkably affect tool-tip temperatures and welding forces, i.e. material inner flow stress. The optimized value of heat transfer coefficient for the workpiece bottom obtained in this study was quite different from values of 0.2 kW/m2/K [15] and 11 kW/m2/K [19] used in other studies. This indicates that the heat transfer coefficient may differ under different experiment settings, so that its value is needed to be optimized for each condition. Slight differences of heat transfer coefficients affected the analysis results of tool temperatures and traverse force. This also indicates that the thermal boundary condition differences while actual welding may affect workpiece and tool internal mechanical phenomena. These suggest that, in actual FSW processes, optimum welding settings are to be determined by considering not only tool rotational and traverse speeds but also heat flows from the workpiece-bottom to the back-plate.
3.3. Validation Result
Validation was conducted with the developed analysis model for a low traverse speed of 1 mm/s to a high speed of 20 mm/s, comparing the analysis data with experiment results. Figure 14 shows the validation result for tool temperatures, and Figure 15 shows for the welding forces. Both experiment and FEM data were acquired 1.0 second after the start of tool traversing. The tool-tip temperatures were almost the same between the experiment and analysis results. The shoulder temperatures of analyses and experiments had the same tendency with the traverse speed, though the analysis results became about 30˚C to 40˚C higher than the experiment ones. The analysis results of traverse force had good agreement with the experiment results. Then, the developed analysis model is able to simulate FSW process at a wide range of tool traverse speeds. This validation indicated that developing a FSW analysis model that calculate reasonable tool temperatures and process parameters needs accuracy of its thermal boundary conditions, especially heat transfer coefficient at the workpiece bottom.
Figure 14. Validation analysis result for tool temperature.
Figure 15. Validation analysis result for traverse direction welding force.
4. Discussion
Internal thermal-mechanical phenomena around the tool/workpiece interface were investigated with the developed analysis model. The analyses for a low traverse speed of 5 mm/s and a high speed of 20 mm/s were performed to investigate the temperatures and stress distributions around the tool/workpiece interfaces. Figure 16 and Figure 17 show the results of 5 mm/s, and Figure 18 and Figure 19 the results of 20 mm/s. The tool rotational speed was fixed to 1000 min−1. Each figure represents values on a cross-section of tool and workpiece along welding line.
Figure 16 and Figure 18 show that temperature inside the probe was around 440˚C at the center of the tip surface and 400˚C at the probe root. Calculated temperature distribution proved that the measured temperature of temp1 is almost the same as temperature at the center of tip surface. On the other hand, the shoulder had a steep temperature gradient along the tool axis, and thus, the temperatures of temp2 and temp3 corresponding to the measurement points inside the shoulder 1.1 mm from the surface were different from the shoulder surface temperatures by about 40˚C to 50˚C. It should be noted that the shoulder
Figure 16. Cross-section showing temperatures (traverse: 5 mm/s, rotation: 1000 min−1).
Figure 17. Cross-section showing stresses (traverse: 5 mm/s, rotation: 1000 min−1).
surface temperatures are higher than temperatures measured with thermocouples.
With regard to workpiece inner temperatures comparing the traverse speed of 5 and 20 mm/s, the tool front temperature distributions were definitely different, showing that the temperature was lower for 20 mm/s than for 5 mm/s. This is because cold workpiece rapidly come close to the tool at higher traverse speed before sufficient temperature increase by plastic and friction work. As for the tool front stresses, the maximum stress was about 70 MPa at 5 mm/s, whilst it exceeded 110 MPa at 20 mm/s. For this reason, a zone of higher stress appeared in front of the tool at 20 mm/s, showing that the flow stress of aluminum alloy 6061-T6 is critically temperature-dependent.
Figure 18. Cross-section showing temperatures (traverse: 20 mm/s, rotation: 1000 min−1).
Figure 19. Cross-section showing stresses (traverse: 20 mm/s, rotation: 1000 min−1).
5. Conclusions
Tool temperatures were measured during FSW, and a thermal boundary condition was optimized with these experiment results to develop an analysis model. The following was learned through this study:
1) The heat transfer coefficient of workpiece bottom was optimized to h = 2.0 kW/m2/K. Slight differences of heat transfer coefficients affected the analysis results for tool temperatures and traverse force, indicating thermal boundary conditions while actual welding may affect the workpiece and tool internal mechanical phenomena.
2) The analysis modeling considered with an optimum value of workpiecebottom heat transfer coefficient led to the reasonable model that can analyze FSW phenomena through a wide range of tool traverse speeds.
3) Analysis results with the developed model show that temperature inside the probe was around 440˚C at the center of the top surface and 400˚C at the probe root. Calculated temperature distribution proved that the measured temperature with the thermocouple at the probe tip is almost the same as the surface temperature. On the other hand, the shoulder had a steep temperature gradient along the tool axis, and thus the measured temperatures inside the tool were different from the surface temperatures by about 40˚C to 50˚C. It showed that the shoulder surface temperatures are higher than temperatures measured with thermocouples.
Acknowledgements
The authors are grateful to YAMAMOTO METAL TECHNOS CO., LTD. for technical support with measurements and tool precision machining.