Next Article in Journal
Evaluation of Corrosion Residual Life Prediction Methods for Metal Pipelines
Previous Article in Journal
Experimental Investigation of Hybrid Beams Utilizing Ultra-High Performance Concrete (UHPC) as Tension Reinforcement
Previous Article in Special Issue
Influence of Inertial Vlasov Foundation Parameters on the Dynamic Response of the Bernoulli—Euler Beam Subjected to A Group of Moving Forces—Analytical Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Effect of the Material Periodic Structure on Free Vibrations of Thin Plates with Different Boundary Conditions

by
Jarosław Jędrysiak
Department of Structural Mechanics, Łódź University of Technology, al. Politechniki 6, 90-924 Łódź, Poland
Materials 2022, 15(16), 5623; https://doi.org/10.3390/ma15165623
Submission received: 22 May 2022 / Revised: 25 July 2022 / Accepted: 4 August 2022 / Published: 16 August 2022

Abstract

:
Thin elastic periodic plates are considered in this paper. Since the plates have a microstructure, the effect of its size on behaviour of the plates can play a crucial role. To take into account this effect, the tolerance modelling method is applied. This method allows us to obtain model equations with constant coefficients, which involve terms dependent of the microstructure size. Using the model equations, not only can formulas of fundamental lower-order vibration frequencies be obtained, but also formulas of higher-order vibration frequencies related to the microstructure. In this paper, the effect of the material periodic microstructure on free vibration frequencies for various boundary conditions of the plates was analysed. To obtain proper formulas of frequencies, the Ritz method is applied. Moreover, some results are compared to the results calculated using the FEM.

1. Introduction

In this paper, thin elastic plates with a periodic structure in planes parallel to the plate midplane are the main object under consideration. Plates of this kind are made of many identical repeatable elements, which are called periodicity cells (marked by dotted lines in Figure 1). Problems of vibrations of these plates can be analysed using partial differential equations, having highly oscillating, periodic, and non-continuous functional coefficients. Because these equations are not good tools to consider special examples of these plates, various simplified approaches are proposed. These methods make it possible to obtain governing equations with constant coefficients. Between them, those based on the asymptotic homogenization should be mentioned, cf. [1]. However, the effect of the microstructure size is usually omitted in the model governing equations.
In order to model various mechanical problems of periodic (or microheterogeneous) structures and composites, other methods are also used. The microlocal parameters approach was applied to describe such plates, e.g., in [2]. In the paper [3], honeycomb sandwich composite shells were analysed using the two-scale asymptotic homogenization method. Euler–Bernoulli beams were considered in [4], where the relationship between the 3D and the homogenisation approach was presented. The stability of multi-cell thin-walled columns of rectangular cross-sections was analysed theoretically, numerically, and experimentally in the work [5]. A dynamic critical load for buckling of columns was considered in [6]. Global and local buckling of sandwich beams and plates were investigated in [7]. A finite element unit-cell method was used in [8] to analyse heterogeneous plates. Brito-Santana et al. [9] applied an asymptotic dispersive method to the problem of shear-wave propagation in a laminated composite. Analytical and numerical methods were used to analyse buckling problems for sandwich beams with variable properties of the core in papers [10,11]. For laminated composited plates, various approaches were applied, e.g., the strong formulation isogeometric analysis in the work [12]; a layer-wise theory combined with a differential quadrature method in [13]. A generalized differential quadrature method was also used for vibrations of nanocomposite functionally graded shells in a series of papers [14,15,16]. An analytical-numerical model based on analytical relations and finite element method was shown to analyse a torsion of auxetic composite beams with a cellular structure in [17].
Unfortunately, the proposed modelling methods for microstructured periodic media usually lead to the model equations without the effect of the microstructure size. On the other side, this effect can play a crucial role in dynamical problems of these media, cf. Brillouin [18], where there were observed macro- and microvibrations related to the macro- and microstructure, respectively. Some special methods were adopted to analyse some similar problems to describe this effect. For instance, a spectral element method was proposed to investigate vibration band gap of periodic Mindlin’s plates in [19]. The centre finite difference method was used in [20] to analyse free flexural vibrations of periodic stiffened thin plates. The simplified super element method was applied for these plates filled with viscoelastic damping material in [21]. Problems of wave bandgaps in periodic plate structures were considered using differential quadrature element in [22].
An alternative approach to analyse various mechanical problems of microstructured (periodic and non-periodic) media is the tolerance modelling method (or the tolerance method), cf. the books Woźniak and Wierzbicki [23], Woźniak et al. (eds.) [24]. This method is general and useful for modelling different problems, which are described by differential equations with highly oscillating, non-continuous functional coefficients. It leads from these exact governing equations to the averaged model equations, which have constant (or slowly-varying) coefficients, some of which explicitly are dependent of the size of the microstructure.
Various applications of this method for different periodic structures were shown in a series of papers. Dynamics of some plane periodic structures were analysed in [25]. Vibrations of one-directional periodic plates, having thickness smaller than the length of the periodic cell, were considered in [26]. Vibrations of periodic wavy-type plates were investigated in [27]. Dynamics of thin plates reinforced by a system of periodic stiffeners was presented in [28]. Vibrations of medium thickness periodic plates were considered in [29]. The buckling of periodic plates interacting with a periodic elastic foundation was investigated in [30]. Vibrations of thin periodic plates with the microstructure size of an order of the plate thickness were described in [31]. Stability of periodic shells was considered in papers [32,33]. Dynamics problems of medium thickness plates on a periodic foundation were described in [34]. Geometrically nonlinear thin periodic plates were considered in [35]. Vibration analysis of periodic geometrically nonlinear beams is shown in [36]. Modelling of vibrations of periodic three-layered plate-type structures were investigated in [37]. Revisiting version of the modelling for dynamics and stability for periodic slender visco-elastic beams on a foundation with damping was shown in [38]. A multi-scale analysis of stress distribution in thin composite periodic plates was considered in [39].
Moreover, the tolerance method was also successfully used to analyse structures with non-periodic micro-heterogeneity, which can be treated as being made of functionally graded materials in the macro-scale. For instance, dynamics for longitudinally graded plates was investigated in [40,41], but stability of similar plates in [42]. Modelling of thin-walled structures with dense system of ribs was shown in [43,44]. Heat conduction of composite cylindrical conductors with non-uniform distribution of constituents was analysed in [45,46]. Effect of microstructure in problems of thermoelasticity for transversally graded laminates was considered in [47]. Dynamics of thin functionally graded plates with one-directional microstructure was investigated in [48]. Free vibrations of medium thickness functionally graded plates were analysed in [49]. Dynamics of thin functionally graded microstructured shells was investigated in [50,51], but dynamics and stability of them in [52]. Stability problems of functionally graded microstructured thin plates on elastic foundation were considered in [53]. However, the above-mentioned works do not cover all the problems that were studied by these authors in the literature on the tolerance modelling method.
The main aim of this contribution is to apply the tolerance and asymptotic models of dynamic problems for thin elastic periodic plates, cf. [26,30], to free vibrations of periodic plate strips with various boundary conditions. The other aim is a certain analysis of the effect of the material periodic cell structure on free vibration frequencies, formulas of them are derived using the Ritz method. Some results are compared and justified by the finite element method.

2. Modelling Foundations

2.1. Preliminaries

Let 0x1x2x3 denote the orthogonal Cartesian co-ordinate system in the physical space and t be the time coordinate. Subscripts α, β, …(i, j, …) run over 1, 2 (1, 2, 3) and indices A, B, … (a, b, …) run over 1, …, N (1, …, n). Summation convention holds for all aforementioned indices. Denote also x ≡ (x1,x2) and zx3. Let the region Ω ≡ {(x,z): −h(x)/2 < z < h(x)/2, x ∈ Π} be the undeformed plate, with Π as the midplane, having along the x1- and x2-axis length dimensions L1 and L2, respectively; and h(x) as the plate thickness.
Plates under consideration are assumed to be periodic in planes parallel to the plate midplane, with periods l1 and l2 along the x1- and x2-axis directions, respectively. Let ∆ ≡ [−l1/2, l1/2] × [−l2/2, l2/2] denote the periodicity basic cell on 0x1x2 plane. The cell size can be described by a parameter l ≡ [(l1)2 + (l2)2]1/2, which satisfies the condition max(h) < l << min(L1,L2), and can be called the microstructure parameter. Denote by (∙), α ≡ ∂/∂xα the partial derivatives with respect to the space co-ordinate, and by the overdots the partial derivatives with respect to the time co-ordinate.
It is assumed that all plate properties can be periodic functions in x, e.g., thickness h(x), elastic moduli aijkl = aijkl(x,z), mass density ρ = ρ(x,z). Moreover, these plate material properties are even functions in z. By aαβγδ, aαβ33, a3333 denote the non-zero components of the elastic moduli tensor, and also introduce: cαβγδaαβγδaαβ33aγδ33(a3333)−1.
Let ui, eij and sij denote displacements, strains, and stresses, respectively; and u ¯ i and e ¯ i j —virtual displacements and virtual strains; and p—loadings (along the z-axis).
Considerations are based on the well-known assumptions of the Kirchhoff-type plate theory, but they are mentioned below.
  • The kinematic assumptions
u α ( x , z , t ) = z α w ( x , t ) , u 3 ( x , z , t ) = w ( x , t ) ,
where w(x,t) is the deflection of the midplane. Similar assumptions are made for virtual displacements:
u ¯ α ( x , z ) = z α w ¯ ( x ) , u ¯ 3 ( x , z ) = w ¯ ( x ) .
  • The strain-displacement relations
e α β = u ( α , β ) .
  • The stress-strain relations
(with the assumption that the plane of elastic symmetry is parallel to the plane z = 0)
s α β = c α β γ δ e γ δ ,
where:
c α β γ δ = a α β γ δ a α β 33 a 33 γ δ / a 3333 , c α 3 γ 3 = a α 3 γ 3 a α 333 a 33 γ 3 / a 3333 .
  • The virtual work equation
    Π h / 2 h / 2 ρ u ¨ i u ¯ i d z d a + Π h / 2 h / 2 s α β e ¯ α β d z d a = Π p u ¯ 3 ( x , h 2 ) d a   ,
Which is satisfied for arbitrary virtual displacements (2), under the assumption that these displacements are neglected on the boundary and sufficiently regular and independent functions; moreover: da = dx1dx2.
The plate properties, i.e., stiffness tensor dαβγδ, inertia properties: μ, j, are periodic functions in x, given by:
d α β γ δ ( x ) = h / 2 h / 2 c α β γ δ ( x , z ) z 2 d z , μ ( x ) = h / 2 h / 2 ρ ( x , z ) d z , j ( x ) = h / 2 h / 2 ρ ( x , z ) z 2 d z .
Combining the above assumptions (1)–(4) together with the virtual work Equation (6), after some manipulations using the divergence theorem and the du Bois-Reymond lemma to (6), the governing equation of thin elastic periodic plates takes the form:
( d α β γ δ w , γ δ ) , α β + μ w ¨ j w ¨ , α α = p .
Periodic plates coefficients of Equation (8) are usually discontinuous and highly oscillating, periodic functions in x. It is very difficult to find solutions to this equation.
Hence, in this note, the original equation is replaced by a system of equations with constant coefficients, which can describe (or not) the information about the microstructure using the tolerance averaging method.

2.2. Foundations of the Tolerance Modelling

2.2.1. Introductory Concepts

Some introductory concepts are applied in the tolerance modelling method. These concepts can be found in the book [24] and they are formulated for beams in [38] or for non-periodic plates in [49,53]. Hence, here they can be only mentioned: the averaging operator <·>, the tolerance parameter, the tolerance-periodic function TP(Π,Δ), the slowly-varying function SV(Π,Δ), the highly oscillating function HO(Π,Δ), and the fluctuation shape function FS(Π,Δ). Below, let us restate the averaging operator.
Let Δ(x) = x + Δ denote a cell at x ∈ ΠΔ, ΠΔ = {x ∈ Π: Δ(x) ⊂ Π}. The averaging operator, one of the fundamental concepts of the tolerance method, is defined by
< f > ( x ) = ( l 1 l 2 ) 1 Δ ( x ) f ( y 1 , y 2 ) d y 1 d y 2 , x Π Δ , y Δ ( x ) ,
for an integrable function f. For function f, being periodic in x its averaged value from (9) is constant.
A very important concept is a fluctuation shape function. It is a function, g(·), being continuous together with gradient 1g and with a piecewise continuous and bounded gradient ∂2g; depending on l as a parameter, which satisfies proper conditions, i.e., g F S ( Π , Δ ) if conditions hold:
( i ) k g O ( l α k )   for   k = 0 , 1 , , α ,   α = 2 ,   0 g g , ( i i ) < g > ( x ) 0   x Π Δ ,
where l is the microstructure parameter. Condition (10ii) can be replaced by <μg>(x) ≈ 0 for every x ∈ ΠΔ, where μ > 0 is a certain tolerance-periodic function.

2.2.2. Tolerance Modelling Assumptions

The fundamental tolerance modelling assumptions in their general form are formulated in the book [24]. For thin plates with a microstructure, they can be found, e.g., in [26,40,49]. Here, these assumptions are written below in the form for thin periodic plates.
One of these assumptions is the micro-macro decomposition, in which the deflection is assumed that can be decomposed in the form:
w ( x , t ) = W ( x , t ) + g A ( x ) V A ( x , t ) , A = 1 , , N ,
with the basic unknowns: W(·,t) called the macrodeflection, VA(·,t) called the fluctuation amplitudes, W ( , t ) , V A ( , t ) S V ( Π , Δ ) ; and the known fluctuation shape functions g A ( ) F S ( Π , Δ ) . Fluctuation shape functions can be derived as solutions to eigenvalue problems posed on the periodicity cell. However, in most cases they can be assumed in an approximate form as trigonometric functions [26,53] or saw-type functions [49].
Similar assumptions to (11) should be introduced for virtual displacements w ¯ ( ) :
w ¯ ( x ) = W ¯ ( x ) + g A ( x ) V ¯ A ( x ) , A = 1 , , N ,
with slowly-varying functions W ¯ ( ) , V ¯ A ( ) S V ( Π , Δ ) .
The second modelling assumption is the tolerance averaging approximation, such that terms O(δ) are negligibly small in the course of modelling, and they can be neglected in the following formulas:
( i ) < ϕ > ( x ) = < ϕ ¯ > ( x ) + O ( δ ) , ( i i ) < ϕ F > ( x ) = < ϕ > ( x ) F ( x ) + O ( δ ) , ( i i i ) < ϕ ( g F ) , γ > ( x ) = < ϕ g , γ > ( x ) F ( x ) + O ( δ ) , x Π ;       γ = 1 , α ; α = 1 , 2 ; 0 < δ < < 1 ; ϕ T P ( Π , Δ ) ,   F S V ( Π , Δ ) ,   g F S ( Π , Δ ) .

2.2.3. Outline of the Modelling Procedure

In the modelling procedure, the above concepts and basic assumptions are applied (11)–(13). The procedure can be divided into a few steps.
The first step stands a substitution of micro-macro decompositions (11) and (12) into the virtual work Equation (6). In the next step the resulting equation is averaged using the averaging operation (9) over a periodicity cell, cf. [26]. In the third step, the tolerance averaged virtual work equation is obtained after using Formula (13) of the tolerance averaging approximation. Introducing the following denotations of averaged parameters, which can be stand averaged constitutive relations:
M α β < h / 2 h / 2 s α β z d z > , M A < g , α β A h / 2 h / 2 s α β z d z > ,
the tolerance averaged virtual work equation takes the form:
Π ( < μ > W ¨ + < μ g B > V ¨ B ) δ W d a + Π ( < μ g A > W ¨ + < μ g A g B > V ¨ B ) δ V A d a Π ( < j > W ¨ , α α + < j g , α B > V ¨ , α B ) δ W d a + Π ( < j g , α A > W ¨ , α + < j g , α A g , α B > V ¨ B ) δ V A d a + + Π M α β , α β δ W d a + Π M A δ V A d a = Π p δ W d a .
At the end, after some manipulations governing equations of the tolerance model are obtained using the divergence theorem and the du Bois-Reymond lemma to Equation (15).

2.3. Governing Equations

2.3.1. Tolerance Model Equations

Let us denote:
D α β γ δ < d α β γ δ > , D α β A < d α β γ δ g , γ δ A > , D A B < d α β γ δ g , α β A g , γ δ B > , m < μ > , m A l 2 < μ g A > , m A B l 4 < μ g A g B > , ϑ < j > , ϑ α A l 1 < j g , α A > , ϑ α β A B l 2 < j g , α A g , β B > , P < p > , P A l 2 < p g A > .  
The tolerance modelling procedure leads to the system of equations for the macrodeflection W and the fluctuation amplitudes of the deflection VA:
( D α β γ δ W , γ δ + D α β A V A ) , α β + m W ¨ + l 2 m A V ¨ A ϑ W ¨ , α α l ϑ α A V ¨ , α A = P , D α β A W , γ δ + D A B V B + l 2 m A W ¨ + l ϑ α A W ¨ , α + l 2 ( l 2 m A B + ϑ α β A B ) V ¨ B = l 2 P A .
The above Equation (17) has constant coefficients and, together with micro-macro decompositions (11) and (12), constitutes the tolerance model of thin elastic periodic plates. This model allows us to describe the effect of the microstructure size on the plate dynamics by terms with the microstructure parameter l. It can be observed that the basic unknowns of Equation (17) have to be slowly varying functions in x, W ( , t ) , V A ( , t ) S V ( Π , Δ ) . Moreover, there have to be determined boundary conditions for the macrodeflection W.

2.3.2. Asymptotic Model Equations

In order to compare results, the asymptotic model equations can be obtained. These equations can be derived using the formal asymptotic modelling procedure. However, below this is made by simple neglecting terms of an order of O(ln), n = 1,2,…, in Equation (17).
Hence, after these manipulations, Equation (17) is simplified to the form of the asymptotic model equations:
( D α β γ δ W , γ δ + D α β A V A ) , α β + m W ¨ ϑ W ¨ , α α = P , D α β A W , γ δ + D A B V B = 0 ,
with all constant coefficients.
Equation (18), together with micro-macro decompositions (11) and (12), constitutes the asymptotic model of thin elastic periodic plates. Similar boundary conditions should be formulated as in the tolerance model, i.e., only for the macrodeflection W. Moreover, it can be observed that Equation (18)2 stand a system of linear algebraic equations for the fluctuation amplitudes VA, A = 1, …, N.

3. Analysis of Free Vibrations of Periodic Plate Strips with Various Boundary Conditions

3.1. Introduction to the Example

Let us consider free vibrations of a thin periodic plate strip having span L along the x1-axis. The loading p is neglected, p = 0. The plate strip has periodic structure along its span, cf. Figure 2. However, all plate properties are independent of the x2-coordinate.
Hence, our considerations can be treated as independent of the x2-coordinate. Let us denote x = x1; z = x3; and ∂ ≡ (·),1; x ∈ [0,L]; z ∈ [−h/2,h/2], with a constant plate thickness h. The time derivatives are still denoted by overdots ( ) . Let Δ [ l / 2 , l / 2 ] be the basic cell in the interval ∧ ≡ [0, L], and l, l << L, be the length of this cell, satisfying the condition l << L. By Δ ( x ) [ x l / 2 , x + l / 2 ] denote a cell with a centre at x ∈ [0, L].
It is assumed that the plate strip is made of two elastic isotropic materials, described by: Young’s moduli E′, E′′, Poisson’s ratios ν′, ν′′ and mass densities ρ′, ρ′′, respectively. These materials are perfectly bonded on interfaces and periodically distributed along the x-axis.
It is assumed that the properties of the plate strip are described by the functions:
E ( y ) = E , for y ( ( 1 γ ) l / 2 , ( 1 + γ ) l / 2 ) , E , for y [ 0 , ( 1 γ ) l / 2 ] [ ( 1 + γ ) l / 2 , l ] ,
ρ ( y ) = ρ , for y ( ( 1 γ ) l / 2 , ( 1 + γ ) l / 2 ) , ρ , for y [ 0 , ( 1 γ ) l / 2 ] [ ( 1 + γ ) l / 2 , l ] ,
ν ( y ) = ν , for y ( ( 1 γ ) l / 2 , ( 1 + γ ) l / 2 ) , ν , for y [ 0 , ( 1 γ ) l / 2 ] [ ( 1 + γ ) l / 2 , l ] ,
with a distribution parameter of material properties γ, cf. Figure 3. Under the above assumptions, the effect of material periodic structure on free vibration frequencies will be considered with using ratios: E″/E′ ∈ [0, 1], ν″/ν′ ∈ [0, 1], ρ″/ρ′ ∈ [0, 1], h/l ∈ (0, 0.1], γ ∈ [0, 1].
Effects of the above parameters will be shown only for the first vibration—lower and higher (for the tolerance model). Hence, the considerations can be restricted to assume only one fluctuation shape function, i.e., A = N = 1. The analysis for these functions can be found in [24]. Denote g ≡ g1, V ≡ V1, ≡ (∙),1. Micro-macro decomposition (11) of the deflection w(x,t) of the plate strip is given by:
w ( x , t ) = W ( x , t ) + g ( x ) V ( x , t ) ,
where W ( , t ) ,   V ( , t ) S V ( Λ , Δ ) for every t ( t 0 , t 1 ) , g ( ) F S ( Λ , Δ ) .
Because the cell is assumed in the form shown in Figure 3, i.e., it has a symmetry axis normal to the midplane, the fluctuation shape function g(x) is assumed as
g ( y ) = l 2 [ cos ( 2 π y / l ) + c ] , y Δ ( x ) , x Λ ,
where the constant c is determined by the condition < μ g > = 0 and has the form:
c = sin ( π γ ) ( ρ ρ ) π [ ρ γ + ρ ( 1 γ ) ]
with γ as the constant distribution parameter of material properties.
Using definitions of periodic Function (7) for the considered plate strip, it can be written:
d ( x ) h 3 12 [ 1 ν ( x ) 2 ] E ( x ) , μ ( x ) h ρ ( x ) , j ( x ) h 3 12 ρ ( x ) .
It can be observed that for the assumed periodic cell, Figure 3, and the assumed fluctuation shape Function (23), term (16)8 neglects:
< j g > = 0 .
Denoting:
D = < d > , D = < d g > , D ¯ = < d g g > , μ = < μ > , μ ¯ = l 4 < μ g g > , ϑ = < ϑ > , ϑ ¯ = l 2 < ϑ g g > ,
tolerance model Equation (17) has the form:
( D W + D V ) + μ W ¨ ϑ W ¨ = 0 , D W + D ¯ V + l 2 ( l 2 μ ¯ + ϑ ¯ ) V ¨ = 0 .
Equation (27) describe free vibrations of the periodic plate strip within the tolerance model.
Using denotations (26) for the plate strip under consideration, Equation (18)1 takes the form
[ ( D D 2 / D ¯ ) W ] + μ W ¨ ϑ W ¨ = 0 .
The above equation describes free vibrations of this plate strip in the framework of the asymptotic model.

3.2. The Using of the Ritz Method

In order to solve the problem of free vibrations and to find free vibration frequencies for different boundary conditions, the known Ritz method is applied, cf. [49,53]. In this method, formulas of the maximal strain energy Y max and the maximal kinetic energy K max should be determined.
The solutions to Equations (28) and (27) for the plate strip under consideration can be assumed as:
W ( x , t ) = A W Ψ ( α x ) cos ( ω t ) , V ( x , t ) = A V Θ ( α x ) cos ( ω t ) ,
where α is a wave number, ω is a free vibration frequency, AW and AV are amplitudes. Functions Ψ(·) and Θ(·) are eigenvalue functions for the macrodeflection and the fluctuation amplitude, respectively, which have to satisfy the given boundary conditions for x = 0, L. To simplify, let us denote the first and second order derivatives of functions Ψ(·) and Θ(·) by:
Ψ ( α x ) α Ψ ˜ ( α x ) , Θ ( α x ) α Θ ˜ ( α x ) , Ψ ( α x ) α 2 Ψ ¯ ( α x ) , Θ ( α x ) α 2 Θ ¯ ( α x ) .
In the following analysis four different cases of boundary conditions of the plate strip will be considered, i.e.:
  • the simply supported plate strip
    Ψ ( 0 ) = Ψ ( 0 ) = Ψ ( L ) = Ψ ( L ) = 0 ;
  • the clamped-hinged plate strip
    Ψ ( 0 ) = Ψ ( 0 ) = Ψ ( L ) = Ψ ( L ) = 0 ;
  • the plate strip clamped on both edges
    Ψ ( 0 ) = Ψ ( 0 ) = Ψ ( L ) = Ψ ( L ) = 0 ;
  • the cantilever plate strip
    Ψ ( 0 ) = Ψ ( 0 ) = Ψ ( L ) = Ψ ( L ) = 0 .
The eigenvalue functions Ψ(·) and Θ(·) from solutions (29) can be assumed as for the homogeneous plate strip as functions satisfying the boundary conditions (31)–(34). In order to describe these solutions, the Krylov–Prager functions are introduced:
S ( α x ) = 1 2 [ cosh ( α x ) + cos ( α x ) ] , T ( α x ) = 1 2 [ sinh ( α x ) + sin ( α x ) ] , W ( α x ) = 1 2 [ cosh ( α x ) cos ( α x ) ] , V ( α x ) = 1 2 [ sinh ( α x ) sin ( α x ) ] .
Hence, the aforementioned eigenvalue functions Ψ(·) and Θ(·) can be written as:
  • the simply supported plate strip
    Ψ ( α x ) = sin ( α x ) , Θ ( α x ) = sin ( α x ) ;
  • the clamped-hinged plate strip
    Ψ ( α x ) = Θ ( α x ) = W ( α x ) coth ( α L ) V ( α x ) ;
  • the plate strip clamped on both edges
    Ψ ( α x ) = Θ ( α x ) = W ( α x ) cosh ( α L ) cos ( α L ) sinh ( α L ) sin ( α L ) V ( α x ) ;
  • the cantilever plate strip
    Ψ ( α x ) = Θ ( α x ) = W ( α x ) sinh ( α L ) sin ( α L ) cosh ( α L ) + cos ( α L ) V ( α x ) .
Now, applying the Ritz method formulas for the maximal strain energy Y max and the maximal kinetic energy K max should be formulated in the framework of both models—the tolerance and the asymptotic. Then, the conditions of the Ritz method have to be used:
( Y max K max ) A W = 0 , ( Y max K max ) A V = 0 ,
from which formulas of free vibrations frequencies can be determined.
Denotations (26), after substituting into them the fluctuation shape function (23), functions of material properties (19)–(21) and eigenvalue functions Ψ(·) and Θ(·), take the form:
D = h 3 12 [ E 1 ν 2 ( 1 γ ) + γ E 1 ν 2 ] 0 L [ Ψ ¯ ( α x ) ] 2 d x , D = π h 3 3 ( E 1 ν 2 E 1 ν 2 ) sin ( π γ ) 0 L Ψ ¯ ( α x ) Θ ( α x ) d x , D ¯ = ( π h ) 3 3 { ( E 1 ν 2 E 1 ν 2 ) [ 2 π γ + sin ( 2 π γ ) ] + 2 π E 1 ν 2 } 0 L [ Θ ( α x ) ] 2 d x , μ = h [ ( 1 γ ) ρ + γ ρ ] 0 L [ Ψ ( α x ) ] 2 d x , ϑ = h 3 12 [ ( 1 γ ) ρ + γ ρ ] 0 L [ Ψ ˜ ( α x ) ] 2 d x , μ ¯ = h 4 π { ( ρ ρ ) [ 2 π γ + sin ( 2 π γ ) ] + 2 π ρ } 0 L [ Θ ( α x ) ] 2 d x + + h π ( ρ ρ ) c [ π c γ 2 sin ( π γ ) ] 0 L [ Θ ( α x ) ] 2 d x + h ρ c 2 0 L [ Θ ( α x ) ] 2 d x , ϑ ¯ = π h 3 12 { ( ρ ρ ) [ 2 π γ sin ( 2 π γ ) ] + 2 π ρ } 0 L [ Θ ( α x ) ] 2 d x .
The maximal strain energy Y max and the maximal kinetic energy K max by the tolerance model can be written as:
Y max T M = 1 2 [ ( D A W 2 α 2 + 2 D A W A V ) α 2 + D ¯ A V 2 ] , K max T M = 1 2 [ A W 2 ( μ + ϑ α 2 ) + A V 2 l 2 ( l 2 μ ¯ + ϑ ¯ ) ] ω 2 .
However, in the framework of the asymptotic model, these formulas have the form:
Y max A M = 1 2 [ ( D A W 2 α 2 + 2 D A W A V ) α 2 + D ¯ A V 2 ] , K max A M = 1 2 A W 2 ( μ + ϑ α 2 ) ω 2 .
Using conditions (40) to relations (42), after some manipulations, the following formulas are obtained:
( ω , + ) 2 l 2 ( l 2 μ ¯ + ϑ ¯ ) α 4 D + ( μ + α 2 ϑ ) D ¯ 2 ( μ + α 2 ϑ ) l 2 ( l 2 μ ¯ + ϑ ¯ ) [ l 2 ( l 2 μ ¯ + ϑ ¯ ) α 4 D ( μ + α 2 ϑ ) D ¯ ] 2 + 4 ( α 2 D ) 2 l 2 ( μ + α 2 ϑ ) ( l 2 μ ¯ + ϑ ¯ ) 2 ( μ + α 2 ϑ ) l 2 ( l 2 μ ¯ + ϑ ¯ ) ,
of the lower  ω and the higher ω +  free vibration frequencies, respectively, in the framework of the tolerance model.
On the other side, the conditions (40) used to relations (43) of the asymptotic model led to the following formula:
ω 2 = α 4 D D ¯ D 2 D ¯ ( μ + ϑ α 2 )
of the lower free vibration frequency  ω .
Under the above analytical results, it can be observed that in the framework of the tolerance model, the effect of the microstructure size of the plate strip can be analysed in the form of higher free vibration frequencies, (44)2, but the asymptotic model leads only to the fundamental lower frequency, (45).

3.3. Results of Calculations

Let us define dimensionless frequency parameters:
Ω 2 ρ E L 2 ω 2 , Ω + 2 ρ E L 2 ω + 2 , Ω 2 ρ E L 2 ω 2 ,
where the free vibration frequencies ω, ω+, and ω are determined by Equations (44) and (45), respectively.
Plots in Figure 4a, Figure 7a, Figure 10a, Figure 14a and Figure 17a are made for the simply supported plate strip, in Figure 4b, Figure 7b, Figure 10b, Figure 14b and Figure 17b for the clamped-hinged plate strip, in Figure 5a, Figure 8a, Figure 11a, Figure 15a and Figure 18a for the plate strip clamped on both edges, and in Figure 5b, Figure 8b, Figure 11b, Figure 15b and Figure 18b for the cantilever plate strip. Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8, Figure 10, Figure 11 and Figure 12, Figure 14, Figure 15, Figure 17 and Figure 18 show results of the lower frequency parameters (46)1,3, but Figure 9, Figure 13, Figure 16 and Figure 19 are the higher frequency parameters (46)2.
The results shown in Figure 4 and Figure 5 and in Figure 9a are in the form of surfaces of the frequency parameters versus pairs of ratios (ρ″/ρ′; E″/E′). In calculations there are assumed: the Poisson’s ratios ν′′ = ν′ = 0.3, ratio l/L = 0.1, ratio h/l = 0.1 and γ = 0.3, 0.5, 0.8.
The lower frequency parameters of the homogeneous plate strips are presented in Figure 4, Figure 5, Figure 10, Figure 11, Figure 14, Figure 15, Figure 17 and Figure 18 by grey planes (γ = 1) and they are related to values Ω: 0.0298656 for the simply supported plate strip (Figure 4a, Figure 10a, Figure 14a and Figure 17a), 0.0466553 for the clamped-hinged plate strip (Figure 4b, Figure 10b, Figure 14b and Figure 17b), 0.0676988 for the clamped-clamped plate strip (Figure 5a, Figure 11a, Figure 15a and Figure 18a), and 0.0106397 for the cantilever plate strip (Figure 5b, Figure 11b, Figure 15b and Figure 18b).
The effect of differences of Young’s moduli (ratios E″/E′) and mass densities (ratios ρ″/ρ′) on the lower frequency parameters can be observed in Figure 4 and Figure 5, i.e., they increase with the increasing of ratio E″/E′; they decrease with the increasing of ratio ρ″/ρ′.
There are pairs of ratios ((ρ″/ρ′)*; (E″/E′)*) such that, cf. Figure 4 and Figure 5: the lower frequency parameters are smaller than this parameter for the homogeneous plate strip (γ = 1) for E″/E′ ≤ (E″/E′)* and ρ″/ρ′ ≥ (ρ″/ρ′)*; they are bigger than this parameter for the homogeneous plate strip (γ = 1) for E″/E′ > (E″/E′)* and ρ″/ρ′ < (ρ″/ρ′)*.
Figure 6 shows traces of surfaces of the lower frequency parameters (46)1,3, from Figure 4 and Figure 5 (pairs of ratios ((ρ″/ρ′)*; (E″/E′)*) as curves), on planes of these frequency parameters for the homogeneous plate strips (γ = 1) (for ν″/ν′ = 1, h/l = 0.1, l/L = 0.1) for all considered boundary conditions, mentioned above. It can be observed that these traces are the same for various boundary conditions of periodic plate strips under consideration, for fixed values of the parameter γ = 0.3, 0.5, 0.8.
However, Figure 7, Figure 8 and Figure 9b present surfaces of dimensionless frequency parameters versus pairs of ratios (γ; h/L), γ ∈ [0, 1], h/l ∈ (0, 0.1]. These calculations are made for: the Poisson’s ratios ν″ = ν′ = 0.3, ratio E″/E′ = 0.5, ratios ρ″/ρ′ = 0.2, 0.5, 0.8.
From Figure 7 and Figure 8, it can be observed that the lower frequency parameters increase with the increasing of h/L and decrease with the increasing of γ for certain values of the ratio ρ″/ρ′, e.g., ρ″/ρ′ = 0.2, but for other values of this ratio they are more varied (it is the parameter γ0 that the frequency parameters decrease for γ ≤ γ0 but increase for γ > γ0). All surfaces of the lower frequency parameters have the same edge for γ = 1.
Similar effects of the ratios ρ″/ρ′, E″/E′ as for the lower frequency parameters presented in Figure 4 and Figure 5 can also be found for the higher frequency parameters in Figure 9a, i.e., they increase with the increasing of ratio E″/E′; they decrease with the increasing of ratio ρ″/ρ′.
The dependency of the higher frequency parameters on the distribution parameter of material properties γ and the ratio of the thickness h/L, can be found in Figure 9b (for fixed ratios: E″/E′ = 0.5; ρ″/ρ′ = 0.2, 0.5, 0.8). It can be observed that the higher frequency parameters increase with the increasing of h/L; the dependency of these parameters on the parameter γ is more complicated, since for fixed ratios E″/E′ = 0.5 and ρ″/ρ′ = 0.2, 0.5, 0.8 there are γ1 and γ2 such that these frequency parameters increase for γ ≤ γ1 or γ > γ2, but decrease for γ1 < γ ≤ γ2. Similarly to lower frequency parameters, cf. Figure 7 and Figure 8, all surfaces of the higher frequency parameters have the same edge for γ = 1.
Surfaces of dimensionless frequency parameters versus pairs of ratios (γ; E″/E′), γ ∈ [0, 1], E″/E′ ∈ [0, 1], are shown in Figure 10, Figure 11 and Figure 13. These calculations are made for: the Poisson’s ratios ν″ = ν′ = 0.3, ratio h/l = 0.1, ratios ρ″/ρ′ = 0.3, 0.5, 0.7.
From surfaces of the lower frequency parameters versus pairs (γ; E″/E′) presented in Figure 10 and Figure 11, it can be observed that the frequency parameters increase with the increasing of E″/E′ and with the increasing of γ. Moreover, all surfaces of the lower frequency parameters have the same edge for γ = 1.
It can be distinguished that for pairs of ratios (γ*; (E″/E′)*) such as cf. Figure 10 and Figure 11, the lower frequency parameters are smaller than this parameter for the homogeneous plate strip (γ = 1) for some values γ* and E″/E′ ≤ (E″/E′)*(γ*), where (E″/E′)* is dependent on γ*; they are bigger than this parameter for the homogeneous plate strip (γ = 1) for E″/E′ > (E″/E′)*(γ*), where (E″/E′)* is dependent on γ*.
In Figure 12, traces of surfaces of the lower frequency parameters are presented (46)1,3, from Figure 10 and Figure 11, on planes related to the homogeneous plate strips (γ = 1) (for ν″/ν′ = 1, h/l = 0.1, l/L = 0.1) for all considered boundary conditions (the pairs of ratios (γ*; (E″/E′)*) from Figure 10 and Figure 11). Similarly to Figure 6, these traces are the same for various boundary conditions of these periodic plate strips, for fixed values of the ratio ρ″/ρ′ = 0.3, 0.5, 0.7.
The higher frequency parameters depend on the distribution parameter of material properties γ and the ratio E″/E′, cf. Figure 13 (for fixed ratios: ν″/ν′ = 1; ρ″/ρ′ = 0.3, 0.5, 0.7). From this figure it can be observed that the frequency parameters increase with the increasing of E″/E′. However, the dependency of these parameters on the parameter γ is more complicated, because for the fixed ratios ν″/ν′ = 1 and ρ″/ρ′ = 0.3, 0.5, 0.7, there are γ1 and γ2 such that these frequency parameters increase for γ ≤ γ1 or γ > γ2 but decrease for γ1 < γ ≤ γ2. Similarly to the lower frequency parameters, all surfaces of the higher frequency parameters have the same edge for γ = 1.
Surfaces of dimensionless frequency parameters versus pairs of ratios (γ; ν″/ν′), γ ∈ [0, 1], ν″/ν′ ∈ [0, 1], are presented in Figure 14 and Figure 15 (for lower frequencies) and in Figure 16 (for higher frequencies). These plots are made for ratios: E″/E′ = 0.5, h/l = 0.1, ρ″/ρ′ = 0.3, 0.5, 0.7.
From surfaces of the lower frequency parameters versus pairs (γ; ν″/ν′) presented in Figure 14 and Figure 15, it can be observed that:
  • the frequency parameters increase with the increasing of ν″/ν′;
  • the dependency of these parameters on the parameter γ is more complicated, since (for the fixed ratio E″/E′ = 0.5):
    -
    for some fixed values of ρ″/ρ′, e.g., ρ″/ρ′ = 0.7, these frequency parameters increase with the increasing of γ,
    -
    for some fixed values of ρ″/ρ′, e.g., ρ″/ρ′ = 0.5, there is γ0 such that these frequency parameters decrease for γ ≤ γ0, but increase for γ0 < γ,
    -
    for some fixed values of ρ″/ρ′, e.g., ρ″/ρ′ = 0.3, there are γ1 and γ2 such that these frequency parameters decrease for γ ≤ γ1 or γ > γ2, but increase for γ1 < γ ≤ γ2;
  • values presented in surfaces of the lower frequency parameters are smaller than this parameter for the homogeneous plate strip (γ = 1) for—the fixed ratio E″/E′ = 0.5 and some fixed ratios ρ″/ρ′, ρ″/ρ′ > 0.3;
  • for some fixed ratios ρ″/ρ′, ρ″/ρ′ ≤ 0.3, there are some values γ* and ν″/ν′ ≥ (ν″/ν′)*(γ*), where (ν″/ν′)* is dependent on γ*, such that the lower frequency parameters are bigger than this parameter for the homogeneous plate strip (γ = 1).
  • Moreover, all surfaces of the lower frequency parameters have the same edge for γ = 1.
Figure 16 shows that the higher frequency parameters increase with the increasing of ν″/ν′. Moreover, the dependency of these parameters on the parameter γ is more complicated, since (for the fixed ratio E″/E′ = 0.5):
-
for some fixed ratios ρ″/ρ′ = 0.5, 0.7 there are γ1 and γ2 such that these frequency parameters increase for γ ≤ γ1 or γ > γ2, but decrease for γ1 < γ ≤ γ2,
-
however, this dependency of these parameters on the parameter γ for the rather small fixed ratio ρ″/ρ′, e.g., ρ″/ρ′ = 0.3, has other form, because there are γ3 and γ4 such that these frequency parameters decrease for γ ≤ γ3 or γ > γ4, but increase for γ3 < γ ≤ γ4.
Similarly to the lower frequency parameters, cf. Figure 14 and Figure 15, all surfaces of the higher frequency parameters have the same edge for γ = 1.
Similar graphs, but versus pairs of ratios (γ; ρ″/ρ′), γ ∈ [0, 1], ρ″/ρ′ ∈ [0, 1], can be found in Figure 17 and Figure 18 (for lower frequencies) and in Figure 19 (for higher frequencies). These results are obtained for parameters: ν″/ν′ = 1, ρ″/ρ′ = 0.5, h/l = 0.1, E″/E′ = 0.3, 0.5, 0.7.
From Figure 17 and Figure 18, it can be concluded that:
  • the frequency parameters decrease with the increasing of ρ″/ρ′;
  • the dependency of these parameters on the parameter γ is more complicated, because (for the fixed ratio ν″/ν′ = 1):
    -
    for some fixed ratios E″/E′, e.g., E″/E′ = 0.5, 0.7, the frequency parameters decrease with the increasing of the parameter γ,
    -
    for some fixed ratios E″/E′, E″/E′ < 0.5, it is γ0 such that these frequency parameters decrease for γ ≤ γ0, but increase for γ0 < γ;
  • the lower frequency parameters are smaller than this parameter for the homogeneous plate strip (γ = 1) for some values γ* and ρ″/ρ′ ≥ (ρ″/ρ′)*(γ*), where (ρ″/ρ′)* is dependent on γ*;
  • the frequency parameters are bigger than this parameter for the homogeneous plate strip (γ = 1) for some γ* and ρ″/ρ′ < (ρ″/ρ′)*(γ*), where (ρ″/ρ′)* depends on γ*.
Moreover, all surfaces of the lower frequency parameters have the same edge for γ = 1.
Figure 19 presents the higher frequency parameters increase with the decreasing of ρ″/ρ′. On the other hand, the dependency of these parameters on the parameter γ is more complicated, since (for the fixed ratio ν″/ν′ = 1):
-
for more fixed ratios E″/E′ < 0.7 (e.g., E″/E′ = 0.3, 0.5) it is γ0 (dependent on ρ″/ρ′) such that these frequency parameters increase for γ ≤ γ0, but decrease for γ0 < γ,
-
however, this dependency of these parameters on the parameter γ for the rather big fixed ratio E″/E′, e.g., E″/E′ = 0.7, has other form, because there are γ1 and γ2 such that these frequency parameters decrease for γ ≤ γ1 or γ > γ2, but increase for γ1 < γ ≤ γ2.
Similarly to the lower frequency parameters, cf. Figure 17 and Figure 18, all surfaces of the higher frequency parameters have the same edge for γ = 1.
Moreover, all calculations are made for the wave number α related to the first form of eigenfunction for every case of the supports, i.e.,: α = π for (31), Figure 4a, Figure 7a, Figure 10a, Figure 14a and Figure 17a; α = 3.9266 for (32), Figure 4b, Figure 7b, Figure 10b, Figure 14b and Figure 17b; α = 4.7300 for (33), Figure 5a, Figure 8a, Figure 11a, Figure 15a and Figure 18a; and α = 1.8751 for (34), Figure 5b, Figure 8b, Figure 11b, Figure 15b and Figure 18b.
It can be observed that the tolerance and the asymptotic models lead to nearly identical values of lower frequency parameters, cf. Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 10, Figure 11, Figure 12, Figure 14, Figure 15, Figure 17 and Figure 18. Moreover, plots of surfaces of higher frequency parameters are the same for various boundary conditions of periodic plate strips under consideration.

4. Some Final Remarks

Vibrations of thin periodic plates for various boundary conditions have been modelled here. Using the tolerance modelling, the known differential equation of Kirchhoff-type plates has been averaged. This method leads from the governing differential equation with non-continuous, periodic coefficients to the system of differential equations with constant coefficients. Using this method, the effect of the microstructure size on the overall dynamic behaviour of the plates under consideration is taken into account in the derived tolerance model equations. On the other hand, a simplified averaged model—the asymptotic model—is described by the governing equation without this effect.
Under the presented considerations, the following general remarks, common to the application of the tolerance modelling method for different microheterogeneous objects, can be formulated.
  • The tolerance model allows us to analyse the effect of the microstructure size on dynamic problems of thin periodic plates under consideration, e.g., the “higher order” vibrations, related to the plate microstructure.
  • A certain a posteriori criterion of physical reliability for the model is that the basic unknowns W, VA, A = 1, …, N, have to be slowly-varying functions. Moreover, under these conditions, the governing equations of the tolerance model have a physical sense.
  • Using the asymptotic model of periodic plates, lower order (fundamental) vibrations can be only analysed.
Summarising the remarks of the results of calculations from Section 3, some additional, more general comments can be formulated:
  • Lower free vibration frequencies (also called fundamental frequencies) can be analysed within both the models—the tolerance and the asymptotic.
  • Values of lower free vibration frequencies of the considered periodic plate strips also depend on boundary conditions of these strips and they change from the highest values to the smallest as for the homogeneous plate strips with identical boundary conditions: the clamped-clamped plate strip, the clamped-hinged plate strip, the simply supported plate strip, the cantilever plate strip.
  • Higher free vibration frequencies, related to the periodic microstructure of the plate, can be investigated only within the tolerance model.
  • Values of higher free vibration frequencies of the considered periodic plate strips do not depend on boundary conditions of these plates, but only on material and geometrical properties of the plates.
  • The effects of differences between material or geometrical parameters, such that Young’s moduli (the ratio E″/E′), mass densities (the ratio ρ″/ρ′), Poisson’s ratios (the ratio ν″/ν′), the plate thickness (the ratio h/L), on free vibration frequencies are similar for both kinds of them—lower and higher.
  • The effect of the distribution parameter of material properties γ on these frequencies is more complicated to describe and is slightly different for both kinds of frequencies—lower and higher.
  • The effect of the parameter γ on the frequencies is related to the material properties Young’s moduli (the ratio E″/E′), mass densities (the ratio ρ″/ρ′), and Poisson’s ratios (the ratio ν″/ν′).
  • Sections by the planes corresponding to the fundamental free vibration frequency of homogeneous plate strips of surfaces for the lower frequencies of periodic plate strips are identical for all boundary conditions. Therefore, created in this way, traces of these surfaces of lower frequencies on the relevant planes of frequencies are the same for all considered boundary conditions.
In this paper, an application of the tolerance model to the analysis free vibration frequencies—lower and higher—are shown for periodic plate strips with various boundary conditions. Selected results were also compared with the results obtained by the finite element method, cf. Appendix A.
We presented a wide analysis of issues of free vibrations of periodic plates, which allows us to notice that the tolerance model seems to be a good tool in the study of this kind of problem. Moreover, this model can be used in the consideration of optimisation of periodic plates.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

Here, obtained results for some selected cases of the simply supported plate strips are compared with results calculated using a commercial computer program of the finite element method (FEM)—Abaqus, with a type of elements S8R5. These results are shown in Table A1. The comparisons are made only for the fundamental lower frequency parameters, because they can only be calculated in the framework of the commercial computer programs of FEM. The parameter of differences between the results of the FEM (Ω0) and the tolerance model (Ω) is denoted by:
ε = | Ω 0 Ω _ Ω _ | 100 % .
Results are obtained for these periodic plate strips for ratios or parameters: Poisson’s ratios ν″ = ν′ = 0.3; E″/E′ = 0.25, 0.5, 0.75; ρ″/ρ′ = 0.2, 0.5, 0.8; γ = 0.2, 0.5, 0.8. Additionally, the first lower frequency parameter for the proper homogeneous plate strip is shown with these boundary conditions (ν″ = ν′ = 0.3; E″/E′ = 1; ρ″/ρ′ = 1; γ = 1) calculated using the tolerance model (Ω = 0.02987), the classical analytical solution (ΩC = 0.02987), and FEM (Ω0 = 0.02986) [%].
Table A1. Frequency parameters Ω, Ω0 for lower free vibration frequencies for the simply supported plate strip (α = π, ν = 0.3, l/L = 0.1, h/l = 0.1).
Table A1. Frequency parameters Ω, Ω0 for lower free vibration frequencies for the simply supported plate strip (α = π, ν = 0.3, l/L = 0.1, h/l = 0.1).
γ E’’/E’ = 0.25E’’/E’ = 0.5E’’/E’ = 0.75E’/E’’ = 1
ρ’’/ρ’ΩΩ0ε [%]ΩΩ0ε [%]ΩΩ0ε [%]ΩΩ0ε [%]
0.20.20.028300.027383.360.037710.038111.050.044380.044450.160.049780.049770.02
0.50.021920.021203.400.029210.029511.020.034370.034430.170.038560.038550.03
0.80.018530.017923.400.024690.024941.000.029050.029100.170.032590.032580.03
1.00.016980.016423.410.022630.022861.010.026630.026670.150.029870.029860.03
0.50.20.025650.026623.640.031850.031730.380.035770.035700.200.038560.038550.03
0.50.022940.023803.610.028490.028380.390.031990.031930.190.034490.034480.03
0.80.020940.021723.590.026010.025900.420.029200.029140.210.031480.031480.00
1.00.019870.020603.540.024670.024570.410.027700.027650.180.029870.029860.03
0.80.20.029080.029621.820.030550.030121.430.031680.031550.410.032590.032580.03
0.50.028100.028611.780.029510.029101.410.030610.030480.000.031480.031480.00
0.80.027200.027701.810.028570.028181.380.029630.029510.410.030480.030480.00
1.00.026650.027141.810.028000.027611.410.029040.028910.450.029870.029860.03
The curves of the parameter of differences ε versus the ratios ρ″/ρ′ are shown in Figure A1.
Figure A1. The parameter of differences ε (A1) versus the ratios ρ″/ρ′ (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1; ν″ = ν′ = 0.3; E″/E′ = 0.25, 0.5, 0.75; γ = 0.2, 0.5, 0.8).
Figure A1. The parameter of differences ε (A1) versus the ratios ρ″/ρ′ (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1; ν″ = ν′ = 0.3; E″/E′ = 0.25, 0.5, 0.75; γ = 0.2, 0.5, 0.8).
Materials 15 05623 g0a1
Analysing results shown in Table A1 and Figure A1, it can be observed that:
  • Using the finite element method, only the lower frequency parameters, obtained in the framework of the tolerance model, can be compared and justified.
  • Differences between values of these parameters calculated using both methods are smaller than about 4% (3.64%) for the assumed geometrical and material parameters.
  • Smaller differences between obtained results (below 1.5%) are for these plate strips, which are described by the ratio E″/E′ ≥ 0.5 (E″/E′ = 0.5, 0.75, 1) and various values of the ratios—γ = 0.2, 0.5, 0.8; ρ″/ρ′ = 0.2, 0.5, 0.8, 1.
  • The biggest differences between calculated lower frequency parameters (above 3%) are for smaller values of the ratio E″/E′ (e.g., E″/E′ = 0.25) and some smaller values of the distribution parameter of material properties γ (e.g., γ = 0.2, 0.5).

References

  1. Kohn, R.V.; Vogelius, M. A new model of thin plates with rapidly varying thickness. Int. J. Solids Struct. 1984, 20, 333–350. [Google Scholar] [CrossRef]
  2. Matysiak, S.J.; Nagórko, W. Microlocal parameters in the modelling of microperiodic plates. Ing. Arch. 1989, 59, 434–444. [Google Scholar] [CrossRef]
  3. Saha, G.C.; Kalamkarov, A.L.; Georgiades, A.V. Effective elastic characteristics of honeycomb sandwich composite shells made of generally orthotropic materials. Compos. Part A Appl. Sci. Manufact. 2007, 38, 1533–1546. [Google Scholar] [CrossRef]
  4. Dallot, J.; Sab, K.; Foret, G. Limit analysis of periodic beams. Eur. J. Mech.-A/Sol. 2009, 28, 166–178. [Google Scholar] [CrossRef]
  5. Królak, M.; Kowal-Michalska, K.; Mania, R.J.; Świniarski, J. Stability and load carrying capacity of multi-cell thin-walled columns of rectangular cross-sections. J. Theor. Appl. Mech. 2009, 47, 435–456. [Google Scholar]
  6. Teter, A. Dynamic critical load based on different stability criteria for coupled buckling of columns with stiffened open cross-sections. Thin-Walled Struct. 2011, 49, 589–595. [Google Scholar] [CrossRef]
  7. Jasion, P.; Magnucka-Blandzi, E.; Szyc, W.; Magnucki, K. Global and local buckling of sandwich circular and beam-rectangular plates with metal foam core. Thin-Walled Struct. 2012, 61, 154–161. [Google Scholar] [CrossRef]
  8. Schmitz, A.; Horst, P. A finite element unit-cell method for homogenised mechanical properties of heterogeneous plates. Compos. Part A Appl. Sci. Manufact. 2014, 61, 23–32. [Google Scholar] [CrossRef]
  9. Brito-Santana, H.; Wang, Y.S.; Rodríguez-Ramos, R.; Bravo-Castillero, J.; Guinovart-Díaz, R.; Tita, V. A dispersive nonlocal model for shear wave propagation in laminated composites with periodic structures. Eur. J. Mech.-A/Sol. 2015, 49, 35–48. [Google Scholar] [CrossRef]
  10. Grygorowicz, M.; Magnucki, K.; Malinowski, M. Elastic buckling of a sandwich beam with variable mechanical properties of the core. Thin-Walled Struct. 2015, 87, 127–132. [Google Scholar] [CrossRef]
  11. Grygorowicz, M.; Magnucka-Blandzi, E. Mathematical modeling for dynamic stability of sandwich beam with variable mechanical properties of core. Appl. Math. Mech. 2016, 37, 361–374. [Google Scholar] [CrossRef]
  12. Fantuzzi, N.; Tornabene, F. Strong Formulation Isogeometric Analysis (SFIGA) for laminated composite arbitrarily shaped plates. Compos. Part B Eng. 2016, 96, 173–203. [Google Scholar] [CrossRef]
  13. Liu, B.; Ferreira, A.J.M.; Xing, Y.F.; Neves, A.M.A. Analysis of composite plates using a layerwise theory and a differential quadrature finite element method. Compos. Struct. 2016, 156, 393–398. [Google Scholar] [CrossRef]
  14. Sobhani, E.; Masoodi, A.R. Natural frequency responses of hybrid polymer/carbon fiber/FG-GNP nanocomposites paraboloidal and hyperboloidal shells based on multiscale approaches. Aerosp. Sci. Technol. 2021, 119, 107111. [Google Scholar] [CrossRef]
  15. Sobhani, E.; Masoodi, A.R. On the circumferential wave responses of connected elliptical-cylindrical shell-like submerged structures strengthened by nano-reinforcer. Ocean. Eng. 2022, 247, 110718. [Google Scholar] [CrossRef]
  16. Sobhani, E.; Moradi-Dastjerdi, R.; Behdinan, K.; Masoodi, A.R.; Ahmadi-Pari, A.R. Multifunctional trace of various reinforcements on vibrations of three-phase nanocomposite combined hemispherical-cylindrical shells. Compos. Struct. 2022, 279, 114798. [Google Scholar] [CrossRef]
  17. Jopek, H.; Strek, T. Torsion of a two-phased composite bar with helical distribution of constituents. Phys. Status Solidi 2017, 254, 1700050. [Google Scholar] [CrossRef]
  18. Brillouin, L. Wave Propagation in Periodic Structures; Dover Pub. Inc.: Dover, UK, 1953. [Google Scholar]
  19. Wu, Z.-J.; Li, F.-M.; Wang, Y.-Z. Vibration band gap properties of periodic Mindlin plate structure using the spectral element method. Meccanica 2014, 49, 725–737. [Google Scholar] [CrossRef]
  20. Zhou, X.Q.; Yu, D.Y.; Shao, X.; Wang, S.; Tian, Y.H. Band gap characteristics of periodically stiffened-thin-plate based on center-finite-difference-method. Thin-Walled Struct. 2014, 82, 115–123. [Google Scholar] [CrossRef]
  21. Zhou, X.Q.; Yu, D.Y.; Shao, X.; Wang, S.; Zhang, S.Q. Simplified-super-element-method for analyzing free flexural vibration characteristics of periodically stiffened-thin-plate filled with viscoelastic damping material. Thin-Walled Struct. 2015, 94, 234–252. [Google Scholar] [CrossRef]
  22. Cheng, Z.B.; Xu, Y.G.; Zhang, L.L. Analysis of flexural wave bandgaps in periodic plate structures using differential quadrature element method. Int. J. Mech. Sci. 2015, 100, 112–125. [Google Scholar] [CrossRef]
  23. Woźniak, C.; Wierzbicki, E. Averaging Techniques in Thermomechanics of Composite Solids. Tolerance Averaging Versus Homogenization; Publishing House of Częstochowa University of Technology: Częstochowa, Poland, 2000. [Google Scholar]
  24. Woźniak, C.; Michalak, B.; Jędrysiak, J. (Eds.) Thermomechanics of Microheterogeneous Solids and Structures. Tolerance Averaging Approach; Publishing House, Łódź University of Technology: Lodz, Poland, 2008. [Google Scholar]
  25. Wierzbicki, E.; Woźniak, C. On the dynamics of combined plane periodic structures. Arch. Appl. Mech. 2000, 70, 387–398. [Google Scholar] [CrossRef]
  26. Jędrysiak, J. On vibrations of thin plates with one-dimensional periodic structure. Int. J. Eng. Sci. 2000, 38, 2023–2043. [Google Scholar] [CrossRef]
  27. Michalak, B. The meso-shape functions for the meso-structural models of wavy-plates. ZAMM 2001, 81, 639–641. [Google Scholar] [CrossRef]
  28. Nagórko, W.; Woźniak, C. Nonasymptotic modelling of thin plates reinforced by a system of stiffeners. Electr. J. Polish Agric. Univ. Civil Eng. 2002, 5, 8. [Google Scholar]
  29. Baron, E. On dynamic behaviour of medium-thickness plates with uniperiodic structure. Arch. Appl. Mech. 2003, 73, 505–516. [Google Scholar] [CrossRef]
  30. Jędrysiak, J. The length-scale effect in the buckling of thin periodic plates resting on a periodic Winkler foundation. Meccanica 2003, 38, 435–451. [Google Scholar] [CrossRef]
  31. Mazur-Śniady, K.; Woźniak, C.; Wierzbicki, E. On the modelling of dynamic problems for plates with a periodic structure. Arch. Appl. Mech. 2004, 74, 179–190. [Google Scholar] [CrossRef]
  32. Tomczyk, B. A non-asymptotic model for the stability analysis of thin biperiodic cylindrical shells. Thin-Walled Struct. 2007, 45, 941–944. [Google Scholar] [CrossRef]
  33. Tomczyk, B. Dynamic stability of micro-periodic cylindrical shells. Mech. Mech. Eng. 2010, 14, 137–150. [Google Scholar]
  34. Jędrysiak, J.; Paś, A. Dynamics of medium thickness plates interacting with a periodic Winkler’s foundation: Non-asymptotic tolerance modeling. Meccanica 2014, 49, 1577–1585. [Google Scholar] [CrossRef]
  35. Domagalski, Ł.; Jędrysiak, J. On the tolerance modelling of geometrically nonlinear thin periodic plates. Thin-Walled Struct. 2015, 87, 183–190. [Google Scholar] [CrossRef]
  36. Domagalski, Ł.; Jędrysiak, J. Geometrically nonlinear vibrations of slender meso-periodic beams. The tolerance modeling approach. Compos. Struct. 2016, 136, 270–277. [Google Scholar] [CrossRef]
  37. Marczak, J.; Jędrysiak, J. Some remarks on modelling of vibrations of periodic sandwich structures with inert core. Compos. Struct. 2018, 202, 752–758. [Google Scholar] [CrossRef]
  38. Jędrysiak, J. Modelling of Vibrations and Stability for Periodic Slender Visco-Elastic Beams on a Foundation with Damping. Revisiting. Materials 2020, 13, 3939. [Google Scholar] [CrossRef]
  39. Marczak, J.; Michalak, B.; Wirowski, A. A multi-scale analysis of stress distribution in thin composite plates with dense system of ribs in two directions. Adv. Eng. Softw. 2021, 153, 102960. [Google Scholar] [CrossRef]
  40. Michalak, B.; Wirowski, A. Dynamic modelling of thin plate made of certain functionally graded materials. Meccanica 2012, 47, 1487–1498. [Google Scholar] [CrossRef]
  41. Wirowski, A.; Michalak, B.; Gajdzicki, M. Dynamic modelling of annular plates of functionally graded structure resting on elastic heterogeneous foundation with two modules. J. Mech. 2015, 31, 493–504. [Google Scholar] [CrossRef]
  42. Perliński, W.; Gajdzicki, M.; Michalak, B. Modelling of annular plates stability with functionally graded structure interacting with elastic heterogeneous subsoil. J. Theor. Appl. Mech. 2014, 52, 485–498. [Google Scholar]
  43. Michalak, B. 2D tolerance and asymptotic models in elastodynamics of a thin-walled structure with dense system of ribs. Arch. Civ. Mech. Eng. 2015, 15, 449–455. [Google Scholar] [CrossRef]
  44. Rabenda, M.; Michalak, B. Natural vibrations of prestressed thin functionally graded plates with dense system of ribs in two directions. Compos. Struct. 2015, 133, 1016–1023. [Google Scholar] [CrossRef]
  45. Ostrowski, P.; Michalak, B. The combined asymptotic-tolerance model of heat conduction in a skeletal micro-heterogeneous hollow cylinder. Compos. Struct. 2015, 134, 343–352. [Google Scholar] [CrossRef]
  46. Ostrowski, P.; Michalak, B. A contribution to the modelling of heat conduction for cylindrical composite conductors with non-uniform distribution of constituents. Int. J. Heat Mass Transf. 2016, 92, 435–448. [Google Scholar] [CrossRef]
  47. Pazera, E.; Jędrysiak, J. Effect of microstructure in thermoelasticity problems of functionally graded laminates. Compos. Struct. 2018, 202, 296–303. [Google Scholar] [CrossRef]
  48. Jędrysiak, J. Tolerance modelling of free vibration frequencies of thin functionally graded plates with one-directional microstructure. Compos. Struct. 2017, 161, 453–468. [Google Scholar] [CrossRef]
  49. Jędrysiak, J. Tolerance modelling of free vibrations of medium thickness functionally graded plates. Compos. Struct. 2018, 202, 1253–1262. [Google Scholar] [CrossRef]
  50. Tomczyk, B.; Szczerba, P. Tolerance and asymptotic modelling of dynamic problems for thin microstructured transversally graded shells. Compos. Struct. 2017, 162, 365–373. [Google Scholar] [CrossRef]
  51. Tomczyk, B.; Szczerba, P. Combined asymptotic-tolerance modelling of dynamic problems for functionally graded shells. Compos. Struct. 2018, 183, 176–184. [Google Scholar] [CrossRef]
  52. Tomczyk, B.; Szczerba, P. A new asymptotic-tolerance model of dynamic and stability problems for longitudinally graded cylindrical shells. Compos. Struct. 2018, 202, 473–481. [Google Scholar] [CrossRef]
  53. Jędrysiak, J.; Kaźmierczak-Sobińska, M. Theoretical Analysis of Buckling for Functionally Graded Thin Plates with Microstructure Resting on an Elastic Foundation. Materials 2020, 13, 4031. [Google Scholar] [CrossRef]
Figure 1. A fragment of a thin periodic plate.
Figure 1. A fragment of a thin periodic plate.
Materials 15 05623 g001
Figure 2. A fragment of a thin periodic plate strip.
Figure 2. A fragment of a thin periodic plate strip.
Materials 15 05623 g002
Figure 3. A periodicity cell of the plate strip under consideration.
Figure 3. A periodicity cell of the plate strip under consideration.
Materials 15 05623 g003
Figure 4. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (ρ″/ρ′; E″/E′) (for h/l = 0.1, l/L = 0.1) for: (a) a simply supported plate strip, (b) a clamped-hinged plate strip.
Figure 4. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (ρ″/ρ′; E″/E′) (for h/l = 0.1, l/L = 0.1) for: (a) a simply supported plate strip, (b) a clamped-hinged plate strip.
Materials 15 05623 g004
Figure 5. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (ρ″/ρ′; E″/E′) (for h/l = 0.1, l/L = 0.1) for: (a) a clamped-clamped plate strip, (b) a cantilever plate strip.
Figure 5. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (ρ″/ρ′; E″/E′) (for h/l = 0.1, l/L = 0.1) for: (a) a clamped-clamped plate strip, (b) a cantilever plate strip.
Materials 15 05623 g005
Figure 6. Traces of surfaces of dimensionless lower frequency parameters (46)1,3 from Figure 4 and Figure 5 on planes of these frequency parameters for the homogeneous plate strips (γ = 1) (for h/l = 0.1, l/L = 0.1) for all boundary conditions under consideration.
Figure 6. Traces of surfaces of dimensionless lower frequency parameters (46)1,3 from Figure 4 and Figure 5 on planes of these frequency parameters for the homogeneous plate strips (γ = 1) (for h/l = 0.1, l/L = 0.1) for all boundary conditions under consideration.
Materials 15 05623 g006
Figure 7. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; h/L) (for E″/E′ = 0.5, l/L = 0.1) for: (a) a simply supported plate strip, (b) a clamped-hinged plate strip.
Figure 7. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; h/L) (for E″/E′ = 0.5, l/L = 0.1) for: (a) a simply supported plate strip, (b) a clamped-hinged plate strip.
Materials 15 05623 g007
Figure 8. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; h/L) (for E″/E′ = 0.5, l/L = 0.1) for: (a) a clamped-clamped plate strip, (b) a cantilever plate strip.
Figure 8. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; h/L) (for E″/E′ = 0.5, l/L = 0.1) for: (a) a clamped-clamped plate strip, (b) a cantilever plate strip.
Materials 15 05623 g008
Figure 9. Dimensionless higher frequency parameters (46)2 versus: (a) pairs of ratios (ρ″/ρ′;E″/E′) (for h/l = 0.1, l/L = 0.1), (b) pairs of ratios (γ; h/L) (for E″/E′ = 0.5, l/L = 0.1).
Figure 9. Dimensionless higher frequency parameters (46)2 versus: (a) pairs of ratios (ρ″/ρ′;E″/E′) (for h/l = 0.1, l/L = 0.1), (b) pairs of ratios (γ; h/L) (for E″/E′ = 0.5, l/L = 0.1).
Materials 15 05623 g009
Figure 10. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; E″/E′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1) for: (a) a simply supported plate strip, (b) a clamped-hinged plate strip.
Figure 10. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; E″/E′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1) for: (a) a simply supported plate strip, (b) a clamped-hinged plate strip.
Materials 15 05623 g010
Figure 11. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; E″/E′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1) for: (a) a clamped-clamped plate strip, (b) a cantilever plate strip.
Figure 11. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; E″/E′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1) for: (a) a clamped-clamped plate strip, (b) a cantilever plate strip.
Materials 15 05623 g011
Figure 12. Traces of surfaces of dimensionless lower frequency parameters (46)1,3 from Figure 10 and Figure 11 on planes of these frequency parameters for the homogeneous plate strips (γ = 1) (for ν″/ν′ = 1, h/l = 0.1, l/L = 0.1) for all boundary conditions under consideration.
Figure 12. Traces of surfaces of dimensionless lower frequency parameters (46)1,3 from Figure 10 and Figure 11 on planes of these frequency parameters for the homogeneous plate strips (γ = 1) (for ν″/ν′ = 1, h/l = 0.1, l/L = 0.1) for all boundary conditions under consideration.
Materials 15 05623 g012
Figure 13. Dimensionless higher frequency parameters (46)2 versus pairs of ratios (γ; E″/E′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1).
Figure 13. Dimensionless higher frequency parameters (46)2 versus pairs of ratios (γ; E″/E′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1).
Materials 15 05623 g013
Figure 14. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; ν″/ν′) (for E″/E′ = 0.5, l/L = 0.1, h/l = 0.1) for: (a) a simply supported plate strip, (b) a clamped-hinged plate strip.
Figure 14. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; ν″/ν′) (for E″/E′ = 0.5, l/L = 0.1, h/l = 0.1) for: (a) a simply supported plate strip, (b) a clamped-hinged plate strip.
Materials 15 05623 g014
Figure 15. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; ν″/ν′) (for E″/E′ = 0.5, l/L = 0.1, h/l = 0.1) for: (a) a clamped-clamped plate strip, (b) a cantilever plate strip.
Figure 15. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; ν″/ν′) (for E″/E′ = 0.5, l/L = 0.1, h/l = 0.1) for: (a) a clamped-clamped plate strip, (b) a cantilever plate strip.
Materials 15 05623 g015
Figure 16. Dimensionless higher frequency parameters (46)2 versus pairs of ratios (γ; ν″/ν′) (for E″/E′ = 0.5, l/L = 0.1, h/l = 0.1).
Figure 16. Dimensionless higher frequency parameters (46)2 versus pairs of ratios (γ; ν″/ν′) (for E″/E′ = 0.5, l/L = 0.1, h/l = 0.1).
Materials 15 05623 g016
Figure 17. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; ρ″/ρ′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1) for: (a) a simply supported plate strip, (b) a clamped-hinged plate strip.
Figure 17. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; ρ″/ρ′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1) for: (a) a simply supported plate strip, (b) a clamped-hinged plate strip.
Materials 15 05623 g017
Figure 18. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; ρ″/ρ′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1) for: (a) a clamped-clamped plate strip, (b) a cantilever plate strip.
Figure 18. Dimensionless lower frequency parameters (46)1,3 versus pairs of ratios (γ; ρ″/ρ′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1) for: (a) a clamped-clamped plate strip, (b) a cantilever plate strip.
Materials 15 05623 g018
Figure 19. Dimensionless higher frequency parameters (46)2 versus pairs of ratios (γ; ρ″/ρ′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1).
Figure 19. Dimensionless higher frequency parameters (46)2 versus pairs of ratios (γ; ρ″/ρ′) (for ν″/ν′ = 1, l/L = 0.1, h/l = 0.1).
Materials 15 05623 g019
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jędrysiak, J. The Effect of the Material Periodic Structure on Free Vibrations of Thin Plates with Different Boundary Conditions. Materials 2022, 15, 5623. https://doi.org/10.3390/ma15165623

AMA Style

Jędrysiak J. The Effect of the Material Periodic Structure on Free Vibrations of Thin Plates with Different Boundary Conditions. Materials. 2022; 15(16):5623. https://doi.org/10.3390/ma15165623

Chicago/Turabian Style

Jędrysiak, Jarosław. 2022. "The Effect of the Material Periodic Structure on Free Vibrations of Thin Plates with Different Boundary Conditions" Materials 15, no. 16: 5623. https://doi.org/10.3390/ma15165623

APA Style

Jędrysiak, J. (2022). The Effect of the Material Periodic Structure on Free Vibrations of Thin Plates with Different Boundary Conditions. Materials, 15(16), 5623. https://doi.org/10.3390/ma15165623

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop