Plastic gradient from geometrically necessary dislocation(GND) can affect material behavior significantly. In this research, mechanical behavior of polycrystalline solid is investigated using the finite element method incorporating plastic gradient from long range dislocation or GND effect. Plastic gradient effect is implemented in the analysis model by considering a long range strain term as well as elastic and plastic terms in the multiplicative decomposition. In the model, gradient hardness coefficient and length parameter are used to evaluate the effect of the long range strains and sensitive study is conducted for the parameters. It is confirmed that the GND amplifies hardening response of polycrystals compared with the single crystal.