Next Article in Journal
A Performance Evaluation Approach for Satellite Attitude Control System in Tracking Mode
Next Article in Special Issue
Investigation of Innovative High-Response Piezoelectric Actuator Used as Smart Actuator–Sensor System
Previous Article in Journal
Unveiling the Biomechanical Insights: Motor Control Shifts Induced by Shoe Friction Adjustments and Their Impact on Defensive Slide, Crossover Dribbling, and Full Approach Jump in Basketball
Previous Article in Special Issue
Ionic Hydraulic Fluids and Seal-Material Compatibility
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Underlapped Symmetrically Ported Valve-Controlled Asymmetric Cylinder Drive

College of Agricultural Equipment Engineering, Henan University of Science and Technology, Luoyang 471000, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(7), 2870; https://doi.org/10.3390/app14072870
Submission received: 21 December 2023 / Revised: 3 February 2024 / Accepted: 27 March 2024 / Published: 28 March 2024
(This article belongs to the Special Issue Research Progress on Hydraulic Fluid and Hydraulic Systems)

Abstract

:
The valve-controlled cylinder drive system is the most common type among hydraulic applications. Nonlinear behaviour in such systems is inevitable when the valve spool is around its null position. We utilised the component linking method to investigate the nonlinearities in a Moog valve-controlled asymmetric cylinder drive system by simulation in Fortran, in which a generalised concept is introduced and validated by comparing to the experimental results. An X factor is proposed in the generalised concept to describe the asymmetric cylinder state, which is a constant when the cylinder is extending or retracting, but numerically calculated when the valve spool is in the underlap region. This analytical solution is approximately 200 times more computationally efficient than the numerical solution method. This paper utilises the component linking method to simulate the Moog valve-controlled asymmetric cylinder drive system in Matlab Simulink, and proposes an analytical solution for the X factor when the valve spool is in the underlap region. This analytical solution is approximately 200 times more computationally efficient than the numerical solution method.

1. Introduction

Hydraulic power systems [1] are widely used in industrial manufacture, building construction, port transportation, funfairs, etc. They use fluid to achieve energy transmission without a rigid joint. Hydraulic power is able to produce massive force and torque within a compact and lightweight package, and conventional rigid components like gears and joints can be avoided.
Hydraulic linear actuators are mainly of the symmetric cylinder or asymmetric cylinder type. The symmetric cylinder type has equal piston areas and the asymmetric type (also called single rod) cylinder comes with different piston areas. The symmetric cylinder is usually utilised for high-performance applications, due to its linear performance with a symmetric ported control valve. This type of cylinder is widely applied in aerospace and military contexts for its high performance and dynamic response [2]. However, the asymmetric cylinder [3] drive constitutes 80% of contemporary hydraulic applications [4], such as cranes, injection moulding systems, hydraulic press machines, etc. The asymmetric cylinder drive is usually controlled by a symmetric ported valve, meaning an unequal flow rate will be delivered into the cylinder chambers with the same valve opening in the extending or retracting state, which means system parameters do not remain the same after the valve crosses its null position [5]. Leaney [6] studied these nonlinearities, but how the system performs when the valve is in the underlap region is not explicitly explained. Viersma [7] indicates that pressure jumps will occur when the valve spool is around its null position, causing ‘implosion’ or ‘explosion’ of oil owing to its compressibility, so that smooth operation around the valve’s null position is impossible. Therefore, this article conducted a nonlinear study on valves in an underlap state and proposed an analytical solution for X factor when the valve spool is in the underlap region to improve computational efficiency.

2. Theory

A test rig of a four-way servovalve-controlled asymmetric cylinder hydraulic system was constructed by Leaney [8]. The system is a typical symmetric ported valve-controlled asymmetric cylinder system, which is depicted in Figure 1. The Moog series 76 model 102 four-way valve is used to control the asymmetric cylinder system to achieve the desired output.
Leaney [6] simulated his test rig with the ‘Component Linking’ method, which models each component such as the valve, actuator/cylinder, and load, then combines them to perform the whole-system simulation. This method aims at simplifying the ‘Small Perturbation’ model [9], which utilised a single transfer function to describe a whole hydraulic system. However, as not all parameters of a hydraulic system remain constant during operation and a single transfer function is mainly designed for a linear time-invariant (LTI) system, a single transfer function [10] is inadequate to capture the behaviours of a hydraulic system. Modelling the components separately is easier for the analysis and adjustment of the system parameters.
The components of the asymmetric cylinder system are connected by the concept of ‘power bonding’, in other words, the components are connected by the power transfer; such a concept is still utilised nowadays for the hydraulic system of a 50-wheel loader [11]. In Leaney’s paper [8], three components are the major concerns, which are the valve, the actuator, and the load model. Their power transfers are indicated in Figure 2. The voltage is input to the valve model to control the valve opening, after which the flow passes through the valve and drops to a certain value. Pressure and flow are delivered to the actuator model, which outputs the force/torque and velocity. Finally, the load model outputs the velocity and force.

2.1. Valve Model

The control valve (electrohydraulic valve) is a Moog 76 model 102 flow control four-way servovalve. If current is applied to the armature coil in the torque motor, the T-shaped flapper rotates slightly. Assuming the flapper rotates clockwise, the offset between the left nozzle and flapper decreases, which leads to the decrease of flow passing through the left nozzle and causes the pressure to increase in the left control line. A similar process occurs in the right control line, leading to a decrease of pressure.
The pressure difference in the control line moves the spool, therefore bends the feedback spring to balance the torque generated by the torque motor, so that the spool displacement is dependent on the rotation of the torque motor; furthermore, the spool displacement is dependent on the current.
Appropriate transfer functions, Equations (1) and (2), are given by Moog [12]. The expressions are linear, empirical relationships that approximate the responses of the actual servovalve. For the frequency range below 50 Hz, the valve dynamic can be approximated by a first-order transfer function given below [12]:
Q ( s ) I ( s ) = K 1 1 + τ s s
where K is the flow gain m 3 s   ·   amp and τ s is the servovalve time constant sec .
A second-order transfer function is necessary to represent a wider frequency range response as given below [12]:
Q s I s = K 1 1 + 2 ζ ω n s + s ω n 2
where ω n is the natural frequency [ rad / s ], and ζ is the damping ratio. In a limit cycle oscillation study of a low-frequency phenomenon, a higher-order valve model does not show significant difference compared to a lower-order model [13]. Hisa [14] stated that a reduced-order transfer function can be used to approximate a higher-order transfer function, and a first-order model is used as a reference model by Persson et al. [15].
As all of the valve operational frequencies in the industrial sector are much lower than 50 Hz, the first-order transfer function Equation (1) is chosen for the valve model. As the displacement of the spool, X s , controls the flow rate of the valve, Equation (1) can be modified as follows:
X s s I s = K s 1 1 + τ s s
where K s is the gain of the system [ m amp ].
The current is generated by an amplifier. Assuming the amplifier is ideal, then current is proportional to the input voltage as given by Equation (4):
I s V v s = K A
where K A is the gain of voltage to current [ amp volt ]. Hence the valve dynamics are represented by:
X s s V v s = K A K s 1 + τ s s
A block diagram for this valve model is depicted in Figure 3. Valve model block diagram.

2.2. Flow Gain

The flow gain is the gain that describes the relationship between the valve opening and the oil flow across it. In order to interpret it clearly, a physical layout of the four-way valve-controlled [16] asymmetric cylinder is presented in Figure 4. Physical layout of the four-way valve-controlled asymmetric cylinder [6].
Q 1 is the flow rate that is delivered into the piston side chamber, V 1 is the volume of fluid that is trapped in the piston chamber, A 1 is the piston side area of the cylinder, and the subscript 2 indicates corresponding parameters in the rod side chamber. P s is the supply pressure, and P T is the return line pressure of the oil tank, which is zero. A V is the opening area of the corresponding port of the valve, and subscripts 1, 2, 3, and 4 indicate different ports of the valve which are depicted in Figure 4. Physical layout of the four-way valve-controlled asymmetric cylinder [6].
x P is the cylinder displacement. Q P 1 is the flow passing through the A V 1 port, Q 1 T is the flow passing through the A V 2 port, Q 2 T is the flow passing through the A V 3 port, and Q P 2 is the flow passing though the A V 4 port.
The flow delivered into the chambers can be described as:
Q 1 = Q P 1 Q 1 T = C D 2 ρ A V 1 P s P 1 1 2 A V 2 P 1 1 2
Q 2 = Q 2 T Q P 2 = C D 2 ρ A V 3 P 2 1 2 A V 4 P s P 2 1 2
For a symmetric ported valve, A V 1 = A V 3 and A V 2 = A V 4 , so that Equations (6) and (7) can be simplified as:
Q 1 = Q P 1 Q 1 T = C D 2 ρ A V 1 P s P 1 1 2 A V 2 P 1 1 2
Q 2 = Q 2 T Q P 2 = C D 2 ρ A V 1 P 2 1 2 A V 2 P s P 2 1 2
These two equations can be applied when the valve is in the extending, retracting, and underlap region, the underlap region being the valve spool position where all port opening areas are non-zero. Even for a zero-lapped valve, due to the manufacturing tolerances [17], there will be a small underlap when the spool is placed in the null position. Figure 5 depicts the theoretical and actual situation for a nominally zero-lapped valve.
The manufacturing tolerances can be identified as:
i.
Radial clearance between spool and sleeve.
ii.
Tolerance on lap.
Such tolerance [1] will lead to variation of flow gain when the valve is in its null region. The flow gain model block diagram is depicted in Figure 6.

2.3. The Actuator Model

The actuator model connects the valve and load model, which accepts the information of flow and pressure as input, then outputs force to the load model. The types of hydraulic actuators are mainly symmetric cylinder, asymmetric cylinder, and hydraulic motor. Due to the focus of this research, the asymmetric cylinder is the only one under analysis.
Whether the asymmetric cylinder is in the extending or retracting state, the flow in its chamber consists of input flow, displacement flow, and compressibility flow. However, under different operational states, the flow balance equations in the different chambers are not the same.
For the cylinder extending state:
Q 1 = Q D 1 + Q C 1
Q 2 = Q D 2 Q C 2
For the cylinder retracting state:
Q 2 = Q D 2 + Q C 2
Q 1 = Q D 1 Q C 1
where
Q D 1 = v   ·   A 1       Q D 2 = v   ·   A 2
Q C 1 = s P 1   ·   V 1 B       Q C 2 = s P 2   ·   V 2 B
A is the area of the cylinder piston side area or rod side area [ m 2 ], v is the cylinder velocity [ m s ], V is the cylinder chamber volume [ m 3 ], P is the cylinder chamber pressure [ N m 2 ] and B is the fluid Bulking Modulus [ N m 2 ]. The subscript 1 and 2 indicate the piston side chamber and rod side chamber, respectively.
A generalised definition of these equations is proposed by Leaney [8], which uses effective parameters to describe the system behaviours rather than analyse them separately in different cylinder chambers. The generalised equations are depicted as follows:
P L E = P 1 A 1 P 2 A 2 A E       Q L E = Q 1 + Q 2 2 A E = A 1 + A 2 2       Q D E = Q D 1 + Q D 2 2 Q C E = Q C 1 + Q C 2 2 = Q L E Q D E
Substituting these effective parameters into the flow and pressure equations reveals the following:
For cylinder extending ( A V 2 = 0 and A V 1 > 0), the service line pressures can be rewritten as follows [8]:
P 1 e x t = P s + γ 2 1 2 1 + γ P L E 1 + γ 3
P 2 e x t = γ P s 1 2 1 + γ P L E 1 + γ 3
For cylinder retracting ( A V 2 > 0 and A V 1 = 0), the service line pressures can be rewritten as follows [8]:
P 1 r e t = γ 2 P s + γ 2 1 2 1 + γ P L E 1 + γ 3
P 2 r e t = γ 3 P s 1 2 1 + γ P L E 1 + γ 3
where γ is the area ratio of the cylinder. Leaney [8] introduced the X factor to unify the expression for P 1 , P 2 and utilised generalised concepts in the component linking method. The X factor value depends on the region of the spool operation: X = 1 for cylinder extending state and X = γ 2 for cylinder retracting state. The service line pressure P 1 and P 2 can be written as follows:
P 1 = X P s + γ 2 1 2 1 + γ P L E 1 + γ 3
P 2 = X γ P s 1 2 1 + γ P L E 1 + γ 3  
P L E = 4 B s V T E Q C E
where γ is the area ratio, V T E is effective total volume [ m 3 ], s is the Laplace operator, and PLE is average pressure. When X = 1 for cylinder extending state, X = γ2 for cylinder retracting state.
V T E = 2 A E A 1 V 1 + A 2 A 1 2 V 2 1 + A 2 A 1 3
The X factor value is known and clear when the valve is outside the underlap region. However, when the valve spool is in the underlap region, the X value is difficult to find explicitly; it can be numerically solved by an implicit Equation (25) [8] as shown below:
A v 1 1 + γ 3 X 1 2 A v 2 X 1 2 = γ A v 1 X γ 1 2 A v 2 1 + γ 3 X γ 1 2
Now the effective load flow can be revealed [8]:
Q L E = C D 2 ρ 1 + γ 2 1 + γ 3 A v 1 X γ P s 1 2 1 + γ P L E 1 2 A v 2 1 + γ 3 X γ P s + 1 2 1 + γ P L E 1 2
The actuator model block diagram is depicted in Figure 7.

2.4. The Load Model

The load model varies depending on the test rig setup. In order to simplify the analysis afterwards, the load setup is a simple asymmetric cylinder coupled with a load, which is placed on a surface, so that the load force balance can be depicted by Equation (27) and Figure 8:
P 1 A 1 P 2 A 2 = P L E A E = m   ·   a + f r i c t i o n   f o r c e + F e x t e r n a l
where m is the load mass [ kg ] and a is the acceleration [ m / s 2 ].

3. Materials and Methods

3.1. Numerical Solution of Chamber Pressure Behaviour in Underlap Region

Leaney [8] stated that the cylinder chamber pressures P 1 and P 2 must be numerically solved when the valve is in the underlap region, due to the derived implicit Equation (25). The chamber pressures are the power source to drive the cylinder, so that the system behaviours are decided by chamber pressures. The way the pressures perform when the valve is outside the underlap region is explicit. The chamber pressures’ behaviour in the underlap region are solved numerically as shown in Figure 9.
An input signal is set to drive the valve to move across the underlap region in a small time interval with a constant speed. The supply pressure is 70 bar and all the modelling parameters are the same as those in the Moog valve-controlled asymmetric cylinder system in the last section.
The chamber pressures of the asymmetric cylinder when the valve is in the underlap region are modelled, and they are numerically solved. The key part is the X factor Equation (25) from [8]. The solving command in Matlab is fsolve and it uses the rust-region algorithm [14] to find the root of factor X with given valve opening areas A v 1 , A v 2 and area ratio γ . Numerically solving Equation (25) takes a lot of computing time; if the X value can be analytically solved, simulation time cost can be massively reduced. A trial is given in the rest of this section.

3.2. Derivation of Analytical Solution for the Value of X Factor

Squaring both sides of Equation (25) leads to:
A V 1 2 γ 2 1 + γ 3 X + A V 2 2 X γ 2 2 A V 1 A V 2 γ 2 X + γ 3 X X 2 1 2 = A V 1 2 X γ + A V 2 2 1 + γ 3 X γ 2 A V 1 A V 2 X γ + X γ 4 X 2 γ 2 1 2
Rearranged:
A V 1 2 A V 2 1 + γ 3 A V 2 2 A V 1 γ 2 + γ 5 A V 1 2 A V 2 1 + γ 3 A V 2 2 A V 1 1 + γ 3 X = 1 + γ 3 X X 2 1 2 γ 2 γ + γ 4 X γ 2 X 2 1 2 .
Let m be given as follows:
m = A V 1 2 A V 2 1 + γ 3 A V 2 2 A V 1 γ 2 + γ 5 = 1 + γ 3 A V 1 2 A V 2 A V 2 γ 2 2 A V 1
Let n be given as follows:
n = A V 1 2 A V 2 1 + γ 3 A V 2 2 A V 1 1 + γ 3 = 1 + γ 3 A V 1 2 A V 2 A V 2 2 A V 1
Substituting m and n into Equation (31) reveals the following:
m n X = 1 + γ 3 X X 2 1 2 γ 2 γ + γ 4 X γ 2 X 2 1 2
Squaring both sides of Equation (32) leads to:
m 2 + n 2 X 2 2 m n X = 1 + γ 3 X X 2 + γ 5 1 + γ 3 X γ 6 X 2   2 γ 2 γ 1 + γ 3 2 X 2 γ 1 + γ 3 X 3 γ 2 1 + γ 3 X 3 + γ 2 X 4 1 2
Rearranging Equation (33) reveals the following:
n 2 + 1 + γ 6 X 2 2 m n + 1 + γ 5 1 + γ 3 X + m 2   = 2 γ 2 γ 1 + γ 3 2 X 2 γ 1 + γ 3 X 3 γ 2 1 + γ 3 X 3 + γ 2 X 4 1 2
Let
p = n 2 + 1 + γ 6 q = 2 m n + 1 + γ 5 1 + γ 3 h = m 2
Equation (34) becomes the following:
p X 2 q X + h = 2 γ 2 γ 1 + γ 3 2 X 2 γ 1 + γ 3 X 3 γ 2 1 + γ 3 X 3 + γ 2 X 4 1 2
Squaring both sides of Equation (36) reveals the following:
p 2 X 4 + q 2 X 2 + h 2 p q X 3 + p h X 2 p q X 3 q h X + p h X 2 q h X   = 4 γ 6 X 4 4 γ 5 1 + γ 1 + γ 3 X 3 + 4 γ 5 1 + γ 3 2 X 2
Rearranging the above equation reveals the following:
p 2 4 γ 6 X 4 + 4 γ 5 1 + γ 1 + γ 3 2 p q X 3 + 2 p h + q 2 4 γ 5 1 + γ 3 2 X 2   2 q h X + h 2 = 0
Let
a 4 = p 2 4 γ 6 a 3 = 4 γ 5 1 + γ 1 + γ 3 2 p q a 2 = 2 p h + q 2 4 γ 5 1 + γ 3 2 a 1 = 2 q h a 0 = h 2
Thus Equation (38) becomes the following:
a 4 X 4 + a 3 X 3 + a 2 X 2 + a 1 X + a 0 = 0
Let
a = a 3 a 4 , b = a 2 a 4 c = a 1 a 4 , d = a 0 a 4
As a result, Equation (40) becomes the following:
X 4 + a X 3 + b X 2 + c X + d = 0
The value of the X factor is a problem of finding the root of a quartic Equation (42).

3.3. Root Form and Identification

A general form of roots can be depicted as below in Equation (43):
X 1 = a 4 + U + V 1 2 2 2 U V W 1 2 2 X 2 = a 4 U + V 1 2 2 2 U V + W 1 2 2 X 3 = a 4 + U + V 1 2 2 + 2 U V W 1 2 2 X 4 = a 4 U + V 1 2 2 + 2 U V + W 1 2 2
where
U = a 2 4 2 b 3 V = L 3 + Z 0 3 L W = a 3 4 a b + 8 c 4 U + V 1 2
where
Z 0 = b 2 3 a c + 12 d Z 1 = 2 b 3 9 a b c + 27 a 2 d + 27 c 2 72 b d L = Z 1 + Z 1 2 4 Z 0 3 1 2 2 1 3
Actually, there is only one value of factor X with a given set of parameters. However, there are four roots from a quartic equation and the value of X is one of them.
To identify the correct root, the method is depicted as the following steps:
(1)
Numerically solving the X factor in the underlap region.
(2)
Identify one of the roots in a certain valve displacement in the underlap region corresponding to the same value from the numerical solution.
(3)
Repeat step 1 and 2 with different area ratio γ
(4)
Repeat step 1, 2, and 3 with different underlap regions.
(5)
Estimate the final analytical solution of the X factor.
Starting with the parameter from Moog’s valve-controlled asymmetric cylinder drive system, for the area ratio γ = 11.43/5.8 = 1.9707 and underlap region from 0.0045 mm to −0.0045 mm, root identification is listed in Table 1.
A new phenomenon can be noticed: the X factor values in the underlap region contain two roots of Equation (25), X 1 and X 4 . When valve displacement is from −0.0045 mm to −0.00075 mm, the value of the X factor is the root X 4 , and when valve displacement is from −0.00074 mm to 0.0045 mm, the value of the X factor is the root X 1 .
All four root values in the underlap region are described as in Figure 10a, different line styles corresponding to different roots. The values of the roots in the underlap are compared with the numerically solved X values, as in Figure 10b. It can be noticed that a part of X 4 and a part of X 1 are merged into the numerically solved curve, and they meet at some specific point.
Another interesting phenomenon can be observed: that some parts of the root curve can be combined with parts of other root curves to generate a smoothed curve. The point where X 4 and X 1 curves meet may change with some parameters. We can change the form of valve opening area A v 1 and A v 2 in Equation (25) into the equation below.
π D x s + 0.0000045 1 + γ 3 X 1 2 π D x s 0.0000045 X 1 2 = γ π D x s + 0.0000045 X γ 1 2 π D x s 0.0000045 1 + γ 3 X γ 1 2
where A v 1 = π D x s + 0.0000045 , A v 2 = π D x s 0.0000045 [ m 2 ].
Cancelling the π D from the above equation, and dividing both sides by 0.0000045 reveals
1 + x s 0.0000045 1 + γ 3 X 1 2 1 x s 0.0000045 X 1 2 = γ 1 + x s 0.0000045 X γ 1 2 1 x s 0.0000045 1 + γ 3 X γ 1 2
The x s 0.0000045 in Equation (47) can be regarded as normalized spool travel from −1 to +1 in underlap. It can be observed that for a certain area ratio γ, the X factor is only affected by the travel percentage. An assumption can be made that the intersection point in Figure 10a is affected by the area ratio; different area ratios γ are applied into Equation (47) to verify this. After running the simulation with different area ratios γ in the underlap region, the root X 4 and X 1 distribution is depicted in Table 2.
Inputting these data into the Matlab CFTOOL (https://www.mathworks.com/products/curvefitting.html) reveals the relationship for X 1 percentage with different area ratio γ as shown in Equation (48).
X 1   P e r c e n t a g e = 0.00999 γ 3 0.09352 γ 2 + 0.3228 γ + 0.2328
The percentage of root X 4 is simply obtained by 1 X 1   Percentage , so that the X factor root distribution when the valve is in the underlap region is obtained.

3.4. Validation of Analytical Solving Method

The overall component linking model block diagram for this Moog valve-controlled asymmetric cylinder drive is depicted in Figure 11.
The overall block model of the Moog valve-controlled asymmetric cylinder drive system is simulated in Matlab Simulink as shown in Figure 12.
The Simulink model with analytical solution simulation results is compared with the numerical solution simulation results in Fortran and the experimental test results in Figure 13 and Figure 14.
The modelling results from Simulink show consistency with the modelling results from Fortran and experimental results [8], which indicates that the analytical solving method is validated.

4. Results and Discussion

The analytical solution for the X factor in the component linking method is able to save computing time during simulation. The xs valve stroke keeps changing during operation, which indicates the computer must solve Equation (25) with different values of xs. If the step size of the simulation is small, for instance 1∗10−4 s, the computer will need to solve thousands of implicit equations in one second; this process will consume a lot of computing power. The analytical solution in the second section improves this situation.
Simulation tests of the Moog valve-controlled asymmetric cylinder drive system are carried out to compare the time consumed between the implicit solving and analytical solving methods, including a square wave test, a sawtooth test, and a sine wave test.
There is no observable difference between the two results, which indicates both solutions have the same performance. The numerical solution takes 5.49 s, while the analytical solution takes 1.07 s, therefore the numerical solution consumes nearly five times as much time to finish the simulation. However, the valve in the square wave test in Figure 15 is operating outside the underlap region most of the time, therefore the advantage of the analytical solution is not highlighted here.
The sawtooth command test is employed to compare the time cost between the numerical solution and the analytical solution when the valve spool sweeps across its underlap region. A total of 2.31% of the valve spool’s travel in the underlap region is considered in the sawtooth test, and the simulation results of the sawtooth command are depicted in Figure 16.
The analytical solution in the sawtooth test produces the same results as the numerical solution, but the analytical solution costs 1.23 s to finish and the numerical method costs 8.26 s to finish. Nonlinear behaviours can be observed in the results in Figure 16, which include the gradient change of the velocity curve after the system switches its state and some observed oscillations when the valve is in the underlap region around its null position.
Now, assuming the underlap region of the Moog four-way is 20%, a sine wave motion command is sent to the four-way valve so that it only operates in the underlap region. Its pressure behaviour under numerical and analytical simulation results are depicted in Figure 17.
The modelling time costs of calculating chamber pressures when valve travel is in the underlap region are compared in Figure 17. The numerically solved simulation used 1647.17 s to finish, while the analytical solution only uses 8.50 s, which is a significant reduction by a factor of 200.
The time consumptions of the three tests are depicted in Table 3.

5. Conclusions

The component linking method is validated in Matlab Simulink (https://www.mathworks.com/products/simulink.html), and its simulation results show consistency with the published results. Therefore, it can be utilised to achieve further research when the valve of the asymmetric cylinder drive system is in the underlap region. How the chamber pressures behave when the valve is in the underlap region is revealed by Simulink modelling, and the inevitable pressure jump process when the valve switches its state is numerically solved with the X factor method. To reduce the computing power required, an analytical solution for the X factor is developed. Equation (47) proved the X factor value is only related to the percentage of valve spool travel in the underlap region, and offered a solution to identify the root distribution percentage in the underlap region. This analytical solution is approximately 200 times more computationally efficient than the numerical solution method.
Underlapped symmetrically ported valves are widely used in systems that require precise flow control, such as chemical processing, liquid distribution systems, etc. By adjusting the opening of the valve, the flow rate of the medium can be precisely controlled. By testing the pressure process formula, it is possible to infer how much displacement the valve core has produced, thereby correcting the valve’s lifespan status.

Author Contributions

Conceptualization, methodology, resources, investigation, supervision, and formal analysis, H.W.; software, M.X.; validation, H.W., M.X. and Z.C.; data curation, M.X. and Z.C; writing—original draft preparation and writing—review and editing, H.W. and M.X; visualization, Z.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Manring, N.D.; Fales, R.C. Hydraulic Control Systems; John Wiley & Sons: Hoboken, NJ, USA, 2019. [Google Scholar]
  2. Burston, M.; Ranasinghe, K.; Gardi, A.; Parezanović, V.; Ajaj, R.; Sabatini, R. Design principles and digital control of advanced distributed propulsion systems. Energy 2022, 241, 122788. [Google Scholar] [CrossRef]
  3. Pan, L.P.; Mi, Z.N.; Ju, S.J. Dynamic Characteristics Analysis on Symmetric Control of Asymmetric Heavy-duty Hydraulic Cylinder. Appl. Mech. Mater. 2013, 336, 924–927. [Google Scholar] [CrossRef]
  4. Vacca, A.; Franzoni, G. Hydraulic Fluid Power: Fundamentals, Applications, and Circuit Design; John Wiley & Sons: Hoboken, NJ, USA, 2021. [Google Scholar]
  5. Ba, K.X.; Yu, B.; Kong, X.D.; Zhao, H.-L.; Zhao, J.-S.; Zhu, Q.-X.; Li, C.-H. The dynamic compliance and its compensation control research of the highly integrated valve-controlled cylinder position control system. Int. J. Control Autom. Syst. 2017, 15, 1814–1825. [Google Scholar] [CrossRef]
  6. Leaney, P.G. The Modelling and Computer Aided Design of Hydraulic Servosystems. Ph.D. Thesis, Loughborough University, Loughborough, UK, 1986. [Google Scholar]
  7. Aboelela, M.A.S.; Essa, M.E.S.M.; Hassan, M.A.M. Modeling and identification of hydraulic servo systems. Int. J. Model. Simul. 2018, 38, 139–149. [Google Scholar] [CrossRef]
  8. Leaney, P.G. Component Oriented Modelling of a Valve Controlled Asymmetric Cylinder Drive Using a Generalized Formulation for Model Linking. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 1991, 205, 285–298. [Google Scholar] [CrossRef]
  9. Gao, F.; Luo, Y.; Pan, Z.; Yang, R. Research on Clamping Force Control of Hydraulic Power Chuck of Numerical Control Machine Tool. In Proceedings of the 2020 IEEE 6th International Conference on Control Science and Systems Engineering (ICCSSE), Beijing, China, 17–19 July 2020; IEEE: New York, NY, USA, 2020; pp. 118–123. [Google Scholar]
  10. Colombo, M.; Cammi, A.; Memoli, V.; Papini, D.; Ricotti, M.E. Transfer function modelling of the Lead-cooled Fast Reactor (LFR) dynamics. Prog. Nucl. Energy 2010, 52, 715–729. [Google Scholar] [CrossRef]
  11. Jing, B.; Luo, Z.; Zhu, X.; Dong, L. Shift hydraulic system modeling and simulation base on power bond graph. In Proceedings of the International Conference on Computer, Mechatronics, Control and Electronic Engineering, Changchun, China, 24–26 August 2010; Volume 5, pp. 111–114. [Google Scholar]
  12. Thayer, W.J. Transfer functions for Moog servovalves. Moog Tech. Bull. 1965, 103, 1–11. [Google Scholar]
  13. Albaidaq, H.; Bjørnarå, O.P.T. Design, Modelling and Simulation of Electro-Hydraulic Self-Contained Cylinders based on Digital Hydraulics. Master’s Thesis, University of Agder, Kristiansand, Norway, 2018. [Google Scholar]
  14. Gao, H.; Lam, J.; Wang, C. Model simplification for switched hybrid systems. Syst. Control Lett. 2006, 55, 1015–1021. [Google Scholar] [CrossRef]
  15. Persson, L.; Plummer, A.; Bowen, C.; Elliott, P. Non-linear Control of a Piezoelectric Two Stage Servovalve. In Proceedings of the 15th Scandinavian International Conference on Fluid Power, Linköping, Sweden, 7–9 June 2017. [Google Scholar]
  16. Shi, A.-Z. Dynamic analysis of symmetrical four-way valve controlled non-symmetrical hydraulic cylinder. Hoisting and Conveying Machinery, 2006.
  17. Zheng, C.; Ma, C.; Yan, X.; Cui, Z.-C.; Zhang, D.-M. Manufacture tolerance analysis and control for a polymer-on-silicon Mach-Zehnder-interferometer-based electro-optic switch. Optoelectron. Lett. 2011, 7, 101–104. [Google Scholar] [CrossRef]
Figure 1. Schematic circuit of the Moog valve-controlled asymmetric cylinder system [8].
Figure 1. Schematic circuit of the Moog valve-controlled asymmetric cylinder system [8].
Applsci 14 02870 g001
Figure 2. Power transfer route.
Figure 2. Power transfer route.
Applsci 14 02870 g002
Figure 3. Valve model block diagram.
Figure 3. Valve model block diagram.
Applsci 14 02870 g003
Figure 4. Physical layout of the four-way valve-controlled asymmetric cylinder [6].
Figure 4. Physical layout of the four-way valve-controlled asymmetric cylinder [6].
Applsci 14 02870 g004
Figure 5. Theoretical and actual situation of a zero-lapped valve [6].
Figure 5. Theoretical and actual situation of a zero-lapped valve [6].
Applsci 14 02870 g005
Figure 6. Flow gain model block diagram.
Figure 6. Flow gain model block diagram.
Applsci 14 02870 g006
Figure 7. The actuator model block diagram.
Figure 7. The actuator model block diagram.
Applsci 14 02870 g007
Figure 8. The load model block diagram.
Figure 8. The load model block diagram.
Applsci 14 02870 g008
Figure 9. Chamber pressures when the valve is in underlap state.
Figure 9. Chamber pressures when the valve is in underlap state.
Applsci 14 02870 g009
Figure 10. Four X roots in underlap region (a) and numerically solved X factor (b).
Figure 10. Four X roots in underlap region (a) and numerically solved X factor (b).
Applsci 14 02870 g010
Figure 11. Overall component linking model block diagram.
Figure 11. Overall component linking model block diagram.
Applsci 14 02870 g011
Figure 12. Component linking model of Moog valve-controlled system in Matlab Simulink.
Figure 12. Component linking model of Moog valve-controlled system in Matlab Simulink.
Applsci 14 02870 g012
Figure 13. 90 mm/s cylinder velocity simulation results vs. experimental results: (a) Simulink results; (b) FORTRAN simulation results [8]; (c) experimental results [8].
Figure 13. 90 mm/s cylinder velocity simulation results vs. experimental results: (a) Simulink results; (b) FORTRAN simulation results [8]; (c) experimental results [8].
Applsci 14 02870 g013
Figure 14. 90 mm/s load pressure simulation results vs. experimental results: (a) Simulink results; (b) FORTRAN simulation results [8]; (c) experimental results [8].
Figure 14. 90 mm/s load pressure simulation results vs. experimental results: (a) Simulink results; (b) FORTRAN simulation results [8]; (c) experimental results [8].
Applsci 14 02870 g014
Figure 15. (a) Analytical solution and (b) numerical solution results of square wave velocity motion.
Figure 15. (a) Analytical solution and (b) numerical solution results of square wave velocity motion.
Applsci 14 02870 g015
Figure 16. (a) Analytical solution and (b) numerical solution of sawtooth velocity simulation results.
Figure 16. (a) Analytical solution and (b) numerical solution of sawtooth velocity simulation results.
Applsci 14 02870 g016
Figure 17. (a) Numerical solution results and (b) analytical solution results of chamber pressures when valve travel is within underlap region.
Figure 17. (a) Numerical solution results and (b) analytical solution results of chamber pressures when valve travel is within underlap region.
Applsci 14 02870 g017
Table 1. Identifying the root from numerically solved results when the area ratio is 1.9707.
Table 1. Identifying the root from numerically solved results when the area ratio is 1.9707.
Valve Disp Xs (mm) X Factor Root
−0.0045 3.883617722 X 4
−0.004404 3.867534491 X 4
−0.004296 3.848941124 X 4
−0.004101 3.814000308 X 4

−0.0009 2.963117621 X 4
−0.000753 2.911892307 X 4
−0.000735 2.905559937 X 1
0 2.637188916 X 1
0.0012 2.172727301 X 1
0.00255 1.65272018 X 1

0.004101 1.119970977 X 1
0.0042 1.0894822 X 1
0.004302 1.058563937 X 1
0.00441 1.026380398 X 1
0.0045 1 X 1
Table 2. Root distribution of X4 and X1 with different area ratios.
Table 2. Root distribution of X4 and X1 with different area ratios.
Area Ratio γ Percentage of x4Percentage of x1
1.5 46% 54%
1.9707 41.67% 58.33%
2.5 39% 61%
3 37% 63%
4 33.3% 66.7%
Table 3. Simulation time cost of numerical, analytical, and fundamental methods.
Table 3. Simulation time cost of numerical, analytical, and fundamental methods.
Square Wave Simulation
Numerical solved time cost Analytical time cost
5.49 s 1.07 s
Underlap Region Simulation
Numerical solved time cost Analytical time cost
1647.17 s 8.50 s
Sawtooth Command Simulation
Numerical solved time cost Analytical time cost
8.26 s 1.23 s
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, H.; Xu, M.; Cao, Z. Analysis of Underlapped Symmetrically Ported Valve-Controlled Asymmetric Cylinder Drive. Appl. Sci. 2024, 14, 2870. https://doi.org/10.3390/app14072870

AMA Style

Wang H, Xu M, Cao Z. Analysis of Underlapped Symmetrically Ported Valve-Controlled Asymmetric Cylinder Drive. Applied Sciences. 2024; 14(7):2870. https://doi.org/10.3390/app14072870

Chicago/Turabian Style

Wang, Huankun, Man Xu, and Zijian Cao. 2024. "Analysis of Underlapped Symmetrically Ported Valve-Controlled Asymmetric Cylinder Drive" Applied Sciences 14, no. 7: 2870. https://doi.org/10.3390/app14072870

APA Style

Wang, H., Xu, M., & Cao, Z. (2024). Analysis of Underlapped Symmetrically Ported Valve-Controlled Asymmetric Cylinder Drive. Applied Sciences, 14(7), 2870. https://doi.org/10.3390/app14072870

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