Forming limit diagram (FLD) is one of the formability criteria which is a plot of major strain versus minor strain. In the present study, Gurson-Tvergaard-Needleman (GTN) model is used for FLD prediction of aluminum alloy 6061. Whereas correct selection of GTN parameters’ is effective in the accuracy of this model, anti-inference method and numerical simulation of the uniaxial tensile test is used for identification of GTN parameters. Proper parameters of GTN model is imported to the finite element analysis of Nakazima test for FLD prediction. Whereas FLD is dependent on forming history and strain path, forming limit stress diagram (FLSD) based on the GTN damage model is also used for forming limit prediction in the numerical method. Numerical results for FLD, FLSD and punch’s load-displacement are compared with experimental results. Results show that there is a good agreement between the numerical and experimental results. The main drawback of numerical results for prediction of the right-hand side of FLD which was concluded in other researchers’ studies was solved in the present study by using GTN damage model.