Open access peer-reviewed chapter

# Distinctive Characteristics of Cosserat Plate Free Vibrations

Written By

Lev Steinberg and Roman Kvasov

Submitted: May 1st, 2019 Reviewed: May 24th, 2019 Published: July 10th, 2019

DOI: 10.5772/intechopen.87044

From the Edited Volume

## Dynamical Systems Theory

Edited by Jan Awrejcewicz and Dariusz Grzelczyk

Chapter metrics overview

View Full Metrics

## Abstract

In this chapter, we present the theoretical analysis of the distinctive characteristics of Cosserat plate vibrations. This analysis is based on the dynamic model of the Cosserat plates, which we developed as an extension of the Reissner plate theory. Primarily, we describe the validation of the model, which is based on the comparison with three-dimensional exact solutions. We present the results of the computer simulations, which allow us to identify different characteristics of the plate vibrations. Particularly, we illustrate and discuss the detection and the classification of the additional high resonance frequencies of a plate depending on the shape and orientation of microelements incorporated into the Cosserat plates.

### Keywords

• variational principle
• Cosserat plate vibrations
• frequencies of micro-vibrations

## 1. Introduction

The theory of asymmetric elasticity introduced in 1909 by the Cosserat brothers [1] gave rise to a variety of Cosserat plate theories. In 1960s, Green and Naghdi specialized their general theory of Cosserat surface to obtain the linear Cosserat plate [2], while independently Eringen proposed a complete theory of plates in the framework of Cosserat elasticity [3]. Numerous plate theories were formulated afterwards; for the review of the latest developments in the area of Cosserat plates we recommend to turn to [4].

The first theory of Cosserat plates based on the Reissner plate theory was developed in [5] and its finite element modeling is provided in [6]. The parametric theory of Cosserat plate, presented by the authors in [7], includes some additional assumptions leading to the introduction of the splitting parameter. This provided the highest level of approximation to the original three-dimensional problem. The theory provides the equilibrium equations and constitutive relations, and the optimal value of the minimization of the elastic energy of the Cosserat plate. The paper [7] also provides the analytical solutions of the presented plate theory and the three-dimensional Cosserat elasticity for simply supported rectangular plate. The comparison of these solutions showed that the precision of the developed Cosserat plate theory is similar to the precision of the classical plate theory developed by Reissner [8, 9].

The numerical modeling of bending of simply supported rectangular plates is given in [10]. We developed the Cosserat plate field equations and a rigorous formula for the optimal value of the splitting parameter. The solution of the Cosserat plate was shown to converge to the Reissner plate as the elastic asymmetric parameters tend to zero. The Cosserat plate theory demonstrates the agreement with the size effect, confirming that the plates of smaller thickness are more rigid than is expected from the Reissner model. The modeling of Cosserat plates with simply supported rectangular holes is also provided. The finite element analysis of the perforated Cosserat plates is given in [11].

The extension of the static model of Cosserat elastic plates to the dynamic problems is presented in [12]. The computations predict a new kind of natural frequencies associated with the material microstructure and were shown to be compatible with the size effect principle reported in [10] for the Cosserat plate bending.

This chapter represents an extension of the paper [12] for different shapes and orientations of micro-elements incorporated into the Cosserat plates. It is based on the generalized variational principle for elastodynamics and includes a non-diagonal rotatory inertia tensor. The numerical computations of the plate free vibrations showed the existence of some additional high frequencies of micro-vibrations depending on the orientation of micro-elements. The comparison with three-dimensional Cosserat elastodynamics shows a high agreement with the exact values of the eigenvalue frequencies.

## 2. Cosserat linear elastodynamics

### 2.1 Fundamental equations

The Cosserat linear elasticity balance laws are

σ ji , j = p i t , E1
ε ijk σ jk + μ ji , j = q i t , E2

where the σ ji is the stress tensor, μ ji the couple stress tensor, p i = ρ u i t and q i = J ji ϕ j t are the linear and angular momenta, ρ and J ji are the material density and the rotatory inertia characteristics, ε ijk is the Levi-Civita tensor.

We will also consider the constitutive equations as in [13]:

σ ji = μ + α γ ji + μ α γ ij + λγ kk δ ij , E3
μ ji = γ + ε χ ji + γ ε χ ij + βχ kk δ ij , E4

and the kinematic relations in the form

γ ji = u i , j + ε ijk ϕ k and χ ji = ϕ i , j , E5

Here u i and ϕ i represent the displacement and microrotation vectors, γ ji and χ ji represent the strain and bend-twist tensors, μ , λ are the Lamé parameters and α , β , γ , ε are the Cosserat elasticity parameters.

The constitutive Eqs. (3)(4) can be written in the reverse form [5].

γ ji = μ + α σ ji + μ α σ ij + λ σ kk , E6
χ ji = γ + ε μ ji + γ ε μ ij + β μ kk , E7

where μ = 1 4 μ , α = 1 4 α , γ = 1 4 γ , ε = 1 4 ε , λ = λ 6 μ λ + 2 μ 3 and β = β 6 μ β + 2 γ 3 .

We will consider the boundary conditions given in [12].

u i = u i 0 , ϕ i = ϕ i 0 , on G 1 t = B 0 \ B σ × t 0 t , E8
σ ji n i = σ j 0 , μ ji n i = μ j 0 on G 2 t = B σ × t 0 t , E9

and initial conditions

u i x 0 = U i 0 , ϕ i x 0 = Φ i 0 , in B 0 , E10
u ̇ i x 0 = U ̇ i 0 , ϕ ̇ i x 0 = Φ ̇ i 0 , in B 0 , E11

where u i 0 and ϕ i 0 are prescribed on G 1 , σ j 0 and μ j 0 on G 2 , and n i is the unit vector normal to the boundary B 0 of the elastic body B 0 .

### 2.2 Cosserat elastic energy

The strain stored energy U C of the body B 0 is defined by the integral [13]:

U C = B 0 W γ χ dv , E12

where

W γ χ = μ + α 2 γ ij γ ij + μ α 2 γ ij γ ji + λ 2 γ kk γ nn + γ + ε 2 χ ij χ ij + γ ε 2 χ ij χ ji + β 2 χ kk χ nn , E13

is non-negative. The relations Eqs. (3)(4) can be written in the form [12]:

σ = γ W and μ = χ W . E14

The stress energy is given as

U K = B 0 Φ σ μ dv , E15

where

Φ σ μ = μ + α 2 σ ij σ ij + μ α 2 σ ij σ ji + λ 2 σ kk σ nn + γ + ε 2 μ ij μ ij + γ ε 2 μ ij μ ji + β 2 μ kk μ nn , E16

and the relations Eqs. (6)(7) can be written as [12].

γ = ∂Φ σ , and χ = ∂Φ μ . E17

We consider the work done by the stresses σ and μ over the strains γ and χ as in [13].

U = B 0 σ γ + μ χ dv E18

and

U = U K = U C E19

Here σ γ = σ ji γ ji and μ χ = μ ji χ ji .

The stored kinetic energy T C is defined as

T C = B 0 ϒ C dv = 1 2 B 0 ρ u t 2 + J ϕ t 2 dv , E20

The kinetic energy T K is given as

T K = B 0 ϒ K p q dv = 1 2 B 0 p 2 ρ 1 + q 2 J 1 dv , E21

where

p = ϒ C u ̇ = ρ u t and q = ϒ C ϕ ̇ = J ϕ t , E22

and

u t = ϒ K p = p ρ 1 and ϕ t = ϒ K q = qJ 1 , E23

The work T W done by the inertia forces over displacement and microrotation is given as in [12].

T W = B 0 ϒ W dv = B 0 p t u + q t ϕ dv E24

Keeping in mind that the variation of p u , q , ϕ , δ u , and δ ϕ is zero at t 0 and t k we can integrate by parts

t 0 t k T K dt = t 0 t k T I dt = 1 2 B 0 p u + q ϕ dv t 0 t k t 0 t k T W dt E25
δ t 0 t k T K = δ t 0 t k T W E26

or

δ T C = δ T K = δ T W E27

and therefore

t 0 t k B 0 p δ u t + q δ ϕ t dvdt = t 0 t k B 0 p t δ u + q t δ ϕ dvdt E28

### 2.3 Variational principle for elastodynamics

We modify the HPR principle [14] for the case of Cosserat elastodynamics in the following way: for any set A v of all admissible states s = u ϕ γ χ σ μ that satisfy the strain-displacement and torsion-rotation relations Eq. (5), the zero variation

δ Θ s = 0

of the functional

Θ s = t 0 t k U K + T C B 0 σ γ + μ χ + p u t + q ϕ t dv dt + t 0 t G 1 σ n u u 0 + μ n ϕ ϕ 0 dadt + t 0 t G 2 σ 0 u + μ 0 ϕ dadt E29

at s A is equivalent of s to be a solution of the system of equilibrium Eqs. (1)(2), constitutive relations Eqs. (6)(7), which satisfies the mixed boundary conditions Eqs. (8)(9).

Proof of the variational principle for elastodynamics

Let us consider the variation of the functional Θ s :

δ Θ s = t 0 t k δ U K + δ T C dt t 0 t k B 0 δ σ γ + σ δ γ + δ μ χ + μ δ χ + u t δ p + p δ u t + δ q ϕ t + q δ ϕ t dvdt + t 0 t G 1 δ σ n u u 0 + σ n δ u + δ μ n ϕ ϕ 0 + μ n δ ϕ dadt + t 0 t G 2 σ 0 δ u + μ 0 δ ϕ dadt

Taking into account Eq. (5) we can perform the integration by parts

B 0 σ δ γ dv = B 0 σ n δ u da B 0 δ u div σ dv + B 0 εσ δ ϕ dv B 0 μ δ χ dv = B 0 μ n δ ϕ da B 0 δ ϕ div μ dv

and based on Eqs. (17)(23)

δ Φ = ∂Φ σ δ σ + ∂Φ μ δ μ , δ ϒ C = ϒ C u t δ u t + ϒ C ϕ t δ ϕ t .

Then keeping in mind that δ T K = δT and Eq. (28) we can rewrite the expression for the variation of the functional δ Θ s in the following form

δ Θ s = t 0 t B 0 ∂Φ σ γ δ σ dvdt + t 0 t B 0 ∂Φ μ χ δ μ dvdt + t 0 t B 0 ρ u t p δ u t dvdt + t 0 t B 0 J φ t q δ φ t dvdt + t 0 t B 0 div σ p t δ u dvdt + t 0 t B 0 div μ + ε σ q t δ u dvdt + t 0 t G 1 u u 0 δ σ n dadt + t 0 t G 1 ϕ ϕ 0 δ μ n dadt + t 0 t G 2 σ n σ 0 δ u dadt + t 0 t G 2 μ n μ 0 δ ϕ dadt

## 3. Dynamic Cosserat plate theory

In this section we review our stress, couple stress and kinematic assumptions of the Cosserat plate [7]. We consider the thin plate P , where h is the thickness of the plate and x 3 = 0 represents its middle plane. The sets T and B are the top and bottom surfaces contained in the planes x 3 = h / 2 , x 3 = h / 2 respectively and the curve Γ is the boundary of the middle plane of the plate.

The set of points P = Γ × [ h 2 h 2 ] T B forms the entire surface of the plate and Γ u × [ h 2 , h 2 ] is the lateral part of the boundary where displacements and microrotations are prescribed. The notation Γ σ = Γ \ Γ u of the remainder we use to describe the lateral part of the boundary edge Γ σ × [ h 2 , h 2 ] where stress and couple stress are prescribed. We also use notation P 0 for the middle plane internal domain of the plate.

In our case we consider the vertical load and pure twisting momentum boundary conditions at the top and bottom of the plate, which can be written in the form:

σ 33 x 1 x 2 h / 2 t = σ t x 1 x 2 t , σ 33 x 1 x 2 h / 2 t = σ b x 1 x 2 t , E30
σ 3 β x 1 x 2 ± h / 2 t = 0 , E31
μ 33 x 1 x 2 h / 2 t = μ t x 1 x 2 t , μ 33 x 1 x 2 h / 2 t = μ b x 1 x 2 t , E32
μ 3 β x 1 x 2 ± h / 2 t = 0 , E33

where x 1 x 2 P 0 .

We will also consider the rotatory inertia J in the form

J = J 11 J 12 0 J 12 J 22 0 0 0 J 33

Let A denote the set of all admissible states that satisfy the Cosserat plate strain-displacement relation Eq. (5) and let Θ be a functional on A defined by

Θ s η = U K S + T C S P 0 S E + P U t P ̂ W + v Ω 3 0 da + Γ σ S o U U o ds + Γ u S n U ds , E34

for every s = U E S A . Here P ̂ = p ̂ 1 p ̂ 2 and W = W W , p ̂ 1 = ηp and p ̂ 2 = 2 3 1 η p

Here the plate stress and kinetic energy density by the formulas

U K S = P 0 Φ S da , T K S = P 0 ϒ C U t da E35

where P 0 is the internal domain of the middle plane of the plate.

Φ S = 3 λ M αα M ββ h 3 μ 3 λ + 2 μ + 3 α + μ M αβ 2 2 h 3 αμ + 3 α + μ 160 h 3 αμ 8 Q ̂ α Q ̂ α + 15 Q α Q ̂ α + 20 Q ̂ α Q α + 8 Q α Q α + 3 α μ M αβ 2 2 h 3 αμ + α μ 280 h 3 αμ 21 Q α 5 Q ̂ α + 4 Q α γ ε 160 hγε 24 R αα 2 + 45 R αα + 60 R αβ R αβ + 48 R 12 R 21 + 3 γ + ε S α S α 2 h 3 γε + γ + ε 160 h 3 γε 8 R αβ 2 + 15 R αβ R αβ + 20 R αβ R αβ 3 β 80 3 β + 2 γ 8 R αα R ββ 15 R αα R ββ 20 R αα R αα β 4 γ 3 β + 2 γ 2 R αα + 3 R αα t h V 2 + T 2 + λ 560 3 λ + 2 μ 5 + 3 η 1 + η pM αα + λ + μ h 840 μ 3 λ + 2 μ 140 + 168 η + 51 η 2 4 1 + η 2 p 2 + λ + μ h 2 μ 3 λ + 2 μ σ 0 2 + εh 12 hγε 3 T 2 + V 2 E36

and

ϒ C U t = 2 W t 2 + 4 15 W t 2 + 2 3 W t W t + h 3 ρ 24 Ψ α t 2 + 4 hJ αβ 15 Ω α 0 t Ω β 0 t + hJ αβ 2 Ω ̂ α 0 t Ω ̂ β 0 t + 2 hJ αβ 3 Ω α 0 t Ω ̂ β 0 t + hJ 33 6 Ω 3 t 2 .

S , U and E are the Cosserat plate stress, displacement and strain sets

S = M αβ Q α Q α Q ̂ α R αβ R αβ S β , E37
S n = M ˇ α Q ˇ Q ˇ α ̂ R ˇ α R ˇ α S ˇ , E38
S o = Π Π o 3 Π o 3 M M M o 3 , E39
U = Ψ α W Ω 3 Ω α 0 W Ω α 0 , E40
E = e αβ ω β ω α ω ̂ α τ 3 α τ αβ τ αβ , E41

where M αβ n β = Π , R αβ n β = M , Q α n α = Π o 3 , S α n α = M o 3 , Q ̂ α n α = Π o 3 , R αβ n β = M , M ˇ α = M αβ n β , Q ˇ = Q β n β , R ˇ α = R αβ n β , S ˇ = S β n β , Q ˇ ̂ = Q ˇ β ̂ n β , R ˇ α = R ˇ αβ n β . ( n β is the outward unit normal vector to Γ u ).

The plate characteristics provide the approximation of the components of the three-dimensional tensors σ ji and μ ji

σ αβ = 6 h 2 ζ M αβ x 1 x 2 t , E42
σ 3 β = 3 2 h 1 ζ 2 Q β x 1 x 2 t , E43
σ β 3 = 3 2 h 1 ζ 2 Q β x 1 x 2 t + 3 2 h Q ̂ β x 1 x 2 t , E44
σ 33 = 3 4 1 3 ζ 3 ζ p 1 x 1 x 2 t + ζ p 2 x 1 x 2 t + σ 0 x 1 x 2 t , E45
μ αβ = 3 2 h 1 ζ 2 R αβ x 1 x 2 t + 3 2 h R αβ x 1 x 2 t , E46
μ β 3 = 6 h 2 ζ S β x 1 x 2 t , E47
μ 3 β = 0 , E48
μ 33 = ζV x 1 x 2 t + T x 1 x 2 t , E49

where

p x 1 x 2 t = σ t x 1 x 2 t σ b x 1 x 2 t , E50
σ 0 x 1 x 2 t = 1 2 σ t x 1 x 2 t + σ b ( x 1 x 2 t ) , E51
V x 1 x 2 t = 1 2 μ t x 1 x 2 t μ b ( x 1 x 2 t ) , E52
T x 1 x 2 t = 1 2 μ t x 1 x 2 t + μ b ( x 1 x 2 t ) . E53

The pressures p 1 and p 2 are chosen in the form

p 1 x 1 x 2 t = ηp x 1 x 2 t , E54
p 2 x 1 x 2 t = 1 η 2 p x 1 x 2 t . E55

and η R is called splitting parameter.

The three-dimensional displacements u i and microrotations ϕ i

u α = h 2 ζ Ψ α x 1 x 2 t , E56
u 3 = W x 1 x 2 t + 1 ζ 2 W x 1 x 2 t , E57
ϕ α = Ω α 0 x 1 x 2 t 1 ζ 2 + Ω ̂ α x 1 x 2 t , E58
ϕ 3 = ζ Ω 3 x 1 x 2 t , E59

and the three-dimensional strain and torsion tensors γ ji and χ ji

γ αβ = 6 h 2 ζ e αβ x 1 x 2 t , E60
γ 3 β = 3 2 h 1 ζ 2 ω β x 1 x 2 t , E61
γ β 3 = 3 2 h 1 ζ 2 ω β x 1 x 2 t + 3 2 h ω ̂ β x 1 x 2 t , E62
χ αβ = 3 2 h 1 ζ 2 τ αβ x 1 x 2 t + 3 2 h τ αβ x 1 x 2 t , E63
χ 3 β = 6 h 2 ζτ β x 1 x 2 t , E64

where ζ = 2 x 3 h .

Then zero variation of the functional

δ Θ s η = 0

is equivalent to the plate bending system of equations (A) and constitutive formulas (B) mixed problems.

A. The bending equilibrium system of equations:

M αβ , α Q β = I 1 2 Ψ β t 2 , E65
Q α , α + p ̂ 1 = I 2 2 W t 2 , E66
R αβ , α + ε 3 βγ Q γ Q γ = I αβ 2 Ω α 0 t 2 , E67
ε 3 βγ M βγ + S α , α = I 3 2 Ω 3 t 2 , E68
Q ̂ α , α + p ̂ 2 = I 2 2 W t 2 , E69
R αβ , α + ε 3 βγ Q ̂ γ = I αβ 0 2 Ω ̂ α 0 t 2 , E70

where I 1 = h 3 12 ρ , I 2 = 2 h 3 ρ , I αβ = 5 h 6 J αβ , I 3 = h 2 6 J 33 , I αβ 0 = 2 h 3 J αβ , p ̂ 1 = η opt p , and p ̂ 2 = 2 3 1 η opt p , with the resultant traction boundary conditions:

M αβ n β = Π , R αβ n β = M , E71
Q α n α = Π o 3 , S α n α = ϒ o 3 , E72

at the part Γ σ and the resultant displacement boundary conditions

Ψ α = Ψ , W = W o , Ω α 0 = Ω 0 , Ω 3 = Ω o 3 , E73

at the part Γ u .

B. Constitutive formulas in the reverse form:1

M αα = μ λ + μ h 3 3 λ + 2 μ Ψ α , α + λμ h 3 6 λ + 2 μ Ψ β , β + 3 p 1 + 5 p 2 λ h 2 30 λ + 2 μ , E74
M βα = μ α h 3 12 Ψ α , β + μ + α h 3 12 Ψ β , α + 1 α ' α h 3 6 Ω 3 , E75
R βα = 5 γ ε h 6 Ω β , α 0 + 5 γ + ε h 6 Ω α , β 0 , E76
R αα = 10 β + γ 3 β + 2 γ Ω α , α 0 + 5 hβγ 3 β + 2 γ Ω β , β 0 , E77
R βα = 2 γ ε h 3 Ω ̂ β , α + 2 γ + ε h 3 Ω ̂ α , β , E78
R αα = 8 γ γ + β h 3 β + 2 γ Ω ̂ α , α + 4 γβh 3 β + 2 γ Ω ̂ β , β , E79
Q α = 5 μ + α h 6 Ψ α + 5 μ α h 6 W , α + 2 μ α h 3 W , α + 1 β 5 3 Ω β 0 + 1 β 5 3 Ω ̂ β , E80
Q α = 5 μ α h 6 Ψ α + 5 μ α 2 h 6 μ + α W , α + 2 μ + α h 3 W , α + 1 α 5 3 Ω β 0 + μ α μ + α Ω ̂ β , E81
Q ̂ α = 8 αμh 3 μ + α W , α + 1 α 8 αμh 3 μ + α Ω ̂ β , E82
S α = 5 γε h 3 3 γ + ε Ω 3 , α , E83

and the optimal value η opt of the splitting parameter is given as in [10]

η opt = 2 W 00 W 10 W 01 2 W 11 + W 00 W 10 W 01 . E84

where

W ij = S η = i E η = j .

We also assume that the initial condition can be presented in the form

U x 1 x 2 0 = U 0 x 1 x 2 , U t x 1 x 2 0 = V 0 x 1 x 2

## 4. Cosserat plate dynamic field equations

The Cosserat plate field equations are obtained by substituting the relations Eqs. (74)(83) into the system of Eqs. (65)(70) similar to [10]:

L U = K 2 U t 2 + F η , E85

where

L = L 11 L 12 L 13 L 14 0 L 16 kL 13 0 L 16 L 12 L 22 L 23 L 24 L 16 0 kL 23 L 16 0 L 13 L 23 L 33 0 L 35 L 36 L 77 L 38 L 39 L 41 L 42 0 L 44 0 0 0 0 0 0 L 16 L 38 0 L 55 L 56 kL 35 L 58 0 L 16 0 L 39 0 L 56 L 66 kL 36 0 L 58 L 13 L 14 L 73 0 L 35 L 36 L 77 L 78 L 79 0 L 16 L 78 0 L 85 L 56 kL 35 L 88 kL 56 L 16 0 L 79 0 L 56 L 55 kL 36 kL 56 L 99 ,
K = h 3 12 ρ 0 0 0 0 0 0 0 0 0 h 3 12 ρ 0 0 0 0 0 0 0 0 0 2 h 3 ρ 0 0 0 0 0 0 0 0 0 h 2 6 J 33 0 0 0 0 0 0 0 0 0 5 h 6 J 11 5 h 6 J 12 0 0 0 0 0 0 0 5 h 6 J 12 5 h 6 J 22 0 0 0 0 0 0 0 0 0 2 h 3 ρ 0 0 0 0 0 0 0 0 0 2 h 3 J 11 2 h 3 J 12 0 0 0 0 0 0 0 2 h 3 J 12 2 h 3 J 22 ,
U = Ψ 1 Ψ 2 W Ω 3 Ω 1 0 Ω 2 0 W Ω 1 0 Ω 2 0 T ,
F η = 3 h 2 λ 3 p 1 , 1 + 5 p 2 , 1 30 λ + 2 μ , 3 h 2 λ 3 p 1 , 2 + 5 p 2 , 2 30 λ + 2 μ , p 1 , 0 , 0 , 0 , h 2 3 p 1 + 4 p 2 24 , 0 , 0 T ,
p 1 = ηp , p 2 = 1 η 2 p

The operators L ij are given as follows

L 11 = c 1 2 x 1 2 + c 2 2 x 2 2 c 3 , L 12 = c 1 c 2 2 x 1 x 2 , L 13 = c 11 x 1 , L 14 = c 12 x 2 , L 16 = c 13 , L 17 = k 1 c 11 x 1 , L 22 = c 2 2 x 1 2 + c 1 2 x 2 2 c 3 , L 23 = c 11 x 2 , L 24 = c 12 x 1 , L 33 = c 3 2 x 1 2 + 2 x 2 2 , L 35 = c 13 x 2 , L 36 = c 13 x 1 , L 38 = c 10 x 2 , L 39 = c 10 x 1 , L 41 = c 12 x 2 , L 42 = c 12 x 1 , L 44 = c 6 2 x 1 2 + 2 x 2 2 2 c 12 , L 55 = c 7 2 x 1 2 + c 8 2 x 2 2 2 c 13 , L 56 = c 7 c 8 2 x 1 x 2 , L 58 = c 9 , L 66 = c 8 2 x 1 2 + c 7 2 x 2 2 2 c 13 , L 73 = c 5 2 x 1 2 + 2 x 2 2 , L 77 = c 4 2 x 1 2 + 2 x 2 2 , L 78 = c 14 x 2 , L 79 = c 14 x 1 , L 85 = c 7 2 x 1 2 + c 8 2 x 2 2 2 c 13 , L 88 = c 7 2 x 1 2 + c 8 2 x 2 2 c 15 , L 99 = c 8 2 x 1 2 + c 7 2 x 2 2 c 15 .

The coefficients c i are given as

c 1 = h 3 μ λ + μ 3 λ + 2 μ , c 2 = h 3 α + μ 12 , c 3 = 5 h α + μ 6 , c 4 = 5 h α μ 2 6 α + μ , c 5 = h 5 α 2 + 6 αμ + 5 μ 2 6 α + μ , c 6 = h 3 γε 3 γ + ε , c 7 = 10 β + γ 3 β + 2 γ , c 8 = 5 h γ + ε 6 , c 9 = 10 h α 2 3 α + μ , c 10 = 5 α μ 3 α + μ , c 11 = 5 h α μ 6 , c 12 = h 3 α 6 , c 13 = 5 3 , c 14 = 5 α + 3 μ 3 α + μ , c 15 = 2 5 α + 4 μ 3 α + μ .

## 5. Numerical validation

For the validation purposes we provide the algorithm and computation results for the three-dimensional Cosserat elastodynamics. We also present the analysis of the numerical results based on the plate theory for the microelements of different shapes and orientations incorporated into the Cosserat plate.

### 5.1 Analysis of Cosserat plate vibrations based on the three-dimensional theory

In our computations we consider the plates made of polyurethane foam—a material reported in the literature to behave Cosserat like—and the values of the technical elastic parameters presented in [15]: E = 299.5 MPa , ν = 0.44 , l t = 0.62 mm , l b = 0.327 mm , N 2 = 0.04 . Taking into account that the ratio β / γ is equal to 1 for bending [15], these values of the technical constants correspond to the following values of Lamé and Cosserat parameters: λ = 762.616 MPa , μ = 103.993 MPa , α = 4.333 MPa , β = 39.975 MPa , γ = 39.975 MPa , ε = 4.505 MPa . We consider a low-density rigid foam usually characterized by the densities of 24–50 kg/m3 [16]. In all further numerical computations we used the density value ρ = 34 kg/m3 and different values the rotatory inertia J .

Let us consider the plate B 0 being a rectangular cuboid 0 a × 0 , a ] × h 2 h 2 . Let the sets T and B be the top and the bottom surfaces contained in the planes x 3 = h 2 and x 3 = h 2 respectively, and the curve Γ = Γ 1 Γ 2 be the lateral part of the boundary:

Γ 1 = x 1 x 2 x 3 : x 1 0 a x 2 0 a x 3 h 2 h 2 , Γ 2 = x 1 x 2 x 3 : x 1 0 a x 2 0 a x 3 h 2 h 2 ,

We solve the three-dimensional Cosserat equilibrium Eqs. (1)(2) accompanied by the constitutive Eqs. (3)(4) and strain-displacement and torsion-rotation relations Eq. (5) complemented by the following boundary conditions:

Γ 1 : u 2 = 0 , u 3 = 0 , φ 1 = 0 , σ 11 = 0 , μ 12 = 0 , μ 13 = 0 ; E86
Γ 2 : u 1 = 0 , u 3 = 0 , φ 2 = 0 , σ 22 = 0 , μ 21 = 0 , μ 23 = 0 ; E87
T : σ 33 = p x 1 x 2 , μ 33 = 0 ; E88
B : σ 33 = 0 , μ 33 = 0 . E89

where the initial distribution of the pressure is given as p = sin π x 1 a sin π x 2 a sin ωt and the rotatory inertia tensor J is assumed to have a diagonal form

J = J x 0 0 0 J y 0 0 0 J z . E90

Using the method of separation of variables and taking into account the boundary conditions Eqs. (86)(87), we express the kinematic variables in the form:

u 1 = cos π x 1 a sin π x 2 a z 1 x 3 sin ωt , E91
u 2 = sin π x 1 a cos π x 2 a z 2 x 3 sin ωt , E92
u 3 = sin π x 1 a sin π x 2 a z 3 x 3 sin ωt , E93
ϕ 1 = sin π x 1 a cos π x 2 a z 4 x 3 sin ωt , E94
ϕ 2 = cos π x 1 a sin π x 2 a z 5 x 3 sin ωt , E95
ϕ 3 = cos π x 1 a cos π x 2 a z 6 x 3 sin ωt , E96

where the functions z i x 3 represent the transverse variations of the kinematic variables.

If we substitute the expressions Eqs. (91)(96) into Eqs. (3)(4) and then into Eqs. (1)(2), we will obtain the following eigenvalue problem

B z = ω 2 A z E97

where

B = b 1 L 2 + b 2 L 0 b 3 L 0 b 4 L 1 0 b 5 L 1 b 6 L 0 b 3 L 0 b 1 L 2 + b 2 L 0 b 4 L 1 b 5 L 1 0 b 6 L 0 b 4 L 1 b 4 L 1 b 7 L 2 b 6 L 0 b 6 L 0 0 0 b 5 L 1 b 6 L 0 b 9 L 2 + b 10 L 0 b 11 L 0 b 12 L 1 b 5 L 1 0 b 6 L 0 b 11 L 0 b 9 L 2 + b 10 L 0 b 12 L 1 b 6 L 0 b 6 L 0 0 b 12 L 1 b 12 L 1 b 13 L 2 + b 2 L 14 , E98
A = a 2 ρ 0 0 0 0 0 0 a 2 ρ 0 0 0 0 0 0 a 2 ρ 0 0 0 0 0 0 a 2 J x 0 0 0 0 0 0 a 2 J y 0 0 0 0 0 0 a 2 J z , E99
z = z 1 , z 2 , z 3 , z 4 , z 5 , z 6 T , E100

and the differential operators L i are defined as

L 0 = I , L 1 = d dx 3 , L 2 = d 2 dx 3 2

and the coefficients b i are defined as

b 1 = a 2 μ + α , b 2 = π 2 α + λ + 3 μ , b 3 = π 2 λ + μ α , b 4 = λ + μ α , b 5 = 2 a 2 α , b 6 = 2 aπα , b 7 = a 2 2 μ + λ , b 8 = 2 π 2 α + μ , b 9 = a 2 γ + ε , b 10 = π 2 β + ε + 3 γ , b 11 = π 2 β + γ ε , b 12 = β + γ ε , b 13 = a 2 β + 2 γ , b 14 = 2 π 2 γ + ε 4 a 2 α

The system of differential Eq. (97) is complemented by the following boundary conditions D z = D 0 for x 3 = h 2 and D z = 0 for x 3 = h 2 .

D = d 1 L 1 0 d 2 L 0 0 d 3 L 0 0 0 d 1 L 1 d 2 L 0 d 3 L 0 0 0 d 4 L 0 d 4 L 0 d 5 L 1 0 0 0 0 0 0 d 6 L 1 0 d 7 L 0 0 0 0 0 d 6 L 1 d 7 L 0 0 0 0 d 8 L 0 d 8 L 0 d 9 L 1 , E101
D 0 = 0 , 0 , a , 0 , 0 , 0 T , E102

and the coefficients d i are defined as

d 1 = a μ + α , d 2 = π μ α , d 3 = 2 , d 4 = a λ + 2 μ , d 5 = πλ , d 6 = a γ + ε , d 7 = a γ ε , d 8 = πβ , d 9 = a β + 2 γ .

The idea for the solution of the eigenvalue problem Eq. (97) is based on the following algorithm:

Step 1. Fix certain frequency value.

We fix certain value of the frequency ω and force the Cosserat body to vibrate at this frequency.

Step 2. Solve the three-dimensional Cosserat system of equations.

Mathematically, fixing certain value of ω implies that three-dimensional system of Eq. (97) has a constant right-hand side and therefore can be solved for the kinematic variables as a static system of equations. We solve the system Eq. (97) using the high-precision Runge-Kutta method incorporated in Mathematica software similar to how it was done in [7].

Step 3. Find large amplitudes of the kinematic variables.

We run ω through an interval of positive real values and take note where the solution changes its sign and the amplitude of the solutions starts to grow indefinitely. This corresponds to the oscillation of the Cosserat body at its resonant frequency. Thus, when the frequency ω coincides with the natural frequency of the plate the resonance will occur and the large amplitude linear vibrations can be observed ( Figure 1 ).

The comparison of the eigenfrequencies of the Cosserat plate with the eigenfrequencies of the three-dimensional Cosserat elasticity is given in the Table 1 . The rotatory inertia principle moments used are J x = 0.001 , J y = 0.001 , J z = 0.001 , which represent a ball-shaped microelement ( Figure 2 ). The relative error of the natural macro frequencies associated with the rotation of the middle plane and the flexural motion is less than 1%.

ω 1 , ω 2 ω 3 , ω 7 ω 4 ω 5 , ω 8 ω 6 , ω 9
Plate theory 0.310 17.881 501.13 205.62 338.95
D Cosserat elasticity 0.309 17.763 530.82 211.98 317.87

### Table 1.

Comparison of the eigenfrequencies ω i (Hz) with the exact values of the 3D Cosserat elasticity.

### 5.2 Analysis of Cosserat plate vibrations based on the plate theory

We consider a plate a × a of thickness h with the boundary G = G 1 G 2

G 1 = x 1 x 2 : x 1 0 a x 2 0 a
G 2 = x 1 x 2 : x 2 0 a x 1 0 a

and the following hard simply supported boundary conditions [7]:

G 1 : W = 0 , W = 0 , Ψ 2 = 0 , Ω 1 0 = 0 , Ω ̂ 1 0 = 0 , Ω 3 = 0 , Ψ 1 n = 0 , Ω 2 0 n = 0 , Ω ̂ 2 0 n = 0 ;
G 2 : W = 0 , W = 0 , Ψ 1 = 0 , Ω 2 0 = 0 , Ω ̂ 2 0 = 0 , Ω 3 = 0 , Ψ 2 n = 0 , Ω 1 0 n = 0 , Ω ̂ 1 0 n = 0 .

Similar to [12] we apply the method of separation of variables for the eigenvalue problem Eq. (85) to solve for the kinematic variables Ψ α , W , Ω 3 , Ω α 0 , W and Ω α 0 . The kinematic variables can be further expressed in the following form

Ψ 1 nm = A 1 cos x 1 a sin x 2 a sin ωt + B 1 sin x 1 a cos x 2 a sin ωt , Ψ 2 nm = A 2 sin x 1 a cos x 2 a sin ωt + B 2 cos x 1 a sin x 2 a sin ωt , W nm = A 3 sin x 1 a sin x 2 a sin ωt + B 3 cos x 1 a cos x 2 a sin ωt , Ω 3 nm = A 4 cos x 1 a cos x 2 a sin ωt + B 4 sin x 1 a sin x 2 a sin ωt , Ω 1 0 , nm = A 5 sin x 1 a cos x 2 a sin ωt + B 5 cos x 1 a sin x 2 a sin ωt , Ω 2 0 , nm = A 6 cos x 1 a sin x 2 a sin ωt + B 6 sin x 1 a cos x 2 a sin ωt , W , nm = A 7 sin x 1 a sin x 2 a sin ωt + B 7 cos x 1 a cos x 2 a sin ωt , Ω ̂ 1 0 , nm = A 8 sin x 1 a cos x 2 a sin ωt + B 8 cos x 1 a sin x 2 a sin ωt , Ω ̂ 2 0 , nm = A 9 cos x 1 a sin x 2 a sin ωt + B 9 sin x 1 a cos x 2 a sin ωt ,

where A i and B i are constants.

We solve an eigenvalue problem by substituting these expressions into the system of Eq. (85). The obtained nine sequences of positive eigenfrequencies ω i nm are associated with the rotation of the middle plane ( ω 1 nm and ω 2 nm ), flexural motion and its transverse variation ( ω 3 nm and ω 7 nm ), micro rotatory inertia ( ω 4 nm , ω 5 nm and ω 6 nm ) and its transverse variation ( ω 8 nm and ω 9 nm ) [12].

We perform all our numerical simulations for a = 3.0 m and h = 0.1 m. We consider different forms of micro elements: ball-shaped elements, horizontally and vertically stretched ellipsoids (see Figure 2 ). For simplicity we will use the notation ω i for the first elements ω i 11 of the sequences ω i nm . The results of the computations are given in the Table 2 . The shape of the micro-elements does not effect the natural macro frequencies ω 1 and ω 2 associated with the rotation of the middle plane and ω 3 and ω 7 associated with the flexural motion and its transverse variation. The ellipsoid elements have higher micro frequencies associated with the micro rotatory inertia ( ω 4 , ω 5 and ω 6 ) and its transverse variation ( ω 8 and ω 9 ), than the ball-shaped elements.

ω 1 , ω 2 ω 3 , ω 7 ω 4 ω 5 , ω 8 ω 6 , ω 9
Shape J x J y J z ω 1 , ω 2 ω 3 , ω 7 ω 4 ω 5 , ω 8 ω 6 , ω 9
Ball 0.001 0.001 0.001 17.88 0.31 501.13 205.62 338.95
Vertical ellipsoid 0.001 0.001 0.0001 17.88 0.31 501.13 650.22 338.95
Horizontal ellipsoid 0.0001 0.001 0.001 17.88 0.31 1363.01 205.62 394.08

### Table 2.

Eigenfrequencies ω i 11 (Hz) for different shapes of micro-elements.

Let J x , J y and J z be the principal moments of inertia of the microelements corresponding to the principal axes of their rotation. We assume that the quantities J x , J y and J z are constant throughout the plate B 0 . If the microelements are rotated around the z -axis by the angle θ the rotatory inertia tensor J can be expressed as

J = J x cos 2 θ + J y sin 2 θ J x J y sin 2 θ 0 J x J y sin 2 θ J x sin 2 θ + J y cos 2 θ 0 0 0 J z E103

The eigenfrequencies for different angles of microrotation of the microelements are given in the Table 3 and the Figure 3 . The rotatory inertia principle moments used are J x = 0.002 , J y = 0.001 , J z = 0.0001 , which represent a horizontally stretched ellipsoid microelement. The case when the microelements are not aligned with the edges of the plate the model predicts some additional natural frequencies related with the microstructure of the material.

Angle θ ω 1 ω 2 ω 3 ω 7 ω 4 ω 5 ω 8 ω 6 ω 9
0 17.88 17.88 0.31 0.31 650.221 265.37 265.37 450.61 450.61
10 17.88 17.88 0.31 0.31 650.221 255.59 279.40 429.89 469.93
20 17.88 17.88 0.31 0.31 650.221 247.75 295.33 406.70 484.79
30 17.88 17.88 0.31 0.31 650.221 242.57 313.65 382.94 495.14
40 17.88 17.88 0.31 0.31 650.221 239.99 333.10 360.57 500.46
45 17.88 17.88 0.31 0.31 650.221 239.68 338.95 354.35 501.13
50 17.88 17.88 0.31 0.31 650.221 239.99 333.10 360.57 500.46
60 17.88 17.88 0.31 0.31 650.221 242.57 313.65 382.94 495.14
70 17.88 17.88 0.31 0.31 650.221 247.75 295.33 406.70 484.79
80 17.88 17.88 0.31 0.31 650.221 255.59 279.40 429.89 469.93
90 17.88 17.88 0.31 0.31 650.221 265.37 265.37 450.61 450.61

### Table 3.

Eigenfrequencies ω i 11 (Hz) for different angles of rotation of horizontal ellipsoid micro-elements.

## 6. Conclusions

In this chapter, we presented a mathematical model of Cosserat plate vibrations. The dynamic model of the plates has been developed as a dynamic extension of the Reissner plate theory. The equations has been presented in both tensorial and the matrix forms. We also described the validation of the model, which is based on the comparison with the three-dimensional Cosserat elastodynamics exact solutions. Based on the presented results of the computer simulations we were able to detect and classify the additional high resonance frequencies of a plate. We have shown that the frequencies depend on the shape and orientation of microelements (ball-shaped elements, horizontally and vertically stretched ellipsoids) incorporated into the Cosserat plates. We also have been able to identify that micro frequencies associated with the micro rotatory inertia and its transverse variation of the ellipsoid elements have higher micro frequencies than the ball-shaped elements. We also showed the dependence of the eigenfrequencies on the angles of rotation of the horizontal ellipsoid micro-elements. These results can be used to identify the characteristics of the plate micro-elements.

## A.1 Conventions

We use the following notation convention:

1. the values of the Latin subindex i take values in the set 1,2,3

2. the values of the Greek indices α and β take values in the set 1 2

3. the Einstein summation notation is used throughout the chapter

## A.2 Notations

x i

artesian coordinates

P

Cosserat thin plate

h

plate thickness

μ , λ

Lamé parameters

α , β , γ , ε

Cosserat elasticity parameters

ρ

material density

J ji or J

rotatory inertia

σ ji or σ

the stress tensor

μ ji or μ

the couple stress tensor

γ ji or γ

strain tensor

χ ji or χ

bend-twist tensor

u i or u

displacement vector

ϕ i or ϕ

microrotation vector

p i or p

linear momentum

q i or q

angular momentum

ε ijk

Levi-Civita tensor

U C

strain stored energy

U K

stress energy

T C

stored kinetic energy

T W

work of inertia forces

S

Cosserat plate stress set

U

Cosserat plate displacement set

E

Cosserat plate strain set

η

splitting parameter

p

pressure

ω

natural frequency of plate vibration

θ

angle of microelement orientation

M αβ

bending and twisting moments

Q α

shear forces

Q α , Q ̂ α

transverse shear forces

R αβ

micropolar bending moments

R αβ

micropolar twisting moments

S α

micropolar couple moments

Ψ α

rotations of the middle plane around x α axis

W , W

vertical deflections of the middle plate

Ω α 0

microrotations in the middle plate around x α axis

Ω 3

rate of change of the microrotation

## References

1. 1. Cosserat E, Cosserat F. Theorie des corps deformables [Theory of Deformable Bodies]. Paris: A. Hermann et fils; 1909
2. 2. Green A, Naghdi P. The linear theory of an elastic Cosserat plate. Mathematical Proceedings of the Cambridge Philosophical Society. 1966;63:537-550
3. 3. Eringen A. Theory of micropolar plates. Journal of Applied Mathematics and Physics. 1967;18:12-31
4. 4. Altenbach H, Eremeyev V. On the theories of plates based on the Cosserat approach. Advances in Applied Mathematics and Mechanics. 2010;21:73-92
5. 5. Steinberg L. Deformation of micropolar plates of moderate thickness. International Journal of Applied Mathematics and Mechanics. 2010;6(17):1-24
6. 6. Kvasov R, Steinberg L. Numerical modeling of bending of Cosserat elastic plates. In: Proceedings of the 5th Computing Alliance of Hispanic-Serving Institutions; 2011. pp. 67-70
7. 7. Steinberg L, Kvasov R. Enhanced mathematical model for Cosserat plate bending. Thin-Walled Structures. 2013;63:51-62
8. 8. Reissner E. On the theory of elastic plates. Journal of Mathematics and Physics. 1944;23:184-191
9. 9. Reissner E. The effect of transverse shear deformation on the bending of elastic plates. Journal of Applied Mechanics. 1945:69-77
10. 10. Kvasov R, Steinberg L. Numerical modeling of bending of micropolar plates. Thin-Walled Structures. 2013;69:67-78
11. 11. Kvasov R, Steinberg L. Modeling of size effects in bending of perforated Cosserat plates. Modelling and Simulation in Engineering. 2017;2017:1-19
12. 12. Steinberg L, Kvasov R. Analytical modeling of vibration of micropolar plates. Applied Mathematics. 2015;6:817-836
13. 13. Nowacki W. Theory of Asymmetric Elasticity. Oxford, New York: Pergamon Press; 1986
14. 14. Gurtin M. The linear theory of elasticity. In: Truesdell C, editor. Handbuch der Physik. Vol. VIa/2. Berlin: Springer-Verlag; 1972. pp. 1-296
15. 15. Lakes R. Experimental methods for study of Cosserat elastic solids and other generalized elastic continua. In: Mühlhaus H, editor. Continuum Models for Materials with Microstructures. New York: Wiley J; 1995. pp. 1-22
16. 16. Singh S. Blowing agents for polyurethane foams. Rapra Review Report. 2002;12(10)

## Notes

• In the following formulas a subindex β = 1 if α = 2 and β = 2 if α = 1.

Written By

Lev Steinberg and Roman Kvasov

Submitted: May 1st, 2019 Reviewed: May 24th, 2019 Published: July 10th, 2019