Reduced Linear Constrained Elastic and Viscoelastic Homogeneous Cosserat Media as Acoustic Metamaterials

: We consider the reduced constrained linear Cosserat continuum, a particular type of a Cosserat medium, for three different material behaviors or symmetries: the isotropic elastic case, a special type of elastic transversely isotropic case, and the isotropic viscoelastic case. Such continua, in which stresses do not work on rates of microrotation gradients, behave as acoustic metamaterials for the (pure) shear waves and also for one branch of the mixed wave in the considered anisotropic material case. In elastic media, those waves do not propagate for frequencies exceeding a certain threshold, whence these media exhibit a single negative acoustic metamaterial behavior in this range. In the isotropic viscoelastic case, dissipation destroys the bandgap and favors wave propagation. This curious effect is, probably, due to the fact that the bandgap is associated not with the dissipation, but with the wave localization which can be destroyed by the viscosity. The dispersion curve is now decreasing in some part of the former bandgap, above a certain frequency, whence the medium is a double negative acoustic metamaterial. We prove the existence of a boundary wavenumber in the viscoelastic case and estimate its value. Below the characteristic frequency corresponding to the boundary of the elastic bandgap, the wave attenuation (logarithmic decrement) is a growing function of the viscous dissipation parameter. Above this frequency, the attenuation decreases as the viscosity increases. of linear constrained reduced Cosserat media: elastic isotropic,


Introduction
Three-dimensional solid continua, in which both translational and rotational motions of body points are considered, and respectively characterized by displacement field u and microrotation vector field θ, were introduced by the Cosserat brothers [1]. Their theory then sank into oblivion for about 50 years before it was systematically revisited by [2]. Among recent books on the subject, one can mention [3].
The present study focusses on acoustic or seismic wave propagation in reduced Cosserat continua, a special class of Cosserat media in which internal stresses do not work on the microrotation gradient velocity. This entails special symmetry conditions on the strain energy, which depends on the Cosserat deformation tensor but not on the wryness tensor. Possible strain measures for continua with rotational degrees of freedom with various types of symmetry are discussed in [4][5][6]. Such models may be applied to granular materials [7][8][9][10], media with high contrasts, geomaterials [11][12][13][14][15][16], in which non-classical wave propagation is observed [17]. One of the first applications of an elastic reduced isotropic Cosserat model to wave propagation in granular materials was suggested in [7].
Media with effective negative parameters (density and/or the relevant elastic modulus), as regards wave propagation, are called acoustic metamaterials [18]. This phenomenon takes place in a certain frequency domain. If both elastic and inertial effective frequency-dependent parameters are negative, one has a double negative acoustic metamaterial. It exhibits negative refraction, i.e., a decreasing dispersion curve. If only one of the effective frequency-dependent parameters is negative, this is a single negative acoustic metamaterial, i.e., the dispersion graph for the corresponding harmonic wave has a band gap [19]. These media can be used, for instance, to control wave beams, energy fluxes, or to conceal objects, rendered acoustically undetectable.
We can consider more general elastic continua, in which the strain energy depends on a certain vector field but not on its gradient. In reduced Cosserat media this is the microrotation vector. One can intuitively view a reduced complex continuum as a hosting continuum, with all elastic connections present, including isolated resonators that are not directly coupled to one another.
Such reduced generalized elastic media were shown to behave as acoustic metamaterials [20]. In particular, reduced elastic Cosserat media are single negative acoustic metamaterials, for which harmonic shear waves are evanescent in a certain frequency domain. Lots of works consider acoustic metamaterials with internal (often rotational) degrees of freedom, on the discrete and continuum level [21][22][23][24][25][26][27].
The possible interpretation of a time-dispersive material as an "extended" material with hidden internal degrees of freedom is suggested, for instance, in [28]. Viscosity may drastically change the wave properties of an elastic metamaterial [13], destroying the bandgap and, in some cases, leading to a decreasing dispersion curve in some range. Relaxed linear micromorphic metamaterials (with band gaps) are discussed in [29,30], where a certain kinematical constraint is accepted, as well as in the full constrained Cosserat medium [31]. Other materials with microstructure having band gaps are investigated in [32,33], the latter one being an example of a hosting continuum enhanced by distributed resonators, discussed in [20]. In linear constrained Cosserat media the microrotation θ is identically equal to the vortex deformation (macrorotation) of the neighbourhood ∇ × u/2. These models are popular in the scientific community [34,35], possibly they could be used for the description of blocky media.
Specialists in the rotational seismology, a science which has developed during the past decades, agree that rotational degrees of freedom of blocks and heterogeneities in the geomedium must be taken into account [36][37][38]. However, there is a strong discussion on the most appropriate model. There are works that insist taking into account only the macro-rotation of the medium, see, for instance, [39]. To satisfy the principle of material objectivity, in this case, it is necessary to accept the kinematical hypothesis that the macro-rotation is equal to the micro-rotation, i.e., to consider the constrained Cosserat theory. On the other hand, the presence of slow waves [40] is an argument in favor of reduced Cosserat continuum or another acoustic metamaterial as an appropriate model.
In this work, we have considered elastic and viscoelastic models which combine features of the constrained Cosserat medium and reduced Cosserat medium. The viscoelastic model is taken in its simplest isotropic variant, which can be found in [41] for the full linear Cosserat medium; we consider a particular case, when it is reduced (i.e., couple stresses are zero) and constrained (i.e., a kinematical hypothesis "microrotation equals macrorotation" is accepted). The investigation of wave propagation in such media (elastic isotropic, elastic with anisotropic coupling and viscoelastic isotropic) will give some information for experimentalists to judge if it is an adequate theory for a real granular or blocky medium. We have chosen this model also since it has a very few parameters and therefore we can clearly see in this case how the viscosity can change the wave propagation in a rotational metamaterial.

Model and Basic Laws
We consider a linear micropolar medium each point of which is characterized by an infinitesimal translational displacement u and an infinitesimal microrotation vector θ. In the linear medium, actual and initial configurations coincide as well as nabla operators in both of them. Since we consider a reduced Cosserat medium, we suppose that no stresses work on ∇θ, see [7,42]. Since we consider a constrained Cosserat medium, we accept a kinematical hypothesis as in [31]. Each body point of the medium possesses density ρ and density of the tensor of inertia ρI. We will shortly recall here dynamic laws for a constrained reduced Cosserat medium of any type, which can be found in [13], and write down its final form for the general case. Balance of force in its local form looks as follows: where τ is the stress tensor, ρK is the density of the external body force, ()˙is the time derivative. We can write this law, representing τ as Here τ S = (τ mn + τ nm )i m i n /2 and τ A = (τ mn − τ nm )i m i n /2 are the symmetric and antisymmetric parts of τ, respectively, and τ × def = τ mn i m × i n , i m being a certain basis, E is the identity tensor. Then (2) takes form Balance of moment in its local form is where ρL is the density of the external body moment. Substituting the kinematical hypothesis (1) into the balance of moment, we obtain Then we substitute this expression into the balance of force (4) and obtain Note that if τ S obeys the equations of a local theory, we still obtain a local theory, although with a sophisticated dynamical term that has mixed spatial-temporal derivatives. This is contrary to the case of the constrained full Cosserat medium, where the kinematical constraint yields in equivalent second gradient theory. This difference is related to the fact that the strain energy in the reduced model does not depend on ∇θ.
The law of balance of energy looks as (heat transfer is not considered). Equation (7) is valid for any type of the linear constrained reduced Cosserat medium. Symmetric part of the stress τ S is determined by the constitutive equations.
We see that antisymmetric part of τ is determined by balance of moment (6) and not by a constitutive equation. This is due to the fact that in the linear non-constrained model τ × works on the strain rate ∇ ×u/2 −θ. Indeed, if U nc is the mass density of the strain energy for a linear non-constrained reduced model, then [42] Note that Constitutive equations for the non-constrained elastic medium are obtained from the law of balance of energy both for the symmetric and antisymmetric part of τ. In the constrained model, θ − ∇ ×u/2 ≡ 0, therefore τ A cannot be obtained from this law, but it is determined by the dynamic law (6).

Constitutive Equations for the Linear Elastic Constrained Reduced Cosserat Media
We will consider the elastic medium in adiabatic or isothermal processes. The second law of thermodynamics is satisfied identically for this case. Since the medium is linear, U a quadratic form of ∇u S : the same as for the classical medium, where C = C mnkl i m i n i k i l is the elastic tensor of the fourth rank. As in the classical elasticity, it satisfies the symmetry conditions C mnkl = C nmkl = C klmn . The law of balance of energy (8) thus gives us the constitutive equation However, since we have rotational degrees of freedom and the kinematical constraint θ = ∇ × u/2, the medium is not classical, and τ A = 0 but determined from dynamic Equation (6).
Concrete form of C is determined by the symmetry of the material: its group of symmetry includes the group of symmetry of the medium, according to the Curie-Neumann's principle [43]. We can also obtain the same equations in terms of modified strain energy with Lagrange multipliers. Consider zero external loads. Indeed, Lagrange equations for the continuum look as where q is a generalised vector consisting of generalised co-ordinates (components of u and θ in our case), Lagrangian density L = ρ(U − K), ρK being the density of the kinetic energy, ρU the density of the strain energy of the continuum. This is true in case when there are no kinematical constraints. When they exist, we have to use the density of the modified strain energy U m instead of U. In our (linear) case U is defined by (11), and where Λ is a Lagrange multiplier, an antisymmetric tensor of the second rank, vector λ = −Λ × . In fact, Λ has a physical sense of τ A .

Note that
This allows us to write down Equation (15) where Λ = Λ A . We see that up to the notation Λ = τ A these equations coincide with balance of force (2), balance of moment (5) for zero loads provided the constitutive relation (12).

Isotropic Elastic Model
First we choose an isotropic homogeneous elastic model. Then the elastic tensor C is also isotropic and looks as where λ, µ are elastic Lamé constants. Then the constitutive equation for the symmetric part of the stress tensor takes form Substituting this expression to the dynamic law (7), and considering zero external loads and spherical density of the inertia tensor I = IE, we obtain equations in displacements: Let us look for the plane wave solution: u = u 0 e i(k·r+ωt) . Then we obtain Denotek = k/k, where k = |k|. We can rewrite this equation as Compression plane wave is separated and it is not influenced by the rotational degrees of freedom. Indeed, if u 0 = u 0k , we obtain the classical dispersion relation for the longitudinal wave: Denote ω 2 1 = 4µ I (this will be a certain characteristic frequency of the medium), c 2 s = µ ρ (wave velocity). The dispersion relation for the shear-rotational wave (u 0 ·k = 0), given by (26), is We see that at low frequency this dispersion branch starts as the classical shear wave, but for higher frequencies dispersion becomes more pronounced, and the curve has a horizontal asymptote ω = ω 1 (boundary frequency). Above this frequency, no shear wave propagates (see Figure 1).
In the reduced non-constrained linear isotropic elastic Cosserat medium (see [42]) we do not accept kinematical hypothesis (1). Instead of that there exist a rotational elastic constant α, and the strain energy has an additional term 1 2 α(∇ × u/2 − θ) 2 . In this case there are two branches of this curve with a bandgap between them. The low branch for the non-constrained model is analogous to the curve presented in Figure 1, but instead of ω 1 for the constrained medium we have a boundary frequency ω 1 = 2(α/(I(1 + α/µ))) 1/2 . The upper branch for the non-constrained medium starts at a cut-off frequency ω 0 = 2(α/I) 1/2 . When the resistance to the strain θ − ∇ × u/2 tends to infinity (α → ∞), the cut-off frequency tends to infinity also, and the boundary frequency tends to ω 1 = 2(µ/I) 1/2 , i.e., we obtain the dispersion graph for the constrained linear isotropic elastic reduced Cosserat medium. This is to expect. This also tells us that the theory under consideration could be possibly applicable for blocky media.

Elastic Model with the Simplest Anisotropic Coupling Term (Axial Symmetry with Axis n)
Now consider the influence of anisotropy on the wave propagation. The anisotropic term of most interest would be a coupling between shear and compressional deformations, a feature which is present in granular-like materials. Such a coupling does not exist in isotropic theories. We consider a simplest variant of such an anisotropic medium: a transversally isotropic medium with the unit vector of the axis of isotropy n, where the anisotropy is only present in this coupling and has a particular form, coupling volumetric and shear deformations. Namely, let the density of the strain energy be where ρU i is the strain energy for the isotropic case. The constitutive equations are Substituting them into dynamic law (7), considering zero external loads K = 0, L = 0 and spherical symmetry for the tensor of inertia (I = IE), we obtain dynamic equations in displacements: Looking for the plane wave propagation u = u 0 e i(ωt+k·r) , we obtain the spectral problem: Denotek = k/k, then n = n ·kk +ñ, whereñ def = (E −kk) · n. Then the spectral problem can be rewritten as: Due to the specific type of coupling there exist a purely shear wave, the same as in the isotropic case, if u || (n ×k). The last and the first term in (33) then vanish and the corresponding dispersion relation the same form (28) as in the isotropic case. However, the compression wave, generally speaking, does not exist, it is coupled to the shear wave, with the exception of some special directions of wave propagation.

Special Directions of Wave Propagation
Consider waves propagating along or orthogonal to the axis of isotropy n. In both cases the last term in parentheses in (33) disappears, since eitherñ = 0 ifk||n, ork · n = 0. The spectral problem for compression and shear wave can be decoupled, as in the isotropic case, and we have the same dispersion relations, with the only difference that for k||n the wave velocity for the compressional wave is influenced by the anisotropic constant, the dispersion relation looks as Thus, if n ·k = ±1 or n ·k = 0, anisotropic term for the shear wave disappears, and its dispersion relation is the same as for the isotropic medium, longitudinal wave is non-dispersive, with constant velocity c p ,

Mixed Wave
If neither n ·k = ±1 nor n ·k = 0, longitudinal wave and shear wave with u not parallel ton ×k do not exist. Waves become mixed. Separating the "isotropic-like" shear wave, we obtain from (33) the dispersion relation for the mixed wave: where It can be rewritten as We can find eigenvectors for the mixed wave from this relation and Equation (33). Omitting some calculations, we give the result: An example of dispersion curves is given in Figure 2. Both curves are dispersive but for the upper one we hardly notice this.
The branch with "+" of Equation (38) at large k tends to c p k, and at small k These two velocities are relatively close if we require the strain energy to be positively definite. For instance, for the parameters shown in Figure 2, if we take c s = 1, we obtain the velocity at low frequencies approximately equal to 2.92, and at large frequencies approximately 2.83. The lower ("−") branch of Equation (38) has a horizontal asymptote ω = ω 1a < ω 1 , analogously to the isotropic case, where ω 2 1a = (1 − c 4 n (k ·n) 2 |k ×n| 2 /(c 2 p c 2 s ))ω 1 2 (for instance, ω 1a ≈ 0.709ω 1 for parameters in Figure 2), and it is highly dispersive. At low frequencies it behaves as and its velocity is less than c s .

Summary for the Elastic Models
We see that for the linear elastic constrained reduced Cosserat media (isotropic and special transversally isotropic with a coupling between volumetric and shear deformations) • there is an infinite bandgap above the boundary frequency ω 1 for the shear waves in both models and above a lower frequency ω 1a , which depends also on the anisotropy and propagation direction, for the slower mixed wave in the considered anisotropic case. • compression wave in the isotropic case is classical, non-dispersive • the upper branch of the mixed wave in the considered anisotropic case is weakly dispersive and has no band gap; its effective moduli, in particular, depend on the direction of the wave propagation and anisotropy parameter

Isotropic Viscoelastic Model
Let us consider how viscosity in case of the simplest viscoelastic model will influence the wave propagation. For the reduced isotropic Cosserat medium it yields in drastic change [13]. As we discussed before, for any constrained linear reduced Cosserat medium τ S is determined from the constitutive equations. We choose an isotropic homogeneous viscoelastic Kelvin-Voight model This constitutive equation satisfies the 2nd law of thermodynamics if κ, ν 0. Consider zero external body loads and spherical tensor of inertia, as in previous sections. This equation together with the balance of force (7) lets us to write down equations in displacements:

Plane Shear Wave
Taking the divergence of (43), we see that as in the purely elastic isotropic case, considered in Section 2.2, compression wave is not influenced by rotational degrees of freedom. It is separated from the shear wave. We will not discuss it here, since it is the same as in the classical viscoelasticity.
We look for the plane wave solution of (43): u = u 0 e i(k·r+ωt) and consider the case of the shear wave (u · k = 0). Shear wave obeys the following equation in the Fourier domain: It gives us the following dispersion relation: Let us investigate how the viscosity will influence the band gap which exists in the elastic case. Let us denote c s Taking the real and imaginary part of the dispersion relation (45), we obtain There is one positive root for K 2 , since n > 0: Figure 3 shows us an example of numerical graphs for the dispersion relation: In this figure we see that both k and | k| first increase with ω, then they reach their maximal values (not at the same frequency), and after this decrease with ω. This is a quite untypical behavior. For all parameters checked numerically k(ω) has one maximum, and | k|(ω) has one maximum for n which is small enough. Approximately at n > 1.5 the extremum for the imaginary part of the wavenumber disappears.

Asymptotics for Various Domains of Frequencies
Let us consider the properties of the dispersion relation (49). First, we can see from this equation, that in the viscoelastic reduced linear isotropic constrained Cosserat continuum there are no band gaps. Indeed, there is one positive root for k 2 , given (up to multiplication by a constant ω 1 /c s ) by Equation (49) for any frequencies. Consider its asymptotics at low frequencies.
At Ω → 0 we obtain almost classical shear wave. Indeed, (45) in this limit gives us In the vicinity of the boundary frequency ω = ω 1 (Ω = 1) So we see that if n < 2 dispersion curve decreases at Ω = 1. The decreasing part of the dispersion curve corresponds to the behavior of a double negative acoustic metamaterial. At n > 2 it increases at this point, and at n = 2 it has a maximum Ω = 1.
At large frequencies our medium is always a double negative acoustic metamaterial, since the dispersion curve decreases: if Ω → ∞ (Ω >> n), then We see that Figure 3 confirms this asymptotics.

Detailed Analysis of the Real Part of the Wavenumber. Influence of the Dissipation Parameter n on Its Behaviour
We see that k(ω) is a differentiable function, increasing at small ω and decreasing at large ω. Therefore, it has at least one maximum and it is limited by a certain value. We will estimate this value from above. Note that Estimating from above each term in the right part of this equation (see Appendix A), we obtain In all numerical examples K increases with Ω until a certain value Ω max where it has a maximum and then decreases, tending to 0 at large Ω, as it follows from the analytics obtained above. Numerical examples of the dispersion curve for the real part of the wave number (Figures 4-7) show us how it changes when the dissipation parameter n increases. We see that K(Ω max ), corresponding to the minimal wavelength, decreases with n. The medium exhibits features of a double negative acoustic metamaterial above Ω max , which can be controlled via the dissipation. However, Ω max is not a monotone function of n. Numerical examples show us that first, it decreases from 1 as n increases, and then increases, overpassing 1 again at n = 2. We can analytically find how the value of k at a given ω changes if we change the dissipation parameter n. If we differentiate Equation (54) with respect to n 2 and find its extremuma, we obtain that • at 0 < ω < ω 1 , i.e., in the zone where the shear wave in the elastic medium propagates, as well as at ω = ω 1 , the value of k decreases when n increases; • at ω > ω 1 , i.e., in the zone where the shear wave in the elastic medium does not exist, k for a given ω first increases as n increases, then reaches its maximal value at n = √ 3(ω 2 −ω 2 1 ) ωω 1 and then decreases as n increases. The corresponding maximal value of ( k) 2 is given by Logarithmic Decrement.
It is curious that the logarithmic decrement at high frequencies behaves as 4πΩ/n = 4πω/ν, i.e., the larger dissipation is, the better wave passes. This peculiar property is related to the fact that the band gap in the elastic medium was related not to the dissipation, but to the localization, which appears to be broken by presence of viscosity. At low frequencies, on the contrary, it is equal to πnΩ = πνω, i.e., it is proportional to the dissipation. At n = 2, the same parameter at which K reaches its maximal value at the characteristic frequency ω = ω 1 , asymptotics for logarithmic decrease at small and large frequencies coincide.
If we differentiate (57) with respect to n, we find that it is a monotone function of n at Ω < 1 (increasing with n) and at Ω > 1 (decreasing with n), and it does not depend on n at Ω = 1 (equals 2π). Numerical examples show a non-monotonic behavior of the logarithmic decrement as a function of frequency for relatively small n near the characteristic frequency ω 1 (see Figure 8).

Discussion
We see that the constrained reduced Cosserat medium in its elastic variant has properties which are similar to the non-constrained reduced Cosserat continuum. It is also a single negative acoustic metamaterial, but the bandgaps for the shear waves in both considered elastic constrained models and for the slower mixed wave in the considered anisotropic case are limited only from below (see Figure 1, Figure 2). The isotropic variant can be obtained as a limit case from the reduced non-constrained isotropic Cosserat model, if we consider in the latter an infinite elastic resistance to the strain θ − ∇ × u/2. The anisotropic model is not the limit case for an anisotropic reduced non-constrained Cosserat theory with a coupling term between volumetric and rotational deformations, considered in [44]. We cannot introduce a similar coupling because of the kinematical hypothesis which we accept in the constrained model. Indeed, in the non-constrained case, it was a term coupling ∇u S and θ − ∇ × u/2. Therefore, it is not strange that one of the branches in the constrained medium is very little dispersive and has no bandgap, unlike the non-constrained continuum. However, both slower waves have some similar to the non-constrained model features (boundary frequencies), due to the fact that both media are reduced, whence some waves cannot propagate.
The constrained reduced viscoelastic Cosserat model demonstrates drastical change in the wave behavior comparing to the elastic case (Figures 3-8). We have seen similar effects in the non-constrained theory [13], but in the model considered here, due to its simplicity, it has been possible to analytically investigate the influence of the dissipation in reduced micropolar continua in more detail. We see that in the viscoelastic constrained isotropic reduced Cosserat medium the bandgap Ω > 1 does not exist anymore, but it converts, partially or completely, into a decreasing part of the dispersion curve. Thus the medium, which was in a certain domain of frequencies a single negative acoustic metamaterial in the elastic case, becomes a double negative acoustic metamaterial in the viscoelastic case, though, generally speaking, in another domain of frequencies.
At high frequencies, the real part of the wavenumber tends to zero, and the group velocity to infinity. The imaginary part tends to −ω 1 /c s at high frequencies (see Figures 3 and 8). This corresponds to a very fast backward wave whose amplitude decreases in space as e −(ω 1 /c s )k·r . These waves become infinitely fast as their wavelength tends to infinity. This peculiar behavior, in our opinion, is related to the existence of the kinematical hypothesis, which introduces a rigid constraint into the medium.
The real part of the wavenumber is limited by a certain value for all frequencies. In all numerical examples k(ω) has only one maximum (Figures 3-8). It does not coincide with the extremum for the imaginary part, when the latter one exists (Figures 3 and 8). Let us discuss the physical meaning of this. All the plane waves existing in the medium have a wavelength that is larger than a certain minimal value, reached at a specific frequency. Above this frequency we have a decreasing part of dispersion curve (a zone of anomalous dispersion), i.e., the medium is a double negative acoustic metamaterial. Thus, introducing the viscosity, we change the type of acoustic metamaterial. For each wavelength larger than this boundary-value, two plane waves co-exist in the continuum. One of them, a forward wave, whose phase velocity has the same sign as the group velocity, has a lower frequency. Another one, with a higher frequency, is a backward wave: its direction of propagation of wave energy, determined by the group velocity, is opposite to the phase velocity. The longer waves are, the faster is the backward wave, and the forward wave is closer to the classical shear wave.
The frequency which separates two zones, of anomalous and normal dispersion, depends on the dissipation parameter, but, as we can suppose from the numerical experiment, cannot be made arbitrarily small in comparison to the characteristic frequency ω 1 for the undamped medium. The group velocity at this point is infinite. Some authors believe that such a point on a dispersion graph corresponds to instabilities in the real medium [45]. Numerical examples (Figures 4-7) show that this frequency depends on the dissipation parameter n in a non-monotonous way.
Logarithmic decrement has, generally speaking, a complex frequency-dependent behavior ( Figure 8), but we know that it decreases with dissipation n for frequencies above the characteristic one (i.e., in the former band gap) and increases with n below it. This, probably, can be explained as follows: in the elastic medium the existence of the bandgap was related not to the dissipation, but to the localization, and the viscosity breaks condition for its existence. At low and high frequencies the logarithmic decrement is proportional to the frequency, but in the middle, it may have (or not have) extremuma, depending on the parameters.
One can observe similar effects also in media of other types. For instance, non-polar dielectrics, if to neglect dissipation, have a bandgap. Dispersion in the low-frequency transparency window is the same as in the elastic constrained reduced isotropic Cosserat medium. The high-frequency transparency window is the same as for the non-constrained elastic reduced isotropic Cosserat medium, but it is absent in our case. When dissipation is taken into account, the bandgap is not present anymore, but a zone of anomalous dispersion appears [45].
In classical viscoelastic metamaterials, considered in [46], viscosity produces similar effects: "transmission is generally worsened in passing bands, while it is enhanced inside locally resonant band gaps". The metamaterial considered in the work mentioned above is a periodic structure of circular metal bars coated with rubber, embedded into a viscoelastic matrix. It is difficult to make a precise comparison with our case, since we do not consider in this work a discrete structure of metamaterial, as it is done in [46]. A similar model is investigated in [47]. It is also concluded that viscosity contracts existing band gaps, and elasticity enhances them. An important difference with the constrained reduced Cosserat medium is that in our case the dissipation destroys completely the bandgap existing in the elastic case, while for the case considered in these two works band gaps are still present in the viscoelastic case. Perhaps the reason is that the considered discrete structure is much more complex than the reduced Cosserat model.
We could expect that these kinds of models may be useful for the description of blocky media and for design of artificial acoustic metamaterials with desired properties.

Materials and Methods
This is a theoretical work with methods based upon fundamental laws of continuum mechanics of complex media in their local form, direct tensor calculus and Fourier transform. Numerical graphs are plotted with help of scripts in python which can be found in Appendix B. Numerical results validate analytical formulae for the obtained dispersion curves.
We describe all the material types considered above and their properties in Table 1.

Medium Wave Properties Classification
Isotropic Classical compression wave. Highly dispersive shear-single negative acoustic elastic rotational wave, which does not propagate above metamaterial at ω > ω 1 boundary frequency ω 1 (large wavenumber limit) Elastic with a Shear wave (same as in isotropic case), non-dispersive single negative acoustic simplest compression wave (for k ⊥ n with classical velocity or metamaterial with respect axisymmetric k||n with larger velocity), and two mixed waves to the shear wave at ω > ω 1 coupling coexist. One mixed wave is weakly, the other is and with respect to the N∇ · u(n · ∇u · n) strongly dispersive with a boundary frequency slower mixed wave between shear and ω 1a < ω 1 depending on the propagation direction at ω > ω 1a volumetric strains Viscoelastic Classical compression wave. Highly dispersive shear-Double negative acoustic isotropic rotational wave. Its dispersion branch has a decreasing metamaterial with respect part and large velocity at large ω; there is a minimal to the shear wave wavelength, k(ω) < k max (n) has at least one above a certain frequency maximum; at small ω the wave tends to classical; band gap disappeared, in its domain dissipation enhances wave propagation, at ω < ω 1 attenuates the wave, at ω 1 does not influence the attenuation Logarithmic decrement is proportional to ω at low and high frequencies (proved analytically and checked numerically).
k(ω) has only one maximum, the group velocity there is infinite, this point separates zones of normal and anomalous dispersion. k has one maximum (at small n) or does not have, it depends on n in a complex way (checked numerically).

Conclusions
We considered three types of linear constrained reduced Cosserat media: elastic isotropic, elastic with specific transversely isotropic coupling between volumetric and shear deformations, and viscoelastic isotropic.
Elastic models present an infinite bandgap above a certain characteristic frequency for some type of waves: shear plane wave and, in case of anisotropic theory, also one of the mixed waves. Therefore they should be classified as single negative acoustic metamaterials. They have also a branch without a bandgap, non-dispersive (in the isotropic case, compressional wave) and weakly dispersive (model with anisotropic coupling, the upper branch of the mixed wave). In the anisotropic case dispersion curves essentially depend on the propagation direction. The isotropic viscoelastic model possesses a classical compressional wave, but the shear wave changes drastically. The bandgap disappears. At high frequencies, we observe a decreasing dispersion curve in some frequency range (i.e., the medium is a double negative acoustic metamaterial in this domain), with a diverging group velocity as the frequency increases. Depending on dissipation, this frequency domain may cover a part of the former band gap ω > ω 1 or extend to the whole former bandgap and occupy also some domain below ω 1 . Above the characteristic frequency ω 1 the logarithmic decrement decreases as the dissipation increases (the viscosity favors wave propagation), and below it, the logarithmic decrement increases with the viscosity increases (which attenuates the wave in this case). Viscous dissipation appears to destroy the wave localization, related to the existence of the bandgap in the elastic case. Shear waves with a wavelength shorter than a certain threshold cannot propagate in this medium. One may thus control the properties of the acoustic metamaterial by means of viscosity.