Based on one-dimensional seepage and deformation governing equations, PDE software is used to analyze the coupled governing equations. The method is available to arbitrary functions of SWCC (soil-water characteristic curve), and it considers that both the coefficient of permeability at saturation and the porosity change with soil strain. Compared with analytical solution, the numerical solution shows high precision and it can effectively solve coupling problems. It is found that coupling effect of seepage and deformation in unsaturated soils plays an important role in the pore-water pressure profiles, and that the coupling effect should be considered during rainfall infiltration. At early stage of rainfall infiltration, settlement obviously increases with time, and the pore-water pressure considering the coupling moves more slowly than that without considering coupling due to positive H. The settlement is related with initial pore-water pressure profiles and H value. The coefficient of permeability at saturation and the porosity have an effect on rainfall infiltration and steady-state seepage, but their change has a weak influence on deformation in unsaturated soils..