次ぎもひねった問題。Rの値によらずE0とEの位相が変わらない条件を導けというもの。



以下の関係が成り立つ

(R+jωL1)*I1+jωM*I2=E0

jωL2*I2+jωM*I1=E0

R*I1=E

K*E=E0

これらをI1,I2,E,Kに関する４元連立方程式として解くと

(%i54) e1:(R+%i*o*L1)*I1+%i*o*M*I2=E0;
(%o54) I1*(R+%i*o*L1)+%i*o*I2*M=E0
(%i55) e2:%i*o*L2*I2+%i*o*M*I1=E0;
(%o55) %i*o*I1*M+%i*o*I2*L2=E0
(%i56) e3:R*I1=E;
(%o56) I1*R=E
(%i57) e4:K*E=E0;
(%o57) E*K=E0
(%i58) solve([e1,e2,e3,e4],[I1,I2,E,K]);
(%o58) [[I1=-(%i*E0*M-%i*E0*L2)/(%i*L2*R+o*M^2-o*L1*L2),I2=-(%i*E0*R+o*E0*M-o*E0*L1)/(o*L2*R-%i*o^2*M^2+%i*o^2*L1*L2),E=-((%i*E0*M-%i*E0*L2)*R)/(%i*L2*R+o*M^2-o*L1*L2),K=-
(%i*L2*R+o*M^2-o*L1*L2)/((%i*M-%i*L2)*R)]]
(%i59) factor(%);
(%o59) [[I1=-(%i*E0*(M-L2))/(%i*L2*R+o*M^2-o*L1*L2),I2=-(E0*(%i*R+o*M-o*L1))/(o*(L2*R-%i*o*M^2+%i*o*L1*L2)),E=-(%i*E0*(M-L2)*R)/(%i*L2*R+o*M^2-o*L1*L2),K=
(%i*(%i*L2*R+o*M^2-o*L1*L2))/((M-L2)*R)]]

E0とEが同相であるためにはKは実数でなければならない。従ってKの虚数部が０である条件

ω*M^2-ω*L1*L2=0

すなわち

∴M^2=L1*L2

これはL1とL2の間の結合係数が1でなければならないことを示す。

またこの条件でKの実数部は正でなければならないので

L2*R/((L2-M)*R) > 0

従って分母が正の値でなければならず

∴L2 > M

L2がMよりも大きいので

∴L1 < L2

でなければならない。

こうした回路は磁気漏洩の少ないトロイダルコアを用いた高周波回路用伝送線路トランスで良く見られる。

M^2=L1*L2でなければならないので結合係数は1でなければならず、バイファイラ巻きとかトリファイラ巻きとか巻き線間の誘導結合が最も強くなるような独特の巻き方が使われるのはこのためである。

バランやマッチングトランスなどの様々な応用がある。

今度も二つぐらいひねった問題。AB端のインピーダンスを求めるのと、CD端を短絡してもインピーダンスが変わらない条件を導けというもの。



以下の式が成り立つ

(R1+jωL1+R3)*I1+jωM*I2=E

(R2+jωL2+R4)*I2+jωM*I1=E

Z*(I1+I2)=E

これをI1,I2,Zに関する３元連立方程式として解くと

(%i60) e1:(R1+%i*o*L1+R3)*I1+%i*o*M*I2=E;
(%o60) I1*(R3+R1+%i*o*L1)+%i*o*I2*M=E
(%i61) e2:(R2+%i*o*L2+R4)*I2+%i*o*M*I1=E;
(%o61) I2*(R4+R2+%i*o*L2)+%i*o*I1*M=E
(%i62) e3:Z*(I1+I2)=E;
(%o62) (I2+I1)*Z=E
(%i63) solve([e1,e2,e3],[I1,I2,Z]);
(%o63) [[I1=(E*R4+E*R2-%i*o*E*M+%i*o*E*L2)/((R3+R1+%i*o*L1)*R4+(R2+%i*o*L2)*R3+(R1+%i*o*L1)*R2+%i*o*L2*R1+o^2*M^2-o^2*L1*L2),I2=(
(E*R3^2+(2*E*R1-%i*o*E*M+2*%i*o*E*L1)*R3+E*R1^2+(2*%i*o*E*L1-%i*o*E*M)*R1+o^2*E*L1*M-o^2*E*L1^2)*R4^2+(E*
R3^3+(2*E*R2+3*E*R1-3*%i*o*E*M+2*%i*o*E*L2+3*%i*o*E*L1)*R3^2+((4*E*R1-2*%i*o*E*M+4*%i*o*E*L1)*R2+3*E*
R1^2+(-6*%i*o*E*M+4*%i*o*E*L2+6*%i*o*E*L1)*R1-o^2*E*M^2+(2*o^2*E*L2+6*o^2*E*L1)*M-4*o^2*E*L1*L2-3*o^2*E*L1^2)*
R3+(2*E*R1^2+(4*%i*o*E*L1-2*%i*o*E*M)*R1+2*o^2*E*L1*M-2*o^2*E*L1^2)*R2+E*R1^3+
(-3*%i*o*E*M+2*%i*o*E*L2+3*%i*o*E*L1)*R1^2+(-o^2*E*M^2+(2*o^2*E*L2+6*o^2*E*L1)*M-4*o^2*E*L1*L2-3*o^2*E*L1^2)*R1
-%i*o^3*E*M^3-%i*o^3*E*L1*M^2+(2*%i*o^3*E*L1*L2+3*%i*o^3*E*L1^2)*M-2*%i*o^3*E*L1^2*L2-%i*o^3*E*L1^3)*R4+
(E*R2+%i*o*E*L2)*R3^3+(E*R2^2+(3*E*R1-3*%i*o*E*M+2*%i*o*E*L2+3*%i*o*E*L1)*R2+3*%i*o*E*L2*R1+o^2*E*M^2+3*
o^2*E*L2*M-o^2*E*L2^2-3*o^2*E*L1*L2)*R3^2+((2*E*R1-%i*o*E*M+2*%i*o*E*L1)*R2^2+(3*E*R1^2+
(-6*%i*o*E*M+4*%i*o*E*L2+6*%i*o*E*L1)*R1-o^2*E*M^2+(2*o^2*E*L2+6*o^2*E*L1)*M-4*o^2*E*L1*L2-3*o^2*E*L1^2)*R2+
3*%i*o*E*L2*R1^2+(2*o^2*E*M^2+6*o^2*E*L2*M-2*o^2*E*L2^2-6*o^2*E*L1*L2)*R1-3*%i*o^3*E*M^3+
(2*%i*o^3*E*L1-%i*o^3*E*L2)*M^2+(%i*o^3*E*L2^2+6*%i*o^3*E*L1*L2)*M-2*%i*o^3*E*L1*L2^2-3*%i*o^3*E*L1^2*L2)*R3+
(E*R1^2+(2*%i*o*E*L1-%i*o*E*M)*R1+o^2*E*L1*M-o^2*E*L1^2)*R2^2+(E*R1^3+(-3*%i*o*E*M+2*%i*o*E*L2+3*%i*o*E*L1)*
R1^2+(-o^2*E*M^2+(2*o^2*E*L2+6*o^2*E*L1)*M-4*o^2*E*L1*L2-3*o^2*E*L1^2)*R1-%i*o^3*E*M^3-%i*o^3*E*L1*M^2+
(2*%i*o^3*E*L1*L2+3*%i*o^3*E*L1^2)*M-2*%i*o^3*E*L1^2*L2-%i*o^3*E*L1^3)*R2+%i*o*E*L2*R1^3+
(o^2*E*M^2+3*o^2*E*L2*M-o^2*E*L2^2-3*o^2*E*L1*L2)*R1^2+(-3*%i*o^3*E*M^3+(2*%i*o^3*E*L1-%i*o^3*E*L2)*M^2+
(%i*o^3*E*L2^2+6*%i*o^3*E*L1*L2)*M-2*%i*o^3*E*L1*L2^2-3*%i*o^3*E*L1^2*L2)*R1-2*o^4*E*M^4+(o^4*E*L2+3*o^4*E*L1)*M^3+
(o^4*E*L1*L2-o^4*E*L1^2)*M^2+(-o^4*E*L1*L2^2-3*o^4*E*L1^2*L2)*M+o^4*E*L1^2*L2^2+o^4*E*L1^3*L2)/(
(R3^2+(2*R1+2*%i*o*L1)*R3+R1^2+2*%i*o*L1*R1-o^2*L1^2)*R4^3+(R3^3+
(3*R2+3*R1-2*%i*o*M+3*%i*o*L2+3*%i*o*L1)*R3^2+
((6*R1+6*%i*o*L1)*R2+3*R1^2+(-4*%i*o*M+6*%i*o*L2+6*%i*o*L1)*R1+2*o^2*M^2+4*o^2*L1*M-6*o^2*L1*L2-3*o^2*L1^2)*
R3+(3*R1^2+6*%i*o*L1*R1-3*o^2*L1^2)*R2+R1^3+(-2*%i*o*M+3*%i*o*L2+3*%i*o*L1)*R1^2+
(2*o^2*M^2+4*o^2*L1*M-6*o^2*L1*L2-3*o^2*L1^2)*R1+2*%i*o^3*L1*M^2+2*%i*o^3*L1^2*M-3*%i*o^3*L1^2*L2-%i*o^3*L1^3)*R4^2
+((2*R2+2*%i*o*L2)*R3^3+
(3*R2^2+(6*R1-4*%i*o*M+6*%i*o*L2+6*%i*o*L1)*R2+6*%i*o*L2*R1+2*o^2*M^2+4*o^2*L2*M-3*o^2*L2^2-6*o^2*L1*L2)*R3^2
+((6*R1+6*%i*o*L1)*R2^2+
(6*R1^2+(-8*%i*o*M+12*%i*o*L2+12*%i*o*L1)*R1+4*o^2*M^2+8*o^2*L1*M-12*o^2*L1*L2-6*o^2*L1^2)*R2+6*%i*o*L2*R1^2
+(4*o^2*M^2+8*o^2*L2*M-6*o^2*L2^2-12*o^2*L1*L2)*R1-4*%i*o^3*M^3+(4*%i*o^3*L2+4*%i*o^3*L1)*M^2+8*%i*o^3*L1*L2*M-6*
%i*o^3*L1*L2^2-6*%i*o^3*L1^2*L2)*R3+(3*R1^2+6*%i*o*L1*R1-3*o^2*L1^2)*R2^2+(2*R1^3+
(-4*%i*o*M+6*%i*o*L2+6*%i*o*L1)*R1^2+(4*o^2*M^2+8*o^2*L1*M-12*o^2*L1*L2-6*o^2*L1^2)*R1+4*%i*o^3*L1*M^2+4*%i*o^3*
L1^2*M-6*%i*o^3*L1^2*L2-2*%i*o^3*L1^3)*R2+2*%i*o*L2*R1^3+(2*o^2*M^2+4*o^2*L2*M-3*o^2*L2^2-6*o^2*L1*L2)*R1^2+
(-4*%i*o^3*M^3+(4*%i*o^3*L2+4*%i*o^3*L1)*M^2+8*%i*o^3*L1*L2*M-6*%i*o^3*L1*L2^2-6*%i*o^3*L1^2*L2)*R1+o^4*M^4+4*o^4*
L1*M^3+(-4*o^4*L1*L2-2*o^4*L1^2)*M^2-4*o^4*L1^2*L2*M+3*o^4*L1^2*L2^2+2*o^4*L1^3*L2)*R4+
(R2^2+2*%i*o*L2*R2-o^2*L2^2)*R3^3+(R2^3+(3*R1-2*%i*o*M+3*%i*o*L2+3*%i*o*L1)*R2^2+
(6*%i*o*L2*R1+2*o^2*M^2+4*o^2*L2*M-3*o^2*L2^2-6*o^2*L1*L2)*R2-3*o^2*L2^2*R1+2*%i*o^3*L2*M^2+2*%i*o^3*L2^2*M-%i*o^3*
L2^3-3*%i*o^3*L1*L2^2)*R3^2+((2*R1+2*%i*o*L1)*R2^3+
(3*R1^2+(-4*%i*o*M+6*%i*o*L2+6*%i*o*L1)*R1+2*o^2*M^2+4*o^2*L1*M-6*o^2*L1*L2-3*o^2*L1^2)*R2^2+(6*%i*o*L2*R1^2+
(4*o^2*M^2+8*o^2*L2*M-6*o^2*L2^2-12*o^2*L1*L2)*R1-4*%i*o^3*M^3+(4*%i*o^3*L2+4*%i*o^3*L1)*M^2+8*%i*o^3*L1*L2*M-6*
%i*o^3*L1*L2^2-6*%i*o^3*L1^2*L2)*R2-3*o^2*L2^2*R1^2+(4*%i*o^3*L2*M^2+4*%i*o^3*L2^2*M-2*%i*o^3*L2^3-6*%i*o^3*L1*L2^2)*
R1+o^4*M^4+4*o^4*L2*M^3+(-2*o^4*L2^2-4*o^4*L1*L2)*M^2-4*o^4*L1*L2^2*M+2*o^4*L1*L2^3+3*o^4*L1^2*L2^2)*R3+
(R1^2+2*%i*o*L1*R1-o^2*L1^2)*R2^3+(R1^3+(-2*%i*o*M+3*%i*o*L2+3*%i*o*L1)*R1^2+
(2*o^2*M^2+4*o^2*L1*M-6*o^2*L1*L2-3*o^2*L1^2)*R1+2*%i*o^3*L1*M^2+2*%i*o^3*L1^2*M-3*%i*o^3*L1^2*L2-%i*o^3*L1^3)*R2^2
+(2*%i*o*L2*R1^3+(2*o^2*M^2+4*o^2*L2*M-3*o^2*L2^2-6*o^2*L1*L2)*R1^2+
(-4*%i*o^3*M^3+(4*%i*o^3*L2+4*%i*o^3*L1)*M^2+8*%i*o^3*L1*L2*M-6*%i*o^3*L1*L2^2-6*%i*o^3*L1^2*L2)*R1+o^4*M^4+4*o^4*
L1*M^3+(-4*o^4*L1*L2-2*o^4*L1^2)*M^2-4*o^4*L1^2*L2*M+3*o^4*L1^2*L2^2+2*o^4*L1^3*L2)*R2-o^2*L2^2*R1^3+
(2*%i*o^3*L2*M^2+2*%i*o^3*L2^2*M-%i*o^3*L2^3-3*%i*o^3*L1*L2^2)*R1^2+
(o^4*M^4+4*o^4*L2*M^3+(-2*o^4*L2^2-4*o^4*L1*L2)*M^2-4*o^4*L1*L2^2*M+2*o^4*L1*L2^3+3*o^4*L1^2*L2^2)*R1-2*%i*o^5*M^5+
(%i*o^5*L2+%i*o^5*L1)*M^4+4*%i*o^5*L1*L2*M^3+(-2*%i*o^5*L1*L2^2-2*%i*o^5*L1^2*L2)*M^2-2*%i*o^5*L1^2*L2^2*M+%i*o^5*
L1^2*L2^3+%i*o^5*L1^3*L2^2),Z=((%i*R3+%i*R1-o*L1)*R4+(%i*R2-o*L2)*R3+(%i*R1-o*L1)*R2-o*L2*R1+%i*o^2*M^2-%i*o^2*L1*L2)/(%i*R4+%i*R3+%i*R2+%i*R1+2*o*M-o*L2-o*L1)
]]

異常にI2の式が長いのだが大丈夫かMaxima？

Zに関して直交形式に整理すると

(%i68) factor(Z=((%i*R3+%i*R1-o*L1)*R4+(%i*R2-o*L2)*R3+(%i*R1-o*L1)*R2-o*L2*R1+%i*o^2*M^2
-%i*o^2*L1*L2)/(%i*R4+%i*R3+%i*R2+%i*R1+2*o*M-o*L2-o*L1));
(%o68) Z=(%i*R3*R4+%i*R1*R4-o*L1*R4+%i*R2*R3-o*L2*R3+%i*R1*R2-o*L1*R2-o*L2*R1+%i*o^2*M^2-%i*o^2*L1*L2)/(%i*R4+%i*R3+%i*R2+%i*R1+2*o*M-o*L2-o*L1)
(%i69) rectform(%);
(%o69) Z=
((2*o*M-o*L2-o*L1)*(-o*L1*R4-o*L2*R3-o*L1*R2-o*L2*R1)-(-R4-R3-R2-R1)*(R3*R4+R1*R4+R2*R3+R1*R2+o^2*M^2-o^2*L1*L2))/((R4+R3+R2+R1)^2+(2*o*M-o*L2-o*L1)^2)
+(%i*((2*o*M-o*L2-o*L1)*(R3*R4+R1*R4+R2*R3+R1*R2+o^2*M^2-o^2*L1*L2)+(-R4-R3-R2-R1)*
(-o*L1*R4-o*L2*R3-o*L1*R2-o*L2*R1)))/((R4+R3+R2+R1)^2+(2*o*M-o*L2-o*L1)^2)

Zについて直交形式に整理すると

Z=((2*ω*M-ω*L2-ω*L1)*(-ω*L1*R4-ω*L2*R3-ω*L1*R2-ω*L2*R1)-(-R4-R3-R2-R1)*(R3*R4+R1*R4+R2*R3+R1*R2+ω^2*M^2-ω^2*L1*L2))/((R4+R3+R2+R1)^2+(2*ω*M-ω*L2-ω*L1)^2)+(j*((2*ω*M-ω*L2-ω*L1)*(R3*R4+R1*R4+R2*R3+R1*R2+ω^2*M^2-ω^2*L1*L2)+(-R4-R3-R2-R1)*(-ω*L1*R4-ω*L2*R3-ω*L1*R2-ω*L2*R1)))/((R4+R3+R2+R1)^2+(2*ω*M-ω*L2-ω*L1)^2)
=(ω^2*(L1+L2-2*M)*(L1*(R2+R4)+L2*(R1+R3))+(R1+R2+R3+R4)*((R1+R3)*(R2+R4)+ω^2*(M^2-L1*L2))/((R1+R2+R3+R4)^2+ω^2*(L1+L2-2*M))+j*ω*((R1+R2+R3+R4)*(L1*R4+L2*R3+L1*R2+L2*R1)-(L1+L2-2*M)*((R1+R3)*(R2+R4)+ω^2*(M^2-L1*L2)))/(R1+R2+R3+R4)^2+ω^2*(L1+L2-2*M)^2)

ということになる。

同じ式でI1,I2だけ解くと

(%i131) solve([e1,e2],[I1,I2]);
(%o131) [[I1=-(E*(-R4-R2-%i*o*L2)+%i*o*E*M)/(L1*(%i*o*(R4+R2)-o^2*L2)+R3*(R4+R2+%i*o*L2)+R1*(R4+R2+%i*o*L2)+o^2*M^2),I2=
(E*R3+E*R1-%i*o*E*M+%i*o*E*L1)/(L1*(%i*o*(R4+R2)-o^2*L2)+R3*(R4+R2+%i*o*L2)+R1*(R4+R2+%i*o*L2)+o^2*M^2)]]
(%i132) factor(%);
(%o132) [[I1=(E*(R4+R2-%i*o*M+%i*o*L2))/(R3*R4+R1*R4+%i*o*L1*R4+R2*R3+%i*o*L2*R3+R1*R2+%i*o*L1*R2+%i*o*L2*R1+o^2*M^2-o^2*L1*L2),I2=
(E*(R3+R1-%i*o*M+%i*o*L1))/(R3*R4+R1*R4+%i*o*L1*R4+R2*R3+%i*o*L2*R3+R1*R2+%i*o*L1*R2+%i*o*L2*R1+o^2*M^2-o^2*L1*L2)]]

今度は至ってシンプルな式が出る。

ここでAC間の電圧をEac,AD間の電圧をEadとするとCD間の電圧Ecdは

Eac=(R1+jωL1)*I1+jωM*I2

Ead=R2*I2

Ecd=Eac-Ead
=(R1+jωL1)*I1+jωM*I2-R2*I2
=(R1+jωL1)*I1+(jωM-R2)*I2

これに先のI1,I2の式を代入すると

(%i149)
subst((E*(R4+R2-%i*o*M+%i*o*L2))/(R3*R4+R1*R4+%i*o*L1*R4+R2*R3+%i*o*L2*R3+R1*R2+%i*o*L1*R2
+%i*o*L2*R1+o^2*M^2-o^2*L1*L2), I1, Ecd=I2*(%i*o*M-R2)+I1*(R1+%i*o*L1));
(%o149) Ecd=(E*(R1+%i*o*L1)*(R4+R2-%i*o*M+%i*o*L2))/(R3*R4+R1*R4+%i*o*L1*R4+R2*R3+%i*o*L2*R3+R1*R2+%i*o*L1*R2+%i*o*L2*R1+o^2*M^2-o^2*L1*L2)+I2*
(%i*o*M-R2)
(%i150)
subst((E*(R3+R1-%i*o*M+%i*o*L1))/(R3*R4+R1*R4+%i*o*L1*R4+R2*R3+%i*o*L2*R3+R1*R2+%i*o*L1*R2
+%i*o*L2*R1+o^2*M^2-o^2*L1*L2), I2, Ecd=(E*(R1+%i*o*L1)*(R4+R2-%i*o*M+%i*o*L2))/(R3*R4
+R1*R4+%i*o*L1*R4+R2*R3+%i*o*L2*R3+R1*R2+%i*o*L1*R2+%i*o*L2*R1+o^2*M^2-o^2*L1*L2)
+I2*(%i*o*M-R2));
(%o150) Ecd=(E*(R1+%i*o*L1)*(R4+R2-%i*o*M+%i*o*L2))/(R3*R4+R1*R4+%i*o*L1*R4+R2*R3+%i*o*L2*R3+R1*R2+%i*o*L1*R2+%i*o*L2*R1+o^2*M^2-o^2*L1*L2)+
(E*(%i*o*M-R2)*(R3+R1-%i*o*M+%i*o*L1))/(R3*R4+R1*R4+%i*o*L1*R4+R2*R3+%i*o*L2*R3+R1*R2+%i*o*L1*R2+%i*o*L2*R1+o^2*M^2-o^2*L1*L2)
(%i151) factor(%);
(%o151) Ecd=(E*(R1*R4+%i*o*L1*R4-R2*R3+%i*o*M*R3+%i*o*M*R2+%i*o*L2*R1+o^2*M^2-o^2*L1*L2))/(R3*R4+R1*R4+%i*o*L1*R4+R2*R3+%i*o*L2*R3+R1*R2+%i*o*L1*R2+%i*o*L2*R1+o^2*M^2-o^2*L1*L2)

従ってCD間の電圧Ecdが０となるためには分子が０となる条件

R1*R4+j*ω*L1*R4-R2*R3+j*ω*M*R3+j*ω*M*R2+j*ω*L2*R1+ω^2*M^2-ω^2*L1*L2=0

が成り立つ必要がある。このためには左辺の実数部と虚数部が共に０でなければならないので直交形式に直すと

R1*R4-R2*R3+ω^2*M^2-ω^2*L1*L2+j*ω*(L1*R4+M*R3+M*R2+L2*R1)
=R1*R4-R2*R3+ω^2*(M^2-L1*L2)+j*ω*(L1*R4+L2*R1+M*(R2+R3))
=0

従って実数部が

R1*R4-R2*R3+ω^2*(M^2-L1*L2)=0

虚数部が

L1*R4+L2*R1+M*(R2+R3)=0

でなければならない。

実数部の式より

∴R1*R4-R2*R3=ω^2*(L1*L2-M^2)

虚数部の式より

∴M=-(L1*R4+L2*R1)/(R2+R3)

ということになる。

当初式をたてる際にはMは正だという前提だったが、実際には負でないといけないということになる。それと特定のωに限定されるので、CD間の電圧を測定すればその周波数の時に０となることがわかる。この回路はある種のブリッジ回路である。

ちなみに著者のZの式には誤植がある。虚数部の分母の式で+ω^2とあるべきところが-ω^2となってしまっている。

次ぎはちょっと変わった相互誘導回路に流れる電流を求めよというもの。



以下の関係が成り立つ

(jωL1+jωL4)*I1+(jωM1+jωM2)*I2=E

(jωL2+R+jωL3)*I2+(jωM1+jωM2)*I1=0

これをI1,I2に関して解くと

(%i24) e1:(%i*o*L1+%i*o*L4)*I1+(%i*o*M1+%i*o*M2)*I2=E;
(%o24) I2*(%i*o*M2+%i*o*M1)+I1*(%i*o*L4+%i*o*L1)=E
(%i25) e2:(%i*o*L2+R+%i*o*L3)*I2+(%i*o*M1+%i*o*M2)*I1=0;
(%o25) I2*(R+%i*o*L3+%i*o*L2)+I1*(%i*o*M2+%i*o*M1)=0
(%i26) solve([e1,e2],[I1,I2]);
(%o26) [[I1=(E*(R+%i*o*(L3+L2)))/(L4*(%i*o*R+o^2*(-L3-L2))+L1*(%i*o*R+o^2*(-L3-L2))+o^2*M2^2+2*o^2*M1*M2+o^2*M1^2),I2=-
(%i*E*M2+%i*E*M1)/(L4*(%i*R+o*(-L3-L2))+L1*(%i*R+o*(-L3-L2))+o*M2^2+2*o*M1*M2+o*M1^2)]]

整理すると

I1=E*(R+jω*(L3+L2))/(L4*(jωR+ω^2*(-L3-L2))+L1*(jωR+ω^2*(-L3-L2))+ω^2*M2^2+2*ω^2*M1*M2+ω^2*M1^2)
=E*(R+jω*(L2+L3))/(ω^2*((M1+M2)^2-(L1+L4)*(L2+L3))+jωR*(L1+L4))

I2=-(j*E*M2+j*E*M1)/(L4*(jR+ω*(-L3-L2))+L1*(jR+ω*(-L3-L2))+ω*M2^2+2*ω*M1*M2+ω*M1^2)
=-E*j*(M1+M2)/(ω*((M1+M2)^2-(L1+L4)*(L2+L3))+jR*(L1+L4))

実効値はそれぞれの絶対値になるので

|I1|=E*sqrt((R^2+ω^2*(L2+L3)^2)/(ω^4*((M1+M2)^2-(L1+L4)*(L2+L3))^2+ω^2*R^2*(L1+L4)^2))

|I2|=E*sqrt((M1+M2)^2/(ω^2*((M1+M2)^2-(L1+L4)*(L2+L3))^2+R^2*(L1+L4)^2))

ということになる。

著者の提供している回路図には相互インダクタンスがどちらもMと記載されているが、解答を見るとM1,M2と異なっていて答えがまったく合わず焦った。出題問題と正解が違っている典型例。こんなんで採点されたら受験者は不幸としかいいようが無い。

図は著者の解答に合わせてM1,M2と書き直したもの。

次ぎは前問と似て非なる相互誘導回路。L4に流れる電流が０となる電源周波数fを導けというもの。



以下の関係が成り立つ

(jωL1-j/(ωC))*I1+jωM1*I2+(j/(ωC))*I3=E

(jωL2+jωL3)*I2+jωM1*I1+jωM2*I3=0

(jωL4-j/(ωC))*I3+jωM2*I2+(j/(ωC))*I1=0

これをI1,I2,I3に関して解くと

(%i35) e1:(%i*o*L1-%i/(o*C))*I1+%i*o*M1*I2+(%i/(o*C))*I3=E;
(%o35) %i*o*I2*M1+I1*(%i*o*L1-%i/(o*C))+(%i*I3)/(o*C)=E
(%i36) e2:(%i*o*L2+%i*o*L3)*I2+%i*o*M1*I1+%i*o*M2*I3=0;
(%o36) %i*o*I3*M2+%i*o*I1*M1+I2*(%i*o*L3+%i*o*L2)=0
(%i37) e3:(%i*o*L4-%i/(o*C))*I3+%i*o*M2*I2+(%i/(o*C))*I1=0;
(%o37) %i*o*I2*M2+I3*(%i*o*L4-%i/(o*C))+(%i*I1)/(o*C)=0
(%i38) solve([e1,e2,e3],[I1,I2,I3]);
(%o38) [[I1=
(E*(o^2*C*(M2^2-L3*L4-L2*L4)+L3+L2))/(L1*(%i*o^3*C*(M2^2-L3*L4-L2*L4)+%i*o*(L3+L2))+%i*o*(-M2^2+L3*L4+L2*L4)-2*%i*o*M1*M2+(%i*o^3*C*L4-%i*o)*M1^2),I2=
(E*(o^2*C*L4-1)*M1-E*M2)/(L1*(%i*o^3*C*(M2^2-L3*L4-L2*L4)+%i*o*(L3+L2))+%i*o*(-M2^2+L3*L4+L2*L4)-2*%i*o*M1*M2+(%i*o^3*C*L4-%i*o)*M1^2),I3=-
(o^2*C*E*M1*M2+E*(-L3-L2))/(L1*(%i*o^3*C*(M2^2-L3*L4-L2*L4)+%i*o*(L3+L2))+%i*o*(-M2^2+L3*L4+L2*L4)-2*%i*o*M1*M2+(%i*o^3*C*L4-%i*o)*M1^2)]]

I3の式を整理すると

I3=-(ω^2*C*E*M1*M2+E*(-L3-L2))/(L1*(jω^3*C*(M2^2-L3*L4-L2*L4)+j*ω*(L3+L2))+j*ω*(-M2^2+L3*L4+L2*L4)-2*j*ω*M1*M2+(j*ω^3*C*L4-j*ω)*M1^2)
=E*((L2+L3)-ω^2*C*M1*M2)/(j*ω*(L1*(ω^2*C*(M2^2-L4*(L2+L3))+(L2+L3))+(L4*(L2+L3)-M2^2)-2*M1*M2+(ω^2*C*L4-1)*M1^2)

従ってI3が０になるためには

(L2+L3)-ω^2*C*M1*M2=0

でなければならず、上記の式をωについて解くと

ω=sqrt((L2+L3)/(C*M1*M2))

ω=2πfを代入すると

2πf=sqrt((L2+L3)/(C*M1*M2))

fについて解くと

f=sqrt((L2+L3)/(C*M1*M2))/2π

ということになる。

次ぎは超難問。トランスの一次側と二次側が容量結合している回路で、二次側に電流が流れない条件を結合係数kを用いて導けというもの。



現実のトランスは理想的なトランスと違って一次側と二次側の巻き線間に容量結合が存在する。今までの問題ではそれは無いものとして扱ってきたが現実には歴然と存在し、低周波領域では無視できても高周波領域では無視できなくなる。

問題はこの回路図ではCは閉回路を構成していないので回路に影響を与えないように見えてしまう点である。これは電気の不思議である。実際には流れるのであるが、理屈をほじくりだすと電磁気学の世界に入ってしまうので、電気回路理論では等価回路で置き換えて考えると都合が良い。今までの問題では等価回路を持ち出さなくても回路方程式を立てることができたが、この問題ではさすがに等価回路に置き換えないと式が立てられない。



これでやっとCを含む閉回路が現れる。以下の関係が成り立つ

(jω*(L1-M)+jωM)*I1-jωM*I3-jω*(L1-M)*I2=E

(jω*(L1-M)+jω*(L2-M)-j/(ωC))*I2-jω*(L1-M)*I1-jω*(L2-M)*I3=0

(jω*(L2-M)+R+jωM)*I3-jωM*I1-jω*(L2-M)*I2=0

これをI1,I2,I3について解くと

(%i51) e1:(%i*o*(L1-M)+%i*o*M)*I1-%i*o*M*I3-%i*o*(L1-M)*I2=E;
(%o51) I1*(%i*o*M+%i*o*(L1-M))-%i*o*I3*M-%i*o*I2*(L1-M)=E
(%i52) e2:(%i*o*(L1-M)+%i*o*(L2-M)-%i/(o*C))*I2-%i*o*(L1-M)*I1-%i*o*(L2-M)*I3=0;
(%o52) -%i*o*I3*(L2-M)-%i*o*I1*(L1-M)+I2*(%i*o*(L2-M)+%i*o*(L1-M)-%i/(o*C))=0
(%i53) e3:(%i*o*(L2-M)+R+%i*o*M)*I3-%i*o*M*I1-%i*o*(L2-M)*I2=0;
(%o53) I3*(R+%i*o*M+%i*o*(L2-M))-%i*o*I1*M-%i*o*I2*(L2-M)=0
(%i54) solve([e1,e2,e3],[I1,I2,I3]);
(%o54) [[I1=(E*(%i*o^2*C*L2*R-%i*R+o*L2)+E*L1*(%i*o^2*C*R-o^3*C*L2)-2*%i*o^2*C*E*M*R+o^3*C*E*M^2)/(L1*(-o^3*C*L2*R+o*R+%i*o^2*L2)+M^2*(o^3*C*R-%i*o^2)),I2=
(E*L1*(%i*o*C*R-o^2*C*L2)-%i*o*C*E*M*R+o^2*C*E*M^2)/(L1*(-o^2*C*L2*R+R+%i*o*L2)+M^2*(o^2*C*R-%i*o)),I3=(o^2*C*E*M^2+E*M-o^2*C*E*L1*L2)/(L1*(-o^2*C*L2*R+R+%i*o*L2)+M^2*(o^2*C*R-%i*o))]]
(%i55) factor(%);
(%o55) [[I1=-(E*(2*%i*o^2*C*M*R-%i*o^2*C*L2*R-%i*o^2*C*L1*R+%i*R-o^3*C*M^2+o^3*C*L1*L2-o*L2))/(o*(o^2*C*M^2*R-o^2*C*L1*L2*R+L1*R-%i*o*M^2+%i*o*L1*L2)),I2=-
(o*C*E*(%i*M*R-%i*L1*R-o*M^2+o*L1*L2))/(o^2*C*M^2*R-o^2*C*L1*L2*R+L1*R-%i*o*M^2+%i*o*L1*L2),I3=(E*(o^2*C*M^2+M-o^2*C*L1*L2))/(o^2*C*M^2*R-o^2*C*L1*L2*R+L1*R-%i*o*M^2+%i*o*L1*L2)]]

いつもながらMaximaの因数分解は詰めが甘い。

I3について整理すると

I3=(E*(ω^2*C*M^2+M-ω^2*C*L1*L2))/(ω^2*C*M^2*R-ω^2*C*L1*L2*R+L1*R-j*ω*M^2+j*ω*L1*L2)
=(E*(ω^2*C*(M^2-L1*L2)+M))/(ω^2*C*R*(M^2-L1*L2)-j*ω*(M^2-L1*L2)

従ってI3が０になる条件は分子が０となる

ω^2*C*(M^2-L1*L2)+M=0

すなわち

ω^2=-M/(C*(M^2-L1*L2))
=M/(C*(L1*L2-M^2))

ここでMを結合係数の定義

M=k*sqrt(L1*L2)

で置き換えると

ω^2=k*sqrt(L1*L2)/(C*(L1*L2-k^2*L1*L2))
=k*sqrt(L1*L2)/(C*L1*L2*(1-k^2))
=(k/(1-k^2))/(C*sqrt(L1*L2))

ということになる。理想的なC=0やk=1の条件ではω=∞となる。Cが大きくなればなるほどωは低くなる。Cが0でない時はL1,L2が大きければ大きいほどωは低くなる。RやL4に流れる電流とは無関係。

等価回路に置き換えた場合、網目電流法では相互インダクタンスによる電圧降下の極性に注意しないとついつい間違ってしまう。

この問題を等価回路を用いないで解く方法は存在するのだろうか？

P.S

これと良く似た回路がアマチュア無線のSWR計とかに使われているCM形電力計にみられる。Cは容量結合、Mは相互インダクタンスの意味らしい。CM形方向性結合器と呼ばれる。昔は特許で押さえられていたらしい。検索しても詳しく解析したものは見あたらない。なにか秘密があるに違いない。

あと残すところ５問で相互誘導回路は終わり、がんばって突破しよう。

次ぎは相互誘導回路の負荷の電流と電圧が一定の場合、電源と負荷の力率の関係を導けというもの。少々わかりにくい。



今までの問題と逆で、負荷側のインピーダンスとそれに流れる電流が既知で電源電圧が未知数だという点。

負荷の力率をcosφとした場合、以下の関係が成り立つ

(jωL1+|ZL|*(cosφ+jsinφ))*IL+jωM*(I-IL)=E

jωL2*(I-IL)+jωM*IL=E

これらをI,Eに関して解くと

(%i1) e1:(%i*o*L1+abs(ZL)*(cos(p)+%i*sin(p)))*IL+%i*o*M*(I-IL)=E;
(%o1) IL*((%i*sin(p)+cos(p))*abs(ZL)+%i*o*L1)+%i*o*(I-IL)*M=E
(%i2) e2:%i*o*L2*(I-IL)+%i*o*M*IL=E;
(%o2) %i*o*IL*M+%i*o*(I-IL)*L2=E
(%i3) solve([e1,e2],[I,E]);
(%o3) [[I=-((%i*sin(p)+cos(p))*IL*abs(ZL)-2*%i*o*IL*M+%i*o*IL*L2+%i*o*IL*L1)/(%i*o*M-%i*o*L2),E=
(IL*(sin(p)*L2-%i*cos(p)*L2)*abs(ZL)-o*IL*M^2+o*IL*L1*L2)/(%i*M-%i*L2)]]
(%i4) abs(%);
(%o4) [[abs(I)=sqrt((-(sin(p)*IL*abs(ZL))/(o*M-o*L2)+(2*o*IL*M)/(o*M-o*L2)-(o*IL*L2)/(o*M-o*L2)-(o*IL*L1)/(o*M-o*L2))^2+(cos(p)^2*IL^2*ZL^2)/(o*M-o*L2)^2),abs(E)=
sqrt((-(sin(p)*IL*L2*abs(ZL))/(M-L2)+(o*IL*M^2)/(M-L2)-(o*IL*L1*L2)/(M-L2))^2+(cos(p)^2*IL^2*L2^2*ZL^2)/(M-L2)^2)]]
(%i5) factor(%);
(%o5) [[abs(I)=(abs(IL)*sqrt((-4*o*sin(p)*M+2*o*sin(p)*L2+2*o*sin(p)*L1)*abs(ZL)+sin(p)^2*ZL^2+cos(p)^2*ZL^2+4*
o^2*M^2+(-4*o^2*L2-4*o^2*L1)*M+o^2*L2^2+2*o^2*L1*L2+o^2*L1^2))/(abs(o)*abs(M-L2)),abs(E)=(abs(IL)*sqrt((2*o*sin(p)*L1*L2^2-2*o*sin(p)*L2*M^2)*abs(ZL)+sin(p)^2*L2^2*ZL^2+cos(p)^2*L2^2*ZL^2+o^2*M^4-2*o^2*L1*L2*M^2+o^2*L1^2*L2^2))/abs(M-L2)]]

整理すると

|E|=(|IL|*sqrt((2*ω*sinφ*L1*L2^2-2*ω*sinφ*L2*M^2)*|ZL|+sinφ^2*L2^2*ZL^2+cosφ^2*L2^2*ZL^2+ω^2*M^4-2*ω^2*L1*L2*M^2+ω^2*L1^2*L2^2))/|M-L2|
=(|IL|*sqrt(2*ω*L2*sinφ*(L1*L2-M^2)*|ZL|+L2^2*ZL^2*(sinφ^2+cosφ^2)+ω^2*(M^4-2*L1*L2*M^2+L1^2*L2^2)))/|M-L2|
=(|IL|*sqrt(2*ω*L2*sinφ*(L1*L2-M^2)*|ZL|+L2^2*ZL^2+ω^2*(M^2-L1*L2)^2))/|M-L2|
=(|IL|*sqrt(L2^2*ZL^2+ω*(L1*L2-M^2)*(2*L2*sinφ*|ZL|+ω*(L1*L2-M^2))/|M-L2|

ということになる。

次ぎは一次側と二次側共に共振回路となっている問題。



一次側と二次側のそれぞれの共振周波数をω1,ω2とした場合の回路全体の共振周波数ω0を結合係数を用いて導けというもの。

未知数は共振周波数と一次側と二次側に流れる電流である。

以下の関係式が成り立つ

(jωL1-j/(ωC1))*I1+jωM*I2=E

(jωL2-j/(ωC2))*I2+jωM*I1=0

(%i6) e1:(%i*o*L1-%i/(o*C1))*I1+%i*o*M*I2=E;
(%o6) %i*o*I2*M+I1*(%i*o*L1-%i/(o*C1))=E
(%i7) e2:(%i*o*L2-%i/(o*C2))*I2+%i*o*M*I1=0;
(%o7) %i*o*I1*M+I2*(%i*o*L2-%i/(o*C2))=0
(%i8) solve([e1,e2],[I1,I2]);
(%o8) [[I1=(C1*E*(%i*o^3*C2*L2-%i*o))/(o^4*C1*C2*M^2+C1*L1*(o^2-o^4*C2*L2)+o^2*C2*L2-1),I2=-(%i*o^3*C1*C2*E*M)/(o^4*C1*C2*M^2+C1*L1*(o^2-o^4*C2*L2)+o^2*C2*L2-1)]]

回路全体が共振状態の場合には一次側の電流が∞になる時なので、I1の分母が０となる条件

ω^4*C1*C2*M^2+C1*L1*(ω^2-ω^4*C2*L2)+ω^2*C2*L2-1=0

ここで

M=k*sqrt(L1*L2)

ω1=1/sqrt(C1*L1)

∴C1*L1=1/ω1^2

ω2=1/sqrt(C2*L2)

∴C2*L2=1/ω2^2

で書き換えると

ω^4*C1*C2*M^2+C1*L1*(ω^2-ω^4*C2*L2)+ω^2*C2*L2-1
=ω^4*C1*C2*k^2*L1*L2+ω^2*(1-ω^2/ω2^2)/ω1^2+ω^2/ω2^2-1
=ω^4*k^2/(ω1^2*ω2^2)+ω^2*(1-ω^2/ω2^2)/ω1^2+ω^2/ω2^2-1
=ω^2*(ω^2*k^2/(ω1^2*ω2^2)+(1/ω1^2-ω^2/(ω1^2*ω2^2))+1/ω2^2)-1
=ω^2*(ω^2*(k^2-1)/(ω1^2*ω2^2)+1/ω1^2+1/ω2^2)-1
=0

両辺にω1^2*ω2^2を乗じると

ω^2*(ω^2*(k^2-1)+ω1^2+ω^2)-ω1^2*ω2^2
=ω^4*(k^2-1)+ω^2*(ω1^2+ω^2)-ω1^2*ω2^2
=0

これをωについて解くとω0は

(%i52) solve([o^2*(o^2*(k^2-1)/(o1^2*o2^2)+1/o1^2+1/o2^2)-1], [o]);
(%o52) [o=-sqrt(sqrt(o2^4+4*k^2*o1^2*o2^2-2*o1^2*o2^2+o1^4)/(k^2-1)-o2^2/(k^2-1)-o1^2/(k^2-1))/sqrt(2),o=
sqrt(sqrt(o2^4+4*k^2*o1^2*o2^2-2*o1^2*o2^2+o1^4)/(k^2-1)-o2^2/(k^2-1)-o1^2/(k^2-1))/sqrt(2),o=-sqrt(sqrt(o2^4+4*k^2*o1^2*o2^2-2*o1^2*o2^2+o1^4)+o2^2+o1^2)/(sqrt(2)*sqrt(1-k^2)),o=
sqrt(sqrt(o2^4+4*k^2*o1^2*o2^2-2*o1^2*o2^2+o1^4)+o2^2+o1^2)/(sqrt(2)*sqrt(1-k^2))]
(%i53) factor(%);
(%o53) [o=-sqrt((sqrt(o2^4+(4*k^2-2)*o1^2*o2^2+o1^4)-o2^2-o1^2)/(k^2-1))/sqrt(2),o=sqrt((sqrt(o2^4+(4*k^2-2)*o1^2*o2^2+o1^4)-o2^2-o1^2)/(k^2-1))/sqrt(2),o=-
sqrt(sqrt(o2^4+(4*k^2-2)*o1^2*o2^2+o1^4)+o2^2+o1^2)/(sqrt(2)*sqrt(1-k^2)),o=sqrt(sqrt(o2^4+(4*k^2-2)*o1^2*o2^2+o1^4)+o2^2+o1^2)/(sqrt(2)*sqrt(1-k^2))]
(%i54) solve([o^4*(k^2-1)+o^2*(o1^2+o2^2)-o1^2*o2^2], [o]);
(%o54) [o=-sqrt(sqrt(o2^4+4*k^2*o1^2*o2^2-2*o1^2*o2^2+o1^4)/(k^2-1)-o2^2/(k^2-1)-o1^2/(k^2-1))/sqrt(2),o=
sqrt(sqrt(o2^4+4*k^2*o1^2*o2^2-2*o1^2*o2^2+o1^4)/(k^2-1)-o2^2/(k^2-1)-o1^2/(k^2-1))/sqrt(2),o=-sqrt(sqrt(o2^4+4*k^2*o1^2*o2^2-2*o1^2*o2^2+o1^4)+o2^2+o1^2)/(sqrt(2)*sqrt(1-k^2)),o=
sqrt(sqrt(o2^4+4*k^2*o1^2*o2^2-2*o1^2*o2^2+o1^4)+o2^2+o1^2)/(sqrt(2)*sqrt(1-k^2))]
(%i55) factor(%);
(%o55) [o=-sqrt((sqrt(o2^4+(4*k^2-2)*o1^2*o2^2+o1^4)-o2^2-o1^2)/(k^2-1))/sqrt(2),o=sqrt((sqrt(o2^4+(4*k^2-2)*o1^2*o2^2+o1^4)-o2^2-o1^2)/(k^2-1))/sqrt(2),o=-
sqrt(sqrt(o2^4+(4*k^2-2)*o1^2*o2^2+o1^4)+o2^2+o1^2)/(sqrt(2)*sqrt(1-k^2)),o=sqrt(sqrt(o2^4+(4*k^2-2)*o1^2*o2^2+o1^4)+o2^2+o1^2)/(sqrt(2)*sqrt(1-k^2))]

ω0は正の実数でなければならないので

ω0=sqrt((sqrt(ω2^4+(4*k^2-2)*ω1^2*ω2^2+ω1^4)-ω2^2-ω1^2)/(k^2-1))/sqrt(2)
=sqrt((sqrt(4*k^2*ω1^2*ω2^2+(ω1^2-ω2^2)^2)-(ω1^2+ω2^2))/(2*(k^2-1)))
=sqrt(((ω1^2+ω2^2)-sqrt(4*k^2*ω1^2*ω2^2+(ω1^2-ω2^2)^2))/(2*(1-k^2)))

ω0'=sqrt(sqrt(ω2^4+(4*k^2-2)*ω1^2*ω2^2+ω1^4)+ω2^2+ω1^2)/(sqrt(2)*sqrt(1-k^2))
=sqrt(sqrt(4*k^2*ω1^2*ω2^2+(ω1^2-ω2^2)^2)+(ω1^2+ω2^2))/(sqrt(2)*sqrt(1-k^2))
=sqrt(((ω1^2+ω2^2)+sqrt(4*k^2*ω1^2*ω2^2+(ω1^2-ω2^2)^2))/(2*(1-k^2)))

ということになる。

P.S

こうした共振回路を誘導結合したトランスは受信機の中間周波トランスとして応用されている。式の上では２つのピークを持つが、結合係数を調整すると２こぶになったり台形になったり、山が重なったりする。

Maximaを使ってある回路定数で結合係数を変えてアドミッタンスの周波数特性をプロットしてみると太平洋海底のハワイ列島のような軌跡が見える。

C1=10pF,C2=20pF,L1=L2=1uH

plot3d((500000*%pi*abs(f)*abs(%pi^2*f^2-12500000000000000))/abs(%pi^4*f^4*k^2-%pi^4*f^4
+37500000000000000*%pi^2*f^2-312500000000000000000000000000000), [k,0,1], [f,10^3,5*10^8],[grid,50,50])$



もともと２つの共振回路の共振点は隣接しているが、ある結合係数では山が完全に２つにスプリットしているのがわかる。

演習問題では一次側は直列共振だったが、中間周波トランスと同じように並列共振にしたらどうなるかとか一次側と二次側の出力比(ゲイン)がどうなるか研究してみると面白いかもしれない。

次ぎは一次側と二次側にそれぞれ負荷が直列に接続された回路。回路全体のインピーダンスを求め、それぞれの負荷が抵抗負荷となった場合に、全体で消費される電力は負荷で消費された電力の合計となることを示せというもの。



著者はここにきて等価回路を使わずに回路方程式で解いている。

別解法として等価回路を使って解いてみよう。



全体のインピーダンスは直並列混成回路なので

Z=Z1+jω*(L1-M)+1/(1/(jωM)+1/(jω*(L2-M)+Z2))
=Z1+jω*(L1-M)+1/((jωM+jω*(L2-M)+Z2)/((jωM)*(jω*(L2-M)+Z2)))
=Z1+jω*(L1-M)+(jωM)*(Z2+jωL2-jωM)/(Z2+jωL2)
=Z1+jω*(L1-M)+jωM+(jωM)*(-jωM)/(Z2+jωL2)
=Z1+jωL1+ω^2*M^2/(Z2+jωL2)

ということになる。

続いてZ1,Z2をそれぞれR1,R2とした場合、インピーダンスの式は

Z=R1+jωL1+ω^2*M^2/(R2+jωL2)
=R1+jωL1+ω^2*M^2*(R2-jωL2)/(R2^2+ω^2*L2^2)
=R1+ω^2*M^2*R2/(R2^2+ω^2*L2^2)+jω*(L1-ω^2*M^2*L2/(R2^2+ω^2*L2^2))

ということになる。従って全体で消費される電力は実効抵抗で消費される電力であるため

Pa=(R1+ω^2*M^2*R2/(R2^2+ω^2*L2^2))*|I1|^2

ということになる。

I1はR1に流れる電流である。R2に流れる電流は分流の法則によって簡単に求めることが出来る。

I2=(jωM/(R2+jω*(L2-M)+jωM))*I1
=(jωM/(R2+jωL2))*I1

従って

|I2|=(ω^2*M^2/(R2^2+ω^2*L2^2))*|I1|

∴Pa=R1*|I1|^2+R2*|I2|^2

である。

確かに方程式をたてるより等価回路を使ったほうが手計算で出来る。

次ぎの問題は一次側と二次側が共に共振回路を構成し、共振周波数がω0で同一の場合、結合係数を変化させて二次側に流れる電流が最も大きくなる結合係数とその時の電流の式を導けというもの。



これも等価回路に置き換えて解いてみよう。



全体のインピーダンスは

Z=R1-j/(ωC1)+jω*(L1-M)+1/(1/(jωM)+1/(jω*(L2-M)+R2-j/(ωC2)))
=R1+j*(ω*(L1-M)-1/(ωC1))+jωM*(R2+j*(ω*(L2-M)-1/(ωC2)))/(R2+j*(ωL2-1/(ωC2)))
=R1+j*(ω*(L1-M)-1/(ωC1))+jωM*(1-jωM/(R2+j*(ωL2-1/(ωC2))))
=R1+j*(ωL1-1/(ωC1))+ω^2*M^2/(R2+j*(ωL2-1/(ωC2)))
=R1+j*(ωL1-1/(ωC1))+ω^2*M^2*(R2-j*(ωL2-1/(ωC2)))/(R2^2+(ωL2-1/(ωC2))^2)
=R1+ω^2*M^2*R2/(R2^2+(ωL2-1/(ωC2))^2)+j*(ωL1-1/(ωC1)-ω^2*M^2*(ωL2-1/(ωC2))/(R2^2+(ωL2-1/(ωC2))^2))

ここで共振角周波数ω0では

ω0=1/sqrt(C1*L1)

∴ω0*L1-1/(ω0*C1)=0

ω0=1/sqrt(C2*L2)

∴ω0*L2-1/(ω0*C2)=0

を適用すると

Z0=R1+ω0^2*M^2/R2



従って一次側を流れる電流I1は

I1=E/Z0
=E/(R1+ω0^2*M^2/R2)
=E*R2/(R1*R2+ω0^2*M^2)

実効値は

|I1|=E*R2/(R1*R2+ω0^2*M^2)

また二次側を流れる電流I2は分流則により

I2=(jωM/(jω(L2-M)+R2-j/(ωC2)+jωM))*I1
=(jωM/(R2+j*(ωL2-1/(ωC2))))*I1

ω0の時はω0*L2-1/(ω0*C2)=0なので

I2=(jω0*M/R2)*I1

従って実効値は

|I2|=(ω0*M/R2)*|I1|

|I1|の式を代入すると

|I2|=(ω0*M/R2)*E*R2/(R1*R2+ω0^2*M^2)
=E*ω0*M/(R1*R2+ω0^2*M^2)

ということになる。

ここで誘導係数の定義

M=k*sqrt(L1*L2)

を代入すると

|I2|=E*ω0*k*sqrt(L1*L2)/(R1*R2+ω0^2*k^2*L1*L2)

|I2|が最大になる誘導係数を求めるにはkで微分すると

(%i69) diff(o0*k*sqrt(L1*L2)/(R1*R2+o0^2*k^2*L1*L2),k);
(%o69) (o0*sqrt(L1*L2))/(R1*R2+k^2*o0^2*L1*L2)-(2*k^2*o0^3*L1*L2*sqrt(L1*L2))/(R1*R2+k^2*o0^2*L1*L2)^2
(%i70) factor(%);
(%o70) (o0*sqrt(L1*L2)*(R1*R2-k^2*o0^2*L1*L2))/(R1*R2+k^2*o0^2*L1*L2)^2

従って微分係数が０となる条件は

R1*R2-k^2*ω0^2*L1*L2=0

これをkについて解くと

(%i71) solve([R1*R2-k^2*o0^2*L1*L2], [k]);
(%o71) [k=-sqrt((R1*R2)/(L1*L2))/o0,k=sqrt((R1*R2)/(L1*L2))/o0]

kは正の値でなければならないので|I2|が最大値を取るkは

kmax=sqrt(R1*R2/L1*L2)/ω0

これを|I2|の式に代入すると

(%i80) subst(sqrt(R1*R2/(L1*L2))/o0, k, (k*o0*E*sqrt(L1*L2))/(R1*R2+k^2*o0^2*L1*L2));
(%o80) (E*sqrt(L1*L2)*sqrt((R1*R2)/(L1*L2)))/(2*R1*R2)

またしてもMaximaの因数分解は中途半端だ

|I2max|=E*sqrt(L1*L2)*sqrt((R1*R2)/(L1*L2))/(2*R1*R2)
=E*sqrt(R1*R2)/(2*R1*R2)
=E/(2*sqrt(R1*R2))

ということになる。結果的にL1,L2は消え、L1,L2,C1,C2に無関係にR1,R2のみによって定まる。すなわち最大電流値がコイルの導体抵抗値のみによって決まるということになる。

実際には共振点は誘導結合が無い時に一緒でも、誘導結合し始めると前前問の解の式よりR1,R2が無視できる程小さい場合には

ω0と

ω0'=ω0*sqrt((1+k)/(1-k))

の２つの山に分かれることになる。



上のプロットではR1,R2を1としているので結合係数が大きくなると片方の山が少し低い周波数へゆっくり移動していく。

いよいよ相互誘導回路の問題の最後。２つの交流電源が供給される回路で消費される全電力の瞬時値と回路の磁気エネルギーを導けというもの。



著者の解を見るといままでの問題と題意は違ったところにあるように見える。単純に考えれば回路全体の瞬時値電力は電源E1,E2とそれから供給される電流I1,I2のそれぞれの積の合計である。

p=E1*I1+E2*I2

またE1,E2に関して以下の関係式が成り立つ

(R1+jωL1)*I1+jωM1*I=E1

(R2+jωL2)*I2+jωM2*I=E2

従ってこれらを先の瞬時値電力の式に代入すると

p=E1*I1+E2*I2
=(R1+jωL1)*I1^2+jωM1*I*I1+(R2+jωL2)*I2^2+jωM2*I*I2

ここで

(jωL3+jωL4)*I+jωM1*I1+jωM2*I2=0

であることからIについて解くと

I=-(jωM1*I1+jωM2*I2)/(jωL3+jωL4)
=-(M1*I1+M2*I2)/(L3+L4)

これを代入すると

p=(R1+jωL1)*I1^2-jωM1*(M1*I1+M2*I2)*I1/(L3+L4)+(R2+jωL2)*I2^2-jωM2*(M1*I1+M2*I2)*I2/(L3+L4)
=(R1+jωL1)*I1^2+(R2+jωL2)*I2^2-jω*(M1^2*I1^2+M1*M2*I1*I2+M1*M2*I1*I2+M2^2*I2^2)/(L3+L4))
=R1*I1^2+R2*I2^2+jω*((L1-M1^2/(L3+L4))*I1^2+(L2-M2^2/(L3+L4))*I2^2-2*M1*M2*I1*I2/(L3+L4))

ということになる。

磁気エネルギーは回路の実効インダクタンスと電源の間で交換される瞬時値電力から抵抗で消費される有効電力を差し引いた無効電力を2*jωで除したものとなる。

wL=(L1-M1^2/(L3+L4)*I1^2/2+(L2-M2^2/(L3+L4))*I2^2/2-M1*M2*I1*I2/(L3+L4)

著者の解は最初の瞬時値電力の式を求める際に致命的なミスを犯している。(M1*i1+M2*i2)*di/dtを展開する際にdi/dtの符号が負なので最終的にM1*M2の項は負とならなければならないが正になってしまっている。

またM1*M2の項が負になるのは、結合係数が1となって漏洩インダクタンスが無くなった場合に磁気エネルギーが相殺されなければならない観点からして理にかなっている。M1*M2の項が正だと結合係数が1になっても実効インダクタンスが決して０にならない式となるため矛盾する。