Článek se zabývá tvorbou numerických metod popisujících chování zemin jako třífázového prostředí. Z praktického hlediska je problém redukován na řesení nenasycených a plně nasycených zemin. Při řešení je akceptována Fillunger-Terzaghiho koncepce efektivních napětí. K popisu plastického chování skeletu zeminy jsou představeny modely Mohr-Coulomb, Drucker-Prager, CAM clay model a model kritického stavu. Proudění tekutiny se řídí Darcyho zákonem, v článku je uvedena experimentálně ověřená proměna součinitele permeability. Časová diskretizace je prováděna zobecnělou metodou středního bodu a prostorová diskretizace je vedena v duchu metody konečných prvků. Pro úspěšnost řešení je obecně zvolena přírustková formulace sdruženého problému