A fundamental problem in mechanics of materials is the computation of the macroscopic response of polycrystalline aggregates from the properties of their constituent single-crystal grains and the microstructure. In this work, the nonlinear homogenization method of deBotton and Ponte Castañeda (1995) was used to compute “variational” self-consistent estimates for the effective behavior of different types of viscoplastic polycrystals, including a two-dimensional model, as well as cubic and hexagonal polycrystals. In contrast with the “incremental” and “tangent” self-consistent estimates, the new results are found to satisfy all known bounds, even in the strongly nonlinear, rate-insensitive limit. The new results also exhibit a more realistic ...