Answers and Discussions for
Contents
Problem 1, part b
Recycled Information from part a
The voltage equation and flux linkage equation for these five identical, magnetically coupled inductors in matrix form is reused(resistance is neglected and linear magnetics assumed).
$ \begin{align} \vec{v}_{12345s} &= p \vec{\lambda}_{12345s} \\ \vec{\lambda}_{12345s} &= \mathbf{L}_{ss} \vec{i}_{12345s} \\ \begin{bmatrix} \lambda_1 \\ \lambda_2 \\ \lambda_3 \\ \lambda_4 \\ \lambda_5 \end{bmatrix} &= \begin{bmatrix} L_s & M_{ns} & M_{fs} & M_{fs} & M_{ns} \\ M_{ns} & L_s & M_{ns} & M_{fs} & M_{fs} \\ M_{fs} & M_{ns} & L_s & M_{ns} & M_{fs} \\ M_{fs} & M_{fs} & M_{ns} & L_s & M_{ns} \\ M_{ns} & M_{fs} & M_{fs} & M_{ns} & L_s \end{bmatrix} \begin{bmatrix} i_1 \\ i_2 \\ i_3 \\ i_4 \\ i_5 \end{bmatrix} \end{align} $
Because the inductors are identical to each other, only three unique, constant entries exist in the inductance matrix $ \mathbf{L}_{ss} $: the self-inductance $ L_s > 0 $, the mutual inductance between inductors $ \frac{2\pi}{5} \, \text{rad} $ apart $ M_{ns} > 0 $, and the mutual inductance between inductors $ \frac{4\pi}{5} \, \text{rad} $ apart $ M_{fs} < 0 $.
The results from the previous current excitation must also be retained. Based on the symmetry of the situation under this operating condition alone, it may be deduced from the flux linkage equation that $ v_5(t) = v_2(t) $ and $ v_3(t) = v_4(t) $.
$ \begin{align} \begin{bmatrix} i_1(t) \\ i_2(t) \\ i_3(t) \\ i_4(t) \\ i_5(t) \end{bmatrix} &= \begin{bmatrix} 10 \cos\left(100 \, \frac{\text{rad}}{\text{s}} t\right) \\ 0 \\ 0 \\ 0 \\ 0 \end{bmatrix} \, \text{A} \\ \begin{bmatrix} v_1(t) \\ v_2(t) \\ v_3(t) \\ v_4(t) \\ v_5(t) \end{bmatrix} &= \begin{bmatrix} 10^3 L_s \cos\left(100 \, \frac{\text{rad}}{\text{s}} t - \frac{\pi}{2} \, \text{rad}\right) \\ 10^3 M_{ns} \cos\left(100 \, \frac{\text{rad}}{\text{s}} t - \frac{\pi}{2} \, \text{rad}\right) \\ 10^3 |M_{fs}| \cos\left(100 \, \frac{\text{rad}}{\text{s}} t + \frac{\pi}{2} \, \text{rad}\right) \\ 10^3 |M_{fs}| \cos\left(100 \, \frac{\text{rad}}{\text{s}} t + \frac{\pi}{2} \, \text{rad}\right) \\ 10^3 M_{ns} \cos\left(100 \, \frac{\text{rad}}{\text{s}} t - \frac{\pi}{2} \, \text{rad}\right) \end{bmatrix} \, \text{V} = \begin{bmatrix} 40 \cos\left(100 \, \frac{\text{rad}}{\text{s}} t - \frac{\pi}{2} \, \text{rad}\right) \\ 10 \cos\left(100 \, \frac{\text{rad}}{\text{s}} t - \frac{\pi}{2} \, \text{rad}\right) \\ 1 \cos\left(100 \, \frac{\text{rad}}{\text{s}} t + \frac{\pi}{2} \, \text{rad}\right) \\ 1 \cos\left(100 \, \frac{\text{rad}}{\text{s}} t + \frac{\pi}{2} \, \text{rad}\right) \\ 10 \cos\left(100 \, \frac{\text{rad}}{\text{s}} t - \frac{\pi}{2} \, \text{rad}\right) \end{bmatrix} \, \text{V} \end{align} $
Field Energy Calculation
The inductance matrix entries are solved based on the previous results: $ L_s = +0.04 \, \text{H} $, $ M_{ns} = +0.01 \, \text{H} $, and $ M_{fs} = -0.001 \, \text{H} $. The new current excitation is written in matrix form.
$ \begin{equation} \begin{bmatrix} i_1(t) \\ i_2(t) \\ i_3(t) \\ i_4(t) \\ i_5(t) \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \\ 10 \\ 0 \\ 2 \end{bmatrix} \, \text{A} \, \forall t \end{equation} $
Coenergy will first be calculated in steps using a sequential ramping process. Before any steps can be completed, the contribution to coenergy of fixing the mechanical system should be documented (zero unless the mechanical system can store energy in its position). This coenergy contribution is trivially zero since the coupled inductor circuit lacks a mechanical system.
$ \begin{equation} W_{c,0} = 0 \end{equation} $
The first step in the coenergy calculation will ramp dummy variable $ i_{3}' $ from $ 0 $ to its final value of $ i_{3} $ while dummy variable $ i_{5}' = 0 $ is held at its initial value. Many terms in the expanded equation have zero value and will be omitted.
$ \begin{align} W_{c,1} &= \left. \sum_{j = 1}^5 \int_{i_{3}' = 0}^{i_{3}'} \lambda_{j}\left(\vec{i}_j\right) \, di_{3}'\right|_{i_{5}' = 0} \\ W_{c,1} &= \int_{i_{3}' = 0}^{i_{3}} \left[M_{fs} i_{3}'\right] \, di_{3}' + \int_{i_{3}' = 0}^{i_{3}} \left[M_{ns} i_{3}'\right] \, di_{3}' + \int_{i_{3}' = 0}^{i_{3}} \left[L_s i_{3}'\right] \, di_{3}' + \int_{i_{3}' = 0}^{i_{3}} \left[M_{ns} i_{3}'\right] \, di_{3}' + \int_{i_{3}' = 0}^{i_{3}} \left[M_{fs} i_{3}'\right] \, di_{3}' \\ W_{c,1} &= \frac{1}{2} M_{fs} i_{3}^2 + \frac{1}{2} M_{ns} i_{3}^2 + \frac{1}{2} L_s i_{3}^2 + \frac{1}{2} M_{ns} i_{3}^2 + \frac{1}{2} M_{fs} i_{3}^2 \end{align} $
