The problems of the propagation of viscoelastic waves in extended laminate and layers of variable thickness are considered in    . In these papers, the boundaries in the formation of the structure of the wave field, as a spectrum of Eigen frequencies and Eigen modes, are revealed in a series for simple problems (for unchanged boundaries), and the boundary change is consistently accompanied by increasing difficulties. We also consider the occurrence of local singularities in wave fields.
In this paper, in contrast to the above, the propagation of waves along the z axis in an infinite viscoelastic cylinder with a radial crack is considered, which is a wedge with some angle ..
2. Statement of the Problem and Methods of Solution
The basic equations of motion of a deformable cylinder (with a radius R) with a longitudinal crack, which at , Case describes a wedge. They are given with three groups of relations. The system of equations of motion of a wedge in a cylindrical coordinate system ( ) it takes the form
Central differences (10), (11) are used for internal straight lines ( ), the left and right differences (10), (11) allow us to take into account the boundary conditions for φ. In the first case, the derivative of φ on the right sides of Equation (8) can be expressed by the formulas  :
Boundary conditions at is taken into account in the equations cor-
responding to the straight lines .
For the main unknowns that do not enter into the boundary conditions, wi, vi, ui use the right difference (10):
For a variables using central difference
The first and third of the conditions (9,c) are taken into account when approximating the derivatives of the function sj, B by j
Similarly, derivatives for the straight line with the number i = N, taking into
account the boundary conditions at . The only difference is the replace-
ment of the right finite difference left:
The number of straight lines can be reduced by using the ant symmetry conditions for the transverse vibrations of a plate at
The corresponding difference relations that take the conditions (16) can be written in the form
The resolving system of ordinary differential equations according to (8) has the form:
In Equation (18), the expressions for the derivatives are chosen from the relations (13)-(17), depending on the boundary conditions with respect to the coordinate φ.
The free surface conditions, the equivalent conditions (9,a) and forming the boundary value problem together with the Equation (18), is obtained in the form
Thus, the original spectral problem (8), (9) is reduced to the canonical problem (18) and (19) by discrimination of the coordinate φ by the method of lines, for the solution of which we apply the method of orthogonal sweep.
3. Numerical Results
The results of calculations are obtained on dimensionless parameters , , , A = 0.048, β = 0.05, α = 0.1, where is the wavelength of the wave, the density ( )   . For numerical realization of the problem a software tool МAPLE 9.5.
Comparisons of the obtained numerical results (without viscosity Math_84#) with the known results   are given in Table 1.
The limiting values of the phase velocity of the first edge mode, depending on the wedge angle at the vertex, found for a material with a Poisson’s ratio υ = 0.25 and , are given in Table 1 of Column 2. Also the results obtained in  for a wedge based on plates Variable thickness are given in Table 1 of column 3 (according to the theory of Kirchhoff-Love plates) and column 4 (according to the theory of Timoshenko’s plates). Column 5 corresponds to a calculation variant with three internal straight lines (N = 3) and boundary conditions (9,a), column 6 corresponds to the boundary conditions:
with the same number of straight lines. From the results, the variants of calculating variable-section plates (according to the Kirchhoff-Love and Tymoshenko theory) and the three-dimensional theory (the proposed method) show that they agree within 7% for wedge angles not exceeding 28˚ (wedge angle ).
In the framework of the work described in this paper (based on the calculation procedure for a three-dimensional wedge), taking into account the viscous properties of the material , they are shown in Figure 1.
Таble 1. Comparisons of the obtained numerical results.
Figure 1. Change in real and imaginary parts of the phase velocity Fluctuations depending onк.
In Figure 1 shown the real (curves CR1, CR2) and imaginary (curves CI1, CI2) parts of the dispersion curves of the first two modes of oscillations in an infinite viscous elastic wedge with angle ( ).
For comparison, the dependence on the wave number of the phase velocity of the first flexural mode of oscillations of a solid cylinder, founded by Pochhomer and Krie, using special functions (curve CR3)   and curve CI3  is shown in the same figure. We note the characteristic features of the curve CR3: at zero, the phase velocity is zero, and at infinity, tends to the Rayleigh wave velocity for a half-space. In the case of a viscoelastic wedge, the first mode has a cutoff frequency, and the phase velocity tends to infinity. At large wave numbers, the limiting phase velocity of this mode also coincides with the velocity of the Rayleigh wave.
Thus, unlike waveguides with a rectangular cross section in wedge-shaped waveguides with a sufficiently small wedge angle, in the analysis of the dispersion dependences of the first mode, it is permissible to use the theory of plates of variable cross section Kirchhoff-Love and Timoshenko. The established fact is explained by the phenomenon of localization of the shape of oscillations near the acute angle of the wedge, described in  .
Based on the results obtained, the following conclusions are drawn:
・ the results of calculating the limiting real velocity part (c = CR) of the propagation of the first mode of a wedge-shaped waveguide in the theory of changing the Kirchhoff-Love plates   and according to the dynamic theory of elasticity differ by no more than 6% for corners of the wedge top not exceeding 28˚;
・ in the wedge-shaped waveguides, the first mode has a cut-off frequency, and the phase velocity tends to infinity.
The work was supported by the Fund for Fundamental Research F-4-14 of the Republic of Uzbekistan.
 Safarov, I.I., Akhmedov, M.Sh. and Boltaev, Z.I. (2016) Ducting in Extended Plates of Variable Thickness. Global Journal of Science Frontier Research: F Mathematics & Decision Sciences, 16, 33-66.
 Bozorov, M.B., Safarov, I.I. and Shokin, Yu.I. (1966) Numerical Simulation of Oscillations of Dissipatively Homogeneous and Inhomogeneous Mechanical Systems. Siberian Branch of the Russian Academy of Sciences, Novosibirsk, 188 p.
 Safarov, I.I., Akhmedov, M.Sh. and Boltaev, Z.I. (2015) Dissemination Sinusoidal Waves in of a Viscoelastic Strip. Global Journal of Science Frontier Research: F Mathematics & Decision Sciences, 15, 39-60.