Time history analysis is performed on geotechnical engineering structures under earthquake actions using 3D nonlinear FEM. The distribution of internal force on slide surfaces has been interpolated from the stress field of FEM with 3D Multi-grid Method. The safety factor with time history of sliding block can be obtained by limit equilibrium analysis and its dynamic stability can be evaluated. For blocks sliced by single slide surface and wedge-shaped double slide surfaces, general formulas of safety factor are given and their applicability have been analyzed. Analysis and evaluation of slide block aseismic stability are performed on high slope on the right bank of Dagangshan arch dam.