In this paper, analytical solutions for one-dimensional nonlinear steady-state vertical flux through unsaturated soils are presented. The analytical solutions are applicable for homogeneous soils for which the hydraulic conductivity function is a generalized form of the well-known Brooks-Corey and rational models. The adopted generalized soil–water retention model involves two additional parameters which reflect the curvature of hydraulic functions near saturation. The main advantage of the used generalized hydraulic conductivity function is that it fits better experimental data. The soil domain is a finite-depth medium overlying a fixed groundwater level. The derived analytical solutions are obtained using the Maclaurin series expansion technique. Analytical solutions for both infiltration and evaporation are developed. For evaporation from a fixed groundwater level, the derived solutions can be used to predict the existence of a vaporization plane below the soil surface. For this case, analytical expression of the position of drying front separating liquid and gas regions is derived. The analytical results are compared to numerical ones for various infiltration and evaporation examples and excellent agreements are obtained. The generalized model is also used to fit experimental data obtained for two soils. The results show that the use of the standard Brooks-Corey model may lead to an overestimation of the position of the drying front.