{"id":76,"date":"2021-02-04T14:28:01","date_gmt":"2021-02-04T19:28:01","guid":{"rendered":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/?post_type=chapter&p=76"},"modified":"2024-03-28T17:47:29","modified_gmt":"2024-03-28T21:47:29","slug":"miscellaneous-topics","status":"publish","type":"chapter","link":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/chapter\/miscellaneous-topics\/","title":{"raw":"Miscellaneous Topics","rendered":"Miscellaneous Topics"},"content":{"raw":"
\r\n

11.1\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Overview<\/a><\/h1>\r\nThis chapter discusses several topics pertinent to modelling and simulation of systems dynamics with focus on BG method. The reader can consult with these sections as independent topics or as supplementary to previously discussed topics.\r\n

11.2\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Energy and Power Conjugate Variables<\/h1>\r\nThe product of conjugate variables related to a physical quantity should, by definition, give the corresponding quantity. In chapter 3, sections 3.2<\/a> and 3.5<\/a>, we discussed the power variables (effort $e$ and flow $f$) and their relations to the state variables (momentum $p$ and displacement $q$). The product ($e.f$) is power\u2014so-called power variables\u2014and integrating them, with respect to time, gives $p$ and $q$, respectively. Hence the $p$ and $q$ are called energy variables. However, further attention to the dimension of the product of the variables, i.e., ($p.q$), indicates that the dimension of this quantity is not equivalent to that of energy. For example, in the SI system of units, we get $[p.q] = [N.s.m] = [J.s]$. Therefore, we can call $e$ and $f$ conjugate power variables because the product ($e.f$) is power. But this definition does not apply to the product ($p.q$). Instead, the product ($\\dot p . q$) or product ($p. \\dot q$) can be defined as true conjugate energy variables.\r\n\r\nOne can ask the question, what is the relationship among these variables? To investigate further and provide a possible answer, we consider the system energy function $\\mathcal{E} = \\mathcal{E} (p,q)$. This functional form is legitimate, since we know that the energy stored, in storage elements $I$ and $C$, is uniquely defined by state variables (see section 3.5<\/a>). The total differential change of energy is then $d \\mathcal{E} = \\dfrac{\\partial \\mathcal{E}}{\\partial p} dp + \\dfrac{\\partial \\mathcal{E}}{\\partial q} dq$. In principle, $\\mathcal{E} = \\mathcal{E}(p,q,t)$. But the partial derivative $\\dfrac{\\partial \\mathcal{E}}{\\partial t} |_{p,q}$ can be ignored at equilibrium state of the system, or when $p$ and $q$ are kept constant and energy is not changing with time[footnote]$d \\mathcal{E} (p,q,t) = \\dfrac{\\partial \\mathcal{E}}{\\partial p} dp + \\dfrac{\\partial \\mathcal{E}}{\\partial q} dq + \\dfrac{\\partial \\mathcal{E}}{\\partial t} dt$[\/footnote]. Now, the total differential change of power $d \\mathcal{P}$ can be written as the time derivative of energy, or $d \\mathcal{P} = \\dfrac{d \\mathcal{E}}{dt} = \\dfrac{\\partial \\mathcal{E}}{\\partial p} \\dfrac{dp}{dt} + \\dfrac{\\partial \\mathcal{E}}{\\partial q} \\dfrac{dq}{dt}$. But $\\dfrac{dp}{dt} = e$ and $\\dfrac{dq}{dt} = f$, the effort and flow in BG method, respectively. Therefore, we can write\r\n\r\n\\begin{equation}\r\n\r\nd \\mathcal{P} = e \\dfrac{\\partial \\mathcal{E}}{\\partial p} + f \\dfrac{\\partial \\mathcal{E}}{\\partial q}\r\n\r\n\\tag{11.1}\r\n\r\n\\end{equation}\r\n\r\nEquation (11.1) clearly defines the relationship among power variables (effort $e$ and flow $f$) and state variables (momentum $p$ and displacement $q$) with respect to the power ($\\mathcal{P}$) and energy ($\\mathcal{E}$) of the system. Note that for a system, we can write this equation for each component and sum it up to get the relation for the system.\r\n\r\nFrom Equation (11.1), we can conclude that by dimensional analysis, $\\dfrac{\\partial \\mathcal{E}}{\\partial p}$ has the dimension of flow and $\\dfrac{\\partial \\mathcal{E}}{\\partial q}$ the dimension of effort, according to BG terminology. We define these quantities as system flow and effort, or (note that the partial derivative definition is explicitly shown here)\r\n\r\n\\begin{equation}\r\n\r\n\\begin{dcases}\r\n\r\n\\dfrac{\\partial \\mathcal{E}}{\\partial p} |_{q=constant} \u00a0, \\: system \\: flow \\\\\r\n\r\n\\dfrac{\\partial \\mathcal{E}}{\\partial q} |_{p=constant} \u00a0, \\: system \\: effort\r\n\r\n\\end{dcases}\r\n\r\n\\tag{11.2}\r\n\r\n\\end{equation}\r\n\r\nFor example, considering a mass, represented by an $I$-element, its kinetic energy can be written as $\\dfrac{p^2}{2m}$, but $\\dfrac{\\partial}{\\partial p} \\left( \\dfrac{p^2}{2m} \\right) = \\dfrac{p}{m} = \\dfrac{mv}{m} = v$, the velocity of the mass or the flow associated with the $I$-element. Similarly, considering a mechanical spring with stiffness $k$, represented by a $C$-element, its potential energy can be written as $k \\frac{q^2}{2}$, but $\\dfrac{\\partial}{\\partial q} \\left( k \\dfrac{q^2}{2} \\right) = kq = kx = F$, the force acting on the spring or the effort associated with the $C$-element. Equation (11.2) applies analogously to mechanical, electrical, hydraulic, or thermal systems. Note that the energy function $\\mathcal{E}$ relates to the Hamiltonian of the system, [13]<\/strong><\/a>.\r\n

11.3\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Including Gravity in BG Models<\/h1>\r\nForces due to gravitational acceleration $g$ can be implemented in a BG model using effort source element $S_e$. But it feeds power to the system when the assigned positive direction of displacement\/velocity and the gravitational acceleration vector are in the same direction. Otherwise, the system feeds power into the gravity field, or loses power. In manually drawn BG, the direction of the power bond is drawn from $S_e$ toward the system when the gravity and displacement directions are matched; otherwise, the power bond should be toward the $S_e$ element. In 20-sim, however, it is not possible to draw a power bond pointing toward an $S_e$ element because it is defined as a source element. Hence, a negative value of $g$ is assigned to the associated parameter in the corresponding equation model, i.e., \u22129.81 $m\/s^2$ (see Figure 11-1<\/a>). Note that $g$ is interpreted as gravitational force per unit mass, here.\u00a0<\/a>\r\n\r\n[caption id=\"attachment_1235\" align=\"aligncenter\" width=\"602\"]\"\"<\/a> Figure 11-1 Implementing gravity force for a BG model in 20-sim.[\/caption]\r\n

<\/a>11.4 \u00a0 \u00a0 \u00a0 \u00a0 Extracting System Equations from BG Models<\/h1>\r\nAlthough graphics of BG models can give us some insights into the dynamics of systems before actually solving the systems equations, having the mathematical model consisting of a system of equations is the ultimate goal resulting from a BG model, or for that matter, any model. In the previous sections, we discussed that for a BG model with integral causalities, we can derive a system of ODEs in terms of state variables, i.e., $p$ and $q$ (see section 3.5<\/a>). For such a system, the number of equations, $n$ is equal to the sum of number of $p$ associated with $I$- elements and number of $q$ associated with the $C$- elements. In other words, the system equations are coupled first-order ODEs explicitly given in terms of state variables and can be solved simultaneously using numerical schemes. The solution of these coupled equations is, obviously, easier to obtain for linear systems rather than it is for non-linear systems.\r\n\r\nIn this section we discuss the procedure and guidelines of how to extract system equations from a BG model. For now, we assume that all storage elements (i.e., $I$- and $C$- elements) in the system are assigned with integral causality strokes. In the next section, we discuss cases when derivative causality and\/or algebraic-loop situations may exist in a BG model, along with the definitions and resulting consequences for such systems equations.\r\n\r\nReferring to the energy management of a system, we have already discussed that the energy input to a system is partially stored in storages and the rest is dissipated through dissipaters. Therefore, the storage elements are key elements for studying the dynamics of any system. This principle also applies to the BG models, and we use it to extract system equations.\r\n\r\nThe system equations can be given in a general form as $\\dot x_i = \\mathcal{F}_i (x_i, u_j, t)$ where $\\dot x_i$, the vector representing time derivative of state variables, $u_j$ , the vector representing the inputs, with indices $i = 1, \\cdots , n,$ and $j = 1, \\cdots, r$. One can expand this functional to get\r\n\r\n\\begin{equation}\r\n\r\n\\begin{dcases}\r\n\r\n\\dot x_1 = \\mathcal{F}_1 (x_1, x_2, \\dots, x_n, u_1, u_2, \\dots, u_r, t) \\\\\r\n\r\n\\dot x_2 = \\mathcal{F}_2 (x_1, x_2, \\dots, x_n, u_1, u_2, \\dots, u_r, t) \\\\\r\n\r\n\\quad \\quad \\quad \\quad \\quad \\quad \\quad \\quad \\vdots \\\\\r\n\r\n\\dot x_n = \\mathcal{F}_n (x_1, x_2, \\dots, x_n, u_1, u_2, \\dots, u_r, t) \\\\\r\n\r\n\\end{dcases}\r\n\r\n\\tag{11.3}\r\n\r\n\\end{equation}\r\n\r\nAfter building the BG model for a system, we ask the following two key questions, [21]<\/strong><\/a>:\r\n\r\nQ1: What does each component\/element send to the system?<\/em>\r\n\r\nQ2: What does the system send back to the storage components\/elements?<\/em>\r\n\r\nEach question guides us to list the corresponding relations involving the state variables and performing some manipulations to write the system equations in terms of the state variables, using the general guideline of $p$ on $I$ and $q$ on $C$. The following steps can be used to extract the system equations:\r\n
    \r\n \t
  1. Build the BG model and assign causalities, with preferred integral causalities.<\/li>\r\n \t
  2. Simplify the BG model, as far as possible.<\/li>\r\n \t
  3. Assign labels as numbers to each power bond in the BG model. The order of numbering is arbitrary.<\/li>\r\n \t
  4. Answer Q1, in terms of all efforts or flows input into the system.<\/li>\r\n \t
  5. Answer Q2, in terms of all $p$\u2019s on $I$\u2019s and all $q$\u2019s on $C$\u2019s, or momenta (or their derivatives) associated with the inertia elements and displacements (or their derivatives) associated with the capacity elements.<\/li>\r\n \t
  6. Apply constraints implemented by 1- and 0- junction elements and perform algebraic manipulations to write the system equations only in terms of state variables, as independent variables with inputs and components\u2019 properties as parameters.<\/li>\r\n<\/ol>\r\nNow, we present some worked-out examples to demonstrate the application of the above mentioned procedure.\r\n

    11.4.1\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Example: System Equations for a Mechanical System<\/h2>\r\nFor this example, we consider a mechanical system as shown in Figure 11-2<\/a>. This systems has three DOF associated with three masses. The number of state variables is seven, associated with the storage elements, three masses, and four springs. The positive displacements are considered as shown, and tension forces are positive (+T).<\/a>\r\n\r\n[caption id=\"attachment_844\" align=\"aligncenter\" width=\"463\"]\"\" Figure 11-2 A mechanical system sketch[\/caption]\r\n\r\nSolution:\r\n\r\nWe first build the BG model, either manually or using the tools available through 20-sim, with all power bonds labelled, as shown in Figure 11-3<\/a>. As mentioned, the labels are arbitrary; they merely help bookkeeping of the variables involved without affecting final solution of the system equations. Therefore, the state variables are: $p_2, p_6, p_{13}$ associated with the momenta of the masses and $q_4, q_{10}, q_{15}, q_{19}$ displacements of the mechanical springs. The corresponding power bonds to these elements are colour coded (see Figure 11-3<\/a>). The system equations are ODEs consisting of these state variables. To extract the system equations, we ask two key questions (see section 11.4<\/a>) and list their corresponding answers, using the power bond labels and causality assignments in terms of BG symbols $e$ and $f$, as follows:<\/a>\r\n\r\n \r\n\r\n[caption id=\"attachment_848\" align=\"aligncenter\" width=\"489\"]\"\" Figure 11-3 BG model for the mechanical system with labelled power bonds[\/caption]\r\n\r\nQ1: What does each component\/element send to the system?<\/em>\r\n\r\nThe inputs to the system are from elements at the boundary of the system associated with bonds 1, 2, 4, 6, 10, 11, 13, 15, 17, and 19, as listed in Table 11-1.\r\n
    \r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
    Table 11-1 Inputs to the mechanical system by all elements\u2014Q1 question<\/caption>\r\n
    Power bond label #<\/strong><\/td>\r\nSystem input\/relation<\/strong><\/td>\r\n<\/tr>\r\n<\/thead>\r\n
    1<\/td>\r\n$e_1 = F$<\/td>\r\n<\/tr>\r\n
    2<\/td>\r\n$f_2 = \\dot q_2 = p_2\/m_3$<\/td>\r\n<\/tr>\r\n
    4<\/td>\r\n$e_4 = k_4q_4$<\/td>\r\n<\/tr>\r\n
    6<\/td>\r\n$f_6 = \\dpt q_6 = p_6\/m_2$<\/td>\r\n<\/tr>\r\n
    10<\/td>\r\n$e_{10} = k_3q_{10}$<\/td>\r\n<\/tr>\r\n
    11<\/td>\r\n$e_{11} = b_1f_{11}$<\/td>\r\n<\/tr>\r\n
    13<\/td>\r\n$f_{13} = \\dot q_{13} = p_{13}\/m_1$<\/td>\r\n<\/tr>\r\n
    15<\/td>\r\n$e_{15} = k_2q_{15}$<\/td>\r\n<\/tr>\r\n
    17<\/td>\r\n$f_{17} = V_{wall}$<\/td>\r\n<\/tr>\r\n
    19<\/td>\r\n$e_{19} = k_1q_{19}$<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<\/div>\r\nNext, we list the relations to answer the second question:\r\n\r\nQ2: What does the system send back to the storage components\/elements?<\/em>\r\n\r\nHere, we are only interested in storage elements, as shown by colour-coded bonds in Figure 11-3<\/a>. For example, considering power bond number 2, the system sends an effort signal (or rate of momentum) to the $I$-element representing the mass $m_3$. Therefore, we can write $e_2 = \\dot p_2 = m_3 \\dot f_2$. Similarly, considering power bond number 15, the system sends a flow signal (or rate of displacement) to the $C$-element representing the spring $k_2$. Therefore, we can write $f_{15} = \\dot q_{15}$. Table 11-2 lists all relations corresponding to storage elements.\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
    Table 11-2 Outputs from the mechanical system to all storage elements\u2014Q2 question<\/caption>\r\n
    Power bond label #(element)<\/strong><\/td>\r\nStorage element\/relation<\/strong><\/td>\r\nRelations, using constraints<\/strong><\/td>\r\n<\/tr>\r\n<\/thead>\r\n
    2 ($I:m_3$)<\/td>\r\n$e_2 = \\dot p_2$<\/td>\r\n$e_2 = \\pmb{\\dot p_2} = e_1 - e_3 = e_1 - e_4 = \\pmb{F - k_4q_4}$<\/td>\r\n<\/tr>\r\n
    4 ($C:k_4$)<\/td>\r\n$f_4 = \\dot q_4$<\/td>\r\n$f_4 = \\pmb{\\dot q_4} = f_3 - f_5 = f_2 - f_6 = \\pmb{\\dfrac{p_2}{m_3} - \\dfrac{p_6}{m_2}}$<\/td>\r\n<\/tr>\r\n
    6 ($I:m_2$)<\/td>\r\n$e_6 = \\dot p_6$<\/td>\r\n$e_6 = \\pmb{\\dot p_6} = e_5 - e_7 - e_8 = \\pmb{k_4q_4 - k_3 q_{10}} - b_1 f_{11} - \\pmb{k_1 q_{19}}$<\/td>\r\n<\/tr>\r\n
    10 ($C:k_3$)<\/td>\r\n$f_{10} = \\dot q_{10}$<\/td>\r\n$f_{10} = \\pmb{\\dot q_{10}} = f_9 = f_7 - f_{12} = f_6 - f_{13} = \\pmb{\\dfrac{p_6}{m_2} - \\dfrac{p_{13}}{m_1}}$<\/td>\r\n<\/tr>\r\n
    13 ($I:m_1$)<\/td>\r\n$e_{13} = \\dot p_{13}$<\/td>\r\n$e_{13} = \\pmb{\\dot p_{13}} = e_{12} - e_{14} = e_9 - e_{15} = (e_{10} + e_{11}) - e_{15} = \\pmb{k_3q_{10}} + \\underbrace{b_1f_{11}}_{=b_1f_{10}} - \\pmb{k_2q_{15} = k_3q_{10}} + b_1 \\pmb{ \\left( \\dfrac{p_6}{m_2} - \\dfrac{p_{13}}{m_1} \\right) - k_2q_{15}}$<\/td>\r\n<\/tr>\r\n
    15 ($C:k_2$)<\/td>\r\n$f_{15} = \\dot q_{15}$<\/td>\r\n$f_{15} = \\pmb{\\dot q_{15}} = f_{14} - f_{16} = f_{13} - f_{17} = \\pmb{\\dfrac{p_{13}}{m_1} - V_{wall}}$<\/td>\r\n<\/tr>\r\n
    19 ($C:k_1$)<\/td>\r\n$f_{19} = \\dot q_{19}$<\/td>\r\n$f_{19} = \\pmb{\\dot q_{19}} = f_8 - f_{18} = f_6 - f_{17} = \\pmb{\\dfrac{p_6}{m_2} - V_{wall}}$<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nIn the relations listed in the third column of Table 11-2, we used the constraints resulted from 1- and 0-junctions and elements\u2019 constitutive equations, and included the relations from Table 11-1. The final relations are given in bold. Finally, after some manipulations, we can write the system equations in terms of seven state variables in matrix form, as given in Equation (11.4).\r\n\r\n\\begin{equation}\r\n\r\n\\begin{matrix}\r\n\r\n\\begin{Bmatrix} \\dot p_2\\\\ \\dot q_4 \\\\ \\dot p_6 \\\\ \\dot q_{10} \\\\ \\dot p_{13} \\\\ \\dot q_{15} \\\\ \\dot q_{19} \\end{Bmatrix}\r\n\r\n= \\begin{bmatrix} 0 & -k_4 & 0 & 0 & 0 & 0 & 0 \\\\ 1\/m_3 & 0 & -1\/m_2 & 0 & 0 & 0 & 0 \\\\ 0 & k_4 & -b_1\/m_2 & -k_3 & b_1\/m_1 & 0 & -k_1 \\\\ 0 & 0 & 1\/m_2 & 0 & -1\/m_1 & 0 & 0 \\\\ 0 & 0 & b_1\/m_2 & k_3 & -b_1\/m_1 & -k_2 & 0 \\\\ 0 & 0 & 0 & 0 & 1\/m_1 & 0 & 0 \\\\ 0 & 0 & 1\/m_2 & 0 & 0 & 0 & 0 \\end{bmatrix} \\\\\r\n\r\n\\\\\r\n\r\n\\begin{Bmatrix} p_2\\\\ q_4 \\\\ p_6 \\\\ q_{10} \\\\ p_{13} \\\\ q_{15} \\\\ q_{19} \\end{Bmatrix} + \\begin{Bmatrix} F\\\\ 0 \\\\ 0 \\\\ 0 \\\\ 0 \\\\ -V_{wall} \\\\ -V_{wall} \\end{Bmatrix}}\r\n\r\n\\end{matrix}\r\n\r\n\\tag{11.4}\r\n\r\n\\end{equation}\r\n

    11.4.2 \u00a0 \u00a0 \u00a0 \u00a0 Example: System Equations for an Electrical System<\/h2>\r\nFor this example, we consider the electrical system and BG model given in section 7.6<\/a> and shown in Figure 7-8<\/a>. The system has two loops. The number of state variables is three, associated with the storage elements, two inductors and one capacitor.\r\n\r\nSolution:\r\n\r\nFigure 11-4<\/a> shows the BG model with labelled power bonds for this system. The power bonds associated with storage elements, hence state variables, are colour coded.<\/a>\r\n\r\n \r\n\r\n[caption id=\"attachment_850\" align=\"aligncenter\" width=\"423\"]\"\" Figure 11-4 BG model for the electrical system with labelled power bonds[\/caption]\r\n\r\n \r\n\r\nAs mentioned, the labels are arbitrary and they merely help bookkeeping of the variables involved without affecting final solution of the system equations. The state variables are: $p_2,p_5,$ associated with the momenta of the inductors, and $q_4$, electrical charge of the capacitor. The system equations are ODEs consisting of these variables. To extract the system equations, we ask two key questions (see section 11.4<\/a>) and list their corresponding answers, using the power bond labels and causality assignments in terms of BG symbols e and f, as follows:\r\n\r\nQ1: What does each component\/element send to the system?<\/em>\r\n\r\nThe inputs to the system are from elements at the boundary of the system associated with bonds\r\n1, 2, 3, 4, 5, and 6, as listed in Table 11-3.\r\n
    \r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n\r\n
    Table 11-3 Inputs to the electrical system by all elements\u2014Q1 question<\/caption>\r\n
    Power bond label<\/strong><\/td>\r\nSystem input\/relation<\/strong><\/td>\r\n<\/tr>\r\n<\/thead>\r\n
    1<\/td>\r\n$e_1 = V$<\/td>\r\n<\/tr>\r\n
    2<\/td>\r\n$f_2 = \\dot q_2 = p_2\/L_2$<\/td>\r\n<\/tr>\r\n
    3<\/td>\r\n$e_3 = R_2f_3$<\/td>\r\n<\/tr>\r\n
    4<\/td>\r\n$e_4 = q_4\/c$<\/td>\r\n<\/tr>\r\n
    5<\/td>\r\n$f_5 = \\dot q_5 = p_5\/L_1$<\/td>\r\n<\/tr>\r\n
    6<\/td>\r\n$e_6 = R_1f_6$<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<\/div>\r\nNext, we list the relations to answer the second question:\r\n\r\nQ2: What does the system send back to the storage components\/elements?<\/em>\r\n\r\nHere, we are only interested in storage elements, as shown by colour-coded bonds in Figure 11-4<\/a>. For example, considering power bond number 2, the system sends an effort signal (or rate of momentum\/flux linkage) to the $I$-element representing the inductor $L_2$. Therefore, we can write $e_2 = \\dot p_2 = L_2 \\dot f_2$. \u00a0Similarly, considering power bond number 4, the system sends a flow signal (i.e., current or rate of charge) to the $C$-element representing the capacitor $c$. Therefore, we can write $f_4 = \\dot q_4$. Table 11-4 lists all relations corresponding to storage elements.\r\n\r\n\r\n\r\n\r\n\r\n\r\n
    Table 11-4 Outputs from the electrical system to all storage elements\u2014Q2 question<\/caption>\r\n
    Power bond label (element)<\/strong><\/td>\r\nStorage element\/relation<\/strong><\/td>\r\nRelations, using constraints<\/strong><\/td>\r\n<\/tr>\r\n<\/thead>\r\n
    2 ($I:L_2$)<\/td>\r\n$e_2 = \\dot p_2$<\/td>\r\n[latex]e_2 = \\dot p_2 = e_9 - e_3 = e_7 - e_3 = e_1 - e_6 - e_3 = V - R_1 \\underbrace{f_6}_{f_8 + f_9} - R_2 \\underbrace{f_3}_{=f_2} = V - R_1 \\left( \\dfrac{p_5}{L_1} + \\dfrac{p_2}{L_2} \\right) - R_2 \\dfrac{p_2}{L_2}[\/latex]<\/td>\r\n<\/tr>\r\n
    4 ($C:c$)<\/td>\r\n$f_4 = \\dot q_4$<\/td>\r\n$f_4 = \\dot q_4 = f_5 = \\dfrac{p_5}{L_1}$<\/td>\r\n<\/tr>\r\n
    5 ($I:L_1$)<\/td>\r\n$e_5 = \\dot p_5$<\/td>\r\n$e_5 = \\dot p_5 = \\underbrace{e_8}_{=e_7} - e_4 = e_1 - e_6 - \\dfrac{q_4}{c} = V- R_1 \\underbrace{f_6}_{f_8+f_9} - \\dfrac{q_4}{c} = V - R_1 \\left( \\dfrac{p_5}{L_1} + \\dfrac{p_2}{L_2} \\right) - \\dfrac{q_4}{c}$<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nNote that $f_6 = f_7 = f_8 + f_9 = f_5 + f_2 - \\dfrac{p_5}{L_1} + \\dfrac{p_2}{L_2}$. In the relations listed in the third column of Table 11-4, we use the constraints resulted from 1- and 0-junctions and elements\u2019 constitutive equations, and included the relations from Table 11-3. The final relations are given in bold<\/strong>. Finally, we can write the system equations, after some manipulations, in terms of three state variables in matrix form, as given in Equation (11.5).\r\n\r\n\\begin{equation*}\r\n\r\n \r\n\r\n\\begin{Bmatrix} \\dot p_2\\\\ \\dot q_4 \\\\ \\dot p_5\u00a0 \\end{Bmatrix} = \\begin{bmatrix} - \\left( \\dfrac{R_1}{L_2} + \\dfrac{R_2}{L_1} \\right) & 0 & - \\dfrac{R_2}{L_1}\\\\ 0 & 0 & \\dfrac{1}{L_1} \\\\ - \\dfrac{R_1}{L_2} & - \\dfrac{1}{c} & - \\dfrac{R_1}{L_1} \\end{bmatrix} \\begin{Bmatrix} p_2 \\\\ q_4 \\\\ p_5 \\end{Bmatrix} + \\begin{Bmatrix} V \\\\ 0 \\\\ V \\end{Bmatrix}\r\n\r\n \r\n\r\n\\tag{11.5}\r\n\r\n\\end{equation*}\r\n

    11.5 \u00a0 \u00a0 \u00a0 \u00a0 Derivative Causality and Algebraic Loop: Implicit System Equations<\/h1>\r\nWhen assigning causality strokes to a BG model, three scenarios may occur [20]<\/strong><\/a>, [21]<\/strong><\/a>. These scenarios are:\r\n
      \r\n \t
    1. all causality strokes are possible to be assigned as preferred integral causalities;<\/li>\r\n \t
    2. at least one storage element ($I$- or $C$- element) cannot be assigned with integral\r\ncausality; instead derivative causality<\/em> is forced upon the element; or<\/li>\r\n \t
    3. there is more than one option for having a BG model with assigned integral causalities for the given system, a so-called algebraic loop<\/em>.<\/li>\r\n<\/ol>\r\nAll of these scenarios are legitimate in terms of BG modelling rules, and with red colour-coded causality strokes, 20-sim provides warnings, not errors, for scenarios 2 and 3 as listed above. As mentioned, scenario 1 is desirable and preferred. For this scenario, the system equations can be derived explicitly in terms of state variables as a system of coupled first-order ODEs. The solution of the system of equations provides the answers for the generalized momenta and displacements associated with storage elements. After having the solutions for system equations, we can calculate other desired quantities related to the system.\r\n\r\nFor the second scenario, the system equations are not independent, i.e., implicit. In other words, the equations related to the storage elements with integral causalities are independent; pertinent state variables can be uniquely calculated by solving these equations, simultaneously. But the equation related to the storage element with derivative causality is not independent; it relates algebraically to the independent equations. Therefore, the related state variable for the element with derivative causality can be calculated using the solution of those of the independent variables. The mathematical dependence of the equation for the derivatively causalled element indicates that the dynamics of the element are related to and defined by the other storage elements in the system. In addition, for the derivative causality case, we may have to force more than one element having non-integral causalities.\r\n\r\nFor the third scenario, after assigning integral causalities for some of the elements in a definite manner, we face one or more options and have to make choices for one or more elements and arbitrarily assign them with causalities. The resulting system equations are again independent coupled ODEs and can be simultaneously solved to find their solutions. However, during the manipulation, we encounter an implicit equation (or equations) for a particular state variable (or its derivative) which is a function of inputs and itself as well\u2014an algebraic loop. For linear systems\/elements, the algebraic-loop situation does not pose a problem; however, for non-linear systems\/components, this makes it more difficult to manipulate the equations to find their solutions.\r\n\r\nNote that an algebraic loop may occur at several levels and, hence, make it harder to manipulate the equations. This is when, after selecting an assigning causality stroke for an element arbitrarily, further selection(s) should be made to proceed.\r\n\r\nFor both scenarios, derivative causalities and algebraic loop, the fact that they appear in a BG model is useful information about the system and\/or assumptions made for the model, even before solving the resulting equations. In many cases, the modeller can improve the system with additional components, usually of the types of $R$- element.\r\n\r\nAs mentioned, 20-sim sends warning messages, but no error messages when derivative causalities and algebraic loop exist. The resulted system equations can be solved by solvers available, and the system simulation and design can proceed as usual.\r\n\r\nIn further sections, we present some examples to demonstrate the derivative causality and algebraic-loop cases.\r\n

      11.5.1 \u00a0 \u00a0 \u00a0 \u00a0 Example: BG Model with Derivative Causality<\/h2>\r\nWe consider the BG model shown in Figure 11-5<\/a>. The colour-coded bonds associated with two inertia $I$-elements are forced, having non-integral causalities or derivative causality. This can be vividly explored by starting assigning causality from $S_f$-element sending a flow signal to the adjacent 1-junction. Therefore, all other connecting bonds to this 1-junction should have their causality strokes at the ports connecting to this 1-junction. This requirement forces the associated $I$-element to have a derivative causality, or sending effort signal, instead of receiving it for being integrally causalled, to the 1-junction. Similarly, the $I$-element $I_1$ is forced to have derivative causality. As a result, the momenta of these two $I$-elements (i.e. $I$ and $I_1$ ) will depend on the remaining storage elements, i.e., $I$-element $I_2$ and $C$-element $c$, and can be calculated having the corresponding solutions for $I_2$- and $C$- elements. For further reading, consult with references cited as [5]<\/strong><\/a>, [20]<\/strong><\/a>, [21]<\/strong><\/a>, and [30]<\/strong><\/a>.\r\n\r\nAs mentioned, in 20-sim a warning message appears on the screen when derivative causality exists in the BG model. But the software solver tools take care of this and solve the system equations. This capability is useful for practical applications in design of the systems.<\/a>\r\n\r\n[caption id=\"attachment_851\" align=\"aligncenter\" width=\"488\"]\"\" Figure 11-5 A BG model with derivative causality, colour coded[\/caption]\r\n

      11.5.2 \u00a0 \u00a0 \u00a0 \u00a0 Example: BG Model with Algebraic Loop<\/h2>\r\nWe consider the BG model shown in Figure 11-6<\/a>. Starting with the process of assigning causalities, we realize that the colour-coded bonds can\u2019t be integrally assigned and present options. In other words, to proceed further with assigning causalities to the whole system\u2019s BG model, we require to make a selection arbitrarily. Usually, the selection can be more easily made with $R$-elements, since it can be assigned with neutral causality. For example, if we select the element $R:R3$ and assign causality stroke to it such that it sends a flow signal to the adjacent 1-junction, then the remaining bonds\u2019 causality strokes can be assigned, as shown colour coded in Figure 11-7<\/a>. Alternatively, element $R:R2$ or $R:R1$ could be selected as an option. The corresponding BG models and related assigned causalities are shown in Figure 11-8<\/a> and Figure 11-9<\/a>, respectively.\r\nFor further reading, consult with references cited as [5]<\/strong><\/a>, [20]<\/strong><\/a>, [21]<\/strong><\/a>, and [30]<\/strong><\/a>.<\/a>\r\n\r\n[caption id=\"attachment_852\" align=\"aligncenter\" width=\"418\"]\"\" Figure 11-6 A BG model with algebraic-loop causality[\/caption]\r\n\r\n \r\n\r\n[caption id=\"attachment_853\" align=\"aligncenter\" width=\"418\"]\"\" Figure 11-7 The BG model with removed algebraic loop\u2014selecting R3[\/caption]\r\n\r\n \r\n\r\n[caption id=\"attachment_854\" align=\"aligncenter\" width=\"392\"]\"\" Figure 11-8 The BG model with removed algebraic loop\u2014selecting R2[\/caption]\r\n\r\n \r\n\r\n[caption id=\"attachment_855\" align=\"aligncenter\" width=\"390\"]\"\" Figure 11-9 The BG model with removed algebraic loop\u2014selecting R1[\/caption]\r\n\r\nAs mentioned, in 20-sim, a warning message appears on the screen when algebraic loop causality exists in the BG model. But the software solver tools take care of this and solve the system equations. This capability is useful for practical applications in design of the systems.\r\n

      11.6 \u00a0 \u00a0 \u00a0 \u00a0 Thermal Systems: Pseudo Bond Graph<\/h1>\r\nThermal systems are unique in terms of BG modelling. They are closely analogous to electrical systems. In practice, to analyze thermal systems, we use this analogy to build models as thermal networks, similar to electrical networks, [39]<\/strong><\/a>, [40]<\/strong><\/a>.\r\n\r\nConsidering heat conduction through a solid slab, e.g., the temperature difference $\\Delta T$ (analogous to voltage) and heat flow rate $\\dot Q_{th}$ (analogous to electrical current) satisfy $\\Delta T = R_{th} \\dot Q_{th}$. The quantity $R_{th} = \\dfrac{\\Delta L}{k_TA}$ (analogous to electrical resistance) is thermal resistance with dimension [K\/W], in $SI$, where, $\\Delta L$ is the length across which the $\\Delta T$ applies; $A$ is the cross-sectional area for heat flow, and $k_T$ is the thermal conductivity of the material. Intuitively, this analogy suggests to consider temperature as the effort and the heat flow rate as the flow when building BG models for thermal systems. With this assumption at hand, the product ($T \\dot Q_{th}$) would be considered as power variables, or the dimension of this quantity should be as of that for power. However, close examination of the dimension for the product ($T \\dot Q_{th}$) gives its dimension as $[(T\\dot Q_{th})] = K \\cdot W$, or $Kelvin \\times Watt$. Therefore, according to BG modelling rules, we cannot accept temperature and heat flow rate as conjugate power variables. This discrepancy leads us to call such a BG model a pseudo bond graph<\/em> (pseudo BG), for which temperature is the effort and heat flow rate is the flow. This is acceptable as far as the thermal system is the only system involved and the modeller is aware of the fact that the \u201cpower\u201d variables are not defined fully and correctly in the related thermal pseudo BG model. However, for multi-domain systems, the exchange of power between a thermal sub-system part and other domains of the system becomes problematic.\r\n\r\nThe solution to this challenge comes from the second law of thermodynamics, which guides us to consider the flow in the BG model as the entropy rate, instead of heat flow rate. Recall that entropy is a state function, and for a reversible system, we have the relationship $\\Delta Q_{th} = T \\Delta S$ for exchange of heat $\\Delta Q_{th}$ and entropy $\\Delta S$ using temperature $T$, [41]<\/strong><\/a>. Therefore, the heat flow rate is $\\Delta \\dot Q_{th} = T \\Delta \\dot S$. The dimension of $\\Delta \\dot S$, the entropy rate, is then [$W\/K$]. Therefore the dimension of the power variables, or ($T \\Delta \\dot S$) is [$W$], as it should be; hence, we can accept the temperature as the effort and entropy rate as the flow (or entropy as displacement) in BG model for a thermal system. Table 11-5 shows the variables involved in Pseudo BG and BG models for thermal systems.\r\n\r\n\r\n\r\n\r\n\r\n
      Table 11-5 Definitions for effort and flow for BG and pseudo BG models, thermal systems<\/caption>\r\n
      Model<\/th>\r\nEffort<\/strong><\/td>\r\nFlow<\/strong><\/td>\r\nPower<\/strong><\/td>\r\nSystem<\/strong><\/td>\r\n<\/tr>\r\n<\/thead>\r\n
      Pseudo BG<\/th>\r\ntemperature<\/td>\r\nheat flow rate<\/td>\r\nill-defined<\/td>\r\nsingle thermal domain<\/td>\r\n<\/tr>\r\n
      BG<\/th>\r\ntemperature<\/td>\r\nentropy rate<\/td>\r\ndefined<\/td>\r\n\u00a0 \u00a0multi-domain<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<\/div>","rendered":"
      \n

      11.1\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Overview<\/a><\/h1>\n

      This chapter discusses several topics pertinent to modelling and simulation of systems dynamics with focus on BG method. The reader can consult with these sections as independent topics or as supplementary to previously discussed topics.<\/p>\n

      11.2\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Energy and Power Conjugate Variables<\/h1>\n

      The product of conjugate variables related to a physical quantity should, by definition, give the corresponding quantity. In chapter 3, sections 3.2<\/a> and 3.5<\/a>, we discussed the power variables (effort \"e\" and flow \"f\") and their relations to the state variables (momentum \"p\" and displacement \"q\"). The product (\"e.f\") is power\u2014so-called power variables\u2014and integrating them, with respect to time, gives \"p\" and \"q\", respectively. Hence the \"p\" and \"q\" are called energy variables. However, further attention to the dimension of the product of the variables, i.e., (\"p.q\"), indicates that the dimension of this quantity is not equivalent to that of energy. For example, in the SI system of units, we get \"[p.q] = [N.s.m] = [J.s]\". Therefore, we can call \"e\" and \"f\" conjugate power variables because the product (\"e.f\") is power. But this definition does not apply to the product (\"p.q\"). Instead, the product (\"\dot p . q\") or product (\"p. \dot q\") can be defined as true conjugate energy variables.<\/p>\n

      One can ask the question, what is the relationship among these variables? To investigate further and provide a possible answer, we consider the system energy function \"\mathcal{E} = \mathcal{E} (p,q)\". This functional form is legitimate, since we know that the energy stored, in storage elements \"I\" and \"C\", is uniquely defined by state variables (see section 3.5<\/a>). The total differential change of energy is then \"d \mathcal{E} = \dfrac{\partial \mathcal{E}}{\partial p} dp + \dfrac{\partial \mathcal{E}}{\partial q} dq\". In principle, \"\mathcal{E} = \mathcal{E}(p,q,t)\". But the partial derivative \"\dfrac{\partial \mathcal{E}}{\partial t} |_{p,q}\" can be ignored at equilibrium state of the system, or when \"p\" and \"q\" are kept constant and energy is not changing with time[1]<\/sup><\/a>. Now, the total differential change of power \"d \mathcal{P}\" can be written as the time derivative of energy, or \"d \mathcal{P} = \dfrac{d \mathcal{E}}{dt} = \dfrac{\partial \mathcal{E}}{\partial p} \dfrac{dp}{dt} + \dfrac{\partial \mathcal{E}}{\partial q} \dfrac{dq}{dt}\". But \"\dfrac{dp}{dt} = e\" and \"\dfrac{dq}{dt} = f\", the effort and flow in BG method, respectively. Therefore, we can write<\/p>\n

      (11.1) <\/span>   <\/span>\"\begin{equation*}<\/p>\n

      Equation (11.1) clearly defines the relationship among power variables (effort \"e\" and flow \"f\") and state variables (momentum \"p\" and displacement \"q\") with respect to the power (\"\mathcal{P}\") and energy (\"\mathcal{E}\") of the system. Note that for a system, we can write this equation for each component and sum it up to get the relation for the system.<\/p>\n

      From Equation (11.1), we can conclude that by dimensional analysis, \"\dfrac{\partial \mathcal{E}}{\partial p}\" has the dimension of flow and \"\dfrac{\partial \mathcal{E}}{\partial q}\" the dimension of effort, according to BG terminology. We define these quantities as system flow and effort, or (note that the partial derivative definition is explicitly shown here)<\/p>\n

      (11.2) <\/span>   <\/span>\"\begin{equation*}<\/p>\n

      For example, considering a mass, represented by an \"I\"-element, its kinetic energy can be written as \"\dfrac{p^2}{2m}\", but \"\dfrac{\partial}{\partial p} \left( \dfrac{p^2}{2m} \right) = \dfrac{p}{m} = \dfrac{mv}{m} = v\", the velocity of the mass or the flow associated with the \"I\"-element. Similarly, considering a mechanical spring with stiffness \"k\", represented by a \"C\"-element, its potential energy can be written as \"k \frac{q^2}{2}\", but \"\dfrac{\partial}{\partial q} \left( k \dfrac{q^2}{2} \right) = kq = kx = F\", the force acting on the spring or the effort associated with the \"C\"-element. Equation (11.2) applies analogously to mechanical, electrical, hydraulic, or thermal systems. Note that the energy function \"\mathcal{E}\" relates to the Hamiltonian of the system, [13]<\/strong><\/a>.<\/p>\n

      11.3\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Including Gravity in BG Models<\/h1>\n

      Forces due to gravitational acceleration \"g\" can be implemented in a BG model using effort source element \"S_e\". But it feeds power to the system when the assigned positive direction of displacement\/velocity and the gravitational acceleration vector are in the same direction. Otherwise, the system feeds power into the gravity field, or loses power. In manually drawn BG, the direction of the power bond is drawn from \"S_e\" toward the system when the gravity and displacement directions are matched; otherwise, the power bond should be toward the \"S_e\" element. In 20-sim, however, it is not possible to draw a power bond pointing toward an \"S_e\" element because it is defined as a source element. Hence, a negative value of \"g\" is assigned to the associated parameter in the corresponding equation model, i.e., \u22129.81 \"m/s^2\" (see Figure 11-1<\/a>). Note that \"g\" is interpreted as gravitational force per unit mass, here.\u00a0<\/a><\/p>\n

      \"\"<\/a>
      Figure 11-1 Implementing gravity force for a BG model in 20-sim.<\/figcaption><\/figure>\n

      <\/a>11.4 \u00a0 \u00a0 \u00a0 \u00a0 Extracting System Equations from BG Models<\/h1>\n

      Although graphics of BG models can give us some insights into the dynamics of systems before actually solving the systems equations, having the mathematical model consisting of a system of equations is the ultimate goal resulting from a BG model, or for that matter, any model. In the previous sections, we discussed that for a BG model with integral causalities, we can derive a system of ODEs in terms of state variables, i.e., \"p\" and \"q\" (see section 3.5<\/a>). For such a system, the number of equations, \"n\" is equal to the sum of number of \"p\" associated with \"I\"– elements and number of \"q\" associated with the \"C\"– elements. In other words, the system equations are coupled first-order ODEs explicitly given in terms of state variables and can be solved simultaneously using numerical schemes. The solution of these coupled equations is, obviously, easier to obtain for linear systems rather than it is for non-linear systems.<\/p>\n

      In this section we discuss the procedure and guidelines of how to extract system equations from a BG model. For now, we assume that all storage elements (i.e., \"I\"– and \"C\"– elements) in the system are assigned with integral causality strokes. In the next section, we discuss cases when derivative causality and\/or algebraic-loop situations may exist in a BG model, along with the definitions and resulting consequences for such systems equations.<\/p>\n

      Referring to the energy management of a system, we have already discussed that the energy input to a system is partially stored in storages and the rest is dissipated through dissipaters. Therefore, the storage elements are key elements for studying the dynamics of any system. This principle also applies to the BG models, and we use it to extract system equations.<\/p>\n

      The system equations can be given in a general form as \"\dot x_i = \mathcal{F}_i (x_i, u_j, t)\" where \"\dot x_i\", the vector representing time derivative of state variables, \"u_j\" , the vector representing the inputs, with indices \"i = 1, \cdots , n,\" and \"j = 1, \cdots, r\". One can expand this functional to get<\/p>\n

      (11.3) <\/span>   <\/span>\"\begin{equation*}<\/p>\n

      After building the BG model for a system, we ask the following two key questions, [21]<\/strong><\/a>:<\/p>\n

      Q1: What does each component\/element send to the system?<\/em><\/p>\n

      Q2: What does the system send back to the storage components\/elements?<\/em><\/p>\n

      Each question guides us to list the corresponding relations involving the state variables and performing some manipulations to write the system equations in terms of the state variables, using the general guideline of \"p\" on \"I\" and \"q\" on \"C\". The following steps can be used to extract the system equations:<\/p>\n

        \n
      1. Build the BG model and assign causalities, with preferred integral causalities.<\/li>\n
      2. Simplify the BG model, as far as possible.<\/li>\n
      3. Assign labels as numbers to each power bond in the BG model. The order of numbering is arbitrary.<\/li>\n
      4. Answer Q1, in terms of all efforts or flows input into the system.<\/li>\n
      5. Answer Q2, in terms of all \"p\"\u2019s on \"I\"\u2019s and all \"q\"\u2019s on \"C\"\u2019s, or momenta (or their derivatives) associated with the inertia elements and displacements (or their derivatives) associated with the capacity elements.<\/li>\n
      6. Apply constraints implemented by 1- and 0- junction elements and perform algebraic manipulations to write the system equations only in terms of state variables, as independent variables with inputs and components\u2019 properties as parameters.<\/li>\n<\/ol>\n

        Now, we present some worked-out examples to demonstrate the application of the above mentioned procedure.<\/p>\n

        11.4.1\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Example: System Equations for a Mechanical System<\/h2>\n

        For this example, we consider a mechanical system as shown in Figure 11-2<\/a>. This systems has three DOF associated with three masses. The number of state variables is seven, associated with the storage elements, three masses, and four springs. The positive displacements are considered as shown, and tension forces are positive (+T).<\/a><\/p>\n

        \"\"
        Figure 11-2 A mechanical system sketch<\/figcaption><\/figure>\n

        Solution:<\/p>\n

        We first build the BG model, either manually or using the tools available through 20-sim, with all power bonds labelled, as shown in Figure 11-3<\/a>. As mentioned, the labels are arbitrary; they merely help bookkeeping of the variables involved without affecting final solution of the system equations. Therefore, the state variables are: \"p_2, p_6, p_{13}\" associated with the momenta of the masses and \"q_4, q_{10}, q_{15}, q_{19}\" displacements of the mechanical springs. The corresponding power bonds to these elements are colour coded (see Figure 11-3<\/a>). The system equations are ODEs consisting of these state variables. To extract the system equations, we ask two key questions (see section 11.4<\/a>) and list their corresponding answers, using the power bond labels and causality assignments in terms of BG symbols \"e\" and \"f\", as follows:<\/a><\/p>\n

         <\/p>\n

        \"\"
        Figure 11-3 BG model for the mechanical system with labelled power bonds<\/figcaption><\/figure>\n

        Q1: What does each component\/element send to the system?<\/em><\/p>\n

        The inputs to the system are from elements at the boundary of the system associated with bonds 1, 2, 4, 6, 10, 11, 13, 15, 17, and 19, as listed in Table 11-1.<\/p>\n

        \n\n\n\n\n\n\n\n\n\n\n\n\n\n\n
        Table 11-1 Inputs to the mechanical system by all elements\u2014Q1 question<\/caption>\n
        Power bond label #<\/strong><\/td>\nSystem input\/relation<\/strong><\/td>\n<\/tr>\n<\/thead>\n
        1<\/td>\n\"e_1 = F\"<\/td>\n<\/tr>\n
        2<\/td>\n\"f_2 = \dot q_2 = p_2/m_3\"<\/td>\n<\/tr>\n
        4<\/td>\n\"e_4 = k_4q_4\"<\/td>\n<\/tr>\n
        6<\/td>\n\"f_6 = \dpt q_6 = p_6/m_2\"<\/td>\n<\/tr>\n
        10<\/td>\n\"e_{10} = k_3q_{10}\"<\/td>\n<\/tr>\n
        11<\/td>\n\"e_{11} = b_1f_{11}\"<\/td>\n<\/tr>\n
        13<\/td>\n\"f_{13} = \dot q_{13} = p_{13}/m_1\"<\/td>\n<\/tr>\n
        15<\/td>\n\"e_{15} = k_2q_{15}\"<\/td>\n<\/tr>\n
        17<\/td>\n\"f_{17} = V_{wall}\"<\/td>\n<\/tr>\n
        19<\/td>\n\"e_{19} = k_1q_{19}\"<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n

        Next, we list the relations to answer the second question:<\/p>\n

        Q2: What does the system send back to the storage components\/elements?<\/em><\/p>\n

        Here, we are only interested in storage elements, as shown by colour-coded bonds in Figure 11-3<\/a>. For example, considering power bond number 2, the system sends an effort signal (or rate of momentum) to the \"I\"-element representing the mass \"m_3\". Therefore, we can write \"e_2 = \dot p_2 = m_3 \dot f_2\". Similarly, considering power bond number 15, the system sends a flow signal (or rate of displacement) to the \"C\"-element representing the spring \"k_2\". Therefore, we can write \"f_{15} = \dot q_{15}\". Table 11-2 lists all relations corresponding to storage elements.<\/p>\n\n\n\n\n\n\n\n\n\n\n\n
        Table 11-2 Outputs from the mechanical system to all storage elements\u2014Q2 question<\/caption>\n
        Power bond label #(element)<\/strong><\/td>\nStorage element\/relation<\/strong><\/td>\nRelations, using constraints<\/strong><\/td>\n<\/tr>\n<\/thead>\n
        2 (\"I:m_3\")<\/td>\n\"e_2 = \dot p_2\"<\/td>\n\"e_2 = \pmb{\dot p_2} = e_1 - e_3 = e_1 - e_4 = \pmb{F - k_4q_4}\"<\/td>\n<\/tr>\n
        4 (\"C:k_4\")<\/td>\n\"f_4 = \dot q_4\"<\/td>\n\"f_4 = \pmb{\dot q_4} = f_3 - f_5 = f_2 - f_6 = \pmb{\dfrac{p_2}{m_3} - \dfrac{p_6}{m_2}}\"<\/td>\n<\/tr>\n
        6 (\"I:m_2\")<\/td>\n\"e_6 = \dot p_6\"<\/td>\n\"e_6 = \pmb{\dot p_6} = e_5 - e_7 - e_8 = \pmb{k_4q_4 - k_3 q_{10}} - b_1 f_{11} - \pmb{k_1 q_{19}}\"<\/td>\n<\/tr>\n
        10 (\"C:k_3\")<\/td>\n\"f_{10} = \dot q_{10}\"<\/td>\n\"f_{10} = \pmb{\dot q_{10}} = f_9 = f_7 - f_{12} = f_6 - f_{13} = \pmb{\dfrac{p_6}{m_2} - \dfrac{p_{13}}{m_1}}\"<\/td>\n<\/tr>\n
        13 (\"I:m_1\")<\/td>\n\"e_{13} = \dot p_{13}\"<\/td>\n\"e_{13} = \pmb{\dot p_{13}} = e_{12} - e_{14} = e_9 - e_{15} = (e_{10} + e_{11}) - e_{15} = \pmb{k_3q_{10}} + \underbrace{b_1f_{11}}_{=b_1f_{10}} - \pmb{k_2q_{15} = k_3q_{10}} + b_1 \pmb{ \left( \dfrac{p_6}{m_2} - \dfrac{p_{13}}{m_1} \right) - k_2q_{15}}\"<\/td>\n<\/tr>\n
        15 (\"C:k_2\")<\/td>\n\"f_{15} = \dot q_{15}\"<\/td>\n\"f_{15} = \pmb{\dot q_{15}} = f_{14} - f_{16} = f_{13} - f_{17} = \pmb{\dfrac{p_{13}}{m_1} - V_{wall}}\"<\/td>\n<\/tr>\n
        19 (\"C:k_1\")<\/td>\n\"f_{19} = \dot q_{19}\"<\/td>\n\"f_{19} = \pmb{\dot q_{19}} = f_8 - f_{18} = f_6 - f_{17} = \pmb{\dfrac{p_6}{m_2} - V_{wall}}\"<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

        In the relations listed in the third column of Table 11-2, we used the constraints resulted from 1- and 0-junctions and elements\u2019 constitutive equations, and included the relations from Table 11-1. The final relations are given in bold. Finally, after some manipulations, we can write the system equations in terms of seven state variables in matrix form, as given in Equation (11.4).<\/p>\n

        (11.4) <\/span>   <\/span>\"\begin{equation*}<\/p>\n

        11.4.2 \u00a0 \u00a0 \u00a0 \u00a0 Example: System Equations for an Electrical System<\/h2>\n

        For this example, we consider the electrical system and BG model given in section 7.6<\/a> and shown in Figure 7-8<\/a>. The system has two loops. The number of state variables is three, associated with the storage elements, two inductors and one capacitor.<\/p>\n

        Solution:<\/p>\n

        Figure 11-4<\/a> shows the BG model with labelled power bonds for this system. The power bonds associated with storage elements, hence state variables, are colour coded.<\/a><\/p>\n

         <\/p>\n

        \"\"
        Figure 11-4 BG model for the electrical system with labelled power bonds<\/figcaption><\/figure>\n

         <\/p>\n

        As mentioned, the labels are arbitrary and they merely help bookkeeping of the variables involved without affecting final solution of the system equations. The state variables are: \"p_2,p_5,\" associated with the momenta of the inductors, and \"q_4\", electrical charge of the capacitor. The system equations are ODEs consisting of these variables. To extract the system equations, we ask two key questions (see section 11.4<\/a>) and list their corresponding answers, using the power bond labels and causality assignments in terms of BG symbols e and f, as follows:<\/p>\n

        Q1: What does each component\/element send to the system?<\/em><\/p>\n

        The inputs to the system are from elements at the boundary of the system associated with bonds
        \n1, 2, 3, 4, 5, and 6, as listed in Table 11-3.<\/p>\n

        \n\n\n\n\n\n\n\n\n\n\n
        Table 11-3 Inputs to the electrical system by all elements\u2014Q1 question<\/caption>\n
        Power bond label<\/strong><\/td>\nSystem input\/relation<\/strong><\/td>\n<\/tr>\n<\/thead>\n
        1<\/td>\n\"e_1 = V\"<\/td>\n<\/tr>\n
        2<\/td>\n\"f_2 = \dot q_2 = p_2/L_2\"<\/td>\n<\/tr>\n
        3<\/td>\n\"e_3 = R_2f_3\"<\/td>\n<\/tr>\n
        4<\/td>\n\"e_4 = q_4/c\"<\/td>\n<\/tr>\n
        5<\/td>\n\"f_5 = \dot q_5 = p_5/L_1\"<\/td>\n<\/tr>\n
        6<\/td>\n\"e_6 = R_1f_6\"<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n

        Next, we list the relations to answer the second question:<\/p>\n

        Q2: What does the system send back to the storage components\/elements?<\/em><\/p>\n

        Here, we are only interested in storage elements, as shown by colour-coded bonds in Figure 11-4<\/a>. For example, considering power bond number 2, the system sends an effort signal (or rate of momentum\/flux linkage) to the \"I\"-element representing the inductor \"L_2\". Therefore, we can write \"e_2 = \dot p_2 = L_2 \dot f_2\". \u00a0Similarly, considering power bond number 4, the system sends a flow signal (i.e., current or rate of charge) to the \"C\"-element representing the capacitor \"c\". Therefore, we can write \"f_4 = \dot q_4\". Table 11-4 lists all relations corresponding to storage elements.<\/p>\n\n\n\n\n\n\n\n
        Table 11-4 Outputs from the electrical system to all storage elements\u2014Q2 question<\/caption>\n
        Power bond label (element)<\/strong><\/td>\nStorage element\/relation<\/strong><\/td>\nRelations, using constraints<\/strong><\/td>\n<\/tr>\n<\/thead>\n
        2 (\"I:L_2\")<\/td>\n\"e_2 = \dot p_2\"<\/td>\n\"e_2 = \dot p_2 = e_9 - e_3 = e_7 - e_3 = e_1 - e_6 - e_3 = V - R_1 \underbrace{f_6}_{f_8 + f_9} - R_2 \underbrace{f_3}_{=f_2} = V - R_1 \left( \dfrac{p_5}{L_1} + \dfrac{p_2}{L_2} \right) - R_2 \dfrac{p_2}{L_2}\"<\/td>\n<\/tr>\n
        4 (\"C:c\")<\/td>\n\"f_4 = \dot q_4\"<\/td>\n\"f_4 = \dot q_4 = f_5 = \dfrac{p_5}{L_1}\"<\/td>\n<\/tr>\n
        5 (\"I:L_1\")<\/td>\n\"e_5 = \dot p_5\"<\/td>\n\"e_5 = \dot p_5 = \underbrace{e_8}_{=e_7} - e_4 = e_1 - e_6 - \dfrac{q_4}{c} = V- R_1 \underbrace{f_6}_{f_8+f_9} - \dfrac{q_4}{c} = V - R_1 \left( \dfrac{p_5}{L_1} + \dfrac{p_2}{L_2} \right) - \dfrac{q_4}{c}\"<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

        Note that \"f_6 = f_7 = f_8 + f_9 = f_5 + f_2 - \dfrac{p_5}{L_1} + \dfrac{p_2}{L_2}\". In the relations listed in the third column of Table 11-4, we use the constraints resulted from 1- and 0-junctions and elements\u2019 constitutive equations, and included the relations from Table 11-3. The final relations are given in bold<\/strong>. Finally, we can write the system equations, after some manipulations, in terms of three state variables in matrix form, as given in Equation (11.5).<\/p>\n

        (11.5) <\/span>   <\/span>\"\begin{equation*}<\/p>\n

        11.5 \u00a0 \u00a0 \u00a0 \u00a0 Derivative Causality and Algebraic Loop: Implicit System Equations<\/h1>\n

        When assigning causality strokes to a BG model, three scenarios may occur [20]<\/strong><\/a>, [21]<\/strong><\/a>. These scenarios are:<\/p>\n

          \n
        1. all causality strokes are possible to be assigned as preferred integral causalities;<\/li>\n
        2. at least one storage element (\"I\"– or \"C\"– element) cannot be assigned with integral
          \ncausality; instead derivative causality<\/em> is forced upon the element; or<\/li>\n
        3. there is more than one option for having a BG model with assigned integral causalities for the given system, a so-called algebraic loop<\/em>.<\/li>\n<\/ol>\n

          All of these scenarios are legitimate in terms of BG modelling rules, and with red colour-coded causality strokes, 20-sim provides warnings, not errors, for scenarios 2 and 3 as listed above. As mentioned, scenario 1 is desirable and preferred. For this scenario, the system equations can be derived explicitly in terms of state variables as a system of coupled first-order ODEs. The solution of the system of equations provides the answers for the generalized momenta and displacements associated with storage elements. After having the solutions for system equations, we can calculate other desired quantities related to the system.<\/p>\n

          For the second scenario, the system equations are not independent, i.e., implicit. In other words, the equations related to the storage elements with integral causalities are independent; pertinent state variables can be uniquely calculated by solving these equations, simultaneously. But the equation related to the storage element with derivative causality is not independent; it relates algebraically to the independent equations. Therefore, the related state variable for the element with derivative causality can be calculated using the solution of those of the independent variables. The mathematical dependence of the equation for the derivatively causalled element indicates that the dynamics of the element are related to and defined by the other storage elements in the system. In addition, for the derivative causality case, we may have to force more than one element having non-integral causalities.<\/p>\n

          For the third scenario, after assigning integral causalities for some of the elements in a definite manner, we face one or more options and have to make choices for one or more elements and arbitrarily assign them with causalities. The resulting system equations are again independent coupled ODEs and can be simultaneously solved to find their solutions. However, during the manipulation, we encounter an implicit equation (or equations) for a particular state variable (or its derivative) which is a function of inputs and itself as well\u2014an algebraic loop. For linear systems\/elements, the algebraic-loop situation does not pose a problem; however, for non-linear systems\/components, this makes it more difficult to manipulate the equations to find their solutions.<\/p>\n

          Note that an algebraic loop may occur at several levels and, hence, make it harder to manipulate the equations. This is when, after selecting an assigning causality stroke for an element arbitrarily, further selection(s) should be made to proceed.<\/p>\n

          For both scenarios, derivative causalities and algebraic loop, the fact that they appear in a BG model is useful information about the system and\/or assumptions made for the model, even before solving the resulting equations. In many cases, the modeller can improve the system with additional components, usually of the types of \"R\"– element.<\/p>\n

          As mentioned, 20-sim sends warning messages, but no error messages when derivative causalities and algebraic loop exist. The resulted system equations can be solved by solvers available, and the system simulation and design can proceed as usual.<\/p>\n

          In further sections, we present some examples to demonstrate the derivative causality and algebraic-loop cases.<\/p>\n

          11.5.1 \u00a0 \u00a0 \u00a0 \u00a0 Example: BG Model with Derivative Causality<\/h2>\n

          We consider the BG model shown in Figure 11-5<\/a>. The colour-coded bonds associated with two inertia \"I\"-elements are forced, having non-integral causalities or derivative causality. This can be vividly explored by starting assigning causality from \"S_f\"-element sending a flow signal to the adjacent 1-junction. Therefore, all other connecting bonds to this 1-junction should have their causality strokes at the ports connecting to this 1-junction. This requirement forces the associated \"I\"-element to have a derivative causality, or sending effort signal, instead of receiving it for being integrally causalled, to the 1-junction. Similarly, the \"I\"-element \"I_1\" is forced to have derivative causality. As a result, the momenta of these two \"I\"-elements (i.e. \"I\" and \"I_1\" ) will depend on the remaining storage elements, i.e., \"I\"-element \"I_2\" and \"C\"-element \"c\", and can be calculated having the corresponding solutions for \"I_2\"– and \"C\"– elements. For further reading, consult with references cited as [5]<\/strong><\/a>, [20]<\/strong><\/a>, [21]<\/strong><\/a>, and [30]<\/strong><\/a>.<\/p>\n

          As mentioned, in 20-sim a warning message appears on the screen when derivative causality exists in the BG model. But the software solver tools take care of this and solve the system equations. This capability is useful for practical applications in design of the systems.<\/a><\/p>\n

          \"\"
          Figure 11-5 A BG model with derivative causality, colour coded<\/figcaption><\/figure>\n

          11.5.2 \u00a0 \u00a0 \u00a0 \u00a0 Example: BG Model with Algebraic Loop<\/h2>\n

          We consider the BG model shown in Figure 11-6<\/a>. Starting with the process of assigning causalities, we realize that the colour-coded bonds can\u2019t be integrally assigned and present options. In other words, to proceed further with assigning causalities to the whole system\u2019s BG model, we require to make a selection arbitrarily. Usually, the selection can be more easily made with \"R\"-elements, since it can be assigned with neutral causality. For example, if we select the element \"R:R3\" and assign causality stroke to it such that it sends a flow signal to the adjacent 1-junction, then the remaining bonds\u2019 causality strokes can be assigned, as shown colour coded in Figure 11-7<\/a>. Alternatively, element \"R:R2\" or \"R:R1\" could be selected as an option. The corresponding BG models and related assigned causalities are shown in Figure 11-8<\/a> and Figure 11-9<\/a>, respectively.
          \nFor further reading, consult with references cited as           [5]<\/strong><\/a>, [20]<\/strong><\/a>, [21]<\/strong><\/a>, and [30]<\/strong><\/a>.<\/a><\/p>\n

          \"\"
          Figure 11-6 A BG model with algebraic-loop causality<\/figcaption><\/figure>\n

           <\/p>\n

          \"\"
          Figure 11-7 The BG model with removed algebraic loop\u2014selecting R3<\/figcaption><\/figure>\n

           <\/p>\n

          \"\"
          Figure 11-8 The BG model with removed algebraic loop\u2014selecting R2<\/figcaption><\/figure>\n

           <\/p>\n

          \"\"
          Figure 11-9 The BG model with removed algebraic loop\u2014selecting R1<\/figcaption><\/figure>\n

          As mentioned, in 20-sim, a warning message appears on the screen when algebraic loop causality exists in the BG model. But the software solver tools take care of this and solve the system equations. This capability is useful for practical applications in design of the systems.<\/p>\n

          11.6 \u00a0 \u00a0 \u00a0 \u00a0 Thermal Systems: Pseudo Bond Graph<\/h1>\n

          Thermal systems are unique in terms of BG modelling. They are closely analogous to electrical systems. In practice, to analyze thermal systems, we use this analogy to build models as thermal networks, similar to electrical networks, [39]<\/strong><\/a>, [40]<\/strong><\/a>.<\/p>\n

          Considering heat conduction through a solid slab, e.g., the temperature difference \"\Delta T\" (analogous to voltage) and heat flow rate \"\dot Q_{th}\" (analogous to electrical current) satisfy \"\Delta T = R_{th} \dot Q_{th}\". The quantity \"R_{th} = \dfrac{\Delta L}{k_TA}\" (analogous to electrical resistance) is thermal resistance with dimension [K\/W], in \"SI\", where, \"\Delta L\" is the length across which the \"\Delta T\" applies; \"A\" is the cross-sectional area for heat flow, and \"k_T\" is the thermal conductivity of the material. Intuitively, this analogy suggests to consider temperature as the effort and the heat flow rate as the flow when building BG models for thermal systems. With this assumption at hand, the product (\"T \dot Q_{th}\") would be considered as power variables, or the dimension of this quantity should be as of that for power. However, close examination of the dimension for the product (\"T \dot Q_{th}\") gives its dimension as \"[(T\dot Q_{th})] = K \cdot W\", or \"Kelvin \times Watt\". Therefore, according to BG modelling rules, we cannot accept temperature and heat flow rate as conjugate power variables. This discrepancy leads us to call such a BG model a pseudo bond graph<\/em> (pseudo BG), for which temperature is the effort and heat flow rate is the flow. This is acceptable as far as the thermal system is the only system involved and the modeller is aware of the fact that the \u201cpower\u201d variables are not defined fully and correctly in the related thermal pseudo BG model. However, for multi-domain systems, the exchange of power between a thermal sub-system part and other domains of the system becomes problematic.<\/p>\n

          The solution to this challenge comes from the second law of thermodynamics, which guides us to consider the flow in the BG model as the entropy rate, instead of heat flow rate. Recall that entropy is a state function, and for a reversible system, we have the relationship \"\Delta Q_{th} = T \Delta S\" for exchange of heat \"\Delta Q_{th}\" and entropy \"\Delta S\" using temperature \"T\", [41]<\/strong><\/a>. Therefore, the heat flow rate is \"\Delta \dot Q_{th} = T \Delta \dot S\". The dimension of \"\Delta \dot S\", the entropy rate, is then [\"W/K\"]. Therefore the dimension of the power variables, or (\"T \Delta \dot S\") is [\"W\"], as it should be; hence, we can accept the temperature as the effort and entropy rate as the flow (or entropy as displacement) in BG model for a thermal system. Table 11-5 shows the variables involved in Pseudo BG and BG models for thermal systems.<\/p>\n\n\n\n\n\n\n
          Table 11-5 Definitions for effort and flow for BG and pseudo BG models, thermal systems<\/caption>\n
          Model<\/th>\nEffort<\/strong><\/td>\nFlow<\/strong><\/td>\nPower<\/strong><\/td>\nSystem<\/strong><\/td>\n<\/tr>\n<\/thead>\n
          Pseudo BG<\/th>\ntemperature<\/td>\nheat flow rate<\/td>\nill-defined<\/td>\nsingle thermal domain<\/td>\n<\/tr>\n
          BG<\/th>\ntemperature<\/td>\nentropy rate<\/td>\ndefined<\/td>\n\u00a0 \u00a0multi-domain<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<\/div>\n
          1. \"d \mathcal{E} (p,q,t) = \dfrac{\partial \mathcal{E}}{\partial p} dp + \dfrac{\partial \mathcal{E}}{\partial q} dq + \dfrac{\partial \mathcal{E}}{\partial t} dt\" ↵<\/a><\/li><\/ol><\/div>","protected":false},"author":801,"menu_order":11,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":[],"pb_section_license":""},"chapter-type":[],"contributor":[],"license":[],"class_list":["post-76","chapter","type-chapter","status-publish","hentry"],"part":3,"_links":{"self":[{"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/pressbooks\/v2\/chapters\/76","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/wp\/v2\/users\/801"}],"version-history":[{"count":25,"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/pressbooks\/v2\/chapters\/76\/revisions"}],"predecessor-version":[{"id":2462,"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/pressbooks\/v2\/chapters\/76\/revisions\/2462"}],"part":[{"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/pressbooks\/v2\/parts\/3"}],"metadata":[{"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/pressbooks\/v2\/chapters\/76\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/wp\/v2\/media?parent=76"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/pressbooks\/v2\/chapter-type?post=76"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/wp\/v2\/contributor?post=76"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/pressbooks.bccampus.ca\/engineeringsystems\/wp-json\/wp\/v2\/license?post=76"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}