CC BY
9doi: 10.18720/MCE.82.11
Nonlinear parametric oscillations of viscoelastic plate of
variable thickness
Нелинейные параметрические колебания вязкоупругой пластинки переменной толщины
M.M. Mirsaidov,
Tashkent Institute of Irrigation and Agricultural
Mechanization Engineers, Tashkent, Uzbekistan,
R.A. Abdikarimov*,
Tashkent Financial Institute, Tashkent,
Uzbekistan,
N.I. Vatin,
Peter the Great St. Petersburg Polytechnic University, St. Petersburg, Russia V.M. Zhgutov,
"Kitezh" LLC Architecture and Construction Company, St. Petersburg, Russia, D.A. Khodzhaev, B.A. Normuminov,
Tashkent Institute of Irrigation and Agricultural Mechanization Engineers, Tashkent, Uzbekistan,
Д-р техн. наук, заведующий кафедрой М.М. Мирсаидов,
Ташкентский институт инженеров
ирригации и механизации сельского
хозяйства, г.Ташкент, Узбекистан,
д-р физ.-мат. наук, профессор
Р.А. Абдикаримов*,
Ташкентский финансовый институт,
г. Ташкент, Узбекистан,
д-р техн. наук, профессор Н.И. Ватин,
Санкт-Петербургский политехнический
университет Петра Великого,
г. Санкт-Петербург, Россия,
канд. техн. наук, генеральный директор
В.М. Жгутов,
ООО «Архитектурно-строительная компания «Китеж», г. Санкт-Петербург, Россия,
канд. физ.-мат. наук, доцент Д.А. Ходжаев, старший преподаватель Б.А. Нормуминов,
Ташкентский институт инженеров ирригации и механизации сельского хозяйства, г.Ташкент, Узбекистан
Key words: thin-walled structures; plates; shells; Ключевые слова: тонкостенные конструкции; variable thickness; periodic load; dynamic пластины; оболочки; переменная толщина; instability domain периодическая нагрузка; область динамической
неустойчивости
Abstract. Isotropic viscoelastic plates of variable thickness under the effect of a uniformly distributed vibration load applied along one of the parallel sides, resulting in parametric resonance (with certain combinations of eigenfrequencies of vibration and excitation forces) are considered in the paper. It is believed that under the effect of this load, the plates undergo the displacements (in particular, deflections) commensurate with their thickness. Geometrically nonlinear mathematical model of the problem of parametric oscillations of a viscoelastic isotropic plate of variable thickness is developed using the classical Kirchhoff-Love hypothesis. Corresponding nonlinear equations of vibration motion of plates under consideration are derived (in displacements). The technique of the nonlinear problem solution by applying the Bubnov-Galerkin method at polynomial approximation of displacements (and deflection) and a numerical method that uses quadrature formula are proposed. The Koltunov-Rzhanitsyn kernel with three different rheological parameters is chosen as a weakly singular kernel. Parametric oscillations of viscoelastic plates of variable thickness under the effect of an external load are investigated. The effect on the domain of dynamic instability of geometric nonlinearity, viscoelastic properties of material, as well as other physical-mechanical and geometric parameters and factors (initial imperfections of the shape, aspect ratios, thickness, boundary conditions, excitation coefficient, rheological parameters) are taken into account. The results obtained are in good agreement with the results and data of other authors. The convergence of the Bubnov-Galerkin method is verified.
Аннотация. Рассматриваются изотропные вязкоупругие пластинки переменной толщины, находящиеся под действием равномерно распределённой вибрационной нагрузки, приложенной по одной из параллельных сторон, приводящей (при определённых сочетаниях частот собственных колебаний и возмущающей силы) к параметрическому резонансу. Считается, что под воздействием указанной нагрузки пластинки допускают перемещения (в частности, прогибы), соизмеримые с их толщиной. Разработана геометрически нелинейная математическая модель задачи о параметрических колебаниях вязкоупругой изотропной пластинки переменной толщины с использованием классической гипотезы Кирхгофа-Лява. Выведены соответствующие нелинейные уравнения колебательного движения рассматриваемых пластинок (в перемещениях). Предложена методика решения рассматриваемой нелинейной задачи на основе применения метода Бубнова-Галеркина при многочленной аппроксимации перемещений (и прогиба), а также численного метода, использующего квадратурные формулы. В качестве слабо-сингулярного ядра выбрано ядро Колтунова-Ржаницына с тремя различными реологическими параметрами. Исследованы параметрические колебания вязкоупругих пластин переменной толщины под воздействием внешней нагрузки. При этом осуществлялся учёт влияния на области динамической неустойчивости геометрической нелинейности, вязкоупругих свойств материала, а также других физико-механических и геометрических параметров и факторов (начальных несовершенств формы, соотношений сторон, толщины, граничных условий, коэффициента возбуждения, реологических параметров). Полученные результаты хорошо согласуются с результатами и данными других авторов. Сходимость метода Бубнова-Галеркина проверена.
1. Introduction
Construction of the dynamic instability domain of thin-walled elastic systems is a very urgent and important task; a multitude of publications are devoted to this issue, of which, first of all, the studies carried out in [1-6] should be mentioned. Beginning with the article by N.M. Belyaev [1], parametric oscillations become the subject of numerous studies in application to various mechanical systems with distributed parameters, in particular, to rods, plates and shells. In the fundamental monograph V.V. Bolotin [2] has analyzed in detail the occurrence of parametric resonance in elastic thin-walled systems. An extensive monograph by T. Schmidt [6] contains a detailed review of many aspects concerning these issues.
Solution of linear problems of parametric oscillations of thin-walled structures (such as plates and shells) often is reduced to the use of ordinary differential equations with periodic coefficients - to the Mathieu-Hill equations. At certain ratios of coefficients, the solution of this type of equation increases indefinitely, since the corresponding oscillations of thin-walled structures are resonant in nature. The range of parameters that cause these oscillations is called the domain of dynamic instability (DDI). An assessment of the DDI is one of the central problems of parametric oscillations of thin-walled structures.
Solution of the problems connected with the study of resonant phenomena in geometrically nonlinear cases of viscoelastic thin-walled structures of variable thickness has been developed inadequately. The development of a methodology for solving such an important class of problems is of considerable theoretical and practical interest.
Note the following important publications in recent years devoted to the solution of this scientific problem:
1) the A.E. Kalinnikov and N.P. Kuznetsov article [7] presents the algorithm and results of a numerical solution of the problem of transverse oscillations of a straight rod under longitudinal force, harmonically varying in time;
2) the article by Q. Yan, H. Ding, L. Chen [8] contains an analysis of nonlinear parametric oscillations of the Timoshenko viscoelastic beam;
3) the article by L.D. Akulenko, S.A. Kumakshev, S.V. Nesterov [9] is devoted to the study of eigenfrequencies and modes of parametric oscillations of mechanical systems (on the example of small parametric oscillations of a mathematical pendulum of variable length);
4) the article by V.V. Karpov, O.V. Ignat'ev, A.A. Semenov [10] presents an analysis of the stressstrain state of shallow shell structures of double curvature, reinforced from the concave side by a various number of stiffeners. Mindlin-Reissner shell deformation theory is used, which accounts for geometrical nonlinearity and transverse shears, as well as for discrete introduction of stiffeners with contact between the stiffener and the shell along thestrip;
5) in the article by L. Kurpa, O.S. Mazur, Ya.V. Tkachenko [11] a numerical-analytical method for studying parametric oscillations of plates under the effect of static and dynamic loads (bearing a periodic nature) is substantiated;
6) the articles [12] (M. Darabi, R. Ganesan), [13] (H.Q. Huynh, H. Nguyen, H.L.T. Nguyen), [14] (R. Kumar, S.C. Dutta, S.K. Panda) and [15] (R. Kumar, S. Mondal, Sh. Guchhait, R. Jamatia) contain an analysis of results of studies of dynamic stability of plates (of various types) subjected to harmonic loads (with and without geometrical nonlinearity);
7) in the article by I.D. levzerov [16] the stability problems for bars and plates are considered. Variation formulations are used for the stability problem;
8) the article by V.M. Dubrovin and Т.А. Butina [17]; the authors propose a method for calculating dynamic stability of a cylindrical shell under axial compressive dynamic load;
9) in the article by A.A. Mochalin [18] on the basis of the V.Z. Vlasov semi-momentless theory, the solution of the problem of dynamic stability of an isotropic cylindrical shell of variable thickness is considered in the direction of its generatrix under external pressure;
10) in the article [19] (M.K. Usarov) is dedicated to improvement of plate theory in order to take into account forces, moments and bimoments, generated by nonlinear law of displacement distribution in plate cross-sections;
11) the article [20] (T. Dey, L.S. Ramachandra) is devoted to parametric oscillations of shells (of various types) subjected to static and dynamic loads of periodic nature (equations of oscillatory motion are obtained by the authors using the Donnell's theory of shell);
12) in the article [21] the author V.M. Budanov has proposed a special linear transformation for expressing the general solution of a second-order differential equation with periodic coefficients; the author has shown with numerical experiments that there exist the periodic solutions of the auxiliary system outside the instability domains of solutions of the Mathieu equation, and the solutions obtained for the DDI are in good agreement with known results;
13) in the article [22] (B.Kh. Eshmatov) is described the analyses of the nonlinear vibrations and dynamic stability of viscoelastic orthotropic plates. The models are based on the Kirchhoff-Love hypothesis and Reissner-Mindlin generalized theory (with the incorporation of shear deformation and rotatory inertia) in geometrically nonlinear statements. With a combination of the Bubnov-Galerkin and the presented method, problems of nonlinear vibrations and dynamic stability in viscoelastic orthotropic rectangular plates have been solved;
14) in the articles [23] (R.A. Abdikarimov, D.P. Goloskokov) and [24] (R.A. Abdikarimov, B.A. Khudayarov) a problem of nonlinear fluctuations of a viscoelastic plate of a variable thickness is considered. The resolving system of the integrodifferential equations of a problem is received. The numerical decision is constructed by Bubnov-Galerkins method with the subsequent application of quadrature formulas for approximation of integrals;
15) in the articles [25, 26] (B.Kh. Eshmatov), accordingly, the dynamic stability and nonlinear vibrations problems of viscoelastic orthotropic and isotropic plates is considered in a geometrically nonlinear formulation using the generalized Timoshenko theory. The problem is solved by the Bubnov-Galerkin procedure combined with a numerical method based on quadrature formulas. The effect of viscoelastic and inhomogeneous properties of the material on the dynamic stability and vibrations of a plate is discussed;
16) in the article [27] (B.Kh. Eshmatov, Kh. Eshmatov, D.A. Khodzhaev) the problem of flutter of viscoelastic rectangular plates and cylindrical panels with concentrated masses is studied in a geometrically nonlinear formulation. The behavior of viscoelastic rectangular plates and cylindrical panels is studied and the critical flow velocities are determined for real composite materials over wide ranges of physicomechanical and geometrical parameters.
Practical value of studies on the assessment of stability of thin-walled structures consists in consideration of the magnitude of existing loads and the law of their variation over time, as well as of all important features of both the structure itself and the properties of materials for their fabrication; this allows a correct estimation of load-bearing capacity of such structures.
Despite the great number of publications devoted to the study of parametric oscillations of thin-walled structures, the number of papers devoted to the study of dynamic stability of these structures with account of viscoelastic properties of their material is few.
2. Methods
Consider a viscoelastic rectangular plate (with sides a and b) of variable thickness h = h(x, y) , made of homogeneous isotropic material (Figure 1). Assume that the plate is subjected to an external dynamic load acting along the edge a, having a periodic character: P(t) = P0 + P cos0t (here
Pq , P = const; 0 - the excitation frequency). Also assume that the plate has initial deflections.
P(t) =p0+p1,:os€k
Я
k
P(t} =P0+PiSQs6>t
Figure 1. Viscoelastic rectangular plate of variable thickness under action of external periodic
load
In [28] R.A. Abdikarimov and V.M. Zhgutov have shown that in this case nonlinear equations of the oscillatory motion of the plate (with respect to deflection w = w(x,y, t) and longitudinal displacements
d 2 w
u = u(x,y, t) , v = v(x,y, t)) under the influence of force P(t)—— with account of initial deflection (at
dx
appropriate boundary and initial conditions) can be written in the form:
(1 -r* )h
dSx
dx
+ u
dsy 1 -ju dY
Л
dx
+
xy
2 дУ
\ 1 - u dh + ~(Sx +Usy j +
dh ( dx
2 dyYxy
+
+
1 -j
E
Px -ph
1 - jU d2u E ~dtT
= 0
(1 -r* )h
ids
, dSx 1 -UdYxy ].dh{
y + u—x + y
dy dy 2 dx J dy
+—(sy +Usx
1 - u dh
dxYxy
+
+
1 -U
E
Pv -Ph
1 -u2 d2v E dt2
= 0
(1-Г*
DV 4w + 2 dD—V2 w + 2——V2 w + V2 DV2 w -
dx dx
dy dy
-(1 -u)
d2D d w „ d2D d2w d2D d 2w
2 л"
dx2 dy2 dxdy dxdy dy2 dx2 J
- 2
- + -
22
V
■ +
+ h
(
dx
dw dw
— g +--T
^ x ^ xy
dx dy
^ di
+ -
J
dy V dy
dw dw
g„ +--T
dx
+ -
+ -
dh
dw
dw
d2 w
dh dx
d2 w
dw dw
— g +--T
x xy
dx dy
Л
+
J
dy V dy
w w w w
^Ту 1++q-ph^r=0
dx2
dt2
(1)
d
Introduce to equations (1) the following dimensionless parameters:
u u =т; _ v w v = — ; w = — ho ho — Wo Wo =~T ho ; x x a ; y~- •y ; hh = h ; Л ho a = b ' 5 = b ho
q q f b Y - = E h ; px = Px ; p = E Fy Pz. E ; © m cot ; 4'} ; ¿o m =P P ' кр = P_ P кр
Using previous notations, a system of integral-differential equations (IDE) in a dimensionless form is obtained.
To solve the system of nonlinear IDE (1), the sought for functions u(x,y, t), v(x,y, t), w(x,y, t), w0(x, y) are presented as an expansion in known (basic) functions <pnm (x, y), <nm (x, y), y/nm (x, y), W0nm (x, y) satisfying the given boundary conditions:
N M N M
u(x, ^ t) = Z Z unm (t<nm y ) ^ t) = Z Z vnm (t}<nm y)
n=1 m=1 n=1 m=1 (2)
N M N M
W(^ У, t) = Z Z Wnm (tWnm (X, y), W0 (X, y) = Z Z W0nmW 0nm (x, y)
n=1 m=1 n=1 m=1
Where Unm = Unm (t) , Vnm = Vnm (t), Wnm = Wnm (t) - are unknown time functionsi <m (x, y)>
<nm (x,y), ¥nm (x,y) , ^0nm (x,y), n = 1,2,...,N; m = 1,2,...,M - known (basic) coordinate functions
that satisfy the given boundary conditions of the problem (described, for example, in Volmir's monograph [4]).
As a rule, algorithms for solving dynamic problems for viscoelastic thin-walled structures are standard ones. By choosing the basic functions that satisfy certain boundary conditions, the original boundary value problem is reduced to the problem of oscillations of a system with a finite number of degrees of freedom, i.e. to a system of linear (or nonlinear) IDE with one independent variable (of time). As a rule, trigonometric functions or beam functions are used as basic ones. It is known that the choice of such basic functions in many cases limits the class of solved problems to structures of the simplest configurations - beams of constant sections, rectangular plates, cylindrical shells [4].
For elastic structures of a more complex geometric form, it is very tempting to use as basic functions
<nm (x,y), <t>„m (x, y), ¥nm (x, y) , W0nm (x, y) the eigenmodes of oscillations, which are, in a sense,
"native" and allow better consideration of all the features of these structures. For the first time, the expansion of the solution in terms of eigenmodes of oscillations in solving problems of described class was used in [29] (by V.N. Chelomei). Subsequently, the eigenmodes of oscillations in the process of expansion of the sought for solution in the problem of forced oscillations of real elastic and viscoelastic structures (including the shell type) were used in [30] (A.N. Ishmatov, M.M. Mirsaidov) and [31] (T.Z. Sultanov, D.A. Khodzhaev, M.M. Mirsaidov).
In [32] (M.M. Mirsaidov, Ya. Mekhmonov), along with expansion in terms of eigenmodes of oscillations, the problem of vibrations of three-dimensional bodies with associated mass was solved. As is well known, the determination of eigenmodes of oscillations for any structures is a completely independent and a very difficult task. That is why quite often a part of a complete system of known functions in analytical form, that satisfies a particular boundary value problem, is used as a basic function.
Substituting (2) into the system (1) written in dimensionless parameters, and following the Bubnov-Galerkin procedure for the functions unm(t) , vnm(t), Wnm(t) to be determined, a system of nonlinear IDE, i.e. the basic resolving nonlinear IDE are obtained in the form:
NM / \i NM/
Z Z aklnmunm - П1 - Г* J Z Z (
NM
\d1klnmUnm + ^e1klnmVnmj+
n=1 m=1
NM / *\j N m Г i
Z Z bklnmV nm - n2 (1 - Г Л Z ZI ~Jd2klnmUnm + e2klnmVnm
n=1 m=1 ln=1 m=1^
Ln=1 m=1
1 NM
+ T Z Z Sknmy - WbW,
n,i=1 m, j=1
NM
= o,
л
+
N M
Z Z
n=1 m=1
1 N M ( "
+ T Z Z g2klnmij \WnmWj - W0nmW0j,
n,i=1 m, j=1
NM
= o,
NM
CklnmWnm +
Пз ZZ Pklnm (1 - 2Mklnm cos ®t )Wnm -ПзГ* Z Z f
Пз
n=1 m=1
N M
n=1 m=1
NM
Z Z Wnm (1 - Г* hd4klnmijUy + e4klnmijVij J +
n,i=1 m, j=1
N M
r W ij rs
= 12%
4kl ■
U„
NM
+ Z Z gklnmijrsWnm t1 - Г M
n,i,r=1m, j,s=1
(o) = u0 , U (oJ = U0 , v (oJ = v0 , V (oJ = V0 , w (oJ =
j\ / 0nm ' nm\ / 0nm ' nm\ / 0nm ' nm\ / 0nm ' nm\ /
(1 -m2
wn
Wnm (0J = W0nm .
1 1
where aklnm = {j 4nmkldxdy ■ 0 0
(3)
3klnmWnm
1klnm
1 klnm
, xx + hхфпш, x y Wkldxdy ,
00 2 2 У
j jf ^^'(«m, xy + Mh'x(nm,y + ^^''yWnm, x \фkldxdУ, nnV 2 2 У
0 0 1 1
ЯГ 1, , , лЛ+mJ, , „ 1 - m ^ ,
I "j'^nm,xVij,xx + 2 h^nm,yVij,xy + m,xWij,yy +
00
2
2
+ 21'hx ^'nm^jx +Mrh'x ^my^jy + Щу'гШУ/У ^l^y ■ 1 1
bklnm = j j hVnmVHdxdy>
00
d
2klnm
e
2klnm
1 1
0 0 11/
jj(
00
, xy + _ hx^Pnm, y + ßMyfinm,x \(Pkldxdy ,
1 -M
1 -M 1 ,,„,
2Я
2 hCm,xx + h(nm,yy ^ hWnm,x + h'y('nm,y (kldxdy ■
2 Я
11 i + ^ 1 _ ^ Slklnmij = \\(hW'nm ,yyij,yy + 2 hynmyij,xy + 2 hynm , yyij,xx + 00 2Я 2Я
+ \_jr К Ум, У,, y + 2 h'y W'rmyyW'j y + К V'nm, Wj x ■
9 тг2 72 r>
cklnm = ЦhVnmVkldxdy. Pklnm = fldnm _ Pklnm»0 ■ ßklnm =-Г"^ ■
0 0 Pk ln m
fknm = i/К" WU + 2^ + ^,yyyy)+
0 0
h(hx )2+h 2hi y^xx+¿VW^)+6h 2hx (yixxx+¿v™ )-
'2h(h'x ) + h2hxx jyym.xx ) + 6h2hx W,.,,, + ) +
+ 6h2hy (Ууу + ¿>1^)+ "[2h(hy )2 + h2h; ^^y + Zy^)+
.27
+ 6(1 _ [2hhx hy + h% ymxy }ykldxdy , 11 1 _
d4klnmij
00 2 + ^ h'Jgy'J', + + j + XSMy'J' +
11
j -у я »^m, xV'ij,xy + я Ö/Jh'x V'nmMj, y + "y Я »'y W'nm,xtij,x +
0 0
+^ яйК Knm y^j,x + Л4Я Km j + (1 _ ж» упплу^;,, Kidxdy,
11 1 _
Sklnmijrs = 12//ynm,xyij,xyrs,xx + _ Я hynm,xyij ,xyr s.yy + 0 0 2
1 + ß l2 7 f f ff 1 7 ff f f я ß 7 ff f f
+ XhVnm, xWV,yW y + 2 hKnm, У Ц, WrS, x + ~ hxWnm, xWj, y¥rs,y +
1 Д i f i2 f ff 1 ^ » f f я ß i ff ff
+ hy£ Wnmx Wjx Wrsy + 2 h Wn m , xWij, xWrs, x +~ hW nm, xWjyWr^ y +
+ Zh W'nm,yW'iW'rs,y + h W'nmWj ,xW'rs,y + h W'nmWj ,W'rS,x +
Я 7 f f f f Я ß 7 f f f f 1 Д л2 1 f f t t
+ Y hWnn,WijWr,, y hyW nm, Wn, xWrs, x + Я hWnm,yWj, Wrs, y +
+ у Wnm, W'jW'rsy +— h¥'L,yy¥'j W'rs, x +
+(1 _ y'myj,y yy ]vudxdy ■
1 1 352
4kt = qjjwkidxdy. n = -^i.
пЯ2
3T
П
2
Twice integrating the obtained system (3) with respect to time t, it can be written in integral form, as shown in the articles by A.F. Verlan, R.A. Abdikarimov, H. Eshmatov [33].
Assuming that t = ti, ti = iAt, i = 1,2,... (here At is the step of numerical integration) and replacing the integrals with quadrature formulas of trapezoid for calculating the unknowns
uinm = uinm (ti) , vinm = vinm (ti) , winm = winm (ti) , the following recurrence formulas are obtained for the
Koltunov-Rzhanitsyn kernel r(t) = Ae• ta-1, A > 0, J3 > 0 , 0 < a < 1, described in the monograph by M.A. Koltunov, A.S. Kravchuk, V.P. Mayboroda [34] and in A.R. Rzhanitsyn's monograph [35].
N M
N M
p-1
NM
Z ZaklnmUpnm = Z Zaklnm(u0nm + U0nmtpJ+n1 Z Aq(tP -tqJ1 Z Zd1klnm n=1m=1 n=1m=1 q=0 In=1m=1
qnm
-A Z »
a
e~ptzU
z q-z,nm
z=0
NM Г
^ZZ
n=1m=1
1klnm
qnm
Z »
a
e zV z q-z,nm
z=0
1 N M
+ T Z Z g1kl)
n,i=1 m, j=1
NM, . p-1 . . i 1 N mI
H Z Aq Ь - tq Jl T ZZ d
wmmw„y — Z»ze ßzw„mw„
qnm qij
a
z=0
NM
Z Z bklnmVpnm - Z Z bklnm (l
q- z,nm q - z,ij
/I
NM
n=1m=1
klnm VV0nm + V0nmtp p-1
q p q'I я
n=1m=1 q=0 [ n=1m=1
2klnm
qnm
Z »
a n
e H zu z^ "q - z,nm
p-1 / \I1 NM + П2 Z Aqltp - tq J1tZZd2klnm
q=0 Iя n=1m =1
qnm
~Z B*~ßzu
V z=0
A
+
a
z^ "q-z,nm
z=0
/
NM Г
+
ZZ
2klnm
N M
Z Z
n=1m=1
CklnmW pnm
n=1m=1
1 N M
+ T Z Z S 2kl
n,i=1 m, j=1
NM
klnm
v_ ^ZB7e~ßzV.
qnm
a
z q - z,nm
z=0
A
/
Wnnmwnij--Z»ze ßzwq-znmwq
qnm qij
a
q-z,nm q - z,ij
z=0
NM
Z Z Cklnm (w0nm + W0nmtp J- Пз Z Z p^m (1 - 2Mk lnm cos 0tК
n=1m=1 p-1
N M
Z Aq (tp - tq J|ZZ f
q=0
3klnm
I n=1m=1
W,
qnm
n=1m=1
- A Z »
a
-ßtz e w
z^ "q-z,nm
z=0
N M
N M
-Z Zd
n,i=1 m, j =1
N M
-ZZ Ze
n,i=1 m, j =1
Г
4klnmijWqnm
A q U„H - A Z ».
qj
a
ze ZUq - z,ij
z=0
4klnmijWqnm
Vqj ~Z »S^V
a
Jz~ q - z,ij
z=0
- Z Z Ski,
n,i,r=1m, j,s=1
'nmijrsW qnm
A
w„uw„„--z »ze ßzw„^ ;,-w.
qij qrs
a
q - z,ij q - z,rs
z=0
-121
(1 -M2 )
яЯчи
U loJ = U , U loJ = U , V (o) = V , V loJ = Vn
nm \ / 0n^ nm\ / 0n^ nm\ / 0n^ nm\ / 0
.lo) = Wn
W (o J = Wn , W 1
nm\ / 0 nm nm '
0nm "
(4)
1
0 0
+
+
+
+
Here Aq, Bz - are numerical coefficients that do not depend on the choice of integrands and could be of different values depending on the quadrature formulas used.
3. Results and Discussion
Results of computational experiments corresponding to different physical and geometric parameters of a viscoelastic isotropic plate are illustrated by the graphs shown in Figures 3-7. The dependence of the change in thickness is taken in the form
h( x) =1 hG (l -а* x ).
2
Here h0 = h(0) = const, a* - is a parameter characterizing the change in thickness. Note that this law leads to a liniar reduction in thickness of structural element in direction of Ox axis (Figure 2).
a)
0.1
0.2
sap
о 4
0.S
0.6
0.7
0.8
0.9
Figure 2. Changes in plate thickness depending on a* parameter:
a) a" = 0.2 ; b) a" = 0.5
Figures 3 and 4 illustrate the effect of viscoelastic properties of material on dynamic behavior of the
plate.
For example, Figure 3 shows the time t variation of displacements of the midpoint ( x = 0.5 , y = 0.5 ) for different values of rheological parameter A of the relaxation kernel of material. An analysis of results has shown that in all the cases considered, at initial time, the oscillations occur near the equilibrium (initial) state. Over time, the amplitude slowly increases and reaches the value comparable with the one of plate thickness.
Curve 1 ( A = 0 ) corresponds to an elastic solution (without considering viscoelastic properties of material); here the oscillations occur without attenuation. Note that in this case the amplitude of oscillations does not increase, but the beating occurs, which is explained by the proximity of the frequency of natural oscillations of the plate to the frequency of periodic effects.
Curve 2 ( A = 0.5 ) and curve 3 ( A = 0.1 ) are obtained considering the viscoelastic properties of material. It is seen that damped oscillations occur near initial oscillations, but over time, oscillations occur with practically constant amplitude, i.e. become steady. In this case, there is some deviation of oscillations from the neutral axis, which is explained by the presence of initial imperfections in the shape.
Figure 3. Graphs of variation with time t of displacements of the midpoint (x = 0.5 , y = 0.5) of the plate at different values of the relaxation kernel parameter A :
a = 0.25; ( = 0.05 ; ju = 0.3 ;S = 25 ; w0 = 0.01;q = 0 ; 2 = 1 ;a* = 0.5 ;S0 = 0;S1 = 0.5 ; ©1 = 1.1;
A = 0.0 (1); 0.05 (2); 0.1 (3)
Figure 4 shows the graphs of deflection versus time t for the midpoint (x = 0.5 , y = 0.5) of the plate corresponding to different values of the relaxation kernel parameter: at a = 0.1 (curve 1); at a = 0.25 (curve 2); at a = 0.5 (curve 3).
In all considered cases, oscillations occur near the equilibrium state in the proximity of initial time. Over time, the amplitude slowly increases, and reaches the thickness of the plate. Then the pattern described is periodically repeated at values of a = 0.25; 0.5 . At sufficiently large values of a and time,
the amplitude of oscillations decreases insignificantly. The amplitude of oscillations at a = 0.1 at large values of time increases, which is explained by the proximity to the resonance phenomenon, since the difference between the eigenfrequency of vibrations of the plate and the frequency of excitations becomes smaller, which completely corresponds to results described in the monograph by A.E. Bogdanovich [5].
0 10 3t> 30 40 50 £0 70 80 90 t
Figure 4. Graphs of dependence of deflection of the midpoint (x = 0.5 , y = 0.5) of the plate on time t at different values of the relaxation kernel parameter a :
A = 0.01; ( = 0.05 ; ju = 0.3 ;S = 25 ; w0 = 0.01;q = 0; 2 = 1 ;a* = 0.5 ;S0 = 0;S1 = 0.5 ; ©1 = 1.1;
a = 0.1 (1); 0.25 (2); 0.5 (3)
Figure 5 shows the graphs of deflection versus time t of the midpoint (x = 0.5 , y = 0.5) of the plate for different aspect ratios: at 2 = 1.0 (curve 1); at 1.5 (curve 2); at 2 (curve 3). Results of calculations
show that for a square plate (with aspect ratio A = 1.0 ), and for a rectangular plate with A = 1.5 at initial time, oscillations occur near the equilibrium state and, over time, approach harmonic oscillations. At A = 2.0, oscillations at initial period of time occur near initial state W0, and over time, oscillation
amplitudes increase sharply, which is explained by the approach of eigenfrequency of the plate to the frequency of excitation force.
Л /\ 2 1 1
3 v
0 20 40 60 80 100
Figure 5. Graphs of dependence of deflection on time t of the midpoint ( x = 0.5 , y = 0.5 ) of the
plate at different aspect ratios A :
A = 0.05 ; a = 0.25 ; ( = 0.05 ; j = 0.3 ; 5 = 25 ; w0 = 0.001 ; q = 0 ; a = 0.5 ; 50 = 0 ; 51 = 0.5 ;
©1 = 1.1 ; A = 1 (1); 1.5 (2); 2 (3)
Figure 6 shows the graphs of deflection versus time t for the midpoint ( x = 0.5 , y = 0.5 ) of the
*
plate at various values of its thickness parameter a , taking into account viscoelastic properties of
material. Namely: at a = 0 (curve 1); at 0.4 (curve 2); at 0.8 (curve 3). Note that an increase in the
*
parameter of a according to the above law corresponds to a decrease in thickness of the plate along the
length of Ox axis. This causes the resonance regime and an increase in oscillation amplitude of the plate
*
(curve 3). At other values of parameter a (curves 1 and 2), the oscillations become steady.
m io JÏ É 40 so a as » щ юо по ж /
Figure 6. Graphs of dependence of deflection on time t of the midpoint (x = 0.5 , y = 0.5) of the
*
plate at various values of the parameter of its thickness a : A = 0.05; a = 0.25; в = 0.05 ;
ju = 0.3 ; S = 25 ; 1 = 1; w0 = 0.01; q = 0 ;£0 = 0 ;S1 = 0.5 ; ©j = 1.1; a* = 0 (1); 0.4 (2); 0.8 (3)
Figure 7 presents results of studies on vibrations of a viscoelastic plate at various values of external constant loads: at q = 0.2 (curve 1); at q = 0.5 (curve 2); at q = 1 (curve 3). It is seen that as external load increases, oscillations occur near the equilibrium state of the plate (Figure 7). At initial period of time,
the nonlinear parametric oscillations retain the beat character, with the passage of time they become less pronounced, and approach stationary ones. Analogous qualitative results are described in Bogdanovich's monograph [5].
0.15-
0.05-
-0.05- -
200
Figure 7. Graphs of dependence of deflection on time t for the midpoint ( x = 0.5 , y = 0.5 ) of the plate at different values of external load q :
A = 0.05 ; a = 0.25 ; p = 0.05 ; ju = 0.3 ; 5 = 25 ; 2 = 1 ; w0 = 0.01 ; a* = 0.5 ; 50 = 0 ; 5 = 0.5 ;
©! = 1.1 ; q = 0.2 (1); 0.5 (2); 1 (3)
4. Conclusions
1. Mathematical model, technique for its investigation and computational algorithm for the problem of parametric nonlinear oscillations of plates of variable thickness are developed taking into account viscoelastic properties of material.
2. An appropriate computer software package is designed to effectively evaluate dynamic instability of plates in question, depending on mechanical characteristics of their material, external excitation loads, geometric parameters, eigenfrequency of vibration, and excitation forces.
3. As a result of numerical experiments it was established that:
a. the presence of viscoelastic properties of material leads to a steady-state oscillation regime even in resonance modes;
b. the presence of initial imperfections of the form leads to a deviation from the reference plane, which is assumed to be neutral (due to the deviation of the latter);
c. a decrease in thickness of the plate along its length in some cases can lead to a resonance mode of oscillations.
Besides, the dependences of some parameters of stable oscillatory process of the plates on the amplitude and frequency of external transverse load are obtained.
References
1. Belyaev, N.M. Ustoychivost prizmaticheskikh sterzhney pod deystviyem peremennykh prodolnykh sil [Stability of Prismatic Rods under the Influence of Variable Longitudinal Forces]. Engineering Constructions and Structural Mechanics. 1924. P. 149-167. (rus)
2. Bolotin, V.V. Dinamicheskaya ustoychivost uprugikh sistem [Dynamic Stability of Elastic Systems]. Moscow: Gostekhizdat. 1956. 600 p. (rus)
3. Vlasov, V.Z. Obshchaya teoriya obolochek i yeye prilozheniya v tekhnike [General Theory of Shells and its Applications in Engineering]. Moscow-Leningrad: Gostekhizdat. 1949. 784 p. (rus)
Литература
1. Беляев Н.М. Устойчивость призматических стержней под действием переменных продольных сил // Инженерные сооружения и строительная механика. 1924. С. 149-167.
2. Болотин В.В. Динамическая устойчивость упругих систем. М.: Гозтехиздат, 1956. 600 с.
3. Власов В.З. Общая теория оболочек и ее приложения в технике. М.-Л.: Гостехиздат. 1949. 784 с.
4. Вольмир А.С. Нелинейная динамика пластинок и оболочек. М.: Наука. 1972. 432 с.
5. Богданович А.Е. Нелинейные задачи динамики цилиндрических композитных оболочек. Рига: Зинатне, 1987. 295 с.
4. Volmir, A.S. Nelineynaya dinamika plastinok i obolochek [Nonlinear Dynamics of Plates and Shells]. Moscow: Science, 1972. 432 p. (rus)
5. Bogdanovich, A.E. Nelineynyye zadachi dinamiki tsilindricheskikh kompozitnykh obolochek [Nonlinear Problems of the Dynamics of Cylindrical Composite Shells]. Riga: Zinatne, 1987. 295 p. (rus)
6. Schmidt, T. Parametricheskiye kolebaniya [Parametric Oscillations]. Moscow: Mir, 1978. 336 p. (rus)
7. Kalinnikov, A.E., Kuznetsov, N.P. Chislennoye issledovaniye parametricheskikh kolebaniy sterzhnya [Numerical Investigation of Parametrical Vibrations of Beam]. Vestnik M.T. Kalashnikov IZHGTU. 2016. Vol. 69. No. 1. Pp. 3-6. (rus)
8. Yan, Q., Ding, H., Chen, L. Nonlinear dynamics of axially moving viscoelastic timoshenko beam under parametric and external excitations. Applied Mathematics and Mechanics. 2015. Vol. 36. No. 8. Pp. 971-984.
9. Akulenko, L.D., Kumakshev, S.A., Nesterov, S.V. Chislenno-analiticheskiy metod issledovaniya parametricheskikh kolebaniy [A Numerical-analytical Method for Investigating Parametric Oscillations]. Journal of Applied Mathematics and Mechanics. 2015. Vol. 79. No. 2. P. 111-121. (rus)
10. Karpov, V.V., Ignat'ev, O.V., Semenov A.A. The stressstrain state of ribbed shell structures. Magazine of Civil Engineering. 2017. No. 6. Pp. 147-160.
11. Kurpa, L., Mazur, O.S., Tkachenko, Ya.V. Parametric Vibration of Multilayer Plates of Complex Shape. Journal of Mathematical Sciences. 2014. Vol. 203. No. 2. Pp. 165-184.
12. Darabi, M., Ganesan, R. Nonlinear dynamic instability analysis of laminated composite thin plates subjected to periodic in-plane loads. Nonlinear Dynamics. 2018. Vol. 91. No. 1. Pp. 187-215.
13. Huynh, H.Q., Nguyen, H., Nguyen, H.L.T. Non-linear parametric vibration and dynamic instability of laminated composite plates using extended dynamic stiffness method. Journal of Engineering Technology. 2017. Vol. 6. Pp. 170-185.
14. Kumar, R., Dutta, S.C., Panda S.K. linear and non-linear dynamic instability of functionally graded plate subjected to non-uniform loading. Composite Structures. 2016. Vol. 154. Pp. 219-230.
15. Kumar, R., Mondal, S., Guchhait, Sh., Jamatia, R. Analytical approach for dynamic instability analysis of functionally graded skew plate under periodic axial compression. International Journal of Mechanical Sciences. 2017. Vol. 130. Pp. 41-51.
16. Ievzerov, I.D. The stability problems for bars and plates. Magazine of Civil Engineering. 2014. No. 1. Pp. 6-11. doi: 10.5862/MCE.45.2 (rus)
17. Dubrovin, V.M., Butina, Т.А. Modelirovaniye dinamicheskoy ustoychivosti tsilindricheskoy obolochki pri tsiklicheskom osevom vozdeystvii [Simulation of Dynamic Stability of a Cylindrical Shell under Cyclic Axial Impact]. Mathematical Modeling and Computational Methods. 2016. No. 3(11). Pp. 24-32. (rus)
18. Mochalin, A.A. Parametricheskiye kolebaniya neodnorodnoy krugovoy tsilindricheskoy obolochki peremennoy plotnosti pri razlichnykh krayevykh usloviyakh [The Parametric Oscillations of Heterogeneous Circular Cylindrical Shell of Variable Density under Different Boundary Conditions]. Izv. Saratov Univ. Ser. Math. Mech. Inform. 2015. Vol. 15. No. 2. Pp. 210-215. (rus)
19. Usarov, M.K. [Buckling of orthotropic plates with bimoments. Magazine of Civil Engineering. 2015. No. 1. Pp. 80-90. (rus)
20. Dey, T., Ramachandra, L.S. Dynamic stability of simply supported composite cylindrical shells under partial axial
6. Шмидт Т. Параметрические колебания. М.: Мир, 1978. 336 с.
7. Калинников А.Е., Кузнецов Н.П. Численное исследование параметрических колебаний стержня // Вестник ИЖГТУ им. М.Т.Калашникова. 2016. Т. 69. № 1. С. 3-6.
8. Yan Q., Ding H., Chen L. Nonlinear dynamics of axially moving viscoelastic Timoshenko beam under parametric and external excitations // Applied Mathematics and Mechanics. 2015. Vol. 36. № 8. Pp. 971-984.
9. Акуленко Л.Д., Кумакшев С.А., Нестеров С.В. Численно-аналитический метод исследования параметрических колебаний // Прикладная математика и механика. 2015. Т. 79. № 2. С. 163-180.
10. Карпов В.В., Игнатьев О.В., Семенов А.А. Напряженно-деформированное состояние ребристых оболочечных конструкций // Инженерно-строительный журнал. 2017. № 6(74). С. 147-160.
11. Kurpa L., Mazur O.S. Tkachenko Ya.V. Parametric Vibration of Multilayer Plates of Complex Shape // Journal of Mathematical Sciences. 2014. Vol. 203. № 2. Pp. 165-184.
12. Darabi M., Ganesan R. Nonlinear dynamic instability analysis of laminated composite thin plates subjected to periodic in-plane loads // Nonlinear Dynamics. 2018. Vol. 91. № 1. Pp. 187-215.
13. Huynh H.Q., Nguyen H., Nguyen H.L.T. Non-linear parametric vibration and dynamic instability of laminated composite plates using extended dynamic stiffness method // Journal of Engineering Technology. 2017. Vol. 6. Pp. 170-185.
14. Kumar R., Dutta S.C., Panda S.K. Linear and non-linear dynamic instability of functionally graded plate subjected to non-uniform loading // Composite Structures. 2016. Vol. 154. Pp. 219-230.
15. Kumar R., Mondal S., Guchhait Sh., Jamatia R. Analytical approach for dynamic instability analysis of functionally graded skew plate under periodic axial compression // International Journal of Mechanical Sciences. 2017. Vol. 130. Pp. 41-51.
16. Евзеров И.Д. Задачи устойчивости для стержней и пластин // Инженерно-строительный журнал. 2014. № 1(45). С. 6-11.
17. Дубровин В.М., Бутина Т.А. Моделирование динамической устойчивости цилиндрической оболочки при циклическом осевом воздействии // Математическое моделирование и численные методы. 2016. № 3(11). С. 24-32.
18. Мочалин А.А. Параметрические колебания неоднородной круговой цилиндрической оболочки переменной плотности при различных краевых условиях // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15. № 2. С. 210-215.
19. Усаров М.К. Изгиб ортотропных пластин с учетом бимоментов // Инженерно-строительный журнал. 2015. № 1(53). С. 80-90.
20. Dey T., Ramachandra L.S. Dynamic stability of simply supported composite cylindrical shells under partial axial loading // Journal of Sound and Vibration. 2015. Vol. 353. Pp. 272-291.
21. Budanov V.M. A form of the solution to the Mathieu equation // Moscow University Mechanics Bulletin. 2017. Vol. 72. № 1. Pp. 7-12.
22. Eshmatov B.Kh. Nonlinear vibrations and dynamic stability of viscoelastic orthotropic rectangular plates // Journal Sound and Vibration. 2007. Vol. 300. № 3. Pp. 709-726.
23. Абдикаримов Р.А., Голоскоков Д.П. Численное исследование нелинейных колебаний вязкоупругих
loading. Journal of Sound and Vibration. 2015. Vol. 353. Pp. 272-291.
21. Budanov, V.M. A form of the solution to the Mathieu equation. Moscow University, Mechanics Bulletin. 2017. Vol. 72. No. 1. Pp. 7-12.
22. Eshmatov, B.Kh. Nonlinear vibrations and dynamic stability of viscoelastic orthotropic rectangular plates. Journal Sound and Vibration. 2007. Vol. 300. No. 3. Pp. 709-726.
23. Abdikarimov, R.A., Goloskokov, D.P. Chislennoye issledovaniye nelineynykh kolebaniy vyazkouprugikh plastin peremennoy tolshchiny [Numerical research of nonlinear fluctuations of viscoelastic plates of the variable thickness]. The Journal of University of Water Communications. 2011. No. 2(9). Pp. 102-107. (rus)
24. Abdikarimov, R.A., Khudayarov, B.A. Modelirovaniye kolebatelnykh protsessov vyazkouprugikh ortotropnykh plastin s peremennoy zhestkostyu [Modelling of vibration processes of viscoelastic orthotropic plates with variable rigidity]. Mechanics. Proceedings of National Academy of Sciences of Armenia. 2011. No. 4(64). Pp. 30-38. (rus)
25. Eshmatov, B.Kh. Dinamicheskaya ustoychivost vyazkouprugikh plastin pri vozrastayushchikh szhimayushchikh nagruzkakh [Dynamic stability of viscoelastic plates under increasing compressing loads]. Journal of Applied Mechanics and Technical Physics. 2006. No. 2(47). Pp. 289-297. (rus)
26. Eshmatov, B.Kh. Nonlinear vibration analysis of viscoelastic plates based on a refined Timoshenko theory. Journal International Applies Mechanics. 2006. Vol. 42. No. 5. Pp. 596-605.
27. Eshmatov, B.Kh., Eshmatov, Kh., Khodzhaev, D.A. Nelineynyy flatter vyazkouprugikh pryamougolnykh plastin i tsilindricheskikh paneley iz kompozitsionnogo materiala s sosredotochennymi massami [Nonlinear flutter of viscoelastic rectangular plates and cylindrical panels of a composite with a concentrated masses]. Journal of Applied Mechanics and Technical Physics. 2013. No. 4(54). Pp. 578-587. (rus)
28. Abdikarimov, R.A., Zhgutov, V.M. Matematicheskiye modeli zadach nelineynoy dinamiki vyazkouprugikh ortotropnykh plastin i obolochek peremennoy tolshchiny [Mathematical models of problems of nonlinear dynamics of viscoelastic orthotropic plates and shells of variable thickness]. Magazine of Civil Engineering. 2010. No. 6. Pp. 38-47. (rus)
29. Chelomei, V.N. Dinamicheskaya ustoychivost elementov aviatsionnykh konstruktsiy [Dynamic stability of the elements of aircraft structures]. Moscow: Aeroflot, 1939. 80 p. (rus)
30. Ishmatov, A.N., Mirsaidov, M. Nonlinear Vibrations of an Axisymmetric Body Acted upon by Pulse Loads. Soviet Applied Mechanics. 1991. Vol. 27. No. 4. Pp. 68-74.
31. Sultanov, T.Z., Khodzhaev, D.A., Mirsaidov, M.M. Assessment of dynamic behavior of heterogeneous systems taking into account non-linear viscoelastic properties of soil. Magazine of Civil Engineering. 2014. No. 1(45). Pp. 80-89.
32. Mirsaidov, M., Mekhmonov, Ya. Non-axisymmetric vibrations of axisymmetric structures with associated mass and hollows (protrusions). J. Strength of Materials. 1987. No. 3. Pp. 111-116.
33. Verlan, A.F., Abdikarimov, R.A., Eshmatov, H. Chislennoye modelirovaniye nelineynykh zadach dinamiki vyazkouprugikh sistem s peremennoy zhestkostyu [Numerical modeling of nonlinear problems of the dynamics of viscoelastic systems with variable rig id it]. Electronic modeling. 2010. Vol. 32. No. 2. Pp. 3-14.
34. Koltunov, M.A., Kravchuk, A.S., Mayboroda, V.P. Prikladnaya mekhanika deformiruyemogo tverdogo tela
пластин переменной толщины // Журнал университета водных коммуникаций. 2011. № 2(9). С. 102-107.
24. Абдикаримов Р.А., Худаяров Б.А. Моделирование колебательных процессов вязкоупругих ортотропных пластин с переменной жёсткостью // Известия НАН Армении. Механика. 2011. № 4(64). С. 30-38.
25. Эшматов Б.Х. Динамическая устойчивость вязкоупругих пластин при возрастающих сжимающих нагрузках // Прикладная механика и техническая физика. 2006. Т. 47. № 2. С. 165-175.
26. Eshmatov B.Kh. Nonlinear vibration analysis of viscoelastic plates based on a refined Timoshenko theory // Journal International Applied Mechanics. 2006. Vol. 42. № 5. Pp. 596-605.
27. Эшматов Б.Х., Эшматов Х., Ходжаев Д.А. Нелинейный флаттер вязкоупругих прямоугольных пластин и цилиндрических панелей из композиционного материала с сосредоточенными массами // Прикладная механика и техническая физика. 2013. Т. 54. № 4. С. 7485.
28. Абдикаримов Р.А., Жгутов В.М. Математические модели задач нелинейной динамики вязкоупругих ортотропных пластин и оболочек переменной толщины // Инженерно-строительный журнал. 2010. № 6. С. 38-47.
29. Челомей В.Н. Динамическая устойчивость элементов авиационных конструкций. М: Аэрофлот, 1939. 80 c.
30. Ishmatov A.N., Mirsaidov M. Nonlinear vibrations of an axisymmetric body acted upon by pulse loads // Soviet Applied Mechanics. 1991. Vol. 27. № 4. Pp. 68-74.
31. Султанов Т.З., Ходжаев Д.А., Мирсаидов М.М. Оценка динамического поведения неоднородных систем с учетом нелинейно-вязкоупругих свойств грунта // Инженерно-строительный журнал. 2014. № 1(45). С. 80-89
32. Mirsaidov M., Mekhmonov Ya. Non-axisymmetric vibrations of axisymmetric structures with associated masses and hollows (protrusions) // Journal Strength of Materials. 1987. № 3. Pp. 111-116.
33. Верлань А.Ф., Абдикаримов Р.А., Эшматов Х. Численное моделирование нелинейных задач динамики вязкоупругих систем с переменной жесткостью // Электронное моделирование. 2010. Т .32. № 2. С. 3-14.
34. Колтунов М.А., Кравчук А.С., Майборода В.П. Прикладная механика деформируемого твердого тела. М.: Высшая школа, 1983. 349 с.
35. Ржаницын А.Р. Теория ползучести. М.: Стройиздат, 1968. 416 с.
[Applied mechanics of a deformable solid]. Moscow: High School. 1983. 349 p.
35. Rzhanitsyn, A.R. Teoriya polzuchesti [Creep Theory]. Moscow: Stroiizdat, 1968. 416 p.
Mirziyod Mirsaidov, +7(987)237-09-81; theormir@mail.ru
Rustamkhan Abdikarimov*, +7(99890)928-45-54; rabdikarimov@mail.ru
Nikolai Vatin,
+7(921)964-37-62; vatin@mail.ru Vladimir Zhgutov,
+7(931)532-25-22; abc.kitezh@gmail.com
Dadakhan Khodzhaev, +7(99871)237-09-93; khodzhaevda@mail.ru
Bakhodir Normuminov,
+9(9898)307-50-65; bnormuminov1977@mail.ru
Мирзиед Мирсаидович Мирсаидов, +7(987)237-09-81; эл. почта: theormir@mail.ru
Рустамхан Алимханович Абдикаримов*,
+7(99890)928-45-54;
эл. почта: rabdikarimov@mail.ru
Николай Иванович Ватин, +7(921)964-37-62; эл. почта: vatin@mail.ru
Владимир Михайлович Жгутов,
+7(931)532-25-22;
эл. почта: abc.kitezh@gmail.com
Дадахан Акмарханович Ходжаев,
+7(99871)237-09-93;
эл. почта: khodzhaevda@mail.ru
Баходир Ашурович Нормуминов,
+9(9898)307-50-65;
эл. почта: bnormuminov1977@mail.ru
© Mirsaidov, M.M.,Abdikarimov, R.A.,Vatin, N.I.,Zhgutov, V.M.,Khodzhaev, D.A.,Normuminov, B.A., 2018
