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 of the circuit. At each node of an electric circuit one writes the currents coming in and out of the node. Take care however in the MNA 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, whereas the branch currents coming out are given a positive sign.

Step 2

Use the BCE in terms of the node voltages of the circuit to eliminate as many branch currents as possible. Writing the BCE's in terms of the node voltages saves one step. If the BCE's 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 G i,  G = 1/R, is used instead of R. We now proceed as explained above.

RC Circuit
RC Circuit
Element Branch equation
Resistor I_R = GV_R
Capacitor I_C = C\frac{dV_C}{dt}

Step 1

In this case there are two nodes, e_1 and e_2. Also there are three currents: i_{V_s},  i_{R} and i_{C}.

At node e1 the KCL yields:

i_{V_s} + i_R = 0

and at node e2:

-i_R + i_C = 0

Step 2

With the provided BCEs in the table and observing that:


V_s = e_1


V_R = e_1 - e_2


V_C = e_2,

the following equations are the result:

G(e_1 - e_2) + i_{V_S} = 0

C\frac{de_2}{dt} + G(e_2  - e_1) = 0

Step 3

Note that at this point there are two equations but three unknowns. The missing equation comes from the fact that

e_1 = V_s

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 \mathbf{x} = \begin{pmatrix}e_1&e_2&i_{V_S}\end{pmatrix}^T is defined, then the above equations can be put in the form Ex'(t) + Ax(t) = f,

where A = \begin{pmatrix}G & -G&  1\\-G & G & 0\\1 & 0 & 0\end{pmatrix}, E = \begin{pmatrix} 0 & 0 & 0\\0& C& 0\\ 0& 0& 0\end{pmatrix} and f = \begin{pmatrix}0&0&V_s\end{pmatrix}^T.

This is a linear differential algebraic equation (DAE), since E 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] 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

  1. Ho, Ruehli, and Brennan (April 1974). "The Modified Nodal Approach to Network Analysis". Proc. 1974 Int. Symposium on Circuits and Systems, San Francisco. pp. 505509.
  2. 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): 101113. doi:10.1109/TCT.1971.1083223.
  3. Cheng, Chung-Kuan. Lecture Notes for CSE245: Computer-Aided Circuit Simulation and Verification. Spring 2006. Lecture 1.
  4. Tischendorf C. Topological index of DAEs in the Circuit Simulation.
  5. 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.
  6. 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.
  7. Markus Kunze (2000). Non-Smooth Dynamical Systems. Springer Science & Business Media. ISBN 978-3-540-67993-6.

External links