Next Article in Journal
An Artificial Neural Network Based Solution Scheme for Periodic Computational Homogenization of Electrostatic Problems
Previous Article in Journal
Nonlocal FEM Formulation for Vibration Analysis of Nanowires on Elastic Matrix with Different Materials
Previous Article in Special Issue
Some Applications of Eigenvalue Problems for Tensor and Tensor–Block Matrices for Mathematical Modeling of Micropolar Thin Bodies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Determination of Wave Velocities through the Eigenvalues of Material Objects

by
Mikhail U. Nikabadze
1,2,
Sergey A. Lurie
3,4,
Hovik A. Matevossian
3,5 and
Armine R. Ulukhanyan
2,*
1
Faculty of Mechanics and Mathematics, Lomonosov Moscow State University, 119991 Moscow, Russia
2
Department of Computational Mathematics and Mathematical Physics, Bauman Moscow State Technical University, 105005 Moscow, Russia
3
Federal Research Center "Computer Science and Control", Russian Academy of Sciences, Vavilov str., 40, 119333 Moscow, Russia
4
Institute of Applied Mechanics, Russian Academy of Sciences, Leningradsky prt 7, 125040 Moscow, Russia
5
Moscow Aviation Institute (National Research University), Volokolomskoe Shosse, 4, 125993 Moscow, Russia
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2019, 24(2), 39; https://0-doi-org.brum.beds.ac.uk/10.3390/mca24020039
Submission received: 26 January 2019 / Revised: 9 April 2019 / Accepted: 9 April 2019 / Published: 11 April 2019
(This article belongs to the Special Issue Mathematical Modeling in Physical Sciences)

Abstract

:
The statement of the eigenvalue problem for a tensor–block matrix of any order and of any even rank is formulated. It is known that the eigenvalues of the tensor and the tensor–block matrix are invariant quantities. Therefore, in this work, our goal is to find the expression for the velocities of wave propagation of some medias through the eigenvalues of the material objects. In particular, we consider the classical and micropolar materials with the different anisotropy symbols and for them we determine the expressions for the velocities of wave propagation through the eigenvalues of the material objects.

1. Introduction

The theory of wave propagation in structural elements, for example, in rods, plates, and shells, generally speaking, is much more complicated than the theory of waves in an unbounded or semi-bounded medium. Although the individual waves are the same, but being repeatedly reflected from the boundaries, they create a wave, the exact description of which is difficult to realize. The problem of the theory is to build such a model in which an acceptable compromise between the requirements of accuracy and simplicity would be achieved. The fundamentals of the theory of waves in an elastic cylindrical rod were created by Pohgammer and Cree at the end of the 19th century. They found the various forms of eigenwaves in rods. Further, the studies of the unsteady wave propagation in the elements of elastic structures were carried out, as a rule, using approximate equations, which were obtained from the corresponding statics equations.
The addition of inertial terms to these equations made it possible to construct solutions of wave propagation problems, however, some conclusions were in conflict with the results of the theory of elasticity. Thus, the velocity of propagation of perturbations in the dynamic bending of the rod, determined by the Bernoulli–Euler equation, has no upper limit, while, according to the theory of elasticity, it should be limited by the velocity of longitudinal waves in a continuous medium. The mentioned equation does not allow establishing the presence of wave fronts at all. The longitudinal wave velocity, determined by the approximate equation of the longitudinal oscillations of the rod, although limited, but does not correspond with the velocity from the theory of elasticity [1]. There was an obvious need for a stronger justification of approximate equations, the establishment of areas of their applicability, the construction of refined equations, that is, the development of the dynamic theory of plates and shells. The first results were obtained when corrections were introduced into the equations, which made it possible to more fully take into account the main factors that determine the elastic wave propagation (Rayleigh [2], Tymoshenko [3]). A significant contribution was made by Timoshenko, who proposed a refined equation of dynamic bending (and shear) of the rod. As it was later established by Uflyand [4] and others, the Timoshenko equation, in contrast to the Bernoulli–Euler equation, determines the finite wave propagation velocities and gives results close to exact results that follow from the theory of elasticity. Tymoshenko’s equations and their solutions were studied in a number of papers (see [1,5,6] and others).
Note that, in Nematic Liquid Crystals (NLC), it is possible to impose electric fields and consequently the nematic orientation. A simple linearization procedure proves that in this way we can electrically vary the acoustic properties of the solid fluid mixture studied in [7]. Therefore, we expect that the methods develop in the present paper may allow a semi-analytical study of the tangent evolution operator for solid–fluid ways in NLC.
In [8], the propagation at discontinuity surfaces in second gradient 3D continua. When endowing the discontinuity surface with material properties, one can expect linear propagation condition of the kind studied in the present paper. We conjecture that in this way it would be to possible design frequency filters at solid interfaces.
In [9], it is proven that second gradient continua must be used as homogenized model for hexagonal lattices if high frequency waves have to be considered. The mathematical method presented here can be applied to get closed form expressions generalizing the results in Figure 6 of the cited paper.
In the theory of non-stationary waves, the wave velocity usually means the front velocity of the wave. Consequently, in the dynamic theory of continuous media, the current problem is the determination of the number of waves and the wave propagation velocity in media. The main types of waves in the continuous elastic medium are a longitudinal or expansion wave and a transverse wave. The existence of these waves was determined by Poisson and Stokes in 1828. Other types of waves that can propagate along the surface of an elastic half-space were discovered by Rayleigh in 1885. The theory of surface waves, very important for seismology, was further developed in the works of Love, Stoneley, and Sobolev. Smirnov and Sobolev introduced a class of functionally invariant solutions of the wave equation and a representation of the wave field using functions of a complex variable. This has achieved significant progress in the general theory of waves in an elastic half-space [10]. One of the main problems of the theory of waves in the half-space is the problem of Lemb [11], namely, the action of a local source on the border of a half-space. One of the main conclusions obtained from the analysis of non-stationary waves in structural elements is that, in the process of propagation of a non-stationary wave through an elastic layer, as a result, in the reflections from surfaces, we obtain an interference pattern, in general, correctly described by simple approximation equations. However, to discover some details of the asymptotics, as well as to make the equations suitable for studying the process, clarifications are necessary, which increase the order of the equations [6]. Note that various methods for obtaining approximate equations for the theories of thin bodies are described in [12]. It also shows various versions of the refined equations for classical and micropolar theories of thin bodies, which can be used to study the wave propagation processes in thin structures. Here, we obtain expressions for velocities of the wave propagation in some unbounded media through the eigenvalues of material objects. Note that some issues of applications of eigenvalue problems for tensor and tensor–block matrices for mathematical modeling of micropolar thin bodies are also considered in [13], while some issues of the modeling of multilayer thin bodies are given in [14].

2. Statement of Eigenvalue Problem of a Tensor–Block Matrix of Any Even Rank

Find all tensor columns U which satisfy equation
M p U = λ U ,
where λ is scalar, and
U = U 1 U 2 · · · U m , M = A 11 A 12 A 13 A 1 m A 21 A 22 A 23 A 2 m A 31 A 32 A 33 A 3 m A m 1 A m 2 A m 3 A m m ,
M = M i 1 i 2 i p j 1 j 2 j p R i 1 i p R j 1 j p = M i j R i R j , R i 1 i p = r 1 r p , R i p R j = g i j , A k l = A k l , j 1 j 2 j p i 1 i 2 i p R i 1 i 2 i p R j 1 j 2 j p = A k l , · j i · R i R j , U k = U k , i 1 i 2 i p R i 1 i p , k , l = 1 , m ¯ , i 1 , i 2 , , i p , j 1 , j 2 , , j p = 1 , n ¯ , i , j = 1 , N ¯ , N = n p .
Note that this problem is solved for the tensor of any even rank and the tensor–block matrix (TBM) of any even rank consisting of four tensors, as well as for the tensor and the tensor–block matrix of the fourth rank, and published in [15]. Therefore, here we do not dwell on the presentation of this problem with the aim of shortening the letter, but, if necessary, we refer to the work mentioned in the previous sentence. We also note that, solving the eigenvalue problem for a TBM of any even rank consisting of four tensors, there is no difficulty in solving the analogous problem for the TBM M (see (1)). Thus, we assume that the eigenvalue problem for a TBM of any order and of any even rank is solved and we consider some of its applications below. In addition, we note that the eigenvalue problem for the tensor of any even rank is a special case of the eigenvalue problem for a TBM of the same rank. Thus, solving the eigenvalue problem for a TBM of any even rank, we believe that it is solved for a tensor of the same rank [15].

3. Kinematic and Dynamic Conditions on the Surface of a Strong Discontinuity in Micropolar Mechanics

Consider a moving regular surface in an unbounded space, x = ( x 1 , x 2 , x 3 ) R 3 , t 0 , whose equation in a fixed Cartesian coordinate system is given by
ψ ( x , t ) = 0 .
The regularity of Equation (2) means the existence of the unit normal vector n ( x , t ) at each point of the surface in Equation (2) at time t:
n ( x , t ) = x ψ | x ψ | , x ψ = k i i ψ , i = 1 , 2 , 3 ,
as well as the required number of times continuously differentiable functions ψ of x. For our purpose, is sufficient the existence of the unit normal, i.e., ψ C 1 , is sufficient. Here, k i is an orthonormal basis of the Cartesian coordinate system, and C 1 is the set of continuously differentiable functions. We use the usual rules of tensor calculus [15,16,17,18,19,20]. We mainly preserve the notation and conventions of the previous works, capital Latin indices assume the values 1, 2. Over repeated indices there is a summation.
Note that the kinematics of a micropolar medium is described by independent vectors of displacements u and rotations φ .
Definition 1.
If the vectors u and φ are continuous when passing through the surface in Equation (2), and the first derivatives k u , t u , k φ , and t φ undergo a discontinuity such that on each side of this surface they take different finite values, then the surface in Equation (2) is called the surface of a strong discontinuity (or the wave of stress and moment stress).
Definition 2.
If the vectors u , k u , t u , φ, k φ and t φ are continuous when passing through the surface in Equation (2), and the second derivatives of the vectors u and φ with respect to x k and t undergo a discontinuity such that on each side of this surface they take different finite values, then the surface in Equation (2) is called the surface of a weak discontinuity (or the acceleration wave).
We introduce the notation: t = / t , i = / x i . Next, the surface of a strong discontinuity is studied.
As in [1,21,22,23], the velocity of moving an arbitrary point of the surface in Equation (2) in the direction of the normal of this surface at this point, we denote by c:
c = 1 | x ψ | t ψ , ψ = k i i ψ .
If the wavefront (the surface of strong discontinuity) of Equation (2) moves in a medium having a velocity field v ( x , t ) , then the velocity of the wavefront relative to the particles of the medium is determined by the formula
θ = c v n , v n = n · v .

3.1. Kinematic Conditions on the Surface of a Strong Discontinuity

The kinematic conditions on the surface of strong discontinuity can be obtained in the same way, as is done, for example, in [1,21,22,23] for the classical case. In this case, the difference lies in that instead of a motion vector should be considered two independent vectors u ( x , t ) , φ ( x , t ) , and they have the form
c [ i u ] + n i [ t u ] = 0 , c [ i φ ] + n i [ t φ ] = 0 ,
where c is determined by the first formula of Equation (3), n i are the components of the unit normal n to the wave front. Writing [ w ] = w + w , where w = u or w = φ , means a jump in the value of w relative to the wave front, w ( w + ) is the limiting value of w when the arbitrarily chosen point is in front (behind) the wave front to the point at the wave front. The square brackets [ ] are called the jump operator.

3.2. Laws of Conservation of Mass and the Tensor of Moments of Inertia at the Wave Front

Applying the law of conservation of mass and the law of the tensor of moments of inertia
d d t V ρ d V = 0 , d d t V J ˜ d V = 0 ,
where ρ is a material density, and J ˜ is the density of the inertia tensor (a special dynamic characteristic of the medium) of the particles of the medium [24,25,26], to the elementary cylinder isolated in the medium. After simple transformations analogous to the classical case [1,21,22,23], we obtain the required laws
[ θ ρ ] = 0 , [ θ J ˜ ] = 0 ,
where θ is defined in Equation (4).

3.3. Dynamic Conditions on the Wave Front

Dynamic conditions on the wave front can be easily obtained with the help of the law on the change of momentum and of the theorem on the change in the angular momentum of the internal rotational motions of the particles of the medium, which are defined as follows:
d d t V ρ v d V = V ρ F d V + Σ P ( n ) d Σ , d d t V J ˜ · ω d V = V ( ρ m + C 2 P ˜ ) d V + Σ μ ( n ) d Σ ,
where F is is the mass force, m is the mass moment, P ( n ) and μ ( n ) are the voltage and moment stress vectors on the area with a unit vector of the normal n respectively, V is the volume of the body, Σ is the boundary of the body, and v = u ˙ , ω = φ ˙ , where the dot over the letter indicates the time derivative. As in the classical case [1,21,22,23], we apply Equation (7) to the above-mentioned cylinder and, taking into account Equation (6), we obtain the conditions
ρ θ [ v ] = [ P ( n ) ] = n · [ P ˜ ] , θ J ˜ · [ ω ] = [ μ ( n ) ] = n · [ μ ˜ ] ,
where ρ θ = ρ θ = ρ + θ + , θ J ˜ = θ J ˜ = θ + J ˜ + .
We note that the relations in Equation (8), in which the Cauchy formulas are taken into account, P ( n ) = n · P ˜ and μ ( n ) = n · μ ˜ , are also valid for any medium. We note that the conditions in Equation (8) and their analogs for different media are derived in [24].

4. Determination of Wave Propagation Velocities in an Infinite Micropolar Solid

Having kinetic (Equation (5)) and dynamic (Equation (8)) conditions on the front of the wave, it is easy to find an equation for determining the propagation velocities of waves in any infinite micropolar medium, including in an infinite micropolar solid. Consider a micropolar solid body, the constitutive relations of which are represented in the form (see [12,27,28])
P ˜ = A 2 u + B 2 φ A 2 C · φ b ˜ ϑ , μ ˜ = C 2 u + D 2 φ C 2 C · φ β ˜ ϑ ,
where b ˜ = A 2 a ˜ + B 2 d ˜ , β ˜ = C 2 a ˜ + D 2 d ˜ . Here, 2 is an inner 2-product, C = B T , A , D are the material tensors (tensors of the elastic moduli) of the fourth rank, ϑ = T T 0 is the temperature drop, b ˜ and β ˜ are the tensors of thermomechanical properties, and T in the upper corner of the tensor means the sign of transposition.
It is known [24,26] that B is an asymmetric tensor, as mentioned above. In particular, in [29], it is proved that B is a symmetric tensor, and C = B . Here, we derive these relations when C = B T , since it is easy to obtain from them the corresponding relations for the case C = B .
Next, for simplicity, let us consider isothermal processes, i.e., we assume that ϑ = 0 . Then, assuming that the material tensors do not undergo a discontinuity when passing through the front of the wave and applying the jump operator to Equation (9), we have
[ P ˜ ] = A 2 [ u ] + B 2 [ φ ] , [ μ ˜ ] = C 2 [ u ] + D 2 [ φ ] .
Multiplying Equation (5) by k i with subsequent summation over i (hereinafter, we omit the index x of the operator x ), we obtain
[ u ] = 1 c n [ v ] , [ φ ] = 1 c n [ ω ] .
From Equation (10) taking into account Equation (11), we find
[ P ˜ ] = 1 c ( A 2 n [ v ] + B 2 n [ ω ] ) , [ μ ˜ ] = 1 c ( C 2 n [ v ] + D 2 n [ ω ] ) ,
with the help of which, from the dynamic conditions in Equation (8), we arrive at the relations
( n · A 2 n E ˜ ) · [ v ] + ( n · B 2 n E ˜ ) · [ ω ] = c ρ θ E ˜ · [ v ] , ( n · C 2 n E ˜ ) · [ v ] + ( n · D 2 n E ˜ ) · [ ω ] = c θ J ˜ · [ ω ] ,
where E ˜ is the unit tensor of the second rank.
We note that, in view of the laws of conservation of mass and the tensor of moments of inertia, Equation (6), the expressions ρ θ and θ J ˜ in Equation (12) can be replaced by ρ + θ + and θ + J ˜ + , respectively, since in front of the wave front ρ + and J ˜ + , and also v n + can be considered known, but to simplify the recording we do not do this.
We introduce the notation
A ˜ = ( 1 / ρ ) n · A 2 n E ˜ , B ˜ = ( 1 / ρ ) n · B 2 n E ˜ , C ˜ = J ˜ 1 · ( n · C 2 n E ˜ ) , D ˜ = J ˜ 1 · ( n · D 2 n E ˜ ) , M ˜ = A ˜ B ˜ C ˜ D ˜ , [ V ] = [ v ] [ ω ] ,
where M ˜ is TBM, V is a vector column of vectors of linear and angular velocities, and [ V ] is the jump of this vector column. Equation (12) can be represented in the form
A ˜ B ˜ C ˜ D ˜ · [ v ] [ ω ] = λ [ v ] [ ω ] , λ = c θ
or, in short,
M ˜ · [ V ] = λ [ V ] .
The equality in Equation (15) (see also Equation (14)) represents a homogeneous system of six algebraic equations with respect to six unknowns (two vectors [ v ] and [ ω ] having a nontrivial solution). For this system to have a nontrivial solution, it is necessary and sufficient that its determinant be zero. Since the determinant is of the sixth order, and θ + = c v n + , then, from the equality to zero of this determinant, we obtain an algebraic equation of degree 6 with respect to c 2 , which is the desired dispersion equation for determining the velocities of waves in an infinite anisotropic micropolar body and their number in a given direction.
Thus, we grt the eigenvalue problem in the form of Equation (15) for TBM. Obviously, λ = c θ is an eigenvalue, and [ V ] is its corresponding jump of the vector column. By virtue of Equation (15), the dispersion equation (the characteristic equation for M ˜ ) can be written in the form
det ( M ˜ λ E ˜ ) = 0 ,
where E is a unit TBM of the second rank. In the expanded form, Equation (16) can be written as follows:
λ 6 I 1 ( M ˜ ) λ 5 + I 2 ( M ˜ ) λ 4 I 3 ( M ˜ ) λ 3 + I 4 ( M ˜ ) λ 2 I 5 ( M ˜ ) λ + I 6 ( M ˜ ) = 0 , I 6 ( M ˜ ) = det ( M ˜ ) .
It is seen that the dispersion in Equation (17) is an algebraic equation of the sixth degree and must have six roots (eigenvalues), counting each root as many times as its multiplicity. Each multiple root determines the square of the velocity of one wave. Hence, in an arbitrary anisotropic infinite micropolar medium, in the general case, no more than six waves can arise in each direction. Note that, based on the dispersion equation of the form in Equation (17), it is easy to establish the number of waves arising in a micropolar elastic medium for different anisotropy. Thus, it is sufficient to find the invariants of the TBM that appear in Equation (17), and then solve the equation itself. Invariants are easily found through the first invariants of powers M ˜ . We have ([15,17,30,31,32])
S k = I k ( M ˜ ) = 1 k ! s 1 1 0 0 s 2 s 1 0 0 s k 1 s k 2 s 1 k s k s k 1 s 2 s 1 , k = 1 , 6 ¯ ,
S k = I k ( M ˜ ) , s k = I 1 ( M ˜ k ) , k = 1 , 6 ¯ , M ˜ k = M ˜ · M ˜ · · M ˜ k ,
where I k ( M ˜ ) , k = 1 , 6 ¯ , denote the invariants of TBM M ˜ . In this case, the inverse relations to Equation (18) are represented in the form
s k = I 1 ( M ˜ k ) = S 1 1 0 0 2 S 2 S 1 1 0 k S k S k 1 S k 2 S 1 , k = 1 , 6 ¯ .
If the material has a center of symmetry, then B = C T = 0 . Then, B ˜ = 0 and C ˜ = 0 , and the TBM M ˜ becomes a diagonal TBM. For such a matrix, the characteristic equation and eigenvalues (wave velocities) are easily found. We have
det ( M ˜ λ E ˜ ) = det A ˜ λ E ˜ 0 ˜ 0 ˜ D ˜ λ E ˜ = det ( A ˜ λ E ˜ ) det ( D ˜ λ E ˜ ) = 0 ,
where 0 ˜ is the zero tensor of rank 2. Equation (19) is equivalent to the following equations:
λ 3 I 1 ( A ˜ ) λ 2 + I 2 ( A ˜ ) λ I 3 ( A ˜ ) = 0 , λ 3 I 1 ( D ˜ ) λ 2 + I 2 ( D ˜ ) λ I 3 ( D ˜ ) = 0 .
Thus, if the material has a center of symmetry, then in this case to determine the wave velocities we have two cubic equations in Equation (20), which are easily solved. Based on Equation (20), we can conclude that, in an arbitrary anisotropic micropolar elastic medium with a center of symmetry, no more than six waves appear in each direction.
We note that in the case of a classical medium, we have one cubic equation analogous to the first of the equations in Equation (20) provided that A ˜ is determined using the first relation in Equation (13), where A is the elastic modulus tensor. Consequently, in an arbitrary anisotropic classical elastic medium, no more than three waves appear in each direction. Note also that, in the case of a micropolar medium, it is reasonable to call M ˜ the dispersion TBM, in the case that the classical medium A ˜ is the dispersion tensor, since their characteristic equations are the dispersion equations.
Having the dispersion equations derived above for the micropolar (classical) medium, we can find the propagation velocities of the waves in the media under consideration for any anisotropy. Of course, the required velocities can be found both in the traditional representation of material objects (TBM, tensors of elastic moduli) and with the help of their eigenvalues, solving in advance the eigenvalue problems of the corresponding object in the case under consideration. The eigenvalue problems for the tensor and TBM of any even rank are solved in [15,17,30,31], in which the notion of the anisotropy (structure) symbol of the material is introduced and the classifications of micropolar and classical media are given.
Therefore, for each material included in these classifications, on the basis of the corresponding dispersion equation obtained for the micropolar (classical) medium, one can obtain a dispersion equation for the considered micropolar (classical) medium, and then determine the number of waves and their propagation velocity.
Definition 3.
The symbol { α 1 , α 2 , , α k } , where k is the number of different eigenvalues of the TBM (tensor), and α i is the multiplicity of the eigenvalue λ i , i = 1 , 2 , , k , is called the symbol of the anisotropy (structure) of the TBM (tensor).
Note that the anisotropy symbol is defined for the TBM (tensor) of even rank. In this case, the symbol of the anisotropy of the material TBM (material tensor) is also called the symbol of anisotropy of the material.
Next, we consider some particular cases of materials and find the propagation velocities of the waves in them using the eigenvalues of the corresponding tensor objects.

5. Application

5.1. Classical Materials with Anisotropy Symbols { 1 , 5 } and { 5 , 1 }

For materials with anisotropy symbols { 1 , 5 } and { 5 , 1 } , the elastic modulus tensor is represented in the form
A = ( λ 1 λ 2 ) a ˜ 1 a ˜ 1 + λ 2 E , B = μ 1 E ( μ 1 μ 6 ) b ˜ 6 b ˜ 6 .
Here, E = k = 1 6 a ˜ k a ˜ k = k = 1 6 b ˜ k b ˜ k = ( 1 / 2 ) ( C ( 2 ) + C ( 3 ) ) is the unit tensor of rank 4, a ˜ k and b ˜ k , k = 1 , 6 ¯ are complete orthonormal systems of proper tensors for tensors A and B , respectively, C ( 2 ) and C ( 3 ) are isotropic tensors of rank 4, λ 1 , λ 2 and μ 1 , μ 6 are eigenvalues, and a ˜ 1 = a ˜ 1 T and b ˜ 6 = b ˜ 6 T are proper tensors corresponding to eigenvalues λ 1 and μ 6 , respectively.
Since a ˜ 1 ( b ˜ 6 ) is a symmetric tensor satisfying the orthonormality condition a ˜ 1 2 a ˜ 1 = 1 ( b ˜ 6 2 b ˜ 6 = 1 ), in a basis constructed using a basis of an arbitrary coordinate system, it is characterized by five components, and in the main basis for a ˜ 1 ( b ˜ 6 )—two components. Hence, the tensor A ( B ) in a basis formed using a basis to an arbitrary coordinate system is characterized by seven parameters, namely two eigenvalues and five components of the tensor a ˜ 1 ( b ˜ 6 ), and in the basis formed by means of the canonical basis for a ˜ 1 ( b ˜ 6 ), four parameters, namely two eigenvalues and two components of the tensor a ˜ 1 ( b ˜ 6 ).
Proposition 1.
Let a ˜ = a ˜ T . Then, a ˜ a ˜ is an isotropic tensor of rank 4 if and only if a ˜ is a spherical tensor. In this case, if a ˜ 2 a ˜ = 1 , then a ˜ = [ ( ± 3 ) / 3 ] E ˜ , and a ˜ a ˜ = ( 1 / 3 ) C ( 1 ) .
Here, C ( 1 ) = E ˜ E ˜ is the first of three isotropic tensors of rank 4. From this statement, it follows that the tensors A and B of Equation (21) is traditionally isotropic if and only if a ˜ 1 a ˜ 1 = b ˜ 6 b ˜ 6 = ( 1 / 3 ) C ( 1 ) . In this case, they can be written in the form
A = ( 1 / 3 ) ( λ 1 λ 2 ) C ( 1 ) + λ 2 E = λ C ( 1 ) + 2 μ E , B = μ 1 E ( 1 / 3 ) ( μ 1 μ 6 ) C ( 1 ) .
We note that an isotropic material whose properties are characterized by the tensor A ( B ) has a positive (negative) Poisson’s ratio [15,17,31]. In addition, the tensor A of Equation (22) is presented in the traditional form (by the Lamé coefficients λ = 1 / 3 ( λ 1 λ 2 ) , μ = 1 / 2 λ 2 ).
It is easy to see that as tensors A and B of Equation (21), and the corresponding dispersion tensors A ˜ and B ˜ have the same structure and are represented as
A ˜ = 1 ρ n · A · n = a ( E ˜ + n n ) + b u ˜ , B ˜ = 1 ρ n · B · n = f ( E ˜ + n n ) + g v ˜ ,
a = λ 2 / ( 2 ρ ) = μ / ρ > 0 , b = ( λ 1 λ 2 ) / ρ = ( 3 λ ) / ρ > 0 , u ˜ = n · a ˜ 1 n · a ˜ 1 , f = μ 1 / ( 2 ρ ) > 0 , g = ( μ 6 μ 1 ) / ρ < 0 , v ˜ = n · b ˜ 1 n · b ˜ 1 .
Assuming in Equation (23) a ˜ 1 = b ˜ 6 = ( ± 3 / 3 ) E ˜ , we obtain the dispersion tensors A ˜ and B ˜ , corresponding to isotropic tensors in Equation (22),
A ˜ = a E ˜ + b 1 n n , B ˜ = f E ˜ + g 1 n n ,
a = λ 2 2 ρ = μ ρ > 0 , b 1 = 2 λ 1 + λ 2 6 ρ = λ + μ ρ > 0 , u ˜ = n · a ˜ 1 n · a ˜ 1 , f = μ 1 / ( 2 ρ ) > 0 , g 1 = ( μ 1 + 2 μ 6 ) / ( 6 ρ ) > 0 , v ˜ = n · b ˜ 1 n · b ˜ 1 .
Obviously, the tensors in Equation (25), which are special cases of the tensors in Equation (23), similar to the latter have the same structure. In this regard, below we consider the first tensor in Equation (23), and for the rest of tensors in Equations (23) and (25) we obtain the corresponding relations by an appropriate renaming of the coefficients and tensors.
Thus, we find the propagation velocities of the waves in the material { 1 , 5 } whose dispersion tensor has the form of Equation (23), and the characteristic equation is represented as the first equation of Equation (20). To do this, we first find I 1 ( A ˜ k ) , k = 1 , 2 , 3 , and then using Equation (18), which is true for a tensor of the corresponding rank, we find I k ( A ˜ ) , k = 1 , 2 , 3 . After simple calculations, we have
I 1 ( A ˜ ) = 4 a + b I 1 ( u ˜ ) , I 1 ( A ˜ 2 ) = 6 a 2 + 2 a b [ I 1 ( u ˜ ) + I 2 ] + b 2 I 1 2 ( u ˜ ) , I 1 ( A ˜ 3 ) = 10 a 3 + 3 a 2 b [ I 1 ( u ˜ ) + 3 I 2 ] + 3 a b 2 I 1 ( u ˜ ) [ I 1 ( u ˜ ) + I 2 ] + b 3 I 1 3 ( u ˜ ) , I 2 ( A ˜ ) = 5 a 2 + 3 a b I 1 ( u ˜ ) a b I 2 , I 3 ( A ˜ ) = 2 a 3 + 2 a 2 b I 1 ( u ˜ ) a 2 b I 2 , I 1 ( u ˜ ) = n · a ˜ 1 2 · n , I = a ˜ 1 2 n n .
Further, by virtue of the corresponding invariants in Equation (27), constructing the characteristic equation from Equation (20) for the considered tensor A ˜ and solving it, we obtain the following expressions for the roots:
η 1 = a , η 2 , 3 = 1 2 3 a + b I 1 ( u ˜ ) ± [ a b I 1 ( u ˜ ) ] 2 + 4 a b I 2 .
Theorem 1.
The dispersion tensor is positive definite, and the dispersion TBM is positively determined.
It follows from this theorem that the eigenvalues of the dispersion tensor and the dispersion tensor–block matrix are positive. It is easy to prove that the eigenvalues in Equation (28) obtained above are positive. Knowing the roots of the characteristic equation of the tensor A ˜ (see Equation (28)), it is easy to find the propagation velocities of waves in an initially resting medium. We have
c 1 = η 1 , c 2 = η 2 , c 3 = η 3 .
Thus, in an initially resting medium with the structure symbol { 1 , 5 } using Equation (29), we can find the propagation velocities of the waves in an arbitrary direction. In the general case, their number is not more than three. Note that, if the medium does not rest, then to determine the wave velocities instead of Equation (29) we have formulas
c 1 ( c 1 v n + ) = η 1 , c 2 ( c 2 v n + ) = η 2 , c 3 ( c 3 v n + ) = η 3 .
It is seen that each relation from Equation (30) is a square equation with respect to the wave propagation velocity. Solving them we obtain explicit expressions for the propagation velocities of waves in the initially disturbing medium with the structure symbol { 1 , 5 } . It can be assumed that the number of waves arising in such a medium in an arbitrary direction are no more than six, but not fewer than three. To conduct this study, and similar to the above for A ˜ (the first equation of Equation (23)) and for B ˜ (the second equation of Equation (23)) is not difficult. Obviously, in the latter case, it is sufficient in the above formulas for A ˜ , as mentioned above, to replace a, b, A ˜ and u ˜ by f, g, B ˜ and v ˜ , respectively. In connection with this simplicity, we do not dwell on these questions.
Next, consider the first dispersion tensor in Equation (25), which corresponds to the first tensor in Equation (22), characterizing the properties of a traditionally isotropic material and being a particular case of the first tensor in Equation (23). It is easy to see that in this case, analogous to Equation (27), the relations are represented in the form
I 1 ( A ˜ ) = 3 a + b 1 , I 1 ( A ˜ 2 ) = 2 a 2 + ( a + b 1 ) 2 , I 1 ( A ˜ 3 ) = 2 a 3 + ( a + b 1 ) 3 , I 2 ( A ˜ ) = a ( 3 a + 2 b 1 ) , I 3 ( A ˜ ) = a 2 ( a + b 1 ) .
Taking into account the invariants corresponding to Equation (31), from the first equation of Equation (20), we obtain the characteristic equation for the investigated dispersion tensor A ˜ . Solving it, we get
η 1 = η 2 = a = λ 2 2 ρ = μ ρ , η 3 = a + b 1 = λ 1 + 2 λ 2 3 ρ = λ + 2 μ ρ .
The velocities of the propagation of waves in an initially resting infinite isotropic elastic medium by virtue of Equation (32) are determined by formulas
c 1 = η 1 = λ 2 2 ρ = μ ρ , c 2 = η 3 = λ 1 + 2 λ 2 3 ρ = λ + 2 μ ρ .
It can be seen that by Equation (33) the propagation velocities of the waves are expressed both in terms of the eigenvalues and the Lame’s parameters. In the case under consideration, writing out the analogous to Equation (30) relations and investigating them is not difficult, so we do not dwell on this.
From the first tensor of Equation (25), we see that any vector perpendicular to n and located in the tangent plane to the wave surface, is an eigenvector of the dispersion tensor A ˜ , and that n is its eigenvector. Consequently, the system of vectors ( s , l , n ) , where s and l are mutually perpendicular unit tangent vectors to the wave surface, is the complete orthonormal system of eigenvectors of the considered dispersion tensor A ˜ (the first tensor of Equation (25)). The eigenvectors of the tensor A ˜ can be found by solving the system of equations corresponding to these vectors, but, in the case under consideration, there is no such need, since they can be easily guessed. Consequently, the canonical representation of the dispersion tensor A ˜ in view of what has been said above, has the form
A ˜ = η 1 ( s s + l l ) + η 3 n n = η 1 I ˜ + η 3 n n , I ˜ = s s + l l .
It is easy to see that, by virtue of Equation (33), from Equation (34) for the initially resting medium, we obtain
A ˜ = c 1 2 ( s s + l l ) + c 2 2 n n = c 1 2 I ˜ + c 2 2 n n = v s v s + v l v l + v n v n ,
and we also have
V ˜ = A ˜ = c 1 ( s s + l l ) + c 2 n n = c 1 I ˜ + c 2 n n = v s s + v l l + v n n = s v s + l v l + n v n ,
where v s = c 1 s , v l = c 1 l and v n = c 2 n .
From Equations (34), (35) and (36), it is seen that the tensors A ˜ and V ˜ are represented by the sum of two or three orthogonal tensors. In addition, it is seen from Equation (33) that the eigenvalues of the tensor V ˜ are the wave propagation velocities, and the eigenvectors coincide with the eigenvectors of the tensor A ˜ , which is quite natural. In this connection, the tensor V ˜ can be called the velocity tensor. Note that the vector equation for determining the eigen-tensors for A ˜ (the first equation of Equation (25)) can be written in the form
{ ( a η ) I ˜ + ( b 1 + a η ) n n } · [ v ] = 0 .
Consequently, for an arbitrary motion of the medium, the expression [ v ] can be represented in the form
[ v ] = [ v τ ] + [ v n ] , v τ = v τ τ , v n = v n n , τ n .
Taking into account Equation (38), from Equation (37), we have
{ ( a η ) I ˜ } · [ v τ ] + { ( b 1 + a η ) n n } · [ v n ] = 0 .
Note that the following theorem holds, which can be proved by the kinematic conditions (the first equation of Equation (5)).
Theorem 2.
rot u = × u = 0 if and only if [ v ] n , and div u = · u = 0 if and only if [ v ] n .
Note also that, if [ v ] n or [ v ] n , then the wave is called longitudinal or transverse, respectively.
If now, the motion of the medium is such that [ v ] n ( [ v ] n ), i.e., v = v n = v n , ( v = v τ = v τ ), then, from Equation (39), we get
{ ( b 1 + a η ) n n } · [ v n ] = 0 { ( a η ) I ˜ } · [ v τ ] = 0 ,
and hence it is easy to obtain the propagation velocities of longitudinal (transverse) waves. They are written above in Equation (33) (see also Equations (30) and (32)), therefore, we do not write them out again here.
Note that a study analogous to the one above can be carried out in a more general case, but for brevity we do not dwell on this.
Further, before considering the materials of other structures, we note that for a fourth rank tensor we apply the four-index and two-index representations ([15,17,30,31])
A = A i j k l e i e j e k e l = m , n = 1 9 A m n e ˜ m e ˜ n = A m n e ˜ m e ˜ n , i , j , k , l = 1 , 2 , 3 , m , n = 1 , 9 ¯ , e i · e j = δ i j , i , j = 1 , 2 , 3 ; e ˜ s = e s e s , s = 1 , 2 , 3 , e ˜ 4 = 1 2 ( e 1 e 2 + e 1 e 2 ) , e ˜ 5 = 1 2 ( e 2 e 3 + e 3 e 2 ) , e ˜ 6 = 1 2 ( e 3 e 1 + e 1 e 3 ) , e ˜ 7 = 1 2 ( e 1 e 2 e 1 e 2 ) , e ˜ 8 = 1 2 ( e 2 e 3 e 3 e 2 ) , e ˜ 9 = 1 2 ( e 3 e 1 e 1 e 3 ) , e ˜ m 2 e ˜ m = δ m n , m , n = 1 , 9 ¯ .
At the same time, if the components of the tensor A have symmetries A i j k l = A k l i j = A j i k l , then we have
A = A i j k l e i e j e k e l = m = 1 6 n = 1 6 A m n e ˜ m e ˜ n = A m n e ˜ m e ˜ n , i , j , k , l = 1 , 2 , 3 , m , n = 1 , 6 ¯ .
Further, the two-index notation of the fourth rank tensor is mainly applied (see Equations (41) and (42)).

5.2. Classical Material with the Anisotropy Symbol { 1 , 2 , 3 } (Cubic Symmetry)

In this case, we have the following canonical representation of the elastic modulus tensor A :
A = ( λ 1 λ 4 ) W ˜ 1 W ˜ 1 + ( λ 2 λ 4 ) ( W ˜ 2 W ˜ 2 + W ˜ 3 W ˜ 3 ) + λ 4 E ( E = m = 1 6 W ˜ m W ˜ m ) ,
where the eigenvalues and proper tensors are represented in the form
λ 1 = A 11 + 2 A 12 , λ 2 = λ 3 = A 11 A 12 , λ 4 = λ 5 = λ 6 = A 44 , W ˜ 1 = ± 3 3 ( e ˜ 1 + e ˜ 2 + e ˜ 3 ) = ± 3 3 E ˜ , W ˜ 2 = ± 2 2 ( e ˜ 1 e ˜ 2 ) , W ˜ 3 = ± 6 6 ( e ˜ 1 + e ˜ 2 2 e ˜ 3 ) , W ˜ 4 = e ˜ 4 , W ˜ 5 = e ˜ 5 , W ˜ 6 = e ˜ 6 .
Taking into account the expressions for the intrinsic tensors in Equation (44), from Equation (43), we get
A = 1 3 ( λ 1 λ 2 ) C 1 + ( λ 2 λ 4 ) k = 1 3 e ˜ k e ˜ k + λ 4 E , E = 1 2 ( C 2 + C 3 ) .
By virtue of Equation (44), for the dispersion tensor A ˜ , we have the expression
A ˜ = a n n + b E ˜ + d k = 1 3 n k 2 e ˜ k .
Based on Equation (46), it is easy to find I 1 ( A ˜ k ) , k = 1 , 2 , 3 , and then using Equation (18), which are also true for A ˜ , to find I k ( A ˜ ) , k = 2 , 3 . In fact, after simple calculations, we find
I 1 ( A ˜ ) = a + 3 b + d , I 1 ( A ˜ 2 ) = a ( a + 2 b ) + 2 b d + 3 b 2 + ( 2 a d + d 2 ) k = 1 3 n k 4 , I 1 ( A ˜ 3 ) = a ( a 2 + 3 a b + 3 b 2 ) + 3 b 2 ( b + d ) + d [ 3 a ( a + 2 b ) + 3 b d ] k = 1 3 n k 4 + d 2 ( 3 a + d ) k = 1 3 n k 6 , I 2 ( A ˜ ) = 1 2 [ 6 b 2 + 4 a b + 2 a d + 4 b d + d 2 d ( 2 a + d ] k = 1 3 n k 4 ] , I 3 ( A ˜ ) = 1 3 ! [ 6 a b 2 + 3 a d 2 + + 3 b d 2 + + 6 a b d + 6 b 3 + 6 b 2 d + d 3 ( 9 a d 2 + 6 a b d + 3 b d 2 + 3 d 3 ) k = 1 3 n k 4 + + ( 6 a d 2 + 2 d 3 ) k = 1 3 n k 6 ] .
Knowing the invariants I k ( A ˜ ) , k = 1 , 2 , 3 , in view of, for example, the first equation of Equation (20), we make the characteristic equation for the tensor under consideration, which is an algebraic equation of the third degree and always has three positive roots, counting each root as many times as its multiplicity. Thus, in the case under consideration, for any direction, depending on the multiplicity of the roots, one can determine the number of waves and find their propagation velocities. Consider, for example, three mutually perpendicular directions of wave propagation in the medium under consideration, which are determined by the following values of the components of the normal: n i = δ i 1 , n i = δ i 2 n i = δ i 3 , i = 1 , 2 , 3 . Note that for each of these directions the first invariant of the tensor A ˜ does not depend on the direction n , and the second and third invariants have the same values. In particular,
I 1 ( A ˜ ) = a + 3 b + d , I 2 ( A ˜ ) = b 2 ( 2 a + 3 b + 2 d ) , I 3 ( A ˜ ) = b 2 ( a + b + d ) .
By Equation (47), the dispersion equation for each of these directions is the same and has roots
μ 1 = λ 1 + 2 λ 2 3 ρ , μ 2 = μ 3 = λ 4 2 ρ .
Consequently, in each of these directions, in an infinite medium, in the case of cubic symmetry, two waves arise and on the basis of Equation (48) their velocities are given by formulas
c 1 = ( λ 1 + 2 λ 2 ) / ( 3 ρ ) , c 2 = λ 4 / ( 2 ρ ) .

5.3. Classical Material with the Anisotropy Symbol { 1 , 1 , 2 , 2 } (Transversal Isotropy)

In this case, the canonical representation of the elastic modulus tensor A has the form
A = μ 1 w ˜ 1 w ˜ 1 + μ 2 w ˜ 2 w ˜ 2 + μ 3 ( w ˜ 3 w ˜ 3 + w ˜ 4 w ˜ 4 ) + μ 5 ( w ˜ 5 w ˜ 5 + w ˜ 6 w ˜ 6 ) ,
where the eigenvalues are given by formulas
μ 1 = 1 2 ( A 11 + A 12 + A 33 ( A 11 + A 12 A 33 ) 2 + 8 A 13 2 ) , μ 3 = μ 4 = A 11 A 12 , μ 2 = 1 2 ( A 11 + A 12 + A 33 + ( A 11 + A 12 A 33 ) 2 + 8 A 13 2 ) , μ 5 = μ 6 = A 55 ,
and the proper tensors are represented in the form
w ˜ 1 = 2 2 sin α ( e ˜ 1 + e ˜ 2 ) + cos α e ˜ 3 = 2 2 sin α I ˜ + cos α e ˜ 3 , w ˜ 2 = 2 2 cos α ( e ˜ 1 + e ˜ 2 ) + sin α e ˜ 3 = 2 2 cos α I ˜ + sin α e ˜ 3 , w ˜ 3 = 2 2 ( e ˜ 1 e ˜ 2 ) , w ˜ i = e ˜ i , i = 4 , 5 , 6 , tg 2 α = 2 2 A 13 A 11 + A 12 A 33 .
Note that transversely isotropic materials according to the classification adopted in [15,17,31], can be of the following types: { 1 , 1 , 2 , 2 } , { 1 , 2 , 1 , 2 } , { 1 , 2 , 2 , 1 } , { 2 , 1 , 1 , 2 } , { 2 , 1 , 2 , 1 } , { 2 , 2 , 1 , 1 } .
Given Equation (51), the tensor in Equation (50) can be written in the form
A = a 2 C ( 1 ) + ( a 1 a 2 ) E + a 3 ( I ˜ e ˜ 3 + e ˜ 3 I ˜ ) + a 4 e ˜ 3 e ˜ 3 + + 1 2 a 5 ( e 3 I ˜ e 3 + e I e ˜ 3 e I + e I e 3 e I e 3 + e 3 e I e 3 e I ) ,
C ( 1 ) = I ˜ I ˜ , E = 1 2 ( C ( 2 ) + C ( 3 ) ) = e ˜ 1 e ˜ 1 + e ˜ 2 e ˜ 2 + e ˜ 4 e ˜ 4 , I ˜ = e ˜ 1 + e ˜ 2 , a 1 = 1 2 ( μ 1 sin 2 α + μ 2 cos 2 α + μ 3 ) , a 2 = 1 2 ( μ 1 sin 2 α + μ 2 cos 2 α μ 3 ) , a 3 = 2 2 ( μ 2 μ 1 ) sin α cos α , a 4 = μ 1 cos 2 α + μ 2 sin 2 , a 5 = μ 5 .
The dispersion tensor for Equation (52) after simple calculations is written as follows:
A ˜ = A 1 m m + A 2 I ˜ + A 3 ( m e 3 + e 3 m ) + A 4 e ˜ 3 ,
m = n I e I , A 1 = 1 ρ ( a 1 + 1 2 a 2 ) , A 2 = 1 ρ ( 1 2 a 2 m 2 + a 5 n 3 2 ) , A 3 = 1 ρ ( a 3 + a 5 ) n 3 , A 4 = 1 ρ ( a 4 n 3 2 + a 5 m 2 ) .
Then, it’s easy to find A ˜ 2 and A ˜ 3 . By virtue of Equation (54), we have
A ˜ 2 = B 1 m m + B 2 I ˜ + B 3 ( m e 3 + e 3 m ) + B 4 e ˜ 3 , A ˜ 3 = C 1 m m + C 2 I ˜ + C 3 ( m e 3 + e 3 m ) + C 4 e ˜ 3 ,
B 1 = A 1 2 m 2 + 2 A 1 A 2 + A 3 2 , B 2 = A 2 2 , B 3 = A 1 A 3 m 2 + A 2 A 3 + A 3 A 4 , B 4 = A 3 2 m 2 + A 4 2 , C 1 = A 1 3 m 4 + ( 3 A 1 2 A 2 + 2 A 1 A 3 2 ) m 2 + 3 A 1 A 2 2 + 2 A 2 A 3 2 + A 3 2 A 4 , C 3 = A 1 2 A 3 m 4 + ( 2 A 1 A 2 A 3 + A 1 A 3 A 4 + A 3 3 ) m 2 + A 2 2 A 3 + A 2 A 3 A 4 + A 3 A 4 2 , C 2 = A 2 3 , C 4 = A 1 A 3 2 m 4 + ( 2 A 3 2 A 4 + A 2 A 3 2 ) m 2 + A 4 3 .
From Equations (54) and (56), we find I 1 ( A ˜ k ) , k = 1 , 2 , 3 , and, by virtue of Equation (18), we get I m ( A ˜ ) , m = 2 , 3 . We have
I 1 ( A ˜ ) = A 1 m 2 + 2 A 2 + A 4 , I 1 ( A ˜ 2 ) = B 1 m 2 + 2 B 2 + B 4 , I 1 ( A ˜ 3 ) = C 1 m 2 + 2 C 2 + C 4 , I 2 ( A ˜ ) = 1 2 { A 1 2 m 4 + [ 2 A 1 ( 2 A 2 + A 4 ) B 1 ] m 2 + ( 2 A 2 + A 4 ) 2 2 B 2 B 4 } , I 3 ( A ˜ ) = 1 3 ! { A 1 3 m 6 + 3 [ A 1 2 ( 2 A 2 + A 4 ) A 1 B 1 ] m 4 + [ 3 A 1 ( 2 A 2 + A 4 ) 2 3 A 1 ( 2 B 2 + B 4 ) 3 B 1 ( 2 A 2 + A 4 ) + 2 C 1 ] m 2 + ( 2 A 2 + A 4 ) 3 3 ( 2 A 2 + A 4 ) ( 2 B 2 + B 4 ) + 2 ( 2 C 2 + C 4 ) } .
Knowing I m ( A ˜ ) , m = 1 , 2 , 3 (see Equation (58)), from the dispersion equation for any direction n , we find the eigenvalues A ˜ , and then the wave velocities.
Let us find, for example, the wave velocities in the following directions: n i = δ i 1 , n i = δ i 2 and n i = δ i 3 , i = 1 , 2 , 3 . Note that from the expression for I m ( A ˜ ) , m = 1 , 2 , 3 (see Equation (58)) follows that for directions n i = δ i 1 and n i = δ i 2 , i = 1 , 2 , 3 , they will accept the same values. In fact, by virtue of the corresponding Equations (55) and (58) for these directions, we have
I 1 ( A ˜ ) = a 1 + a 2 ρ + a 2 2 ρ + a 5 ρ , I 2 ( A ˜ ) = a 1 + a 2 ρ a 2 2 ρ + a 1 + a 2 ρ a 5 ρ + a 2 2 ρ a 5 ρ , I 3 ( A ˜ ) = a 1 + a 2 ρ a 2 2 ρ a 5 ρ .
From Equation (58), we see that the characteristic equation of the tensor A ˜ has roots
μ 1 = a 1 + a 2 ρ , μ 2 = a 2 2 ρ , μ 3 = a 5 ρ .
Hence, for the wave velocities in the initially resting medium, we have the following values:
c 1 = a 1 + a 2 ρ , c 2 = a 2 2 ρ , c 3 = a 5 ρ .
For the direction n i = δ i 3 , i = 1 , 2 , 3 , similar to Equations (59)–(61), the relations have the form
I 1 ( A ˜ ) = a 4 ρ + 2 a 5 ρ , I 2 ( A ˜ ) = 2 a 4 ρ a 5 ρ + a 5 ρ 2 , I 3 ( A ˜ ) = a 4 ρ a 5 ρ 2 , μ 1 = a 4 ρ , μ 2 = μ 3 = a 5 ρ , c 1 = a 4 ρ , c 2 = a 5 ρ .
Note that the consideration of the material with the symbol of the structure { 1 , 1 , 2 , 2 } of the trigonal syngonies (six essential components) reduces to the previous case.

5.4. Micropolar Material with a Center of Symmetry and the Anisotropy Symbol { 1 , 5 , 3 }

In this case, the properties of the medium are characterized by two tensors of the fourth rank A and D , which have the same structure. Therefore, it suffices to consider one of these tensors, since all the relations obtained for one tensor are obtained in a completely analogous way for the other (see Equation (20)). Consider, for example, the tensor A . In this case, its canonical representation has the form
A = ( λ 1 λ 2 ) u ˜ 1 u ˜ 1 + λ 2 E + ( λ 7 λ 2 ) m = 7 9 u ˜ m u ˜ m , E = C ( 2 ) = m = 1 9 u ˜ m u ˜ m .
By virtue of Equation (62), it is easy to find I 1 ( A ˜ k ) , k = 1 , 2 , 3 , and then by virtue of Equation (18) get I m ( A ˜ ) , m = 2 , 3 . Knowing I m ( A ˜ ) , m = 1 , 2 , 3 , we can write the characteristic equation A ˜ and find its roots, and then determine the desired velocity of the waves. Formulas are cumbersome, but the main thing is that they can be obtained in an explicit form. Note that the case of the medium { 5 , 1 , 3 } is a similar consideration as { 1 , 5 , 3 } . Therefore, in the general case, we do not stop. Next, we consider the case of a micropolar isotropic elastic medium with a center of symmetry and the symbol of structure { 1 , 5 , 3 } . In this case, the tensors A and D have representations
A = a 1 C ( 1 ) + a 2 C ( 2 ) + a 1 C ( 3 ) , D = d 1 C ( 1 ) + d 2 C ( 2 ) + d 1 C ( 3 ) ,
where the coefficients a i , d i , i = 1 , 2 , are determined in terms of the eigenvalues of these tensors by the formulas:
a 1 = 1 3 ( λ 1 λ 2 ) , a 2 = 1 2 ( λ 2 + λ 7 ) , a 3 = 1 2 ( λ 2 λ 7 ) , d 1 = 1 3 ( μ 1 μ 2 ) , d 2 = 1 2 ( μ 2 + μ 7 ) , d 3 = 1 2 ( μ 2 μ 7 ) .
It is easy to see that in the present case the dispersion tensors A ˜ and D ˜ , corresponding to tensors A and D of Equation (63), respectively, have the form
A ˜ = 1 ρ [ ( a 1 + a 3 ) n n + a 2 E ˜ ] , D ˜ = 1 J [ ( d 1 + d 3 ) n n + d 2 E ˜ ] .
From Equation (64), we see that A ˜ and D ˜ are similar to the tensors considered above from Equation (25). Introducing the notation
a = a 2 ρ , b = a 1 + a 3 ρ , f = d 2 J , g = d 1 + d 3 J ,
the tensors in Equation (64) can be written as follows:
A ˜ = a E ˜ + b n n , D ˜ = f E ˜ + g n n .
Then, similar to Equations (32) and (33), the roots of the dispersion equations of the tensors A ˜ and D ˜ and the wave velocities in the initially at rest media are determined by the formulas
η 1 = η 2 = a = λ 2 + λ 7 2 ρ , η 3 = a + b = λ 1 + 2 λ 2 3 ρ = λ + 2 μ ρ , η 4 = η 5 = f = μ 2 + μ 7 2 ρ , η 6 = f + g = μ 1 + 2 μ 2 3 ρ ,
c 1 = λ 2 + λ 7 2 ρ , c 2 = λ 1 + 2 λ 2 3 ρ , c 3 = μ 2 + μ 7 2 ρ , c 4 = μ 1 + 2 μ 2 3 ρ .
Thus, in a micropolar isotropic, infinite, initially resting medium, four waves appear in each direction, the velocities of which are calculated from Equation (67).
Note that, to determine the direction of propagation of these waves, it is necessary to find a complete system of eigenvectors and give a canonical representation of the dispersion tensors, and then conduct a study just as it was done above in the case of classical isotropic material. Of course, such a study can always be carried out for any anisotropic medium, as well as for a more general medium.
Note that, in the traditional form, the wave velocities in classical and micropolar media are researched in [33].

6. Conclusions

The statement of the eigenvalue problem for TBM of any order and any even rank is formulated. We consider materials with the anisotropy symbols { 1 . 5 } and { 5 . 1 } , as well as isotropic materials, and for them we determine the expressions for the velocities of wave propagation. In addition, we obtain expressions for the velocities of wave propagation for materials of cubic syngony with the anisotropy symbol { 1 , 2 , 3 } , hexagonal system (transversal isotropy) with anisotropy symbol { 1 , 1 , 2 , 2 } , trigonal system with anisotropy symbol { 1 , 1 , 2 , 2 } , and tetragonal system with anisotropy symbol { 1 , 1 , 1 , 2 , 1 } . We also obtain the expressions for the velocities of wave propagation for a micropolar medium with the anisotropy symbols { 1 . 5 . 3 } and { 5 . 1 . 3 } , and for an isotropic micropolar material.
Some results of the article were presented at the 6th International Conference on Mathematical Modeling in Physical Sciences (6th ICMSQUARE, August 28–31, 2017, Pafos, Cyprus) and published in [34].

Author Contributions

M.N. wrote the most part of the paper. A.U. took part in obtaining the relations for the tensors and their invariants for the different anisotropic materials, also in typing and presenting this work, and in translating the paper into English. S.L. and H.M. took part in the discussion about the applicability of the proposed method in many issues of mathematics and physics. Authors together discussed the problems in the study of which this method should be used.

Acknowledgments

This work was supported by the Russian Foundation for Basic Research (project no. 19-01-00016-A).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Poruchikov, V.B. Methods of the Dynamic Theory of Elasticity; Nauka: Moscow, Russia, 1986. (In Russian) [Google Scholar]
  2. Strett, D.V. Teoriya Zvuka; Book on Demand Ltd.: Moscow, Russia, 1955; pp. 1–2. (In Russian) [Google Scholar]
  3. Timoshenko, S.P. Kolebaniya v Inzhenernom Dele; Fizmatgiz: Moscow, Russia, 1959. (In Russian) [Google Scholar]
  4. Uflyand, Y.S. Rasprostraneniye voln pri poperechnykh kolebaniyakh sterzhney i plastin. Prikl. Matem. Mekhan. 1948, 12, 3. (In Russian) [Google Scholar]
  5. Novozhilov, V.V.; Slepyan, L.I. On Saint-Venant’s principle in the dynamics of beams. J. Appl. Math. Mech. 1965, 29, 293–315. [Google Scholar] [CrossRef]
  6. Slepyan, L.I. Non-Stationary Elastic Waves; Sudostroenie: Saint-Petersburg, Russia, 1972. (In Russian) [Google Scholar]
  7. Rosi, G.; Placidi, L.; dell’Isola, F. “First” and “slow” pressure waves electrically induced by nonlinear coupling in Biot-type porous medium saturated by a nematic liquid crystal. Zeitschrift für Angewandte Mathematik und Physik 2017, 68, 51. [Google Scholar] [CrossRef]
  8. dell’Isola, F.; Madeo, A.; Placidi, L. Linear plane wave propagation and normal transmission and reflection at discontinuity surfaces in second gradient 3D continua. Zeitschrift für Angewandte Mathematik und Mechanik 2012, 92, 52–71. [Google Scholar] [CrossRef]
  9. Rosi, G.; Auffrey, N. Anisotropic and dispersive wave propagation within strain-gradient framework. Wave Motion 2016, 63, 120–134. [Google Scholar] [CrossRef]
  10. Sobolev, S.L. Nekotoryye voprosy teorii rasprostraneniya kolebaniy. In Differentsialnyye i Integralnyye Uravneniya Matematicheskoy Fiziki; Frank, F., Mises, R., Eds.; ONTI: Leningrad, Moscow, 1937. (In Russian) [Google Scholar]
  11. Lamb, H. On the propagation of tremors over the surface on an elastic solid. Proc. R. Soc. Lond. A 1903, 72, 128–130. [Google Scholar] [CrossRef]
  12. Nikabadze, M.U. Development of the Method of Orthogonal Polynomials in the Classical and Micropolar Mechanics of Elastic Thin Bodies; Moscow University Press: Moscow, Russia, 2014. (In Russian) [Google Scholar]
  13. Nikabadze, M.U.; Ulukhanyan, A.R. Some Applications of Eigenvalue Problems for Tensor and Tensor Block Matrices for Mathematical Modeling of Micropolar Thin Bodies. Math. Comput. Appl. 2019, 24, 33. [Google Scholar] [CrossRef]
  14. Nikabadze, M.U.; Ulukhanyan, A.R. Modeling of multilayer thin bodies. Contin. Mech. Thermodyn. 2019, 31, 1–26. [Google Scholar] [CrossRef]
  15. Nikabadze, M.U. Topics on tensor calculus with applications to mechanics. J. Math. Sci. 2017, 225, 1–194. [Google Scholar] [CrossRef]
  16. Luré, A.I. Nonlinear Elasticity; Nauka: Moscow, Russia, 1980. (In Russian) [Google Scholar]
  17. Nikabadze, M.U. On Some Issues of Tensor Calculus with Applications to Mechanics. Contemp. Math. Fundam. Direct. 2015, 55, 3–194. [Google Scholar]
  18. Pobedrya, B.E. Numerical Methods in the Theory of Elasticity and Plasticity, 2nd ed.; Moscow State University: Moskva, Russia, 1995; 366p. (In Russian) [Google Scholar]
  19. Vekua, I.N. Fundamentals of Tensor Analysis and Covariant Theory; Nauka: Moscow, Russia, 1978. (In Russian) [Google Scholar]
  20. Pobedrya, B.E. Lectures in Tensor Analysis; Moscow State University: Moskva, Russia, 1986. (In Russian) [Google Scholar]
  21. Petrashen, G.I. Basics of the mathematical theory of the propagation of elastic waves. In Issues of Dynamic Propagation of Seismic Waves; Nauka: Moscow, Russia, 1978; 248p. (In Russian) [Google Scholar]
  22. Petrashen, G.I. Propagation of Waves in Anisotropic Elastic Media; Nauka: Leningrad, Russia, 1980; 280p. (In Russian) [Google Scholar]
  23. Sagomonyan, A.Y. Stress Waves in Continuous Media; Moscow University Press: Moscow, Russia, 1985; 416p. (In Russian) [Google Scholar]
  24. Eringen, A.C. Microcontinuum Field Theories I. Foundations and Solids; Springer: New York, NY, USA, 1999; 340p. [Google Scholar]
  25. Kupradze, V.D.; Gegelia, T.G.; Basheleishvili, M.O.; Burchuladze, T.V. Three-Dimensional Problems of the Mathematical Theory of Elasticity and Thermoelasticity; Nauka: Moscow, Russia, 1976; 664p. (In Russian) [Google Scholar]
  26. Nowacki, W. Theory of Elasticity; Mir: Moscow, Russia, 1975; 872p. (In Russian) [Google Scholar]
  27. Nikabadze, M.U. Some issues concerning a version of the theory of thin solids based on expansions in a system of Chebyshev polynomials of the second kind. Mech. Solids 2007, 42, 391–421. [Google Scholar] [CrossRef]
  28. Nikabadze, M.U.; Ulukhanyan, A.R. Analytical Solutions in the Theory of Thin Bodies. In Generalized Continua as Models for Classical and Advanced Materials, Advanced Structured Materials; Altenbach, H., Forest, S., Eds.; Springer: Cham, Switzerland, 2016; Volume 42, pp. 319–361. [Google Scholar]
  29. Baskakov, V.A.; Bestuzheva, N.P.; Konchakova, N.A. Linear Dynamic Theory of Thermoelastic Medium with Microstructure; Voronezh University Press: Voronezh, Russia, 2001; 162p. [Google Scholar]
  30. Nikabadze, M.U. Construction of Eigentensor Columns in the Linear Micropolar Theory of Elasticity. Mosc. Univ. Mech. Bull. 2014, 69, 1–9. [Google Scholar] [CrossRef]
  31. Nikabadze, M.U. Eigenvalue Problems of a Tensor and a Tensor-Block Matrix (TMB) of Any Even Rank with Some Applications in Mechanics. In Generalized Continua as Models for Classical and Advanced Materials, Advanced Structured Materials; Altenbach, H., Forest, S., Eds.; Springer: Cham, Switzerland, 2016; pp. 279–317. [Google Scholar]
  32. Nikabadze, M.U. To the problem of decomposition of the initial boundary value problems in mechanics. J. Phys. Conf. Ser. 2017, 739, 012056. [Google Scholar] [CrossRef]
  33. Ulukhanyan, A.R. The Dispersive Equations and Velocities of Propagation of Waves in the Couple Theory. Available online: https://www.researchgate.net/publication/318300021_Dispersionnye_uravnenia_i_skorosti_rasprostranenia_voln_v_momentnoj_teorii_The_Dispersive_Equations_and_Velocities_of_Propagation_of_Waves_in_the_Couple_Theory (accessed on 10 April 2019). (In Russian).
  34. Matevossian, H.; Nikabadze, M.; Ulukhanian, A. Determination of velocities of wave propagation in some media through the eigenvalues of the material tensors. J. Phys. Conf. Ser. 2018, 1141, 012154. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Nikabadze, M.U.; Lurie, S.A.; Matevossian, H.A.; Ulukhanyan, A.R. On Determination of Wave Velocities through the Eigenvalues of Material Objects. Math. Comput. Appl. 2019, 24, 39. https://0-doi-org.brum.beds.ac.uk/10.3390/mca24020039

AMA Style

Nikabadze MU, Lurie SA, Matevossian HA, Ulukhanyan AR. On Determination of Wave Velocities through the Eigenvalues of Material Objects. Mathematical and Computational Applications. 2019; 24(2):39. https://0-doi-org.brum.beds.ac.uk/10.3390/mca24020039

Chicago/Turabian Style

Nikabadze, Mikhail U., Sergey A. Lurie, Hovik A. Matevossian, and Armine R. Ulukhanyan. 2019. "On Determination of Wave Velocities through the Eigenvalues of Material Objects" Mathematical and Computational Applications 24, no. 2: 39. https://0-doi-org.brum.beds.ac.uk/10.3390/mca24020039

Article Metrics

Back to TopTop