[en] This article presents an analytic method and calculations of strong motion spectra for the energy, displacement, velocity and acceleration based on the physical and geometric ground properties at a site. Although earthquakes occur with large deformations and high stress intensities which necessarily lead to nonlinear phenomena, most analytical efforts to date have been based on linear analyses in engineering seismology and soil dynamics. There are, however, a wealth of problems such as the shifts in frequency, dispersion due to the amplitude, the generation of harmonics, removal of resonance infinities, which cannot be accounted for by a linear theory. In the study, the stress-strain law for soil is taken as tau=G₀γ+G₁γ³+etaγ where tau is the stress, γ is the strain, G₀ and G₁ are the elasticity coefficients and eta is the damping and are different in each layer. The above stress-strain law describes soils with hysterisis where the hysterisis loops for various amplitudes of the strain are no longer concentric ellipses as for linear relations but are oval shapes rotated with respect to each other similar to the materials with the Osgood-Ramberg law. It is observed that even slight nonlinearities may drastically alter the various response spectra from that given by linear analysis. In fact, primary waves cause resonance conditions such that secondary waves are generated. As a result, a weak energy transfer from the primary to the secondary waves takes place, thus altering the wave spectrum. The mathematical technique that is utilized for the solution of the nonlinear equation is a special perturbation method as an extension of Poincare's procedure. The method considers shifts in the frequencies which are determined by the boundedness of the energy