Abstract: A plastic damage constitutive model was proposed with special concerns on simulating the stress-triaxiality-dependent hardening phenomena.The stress triaxiality was defined in terms of stress invariants.A subroutine was developed for validation of the constitutive model at local level.Comparison between the numerical results and experimental data indicates the validity of the proposed model and the related numerical scheme.