lib hydraulic amesim software tutorial
DESCRIPTION
amesim hydraulic model tutorial examples and guide imagine R11 multi physics softwareTRANSCRIPT
-
Hydraulic LibraryVersion 4.2 - September 2004
-
Copyright IMAGINE S.A. 1995-2004
AMESim is the registered trademark of IMAGINE S.A.
AMESet is the registered trademark of IMAGINE S.A.
Ara
AIn
M
NcaNC
P
Uli
Wte
TSan
X
AcoDAMS is a registered United States trademark of Mechanical Dynamics, Incorpo-ted.
DAMS/Solver and ADAMS/View are trademarks of Mechanical Dynamics, corporated.
ATLAB and SIMULINK are registered trademarks of the Math Works, Inc.
etscape and Netscape Navigator are registered trademarks of Netscape Communi-tions Corporation in the United States and other countries. Netscapes logos and etscape product and service names are also trademarks of Netscape Communications orporation, which may be registered in other countries.
ostScript is a trademark of Adobe Systems Inc.
NIX is a registered trademark in the United States and other countries exclusively censed by X / Open Company Ltd.
indows, Windows NT, Windows 2000, Windows XP and Visual C++ are regis-red trademarks of the Microsoft Corporation.
he GNU Compiler Collection (GCC) is a product of the Free Software Foundation. ee the GNU General Public License terms and conditions for copying, distribution d modification in the license file.
windows is a trademark of the Massachusetts Institute of Technology.
ll other product names are trademarks or registered trademarks of their respective mpanies.
-
Hydraulic Library 4.2User Manual
Chapter 1: Tutorial examples
1.1 In
T
Hsomfr
M
Um
S
U
H
Uel
H
Tty1
troduction
he AMESim Hydraulic library consists of:
A collection of commonly used hydraulic components such as pumps, motors, orifices, etc. including special valves.
Submodels of pipes and hoses.
Sources of pressure and flow rate.
Sensors of pressure and flow rate.
A collection of fluid properties.
ydraulic systems in isolation are completely useless! It is necessary to do mething with the fluid and also to control the process. This means that the library ust be compatible with other AMESim libraries. The following libraries are equently used with the Hydraulic library:
echanical library
sed in fluid power application when hydraulic power is translated into echanical power.
ignal, Control and Observer library
sed to control the hydraulic system.
ydraulic component design library
sed to build specialist components from very basic hydraulic and mechanical ements.
ydraulic resistance library
his is a collection of submodels of bends, tee-junctions, elbows etc. It is used pically in low pressure applications such as cooling and lubrication systems.
-
Chapter 1Tutorial examples
2
CreledoC
1.2 E
O
Inthcofr
T
Note: It is possible to use more than one fluid in the Hydraulic library. This is important because you can model combined cooling and lubrication systems of a library.
The hydraulic library assumes a uniform temperature throughout hapter 1 of the manual consists of a collection of tutorial examples. We strongly commend that you do these tutorial examples. They assume you have a basic vel of experience using AMESim. As an absolute minimum you should have ne the examples in Chapter 3 of the AMESim manual and the first example of
hapter 5 which describes how to do a batch run.
xample 1: A simple hydraulic system
bjectives
Construct a very simple hydraulic system
Introduce the simplest pipe/hose submodels
Interpret the results with a special reference to air release and cavitationFigure 1.1: A very simple hydraulic system
this exercise you will construct the system shown in Figure 1.1. This is perhaps e simplest possible meaningful hydraulic system. It is built partly from mponents from the Hydraulic category (which are normally blue) and partly
om the Mechanical category.
he hydraulic part is built up from standard symbols used for hydraulic systems.
the system. If thermal effects are considered to be important, the Thermal Hydraulic and Thermal Hydraulic Component Design libraries should be used.
There are models of cavitation and air release in the hydraulic library. Note also there is a special two-phase flow library. A typical application for this is air conditioning systems.
Pump
Primemover
Tank
Relief valveRotary Load
Motor
-
Hydraulic Library 4.2User Manual
The prime mover supplies power to the pump, which draws hydraulic fluid from a tank. This fluid is supplied under pressure to a hydraulic motor, which drives a rotary load. A relief valve opens when the pressure reaches a certain value. The output from the motor and the relief valve returns to the tank. The diagram shows three tanks but it is quite likely that a single tank is employed.
Tspbeiccom
S
Tasch3
he first category contains general hydraulic components. The seconds contains ecial valves.The hydraulic components used in the model you will build can all found in the first of these Hydraulic categories. If you click on this category on, you will have the dialog box shown in Figure 1.2. First look at the mponents available in this library. Display the title of some components by oving the pointer over the icons.
Figure 1.2: The components in the first hydraulic category.
tep 1: Use File New... to produce the following dialog box.
here are two categories in the Hydraulic library. These have blue the standard color. If you do not have this categories displayed, eck the path list in the Options menu.
-
Chapter 1Tutorial examples
4
Figure 1.3: The hydraulic starter system.
S
1.
2.
3.
4.
SOco
Yyotep 2: Construct the rest of the system and assign submodels
Construct the system with the components as shown in Figure 1.1.
Save it as hydraulic1.
Go to Submodel mode.Notice that the drop, the prime mover, the node and the pipes are not of normal appearance because they do not have submodels associated with them. The easiest way to proceed is:
Click on the Premier submodel button which is situated in the horizontal menu bar.
elect the hydraulic starter circuit libhydr.amt and then click on K. A new system with a fluid properties icon in the top left rner of the sketch will be created.
ou could also have clicked on the New icon in the tool bar but if u do this you will have to add the fluid properties icon yourself.
-
Hydraulic Library 4.2User Manual
Figure 1.4: The line submodels.
Yasthre
Awap
Rof
Inthinbyno
S
1.
2.
5.
6.5
ou get something like Figure 1.4. It is possible that your system may have HL000 sociated with one of the other line runs. These minor variations are dependent on e order in which you constructed the lines. They will not influence the simulation sults.
n important feature to notice is that a line run has a special submodel (HL000) hich is not a direct connection. To emphasize this point the line run has a special pearance.
emember the submodel DIRECT does nothing at all. It is as if the ports at the end the line were connected directly together.
contrast, HL000 computes the net flow into the pipe and uses this to determine e time derivative of pressure. If the net flow into the pipe is positive, pressure creases with time. If it is negative, it decreases with time. The pressure created HL000 is conveyed to the relief valve inlet. The motor inlet is conveyed by the de and submodel DIRECT.
tep 3: Set parameters
Change to Parameter mode.
Set the following parameters and leave the others at their default values:
Press the mouse right-button.
Select Show line labels in the label menu.
Submodel Title Value
HL000 pipe length [m] 4
RL00 coefficient of viscous friction [Nm/(rev/min)]0.02
-
Chapter 1Tutorial examples
6
Figure 1.5: Setting the line submodel HL000 parameters.
3.
4.To display the parameters of a line submodel click the left mouse button with the pointer on or near the appropriate line run.
Part of the dialog box for HL000 is shown in Figure 1.5. The compressibility of the oil and the expansion of the pipe or hose with pressure are taken into account together with the pipe volume. HL000 normally requires the bulk modulus of the hydraulic fluid and pipe wall thickness together with the Youngs modulus of the wall material. From these values an effective bulk modulus of the combined fluid and pipe walls can be calculated. The effective bulk modulus of a hose is normally very much less than that of a rigid steel pipe.
Click on the fluid icon FP04 in the sketch.
A new dialog box as shown in Figure 1.6 is displayed. This shows you the properties of the hydraulic fluid. Currently they are at their default values and the absolute viscosity, bulk modulus, air/gas content and temperature are given in common units.
-
Hydraulic Library 4.2User Manual
Figure 1.6: Parameter for fluid properties submodel FP04.
5.
S
1.
2.
3.7
Click on OK.
tep 4: Run a simulation
Go to Run mode and do a simulation run.
The default values in the Run Parameters dialog box are suitable for this example.
Click on the Start run button.
Click on the pump component to produce the dialog box shown in Figure 1.7.
Some variables such as a pressure have no direction associated with them. A (gauge) pressure of -0.1 bar indicates that the pressure is below atmospheric. In contrast other variables, such as flow rate, do have a direction associated with them. A flow rate of -6 L/min indicates that the flow is in the opposite direction to some agreed standard direction.
Note that the first item in the list is an enumeration integer parameter. A collection of properties of varying complexity are available but for this exercise elementary is satisfactory.
-
Chapter 1Tutorial examples
8
Figure 1.7: The Variable List for PU001.
4.
5.Note that you can use the Replay facility to give you a global picture of the results. Figure 1.8 also shows the flow rates in L/min at a time of 10 seconds.
Figure 1.8: Flow rates displayed in replay.
To plot a variable associated with a line submodel, click on or near the corresponding line run.
Plot pressure at port 1 for HL000.
-
Hydraulic Library 4.2User Manual
Figure 1.9: The pressure in the hydraulic pipe.
Cavitation
W
Tcalo9
Notice how the pressure goes up to just over the relief valve setting (150 bar). During this time the load speeds up rapidly and actually 'over-speeds'. At this point the motor is demanding more hydraulic flow than the pump can supply. The result is that the pressure must drop and the relief valve closes. The pressure continues to drop and falls below zero bar gauge. However, pressure is not like voltage or force. We cannot have a pressure of -100 bar. The absolute zero of pressure is about -1.013 bar gauge. It is time to introduce two terms.
and air release
hen pressure falls to very low levels, two things can happen:
Air previously dissolved in the fluid begins to form air bubbles.
The pressure reaches the saturated vapor pressure of the liquid and bubbles of vapor appear.
hese phenomena are known as air release and cavitation respectively. They can use serious damage. Using the Zoom facility, the graph gives a better view of the wer pressure values:
Figure 1.10: Low pressure in the hydraulic pipe.
-
Chapter 1Tutorial examples
10
All AMESim submodels have hydraulic pressure in bar gauge. The low pressure shown in Figure 1.10 : Low pressure in the hydraulic pipe. is caused by the load speed exceeding its steady state or equilibrium value and it is a highly undesirable behavior as it can result in damage to the real system.
In reality the starting values we have given for the pipe pressure and load speed are nobesy
NT
ATflH
TcatoA
Tspru
1.3 Ep
O
Int very realistic and the prime mover would start from rest or some valve would used to regulate the flow to the motor. However, hydrostatic transmission stems like this often do suffer badly from cavitation and air release problems.
ote that all AMESim submodels have hydraulic volumetric flow rate in L/min. here are two possible interpretations of this flow rate:
The flow rate is measured at the local current hydraulic pressure, or
The flow rate is measured at a reference pressure.
MESim adopts the second alternative with a reference pressure of 0 bar gauge. his means that the volumetric flow rate is always directly proportional to the mass ow rate. In most situations the difference between the two flow rates is negligible. owever, there are three situations when there is a significant difference:
1. There is a very large air content; the pressure drops below the satu-ration pressure for air in the liquid and air bubbles are formed in the liquid.
2. The pressure drops to the level of the saturated vapor pressure of the liquid and cavities of vapor form.
3. Extremely high variations in pressure occur such as in certain types of fuel injection systems.
he first situation is called air release and the second cavitation. If there is vitation or significant air release at the inlet to a pump, the flow rate according the first definition will not be reduced but with the approach adopted by MESim it will be significantly reduced.
he properties of hydraulic fluids vary a great deal. Modeling them is a very ecialist process and the model can be extremely simple or highly complex. The n times are greatly influenced by this level of complexity.
xample 2: Using more complex hydraulic roperties
bjectives:
Use more complex models of fluid properties.
See how air content changes the performance of the system.
the Hydraulic category two special components can modify the fluid properties:
-
Hydraulic Library 4.2User Manual
Figure 1.11: The two fluid properties icons.
Tan
T
1.
Inwth
2.
A collect of
Special model usedto ensure compatibility
simple andcomplex fluidproperties.11
his is an example of a component without ports. We cannot connect this icon to y other.
here are two important thing about FP04.
It has an integer parameter index of hydraulic fluid that is in the range 0 to 100 inclusive. This arrangement means that it is possible to have more than one fluid in an AMESim system.
elementary This is the default and features a constant liquid bulk modulus with absolute viscosity. The treatment of fluid properties under air release and cavitation is done.
simplest This has a constant absolute viscosity. The bulk modulus is constant above the gas saturation pressure and is 1/1000 of this value below the gas saturation pressure. This model is very old but is still used by some AMESim users.It is likely to give the fastest runs.
advanced This gives you access to some cavitation parameters not accessible in the elementary properties.
advanced using tables This is like the advanced option but you install tables of data to give variation of bulk modulus and absolute viscosity with pressure and temperature.
Robert Bosch adiabatic diesel These properties are provided by Robert Bosch GmbH and comprise a number of common types of diesel fuel.
AMESim always use this fluid properties icon. It is associated ith one submodel: FP04. The other icon and its submodel is ere only for backward compatibility.
The characteristics of the fluid properties are determined by its parameters. One of which is an enumeration integer parameter. There are 5 possibilities:
between 4.0 models andearlier. Do not use this one.
-
Chapter 1Tutorial examples
12
Using one of the special fluids
Step 1: Use the Advanced fluid properties.
1. Return to the first example of this manual, add another fluid properties icon.
2. Use Premier submodel and go to Parameter mode. Your sketch should look
3.like this.
Figure 1.12: The sketch with two instances of FP04.
Look at parameters of FP04-2. Change the enumeration integer parameter to advanced. The Change Parameters list should now look like this:
Figure 1.13: The advanced fluid properties.
Change the index of hydraulic fluid in FP04-2 to 1. This is a number in the range 0 to 100. If you look at the other hydraulic components in the system you will find they have index 0 and hence they will still use the fluid properties of FP04-1. We could go into every hydraulic component using this second fluid and set the parameter index of hydraulic fluid to 1. This would be extremely tedious with a big system and there is always the possibility of missing one.
-
Hydraulic Library 4.2User Manual
Step 2: Set all fluid indices to the same value of 1
The best way to do this is to use the common parameters facility.
1. Use Edit Select all.
All the system components will be selected, unselect FP04-1 holding the
2.
3.
S
S
1.
2.
3.
4.13
SHIFT key and clicking on it.
Use Parameters Common parameters.
Figure 1.14 shows the Common parameters dialog box. This is a list of common parameters for selected objects. They occur at least twice. Since there are 3 hydraulic tanks and they all have pressures of 0 bar, this value is displayed. There are a number of submodels that have a parameter index of hydraulic fluid. In FP02 the index of hydraulic fluid is set to 1 whereas in other submodels its value is 0. The value is displayed as ???. Similarly the prime mover and rotary load both have a parameter (strictly speaking variable) with title shaft speed. Since the two values are different, ??? is displayed.
Figure 1.14: Different values for common parameters
Set the parameter index of hydraulic fluid to 1. This will change all the parameters in the system except FP01 (remember we selected Select all except FP01).
tep 3: Run a simulation and plot some variables
You will probably find the results very much the same as in example 1.
tep 4: Organize a batch run to vary the air content
In Parameter mode use Parameters Batch parameters.
Drag and drop the air gas content from FP04-2 to the Batch control parameter setup dialog box.
Set up the batch parameters as in Figure 1.15 so that the air content goes from 0% to 10% in steps of 2%.
Specify a batch run in the Run parameters dialog box and initiate the run.
-
Chapter 1Tutorial examples
14
Figure 1.15: Setting up a batch run varying air content.
5.
6.
7.Plot several graphs of the batch run to compare results with various air contents.
Figure 1.16: Pressure in pipe.
By zooming on the curve in regions where the pressure is below 0 bar you will probably find some, but not a remarkable variation in the results.
Change the saturation pressure in FP04-2 to 400 bar.
Repeat the batch run and update your plot.
-
Hydraulic Library 4.2User Manual
Figure 1.17: Pressure in pipe with saturation pressure 400 bar.
Tchex
NItpooi
Aththsaa
S(Fsa
Hinwhi15
he variation between the runs is now very pronounced. The dynamic aracteristics of the system is completely transformed. A few words of planation are necessary.
ormally the air content of a hydraulic oil is well below 1%, and 0.1% is typical. is normally considered good practice to keep the value as low as reasonably ssible. However, in a few application, such as lubrication oil in gearboxes, the l and air are well mixed up and 2.5% is typical and up to about 10% is possible.
reasonable quantity of air, given time, will completely or partially dissolve in e hydraulic fluid. The lowest pressure at which all the air is dissolved is called e saturation pressure. For very slow systems all the air is dissolved above the turation pressure and partially dissolved below this pressure. Henrys law gives reasonable approximation for the fraction of air that is dissolved in equilibrium.
ome systems are slow enough to stay very close to this equilibrium position igure 1.16). Often classic fluid power systems behave like this. The original turation pressure is better for the current example.
owever, it does take time for the air to dissolve and this time will not be available fast acting systems. Fuel injection systems are a good example of this. Hence ith such systems it may be appropriate to set the saturation pressure artificially gh to allow for significant quantities of air to be undissolved at all pressures.
-
Chapter 1Tutorial examples
16
1.4 Example 3: Using more complex line submodels
Objectives:
Tdethsu
S
1. Use more complex line submodels.
Understand the need for a variety of line submodels.
To understand the importance of setting an appropriate line submodel.
he system for this example is the same as for example 2 (Figure 1.12). We will scribe the modification of the system to use more complex line submodels and e experiments performed. Finally we give a little of the theory behind the bmodels.
tep 1: Change submodels
All the submodels in the current system were selected automatically. We will change some of them manually.
Go to Submodel mode.
You will now change some line submodels.
Before continuing note the following points:
The corners in the pipe runs are not physical but diagrammatic.
There are three hydraulic pipes and they meet at a point which physically will be a tee-junction.
It is necessary to have a large number of hydraulic pipe submodels.
In the present system three submodels are set: DIRECT, DIRECT and HL000.
Figure 1.18: The current line submodels.
This tee-junction in the sketch is described as a 3 port node and it has the submodel H3NODE1. This models the junction as a common pressure with flow rates that give conservation of mass.
-
Hydraulic Library 4.2User Manual
None of these line submodels take friction into account. We will suppose that the relief valve is close to the node but the pump and the motor are at such distances from the node that the pressure drop along the pipes cannot be ignored. We need to select new pipe submodels that take friction into account for the pipe runs from the pump to the node and from the node to the motor.
2.
3.
4.
S
1.17
Click on the line run attached to the pump and select HL03 in the Submodel list.
Figure 1.19: The hydraulic line submodels available.
Why did we not choose a more complex submodel that also included inertia? We answer this question later in this exercise.
For the line from the node to the motor, select the submodel HL01.
For the line between the node and the relief valve, the submodel DIRECT is already selected and this is exactly what we want.
tep 2: Set parameters and run a simulation
Go to Parameter mode and set parameters for HL01 and HL03 so that both pipe lengths are 5 m and pipe diameters are 10 mm.
This can be done one at a time. However, we can do it another way. Press the Shift key on click on the HL03 and HL01 line runs so that they are selected. Use Parameters Common parameters. Figure 1.20 shows the Common parameters dialog box.
Note the brief description of each line submodel. In these descriptions C stands for compressibility, R for resistance (pipe friction) and I for inertia (fluid momentum). HL000 which we used before takes into account compressibility only. HL03 takes into account compressibility and friction. It is modeled like two hydraulic compressible volumes with a resistance between them.
-
Chapter 1Tutorial examples
18
Figure 1.20: The common parameters of the two line submodels.
2.
3.
4.
NsifaNote that ??? indicates that different values are set in the line submodels. Set the index of hydraulic fluid to 1, diameter of pipe to 10 and pipe length to 5.
In FP04-2 reset the saturation pressure (for dissolved air/gas) to 0 bar.
Run a simulation with the default run parameters.
Plot the two pressures in HL03.
Figure 1.21: Pressures at the ends of pipe joining pump to node.
ote that there is a large pressure drop along the line. This could be regarded as a zing problem but in addition it would be bad practice to site the relief valve so r from the high pressure point.
-
Hydraulic Library 4.2User Manual
Step 3: We now investigate other line submodels.
1. Return to Sketch mode and Copy-Paste part of the system as shown:Figure 1.22: Part of the system is duplicated.
2.
3.19
In Submodel mode change the lower two line submodels as follows:Figure 1.23: New line submodels.
This system will enable you to make direct comparisons between results.
Go to Run mode and do a simulation. Plot the pressure at the pump outlet (pressure at port 2).
Figure 1.24: Pressure at pump outlet.
-
Chapter 1Tutorial examples
20
We note that the curves are virtually the same. (Try zooming.) There is absolutely no advantage to using HL04 and HL06 instead of HL01 and HL03. If we separated the two systems and ran then independently we would find run times for the more complex submodels were higher.
4. Change the communication interval in the Run Parameters dialog box to
0.001s and rerun the simulation.
If you have a look at the Warnings/Errors tab of the Simulation run dialog box, you will find that some checks are performed by the line submodels (see Figure 1.25). A similar message is issued for HL03.
Figure 1.25: Messages under the Warning tab.
It is suggested that:
HL01 should be replaced by HL04 and
HL03 should be replaced by HL06.
In other words with this communication interval the lower subsystem is better than the upper. If you replot the pressures at the pump outlets, there are clearly differences. This is what happens if you zoom.
Figure 1.26: Zoomed pressures at pump outlet.
The violent (and unrealistic) start up has created this oscillation in the pressure of about 56 Hz. It is damped out by 0.1 seconds. Why did we get no warning message in the previous run? The answer is that a lot of checks are applied to your submodel choices when the run starts. These take into account the fluid properties, the pipe dimensions and the communication interval. When the communication interval was 0.1 seconds, it would have been impossible to see these oscillations and hence no warnings are issued.
-
Hydraulic Library 4.2User Manual
Some very simple arithmetic gives two important points:
1. If you want to see f Hz you need a communication interval no bigger than about 1/(10f) seconds.
2. A communication interval of x seconds enables you to see frequen-cies down to about 1/(10x) Hz. Thus with 0.1 seconds you can see
1.5 E
O21
down to about 1 Hz.
In the current example we are probably not interested below 1 Hz and 0.1 seconds, HL01 and HL03 are very suitable. The following note is vitally important.
See also Chapter 4:Selecting submodels for Hydraulic Lines.
xample 4: Valves with duty cycles
Figure 1.27: Hydraulic system with servovalve.
bjectives
Introduce valves controlled by duty cycles.
Use the plot manager to plot flow rates against differential pressure.
Note: It is a common mistake in modeling hydraulic systems to always use line submodels of high complexity. The correct procedure is to use the simplest line submodels that will achieve the modeling objective. Be aware of the frequencies you are interested in and the frequencies you can see with the current communication interval.
If you make a bad choice, you may multiply the CPU time by 10 or 100 and force the integrator to compute high frequency phenomena which are of no interest to you and will be invisible with the communication interval set! The messages under the Warning tab are very helpful. Read them.
-
Chapter 1Tutorial examples
22
Step 1: Build the system and set parameters
1. Build the system shown in Figure 1.27 and save it as servovalve.
Note that you have a directional valve that you will use to change the direction of rotation of the load. You will need to use two new components:
2.
3. A 3-position 4-port direction valve found in the first Hydraulic category, and
A duty cycle component found in the Signal, Control and Observers category.
When the new system sketch is complete, use Premier submodel to get the simplest combination of submodels.
Set the parameter values for the duty cycle submodel UD00 as follows:
This gives a signal as follows:
Figure 1.28: The duty cycle controlling the valve.
Title Value
duration of stage 1 [s] 1
output at start of stage 2 [null] 40
output at end of stage 2 [null] 40
duration of stage 2 [s] 3
output at start of stage 3 [null] -40
output at end of stage 3 [null] -40
duration of stage 3 [s] 3
Note: If you do not change the parameters, the valve will not open. The motor and load will not move at all! For simplicity leave the other submodel components with their default settings.
-
Hydraulic Library 4.2User Manual
4. Select the directional valve (Figure 1.27).
You will not change any parameters but an understanding of the parameters of SV00 will help you to set those of UD00. The spool has some state variables which are the first two items in the list. In Parameter mode, their values are the initial values of these state variables. The spool position is a fraction and so is 23
a dimensionless quantity in the range -1 to 1.
Figure 1.29: Parameters of the servo valve
The next 12 items determine the hydraulic flow characteristics of the valve covering the 4 possible flow paths. When the valve is in one extreme position with fractional spool position +1, P is connected to A and T to B. In the other extreme position the spool position is -1 and the connections are A to T and B to P.
When the spool position is 0, there is no flow. To define the flow characteristics of the valve in the extreme positions a flow rate pressure drop pair is used. The default values of these are 1 L/min and 1 bar. These values can normally be found in a manufacturer's catalogue. The parameter critical flow number (laminar -> turbulent) is less important and can be left at its default value. You can find the details for any submodel if you click on the Help button. For SV00, this produces the dialog box shown in Figure 1.30.
-
Chapter 1Tutorial examples
24
Figure 1.30: Help for submodel SV00.
5.
6.
7.
Sre
1.The rating of the valve rated current is set to 40 mA. This means that an input signal of 40 units will produce a fraction spool position of 1. As the spool moves it behaves like a second order system. You can specify the natural frequency and damping ratio.
Enable Discontinuities Printout in the Run Parameters dialog box.
Run a simulation with default run parameters.
Select the relief valve component and plot the three quantities:
Flow at relief valve outlet [L/min]
Pressure at relief valve inlet [bar]
Pressure at relief valve outlet [bar]
on the same plot.
tep 2: Plot the flow rate against the differential pressure for the lief valve.
This is a very common requirement for a 2-port valve and it involves use of the plot manager.
Start the Plot manager by clicking on the button in the tool bar on your plot.
The Plot manager is displayed as in Figure 1.31.
-
Hydraulic Library 4.2User Manual
Figure 1.31: The Plot Manager.
2.
3.
5.
4.25
In the right window we have the three quantities we requested and time. We must create a new variable which will be the pressure drop across the valve.
Click on Add item and construct the new variable.
Click on Add curve and drag and drop the flow rate onto Y and the pressure drop onto X.
Finally click on OK to see the plot.
Curve 1 to 3 are no longer required so select each in turn and then click on Remove curve so that only Curve 4 remains.
-
Chapter 1Tutorial examples
26
Figure 1.32: The relief valve flow rate pressure drop characteristics.
1.6 Ea
Oxample 5: Position control for a hydraulic ctuator
bjectives:
Use a simple proportional control system to achieve a prescribed cycle in a hydraulic system.
Show the consequences of using an unequal area actuator.
Show saturation in a servo-valve.
Study stability and instability in the control system.
-
Hydraulic Library 4.2User Manual
Figure 1.33: The position control system.
T(osecywtoan
S
1.
2.
Position sensorHydraulic actuator or jack
227
he system sketch for this exercise is shown in Figure 1.33. The hydraulic actuator r jack) moves a load and there is control using position feedback. The position nsor is used to convert the actuator displacement to a signal. A position duty cle is specified by a duty cycle submodel. The duty cycle position is compared ith the position indicated by the sensor to produce an error. The error is subjected a gain and the signal transferred to the servo-valve. A further duty cycle supplies external force to the actuator via the position transducer.
tep 1: Build the system and set parameters
Build the new system and save it as actuator.
The position sensor is found in the category labeled Mechanical. A signal port is used to pass the displacement into the feedback loop.
Use the Premier submodel button to select the simplest possible submodels combinations.
Forceconversionblock
Servo valve
Force dutycycle
1 3 4
Hydraulicaccumulator
-
Chapter 1Tutorial examples
28
3. Set parameters for the submodels using the suggested values in the following table:
N
Submodel Number on sketch if
any
Title Valueote:
HJ000 piston diameter [mm] 30
diameter of rod [mm] 20
length of stroke [m] 1
PU001 pump displacement [cc/rev] 35
UD00 1 duration of stage 1 [s] 1
output at end of stage 2 [null] 0.8
duration of stage 2 [s] 3
output at start of stage 3 [null] 0.8
output at end of stage 3 [null] 0.8
duration of stage 3 [s] 1
output at start of stage 4 [null] 0.8
output at end of stage 4 [null] 0.2
duration of stage 4 [s] 3
output at start of stage 5 [null] 0.2
output at end of stage 5 [null] 0.2
UD00 2 output at the end of stage 1 [null] 1000
output at end of stage 1 [null] 1000
SV00 valve natural frequency [Hz] 50
valve damping ratio [null] 1
valve rated current [mA] 200
DT000 gain for signal output [1/m] 10
GA00 3 value of gain [null] 10
GA00 4 value of gain [null] 250
-
Hydraulic Library 4.2User Manual
The parameters of HJ000 give a very unequal area actuator and the plots can demonstrate the consequences of this.
The external force to the right of the actuator is a constant value of 1000 N.
The gain in the displacement sensor converts the jack position that is in
4.
S
1.29
the range 0 to 1 m to a signal in the range 0 to 10. The gain for the submodel GA00 attached to the duty cycle submodel, is also 10. By this means, the duty cycle will directly represent the actuator displacement in m.
When you set the parameters for HJ000, click on the External variables button to call up the dialog box shown in Figure 1.34.
Figure 1.34: External variables of HJ000
This indicates that a positive velocity means the rod is moving to the right. The greater the displacement, the further it is to the right. In the current case, a zero displacement and velocity means that the rod and piston are stationary and the piston is at the extreme left end of the jack.
The meaning of the sign of the acceleration and external force should be clear. A positive external force opposes the other variables i.e. makes a negative con-tribution to the acceleration. Hence it is trying to reduce the velocity and dis-placement.
Remove the dialog box by clicking on Close.
tep 2: Run simulation and plot results
Run a simulation setting a final time of 12 s and a communication interval of 0.05 s.
-
Chapter 1Tutorial examples
30
2. Plot the following graphs:
Actuator displacement and duty cycle output on the same graph.
Flow rate at the two actuator ports on the same graph.
Flow rate at pump outlet and flow rate at relief valve outlet on the same graph.
3.
4.
5.
6. Fractional spool position.
Figure 1.35: The required and the actual displacement.
The first plot (Figure 1.35) gives an idea of how closely the actual performance matches the required duty cycle.
Plot the output from the summing junction (strictly speaking a differencing junction) that gives you the position error in m.
Try changing the gain attached to the servo valve, the servo valve natural frequency and damping ratio.
Include a high gain value that makes the system unstable.
Try introducing a dead band, up to about 10%.
-
Hydraulic Library 4.2User Manual
Figure 1.36: Pump and relief valve flow rates.
7.
8.31
A typical plot for the flow rates from the pump and relief valve outlets is shown in Figure 1.36. If you had chosen the pump inlet flow rate instead of the pump outlet flow rate, negative values would have appeared on the graph. This is easily explained if you click on the External variables button of the Variable List dialog box. For both ports of the pump a positive flow rate indicates flow out of the pump. It follows that the flow rate at the pump inlet must be negative.
Plot the two flow rates in the actuator HJ000.
For the this submodel, flow rate is an input on both flow ports. This means a positive flow rate indicates flow into the component. Figure 1.37 shows typical results. Note how different the magnitudes of the flow rates are due to the unequal areas.
Figure 1.37: Hydraulic actuator flow rates.
Plot the valve spool fractional displacement.
This gives an idea of how close to saturation the valve is during the duty cycle. If a value of +1 or -1 is reached, the valve is saturated.
-
Chapter 1Tutorial examples
32
1.7 Example 6: Simple design exercise for a hydraulic suspension
Objectives:
D
Tdath
S
1.o a simple initial design study for a hydraulic suspension using:
Analytical analysis.
AMESim standard runs.
Batch runs.
Linear analysis.
he system is shown in Figure 1.38. The hydraulic jack with the two orifices is the mper and the accumulator is the spring. It is proposed to use this suspension on e cab of a truck. The load on each suspension strut is 250kg.
Figure 1.38: A simplified hydraulic suspension.
tep 1: Build the system and run a simulation
Build the system using Premier submodel.
Much sizing can be done by simple calculations but simulation can be a great help in rapidly confirming the calculations and adding dynamics to the steady state values. The two ports of the jack are interconnected and in equilibrium. The pressures above and below the jack piston will be the same. Using a force balance in the equilibrium position in terms of the piston area Apist and rod area Arod
It follows that
PApist P Apist Arod( ) 250xg=
PArod 250g=
-
Hydraulic Library 4.2User Manual
From this if we want an operating pressure of about 70 bar the diameter of the rod must be about 22.3 mm. We will use a rod diameter of 20 mm and a piston diameter of 40 mm.
2. Set the parameters of the following table.
3.
P
Submodel Parameter title Value33
Run a dynamic simulation for 10 s.
Figure 1.39: Pressure and displacement plots.
Figure 1.39 shows the system pressure and the displacement.
roblem 1: The starting values are poor.
HJ000
rod displacement [m] 0.15
piston diameter [mm] 40
diameter of rod [mm] 20
angle rod makes with horizontal [degree] 90
total mass being moved 250
-
Chapter 1Tutorial examples
34
Problem 2: The accumulator spring with its precharge pressure of 100 bar is taking no part in this simulation. The only spring involved at the moment is the hydraulic fluid.
Solution to problem 1:
1. In Parameter mode select Parameters Set final values.
2.
3.
S
Tacbe
Tto
wce
S
1.
1.This will give reasonable starting values for state variables. You will find that the piston has dropped slightly from the mid-position.
Reset the following parameters:
Run a simulation again and check that the system is in equilibrium with the rod in mid-position.
olution to problem 2:
he two parameters we can vary are the precharge pressure and volume of the cumulator. For the accumulator to work as a spring, the precharge pressure must lower than the equilibrium pressure.
he volume of fluid in the jack varies according to the piston position. This is due the rod. The difference between the minimum and maximum oil volume is
hich is 0.1 L. The accumulator volume should be a bit bigger than this but rtainly not 10 L.
tep 1: Investigate the spring rate
Set the following values
Do a run and verify that these values do not disturb the equilibrium.
The values should have changed the spring rate but not the equilibrium position. We need now to investigate the spring rate.
Submodel Parameter title Value
HJ000rod displacement [m] 0.15
rod velocity [m/s] 0
Submodel Parameter title Value
HA001gas precharge pressure [bar] 10
accumulator volume [L] 0.5
Arod x stroke
-
Hydraulic Library 4.2User Manual
2. Set the following values
3.
4.
Tcocoen
Toft
Submodel Parameter title Value
output at start of stage 1 [null] 0
output at end of stage 1 [null] 250035
Do a run for 120 s.
Plot graphs of:
rod displacement (HJ000)
pressure at port 1 (HJ000)
against
external force on rod (HJ000)
Figure 1.40: Displacement against force.
he force value of 2500 N pushes down on the suspension with a value rresponding to the weight of the car. The force of -2500 effectively takes the mplete weight off the suspension. The slow evolution of the force duty cycle sures that the system is very close to equilibrium at all times.
he plot of displacement against force (Figure 1.40) shows the non-linear nature the spring. It also shows that the suspension does not bottom out but it does op out.
UDOOduration of stage 1 [s] 40
output at start of stage 2 [null] 2500
output at end of stage 2 [null] -2500
duration of stage 2 [s] 80
-
Chapter 1Tutorial examples
36
Figure 1.41: Force against pressure.
Fba
Wdovo
Hwar
Sth
1.
2.
3.igure 1.41 shows that maximum pressure is 160 bar and the minimum is about 40 r which occurs when the suspension tops out.
e could continue by doing further analytical calculations. Alternatively we could batch runs varying the accumulation pre-charge pressure and accumulator lume and the interested reader could try this.
owever, we will end the exercise by considering the damping of the suspension hich is mainly provided by the two orifices. For simplicity we will assume they e of the same characteristics.
tep 2: Setup a batch run varying the diameters of the orifices with e vehicle subject to a step change in force
Select Parameters Global parameters.
Set up the global parameter shown.
Set the following parameters for BOTH orifices:
Submodel Parameter title Value
OR0001 for pressure drop/flow rate pair 2 for orifice diameter
2
equivalent orifice diameter [mm] DIAM
-
Hydraulic Library 4.2User Manual
4. Set up a duty cycle to give a step increase in force:
5.
6.
7.
T
Submodel Parameter title Value
output at start of stage 1 [null] 0
output at end of stage 1 [null] 037
Select Parameters Batch parameters.
Drag and drop the global parameter into the Batch parameters dialog box and set the following values for a batch run:
Figure 1.42: Batch parameters
Perform a batch run for 10 s and plot the displacement of the piston.
Figure 1.43: Batch run results for rod displacement
he batch run will use orifice diameters of 1 to 6 mm in steps of 0.5mm. Zooming
UD00duration of stage 1 [s] 1
output at start of stage 2 [null] 500
output at end of stage 2 [null] 500
duration of stage 2 [s] 9
-
Chapter 1Tutorial examples
38
in on the plot it becomes clear that 3 mm gives a reasonable degree of damping.
8. Remove the step from UD000 so that there is a constant force of 0 N.
9. Insert a linearization time at 10 s.
10. Repeat the batch run and look at the damping ratio for the oscillatory frequency.
L2.dibegras
Wooking at the eigenvalues selecting the .jac0.1 to .jac0.11 files we see that below 5 mm the system is very highly damped. However, the results for the 1 mm ameter give an oscillatory frequency of about 25 Hz which is curious but could investigated with tools such as modal shapes. For diameters of 2.5 mm and eater there is an oscillatory frequency of about 1.23 Hz and the damping ratio is follows:
e can see the evolution of these eigenvalues in a root locus plot.
Diameter of orifice [mm]
Damping ratio
2.5 0.533
3 0.308
3.5 0.194
4 0.130
4.5 0.091
5 0.067
5.5 0.050
6 0.039
-
Hydraulic Library 4.2User Manual
Figure 1.44: Root locus plot.
Ase39
more refined search between 2.0 and 3.0 mm would be a good idea but 2.5 mm ems reasonable.
-
Chapter 1Tutorial examples
40
-
Hydraulic Library 4.2User Manual
Chapter 2: Theory of fluid properties
W
2.1 D
T
D[kte
T
se
T
w
an
Tth41
e will concentrate mainly on three fluid properties in this chapter:
The density which leads to mass and hence to hydraulic inertia effects.
The viscosity which leads to the hydraulic friction effects.
The compressibility and thus the bulk modulus which leads to the hydraulic system stiffness. Notice that the compressibility effect can be modified by air release, cavitation phenomena and by expansion of a pipe, hose or chamber containing the hydraulic fluid.
ensity and compressibility coefficient
he density is the mass of a substance per unit volume:
ensity has dimensions of [M/L3] and is expressed in kilograms per cubic meter g/m3]. As mentioned previously the density is a function of the pressure and the mperature:
his function can be approximated by the first three terms of a Taylor
ries:
his can also be expressed as:
ith
d
his equation is the linearized state equation for a liquid. Using the definition of e density, the two coefficients and B can also be expressed as:
MV-----=
P T nature of fluid, ,( )=
P P T T+,+( ) P------
TP T------
P
T+ +=
1 PB------- T+
=
B P------
T=
1---
T------
P
=
B V PV------
Tand 1V
--- VT------
P= =
-
Chapter 2Theory of fluid properties
42
B is known as the isothermal bulk modulus or for simplicity the bulk modulus and is known as the cubical expansion coefficient. Since fluid density varies with the applied pressure, this implies that a given mass of fluid submitted to a pressure change changes its volume. This phenomenon leads to the definition of the compressibility coefficient :
wfo
it
U
Mm
TRcoabexthsh
Tai
Entrapped
Waili
1 V here is expressed in units Pa -1 (or m2/N). Considering the relation r a closed hydraulic circuit the mass is constant, and hence:
follows that
sing the definition of the compressibility coefficient we obtain:
ore usually we use the bulk modulus B also known as the volumetric elasticity odulus:
he relation between and B implies mass conservation. This relation must be IGOROUSLY RESPECTED in the calculations. In the modeling and simulation ntext of fluid energy systems, disregarding the relation between and B leads to normal evolutions of pressure in the closed circuit submitted to compression and pansion cycles. This phenomenon is strongly accentuated if aeration occurs in e circuit (when dissolved air in the fluid reappears in the form of bubbles). We all approach this point by examining the phenomena of aeration and cavitation.
he air can also have adverse consequences on a fluid compressibility. In liquid r can be present in two forms: entrapped and dissolved.
air
hen the return pipe is not submersed in the tank the liquid jet can entrain some r bubbles in the tank. Another phenomenon that affects the quantity of air in quid is the leakage.
V--- P------ T=
V M=
d V( ) 0 and Vd dV+ 0= =
d------
dVV
-------=
1---
P-------------=
B P-------------=
-
Hydraulic Library 4.2User Manual
Figure 2.45: Liquid leakage
Dissolved a
43
This air stays in the liquid as cavities and can modify the fluid compressibility. In this context we talk about effective bulk modulus. Figure 2.46 shows the bulk modulus of a diesel fuel at 40 C with 0, 0.01, 0.1, 1, 10% air. The plot is obtained using the system shown. The model of the diesel fuel properties is based on accurate ex-perimental measurements and are designed for use with injection system which are very fast acting. For this reason air is assumed to be entrained rather than dissolved.
Figure 2.46: Air is entrained
ir
Air can also be dissolved in a liquid. A certain amount of air molecule can be part of the liquid. In this case the dissolved air does not significantly change the fluid properties.
-
Chapter 2Theory of fluid properties
44
2.2 Air release and cavitation
Air can be dissolved or entrained in liquids and it is possible for air to change from one of these two forms to the other depending on the conditions to which the fluid is
S(ucadifrdiofnotore
IfpAusco
T
Tnu
Saturation
Vapor pr subjected.
uppose the fluid is in equilibrium with a certain percentage of dissolved gas sually air: nitrogen and oxygen). Lowering the pressure above a critical value lled the saturation pressure induces aeration. This is the process where the ssolved gas forms air bubbles in the liquid until all the dissolved gases or air are ee. The exact point where all the dissolved gas has come out of solution is fficult to pin-point because it depends on the chemical composition and behavior the gas. This is a non-symmetrical dynamic process: the growing process does t have the same dynamics as when air bubbles disappear. In consequence the tal amount of bubbles created when the pressure drops may or may not be dissolved in the liquid when it rises again.
the pressure is dropped further and above another critical value called the vapor ressure, the fluid itself starts to vaporize. It corresponds to a liquid phase change. t some point only fluid vapor and gas exist. In liquid systems the term cavitation ually refers to the formation and collapse of cavities in the liquid even if cavities ntain air or liquid vapor.
o summarize with a sketch what we have introduced see above:
Figure 2.47: Air release and cavitation
he development of a cavity is now recognized as being associated with a cleation center such as microscopic gas particles, wear or wall asperities. When
Time
cavitation Vapor (vaporized liquid)
pressure
Liquid pressure Re dissolution (total or partial)
essure
Air bubbles appearance
Air bubbles
-
Hydraulic Library 4.2User Manual
the liquid is subjected to a tensile stress, cavities do not form as a result of liquid rupture but are caused by the rapid growth of these nuclei.
To understand this, think of beer (or champagne if you prefer) in a bottle, when it is closed you see no air bubbles and the liquid does not look fizzy. The pressure in the bottle is above the saturation pressure of the gas in the liquid. When you open thinatpheftoyore
Tvaph
Bca
Inmphor
Inhianenknai
Inprto45
e bottle suddenly bubbles appear and so the dissolved gas (molecules of gas held the liquid) starts to appear as gas. In fact the liquid is gas saturated and the mospheric pressure is less than the saturation pressure of the liquid. This enomenon is clearly not cavitation but air release (aeration). Considering nuclei fects, bubbles form only at particular places in your glass: around the glass (due small asperities) and round any particles present in the liquid. Theoretically, if ur liquid was perfectly pure and the wall of the system perfectly regular, air lease or cavitation would occur with great difficulty!
he key point about cavitation is that it is a phase change: the liquid changes to por. A comparison can be made between cavitation and boiling. If we look at the ase diagram below:
Figure 2.48: Cavitation and boiling
oiling is a phase change at constant pressure and variable temperature and vitation is a phase change at constant temperature and variable pressure.
any system air release starts first and if the pressure decreases further, cavitation ay occur. This means that, sometimes, people talk about cavitation when the real enomenon is air release. Both phenomena can lead to destruction of the material component.
both cases it is entrained gas that causes the troubles. When cavities encounter gh pressure in the downstream circuit, these bubbles or cavities can be unstable d can collapse implosively. The pressure developed at collapse can be large ough to cause severe mechanical damage in the containing vessel. It is well-own that hydraulic pumps and pipework can be badly damaged by cavitaton and
r release.
all classical hydraulic systems air release and cavitation must be avoided to event material destruction but sometimes it is required like for injection systems prepare the spray formation.
pressure
boiling liquid
solid gas (vapor)
cavitation
temperature
-
Chapter 2Theory of fluid properties
46
2.3 Viscosity
Figure 2.49: Viscosity
Vbooia
Vthmse
Fve
Tdi
wvealviolhe
Tiniscosity is a measure of the resistance of the fluid to flow. This characteristic has th positive and negative effects on fluid power systems. A low viscosity leads to l leaks in the dead zone formed between the mechanical parts in movement, and high viscosity will lead to loss of pressure in hydraulic ducts.
iscosity is a characteristic of liquids and gases and is manifested in motion rough internal damping. Viscosity results from an exchange of momentum by olecular diffusion between two layers of fluid with different velocities. In this nse, the viscosity is a fluid property and not a flow property.
igure 2.49 shows the relation between shearing constraint and difference of flow locity between two layers
he definition of viscosity was first given by Newton. Between two layers of stance dy, the exerted force between these two layers is given by:
here U(y) is the velocity depending on the radial position y and dU/dy the locity gradient. This proportionality expresses the notion of Newtonian fluid and lows the introduction of defined as the dynamic viscosity or the absolute scosity. The dimension of is [ML-1T -1] and the SI unit is kg/m/s or Pa s. The der unit is the Poise, P, which is 0.1 kg/m/s. However, this is very small and nce the milli Poise, mP, is the common unit which is 10-4 kg/m/s.
he dynamic viscosity is the constant of proportionality between a stress and the tensity of shearing between two neighboring layers:
yu + du
u
F AdU y( )dy---------------=
dU y( )dy---------------and hence
dU y( )dy
------------------------------ shear stress
shear rate---------------------------= = =
-
Hydraulic Library 4.2User Manual
However the absolute viscosity is not very often used in fundamental equations. For example the dynamics of the elementary volume between the two layers is expressed as:
an
Inthin
ofviunpa
N
Ngr
Addy------dy AdydU y( )dt---------------=47
d thus using the shear stress calculation:
other formulas (e.g. Navier Stokes) the ratio between the absolute viscosity and e density occurs so often that a new parameter called the kinematic viscosity is troduced
dimension [L2 T -1] and so the SI unit is the m2/s. The older unit of kinematic scosity is the Stoke, St, which is 10-4 m2/s. However, even this is a very small it and hence the centistoke cSt is the common unit with 1 cSt = 10-6 m2/s. This rameter is easily measured with viscometers.
ote that the viscosity varies significantly with the fluid temperature.
Figure 2.50: Viscosity against temperature
ormally in absence of air release and cavitation the variation with pressure is not eat unless the pressure is very extreme.
dU y( )dt
--------------- 1---ddy------ ---
d2U y( )dy2
------------------= =
---=
-
Chapter 2Theory of fluid properties
48
Figure 2.51: Variation with pressure
Viscosity
Ath
Tde
TnuhyFlo influence on the flow
nother important aspect of the viscosity is its influence on the flow conditions of e fluid. We can distinguish two types of flow conditions:
Laminar flow for which the flow lines are parallel and shearing forces create a pressure drop.
Turbulent flow for which the fluid particles have a disordered, random movement leading to a loss of pressure.
hese two conditions can be distinguished using the Reynolds number which is fined as follows:
with
U: average fluid velocity
d: diameter of the duct (hydraulic diameter for others geometries)
: density: dynamic viscosity: kinematic viscosity
he transition between laminar to turbulent flow occurs at the critical Reynolds mber. This is not well defined, there exists always a transition region. In a draulic line, the critical Reynolds number is generally between 1500 to 2000.
or uneven geometries (thin-walled orifices), the critical Reynolds number can be wer than 100.
Re Ud-----------Ud
------- inertia effectsviscous effects-----------------------------------= = =
-
Hydraulic Library 4.2User Manual
For non-circular cross sections, the hydraulic diameter can be used to determine the Reynolds number. Hydraulic diameter is defined as follows:
We now give two examples:
H
Flow thro
Onudeth
w
Twor
CCco
ne
veC
Itex
dh4xcross sectional area
wet perimeter-----------------------------------------------------=49
Circular orifice of diameter:
Rectangular orifice, length L and width l:
ence if .
ugh orifices
rifices (also called restrictions) can be fixed or variable and occur in huge mbers in fluid systems. Not surprisingly in Engineering courses a mathematical scription is presented. This is usually based on Bernoullis equation and leads to e form
here Cq is the flow coefficient. This is variously described as
typically 0.7 or
varying with orifice geometry and Reynolds number.
he second alternative is obviously more correct. If we do take a constant value, e are forced to have the gradient of Q against infinity at the igin! This cannot be and if you try to implement it is a numerical disaster!
learly the flow is laminar for sufficiently small pressure drops which means that q is certainly not constant. One solution is to perform detailed experiments and mpute Cq against Reynolds number. In the context of the orifice (not
cessarily circular) the Reynolds number is where U is a mean
locity and dh the hydraulic diameter. If we take U=Q/A, we end up with the form q =f(Q) and ultimately with
is possible to work with an implicit relationship like this but we would prefer an plicit formula.
dh
4 d2---
2
d------------------ d= =
dh4Ll
2 L l+( )-------------------2LlL l+-----------= =
dh 2l L l
Q CqA2 Pup Pdown( )
---------------------------------------=
P Pup Pdown=
ReUdh
----------=
Q F Q( )=
-
Chapter 2Theory of fluid properties
50
This is provided by introducing another dimensionless number known as the flow number and denoted by [reference 1]. This is defined as
From a modeling point of view contains quantities we know. Using we have
an
toth
Fli
an
W
dh-----
2 Pup Pdown( )---------------------------------------=d provided we have , we have an explicit relationship which is easy
evaluate. There are no more problems to obtain measurements for an for and so the flow number form has many advantages.
or general use the AMESim submodels OR000 and OR002 require and the miting value of Cq. The value of is computed from
d the flow coefficient is calculated as
hen , Cq is about 96% of Cqmax.
Note: Both Pup and Pdown are needed. P is not enough because a pressure drop of 1 bar to 0 bar is not the same as 1001 bar to 1000 bar.
It is not clear which pressure should be used to calculate and . Possibilities are Pup, Pdown, (Pup + Pdown)/2 or the pressure at the vena contracta. AMESim uses (Pup + Pdown)/2.
Tables of can also be compiled using CFD (computational fluid dynamics) software.
For high values of , Cq is approximately constant. The lowest value of at which Cq is approximately
constant is called the critical flow number .
The critical flow number for a thin or sharp edge orifice is about 100 and for a long orifice is about 3000.
For a long edge orifice the constant Cq value is also the maximum value.
For a sharp edge orifice the maximum Cq value can be slightly greater than the constant value and occurs at a value slightly below .
QCqA
dh-----------------=
Cq Cq ( )Cq Cq ( )
Cq Cq Re( )=
Cq Cq ( )
crit
crit
crit
dh-----
2 Pup Pdown( )---------------------------------------=
Cq Cqmax 2
crit-----------
tanh=
crit=
-
Hydraulic Library 4.2User Manual
Figure 2.52 shows a graph of Cq against .Figure 2.52: Graph of Cq against l
Frictional
Sstcr
w
Inoflaviha
Cqmax
D51
drag
ubmodels belonging to this category are used to model resistance to flow in raight tubes and conduit. The pressure losses along a straight tube of constant oss-sectional area are calculated from the Darcy-Weisbach equation:
here:
friction coefficient of the segment of relative length l/Dh=1;
hydraulic or equivalent diameter;
length of flow segment.
For this type of submodel the total friction factor is given by:
straight tubes, the resistance to the motion of a liquid or a gas under conditions laminar flow is due to the force of internal friction. This happens when one yer of the liquid (or gas) has a relative motion compared to the others. These scosity forces are proportional to the first power of the flow velocity. We then ve:
crit
p lD----Q2
2Amin2
-------------- =
h
l
lDh------=
Re( )
-
Chapter 2Theory of fluid properties
52
As the Reynolds number increases, the inertia forces, which are proportional to the velocity squared, begin to dominate. As flow becomes turbulent, there is a significant increase in the resistance to motion. Part of this increase is due to the roughness of the wall surface. Therefore, we have:
w
rr
Tto
T
w
AB
Re rr( , )
Dhere:
is the relative roughness.
he relative roughness is calculated as the ratio of the average height of asperities the tube diameter. See details in Figure 2.53 :
Figure 2.53: Relative roughness
he relative roughness of a pipe is given by:
here:
is the equivalent uniform roughness of the pipe;
is the hydraulic diameter of the pipe.
sampling of absolute pipe roughness for new clean pipes is proposed by inder [2]:
Pipe material Absolute roughness
Drawn brass 1.5
Drawn copper 1.5
Commercial steel 45
Wrought iron 45
l
Dh
rr Dh------=
h
m
m
m
m
-
Hydraulic Library 4.2User Manual
TreN
F
A
Reference
[1In
[2E
Asphalted cast iron 120
Galvanized iron 150
Pipe material Absolute roughness
m
m53
he dependence of the friction coefficient on the Reynolds number and the lative roughness as shown in Figure 2.54 is often known as the harp of ikuradse.
igure 2.54: Evolution of the frictional drag factor with the Reynolds number and with relative roughness
ll lines with friction in the hydraulic category use such a frictionnal drag factor.
s
] McCloy D, Discharge Characteristics of Servo Valve Orifices, 1968 Fluid ternational Conference, pp 43-50.
] R.C. Binder, Fluid Mechanics. 3rd Edition, 3rd Printing. Prentice-Hall, Inc., nglewood Cliffs, NJ. 1956.
Cast iron 260
Wood stave 0.2 to 0.9
Concrete 0.3 to 3
Riveted steel 0.9 to 9
m
mm
mm
mm
Re64
= rr = 0.033
rr = 0.016
rr = 0.008
rr = 0.004 rr = 0.002 rr = 0.001
rr = 0.0
-
Chapter 2Theory of fluid properties
54
-
Hydraulic Library 4.2User Manual
Chapter 3: AMESim Fluid Properties
3.1 In
Afl
YTde
Wan
FP04
el
Ticwpr
Ww55
troduction
MESim allows you to use systems with several fluids in a single sketch. For each uid you use, you need to add a fluid property icon to your sketch.
ou can use each icon to install an index of hydraulic fluid in the range 0 to 100. hese icons give you access to a number of submodels which will now be scribed.
e do not describe here the FPDROP submodel since it is considered as obsolete d it is available only for compatibility with old systems (4.0 and earlier).
ementary
This has the following parameters:
index of hydraulic fluid
density
bulk modulus
absolute viscosity
saturation pressure (for dissolved air/gas)
air/gas content
temperature
polytropic index for air/gas/vapor content
absolute viscosity of air/gas
name of fluid
his submodel which is associated with the on shown has an enumeration parameter hich gives you access to a collection of fluid operties of varying levels of complexity.
e now describe the parameters associated ith each enumeration option.
-
Chapter 3AMESim Fluid Properties
56
This option makes the following assumptions:
1. The bulk modulus of the liquid with zero air/gas content is constant. This means the corresponding density varies exponentially with pressure.
2. The viscosity of the liquid with zero air/gas content is constant.
Tte
ad
T
Wrethbopuprasm
Wreto
Wfo
ad
Tvifohere is an air release and cavitation model included. Note that name of fluid is a xt string (e.g. cooling water) that identifies the fluid.
vanced
his is similar to elementary but there are addition parameters.
hen the pressure reaches the saturation pressure of the fluid, some air/gas is leased. If the pressure continues to decrease, the high saturated vapor pressure of e fluid can be reached and some vapor appears (cavitation: the liquid starts to il). Remember the fluids used in real engineering systems are not chemically re substances. For this reason cavitation is assumed to occur over a range of essures and the low saturation vapor pressure is the pressure at which it is sumed that all liquid has become vapor. All these changes of state strongly odify the fluid characteristics.
ith the elementary option this behavior is taken into account with some asonable fixed cavitation parameters. However, with advanced you are allowed set these values yourself. They are:
high saturated vapor pressure
low saturated vapor pressure
absolute viscosity of vapor
effective molecular mass of vapor
ith the elementary option, advanced user parameters do not appear, but have the llowing constant values:
high saturated vapor pressure: 0.9 bar
low saturated vapor pressure: 0.95 bar
absolute viscosity of vapor: 0.02 cP
effective molecular mass of vapor: 200
vanced using tables
his option has been created to use values for the fluid density, bulk modulus and scosity depending on the current pressure and temperature. The parameters used r this submodel are the following ones:
index of hydraulic fluid
saturation pressure (for dissolved air/gas)
-
Hydraulic Library 4.2User Manual
air/gas content
temperature
polytropic index for air/gas/vapor content
absolute viscosity of air/gas
NTfutafl
T
Y
Ed
T57
(advanced users) high saturated vapor pressure (cavitation)
(advanced users) low saturated vapor pressure (cavitation)
(advanced users) absolute viscosity of vapor
(advanced users) effective molecular mass of vapor
(advanced users) air/gas density at atmospheric pressure 0 degC
name of fluid
name of file specifying fluid properties
ote that density, bulk modulus and viscosity do not appear in the parameters. hey are calculated from values in tables processed by specific functions. These nctions apply interpolation processing to calculate the fluid characteristics from bles. These tables are given in a text file specified by the name of file specifying uid properties parameter of the submodel.
hree samples of such files are supplied in the AMESim installation CD:
tblprop1.txt
tblprop2.txt
tblprop3.txt
ou should be able to copy these files from the directory:
$AME/misc (for Unix) or
%AME%\misc (for Windows)
ach file describes a particular mode of definition of the fluid properties. For the ensity and bulk modulus, three modes are available:
In mode 1, density and bulk modulus are defined from a reference density, a reference pressure and a set of tables of bulk modulus values against pressure. Each table is written for a given temperature (see tblprop1.txt).
In mode 2, density and bulk modulus are defined from a set of tables of density values against pressure. Each table is written for a given temperature (see tblprop2.txt).
In mode 3, density and bulk modulus are defined from a reference density, a reference pressure and a set of tables of speed of sound values against pressure. Each table is written for a given temperature (see tblprop3.txt).
he viscosity of the fluid is also given in these files after the definition of the
-
Chapter 3AMESim Fluid Properties
58
density and the bulk modulus. Two modes are available for the viscosity:
In mode 1 the absolute viscosity is defined from tables of absolute viscosities in cP. Each table is written for a given temperature (see tblprop1.txt).
In mode 2 the absolute viscosity is defined from tables of kinematic
Tarthcaus
R
TB
Itnotere
Tinacviscosities in cSt. Each table is written for a given temperature (see tblprop2.txt).
he best plan if you want to use this facility is to copy these files to a suitable local ea and examine them in an editor. Lines beginning with a '#' are comments and ese comments give further information on how the data is arranged. Then you n select the file that uses the modes you find suitable, and modify it in order to e your own data.
obert Bosch adiabatic diesel
his submodel is for Diesel fluid properties and is kindly supplied by Robert osch GmbH. Its parameters are:
fuel type
index of hydraulic fluid
(advanced users) high saturated vapor pressure
(advanced users) low saturated vapor pressure
(advanced users) effective molecular mass
absolute viscosity of air/gas
(advanced users) absolute viscosity of vapor
air/gas content
temperature
is assumed that these fluids are used in fast acting injection systems and there is time for the air content to dissolve or undissolve. The user sets a fixed
mperature and the local temperature is computed using an approximate lationship for an adiabatic change.
he fuel type is an enumeration teger parameter which gives cess to 9 diesel fuels.
-
Hydraulic Library 4.2User Manual
simplest
This option gives the simplest hydraulic fluid properties. Its parameters are:
index of hydraulic fluid
density
Tduot
3.2 T
Cdidi
Fofa sp
SF
Nth59
bulk modulus
absolute viscosity
saturation pressure
air/gas content
temperature
polytropic index for air/gas content
name of fluid
his submodel can be useful in difficult cases. The integrator has an easier task ring cavitation and air release and so it may be possible to get a solution when her methods are unsuccessful.
utorial example
opy $AME/misc/tblprop1.txt or %AME%\misc\tblprop1.txt in a suitable rectory. Start AMESim an build the system shown in Figure 3.55 in this same rectory
Figure 3.55: A simple system for plotting fluid properties.
P04 is the only submodel available for this icon. Set the index of hydraulic fluid the hydraulic submodels to 1. Change the parameters of the pressure input to get ramp from 0 to 100 bar in 10 seconds. Change the parameter name of file ecifying fluid properties so that it specified your own file tblprop1.txt.
tart a simulation and plot the density, the bulk modulus and the viscosity of the PROP submodel against the pressure:
ow edit the values of your file tblprop1.txt and rerun the simulation. Note how e properties change.
-
Chapter 3AMESim Fluid Properties
60
-
Hydraulic Library 4.2User Manual
Chapter 4: Selecting submodels for Hydraulic Lines
Ttrfo
4.1 In
Tm
Tfo
Lumped p
NwInwsisuO
Capacitan
Lm61
his problem can create a lot of worry for some users and hence in this chapter we y to give some pragmatic rules to help you select an appropriate submodel. In the rmulae below it is assumes that SI units are used throughout.
troduction
he submodels for lines are arranged with the simplest submodel at the top and the ost complex at the bottom. Why are there so many line submodels?
he main problem is the complexity of flow in hydraulic pipes and hoses. The llowing features might be important:
Variation of fluid compressibility and expansion of pipe/hose walls with pressure.
Inertia of fluid.
Variation of bulk modulus with pressure.
Variation of viscosity with pressure.
Laminar, turbulent and transition flow.
Frequency dependent friction.
Air release and cavitation.
arameter and distributed parameter line submodels
ormally it is perfectly satisfactory to use a lumped parameter submodel in hich properties such as pressure are represented by a single representative value. other words, within the pipe we assume there is negligible variation of pressure ith position. However, if the pipe is extremely long or if wave dynamics are gnificant a distributed parameter submodel should be used. For this type of bmodel, there are quantities such as pressure calculated at a number of positions. ften these are stored as an array of values.
ce, resistance and inertia in line submodels
iquids are compressible - not as much as gases but it is essential to recognize in odeling that they are compressible. When subjected to high pressure the density
-
Chapter 4Selecting submodels for Hydraulic Lines
62
of a liquid increases. In addition the pipe or hose containing the liquid expands with pressure. The net result of this is a capacitance (spring) effect.
To cause a hydraulic fluid to travel along a horizontal pipe we must provide a pressure gradient to drive the fluid. This is a resistance effect.
The moving fluid has mass and hence it has inertia.
Zero-dime
Tas
Tth
One-dime
Ifacdim
Ufoeq
Thacommw
Yroimco
Hfogiesnsional line submodels
he simplest line submodels are DIRECT and HL000 and these can be described zero-dimensional.
he DIRECT line submodel assumes that the two ports are very close together and e fluid and pipe in between contributes nothing.
nsional line submodels
we take into account resistance or inertia, we have a one-dimensional (or more curately one spatial dimension) submodel. It is possible to have 2- and 3-mensional submodels but this is the field of computation fluid dynamics. For odeling systems one dimension is normally enough.
ntil AMESim 4.2 all hydraulic line submodels in AMESim employ an implicit mulation which was very stable. This ensures that when the system is close to uilibrium very large integration steps may be taken.
here are other classes of methods which have much more limited stability and ve a strict limitation on step size due to the CFL (Courant Friedrichs Lewy) ndition. AMESim 4.2 contains an experimental version of one of these ethods. These are HLG20, HLG21 and HLG22 and implement the Godunov ethod. They should be used as alternatives to HL020, HL021 and HL022. In other ords for very low viscosity situations.
ou can try these methods but they tend to give slow simulation runs and are less bust that the regular line submodels. They will be replaced soon by an plementation of a different method which, although still restricted by the CFL ndition, is faster and more robust than Godunov.
L000 considers the capacitance only. The length of the line is too small r significant resistance. The fluid velocity and the mass are too small to ve significant inertia. The hydraulic chamber submodel HC00 is sentially the same as HL000.
-
Hydraulic Library 4.2User Manual
4.2 Line submodels occur in threes
AMESim line submodels normally occur in groups of three. The reason for this is the input and output characteristics of external variables of a submodel. If we cofo
Incoarfithli
4.3 T
Aspect ra
Ton
p
p
p63
nnect a pipe to a component, the component submodel normally does one of two llowing things at the connection port:
It calculates the flow rate (output) from the pressure (input) or
It calculates the pressure (output) from the flow rate (input)
each case the pipe submodel must provide the correct variable for the mponent submodel. Figure 4.1 shows the three standard possibilities. The rows indicate the direction of the flow of information. Thus in the left port of the rst case the flow rate is calculated by the attached component submodel whereas e pressure is calculated by the line submodel. The important feature is that most ne submodels are produced in groups of three to cover the three possibilities.
Figure 4.56: Three different causalities.
hree important quantities
tio
he checking algorithm in AMESim issues warning messages when you use a e-dimensional submodel that has an aspect ratio length/diameter ratio less than
pressure
ressure flow rate
flow rate
flow rate
ressure
flow rate pressure
pressure
flow rate flow rate
ressure
-
Chapter 4Selecting submodels for Hydraulic Lines
64
6. This is defined in terms of the length L and diameter D as follows:
Short fat pipes require different submodels than long thin pipes.
For distributed line submodels the line is divided into a collection of cells and the te
Dissipatio
A
Wwin
Ifw
AratioLD----=st is that the cell length/diameter ratio must not be greater than 6.
n number
nother important measure is the dissipation number. This is defined as:
the kinematic viscosity and c the speed of sound.
hen the dissipation number reaches 1, the principal eigenvalues become real and ave effects are not significant. When this is true no models that take into account ertia should be used.
the dissipation number is significant less than 1, it may be important to consider ave effects. This motivates the following table.
Line submodels
Takes into account:
Dissipation number
Use for
HL01HL02HL03
Capacitance + resistance
> 0.8 relatively short pipes with high dissipation number
HL004 HL005 HL006
Capacitance + resistance + inertia+ frequency dependent friction
< 1.2 and > 0.8e-3 relatively short pipes with low dissipation number
HL04HL05HL06
Capacitance + resistance + inertia
< 1.2e-3 relatively short pipes with very low dissipation number
HL02IResistance + inertia
relatively short pipes with very high fluid velocity
Ndiss4LcD2----------=
c B----=
-
Hydraulic Library 4.2User Manual
T
Commun
T
IfseTdi
War
TasZ
WT
HL10HL11HL12
Capacitance + resistance
> 0.8 long pipes with high dissipation number
Capacitance + resistance +
< 1.2 and > 0.8e-3 moderate lengths with low but not very low 65
he result of this test must be qualified by considering the next important number.
ication interval
he time taken for a wave to travel down the pipe is
this time is significantly less than the communication interval, you will never e the waves in plots and so it is not useful to use a wave dynamics submodel. his is why changing the communication interval leads to the appearance/sappearance of warning messages.
e will compare Twave with the communication interval Tcom to determine if we e likely to see waves in results.
he hydraulic volume submodels HC00 and HC01, which are basically the same HL000, are included for completeness. Similarly the zero volume submodel
EROHV is also included.
Figure 4.57: Three other lines.
e are about to display charts which help decide which line submodel to select. hese must be studied bearing in mind the following notes.
HL030 HL031 HL032
inertia+ frequency dependent friction
dissipation number
HL020 HL021 HL022
HLG20
HLG21
HLG22
Capacitance + resistance + inertia
< 1.2e-3 moderate lengths with very low dissipation number
TwaveLc---=
Hydraulic chambers Zero volume
-
Chapter 4Selecting submodels for Hydraulic Lines
66
Note: Since many lines in AMESim are constituted with several segments (1, 5, 10 or 20) it can be noticed that in the below selection process the aspect ratio is compared to 6, 30 or 60 corresponding to segment of length L, L/6, and L/10.
The decision process employed in the charts that follow is very
similar to one employed by AMESim when it checks the suitability of your submodels. If the submodel is regarded as unsuitable, a warning message is issued, values of aspect ratio and dissipation number are given for one segment of the choosen line.
Often the final result from the chart is three submodels such as HL01, HL02, HL03. Since AMESim will check causality, only one submodel, the one that is compatible with adjoining submodels, will be offered to you as a choice.
The charts are intended for general guidance and give a good choice most of the time. However, there are circumstances in which an advanced users may wish to break these rules.
-
Hydraulic Library 4.2User Manual
4.4 The selection process
HL0HL0HL0HL1HL1HL1
1
Y
FrictioDomin67
Aratio > 6?
Yes No
Suitable for one spatialdimension assumption
Zero spatial dimensions
Ndiss
> 1
123012
es
nated
No
Possibly waveeffects important
DIRECTHL000HC00HL02IZEROHV
HL01 HL02 HL03HL04 HL05 HL06HL004 HL005 HL006HL10 HL11 HL12HL020 HL021 HL022
HL030 HL031 HL032
2
3
HLG20 HLG21 HLG22
Aratio > 60
Use distributed submodelHL10HL11HL12
Yes
Use lumped submodelHL01HL02HL03
Friction dominated pipes and hoses
No
1
-
Chapter 4Selecting submodels for Hydraulic Lines
68
T < T
Cannot see wregard as fric
1
Ff
Possibly wave effects important
NoYes
2wave com
aves sotion dominated Aratio
> 30
Use simplestwave submodel
HL04HL05HL06
No
Ndiss > 1.0e-3
Yes
requency dependentriction important
HL030HL031HL032
Yes
HL020HL021HL022
No
Use frequencydependent friction
HL004HL005HL006
Ndiss > 1.0e-3Yes No
HLG20HLG21HLG22
High flowacceleration?
Try to useHL02I
Yes
Try to use HL000(or HC00) orDIRECT
No
Problem withslow simulation. Suspectcaused by pressure in submodel.
Try ZEROHV
Yes
Short pipes
High fluid inertia.
3
-
Hydraulic Library 4.1User Manual
Index
AadvanceAdvancadvanceAir releAspect
BBatch pBatch ruBernoulBulk m
VaCCapacitCavitatiCFL coCommoCommu
ImCompleCompreCompreCourantCubical
DDensityDiesel
FluDIRECDissipaDissolvDistribuDuty cy
EelementEntrapp
FFlow .
LaFlow raFluid coFluid prFP04 .Frequen69
d . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56ed properties . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12d using tables . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56ase . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2, 9, 10, 44, 45, 61ratio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
arameters . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13n . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13lis equation . . . . . . . . . . . . . .