On the cell-dependent vibrations and wave propagation in uniperiodic cylindrical shells

The objects of consideration are thin linearly elastic Kirchhoff–Love-type circular cylindrical shells having a periodically micro-heterogeneous structure in circumferential direction (uniperiodic shells). The aim of this contribution is to study certain problems of micro-vibrations and of wave propagation related to micro-fluctuations of displacement field caused by a periodic structure of the shells. These micro-dynamic problems will be analysed in the framework of a certain mathematical averaged model derived by means of the combined modelling procedure. The combined modelling includes both the asymptotic and the tolerance non-asymptotic modelling techniques, which are conjugated with themselves under special conditions. Contrary to the starting exact shell equations with highly oscillating, non-continuous and periodic coefficients, governing equations of the combined model have constant coefficients depending also on a cell size. Hence, this model takes into account the effect of a microstructure size on the dynamic behaviour of the shells (the length-scale effect). It will be shown that the micro-periodic heterogeneity of the shells leads to cell-depending micro-vibrations and to exponential waves as well as to dispersion effects, which cannot be analysed in the framework of the asymptotic models commonly used for investigations of vibrations and wave propagation in the periodic structures.


Introduction
Thin linearly elastic Kirchhoff-Love-type circular cylindrical shells with a periodically micro-inhomogeneous structure (a periodically varying thickness and/or periodically varying elastic and inertial properties) in circumferential direction are analysed. Shells of this kind are termed uniperiodic. At the same time, the shells have constant structure in axial direction. The shells under consideration are composed of a large number of identical elements, and every such element, called a periodicity cell, can be treated as a thin shell. It means that the period of inhomogeneity is very large compared with the maximum shell thickness and very small as compared to the midsurface curvature radius as well as the length dimension of the shell midsurface in the periodicity direction. As examples, we can mention cylindrical shells with periodically spaced stiffeners as shown in Fig. 1 or made of two kinds of periodically distributed materials as shown in Fig. 2.
It should be noted that in the general case, on the shell midsurface we deal with not periodic but locally periodic structure. By a locally periodic shell we mean a shell which, in small subregions of the shell midsurface, can be approximately regarded as periodic. Hence, a locally periodic shell is made of a large number of not identical, but similar elements. However, for cylindrical shells the Gaussian curvature is equal to zero, and hence, on the developable cylindrical surface we can separate a cell which can be referred to as the representative cell for the whole shell midsurface. It means that on cylindrical surface, we deal not with locally periodic but with a periodic structure.
Dynamic problems of periodic shells are described by partial differential equations with highly oscillating, periodic and non-continuous coefficients. Thus, these equations are too complicated to be applied to the investigations of engineering problems. To obtain averaged equations with constant coefficients, various approximate modelling procedures for shells of this kind have been proposed. Periodic cylindrical shells (plates) are usually described using homogenized models derived by means of asymptotic methods. These models from a formal  [1], where asymptotic modelling of plates, laminates and shells is discussed. Unfortunately, in the models of this kind the effect of a periodicity cell length dimensions (called the length-scale effect) on the overall shell behaviour is neglected.
This effect can be taken into account using the modified couple stress-based theories of continuous media. We mention here paper by Awrejcewicz et al. [2], where mathematical model for the analysis of static and dynamic problems of micro/nano-beams is derived and discussed; the size-dependent model equations are formulated on the basis of the Grigolyuk-Chulkov hypotheses and the modified couple stress theory.
The length-scale effect can be bearing in mind applying the multiscale-multifield models derived from a non-classical (generalized) continuum formulation. These models take into account the microstructure size by means of microdisplacement variables added to the standard macrodisplacements and of material internal length parameters, cf., e.g. Settimi et al. [3], where the length-scale effect on the dynamic properties of a composite microcracked elastic bar is studied; the internal parameters represent here density and length of microcracks.
Some numerical approaches are also proposed to study the size effects in mechanical problems for microheterogeneous structures. As example, we can mention here the paper by Hassani et al. [4], where the sizedependent variational differential quadrature procedure is combined with the finite element method into a new technique. In this paper, by considering several numerical examples, it has been shown that the proposed size-dependent formulation and numerical solution approach have a good performance to study the large deformations of hyperelastic microstructured bodies.
The length-scale effect can be also taken into account using the non-asymptotic-tolerance averaging technique, cf. Woźniak and Wierzbicki [5], Woźniak et al. [6,7]. Some applications of this method to the modelling of mechanical and thermomechanical problems for various periodic structures are shown in many works. The extended list of publications on this topic can be found in [5][6][7]. We mention here monograph by Tomczyk [8], where the length-scale effect in dynamics and stability of periodic cylindrical shells is investigated, paper by Marczak and Jędrysiak [9], where vibrations of periodic three-layered plates with inert core are studied and papers by Tomczyk and Litawska [10,11], where certain extended co-called general tolerance and general asymptotic-tolerance models for the analysis of dynamic problems for periodic cylindrical shells are proposed and discussed. These general models are derived by means of a certain extended version of the tolerance modelling technique presented by Tomczyk and Woźniak [12].
In the last years, the tolerance modelling was adopted for mechanical and thermomechanical problems of functionally graded structures, e.g. for heat conduction in longitudinally graded structures by Ostrowski and Michalak [13], for thermoelasticity of transversally graded laminates by Pazera and Jędrysiak [14], for vibrations of functionally graded thin plates by Wirowski [15], for dynamics of transversally and longitudinally graded thin cylindrical shells by Tomczyk and Szczerba [16][17][18].
The aim of this note is to study certain problems of cell-depending vibrations and of long wave propagation related to micro-fluctuations of displacement field caused by a periodic structure of the shells. Note that we deal with long waves if condition λ/L << 1 holds, where λ is the characteristic length dimension of the cell and Lis the wavelength. These micro-dynamic problems will be analysed in the framework of the combined asymptotic-tolerance model proposed in [8]. Governing equations of this averaged model have constant coefficients depending also on a microstructure size. An important advantage of this model is that it makes it possible to separate the macroscopic description of the modelling problem from its microscopic description. It will be shown that the periodic micro-heterogeneity of the shells leads to vibrations depending on a cell size and to exponential waves as well as to dispersion effects, which cannot be analysed in the framework of the asymptotic models commonly used for investigations of vibrations and wave propagation in the periodic shells under consideration. The new wave propagation speed depending on a cell size will be obtained and analysed.
It should be noted that this article is a certain continuation of papers by Tomczyk and Szczerba [16][17][18] and by Tomczyk and Litawska [10,11], in which some special dynamic problems for thin functionally graded cylindrical shells [16][17][18] and for shells with two-directional periodic structure in directions tangent to the shell midsurface (biperiodic shells) [10,11] are analysed by applying the tolerance modelling technique. Note that in the non-asymptotic-tolerance approach, shells with one-directional periodic structure (uniperiodic shells) being objects of consideration in this work are not special cases of biperiodic shells. Model equations for uniperiodic shells are more complicated than those for biperiodic shells and contain a lot of length-scale terms which do not have counterparts in the equations for biperiodic shells. The occurrence of these terms is strictly related to the fact that the modelling physical reliability conditions for uniperiodic shells are less restrictive than those for biperiodic shells.
The periodic shells being objects of considerations in this paper are widely applied in civil engineering, most often as roof girders and bridge girders. They are also widely used as housings of reactors and tanks. Periodic shells having small length dimensions are elements of air-planes, ships and machines.

Formulation of the problem: starting equations
We assume that x 1 and x 2 are coordinates parametrizing the shell midsurface Min circumferential and axial directions, respectively. We denote x ≡ x 1 ∈ ≡ (0, L 1 ) and ξ ≡ x 2 ∈ ≡ (0, L 2 ), where L 1 , L 2 are length dimensions of M,cf. Figs. 1 and 2. Let Ox 1x 2x 3 stand for a Cartesian orthogonal coordinate system in the physical space R 3 and denotex ≡ ( It means that on M, the orthonormal parametrization is introduced. Note that derivative ∂r/∂ x α , α = 1, 2, should be understood as differentiation of each component ofr, Sub-and superscripts α, β,… run over 1, 2 and are related to x 1 , x 2 , summation convention holds. Partial differentiation related to x α is represented by Let a αβ and b αβ stand for the midsurface first and second metric tensors, respectively. Under orthonormal parametrization introduced on M, a αβ is a unit tensor and components of tensor b αβ are: b 22 . Differentiation with respect to time is represented by the overdot. Let d(x), r stand for the shell thickness and the midsurface curvature radius, respectively.
The basic cell and an arbitrary cell (x) with the centre at point x ∈ are defined by means of: The microstructure length parameter λ satisfies conditions: λ/d max >> 1, λ/r << 1 and λ/L 1 << 1.
Setting z ≡ z 1 ∈ [−λ/2, λ/2], we assume that the cell has a symmetry axis for z = 0. It is also assumed that inside the cell the geometrical, elastic and inertial properties of the shell are described by even functions of argument z.
Denote by u α = u α (x, ξ, t), w = w(x, ξ, t), (x, ξ, t) ∈ × × I, the shell displacements in directions tangent and normal to M, respectively. Elastic properties of the shells are described by shell stiffness tensors D αβγδ (x), B αβγδ (x). Let μ(x) stand for a shell mass density per midsurface unit area. The external forces will be neglected.
The considerations are based on the well-known Kirchhoff-Love theory of thin elastic shells, cf. Kaliski [19].
In the framework of the shell theory under consideration, strain energy function E(x, ξ, t), (x, ξ, t) ∈ × × I, related to midsurface Mhas the form where the membrane ε αβ (x, ξ, t) and curvature κ αβ (x, ξ, t) strain tensors are The kinetic energy function K (x, ξ, t), (x, ξ, t) ∈ × × I, related to midsurface M is given by Let us introduce the action functional with lagrangian L being a highly oscillating function with respect to x. Here, under assumption that the external forces are neglected, lagrangian L have the form where kinetic energy K and strain energy Eare given above. Substituting (1)-(3) into (5) and taking into account that b 22 = b 12 = b 21 = 0 and b 11 = −r −1 , we arrive at Lagrange function (5) in the form Under assumption that ∂ L/∂(∂ β u α ) and ∂ L/∂(∂ αβ w) are continuous, from the principle of stationary action applied to A(u α , w), we obtain the system of Euler-Lagrange equations, which can be written in an explicit form as Equations (7) coincide with the well-known governing equations of Kirchhoff-Love theory of thin elastic shells, cf. Kaliski [19]. For periodic shells, coefficients D αβγδ (x), B αβγδ (x), μ(x) of (6) and (7) are highly oscillating, non-continuous and periodic functions in x. Applying the combined asymptotic-tolerance modelling technique, the averaged model equations with constant coefficients depending also on a cell size were derived in [8]. The combined modelling given in the general form in Woźniak et al. [7] includes both the asymptotic and the tolerance non-asymptotic modelling techniques, which are conjugated with themselves under special conditions. Contrary to the starting exact shell equations with highly oscillating, non-continuous and periodic coefficients, governing equations of the combined model have constant coefficients depending also on a cell size. Hence, this model takes into account the effect of a microstructure size on the dynamic behaviour of the shells (the length-scale effect). Here, the combined model equations formulated in [8] will be used for investigations of two special micro-dynamic problems for cylindrical shell having constant thickness and made of two component linearly elastic isotropic materials periodically distributed in circumferential direction as shown in Fig. 2. The first of these problems deals with cell-depending micro-vibrations caused by a periodic structure of the shell. The second one deals with propagation of the long waves related to micro-fluctuations of axial displacements. It has to be emphasized that these aforementioned special micro-dynamic problems can be studied in the framework of neither the asymptotic models nor the known commercial numerical models based on the finite element method.
To make the analysis more clear, in the next section the asymptotic-tolerance model for the shells under consideration will be reminded, following [8]. Moreover, the basic concepts and assumptions of the tolerance modelling technique and of the consistent asymptotic approach will be outlined, following [7,8].

Modelling procedure: asymptotic-tolerance model
The combined modelling technique under consideration is realized in two steps. The first step is based on the consistent asymptotic modelling procedure [7,8]. The second one is realized by means ofthe tolerance non-asymptotic modelling technique [7,8].

Step 1. Consistent asymptotic modelling
The fundamental concepts of the consistent asymptotic procedure are those of an averaging operation and fluctuation shape functions. In what follows, the above-mentioned concepts will be specified with respect to one-dimensional region ≡ (0, L 1 ) defined in this contribution.
Let f (x) be a function defined in¯ ≡ [0, L 1 ], which is integrable and bounded in every cell (x), x ∈ . The averaging operation of f (·) is defined by It can be seen that if f (·) is -periodic, then < f > is constant. Denote by ∂ k 1 the kth derivative of function defined in . Let h(x) be a λ-periodic, highly oscillating function defined in¯ = [0, L 1 ], which is continuous together with derivatives ∂ k 1 h, k = 1, . . . , R − 1, and has a continuous or piecewise continuous bounded derivative ∂ R 1 h. Function h(·) will be called the fluctuation shape function of the Rth kind, is a shell mass density. Nonnegative integer R is assumed to be specified in every problem under consideration.
The first step of the combined modelling is based on the consistent asymptotic averaging of lagrangian (6). To this end, we shall restrict considerations to displacement fields Unknown functions u 0 α , U α in (9) are assumed to be continuous and bounded in together with their first derivatives. Unknown functions w 0 , W in (9) are assumed to be continuous and bounded in together with their derivatives up to the second order. Unknowns u 0 α , w 0 and U , α W are called macrodisplacements and fluctuation amplitudes, respectively. They are independent of ε. This is the main difference between the asymptotic approach under consideration and approach which is used in the known homogenization theory, cf. Bensoussan et al. [20], Jikov et al. [21].
By (9) are denoted λ-periodic highly oscillating fluctuation shape functions depending on ε. The fluctuation shape functions are assumed to be known in every problem under consideration. They have to satisfy conditions: We substitute the right-hand sides of (9) into (6) and take into account that under limit passage ε → 0, terms depending on ε can be neglected and every continuous and bounded function of argument z ∈ ε (x) tends to function of argument x ∈¯ . Moreover, if ε → 0 then by means of a property of the mean value, cf. [21], the obtained result tends weakly to the function being the averaged form of starting lagrangian (6) under consistent asymptotic decomposition (9). Then, applying the principle of stationary action we obtaingoverning equations of the consistent asymptotic model for the unperiodic shells under consideration. These equations consist of partial differential equations for macrodisplacements u 0 α , w 0 coupled with linear algebraic equations for fluctuation amplitudes U α , W . After eliminating fluctuation amplitudes from the governing equations by means of where where Since displacement fields u α (x, ξ, t), w(x, ξ, t) have to be uniquely defined in × × I, we conclude that u α (x, ξ, t), w(x, ξ, t) have to take the form with U α , W given by (10).
given by (12) are tensors of effective elastic moduli for the considered composite uniperiodic shells.
In contrast to starting equations (7) with discontinuous, highly oscillating and periodic coefficients, the asymptotic model equations (11) have coefficients constant but independent of the microstructure size λ. Hence, the above model is not able to describe the length-scale effect on the overall shell dynamics. That is why, the model derived in the first step of combined modelling is referred to as the macroscopic model for the problem under consideration.
Unknown macrodisplacements u 0 α , w 0 and fluctuation amplitudes U α , W must be continuously bounded in .
The resulting equations (11) are uniquely determined by the postulated a priori periodic fluctuations shape , representing oscillations inside a cell. These functions can be derived from the periodic discretization of the cell using, for example, the finite element method or obtained as exact or approximate solutions to certain periodic eigenvalue problems on the cell describing free periodic vibrations. If the fluctuation shape functions are not derived as solutions to periodic eigenvalue cell problems mentioned above, then the effective moduli (12) of the shell are obtained without specification of the periodic cell problems. This situation is different from that occurring in the known asymptotic homogenization approach, cf., e.g. [20], where only solutions to the periodic cell problems make it possible to define the effective moduli of the structure under consideration.
In the first step of combined modelling, it is assumed that within the asymptotic model, solutions u 0 α , w 0 to the problem under consideration are known. Hence, there are also known functions where U α , W are given by means of (10).

Step 2. Tolerance modelling
The second step of the combined modelling is based on the tolerance modelling technique, cf [7,8].
The fundamental concepts of the tolerance modelling procedure under consideration are those of two tolerance relations between points and real numbers determined by tolerance parameters, slowly varying functions, tolerance-periodic functions, fluctuation shape functions and the averaging operation.
In what follows, some of the above-mentioned concepts and assumptions will be specified with respect to one-dimensional region ≡ (0, L 1 ) defined in this contribution.
Let F(x) be a function defined in¯ = [0, L 1 ], which is continuous, bounded and differentiable in¯ together with their derivatives up to the Rth order. Note that function Fcan also depend on ξ ∈¯ = [0, L 2 ] and time coordinate t as parameters. Let δ ≡ (λ, δ 0 , δ 1 , . . . , δ R ) be the set of tolerance parameters. The first of them is related to the distances between points in , the second one is related to the distances between values of function F(·) and the kth one to the distances between values of the kth derivative of F(·), k = 1, . . . , R. A function F(·) is called slowly varying of the Rth kind with respect to cell and tolerance parameters δ, F ∈ SV R δ ( , ), if and only if the following two conditions are satisfied where symbols " λ ≈" and " Roughly speaking, from (15) and (16) it follows that slowly varying function F(·) can be treated as constant on an arbitrary cell and that the products of derivatives of slowly varying function in periodicity direction andmicrostructure length parameter λ are treated as negligibly small.
An integrable and bounded function f (x) defined in¯ = [0, L 1 ], which can also depend on ξ ∈¯ and time coordinate t as parameters, is called tolerance-periodic with respect to cell and tolerance parameters . For function f (·) being tolerance-periodic together with its derivatives up to the Rth order, we shall write f ∈ T P R δ ( , ), δ ≡ (λ, δ 0 , δ 1 , . . . , δ R ). The concepts of fluctuation shape functions and averaging operation have been explained in Sect. 3.1. The tolerance modelling is based on two assumptions. The first assumption is called the tolerance averaging approximation. The second one is termed the micro-macro decomposition.
Let f (·) be an integrable periodic function defined in¯ = [0, The tolerance averaging approximation has the form and In the course of modelling, terms O(δ) in (17) and (18) are neglected. Let us observe that the slowly varying functions can be regarded as invariant under averaging. Approximations given above are applied in the modelling problems discussed in this contribution. For details, the reader is referred to [5][6][7][8].
The second fundamental assumption, called the micro-macro decomposition, states that the displacements fields occurring in the starting lagrangian under consideration can be decomposed into unknown averaged (macroscopic) displacements being slowly varying functions in x ∈ and highly oscillating fluctuations represented by the known highly oscillating λ-periodicfluctuation shape functions multiplied by unknownfluctuation amplitudes (microscopic variables) slowly varying in x.
In the second step of combined modelling, we introduce the extra micro-macro decomposition superimposed on the known solutions u 0α , w 0 obtained within the macroscopic model.
where fluctuation (microscopic) amplitudes Q α , V are the new slowly varying unknowns, i.e.
are the new periodic, continuous and highly oscillatingfluctuation shape functions which are assumed to be known in every problem under consideration. These functions have to satisfy conditions: We substitute the right-hand sides of (19) into (6). The resulting lagrangian is denoted by L cb . Then, we average L cb over cell using averaging formula (8) and applying the tolerance averaging approximation (17), (18). As a result, we obtain function < L cb > called the tolerance averaging of starting lagrangian (6) in under superimposed decomposition (19). Next, applying the principle of stationary action, under the extra approximation 1 + λ/r ≈ 1, we arrive at the system of Euler-Lagrange equations for Q α , V , which can be written in an explicit form as Equations (20) and (21) together with the micro-macro decomposition (19) constitute the superimposed microscopic model. Coefficients of the derived model equations are constant, and some of them depend on a cell size λ (the underlined terms). The right-hand sides of (20) and (21) are known under assumption that u 0α , w 0 were determined in the first step of modelling. The basic unknowns Q α , V of the model equations must be the slowly varying functions in periodicity directions. This requirement can be verified only a posteriori, and it determines the range of the physical applicability of the model. The boundary conditions for Q α , V should be defined only on boundaries ξ = 0, ξ = L 2 . It can be shown that under assumption that fluctuation shape functions h(x), g(x) of macroscopic model coincide with fluctuation shape functions c(x), b(x) of microscopic model, we can obtain microscopic model equations (20), (21), in which c(x) and b(x) are replaced by h(x) and g(x), respectively, and in which the right-hand sides are equal to zero. Moreover, taking into account a symmetric form of tensor D αβγδ we arrive finally at three equations for unknown fluctuation amplitudes Q 1 (x, ξ, t), Q 2 (x, ξ, t) and V (x, ξ, t), which are not conjugated with themselves Equations (22)-(24) are independent of the solutions u 0α , w 0 obtained in the framework of the macroscopic model. Hence, they describe selected problems of the shell micro-dynamics (e.g. the free micro-vibrations, propagation of waves related to the micro-fluctuation amplitudes) independently of the shell macro-dynamics. Moreover, micro-dynamic behaviour of the shell in the axial and circumferential directions can be analysed independently of its micro-dynamic behaviour in the direction normal to the shell midsurface. Microscopic model equations (22)-(24) also describe certain time-boundary and space-boundary phenomena strictly related to the specific form of initial and boundary conditions imposed on unknown fluctuation amplitudes Q α , V . That is why, these equations are referred to as the boundary layer equations, where the term "boundary" is related both to time and space.
Since equations (22)-(24) are not conjugated with themselves, the micro-dynamic behaviour of the shells in the axial, circumferential and normal directions can be investigated independently of each other.

Combined asymptotic-tolerance model
Summarizing results obtained in Sects. 3.1 and 3.2, we conclude that the combined asymptotic-tolerance model of selected dynamic problems for the uniperiodic shells under consideration presented here following Tomczyk [8] is represented by: where functions u 0 α , U α , w 0 , W have to be obtained in the first step of combined modelling, i.e. in the framework of the consistent asymptotic modelling.
Coefficients of all equations derived in the framework of combined modelling are constant in contrast to coefficients in starting Eq. (7) which are discontinuous, highly oscillating and periodic in x. Moreover, some of them depend on a cell size λ. Thus, the combined model can be applied to the analysis of many phenomena caused by the length-scale effect.
Under special conditions imposed on the fluctuation shape functions, we can obtain microscopic model equations (22)-(24), which are independent of the solutions obtained in the framework of the macroscopic model. It means that an important advantage of the combined model is that it makes it possible to separate the macroscopic description of some special dynamic problems from the microscopic description of these problems.
For details, the reader is referred to Tomczyk [8].
It should be noted that the combined asymptotic-tolerance model of dynamic problems for cylindrical shells with periodic structure in both circumferential and axial directions (biperiodic shells) proposed in Tomczyk and Litawska [11] cannot be applied for analysis of dynamic problems for uniperiodic shells considered here. Model presented in [11] is derived in the framework of the extended version of the tolerance modelling technique based on a new notion of weakly slowly varying function, cf. Tomczyk and Woźniak [12]. For this function, restrictive condition (16) and approximations (18) do not hold. Moreover, in the non-asymptotic-tolerance approach, the uniperiodic shells are not special cases of shells with two-directional periodic structure.
It should be also noted that the combined asymptotic-tolerance models for functionally graded cylindrical shells are presented by Tomczyk and Szczerba in [17,18]. Coefficients of governing equations of these models are not constant. They are smooth and slowly varying either in circumferential direction [17] or in the axial one [18].
Some applications of micro-dynamic Eqs. (22)-(24) will be shown in the next section.

Examples of applications
In this section, we shall investigate two special micro-dynamic problems applying Eqs. (22)-(24). The first of them deals with free cell-depending micro-vibrations. The second one deals with propagation of the waves related to micro-fluctuations of axial displacements. It has to be emphasized that these aforementioned special micro-dynamic problems can be studied in the framework of neither the asymptotic models nor the known commercial numerical models for the periodic shells under consideration.

Formulation of the problem
The object of considerations is a thin cylindrical shell with L 1 , L 2 , r , d as its circumferential length, axial length, midsurface curvature radius and constant thickness, respectively. The shell has a periodically heterogeneous structure along circumferential direction and constant structure in the axial direction. It is assumed that the shell is made of two homogeneous elastic isotropic materials, which are perfectly bonded on interfaces, cf. Fig. 2. The free micro-vibration problem will be studied for an open simply supported shell, i.e. for a shell with hinged edges and with supports free to move, cf. [19]. The wave propagation problem will be investigated for a closed shell (obviously, in this case L 1 = 2πr ). Moreover, we assume that L 2 ≥ L 1 .
The shell's mass density per midsurface unit area μ(x) and stiffness tensors D αβγδ (x), B αβγδ (x) are described by functions λ-periodic in x and independent of ξ.
The basic cell is defined by: Figs. 2 and 3. We recall that the microstructure length parameter λ has to satisfy conditions: λ/d >> 1, λ/r << 1 and λ/L 1 << 1. Setting z ≡ z 1 ∈ [−λ/2, λ/2], we assume that the cell has a symmetry axis for z = 0. Inside the cell, the geometrical, elastic and inertial properties of the shell are described by symmetric (i.e. even) functions of argument z.
Properties of the component materials are described by: Young's moduli E 1 , E 2 , Poisson's ratios ν 1 , ν 2 and mass densities ρ 1 , ρ 2 , cf. Fig. 3. It is assumed that elastic E(·) and inertial ρ(·) properties of the composite shell are periodic functions in x, but Poisson's ratio ν ≡ ν 1 = ν 2 is constant. Inside the cell, functions E(z), ρ(z) take the form The fluctuation shape functions h(z) ∈ F S 1 ( , ), g(z) ∈ F S 2 ( , ) describe the expected form of displacement disturbances caused by a periodic structure of the shell. It means that they should approximate the expected principal modes of the shell's free vibrations. These modes have to be λ-periodic, and their mean values in every cell must be equal to zero. On the basis of knowledge of the physically reasonable approximations of principal modes of free vibrations in thin Kirchhoff-Love-type periodic cylindrical shells, cf. Tomczyk [8], and also in thin Kirchhoff-type periodic plates, cf. Jędrysiak [22], in the problem under consideration the fluctuation shape functions can be taken as: h(z) = λ sin(2πz/λ), g(z) = λ 2 [cos(2πz/λ) + c], z ∈ (x), x ∈ , where constant c, calculated from condition < μg >= 0, is equal to c = −(ρ 1 −ρ 2 ) sin(ηπ)[π(ηρ 1 +(1−η)ρ 2 )] −1 .
The subsequent analysis will be based on Eqs. (22)-(24) describing the shell micro-dynamics.

Discussion of computational results
On the basis of results shown in Figs. 4, 5, 6, 7, 8 and 9, the following conclusions can be formulated: 1. Values of the dimensionless frequencies decrease with the increasing of ratio λ/L 1 , i.e. with the decreasing of differences between period length λ and the length dimension L 1 of the shell midsurface in periodicity direction, cf. Figs. 4, 6, 8.  Fig. 5, , cf. Fig. 7, , cf. Fig. 9, are obtained for pair of ratios (E 2 /E 1 = 1, ρ 2 /ρ 1 = 0.01), i.e. for a periodic shell with a very strong inertial heterogeneity and with   Fig. 9 Diagrams of dimensionless frequencies of transversal free micro-vibrations versus: a ratio E 2 /E 1 and b ratio ρ 2 /ρ 1 ; λ/L 1 = 0.1

Long wave propagation problem
Now, let the shell under consideration be closed and unbounded along the axial coordinate ξ. We shall analyse the long wave propagation problem. We deal with long waves if condition λ/L << 1 holds, where λ is the characteristic length dimension of the cell and L is the wavelength. The waves related to micro-fluctuation amplitude Q 2 are taken into account. Hence, Eq. (23) describing the shells' micro-dynamics in an axial direction will be applied. Let the investigated problem be rotationally symmetric with a period λ/r ; thus, unknown Q 2 (·, t) of (23) is independent of x. We look for solution to Eq. (23) in the form Q 2 (ξ, t) = F(ξ − ct), where c is the wave propagation velocity. Settingh = λ −1 h, from Eq. (23) we obtain where speedsc andc are defined bỹ Equation (35) where under assumption that L >> λ the following condition holds kλ = 2πλ/L << 1.
The above equation describes the effect of dispersion. It can be seen that for kλ → 0, the dispersion effect disappears. From Eq. (37), it follows that the dispersive long waves related to micro-fluctuation amplitude Q 2 (ξ, t) can propagate across the unbounded uniperiodic shells under consideration with propagation speed depending on microstructure size λ. Note that for homogeneous isotropic shells, expression (38) leads to the well-known results c 2 = D/μ, D = Eδ/(1 − ν 2 ), where E, ν, δ, μ are Young's modulus, Poisson's ratio, the shell thickness and mass density of the shell material, respectively, cf. Kaliski [19].

Numerical results
Let us define the following dimensionless wave propagation speed where speed c is determined by formula (38).

Discussion of analytical and computational results
It was shown that the tolerance-periodic heterogeneity of the shells leads to exponential waves and to dispersion effects, which cannot be analysed in the framework of the asymptotic models for periodic shells. Moreover, the new wave propagation speed depending on the microstructure size has been obtained, cf. formula (38).
On the basis of results shown in Figs. 10, 11 and 12, the following conclusions can be formulated: i.e. with the decreasing of differences between elastic properties of the shell component materials, cf. Fig. 11, but they decrease with the increasing of ratio ρ 2 /ρ 1 ∈ [0.01, 1], i.e. with the decreasing of differences between inertial properties of the component materials, cf. Fig. 12. 4. The highest values of dimensionless speed C, cf. Figs. 11 and 12, are obtained for pair of ratios (E 2 /E 1 = 1, ρ 2 /ρ 1 = 0.01), i.e. for a periodic shell with a very strong inertial heterogeneity and with elastic homogeneous structure. The smallest values of this speed are obtained for pair of ratios (E 2 /E 1 = 0.01, ρ 2 /ρ 1 = 1), i.e. for a periodic shell with a very strong elastic heterogeneity and with inertial homogeneous structure. 5. For a homogeneous isotropic shell, expression (38) leads to result: For ν = 0.3, the value of speed C hom is equal 1.05. Comparing this result with results shown in Figs. 11 and 12, we conclude that in the unbounded homogeneous isotropic shell, the displacement wave propagates along axial direction with speed which is much smaller, i.e. about 70 times smaller, than the smallest velocity obtained for the periodic shell with a very strong elastic heterogeneity (E 2 /E 1 = 0.01) and at the same time with inertial homogeneous structure (ρ 2 /ρ 1 = 1).

Final remarks and conclusions
The following remarks and conclusions can be formulated: • Thin linearly elastic Kirchhoff-Love-type circular cylindrical shells having a periodic microstructure in circumferential direction (uniperiodic shells) are objects of consideration, cf. Figs. 1 and 2. At the same time, the shells have constant geometrical and material properties in axial direction. • The new averaged combined asymptotic-tolerance model for the analysis of selected dynamic problems for the uniperiodic cylindrical shells under consideration was derived in Tomczyk [8]. Here, the governing equations of this model are recalled and applied for investigations of certain micro-dynamic problems for the shells under consideration. The aforementioned model equations consist of macroscopic (asymptotic) model equations (11) for macrodisplacements u 0 α (x, ξ, t), w 0 (x, ξ, t), (x, ξ) ∈ × , t ∈ I, derived by means of the consistent asymptotic procedure, cf. Woźniak et al. [7], and of microscopic tolerance (nonasymptotic) model equations (20), (21) for fluctuation amplitudes Q α (x, ξ, t), V (x, ξ, t) formulated by applying the tolerance modelling technique, cf. Woźniak et al. [7]. The tolerance modelling is based on the concept of tolerance relations between points and real numbers related to the accuracy of the performed measurements and calculations. The tolerance relations are determined by the tolerance parameters. Macroand microscopic models are combined together under assumption that in the framework of the asymptotic model the solutions (14) to the problem under consideration are known. Contrary to the starting wellknown governing Eq. (7) of Kirchhoff-Love theory with highly oscillating, non-continuous and periodic coefficients, equations of the asymptotic-tolerance model have constant coefficients depending also on a microstructure size. Hence, this model allows us to describe the effect of a length scale on the dynamic shell behaviour. The resulting combined model equations are uniquely determined by the highly oscillating periodic fluctuation shape functions describing oscillations inside the cell. These functions have to be known in every problem under consideration. Under special conditions imposed on the fluctuation shape functions, we can derive microscopic equations (22)-(24), which are independent of solutions (14) obtained within the macroscopic model. It means that an important advantage of this model is thatit makes it possible to separate the macroscopic description of some special problems from their microscopic description. Moreover, Eqs. (22)-(24) involve terms with time and spatial derivatives of fluctuation amplitudes. Hence, these equations describe certain time-boundary-layer and space-boundary-layer phenomena strictly related to the specific form of initial and boundary conditions imposed on the unknown fluctuation amplitudes.
• The main aim of this contribution was to apply micro-dynamic equations (22)-(24), proposed in [8] and recalled here, to study two special micro-dynamic problems for a certain cylindrical shell made of two homogeneous elastic isotropic component materials densely and periodically distributed in circumferential direction, cf. Fig. 2. The first of these problems deals with cell-dependent micro-vibrations. The second one deals with propagation of the long waves related to micro-fluctuations of axial displacements. • The free micro-vibration frequencies have been determined, cf. Eqs. (31)-(33) and investigated. These frequencies depend on a periodicity cell size. Hence, they cannot be obtained in the framework of the asymptotic models commonly used for investigations of dynamics of periodic shells. The influence of the shell elastic, inertial and geometrical properties on the free micro-vibration frequencies has been analysed. From the numerical example, it follows that the free micro-vibration frequencies decrease with the decreasing of differences between inertial properties of the component materials, i.e. with the increasing of ratio ρ 2 /ρ 1 ∈ [0.01, 1], cf. Figs. 5b, 7b, 9b, but they increase with the decreasing of differences between elastic properties of the shell material components, i.e. with the increasing of ratio E 2 /E 1 ∈ [0.01, 1], cf. Figs. 5a, 7a, 9a. They also decrease with the decreasing of differences between the period length λ and the length dimension L 1 of the shell midsurface in periodicity direction, i.e. with the increasing of ratio λ/L 1 , cf. Figs. 4, 6, 8. • Some new important results have been obtained analysing the long wave propagation problem related to micro-fluctuations in axial direction. We deal with long waves if condition λ/L << 1 holds, where λ is the characteristic length dimension of the cell and L is the wavelength. It was shown that the tolerance-periodic heterogeneity of the shells leads to exponential waves and to dispersion effects, which cannot be analysed in the framework of the asymptotic models for periodic shells. Moreover, the new wave propagation speed (38) depending on the microstructure size has been obtained and investigated. The influence of the shell elastic, inertial and geometrical properties on this cell-dependent speed has been analysed. From the numerical example, it follows that the values of the wave propagation velocity increase with the decreasing of differences between elastic properties of the shell component materials, cf. Fig. 11, but they decrease with the decreasing of differences between inertial properties of the component materials, cf. Fig. 12. Values of the wave propagation speed decrease with the decreasing of differences between period length λ and the wavelength L ≡ L 2 , cf. Fig. 10. The strongest decrease in the speed takes place for ε ≡ λ/L ∈ [0.0001, 0.001].
Some others applications of the combined asymptotic-tolerance model will be shown in forthcoming papers.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.