The solution of boundary-value problem in mechanics of polycrystalline materials is represented as perturbation series upon intergrains interaction like it is made in quantum fields theory. The terms of this series satisfy the infinite set of integral equations that reduced to systems of linear algebraic equations if to neglect the heterogeneities of strains within domain of individual grain. The method allows to take into account a microstructure of materials. Few examples of method application are demonstrated.