In the field of geotechnical engineering, it is of great significance to analyze the law and mechanism of soil hydraulic fracturing, which can provide a basis for the study of the diffusion and distribution of splitting grouting veins and guide practical engineering. In this paper, a three-dimensional finite element program is developed by using block diagonal preprocessing, PSQMR iterative method and improved sparse matrix vector multiplication parallel algorithm, which can realize large-scale Biot consolidation finite element calculation. The finite element solution based on CPU serial computing platform and GPU parallel computing platform is realized. The calculation scale is greatly improved, and large scale finite element calculation can be carried out on personal computer. By comparing the numerical simulation results, experimental simulation results and theoretical analysis with the calculation results, the correctness of the calculation method is verified, which provides a numerical simulation tool for the study of hydraulic fracturing in soil.