8 Bond Graph Models for Hydraulic Systems

8.1       Overview

The generalized BG elements and relations apply to the modelling of dynamics of hydraulic systems in a similar way that the mechanical or electrical systems were treated; i.e.; they are analogous (see Table 3‑1). In this chapter, we define the effort and flow for hydraulic systems and derive the relations for hydraulic capacitance, inertance, and resistance corresponding to BG elements C, I, and R, respectively. Note that the complexity of fluid behaviour in static or dynamic flow conditions require us to pay more attention to identify these quantities and relations as compared to those for mechanical and electrical systems.

For modelling hydraulic systems, we are usually interested in having a relationship between pressure and fluid volume in static conditions and between pressure drop and fluid volume flow rate in dynamic conditions. For example, we might be interested to know the pressure drop for a given flow rate in a pipe, or we might want to know the pressure at a given depth in a storage tank, as sketched in Figure 8‑1.

Figure 8-1 Sketches for pressure drop in a pipe and in a storage tank

8.2       Definitions of Effort, Flow, and Momentum for Hydraulic Systems

Recall that power is the quantity of interest in BG method (see section 3.2). Considering a control volume \mathbb{V}  of an incompressible fluid flowing under pressure \mathbb{p} , we can write the power \mathcal{P} as the product of the force F  exerting on the fluid, resulted from applied pressure, and the velocity v of the fluid flowing through the volume, or \mathcal{P}=Fv. But the velocity of the fluid can be written as v={\dot{\mathbb{V}}}/{\mathbb{A}}, using the continuity relation, where \dot{\mathbb{V}}  is volume flow rate of the fluid and \mathbb{A}  is the cross-sectional area of the control volume. Therefore, after substitution, we get \mathcal{P}  =  F \dfrac{\dot{\mathbb{V}}}{\mathbb{A}} = \dfrac{F}{A} \dot{\mathbb{V}}  =  \mathbb{p}\dot{\mathbb{V}}, or equivalently rate of energy \mathbb{pV} . Comparing the relation \mathcal{P} =  \mathbb{p}\dot{\mathbb{V}} with the BG generalized relation for power, i.e., \mathcal{P} = ef , we can write e  \equiv  \mathbb{p}  and f  \equiv  \dot{\mathbb{V}}. In other words, for hydraulic systems, pressure is equivalent to BG effort, and fluid volume flow rate is the BG flow. Similarly, we can write the generalized BG displacement q  as the volume of the fluid, or q = \int f dt = \int \dot{\mathbb{V}} dt = \mathbb{V}.

In BG method, the generalized momentum is the integral of effort. Therefore, we can write the fluid momentum p_{\mathbb{P}} as the integral of pressure, or p_{\mathbb{p}}  =  \int \mathbb{p}  dt . Summarizing these relations, we have

    \begin{equation*} \begin{dcases} e \equiv \mathbb{p}, pressure \\ f \equiv \dot{\mathbb{V}} , volume \: flow \: rate \\ e \cdot f \equiv \mathbb{p} \dot{\mathbb{V}}, power \\ q \equiv \mathbb{V} = \int \dot{\mathbb{V}}dt, volume \\ p \equiv p_{\mathbb{p}} = \int \mathbb{p} dt, pressure \:  momentum \\ \end{dcases} \end{equation*}

For hydraulic systems, fluid pressure is equivalent to generalized BG effort, and fluid volume flow rate is the generalized BG flow.

 

8.3       Fluid Compliance: C-element

Fluid compliance or hydraulic capacitance describes potential energy storage with a fluid, e.g., the height of fluid in a tank. It is equivalent to mechanical spring compliance or electrical capacitor capacitance. For a C-element in BG method, we have e=q/c . Using equivalent quantities for hydraulic systems, we can write \mathbb{p} = \mathbb{V}/c, or fluid compliance is volume change per unit of pressure acting on the fluid volume. For an incompressible fluid with density \rho, the hydrostatic pressure at depth h is \mathbb{p} = \rho gh and the volume of the fluid is \mathbb{V} = \mathbb{A}h . After substitution, we get \mathbb{p} = \mathbb{A}h/c , or after rearranging and simplifying, the hydraulic capacitance for incompressible fluid is

(8.1)   \begin{equation*} c = \dfrac{\mathbb{A}}{{\rho}g} , incompressible \end{equation*}

where g is the gravitational acceleration. The dimension of fluid compliance can be worked out as  [\dfrac{m^4s^2}{kg}] = [\dfrac{m^5}{N}] = \dfrac{volume}{pressure}.

Note that the pressure could be replaced by total dynamic pressure for fluid in motion.

If the fluid is compressible, we use the bulk modulus of elasticity \beta for calculating the change in volume. By definition, \beta is pressure needed to change fluid volume per unit of volume, or \beta = \vartriangle \mathbb{p}/(\vartriangle\mathbb{V}/\mathbb{V}).  Therefore, c = \dfrac{\vartriangle\mathbb{v}}{\vartriangle\mathbb{p}} = \dfrac{\mathbb{v}}{\beta}.

 

(8.2)   \begin{equation*} c = \dfrac{\mathbb{V}}{\beta} , compressible \end{equation*}

 

For more complex flow and non-uniform, flexible tubes, consult with chapter 4 of Dean, Karnopp, Margolis, and Rosenberg [20].

Having the hydraulic capacitance, we can write the relation between the flow rate and the pressure as \vartriangle\dot{\mathbb{V}} = \dot{\mathbb{V}}_2 - \dot{\mathbb{V}}_1 = c\dfrac{d\mathbb{p}}{dt}, useful to calculate the flow rate for given pressure. Similarly, we can write \mathbb{p} = \dfrac{1}{c}\int\vartriangle\overset{.}{\mathbb{V}}dt = \dfrac{1}{c}\int( \dot{\mathbb{V}}_2 - \dot{\mathbb{V}}_1)dt, useful for calculating pressure for given flow rates.

 

(8.3)   \begin{equation*} \begin{dcases} \vartriangle\dot{\mathbb{V}} = \dot{\mathbb{V}}_2 - \dot{\mathbb{V}}_1 = c\dfrac{d\mathbb{p}}{dt}\\ \mathbb{p} = \dfrac{1}{c}\int( \dot{\mathbb{V}}_2 - \dot{\mathbb{V}}_1)dt\\ \end{dcases} \end{equation*}

Note the similarity between relations given by Equation (8.3) and those given for mechanical spring (when pressure is replaced by force and fluid volume flow rate by velocity) and electrical capacitance (when pressure is replaced by voltage and fluid volume flow rate by current).

8.4       Fluid Inertia: I-element

Fluid inertia, or hydraulic inertance, describes kinetic energy storage with a fluid or the inertia, e.g., of a fluid flowing in a pipe. It is equivalent to inertia related to mass in mechanical or inductance in electrical systems. For an I-element in BG method, p = If describes the relation between generalized momentum and flow. Using equivalent quantities for hydraulic systems, we can write p_\mathbb{p} = I\dot{\mathbb{V}}, or fluid pressure momentum is the product of fluid inertia by its volume flow rate. To derive the relation for I, we require to have the relationship between the momentum and volume flow rate of the fluid flow. For derivation, we consider a control volume with length \ell  and cross-sectional area \mathbb{A} of the fluid with density \rho , as shown in Figure 8-2

Figure 8-2 Sketch for a control volume of flowing fluid in a pipe

Assuming a pressure difference \vartriangle\mathbb{p} between two ends of the control volume acting on the fluid, we can write Newton’s second law for the fluid motion as \underbrace{\vartriangle\mathbb{pA}}_{force} = \underbrace{\rho\mathbb{V}}_{mass}\dfrac{d}{dt} \underbrace{\left(\dfrac{\dot{\mathbb{V}}}{\mathbb{A}}\right)}_{velocity}. Rearranging the terms and integrating the pressure, we get the pressure momentum p_{\mathbb{p}} = \int {\mathbb{p}dt = \dfrac{\rho{\mathbb{V}}}{\mathbb{A}}\dfrac{\dot{\mathbb{V}}}{{\mathbb{A}}}. But volume is \mathbb{V} = \mathbb{A}\ell. After substituting, we get the relationship between pressure momentum and the volume flow rate as

(8.4)   \begin{equation*} p_\mathbb{p} = \dfrac{\rho\ell}{\mathbb{A}}\dot{\mathbb{V}} \end{equation*}

Comparing Equation (8.4) with the generalized momentum equation for I-element, p = If we can write the fluid inertia as

(8.5)   \begin{equation*} I= \dfrac{\rho\ell}{\mathbb{A}}, incompressible \end{equation*}

The dimension of fluid inertia can be worked out as [\,\dfrac{kg}{m^4}]\, \equiv \dfrac{pressure\: momentum}{volume\: flow \:rate}.

From Equation (8.5), we can conclude that a fluid has larger inertia when flowing in small diameter tubes, compared to in larger tubes because I is inversely proportional to \mathbb{A}. This effect is counterintuitive and is sometimes misinterpreted with the wrong assumption that large-size tubes should exhibit larger inertia effects. Note that we consider only force due to pressure, and not friction due to viscosity, assuming an ideal fluid.

For more complex flow and non-uniform and/or flexible tubes, consult with chapters 4 and 12 of Dean, Karnopp, Margolis, and Rosenberg [20]. For example, if the cross-section of the pipe and the density of fluid change along its x-axis, then we get

(8.6)   \begin{equation*} I = \int\dfrac{pdx}{{\mathbb{A}}(x)} \end{equation*}

Having the hydraulic inertance, we can write the relation between the flow rate and the pressure as  \vartriangle\dot{\mathbb{V}} = \dot{\mathbb{V}}_2 - \dot{\mathbb{V}}_1 = \dfrac{1}{I}\int\mathbb{p}dt = \dfrac{\mathbb{A}}{\rho\ell}\int\mathbb{p}dt, useful to calculate the flow rate for given pressure. Similarly, we can write \vartriangle\mathbb{p} = I\dfrac{d\dot{\mathbb{V}}}{dt} = \dfrac{\rho\ell}{\mathbb{A}}\dfrac{d\dot{\mathbb{V}}}{dt}, useful for calculating pressure for given flow rates.

(8.7)   \begin{equation*} \begin{dcases} \vartriangle\dot{\mathbb{V}} = \dfrac{1}{I}\int\mathbb{p}dt\\ \vartriangle\mathbb{p} = I\dfrac{d\dot{\mathbb{V}}}{dt}\\ \end{dcases} \end{equation*}

Note the similarity between relations given by Equation (8.7) and those given for mechanical systems (when pressure is replaced by force, fluid flow rate by velocity, and inertance by mass) and electrical systems (when pressure is replaced by voltage, fluid flow rate by current, and inertance by inductance).

8.5        Fluid Resistance: R-Element

Gotthilf Hagen, (1797–1884)

Fluid, or hydraulic, resistance describes energy dissipation with a fluid, e.g., friction of a fluid flowing in a pipe. Fluid resistance is equivalent to dampers in mechanical or resistors in electrical systems. For an R-element in BG method, we have e=Rf. Using equivalent quantities for hydraulic systems, we can write

(8.8)   \begin{equation*} \mathbb{p} = R\dot{\mathbb{V}} \end{equation*}

or fluid resistance is equal to pressure change per unit volume flow rate. This relationship depends on the state of the flow (e.g., laminar, turbulent) and the fluid properties (e.g., ideal, viscous,) [26], [27], [28]. Note the similarity between relations given by Equation (8.8) and those given for mechanical systems (when pressure is replaced by force and fluid flow rate by velocity) and electrical systems (when pressure is replaced by voltage and fluid flow rate by current).

Jean-Léonard-Marie Poiseuille, (1797–1869)

To demonstrate the derivation of the relation for R, we consider a laminar flow of a viscous incompressible fluid in a pipe (so-called Hagen-Poiseuille flow) and write Newton’s second law for a cylindrical differential control volume of the fluid with length dz along the pipe axis and a cross-section with radius r , as shown in Figure 8‑3. This flow is axisymmetric, and the velocity profile changes along the radius related to a cylindrical coordinate system (r, \theta, z).

Figure 8-3 Sketch of velocity profile for a Hagen-Poiseuille flow in a pipe

For a steady flow (i.e., non-transient), we have \sumF_z = 0. But the forces applied on the fluid are due to pressures \mathbb{p}_1 at point z_1 and \mathbb{p}_2 = (\mathbb{p}_1 - d\mathbb{p}) <  \mathbb{p}_1 at point z_2 and the viscous-induced shear stress \tau . Hence, (\mathbb{p}_1 - \mathbb{p}_1 + d\mathbb{p}) \pi r^2  - \tau(2\pi{r})dz = 0 where r is the radial dimension in the x-y plane parallel to the pipe cross-section. We need the relation for fluid friction effect due to viscosity. According to Newton’s law for a viscous fluid, we have \tau = \mu\dfrac{du(r)}{dr}, assuming the shear stress due to the fluid’s viscosity be proportional to the velocity gradient along the radius with the proportionality constant being the dynamic viscosity \mu . Note that velocity profile at any cross-section of the pipe is only a function of radius, or velocity vector is \overrightarrow{u} = (0, 0, u(r))\equiv (0, 0, u) . After substitution, we have \pi r^2 d\mathbb{p} = \mu\dfrac{du}{dr} (2\pi r)dz . After simplifying and rearranging the terms, we have \dfrac{r}{2\mu}d\mathbb{p} = \dfrac{du}{dr}dz . Integrating the latter relation, noting that velocity is not a function of z, gives \dfrac{r}{2\mu}\displaystyle\int_{\mathbb{p}_1}^{\mathbb{p}_2}d\mathbb{p} = \dfrac{du}{dr}\displaystyle\int_{z_1}^{z_2}dz , or \dfrac{r}{2\mu}({\mathbb{p}_2}-{\mathbb{p}_1}) = \dfrac{du}{dr}(z_2 - z_1) . Now, we rearrange the terms and let z_2 - z_1 = \ell, the length of the control volume, and use the pressure difference as a positive constant quantity \mathbb{p} = (\mathbb{p}_1 - \mathbb{p}_2) > 0 in the direction of the fluid flow. Hence, du = -\dfrac {rdr}{2\mu\ell}\mathbb{p}. Integrating both sides (the left-hand side with respect to u and the right-hand side with respect to r ) gives u = -\dfrac{r^2}{4\mu\ell}\mathbb{p} + constant . The constant of integration can be obtained using the information at the boundary of the pipe assuming the no-slip condition, or u|_{r=D/2} = 0 , where D  is the pipe diameter. Therefore, constant = \dfrac{\mathbb{p}D^2}{16\mu\ell} . Hence, after back substitution, we get the velocity profile u = - \dfrac{\mathbb{p}r^2}{4\mu\ell} + \dfrac{\mathbb{p}D^2}{16\mu\ell} . This relation is the famous parabolic velocity profile for the flow in a pipe and can be written in its functional form as

(8.9)   \begin{equation*} u(r) = \dfrac{\mathbb{p}}{4\mu\ell}\left({\dfrac{D^2}{4} - r^2}\right) \end{equation*}

Using Equation (8.9), we can calculate the velocity for any given value of r, e.g., u|_{r=0} = \dfrac{\mathbb{p}D^2}{16\mu\ell} at the centre-line of the pipe or u|_{r=D/2} = 0 at the interior wall of the pipe.

Now, to find the volume flow rate, we integrate the velocity u(r) over the whole cross-section of the pipe using a differential area element d\mathbb{A} = rdrd\theta. Or, the volume of the fluid passing through the whole cross-section of the pipe per unit of time is given by

\dot{\mathbb{V}}=\iint ud{\mathbb{A}}=\int_o^{2\pi} d\theta \int_0^{D/2}urdr = 2\pi \dfrac{\mathbb{p}}{4\mu\ell}\int_o^{D/2}\left({\dfrac{D^2}{4} - r^2}\right)rdr = \dfrac{\pi{D^4}}{128\mu\ell}\mathbb{p}, or

(8.10)   \begin{equation*} \dot{\mathbb{V}} = \dfrac{\pi{D^4}}{128\mu\ell}\mathbb{p} \end{equation*}

Comparing Equation (8.10) with Equation (8.8), we can write the fluid resistance as

(8.11)   \begin{equation*} R = \dfrac{128\mu\ell}{\pi{D^4}} \end{equation*}

The fluid resistance can be interpreted as the amount of pressure drop per unit of volume flow rate of the fluid in the pipe. The dimension of fluid resistance can be worked out as \Big[\dfrac{kg}{s.m^4}\Big] = \Big[\dfrac{N.s}{m^5}\Big]\equiv \dfrac{pressure}{volume\: flow\: rate}.

Other BG elements for hydraulic systems are sources of flow (e.g., centrifugal pumps) and efforts (e.g., reservoirs, tanks, displacement pumps). Pumps provide flow of a fluid at a certain flow rate according to their types and specifications. Reservoirs or pressure chambers provide certain pressure to the system as an effort source. The transformers elements are those like piston-cylinder (plunger), and gyrators are those elements like reaction turbines or hydraulic motors. Sketches below show some related elements.

Typical hydraulic components are shown in Table 8‑1.

Table 8-1 Typical hydraulic components and their corresponding BG elements
R-element
(valve)
C-element
(storage)
I-element
(fluid mass)
TF-element
(plunger)
GY-element
(pump)

8.6        Sign Convention for BG Modelling of Hydraulic Systems

The sign convention for hydraulic systems can be defined by specifying the relative high/low pressure points in the system and, hence, the positive fluid flow direction along the pressure drop. The pressure reference is commonly taken to be the atmospheric pressure (i.e., one atm for absolute and zero for gauge pressures). For BG modelling, it is recommended to have all pressures in gauge and define a zero-pressure point for reference atmospheric pressure. If the results are required in absolute pressure units, then one unit of atmospheric pressure can be added to the obtained values from the BG model.

8.7        Guidelines for Drawing BG for Hydraulic Systems

As mentioned in chapter 4, the general guidelines for drawing BG models can be applied to hydraulic systems, along with causality assignment rules. For hydraulic systems, we follow the guidelines given for electrical systems (see section 7.3) as described in the following steps:

1) Assign sign convention for fluid flow directions.

2) Assign 0-junction for each distinct pressure point in the system.

3) Assign 1-junction for each element in the system. This is for taking care of relative pressure drops related to each element located between two adjacent 0-junctions, since 1-junction is effort summator.

4) Select a node in the system as a reference, i.e., the atmospheric pressure point, and assign a 0-junction element to it. If gauge pressures are used, then this 0-junction and all connected power bonds can be eliminated to simplify the model.

5) Assign C-element for storage/capacitors, R-element for friction, I-element for fluid mass, and S_e for pressure and S_f for flow sources.

6) Assign TF-element for hydraulic transformers and GY-element for hydraulic gyrators.

7) Connect the elements with power bonds and assign causalities. Simplify by neglecting the bonds and the 0-junction which are connected to the 0-junction representing the atmospheric pressure.

Similarly, a 1-junction-based approach can be used for distinct flow rates and hence simplifying the BG model, as we demonstrated in the previous chapter with electrical systems.

8.8        Example: Hydraulic Reservoir-Valve System

Figure 8‑4 shows a hydraulic system consisting of two tanks, pipes, and valves. Build a BG model for this system.

Figure 8-4 Sketch for a two-tank hydraulic system

Solution:

We follow the systematic procedure for building the BG model, listed in section 8.7. For this system, we can easily identify two pressure points located at the bottom of tanks. We assign two 0-junctions for each. For flow input, we assign a flow source element, and for the output, an effort source element to define the atmospheric pressure at that location. For the tanks, we only consider capacitance, assuming slow fluid motion and neglect inertia and friction (i.e., no inertance nor resistance). For the pipe sections, we consider inertance and resistance. As well, we assign 1-junctions for flows in the pipes that represent the pressure changes for these components. Figure 8-5 shows the resulting BG model.

Figure 8-5 Bond graph model for a two-tank hydraulic system

8.9        Example: Hydraulic Reservoir-Valve System Simulation

In this example, we use the BG model developed in section 8.8, along with data assigned to parameters for simulation. Considering water as the fluid (\rho = 10^3kg/m^3 , \mu = 10^-3Pa.s) and the data given in Table 8‑2, we can calculate the related C, R, and I of the elements in the system. The diameter of the pipes is 15 cm, and g = 9.81\dfrac{m}{s^2}.

The following video shows how to build and run the model for this example in 20-sim.

Table 8‑2 Data for system components in example 8-9
Component X-section area
[m^2]
Length
[m]
C
[m^5/N]
Eq. (8.1)
I
[kg/m^4]
Eq. (8.5)
R
[N.s/m^5 ]
Eq. (8.11)
Storage tanks 2 2.04 X 10^4
Pipe1 0.01767 4 226372.4 322
Pipe2 0.01767 2 113186.2 161

 

8.10       Example: Hydraulic Pump-Reservoir-Valve System

In this example, we use the BG model developed in section 8.8, adding a pump to the system as shown in Figure 8‑6. In this example we discuss in more detail the BG model of a pump. For further details related to BG modelling of pumps, consult with references cited as [21] and [29].

Figure 8-6 A hydraulic system with a pump

The following video shows how to build and run the model for this example in 20-sim. The resulting BG model is shown in Figure 8‑7.

 

Figure 8-7 Bond graph model for the pump-reservoir-valve hydraulic system

8.11        Example: A Hydraulic Lift System

Blaise Pascal (1623-1662)

In this example, we consider a hydraulic lift, as sketched in Figure 8‑8. We build a BG model for this hydraulic system. The continuity relation applies to the fluid motion and Pascal’s law defines the pressure distribution of the fluid in the cylinders. Two transformer elements are used in the BG model to convert linear velocities of the pistons to/from volume flow rate and convert forces to pressures (see Figure 8‑9). The transformers’ parameters are explained in the video clip.

 

Figure 8‑8 A hydraulic lift system

The following video shows how to build and run the model for this example in 20-sim.

 

The BG model is shown in Figure 8‑9, along with the detail of the transformers’ inputs and outputs.

Figure 8-9 Bond graph model for hydraulic lift system

 

Exercise Problems for Chapter 8

Exercises

  1. Build the bond graph for a two-way safety valve.
  2. Repeat the example in section 8-9 and perform a parametric sweep for some parameters in the simulation, for example pipe diameters and lengths.
  3. Expand the BG model given in section 8.10 with running simulation with some data for the system parameters, similar to those given in section 8.9. Also, expand the model of the pump using some pump-chart (H-Q).
  4. Use some data and run simulation for the example given in section 8.11, the hydraulic lift.

Media Attributions

License

Icon for the Creative Commons Attribution-NonCommercial 4.0 International License

Engineering Systems Dynamics, Modelling, Simulation, and Design Copyright © 2021 by Mehrzad Tabatabaian is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License, except where otherwise noted.

Share This Book