The micromechanical and macromechanical behavior of idealized granular assemblies, made by linearly elastic, frictionless, polydisperse spheres, are studied in a periodic, triaxial box geometry, using the discrete element method. Emphasis is put on the effect of polydispersity under purely isotropic loading and unloading, deviatoric (volume conserving), and uniaxial compression paths. We show that scaled pressure, coordination number and fraction of rattlers behave in a very similar fashion as functions of volume fraction, irrespective of the deformation path applied. Interestingly, they show a systematic dependence on the deformation mode and polydispersity via the respective jamming volume fraction. This confirms that the concept of a single jamming point has to be rephrased to a range of variable jamming points, dependent on microstructure and history of the sample, making the jamming volume fraction a state-variable. This behavior is confirmed when a simplified constitutive model involving structural anisotropy is calibrated using the purely isotropic and deviatoric simulations. The basic model parameters are found to depend on the polydispersity of the sample through the different jamming volume fractions. The predictive power of the calibrated model is checked by comparison with an independent test, namely uniaxial compression. The important features of the uniaxial experiment are captured and a qualitative prediction for the evolution of stress and fabric is shown involving a “softening” regime in both stress and fabric – stronger for the latter – that was not prescribed into the model a priori.