An analytical chemo-hydromechanical model of residual soils was presented in this study. The modeling procedure utilized a unique phenomenological approach, which highlighted some of the macro-mechanical influences responsible for the behavior of residual soils. The macromechanical analysis yielded an extended poroelastic theory. The developed model was applied in solving a typical civil engineering problem of soil consolidation. For a chemical loading treatment, the problem was using requisite boundary conditions. Thereafter, the mathematical software Mathematica was used for solving the ensuing diffusion equations using physico-chemical properties typical of a residual soil sample. Even though the chemical loading was instantaneous, it took some time before its effect could be felt everywhere according to the poroelastic theory. From obtained results the evolution of the generated pore pressure showed increasing segments of the soil returning to their initial state after the passage of the pore pressure front. The vertical displacement showed a slight increase or elevation of the soil surface. The flux of water flow in the soil was initially positive before turning negative, which can be explained by initial successful penetration of the infiltrating liquid until it met increasingly tortuous paths as result of lower permeability. The volume (per unit area) of water leaving the soil after chemical loading demonstrated that a portion of the water flux associated with the infiltrating liquid ends up leaving the soil through the sides. The developed model was later compared with a similar porochemoelastic model found in literature with reasonable agreement.