In this paper, the process of simulating the occurrence of an earthquake on the fault plane for the 2016 Kumamoto, Japan earthquake is investigated dynamically using finite element methods. The two-dimensional finite element method is used by concentrating friction at the finite points. The used element is the elastic element, whose properties are determined based on the shear wave velocity and longitudinal wave. In order to verify the proposed method, the acceleration and velocity response spectra, Fourier spectrum, and other main parameters of the record obtained from the simulation were compared with the Kik-net data. From the obtained results, the role of displacement and velocity of its application to the fault plate in PGA can be pointed out, so that with a faster displacement of the fault plane, a larger PGA is expected. Friction and shear stiffness between the two fault planes and the dependence of friction on velocity are the main and influential factors in the characteristics of the produced earthquake. In high-frequency earthquakes, friction is dependent on velocity, and when friction is not dependent on velocity, static plays a major role, and no drastic changes in the response spectrum are observed. [ABSTRACT FROM AUTHOR]