a b s t r a c tInitial and subsequent yield surfaces for 6061 aluminum, determined by a method of automated yield stress probing, are presented in the 2D (r zz À r hz ) and 3D (r hh À r zz À r hz ) stress spaces. In the (r zz À r hz ) space, yield surfaces at small pre-strains show the noses and unapparent cross effect. At larger pre-strains, they become ellipses with positive cross effect. In the (r hh À r zz À r hz ) space, the initial yield surface is not well described by von Mises yield criterion due to material anisotropy. The yield surfaces of various torsional pre-strains show obvious rotation around the r zz axis but they do not rotate when subjected to axial pre-strains. Therefore, the rotation behavior of yield surface is pre-strain path dependent. The rotation of yield surfaces in the 3D space is the emphasis of the present paper. Coupled axial-torsional behavior subjected to torsion after axial pre-strain are also presented for the same material that is used to determine the yield surfaces. This information is useful for verification of constitutive models.