This paper confirms that in the u-p coupling model, the calculation results are consistent with the results of the analytical solution. This paper confirms the rationality and feasibility of the u-p-U model solving technology in the u-p-U coupling model. This paper confirms that the calculation results of the two forms are almost the same whenk=10-4m/s and k=10-5m/s. When k=10-3m/s the results of the two calculation forms are slightly different. The finite incremental calculus stabilization method can effectively eliminate the phenomenon of numerical oscillation in the region where the pore pressure gradient changes drastically. The finite incremental calculus stability method has a good stability effect in dealing with the dynamic equation of saturated soil which is difficult to calculate.