Mathematical Problems in Engineering
Volume 2017 (2017), Article ID 5676317, 11 pages
https://doi.org/10.1155/2017/5676317

Research Article

Stability Analysis of a Direct-Operated Seawater Hydraulic Relief Valve under Deep Sea

Shan Wu,1,2 Chao Li,² and Yipan Deng²

¹Department of Public Basic Course, Wuhan Technology and Business University, Wuhan, Hubei 430065, China
²School of Mechanical Science and Engineering, Huazhong University of Science and Technology, Wuhan 430074, China

Correspondence should be addressed to Chao Li; lichaobxd@hust.edu.cn

Received 16 February 2017; Revised 10 April 2017; Accepted 16 April 2017; Published 16 May 2017

Academic Editor: Oleg V. Gendelman

Copyright © 2017 Shan Wu et al. 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.

Abstract.

Taking into account the deformation of a designed direct-operated seawater hydraulic relief valve in deep sea, which might have a great influence on the stability of the valve, a mathematic model of the relief valve was established and stability analysis was conducted. As the fitting clearances between the damping sleeve and the damping bar play a key role in the performance of the relief valve, the fitting clearances after deformation under pressure of different ocean depths were obtained using finite element method. Applying the deformation data to the relief valve model, the stability and relative stability could be analyzed quantitatively through both the frequency domain analysis method and the time domain analysis method to detect the influence of the fitting clearance after deformation. The simulation results show that the seawater relief valve has a stable performance within 4000 meters deep under the sea.

1. Introduction

The traditional hydraulic technology takes mineral oil as its working medium. By contrast, the seawater hydraulic technology which takes seawater directly as a medium has a lot of advantages, including no pollution of leakage, simple structure, and easy maintenance, and thus it has a wide range of applications in deep-sea equipment. In recent years, seawater hydraulics have been successfully applied to walking/propulsion systems [1], buoyancy adjusting systems [2, 3], operating systems [4], and other key subsystems of marine equipment, and the researches focus mostly on matching materials [5–8] and water hydraulic components [9]. As a consequence, water hydraulics have become internationally the focus of attention in the field of fluid power transmission and control.

Almost every hydraulic system is equipped with a pressure relief valve to maintain the working pressure of the system at a predetermined level and protect the system from overpressure. There are two kinds of relief valves that are available: direct-operated (one-stage) [10] and pilot type [11]. Compared with the pilot-operated relief valve, the structure of the direct-operated type is simpler, and the operation is more stable and reliable [12]. Besides, for the direct-operated type, the leakage is less and the sealing is relatively easier. Therefore, such direct-operated relief valve using seawater as the working medium is frequently used [13]. However, the features of both physical and chemical performances of water are totally different from oil. The density of water is larger than that of oil. When the temperature is 20°C, the viscosity of seawater is 1.0 × 10⁻³ Pa·s, about 1/30–1/50 of the viscosity of commonly used hydraulic oil. The elastic modulus of water is about 1.2 times as much as that of oil [14, 15]. Due to these factors, the vibration and impact problem are highlighted for the water relief valve [16]. Adding a damping device is often used to reduce vibration and impact [10, 16, 17].

There are two common methods to analyze the dynamic characteristics of the hydraulic valves, that is, time domain analysis and frequency domain analysis. The results are more intuitionistic to use time domain analysis to simulate a hydraulic valve. Darby [18] and Hos et al. [19] conducted time response analysis for pressure relief valves in gas service, which can clearly show the dynamic response of a relief valve over time. In addition, as a method of frequency domain analysis, the bond diagram is usually used to analyze the characteristics of the hydraulic valves, such as Dasgupta’s research on a proportional solenoid controlled piloted relief valve [20] and Zanj’s research on an indirect hydrocontrol valve [21]. However, the deep-sea environment is greatly different from the environment of sea surface. It is characterized by high hydrostatic pressure (pressure increases by 1 atm for every 10 m in depth). This will make hydraulic components produce large deformation [22], which has a serious effect on the dynamic characteristics of the relief valve. In order to analyze the stability of the direct-operated seawater hydraulic relief valve under different ocean depths, this paper first uses a FEM (finite element method) technique to obtain the fitting clearance deformation under different hydrostatic pressures. Based on the deformation data, the state-space method and modern control theory, which can determine the stability of a linear system, are used to predict the stable operation of the direct-operated seawater hydraulic relief valve. In the end, the time response analysis is conducted to detect the dynamic performance of the relief valve.

Figure 1 shows the structure of the direct-operated seawater hydraulic relief valve studied in this paper. A damping cavity consists of a damping bar and a damping sleeve, which can minimize the vibration of plate spool. Its main structural parameters are shown in Table 1. While working in the submarine environment, the pressure of large ocean depth will lead to deformation of fitting clearances between the damping bar and the damping sleeve and will have an influence on the stability of the relief valve. In this paper, the 3D model of the valve is established to calculate the fitting clearance deformation under the ocean environment with FEM. On the basis of the deformation data, the state-space equation and TF (transfer function) are used to quantitatively analyze the stability and relative stability of the valve operating under different ocean depths.

Table 1: Main structure parameters of the relief valve.

Physical quantity
Sign
Value selection

Valve port flow coefficient

0.8

Seawater dynamic viscosity (Pa⋅s)

0.001

Elasticity modulus of seawater bulk (MPa)

1780

Diameter of front section of damping bar (mm)

10

Initial fitting clearance of front section of damping bar (mm)

0.02

Matched length of front section of damping bar (mm)

6.5

Initial length of damping cavity (mm)

3

Pressure value of working point (MPa)

14

Precompressed spring length (mm)

10.28

Pipeline volume from pump outlet to valve inlet (L)

2

Valve port speed coefficient

0.8

Seawater density (kg⋅m⁻³)

1025

Diameter of valve port (mm)

6

Contact surface diameter of back section of damping bar (mm)

6

Initial fitting clearance of back section of damping bar (mm)

0.02

Matched length of back section of damping bar (mm)

8

Equivalent mass of moving parts (g)

40

Stiffness coefficient of spring (N⋅mm⁻¹)

34.59

Working point spring compression variation (mm)

0.20

Pump outlet flow (L⋅min⁻¹)

30

Figure 1: Structure schematic diagram of direct-operated seawater relief valve.

2. Mathematical Model

2.1. Basic Model of Relief Valve

In order to simplify the analysis, the following assumptions are made:

(1)When the spool valve is moving, there are no considerations on friction and instantaneous hydraulic dynamics.(2)The influence of fluid mass is ignored.(3)The ambient pressure is considered as zero to analyze the performance. This is because the pressure compensated method is often taken in the deep-sea hydraulic system. When the relief valve is used in deep-sea equipment, the pressure difference of the valve between the input and the output is the same as that in land. Consequently, the force acting on the spool due to ambient pressure of the input and output is balanced. Thus, the ambient pressure does not have an influence on the analysis.(4)The gap width between the damping bar and the damping sleeve is uniform.(5)The transient flow force is generally disregarded in relief valves because of a much smaller part in various forces, and it is only considered in servo valves and proportional valves.(6), , and are taken as constant. This is because the stroke of the spool is relatively smaller than the matched length. Thus, the change of matched length will have a little influence on the performance of the valve.

Under the condition of the above assumptions, the mathematical model of seawater hydraulic relief valve (as shown in Figure 2) can be expressed by the following equations [10, 11, 16, 20, 23]:

Figure 2: Schematic diagram of direct-operated seawater relief valve with damping cavity.

(1) Flow equation of valve port:

where is the flow through the valve port and is the inlet pressure of valve port.

(2) Differential equation of valve spool movement:

where is the pressure of damping cavity; , , and are, respectively, the flow area of the valve seat, sectional area of the back end of the damping bar, and sectional area of the front end of the damping bar; and are the damping coefficients of movement, respectively, for the back and front end of the damping bar; that is, (); denotes the steady-state flow force generated by the change of speed and direction of the flow through the valve port, and [24].

(3) Continuity equation of valve port:

where is the pump output flow.

(4) Flow equation of damping cavity [25]:

where is the clearance flow of the back end of the damping bar and is the clearance flow of the front end of the damping bar.

(5) Continuity equation of damping cavity:

where is the initial volume of damping cavity; that is, .

2.2. Fitting Clearance Deformation of the Damping Bar under the Pressure of Large Depth

Under the high pressure of large depth, the fitting clearance between the damping bar and the damping sleeve will become smaller, which has a great influence on the dynamic characteristics of the seawater relief valve. To get fitting clearances at different ocean depths, structural static analyses of the damping sleeve and bar are made by ANSYS Workbench.

Figure 3 shows the finite element simulation model of the damping sleeve and bar, whose materials are, respectively, PEEK and ANSI 316 Stainless Steel. According to the actual operation condition (as shown in Figure 1), both ends of the damping sleeve and the bottom of the damping bar are fixed. These two parts are wholly steeped in seawater, so ocean depth pressure is added to all faces of the two parts. Then, the fitting clearances at different depths can be obtained.

(a) Simulation model of the damping sleeve

(b) Simulation model of the damping bar

Figure 3: Simulation model of fitting clearance deformation of the damping sleeve and bar.

Figure 4 shows the deformation cloud on the damping sleeve and bar at the depth of 4000 m, respectively. As can be seen from the picture, the deformation of the damping bar is relatively smaller than that of the damping sleeve. Thus, the deformation of the damping sleeve is the major part of the fitting clearance deformation.

(a) Deformation of the damping sleeve

(b) Deformation of the damping bar

Figure 4: Deformation of the damping sleeve and damping bar at the depth of 4000 m.

The predetermined clearances in land between the damping sleeve and both ends of the damping bar are  mm and  mm, respectively. And the variation of fitting clearances (as shown in Figure 5) of the front end and back end at different ocean depths can be obtained by finite element simulation. When the ocean depth is 4000 m, the fitting clearance of the damping sleeve and the front end of the damping bar decreases to 0 mm under huge seawater pressure, which means failure for the valve. Therefore, the discussion of the stability in this paper is all within the ocean depth of 4000 m.

(a) Fitting clearances of back end

(b) Fitting clearances of front end

Figure 5: Variation of damping fitting clearances at different ocean depths.

3. Analysis of Stability Characteristics

Taking the deformation at different ocean depths into account, qualitative and quantitative analysis of stability and relative stability for this relief valve will be carried out, respectively.

3.1. Stability Analysis

In order to facilitate analysis, the following variables are defined:

Defining , , , and and expanding the system around the normal working point , the model above can be expressed as

And then, to order , , a state-space equation of the relief valve can obtained:

where

The model is a system with two inputs and one output, and the stability of the system is merely related to matrix . Furthermore, by putting fixed parameters shown in Table 1 into formula (10) and taking the determinant of , the characteristic equation of the system can be obtained:

where

According to the Routh Stability Criterion, if all coefficients in the characteristic equation are positive and the signs of all elements in the first column of Routh Table are positive, the system is stable. With formulae (14)-(15), it can be known that every coefficient of the characteristic equations is positive if and . As a fourth-order equation for formula (14), the relief valve will be stable under the condition that signs of and in the first column of Routh Table are positive. To examine the stability of the relief valve at the ocean depth of 0–4000 m, we need to put fitting clearances at different ocean depths into formula (15). As shown in Figure 6, the coefficients in the first column of and in the Routh Table have changed over various ocean depths, and two coefficients within the ocean depth of 0–4000 m are all positive values. Therefore, it can be confirmed that this relief valve is stable at ocean depths of 0–4000 m.

(a) Coefficient value in the first column of

(b) Coefficient value in the first column of

Figure 6: Values of the first column coefficient in Routh Table with different ocean depths.

3.2. Relative Stability Analysis

The relative stability analysis gives the degree of stability of a stable system, which shows the margin when the system gets unstable. The relative stability analysis of the relief valve is also conducted.

In a practical system, a relief valve is often installed at the output of a pump, which means that the output flow of the pump is the input flow of the valve when all the flow of the pump is relieved by the relief valve. In formula (12), and are, respectively, seen as equivalent precompressed spring length and input flow of pump. Because the precompressed spring length is a definite value without affecting the stability of the system, the system can be simplified to a single-input one as follows:

The transfer function of the system can be obtained as follows:

According to the transfer function, we can obtain the BODE diagram of the relief valve in the initial conditions (i.e., 0 m ocean depth, fitting clearance without deformation) shown in Figure 7. It can be seen from the figure that the values of phase and magnitude margin of the relief valve are positive in initial conditions and also with ample margins, which shows that this valve operating on the surface of sea has a good relative stability.

Figure 7: BODE diagram of the relief valve in initial conditions.

Putting deformation data of the damping bar and damping sleeve into formula (17), the BODE diagram of the relief valve at different ocean depths can be obtained in Figure 8. From the figure, it can be seen that the peak gain of the relief valve gradually increases with greater ocean depth. Since the 0 Hz gains of the relief valve at different ocean depths are similar to each other, it can be concluded that the overshoot of the relief valve will increase with greater ocean depths.

(a) 400 m ocean depth

(b) 1300 m ocean depth

(c) 2500 m ocean depth

(d) 3900 m ocean depth

Figure 8: BODE diagrams of the relief valve under environment of different ocean depths.

Moreover, the value of phase and magnitude margin of the relief valve can be positive and simultaneously with ample margins at the ocean depths shown in Figure 8, so the fitting clearance changes of the damping bar and sleeve caused by huge pressure have a little influence on the relative stability of the relief valve.

4. Performances Simulation of Relief Valve under Different Ocean Depths

4.1. Dynamic Characteristics Simulation with Step Flow Input

On the basis of the structural parameters given above, the dynamic characteristic of the relief valve on the initial condition (i.e., 0 m ocean depth) can be obtained as shown in Figure 9. From the stimulation data, it can be calculated that the rise time  s, peak time  s, overshoot , setting time  s, and the value of control precision is 0.95.

Figure 9: Dynamic characteristics P-T curve of seawater relief valve in initial conditions.

With the ocean depth increasing, the fitting clearances of damping bar and sleeve will gradually decrease, and the changes of clearances are bound to have great influences on the dynamic characteristics of the valve. Putting the fitting clearances after deformation into the simulation model, the dynamic characteristics at different depths can be gotten as shown in Table 2. From the simulation results, the rise time of the valve will decrease slightly with the increase of ocean depths, but the overshoot will increase, which is up to 75.3% when the valve works at 2500 m ocean depth. Moreover, excessive amount of overshoot will increase the adjustment time of the relief valve to achieve stable pressure with the increase of ocean depth.

Table 2: Dynamic characteristics of seawater relief valve under different ocean depths.

Depth (m)
Rise time (s)
Overshoot (%)
Peak time (s)
Setting time (s)

0
0.0034
16.8%
0.0077
0.0130

400
0.0033
23.3%
0.0078
0.0141

800
0.0032
30.8%
0.0084
0.0228

1300
0.0031
41.2%
0.0083
0.0276

2500
0.00295
75.3%
0.01058
0.0500

4.2. Dynamic Characteristics Simulation with Pulsating Flow Input

The valve has been successfully applied to an underwater tool system. In this system, the valve is placed at the pump outlet. To detect the dynamic performance of the valve, the output flow of the pump is regarded as the input flow of the valve. Therefore, considering the flow ripple of the plunger pump, the real-time output flow of a plunger pump can be obtained as

where is the displacement of the plunger pump (40 mL/r); is the number of plungers of the pump (5); is the rated speed of the pump (750 r/min); is the central angle between two plungers of the pump, and .

Ignoring the influences of friction, corrosion, cavitations, strain, and other factors of the valve, the simulation results with design clearance are shown in Figure 10. Compared with Figure 9, it can be seen that the main dynamic characteristics of the system under pulsating flow input are in accordance with the performance under ideal input condition, while pressure response will be accompanied by a small pulsation at last. Fortunately, simulation results show that the valve has a good anti-interference ability.

Figure 10: Dynamic P-T curve of seawater relief valve under pulsating flow input (ocean depth of 0 m).

Figure 11 is the pressure response with pulsating flow at different ocean depths. Compared with the response with constant flow (shown in Figure 9), it can be figured out that the decrease of fitting clearance resulting from pressure of ocean depth does not reduce the anti-interference ability of the seawater relief valve. At the same time, the pressure pulsation of the valve caused by flow pulsation is similar to the performance of the initial clearance, without the increasing oscillation or divergence. Thus, simulation results give evidence to the stability of the relief valve.

(a) 400 m ocean depth

(b) 800 m ocean depth

(c) 1300 m ocean depth

(d) 2500 m ocean depth

Figure 11: Dynamic P-T curves of seawater relief valve at different ocean depths under pulsating flow input.

5. Conclusions

In this paper, a direct-operated seawater hydraulic relief valve with a damping device is introduced. The mathematical model of the valve is derived, in which the fitting clearance between the damping sleeve and damping bar () has a great effect on the dynamic stability. Considering the deformation of the damping sleeve and bar due to huge environmental pressure, the deformation simulation of the damping sleeve and bar at different ocean depths is achieved and different fitting clearances over varying depths are obtained. Based on the results obtained from the deformation simulation, the following conclusions can be made:

(1)Because of the deformation of the damping sleeve and bar, the relief valve has an effective working depth and it will fail to work at a greater depth. This is because the fitting clearance, that is, the clearance between the damping sleeve and bar, will decrease with the increase of the working depth. If the deformation is larger than the initial clearance, the moving parts of the valve will be stuck. Actually, this is a common problem for the components with precise clearance of matching pairs, such as piston/sleeve in a hydraulic pump, spool/sleeve in a hydraulic valve, and bearings in an electric motor, which must be considered in design.(2)For the relief valve above, the damping device is an important part and the clearance between the damping sleeve and bar is a depth sensitive parameter. When the depth is 4000 m, the clearance is 0 mm and the valve is invalid. Considering a certain safe margin, the working depth is 3500 m or less.(3)Although the fitting clearance becomes less with the increase of working depth, the relief valve is still stable according to the stability analysis.(4)The resonance peak of the relief valve in the relative stability analysis, which is positively correlated to the overshoot of the relief valve, increases with greater and greater ocean depth. And the relative stability analysis shows that relief valve always keeps an ample positive phase and magnitude margin within the depths of 0–4000 m.(5)The overshoot of the relief valve of step response increases with the increase of the ocean depth, and larger overshoot results in a longer setting time of the relief valve. Considering the 30% maximum overshoot of the criterion for land used relief valves, the maximum ocean depth of the relief valve suited for normal work is approximately 800 m. Irrespective of the 75.3% overshoot of the relief valve at 2500 m ocean depth, the relief valve still can achieve the desired stable pressure with a very short time.(6)It can be figured out that the decrease of fitting clearance resulting from the pressure of ocean depth does not reduce the anti-interference ability of the seawater relief valve. The pressure pulsation of the relief valve in deep sea, which is caused by the flow pulsation of piston pump, is similar to the situation with the initial clearance. And all the situations do not exhibit increasing oscillation or divergence.

The conclusions drawn in this paper not only offer guidelines for the selection or design of a direct-operated seawater hydraulic relief valve used in deep sea, but also provide the references for analyzing other types of hydraulic control components, such as flow pressure compensated control valves, pressure relief valves, or pilot relief valves, since those valves all can be abstracted as a coupling system of spring-mass and liquid accompanied by the stable problem. Moreover, this paper has also provided an analysis method, which puts deformation data of key fitting clearances of hydraulic components at different ocean depths into a state-space model for underwater hydraulic system.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the publication of this paper.

Acknowledgments

This paper is supported by the Program for Science Research from Education Department of Hubei Province (B2014168) and Wuhan Technology and Business University (B2013005), New Century Excellent Talents in University (NCET-13-0231), National Sci-Tech Support Plan (2014BAB12B00), and National Key R&D Plan (2016YFC0304800&2016YFC0300600).

References

B. Xin, L. Xiaohui, S. Zhaocun, and Z. Yuquan, “A vectored water jet propulsion method for autonomous underwater vehicles,” Ocean Engineering, vol. 74, pp. 133–140, 2013.

Y. Liu, X. Zhao, D. Wu, D. Li, and X. Li, “Study on the control methods of a water hydraulic variable ballast system for submersible vehicles,” Ocean Engineering, vol. 108, pp. 648–661, 2015.

X. Zhao, Y. Liu, M. Han, D. Wu, and D. Li, “Improving the performance of an AUV hovering system by introducing low-cost flow rate control into water hydraulic variable ballast system,” Ocean Engineering, vol. 125, pp. 155–169, 2016.

Y. Liu, D. Wu, D. Li, X. Zhao, and X. Li, “Applications of seawater hydraulics in deep-sea equipment,” Journal of Mechanical Engineering, vol. 50, no. 2, pp. 28–35, 2014.

D. Wu, Y. Liu, D. Li, X. Zhao, and Y. Liu, “Tribo-corrosion properties of WC-10Co-4Cr coating in natural silt-laden waters when sliding against Si3N4,” International Journal of Refractory Metals and Hard Materials, vol. 58, pp. 143–151, 2016.

D. Wu, Y. Liu, X. Zhao, D. Li, and X. Ren, “The tribological behaviors of different mass ratio Al2O3-TiO2 coatings in water lubrication sliding against Si3N4,” Tribology Transactions, vol. 59, no. 2, pp. 352–362, 2016.

D. Wu, Y. Liu, D. Li, X. Zhao, and X. Ren, “The applicability of WC-10Co-4Cr/Si3N4 tribopair to the different natural waters,” International Journal of Refractory Metals & Hard Materials, vol. 54, pp. 19–26, 2016.

D. Wu, Y. Liu, D. Li, and X. Zhao, “Effect of materials on the noise of a water hydraulic pump used in submersible,” Ocean Engineering, vol. 131, pp. 107–113, 2017.

J. Liang, X. Luo, Y. Liu, X. Li, and T. Shi, “A numerical investigation in effects of inlet pressure fluctuations on the flow and cavitation characteristics inside water hydraulic poppet valves,” International Journal of Heat and Mass Transfer, vol. 103, pp. 684–700, 2016.

K. Dasgupta and R. Karmakar, “Modelling and dynamics of single-stage pressure relief valve with directional damping,” Simulation Modelling Practice and Theory, vol. 10, no. 1-2, pp. 51–67, 2002.

K. Dasgupta and R. Karmakar, “Dynamic analysis of pilot operated pressure relief valve,” Simulation Modelling Practice and Theory, vol. 10, no. 1-2, pp. 35–49, 2002.

C. Bazsó and C. J. Hos, “An experimental study on the stability of a direct spring loaded poppet relief valve,” Journal of Fluids and Structures, vol. 42, pp. 456–465, 2013.

K. Suzuki and E. Urata, “Development of a direct pressure-sensing water hydraulic relief valve,” International Journal of Fluid Power, vol. 9, no. 2, pp. 5–13, 2008.

S. Hayashi, A. Shirai, N. N. Guo, and T. Hayase, “Static and dynamic characteristics of a pressure-compensated flow control,” in Proceedings of the 5th International Conference on Fluid Power Transmission and Control (ICFP '01), pp. 425–429, Hangzhou, China, 2001.

B. Wolfgang, “Water~or~oil hydraulics in the future,” in Proceedings of 6th Scandinavian International Conference on Fluid Power, pp. 51–64, Tampere, Finland, 1999.

R. D. Eyres, A. R. Champneys, and N. A. J. Lieven, “Modelling and dynamic response of a damper with relief valve,” Nonlinear Dynamics, vol. 40, no. 2, pp. 119–147, 2005.

X. Luo, X. He, S. Cao, and X. Ba, “Theoretical and experimental analysis of a one-stage water hydraulic relief valve with a one-way damper,” Journal of Pressure Vessel Technology, vol. 135, no. 6, Article ID 061210, 2013.

R. Darby, “The dynamic response of pressure relief valves in vapor or gas service, part I: mathematical model,” Journal of Loss Prevention in the Process Industries, vol. 26, no. 6, pp. 1262–1268, 2013.

C. J. Hos, A. R. Champneys, K. Paul, and M. McNeely, “Dynamic behaviour of direct spring loaded pressure relief valves in gas service: II reduced order modelling,” Journal of Loss Prevention in the Process Industries, vol. 36, pp. 1–12, 2015.

K. Dasgupta and J. Watton, “Dynamic analysis of proportional solenoid controlled piloted relief valve by bondgraph,” Simulation Modelling Practice and Theory, vol. 13, no. 1, pp. 21–38, 2005.

A. Zanj, H. Karimi, A. J. Gholi, and M. Shafiee, “Dynamic modeling of indirect hydro-control valve—Bondgraph approach,” Simulation Modelling Practice and Theory, vol. 28, pp. 65–80, 2012.

F. Yuan and X. Wu, “Hydrostatic pressure effects on deformation mechanisms of nanocrystalline fcc metals,” Computational Materials Science, vol. 85, pp. 8–15, 2014.

Q. Ye and J. Chen, “Dynamic analysis of a pilot-operated two-stage solenoid valve used in pneumatic system,” Simulation Modelling Practice and Theory, vol. 17, no. 5, pp. 794–816, 2009.

L. Yinshui, R. Xiaojun, W. Defa, L. Donglin, and L. Xiaohui, “Simulation and analysis of a seawater hydraulic relief valve in deep-sea environment,” Ocean Engineering, vol. 125, pp. 182–190, 2016.

G. E. Totten and V. J. De Negri, Handbook of Hydraulic Fluid Technology, CRC Press, 2nd edition, 2012.