On the equations of the analytical dynamics of the quasi-3 D plate theory of I . N . Vekua type and some their solutions

The plate theory of th order is constructed on the background of the Lagrangian variational formalism of analytical dynamics of continuum systems and the dimensional reduction approach of I. N. Vekua – A. A. Amosov. The plate model is defined within the configuration space, the set of field variables, and the Lagrangian density. The field variables are determined by the coefficients of the biorthogonal expansion of the spatial displacement vector field with respect to the dimensionless normal coordinate. The dynamic equations are derived as Lagrange equations of the second kind of the two-dimensional continuum system. The dynamics of the plane elastic layer is considered as an example, the normal wave propagation is described on the basis of refined plate theories of various orders, and the convergence of approximate solutions to the exact solution of the three-dimensional elastodynamics problem is analyzed for different wavenumbers.


Introduction
A plate is used nowadays as a mathematical model of many modern devices in machine industry.It is to be noted that in the high-frequency dynamics of composite plates the classical Kirchhoff's and even the refined first-order shear deformation plate theories fail whereas they are quietly consistent with the low-frequency dynamics of thin-walled structures [1][2][3][4][5][6].Moreover, the waveguide modeling requires refined plate theories; the authors of the paper [7] note that "more reliable 2D models are needed for high-frequency vibrations, wave propagation etc.".In general, the refinement of plate models consists in the accounting of supplementary degrees of freedom in addition to the translation and rotation of the middle surface point in the plate kinematics [1,[8][9][10][11].Such plate theories so-called "quasi-3D models" are used as well for highly anisotropic composite plates and functionally graded thin-walled structures [6,[12][13][14], as for the investigation of boundary layers in dynamics [15], boundary and edge waves [16][17][18], etc.The refined plate and shell theories can be also useful in problems of interaction of acoustic waves and thin-walled structures based on approximate diffraction models (e. g. see [19][20][21], and others).
Many methods of construction of refined plate and shell models can be used.The asymptotic integration approach [15] seems to be powerful and efficient method of the qualitative analysis of the plate and shell dynamics (for instance, see [22]).On the other hand, the asymptotic method does not allow one to construct the full hierarchy of solutions [23] approximating the three-dimensional solution in various norms [24].At the same time the formal series expansion of the displacement vector, stress tensor etc. offers some features of numerical algorithm construction if such a plate model is used together with finite element [24] or meshless numerical simulation [25].As well power series can be used [1,4,26,27] as special function expansions [12,13,28].One of the most powerful and universal approaches is based on generalized Fourier series [6,9,10,24,26,[29][30][31].
Here the higher-order plate theory based on the Lagrangian formalism of analytical mechanics of continua combined with the dimensional reduction approach [9] is used.The plate model interpreted as a two-dimensional continuum consists in the configuration space, the set of field variables being the biorthogonal expansion coefficients of the three-dimensional displacement vector field with respect to the thickness coordinate, and the Lagrangian density defined on the two-dimensional area corresponding to the plate middle surface (for more details, see [11,31,32]).This kind of theory so-called "elementary" allows one to satisfy the boundary conditions on the faces of a plate approximately as a result of the convergence of the sequence of two-dimensional solutions.To satisfy the boundary conditions exactly, for instance in the case of contact interaction, the mixed formulation [29,30] can be used, or the "extended" theory [9] or can be constructed.In particular, the boundary conditions shifted from the faces onto the middle surface become constraint equations, and the constrained variational problem is solved by the Lagrange multipliers method [33]; this approach allows one to obtain consistent low-order approximations [27] but seems to be a bit too complex when the order of the theory rises.
The well investigated problem of normal wave propagation in the plane elastic layer [34] can be used to analyse the properties of the constructed theories hierarchy and the convergence of the two-dimensional solutions [35][36][37][38][39][40].Here the normal waveforms corresponding to some specific wavenumbers are analysed; some results unpublished in the cited articles are presented.

Equations of the analytical dynamics of a plate as a two-dimensional Lagrangian continuum system
Let the plate be a three-dimensional elastic body: ⊂ , = ∪ , = ± ⊕ , with the faces ± and a lateral surface [11,[31][32][33]: where is the smooth base surface and 2ℎ is the plate thickness [11,31].
The mathematical model of a plate consists in the two-dimensional manifold , ̅ = ∪ ( = ∩ ) with the curvilinear chart , [11]: The covariant base vectors defined as = , ≡ / allow one to determine the metrics = ⋅ ; = ( × )/ √ = const is the normal unit vector, and = det .The linear dynamics problem statement for a plate can be based on the Hamilton principle [11,32]: The volumetric and surface Lagrangian densities can be written as follows [36]: where = + is the spatial distribution of the displacement vector field.The Lagrangian formalism of analytical dynamics allows one to introduce the configuration manifold Ω with the generalized coordinates ( ) : × → , ∈ ∪ {0} [11,31], so that = ( ( ) ), and the tangent fibration Ω can be defined with the base vectors ( ) = / ( ) .For linear systems such as the Eq. ( 1) the configuration manifold Ω becomes a Euclidian space; thus, the reduction of the three-dimensional elastodynamics problem given by Eq. ( 1) consists in EKATERINA L. KUZNETSOVA, ELENA L. KUZNETSOVA, LEV N. RABINSKIY, SERGEY I. ZHAVORONOK the projection of Ω onto its subspace Ω ( = 0,1 … ) [11,36,38].To construct a plate theory, the biorthogonal system ( ) ( ) , ( ) ( ) is used, therefore ( ) can be interpreted as field variables of the first kind [11]: ( ) = ( , ( ) ) .Here ( , ) is the scalar product [11,32], and the vector components ( , ) are supposed to be square integrable over [-1,1]∋ .The density of Lagrangian can be now defined on ̅ as follows [32,38]: ( Here ∇ denotes the covariant derivative on , and the following linear operators are defined (see [11]): The generalized stiffness for arbitrary anisotropic plates are defined as follows [11,32,38]: Here are contravariant components of the elastic constants tensor .Thus, the two-dimensional continuum system is defined on ̅ within the finite-dimensional configuration space Ω , + 1 field variable ( ) , and the Lagrangian densities , .The dynamic equations of the generalized plate theory of th order can be obtained as Lagrange equations of the second kind [31] for the two-dimensional continuum system given by the Eq. ( 2) (see also [32,38]): Their natural boundary conditions can be represented in the following notation [32,38]: Finally, the initial conditions are represented as follows: The initial-boundary value problem statement, Eq. (3-7), corresponds to the so-called "elementary" theory of plates of the th order; in other words, the boundary conditions on the faces ± are considered implicitly and can be satisfied only after solution's convergence at the point = ± 1.This theory allows one to obtain the simplest equations system and is preferable for higher orders.To satisfy the boundary conditions exactly we can formulate the "extended" theory of plates (for more details see [33]).

Modeling of the normal waves in the plane elastic layer on the background of the th order plate theory
Let us consider the plane problem for an elastic layer of thickness 2ℎ; here and below denotes shear modulus, denotes the Poisson ratio, and is the mass density.Let the normal waves be propagating along the axis .Finally, let us introduce the following dimensionless variables (see also [11,[35][36][37][38]): where = ⁄ ⁄ denotes the shear wave velocity.Thus, accounting for the material isotropy and homogeneity, for the homogeneous boundary conditions ± = 0 and absence of bulk forces ( ( ) = 0), and for the definition of the dilatation wave velocity = ( + 2 ) ⁄ ⁄ , we transform the general dynamic Eqs. ( 3), (4) to the following dimensionless formulation derived in [36,37]: The displacement field in the normal wave propagating along is defined as follows [34]: ̃ = ⁄ denotes the dimensionless phase velocity depending on the dimensionless wavenumber = ℎ, ( ) is the amplitude vector.Substituting Eq. ( 10) into Eq.( 9) we obtain the spectral problem defined by the Eq. ( 11) analogous to the one described in [35]: The eigenvalues of the matrix , Eq. ( 12), depending on form the spectra of longitudinal and bending waves.The corresponding dispersion curves for the phase velocities are shown on Figs. 1 and 2 (see also [11] and [35] for the dispersion curves corresponding to the dimensionless phase frequencies).It can be seen that the curves corresponding to the phase velocities of the lowest-order longitudinal and bending modes tend to the Rayleigh wave velocity as the wavenumber rises, → ∞.All dispersion curves except the 1st one have the characteristic cross points with the line corresponding to the dilatation velocity, ( ) = , > 0; the exact solution of the three-dimensional elastodynamics problem [34] gives the characteristic wavenumber * = ( − 1) / .The phase velocities computed on the basis of the spectral problem given by the Eqs.( 11), (12) for the plate theories of the order = 1, 2, 3, …, 10 are shown below in the Table 1.The same results for the phase frequencies are presented in [36].Analyzing it, we can conclude that the solution at the point * converges slowly that the at the point = 0 corresponding to the locking frequencies of normal waves (see [11,35]).
For the longitudinal waves ( = (2 − 1)/2) we have the following exact solution [34]: The normalized exact and approximate waveforms corresponding to these modes are shown on the Fig. 5.The relative mean square error of the waveform approximation given by the th order plate theory can be introduced as follows [37]: This error definition will be used below to estimate minimum orders of plate theory that allow one to approximate the exact solution [34].

Conclusions
Thus, the convergence of the approximation given by the sequence of solutions based on the plate theories of various orders depends significantly on the wavenumber especially in the phase velocity domain ̃> .The phase velocities converge slowly near the point = .It can be seen that five lowest phase velocities obtained on the background of the quasi-3D plate theory tends to the dilatation velocity at following plate model's orders (Table 2).The same situation can be observed after constructing the normal waveforms at different nonzero wavenumbers (for more results see [37]).
The convergence of the waveforms computed at the point = ( − 1) / (Goodier and Bishop modes) as well as at the point = (2 − 1)/2 (Lamé modes) on the groundwork of the quasi-3D plate theory of th order is estimated by the mean square error ∆ , Eq. ( 16).
Let us define the numerical convergence by the threshold level ∆ ≤ 0,05 (this depend on the problem, but this level seems to be satisfying in most applications).The minimum model orders allowing one to obtain such a convergence of the waveforms corresponding to the lowest propagating modes are shown in the Table 3.It can be also shown that the relative mean square error of the approximate normal modes become almost stable in the wavenumber domain ∈ [ , 15] where = is the phase frequency, and can be roughly estimated by the error ∆ ( * ), * = ( − 1) / [37].The obtained solution can be useful in the transient waveguide dynamics analysis based on the approximate models.It can be noted that these approximations give the upper estimate of the phase velocities of the normal waves due to the reduced number of degrees of freedom (field variables), but the obtained overestimation depends not only on the model order.
It is shown that the approximate satisfaction of the boundary conditions on the faces in terms of the "elementary" plate theory results in the overestimated longitudinal and bending stiffness especially for the lowest-order models.This drawback can be eliminated on the basis of the "extended" plate theories [33] that allow one to take into account the boundary conditions on the plate faces that are "lost" in the "elementary" theories and to account their effect on the tangent components of the stiffness tensors ( ) .

2841 .
ON THE EQUATIONS OF THE ANALYTICAL DYNAMICS OF THE QUASI-3D PLATE THEORY OF I. N. VEKUA TYPE AND SOME THEIR SOLUTIONS.

Table 2 .
Convergence of the 1, 2, …, 5th phase velocities given by the th order plate theory: