
はじめに
RLC直列交流回路の電流値を求める際は、単に電流の大きさだけでなく、位相が電圧に関して進むのか遅れるのかまでチェックしなければならない。一般的に位相差を求めるためには、三角関数の逆関数や三角関数の合成公式を使用する必要性があるため、その部分をPythonによる数値解析に置き換えて、電流と電圧の特性を調べることを本記事の目的とした。
また、最後に三平方の定理を用いた電力計算についても言及した。
電流値の導出
回路の抵抗、インダクタンス、キャパシタンスをそれぞれ以下の様にする。
この場合、RLC直列回路に流れる電流は、以下の様に表すことが出来る。
ただし、電圧は以下の様に表されるとする。
また、αは以下の様な三角形の底角とする。
つまり、
Pythonプログラム
このことに留意してPythonでプログラミングを書くと以下の様になる。
import numpy as np import matplotlib.pyplot as plt import math R=int((input("抵抗値を入力してください"))) X_L=int((input("インダクタンスを入力してください"))) X_C=int((input("キャパシタンスを入力してください"))) V=100 a=float((X_L-X_C)/R) x = np.arange(-5, 5, 0.1) y1 = V*np.sin(x) y2 = V*np.sin(x+math.atan(a))/(R**2+(X_L-X_C)**2)**0.5 plt.title("交流:直列RCL回路における電流特性", fontname="MS Gothic")#日本語化 plt.xlabel("time") plt.ylabel("size") plt.plot(x, y1, label="V",) # labelに凡例に表示する文字を指定 plt.plot(x, y2, label="I") plt.legend() # 凡例を表示 plt.xlim(-np.pi, np.pi) plt.ylim(-1.2*V, 1.2*V) plt.show()
ここでは、R=1Ω,XL=2Ω,XC=1Ωの時の遅れ位相のグラフを以下に示す。
このように、電流の位相が基準となる電圧の位相に比べて遅れてしまうことが分かる。
次に、R=2Ω,XL=2Ω,XC=2Ωの時のグラフを以下に示す。
この様に、電流と電圧が同相になることが分かる。
最後に、電流が進みになる場合のグラフを以下に示す。
(R=1Ω,XL=2Ω,XC=3Ω)
電力計算
インピーダンス計算、電圧計算、電力計算まで考慮に入れたプログラムを以下に示す。
import numpy as np import matplotlib.pyplot as plt import math R=int((input("抵抗値を入力してください"))) X_L=int((input("インダクタンスを入力してください"))) X_C=int((input("キャパシタンスを入力してください"))) V=100 a=float((X_L-X_C)/R) x = np.arange(-5, 5, 0.1) y1 = 2**0.5*V*np.sin(x) y2 = 2**0.5*V*np.sin(x+math.atan(a))/(R**2+(X_L-X_C)**2)**0.5 plt.title("交流:直列RCL回路における電流特性", fontname="MS Gothic")#日本語化 plt.xlabel("time") plt.ylabel("size") plt.plot(x, y1, label="V",) # labelに凡例に表示する文字を指定 plt.plot(x, y2, label="I") plt.legend() # 凡例を表示 plt.xlim(-np.pi, np.pi) plt.ylim(-1.2*V, 1.2*V) plt.show() #以下力率と電力計算 Z=(R**2+(X_L-X_C)**2)**0.5 #インピーダンス I=V/Z V_R=R*I V_C=X_C*I V_L=X_L*I print("電圧の三角形が成立すれば0になるはず") print(V**2-(V_R**2+(V_L-V_C)**2)) P=V_R*I Q=(V_L-V_C)*I S=V*I print("電力の三角形が成立すれば0になるはず") print(S**2-Q**2-P**2) print("力率") theta = P/S print(P/S)
この様に、RLC直列回路の交流計算は、三平方の定理を用いると上手く計算することができるということがわかる。
まとめ
今回はPythonを用いてRLC直列回路の位相解析を行った。これにより、リアクトルがキャパシタンスよりも優位となる一般的な状態の場合は、電流が電圧よりも遅れの位相を取ることが分かった。一方で、インダクタンスとキャパシタンスが同じ場合は、位相が電圧と電流で同相になり、キャパシタンス優位の場合は電流の位相が電圧の位相に対して進みになることが確認できた。また、電力計算では、電力の三角形を書くことが重要であることが分かった。

電材館 第二種電気工事士 技能試験セット 「 電線2回分 」 「 ホーザン 工具 」 「 動画解説・テキスト付き 」 2023年
参加者500人以上
電験三種や電気工事士など
資格合格を目指す人が集まる
無料のオンラインサロン
【みんなのデンキ塾】
コメント