Bending and Vibrations of a Thick Plate with Consideration of Bimoments ()
Received 1 July 2016; accepted 27 August 2016; published 30 August 2016

1. Introduction
Theory of plates and shells has a special place in design of structural elements. Specified theories of plates are built by many authors. All existing specified theories of plates are developed on the basis of a number of simplifying hypotheses. An overview of the main statements and common methods of constructing an improved theory of plates and shells can be found in the works of S. A. Ambartsumyan [1] , K. Z. Galimov [2] , Sh. K. Galimov [3] , Kh. M. Mushtari [4] and others. Static problem of the bending of a thick isotropic plate in three- dimensional theory of elasticity is considered by B.F. Vlasov in [5] , which gives an exact analytical solution in trigonometric series. Monograph by E.N. Baida [6] is devoted to solving the problem of bending of orthotropic plates in trigonometric series. Numerical results of displacements and stresses are obtained.
The authors in [7] - [10] deal with dynamic problems of plates with anisotropic properties. Karamooz Ravari M.R., Forouzan M.R. [7] have studied the problems of plates oscillations. Frequency equations of orthotropic circular ring plate were obtained for general boundary conditions in oscillation plane. In [8] the solution of transition oscillations of rectangular viscous-elastic orthotropic plate are given for concrete strain models according to Flugge and Timoshenko-Mindlin’s theories. The paper [9] is devoted to analytical solution of the problem of forced steady-state vibrations of orthotropic plate. By the method of superposition the problem is reduced to a quasi-regular infinite system of linear equations. In [10] an analytical method of solution of spatial problem of bending of orthotropic elastic plates subjected to external loads on upper and lower edges is developed. In [11] a problem is considered of a bending of orthotropic rectangular plate laying on two-parameter elastic foundation. Research in the field of thick plates has shown that in the case of spatial deformation of a plate along its thickness there occurs the nonlinear laws of displacements distribution and the hypothesis of plane sections is violated. In the cross-sections of the plate except for the tensile and shear forces, bending and torsional moments, there appear the additional force factors, the so-called bimoments. The author of the article addresses the problem of bending and vibrations of thick plates based on bimoment theory of plates [12] - [15] , built as a part of three-dimensional theory of elasticity, using the method of displacements decomposition in one of the spatial coordinates in Maclaurin infinite series.
This paper gives a brief description of the technique of constructing a theory of plates with consideration of bimoments generated due to displacements distribution of cross-section points by a non-linear law. Here the equations of bimoments are built with the equation of three-dimensional dynamic theory of elasticity, described on face surfaces of the plate. The bimoments are introduced in stress dimensions and are characterized by the intensity of generated bimoments. We would use the designations and determinant correlations of forces, moments, bimoments and equations of motion relative to these force factors.
Unlike bimoment theory in [14] and [15] , here the bimoment equations are built with the equation of three- dimensional dynamic theory of elasticity, described on face surfaces of the plate. Bimoments are introduced in stress dimensions, and they characterize the intensity of generated bimoments.
Determinant relationships of forces, moments, bimoments and equations of motion relative to these force factors are given.
2. Statement of the Problem
Consider an orthotropic thick plate of constant thickness
and dimensions
in plane. Introduce the designations:
―elasticity moduli;
―shear moduli;
―Poisson ratio of plate material.
When building an equation of motion the plate is considered as a three-dimensional body and all components of stress and strain tensors:
are taken into consideration. The components of displacement vector are the functions of three spatial coordinates and time
.
The components of strain tensor
are determined from Cauchy relation as:
(1.а)
(1.b)
For orthotropic plate, the Hooke’ law, in a general case, is written as:
(2.a)
(2.b)
(2.c)
(2.d)
where
are the elastic constants, determined through Poisson ratio and the moduli of elasticity in the form [14] [15] .
As an equation of motion of a plate we would use three-dimensional equations of dynamic theory of elasticity:
(3.а)
(3.b)
(3.c)
where
is a density of plate material.
Boundary conditions on lower and upper face surfaces of the plate
and
are:
(4.a)
(4.b)
Here
and
are distributed external loads, applied to upper and lower face surfaces of the plate
and
along the direction of
coordinates axes.
3. Method of Solution
The methods of building the bimoment theory of plates are based on Cauchy relation (1), generalized Hooke’s law (2), three-dimensional equations of the theory of elasticity (3), boundary conditions on face surfaces (4). A proposed bimoment theory of plates is also described by two non-connected problems, each of which is formulated on the basis of six two-dimensional equations of motion with corresponding boundary conditions.
The components of displacement vector are expanded into Maclaurin infinite series in the form:
(5.а)
(5.b)
Here
are unknown functions of two spatial coordinates and time
,
. In a general case, these functions are determined according to the formulae:
![]()
The displacements in stresses in upper
and lower points
in plate fibers we would designate as
and
.
The first problem of bimoment theory describes tension-compression and transverse reduction of the plate, and the second one―the bending and transverse shear of the plate. Determinant relationships and corresponding equations of motion of the plate in the first and second problems are briefly described below.
The first problem is described by the forces and bimoments with six generalized functions
, which are determined by relationships:
(6)
(7)
Introduce the external loads for the first problem
(8)
The expressions of longitudinal and tangential forces are written as [12] - [15] :
,
(9.а)
(9.b)
The intensities of the bimoments
from tangential stresses
have the expressions
(10.а)
The intensity of the bimoment
from normal stress
is written in the form:
(10.b)
The equations of motion relative to longitudinal and tangential forces and bimoments from tangential and normal stresses have the form [12] - [15] :
(11)
(12)
Note, that the expressions of force factors (9), (10), and hence, the equations of motion of the system (11), (12) is rigorously built. This system consists of three equations relative to six unknown functions
. As could be seen, three equations are missed. If in expressions (9.а) the terms
are omitted, then we would obtain two equations of motion of classic theory of plates in the form (11), since the equation of motion (12) becomes isolated and fail.
The second problem of bimoment theory consists of the equations for bending moments, torsional moments, shear forces relative to six kinematic functions
, determined by formulae:
(13)
(14)
Introduce the generalized external loads for the second problem
(15)
Bending, torsional moments and shear forces, which are rigorously built, have the form [12] - [15] :
(16.а)
(16.b)
(16.c)
(16.d)
The system of equations of motion of the second problem consists of two equations relative to bending, torsional moments and one equation relative to shear force and it is written in the form [12] - [15] :
(17)
(18)
Note, that the expressions of forces and moments (16), hence, the equations of motion of the system (17), (18) are rigorously built. Similar to the first problem, here three equations are missed. The system of equations of motion (17), (18) consists of three equations relative to six unknown functions
. If in expressions of forces and moments, into the equations of motion (17) and (18) conventionally introduce
and
, and the shear modulus
substitute for
, (where
is a shear coefficient), then an equation of motion of plates could be obtained according to Timoshenko’s theory.
To complete the systems (11), (12) and (17) and (18) it is necessary to build two more systems, with three equations in each. Write down three equations of motion of the theory of elasticity (3) on face surfaces of the plate
and
. Adding and subtracting the equations of the theory of elasticity (3) on face surfaces of the plate
and
, and taking into account the Hooke’s law (2), surface conditions (4) and designations (6), (7) and (13), (14), two independent systems with three equations in each could be obtained. The first of these systems describes the first problem and has the form:
(19)
(20)
Here the intensities of the bimoments
―under transverse reduction and tension-compression of the plate, generated due to
are:
(21)
are the intensities of the bimoments generated due to transverse stresses
:
(22)
The second system of equations obtained from the equations of the theory of elasticity (3) is written in the form:
(23)
(24)
Here
are the intensities of the bimoments under transverse bending and shear for the second problem generated due to the stresses
:
(25)
The intensities of the bimoments
, generated due to the stresses
, under transverse shear and bending are written in the form:
(26)
The intensities of the bimoments
are determined from Hooke’s law (2) with consideration of the conditions on face surfaces
and
(4) as:
(27)
(28)
Here
.
The expressions of the intensities of the bimoments
are determined by the solution of the system of linear algebraic equations relative to coefficients of Maclaurin series
, which are obtained by the substitution of the series (5) into the conditions on face surfaces at
and
(4) and designations (6), (7).
(29.а)
(29.b)
(30)
The expressions of the intensities of the bimoments
are determined by the solution of the system of linear algebraic equations relative to coefficients of Maclaurin series
, which are obtained by the substitution of the series (5) into the conditions on the face surfaces at
and
(4) and designations (13), (14).
(31.а)
(31.b)
(32)
Write down the formulae to determine the displacements on the face surfaces of the plate
and
:
(33)
Formulae for stresses on the face surfaces of the plate
and
have the form:
(34)
Maximum values of displacements and stresses of the plate are reached on the face surfaces of the plate and are determined by the solutions of the first and second problems by the formulae (33) and (34).
Note, that the expressions of intensities of the bimoments (10), (27), (28), (29), (30), (31) and (32) are built for the first time and are new in the theory of plates.
Consider the boundary conditions of a discussed problem for the thick plates.
1) On the border of the plate the displacements are zero. On the edges of the plate
and
the conditions should be as follows:
(35)
(36)
2) On the border
the plate is supported. The following conditions should be satisfied:
(37)
(38)
3) On the border
the plate is free of supports. The following conditions should be satisfied
(39)
(40)
Boundary conditions on the border
are similarly written.
When studying the problem of transverse bending and shear it is enough to consider only the second problem with the equations of motion (17), (18), (23), (24) and boundary conditions (35)-(40).
4. Solution of Tests Problem
As an example, consider the forced harmonic vibrations of a cantilever rectangular plate fixed on both ends under the effect of harmonic periodic external load:
(41)
where
is an amplitude, frequency and the mode of vibration of an external load, respectively. Note, that if
, we obtain the problem of static bending of the plate.
Substituting (41) into (8) and (15) determine the load terms of the equation of motion. For a plate fixed on both ends the boundary conditions are written in the form (35) and (36).
5. Numeric Results
First determine eigenfrequencies of the plate. After dividing the variables by spatial coordinates and time, the problem is solved by Finite Difference Method. The step in spatial coordinates is
. In calculations, for isotropic plates
are given as an initial data.
For square plates with dimensions
the value of eigenfrequency is
. With increasing dimensions of the plate up to
the value of eigenfrequency is
. For square plates with dimensions
the value of eigenfrequency is
.
Table 1 shows the results obtained for the displacements, moments and forces in fixed square plates
under different values of dimensionless frequency
. When the value of the frequency
of external effect
approaches the eigenfrequency
the values of the displacements, forces and moments dramatically increase; this testifies of gradual transition of the motion of plate points into resonant mode. As seen, an abrupt increase in the values of displacements, forces and moments could be observed.
Table 2 and Table 3 show numeric values of displacements, moments and forces, calculated for the fixed square plates with dimensions
and
, respectively, for different values of dimensionless frequency
.
Calculations show that when the value of the frequency of external effect
approaches eigenfrequency, an increase in the values of displacements, forces and moments is observed; this testifies of gradual transition of the motion of plate points into resonant mode.
![]()
Table 1. Displacements, forces and moments at
.
![]()
Table 2. Displacements, forces and moments at
.
![]()
Table 3. Displacements, forces and moments at
.
Calculations show that the equations of motion of a plate (23) may be substituted by kinematic conditions relative to tangential stresses:
(26)
Kinematic equations serve to determine the generalized displacements
.
The equations (26) are determined by the solution of the system of linear algebraic equations relative to coefficients of the series (5)
, which are obtained by the substitution of the series (5) into the conditions on the face surfaces at
and
(4) and designations (13), (14).
6. Conclusion
Based on these studies, we would note that using the method of expansion in a series as part of three-dimen- sional dynamic theory of elasticity, a two-dimensional bimoment theory of orthotropic thick plates was developed and the equations of motion of the plate relative tot forces, moments and bimoments were built. It is shown that the problem in the general case is reduced to the definition of twelve unknown functions of two spatial coordinates and time. New expressions to determine the forces, moments and bimoments of the plates were built, as well as the methods for solving the problems of free and forced vibrations of plates based on Finite Difference Method.