Computational Modeling of Heat Transfer and Visco-Plastic Flow in Friction Stir Welding

With a focus to develop a quantitative understanding of the FSW process, a comprehensive three dimensional heat transfer and plastic flow model is developed. The model can predict variables such as temperature and velocity fields and torque based on the given welding parameters like weld velocity, tool rotational speed and axial pressure. It considers tool design dependent spatially variable heat generation rates, deformational work, non-Newtonian viscosity as a function of local strain rate, temperature and the nature of the material and temperature dependent thermal conductivity, specific heat capacity and yield stress. It is shown that the temperature fields, cooling rates, the plastic flow fields and the geometry of the thermo-mechanically affected zone (TMAZ) can be adequately described by solving the equations of conservation of mass, momentum and energy in three dimensions with appropriate boundary conditions and constitutive equations for viscosity. The model is tested for four different alloys: (1) AA 6061-T6, (2) 1018 Mn steel, (3) 304L stainless steel and (4) Ti-6Al-4V which have widely different thermophysical and rheological properties. Numerically computed temperature fields, variations of peak temperatures with FSW variables and TMAZ geometry were compared with the experimental results. Currently, due to unknown parameters in existing transport phenomena based models, the computed temperature and velocity fields and torque may not always agree with the corresponding experimentally determined values and may not show the same trend as experimental results for a range of welding variables. Here, it is shown that this problem can be solved by combining the rigorous phenomenological process sub-model with a multivariable optimization scheme called Differential Evolution. The values of the uncertain model input parameters from a limited volume of independent experimental data which includes temperature measurements obtained using thermocouples and torque measured using dynamometers. This approach resulted in agreement between the phenomenological model and the experimental results with a greater degree of certainty. It is tested for FSW of: (1) dissimilar AA 6061-T6 to AA 1200, (2) 1018 Mn steel and (3) Ti-6Al-4V. Independent thermocouple and dynamometer measurements are also used for validation and verification of results. Improvement in the reliability of the numerical model is an important first step towards increasing its practical usefulness. Also, one of the reasons why current models do not find extensive applications is because they cannot be used to tailor weld attributes. The aim of the present research is to develop a reliable bi-directional model which can find wide use in manufacturing and process control. It is shown that by coupling a reliable model with an evolutionary search algorithm, we can find multiple sets of welding parameters to achieve a target peak temperature and cooling rate in welds. The model is tested for dissimilar welds of AA 6351 and AA 1200. FSW is being increasingly used for dissimilar metal joining. Models are needed to calculate the redistribution of alloying elements when two alloys with dissimilar alloying element contents are joined. The transport and mixing of magnesium from Mg-rich AA 6061 alloy into a commercially pure aluminum AA 1200 was examined experimentally and numerically at various locations in the welded workpiece. The concentration of the solute is measured in transverse cross-sections across the weld-center line at various depths from the top surface of the workpiece. The measurement was done using electron probe micro-analysis (EPMA) of polished transverse-cut...