# Mathematical model of proportional spool valve.

1. IntroductionDifferent type membranate elements, such as pneumatic artificial muscles (PAM) and bellow cylinders, have a great potential to be used as power actuators in robotics and automation due to their advantages as lightness, easy replacement, safe operation, etc. [1]. For motion control of membranate actuator it is necessary to know not only mathematical model of such element, but also mathematical model of a pneumatic valve. There are number of papers to date in which were designed mathematical models of PAMs [2-7] and bellows [8-9]. However many of them include simplified mathematical model of the valve. Should be mentioned that there are some papers that present mathematical models of pneumatic proportional valves used for precise pressure control [10]. Pneumatic proportional spool valve (fig. 1) presented in this paper is used for fast and precise gas flow control. In combination with an external position and pressure sensors and controllers, a precise and fast-acting pneumatic system could be created.

In this paper is presented improved mathematical model of flow rate of pneumatic proportional valve Festo MPYE-5-1/8LF-010-B. Spool dynamics and electro-magnetic part are not included in this model (see in previous papers of authors [11-12]). In previous research of the authors were not included in the model relation with varying average temperature of gas flow through the valve as well as were given only general recommendations of choosing of friction factor coefficient.

2. Description of air flow through the valve

Pneumatic valve usually mathematically described as throttle mass flow rate through which one can calculate [13]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII];

where [G.sup.*]--means critical mass flow rate (Ma [approximately equal to] 1);

[p.sub.a]--means constant input pressure;

[p.sub.b]--means constant output pressure;

[mu]--means flow coefficient of throttle aperture;

b--means aperture width;

[[theta].sub.a]--means input gas temperature.

After air flow through the valve analysis, was developed equation system which describes mass flow rate through the valve [12]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII];

where [G.sub.1] and [G.sub.2]--means mass flow rates through first and second throttle apertures; [A.sub.ch]--means cross sectional area of the channel.

Valve inner space longitudinal section is shown in Figure 2. With numbers 1 and 2 are marked surfaces of air flow parameters (temperature, pressure, density, etc.) after and before corresponding throttle aperture. Then the solution of equation system (4) is the following:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

3. Estimation of head loss coefficient

Let begin with estimation of head loss coefficient [xi]. This coefficient consist of following three coefficients: coefficient of so-called "minor losses" because of air flow turn and collision after the first throttle aperture pair [[xi].sub.1], coefficient of friction in valve channel [[xi].sub.2] and coefficient of losses because of air flow turn and stopping before the second pair of throttle apertures [[xi].sub.3]. According to the parameters of channel length L and hydraulic diameter D air flow turn is about 45[degrees], and losses coefficients [[xi].sub.1] = [[xi].sub.3] = 1 [14]. The most complicate is to evaluate friction coefficient [[xi].sub.2] which could be calculated with Darcy-Weisbach equation [13]:

[xi] = 4f L/D

where f--means fanning friction factor;

L--means channel length;

D--means channel hydraulic diameter.

For the given valve channel length and hydraulic diameter rate is the following:

L/D = 10 mm / 4 mm = 2.5;

Let evaluate fanning friction factor. In general, this coefficient depends on channel parameters, flow character and air flow velocity. Should be mentioned that air flow in channel is turbulent because two air flows after passing radial located throttle apertures collide (in figure 3 shown turbulence kinetic energy).

Analyzing modelling results one can see that supersonic flow occurs mainly in downstream when small pressure ratio. Maximal air flow velocity in the valve channel is v = 170 m/s. Then Reynolds numbers for air flow in channel in different cases are in a range of 1000 to 45000. So far as air flow in the channel is turbulent and channel section is a ring, to calculate fanning friction factor one can use Blasius equation [13]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII];

Values of friction coefficient [[xi].sub.2] have the following limits:

0.054 [less than or equal to] [[xi].sub.2] [less than or equal to] 0.14'

Because of relatively small values of friction coefficient and to simplify calculations, one can take [[xi].sub.2] = 0.1. As a result, a value of head loss coefficient could be taken constant and [xi] = 2.1.

4. Evaluation of air flow temperature

The main complexity in air flow temperature evaluation is that for precise solution it is necessary to solve Navier-Stokes equations. At the same time, air flows through the valve pretty fast, therefore for calculation of temperature could be used equations described adiabatic process, but pressure values p1 and p2 still unknown. For this reason and to simplify calculation procedure will be use empirical equations.

In figure 4 are shown contours of temperature when different values of input and output pressure and opening of the valve, besides in case b) take place supercritical flow.

By modelling results was found out that average value of temperature [[theta].sub.1] mostly depends on input/output pressure ratio. This relation is not highly non-linear, so temperature value [[theta].sub.1] could be calculated with a help of the following linear empirical equation:

[[theta].sub.1] = [[theta].sub.i] (1 - 0.075 (1 - [p.sub.o]/[p.sub.i])) (2)

Average temperature value [[theta].sub.2] depends not only on input/output pressure ratio, but also on valve opening. There was developed the following equation to calculate temperature [[theta].sub.2] value:

[[theta].sub.2] = [[theta].sub.i] - 45 x [([p.sub.o]/[p.sub.i] - 1).sup.2] x/[x.sub.m] (3)

5. Results

In previous sections were evaluated values of average flow temperature [[theta].sub.1] and [[theta].sub.2] with equations (6)-(7) and head loss coefficient [xi]. All other constants from equation (5) were found by authors in previous research [10]. Values of mass flow rate through the valve received by modelling results and with an equation (5) with different values of input and output pressure are shown in figures 5 and 6.

Comparison of simulation and theoretical calculation shows us that the model works accurate (error is less than 6%), especially when pressure ratio is small. The main reason of deviation is probably the fact that there were developed simplified equations for flow average temperature and head loss coefficient because of very complex flow process in the valve channel. Nevertheless introduction of these simplified equations decrease model error from 10-15% to 6% in comparison with previous research.

As a result was developed accurate mathematical model of a proportional pneumatic spool valve. The main advantage of this model is that the model based not on experimental data of particular valve, but on general gas dynamics theory. Consequently this model could be used for a class of pneumatic valves.

Finally use of developed improved mathematical model of valve flow rate together with membranate actuators models makes possible a design of fast and precise robotics and automation systems.

6. Conclusion

Improved mathematical model of flow rate of proportional pneumatic spool valve was presented in this paper. The main idea of creation of such model was consideration of air flow temperature average values variations and head loss coefficients variations to decrease model error. Development of the model included:

* General description of the air flow through the valve

* Evaluation of head loss coefficient

* Evaluation of air flow temperature variations

* Verifying results using modelling

Emphasis was made on developing of accurate easy-to-calculate universal mathematical model which could be used for a class of pneumatic valves. This model of the valve could be used together with the model of the pneumatic muscle or bellow cylinder to simulate different types of closed-loop systems, for example position system or active damping system.

The logical next step would be verification of received results by experimental data. Further research will be focused on verifying received results with help of experiment.

DOI: 10.2507/26th.daaam.proceedings.085

7. References

[1]. Daerden, F., Lefeber, D.: Pneumatic Artificial Muscles: actuators for robotics and automation. European Journal of Mechanical and Environmental Engineering. Citeseer, 2002, issue 1, vol. 47, pp. 11 -21

[2]. Loshitskiy, P.A., Sharovatov, V.T., Mathematical model of dual-acting force mambranate rodeless pneumocylinder, Mechatronics, Automation, Control, Moscow, 2012, No. 4, pp. 24-30.

[3]. J. Pitel, R. Neydorf, J. Borzikova, Arm position simulation of PAM based actuator, Annals of DAAAM for 2011 & Proceedings of the 22nd International DAAAM Symposium, Katalinic, B. (ed.), DAAAM International Vienna, ISBN 978-3-901509-83-4, Vienna, 2011, pp. 0145-0146.

[4]. K. C. Wickramatunge, T. Leephakpreeda, Empirical modeling of dynamic behaviors of pneumatic artificial muscle actuators, ISA Transactions, Vol. 52, Issue 6, 2013, pp. 825-834.

[5]. S. Ganguly, A. Garg, A. Pasricha, S.K. Dwivedy, Control of pneumatic artificial muscle system through experimental modelling, Mechatronics, Vol. 22, Issue 8, 2012, pp. 1135-1147.

[6]. M. Tothova, A. Hosovsky, Dynamic Simulation Model of Pneumatic Actuator with Artificial Muscle, SAMI 2013, IEEE 11th International Symposium on Applied Machine Intelligence and Informatic, 2013, pp. 47-51.

[7]. M. Tothova, J. Pitel', Dynamic Simulation of Pneumatic Muscle Actuator in Matlab/Simulink Environment, SISY 2014, IEEE 12th International Symposium on Intelligent Systems and Informatics, 2014, pp. 209-213.

[8]. V. Sharovatov, P. Chernus, Mathematical Model of a Strength Part of a Push Type Single-Acting Rodless Pneumatic Cylinder, Mechatronics, Automation, Control, No. 9, 2014, pp. 30-36.

[9]. Chernus Pavel, Sharovatov Valery. Dynamic Mathematical Model of Two-way Bellow Actuator. // ProcediaEngineering, Vol. 100, 2015, pp. 1040-1045.

[10]. Varga, Z., Honkola, P.-K., Mathematical model of pneumatic proportional valve, Journal of applied science in thermodynamics and fluid mechanics Vol. 1, No. 1/2012.

[11]. P.P. Chernus, V.T. Sharovatov, Accounting for compressed gas properties in the mathematical model of a spool valve // Proceedings of higher educational institutions. Machine building, No. 12, 2014, pp. 22-29

[12]. Petr Chernus, Valery Sharovatov, Consideration of influence of peculiarities of compressed gas on mathematical model parameters of spool valve. // Procedia Engineering , Vol. 100, 2015, pp. 1046-1054.

[13]. Popov, D.N., Dynamics of hydraulic and pneumatic drives, Moscow, 2002.

[14]. Herz, E.B., Pneumatic devices and systems. Reference book, Moscow, 2010.

Petr Chernus, Valery Sharovatov, Pavel Chernus

Baltic state technical university "VOENMEH" named after D.F. Ustinov, 1st

Kranoarmeyskaya str. 1, 190005, St

-Petersburg, Russia

Caption: Fig. 1. Proportional pneumatic valve. 1--housing, 2--valve spool, 3 --housing for electronics

Caption: Fig. 2. Valve inner space longitudinal section

Caption: Fig. 3. Turbulence kinetic energy

Caption: Fig. 4. Air flow temperature. a) valve open on 40%, [p.sub.i] = 5 bar, [p.sub.o] = 1 bar; b) valve open on 80%, [p.sub.i] = 5 bar, [p.sub.o] = 1 bar; c) valve open on 60%, [p.sub.i] = 5 bar, [p.sub.o] = 3 bar; d) valve open on 100%, [p.sub.i] = 5 bar, [p.sub.o] = 3 bar.

Caption: Fig. 5. Mass flow rate when pi = 5 bar, po = 1 bar

Caption: Fig. 6. Mass flow rate when pi = 5 bar, po = 3 bar

Printer friendly Cite/link Email Feedback | |

Author: | Chernus, Petr; Sharovatov, Valery; Chernus, Pavel |
---|---|

Publication: | Annals of DAAAM & Proceedings |

Article Type: | Technical report |

Date: | Jan 1, 2015 |

Words: | 1864 |

Previous Article: | Modelling of two-way Bellow actuator positioning. |

Next Article: | Robot based machine hammer peening using an electromagnetic driven hammering device. |

Topics: |