Next Article in Journal
A Procedure for Tracing Supply Chains for Perishable Food Based on Blockchain, Machine Learning and Fuzzy Logic
Next Article in Special Issue
Usage of Mathematical Modeling and Optimization in Development of Hydrogel Medical Dressings Production
Previous Article in Journal
A Hybrid Prognostics Deep Learning Model for Remaining Useful Life Prediction
Previous Article in Special Issue
The Sentiment Analysis Model of Services Providers’ Feedback
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Methodology for Increasing the Efficiency of Dynamic Process Calculations in Elastic Elements of Complex Engineering Constructions

1
Scientific and Technical Library, Lviv Polytechnic National University, 79007 Lviv, Ukraine
2
Department of Transport Technologies, Institute of Engineering Mechanics and Transport, Lviv Polytechnic National University, 79007 Lviv, Ukraine
3
Department of Social Communication and Information Activities, Institute of the Humanities and Social Sciences, Lviv Polytechnic National University, 79007 Lviv, Ukraine
4
Department of Philology, Interpreting and Translation, Institute of Humanities and Public Administration, Ivano-Frankivsk National Technical University of Oil and Gas, 76019 Ivano-Frankivsk, Ukraine
5
Department of Information Systems, Faculty of Management, Comenius University in Bratislava, 820 05 Bratislava, Slovakia
*
Author to whom correspondence should be addressed.
Electronics 2021, 10(1), 40; https://doi.org/10.3390/electronics10010040
Submission received: 27 November 2020 / Revised: 17 December 2020 / Accepted: 22 December 2020 / Published: 29 December 2020
(This article belongs to the Special Issue Electronization of Businesses - Systems Engineering and Analytics)

Abstract

:
This study deals with a methodology for increasing the efficiency of dynamic process calculations in elastic elements of complex engineering constructions. We studied the complex dynamic processes in a simple engineering construction, a mechanical system of an elastic body–continuous flow of homogeneous medium. The developed methodology is based on the use of a priori information on some of the vibrations forms, the construction of a “simplified” mathematical model of system dynamics, and the obtaining of an analytical relationship that describe the overall range of factors on the elastic vibrations of system. The methodology is used for cases of complex vibrations of elastic bodies, and the obtained results can serve as a basis for choosing the main technological and operational parameters of elastic elements of mechanisms and machines that perform complex vibrations. The results obtained in this work are the basis for calculating the blast effect on the elements of protective structures in order to increase their protective capacity by improving the method of their attachment or by using additional reinforcement, buff load effects on the elements of drilling strings and dynamic processes that occur during surface strengthening by work hardening in order to avoid resonance phenomena, and technological processes of vibration displacement or vibration separation of granular media.

1. Introduction

A methodology for researching complex vibrations of the structural elements of the simplest engineering constructions has been developed. It is based on the usage of a priori information about one of the vibration forms of an elastic body, its processing and usage during the construction of a mathematical model of another form of vibrations, and development or usage of existing analytical methods for constructing a solution for the latter one. The methodology, however, simplifies the mathematical model of process dynamics and thus simplifies the procedure of finding and researching its solution. Using this methodology enables one to establish a number of features of the structural elements’ dynamics, in particular the conditions of external and internal resonances and their initial peculiarities. The main idea of the methodology is presented using the example of an elastic body that performs bending and longitudinal vibrations, and in addition to that a continuous flow of homogeneous medium (CFHM) of zero stiffness that moves along the body, which greatly complicates the problem of constructing a solution for the mathematical model of elastic body dynamics. For this reason, the following issues were predetermined:
  • even movement at a constant speed along the elastic body of the CFHM affects the natural frequencies of the component vibrations;
  • the most dangerous (resonant) vibrations of the specified type of systems can be caused both by external periodic forces and mutual influence of certain forms of vibrations on the others (internal resonances);
  • internal resonances of an elastic body that performs bending and longitudinal vibrations can exist with even and odd modes of bending vibrations; with even modes they are caused by its longitudinal vibrations and with odd modes by its longitudinal vibrations and continuous flow of homogeneous medium motion;
  • the amplitude of bending vibrations of the transition through the internal resonance at higher frequencies of longitudinal vibrations is less than at the second frequency;
  • the amplitude of the transition through the main resonance for larger values of the amount of motion of the continuous flow of homogeneous medium is smaller than the amplitude of the transition during the rapid transition through the resonance.
The methodology can be generalized in other cases of complex vibrations of elastic bodies, and the obtained results can serve as the basis for choosing the main technological and operational parameters of elastic elements of mechanisms and machines that perform complex vibrations. The reliability of the obtained calculated dependences is confirmed by the boundary transition from which follow the results known from scientific publications concerning the one form of vibrations of elastic bodies.

2. Problem Statement and Research Relevance

Analysis of dynamic processes in the elements of machines and structures, which are elastic bodies, shows that in most cases they carry out such complex oscillations through a combination of longitudinal forces, bending, and torsion. Isolation during research of some of them, even the largest amplitude with simultaneous neglect of the influence of others (even small amplitude) will lead to significant inaccuracies in describing the dynamic process of the element and the system as a whole. The value of the error increases significantly in cases where the natural frequencies of the oscillations are close in value or related by a rational ratio, and the elastic properties of the object under study are described by nonlinear relations.
The analytical study of complex nonlinear dynamic processes in elastic bodies (construction elements) based on mathematical models that take into account a wide range of external and internal factors is associated with significant mathematical difficulties. This is the integration of systems of nonlinear non-autonomous boundary value problems for equations with partial derivatives. Therefore, there is an actual problem of development of such an analytical methodology that takes into account a broad range of factors that impact on the object and receive its calculation-dependent base that would be available for engineering calculations at the design and upgrade stage. A partial solution for these problems is illustrated with the case of an elastic element mechanism (construction), which provides both longitudinal and bending vibrations. In addition, the elastic element is characterized by longitudinal frequencies, or it moves along the continuous flow of an incompressible homogeneous environment.
The study of even its monoform nonlinear oscillations is a difficult mathematical problem, because it is impossible to use existing analytical methods of partial differential equations to solve it.
To partially resolve the problem, an analytical and empirical method is developed. Its main idea includes the following tasks: based on empirical information on one of the forms of oscillations (usually of smaller amplitude), the approximate analytical relations of their amplitude–frequency characteristics are obtained; the obtained dependences are taken into account in the refined mathematical model of another form of oscillations. The latter non-autonomous type is much simpler to study than the “idealized model”. One of the most convenient analytical methods is used for this model, or a new approach is developed, for example, for systems with strong nonlinearity.

3. Related Works

The most effective methods of analytical study of nonlinear oscillatory processes of systems with concentrated masses and distributed parameters are those based on the main idea of perturbation methods [1,2]. They are effectively used to study single or multifrequency oscillations of elastic bodies, as well as systems with distributed parameters [3,4]. However, the operation of existing structures is usually accompanied by complex oscillations of their elements. For example, turbine rotors perform simultaneous bending and torsional oscillations [5,6], columns for drilling oil and gas wells perform both longitudinal bending and torsional oscillations [7,8], etc. Some experimental and theoretical studies show that components of vibrations of structural elements, in some cases interacting with each other, can lead to significant dynamic loads in them. These loads increase rapidly in magnitude in the so-called resonant cases, when the frequencies of the component oscillations (bending, longitudinal, torsional) or external and natural frequencies are related by rational relations [2,9] (i.e., oscillations in bodies with material properties that are or are nearly linearly elastic). To prevent them means to extend the service life of elastic bodies, i.e., structures as a whole, only on the basis of mathematical models adequate to the dynamic process and their solutions. If such phenomena under the action of external perturbation for both quasilinear and nonlinear systems with concentrated masses have found a relatively proper development (see, for example [1,9]), the issues of internal resonances for systems with distributed parameters were considered only in some works [9,10]. To prevent such phenomena, it is necessary to know in advance about the reaction of the object under study to a particular type of disturbance, as well as the result of the interaction of certain types of oscillations. However, the development of analytical methods for appropriate mathematical models that take into account a wide range of external and internal factors is a complex problem. The use of the same for the latter methods of numerical simulation in many cases does not give the desired results. The problem is much more complicated for elastic elements whose elastic properties of the material are described by nonlinear relations [11,12,13,14,15] or are characterized by longitudinal velocity [16,17,18,19,20,21,22]. The peculiarity of the dynamics of the last class of systems is an example of this fact. The velocity of the relative motion of the medium along the elastic body (or the body itself) affects the quantitative characteristics of the components of oscillations (longitudinal or bending), and thus the second component of oscillations (torsional). On the other hand, taking into account the constant velocity of longitudinal motion in the mathematical model of monoform oscillations of elastic bodies leads to a qualitatively new form for which it is not possible to apply known classical methods of integrating equations with partial derivatives [16,17,20,21]. Taken together, it requires the development of simplified methods for studying this class of systems and obtaining relatively simple analytical dependencies that could predict in advance such undesirable conditions for the resonant phenomena of the class of bodies (systems).
Thus, the paper attempts to partially solve such a little-covered problem in the literature as the development of methods for analytical research of complex nonlinear oscillations of the simplest elements of machines and mechanisms. Its main idea is illustrated by the example of nonlinear elastic bodies that carry both bending and longitudinal vibrations besides the body that moves along a continuous flow in a homogenous environment. Such systems include, in particular, columns for drilling wells along which a continuous flow of fluid moves, elements of pipelines for transporting liquids, working bodies of screw conveyors for moving bulk media [7,19,23], and others. They are the main focus of the work. Regarding the external perturbation of the oscillations of an elastic body, the cases of continuous and discrete action are considered in the work. The latter is analytically represented by delta functions [24,25,26,27,28,29,30,31]. The resonance condition for the discrete action of pulsed perturbation is obtained, and on the basis of generalization of the main idea of asymptotic methods [29,30,31] of nonlinear mechanics on systems with pulsed perturbation, the laws of change of basic parameters of resonant oscillations are described.

4. Materials and Methods

The main idea of the developed methodology for studying complex nonlinear elastic body oscillations is presented below on the example of the elastic body–CFHM system, which moves along a body with a constant velocity V . The elastic body under the action of external factors performs complex vibrations, known as bending and longitudinal vibrations. Partial information about one of the forms of vibrations, in this case longitudinal, which is known from experimental studies, is considered in experimental studies based on the processing of information obtained by the use of sensors. Regarding these specified data, the amplitude of these oscillations is much smaller than the bending. Despite this fact, the following vibrations under certain conditions can lead to a significant increase in the amplitude of another form of bending vibrations of the body and thus reduce the service life. The task is to:
(a)
develop methods for processing experimental information about longitudinal vibrations and bringing them to a form convenient for consideration when building a mathematical model of bending vibrations;
(b)
construct a mathematical model of the dynamics of the elastic body–CFHM system, which would take into account the above a priori information and a wide range of internal and external factors of continuous and discrete action;
(c)
develop methods for constructing the solution of the above mathematical model;
(d)
obtain analytical relations that describe the defining parameters of bending vibrations of an elastic body depending on external and internal factors of the system on the basis of the above method. That would enable calculations of specific elements of machines or systems at the design stage or during their operation.
Methods of processing experimental information about one of the forms of vibrations and bringing it to a form convenient for building a mathematical model of another form of vibrations. The peculiarity of the elastic body dynamics is the fact that its material satisfies the quasilinear elastic properties, and therefore the frequency of the dynamic process of such an elastic body depends on the amplitude. It is proposed to use a priori information about the amplitude and period of longitudinal vibrations of an elastic body to establish this relationship with respect to longitudinal vibrations [14]. Thus, it is believed that it is determined by a set of successive values of amplitudes and periods of vibrations. These characteristics of the longitudinal vibrations of a rigid body determine the following approximate ratios:
Δ b τ = d b d t ,   Δ τ τ = d τ d t
in which Δ b = b j + 1 b j , Δ τ = τ j + 1 τ j , Δ τ —length of j -scope, b j -magnitude.
Logarithms of amplitude change δ b = ln b j b j + 1 and period δ τ = ln τ j τ j + 1 of vibrations with a sufficient degree of accuracy can be approximated by polynomials
δ b = β 1 ( b j + 1 2 , ε ) ,   δ τ = [ ρ ¯ + β 2 ( b j + 1 2 , ε ) ] β 1 ( b j + 1 2 , ε )
In the above dependences, β 1 ( b j + 1 2 , ε ) and β 2 ( b j 1 2 , ε ) polynomials have the form of
β 1 ( b j + 1 2 , ε ) = n k ε n β 1 n k b j + 1 2 2 k
β 2 ( b j + 1 2 , ε ) = n k ε n β 2 n k b j + 1 2 2 k
It should be highlighted that in relations (3) and (4), b j + 1 2 = b j + b j + 1 2 , ρ ¯ is a constant, which is expressed in terms of physical and mechanical properties of the elastic body, the method of its fixation and for the quasilinear properties of the latter is zero, a small parameter indicating a slow change and time of amplitude and period of vibrations.
From the expressions for the logarithms of the change in the amplitude and period of vibrations on the basis of the mean theorem we obtain
τ j τ j + 1 d d t ( β 1 ( b j + 1 2 , ε ) ) d t = τ j ( d d t ln b ) | b = b j + 1 2
τ j τ j + 1 d d t ( β 2 ( b j + 1 2 , ε ) ) d t = τ j ( d d t ln τ ) | τ = j + 1 2 ,
2 π = ψ j + 1 ψ j = τ j ( d ψ d t ) | τ = τ j + 1 2
where ψ is phase of the longitudinal component of vibrations.
Regarding the obtained expressions to the characteristics of the vibrating process at a given time, we obtain
d b d t = n ε n b A ^ n ( b ) , d ψ d t = ϖ b ρ ¯ + n ε n B ^ n ( b ) ,
where A ^ n ( b ) , B ^ n ( b ) are known functions.
According to the above, the longitudinal movement of the cross section of the elastic body with the coordinate y for the first approximation can be represented as u ( y , t ) = b r ( ε t ) sin r π l y cos θ r ,   θ r = ( ϖ r t + ϕ r ) (l is length of the elastic body and the parameters r , ϖ are determined by the method of fixing the body, its physical and mechanical properties and cross-sectional shape).
Mathematical model of elastic body dynamics. As it follows from the formulation of the problem to be solved in the work, an arbitrary conditionally-isolated element of an elastic body, as well as SPOS, are in complex motion. To describe the kinematic parameters of CFHM, portable motion is considered to be the motion of an elastic body. The movement of the continuous flow of homogeneous medium is relative to the body and will be nothing more than a relative motion. As for the elastic body, which performs complex vibrations, bending and longitudinal ones, the first of them will be considered portable, then the longitudinal deformation of the thallus relative to the transversely deformed elastic body determines the relative motion. Assuming that in the undeformed position of the body its neutral axis occupies a position that coincides with the vertical axis OY, provided that the bending vibrations of the body occur in one plane (XOY plane), the portable transverse displacement of the elastic body is uniquely determined by the movement perpendicular to OY. Thus, the transverse vibrations are described by the function—the transverse movement of the neutral axis of the body with the coordinate at any time t; the relative displacement of the elastic body (longitudinal vibrations) is unambiguously described by the displacement of the specified cross section of the body along the deformed elastic axis. The specified value is determined, as mentioned above, by the function u ( y , t ) .
A mathematical model of the dynamics of the elastic body system—continuous flow of CFHM can be obtained based on the “dynamic equilibrium” of the conditionally selected element with the coordinate of its middle y and length dy.
In the case when the elastic body continuous flow of homogeneous medium performs one-form vibrations (bending) at hinged ends, its mathematical model has the following form [11]:
( m t + m s ) 2 w t 2 + 2 m s 2 w y t V + m s V 2 2 w y 2 + E I 4 w y 4 = ε f 1 ( w , w t , , 3 w y 3 , γ )
w ( y , t ) | y = j = 0 ,   2 w y 2 | y = j = 0 ,   j = 0 , l
in which l is the distance between the attachment points of the elastic body, m t ,   m s is accordingly the mass of the unit length of the elastic body and CFHM, E I is its bending stiffness ( E –E-modulus, moment of inertia of the body cross-section), f 1 ( w , w t , , 3 w y 3 , γ ) is a function that describes the resistance forces, the nonlinear component of the reducing force, the external 2 π -periodic γ perturbation with frequency μ ( γ = μ t ) and other forces whose maximum value is significantly less than the maximum value of the last term of the left side, as indicated by the small parameter ε .
As for the considered much more complex case, the mathematical model should additionally take into account the inertial forces of relative motion and the Coriolis force of an element of an elastic body and CFHM. In accordance with the above-accepted components of the movements of the system under consideration, the relative ( a r t ) and the Coriolis ( a c o r t ) acceleration of a point that coincides with the middle of an element of an elastic body (a point with a coordinate y ) is determined by the dependences
a r t = i 2 u t 2 sin α + j 2 u t 2 cos α ,
a c o r t = 2 ( 2 w y t ) k × ( i u t sin α 0 + j u t cos α 0 ) = 2 2 w y t { u t cos α 0 i + u t sin α 0 j }
where i , j , k unit vectors are directed along the axes OX, OY, and OZ (the OZ axis is directed perpendicular to the XOY plane); and α 0 is the angle that is formed by the tangent to the curved axis of the elastic body with the OY axis at the point with the coordinate y . For small vibrations of the system (such vibrations are considered in this paper), the condition for determining the specified angle is the correlation sin α 0 = α 0 = w y ,   cos α = 1 , and therefore in accordance with the specified 2 w y t is nothing more than the angular velocity of the transpositional motion of elastic body. As for the acceleration of CFHM, for its movement we have the following: the relative acceleration is equal to a e s = V 2 ρ п , where п is the unit vector of the internal normal to the curved axis of an elastic body, and that is why п = i cos α 0 + j sin α 0 , ρ -radius of curvature of the neutral axis of the body at a point with an ordinate y . Projections of the specified vector on the axis of the fixed system XOYZ are determined in accordance with the correlation
a e s = 2 w t 2 i + 2 u t 2 sin α 0 i + 2 u t 2 cos α 0 j
Similarly, the Coriolis acceleration vector continuous flow of homogeneous medium is determined by
a c o r s = 2 ( 2 w y t + 2 u y t ) k × V ( i sin α 0 + j cos α 0 ) = 2 V ( 2 w y t + 2 u y t ) ( i cos α 0 + j sin α 0 )
in which 2 u y t k is the vector of angular velocity of an elastic body in relative motion.
Taken together, the above allows us to write down the differential equation of bending vibrations of an elastic body, taking into account longitudinal vibrations and motion along CFHM in the form of
( m t + m s ) 2 w t 2 + E I 4 w y 4 + m s { 2 2 w y t V + V 2 2 w y 2 + 2 V ( 2 w y t + 2 u y t ) + 2 u t 2 w y ) } + m t { 2 u t 2 w y 2 2 w y t u t } = f 1 ( w , w t , , 3 w y 3 , γ ) .
Note 1:
(a) Equation (15) takes into account that for the case of small vibrations, it is true that sin α 0 = α 0 = w y ,   cos α = 1 , and the deplanation of flat sections is neglected;
(b) longitudinal vibrations do not change the curvature of the neutral axis of the elastic body;
(c) the differential relation of the “dynamic equilibrium” of the selected element in the projection on the OY axis can, together with Equation (15), serve as a basis for determining the dynamic response of anchorages.
Methodology for studying nonlinear bending vibrations of an elastic body, non-resonant case. As noted above, the main idea of the methodology is based on the use of a priori information for longitudinal vibrations, so the function u ( t , y ) and its partial derivatives are determined in accordance with the relation obtained on its basis u ( y , t ) = b r ( ε t ) sin r π l y cos θ r ) . Restrictions on the CFHM and magnitude of the amplitude of longitudinal vibrations allow us to represent differential Equation (15) in the form
( m t + m s ) 2 w t 2 + E I 4 w y 4 = ε f ( w , w t , , 3 w y 3 , γ ) + g ( w y , 2 w y 2 , 2 w y t , u t ) ,
where
ε f 1 ( w , w t , , 3 w y 3 , γ ) = ε f ( w , w t , , 3 w y 3 , γ ) m s ( 2 2 w y t V + V 2 2 w y 2 ) ,
g ( w y , 2 w y 2 , 2 w y t , u t , 2 u t 2 ) = m s { 2 V ( 2 w y t + 2 w y t ) + 2 u t 2 w y ) } m t { 2 u t 2 w y 2 2 w y t u t } expresses the effect on the bending vibrations of the longitudinal and movements of the CFHM.
It is obvious that the maximum value of the right-hand side of the obtained Equation (16) is a small value compared to the maximum value of the term E I 4 w y 4 . This is the basis for applying general ideas for constructing solutions to quasilinear equations with a small right-hand side [15,16,17,18] for the boundary value problem (16) and (10). Its peculiarity is that its right part is 2 π -periodic according to two arguments, θ (takes into account the influence of information known on the basis of experimental information, the influence of longitudinal vibrations on bending ones) and γ (takes into account the influence of external perturbation). The solution of Equation (17), more precisely its first approximation in the s -form of “dynamic equilibrium”, is presented in the form
w ( t , y ) = a s ( t ) ( sin ( χ s y + ψ s ( t ) ) + sin ( χ s y ψ s ( t ) ) ) + ε U 1 ( a s , y , ψ s , γ , θ r ) .
in which χ s = s π l , and U ( a s , x , ψ s , γ , θ r ) is an unknown 2 π -periodic function by arguments ψ s , γ , θ , which takes into account the influence of external and internal factors on changing the vibrations. As for the law of changes in its parameters a s ( t ) and ψ s ( t ) , their value significantly depends on the correlation between the frequency of natural bending vibrations ω s = ( s π l ) 2 E I m t + m s on the one hand and the frequency of the external periodic perturbation μ or longitudinal vibrations ϖ on the other hand. In the case when the correlations are justified, q ω s p μ тa and q ω s p ϖ r ( p ,   q are coprime numbers and this case will be called non-resonant) for bending vibrations of the function a s ( t ) and ψ s ( t ) in the asymptotic representation (18) are related by ordinary differential equations
d a s d t = ε A 1 ( a s ) ,   d ψ s d t = ω s ( χ s ) 2 m s m t + m s V 2 ω 1 + ε B 1 ( a s )
The right-hand sides of relations (19), more precisely unknown functions A 1 ( a s ) and B 1 ( a s ) , must be defined so that the asymptotic representation (18), taking into account (19), satisfies the original Equation (17) with the considered accuracy. For an unambiguous definition of the functions A 1 ( a s ) and B 1 ( a s ) on U ( a s , y , ψ s , γ , θ ) , we impose an additional condition—it does not contain the first mods in the argument scenario. This is equivalent to the following relation:
0 2 π U ( a s , y , ψ s , γ , θ r ) { sin ψ s cos ψ s } d ψ s = 0
In the physical interpretation, this condition is equivalent to the choice of the s mode amplitude for the bending vibrations amplitude. To find the indeterminate functions A 1 ( a s ) and B 1 ( a s ) , which describe the laws of variation of the basic parameters of bending vibrations of an elastic body by differentiating the asymptotic representation of the solution taking into account the relations (19), we obtain
w t = ε A 1 ( a s ) ( cos ( χ s y + ψ s ) cos ( χ s y + ψ s ) ) a s ( ω + ε B 1 ( a s ) ) ( sin ( χ s y + ψ s ) +   s i n ( χ s y + ψ s ) ) + ε ( ω s U 1 ( a s , x , ψ s , γ , θ r ) ψ s + U 1 ( a s , x , ψ s , γ , θ r ) θ r ϖ r + U 1 ( a s , x , ψ s , γ , θ r ) γ μ ) ,
2 U t 2 = a s ω s 2 ( sin ( χ s y + ψ s ) + sin ( χ s y + ψ s ) ) 2 ε ( A 1 ( a s ) + a s B 1 ( a s )   ω s ) ( cos ( χ s y + ψ s ) cos ( χ s y + ψ s ) ) + ε ( 2 U ψ s 2 ω 1 2 + 2 U θ r 2 ϖ r 2 + 2 U γ 2 μ 2 + 2 2 U ψ s θ r ω s ϖ r + 2 2 U ψ s γ ω s μ + 2 2 U γ θ r μ ϖ r ) ,
2 u y 2 = a s χ s 2 ( cos ( χ s y + ψ s ) cos ( χ s y + ψ s ) ) + ε 2 U 1 y 2 u y = a s χ s ( sin ( χ s y + ψ s ) sin ( χ s y ψ s ) ) + ε U 1 y 3 u y 3 = a s χ s 3 ( sin ( χ s y + ψ s ) + sin ( χ s y + ψ s ) ) + ε 3 U 1 y 3 , 4 u y 4 = a χ s 4 ( cos ( χ s y + ψ s ) cos ( χ s y + ψ s ) ) + ε 4 U 1 y 4 .
Thus, we obtain the differential equation for the connection between the indeterminate functions U 1 ( a 1 , y , ψ 1 , γ , θ ) and A 1 ( a 1 ) and B 1 ( a 1 ) for the first mode of vibrations ( s = 1 ).
( ω 1 ) 2 2 U 1 ψ 1 2 + ϖ r 2 2 U 1 θ r 2 + μ 2 2 U 1 γ 2 + 2 ω 1 μ 2 U 1 ψ θ r + 2 ω 1 ϖ r 2 U 1 ψ θ r + 2 μ ϖ r 2 U 1 γ θ r + E I m t + m s 4 U 1 x 4 = 1 m s + m t F 1 ( a 1 , y , ψ 1 , γ , θ r ) + G 1 ( a 1 , y , ψ 1 , θ r ) + 4 A ( a 1 ) ω 1 sin χ 1 y sin ψ 1 + 4 B ( a 1 ) a 1 ω 1 sin χ 1 y cos ψ 1
where F ¯ 1 ( a s , y , ψ s , γ ) corresponds to the value f ( w , w t , , 3 w y 3 , γ ) m s ( 2 2 w y t V + V 2 2 w y 2 ) if the function w ( y , t ) in it and its derivatives is determined according to the main part in representation (5). As for the function G 1 ( a s , y , ψ s , θ r ) , it, in accordance with the agreed notation, takes approximate values for the considered value as follows. It allows one to obtain the relations to determine the required functions
G 1 ( a 1 , y , ψ 1 , θ r ) 1 m s + m t { 2 m s V ( b r ϖ r r π l cos r π l y sin θ r ) 2 a 1 m t b r χ 1 ϖ r 2 sin r π l y cos θ r cos χ s y cos ψ s 2 a 1 [ m t b r ϖ r sin r π l y sin θ r ] cos χ 1 y sin ψ 1 } d a 1 d t = ε 6 π 2 ( m t + m s ) l ω 1 0 l 0 2 π 0 2 π F 1 ( a 1 , y , ψ 1 , γ ) sin χ 1 y sin ψ 1 d y d ψ 1 d γ + ε 6 π 2 ( m t + m s ) l ω 1 0 l 0 2 π 0 2 π G 1 ( a 1 , y , ψ 1 , θ r ) sin χ 1 y sin ψ 1 d y d ψ 1 d θ r , d ψ 1 d t = ω 1 ε 6 π 2 ( m t + m s ) l a 1 ω 1 0 l 0 2 π 0 2 π F 1 ( a 1 , y , ψ 1 , γ ) sin χ 1 y cos ψ 1 d y d ψ 1 d γ + + ε 6 π 2 ( m t + m s ) l a 1 ω 1 0 l 0 2 π 0 2 π G 1 ( a 1 , y , ψ 1 , θ r ) sin χ 1 y cos ψ 1 d y d ψ 1 d θ r .
From the form of the function G 1 ( a 1 , y , ψ 1 , θ r ) , it can be seen that
0 l 0 2 π 0 2 π { cos ψ s sin ψ s } G 1 ( a s , y , ψ s , θ r ) sin χ s y d y d ψ s d θ r 0
The obtained relation allows us to draw the following conclusion: for the first non-resonance approximation, the longitudinal vibrations of the elastic body do not affect the law of change of the single-frequency bending vibrations amplitude. As for the CFHM, it follows from the specified properties of the system that describe the shape of the defected axis of the elastic body, the frequency and amplitude of vibrations. Taken into consideration all the above, the system of differential equations for the first non-resonance approximation (25) can be transformed into
d a 1 d t = ε 8 π 2 ( m t + m s ) l ω 1 s 0 2 π 0 2 π 0 l F 1 ( a 1 , y , ψ 1 , γ ) sin χ 1 y sin ψ 1 d y d ψ 1 = ε A 1 ( a 1 ) d ψ 1 d t = ω 1 m s ( m t + m s ) ω 1 V 2 ( χ 1 ) 2 ε 8 π 2 ( m t + m s ) l a ω 1 0 2 π 0 2 π 0 l F 1 ( a 1 , y , ψ 1 , γ ) sin χ 1 y cos ψ 1 d y d ψ 1 d γ = ω 1 m s ( m t + m s ) ω s V 2 ( χ 1 ) 2 + ε B 1 ( a 1 ) .
For the case when the material of the elastic body satisfies the nonlinear technical law of elasticity [6], and the resistance force is proportional to the speed in the degree s 1 , the right part of Equation (19) takes the form for the case of the main mode of vibrations
A 1 ( a ) = υ a s 1 , ε B 1 ( a ) = 3 k 2 ε 32 π 2 l 2 a 2 ω 1 .
According to the relations (19), taking into account (28) in Figure 1, the laws of change of the natural bending vibrations frequency from the amplitude and speed of the CFHM relative motion (Figure 1a) and the CFHM bulk mass and the length of the elastic element (Figure 1b) at the constant speed of the CFHM relative motion at such parameters of the system l = 8   m , m t = 30   kg / m ,   m s = 20   kg / m ,   E = 2 × 10 12   H / m 2 , I = 6 × 10 6   m 4 are presented for different values of the CFHM bulk mass.
Thus, the influence of nonlinear elastic characteristics of the elastic body material and the motion of a continuous flow of medium along it is manifested as a function of its natural frequency of vibrations on the amplitude and speed of the latter. In particular, the natural frequency of nonlinear vibrations of an elastic body:
for larger values of the relative flow rate of the continuous medium is less;
decreases with increase of elastic body length;
is smaller for a continuous flow of medium of smaller specific gravity.
The above facts play a significant role in the study of resonance vibrations of the system of elastic body–CFHM continuous flow. These vibrations are the subject of further research.
Methods for studying nonlinear bending vibrations of an elastic body, resonance case. Resonance vibrations almost always play a negative role for mechanical systems, because during resonance they significantly increase the amplitude of vibrations, which means that dynamical loads increase. This reduces the service life of the systems. The peculiarity of the resonance vibrations of the considered elastic body is the fact that they can be caused not only by external forces but also by longitudinal vibrations; in addition, they are affected by the CFHM motion.
In order to find their influence on the resonance values of the amplitude, as in the non-resonance case, we will look for the solution of the problem in the asymptotic representation (18) with the difference that in the resonance case, the pipe vibration’s amplitude significantly depends on the phase difference of natural and forced vibrations. First of all, consider the internal resonance. As noted above, it occurs due to longitudinal vibrations, the main parameters of which are known from the processing of experimental studies, and they can be interpreted as external (in relation to bending) action. The phenomenon of resonance in nonlinear systems depends significantly on the phase difference of their natural and coercive forces. For this purpose, we introduce the specified parameter according to the relation ϕ = ψ 1 θ r , and therefore in the asymptotic representation (18), the amplitude parameter is determined by a more complex relation, namely
d a d t = ε A 1 ( a , ϕ ) ,   d ϕ d t = ω 1 ( χ 1 ) 2 m s m t + m s V 2 ω 1 q p ϖ r + ε B 1 ( a , ϕ ) .
Here and below, the index in the amplitude of the main mode of vibrations is omitted for simplicity. The indeterminate functions A 1 ( a , ϕ ) , B 1 ( a , ϕ ) are located in such a way that the asymptotic representation (18), taking into account (28), satisfies the original equation with the considered degree of accuracy. To do this, by differentiating (18), taking into account the above, we have
u t = ε   A 1 ( a , ϕ ) ( cos ( χ 1 y + ϕ + θ r ) cos ( χ 1 y ( ϕ + θ r ) ) ) a ( ω + μ   B 1 ( a , ϕ ) ) ( sin ( χ 1 y + ϕ + θ r ) + sin ( χ 1 y ( ϕ + θ r ) ) ) + q p ϖ r ε U 1 θ r + ω 1 ε U 1 ψ 1 2 u t 2 = ε A ( a , ϕ ) ϕ ( ω 1 q p ϖ r ) ( cos ( χ 1 y + ϕ + θ r ) cos ( χ 1 y ( ϕ + θ r ) ) ) ε ω A 1 ( a , ϕ ) ( sin ( χ 1 y + ϕ + θ r ) + sin ( χ 1 y ( ϕ + θ ) ) ) a ( ω 2 + 2 ε ω B 1 ( a , χ 1 y + ϕ + θ r ) ) ( cos ( χ 1 y + ϕ + θ r ) cos ( χ 1 y ( ϕ + θ r ) ) ) ε a B 1 ( a , ϕ ) ϕ ( ω q p ϖ r ) ( sin ( χ 1 y + ϕ + θ r ) + sin ( ( χ 1 y ( ϕ + θ r ) ) ) + ε ( 2 U 1 ψ 1 2 ω 1 2 + 2 U 1 θ 2 ( q p ϖ r ) 2 + 2 U 1 γ 2 μ 2 + 2 2 U 1 θ ψ 1 q p ϖ r ω 1 + 2 2 U 1 θ γ q p ϖ r μ + 2 2 U 1 γ ψ 1 μ ω 1 )
The obtained relations allow us to write a differential equation, which connects the required functions in the form
ω 1 2 2 U 1 ψ 1 2 + ( q p ϖ r ) 2 2 U 1 θ r 2 + μ 2 2 U 1 γ 2 + 2 q p ω 1 ϖ r 2 U 1 ψ 1 θ r + 2 ω 1 μ 2 U 1 ψ γ + 2 q p μ ϖ r 2 U 1 γ θ r + E I m t + m s 4 U 1 y 4 =   F ( y , a , ψ , γ ) + a m s m t + m s V 2 ( π l ) 2 sin χ π l x cos ψ 1 +   G ( y , a , ψ , θ r ) + sin π y l × ( cos ψ ( A ( a , ϕ ) ϕ ( ω 1 q p ϖ r ) + 2 a ω 1 B ( a , ϕ ) ) + sin ψ ( a B ( a , ϕ ) ϕ ( ω q p ϖ r ) + 2 A ( a , ϕ ) ω ) ) .
An indeterminate function will satisfy homogeneous boundary conditions if it can be represented as a series
U 1 ( y , a , ψ , γ , θ , ) = m sin m π l y U 1 m ( a , θ , ψ 1 , γ ) .
In this case, the functions A 1 ( a , ϕ ) , B 1 ( a , ϕ ) and their derivatives are connected by a relation
a V 2 m s m t + m s ( π 2 2 l ) cos ψ 1 + 1 ( m t + m s ) l ω s 0 l F ( a , y , ψ 1 , θ ) sin π l y d y + 1 l 0 l G ( a , y , ψ 1 , θ ) sin π l y d y        + ( cos ψ ( A ( a , ϕ ) ϕ ( ω q p ϖ r ) + 2 a ω B ( a , ϕ ) )        + sin ψ ( a B ( a , ϕ ) ϕ ( ω q p ϖ r ) + 2 A ( a , ϕ ) ω ) ) = 0
Relation (33) is the basis for determining the right-hand members of differential Equation (29), which describe the change in the basic parameters of bending vibrations of an elastic body in the resonant zone. Indeed, if U 1 ( a , x , ψ 1 , γ , θ r ) is 2 π periodically with an argument ψ 1 , and does not contain in the splitting of its first modes of the splitting, then the coefficients of its splitting for the specified mode of the function U 1 m ( a , θ , ψ 1 , γ ) have the same properties. Thus, we obtain a system of equations from relation (33).
a B ( a , ϕ ) ϕ ( ω q p ϖ r ) + 2 A ( a , ϕ ) ω = f ¯ 1 ( a ) 1 π 0 l 0 2 π G ( a , y , ϕ + θ r , θ r ) sin π l y sin ( ϕ + θ r ) d θ r d y , A ( a , ϕ ) ϕ ( ω q p ϖ r ) 2 a ω B ( a , ϕ ) = f ˜ 1 ( a ) + a V 2 m s m t + m s ( π 2 2 l ) 1 π l 0 l 0 2 π G ( a , y , ϕ + θ r , θ r ) sin π l y cos ( ϕ + θ r ) d θ r d y ,
where
f ¯ 1 ( a ) = 1 4 π 2 ( m t + m s ) l 0 l 0 2 π 0 2 π F ( a , y , ψ 1 , θ r ) sin π l y sin ψ d ψ 1 d θ r d y , f ˜ 1 ( a ) = 1 4 π 2 ( m t + m s ) l 0 l 0 2 π 0 2 π F ( a , y , ψ 1 , θ r ) sin π l y cos d ψ 1 d θ r d y .
The system of differential Equation (34) allows us to represent equation of the first approximation of resonance vibrations caused by the CFHM motion and longitudinal vibrations in the form
d a d t = ε f ¯ 1 ( a ) + 1 π ω 1 1 π 0 l 0 2 π G ( a , y , ϕ + θ r , θ r ) sin π l y sin ( ϕ + θ r ) d θ r d y d ϕ d t = ω 1 q p ϖ r a V 2 m s m t + m s ( π 2 2 l ) + ε f ˜ 1 ( a ) + 1 π a ω 1 1 π 0 l 0 2 π G ( a , y , ϕ + θ r , θ r ) sin π l y cos ( ϕ + θ r ) d θ r d y
here
f ¯ 1 ( a ) = 1 4 π 2 ( m t + m s ) l 0 l 0 2 π 0 2 π F ( a , y , ψ 1 , θ r ) sin π l y sin ψ d ψ 1 d θ r d y f ˜ 1 ( a ) = 1 4 π 2 ( m t + m s ) l 0 l 0 2 π 0 2 π F ( a , y , ψ 1 , θ r ) sin π l y cos d ψ 1 d θ r d y .
Differential equations describe the change in the amplitude of transverse vibrations of an elastic body during the transition through resonance at the main frequency of longitudinal vibrations at different speeds of motion of CFHM, provided that the nonlinear elastic properties of the body material satisfy the nonlinear technical law of elasticity, and the resistance force is proportional to the speed of portable motion in the power of s 1 :
d a d t = k ¯ 1 ( a ) s 1 b 4 ( m t + m s ) m s V l ϖ κ π 2 l sin ϕ d ϕ d t = ( π l ) 2 E I m t + m s ϖ 2 m s V 2 m t + m s ( π l ) 2 3 k ¯ 2 32 a 2 + b 4 ( m t + m s ) a m s V l ϖ κ π 2 l cos ϕ
where a is the amplitude of the main internal resonant vibrations.
As for the subresonances that occur when the condition is met κ = 2 s , they are described by dependencies.
d a d t = 1 4 π ( m t + m s ) l ω s { ε 0 2 π 0 2 π 0 l F 1 ( a , y , ψ s , γ ) sin χ s y sin ψ s d y γ d ψ s + + 8 l a m s b χ s ϖ 2 m t + m s p q 4 p 2 q 2 cos p π q [ c o s p π q sin 2 p π q 2 c o s p π q sin 2 p π q cos 2 ϕ + + 1 2 s i n p π q ( 1 4 cos 2 p π q ) s i n 2 ϕ ] a 4 p q l 4 p 2 q 2 m t b ϖ cos p π q [ c o s p π q ( 2 sin 2 ϕ cos 2 p π q 2 sin 2 ϕ + 4 cos 2 ϕ cos p π q sin p π q sin 2 p π q ) s i n p π q ( 2 cos 2 ϕ 1 ) ] } d ϕ d t = 2 ω s p q ϖ + 1 4 π ( m t + m s ) l a ω s { ε 0 l 0 2 π 0 2 π F 1 ( a s r i s , y , ψ s , γ ) sin χ s y cos ψ s d y d ψ s d γ 2 m s V 2 ( χ s ) 2 d ϕ d t 2 l a m s b χ s ϖ 2 m t + m s 4 p q 4 p 2 q 2 cos m π n [ sin p π q ( 1 2 cos 2 p π q ) 2 cos p π q sin 2 p π q sin 2 ϕ 2 sin p π q ( 1 2 cos 2 p π q ) cos 2 ϕ ] a 4 p q l 4 p 2 q 2 m t b ϖ cos p π q [ cos p π q sin 2 p π q ( 1 2 cos 2 ϕ ) 1 2 sin p π q ( 1 4 cos 2 p π q ) sin 2 ϕ ] }
where ϕ = 2 ψ p q θ r .
In Figure 2 is presented, in accordance with relations (34), the change in the amplitude of bending vibrations caused by longitudinal vibrations for the case when the material of elastic body satisfies the nonlinear technical law of elasticity, and the resistance force is proportional to the velocity for s 1 = 3 for different values of the relative amount of motion.
In Figure 3a, a change in the amplitude of bending vibrations caused by longitudinal vibrations at different values of amplitudes can be observed, while in Figure 3b is shown the motion of CFHM at different values of the speed of its movement.
By analyzing the analytical dependences obtained above, it is established that internal resonances for an elastic body that performs bending and longitudinal vibrations can exist on even and odd modes of bending vibrations, and on even modes they are caused by longitudinal vibrations, and on odd modes by the motion CFHM. As for the quantitative characteristics of the internal resonance, the given graphical dependences show that the amplitude of the transition through the internal resonance is less:
at higher frequencies of longitudinal vibrations than at the main one;
for large values of the relative amount of motion CFHM;
with a faster transition through resonance.
The above represented results have a practical component; in order to reduce dynamic loads in an elastic body that performs longitudinal and bending vibrations and along which a continuous stream motion of CFHM, it is necessary to move the latter at the highest technologically possible speed. At the same time, a special case mentioned above, m s = 0 , is the result concerning resonant bending vibrations caused only by longitudinal vibrations. They, as follows from relations (34), exist only on paired modes of bending vibrations.
The main idea of studying internal resonant phenomena in the elastic body–CFHM system described above has been transferred to the study of external resonant phenomena. Below, we give only the basic relations concerning the continuous action of an external periodic perturbation, and the case of a discrete periodic (pulsed) action on the system under consideration is considered in more detail, because it requires some additional mathematical calculations.
Therefore, the condition for the main external resonance under the continuous action of an external perturbation is ω 1 μ The relations that describe the resonant process at the main frequency of natural vibrations take the form
a B ¯ ( a , ϕ ¯ ) ϕ ( ω 1 μ ) + 2 A ¯ ( a , ϕ ¯ ) ω 1 = 1 2 ( m t + m s ) l π 0 l 0 2 ϕ F ( a , y , ϕ ¯ + θ , θ ) F sin π l y sin ( ϕ ¯ + θ ) d y d θ + 1 2 ( m t + m s ) l π 0 l 0 2 ϕ F ( a , y , ϕ ¯ + θ , θ ) F sin π l y cos ( ϕ ¯ + θ ) d y d θ ,
where ϕ ¯ is the relations that describe the resonant process at the main frequency of natural vibrations that take the form of the phase difference between natural and external forced vibrations, i.e., ϕ ¯ = ψ θ , so functions A ¯ ( a , ϕ ¯ ) and B ¯ ( a , ϕ ¯ ) describe the laws of change in the resonant parameters of the system under study using relations similar to (29), with the only difference that the resonant process is considered at the frequency of external perturbation.
In the specific case considered for non-resonant vibrations under the condition of external periodic perturbation, the system of Equation (35) is transformed to the form
a B ¯ ( a , ϕ ¯ ) ϕ ( ω 1 μ ) + 2 A ¯ ( a , ϕ ¯ ) ω 1 = υ a s 1 + H 2 ( m t + m s ) ω l π cos ϕ ¯ , A ¯ ( a , ϕ ¯ ) ϕ ( ω μ ) 2 a ω   B ¯ ( a , ϕ ¯ ) = a V 2 m s m t + m s ( π 2 l ) 3 ε 32 π 2 k 2 ¯ l 2 a 2 ω 1 H 2 ( m t + m s ) a ω l π sin ϕ ¯ ,
where H is the value of the external sinusoidal periodic perturbation.
The differential dependences obtained above also allow us to obtain relations for determining the amplitude of stationary resonant vibrations
β m + m 1 ( ω ¯ ) s 1 1 a s 1 2 ε H ¯ π ( ω ¯ + ν ( t ) ) cos ϕ = 0 , ω ¯ μ ( t ) ( π l ) 2 m s m t + m s V 2 8 ω = ε ( 2 H ¯ π ( ω ¯ + μ ( t ) ) a sin ϕ + 3 k ¯ 2 32 π 2 l 2 a 2 ω )
By simple transformations from dependencies (37), we obtain a resonant curve that describes the dependence of the amplitude of stationary resonant vibrations on the unbalance of the frequencies of natural and forced vibrations and the parameters of the system
( 2 ε H ¯ π ( ω + μ ( t ) ) ) 2 = ( β m t + m s ( ω ¯ ) s 1 1 a s 1 ) 2 + a 2 ( ω ¯ μ ( t ) ( π l ) 2 m s m t + m s V 2 8 ω ¯ 3 k ¯ 2 ε 32 π 2 l 2 a 2 ω ) 2
Based on ratio (38), Figure 4 presents the dependencies of the amplitude of resonant stationary vibrations of an elastic body on the unbalance α = ω μ of the frequencies of natural and forced vibrations at different speeds of motion of CFHM.
It follows from the above that when the sign of the nonlinear component of the reducing force changes, the slope of the resonant curves changes. Comparing the graphical dependencies for different models of the resistance force, it can be argued that an increase in the degree of its nonlinearity leads to a decrease in the resonant value of the vibration amplitude and a narrowing of the resonant zone.
As for the unknown function U 1 ( a 1 , y , ψ 1 , γ , θ ) , and for non-resonant or resonant cases, it partially affects only the change in the shape of bending vibrations and can be found without much difficulty by decomposing into multiple Fourier series, followed by equating the coefficients in the right and left parts of the relations that follow from (24) and (25), for the non-resonant case, and (31) and (36) for the resonant one. The case of pulsed (discrete) action of an external disturbance is much more important and at the same time much more difficult to develop a research methodology. Below, we consider the case when the specified perturbation acts only at a fixed point of the elastic body and is repeated at regular intervals. As for the limitations concerning the specified external perturbation of the vibrations, the elastic body–continuous flow of CFHM, as for continuous periodic perturbation, the maximum value of the pulse perturbation is small compared to the maximum value of the reduced transverse force. Thus, the differential equation of the bending component of the elastic oscillations takes the following form:
( m t + m s ) 2 w t 2 + E I 4 w y 4      = ε { f 1 ( w , w t , , 3 w y 3 ) + δ ( y y 0 ) j = 1 f ^ j ( w , w t , , 3 w y 3 ) δ ( t j T ) }      + g ( w y , 2 w y 2 , 2 w y t , u t ) ,
in which y 0 is the coordinate of the elastic body in which the periodic impulse force acts with the period T , and δ ( .... ) is a function of the corresponding argument. Regarding the following forms, f 1 ( w , w t , , 3 w y 3 ) and f ^ j ( w , w t , , 3 w y 3 ) , we can consider them to be analytical. The form g ( w y , 2 w y 2 , 2 w y t , u t ) takes into account the effect of longitudinal oscillations and CFHM on bending.
From the properties of delta functions [19,20,21], it follows that the maximum values of the periodic action of pulsed perturbation at time points j T are equal f ^ j ( w ( t , y ) , w ( t , y ) t , , 3 w ( t , y ) y 3 ) | t = j T , y = y 0 . To generalize the above method of studying the effect of longitudinal vibrations and dynamics of CFHM on bending, primarily on the basis of the properties of completeness and orthonormal system of functions that describe the forms of natural bending vibrations of an elastic body, delta function from a linear variable δ ( y y 0 ) will be as follows: δ ( y y 0 ) = j = 1 c m sin m π l y , in which the modulus c m is as follows c m = 2 l 0 l δ ( y y 0 ) sin m π l y d y = 2 l sin m π l y 0 . As for the time component of the pulse perturbation, the most interesting case of the considered drilling is the case when the perturbation has a period close to the period of its own bending oscillations. Hence, in the form presented below we will consider T 2 π ω 1 pulsed perturbations causing resonant bending vibrations. In this case, the temporal component of its action can be precisely represented [20] in the following form: δ ( t j T ) = cos ( ω 1 t ) δ ( t j 2 π ω 1 ) = cos θ δ ( θ θ j ω 1 ) , θ j = 2 π j . For simplicity, it is assumed that the first impulse action coincides with the initial moment of motion. Thus, the conducted partial sampling allows to present the initial mathematical model in a simpler form:
( m t + m s ) 2 w t 2 + E I 4 w y 4 = g ( w y , 2 w y 2 , 2 w y t , u t ) + ε { f 1 ( w , w t , , 3 w y 3 ) + 2 l m sin m π l y 0 sin m π l y j = 1 f ^ j ( w , w t , , 3 w y 3 ) cos θ δ ( θ θ j ω 1 ) } .
Note that the procedures performed with delta functions do not violate the accuracy of the obtained mathematical model in the form of relation (40) compared to the original (differential Equation (39)); in addition, to operate with the right part of the latter equation is much easier than the original. The use of the developed methodology for the study of resonant oscillations due to external periodic action for the boundary value problem (40), (10) allows the relationship of the main resonance to be represented as follows:
a B ¯ ( a , ϕ ¯ ) ϕ ( ω 1 2 π T ) + 2 A ¯ ( a , ϕ ¯ ) ω 1 = f ¯ 1 ( a ) + 1 2 ( m t + m s ) l π m sin m π l y 0 j 0 l 0 2 ϕ sin m π l y F j ( a , y , ϕ ¯ + θ ) cos θ δ ( θ θ j ω 1 ) sin π l y sin ( ϕ ¯ + θ ) d y d θ    A ¯ ( a , ϕ ¯ ) ϕ ( ω 1 2 π T ) 2 a ω 1 B ¯ ( a , ϕ ¯ )      = f ˜ 1 ( a ) + a V 2 m s m t + m s ( π 2 l )      + 1 2 ( m t + m s ) l π m sin m π l y 0 j 0 l 0 2 ϕ sin m π l y F j ( a , y , ϕ ¯      + θ ) cos θ δ ( θ θ j ω 1 ) sin π l y cos ( ϕ ¯ + θ ) d y d θ
The system of differential equations is the basis for the study in the general case of the influence of impulse perturbation on the bending vibrations of the system continuous flow of CFHM. It is greatly simplified in the case where the impulse perturbation is constant, i.e., f j ( w , w t , , 3 w y 3 ) F 0 = c o n s t . For this case, the basic equations for the main resonance in the case when the elastic properties of the body are described by the nonlinear technical law of elasticity, and the force of resistance is proportional to the velocity to the degree they take the form
a B ¯ ( a , ϕ ¯ ) ϕ ( ω 1 2 π T ) + 2 A ¯ ( a , ϕ ¯ ) ω 1 = υ a s 1 + T F 0 ( m t + m s ) ω l sin π y 0 l j = 1 cos j T cos ϕ ¯ ¯ , A ¯ ( a , ϕ ¯ ) ϕ ( ω 1 2 π T ) 2 a ω B ¯ ( a , ϕ ¯ ) = a V 2 m s m t + m s ( π 2 l ) 3 k 2 ε 32 π 2 l 2 a 2 ω 1 T F 0 2 ( m t + m s ) a ω l π sin π y 0 l j = 1 cos j T sin ϕ ¯ ¯ ,
For the first resonant approximation, the system of differential Equation (42) describes the main parameters of resonant vibrations under the action of external periodic pulse perturbation by the following dependence:
d a d t = = υ a s 1 + T F 0 ( m t + m s ) ω l sin π y 0 l j = 1 cos j T cos j T cos ϕ ¯ ¯ , d ϕ ¯ d t = ω 1 2 π T a V 2 m s m t + m s ( π 2 l ) 3 k 2 ε 32 π 2 l 2 a 2 ω 1 T F 0 2 ( m t + m s ) a ω l π sin π y 0 l j = 1 cos j T sin ϕ ¯ ¯ .
According to the above dependence, Figure 5 shows the change in the amplitude of bending vibrations during the transition through resonance at different lengths of the elastic body, the quantities of the relative motion of the CFHM, and the points of application of impulse perturbation.
The presented graphical dependencies show the amplitude of the passage through the main resonance of the bending vibrations of the elastic body, which is due to the periodic pulse action. It becomes more important in the case when the point of application of this action is closer to the middle of the elastic body. As for the continuous flow of a homogeneous medium, for larger values of the relative amount of motion, it not only reduces the natural frequency of the elastic body but also slightly reduces the maximum value of the amplitudes of the system through the resonance.
The above results can be generalized without much difficulty in the case of the action of impulse perturbation at several points of an elastic body.

5. Conclusions

A methodological study of complex vibrations has been developed for the mechanical system of an elastic body–continuous flow of homogeneous medium (CFHM). The methodology is based on:
(a)
the use of a priori information on some of the vibrations forms;
(b)
the construction of a “simplified” mathematical model of system dynamics, which provides a wide range of external and internal factors and information on the longitudinal vibrations of an elastic body;
(c)
obtaining analytical relationship that describe the overall range of factors on the elastic vibrations of system.
Considering the above given information, it allows us to conclude the following:
  • CFHM motion affects both the natural frequency of bending and longitudinal vibrations;
  • the most dangerous resonance vibrations of the specified systems type can be caused both by external periodic forces and mutual influence of one form of vibrations on the others (internal resonances);
  • internal resonances for an elastic body that performs bending and longitudinal vibrations can exist in even and odd modes of bending vibrations, and they are caused by longitudinal vibrations in even modes, and by continuous flow of homogeneous medium motion in odd modes;
  • the amplitude of the bending vibrations of the transition through the internal resonance at higher longitudinal vibrations frequencies is less than at the second frequency;
  • the amplitude of the transition through the main resonance for larger values of the amount of CFHM motion is smaller than the amplitude of the transition during the rapid transition through the resonance.
The methodology can be used for other cases of complex vibrations of elastic bodies, and the obtained results can serve as a basis for choosing the main technological and operational parameters of elastic elements of mechanisms and machines that perform complex vibrations.
Practical use of the obtained results. The theoretical results obtained in this work can be the basis for calculating the following:
  • the blast effect on the elements of protective structures in order to increase their protective capacity by improving the method of their attachment or by using additional reinforcement;
  • buff load effects on the elements of drilling strings and dynamic processes that occur during surface strengthening by work hardening in order to avoid resonance phenomena;
  • technological processes of vibration displacement or vibration separation of granular media;
  • longitudinal frequencies of movement of rope lifts in order to avoid resonant phenomena in them, and thus to extend the resource to make the object;
  • reducing dynamic performance bearings for mounting rotors;
  • assessing structural strength and vibration reduction of mechanical systems in energy, taking into account the optimal parameters (motion, mass, stiffness, etc.), which are found on the basis of empirical information and methodology developed in the work using interpretation and synthesis.

Author Contributions

Study analysis, A.A., S.F., and M.S.; methodology, A.A., S.F., and M.S.; resources, S.F., A.A., and M.S.; validation, S.F., A.A., and M.S.; writing—original draft preparation, Y.S., S.F., A.A., Y.K., and M.S.; writing—review and editing, Y.S., S.F., A.A., T.P. and M.S.; visualization, S.F., A.A. and M.S.; project administration, S.F., T.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Slovak Research and Development Agency under the contract No. APVV-19-0581 and by of the Faculty of Management of Comenius University in Bratislava, Slovak Republic.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bogolyubov, N.N.; Mitropolsky, Y.A. Asymptotic Methods in the Theory of Nonlinear Oscillations; Gordon and Breach Publishers: London, UK, 1958; p. 408. [Google Scholar]
  2. Mitropolskii, Y.A.; Moseenkov, B.I. Asymptotic Solutions of Partial Differential Equations; Vyshcha Shkola: Kyiv, Ukraine, 1976. [Google Scholar]
  3. Pukach, P.Y.; Kuzio, I.V. Resonance phenomena in quasi-zero stiffness vibration isolation systems. Sci. Bull. Natl. Min. Univ. 2015, 3, 62–67. [Google Scholar]
  4. Yarka, U.; Fedushko, S.; Veselý, P. The Dirichlet Problem for the Perturbed Elliptic Equation. Mathematics 2020, 8, 2108. [Google Scholar] [CrossRef]
  5. Goldin, A.S. Vibration of Rotary Machines; Mechanical Engineering: Moscow, Russia, 2000; p. 344. [Google Scholar]
  6. Peleschyshyn, A.; Syerov, Y.; Tarasov, D.; Zhezhnych, P. Intelligent system of information gathering and analysis from news Internet-resources. In Proceedings of the Advanced Computer Systems and Networks: Design and Application (ACSN-2005), Lviv, Ukraine, 21–23 September 2005; pp. 48–52. [Google Scholar]
  7. Ohorodnikov, P. Investigation of the relationship between longitudinal and torsional vibrations of the drill string. Oil Ind. Ukr. 2014, 2, 6–9. [Google Scholar]
  8. Huliaiev, V.I.; Borshch, O.I. Spiral waves in swirling elastic quills rotating with internal fluid flow. Acoust. Bull. 2007, 10, 12–18. [Google Scholar]
  9. Andrukhiv, A.; Sokil, B.; Sokil, M. Resonant phenomena of elastic bodies that perform bending and torsion vibrations. Ukr. J. Mech. Eng. Mater. Sci. 2018, 3, 65–73. [Google Scholar] [CrossRef]
  10. Andrianov, I.V.; Danishevskyi, V.V.; Ivankov, A.O. Asymptotic Methods in the Theory of Vibrations of Beams and Plates; Dnipropetrovsk State Academy of Civil Engineering and Architecture: Dnipro, Ukraine, 2010; p. 217. [Google Scholar]
  11. Gendelman, O.; Vakakis, A.F. Transitions from localization to nonlocalization in strongly nonlinear damped oscillators. Chaos Solitons Fractals 2000, 11, 1535–1542. [Google Scholar] [CrossRef]
  12. Cveticanin, L. Period of vibration of axially vibrating truly nonlinear rod. J. Sound Vib. 2016, 374, 199–210. [Google Scholar] [CrossRef]
  13. Cveticanin, L.; Pogany, T. Oscillator with a sum of non-integer order non-linearities. J. Appl. Math. 2012, 649050. [Google Scholar] [CrossRef]
  14. Pukach, P.Y.; Kuzio, I.V.; Nytrebych, Z.M.; Ilkiv, V.S. Analytical methods for determining the effect of the dynamic process on the nonlinear flexural vibrations and the strength of compressed shaft. Sci. Bull. Natl. Min. Univ. 2017, 5, 69–76. [Google Scholar]
  15. Pukach, P.Y.; Kuzio, I.V.; Nytrebych, Z.M.; Ilkiv, V.S. Asymptotic method for investigating resonant regimes of non-linear bending vibrations of elastic shaft. Sci. Bull. Natl. Min. Univ. 2018, 1, 68–73. [Google Scholar]
  16. Chen, L.Q.; Zu, J.W.; Wu, J. Transverse vibrations of an axially accelerating viscoelastic string with geometric nonlinearity. J. Eng. Math. 2004, 48, 171–182. [Google Scholar] [CrossRef]
  17. Chen, L.Q. Analysis and Control of Transverse Vibrations of Axially Moving Strings. Appl. Mech. Rev. 2005, 58, 91–116. [Google Scholar] [CrossRef]
  18. Gevko, I.; Lyashuk, O.; Djachun, A.; Dovbush, A. Interpretation of the choice of conveyers with improved technological characteristics. MOTROL. In Commission of Motorization and Energetics in Agriculture; Polska Akademia Nauk: Lublin–Rzeszow, Poland, 2015; Volume 17, pp. 107–116. [Google Scholar]
  19. Lyashuk, O.; Vovk, Y.; Sokil, B.; Klendii, V.; Ivasechko, R.; Dovbush, T. Mathematical model of a dynamic process of transporting a bulk material by means of a tube scraping conveyor. Agric. Eng. Int. CIGR J. 2019, 21, 74–81. [Google Scholar]
  20. Sokil, M. Nonlinear vibrations of flexible tubular bodies, along which a continuous flow of medium moves. Collection of scientific and technical works. Sci. Bull. Natl. Min. Univ. 2014, 24, 351–356. [Google Scholar]
  21. Sokil, B.; Senyk, A.; Sokil, M.; Andrukhiv, A.; Kovtonyuk, M.; Gromaszek, K.; Ziyatbekova, G.; Turgynbekov, Y. Mathematical models of dynamics of friable media and analytical methods of their research. Prz. Elektrotechniczny 2019, 95, 74–78. [Google Scholar] [CrossRef]
  22. Kharchenko, Y.V.; Sokil, M.B. Vibrations of Moving Nonlinearly Elastic Media and Asymptotic Method in Their Study; Science Bulletin; Coll. Scientific and Technical Works: Lviv, Ukraine, 2006; Volume 16, pp. 134–138. [Google Scholar]
  23. Delta Function. Mathematics. Available online: https://mathworld.wolfram.com/DeltaFunction.html (accessed on 17 July 2020).
  24. Markovych, B.M. Equations of Mathematical Physics; Lviv Polytechnic Publishing House: Lviv, Ukraine, 2010; p. 384. [Google Scholar]
  25. Oleynik, O.A. Lectures on Partial Differential Equations; Binomial: Moscow, Russia, 2005. [Google Scholar]
  26. Sokil, B.I.; Khytriak, O.I. Inverse Problems of Dynamics on Vibrations of Media Described by a Perturbed Breserton Equation. Scientific Herald of the NLTU of Ukraine; Coll. Scientific and Technical Works: Lviv, Ukraine, 2010; Volume 20, pp. 290–296. [Google Scholar]
  27. Goroshko, A.V.; Roizman, V.P. Statistical Methods for Providing the Stability of the Solutions of Inverse Problems and their Application to Decrease Rotor Vibroactivity. J. Mach. Manuf. Reliab. 2015, 44, 232–238. [Google Scholar] [CrossRef] [Green Version]
  28. Ryabchuk, N.; Grishko, N.; Grishko, V.; Rudenko, A.; Petryk, V.; Bapiyev, I.; Fedushko, S. Artificial Intelligence Technologies Using in Social Engineering Attacks. Available online: http://ceur-ws.org/Vol-2654/paper43.pdf (accessed on 19 August 2020).
  29. Sokolov, G.; Syerov, Y. Building a Computer Model of an Acoustic Signal Recognition Device. Available online: http://ceur-ws.org/Vol-2654/paper37.pdf (accessed on 19 August 2020).
  30. Goroshko, A.V.; Royzman, V.P.; Barmina, O.V. Study of the Stability and Accuracy of the Eccentricity Identification Algorithm in Fast-Revolving Rotor Balancing Problems. J. Mach. Manuf. Reliab. 2016, 45, 227–232. [Google Scholar] [CrossRef] [Green Version]
  31. Royzman, V.; Goroshko, A. Multiple inverse problem. J. Vibroeng. 2012, 14, 1417–1424. [Google Scholar]
Figure 1. (a) Relation of the frequency of nonlinear bending vibrations of the elastic body on the amplitude a and continuous flow V of homogeneous medium (CFHM) relative speed V and (b) the length of the elastic body and CFHM bulk mass at a constant relative speed of the latter.
Figure 1. (a) Relation of the frequency of nonlinear bending vibrations of the elastic body on the amplitude a and continuous flow V of homogeneous medium (CFHM) relative speed V and (b) the length of the elastic body and CFHM bulk mass at a constant relative speed of the latter.
Electronics 10 00040 g001
Figure 2. Laws of changes in the amplitude of bending vibrations during the transition through internal resonance at different values of the relative amount of motion CFHM ( k = m s V ): (a) - - - -   k = 200   kg s m , –– –– –– k = 150   kg s m , –––– k = 50   kg s m ,   m t = 30 ; (b) - - - -   k = 400   kg s m , –– –– –– k = 300   kg s m , –––– k = 100   kg s m ,   m t = 30 ; (c) - - - -   k = 400   kg s m , –– –– –– k = 300   kg s m , –––– k = 100   kg s m ,   m t = 20 ; (d) - - - -   k = 600   kg s m , –– –– –– k = 450   kg s m , –––– k = 10   kg s m ,   m t = 20 .
Figure 2. Laws of changes in the amplitude of bending vibrations during the transition through internal resonance at different values of the relative amount of motion CFHM ( k = m s V ): (a) - - - -   k = 200   kg s m , –– –– –– k = 150   kg s m , –––– k = 50   kg s m ,   m t = 30 ; (b) - - - -   k = 400   kg s m , –– –– –– k = 300   kg s m , –––– k = 100   kg s m ,   m t = 30 ; (c) - - - -   k = 400   kg s m , –– –– –– k = 300   kg s m , –––– k = 100   kg s m ,   m t = 20 ; (d) - - - -   k = 600   kg s m , –– –– –– k = 450   kg s m , –––– k = 10   kg s m ,   m t = 20 .
Electronics 10 00040 g002
Figure 3. Change in the amplitude of bending vibrations during the transition through internal resonance at different values of the amplitude of longitudinal vibrations (a) (- - - b = 0.0015 m, – – – b = 0.001 m, ––– b = 0.0125 m); and different speeds of motion CFHM (b) (- - -V = 1 m/s, – – – V = 1.25 m/s, ––– V = 5 m/s).
Figure 3. Change in the amplitude of bending vibrations during the transition through internal resonance at different values of the amplitude of longitudinal vibrations (a) (- - - b = 0.0015 m, – – – b = 0.001 m, ––– b = 0.0125 m); and different speeds of motion CFHM (b) (- - -V = 1 m/s, – – – V = 1.25 m/s, ––– V = 5 m/s).
Electronics 10 00040 g003
Figure 4. Dependences of the resonant amplitude of vibrations of an elastic body on the frequency unbalance at different values of the nonlinear component of the reducing force: k 2 < 0 , k 2 > 0 .
Figure 4. Dependences of the resonant amplitude of vibrations of an elastic body on the frequency unbalance at different values of the nonlinear component of the reducing force: k 2 < 0 , k 2 > 0 .
Electronics 10 00040 g004
Figure 5. Time change of the amplitude of resonant vibrations due to periodic pulse action at different points of application. (a) ----   k = 400   kg s m , –– –– –– k = 300   kg s m , –––– k = 100   kg s m ,   y 0 = 1 4 , l = 8   m ; (b) ----   k = 400   kg s m , –– –– –– k = 300   kg s m , –––– k = 100   kg s m ,   y 0 = 1 4 , l = 6   m ;   (c) ----   k = 600   kg s m , –– –– –– k = 450   kg s m , –––– k = 150   kg s m ,   y 0 = 1 6 , l = 8   m ; (d) ----   k = 600   kg s m , –– –– –– k = 450   kg s m , –––– k = 150   kg s m ,   y 0 = 1 6 , l = 6   m ; (e) ----   k = 600   kg s m , –– –– –– k = 450   kg s m , –––– k = 150   kg s m ,   y 0 = 1 12 , l = 8   m ; (f) ---   k = 600   kg s m , –– –– –– k = 450   kg s m , –––– k = 150   kg s m ,   y 0 = 1 12 , l = 6   m ).
Figure 5. Time change of the amplitude of resonant vibrations due to periodic pulse action at different points of application. (a) ----   k = 400   kg s m , –– –– –– k = 300   kg s m , –––– k = 100   kg s m ,   y 0 = 1 4 , l = 8   m ; (b) ----   k = 400   kg s m , –– –– –– k = 300   kg s m , –––– k = 100   kg s m ,   y 0 = 1 4 , l = 6   m ;   (c) ----   k = 600   kg s m , –– –– –– k = 450   kg s m , –––– k = 150   kg s m ,   y 0 = 1 6 , l = 8   m ; (d) ----   k = 600   kg s m , –– –– –– k = 450   kg s m , –––– k = 150   kg s m ,   y 0 = 1 6 , l = 6   m ; (e) ----   k = 600   kg s m , –– –– –– k = 450   kg s m , –––– k = 150   kg s m ,   y 0 = 1 12 , l = 8   m ; (f) ---   k = 600   kg s m , –– –– –– k = 450   kg s m , –––– k = 150   kg s m ,   y 0 = 1 12 , l = 6   m ).
Electronics 10 00040 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Andrukhiv, A.; Sokil, M.; Fedushko, S.; Syerov, Y.; Kalambet, Y.; Peracek, T. Methodology for Increasing the Efficiency of Dynamic Process Calculations in Elastic Elements of Complex Engineering Constructions. Electronics 2021, 10, 40. https://doi.org/10.3390/electronics10010040

AMA Style

Andrukhiv A, Sokil M, Fedushko S, Syerov Y, Kalambet Y, Peracek T. Methodology for Increasing the Efficiency of Dynamic Process Calculations in Elastic Elements of Complex Engineering Constructions. Electronics. 2021; 10(1):40. https://doi.org/10.3390/electronics10010040

Chicago/Turabian Style

Andrukhiv, Andriy, Maria Sokil, Solomiia Fedushko, Yuriy Syerov, Yaryna Kalambet, and Tomas Peracek. 2021. "Methodology for Increasing the Efficiency of Dynamic Process Calculations in Elastic Elements of Complex Engineering Constructions" Electronics 10, no. 1: 40. https://doi.org/10.3390/electronics10010040

APA Style

Andrukhiv, A., Sokil, M., Fedushko, S., Syerov, Y., Kalambet, Y., & Peracek, T. (2021). Methodology for Increasing the Efficiency of Dynamic Process Calculations in Elastic Elements of Complex Engineering Constructions. Electronics, 10(1), 40. https://doi.org/10.3390/electronics10010040

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop