在岩土工程领域，分析土体水力劈裂规律和机理具有重要意义。通过使用块对角预处理、PSQMR迭代法及改进的稀疏矩阵矢量乘并行算法，开发土体中水力劈裂三维有限元程序，实现了基于CPU串行计算和GPU并行计算平台的有限元求解，大幅提高计算规模，在个人计算机上可进行工程尺度的大规模有限元计算。通过数值模拟结果、试验模拟结果和理论分析与计算结果进行的对比，验证了所述计算方法的正确性，为研究土体中水力劈裂提供了有力的数值模拟工具。

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.