In this article, the authors numerically analyzed the aerostatic thrust bearing with porous restrictors using the finite-element method. After the existing methods were considered, a new analysis procedure is developed for improving the efficiency and application of the finite-element analysis program. In this procedure, HYPERMESH is used for preprocessing and postprocessing with the solver programed with C. The one-dimensional (1D)/three-dimensional (3D) source flow model derived from D'Arcy's law and the partial porous-bearing model are both supported by the solver. The matrix compress storage arithmetic and INTEL MATH KERNEL LIBRARY linear sparse solvers are used for improving the utilization ratio of memory and the solution efficiency. Therefore, the solver can be used to analyze the different applications of the porous restrictors in engineering. The program is validated by comparing the results obtained by different methods. The calculation efficiency and results of the 1D/3D source flow model are analyzed and compared with experimental results. The 1D model is shown to be suitable for engineering issues.