Deformation quantization of field models: iterated variations in the star-products