The Gurson-Tvergaard-Needleman (GTN) damage model was developed basing on anisotropic yield criterion to predict the damage evolution for anisotropic voided ductile materials. Hill’s quadratic anisotropic yield criterion (1948) and Barlat’s 3-component anisotropic yield criterion (1989) were used to describe the anisotropy of the matrix. User defined subroutines were developed using the above models. Taking the benchmark of NUMISHEET’93 square cup deep drawing as an example, the effect of matrix plastic anisotropy on a ductile material was studied. The predicted result by Barlat’89-GTN model has a better agreement with the experimental data than that by Hill’48-GTN and the original GTN model.