Bond graph

Short description: Graphical representation of a dynamic system

A simple mass–spring–damper system, and its equivalent bond-graph form

A bond graph is a graphical representation of a physical dynamic system. It allows the conversion of the system into a state-space representation. It is similar to a block diagram or signal-flow graph, with the major difference that the arcs in bond graphs represent bi-directional exchange of physical energy, while those in block diagrams and signal-flow graphs represent uni-directional flow of information. Bond graphs are multi-energy domain (e.g. mechanical, electrical, hydraulic, etc.) and domain neutral. This means a bond graph can incorporate multiple domains seamlessly.

The bond graph is composed of the "bonds" which link together "single-port", "double-port" and "multi-port" elements (see below for details). Each bond represents the instantaneous flow of energy (dE/dt) or power. The flow in each bond is denoted by a pair of variables called power variables, akin to conjugate variables, whose product is the instantaneous power of the bond. The power variables are broken into two parts: flow and effort. For example, for the bond of an electrical system, the flow is the current, while the effort is the voltage. By multiplying current and voltage in this example you can get the instantaneous power of the bond.

A bond has two other features described briefly here, and discussed in more detail below. One is the "half-arrow" sign convention. This defines the assumed direction of positive energy flow. As with electrical circuit diagrams and free-body diagrams, the choice of positive direction is arbitrary, with the caveat that the analyst must be consistent throughout with the chosen definition. The other feature is the "causality". This is a vertical bar placed on only one end of the bond. It is not arbitrary. As described below, there are rules for assigning the proper causality to a given port, and rules for the precedence among ports. Causality explains the mathematical relationship between effort and flow. The positions of the causalities show which of the power variables are dependent and which are independent.

If the dynamics of the physical system to be modeled operate on widely varying time scales, fast continuous-time behaviors can be modeled as instantaneous phenomena by using a hybrid bond graph. Bond graphs were invented by Henry Paynter.

Systems for bond graph

Many systems can be expressed in terms used in bond graph. These terms are expressed in the table below.

Conventions for the table below:

• $P$ is the active power;
• $\hat{X}$ is a matrix object;
• $\overrightarrow{x}$ is a vector object;
• $x^{†}$ is the Hermitian of x; it is the complex conjugate of the transpose of x. If x is a scalar, then the Hermitian is the same as the complex conjugate;
• $D_t^n$ is the Euler notation for differentiation, where: $$D_t^{n}f(t)=\{\begin{array}{ll}\underset{-\mathrm{\infty }}{\overset{t}{\int }}f(s)ds,& n<0\\ f(t),& n=0\\ \frac{{\partial }^{n}f(t)}{\partial t^n},& n>0\end{array}$$
• $\{\begin{array}{l}⟨x{⟩}^{\alpha }:=|x{|}^{\alpha }\mathrm{sgn}(x)\\ ⟨a⟩=k⟨b{⟩}^{\beta }⟹⟨b⟩={(\frac{1}{k}⟨a⟩)}^{1/\beta }\end{array}$
• Vergent-factor: $${\varphi }_L=\{\begin{array}{l}\text{<mrow data-mjx-texclass="ORD">Prismatic</mrow>}:\frac{\text{<mrow data-mjx-texclass="ORD">length</mrow>}}{\text{<mrow data-mjx-texclass="ORD">cross<mo stretchy="false">−</mo>sectional</mrow>}\text{<mrow data-mjx-texclass="ORD">area</mrow>}}\\ \text{<mrow data-mjx-texclass="ORD">Cylinder</mrow>}:\frac{\ln \left(\frac{\mathrm{r}\mathrm{a}\mathrm{d}\mathrm{i}\mathrm{u}{\mathrm{s}}_{\mathrm{o}\mathrm{u}\mathrm{t}}}{\mathrm{r}\mathrm{a}\mathrm{d}\mathrm{i}\mathrm{u}{\mathrm{s}}_{\mathrm{i}\mathrm{n}}}\right)}{2\pi \cdot \text{<mrow data-mjx-texclass="ORD">length</mrow>}}\\ \text{<mrow data-mjx-texclass="ORD">Sphere</mrow>}:\frac{1}{4\pi (\mathrm{r}\mathrm{a}\mathrm{d}\mathrm{i}\mathrm{u}{\mathrm{s}}_{\mathrm{i}\mathrm{n}}\parallel \mathrm{-}\mathrm{r}\mathrm{a}\mathrm{d}\mathrm{i}\mathrm{u}{\mathrm{s}}_{\mathrm{o}\mathrm{u}\mathrm{t}})}\end{array}$$

	Generalized flow	Generalized displacement	Generalized effort	Generalized momentum	Generalized power (in watt for power systems)	Generalized energy (in joule for power systems)
Name	$\overrightarrow{f}(t)$	$\overrightarrow{q}(t)$	$\overrightarrow{e}(t)$	$\overrightarrow{p}(t)$	$P=\overrightarrow{f}(t{)}^{†}\overrightarrow{e}(t)$	$E=\overrightarrow{q}(t{)}^{†}\overrightarrow{e}(t)$
Description	Time derivative of displacement	A quality related to static behaviour.	The energy per unit of displacement	Time integral of effort	Transformation of energy from one to another form	Conserved quantity in closed systems
Elements
Name	Hyperance $H$ Hyperrigitance $P=H^{-1}$	Compliance $C$ Rigitance $K=C^{-1}$	Resistance $R$	Inertance $I$ (or $L$)	Abrahance $A$	Magnance $M$
Properties	Power dissipative element	Charge storage element (State variable: displacement) (Costate variable: effort)	Power dissipative element	Momentum storage element (State variable: momentum) (Costate variable: flow)	Power dissipative element	Power dissipative element
Quantitative behaviour	For 1-dimension systems (linear): $$P=H\cdot {(D_t^{0}q(t))}^2$$ For 1-dimension systems: $$e(t)=P\gamma [D_t^{-1}q(t)]$$ Impedance: $$Z(s)=\frac{1}{s^{2}H}=\frac{1}{s^2}P$$	Potential energy for N-dimension systems: $$\begin{array}{lcl}V& =& \frac{1}{2}\overrightarrow{q}(t{)}^{†}\overrightarrow{e}(t)\\ \overrightarrow{q}(t)& =& \hat{C}\overrightarrow{e}(t)\end{array}$$ Potential energy: $$V=\underset{0}{\overset{q}{\int }}e(q)dq$$ Potential coenergy: $$\bar{V}=\underset{0}{\overset{e}{\int }}q(e)de$$ For 1-dimension systems: $$f(t)=C\cdot \frac{de}{dt}+e\frac{dC}{dt}$$ Impedance: $$Z(s)=\frac{1}{sC}=\frac{1}{s}k$$	For 1-dimension systems (linear): $$P=R\cdot {(D_{t}q(t))}^2$$ Power for 1-dimension non-linear resistances ($e(f)$ is the effort developed by the element): $$P=e(f)f$$ Rayleigh Power: $$\mathfrak{\Re }=\frac{1}{2}R\cdot f(t{)}^2$$ Rayleigh Power for non-linear resistances: $$\mathfrak{\Re }=\underset{0}{\overset{f}{\int }}e(f)df$$ Rayleigh Effort: $$e_{\mathfrak{\Re }}=\frac{d\mathfrak{\Re }}{df}=e(f)$$ For N-dimension systems: $$\begin{array}{lcr}P& =& \overrightarrow{f}(t{)}^{†}\overrightarrow{e}(t)\\ \overrightarrow{e}(t)& =& \hat{R}\overrightarrow{f}(t)\end{array}$$ For 1-dimension systems: $$e(t)=R\cdot \gamma [D_t^{1}q(t)]$$ Impedance: $$Z(s)=R$$	Kinetic energy for N-dimension systems: $$\begin{array}{l}T=\frac{1}{2}\overrightarrow{\rho }(t{)}^{†}\overrightarrow{f}(t)\\ \overrightarrow{\rho }(t)=\hat{L}\overrightarrow{f}(t)\end{array}$$ Kinetic energy: $$T=\underset{0}{\overset{\rho }{\int }}f(\rho )d\rho$$ Kinetic coenergy: $$\bar{T}=\underset{0}{\overset{f}{\int }}\rho (f)df$$ For 1-dimension systems: $$e(t)=L\cdot \frac{df}{dt}+f\cdot \frac{dL}{dt}$$ Impedance: $$Z(s)=sL$$	For 1-dimension systems (linear): $$P=A\cdot {(D_t^{2}q(t))}^2$$ For 1-dimension systems: $$e(t)=A\cdot \gamma [D_t^{3}q(t)]$$ Impedance $$Z(s)=s^{2}A$$	For 1-dimension systems (linear) $$P=A\cdot {(D_t^{3}q(t))}^2$$ For 1-dimension systems $$e(t)=M\cdot \gamma [D_t^{5}q(t)]$$ Impedance $$Z(s)=s^{4}M$$
Generalized behaviour	Energy from active effort sources: $W=\underset{0}{\overset{q}{\int }}{e}_{\text{source}}dq$ Lagrangian: $\mathfrak{𝔏}=T-(V-W)$ Hamiltonian: $\mathfrak{ℌ}=T+(V-W)$ Hamiltonian Effort: $e_{\mathfrak{ℌ}}=\frac{d\mathfrak{ℌ}}{dq}$ Lagrangian Effort: $e_{\mathfrak{𝔏}}=\frac{d\mathfrak{𝔏}}{dq}$ Passive Effort: $e_{\mathfrak{𝔏}\mathfrak{ℌ}}=\frac{d}{dt}\frac{d\mathfrak{𝔏}}{df}$ Power Equation: $\frac{d\mathfrak{𝔏}}{dt}+\frac{d\mathfrak{ℌ}}{dt}=0$ Effort Equation: $e_{\mathfrak{𝔏}}+e_{\mathfrak{ℌ}}=0$ Lagrangian Equation: $e_{\mathfrak{\Re }}+e_{\mathfrak{𝔏}\mathfrak{ℌ}}=e_{\mathfrak{𝔏}}$ Hamiltonian Equation: $e_{\mathfrak{\Re }}+e_{\mathfrak{𝔏}\mathfrak{ℌ}}=-e_{\mathfrak{ℌ}}$ If $\bar{W}$ is the coenergy, $W$ is the energy, $S$ is the state variable and $\bar{S}$ is the costate variable, $\begin{array}{r}\bar{W}+W=S\cdot \bar{S}\\ \bar{W}=\underset{0}{\overset{\bar{S}}{\int }}S(\bar{S})d\bar{S}\\ W=\underset{0}{\overset{S}{\int }}\bar{S}(S)dS\end{array}$ For linear elements: $\bar{W}=W=\frac{1}{2}S\cdot \bar{S}$

Longitudinal mechanical power system

Flow-related variables	$D_t^{6}x$: Pounce / Pop $\mathrm{[}\mathrm{m}/{\mathrm{s}}^{\mathrm{6}}\mathrm{]}$	$D_t^{5}x$: Flounce / Crackle $\mathrm{[}\mathrm{m}/{\mathrm{s}}^{\mathrm{5}}\mathrm{]}$	$D_t^{4}x$: Jounce / Snap $\mathrm{[}\mathrm{m}/{\mathrm{s}}^{\mathrm{4}}\mathrm{]}$	$D_t^{3}x$: Jerk $\mathrm{[}\mathrm{m}/{\mathrm{s}}^{\mathrm{3}}\mathrm{]}$
	$D_t^{2}x$: Acceleration $(x_{2t})\mathrm{[}\mathrm{m}/{\mathrm{s}}^{\mathrm{2}}\mathrm{]}$	$D_t^{1}x$: Velocity $(x_t)\mathrm{[}\mathrm{m}/\mathrm{s}\mathrm{]}$ (Flow)	$D_t^{0}x$: Displacement $(x)\mathrm{[}\mathrm{m}\mathrm{]}$ (Displacement)	$D_t^{-1}x$: Absement $\mathrm{[}\mathrm{m}\mathrm{\cdot }\mathrm{s}\mathrm{]}$
	$D_t^{-2}x$: Absity $\mathrm{[}\mathrm{m}\mathrm{\cdot }{\mathrm{s}}^{\mathrm{2}}\mathrm{]}$	$D_t^{-3}x$: Abseleration $\mathrm{[}\mathrm{m}\mathrm{\cdot }{\mathrm{s}}^{\mathrm{3}}\mathrm{]}$	$D_t^{-4}x$: Abserk $\mathrm{[}\mathrm{m}\mathrm{\cdot }{\mathrm{s}}^{\mathrm{4}}\mathrm{]}$
Effort-related variables	$D_t^{1}F$: Yank $\mathrm{[}\mathrm{N}/\mathrm{s}\mathrm{]}$	$D_t^{0}F$: Force $(P_x^t)\mathrm{[}\mathrm{N}\mathrm{]}$ (Effort)	$D_t^{-1}F$: Linear momentum $(P_x^{2t})\mathrm{[}\mathrm{k}\mathrm{g}\mathrm{\cdot }\mathrm{m}/\mathrm{s}\mathrm{]}$ (Momentum)
Passive elements
Compliance (C)	Resistance (R)	Inertance (I)	Abrahance (A)	Magnance (M)
Spring $$⟨P_x^t⟩=k⟨x⟩$$ where $k$ is the spring Stiffness	Damper $$⟨P_x^t⟩=b⟨x_t⟩$$ where $b$ is the damper parameter	Mass $$⟨P_x^t⟩=m⟨x_{2t}⟩$$ where $m$ is the mass	Abraham–Lorentz force $$⟨P_x^t⟩=\frac{{\mu }_0{q}^2}{6\pi c}⟨x_{3t}⟩$$ where ${\mu }_0$ is the permeability $q$ is the electric charge $c$ is the speed of light	Magnetic radiation reaction force $$⟨P_q^t⟩=\frac{{\mu }_0{q}^{2}R}{24\pi c^3}⟨x_{5t}⟩$$ where ${\mu }_0$ is the permeability $q$ is the electric charge $c$ is the speed of light $R$ is the radius of the magnetic moment
Cantilever $$⟨P_x^t⟩=\frac{3EI}{L^3}⟨x⟩$$ $L$: Cantilever length $E$: Young's modulus $I$: Second moment of area	Cyclotron radiation resistance $$⟨P_x^t⟩=\frac{{\sigma }_t{B}^2}{c{\mu }_0}⟨x_t⟩$$ where ${\sigma }_t$: Thompson cross-section $c$: Speed of light ${\mu }_0$: Permeability $B$: Magnetic field density
Prismatic floater in a wide waterbody $$⟨P_x^t⟩={\rho }_{L}gA⟨x⟩$$ where ${\rho }_L$: Liquid density $A$: Area $g$: Gravitational acceleration	Viscous friction $$⟨P_x^t⟩=b⟨x_t⟩$$ where $b$ is the viscous friction paramete
Elastic Rod $$⟨P_x^t⟩=\frac{EA}{L}⟨x⟩$$ where $E$: Young's modulus $A$: Area $L$: Rod length	Inverse of kinetic mobility $$⟨P_x^t⟩=\frac{1}{\mu }⟨x_t⟩$$ where $\mu$ is the kinetic mobility
Newton's law of gravitation $$⟨P_x^t⟩=GMm⟨x{⟩}^{-2}$$ where $G$: Gravitational constant $M$: Mass of body 1 $m$: Mass of body 2	Geometry interacting with air (e.g. drag) $$⟨P_x^t⟩=\frac{1}{2}c\rho A⟨x_t{⟩}^2$$ where $A$: Contact area $c$: Aerodynamic geometry coefficient $\rho$: Fluid density
Coulomb's law $$⟨P_x^t⟩=\frac{Qq}{4\pi {\epsilon }_0}⟨x{⟩}^{-2}$$ where ${\epsilon }_0$: Permittivity $Q$: Charge of body 1 $q$: Charge of body 2	Absquare damper $$⟨P_x^t⟩=B⟨x_t{⟩}^2$$ where $B$ is the Absquare damper parameter
Casimir force $$⟨P_x^t⟩=\frac{A\hslash c{\pi }^2}{240}⟨x{⟩}^{-4}$$ where $A$: Area of plate $\hslash$: Planck's reduced constant $c$: Speed of light	Dry friction $$⟨P_x^t⟩=\mu F_n⟨x_t{⟩}^0$$ where $\mu$: Kinetic friction coefficient $F_n$: Normal force
Biot-Savart Law $$⟨P_x^t⟩=\frac{{\mu }_0{I}_1{I}_{2}l}{2\pi }⟨x{⟩}^{-1}$$ where ${\mu }_0$: Permeability $I_1$: Current in wire 1 $I_2$: Current in wire 2 $l$: Length of wires
Piston pressing fluid inside an adiabatic chamber $$P_x^t=A{P}_0{(1+\frac{x}{x_0})}^{-\gamma }$$ where $A$: Area of piston $P_0$: Initial pressure inside $x_0$: Initial position of piston $\gamma$: Poisson coefficient

Angular mechanical power system

Flow-related variables	$D_t^3\theta$: Angular jerk $\mathrm{[}\mathrm{r}\mathrm{a}\mathrm{d}/{\mathrm{s}}^{\mathrm{3}}\mathrm{]}$	$D_t^2\theta$: Angular acceleration $({\theta }_{2t})\mathrm{[}\mathrm{r}\mathrm{a}\mathrm{d}/{\mathrm{s}}^{\mathrm{2}}\mathrm{]}$
	$D_t^1\theta$: Angular velocity $({\theta }_t)\mathrm{[}\mathrm{r}\mathrm{a}\mathrm{d}/\mathrm{s}\mathrm{]}$ (Flow)	$D_t^0\theta$: Angular displacement $(\theta )\mathrm{[}\mathrm{r}\mathrm{a}\mathrm{d}\mathrm{]}$ (Displacement)
Effort-related variables	$D_t^1\tau$: Rotatum $\mathrm{[}\mathrm{W}/\mathrm{r}\mathrm{a}\mathrm{d}\mathrm{]}$	$D_t^0\tau$: Torque $(P_{\theta }^t)\mathrm{[}\mathrm{J}/\mathrm{r}\mathrm{a}\mathrm{d}\mathrm{]}$ (Effort)
	$D_t^{-1}\tau$: Angular momentum $(P_{\theta }^{2t})\mathrm{[}\mathrm{J}\mathrm{\cdot }\mathrm{s}/\mathrm{r}\mathrm{a}\mathrm{d}\mathrm{]}$ (Momentum)
Passive elements
Compliance (C)	Resistance (R)	Inertance (I)
Inverse of the angular spring constant $$⟨P_{\theta }^t⟩=k_r⟨\theta ⟩$$ where $k_r$ is the angular spring constant	Angular damping $$⟨P_{\theta }^t⟩=R⟨{\theta }_t⟩$$ where $R$ is the damping constant	Mass moment of inertia Type $⟨P_{\theta }^t⟩=I⟨{\theta }_{2t}⟩$ where $I$ is the mass moment of inertia
Rod Torsion $$⟨P_{\theta }^t⟩=\frac{GJ}{L}⟨\theta ⟩$$ where $G$: Shear modulus $J$: Polar moment of area $L$: Rod length	Governor (e.g. used in music boxes) $$⟨P_{\theta }^t⟩=R⟨{\theta }_t{⟩}^2$$ where $R$ is the Governor constant
Bending moment (cantilever) $$⟨P_{\theta }^t⟩=\frac{EI}{L}⟨\theta ⟩$$ where $E$: Young modulus $I$: Second moment of area $L$: Beam length
Parallel force field $$⟨P_{\theta }^t⟩=⟨P_x^t⟩\sin (\alpha -\theta )$$ where $\alpha$: Angle (counter-clockwise positive) between the polar axis and the force field $\theta$: Current angle of the object (e.g. pendulum)

Electric power system

Flow-related variables	$D_t^{2}q$: Electric Inertia $\mathrm{[}\mathrm{A}/\mathrm{s}\mathrm{]}$	$D_t^{1}q$: Electric current $(q_t)\mathrm{[}\mathrm{A}\mathrm{]}$ (Flow)	$D_t^{0}q$: Electric charge $(q)\mathrm{[}\mathrm{C}\mathrm{]}$ (Displacement)
Effort-related variables	$D_t^{1}V$: Distension $\mathrm{[}\mathrm{V}/\mathrm{s}\mathrm{]}$	$D_t^{0}V$: Voltage $(P_q^t)\mathrm{[}\mathrm{V}\mathrm{]}$ (Effort)	$D_t^{-1}V$: Flux linkage $(P_q^{2t})\mathrm{[}\mathrm{V}\mathrm{\cdot }\mathrm{s}\mathrm{]}\text{ or }\mathrm{[}\mathrm{W}\mathrm{b}\mathrm{\cdot }\mathrm{t}\mathrm{u}\mathrm{r}\mathrm{n}\mathrm{]}$ (Momentum)
Elements
Hyperance (H)	Compliance (C)	Resistance (R)	Inertance (I)	Abrahance (A)
Frequency Dependent Negative Resistor (FDNR) $$⟨P_q^t⟩=\frac{1}{H}⟨q^t⟩$$	Linear capacitor $$⟨P_q^t⟩=\frac{1}{\epsilon }{\varphi }_L⟨q⟩$$ where $\epsilon$: Permittivity ${\varphi }_L$: Vergent-factor	Linear resistor $$⟨P_q^t⟩=\rho {\varphi }_L(1+\alpha (T-T_0))⟨q_t⟩$$ where $\rho$: Resistivity ${\varphi }_L$: Vergent-factor $\alpha$: Thermal coefficient $T$: Current temperature $T_0$: Reference Temperature	Linear inductor (solenoid) $$⟨P_q^t⟩=\frac{{\mu }_0{N}^{2}A}{L}⟨q_{2t}⟩$$ where ${\mu }_0$: Permeability $N$: Number of turns $A$: Area $L$: Length	Frequency Dependent Negative Conductance (FDNC) Type $${\gamma }^1⟺V=R{D}_t^{2}i[H\cdot s]$$
		Diode $$P_q^t=n{V}_T\ln (1+\frac{q_t}{I_s})$$ where $n$: Ideality factor $V_T$: Thermal voltage $I_s$: Leakage current	Toroid $$⟨P_q^t⟩=\frac{\mu N^{2}A}{2\pi r}⟨q_{2t}⟩$$ where $\mu$: Permeability $N$: Number of turns $A$: Cross-sectional area $r$: Toroid radius to centerline
Intra-Gyrator			Inter-Gyrator
Compliant Gyrator	Resistive Gyrator	Inertant Gyrator
	Hall effect device $$\begin{array}{r}{e}_x=R_H\frac{B_z}{t_z}{i}_y\\ e_y=R_H\frac{B_z}{t_z}{i}_x\end{array}$$ where $R_H$: Hall coefficient $B_z$: Vertical magnetic induction field $t_z$: Vertical thickness	Induction motor $$\begin{array}{r}{e}_r^1=M\cos (\theta )D_t{f}_s^1+f_s^1{D}_t[M\cos (\theta )]\\ e_s^1=M\cos (\theta )D_t{f}_r^1+f_r^1{D}_t[M\cos (\theta )]\end{array}$$ where $\theta$: Electrical angle $e_r^1$: Voltage at rotor bar $e_s^1$: Voltage at stator bar $f_s^1$: Flow through stator bar $f_r^1$: Flow through rotor bar	DC motor $$\begin{array}{r}{\tau }_e=k_a{\varphi }_a(i_e)i_a\\ {\omega }_e=\frac{1}{k_a{\varphi }_a(i_e)}{e}_a\end{array}$$ where ${\tau }_e$: Electromagnetic torque ${\omega }_e$: Axis angular velocity $i_e$: Field current $i_a$: Armature current	Faraday gyrator $$\begin{array}{rl}F& =Bli\\ V& =\frac{1}{Bl}v\end{array}$$ where $F$: Force $v$: Rod velocity $V$: Rod voltage $B$: Magnetic induction field $l$: Rod length
			Faraday Disk $$\begin{array}{rl}V& =\frac{1}{2}B{r}^2\omega \\ \tau & =\frac{1}{2}B{r}^{2}i\end{array}$$ where $B$: Magnetic induction field $r$: Disk radius
Intra-Transformer			Inter-Transformer
Electrical Transformer (Only for AC signals) $$\begin{array}{rl}{V}_2& =\frac{N_2}{N_1}{V}_1\\ f_2& =\frac{N_1}{N_2}{f}_1\end{array}$$

Hydraulic / pneumatic power system

Flow-related variables	$D_t^{1}V$: Volumetric flow rate $(V_t)\mathrm{[}{\mathrm{m}}^{\mathrm{3}}/\mathrm{s}\mathrm{]}$ (Flow)	$D_t^{0}V$: Volume $(V)\mathrm{[}{\mathrm{m}}^{\mathrm{3}}\mathrm{]}$ (Displacement)
Effort-related variables	$D_t^{0}P$: Pressure $(P_V^t)\mathrm{[}\mathrm{P}\mathrm{a}\mathrm{]}$ (Effort)	$D_t^{-1}P$: Fluid momentum $(P_V^{2t})\mathrm{[}\mathrm{P}\mathrm{a}\mathrm{\cdot }\mathrm{s}\mathrm{]}$ (Momentum)
Elements
Compliance (C)	Resistance (R)	Inertance (I)
Pipe elasticity $$⟨P_V^t⟩=\left(\frac{t_{W}E}{2r_0{V}_0}\right)⟨V⟩$$ where $r_0$: Nominal pipe radius $V_0$: Pipe volume (without stress) $t_W$: Wall thickness $E$: Young's modulus	Darcy sponge $$⟨P_V^t⟩=\left(\frac{\mu }{k}{\varphi }_L\right)⟨V_t⟩$$ where $\mu$: Dynamic viscosity $k$: Permeability ${\varphi }_L$: Vergent-factor	Fluid inertia in pipes $$⟨P_V^t⟩=\left(\rho {\varphi }_L\right)⟨V_{2t}⟩$$ where $\rho$: Fluid density ${\varphi }_L$: Vergent-factor
Compressible fluid (approximation) $$⟨P_V^t⟩=\left(\frac{B}{V_0}\right)⟨V⟩=\underset{\begin{array}{c}\text{Acoustic}\\ \text{Approximation}\end{array}}{\underbrace{\left(\frac{{\rho }_0{c}^2}{V_0}\right)}}⟨V⟩$$ where $V_0$: Pipe Volume (without stress) $B$: Bulk Modulus ${\rho }_0$: Gas density at reference pressure $c$: Speed of sound	Valve $$⟨P_V^t⟩=\left(\frac{\rho }{2C_d^2{A}_0^2}\right)⟨V_t{⟩}^2$$ $C_d$: Discharge coefficient $\rho$: Fluid density $A_0$: Smallest area for fluid passage
Tank with area $A\frac{1}{[m{]}^{n-1}}{(h+h_0)}^{n-1}$: $$P_V^t=\rho g{h}_0[{(1+\frac{nV}{A\frac{1}{[m{]}^{n-1}}{h}_0^n})}^{1/n}-1]$$ where $h$: Height with respect to the ground $h_0$: Inlet height $\rho$: Fluid density $A$: Curve scaling (dimension of area) $[m{]}^{n-1}$: Correction of dimension $n$: Curve parameter $g$: Acceleration of gravity Case $n=1$ (prismatic case): $$P_V^t=\frac{\rho g}{A}V$$ Case $n=0$ (current diode case): $$P_V^t=\rho g{h}_0(e^{\frac{V}{A[m]}}-1)$$	Poiseuille resistance for cylinders $$⟨P_V^t⟩=\left(\frac{8\mu L}{\pi R^4}\right)⟨V_t⟩$$ where $\mu$: Dynamic viscosity $L$: Pipe length $R$: Pipe radius
Isothermal chamber $$⟨P_V^t⟩=k⟨V_t{⟩}^{-1}$$ where $k$ is Chamber's constant	Turbulence resistance $$⟨P_V^t⟩=a_t⟨V_t{⟩}^{\frac{7}{4}}$$ where $a_t$ is an empirical parameter
Compressible fluid $$P_V^t=-B\ln V$$ where $B$ is the bulk modulus	Nozzle $$⟨P_V^t⟩=\frac{\rho }{2(A_{\text{out}}^2\parallel -A_{\text{in}}^2)}⟨V_t{⟩}^2$$ where $\rho$: Fluid density $A_{\text{out}}$: Outlet area $A_{\text{in}}$: Inlet area
Adiabatic bladder $$P_V^t=P_{V,0}^t{(1-\frac{V}{V_0})}^{-\gamma }$$ where $P_{V,0}^t$: Reference pressure $V_0$: Reference volume $\gamma$: Poisson coefficient	Check valve $$P_V^t=k\ln (1+\frac{V_t}{V_{t,0}})$$ where $k$: Empirical constant $V_{t,0}$: Reverse bias absolute flow

Gyrator-capacitor power system

Flow-related variables	$D_t^0\phi$: Magnetic Flux $(\phi )\mathrm{[}\mathrm{W}\mathrm{b}\mathrm{]}$ (Displacement)
Effort-related variables	$D_t^0{ℱ}$: Magnetomotive force $({ℱ})\mathrm{[}\mathrm{A}\mathrm{\cdot }\mathrm{t}\mathrm{u}\mathrm{r}\mathrm{n}\mathrm{]}$ (Momentum)
Elements
Compliance (C)	Resistance (R)	Inertance (I)
Permeance (${𝒫}$) $$⟨{ℱ}⟩=\frac{1}{\mu }{\varphi }_L⟨\phi ⟩\mathrm{[}\mathrm{H}/\mathrm{t}\mathrm{u}\mathrm{r}{\mathrm{n}}^{\mathrm{2}}\mathrm{]}$$ where $\mu$: Permeability ${\varphi }_L$: Vergent-factor	Magnetic Complex Impedance ($Z_M$) $${ℱ}=Z_M{\phi }_t\mathrm{[}{\mathrm{t}\mathrm{u}\mathrm{r}\mathrm{n}}^{\mathrm{2}}/\mathrm{\Omega }\mathrm{]}$$	Magnetic Complex Inductance ($L_M$) $${ℱ}=L_M{\phi }_{2t}\mathrm{[}{\mathrm{t}\mathrm{u}\mathrm{r}\mathrm{n}}^{\mathrm{2}}\mathrm{F}\mathrm{]}$$

Gravitational power system

Flow-related variables	$D_t^0{i}_g$: Gravitational Current Loop $(i_g=2{\epsilon }_g{v}_{\text{orbit}}^3)\mathrm{[}\mathrm{k}\mathrm{g}/\mathrm{s}\mathrm{]}$ (Flow)	$D_t^{-1}{i}_g$: Gravitational Charge $(M)\mathrm{[}\mathrm{k}\mathrm{g}\mathrm{]}$ (Displacement)
Effort-related variables	$D_t^0{V}_g$: Gravitational Voltage $(V_g=\frac{1}{2}\frac{v_{\text{orbit}}^4}{c^2})\mathrm{[}{\mathrm{m}}^{\mathrm{2}}/{\mathrm{s}}^{\mathrm{2}}\mathrm{]}$ (Effort)	$D_t^{-1}{V}_g$: Gravitational Momentum $({\varphi }_g=\frac{\pi v_{\text{orbit}}^{3}r}{c^2})\mathrm{[}{\mathrm{m}}^{\mathrm{2}}/\mathrm{s}\mathrm{]}$ (Momentum)
Elements
Compliance (C)	Resistance (R)	Inertance (I)
Gravitational Capacitance Type ${\gamma }^1⟺M_g=C{V}_g\mathrm{[}\mathrm{k}\mathrm{g}/{\mathrm{m}}^{\mathrm{2}}\mathrm{]}$ $C_g=\frac{2c^2{r}^2}{GM}$	Gravitational Orbital Resistance Type ${\gamma }^1⟺V_g=R_g{i}_g\mathrm{[}{\mathrm{m}}^{\mathrm{2}}/\mathrm{k}\mathrm{g}\mathrm{\cdot }\mathrm{s}\mathrm{]}$ $R_g=\frac{{\mu }_g}{4}{v}_{\text{orbit}}$	Gravitational Inductance Type ${\gamma }^1⟺{\varphi }_g=I{i}_g\mathrm{[}{\mathrm{m}}^{\mathrm{2}}/\mathrm{k}\mathrm{g}\mathrm{\cdot }{\mathrm{s}}^{\mathrm{2}}\mathrm{]}$ $L_g=\frac{2{\pi }^{2}G}{c^2}r$

Electric field volumetric power density system

Flow-related variables	$D_t^{0}J$: Current Density $(J)\mathrm{[}\mathrm{A}/{\mathrm{m}}^{\mathrm{2}}\mathrm{]}$ (Flow)	$D_t^{-1}J$: Electric Displacement Field $(D)\mathrm{[}\mathrm{C}/{\mathrm{m}}^{\mathrm{2}}\mathrm{]}$ (Displacement)
Effort-related variables	$D_t^{0}E$: Electric Field $(E)\mathrm{[}\mathrm{V}/\mathrm{m}\mathrm{]}$ (Effort)	$D_t^{-1}E$: Magnetic Potential Vector $(A)\mathrm{[}\mathrm{V}\mathrm{\cdot }\mathrm{s}/\mathrm{m}\mathrm{]}$ (Momentum)
Elements
Compliance (C)	Resistance (R)	Volumetric Power Density Inertance ($I_V$)
Electrical Permittivity Type ${\gamma }^1⟺D={ϵ}_{0}E\mathrm{[}\mathrm{F}/\mathrm{m}\mathrm{]}$	Electrical Resistivity Type ${\gamma }^1⟺E=\rho J\mathrm{[}\mathrm{\Omega }\mathrm{m}\mathrm{]}$	Magnetic Permeability ${\mu }_0\mathrm{[}\mathrm{H}/\mathrm{m}\mathrm{]}$

Magnetic field volumetric power density system

Flow-related variables	$D_t^{0}B$: Magnetic Flux Density $(B)\mathrm{[}\mathrm{T}\mathrm{]}\text{ or }\mathrm{[}\mathrm{W}\mathrm{b}/{\mathrm{m}}^{\mathrm{2}}\mathrm{]}$ (Displacement)
Effort-related variables	$D_t^{0}H$: Magnetic Field Strength $(H)\mathrm{[}\mathrm{A}\mathrm{\cdot }\mathrm{t}\mathrm{u}\mathrm{r}\mathrm{n}/\mathrm{m}\mathrm{]}$ (Effort)
Elements
Compliance (C)	Resistance (R)	Volumetric Power Density Inertance ($I_V$)
Magnetic Permeability for Magnetic Circuits Type ${\gamma }^1⟺B=\mu H\mathrm{[}\mathrm{H}/\mathrm{(}\mathrm{m}\mathrm{\cdot }{\mathrm{t}\mathrm{u}\mathrm{r}\mathrm{n}}^{\mathrm{2}}\mathrm{)}\mathrm{]}$

Gravitoelectric volumetric power density system

Flow-related variables	$D_t^0{J}_g$: Flux of Mass $(J_g)\mathrm{[}\mathrm{k}\mathrm{g}/{\mathrm{m}}^{\mathrm{2}}\mathrm{s}\mathrm{]}$ (Flow)	$D_t^{-1}{J}_g$: Accumulated Flux of Mass $(D_g)\mathrm{[}\mathrm{k}\mathrm{g}/{\mathrm{m}}^{\mathrm{2}}\mathrm{]}$ (Displacement)
Effort-related variables	$D_t^{0}g$: Acceleration of Gravity $(g)\mathrm{[}\mathrm{m}/{\mathrm{s}}^{\mathrm{2}}\mathrm{]}$ (Effort)
Elements
Compliance (C)	Resistance (R)	Volumetric Power Density Inertance ($I_V$)
Gravitational Permittivity Type ${\gamma }^1⟺D_g={ϵ}_{g}g\mathrm{[}\mathrm{k}\mathrm{g}\mathrm{\cdot }{\mathrm{s}}^{\mathrm{2}}/{\mathrm{m}}^{\mathrm{3}}\mathrm{]}$ ${\epsilon }_g=\frac{1}{4\pi G}$		Gravitational Permeability $\mathrm{[}\mathrm{m}/\mathrm{k}\mathrm{g}\mathrm{]}$$$I_V={\mu }_g=\frac{4\pi G}{c^2}$$

Gravitomagnetic volumetric power density system

Flow-related variables	$D_t^0{B}_g$: Gravitomagnetic Field $(B_g={\omega }_{\text{orbit}}{\left(\frac{v_{\text{orbit}}}{c}\right)}^2)\mathrm{[}\mathrm{H}\mathrm{z}\mathrm{]}$ (Displacement)
Effort-related variables	$D_t^0{H}_g$: Gravitomagnetic Field Strength $(H_g)\mathrm{[}\mathrm{P}\mathrm{a}\mathrm{\cdot }\mathrm{s}\mathrm{]}$ (Effort)
Elements
Compliance (C)	Resistance (R)	Volumetric Power Density Inertance ($I_V$)
Gravitational Permeability $-\frac{4\pi G}{c^2}\mathrm{[}\mathrm{m}/\mathrm{k}\mathrm{g}\mathrm{]}$

Thermal power-temperature system

Flow-related variables	$D_t^{1}Q$: Heat-rate $({\psi }_t)\mathrm{[}\mathrm{W}\mathrm{]}$ (Flow)	$D_t^{0}Q$: Total heat $(\psi )\mathrm{[}\mathrm{J}\mathrm{]}$ (Displacement)
Effort-related variables	$D_t^{0}T$: Temperature $(T)\mathrm{[}\mathrm{K}\mathrm{]}$ (Effort)
Elements
Compliance (C)	Resistance (R)	Inertance (I)
Isobaric heat $$T=\frac{1}{\rho V{c}_p}\psi =\frac{1}{C_P}\psi$$ where $\rho$: Object mass density $V$: Volume of object $c_p$: Pressure specific heat capacity	Conduction Resistance $$T=\frac{1}{k}{\varphi }_L{\psi }_t$$ where $k$: Thermal conductivity ${\varphi }_L$: Vergent-factor
Isocoric heat $$T=\frac{1}{C_V}\psi$$ where $C_V$ is the constant-volume heat capacitance	Convection Resistance $$T=\frac{1}{hA}{\psi }_t$$ where $h$: Convection coefficient $A$: Interface area
Isothermal heat $$T=\frac{1}{nR\ln \left(\frac{V_f}{V_i}\right)}\psi$$ where $R$: Universal gas constant $n$: Number of moles $V_f$: Final volume $V_i$: Initial volume	Stefan-Boltzmann law $$⟨T⟩={\left(\frac{1}{e\sigma A}\right)}^{0.25}⟨{\psi }_t{⟩}^{0.25}$$ where $e$: Emissivity $\sigma$: Stefan-Boltzmann constant $A$: Interface area

Continuum mechanics volumetric power density system

Flow-related variables	$D_t^1\epsilon$: Strain-rate $(\dot{\epsilon })\mathrm{[}\mathrm{H}\mathrm{z}\mathrm{]}$ (Flow)	$D_t^0\epsilon$: Strain $(\epsilon )\mathrm{[}\mathrm{1}\mathrm{]}$ (Displacement)
Effort-related variables	$D_t^0\sigma$: Stress $(\sigma )\mathrm{[}\mathrm{P}\mathrm{a}\mathrm{]}$ (Effort)
Elements
Compliance (C)	Resistance (R)	Volumetric Power Density Inertance ($I_V$)
Inverse of Rigidity Type ${\gamma }^1⟺\epsilon =C\sigma \mathrm{[}\mathrm{P}{\mathrm{a}}^{\mathrm{-}\mathrm{1}}\mathrm{]}$ $C=\frac{1}{K}$	Viscosity Type ${\gamma }^1⟺\sigma =R\dot{\epsilon }\mathrm{[}\mathrm{P}\mathrm{a}\mathrm{\cdot }\mathrm{s}\mathrm{]}$	Power Density Inertance: Density of Material $\rho \mathrm{[}\mathrm{k}\mathrm{g}/{\mathrm{m}}^{\mathrm{3}}\mathrm{]}$

Other systems:

• Thermodynamic power system (flow is entropy-rate and effort is temperature)
• Electrochemical power system (flow is chemical activity and effort is chemical potential)
• Thermochemical power system (flow is mass-rate and effort is mass specific enthalpy)
• Macroeconomics currency-rate system (displacement is commodity and effort is price per commodity)
• Microeconomics currency-rate system (displacement is population and effort is GDP per capita)

Tetrahedron of state

Tetrahedron of state

The tetrahedron of state is a tetrahedron that graphically shows the conversion between effort and flow. The adjacent image shows the tetrahedron in its generalized form. The tetrahedron can be modified depending on the energy domain.

Using the tetrahedron of state, one can find a mathematical relationship between any variables on the tetrahedron. This is done by following the arrows around the diagram and multiplying any constants along the way. For example, if you wanted to find the relationship between generalized flow and generalized displacement, you would start at the f(t) and then integrate it to get q(t). More examples of equations can be seen below.

Relationship between generalized displacement and generalized flow.

$$q(t)=\int f(t)dt$$

Relationship between generalized flow and generalized effort.

$$f(t)=\frac{1}{R}\cdot e(t)$$

Relationship between generalized flow and generalized momentum.

$$f(t)=\frac{1}{I}\cdot p(t)$$

Relationship between generalized momentum and generalized effort.

$$p(t)=\int e(t)dt$$

Relationship between generalized flow and generalized effort, involving the constant C.

$$e(t)=\frac{1}{C}\int f(t)dt$$

All of the mathematical relationships remain the same when switching energy domains, only the symbols change. This can be seen with the following examples.

Relationship between displacement and velocity.

$$x(t)=\int v(t)dt$$

Relationship between current and voltage, this is also known as Ohm's law.

$$i(t)=\frac{1}{R}V(t)$$

Relationship between force and displacement, also known as Hooke's law. The negative sign is dropped in this equation because the sign is factored into the way the arrow is pointing in the bond graph.

$$F(t)=kx(t)$$

For power systems, the formula for the frequency of resonance is as follows: $$\omega =\sqrt{\frac{1}{LC}}$$

For power density systems, the formula for the velocity of the resonance wave is as follows: $$c=\sqrt{\frac{1}{LC}}$$

Components

If an engine is connected to a wheel through a shaft, the power is being transmitted in the rotational mechanical domain, meaning the effort and the flow are torque (τ) and angular velocity (ω) respectively. A word bond graph is a first step towards a bond graph, in which words define the components. As a word bond graph, this system would look like: $$\text{engine}\stackrel{\tau }{\underset{\omega }{-----}}\text{wheel}$$ A half-arrow is used to provide a sign convention, so if the engine is doing work when τ and ω are positive, then the diagram would be drawn: $$\text{engine}\stackrel{\tau }{\underset{\omega }{---\rightharpoondown }}\text{wheel}$$ This system can also be represented in a more general method. This involves changing from using the words, to symbols representing the same items. These symbols are based on the generalized form, as explained above. As the engine is applying a torque to the wheel, it will be represented as a source of effort for the system. The wheel can be presented by an impedance on the system. Further, the torque and angular velocity symbols are dropped and replaced with the generalized symbols for effort and flow. While not necessary in the example, it is common to number the bonds, to keep track of in equations. The simplified diagram can be seen below.

$$S_e\stackrel{e_1}{\underset{f_1}{---\rightharpoondown }}\text{I}$$

Given that effort is always above the flow on the bond, it is also possible to drop the effort and flow symbols altogether, without losing any relevant information. However, the bond number should not be dropped. The example can be seen below.

$$S_e\stackrel{{}_1}{\underset{}{---\rightharpoondown }}\text{I}$$

The bond number will be important later when converting from the bond graph to state-space equations.

Association of elements

Series association

Suppose that an element has the following behavior: $$e(t)=\alpha g(q(t))$$ where $g(x)$ is a generic function (it can even differentiate/integrate its input) and $\alpha$ is the element's constant. Then, suppose that in a 1-junction you have many of this type of element. Then the total voltage across the junction is: $$e(t)=\left(\sum _i{\alpha }_i\right)g(q(t))⟹\begin{array}{c}{\alpha }_{\text{eq}}=\sum _{i=1}^N{\alpha }_i\end{array}$$

Parallel association

Suppose that an element has the following behavior: $$e(t)=g(\alpha q(t))$$ where $g(x)$ is a generic function (it can even differentiate/integrate its input) and $\alpha$ is the element's constant. Then, suppose that in a 0-junction you have many of this type of element. Then it is valid:

$$g^{-1}(e(t))={\alpha }_i{q}_i(t)⟹\frac{1}{{\alpha }_i}{g}^{-1}(e(t))=q_i(t)⟹\left(\sum _i\frac{1}{{\alpha }_i}\right)g^{-1}(e(t))=q(t)⟹g(g^{-1}(e(t)))=g(\frac{1}{\sum _i\frac{1}{{\alpha }_i}}q(t))⟹\begin{array}{c}{\alpha }_{\text{eq}}=\underset{i=1}{\overset{N}{\parallel }}{\alpha }_i\end{array}$$

Single-port elements

Single-port elements are elements in a bond graph that can have only one port.

Sources and sinks

Sources are elements that represent the input for a system. They will either input effort or flow into a system. They are denoted by a capital "S" with either a lower case "e" or "f" for effort or flow respectively. Sources will always have the arrow pointing away from the element. Examples of sources include: motors (source of effort, torque), voltage sources (source of effort), and current sources (source of flow).

$$S_e\stackrel{}{\underset{}{---\rightharpoonup }}J\text{and}{S}_f\stackrel{}{\underset{}{---\rightharpoonup }}J$$ where J indicates a junction.

Sinks are elements that represent the output for a system. They are represented the same way as sources, but have the arrow pointing into the element instead of away from it.

$$J\stackrel{}{\underset{}{---\rightharpoonup }}{S}_e\text{and}J\stackrel{}{\underset{}{---\rightharpoonup }}{S}_f$$

Inertia

Inertia elements are denoted by a capital "I", and always have power flowing into them. Inertia elements are elements that store energy. Most commonly these are a mass for mechanical systems, and inductors for electrical systems.

$$J\stackrel{}{\underset{}{---\rightharpoonup }}I$$

Resistance

Resistance elements are denoted by a capital "R", and always have power flowing into them. Resistance elements are elements that dissipate energy. Most commonly these are a damper, for mechanical systems, and resistors for electrical systems.

$$J\stackrel{}{\underset{}{---\rightharpoonup }}R$$

Compliance

Compliance elements are denoted by a capital "C", and always have power flowing into them. Compliance elements are elements that store potential energy. Most commonly these are springs for mechanical systems, and capacitors for electrical systems.

$$J\stackrel{}{\underset{}{---\rightharpoonup }}C$$

Two-port elements

These elements have two ports. They are used to change the power between or within a system. When converting from one to the other, no power is lost during the transfer. The elements have a constant that will be given with it. The constant is called a transformer constant or gyrator constant depending on which element is being used. These constants will commonly be displayed as a ratio below the element.

Transformer

A transformer applies a relationship between flow in flow out, and effort in effort out. Examples include an ideal electrical transformer or a lever.

Denoted $$\begin{array}{c}\stackrel{{}_1}{\underset{}{---\rightharpoonup }}TR\stackrel{{}_2}{\underset{}{---\rightharpoonup }}\\ {}^{r:1}\end{array}$$ where the r denotes the modulus of the transformer. This means $$f_{1}r=f_2$$ and $$e_{2}r=e_1$$

Gyrator

A gyrator applies a relationship between flow in effort out, and effort in flow out. An example of a gyrator is a DC motor, which converts voltage (electrical effort) into angular velocity (angular mechanical flow).

$$\begin{array}{c}\stackrel{{}_1}{\underset{}{---\rightharpoonup }}GY\stackrel{{}_2}{\underset{}{---\rightharpoonup }}\\ {}^{g:1}\end{array}$$ meaning that $$e_2=g{f}_1$$ and $$e_1=g{f}_2.$$

Multi-port elements

Junctions, unlike the other elements can have any number of ports either in or out. Junctions split power across their ports. There are two distinct junctions, the 0-junction and the 1-junction which differ only in how effort and flow are carried across. The same junction in series can be combined, but different junctions in series cannot.

0-junctions

0-junctions behave such that all effort values (and its time integral/derivative) are equal across the bonds, but the sum of the flow values in equals the sum of the flow values out, or equivalently, all flows sum to zero. In an electrical circuit, the 0-junction is a node and represents a voltage shared by all components at that node. In a mechanical circuit, the 0-junction is a joint among components, and represents a force shared by all components connected to it.

$$\text{all }e\text{'s are equal}$$ $$\sum f_{\text{in}}=\sum f_{\text{out}}$$

An example is shown below.

$$\stackrel{{}_1}{\underset{}{---\rightharpoondown }}\stackrel{\stackrel{{}_2}{\upharpoonright }}{0}\stackrel{{}_3}{\underset{}{---\rightharpoondown }}$$

Resulting equations: $$e_1=e_2=e_3$$ $$f_1=f_2+f_3$$

1-junctions

1-junctions behave opposite of 0-junctions. 1-junctions behave such that all flow values (and its time integral/derivative) are equal across the bonds, but the sum of the effort values in equals the sum the effort values out, or equivalently, all efforts sum to zero. In an electrical circuit, the 1 junction represents a series connection among components. In a mechanical circuit, the 1-junction represents a velocity shared by all components connected to it.

$$\text{all }f\text{'s are equal}$$ $$\sum e_{\text{in}}=\sum e_{\text{out}}$$

An example is shown below.

$$\stackrel{{}_1}{\underset{}{---\rightharpoondown }}\stackrel{\stackrel{{}_2}{\upharpoonright }}{1}\stackrel{{}_3}{\underset{}{---\rightharpoondown }}$$

Resulting equations: $$f_1=f_2=f_3$$ $$e_1=e_2+e_3$$

Causality

Bond graphs have a notion of causality, indicating which side of a bond determines the instantaneous effort and which determines the instantaneous flow. In formulating the dynamic equations that describe the system, causality defines, for each modeling element, which variable is dependent and which is independent. By propagating the causation graphically from one modeling element to the other, analysis of large-scale models becomes easier. Completing causal assignment in a bond graph model will allow the detection of modeling situation where an algebraic loop exists; that is the situation when a variable is defined recursively as a function of itself.

As an example of causality, consider a capacitor in series with a battery. It is not physically possible to charge a capacitor instantly, so anything connected in parallel with a capacitor will necessarily have the same voltage (effort variable) as that across the capacitor. Similarly, an inductor cannot change flux instantly and so any component in series with an inductor will necessarily have the same flow as the inductor. Because capacitors and inductors are passive devices, they cannot maintain their respective voltage and flow indefinitely—the components to which they are attached will affect their respective voltage and flow, but only indirectly by affecting their current and voltage respectively.

Note: Causality is a symmetric relationship. When one side "causes" effort, the other side "causes" flow.

In bond graph notation, a causal stroke may be added to one end of the power bond to indicate that this side is defining the flow. Consequently, the side opposite from the casual stroke controls the effort.

Sources of flow ($S_f$) define flow, so they host the causal stroke: $$S_f|---\rightharpoonup$$ Sources of effort ($S_e$) define effort, so the other end hosts the causal stroke: $$S_e---\rightharpoonup |$$

Consider a constant-torque motor driving a wheel, i.e. a source of effort ($S_e$). That would be drawn as follows: $$\begin{array}{r}\text{motor}\\ S_e\end{array}\stackrel{\tau }{\underset{\omega }{---\rightharpoonup |}}\text{wheel}$$

Symmetrically, the side with the causal stroke (in this case the wheel) defines the flow for the bond.

Causality results in compatibility constraints. Clearly only one end of a power bond can define the effort and so only one end of a bond can (the other end) have a causal stroke. In addition, the two passive components with time-dependent behavior, $I$ and $C$, can only have one sort of causation: an $I$ component determines flow; a $C$ component defines effort. So from a junction, $J$, the preferred causal orientation is as follows: $$J\stackrel{}{\underset{}{---\rightharpoonup |}}I\text{and}J\stackrel{}{\underset{}{|---\rightharpoonup }}C$$

The reason that this is the preferred method for these elements can be further analyzed if you consider the equations they would give shown by the tetrahedron of state.

$$f(t)=\frac{1}{I}\int e(t)dt\text{and}e(t)=\frac{1}{C}\int f(t)dt$$

The resulting equations involve the integral of the independent power variable. This is preferred over the result of having the causality the other way, which results in derivative. The equations can be seen below.

$$e(t)=I\dot{f}(t)\text{and}f(t)=C\dot{e}(t)$$

It is possible for a bond graph to have a causal bar on one of these elements in the non-preferred manner. In such a case a "causal conflict" is said to have occurred at that bond. The results of a causal conflict are only seen when writing the state-space equations for the graph. It is explained in more details in that section.

A resistor has no time-dependent behavior: apply a voltage and get a flow instantly, or apply a flow and get a voltage instantly, thus a resistor can be at either end of a causal bond: $$J\stackrel{}{\underset{}{---\rightharpoonup |}}R\text{and}J\stackrel{}{\underset{}{|---\rightharpoonup }}R$$

Transformers are passive, neither dissipating nor storing energy, so causality passes through them: $$\stackrel{}{\underset{}{-----|}}TF\stackrel{}{\underset{}{-----|}}\text{or}\stackrel{}{\underset{}{|-----}}TF\stackrel{}{\underset{}{|-----}}$$

A gyrator transforms flow to effort and effort to flow, so if flow is caused on one side, effort is caused on the other side and vice versa: $$\stackrel{}{\underset{}{|-----}}GY\stackrel{}{\underset{}{-----|}}\text{or}\stackrel{}{\underset{}{-----|}}GY\stackrel{}{\underset{}{|-----}}$$

Junctions

In a 0-junction, efforts are equal; in a 1-junction, flows are equal. Thus, with causal bonds, only one bond can cause the effort in a 0-junction and only one can cause the flow in a 1-junction. Thus, if the causality of one bond of a junction is known, the causality of the others is also known. That one bond is called the 'strong bond' $$\text{strong bond}\to ⊣\stackrel{\mathrm{\top }}{\underset{\mathrm{\perp }}{0}}⊣\text{and}\text{strong bond}\to ⊢\stackrel{\mathrm{\perp }}{\underset{\mathrm{\top }}{1}}⊢$$ In a nutshell, 0-junctions must have a single causal bar, 1-junctions must have all but one causal bars.

Determining causality

In order to determine the causality of a bond graph certain steps must be followed. Those steps are:

1. Draw Source Causal Bars
2. Draw Preferred causality for C and I bonds
3. Draw causal bars for 0 and 1 junctions, transformers and gyrators
4. Draw R bond causal bars
5. If a causal conflict occurs, change C or I bond to differentiation

A walk-through of the steps is shown below. $$\begin{array}{ccccccccc}{S}_f& \stackrel{}{\underset{}{---\rightharpoonup }}& 0& \stackrel{}{\underset{}{---\rightharpoonup }}& TR& \stackrel{}{\underset{}{---\rightharpoonup }}& 0& \stackrel{}{\underset{}{---\rightharpoonup }}& C_5\\ & & \downharpoonleft & & {}^{r:1}& & \downharpoonleft & & \\ & & C_2& & & & R_6& & \end{array}$$

The first step is to draw causality for the sources, over which there is only one. This results in the graph below. $$\begin{array}{ccccccccc}{S}_f& \stackrel{}{\underset{}{|---\rightharpoonup }}& 0& \stackrel{}{\underset{}{---\rightharpoonup }}& TR& \stackrel{}{\underset{}{---\rightharpoonup }}& 0& \stackrel{}{\underset{}{---\rightharpoonup }}& C_5\\ & & \downharpoonleft & & {}^{r:1}& & \downharpoonleft & & \\ & & C_2& & & & R_6& & \end{array}$$

The next step is to draw the preferred causality for the C bonds. $$\begin{array}{ccccccccc}{S}_f& \stackrel{}{\underset{}{|---\rightharpoonup }}& 0& \stackrel{}{\underset{}{---\rightharpoonup }}& TR& \stackrel{}{\underset{}{---\rightharpoonup }}& 0& \stackrel{}{\underset{}{|---\rightharpoonup }}& C_5\\ & & \overline{\downharpoonleft }& & {}^{r:1}& & \downharpoonleft & & \\ & & C_2& & & & R_6& & \end{array}$$

Next apply the causality for the 0 and 1 junctions, transformers, and gyrators. $$\begin{array}{ccccccccc}{S}_f& \stackrel{}{\underset{}{|---\rightharpoonup }}& 0& \stackrel{}{\underset{}{|---\rightharpoonup }}& TR& \stackrel{}{\underset{}{|---\rightharpoonup }}& 0& \stackrel{}{\underset{}{|---\rightharpoonup }}& C_5\\ & & \overline{\downharpoonleft }& & {}^{r:1}& & \underset{\_}{\downharpoonleft }& & \\ & & C_2& & & & R_6& & \end{array}$$

However, there is an issue with 0-junction on the left. The 0-junction has two causal bars at the junction, but the 0-junction wants one and only one at the junction. This was caused by having $C_2$ be in the preferred causality. The only way to fix this is to flip that causal bar. This results in a causal conflict, the corrected version of the graph is below, with the $\star$ representing the causal conflict.

$$\begin{array}{ccccccccc}{S}_f& \stackrel{}{\underset{}{|---\rightharpoonup }}& 0& \stackrel{}{\underset{}{|---\rightharpoonup }}& TR& \stackrel{}{\underset{}{|---\rightharpoonup }}& 0& \stackrel{}{\underset{}{|---\rightharpoonup }}& C_5\\ & & \underset{\_}{\downharpoonleft }\star & & {}^{r:1}& & \underset{\_}{\downharpoonleft }& & \\ & & C_2& & & & R_6& & \end{array}$$

Converting from other systems

One of the main advantages of using bond graphs is that once you have a bond graph it doesn't matter the original energy domain. Below are some of the steps to apply when converting from the energy domain to a bond graph.

Electromagnetic

The steps for solving an Electromagnetic problem as a bond graph are as follows:

1. Place an 0-junction at each node
2. Insert Sources, R, I, C, TR, and GY bonds with 1 junctions
3. Ground (both sides if a transformer or gyrator is present)
4. Assign power flow direction
5. Simplify

These steps are shown more clearly in the examples below.

Linear mechanical

The steps for solving a Linear Mechanical problem as a bond graph are as follows:

1. Place 1-junctions for each distinct velocity (usually at a mass)
2. Insert R and C bonds at their own 0-junctions between the 1 junctions where they act
3. Insert Sources and I bonds on the 1 junctions where they act
4. Assign power flow direction
5. Simplify

These steps are shown more clearly in the examples below.

Simplifying

The simplifying step is the same regardless if the system was electromagnetic or linear mechanical. The steps are:

1. Remove Bond of zero power (due to ground or zero velocity)
2. Remove 0 and 1 junctions with less than three bonds
3. Simplify parallel power
4. Combine 0 junctions in series
5. Combine 1 junctions in series

These steps are shown more clearly in the examples below.

Parallel power

Parallel power is when power runs in parallel in a bond graph. An example of parallel power is shown below.

Parallel power can be simplified, by recalling the relationship between effort and flow for 0 and 1-junctions. To solve parallel power you will first want to write down all of the equations for the junctions. For the example provided, the equations can be seen below. (Please make note of the number bond the effort/flow variable represents). $$\begin{array}{ccc}{f}_1=f_2=f_3& & e_2=e_4=e_7\\ e_1=e_2+e_3& & f_2=f_4+f_7\\ & & \\ e_3=e_5=e_6& & f_7=f_6=f_8\\ f_3=f_5+f_6& & e_7+e_6=e_8\end{array}$$

By manipulating these equations you can arrange them such that you can find an equivalent set of 0 and 1-junctions to describe the parallel power.

For example, because $e_3=e_6$ and $e_2=e_7$ you can replace the variables in the equation $e_1=e_2+e_3$ resulting in $e_1=e_6+e_7$ and since $e_6+e_7=e_8$, we now know that $e_1=e_8$. This relationship of two effort variables equaling can be explained by an 0-junction. Manipulating other equations you can find that $f_4=f_5$ which describes the relationship of a 1-junction. Once you have determined the relationships that you need you can redraw the parallel power section with the new junctions. The result for the example show is seen below.

Examples

Simple electrical system

A simple electrical circuit consisting of a voltage source, resistor, and capacitor in series.

The first step is to draw 0-junctions at all of the nodes: $$\begin{array}{ccccc}& 0& & 0& \\ & & & & \\ & & & & \\ & 0& & 0& \end{array}$$

The next step is to add all of the elements acting at their own 1-junction: $$\begin{array}{ccccccccc}& & & & R& & & & \\ & & & & |& & & & \\ & & 0& -& 1& -& 0& & \\ & & |& & & & |& & \\ S_e& -& 1& & & & 1& -& C\\ & & |& & & & |& & \\ & & \underset{\_}{0}& -& -& -& 0& & \end{array}$$

The next step is to pick a ground. The ground is simply an 0-junction that is going to be assumed to have no voltage. For this case, the ground will be chosen to be the lower left 0-junction, that is underlined above. The next step is to draw all of the arrows for the bond graph. The arrows on junctions should point towards ground (following a similar path to current). For resistance, inertance, and compliance elements, the arrows always point towards the elements. The result of drawing the arrows can be seen below, with the 0-junction marked with a star as the ground.

Now that we have the Bond graph, we can start the process of simplifying it. The first step is to remove all the ground nodes. Both of the bottom 0-junctions can be removed, because they are both grounded. The result is shown below.

Next, the junctions with less than three bonds can be removed. This is because flow and effort pass through these junctions without being modified, so they can be removed to allow us to draw less. The result can be seen below.

The final step is to apply causality to the bond graph. Applying causality was explained above. The final bond graph is shown below.

Advanced electrical system

A more advanced electrical system with a current source, resistors, capacitors, and a transformer

Following the steps with this circuit will result in the bond graph below, before it is simplified. The nodes marked with the star denote the ground.

Simplifying the bond graph will result in the image below.

Lastly, applying causality will result in the bond graph below. The bond with star denotes a causal conflict.

Simple linear mechanical

A simple linear mechanical system, consisting of a mass on a spring that is attached to a wall. The mass has some force being applied to it. An image of the system is shown below.

For a mechanical system, the first step is to place a 1-junction at each distinct velocity, in this case there are two distinct velocities, the mass and the wall. It is usually helpful to label the 1-junctions for reference. The result is below. $$\begin{array}{ccc}& & \\ & & \\ 1_{\text{mass}}& & \\ & & \\ & & \\ & & \\ 1_{\text{wall}}& & \end{array}$$

The next step is to draw the R and C bonds at their own 0-junctions between the 1-junctions where they act. For this example there is only one of these bonds, the C bond for the spring. It acts between the 1-junction representing the mass and the 1-junction representing the wall. The result is below. $$\begin{array}{ccc}& & \\ & & \\ 1_{\text{mass}}& & \\ |& & \\ 0& -& C:\frac{1}{k}\\ |& & \\ 1_{\text{wall}}& & \end{array}$$

Next you want to add the sources and I bonds on the 1-junction where they act. There is one source, the source of effort (force) and one I bond, the mass of the mass both of which act on the 1-junction of the mass. The result is shown below. $$\begin{array}{ccc}{S}_e:F(t)& & \\ |& & \\ 1_{\text{mass}}& -& I:m\\ |& & \\ 0& -& C:\frac{1}{k}\\ |& & \\ 1_{\text{wall}}& & \end{array}$$

Next power flow is to be assigned. Like the electrical examples, power should flow towards ground, in this case the 1-junction of the wall. Exceptions to this are R, C, or I bond, which always point towards the element. The resulting bond graph is below.

Now that the bond graph has been generated, it can be simplified. Because the wall is grounded (has zero velocity), you can remove that junction. As such the 0-junction the C bond is on, can also be removed because it will then have less than three bonds. The simplified bond graph can be seen below.

The last step is to apply causality, the final bond graph can be seen below.

Advanced linear mechanical

A more advanced linear mechanical system can be seen below.

Just like the above example, the first step is to make 1-junctions at each of the distant velocities. In this example there are three distant velocity, Mass 1, Mass 2, and the wall. Then you connect all of the bonds and assign power flow. The bond can be seen below.

Next you start the process of simplifying the bond graph, by removing the 1-junction of the wall, and removing junctions with less than three bonds. The bond graph can be seen below.

There is parallel power in the bond graph. Solving parallel power was explained above. The result of solving it can be seen below.

Lastly, apply causality, the final bond graph can be seen below.

State equations

Once a bond graph is complete, it can be utilized to generate the state-space representation equations of the system. State-space representation is especially powerful as it allows complex multi-order differential system to be solved as a system of first-order equations instead. The general form of the state equation is $$\dot{\mathbf{𝐱}}(t)=\mathbf{𝐀}\mathbf{𝐱}(t)+\mathbf{𝐁}\mathbf{𝐮}(t)$$ where $\mathbf{𝐱}(t)$ is a column matrix of the state variables, or the unknowns of the system. $\dot{\mathbf{𝐱}}(t)$ is the time derivative of the state variables. $\mathbf{𝐮}(t)$ is a column matrix of the inputs of the system. And $\mathbf{𝐀}$ and $\mathbf{𝐁}$ are matrices of constants based on the system. The state variables of a system are $q(t)$ and $p(t)$ values for each C and I bond without a causal conflict. Each I bond gets a $p(t)$ while each C bond gets a $q(t)$.

For example, if you have the following bond graph

you would have the following $\dot{\mathbf{𝐱}}(t)$, $\mathbf{𝐱}(t)$, and $\mathbf{𝐮}(t)$ matrices:

$$\dot{\mathbf{𝐱}}(t)=\left[\begin{array}{c}{\dot{p}}_3(t)\\ {\dot{q}}_6(t)\end{array}\right]\text{and}\mathbf{𝐱}(t)=\left[\begin{array}{c}{p}_3(t)\\ q_6(t)\end{array}\right]\text{and}\mathbf{𝐮}(t)=\left[\begin{array}{c}{e}_1(t)\end{array}\right]$$

The matrices of $\mathbf{𝐀}$ and $\mathbf{𝐁}$ are solved by determining the relationship of the state variables and their respective elements, as was described in the tetrahedron of state. The first step to solve the state equations is to list all of the governing equations for the bond graph. The table below shows the relationship between bonds and their governing equations.

Element type	Bond Name	Bond with causality	Governing equation(s)
Single-port elements	Source/ Sink, S	$S_e\stackrel{}{\underset{}{---\rightharpoonup |}}$	$\text{input}=e(t)$
		$S_f\stackrel{}{\underset{}{|---\rightharpoonup }}$	$\text{input}=f(t)$
	Resistance, R: Dissipated Energy	$\stackrel{}{\underset{}{---\rightharpoonup |}}R$	$f(t)=\frac{1}{R}e(t)$
		$\stackrel{}{\underset{}{|---\rightharpoonup }}R$	$e(t)=Rf(t)$
	Inertance, I: Kinetic Energy	♦$\stackrel{}{\underset{}{---\rightharpoonup |}}I$	$f(t)=\frac{1}{I}\int e(t)dt$
		$\stackrel{}{\underset{}{|---\rightharpoonup }}I$	$e(t)=I\dot{f}(t)$
	Compliance, C: Potential Energy	$\stackrel{}{\underset{}{---\rightharpoonup |}}C$	$f(t)=C\dot{e}(t)$
		♦$\stackrel{}{\underset{}{|---\rightharpoonup }}C$	$e(t)=\frac{1}{C}\int f(t)dt$
Double-port elements	Transformer, TR	$\begin{array}{c}\stackrel{}{\underset{}{---\rightharpoonup }|}TR\stackrel{}{\underset{}{---\rightharpoonup }}|\\ {}^{r:1}\end{array}$	$f_1=\frac{1}{r}{f}_2$ $e_2=\frac{1}{r}{e}_1$
		$\begin{array}{c}|\stackrel{}{\underset{}{---\rightharpoonup }}TR|\stackrel{}{\underset{}{---\rightharpoonup }}\\ {}^{r:1}\end{array}$	$e_1=r{e}_2$ $f_2=r{f}_1$
	Gyrator, GY	$\begin{array}{c}|\stackrel{}{\underset{}{---\rightharpoonup }}GY\stackrel{}{\underset{}{---\rightharpoonup }}|\\ {}^{g:1}\end{array}$	$e_1=g{f}_2$ $e_2=g{f}_1$
		$\begin{array}{c}\stackrel{}{\underset{}{---\rightharpoonup }}|GY\stackrel{}{\underset{}{|---\rightharpoonup }}\\ {}^{g:1}\end{array}$	$f_1=\frac{1}{g}{e}_2$ $f_2=\frac{1}{g}{e}_1$
Multi-port elements	0 junction	One and only one causal bar at the junction	$\text{all }e(t)=$
			$\sum f(t{)}_{\text{in}}=\sum f(t{)}_{\text{out}}$
	1 junction	one and only one causal bar away from the junction	$\sum e(t{)}_{\text{in}}=\sum e(t{)}_{\text{out}}$
			$\text{all }f(t)=$

"♦" denotes preferred causality.

For the example provided,

the governing equations are the following.

1. $e_1=\text{input}$
2. $e_3=e_1-e_2-e_4$
3. $f_1=f_2=f_4=f_3$
4. $e_2=R_2{f}_2$
5. $f_3=\frac{1}{I_3}\int e_{3}dt=\frac{1}{I_3}{p}_3$
6. $f_5=f_4\cdot r$
7. $e_4=e_5\cdot r$
8. $e_5=e_7=e_6$
9. $f_6=f_7-f_5$
10. $e_6=\frac{1}{C_6}\int f_{6}dt=\frac{1}{C_6}{q}_6$
11. $f_7=\frac{1}{R_7}{e}_7$

These equations can be manipulated to yield the state equations. For this example, you are trying to find equations that relate ${\dot{p}}_3(t)$ and ${\dot{q}}_6(t)$ in terms of $p_3(t)$, $q_6(t)$, and $e_1(t)$.

To start you should recall from the tetrahedron of state that ${\dot{p}}_3(t)=e_3(t)$starting with equation 2, you can rearrange it so that $e_3=e_1-e_2-e_4$. $e_2$ can be substituted for equation 4, while in equation 4, $f_2$ can be replaced by $f_3$due to equation 3, which can then be replaced by equation 5. $e_4$ can likewise be replaced using equation 7, in which $e_5$can be replaced with $e_6$ which can then be replaced with equation 10. Following these substituted yields the first state equation which is shown below.

$${\dot{p}}_3(t)=e_3(t)=e_1(t)-\frac{R_2}{I_3}{p}_3(t)-\frac{r}{C_6}{q}_6(t)$$

The second state equation can likewise be solved, by recalling that ${\dot{q}}_6(t)=f_6(t)$. The second state equation is shown below.

$${\dot{q}}_6(t)=f_6(t)=\frac{r}{I_3}{p}_3(t)-\frac{1}{R_7\cdot C_6}{q}_6(t)$$

Both equations can further be rearranged into matrix form. The result of which is below.

$$\left[\begin{array}{c}{\dot{p}}_3(t)\\ {\dot{q}}_6(t)\end{array}\right]=\left[\begin{array}{cc}-\frac{R_2}{I_3}& -\frac{r}{C_6}\\ \frac{r}{I_3}& -\frac{1}{R_7\cdot C_6}\end{array}\right]\left[\begin{array}{c}{p}_3(t)\\ q_6(t)\end{array}\right]+\left[\begin{array}{c}1\\ 0\end{array}\right]\left[\begin{array}{c}{e}_1(t)\end{array}\right]$$

At this point the equations can be treated as any other state-space representation problem.

International conferences on bond graph modeling (ECMS and ICBGM)

A bibliography on bond graph modeling may be extracted from the following conferences :

• ECMS-2013 27th European Conference on Modelling and Simulation, May 27–30, 2013, Ålesund, Norway
• ECMS-2008 22nd European Conference on Modelling and Simulation, June 3–6, 2008 Nicosia, Cyprus
• ICBGM-2007: 8th International Conference on Bond Graph Modeling And Simulation, January 15–17, 2007, San Diego, California, U.S.A.
• ECMS-2006 20TH European Conference on Modelling and Simulation, May 28–31, 2006, Bonn, Germany
• IMAACA-2005 International Mediterranean Modeling Multiconference
• ICBGM-2005 International Conference on Bond Graph Modeling and Simulation, January 23–27, 2005, New Orleans, Louisiana, U.S.A. – Papers
• ICBGM-2003 International Conference on Bond Graph Modeling and Simulation (ICBGM'2003) January 19–23, 2003, Orlando, Florida, USA – Papers
• 14TH European Simulation symposium October 23–26, 2002 Dresden, Germany
• ESS'2001 13th European Simulation symposium, Marseilles, France October 18–20, 2001
• ICBGM-2001 International Conference on Bond Graph Modeling and Simulation (ICBGM 2001), Phoenix, Arizona U.S.A.
• European Simulation Multi-conference 23-26 May, 2000, Gent, Belgium
• 11th European Simulation symposium, October 26–28, 1999 Castle, Friedrich-Alexander University, Erlangen-Nuremberg, Germany
• ICBGM-1999 International Conference on Bond Graph Modeling and Simulation January 17–20, 1999 San Francisco, California
• ESS-97 9TH European Simulation Symposium and Exhibition Simulation in Industry, Passau, Germany, October 19–22, 1997
• ICBGM-1997 3rd International Conference on Bond Graph Modeling And Simulation, January 12–15, 1997, Sheraton-Crescent Hotel, Phoenix, Arizona
• 11th European Simulation Multiconference Istanbul, Turkey, June 1–4, 1997
• ESM-1996 10th annual European Simulation Multiconference Budapest, Hungary, June 2–6, 1996
• ICBGM-1995 Int. Conf. on Bond Graph Modeling and Simulation (ICBGM’95), January 15–18, 1995, Las Vegas, Nevada.

See also

• 20-sim simulation software based on the bond graph theory
• AMESim simulation software based on the bond graph theory
• Hybrid bond graph
• Coenergy

Further reading

• Kypuros, Javier (2013). System dynamics and control with bond graph modeling. Boca Raton: Taylor&Francis. doi:10.1201/b14676. ISBN 978-1-4665-6075-8. 
• Paynter, Henry M. (1960). Analysis and design of engineering systems. M.I.T. Press. ISBN 0-262-16004-8. 
• Karnopp, Dean C.; Margolis, Donald L.; Rosenberg, Ronald C. (1990). System dynamics: a unified approach. New York: John Wiley & Sons. ISBN 0-471-62171-4. https://archive.org/details/systemdynamicsun0000karn. 
• Thoma, Jean Ulrich (1975). Bond graphs: introduction and applications. Oxford: Pergamon Press. ISBN 0-08-018882-6. 
• Gawthrop, Peter J.; Smith, Lorcan P. S. (1996). Metamodelling: bond graphs and dynamic systems. London: Prentice Hall. ISBN 0-13-489824-9. 
• Brown, Forbes T. (2007). Engineering system dynamics – a unified graph-centered approach. Boca Raton: Taylor & Francis. ISBN 0-8493-9648-4. 
• Mukherjee, Amalendu; Karmakar, Ranjit (2000). Modelling and simulation of engineering systems through bondgraphs. Boca Raton: CRC Press. ISBN 978-0-8493-0982-3. 
• Gawthrop, P.J.; Ballance, D.J. (1999). "Chapter 2: Symbolic computation for manipulation of hierarchical bond graphs". in Munro, N.. Symbolic Methods in Control System Analysis and Design. London: Institution of Electrical Engineers. pp. 23-52. ISBN 0-85296-943-0. https://archive.org/details/symbolicmethodsc00munr. 
• Borutzky, Wolfgang (2010). Bond Graph Methodology. London: Springer. doi:10.1007/978-1-84882-882-7. ISBN 978-1-84882-881-0. 
• http://www.site.uottawa.ca/~rhabash/ESSModelFluid.pdf Explains modeling the bond graph in the fluid domain
• http://www.dartmouth.edu/~sullivan/22files/Fluid_sys_anal_w_chart.pdf Explains modeling the bond graph in the fluid domain

External links

• Simscape Official MATLAB/Simulink add-on library for graphical bond graph programming
• BG V.2.1 Freeware MATLAB/Simulink add-on library for graphical bond graph programming