Abstract：Multigrid method solver is of high numerical efficiency when used in solving linear equations derived from boundary-value problems of partial differential equations. There are some shortages in geometrical multigrid method which restricts its application area. The algebraic multigrid method is used to solve finite element linear equations which are derived from three-dimensional finite element analysis of rock mechanics and engineering. The three-dimensional coarse-grid selection method based on element agglomeration and the three-dimensional interpolation operator are briefly introduced. By using the newly developed three-dimensional finite element program based on algebraic multigrid method，four different numerical experiments are designed，and carried out to validate its convergence character，numerical efficiency and practical application to modeling excavation problem of rock engineering. The numerical experiments show that the algebraic multigrid method is of better stability，good convergence character and better adaptability，with much higher numerical efficiency with increasing of the number of linear equations and much less computer memory compared with direct method. Increasing. The algebraic multigrid method has much better numerical efficiency. The algebraic multigrid method is suitable and efficient for three-dimensional finite element modelling of large-scale geomechanical engineering.
