The simultaneous action of external excitations – features of the dynamical damping of vibrations
Sergey Eliseev^{1} , Alexey Orlenko^{2} , Andrey Eliseev^{3}
^{1, 3}Irkutsk State Transport University, Irkutsk, Russia
^{2}Krasnoyarsk Institute of Railway Transport – Branch of Irkutsk State Transport University, Krasnoyarsk, Russia
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 18, 2018, p. 1-7.
https://doi.org/10.21595/vp.2018.19893
Received 7 April 2018; accepted 18 April 2018; published 22 May 2018
The article considers new dynamic effects arising in problems of vibration protection of objects under external disturbances, on the side of support surfaces. A method for constructing mathematical models for a vibration protection system is developed, in which the protection object has two degrees of freedom. It is demonstrated that is possible to introduce additional constraints in the form of motion transformation devices. The article estimates the possibilities of using screw non-locking mechanisms with flywheel nuts, the given moment of inertia of which can be regulated.
Keywords: frequency diagram, simultaneous vibration damping modes, transfer functions.
1. Introduction
In the problems of the dynamics of technological machines and vehicles, much attention is paid to the issues of vibrational interactions of elements of mechanical oscillation systems (MOS) [1-4]. Effects of dynamic damping of oscillations are widely used in various branches of technology [5-8], which is associated with the study of the features of the dynamic properties of objects [9-11], most often selected as a mass-and-inertia element, whose dynamic state is estimated by one coordinate.
In works concerned with controllable systems as part of the MOS, servo drives, elements of pneumatic, hydro-, and electrical automation are used [12-14]. Although the problems of analysis and dynamic synthesis of the MOS have been in progress, the methods of structural mathematical modeling have been developed to a lesser extent, oriented toward using the analytical tools of the automatic control theory, which provides great possibilities for estimating the dynamic states of objects having several degrees of freedom [2, 3, 9, 15]. The main provisions of structural mathematical modeling are given in [8, 13, 16-18].
The proposed article develops a methodological basis for constructing mathematical models for vibration protection systems with an object of protection in the form of a solid body that performs a plane motion under the influence of a kinematic perturbation, taking into account that it’s possible to implement dynamic damping modes simultaneously in two coordinates.
2. Some general provisions. Statement of the research task
We consider a vibration protection platform for the installation of instrumentation, which is subject to vibrations on the side of the supporting surface (Fig. 1). The object in the form of a solid body with mass $M$ and moment of inertia $J$ is supported by elastic elements with stiffnesses ${k}_{1}$ and ${k}_{2}$ with parallel-mounted motion transformation devices (MTD) having reduced masses ${L}_{1}$ and ${L}_{2}$, respectively [11, 17].
The system under consideration has linear properties and performs small oscillations with respect to the static equilibrium position. The position of the center of gravity (p. O) is determined by the lengths of the arms ${l}_{1}$ and ${l}_{2}$. The ${y}_{1}$ and ${y}_{2}$ coordinate system associated with the fixed base is used. The frictional forces in the system are not taken into account, the external action $z\left(t\right)$ is a harmonic function.
Fig. 1. Schematic diagram of the vibration protection platform in the form of a rigid body with elastic-inertial supports with kinematic perturbation $z\left(t\right)$
To obtain a mathematical model, using the Lagrange equations of the second kind, an expression is determined for the kinetic and potential energies:
To carry out the appropriate calculations, we take into account the following relations:
$a=\frac{{l}_{2}}{{l}_{1}+{l}_{2}},b=\frac{{l}_{1}}{{l}_{1}+{l}_{2}},c=\frac{1}{{l}_{1}+{l}_{2}}.$
Using the Eq. (1)-(3), on the basis of known methods we obtain a mathematical model in the form of a system of two linear inhomogeneous differential equations with constant coefficients. After the Laplace transform under zero initial conditions, these equations in the operator form become:
where $p=j\omega $ is the complex variable ($j=\sqrt{-1}$); icon “–” corresponds to the Laplace transform of the variable.
On the basis of Eqs. (4), (5) a structural mathematical model can be constructed, as shown in Fig. 2; the system has an external perturbation simultaneously at the inputs of two partial systems.
Fig. 2. Structural mathematical model of the system by Fig. 1 in the form of a block diagram of an equivalent (dynamically) automatic control system
Using the block diagram in Fig. 2, we write the transfer functions of the system:
where:
is the characteristic frequency equation.
We note that ${\stackrel{-}{y}}_{1}$ and ${\stackrel{-}{y}}_{2}$ can take zero values if condition:
At a frequency determined by the relation ${k}_{1}{L}_{2}={k}_{2}{L}_{1}$, it is possible to simultaneously “nullify” the coordinates ${y}_{1}$ and ${y}_{2}$, which can be considered as a mode of dynamic damping of the oscillations of the protection object simultaneously in the two coordinates ${y}_{1}$ and ${y}_{2}$ under the kinematic perturbation from the support surface side.
The task of the research is to study the possibilities of forming modes of dynamic damping of the object’s oscillations simultaneously in two coordinates, using data on the available forms of a purposeful change in the ratio of the parameters of the systems for configuring the parameters of the dynamic state of the object.
3. Construction of mathematical models of interactions
Using relations Eq. (9), assuming that ${L}_{2}/{L}_{1}=\alpha $, we write the transfer functions Eqs. (6), (7) in the form:
where:
We note that in the construction of transfer functions of the system, simultaneous transmission of influences to both inputs of the system is taken into account, while $\alpha $ is regarded as an adjustable connectivity coefficient of the parameters of the reduced masses of the MTD.
The partial frequencies of the system are determined by the expressions, respectively:
Note that the partial frequency ${n}_{1}^{2}$ does not depend on α.
In the MOS under consideration (Fig. 1), it is possible to create dynamic oscillation damping modes in two coordinates, which is determined by the possibility of “nullifying” the numerators of the transfer functions Eqs. (10), (11).
For the coordinate ${\stackrel{-}{y}}_{1}$, we can write the following equation for finding the frequencies of the dynamic damping of oscillations under the kinematic perturbation $\stackrel{-}{z}$:
$+{k}_{2}\left({L}_{1}+J{c}^{2}-Mab\right)+{k}_{1}{k}_{2}=0..$
Eq. (15) can be reduced to the form of a biquadratic equation, assuming that ${p}^{2}=-{\omega}^{2}=-\gamma $, then Eq. (15) takes the form:
where:
${R}_{2}={k}_{1}\left(M{b}^{2}+J{c}^{2}+\alpha {L}_{1}\right)+{k}_{2}\left({L}_{1}+J{c}^{2}\u2013Mab\right),{R}_{3}={k}_{1}{k}_{2}.$
The solution of Eq. (15) can have, with a certain set of parameters of the system, two real positive roots:
In turn, from the coordinate we obtain that:
$+{k}_{1}(\alpha {L}_{1}+J{c}^{2}-Mab)+{k}_{1}{k}_{2}=0.$
where:
${R}_{2}^{\text{'}}={k}_{2}\left(M{a}^{2}+J{c}^{2}+{L}_{1}\right)+{k}_{1}\left(\alpha {L}_{1}+J{c}^{2}\u2013Mab\right),{R}_{3}^{\text{'}}={k}_{1}{k}_{2}.$
The introduction of the connectivity coefficient α between the parameters ${L}_{1}$ and ${L}_{2}$ changes the transfer functions Eqs. (6), (7), and hence the values of the frequency of the dynamic damping of the oscillations determined by Eqs. (16) and (19). In this case, the characteristic frequency equation also changes, that is, the frequencies of proper oscillations also change accordingly.
Taking into account the connectivity coefficient $\alpha $, we shall use, using Eqs. (15) and (18), an equation for determining the frequencies of dynamic damping, ${\omega}_{1din}^{2}$, ${\omega}_{2din}^{2}$.
The frequencies of the natural oscillations ${\omega}_{1sob}$, ${\omega}_{2sob}$ can be found from the solution of the characteristic Eq. (12), which can be represented in the form of:
$+{p}^{2}\left[{k}_{2}\right(M{a}^{2}+J{c}^{2}+{L}_{1})+{k}_{1}(M{b}^{2}+J{c}^{2}+\alpha {L}_{1}\left)\right]+{k}_{1}{k}_{2}=0.$
We note that the frequencies of the natural oscillations in this case will depend on the connectivity coefficient $\alpha $.
4. Discussion of the results
Under the kinematic perturbation of the initial system (Fig. 1), as shown in the structural scheme shown in Fig. 2, an external disturbance is distributed simultaneously over two inputs. The use of the principle of superposition results in a change in the form of transfer functions in comparison with conventional approaches, when only one input is perturbed.
Under the action of one external perturbation in a system with two degrees of freedom, the numerator of the transfer function along the coordinate of the application of force. It is formed on the basis of the parameters of the corresponding partial system. According to another coordinate, the numerator of the transfer function is created due to the parameters of the link that implements interpartial bonds, which, on the whole, depends on the choice of the coordinate system. In the standard situation, the MOS in Fig. 1 in the coordinates ${y}_{1}$, ${y}_{2}$ for a force perturbation with respect to the coordinate ${y}_{1}$ will have one mode of dynamic damping of the oscillations, at a frequency determined by the expression:
In the second coordinate ${y}_{2}$ with force perturbation, the dynamic damping mode is not assumed, but due to the specificity of the transfer function of the link of the interpartial bond with the transfer function:
the coordinate ${y}_{2}$ can acquire zero values when the condition $J{c}^{2}=Mab$ is fulfilled.
Using formulas for determining the frequencies of partial and natural oscillations, as well as the frequencies of dynamic damping of oscillations, it is possible to construct a frequency diagram that reflects the features of the dynamic interactions of the initial MOS, as shown in Fig. 3.
To specify the representations, a model problem with parameters $M=$ 10 kg was solved; $J=$ 5 kg.m^{2}; $a=$ 0.4 m; $b=$ 0.6 m; $c=$ 1; ${k}_{1}=$ 5000 N/m; ${k}_{2}=$ 10,000 N/m; ${L}_{1}=$ 10 kg; ${L}_{2}=\alpha {L}_{1}$.
Fig. 3. Diagram of dynamic interaction modes ((${L}_{1}$, ${L}_{2}=\alpha {L}_{1}$): the notation of the dependencies ${\omega}_{1din}^{2}\left(\alpha \right)$, ${\omega}_{2din}^{2}\left(\alpha \right)$, ${\omega}_{1sob}^{2}\left(\alpha \right)$, ${\omega}_{2sob}^{2}\left(\alpha \right)$, ${n}_{1}^{2}\left(\alpha \right)$ and ${n}_{2}^{2}\left(\alpha \right)$ are given directly in the text
In the diagram (Fig. 3) the solid line (––) indicates the dependency ${\omega}_{1din}^{2}\left(\alpha \right)$ graphs, the numerator of the transfer function Eq. (10) is used to construct the ${\omega}_{1din}^{2}\left(\alpha \right)$ graph. The value of the frequencies ${\omega}_{1din}^{2}$ depends on the connectivity coefficient $\alpha $. Since the frequency of the dynamic damping of the oscillations along the ${y}_{1}$ coordinate is determined from the biquadratic frequency equation, the graph is represented in the diagram (Fig. 3) by two fragments (or branches): the lower branch of the graph ${\omega}_{1din}^{2}\left(\alpha \right)$ corresponds to low-frequency forms of motion. In turn, the plot of the frequency of the dynamic damping of the oscillations ${\omega}_{2din}^{2}\left(\alpha \right)$ along the ${y}_{1}$ coordinate is determined from the frequency equation formed by the numerator of the transfer function Eq. (11). In the diagram (Fig. 3), both branches of the graph ${\omega}_{2din}^{2}\left(\alpha \right)$ are indicated by a dashed line (- - -).
In the case under consideration, the partial frequencies of the system determined by Eqs. (13), (14) are also represented in Fig. 3 by graphs ${n}_{1}^{2}\left(\alpha \right)$, ${n}_{2}^{2}\left(\alpha \right)$ and are denoted by dash-dotted lines (∙-∙-). The graph of the dependence ${n}_{1}^{2}\left(\alpha \right)$ is a straight line parallel to the abscissa axis, since ${n}_{1}$, as it follows from Eq. (13), and does not depend on $\alpha $. The graphs of the dependences ${\omega}_{1sob}^{2}\left(\alpha \right)$ and ${\omega}_{2sob}^{2}\left(\alpha \right)$ can be constructed based on the use of the frequency characteristic Eq. (12). In Fig. 3 graphs ${\omega}_{1sob}^{2}\left(\alpha \right)$ and ${\omega}_{2sob}^{2}\left(\alpha \right)$ are denoted by dotted lines (∙∙∙∙).
5. Conclusions
Technically, such capabilities can be implemented with the help of special mechanisms controlled on the basis of measuring the current parameters of the dynamic state of the object of vibration protection. The authors suggest a method for constructing mathematical models of the system, based on the use of transfer functions of the system. The proposed analytical tools rely on the method of structural modeling, in which the mechanical oscillatory system is compared with a dynamically equivalent structural diagram of the automatic control system.
References
- Chelomei V. N. Vibrations in Technology: a Handbook in 6 Volumes. Protection from Vibration and Shocks, Vol. 6, Mashinostroenie Publishers, Moscow, 1981, p. 456, (in Russian). [Search CrossRef]
- Harris S. M., Piersol A. G. Shock and Vibration Handbook. McGraw-Hill Book Co, New York, 2002, p. 1457 [Search CrossRef]
- De Silva C. W. Vibration. Fundamentals and Practice. CRC Press, Boca Raton, London, New York, Washington, 2000, p. 957, (in Russian). [Search CrossRef]
- Belokobylsky S. V., Eliseev S. V., Kashuba V. B. Applied Problems of the Structural Theory of Vibration Protection Systems. Polytekhnika, St. Petersburg, 2013, p. 363, (in Russian). [Search CrossRef]
- Eliseev S. V., Nerubenko G. P. Dynamic Vibration Dampers. Nauka Publishers, Novosibirsk, 1982, p. 144, (in Russian). [Search CrossRef]
- Korenev B. G., Reznikov L. M. Dynamic Vibration Dampers: Theory and Technical Applications. Nauka Publishers, Moscow, 1988, p. 304, (in Russian). [Search CrossRef]
- Karamyshkin V. V. Dynamic Damping of Oscillations. Mashinostroenie, Leningrad, 1988, p. 108, (in Russian). [Search CrossRef]
- Eliseev S. V., Khomenko A. P. Dynamic Damping of Vibrations: the Concept of Feedback and Structural Methods of Mathematical Modeling. Nauka Publishers, Novosibirsk, 2014, p. 357, (in Russian). [Search CrossRef]
- Kolovsky M. Z. Automatic Control of Vibration Protection Systems. Nauka Publishers, Moscow, 1976, p. 320, (in Russian). [Search CrossRef]
- Eliseev S. V. Structural Theory of Vibration Protection System. Nauka Publishers, Novosibirsk, 1978, p. 221, (in Russian). [Search CrossRef]
- Eliseev S. V., Lukyanov A. V., Reznik Yu N., Khomenko A. P. Dynamics of Mechanical Systems with Additional Ties. Irkutsk State University, Irkutsk, 2006, p. 315, (in Russian). [Search CrossRef]
- Eliseev S. V., Kuznetsov N. K., Lukyanov A. V. Controlling the Vibrations of Robots. Nauka Publishers, Novosibirsk, 1990, p. 320, (in Russian). [Search CrossRef]
- Eliseev Reznik Yu S. V. N., Khomenko A. P. Mechatronic Approaches to the Dynamics of Mechanical Oscillation Systems. Nauka Publishers, Novosibirsk, 2011, p. 384, (in Russian). [Search CrossRef]
- Kuznetsov N. K. Dynamics of Controlled Machines with Additional Constraints. IrGTU, Irkutsk, 2009, p. 290, (in Russian). [Search CrossRef]
- Chernousko F. L., Akulenko L. D., Sokolov B. N. Control of Oscillations. Nauka Publishers, Moscow, 1980, p. 384, (in Russian). [Search CrossRef]
- Khomenko A. P., Eliseev S. V., Ermoshenko Yu V. System Analysis and Mathematical Modeling in the Mechatronics of Vibration Protection Systems. ISTU, Irkutsk, 2012, p. 288, (in Russian). [Search CrossRef]
- Eliseev S. V., Artyunin A. I. Applied Theory of Oscillations in Problems of the Dynamics of Linear Mechanical Systems. Nauka Publishers, Novosibirsk, 2016, p. 459, (in Russian). [Search CrossRef]
- Karnovsky I. A., Lebed E. Theory of Vibration Protection. Springer, Switzerland, 2016, p. 708. [Publisher]