Abstract We present a numerical scheme for obtaining guaranteed (reliable) and arbitrarily close two sided bounds to effective (homogenized) parameters of the linear elasticity problem. For upper bounds, we use standard finite element (FE) discretization so‐called primal problem with preconditioning based on fast discrete Fourier transformation (FFT). lower dual formulation some smoother FE app...