Engineering:Modified nodal analysis
In electrical engineering, modified nodal analysis[1] or MNA is an extension of nodal analysis which not only determines the circuit's node voltages (as in classical nodal analysis), but also some branch currents. Modified nodal analysis was developed as a formalism to mitigate the difficulty of representing voltage-defined components in nodal analysis (e.g. voltage-controlled voltage sources). It is one such formalism. Others, such as sparse tableau formulation,[2] are equally general and related via matrix transformations.
Method
The MNA uses the element's branch constitutive equations or BCE, i.e., their voltage - current characteristic and the Kirchhoff's circuit laws. The method is often done in four steps,[3] but it can be reduced to three:
Step 1
Write the KCL equations of the circuit. At each node of an electric circuit, write the currents coming into and out of the node. Take care, however, in the MNA method, the current of the independent voltage sources is taken from the "plus" to the "minus" (see Figure 1). Also, note that the right hand side of each equation is always equal to zero, so that the branch currents that come into the node are given a negative sign and those that go out are given a positive sign.
Step 2
Use the BCEs in terms of the node voltages of the circuit to eliminate as many branch currents as possible. Writing the BCEs in terms of the node voltages saves one step. If the BCEs were written in terms of the branch voltages, one more step, i.e., replacing the branches voltages for the node ones, would be necessary. In this article the letter "e" is used to name the node voltages, while the letter "v" is used to name the branch voltages.
Step 3
Finally, write down the unused equations.
Example
The figure shows a RC series circuit and the table shows the BCE of a linear resistor and a linear capacitor. Note that in the case of the resistor the admittance [math]\displaystyle{ G }[/math] i, [math]\displaystyle{ G = 1/R }[/math], is used instead of [math]\displaystyle{ R }[/math]. We now proceed as explained above.
| Element | Branch equation |
|---|---|
| Resistor | [math]\displaystyle{ I_R = GV_R }[/math] |
| Capacitor | [math]\displaystyle{ I_C = C\frac{dV_C}{dt} }[/math] |
Step 1
In this case there are two nodes, [math]\displaystyle{ e_1 }[/math] and [math]\displaystyle{ e_2 }[/math]. Also there are three currents: [math]\displaystyle{ i_{V_s} }[/math], [math]\displaystyle{ i_{R} }[/math] and [math]\displaystyle{ i_{C} }[/math].
At node e1 the KCL yields:
[math]\displaystyle{ i_{V_s} + i_R = 0 }[/math]
and at node e2:
[math]\displaystyle{ -i_R + i_C = 0 }[/math]
Step 2
With the provided BCEs in the table and observing that:
[math]\displaystyle{ V_s = e_1 }[/math]
[math]\displaystyle{ V_R = e_1 - e_2 }[/math]
[math]\displaystyle{ V_C = e_2, }[/math]
the following equations are the result:
[math]\displaystyle{ G(e_1 - e_2) + i_{V_S} = 0 }[/math]
[math]\displaystyle{ C\frac{de_2}{dt} + G(e_2 - e_1) = 0 }[/math]
Step 3
Note that at this point there are two equations but three unknowns. The missing equation comes from the fact that
[math]\displaystyle{ e_1 = V_s }[/math]
and so finally we have three equations and three unknowns, that leads to a solvable linear system.
Modified Nodal Analysis and DAEs
If the vector [math]\displaystyle{ \mathbf{x} = \begin{pmatrix}e_1&e_2&i_{V_S}\end{pmatrix}^T }[/math] is defined, then the above equations can be put in the form [math]\displaystyle{ Ex'(t) + Ax(t) = f, }[/math]
where [math]\displaystyle{ A = \begin{pmatrix}G & -G& 1\\-G & G & 0\\1 & 0 & 0\end{pmatrix} }[/math], [math]\displaystyle{ E = \begin{pmatrix} 0 & 0 & 0\\0& C& 0\\ 0& 0& 0\end{pmatrix} }[/math] and [math]\displaystyle{ f = \begin{pmatrix}0&0&V_s\end{pmatrix}^T }[/math].
This is a linear differential algebraic equation (DAE), since [math]\displaystyle{ E }[/math] is singular. It can be proved that such a DAE coming from the Modified Nodal Analysis will have differentiation index less or equal than two as long as only passive RLC components are used.[4][full citation needed] When using active components, such as operational amplifiers, the differentiation index can be arbitrarily high.[5]
Non-smooth analysis
DAEs assume smooth characteristics for individual components; for example, a diode can be modeled/represented in a MNA with DAEs via the Shockley equation, but one cannot use an apparently simpler (more ideal) model where the sharply exponential forward and breakdown conduction regions of the curve are just straight vertical lines. Circuit analysis (including MNA) with the latter kind of equations is actually more involved (than using DAEs) and is the topic of non-smooth dynamical systems (NSDS) analysis, which relies on the theory of differential inclusions.[6][7]
References
- ↑ Ho, Ruehli, and Brennan (April 1974). "The Modified Nodal Approach to Network Analysis". pp. 505–509. doi:10.1109/TCS.1975.1084079.
- ↑ Hachtel, G., Brayton, R, and Gustavson, F. (January 1971). "The Sparse Tableau Approach to Network Analysis and Design". IEEE Transactions on Circuit Theory 18 (1): 101–113. doi:10.1109/TCT.1971.1083223.
- ↑ Cheng, Chung-Kuan. Lecture Notes for CSE245: Computer-Aided Circuit Simulation and Verification. Spring 2006. Lecture 1.
- ↑ Tischendorf C. Topological index of DAEs in the Circuit Simulation.
- ↑ K. E. Brenan; S. L. Campbell; L. R. Petzold (1996). Numerical Solution of Initial-value Problems in Differential-algebraic Equations. SIAM. pp. 173–177. ISBN 978-1-61197-122-4.
- ↑ Vincent Acary; Olivier Bonnefon; Bernard Brogliato (2010). Nonsmooth Modeling and Simulation for Switched Circuits. Springer Science & Business Media. pp. 3–4 (for the diode example). ISBN 978-90-481-9681-4.
- ↑ Markus Kunze (2000). Non-Smooth Dynamical Systems. Springer Science & Business Media. ISBN 978-3-540-67993-6. https://archive.org/details/nonsmoothdynamic0000kunz.
External links
- "Modified Nodal Analysis (DC algorithm description in Qucs technical documentation)". http://qucs.sourceforge.net/tech/node14.html.
 |  |
