There is currently a need for an efficient numerical optimization strategy for the quality of friction stir welded (FSW) joints. However, due to the computational complexity of the multi-physics problem, process parameter optimization has been a goal that is out of reach of the current state-of-the-art simulation codes. In this work, we describe an advanced meshfree computational framework that can be used to determine numerically optimized process parameters while minimizing defects in the friction stir weld zone. The simulation code, SPHriction-3D, uses an innovative parallelization strategy on the graphics processing unit (GPU). This approach allows determination of optimal parameters faster than is possible with costly laboratory testing. The meshfree strategy is firstly outlined. Then, a novel metric is proposed that automatically evaluates the presence and severity of defects in the weld zone. Next, the code is validated against a set of experimental results for ½” AA6061-T6 butt joint FSW joints. Finally, the code is used to determine the optimal advancing speed and rpm while minimizing defect volume based on the proposed defect metric.

Recognition of the friction stir welding process is growing in the aeronautical and aero-space industries. To make the process more available to the structural fabrication industry (buildings and bridges), being able to model the process to determine the highest speed of advance possible that will not cause unwanted welding defects is desirable. A numerical solution to the transient two-dimensional heat diffusion equation for the friction stir welding process is presented. A non-linear heat generation term based on an arbitrary piecewise linear model of friction as a function of temperature is used. The solution is used to solve for the temperature distribution in the Al 6061-T6 work pieces. The finite difference solution of the non-linear problem is used to perform a Monte-Carlo simulation (MCS). A polynomial response surface (maximum welding temperature as a function of advancing and rotational speed) is constructed from the MCS results. The response surface is used to determine the optimum tool speed of advance and rotational speed. The exterior penalty method is used to find the highest speed of advance and the associated rotational speed of the tool for the FSW process considered. We show that good agreement with experimental optimization work is possible with this simplified model. Using our approach an optimal weld pitch of 0.52 mm/rev is obtained for 3.18 mm thick AA6061-T6 plate. Our method provides an estimate of the optimal welding parameters in less than 30 min of calculation time.