ECE PhD QE PE1 2013 Problem1 - Rhea

Problem 1

Coenergy Calculation

Coenergy will 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).

$ \begin{equation} W_{c,0} = 0 \end{equation} $

The first step in the coenergy calculation will ramp dummy variable $ i_{1}' $ from $ 0 $ to its final value of $ i_{1} $ while dummy variable $ i_{2}' = 0 $ is held at its initial value. Recall that $ \int \frac{au}{\sqrt{b + u^2}} \, du = a \sqrt{b + u^2} + C $ where $ C $ is an arbitrary constant.

$ \begin{align} W_{c,1} &= \left.\int_{i_{1}' = 0}^{i_{1}} \lambda_{1}\left(i_{1}', i_{2}', \theta_{rm}\right) \, di_{1}'\right|_{i_{2}' = 0} + \left.\int_{i_{1}' = 0}^{i_{1}} \lambda_{2}\left(i_{1}', i_{2}', \theta_{rm}\right) \, di_{1}'\right|_{i_{2}' = 0} \\ W_{c,1} &= \int_{i_{1}' = 0}^{i_{1}} \left[10i_{1}' + 2 \left(5 + 2\sin(2\theta_{rm})\right) \frac{2i_{1}' + \cancelto{0}{i_{2}'}}{\sqrt{1 + (2i_{1}' + \cancelto{0}{i_{2}'})^2}}\right] \, di_{1}' + \int_{i_{1}' = 0}^{i_{1}} \left[5\cancelto{0}{i_{2}'} + \left(5 + 2\sin(2\theta_{rm})\right) \frac{2i_{1}' + \cancelto{0}{i_{2}'}}{\sqrt{1 + (2i_{1}' + \cancelto{0}{i_{2}'})^2}}\right] \, di_{1}' \\ W_{c,1} &= 5 i_{1}^2 + 2 \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1})^2} - 2 \left(5 + 2\sin(2\theta_{rm})\right) + \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1})^2} - \left(5 + 2\sin(2\theta_{rm})\right) \\ W_{c,1} &= 5 i_{1}^2 + 3 \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1})^2} - 3 \left(5 + 2\sin(2\theta_{rm})\right) \end{align} $

The evaluation of the previous integral has terms for $ i_{1}' = 0 $ (lower bound evaluation) that must not be forgotten. The second step in the coenergy calculation will ramp dummy variable $ i_{2}' $ from $ 0 $ to its final value of $ i_{2} $ while dummy variable $ i_{1}' = i_{1} $ is held at its final value.

$ \begin{align} W_{c,2} &= \left.\int_{i_{2}' = 0}^{i_{2}} \lambda_{1}\left(i_{1}', i_{2}', \theta_{rm}\right) \, di_{2}'\right|_{i_{1}' = i_{1}} + \left.\int_{i_{2}' = 0}^{i_{2}} \lambda_{2}\left(i_{1}', i_{2}', \theta_{rm}\right) \, di_{2}'\right|_{i_{1}' = i_{1}} \\ W_{c,2} &= \int_{i_{2}' = 0}^{i_{2}} \left[10i_{1} + 2 \left(5 + 2\sin(2\theta_{rm})\right) \frac{2i_{1} + i_{2}'}{\sqrt{1 + (2i_{1} + i_{2}')^2}}\right] \, di_{2}' + \int_{i_{2}' = 0}^{i_{2}} \left[5i_{2}' + \left(5 + 2\sin(2\theta_{rm})\right) \frac{2i_{1} + i_{2}'}{\sqrt{1 + (2i_{1} + i_{2}')^2}}\right] \, di_{2}' \\ W_{c,2} &= 10 i_{1} i_{2} + 2 \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1} + i_{2})^2} - 2 \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1})^2} + \frac{5}{2} i_{2}^2 + \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1} + i_{2})^2} - \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1})^2} \\ W_{c,2} &= 10 i_{1} i_{2} + \frac{5}{2} i_{2}^2 + 3 \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1} + i_{2})^2} - 3 \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1})^2} \end{align} $

The evaluation of the previous integral has terms for $ i_{2}' = 0 $ (lower bound evaluation) that must not be forgotten. The final step is the sum all the individual coenergy contributions $ W_c(i_{1}, i_{2}, \theta_{rm}) = \sum_{k = 0}^{2} W_{c,k}(i_{1}, i_{2}, \theta_{rm}) $.

$ \begin{align} W_c(i_{1}, i_{2}, \theta_{rm}) &= W_{c,0} + W_{c,1} + W_{c,2} \\ W_c(i_{1}, i_{2}, \theta_{rm}) &= \begin{split} &{}0 + 5 i_{1}^2 + 3 \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1})^2} - 3 \left(5 + 2\sin(2\theta_{rm})\right) \\ &{}+ 10 i_{1} i_{2} + \frac{5}{2} i_{2}^2 + 3 \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1} + i_{2})^2} - 3 \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1})^2} \\ \end{split} \\ W_c(i_{1}, i_{2}, \theta_{rm}) &= 5 i_{1}^2 + 10 i_{1} i_{2} + \frac{5}{2} i_{2}^2 + 3 \left(5 + 2\sin(2\theta_{rm})\right) \sqrt{1 + (2i_{1} + i_{2})^2} - 3 \left(5 + 2\sin(2\theta_{rm})\right) \end{align} $

Electromechanical Torque Calculation

Electromechanical torque is just the partial derivative of coenergy with respect to mechanical rotor position.

$ \begin{align} T_e &= \frac{\partial W_c(i_{1}, i_{2}, \theta_{rm})}{\partial \theta_{rm}} \\ T_e &= 0 + 0 + 0 + 12\cos(2\theta_{rm}) \sqrt{1 + (2i_{1} + i_{2})^2} - 12\cos(2\theta_{rm}) \end{align} $

Thus, the electromechanical torque equation is obtained for this device.

$ \begin{equation} \boxed{T_e = 12\cos(2\theta_{rm}) \left[\sqrt{1 + (2i_{1} + i_{2})^2} - 1\right]} \end{equation} $

Discussion

Back to PE-1, August 2013