Effect of Bulk Modulus On Performance of A Hydrost
Effect of Bulk Modulus On Performance of A Hydrost
Effect of Bulk Modulus On Performance of A Hydrost
net/publication/225921378
CITATIONS READS
52 1,735
1 author:
SEE PROFILE
All content following this page was uploaded by Ali Volkan Akkaya on 15 May 2014.
1. Introduction
543
544 Ali Volkan Akkaya
the hydraulic systems (McCloy & Martin 1980; Watton 1989). Due to temperature variations
and air entrapment, the bulk modulus may vary during the operation of the hydraulic sys-
tems (Eryilmaz & Wilson 2001). A little entrapped air is enough to reduce the bulk modulus
significantly (Merrit 1967; Tan & Sepehri 2002). Moreover, system pressure plays an impor-
tant role on the bulk modulus value (Wu et al 2004). Some effects of instabilities induced by
bulk modulus nonlinearities such as pressure oscillations in the form of pressure waves can
be detrimental to operation of hydraulic systems and may result in reduced component life,
loss of performance, disturbance in control systems, reduced efficiency and increased acous-
tic noise. In spite of these adverse effects, there are few studies about bulk modulus within
hydrostatic transmission systems. Yu et al (1994) developed an on-line parameter identifica-
tion method, determining the effective oil bulk modulus within an actual hydraulic system by
measuring the propagation of a pressure wave through a long pipe. Marning (1997) devel-
oped a linear relation between oil bulk modulus and pressure for a HST system. However, to
date, nothing has appeared in the literature that addresses the effect of bulk modulus dynam-
ics incorporated into a hydrostatic transmission model on control design process of the HST
system. In fact, models of hydrostatic transmission systems with variable bulk modulus have
more complex dynamic behaviour than normal. Moreover, having servo control of the sys-
tem, dynamics of bulk modulus becomes more important because the closed-loop system
itself raises the issue of stability.
Bulk modulus cannot be determined directly and hence needs to be estimated. Based on
this estimation, corrective actions may be taken in control applications for HST systems. The
complex dynamic interactions between variable bulk modulus and the control action is inves-
tigated using modelling and simulation analysis. Simulation tests are particularly beneficial
when preparing a model of a real system is complicated and time-consuming. A servo hydro-
static transmission control system is a good example for this issue. The determination of static
and dynamic behaviours using simulation tests is possible without expensive prototypes. The
simulation also makes a shorter product-designing cycle possible.
This study focuses on control performance of a typical HST system. A nonlinear model
of the system is studied by means of Matlab-Simulink software. The system model is a
combination of each individual component model consisting of pump, valve, hydraulic hose
and hydraulic motor. In addition, the variable bulk modulus is presented to describe the
effects of this phenomenon on system dynamics and control algorithm. For this purpose, two
different hydraulic hose Simulink models are incorporated separately into the system model.
In addition, the models are utilized in the control design process. The control of the angular
velocity of the hydraulic motor coupled with load is achieved by PID (proportional integral
derivative) and fuzzy types of controller. In the first model, bulk modulus is assumed to have
a fixed value and angular velocity control of the HST system is carried out with the classical
PID control algorithm. In the second model, bulk modulus is defined as a variable parameter
depending on entrapped air and system pressure. This new model is applied on velocity control
of the HST system under the same PID control parameters. In the following, fuzzy controller
is implemented in this new model in order to judge its capability against variable bulk modulus
nonlinearity. The simulation results of two control approaches are then compared to analyse
the differences in the performance of the HST system in terms of bulk modulus dynamics.
2. Mathematical model
The physical model of the HST system considered for this study is shown in figure 1. The
variable displacement pump driven by an induction motor supplies hydraulic power to a fixed
Effect of bulk modulus on performance of a transmission control system 545
displacement hydraulic motor for driving load. To protect the system from excessive pressure,
a pressure relief valve is used.
From a research objective point of view, the descriptions of a system mathematical model
should be as simple as possible. At the same time, it must include important characteristics of
the real event. One way to understand the system is to separate the system into components
for the purpose of modelling. Using a fundamental knowledge of physics, for instance the
moment equilibrium and continuity equation, a model that represents the dynamics behaviour
of each component can be derived at the component levels. Having understood each individual
component, we can understand the overall system by interconnecting the components together
to obtain an overall system model (Prasetiawan 2001). In this paper, the model of each
component used for the HST system is developed using earlier methods (Jedrzykiewicz et al
1997, 1998).
where, Qp is pump flow rate (m3 /s), α is displacement angle of swashplate (◦ ), kp is pump
coefficient (m3 /s), ηvp is pump volumetric efficiency (−) which is assumed not to depend on
pump rotation angle.
where, kv is slope coefficient of valve static characteristic (m5 /Ns), P is system pressure (Pa)
and Pv is valve opening pressure (Pa).
2.3a Variable bulk modulus Fluid is an important element of hydrostatic systems and enables
power transmission, hence it can influence the dynamic behaviours of the system and the
control system. The bulk modulus of non-aerated hydraulic oil depends on temperature and
pressure, for mineral oils with additives its value ranges from 1200 to 2000 MPa. Moreover,
system pressure and entrapped air affect the bulk modulus value. If a hydraulic hose is used
rather than a steel pipe, the bulk modulus of this section may be considerably reduced. Owing
to these reasons, the parameters influencing bulk modulus value must be included in the HST
model for more accurate system dynamics.
The equation which gives the variable bulk modulus of fluid-air mixture in a flexible
container is as follows (McCloy & Martin 1980):
1 1 1 Va 1
= + + · , (6)
βv βf βh Vt βa
where, the subcripts a, f and h refer to air, fluid, and hose respectively. It is assumed that the
initial total volume Vt = Vf + Va , and that βf βa . Thus bulk modulus will be less than any
βf , βh , or Vt /Va βa . The bulk modulus of the fluid βf is obtained from the manufacturers
data. The adiabatic bulk modulus used for air is (Cp /Cv )P = 1·4P . With these assumptions,
(6) can be rewritten as in,
1 1 1 s
= + + , (7)
βv βf βh 1·4 · P
where, s is entrapped air percent in the total volume (s = Va /Vt ).
where, kmt is motor torque coefficient (m3 ), P is pressure drop in hydraulic motor (Pa)
and ηmm is mechanical efficiency of hydraulic the motor (−). The torque produced in the
hydraulic motor (Nm) is equal to the sum of the moments from the motor loads and can be
given as,
Mm = MI + MB + Mo , (10)
where, MI , MB and Mo are the moments resulting from load inertia, friction force and machine
operation respectively. These moments can be denoted as
Mm = Im (dω/dt) + Bω + Mo , (11)
where, Im is the inertia of the hydraulic motor shaft (Nms2 ), B is viscous friction coefficient
of motor and its shaft (Ns/m), and ω is angular velocity of motor shaft (1/s). Equation (11)
can be used to determine the angular velocity of the hydraulic motor shaft. This equation is
rearranged for angular velocity as
Qp = Qm + Qc + Qv , (13)
moment:
Mm = MI + MB + Mo . (14)
3. Control applications
Most publications related to the HST control are related to the speed control of the hydraulic
motor connected to the load. In order to achieve this goal, different closed-loop control design
strategies can be used. However, Lee & Wu (1996) showed that using only pump displacement
to regulate load speed is the most effective of all the methods they tested. In addition, Re et al
(1996) concluded that the sole use of pump displacement actuation to control one load speed
of a system with variable-displacement pump and motor is the most efficient, and should be
always preferred whenever possible. For this reason, in the HST systems being considered in
this study, the output angular velocity is controlled by the flow rate supplied to the hydraulic
motor, and this flowrate is adjusted by the swashplate angle of the variable-displacement
pump. Swashplate dynamics are not taken into consideration in the control application in
this study for the sake of simplicity. In addition, the swashplate control system usually has
faster dynamics than the rest of the system, and therefore neglecting its dynamics is justified
(Watton 1989).
To precisely control the angular velocity of the hydraulic motor in hydrostatic transmission
control systems, an appropriate controller must be designed in advance. In industrial appli-
cations, classical control methods such as PI, PID are being used for velocity control of HST
systems. It is crucial to determine controller parameters accurately because PID control meth-
ods have linear characteristics. They are sometimes insufficient to overcome nonlinearities
which exist in the nature of the HST systems for high precision applications (Tikkanen et al
1995; Prasetiawan 2001). In particular, the bulk modulus ought to be regarded as a source of
significant nonlinearity for this type of controller. Thus, the controller has to be very robust
to account for such wide variation. Use of knowledge-based systems in process control is
increasing, especially in the fields of fuzzy control (Tanaka 1996). Unlike classical control
methods, the fuzzy controller is designed with linguistic terms to cope with the nonlineari-
ties. Therefore, this control method is also applied to judge its capacity to reduce the adverse
effect of variable bulk modulus.
determined for the reference angular velocity. Figure 3 shows the Simulink model of the
PID-controlled HST system.
de(t) Kp
uv(t) = Kp · e(t) + Kp · τd · + e(t) · dt, (17)
dt τi
e(t) = ωr − ω. (18)
The inference process consists of forty nine rules driven by the linguistic values of the input
signals. These fuzzy rules written as a rule base are shown given in table 1. The rule base is
developed by heuristics with error in motor angular velocity and derivation of error in this
velocity. For instance, one of the possible rules is: IF e = PS and de = NB THAN u = NM.
This rule can be explained as in the following: If the error is small, angular velocity of hydraulic
motor is around the reference velocity. Significantly big negative value of derivation of error
shows that the motor velocity is rapidly approaching the reference position. Consequently,
controller output should be negative middle to prevent overshoot and to create a brake effect.
As a rule-inference method, the Mamdani Method is selected because of its general acceptance
(Tanaka 1996).
Defuzzification transforms the control linguistic variables into the exact control output. In
defuzzification, the method of centre of gravity is implemented (Tanaka 1996), as
n n
u= Wi Bi / Wi (19)
i=1 i=1
de\e NB NM NS ZE PS PM PB
NB NB NB NB NM NM NS ZE
NM NB NB NM NS NS ZE PS
NS NB NM NS NS ZE PS PM
ZE NM NS NS ZE PS PS PM
PS NM NS ZE PS PS PM PB
PM NS ZE PS PS PM PB PB
PB ZE PS PM PM PB PB PB
where, u is the output signal of the fuzzy controller, Wi is the degree of the firing of the ith
rule, Bi is the centroid of the consequent fuzzy subset of ith rule. Real values of control output
signal (uv) are determined by the scaling factor of Guv. As a result, the fuzzy controller
built-in Fuzzy Logic Toolbox of the Matlab program has been added to the Simulink model
of hydrostatic transmission system for simulation analysis (figure 6).
The validity of the influence of bulk modulus dynamics on HST control system has been
tested in computer simulations. In order to carry out simulation, some physical and simulation
parameters corresponding to HST system are taken from work of McCloy & Martin (1980)
and Jedrzykiewicz et al (1997, 1998), and other control parameters are as given in table 2.
Open loop pressure and angular velocity responses of the HST system are given in figures 7a
and b respectively, under fixed bulk modulus and variable bulk modulus. Comparing the sim-
ulation results shows that the model including variable bulk modulus shows flexible dynam-
ics and decreasing system stiffness (figure 7a). Moreover, a degree of aeration less than 1%
brings about considerable changes of velocity and pressure responses because the aeration of
the working fluid results in decrease of fluid bulk modulus and changes its characteristics as
a function of pressure.
The dynamic behaviour patterns in figure 8 are obtained from PID control. It is observed
from figure 8a that the system model including the fixed bulk modulus is in good agreement
with reference velocity. In contrast to this, the simulation result of the model with variable
bulk modulus show oscillations in the transient regime. The reason is that the bulk modulus
Figure 7. Open loop responses of HST system under fixed and variable (a) system pressure, (b) angular
velocity.
Effect of bulk modulus on performance of a transmission control system 553
Figure 8. PID control responses of HST system for fixed and variable bulk modulus (a) angular
velocity, (b) system pressure, (c) variable bulk modulus.
value becomes lower than that of the fixed one (figure 8b) and changes with system pressure
(figure 8c). The same characteristic is seen for increasing load moment at 0·03 s. The above
results indicate that with change in the bulk modulus, peak pressure as well as fluctuation of
the fluid pressure increase in closed loop control applications. This increases the settling time
of the system’s responses. This may cause to failure of the stability of the system. Therefore,
it is necessary to revise the control parameters or apply a more robust controller in terms of
variable bulk modulus.
Fuzzy and PID control responses of HST system under the variable bulk modulus are
depicted in figure 9. Simulation results of fuzzy controller show good performance compared
with PID controller in tracking referenced velocity (figure 9a). In addition, fuzzy controller
rejects the effect of loading on pressure dynamics (figure 9b). Figure 9c indicates that fuzzy
controller minimizes the fluctuations of bulk modulus value. This is the reason why the
fuzzy response is more robust. Figure 9d shows that PID controller can work safely up to
200 Hz, whereas a fuzzy controller is capable of enduring 400 Hz. Therefore, a fuzzy controller
554 Ali Volkan Akkaya
Figure 9. Fuzzy and PID control responses of HST system under variable bulk modulus. (a) Angular
velocity, (b) system pressure, (c) variable bulk modulus, (d) velocity bode diagram.
5. Conclusion
demonstrated that omitting the bulk modulus dynamics in hydrostatic transmission control
systems may lead to major errors in system response and have implications on the safety of
operation. Therefore, bulk modulus should be considered as a variable parameter to obtain
a more realistic overall model and to determine more accurate control parameters in PID
controller. Analysis including bulk modulus dynamics in an HST-control system model with
this control design feature has not been described in the literature to date. Therefore, it may
be useful for early design of an HST system used for PID control application. In addition, it
is clearly seen that a fuzzy controller has the capability of eliminating the adverse effects of
variable bulk modulus. This will also benefit the control design process in terms of developing
a robust controller. For future research, model development will be expanded to include
swashplate dynamics, valve dynamics and more complex flow and torque models of the
pump and the motor. Furthermore, an adaptive control method will be applied for changeable
velocity reference and load moment.
List of symbols
References
Dasgupta K 2000 Analysis of a hydrostatic transmission system using low speed high torque motor.
Mech. Mach. Theory 35: 1481–1499
Dasgupta K, Chattapadhyay A, Mondal S K 2005 Selection of fire-resistant hydraulic fluids through
system modelling and simulation. Simul. Model. Pract. Theory 13: 1–20
Eryilmaz B, Wilson B H 2001 Improved tracking control of hydraulic systems. Trans. ASME: J. Dyn.
Syst. Meas. Control 123: 457–462
Huhtala K 1996 Modelling of hydrostatic transmission – steady state, linear and nonlinear models.
Acta Polytech. Sci. Me. 123:
Jedrzykiewicz Z, Pluta J, Stojek J 1997 Research on the properties of a hydrostatic transmission for
different efficiency models of its elements. Acta Montanistica Slovaca 2: 373–380
Jedrzykiewicz Z, Pluta J, Stojek J 1998 Application of the Matlab-Simulink package in the simulation
tests on hydrostatic systems. Acta Montanistica Slovaca Rocnik 3: 29–36
Kugi A, Schlacher K, Aitzetmüller H, Hirmann G 2000 Modelling and simulation of a hydrostatic
transmission with variable-displacement pump. Math. Comput. Simul. 53: 409–414
Lee C B, Wu H W 1996 Self-tuning adaptive speed control for hydrostatic transmission systems. Int.
J. Comput. Appl. Technol. 9: 18–33
Lennevi J, Palmberg J O 1995 Application and implementation of LQ design method for the velocity
control of hydrostatic transmissions. Proc. Inst. Mech. Eng., J. Syst. Control Eng. 209: 255–268
Manring N D 1997 The effective fluid bulk modulus within a hydrostatic transmission. Trans. ASME:
J. Dyn. Syst. Meas. Control 119: 462–466
Manring N D, Luecke G R 1998 Modelling and designing a hydrostatic transmission with a fixed-
displacement motor. Trans. ASME: J. Dyn. Syst. Meas. Control 120: 45–49
McCloy D, Martin H R 1980 Control of fluid power, analysis and design (New York: John Wiley &
Sons)
Merrit H E 1967 Hydraulic control systems (New York: John Wiley & Sons)
Ogata K 1990 Modern control engineering (Englewood Chiffs, NJ: Prentice-Hall)
Piotrowska A 2003 The control of the rotational speed of hydraulic engine in hydrostatic transmission
by use of the module DSP. 28th ASR Seminar, Instruments and Control (Ostrava: VS̆B-TU) pp. 291–
297
Prasetiawan E A 2001 Modelling, simulation and control of an earthmoving vehicle powertrain simu-
lator. M Sc thesis, Mechanical Engineering in Graduate College, University of Illinois, Urbana, Il
Re L, Goransson A, Astolfi A 1996 Enhancing hydrostatic gear efficiency through nonlinear optimal
control strategies. Trans. ASME: J. Dyn. Syst. Meas. Control 118: 727–732
Tan H Z, Sepehri N 2002 Parametric fault diagnosis for electrohydraulic cylinder drive units. IEEE
Trans. Ind. Electron. 49: 96–106
Tanaka K 1996 Introduction to fuzzy logic for engineering application (Berlin Springer)
Tikkanen S, Huhtala K, Vilenius M 1995 Fuzzy controllers in hydrostatic transmission. IEE Collo-
quium on Innovative Actuators for Mechatronic Systems (London: Inst. Elec. Eng.) 15/1–15/3
Watton J 1989 Fluid power systems: Modelling, simulation, analog and microcomputer control (Engle-
wood Chiffs, NJ: Prentice-Hall)
Wu K, Zhang Q, Hansen 2004 Modelling and identification of a hydrostatic transmission hardware-
in-the-loop simulator. Int. J. Vehicle Des. 34: 63–75
Yu J, Chen Z, Lu Y 1994 The variation of oil effective bulk modulus with pressure in hydraulic
systems. Trans. ASME: J. Dyn. Syst. Meas. Control 116: 146–150
Zadeh L 1965 Fuzzy sets. Inf. Control 8: 338–353