How can limit the deformation in initial stress condition-inistate

[jiligeo]

Dear friends, I am modeling a soil mass using ANSYS classic v.11.
EDP material model were used for soil.
The initial stresses were defined by inistate command.
My problem is: the max deformation of model is too much (54 cm) and it is not true (the real value is zero or every small one).
The calculated values of initial stresses are true.
Is it possible to limit the deformations in initial stress condition? Or set the deformation of model to zero before to start the next load step?

Please find my APDL code for this calculation as follows.

Regards
Jalil

finish
/clear
/vup,1,-z
/view,1,1,1,+1
/plopts,logo,on
/triad,on
/efacet,4
seltol,1e-9
/title,Initial stress in soil mass

! ************************************************
! * main parameters *
! ************************************************
gg=9.81
*afun,deg

! material - m KPa ton/m3
dpEs1=40e3 $dpNus1=0.35 $dpRhos1=2.00
dpCs1=10 $dpPhis1=25 $dppsi1=10
dpkapa1dp=6*dpCs1*cos(dpPhis1)/(sqrt(3)*(3-sin(dpPhis1)))
dpalfah1dp=2*sin(dpPhis1)/(sqrt(3)*(3-sin(dpPhis1)))
dpalfap1dp=2*sin(dpPsi1)/(sqrt(3)*(3-sin(dpPsi1)))

dpalfah1edp=3*sqrt(3)*dpalfah1dp
dpkapa1edp=sqrt(3)*dpkapa1dp
dpalfap1edp=3*sqrt(3)*dpalfap1dp

/prep7
! ************************************************
! * element table *
! ************************************************
et,1,solid186 !3-D 20-Node Structural Solid

block,0,0.4,0,0.4,0,60

lesize, all, 0.2
vatt,1,1,1,0
MSHAPE,0,3D
MSHKEY,1
MSHMID,0
vmesh,all
uimp,1,ex,dens,nuxy,dpEs1,dpRhos1,dpNus1
tb,EDP,1, , ,lyfun
tbdata,1,dpalfah1edp, dpkapa1edp
!tbdata,1,1., 21
tb,EDP,1, , ,lfpot
tbdata,1,dpalfap1edp, dpkapa1edp
!tbdata,1,1., 21
!**************************
/solu
antype,static,new
allsel,all
inistate,write,1,,,,
acel,0,0,-gg

nsel,s,loc,x,0 $ d,all,ux,0
nsel,s,loc,x,0.4 $ d,all,ux,0
nsel,s,loc,y,0 $ d,all,uy,0
nsel,s,loc,y,0.4 $ d,all,uy,0
nsel,s,loc,z,60 $ d,all,uz,0

allsel,
nropt,
solcontorl,on,on
nlgeom,on
nsubst,10,10,2
neqit,30
time,1
solve
lsclear,all
finish

/solu

allsel,

inistate,read,,ist

nsel,s,loc,x,0 $ d,all,ux,0
nsel,s,loc,x,0.4 $ d,all,ux,0
nsel,s,loc,y,0 $ d,all,uy,0
nsel,s,loc,y,0.4 $ d,all,uy,0
nsel,s,loc,z,60 $ d,all,uz,0
!*********************
acel,0,0,-gg
! ***************************
allsel
nropt,
nlgeom,on
solcontorl,on,on
nsubst,10,10,2
neqit,30
time,1
solve
_________________