Hydraulics
Hydraulics
Hydraulics
Control and
Information Sciences
Edited by A.V. Balakrishnan and M.Thoma
33
Peter Dransfield
Springer-Verlag
Berlin Heidelberg New York 1981
Series Editors
h~ V. Balakrishnan - M. Thoma
A d v i s o r y Board
I D. Davisson A. G. J. MacFarlane. H. Kwakernaak
J. I Massey Ya. 7_ Tsypkin A. J. Viterbi
Author
Peter Dransfield, Ph.D.,
Department of Mechanical Engineering
Monash University, Australia
This work is subject to copyright. All rights are reserved, whether the whole
or part of the material is concerned, specifically those of translation, re-
printing, re-use of illustrations, broadcasting, reproduction by photocopying
machine or similar means, and storage in data banks.
Under 54 of the German Copyright Law where copies are made for other
than private use a fee is payable to 'Verwertungsgesellschaft Wort', Munich.
Springer-Verlag Berlin Heidelberg 1981
Printed in Germany
Printing and binding: Beltz Offsetdruck, Hemsbach/Bergstr.
906113020-543210
PREFACE
The Text centres around the ideas that for hydraulic control systems
Among these, the use of bond graphs (abbreviation for power bond graphs) is paramount
to the Text. Therefore, bond graphs are introduced, explained, demonstrated, and
utilized in substantial detail. The features of bond graphs which make them
desirable in the hydraulic control system design situation include their structural
affinity with the real system and its components, the formality in the assembly of a
bond graph structure, their modular and re-usable nature, and the provision of a
structure from which the set of equations which become the model can be formally
prepared.
It is assumed that the reader has a general appreciation of the nature and
operation of hydraulic control systems and their major components. The Text is
written mainly for those mechanical and control engineers who would like to be able
to include predictive dynamic analysis among their techniques for the design of
hydraulic control systems but who find existing approaches unsatisfying or
inappropriate. The Text should be suitable also for up to twenty hours in
appropriate graduate or senior undergraduate courses, and for specialized extension-
type courses.
For full appreciation of the Text, the reader will need to make a "break through"
by thoroughly understanding and using the symbols and structures of bond graphs.
Such a break through was literally forced on me by a graduate student whom I was
supervising for a research Master degree. I gladly acknowledge my debt to
Bevis Barnard, Senior Lecturer at Caulfield Institute of Technology, for the gentle
persuasion emanating from his own discovery and conviction that bond graphs offered
something unique to the modelling, simulation and design of powered control systems.
I acknowledge also the conceptions of Henry Paynter, the original and sustained
developmental work of Dean Karnopp and Ron Rosenburg, the work of bond graph
entrepreneur Jean Thoma, the assistance of colleague Jacek Stecki, and the
contributing work of a small but steady stream of undergraduate and graduate students
including Rob Winton, M.K. Teo, Roger LaBrooy, Kevin Duke, S. Ramachandran, V.K.L. Mai,
and Merren Cliff.
Many people contribute to the production of a text. I single out two for
special thanks and appreciation; Lorry Ryan and John Millar, both of Monash
University, for their excellent work in typing and drafting respectively.
And there is an individual without whose encouragement this Text would not have
been completed my wife Nevillie whose support was sustained, inspirational, and
essential.
To those willing to learn, and to those willing to try a new way, welcome aboard.
CHAPTER 1 INTRODUCTION 1
1.1 Introduction 1
1.2 A Design Philosophy 2
1.3 Dynamic Modelling 7
1.4 A Desirable Background ii
1.5 Arrangement of Text ii
1.6 Conclusion 12
CHAPTER 7 (contd.)
7.9 4-Way Control Valves 125
7.9.1 Introduction 125
7.9.2 The Basic Relationship 126
7.9.3 Valve Flow Nomenclature 129
7.9.4 The Closed-Centre Control Valve 129
7.9.5 The Open-Centre Control Valve 132
7.9.6 The Tandem-Centre Control Valve 136
7.9.7 Summary 138
7.10 Actuators 139
7.10.1 Introduction 139
7.10.2 Linear Actuator 139
7.10.3 Rotary Actuator (Hydraulic Motor) 142
7.10.4 Equations and Coefficients 142
7.11 Some Other Common Components 144
7.11.1 Hydraulic Lines 144
7.11.2 Filters 146
7.11.3 Accumulators 146
7.11.4 Some Other Valves (Check Valve, 148
Counterbalance Valve)
7.12 Conclusion 151
7.13 References 152
CHAPTER i0 (contd.)
10.2 cavitation 197
i0.2.1 General Discussion 197
10.2.2 Cavitation in Modelling and Simulation 200
10.3 Hydraulic Backlash 201
10.4 F l o w Forces in V a l v e s 204
i0.5 Hydraulic L o c k 205
i0.6 Contaminated Fluid 210
i0.7 Conclusion 212
10.8 References 212
C H A P T E R ii CONCLUSION 215
INDEX 226
CHAPTER 1
INTRODUCTION
i.i INTRODUCTION
Hydraulic control systems are used to control the position or speed of resisting
loads. Final drive is usually provided by a hydraulic actuator either a linear-
motion hydraulic cylinder or a rotary-motion hydraulic motor. The actuator develops
its force or torque by receiving liquid from a positive displacement pump at a
relatively high pressure, usually between 7 and 35 MPa. The actuator develops its
motion by receiving flowrate of liquid from the pumping system. Fig. i.i illustrates
several simple hydraulic control systems.
I I HydraulicCylinder7
i .
~-----~ilter I ~ InertiaLad~
I~ ~ ~ %---4-WayManualControlValve
II ~-Pressure
ReliefValve
I] II Reservoir
~NN
f
i:
It is not sufficient for the designer of a hydraulic system to know that his
proposed system WiZ~ move the driven load from one state to another. He should know
also how the load will move. He needs to appreciate not only the initial and final
states of the responses, but also the time-domain path between these states (Fig. 1.2).
He should know if system response is stable, if it is fast enough or too fast, if it
is oscillatory or aperiodic, etc.
corrective in nature, the objective being to find out why an existing system
is not performing satisfactorily and to deduce what can be done to improve the
situation,
In the latter situation, the proposed system's dynamics are designed as part of the
system design process. The present Text is mostly concerned with this approach,
though it is relevant also to corrective analyses.
Fig. 1.3 illustrates the stages in a dynamic design philosophy. Specify Task
does not imply use of any particular system or technology - - simply a definition of
the task to be accomplished. Define Required Response is also a consideration
independent of the particular system or technology to be employed ....... typically this
would be a time-based curve of the desired system response form, with time and response
magnitudes as explicit scales, such as is shown on Fig. 1.2.
past experience,
the maximum pressure to be allowed,
static power and overall motion considerations,
the available hardware components, and their cost.
a. Linear Drive b. Rotary Drive
200~
time, s I time, s
Close(
!
0 0.5 ~v 0 0.5
i m. 50~
0.5m-
me, s . . . . . . . Re; s
0 200
0 0.5 1 0 0.5
Dynamic performance analyses are often carried out for the high-cost hydraulic
systems associated with aircraft, military, space, and some machine tool applications.
They are less often associated with the design of industrial systems. One result is
that while many such systems are accepted by the user, their performance rarely
approaches a considered optimum.
Relatively few mechanical (or other') engineers are confident and competent with
predictive dynamic analysis of systems. Predictive dynamic analysis hinges on the
formation of a mathematical representation of the system. ~is mathematical model
is used to extract information on the response characteristics of the system. If an
adequate mathematical model can't be formed quickly and with confidence, dynamic
analysis is not an attractive proposition to most designer-analysts.
Few designers of hydraulic control systems are able to confidently form dynamic
models of proposed systems. On the other hand, given a mathematical model, most
engineers can arrange for its solution, which provides the required prediction of
dynamic performance. Generally, solution of anything other than the simplest models
is best effected by digital computer. Digital simulation procedures, widely
available in easy-to-use forms, are of particular relevance. There is no need for
the designer-analyst to be himself a computer programmer. Given the means to solve
(simulate) his model, the designer-analyst can change one or more of the model's
parameters in order to seek a desired form of system response. Such optimization of
response can be sought through repetitive solutions of the model as parameters are
changed. Automated optimization procedures can also be used.
~'r-t
0 ~ t2 tr
Input Expected Response R e q u i r e d
~ 3. Propose System, I
and its Hardware, I
to Accomplish 1
Task I
4- Develop Dynamic
Model
I l
~ ~esponso L 7. Judge Response I J9 Promising;
Response I
I Hopeless.
Redesign System I TM
Quality; Decide 'i
Between 8, 9, I0 vI Optimize it
Acceptable.
Finalize Design
transfer functions,
block diagrams and signal flow graphs,
sets of equations not having a unique formal arrangement,
vector-matrix (state space) expressions,
have been around for a long time. Journals contain many examples of successful
applications of each of them. Yet the lack of general use of these model forms by
designers of hydraulic systems suggests that they are insufficiently attractive.
the need for reliable re-useable models of common hydraulic system components.
Hydraulic control systems are of modular construction, being made up by
assembling and interconnecting a set of well-defined off-the-shelf hardware
devices. A particular type and size of pump, valve, actuator, etc. is used
in many different systems. The designer needs to have available a set of
dynamic models, one for each common hydraulic system component, which is
capable of being interconnected to form a system model in a manner reflecting
the formation of the hardware system itself,
the need to be able to change easily the model of a proposed system to reflect
conceivable variations in the proposed system, aimed at ensuring choice of the
best design for the particular system.
They are inherently non-linear due to the presence of turbulent orifice flow,
friction, valve lap, pressure limiting actions, and other effects;
They usually are expected to provide large scale response. Quite normally
they are commanded to move the controlled load from one extreme position to
the other;
Power flow modelling, and in particular its recent extension POWER BOND GRAPHS,
provides the designer of powered control systems with a dynamic modelling procedure
which is more fundamental to, and more useable in, the system design environment.
It is intended to demonstrate this by the end of the Text. It is intended also to
leave the reader with a library of power bond graph models of the more common hydraulic
system components, and with the ability to develop power bond graph models of any
additional components he might wish to use.
Fig. 1.4 (a) illustrates a simple hydraulic system, Fig- 1.4 (b) a power bond graph
structure for the system, and Fig. 1.4 (c) the set of equations which is a dynamic
model of the system. The bond graph structure is the key for linking the system to
its dynamic model. The bond graph is a type of signal flow diagram which to a person
familiar with bond graphs bears a striking affiliation with the hardware system
diagram. The set of equations is readily formed from the bond graph structure.
Each line in the horizontal chain of Fig. 1.4(b) represents a physical component
through which power flows between the adjacent system components. The product of
a. The S y s tern
A ....... Fa Ira
r
I / /,t,,,,,ILr
............ ~
I Friction
Ca Im
C. Dynamic Equations
Ps = constant (assumption)
Qs = K - Xi - ~Pv ~ (for Xv positive)
Qc = Os - Q a
the two variables written against each line describes the power flowing in the line.
Thus Ps - Qs (supply pressure times supply flowrate) describes the power flowing
into the control valve. Similarly, Fa - Xm (force applied to the mass times
velocity of the mass) describes the power flowing into the driven load. The R ,
C , and I symbols appearing at the ends of the spur lines of the diagram indicate
respectively the presence of resistive (R) , capacitive (C) , and inertive (I)
effects which affect the system's dynamic response. The half-arrows on each line
indicate the power flow direction in each power bond, either real or assumed. The
short transverse bar at one end of each line indicates a cause-effect decision aimed
at arranging equations, to be formed later from the bond graph structure, in a form
most suitable for simulation. These various points will be developed and explained
in the Text. The Text's nomenclature and units are discussed in Chapter 2.
Formation of the model structure, i.e. of the power bond graph, decisions on the
precise relationships to be used in the equations, and formation of and arrangement
of the equation set, are separate and progressive steps in the dynamic model derivation
procedure. In conventional modelling procedures these steps become jumbled together,
often creating confusion in the minds of part-time dynamic analysts. It will be
demonstrated in later chapters that power flow models allow the formality and
flexibility required in the system design situation, particularly in that
component models are useable in a variety of system models, as are the hardware
components themselves useable in a wide variety of systems,
Power flow and power bond graph modelling procedures can be used with all types
of powered control s y s t e m s - - e l e c t r i c a l , pneumatic, mechanical, etc. The present
Text however concentrates on their use with hydraulic systems, although this
encompasses electric motors, electrohydraulic devices, and a variety of mechanical
sub-systems. The ideas of power flow and power bond graph modelling have been
developed by others, and references will be provided in appropriate places. The
present Text simply utilizes power flow techniques and orientates them specifically
towards use in the design and analysis of hydraulic control systems.
Beyond a review in Appendix i, with References, the present Text does not provide
instruction on how to develop a basic design for a hydraulic system for a specific
I!
familiarity with basic hydraulic control system components, systems, and the
symbols used to represent hydraulic components in system diagrams of the type
of Fig. i.i,
Chapter 1 provides the setting for the Text and makes its case;
Chapter 5 introduces power bond graphs, the key feature of the Text, which
is offered as a major step towards making predictive dynamic analysis viable
for more designers of hydraulic control systems;
Chapter I0 briefly introduces some factors which can affect the performance
of hydraulic systems, and about which the system designer should be aware;
1.6 CONCLUSION
In conclusion, the orientation attempted with the Text can be stated: "Given a
proposed hydraulic control system, simulate it to check and to optimize its dynamic
performance before any final commitment is made to expensive hardware".
CHAPTER 2
2.1 INTRODUCTION
What a help if everyone would use the same notation and units' This would allow
immediate reCognition of symbols, and of the significance of the magnitudes of
coefficients and numbers appearing in relationships. The value of this kind of
recognition in the system design environment should not be under-estimated. The
designer needs to be able to quickly utilize established techniques and data, rather
than seek out new unproven and possibly ambiguous materlal.
The notation, the way of writing relationships, and the units used in this Text
seek standardization.
2.2 NOTATION
2.2 .I Symbols
Upper case English letters are used generally for system variables and
coefficients, each being reserved for particular quantities; some widely-
14
A -- a r e a a - - angular acceleration
B - bulk modulus of oil b - B-
C - capacitance, or capacitive c - y -
effect, or compliance d- -
D - operator d/dt e - E -
E - general symbol for effort f- - damping ratio
(i.e. potential) variable; g - gravitational - efficiency
also voltage acceleration 8 - angular displacement
F - force h - I
G - i - electrical K
GY - gyration power transformation current
H - j- - coefficient of
I - inertive or inductive effect; k - friction
~ -
inertia or inductance v - kinematic viscosity
J - m -
C u r r e n t l y a c c e p t e d u s a g e and p h y s i c a l a s s o c i a t i o n o f symbols w i t h p a r t i c u l a r
q u a n t i t i e s are r e t a i n e d as far as is feasible;
8a d e n o t e s the a n g u l a r d i s p l a c e m e n t o f a h y d r a u l i c m o t o r (actuator).
2.2.3 Operations
M u l t i p l i c a t i o n is i n d i c a t e d b y a central dot: N P ; Np Ck ; 13 Xm ;
J J
2.2.4 Equations
QEp = Pp/REp
Qp = Kv Xi dp
= ~m(0} + z / m - I~ dt
2.2.5 Comments
A little usage will show that this scheme is easy to write, easy to type,
easy to understand, and is easily associated with digital computer programming
and print out.
Some of the terms included in Fig. 2.1 apply specifically to power bond graph
terminology to be introduced later.
Only symbols used in this Text are designated in Fig. 2.1. There are plenty
of gaps for future assignment or general use.
2.3 UNITS
Syst~mes Internationales (SI) units are used without exception. Apart from its
substantial and growing international acceptance, the SI system overcomes most of the
ambiguity associated with the use in engineering and science of the British and the
standard metric systems of units.
The SI units which appear in this Text are shown in Fig. 2.2. Appendix 2
Conversational
Physical Quantity Name of Unit Symbol of Unit
Alternatives
Basic Units
Length metre m centimetre (cm) ,
millimetre (mm)
Mass kilogram kg tonne (E3 kg)
Time second s minute s
Electric current ampere A
Temperature degree celsius C
Derived Units
Area metre 2 m2
Volume metre 3 m3 litre (E-3 m 3)
Density kilogram/me tre 3 kg/m 3
Force (also weight) newton N kilonewton (kN)
Torque newton metre N - m
Rotating inertia kilogram metre 2 kg - m 2
Pressure pascal (N/m 2) Pa kPa, MPa, bar(E5 Pa)
Energy (work) joule J
Power watt W kilowatt (kW)
Electric potential volts E
(electric voltage)
Electric resistance ohm
Velocity metre/second m/s
Acceleration metre/second 2 m/s 2
Flowrate metre 3/second m3/s litre/sec or
litre/minute
Angle radian rad revolution
Angular velocity radian/second rad/s revolution/sec or
minute
Frequency hertz Hz cycle per second
Dynamic viscosity newton second/metre 2 N s/m 2, centipoise (cP)
or pascal second or Pa o s
Kinematic viscosity metre2/second m2/s centistoke (cSt)
Common multiples: micro (u) E-6 ; milli (m) E-3 ; centi (c) E-2 ;
kilo (k) E3 ; mega (M) E6
2.3 .I Comments
or w r i t e descriptively of i0 M P a , b u t should w r i t e it as i0 E6 P a w h e n it
comes to equations.
kg -m 2 , not k g m 2 , for k i l o g r a m m e t r e 2.
2.4 REFERENCES
2.4 The A B C ' s o f SI, National Fluid Power Association, USA, 1974.
CHAPTER 3
3.1 INTRODUCTION
Transfer functions,
block diagrams and signal flow graphs,
vector-matrix (state-space) equations,
informally arranged sets of dynamic and algebraic equation,
are regarded in this Text as conventional forms for dynamic models. They are
conventional in that each is widely used in hydraulic system analyses readily
available in texts and articles.
a. System Schematic
Xi
supply Oil
const, press, ps
XV
.......... , ......
Kvl, Kv2--
\\ L2
Xm
Rf
xv = (L2 - L I ) / L 2 Xi - L I / L 2 Xm 3.1
Qs = Kvl " Xv - Kv2 - Pa 3.2
Qs = A - ~m+v/B ~a 3.3
Im ~m = Pa " A - Rf - Xm - Ck - Xm 3.4
c. Transfer Function
X m = K1
Xi K2 D3 + K3 D2 K4 D + 1
e. Block Diagram
I D ...........
B~~/(V-D) -A
Xi'-- ( L 2 - L I I / L 2 . Kvl ~ __
-LI/L2
FIG. 3.1 SOME CONVENTIONAL MODEL FORMS FOR SYSTEM OF FIG. 3.1a
21
Xi as the controlling input and load position Xm as the controlled output. Supply
pressure Ps is assumed to remain constant. The lever system L1 , L2 is such as
to allow small control valve motions Xv and much larger load motions Xm . Xi is
a m a l l motion. This system will be used to demonstrate conventional model forms.
Depending on the analyst the equations 3.1 to 3.4 could be differently arranged,
or even different in content. The present equations are all linear, whereas several
non-linearities can readily be included if the analyst considers them to be important
to his response predictions.
All approaches to modelling require insight of the physical laws and relationships
which govern the b e h a v i o u r o f the various devices used in the system components,
together with an appreciation of the interaction of the components and of the structure
of the system. The successful formation of dynamic models is an art as much as a
science, and experience and experiment are important assets.
2. Select that response variable (or set of variables) study of which will
provide clearest appreciation of system response;
3. Designate an input (or set of inputs) which will cause the system to respond
in a realistic operational mode. It is not necessary to specify the form
of the input at this stage, merely its physical variable and location;
4. List ass~nptions which are valid in the particular circumstances, and which
will render modelling easier but still adequately valid. Typical
assumptions which might be valid include: input drive speed constant;
pressure relief valve does not lift; all driven inertias can be lumped in
with the designated load inertia; all friction is viscous; pressure in the
return lines is negligibly small (zero); motor-load coupling is rigid; the
model is for small scale response predictions only; etc. This is an area
requiring considerable thought for each modelling attempt. However, there
is a lot of confirmed past experience available in the literature;
3.4 S E L E C T I O N OF R E L A T I O N S H I P S
V = V(0) + A X m 3.7
3.5 S E L E C T I O N OF P A R A M E T E R V A L U E S
decided to regard the input drive speed to the pump (~p) , as constant
(i.e. as a source). This eliminates primary drive (electric motor and
shaft to pump) dynamics, and the inertia of the pump rotating assembly,
from the analysis;
decided to lump the internal leakages (case drains) of pump and motor and
account for them via the leakage coefficient KZ ;
decided that the resistance of the short hydraulic lines joining pump and
motor do not significantly affect the response ~j(t) ; also that the
return line pressure is so low relative to the drive pressure Pa that it
can be neglected for the dynamic analysis;
lumped the hydraulic motor and its shaft inertias in with the load inertia
~o give the total Ij ;
Qp = Kq up #i 3.8
The volume flowrate of oil effectively used at the hydraulic motor, Qa m3/s ,
is given by,
Qa = Qp- Q~ - ~ c
=
Qp - K " Pa - V/B Pa 3.9
/t , Lino / "
~D'"Pump Make~-up ZHyd , Motor /__
v~ m 3 Ka, rad/m 3
*~' K, N.m/Pa Drive Friction
Rf, N.m
Qa
B is b u l k m o d u l u s o f oil, Pa .
Motor s p e e d mj due to Qa is g i v e n b y
ej = Ka Qa 3 .i0
The m o t o r p r e s s u r e - t o r q u e r e l a t i o n s h i p is
La = K Pa 3 .ii
A p p l y i n g N e w t o n ' s l a w to the d r i v e n i n e r t i a ,
L a - Rf ~j = Ij ~9 3.12
~_i
~i = K2 " D z + K3 D + i 3.13
3.15
and a corresponding change in the block diagram of Fig. 3.2 b . This formulation
makes sense if swashplate angle ~i is controlled from neutral, through a desired
input cycle, and back to neutral, causing the load to move from one equilibrium
position to another.
3.7 CONCLUSION
3.8 REFERENCES
3.1 Blackburn, Reetof, and Shearer, Fluid Power Control, Tech. Press M.I.T.,
1960.
3.2 Lewis and Stern, Design of Hydraulic Control Systems, McGraw-Hill, 1962.
3.4 Walters, R., Hydraulic and Electrohydraulic Servo Systems, Iliffe, 1967.
3.6 Aizerman, M., Pneumatic and Hydraulic Control Systems, Pergamon, 1968.
3.8 Kellar, G.R., Hydraulic System Analysis, Indust. Pub. Co., 1970.
3.9 McCloy, D. and Martin, H., Control of Fluid Power, Longman, 1973.
4.1 INTRODUCTION
to discuss the case for modelling hydraulic control systems via study of the
instant (dynamic) power flows taking place between the components of the
systems;
Those readers familiar with these concepts may prefer to pass over or scan this
Chapter and move directly into power bond graph formulation. Those for whom power
flow modelling is relatively unknown should study this Chapter as an introduction to
the more sophisticated power bend graph material.
Power
i Electric PowerTransformers Power
Absorbing
Supply Load
Eieetric Hydraulic
Pum
Hydraulic
C~
I
Inertia
Load
IElec.Power,~EM,Mech.P o w ~ H v d . PowerL, HC ,Mech-
l(Hyd'iLine~
I
PowerS"1
[ I IPwerLines)~I(Shaft)-I
PowerTransmitters
energy. A low energy-utilization rate means a slow motion imparted to the load; a
higher energy-utilization rate causes faster motion. Thus, the dynamic response of
a hydraulic control system is dictated by the dynamic energy flows (i.e. energy flow-
rates) taking place.
Accepting the idea of power flow, it follows that system dynamic response is
dictated by dynamic power flow. From this arises the concept that dynamic response
should be predicted through study of dynamic power flow in the system. This idea is
reinforced by realization that a hydraulic system ~ o e 8 function via dynamic power
interchange between its hardware components. The need then becomes to be able to
describe the instant power flows taking place between the assembly of components
which together comprise the system. The components of the system of Fig. 4.1 can be
considered to be a power source (the mains supply); power transformers (the electric
motor which accepts electric power and puts out mechanical power, the hydraulic pump
which transforms mechanical input power to hydraulic power, the power cylinder which
accepts hydraulic power and puts out mechanical power to the load); power transmitters
(the shafting connecting pump to electric motor, the hydraulic line connecting cylinder
to pump); and a power absorber (the driven load).
Further, hydraulic control systems are controlled by very low powered manual or
automatic command signals, and can consequently be regarded as power amplifying in
addition to power transforming and power transmitting systems.
Let us accept the natural prima facie case for investigating the dynamic response
of hydraulic control systems through study of the power flows which take place
following a command input to the system. The benefits of doing this will emerge in
32
Each component of a hydraulic control system has only a limited number of ways
through which power can be taken in or given out. Neglecting all losses for the
moment, an electric motor can take in power from the mains supply, and can give it
out to its output shaft. That is, it has only two physical power flow connection
points. These are known as power ports, or simply as ports. Fig. 4.2a illustrates
the motor as a schematic, and as a block having two power ports (column 3).
A hydraulic accumulator has only one port through which it can accept or give
out power. It is represented in Fig. 4.2b as a schematic, and as a block having a
single power port.
A hydraulic pump can also be regarded as essentially a 2-port device, the ports
being the drive shaft and the discharge port. However, the internal leakage (case
drain) of a pump affects its power output and consequently the dynamic response of
the system. It should be accounted for in the power port structure as a third port,
through which power is lost, as illustrated on Fig. 4.2c, column 3.
The 4-way control valve is a 4-port device, with power able to flow into and out
from it via the supply port, the two control ports, and the exhaust port. Fig. 4.2d
illustrates this.
A well designed and well placed system reservoir (tank) rarely affects a
system's dynamic response. However, it receives power from the pump leakage, and
from the return line from the control valve. A reservoir can be shown as a 2-port
device receiving these two power flows (Fig. 4.2h).
Power Port Schematic
Co.,~onent
Ideal With Losses
a. Induction
elec. mech.
Electric Motor
i:~wer--~power
in out
T, friction, windage
b. Hydraulic
~ thermal
Accumulator
loss
hydraulicpower,
in and out
c. Hydraulic Pump pump friction,
mech. ~
power - - ~
in ~
hyd.
power
out
~ nert/a
losses
case drain
leakage
d. Control Valve control
supply-~ port 1
tank control
port 2
g. Hydraulic
Cylinder supply-~--~loa d
tank--4.-~
friction,
entry losses
h. Reservoir return flows
I I
i. Inertia Load friction loss
me~.
wer--D
in
A driven load of the kind illustrated on Fig. 4.1 is ideally a 1-port component,
this being the mechanical connection whereby the actuator is connected to the load
(Fig. 4.2i).
Basic power port designations are readily expanded to include loss and other
effects which are thought, in a particular case, to affect a system's dynamic
response. For example, an electric induction motor loses power to speed slip and
should be regarded as a 3-port component. It may also lose significant power to
internal friction and windage and could then be regarded as a 4-port component
(Fig. 4.2a, column 4). Similarly, the n-port designation of a hydraulic pump
could include an environmental port allowing for friction loss (Fig. 4.2c, column 4).
Fig. 4.2, column 4 includes some power loss allowances for other devices. The
designer-analyst should not feel compelled to use the most complete power pert
representation. Rather, he should include only those power flows which he feels
will significantly affect the dynamic response of his particular system. The
simplest model adequate for the purpose in hand should always be sought.
Power port blocks can be used to form a power flow schematic of a system.
Simply couple the blocks as adjacent components are themselves coupled. Fig. 4.3
shows a power flow schematic for the system of Fig. 4.1, expanded to include a
control valve and the reservoir, both of which manipulate or receive power in this
system. A power port schematic provides a structure from which the set of equations
which is the model of the system can be prepared.
considering the power flowing into and from each component to be used in the
system, and forming a set of power port blocks for the components. These
component blocks are quite reuseable;
simple joining of the component blocks to form the system power flow
schematic. Such joining is directly analogous to physical connection of
the hardware components to form the system. For example, the power output
end of the electric motor power port is made one with the power input end of
the hydraulic pump block. The line joining the two blocks is precisely
analogous to the shaft which connects the two physical components. In both
cases the connection is a power bond.
35
a. Detailed S c h e m a t i c
constant
Voltage
Source
b. Shaft C o m p l i a n c e R e c o g n i z e d
EM -F---
..... Iv
C. Pump Speed A s s u m e d C o n s t a n t
Constant
Speed Line ~ m m
Source
~I Reservoir
d. C o n s t a n t P r e s s u r e Supply to C o n t r o l V a l v e
|
Line
Control
HC Im
Sourj Valve
~ Hyd.
- I sine I~ I
llj
FIG. 4.3 POWER PORT SCHEMATICS
38
The power port schematic is readily expanded to allow for effects thought to
affect system response. For example, if it is required to include the compliance of
the shaft between electric motor and pump, place a block labelled SHAFT in the line
joining the motor and pump blocks (Fig. 4.3b).
Mechanical power flow is Net force (or torque) Velocity (or angular
velocity)
Note that
Fig. 4.4a illustrates a simple hydraulic control system; Fig. 4.4b includes the
37
power flow product variables, with the potential variables above and the flowrate
variables below the component symbols. Arrows indicating the natural flow direction
of power are included - - p o w e r is drawn from the electric mains supply (source of
voltage) and is delivered to the driven load- Neglecting all losses, power is
conserved and
The relationship between the potential and the flowrate variables of a particular
power product is of a cause-effect form. For example, force is the cause w h i c h
induces velocity of the load mass. Voltage applied to an electric motor is a cause,
and the resulting current drawn from the supply is an effect induced by the load
impedance seen by the motor.
At first sight it might appear that the potential variable is always cause and
the flowrate variable is always effect. This is not the case. It is not always
obvious which variable of a power product is cause, and which is effect. Appreciation
of cause-effect is a great aid to the understanding of a system's dynamic behaviour.
It will be demonstrated in later Chapters that assignment of causality (i.e. decision
on which variable of a power product is cause and which is effect) can be a key link
in obtaining from a power flow diagram the correct set of dynamic equations for
simulation.
Consider the system of Fig. 4.4a and b. It is proposed to decide which variable
in each power product is cause and which is effect. We shall rearrange Fig. 4.4b to
show cause variables above and effect variables below, rather than the Fig. 4.4b
arrangement of potential variables above and flowrate variables below. Further, we
shall designate the cause variable by an arrow in the direction of desired power flow,
and the effect variable by an arrow directed against the power flow. This is a
convention. However, it has physical significance, as the effect variables are back
effect8 dictated by the dynamic nature of the system. For example, applied voltage
is a source or constant and can only be designated as a cause as its magnitude is not
affected by the load effects of the system. Designation of voltage as cause
automatically requires current to be designated as an effect. However, this is also
physically the case. The current drawn by the electric motor depends on the load
applied to the motor by t~9 8yst~n, i.eo the current i8 a back effect:
This leaves cause-effect decisions to be made for the hydraulic pump and for the
38
a. Power Flow
Potential Variables
/ \
Voltage Torque Pressure Force
HC ~ Im
I Shaft - Hyd. Line-- Cyl. Rod
\ ...................
Flowrate Variables
Power I HP HC ...... ~ Im
Blow ~ EM
...i
.
\ /
Effect Variables
Check to see that physical sense is not violated. Voltage is applied to the
electric motor, whose motion is resisted by the system load; this requires current
to be drawn from the supply main. Resultant speed of the motor is transmitted
directly to the hydraulic pump, which, being of positive displacement, causes positive
flowrate of liquid to the cylinder. Liquid pumped positively into the hydraulic
cylinder, whose piston motion is resisted by the inertia load, creates rapidly
increasing pressure. This pressure acts on the piston to generate a force, which
acts on the inertia load. The inertia load is put into motion. The magnitude of
the inertia load dictates its acceleration and consequently its velocity, which is
also the velocity of the piston. Piston velocity affects the rate of pressure
change. Pressure dictates the torque required to drive the pump and thus the torque
output required of the electric motor. The torque demanded of the motor dictates
the current drawn by it. These actions do not take place in complete time sequence
as implied here. They are occurring simultaneously, but are out of phase in the
time domain. That is, it takes time for an action to fully develop following an
applied cause. Predicting the sum total of these time effects is what dynamic
analysis is all about.
Note that the effect variable does not necessarily result directly from its own
cause variable. In the present example, this correlation occurs only for the inertia
load, for which velocity is the effect of the causing force. Certainly current is
not caused directly by the applied voltage; it is dictated by the shaft torque
demand.
So far, we have reached the stage where a power flow schematic can be proposed.
Thinking this out and making the associated decisions provides insight into those
aspects of its components which will affect a system's dynamic performance. The
ultimate objective remains the production of a set of e Uuations (the dynamic model)
which can be ~imulated (solved) to obtain the predicted behaviour of the significant
system variables. It will become apparent that a power flow schematic can provide
the structure from which a complete and suitably arranged set of equations can be
written.
For example, consider the simple situation of Fig. 4.5a, for which supply
pressure Ps is regarded as constant (a source). When the variable position 4-way
closed-centre control valve is actuated by the displacement input signal Xv , the
system will respond dynamically to move the inertia load. Fig. 4 . 5 b s h o w s the
initial power port structure proposed to describe dynamic power flow following
application of Xv . Hydraulic line resistance has been neglected. If power loss
in the return line from cylinder through valve to tank is considered to be negligible
as far as load dynamic response is concerned, the power flow schematic can be
simplified to that of Fig. 4o5C. Each power bond (line) on Fig. 4.5c is allocated
its power product variables, as follows:
~.~m describes power absorbed by the load (i.e. flowing from cylinder to
load.
Fm is the net force applied to the driven load, and ~ the velocity
developed by the load.
Also, we propose to account for a power loss due to friction associated with the
41
a. The System
Valve
Friction
Source
Ps ~ Cont. i
Valve |~
iTM
Sourcel P~
Ps I Qs
driven load, and which is described by Ff " ~ , where Ff is the force required to
overcome friction.
Being satisfied that the power flow schematic of Fig. 4.5c is adequate for our
purpose, we now require to prepare from it an appropriate and complete set of
equations. This requires an equation for each power variable, expressed in terms of
only other system power variables and system coefficients. Further equations which
might be required to describe variables other than these po~Jez~ 8tatg UGI~GbZe8 can be
added if desired. Environmental power loss effects, such as load friction in the
present case, can be described as separate equations, or can be built into the power
state equations, whichever the analyst may prefer. For the present case, we can
write the six power state equations:
to specify any significant response variables other than the power state
variables,
xm -- I ~ at + Xm(O) 4.7
V = V(O) + A Xm 4.8
( V is volume of oil under compression)
Ff = to be decided 4.9
(friction force)
43
Xv = to b e s p e c i f i e d 4 .I0
(control valve d i s p l a c e m e n t input)
Ff = Kf Xm 4.9a
Ff = s g n Xm - Fu + Kf Xm 4.9b
4.5 AN EXAMPLE
Fig. 4.6a illustrates a hydrostatic drive system whose dynamic response and
performance characteristics are required prior to purchasing the components and
constructing the system. Fig. 4.6b shows a power port diagram prepared for the
system. Each line between blocks is a power bond, representing a physical power-
transmitting connection. Each line has two power state variables associated with
it, their product describing the power flowing. At first sight it appears that
there will be 2n power state variables, where n is the number of power bonds
present. However some bonds share a power state variable. For example, Pp is
the pressure on the bond between pump and relief valve, and also on the bond between
relief valve and filter. Also, ~a is the velocity (flowrate variable) on both the
hydraulic motor-load bond and the load friction power dissipation bond. Examination
of the power port diagram reveals the following set of power state variables, starting
from the assumptions that the electric motor runs at constant speed, and that the
charge pump power requirement does not affgct the dynamics of the main system:
a. Circuit Diagram
~'tric Motor
Pressure Compensated Pump
From Tank
b. Power Port
IElectric I
Motor I
T S~ Tfa ~a
Pal pressure developed in line between control valve and hydraulic motor
Qvl flowrate from control valve towards motor
Qal flowrate received by hydraulic motor
Pa2 pressure developed in return line between motor and control valve
Qa2 flowrate from motor towards control valve
Qv2 return flowrate entering control valve
Qe discharge flowrate from control valve to reservoir
Qre flowrate from crossline relief block to reservoir
Prc pressure in charge line to crossline relief block
Qrc flowrate in charge line to crossline relief block
Qia leakage flowrate from hydraulic motor
Ta torque developed in motor-load shaft
~a angular velocity of motor and load
Tfa friction torque
To form the model, an equation for each of these 24 power state variables is
required. If numbers can be provided for the various coefficients, the set of
equations can be simulated to provide the sought-after dynamic response information.
The system can be controlled or disturbed by adjusting the state of the variable
delivery pump, by adjusting the 4-way control valve, or by adjusting load inertia or
friction. This example and its equations are developed more fully via bond graphs
in Chapter 8.
4.6 SUMMARY
Power can enter or leave a component through only a small number of readily
identified and physically meaningful power ports. Adjacent components are
connected together functionally only via these power ports.
Each power variable is a power state variable which is the object of one of the
set of power state equations which comprise the basic dynamic model of the system.
Using preferred causality, the power state equations can be arranged in algebraic
and integral forms suitable for computer solution or simulation.
System power port structures are modular and closely resemble the modular
construction of the hydraulic system itself. This allows components to be
readily changed, added, or deleted from a system's power port block diagram.
Dynamic coupling effects are inherently included in mQdels derived from power
port structures.
The development of power port models will not be pursued further, as the Power
Bond Graph techniques to be developed in the next Chapters are a more advanced and
formal development of power flow modelling.
4.7 REFERENCES
4.1 Unruh, D., A Standard Format for Mathematical Models of Fluid Power
Systems, Proc. N.C.F.P., V26, 1972.
4.3 Stecki, J.S. and Cuthbertson, S.C., The Effects of Changeover to Fire
Resistant Fluilds on Hydraulic Performance - Part 2: System Criteria,
FRH Jnl., VI, N2, pp 167-173, FPRC, Oklahoma State University, 1981.
CHAPTER 5
5.1 INTRODUCTION
that conventional model forms are not particularly suitable for use in the
system design situation,
that dynamic power flows are particularly descriptive of the physical actions
dictating the dynamic response of hydraulic control systems,
that power flow models offer the system designer modular reuseable model
forms which are close to direct representation of the similarly modular
hydraulic system,
that power flow modelling allows increased separation of the essential steps,
and consequent increased formality, in forming a set of dynamic equations
(the model) suitable for solution by computer.
The present Chapter introduces the power bond graph technique for developing
dynamic models of powered control systems. It has become customary to refer to
power bond graphs simply as bond graphs. Although the present author considers this
to be an undesirable abbreviation, it is adopted for the Text. Bond graphs have
49
A bond graph is a dual-signal flow diagram formed from a limited set of symbols.
From it a properly ordered and arranged set of dynamic equations can be prepared in a
form suitable for computer solution. A new but very limited set of terms and
symbols has to be learned. In effect, the multiport power flow approach of
Chapter 4 is made more general, more concise, and more formal through the (power)
bond graph symbols. We shall begin by introducing and discussing the basic terms
and symbols which w i l l b e used in subsequent development of bond graph derived models.
The reasons why bond graph techniques are being proposed in this Text as superior in
the design-orientated analysis area will be presented, discussed, and demonstrated
after the reader has gained an appreciation of what bond graphs are.
to allow ready formation of, and subsequent changes in, the structure as is
required in the creative design situation;
used as the general symbols for effort and flow respectively. In specific situations,
specific symbols will be used; for example p for effort (pressure) and Q for f l ~
5.2.2 Sources
One or more of the power variables associated with a dynamic system can be
constant or can be assumed to be constant for purposes of analysis. For example,
some hydraulic control systems are analysed on the basis that supply to the
control valve is at constant pressure,
In bond graph terms, each of these constants is called a 8ou~de, and is designated by
the symbol S - In general, in the present Text,
although specific effort and flow variables will be specifically described. For
example:
Designation of a power variable as a constant does not imply that the associated
power is constant. In fact, in a dynamic situation power will rarely be constant.
The product variable associated with the designated source can vary, thus allowing
power to vary. For the electric motor, voltage (effort) is constant b u t current
(flow) varies as dictated b y the load on the motor. For a constant supply pressure
hydraulic system, flowrate can vary and power flow will vary accordingly. For power
flow to be constant, either both effort and flow variables must be sources, or both
must be allowed to vary such that their product remains constant. The former case
is not compatible with dynamic states, though the latter may be sought in some
control system designs.
A pc~er state variable which is independent of the system, though not necessarily
constant, is also regarded as a source.
E
Q
P
Q
A s d i s c u s s e d in C h a p t e r 4, m a n y o f the m a j o r c o m p o n e n t s in a h y d r a u l i c control
system are p r i m a r i l y p o w e r transforming devices. Thus, an e l e c t r i c m o t o r c o n v e r t s
electric p o w e r to m e c h a n i c a l power; a h y d r a u l i c p u m p c o n v e r t s m e c h a n i c a l p o w e r to
hydraulic power; a h y d r a u l i c a c t u a t o r c o n v e r t s h y d r a u l i c p o w e r back to m e c h a n i c a l
power. P o w e r t r a n s f o r m a t i o n is r e p r e s e n t e d in b o n d g r a p h s as the symbol TF in the
centre of a p o w e r bond, as follows:
- - TF
and,
E1 E2
Q1 Q2
which m e a n s t h a t
E1 t r a n s f o r m s to E2
Q1 t r a n s f o r m s to Q2
E1 Q1 = E 2 Q2
P F
TF
Q A
P A = F
Note t h a t p o w e r is conserved,
(P - A) - Q/A = P - Q = F - X .
T P
Q
vp-i
which means
T V p -I = P
Vp = Q
T w = P - Q
Q = t0 - v p
T = P - Vp
T P
- - MTF - -
Q
V p -I ()
Losses and dynamic effects which occur during power transformations are accounted
for separately from the actual TF operation, which is treated as ideal; represent-
ation of losses and dynamic effects in bond graphs will be discussed and demonstrated
shortly.
Gyration
E1 E2
GY
Q1 K Q2
E1 K = Q2
QI/K = E2
rather than
E1 K = E2
QI/K = Q2
E T
GY
i
K
which means
E " K =
i/K = T
(E K) i/K = E i = T
Resistive Power D i s s i p a t i o n
C a p a c i t i v e Power Storage
Inertive (or Inductive) Power Storage
R C I
Ill
That is, a p o w e r b o n d l e a d i n g to one of the three b a s i c d y n a m i c a c t i o n classifications.
E a c h c l a s s i f i c a t i o n will n o w b e d e s c r i b e d through examples.
The examples all illustrate that power dissipation actions are resistive in
nature. Further, they cause irreversible power loss. The bond graph symbol
recognizing a power dissipation action is a power bond with the letter R (for
Eesistive effect), at its end. The R indicates that power transmission in the
region of the b o n d is dictated or affected by a resistive effect. Power flow in the
bond is described by the product o f an effort variable and a flow variable. Thus
the symbol becomes
Rh
dPh Qh
If, o n the other hand, flow in the line w a s turbulent, the Rh effect cannot be
accounted for via a simple coefficient. The pressure drop-flowrate relationship is
better described b y the non-linear equation
I r r e s p e c t i v e of w h a t Rh r e l a t i o n s h i p is e v e n t u a l l y chosen to relate Qh to
APh , the b o n d g r a p h symbol r e m a i n s as shown. The symbol r e c o g n i z e s the action, not
the s p e c i f i c e q u a t i o n to b e used. The symbol i n d i c a t e s that f l o w r a t e is a function
o f p r e s s u r e drop, or vice versa. E q u a t i o n s e l e c t i o n c o m e s later.
Rf
which means
Rf
Ff
where Ff is f r i c t i o n force r e s i s t i n g m o t i o n X .
Ff
s l o p e Rf Ff = Rf X 5.4
Ff Ff
mechanical springs,
drive shafts in torsion,
mechanical compliance in general,
hydraulic accumulators,
electric capacitors,
in general, oil-containing volumes in hydraulic control systems
(filters, hydraulic hoses, cylinder end volumes, etc.).
In general, these devices can store and give up potential energy in manners
affecting dynamic performance of the system of which they are part. The bond graph
symbol for a capacitive effect is:
t
the line (power bond) representing the path through which power can flow to or from
the particular capacitive effect recognized by the C at the end of the bond.
Implied on the bond are the effort and flow variables associated with the particular
C being considered,
E = C(Q) or Q = C(E)
Note that equation 5.7 is the dynamic form of the static spring relationship
Ff = 1/Ck X
relationships will be demonstrated later in the Text. Note that w~ have, in each
case, expressed the potential variable as a function of the flowrate variable. This
satisfies the b a s i c causality requirement that the effect is a result of the cause.
with the power bond (the line) indicating the means through which power may be given
to the inertia, or given up by it. Adding the ~ffort and the flow variable
associated with the possibilities
leads to
Im ij
Fm ~ or TjI~J
whe re Fm is the net force applied to mass Im to cause dynamic linear velocity
Xm of the mass;
The relationship implied on the I power bond is, of course, Newton's law.
However, it is best expressed as
= I(F) or ~j = I(Tj)
Fm = I(X) or Tj = I(~j)
In fact, when we realize that force (torque) is cause, and velocity (angular velocity)
is effect, it is more appropriate to express Newton's law as:
Xm = I(Fm) or mj = I(Tj)
or
f
~j = i/Ij )Tj dt + ~j(O) 5 .i0
E l e c t r i c a l i n d u c t a n c e is a n o t h e r d y n a m i c I effect. It is d e s c r i b e d b y a
d y n a m i c e q u a t i o n similar in form and nature to N e w t o n ' s law, i.e.
E = L di/dt , or i = 1/L IE dt
E Q2
I
E
Ql 0
Q3
Q1 + Q2 + Q3 = 0
E = same in all three b o n d s .
T
Q1 * P ~ Q3
61
P P
Q1 0 Q3
which means
Q1 + Q2 + Q3 = 0
p = same in all three bonds
E2 Q
E1 E3
1
Q Q
which means
E1 + E2 + E3 = 0
Q = same in all three bonds.
A typical use of the 1-junction occurs when we wish to allow for a frictional
resisting force opposing the motion of a hydraulic cylinder-driven load, as follows:
~--~-- Xm
Fa ~ I Im I
J Ff
7/////////
where Fa is pressure induced force acting on mass Im to induce displacement Xm ;
Ff is friction force.
Fa Fm
1
is i n t e n d e d to e x p r e s s
Fa - Ff = Fm ,
the net force a c t i n g o n the mass. It is a p p a r e n t enough that the p i s t o n rod, the
mass, and the friction force all e x p e r i e n c e the same v e l o c i t y , leading to the
c o m p l e t e d 1-junction,
Rf
Fa Fm
1
w h i c h means, algebraically,
F a + F f + Fm = 0
Xm = same for all b o n d s
Ff = Rf(~) or vice versa.
5.3 F O R M I N G P O W E R BOND G R A P H S T R U C T U R E S
SE
SE E f f o r t is constant, or independent.
Source
Q
E
SQ SQ Flow is c o n s t a n t , or i n d e p e n d e n t .
E1 E2 E1 t r a n s f o r m s to E 2 , Q1
TF - TF
Q1 Q2 t r a n s f o r m s to Q2 , such t h a t
Transformer EI-QI = E2.Q2 .
E1 E2 E1 transforms to Q2 , Q1
GY ~ GY-
Q1 Q2 t r a n s f o r m s to E2 , such t h a t
E1.QI = E2-Q2 .
Resistive Effect R R Recognizes a resistive effect.
Capacitive Effect C - - C Recognizes a capacitive effect.
Inertive E f f e c t Recognizes an i n e r t i v e or i n d u c t i v e
I I
effect.
E E
0-Junction
0 A l g e b r a i c s u m m i n g of flow v a r i a b l e s
a t c o n s t a n t e f f o r t , i.e.
Q1 + Q2 + Q3 = 0 , E = c o n s t a n t .
EIQ2
E1 E3
1-Junction Algebraic summing of effort variables
Q
a t c o n s t a n t flow, i.e.
E1 + E2 + E3 = 0 , Q = c o n s t a n t .
E21Q
Fa
(e) Fm
Ff Im
a c t s on load i n e r t i a
(mass) .
compl e te d.
Structure l
I
!
Fm Im
Xm
F.~ rr I F.m Im
X X
FIG. 5.2 DEVELOPMENT OF BOND GRAPH FOR INERTIA LOAD WITH FRICTION
84
Step (c): Fm acts to accelerate the load inertia; this is recognized b y placing
Im at the end of the Fm bond.
Step (d): the two power symbols appropriate to each bond are included on the
structure.
Formation of equations from the bond graph will be dealt with later in the Text.
Step (a) : it is recognized that the cylinder will receive power as the product of
pressure and flowrate of oil, P Q .
Step (b) : it is recognized that the cylinder will provide power as the product of
force Fa and velocity Xa . The implication is that the cylinder
w4ZZ be connected to a resisting load.
Step (c) : it is considered that the capacitive effect of the cylinder chamber
receiving flow is dynamically significant. A certain (small) portion
Qc of the inflowing oil is absorbed in compression of the oil in the
cylinder. Thus, the effect requires summation at an 0-junction and
is designated as a capacitive e f f e c t Ca .
step (e): the presence of a frictional force Ffa opposing motion of the ram,
and due to piston seals and rod glands, is recognized and allowed for
Step (f): the inertia la of the piston and rod will absorb some power as it is
accelerated. This is allowed for by the Ia bond; it leads from a
85
mPQ Xa
I "
J
Fa
~a
I i
for by Ia at a
1-junction. This
completes the structure.
0~~~1 ~1
Ca Rfa Ia
(g) Power variables added to
PIQc FfalXa F21Xa
each bond.
....... P ....
Q
0 ~ P
Qa
~
A
"-~1 "-1Fo
~a Xa ~a
Ca
(h) Simplification if cylinder
inertia Ia and friction PI Qc
Rfa effects are neglected.
P 0 P Fa
FIG. 5.3 D E V E L O P M E N T OF A B O N D G R A P H S T R U C T U R E
FOR A HYDRAULIC CYLINDER
1-junction as it is force (effort variable) F2 which is being allowed
for. The designated effects have all been accounted for. An output
bond is joined to the last 1-junction to indicate the actuator will be
joined to a load of some kind. ~his completes the structure.
Step (g) : Appropriate power product variables are added to each bond.
The bond graph structure developed for the cylinder is reasonably detailed. It
can be simplified by neglecting the cylinder rod inertia and friction (they can often
be lumped with the inertia a n d friction of the load attached to the rod). Such
simplification means neglect of the Rfa and Ia effects on Fig. 5.3g in which case
the bond graph structure becomes Fig. 5.3h
Fig. 5.4 shows the steps in development of a bond graph structure, starting with
a constant voltage (SE) electric power supply, and ending with a torgue (To) and
speed ( ~ o ) output power product. Points to note are:
Step (b): it is recognized that the motor gives out power as the mechanical
product torque (To) speed (co) .
Step (c)~ the gyration (GY) of constant voltage to constant motor field speed
(Source S~) is recognized.
Step (d): the speed slip essential to the action of an induction motor is
recognized by the P~ bond from an O-junction. It is a resistive
effect.
Step (e): torgue loss due to friction is allowed for by the resistive Rf bond
leading from a 1-junction.
Step (f): absorption of some torgue to accelerate the motor rotor and shaft
assembly is allowed for by the Ix bond leading from a 1-junction.
Step (g): the torsional compliance of the output shaft is allowed for with the
Ck bond leading from an 0-junction, recognizing that some speed
67
R~ Rf Ix
(f) Torque absorbed
accelerating motor-
shaft inertia Ix .
GY ~ 0
I { +
,, 1 ,a, 1
To
OJX
To
~00
R~ Rf IX Ck
(g) Speed variation due
to compliance of
motor-shaft,
This completes
Ck .
I { I +
structure. SE GY
~ 0 ~ I I To 0 To
l
Rf Ix Ck
(h) Structure with
power state
variable symbols
added.
Ti law Tf ] ~x
SE
l
GY Ti
S~
0 Ti
~x
I 1 ~To 0 -~
(i) Simplified
structure, if shaft Rm
compliance Ck ,
rotor inertia Ix , TO lAw
rotor friction Rf
are all neglected,
and motor field
~-----0
S~
To
~o
speed ~ is
regarded as constant
(i.e. as a source
S~ ).
FIG. 5.4 DEVELOPMENT OF B O N D G R A P H F O R I N D U C T I O N ELECTRIC MOTOR
RECOGNIZING S L I P (R~), F R I C T I O N (Rf), I N E R T I A (Ia), A N D
SHAFT COMPLIANCE (Ck).
68
Step (h): dual variables appropriate to each bond have been added. These are
the variables which will appear in the dynamic equations.
It can be seen that the bond graph comprises a main power flow sequence from
left to right consisting of power bonds, transformer, 0-junctions and 1-junctions.
From the main stream, power bonds lead from the junctions to the various R , C and
I effects judged to affect the motor dynamics.
Consider the hydraulic cylinder driven load combination shown on Fig. 5.5a.
The bond graph structures for each of the two major components, previously developed,
are shown in Fig. 5.5b and c, directly beneath the appropriate hardware components.
Actuator inertia and friction are neglected, knowing that they can be lumped with the
load's values in this rigidly connected situation.
The system bond graph structure is formed by simple joining of the component
bond graphs as shown in Fig. 5.5d. The Fa - ~ bond directly represents the
cylinder rod which connects cylinder piston to load. The ease with which the
component structures are joined to form the system structure is a characteristic of
the power bond graph technique.
89
a. The Sy s tern
Xm
Pa i >
Qs
W f ~ r - f r f
Friction
b. Cylinder Bond Graph (Inertia,
Friction, and Return Line c. Load Bond Graph
Effects Neglected)
, pa ~ pa TF Fa Fa I Fm
es ~ Qa A Xm ~m ~m
PalQc FflXm
ca Rf
Pa ! Qc FflXm
Ca Rf
e. System Bond Graph with Power Flow Direction and Causality Assigned
(Ref. Section 5.4.6)
~ Im
Ca Rf
More detailed bond graph structures of the linear actuator driven load will be
developed in Chapters 7 and 8. Later structures will allow for return line effects,
compliance in the actuator-load connection, internal cylinder leakage, and close-
centre, open-centre, and tandem-centre valve control.
Fig. 5.5e shows the bond graph with power flow directions (half arrow on each
bond) and causality (short transverse bar across one end of each bond) decided. The
meanings of these decisions will be presented next.
5.4.1 Introduction
So far in Chapter 5, we have described the basic power bond graph terms and
symbols, and have used them to form structures for some components or sub-systems of
the kinds included in typical hydraulic control systems. Later, we shall develop
more sub-system structures, and also structures for complete systems. A power bond
graph structure recognizes the presence of the various R , C and I effects which
dictate the dynamic response of the system. The structure does not indicate the
relationships by which the dynamic effects are described. That comes later.
However, the structure is useful only if it is leading towards eventual preparation
of a set of dynamic equations which describe the system dynamics. The assigning of
power flow directions, and of causality, are steps preparing the bond graph to yield
the appropriate set of equations. Assigr~ent of power flow directions provides the
signs of the terms in the algebraic s~um~ing equations to be eventually written for
each 0-junction and 1-junction. Assignment of causality reflects decisions on which
power state variable is to be cause, and which is to be effect, on each R , C and
I bond and its equation.
decide signs b y the physical action which takes place in the system
immediately after it is activated.
R C I
The directions of power flow in the remaining bonds of a bond graph follows physical
common sense - - p o w e r flows from the source to the load, usually from left to right
across the bond graph.
Consider the bond graph structure developed for the hydraulic cylinder in
Fig. 5.3, and shown again on Fig. 5.6. Obviously, it is intended that power flows
from left to right, i.e. from the hydraulic supply P Q towards the driven load.
Hence, directional half arrows are added on the m a i n (horizontal) power flow path
from left to right. For Ca , Rf and Ia , show the half arrows pointing into the
effect, utilizing the convention that power flows into these elements.
It is usually a simple matter to designate the power flow direction on each bond.
For the present, this m e r e l y helps in appreciating what is going on in the system.
Later w e shall use power flow directions in deciding the signs in the equations being
formed from the structure.
5.4.3 Causality
Ca Rf Ia
Fa L~
~a
Ca Rf Ia
i ~0
I 1 T
FIG. 5.7 BOND G R A P H OF H Y D R A U L I C CYLINDER (EX. FIG. 5.6),
WITH CAUSAL BARS ADDED.
73
variable as a cause can be extended and applied to a power bond graph structure
so as to show how the equations to b e prepared for R , C and I effects are to be
arranged. Basically, it is m o s t appropriate to write effect = f(cause) ; for
subsequent digital simulation, it is preferred that a differential relationship be
expressed in integral rather than derivative form.
A causal decision is indicated b y a short transverse bar, called the causaZ bar,
drawn across one end of a power bond. For example,
with the causal bar drawn at the R end of the power bond, indicates that the R
expression is to be written in the form
Q = R(E)
where R() indicates "a resistive function of (what is in parenthesis)". That is,
it has b e e n decided that E is the cause and Q is the effect for this particular
R element.
If the causal bar is placed on the power bond at the end away from R ,
E = R(Q)
Causal assignment should reflect reality both as regards to how the system
components behave and h o w particular equations should be arranged. For example,
computer solution (simulation) procedures prefer equations in integration forms
rather than in differential forms. Thus, capacitive equations should be expressed:
74
Q = C(E) = (possibly) C
(or Fnet = Im X m )
Im
means that force F is directed at the mass Ira and velocity X is a resulting
(back) effect (which is not to say that velocity is backwards, or negative:).
Conversely, for a hydraulic accumulator,
For example, there can be only one effort variable at an 0-junction, as previously
demonstrated, indicating that there can be only one causal bar at an 0-junction,
!
lo I
If, as in the above case, a C element occurs at an O-junction, the other bonds at
the junction ~ s t have causal bars away from the junction.
On the other hand, there can be only one flow variable at a 1-junction, as
previously demonstrated, and this means that there, can be only one bond at a 1-junction
without a causal bar at the junction. Representation of an I element demonstrates
this
--Ill
with the preferred causal bar against the I end of the I bond, the remaining two
power bonds nTust have causal bars at the 1-junction if the symbolic representation is
to make physical sense.
flow directions.
As a more complete example, consider the bond graph structure for the hydraulic
cylinder from Section 5.3.2 and Figs. 5.3 and 5.6. Fig. 5.7 shows the bond graph
with causality added, the required steps being:
Fig. 5.4 showed the development of a bond graph structure for an induction
electric motor. Fig. 5.8a shows it with power flow directions assigned. Fig. 5.8b
shows the bond graph with causality assigned. Considering Fig. 5.8a,
2. The causal bars of the other two bonds at the Ck 0-junction must now be
placed at the bond ends away from the 0-junction.
3. The causal bar of the left bond at the Ix 1-junction now must be at the
1-junction.
4. Because the right-hand bond at the Rf 1-junction cannot now have a causal
bar at the junction, the remaining two bonds must have causal bars at the
junction.
This leaves the gyrator bonds and the Rm 0-junction as causally undecided,
except that the right-hand bond from the 0-junction cannot have a causal bar
at its junction end. Rational thought lets us resolve this.
77
K~ Rf Ix Ck
SE
i
~ GY
'~0
] "-1
] ] ~1 "0
1 To
~o
b. Causality Decided
P~ Rf Ix Ck
sE ~,
i GY I
,
1171
0 "-', 1 ", 1 : ~ 0 :o .:
T = R(~)
by placing the Re bond's causal bar at its 0-junction end. This immediately
locates the causal bar for the adjacent left bond at its GY end, leaving only
the SE/i bond causality undecided.
The causal assigr.nent is now complete, and the bond graph is ready to have the
dynamic equation set prepared from it.
Consider the system of previous Fig. 5.5 whose bond graph structure was shown on
Fig. 5.5d. Using only the criteria
the power flow directions and causal bars of Fig. 5.5e are simple to complete.
5.4.7 Summary
The reasons and the techniques for assigning power flow directions and causality
have b e e n presented and discussed. Their significance will be confirmed in Section
5.5, where equations will be developed from power bond graph structures. However,
patterns which reflect b o t h computational requirements and physical reality have
emerged, and will be summarised here.
Power flow direction and causality decisions are intended to help in the
processing of a rational set of dynamic equations from a power bond graph
structure.
79
Given a bond graph structure, allot each C and I element its preferred
causality; c effects have the causal bar placed away from the C ; I
effects have the causal bar placed against the I .
a 1-junction can have only one bond without a causal bar adjacent to the
junction,
If all remaining causalities do not flow automatically from the above steps,
seek the most physically obvious causal decision(s) to allow completion of
causal assignments.
Each power bond must have a causal bar, at one end only.
ie or {
For GY power transformation the opposite holds, i.e. there must be either
80
SE E f f o r t is constant, or i n d e p e n d e n t
SE
Source
Q o f the system.
E F l o w is constant, or i n d e p e n d e n t of
sQ
SQ the system.
E1 E2 E1 t r a n s f o r m s to E2 , and Q1
TF -TF
Q1 Q2 t r a n s f o r m s t o Q2 such that
Transformer EI oQI = E 2 - Q 2 .
E1 GY E2 E1 t r a n s f o r m s to Q2 , and Q1
GY
Q1 Q2 t r a n s f o r m s to E2 such that
E I . Q I = E2-Q2 .
O-Junction
0 7 E0__ Q3 A l g e b r a i c s u m m i n @ o f flow variables
at same E , i.e. Q1 + Q2 + Q3 = 0 .
E [Q2
1-Junction E1 I E3 A l g e b r a i c s u m m i n g o f e f f o r t variables
- -Q -7 at same Q , i.e. E1 + E2 + E3 = 0 .
E2JQ
P o w e r F l o w Directioni s h o w s d i r e c t i o n o f p o w e r flow i n each
(half arrow bond.
o n one e n d
E a s i l Y d e c i d e d for m a i n p o w e r flow
o f each
se que nce.
bond)
P o w e r flows into R , C , I
e l e m e n t s (a convention).
G i v e s s i g n c o n v e n t i o n for equations
to b e formed for junctlons.
ie - - ~ I or I GY--
b u t not I GY - - I nor I GY l - -
Fig. 5.9 shows t a b u l a t e d the c o m m o n b o n d graph terms and symbols. Fig. 5.9 is
similar to Fig. 5.1, b u t w i t h p o w e r flow d i r e c t i o n a n d causality included.
5.5.1 Introduction
3. Fa = ? to be s p e c i f i e d if F a is an e x t e r n a l input,
or to b e d e s c r i b e d from b o n d g r a p h structure
o f the c o m p o n e n t d r i v i n g the load, as
d i s c u s s e d i n the n e x t example.
Xm
Im
Fa I/2////z
Ff
Rf
a. System
xm
Pa
m
Qs
r
Im
Friction
l Pa Pa ~ TF Im
0 Qa " A Fa ~ 1 1 ~ x m
xm
Ca Rf
5.5.3 An E l e m e n t a r y S y s t e m E x a m p l e
4. Fa = Pa A --- t r a n s f o r m a t i o n o f c y l i n d e r p r e s s u r e Pa to
force Fa via m o d u l u s A .
7. Qs = ? --- d e t e r m i n e d b y n e x t u p s t r e a m component.
It is informative and significant that the causal bar had to be placed at the
input end of the Pa Qs bond. It conveys that the component directing flow to
the actuator is essentially a f!owrate provider, and that pressure is a consequent
back e f f e c t dictated b y the actuator and load system. The bond graph reveals the
fallacy of thinking that pressure can be supplied to an actuator.
that the pumping system supplies oil to the control valve at constant pressure
and that pump supply is adequate. This allows neglect of the dynamic effects
of all components upstream of the control valve, reducing the hardware to be
considered to that shown on Fig. 5.12b. Further, supply pressure to the
valve is designated as a source, SP ,
that the flow in the return lines (cylinder to valve to tank) has a
negligible effect on load motion response Xm following valve actuation,
The assumptions allow reduction of the load and cylinder bond graphs to those, used
previously, shown on Fig. 5.12c under the respective components of Fig. 5.12b. We
have not previously required a b o n d graph for a 4-way valve.
When return line effects are neglected, only the supply port and the pressure-
side control port need be considered. Further, the ports are physically close, and
together with the internal passage connecting them, they can be lumped as a single
resistance. It is of course pressure that is "lost" during flow through a restriction.
Thus, for dynamic analysis purposes, the left-hand bond graph structure of Fig. 5.12c
summarises the contribution of the valve to the present system's dynamic action. What
we have is
Rv
APvIQs
SP Pa
1
Qs Qs
with causality for the other two bonds to be decided. To resolve this, consider the
physical reality of either bond. For the Rv bond, for example, it is usual to
express flowrate through the valve as a function of pressure drop
Qs = f(APv, Xv)
Thus, the causal bar on the Rv bond needs to be at the bond's Rv end. This means
in turn that the causal bar on the SP bond must be at its 1-junction end.
a. The Sy stem
r- -I
xm
iI I, I II I, t ii ///////"
APviQs
sp I Pa 0 Pa Fa Fa I Fm , Im
Qs Qs ~
Pa IQc
Ca
d. System Bond Graph
Rv Rf
Apv~Qs
Qs Qa A Xm Xm
Ca
System's Equations
There are ten d i f f e r e n t p o w e r variables o n Fig. 5.12d, and so we need ten p o w e r
state equations. We shall p r e p a r e them in the o r d e r
Sources
I effects
C effects
R e ffects
TF
Junctions (0 to 1 from l e f t to right), giving,
Source : SP = constant
TF: Fa
= Pa A [ a r r a n g e d to h a v e p r e v i o u s l y u n u s e d o b j e c t
Qa = Xm - A J variables
Xv = (input to be specified)
Xm = I~ dt (desired p r i m a r y r e s p o n s e variable).
J
5.5.5 Summary
The examples given can be worked back to include a pressure relief valve,
hydraulic pump, an electric motor and other components. However, this will be left
to future examples and case studies. The purpose at present is to illustrate how
the equation set is formed from a b o n d graph structure.
There is not an exclusige order of the equations. However, the bond graph
novice is advised to seek equations in the following sequence. Having obtained a
bond graph with power flow directions and causality assigned,
write the effort equation for each TF and GY , at least in functional form,
using as the object a variable not yet the object of an equation,
write the flow equation for each TF and GY , again using as the object a
previously unused variable,
write a summing equation for each 0 and 1-junction, using power flow
directions to dictate signs, and choosing as the object of each equation a
power variable not already the object of an equation,
check that there is an equation for each power state variable, and that no
such variable is without an equation,
write an equation for each of any system variables required for study which
are not among the power state variables,
With some experience gained, an analyst m a y prefer to prepare the equations from
a bond graph in a somewhat different order; typically to start at the first (or last')
summing junction. For the systems of Fig. 5.12d, this would result in
QS = R(APv) = K Xv sgn Pv Pv ~
APv = SP - Pa
Qc = Qs - Qa
and so on. With this approach, the choice of object for each summing equation is
less certain, and may require subsequent re-arrangement of the equations.
The power bond graph technique results in a set of relatively simple equations,
readily checked for errors, and readily programmed for computer solution. The same
forms of equation appear in example after example, a n d those relationships relevant
to hydraulic control systems quickly become ingrained and are easily recalled and
used.
advantage. Indeed, the structural rigidity required for transfer function and
state-variable model forms may be a positive hindrance. Suffice it to keep in mind
that models in conventional forms can be formed from bond graphs if there is a
compelling reason.
Bond graph methods well beyond w h a t appears in the present Text have been
developed and used. The present Text concentrates only on those aspects of bond
graph methods which the author has found to be useful and appropriate to the dynamic
analysis and performance prediction of industrial hydraulic control systems. Readers
wishing to delve further into bond graph techniques can do so via References 5.1 to
5.4. However, there are several aspects of b o n d graphs not discussed in detail in
the present Text which should be b r o u g h t to the attention of readers.
90
In some systems, one or more effects may be modulated during system action.
The modulation can be independent of the system, or can be made via automatic feed-
back. In either case, it is usual that the modulation action does not utilize
significant system power. It is usual that one system variable, rather than a power
pair, is utilized to effect any internal system feedback. For example if a variable
displacement pump was adjusted manually during operation, no system power is utilized.
If a solenoid-operated hydraulic valve is actuated b y a system pressure, no power is
drawn from the hydraulic system itself. In general, logic, control, and information
signals are sensed, processed, and conveyed w i t h o u t loss of power from the powered
control system either because they are v e r y low-powered signals, or because external
power is utilized.
There are several ways modulation signals, if they exist, can be highlighted on
bond graphs. For example, in Fig. 5.12 it is intended that the control valve be
manually operated such that its resistance Rv is Xv dependent. This can be
illustrated as
R v (Xv) Rv --4~-Xv
o a T
For example, consider a pilot-operated pressure relief valve where the relief
valve is situated at the pump but its o p e r a t i n g pressure is sensed well downstream
at the actuator. Fig. 5.13a illustrates this situation, and includes the pilot
actuation line as a b r o k e n line. Fig. 5.13b shows the corresponding part of a bond
graph representation of the system, and includes the modulating feedback effect as a
broken line w i t h a full arrowhead on it. Thus, the modulation bond is distinguished
from power bonds b y being a b r o k e n line of light construction, and b y having a full
arrowhead on it in the direction of signal flow. Note that the modulation bond,
referred to in bond graph terminology as an a~tiue or a~tivatsd bond, has only one
power variable associated with it, in this case pressure Pa . The implication is
that the other power variable is zero and therefore that no power flows in the active
bond.
a. The S y s t e m
@
I
I
u_ Pilot L i n e J
b. Bond Graph
Pilot Signal
J %
/ %
/
I
Rvl %% Cal
!
Rr
/
/
1 ",1
~--,- Vp
~,-.,~ 0 ~-.,. 0 ,--- 0 1 ~-~ ~
Pa2 ~ / /
1
Rp
T
Cp Ka
J
Rv2 Ca2
~.,, / \ / k I \ /
Pump, Relief Valve, Filter, Line 4-Way Valve Hydraulic Motor Load
I I1~ MTFJ
i i
Vp
5.6.3 F i e l d s and J u n c t i o n S t r u c t u r e s
5.6.4 Simplifications
Rf Rf
I b~
I
Ia Ia
I Q2 E2I
0 ~ or EZ ~ I E3___E!__L
Q1 E Q3 Q
5.7 CONCLUSION
Bond graph terms and symbols have been introduced and their uses demonstrated.
Only those aspects of bond graph techniques which have been found to be relevant to
the present Text have been discussed in detail. For further reading in bond graph
techniques and their uses in a wide variety of dynamic systems, Refs. 5.1 to 5.4
are recc~unended.
5.8 REFERENCES
5.2 Thoma, J., Introduction to Bond Graphs and their Applications, Pergamon,
1975.
5.3 Transactions ASME J. of Dyn. Syst. Meas. and Control, V.94, Series G,
N.3, 1972 (whole issue).
5.4 Journal of the Franklin Institute, V.308, N.3, 1979 (whole issue).
CHAPTER 6
6.1 INTRODUCTION
Power flow modelling procedures provide the system model as a set of low-order
algebraic and differential equations. A detailed model of a system of some substance
could be a set containing a substantial number of such equations. Digital simulation
procedures are most compatible with this situation. There is no n e e d - - i n fact it
is undesirable - - to run the low-order equation groups together to form more complex
groupings of variables and coefficients. In the system design situation, in which
the sensitivity of system performance to changes of some parameters is likely to be
investigated, the simple presentation of coefficients and variables resulting from
power flow modelling is positively an advantage.
9S
Thus, power flow modelling and digital simulation are highly compatible. The
existence of the digital simulation solution procedure has made it possible to realize
the advantages of the power flow modelling procedure. Digital simulation is the only
solution method considered to be relevant to the present Text.
The digital simulation procedures relevant to the present Text are generally
classified as Continuous System Simulation Languages ICSSL's). A continuous system
is regarded as one which can have a continuous time history associated with its
dynamic performance, rather than one intended or known to have discrete states in
time.
Fig. 6.1 illustrates the relationship of the CSSL's among the dynamic problem
solving languages con~nonly used in digital computing. To help orientation of the
reader, some of the specific languages associated with each grouping are included.
The high level problem solving languages such as Fortran can be utilized directly to
prepare a model for solution. However, the CSSL's (for continuous systems) and the
discrete system languages provide the bridge between the user and the computer's
basic compiling language (Fortran or alternative).
The CSSL's fall into t~o categories. Those in the block-orientated group are
intended to fairly directly rival analogue ccmputing, in which a block diagram type
of model is most appropriate. The expression-orientated CSSL's anticipate receiving
the model in the form of a set of equations. It is this group which is of primary
interest in the present Text, and only this group will be considered further.
DIGITAL
SIMULATION
FORTRAN GASP
ALGOL GERT
COBOL SIMULA
BASIC etc.
etc.
Continuous System
Simulation Procedures
(CSSL's)
Expression-Orientated Block-Orientated
CSSL's CSSLfs
MIMIC DYNAMO
CSMP MIDAS
CSSL PACTOLUS
ACSL etc.
etc.
The CSSL's date back to the mid-1950's. Ref. 6.1 contains a survey of their
98
MIMIC (Control Data Corporation) and DSL/90 (IBM) were perhaps the major links
drawing the early digital simulations towards the standard CSSL concept. Both were
widely and generally used. Control Data Corporation, influenced by the standardization
move, now has CSSL 3 [6.4]. Similarly, IBM's DSL/90 has been developed into CSMP 3
(Continuous System Modelling Program), [6.5]. Xerox Data Systems SL-I [6.6] is
particularly close to the standard CSSL. Mitchell and Gauthier [6.7] have produced
ACSL, pronounced "axle" and referring to Advanced Continuous Simulation Language, for
use with most computers. Most CSSL's have been produced for use with a particular
company's computers. Thus, ACSL represents a major step forward.
Apart from the general purpose CSSL's, there exist many special purpose orientated
digital simulation procedures. For example, HYDSIM [6.8] is especially designed for
hydraulic control system response analysis using power port models of components.
ENPORT [6.9] has been especially developed to accept a model in bond graph form. It
has the disadvantage (very serious in the present Text) of being able to handle only
linear expressions. Both HYDSIM and ENPORT have the significant feature of being
able to store "standard" models of components likely to be used and re-used in various
control systems. Some other CSSL's are orientated to incorporate response optimization
schemes.
6.4 CONCLUSION
Programming and running time can be costly if inexpertly used. Round-off error and
precise accuracy may be difficult to estimate. Numerical instability can be a problem
where different parts of a model have widely different response speeds. This
situation may be reflected by the presence in the model of widely different natural
frequencies. These kinds of model are referred to as stiff. Numerical stability
in the solution of stiff models is being better overcome with the development of more
flexible integration procedures which do not necessarily require formal explicitness
of all of the relationships of the model. ACSL edges into this category. It may
require some experience in particular problem areas before the best integration
procedure and step size specification can be made with confidence. A ~odern CSSL
may have as options several Runge-Kutta integration procedures plus several predictor-
corrector types (i.e. Milne). Some are more complex than others, requiring more
computer capability. Thus a simple self-starting Runge-Kutta 4th order procedure
should be preferred to the more complex Milne 5th order form (needs starter, complex
but very accurate) unless it is known to be inadequate to a particular set of
equations or relationship.
CSSL's are very compatible with power flow models. They are also very compatible
with the equation sets and relationships which are relevant to hydraulic control system
response, including response optimization procedures.
6.5 REFERENCES
6.1 Branden, D., Burwen, M., and Wachouski, A., Digital Simulation of Physical
Systems, Assoc. for Computing Machinery.
6.2 Strauss, J.C., The SCI Continuous System Simulation Language (CSSL),
Simulation, V.9, N.6, Dec. 1967, pp.281-303.
6.5 Continuous System Modelling Program 3 (CSMP 3) and Graphic Feature, IBM
Canada, Ontario, Canada, 1972.
6.6 SL-I Reference Manual, Xerox Data Systems, E1 Sequndo, California, 1970.
6.8 Smith, C.K., HYDSIM User's Manual, School of Engineering, Oklahoma State
University, 1973.
100
7.1 INTRODUCTION
The bond graph material of the Text to here has emphasized formation of the
structure of a system's dynamic model rather than the actual equations whose solution
will provide dyn~m~ic response. Equations have b e e n presented in some examples, b u t
with little discussion.
The present Chapter brings together physical relationships which are relevant to
the Text. Where ambiguity or choice exists, some discussion is offered. Selection
of values for those coefficients which are not straightforward is discussed, and some
values for coefficients are included.
7.2 COMPLIANCE
7.2.1 Introduction
Fig. 7.1 illustrates the idea that for a particular hydraulic system and its
driven load, the pressure development at the actuator following a control command at
the control valve might look like
Pressure
l i
f {
II Ii f~% /r-- 3. Correct Bulk Modulus
' t / ..~.
# i ; i .i., ,
; ~ li ,-
i /___..l
Time
Pl
,2 I i I Im
B = - V AP/AV 7.1(a)
If v o l u m e V of l i q u i d is s u b j e c t e d to ~ n c r e a s e d p r e s s u r e AP , V dew,Base8 by
A V , i.e. AP/AV is negative, and the n e g a t i v e sign o f e q u a t i o n 7.1(a) is u s e d to
make bulk modulus B a p o s i t i v e value.
E q u a t i o n 7.1(a) c a n be e x p r e s s e d
B = - V dP/dV 7.1(b)
giving
B = V P/Qc 7.1(d)
~ene
C = V/B 7.1(f)
E q u a t i o n 7.1(a) can be e x p r e s s e d
B
AP . . . . AV 7.1(g)
V
V = v(o) + Xa Aa 7.2
Note that for single rod actuators, Aa is one v a l u e for the h e a d side a n d a n o t h e r
(smaller) v a l u e for the rod side. The p r o b l e m o f v a r y i n g volume does not arise with
rotary drives, the volume of oil u n d e r c o m p r e s s i o n o n e i t h e r side o f the r o t a r y
actuator r e m a i n i n g v i r t u a l l y constant.
The bulk modulus of a particular liquid is affected by its temperature and its
pressure, though the effects are not significant in the working pressure and
temperature ranges associated normally with hydraulic control systems fluids.
Aeration of the fluid, prevalent in working hydraulic systems, is a much more
significant factor. A typical mineral hydraulic fluid in its quiescent state has a
bulk modulus of around 1750 MPa. Under the influence of turbulence and cavitation in
working hydraulic systems, bulk modulus might drop locally to 10% of this value, or
even to zero locally if cavitation is fully developed. The lack of drive stiffness
caused by low bulk modulus can seriously affect the response and performance of
hydraulic control systems. In its cycle through a particular system the oil's bulk
modulus can vary from a high value in the reservoir, to a lower value after pumping,
to a lower value again after expansion through the control valve, and so on. In a
well designed reservoir, the oil has a chance to regain or approach its quiescent
bulk modulus before again being drawn into the working system. Merritt [7.1]
discusses and provides typical values for bulk modulus at different parts of working
hydraulic systems. A detailed discussion of bulk modulus and the variations in its
definition is provided in Ref. 7.2.
the shaft between drive and pump is subjected to torsional deflection when
under torque,
Most hydraulic systems are designed such that the mechanical components are very
stiff and will experience minimal deflection under load. Therefore, mechanical
compliance
A~ 7.3
where Xl , ~2 are the instant velocities of the two sections of the compliance
being considered
is the velocity difference between X1 and X2
Ck is the compliance of the elastic mechanical component
Kk = i/Ck is the "stiffness" of the elastic mechanical component
F is the force developed in and transmitted b y the compliance
AX is the instant deflection between the two sections considered.
Note that equation 7.4 is in the form required for power flow modelling whereas
equation 7.4(a) is the conventional static form for the force-deflection relationship
for a linear spring.
Equations 7.3, 7.4, 7.4(a) have the equivalents for rotation motion
= ml - ~2 7.5
where el , ~2 are the instant angular velocities of the two sections of the
compliance being considered
is the difference between el and ~2
is the torque developed in and transmitted by the compliant
component
Ck is the compliance coefficient of the component
= 1/Ck is the "stiffness" coefficient of the component
Ck
! F ~ 0 F ~!
~2 '
ope K f T
/
X
j /
/q /
FIG. 7.4 REPRESENTING FRICTION
109
7.3 FRICTION
Friction in these various circumstances can range between the two limits:
i. V{SaOU8 free,on, where good film lubrication is implied. That is, the mating
mechanical parts are separated by a film of lubricant, so that only fluid shear is
present to resist motion. Under these circumstances
Ff ~ ~ 7.V(a)
Ff = Kf " X 7.7
2. CouZomb J'~ei~on, where the mating mechanical components are in direct contact,
with incomplete or no effective lubrication. In this situation static friction has
to be overcome before relative motion can start, and the static friction force F~
is described by
F~ = ~ Fn 7.8
It is well known that when coulomb friction is overcome and motion starts, the
friction force drops sharply. Fig. 7.4b illustrates coulomb friction.
Many of the components in a hydraulic control system are flooded with the
hydraulic fluid. One of the properties desired of a hydraulic fluid is that it be a
good lubricant. This is certainly realized if a mineral oil is used, as is the case
for most hydraulic systems. Thus, good lubrication can be expected for the internal
110
moving parts of pumps, motors, actuators, and valves. However, the presence of
generous quantities of oil does not absolutely guarantee that no direct contact of
parts takes place. This is particularly so if the parts have stationary periods,
such as is usual with linear actuators, spool valves, and hydraulic motors. A film
of lubricant developed and stabilized when parts are in motion can be squeezed out to
some degree when motion stops. It is preferable to allow for this effect when
considering friction of even well-lubricated parts, by using a friction force-velocity
characteristic of the type illustrated on Fig. 7.4c.
Ff = F~-sg~~f-~ 79
where sgn means "sign of", allowing that X can be positive or negative.
Frictional forces and torques always oppose motion and irreversibly dissipate
(waste) power. However, it should not be presumed that friction should be eliminated
or even minimized. Friction is often the only mechanism providing damping of load
response to sharp inputs. A minimum-friction actuator driving a minimum-friction
load may well prove too oscillatory, or even result in limit cycling which borders on
instability.
Numerical values for friction coefficients are included in some of the examples
of Chapter 8.
structural compliance, and possibly of external load forces, is the task faced by the
majority of hydraulic control systems. The designer-analyst needs to understand the
basic load arrangements he m a y have to choose between or adapt frc~. For example to
control the rudder of an aircraft, the control system needs to overcome inertia,
friction, structural compliance, and the aerodynamic forces on the rudder. The
hydraulic system of a large press has to cope with inertia of the press head assembly,
friction in the press guides and bearing, deflections in the press frame, and the
large load induced when the press head comes in contact with the work piece.
Fig. 7.5 illustrates eight basic load arrangement models, together with their
bond graphs, causal requirements, and probable (not all exclusive) equations. The
first three entries are straightforward. The force Fa developed in the actuator
is the causal input variable, load-actuator velocity is the resulting effect. There
is no compliance effect between actuator and load, and this is often a reasonable
assumption.
Entry 6 illustrates a compliant drive between actuator and load, with damping or
friction also present in the connection. This might be the situation if actuator
friction was considered to be significant or even necessary to damp down oscillations
in Xa and Xm due to the presence of compliance Ck . Again, and for the same
reasons as for Entry 4, actuator force Fa is a back effect resulting from causal
velocity Xa . The 0-junction on the bond graphs tells us that the compliant drive
causes Xm and Xa to be different by the amount AX (and consequently Xm and
Xa to be different by ~X ). The 1-junction allows that both Ck and Rf affect
, and that they act in parallel with each other.
Entry 8 shows the situation of Entry 3 operating in the vertical plane, when
gravitational force acting on the mass must be included. This is shown as the
constant force SFg (Source, Force, gravitationai). SFg could have been added to
the single 1-junction of Entry 3.
Most hydraulic control system loads include inertias. Many also include other
forces. For example, a machine tool drive or tool positioner will experience cutting
forces. The tail-plane control on an aircraft will experience an aerodynamic force
which is dependent on both the speed of the airflow over the surface and on the angle
of attack of the surface to that airflow. A sensitive model for an eleotrohydraulic
servovalve may have to account for both the inertia of the valve spool and the flow-
forces generated on it.
Consider the case of an aircraft's rudder control. Both the inertia of the
rudder, and the aerodynamic forces on it should be considered when assessing its
position response to input cen~ands. Whatever else may have been included in a
bond graph for the control system and its load, the structure should end with the
inclusion of an R8 , vel bond from a 1-junction and with causal bar at the junction,
indicating recognition of a force which is dependent on rudder angle (8) and on
airspeed (vel) . It is of no consequence whether Im is shown projecting from the
same 1-junction. It can project horizontally as previously, or at any other angle
for that matter.
7.5.1 Relationships
The working parts of pumps, valves, and actuators must have clearances. For
example, the pistons in their bores in piston pumps and actuators, the spool in its
bore in spool valves, hhe gears and the side plates in gear pumps, etc. These
clearances, though very small, will act as leakage paths when under high pressure
drops. Normal hydraulic fluids are sufficiently viscous, and working clearances are
sufficiently small, that leakage flows can be assumed to be laminar. That is
QZ ~ dP 7.10(a)
Then,
Q = I/R ~P 7.10
113
EQU&TIONS C ~ v ~ HTS
1. Mall only
ga - causal f~rce Fa derlwas f z ~ dxlwi~
ece~orent.
Fa ILl
9= - iln~.[rad~+~(o) Wew~n'm law.
"J--CLt] xm - I~m d t + m . ( o ) Auxiliary equation.
l/~-I~ ~t+~(o)
R(| d l n o ~ l a " a r l l i l t i v e
function of'.
- K~,~m for visc.frlc K is visc. frlc. colff.
Fm - Fa-FZ Net fozol on I~ .
C~ Fa kl Fm ~! Im Ff " RlXm)
Fa ~ ~ ~= ,1 9= rk - c19=) CO d e n o t e s "a ~ a t i t i ~
f u n c t i o n of'.
For l i n e a r a p r i n g , o f
-I c~ff. ~k .
ok Fm = Fa-Ff-Fk
1 ~&
~ -: m
~
ra
F~
-
-
-
~.-~
c(~)
Fa-rf
Ff RiAX) Use d e s i r e d e q u a t i o n foz
friction.
6. Mass, w i t h compliance and f~ictlon ~a - causal input AS for 4 and 5,
i n th e d r i ~ Fa ~ Fa = : I m i.e. for ~ p l l ~ t drives
xm = i/Zm-[Fa dt+Xm(o) 4, 5, 6, ~elocity is the
cause varlable and force
the resulting eZfect
variable.
Ff - R(AX)
Interesting, isn't it?
rk = c(al)
Rf ~ ~,,I~ ~ ~
Pa - Ff+Fk
~. Ks f o r 6, p l u s external f o r ~ Fi
xeslst/ng motion
~a - causal input
- c(~x)
Ill
S. ~ f o r 3, b u t s u s p e n d e d v e r t i c a l l y
~f
Fa
~m sFq i s c o n s t a n t f o r c e
As f o r l . e x e r t e d o n ma~s due t o
rf !
Fa I gravity.
~he ~ o a d j a c e n t l - J u ~ t i ~ s
Fro' - Fa-F f - F k can be combined.
Pm - nn'-SFg
Fa (~k SFg
AP will be a variable in the power state equations and the simulation. Usu~ID
leakage flows back to tank, the discharge side pressure can be taken as zero, and
approximates the pressure at the entry to the clearance. Thus, leakage from a p~mp
can be expressed
Qp = i/Rp Pp 7.10(h)
Q = ~ D H3/(96 ~ L)J"
~ Ap 7.11
However, the attitude of a piston working in its bore can rarely be guaranteed.
The piston can tilt, and can become eccentric in the bore as illustrated on Figs.
7.6b, c, respectively. For the case of full eccentricity, Fig. 7.6d, it is readily
shown that QZ will be 2.5 times the value calculated from equation 7.11. That is,
for full eccentricity, R drops to 0.4 of its value for the concentric case of
Fig. 7.6a. It is worth observing that leakage is proportional to clearance cubed.
That is, if clearance is doubled, the resistance is decreased to an eighth, and the
leakage flowrate goes up eight times. It is also w o r t h noting that clearance is
affected - - increased ....b y pressure expansion effects, increasing leakage from the
value calculated using manufactured clearance. Leakage in the clearances of piston
components is discussed more fully in Ref. 7.3.
7.5.2 Coefficients
a. Concentric
2. i ~," / I I ! I :
T
b. Tilted Piston
J / / J / J J / J J
/ I / / / I I'////
C. Eccentric Piston
Ll / L / / i / / / /
/ / / / / / / / / / /
d. Fully Eccentric
H q _/ / / / / / / / / / /
I
/ J / / / / ," / / / / /
7.6 R E L I E F VALVE F L O W R A T E S
7.6 .I Relationships
Qr = f(p) 7.13(a)
Qr = 0 if P .< P s e t
7.13
= i/Rr (P - Pset) if P > Pset
where i/Rr is the slope ~Q/AP of the Qr(P) line for P > Pset .
117
slope 1/Rr
P
Q1
~
PIQr
~%
U
P
Q2
Qr = f (P)
Q2 = Q1 - Qr
! ~-P
Pset
FIG. 7.7 PRES SURE-FLOWRATE FIG. 7.8 REPRESENTING A
CHARACTERISTIC OF A PRESSURE RELIEF VALVE
PRESSURE RELIEF VALVE BY BOND GRAPH
P~ Rf Ix Ck
T~
i| OperatingRegion
!~ Shaft Speed, ~x
Speed at Max. Operating Torque ~ I
Synchronous Fie'Id
Speed, S~
analyst can expect that a relief valve manufacturer can supply it. The P > Pset
characteristic is not always a perfect straight line. It may have a curve to it, in
which case the curve can be approximated by a straight line, described b y a non-line~
relationship, or represented by data points. Alternatively, experimental data for
Qr(P) can be used in a simulation. Ref. 7.5 describes in more detail the case for
representing relief valves by equation 7.13 in system simulations. It includes the
establishing of Rr values for a number of different valves, both direct-acting and
pilot-operated. It describes also the results of simulations using the equation 7.13
representation of relief valves, and their support by experiment.
In the occasional event that an analyst may wish to study the dynamics of the
relief valve, for example to study chatter, a number of references are available
(Ref. 7.6 for example).
In stm~ary, for general system analysis a pressure relief valve can be represen~d
as an R-effect leading from an 0-junction in a bond graph, as illustrated on Fig. 7.~
with the Rr effect to be described by equation 7.13. Note that the single causal
bar allowed at the 0-junction can be o n either of the main power flow bonds, d e p e n ~
on the particular system the relief valve is used in.
7.6.2 Coefficients
As an example of Rx values, the system of Ref. 7.5 used six different relief
valves, all rated in the region 4 to 20 MPa set pressure, and 60 to 180 /min. Values
for Rr ranged from approximately 0.5 to 5 x 109 pa/(m3/s) depending upon the
particular valve and the particular set pressure. Values for Rr should be available
from the relief valve manufacturer, or should be assessable from his valve's perform~ce
data.
7.7 .I Relationships
Many hydraulic pumps are driven by induction electric motors. Often the motor's
output speed c a n be regarded as c o n s t a n t , in w h i c h case the motor's internal
characteristics and dynamics can be ~eglected and the input speed to the pump can be
taken as a source (constant) . However, b y its very nature, an induction motor's
output speed varies with load. Sometimes this "speed slip" will b e significant to a
particular hydraulic system response analysis, and it is then necessary to have a
dynamic model of the induction motor. In Chapter 5, Section 5.3.3 and Fig. 5.4, a
detailed bond graph was developed for this situation. Fig. 7.9 is a reproduction of
Fig. 5.4b w i t h the gyration from source voltage to source motor field speed omitted
119
friction, as the e f f e c t Rf
T~ = R~ /~0 7.14
Tf = Kf ~x 3 7.15
Kf is a f r i c t i o n coefficient
~x is m o t o r shaft speed
where ~x is m o t o r speed
7.7.2 Model
Source : Sm = constant
R~: T~ --- R e ~
Rf: Tf = Kf mx 3
1-junction: Ty = T~ - Tf
1-junction: Tx = Ty- To
0-junction: d~k = ~x - ~o
7.7.3 Coefficients
Sm = 1500 r p m
Rm = 53 N m/(rad/s)
Kf = 0.316 - 10 -3 N - m / ( r a d / s ) 2
Ix = 16.5 10 -3 k g m2
Ck = i0 10 -6 r a d / ( N m)
121
7.7.4 Conclusion
The electric induction motor has been modelled in some detail. The model
developed and the particular coefficients for it will be utilized in an example in
Chapter 8 where motor dynamics were considered to be significant to the analysis.
It is emphasized, again, that it is rare that electric motor dynamics are of primary
significance in the dynamic performance analysis of hydraulic control systems.
7.8 .i Relationships
For a particular drive speed ep and swept volume rate Vp the bond graph can
be simplified to that of Fig. 7.11b.
It is of interest and significance that the causal bar on the pump's output bond
is at the left, indicating that pressure is a back effect, to be induced by whatever
impedance the pump's discharge flowrate experiences. Thus, the bond graph indicates
unambiguously that a positive displacement pump is primarily a generator of flowrate,
and not of pressure.
Sup Vp
c As for a, but allowing for pump d. Allowing for the
friction (Rfp) and pump rotating capacitance Effect After
parts inertia (Ip) Pump Discharge
Rfp Ip Rp cp
up mp Vp Pp
Ip cp
~ Qp falls dramatically if
Pp 9 Ppset
J| ~pp
Set Pressure
For a variable displacement pump, the displacement per radian of the pump is
modulated by a suitable mechanism. For a swashplate pump for example Vp on Fig.
7.11a and c would have to be replaced by
vp = Vp' - e
where Vp' is the volumetric displacement of the pump per radian of pump rotation
per radian of swashplate angle measured from zero (swashplate vertical,
piston stroke zero)
is the swashplate angle measured from zero.
If allowing for variable displacement, the inertia of the swashplate and its
driving mechanism may have to be considered, whether or not an assumption of constant
pump rotational speed is assumed.
A pump invariably discharges into some kind of a volume. The volume will
include the pump's own porting and outlet, and can include the volumes of connecting
hydraulic lines, of a high-pressure filter, and an accumulator. The capacitance
effect of the pump's discharge side volumes can be recognized as a C effect, as
illustrated on Fig. 7.11d, where
The majority of hydraulic systems include a pressure relief valve adjacent to the
pump, and this can be allowed for as an Rr effect as shown in Fig. 7.11e, where
I n c l u d i n g i n t e r n a l l e a k a g e t a p r e s s u r e - c o m p e n s a t e p u m p has a Qp - P p relation-
s h i p o f the form i d e a l i z e d on Fig. 7.12. A t c o n s t a n t p u m p speed, Q~p is insensitive
to pressure, in the m a n n e r o f a n o r m a l p o s i t i v e d i s p l a c e m e n t pump, u n t i l a set pressurl
Pset is reached. P r e s s u r e in e x c e s s of Pset is u s e d to activate a m e c h a n i s m which
c a u s e s p u m p d i s p l a c e m e n t and c o n s e q u e n t l y Qp to d e c r e a s e v e r y sharply. When Pp
falls below Pset , the p o s i t i v e d i s p l a c e m e n t p u m p c h a r a c t e r i s t i c is resumed. The
r e l a t i o n s h i p of Fig. 7.12 is r e a d i l y d e s c r i b e d as
P u m p f r i c t i o n a r i s e s from
shaft bearings
m o v e m e n t o f control m e c h a n i s m s i n the p u m p
P r e c i s e c o m p o n e n t - b y - c o m p o n e n t a s s e s s m e n t of friction i n p u m p s is a c o m p l i c a t e d
125
process. However, it can be appreciated that the torque required to overcome friction
is affected by both pump speed, and pump pressure, lhe moving parts will be i~nersed
in oil, providing good lubrication. Pump speed is likely to induce a viscous friction
effect, i.e. the higher the speed the higher the frictional torque, pump pressure
will increase the normal loading between mated mechanical parts in relative motion,
and thus increase the effort needed to overcome friction. That is, in general,
Normally, a pump's speed will not vary much, and thus the ~p component of Tfp will
be fairly constant. Pump pressure can vary widely and quickly as a system operates.
Thus, the pp component of Tfp can be dynamic. All things considered, if pump
friction needs to be included in a particular dynamic analysis, it can be represented
7.8.2 Coefficients
7.9 .i Introduction
Because of
The 4-way valve has four ports through which power can enter or leave the valve.
One port accepts power as an input from the pump system and is designated the supply
port. Another port is connected to exhaust (tank). The remaining two, the control
ports, typically but not exclusively are connected to the two hydraulic ports of a
hydraulic actuator. The direction in which a driven actuator moves is dictated by
126
the direction in which the control valve is actuated. The speed at which a driven
actuator moves is dictated b y the amount b y which the control valve is opened.
closed centre: all four ports are closed and are isolated from each other,
Q = f(Am, Ap n) 7.20(a)
A = K ' X
127
Y Z
a. Closed Centre
I
S
I
E
Y Z S Supply Port
IIII i
E Exhaust Port
b. Open Centre
Y" Control Port
I I Z Control Port
S E
Y Z
i I
C. Tandem Centre
S E
Spool Lan
J-
FIG. 7.14 ILLUSTRATING SQUARE-PORTED SPOOL VALVE
For line-on-line, W = Wx
For overlap, Wx > W and overlap = 0.5(Wx-W)
For underlap, Wx < W and underlap = 0.5(W-Wx)
(Assuming axial symmetry)
128
m = 1
giying
Q = f(K X AP ~) 7.20(b)
Q = Kv X A~ % 7.20
s u c h that
Q = I / R ~P 7.22
of Fig. 7.13.
129
Allowing that the spool can be moved either left or right of centre, there are
four flow controlling restrictions available in the valve, as illustrated on Fig.
7.15c. If the spool is central, each of the four possible flow paths have infinite
resistance, and no flow takes place. If the spool is moved left of centre, Xv
positive Fig. 7.15d, the resistance of openings Rsy and Rze becomes finite,
allowing flow, while the resistances of openings Rye , Rsz remain infinite,
preventing flow. Thus, flow can take place from supply through control port Y ,
and through control port Z to exhaust. Had the spool been moved right of centre
130
( Xv negative), Rsz and Rye would become finite, allowing flow from s u p p l y to
port Z and from port Y to e x h a u s t .
Fig. 7.15e illustrates the flow paths of the 4-way closed-centre valve in a
Effect Equation
Rsy; Qsy = f(APsy): Qsy = Kv Xv APsy 7.23(a)
The junction equations, (a) to (), are general. However, the ~ ~<~uations,
(a) to (d), are in p r a c t i c e discontinuous. Each needs to b e e x p a n d e d to e m b r a c e
physical reality, as f o l l o w s
Qsy = 0 if xv is 0
Qsz = 0 if Xv is 0
= 0 if Xv is positive
7.23(n)
= Kv IXvl APsz ~ if Xv is negative a n d ~ 97
S ~ ~ Y
~ Xv
E Z
,
Flow Qsy _ ~
through Rsy through Rze
e. Pressure-Flow- f. General Bond Graph
Resistance Circuit
QYel Apye ~,.Rye
Qye
PS
~sy P
s Be
Qs
'Rsz
1
Qs~l ~ R~
f Qz
Rze
pe~ Ps
E Qe " ~ze S Qs LlOPs
g, Bond Graph for Xv
Positive Only
Rsy h. Overlap E~L~ 0p e Qszl APszNR~z
~Psy~ I wx
pz01 Pz ~ Z
__--=--% Qze Qz
Q~e = 0 if xv is 0 |
= 0 if xv is positive
7.23(0)
= Kv IXvl APye ~ if Xv is negative and ~ W
= Kv W APye % if Xv is n e g a t i v e and > W
Qze = 0 if Xv is 0
= K v Xv APze % if Xv is positive and ~ w
7.23 (p)
= Kv W APze if xv is p o s i t i v e and > w
= 0 if Xv is negative
a d e q u a t e to the p a r t i c u l a r task".
plus the R effect equations appropriate to open-centre with uniform underlap x(o)
and port width W
Equations 7.25(a) to (j) provide a model allowing for all degrees Df valve
opening, positive and negative and with port opening saturation. A particular
analysis might not require comprehensive description of the valve characteristics,
and the bond graph and the model could be simplified accordingly, as indicated with
the closed-centre valve.
s
/ / /////
E Z
l
QYel Apye ~, Rye
c. Fl~;-Resistance Circuit
(same as for C l o s e d -
C e n t r e Valve)
PY 0 1 PY ~ y
QY
Qye PY
~e
1
Qy Y Q~Yl ap~v .:~y
Rs
~Qs,
s Ps
Qs J ": 0 P~
Pz Z
Qz
E~OPe ~szI ~Psz~,R~z
E Pe ~ -~ Rze~
]-Qz
Qe Qze
Pz ~ z
Pz : Qz
Qze
Qzel APze L: Rze
Rye , Rze , of the closed and open-centre valves, plus a fifth o r i f i c e Rse through
which flow can p a s s from supply to exhaust.
APse = Ps - Pe 7.26(i)
137
xv
E
d. Bond Graph
Qye I APye~ ~
~ Rye
Osy I ~psy ~,Rsy
dRsy
S Ps~ s Ps 3 ~ N P~
Qs , Rsz
--~----~ Z Pe I~ ! T ~ A
~'Rze Qz
InPse T
E Pe Rse ~sz
Qe ~
Pz01 PZ ~ Z
Qz
Qze ~Qze
Qsz = 0 if Xv is 0
= 0 if Xv is positive
7 2 6 kl
= Kv-lXvl.APsz % if Xv is negative but < W
= Kv-W-~Psz ~ if XV is negative but ~ W
Qye = 0 if Xv is 0 I
= 0 if Xv is positive
= ~v Ixvl'~Ye if Xv is negative but < W 7.26()
= Kv W APye if Xv is negative but ~ W
Qze = 0 if Xv is 0
= Kv. Xv.APze ~ if Xv is positive but < W
7.26 (m)
= Kv.W.APze ~ if Xv is positive but 9 W
= 0 if Xv is negative
7.9.7 Summary
Consideration of the bond graphs of Figs. 7.15, 7.16, and 7.17 reveals that
Fig. 7.17d could be regarded as a general bond graph for all three centre configurations
considered. For the closed and open-centre arrangements, it would be necessary only
to specify that the resistive effect Rse is infinitely high for all spool displacement
Xv , allowing no flow Qse . The equation set derived from the general bond graph
would be different for each of the three cases only in the Q = f(AP) equations.
Fig. 7.17d can be put into the circular form of Fig. 7.18a, or in the form of
Fig. 7.18b. qhis latter form might be the most handy, and will be used in future
examples.
Equations 7.20 to 7.26 show the types of equations and equation sets required to
describe the valve flows. Specification of valve opening X for each flow equation,
in terms of Xv (magnitude and sign), and of particular valve port and lap conditions,
is the major challenge. All flows, if finite, are unidirectional except for Qy and
Qz , the external flows at the control ports. If Qy is positive, Qz must be
negative, and vice versa, and the relevant equations must reflect this.
7 .i0 ACTUATORS
7. I0.1 Introduction
Bond graph structures for a linear actuator were developed in some detail in
Chapter 5, Section 5.3.2. Fig. 5.3 showed the developments for the conditions
Inertia of the actuator's moving assembly (la) and its friction (Rfa) were
included in Fig. 5.3g, and neglected in Fig. 5.3h. In both cases the capacitance
effect of oil flowing into the cylinder is included. A linear actuator is included
in the system whose bond graph and equations were developed in Section 5.5.4, with
the same limitations associated with Fig. 5.3.
In the present Section, bond graphs for both linear and rotary actuators will be
further developed.
The back pressure developed in the discharge size of an actuator as its piston
is driven in one direction generates a force which opposes motion. This return-side
t40
a. A Circular Form
PY ~ y
Qy
Ps
Qs '~.0 ..... ":I I Q~e " O: Pe
Qe
~ E
Q~I ~ I/Qze
l Pz ~ Z
Qz
b. A Parallel Form
.................. PY ~ y
Rsy Rye
APsy~ Apy
IF---~O ~,II
APsz~ I APze~
Rsz Rze
............... Pz a,.- Z
Qz
effect can be included in a bond graph via the 1-junction shown on Fig. 7.19a, which
sums algebraically the forward-side pressure-generated force Fay and the return-side
pressure force Faz . The 1-junction merely describes the net pressure force Fa as
Note that the back pressure effect causes a power loss through a loss of effective
force.
Fig. 7.19b shows the development of Fig. 5.3g to allow for the back pressure
effect. Pressure Paz is generated because flow is being forced from the discharge
side of the actuator into the return line by forward motion of the actuator piston.
This is clearly a capacitive effect (flow into a C causes effort build-up in the C )
and this is recognized by including Caz in the return line bond graph. It is also
the reason why back pressure Paz is causal at the TF end of its bond.
If it was decided that actuator friction (Rfa) and inertia (Ia) could be
neglected in a particular analysis, or combined with the friction and inertia of a
driven load, they and their 1-junctions are simply deleted from Fig. 7.19b.
To allow for the effects of internal leakage across the actuator piston, the bond
graph can be modified to that illustrated on Fig. 7.19c. A bond from the 0-junction
of Cay , or from a separate 0-junction downstream of it, allows that some of the
incoming flow Qay is lost as leakage Qa . A bond into the 0-junction of Caz ,
or into a separate O-junction upstream of it, allows Qa to he added to Qaz . The
l-junction placed in the Qa flow bonds recognizes that it is pressure drop
(Pay - Paz) across the piston that induces leakage.
Fig. 7.19c is a quite detailed model structure for a linear actuator, and is
suitable for study of the dynamic performance of the actuator itself. It is suitable
als0 for connection to the bond graphs of inertia loads having compliance between
actuator and inertia such as were discussed with Fig. 7.5, Entries 4 to 7. Entry
and/or exit pressure loss can be included by adding 1-junctions at the left.
It is important to be aware that Fig. 7.19c and its causalities are not unique.
The actuator bond graph can depend also on the type of load it is to drive. For
example if it is to drive a simple inertia without compliance in the connection, the
causality of the force output bond on Fig. 7.19c will be incorrect. The inertia load,
Fig. 7.5 Entry i, requires the causal bar at the right end of the Fa bond of
Fig. 7.19c. A little thought resolves this conflict. If there is no compliance
between Ia and Im they should be lumped, reflecting the rigidity of connection.
The Ia bond of Fig. 7.19c becomes the actuator's output bond and the former Fa bond
is deleted as irrelevant, Fig. 7.19d. The same arrangement applies for the same
142
reason to the load systems represented on Fig. 7.5, Entries 2, 3, and 8. Actuator
friction and load friction also can be lumped together if there is no compliance
between them, i.e. if their bonds share the same velocity.
Note that the effects Cay and Caz cannot be represented by simple coefficient,
except for small perturbation response analysis. In general the cylinder end volumes
vary with their displacement as discussed in Section 7.2, and the formulation of the
Ca effects should reflect this.
The hydraulic motor can be viewed as the rotational equivalent of the linear
actuator. Its moving parts have inertia and friction, its chambers capacitance,
and it has internal leakage. Fig. 7.19 can be quickly modified to its rotational
equivalent by replacing force by torque, velocity by angular velocity, and the TF
modulator coefficients A by Va , where Va is the swept volume per radian of
rotation of the motor. The TF equations become
Pay Ha = T 7.28
Qay Va -I = ~ 7~29
The equation forms associated with Pig. 7.19 and its rotational equivalent have
all been discussed and will not be re-presented here. Attention is drawn to the
variation of the end volumes as the piston of a linear actuator moves. This affects
the capacitances of the linear actuator, as discussed in Section 7.2. The varying
capacitance effect does not occur to anywhere near the same extent with hydraulic
motors and can usually be ignored.
Loaded hydraulic actuators usually have low damping, with damping ratios of the
order of 0.01 to 0.2 depending on the conditions of bearings, seals, and glands.
143
pqy Fa, Xa
ay ~! / pa;I mv' Supply I Pay.:
Fa
Paz ~ TF --~- Fay
Qay Ay Xa~
"1 ---~I---~I ~ ~a
Fa Return ~ ~J
1 ~a ~-- 0 Qa~Az
. / Ia
A Paz TF ~ a z
Qaz Az Caz
AQay
~Qa
Pe
"~AQaz Az Ia
Caz
I ": (~m + ~
Az
~I:
C. The I n c o m p a t i b i l i t y of S h o w i n g
Adjacent C Effects. (Where
b. Including Hydraulic Line can the c a u s a l b a r go on the
Capacitance (Ch) c e n t r e bond?)
Cl c2
~Qh
"0 ~, b a d
Rh Ra
Rb Cb
Pl P2 P2
Accurate dynamic modelling of long fluid lines is a complex subject beyond the
scope of the present Text. Ref. 7.7 provides a useful and relevant review.
7.11.2 Filters
Apb = f(Ql)
Pl = P2 + APb
P2 = f (AQb)
AQb = Q1 - Q2 7.30
7.11.3 Accumulators
to provide supply flowrate, and hence power additional to the pumps, for
peak load demands o f short duration,
to run continuously,
where C = V/B
could be used. It should be borne in m i n d that V , the volume of oil under pressure
in and around the accumulator, will vary depending on the state of charge of the
accumulator. It should be borne in m i n d also that for most system analyses, the
primary requirement will be to study the dynamics of the driven load, and not those
of the accumulator. Usually a reasonable approximation of the accumulator dynamics
is adequate. If a more accurate model of the accumulator dynamics is required, V
may have to be appropriately modulated, and R may have to be described by a square
root of pressure drop relationship.
Pressure relief valves and directional control valves have been discussed in
Sections 7.6 and 7.9 respectively. It is not practicable to present here the bond
graphs of the large number of valve variations available commercially. Several will
be presented now, and it is left to the reader to form bond graphs for other valves
when they are required.
Check Valve
The check valve, or non-return valve, is intended to allow free flow in one
direction and to prevent totally flow in the other direction. Most check valves
Will have a flow-pressure drop characteristic of the kind idealized on Fig. 7.23a.
Negative pressure drop results in no flow. A small positive pressure drop is
usually present when flow is positive and this pressure drop can increase somewhat
with flowrate. Pressure drop in the presence of flow means that there is power
loss. Just how this is best accounted for on a bond graph depends on the location
of the check valve relative to its adjacent components.
For example, consider the system of Fig. 7.23b. There will be three pressures
involved when the load is being driven, these being supply Ps , pressure Px
between control valve and check valve, and pressure Pa downstream of the check
valve. To accommodate this the bond graph of Fig. 7.23o could be used. It
includes capacitance Cx between control valve and check valve. The causality of
the check valve Rx effect allows that flowrate through the check valve depends on
the direction and size of the pressure drop ~Px across it. Some appropriate
equations are included with Fig. 7.23c.
APx
QI, Px Pa
Ps
~1: '~ 0 ~:I: ~ 0 Q---'-~ TF
Ql Px Q2 Pa
Px = f (Qcx) Qcx = Q1 - Q2
Q2 = f (APx) APx = Px - Pa
Pa = f(Qca) Qca = Q2-Qa
Rv RX Ca Rv~ _Rx Ca
i iQca Ps
v SP/ Qca
Zz---~l ~----"'- 1 : - 0 --'~'a-" ' or ~'~ 1 : ~ 0 . . .Qa
.
Q1 Q1 Pa Q1 pa
~-------"
PY
Qy
"*'
I:
"'-- - ,TF ---~ k
SFg
Ay %, 4
\ k
k~& \ I
Cay \ .~
Y- Fa
~1~a-~ I
.... ~ Im
4 caz ~a
T
SPx I / ,, // I
k Ay ~k i
k~Cay k%
wl ~,Caz i )~ Xa 3- Xa
SP // / I
/
'I" Qz
~= 1 ," !-)~
paz . Qaz
. . . ~TF "- --""
Az
Rf
\ /
Pressure-actuated Valve Qaz = 0 if Paz < SPw
of Counterbalance Valve = Xa - Az if Paz > SPw
Pz = Paz - SPw
SPw = source pressure
down, the pressure build-up in the upper actuator line causes increased pressure in
the lower side and the counterbalance valve opens to allow discharge flow, and hence
downward motion of the load.
Fig. 7.24c shows a bond graph for the counterbalance valve for the condition
~at the load is being lifted. The check valve structure of Fig. 7.23d has been
adopted. That is, the check valve is seen as providing a particular form of
pressure source SPx . Supply Pz , Qz can be regarded as having come from the
Z port of control valve structure of the form of Fig. 7.18, while discharge Py ,
Qy is feeding into the control valve's Y port. The various equations appropriate
to Fig. 7.24c have all been developed previously.
Fig. 7.24d shows a bond graph for the counterbalance valve for the condition
that the load is being driven downwards. The pressure-actuated valve is set at a
certain threshold pressure. Thus, it can be shown on a bond graph with a symbol and
causality similar to that used for a check valve. In this case the R effect is
designated as the pressure source SPw . SPw will have a particular form which is
substantially independent of the system in which the valve is used.
It is probably possible to combine Figs. 7.24c and d into a single bond graph.
B0wever, the algebra of the equations would become confusing. It is probably better
to use the separate diagrams to prepare the joint equation set using the notions that
Fig. 7.24c applies when control Xv is negative, and Fig. 7.24d when Xv is
positive.
The more common components, actions, and phenomena relating to the dynamic
response of hydraulic control systems have been considered. The relationships
adopted or proposed are not necessarily unique. A particular analyst may prefer
different equation forms. Particular systems may require the use of additional
equations or of different equations to describe adequately what the analyst feels
will occur. With some experience in modelling, simulation, and in observation of
hardware system performance, the designer-analyst can develop modified bond graphs
and relationships. The objective of the present Chapter - - indeed of the present
Text - - is to provide an approach to modelling, enough bond graph material, and
enough equation forms, for simulation as a system design aid to be viable.
152
Quite apart from the deciding of equation forms, the selection of numerical
values for coefficients in the equations is not always easy. Component manufacturers
may have to be spec~fically solicited for some data.
7.13 REFERENCES
7.2 Davis, D. and Stecki, J.S., Effective Bulk Modulus in Hydraulic Systems,
Proc. Fluid Power Conf., FPS, Sydney, 1978.
7.3 Dransfield, P. and Bruce, D.M., Leakage Flow Rata Past Pistons of Oil
Hydraulic System Components, J. of Aircraft, AIAA, VS, N2, 1968.
7.5 Chong, F.K. and Dransfield, P., The Effect of Choice of Relief Valve on
the Response of a Hydraulic Control System, Proc. ist Aust. Conf. on
Control Eng., I.E. Aust., 1979.
7.6 Steber, G.R. and Romer, I.C., Simulation of a Pressure Relief Valve,
Simulation, V12, N5, 1969.
7.9 Svoboda, J., Analogue and Digital Modelling in the Design of a Hydraulic
Vehicular Transmission, Proc. I.Mech.E., V193, 1979.
CHAPTER 8
8.1 INTRODUCTION
Chapters 4 and 5 introduced the ideas of power flow and power bond graph
modelling respectively. Chapter 6 outlined the solution procedure once a model is
obtained. Chapter 7 presented and discussed bond graph structures for components
commonly used in hydraulic control systems. Application of the techniques will now
be demonstrated more adequately using a series of case studies escalating in complete-
ness and complexity. Where judged necessary, explanatory detail is provided.
Fig. 8.1a illustrates the system, which was used previously as an example in
Chapter 5 (Section 5.5.4, Fig. 5.12) . At that stage return line effects were
neglected, the model was restricted to positive spool displacement only, and it was
assumed that fluid was supplied to the valve at constant pressure. The example will
now be extended to include return line effects and to allow spool displacement in
both directions and thus load motion in either direction. The assumption of constant
pressure supply will be retained, though its validity is questionable.
Fig. 8.1b shows the components considered to be relevant to the present analyis.
Fig. 8.1c shows the bond graph, formed by direct coupling of the bond graphs of the
major components, taken from Figs. 7.18 for the valve (with Rse neglected), 7.19b
for the actuator (with Rfa and la deleted as they will be lumped with the load,
the actuator and load being rigidly coupled; and internal leakage Qa neglected) ,
and 7.5(2) for the load. Rsy , Rsz , Rye , and Rze have each been shown as Xv
a. The System
Xm
S xv i
//III~I
C. ~e Bond Graph PY
Rf
/ ~xv) Rye(xv) ~
sPs
Qs ~ 0 sPs sPe0:~ ~ ~m I
Xm
~z(xv) ~e(xv) / F/
Pz
I \ I
4-Way Control Valve Linear Actuator Driven I~ad
Cy and Cz b e i n g s h o w n as Xm d e p e n d e n t v i a d i r e c t feedback of Xm .
Fm = Fy - Fz - F f - Fk - SFe 8.2.1
the Xv e q u a t i o n 8.1.30 b e c o m e s
the i n p u t e q u a t i o n b e c o m e s
Xi = to be s p e c i f i e d 8.2.3
a. The System
Pressure
Com~n
st~pated I'l ~lI
b. Simplified System
supply Oil Xi
ConstantPressure
! !
., I I
Ck
Xm
Rf
Co The P r o p o s e d Bond Graph
Supply Control Valve Actuator Load
/ ~ / \
~------Xi, ~ Py
~ o-.~ APye
~e/~'
I
Qs--~0sPs
$7~~-~ 7-,_t,......~n ~ -~
\---Xi, >~ Qz Pz
FIG. 8.2 BOND GRAPH FOR LOADED HYDRAULIC SERVOSYSTEM
159
Many systems also include a high pressure filter near the pump discharge. Some
include an accumulator. It is reasonable to think of a power unit as being versatile
in that it could be used to provide power to any of a number of actuation systems.
Hence, for the system design situation it is attractive to have on hand dynamic models
for typical power supply sub-systems.
Fig. 8.3a shows one typical pumping system. The variable restriction is
included to represent the impedance load the sub-system will "see" when it is
connected to an actuator sub-system. The intention now is to prepare a fairly
detailed model which could be useful in a variety of applications of the pumping
system. Such a model could also be used to study the internal dynamics of the
pumping system itself.
Previous Figs. 5.4, 5.8, 7.9 developed a bond graph for an induction electric motor,
including its output shaft compliance. Fig. 7.9 is repeated as Fig. 8.3b. The
coupling and pump shaft unit can be included by lumping its inertia with rotor inertia
Ix , its friction with Rf and its compliance with Ck . From previous Fig. 7.11
the bond graph of Fig. 8.3c can be proposed for the pump. Fig. 8.3d shows a bond
graph encompassing the relief valve (Rr) , the resistance of the line plus filter
{Rb) , the capacitance of the line plus filter (Cb) , and the resistance of the load
valve (Rv) .
Fig. 8.3e shows the system bond graph formed simply by coupling the bond graphs
of Figs. 8.3b, c and d. Adjacent 1-junctions and adjacent 0-junctions have been
combined. The direction half-arrows and the causal bars were readily added using
physical sense and causal preferences. Combining the component sub-systems did not
upset their power flow directions and causalities.
The system equations written using the bond graph are given in Table 8.2. The
terms and symbols are as previously used.
Linear relationships were assumed for most of the dynamic effects in the model.
This assumption allows effects such as Ix , Ip , Ck , Cb , Rm , Rp , Rr ,
Rb , and Rv to be taken as linear coefficients. This is by no means always
160
The friction equations used for Tf and Tfp were discussed in Chapter 7.
Coulomb friction could be added.
pump torque = volumetric displacement per radian (Vp) - pump pressure (Pp)
pump flow = volumetric displacement per radian (Vp) pump speed (up)
For a variable displacement pump, volumetric displacement per radian of rotation pez
radian of swashplate angle ( V ~ ) multiplied by the swashplate angle (8p) in radi~s
gives the required pump volumetric displacement (equation 8.3.24).
The model, equations 8.3.1 to 25, was applied to a specific system, which was
also built and tested experimentally (Refs. 8.1, 8.2). A fixed restriction Rv was
used, and the system was activated for dynamic response by moving swashplate angle
8p from one position to another, increasing the flowrate demand of the system.
Some of the coefficients needed for the model were easily obtained, others
required experimental investigation or estimation. The values and how they were
obtained are indicated in Table 8.3.
Fig. 8.4 shows the results of a simulation of the model effected using CSMP on a
B6700 computer. Fig. 8.4a shows the input used, a change in swashplate angle.
Fig. 8.4b shows the consequent development of system pressure Pv(t) , together with
a corresponding experimental measurement. Fig. 8.4c shows simulation results using
a variety of values for the bulk modulus of the oil. Once a simulation is working
it is a simple matter to make this sort of excursion to investigate the sensitivity
161
a. A Pumping System
Variable Delivery Pump Restriction
Ele Ill I
~~ PressurFil~iief Val~e I .... I
TfPI T Qp
o ='I: ~ " ~ ~ 0 : ~ .
Ty~ Vp pp
Ip
~Qr ~ I T
~--~0:
Pp
"I,
Qpy
~0
Pv
'*:1: ~
Qv
'=
~-~0--~I~-L0--~I~ ~ ~K-L0~--~I~-~0----~I: ~
Ix Ip Rr Cb
/
Iml
0.i
time
I I I : -~
0.i 0.2 0.3 second
b. Response of Pv
Pv, MPa
2O
15
I0
Experimental
Initial Simulation
time
. i i I I
0.2 0.4 0.6 0.8 1 second
lO
J - - B = 1380 MPa
5 / / S/ .... B = 1035 MPa
p u m p d i s p l a c e m e n t can b e zero or p o s i t i v e o n l y
Fig. 8.5b shows a b o n d graph for the system. It u t i l i z e s a structure for the
p u m p m o d i f i e d from Fig. 7.11g b y n e g l e c t i n g p u m p f r i c t i o n (Rfp) and i n e r t i a (Ip) .
It i n c l u d e s d e l i v e r y line r e s i s t a n c e Rh and h y d r a u l i c m o t o r leakage Ra . It
n e g l e c t s the e f f e c t s o n d y n a m i c response of the r e t u r n line flow. Pump, d e l i v e r y
line, and m o t o r c a p a c i t a n c e s are l u m p e d as Cpha . The h y d r a u l i c m o t o r is regarded
as a f i x e d - m o d u l u s p o w e r transformer, w h i l e the d r i v e n a s s e m b l y is r e p r e s e n t e d as a
r o t a r y i n e r t i a load w i t h friction. The set of p o w e r state e q u a t i o n s p r e p a r e d using
the b o n d g r a p h is shown on Table 8.4, the nomenclature b e i n g i n c l u d e d w i t h Fig. 8.4.
a. S c h e m a t i c of S y s t e m
S~op
Rp Ra Rfj
~ Q E P / Q r a
I Qa TfJ1
Pp
Sup
Pa Va mj
/
/ Q AP
/
Xi Cpha
c. D e f i n i t i o n of Symbols
Effect Equation
P~er S t a t e Equations:
S u p p l y source Sago = constant 8.4.1
Driven inertia Ij ~j = i/Ij JTj at + ~j(o) 8.4.2
J
R e l i e f valve Rr Qr = 0 if Pp ~ P s e t 8.4.5
= Kr (Pp - Pset) if Pp > P s e t
Line r e s i s t a n c e Rh dPh = Kh Qpa 8.4.6
M o t o r leakage Ra Qa = K % a Pa 8.4.7
Load friction Rfj Tfj = Kfj ~j 8.4.8
Pump transformer Vp Tp = Vp - Pp 8.4.9
Qp = vp. s~p 8.4.10
Motor transformer Va Ta = Va - Pa 8.4.11
Qa = Va - ~j 8.4.12
Auxiliary Equations :
C a p a c i t a n c e of pump, hose
Cpha = Vpha/B 8.4.17
and actuator
P u m p displ a c e m e n t / r a d Vp = Kp Xi 8.4.18
Control i n p u t Xi = i n p u t to b e s p e c i f i e d 8.4.19
The b o n d g r a p h i m p l i e s a q u i t e d e t a i l e d d y n a m i c m o d e l , as it includes
a p r e s s u r e c o m p e n s a t e d p u m p w i t h its i n p u t d r i v e
a r e l i e f valve
a filter
a 4 - w a y d i r e c t i o n a l control valve
a c r o s s - p o r t r e l i e f valve system, w i t h m a k e - u p supply from th~ charge pump
a hydraulic motor
the d r i v e n rotary i n e r t i a load.
167
-[ TI
c. non4 Grsph
-,,T/ ~
,,
Lj
i -~,i b0 "I' ~
~ . . . . 'P I
## ~3 II" ( w )
L; O~
I/
O"
ill
\ .............. 1%_ _ - J %. . . . . . . . . . . . . . . . . / % . . . . . . . . . . . . . . . . . . /
p~l~ ~ , ~ / control c~n Port ~ e t sylz~ S~li'luZ~= ~ = ~ ~,la
The additional information and model development on a bond graph relative to that
on a power port schematic is clearly illustrated. Each has a block representing the
system's c~mponents or sub-systems. Each shows the main power flow bonds. But that
is all that the power port schematic shows, while the bond graph goes on to expose
those effects in each block which are necessary, or that the analyst decides are
significant, to the system's dynamic performance. Further, the decisions on
causality provide both a technique for checking the appropriateness of a developing
model, and a framework for arrangement of the equations. The full arrows on the
power port diagram and the half arrows on the bond graph both indicate the presumed
directions of power flow in the system, and provide equation sign conventions.
Equations for most of the system's components have been presented in previous
examples and will not be repeated here. The cross-port relief block has not been
detailed previously. Using the bond graph its actions can be described by the
equations of Table 8.5.
Effect Equation
F l o w t h r o u g h check v a l v e Qyl = 0 if &Pl p o s i t i v e ~ 8.5.1
R1 from Y line, and = K1 API if AP1 n e g a t i v e J
the a s s o c i a t e d 0- a n d
API = Py - Prl 8.5.2
1-junctions
Qy = Q y - Q y l + Qy3 8.5.3
F l o w t h r o u g h check v a l v e Qz4 = 0 if p z ~ pc
8.5.18
R4 into Z line, a n d = K4 . AP4 if Pz < Pc J
the a s s o c i a t e d i- a n d
AP4 = sPc - Pz 8.5.19
0- j u n c t i o n s
8.7.1 Introduction
were
170
control is to be manual;
Fig. 8.7c shows the main components of the hydraulic system proposed. It was
decided to use an open-centre control valve. The high gain between valve motion and
control flow in the valve's central position allows good control of load speed with-
out the need for constant pressure supply. The counterbalance valve in the lower
control line is needed to hold the load against gravity. The pressure relief valve
in the upper control line is used to limit cylinder pressure (and hence force and
acceleration) when the load is to be driven downwards, gravity providing much of the
force needed for the downward stroke.
to select the pump and cylinder ....... this was associated with the selection
of beam lever ratio, with requirement to use system pressure maximum of
13.8 MPa, and with the desire to use the smallest pump and cylinder which
would provide the speed of load motion required (overall 0.6 m/s) ;
,,
(
~ i Actuato~
y~aullc
0.15 O.B5 1 o
/ f
ELI/
~13.8MPa~ Cont~oI
- Va/ve
ntxe
d. Bond Graph
C~
~'
~Y ......
I
h0 "~' "
LS ~
3 ~ ~ = ~-~0 ~-~-~0 ~ " 0 ~
~z sFg
T "
~Az TM
FIG. 8.7 LOAD L I F T I N G SYSTF~4 (EXAMPLE 8.7, TABLES 8.6 AND 8.7)
172
main system pressure relief valve setting (as specified) 13.8 MPa
This brought the system design to the stage where a reasonably accurate assessment of
its dynamic response was needed to see if it satisfied the speed of response, and the
acceleration limiting specifications.
Fig. 8.7d shows a bond graph for the system developed on the bases that
the resistance of the short hydraulic lines have a negligible effect on load
position response
the capacitances of the lines can be lumped with Cp for the pump-valve
lines, and with Cy and Cz for the valve-actuator lines
the beam can be assumed to supply a constant lever ratio between actuator
displacement and vertical displacement of the load Xm .
The main purpose of the model is to allow prediction of the time response motion
of the driven load following cca~mand inputs applied to the control valve. The bond
graph of Fig. 8.7d allows development of equations for modelling the system for both
lifting and lowering of the load. AS seen in Section 7.9.5, many equations are
required for comprehensive modelling of the open-centre valve. In the present
system, lowering of the load is greatly assisted by gravity with little hydraulic
power requirement, the main need being to control the speed of descent. Thus
lifting of the load with hydraulic power is the main dynamic aspect. The equations
of the model are given in Table 8.6. They are written for lifting only.
8.7.4 Simulation
Fig. 8.8 shows the input and the responses of some of the system variables when
the model was simulated with MIMIC on a CDC Cyber 6400 computer. Load displacement
Xm(t) is directly related to piston displacement Xa(t) by the beam lever ratio.
Corresponding results obtained from the hardware system, which was assembled
8u]~sequent to the simulation, are included on Fig. 8.8. Considering the ambiguity
of some model parameter coefficients, the correspondence is good. It is a simple
matter to improve the correspondence between simulated and actual system performance
by adjustment of some of the model's parameters.
Comparison between Figs. 8.7b and 8.8b shows that the specified response is
approached but not achieved. The main discrepancy is the hump near t = 0.5 s .
Some parameter adjustment, possibly using optimization procedures of the type to be
discussed in Chapter 9, is called for. with the simulation working, it is a simple
matter to study the effect on response of adjusting a parameter of the model.
8.7.5 Conclusion
8.8.1 Introduction
Hydraulic control systems are commonly used where the precise positioning of
large resisting loads is required. High pressure hydraulics has little logic
capability, being a power transfer technology. Where logic or feedback is required
the hydraulic power system usually has to be interfaced with manual, mechanical,
pneumatic, electronic, or electrical control logic. Electrohydraulic servo-systems,
174
Internal leakage, Rp
QZp = KZp - Pp 8.7.3
effect
Relief valve, Rr effect Qr = 0 if Pp ~ P s e t
8.7.4
= Kr - (Pp - Pset) if Pp > Pset
Tp = Pp vp 8.7.6
Control Valve
F o r all Xv , 0 to -+ x v m a x
Xv is spool d i s p l a c e m e n t
f r o m centre. ~ = Qsy + Qsz 8.7.8
Fy = Py Ay 8.7.27
Qzl = Xa Az 8.7.30
Note i: W h e n the c o u n t e r -
b a l a n c e v a l v e opens, the
Rz effect becomes
negligible, Pzl and
Pz identical, and Czl
and Cz2 can be l u m p e d
as Cz .
Note 2: Equations 8.7.28
to 35 h o l d also for Xv
n e g a t i v e so l o n g as the
c o u n t e r b a l a n c e valve is
s e e n to b e e i t h e r full
open or full shut
d e p e n d i n g on the v a l u e
of Pzl and w h e n o p e n
h a s n e g l i g i b l e resist-
ance to flow.
Auxiliary Equations
Xm d t + Xm(o) 8.7.43
a. Input xv(t)
o V
I ! I I ~ I t I I ~ Q a
0 i
0.15
0.1
Xa, m
measured
0.05
simulated ..+.........
f
- tim, s
0 1
time, s
0 1
14
Pzl, MPa
;:.IJ -- time, S
0 0.2 0.4 0.6 0.8 I
long used in military and space technology, are now widely utilized in industrial
applications such as machine tool control, flight simulators, structure fatigue
testing rigs, computer-aided manufacture, and industrial robots. With increasing
use of microprocessors and computer control it can be expected that electrohydraulic
servovalves will be used increasingly in industrial and manufacturing automation.
Electrohydraulic drives are often used where the dynamic response requirements
of the driven load(s) in terms of speed, cycling, accuracy, and stability, are
particularly severe. One such application, the subjec~ of the present Example, is
the control and drive system associated with a 28 tonne seismic block facility.
The intention is that the electrohydraulic system shakes, at controlled force and
displacement amplitudes and frequencies, a structure or machine undergoing
vibrations analysis. The present system is intended to be useful to 300 Hz.
Fig. 8.9a illustrates the shaker system layout and Fig. 8.9b the essentials of the
hydraulic system. The hydraulic power supply centres on a 30 kW Denison piston
pump operating at up to 21 MPa set pressure and up to 60 /min flowrate.
A 30 kW induction motor drives the axial piston pump at i000 rpm. The main
pressure control valve is set at between 2.76 and 20.7 MPa, depending on require-
ments, to provide constant pressure at that value to the servovalve. The return
line pressure control valve is set at between 2.1 and 3.45 MPa, depending on the
system performance requirements, to prevent the possibility of cavitation with
consequent dynamic performance uncertainty. Both pressure control valves are
pilot operated by remote control from the servo-controller. The electrohydraulic
servovalve receives its position feedback modulated input signal from the servo-
controller. The water-cooled oil cooler enables oil temperature to be controlled,
thus controlling oil viscosity, an important factor in high performance hydraulics.
The long hydraulic hoses provide flexibility so that the actuator can be taken to
any desired position on or around the seismic block (approximate dimensions 3 m x
4.5 m area x 0.65 m depth 28 tonnes mass) . The accumulator package provides the
smoothing of pressure transients and the capacitive effects needed for smooth load
position control. The electrohydraulic system has various override, safety, and
functional features not included on Fig. 8.9 as they have little effect on dynamic
performance under normal system operation.
In order to better understand the dynamic nature and performance of the system
beyond the rather meagre information supplied by the manufacturer, a simulation of
the shaker system was required. Because the main interest lay in employing the
upper limits of dynamlc response, the model for simulation needed to be detailed,
accurate, and versatile.
179
a. ~hmlat~c ~a~t
+" +" I l
I,+'.I)Ut SI~nal'*---ll,,1- ~ } ~$I~II I . l
I - I
1 ] 'l
+.- . . . . . . . . . . I /
_/
b. Ellentlal# o f H y d ~ a u l i c Sylte~
CICU~K3~tI~I~
S*+~I++
. ]~:.z~. C,aph
/
mp
"+-/
F/ire:
\ 2-n . . . . . ,,
f
.1.. F~ 1' n....++./ Roe 1
\I
lllltxoh~*mL~+c T m l m l l l . In~ toId
Py
"\2 I +\,/
+T ? ? - i
r ' ~ l ~+'+'~ 0 - " * + I + - - ' ,
r~
1!_! + i~ / + , + +\
*/
]
0 +~+-i 1 t.,..-- 0 J.,,L~ I p.+.-- 0 }+..--+ O .,+--.i 1
p.iJr.) Itrm(R)
, ~ J+r,
X#
; j,J,~t +,+F+.'~+++++ j,
++
cooler
control
Fig. 8.9c shows the power bond graph structure developed for the system, the
various hardware components being indicated adjacent to the appropriate sections of
the bond graph. Some of the R and C effects included will not be simple
coefficients. Discontinuous and other non-linearities are present.
The actuator and driven load sections of the bond graph are similar to those
used in a number of earlier examples. The return line bond graph includes
recognition of its hoses, its accumulator, and its pressure control valve (Rr2) .
The return line from the EHSV joins with the discharge from the supply line pressure
control (Rrl) , and the total discharge passes through the oil cooler to return to
tank. The capacitance and resistance effects of the oil cooler and its associated
lines are included on the bond graph as Cx and Rx .
181
AYE
I I !
I
Signals ---P"
I "0~
x'/
",0 01 Xa I ~ Im
~I:
I T I
"0 ~I' I
/
FIG. 8.10 DIRECT INCLUSION ON BOND GRAPH OF MODULATING SIGNALS
8.8.3 Conclusion
Further explanation of bond graph development, the equations derived using the
bond graph, and results of a simulation of the system are given in Refs. 8.6 and 8.7.
There are 77 equations in the model, most of forms already used in this Text. It
should not be thought that the system bond graph of Fig. 8.9c came easily, without
effort. To form and link the component bond graph structures took considerable
effort in satisfying physical reality, adequate completeness, and causality.
However, development of the bond graph provided the key to orderly and progressive
model development in a situation which would have been daunting if approached by
conventional modelling procedures.
8.9 CONCLUSION
A number of examples has been given of the use of bond graphs to provide a
link between a hydraulic control system, real or proposed, and its simulation for
analysis and synthesis purposes.
182
8.9 REFERENCES
8.2 Dransfield, P. and Winton, R.D., Applying Power Bond Graphs to Hydraulic
Control Systems, Proc. III Congresso Brasiliero de Engenharia Mecanica,
Rio de Janeiro, 1975.
8.3 Stecki, J.S. and Cuthbertson, S.C., already cited as Ref. 4.3.
8.5 Barnard, B.W. and Dransfield, P., already cited as Ref. 7.4.
9.1 INTRODUCTION
The ethos of this Text is that ensuring good dynamic response should be an
integral part of the design of a proposed hydraulic control system. Further, that
the securing of this aim is b e s t effected via mathematical modelling of the proposed
system and digital simulation of the model. Further still that bond graph techniques
provide a dynamic modelling procedure particularly suited to the modelling of hydraulic
control systems.
To here the Text has concentrated on developing the bond graph procedures needed
to model hydraulic control systems. It has been implicit that such models are best
simulated by using digital simulation.
As seen in the examples of Chapter 8, a system's model will have many coefficients.
To retain satisfactory power and motion characteristics m o s t of the coefficients will
not be able to b e varied much w i t h o u t basic re-design o f the system. However, there
are usually a few coefficients for which the designer-analyst has a reasonable choice
of values. The idea becomes to seek that combination of values for these parameters
which yields optimum response.
In the simplest case only one parameter is varied. A set of discrete solutions
of the system model, with the chosen parameter varied over its allowed range, quickly
leads to the optimum value of the parameter; that is, that value w h i c h produces the
desired system response, or the possible response closest to the desired response.
In the more likely case where two or more parameters can be varied simultaneously,
the number of possible combinations of values of the chosen parameters is so formidable
that computer techniques are required for the search for the optimal combination of
parameter values .....again, that cc~bination w h i c h provides response closest to that
desired.
This Chapter describes a response optimization scheme that has been developed
and used with hydraulic control system response (Refs. 9.1, 9.2). There is no
attempt to generally survey control system response optimization.
t
~E =
i o
li~ dt
L a B r o o y [9.1] c o n s i d e r e d m o s t of these in r e l a t i o n to r e s p o n s e s t y p i c a l l y r e q u i r e d
of h y d r a u l i c control systems. It was confirmed o r c o n c l u d e d that
/ x
f l /
~11 \Desired
- ..~
r o r E{t)
response xl(t)
time
ISE
Its E " E dt
o
ItSE
I ts t E E dt
o
IStSE
I ts t t E E dt
o
z~
I ts IEI dt
o
ItAE
I ts t " IIE d t
o
ts 2Kt
IXSE
i o
e E E dt
IMXSE C[1 o
- e -2K E E dt
I d e n o t e s 'integral of'
S d e n o t e s ' square of'
E denotes error
t d e n o t e s time
A d e n o t e s 'absolute v a l u e of'
X d e n o t e s 'exponential w e i g h t i n g of'
M d e n o t e s 'modified v a l u e of'
ts is s i m u l a t i o n time (> s e t t l i n g time of m o d e l response)
e is the e x p o n e n t i a l
K is the a r b i t r a r y c o n s t a n t u s e d w i t h IXSE and IMXSE
187
a Unimodal
Value of integral
b. Mul timodal error function
c. Non-modal
Value o~ parameter
Hydraulic control systems are stiff and are u s u a l l y h e a v i l y loaded. They are
o t h e r error functions.
9.4 .i Introduction
designer-analyst would rarely know which of the three s i t u a t i o n s o f Fig. 9.2 he has.
The Golden Section Search Procedure for 1-dimensional optimization does not
effectively extend to multi-parameter use [9.5]. More practical search techniques
are based on the evolutionary operation method of Box [9.6].
The Simplex method [9.7] is one EVOP approach. A Simplex is defined generally
as a shape formed by (n + 2) points distributed in (n + i) dimensional space, n
being the number of optimization parameters. Error function evaluations are made
for each of the (n + 2) points of the Simplex. The point producing the worst
evaluation, i.e. that point furthest from the point producing the smallest evaluation,
is discaxded and replaced by a point located in a more favourable position. This
point is located using a set of mechanical rules intended to help the Simplex converge
towards the point of minimum possible error evaluation. It has been reported [9.8]
that Simplex is most satisfactory when the error function is linear and when no
constraints are placed on the range of each optimization parameter. Also that
Simplex is susceptible to collapse into a sub-space to yield a local rather than the
global optimtun. For these reasons, Simplex seems to be inadequate for optimization
of the response of hydraulic control systems. It is recognized that Simplex
techniques might have been improved to accommodate its earlier limitations.
The Complex search method proposed by Box [9.9] and improved by Guin [9.10]
appears to be well suited to the present requirements. Complex is heuristic (self-
learning), is capable of handling non-linear error functions with unequally constrained
optimization parameters, can negotiate hard constraint boundaries, and can turn
corners in the space being scanned. It is unlikely to collapse into a sub-space to
yield a local minimum, as a single point existing in a "better" region (i.e. closer
to the global minimum) can cause the complex to be re-directed towards this point.
190
9.4.4 Complex
The random nature of the casting of the initial complex of points allows a
probability that a point may be in the region of the global error function minimum
which locates the optimal combination of optimization parameters. In this case, the
Complex procedure will cause the complex to collapse towards the optimal point. If
no point of the initial complex lies in the region of optimality, the heuristic rules
governing the motion of the complex cause eventual scanning of the optimal region.
Like simplex, the rules of Complex generate new points to replace "worst" points with
the objective of collapsing the complex towards the global optimum point. Like most
search procedures, a wide choice of error function is allowed.
Fig. 9.3 illustrates, with a simple 2-parameter example, the nature of Complex.
The ordinates show the boundaries of the two optimization parameters N1 and N2 .
The small-chain broken lines illustrate contours of constant values of the chosen
error function as N1 and N2 are varied. Adopting sigma = 3 (i.e. at least
n + 1 ), the three points of the complex are shown as N1 , N2 , N3 , and their
centroid is at N . To minimize the value of error function, the complex must be
directed from the i0 contour towards the 4 contour. This is effected by reflecting
the worst point, N1 , through and beyond the centroid to NI' The next-stage
complex becomes N2 , N3 , NI' It can be seen that N3 would suffer the next
reflection through the new centroid.
Fig. 9.4 shows the flow chart developed b y La Brooy [9.1] for compiling the
optimization algorithm COMPLEX which he used with MIMIC to provide an automated
optimization procedure. The flow chart was modified from that of Box [9.9] for
increased efficiency. A system model simulated in MIMIC calls on the sub-program
COMPLEX which returns a parameter combination via multiple RETURN statements for
MIMIC to effect a response solution. MIMIC then evaluates the error function for
the particular parametric combination and recalls COMPLEX, passing the error function
value and receiving another parameter combination in return. The sole purpose of
the digital simulation procedure is to provide objective error function evaluations.
COMPLEX implements the logic of Fig. 9.4 to enable the Complex procedure to function.
It is a discrete Fortran optimization sub-program containing several parameters such
191
N2
N2max
%.,. __ .~, ~ ' ~ " NI' ,~/7
N2 / ~'
N2min I l- N1
Nlmin Nlmax
9.5 EXAMPLE
K1 = 0.024
K2 = 0.3 E-3
was accepted as o p t i m i z a t i o n .
values respectively.
9.6 C O N C L U S ION
,,, ~ ,
81 ptlmum
/W##~ desired 1
; ~ t SeCS , : ~ t secs
0 0.2 0.3 0 0.I 0.2 0.3
0.i
e R e s p o n s e s for K1 a n d K2 a t T o p
and Bottom of Their Ranges
xo
8!
4-
//
/
ldesired
/
l
o.1
~<i -- o . o i ,
l
0.2
K2 = o.oooa
~ t secs
o:a'-
f.
Xo
Optimum response
Response Sought
F
if F a s t e s t
4 //~///~/
I/// top range response
Xo
f" , Lt secs
s~ 0 011 0:2 02a-
K1 = 0.05, K2 = 0.0004
| .... | ; ~- t secs
0.i 0.2 O.3
scheme for optimizing the response of hydraulic control systems has been described
and demonstrated. The component steps and procedures are:
9.7 REFERENCES
9.7 Nelder, J.A. and Mead, R., A Simplex Method for Function Minimisation,
Computer Jrl., V7, 1964-65.
195
9.8 Jacoby, S.L., Kowalik, J.S. and Pizzo, J.T., Iterative Methods of Non-
linear Optimisation Problems, Prentice Hall, 1972.
9.9 Box, M.J., A New Method for Constained Optimisation and a Comparison with
Other Methods, Computer Jrl., VS, 1965-66.
9.11 Dransfield, P. and LaBrooy, R., Hydraulic Systems with Precision Reflexes,
Machine Design, May 20, 1976. The example developed in this article came
from Fluid Power Control, Blackburn, Reethof, and Shearer.
CHAPTER 10
i0.1 INTRODUCTION
Phenomena which are inherent in working hydraulic control systems to the extent
that they are always present have been included in the modelling procedures dealt
with in this Text. Included in this inherent phenomena category are
Some can be allowed for to some degree, if necessary, in system simulation models.
It is better that components and systems be designed and operated so that parasitic
phenomena do not appear. To achieve this requires an understanding of the phenomena
and an appreciation of how each can be avoided, or its affects on system response
minimized.
197
With the exception of coulomb friction, which has been considered already, each
of the mentioned parasitic phenomena will be discussed briefly. References for
further study are provided.
i0.2 CAVITATION
To provide these characteristics, the fluid used in the system must retain a high bulk
modulus as it passes around the operating system.
If the pressure of a quiescent oil is changed slowly, air is absorbed and dissolved
(for increasing pressure) or comes out of solution and is rejected (for decreasing
pressure), as the natural saturation state appropriate to the particular pressure is
sought. The mechanism of rejection of previously dissolved air as pressure is
lowered is for bubbles of air to form and to float through the surface of the fluid.
198
This takes time, particularly with viscous fluids. Air is vastly more compressible
than is oil, and consequently lowering of pressure has far more effect on the volume
of the air component of the fluid than on the oil component.
The presence of dissolved air up to the natural saturation level has little
effect on the bulk modulus of an oil. Thus, air which remains dissolved in the
fluid will not affect the response of a hydraulic control system. However, the
presence of bubbles of air in the oil will render the fluid more compressible, as any
increase in pressure will quickly compress the bubbles of air. In these circumstances,
system response can be adversely affected. Taken to the extreme, if a large quantity
of air bubbles exists in the oil, the fluid suffers a drastic reduction in its bulk
modulus.
The fluid in a working hydraulic control system experiences frequent, large, and
sudden changes of pressure, in both directions. Thus, the saturation level of
dissolved air will be changing rapidly and the actual level of dissolved air will be
trying to adjust itself. This does not matter much, so far as system response is
concerned, while pressures are increasing, as the fluid will be sub-saturated.
However, when pressure is decreased rapidly, such as between a reservoir and a pump's
suction port, or downstream of the orifice of a control valve, the fluid can be super-
saturated and air coming out of solution as bubbles is a likely consequence. When
the pressure drop is from one positive value to another positive value well above
atmospheric pressure, this condition should not be serious. Nor should it seriously
affect system performance. The bulk modulus will be transiently lowered but will
remain high. If however the low pressure approaches atmospheric pressure or is below
it, the volume of air coming out of solution can be large, the bulk modulus of the
expanding fluid can approach zero, and system performance can be grossly degraded.
The formation of bubbles in the fluid due to air or gases coming out of solution
due to reducing pressure is known as gaseou~ er~uitation. In effect, the bubbles are
cavities in the liquid. Such cavities are highly compressible and cause reduction
of the fluid's bulk modulus, Once formed, gas cavities don't inu~ediately dissolve
if system pressure is again increased. It takes time for re-dissolving t o occur.
Thus a reduction in bulk modulus arising from gaseous cavitation can p e r s i s t
temporarily even if system pressure rises.
Quite apart from the natural gaseous cavitation described, a rapidly moving fluid
undergoing pressure changes, directional changes, and temperature changes in the
presence of air or other gases can entrain bubbles of the gas. This is likely to
occur for example if return lines to a reservoir are not drowned, or if there are air
leaks, n o matter how minor, in uncharged pump suction lines. Entrained bubbles
affect the fluid just as do bubbles of previously dissolved air. Thus the term
199
Air entrained at low pressure can dissolve into the oil when pressure is increased,
raising its air content level towards the saturation state. This is highly
undesirable as the potential for gaseous cavitation, when pressure is inevitably
lowered at a later part of the hydraulic circuit, is increased. The oil in a well-
designed, well-operated, well-maintained hydraulic control system will tend to become
sub-saturated. Any air entrainment will negate this desirable condition.
Like any liquids, hydraulic system fluids will boil (vapourize) if either the
temperature is high enough or the pressure is low enough. Vapourization is
characterized by the development in the liquid of vaporous bubbles which in effect are
cavities in the liquid. Each liquid has a pressure, at which vapourization will
develop. If pressure is reduced to or below this vapour pressure, vapourization will
occur. A liquid's vapour pressure is temperature dependent, increasing as temperature
is increased. Thus, fresh water boils at 100C at sea-level atmospheric pressure
0.i MPa. That is, at 100C, its vapour pressure is 0.i MPa. At temperatures below
100C, the water's vapour pressure is below 0.i MPa and if atmospheric pressure is
maintained water re/nains in its liquid state. Water at 21C will vapourize if its
pressure is reduced to 0.0024 MPa, a strong vacuum state. Hydraulic fluids exhibit
similar characteristics. Mineral oils have vapour pressure characteristics lower
than those of water and are consequently less susceptible to cavitation. Water-based
hydraulic fluids have vapour pressure characteristics more like water.
P - Pv
Nc = y p ~2
Low Nc indicates a high tendency to cavitate and vice versa. Hence high Nc ,
which results from low Pv , low X and low p , is desirable.
pressure of the oil falling below its vapour pressure, whereon the oil
boils (vaporous cavitation).
In working hydraulic control systems, gaseous and vaporous cavitation can co-exist.
Gaseous cavities are likely to form as soon as system pressure is reduced below the
appropriate saturation value~ which can be a quite high pressure. Vapour cavities
will not form until pressure drops below the appropriate vapour pressure which is
normally around or below atmospheric pressure.
The presence of cavities in the working fluid reduces its bulk modulus. If the
reduction is dramatic, system response will be degraded.
The presence of dissolved air, and its coming and going as pressure changes, is a
natural phenomena. Little can be done about it. Unless accompanied by vaporous
cavitation, it is not usually of serious consequence. Entrained air is always un-
desirable and can seriously degrade system performance. Good reservoir design
and good suction line design can minimize the presence and effect of entrainment.
Vaporous cavitation can be avoided by keeping pressures at all parts of the system
above the fluid's vapour pressure. If necessary, reservoirs can be pressurized
and/or main pump suction lines charged.
Cavitation not only degrades system response and performance. It can also damage
the surfaces of components due to very high instantaneous pressures generated as
cavities implode. It can also damage the system fluid by oxidation due to the
very high local temperature generated as bubbles of air or gas or vapour are
compressed when pressure increases. Cavitation also degrades a system's
performance by generating offensive noise.
Refs. i0.I to 10.3 provide more detailed discussions of cavitation and its
relevance to hydraulic control systems.
Fully-developed vaporous cavitation c~n cause the bulk modulus of the system
fluid in the region of cavitation to be zero. Under these conditions the pressure
in the cavitating region is effectively zero and no power will be being transmitted
from the region. The zero bulk modulus, zero pressure, zero power condition applies
only to the region of cavitation and is usually highly transient.
201
associated with the computed negative pressures. That is, there is no recognition
of cavity formation. This could lead to errors in load response predictions, though
the errors would probably be transient in most situations.
McCloy [10.4] provides an approach for allowing for cavity formation in simulation
models of hydraulic cylinders. He also suggests that the model complication required
is rarely justified on the basis that its affects on cylinder response are minimal.
Rogers [10.5] found it necessary to include modelling of cavity formation in simulation
studies of hydraulic backlash, to be discussed in Section 10.3.
One of the most common applications for hydraulic control systems requires the
moving of an inertia load from one state (position or speed) to another in a
specified range. There is usually the requirement or expectation that the driven
load will be restrained accurately and stiffly in its new state prior to the next
instruction, which might be milliseconds, seconds, minutes, or hours later. For the
case of linear positioning, this objective is usually sought by using a hydraulic
actuator whose output is rigidly connected to the driven load and which receives oil
from a closed-centre or tandem-centre control valve. The intention is that when the
control valve is centred, both hydraulic lines joining the valve to the actuator will
be, and will remain, full of oil in a high bulk modulus state, thus stiffly locking
the load in its required position. With rotary drives, the intention is that both
the supply line and the discharge line of the rotary actuator be filled with oil of
a high bulk modulus and which is flowing at a constant rate.
202
Unfortunately, these states are not always achieved, and the driven load can
sometimes be easily moved about its desired state. This effect is called h y ~ l i c
backlash. Its basic characteristic is illustrated on Fig. i0.i. Ideally, any
attempt to move the located load should immediately induce a high resisting force or
torque from the hydraulic drive. If an atten~t to displace the load meets zero or
little resistance, at least over a band of displacement, hydraulic backlash exists.
In the circumstances described, hydraulic backlash can exist only if the hydraulic
fluid in one or both of the hydraulic lines of the actuator has a zero or very low
bulk modulus. It is assumed that the hydraulic components have normal resistance to
internal leakage and that the mechanical cc~ponents of the system and load are of
normal rigidity.
Not all hydraulic control systems are susceptible to the development of hydraulic
backlash. The phenomenon will develop in susceptible systems only under some
operating conditions, notably during rapid braking of inertia loads by valve or pump
displacement control. It tends to be transitory in nature and usually decays in
time after a load reaches equilibrium, though this time can be of significant duration.
, /
/
!
Displacement
applied to load
/
W
rl
L.~l
tl-,,~
.J
It
Q = Kv - Xv P
Qr = 0 if P ~ Pset
= Kr - (P - Pset) if P > Pset
The thrust of the Text has been not the accurate simulation of individual hydraulic
components, but rather the use of simplified representations of them which allow
reasonable predictions of the major system response characteristics such as response
of the driven load, and pressure developments in the system.
valves experiencing flow develop force reactions on their moving component due
to the flow. These parasitic flow forces are undesirable and can be minimized by
good design. As they may cause or account for unexpected system performance they
should be appreciated.
The origins and nature of flow forces can be appreciated by considering the
simple spool configuration of Fig. 10.3. In the centred state, the supply pressure
acts on the two internal spool land faces equally and the pressure forces on the spool
sum to zero. If the spool is moved to the right, as shown, the left control port
remains blocked while the right control port experiences flow. The left internal
spool land face continues to experience stagnant fluid at uniform pressure. However,
the fluid flowing at relatively high velocity through the right-hand port gains
kinetic energy, and this can only be achieved at the expense of pressure (potential
energy). Thus, in the region of the jet formed by the partly-opened orifice, the
pressure acting on the internal spool land face is reduced as shown. The axial
force due to pressure on the right-hand spool land is reduced below that on the
left-hand and a net axial force exists on the spool. The axial flow force acts to
close the valve. If not resisted, it will do so. With manually controlled valves,
the manual operator, possibly aided by a detente mechanism, usually can overcome the
flow force and hold the valve open. It may become a problem with large valves.
Flow forces axe more significant with electrohydraulic valves, as their axial
205
Most of the major texts on hydraulic control systems provide details (e.g.
Refs. 10.7, 10.8) allowing calculation of steady-state flow forces. Dynamic flow
forces are less readily assessed.
Piston pumps and motors and spool-type valves rely only on the close fit of
piston to bore to minimize leakage past the piston due to axial pressure drop. The
fine clearances needed demand near-perfect geometry of both piston and bore, which is
a main reason why hydraulic piston and spool components are expensive to produce.
viscous drag force, in their axial direction, due to viscous leakage flow
in the clearance between piston and bore, and due to viscous shear of the
oil film due to axial motion of the pistQn,
mechanical friction due to any influence which causes the piston and bore
to be in contact instead of the piston "floating" in the clearance oil film.
The first four of these effects are inherent and are well understood. Single-
ended piston devices (pumps, motors, some pressure control valves) experience large
but expected end pressure forces which are mechanically reacted in the device. Spool
pistons axe usually symmetrical so that the net axial pressure force is zero. These
devices have the potential to be operated by a miniscule axial force, which has only
to overcome the inertia and the viscous drag forces. A piston 1 cm diameter, 2 cm
long, with diametral clearance 0.01 ram, and with axial pressure drop of 20 MPa
experiences a drag force of only about 2 N . Valve flow can induce axial and
radial forces on valve spools, inhibiting their performance and response, as
206
Fx
i?
,x ...i,...
if/M////,////
[ ~
#
~, -
I
,
I
Fflow
/i//-
I / / I / / / / ///
////////
C. Bent Spool
Lateral force
FIG. 10.4 SOME CAUSES OF UNSYM-M~TRICAL CLEARANCES
IN WORKING PISTON COMPONENTS
207
Mechanical friction requires that a lateral force exists to force the piston into
contact with its bore. Then, any attempt to move the piston axially must overcome
the resulting friction, and the force required to induce axial motion would include a
component
Ffx = p Fz
with piston pumps and motors, the eccentric nature of the drive acts to induce
piston-bore contact as the piston is forced against the bore (Fig. 10.4a). With
pumps, the high speed reversing motion of the piston allows little time for lateral
and tilting displacement of the piston to fully develop, and also encourages good
hydrodynamic lubrication to be sustained. It takes time to squeeze out a film of
mineral oil. In any case, the large power input to a p u m p completely overwhelms its
piston friction effects, and piston friction is reduced to a negligible effect as far
as dynamic response of the driven load is concerned. This is not so with hydraulic
motors unless they too are run at reasonably high and constant speed. The well-
known tendency of hydraulic piston motors to stick or to be jerky at stall or slow
speeds is at least partly due to friction deriving from lateral piston displacement
and tilt arising from the eccentricity of the drive.
Fig. 10.Sa shows a piston lying centrally in its bore with a high pressure Pl
at one end and a lower pressure P2 at the other end. Ideally, one would imagine
that the pressure distribution in the clearance is linear as shown, and that it is
the same along all axial generators on the piston surface. That is, the axial
pressure distribution at 8 = 0 is the same as that at 8 = w and at all other
values of @ , and the radial forces generated on the piston surface sum to zero,
leaving the piston floating freely in its surrounding oil film.
208
o
e I
/ / Z z ] ~' f / W / W / / / / / /
P2 P1
Fz = 0 FZ = 0
f / / / / / i /
/ / / / / / / ~ / / / / / / / / / /
/ / / / / / / f / / I / / / / / / / / /
// / / ~" / / // /
// W z" / // / / /
/ - / / / / / / / /," / / / / / / / / / /
Friction between piston and bore arising frcm uneven axial pressure distribution
around a piston is known as hydraulic lock.
Under what conditions can unequal axial pressure distributions ~ radian apart
arise? A n d under what conditions will the n e t lateral force of such pressure
distributions cause decentring of the piston towards contact with the bore? These
questions have been studied fairly extensively, for example in Refs. 10.9, 10.10,
I0.II. If for any reason, the actual axial clearance between a piston and its bore
is diverging in the direction of axial pressure drop, hydraulic lock can develop.
The attitude of a piston in its bore can rarely be guaranteed. Even in the
ideal case of Fig. 10.5a, why should the piston be concentric with its bore? Its
own mass can cause it to sink into the clearance oil film. The piston's supporting
mechanism m a y cause it to tilt (Fig. 10.4a and b). A slight banana bend in a long
spool can cause the situation of Fig. i0.4c. A geometrically perfect piston might
have any of the attitudes shown in Fig. 10.5a to d. Also, geometric perfection of
piston and bore (soundness, parallelism, straightness) is impossible to achieve.
Imperfections contribute to uneven clearances and possible uneven axial pressure
distributions. Lack of parallelism, in either or both components, can generate the
clearance situations illustrated on Fig. 10.5e and f.
Fig. 10.5 includes the pressure distributions likely to be associated with various
piston-bore clearance ccnfigurations. Each indicates the net lateral force Fz which
will be generated by the pair of diametrically opposite pressure distributions shown.
For the parallel concentric and parallel eccentric cases (a and b) Fz is zero. For
the tilted piston cases (c and d) Fz acts to decentre the piston towards contact
with the bore. With the convergently-tapered clearance situation (e) Fz is finite
but acts to centre the piston in its bore, a highly desirable situation w h i c h holds
whether the piston is tilted or parallel in its bore. with the divergently-tapered
clearance (f) Fz acts to decentre the piston whether or not the piston is tilted.
The decentring lateral force necessary for hydraulic lock to develop will be present
if there is a diverging axial clearance due to piston tilt or component tapering.
210
Physically cc~%taminated fluid can erode and wear ccrnponents, degrading the
characteristics of control orifices and increasing internal leakage. It can also
cause direct physical malfunction of sensitive parts. For example, a sensitive
electrical control valve can silt up and jam if left stationary under pressure in
the presence of excessively cc~taminated oil.
materials of the particulate are also important. So too are the materials and
clearances of the operating system components.
Fig. 10.6 illustrates crudely some mechanisms which could lead to performance
aberrations in piston components. Hard particles of sizes of the order of the
diametral clearance, and of half the diametral clearance, would be highly undesirable
as they could enter the clearance and lodge in it, especially if pressure was reduced
after entry. Silting at the entry to a fine clearance can affect pressure distribution
in the clearance and so give rise to hydraulic lock in sensitive components.
a. Particle of Size Half c. Tilting Due to Lodged P a r t i c l e
Diameter Clearance
i Ps
b. P a r t i c l e of Size D i a m e t e r
Clearance
The person concerned with the dynamic response and performance of hydraulic
control systems should always be on the alert for aberrations in system performance
which may be attributable to the effects of fluid-borne contamination. Such
contamination will always be present. The idea is to keep it to levels at which the
system performance will not be affected, and which will allow long and reliable life.
212
10.7 CONCLUSION
Some phenomena which might or might not manifest themselves in working hydraulic
control systems have been discussed at an introductory level. There are other
parasitic phenomena. Thermal lock is one. Thermal lock occurs if working
components expand under the influence of temperature such that running clearances
become impaired. If manifested, parasitic phenomena will affect system performance
and response in a degradatory way. Most are not readily allowed for in dynamic
system models. If a system's actual response does not correlate well with a predicted
desired response, it could be due to one or more parasitic phenomena. For this
reason, the designer-analyst of hydraulic control systems should have familiarity
with them.
The phenomena discussed have been related mostly to piston configurations. They
may apply also to non-piston components.
10.8 REFERENCES
10.5 Rogers, K.J., On the Origins and Effects of Hydraulic Backlash in Tandem
Centre Valve Controlled Hydraulic Systems, Ph.D. thesis, Monash
University, 1976.
10.8 McCloy, D. and Martin, H., Control of Fluid Power, Longman, 1973.
10.11 Dransfield, P., Bruce, D.M., and Wadsworth, M., A General Approach to
Hydraulic Lock, Proc. I.Mech.E., V182, Part i, N27, 1967-68.
10.13 Tessman, R.K. and Fitch, E.C., Contaminant Induced Wear Debris for
Fluid Power Components, Proc. I.Mech.E. Technology Conf., Swansea,
U.K., 1978.
CONCLUSION
It is hoped that the Text has provided the reader with a reasonably deep
appreciation of the potential and the use of bond graphs as a key link in improving
his (read his or her) ability to predict the dynamic response of hydraulic control
systems. I encourage the reader with the desire but not the confidence to effect
such analyses, to try the bond graph approach. The reader already confident of his
ability to perform such analyses on proposed systems on the drawing board, may well
see bond graphs as just another modelling technique. I challenge such readers to
dig in sufficiently to be able to properly try the bond graph approach to simulation.
At the very least you will find it enlightening.
that bond graphs provide a key link in the confident development of dynamic
models,
that there are phenomena parasitic to the response of hydraulic systems and
that the system designer should be aware of these and of their possible
effects.
It is not claimed that the Text is exhaustive. There are bond graph techniques,
inherent and parasitic phenomena, optimization techniques, and related matters which
have been barely referred to or neglected. There are however sufficient material
and e x a ~ l e s to allow the hydraulic system designer-analyst to investigate many
systems, and References from which he can extend his knowledge when ready. The
Author claims no omnipotence in the area of dynamic performance of hydraulic control
systems - - simply a reasonable knowledge and experience of it which has been
significantly enhanced by his adoption of bond graphs.
APPENDIX 1
INTRCDUCTION
A considerable number of books provide information and techniques for the design
of hydraulic control systems. Some are included in the References of this Appendix.
The designer rarely needs to design a major component. His task is to conceive a
strategy appropriate to the particular application, to represent this strategy on a
circuit diagram, and to select the components for his system from a wide range of
commercial stock most of which will be available on-the-shelf. There may be delays
in supply if the shelf happens to be in another country~
It has been assumed that the readers of this Text are reasonably familiar with
hydraulic control systems and to some degree with their design. This Appendix will
only outline a basic approach for arriving at a hydraulic system capable of
performing a specific task.
POINTS TO BE CONSIDERED
Safety requirements
protection of operators and others;
protection of the system;
adjacent machinery or equipment.
Reliability
life of system;
likely maintenance requirements;
access to systems;
ease of maintenance;
availability of maintenance expertise and facilities;
operator ability.
Other factors
will operating noise be a problem?
does the client have any personal whims or requirements? Layout, preferred
maker of components, compatibility with existing system, etc.;
Size, space, mass, or layout limitations.
AN APPROACH
2. Propose a hydraulic circuit diagram for the system showing only major components
(actuators, control valves, pumps). This diagram constitutes a design strategy.
3. Appreciate any basic constraints placed on the system, e.g. maximum system
pressure, pump drive speed.
4. Specify the maximum strokes required of linear actuators and/or the maximum
speed required of hydraulic motors.
5. Assess the mauximum force or torque required of each actuator, allowing for
inertias, friction, and external forces or torques the drive will be expected
to overcome.
219
9. Determine the approximate pump discharge rate required on the bases of:
i0. Select commercial actuator(s) and control valve (s). Refine approximate pump
discharge rate requirement accordingly.
13. Decide if any intermediate valving or control circuits are required and select
suitable components.
14. Decide hydraulic line or hose sizes (length, internal diameter, strength).
16. Estimate the maximum temperature rise in the system oil due to working of the
system.
19. Allowing for pressure drops, internal leakages, and for the actual sizes and
ratings of the commercial items, refine all calculations.
20. Check that the system meets the action and performance requirements.
21. Estimate the cost of system components and of building the system.
22. Optimize the system and component selection on the bases of:
minimum cost;
adequacy of performance;
minimum number of components.
23. Finalize the design, component selection, and the estimate of cost.
Yeaple, F.D. (Ed.) , Hydraulic and Pneumatic Power and Control, McGraw-Hill, 1966.
Kibble, J.D., Basic Factors in the Design of Hydraulic Circuits, Proc. I.Mech.E.,
V.184, Part I, N.I, 1969-70.
SI CONVERSION FACTORS
T h e U S A is t h e m a j o r i n d u s t r i a l n a t i o n n o t y e t f u l l y u s i n g the SI s y s t e m . The
f o l l o w i n g t a b l e a l l o w s c c ~ v e r s i o n o f u n i t s co~muonly u s e d in t h a t c o u n t r y ' s f l u i d
p o w e r i n d u s t r i e s to b e c o n v e r t e d to SI.
Basic Common
Quantity Basic Conversion Alternatives
SI U n i t USA Unit
Basic Units
Derived Units
Area m2 in 2 in2/1550 = m 2 in 2 , 6 . 4 5 2 = c~n2
Volume m3 US g a l l o n US gal/264 = m 3 US g a l 3 . 7 8 5 =
Density kg/m 3 Ibm/ft 3 lhm/ft3,16.02 = kg/m 3
F o r c e (also newton (N) ibf ibf 4.448 = N
weight)
Torque N -m i b f " in (Ibf-in)/8.85 = N m
Rot.inertia kg - m 2 i b m - in 2 ( i b m ' i n 2 ) / 3 4 1 3 = kg-m;
Pressure pascal (Pa) i b f / i n 2 (psi) psi - 6895 = Pa )si/145 = M P a
Pa = N/m 2 psi/14.5 = bar
Energy(also joule ft " ibf (ft'ibf) 1.356 = J
work) (N - m)
Note: i0 n is s h o w n as En
i 0 -n is s h o w n as E-n
APPENDIX 3
There is a remarkably large number of English language papers concerned with the
dynamic response and performance of hydraulic control systems. They appear in a wide
variety of journals, with fairly regular contributions in the Proceedings of the
Institution of Mechanical Engineers, the Proceedings of BHRA's International Fluid
Power Symposia, the Bulletin of the Japan Society for Mechanical Engineers, and in a
number of the ASME journals. The present Bibliography is quite limited, choices for
it being based on whether a papex is directly referenced in the Text, whether it is
available in accessible literature, whether it is clearly relevant to the work of the
Text, and an attempt to provide the names of persons working in the area. Some
excellent older references are not included as they are referred to in later papers.
Books
Blackburn, J.F., Reethof, G., and Shearer, J.L., Fluid Power Control, M.I.T. Press,
1960.
Guillon, M., Hydraulic Servo Systems: Analysis and Design, Butterworths, 1969
(translation from French).
McCloy, D. and Martin, H.R., The Control of Fluid Power, Longman, 1973.
Wang, P.K.C., Mathematical Models for the Time Domain Design of Electrohydraulic
Servomechanisms, A.I.E.E. Trans. Applic. and Industry, VS0, 1961.
Wang, P.K.C. and Ma, J.T.S., Cavitation in Valve Controlled Hydraulic Actuators,
Trans. A.S.M.E. Jrl. App. Mech., V30, 1963.
McCloy, D. and Martin, H.R., Some Effects of Cavitation and Flow Forces in the
Eleetrohydraulic Servomechanism, Proc. I.Mech.E., V178, Pt.l, 1963-64.
Dransfield, P. and Bruce, D.M., Leakage Flow Rate Past Pistons of A Hydraulic
System Components, J. of Aircraft, A.I.A.A., V5, N2, 1968.
Lepp, R.M., A Study of the Time Optimized Hydraulic Servomechanism Operating Under
Cavitation Conditions, Ph.D. Thesis, U. of Saskatchewan, Canada, 1968.
Takenaka, T. and Urata, E., The Dynamic Characteristics of oil Hydraulic Control
Valves, Bull. J.S.M.E., VI2, 1969.
Burrows, C.R. and Webb, C.R., The Design of Fluid Power Systems Using Analogue
Computer, Trans. Inst. Meas. & Control, V2, 1969.
Haines, D.F. and Davies, R.M., Harmonic and Transient Behaviour of an Electro-
hydraulic Servomechanism Test Rig, Control, VI3, 1969.
Nikiforuk, P.N. and Tsai, S., Detailed Analysis of a Two-Stage Four-Way Electro-
hydraulic Control Valve, J. Mech. Eng. Sci., VII, N2, 1969.
Urata, E., The Ramp Response of a Loaded Hydraulic Servomechanism, Bull. J.S.M.E.,
VI3, 1970.
Green, W.L. and Crossley, T.R., An Analysis of the Control Mechanism Used in
Variable Delivery Hydraulic Pumps, Proc. I.Mech.E., V185, 1970-71.
Ulrich, H.J., Some Factors Influencing the Natural Frequency of Linear Hydraulic
Actuators, Intl. of Mach. Tool Design and Research, VII, 1971.
Armstrong, P.J. and McCloy, D., Optimization of a Hydraulic Servo Using On-Off
Controllers, Proc. 2nd Intl. Fluid Power Symp., BHRA, 1971.
Knight, G.C., McCallion, H., and Dudley, B.R., Connection Capacitance Effects in
Hydrostatic Transmission Systems and Their Prediction by Mathematical Model, Proc.
I.Mech.E., V86, 1972.
Keating, T and Martin, H.R., Mathematical Models for the Design of Hydraulic
Actuators, Trans. Inst. Soc. of Am., 1972.
Unruh, D., A Standard Format for Mathematical Models of Fluid Power Systems, Proc.
26th N.C.F.P., 1972.
Armstrong, P.J. and McCloy, D., Design of High Speed Hydraulic Position Servo-
mechanisms Using Optimization Techniques, Proc. I.Mech.E., V187, 1973.
Barnard, B.W., The Dynamic Response of a Hydraulic System, M.Eng.Sc. Thesis, Monash
University, Australia, 1973.
Martin, K.F., Design Approximations for Hydraulic Servos with Non-Symmetrical Oil
Volumes, Meast. & Control, V6, 1973.
Muto, T. and Hattori, T., A Study on the Unstable Phenomena in Hydraulic Driving
System and Direct Servosystem, Bull. J.S.M.E., VI7, 1974.
Helgestad, B.O., Foster, K., and Bannister, F.K., Pressure Transients in an Axial
Piston Hydraulic Pump, Proc. I.Mech.E., V188, 1974.
Dreymuller, J., Pilot-Operated and Directly Actuated Pressure Control with variable
Delivery Axial Piston Pumps, Proc. 4th Intl. Fluid Power Symp., BHRA, 1975.
Dransfield, P. and La Brooy, R., Hydraulic Systems with Precision Reflexes, Machine
Design, 20 May 1976.
Martin, K.F., Ong, C.S., Tee, B.L., and Little, E., Dynamic Analysis of an
Electrically-Operated Pressure Control Valve, Proc. 5th I.F.P.S., BHRA, 1978.
Svoboda, J., Analogue and Digital Modelling in the Design of a Hydraulic Vehicular
Transmission, Proc. I.Mech.E., V193, 1979 .
Akashi, H., Nakagawa, T., and O~umi, T., Analysis of Vibration in Hydraulic Drive
System, Hull. J.S.M.E., V22, 1979.
Stecki, J.S. and Cuthbertson, S.C., The Effects of Changeover to Fire Resistant
Fluids on Hydraulic Performance, Fire Resistant Hydraulics Jrl., VI, N2, pp 167-173,
FPRC, Oklahoma State University, 1981.
Backe, W., DSH Program System for Digital Simulation of Hydraulic Systems, Proc. 6th
Intl. Fluid Power Syrup., BHRA, 1981.
Bowns, D.E., Tomlinson, S.P., and Bonson, L.A., Progress Towards a General Purpose
Hydraulic System Simulation Language, Proc. 6th Intl. Fluid Power Syrup., BHRA, 1981.