かご型誘導機の二次元解析におけるロータバーとエンドリングの取り扱い

かご型誘導電動機（かご型誘導モータ）の磁場解析を行うためには、回転子内のロータバーとエンドリング（概念図Fig.1）を有限要素に分割し、三次元解析する方法が考えられます。しかし、三次元解析では多大の計算時間を必要とし、さらに一般に誘導モータの解析は定常に達するのが遅いことから、端部の影響を考慮しない二次元解析が多く行われております。EMSolutionでもその解析例を「誘導電動機の解析」に示しています。エンドリングを解析に含めず、その効果（主とて抵抗によるもの）を、ロータバーの導電率を小さくすることで近似的に扱っております。この場合、ロータバーの軸方向電流の面積分値はゼロとなるはずですので、SUFCUR（面流入電流ソース）を用いて、その軸方向電流面積分値のトータルをゼロになるように設定しております。二次元解析の場合、このような拘束条件を入れなければ、一般的に軸方向電流積分値がゼロになりません。

二次元解析で、より現実的な近似法として、エンドリングを抵抗で模擬して電気回路に組み込み、有限要素法と連成して解析する方法が考えられます。Fig.2にその等価回路を示します。有限要素解析領域の点線で示した四角形内はロータバーを表します。また、抵抗はエンドリングを表し、上下で隣接するロータバーを連結しています。ロータバーが等間隔に配置されていると仮定しますと、抵抗値は全て等しくなります。さらに、ロータバーの端部の抵抗を加えたり、インダクティブな効果を入れたりするためのインダクタンスを挿入することも考えられますが、ここでは簡単のため抵抗値のみを扱います。また、上部と下部は同じ電流が流れますので、下部は節点（Fig.2内の黒丸）間を抵抗ゼロで連結し、上部の抵抗値を倍にしても構いません。
Fig.2のような電気回路は、NETWORKモジュールを使用すれば容易に作成でき、有限要素法との連成解析が可能となります。回転周期対称性を考慮しない360度モデルの場合、Fig.2のようにすべてのロータバーとエンドリングを考慮すればよいのですが、多くの場合対称性がありますので、その場合はFig.3のように等価的に分割して有限要素領域と合わせ、その周期分だけ（網目部）をモデル化すれば解析することができます。このとき、エンドリングに流れる電流は、REGION_FACTOR倍とはなりませんので、エンドリングの抵抗をREGION_FACTOR倍として設定し、全周モデルで流れる電流値と等しくします。REGION_FACTORを1としますと、抵抗値は全周モデルの場合と変わらず設定できますが、入力電圧、抵抗等は1/2と設定することになります。なお、ロータバーの導電率はREGION_FACTORによらず物性データそのものの値を設定します。

「誘導電動機の解析」のモデルにエンドリングを抵抗として近似した例を示します。Fig.4に有限要素メッシュを示します。回転子のロータバーは、全て異なるプロパティ番号を設定する必要があります。また、周期対称面のロータバーは回転対称となりますので、同じプロパティ番号を設定しています。「誘導電動機の解析」では電源回路はCIRCUITで入力し、前述したようにロータバーの導電率をエンドリングの効果を考慮して設定していますが、ここで設定するNETWORKでは、エンドリングの抵抗値を別々に設定できます。
三相交流電源を設定するには、三相独立に回路を設定することもできますが、ここではY字結線することを考えます。CIRCUITで設定する場合、CONNECTION_MATRIXを入力する必要があり、直感的に設定することが容易ではありません。NETWORKを使用すると、実際の電源回路そのままを設定することができます。ここでは、相間電圧が100/√3V（実効値）となるようにY字結線しております。電流電源入力とする場合、中性点を大きな抵抗で結ぶ必要があります。「誘導電動機の解析」ではΔ結線をしており、相間電圧が100V（実効値）となります。この場合の結線方法と入力例、および三相交流電源の結線の注意点は、「Δ結線とY結線」をご参照下さい。

現実的には、エンドリングの抵抗は表皮効果による周波数依存性もあり、またエンドリングの形状にも依存すると考えられますので、一概に決定することは難しいと思われます。むしろ、実験結果に良く合うように、抵抗値を調整するのが現実的かと思われます。

＜交流定常解析の入力例＞

使用データ：inputNETWORK_Y_AC：交流定常解析
　　　　　　　inputNETWORK_Y_1123rpm：過渡解析
　　　　　　　pre_geom2D.neu：固定子メッシュデータ
　　　　　　　rotor_mesh2D.neu：回転子メッシュデータ