A constitutive model for a porous metal subjected to general three-dimensional finite deformations is presented. The model takes into account the evolution of porosity and the development of anisotropy due to changes in the shape and the orientation of the voids during deformation. The pores are initially ellipsoidal and distributed randomly in an elastic-plastic matrix (metal). Under finite plastic deformation, the voids are assumed to remain ellipsoids, and to change their volume, shape, and orientation. At every point in the homogenized continuum, a “representative” ellipsoid is considered with principal axes defined by the unit vectors n⁽¹⁾,  n⁽²⁾,  n⁽³⁾  = n⁽¹⁾  x n⁽²⁾  and corresponding principal lengths a, b and c. The homogenized continuum is locally orthotropic, with the local axes of orthotropy coinciding with the principal axes of the representative ellipsoid. The basic “internal variables” characterizing the state of the microstructure at every point in the homogenized continuum are given by the local equivalent plastic strain , the local void volume fraction or porosity , the two aspect ratios of the local representative ellipsoid (w₁=c/a and w₂=c/b) and the orientation of the principal axes of the ellipsoid (n⁽¹⁾, n⁽²⁾, n⁽³⁾). A methodology for the numerical integration of the elastoplastic constitutive model is developed. The problem of ductile fracture near the tip of a blunt crack is studied by using the finite element method and comparisons with traditional constitutive models that assume isotropic behavior are made. A “plastic strain-gradient” version of the model will be also presented and issues associated with its numerical implementation will be discussed.