Software for Fluid Power Technology


From Editor

The purpose of the Software Review section of the Journal is to present information to the reader about engineering software, including simulation programs, to highlight their specific features and their "fitness to purpose" in the unique field of fluid power and motion control. It is, of course, impossible to establish evaluation criteria matching the needs of all readers, therefore readers should not look for absolute ratings but more or less "fuzzy" opinions of the reviewer. A software program is like a wrench, just a tool to solve problems. It is good to solve some problems and not so good for others and this depends on both the nature of the problem and the users' attitude - and generally when we review software we do not know either. A software tool can be highly specialised and great for a some applications but not so well suited for others, on the other hand another software tool can be more flexible and generally applicable but without outstanding features. It is impossible, and even misleading, to say which one is better. What we hope to accomplish is to give the reader information necessary to take his/her own decision.

CoCO-SIM Object-Oriented Multi-Pole Modelling and Simulation Environment for Fluid Power Systems Part 2: Modelling and simulation of Hydraulic-Mechanical Load-Sensing System

Contact Information


   
Vendor:  Tallinn University of Technology, Institute of Machinery and Institute of Cybernetics
Authors of the system: Gunnar Grossschmidt1, Mait Harf
Address: 1) Ehitajate tee 5, 19086 Tallinn, Estonia
               2) Akadeemia tee 21, 12618 Tallinn, Estonia
Telephone Number: 1) +372 519 89785, 2) +372 620 4215 
Fax Number: 1) +372 620 3264, 2) +372 620 4151
Email: gunnar.grossschmidt@ttu.ee, mait@cs.ioc.ee
Web: http://www.cs.ioc.ee/cocovila/
Platform: Java based, platform independent (tested on Windows, Linux, Mac OS)

1. Introduction


        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.

An output coupled power-split transmission

Fig. 1: Scheme of the hydraulic–mechanical load-sensing system

An output coupled power-split transmission

Fig. 2: Scheme of a hydraulic-mechanical controller

An output coupled power-split transmission

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 m3/rev has been taken. An axial-piston variable displacement pump with a maximum working volume Vmax = 63 x 10-06 m3/rev, a nominal rotation frequency nnom = 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 pmax = 250 bar, pressure drop in measuring valve ΔpIDW = ~5...7 bar, pressure drop in measuring    valve    with    pressure    compensator ΔpIDVW =~14...19 bar, pump control system feeding pressure p0 = 60 bar and pump control pressure pCP = ~4...19 bar were used.

An output coupled power-split transmission

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 yV = 0.003 m, maximum displacement of the directional valve yVW = 0.007 m and maximum displacement of the pressure compensator yIDW = 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 = π D2 /4;     F   = Ff0    + k   (p1+ p2) /2.

Difference of control valve velocity v

dv=(Δt/m)((p1+ p2)A-(y+y0) c – F sign(v)–hv)    (2)

Difference of control valve displacement y
dy=Δtv (3)


Parameters: D = 0.008 m, Ff 0 = 0.3 N, k= 1E-7 m2, 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 / (μ2A12) + 1 / (μ2A22)
A1 - longitudinal area of two identical slots,
A2 - sectional area of two identical slots.
Scheme of the meter-in throttle edge of the control valve is shown in Fig. 5.

An output coupled power-split transmission

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,

An output coupled power-split transmission

A2 = 2 (5e-08 + dg)

Output volumetric flow Q1 = Q2.    (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 p2 = 1.1806E6 Pa.

4.3 RVT - Meter-out Throttle Edge of the Control Valve

Output volumetric flows

An output coupled power-split transmission


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


An output coupled power-split transmission

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

An output coupled power-split transmission


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

An output coupled power-split transmission


4.7    ZV - Positioning Cylinder with Swash Plate


Static displacement of the swash plate spring

An output coupled power-split transmission

4.8    ME - Electric Motor


Angular velocity for steady state conditions

An output coupled power-split transmission

An output coupled power-split transmission

4.9 CJh – Clutch with Rotor of the Pump


For steady state conditions

An output coupled power-split transmission

4.10 PV - Variable Displacement Axial Piston Pump



An output coupled power-split transmission

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).

An output coupled power-split transmission

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

An output coupled power-split transmission

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

An output coupled power-split transmission



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).

An output coupled power-split transmission

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  p1 = 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 m3/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.00069m3/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 m3/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:
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.

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.



BACK



©2009 Fluid Power Net All rights reserved. Disclaimer