The present research studies the numerical modeling of the mechanical behavior of residual soils, by a micro-mechanical analysis, using as modeling tool the discrete elements method (DEM), which allows the representation of the soil as 2D particle assembly of circular discrete elements. There are several micro and macro structural phenomena present in residual soils, which determine its mechanical response: double structure and double porosity of soil mass, the presence of a clay-like matrix, partially saturated condition, grain cementation, grain crushing, and the presence of small blocks of not-weathered parental rock. Based on the study of these phenomena, mathematical formulations have been proposed in order to describe them, which have been implemented in the DEMlib discrete elements code. After calibration and validation of the implemented analysis program, the influence of each phenomena in soil response was studied by simulating edometrical and direct shear tests in discrete samples with several densities, criticizing the results and getting conclusions. Afterwards, these phenomena were studied in a combined way in order to reproduce the response of some young and mature residual soils. To finish up this work the conclusions obtained have been summarized and suggestions for deepening and further developing of the implemented research line have been stated.