Of Own and Forced Vibrations of Dissipative Inhomogeneous Mechanical Systems

Show more

1. Introduction

The use of damping vibrations of dynamic viscoelastic mechanical systems with different rheological properties, for all the studies of the problem [1] [2] , is rarely considered in the scientific literature. At the same time, modern machine building is characterized by a wide use of polymeric and metallic materials with various viscoelastic properties [3] [4] [5] [6] . A mechanical system, consisting of rigid and deformable bodies, connected to each other and to a base by deformable (elastic or viscoelastic) elements, is studied.

2. Statement of the Problem and Methods of Solution

Let us consider a composite structure occupying a volume $V={\displaystyle \underset{n=1}{\overset{N}{\sum}}{V}_{n}}$ , bounded

by a surface
$\sum ={\sum}_{u}+{\sum}_{p}$ . Each of N volumes
${V}_{n}$ is filled with a viscoelastic medium, the parameters of which depend on the number n. On a part of the surface
${\sum}_{u}$ ―zero motions, on
${\sum}_{p}$ ―zero stresses; mass forces are absent. The frequencies and damping indices of the natural vibrations of the structures are to be determined. The physical properties of the material of the n^{th} volume of the structure, are described by the relations:

${\sigma}_{ij}={\stackrel{\u02dc}{\lambda}}_{n}{\epsilon}_{kk}{\delta}_{ij}+2{\stackrel{\u02dc}{\mu}}_{n}{\epsilon}_{ij}\left(i,j,k=1,2,3;n=1,2,\cdots ,N\right)$ ,

where ${\sigma}_{ij},{\epsilon}_{ij}$ ―components of stress and strain tensors , ${\stackrel{\u02dc}{\lambda}}_{n},{\stackrel{\u02dc}{\mu}}_{n}$ ―Volterra operators

$\begin{array}{l}{\stackrel{\u02dc}{\lambda}}_{n}\phi \left(t\right)={\lambda}_{n}\left[\phi \left(t\right)-{\displaystyle \underset{0}{\overset{t}{\int}}{R}_{\lambda n}\left(t-\tau \right)\phi \left(\tau \right)\text{d}\tau}\right];\\ {\stackrel{\u02dc}{\mu}}_{n}\phi \left(t\right)={\mu}_{n}\left[\phi \left(t\right)-{\displaystyle \underset{0}{\overset{t}{\int}}{R}_{\mu n}\left(t-\tau \right)\phi \left(\tau \right)\text{d}\tau}\right]\end{array}$ (1)

here ${\lambda}_{n},{\mu}_{n},{R}_{\lambda n},{R}_{\mu n}$ ―Lame parameters and the relaxation core of the medium occupying the volume ${V}_{n}$ ; $\phi \left(t\right)$ ―an arbitrary function of time.

We assume, that the integral terms in (1) are small. Let the functions $\phi \left(t\right)$ has the form $\phi \left(t\right)=\psi \left(t\right){\text{e}}^{-i{\omega}_{R}t}$ , where $\psi \left(t\right)$ ―slowly varying function of time, ${\omega}_{R}$ ―real constants. Further, applying the freezing procedure, which has been widely studied [7] , to the place (1), we can obtain approximate relations:

$\begin{array}{l}{\stackrel{\u02dc}{\lambda}}_{n}\phi \approx {\stackrel{\xaf}{\lambda}}_{n}\phi ={\lambda}_{n}\left[1-{\Gamma}_{\lambda n}^{\u0421}\left({\omega}_{R}\right)-i{\Gamma}_{\lambda n}^{S}\left({\omega}_{R}\right)\right]\phi ;\\ {\stackrel{\u02dc}{\mu}}_{n}\phi \approx {\stackrel{\xaf}{\mu}}_{n}\phi ={\mu}_{n}\left[1-{\Gamma}_{\mu n}^{\u0421}\left({\omega}_{R}\right)-i{\Gamma}_{\mu n}^{S}\left({\omega}_{R}\right)\right]\phi ;\\ {\Gamma}_{\lambda n}^{c}\left({\omega}_{R}\right)={\displaystyle \underset{0}{\overset{\infty}{\int}}{R}_{\lambda n}\left(\tau \right)\mathrm{cos}{\omega}_{R}\tau \text{\hspace{0.17em}}\text{d}\tau};\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\Gamma}_{\mu n}^{c}\left({\omega}_{R}\right)={\displaystyle \underset{0}{\overset{\infty}{\int}}{R}_{\mu n}\left(\tau \right)\mathrm{cos}{\omega}_{R}\tau \text{\hspace{0.17em}}\text{d}\tau};\\ {\Gamma}_{\lambda n}^{s}\left({\omega}_{R}\right)={\displaystyle \underset{0}{\overset{\infty}{\int}}{R}_{\lambda n}\left(\tau \right)\mathrm{sin}{\omega}_{R}\tau \text{\hspace{0.17em}}\text{d}\tau};{\Gamma}_{\mu n}^{s}\left({\omega}_{R}\right)={\displaystyle \underset{0}{\overset{\infty}{\int}}{R}_{\mu n}\left(\tau \right)\mathrm{sin}{\omega}_{R}\tau \text{\hspace{0.17em}}\text{d}\tau}.\end{array}$

Physical relationships for deformable elements of zero volume

${F}_{e}=-{\stackrel{\xaf}{C}}_{e}\Delta l=-{C}_{e}\left[1-{\Gamma}_{e}^{\u0421}\left({\omega}_{R}\right)-i{\Gamma}_{e}^{S}\left({\omega}_{R}\right)\right]\Delta l,$

where $\Delta l$ ―lengthening this element; ${C}_{e}$ ―instant stiffness of the massless element. When posing the problem of proper and forced oscillations of a dissipative mechanical system, consisting of rigid and deformable bodies, the principle of possible displacements is used, according to which the sum of all active forces, acting on the system, including inertia forces, is zero:

$\delta \u0410=\delta {\u0410}_{F}+\delta {\u0410}_{u}+\delta {\u0410}_{I}=0.$ (2а)

where,

$\begin{array}{l}\delta {A}_{F}=-{\displaystyle \underset{n=1}{\overset{{s}_{2}}{\sum}}{\displaystyle \underset{{V}_{n}}{\int}{\sigma}_{ij}}}\delta {\epsilon}_{ij}\text{d}V-{\displaystyle \underset{e=1}{\overset{{s}_{1}}{\sum}}{F}_{e}}\Delta l;\\ \delta {A}_{u}=-{\displaystyle \underset{n=1}{\overset{{s}_{2}}{\sum}}{\displaystyle \underset{{V}_{n}}{\int}{\rho}_{u}\frac{{\partial}^{2}u}{\partial {t}^{2}}\delta u\text{\hspace{0.05em}}\text{d}V}}-{\displaystyle \underset{k=1}{\overset{n}{\sum}}{m}_{k}\frac{{\text{d}}^{2}u}{\text{d}{t}^{2}}\delta {u}_{k}}-{\displaystyle \underset{k=1}{\overset{n}{\sum}}{I}_{k}\frac{{\text{d}}^{2}u}{\text{d}{t}^{2}}\delta {\phi}_{k}};\\ {A}_{I}=-{\displaystyle \underset{n=1}{\overset{{s}_{2}}{\sum}}{\displaystyle \underset{{V}_{n}}{\int}{\rho}_{n}f}}\delta u\text{d}V+{\displaystyle \underset{n=1}{\overset{{s}_{1}}{\sum}}{\displaystyle \underset{{V}_{n}}{\int}f}}\delta u\text{d}V+{\displaystyle \underset{n=1}{\overset{N}{\sum}}{F}_{m}}\delta {u}_{n}+{\displaystyle \underset{k=1}{\overset{n}{\sum}}{m}_{k}\delta {\phi}_{k}};\end{array}$

here,
$\delta {\epsilon}_{ij}$ ,
$\delta \Delta l$ ―the variation of the deformations of distributed and elongated lumped elements;
${\rho}_{n}$ ―material density n-th concentrated element;
${m}_{k}$ ― weight к^{th} rigid body;
${u}_{1},{u}_{k},\delta {u}_{1},\delta {u}_{k}$ ―vectors of displacements of points of distributed elements and centers of mass of rigid bodies and their variations;
$\phi ,\delta \phi $ ―vectors of angular displacements of rigid bodies and their variations;
${V}_{n},{\sum}_{n}$ ―volume and surface n-th distributed item;
${I}_{n}$ ―central inertia tensor n-th rigid body;
${F}_{m},{M}_{k}$ ―the principal vector and principal moment of the forces, applied to к-hard body;
${s}_{1}$ ―otherwise deformability elements;
${s}_{2}$ ― otherwise the wire deformability.

If the dissipative mechanical system consists of deformable bodies, then instead of (2a) the uniform equations of the movetion are used:

${\stackrel{\u02dc}{\mu}}_{n}\frac{{\partial}^{2}{u}_{i}}{\partial {x}_{j}\partial {x}_{j}}+\left({\stackrel{\u02dc}{\lambda}}_{n}+{\stackrel{\u02dc}{\mu}}_{n}\right)\frac{{\partial}^{2}{u}_{j}}{\partial {x}_{i}\partial {x}_{j}}-{\rho}_{n}\frac{{\partial}^{2}{u}_{i}}{\partial {t}^{2}}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}x\in {V}_{n}$ (2b)

and boundary conditions

${u}_{i}=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{at}\text{\hspace{0.17em}}x\in {\sum}_{u},$

$\left[{\stackrel{\u02dc}{\lambda}}_{n}\frac{\partial {u}_{k}}{\partial {x}_{k}}{\delta}_{ij}+{\stackrel{\u02dc}{\mu}}_{n}\left(\frac{\partial {u}_{i}}{\partial {x}_{j}}+\frac{\partial {u}_{j}}{\partial {x}_{i}}\right)\right]{\nu}_{j},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{at}\text{\hspace{0.17em}}x\in {\sum}_{p},$ (3)

were i, j, k = 1, 2, 3;
$n=1,2,3,\cdots ,N$ ;
${\rho}_{n}$ ―material density n^{th} volume;
${\nu}_{j}$ ― components of the normal to the surface
${\sum}_{p}$ ;
$x=\left({x}_{1},{x}_{2},{x}_{3}\right)$ ―radius vector of the point of the body;
${u}_{i}$ ―components of the displacement vector. At the volume boundaries
${V}_{n}$ we assume continuity of the displacement vector components, normal and tangent to the stress interface. We seek the solution of problem (2) and (3) in the form

${u}_{i}\left({x}_{1},{x}_{2},{x}_{3},t\right)={U}_{i}\left({x}_{1},{x}_{2},{x}_{3}\right){\text{e}}^{-i\omega t},$ (4)

were $\omega ={\omega}_{R}+i{\omega}_{I}$ ―complex natural frequency; ${U}_{i}$ ―complex eigenmode of os- cillations. Values $\omega $ and ${U}_{i}$ we define solving the problem on eigenvalues:

$\begin{array}{l}x\in {V}_{n}:\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\stackrel{\xaf}{\mu}}_{n}\left({\omega}_{R}\right)\frac{{\partial}^{2}{U}_{i}}{\partial {x}_{j}\partial {x}_{j}}+\left({\stackrel{\xaf}{\lambda}}_{n}\left({\omega}_{R}\right)+{\stackrel{\xaf}{\mu}}_{n}\left({\omega}_{R}\right)\right)\frac{{\partial}^{2}{U}_{j}}{\partial {x}_{i}\partial {x}_{j}}+{\rho}_{n}{\omega}^{2}{U}_{i}=0;\\ x\in {\sum}_{u}:\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{U}_{i}=0;\\ x\in {\sum}_{p}:\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left[{\stackrel{\xaf}{\lambda}}_{n}\left({\omega}_{R}\right)\frac{\partial {U}_{k}}{\partial {x}_{k}}{\delta}_{ij}+{\stackrel{\xaf}{\mu}}_{n}\left({\omega}_{R}\right)\left(\frac{\partial {U}_{i}}{\partial {x}_{j}}+\frac{\partial {U}_{j}}{\partial {x}_{i}}\right)\right]\nu =0\end{array}$ (5)

We consider a construction with a finite number of degrees of freedom, consisting of a finite number of material points and absolutely rigid bodies, connected by massless viscoelastic elements. If the rheological characteristics of the elements are the same, then this system is a finite-dimensional model of a dissipative homogeneous system, under various rheological characteristics of the elements, the system under consideration simulates a dissipative inhomogeneous system with a finite number of degrees of freedom. If some elements of the mechanical system are deformed, then the designs represent a system with distributed parameters. For natural oscillations, it is necessary to determine the natural frequencies and damping coefficients of the natural oscillations. The Lagrange equations for the system under consideration have the form:

$\underset{k=1}{\overset{n}{\sum}}\left({a}_{jk}{\stackrel{\xa8}{q}}_{k}+{\stackrel{\u02dc}{c}}_{jk}{q}_{k}\right)=0},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(j=1,2,3,\cdots ,n\right)$ (6)

were ${\u0430}_{jk}$ ―symmetric positive definite matrix of generalized masses; ${q}_{k}$ ―ge- neralized coordinates; n―number of degrees of freedom; ${\stackrel{\u02dc}{\u0441}}_{jk}$ ―a symmetric matrix whose elements are the operator:

${\stackrel{\u02dc}{\u0441}}_{jk}\phi \left(t\right)={c}_{jk}\left[\phi \left(t\right)-{\displaystyle \underset{0}{\overset{t}{\int}}{R}_{jk}\left(t-\tau \right)\phi \left(\tau \right)\text{d}\tau}\right].$ (7)

here с_{jk}―elements of a symmetric positive definite matrix of instantaneous generalized rigidities; R_{jk}―elements of a nonnegatively definite matrix of generalized relaxation kernels. For a homogeneous system (all elements are made of a single viscoelastic material and described by the same relaxation nuclei, but different instantaneous stiffnesses due to different sizes), all the relaxation nuclei (7) are the same. Consequently, the matrix of generalized operator stiffnesses is a positive definite real number matrix, multiplied by the operator scalar.

Thus, in normal coordinates ${\theta}_{\u043a}$ the elastic system (6) takes the form

${{\theta}^{\u2033}}_{\u043a}+{\Omega}_{\u043a}^{2}\left[{\theta}_{\u043a}-{\displaystyle \underset{0}{\overset{t}{\int}}R\left(t-\tau \right){\theta}_{k}\left(\tau \right)\text{d}\tau}\right]=0,$ (8)

were, ${\Omega}_{\u043a}$ ―natural oscillation frequency of the elastic system. The system (8) splits into separate equations, that from the mechanical point vision, means the independence of the variation of various normal coordinates with free oscillations. For an inhomogeneous system (its deformable elements have different rheological characteristics, in particular, some of them may be elastic), the operator coefficients ${\stackrel{\u02dc}{\u0441}}_{jk}$ in (6) is the sum of two matrices-the numerical and the operator matrix:

${\stackrel{\u02dc}{\u0441}}_{jk}={\u0441}_{jk}+{\stackrel{\u02dc}{\u0412}}_{jk};\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\stackrel{\u02dc}{\u0412}}_{jk}\phi ={\u0441}_{jk}+{\displaystyle \underset{0}{\overset{t}{\int}}{R}_{jk}}\left(t-\tau \right)\phi \left(\tau \right)\text{d}\tau $ .

Three symmetric matrices ${\u0430}_{jk},{R}_{jk},{\u0441}_{jk}$ can not, in the general case, be reduced to a diagonal form by a single no degenerate transformation; therefore, for a dissipative inhomogeneous system of Lagrange’s equation in normal coordinates ${\theta}_{\u043a}$ does not decay, but takes the form:

${{\theta}^{\u2033}}_{\u043a}+{\Omega}_{\u043a}^{2}\left[{\theta}_{\u043a}-{\displaystyle \underset{j=1}{\overset{n}{\sum}}{\displaystyle \underset{0}{\overset{t}{\int}}{\Theta}_{kj}\left(t-\tau \right){\theta}_{k}\left(\tau \right)\text{d}\tau}}\right]=0,$ (9)

were ${\Theta}_{kj}$ ―symmetric positive definite matrix of generalized relaxation kernels in the normal coordinates of the elastic problem. The difference between systems (8) and (9) can be interpreted as a fact of the mutual influence of the normal coordinates of a dissipative inhomogeneous system with free oscillations. This mutual influence, which, as a rule, is neglected, can be of fundamental importance. The problems of the natural oscillations of the system, described by Equations (6), can be solved in the following way. We replace the operators ${\u0441}_{jk}$ in

(6), by complex generalized rigidities: ${\stackrel{\xaf}{\u0441}}_{jk}={c}_{jk}\left[1-{\Gamma}_{jk}^{\u0421}\left({\omega}_{R}\right)-i{\Gamma}_{jk}^{S}\left({\omega}_{R}\right)\right].$

A particular solution of system (6) is sought in the form:

${q}_{k}={U}_{k}{\text{e}}^{-i\omega t},$ (10)

were $\omega ={\omega}_{R}+i{\omega}_{I}$ ―desired complex natural frequency; ${U}_{k}$ ―components of the desired complex eigenforms. For the quantities $\omega $ and ${U}_{k}$ one can write a homogeneous linear algebraic system with a nonlinearly entering complex parameter. The characteristic equation of the eigenvalue problem has the form

${c}_{jk}\left[1-{\Gamma}_{jk}^{\u0421}\left({\omega}_{R}\right)-i{\Gamma}_{jk}^{S}\left({\omega}_{R}\right)\right]-{\omega}^{2}{a}_{jk}=0$ (11)

The roots of Equation (11) can be determined by the method of Mueller assuming the eigenfrequencies of the elastic problem as the initial approximation. The left-hand side of Equation (11) is calculated at each iteration by the method of Gauss with the separation of the principal element. Thus, to solve equation (11) we do not need to disclose the determinant on the left-hand side of.

As an example, consider a system with two degrees of freedom, consisting of two bodies of masses ${m}_{1}=1$ , ${m}_{1}=1$ , and three deformable elements with operator stiffnesses ${\stackrel{\u02dc}{\u0441}}_{1},{\stackrel{\u02dc}{\u0441}}_{2}={\stackrel{\u02dc}{\u0441}}_{3}$ (Figure 1).

The dependence of the natural frequencies and damping factors on the instantaneous stiffness ${\u0441}_{2}$ at fixed values ${m}_{1},{m}_{1},{c}_{1}$ and relaxation nucleus ${R}_{1}$ , ${R}_{2}$ . We consider two variants of the problem. For the first variant we have a homogeneous system ${R}_{1}={R}_{2}=A{\text{e}}^{-\beta t}/{t}^{1-\alpha},$ were $\u0410=0.01$ ;

$\beta =1;\alpha =0.1;{\u0441}_{1}=1;$ instant stiffness ${\u0441}_{2}={10}^{-3}$ which considered by Koltunov M.A., Mayboroda V.P., Zubcheninov V.G. [4] . The results of the calculations are shown in Figure 1. Dependence of natural frequencies and damping factors on rigidity ${\u0441}_{2}$ monotonic, and the character of the dependence is the same for frequencies and damping coefficients. In the second variant, a dissipative inho-

Figure 1. The change in complex frequencies as a function of ${\u0441}_{2}$ , for dissipative homogeneous mechanical systems.

mogeneous system is considered: the first deformable element is elastic, $\left({R}_{1}=0\right)$ the other parameters are the same as for the first option. The results of the calculations are shown in Figure 2. Dependence of natural frequencies ${\omega}_{Rk}$ from ${c}_{2}$ the same as in the case of a homogeneous system, the corresponding curves coincide with an accuracy of up to 5%. As for the coefficients of damping, their behavior varies radically: the dependence ${\omega}_{Ik}$ from ${c}_{2}$ becomes no monotonic. Of particular interest is the minimum value of the damping coefficient for a fixed value ${c}_{2}$ damping factor:

${\delta}_{\omega}=\underset{k}{\mathrm{min}}\left(-{\omega}_{Ik}\right)$

Value $\delta $ determines the damping properties of the system as a whole. In the case of a homogeneous system, the quantity ${\delta}_{\omega}$ (we call it the global damping factor) is entirely determined by the imaginary part of the smallest modulo complex eigenfrequency. In the case of an inhomogeneous system, as the global damping factor as a function of the magnitude ${c}_{2}$ the imaginary parts of both the first and second eigenfrequencies appear. “Change of roles” occurs with a characteristic value of the value ${c}_{2}$ , when the real parts of the first and second eigenfrequencies are closest. The global damping factor at the specified characteristic value ${c}_{2}$ has a pronounced maximum. This circumstance represents, in our opinion, a new mechanical effect, which can be formulated as follows: oscillations of the eigenmodes of a dissipative inhomogeneous viscoelastic system with close frequencies mutually cancel each other. Instant stiffness ${c}_{2}$ is a geometric parameter, determined by the dimensions of the element, and not by the physical properties of the material. The main feature of the observed effect is the qualitative dependence of the dissipative properties of the system on its geometric parameters.

Thus, the results obtained for the dissipative inhomogeneous viscoelastic design under consideration are completely consistent with the solutions of the problem of free damping oscillations and confirm the fact of a sharp increase in the intensity of dissipative processes when the fundamental frequencies approach

Figure 2. The dependence of complex frequencies on ${\u0441}_{2}$ .

in inhomogeneous viscoelastic systems. In this case, the role of rheology reduces both to the damping of oscillations and to the mutually increasing interaction of oscillations and different modes, which substantially increases the dissipative properties of the system as a whole. This effect is called “the effect of Troyanovsky-Safarov”.

The effect of interaction of various forms of motion of solid bodies has a principal perspective for the synthesis of dissipative heterogeneous engineering structures, construction products, damping compounds, materials and composites of various vibration protection systems and devices that are optimal for dissipative properties and material consumption.

Let us consider the natural oscillations of a system with two degrees of freedom (Figure 3). The following parameter values [2] are accepted: $A=0.048;\beta =0.05;\alpha =0.1;{\u0441}_{1}=1;M=1$ , instant stiffness ${\u0421}_{2}$ varies within $0.001{\scriptscriptstyle \frac{\u2022}{\u2022}}1.5$ , which was later contradicted by Bozorov M.B., Safarov I.I., Shokin Yu.I. [8] .

Each motion of an inhomogeneous system is a superposition of interacting oscillations of several normal coordinates and this interaction of different normal coordinates, the most intense at close natural frequencies, leads to an intensification of the dissipative processes in the system. Let the dissipative mechanical system consist of layered bodies. As an example, consider a hollow cylinder of length l, consisting of coaxial viscoelastic layers. We introduce a cylindrical coordinate system: $r,\phi ,z$ , which was later contradicted by Safarov I.I, Akhmedov M. Sh., Boltaev Z.I. [9] . Internal $r={R}_{0}$ and the outer surface is free of stress ${\sigma}_{rr}={\sigma}_{r\phi}={\sigma}_{rz}=0$ ; at the ends $z=0,z=l$ offsets ${u}_{z}=0$ and shearing stresses: ${\sigma}_{rz}={\sigma}_{\phi z}=0$ ; on the boundaries of layers: $r={r}_{1}{r}_{2}\cdots {r}_{N-1}$ offsets: ${u}_{r},{u}_{\phi},{u}_{z}$ and stresses ${\sigma}_{rr},{\sigma}_{r\phi},{\sigma}_{rz}$ are continuous. As the required task variables, we take the ${u}_{r},{u}_{\phi},{u}_{z}$ and stresses ${\sigma}_{rr},{\sigma}_{r\phi},{\sigma}_{rz}$ . Equations for these variables are obtained by eliminating from the Cauchy relations, Hooke’s law and the equations of motion of strains and stresses ${\sigma}_{\phi \phi},{\sigma}_{zz},{\sigma}_{\phi z}$ ( ${r}_{n-1}<r<{r}_{n},n=1,\cdots ,N$ ), which was later contradicted by Safarov I.I, Akhmedov M. Sh., Boltaev Z.I. [10] :

$\begin{array}{l}\frac{\partial {u}_{r}}{\partial r}=-\frac{{\stackrel{\xaf}{\nu}}_{n}}{1-{\stackrel{\xaf}{\nu}}_{n}}\left(\frac{{u}_{r}}{r}+\frac{1}{r}\frac{\partial {u}_{\phi}}{\partial \phi}+\frac{\partial {u}_{z}}{\partial z}\right)+\frac{\left(1+{\stackrel{\xaf}{\nu}}_{n}\right)\left(1-2{\stackrel{\xaf}{\nu}}_{n}\right)}{{\stackrel{\xaf}{E}}_{n}\left(1-{\stackrel{\xaf}{\nu}}_{n}\right)}{\sigma}_{rr};\\ \frac{\partial {u}_{\phi}}{\partial r}=-\frac{1}{r}\frac{\partial {u}_{r}}{\partial \phi}+\frac{{u}_{\phi}}{r}+\frac{2\left(1+{\stackrel{\xaf}{\nu}}_{n}\right)}{{\stackrel{\xaf}{E}}_{n}}{\sigma}_{r\phi};\\ \frac{\partial {u}_{z}}{\partial r}=-\frac{\partial {u}_{r}}{\partial z}+\frac{2\left(1+{\stackrel{\xaf}{\nu}}_{n}\right)}{{\stackrel{\xaf}{E}}_{n}}{\sigma}_{rz};\end{array}$ (12)

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

where

(a)(b)

Figure 3. The dependence of complex frequencies on ${c}_{2}$ ; (a) Dissipatively homogeneous mechanical system; (b) Dissipatively homogeneous mechanical system.

$\begin{array}{l}{\sigma}_{\phi \phi}=\frac{{\stackrel{\xaf}{E}}_{n}}{2\left(1+{\stackrel{\xaf}{\nu}}_{n}\right)}\left(\frac{{u}_{r}}{r}+\frac{1}{r}\frac{\partial {u}_{\phi}}{\partial \phi}+{\stackrel{\xaf}{\nu}}_{n}\frac{\partial {u}_{z}}{\partial z}\right)+\frac{{\stackrel{\xaf}{\nu}}_{n}}{1-{\stackrel{\xaf}{\nu}}_{n}}{\sigma}_{rr};\\ {\sigma}_{\phi z}=\frac{{\stackrel{\xaf}{E}}_{n}}{2\left(1+{\stackrel{\u02dc}{\nu}}_{n}\right)}\left(\frac{\partial {u}_{\phi}}{\partial z}+\frac{1}{r}\frac{\partial {u}_{z}}{\partial \phi}\right);\end{array}$

$\begin{array}{l}{\sigma}_{zz}=\frac{{\stackrel{\xaf}{E}}_{n}}{1-{\stackrel{\xaf}{\nu}}_{n}^{2}}\left({\stackrel{\xaf}{\nu}}_{n}\frac{{u}_{r}}{r}+\frac{{\stackrel{\xaf}{\nu}}_{n}}{r}\frac{\partial {u}_{\phi}}{\partial \phi}+\frac{\partial {u}_{z}}{\partial z}\right)+\frac{{\stackrel{\xaf}{\nu}}_{n}}{1-{\stackrel{\xaf}{\nu}}_{n}}{\sigma}_{rr};\\ {\stackrel{\u02dc}{E}}_{n}=\frac{{\stackrel{\xaf}{\mu}}_{n}\left(3{\stackrel{\xaf}{\lambda}}_{n}+2{\stackrel{\xaf}{\mu}}_{n}\right)}{{\stackrel{\xaf}{\mu}}_{n}+{\stackrel{\xaf}{\lambda}}_{n}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\stackrel{\xaf}{\nu}}_{n}=\frac{{\stackrel{\xaf}{\lambda}}_{n}}{2\left({\stackrel{\xaf}{\mu}}_{n}+{\stackrel{\xaf}{\lambda}}_{n}\right)}.\end{array}$

The boundary conditions of the problem have the form:

$\begin{array}{l}r={R}_{0},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}r={r}_{n}\text{}:\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\sigma}_{rr}={\sigma}_{r\phi}={\sigma}_{rz}=0;\\ z=0,\text{\hspace{0.17em}}z=l:\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\sigma}_{rz}={\sigma}_{\phi z}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}{u}_{z}=0\end{array}$ (13)

The problem (12)-(13) does not degenerate when ${\stackrel{\u02dc}{\nu}}_{n}=0.5$ (incompressible medium).

We seek the solution of the problem in the form:

$\left(\begin{array}{c}{u}_{r}\\ {u}_{\phi}\\ {u}_{z}\end{array}\right)=\left(\begin{array}{c}{U}_{r}\\ i{U}_{\phi}\\ i{U}_{z}\end{array}\right){\text{e}}^{-i\left(\frac{k\text{\pi}}{l}z+\omega t+m\phi \right)}$ , (14a)

$\left(\begin{array}{c}{\sigma}_{rr}\\ {\sigma}_{r\phi}\\ {\sigma}_{rz}\end{array}\right)=\left(\begin{array}{c}{\tau}_{rr}\\ i{\tau}_{r\phi}\\ i{\tau}_{rz}\end{array}\right){\text{e}}^{-i\left(\frac{k\text{\pi}}{l}z+\omega t+m\phi \right)}$ , (14b)

where $m,k$ ―positive integers; ${U}_{r},{U}_{\phi},{U}_{z},{\tau}_{rr},{\tau}_{r\phi},{\tau}_{rz}$ ―complex amplitudes de- pending on the radius r. The problem reduces to finding the eigenvalues for a system of six ordinary differential equations with respect to complex amplitudes ( ${r}_{n-1}<r<{r}_{n},\text{\hspace{0.17em}}n=1,\cdots ,N$ ):

#Math_132# (15a)

$\begin{array}{l}\frac{\text{d}{\tau}_{rr}}{\text{d}r}=i\frac{m}{r}{\tau}_{r\phi}+i\frac{k\text{\pi}}{l}{\tau}_{rz}-\frac{{\tau}_{rr}-{\tau}_{r\phi}}{r}-{\rho}_{n}{\omega}^{2}{U}_{r};\\ \frac{\text{d}{\tau}_{r\phi}}{\text{d}r}=i\frac{m}{r}{\tau}_{\phi \phi}+i\frac{k\text{\pi}}{l}{\tau}_{\phi z}-\frac{2{\tau}_{r\phi}}{r}-{\rho}_{n}{\omega}^{2}{U}_{\phi};\\ \frac{\text{d}{\tau}_{rz}}{\text{d}r}=i\frac{m}{r}{\tau}_{r\phi}+i\frac{k\text{\pi}}{l}{\tau}_{rz}-\frac{{\tau}_{rz}}{r}-{\rho}_{n}{\omega}^{2}{U}_{z}.\end{array}$ (15b)

here

$\begin{array}{l}{\tau}_{\phi \phi}=\frac{{\stackrel{\xaf}{E}}_{n}}{2\left(1+{\stackrel{\xaf}{\nu}}_{n}\right)}\left(\frac{{U}_{r}}{r}-i\frac{m}{r}{U}_{\phi}-i{\stackrel{\xaf}{\nu}}_{n}\frac{k\text{\pi}}{l}{U}_{r}\right)+\frac{{\stackrel{\xaf}{\nu}}_{n}}{1-{\stackrel{\xaf}{\nu}}_{n}}{\tau}_{rr};\\ {\tau}_{\phi z}=-\frac{i{\stackrel{\xaf}{E}}_{n}}{2\left(1+{\stackrel{\xaf}{\nu}}_{n}\right)}\left(\frac{k\text{\pi}}{l}{U}_{\phi}+\frac{m}{r}{U}_{z}\right);\\ {\tau}_{zz}=\frac{{\stackrel{\xaf}{E}}_{n}}{1-{\stackrel{\xaf}{\nu}}_{n}^{2}}\left({\stackrel{\xaf}{\nu}}_{n}\frac{{U}_{r}}{r}-i\frac{{\stackrel{\xaf}{\nu}}_{n}m}{r}{U}_{\phi}-i\frac{k\text{\pi}}{l}{U}_{z}\right)+\frac{{\stackrel{\xaf}{\nu}}_{n}}{1-{\stackrel{\xaf}{\nu}}_{n}}{\tau}_{rr};\end{array}$

when

$r={R}_{0},\text{\hspace{0.17em}}r={r}_{n}:\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\tau}_{rr}={\tau}_{r\phi}={\tau}_{rz}=0;\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}z=0\text{\hspace{0.17em}},z=l:\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\tau}_{rz}={\tau}_{\phi z}=0,\text{\hspace{0.17em}}{U}_{z}=0.$

The characteristic equation of the problem is constructed by the method of orthogonal sweep, the roots of the characteristic equation are determined by the Muller method. Figure 4 shows the results of calculations for the case of plane oscillations ( $m=2,\text{\hspace{0.17em}}k=0$ ) cylinder, consisting of an outer thin elastic layer ( ${r}_{1}=1.5;{r}_{2}=1.7$ and an internal incompressible viscoelastic layer). Inner radius ${R}_{0}$ , which varied from 0.6 to 1.3. The following parameter values are accepted: ${E}_{1}={10}^{-4}{\nu}_{1}=0/5$ , ${R}_{\mu 1}=A{\text{e}}^{-\beta t}/{t}^{1-\alpha}$ ,#Math_141#; $\beta =1$ ; $\alpha =0.1$ ; ${E}_{2}=1$ ; ${\nu}_{2}=0.3$ ; ${R}_{2}=0$ . Figure 4 shows the dependences of the real and imaginary parts of the first two natural frequencies on the inner radius ${R}_{0}$ . The above dependences confirm the effect found in the analysis of the natural oscillations of a system

Figure 4. The change in complex frequencies as a function of ${c}_{2}$ , for dissipative inhomogeneous mechanical systems.

with a finite number of degrees of freedom.

3. Forced Oscillations of Dissipative Mechanical Systems

For a system with a finite number of degrees of freedom, the variational problem (2a) reduces to a system of linear Lagrange equations of the second kind with complex generalized rigidity:

$\underset{k=1}{\overset{N}{\sum}}\left({a}_{jk}{\stackrel{\xa8}{q}}_{k}+{\stackrel{\u02dc}{\u0441}}_{jk}{q}_{k}\right)={f}_{j}{\text{e}}^{-i\lambda t}},\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(j=1,2,3,\cdots ,N\right)$ (16)

where

${\stackrel{\u02dc}{c}}_{jk}\phi \left(t\right)={c}_{jk}\left[\phi \left(t\right)-{\displaystyle \underset{-\infty}{\overset{t}{\int}}{R}_{jk}\left(t-\tau \right)\phi \left(\tau \right)\text{d}\tau}\right]$ ,

${a}_{ij}$ ―components of the real symmetric matrix of generalized masses; ${q}_{k}$ ―com- plex generalized coordinates; ${f}_{j}$ ―complex amplitudes of generalized forces; $\lambda $ ―actual frequency of external influences. The solution of the problem of forced oscillations of system (16) is sought in the form:

${q}_{j}={A}_{j}{\text{e}}^{-i\lambda t}$ , (17)

where ${A}_{j}$ desired complex amplitudes. The problem of forced oscillations, reduces to solving an inhomogeneous algebraic system:

$\underset{j=1}{\overset{N}{\sum}}\left({c}_{lk}\left(\lambda \right)-{\lambda}^{2}{a}_{jk}\right){A}_{k}}={f}_{j$ , (18)

the solution, which is carried out in the Gauss method. In the presence of a medium of deformable elements of a system of at least one with a distributed mass, the number of its degrees of freedom becomes infinite (countable). In this paper we propose two approaches to solving an infinite-dimensional version of the variational problem (2a). The first approach reduces to constructing an infinite system of Lagrange equations of the second kind, which after truncation turns into a finite system of the form (16). As the generalized coordinates, the components of the displacements of the centers of mass of rigid bodies, and their small angles of rotation with respect to the coordinate axes (or principal directions of inertia) are taken, as well as, the components of the expansion of the displacements of the massive deformable elements from their elastic self-oscillations (here we mean the eigenmodes of the vibrations of an elastic isolated deformable element). In this case the variational problem (2a) reduces to a system of partial differential equations, which has been widely studied [10] :

$\underset{k=1}{\overset{{s}_{2}}{\sum}}{L}_{jk}}{W}_{j}\pm {\rho}_{j}\frac{{\partial}^{2}{W}_{j}}{\partial {t}^{2}}=\rho {f}_{j}{\text{e}}^{-i\lambda t$ , (19)

wherе
${W}_{j}$ ―components of the displacement vector of points j^{th} massive of the deformable element;
${\rho}_{j}$ ―its linear density, depending on the coordinates;
${f}_{j}$ ― components of complex amplitudes of mass forces, including portable inertia forces for a given vibration of a rigid base;
${L}_{jk}$ ―linear combinations of differentiation operators with respect to the spatial coordinate of the k^{th} rod with complex coefficients. These coefficients in the case of natural oscillations depend on the real parts of the desired complex eigenfrequencies, and for the steady- state ones, on the actual frequency of the external action. The system (19) does not include masses and moments of inertia of absolutely rigid bodies and complex stiffnesses of massless deformable elements of zero volume. These parameters in the approach are taken into account, both in the boundary conditions and in the internal (depending on the location of the attached masses) points of the j^{th} rod. As shearing forces and bending moments, the elastic forces of concentrated deformed elements, forces and moments of inertia forces (relative portable ones) applied to the absolute rigid body of the system will enter into the specified boundary conditions. The boundary conditions, for the system (19), in the general case are not written out in view of the extreme variety. Let us determine the dependence of the resonance amplitudes
${A}_{jk}$ (j^{th} number of the generalized coordinate, k is the number of the resonant frequency ) (18) from the system parameters. The algorithm for constructing this dependence includes the construction of amplitude-frequency characteristics, for each generalized coordinate
${q}_{j}$ and finding the maximums
${A}_{jk}$ for each
${\lambda}_{j}$ . To find the maxima and minima on the resonance curve, the Mueller method is applied, which has been widely studied [10] . The amplitude-frequency response, this curve, is described by equation
${A}_{j}\left(\lambda \right)=\left|{q}_{j}\left(\lambda \right)\right|$ . Its maxima and minima satisfy equation

$\frac{\partial {A}_{j}\left(\lambda \right)}{\partial \lambda}=0$ (20)

It is necessary to have the roots of this equation (by the Muller method).

To do this, we must calculate the left-hand sides of (20) for given values $\lambda ={\lambda}_{1},{\lambda}_{2},\cdots ,{\lambda}_{n}$ . Having chosen $\Delta \lambda ={10}^{-4}\text{\hspace{0.17em}}\text{-}\text{\hspace{0.17em}}{10}^{-6}$ , we define the derivatives by the formula:

$\frac{\partial {A}_{j}\left(\lambda \right)}{\partial \lambda}=\frac{{A}_{j}\left({\lambda}_{j}+\Delta \lambda \right)-{A}_{j}\left({\lambda}_{j}-\Delta \lambda \right)}{2\Delta \lambda}=0$ .

The first root is known #Math_173# , so you need to start looking for the second one. As an example, we consider a system with two degrees of freedom (Figure 3) , which has been widely studied [11] . The values are accepted, the other parameters coincide with those adopted in the second task. Two variants of mechanical systems are considered. The results of the calculation of the first variant are shown in Figure 5(a). Addiction
${A}_{jk}$ (j = 1, 2; k = 1, 2) of the parameter
${\u0441}_{2}$ monotonous. In the second variant, the first deformable element is elastic:
${R}_{1}=0$ , the other parameters are the same as those used above. The results of the calculation are shown in Figure 5(b), according to which when the natural frequencies approach the corresponding amplitudes А_{11} and А_{12} becomes equal. Addiction
${A}_{jk}$ from
${\u0441}_{2}$ no monotonic. The damping properties of the system as a whole for forced oscillations are determined by the maximum resonance amplitude (we call it the global resonant amplitude). The intensity of dissipative processes in the system is higher the lower the global resonant amplitude (GRA):
${\delta}_{\u0410}=\underset{\lambda}{\mathrm{max}}\left({A}_{jk}\right)$ .

In a dissipative homogeneous system, the role of the global resonant amplitude is fulfilled for all values of the parameter by the first resonance amplitude. In a dissipative inhomogeneous system, both the first and second resonant amplitudes act as the GRA in terms of the magnitude of the parameter. “Change of roles” also occurs in the case of a global damping factor ${\delta}_{\u0410}$ with a characteristic value of the parameter at which the real parts of the natural frequencies are closest.

At this value of the parameter, the global resonant amplitude (GRA) is minimal and, consequently, the dissipative processes in the system proceed most intensively, and the global damping coefficient has a pronounced maximum. This effect is a continuation of the effect of Troyanovsky-Safarov. To clarify the physical nature of the detected effect, we write the equations of motion of a system with n degrees of freedom. In the case of a homogeneous system, all the relaxation nuclei ${R}_{ij}$ is the same: ${R}_{ij}=R$ , since the matrix of generalized complex rigidity ${\u0441}_{ij}$ is a positive definite real matrix ${\u0441}_{ij}$ , multiplied by a complex scalar:

${\stackrel{\xaf}{\u0441}}_{jk}={c}_{jk}\left[1-{\Gamma}_{jk}^{\u0421}\left({\omega}_{R}\right)-i{\Gamma}_{jk}^{S}\left({\omega}_{R}\right)\right].$

In the normal coordinates of the viscoelastic problem, the system (16) takes the form:

${\stackrel{\xa8}{\theta}}_{n}+{\Omega}_{n}^{2}{\theta}_{n}\left(1-{\Gamma}^{c}-i{\Gamma}^{s}\right)={\Psi}_{n}$ (21)

(a)(b)

Figure 5. Dependence of the resonance amplitude on с_{2}; (a) Dissipatively homogeneous mechanical system; (b) Dissipatively homogeneous mechanical system.

where
$\Omega $ ―complex own frequency of elastic system;
${\Psi}_{n}$ ―generalized force corresponding to n^{th} normal coordinate. The system (21) is divided into n separate equations. This means that, the motion of a mechanical viscoelastic system is a superposition of independent normal vibrations decay, and the forced ones have a finite resonant amplitude. The main property of conservative systems― the possibility of excitation of the vibration of one normal coordinate without excitation of the others―is completely preserved in the case of a homogeneous viscoelastic system. Therefore, in the case of an inhomogeneous system, the Lagrange equation in normal cords of an elastic system has the form

${\stackrel{\xa8}{\theta}}_{n}+{\Omega}_{n}^{2}{\theta}_{n}-{\Omega}_{n}^{2}{\displaystyle \underset{j=1}{\overset{N}{\sum}}\left({\theta}_{nj}^{c}+{\theta}_{nj}^{s}\right){\theta}_{j}={\Psi}_{k}}$ (22)

where ${\theta}_{nj}^{c},\text{\hspace{0.17em}}{\theta}_{nj}^{s}$ ―non-negative definition of a real matrix. The system (22) consists of related equations. This mechanical coupling means that it is impossible to excite the oscillation of an individual normal coordinate. This effect has been widely studied [12] [13] [14] . Each motion of an inhomogeneous system is a superposition of interacting oscillations of several normal coordinates, and this interaction of different normal coordinates, the most intense at close natural frequencies, leads to an intensification of the dissipative processes in the system.

References

[1] Safarov, I.I. (1992) Vibrations of a Wave in Dissipatively Inhomogeneous Media and Structures. Tashkent: Fan, 252 p.

[2] Bozorov, M.B., Safarov, I.I. and Shokin, Yu.I. (1996) Numerical Simulation of Oscillations of Dissipatively Homogeneous and Inhomogeneous Mechanical Systems. Sibirian Branch of the Russian Akademy of Science, Novosibirsk, 189 p.

[3] Safarov, I.I., Teshaev, M.Kh. and Majidov, M. (2012) Damping of Oscillations of Dissipative-Inhomogeneous Mechanical Systems. Foundations, Concepts, Methods. LAP, Lambert Academic Publishing (Germany), 8 p.

[4] Koltunov, M.A., Mayboroda, V.P. and Zubcheninov, V.G. (1983) Strength Calculations of Products Made of Polymer Materials. Mashinostroenie, Moscow, 239 p.

[5] Koltunov, M.A. (1976) Creep and Relaxation. Higher Shkola, Moscow, 277 p.

[6] Leibenzon, L.S. (1943) Variational Methods for Solving Problems in the Theory of Elasticity. Gostekhizdat, Moscow, 286 p.

[7] Mayboroda, V.P., Troyanovsky, I.E. and Safarov, I.I. (1983) Free and Forced Oscillations of Systems of Solids on Inhomogeneous Viscoelastic Shock Absorbers. Journal of the USSR Academy of Sciences, 3, 71-77.

[8] Safarov, I.I., Teshaev, M.Kh. and Madjidov, M. (2014) Natural Oscillations of Viscoelastic Lamellar Mechanical Systems with Point Communications. Applied Mathematics, 5, 3018-3025.

https://doi.org/10.4236/am.2014.519289

[9] Safarov, I.I., Akhmedov, M.Sh. and Boltaev, Z.I. (2014) Loose Waves in Viscoelastic Cylindrical Wave Guide with Radial Crack. Applied Mathematics, 6, 214-225.

https://doi.org/10.4236/am.2014.521329

[10] Safarov, I.I., Akhmedov, M.Sh. and Boltaev, Z.I. (2015) Setting the Linear Oscillations of Structural Heterogeneity. Viscoelastic Lamellar Systems with Point Relations. Applied Mathematics, 6, 225-234.

https://doi.org/10.4236/am.2015.62022

[11] Safarov, I.I., Akhmedov, M.Sh. and Boltaev, Z.I. (2015) Natural Oscillations of Cylindrical Bodies with External Friction on the Boundary. Applied Mathematics, 6, 629-645.

https://doi.org/10.4236/am.2015.63057

[12] Safarov, I.I., Boltaev, Z.I. and Umarov, A.O. (2013) Forced Oscillations of Cylindrical Bodies with External Friction at the Boundary. Interuniversity Collection of Scientific Papers “Problems of Mechanics and Management. Nonlinear Dynamic Systems. Permian, 45, 114-122.

[13] Safarov, I.I., Teshaev, M.K.H. and Boltaev, Z.I. (2016) Waves in a Cylindrical Shell with a Viscous Liquid. Vestnik of Perm University. Mathematics. Mechanics. Computer Science, Perm, 3, 82-93.

[14] Safarov, I.I., Boltaev, Z.I. and Axmedov, M.Sh. (2016) Ford Hormonal Vibrations in a Visocoelastic Layered Bodies Lying Deformable Nalf-Plane. International Journal of Applied Research, 2, 138-147.