表皮効果とベッセル関数

表皮効果とベッセル関数

はじめに

一般的に送電においては円筒状の銅線が用いられるが、高周波電流を流す場合は、渦電流により表面に電流が流れやすくなる表皮効果の影響により、銅線の抵抗値が増加してしまう。今回は、天下り的ではあるが、ベッセル関数を用いて表皮効果の影響についてPythonを用いてグラフ化することを試みる。

表皮効果の定量化

まず、モデルとしては、円筒座標を用いて解析するわけだが、ここでは以下の記事の解のみを引用する。

電流密度Jz(r)はベッセル関数を使い厳密に解くと、以下の様に表せる。

ただし、σ:電気伝導率,μ:透磁率,ω:角速度,I:総電流,r:半径,r0:総半径

これをpythonプログラムとして書くと以下の様になる。

Python


import matplotlib.pyplot as plt
import math 
import numpy as np
#ベッセル関数の読み込み
from scipy.special import jv

#円筒銅線の材料パラメーター

#電気伝導率σ
sigma = 1.0
#透磁率μ
mu = 1.0

#角周波数
omega = 10
#電流値
I=1.0
#銅線の半径
r_0=100
k= (-1j*omega*mu*sigma)**0.5
r = np.linspace(-1*r_0,r_0,100)
#電流密度(rの関数)
y = (k*I/(2*math.pi))*jv(0,k*r)/jv(1,k*r_0)

plt.plot(r,y)
plt.show()

結果のグラフ表示

ωを1から100まで動かしたときの結果をグラフ表示すると以下の様になる。

ω=100のとき

表面付近に電流が集中していることが分かる。

ω=10のとき

ω=1のとき

表面付近に電流は集中はしているが影響は小さい。

まとめ

周波数が大きい場合は、電流が銅線の表面付近に流れやすくなる表皮効果の影響が無視できなくなってしまうということが、ベッセル関数を用いた解析モデルによって分かった。ゆえに、流れる電流路の断面積が小さくなるので抵抗値が増加してしまうので、交流機器を作る時はそのような影響も考慮する必要性がある。

参考:円筒表示

計算量は格段に大きくなるのだが、円筒座標にして対称性を考慮せずにコンター図を作成するプログラムを書くと以下の通りである。

import matplotlib.pyplot as plt
import math 
import numpy as np
#ベッセル関数の読み込み
from scipy.special import jv
#3dデータを上手く加工する
from scipy.interpolate import griddata



#円筒銅線の材料パラメーター

#電気伝導率σ
sigma = 1.0
#透磁率μ
mu = 1.0

#角周波数
omega = 10
#電流値
I=1.0
#銅線の半径
r_0=10

k= (-1j*omega*mu*sigma)**0.5

#プロット

NUM=10000                          # サンプリング数
x=0                                     # 点のx値
y=0                                     # 点のy値                                  
#zを求めるために使用
i=0   



#a<r<b つまり -r_0<r<r_0の範囲で乱数を生成する
a=-1*r_0
b=r_0

x_a=[]
y_a=[]
z_a=[]
                                  # ループカウンタ
while i<NUM:                            # NUM回ループ
    x=(b - a) * np.random.rand() + a                  # ランダムな点の位置x a<x<bにおけるランダムな点を作成
    y=(b - a) * np.random.rand() + a                  # ランダムな点の位置y
    if x**2+y**2<r_0**2:                      # 円筒内部ならば
 
        z= (k*I/(2*math.pi))*jv(0,k*(x**2+y**2)**0.5)/jv(1,k*r_0)
        x_a.append(x)
        y_a.append(y)
        z_a.append(z)
    else:                               # 円の外
        #plt.scatter(x,y,color='blue')
        z=0   
    i+=1                                # ループカウンタ増


#データの整理
x_new, y_new = np.meshgrid(np.unique(x_a), np.unique(y_a))
z_new = griddata((x_a, y_a), z_a, (x_new, y_new))


#コンター図で表示する
plt.contourf(x_new, y_new, z_new)
#カラーバー
plt.colorbar()
plt.show()

これにより作成したコンター図は以下の様になる。

計算量は大きくなるが、上手く表皮効果のイメージをつかむことが出来た。

一方で、コンター図でなく散布図で表すと以下の様なグラフになる。

こちらも、銅線表面付近に行くほど電流が集中していることが分かる。

上記は、ω=10の場合なので、ω=1,100の場合をそれぞれ示して締めようと思う。

(1)ω=100の場合について

この様に、ω=10の場合と比較して表面により電流が集中していることが分かる。なので、高周波を扱う際には表面付近での加熱による熱損失や熱破壊による誤動作などが問題となってくる。

(2)ω=1の場合について

電流は、上2つの場合と比較してなだらかで中央付近にもより多く分布しているということが分かった。



電材館 第二種電気工事士 技能試験セット 「 電線2回分 」 「 ホーザン 工具 」 「 動画解説・テキスト付き 」 2023年


参加者500人以上
電験三種や電気工事士など
資格合格を目指す人が集まる
無料のオンラインサロン
【みんなのデンキ塾】

参加申請は公式ツイッターから



関連記事

コメント

  1. この記事へのコメントはありません。

コメントするためには、 ログイン してください。