
,
which forms the boundaries in the space of system parameters within the first zone of the parametric instability.
 From the physical viewpoint, one can see that the parametric excitation of bending waves appears as a degenerated case of nonlinear wave interactions. It means that the study of resonant properties in nonlinear elastic systems is of primary importance to understand the nature of dynamical instability, even considering free nonlinear oscillations.
 Normal forms
 The linear subset of eqs. (0) describes a superposition of harmonic waves characterized by the dispersion relation
   
,
where 
  
 refer the 
 
 branches of the natural frequencies depending upon wave vectors 
 
. The spectrum of the wave vectors and the eigenfrequencies can be both continuous and discrete one that finally depends upon the boundary and initial conditions of the problem. The normalization of the first order, through a special invertible linear transform
  
leads to the following linearly uncoupled equations
   
,
where the 
  
 matrix 
 
 is composed by 
 
-dimensional polarization eigenvectors 
 
 defined by the characteristic equation
  
;
  
 is the 
 
 diagonal matrix of differential operators with eigenvalues 
 
; 
 
 and 
 
 are reverse matrices.
The linearly uncoupled equations can be rewritten in an equivalent matrix form [5]
 (12) 
  
 and 
 
,
using the complex variables 
  
. Here 
 
 is the 
 
 unity matrix. Here 
 
 is the 
 
-dimensional vector of nonlinear terms analytical at the origin 
 
. So, this can be presented as a series in 
 
, i. e.
 
, where 
  
 are the vectors of homogeneous polynomials of degree 
 
, e. g.
 
 Here 
  
and 
 
 are some given differential operators. Together with the system (12), we consider the corresponding linearized subset
(13) 
  
 and 
 
,
whose analytical solutions can be written immediately as a superposition of harmonic waves
   
,
where 
  
 are constant complex amplitudes; 
 
 is the number of normal waves of the 
 
-th type, so that 
 
 (for instance, if the operator 
 
 is a polynomial, then 
 
, where 
 
 is a scalar, 
 
 is a constant vector, 
 
 is some differentiable function. For more detail see [6]).
A question is following. What is the difference between these two systems, or in other words, how the small nonlinearity is effective?
 According to a method of normal forms (see for example [7,8]), we look for a solution to eqs. (12) in the form of a quasi-automorphism, i. e.
 (14) 
 
 where 
  
 denotes an unknown 
 
-dimensional vector function, whose components 
 
 can be represented as formal power series in 
 
, i. e. a quasi-bilinear form:
(15) 
 
, for example
 
 
 where 
  
 and 
 
 are unknown coefficients which have to be determined.
By substituting the transform (14) into eqs. (12), we obtain the following partial differential equations to define 
  
:
(16) 
 
. It is obvious that the eigenvalues of the operator 
  
 acting on the polynomial components of 
 
 (i. e. 
 
) are the linear integer-valued combinational values of the operator 
 
 given at various arguments of the wave vector 
 
.
In the lowest-order approximation in 
  
 eqs. (16) read
 
. The polynomial components of 
  
 are associated with their eigenvalues 
 
, i. e. 
 
, where
 
  
or 
  
,
while 
  
 in the lower-order approximation in 
 
.