In the finite element numerical simulation, it is an urgent challenge to determine the ground motion input for the layered foundation under SV wave of oblique incidence over critical angle. Exact dynamic stiffness matrix in the frequency domain (i.e., frequency domain stiffness matrix method) of layered foundation is used to derive the formula for calculating the equivalent node force of ground motion input under arbitrary angular oblique incidence of SV wave, and the validity and accuracy of the ground motion input applied to ABAQUS numerical simulation are demonstrated by simulating seismic wave field of homogeneous half-space and layered foundation under oblique incidence of SV wave through ABAQUS software. It is shown that the inclined SV wave incidence with an arbitrary angle can be realized through ABAQUS simulation by means of stiffness matrix method, and the method is of high accuracy, especially in the case of SV wave incidence with an inclined angle greater than the critical one. The oval like motion trails of surface particles in uniform half space and the wave motion characteristics of layered half space are simulated commendably by finite element numerical modelling. On the above basis, the frequency domain stiffness matrix method is further combined with equivalent linearization method, which solves the ground motion input problem of two-dimensional layered half space considering soil nonlinearity for the inclined incidence with arbitrary angle.