Low-Frequency Vibration Modes of Strongly Inhomogeneous
Elastic Laminates
J. D. Kaplunov, L. A. Prikazchikova
Julius D. Kaplunov, https://orcid.org/0000-0001-7505-4546, School of Computing and Mathematics, Keele University, Keele, Staffordshire, ST5 5BG, UK, j.kaplunov@keele.ac.uk
Ludmila A. Prikazchikova, https://orcid.org/0000-0001-9051-2103, School of Computing and Mathematics, Keele University, Keele, Staffordshire, ST5 5BG, UK, l.prikazchikova@keele.ac.uk
The dynamic behaviour of thin multi-layered structures, composed of contrasting "strong" and "weak" layers, is considered. An asymptotic procedure for analysing the lowest cutoffs is developed. A polynomial frequency equation is derived, along with the linear equations for the associated eigenforms corresponding to displacement variation across the thickness. For a five-layered laminate with clamped faces two term expansions for eigenfrequencies and eigenforms are compared with those obtained from the exact solution of the original problem for thickness resonances.
Key words: asymptotic, high-contrast, thin elastic structure, thickness fluctuations, low frequency, cutoff. DOI: https://doi.org/10.18500/1816-9791-2018-18-4-447-457
INTRODUCTION
The classical theory for thin homogeneous elastic structures governing bending, extensional and torsional-shear vibration modes, e.g. see [1-3], is asymptotically consistent over the low-frequency range. Obviously, it does not take into consideration any of high-frequency modes with the cut-off frequencies related to thickness resonances. At the same time, currently available refined formulations incorporating the lowest high-frequency modes ( [4, 5] and reference therein) do not appear to be uniformly asymptotically valid, as noted for example in [6]. These observations are also true for layered structures not demonstrating a substantial contrast in material and geometric properties of the layers, e.g. [7]. However, multi-parametric analysis of strongly vertical inhomogeneous three-layered plates in [6], see also [8] and [9], indicates that for certain combination of problem parameters the lowest thickness shear resonance frequency becomes asymptotically small, resulting in an extra low-frequency vibration mode in comparison with the traditional layout. Four scenarios studied in [6] correspond to the practically important setups of sandwich structures, laminated glass, photovoltaic panels and precipitator plates utilised in gas filters, see [10-13].
In this paper we extend the framework of the cited publication [6] on thin strongly inhomogeneous plates to multi-layered structures with arbitrary number of layers. For the sake of definiteness, the upper and lower faces are supposed to be clamped. For homogeneous plates, such boundary conditions do not support low-frequency vibrations [14,15]. Below we consider layers of two different types, namely, "strong" and "weak" ones, adapting the ratio of their stiffnesses and densities as a single small parameter, such that the wave speeds in strong and weak components are of the same order. In addition, we do not impose any special restriction on the thicknesses of the layers, apart from the assumption that all of them are also of the same asymptotic order. The consideration is restricted to the evaluation of thickness resonances.
The paper is organised as follows. The problem is formulated in Section 1. The general asymptotic procedure is developed in Section 2. A polynomial equation is derived for the
leading order values of the thickness resonances along with linear algebraic equations for the associated eigenforms. As might be expected, the order of the aforementioned equation coincides with the number of strong layers undergoing almost rigid body motion [9]. In this case the weaker layers are subject to nearly homogeneous thickness deformation. The results of Section 2 are specified in the next Section 3 for a five-layered laminate. A correction to the leading order approximation is also calculated, in order to illustrate small deviations from rigid body motions and homogeneous deformations. Numerical results are discussed in the last section.
1. STATEMENT OF THE PROBLEM
Consider a thin elastic laminate composed of n alternating strong and weak layers of thickness hi, i = 1,...,n, see Fig. 1. Denote Young moduli, Poisson ratios and densities of the layers by Em, vm and pm with m = 1 and m = 2 corresponding to weak and strong layers, respectively. Thus, Ei ^ E2 and also, for the sake of definiteness, we
assume that E1 /E2 ~ p1 /p2.
The consideration below is oriented to harmonic vibrations with angular frequency u of a laminate with fixed faces over the low frequency range
hn
hn-i
h4
h3
h2
hi
и < -1. h,
Em Pm
i = 1,..., n, m = 1, 2.
Fig. 1. Multi-layered laminate with n alternating weak and strong layers
We start from the equations [2]
(i)
The main focus is on the effect of elastic contrast, since low frequency vibrations is not a feature of homogeneous structures with fixed faces, see [14]. We restrict ourselves to the evaluation of the cutoffs of low-frequency modes arising near the thickness resonances of a flat plate clamped along the faces of weak outer layers.
Em
d2 v
2(1 + Vm)Km dz2
3, , 2
+ Pm U V3i
0,
1, .. , n,
m
1, 2,
where
and
Km
'1 - 2vr
2- 2v
Em
d2Vfc, 2
+ PmU Vki
0,
k = 1, 2,
2(1 + Vm ) dz2
for stretch and shear thickness vibrations, respectively, with vki and v3i denoting tangential and transverse displacements. The continuity conditions along interfaces are written as
V3i = V3j,
Ei
dv
3i
E2
dv3j
2(1 + v1 )k1 dz 2(1 + v2)k2 d;
or
vki = vkj, k = 1, 2,
Ei dvki = E2 dvkj 2(1 + v1) dz = 2(1 + v2) dz '
where j = i + 1 for a "weak-strong" interface and j = i — 1 for a "strong-weak" one. In addition, along fixed faces z = 0 and z = hi + h2 + ... + hn
vki = V31 = vkn = V3n = 0, k = 1, 2. The equations above can be re-written in the form
d2 m
Dm^—7 + Pm V2 Ui = 0, (2)
dz2
with
Ui = Uj ,
D dui D duj (3) dz dz
and
u1 = un = 0 (4)
E
at z = 0 and z = h1 + h2 + ... + hn, where ui = v3i and Dm = —-^- for stretch
2(1 + Vm )Km
Ei
vibrations and ui = vki and Dm = ——i—- for shear vibrations.
2(1 + Vi)
The formulated boundary value problem can be treated asymptotically over low-frequency range (1), due to the contrast properties of the layers.
2. ASYMPTOTIC PROCEDURE
First, re-write the equations in the previous section in dimensionless local coordinates Zi = z/hi and frequencies Qi = vhi/cm with cm = \JDm/pm, i = 1, 2,n and m = 1,2, having
d2 Ui
dZ2 u = 0. (5)
The frequency parameters corresponding to both weak or strong i-th and j-th layers are related to each other as
Qi = Lj Qj, (6) where Lj = hi/hj. At the same time, for i-th weak layer and j-th strong one
Qi = cLj Qj, (7) where c = c2/ci. Boundary conditions (4) become
u11 Zi =0 = |Zn =bn+1 = 0, (8)
where
1 i-1
bi =0 and bi = , i = 2,...,n,
t hi
i n=0
such that bi < Zi < bi + 1.
Continuity conditions (3) along interfaces take the form
ui \ Zi =bi+1 = ui+1 \zi+1 =bi+1
and
du,
dZi
= £p L
du
i+1
Zi=bi+1
i+1 dZ
i+1
(9)
(10)
Zi+1 =bi+1
where p = 1 if i-th layer is strong, and p = —1 if it is weak. Here small parameter s is introduced as the ratio
D1
s = D «1
As it has been already mentioned, we also assume that p1 /p2 ~ s. In this case we have 0 ~ 02 ~ ... ~ ~ s over low frequency range (1). Similarly to [9], we expand frequencies and displacements ui in the asymptotic series
and
0 = £ (nj, + en?1 + ...)
Ui = Ui0 + £Ui1 + ...
At leading order, we have from equation (5) for strong layers
d2 Ui0
dZ2 "
subject to the Neumann boundary conditions
0,
(11)
dUiO
dZ,
du
iO
Zi =bi
dZi
= 0.
Zi =bi+1
As a result, we arrive at uniform variation across the thickness
Ui0 = Ci,0 = const,
(12)
corresponding to rigid body motion.
Next, we proceed with equations (11) for weak layers. For inner weak layers we have
UiO |
Zi =bi
Ci-1,0 5 ui0 1
Zi =bi+1
Ci + 1,0 5
while for outer ones
U10 |Zi =0 = Un0 lzn =bn+1 =
Thus, we obtain for eigenforms
U10 = C2,0 Z1,
ui0 = Ci-1,0 + (Ci+i,0 — Ci-1,0) (Zi — bi)
(13)
un0 = 1,0 (bn + 1 — Zn) •
At next order, we start again with the equations for strong layers following from (5). They are
dff + Uo = 0. (14)
The associated boundary conditions, derived from (10) and (12), become
du^i
dZi dui1
dZi
Zi =bi
Zi =bi+1
Li-1 (Ci,Q — 2,0 ),
Li+1 (Ci+2,Q — Ci,Q)-
(15)
The compatibility of equations (14) and boundary conditions (15) results in the relations
C2,0 ^20 = L1 C2,0 — L3 (C4,0 — C2,0) ,
Ci,Q ^¿Q — Li-1 (Ci,Q — 2,0) — Li+1 (Ci+2,Q — Ci,0) ,
(16)
1,0 0 = Ln 1 1,0 + LJJ-2 (Cn-1,0 — 3,0) •
The latter together with formulae Oi0 = LjOj0, see (6), can be used to determine unknown constants Ci,0 and, therefore, leading order eigenfrequencies as will be demonstrated below for a five-layered laminate. In addition, next order corrections to eigenfrequencies and eigenforms will be derived.
3. FIVE-LAYERED LAMINATE
Consider a five-layered laminate with alternating weak and strong layers clamped along the faces of outer weak layers, see Fig. 2.
First, we have from (16) for two strong components
C2,0 °20 = L1 C2,0 — L3 (C4,0 — C2,0) 5 C4 ,0 °40 = L5C4,0 + L4 (C4,0 — C2,0) 5
"ШШШШ^ШШШШШШШШ^
Fig. 2. Five-layered laminate with clamped faces composed of two strong and three weak layers
with 04Q — ¿4 02Q. These equations yield
i4 — (l2)2 (L4 + ¿4 (1 - k)),
together with the relation with parameter k defined as
C2,q — k C4,q,
(17)
(18)
k'1,2 —
-b ±y/b2 + 4L4 (L4)
2L|
b — L4 l4 + L4 l4 - L3 - ¿4.
2
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2018. Т. 18, вып. 4 Then, we obtain from (13) and (18) for leading order eigenforms
uio = 64,0kZi, W20 = 64,0k, M30 = 64,0 (k + (1 - k) (Z3 - 63)), (19)
U40 = 64,0, U50 = 64,0 (1 + 65 — Z5).
Next, we have from boundary value problem (14)-(15) for the first order correction to strong layer displacements
uii = Ai,i Z2 + Bi,i Zi + 6^,1, i = 2,4, (20)
where
A2,i = — 264,0k ^05 B2,i = 64,0k(L2 + ^062),
A4,i = — 264,0^0 (L4) 5 B4,i = 64,0(^20 (L4) (b4 + 1) — L5 )•
For weak layers the first order correction is determined from equations (14) subject to the boundary conditions
uii U =0 = 0 uii U=i = U2i |Z2 =b2 5
u3i |Z3 =b3 = U2i |Z2 =b2+i 5 u3i |Z3=b3+i = U4i U4=64 ' u5i |Z5 =65 = u4i |Z4 =64 + i 5 u5i |Z5 =65 + i = 0
As a result, we arrive at
Uii = Fi,iZ3 + Gi,iZ2 + Hi,iZi + Ki,i5 i = 1 35 55 (21)
where
Fi,i = — 164,0^0 (L2)2 c2k Gi,i = 05
F3,i = 6C4,0^0 (L2)2 c2(k - 1), G3,i = 2C4,0^0 (L2)2 c2 (63(1 - k) - k), F5,1 = ^C4,0^0 (L5)2 c2, = -1 C4,0^0 (L5)2 c2(65 +1).
Then, setting one of the constants equal to zero, namely C4)1 = 0, and applying Dirichlet boundary conditions for weak layers we obtain for the rest of the constants
#1,1 = A2,1 62 + B2,1 62 — F1,1 + C2,1, #3,1 = A4,1 b4 + B4,1 64 - A2,1 (b2 + 1)2 - B2,1 (62 + 1) - Fj,1 (3&2 + 3b3 + 1) -
-G3,1(2b3 + 1) -C2,1, #5,1 = A-4,1 (64 + 1)2 - ^4,1(64 + 1) - F5,1 (362 + 365 + 1) - G5,1 (265 + 1)
and
Ki,i = 0, Кз,1 = A4,162 + B4164 - F3,1 (63 + 1)3 - G3,1 (63 + 1)2 - Нз,1 (63 + 1). K5,1 = -F5,1 (65 + 1)3 - G5,1 (65 + 1)2 - H5,1
where C2,1 is still unknown.
Now, we consider vibrations of strong layers at second order, starting from the equations
dzf
subject to the boundary conditions
Y + fis, Uii + Щ Ui0 = 0, i = 2,4,
(22)
du
i2
dZi
L
dui—
i—11
i-1
Zi =bi
dZi—1
Zi_i =bi_i+1
du,
i2
dZi
L
du,
i+11
i+1
dZ
i+1
(23)
Zi+i =bi+i
Zi =bi+1
The compatibility condition for the last boundary value problem leads to
Q2
-k ( - 3b2 fi20 + 3L1b2 + 3L3 b2 - 3b2 fi20 + 3L3 + 6l2 bs - ^20) --6A4,1 l2b2 + 3B2,1 ( - 2b2fi2o + 2L1 b2 + 2l2b2 - fi2o + 2L3) - 6B4,1 l2b4+
+6С2д ( - fi20 + L1 + L2) + 12Е1ДL2 + 6^L2 (3Ьз + ) + 6G3,1 L2), (24)
where
C * C2,1
C
2,1
1
1
C
4,0
6 L4L4 + L2L3 + L3k - (L4)2^0
4 ) 2
X
x ^ 3L4 L4b2 + 3L4 l4 b2 - 3b2 ^ (l4)2 - 3b2 fi20 (L )2 + -l4 l4 b2 -
-fi20 (l2)2 + 3l3 b2 k + 3l2 l4 + 6l3 b2k + 3l4 k) + 2A*, 1 ( 3b|kfi20 (L )2 + +3b4 k fi20 (L )2 - 3l4 l4 b4 + ш20 (l2)2 - 3l4 b4 k - 3L4b4 k - 6l4 b4k - 3l5 k) -
2b2 fi20 (l4 )2 + 2l4l4 b2 + 6l2 l4 b2 - fi20 (l4)2 + 2l2 l4 + 2l4 b2 k + 2l4 k)+
+3В*д (kfi20 (L)2 + 2b4kfi20 (L4)2 - 2L4L3bA - 2L4bAk - 2L5b4k - 2L5k) + +12Е*ДL4L4 + 6Е3*д (3L2L4b3 + L2L3 - 3L4b3k - 2L4k) - 6е5*д (3L4bsk + L4k) +
+6G3,^L4L3 - L3k) - 6G*,1 L4k
and all the constants with ★ are obtained from their counterparts through division by C4,0.
4. NUMERICAL RESULTS
Consider a five-layered laminate studied in the previous section with the layers of same thickness h, i.e. h = h, i = 1,2,..., 5. Then, we have from (17) and (24) for eigenfrequencies, and from (19), (20), and (21) for eigenform at k = 1 and k = —1, respectively
(25)
п2
e (i _ - (5c2 + 2)+ ...),
with
ui = Zi + 6Zi (-Z?c2 + c2 + 24) + ...,
U2 = 1 + 2 (-Z? + 4Z2 + 5) + ..., 2
из
U4
1 + - (-c2Z32 + 5c2Z3 - 6c2 + 9) + ...,
(26)
1 - 2Z4 (Z4 - 6) + ...,
И5 = -Z5 + 5 + 6 (c2Z53 - 15c2Z52 - 24Z5 + 74Z5c2 - 120c2 + 120) + ...
and with
02 = - (3 - 2 (3c2 + 2) +...),
(27)
ui
-Zi + 2Zi (Z?c2 - c2 - 40) + ...,
U2 = -1 + 2 (3Z22 - 8Z2 - 35) + ..., 2
из = 2Z3 - 5 + 2 (-2c2Z33 + 15c2Z32 - 37c2Z3 + 78Z3 + 30c2 - 195) + ..., 2
U4 = 1 - 2Z4 (3Z4 - 22) + ..., 2
(28)
u5 = -Z5 + 5 + - (c2- 15c2Zf + 74Z5c2 - 40Z5 - 120c2 + 200) + ... 2
The eigenfrequencies, calculated from one-term and two-term asymptotic formulae, see (25) and (27), versus their exact values, for which detM = 0 in (30), are presented in the Table. Here and below e = 0.01 and c =1.
Comparison of exact and approximate eigenfrequencies
k 02 = -02o 02 = - (O2o + -02i ) exact value 02
1 0.01 0.009883 0.009884
-1 0.03 0.02975 0.02975
Fig. 3 demonstrates a good agreement between two-term expansions (26) and (28) and exact solution (29) for both k = ±1. In this figure the exact eigenforms are normalized by constant U24 (U24 « C40), see Appendix.
u/C4fi] u/C4,Q _ 1
0.8 0.6 0.4 0.2 0
3 4 z/h
0.5
-0.5
-1
1— —H
1 \\
1 \\
1 \\
\\ \\
\ 1 2 3 4 z/h
\
\ //
V/
b
Fig. 3. Asymptotic expansions (26), k = 1 (a) and (28), k = -1 (b) of the eigenforms of a regular five-layered laminate at e = 0 (dashed line) and e = 0.01, c =1 (dotted line) together with exact solution (29) (solid line)
1
5. CONCLUDING REMARKS
The developed methodology seems to be the initial stage in analysing low-frequency vibrations of strongly inhomogeneous multi-layered structures. The next step, following the evaluation of low-frequency thickness resonances, is concerned with the derivation of multi-mode polynomial approximations of the original dispersion relations, similarly to [6]. The final expected outcome should involve two-dimensional equations of motion supplied with appropriate boundary conditions at the edges. The proposed approach may be easily extended to layered anisotropic and pre-stressed structures, previously investigated only for homogeneous and non-contrast configurations, e.g. see [16-18]. Certainly, various generalisations dealing with another scaling for problem parameters are also possible.
Acknowledgements: This work was supported by the grant J2-9224 from the Slovenian Research Agency and also by Program OP20.00362 EVA4GREEN. L. A. Prikazchikova acknowledges support by the Faculty Award by Keele University, UK.
APPENDIX
The solution of equations (5) for a five-component laminate can be written as
ui = Uii sin (OZi) + Ui2 cos (OiZi), i = 1,..., 5, (29)
where Ui1 and Ui2 are constants. On substituting (29) into continuity conditions (9), (10) and boundary conditions (8), we arrive at the eigenvalue problem
M • U = 0, (30)
where U = (U11? U12, U21, U22, U31, U32, U41, U42, U51, U52)T, is an eigenvector and M is a 10 x 10 matrix with non-zero components given by
M12 = 1, M21 = sin(01), M22 = COs(01), M23 = - sin(02&2), M24 = - COs(02&2), M33 = sin(02(b2 + 1)), M34 = cos(02(b2 + 1)), M35 = - sin(03b3), M36 = - cos(03b3), M45 = sin(03(b3 + 1)), M46 = cos(03(b3 + 1)), M47 = - sin(04b4), M48 = - cos(04b4), M57 = sin(04 (b4 + 1)), M58 = cos(04 (b4 + 1)), M59 = - sin(05 b5), M510 = - cos(05 b5),
M69 = sin(05 (b5 + 1)), M610 = cos(05 (b5 + 1)), M71 = 01 cos(01 )e, M72 = -01 sin(01)e, M73 = -¿2 cos(02b2)02, M74 = ¿2 sin(02b2)02, M83 = O2 cos(02(b2 + 1)), M84 = -O2 sin(0(b2 + 1)), M85 = -¿3 cos(0b3)0e, M86 = ¿3 sin(03b3)03e, M95 = 03 cos(03(b3 + 1))e, M96 - 03 sin(03(b3 + 1))e,
M97 = -L3 cos(04 b4 )04, M98 = L4 sin(04 b4 )04, M107 = 04 cos(04 (b4 + 1)), M108 = -04 sin(04 (b4 + 1)), M109 = -cos(05 b5 )05 e, M1010 sin(05 b5 )05 e.
References
1. Kossovich L. Yu. Nestacionarnye zadachi teorii uprugih tonkih obolochek [Nonstationary problems in the theory of elastic thin shells]. Saratov, Saratov Univ. Press, 1986 (in Russian).
2. Kaplunov J. D., Kossovich L. Y., Nolde E. V. Dynamics of thin walled elastic bodies. Academic Press, 1998.
3. Beresin V. L., Kossovich L. Y., Kaplunov J. D. Synthesis of the dispersion curves for a cylindrical shell on the basis of approximate theories. Journal of sound and vibration, 1995, vol. 186, no. 1, pp. 37-53. DOI: https://doi.org/10.1006/jsvi.1995.0432
4. Le K. C. Vibrations of shells and rods. Berlin, Springer, 1999.
5. Berdichevsky V. Variacionnye principy mekhaniki sploshnoj sredy [Variational Principles of Continuum Mechanics]. Moscow, Nauka, Glav. red. fiz.-mat. lit., 1983 (in Russian).
6. Kaplunov J., Prikazchikov D. A., Prikazchikova L. A. Dispersion of elastic waves in a strongly inhomogeneous three-layered plate. International Journal of Solids and Structures, 2017, vol. 113-114, pp. 169-179. DOI: https://doi.org/10.1016/j-.ijsolstr.2017.01.042
7. Kossovich L. Yu., Shevtsova Yu. V. Asymptotic approximations of three-dimensional dynamic equations of elasticity theory in case of two-layered plates. Problems of strenght and plasticity, 2005, vol. 76, pp. 102-111 (in Russian). DOI: https://doi.org/10.32326/1814-9146-2005-67-1-102-110
8. Prikazchikova L., Ece Aydin Y., Erbas B., Kaplunov J. Asymptotic analysis of an anti-plane dynamic problem for a three-layered strongly inhomogeneous laminate. Mathematics and Mechanics of Solids, 2018. DOI: https://doi.org/10.1177/1081286518790804
9. Kaplunov J., Prikazchikov D., Sergushova O. Multi-parametric analysis of the lowest natural frequencies of strongly inhomogeneous elastic rods. Journal of Sound and Vibration, 2016, vol. 366, pp. 264-276. DOI: https://doi.org/10.1016/jjsv.2015.12.008
10. Vinson J. R. The behavior of sandwich structures of isotropic and composite materials. CRC Press, 1999.
11. Ivanov I. V. Analysis, modelling, and optimization of laminated glasses as plane beam. International Journal of Solids and Structures, 2006, vol. 43, no. 22-23, pp. 6887-6907. DOI: https://doi.org/10.1016/j.ijsolstr.2006.02.014
12. Schulze S.-H., Pander M., Naumenko K., Altenbach H. Analysis of laminated glass beams for photovoltaic applications. International Journal of Solids and Structures, 2012, vol. 49, pp. 2027-2036. DOI: https://doi.org/10.1016/j.ijsolstr.2012.03.028
13. Lee P., Chang N. Harmonic waves in elastic sandwich plates. Journal of Elasticity, 1979, vol. 9, pp. 51-69. DOI: https://doi.org/10.1007/BF00040980
14. Kaplunov J. D. Long-wave vibrations of a thinwalled body with fixed faces. The Quarterly Journal of Mechanics and Applied Mathematics, 1995, vol. 48, no. 3, pp. 311-327. DOI: https://doi.org/10.1093/qjmam/48.3.311
15. Kaplunov J. D., Nolde E. V. Long-wave vibrations of a nearly incompressible isotropic plate with fixed faces. The Quarterly Journal of Mechanics and Applied Mathematics, 2002, vol. 55, no. 3, pp. 345-356. DOI: https://doi.org/10.1093/qjmam/55.3.345
16. Kaplunov J. D., Kossovich L. Yu., Rogerson G. A. Direct asymptotic integration of the equations of transversely isotropic elasticity for a plate near cut-off frequencies. Quarterly Journal of Mechanics and Applied Mathematics, 2000, vol. 53, no. 2, pp. 323-341.
17. Nolde E. V., Rogerson G. A. Long wave asymptotic integration of the governing equations for a pre-stressed incompressible elastic layer with fixed faces. Wave Motion, 2002, vol. 36, no. 3, pp. 287-304. DOI: https://doi.org/10.1016/S0165-2125(02)00017-3
18. Rogerson G. A., Sandiford K. J., Prikazchikova L. A. Abnormal long wave dispersion phenomena in a slightly compressible elastic plate with non-classical boundary conditions. International Journal of Non-Linear Mechanics, 2007, vol. 42, no. 2, pp. 298-309. DOI: https://doi.org/10.1016/j.ijnonlinmec.2007.01.005
Cite this article as:
Kaplunov J. D., Prikazchikova L. A. Low-Frequency Vibration Modes of Strongly Inhomogeneous Elastic Laminates. Izv. Saratov Univ. (N. S.), Ser. Math. Mech. Inform., 2018, vol. 18, iss. 4, pp. 447-457. DOI: https://doi.org/10.18500/1816-9791-2018-18-4-447-457
УДК 539.3
ДЛИННОВОЛНОВЫЕ МОДЫ КОЛЕБАНИЙ СИЛЬНО НЕОДНОРОДНЫХ УПРУГИХ СЛОИСТЫХ КОНСТРУКЦИЙ
Ю. Д. Каплунов, Л. А. Приказчикова
Каплунов Юлий Давидович, доктор физико-математических наук, профессор, Школа компьютерных наук и математики, Университет г. Киль, Киль, Стаффордшир, ЭТ5 БЭЭ, Великобритания, j.kaplunov@keele.ac.uk
Приказчикова Людмила Анатольевна, кандидат физико-математических наук, Школа компьютерных наук и математики, Университет г. Киль, Киль, Стаффордшир, ЭТ5 5Б0, Великобритания, l.prikazchikova@keele.ac.uk
В статье изучается динамика тонких многослойных конструкций с контрастными «мягкими» и «жесткими» слоями. Разработана асимптотическая процедура для анализа малых частот среза. Получено полиномиальное частотное уравнение, а также линейные уравнения для собственных форм тол-щинных колебаний. В случае пятислойной пластины с зажатыми лицевыми поверхностями выведены двухчленные асимптотические разложения для собственных частот и форм, которые сопоставляются с точным решением исходной задачи о свободных толщинных колебаниях.
Ключевые слова: асимптотика, контраст, слоистые упругие конструкции, толщинные колебания, малые частоты, срез.
Благодарности. Статья выполнена при финансовой поддержке Словенского исследовательского агентства (проект № Л2-9224) и программы ОР20.00362 ЕУЛ40кЕЕМ. Работа Л. А. Приказчиковой выполнена при финансовой поддержке премии факультета Университета Киля, Великобритания.
Образец для цитирования:
Kaplunov J. D., Prikazchikova L. A. Low-Frequency Vibration Modes of Strongly Inhomogeneous Elastic Laminates [Каплунов Ю. Д., Приказчикова Л. А. Длинноволновые моды колебаний сильно неоднородных упругих слоистых конструкций] // Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2018. Т. 18, вып. 4. С. 447-457. DOI: https://doi.org/10.18500/1816-9791-2018-18-4-447-457