# The Analytical Form of the Dispersion Equation of Elastic Waves in Periodically Inhomogeneous Medium of Different Classes of Crystals

Authors: Nurlybek A. Ispulov; Abdul Qadir; Marat K. Zhukenov; Talgat S. Dossanov; Tanat G. Kissikov
Journal: Advances in Mathematical Physics
 (2017)

---

## Abstract

The investigation of thermoelastic wave propagation in elastic media is bound to have much influence in the fields of material science, geophysics, seismology, and so on. The heat conduction equations and bound equations of motions differ by the difficulty level and presence of many physical and mechanical parameters in them. Therefore thermoelasticity is being extensively studied and developed. In this paper by using analytical matrizant method set of equation of motions in elastic media are reduced to equivalent set of first-order differential equations. Moreover, for given set of equations, the structure of fundamental solutions for the general case has been derived and also dispersion relations are obtained.


---

## Body

## 1. Introduction

The theory of thermoelasticity deals with the study of mutual interactions of thermal and mechanical fields in elastic bodies [1, 2]. It has vast applications in the various branches of Physics as well as in engineering, like materials engineering, mechanical engineering, nuclear engineering, and so forth. Theory of thermoelasticity is based on assumption that temperature distribution in an elastic object is governed by hyperbolic type parabolic-type partial differential equation as described by Fourier law of heat conduction [3–5]. According to Fourier law any thermal impulse is felt everywhere instantly in an object. Obviously it raised some serious concerns due to its unrealistic point of view. In order to circumvent this problem and to make it realistic a generalized theory of thermoelasticity was proposed which takes into account a finite thermal relaxation time. In this theory the temperature distribution is governed by hyperbolic type equations, which results in heat propagation in solids being considered as wave phenomenon instead of diffusion phenomenon.In order to investigate the wave propagation in anisotropic inhomogeneous medium a new method of matrizant was developed. This method allows investigation of wave propagation in anisotropic medium with various physical and mechanical properties [6–8].In 1950 Thompson [9] proposed a matrix method in order to investigate the elastic wave propagation in isotropic stratified media. Haskell also enhanced the method in 1953 [10]. After that major work was carried out by Stroh and others [11, 12]. He analytically investigated the dislocations in anisotropic medium by expressing first-order motion equations using (6×6) matrix. In order to investigate the insulators made up of piezoelectric materials, six-dimensional framework was enhanced to eight-dimensional formalism. Matrix method also paved the way for carrying out numerical simulation in anisotropic media [13, 14]. Various researchers have investigated the ordered structures and layered medium by using the matrix method. In this connection, following papers are of particular interest. Wave propagation has been investigated by matrix algebra method [15, 16], WKB method, and ray method [17–19]. Some investigations have tried to employ matricant, in which infinite product of truncated exponentials of the matrix of system coefficients and also Peano expansion are satisfied [7, 8]. However in case of periodic structure, Peano expansion cannot be fully solved. Therefore development of analytical techniques will open new dimensions to understand wave propagation in periodic structure.In case of layered and periodic medium the dispersion equations have been obtained and also the matricant structure was formed for nonhomogenous isotropic medium [3]. In [4] the matricant was obtained employing Chebyshev-Gegenbauer’s polynomial form, for the case of finite periodic inhomogeneous layer. For such structures, the modified conditions in determining the dispersion relationship having mutual transformation of transverse and longitudinal waves are obtained. In [20, 21] these results have been generalized in case of anisotropic inhomogeneous media.The applications of matrizant method for nondestructive testing and wave propagation in thermoelastic media are considered [22].Periodically heterogeneous media have lot of importance from applied and theoretical perspective. Wave propagation in discrete periodic structures has been extensively studied. While in case of continuous medium, layered homogeneous isotropic periodic structures are well studied. However the investigations of wave propagation in more complex periodically heterogeneous medium are carried out using various numerical methods or approximate analytical method. One of them is matrizant method; it was initially developed in late 70s in the Kazakh scientific school of Jakhan Suleimenuly Erzhanov, to investigate the dynamics of inhomogeneous medium.The method aims at reducing original equations of motion, by using method of separation of variables, to an equivalent system of ordinary differential equations of first order with variable coefficients. After that the resulting system of equations defines the structure of matrizant.


## 2. Elastic Waves

The motion equations in elastic medium and generalized Hooke’s law describe wave propagation in elastic media as follows [21]:(1)∂σij∂xj=ρ∂2ui∂t2σij=cijklεkl,where εkl=1/2∂ui/∂xj+∂uj/∂xi represents the deformation tensor components, ui denotes the mechanical displacement vector, σij are the stress tensor, cijkl represents the elastic parameters of nonisotropic media, and density of medium is represented by ρ.The medium is assumed to be stratified; that is, parameters employed to describe the material depend on space variable alongz-axis.Using the representation of the solution(2)fx,y,z,t=fzexp⁡iωt-ikxx-ikyy,where ω is radial frequency, ki are a projection of wave number. The multiplier exp⁡iωt-ikxx-ikyy is omitted in the following equations for clarity.Taking into consideration propagation direction, derivative of anisotropic medium onz-axis, and using (2) then (1) are reduced to a system of first-order ODEs having variable coefficients. (3)dW→dz=BW→;W→=uz,σzz,ux,σxz,uy,σyzt;the transposition operator is denoted by t.The coefficient matrixB, for the case of triclinic anisotropic medium, takes the following form: (4)B=0b12b13b14b15b16b2100b240b26b24b14b33b34b35b360b13b43b33b45b46b26b16b46b36b55b560b15b45b35b65b55.For the orthorhombic anisotropy:(5)B=0b12b130b150b2100b240b26b2400b340b360b13b430b450b2600b360b560b15b450b650.If the vector, representing the direction of wave propagation, is in(xz) plane of anisotropy orthorhombic medium, the coefficient matrix given in (5) is divided into 4.If the wave propagation direction vector lies in the media matrix (5) splits into (4×4) and (2×2) matrices: (6)B′=0b12b130b2100b24b2400b340b13b430;W→=uz,σzz,ux,σxzt;B′′=0b56b650;W→=uy,σyzt.


## 3. Matricant Structure and Its Implications

### 3.1. Matricant Method

In order to describe the wave propagation in elastic medium various analytical approaches are used, like the formalism proposed by Stroh–Barnett–Lothe [23] for piezoelasticity, the state vector “W” method [13]. The matricant approach is different from others analytical techniques used to investigate wave propagation in elastic medium. In the matricant approach the vector W of coefficient matrix B is chosen in (ui, σij), (θ, qi), pairs, for instance. The selection of pairs depends on the type of wave to be investigated. Coupled waves in the general case and along main crystal axes are suitably described by the use of such notations.Solutions of (5) are written as(7)W→z=Tz,z0W→z0Here, T(z,z0) is the matricant, that is, the normalized fundamental solution matrix of the systems of ODEs. For T(z,z0) and T-1(z,z0), there are representations in the form of the infinite matrix integral series of exponential type as follows [9, 10]:(8)Tz,z0=I+∫z0zBz1dz1+∫z0z∫z0z1Bz1Bz2dz1dz2+⋯(9)T-1z,z0=I-∫0zB1z1dz1+∫0z∫0z1B2z2B1z1dz1dz-⋯;the identity matrix is denoted byI.The expansion in (9) is the alternating-sign series with reverse argument ordering of the integrated product ofB(zi). Note that the matrixB(zi) does not commute. As the initial system of equations are satisfied by the matricant, so the successive approximation methods can be used to obtain (8).(10)dTdz=BT.It follows from substitution of (7) into (10).Similar to (10), inverse matricant T-1 is the solution of the equation(11)dT-1dz=-T-1B.It follows from differentiation that the identity(12)TT-1≡T-1T≡IThe solutions of (2×2) order matrix are well known.(13)B=0b12b210;the matricants have the structure(14)T=t11t12t21t22;T-1=t22-t12-t21t11.For unimodular matrices the above result can be carried out.In case of (4×4) coefficient matrix in (6), its inverse matrix is given by (15)T-1=t22-t12-t42t32-t21t11t41-t31-t24t14t44-t34t23-t13-t43t33.The result given in (15) defines the properties of solutions of systems of first-order ODEs with variable coefficients. This has been obtained with term to term comparison of elements of (8) and (9) and with the help of mathematical induction [5, 24].In the case of the coefficient matrix (5), the (6×6) matricant structure is obtained in the following form: (16)T-1=t22-t12-t42t32-t62t52-t21t11t41-t31t61-t51-t24t14t44-t34t64-t54t23-t13-t43t33-t63t53-t26t16t46-t36t66-t56t25-t15-t45t35-t65t55.It may be noted that elements ofT-1 contain only of elements from matricantT. It is a one to one correspondence of elements of direct and its inverse matricant given in (15) and (16). The conservation laws are contained by invariant relationships; these laws have to be satisfied in wave processes. In 1D inhomogeneous isotropic medium having various crystals symmetry is described by matricant structure. The matricant of coefficient matrix is of order (2n×2n) [5].


### 3.2. Periodic Structures

Suppose the variation in parameters that describes medium isfi(z). The conditionfi(z+h)=fi(z) is satisfied by periodic structure and also by parameters altered by environment. The period of inhomogeneity is denoted by h. Monodromy matrix is matricant of single period of inhomogeneity. It is known for single period inhomogeneity that W→(h)=T(0,h)W→0(0) and the Floquet–Bloch conditions W→z+h=W→0zexp⁡ikzh. From this follows the following equality: (17)T-Iexp⁡ikzhW→0=0;the equivalent condition are obtained by multiplying (17) with T-1eikh.(18)T-1-Iexp⁡-ikzhW→0=0.Using auxiliary matrix and (17) and (18), the modified conditions as given below are derived as follows:(19)p^-Icos⁡kzhW→0=0,p^=12T+T-1.Matrix p^ for cases (6) is given by(20)p^=p110p13p140p11p23p24-p24p14p220p23-p130p22.The dispersion in periodic structure is determined by characteristics of (19).(21)sdet⁡p^-Icos⁡kzh=0.The dispersion equation in case of (4×4) matrices is found by(22)
cos
⁡
k
1
z
h
=
p
~
1
cos
⁡
k
2
z
h
=
p
~
2
=
1
2
p
11
+
p
22
±
1
2
p
11
-
p
22
2
+
4
p
14
p
23
-
p
13
p
24
.The equation of dispersion is obtained. However, the order of characteristic equation is reduced by half when the conditions as laid down in (19) are imposed.Following recurrence relationship can be obtained from (19) as follows: (23)T2=2p^T-I.Matricant representing periodically inhomogeneous layers can be obtained applying (24).In the presence ofn periods of heterogeneity in form H=nh, we can obtain (24)TH=Tmnh=Pnp^Tm-Pn-1p^,where Tm=Th represents the monodromy matrix and Pnp^ denotes matrix polynomials of Chebyshev-Gegenbauer [5, 25]. The results of Brillouin and Parodi [26] are generalized by above equations.


### 3.3. Structuring the Matrizant

Structuring the matrizant (normalized matrix of fundamental solutions) is based on its representation in the form of the exponential matrix series [5, 7], (8) and (9).These matrix series converge absolutely and uniformly on any finite interval. In this case, the following relation is true, (12).For the matrizant the following expressions also hold:(i)
Tz0,z=Tz1,zTz0,z1.(ii)
ln⁡Tz0,z=∫z0zspB(z1)dz1.(iii)
IfB=B0 – constant matrix, then T0=exp⁡B0z-z0.(iv)
dT-1/dz=-T-1B.Matrix series (8) and (9) can be written in summation form(25)T=∑n=0∞Tn,T-1=∑n=0∞Tn-1.The indexn corresponds to the number of multiplied under signs of integral matrices, where B(zi) is the number of integrals of the matrix in each term of the series as given in (1). Moreover, the terms of the series with even and odd values of n are as follows:(26)Teven=∑n=0∞T2n,Todd=∑n=0∞T2n+1Teven-1=∑n=0∞T2n-1,Todd-1=∑n=0∞T2n+1-1,n=0,1,2,3,….In such case constructing the matrizant is basically expressing the relationships between the elements of theT and T-1 matrices, it is based on the element-wise comparison.As a first approximationtij(1)=-tij(1), tij(-1)=∫0zbij(z)dz.Elastic waves propagating in the orthorhombic syngony of the classes mm2 and 222, hexagonal syngony (6,6-, 622, 6mm, 6-m2), tetragonal syngony (class 422), and matrizant structure are built based on the structure of the coefficient matrix, based on the element-wise comparison of matrix T and T-1. The structure of the coefficient matrix in the bulk case is as follows [5]:(27)B=0b12b130b150b2100b240b26b2400b34000b13b430b450b260000b560b15b450b650.See (3).(28)W→=uz,σzz,ux,σxz,uy,σyzt.The coefficients of the matrixbij (27) for the crystals of orthorhombic syngony are equal to(29)b12=1c33;b13=b42=ikxc13c33;b15=b62=ikyc32c33;b21=-ρω2;b24=b31=ikx;b26=b51=iky;b34=1c55;b43=ky2c66-ρω2+kx2c11-c132c33;b45=b63=c66+c12-c13c32c33kxky;b56=1c44;b65=kx2c66-ρω2+c22-c232c33ky2,where cαβ- are elastic parameters, kx,ky- are components of the wave vector, ρ is medium density, and ω isangular frequency.Coefficientsbij for hexagonal syngony (classes 6, 6-, 622, 6mm, 6-m2) have the form(30)b12=1c33;b13=ikxc13c33;b15=ikyc13c33;b21=-ρω2;b24=ikx;b26=iky;b34=b56=1c44;b43=ky2c11-c122+kx2c11-c132c33-ρω2;b45=c12+c11-c122-c132c33kxky;b65=kx2c11-c122+c11-c132c33ky2-ρω2.Coefficientsbij for tetragonal syngony (class 422) have the form(31)b12=1c33;b13=ikxc13c33;b15=ikyc13c33;b21=-ρω2;b24=ikx;b26=iky;b34=b56=1c44;b43=ky2c66+kx2c11-c132c33-ρω2;b45=c12+c66-c132c33kxky;b65=kx2c66+-c132c33ky2-ρω2.As it can be seen from the last relations coefficientsbij of crystals of high and average symmetry differ only in the values ​​of the elastic constants cαβ.Matrix of order 6 describes the propagation of bound elastic one longitudinal and two transverse waves.The second approximation matrizant has the form(32)T2=∫0z∫0z1Bz1Bz2dz1dz2.Inverse matrizant in the second approximation takes the form(33)T2-1=∫0z∫0z1Bz2Bz1dz1dz2.Comparison of the terms of the second approximation gives the following relationship between the elements of matrizantT and T-1:(34)T2-1=t2200t320t520t11t410t6100t14t440t640t2300t330t530t16t460t660t2500t350t552.The elementstij are the elements of the direct matrizant (27).Similarly, elements of the third approximation are compared.(35)T3=∭Bz1Bz2Bz3dz1dz2dz3orT3=T2Bz3T3=0t12t130t150t2100t240t26t3100t340t360t42t430t450t5100t540t560t62t630t6503.Inverse matrizant in the third approximation has the form(36)T3-1=∭Bz3Bz2Bz1dz1dz2dz3orT3-1=Bz3T-12and has the structure(37)T3-1=0t12t420t620t2100t150t51t2400t340t540t13t430t630t2600t360t560t15t450t6503.The elementstij are the elements of the matrizant (35). Endless rows of the matrix can be written as follows [5]:(38)T=Teven+Todd,T-1=Teven-1-Todd-1,where T±1 – corresponds to the sum of odd and even rows (9, 10).Mathematical induction proves that the structure of theTn-1 is preserved for any n.Structure (34) is valid for all even T-1 and the structure of (37) is valid for all odd values of T-1. Generalizing (34) and (37) according to (38), we obtain the structure of the T-1 matrizant as follows [5]:(39)T-1=t22-t12-t42t32-t62t52-t21t11t41-t31t61-t51-t24t14t44-t34t64-t54t23-t13-t43t33-t63t53-t26t16t46-t36t66-t56t25-t15-t45t35-t65t55.Elementstij of matrices Teven-1 and Todd-1 are elements of matrices; Teven and Todd are the elements of direct matrix T, respectively.Thus, the structure of matrizant is a relationship between elements of the forward and inverse matrizant in the form (39), as well as the relationship between the elements T and T-1 as follows from (12).Analytical representation of matrizant of periodically inhomogeneous layer is derived from knowing the structureT-1.


## 3.1. Matricant Method

In order to describe the wave propagation in elastic medium various analytical approaches are used, like the formalism proposed by Stroh–Barnett–Lothe [23] for piezoelasticity, the state vector “W” method [13]. The matricant approach is different from others analytical techniques used to investigate wave propagation in elastic medium. In the matricant approach the vector W of coefficient matrix B is chosen in (ui, σij), (θ, qi), pairs, for instance. The selection of pairs depends on the type of wave to be investigated. Coupled waves in the general case and along main crystal axes are suitably described by the use of such notations.Solutions of (5) are written as(7)W→z=Tz,z0W→z0Here, T(z,z0) is the matricant, that is, the normalized fundamental solution matrix of the systems of ODEs. For T(z,z0) and T-1(z,z0), there are representations in the form of the infinite matrix integral series of exponential type as follows [9, 10]:(8)Tz,z0=I+∫z0zBz1dz1+∫z0z∫z0z1Bz1Bz2dz1dz2+⋯(9)T-1z,z0=I-∫0zB1z1dz1+∫0z∫0z1B2z2B1z1dz1dz-⋯;the identity matrix is denoted byI.The expansion in (9) is the alternating-sign series with reverse argument ordering of the integrated product ofB(zi). Note that the matrixB(zi) does not commute. As the initial system of equations are satisfied by the matricant, so the successive approximation methods can be used to obtain (8).(10)dTdz=BT.It follows from substitution of (7) into (10).Similar to (10), inverse matricant T-1 is the solution of the equation(11)dT-1dz=-T-1B.It follows from differentiation that the identity(12)TT-1≡T-1T≡IThe solutions of (2×2) order matrix are well known.(13)B=0b12b210;the matricants have the structure(14)T=t11t12t21t22;T-1=t22-t12-t21t11.For unimodular matrices the above result can be carried out.In case of (4×4) coefficient matrix in (6), its inverse matrix is given by (15)T-1=t22-t12-t42t32-t21t11t41-t31-t24t14t44-t34t23-t13-t43t33.The result given in (15) defines the properties of solutions of systems of first-order ODEs with variable coefficients. This has been obtained with term to term comparison of elements of (8) and (9) and with the help of mathematical induction [5, 24].In the case of the coefficient matrix (5), the (6×6) matricant structure is obtained in the following form: (16)T-1=t22-t12-t42t32-t62t52-t21t11t41-t31t61-t51-t24t14t44-t34t64-t54t23-t13-t43t33-t63t53-t26t16t46-t36t66-t56t25-t15-t45t35-t65t55.It may be noted that elements ofT-1 contain only of elements from matricantT. It is a one to one correspondence of elements of direct and its inverse matricant given in (15) and (16). The conservation laws are contained by invariant relationships; these laws have to be satisfied in wave processes. In 1D inhomogeneous isotropic medium having various crystals symmetry is described by matricant structure. The matricant of coefficient matrix is of order (2n×2n) [5].


## 3.2. Periodic Structures

Suppose the variation in parameters that describes medium isfi(z). The conditionfi(z+h)=fi(z) is satisfied by periodic structure and also by parameters altered by environment. The period of inhomogeneity is denoted by h. Monodromy matrix is matricant of single period of inhomogeneity. It is known for single period inhomogeneity that W→(h)=T(0,h)W→0(0) and the Floquet–Bloch conditions W→z+h=W→0zexp⁡ikzh. From this follows the following equality: (17)T-Iexp⁡ikzhW→0=0;the equivalent condition are obtained by multiplying (17) with T-1eikh.(18)T-1-Iexp⁡-ikzhW→0=0.Using auxiliary matrix and (17) and (18), the modified conditions as given below are derived as follows:(19)p^-Icos⁡kzhW→0=0,p^=12T+T-1.Matrix p^ for cases (6) is given by(20)p^=p110p13p140p11p23p24-p24p14p220p23-p130p22.The dispersion in periodic structure is determined by characteristics of (19).(21)sdet⁡p^-Icos⁡kzh=0.The dispersion equation in case of (4×4) matrices is found by(22)
cos
⁡
k
1
z
h
=
p
~
1
cos
⁡
k
2
z
h
=
p
~
2
=
1
2
p
11
+
p
22
±
1
2
p
11
-
p
22
2
+
4
p
14
p
23
-
p
13
p
24
.The equation of dispersion is obtained. However, the order of characteristic equation is reduced by half when the conditions as laid down in (19) are imposed.Following recurrence relationship can be obtained from (19) as follows: (23)T2=2p^T-I.Matricant representing periodically inhomogeneous layers can be obtained applying (24).In the presence ofn periods of heterogeneity in form H=nh, we can obtain (24)TH=Tmnh=Pnp^Tm-Pn-1p^,where Tm=Th represents the monodromy matrix and Pnp^ denotes matrix polynomials of Chebyshev-Gegenbauer [5, 25]. The results of Brillouin and Parodi [26] are generalized by above equations.


## 3.3. Structuring the Matrizant

Structuring the matrizant (normalized matrix of fundamental solutions) is based on its representation in the form of the exponential matrix series [5, 7], (8) and (9).These matrix series converge absolutely and uniformly on any finite interval. In this case, the following relation is true, (12).For the matrizant the following expressions also hold:(i)
Tz0,z=Tz1,zTz0,z1.(ii)
ln⁡Tz0,z=∫z0zspB(z1)dz1.(iii)
IfB=B0 – constant matrix, then T0=exp⁡B0z-z0.(iv)
dT-1/dz=-T-1B.Matrix series (8) and (9) can be written in summation form(25)T=∑n=0∞Tn,T-1=∑n=0∞Tn-1.The indexn corresponds to the number of multiplied under signs of integral matrices, where B(zi) is the number of integrals of the matrix in each term of the series as given in (1). Moreover, the terms of the series with even and odd values of n are as follows:(26)Teven=∑n=0∞T2n,Todd=∑n=0∞T2n+1Teven-1=∑n=0∞T2n-1,Todd-1=∑n=0∞T2n+1-1,n=0,1,2,3,….In such case constructing the matrizant is basically expressing the relationships between the elements of theT and T-1 matrices, it is based on the element-wise comparison.As a first approximationtij(1)=-tij(1), tij(-1)=∫0zbij(z)dz.Elastic waves propagating in the orthorhombic syngony of the classes mm2 and 222, hexagonal syngony (6,6-, 622, 6mm, 6-m2), tetragonal syngony (class 422), and matrizant structure are built based on the structure of the coefficient matrix, based on the element-wise comparison of matrix T and T-1. The structure of the coefficient matrix in the bulk case is as follows [5]:(27)B=0b12b130b150b2100b240b26b2400b34000b13b430b450b260000b560b15b450b650.See (3).(28)W→=uz,σzz,ux,σxz,uy,σyzt.The coefficients of the matrixbij (27) for the crystals of orthorhombic syngony are equal to(29)b12=1c33;b13=b42=ikxc13c33;b15=b62=ikyc32c33;b21=-ρω2;b24=b31=ikx;b26=b51=iky;b34=1c55;b43=ky2c66-ρω2+kx2c11-c132c33;b45=b63=c66+c12-c13c32c33kxky;b56=1c44;b65=kx2c66-ρω2+c22-c232c33ky2,where cαβ- are elastic parameters, kx,ky- are components of the wave vector, ρ is medium density, and ω isangular frequency.Coefficientsbij for hexagonal syngony (classes 6, 6-, 622, 6mm, 6-m2) have the form(30)b12=1c33;b13=ikxc13c33;b15=ikyc13c33;b21=-ρω2;b24=ikx;b26=iky;b34=b56=1c44;b43=ky2c11-c122+kx2c11-c132c33-ρω2;b45=c12+c11-c122-c132c33kxky;b65=kx2c11-c122+c11-c132c33ky2-ρω2.Coefficientsbij for tetragonal syngony (class 422) have the form(31)b12=1c33;b13=ikxc13c33;b15=ikyc13c33;b21=-ρω2;b24=ikx;b26=iky;b34=b56=1c44;b43=ky2c66+kx2c11-c132c33-ρω2;b45=c12+c66-c132c33kxky;b65=kx2c66+-c132c33ky2-ρω2.As it can be seen from the last relations coefficientsbij of crystals of high and average symmetry differ only in the values ​​of the elastic constants cαβ.Matrix of order 6 describes the propagation of bound elastic one longitudinal and two transverse waves.The second approximation matrizant has the form(32)T2=∫0z∫0z1Bz1Bz2dz1dz2.Inverse matrizant in the second approximation takes the form(33)T2-1=∫0z∫0z1Bz2Bz1dz1dz2.Comparison of the terms of the second approximation gives the following relationship between the elements of matrizantT and T-1:(34)T2-1=t2200t320t520t11t410t6100t14t440t640t2300t330t530t16t460t660t2500t350t552.The elementstij are the elements of the direct matrizant (27).Similarly, elements of the third approximation are compared.(35)T3=∭Bz1Bz2Bz3dz1dz2dz3orT3=T2Bz3T3=0t12t130t150t2100t240t26t3100t340t360t42t430t450t5100t540t560t62t630t6503.Inverse matrizant in the third approximation has the form(36)T3-1=∭Bz3Bz2Bz1dz1dz2dz3orT3-1=Bz3T-12and has the structure(37)T3-1=0t12t420t620t2100t150t51t2400t340t540t13t430t630t2600t360t560t15t450t6503.The elementstij are the elements of the matrizant (35). Endless rows of the matrix can be written as follows [5]:(38)T=Teven+Todd,T-1=Teven-1-Todd-1,where T±1 – corresponds to the sum of odd and even rows (9, 10).Mathematical induction proves that the structure of theTn-1 is preserved for any n.Structure (34) is valid for all even T-1 and the structure of (37) is valid for all odd values of T-1. Generalizing (34) and (37) according to (38), we obtain the structure of the T-1 matrizant as follows [5]:(39)T-1=t22-t12-t42t32-t62t52-t21t11t41-t31t61-t51-t24t14t44-t34t64-t54t23-t13-t43t33-t63t53-t26t16t46-t36t66-t56t25-t15-t45t35-t65t55.Elementstij of matrices Teven-1 and Todd-1 are elements of matrices; Teven and Todd are the elements of direct matrix T, respectively.Thus, the structure of matrizant is a relationship between elements of the forward and inverse matrizant in the form (39), as well as the relationship between the elements T and T-1 as follows from (12).Analytical representation of matrizant of periodically inhomogeneous layer is derived from knowing the structureT-1.


## 4. Dispersion Equations for the Elastic Anisotropic Mediums

We introduce the following matrix [5]:(40)p^=12T+T-1;substituting values of T and T-1 in (40), we obtain(41)2p^=p10p13p14p15p160p1p23p24p25p26-p24p14p30p35p36p23-p130p3p45p46-p26p16p46-p36p50p25-p15-p45p350p5.The matrixp as given in (40) which is important for the regular structures gives the recurrence relation as in (23).Consistent application of (23) gives possibility of representing Tn in the form (42)Tn=PnpT-Pn-1p,where Pn(p) is Chebyshev-Gegenbaur matrix polynomials of the second kind.Equation (42) allows obtaining Tn in an explicit analytic form(43)Tn=∑i=13PiPnp~iT-Pnp~iI,where Pi=1/p~i-p~jp~i-p~kp^-p~jIp^-p~kI, i, j,k-1,2,3;i≠j,j≠k,i≠k.p ~ 1 , p ~ 21 , p ~ 31are the roots of the characteristic equation, satisfying the following condition (23).The matrix method of matrizant allows twice lowering the degree of the characteristic equation, which in the end has the form(44)λ3+aλ2+bλ+c2=0λ3+aλ2+bλ+c=0,where (45)λ=cos⁡k~hia=-p1+p2+p5b=p1p3+p1p5+p3p5-p14p23+p13p24-p16p25+p15p26+p35p46c=-p1p3p5+p36p45-p35p46+p3p16p25-p15p26+p5p14p23-p13p24--p16p23p35+p24p45+p13p26p35-p25p36+p15p23p36+p24p46+p14p26p45-p25p46.The solution of the characteristic equation (44) gives three roots, which have the following form:(46)cos⁡k~h1=161.58δ3-2a+2.52a2-3bδ3cos⁡k~h2=112-3.17δ3-4a-4-23a2-3bδ3cos⁡k~h3=112-3.17-13δ3-4a+5a2-3bδ3,where δ=-2a3+9ab-27c+2a3-9ab+27c2-4a2-3b3.Relations (46) determine the dispersion equations of elastic waves in above-mentioned crystals. The difference lies in the roots of the difference between the coefficients bij in the matrix (27).


## 5. Conclusion

In this paper we have developed the structure of matrizant and from it obtained invariant relations which reflects the inner symmetry of inner equations and contains conservation laws. Also we have derived an analytical representation of matrizant of periodically inhomogeneous layer using Chebyshev-Gegenbauer polynomials and obtained a separate dispersion equation. Finally, by using different crystals systems the analytical solution of equations of motion for a wide class of homogenous anisotropic medium has been obtained.



---

*Source: 5236898-2017-01-29.xml*