Earthquake preparation and occurrence is a complex physical process. The use of numerical analysis, which constructs three-dimensional dynamic models of the crust and upper mantle to simulate the rock deformation process, is currently one of the most effective methods to study the crustal energy transfer and accumulation. The simulation result of current crustal deformation is consistent with the existing GPS data around the Eastern Himalayan Syntaxis and its surrounding areas. The numerical simulation results of deformation process in the Capital area from 1989 to 1998 clearly show that the high strain energy concentration region shifted from Datong area where 1989 earthquake (MS 5.8) occurred to Zhangbei area where 1998 earthquake happened.