簡易版時間周期法による誘導電動機の定常解析　　　

　かご型誘導電動機の定常状態の近似解は交流定常解析によっても求めることが可能であり，それを初期値として時間ステップ解析を行えば素早く定常解を求めることができることを誘導電動機の解析において示しました。交流定常解析は線形計算を行うため，適当な透磁率を与えて計算する必要があります。 　しかし，非線形性が強いとき，例えばソリッドロータの誘導電動機のような場合など，線形定常解析結果は定常解から大きく異なってしまう場合が多いと思われます。このような場合，交流定常解を初期値として時間ステップ解析を行っても定常解に達するには多大のステップを要してしまいます。 　 時間周期問題の定常解への高速収束で紹介したように，同期回転機の場合，SD-EEC時間周期有限要素法は定常解を高速に得ることができます。しかし，この方法は固定子，回転子共に補正を行うため，誘導機のようなすべりがある場合は適用できません。
　誘導機ではステータ側である電機子電流の一周期とロータ側のすべり周期が一致しませんが，定常状態においては固定子のみで見た場合，電機子電流と磁場は近似的には半周期性（一周期の半分で符号が反転する）を有すると考えられます。回転子のロータバーとエンドリングに流れる電流の時定数は，固定子巻線の時定数に比べて小さいと考えられ，また回転子による非周期性も小さいとして，固定子のみの半周期の時間周期法を適用することが可能であると考えられます。そこでこの半周期性を利用し，固定子のみに時間周期法を適用する機能を開発しましたので報告します。

　適用する方法は簡易版時間周期法と呼ばれるもので，SD-EEC法では補正方程式を導出して補正量を求め，半周期もしくは一周期毎に補正を行うものですが，簡易版時間周期法の場合，半周期性に対して次の補正を行います。時間ステップ解析において，ある時点の方程式の解x₀ と半周期後の解 x₁に対して、補正量 pを次式で求めます。
　

　

　　誘導機の解析で示したモデルに，ロータバーとエンドリングを等価回路としてモデル化したかご型誘導機の二次元解析におけるロータバーとエンドリングの取り扱いのモデルに簡易版時間周期法を適用した例を示します。ロータバーは別々のプロパティで作成し，それぞれSUFCURとして定義し，NETWORKによりSUFCURをエンドリングを模擬した仮想抵抗で接続した等価回路としています。上記方法の威力を見るため，初期値はゼロとして時間ステップ解析を行います。すべりは0.25で，時間ステップはやや粗いですが一周期34ステップとし，半周期である17ステップ目で補正を行っています。
　　通常の方法（Conventional）では過渡状態が次第に減衰していきますが，未だ定常には達していないのに対し，簡易版時間周期法（Simplified EEC）では，補正4回ほどでほぼ定常解が得られていることがわかります。

Fig.2　回転数 1125rpm，		Fig.3　回転数 1125rpm，
すべり S=0.25 時の固定子一次巻き線の電流変化		すべり S=0.25 時のロータバーの渦電流損失

Fig.4　回転数 1125rpm，		Fig.5　回転数 1125rpm，
すべり S=0.25 時の回転子（半分）に働くトルク		すべり S=0.25 時のロータバー電流変化

　これより，誘導機の定常解析に簡易版時間周期法が有効であることを示せたと思います。誘導機に適用する場合，先にも説明しましたように，この方法は固定子のみに補正を行う方法であり，固定子の電機子電流と磁場を半周期性であるとして解くものですが，半周期ごとにステータとロータバーの位置関係が必ずしも一致しないため，かなり定常解には近くなっていると思いますが厳密には定常解とはなりません。そのため，周期的に印加される磁場に対する解の収束性と同様の考え方で，最初は簡易版時間周期法で定常解に近いところまで計算し，それを初期値としてリスタート解析を行って定常解を得る，といったようにお使い頂ければと思います。
　またSD-EEC法では，発電機や大型の界磁巻線型同期機でかご（短絡環）のロータバーをSUFCURとして定義して解析することはできませんが，上記誘導機の例で示しましたように，簡易時間周期法であれば解析することができます。

使用法
　inputファイルにおいてSD-EEC時間周期法のデータの後にTP_EEC_OPTION=1と入力するのみで簡易版時間周期法が適用されます。解析は半周期性を持つ過渡解析(TRANSIENT)計算のみで使用できます。時間周期法の設定は，時間周期問題の定常解への高速収束と同様です。

入力例


　また，時間周期問題の定常解への高速収束と同様に，特にNETWORKを使用する場合はTHETAとTHETA_NETWORKを，さらにMOTIONを使用する場合はTHETA_MOTIONを同じ1とした方が安定します。

入力例