distinct target of this paper is development of 'implicit discrete element method (dem)' in three dimensions for estimating the ultimate bearing capacity of shallow foundations on layered (non-homogeneous) soils consisting of sand overlying sand or clay overlying sand and vice versa with different shear resisting parameters. in this method, shear resisting parameters and density of ea...