Mathematical modeling of forced oscillations of semidefinite vibroimpact system sliding along rough horizontal surface
Vitaliy Korendiy^{1} , Volodymyr Gursky^{2} , Oleksandr Kachur^{3} , Volodymyr Gurey^{4} , Oleksandr Havrylchenko^{5} , Oleh Kotsiumbas^{6}
^{1, 2, 3, 4, 5}Department of Robotics and Integrated Mechanical Engineering Technologies, Lviv Polytechnic National University, Lviv, Ukraine
^{6}Department of Motor Vehicle Transport, Lviv Polytechnic National University, Lviv, Ukraine
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 39, 2021, p. 164169.
https://doi.org/10.21595/vp.2021.22298
Received 2 November 2021; received in revised form 15 November 2021; accepted 21 November 2021; published 2 December 2021
53rd International Conference on Vibroengineering, Colombo, Sri Lanka and Udaipur, India, December 2, 2021
The paper considers the motion conditions of a semidefinite vibratory system placed upon a rough horizontal surface. Such systems are sometimes called unrestrained or degenerate ones, and are usually used in various vibrationdriven robots and capsules. Unlike the numerous existent investigations dedicated to a similar subject, the novelty of the present paper consists in the implementation of a crank mechanism for exciting oscillations of a doublemass vibroimpact system setting into planar locomotion a robot’s movable body. A general design diagram of the improved semidefinite vibroimpact system is proposed, and the corresponding mechanical diagram is considered. The differential equations describing the system sliding (planar locomotion) along a rough horizontal surface are derived. A thorough analysis of the main inertiastiffness, design, and excitation parameters influencing the system motion conditions is carried out. Performing the numerical modeling in MathCad software, the dynamic behavior of the robot’s movable body is studied under the specified system’s parameters and operational conditions.
 A general design diagram of the improved semidefinite vibroimpact system is proposed, and the corresponding mechanical diagram is considered.
 The differential equations describing the system sliding (planar locomotion) along a rough horizontal surface are derived.
 A thorough analysis of the main inertiastiffness, design, and excitation parameters influencing the system motion conditions is carried out.
 The dynamic behavior of the robot’s movable body is studied under different system’s parameters and operational conditions.
Keywords: vibrationdriven robot, crank mechanism, design diagram, differential equations, excitation, numerical modeling, dynamic behavior, operational conditions, optimization.
1. Introduction
It is wellknown that the motion of a semidefinite vibratory system is composed of translation and vibration. In such a case, one of the system’s natural frequencies is equal to zero; this means that the system doesn’t oscillate, and translationally moves as a single rigid body without any relative movements of the masses forming this system [1]. That’s why such systems are usually used in various vibrationdriven robots and capsules.
The vibrationdriven system consisting of a movable rigid body and two internal masses set into oscillatory motion by orthogonally directed periodical forces is studied in [2]. The similar system driven by two orthogonally oscillating internal masses and sliding along a rough horizontal plane is investigated in [3]. In [4], there are considered two types of vibroimpact capsule systems with different types of constraints. The paper [5] is dedicated to the analysis of a stickslip motion of a 2DOF vibrationdriven system.
In [6], the authors studied the dynamic behavior of a semidefinite vibratory system with a rotating internal mass. The improved design diagram of a vibrationdriven system (vibration manipulator) is proposed in [7]. In [8], there is analyzed the dynamic behavior of a vibroimpact locomotion system subjected to the action different friction levels. The paper [9] presents a thorough comparative analysis of two vibrationdriven locomotion systems: a common vibratory system and a vibroimpact one.
The novelty of the present paper consists in the implementation of a crank mechanism for exciting oscillations of a doublemass vibroimpact system setting into planar locomotion a robot’s movable body. Some aspects of this topic regarding the sliding motion of the doublemass vibratory system were presented in the paper [10]. The possibilities of implementation of the optimized vibroimpact systems operating under the resonance conditions were considered in [11]. The improved design of a crank mechanism used for exciting oscillations of vibratory equipment was presented in [12]. The present paper is aimed at combining the controllable crank excitation mechanism with a doublemass vibroimpact system in order to provide the effective conditions of the robot motion.
2. Design diagram of semidefinite vibratory system
The proposed design diagram of a semidefinite vibroimpact system is presented in Fig. 1(a). The electric motor 1 with changeable rotor speed sets into motion the eccentric disk 2. The latter is hingedly joined with the connecting rods 3 and 4. The sliding shafts 7 and 8 are hinged to the rods 3 and 4, respectively, and translationally move in the guide (pilot) bearings 5 and 6. The springs 9 and 10 are connected to the free ends of the sliding shafts 7 and 8, and, at the same time, to the impact bodies 11 and 12 sliding along the guide rods 17 and 18. The impact plates 13 and 14 are connected to the robot’s body 19 with the help of the springs 15 and 16.
The main idea of the proposed mobile vibratory system consists in setting the robot’s body into the translational motion by implementing the vibroimpact operation modes. The steadystate rotation of the disk 2 causes the straightline oscillations of the shafts 7 and 8. The corresponding shafts influence the springs 9 and 10 exciting the oscillations of the bodies 11 and 12. The latter ones can impact the plates 13 and 14 under specific operational conditions. The kinetic energy being released during the impact of the body 11 and the plate 13 causes the reduction of the vertical force exerted by the robot’s body upon the supporting surface. Similarly, the kinetic energy being released during the impact of the body 12 and the plate 14 causes both the reduction of the vertical force exerted by the robot’s body upon the supporting surface, and the sliding of the robot’s body along the surface. The mentioned impacts are shifted in time (displaced in phase) and can be easily controlled by changing the gaps between the impact bodies 11, 12 and the corresponding impact plates 13, 14. The angle $\phi $ between the sliding shafts 7 and 8 is also a controllable parameter effecting the system’s operation.
Fig. 1. The proposed semidefinite vibroimpact system: a) general design; b) kinematic diagram
a)
b)
While designing and implementing such vibroimpact systems, the special attention must be paid to substantiation of the inertial, stiffness, design and excitation parameters in accordance with the specific operational conditions in order to ensure the maximal locomotion speed and the minimal energy consumption. To solve the stated problem let us consider the corresponding kinematic diagram of the semidefinite vibroimpact system (Fig. 1(b)). The driving mechanism consists of the crank DA joined with the connecting rods AB and AC. The angular velocity $\mathrm{\omega}$ of the crank is considered to be constant. The connecting rods AB and AC are hinged to the sliders B and C. The corresponding springs are characterized by the stiffness coefficients ${k}_{1}$, ${k}_{2}$, ${k}_{3}$, ${k}_{4}$, and damping coefficients ${c}_{1}$, ${c}_{2}$, ${c}_{3}$, ${c}_{4}$. The impact bodies are modelled as springloaded sliders (oscillating masses) ${m}_{2}$, ${m}_{3}$. The robot’s body is characterized by the mass ${m}_{1}$.
3. Differential equations describing the system motion
The studied semidefinite vibroimpact system is characterized by three degrees of freedom. The robot’s body (mass ${m}_{1}$) is considered to move translationally along a rough horizontal plane in the direction of the $Ox$ axis. The sliders B and C are driven by the crank mechanism rotating at the constant angular velocity $\omega $. The sliders set the masses ${m}_{2}$, ${m}_{3}$ into the oscillatory motion along the axes $D\zeta $ and $D\xi $, respectively. The system of differential equations describing the motions of all the bodies is following:
${m}_{3}\bullet {\ddot{\xi}}_{3}+{c}_{3}\bullet \left({\dot{\xi}}_{3}{\dot{\xi}}_{C}\right)+{k}_{3}\bullet \left({\xi}_{3}{\xi}_{C}\right)+{c}_{1}^{*}\bullet \left({\dot{\xi}}_{3}{\dot{\xi}}_{H}\right)+{k}_{1}^{*}\bullet \left({\xi}_{3}{\xi}_{H}\right)=0,$
$\left({m}_{1}+{m}_{2}+{m}_{3}\right)\bullet {\ddot{x}}_{1}+{c}_{3}\bullet \left({\dot{\xi}}_{3}{\dot{\xi}}_{C}\right)\bullet \mathrm{sin}\phi +{k}_{3}\bullet \left({\xi}_{3}{\xi}_{C}\right)\bullet \mathrm{sin}\phi +{c}_{1}^{*}\bullet \left({\dot{\xi}}_{3}{\dot{\xi}}_{H}\right)\bullet \mathrm{sin}\phi +{k}_{1}^{*}\bullet \left({\xi}_{2}{\xi}_{H}\right)\bullet \mathrm{sin}\phi ={F}_{fr},$
where ${\zeta}_{2}$, ${\dot{\zeta}}_{2}$, ${\ddot{\zeta}}_{2}$ are the coordinate, velocity and acceleration of the mass ${m}_{2}$ moving along the axis $D\zeta $; ${\xi}_{3}$, ${\dot{\xi}}_{3}$, ${\ddot{\xi}}_{3}$ are the coordinate, velocity and acceleration of the mass ${m}_{3}$ moving along the axis $D\xi $; ${\zeta}_{B}$, ${\dot{\zeta}}_{B}$ are the coordinate and velocity of the slider B moving along the axis $D\zeta $; ${\xi}_{C}$, ${\dot{\xi}}_{C}$ are the coordinate and velocity of the slider C moving along the axis $D\xi $; ${\zeta}_{E}$, ${\dot{\zeta}}_{E}$ are the coordinate and velocity of the plate E moving along the axis $D\zeta $; ${\xi}_{H}$, ${\dot{\xi}}_{H}$ are the coordinate and velocity of the plate H moving along the axis $D\xi $; ${\ddot{x}}_{1}$ is the acceleration of the robot’s body moving along the axis $Ox$; ${c}_{1}^{*}$, ${c}_{4}^{*}$ are the functions describing the change of damping coefficients of the springs 16 and 15 during the impacts of the bodies 12–14 and 11–13, respectively (see Fig. 1(a)); ${k}_{1}^{*}$, ${k}_{4}^{*}$ are the functions describing the change of stiffness coefficients of the springs 16 and 15 during the impacts of the bodies 12–14 and 11–13, respectively (see Fig. 1(a)); ${F}_{fr}$ is the dry friction force acting upon the robot’s body during its sliding along a rough horizontal plane.
The coordinate of the slider B moving along the axis $D\zeta $ and driven by the crank DA rotating at the constant angular velocity $\omega $ can be determined as follows:
where ${l}_{DA}$, ${l}_{AB}$ are the lengths of the rods DA and AB, respectively.
The coordinate of the slider C moving along the axis $D\xi $ and driven by the crank DA rotating at the constant angular velocity $\omega $ can be determined as follows:
where ${l}_{AC}$ is the length of the rod AC.
The functions describing the change of damping and stiffness coefficients of the springs 16 and 15 during the impacts of the bodies 12, 14, and 11, 13, respectively, are following:
${c}_{4}^{*}\left(t\right)=\left\{\begin{array}{c}{c}_{4},\left({\zeta}_{2}\left(t\right){\zeta}_{E}\right)\ge 0,\\ 0,\left({\zeta}_{2}\left(t\right){\zeta}_{E}\right)0,\end{array}\right.{k}_{4}^{*}\left(t\right)=\left\{\begin{array}{c}{k}_{4},\left({\zeta}_{2}\left(t\right){\zeta}_{E}\right)\ge 0,\\ 0,\left({\zeta}_{2}\left(t\right){\zeta}_{E}\right)0.\end{array}\right.$
The dry friction force acting upon the robot’s body during its sliding along a rough horizontal plane can be determined as follows:
where ${F}_{k.fr}\left(t\right)$ is the function describing the change of the kinetic friction force acting upon the robot’s body during its sliding along a rough horizontal plane; ${F}_{in}\left(t\right)$ is the function describing the change of all the internal inertial forces acting upon the robot’s body in the $Ox$ axis direction; ${\dot{x}}_{1}$ is the velocity of the robot’s body (mass ${m}_{1}$).
The kinetic friction force can be described by the following function:
${c}_{4}^{*}\left(t\right)\bullet \left({\dot{\zeta}}_{2}\left(t\right){\dot{\zeta}}_{E}\right){k}_{4}^{*}\left(t\right)\bullet \left({\zeta}_{2}\left(t\right){\zeta}_{E}\right){c}_{3}\bullet \left({\dot{\xi}}_{3}\left(t\right){\dot{\xi}}_{C}\left(t\right)\right)\bullet \mathrm{cos}\phi $
${k}_{3}\bullet \left({\xi}_{3}\left(t\right){\xi}_{C}\left(t\right)\right)\bullet \mathrm{cos}\phi {c}_{1}^{*}\left(t\right)\bullet \left({\dot{\xi}}_{3}\left(t\right){\dot{\xi}}_{H}\right)\bullet \mathrm{cos}\phi $
${k}_{1}^{*}\left(t\right)\bullet \left({\xi}_{2}\left(t\right){\xi}_{H}\right)\bullet \mathrm{cos}\phi )\bullet {k}_{fr},$
where ${k}_{fr}$ is the coefficient of kinetic (sliding) friction.
The function describing the change of all the internal inertial forces acting upon the robot’s body in the $Ox$ axis direction can be written as follows:
${c}_{1}^{*}\left(t\right)\bullet \left({\dot{\xi}}_{3}\left(t\right){\dot{\xi}}_{H}\right)\bullet \mathrm{sin}\phi {k}_{1}^{*}\left(t\right)\bullet \left({\xi}_{2}\left(t\right){\xi}_{H}\right)\bullet \mathrm{sin}\phi .$
Substituting Eq. (2)–(7) into Eq. (1), the mathematical model describing forced oscillations of the semidefinite vibroimpact system can be developed.
4. Results of numerical modeling
In order to verify the correctness of the developed mathematical model let us carry out simulation of the system motion in MathCad software. In accordance with the 3D model of the vibroimpact system designed in SolidWorks software, let us specify the following parameters: ${m}_{1}=$ 1 kg, ${m}_{2}=$0.1 kg, ${m}_{3}=$0.08 kg, ${l}_{DA}=$0.02 m, ${l}_{AB}={l}_{AC}=$0.06 m, $\phi =$50° = 0.873 rad, ${k}_{1}={k}_{4}=$ 10^{7} N/m, ${k}_{2}={k}_{3}=$ 10^{3} N/m, ${c}_{1}={c}_{2}={c}_{3}={c}_{4}=$ 0 (N∙s)/m; $\omega =$ 100 rad/s, ${\xi}_{H}={\zeta}_{E}=$ 0.06 m, ${k}_{fr}=$ 0.2.
Fig. 2(a) presents the time dependencies of displacements of the impact bodies (masses ${m}_{2}$ and ${m}_{3}$); the masses are set into oscillatory motion due to the crank rotation at a constant angular velocity. The operation of the vibroimpact system causes the translational motion of the robot’s body whose displacement changes in time according to the dependences presented in Fig. 2(b); the simulation was carried out for different values of the angle $\phi $: 40°, 50°, 60°, 70°, 80°, 90°. Based on the obtained results, it can be concluded that the average speed of the robot’s body motion reaches its maximal value at the angle $\phi =$ 80°. Further increase of the inclination angle over 85…90° causes reduction of the average translational speed.
Fig. 2. Dynamic behavior of the semidefinite vibroimpact system under the specified operational conditions and design parameters: a) displacements of the impact bodies; b) displacements of the robot’s body at different angles $\phi $
a)
b)
Fig. 3. Results of modelling the system motion under different excitation conditions and design parameters: a) the case of changing the impact mass; b) the case of changing the crank length
a)
b)
During further optimization of the studied system it is necessary to consider the possibilities of maximization of the robot translational speed by changing the design and excitation parameters. As the first example, let us consider the change of the impact mass ${m}_{3}$ at the optimal angle $\phi =$ 80° (Fig. 3(a)). The gradual increase of the mass ${m}_{3}$ from 0.08 kg to 0.25 kg causes the reduction of the robot average translational speed, whereas the decrease of the mass ${m}_{3}$ to 0.02 kg ensures the maximal value of the speed. One more parameter effecting the robot speed is the driving mechanism eccentricity value, or the length of the driving crank ${l}_{DA}$. The simulation showed that the optimal values $\phi =$ 80°, ${m}_{3}=$ 0.02 kg provide the quickest motion of the robot’s body under the condition of ${l}_{DA}=$ 0.06 m (Fig. 3(b)). The larger value of eccentricity is ensured, the larger robot average speed can be reached.
Further investigations on the subject of the paper should include the analysis and solution of the system optimization problem. The latter consists in maximizing the robot’s body average speed ${V}_{aver.}$ and minimizing the power consumption ${P}_{cons.}$:
5. Conclusions
The paper substantiates the possibilities of implementation of a crank mechanism for exciting oscillations of a doublemass vibroimpact system setting into planar locomotion a robot’s movable body. The improved design of a semidefinite vibroimpact system is proposed, and the corresponding kinematic diagram is developed. The system consists of two impact bodies (internal masses); one mass moves along a vertical axis, whilst the other one moves along an inclined axis. The masses are springloaded and set into the oscillatory motion be the crank excitation mechanism. The maximal displacements of the masses are limited by the springloaded impact plates. The kinetic energy being released during the impact of the bodies and the plates causes both the reduction of the vertical force exerted by the robot’s body upon the supporting surface, and the sliding of the robot’s body along the surface.
The studied semidefinite vibroimpact system is characterized by three generalized coordinates (degrees of freedom): two of them describe the translational motions of the impact bodies along the corresponding guide axes, and the last one considers the translational motion of the robot’s body along a rough horizontal surface. Based on the developed mathematical model of forced oscillation of the studied semidefinite vibroimpact system, the simulation of the impact masses and the robot’s body motions is carried out in MathCad software. All the input parameters are specified in accordance with the 3D model of the vibroimpact system designed in SolidWorks software. The results of simulation of the system dynamic behavior allowed to draw the conclusions about the certain optimal design and excitation parameters: the optimal angle between the guide axes of the impact bodies $\phi =$ 80°, the optimal value of the inclined impact mass ${m}_{3}=$ 0.02 kg, the optimal driving mechanism eccentricity value (the length of the driving crank) ${l}_{DA}=$ 0.06 m.
References
 Rao S. S., Mechanical Vibrations. Sixth Edition in SI Units. Harlow, United Kingdom: Pearson Education Limited, 2018. [Search CrossRef]
 X. Zhan, J. Xu, and H. Fang, “Planar locomotion of a vibrationdriven system with two internal masses,” Applied Mathematical Modelling, Vol. 40, No. 2, pp. 871–885, Jan. 2016, https://doi.org/10.1016/j.apm.2015.06.016 [Publisher]
 Fritzkowski P. and Starosta R., “Dynamics of an inertially driven robot,” Vibrations in Physical Systems, Vol. 29, pp. 1–9, 2018. [Search CrossRef]
 Y. Yan, Y. Liu, and M. Liao, “A comparative study of the vibroimpact capsule systems with onesided and twosided constraints,” Nonlinear Dynamics, Vol. 89, No. 2, pp. 1063–1087, Jul. 2017, https://doi.org/10.1007/s1107101735007 [Publisher]
 P. Li and Z. Jiang, “Bifurcation analysis of stickslip motion of the vibrationdriven system with dry friction,” Mathematical Problems in Engineering, Vol. 2018, No. 3, pp. 1–10, Aug. 2018, https://doi.org/10.1155/2018/2305187 [Publisher]
 B. Bardin and A. Panev, “On dynamics of a rigid body moving on a horizontal plane by means of motion of an internal particle,” Vibroengineering PROCEDIA, Vol. 8, pp. 135–141, Oct. 2016, https://doi.org/10.1063/1.5034582 [Publisher]
 K. Ragulskis, B. Spruogis, M. Bogdevičius, A. Matuliauskas, V. Mištinas, and L. Ragulskis, “Motion of vibration manipulators with self stopping device in one direction with interactions of two non deformable elements in one direction,” Mechanics, Vol. 26, No. 6, pp. 526–531, Dec. 2020, https://doi.org/10.5755/j01.mech.26.6.25527 [Publisher]
 K.T. Nguyen, N.T. La, K.T. Ho, Q.H. Ngo, N.H. Chu, and V.D. Nguyen, “The effect of friction on the vibroimpact locomotion system: modeling and dynamic response,” Meccanica, Vol. 56, No. 8, pp. 2121–2137, Aug. 2021, https://doi.org/10.1007/s1101202101348w [Publisher]
 N.T. La, T.T. Nguyen, and V.D. Nguyen, “A comparative study on the two vibration driven locomotion systems in various friction levels,” Vietnam Journal of Mechanics, Vol. 43, pp. 121–137, Apr. 2021, https://doi.org/10.15625/08667136/15662 [Publisher]
 V. Korendiy, “Dynamics of twomass mobile vibratory robot with electromagnetic drive and vibroimpact operation mode,” Ukrainian Journal of Mechanical Engineering and Materials Science, Vol. 4, No. 2, pp. 80–93, 2018, https://doi.org/10.23939/ujmems2018.02.080 [Publisher]
 V. Gursky, I. Kuzio, and V. Korendiy, “Optimal synthesis and implementation of resonant vibratory systems,” Universal Journal of Mechanical Engineering, Vol. 6, No. 2, pp. 38–46, Mar. 2018, https://doi.org/10.13189/ujme.2018.060202 [Publisher]
 O. Lanets, O. Kachur, V. Korendiy, and V. Lozynskyy, “Controllable crank mechanism for exciting oscillations of vibratory equipment,” in Lecture Notes in Mechanical Engineering, pp. 43–52, 2021, https://doi.org/10.1007/9783030778231_5 [Publisher]
Cited By
Lecture Notes in Mechanical Engineering
Vitaliy Korendiy, Oleksandr Kachur

2023

Vibroengineering PROCEDIA
Vitaliy Korendiy, Oleksandr Kachur, Volodymyr Gursky, Oleh Kotsiumbas, Petro Dmyterko, Serhij Nikipchuk, Yaroslav Danylo

2022

Vibroengineering PROCEDIA
Vitaliy Korendiy, Volodymyr Gursky, Oleksandr Kachur, Petro Dmyterko, Oleh Kotsiumbas, Oleksandr Havrylchenko

2022

Vibroengineering PROCEDIA
Vitaliy Korendiy, Volodymyr Gurey, Volodymyr Borovets, Oleh Kotsiumbas, Vasyl Lozynskyy

2022

Vibroengineering PROCEDIA
Vitaliy Korendiy, Oleksandr Kachur, Volodymyr Gursky, Volodymyr Gurey, Nadiia Maherus, Oleh Kotsiumbas, Oleksandr Havrylchenko

2022

Vibroengineering PROCEDIA
Vitaliy Korendiy, Oleksandr Kachur, Volodymyr Gursky, Volodymyr Gurey, Roman Pelio, Oleh Kotsiumbas

2022

Vibroengineering PROCEDIA
Vitaliy Korendiy, Igor Kuzio, Serhij Nikipchuk, Oleh Kotsiumbas, Petro Dmyterko

2022
