For the soil is a very complex natural material, significant strain gradient effect exist in soil analysis. Based on the "gradient" phenomenon, we add the plastic strain gradient hardening item into the traditional Cambridge yield surface. By using the consistency conditions and associated flow rule, we get the explicit expression of plastic strain gradient stiffness matrix. And the finite element method of plastic strain gradient is also shown in this article. Plastic strain gradient is actually a phenomenological non-local model containing microstructure information of the material. It may overcome the difficulties in simulating the gradient phenomenon by traditional mechanical model.