^{1}

^{2}

^{3}

^{*}

A vibrational statement, method and algorithm to assess the damping capacity of structurally inhomogeneous viscoelastic mechanical systems, consisting of a package of rectangular plates and shells with point relations and concentrated masses at different rheological properties of deformable elements, are proposed in the paper. To describe rheological properties of the material, the linear hereditary Boltzmann-Volter theory was used. To assess the damping capacity of the system, the problem in question in each case was reduced to solving the proper problems of algebraic equations with complex parameters solved by the Muller method. The accuracy of the methods was demonstrated by comparing the calculated results with known published data and a numerical experiment. Complex natural frequency of the system was used to assess the damping capacity of inhomogeneous viscoelastic systems. Various eigenvalue problems have been solved for structurally inhomogeneous mechanical systems consisting of a package of plate and shell systems with concentrated masses and shock absorbers. A number of new mechanical effects have been discovered, related to the manifestation of the damping capacity of mechanical systems under consideration.

Currently, many technical structures are widely used shell and plate structures. Thin-walled pipes (shells), plates and panels in real conditions, as a rule, interact with other structures and bodies, rely on point rigid and (or) elastic, articulated and (or) pinched supports, have point-attached masses. Boundary conditions often differ from homogeneous. If under static loads these features cause a local redistribution of stress and strain fields, then in dynamic problems point connections significantly affect all the dynamic characteristics of a plate or shell—the spectrum of natural frequencies, the forms of natural vibrations, resonant frequencies and amplitudes, etc., for example:

Leissa [

Boay [

Avalos et al. [

Xiang et al. [

Laura and Grossi [

Wu and Luo [

Nicholson and Bergman [

The design of submarine tunnels is inevitably characterized by the presence of attached masses due to structural, strength or operational necessity [

Rod systems have a high degree of vibration isolation. E. Skudrzyk [

To generalize the methods of the above works for a finite number of mechanical systems, consisting of dissipative plates and shells is an impossible task.

The task is complicated if the plate or shell has hereditary properties, which has recently become increasingly relevant, especially in vibration protection systems [

Of great practical interest are the dynamic problems of dissipative mechanical systems, the elements of which have different dissipative properties. In the future, such mechanical systems are called structurally inhomogeneous viscoelastic systems, the elements of which have different rheological properties.

An analysis of the work on the natural vibrations of elastic plate mechanical systems shows that practically no methods have been developed for the dynamic calculation of plate mechanical systems whose elements have different rheological properties (structurally inhomogeneous mechanical system).

Structurally heterogeneous mechanical systems with a finite number of degrees of freedom are addressed by I. I. Safarov and et al. [

In the works of I. I. Safarov and et al. [

In the scientific literature there are a large number of publications on the study of the dynamic characteristics of shells bearing attached masses [

For a package of plates (or shells) with point bonds and concentrated masses, the traditional approach is possible when the directions and coordinates of the active and passive loads must be known.

This is a difficult task. Therefore, the traditional approach is applied for a single plate with concentrated masses. For two or more plates with concentrated masses and point bonds, a variational approach is applied, which allows one to obtain a system of integro-differential equations satisfying boundary conditions and allows one to analyze the dissipative properties of the mechanical system as a whole.

To date, the question of assessing the level of oscillations and dissipative properties of plates and plate systems of new structural elements of various apparatuses is still insufficiently studied and requires extensive research. Therefore, this problem is relevant, requiring its gradual solution.

Consider a mechanical system consisting of N isotropic viscoelastic bodies (a package of rectangular plates or cylindrical shells) occupying volume V n and restricted by surfaces Ω n ( n = 1 , 2 , ⋯ , N ) . It is assumed that one linear size of each body is much smaller than the other two (the class of plates and shells). For each n, homogeneous boundary conditions Ω n b o are set on parts of the surface of the n-th body; kinematic and dynamic relations are imposed on the remaining free surface Ω n f r = Ω n / Ω n b o at a finite number of points: point-like rigid, elastic and (or) viscoelastic hinge-type supports (rigid supports can be the fixed ones), rigid elastic and (or) viscoelastic shock absorbers connecting the bodies (at N > 1 ), concentrated masses M q n ( q = 1 , 2 , ⋯ , Q ) . The arrangement of relations and masses on surfaces Ω n f r is arbitrary.

A special case of such a structurally inhomogeneous viscoelastic system is a system with elastic and viscoelastic elements.

Then, the system under consideration is transformed into a non-conservative one with complex Eigen frequencies, i.e. ω = ω R + i ω I . The real part of the complex natural frequency ω R means the natural frequency of the system, the fictitious part ω I determines the damping rate of vibrations and has the meaning of the coefficient of damping. In the general case, the dissipative (damping) properties of the elements of such systems are different.

Further, in the problems under consideration, it is required to determine the natural vibrations frequencies of the inhomogeneous viscoelastic system, and to estimate its damping capacity.

In mathematical statement, a viscoelastic problem is as follows: let all points of the n-th body obey the harmonic law of vibrations, i.e.

U n j ( x ¯ n , t ) = U n j 0 ( x ¯ n ) e − i ω t , n = 1 , ⋯ , N ; j = 1 , ⋯ , J , (1)

where U n j 0 ( x ¯ n ) is the j-th component of the displacement vector of the n-th body, J is the number of components of the displacements vector, x ¯ n = ( x ¯ 1 n , x ¯ 2 n , x ¯ 3 n ) is the radius vector of the point of the n-th body, ω = ω R + i ω I is the sought-for complex frequency of the system, ω R is the natural frequency and ω I is the coefficient of damping ( ω I < 0 ). Since each component of the displacements vector already has an index n, the latter is not used in what follows to designate the components of the radius vector. For rectangular plates J = 1

U n 1 0 ( x 1 , x 2 ) = W n 0 ( x , y ) ,

and for the shells of revolution J = 3

U n 1 0 ( x 1 , x 2 ) = U n 0 ( x , y ) , U n 2 0 ( x 1 , x 2 ) = V n 0 ( x , y ) , U n 3 0 ( x 1 , x 2 ) = W n 0 ( x , y ) ,

where x, y are the coordinates.

Based on the principle of possible displacements, we equate the sum of the virtual work of all active forces, including inertial ones, on possible displacements δ U n j 0 ( x ¯ , t ) , to zero

δ A σ + δ A a + δ A m = 0 , (2)

where δ A σ , δ A a , δ A m is a virtual work of the internal forces of the bodies of springs, and the forces of inertia, with account for concentrated masses. This work can be represented by the following relations:

δ A σ = − ∑ n = 1 N σ m k n δ ε m k n d V ,

δ A a = − ∑ n = 1 N − 1 ∑ l = 1 L n σ l n δ ε l n − ∑ n = 1 N ∑ l ′ = 1 L ′ n σ l ′ n δ ε l ′ n (3)

δ A m = − ∑ n = 1 N ρ n ∫ V n ( ∑ j = 1 J U ¨ n j ( x ¯ , t ) δ U n j ) d V − ∑ n = 1 N ∑ q = 1 Q n M q n ∑ j = 1 J U ¨ n j ( x ¯ n q , t ) δ U n j

where x ¯ n q = ( x ¯ n 1 q , x ¯ n 2 q , x ¯ n 3 q ) are the coordinates, L n is the number of deformable elements (springs or shock absorbers) between the n-th and (n + 1)-th bodies, Q n is the number of concentrated masses on the n-th body, L ′ n is the number of elastic (viscoelastic) supports on the n-th body, σ m k n , ε m k n , σ l n , ε l n , σ l ′ n , ε l ′ n are the components of stress and strain tensors, respectively, of the n-th body, l-th spring (deformable element or shock absorber), and l ′ -th elastic (viscoelastic) support.

To describe viscoelastic properties of a body material with the linear hereditary Boltzmann-Volter theory, physical relations for the n-th viscoelastic element of the system are defined as [

σ m k n ( t ) = λ ˜ n Θ n ( t ) δ m k + 2 μ ˜ n ε m k n ( t ) (4)

where λ ˜ n , μ ˜ n are the Volter integral operators, replaced below by one operator. The Poisson’s ratio ν n in the proposed statement of the problem is assumed constant. This means that for a structurally homogeneous viscoelastic system, the modes of natural vibrations are equal to the eigenvectors of corresponding elastic problem [

σ m k n ( t ) = E ˜ n 1 + ν n [ ν n 1 − 2 ν n Θ n ( t ) δ m k + ε m k n ( t ) ] (5)

where Е ˜ n is the Volter operator of the form [

E ˜ n φ ( t ) = E 0 n [ φ ( t ) − ∫ 0 t R E n ( t − τ ) φ ( τ ) d τ ] (6)

here E 0 n is the instantaneous modulus of elasticity, and R E n is the kernel of relaxation.

Given (1), the time function in Equality (6) is φ ( t ) = exp ( − i ω t ) at a slowly varying amplitude. Assuming the smallness of the integral term ∫ 0 ∞ R ( τ ) d τ , using the freeze-etch method [

E ˜ n φ ( t ) ≅ E 0 j [ 1 − Γ n c ( ω R ) − i Γ n s ( ω R ) ] φ ( t ) ≡ E ¯ n φ ( t ) , (7)

where

Γ n c ( ω R ) = ∫ 0 ∞ R n ( τ ) cos ( ω R τ ) d τ , Γ n s ( ω R ) = ∫ 0 ∞ R n ( τ ) sin ( ω R τ ) d τ .

This allows eliminating the integral terms and, ultimately, time from the vibrational equation. In symbolic form, it can be represented as

δ G ( U n j 0 ( x ¯ ) , ω 2 ) = 0 . (8)

Write out a specific representation of the functional G, for example, for a package of rectangular plates with point relations:

G [ W n 0 ( x , y ) , ω 2 ] = − 1 2 ∑ n = 1 N D ¯ n ∫ 0 a n ∫ 0 b n [ ( ∂ 2 W n 0 ∂ x 2 + ∂ 2 W n 0 ∂ y 2 ) 2 − 2 ( 1 − ν n ) ( ∂ 2 W n 0 ∂ x 2 ∂ 2 W n 0 ∂ y 2 − ( ∂ 2 W n 0 ∂ x ∂ y ) 2 ) ] d x d y − 1 2 ∑ n = 1 N − 1 ∑ l = 1 L n D ¯ n [ W n 0 ( x n l , y n l ) − W n + 1 0 ( x n l , y n l ) ] 2 − 1 2 ∑ n = 1 N ∑ l ′ = 1 L ′ n C l ′ n ( W n 0 ) 2 ( x n l ′ , y n l ′ ) + ω 2 2 ∑ n = 1 N ρ n h n ∫ 0 a n ∫ 0 b n ( W n 0 ) 2 d x d y + ω 2 2 ∑ n = 1 N ∑ q = 1 Q n M q n ( W n 0 ) 2 ( x n q , y n q ) ,

where h n , a n , b n are the thickness and linear dimensions of the n-th plate, x n q , y n q are the coordinates of the n-th concentrated mass, x n l , y n l are the coordinates of the l-th spring (shock absorber), x n l ′ , y n l ′ are the coordinates of the l ′ -th elastic (viscoelastic) support.

If the n-th plate, l-th spring, and l ′ -th support are viscoelastic, then D ¯ n , C ¯ l n , С ¯ l ′ n are represented by the following formulas:

D ¯ n = D n f n ( ω R ) , C ¯ l n = C l n f l n ( ω R ) , С ¯ l ′ n = С l ′ n f l ′ n ( ω R ) ,

where f ( ω R ) = 1 − Γ c ( ω R ) − i Γ s ( ω R ) is a complex function whose numerical coefficients depend on the parameters of the relaxation kernel of corresponding viscoelastic elements, D n = E n h n 3 / ( 12 ( 1 − ν n 2 ) ) , C l n , C l ′ n is the generalized instantaneous stiffness of the n-th plate, l-th shock absorber and l ′ -th support, respectively. In the elastic case, D ¯ n = D n , C ¯ l n = C l n , С ¯ l ′ n = С l ′ n where D n , С l n , С l ′ n is the generalized stiffness of the n-th plate, l-th shock absorber and l ′ -th support, respectively. A similar functional can be written for the system of the shells of revolution. The components of the displacements vector U n j 0 ( x ) are the sought for functions of vibrational Equation (8) and must satisfy the boundary conditions on the surface Ω n b o ,

L n U n j 0 ( x ¯ ) = 0 , x ¯ ∈ Ω n b o . (9)

It remains to impose on the system rigid point relations that do not perform work under vibrations. The conditions of rigid hinge support of the n-th body in point supports S_{n} are written as

U n j 0 ( x ¯ n s ) = 0 ( s = 1 , ⋯ , S n ; j = 1 , ⋯ , J ) , (10)

where x ¯ n s are the coordinates of the s-the support of the n-the body.

If the part of supports is fixed, then the following conditions are added.

∂ U n j 0 ( x ¯ n s ) ∂ α n s = 0 , ( s = 1 , ⋯ , S n α ; j = 1 , ⋯ , J ) (11)

where α n s is the direction of the unit vector along which the rigid fixing of the body is done at the point x n s .

In the program implementing the algorithm, condition (11) is considered only for the shells of revolution. The presence of rigid posts between the n-th and (n+1)-the bodies at N ≥ 2 is taken into account by the relationships

U n j 0 ( x ¯ n r ) − U n + 1 , j 0 ( x ¯ n r ) = 0 ( r = 1 , ⋯ , R n ; j = 1 , ⋯ , J ) , (12)

where x ¯ r is the coordinate of the r-the post, R n is the number of posts between the n-the and (n + 1)-the bodies. At N = 1 conditions (12) are absent.

Thus, restrictions of the type (10)-(12) are imposed on the displacement vector. The imposition on the system of point relations is considered using the Lagrange multipliers method. Then the vibrational Equation (8) is rewritten as

δ { G ( U n j 0 ( x ¯ ) , ω 2 ) + ∑ n = 1 N ∑ s = 1 S n ∑ j = 1 J λ n j s U n j 0 ( x ¯ n s ) + ∑ n = 1 N ∑ s = 1 S n α ∑ j = 1 J k n j s ∂ U n j 0 ( x ¯ n s ) ∂ α n s + ∑ n = 1 N − 1 ∑ r = 1 R n ∑ j = 1 J μ n j r [ U n j 0 ( x ¯ n r ) − U n + 1 , j 0 ( x ¯ n r ) ] } = 0 , (13)

where λ n j s , k n j s , μ n j s are the Lagrange multipliers. It is necessary to find the spectrum of complex natural frequencies, ω k = ω R k + i ω I k where ω R k are the frequencies, and ω I k are the coefficients of damping of natural vibrations.

Approximate solution of vibrational Equation (13) is sought using an expansion in approximating forms composed of fundamental functions that satisfy the equation and the given geometric boundary conditions on Ω n f r surfaces of each body. It is assumed that functions Φ n j k ( x ¯ ) for such bodies are known (for rectangular plates and circular cylindrical shells this is a fundamental sequence of beam functions). Then the approximating forms can be constructed as a finite expansion in known functions [

U n j 0 ( x ¯ ) = ∑ k = 1 K γ n j k Φ n j k ( x ¯ ) , (14)

where γ n j k is the sought for complex coefficient.

Preliminary Φ n j k ( x ¯ ) can be normalized. The sum (14) automatically satisfies the boundary conditions on Ω n f r by virtue of terms choice. Varying Equation (13) on the generalized coordinates λ n j s , k n j s , μ n j s , γ n j s , a homogeneous system of linear equations is obtained. Dimension of this systems is J ⋅ N ′ × J ⋅ N ′ , where;

N ′ = ∑ n = 1 N ( S n + S n α + R n ) + N ⋅ K ; J is the number of components of the displacements vector U n j 0 ( x ¯ ) . Without specific calculations, we write this system in a matrix form:

( A + ∑ n = 1 N n f n ( ω R ) A n n + ∑ n = 1 N − 1 ∑ l = 1 L n f l n ( ω R ) A l n n + ∑ n = 1 N ∑ l ′ = 1 L n f l ′ n ( ω R ) A l ′ n n − ω 2 B ) ξ ¯ = 0 , (15)

where ξ ¯ is a column vector of generalized coordinates; N n -the number of viscoelastic bodies of the system; B is a symmetric, degenerate matrix of the generalized masses of the system; A n n , A l n n , A l ′ n n -square matrices of dimension J ⋅ N ′ × J ⋅ N ′ consisting of zeros, except for sub matrices of instantaneous stiffness of the n-the viscoelastic body, l-the shock absorber and l ′ -the viscoelastic support, respectively; A is a symmetric matrix (its sub matrix is A^{0} of dimension J ⋅ K × J ⋅ K ); it presents a generalized total stiffness of elastic elements of the system, and the sub matrices A H = A b T take into account kinematic conditions imposed on the system of rigid point relations.

Structurally, the matrices A and B are similar to those described in [

N ″ = J ⋅ N ′ − 2 J ∑ n = 1 N ( S n + S n α + R n ) . Equating the system determinant to zero, we obtain the frequency equation, which, unlike the case of elastic problem [

In engineering, an estimated logarithmic decrement of vibrations is assessed using ω I , and determined [

δ = − 2 π ω I ω R

In this section we will not consider the question of strict convergence of the method from a mathematical point of view, since it is not important for the following reasons. The energy approach used in the statement of the problem is essentially the Ritz method, the convergence of which is proved, for example, in [

In [

ω 1 = 2 π 2 a 2 D ρ h + 4 M / a 2 . (16)

Studies have shown that the first approximation (for supported on the edges of the plate) is quite close to the exact one. So, an increase in the number of members of fundamental functions to K = 16 “clarifies” the main frequency by only 1% - 2%; if the mass of the attached load is commensurate with the weight of the plate, the relative error is slightly larger (up to 6%) for K = 16.

At the edges of the free square plate is supported by four symmetrical supports located diagonally. For such problems, to determine the fundamental natural frequency, a theoretical and experimental formula is given [

ω 1 = ω ¯ a 2 D ρ h , (17)

where is the frequency coefficient determined experimentally depending on the location of the supports on the diagonals. Consider the two extreme cases—bearing on one support in the center and bearing in the corners (see

A square plate with a free circuit, spaced from the corners by a distance of r = 0.5 cm, is supported at four symmetric points:

ω ′ 1 = 246 ; ω ′ 2 = 1290.3 ; ω ′ 3 = 2815 ; ω ′ 4 = 3840.

When K = 25 in the formula (17), the calculated frequencies are as follows [

ω 1 = 250.7 ; ω 2 = 1380.3 ; ω 3 = 2839.8 ; ω 4 = 4077.3.

We reduce the initial problem to the following: a square plate is pivotally supported along the contour, one of its sides has discrete uniformly located

К | 6 | 9 | 16 | 25 | |
---|---|---|---|---|---|

Leaning in the corners | ϖ | 10.9082 | |||

ϖ K | 13.8061 | - | 11.6074 | 11.5811 | |

δ 1 % | 20.1833 | - | 6.9229 | 5.7265 | |

Leaning in the center | ϖ | 6.7806 | |||

ϖ K | - | 7.2712 | 7.2397 | 6.8214 | |

δ 1 % | - | 7.212 | 6.7030 | 0.6573 |

pinch points in the direction perpendicular to this side. We varied the number of fundamental functions of the sum (14) and the number of pinches that must satisfy (for a pivotally supported plate) the boundary conditions. The relative error was calculated by the formula [

δ i j = − ω i j − ω i j T ω i j 100 % ,

where ω i j T -is the natural frequency determined exactly; ω i j -natural frequency, determined approximately; i and j are the number of half-waves of the form along the axes 0X and 0Y, respectively.

Thus, the developed algorithm and the results obtained on the basis of the developed programs are reliable and the solution converges. The physical interpretation of the graph is as follows: the left branch of the graph is the energy of the plate, sufficient for the implementation of close to real (true) forms of vibration; curve minimum-saturation of the mechanical system with bonds; the right branch shows that the total energy of the plate decreases due to an excess of bonds, and, at higher frequencies, the shape of the vibrations is more and more distorted.

The examples given below are mostly theoretical in nature, but they allow us to draw practical conclusions. The problems are solved using the algorithm described in paragraph 3 of this paper.

Problem 1. An inhomogeneous system is considered, which consists of rectangular plates with point relations. The kernel of relaxation for deformable viscoelastic elements (shock absorbers) is chosen in the form of the Rzhanitsyn-Koltunov kernel [

R ( t ) = A e − β t t α − 1

where A , α , β are the kernel parameters [

In the general case, the methods to determine the kernel parameters (16) for various materials are given in [

In this case, in specific calculations, the kernel parameters are as follows:

A = 0.01 , α = 0.1 , β = 0.05 .

In contrast to the elastic problem, the dependence of the two lowest frequencies and corresponding damping coefficients on the magnitude of the instantaneous stiffness of the shock absorber was studied here. The stiffness varied from 10^{−4} to 10^{−1}. On the right, this range is limited to C_{2}, since at C = C_{2}, the second mode changes.

Real illustration of this effect is the presence of the point of graphs intersection of the damping coefficients ω I 1 and ω I 2 at C = C ∗ .

And one more feature: at this point, the difference in damping rates of two modes of structure vibrations changes the sign if the spring is viscoelastic (that is, the system is structurally inhomogeneous). As in the elastic problem, the Eigen frequencies ω R 1 , ω R 2 that are less than ω 1 , ω 2 at this point converge. Analysis of the problems of this type showed that the effect of the interaction of Eigen modes is observed only in structurally inhomogeneous systems (in this case with elastic and viscoelastic elements) and at a noticeable convergence of the real parts of the Eigen frequencies ω R 1 , ω R 2 ; the fictitious parts of the frequencies (the damping coefficients) change dramatically (increases or decreases). The absence of at least one of these conditions eliminates the manifestation of the effect.

Physical explanation of the observed effect should be sought in the nature of dynamic redistribution of energy of the system between two modes of the properties described above.

Now consider what differences exist between structurally homogeneous and inhomogeneous viscoelastic systems.

Structurally homogeneous viscoelastic system (all elements are viscoelastic with the same rheological properties) is characterized by the fact that, in formula (15), firstly, there is no matrix A (sub matrices A_{H} and A_{b} can be transferred to the next matrix) and, secondly, all functions f n ( ω R ) are identical. Then the system of Equations (15) in the matrix form can be rewritten as:

[ f ( ω R ) A n − ω 2 B ] ξ ¯ = 0 , (18)

where A^{n} is a numerical matrix of total instantaneous stiffness of all viscoelastic elements of the system.

After elimination of linearly dependent components from system (18), the transformed matrices of generalized instantaneous stiffness A ˜ n and B ˜ can be written in a canonical form, i.e. by special transformation of the generalized coordinates they are reduced to a matrix of diagonal form. This means that mechanical system is a set of independent partial systems with one degree of freedom. In other words, the Eigen modes of such a system are independent and can be considered and calculated separately.

Another situation is formed for a structurally inhomogeneous viscoelastic system, for which the matrix of generalized stiffness of elastic elements A is added to (18):

[ A ¯ ( ω R ) − ω 2 B ] ξ ¯ = 0 , (19)

where A ¯ ( ω R ) = A + f ( ω R ) .

In the general case, for two matrices А ( ω R ) and B (one of which is functional), after eliminating linearly dependent components, it is impossible to simultaneously select a non-degenerate coordinate transformation, leading to a canonical form. This means that natural modes of such a mechanical system cannot be considered separately from each other, i.e. they are interdependent. Consequently, under free oscillations, energy exchange occurs between the modes. This is seen when the modes have close natural frequencies. Then at the point of approaching of the graphs of natural frequencies ω R 1 , ω R 2 (at the point C = C ∗ ), the fictitious parts of the frequency ω I 1 and ω I 2 intersect.

At the intersection point of the graphs (

From the above, the following practical conclusion can be made: the damping capacity of inhomogeneous viscoelastic systems basically determines the minimum damping coefficient (in this case, the vibrations of this particular mode are damped out last); the system’s global (determining) damping coefficient is ω I 1 up to the intersection point and then ω I 2 . The optimal, in the sense of damping, vibration mode of the structure is at C = C ∗ , when this global damping coefficient has a maximal value.

Problem 2. A similar effect was found for a dissipative inhomogeneous viscoelastic structure consisting of two circular coaxially located cylindrical elastic shells interconnected by a massless viscoelastic element. Parameters of its relaxation kernel are [^{−2} to 10^{2}. Mechanical characteristics of the shells are the same and equal to E = 1 , ρ = 1 , ν = 0.35 . Geometric parameters of the first (inner) shell are L = 10 , R 1 = 1 , h 1 = 0.1 ; of the second (outer) shell: L = 10 , R 1 = 2 , h 1 = 0.2 . On the second shell there is a concentrated mass M = 0.5 attached at a point X M = 5 , Y M = 0 .

The shock absorber connecting the shells has the coordinates of the application X a = 5 , φ a = Y a 1 R 1 = Y a 2 R 2 = π 2 . Shells are hinge supported (Navier conditions) [

Dependence of the frequencies ω R and damping coefficients ω I on the rigidity of the shock absorber С_{0 }for the considered mechanical systems is shown in

This is explained by the fact that the deformable element less effectively damps vibrations of shells, i.e. torsional and longitudinal modes of vibrations remain undammed.

Unlike the previous problem, the global damping coefficient is, in the order of succession, not ω I 2 and ω I 1 , but ω I 1 and ω I 2 , respectively. The stiffness corresponding to the point of intersection of the graphs of damping coefficients is an optimal instantaneous stiffness of the shock absorber [

Problem 3. Consider a mechanical system consisting of two parallel, identical (in geometry and mechanical properties) elastic plates ( E = 2 × 10 11 N / m 2 ,

ρ = 7.8 × 10 3 kg / m 2 , ν = 0.35 , h = 0.001 m ), connected by one massless viscoelastic deformable element. The parameters of its relaxation kernel are A = 0.01 , α = 0.1 , β = 0.05 and instantaneous stiffness is C = 10.

The plates are of square form ( a = b = 0.2 m ), supported along the contour, each having one added mass ( M 1 = M 2 = 0.05 kg ) [

The aim of the study is to determine the nature of dependence of damping capacity of structure on the location of shock absorber and the masses on the area of the plates.

In this problem, the following option is investigated: a shock absorber is located in the center of the plate, the added masses are on the central axis of the structure, M_{1} is attached to the first (lower) plate at 0.04 m from the shock absorber; on the second (upper) plate a shock absorber is located in the center and positions of mass М_{2} on the second plate are varied along the central axis of the structure.

_{2}. Calculations showed a strong dependence of the damping coefficient of the first global mode ω I 1 on the location of the coordinates of mass М_{2}. The damping coefficient of the second mode ω I 2 and natural frequencies ω R 1 , ω R 2 remained practically unchanged.

The presence of two points on the x-axis, where the first global mode of the structure does not damp, is explained by the fact that the wave modes of the

plate are indistinguishable in this case. When mass М_{2} is located at these points, the plates become identical in terms of inertial and stiffness characteristics [

If to move the mass М_{2} on the upper plate in other directions, then the points could be found at which the damping coefficient ω I 1 is zero. Thus, for mass М_{2}, there are an infinite number of positions (at fixed positions of the shock absorber and mass M_{1}) in which the first global mode does not damp. These points on the second plate form a closed curve (close to a circle).

Analysis of the results (

The effect described shows that the magnitude of dissipative energy of the system depends on damping (rheological) properties of the material and on structure geometry as a whole. This effect does not manifest itself if the viscoelastic system is dissipative homogeneous [

Problem 4. Consider the structure with the same initial data as in Problem 3. The masses are located at distances of 0.04 m on both sides of the massless deformable element. The Eigen frequencies of two modes of vibrations are almost unchanged. Consider the dependence of frequency and damping characteristics of the system on simultaneous displacement of a massless deformable element and mass.

As in the previous example, there is a point (when the shock absorber is in the center) where ω I 1 = 0 . At this point, the second damping coefficient is a maximal one, i.e. ω I 2 . In the presence of concentrated masses the dependence of damping properties of the system is so strong that, in the absence of them, the damping coefficient ω I 1 (global one) is zero for any location of the shock absorber (the corresponding graph is shown in

Knowledge of the law of variation determining damping coefficient ω I 1 made it possible to make various recommendations of an optimization nature. The optimization of the damping properties of structure did not affect the spectrum of the lowest frequencies. The instantaneous stiffness of the shock absorber, the magnitude of the added mass and viscosity did not vary, since these parameters did not (qualitatively) affect the results obtained. Only the quantitative characteristics of the system varied, which is not fundamental.

The effect revealed in this problem does not manifest itself if to consider a dissipatively homogeneous viscoelastic system.

After analyzing problems 3 and 4, we can conclude that for given damping coefficients of the system, it is possible to determine the necessary rheological properties of the shock absorber, the size and location of the added masses, its instantaneous stiffness and the location of masses and shock absorber in the plane of the plates [

Problem 5. Consider a viscoelastic shell, the rheological properties of which are described by a kernel of the form given in (17); the parameters are: A = 0.01 , α = 0.1 , β = 0.05 . Geometric and mechanical characteristics are: E = 1 , ρ = 1 , ν = 0.35 , R = 1 , L = 7 , h = 0.01 . Coordinates of rigid support are: x = 2 , y = 3 π / 2 . Shell thickness h varied from 0.01 to 0.2. Results of calculations for the two lower frequencies ω R 1 and ω R 2 of corresponding damping coefficients ω I 1 , ω I 2 are shown in

For comparison, the first two Eigen frequencies of elastic shell ω 1 , ω 2 are given. A twentyfold increase in the shell thickness increases the first and second Eigen frequencies and the damping coefficient by almost 4 times.

1) The mathematical formulation and methods for solving the problem of natural vibrations of structurally inhomogeneous dissipatively homogeneous and heterogeneous mechanical systems consisting of a package of plates (or shells) with point supports and attached masses are formulated.

2) The convergence of solution is numerically proved depending on the terms taken in the sought-for solution (i.e., in solution expansion in fundamental functions) and on various parameters of mechanical systems.

3) The problems of natural vibrations of dissipative homogeneous and inhomogeneous lamellar systems with concentrated supports and added masses are solved.

h | ω 1 | ω R 1 | − ω I 1 | ω 2 | ω R 2 | − ω I 2 |
---|---|---|---|---|---|---|

0.01 | 0.0927 | 0.0872 | 0.4 × 10^{−4} | 0.142 | 0.135 | 0.94 × 10^{−4} |

0,1 | 0.274 | 0.26 | 0.34 × 10^{−3} | 0.281 | 0.27 | 0.44 × 10^{−3} |

0.2 | 0.339 | 0.322 | 0.53 × 10^{−3} | 0.502 | 0.48 | 0.11 × 10^{−2} |

4) To describe dissipative properties of the system as a whole, the concept of “global damping coefficient” is introduced. In the case of a structurally homogeneous mechanical system, the global damping coefficient is determined by the fictitious part of the first (in modulus) complex Eigen frequency. In the case of a structurally inhomogeneous mechanical system, the fictitious parts of the first and second Eigen frequencies play the role of a global damping coefficient, depending on the stiffness magnitude of a shock absorber.

5) It is established that optimal damping of vibrations of inhomogeneous systems occurs when the values of close frequencies of vibration modes converge, and the fictitious parts of these frequencies are equal to each other (i.e., at C=C*). In this case, both modes of vibrations provide the same energy dissipation.

6) It is revealed that for given damping coefficients of the system, it is possible to determine the necessary rheological properties of a shock absorber, the size and location of the added masses in the plane of the plates, to ensure effective control of system vibrations with maximum damping capacity.

1) To reduce the amplitude of the interference of radio electronic equipment in resonance mode, we recommend that we use the developed methodology and algorithm of structurally inhomogeneous dissipative mechanical systems, which effectively reduces the amplitudes of mixing and other force factors to 60%.

2) Significant intensification of dissipative processes in structurally heterogeneous dynamical systems occurs when the corresponding natural frequencies approach each other.

3) The effect of the interaction of various forms of motion of solid bodies has a fundamental perspective for the synthesis of structurally heterogeneous machine-building structures that are optimal in dissipative properties and material consumption.

4) Analytically, to obtain several eigenfrequency modes of structurally heterogeneous dynamic mechanical systems, depending on the physical and geometric parameters, it remains open.

The authors declare no conflicts of interest regarding the publication of this paper.

Mirsaidovich, M.M., Ibrokhimovich, S.I. and Khudoyberdievich, T.M. (2019) Dynamics of Structurally Inhomogeneous Lamellar and Shell Mechanical Systems. Part 1. Journal of Applied Mathematics and Physics, 7, 2283-2302. https://doi.org/10.4236/jamp.2019.710155