Statistics of Using local second gradient model and shear strain localisation to model the excavation damaged zone in unsaturated claystone