The rolled aluminum alloy sheets usually have anisotropic properties that affect its ductile fracture behavior. In order to study the plastic anisotropy and ductile fracture behavior of aluminum alloy 6016 sheet, the Hill48 anisotropic yield criterion is introduced into the Gurson–Tvergaard–Needleman (GTN) model forming an anisotropic GTN model. In this model, two groups of Hill48 anisotropic yield constants, which are respectively calculated by Lankford's coefficient in different directions and the combination of the Lankford's coefficient and yield stress in different directions, are separately used to describe anisotropic mechanical behaviors. The failure void volume fraction is identified through microstructure analysis and the other damage parameters are calibrated through a finite element (FE) inverse calibration method. The abilities of these two groups of Hill48 anisotropic yield constants in predicting anisotropic mechanical behaviors of aluminum alloy 6016 sheet are examined by performing FE simulations of 0°, 45°, and 90° plate tensile tests, deep-drawing tests and Erichsen cupping tests. Results show that the Hill48 anisotropic yield constants calculated by the Lankford's coefficient in different directions have a better ability to predict local anisotropic plastic deformation especially the earing effect in the deep-drawing test, but the Hill48 anisotropic yield constants calculated by the combination of the Lankford's coefficient and yield stress in different directions have a better prediction ability in predicting global anisotropic mechanical behaviors of the aluminum alloy 6016 sheet.