Phenomenological yield criteria are generally described by many material parameters. A technique to identify these parameters is required to find the best fit to the results of the mechanical tests. The parameter identification by the classical simulated annealing technique is presented in this paper. This algorithm, based on works by Metropolis et al, is a global optimization method that distinguishes between different local optima to reach the global optimum. The anisotropic model used in this study is the one proposed by Cazacu et al. To prove the efficiency of the proposed algorithm, the material parameters of Ti6Al4V titanium alloy are identified and compared with those obtained using different identification procedures and the same experimental data.