Some recommendations for the calculation of pneumatic engines taking into account the smooth stop

Raevsky, Vladimir; Nasonov, Dmitry; Ilichev, Vladimir

doi:10.21595/vp.2021.22064

Vibroengineering Procedia

Browse Procedia

Published: 28 June 2021

Check for updates

Some recommendations for the calculation of pneumatic engines taking into account the smooth stop

Vladimir Raevsky¹

Dmitry Nasonov²

Vladimir Ilichev³

¹Department of Informatics and Information Technologies, Kaluga State University named after K.E. Tsiolkovski, Kaluga, Russian Federation

²Department of Machines’ Vibroacoustics, Mechanical Engineering Research Institute of the Russian Academy of Sciences, Moscow, Russian Federation

³Department of Heat Engines and Hydraulic Machines, Bauman State Technical University (Kaluga Branch), Kaluga, Russian Federation

Corresponding Author:

Dmitry Nasonov

Cite the article Download PDF

Downloads 689

Abstract

The purpose of this paper is to clarify the recommendations and to develop a methodology for calculating multi-piston pneumatic motors (pneumatic positioners) for analyzing the movement parameters of a pneumatic drive and taking into account the need for a smooth stop of the output rod. A system of differential equations is obtained, the solution of which gives a description of the dynamics processes occurring during the operation of a multi-piston pneumatic actuator. The software has been developed that allows to obtain a number of graphical dependencies, on the basis of which the designer decides on the expediency of using braking by the method of throttling the inputs and outputs of the pneumatic motor chambers.

1. Introduction

Pneumatic industrial robots (PIR) are still widely used to perform a wide range of lifting and transport operations, although they are currently being replaced by robots with an electric drive [1, 2]. It is believed that PIR are reliable and easy-to-operate, high-speed and low-cost robotic devices. PIR with a cyclic control system are characterized by high accuracy of the output to the desired point, but have a small number of positioning points. PIR with a positional system have a larger number of positioning points, but have less accuracy.

One of the ways to increase the number of positioning points in the cycle while maintaining high accuracy of the point output is the use of multi-piston pneumatic motors (pneumatic positioners). Multi-piston positioners are characterized by the need to dampen the movement of the output link of the drive at the end of the stroke (shock-free stop of the pistons at the end of the stroke), since at high speeds of movement of the output link when approaching the stops, strong impacts of the working body of the robot on the stops are possible. This worsens the dynamic characteristics of the drive, increases its operating time, and causes premature wear of the structure.

At the same time, the calculation of multi-piston pneumatic positioners has certain difficulties due to the fact that the internal pistons of the pneumatic motor do not have the properties of a rigid stop due to the pliability of the air, so the recommendations for the calculation of such pneumatic motors are of great practical importance.

In papers [3-10], methods for calculating robot pneumatic motors are given, in papers [11, 12] – methods for calculating the smooth stop of the piston at the end of the stroke for pneumatic motors with a single piston. However, the calculation of the smooth stop mode of a multi-piston engine has some features caused by the variable load on the pistons and many different combinations of movements of the piston groups, which cause difficulties in determining the parameters of the pneumatic drive.

2. Analytical dependencies for calculation

Consider a discrete-action $N$ -piston drive operating on a differential principle. The scheme of such a drive for $N =$ 3 is shown in Fig. 1.

Fig. 1N-discrete-action piston drive (N= 3)

Let’s introduce the state space $C_{1}, C_{2}, \dots, C_{N},$ that characterizes the position of the pneumatic drive pistons:

1

C_{i} = \{\begin{matrix} 1, \\ 0, \end{matrix}

where: $C_{i} = 1$ , if the $i$ -th cavity is connected to the high-pressure line, the $i$ -th piston is extended; $C_{i} = 0$ , if the $i$ -th cavity is connected to the atmosphere, the $i$ -th piston is retracted.

The variables $C_{i}$ of the state space establish a one-to-one correspondence between the state of the distributor valves and the position of the pneumatic positioner pistons for the time periods when the transients in the drive have ended. The values of the variables $C_{i}$ are assigned depending on the desired state of the motor output rod, which is conveniently set using binary code. Then the coordinate of the output rod $x_{ξ}$ measured along the $O X$ axis (Fig. 1) is determined by the expression:

2

x_{ξ} = \sum_{i = 1}^{N} C_{i} x_{0 i},

where: $x_{0 i}$ is the stroke of the $i$ -th piston.

To describe the relative motion of pistons and piston groups, let’s introduce the variables $u_{i}$ :

3

u_{i} = \{\begin{matrix} 1, \\ 0, \end{matrix}

where: $u_{i} = 1$ , if the $i$ -th piston is moving relatively to the $(i - 1)$ -th; $u_{i} = 0$ if the $i$ -th piston is stationary relatively to the $(i - 1)$ -th.

The values of the variables $u_{i}$ are determined by analyzing the forces acting on the pistons.

The system of equations describing the movements of the pneumatic drive pistons has the form:

4

{\ddot{x}}_{i} = u_{i} \frac{p_{i} S_{i} - p_{n} S_{n} - \sum_{j = 1}^{n} F_{t r j}}{\sum_{j = 1}^{n} m_{j}},

where: ${\ddot{x}}_{i}$ is the acceleration of the $i$ -th piston relatively to the $(i - 1)$ -th; $p_{i}$ is the air pressure in the $i$ -th cavity; $S_{i}$ is the area of the $i$ -th piston from the $i$ -th cavity; $F_{t r j}$ is the friction force acting on the $j$ -th piston; $m_{j}$ is the mass of the $j$ -th piston; $n$ is the number of the last piston in the group of pistons moved by the $i$ -th piston.

The pressure changes in the cavities of the pneumatic motor are described by the following equation [11, 12], modified for a multi-piston device:

5

{\dot{p}}_{i} = α_{i} ((C_{i} μ_{i 1} φ (σ_{1}) p_{M} - (1 - C_{i}) μ_{i 2} φ (σ_{2}) p_{i}^{β_{1}} p_{M}^{{- β}_{2}}) S_{i^{- 1}} + (1 - 2 C_{i}) k p_{i} x_{i}),

where: $α_{i} = \frac{k K \sqrt{R T_{M}} f_{i}}{x_{n i} + x_{i}}$ , $i = 1, 2, \dots, N$ ; $k =$ 1,4 – adiabatic index; $K = \sqrt{\frac{2 g k}{k - 1}} =$ 8,283 m^1/2/s at $g = 9,8$ m/s²; $R =$ 29,27 kgf∙m/ kgf °С – gas constant; $T_{M}$ – air temperature in the main line; $f_{i}$ –cross-sectional area of the $i$ -th entrance hole; $x_{n i} = \frac{V_{n i}}{S_{i}}$ – the given initial coordinate of the position of the $i$ -th piston at $i = 1, 2, \dots, N$ ; $x_{i}$ – the current position coordinate of the $i$ -th piston relatively to the $(i - 1)$ -th; $V_{n i}$ – initial volume of the $i$ -th cavity; $p_{M}$ – air pressure in the main line; $β_{1} = 3 k - 1 / 2 k$ ; $β_{2} = k - 1 / k$ ; $σ_{1} = p_{i} / p_{M}$ ; $σ_{2} = p_{A} / p_{i}$ ; $p_{A}$ – atmosphere pressure.

The correction factor of the air flow rate was determined by the dependence:

6

μ_{i j} = {φ (σ_{j})}^{- 1} \sqrt{\frac{k - 1}{2 k} \frac{1 - σ_{j}^{2}}{ε_{c} - 2 ln (σ_{j})},}

where: $j = 1, 2$ ; $ε_{c}$ – system resistance coefficient, and the air flow function is:

7

φ (σ_{j}) = \{\begin{matrix} \sqrt{{σ_{j}}^{\frac{2}{k}} - {σ_{j}}^{\frac{k + 1}{k}}}, {0,528 < σ}_{j} \leq 1, \\ 0,2588, {0 < σ}_{j} \leq 0,528 . \end{matrix}

The friction force is defined as the sum of the viscous friction force and the friction in the seals:

8

F_{t r j} = f_{v} {\dot{x}}_{A j} + π l D F_{t r} Δ p,

where: $f_{v}$ – viscosity index; ${\dot{x}}_{A j}$ the velocity of the $j$ -th piston in a fixed coordinate system; $D$ – the diameter of a ring seal; $l$ – width of the sealing surface of the ring; $F_{t r}$ – coefficient of the rubber friction on the cylinder wall; $Δ p$ – pressure drop between adjacent cavities.

The collisions of the pistons during the movement were considered as central impacts [13], the velocities of the pistons after the impact ${\dot{x}}_{i +}$ , ${\dot{x}}_{i + 1 +}$ were determined from the expressions:

9

{\dot{x}}_{i +} = \frac{(m_{i} - b m_{i + 1}) {\dot{x}}_{i -} - m_{i + 1} (1 + b) {\dot{x}}_{i + 1 -}}{m_{i} + m_{i + 1}},

10

{\dot{x}}_{i + 1 +} = \frac{m_{i} (1 + b) {\dot{x}}_{i -} + (m_{i + 1} - b m_{i}) {\dot{x}}_{i + 1 -}}{m_{i} + m_{i + 1}} .

where: ${\dot{x}}_{i -}$ , ${\dot{x}}_{i + 1 -}$ the speeds of the advising pistons before impact.

Speed of movement of the output rod of the pneumatic positioner is:

11

{\dot{x}}_{ξ} = \sum_{i = 1}^{n} {\dot{x}}_{i},

and the current coordinate is:

12

x_{ξ} = \sum_{i = 1}^{n} x_{i} .

Combining the initial conditions and integrating the system of Eqs. (4, 5, 8-12), taking into account Eqs. (6, 7), it is possible to obtain pictures of the processes occurring during the operation of the pneumatic motor when its output rod is extended to a different length, determined by the value of the strokes of the pistons. The integration of the system continues until all the pistons of the drive do not work out the specified positions. The position of the piston is considered exhausted if its current coordinate differs from the set one by an amount not exceeding the robot positioning error $Δ x$ .

In the environment of the Python 3 programming language (Scipy, NumPy, Matplotlib libraries), the software that implements the solution of the above dependencies has been created.

3. Discussion of the results

Fig. 2 shows the dependence of the speed of movement of 10 kg mass on time, obtained as a result of solving a system of equations in the developed software. The following parameters are set: piston diameter $D =$ 0,02 m, $D_{t r}$ = 0,01 m, line pressure $p_{M}$ = 4∙10⁵ Pa, $x_{01}$ = 0,016 m, $f_{1}$ = 0,002 m, system drag coefficient during acceleration $ε_{c}$ = 40, during braking $ε_{c}$ = 800. Fig. 3 shows the vibration of the piston after impact. Solid curves correspond to the mode of operation of the drive without braking, dotted ones-with braking by throttling the input when the pistons are extended and the output when they are retracted.

Fig. 210 kg mass movement speed

Fig. 3Piston vibrations during the 10 kg mass movement

As can be seen from the graphs shown in Fig. 3, the working time of the piston of the specified position of 0,016 m without braking is 0,103 s, with braking it is 0,195 s. However, in the first case, the time when the vibration attenuation to a value not exceeding the positioning error of 0,5 mm was 0,246 s, and in the second case it was 0. Thus, the value of the impact velocity of the piston at the end of its stroke significantly affects the time of working out the specified coordinate by the piston.

Calculations have shown that the main operating time of the drive is occupied by the movement of pistons with a long stroke length. This makes it possible to reduce the speed of movement of the pistons with a small stroke without increasing the operating time of the drive. So, with the area ratio $S_{N} / S_{N + 1} =$ 0,25 and movement of the 2 kg mass for the stroke length of 1, 2, 4 and 8 mm, the movement time was 0,04, 0.073, 0,15 and 0,3 s, respectively, while the speed of the piston at the end of the stroke did not exceed 0,04 m/s. To move the same mass to the length of the stroke of the piston 64 and 128 mm, a braking scheme with throttling of the inlet was used. The braking distance was set constant, equal to 3/4 of the stroke of the piston, the coefficients $ε_{i}$ were calculated, providing the necessary braking intensity and speed at the end of the stroke of the piston, which did not exceed 0,04 m/s. At a steady speed of 0,6 and 0,66 m/s, braking was provided at $ε_{c}$ = 600 and $ε_{c}$ = 400, respectively, and the travel time in this case was 0,22 and 1,2 seconds. Thus, the travel time of the pistons with a stroke of 1 and 2 mm can be neglected in comparison with the travel time of the piston with a stroke of 128 mm. The question of whether braking is necessary for pistons with a stroke length of 4 and 8 mm can be resolved depending on whether the overall speed of the drive is provided with an increase in the travel time of the pistons.

Fig. 4 shows the range of values of the drag coefficients $ε_{c}$ , at which the speed of the piston at the end of the stroke does not exceed 0,1 m/s with the ratio of the areas of the pistons $S_{N} / S_{N + 1} =$ 0,75 for the moving masses 2,6 and 10 kg. The values of $ε_{c} <$ 500 correspond to the unsteady movement during acceleration of the piston, but the time of movement of the pistons is minimal. At $ε_{c} >$ 500, the steady-state speed of the pistons does not exceed 0,1 m/s for different moving masses and different piston strokes (1,..., 128 mm).

Fig. 4The range of values of the resistance coefficients εс

Fig. 5Dependence of the time of movement of the piston with braking at the end of the stroke on the stroke length

Fig. 5 shows the dependence of the time of movement of the piston with braking at the end of the stroke of the piston on the length of its stroke. As the calculations show, with an increase in the speed of movement and the implementation of the braking mode, the time of movement of the piston weighing 2 kg for a stroke length of 16 mm decreased from 0,21 to 0,12 s, which indicates the feasibility of using the braking mode in this case.

4. Conclusions

The practical value of the work is as follows.

Analytical dependences are obtained that allow us to refine the dynamic characteristics of the movement of a pneumatic multi-piston drive, taking into account the need for a smooth stop of the output rod. The method is software-implemented, the correctness of the results obtained is confirmed by test calculations.

The novelty of the work is as follows.

As a result of research, it was found that when determining the operating time of the drive, it is necessary to take into account not only the time of movement of the pistons, but also the time of attenuation of the amplitude of their oscillations after stopping.

Reducing the speed of movement of the pistons with small strokes practically does not affect the operating time of the drive.

The using of braking modes for pistons with a long stroke length provides lower piston stroke speeds, while the drive operation time does not significantly increase.

The use of braking at the end of the stroke of the pneumatic drive pistons allows you to improve its dynamic characteristics without significant hardware costs.

References

Gerasimov I., Belyaev D., Kazachkov A. There is an opinion. Electric drive and robotics. Control Engineering Russia, Vol. 4, Issue 46, 2013, p. 18-19, (in Russian).

Search CrossRef
Sarvarov A., Vasiliev A., Danilenko K., Menshchikova E. Comparative analysis of mechatronic system drives. Electrotechnical Systems and Complexes, Vol. 4, Issue 25, 2014, p. 21-25, (in Russian).

Search CrossRef
Rachkov M. Optimal control of the position of a pneumatic manipulator with the use of an observation loop. Journal of Computer and Systems Sciences International, Vol. 40, Issue 3, 2001, p. 504-509.

Search CrossRef
Yu X. Study on servo control strategy for pneumatic manipulator. International Conference on Measuring Technology and Mechatronics Automation, 2010, p. 641-644.

Publisher
Shang C., Tao Meng G. D. Adaptive robust trajectory tracking control of a parallel manipulator driven by pneumatic cylinders. Advances in Mechanical Engineering, Vol. 8, Issue 4, 2016, p. 1-15.

Publisher
Yin Y., Wang J. A nonlinear feedback tracking control for pneumatic cylinders and experiment study. American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA, 2009, p. 3476-3481.

Search CrossRef
Wang J., Wang D. Modelling study, analysis and robust servo control of pneumatic cylinder actuator systems. IEE Proceedings. Control Theory and Applications, Vol. 148, Issue 1, 2001, p. 35-42.

Publisher
Meng D., Tao G., Zhu X. Adaptive robust motion trajectory tracking control of pneumatic cylinders. Journal of Central South University, Vol. 20, 2013, p. 3445-3460.

Publisher
Munasypov R., Tselishchev O. Electropneumatic tracking drive in the adaptive springing system of a mobile robotic complex. Bulletin of the Ufa State Aviation Technical University, Vol. 3, Issue 56, 2016, p. 194-200, (in Russian).

Search CrossRef
Gribkov A., Shilin D., Filin I. Test of kinematic model by pneumatic manipulator. Bulletin of the Samara State Aerospace University, Vol. 3, Issue 34, 2012, p. 203-208, (in Russian).

Search CrossRef
Hertz E. V., Kreinin G. V. Calculation of Pneumatic Actuators. Reference Manual. Mashinostroenie, Moscow, 1975, p. 272, (in Russian).

Search CrossRef
Belyaev N., Uvarov Stepanchuk E. E. Pneumohydraulic Systems. Calculation and Design. Higher School, Moscow, 1988, p. 271, (in Russian).

Search CrossRef
Panovko Ya Introduction to the Theory of Mechanical Shock. Nauka, Moscow, 1977, p. 224, (in Russian).

Search CrossRef

About this article

Received

17 May 2021

Accepted

09 June 2021

Published

28 June 2021

SUBJECTS

System dynamics in manufacturing system modeling

DOI

https://doi.org/10.21595/vp.2021.22064

Keywords

pneumatic drive

multi-piston actuator

smooth rod stop

positioning dynamics

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Previous article in issue Previous Next article in issue Next

Research article

2019 08 15

Multi-objective optimization of active suspension predictive control based on improved PSO algorithm

Qiang Zhao, Baoquan Zhu

Research article

2017 08 15

Active control for non-autonomous diaphragm-type pneumatic isolation system by using an augmented adaptive sliding-mode controller

Hung-Yi Chen, Jin-Wei Liang

Research article

2016 09 30

Improvement of ride performance with an active suspension based on fuzzy logic control

Yanghai Nan, Wei Shi, Pan Fang

Research article

2016 02 15

Sliding mode force tracking control for active hydro-pneumatic suspension

Jifu Guan, Dongdong Huo, Xinxia Ji

V. Raevsky, D. Nasonov, and V. Ilichev, “Some recommendations for the calculation of pneumatic engines taking into account the smooth stop,” Vibroengineering PROCEDIA, Vol. 38, pp. 166–171, Jun. 2021, https://doi.org/10.21595/vp.2021.22064

Copy Extrica

Copied to clipboard!

TY  - JOUR
DO  - 10.21595/vp.2021.22064
UR  - https://doi.org/10.21595/vp.2021.22064
TI  - Some recommendations for the calculation of pneumatic engines taking into account the smooth stop
T2  - Vibroengineering PROCEDIA
AU  - Nasonov, Dmitry
AU  - Raevsky, Vladimir
AU  - Ilichev, Vladimir
PY  - 2021
DA  - 2021/06/28
PB  - JVE International Ltd.
SP  - 166-171
VL  - 38
SN  - 2345-0533
SN  - 2538-8479
ER  - 

Copy Ris

Copied to clipboard!

@article{Nasonov_2021,
	doi = {10.21595/vp.2021.22064},
	url = {https://doi.org/10.21595/vp.2021.22064},
	year = 2021,
	month = {jun},
	publisher = {{JVE} International Ltd.},
	volume = {38},
	pages = {166--171},
	author = {Dmitry Nasonov and Vladimir Raevsky and Vladimir Ilichev},
	title = {Some recommendations for the calculation of pneumatic engines taking into account the smooth stop},
	journal = {Vibroengineering {PROCEDIA}}
}

Copy Bibtex

Copied to clipboard!

[1]D. Nasonov, V. Raevsky, and V. Ilichev, “Some recommendations for the calculation of pneumatic engines taking into account the smooth stop,” Vibroengineering PROCEDIA, vol. 38, pp. 166–171, Jun. 2021, doi: 10.21595/vp.2021.22064.

Copy IEEE

Copied to clipboard!

Nasonov, Dmitry, Vladimir Raevsky, and Vladimir Ilichev. “Some Recommendations for the Calculation of Pneumatic Engines Taking into Account the Smooth Stop.” Vibroengineering PROCEDIA 38 (June 28, 2021): 166–71. https://doi.org/10.21595/vp.2021.22064.

Copy Chicago

Copied to clipboard!

Some recommendations for the calculation of pneumatic engines taking into account the smooth stop

Abstract

1. Introduction

2. Analytical dependencies for calculation

3. Discussion of the results

4. Conclusions

References

About this article

Related Articles