Friction stir welding (FSW) process of aluminum alloys was investigated using a two-dimensional Eulerian formulation coupling viscoplastic flow and heat transfer and strain hardening. The thermal equation for the temperature was modified to stabilize temperature distribution using a Petrov-Galerkin method. The evolution equation for strength was calculated using a streamline integration method. Predicted strength was compared with experiments. Based on crystal plasticity, texture evolution was predicted during FSW of AA6061.