Distribution Free Waves in Viscoelastic Wedge with an Arbitrary Angle Tops

Author(s)
Safarov Ismail Ibrahimovich^{1},
Teshaev Mukhsin Khudoyberdiyevich^{1},
Akhmedov Maqsud Sharipovich^{2},
Boltayev Zafar Ixtiyorovich^{1}

Affiliation(s)

^{1}
Bukhara Technological Institute of Engineering, Bukhara, Republic of Uzbekistan.

^{2}
Institute of Engineering, Bukhara, Republic of Uzbekistan.

Abstract

In work questions of distribution of waves in a viscoelastic wedge with any corner of top is considered. The elastic cylinder with a radial crack is a wedge corner. The regional task for system of the differential equations in private derivatives is decided by means of a method of straight lines that allows using a method of orthogonal prorace.

In work questions of distribution of waves in a viscoelastic wedge with any corner of top is considered. The elastic cylinder with a radial crack is a wedge corner. The regional task for system of the differential equations in private derivatives is decided by means of a method of straight lines that allows using a method of orthogonal prorace.

Keywords

The Cylinder, The Differential Equations, The Orthogonal Prorace, Approximating Formulas, A Wedge of Corner, A Wave Guide, Dispersive Dependence

The Cylinder, The Differential Equations, The Orthogonal Prorace, Approximating Formulas, A Wedge of Corner, A Wave Guide, Dispersive Dependence

1. Introduction

The problems of the propagation of viscoelastic waves in extended laminate and layers of variable thickness are considered in [1] [2] [3] . 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 $\left|\phi \right|<{180}^{\circ}$ ..

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 $\phi =\left|{\phi}_{0}\right|<{180}^{\circ}$ , 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 ( $r,\phi ,z$ ) it takes the form

$\begin{array}{l}\rho \frac{\partial {u}_{r}}{d{t}^{2}}=\frac{\partial {\sigma}_{rr}}{dr}+\frac{{\sigma}_{rr}-{\sigma}_{r\phi}}{r}+\frac{1}{r}\frac{\partial {\sigma}_{r\phi}}{\partial \phi}+\frac{\partial {\sigma}_{rz}}{\partial z};\\ \rho \frac{{\partial}^{2}{u}_{\phi}}{\partial {t}^{2}}=\frac{1}{r}\frac{\partial {\sigma}_{\phi \phi}}{\partial \phi}+\frac{2{\sigma}_{r\phi}}{r}+\frac{\partial {\sigma}_{r\phi}}{\partial r}+\frac{\partial {\sigma}_{z\phi}}{\partial z},\\ \rho \frac{{\partial}^{2}{u}_{z}}{\partial {t}^{2}}=\frac{\partial {\sigma}_{zz}}{\partial z}+\frac{\partial {\sigma}_{rz}}{\partial r}+\frac{{\sigma}_{zz}}{r}+\frac{1}{r}\frac{\partial {\sigma}_{z\phi}}{\partial \phi}.\end{array}$ (1)

Here

$\begin{array}{l}{\epsilon}_{rr}=\frac{\partial {u}_{r}}{\partial r};\text{\hspace{0.17em}}{\epsilon}_{zz}=\frac{\partial {u}_{z}}{\partial z};\text{\hspace{0.17em}}{\epsilon}_{\phi \phi}=\frac{1}{r}\frac{\partial {u}_{r}}{\partial \phi}+\frac{{u}_{r}}{r};\\ \\ \})\text{T},\text{\hspace{0.17em}}i=\stackrel{\xaf}{1,N}\end{array}$ (11)

Central differences (10), (11) are used for internal straight lines ( $1<i<N$ ), 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 [4] :

$1<i<N$

$\begin{array}{l}{w}_{i,\phi}=\left({w}_{i+1}-{w}_{i-1}\right)/2\Delta ;\\ {u}_{i,\phi}=\left({u}_{i+1}-{u}_{i-1}\right)/2\Delta ;\\ {v}_{i,\phi}=\left({v}_{i+1}-{v}_{i-1}\right)/2\Delta ;\end{array}$ (12)

$\begin{array}{l}{\tau}_{\phi i,\phi}=\left({\tau}_{\phi \left(i+1\right)}-{\tau}_{\phi \left(i-1\right)}\right)/2\Delta ;\\ {\sigma}_{\phi i,\phi}=a\left({\sigma}_{i+1}-{\sigma}_{i-1}\right)/2\Delta +\frac{b}{r}\left[\left({v}_{i+1}-2{v}_{i}+{v}_{i-1}\right)/{\Delta}^{2}+{w}_{i,\phi}\right]+ck{u}_{i,\phi};\end{array}$

${B}_{i}=\left({u}_{i+1}-2{u}_{i}+{u}_{i-1}\right)/{\Delta}^{2}/k-k{v}_{i,\phi}.$

Boundary conditions at $\phi =-\frac{{\phi}_{0}}{2}$ is taken into account in the equations cor-

responding to the straight lines $i=1$ .

For the main unknowns that do not enter into the boundary conditions, w_{i}, v_{i}_{,} u_{i} use the right difference (10):

$i=1:$

$\begin{array}{l}{w}_{i,\phi}=\left(-3{w}_{1}+4{w}_{2}-{w}_{3}\right)/2\Delta ;\\ {v}_{i,\phi}=\left(-3{v}_{1}+4{v}_{2}-{v}_{3}\right)/2\Delta ;\\ {u}_{i,\phi}=\left(-3{u}_{1}+4{u}_{2}-{u}_{3}\right)/2\Delta .\end{array}$ (13)

For a variables ${\tau}_{\phi}$ using central difference

${\tau}_{{\phi}_{i},\phi}\approx \left({\tau}_{{\phi}_{0}}-{\tau}_{{\phi}_{2}}\right)/2\Delta =-{\tau}_{{\phi}_{2}}/2\Delta $ (14)

The first and third of the conditions (9,c) are taken into account when approximating the derivatives of the function s_{j}, B by j

$\begin{array}{c}{\sigma}_{{\phi}_{1},\phi}=\left({\sigma}_{{\phi}_{2}}-{\sigma}_{{\phi}_{0}}\right)/2\Delta ={\sigma}_{{\phi}_{2}}/2\Delta \\ =\left(a{\sigma}_{{r}_{2}}+\frac{b}{r}\left[\left({v}_{3}-{v}_{1}\right)/2\Delta +{w}_{2}\right]-ck{u}_{2}\right)/2\Delta \end{array}$

${B}_{1,\phi}\approx \left({B}_{2}-{B}_{0}\right)/2\Delta ={B}_{2}/2\Delta =\left[\left({u}_{3}-{u}_{1}\right)/2\Delta /r-k{v}_{2}\right]/2\Delta .$

Similarly, derivatives for the straight line with the number i = N, taking into

account the boundary conditions at $\phi =\frac{{\phi}_{0}}{2}$ . The only difference is the replace-

ment of the right finite difference left:

$i=N:$

$\begin{array}{l}{w}_{i,\phi}=\left(3{w}_{N}-4{w}_{N-1}+{w}_{N-2}\right)/2\Delta ;\\ {v}_{i,\phi}=\left(3{v}_{N}-4{v}_{N-1}+{v}_{N+1}\right)/2\Delta ;\\ {u}_{i,\phi}=\left(3{u}_{N}-4{u}_{N-1}+{u}_{N+1}\right)/2\Delta ;\end{array}$ (15)

$\begin{array}{l}{\tau}_{\phi \text{\hspace{0.17em}}i,\phi}=-{\tau}_{\phi \left(N-1\right)}/2\Delta ;\\ {\sigma}_{i,\phi}=-\left(a{\sigma}_{N-1}+\frac{b}{r}\left[\left({v}_{N}-{v}_{N-2}\right)/2\Delta +{w}_{N-1}\right]+ck{u}_{N-1}\right)/2\Delta =-\frac{{\sigma}_{N-1}}{2\Delta}\end{array}$

${B}_{i,\phi}=-\left[\left({u}_{N}-{u}_{N-2}\right)/2\Delta /r-k{v}_{N-1}\right]/2\Delta =-\frac{{B}_{N-1}}{2\Delta}.$

The number of straight lines can be reduced by using the ant symmetry conditions for the transverse vibrations of a plate at

$\phi =0:w=v={\sigma}_{\phi}=0.$ (16)

The corresponding difference relations that take the conditions (16) can be written in the form

$i=N:$

$\begin{array}{l}{w}_{i,\phi}=-{w}_{N-1}/2\Delta ;\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{u}_{i,\phi}=-{u}_{N-1}/2\Delta ;\\ {v}_{i,\phi}=\left(3{v}_{N}+4{v}_{N-1}-{v}_{N-2}\right)/2\Delta .\end{array}$ (17)

$\begin{array}{l}{\tau}_{{\phi}_{i},\phi}=\left(3{\tau}_{\phi N}-4{\tau}_{\phi \left(N-1\right)}+{\tau}_{\phi \left(N-2\right)}\right)/2\Delta ;\\ {\sigma}_{i,\phi}=-\left(a{\sigma}_{N-1}+\frac{b}{r}\left[\left({v}_{N}-{v}_{N-2}\right)/2\Delta +{w}_{N-1}\right]+ck{u}_{N-1}\right)/2\Delta =-\frac{{\sigma}_{N-1}}{2\Delta}\end{array}$

${B}_{i,\phi}=\left(-2{u}_{N}+{u}_{N-1}\right)/{\Delta}^{2}/r-k{v}_{i,\phi}.$

The resolving system of ordinary differential equations according to (8) has the form:

$\begin{array}{l}{{w}^{\prime}}_{i}={\sigma}_{i}/k-a\left(k{u}_{i}+\left({w}_{i}+{v}_{i,\phi}\right)/R\right);\\ {{v}^{\prime}}_{i}={\tau}_{{\phi}_{i}}+\left({v}_{i}-{w}_{i,\phi}\right)/R;\\ {{u}^{\prime}}_{i}={\tau}_{zi}+k{w}_{i};\\ {{\sigma}^{\prime}}_{i}=-{\omega}^{2}{w}_{i}+\left[2\left(\left({w}_{i}+{v}_{i,\phi}\right)/R-{{w}^{\prime}}_{i}\right)-{\tau}_{{\phi}_{i},\phi}\right]/r-k{\tau}_{zi};\\ {{\tau}^{\prime}}_{\phi i}=-{\omega}^{2}{v}_{i}+\left({\sigma}_{i,\phi}+2{\tau}_{\phi i}\right)/r-k\left({u}_{i,\phi}/R-k{v}_{i}\right);\\ {{\tau}^{\prime}}_{zi}=-{\omega}^{2}{u}_{i}+\left({B}_{i,\phi}+{\tau}_{zi}\right)/r+k\left({\sigma}_{i}+2\left(k{u}_{i}-{{w}^{\prime}}_{i}\right)\right).\end{array}$ (18)

In Equation (18), the expressions for the derivatives ${w}_{i,\phi},{v}_{i,\phi},{u}_{i,\phi},\text{\hspace{0.17em}}{\tau}_{\phi i,\phi}\text{\hspace{0.17em}},{\sigma}_{i,\phi},\text{\hspace{0.17em}}{B}_{i,\phi}$ 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

${\tau}_{\phi i}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\sigma}_{\phi i}=0\text{\hspace{0.17em}},\text{\hspace{0.17em}}{B}_{i}=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(i=1,\cdots ,N\right)$ (19)

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 $\nu =0.25$ , $\u0441/{c}_{s}$ , $R/\lambda $ , A = 0.048, β = 0.05, α = 0.1, where $\lambda =2\text{\pi}/k$ is the wavelength of the wave, the density $\rho =\rho /{\rho}_{0}$ ( $\lambda =\rho \left({\u0441}_{\u0440}^{2}-2{\u0441}_{s}^{2}\right),\mu ={c}_{s}^{2}\rho $ ) [9] [10] . For numerical realization of the problem a software tool МAPLE 9.5.

Comparisons of the obtained numerical results (without viscosity ${R}_{\lambda}\left(t\right)=$ Math_84#) with the known results [1] [3] 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 $r=R=1$ , ${R}_{\lambda}\left(t\right)={R}_{\mu}\left(t\right)=0$ are given in Table 1 of Column 2. Also the results obtained in [11] 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:

$\phi =-\frac{{\phi}_{0}}{2}:\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\sigma}_{\phi}={\sigma}_{\phi r}={\sigma}_{\phi z}=0;\text{\hspace{0.17em}}\phi =\frac{{\phi}_{0}}{2}:\text{\hspace{0.17em}}{u}_{r}={u}_{z}={\sigma}_{\phi \phi}=0$ (20)

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 ${\phi}_{0}\le {28}^{\circ}$ ).

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 ${R}_{\lambda}\left(t\right)={R}_{\mu}\left(t\right)=A{\text{e}}^{-\beta t}/{t}^{1-\alpha}$ , 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 C_{R}_{1}, C_{R}_{2}) and imaginary (curves C_{I}_{1}, C_{I}_{2}) parts of the dispersion curves of the first two modes of oscillations in an infinite viscous elastic wedge with angle (
${\phi}_{0}={45}^{\circ}$ ).

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 C_{R}_{3}) [12] [13] and curve C_{I}_{3} [14] is shown in the same figure. We note the characteristic features of the curve C_{R}_{3}: 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.

4. Conclusions

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 [3] .

Based on the results obtained, the following conclusions are drawn:

・ the results of calculating the limiting real velocity part (c = C_{R}) of the propagation of the first mode of a wedge-shaped waveguide in the theory of changing the Kirchhoff-Love plates [10] [11] 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.

Cite this paper

Ibrahimovich, S. , Khudoyberdiyevich, T. , Sharipovich, A. and Ixtiyorovich, B. (2017) Distribution Free Waves in Viscoelastic Wedge with an Arbitrary Angle Tops.*Applied Mathematics*, **8**, 736-745. doi: 10.4236/am.2017.85058.

Ibrahimovich, S. , Khudoyberdiyevich, T. , Sharipovich, A. and Ixtiyorovich, B. (2017) Distribution Free Waves in Viscoelastic Wedge with an Arbitrary Angle Tops.

References

[1] Safarov, I.I., Akhmedov, M.Sh. and Boltaev, Z.I. (2015) Distribution of the Natural Waves of Extended Lamellar Viscoelastic Bodies of Variable Thickness. LAP, Lambert Academic Publishing, 112 р.

[2] Safarov, I.I., Akhmedov, M.Sh. and Qilichov, O. (2016) Dynamics of Underground Hiheline from the Flowing Fluid. Lambert Academic Publishing, 345 р.

[3] Safarov, I.I., Teshaev, M.H. and Boltaev, Z.I. (2013) Mathematical Modeling of the Wave Process in a Mechanical Waveguide with Allowance for Internal Friction. LAP, Germany, 243 p.

[4] Godunov, S.K. (1977) On Numerical Solutions of Boundary Value Problems for Systems of Linear Ordinary Differential Equations. ELIBRARY.RU. Uspekhi Matematicheskikh Nauk, 1061, 16, 171-174.

[5] Ilyushin, A.A. and Pobedrya, B.E. (1970) Fundamentals of the Mathematical Theory of Thermal Viscoelasticity. Science, Moscow, 280 p.

[6] Filatov, A.R. and Sunchaliev, P.M. (1972) On Some Methods of Studying Nonlinear Problems in the Theory of Viscoelasticity. Reports of the Academy of Sciences of the USSR, 206, 201-203.

[7] Ewing, W.M., Jardetzku, W.S. and Press, F. (1962) Elastic Waves in Layered Media. McGraw-Hill, New York.

[8] Koltunov, M.A. (1976) Creep and Relaxation. Higher School, Moscow, 276 p.

[9] 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.

[10] 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.

[11] 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.

[12] Pochhammer, L. (1876) über die Fortpflanzung segechwindigkeiten schwingungen in einem unbergrawzten isotropen kreiscylinder. Y.reine und angew. Math., 81, 324-336.

[13] Chree, C. (1886) Longitudinal Vibrations of a Corcablar Bar. The Quarterly Journal of Pure and Applied Mathematics, 21, 287-298.

[14] Safarov, I.I., Akhmedov, M.Sh. and Boltaev, Z.I. (2016) Vibrations and Diffraction of Waves on a Cylindrical Body by a Viscoelastic Medium. Lambert Academic Publishing, Germany, 262 р.