Mathematical model and motion analysis of a wheeled vibroimpact locomotion system
Vitaliy Korendiy^{1} , Volodymyr Gursky^{2} , Oleksandr Kachur^{3} , Petro Dmyterko^{4} , Oleh Kotsiumbas^{5} , Oleksandr Havrylchenko^{6}
^{1, 2, 3, 4, 6}Department of Robotics and Integrated Mechanical Engineering Technologies, Lviv Polytechnic National University, Lviv, Ukraine
^{5}Department of Motor Vehicle Transport, Lviv Polytechnic National University, Lviv, Ukraine
^{1}Corresponding author
Vibroengineering PROCEDIA, Vol. 41, 2022, p. 7783.
https://doi.org/10.21595/vp.2022.22422
Received 31 January 2022; received in revised form 15 February 2022; accepted 25 February 2022; published 21 April 2022
55th International Conference on Vibroengineering, Almaty, Kazakhstan, April 21, 2022
The paper is aimed at investigating the motion conditions of the wheeled vibroimpact locomotion system equipped with the twin crankslider excitation mechanism and the additional braking mechanisms allowing only oneway rotation of the wheels. The novelty of the present research consists in the improved mathematical model describing the motion conditions of the vibroimpact system and the proposed parameters optimization technique that allows for maximizing the average translational velocity of the wheeled platform. The main idea of this technique is to provide the maximal velocities of internal bodies when they get in contact with the corresponding impact plates. The numerical modeling results describing the dynamic behavior of the vibroimpact system are obtained in Mathematica software and substantiate the correctness of the developed mathematical model and of the proposed parameters optimization technique. The paper can be of significant practical and scientific interest for researchers and engineers studying and improving the vibratory locomotion systems, e.g., for inspecting and cleaning the pipelines.
 The improved design of the wheeled vibroimpact locomotion system equipped with the twin crankslider excitation mechanism and the braking mechanism allowing the wheels oneway rotation is proposed.
 The mathematical model describing the vibroimpact system motion is derived and the parameters optimization technique that allows for maximizing the system average translational velocity is proposed.
 Providing the vibroimpact operational conditions, the stepwise (discontinuous) translational motion of the wheeled platform is studied by means of numerical modeling in the Mathematica software.
 The obtained results were used while designing and implementing the improved vibratory locomotion system for inspecting and cleaning the pipelines.
Keywords: twin crankslider mechanism, braking mechanism, motion conditions, mathematical model, optimization technique, numerical modeling, dynamic behavior.
1. Introduction
The investigations dealing with mathematical modeling and analyzing the dynamic behavior of the vibrationdriven locomotion systems are currently of significant interest among researchers and engineers all over the world. One of the most important problems occurring during designing and implementation of such systems consists in providing their most efficient operational conditions. Considering the wheeled locomotion systems, the average velocity factor is of the most significant ones that should be maximized under the other constant conditions. This paper is focused on motion modeling and testing the parameters optimization technique of the wheeled vibroimpact locomotion system with the twin crankslider excitation mechanism.
Different types of the vibroimpact locomotion systems are thoroughly investigated in numerous scientific papers, e.g. [1][4]. The experimental investigations on the selfpropelled vibroimpact locomotion system operating under anisotropic friction are presented in [1]. The influence of dry and isotropic friction on the dynamic behavior of a vibroimpact locomotion system is studied in [2]. The paper [3] is devoted to investigating the locomotion conditions of a selfpropelled vibrationdriven capsule system equipped with the pushtype solenoid. In [4], the authors considered the basic approaches to improving the locomotion performance of vibrationdriven systems. The paper [5] is dedicated to the parameters optimization problem focused on maximizing the average velocity of the multimodule vibrationdriven locomotion robot. The wheeled vibrationdriven robot equipped with the inertial vibration exciter (unbalanced rotor) is studied in the papers [6][8]. The robot dynamic behavior is thoroughly studied in [6], [7], and the problems of its vibratory system’s parameters optimization are considered in [8].
The analysis of numerous scientific papers (e.g. [1][8]) has shown that the problems of implementing the crank excitation mechanisms in the multimass vibroimpact locomotion systems are not thoroughly investigated. Therefore, the authors of the present paper have previously studied the improved design of the controllable eccentrictype excitation mechanism [9], investigated the motion conditions of the threemass vibroimpact system with the crankslider excitation mechanism [10], and analyzed the possibilities of implementing the proposed excitation mechanism in vibratory compacting equipment [11]. The novelty of the present paper consists in improving the mathematical model and parameters optimization technique of the wheeled vibroimpact locomotion system equipped with the twin crankslider excitation mechanism and the additional braking mechanisms allowing only oneway rotation of the wheels. The proposed wheeled platform can be implemented in practice for inspecting and cleaning the pipelines.
2. Research methodology
2.1. Dynamic diagram of the wheeled vibroimpact locomotion system
The dynamic diagram of the threemass vibroimpact system is presented in Fig. 1. The inertial reference frame $xOy$ is adopted to study the translational locomotion of the wheeled platform of the mass ${m}_{1}$. For this purpose, the generalized coordinate ${x}_{1}$ is introduced. The twin crankslider mechanism “ABCADE” is used to excite the oscillations of the internal bodies of the masses ${m}_{2}$ and ${m}_{3}$. The angular velocity $\omega $ of the crank BAD is considered to be constant, and the initial positions of the crank and the connecting rods BC, DE are the horizontal ones. The relative displacements of the oscillating bodies ${m}_{2}$, ${m}_{3}$ with respect to the wheeled platform are described by the corresponding generalized coordinates ${x}_{2}$ and ${x}_{3}$. The rotation of the crank BAD is transformed into the translational reciprocating motions of the sliders C and E, which are connected with the internal bodies ${m}_{2}$, ${m}_{3}$ by the springs of the stiffnesses ${k}_{1}$, ${k}_{2}$, respectively. The potential energies accumulated in the springs during the sliders reciprocating motions are transformed into the kinetic energies of the internal bodies oscillations.
In order to ensure the translational motion of the wheeled platform, the vibroimpact operational conditions are used. The relative displacements of the internal bodies (masses ${m}_{2}$, ${m}_{3}$) with respect to the wheeled platform (mass ${m}_{1}$) are limited by the corresponding impact plates and springs ${k}_{3}$, ${k}_{4}$. When the internal bodies come into contact with the impact plates, their kinetic energies reduce causing the increase in the potential energies of the impact springs. The energies accumulated in the springs are transformed into the kinetic energy of the wheeled platform. In such a manner, the stepwise (discontinuous) translational motion of the platform is performed to the right; the latter is adopted as the forward motion direction. The backward motion (directed to the left) is constrained due to the use of the special braking mechanisms in the wheels allowing only oneway rotation. For example, these mechanisms can be designed as unidirectional (overrunning) mechanical clutches (Fig. 1). While carrying out further parametric synthesis of the vibroimpact system, it is necessary to maximize the average velocity of the wheeled platform.
Fig. 1. Dynamic diagram of the vibroimpact locomotion system
2.2. Mathematical model of the system locomotion
In order to develop the mathematical model of the considered vibroimpact system locomotion, the corresponding Lagrange equation of the second kind was derived for each generalized coordinate (${x}_{1}$, ${x}_{2}$, ${x}_{3}$). At the present stage of the investigations, all the energy losses are neglected. Three differential equations describing the system locomotion are following:
where ${x}_{C}$, ${x}_{E}$ are the functions describing the relative displacements of the sliders C, E with respect to their initial positions; ${k}_{3}^{*}$, ${k}_{4}^{*}$ are the functions describing the change of the stiffness coefficients during the impacts of the corresponding bodies; ${\delta}_{2}^{*}$, ${\delta}_{3}^{*}$ are the initial distances between the oscillating bodies and the corresponding impact plates (i.e., the initial impact gaps); ${F}_{br}$ is the braking force acting on the wheels when they attempt to rotate in a counterclockwise direction. All the functions and values mentioned above can be expressed as follows:
where ${l}_{AB}$, ${l}_{BC}$, ${l}_{AD}$, ${l}_{DE}$ are the lengths of the corresponding rods AB, BC, AD, DE forming the crankslider excitation mechanism; $\omega $ is the angular velocity of the crank BAD; ${\delta}_{2}$, ${\delta}_{3}$ are the smallest distances between the internal bodies and the corresponding impact plates at the state of rest when the crank angular position is 0 (for the mass ${m}_{2}$) and $\pi $ (for the mass ${m}_{3}$).
2.3. Analyzing the basic stages of the internal bodies oscillations
Substituting Eq. (4)(9) into Eq. (2)(3), let us solve the obtained system of two differential equations in the Mathematica software considering the following stages. At the first stage, the second body (mass ${m}_{2}$) moves from its left initial position to the position where it comes into contact with the impact plate: ${x}_{2\left(1\right)}\left({t}_{2\left(1\right)}\right)={\delta}_{2}^{*}=2{l}_{AB}+{\delta}_{2}$. The third body (mass ${m}_{3}$) moves from its right initial position to the position where its velocity becomes zero: ${\dot{x}}_{3\left(1\right)}\left({t}_{3\left(1\right)}\right)=0$. For this stage, adopting zero values of the initial conditions $\left({x}_{2}\left(0\right),{\dot{x}}_{2}\left(0\right),{x}_{3}\left(0\right),{\dot{x}}_{3}\left(0\right)=0\right)$, the laws of motion of the corresponding internal bodies are following:
At the second stage, the second body (mass ${m}_{2}$) moves from the position where it got in contact with the impact plate to the position where its velocity becomes zero: ${\dot{x}}_{2\left(2\right)}\left({t}_{2\left(2\right)}\right)=0$. The third body (mass ${m}_{3}$) moves from its left position where its velocity became zero to the position where it comes into contact with the impact plate: ${x}_{3\left(2\right)}\left({t}_{3\left(2\right)}\right)={x}_{3\left(1\right)}\left({t}_{3\left(1\right)}\right)+{\delta}_{3}^{*}={x}_{3\left(1\right)}\left({t}_{3\left(1\right)}\right)+{\delta}_{3}$. For the considered stage, the initial conditions can be presented as follows: ${x}_{2}\left(0\right)={x}_{3}\left(0\right)={\dot{x}}_{3}\left(0\right)=0$, ${\dot{x}}_{2}\left(0\right)={\dot{x}}_{2\left(1\right)}\left({t}_{2\left(1\right)}\right)$. The laws of motion of the internal bodies are following:
$+\sqrt{{m}_{2}/\left({k}_{1}{k}_{3}\right)}{\dot{x}}_{2\left(1\right)}\left({t}_{2\left(1\right)}\right)\mathrm{sinh}\left(\sqrt{\left({k}_{1}{k}_{3}\right)/{m}_{2}}t\right)$
$+\frac{\left({m}_{2}\left({k}_{1}{l}_{AB}+{k}_{3}\left(2{l}_{AB}+{\delta}_{2}\right)\right){\omega}^{2}{k}_{3}\left({k}_{1}+{k}_{3}\right)\left(2{l}_{AB}+{\delta}_{2}\right)\right)\mathrm{cosh}\left(\sqrt{\frac{{k}_{1}{k}_{3}}{{m}_{2}}}t\right)}{\left({k}_{1}+{k}_{3}\right)\left({k}_{1}+{k}_{3}{m}_{2}{\omega}^{2}\right)},$
There are much more stages of the internal bodies oscillations. Due to the limited length of the present paper, let us pay major attention to the first two ones describing the starting conditions of the wheeled vibroimpact locomotion system.
3. Results and discussion
3.1. Stating and solving the parameters optimization problem
In the present research, let us consider the optimization problem consisting in maximizing the average translational velocity ${\dot{x}}_{1aver.}$ of the wheeled platform. The largest velocity can be reached in the case when the velocities of the internal bodies reach their maximal values during the contacts of the corresponding impact bodies and plates. Considering the stages of the internal bodies motions analyzed above, the optimization problem can be stated as follows:
The duration ${t}_{2\left(1\right)}$ of the first stage of the mass ${m}_{2}$ oscillations can be determined by substituting $x={x}_{2\left(1\right)}\left({t}_{2\left(1\right)}\right)=2{l}_{AB}+{\delta}_{2}$ into Eq. (11) and solving it with respect to time $t$. Then, by differentiating Eq. (11) with respect to time $t$, the velocity expression of the mass ${m}_{2}$ can be derived. Substituting the obtained value of the stage duration ${t}_{2\left(1\right)}$ into the derived expression ${\dot{x}}_{2\left(1\right)}$, the optimization function ${\dot{x}}_{2\left(1\right)}\left({l}_{AB},{k}_{1},{m}_{2},\omega ,{\delta}_{2}\right)$ can be deduced.
The duration ${t}_{3\left(1\right)}$ of the first stage of the mass ${m}_{3}$ oscillations can be determined after deriving the velocity expression by differentiating Eq. (12) with respect to time $t$ and substituting $\dot{x}={\dot{x}}_{3\left(1\right)}\left({t}_{3\left(1\right)}\right)=0$. The solution of the obtained equation with respect to time $t$ gives the corresponding duration ${t}_{3\left(1\right)}$. Substituting the obtained value ${t}_{3\left(1\right)}$ into Eq. (12), the maximal displacement ${x}_{3\left(1\right)}\left({t}_{3\left(1\right)}\right)$ can be determined. The duration ${t}_{3\left(2\right)}$ of the second stage of the mass ${m}_{3}$ oscillations can be determined by substituting $x={x}_{3\left(2\right)}\left({t}_{3\left(2\right)}\right)={x}_{3\left(1\right)}\left({t}_{3\left(1\right)}\right)+{\delta}_{3}$ into Eq. (14), and solving it with respect to time $t$. Then, by differentiating Eq. (14) with respect to time $t$, the velocity expression of the mass ${m}_{3}$ can be derived for the second stage of its oscillations. Substituting the obtained value of the stage duration ${t}_{3\left(2\right)}$ into the derived expression ${\dot{x}}_{3\left(2\right)}$, the optimization function ${\dot{x}}_{3\left(2\right)}\left({l}_{AD},{k}_{2},{m}_{3},\omega ,{\delta}_{3}\right)$ can be deduced.
The parameters optimization technique considered above has been programmed in the Mathematica software using the numerical maximization algorithm defining the global maximum of the functions ${\dot{x}}_{2\left(1\right)}\left({l}_{AB},{k}_{1},{m}_{2},\omega ,{\delta}_{2}\right)$, ${\dot{x}}_{3\left(2\right)}\left({l}_{AD},{k}_{2},{m}_{3},\omega ,{\delta}_{3}\right)$ subject to the given constraints: 1 kg $<{m}_{2}<$4 kg, 1 kg $<{m}_{3}<$4 kg, 10^{3}^{}N/m $<{k}_{1}<$ 10^{7}^{}N/m, 10^{3} N/m$<{k}_{2}<$10^{7}^{}N/m. Specifying the crank lengths ${l}_{AB}=\text{0.03m}$, ${l}_{AD}=\text{0.03m}$, the angular velocity $\omega =\text{314}\text{rad}/\text{s}$, and the impact gaps ${\delta}_{2}=\text{0.01m}$, ${\delta}_{3}=\text{0.01m}$, the optimal values of the masses of the internal bodies and the stiffness coefficients of the corresponding springs have been calculated: ${m}_{2}=\text{3.04kg}$, ${m}_{3}=\text{2.99kg}$, ${k}_{1}=$9.98∙10^{5} N/m, ${k}_{2}=$9.84∙10^{5}^{}N/m.
4. Numerical modeling of the system locomotion in Mathematica software
In order to perform further numerical modeling of the wheeled vibroimpact system locomotion in the Mathematica software, let us introduce the mass of the wheeled platform ${m}_{1}=$ 30 kg, and the stiffness coefficients of the impact springs ${k}_{3}=$ 10^{8}^{}N/m, ${k}_{4}=$ 10^{8} N/m.
The numerical modeling results of the system locomotion under the following controlled parameters ${\delta}_{2}={\delta}_{3}=$ 0.01 m (${\delta}_{2}^{*}=$ 0.07 m, ${\delta}_{3}^{*}=$ 0.01 m) are presented in Fig. 3. All the changes in the position of the wheeled platform are caused by the impacts (${x}_{2}>{\delta}_{2}^{*}$) of the mass ${m}_{2}$ at $t=$ 0.008 s, 0.027 s, 0.034 s, 0.053 s, and by the impacts (${x}_{3}>{\delta}_{3}^{*}$) of the mass ${m}_{3}$ at $t=$ 0.018 s, 0.037 s, 0.044 s, 0.056 s. These results allow for substantiating the correctness of the proposed mathematical model describing the system locomotion.
Fig. 2. Time dependences of the system’s bodies: a) displacements and b) velocities
a)
b)
In order to check the correctness of the optimization technique considered above, let us adopt the distances ${\delta}_{2}$, ${\delta}_{3}$ characterizing the corresponding impact gaps as the changeable parameters effecting the kinematic and dynamic characteristics of the vibroimpact system. Fig. 3 presents the dependencies of the wheeled platform displacements ${x}_{1}$ and average velocity ${\dot{x}}_{1aver.}$ on the impact gaps values ($0\le {\delta}_{2}={\delta}_{3}\le $ 0.03 m). The optimization problem defining the system’s inertiastiffness parameters (${m}_{2}$, ${m}_{3}$, ${k}_{1}$, ${k}_{2}$) was solved above for the values ${\delta}_{2}=$ 0.01 m, ${\delta}_{3}=$ 0.01 m. The obtained plots substantiate the maximal values of the wheeled platform average velocity of about 7 m/s in the impact gaps range 0.01 m $\le {\delta}_{2}={\delta}_{3}\le $ 0.015 m.
Fig. 3. Dependencies of the platform a) displacement and b) average velocity on the impact gaps values
a)
b)
5. Conclusions
The paper considers the improved design of the wheeled vibroimpact locomotion system (Fig. 1), which can be used for inspecting and cleaning the pipelines. It is equipped with the additional braking mechanisms that allows oneway rotation of the wheels. The vibrations are excited by the twin crankslider mechanism setting two internal bodies into the translational oscillatory motion. Providing the vibroimpact operational conditions, the stepwise (discontinuous) translational motion of the wheeled platform is studied. The corresponding mathematical model describing the system motion conditions is developed, and the analysis of the basic stages of the internal bodies oscillations is carried out. The analytical solutions of the derived differential equations are presented for the first two stages of the internal bodies oscillations. Based on the obtained solutions, the system’s parameters optimization problem is stated and solved using the Mathematica software. The numerical modeling of the system locomotion allows for substantiating the correctness of the developed mathematical model and of the proposed parameters optimization technique. All the changes in the position of the wheeled platform are caused by the impacts of the internal masses (Fig. 2), and the maximal values of the wheeled platform average velocity of about 7 m/s are located in the impact gaps range 0.01 m$\le {\delta}_{2}={\delta}_{3}\le $ 0.015 m (Fig. 3), which were initially used for solving the optimization problem. The increase in the impact gaps over 0.03 m causes the nonimpact operational conditions characterized by the lowest values of the wheeled platform average velocity.
References
 V.D. Nguyen, K.T. Ho, N.T. La, Q.H. Ngo, and K.T. Nguyen, “An experimental study on the selfpropelled locomotion system with anisotropic friction,” in Lecture Notes in Mechanical Engineering, Singapore: Springer Singapore, 2022, pp. 537–545, https://doi.org/10.1007/9789811632396_40 [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]
 Y. Yan, Y. Liu, J. Páez Chávez, F. Zonta, and A. Yusupov, “Proofofconcept prototype development of the selfpropelled capsule system for pipeline inspection,” Meccanica, Vol. 53, No. 8, pp. 1997–2012, Jun. 2018, https://doi.org/10.1007/s1101201708013 [Publisher]
 J. Xu and H. Fang, “Improving performance: recent progress on vibrationdriven locomotion systems,” Nonlinear Dynamics, Vol. 98, No. 4, pp. 2651–2669, Dec. 2019, https://doi.org/10.1007/s1107101904982y [Publisher]
 B. Diao, X. Zhang, H. Fang, and J. Xu, “Biobjective optimization for improving the locomotion performance of the vibrationdriven robot,” Archive of Applied Mechanics, Vol. 91, No. 5, pp. 2073–2088, May 2021, https://doi.org/10.1007/s00419020018705 [Publisher]
 I. A. Loukanov, V. G. Vitliemov, and I. V. Ivanov, “Dynamics of a vibrationdriven oneway moving wheeled robot,” IOSR Journal of Mechanical and Civil Engineering, Vol. 13, No. 3, pp. 14–22, 2016, https://doi.org/10.9790/16841303051422 [Search CrossRef]
 I. A. Loukanov, V. G. Vitliemov, and I. V. Ivanov, “Dynamics of a mobile mechanical system with vibration propulsion (VibroBot),” International Journal of Engineering Research and General Science, Vol. 4, No. 6, pp. 2320–9364, Jul. 2016. [Search CrossRef]
 I. A. Loukanov, “Multicriteria identification of VibroBot dynamic characteristics,” IOSR Journal of Engineering, Vol. 6, No. 7, pp. 26–35, Jul. 2016, https://doi.org/10.9790/3021067012635 [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, Cham: Springer International Publishing, 2021, pp. 43–52, https://doi.org/10.1007/9783030778231_5 [Publisher]
 V. Korendiy, V. Gursky, O. Kachur, V. Gurey, O. Havrylchenko, and O. Kotsiumbas, “Mathematical modeling of forced oscillations of semidefinite vibroimpact system sliding along rough horizontal surface,” Vibroengineering Procedia, Vol. 39, pp. 164–169, Dec. 2021, https://doi.org/10.21595/vp.2021.22298 [Publisher]
 V. Korendiy et al., “Kinematic and dynamic analysis of threemass oscillatory system of vibroimpact plate compactor with crank excitation mechanism,” Vibroengineering Procedia, Vol. 40, pp. 14–19, Feb. 2022, https://doi.org/10.21595/vp.2022.22393 [Publisher]
Cited By
Vitaliy Korendiy, Oleksandr Kachur, Volodymyr Gurskyi, Pavlo Krot 
2022

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

2022
