Abaqus - Kirsch equation for stress around a hole?

Hello,

I have defined a basic 2D rectangular model with a hole in the center and I am trying to apply initial geostatic stress to this model and see if the results from Abaqus are same as the analytical solution of Kirsch equation. When I use "predefined field>initial>stress" to define geostatic stress (image 35, below), which is the method used in ABAQUS sample example file "hydrfract_xfem_c3d8rp" to define geostatic stress, the results are not same as the Kirsch equation but when I define the stresses as "load>pressure" (image 36), the results are same with the Kirsch equation. But my problem is that I have seen in some videos that I shouldn't use pressure load in geostatic step and also when I use pressure to define geostatic stress, I can not run the *cflow in soils step and it gives errors, while if I use predefined stress, I don't get any errors. How can I define geostatic stress properly that the stress around the hole can agree with the Kirsch equation? And there is another option in "predefined field>geostatic stress" that I read the user assistance but I don't know how to use it. Any help is appreciated

Thank you