Next Article in Journal
Energy and Exergy Analyses of Geothermal Organic Rankine Cycles Considering the Effect of Brine Reinjection Temperature
Previous Article in Journal
Modeling of a Hybrid Fuel Cell Powertrain with Power Split Logic for Onboard Energy Management Using a Longitudinal Dynamics Simulation Tool
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stator Curvature Optimization and Analysis of Axial Hydraulic Vane Pumps

School of Mechanical Power Engineering, Harbin University of Science and Technology, Harbin 150008, China
*
Author to whom correspondence should be addressed.
Energies 2022, 15(17), 6229; https://doi.org/10.3390/en15176229
Submission received: 30 July 2022 / Revised: 22 August 2022 / Accepted: 22 August 2022 / Published: 26 August 2022
(This article belongs to the Topic Exergy Analysis and Its Applications – 2nd Volume)

Abstract

:
Aiming at the problem of large vibration of a high-subside stator inner cavity curve vane pump, the force analysis of the vane at the transition curve is carried out, and the functional relationship between the vane turning angle θ, the large arc radius R, the small arc radius r and the sliding friction is established. The particle swarm algorithm is used to optimize the solution of the objective function, and the optimized parameter values are brought into the MATLAB simulation program to obtain the optimized stator curve profile diagram. The dynamic performance of the vane pump before and after optimization is simulated using ADAMS. The results show that: the acceleration of the vane pump slide is significantly reduced; the friction between the slide and the slide groove is significantly reduced; the contact force between the slide and the stator is significantly reduced; the impact vibration of the optimized vane pump is significantly reduced; and the dynamic performance of the vane pump is improved.

1. Introduction

A vane pump is a type of displacement pump that works by regularly changing the volume of a sealed working chamber. Compared to other structured pump types, the vane pump has the benefits of a compact structure, uniform flow rate and low rotor inert. It is widely employed in petrochemical, aerospace, metallurgical and other fields [1,2]. Figure 1 depicts the structure of a double-acting vane pump. The pump rotor is concentrically attached to the stator; the vanes are radially installed in the rotor slots; and the rotor glides back and forth in the places as it rotates with the main shaft. The cavity is formed by the rotor, vane, stator inner cavity surface and pump end face. In the first half of the rotation of the rotor, a vacuum is formed in the pump and the liquid is sucked in; when the volume of the cavity is reduced from the maximum value to the minimum value, the liquid medium is pressed out. The media transfer is completed by periodic adjustments in the volume cavity [3].
With the trend toward low noise and low pulsation in vane pumps, the classic vane pump’s long-term vibration, noise, durability and other lingering difficulties must be addressed [4], particularly noise and vibration. The pulsation of the pump output and pressure is propagated to all components of the hydraulic system, including the hydraulic hose and valves [5,6].
The noise is mainly created by the pressure pulsation, the stator curve structure and the working chamber pressure rise in three aspects [7], whereas the vibration is primarily caused by the flow pulsation caused by leakage and the impact vibration between the vane and the stator curve’s two elements. As a result, a proper design of the stator curve structure and the construction of the vane pump’s force model are critical to reducing noise and vibration and increasing the pump’s dynamic performance. The dynamic performance of the pump has been explored in a few related articles.
One of the origins of vibration in vane pumps is the imbalance forces, which increase the amplitude as the manufacturing errors rise [8]. Numerical approaches have been demonstrated to predict pump vibrations, particularly those generated by hydrodynamic forces. This technique, which may be used in the initial design phase [9], can reduce the computational effort and enhance development efficiency. Pump vibration can also be caused by the guiding vane of an axial flow pump. It was discovered that an appropriate number of guide vanes has a role in suppressing the pressure fluctuations and vibroacoustics of axial pumps by comparing the radiated noise field stimulated by diffusers with the varied numbers of vanes [10]. Experiments have also revealed that different vane leading edge forms alter vane pump vibrations, and a method for determining the degree of cavitation behavior has been developed. As a result, a measurement system with a pump and sensors was designed to successfully eliminate the cavitation phenomena when the pump is operated [11]. To forecast the vibration characteristics of a multistage rotor system, a novel axial dynamic model integrating the transient force of the balance disk was presented, which can serve as a reference for the axial vibration model of a multistage pump rotor system [12]. Meanwhile, a two-way iterative fluid–structure coupling-based technique is presented to analyze the vibration characteristics of axial flow pumps. The fluid–structure coupling interaction has received less attention. The apparent time-frequency law of fluid pressure pulsation and structure vibration is given [13].
Meanwhile, other researchers optimized pump vibration and noise, leading to developing a multi-objective experimental and response surface-based optimization strategy for improving centrifugal pump hydroacoustic performance. The ideal combination of parameters for the synergistic optimization of hydraulic performance and noise was established using a multiple regression response surface model, which might result in a 3.03 percent gain in hydraulic efficiency and a 3.01 percent reduction in noise [14]. The results show that the radiation noise has typical dipole characteristics, and its directionality varies with the flow rate. A hybrid numerical method was proposed to obtain the streamwise noise source for the radiation noise of centrifugal pumps; the finite element method was used to obtain the acoustic simulation results, and the numerical method was verified by experiments [15]. Simulating the force situation of the components under various dynamic conditions using computational fluid dynamics (CFD) software COMSOL, some scholars proposed a new type of cylindrical vane pump by deriving the displacement calculation formula based on the design structure of the cylindrical vane pump, deriving the theoretical stator internal surface equation and establishing a three-dimensional model and using CFD software COMSOL to simulate the force situation of the components under various dynamic conditions [16]. The candidate transition curve functions for the stator with triangular and polynomial forms have also been discussed, and the V-A-J method (which converts the shape optimization task into a parameter value optimization task) has been used to find the best candidate for the “velocity criterion” and “acceleration criterion” curve equations. Without considering the effect of fluid movement, the dynamic features of the transition curves were investigated, and the findings revealed that different stator curves are acceptable for various applications [17].
Although these studies aid in understanding the source of vibration and noise in the vane pump, there is less research in the form of a high-dimensional curve stator curve of the vane pump generated by vibration optimization. This work is based on a balanced eight-dimensional stator curve vane pump, the vane of the vane pump force analysis. Few people consider the effect of the ratio of large to small arc radius of the stator cavity on the performance of the transition curve. Especially in the process of constructing the derivation of higher order equations, the difference between the radius of large and small arcs directly affects the dynamic characteristics of the curve as the coefficient of the transition curve, which in turn affects the motion form of the vane, the contact state. Therefore, a reasonable stator curve is crucial to the performance of the vane pump. In this paper, we will combine the vane force analysis and consider the design requirements to optimize the large and small arc radius ratio of the stator inner cavity curve of double acting sliding vane pump by intelligent algorithm to get the optimal large and small arc radius ratio under the constraint condition.
The vibration noise during the operation of the vane pump and the stator inner cavity transition curve part of the wear are significant, as shown in Figure 2, and the size of the arc part has almost no scratch, according to the experiments. Each force on the slider is related to the octant of the stator curve, and the boundary condition of the octant is set with Rr and the angle δ of the transition curve section; in this vane pump design, R/r = 39/32 = 1.21875, and the stator inner cavity transition curve part of the wear is serious. The stator inner cavity length and short radius ratio must be optimized, as indicated by the red outlined area in Figure 2.

2. Establishment of the Blade Force Model’s Objective Function

Because there are apparent cut marks on the stator inner cavity surface transition curve of the pump after the experiment, consider the sliding friction between the vane and the wall as the objective function of optimization, and the value of R and r is the optimized result when the friction force is reduced to the minimum value.
The force analysis of the blade in the transition curve section is shown in Figure 3 [18].
In the transition curve section, the blades are exposed to centrifugal force Fi, spring force Ft, inertia force Fm, friction force Ff between the stator cavity and the stator, Coe-type inertia force Fk, viscous friction force Fη, contact reaction force Fn1, Fn2 of the rotor and hydraulic pressure Fp1, Fp2, Fp3, Fp4, Fp5. Each force’s angular function is illustrated below.
(1)
The centrifugal force of the blade as a function of the rotation angle [19] is shown in Equation (1).
F i ( θ ) = { m v ω 2 ( R λ C )         θ [ 0 , 20 ° ] & [ 160 ° , 200 ° ] & [ 340 ° , 360 ° ] m v ω 2 ( ρ ( θ ) λ C )         θ [ 20 ° , 70 ° ] & [ 110 ° , 160 ° ] & [ 200 ° , 250 ° ] & [ 290 ° , 340 ° ] m v ω 2 ( r λ C )         θ [ 70 ° , 110 ° ] & [ 250 ° , 290 ° ]
where λc is the blade center of mass to the stator cavity distance, mv is the rotating blade center mass, ω is the spindle angular velocity, R is the stator cavity large arc radius, r is the stator cavity small arc radius, ρ(θ) is the vector diameter of the blade.
(2)
The spring force on the blade as a function of the angle of rotation is shown in Equation (2).
F t ( θ ) = { k ( μ 15 )         θ [ 0 , 20 ° ] & [ 160 ° , 200 ° ] & [ 340 ° , 360 ° ] k ( μ ρ ( θ ) + 24 )         θ [ 20 ° , 70 ° ] & [ 110 ° , 160 ° ] & [ 200 ° , 250 ° ] & [ 290 ° , 340 ° ] k ( μ 8 )         θ [ 70 ° , 110 ° ] & [ 250 ° , 290 ° ]
where k is the spring stiffness factor, and μ is the spring free length.
(3)
The inertia force on the blade as a function of the angle of rotation is shown in Equation (3).
F m ( θ ) = { m ω 2 R = m π 2 n 2 900 R         θ [ 0 , 20 ° ] & [ 160 ° , 200 ° ] & [ 340 ° , 360 ° ] m π 2 n 2 900 a ( θ )         θ [ 20 ° , 70 ° ] & [ 110 ° , 160 ° ] & [ 200 ° , 250 ° ] & [ 290 ° , 340 ° ] m ω 2 r = m π 2 n 2 900 r         θ [ 70 ° , 110 ° ] & [ 250 ° , 290 ° ]
where m is the blade mass, n is the spindle speed, a(θ) is the radial acceleration of the blade in the transition curve section.
(4)
The frictional force between the blade and the inner cavity of the stator as a function of the rotation angle is shown in Equation (4).
F trip ( θ ) = { m ω 2 [ R L 2 ] · cos ( 90 ° α ) · f   θ [ 0 , 20 ° ] & [ 160 ° , 200 ° ] & [ 340 ° , 360 ° ] m ω 2 [ ρ ( θ ) L 2 ] · cos ( 90 ° α ) · f   θ [ 20 ° , 70 ° ] & [ 110 ° , 160 ° ] & [ 200 ° , 250 ° ] & [ 290 ° , 340 ° ] m ω 2 [ r L 2 ] · cos ( 90 ° α ) · f   θ [ 70 ° , 110 ° ] & [ 250 ° , 290 ° ]
where L is the blade length, α is the pressure angle formed by the blade and the stator cavity surface, and ƒ is the friction coefficient.
(5)
The Coe-type inertia force oonthe blade as a function of the rotation angle [20] is shown in Equation (5).
F k ( θ ) = { 2 m ω 2 [ R L 2 ]     θ [ 0 , 20 ° ] & [ 160 ° , 200 ° ] & [ 340 ° , 360 ° ] 2 m ω V ( θ ) · sin [ ω , V ( θ ) ]     θ [ 20 ° , 70 ° ] & [ 110 ° , 160 ° ] & [ 200 ° , 250 ° ] & [ 290 ° , 340 ° ] 2 m ω 2 [ r L 2 ]     θ [ 70 ° , 110 ° ] & [ 250 ° , 290 ° ]
where V(θ) is the radial velocity of the blade in the transition curve section.
(6)
The viscous friction force on the blade as a function of the rotation angle is shown in Equation (6).
F η ( θ ) = { η ω ( R L 2 ) [ A v r 1 h v r A v s h v s ]     θ [ 0 , 20 ° ] & [ 160 ° , 200 ° ] & [ 340 ° , 360 ° ] η V ( θ ) [ A v r h v r A v s h v s ]     θ [ 20 ° , 70 ° ] & [ 110 ° , 160 ° ] & [ 200 ° , 250 ° ] & [ 290 ° , 340 ° ] η ω ( r L 2 ) [ A v r 2 h v r A v s h v s ]     θ [ 70 ° , 110 ° ] & [ 250 ° , 290 ° ]
where η is the oil viscosity, Ab is the actual contact area between the bottom of the blade and the oil according to the design size of the blade, Avs is the actual contact area between the blade and the end plate, Avr1 is the actual contact area between the large circular arc section blade and the rotor, Avr2 is the actual contact area between the small circular arc section blade and the rotor, the actual contact area between the blade and the rotor in the stator transition curve section A v r = 65 [ 46 ρ ( θ ) ] . hvr is the clearance between the blade and the rotor, hvs is the clearance between the blade and the end plate.
(7)
The contact reaction force oonthe blade by the rotor is shown in Equation (7) as a function of the rotation angle.
F N ( θ ) = { F N 1 = ξ · L 1 ( f cos α sin α ) + ( L 1 L / 2 ) F k ( θ ) L L 1 F N 2 = ξ · L ( f cos α sin α ) + L / 2 · F k ( θ ) L L 1 θ [ 0 , 360 ° ]
where L1 is the length of the blade extending outside the rotor slot, L is the length of the blade, and ξ is the coefficient of the contact reaction force on the blade by the rotor.
(8)
The hydraulic pressure on the blade as a function of the turning angle is shown in Equation (8).
F p ( θ ) = { F p 1 = a 1 · sin ( b 1 θ + C 1 ) +                                     a 2 · sin ( b 2 θ + C 2 ) + a 3 · sin ( b 3 θ C 3 ) F p 2 = a 4 · sin ( b 4 θ C 4 ) +                                     a 5 · sin ( b 5 θ + C 5 ) + a 6 · sin ( b 6 θ + C 6 ) F p 3 = a 7 · sin ( b 7 θ C 7 ) +                                     a 8 · sin ( b 8 θ C 8 ) + a 9 · sin ( b 9 θ + C 9 ) F p 4 = a 10 · sin ( b 10 θ + C 10 ) +                                     a 11 · sin ( b 11 θ + C 11 ) + a 12 · sin ( b 12 θ + C 12 ) F p 5 = 41840158.12 N θ [ 20 ° , 70 ° ] & [ 110 ° , 160 ° ] & [ 200 ° , 250 ° ] & [ 290 ° , 340 ° ]
where the blade is subjected to the various hydraulic component forces by design of experiment (DOE) simulation results for fitting the way to establish the individual component forces as follows. The function equation of the hydraulic pressure on the slide is detailed in Appendix A.
The blade rotation angle as a function of the hydraulic pressure at the wall of the pressure chamber is shown in Equation (9).
F p 1 ( θ ) = a 1 · sin ( b 1 · θ + c 1 ) + a 2 · sin ( b 2 · θ + c 2 ) + a 3 · sin ( b 3 · θ + c 3 )
a 1 = 9.277 × 10 5   ( 4.131 × 10 8 , 4.15 × 10 8 ) b 1 = 0.0597   ( 8.295 , 8.415 ) c 1 = 0.255   ( 403.3 , 403.8 ) a 2 = 5.412 × 10 5   ( 4.329 × 10 8 , 4.34 × 10 8 ) b 2 = 0.088   ( 11.26 , 11.43 ) c 2 = 2.063   ( 482.6 , 486.7 ) a 3 = 4.692 × 10 4   ( 1.088 × 10 6 , 1.182 × 10 6 ) b 3 = 0.2742   ( 1.396 , 1.944 ) c 3 = 1.033   ( 52.72 , 50.65 )
The blade rotation angle as a function of the hydraulic pressure on the circular surface of the oil pressure chamber is shown in Equation (10).
F p 2 ( θ ) = a 4 · sin ( b 4 · θ + c 4 ) + a 5 · sin ( b 5 · θ + c 5 ) + a 6 · sin ( b 6 · θ + c 6 )
a 4 = 4.109 × 10 6   ( 1.808 × 10 11 , 1.808 × 10 8 ) b 4 = 0.1298   ( 443.2 , 443.4 ) c 4 = 1.872   ( 1.493 × 10 4 , 1.492 × 10 4 ) a 5 = 9.939 × 10 6   ( 4.512 × 10 11 , 4.512 × 10 11 ) b 5 = 0.1598   ( 1068 , 1068 ) c 5 = 0.2701   ( 3.534 × 10 4 , 3.534 × 10 4 ) a 6 = 6.122 × 10 6   ( 6.32 × 10 11 , 6.32 × 10 11 ) b 6 = 0.177   ( 625.6 , 626 ) c 6 = 2.845   ( 2.063 × 10 4 , 2.063 × 10 4 )
The blade rotation angle as a function of the liquid pressure at the wall of the suction chamber is shown in Equation (11).
F p 3 ( θ ) = a 7 · sin ( b 7 · θ + c 7 ) + a 8 · sin ( b 8 · θ + c 8 ) + a 9 · sin ( b 9 · θ + c 9 )
a 7 = 1.615 × 10 6   ( 1.544 × 10 10 , 1.545 × 10 10 ) b 7 = 0.1021   ( 85.82 , 86.03 ) c 7 = 1.336   ( 3091 , 3089 ) a 8 = 3.58 × 10 6   ( 1.156 × 10 11 , 1.156 × 10 11 ) b 8 = 0.131   ( 370.9 , 371.2 ) c 8 = 0.8403   ( 1.161 × 10 4 , 1.161 × 10 4 ) a 9 = 2.279 × 10 6   ( 1.311 × 10 11 , 1.311 × 10 11 ) b 9 = 0.1435   ( 274.4 , 274.7 ) c 9 = 3.598   ( 8224 , 8231 )
The blade rotation angle as a function of the fluid pressure at the circular arc of the suction chamber is shown in Equation (12).
F p 4 ( θ ) = a 10 · sin ( b 10 · θ + c 10 ) + a 11 · sin ( b 11 · θ + c 11 ) + a 12 · sin ( b 12 · θ + c 12 )
a 10 = 7.522 × 10 5   ( 2.657 × 10 5 , 1.239 × 10 6 ) b 10 = 0.03875   ( 0.002936 , 0.08044 ) c 10 = 0.3424   ( 1.531 , 2.215 ) a 11 = 1.541 × 10 5   ( 1.575 × 10 6 , 1.883 × 10 6 ) b 11 = 0.219   ( 0.8047 , 1.243 ) c 11 = 0.309   ( 29.13 , 29.75 ) a 12 = 2.349 × 10 5   ( 1.531 × 10 6 , 2.001 × 10 6 ) b 12 = 0.3297   ( 0.2796 , 0.9391 ) c 12 = 0.2197   ( 18.26 , 18.7 )
The contact angle and rotation angle function equations are as follows. MATLAB was used to fit the function equation of the blade contact angle α and the rotation angle θ by sampling ten coordinate points [21]. After comparing various types of suitable approximation methods, it was discovered that the Gaussian approximation method was the closest to the resultant data, so the function equation of the blade contact angle and the rotation angle was fitted by this method as follows. The functional relationship between the contact angle of the arc top profile of the sliding plate and the rotation angle of the sliding plate is described in detail in Appendix B.
The Gaussian approximation method generic function model is shown in Equation (13).
f ( x ) = a 1 · exp ( ( ( x b 1 ) / c 1 ) 2 )
The coefficient (95% confidence level) is as follows.
a 1 = 48.11           ( 46.15 , 50.07 ) b 1 = 27.1                 ( 26.64 , 27.57 ) c 1 = 13.89             ( 13.27 , 14.55 )
The results are showed in Equation (14).
α = 48.11 · exp ( ( ( θ 27.1 ) / 13.89 ) 2 )
To make it easier to develop the equation of the force on the slider as a function of the slider rotation angle, the relationship between the slider rotation angle θ and the contact angle α of the stator is determined.
As the blade and stator cavity mutual wear situation is serious, fully consider the blade force situation, the blade and stator inner wall between the minimum value of the sliding friction as the optimization objective function. Additionally, the calculation formula of sliding friction is F f = f · F n , according to Figure 4 force analysis, to obtain the force balance sketch in the direction of the common normal, combined with the actual force situation of the vane pump to obtain the blade by each force, according to the sliding friction formula to find the reverse support force Fn on the wall surface of the blade.
The equation of force balance in the direction of the common normal is shown in Equation (15).
F n + F η · cos α + ( F n 1 + F k + F p 3 ) · cos ( 90 ° α ) ( F n 2 + F p 1 ) · cos ( 90 ° α )                 ( F p 5 + F k + F i ) · cos α + F p 4 · cos ( 90 ° α 2 ) + F p 2 · cos ( 90 ° + α 2 ) = 0
The collation yields Equation (16).
F n = ( F p 5 + F k + F i F η ) · cos α + ( F n 2 + F p 1 F p 3 F n 1 F k ) · cos ( 90 ° α )                                       F p 4 · cos ( 90 ° α 2 ) F p 2 · cos ( 90 ° + α 2 )
The resulting sliding friction relationship between the slider and the stator cavity is shown in Equations (17) and (18).
F f = f · F n = f · [ ( F p 5 + F k + F i F η ) · cos α + ( F n 2 + F p 1 F p 3 F n 1 F k ) · cos ( 90 ° α ) F p 4 · cos ( 90 ° α 2 ) F p 2 · cos ( 90 ° + α 2 ) ]
F f = f · F n = f · [ ( F p 5 + k · ( μ ρ ( θ ) + 24 ) + m · ω 2 ( ρ ( θ ) λ c ) η V ( θ ) [ A v r h v r A v s h v s ] ) · cos α + cos ( 90 ° α ) · ( ξ · [ L · ( f · cos α sin α ) + ( L / 2 ) · F k ( θ ) ( L 1 · ( f · cos α sin α ) + ( L 1 L / 2 ) · F k ( θ ) ) ] L L 2 + F p 1 F p 3 2 m ω · V ( θ ) ) F p 4 · cos ( 90 ° α 2 ) F p 2 · cos ( 90 ° + α 2 ) ]
which A v r = 65 [ 46 ρ ( θ ) ] and ξ = ( F p + F i F m ) ( L L 1 ) L 1 f F k ( θ ) [ ( f 2 + 1 ) L + ( f 2 1 ) L 1 ] cos α 2 L 1 f fsin α .
The remaining values of each parameter are shown in Table 1.
The comprehensive values above can be obtained and will be the slide and the stator cavity of the minimum value of the sliding friction of the objective function relationship, as shown in Equation (19).
min F f = f · F n = f · [ ( F p 5 + k · ( μ ρ ( θ ) + 24 ) + m · ω 2 ( ρ ( θ ) λ c ) η V ( θ ) [ A v r h v r A v s h v s ] ) · cos α + cos ( 90 ° α ) · ( ξ · [ L · ( f · cos α sin α ) + ( L / 2 ) · F k ( θ ) ( L 1 · ( f · cos α sin α ) + ( L 1 L / 2 ) · F k ( θ ) ) ] L L 2 + F p 1 F p 3 2 m ω · V ( θ ) ) F p 4 · cos ( 90 ° α 2 ) F p 2 · cos ( 90 ° + α 2 ) ]

3. Stator Inner Cavity Curve Length and Radius Optimization

Eberhart and Kennedy proposed that the particle swarm optimization (PSO) algorithm consists of a series of iterative forms shown in Equations (20) and (21), iterating the velocity magnitude and direction of travel at times until the termination condition is reached [22].
V i , j ( t + 1 ) = ω V i , j ( t ) + C 1 R   and   [ P i , j x i , j ( t ) ] + C 2 R   and   [ P i , j x i , j ( t ) ]
x i , j ( t + 1 ) = x i , j ( t ) + V i , j ( t + 1 )       ( i , j   =   1 ,   2 ,   ,   d )
where ω is the inertia weight factor, C1 and C2 are learning elements with values in the range between 0 and 2, and rand and Rand are two mutually independent random functions following U(0,1). C1 is the particle’s distance from its optimal position step, C2 is the particle’s distance from the global optimal position step, and C1 = C2 = 1.49445 is taken in the requested objective function.
Constraints on the optimization of the stator length and radius parameters. The objective function of minimizing the sliding friction between the slider and the stator cavity is solved by the intelligent particle swarm algorithm, and some constraints still need to be set to make the solution converge more easily and find the optimal solution according to the design requirements and size limitations [23]. The constraints are shown in Table 2.
We solve the objective function with constraints according to the particle swarm algorithm.
(1)
When Rr = 1, the minimum value of the objective function is 4.8298 × 106, at which time R = 33, r = 32, θ = 26.65°.
(2)
When Rr = 2, the minimum value of the objective function is 4.4148 × 106, when R = 34, r = 32, θ = 24.56°.
(3)
When Rr = 3, the minimum value of the objective function is 4.0240 × 106, when R = 35, r = 32, θ = 23.71°.
(4)
When Rr = 4, the minimum value of the objective function is 3.6583 × 106, when R = 36, r = 32, θ = 22.90°.
(5)
When Rr = 5, the minimum value of the objective function is 3.3217 × 106, when R = 37, r = 32, θ = 23.13°.
(6)
When Rr = 6, the minimum value of the objective function is 3.4313 × 106, when R = 38, r = 32, θ = 22.14°.
The objective function value of the sliding friction between the slide and the stator cavity is minimized for R = 37, r = 32 and θ = 23.13°, as calculated by the particle swarm algorithm analysis. The simulation results are shown in Figure 5 and Figure 6.
The optimized results R = 37, r = 32 are brought into the octant curve simulation program. The comparison results before and after the optimization of the stator cavity transition curve can be obtained by the simulation program, as shown in Figure 7.

4. Results

Based on the stator model before and after optimization, the ADAMS dynamics simulation of the virtual prototype before and after stator optimization is carried out to predict the effect of the optimized stator on pump kinematics and dynamics performance.
We created the simulation model in ADAMS/View, set up the working environment and set up the gravity, coordinate system, etc., as shown in Figure 8.
Spring force is applied between the blade and the rotor blade slot, friction force between the edge and the inner wall of the rotor blade slot. Solid-to-solid contact force is applied between the blade and the stator cavity and a force acting on the bottom and both sides of the edge, applying a given speed to the spindle, setting the simulated spindle speed to 300 rpm.
Sliding vane pump simulation results are compared and analyzed as follows.
(1)
Spindle torque curve before and after stator cavity length to radius ratio optimization.
Comparing Figure 9 and Figure 10, it can be seen from the simulation curve that there is a prominent peak in the shaft torque before the optimization of the central stator cavity. This analysis is put forward because when the difference of Rr is significant, a curve section with large curvature will be generated at the transition curve. Therefore, an ample torque is developed in this curve section when the main shaft rotates. The optimized simulation results show that the maximum value of the spindle torque is controlled within 5800 N·mm. Compared with the maximum torque peak of 16,000 N·mm before optimization, there is no prominent peak of the spindle torque in the whole rotation cycle. The optimization of the length and radius of the stator cavity has a good effect on reducing the spindle torque.
(2)
Slider acceleration curves before and after optimizing stator cavity length to radius ratio.
As shown in Figure 11 and Figure 12, the simulation results before the optimization of the stator cavity show that the acceleration of the slide has a sudden change of the great peak; the reason is that when the difference of Rr is significant, it will produce a curve section with large curvature at the transition curve. Therefore, when the slide slides on the stator transition curved surface, it will significantly impact this curve surface, resulting in severe wear and vibration. The simulation results of stator cavity optimization show that the maximum value of slide acceleration is controlled within 3.0 × 105 mm/s2 compared with the peak value of 5.0 × 106 mm/s2 before optimization; there is no significant peak value of slide acceleration in the whole rotation cycle, which shows that the optimization of the length and radius of stator cavity has a more significant effect on reducing the slide acceleration impact.
(3)
Friction curves between the slider and rotor slider slot before and after optimizing the stator cavity length and short radius ratio.
As shown in Figure 13 and Figure 14, the simulation results before optimizing the stator cavity show that the friction between the slide and the rotor slide slot changed abruptly to a great peak. The reason is that when the difference of Rr is significant, a curve section with large curvature will be generated at the transition curve. Therefore, when the slider slides on the stator transition curved surface, it will generate a large reaction force on the slider arc side, which will increase the pressure between the slider and the rotor slider slot wall, and the sliding friction will increase when the slider slides radially, resulting in deeper scratches. From the stator cavity optimization simulation results, it can be seen that the friction force between the slider and the rotor slider slot is controlled within 0.0032 N. Compared with the peak friction force of 50 N between the slider and the rotor slider slot before optimization, there is no large peak friction force between the slider and the rotor slider slot during the whole rotation cycle.
(4)
Slider–stator contact force curve before and after optimizing the stator cavity length to short radius ratio.
As shown in Figure 15 and Figure 16, the simulation results before stator cavity optimization show that the contact force between the slide and the stator has a sudden and extreme peak because when the difference of Rr is significant, a curve section with large curvature will be generated at the transition curve. As a result, when the slide slides on the stator transition curved surface, it will produce an excellent acceleration impact on this curved surface. Then, the slide and the stator will have a great elastic deformation, resulting in a sudden change in the contact force and causing a severe wear phenomenon. The simulation results of the stator cavity optimization show that the maximum value of the contact force between the slider and the stator is controlled within 21 N. Compared with the peak value of the contact force between the slider and the stator before optimization, which reaches 75 N, the peak value of the contact force between the slider and the stator is significantly reduced during the whole rotation cycle, which indicates that the optimization of the length and radius of the stator cavity has a beneficial effect on reducing the contact force between the slider and the stator.
From the ADAMS simulation results, it can be seen that the spindle torque, slider acceleration, blade and rotor friction, and blade and stator contact force numerical magnitude are reduced before and after stator optimization, indicating that the optimization of the length and radius of the stator inner cavity has a beneficial effect on lowering the blade and stator contact force.

5. Conclusions

With algorithm optimization and dynamics simulation of the transition curve of the stator cavity of the vane pump, the following research results were obtained.
  • The explicit expression of the blade contact angle as a function of the blade rotation angle was obtained by fitting the values of the sampling points. The analysis showed that the blade contact angle was maximum in the middle of the transition curve and that the Gaussian approximation method had a good approximation effect.
  • The optimization model of the long–short radius ratio using an intelligent particle swarm algorithm. According to the force balance equation, the expression of the sliding friction between the blade and the stator was derived. To minimize the sliding friction force on the edge, the optimized solution of the long–short radius ratio of the inner cavity of the stator was calculated.
  • The dynamics of the vane pump was simulated using Adams separately, and the main shaft torque, vane acceleration, vane and rotor friction, and vane and stator contact force before and after optimization were compared and analyzed. The results showed that the optimized stator performance was better.

Author Contributions

Conceptualization, D.X.; methodology, D.X.; software, S.L.; validation, D.X., J.W. and X.B.; formal analysis, D.X.; investigation, S.L.; data curation, D.X.; writing—original draft preparation, D.X.; writing—review and editing, D.X. and Y.S.; visualization, J.W. and X.B.; supervision, Y.S.; project administration, Y.S.; funding acquisition, Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

The research is supported by Heilongjiang Province to support the transformation of major scientific and technological achievements project (Grant No. CG20A006); National Defense Basic Research Program Projects (Grant No. JCKY2019412D004).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no potential conflict of interest with respect to the research, authorship and publication of this article.

Abbreviations

FiCentrifugal force
FtSpring force
FmInertia force
FkCoe-type inertia force
FηViscous friction force
Fn1, Fn2Contact reaction force of the rotor
mvthe rotating blade center mass
ωthe spindle angular velocity
Rthe stator internal cavity large arc radius
rthe stator internal cavity small arc radius
ρ(θ)the blade vector diameter
kthe spring stiffness coefficient
μthe spring free length
mthe blade mass
nthe spindle speed
FfFriction force between the stator cavity and the stator
λcthe distance between the vane center of mass and the stator inner cavity
a(θ)the radial acceleration of the blade in the transition curve section
Lthe blade length
αthe pressure angle formed by the blade and the stator cavity surface
ƒthe friction coefficient
V(θ)the radial velocity of the blade in the transition curve section
ηthe oil viscosity
Abthe actual contact area between the bottom of the blade and the oil fluid according to the design size of the blade
Avsthe actual contact area between the blade and the endplate
Avr1the actual contact area between the blade and the rotor in the sizeable circular arc section
Avr2the actual contact area between the blade and the rotor in the small circular arc section
Avrthe actual contact area between the blade and the rotor in the stator transition curve section
hvrthe clearance between the blade and the rotor
hvsthe clearance between the blade and the endplate
L1the length of the blade extending outside the rotor slot
ξthe coefficient of the contact reaction force of the edge by the rotor on the edge
Fp1, Fp2, Fp3, Fp4, Fp5Hydraulic pressure

Appendix A

Since the hydraulic pressure on the suction cavity, pressure cavity wall surface and circular arc surface of the slide is changing with the rotation angle of the slide at all times during the rotation of the slide, the hydraulic pressure on the slide must be clarified as a function of the rotation angle during the force analysis of the slide and the stator cavity surface. Here, the finite element simulation method is used to derive the slider as a function of the change in the rotation angle.
The finite element solution sets the time step according to the different speed, and every 6° is 1 step—total 10 steps. The transition curve cavity slide fluid pressure results are analyzed as follows.
The force cloud diagram of the slide’s hydraulic pressure in the transition curve cavity section is shown in Figure A1.
Figure A1. Fluid pressure nebulae of the model with different slide angles.
Figure A1. Fluid pressure nebulae of the model with different slide angles.
Energies 15 06229 g0a1aEnergies 15 06229 g0a1b
The force cloud diagram of the hydraulic pressure section of the slide in the transition curve cavity section is shown in Figure A2.
Figure A2. Fluid pressure nephogram under different slide angles.
Figure A2. Fluid pressure nephogram under different slide angles.
Energies 15 06229 g0a2aEnergies 15 06229 g0a2b
(1)
Material mechanics analysis of slide section
According to the principle of material mechanics, the slider hydraulic pressure cloud can be simplified to a cantilever beam structure for linear load distribution and uniform load analysis of the circular part of the slider. As shown in Figure A3.
Figure A3. Cantilever beam load distribution of slide plate.
Figure A3. Cantilever beam load distribution of slide plate.
Energies 15 06229 g0a3
Hydraulic pressure analysis of the oil pressure chamber slider.
When the slider angle is θ .
F max θ = P max × S
where P max is the maximum pressure point on the wall of the slide in the pressure chamber, S is the unit area, F max θ denotes the maximum hydraulic pressure on the wall of the slide in the pressure chamber when the angle of rotation of the slide is θ. From the fluid simulation measurement results, it can be seen that the hydraulic pressure on the wall at the interface between the wall of the slide extension and the top profile arc is the maximum.
Length of rotor protruding from the slider: L 1
At this time, the equivalent liquid pressure at the wall of the pressure oil chamber slide
F p 1 = F max θ · L 1 2
where F p 1 is the equivalent fluid pressure at the wall of the pressure chamber slide.
Average hydraulic pressure in the circular part of the slider
q = P ¯ × S
where P ¯ is the average pressure on the circular part of the slide; q indicates the average hydraulic pressure on the circular part of the slide, measured by the fluid simulation results of the root of the arc to the stator contact point. The hydraulic pressure value gradually increases, and the assumption is that the average hydraulic pressure is the load on the circular part of the slide.
The arc length of the circle is
L 2 = R d · ( π 2 + α ) , R d = 3   mm
where L 2 is the length of the arc on the side of the oil pressure chamber at the contact point of the slide.
At this time, the pressure oil chamber slide arc equivalent fluid pressure
F p 2 = q · L 2 2
where F p 2 is the equivalent hydraulic pressure of the pressure chamber slider arc.
(2)
Suction chamber slider fluid pressure analysis.
When the angle of rotation of the slide is θ
F max θ = P max θ × S
where P max is the maximum pressure point on the wall of the suction chamber slide, S is the unit area, F max θ is the maximum hydraulic pressure on the wall of the suction chamber slide when the angle of rotation of the slide is θ. The fluid simulation measurement results show that the hydraulic pressure on the wall of the slide is highest at the point where the wall of the slide extends, and the top contour arc meets.
The equivalent liquid pressure at the wall of the suction chamber slide at this time
F p 3 = F max θ · L 1 2
where L 1 is the length of the slide extending out of the rotor at this rotation angle; Fp3 is the equivalent liquid pressure at the wall of the slide in the suction chamber.
The average hydraulic pressure in the circular part of the slider
q = P ¯ × S
where p ¯ is the average pressure applied to the circular part of the slide; q is the average fluid pressure in the circular part of the slide, measured by fluid simulation, which results in the root of the arc to the stator contact point. The fluid pressure value gradually increases, and the assumption is that the average fluid pressure is the load on the circular part of the slide.
Arc length of the circle
L 3 = R d · ( π 2 α ) , R d = 3   mm
where L 3 is the length of the arc on the suction chamber side of the slide contact point.
At this time, the suction chamber slider arc equivalent liquid pressure
F P 4 = q · L 3 2
where F p 4 is the suction chamber slider arc equivalent pressure.
Combined with the above transition curve cavity slide different angle of the case of the hydraulic pressure on the wall of the slide pressure cavity, the hydraulic pressure on the circular surface of the slide pressure cavity, the hydraulic pressure on the wall of the slide suction cavity and the hydraulic pressure on the circular surface of the slide suction cavity are shown in Table A1.
Table A1. Equivalent stress table of alternating fluid pressure on slide plate.
Table A1. Equivalent stress table of alternating fluid pressure on slide plate.
θ F max θ
(kN)
F p 1
(kN)
q
(kN)
L 2
(mm)
F p 2
(kN)
F max θ
(kN)
L 1
(mm)
F p 3
(kN)
q
(kN)
L 3
(mm)
F p 4
(kN)
207670.681984.710466.29286.4890.721474.71346.18
260842.43125.024783.74966.48311.041854.39406.63
12°250807.52765.432749.64996.46319.772103.98418.71
18°197602.82175.56.411562.561186.12361.081453.01218.22
24°162413.1178.57.049629.151325.10336.62452.37290.41
30°187328.181777.130631.271693.51296.591712.28195.52
36°197181.241986.421635.701971.84181.241982.99296.87
42°20454.061885.453512.592060.5354.59199.51.32131.86
48°18821.621087.693415.421980.2322.771711.72147.65
54°13433.51124.720264.341890.5047.251034.69242.02
The hydraulic pressure at the bottom of the slide has little effect on the overall force on the slide and can be abstracted as a constant value, taking the average hydraulic pressure at the bottom of the slide. According to the CFD results, the average pressure P5 at the bottom of the slide corresponding to the transition curve cavity is 0.122 MPa, and the area S3 at the bottom of the slide is 340.58 mm2, and the hydraulic pressure Fp5 at the bottom of the slide can be calculated according to Equation (A10) as 41,840,158.12 N.
F p 5 = P 5 × S 3
The function approximation method was used to obtain the function of liquid pressure in the transition curve cavity. According to the results data of the liquid pressure analysis, the function fitting was carried out using Matlab’s cftool toolbox to obtain the fitted function relationship equations of liquid pressure under different turning angles. The obtained function relations are shown in Equations (9)–(12) in the paper.

Appendix B

As shown in Figure A4, the contact angle between the top profile of the slide arc and the slide rotation angle. By finding the equation of the relationship between the two angles and sampling points to draw the curve, the function fitting method is applied to the curve to fit the function of the contact angle of the top profile of the slider arc and the slider rotation angle.
The equation of the polar coordinates of point A is show in Equation (A12).
{ x = ρ ( θ ) cos β y = ρ ( θ ) sin β
Figure A4. Stator contact angle diagram of slide piece.
Figure A4. Stator contact angle diagram of slide piece.
Energies 15 06229 g0a4
Set the slope of the tangent direction at the contact point A between the transition curve of the inner cavity of the sub and the top profile of the arc as k1.
k 1 = d y d x = d y d β · d β d x = ρ ( β ) sin β + ρ ( β ) cos β ρ ( β ) cos β ρ ( β ) sin β
Setting the direction of the common normal at contact point A, i.e., the slope k2 at n − n in the figure, it is calculated that
k 2 = 1 k 1 = ρ ( β ) sin β ρ ( β ) cos β ρ ( β ) sin β + ρ ( β ) cos β = f 1 ( β )
Then, the equation of the line O1A can be expressed as
y ρ ( β ) sin β = f 1 ( β ) ( x ρ ( β ) cos β )
Let the coordinates of the center O1 of the top profile arc of the slide be (x1, y1), and since the point (x1, y1) is on the line O1A, the following equation can be obtained.
y 1 x 1 = tan θ
The collation yields Equation (A17).
y 1 = x 1 · tan θ
Substituting y1 into the equation of the line O1A.
x 1 tan α ρ ( α ) sin β = f 1 ( β ) ( x 1 ρ ( β ) cos β )
The collation yields Equation (A19).
x 1 = f 1 ( β ) ρ ( β ) cos β + ρ ( β ) sin β tan θ f 1 ( β )
Let the radius of the top profile of the blade be Rd.
R d = ( y y 1 ) 2 + ( x x 1 ) 2 = f 3 ( θ , β )
Rd = 3 mm in this design model, and the derived Equation (A21) is obtained.
9 = ρ 2 ( β ) + ( sec 2 ( θ ) 2 ρ ( β ) ( sin β + cos β ) ) [ ρ ( β ) sin β ρ 2 ( β ) cos β sin β ρ ( β ) V ( β ) cos 2 β V ( β ) sin β + ρ ( β ) cos β tan α ρ ( β ) sin β V ( β ) cos β V ( β ) sin β + ρ ( β ) cos β ]
Since Rd is constant, θ is a function of β. β = f 4 ( θ ) .
The following geometric relationship can be obtained from Figure A4.
{ sin ( θ β ) = | A B | ρ ( β ) sin α = | A B | R d
The collation yields Equation (A23).
α = sin 1 ( ρ ( β ) sin ( θ β ) 3 )
Sampling 10 coordinate points using Matlab’s cftool toolbox for function fitting, comparing multiple types of fitting approximation methods, it was found that the Gaussian approximation method was the closest to the resultant data, fitting the function relationship between the slide contact angle α and the rotation angle θ Equation (A24).
Gaussian approximation method general function model.
f ( x ) = a 1 · exp ( ( ( x b 1 ) / c 1 ) 2 )
The coefficient (95% confidence level) is as follows.
a 1 = 48.11 ( 46.15 , 50.07 ) b 1 = 27.1 ( 26.64 , 27.57 ) c 1 = 13.89 ( 13.27 , 14.55 )
The result is as follows.
α = 48.11 · exp ( ( ( θ 27.1 ) / 13.89 ) 2 )

References

  1. Ma, W.C. Research status of high-pressure vane pump and its outlook. Fluid Drive Control 2006, 6, 1–6. [Google Scholar]
  2. Liu, J.G. Development of double-acting variable vane pump. J. Yanshan Univ. 2002, 3, 276–279. [Google Scholar]
  3. Jing, Y. Sliding vane pump features and development trend. China Storage Transp. 2008, 12, 116–117. [Google Scholar]
  4. Huang, S.B. Innovative Design and Simulation for a Novel Vane Pump. In Proceedings of the International Conference on Mechatronics and Materials Processing, Guangzhou, China, 18–20 November 2011; pp. 354–359. [Google Scholar]
  5. Karpenko, M.; Prentkovskis, O.; Šukevičius, Š. Research on high-pressure hose with repairing fitting and influence on energy parameter of the hydraulic drive. Eksploat. I Niezawodn. Maint. Reliab. 2022, 24, 25–32. [Google Scholar] [CrossRef]
  6. Stosiak, M. Ways of reducing the impact of mechanical vibrations on hydraulic valves. Arch. Civ. Mech. Eng. 2015, 15, 392–400. [Google Scholar] [CrossRef]
  7. Bai, J.J.; Xu, M.M. Research Status and Development Trend of Vane Pump. In Proceedings of the International Conference on Mechatronics Engineering and Computing Technology, Shanghai, China, 9–10 April 2014; pp. 1143–1146. [Google Scholar]
  8. Yoshida, Y.; Tsujimoto, Y.; Kawakami, T.; Sakatani, T. Unbalance hydraulic forces caused by geometrical manufacturing deviations of centrifugal impellers. J. Fluids Eng. 1998, 120, 531. [Google Scholar] [CrossRef]
  9. Ravindra, B.; Appasaheb, K. Prediction of flow-induced vibration due to impeller hydraulic unbalance in vertical turbine pumps using one-way fluid-structure interaction. J. Vib. Eng. Technol. 2019, 8, 417–430. [Google Scholar]
  10. Li, G.P.; Chen, E.Y.; Yang, A.L.; Xie, Z.B.; Zhao, G.P. Effect of Guide Vanes on Flow and Vibroacoustic in an Axial-Flow Pump. Math. Probl. Eng. 2018, 2018, 3095890. [Google Scholar] [CrossRef]
  11. Christopher, S.; Kumaraswamy, S. Identification of Critical Net Positive Suction Head Form Noise and Vibration in a Radial Flow Pump for Different Leading Edge Profiles of the Vane. J. Fluids Eng. Trans. ASME 2013, 135, 121301. [Google Scholar] [CrossRef]
  12. Zhou, W.J.; Cao, Y.H.; Zhang, N.; Gao, B. A novel axial vibration model of multistage pump rotor system with dynamic force of balance disc. J. Vib. Eng. Technol. 2020, 8, 673–683. [Google Scholar] [CrossRef]
  13. Zhang, L.J.; Wang, S.; Yin, G.J.; Guan, C.N. Fluid-structure interaction analysis of fluid pressure pulsation and structural vibration features in a vertical axial pump. Adv. Mech. Eng. 2019, 11, 1687814019828585. [Google Scholar] [CrossRef]
  14. Liu, H.L.; Cheng, Z.M.; Ge, Z.P. Collaborative improvement of efficiency and noise of bionic bane centrifugal pump based on multi-objective optimization. Adv. Mech. Eng. 2021, 13, 1687814021994976. [Google Scholar] [CrossRef]
  15. Si, Q.R.; Wang, B.B.; Yuan, J.P. Numerical and Experimental Investigation on Radiated Noise Characteristics of the Multistage Centrifugal Pump. Processes 2019, 7, 793. [Google Scholar] [CrossRef]
  16. Cheng, Y.Q.; Wang, X.H.; Chai, H. The Theoretical Performance Analysis and Numerical Simulation of the Cylindrical Vane Pump. Arab. J. Sci. Eng. 2021, 46, 2947–2961. [Google Scholar] [CrossRef]
  17. Liu, L.; Ding, C.; Wang, P.F. Design and Optimization of the Transition Curve in the Novel Profiled Chamber Metering Pump. J. Multi-Body Dyn. 2020, 234, 435–446. [Google Scholar] [CrossRef]
  18. Fan, J.; Wang, Z.R. Calculation and analysis of single-acting vane pump stator radial fluid pressure for high-powered remote water supply fire truck water collection system. Fire Technol. Prod. Inf. 2013, 11, 47–50. [Google Scholar]
  19. Cho, I.S. Behavioral characteristics of the vane of a hydraulic vane pump for power steering systems. J. Mech. Sci. Technol. 2015, 29, 4483–4489. [Google Scholar] [CrossRef]
  20. Lei, H.; Hu, H.J.; Lu, Y. A dynamic analysis on the transition curve of profiled chamber metering pump. J. Dyn. Syst. Meas. Control 2016, 138, 071003. [Google Scholar] [CrossRef]
  21. Gao, P. Application of Matlab’s CFtool toolbox to flotation tailings ash and image grayness curve fitting. Coal Process. Technol. 2015, 1, 67–70. [Google Scholar]
  22. Hou, Z.R.; Lu, Z.S. MATLAB-based particle swarm optimization algorithm and its application. Comput. Simul. 2003, 10, 68–70. [Google Scholar]
  23. Wang, L.J.; Jiang, S.F.; Xu, F.X. Analysis of particle swarm algorithms and their covariates. J. Qiqihar Univ. 2019, 35, 14–17. [Google Scholar]
Figure 1. Vane pump cross-sectional configuration.
Figure 1. Vane pump cross-sectional configuration.
Energies 15 06229 g001
Figure 2. Stator wear location diagram.
Figure 2. Stator wear location diagram.
Energies 15 06229 g002
Figure 3. Blade transition cavity force analysis diagram.
Figure 3. Blade transition cavity force analysis diagram.
Energies 15 06229 g003
Figure 4. Slider transition chamber force analysis diagram.
Figure 4. Slider transition chamber force analysis diagram.
Energies 15 06229 g004
Figure 5. Optimized minimum value of sliding friction.
Figure 5. Optimized minimum value of sliding friction.
Energies 15 06229 g005
Figure 6. Optimized values of R, r and θ.
Figure 6. Optimized values of R, r and θ.
Energies 15 06229 g006
Figure 7. Stator cavity transition curve optimization comparison chart.
Figure 7. Stator cavity transition curve optimization comparison chart.
Energies 15 06229 g007
Figure 8. ADAMS model of sliding vane pump.
Figure 8. ADAMS model of sliding vane pump.
Energies 15 06229 g008
Figure 9. Comparison of spindle torque before stator optimization.
Figure 9. Comparison of spindle torque before stator optimization.
Energies 15 06229 g009
Figure 10. Comparison of spindle torque after stator optimization.
Figure 10. Comparison of spindle torque after stator optimization.
Energies 15 06229 g010
Figure 11. Comparison of slider acceleration before stator optimization.
Figure 11. Comparison of slider acceleration before stator optimization.
Energies 15 06229 g011
Figure 12. Comparison of slider acceleration after stator optimization.
Figure 12. Comparison of slider acceleration after stator optimization.
Energies 15 06229 g012
Figure 13. Comparison of friction between slider and rotor slider slot before stator optimization.
Figure 13. Comparison of friction between slider and rotor slider slot before stator optimization.
Energies 15 06229 g013
Figure 14. Comparison of friction between slider and rotor slider slot after stator optimization.
Figure 14. Comparison of friction between slider and rotor slider slot after stator optimization.
Energies 15 06229 g014
Figure 15. Comparison of slide and stator contact force before stator optimization.
Figure 15. Comparison of slide and stator contact force before stator optimization.
Energies 15 06229 g015
Figure 16. Comparison of slide and stator contact force after stator optimization.
Figure 16. Comparison of slide and stator contact force after stator optimization.
Energies 15 06229 g016
Table 1. Table of parameter values.
Table 1. Table of parameter values.
Ab
(mm2)
Avs
(mm2)
Avr1
(mm2)
Avr2
(mm2)
η
(mm2/s)
hvr
(mm)
hvs
(mm)
340.5880.03423.54910300.030.30
Table 2. Constraints.
Table 2. Constraints.
aω
(rad/s)
ƒαr
(mm)
R
(mm)
R/r
(Rr)
50°1000.13[0, 90°][32, 34][32, 40](1, 1.2603]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sun, Y.; Xue, D.; Liu, S.; Wu, J.; Bai, X. Stator Curvature Optimization and Analysis of Axial Hydraulic Vane Pumps. Energies 2022, 15, 6229. https://doi.org/10.3390/en15176229

AMA Style

Sun Y, Xue D, Liu S, Wu J, Bai X. Stator Curvature Optimization and Analysis of Axial Hydraulic Vane Pumps. Energies. 2022; 15(17):6229. https://doi.org/10.3390/en15176229

Chicago/Turabian Style

Sun, Yongguo, Dong Xue, Shisheng Liu, Jinghang Wu, and Xingyu Bai. 2022. "Stator Curvature Optimization and Analysis of Axial Hydraulic Vane Pumps" Energies 15, no. 17: 6229. https://doi.org/10.3390/en15176229

APA Style

Sun, Y., Xue, D., Liu, S., Wu, J., & Bai, X. (2022). Stator Curvature Optimization and Analysis of Axial Hydraulic Vane Pumps. Energies, 15(17), 6229. https://doi.org/10.3390/en15176229

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop