Dynamics of Structurally Inhomogeneous Lamellar and Shell Mechanical Systems. Part 1

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.

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 [1] [2] developed exact solutions for problems with free vibrational vibrations of rectangular plates.
Boay [3] analyzed the natural frequencies of plates with and without concentrated masses using the Rayleigh energy method.
Avalos et al. [4] dealt with solving vibration problems using a simple mounted concentrated mass using the well-known normal mode.
Xiang et al. [5] used the Ritz method in combination with the variation method to solve the problem of vibration of rectangular plates, which rely on elastic edge supports.
Laura and Grossi [6] [7] calculated the fundamental frequency coefficient for a rectangular plate with edges elastically held against both displacements and rotations using polynomial coordinates and the functions of the Rayleigh-Ritz method.
Wu and Luo [8] solved the problem of eigenfrequencies and corresponding mode shapes of a uniform rectangular plane of a plate carrying any number of point masses using analytical and combined numerical-analytical methods.
Nicholson and Bergman [9] used the Green's function in describing natural modes for damping oscillatory systems. Gorman [10] solved the problem of free vibrations of deformable systems during shear of plates based on homogeneous elastic substrates using the modified Galerkin method. This work is devoted to the application of the differential quadrature method to the vibration of plates resting on elastic bases and bearing any number of spring-loaded masses.
The design of submarine tunnels is inevitably characterized by the presence of attached masses due to structural, strength or operational necessity [11] [12] [13]. Such a system "Shell-mass" in operating conditions is subjected to intense dynamic, in particular, periodic loads, which requires increased accuracy of calculations [14].
Rod systems have a high degree of vibration isolation. E. Skudrzyk [15] considered vibration isolation of a structure consisting of masses and rods connected in series. Lamellar mechanical systems with point supports and lumped masses are also dedicated to the work [16]- [24].
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 [25] [26]. Free vibrations of a dissipative system, which are of a damped nature, were studied in [27] [28] [29] [30] [31]. The amplitudes of the vibration modes decrease over time, therefore, such a process, strictly speaking, is not periodic. But the frequencies of the corresponding forms remain constant, and in this sense, the dissipative system can be investigated as a system with its own vibrations.
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. [32].
In the works of I. I. Safarov  However, to date, theoretical results are not always consistent with numerical and experimental data [38]. For example, in [39] it was found that the decrease in the smaller of the split eigenfrequencies is stronger than that predicted by the traditional theory [40], but also depends on the geometric parameters of the shell [41]. This circumstance requires clarification of the traditional mathemati- 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. 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. 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.

Vibrational Statement of the Problem of Natural Vibrations for Structurally Inhomogeneous Plate and Shell Mechanical Systems with Point Relations and Concentrated Masses
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.
( ) ( ) 0 , e , 1, , ; 1, , 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 where , , 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:  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 [14] ( ) ( ) ( ) , instead of (4) we get where n Е  is the Volter operator of the form [42] [43]: here 0n E is the instantaneous modulus of elasticity, and En R is the kernel of relaxation.
Given (1), the time function in Equality (6) using the freeze-etch method [38], the relation (6) is substituted by the approximate one: where This allows eliminating the integral terms and, ultimately, time from the vibrational equation. In symbolic form, it can be represented as Write out a specific representation of the functional G, for example, for a package of rectangular plates with point relations: is a complex function whose numerical coefficients depend on the parameters of the relaxation kernel of corresponding viscoelastic elements, If the part of supports is fixed, then the following conditions are added.
where s n α is the direction of the unit vector along which the rigid fixing of the body is done at the point s n x .
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 2 N ≥ is taken into account by the relationships where r x is the coordinate of the r-the post, n R is the number of posts between the n-the and (n + 1)-the bodies. At 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 where , ,

Algorithm for the Implementation of the Vibrational Method to Solve Viscoelastic Problem 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 fr n Ω surfaces of each body. It is assumed that functions 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 [35] [36]: where k nj γ is the sought for complex coefficient.
can be normalized. The sum (14) automatically satisfies the boundary conditions on 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; , , n n n n ln l 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 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 [17] for elastic mechanical systems. The generalized stiffness and masses of the n-th elastic element are calculated using the method given in [17], if the elements of the system are rectangular plates. If all viscoelastic elements of the system have the same rheological properties, then ( ) ( ) 1 2 R R f f ω ω = =  so, the second, third and fourth terms in (15) are replaced by one matrix of total instantaneous stiffness of all viscoelastic elements (in case of structurally homogeneous viscoelastic system). The degeneracy of the matrix B, as in the elastic problem [17], is due to the introduction of additional point relations (rigid supports and posts) into the system. By the method described in [17], we eliminate linearly dependent coordinates and bring system (15)

Evaluation of the Practical Convergence of the Algorithm and the Reliability of Numerical Results on the Problems of Natural Vibrations of Dissipative Mechanical Systems
In this section we will not consider the question of strict convergence of the me-thod 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 [25] [28]. Accounting with the help of Lagrange multipliers of point bonds superimposed on plates is also a well-known method for finding a conditional extremum. Superimposed point relationships affect the rate of convergence, but not the convergence itself. To illustrate the convergence of the method, we compare the data available in the literature (theoretical and experimental) with the results of calculations by the described method that we obtained.
In [25] [28], a relation was obtained that allows (in a first approximation) to calculate the first (main) frequency of a square elastic plate supported at the edges, with an attached mass M in the center: 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 [25] 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 Table  1). Data analysis indicates satisfactory convergence of the method and the dependence of the convergence rate on the number of supports is observed. The more point dependencies, the more the convergence rate worsens, i.e. the speed of convergence depends on the number of supports. The more point bonds are superimposed on the plate, at a fixed row length (14), the greater the error. In [28] for symmetric natural vibrations of a duralumin square plate (a = 13 cm, h = 0.193 cm) experimental results of the first frequency.
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: 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 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 [28]: where T ij ω -is the natural frequency determined exactly; ij ω -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.

Solution and Analysis of the Problem of Natural Vibrations of Structurally Inhomogeneous Viscoelastic Systems
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 [26] ( ) where , , A α β are the kernel parameters [24]. The viscosity of the shock absorber is taken such that its creep strain during a quasistatic process is a small fraction (~12%) of the total strain.
In the general case, the methods to determine the kernel parameters (16) for various materials are given in [26], and in [33], using these methods, specific parameters for some materials are given.
In this case, in specific calculations, the kernel parameters are as follows: 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. Figure 1 shows the dependence of the first two frequencies 1 2 ,

R R
ω ω and corresponding damping coefficients on the magnitude of instantaneous stiffness of the shock absorber C. From the analysis of graphs, it follows that the dissipative properties of this system are determined not only by the rheology of its elements, but greatly depend on the interaction of natural modes of vibration. This effect is expressed in the fact that under certain conditions and to some value of the rigidity of the shock absorber, the mode energetically more capacious (in this case, the second one) dissipates less energy than the mode less energy-intensive. Then, beginning from some value of instantaneous stiffness of the shock absorber C * (in this case,  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 1 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 1 2 , R R ω ω ; 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 ( ) n R f ω are identical. Then the system of Equations (15) in the matrix form can be rewritten as: 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 n A  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): 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 1 2 ,

R R
ω ω (at the point C C * = ), the fictitious parts of the frequency 1 I ω and 2 I ω intersect.
At the intersection point of the graphs (Figure 1) of damping coefficients 1 I ω and 2 I ω both modes dissipate energy in the same way, although they differ from each other (up to a phase). Up to the point C C * = there is a "transfer" of energy from the second mode to the first one, therefore the latter dissipates energy more intensively. Past the intersection point, the difference between the first Eigen frequencies increases, the interaction of corresponding modes decreases, and their dissipative properties take on an ordinary character.
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 1 I ω , respectively. The stiffness corresponding to the point of intersection of the graphs of damping coefficients is an optimal instantaneous stiffness of the shock absorber [27].
Problem 3. Consider a mechanical system consisting of two parallel, identical (in geometry and mechanical properties) elastic plates (  [29]. 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. Figure 3 shows  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 [30]. As in the previous example, there is a point (when the shock absorber is in the center) where 1 0 I ω = . At this point, the second damping coefficient is a maximal one, i.e. 2 I ω . 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 1 I ω (global one) is zero for any location of the shock absorber (the corresponding graph is shown in Figure 3).
Knowledge of the law of variation determining damping coefficient 1 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 [27].
Problem 5. Consider a viscoelastic shell, the rheological properties of which are described by a kernel of the form given in (17) ω ω are shown in Table 2.
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.

Conclusions
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. Journal of Applied Mathematics and Physics  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.

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