Abstract: The present paper uses the mathematics model for consolidation of unsaturated soil developed in ref.[I] to solve boundary value problems.The analytical solutions for one-dimensional consolidation problem are gained by making use of Laplace transform and finite Fourier transform.The displacement and the pore water pressure as well as the pore gas pressure are found from governing equations simultaneously.The theoretical formulae of coefficient and degree of consolidation are also given in the paper.With the help of the method of Galerkin Weighted Residuals, the finite element equations for two-dimensional consolidation problem are derived.A FORTRAN program named CSU8 using 8-node isoparameter element is designed.A plane strain consolidation problem is solved using the program, and some distinguishing features on consolidation of unsaturated soil and certain peculiarities on numerical analysis are revealed.These achievements make it convenient to apply the theory proposed by the author in engineering practice.

Key words: pore water pressure, pore gas pressure, coefficient of consolidation, degree of consolidation, initial settlement, finite element