-□ □-
Розглянуто тривимiрну задачу про напружено-деформований стан плити, у яког характеристики пружностi матер^ алу змтюються за товщиною та опису-ються довшьними ттегровними функщя-ми. Показано, що точт, в сенЫ Сен-Венана, аналтичн розв'язки крайовог задачi можна отримати, якщо розподш навантажень задовольняе двовимiрному полкармошч-ному рiвнянню. Сформульовано наближе-н теори для iнженерних розрахун^в згину неодноридних пластин
Ключовi слова: теорiя пружностi, iзотропнi тша, неоднор^дн матерiали, згин
плит, напруження та деформащг
□-□
Рассмотрена трехмерная задача о напряженно-деформированном состоянии плиты, в которой характеристики упругости материала изменяются по толщине и описываются произвольными интегрируемыми функциями. Показано, что точные, в смысле Сен-Венана, аналитические решения краевой задачи можно получить, если распределение нагрузок удовлетворяет двумерному полигармоническому уравнению. Сформулированы приближенные теории для инженерных расчетов изгиба неоднородных пластин
Ключевые слова: теория упругости, изотропные тела, неоднородные материалы, изгиб плит, напряжения и деформации -□ □-
UDC 539.3
|doi: 10.15587/1729-4061.2016J5052|
ANALYTICAL STUDY OF THE BENDING OF ISOTROPIC PLATES, INHOMOGENEOUS IN THICKNESS
V. Plevako
Doctor of Technical Sciences, Professor** E-mail: V_Plev@ukr.net V. Potapov Doctor of Technical Sciences, Professor, Head of Department* E-mail: potapov@bigmir.net V. Kucenko PhD, Associate Professor* E-mail: vakutsenko@mail.ua I. Lebedinec PhD, Associate Professor* E-mail: igor13lebedinec@gmail.com I. Pe d o ry c h Senior Lecturer* E-mail: pedorich_i@mail.ru *The department of refrigerating and trading equipment and applied mechanics** **Kharkiv State Universiti of Food Technology and Trade Klochkovskaya str., 333, Kharkiv, Ukraine, 61051
1. Introduction
The heterogeneity is a characteristic feature of virtually all materials used in engineering and construction. It is due to a number of factors, which conditionally can be divided into three groups:
1. Action of environment (temperature fields, radioactive radiation, uneven humidity, etc.).
2. Manufacturing techniques peculiarities (rolling, forging, hardening of casting and concrete etc.).
3. Project implementation plan (availability of the reinforcement, layers of other materials, etc.), which is an important source of reduction in weight, size and cost of the projected designs.
The improvement of calculations of machine elements, constructions and structures is associated with taking into account the impact of heterogeneity of actual materials on the stressed and deformed states of elastic bodies. Therefore, a relevant task is the search for the methods of determining the displacements and stresses in three-dimensional bodies of a relatively simple form - slabs and plates made of inhomogeneous materials.
2. Analysis of scientific literature and the problem setting
Nowadays most of the studies in this direction consider a fairly limited number of the types of heterogeneities, for which
analytical solutions to the problems of elasticity theory were obtained [1, 2]. This is due to the fact that the difficulties faced during consideration of specific problems are relatively more complex than the similar problems of the classical theory of elasticity, since variable coefficients appear in its fundamental equations. These characteristic features offer explanation to the fact that analytical solutions to a limited number of tasks for the bodies of the simplest geometric shapes have been obtained until now: rectangular [1, 5, 6] and circular plate [2], structure [3] and dimensions [4, 6] and elementary dependencies of the characteristics of materials elasticity on the coordinates of the points of the body [5].
The analysis of the above mentioned publications demonstrated that the exact solutions to variable thickness of inho-mogeneous rectangular plates were considered, the influence of heterogeneity degree, the ratios of sides [1] and the heterogeneous boundary conditions [4] were taken into account. As result of the research, it was found that the influence of elasticity module on the plate of (FG) and (EG) materials rested on the elastic foundation is produced according to the exponential law [5, 6], and the influence of material properties and the strength of the elastic foundation on mechanical properties of the plate was studied.
The bodies of the type of sphere, half-space and the layered system were most often studied. The bodies with exponential [7, 8] or degree [9, 10] Young's law of the module change and stable [4], and in some cases even the variable
©
Poisson coefficient [11, 12] were considered. It was proven that such heterogeneities usually significantly affect both stressed and deformed state of bodies. The works [13, 14] are devoted to more complex dependencies.
Therefore, one can state that only quite a limited number of types of heterogeneities have been studied until now, for which analytical solutions to the problems of the theory of elasticity were obtained. At the same time, there are no exact, in the sense of Saint Venant, solutions to the boundary problems for the plates, inhomogeneous in thickness, exposed to the actions of the mass forces and the surface loads in a general case. Therefore, various approximate methods of solving the problems of mechanics of a deformed solid body have been developed over the recent decades, which allow receiving numeric solutions relatively easy.
But the vast majority of researchers feel the lack of a sufficient number of analytical solutions to the problems, which might serve as reliable test examples for approximate methods. Therefore, this research is devoted to searching for the new exact, in the sense of Saint Venant, solutions to the boundary problems for the plates, inhomogeneous in thickness, as well as the ways of building up approximate solutions that have the set degrees of accuracy.
3. The aim and the tasks of the research
The purpose of the research is to build up analytical solutions to the boundary problems for the plates, inhomo-geneous in thickness, exposed to the action of mass forces and surface loads.
To achieve this aim, the following tasks were set:
- to find the criteria, to which the distribution of the forces inside and on the body surface should correspond, if the material's heterogeneity is described by arbitrary integrated functions of one Cartesian coordinate;
- to obtain the exact, in the sense of Saint Venant, solutions of boundary problems for the plates, inhomogeneous in thickness, that can be used in engineering calculations, as well as when testing the existing approximate theories;
- to develop methods for building up approximate solutions that have a specified degree of accuracy.
4. Analytical solution of the problem about the bending of a plate, inhomogeneous in thickness, and its discussion
Let us assume that the plate, inhomogeneous in thickness, is loaded with mass forces X, Y, Z and is limited from the sides by the cylindrical surface r (Fig. 1) (the load of the plate surface is considered as a special case).
Fig. 1. Scheme of the plate, inhomogeneous in thickness (external forces are not shown)
We will try to set the conditions that the mass forces should meet for the problem on stressed and deformed states of a studied elastic body to have the solution, exact in the sense of Saint Venant, i. e. the plate thickness h will be considered a relatively small magnitude.
We will formulate the ways of building up so-called "technical" (approximate) theories of the bending of plates.
1. To solve the problem, we will connect the body with the Cartesian coordinate system, the O origin of which lies inside the plate and the axes Ox, Oy are directed by parallel bordering planes z=-h1 and z=h2, that is, the plate thickness h=h1+h2 (position of the point O by height will be chosen later from the condition of the simplest recording of final formulas). With this choice of axes, the material's heterogeneity will be described by the functions of only one coordinate z, which will be considered arbitrary and integrable.
As it is shown in the paper [14], the solution for the problem formulated above comes down to the search for the solutions of two linear differential equations with variable coefficients, one of which is of the fourth, and the other is of the second order. The equations have the form:
( 1 92L
dz2
V a4 dz2
dz2
—L
d ( 9L
—I at —
dz V dz
a^ dL
a4 dz2
AN + G,
A 2L = R,
If1 dN 1 = -G AfAi
dzl G dz
dz I G
(1) (2)
where
R _ d2
A = -
9x2
1 f^L-a, a4V dz 2
d-2_ dy2 '
d / . \ a2 fdA *
dz(aA )-a: lizr-a2
where A4, A2, A3 are the special solutions for these problems:
A 2 A -7 A A df2 9f3
A2Aj = Z, AA2 =—+—, dx dy
dy dx '
In this case, = a^z) (i=1, 2, 3, 4) may be expressed through the material's elasticity parameters like this:
AA3 --(Af2=Xi, AÍ3=y),
(3)
_ 1 _ a — , a0—
1 G 2
1 -
2v2 1 -v
V
E' a4" 1 -v2'
(4)
where G=G(z) is the module of the plate material's displacement; E=E(z) is the Young module; v=v(z) is the Poisson ratio.
If the functions Si(i=1, 2, 3, 4), Tj(j=1,2) are introduced for consideration with the help of the following dependencies:
S,
a„
■±(aL +
dz I dz
+ A, + a4AL-A2
_ JdL .) as,
S3 = A2L, S4 —A^+Ail,
T=-Gif+Aii- T=AN
(5)
then for determining the displacements ux, uy, uz and the stresses ox, oy,..., oxy, we will receive the following formulas:
9S, 9T 9S, 9T
ux =—1 +--1, uy =—1--1, uz = S2,
dx dy dy dx
G„ = -
aiA + 2G-T dy
a4A+
dx
Si — + 2G
S^ 2G
92Tj 9x9y'
92Tt 9x9y'
o = S3, T A + dTLi T =dS4 +
dx dy dy dx
'V
-— ^^ 2 )
—L''
n—1
-[ai (Ln—1 + «in ^)] — — x
x[Ln'—1 + Sin (^— eX2)] —(1 — Sin)
Ml I.
1 1 1
—N' = -—Nm G m I G
(9)
(10)
T xy = 2G - G dxdy
a2Tt a2Tt
9x2 9y2
(6)
Thus, when solving any problem of the theory of elasticity for bodies of one-dimensional heterogeneity, it is necessary in general case to find the solutions for the differential equations (1) and (2) by the corresponding boundary conditions. In some cases, it is enough to examine only one of them.
2. To find the solutions for the equations (1) and (2), let us first proceed to the dimensionless coordinates:
t x y J, z
s=x n=y z=h
and designate:
D, =—, D2 =—, D2 = D2 + D2,
1 2 dn 1 2
1 -Ai 1 -Al 1 = al
h h h2 f2 f3 e = J' K1 = p K2 = "¡f' ^ = jl' ^3 = jl".
Here l is the characteristic plate dimension on the plan. Now, the dependencies (3) take the form:
DX = Z, D2A2 = ^+ ^3, d2a3 = ^
d%
d%
(D2^2 = X, = y).
L = hl4 j^L0 + (eD)2 Li + (eD)4 L2 + ...J, N = hl3 [ N0 + (eD)2 N + (eD)4 N2 +...],
The formal solution for the equations (1) and (2) will be searched for in the form of the series, similar to how it was made in the paper [15]
(7)
(8)
where Ln, Nm(m, n=0,1,...) are the functions to be defined.
Let us substitute L and N to the equations (1) and (2) and collect similar members by the same degrees eD. Then, equating them to zero, we will obtain two recurrent sequences of ordinary differential equations:
We will enter the symbol here:
S.. = f°' i * j, ,j U i=j.
In these formulas, the differentiation by the coordinate Z. is marked with a stroke.
The boundary conditions for the functions Ln, Nn follow from the boundary conditions of the problem.
We agree to designate the value of any function by Z = -k4, with the top index in the form of a degree, and the value Z = k2 is designated in the form of an asterisk. First consider a separate case of the plate's load. Let us assume that the lower and the upper surfaces of the plate Z = -k4 and Z = k2 are free from the action of forces, and only mass forces affect the body.
Then the boundary conditions on the planes Z = -k4 and Z = k2 will be written down as follows:
o° = so = o, lTOx = D4S4+D2T2°=o,
k; = D2S0 - D4T2° = 0, (11)
o; = S3 = 0, lTZx = DiS,i+d2t;=0,
lTZy = D2s; - DjT* = 0.
Thus, S4 and T2o, as well as SJ and T2*, are interrelated by the Cauchy-Riemann conditions and, therefore, are two-dimensional harmonic functions ra° and ra*, i. e.
S4 = D4ra°, T2° = D2ra°,
S*4 = Djra*, T2 = D2ra* (D2ra°*= 0).
However, it is easy to prove that the influence of the functions ra° and ra* on the stressed and deformed states of the plate can be taken into account by using the arbitrary elements in the selection of the functions X2 and X3. So we put ra° = ra* = 0. Then the boundary conditions will take the form:
so _ s* _ so - S* - To - T* - 0
They are automatically satisfied if the condition is set:
l n = l; = 0, L' n = -Sonx
L'n =-Son^*, Nm = Nm = 0 (n,m = 0,1,...).
We will consider the following geometric characteristics of cross section of the plate:
Z Z
F = -J a4dZ; S =J a4ZdZ;
z z
I = J SdZ; F2 = J Gdz
Lo = eQ2-Qi -(eQ 2-Qî)-Vi-eQ
where
z z z
Qi = JMZ, Q2 =J dzj^2dz,
Vi
= I*- V2 = F* (S + ZF1 ) K2 Vi.
Solutions of other equations are determined by the following recurrent dependencies:
Ln = Mn - m; - M'n> 2,
In this case,
z z
Mn = - J dZ J [a2Ln-4 + a4 (Mln - M2n )] dZ,
-Kl -Kl
z z
Mln = J dzj] ^ [Lf.! + sln 2 )]+
-(1 -Sin )
—K1 -K Ia 4 2
Ln—2 \dZ,
M2n = J a, (L+ 8ln)dZ.
By solving equations of the sequence (10), we obtain: No = -Qs + Q 3x, Nm = -Rm + Rm x,
l4
S
A — -, O
[(eD)2 X +]T(eD)2n+2L n ],
(17)
Ti = l31(£D)2mU2
m=0
Here,
Uio =-qÎy* + e
T2 = - ¿(eD)2m+2Nm
h m=0
Uln =Mln -M2„ -
(12)
and choose the position of the origin of the coordinate O from the condition S*=0.
Then the solution of the first equation of the sequence (9) can be recorded as:
(Q2 -KQ 2" ) f - Q
n-Ml+(Mn-K2M'n- )4,
n Fi< v n 2 n <T
Wo=[q; -e(q; -k2q 2- )] 1 Wn =-M1n -(Mn -K2Mn* )4 ,
(13)
U2m =
Q3
r m _ r*
G F*
(14)
(15)
Q.3 = J MZ, Rm = J GdZj Nm-1dZ, x = F. (16) Pk =î*ï(№
-k. -k. -ki G F2 i+j=0
Thus, the stress and deformations have the form of the series, the members of which grow by degrees of the operator D2. Therefore, if the functions X, Y, Z are polyharmonic, i. e. satisfy a two-dimensional equation:
D2k||x(4, n, Z), Y(4, n, Z), Z(£,, n, Z)|| = 0(k = 1,2,3...)
(the variable Z plays the role of a parameter), the functions Ln and Nn will be also polyharmonic. As a result, the series (7), (8) and (17) will break and by using the formulas (17) and (6) we will find the exact solution to the problem. Let us remind that in this case the boundary conditions on the planes Z = -k4 and Z = k2 will be satisfied exactly.
3. Let us focus in more detail on the possible forms of setting the functions X, Y and Z in the search for the exact solutions to the problem on the stressed and deformed states of the plate.
It is easy to prove that the following features are characteristic for the polyharmonic functions:
if D2l> = 0, then D2(m+1)|n, P2||^ = 0 (p = ^2 + n2).
So, if there is any polyharmonic function ^ with arbitrary m, then the new solution O of the polyharmonic equation with a larger m can be built up by using the formula:
O = ^Pk,
In this case k
If we substitute series (7), (8) to the formulas (5) for S1, Ti, T2, S2, S3, S4, we will find:
Si = T" Uln, S2 = ^lO^fX,
h n=0 h n=0
14 »
S3 = ¿T K^fX,
Pk = 1« i(Z)p2i.
i=0
Here a¡j(Z) and a¡(Z) are the arbitrary integrable functions of the coordinate z.
We will receive additional solutions if we take: a) linear combinations of the known solutions;
2
or
b) the derivative of an arbitrary order by the parameter S
n Z, s).
ds .
c) the integral by the parameter with the weight function dependent on it
Jn, Z, s)f(s)ds.
As an example, here are a few more types of polyhar-monic functions, by which it is convenient to approximate volumetric loads during the solution of the problems of the theory of elasticity in Cartesian coordinates:
O =
O =
sin n£, cos n£,
sin nn
cos nn
Si lh
z z
Z-(eD)2 J dZj a2ZdZ-J GSdZ
-Ki -Ki
Zl -
1 -(eD)2
J |rdZ+ J dZj a2dZ
f;'
1 -(eD)2 J a2ZdZ
-Kj
S
^-(eD)2 J a2dZ^^f,
S3 = 0, -f = -(eD)
sy, - 41 f; 2
T= lh
1 -(ed)2 j |dz
^ T F
, = 4-(eD)2
l F,*1 ' 3
We will note that the problem of searching for the "homogeneous" solutions for a non-homogeneous plate, but from the electro elastic material, was studied earlier in [16].
5. The solution to the problem on the balance of the plate subject to the actions of the surface polyharmonic loads may be obtained from the solution given above.
Without limiting the studies generality, we will demonstrate this with the example of loading the upper surface Z = k2 with the normal efforts oZZ =o(%, n).
Put
Z = S(Z-k2) c(%, n).
(19)
4. The solution in the form of (17) is built in such a way that the boundary conditions on the planes Z = -k4 and Z = k2 are satisfied exactly. The conditions on the lateral area D have to be satisfied by using arbitrariness in choosing the functions A4, X2 and X3. The so-called "homogeneous solutions", which determine the stressed and deformed states of the plate with the load of the lateral side D, correspond to this arbitrariness.
Homogeneous solutions are easy to find if you put X=Y=Z=0. Then the above mentioned formulas imply that Ai, X2 and X3 are biharmonic functions.
Now, if we find Ln, and Nm from the dependencies (15) and (16) and substitute them to the formulas (17), we will find functions Sj, and Tj:
In this case, S(Z-k2) is the asymmetric impulse function [17], and we will first obtain the solution to the problem in the case when the plate is exposed to the mass forces Z in the form (19), that is, we might consider that the load is applied inside the plate under the surface Z = k2.
We have the equation for defining the functions X4, X2 and A3:
D% = s(z-k2)o(£,,n), D2X2 = D^2 + D2^3,
D2X3 = D2^2 - D^3, D2^2 = 0, D2^3 = 0. (20)
Obviously, it is sufficient to find any special solutions to these equations, whereas the general homogeneous solutions, corresponding to them, are included in the formulas (18). So we take:
X, = s(z-k2)f(^,t|), X2 = X3 = 0 (D4F = 0).
From the formulas (13), (14) and (16) we have:
Lo =|>!-e(C-K2)]F(%,n), N0 = 0,
where
e(Z-K2) =
z<k 2, Z = K2.
(18)
Here symbols ^123 designate the two-dimensional biharmonic functions. Between the two of them - and there is a link: D2x¥2 and D2^3 are adjoint harmonic functions.
Now using the dependencies (6), it is easy to find the formulas for determining the displacements and stresses in the plate. They, along with (18), allow receiving the solution to the problem on the balance of the plate, inhomogeneous in thickness, loaded with mass forces of polyharmonic type. Note that the boundary conditions may be satisfied exactly on the planes Z = -k1 and Z = K2, and on the lateral side r, they are "softened" in the sense of Saint Venant.
Hence, using the recurrent dependency (15), it is easy to find the required number of functions Ln (according to (16), Nm=0 for all m), and when substituting them to the formulas (5) and (6), one can find components of the displacement vector and the stress tensor.
In this case, the boundary conditions on the planes z = —k1 and Z = k2 of the plate will have the form (11). Now, to proceed to the problem of the influence of the surface forces, it is necessary to exclude the boundary area Z = K2 from the area belonging to the plate, and thereby, to "bring" the load o(£, ,n) onto the surface Z = K2.
In this case, for example, the function L0 will take the form:
Lo =v№, n).
If now we use the homogeneous solutions in the form of (18), then it is possible in principle to satisfy the "softened" boundary conditions on the surface r and by so doing, to solve exactly in the sense of Saint Venant the problem of balance of a non-homogeneous plate under the action of the surface load o on it.
Hence, if the external forces are represented by the poly-harmonic functions with the variables X, h, then the problem
S
of balance of a non-homogeneous plate may be solved exactly in the sense of Saint Venant.
The thinner the plate, the more exact the solution is. All the series introduced for consideration, break, and the result contains the finite number of members.
6. It follows from the formulas (17) that the order of each successive summand is proportional to the value e2 relative to the previous one. Therefore, if e is small, it is possible in the series (17) to be limited by taking into consideration only the first few members, while the others may be discarded; an excessive precision of formulas in the theory, which permits using the principle of Saint-Venant, is hardly appropriate.
If the external loads are not polyharmonic functions, they may be in a sense approximated by such functions and, thereby, a solution may be found. Moreover, if e<<1, then the precision, sufficient for practical purposes, can be obtained by leaving only one summand in the obtained formulas for calculating stresses and displacements.
However, a similar result may be found when the loads are not approximated by the polyharmonic functions. It is only necessary to break the series in the formulas.
The elementary «technical» theory of the bending of non-homogeneous plates may be obtained by leaving only one summand in the series. So, if X=Y=0, it is easy to find
ux = -e^DjW, uy = -eZD2w, uz = w,
o = —
o. =h
2eÇG /ri2 ^ _ 2eÇG ^2 , ^
, (D2 + vD2)w, oy = -l(1 — v) 1 27 y l(1 -v)
(vD2 + D2)w,
-L J MÇ-J ZdZ
2et
, Txy =-TÇGD1D2w,
t =--SD,D2w, t =--SD2D2w.
zx 1 1 ' zy 1 2
Here w is the solution to a biharmonic equation
14 K2
D4w = q, q = hi* J ZdZ.
-Ki
To simplify these dependencies, it makes sense to proceed from the dimensionless coordinates to the ordinary ones.
In a separate case when the material is homogeneous, the above given formulas will in fact coincide with those used in the theory of bending thin plates. The classical theory does not take into account only the stress oz.
If we leave only two summands in the series for ux, uy, uz, we will obtain a more precise method of calculation of the bending of non-homogeneous plates.
ux = -eDi
uy = -eD2
Zw - (eD/U^
Zw - (ed)2^^
uz = w+(eD)2W1hj.
We will note that it follows from the obtained results that the results, exact in the Saint Venant approximation, may be obtained in the case when the functions describing the distribution of loads inside the plate and on its surface satisfy two-dimensional polyharmonic equations (harmonic, biharmonic, etc.). This somehow limits the class of functions, which the real loads in a plate can correspond to, but allows receiving analytical solutions, exact in the Saint Venant approximation, which makes it possible, for example, to test numerical methods of solving the problems of mechanics of a deformed solid body.
5. Conclusions
A three-dimensional problem of bending the plate, in which the parameters of elasticity of the material vary by thickness and are arbitrary integrable functions, was examined. And the plate itself is exposed to the action of mass forces while the action of the surface loads is studied as a separate case:
1. We obtained analytical solution to the boundary problem by the operator methods in a case when the boundary conditions are satisfied exactly on the flat surfaces of the plates, and on the lateral surface - in the Saint Venant approximation.
2. It was theoretically proved that the exact, in the sense of Saint Venant, analytical solutions may be obtained if the plate is exposed to the action of mass and surface forces, distributed on the plate and on its surface by the two-dimensional polyharmonic law. In this case, the thinner the plate, the more exact the solution will be, since the corresponding solutions represent the series that contain a finite number of members.
3. The obtained solutions allow using them as an approximate, "technical" theory for engineering calculations of the stressed and deformed state of non-homogeneous plates.
4. It was demonstrated that the obtained formulas for the calculation of the bending of thin plates in the case of homogeneous material, transfer to the classic formulas of the theory of bending thin plates.
4
References
1. Zenkour, A. M. Exact solutions for variable-thickness inhomogeneous elastic plates under various boundary conditions [Text] / A. M. Zenkour, D. S. Mashat // Meccanica. - 2008. - Vol. 44, Issue 4. - P. 433-447. doi: 10.1007/s11012-008-9181-z
2. Li, X. Pure bending of simply supported circular plate of transversely isotropic functionally graded material [Text] / X. Li, H. Ding, W. Chen // Journal of Zhejiang University SCIENCE A. - 2006. - Vol. 7, Issue 8. - P. 1324-1328. doi: 10.1631/jzus.2006.a1324
3. Zhang, Z. M. An explicit solution for the elastic quarter-space problem in matrix formulation [Text] / Z. M. Zhang, W. Wang, P. L. Wong // International Journal of Solids and Structures. - 2013. - Vol. 50, Issue 6. - P. 976-980. doi: 10.1016/ j.ijsolstr.2012.12.001
4. Gao, Y. A refined theory of elastic thick plates for extensional deformation [Text] / Y. Gao, B.-S. Zhao // Archive of Applied Mechanics. - 2009. - Vol. 79, Issue 1. - P. 5-18. doi: 10.1007/s00419-008-0210-0
5. Salehipour, H. Closed-form elasticity solution for three-dimensional deformation of functionally graded micro/nano plates on elastic foundation [Text] / H. Salehipour, H. Nahvi, A. R. Shahidi // Latin American Journal of Solids and Structures. - 2015. - Vol. 12, Issue 4. - P. 747-762. doi: 10.1590/1679-78251398
6. Zenkour, A. M. Thermo-Mechanical Bending Response of Exponentially Graded Thick Plates Resting on Elastic Foundations [Text] / A. M. Zenkour, M. N. M. Allam, A. F. Radwan, H. F. El-Mekawy // International Journal of Applied Mechanics. - 2015. -Vol. 07, Issue 04. - P. 1550062. doi: 10.1142/s1758825115500623
7. Chondhury, P. Stresses in an elastic layer with varying modulus of elasticity [Text] / P. Chondhury // Bulletin ot the Calcutta Mathematical Society. - 1957. - Vol. 49, Issue 2. - P. 99-104.
8. Kogan, B. N. Stresses and strains in the coating with a continuously varying modulus [Text] / B. N. Kogan // Tr. Khark. avtomob.-dor. yn-ta. - 1957. - Issue 19. - P. 53-66.
9. Gibson, R. E. Some results concerning displacements in a non-homogeneous elastic layer [Text] / R. E. Gibson, P. T. Brown, K. R. F. Andrews // Zeitschrift für angewandte Mathematik und Physik ZAMP. - 1971. - Vol. 22, Issue 5. - P. 855-864. doi: 10.1007/bf01591813
10. Butter, H. Elastische Schacht und elastischer Halbraum bci mit der Tiete stetig abnehmenden Elastizitatmodul [Text] / H. Butter // Ingenieur-Archiv. - 1963. - Vol. 32, Issue 6. - P. 430-436. doi: 10.1007/bf00531790
11. Awojobi, A. O. Plane strain and axially symmetric problems of a linearly nonhomogeneous elastic half-space [Text] / A. O. Awojobi, R. E. Gibson // The Quarterly Journal of Mechanics and Applied Mathematics. - 1973. - Vol. 26, Issue 3. - P. 285-302. doi: 10.1093/qjmam/26.3.285
12. Awojobi, A. O. On the hyperbolic variation of elastic moduls in a non-homogeneous stratum [Text] / A. O. Awojobi // International Journal of Solids and Structures. - 1976. - Vol. 12, Issue 11. - P. 739-748. doi: 10.1016/0020-7683(76)90039-1
13. Puro, A. E. On the solution of equations of the theory of an inhomogeneous medium elasticit [Text] / A. E. Puro // Journal of applied mechanics. - 1975. - Vol. 11, Issue 3. - P. 50-55.
14. Plevako, V. P. Common solutions in problems of the theory of elasticity of inhomogeneous media [Text] / V. P. Plevako. - Kharkov: Osnova, 1997. - 159 p.
15. Plevako, V. P. Plane problems bending inhomogeneous bodies with the polynomial loadings [Text] / V. P. Plevako // Journal of applied mechanics. - 1978. - Vol. 14, Issue 3. - P. 59-65.
16. Ustynov, Y. A. Homogeneous solutions and the problems of transition from three-dimensional to two-dimensional problems for the boards of elektrouprugih materials with variable thickness properties [Text]: conference / Y. A. Ustynov. - Tbylysy: Yzd-vo «Metsnyereba», 1975.
17. Korn, H. Mathematical Handbook for Scientists and Engineers [Text] / H. Korn, T. Korn. - Moscow: Nauka, 1973. - 832 p.