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

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.


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][4][5][6][7][8][9][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.

Analytical dependencies for calculation
Consider a discrete-action -piston drive operating on a differential principle. The scheme of such a drive for 3 is shown in Fig. 1. Let's introduce the state space , , … , , that characterizes the position of the pneumatic drive pistons: where: 1 , if the -th cavity is connected to the high-pressure line, the -th piston is extended; 0, if the -th cavity is connected to the atmosphere, the -th piston is retracted. The variables 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 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 measured along the axis ( Fig. 1) is determined by the expression: , where: is the stroke of the -th piston. To describe the relative motion of pistons and piston groups, let's introduce the variables : where: 1, if the -th piston is moving relatively to the − 1 -th; 0 if the -th piston is stationary relatively to the − 1 -th.
The values of the variables 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: where: is the acceleration of the -th piston relatively to the − 1 -th; is the air pressure in the -th cavity; is the area of the -th piston from the -th cavity; is the friction force acting on the -th piston; is the mass of the -th piston; is the number of the last piston in the group of pistons moved by the -th piston.
The pressure changes in the cavities of the pneumatic motor are described by the following VIBROENGINEERING PROCEDIA. JUNE 2021, VOLUME 38 equation [11,12], modified for a multi-piston device: -atmosphere pressure. The correction factor of the air flow rate was determined by the dependence: where: = 1, 2; -system resistance coefficient, and the air flow function is: The friction force is defined as the sum of the viscous friction force and the friction in the seals: where: -viscosity index; the velocity of the -th piston in a fixed coordinate system;the diameter of a ring seal; -width of the sealing surface of the ring; -coefficient of the rubber friction on the cylinder wall; Δ -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 , were determined from the expressions: where: , the speeds of the advising pistons before impact. Speed of movement of the output rod of the pneumatic positioner is: and the current coordinate is: Combining the initial conditions and integrating the system of Eqs. (4,5,(8)(9)(10)(11)(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 Δ .
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.  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.

Discussion of the results
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 ⁄ = 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 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 = 600 and = 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 , 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 ⁄ = 0,75 for the moving masses 2,6 and 10 kg. The values of < 500 correspond to the unsteady movement during acceleration of the piston, but the time of movement of the pistons is minimal. At > 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).

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.