The dynamics of eccentric vibration mechanism (Part 1)
A. Leonid Igumnov1 , S. Vladimir Metrikin2 , V. Irina Nikiforova3
1, 2Research Institute for Mechanics, National Research Lobachevsky State University of Nizhni Novgorod, Nizhny Novgorod, Russia
3National Research Lobachevsky State University of Nizhni Novgorod, Nizhny Novgorod, Russia
Journal of Vibroengineering, Vol. 19, Issue 7, 2017, p. 4854-4865.
Received 14 March 2017; received in revised form 23 August 2017; accepted 31 August 2017; published 15 November 2017
This paper presents a mathematical model of a new multi-striker eccentric shock-vibration mechanism (ESVM) with a crank-sliding bar vibration exciter (CSVE) and an arbitrary number of pistons. Analytical solutions for the parameters of the model are obtained to determine the regions of existence of stable periodic motions. Under the assumption of an absolutely inelastic collision of the piston, we derive equations that single out a bifurcational unattainable boundary in the parameter space, which has a countable number of arbitrarily complex stable periodic motions in its neighborhood. We present results of numerical simulations, which illustrate the existence of periodic and stochastic motions. The methods proposed in this paper for investigating the dynamical characteristics of the new crank-type conrod mechanisms allow practitioners to indicate regions in the parameter space, which allow tuning these mechanisms into the most efficient periodic mode of operation, and to effectively analyze the main changes in their operational regimes when the system parameters are changed.
Keywords: point mapping, stability, bifurcation diagram, coefficient of restitution, impact interaction, unreachable boundary.
Recently, alongside with shock-vibration mechanisms with unbalance vibration exciters, eccentric shock-vibration mechanisms (ESVM) with a crank-sliding bar vibration exciter (CSVE) [1, 2] have come to be widely used in civil engineering. This innovative design was based on the principle of “the inverted vibrator”, with its working unit, unbalanced mass, hinged on an eccentric shaft and balanced in rotation by the unbalanced mass. The loading pulse transmitted onto the surface (ground, piles etc.) is produced both due to a spreader with a shoulder of the eccentric shaft and due to the kinetic energy of the falling of the working unit. The effectiveness of compacting and immersing machines not only depends on the quantity of the energy transmitted to the medium being processed, but, to the larger extent, on the character of transmitting this energy – the pulse “form” which should be varied by rearranging certain dynamic factors of a single loading pulse. It is evident that a dense and, at the same time, strong structure of the ground can be achieved only when specific pressure on the surface of the contact of the working unit and the ground in the process of compaction is increased by degrees, its lower limit being determined by the physical properties of the ground in its initial state (before compaction), and the upper one by the ultimate strength of the ground or by technological conditions. Thus, parameters of such machines and mechanisms must be defined from the conditions close to quasi-plastic interaction. Besides, the pulse rate in each single cycle must be such that it precludes possible development of elastic aftereffects of the processed medium in the intervals between the pulses. Such multi-pulse loading mode can be realized using multi-striker ESVM with CSVE designed to easily regulate working modes by changing the geometry of cinematic connections and to solve the problem of compaction of soils in strained conditions of industrial and civil engineering processes. The variety of problems that require the study of dynamic models with discontinuous or piecewise-continuous non-linearities is enormous. They include dry friction systems, shock-vibration and relay systems, variable-structure control systems, numerous mechanic, electric- and radio-technical systems with non-smooth characteristics of certain elements, etc. The investigators’ interest for such substantially nonlinear systems remains alive. Thus, nonlinear dynamic problems of vibration systems with impact interaction of the elements are currently being analyzed by scientific teams headed by B. I. Van de Vrande, D. N. Van Campen, P. I. Leine, G. W. Luo, V. K. Astashov, V. I. Babitsky, М. I. Feygin et al. [3-18]. These works analyzed the complex dynamics of vibration systems with impact interactions of various nature. Periodical motion modes in non-self-supporting systems with finite (mostly with one) and infinite (quasi-plastic) number of impacts during one period of the periodic motion have been studied. Where possible, the mathematical apparatus of the point mapping method [17-29] was used in the studies. The systems and mechanisms used as vibro-impact rams for various media (soils, concrete, loose materials etc.) as well as for driving in or extracting piles, were designed in the form of unbalance mechanisms. The designs of such mechanisms, especially when being used as vibro-impact rams, often could not cope with their tasks for various reasons. One of such reasons is their being unsuited for working in conditions of limited space, as well as being highly power-consuming.
The innovative design of a shock-vibration mechanism with a random (finite) number of sliding bar-strikers  can provide ramming of, for instance, soil (concrete) with higher quality and lower power consumption in spatially limited conditions. At the same time, the introduction of additional impact couples makes the behavior of the mechanism more complicated. That is why the dynamic analysis of such systems and the study of their properties, accounting for the impact interaction with the processed media, is an important issue of both scientific and practical interest. However, the scientific literature available gives no account of certain important aspects of the behavior of such systems. In the first place, it concerns bifurcations related to the expansion of the phase trajectory of a periodical solution stitched up of separate portions beyond the sub-regions of its definition and having no counterparts in analytical systems.
In this connection, the work presents a mathematical model of multi-striker ESVM with CSVE with and without account for the motion of the processed medium, describes the structure of the phase space of model and, based on this and using the point mapping method, determines analytical relations for finding in the parameter space the boundaries of the regions of existence and stability of periodical motion modes. Additionally, relations for the system parameters are given that make it possible to do preliminary tuning of the mechanism onto to main working mode. Using numerical computations with the help of a software complex developed in the Borland C++ Builder 6, bifurcation diagrams were obtained which make it possible to monitor the main changes of the motion modes of the mechanisms.
2. Problem setting
The impact-vibration mechanism considered in the paper (Fig. 1) consists of the frame 1 encasing the eccentric shaft 2 with flywheels at its ends. On the shaft, there are eccentric mechanisms 3, each consisting of two eccentrics nested into each other with a possibility to change the position of the washers, which makes it possible to regulate the value of eccentricity and phase shift between them ( 1, 2, ..., ). Slides-strikers (SS) 4 are hinged at the free ends of the connecting rods. The eccentric mechanisms in combination with the connecting rods and SS convert the rotary motion of the of the shaft with the constant cyclic rate of rotation of the flywheel into reciprocating motion of the frame relative to stands 5. Each SS strikes its own anvil.
The processed medium (soil, pile, etc.) is represented as an elastically fastened mass with elasticity coefficient . Energy losses in the slide-bars of the frame and the medium is accounted for in the form of viscous friction with equivalent damping coefficients and , respectively.
As the eccentrics can be phase-shifted by angles relative to each other and have, generally speaking, different lengths, the body oscillates relative the mechanism with the largest sum of projections of the geometrical dimensions onto the vertical axis at the moment. The impact on the anvil of one or several SS’s takes place either during the change of SS’s interacting with the anvil, or after the breakage of the body together with the eccentric mechanisms.
Fig. 1. The impact-vibration mechanism
Neglecting the masses of SS’s, connecting rods and cranks, the equation of free (without impact) motion of the system for can be written in the following form:
where – coordinate of the mass centre of the body counted off from the equilibrium position of mass for a non-deformed spring , is deviation value of mass from the same initial point, are values characterizing deviation of bases of the th SS, is free-fall acceleration.
When one of the SS contacts the medium there takes place either a momentary impact interaction, such that if , and , are velocities of the th SS and the medium immediately before and after the impact interaction, respectively, then:
Or simultaneous ( 0) motion of the shock-vibration mechanism with the medium described by the following equation:
In Eqs. (2) the following definitions are used: 01 is impact velocity recovery coefficient, .
The position of the eccentricities of length will be measured by angles , counted off the vertical axis. Then it follows that:
where is distance from the fixture point of the th connecting rod to the base of the th SS, is length of the th connecting rod.
When changing in system (1-4) to dimensionless time , coordinates , of the body and the medium, parameters:
As well as introducing function , for the following equations describing impact-oscillatory motions of the mechanism, accounting for the mobility of the processed medium, are obtained:
-periodical function on Eqs. (6) describing simultaneous motion of the striking mechanism with the medium after the impact of the th SS has the following form:
3. Solution method
Assuming 0, , 0 in Eqs. (5-7), the equations of motion with momentary stops  of one or several SS’s can be rewritten as:
The phase space of system Eq. (8) in coordinates , , is truncated along , and, besides, , (Fig. 2).
Fig. 2. The phase space
All the phase trajectories are situated either along the surface ore above it, which represents (generally speaking) of intersecting cylindrical surfaces ( 1, 2, ..., ). The case of corresponds to free motion of the falling mechanism, and corresponds to impact interaction of one of the SS’s with the stopper.
In the case of impact interactions of SS with the stopper for 0, the representing point in the phase space moves as follows: from region , the representing point gets on surface into a point (see Fig. 2), and then, moving along surface into point , ), leaves it and gets into subspace , where it stays until it gets once again onto surface into a point , etc. If at a time , when the representing point gets onto surface into point , the value before the impact velocity is equal to the velocity of change at that moment, or a limiting case of the value of the velocity recovery coefficient is considered for the impact 0, then the representing point displaces from to point and then will slide along surface until, at time , it reaches boundary 0 of the sliding potion plate. These time intervals correspond to prolonged contacts of SS with the stopper.
In what follows, the main attention will be paid to periodical motion modes with successive impacts of each SS during period . It is evident that such a mode is only possible under the condition of paired intersection of two successive surfaces , . For this to be possible, the following relations for the system parameters must hold:
If for any pair , this inequality does not hold (no intersection of the surfaces) then, in the technological process of the motion of the mechanism, the mode with the stroke of the th or the th SS is skipped. To this end, conditions Eq. (9) are of paramount importance when adjusting a specific mechanism to the mode with alternating strokes of each of the SS’sN.
4. Constructing a point map
It follows from the description of the structure of phase space in system Eq. (8) that the dynamics of the mathematical model considered can be investigated by studying the properties of the point transform [18-19, 24-29] of surface into itself.
Let us introduce into the analysis the simplest point transforms transforming surface points into surface points , which, using (8), can be written in the following form ( 0):
The parameters in Eq. (10) have to satisfy relations Eq. (9) as well as the conditions of existence of point transforms , namely:
Point transform of point into point under conditions Eq. (11) is defined as:
5. Existence and stability conditions for fixed points at a point map t corresponding to periodic two-impact motion modes
If is the number of strokes of the th SS against the stopper during one period , and is multiplicity of period Г of the periodic motion mode to period , then the periodic motion will be characterized by numbers , 1, 2,..., . It means that the problem of analyzing the periodic motion mode with a single stroke of each of the SS’s can be reduced to analyzing the properties of point mapping and the conditions of existence of the latter (-bifurcations), analogous to analyzing systems with discontinuous nonlinearities [18, 19]. Thus, the coordinates of a fixed point corresponding to the periodic motion mode with successive strokes of each of the SS’s are determined from system of Eq. (8) supplemented by periodicity conditions:
Using simple transformations, it can be found from system (8) and (13) that:
where the components of -dimensional vector-function are independent of , but are functions of , and the parameters of the system:
Times , 1, 2, 3,..., are determined by solving the system of nonlinear equations of the form:
Stability in the small of the main of the periodic motion modes depends on the value of the roots of the characteristic equation:
where coefficients 0, 1, 2 are defined, as usual , after the linearization of the equations of point transform of in the vicinity of a fixed point. After some simple transformations one obtains:
where , , , are square matrices, with non-zero entries having the following form:
Let us consider in a more detail the dynamics of the mechanism for various values of of SS’s.
6. The dynamics of single piston mechanism
1. Using Eqs. (14-17), the coordinates of the fixed point are determined through the parameters of the system as:
where , and the coefficients of the characteristic equation are.
The regions of existence and stability of periodic single-stroke motion modes are limited by surfaces , , whose equations in the parametric form have, respectively, the following form:
Fig. 3 shows subdivision of the parameter plane , into the planes (between the upper and the lower boundaries) of existence and stability of a periodic single stroke -tuple motion mode.
Fig. 3. The boundaries of the region of existence and stability
In Fig. 4(a) bifurcation diagram for parameter is shown for and 0, 12 wherefrom it follows that a single-stroke motion mode of the mechanism does exist on the interval for the frequency parameter presented in Fig. 3. In the region of the parameters that do not belong to this interval, there exist, as was expected, periodic motions with two (Feigenbaum doubling bifurcation ), three, etc. strokes of SS for the period of motion.
Fig. 4. The bifurcation diagram for parameter
It is to be noted that for 1°, is an analytical function and, thus, by substituting variable in Eqs. (8), a mathematical model of a conventional impact oscillator is obtained. Its dynamics has been investigated in numerous works (see [4, 5] and the references thereof).
It is known  that for 0 point mapping Eq. (10) is reduced to analyzing the point mapping of a circle into a circle. In what follows, for a non-monotonous point transform of a circle into a circle, the structure of the boundaries of the existence and stability regions in the large is studied, and the theorem of un-attainability of these boundaries is proved [21, 22].
Consider a point mapping of a circle into a circle:
where nonlinear 1-periodic function , satisfies the following conditions:
А. Derivative exists everywhere except, perhaps, , where existence of one-sided derivatives is assumed, as well as in points 1.
В. , , , , , .
As function is non-monotonous, the mapping of Eq. (19) for can have simple, only one-turn fixed points and , the first one always being stable due to conditions B and C, whereas the second one always being unstable. Further, if not specially stated, the case of will be considered. It is to be noted that mapping Eq. (19) is monotonous if (see condition B), and, as derivative has in point 0. Equation –1 defines in the plane of parameters , boundary of monotonous transform .Define the region of the values of parameters , , where is monotonous through .
Coordinates , of simple fixed points of mapping Eq. (19) are found from condition and their boundaries of existence are:
It is evident that equations 1, 1 define in the plane of parameters the boundaries of birth and disappearance of the fixed points, respectively. Let signify the region of existence of fixed mapping points Eq. (19), limited in plane , by Eq. (20).
Based on the above-mentioned properties of function and transform , it can be shown that
1. For , despite function being non-monotonous, all the points of a unit segment are transformed into a stable fixed point , except for .
2. For , counting number of the points of a unit segment is transformed into an unstable fixed point.
3. Equation in the plane of parameters , defines the bifurcation curve:
Depicting in the region of stability in the large. In Eq. (21), the value of is determined from 1, where .
Boundary Eq. (21), defined as , is shown to be an unreachable boundary, that is, in any of its vicinities there are bifurcation boundaries corresponding to creation of transform points of .
Define from equation 0, i.e. (under conditions (A, B, C) and for the above mentioned properties of point transform, is defined uniquely).
Consider the sequence of inverse mappings with the values of on the interval . As function is monotonous on this interval, sequence is defined uniquely.
Let . From the definition of it follows that , as well as the evident equalities:
Using relations Eq. (22), as well as conditions A, B, C it is shown that the following lemma holds:
Lemma. If , then
1. , 0, 1, ..., , 0.
2. In points the function has local maximum.
Proof. 1. The first assumption directly follows from the equalities:
2. Let , .
Then, for the values of containing a derivative of function , we have:
Let us assume that in points has a local maximum ( may not exist in these points). As if the conditions of lemma hold, then, due to (A), . Hence, are points of local maximum for transform .
Let . Then 0, i.е. changes its sign after crocking this point. But, in this case because for there is .
Thus, if has local maximums in points then also has local maximums in points . Keeping in mind that has a local maximum in point 0, induction provides a proof of the lemma.
Defines on the plane of parameters a bifurcation curve corresponding to the creation of fixed points of mapping , i.е. to the creation of a cycle of -tuple mapping points .
The bifurcation curve on the plane of parameters that corresponds to the appearance of fixed mapping points of under condition Eq. (23) will be denoted by . The equation of this curve is:
Comparing the equations of bifurcation curves and of curve and also bearing in mind relation Eq. (23), one can postulate.
Theorem. For any , there exists natural , such that, for any , bifurcation curves lie in the plane of parameters , in the vicinity of curve . – unreachable boundary.
Comment 1. An unreachable boundary exists even for more general assumptions, in particular, without limitations on the sign of and without condition .
Comment 2. If is an analytical function then, after crossing bifurcation boundary two cycles of -tuple points are created: a stable one and a non-stable one, i.е. in this case any vicinity of curve comprises an even number of stability regions of cycles of multiple points.
Comment 3. However, in realistic systems, is normally a smooth nonlinear function. Then, evidently, boundary will also be unreachable, and after crossing each of cycles of -tuple points will be created, one of them being stable. In this case, is boundary of capture [24, 25], region is region of stability in the large of the main mode, for adjusting and maintaining of which parameters , are to be chosen in such a way that, if they change, conditions of technological stability was provided [26, 27], or distance of the working point from was considerably larger than . Boundary in this case is dangerous [28, 29].
Comment 4. It is to be noted that for certain conditions imposed on one can absolutely analogously consider regions of existence of m-revolution simple fixed points, investigate their stability and write an equation of the curve analogous to . In doing so it is easily shown that boundaries , 1, 2, …. are unreachable.
The work describes:
– a model of a new multi-striker eccentric shock-vibration mechanism with a crank-sliding bar vibration exciter and an arbitrary number of pistons;
– the phase space of the mathematical model;
– equations of a general (for a random number of SS) point mapping of a two-dimensional non-analytical intersecting surface into itself;
– equations (in the parametric form) of the coordinates of the fixed points corresponding to periodic motions with alternating strokes of each of the SS’s;
– bifurcation boundaries of existence of the above mentioned equations of point mappings.
The main results achieved:
– a methodology and algorithms of constructing point mapping equations of non-analytical secant surfaces truncated along phase coordinates have been developed;
– the parameter space has been subdivided into regions of stability of periodic motion modes;
– the study of dynamics of a single-piston mechanism undergoing a perfectly inelastic collision has been formulated as a Poincare map of a circle into itself;
– the existence of a denumerable number of stable periodic modes of motion in the neighborhood of unreachable boundaries in the extended phase space has been proven;
– bifurcation values of the parameters have been demonstrated, for which chaotic motion modes exist. The scenario of the appearance of chaotic motion modes is very similar to the scenario of doubling of Feigenbaum period [23-28].
The work was financially supported by Russian Scientific Foundation (Grant No. 16-19-10237).
- Bobylev L. М. Compaction of Soils of Backfilling in Strained Working Conditions. Stroyizdat, 1980, p. 252, (in Russian). [Search CrossRef]
- Shilkov V. А., Savalyuk А. D., Metrikin V. S., Polyakov А. А., Shabardin А. K., Alyokhin А. I., Omenenkoi Ya А.с. 1020479 USSR, A Vibro-ram, No. 3376593/29-33, 2005, (in Russian). [Search CrossRef]
- Kobrinskiy A. A., Kоbrinskiy А. Е. Vibro-Impact Systems. Science, 1973, p. 592, (in Russian). [Search CrossRef]
- Babitskiy V. I. The Theory of Vibro-Impact Systems. Science, 1978, p. 352, (in Russian). [Search CrossRef]
- Babitskiy V. I., Krupenin V. L. Oscillations in Strongly Nonlinear Systems. Science, 1985, p. 320, (in Russian). [Search CrossRef]
- Biderman V. L. The Applied Theory of Mechanical Oscillations. High School, 1972, p. 416, (in Russian). [Search CrossRef]
- Bogodukhov S. I., Grebenyuk V. F., Proskurin А. D. Processing of strengthened surfaces in mechanical engineering a maintenance. Mechanical Engineering, 2005, p. 256, (in Russian). [Search CrossRef]
- Vagapov I. K., Ganiyev М. М., Shinkaryov A. S. Theoretically and experimentally investigating the dynamics of an ultrasonic vibro-impact system with an intermediate striker. Mechanical Engineering, Vol. 5, 2008, p. 3-24, (in Russian). [Search CrossRef]
- Zakrzhevskiy М. V. Oscillations of a Substantially Nonlinear Substantially-Nonlinear Mechanical Systems. Zinatne, Riga, 1980, p. 190. [Search CrossRef]
- Babitsky V. I. Theory of Vibro-Impact Systems. Springer, New York, 1998, (in Russian). [Publisher]
- Cveticanin L. The motion of a two-mass system with non-linear connection. Journal of Sound and Vibration, Vol. 252, Issue 2, 2002, p. 361-369. [Search CrossRef]
- Pavloaskaia E., Wiercigroch M. Periodic solution finder for an impact oscillator with a drift. Journal of Sound and Vibration, Vol. 267, Issue 4, 2003, p. 893-911. [Search CrossRef]
- Vagapov I. K., Ganiev M. M., Shinkarev A. S. Forced vibrations of two nonlinearly connected solid waveguides under static load. Journal of Sound and Vibration, Vol. 302, 2007, p. 425-141, (in Russian). [Search CrossRef]
- Luo G. W., Ma X. H. L. Periodic – impact motions and bifurcations in dynamics of a plastic impact oscillator with a frictional slider. Journal Mechanics A/Solids, Vol. 27, 2008, p. 1088-1107. [Search CrossRef]
- Leine R. I., Heimsch T. F. Global uniform symptotic attractive stability of the non-autonomous bouning ball system. Journal Physica, Vol. 241, 2012, p. 2029-2041. [Search CrossRef]
- Goebel R., Sanfelice R. G., Teel A. R. Invariance principles for switching systems via hybrid systems techniques. Systems and Control Letters, Vol. 57, 2008, p. 980-986. [Search CrossRef]
- Zheleztsov N. А. The point transform method and the problem of forced oscillations of an oscillator with “combined” friction. PММ, Vol. 1949, Issue 13, 1, p. 3-40, (in Russian). [Search CrossRef]
- Feygin M. I. Forced Oscillations of Systems with Discontinuous Nonlinearities. Science, 1994, p. 285, (in Russian). [Search CrossRef]
- Bespalova L. V. On the theory of a vibro-impact mechanism. ОТN, Vol. 5, 1957, p. 1-14, (in Russian). [Search CrossRef]
- Bykhovsky I. I., Dorokhova А. D., Zaretskiy L. B., Lukomskiy S. I. On certain periodic motions and the structure of the phase space of a vibro-impact system with a constant righting force. Mechanics and Mechanical Engineering, Vol. 2, 1964, p. 161-165, (in Russian). [Search CrossRef]
- Butenina N. N., Metrikin V. S. On the Unreachable Boundaries of the Point Mapping of a Circle Onto Itself in the Collection: Differential and Integral Equations. Gorky University Publishing, 1981, p. 80-84, (in Russian). [Search CrossRef]
- Butenina N. N., Metrikin V. S. Investigating stability of a class of systems of pulse-phase auto-adjustment of the frequency. University News, Instrument Engineering, Vol. 6, 1986, p. 24-28. [Search CrossRef]
- Shuster G. Deterministic Chaos. World, 1988, p. 237, (in Russian). [Search CrossRef]
- Anishhenko V. S., Vadivasova T. E. Lecture Notes on Nonlinear Dynamics. NIC, Regular and Chaotic Dynamics, 2011, p. 516, http://mirknig.com/2012/04/18/lekcii-po-nelineynoy-dinamike.html, (in Russian) [Search CrossRef]
- Nekorkin V. I. Lectures on the Foundations on the Theory of Oscillations. Nizhnij Novgorod, Preprint of the Nizhniy Novgorod State University 2012, p. 311, (in Russian). [Search CrossRef]
- Nejmark Ju I. Dynamical Systems and Controllable Processes. Second Edition, Librokom, 2010, p. 336, (in Russian). [Search CrossRef]
- Kuznecov A. P., Savin A. V., Tjurjukina L. V. Introduction to Physics of Nonlinear Maps. Saratov, 2010, p. 134, http://www.studmed.ru/kuznecov-ap-vvedenie-v-fiziku-nelineynyh-otobrazheniy _4bd0dfbdedc.html, (in Russian). [Search CrossRef]
- Danilov Ju A. Lecture Notes on Nonlinear Dynamics: Basic Introduction. Second Edition, KomKniga, 2006, p. 208, (in Russian). [Search CrossRef]
Articles Citing this One
Determination of inertia-stiffness parameters and motion modelling of three-mass vibratory system with crank excitation mechanism
Vibroengineering PROCEDIAVitaliy Korendiy, Oleksii Lanets, Oleksandr Kachur, Petro Dmyterko, Roman Kachmar