Non-linear regression combined with Markov chain Monte Carlo method is utilized. A combination of Voce-modified JC can predict hardening behavior of DP600 accurately. Rousselier model is calibrated to simulate uniaxial tension at different strain-rates. The onset of rapid void growth depends on the hardening rate at large strains. The predicted geometry of the fracture surface is dependent on the hardening function.