# Axisymmetric Longitudinal-Bending Waves in a Cylindrical Shell Interacting with a Nonlinear Elastic Medium

Authors: Alexander I. Zemlyanukhin; Andrey V. Bochkarev; Lev I. Mogilevich; Ekaterina G. Tindova
Journal: Modelling and Simulation in Engineering
 (2016)

---

## Abstract

A nonlinear differential equation is derived which describes the propagation of axisymmetric stationary longitudinal-bending waves in infinite cylindrical shell of Timoshenko type, interacting with the external nonlinear elastic medium. A modified perturbation method based on the use of diagonal Pade approximants was applied to build exact solitary-wave solutions of the derived equation in the form of traveling front and the traveling pulse. Numerical solutions of the equation, obtained by means of finite difference method, are in good agreement with the corresponding exact analytical ones.


---

## Body

## 1. Introduction

The shell theory is a well-developed section of deformable solids mechanics and has numerous technical applications. The interest in nonclassical variants of shell theory has grown recently due to the need of analysing nanoobjects elastic properties, carbon nanotubes (CNTs), in particular [1–4].Nonlinearity and dispersion are the main determinants of wave processes in elastic thin-walled structures. Longitudinal waves in rectilinear rods and plates interact with bending waves only in nonlinear approximation. Thus, longitudinal waves parametrically affect bending waves, which serve, in turn, as a nonlinear source for longitudinal waves. As for the shells the longitudinal and normal components of displacement are already connected in the linear approximation [5]. For this reason, the dynamic processes in shells are significantly more complex than similar processes in rods and plates.Most of the tasks of wave propagation in elastic shells are considered in the linear formulation [6]. Analytical study of coupled equations system of shell element motion, written in displacements, is task of some difficulty; therefore researchers need to simplify it. Rather often, the longitudinal displacements’ inertia is neglected; that is, the shell middle surface is considered to be unstretchable. For the axisymmetric case this allows us to integrate the first equation of the system and reduce the problem to solving a single equation for the normal displacement. This approach is traditional for one-dimensional deformable systems and is often referred to as the Kirchhoff hypothesis [7]. However, the wave process in the shells cannot be adequately described with the use of such simplifications.Equations of shell theory contain several small dimensionless parameters, the ratio of the thickness to the radius of curvature, in particular. Therefore, approximate solutions to these equations can naturally be obtained by means of the perturbation method [8], based on the small-parameter expansion of an unknown solution. The perturbation method is widely used to solve various problems of statics and dynamics of deformable structures. For example, in [9], the authors investigate the influence of small first-order and second-order perturbations of structural parameters such as stiffness and mass on the perturbations of structural vibration eigenvalues and eigenvectors. Traditionally, only the first few terms of the perturbation series are retained to find the approximate solution of the problem. However, the requirement of geometricity for this series allows us to find the sum of the series, which gives the problem exact solution [10].The article aims at deriving and analysing the equation modeling the propagation of axisymmetric longitudinal-bending waves in infinite cylindrical shell, interacting with a external nonlinear elastic medium [11, 12]. The assumption of the infinity of the shell is valid if its edges are damped in such an optimal way that incident waves are not reflected [13]. Then, the influence of the boundary conditions can be ignored and vibrations propagation through the shell can be considered as traveling elastic waves. To find exact solutions to the derived equation we use the perturbation method with the Pade approximants [14]. In this paper we will restrict our search to solitary-wave (soliton-like) solutions.


## 2. The Equation Derivation

We start from the model of Timoshenko-type shell [15], which, in contrast to the Kirchhoff–Love model, takes into account the additional deformations associated with transverse forces and rotational inertia.The equations of motion of cylindrical shell element under the assumption of independence of the wave process from the circumferential coordinate have the form(1)
∂
2
u
∂
x
2
-
μ
k
y
∂
w
∂
x
-
γ
g
1
-
μ
2
E
∂
2
u
∂
t
2
=
0
,
k
2
1
-
μ
2
∂
2
w
∂
x
2
+
∂
ψ
∂
x
+
μ
k
y
∂
u
∂
x
-
k
y
2
w
-
γ
g
1
-
μ
2
E
∂
2
w
∂
t
2
=
k
1
w
+
k
2
w
3
,
∂
2
ψ
∂
x
2
-
γ
g
1
-
μ
2
E
∂
2
ψ
∂
t
2
-
6
k
2
1
-
μ
h
2
∂
w
∂
x
+
ψ
=
0
,where x, z are longitudinal and radial coordinates; u
,
w are displacements of points of the shell’s middle surface in the directions x, z; ψ is the angle of rotation of a normal to the middle surface, E is Young’s modulus, μ is Poisson’s ratio, k
y
=
1
/
R is the parameter of curvature, R is the radius of curvature of the shell, γ is the specific weight of the shell material, g is the gravitational acceleration, t is time, h is the shell thickness, and k is the correction factor of the Timoshenko model; and k
1, k
2 are the coefficients characterizing the influence of external nonlinear elastic medium (Figures 1 and 2).Figure 1
Geometry of cylindrical shell.Figure 2
Displacements and rotation of shell element.Solutions of system (1) will be sought in the class of stationary waves:(2)
u
=
u
ξ
,
w
=
w
ξ
,
ψ
=
ψ
ξ
,where ξ
=
x
-
V
t is “running coordinate” and V is the constant wave speed. Equations (1) take the form(3)
1
-
k
4
d
2
u
d
ξ
2
-
μ
k
y
d
w
d
ξ
=
0
,
(4)
k
3
d
2
w
d
ξ
2
+
d
ψ
d
ξ
+
μ
k
y
d
u
d
ξ
-
k
y
2
w
-
k
4
d
2
w
d
ξ
2
=
k
1
w
+
k
2
w
3
,
(5)
1
-
k
4
d
2
ψ
d
ξ
2
-
6
k
2
1
-
μ
h
2
d
w
d
ξ
+
ψ
=
0
,where k
3
=
k
2
(
1
-
μ
/
2
) and k
4
=
V
2
(
γ
/
g
)
(
1
-
μ
2
/
E
).The relationship between normal and longitudinal displacements is determined from (3): (6)
d
u
d
ξ
=
μ
k
y
1
-
k
4
w
.By substituting (6) into (4), we express the derivative of the normal’s rotation angle d
ψ
/
d
ξ through the component of normal displacement w:(7)
d
ψ
d
ξ
=
-
1
k
3
k
3
-
k
4
d
2
w
d
ξ
2
+
μ
2
k
y
2
1
-
k
4
-
k
y
2
-
k
1
w
-
k
2
w
3
.By differentiating equation (5) by ξ and substituting the derivative (7) into the resulting expression, we finally get a nonlinear ordinary differential equation for the component of normal displacement(8)
d
4
w
d
ξ
4
+
a
1
d
2
w
d
ξ
2
+
a
2
d
2
w
3
d
ξ
2
+
a
3
w
+
a
4
w
3
=
0
,where(9)
a
1
=
-
1
k
3
-
k
4
-
h
2
μ
2
k
y
2
+
12
k
3
k
4
h
2
1
-
k
4
+
k
y
2
+
k
1
,
a
2
=
-
k
2
k
3
-
k
4
,
a
3
=
12
k
3
1
-
k
4
k
3
-
k
4
h
2
-
μ
2
k
y
2
1
-
k
4
+
k
y
2
+
k
1
,
a
4
=
12
k
2
k
3
1
-
k
4
k
3
-
k
4
h
2
.Note that the equation similar to (8) was previously obtained in [16] for bending waves in Timoshenko-type beam lying on a nonlinear elastic foundation. The qualitative analysis of this equation was carried out for the case a
3
=
a
4
=
0. The full equation is considered in this paper.


## 3. Solitary-Wave Solutions

To find exact solitary-wave solutions of equation (8) we will use the direct perturbation method with the Pade approximants [14], proposed in [10]. The basic idea of the method consists in establishing the relationships between the equation coefficients and parameters of the sought solution in which the perturbation series becomes geometric. In this case the sum of the series gives the expression for the sought solution. The criterion of geometricity of the series is the equality of sequential diagonal Pade approximants [
Q
/
Q
]
,
[
Q
+
1
/
Q
+
1
]
,
[
Q
+
2
/
Q
+
2
]
,
…, the minimum order Q of which is determined by the pole order of the sought solution.The pole order of the solution is determined from the balance of the dominant terms of equation [17]. By substituting the function w
=
a
z
-
p, where a and p are constant, in the left part of (8), we obtain the expression(10)
a
p
p
+
1
p
+
2
p
+
3
z
-
p
-
4
+
3
a
2
a
3
p
3
p
+
1
z
-
3
p
-
2
+
a
1
a
p
p
+
1
z
-
p
-
2
+
a
3
a
z
-
p
+
a
4
a
3
z
-
3
p
,the first two terms of which are dominant. The balance is achieved by the equality of degrees -
p
-
4
=
-
3
p
-
2 from which we obtain p
=
1. Thus, the solution has a pole of the 1st order and it is necessary to choose Q
=
1 for the Pade approximants.According to the perturbation method [8], we seek a solution of equation (8) in the form of functional series(11)
w
=
∑
n
=
1
∞
ε
n
w
n
ξ
,where w
n
ξ are functions to be determined and ε is a formal parameter. By substituting (11) into (8) and grouping in powers of ε, we obtain the infinite system of equations for the functions w
n
ξ: (12)
ε
1
:
d
4
w
1
d
ξ
4
+
a
1
d
2
w
1
d
ξ
2
+
a
3
w
1
=
0
,
ε
2
:
d
4
w
2
d
ξ
4
+
a
1
d
2
w
2
d
ξ
2
+
a
3
w
2
=
0
,
ε
3
:
d
4
w
3
d
ξ
4
+
a
1
d
2
w
3
d
ξ
2
+
a
3
w
3
=
-
a
4
w
1
3
-
3
a
2
w
1
w
1
d
2
w
1
d
ξ
2
+
2
d
w
1
d
ξ
2
,
ε
4
:
d
4
w
4
d
ξ
4
+
a
1
d
2
w
4
d
ξ
2
+
a
3
w
4
=
-
3
a
4
w
1
2
w
2
-
3
a
2
2
w
1
w
2
d
2
w
1
d
ξ
2
+
w
1
2
d
2
w
2
d
ξ
2
+
4
w
1
d
w
1
d
ξ
d
w
2
d
ξ
+
2
w
2
d
w
1
d
ξ
2
,
⋮The first equation of system (12) has particular solution w
1
=
e
ξ under the condition of(13)
a
1
=
-
a
3
-
1
.The equality w
2
=
0 must be chosen as a solution to the second equation of system (12). This allows us to obtain the solution of the nth equation of the system in the form of w
n
=
K
n
w
1
n. By sequentially solving the other equations and determining the constants K
n for decomposition (11), we have(14)
w
=
ε
e
ξ
+
9
a
2
+
a
4
8
a
3
-
9
ε
e
ξ
3
+
9
a
2
+
a
4
25
a
2
+
a
4
64
a
3
-
9
a
3
-
25
ε
e
ξ
5
+
9
a
2
+
a
4
49
a
2
+
a
4
17
a
2
a
3
+
a
3
a
4
-
225
a
2
-
17
a
4
512
a
3
-
9
2
a
3
-
25
a
3
-
49
ε
e
ξ
7
+
⋯
.The right part of equality (14) contains only odd degree of product ε
e
ξ. By multiplying both parts of the equality by ε
e
ξ and denoting z
=
ε
e
ξ
2, we obtain a power series containing all powers of the variable z:(15)
w
~
=
z
+
9
a
2
+
a
4
8
a
3
-
9
z
2
+
9
a
2
+
a
4
25
a
2
+
a
4
64
a
3
-
9
a
3
-
25
z
3
+
9
a
2
+
a
4
49
a
2
+
a
4
17
a
2
a
3
+
a
3
a
4
-
225
a
2
-
17
a
4
512
a
3
-
9
2
a
3
-
25
a
3
-
49
z
4
+
⋯
.By calculating the diagonal Pade approximants [1/1] and [2/2] for series (15), we perform the factorization of their difference. The numerator of the resulting expression contains the multipliers(16)
a
3
-
49
9
a
2
+
a
4
a
2
a
3
+
a
4
.The first multiplier (16) is in the denominators of (15) and cannot go to zero. The second multiplier is included in the numerators of all terms (15), except the first, and leads to a trivial solution w
~
=
z. Finally, the third multiplier is reset to zero under the condition of(17)
a
3
=
-
a
4
a
2
,which transforms the series (15) in the geometric one(18)
w
~
=
z
-
a
2
8
z
2
+
a
2
2
64
z
3
-
a
2
3
512
z
4
+
⋯
.The sum of the series (18) coincides with approximants [1/1], [2/2] and after returning to the variables ξ and w, we have an equation(19)
w
=
8
ε
e
ξ
a
2
ε
2
e
ξ
2
+
8
=
±
2
a
2
sech
⁡
ξ
+
ln
⁡
±
ε
a
2
8
,in which the signs are chosen the same as the sign of ε. The expression (19) contains two arbitrary parameters ε, a
2 and is an exact solution of equation (8) whose coefficients are connected by the relations (13) and (17). Under the assumption that a
2
>
0, the solution has the form of a traveling pulse. Note, that unlike the classical perturbation method, where ε is a small parameter, in this case, ε is an arbitrary real number.Strictly speaking, a power series would be geometric or becomes so after regrouping the terms, if all diagonal Pade approximants[
Q
/
Q
]
,
[
Q
+
1
/
Q
+
1
]
,
[
Q
+
2
/
Q
+
2
]
,
… coincide with each other. Condition (17) guarantees equality of only the first two approximants of the specified sequence. However, the substitution in (8) shows that (19) is the exact solution. Otherwise, we would require the equality of the two following approximants [
Q
+
1
/
Q
+
1
] and [
Q
+
2
/
Q
+
2
] to obtain the second equation for the coefficients entering the decomposition (15) and to check whether these approximants give exact solution. Such steps are repeated until an exact solution is found, either the system of equations for the coefficients will be inconsistent. In the latter case, the solution by this method is impossible to obtain.To find a kink-type solution we need to substitute the expression(20)
w
=
E
+
∑
n
=
1
∞
ε
n
w
n
ξ
,where E is an arbitrary constant, into equation (8). After grouping by powers of ε we have(21)
ε
0
:
E
3
a
4
+
E
a
3
=
0
,
ε
1
:
d
4
w
1
d
ξ
4
+
3
E
2
a
2
+
a
1
d
2
w
1
d
ξ
2
+
3
E
2
a
4
+
a
3
w
1
=
0
,
ε
2
:
d
4
w
2
d
ξ
4
+
3
E
2
a
2
+
a
1
d
2
w
2
d
ξ
2
+
3
E
2
a
4
+
a
3
w
2
=
-
6
E
a
2
d
d
ξ
w
1
d
w
1
d
ξ
-
3
E
a
4
w
1
2
,
ε
3
:
d
4
w
3
d
ξ
4
+
3
E
2
a
2
+
a
1
d
2
w
3
d
ξ
2
+
3
E
2
a
4
+
a
3
w
3
=
-
6
E
a
2
d
2
d
ξ
2
w
1
w
2
-
3
a
2
d
d
z
w
1
2
d
w
1
d
z
-
a
4
w
1
6
E
w
2
+
w
1
2
,
⋮From the first equation of system (21) we find(22)
a
4
=
-
a
3
E
2
;the second equation (21) has a particular solution w
1
=
e
ξ if(23)
a
1
=
-
3
E
2
a
2
+
2
a
3
-
1
;the remaining equations (21) have the solutions in the form of w
n
=
K
n
w
1
n. By sequentially determining constants K
n
,
n
=
2,3
,
…, we can represent the series from the right part of equality (20) in the form of(24)
∑
n
=
1
∞
ε
n
w
n
=
z
-
4
E
2
a
2
-
a
3
2
E
a
3
+
2
z
2
+
6
E
2
a
2
-
2
a
3
-
1
9
E
2
a
2
-
a
3
4
E
2
a
3
+
2
2
a
3
+
9
z
3
-
28
E
4
a
2
2
a
3
+
72
E
4
a
2
2
-
16
E
2
a
2
a
3
2
-
58
E
2
a
2
a
3
-
36
E
2
a
2
+
2
a
3
3
+
9
a
3
2
+
8
a
3
16
E
2
a
2
-
a
3
8
E
3
a
3
+
2
2
a
3
+
8
2
a
3
+
9
z
3
+
⋯
,where z
=
ε
e
ξ.The numerator of the difference of the diagonal Pade approximants [1/1] and [2/2] for the series (24) contains the multipliers(25)
a
3
+
8
11
E
2
a
2
a
3
-
18
E
2
a
2
-
4
a
3
2
+
2
a
3
2
E
2
a
2
+
1
,the first of which cannot be zeroed, the second does not lead to exact solutions, and the third gives the condition(26)
a
2
=
-
1
2
E
2
,under which the series (24) becomes geometric(27)
∑
n
=
1
∞
ε
n
w
n
=
z
+
z
2
2
E
+
z
3
4
E
2
+
z
4
8
E
3
+
⋯and has the sum(28)
2
E
z
2
E
-
z
.After substituting (28) instead of the series into (20) and returning to the variable ξ we obtain the expression, being the exact solution of equation (8) under conditions (22), (23), and (26):(29)
w
=
E
+
2
E
ε
e
ξ
2
E
-
ε
e
ξ
=
-
E
tanh
⁡
1
2
ξ
+
ln
⁡
-
ε
2
E
.The solution (29), containing two arbitrary constants E and ε, is bounded and has the form of a traveling front (kink) if E
ε
<
0.


## 4. Numerical Results

In this chapter, the verification is performed to ensure that the proposed analytical method is correct.We use finite difference method [18] to obtain numerical solutions of the nonlinear ordinary differential equation (8). Difference scheme for solving the boundary value problem can be represented by using central differences:(30)
w
j
+
2
-
4
w
j
+
1
+
6
w
j
-
4
w
j
-
1
+
w
j
-
2
h
4
+
a
1
w
j
+
1
-
2
w
j
+
w
j
-
1
h
2
+
a
2
w
j
+
1
3
-
2
w
j
3
+
w
j
-
1
3
h
2
+
a
3
w
j
+
a
4
w
j
3
=
0
.To linearize the difference equation (30), we use the following approximate replacement for the terms with cubic nonlinearity:(31)
w
k
+
1
3
=
w
k
+
1
3
-
w
k
3
+
w
k
3
=
w
k
+
1
-
w
k
·
w
k
+
1
2
+
w
k
+
1
w
k
+
w
k
2
+
w
k
3
≈
w
k
+
1
-
w
k
·
3
w
k
2
+
w
k
3
=
3
w
k
+
1
w
k
2
-
2
w
k
3
,where k denotes the iteration number. Substituting (31) into (30), we get the difference problem to be linear in the grid variable w
k
+
1:(32)
w
0
k
+
1
=
w
~
-
ξ
0
,
w
1
k
+
1
=
w
~
-
ξ
0
+
h
,
w
j
+
2
k
+
1
-
4
w
j
+
1
k
+
1
+
6
w
j
k
+
1
-
4
w
j
-
1
k
+
1
+
w
j
-
2
k
+
1
h
4
+
a
1
w
j
+
1
k
+
1
-
2
w
j
k
+
1
+
w
j
-
1
k
+
1
h
2
+
a
2
3
w
j
+
1
k
+
1
w
j
+
1
k
2
-
2
w
j
+
1
k
3
-
2
3
w
j
k
+
1
w
j
k
2
-
2
w
j
k
3
+
3
w
j
-
1
k
+
1
w
j
-
1
k
2
-
2
w
j
-
1
k
3
h
2
+
a
3
w
j
k
+
1
+
a
4
3
w
j
k
+
1
w
j
k
2
-
2
w
j
k
3
=
0
,
j
=
2
,
N
-
3
¯
,
w
N
-
2
k
+
1
=
w
~
ξ
0
-
h
,
w
N
-
1
k
+
1
=
w
~
ξ
0
,where w
n is a grid variable set on a uniform grid ξ
n
,
ξ
n
=
-
ξ
0
+
n
h
,
n
=
0
,
N
-
1
¯; N
=
1
+
2
ξ
0
/
h is a number of grid nodes; -
ξ
0
,
ξ
0 is an integration interval;h is a grid spacing;k and k
+
1 are the numbers of the previous and current iterations, respectively; and w
~
ξ is an exact solution of equation (8).The problem (32) is solved by means of fixed-point iteration method. The calculation process is stopped when L
2-norm of the difference between the solution from the current iteration and the solution from previous iteration is less than a specified small number α:(33)
w
k
+
1
-
w
k
2
≡
∑
n
=
0
N
-
1
w
n
k
+
1
-
w
n
k
2
1
/
2
<
α
.Deviations of the numerical solution from the exact solution for different values of the coefficients of equation (8) are shown in Tables 1 and 2. The calculations are performed under the conditions ξ
0
=
10
,
α
=
1.0
e
-
6.Table 1
Deviations of the numerical solutionw
n from the exact solution w
~
ξ
n (19).
h
1.00
e
-
01
5.00
e
-
02
2.50
e
-
02
1.25
e
-
02
6.25
e
-
03
a
1
=
1.0
,
a
2
=
1.0
,
a
3
=
-
2.0
,
a
4
=
2.0
,
ε
=
2.0
w
n
-
w
~
ξ
n
2
5.83
e
-
03
2.05
e
-
03
7.24
e
-
04
2.54
e
-
04
8.95
e
-
04
a
1
=
0.0
,
a
2
=
1.0
,
a
3
=
-
1.0
,
a
4
=
1.0
,
ε
=
2.0
w
n
-
w
~
ξ
n
2
5.89
e
-
03
2.07
e
-
03
7.30
e
-
04
2.71
e
-
04
4.27
e
-
04
a
1
=
-
0.5
,
a
2
=
2.0
,
a
3
=
-
0.5
,
a
4
=
2.0
,
ε
=
1.0
w
n
-
w
~
ξ
n
2
4.16
e
-
03
1.46
e
-
03
5.15
e
-
04
1.90
e
-
04
4.90
e
-
04Table 2
Deviations of the numerical solutionw
n from the exact solution w
~
ξ
n (29).
h
1.00
e
-
01
5.00
e
-
02
2.50
e
-
02
1.25
e
-
02
6.25
e
-
03
a
1
=
4.5
,
a
2
=
0.5
,
a
3
=
2.0
,
a
4
=
-
2.0
,
E
=
1.0
,
ε
=
-
1.0
w
n
-
w
~
ξ
n
2
9.19
e
-
03
3.23
e
-
03
1.14
e
-
03
3.96
e
-
04
1.57
e
-
05
a
1
=
4.5
,
a
2
=
-
0.5
,
a
3
=
2.0
,
a
4
=
-
2.0
,
E
=
1.0
,
ε
=
-
0.5
w
n
-
w
~
ξ
n
2
5.14
e
-
03
1.86
e
-
03
6.69
e
-
04
2.18
e
-
04
7.97
e
-
05
a
1
=
4.5
,
a
2
=
-
0.125
,
a
3
=
2.0
,
a
4
=
-
0.5
,
E
=
2.0
,
ε
=
-
1.5
w
n
-
w
~
ξ
n
2
1.36
e
-
02
4.83
e
-
03
1.72
e
-
03
6.60
e
-
04
1.38
e
-
03As shown in Figure3, the curves of the exact solution and the numerical solution practically coincide: the maximum value of the absolute difference of the solutions does not exceed 0
.
83
e
-
4.Figure 3
The exact solution (solid line) and error of the numerical solution multiplied by1
0
4 (dashed line).
(a)
Solution (19) of equation (8) at a
1
=
1.0
,
a
2
=
1.0
,
a
3
=
-
2.0
,
a
4
=
2.0, and ε
=
2.0
(b)
Solution (29) of equation (8) at a
1
=
4.5
,
a
2
=
0.5
,
a
3
=
2.0
,
a
4
=
-
2.0
,
E
=
1.0, and ε
=
-
1.0


## 5. Conclusion

The analysis showed that the search for solutions to the equations of motion of an element of a cylindrical shell of Timoshenko type, interacting with a nonlinear elastic medium, in the class of stationary longitudinal-bending waves, leads to a nonlinear ordinary differential equation (8) for the normal displacement. The exact solitary-wave solutions in the form of a traveling pulse and traveling front are found for this equation. The numerical solutions found by means of finite difference method are in good agreement with the corresponding exact ones.



---

*Source: 6596231-2016-11-27.xml*