Fluid power systems, in which working pressure
(pressure in pump output) is kept proportional to load, are called
hydraulic load-sensing systems. Such systems are mainly used in
mechanisms containing numerous drives to run with the purpose to save
energy. Hydraulic load-sensing systems are automatically regulating
systems with a number of components and several feedbacks. It is
difficult to find optimal solutions for their static, steady-state
motion and dynamics. A very precise parameter setting, especially for
resistances of hydraulic valve spools and for spring characteristics,
is required to make the system function.
As a result of current research, a simulation system is proposed that
enables one to perform computer experiments at the first stage of
design.
Theoretical fundamentals of the approach have been described in (Grossschmidt et al., 2009).
2. Hydraulic-MechanicalLoad-SensingSystem
At the first stage of the research the initial scheme of the hydraulic
load-sensing system of Bosch GmbH has been modelled and simulated.
Examples of modelling and simulation of this load-sensing system were
discussed in (Grossschmidt et al., 2006). In the initial scheme feeding
pressure of the controller was taken from the output of a variable
displacement hydraulic pump. But the pressure in the pump output
depends on the load and resistance of the system. Therefore the
behaviour of the load-sensing system depends on the output pressure of
the hydraulic pump. To get better results, it was necessary to improve
the model.
A modified scheme of load-sensing system is proposed (Fig. 1), in which
the controller has an independent constant pressure feeding. Feedback
pressures have been taken directly from the measuring valve with
pressure compensator RIDVW.
The scheme of hydraulic-mechanical controller is shown in Fig. 2 and the section of the valve block in Fig. 3.
Fig. 1: Scheme of the hydraulic–mechanical load-sensing system
Fig. 2: Scheme of a hydraulic-mechanical controller
Fig. 3: Section of the valve block
In Fig. 1, the variable
displacement axial piston pump is driven by an electric motor M.
Hydraulic- mechanical control of the pump volumetric flow is performed
by a control valve and positioning hydraulic cylinder of the swash
plate. The feeding chain of the hydraulic motor RVerbr contains tube
RL-zu, pressure compensator RIDW, measuring valve RVW, check valve,
meter-in throttle edge RSK-zu and connection elements. The output chain
of the hydraulic motor RVerb contains a meter-out throttle edge RSK-r,
and tube RL-ab. The device contains load-sensing pressure feedback.
The scheme of a hydraulic-mechanical controller (Fig. 2) contains a
control valve (effective area AV) with meter-in and meter-out throttle
edges, local hydraulic resistor (volumetric flow QDr), positioning
cylinder (effective area AZ), and swash plate with spring.
The main part of the valve block (Fig. 3) is a directional valve with
measuring throttle edge RVW, meter-in throttle edge RVS-zu and
meter-out throttle edge RSK-r of hydraulic motor. The valve block also
contains pres- sure compensator with throttle edge RIDW and check valve
causing pressure drop 0.5 bars.
3. Composing the Multi-pole Model
To build up the multi-pole model it is necessary to decompose the load-sensing system into subsystems and components.
The hydraulic-mechanical controller subsystem includes control valve,
meter-in throttle edge, meter-out throttle edge, local hydraulic
resistor and positioning cylinder with swash plate.
The valve block subsystem includes a measuring valve with pressure
compensator, check valve, meter-in throttle edge and meter-out throttle
edge of the hydraulic motor.
The hydraulic pump subsystem includes a variable displacement axial
piston pump together with electric motor and clutch. The hydraulic
motor subsystem includes a hydraulic motor together with meter-in and
meter-our throttle edges and clutch with drive mechanism.
The system includes tubes and tube connection elements as well.
Multi-pole models of all the
components must be composed. The multi-pole model of the load-sensing
system is built up using the components models. First, necessary
components must be connected through poles. Second, variables of
connection poles must be defined as inputs or outputs for every
component depending on required oriented causalities (Grossschmidt et
al., 1998, 2000 and 2003).
All the multi-pole models of the
load-sensing system components have been described as CoCoViLa
(Grigorenko et al., 2005 and 2006) classes together with their icons
and images. Relations (formulae) between variables of multi-pole models
are described as equations or as user-defined Java methods. In
relations state variables, inner iterations, logical conditions and
CoCoViLa subtasks are used for describing the behaviour of the
component.
The multi-pole model (Fig. 4) represents the scheme of the load-sensing system (Fig. 1).
The next step after configuring the hydraulic load- sensing system is parameters specification.
Firstly, the hydraulic motor and hydraulic pump parameters must be
chosen. An axial piston hydraulic motor with a working volume V = 40 x
10-06 m
3/rev has been taken. An axial-piston variable displacement pump with a maximum working volume V
max = 63 x 10-06 m
3/rev, a nominal rotation frequency n
nom = 1475 min-1 and a maximum displacement angle of the pump swash plate α
max = 0.3264 rad have been taken.
Secondly, the fluid and its properties must be chosen. The oil HLP46 was chosen.
Third, initial approximate values
of pressures and pressure drops for pump control must be set. Maximum
working pressure p
max = 250 bar, pressure drop in measuring valve Δp
IDW =
~5...7 bar, pressure drop in measuring
valve with
pressure compensator Δp
IDVW =~14...19 bar, pump control system feeding pressure p
0 = 60 bar and pump control pressure p
CP = ~4...19 bar were used.

Fig. 4: Multi-pole model of the hydraulic-mechanical load-sensing system
Fourth, maximum displacements of the valves must be set. Maximum displacement of the pump control valve y
V = 0.003 m, maximum displacement of the directional valve y
VW = 0.007 m and maximum displacement of the pressure compensator y
IDW = 0.007 m have been taken.
4. Mathematical Models
In this section mathematical
models implemented as multi-pole models of the components of the load
sensing system are presented. Values of parameters presented in the
section were obtained from performed simulations. The initial values
for iteration calculations, which have been found from stationary
calculations, correspond to displacement of the direction valve 0.0045
m and load moment 65 Nm.
4.1 VP - Displacement of the Control Valve
Static displacement of the control valve
y = ((p1+ p2) A - Ff) / c - y0 (1)
where A = π D
2 /4; F
f = F
f0 + k
f (p
1+ p
2) /2.
Difference of control valve velocity v
dv=(Δt/m)((p1+ p2)A-(y+y0) c – Ff sign(v)–hv) (2)
Difference of control valve displacement y
dy=Δtv (3)
Parameters: D = 0.008 m, F
f 0 = 0.3 N, k
f = 1E-7 m
2, c = 8377.5 N/m, y0 = 0.00846 m, m = 0.02 kg, h = 1200 Ns/m. Initial value: inity = 1.5819E-3 m.
4.2 RVP - Meter-in Throttle Edge of the Control Valve
Output pressure
p2=p1 -Q2 Q2 TR ρ/2, (4)
where TR = 1 / (μ
2A
12) + 1 / (μ
2A
22)
A
1 - longitudinal area of two identical slots,
A
2 - sectional area of two identical slots.
Scheme of the meter-in throttle edge of the control valve is shown in Fig. 5.

Fig. 5: Scheme of the meter-in throttle edge of the control valve
Geometrical relations: yS = y+l0,
h = yS (if yS <= hmax), hmax= R (1 – cos ((90 - β) π / 360)) l = Rα lmax = Rαmax,

A
2 = 2 (5e-08 +
dg)
Output volumetric flow Q
1 = Q
2. (5)
Parameters: D = 0.008 m, μ = 0.8, R = 1.1E-4 m, β = 0.30, l0 = 5.6E-4 m, g = 0.001 m.
Initial value: init p
2 = 1.1806E6 Pa.
4.3 RVT - Meter-out Throttle Edge of the Control Valve
Output volumetric flows

4.4 ResG, ResH, ResY – Local Hydraulic Resistors Output pressure of local hydraulic resistor ResG

4.5 IEH1-3 - Hydraulic Interface Element
The hydraulic interface element expresses:
• equality of output pressures, if pressure of one pole is given as input;
• the equation of continuity of volumetric flows

4.6 IEH11-2, IEH12-2, IEH2-12 – HydraulicInter- face Elements
The elements express equality of output pressures, if pressure of one
pole is given as input. The equations of continuity of volumetric flows
are as follows

4.7 ZV - Positioning Cylinder with Swash Plate
Static displacement of the swash plate spring

4.8 ME - Electric Motor
Angular velocity for steady state conditions


4.9 CJh – Clutch with Rotor of the Pump
For steady state conditions

4.10 PV - Variable Displacement Axial Piston Pump

4.11 TubeG, TubeH –Tubes form G and H
The tube model form G contains the
lumped inertia L, resistance R and volume elasticity C in sequence L
– R – C (Fig. 6a).

Fig. 6a: The structure of the four-pole model of form G

Fig. 6b: The structure of the four-pole model of form H
The tube model form H contains the
lumped vol- ume elasticity C, resistance R and inertia L in sequence C
– R – L (Fig. 6b).
For steady state conditions of TubeG

and volume elasticity C for achieving the same free oscillation
frequency as the model of tube with distributed parameters
For calculation of fluid properties and values Rl, Rt,
L, C of tubes a Java method is used. Kinematic viscosity, density and
compressibility factor depending on the pressure are calculated at each
time step.
4.12 Silencer with Interface Element
For damping oscillations in tubes,
silencers have been used. The silencer consists of a local hydraulic
resistor (ResG or ResH) with a deadlock tube (Tubeg or Tubeh).
Silencers are connected to the ends of tubes using interface elements
IEH11-2, IEH12-2, IEH2-12 and IEH1-3.
Parameters and initial values:

4.13RIDVW - Measuring Valve with Pressure Compensator
Measuring valve contains two identical slots. Slot profiles have been
designed to consist of 6 linear fragments of different shape in order
to achieve almost linear dependence of pressure drop on the
displacement of directional valve.
The pressure compensator contains four identical slots. Linear
dependence of the sum of pressure drops in the measuring valve and the
pressure compensator on the directional valve displacement must be
achieved. Slot profiles of the pressure compensator have been designed
to contain 15 linear fragments of different shape. It was quite
difficult to achieve linear depend- ence in the cases of minor
displacements of the directional valve (Fig. 7).

Fig. 7: Simulated pressure drop in measuring valve with pressure compensator depending on the displace- ment of the directional valve
However, in order to simplify the model of the load-sensing system,
model of the measuring valve with pressure compensator RIDVW has been
replaced by linear model RIDVWlin as follows.
For steady state conditions

4.14 RSKZ - Meter-in Throttle Edge of Hydraulic Motor and Check Valve
Meter-in throttle edge of hydraulic motor is as a cir- cular slot of
the directional valve. The pressure drop of the check valve is 5E4 Pa.
Output pressure

where TR = μ ((5e-08)+ π D y). Parameters: D = 0.018 m, μ = 0.7. Initial value:
init p
1 = 1.2008E7 Pa.
4.15 MH – Hydraulic Motor
Initial values: initp1 = 1.1892E7 Pa, initω = 107.98 rad/s, init Q1 = 7.0810E-4 m3/s.
4.16 RSKA – Meter-out Throttle Edge of Hydraulic Motor
The meter-out throttle edge of hydraulic motor is as a slot, which is formed by conical part of the directional valve.
4.17 CJhM – Clutch with Rotor of the Hydraulic Motor and Drive Mechanism
4.17 CJhM – Clutch with Rotor of the Hydraulic Motor and Drive Mechanism
Effiency coefficient of the hydraulic pump
eP = PPV/PME (64)
Effiency coefficient of the load-sensing system without pump
eHS = PMH/PPV (65)
Effiency coefficient of the load-sensing system
eG = PMH/PME (66)
5. Computing Process Organization
Using visual specifications of
described multi-pole models of fluid power system components one can
graphically compose models of various fluid power systems.
It is possible to solve a number of various computing tasks on each
fluid power system model evaluating some components as inputs and
computing some other components as outputs. When solving specific
simulation problem, the model has to be adjusted by evaluating
different parameters of the elements and adding sources to elements of
the model that describe disturbances of necessary shape and value.
When simulating steady state conditions behaviour of some parameter of
the fluid power system is calculated depending on the change of some
other parameter of the system.
When simulating dynamic behaviour, transient responses caused by
disturbances are calculated. Disturbances are considered as changes of
input parameters in certain points of the hydraulic system (pressures,
volumetric flows, load forces or moments, control signals, etc.).
Disturbances of different shapes are considered, such as constant,
step, impulse, sine, linear, etc. Calculations of dynamic responses are
performed in time. Time step length and number of steps are to be
specified. In all dynamic calculations the fourth-order classical
Runge–Kutta method has been used.
Steady state and dynamic computing processes are organized by
corresponding process classes. To follow the system behaviour in time,
the concept of state is invoked. State variables are introduced for
each functional element to characterize the element at the current
simulation time step. The simulation process starts from the initial
state and includes calculation of following state (nextstate) from
previous states (usually from oldstate and state). Final state
(finalstate) is computed as a result of simulation.
For solving a specific simulation task on the fluid power system model
the CoCoViLa program synthesizer is used for automatic construction of
the problem solving program.
We use a special technique for calculating variables in loop
dependences. We split the variable, assign it an initial approximate
value and try iteratively recompute the variable. Recomputing algorithm
can be synthesized by the CoCoViLa program synthesizer. Splitting the
variable, assigning an approximate value and asking CoCoViLa to find
algorithm for recomputing must be described in the multi-pole model of
the fluid power system component.
Variables with index “e” are obtained as results of iterations.
6. Modelling and Simulation Stages
The following steps must be performed for simulation both steady state conditions and dynamic transient responses.
First, all the models of components must be tested separately. For
every component the simulation task must be composed and input signals
must be chosen. Behaviour of the component must be simulated. Initial
approximate components parameters values (e.g. etc.) must be refined as
a result of simulation. Stiffness of springs, geometry of valves
working slots, etc. must be refined for steady state conditions.
Dynamic resistances, time constants, damping constants, etc. must be
refined for dynamics.
Second, the separately tested
components models must be connected into sub-systems, tested in
behaviour and adjusted.
Third, the whole load-sensing system must be built up, simulation tasks must be solved and the parameters must be adjusted.
7. Simulating Steady State Conditions
7.1 SimulatingSubsystems
7.1.1Hydraulic-mechanical Controller
The hydraulic-mechanical controller includes constant pressure feeding
that enables one to make the hydraulic pump control independent of
pressure in pump output. The simulation task description is shown in
Fig. 8.
Fig. 8: Simulation task description of hydraulic- mechanical controller
Notations: VP
- Displacement of the control valve; constant Source - Constant input
values of p1, p2 and F; static Source - Input pressure p1 of VP; p1p2 -
Calculation of pressure difference p1-p2; RVP - Meter-in throttle edge
of the control valve; IEH1-3 - Hydraulic interface element; ZV -
Positioning cylinder with swash plate; ResY - Local hydraulic resistor;
RVT - Meter-out throttle edge of the control valve.
Simulation results are shown in Fig. 9.
The hydraulic-mechanical controller has been simulated in case the
pressure difference acting on the control valve changes in the interval
14.1...19.1 bars.
Fig. 9: Simulated graphs of the hydraulic-mechanical controller (dependences on the pressure difference acting to the control valve)
Control valve (1) moves 0 ... 0.003 m. Volumetric flow to the control valve (2) changes in the interval 2.28E-5...4.25E-5 m
3/s.
Pump control pressure (3) interval, required for pump volumetric flow
regulation from max to min is 4...19 bars. Angular position of the pump
swash plate (4) decreases from maximum 0.3264 to 0 rad.
7.1.2Hydraulic Motor Subsystem
The simulation task description for the steady state conditions of the hydraulic motor subsystem is shown in Fig. 10.
Fig. 10: Simulation task of hydraulic motor subsystem
Notations:
RSKZ - Meter-in throttle edge for hydraulic motor; MH - Hydraulic
motor; RSKA - Meter-out throttle edge for hydraulic motor; TubeH -
Tubes; constant Source - Constant input for pressure p2 at the right
end of the subsystem and for load moment of the hydraulic motor M;
static Source - Input volumetric flow Q1 to RSKZ and input displacement
y of the directional valve.
Simulation results are shown in Fig. 11. Dependen- ces on the
displacement of the directional valve in the interval 0.0025...0.00685
m have been calculated for the case the load moment of the hydraulic
motor has constant value 65 Nm. Pressure at hydraulic motor inlet (1)
increases from 116.5 to 121.0 bar. Pressure at meter-in throttle edge
inlet (2) increases from 117 to 123.5 bar. Angular velocity of the
hydraulic motor (3) linearly increases from 0 to 225 rad/s.
Fig. 11: Simulated dependences of hydraulic motor subsys- tem on the displacement of the direction valve (load moment 65 Nm)
7.2 Simulating Steady State Conditions of the Load-Sensing System
A simulation task description for calculating steady
state conditions of the load-sensing system is shown in Fig. 12.
When simulating the steady state conditions of the hydraulic-mechanical
load-sensing system, two types of tasks are considered. First, we
consider dependences on the displacement of the directional valve in
the interval 0.0025...0.00685 m if the load moment of the hydraulic
motor is of constant value 65 Nm. Second, we consider dependences on
the load moment of the hydraulic motor in the interval 0...130 Nm, if
the displacement of the directional valve has a constant value 0.0045 m.
The task descriptions for simulations of both types differ only by
parameters of “static Sources” that specify simulation
input values.
Simulated dependences of the task of the first type are shown in Fig. 13 and Fig. 14.
Fig. 12: Simulation task description for the steady state conditions of the load-sensing system
Notations:
VP - Displacement of the pump control valve; RVP - Meter-in throttle
edge of the pump control valve; IEH1-3, IEH2-2m - Hydraulic interface
elements; ZV - Positioning cylinder with swash plate; ResY - Hydraulic
resistor; RVT - Meter-out throttle edge of the pump control valve; PV -
Variable displacement pump; ME - Electric motor; WG - Efficiency
coefficients calculator; TubeH - Tubes; RIDVWlin - Measuring valve with
pressure compensator; RSKZ - Meter-in throttle edge of the hydraulic
motor; MH - Hydraulic motor; RSKA - Meter-out throttle edge of the
hydraulic motor; constant Source – Constant input values of pump
control system feeding pressure p1 and outlet pressure p2; static
Source - Input values of the hydraulic motor load moment M and the
directional valve displacement y; static Process – Simulation
manager.
Fig. 13: Simulated dependences on the displacement of the directional valve (load moment 65 Nm)
In the ideal case the
dependences of pump control pressure and position angle of the swash
plate (Fig. 13) would be linear. As the shape of the graphs depends
mostly on the passage areas of the throttle edges of the valves it is
very difficult to achieve exact linearity. Displacement of the control
valve (1) linearly drops from 0.003 to 0 m. Pump control pressure (2)
drops from 19 to 4 bars and the volumetric flow of the pump (3)
increases from 0 to 0.00148 m3/s.
In Fig. 14 the efficiency coefficient of the whole load-sensing system
(1) increases from 0.63 to 0.67 if displacement of the directional
valve is higher than 0.003 m.

Fig. 14: Simulated dependences on the displacement of the directional valve (load moment 65 Nm)
In cases of less displacement, the
efficiency coefficient falls rapidly. Angular velocity of the hydraulic
motor (2) increases from 0 to 224 rad/s.
Resulting graphs of simulating task of the second type are shown in Fig. 15.

Fig. 15: Simulated dependences on the load moment of the hydraulic motor (displacement of the directional valve 0.0045 m)
The efficiency coefficient of the
pump (1) is 0.903 if the load moment of the hydraulic motor is zero.
The efficiency coefficient of the pump is maximal 0.944 if the load
moment is 45 Nm. If the load moment increases the efficiency
coefficient slightly falls to 0.927, because of increasing the pump
volumetric losses. The efficiency coefficient of the load-sensing
system (2) increases in the interval 0...0.70 in the full interval of
change of the hydraulic motor load moment.
8. Simulating Dynamics of the Load-Sensing System
Initially simulations of dynamic
transient responses were performed on the model of the load sensing
system (see Fig. 4) (Grossschmidt et al., 2008). But the results of
simulations showed that the system turned to be unstable. The model of
the dynamics of the load- sensing system had to be improved. First,
hydraulic damping resistors were inserted to pressure inputs of pump
control valve VP and to input of the positioning cylinder with swash
plate ZV. Second, silencers were inserted to the ends of connecting
tubes.
A simulation task description for calculating transient responses of
the improved load-sensing system is shown in Fig. 16. Below three
examples are considered.
In all the examples simulations have been performed for 1.3 s, using
time step Δt = 10 μs. Results have been calculated for 130 000
points.
In the first example, a constant displacement 0.0045 m of the
directional valve was taken. A step change +10 Nm (during 0.01 s) was
applied to the load moment of the drive mechanism 65 Nm.
In Fig. 17 and Fig. 18 graphs of simulating the dynamic responses are shown.
Fig. 16: Simulation task description for the dynamic transient responses of the load-sensing system
Notations:
VP - Displacement of the pump control valve; RVP - Meter-in throttle
edge of the pump control valve; ZV - Positioning cylinder with swash
plate; ResG, ResH, ResY –Local hydraulic resistors; RVT -
Meter-out throttle edge of the pump control valve; PV - Variable
displacement pump; ME - Electric motor; RIDVWlin - Measuring valve with
pressure compensator; RSKZ - Meter-in throttle edge of the hydraulic
motor; MH - Hydraulic motor; RSKA - Meter-out throttle edge of the
hydraulic motor; CJh – clutch; CJhM – clutch with drive
mechanism; IEH1-3, IEH11-2, IEH12-2, IEH2-12 - Hydraulic interface
elements; TubeH, TubeG, Tubeg, Tubeh - Tubes; constant Source –
Constant inputs; dynamic Source - Input values of the hydraulic motor
load moment M and the directional valve displacement y; dynamic Process
– Simulation manager.
Fig. 17: Graphs considering hydraulic pump (first example)
In Fig. 17 all the transient responses damp in 1.3 s. Angular velocity
of the hydraulic pump (1) stabilizes at the initial level 153.15 rad/s.
Volumetric flow of the hydraulic pump (2) drops from 0.00075 to 0.00069m
3/s. Output pressure of the hydraulic pump (3) stabilizes at the level 1.54E7 Pa.

Fig. 18:Graphs considering the system outputs (first example)
In Fig. 18, the angular velocity of
the drive mechanism (1) oscillates and stabilizes at the value 102
rad/s in 1.3 s. The graph of pressure at the outlet of the hydraulic
motor (2) has an almost similar shape.
In the second example, a constant load moment of the drive mechanism,
65 Nm, was taken. A step change -0.0015 m (during 0.01 s) was applied
to the displacement of the directional valve, 0.0045 m.
In Fig. 19 to Fig. 22 graphs of simulating the dynamic responses are shown.

Fig. 19: Graphs considering hydraulic pump (second example)
In Fig. 19 angular velocity of the hydraulic pump (1) increases from
153.2 to 154.0 rad/s. Volumetric flow of the hydraulic pump (2) drops
from 0.00075 to 0.00027 m
3/s. Output pressure of the hydraulic pump (3) oscillates and damps to the initial level.

Fig. 20: Graphs considering volumetric flows in hydraulic pump control device (second example)
In Fig. 20 volumetric flow of the meter-in throttle edge of the control valve (1) increases from 3.94E-5 to 4.225E-5 m3/s, and volumetric flow of the meter-out throttle edge (2) drops from 1.2E-5 to 0.9E-5 m3/s. Volumetric flow of the constant hydraulic resistor (3) increases from 2.75E-5 to 3.3E-5 m3/s. Volumetric flow of the positioning cylinder (4) oscillates and damps at the initial level 0 m3/s.

Fig. 21: Graphs considering hydraulic-mechanical controller (second example)
In Fig. 21 pump regulating pressure
(1) increases from 1.1E6 to 1.6E6 Pa, position angle of the pump swash
plate (2) drops from 0.17 rad to 0.06 rad and displacement of the pump
control valve (3) increases from 0.0014 to 0.0026 m.

Fig. 22: Graphs considering the system outputs (second example)
In Fig. 22, the angular velocity of the drive mechanism (1) drops
asymptotically from the initial value 115 rad/s to 40 rad/s. The graph
of pressure at the outlet of the hydraulic motor (2) carries higher
frequency oscillations in the initial phase of the transient response.
These oscillations are caused by absence of silencer at the tube
between hydraulic motor and meter-out throttle edge of the hydraulic
motor. There is no need for silencer as it does not affect the
behaviour of the whole load sensing system.
In the third example, both a step change +10 Nm to the load moment of
the drive mechanism 65 Nm and step change -0.0015 m to the displacement
of the directional valve 0.0045 m were applied simultaneously (during
0.01 s).
In Fig. 2 to 25 graphs of simulating the dynamic responses are shown.

Fig. 23: Graphs considering hydraulic-mechanical controller (third example)
In Fig. 23, the hydraulic pump
regulating pressure (1) increases from 1.1E6 to 1.6E6 Pa, the position
angle of the pump swash plate (2) drops from 0.17 rad to 0.06 rad and
the displacement of the pump control valve (3) increases from 0.0015 to
0.0026 m.

Fig. 24: Graphs considering hydraulic pump (third example)
In Fig. 24, the angular velocity of
the pump (1) increases from 153.2 to 153.9 rad/s and the volumetric
flow of the hydraulic pump (2) drops from 0.00075 to 0.00027 m3/s in
1.3 s. Output pressure of the hydraulic pump (3) oscillates and damps
at the level 1.55E7 Pa.

Fig. 25: Graphs considering the system outputs (third ex- ample)
In Fig. 25, the angular velocity of
the drive mechanism (1) drops from the initial value 110 to 40 rad/s.
The graph of pressure at the outlet of the hydraulic motor (2) has an
interesting shape, because no silencer is used at the tube between the
hydraulic motor and the meter-out throttle edge of the hydraulic motor.
9. Size and Complexity
The simulating task for calculating transient responses of the load-sensing system considered in the simulation contains:
- 43 classes, including 27 component classes;
- more than 1200 parameters;
- 24 variables that have to be iterated during the computations;
- 112 links between system components.
The automatically constructed Java code for solving the simulation task
for calculating transient responses of the load-sensing system contains
6967 lines of Java source code and includes 4 algorithms for solving
different subtasks.
10. Recommendations
The following recommendations concerning the load-sensing system can be pointed out as a result of the performed simulations.
- The hydraulic-mechanical controller must have independent feeding with constant pressure.
- Feedback pressures have been taken directly from the measuring valve with pressure compensator.
- A very precise parameter setting, especially for spring
characteristics and for slot geometry of hydraulic valves, is required
to make the system function.
- All the system components must react with optimal delay and with optimal damping.
- The clutches must have elasticity of optimal values.
- Diameters of the connecting tubes must be optimal and lengths must be different to avoid resonance.
- To achieve dynamic stability of the system hydraulic resistors
for damping of pump control valve and positioning cylinder must be
invoked. Silencers must be invoked into connecting tubes as well.
11. Using Simulation in Fluid Power System Design
Simulations must be performed at
the first stage of design. The results of simulations can be used as a
starting point when building trial versions of real load- sensing fluid
power system components and subsystems.
As a result of tests of trial
versions of the fluid power system components and subsystems
mathematical models and parameters must be refined.
Simulations must be performed once again to prove correctness of used
solutions. Finally the designed system must be experimentally refined
and adjusted. In this way we can achieve a good performance of the
designed fluid power system.
12. Conclusions
As a result of the current research, a modelling and
simulation system is proposed that enables one user- friendly perform
computer experiments in design stage of large and complicated fluid
power systems.
An example of multi-pole modelling and simulation of steady-state
conditions and dynamics of a hydraulic-mechanical load-sensing system
has been considered.
The visual programming environment CoCoViLa supporting declarative
programming in a high-level language and automatic program synthesis
has been used.
The main features of the approach proposed are as follows:
• Mathematical models of the functional elements
are composed as multi-pole models taking into account signal
propagation in both directions.
• Used multi-pole models can have various oriented causalities.
• The mathematical model of the fluid-power
system includes models of functional elements, carries the full
information about connections of in- put/output variables, expresses
the considered mathematical oriented causalities and thus guarantees
the completeness of the model.
• Simulation is performed step by step, starting
from simulation of components and moving to more complicated subsystems.
• Calculations are performed separately for each
multi-pole model. Iteration methods are used in cases of loop
dependencies that may appear between component models when they are
connected together into more complicated ones.
Nomenclature
As a result of the
current research, a modelling and simulation system is proposed that
enables one user- friendly perform computer experiments in design stage
of large and complicated fluid power systems.



Acknowledgement
This research was supported by the Estonian Sci- ence Foundation (Grant No. 7091).
References
Grigorenko, P.,
Saabas, A. and
Tyugu,
E. 2005. COCOVILA – Compiler-Compiler for Visual Languages. Proc.
of the 5th Workshop on Language Descriptions, Tools and Applications,
(Edinburgh, 2005), v. 141, n. 4 of Electron. Notes in Theor. Comput.
Sci., pp. 137-142. Elsevier.
Grigorenko, P. and
Tyugu,
E. 2006. Deep Semantics of Visual Languages. In: E. Tyugu, T. Yamaguchi
(eds.) Knowledge-Based Software Engineering. Frontiers in Artificial
Intelligence and Applications, Vol. 140. IOS Press, pp. 83 - 95.
Grossschmidt, G. and
Vanaveski,
J. 1998. Causality of Mathematical Models of Technical Systems. Proc.
of the 12th European Simulation Multiconference ESM ́98, as part of the
50th anniversary celebrations of the University of Manchester,
Simulation - Past, Present and Future, June 16-19, Manchester, United
Kingdom, pp. 191-195.
Grossschmidt, G. and
Vanaveski,
J. and Harf, M. 2000. Simulation of Hydraulic Chains using Multi-pole
Models in the NUT Programming Environment. Proc. of the 14th European
Simulation Multiconference ESM’2000, Simulation and Modelling:
Enables for a better quality of life, May 23-26, Ghent, Belgium, pp.
709-713.
Grossschmidt, G. and
Harf,
M. 2003. Multi-pole Modelling and Simulation of an Electro-hydraulic
Servo-system with a Non-linear Regulator in NUT Environment.
Fünftes Deutsch-Polnisches Seminar “Innovation und
Fortschritt in der Fluidtechnik”, 18. -19. September, Warszawa,
pp. 146-161.
Grossschmidt, G. and
Harf,
M. and Sallaste, T. 2006. Modelling and Simulation of Fluid Power
Systems in Object-Oriented Programming Environment. 8th Biennial ASME
Conference on Engineering Systems Design and Analysis ESDA2006, July
4-7, Turin, Italy, ASME, paper ESDA 2006-95387, pp. 1- 10.
Grossschmidt, G. and
Harf,
M. 2008. Multi-pole modelling and simulation of a hydraulic mechanical
load-sensing system using the CoCoViLa programming environment.
Proceedings of 6th International Fluid Power Conference "Fluid Power in
Motion", April 1-2, Dresden, Germany, Vol. 2, pp. 553-568.
Grossschmidt, G. and
Harf,
M. 2008. Modelling and simulation of fluid power systems in an
intelligent programming environment. Proceedings of the 6th Industrial
Simulation Conference, June 9-11, Lyon, France, pp. 224-230.
Grossschmidt, G. and
Harf,
M. 2009. COCO-SIM Object-Oriented Multi-Pole Modelling and Simulation
Environment for Fluid Power Systems, Part 1: Fundamentals.
International Journal of Fluid Power, August Issue.