Розглянуто задачi статики трансверсально-iзотропних пластин, як знаходяться за ди зосе-редженог сили. Проведено аналiз фундаменталь-них розв'яз^в, отриманих iз використанням рiзних наближених теорш. За допомогою розглянутих теорш тривимiрна задача теори пружностi зведена до двовимiрног. Розглянуто безмоментний пружний стан i стан вигину для трансверсальнонзотропних пластин. Проаналiзовано фундаментальт розв'яз-ки рiвнянь статики трансверсально^зотропних пластин, що побудован на базi узагальненог теори {т,п}-аппрокымацп для рiзних наближень
Ключовi слова: {т,п}-аппроктмащя, силова дiя, рiвняння статики, трансверсальнонзотропна пластина, полноми Лежандра
□-□
Рассмотрены задачи статики трансверсально-изотропных пластин, которые находятся под действием сосредоточенной силы. Проведен анализ фундаментальных решений, полученных с использованием различных приближенных теорий. С помощью рассмотренных теорий трехмерная задача теории упругости сведена к двумерной. Рассмотрены безмоментное напряженное состояние и состояние изгиба для трансверсально-изотропных пластин. Проанализированы фундаментальные решения уравнений статики трансверсально-изотропных пластин, построенные на базе обобщённой теории {т,п}-аппроксимации для разных приближений
Ключевые слова: {т,п}-аппроксимация, силовое воздействие, уравнения статики, трансверсально-и-
зотропная пластина, полиномы Лежандра -□ □-
UDC 539.3
|DOI: 10.15587/1729-4061.2017.96508]
ANALYSIS OF FUNDAMENTAL SOLUTIONS TO THE EQUATIONS OF STATICS CONSTRUCTED FOR TRANSVERSAL-ISOTROPIC PLATES
I. Bo kov
Postgraduate Student* E-mail: bokov.dev@gmail.com N. Bondarenko
PhD*
E-mail: Bondarenko.Natalya.Sergeevna@gmail.com
E. Strelnikova
Doctor of Technical Sciences, Professor, Senior Researcher* E-mail: elena15@gmx.com *Department of General Research in Power Engineering A. N. Podgorny Institute of Mechanical Engineering Problems NAS of Ukraine Pozharskoho str., 2/10, Kharkiv, Ukraine, 61046
1. Introduction
The thin-walled elements of structures (plates and shells) are the basic bearing elements of many contemporary designs and articles (vessels, aircraft, pipelines, cisterns, etc.). The action of the concentrated or local force influence on a thin-walled element causes dangerous stress concentration, which can lead to the destruction of the structure. That is why examining the stressed state of plates and shells, exposed to considerable force influences, has been, and remains, a theoretically and practically relevant task.
Until recently, the strength calculations of the thin-walled elements of structures have mostly employed the classical theory, based on the hypothesis about undeformable normal elements. In recent years, new composite materials have been utilized in engineering intensively in order to create protective coatings on the friction surfaces, as well as to fabricate different parts of equipment. Such materials typically possess low shear rigidity, which is why the hypotheses of Kirchhoff theory do not hold for them and the classical theory yields significant error in the calculations.
The use of contemporary composite materials leads to the need of constructing the refined theories of plates and shells, which consider phenomena that are connected to transverse
dislocations and compression. However, solving the problems on the theory of elasticity based on the refined theories in a three-dimensional setting is linked with considerable mathematical difficulties. That is why an issue about constructing the refined theories is closely connected to the problem of bringing three-dimensional tasks to the two-dimensional ones.
Thus, examining the SSS of plates under the action of concentrated force influences based on the refined theories is an important and relevant scientific and technical task.
2. Literature review and problem statement
The analytical solution of the axisymmetrical mixed problem for the isotropic half-space with the surface, elas-tically fixed beyond the circular region of the application of distributed load, is given in article [1]. When solving the given problem, the procedure of transition from the distributed load to the concentrated force was substantiated. They obtained the compact form of precise analytical solution of the problem on the concentrated force, applied to a halfspace with elastic surface.
Paper [2] explores the problem on the bend of a rectangular plate with rotationally fixed edges under the action of
©
concentrated force. The focus is on determining the angular forces and deflections.
The problem on the concentrated force that acts in the inner part of an infinite plate is solved in article [3]. This publication examines an isotropic plate of arbitrary thickness, in which the moduli of elasticity are any assigned functions. Solution of the problem presented is based on the classical solution for the concentrated force in a thin elastic plate.
Paper [4] addresses the displacement of a shallow spherical thin shell under the action of concentrated load. Deflections under the point of the application of concentrated loads on the spherical thin shells are calculated. This work also analyzes and compares the methods, which make it possible to obtain the refining calculations of shell deflections.
The influence of surface stresses on the stressed-strained state of elastic shell in is investigated in article [5]. In this study, the surface stresses are represented in the form of a static load, localized in the ultrathin layers of shell near its surface. The three-dimensional equations of elasticity are analyzed by the asymptotic method with the use of several asymptotic parameters.
Paper [6] presents results of analysis of a flexible plate of different thickness in accordance with the Kirchhoff model, the model of first order shift and the three-dimensional elastic model. Numerical algorithms of the method of difference energy and the method of finite elements are used to solve the problem. Solution of the contact problem and comparative evaluation of the obtained results are given.
An analysis of the publications over recent time allows us to conclude that the refined theory of the {m,n}-approxi-mation has not been commonly applied to study the SSS of plates and shells, which are exposed to the action of concentrated force.
This is linked to the fact that the solutions of problems on the concentrated influence, employing the theory of {m,n}-approximation, is connected with the larger mathematical difficulties than the solution of analogous problems based on classical theory. The mathematical apparatus, which is based on the application of a specialized G-function [7], makes it possible to overcome the indicated difficulties.
Among the publications, which apply a generalized theory of the {m,n}-approximation and the procedure of mapping, built with the help of specialized G-function, those articles should be noted that deal with the construction of fundamental solutions of the equations of statics for the transversal-isotropic plates based on the {1,0}-approximation for the state of bending [8] and based on the {1,2}-approximation for the zero spin stressed state and the state of bending [9, 10]. Worth noting are also the publications, in which authors obtained fundamental solutions of differential equations on thermoelasticity of the {1,0}- and {1,2}-approximation for the transversal-isotropic plates [11, 12].
3. The aim and tasks of the study
The aim of present study is the analysis of fundamental solutions of the equations of statics of the transversal-iso-tropic plates, built on the basis of a generalized theory of
the {m,n}-approximation for different approximations. This will make it possible to conduct strength calculations of the thin-walled elements of structures that contain dangerous stress concentrators in the form of concentrated and local force impacts.
To achieve the aim, the following tasks were set:
- to analyze the methods for reducing the three-dimensional problems of the theory of elasticity for plates to the two-dimensional ones;
- to compare fundamental solutions, obtained on the basis of the {1,0}- and {1,2}-approximation.
4. Materials and the methods for examining the decrease in order of the three-dimensional problems on the theory of elasticity
4. 1. System of equations of the theory of elasticity for transversal-isotropic bodies
The system of equations of the theory of elasticity for the transversal-isotropic bodies includes:
1) the Cauchy ratios [13]
e e —dv e — duz;
x 9x ' y 9y z 9z
9u 9u„ du du 9u„ 9u
ex ——-+—y; exz — —- + —-; ez ——y+—-; (1)
xy dy 9x xz dz 9x yz dz 9y
2) Hooke's law for the transversal-isotropic bodies [14]
Ox = E (ex + vey ) + Xoz ; oy = E (ey + vex ) + Xoz ;
o =
E'[ez+%+ey )] ;
t = Ge ; t = G'e ; t = G'e ,
xy xy ' xz xz ' yz yz '
where
EE= (1-v)E'
1 -v2' 1 -v-2 (v')2E/E''
v' E
—; g—r, 1 - v E 2(1 + v)
where E, E' are the Young moduli for the directions in the planes of isotropy and perpendicular to it; v, v', G, G' are the Poisson coefficients and the shear moduli, which correspond to these directions;
3) the equations of equilibrium [14]
dox 9T 9t
—-+——+—— + Fx = 0;
dx 9y dz x
9t do 9t —— + ^ + —^ + Fy = 0; dx 9y dz y
9t doz „ _
——+—— + —- + F = 0,
9x 9y dz where F = (Fx ,Fy ,Fz ) is the volume force vector.
(2)
4. 2. Description of methods for reducing the three-dimensional problems of the theory of elasticity for plates to the two-dimensional ones
Precise solution of the problems of the theory of elasticity for plates in a three-dimensional setting is associated with significant mathematical difficulties, in order to overcome which, the methods of approximation are employed. The method of decreasing the quantity of independent variables is one of them. In this case, a transition from the three-dimensional problems to the two-dimensional ones is accomplished in different ways.
The first group includes the approaches based on the acceptance of various hypotheses. The most popular are the hypotheses of classical Kirchhoff theory for plates. The same group includes an S. P. Timoshenko model, proposed in 1916 for the problems on the dynamics of rods.
Another technique for the reduction of three-dimensional equations to the two-dimensional ones consists in the expansion of the desired functions into exponential or functional series along the normal z coordinate and in the retention of specified sections of these series depending on the required accuracy and character of the problem. This approach is used, in particular, by the refined theory of the {m,n}- approximation [14]. Within the framework of the given theory, when deriving the two-dimensional equations, they employ the method of expansion of the desired and assigned functions of three-dimensional equations (1), (2) into the Fourier series by the Legendre polynomials Pk — Pk (z/h). Moreover, in the theory of {m,n}-approximation, m characterizes a quantity of retained members of the Fourier series for the components of SSS in the xOy plane; n is the number of retained terms of series for the desired and assigned functions in the direction of the Oz axis. In particular, the representations of components of the strain tensor within the {m,n}-approximation take the form [14]
where u, v, w0 are the analogs of displacements of the points of median surface of the plate; gx, gy are the analogs of the rotation angles of the normal; Nx, Ny, S are the analogs of membrane forces; Mx, My, H are the analogs of the bending moment and the torque; Qx, Qy are the analogs of transverse forces.
Components of the volume force vector are represented as [14]
F _ aP-+f« (x-_
(4)
The equations of statics in the case of {1,0}-approxi-mation for the transversal-isotropic plates, recorded in the dimensionless coordinate system (x1=x/h, x2=x/h, x3=x/h), contain [14]
- equations of Hooke's law:
N = bc
9u 9v
-+ v—
9x1 dx2 y
; N2 =B0
dv du
-+ v-
^Эх2 9x1
(5)
S = —B0 2 0
du dv Эх2 9x1
; h=i-VDo
' 2 0
dY1 + Эу2 Эх2 Эх1
Ml = Do
Qi = Ло where
dYi I vdY 2 dx1 Эх2 y
9w,
; M2 = Do
dY 2 +vdYi
V dx2
9x1
Yr
dx1
; Q2=Ло
Y2
dwo 9x.
2
B = JL-. D - _D_ - 2 1 ; л _
B _ • T ; D0 _ т-ч 2 _ ^ л 2 ; л о _
1 -v2
Eh2 31 -v2
3E/G''
ex _jCexkPk (x - y); exy exykpk;
k_0 k_0
n+2 n+1 exZ exzkpk (x—y); ez ezkpk-
k_0
where Pk = Pk (z/h) are the Legendre polynomials.
4. 3. Fundamental relationships and mathematical statement of the problem
We consider a transversal-isotropic plate of thickness 2h in the rectangular Cartesian coordinate system x, y, z. Concentrated force FC, applied at the origin of coordinates (singular point), acts on the plate.
In the case of {1,0}-approximation, the components of displacement vector and stress tensor are represented as:
ux _ uP0 + hYxp1;; uy _ vP0 + hYyP1;; uz _ w0P0;;
ox _ N7P0+3MfP1 (x - y); x 2h 0 2h2 1
S 3H °z _0; Txy _2hP0 + 2h2P1;
Qx
(3)
Txz _ (P0 -P2) (x - y),
E/G' is the shear compliance parameter; - equation of equilibrium
9N1 9S „ 9N2 9S —1 + —+ q1 _ 0; —± + — + q2 _ 0;(6)
9x1 dx2 1 ox2 dx1
ЭМ1 ЭН ЭМ2 ЭН
--Q1 + m1 _0; —± + ---Q2 + m2 _0;
ox1 ox2 ox2 ox1
dQi 9Q2 —^ + —^ + q3 — 0. 9x1 9x2
Let us note that the system of equations, which describes the zero spin stressed state in the case of {1,0}-approximation (first three equations (5), first two equations (6)), is similar to the system of equations, which describes the flat stressed state of isotropic plates within the framework of classical theory.
Within the framework of the {1,2}-approximation, the representation of components of displacement vector ux, uy, of stress tensor ax, ay, Txy and the components of volume force vector Fx, Fy takes the same form as in the case of {^^approximation (formulas (3), (4)), while the remaining components are determined as [14]:
uz — w0P0 + w1P1 + w2P2;
Txz=Qf (P. - P4 )+^t (P. - P3 )+5Qf (P- - P4 ) (x - y);
25 2 1 3 E
Q. =—Q0; D. =--; A. = — A.; E* = — ;
. 21 . . 31-v2 . 4 . E'
=I (P. - P2 )+3Rl (P. - P3 );
2h
3qz^ , 5qz:
Fz = —P. ^^ziP. + ^^P2.
2h . 2h
2h
equation of equilibrium
dN. dS „ 9N2 dS
—1+-+q. = O; —2+-+ q2 = 0;
dx. dx2 dx2 dx.
The system of equations for transversal-isotropic plates in the case of {1,2}-approximation, recorded in the dimen-sionless coordinate system, contains - equations of Hooke's law:
N. = B.
N2 = B.
du dv
-+v—
dx. dx2 y
dv du -+v-
dx2 dx.
+ X.R.;
+ X.R.;
S = —B. 2 0
dv du -+ —
dx. dx2 y
; h=^D0 2 0
dY 2 + dx. dx2 y
M. = D.
M2 = D.
Q.0 =A0
Q20 =A0
dY^ vdY 2
dx. dx
2
dY 2 + dx2 dx.
-XoR.;
Y + I w0
Y2 +--1 w. -
'2 dx21 0 .4
(7)
3M. dH dM2 dH
— Q.O + m. = --Q20+m2 = 0;
dx. dx2 dx2 dx.
9Q20 + = 0;9Q.^ dQ2.
dx. dx2
+ ^ - R. + q4 = 0;
dx. dx.
^ + dQ2 - 3R. + q5 = 0.
dx. dx2
The moments in relationships (5)-(8) are determined with accuracy to the magnitude Eh2, and the generalized efforts - to Eh.
A fundamental solution of system (5), (6) and system (7), (8) has a certain mechanical sense - this is the solution of the problem on the action of concentrated force on a plate. Dirac's delta function [.5] is located in the place of load functions in the right sides of resolving equations. Therefore, if the concentrated force FC, acts on the plate, then the volume force vector must be taken in the form
F (? ) = F Cs( ? ),
where r = (x,y,z) is the radius-vector of the force application point in the initial (dimensional) coordinate system; FC = (FxC,FyC,FzC ).
Then we obtain the expressions in the right sides of equations of equilibrium (6), (8) in the dimensionless coordinate system
Q.. = A. ^ Q2. = A0 ^
dx. dx2
Q.2 =-a-.4
Q22 = -Al ^22 .4
Y. (w.- 2w2 )
Y 2 (w.- 2w2 )
R0 =Q0 R. =Q.
where
w. + ^0
du dv -+-
dx. dx2 y
dY. , dY2^ dx. dx2 y
2 1 v' * 28 .
B„ =-' K=-E; A„ =--
Bo .-v2' 0 !-v 0 .5 E/G'
Q=7 (.-v)/E* .
mj = m*8(x.,x2 ); qj = q* S(x.,x2 ), (9)
where S(x.,x2 ) is the two-dimensional Dirac's delta function.
5. Results of examining the fundamental solutions, built with the use of refined theories
Let us examine the fundamental solutions, built on the basis of theory of the {m,n}-approximation. The fundamental solution of the equations of statics at bending in the case of {^^approximation is obtained in article [8]. The statement of a flat problem within the framework of this approximation coincides with the formulation of analogous problem for isotropic plates based on the Kirchhoff theory. The fundamental solution of the specified problem is built in paper [16]. The fundamental solutions for a flat problem and bending in the case of {^^approximation are received in articles [11, 12].
The expressions of inner force factors, obtained with the use of the generalized theory, take the form: - for the {1,0}-approximation
5.-v-2(v')2E* '
N. = -2n[q* fr (x.,x2 ) + (v + 2)^2 (x2,x. )} +
+ q2 jv¥. (x2,x. )-^2 (x.,x2 )}] ; a°)
N2 —-¿[q* (x1,x2 )-^2 (x2,x1)} +
+q2 {^1 tax ) + (v+2)^2 (x1,x2)}];
S — - ^^[q* {^1 (x2,x1 )-V^2 (x1,x2 )} + +q2 {^1 (x1,x2 )-v^2 (x2,x1)}];
M1 — - 2n[m* {^1 (x1,x2 ) + V^2 (x2,x1)+ 2^3 (x2,x1)} +
+m*2 {V^1 (x2,x1 )+^2 (x1,x2 )-2XP3 (x1,x2 )}-
-q {^ (x1,x2 ) + v^5 (x2,x1)}];
M2 — - 2n[m* K (x1,x2 ) + ^2 (x2,x1)-2^3 (x2,x1)} +
+m*2 {^1 (x2 ,x1) + V^2 (x1 ,x2 ) + 2XP3 (x1,x2 )} -
-q3 {v^5 (x1,x2) + ^5 (x2,x1)}];
H — ^¿[m* {(1-v)^2 (x1,x2 )-^3 (x1,x2 )+^4 (x2,x1)} +
+m*2 {(1-v)^2 (x2,x1 )-^3 (x2,x )+^4 (^ )}-
-(1-V) q3^6 (x1,x2)];
2* Q1 — [m*^7 (x2,x1) - m2^8 (x1,x2)] - J^;
Q2 — 2- [m2 ^7 (x1 ,x2) - m* ^ (x ,x2)] - J^f,
where
(x1,x2 ) —
x1 (x2 + 3xi;)
2r4
^3 (x1,x2) — 2r2G0,1 (ar) + G12 (ar);
^2 (x1,x2) —-
2r4
^ (x x)—- 1lnYr- 1x2 - x2.
(x1,x2)— 2ln2 4 r2 ;
^4 (x1,x2) — §fG0,1 (ar) + fiMlM G1,2 (ar);
(x1,x2 )— 2'r22 ;
2r4
GnV(z) is the specialized G-function [7]; - for the {1,2}-approximation
N —-— 1 2n
B B
qy °1(x1,x2 )+q2"AL O (x1,x2)+
+q**a2®3 (x1,x2)+q2a0° (x1,x2)+q3^0a2®5 (x1,x2)+
+2q1®6 (x2,x1)-2q2®6 (x1,x2)+
B v B v
+q**-A O (x2,x1)+q2^ ®1(x2,x1)+q**a0v®4 (x2,x0+ +q2a2 V0 (x2,x1)+q3V0V0 (x2,x0]-
' q**®7(x1,x2) + q2®7(x2,x1)-(x1,x2)
2nA
N2 —-— 2 2n
(11)
B B
q1^A 0 2 (x2,x1)+q2^° O (x2,x1)+
+q*a2®4 (x2,x1)+q2a2° 3 (x2,x1)+q3V2®5 (x2,x1 )-
-2q1®6 (x2,x1)+2q2®6 (x1,x2)+
* b0v , , \ . B0V / v
+q* it- °1 (x1,x2)+q^^ °2 (x1,x2)+
+q*a2 V0 (x1,x2 )+q2a2v®4 (x1,x2)+q3 V2v05 (x1,x2)]-
A 2q'
2nA
q**° 7 (x1,x2) + q2® 7 (x2,x1)-°9 (x1,x2)
S —-
1 -V 2n
B B
q**^0 2 (x1,x2)+q2B°L O (x2,x1)+
+q1a2® 4 (x1 ,x2) + q2a00 4 (x2 ,x1) + q3V2°8 (x1 ,x2)-
1 . 1
-q* v-[ °10 (x2,x1)+q2 v-1 0 (x2,x1)+
1 ,1
+q* v-[ °6 (x1,x2)-q2 v-1 (x1,x2)
M1 — - 2n '{miai2°11 (x1,x2 ) + m2a2°12 (x1,x2 ) +
-q5 o13 (x1,x2,a1)-q4 [ D 45 3 131 1 2 " 44 A1 42
013 (x1,x2,a1) +
+m** 014 (x1,x2,a1) + m2 015 (x1,x2,a1) -A1 A1
-q4ai2°16 (x1,x2 ) + 2mi*°15 (x2,x1,b1)-
2m2°15 (x1,x2,b1) +
+m1*a2!v®12 (x2,x1) + m2a2!v®11 (x2,x1) +
+°13 (x2,x1,a1)-q4v
D0
A1 42
013 (x2,x1,a1) +
* vD0 , \ * vD0 , v
+m*~T"L015 (x2,x1,a1) + m2 7 0°14 (x2,x1,a1) -A1 A1
-q4a22v°16 (x2,x1)}+
1+ A1a2
14A„Q„
017 (x1,x2 ) +
^ (x1,x2) — 1 |g00 (ar) + ^^Gu (ar
__Aa2 1
(x1,x2) — xf2G1j1 (ar); a2 — 3(1 + v)A0; r— Vr2^; -q5 ^^°17(x1,x2)-m**°18(x1,x2)-m2°18(x2,x1;
M2 = -2~{mia2012 (x2,x1) + m2a2®11 (x2,x1) +
2n
+q5 013 (x2,x1,a1)-q4
Do
A1 42
$13 (x2,x1,a1) +
+m1Dl $15 (x2,x1,a1) + m2D0 $14 (x2,x1,a1 ) -A1 A1
-q4a2®16 (x2,x1) - 2m1$15 (x2,x1,b1) + 2m2$15 (x1,x2,b1) +
+m1a2 VO11 (x1,x2) + m2a2v$12 (x1,x2) +
Do V2
+q5^3^ $13 (x1,x2,a1)-q4V
A1 42
* vD0 / \ * vD0 , v
+m1~T"°$14 (x1,x2,a1) + m2 7 0$15 (x1,x2,a1) -A1 A1
—q4ai v$16 (x1,x2)}+^ ]q4
1+
AA2
14A„Q„
$17 (x1,x2 ) +
+q5 $17 (x1,x2 )-mi*$18 (x1,x2) - m2$18 (x2,x1)};
1
Q22 = 28n {m'*b2$19 (x1,x2,b1) - m2b2$13 (x1,x2,b1) + * 14A°Q° / v
+m* ,0 0 $19 (x1,x2,a1) + A1
+m2 14^°Q° $13 (x2,x1,a1) - 14q5$18 (x2,x1) + A1
+q4
14A„Q„
$21 (x2,x1)i;
$13 (x1,x2,a1)+ R = ^2Q0
2nA
A1
q1$7 (x1,x2 )+q2$7 (x2,x1)-q3^H^ $9 (x1,x2)
^An
$17 (x1,x2 ) +
= Iq4 1 + ^L
1 6nA1 [ 4 14A°Q° , Aa2
+q5 FoL$17 (x1,x2 ) - mi$18 (x1,x2 ) - m2$18 (x2,x1,
where
H = - 12—n—{m^a22$12 (x1,x2 ) + m2al2$12 (x2,x1) +
+q5 ^Oa1- $19 (x1, x2, a1) q4
A1 42
$19 (x1,x2,a1) +
+m1D° $15 (x1,x2,a1) + m2D° $15 (x2,x1,a1 ) -
A1 A1 1
-q4ai$20 (x1,x2) - m.* $14 (x2,x1,b1) --m2 v—1 $14 (x1,x2,b1) + +m* v-1 $15 (x1,x2,b1)+ m2 v— $15 (x2,x1,b1 )}; 1
Q1° = 2n{m1b2$13 (x2,x1,b1 )-m2b2$19 (x1,x2,b1 )--q4$18 (x1,x2)- q4a2$21 (x1,x2)};
Q2° = 2n{-mi*b2$19 (x1,x2,b1) + m2b2$13 (x1,x2,b1) --q4$18 (x2,x1)- q4a2$21 (x2,x1)};
Q11=
2nA
q1$5 (x1,x2)+q2$ 8 (x1,x2)-q3 777$ 7 (x1,x2
a2 = a2 = 196D°Q°; b2 = 2A° ;
0 AA°; 1 9A1A0 ; 1 D°(1 -v);
1
A =B° + A^iOO; A1 = D°+ -A2Qo;
^ i \ , V x1 (x22 - 3x2) , ,
$1 (xx1,:x2 ) = 22~G°1 (a°r ) + q„4-- G1,2 (a°r );
2r4
$2(x1,x2)=2r2G01 (a°r)+ 2(2r2—~Gn(a°r);
/ \ 3 / \ x1 (x2 3x2) f v $3 (x1,x2 ) = - 8x1G1,° (a°r)--^- G2,1 (a°r);
t \ 1 I \ x2 (Sx2 - x2)
$4 (x1,x2 ) = -8x2G10 (a°r)--^- G2,1 (a°r);
$5 (x1,x2 $6 (x1,x2
x2 x2
G0° (a°r)+ 1r2 2G11 (a°r)
) = - ^^^ $7 (x1,x2 ) = 2^01 (a°r); 2r4 r
Q21 =
2nA
A
q.$8 (x1,x2)+q2$5 (x2,x1)-q3 t-^7$7 (x2,x1)
1
Q12 = 28n {-m1b2$13 (x2,x1,b1) + m2b2$19 (x1,x2,b1) +
* 14A °Q° / \ * 14A °Q° / v
+m1 ° $13 (x1,x2,a1) + m2 ° $19 (x1,x2,a1) -
A1 A1
14q5$18 (x1,x2)+q4
14A„Q„
CLA T
1 A1
$21 (x1,x2)i;
$8 (x1,x2 ) = xrx2G11 (a°r); $9 (x1,x2) = G°,° (a°r);
x1 (x2 + 3x2)
$10 (x1,x2 )= 2r
$11 (x1,x2 ) = -|x1G1,0 (a1r)--1 ( 8- 2 ) G2,1 (a1r);
/ \ 1 / \ x2 (3x2 - x2) / \ $12 (x1,x2 )=- 8x2G1,° (a1r )--- 2,1 (a1r);
1
ф13 (x1,x2,c) = ö
G0,0 (cr)+ X1 2X2 G1,1 (cr)
л i \ 3x^ , ч x1 (x2 - 3x2) , ч Ф14 (X1,X2,c) = 22 G0,1 (cr)+ 2r4-"G1,2 (cr);
Ф15 (X1, X2, c) = 7XjG0,1 (cr)+ X2 ^ X2 ) G1,2 (cr);
ф1б (x1, x2 )=-21yln 2 " -"212G0,0 (a1r) 1 8 2 G2,0 (a1r);
2x
ф17 (X1,X2)= G0,0 (a1r); ф18 (X1,X2(a1r); ф19 (x1,x2,c)=x^G1,1 (cr);
Ф20 (X1,X2 ) =--4-^2,0 (a1r);
021 (x1,x2) — -"2LG1,0 (a1r).
Graphs of internal force factors (10), (11) (Fig. 1-8) are constructed along the X-axis in the dimensionless Cartesian coordinate system x1, x2, determined with accuracy to the half-thickness of plate h. Fig. 1-3 show the generalized membrane forces for a flat problem. Fig. 4-6 show generalized bending and twisting moments; Fig. 7-8 - generalized transverse forces for the state of bending.
Fig. 1. Generalized membrane effort N,: curve of green color — generalized effort, obtained within the framework of {1,0}-approximation; curve of red color — within the framework of {1,2}-approximation
When conducting the numerical studies, components of the expansion of volume force vector were determined from formulas (9), in which
m**— m2— q1— q2— q3— q4—q5— 1.
The computations are performed for the following values of elastic constants of the transversal-isotropic material:
E* = 5; v = 0,3; v' = 0,07; E/G' = 40.
Fig. 2. Generalized membrane effort N2: curve of green color — generalized effort, obtained within the framework of {1,0}-approximation; curve of red color — within the framework of {1,2}-approximation
Fig. 3. Generalized membrane effort S: curve of green color — generalized effort, obtained within the framework of {1,0}-approximation; curve of red color — within the framework of {1,2}-approximation
Fig. 4. Generalized bending moment M,: curve of green color — generalized moment, obtained within the framework of {1,0}-approximation; curve of red color — within the framework of {1,2}-approximation
Fig. 5. Generalized bending moment M2: curve of green color — generalized moment, obtained within the framework of {1,0}-approximation; curve of red color — within the framework of {1,2}-approximation
Fig. 6. Generalized torque H: curve of green color — generalized moment, obtained within the framework of {1,0}-approximation; curve of red color — within the framework of {1,2}-approximation
allowed us to estimate the refinement, introduced by the use of equations of the {1,2}-approximation for transver-sal-isotropic plates, exposed to the action of concentrated force, instead of applying the analogous equations of {1,0}-approximation.
Fig. 7. Generalized transverse forces Q; curve of green color — generalized transverse force, obtained within the
framework of {1,0}-approximation (i=1); curve of red color — within the framework of {1,2}-approximation (i=10)
We carried out numerical calculations, which make it possible to compare the fundamental solutions (10), (11), built on the basis of generalized theory in the variants of {1,0}- and {1,2}-approximation, respectively. These studies
Fig. 8. Generalized transverse forces Qj: curve of green color — generalized transverse force, obtained within the framework of {1,0}-approximation (j=2); curve of red color — within the framework of {1,2}-approximation (j=20)
6. Discussion of results, obtained with the help of approximated theories
We built fundamental solutions of equations of the theory of elasticity on the basis of refined theory of the {m,n}-approximation. The selected theory is the most acceptable for reducing the three-dimensional equations to the two-dimensional ones since it is not based on any hypotheses, but employs the method of I. N. Vekua for the expansion of desired functions into the Fourier series by the Legendre polynomials [14]. This approach makes it possible to examine not only the thin plates, but also the plates of medium and large thickness. In this case, the accuracy of the obtained solutions depends on the number of retained summands in the expansions of the assigned and desired functions. In addition, this theory makes it possible to consider transverse shearing and normal stresses.
The numerical studies conducted have demonstrated that the inner force factors of the zero spin stressed state and the state of bending, obtained with the use of equations of both the {1,0}- and {1,2}-approximation, have identical character, and their numerical values differ insignificantly. This research confirms that the generalized theory in the variant of {1,0}-ap-proximation makes it possible to ensure sufficiently high accuracy of the approximation of a three-dimensional problem of the theory of elasticity to the two-dimensional one.
7. Conclusions
1. We examined the methods of reducing the three-dimensional problems of the theory of elasticity for plates to the two-dimensional ones. An analysis of the approaches to decreasing the quantity of independent variables allowed us to conclude that when calculating the thin-wall elements of structures, made of contemporary composite materials, it is necessary to use the refined theories of plates and shells
for the concentrated force influences. These theories make it possible to evaluate the phenomena, connected to taking account of transverse shifts and compression. Since the classical theory of Kirchhoff-Love does not consider these phenomena, then this theory yields significant error during calculations. The refined theories also make it possible to examine not only the thin plates, but also the plates of medium and large thickness.
2. We analyzed the fundamental solutions of the equations of statics, obtained on the basis of generalized theory in the variants of {1,0}- and {1,2}-approximation for the purpose of determining the refinement, which is introduced by the retention of a large number of terms in the expansion series of the desired functions. The numerical studies are con-
ducted, which confirmed that the expansion of the desired functions into series by the Legendre polynomials from the thickness coordinate and the retention of particular sections of these series makes it possible to obtain a solution of the problem with required accuracy. The order of the resolving system of equations depends on the selection of approximation and at large m and n it can be sufficiently high. The use in practice of such equations is connected with the considerable mathematical complexities. That is why, depending on the specific character of the examined problems, one selects such approximation so that the resolving equations would take the simpler form and, at the same time, would reflect the specific character of mechanical behavior of plate with a sufficient degree of accuracy.
References
1. Khapilova, N. S. The exact solution of the problem on a concentrated-force action on the isotropic half-space with the boundary fixed elastically [Text] / N. S. Khapilova, S. V. Zaletov // St. Petersburg Polytechnical University Journal: Physics and Mathematics. -2015. - Vol. 1, Issue 3. - P. 287-292. doi: 10.1016/j.spjpm.2015.11.004
2. Shi, W. Bending of a rectangular plate with rotationally restrained edges under a concentrated force [Text] / W. Shi, X.-F. Li,
C. Y. Wang // Applied Mathematics and Computation. - 2016. - Vol. 286. - P. 265-278. doi: 10.1016/j.amc.2016.04.029
3. Spencer, A. J. M. Concentrated force solutions for an inhomogeneous thick elastic plate [Text] / A. J. M. Spencer // Zeitschrift fur angewandte Mathematik und Physik. - 2000. - Vol. 51, Issue 4. - P. 573. doi: 10.1007/s000330050018
4. Wang, F. Deflection Solutions for Concentrated Force on Spherical Shell [Text] / F. Wang, X. Wang // Proceedings of the 5th International Conference on Civil Engineering and Transportation 2015. - 2015. doi: 10.2991/iccet-15.2015.50
5. Rogacheva, N. N. The effect of surface stresses on the stress-strain state of shells [Text] / N. N. Rogacheva // Journal of Applied Mathematics and Mechanics. - 2016. - Vol. 80, Issue 2. - P. 173-181. doi: 10.1016/j.jappmathmech.2016.06.011
6. Trushin, S. Numerical Evaluation of Stress-Strain State of Bending Plates Based on Various Models [Text] / S. Trushin,
D. Goryachkin // Procedia Engineering. - 2016. - Vol. 153. - P. 781-784. doi: 10.1016/j.proeng.2016.08.242
7. Khizhnyak, V. Mixed problem in the theory of plates and shells [Text]: tutorial / V. Khizhnyak, V. Shevchenko. - Donetsk: DonGU, 1980. - 128 p.
8. Bokov, I. P. Fundamental solution of static equations of transversely isotropic plates [Text] / I. P. Bokov, E. A. Strelnikova // International Journal of Innovative Research in Engineering & Management. - 2015. - Vol. 2, Issue 6. - P. 56-62.
9. Bokov, I. Construction of fundamental solution of {1,2}-approximation static equations of momentless stress state for transversely-isotropic plates [Text] / I. Bokov, N. Bondarenko, E. Strelnikova // ScienceRise. - 2016. - Vol. 8, Issue 2 (25). - P. 41-48. doi: 10.15587/2313-8416.2016.76534
10. Bokov, I. Investigation of stress-strain state of transversely isotropic plates under bending using equation of statics {1,2}-approximation [Text] / I. Bokov, N. Bondarenko, E. Strelnikova // EUREKA: Physics and Engineering. - 2016. - Vol. 5. - P. 58-66. doi: 10.21303/2461-4262.2016.00159
11. Bondarenko, N. The fundamental solution of differential equations of thermoelasticity {1,0}-approximation for transversely isotropic plates [Text] / N. Bondarenko // Proceedings of the Institute of Applied Mathematics and Mechanics of National Academy of Sciences of Ukraine. - 2009. - Vol. 18. - P. 11-18.
12. Bondarenko, N. Fundamental solution {1,2}-approximation the membrane thermoelastic state transversely isotropic plates [Text] / N. Bondarenko, A. Goltsev, V. Shevchenko // Reports of National Academy of Sciences of Ukraine. - 2009. - Issue 11. - P. 46-52.
13. Gorshkov, A. The theory of elasticity and plasticity [Text] / A. Gorshkov, E. Starovoytov, D. Tarlakovskiy. - Moscow: FIZMATLIT, 2002. - 416 p.
14. Pelekh, B. Laminated anisotropic plates and shells with stress concentrators [Text] / B. Pelekh, V. Lazko. - Kyiv: Science thought, 1982. - 296 p.
15. Vladimirov, V. Generalized functions in mathematical physics [Text] / V. Vladimirov. - Moscow: Science, 1976. - 280 p.
16. Shevchenko, V. Methods of fundamental solutions in the theory of thin shells [Text]: dis. ... Dr. Sc. Phys.-Math. / V. Shevchenko. -Kazan, 1982. - 332 p.