long-wave propagation in multi-layered and multi-component ... · int. j. mech. sci, 131, 171-178....
TRANSCRIPT
![Page 1: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/1.jpg)
Long-wave propagation in multi-layered and multi-component
strongly inhomogeneous waveguides
Julius Kaplunov
(with D. Prikazchikov and L. Prikazchikova)
Keele University, UK
“Wave propagation in complex and microstructured media”
Corsica, August 2019
![Page 2: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/2.jpg)
Outline
• Introduction
• Strongly inhomogeneous rods
• Anti-plane motion of circular layered cylinders
• Layered plates
• Coated half-space
• Periodic structures
![Page 3: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/3.jpg)
1. Introduction
High-contrast layered structures• photovoltaic panels • laminated glass
www.dupont.com www.gscglass.com
![Page 4: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/4.jpg)
Introduction
Sandwich structures
• Foam insulation panels
Rus & Tolley, 2015. Design, fabrication and control of soft robots. Nature, 521(7553), 467.
Soft robots• Classical sandwich plate
Stokes et al. A hybrid combining hard and soft robots. Soft Robotics 1.1 (2014): 70-74
Bio-composites
• (Teeth, bones, etc., contain bothsoft protein/collagen matrix and hard mineral inclusions)
Slesarenko et al., 2017. Understanding the strength of bioinspired soft composites.Int. J. Mech. Sci, 131, 171-178.
![Page 5: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/5.jpg)
2. Low-frequency vibrations of multi-component high-contrast elastic rods
Contrast in• Stiffness
• Density
• Length
Small parameters asymptotic methods
Physical intuition:
Strong components (free ends) - almost rigid body motions
Weak components (fixed b.c.) - almost homogeneous deformations
,j
i
E
E,
j
i
j
i
l
l
J. Kaplunov et al., J. Sound Vib., 445 (2019): 132-147
![Page 6: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/6.jpg)
Toy problem: three-component rod (antisymmetric)
J. Kaplunov et al. J. Sound Vib. 366 (2016): 264-276
![Page 7: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/7.jpg)
![Page 8: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/8.jpg)
![Page 9: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/9.jpg)
.1| ,
;1| ,1
χ|χ
χ|U
Approximate polynomial eigenform
![Page 10: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/10.jpg)
![Page 11: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/11.jpg)
Multi-component high-contrast elastic rods
. , ,~ ,1 2
21
1
2
m
mm
i
jj
i
Ec
l
lL
E
E
,02
2
2
ii
i
i udX
udEquations of motion
Boundary conditions
Problem parameters
.21 ;,1 , 1, , ,1
0
1 ,mnillbbXbc
l
l
xX
i
k
kiiiii
m
ii
i
ii
. ,
11
11
1
11
1
11
iiii
iiii
bXi
ii
i
j
bXi
i
bXibXidX
duL
dX
duuu
,0
101
1
1
nn bXn
n
XdX
du
dX
du
Continuity
Dimensionless scaling
)lyrespective soft,or stiff being component for 1or 1( thi-j
![Page 12: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/12.jpg)
Multi-component high-contrast elastic rods
• Leading order problem
for stiff components
Asymptotic expansions
Global low-frequency regime . ... ,~ 2
1
2
0
22 iiii
, ...10 iii uuu
,02
0
2
i
i
dX
ud.0
1
00
iiii bXi
i
bXi
i
dX
du
dX
du
const 0 ii Cu( almost rigid body motion )
![Page 13: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/13.jpg)
Multi-component high-contrast elastic rods
• Leading order problem for soft components
,02
0
2
i
i
dX
ud
. 1110 iiiiii bXCCCu
( almost homogeneous deformation )
. , 111 ibXiibXi CuCu
iiii
![Page 14: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/14.jpg)
Multi-component high-contrast elastic rods
From solvability of next order for stiff components
.,...,5,3,1 ni
Polynomial equation for frequency!
,1 1
31
2
2
10
C
CL
, 21
2111
2
0
i
ii
i
i
ii
i
i
i
i
iiC
CL
C
CLLL
; 1 21
2
0
n
nn
nnC
CL
• Leading order eigenform – piecewise linear (constant for stiff parts)
• Next order – correction for frequencies, polynomial correction to the eigenform
![Page 15: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/15.jpg)
Example: Five-component rod (free ends)
Bicubic frequency equation
02
102
4
101
6
10 aa
Approximate polynomial eigenform
.2
4or 0 ,
1
2
2
2
111
2
2
10L
aaakkkL
• For k = 0 the solution is exact – rigid body motion!
![Page 16: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/16.jpg)
3. Antiplane motion of concentric circular cylinders
.2,1 , , ,~ ,1 2
21
2
1 mcl
lL
m
mm
i
jj
i
.02
2
2
iii
i
i
i
ii uR
dR
du
dR
udREquations of motion
Boundary conditions
Problem
parameters.,1 , 1, , ,
1
0
1 nillbbRbc
l
l
rR
i
k
kiiiii
m
ii
i
ii
. ,
11
11
1
11
1
11
iiii
iiii
bRi
ii
i
j
bRi
i
bRibRidR
duL
dR
duuu Continuity
)lyrespective soft,or stiff being component for 1or 1( thi-j
. 011
11
nn bRnbR
uu
![Page 17: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/17.jpg)
Antiplane motion of concentric circular cylinders
Summary of the approach
• Global low-frequency perturbation
• Rigid body motions of stiffer components
(at leading order)
• Leading order solution for softer components, involving logarithmic functions
• Solvability of the next order problem
for stiffer componentsPolynomial equation for frequency
![Page 18: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/18.jpg)
Example: Three-layered cylinder
Frequency
Eigenform
![Page 19: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/19.jpg)
Example: Five-layered cylinder
( two roots for k )Frequency
Eigenform
![Page 20: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/20.jpg)
Example: square cylinder with a circular annular inclusion
Frequency
Eigenform
2
2
2
2
1
212
2014
2
bb
L
J. Kaplunov et al. Adv. Struct Mat. 46 (2017): 265-277
![Page 21: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/21.jpg)
4. High-contrast three-layered plates (antisymmetric)
Preliminary remarks
• Rayleigh-Lamb dispersion relation for a single-layered plate
NO CHANCE OF TWO-MODE APPROXIMATIONS!
J. Kaplunov et al. Dynamics of thin walled elastic bodies, Academic Press, 1998
![Page 22: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/22.jpg)
![Page 23: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/23.jpg)
Low-frequency vibrations of high-contrast three-layered plates
Kaplunov et al. Int. J. Solids Struct. 113 (2017): 169-179
![Page 24: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/24.jpg)
![Page 25: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/25.jpg)
NEED OF TWO-MODE MODELS!
No contrast
Effect of contrast
![Page 26: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/26.jpg)
Frequency equation
Frequency
![Page 27: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/27.jpg)
Some three-layered structures satisfying the condition
A) Photovoltaic panels B) Laminated glass
~ 1,~ ,1 h 1~ ,~ ,1 -1/4 h
C) Sandwich structure
2~ ,~ ,1 h
(stiff skin layers and light core layer ) (stiff skin layers and light thin core layer ) (stiff thin skin layers and light core layer )
UNEXPECTEDLY LOW FIRST SHEAR CUT-OFF FREQUENCIES!
![Page 28: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/28.jpg)
![Page 29: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/29.jpg)
Low-frequency dispersion behaviour
A) Photovoltaic panels(stiff skin layers and light core layer )
UNIFORM TWO-MODE APPROXIMATIONS
![Page 30: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/30.jpg)
![Page 31: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/31.jpg)
CLASSICAL KIRCHHOFF-TYPE THEORY IS NOT APPLICABLE!
![Page 32: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/32.jpg)
![Page 33: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/33.jpg)
![Page 34: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/34.jpg)
Low-frequency dispersion behaviour
UNIFORM TWO-MODE APPROXIMATIONS
B) Laminated glass
(stiff skin layers and light light core layer )
J. Kaplunov et al. Proc. Eng. 199 (2017): 1489-1494
1~ ,~ ,1 -1/4 h
![Page 35: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/35.jpg)
Low-frequency dispersion behaviour
COMPOSITE APPROXIMATIONS
C) Sandwich structure
(stiff thin skin layers and light core layer )
![Page 36: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/36.jpg)
NO OVERLAP REGION!
![Page 37: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/37.jpg)
Lightweight structures
Where is transition from uniform approximation to a composite one?
Low thickness shear cut-off frequency
J. Kaplunov et al., to appear
![Page 38: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/38.jpg)
Uniform approximation
3
10 a Composite approximation
1
3
1a
9
1a
3
1a
3
2a
Limiting case
• Two-mode uniform
• Composite
![Page 39: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/39.jpg)
Anti-plane antisymmetric motion
s
①
①
♥ ②r
♦r ②r
Equations of motion
∂σq13∂x1
+∂σq23∂x2
− ρq∂2uq
∂t2= 0, q = c, s,
with
σqi3 = µq∂uq
∂xi, i = 1, 2,
uq are out of plane displacements, σqi3 are shear stresses.
![Page 40: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/40.jpg)
Dispersion relation
Continuity conditions along interfaces x2 = ±hc
σc23 = σs23 and uc = us.
Traction-free boundary conditions
σs23 = 0 at x2 = ±(hc + hs).
Equations of motion
∆uq −1
(cq2)2∂2uq
∂t2= 0, q = c, s.
Dispersion relation
µα1 cosh(α1) cosh(α2h) + α2 sinh(α1) sinh(α2h) = 0,
with
α1 =√
K2 − Ω2, α2 =
√K2 − µ
ρΩ2,
Ω =ωhc
cc2, K = khc, h =
hs
hc, µ =
µcµs, ρ =
ρcρs.
![Page 41: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/41.jpg)
Exact solutions for displacements and stresses
uc = hcsinh(α1ξ2c)
α1, σc13 = iµcK
sinh(α1ξ2c)
α1, σc23 = µc cosh(α1ξ2c),
and
us = hcβ (cosh [α2(hξ2s + 1)]− tanh [α2(h + 1)] sinh [α2(hξ2s + 1)]) ,
σs13 = iµsKβ (cosh [α2(hξ2s + 1)]− tanh [α2(h + 1)] sinh [α2(hξ2s + 1)]) ,
σs23 = µsα2β (sinh [α2(hξ2s + 1)]− tanh [α2(h + 1)] cosh [α2(hξ2s + 1)]) ,
where
β =sinhα1
α1
(coshα2 − sinhα2 tanh[α2(h + 1)]
) .Dimensionless variables
ξ2c =x2
hc, 0 ≤ x2 ≤ hc,
ξ2s =x2 − hc
hs, hc ≤ x2 ≤ hc + hs.
![Page 42: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/42.jpg)
Long-wave low-frequency limit
Polynomial dispersion relation
µ+ γ1K2 + γ2K
4 + γ3K2Ω2 + γ4Ω2 + γ5Ω4 + · · · = 0,
with
γ1 =µ
2
(1 + h2
)+ h,
γ2 =µ
24
(1 + 6h2 + h4
)+
h
6(1 + h2),
γ3 = − µ
12(1 + 3h2)− h
6− µh
12ρ(2 + 3µh)− µh3
12ρ(4 + µh) ,
γ4 = −µ2− µh
ρ
(1 +
µh
2
),
γ5 =µ
24+µh
12ρ(2 + 3µh) +
µ2h3
24ρ2(4 + µh).
![Page 43: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/43.jpg)
Dispersion curves
No contrast
0 1.18 3.63 6.16
0.5K
Ω
µ = 0.232, ρ = 3.0, h = 1.0
Effect of contrast
00.17 3.13 4.60 6.29
0.5K
Ω
µ = 0.014, ρ = 0.03, h = 1.0
No fundamental mode. It appears in case of symmetricmotion.
The lowest cut-off frequency in case of a contrast isΩ = 0.17
Consider two setups of the contrast:A. Photovoltaic panels and B. Sandwich structures
![Page 44: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/44.jpg)
A. Photovoltaic panels. Shortened polynomial dispersion relation
Plate with stiff outer layers and light coreµ 1, h ∼ 1, ρ ∼ µ
γ1 ∼ γ2 ∼ γ3 ∼ γ4 ∼ γ5 ∼ 1.
Shortened dispersion relation
µ
h+ K2 − 1
ρµΩ2 = 0.
Scaled dimensionless frequency and wavenumber
Ω2 = µαΩ2∗ and K2 = µαK2
∗,
where Ω∗ ∼ K∗ ∼ 1 and 0 < α ≤ 1.
α covers the whole long-wave low-frequency band, given byΩ 1, and K 1.
![Page 45: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/45.jpg)
Shortened polynomial dispersion relation
Dispersion relation expressed in Ω∗ and K∗ becomes
Ω2∗ = ρµ
(K2
∗ +µ1−α
h
).
At α < 1 we have Ω∗ ∼√ρµK∗ or ω ∼ cs2k, corresponding to the
short-wave limit for stiffer skin layers.
0 0.17 0.5 1.0 1.5
0.5
1.0
K
Ω
![Page 46: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/46.jpg)
B. Sandwich structure. Shortened polynomial dispersion relation
Plate with stiff outer layers and light coreµ 1, h ∼ µ, ρ ∼ µ2
γ1 ∼ γ2 ∼ µ and γ3 ∼ γ4 ∼ γ5 ∼ 1.
Approximate dispersion relation
µ+ µ
(1
2+ hµ
)K2 − hµ
6ρµK2Ω2 −
(µ
2+
hµρµ
)Ω2 +
hµ6ρµ
Ω4 = 0.
Normalized wavenumber and frequency
K2 = µK2∗ and Ω2 = µΩ2
∗,
we obtain
1 +µ
(1
2+ hµ
)K2
∗−µhµ6ρµ
K2∗Ω2
∗−(µ
2+
hµρµ
)Ω2∗ +µ
hµ6ρµ
Ω4∗ = 0.
![Page 47: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/47.jpg)
Shortened polynomial dispersion relation
Adapt a near cut-off asymptotic expansion in the form
Ω2∗ = Ω2
0 + µΩ21 + · · ·
where
Ω20 =
ρµhµ
and Ω21 =
ρµhµ
(1
3+ hµ
)K2
∗ −1
3
ρ2µh2µ
,
leading to the optimal shortened dispersion relation(hµ +
1
3
)K2 − 1
µ
hµρµ
Ω2 +
(1− µρµ
3hµ
)= 0.
Valid only over a narrow vicinity of the cut-off frequency!
![Page 48: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/48.jpg)
Numerical illustration
0.1 0.14 0.2
0.5
1.0
K
Ω
µ = 0.014, ρ = 0.03, and h = 1.0
![Page 49: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/49.jpg)
Asymptotic formulae for displacements and stresses (setup A)
Leading order displacements and stresses
uc = hcξ2c,
σc13 = iµc√µK∗ξ2c,
σc23 = µc,
and
us = hc,
σs13 = iµs√µK∗,
σs23 = µch
(K2
∗ −Ω2∗
ρµ
)(ξ2s − 1) .
We obtainuq
hc∼ σq23
µc∼ σq13µq√µ, q = c, s.
![Page 50: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/50.jpg)
Normalised displacement and stress σ23 (setup A)
ξ2 = ξ2c, u =uc
hc, and σ23 =
σc23µc
, (0 < ξ2 ≤ 1)
or ξ2 = 1 + ξ2s, u =us
hc, and σ23 =
σs23µc
, (1 < ξ2 ≤ 2)
0 0.5 1.0 1.5
0.5
1.0
1.5
u
ξ2
0 0.5 1.0 1.5
0.5
1.0
1.5
σ23
ξ2
![Page 51: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/51.jpg)
Model construction (setup A)
Scaled longitudinal coordinate and time
x1 =hc√µξ1 and t =
hc
c2c√µτ,
Normalised displacement and stresses
uq = hcvq, σq13 = µq
õSq
13, σq23 = µcSq23, q = c, s.
with all dimensionless quantities assumed to be of order unity.
Core layer
µ∂Sc
13
∂ξ1+∂Sc
23
∂ξ2c− µ∂
2vc
∂τ2= 0,
Sc13 =
∂vc
∂ξ1, Sc
23 =∂vc
∂ξ2c.
Skin layer
∂Ss13
∂ξ1+
1
h
∂Ss23
∂ξ2s− 1
ρµ
∂2vs
∂τ2= 0,
Ss13 =
∂vs
∂ξ1, µhSs
23 =∂vs
∂ξ2s.
![Page 52: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/52.jpg)
Derivation of a shortened equation (setup A)
Continuity and boundary conditions
vc∣∣ξ2c=1
=vs∣∣ξ2s=0
,
Sc23
∣∣ξ2c=1
=Ss23
∣∣ξ2s=0
,
and
Ss23
∣∣ξ2s=1
= 0.
Expand displacements and stresses into asymptotic series as
vq =vq0 + µvq
1 + · · · ,Sqj3 =Sq
j3,0 + µSqj3,1 + · · · , q = c, s and j = 1, 2.
![Page 53: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/53.jpg)
Leading order problem
Sc13,0 =
∂vc0
∂ξ1,
∂Sc23,0
∂ξ2c= 0, Sc
23,0 =∂vc
0
∂ξ2c,
and
∂Ss13,0
∂ξ1+
1
h
∂Ss23,0
∂ξ2s− 1
ρµ
∂2vs0
∂τ2= 0,
Ss13,0 =
∂vs0
∂ξ1,
∂vs0
∂ξ2s= 0,
with
vc0
∣∣∣ξ2c=1
=vs0
∣∣∣ξ2s=0
,
Sc23,0
∣∣∣ξ2c=1
=Ss23,0
∣∣∣ξ2s=0
,
and
Ss23
∣∣ξ2s=1
= 0.
![Page 54: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/54.jpg)
Leading order solution
vs0 = w(ξ1, τ).
The rest of the quantities are expressed in terms of w as
Sc13,0 =ξ2c
∂w
∂ξ1, Sc
23,0 = w, vc0 = ξ2cw,
Ss13,0 =
∂w
∂ξ1, Ss
23,0 = w(1− ξ2s),
with w satisfying the 1D equation
∂2w
∂ξ21− 1
ρµ
∂2w
∂τ2− 1
hw = 0,
which may be presented in the original variables as
∂2us
∂x21
− ρsµs
∂2us
∂t2− µcµshchs
us = 0,
where us(x1, t) ≈ w(x1, t).
![Page 55: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/55.jpg)
Justification of the model
Insert ansatz us = exp i(kx1 − ωt) into the last equation. As aresult, we have the dispersion relation
k2 − ρsµsω2 +
µcµshchs
= 0.
Coincides with the shortened dispersion relation for setup A!
![Page 56: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/56.jpg)
Anti-plane shear of three-layered asymmetric plates
-h3
h2
2
-h2
2
h1
μ1
μ2
μ1
ρ1
ρ2
ρ1
x1
x2
![Page 57: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/57.jpg)
More sophisticated dispersion relation
µα1α2 tanh(hα1) + µ2α22 tanh(α2)+
µα1α2 tanh(h∗α1) + α12 tanh(h∗α1) tanh(α2) tanh(hα1) = 0,
where
α1 =
√K2 − µ
ρΩ2, α2 =
√K2 − Ω2,
with
Ω =ωh2
c(2)2
, K = kh2,
and
h =h1
h2, h∗ =
h3
h2, µ =
µ2µ1, ρ =
ρ2ρ1, c
(i)2 =
√µiρi, i = 1, 2
![Page 58: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/58.jpg)
Effect of contrast
No contrast
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4
0.0
0.2
0.4
0.6
0.8
1.0
Ω
k
Contrast parameters
0.0 0.2 0.4 0.6 0.8 1.0
0.0
0.2
0.4
0.6
0.8
1.0
Ω
k
Two modes in case of high contrast for a scalar problem!
![Page 59: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/59.jpg)
Cut-off frequencies
Frequency equation
√µρ
(tan
(h
õ
ρΩ
)+ tan
(h∗√µ
ρΩ
))+ µρ tan (Ω)− tan
(h
õ
ρΩ
)tan (Ω) tan
(h∗√µ
ρΩ
)= 0.
Lowest cut-off
Ω ≈
√µρ(h + h∗ + ρ)
hh∗µ
![Page 60: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/60.jpg)
A. Photovoltaic panels. Two-mode approximation
Shortened polynomial dispersion relation for two modes
G1K2+G2Ω2+G3K
4+G4K2Ω2+G5Ω4+G6K
4Ω2++G7K2Ω4 = 0
Out[ ]=
0.0 0.2 0.4 0.6 0.8 1.0
0.0
0.2
0.4
0.6
0.8
1.0
ω
k Appr Dispersion Relation
Exact Dispersion Relation
![Page 61: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/61.jpg)
Coated half-space
H.-H. Dai, J. Kaplunov, D.A. Prikazchikov. A long-wave model for the surface elastic wave in a coated half-space. Proc. Roy. Soc. A,
466, 3097–3116 (2010)
Singularly perturbed hyperbolic equation on the interface x3 = 0
ϕ,11 − c−2R ϕ,tt − bhϕ,111 = ARP,
isotropic coating b =µ0
2µB(1− β2
R)
(c2Rc220
(αR + βR)− 4βR
(1− c220
c210
))orthorhombic coating b =
c0662µB
(1− β2R)
(c2Rc260
(αR + βR)− 4βRc2c0c260
),
with c260 = c066/ρ0, c2c0 =(c011c
022 − (c012)2
)/ρ0.
![Page 62: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/62.jpg)
Hyperbolic-elliptic model for the Rayleigh wave on an isotropichalf-space
• Formulation of the problem: (−∞ < x1 <∞, 0 ≤ x3 <∞)
Equations of motion c21∆ϕ− ϕ,tt = 0, c
22∆ψ − ψ,tt = 0.
Boundary conditions σ33
∣∣x3=0
= P (x1, t), σ31
∣∣x3=0
= Q(x1, t).
• Asymptotic model:
Elliptic equations ϕ,33 + α2Rϕ,11 = 0, ψ,33 + β
2Rψ,11 = 0,
where αR =
√1 −
c2Rc21, βR =
√1 −
c2Rc22.
Potentials are related ϕ(x1 − cRt, αRx3) = γψ(x1 − cRt, αRx3).
Hyperbolic equations on the surface x3 = 0
ϕ,11 − c−2R ϕ,tt = ARP, ψ,11 − c
−2R ψ,tt = −ARQ.
where AR =1 + β2
R
2µB, B =
αR
βR
(1 − β2R) +
βR
αR
(1 − α2R) − 1 + β
4R.
Kaplunov, J., Prikazchikov, D.A.: Asymptotic theory for Rayleigh and Rayleigh-type waves. Advances in Applied Mechanics
50, 1–106 (2017)
![Page 63: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/63.jpg)
Clamped surface
Homogeneous half-space
⇒ No surface wave
Inhomogeneous half-space
Rigorous mathematical analysis proves a possibility of localized waves
K.D. Cherednichenko, S. Cooper. On the existence of high-frequency boundary resonances in layered elastic media. Proc. R. Soc.
A, 471(2178), 20140878 (2015)
![Page 64: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/64.jpg)
Coated half-space with fixed surface. Problem statement
Equations of motion
σ±i1,1 + σ±
i3,3 − ρ±u±
i,tt = 0, i = 1, 3.
Constitutive relations
σ±ij = λ±δij(u
±1,1 + u±
2,2 + u±3,3) + µ±(u±
i,j + u±j,i), j = 1, 3.
Boundary conditions at x3 = −h
u−i = 0.
– Novelty!
Continuity conditions at x3 = 0
u−i = u+
i , σ−i3 = σ+
i3.
![Page 65: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/65.jpg)
Coated half-space with fixed surface. Exact solution
Elastic wave potentials
u±1 =
∂ϕ±
∂x1− ∂ψ±
∂x3, u±
3 =∂ϕ±
∂x3+∂ψ±
∂x1.
Wave equations
∆ϕ± − 1
(c±1 )2∂2ϕ±
∂t2= 0, ∆ψ± − 1
(c±2 )2∂2ψ±
∂t2= 0,
where c±1 =
√λ± + 2µ±
ρ±, c±2 =
√µ±
ρ±and ∆ =
∂2
∂x1+
∂2
∂x3.
The sought for wave potentials
ϕ− = [A1 cos(α−kx3) +A2 sin(α−kx3)]eik(x1−ct),
ψ− = [A3 cos(β−kx3) +A4 sin(β−kx3)]eik(x1−ct),
ϕ+ = A5eik(x1−ct)−α+kx3 , ψ+ = A6eik(x1−ct)−β
+kx3 ,
where
α− =
√c2
(c−1 )2− 1, α+ =
√1− c2
(c+1 )2, β− =
√c2
(c−2 )2− 1, β+ =
√1− c2
(c+2 )2.
![Page 66: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/66.jpg)
Analysis of full dispersion relation
Dispersion relation detA = 0, where A is a 6× 6 matrix.
Material parameter (relative stiffness) µ =µ−
µ+.
Dimensionless variables ζ =c
c+2, Ω+ =
ωh
c+2, K = kh.
Numerical values ρ− = 150, ρ+ = 250, ν− = 0.3, ν+ = 0.25.
µ = 0.1 µ = 0.3 µ = 0.5
0 2 4 6 8 100
2
4
6
8
10
Ω+
K
0 2 4 6 8 100
2
4
6
8
10
Ω+
K
0 2 4 6 8 100
2
4
6
8
10
Ω+K
![Page 67: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/67.jpg)
Love waves
Dispersion relation tan(Kβ−) +µβ−
β+= 0.
μ=0.10 2 4 6 8 10
0
2
4
6
8
10
Ω+
K
μ=0.30 2 4 6 8 10
0
2
4
6
8
10
Ω+
Kμ=0.5
0 2 4 6 8 100
2
4
6
8
10
Ω+
K
Existence of surface waves
Existence triangle c2
- < c2+
or
μ-
μ+< ρ-
ρ+
c = c2+
c = c2-
0 2 4 6 8 100
2
4
6
8
10
Ω+
K
Initial pointsc = c+2 ⇒ β+ → 0⇒ tan(Kβ−)→∞
Ω+in = Kin =
πn
2
√ρ
µ− 1
, n = 1, 3, ...
Approximations for K 1c→ c−2 ⇒ β− → 0⇒ tan(Kβ−)→ 0
Kβ− = πm, m = 1, 2, ...
![Page 68: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/68.jpg)
Analysis of full dispersion relation
Initial points
R(Ω+) = 0,
where R(Ω+) = detA at c = c+2 or K = Ω+.
μ=0.31
2
3
0 2 4 6 8 100
2
4
6
8
10
Ω+
K
μ=0.5
1
2
0 2 4 6 8 100
2
4
6
8
10
Ω+
K
μ=0.3
1 2 3
0 2 4 6 8 10
-0.4
-0.2
0.0
0.2
0.4
0.6
0.8
1.0
Ω+
R
μ=0.5
1 2
0 2 4 6 8 10
-0.4
-0.2
0.0
0.2
0.4
0.6
0.8
1.0
Ω+
R
Approximations for K 1
Similarly to the Love waves
Ω+ =c−2c+2
√π2m2 +K2
m = 1, 2, ...
12.0 12.5 13.0 13.5 14.0 14.530
31
32
33
34
Ω+
K
![Page 69: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/69.jpg)
Analysis of full dispersion relation
High contrast (soft layer, stiff half-space)µ = 0.001
Rayleigh wave
Thickness shear resonances
Thickness stretch resonance
0.0 0.1 0.2 0.3 0.4 0.5
0.0
0.1
0.2
0.3
0.4
0.5
Ω+
K
⇓Inspiration for using the asymptotic model for Rayleigh-type waves
The thickness resonances areeigenvalues of the problems
u−3,33 +
(Ω−
κ−
)2
u−3 = 0,
u−1,33 + (Ω−)2u−
1 = 0
with u−i = 0, i = 1, 3 at the faces.
u−i = 0
u−i = 0
stretchshear
![Page 70: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/70.jpg)
Thin coating
Effective boundary⇒ conditions for
stresses at x3 = 0
Dirichlet boundary conditions
u−i = 0, x3 = −h, u−
i = vi, x3 = 0,
where vi are prescribed values.
Small parameter (long-wave approximation)
ε =h
l 1.
Dimensionless variables
ξ1 =x1l, ξ3 =
x3 + h
h, Ω− =
ωh
c−2
Therefore Ω− = Ω+ c+2
c−2= Ω+
√ρ
µwith ρ =
ρ−
ρ+.
![Page 71: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/71.jpg)
Thin coating. Asymptotic procedure
Method of direct asymptotic integration of 3D equations in linear elasticity
Goldenveizer, A.L., Kaplunov, J.D., Nolde, E.V.: On Timoshenko-Reissner type theories of plates and shells. Int. J. Solids Struct.,
30, 675-694 (1993)
Kaplunov, J., Prikazchikov, D., Sultanova, L.: Justification and refinement of Winkler–Fuss hypothesis. Z. Angew. Math. Phys.
69(3), 80 (2018)
Scaling of the displacements and stresses
u−i = hu∗−
i , σ−ij = µ−σ∗−
ij , Ω− ∼ 1, vi = hv∗i .
Note that Ω− ∼ 1 is associated with high-frequency localized wave.
Governing equations
εσ∗−i1,1 + σ∗−
i3,3 + (Ω−)2u∗−i = 0, σ∗−
11 = ε(κ−)2u∗−1,1 + ((κ−)2 − 2)u∗−
3,3,
σ∗−33 = ε((κ−)2 − 2)u∗−
1,1 + (κ−)2u∗−3,3, σ∗−
13 = u∗−1,3 + εu∗−
3,1,
where κ± =c±1c±2
=
√2− 2ν±
1− 2ν±.
Boundary conditions u∗−i = 0, ξ3 = 0, u∗−
i = v∗i , ξ3 = 1.
![Page 72: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/72.jpg)
Thin coating. Leading order
Asymptotic series(u∗−i
σ∗−ij
)=
(u−(0)i
σ−(0)ij
)+ ε
(u−(1)i
σ−(1)ij
)+ ...
Leading order equations
σ−(0)i3,3 + (Ω−)2u
−(0)i = 0, σ
−(0)11 = ((κ−)2 − 2)u
−(0)3,3 ,
σ−(0)33 = (κ−)2u
−(0)3,3 , σ
−(0)13 = u
−(0)1,3 .
Leading order boundary conditions
u−(0)i = 0, ξ3 = 0, u
−(0)i = v∗i , ξ3 = 1.
Leading order stresses
σ−(0)33 =
κ−Ω−v∗3
sin
(Ω−
κ−
) cos
(Ω−
κ− ξ3
), σ
−(0)13 =
Ω−v∗1sin(Ω−)
cos(Ω−ξ3).
Sinusoidal, not polynomial thickness variation!
![Page 73: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/73.jpg)
Thin coating. Asymptotic effective boundary conditions
Asymptotic dynamic effective boundary conditions at the interface x3 = 0
σ+33 = µ−κ−Ω+
√ρ
µ
v3h
cot
(Ω+
κ−
√ρ
µ
),
σ+13 = µ−Ω+
√ρ
µ
v1h
cot
(Ω+
√ρ
µ
).
Non-traditional effective boundary conditions, corresponding tohigh-frequency long-wave phenomena failing at thickness resonances(
sin
(Ω+
κ−
√ρ
µ
)= sin
(Ω+
√ρ
µ
)= 0
).
Kaplunov, J.D., Kossovitch, L.Yu., Nolde, E.V.: Dynamics of thin walled elastic bodies. Academic Press (1998)
Kaplunov, J., Krynkin, A.: Resonance vibrations of an elastic interfacial layer. J. Sound Vib. 294(4-5), 663–677 (2006)
![Page 74: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/74.jpg)
Model for Rayleigh-type waves on a stiff half-space
Kaplunov, J., Prikazchikov, D.A.: Asymptotic theory for Rayleigh and Rayleigh-type waves. Advances in Applied Mechanics 50,
1–106 (2017)
Displacements at the surface x3 = 0
v1 =1− β2
R
2ϕ,1, v3 = −1− β2
R
1 + β2R
ϕ,3.
Hyperbolic equation for the volume wave potential at the interface
ϕ,11 −1
c2Rϕ,tt = A(P + ϑQ) = A ϕ,1.
Here P (x1, t) = σ+33 and Q(x1, t) = σ+
13, and the bar may be interpreted in asense of the Hilbert transform.
![Page 75: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/75.jpg)
Effect of contrast (soft coating)
Considering at the boundary x3 = 0 the potential ϕ = Beik(x1−ct) with k > 0
ζ2
ζ2R−√µργζ − 1 = 0,
where ζ =c
c+2, ζR =
cR
c+2, and
γ = KR
[cot
(Ω+
√ρ
µ
)+ κ− cot
(Ω+
κ−
√ρ
µ
)],
with Ω+ =ωh
c+2.
Asymptotic series ζ = ζ(0) +√µζ(1) + ...
Small parameter for soft coating(µ =
µ−
µ+ 1
)At leading order ζ(0) = ζR.
Correction2ζ(0)ζ(1)
ζ2R−√ργζ(0) = 0 ⇒ ζ(1) =
1
2ζ2R√ργ.
![Page 76: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/76.jpg)
Effect of high contrast (soft coating)
Non-uniform approximation for dimensionless velocity
ζ = ζR +√µ
(√ρ
2γ ζ2R
)+ ...
Note, γ is the source of non-uniformity, since γ →∞ at
Ω+ =
õ
ρπn and Ω+ = κ−
õ
ρπn, n = 0, 1, 2, ...
![Page 77: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/77.jpg)
Numerical results
Numerical values ρ− = 150, ρ+ = 250, ν− = 0.3, ν+ = 0.25.
µ = 0.01 µ = 0.001
![Page 78: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/78.jpg)
6. Periodic structures
R. Craster, L. Joseph, J. Kaplunov, Wave Motion, 2014, 51, 581-588
• Similarity between asymptotic procedures for thin and periodic structures
• Knowledge transfer from theory of plates and shells to homogenisation
• High-frequency homogenisation
![Page 79: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/79.jpg)
Industrial motivation
Functionally graded microstructures
the picture is taken from NASA webpage
![Page 80: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/80.jpg)
Industrial motivation
Civil structures
the picture is taken from http://www.wikipedia.org
![Page 81: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/81.jpg)
Industrial motivation
Photonic crystals
the picture is taken from the review of P Russell, Science 2003
![Page 82: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/82.jpg)
Dynamic homogenization
Two toy problems
The goal is to demonstrate the similarity of the homogenization
procedures for 2D thin functionally graded structures and 1D periodic
structures, see R.V. Craster, L.M. Joseph & J. Kaplunov in Wave Motion
2014.
(A) SH waves in a functionally graded layer (2D problem)
x
h
-h
y
0
(B) Longitudinal waves in a periodic rod (1D problem)
xh-h 0
![Page 83: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/83.jpg)
Dynamic homogenization
Problem A
∂2u
∂x2+
∂2u
∂y2+
ω2
c2(y)u = 0
where u = u(x, y)
Problem B
d2u
dx2+
ω2
c2(x)u = 0
where u = u(x)
traction free faces
∂u/∂y|y=±h = 0
periodicity
c(x) = c(x+ 2h)
Small parameter
ǫ = h/L ≪ 1 (L is typical wavelength along x-axis)
Scaling
X = x/L, ξ = α/h, where
α = y α = x
![Page 84: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/84.jpg)
Dynamic homogenization
Dimensionless equations in u(X, ξ)
uξξ + ǫ2uXX + λ2
C2(ξ)u = 0 uξξ + 2ǫuXξ
︸ ︷︷ ︸
the only difference
+ǫ2uXX + λ2
C2(ξ)u = 0
with λ = ωhc0
and C(ξ) = c(ξ)c0
Classical low frequency limit (λ ∼ ǫ)
u(X, ξ) = u0(X, ξ) + ǫu1(X, ξ) + ǫ2u2(X, ξ) + . . .
and
λ2 = ǫ2(λ20 + ǫλ2
1 + ǫ2λ22 + . . .
)
with Neumann boundary conditions
uiξ|ξ=±1 = 0
with periodicity conditions
ui(X, 1) = ui(X,−1),uiξ(X, 1) = uiξ(X,−1)
![Page 85: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/85.jpg)
Dynamic homogenization
At leading order we get over a microscale
u0ξξ = 0resulting in uniform static variation along thickness or cell
u0(X, ξ) = v0(X).
ξ
u0(X, ξ)
0
v0(X )
1-1
Proceeding to higher orders
u1(X, ξ) = 0, λ1 = 0 and u2ξξ = −v0XX −λ22
C2(ξ)v0
Finally, we arrive at 1D homogenized equation
d2v0dx2
+ω2
〈c〉2v0 = 0, with 〈c〉 =
[1
2h
∫ h
−hc−2(z)dz
]−1/2
![Page 86: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/86.jpg)
Dynamic homogenization
Non-classical high frequency limit (λ ∼ 1)
The so-called high frequency long wave theory for thin elastic structures established some time ago
(e.g. see J.D.Kaplunov, L.Yu.Kossovich & E.V.Nolde, Dynamics of Thin Walled Elastic Bodies,
Academic Press, N.-Y. 1998) inspired a more recently developed high frequency homogenization
procedure (see R.V.Craster, J.Kaplunov & A.V.Pichugin in Proc R Soc A 2010, J.Kaplunov & A.Nobili
in Math Meth Appl Sci 2017, and D.J.Colquitt, V. Danishevsky & J.Kaplunov in Math Mech Solids
2018)
At leading order u0(X, ξ) = v0(X)U0(ξ) and U0ξξ +λ2
0
C2(ξ)U0 = 0
with Neumann boundary conditions
U0ξ|ξ=±1 = 0
with periodicity conditions
U0(X, 1) = U0(X,−1),U0ξ(X, 1) = U0ξ(X,−1)
or antiperiodicity conditions (leading to
periodicity with a double period)
U0(X, 1) = −U0(X,−1),U0ξ(X, 1) = −U0ξ(X,−1)
![Page 87: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/87.jpg)
Dynamic homogenization
Eigenvalues λ0 correspond to
thickness resonances cell resonances
ξ
U0(ξ)
0 1-1
The sought for 1D homogenized equation is
h2Tv0xx + (λ2 − λ20)v0 = 0
T =
∫ h
−hU2
0 (z)C−2(z)dz
∫ h
−hU2
0 (z)dz
T takes slightly more complicated form
(see R.V.Craster, J.Kaplunov & A.V.Pichugin in
Proc R Soc A 2010)
![Page 88: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/88.jpg)
Dynamic homogenization
Floquet-Bloch waves
[u(−1)uξ(−1)
]
= exp(i2κε)
[u(1)uξ(1)
]
where κ - Bloch parameter.
Bloch spectra λ(κ) near edges of stop bands
κ = 0 κ =π
2ε
almost periodic solutions almost anti-periodic solutions
![Page 89: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/89.jpg)
Dynamic homogenization
a) piecewise uniform sound speed (constant coefficients)
ξ
C(ξ)
1
r
−1 0 1
b) Mathieu’s equation (variable coefficients)C−2(ξ) = α− 2θ cos ξ
![Page 90: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/90.jpg)
Dynamic homogenization
Piecewise uniform rod (r = 1/3)
0 1 2 30
2
4
6
8
10
12
14
16
18
20(a)
2ε κ
λ
0 1 2 3 4 5
100
102
104
r
-T0(1
) ,T0(2
)
(b)
0 1 2 3 4 5
100
102
104
r-T
0(3),T
0(4)
(c)
0 1 2 3 4 5
102
104
r
-T0(5
),T
0(6)
(d)
![Page 91: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/91.jpg)
Dynamic homogenization
Mathieu’s equation (α = 1, θ = 1/2)
0 0.5 10
1
2
3
4
5
6
7
8
εκ
λ
(a)
0 1 2 3-1
-0.5
0
0.5
1
ξ
(b)
0 1 2 3
-0.5
0
0.5
1
ξ
(c)
U0(ξ,λ
0
(2))
U0(ξ,λ
0
(3))
W1(ξ,λ
0
(2))
W1(ξ,λ
0
(3))
![Page 92: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/92.jpg)
Dynamic homogenization
High frequency homogenization in 2D(see R.V.Craster, J.Kaplunov & A.V.Pichugin in Proc R Soc A 2010)
L
2l
∇x · [a(x)∇xu(x)] + ω2ρ(x)u(x) = 0
with double periodic a(x) and ρ(x), where x = (x1, x2)
Small parameter ǫ = l/L ≪ 1
Scaling X =x
L, ξ =
x
l
![Page 93: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/93.jpg)
Dynamic homogenization
Asymptotic series
u(X, ξ)=u0(X, ξ) + εu1(X, ξ) + ε2u2(X, ξ) + . . .and
λ2 = λ2
0+ ελ2
1+ ε2λ2
2+ . . . , where λ =
ωl
c0
Double periodicity - antiperiodicity conditions
ui (X;−1, ξ2) = ±ui (X; 1, ξ2)uiξ1 (X;−1, ξ2) = ±uiξ1 (X; 1, ξ2)ui (X; ξ1,−1) = ±ui (X; ξ1, 1)uiξ2 (X; ξ1,−1) = ±uiξ2 (X; ξ1, 1)
At leading order
u0 (X, ξ) = v0 (X)U0 (ξ)and
∇ξ · [a(ξ)∇ξU0] + λ2
0c20ρ(ξ)U0 = 0 (*)
with periodicity - antiperiodicity conditions on cell contour
![Page 94: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/94.jpg)
Dynamic homogenization
The final macroscale equation becomes
l2Tij
∂2v0∂xi∂xj
+(λ2 − λ2
0
)v0 = 0 (i, j = 1, 2) (**)
with Tij expressed through double integrals over the domain −1 ≤ ξ1, ξ2 ≤ 1 containing
double periodic eigenfunction U0 (ξ) and a pair of single periodic functions Vi (ξ),calculated from non-homogeneous boundary value problems for the equation (*).
Remarks.
(i) The equation (**) is valid near edges of stop bands.
(ii) The type of the equation (**) depends on problem parameters.
(iii) Simple explicit expressions for the coefficients Tij are available only in the case of thecheckerboard structures with piece-wise parameters governed by (see R.V.Craster,J.Kaplunov, E.Nolde & S.Guenneau, JOSA 2011)
∂2u
∂x21
+∂2u
∂x22
+ω2
c2[1 + g(x1) + g(x2)]u = 0
where gi(xi) = 0 for − 1 ≤ xi < 0; gi(xi) = r2 for 0 ≤ xi < 1
![Page 95: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/95.jpg)
Dynamic homogenization
Applications of high frequency homogenization theory
for checkerboards in optics
Defect modes
![Page 96: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/96.jpg)
Dynamic homogenization
Ultra-refraction
![Page 97: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/97.jpg)
Dynamic homogenization
All-angle negative refraction
For further detail see R.V.Craster, J.Kaplunov, E.Nolde & S.Guenneau, JOSA 2011
![Page 98: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/98.jpg)
Dynamic homogenization
High frequency homogenization for lattice structures
e.g. see R.V.Craster, J.Kaplunov & J.Postnova, QJMAM 2010 for spring mass structures
and E.Nolde, R.V.Craster & J.Kaplunov, JMPS 2011 for frame and truss structures
(2n,2m) (2n+1,2m)
(2n,2m+1) (2n+1,2m+1)
The square lattice
(2n,2m) (2n+1,2m)
(2n,2m+1) (2n+1,2m+1)
Two-scale approach
u = u (X, ξ)
continuous discrete
(ξ = (0, 0), (0, 1), (1, 0), (1, 1))
By applying Taylor series in X and periodicity - antiperiodicity conditions we arrive at a
matrix-differential problem [4× 4]. It is[
(A0 − λ2M)︸ ︷︷ ︸
+εA1(∂i, λ) + ε2A2(∂i∂j , λ) + ...
]
u(X, ξ) = 0
linear algebra
with ε = 1/N ≪ 1, ∂i = ∂/∂Xi, M = diag (M1,M1,M2,M2)
![Page 99: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/99.jpg)
Dynamic homogenization
Variety of homogenized models
e.g. for periodicity in both directions when u (X, ξ) = v0 (X) [0, 0,−1, 1]T
The result is
l4
4 (M2 −M1)
(∇4
xv0 − 4∂2x1∂2x2v0)+
(
λ2 −4
M2
)
v0 = 0
which is particularly handy for analyzing localized phenomena.
-1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 3.50
1
2
3
4
5
6
Dispersion curves: M1=2, M
2=1
Length along BAC
λ2
B A C B
1
1′
2
3
3′
2′
4
4′
5
5′
-1 0 1
-1
0
1
A
BC
κ1
κ 2
![Page 100: Long-wave propagation in multi-layered and multi-component ... · Int. J. Mech. Sci, 131, 171-178. 2. Low-frequency vibrations of multi-component high-contrast elastic rods ... (stiff](https://reader036.vdocument.in/reader036/viewer/2022071012/5fca1abdb35a447f346265a2/html5/thumbnails/100.jpg)
Concluding remarks
Multi-component high-contrast waveguides require specialised theory
High contrast may lead to unexpectedly low natural frequencies
Stronger components subject to Neumann conditions, perform almost rigid body motions
Two-mode theories for long-wave low-frequency motion of layered plates
(asymptotically uniform or composite)
Deep parallels between long-wave dispersion of thin and periodic structures