a b s t r a c tIn this work, a new and simple numerical approach to simulate nonlinear wave propagation in purely hysteretic elastic solids is presented. Conversely to classical time discretization method, which fully integrates the nonlinear equation of motion, this method utilizes a first-order approximation of the nonlinear strain in order to separate linear and nonlinear contributions. The problem for the nonlinear displacements is then posed as a linear one in which the solid is enforced with nonlinear forces derived from the linear strain. In this manner, a frequency analysis can be easily conducted, leading directly to a well-known frequency spectrum for the nonlinear strain. A mesoscale approach known as Preisach-Mayergoyz space (PM space) is used for the chacterization of the nonlinear elastic region of the solid. A meshless element free Galerkin method is implemented for the discretized equations of motion. Nevertheless, a mesh-based method can be still used as well without loss of generality. Results are presented for bidimensional isotropic plates both in plane stress and in plane strain subjected to harmonic monotone excitation.